==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 79043 10.1038/s41598-020-79043-1 Article Acid resistance of Masson pine (Pinus massoniana Lamb.) families and their root morphology and physiological response to simulated acid deposition Zhou Sijie 12 Zhang Min 12 Chen Shuzhan 12 Xu Wen 12 Zhu Liting 12 Gong Shurui 12 He Xiaoqin 12 Wang Ping 18061713169@163.com 12 1 grid.410625.40000 0001 2293 4910College of Biology and the Environment, Nanjing Forestry University, Nanjing, 210037 China 2 grid.410625.40000 0001 2293 4910Key Laboratory of Forest Genetics and Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, 210037 China 16 12 2020 16 12 2020 2020 10 2206631 1 2020 4 11 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.Pinus massoniana Lamb. is one of the most sensitive species to acid deposition among forest woody plants, but differences in acid resistance among pine families still exist. It is of great significance to study the differences in acid resistance of Masson pine families and to analyze the physiological regulation mechanism of their acid resistance. In this study, the 100-day-old seedlings of 16 Masson pine families were treated with the simulated acid rain (SAR) at different pH levels (5.6, 4.5, 3.5 and 2.5) for 100 days to investigate the plant morphology, chlorophyll content, and root physiological responses. Results showed that pine family No. 35 maintained the good morphology, high chlorophyll content and organic acids secretion, and low plasma membrane permeability exposed to SAR, while family No. 79 presented the opposite. SAR not only increased the root plasma membrane permeability, but also induced an exudation of organic acids from the pine roots, and the test parameters changed sharply when the SAR pH was lower than 4.5. The results indicated that Masson pine could resist to acidic environment (pH 4.5–5.6), and family No. 35 had the acid resistance while the family No. 79 was sensitive to acid stress. The acid resistance diversity of different pine families had close relation with the root physiological processes, including the root plasma membrane permeability and organic acids secretion. For the future research, the natural genetic variation of Masson pine in response to acid stress and its acid resistance mechanism should be further studied. Subject terms EcologyEnvironmental scienceshttp://dx.doi.org/10.13039/501100001809National Natural Science Foundation of China31270680Wang Ping issue-copyright-statement© The Author(s) 2020 ==== Body Introduction The rapid growth of economy and population has resulted in the excessive consumption of fossil energy, and it consequently increased the emission of SOx and NOx, leading to the atmospheric acid deposition. Acid rain pollution has become one of the global top ten serious environmental challenges since the 1970s1. After Europe and North America, China is the third-largest heavy acid deposition region in the world2. Approximately 40% of China suffered from acid precipitation, and the incidence of acid rain in cities reached up to 19.8% in 20163,4. Busch et al.5 (2001) predicted that acid deposition in forest ecosystems would be about 46% higher in 2040 than that was in the 1980s. Acid deposition not only seriously threatens the biochemical cycle of the earth, but also destroys the stability of terrestrial ecosystems6. As the important part of the terrestrial ecosystem, plants can be considered as the main victim of acid deposition pollution7. Previous studies found that the damage mechanism of acid deposition on plants mainly involved destroying the ultrastructure of chloroplasts, inhibiting photosynthetic capacity, inducing the membrane lipid peroxidation, restraining respiration and other physiological metabolisms8–10. The soil acidification and the release of soil base ions caused by acid rain erosion can seriously damage plant roots, such as inhibiting root elongation and reducing the number of lateral roots7, thereby weakening their ability to absorb and transport water and nutrients11. These inhibiting effects will eventually cause the wilting, fading and death of plants. Masson pine (Pinus massoniana Lamb.) is a native tree species in China with a wide distribution range and large forest area. It is characterized by drought resistance, rapid growth, strong adaptability and high economic value and is a pioneer tree species for afforestation and reforestation12. However, large-scale observations showed that the regions that are seriously polluted by acid deposition are concentrated mainly in the south and the southwest of China, especially in city-clusters such as the Yangtze River Delta area, the Pearl River Delta area, and the Sichuan Basin13. Where the annual average pH of precipitation has been generally below 4.5 in recent years, even as low as 3.014. Wu et al.15 (1998) reported that Masson pine was sensitive to acid deposition, and the yield decreased by 43% when exposed to acid rain of pH ≤ 4.0. But Quan and Ding16 (2017) claimed that Masson pine showed the strong adaptability to the acidic soils in south China. Under the long-term natural selection and reproductive isolation, Masson pine had abundant intra species natural genetic variation. He et al.17 (2013) reported that the physiological responses of Masson pine to different levels of phosphorus stress exhibited highly significant variation among 11 elite pine families. But the differences in acid resistance of Masson pine families is not well-understood. Therefore, studying the differences in acid resistance of the pine families and finding scientific ways to alleviate the damage on Masson pine caused by acid deposition is of great significance for the restoration of forest ecosystems eroded by acid deposition. When plants are subjected to stress, the first response is the plasma membrane and the functional protein cells expressed on the plasma membrane18. Consequently, the relative permeability of root plasma membrane can directly reflect the stability of plant root cells internal environment and the adaptability and resistance of plant to external soil environment changes. Moreover, acid rain has also been reported to affect the plant chlorophyll content and photosynthesis as well19. Low molecular weight organic acids (e.g., oxalic acid (OA), malic acid (MA), citric acid (CA), etc.) are important components of root exudate, which play a key role in plant vivo osmosis, pH regulation, ion balance, and other rhizospheric reactions in response to environmental stress20. Root organic acids can also act in reducing the high oxidation state of some metal elements [e.g., Fe (III) and Cr (VI), etc.], so that they can be better absorbed and utilized by plants21. Besides, Zhang et al.22 (2018) found that citric acid and fumaric acid significantly stimulated Hansschlegelia zhihuaiae S113 colonization on cucumber roots, which allowed for improved cucumber growth. However, both the types and contents of organic acids secretion are affected by plant environmental stress and plant species23. Since the genetic variation of Masson pine is complex and diverse, the responses of different families to environmental stress are also different. Hence, studying the organic acids secretion of Masson pine families with different acid resistance under acid deposition is very important to further understand the physiological mechanism of plant acid resistance. In this work, 100-day-old seedlings of 16 Masson pine families were used as the test materials to study the response of plant morphology, photosynthetic traits, and root physiological parameters of Masson pine families to the simulated acid rain with different pH values. The purpose of this study is to investigate the differences in the acid resistance of pine families and to identify the key physiological factors of Masson pine in response to simulated acid deposition. Results Changes in morphology indexes of Masson pine families When exposed to simulated acid rain of different pH, the family No. 35 has always had longer root and shoot length, heavier root dry weight and more primary lateral root number among all 16 pine families, while families Nos. 79, 114 and 116 showed the opposite performance (Fig. 1). Although the length of roots and shoots (Fig. 1A,B), the dry weight of roots (Fig. 1C), and the number of primary lateral roots (Fig. 1D) showed different degrees of decrease with the decrease of pH value of SAR, the pine morphology was slightly effected by the SAR at pH 5.6–4.5, while under the acid rain gradient of pH 2.5–3.5, the pine morphology indexes decreased obviously.Figure 1 The morphology indexes, (A) root length (l/cm), (B) shoot length (l/cm), (C) root dry weight (m/mg), and (D) primary lateral root number of Masson pine families exposed to the SAR at pH of 2.5, 3.5, 4.5 and 5.6. Results are expressed as mean ± standard error (n = 3). Changes of chlorophyll contents for the pine seedlings The contents of chlorophyll a, chlorophyll b and chlorophyll a + b decreased with the increase of acid intensity, however, the value of chlorophyll a/b increased (Fig. 2). In addition, the chlorophyll b, a + b contents of the all families under severe acid stress of pH 2.5 were lower than those under moderate (pH 3.5), mild (pH 4.5) acid treatments and control treatment (pH 5.6), and the chlorophyll a/b value under severe acid stress was significantly higher than those of the other three treatments. Moreover, treated with the SAR of different pH, family No. 35 had higher contents of chlorophyll a, chlorophyll b, chlorophyll a + b and lower value of chlorophyll a/b, while the opposite phenomena were observed in the families Nos. 79, 114, and 116.Figure 2 Contents of (A) chlorophyll a (w/mg g−1), (B) chlorophyll b (w/mg g−1), (C) chlorophyll a + b (w/mg g−1), and (D) the value of chlorophyll a/b for the Masson pine families exposed to the SAR at pH of 2.5, 3.5, 4.5, and 5.6, respectively. Results are expressed as mean ± standard error (n = 3). Changes in the relative permeability of root plasma membrane The permeability and the damage degree of root plasma membrane of the Masson pine families both increased with the increase of acid intensity (Fig. 3). In addition, Masson pine family No. 35 kept low permeability and damage degree of root plasma membrane among all test pine families, while families Nos. 79, 114 and 116 presented opposite performances.Figure 3 (A) The relative permeability (p/%) and (B) damage degree (p/%) of root plasma membrane of 16 Masson pine families exposed to the SAR. The value of membrane damage degree of Masson pine families at pH 5.6 is 0, which is not shown in (B). Results are expressed as mean ± standard error (n = 3). Changes in the root secretion of organic acids and the rhizosphere pH of Masson pine family No. 35 and No. 79 In order to reflect the difference of root organic acids secretion of Masson pine families with different acid resistance strongly and intuitively, the pine family No. 35 and No. 79 were selected as the acid-resistant and acid-sensitive species respectively, to further study their root physiological response to the SAR (Fig. 4). The results of the ANOVA are tabulated in Table 1.Figure 4 Root organic acids secretion, (A) oxalic acid (w/ng L−1), (B) malic acid (w/ng L−1), (C) citric acid (w/ng L−1), and (D) soil pH of Pinus massoniana families Nos. 35 and 79. Based on Duncan’s multiple range tests with P < 0.05, significant differences between Masson pine family No. 35 and No. 79 are represented by lowercase letters. Results are expressed as mean ± standard error (n = 3). Table 1 ANOVA results of root organic acids secretion and soil pH of Masson pine exposed to the SAR (at 95% confidence level) (n = 3). SAR pH Source SS DF MS F P Oxalic acid secretion 2.5 Masson pine family 40,639.740 1 40,639.740 7.153 0.056 Error 22,725.000 4 5681.250 Total 63,364.740 5 3.5 Masson pine family 28,731.840 1 28,731.840 8.820 0.041 Error 13,029.700 4 3257.425 Total 41,761.540 5 4.5 Masson pine family 17,109.360 1 17,109.360 16.050 0.016 Error 4264.100 4 1066.025 Total 21,373.460 5 5.6 Masson pine family 1209.840 1 1209.840 2.100 0.221 Error 2304.020 4 576.005 Total 3513.860 5 Malic acid secretion 2.5 Masson pine family 7434.240 1 7434.240 100.701 0.001 Error 295.300 4 73.825 Total 7729.540 5 3.5 Masson pine family 5063.415 1 5063.415 40.216 0.003 Error 503.620 4 125.905 Total 5567.035 5 4.5 Masson pine family 3484.860 1 3484.860 37.731 0.004 Error 369.440 4 92.360 Total 3854.300 5 5.6 Masson pine family 1117.935 1 1117.935 14.678 0.019 Error 304.660 4 76.165 Total 1422.595 5 7434.240 Citric acid secretion 2.5 Masson pine family 289.815 1 289.815 3.637 0.129 Error 318.760 4 79.690 Total 608.575 5 3.5 Masson pine family 636.540 1 636.540 7.573 0.051 Error 336.200 4 84.050 Total 972.740 5 4.5 Masson pine family 393.660 1 393.660 5.308 0.083 Error 296.680 4 74.170 Total 690.340 5 5.6 Masson pine family 174.960 1 174.960 5.151 0.086 Error 135.860 4 33.965 Total 310.820 5 289.815 Rhizospheric soil pH 2.5 Masson pine family 0.558 1 0.558 9.264 0.038 Error 0.241 4 0.060 Total 0.799 5 3.5 Masson pine family 0.375 1 0.375 3.662 0.128 Error 0.410 4 0.102 Total 0.785 5 4.5 Masson pine family 0.018 1 0.018 0.161 0.709 Error 0.451 4 0.113 Total 0.470 5 5.6 Masson pine family 0.001 1 0.001 0.018 0.899 Error 0.293 4 0.073 Total 0.294 5 Non-rhizospheric soil pH 2.5 Masson pine family 0.228 1 0.228 2.700 0.176 Error 0.338 4 0.085 Total 0.566 5 3.5 Masson pine family 0.240 1 0.240 3.529 0.133 Error 0.272 4 0.068 Total 0.512 5 4.5 Masson pine family 0.029 1 0.029 0.325 0.599 Error 0.362 4 0.090 Total 0.391 5 5.6 Masson pine family 0.005 1 0.005 0.069 0.806 Error 0.315 4 0.079 Total 0.320 5 SS is the sum of squares, DF is the degrees of freedom, and MS is the mean square. The secretion of root organic acids of Masson pine increased with the pH value of SAR decreased from 5.6 to 3.5, and the secretion of oxalic acid and malic acid was obviously higher than that of citric acid according to Fig. 4A–C. No succinic acid and acetic acid were detected (Table S1). The root organic acids (oxalic acid, malic acid, and citric acid) secretion of Masson pine family No. 35 was higher than that of No. 79 exposed to the SAR of different acidity. Specifically, the root oxalic acid secretion of pine family No. 35 was significantly higher than that of No. 79 treated with the SAR at pH 3.5 (P = 0.041 < 0.05, Table 1) and 4.5 (P = 0.016 < 0.05, Table 1). Similarly, the root malic acid secretion of pine family No. 35 was significantly higher than that of No. 79 treated with the SAR at pH 2.5 (P = 0.001 < 0.05, Table 1), 3.5 (P = 0.003 < 0.05, Table 1) and 4.5 (P = 0.004 < 0.05, Table 1). The comparison between the soil background pH (6.65, Table S3) and the soil pH of Masson pine seedlings subjected to SAR (Fig. 4D, Table S2) showed that the rhizosphere and non-rhizosphere pH were lower and decreased with the decrease of SAR pH. Figure 4D also showed that although the soil pH of family No. 35 decreased with the decrease of SAR pH, it was still higher than that of No. 79 treated with SAR at the same acidity level. Particularly, the rhizospheric soil pH of family No. 35 treated with SAR of pH 2.5 was significantly higher than that of No. 79 (P = 0.038 < 0.05, Table 1). Correlation coefficients, principle component analysis (PCA), and cluster analysis Correlations among the morphology indexes, photosynthetic traits, and root physiological parameters of Masson pine families subjected to SAR were revealed by Pearson’s correlation coefficients. The significant negative correlations between oxalic acid secretion and root length (− 0.726*), shoot length (− 0.709*), rhizosphere pH (− 0.924**), and non-rhizosphere pH (− 0.948**) were observed (Table 2). The root plasma membrane permeability was significantly negatively correlated with photosynthetic traits (except chlorophyll a/b), morphology indexes, rhizospheric and non-rhizospheric soil pH (P < 0.01, Table 2). The rhizospheric soil pH was significantly negatively correlated with the chlorophyll a/b (− 0.860**, Table 2).Table 2 Correlation coefficients (r) among the photosynthetic traits, morphology indexes, and root physiological parameters for the pine seedlings of families Nos. 35 and 79, which were exposed to the SAR at pH of 5.6, 4.5, 3.5, and 2.5. Variables Chl. a/b RPRPM OA MA CA Chl. a − 0.726* − 0.865** − 0.551 − 0.099 0.052 Chl. b − 0.959** − 0.966** − 0.693 0.021 0.087 Chl. a + b − 0.899** − 0.959** − 0.661 − 0.026 0.076 Chl. a/b 1 0.928** 0.650 − 0.172 − 0.170 RPRPM 0.928** 1 0.798* 0.077 0.046 RL − 0.935** − 0.982** − 0.726* 0.027 0.049 SL − 0.936** − 0.966** − 0.709* 0.063 0.071 RDW − 0.931** − 0.958** − 0.698 − 0.012 0.009 PLRN − 0.946** − 0.943** − 0.584 0.242 0.286 OA 0.650 0.798* 1 0.529 0.570 MA − 0.172 0.077 0.529 1 0.959** CA − 0.170 0.046 0.570 0.959** 1 RpH − 0.860** − 0.943** − 0.924** − 0.280 − 0.291 NRpH − 0.836** − 0.929** − 0.948** − 0.331 − 0.342 **, and * = significant at 0.01 and 0.05 levels, respectively. The underlined data indicates a significant negative correlation between the two indicators. Chl. a chlorophyll a, Chl. b chlorophyll b, Chl. a + b chlorophyll a + b, Chl. a/b chlorophyll a/b, RPRMP relative permeability of plasma membrane, RL root length, SL shoot length, RDW root dry weight, PLRN primary lateral root number, OA oxalic acid, MA malic acid, CA citric acid, RpH rhizosphere pH, NRpH non-rhizosphere pH. The abbreviations above are the same as the following. In PCA, the principal components 1 and 2 accounted for 94.83% of the observed variables (F1 = 76.33% and F2 = 18.50%, Fig. 5). F1 was characterized by high positive scores for photosynthetic traits (except chlorophyll a/b), morphology indexes, rhizosphere and non-rhizosphere pH, and high negative scores for chlorophyll a/b, relative permeability of root plasma membrane, and oxalic acid secretion. F2 was characterized by high positive scores for malic acid and citric acid secretion.Figure 5 Principle component analysis (PCA) showing the associations among the photosynthetic traits, morphology indexes, and root physiological parameters of Masson pine families Nos. 35 and 79. Cluster analysis (Fig. 6) demonstrated the groups of plant physiological attributes of Masson pine families, which responded by a similar way to SAR treatment. They were grouped into three distinct sections according to their similarities when exposed to the SAR. The morphology indexes (root length, shoot length, root dry weight, primary lateral root number), photosynthetic traits (chlorophyll a, chlorophyll b, chlorophyll a + b), rhizospheric pH and non-rhizospheric pH were affected collectively and were grouped as G-1, which were positively correlated with the SAR pH. And the chlorophyll a/b, root plasma membrane permeability and oxalic acid were grouped as G-3, they were negatively correlated with the SAR pH. Besides, Malic acid and citric acid secretion exhibited a similar pattern and were grouped as G-2.Figure 6 Dendrograms of cluster analysis showing the similarity for the test indicators of Masson pine family No. 35 and No. 79 exposed to the SAR at pH of 5.6, 4.5, 3.5, and 2.5. Discussion Differences in plant morphology, chlorophyll content, and root plasma membrane permeability of the 16 Masson pine families Our results in Fig. 1 showed that the pine seedlings morphology was closely related to the SAR acidity. When treated with SAR at pH below 3.5, the pine seedlings showed the phenomenon of shortened root and shoot length, reduced root dry weight and lateral roots. A large amount of H+ caused by SAR was a major factor that threatened the growth of trees, which was absorbed into plant cells through roots, inducing intracellular free radical and proton effects24, thus adversely affecting the root morphology. Also, acid rain erosion will cause soil acidification, and the leaching of salt-based cations in acidified soil will reduce the ability of soil to supply nutrients25, thus affecting the absorption and utilization of nutrients in Masson pine roots, and then leading to the decrease of seedling morphology indexes26. Thus, we can infer that Masson pine seedlings had a certain resistance to acid rain of pH 4.5, but had poor resistance to acid rain with pH lower than 3.5. When plants are exposed to acid rain, the synthesis of chlorophyll can be inhibited, which accordingly leads to negative effects on leaf photosynthesis and metabolism6,7. Liu et al.27 (2007) studied the changes of chlorophyll in Masson pine under continuous acid rain treatment, indicating that the photochemical efficiency and energy conversion efficiency were decreased in the case of acid rain erosion and strong light. Therefore, plant chlorophyll content can be used as an indicator of plant stress resistance to characterize the physiological situation of plants under environmental stress28. High chlorophyll content and low chlorophyll a/b value are the important traits possessed by excellent Masson pine families27. Generally, the higher the chlorophyll content, the stronger the photosynthesis, accordingly producing more photosynthetic products, which can promote the growth of plants. Liang and Zhang28 (2018) found that chlorophyll content of soybean (Glycine max) seedlings decreased with the application of acid rain (pH 4.5/pH 3.0). Similarly, our results (Fig. 2) showed that treated with SAR, the chlorophyll content in pine needles decreased and the ratio of chlorophyll a/b increased. Under severe (pH 2.5) and moderate (pH 3.5) SAR treatments, the chlorophyll contents of Masson pine were much lower than that of the control. It indicated that acid rain of pH 2.5 and 3.5 inhibited the synthesis of chlorophyll, causing Masson pine to fail to make full use of the absorbed light energy, weakening photosynthesis, and reducing the resistance of Masson pine to acid rain. Plant cell membrane plays an important role in regulating the transport of intracellular and extracellular substances and maintaining the normal physiological metabolism of cells, the permeability of which is closely related to the intensity of stress under adversity29. Wagatsuma30 (2017) proposed that the retardation of aluminum permeation through the plasma membrane lipid bilayer mechanism was a new concept of aluminum tolerance in plants. The results in Fig. 3 showed that the root plasma membrane permeability of Masson pine treated with SAR of pH 5.6 and 4.5 remained low, while the plasma membrane was severely damaged by SAR of pH 3.5 and 2.5. Hence, it can be inferred that SAR can destroy the integrity of plant plasma membrane, but Masson pine seedlings have a certain acid resistance and can maintain the integrity of membrane under mild SAR stress. The similar phenomena were reported by Liang et al.31 (2019) who found that the membrane permeability of rice and soybean leaves did not change exposed to SAR at pH 5.0, 4.5 or 4.0, but which increased when treated with SAR at pH 3.5, 3.0 or 2.5. Differences in plant morphology and physiological traits of the Masson pine families No. 35 and No. 79 The results in Figs. 1 and 2 showed that the Masson pine family No. 35 maintained better phenotypic morphology, higher content of chlorophyll a, chlorophyll b, and chlorophyll a + b, lower chlorophyll a/b value in response to SAR of different pH among 16 Masson pine families. On the contrary, the pine families Nos. 79, 114 and 116 presented the worse plant morphology and lower chlorophyll content subjected to simulated acid deposition compared with others. In general, plants that are resistant to environmental stress exhibit some good phenotypic and physiological traits to a certain extent24. When exposed to adverse environment, plants with good phenotypic morphology and high chlorophyll content can better absorb nutrients from soil and carry out photosynthesis, which is conducive to its growth under environmental stress26,28. Therefore, the results indicated that the acid resistance of Masson pine family No. 35 was stronger than that of families Nos. 79, 114 and 116. When plants are subjected to adversity stress, a large amount of reactive oxygen species are produced to attack the membrane system, leading to the peroxidation of the cytoplasmic membrane and changes in membrane lipid components32. Which directly affecting the fluidity and permeability of membrane, resulting in the leakage of intracellular electrolytes and small molecular organic substances32. According to Fig. 3, the stability and integrity of membrane of Masson pine families Nos. 79, 114 and 116 was damaged more severely than that of family No. 35 exposed to SAR of the same pH level. Serious damage to the root plasma membrane of Masson pine families Nos. 79, 114 and 116 was likely to cause the leakage of intracellular nutrients and electrolytes, accordingly making the Masson pine root physiological metabolism unable to carry out normally. But the plasma membrane integrity of pine family No. 35 maintained well, which was beneficial to its resistance to acid stress. Among studied 16 Masson pine families, family No. 35 and families Nos. 79, 114 and 116 showed the significant differences in phenotypic and physiological indexes responding to acid stress (Figs. 1, 2, 3). The metabolic system, especially root organic acids secretion, plays an important role in regulating the response of plants to abiotic stress. In order to reveal the key physiological factors working in mechanism of resistance to acid stress in Masson pine roots, pine families No. 35 and No. 79 were selected as relatively acid-resistant and acid-sensitive plants respectively to observe the differences in root organic acids secretion and soil pH in response to acid stress. Organic acids are important plant root exudates, which can regulate the osmotic balance of plant cells, intracellular pH and ionic concentration20. Root organic acids can interact with some metal ions, such as cadmium (Cd), aluminum (Al)33, gallium (Ga), copper (Cu), manganese (Mn) and lead (Pb)34, chelated to form organometallic complexes35, thereby reducing the toxicity of metal ions to plants36. Therefore, root secretion of organic acids can be used as a special protective mechanism for plants to cope with acid stress. We observed that induced by SAR, the root secretion of organic acids increased (Fig. 4A–C), that was the same as Ahmed et al.37 (2016) and Wang et al.38 (2006). This is probably because the large amount of organic acids secreted by pine roots can react with H+ and metal ions activated by SAR to alleviate the injury caused by SAR to plants. In addition, root organic acids can regulate the availability of soil nutrients which are difficult to use, making them easy to be absorbed and utilized by plants39. Hence, the increase of root organic acids secretion induced by acid stress is a physiological defense mechanism of plants to against adversity. According to Fig. 4 and Table 1, when treated with SAR of the same acidity, the secretion of organic acids (except citric acid) exuded from pine roots of family No. 35 was obviously higher than that of family No. 79. This is probably because SAR treatment increased the membrane permeability of pine family No. 79, destroyed the integrity of plasma membrane, and adversely affected the root physiological metabolism of Masson pine, resulting in the less secretion of root organic acids of the family No. 79. Thus, the result indicated that the root organic acids secretion also contributed to the acid resistance of Masson pine. Soil buffer system, plant nutrients and metal cations, respiratory CO2 production, relative permeability of roots plasma membrane, root secretion of organic acids and proton can lead to fluctuations in the rhizosphere pH40. Due to the existence of buffer systems such as hydroxyl aluminum and exchangeable cations in the soil, which can react with H+ to reduce soil acidification and increase the soil pH. According to Table S3 and Fig. 4D, the soil pH of Masson pine seedlings subjected to SAR were lower than the soil background pH and decreased with the decrease of SAR pH. Which illustrated that the simulated acid rain could cause soil acidification and affect the rhizosphere environment of Masson pine. Moreover, it was probably that SAR of pH 3.5 and 2.5 input a large amount of H+ into the soil, resulting in consumption of the buffer system and the weakness of the regulation, so that the soil acidification intensified. When subjected to acid rain, a large amount of soil H+ could be absorbed by plant roots through root plasma membrane, and other soil H+ could react with anions (e.g., HCO3−, OH−, etc.) and organic acids (e.g., OA, MA, CA, etc.) secreted by plant roots. Wang et al.41 (2006) found that the rhizospheric pH of wheat treated with aluminum was positively correlated with root organic acids secretion. Figure 4D and Table 1 showed that treated with SAR at the same acidity level, the rhizospheric soil pH of Masson pine family No. 35 was significantly higher than that of family No. 79. In addition, the Pearson’s correlation coefficients between rhizospheric soil pH and chlorophyll a/b (− 0.860**), relative permeability of plasma membrane (− 0.943**), oxalic acid (− 0.924**) indicated that rhizosphere pH was significantly negatively correlated with these physiological properties (P < 0.01, Table 2). Treated with the SAR at the same pH, the root plasma membrane permeability of family No. 35 was lower than that of family No. 79, so that it could prevent more soil H+ from entering the roots and transport more anions to react with the soil H+. Furthermore, the organic acids secreted from family No. 35 was much more than that of family No. 79, so that it could bind to soil H+ better. These could explain why the rhizospheric pH of family No. 35 was higher than that of family No. 79 when subjected to SAR at the same pH. It could be proved by the results that the relative permeability of root plasma membrane and root oxalic acid secretion were instrumental in the resistance of Masson pine responding to acid stress. The results in the correlation analysis (Table 2) of the test indicators indicated that the root plasma membrane permeability and oxalic acid secretion are negatively correlated with the morphology indexes, photosynthetic traits (except chlorophyll a/b), and soil pH of Masson pine families treated with SAR of different pH. Furthermore, the chlorophyll a/b, root plasma membrane permeability, and oxalic acid secretion were characterized as the principal components 1 (F1 = 76.33%, Fig. 5) by high negative scores according to PCA and were grouped as G-3 due to the negative correlations with the SAR pH based on the cluster analysis (Fig. 6). Therefore, the root plasma membrane permeability and root oxalic acid secretion of Masson pine can be considered as the key physiological factors to further study the pine root physiological mechanism resisting to acid stress. Conclusion The effects of SAR on the growth and root physiological processes of the 16 pine families were studied in this paper. Herein, we focus on elucidating the relationship between acid resistance of Masson pine families and root physiological processes. The results evidently proved that Masson pine could tolerate the SAR at pH from 4.5 to 5.6, but there was significant difference in acid resistance among the pine families. Masson pine family No. 35 resisted to acid stress, but families Nos. 79, 114 and 116 showed acid sensitivity. Based on the data obtained, family No. 35 had better morphological characters, much more chlorophyll content and organic acids secretion, lower root plasma membrane permeability, and higher rhizospheric soil pH compared with family No. 79. Additionally, root plasma membrane permeability and organic acids secretion (especially oxalic acid) contribute greatly to the acid resistance of Masson pine. Understanding how plant regulators play a role in acid resistance mechanism of Masson pine and exploring their interaction is of great significance to the cultivation of Masson Pine under acid deposition. In future studies, the natural genetic variation of Masson pine in response to acid stress should be further studied, and its acid resistance mechanism, including root physiological process, metabolic system, enzyme system, as well as genetic regulation needs to be established. Methods Experimental soil The woodland soil was collected at a soil depth of 0–20 cm from Nanjing Laoshan. The soil sample was sifted by 1 mm screen after being air dried, and base fertilizers containing of 0.2 g N (urea), 0.1 g P2O5 (KH2PO4) and 0.1 g K2O (K2SO4) were provided per kilogram of the soil, which were mixed and aged for 1 week and spared to be used. The soil physicochemical properties were shown in Supplementary as Table S3. Preparation of simulated acid rain Based on the monitoring data of acid rain in Nanjing from April to July 2017 and the average concentrations of anions and cations in natural precipitation in Southern China, H2SO4, HNO3 with KF, NaCl, Ca(NO3)2, Mg(NO3)2, and (NH4)2SO4 were used to prepare the SAR by gradual dilution method. The as-prepared ion concentrations were (μmol L−1): K+, 8; Na+, 21; Ca2+, 33; Mg2+, 5; NH4+, 44; Cl−, 21; F−, 8. During 1992–2003, the concentration of SO42− had an obvious decrease in Nanjing rainwater, while NO3− increased markedly42. The molar ratio of H2SO4 to HNO3 in SAR was set as 4:1. Referring to the maximum rainfall acidity (pH 2.1) that has emerged in China43 and the annual average pH value of acid rain in Jiangsu Province (pH 4.93) in 2017, the pH of SAR was set at 3 levels: pH 4.5 (mild acidification), pH 3.5 (moderate acidification), and pH 2.5 (severe acidification); pH 5.6 was used as the control condition. Three Masson pine per family were cultivated in each treatment to conduct three repetitions. Cultivation of Masson pine seedlings and sample collection Sixteen families of Masson pine seeds collected from Jiangsu, Zhejiang, and Sichuan were used with the designated numbers of 35, 38, 76, 78, 79, 80, 88, 89, 90, 114, 115, 116, 117, 119, 147 and 151, respectively44. The seeds were sterilized with 75% ethanol for 10 min, and then washed with deionized water for 5 times with 20 min each. A 24 h soak time was allowed in 40 °C warm water, and subsequently the seeds were seeded on wet cotton and placed in incubator at constant 25 °C for germination. The 10-day healthy seedlings after germination were selected to cultivate with soil in the experimental rhizobox (Supplementary, Fig. S1). The moisture content of cultivation soil was kept at about 60% during the growth of Masson pine. The spray amount, time and frequency of SAR were implemented according to the precipitation status of Nanjing from April to July, 2017. The sowing time of Masson pine seeds was 20 April, 2018, and the times of germination and seedlings transfer were on April 29 and May 15, 2018, respectively. The next day after transplantation, the survival condition was checked and replanted immediately to ensure a survival rate of transplantation was ≥ 99.0%. These seedlings grew in a sunlight chamber. Considering the rapid growth rate of Masson pine in seedling stage and the obvious effect of environmental stress, the experiments were conducted on August 22, 2018, and the samples were collected after 100-day treatments. Three seedlings of each Masson pine family per treatment were harvested and were carefully washed with tap water and deionized water three times, respectively. After that, the pine seedlings were separated into roots, stems, and needles. Rhizospheric and non-rhizospheric soils were excavated down to the same horizon (20 cm) in the rhizobox for separation analyses. Analytical methods Morphology indexes Plant samples were taken from each treatment, the root length and shoot length were measured, and the number of primary lateral roots with the length ≥ 2 cm grown from the main root was counted directly. Root samples were oven-dried at 60 °C until reaching a constant weight to measure the root dry weight. Chlorophyll content The contents of chlorophyll a and chlorophyll b were determined by spectrophotometry at 663 nm and 645 nm, respectively45. Samples from each family were selected randomly and the mid-upper needles on the sunny side of which were mixed. Exact 1 g sample was taken to be measured by UV-2401PC UV–visible spectrophotometer (Japan Shimadzu Company) after being extracted with 100 mL of acetone; repeated 3 times per family. Relative permeability of roots plasma membrane The relative permeability of roots plasma membrane was determined by conductivity method according to Inone and Kinoshita46 (2017). Fresh root of 0.1 g was washed with deionized water and cut into small pieces with length of 1 cm. Then, the root sample was put into a 20 mL stoppered test tube added with 10 mL of deionized water and oscillated in a shaking incubator for 30 min. The DDS-11C digital electric conductivity meter (Shenzhen Shengbang Electrical Appliance Factory, Tianjin, China) was used for determination of electric conductivity E1; then it was treated in boiling water for 10 min and the electric conductivity E2 was determined after cooling down; the background electric conductivity E0 was also measured. The relative plasma membrane permeability was calculated as: p (%) = (E1 − E0)/(E2 − E0) × 100%. The degree of membrane damage was calculated according to the following formula: the degree of membrane damage (p/%) = (PT − PCK) × 100%, where PT was the relative plasma membrane permeability of the treatments; PCK was the relative plasma membrane permeability of the control. Determination of organic acids The ultra-trace organic acids were determined by an accelerated solvent extraction-solid-phase extraction-liquid chromatography with electrospray ionization-tandem mass spectrometry methodology according to Wang et al.47 (2015). The detailed data were shown in Supplementary as Table S1. Soil pH After 100 days of seedling cultivation, the rhizospheric and non-rhizospheric soil samples (Supplementary, Fig. S1) were extracted by 1 mol L−1 KCl (2.5:1 water/soil ratio). The soil pH was measured by a PHS-25 digital acidity meter (Shanghai Precision Instrument Co., Ltd.) according to Yao et al.48 (2019). The detailed data were shown in Supplementary as Table S2. Statistical analysis Data were processed by One-Way ANOVA, Pearson’s correlation coefficients, principle component analysis, and cluster analysis using SPSS 11.5 software. The factor in one-way ANOVA was the Masson pine family and the significance (P < 0.05) of differences among the Masson pine families was evaluated by Duncan’s multiple range tests; tabulation and plotting were performed by Excel software and Origin Pro Version 2017 (OriginLab Corporation, Northampton, MA, USA). Supplementary Information Supplementary Information. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Information The online version contains supplementary material available at 10.1038/s41598-020-79043-1. Acknowledgements This research was financially supported by National Natural Science Foundation of China (grant numbers: 30771696, 31270680) and Jiangsu University Advantages Construction Project Funding Project. Dr. Wang Ping acknowledges support from the Analytical Center of Nanjing Forestry University. Author contributions S.Z. made the draft on main text, figures and tables and collected data and did experiments for this manuscript. M.Z. modified the main text. S.C. and W.X. and L.Z. did some experiments and collected some data for this manuscript. P.W. designed experiments, made the outline, and revised main text for this manuscript. S.G. and X.H. helped to collect samples and recorded some data for this manuscript. Competing interests The authors declare no competing interests. ==== Refs References 1. Reis S Grennfelt P Klimont Z Amann M ApSimon H Hettelingh JP Holland M LeGall AC Maas R Posch M Spranger T Sutton MA Williams M From acid rain to climate change Science 2012 338 1153 1154 10.1126/science.1226514 23197517 2. Wang L Chen Z Shang H Wang J Zhang PY Impact of simulated acid rain on soil microbial community function in Masson pine seedlings Electron. J. Biotechnol. 2014 17 199 203 10.1016/j.ejbt.2014.07.008 3. Wang WX Xu PJ Research progress in precipitation chemistry in China Prog. Chem. 2010 21 266 281 4. Meng Y Zhao YL Li R Li JL Cui LL Kong LD Fu HB Characterization of inorganic ions in rainwater in the megacity of Shanghai: Spatiotemporal variations and source apportionment Atmos. Res. 2019 222 12 24 10.1016/j.atmosres.2019.01.023 5. Busch G Lammel G Beese FO Feichter J Dentener FJ Roelofs GJ Forest ecosystems and the changing patterns of nitrogen input and acid deposition today and in the future based on a scenario Environ. Sci. Pollut. Res. 2001 8 95 102 10.1007/BF02987301 6. Wang Y Zhang XL Hu YB Teng ZY Zhang SB Chi Q Sun GY Phenotypic response of tobacco leaves to simulated acid rain and its impact on photosynthesis Int. J. Agric. Biol. 2019 21 391 398 7. Ramlall C Varghese B Ramdhani S Pammenter NW Bhatt A Berjak P Sershen P Effects of simulated acid rain on germination, seedling growth and oxidative metabolism of recalcitrant-seeded Trichilia dregeana grown in its natural seed bank Physiol. Plant. 2015 153 149 160 10.1111/ppl.12230 24835442 8. Wang XQ Liu Z Niu L Fu B Long-term effects of simulated acid rain stress on a staple forest plant, Pinus massoniana Lamb: A proteomic analysis Trees Struct. Funct. 2013 27 297 309 10.1007/s00468-012-0799-z 9. Tong SM Zhang LQ Differential sensitivity of growth and net photosynthetic rates in five tree species seedlings under simulated acid rain stress Pol. J. Environ. Stud. 2014 23 2259 2264 10.15244/pjoes/24930 10. Wu X Liang CJ Enhancing tolerance of rice (Oryza sativa ) to simulated acid rain by exogenous abscisic acid Environ. Sci. Pollut. Res. 2017 24 4860 4870 10.1007/s11356-016-8219-3 11. Hu WJ Chen J Liu TW Wu Q Wang WH Liu X Shen ZJ Simon M Chen J Wu FH Pei ZM Zheng HL Proteome and calcium-related gene expression in Pinus massoniana needles in response to acid rain under different calcium levels Plant Soil 2014 380 285 303 10.1007/s11104-014-2086-9 12. Luo SP He BH Zeng QP Li NJ Yang L Effects of seasonal variation on soil microbial community structure and enzyme activity in a Masson pine forest in Southwest China J. Mt. Sci. 2020 17 1398 1409 10.1007/s11629-019-5825-9 13. Zhang MY Wang SJ Wu FC Yuan XH Zhang Y Chemical compositions of wet precipitation and anthropogenic influences at a developing urban site in southeastern China Atmos. Res. 2007 84 311 322 10.1016/j.atmosres.2006.09.003 14. Li YF Wang YJ Wang B Wang YQ Response of soil respiration and its components to simulated acid rain in a typical forest stand in the three gorges reservoir area Environ. Sci. 2019 40 1457 1467 10.13227/j.hjkx.201803170 15. Wu, G. Effect of acidic deposition on productivity of forest ecosystem and estimation of its economic losses in southern suburbs of Chongqing China. J. Environ. Sci-China10, 83–88. http://kns.cnki.net/kns/detail/detail.aspx?FileName=HJKB802.010&DbName=CJFQ1998 (1998). 16. Quan WX Ding GJ Root tip structure and volatile organic compound responses to drought stress in Masson pine (Pinus massoniana Lamb.) Acta. Physiol. Plant. 2017 39 258 10.1007/s11738-017-2558-7 17. He YL Liu AQ Mulualem T Wu PF Ma XQ Wang C Per CO Physiological responses of needles of Pinus massoniana elite families to phosphorus stress in acid soil J. For. Res. 2013 24 325 332 10.1007/s11676-013-0356-7 18. DeHayes DH Schaberg PG Hawley GJ Strimbeck GR Acid rain impacts on calcium nutrition and forest health Bioscience 1999 49 789 800 10.2307/1313570 19. Ju SM Wang LP Chen JY Effects of silicon on the growth, photosynthesis and chloroplast ultrastructure of Oryza sativa L. seedlings under acid rain stress Silicon 2020 12 655 664 10.1007/s12633-019-00176-8 20. Ma Y Guo LQ Wang HX Bai B Shi DC Accumulation, distribution, and physiological contribution of oxalic acid and other solutes in an alkali-resistant forage plant, Kochia sieversiana , during adaptation to saline and alkaline conditions J. Plant Nutr. Soil Sci. 2011 174 655 663 10.1002/jpln.201000337 21. Rajniak J Giehl RFH Chang E Murgia I von Wiren N Sattely ES Biosynthesis of redox-active metabolites in response to iron deficiency in plants Nat. Chem. Biol. 2018 14 442 450 10.1038/s41589-018-0019-2 29581584 22. Zhang H Chen F Zhao HZ Lu JS Zhao MJ Hong Q Huang X Colonization on cucumber root and enhancement of chlorimuron-ethyl degradation in rhizosphere by Hansschlegelia zhihuaiae S113 and root exudates J. Agric. Food Chem. 2018 66 4584 4591 10.1021/acs.jafc.8b00041 29672047 23. Chen YT Wang Y Yeh KC Role of root exudates in metal acquisition and tolerance Curr. Opin. Plant Biol. 2017 39 66 72 10.1016/j.pbi.2017.06.004 28654805 24. Yan F Schubert S Mengel K Effect of low root medium pH on net proton release, root respiration, and root growth of corn (Zea mays L.) and broad bean (Vicia faba L.) Plant Physiol. 1992 99 415 421 10.1104/pp.99.2.415 16668900 25. Hu XF Wu AQ Wang FC Chen FS The effects of simulated acid rain on internal nutrient cycling and the ratios of Mg, Al, Ca, N, and P in tea plants of a subtropical plantation Environ. Monit. Assess. 2019 191 99 10.1007/s10661-019-7248-z 30680466 26. Ericsson T Growth and shoot: root ratio of seedlings in relation to nutrient availability Plant Soil 1995 168–169 205 214 10.1007/BF00029330 27. Liu JX Zhou GY Yang CW Ou ZY Peng CL Responses of chlorophyll fluorescence and xanthophyll cycle in leaves of Schima superba Gardn. & Champ. and Pinus massoniana Lamb. to simulated acid rain at Dinghushan biosphere reserve, china Acta Physiol. Plant. 2007 29 33 38 10.1007/s11738-006-0005-2 28. Liang CJ Zhang BJ Effect of exogenous calcium on growth, nutrients uptake and plasma membrane H+ -ATPase and Ca2+ -ATPase activities in soybean (Glycine max ) seedlings under simulated acid rain stress Ecotoxicol. Environ. Safe 2018 165 261 269 10.1016/j.ecoenv.2018.09.019 29. Li XW Li YL Mai JW Tao L Qu M Liu JY Shen RF Xu GL Feng YM Xiao HD Wu LS Shi L Guo SX Liang J Zhu YY He YM Baluska F Shabala S Yu M Boron alleviates aluminum toxicity by promoting root alkalization in transition zone via polar auxin transport Plant Physiol. 2018 177 1254 1266 10.1104/pp.18.00188 29784768 30. Wagatsuma T The membrane lipid bilayer as a regulated barrier to cope with detrimental ionic conditions: making new tolerant plant lines with altered membrane lipid bilayer Soil Sci. Plant Nutr. 2017 63 507 516 10.1080/00380768.2017.1369362 31. Liang CJ Ma YJ Li LR Comparison of plasma membrane H+ -ATPase response to acid rain stress between rice and soybean Environ. Sci. Pollut. Res. 2020 27 6389 6400 10.1007/s11356-019-07285-2 32. Guo Q Liu L Barkla BJ Membrane lipid remodeling in response to salinity Int. J. Mol. Sci. 2019 20 4264 10.3390/ijms20174264 6747501 33. Pellet DM Grunes DL Kochian LV Organic acid exudation as an aluminum-tolerance mechanism in maize (Zea mays L.) Planta 1995 196 788 795 10.1007/BF01106775 34. Wang HH Shan XQ Liu T Xie YN Wen B Zhang SZ Han F van Genuchten MT Organic acids enhance the uptake of lead by wheat roots Planta 2007 225 1483 1494 10.1007/s00425-006-0433-7 17106683 35. Li ZR Wang JX An LZ Tan JB Zhan FD Wu J Zu YQ Effect of root exudates of intercropping vicia faba and arabis alpina on accumulation and sub-cellular distribution of lead and cadmium Int. J. Phytoremediat. 2019 21 4 13 10.1080/15226514.2018.1523867 36. Jia H Hou DY Dai Y Lu HL Yan CL Effects of root exudates on the mobility of pyrene in mangrove sediment water system CATENA 2018 162 396 401 10.1016/j.catena.2017.10.022 37. Ahmed IM Nadira UA Cao FB He XY Zhang GP Wu FB Physiological and molecular analysis on root growth associated with the tolerance to aluminumand drought individual and combined in Tibetan wild and cultivated barley Planta 2016 243 973 985 10.1007/s00425-015-2442-x 26748913 38. Wang P Bi SP Wang S Ding QY Variation of wheat root exudates under aluminum stress J. Agric. Food Chem. 2006 54 10040 10046 10.1021/jf061249o 17177539 39. Yao Y Zhang FY Wang MY Liu F Liu WF Li XW Qin DD Geng XH Huang XX Zhang P Thallium-induced oxalate secretion from rice (Oryza sativa L.) root contributes to the reduction of Tl(III) to Tl(I) Environ. Exp. Bot. 2018 155 387 393 10.1016/j.envexpbot.2018.07.028 40. Javed M Akram M Habib N Tanwir K Ali Q Niazi N Gul H Iqbal N Deciphering the growth, organic acid exudations, and ionic homeostasis of Amaranthus viridis L. and Portulaca oleracea L. under lead chloride stress Environ. Sci. Pollut. Res. 2017 25 2958 2971 10.1007/s11356-017-0735-2 41. Wang P Bi SP Ma LP Han WY Aluminum tolerance of two wheat cultivars (Brevor and Atlas66 ) in relation to the irrhizosphere pH and organic acids exuded from roots J. Agric. Food. Chem. 2006 54 10033 10039 10.1021/jf0611769 17177538 42. Tu J Wang HS Zhang ZF Jin X Li WQ Trends in chemical composition of precipitation in Nanjing, China, during 1992–2003 Atmos. Res. 2005 73 283 298 10.1016/j.atmosres.2004.11.002 43. Liang CJ Wang WM Antioxidant response of soybean seedlings to joint stress of lanthanum and acid rain Environ. Sci. Pollut. Res. 2013 20 8182 8191 10.1007/s11356-013-1776-9 44. Tang, X. R., Li, W. P., Zuo, H. S. & Yin, Y. L. Study on the growth stability of Pinus Massoniana. J. Hunan For. Sci. Technol.29, 20–24, http://kns.cnki.net/kns/detail/detail.aspx?FileName=HLKJ200204005&DbName=CJFQ2002 (2002) (in Chinese). 45. Jia XM Wang H Svetla S Zhu YF Hu Y Cheng L Zhao T Wang YX Comparative physiological responses and adaptive strategies of apple Malus halliana to salt, alkali and saline-alkali stress Sci. Hortic. Amsterdam 2019 245 154 162 10.1016/j.scienta.2018.10.017 46. Inoue S Kinoshita T Blue light regulation of stomatal opening and the plasma membrane H+ -ATPase Plant Physiol. 2017 174 531 538 10.1104/pp.17.00166 28465463 47. Wang SL Fan CNQ Wang P Determination of ultra-trace organic acid in Masson pine (Pinus massoniana L.) by accelerated solvent extraction and liquid chromatography-tandem mass spectrometry J. Chromatogr. B 2015 981–982 1 8 48. Yao YW Ren BL Yang Y Huang CJ Li MY Preparation and electrochemical treatment application of Ce-PbO2 /ZrO2 composite electrode in the degradation of acridine orange by electrochemical advanced oxidation process J. Hazard. Mater. 2019 361 141 151 10.1016/j.jhazmat.2018.08.081 30179785