- Open Access
Plant sex affects plant-microbiome assemblies of dioecious Populus cathayana trees under different soil nitrogen conditions
Microbiome volume 10, Article number: 191 (2022)
Dioecious plants have coevolved with diverse plant microbiomes, which are crucial for the fitness and productivity of their host. Sexual dimorphism in morphology, physiology, or gene expression may relate to different microbial compositions that affect male and female fitness in different environments. However, sex-specific impacts on ecological processes that control the microbiome assembly are not well known. In this study, Populus cathayana males and females were planted in different nitrogen conditions. It was hypothesized that males and females differently affect bacterial and fungal communities in the rhizosphere soil, roots, old leaves, and young leaves. Physiological traits and transcriptome profiles of male and female plants were investigated to reveal potential mechanisms that control the microbiome assembly.
Our results showed strong niche differentiation that shapes microbial communities leading to a rapid loss of diversity along a decreasing pH gradient from the rhizosphere soil to leaves. Sex had different impacts on the microbial assembly in each niche. Especially fungal endophytes showed great differences in the community structure, keystone species, and community complexity between P. cathayana males and females. For example, the fungal co-occurrence network was more complex and the alpha diversity was significantly higher in young female leaves compared to young male leaves. Transcriptome profiles revealed substantial differences in plant-pathogen interactions and physiological traits that clearly demonstrated divergent internal environments for endophytes inhabiting males and females. Starch and pH of young leaves significantly affected the abundance of Proteobacteria, while tannin and pH of roots showed significant effects on the abundance of Chloroflexi, Actinobacteria, and Proteobacteria, and on the bacterial Shannon diversity.
Our results provided important knowledge for understanding sexual dimorphism that affects microbial assemblies, thus advancing our understanding of plant-microbiome interactions.
Dioecious angiosperms, including 15,600 species in 987 genera and 175 families, have evolved independently from hermaphroditic ancestors . Genetic sex determination mechanisms that lead to separate male and female plants have been identified . Males and females have sex-specific adaptions to different environments resulting in biased sex ratios in populations [3,4,5]. Sexual dimorphism in morphology, physiology, and gene expression, and biased sex ratios are due to selective pressures acting differently on males and females [6,7,8]. Populus species are widely distributed, and they are suitable woody model species to study sex-specific responses of dioecious plants to different environments . Previous studies have provided convincing evidence of male poplars having a better capacity to forage limited nutrients in infertile soil and a higher resistance to other abiotic stresses when compared to female poplars [5, 7, 8].
Plant microbiota or microbiome comprises diverse and complex microbial communities containing bacteria, archaea, and fungi living on the plant surface (e.g., rhizosphere and phyllosphere) as well as inside plants (endosphere) [9,10,11]. Plants and their microbial inhabitants have coevolved over millions of years and are linked to plant fitness as a “holobiont” . Recent studies suggest that host health, growth, nutrient acquisition, and plant microbiota assembly are largely impacted by complex interactions between plants, microbes, and environments [11,12,13]. However, the underlying ecological processes in regulating plant-microbiome-environment interactions are still poorly understood.
Studies on dioecious plants have reported sex-specific responses to arbuscular mycorrhizal fungi suggesting that female growth and reproduction benefit more from arbuscular mycorrhizal fungi when compared to males . However, a higher colonization rate or hyphal biomass of arbuscular mycorrhizal fungi in male plants have been suggested to be crucial for males for the maintenance of a higher resistance or tolerance to different abiotic stresses [7, 15]. For example, a closer relationship between arbuscular mycorrhizal fungi and the root system of P. cathayana males has been suggested to enhance the capacity of phosphorus forage compared to females . A selective enrichment of a specific microbiome that contributes to stress resistance is driven by variation in plant metabolism and immune-related characters [16, 17]. However, sex-specific plant microbiomes have not been investigated, whereas differences in defense, growth, and reproduction between male and female plants have been well studied [6, 18,19,20].
The microbial community assembly of the rhizosphere is largely determined by different plant species or genotypes through sensing and responding to root-derived signals . A recent study revealed that Populus euphratica males and females harbor sex-specific bacterial and fungal communities in the rhizosphere soil of natural forests . It has been reported that a rapid loss of diversity from soil to roots and then to leaves and flowers indicate a strong selective pressure along the soil-plant continuum [10, 11, 23]. The microorganisms overcome or escape the host’s immune system and successfully colonize plant compartments (endosphere), such as roots, stems, or leaves without causing any disease. The community assembly of endophytes largely depends on the plant species, genotype, and development stage because of differences in plants’ physiology and metabolism [11, 19, 24]. Sexual differences in morphology, physiology, and metabolism between males and females are supposed to lead to distinct plant endophytes. The metabolism of different plant compartment niches generates chemically distinct environments and each host tissue provides a stable availability of metabolites for associated endophytes to maintain the biosynthesis of diverse compounds [23, 25]. Changes in plant metabolism caused by varied environmental conditions are crucial factors that affect microbial communities within plants.
As an important dioecious plant, P. cathayana has been extensively studied for its responses to different environmental conditions. It has been discovered that the sexual differences are greatest in poor environments . P. cathayana males have higher carbon fixation and transportation from a carbon source (leaf) to sinks, whereas P. cathayana females show higher secondary metabolic activities and glycolysis, but a lower expression of resistance genes under soil nitrogen-limited conditions [8, 20, 26]. In internal plant compartments, such as root and leaf tissues, a fine-tuning and adaptation of the microbiome is clearly evident , but the specific factors controlling the local microbial assembly and stability are still unclear. In the present study, the effects of internal traits on the microbial assembly and functions in P. cathayana males and females were explored. Three soil nitrogen levels were used to cultivate plant materials, and sampled rhizosphere soils, roots, old leaves, and young leaves. We firstly hypothesized that niche-dependent bacterial and fungal communities occur along the soil-plant continuum. Then, sex would impose distinct impacts on the bacterial and fungal communities in each niche. As previously shown, P. tremula females prioritize nutrient acquisition, and flavonoid and condensed tannin production, while males allocate more to the growth of aerial parts under limited nutrient conditions . Therefore, the sex effect on the microbial assembly of endophytes was supposed to be stronger when soil nitrogen level is limited. If the nitrogen level was recovered to a normal level, the sex effect would be expected to decrease.
Materials and methods
Populus cathayana forests widely distribute in Qinghai province, China. The original cuttings of P. cathayana males and females were sampled from five populations as detailed in a previous study . Five populations were separated by at least 50 km for each other. All cuttings were then cultivated in a greenhouse at the Hangzhou Normal University, Zhejiang Province, China. The whole experiment was conducted in this greenhouse with the temperature ranging between 21 and 25 °C during the day and 15 and 18 °C at night, with 12–14 h photoperiod . Soil with a soil organic matter content of 4 g kg−1, total nitrogen content 1.5 g kg−1, and pH 7.1 was collected in the garden and homogenized. Then, the collected soil, sand and vermiculite in a 1:1:1 ratio were homogenized to set “nutrient-limited” soil. We sampled 30 male and 30 female cuttings, and planted them in plastic pots (30 cm diameter and 21 cm height), which were filled with nutrient-limited soil at the beginning of March, 2019. Each pot contained one cutting. All plants were irrigated once a week with 150 ml Long Ashton solution containing 2 mM NH4NO3 until the beginning of April when the cuttings started root formation  (see Supplementary methods S1).
In all, 24 male and 21 female cuttings successfully sprouted, of which 15 male and 15 female individuals with a similar height for further experiments were evenly chosen from the five original populations. Three nitrogen treatments through controlling the N content of the Long Ashton solution were designed. The control-nitrogen treatment, re-irrigated with control-nitrogen solution treatment, and limited-nitrogen treatment were named as N1, N2, and N3, respectively. Thus, the experimental setup with two sexes (males and females) and three N treatments (N1, N2, and N3) included six treatments in total. Each treatment contained five replicates. The control-nitrogen treatment was irrigated with a solution containing 2 mM NH4NO3 (considered as a normal N supply) twice a week till the end of the experiment. Similarly, the limited-nitrogen treatment was controlled by excluding the N source from the solution when irrigated twice a week till end of the experiment. Two steps were used to construct the re-irrigated with control-nitrogen solution treatment. Firstly, 5 male and 5 female individuals were irrigated with 200 ml solution without N source twice a week until mid-July and then re-irrigated with the solution containing 2 mM NH4NO3 twice a week until end of the experiment.
Sample collection and soil traits
Rhizosphere soil, roots, and old and young leaves were collected at the beginning of September, 2019. Rhizosphere soil was defined as the soil tightly attached to plant roots. Roots of each plant were sampled. After the N treatment, 2 or 3 new emerging leaves of each plant were labelled and named them as old leaves since they began to be senescent by the end of the experiment. Recent fully expanded leaves of each plant present before sampling were named as young leaves. Both old and young leaves were collected. Soil, roots, old leaves, and young leaves were referred as “microhabitats” in this study. All samples were first kept on ice and stored at −80 °C before extracting DNA in the laboratory. The soil chemical traits including pH, total nitrogen (TN), total phosphorus (TP), available phosphorus (AP), NH4+, and NO3−, were measured according to standard protocols . Activities of β-1,4-glucosidase (BG) related to labile-C-cycling and β-1,4-N-acetylglucosaminidase (NAG) related with N-cycling were tested by Elisa kits (Shanghai, China). Soil microbial carbon and nitrogen biomass were determined by the chloroform fumigation extraction method .
DNA extraction and amplification
For the rhizosphere soil, DNA was extracted from 0.5 g soil using the E.Z.N.A.® soil DNA Kit (Omega Bio-tek, Norcross, GA, USA). The primers used for bacterial 16S rRNA (V3-V4) gene amplifications were 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) . Amplifications of the internal transcribed spacer region ITS2 in fungi were conducted using the primers ITS3F (5′-GCATCGATGAAGAACGCAGC-3′) and ITS4R (5′-TCCTCCGCTTATTGATATGC3′) . Before endophytic DNA extractions, leaves or roots were washed with sterile H2O and then treated with 70% ethanol for 5 min. After that, the samples were immersed in 5.25% NaClO solution for 5 min, and in 70% ethanol for 30 s. Finally, the samples were washed with sterile H2O . A two-round amplification process was used to amplify the 16S rRNA amplicon of bacteria inhabiting roots and leaves. The first round of PCR was performed using the primers 799F (5′-AACMGGATTAGATACCCKG-3′) and 1392R (5′-ACGGGCGGTGTGTRC-3′) . The second round of PCR was conducted with the primers 799F (5′-AACMGGATTAGATACCCKG-3′) and 1193R (5′-ACGTCATCCCCACCTTCC-3′) targeting the V5–V7 region . The amplification of endophytic fungi was also performed using the ITS3F/ITS4R primers. All samples were sequenced on an Illumina MiSeq platform (Illumina, San Diego, USA) with a paired end sequencing (2 × 300). Low-quality read ends and primer sequences were trimmed. Paired end sequences were merged to a single sequence by FLASH version 1.2.11 . Operational taxonomic units (OTUs) with 97% similarity cutoff were clustered using UPARSE version 7.1  after removing chimeric sequences and singletons. Representative sequences were classified using RDP Classifier version 11.5 against the SILVA reference database (v. 132) and UNITE (v. 8.0) to classify each OTU into bacteria and fungi, respectively. The raw data of bacterial and fungal sequences are available at the National Center for Biotechnology Information under BioProject ID PRJNA797907.
Physiological traits and transcriptome sequencing
Physiological traits of roots and young leaves, including pH, total phenols, tannin, NH4+, NO3−, and starch, were measured for four replicates per treatment. The remaining leaves and roots after DNA extraction were separated into two parts. One part was dried at 75 °C for 72 h and ground into fine powder, for which measurements of pH and starch were conducted. Mixtures of dried samples and deionized water (volume ratio 1: 8) were shaken at 250 rpm for 1 h and then centrifuged at 13,000 rpm for 5 min to test pH . Fine powder was extracted by 80% ethanol (v/v) at 80 °C for 30 min, and the supernatant was removed. The residuals were hydrolyzed by 9.2 M HClO4, and an anthrone reagent was used to measure starch . Fresh samples (the second part) of roots and young leaves were ground into powder to determine NH4+ and NO3− . Fresh roots or young leaves (0.5 g) were extracted using 95% ethanol for 2 h at 60 °C. A folin reagent and 20% Na2CO3 (g/L) were added into the extract for 0.5 h at 30 °C and then tested for total phenols at 760 nm.
To compare males and females, a transcriptome analysis of young leaves was performed between control- and limited-nitrogen treatments (N1 and N3). Three replicates per treatment were analyzed. Total RNA was extracted using Plant RNA Purification Reagent, and genomic DNA was removed using DNase I (TaKara). A detailed description of the RNA quality determination, library preparation, and Illumina Hiseq sequencing are available in Methods S2. The clean reads were separately aligned using the reference genome Populus_trichocarpa (http://plants.ensembl.org/Populus_trichocarpa/Info/Index). RSEM was used to quantify gene abundances (http://deweylab.biostat.wisc.edu/rsem/). The expression level of each transcript was calculated according to the transcript per million reads method. Heatmaps of the expression profiles were generated using the heatmap package in R. A differential expression analysis was performed to identify significantly differently expressed genes (DEGs) using DESeq2 (P-adjust value < 0.05, |log2FC|>1) . A functional-enrichment analysis (KEGG) was performed to identify which DEGs were significantly enriched at the Bonferroni-corrected P-value <0.05 compared with the whole-transcriptome background.
The effects of sex and N treatment on soil chemicals, plant physiological traits, and Shannon diversity were analyzed by a two-way ANOVA analysis. The beta diversity of fungi and bacteria was assessed using nonmetric multidimensional scaling (NMDS) based on Bray-Curtis distances. The significance of sex, N treatment, and their interaction effects on community dissimilarity was tested with PERMANOVA by using the “adonis” function in vegan package in R . A hierarchical clustering analysis was performed by using the method “Ward.D2” with “hclust ()” function in R. The Venn diagram analysis was used to define bacterial genera and fungal OTUs shared or not shared among rhizosphere soil, roots, old leaves, and young leaves and visualized in Cytoscape (v. 3.5). Interactions and connections of physiological traits with the alpha diversity and dominant bacterial/fungal taxa were assessed by Mantel test with vegan package in R . The co-occurrence networks were visualized in Gephi and performed by Conet in Cytoscape (v. 3.5) to show the network topological attributes of bacterial and fungal communities . Only significant (P < 0.01) and robust Spearman’s correlations (r > 0.6 or r < 0.6) were kept in the networks. The linear discriminant analysis effect size (LEfSe) with logarithmic LDA> 3 (Wilcoxon P < 0.05) was applied to identify biomarker taxa of different niches to show the effects of sex. In each niche, the EdgeR’ was used to reveal different abundances of OTUs. OTUs present in at least 60% of samples and with a relative abundance ≥ 0.1% were defined as dominant taxa. The phylogenetic tree of roots and young leaves was annotated and visualized in iTOL software . The dominant taxa, biomarker taxa (LEfSe), and hub nodes among networks were defined as keystone taxa.
The differences of bacterial communities in the rhizosphere soil, roots, and young leaves were significantly affected by soil N (Fig. 1). However, sex imposes stronger effects on bacterial communities than soil N levels with clear differences between the roots, old leaves, and young leaves of males and females. The beta diversity of the bacterial community of old leaves was significantly affected only by sex (Fig. 1c). Similarly, sex was the main factor in driving differences of fungal communities in roots, old leaves, and young leaves (Fig. 1f–h).
The hierarchical clustering analysis revealed three clusters in bacterial communities: rhizosphere soil, roots, and leaves (old and young) were well separated at the genus level (Fig. 2a). The fungal communities of the rhizosphere soil, roots, and leaves (old and young) were well separated, as also the fungal communities of male and female roots (Fig. 2b). The relative abundance of Acidobacteria and Chloroflexi declined in roots, old leaves, and young leaves relative to the rhizosphere soil (Fig. 2c). Unclassified genera belonging to Ascomycota were enriched in roots, old leaves, and young leaves (Fig. 2c). The relative abundance of Actinobacteria and Geopora was significantly higher in roots than in other parts in different treatments (Fig. S1). The alpha diversity of bacterial and fungal communities was mainly affected by the niche, the rhizosphere soil showing the highest diversity (Fig. 2d). The alpha diversity of fungal communities was significantly impacted by sex as well (Fig. 2d, Fig. S2b).
It was found that the four niches harbored 280 bacterial genera and 22 fungal OTUs in common, the majority of them belonging to Proteobacteria and Ascomycota (Fig. 3). The rhizosphere soil, roots, old leaves, and young leaves had 498, 63, 59, and 104 unique bacterial genera, and 257, 3, 1, and 25 unique fungal OTUs, respectively (Fig. 3). The numbers of bacterial genera belonging to Bacteroidetes and Firmicutes were much higher in old and young leaves than in roots and soil (Fig. 3a). The differential abundance analysis indicated that more fungal OTUs were absent in the endophytes of females compared to males (Fig. S3).
The LDA effect size of LEfSe clearly indicated the presence of bacterial and fungal biomarker taxa in different niches (Fig. 4). The most significant bacterial biomarker taxon in the rhizosphere soil of P. cathayana females was Actinobacteria. The most significant bacterial and fungal biomarker taxa in the root endosphere of females were Actinobacteria and the genus Geopora, respectively, while in the root endosphere of males they were Proteobacteria and a genus belonging to Ascomycota, respectively. Other significant bacterial and fungal biomarker taxa were found in the old leaf endosphere. Young leaves of P. cathayana males and females also harbored distinct bacterial and fungal biomarker taxa.
To further characterize the niche and sex effect on plant microbiomes, particularly on plant endophytes, the co-occurrence patterns of bacterial and fungal communities were assessed along the soil-plant continuum. Bacterial communities of the rhizosphere soil and young leaves showed a higher network complexity than those of roots or old leaves with a higher average degree and clustering coefficient (Fig. S4a; Table S1). However, fungal communities of young leaves showed the highest network complexity with the highest average degree (40.25) and clustering coefficient (0.75) (Fig. S4b; Table S1). The bacterial taxonomic composition of “hub nodes” in the network differed between the rhizosphere soil and plant compartments (Fig. 5a). Positive edges and the average degree of bacterial networks were higher in the rhizosphere soil, roots, old leaves, and young leaves of females when compared to males, while positive edges and the average degree of fungal communities in young leaves of P. cathayana females were much higher than those of other niches (Table 1). OTUs belonging to Bacteroidetes and Firmicutes in young leaves of P. cathayana males, and OTUs belonging to Proteobacteria and Firmicutes in young leaves of females were more connected with others (Fig. 5b). OTUs belonging to Ascomycota in young female leaves were more connected with other OTUs in the network (Fig. 5b).
In microbial communities of the endosphere, the dominant bacterial OTUs mainly belonged to Gammaproteobacteria, Alphaproteobacteria, and Actinobacteria at the class level, while dominant fungal OTUs were affiliated with Ascomycota at the genus level (Fig. 6, Fig. S5). A total of 114 and 20 dominant OTUs of bacteria and fungi, respectively, were identified. P. cathayana males and females harbored 30 and 25 distinct dominant bacterial OTUs, while 5 and 4 unique and dominant fungal OTUs, respectively. Unique dominant bacterial OTU1204, OTU1621, and OTU 825 in the female root endosphere showed a higher relative abundance than other root endophytes (Fig. 6a). The shared OTU4235, OTU4835, and OTU2362 between males and females had higher relative abundances than other leaf (both old and young) endophytes. Dominant fungal OTU229 and OTU125 with much higher relative abundances than other dominant OTUs were shared among roots, old leaves, and young leaves in both males and females (Fig. S5).
Sex significantly affected soil total nitrogen, total phosphorus, NH4+, NO3−, microbial carbon, microbial nitrogen, and NAG (Table S2). Nitrogen and sex showed significant effects on plant growth (Table S3). Concentrations of phenols, tannins, NH4+ and NO3− of young leaves, concentrations of phenols, tannin, and starch, and the pH level of roots were significantly affected by sex (Fig. 7). The soil N treatment also showed significant effects on tannins, pH, NH4+ and NO3− in young leaves, and on pH and of NO3− in roots (Fig. 7). Young male and female leaves showed great differences in gene expression, females being more sensitive to the N treatment, as revealed by transcriptomic data (Fig. S6). A total of 37,008 genes were identified. The gene expression analysis revealed 272 genes that were significantly upregulated in males, and 275 genes that were significantly downregulated in females in the control treatment. The N-limited treatment caused much greater differences in gene expression between males and females (Fig. S6b). The KEGG pathway enrichment showed that more genes related to plant-pathogen interaction functions were significantly different between males and females in control and N-limited conditions (Fig. S6c).
The pH value was negatively related with tannins in young leaves and roots (Fig. 8). Starch and pH of young leaves had significant effects on the abundance of Proteobacteria, while tannins and pH of roots showed significant effects on the abundance of Chloroflexi, Actinobacteria, and Proteobacteria, and bacterial Shannon diversity (Fig. 8a, b). The abundance of unclassified genera belonging to Ascomycota was significantly affected by starch, while the fungal Shannon diversity was affected by tannins in young leaves (Fig. 8c). The abundance of Geopora was affected by starch, while the fungal Shannon diversity was affected by the NH4+ concentration of roots (Fig. 8d).
Our initial hypothesis was that sex is an important factor influencing the properties (including microbial communities) and functions of the rhizosphere soil. To adapt to varied environments, different plant species or genotypes are able to recruit beneficial microbes and vary interactions among microbes in the rhizosphere soil through changing chemical conditions and by releasing root-derived compounds [13, 22, 38]. The contents of total nitrogen, total phosphorus, NH4+, and NO3− were lower in the rhizosphere soil of females compared to males, probably because of higher nutrient demands of P. cathayana females for growth (Tables S2 and S3).
Different chemical traits of the rhizosphere soil have been reported as being important factors influencing bacterial and fungal communities and functions in P. euphratica males and females in natural forests . Another important factor is the microbial substrate preferences to different root exudates , while the quality and amount of photosynthesis-derived exudates in P. cathayana males and females are likely to be different because of differences in carbon fixation and allocation, nutrient uptake, and utilization, and connections with arbuscular mycorrhizal fungi [7, 8, 26]. In this study, sex greatly altered interactions among hub nodes of different OTUs as indicated by the co-occurrence network, the biomarker taxa, and the microbial carbon biomass of the rhizosphere soil (Table 1, S2; Figs. 4 and 5). Microbial enzyme activities dynamically respond to alterations in the nutrient and substrate availability and they represent the biological metabolism of the soil microbiome . Sex significantly affected the NAG activity related to chitin degradation (Table S2), which reflected variations in soil properties and microbial communities imposed by P. cathayana males and females.
Differences have been reported in poplars [9, 40, 41] and other species [10, 11, 23] in the plant microbiome among soil, root, and other microhabitats. This study also demonstrated a strong niche differentiation, such as a large number of unique taxa in the rhizosphere soil, a decreased alpha diversity, and varied abundance of major taxa along the soil-P. cathayana continuum. Extreme pH conditions impose stronger selective pressures for the survival and fitness of microbes, which suggests a stochastic process in neutral conditions but a deterministic assembly in acidic or alkaline pH conditions . This study revealed a clear decreasing pH gradient from neutral (rhizosphere soil) to acidic conditions (leaf) (Table S1, Fig. 7), which demonstrated an ecological filter during community assembly processes along the soil-root-leaf continuum. For example, the rapid decrease of Basidiomycota within plant compartments, particularly in young leaves, was likely to be related with the decreasing pH because most Basidiomycota taxa perform better in higher pH conditions . Fine-tuning and adaptation of the microbial assembly and stability, e.g., in roots and leaves, are clearly evident [9, 11, 23]. Besides different pH conditions, phenol and tannin levels of young leaves were higher than those of roots, which indicated that the microbes needed to have an ability to deal with defense chemicals. However, the sex impacts on the endophyte composition and abundance in each niche were not amplified by the great differences in the transcriptome and physiological traits between males and females under nitrogen-limited conditions (Fig. S1). Particularly, the fungal community inhabiting the plant endosphere stayed stable under different nitrogen treatments, which demonstrated a stronger niche-depended selection instead of external environmental effects.
A successful colonization of bacteria and fungi in roots and leaves requires that these microbes have the capacity to overcome the host’s immune system in different microhabitats [9,10,11]. Aerial plant parts impose much stronger selective pressures than belowground parts on the microbiome assembly [11, 44]. Indeed, the composition and structure of microbes residing in leaves (both old and young) were very different compared to those belowground, while young and old leaves had a similar microbe composition, relative abundance, and alpha and beta diversity (Figs. 2 and 6, S1, S5). Fungal taxa from the phylum Ascomycota enrich in live and senescent leaves of trees . For most microbes in this study, young leaves did not impose a stronger inhibitory influence on the colonization or growth of the endophytes than old leaves. It was evident that these bacterial and fungal microbes colonize leaves shortly after emergence and reach and keep a constant level before leaf abscission. However, some fungal taxa with a relatively low abundance (classified as “others” in this study) successfully colonized old leaves (Fig. S1b), probably as a consequence of the declined action of protective mechanisms resulting in a rapid proliferation of opportunistic fungi during leaf senescence [45, 46].
One interesting finding was that the network hubs in young leaves were more connected with other taxa in the fungal community than those in old leaves (Table 1; Fig. 5). Keystone species belonging to Basidiomycota were more connected with other taxa in young leaves than in old leaves (Fig. 5). Positive edges in the fungal network greatly declined from young to old leaves, implying substantially decreased competition among fungal species [11, 47]. A loosely assembled fungal network reflects subtle changes in the host status; old leaves are suggested to decompose under saprophytic fungi from Ascomycota before abscission [45, 48].
Our results further confirmed that the selective pressure imposed by sex on endophytes in different plant parts is much stronger than that on microbial communities of the rhizosphere soil. The community structure, biomarker/core taxa, and interactions among bacterial and fungal taxa in roots, old leaves, and young leaves were different between P. cathayana males and females. Key microbial taxa typically enriched in different parts of P. cathayana males and females, such as fast-growing Alphaproteobacteria and Gammaproteobacteria (Fig. 6), reflecting adaptation of these plant-associated microbes to interior environments of plants by acquiring carbohydrates and nutrients from the host . For example, the Gammaproteobacteria species are capable of colonizing and then dominating diverse niches, including rhizosphere, and root and leaf endosphere, with crucial roles in regulating host fitness, pathogen suppression, and tolerance to varied stresses [49, 50]. However, different relative abundances of these taxa and distinct taxa in each plant compartment indicated sex specificity in the structure and function of entophytic communities (Fig. 6). Streptomycetes and Geopora were biomarker taxa and enriched in the female root endosphere (Fig. 4). Streptomycetes species colonizing root systems are considered beneficial due to their capacity to promote plant growth and induce defense responses against plant pathogens [51, 52]. Geopora had a higher relative abundance in the root endosphere than in the rhizosphere soil (Fig. 2), and the enrichment was greater in female roots (Fig. S1b).
Based on the limited knowledge on how male and female plants, especially in Populus species, resist biotic stresses, the enriched Streptomycetes species have the potential ability to enhance pathogen resistance in females. The species of the genus Geopora are important mycobionts of ectomycorrhizal fungi [53, 54]. Female roots probably establish a closer mutualistic relationship with Geopora species than male roots, which contributes to meeting the higher nitrogen demand of P. cathayana females during growth . Alternaria species are mainly saprophytic fungi . In this study, they were identified as biomarker taxa by LEfSe and also as dominant taxa in old leaves of P. cathayana females, indicating sex-specific differences in recruiting microbes related with the decomposition process before leaf abscission.
This study further found distinct bacterial and fungal biomarker and dominant taxa in young leaves of P. cathayana males and females (Fig. 4, S5). Interestingly, the fungal co-occurrence network was much more complex and the alpha diversity of young female leaves was significantly higher compared to male leaves (Fig. 5b, S2). The more connected network among OTUs from Ascomycota, Basidiomycota, Zoopagomycota, and other taxa suggested a closer relationship between female leaves and its associated microbes, which implied that the fitness of the female host was more dependent on the fungal community compared to males. Hub species in microbial networks are recognized as mediators between the host and its associated microbiome. Through network hubs, a host can selectively influence the assembly of their associated microbiome by regulating microbe-microbe interactions . Variations in the constituting microbial taxa and in the number of hub microorganisms have a significant influence on the assembly of microbial communities, and this effect is suggested to be independent from external environmental factors .
Furthermore, the N levels in this study clearly caused sexually different gene expressions in various functions by transcriptomic data (Fig. S6). Previous studies also reported greatly differences in gene expressions and physiological traits at similar N gradient [20, 26, 27]. In the present study, substantial differences between males and females in gene expression related to plant-pathogen interactions were found (Fig. S6). However, these differences did not lead to great changes in the composition and structure of leaf endophyte (Figs. S1 and S2). The composition and abundance of endophytic communities in the re-irrigated nitrogen treatment were similar to those of other two N treatments. Instead, there was convincing evidence that a stronger control by a host rather than by the external N dose shapes specific endophytic communities, probably due to different traits in males and females. Plant-associated microbes differ in their ability to successfully colonize different plant parts due to their differences in overcoming cell wall, in resisting defense-related chemicals and in carbohydrate metabolic functions of the host [25, 57, 58]. Leaves of P. cathayana males have a higher mesophyll palisade tissue thickness but a lower spongy cell density , while males have a smaller root system compared to females . These differences in the structure are probably important forces in selecting endophytes inhabiting leaves or roots in P. cathayana males and females.
Physiological traits of roots and young leaves were used to reflect the environment of the plant endosphere, and most of these traits were significantly affected by sex (Fig. 7). Leaf carbon is suggested being the main driving force for changes in the diversity, richness, and composition of foliar fungal endophytes , while many studies have emphasized that secondary metabolites of plants, particularly those related to plant defense traits, influence the plant microbiome . However, after exploring twelve P. trichocarpa genotypes, which produced different amounts of salicylic acid-related secondary metabolites, Veach et al.  suggested that plant defense strategies partially drive colonization and the assembly of taxon-specific microorganisms. In addition, increased primary metabolites facilitate the colonization of suitable microbes within roots . Our results were unable to explain the specific biological mechanisms that control a host plant’s microbiome. However, this study demonstrated that the microbial assembly was affected or driven by diverse conditions within plant tissues, such as pH, carbon supply (starch), nitrogen status (NH4+), or defense traits (total phenols and tannins) (Fig. 8). Previous empirical evidence has showed that plant hosts exert strong effects on plant microbiomes by deterministic selection during plant growth and development . We found that P. cathayana males and females have sex-specific selection pressures or specialized endosphere habitats that microbes select to form distinct microbial compositions.
In this study, bacterial and fungal community assemblies across a plant-soil continuum were mainly influenced by niches, which suggested that ecological filters, for example, a decreasing pH gradient from the rhizosphere soil to leaves, affect the community assembly processes along the soil-root-leaf continuum. Endophytic fungal communities showed a higher stability than bacterial communities in each niche in response to different soil N treatments, which demonstrated the presence of a powerful control of a host in selecting fungi. Our results suggested potentially different functions in N-cycling, nutrient acquisition, and pathogen resistance due to divergent microbial communities in the rhizosphere soil, roots, and leaves. All results provided important new knowledge for understanding the effects of sexual dimorphism in microbial assemblies and plant-microbiome interactions.
Availability of data and materials
All bacterial and fungal databases in this study are available at National Center for Biotechnology Information under BioProject ID PRJNA797907.
Nonmetric multidimensional scaling
Total nitrogen concentration
Available phosphorus concentration
Total phosphorus concentration
Linear discriminant analysis effect size
Renner SS. The relative and absolute frequencies of angiosperm sexual systems: dioecy, monoecy, gynodioecy, and an updated online database. Am J Bot. 2014;101:1588–96.
Xue L, Wu H, Chen Y, Li X, Yin T. Evidences for a role of two Y-specific genes in sex determination in Populus deltoides. Nat Commun. 2020;11:5893.
Bertiller MB, Sain CL, Bisigato AJ, Coronato FR, Aries JO, Graff P. Spatial sex segregation in the dioecious grass Poa ligularis in northern Patagonia: the role of environmental patchiness. Biodivers Conserv. 2002;11:69–84.
Hultine KR, Grady KC, Wood TE, Shuster SM, Stella JC, Whitham TG. Climate change perils for dioecious plant species. Nat Plants. 2016;2:16109.
Liu M, Korpelainen H, Li CY. Sexual differences and sex ratios of dioecious plants under stressful environments. J Plant Ecol. 2021;14:920–33.
Avila-Sakar G, Romanow CA. Divergence in defence against herbivores between males and females of dioecious plant species. Int J Evol Biol. 2012;23:897157.
Xia ZC, He Y, Yu L, Lv RB, Korpelainen H, Li CY. Sex-specific strategies of phosphorus (P) acquisition in Populus cathayana as affected by soil P availability and distribution. New Phytol. 2020;225:782–92.
Wu XY, Liu JT, Meng QQ, Fang SY, Kang JY, Guo QX. Differences in carbon and nitrogen metabolism between male and female Populus cathayana in response to deficient nitrogen. Tree Physiol. 2021;41:119–33.
Beckers B, Beeck M, Weyens N, Boerjan W, Vangronsveld J. Structural variability and niche differentiation in the rhizosphere and endosphere bacterial microbiome of field-grown poplar trees. Microbiome. 2017;5:25.
Bernard J, Wall CB, Costantini MS, Rollins RL, Amend AS. Plant part and a steep environmental gradient predict plant microbial composition in a tropical watershed. ISME J. 2021;15:999–1009.
Xiong C, Singh BK, He JZ, Han YL, Li PP, Wan LH, et al. Plant developmental stage drives the differentiation in ecological role of the maize microbiome. Microbiome. 2021;9:171.
Trivedi P, Leach JE, Tringe SG, Sa TM, Singh BK. Plant-microbiome interactions: from community assembly to plant health. Nat Rev Microbiol. 2020;18:607–21.
Guo QX, Yan LJ, Korpelainen H, Niinemets Ü, Li CY. Plant-plant interactions and N fertilization shape soil bacterial and fungal communities. Soil Biol Biochem. 2019;127:127–38.
Vega-Frutis R, Munguía-Rosas MA, Varga S, Kytöviita MM. Sexspecific patterns of antagonistic and mutualistic biotic interactions in dioecious and gynodioecious plants. Perspec Plant Ecol Evol Syst. 2013;15:45–55.
Wu QP, Tang Y, Dong TF, Liao YM, Li DD, He XH, et al. Additional AM fungi inoculation increase Populus cathayana intersexual competition. Front Plant Sci. 2018;9:607.
Stringlis IA, Yu K, Feussner K, Jonge RD, Bentum SV. MYB72-dependent coumarin exudation shapes root microbiome assembly to promote plant health. Proc Natl Acad Sci USA. 2018;115:5213–22.
He Y, Zhu ZD, Guo QX, Xia ZC. Sex-specific interactions affect foliar defense compound accumulation and resistance to herbivores in Populus cathayana. Sci Total Environ. 2021;774:145819.
Randriamanana TR, Nybakken L, Lavola A, Aphalo PJ, Nissinen K, Julkunen-Tiitto R. Sex-related differences in growth and carbon allocation to defence in Populus tremula as explained by current plant defense theories. Tree Physiol. 2014;34:471–87.
Vega-Frutis R, Luna-Esquivel G, Figueroa-Esquivel EM. Land-use change impact on mycorrhizal symbiosis in female and male plants of wild Carica papaya (Caricaceae). Symbiosis. 2018;76:1–11.
Song HF, Cai ZY, Liao J, Tang DT, Zhang S. Sexually differential gene expressions in poplar roots in response to nitrogen deficiency. Tree Physiol. 2019;39:1614–29.
Mushinski RM, Payne ZC, Raff JD, Craig ME, Pusede SE, Rusch DB, et al. Nitrogen cycling microbiomes are structured by plant mycorrhizal associations with consequences for nitrogen oxide fluxes in forests. Global Change Biol. 2020;27:1068–82.
Guo QX, Liu JT, Yu L, Korpelainen H, Li CY. Different sexual impacts of dioecious Populus euphratica on microbial communities and nitrogen cycle processes in natural forests. Forest Ecol Manage. 2021;496:119403.
Hou SJ, Thiergart T, Vannier N, Mesny F, Ziegler J, Pickel B, et al. A microbiota-root-shoot circuit favours Arabidopsis growth over defence under suboptimal light. Nat Plants. 2021;7:1078–92.
Lebeis SL, Paredes SH, Lundberg DS, Breakfield N, Gehring J, McDonald M, et al. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science. 2015;349:860–4.
Cole BJ, Feltcher ME, Waters RJ, Wetmore KM, Visel A. Genome-wide identification of bacterial plant colonization genes. PLoS Biol. 2018;15:e2002860.
Zhang S, Jiang H, Zhao HX, Korpelainen H, Li CY. Sexually different physiological responses of Populus cathayana to nitrogen and phosphorus deficiencies. Tree Physiol. 2014;34:343–54.
Luo J, Zhou J, Li H, Shi W, Polle A, Lu M, et al. Global poplar root and leaf transcriptomes reveal links between growth and stress responses under nitrogen starvation and excess. Tree Physiol. 2015;35:1283–302.
Schindlbacher A, Rodler A, Kuffner M, Kitzler B, Sessitsch A, Zechmeister-Boltenstern S. Experimental warming effects on the microbial community of a temperate mountain forest soil. Soil Biol Biochem. 2011;43:1417–25.
Caporaso JG, Lauber CL, Walters WA, Berglyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108:4516–22.
White TJ, Bruns TD, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfaud DH, Sninsky JJ, White TJ, editors. PCR protocols: a guide to methods and applications. London: Academic; 1990. p. 315–22.
Samad A, Trognitz F, Compant S, Antonielli L, Sessitsch A. Shared and host-specific microbiome diversity and functioning of grapevine and accompanying weed plants. Environ Microbiol. 2017;19:1407–24.
Wang F, Men X, Zhang G, Liang KC, Xin YH, Wang J, et al. Assessment of 16S rRNA gene primers for studying bacterial community structure and function of aging flue-cured tobaccos. AMB Expr. 2018;8:182.
Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.
Cornelissen J, Quested H, van Logtestijn R, Pérez-Harguindeguy M, Gwynn-Jones D, Díaz S, et al. Foliar pH as a new plant trait: can it explain variation in foliar chemistry and carbon cycling processes among subarctic plant species and types? Oecologia. 2006;147:315–26.
Oksanen J, Kindt R, Legendre P, O’Hara B, Stevens MHH, Oksanen MJ, et al. The vegan package. Community Ecol Package. 2007;10:631–7.
Sunagawa S, Coelho LP, Chaffron S, et al. Structure and function of the global ocean microbiome. Science. 2015;348:1261359.
Zhalnina K, Louie KB, Hao Z, Mansoori N, da Rocha UN, Shi SJ, et al. Dynamic root exudate chemistry and microbial substrate preferences drive patterns in rhizosphere microbial community assembly. Nat Microbiol. 2018;3:470–80.
Chen Y, Liu X, Hou YH, Zhou SR, Zhu B. Particulate organic carbon is more vulnerable to nitrogen addition than mineral-associated organic carbon in soil of an alpine meadow. Plant Soil. 2021;458:93–103.
Cregger MA, Veach AM, Yang ZK, Crouch MJ, Vilgalys R, Tuskan AG. The Populus holobiont: dissecting the effects of plant niches and genotype on the microbiome. Microbiome. 2018;6:31.
Bonito G, Reynolds H, Robeson MS, Nelson J, Hodkinson BP, Tuskan G, et al. Plant host and soil origin influence fungal and bacterial assemblages in the roots of woody plants. Mol Ecol. 2014;23:3356–70.
Tripathi BM, Stegen JC, Kim M, Dong K, Adams JM, Lee YK. Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. ISME J. 2018;12:1072–83.
Xu F, Cai TJ, Yang X, Sui WZ. Soil fungal community variation by large-scale reclamation in Sanjiang plain, China. Ann Microbiol. 2017;67:679–89.
Wagner M, Lundberg D, del Rio T, Tringe SG, Dangl JL, Mitchell-Olds T. Host genotype and age shape the leaf and root microbiomes of a wild perennial plant. Nat Commun. 2016;7:12151.
Voříšková J, Baldrian P. Fungal community on decomposing leaf litter undergoes rapid successional changes. ISME J. 2013;7:477–86.
Pang ZQ, Chen J, Wang TH, Gao CS, Li ZM, Guo LT, et al. Linking plant secondary metabolites and plant microbiomes: a review. Front Plant Sci. 2021;12:621276.
van der Heijden MGA, Hartmann M. Networking in the plant microbiome. PLoS Biol. 2016;14:e1002378.
Guerreiro MA, Brachmann A, Begerow D, Peroh D. Transient leaf endophytes are the most active fungi in 1-year-old beech leaf litter. Fungal Divers. 2017;89:1–15.
Mendes R, Kruijt M, de Bruijn I, Dekkers E, van der Voort M, Schneider JH, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science. 2011;332:1097–100.
Alvarez-Perez JM, Gonzalez-Garcia S, Cobos R, Olego MA, Ibanez A, Diez-Galan A, et al. Use of endophytic and rhizosphere Actinobacteria from grapevine plants to reduce nursery fungal graft infections that lead to young grapevine decline. Appl Environ Microb. 2017;83:e01564–17.
Salla TD, Thanise R, Astarita LV, Santarém ER. Streptomyces rhizobacteria modulate the secondary metabolism of Eucalyptus plants. Plant Physiol Biochem. 2014;85:14–20.
Tian X, Cao L, Tan H, Han W, Chen M, Liu Y, et al. Diversity of cultivated and uncultivated actinobacterial endophytes in the stems and roots of rice. Microb Ecol. 2007;53:700–7.
Patterson A, Flores-Rentería L, Whipple A, Whitham T, Gehring C. Common garden experiments disentangle plant genetic and environmental contributions to ectomycorrhizal fungal community structure. New Phytol. 2019;221:493–502.
Shemesh H, Boaz BE, Millar CI, Bruns TD. Symbiotic interactions above treeline of long-lived pines: Mycorrhizal advantage of limber pine (Pinus flexilis) over Great Basin bristlecone pine (Pinus longaeva) at the seedling stage. J Ecol. 2019;108:908–16.
Jiang H, Shi YT, Zhou ZX, Yang C, Chen YJ, Chen LM, et al. Leaf chemistry and co-occurring species interactions affecting the endophytic fungal composition of Eupatorium adenophorum. Ann Microbiol. 2011;61:655–62.
Agler MT, Ruhe J, Kroll S, Morhenn C, Kim ST, Weigel D, et al. Microbial hub taxa link host and abiotic factors to plant microbiome variation. PLoS Biol. 2016;14:e1002352.
Levy A, Salas Gonzalez I, Mittelviefhaus M, Clingenpeel S, Paredes SH, Miao JM, et al. Genomic features of bacterial adaptation to plants. Nat Genet. 2018;50:138–50.
Veach AM, Morris R, Yip DZ, Yang ZK, Engle NL, Cregger MA, et al. Rhizosphere microbiomes diverge among Populus trichocarpa plant-host genotypes and chemotypes, but it depends on soil origin. Microbiome. 2019;7:1–15.
Liu M, Liu XC, Du XH, Korpelainen H, Niinemets U, Li CY. Anatomical variation of mesophyll conductance due to salt stress in Populus cathayana females and males growing under different inorganic nitrogen sources. Tree Physiol. 2021;41:1462–78.
Yang T, Weisenhorn P, Gilbert JA, Ni Y, Sun R, Shi Y, et al. Carbon constrains fungal endophyte assemblages along the timberline. Environ Microbiol. 2016;18:2455–69.
Mangeot-Peter L, Tschaplinski TJ, Engle NL, Veneault-Fourrey C, Martin F, Deveau A. Impacts of soil microbiome variations on root colonization by fungi and bacteria and on the metabolome of Populus tremula × alba. Phytobiomes J. 2020;4:142–55.
We acknowledge the help of Dr. Lin Zhang with data analysis, and of Ms Xiaoyi Wu with experimental work.
This research was supported by the Natural Science Foundation of China (31800508).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary methods. Table S1. Bacterial and fungal co-occurrence network characteristics in each niche. Table S2. The physiological traits and microbial functions of rhizosphere soil affected by sex and nitrogen treatment. Table S3. The growth characteristics of P. cathayana males and females under different N treatments. Figure S1. The composition and relative abundance of bacterial (a) and fungal (b) communities in the rhizosphere soil, roots, old leaves and young leaves of P. cathayana males and females in different N treatments. Figure S2. The Shannon diversity of bacterial and fungal communities in different N treatments. From left to right: rhizosphere soil, roots, old leaves and young leaves. Figure S3. The volcano plot illustrating the enrichment and depletion patterns of bacterial and fungal microbiomes in each niche (male vs female). Figure S4. Different niche effects on the networks of bacterial (a) and fungal (b) communities. Figure S5. Phylogenetic tree, taxonomic composition, and distribution patterns of dominant fungal taxa living in roots, old leaves and young leaves. (a) Dominant fungal taxa in P. cathayana females and males, respectively. (b) The distribution pattern of dominant fungal taxa in P. cathayana females. (c) The number of dominant fungal taxa, the shared and specific bacterial OTUs in P. cathayana females and males. (d) The distribution pattern of dominant fungal taxa in P. cathayana males. Figure S6. Heatmaps of the expression profiles (a), differential expression analysis (DEGs) (b) and functional-enrichment analysis (KEGG) (c). MC, young male leaves in control-nitrogen treatment; MN, young male leaves in limited-nitrogen treatment; FC, young female leaves in control-nitrogen treatment; FN, young female leaves in limited-nitrogen treatment.
About this article
Cite this article
Guo, Q., Liu, L., Liu, J. et al. Plant sex affects plant-microbiome assemblies of dioecious Populus cathayana trees under different soil nitrogen conditions. Microbiome 10, 191 (2022). https://doi.org/10.1186/s40168-022-01387-9
- Dioecious plants
- Plant-soil continuum
- Plant microbiome
- Plant-microbe interactions