Pinus massoniana is the major afforestation and vegetation restoration tree in southern China, and it plays an important role in the sustainable development of plantations. However, long-term single planting of P. massoniana has resulted in the decline of soil quality and forest productivity, and a soil fertility assessment is urgently needed. We selected P. massoniana plantations of four age stages for plot investigation and sampling to determine the soil physicochemical properties, microbial diversity and composition, and enzyme activities at different soil depths. The results showed that soil total phosphorus (TP) and available phosphorus (AP) decreased with the increase of age, especially low C/N ratio and high C/P and N/P ratio in the 30-year and 36-year stands, leading to P limitation. Meanwhile, the bacterial Shannon index also decreased with the increase of age and was positively correlated with AP, NO₃⁻-N, and pH. However, the fungal Shannon index decreased first and then increased with the increase of age; soil acid phosphatase (S-ACP) and urease activities showed a similar trend. Correlation analysis demonstrated that the increase of total organic carbon (TOC) and total nitrogen (TN) promoted the increase of fungal Shannon index, which was beneficial to the secretion of more enzymes. We found that soil physicochemical properties, microbial diversity, and enzyme activity decreased simultaneously when soil depths increased. Moreover, Acidobacteria and Basidiomycota were the most abundant bacterial and fungal communities, respectively, followed by Proteobacteria and Actinobacteria for bacteria and Ascomycota for fungi, and these microbial taxa were significantly affected by soil water content (SWC), TOC, AP, and C/P. In conclusion, this work reveals the potential correlation among soil physicochemical properties, microbial diversity and composition, and enzyme activities, and revealed potential correlations among them which will help to improve understanding of soil conditions and provide a reference for rational management of soil resources.

Background: Rhododendron delavayi Franch. var. delavayi is a wild ornamental plant species in Guizhou Province, China. The lack of its plastid genome information seriously hinders the further application and conservation of the valuable resource. Methods: The complete plastid genome of R. delavayi was assembled from long sequence reads. The genome was then characterized, and compared with those of other photosynthetic Ericaceae species. Results: The plastid genome of R. delavayi has a typical quadripartite structure, and a length of 202,169 bp. It contains a large number of repeat sequences and shows preference for codon usage. The comparative analysis revealed the irregular recombination of gene sets, including rearrangement and inversion, in the large single copy region. The extreme expansion of the inverted repeat region shortened the small single copy, and expanded the full length of the genome. In addition, consistent with traditional taxonomy, R. delavayi with nine other species of the same family were clustered into Ericaceae based on the homologous protein-coding sequences of the plastid genomes. Thus, the long-read assembly of the plastid genome of R. delavayi would provide basic information for the further study of the evolution, genetic diversity, and conservation of R. delavayi and its relatives.

Four evergreen broadleaf Rhododendron spp. (Rhododendrons), namely, Rhododendron aganniphum, R. nyingchiense, R. wardii, and R. triflorum, occur in harsh subalpine habitats in the southwest Qinghai-Tibet Plateau (QTP), China. Considering that the four Rhododendrons cannot escape their unique environment, they must evolve a set of adaptations to survive, but the information is lacking. To uncover their physiological adaptation characteristics, in the present study, we monitored their physiological characteristics by determination of their seasonal variation in antioxidant enzyme activity, osmotic adjustment substrates, and carbohydrate contents, and their pigment content and photosynthetic efficiency. The results showed that superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) activities and proline content of four Rhododendrons had a significant difference in autumn and were insignificant in summer. Specifically, R. aganniphum had the maximum protective enzyme activity and proline content in winter as well as chl a, b, and car contents. The values of maximal quantum yield (Fv/Fm), photochemical efficiency (ΦPSII), and non-photochemical quenching (NPQ) of four Rhododendrons were significantly higher in summer than in other seasons. The lower qP indicated the four Rhododendrons were susceptible to photoinhibition. Overall, the four Rhododendrons had similar physical characteristics in subalpine habitats. The parameters of the maximum quantum yield of photosystem II (PSII), the actual quantum yield of PSII, the non-photochemical chlorophyll fluorescence quenching, and chlorophyll a content increased in summer. Meanwhile, the protective enzyme activity and total soluble sugar content, proline content, and carotenoid content increased in spring, autumn, and winter. These results suggested that the four Rhododendrons can adapt to subalpine habitats by heat dissipation to avoid the damage of excessive radiation during the warm season while scavenging reactive oxygen and increasing the intracellular fluid concentration to avoid damage caused by chilling temperatures during the cold seasons. These findings would provide a reference for the conservation and application of these valuable ornamental evergreen broadleaf Rhododendrons, and enrich theory of plant eco-physiology in the high altitudes of the QTP.

Flower color is one of the crucial traits of ornamental plants. Rhododendron delavayi Franch. is a famous ornamental plant species distributed in the mountain areas of Southwest China. This plant has red inflorescence and young branchlets. However, the molecular basis of the color formation of R. delavayi is unclear. In this study, 184 MYB genes were identified based on the released genome of R. delavayi. These genes included 78 1R-MYB, 101 R2R3-MYB, 4 3R-MYB, and 1 4R-MYB. The MYBs were divided into 35 subgroups using phylogenetic analysis of the MYBs of Arabidopsis thaliana. The members of the same subgroup in R. delavayi had similar conserved domains and motifs, gene structures, and promoter cis-acting elements, which indicate their relatively conserved function. In addition, transcriptome based on unique molecular identifier strategy and color difference of the spotted petals, unspotted petals, spotted throat, unspotted throat, and branchlet cortex were detected. Results showed significant differences in the expression levels of R2R3-MYB genes. Weighted co-expression network analysis between transcriptome and chromatic aberration values of five types of red samples showed that the MYBs were the most important TFs involved in the color formation, of which seven were R2R3-MYB, and three were 1R-MYB. Two R2R3-MYB (DUH019226.1 and DUH019400.1) had the highest connectivity in the whole regulation network, and they were identified as hub genes for red color formation. These two MYB hub genes provide references for the study of transcriptional regulation of the red color formation of R. delavayi.