Histopathological data for each tissue type with observable lesions were analyzed using bivariable ordered logistic regression models, with treatment (high or low NO3-N) as the independent variable in each model and specific lesion score as the ordinal dependent variable