ResultsHeterogeneous nutrient environments, compared with a homogeneous environment, significantly increased the dry matter accumulation and root morphology indexes of the root system of F. hodginsii, which proliferated in nutrient-rich patches, and the P heterogeneous environment had the most significant enhancement effect, with dry matter accumulation 70.2%, 7.0%, and 27.0% higher than that in homogeneous and N and K heterogeneous environments, respectively. Homogeneous environments significantly increased the specific root length and root area of the root system; the dry matter mass and morphological structure of the root system of F. hodginsii with a heterospecific neighbor were higher than those under conspecific neighbor and single-plant treatments, and the root area of the root system under the conspecific neighbor treatment was higher than that under the heterospecific neighbor treatment, by 20% and 23%, respectively. Moreover, the root system under heterospecific neighbor treatment had high sensitivity; the heterogeneous nutrient environment increased the mean diameter of the fine roots of the seedlings of F. hodginsii and the diameter of the vascular bundle, and the effect was most significant in the P heterogeneous environment, exceeding that in the N and K heterogeneous environments. The effect was most significant in the P heterogeneous environment, which increased fine root diameter by 20.5% and 10.3%, respectively, compared with the homogeneous environment; in contrast, the fine root vascular ratio was highest in the homogeneous environment, and most of the indicators of the fine root anatomical structure in the nutrient-rich patches were of greater values than those in the nutrient-poor patches in the different heterogeneous environments; competition promoted most of the indicators of the fine root anatomical structure of F. hodginsii seedlings. According a principal component analysis (PCA), the N, Pm and K heterogeneous environments with heterospecific neighbors and the P heterogeneous environment with a conspecific neighbor had higher evaluation in the calculation of eigenvalues of the PCA.