Abstract
Introduction:
Microbiome within plant tissues is pivotal for co-evolution with host plants. This microbiome can colonize the plant, with potential transmission via seeds between parents and offspring, affecting seedling growth and host plant adaptability to the environment.
Methods:
We employed 16S rRNA gene amplicon analysis to investigate the vertical distribution of core microbiome in cotton seeds across ecological niches [rhizosphere, root, stem, leaf, seed and seed-P (parental seed)] of the three cotton genotypes.
Results:
The findings demonstrated a significant decrease in microbiome diversity and network complexity from roots, stems, and leaves to seeds. The microenvironment exerted a more substantial influence on the microbiome structure of cotton than the genotypes. The core endophytic microorganisms in cotton seeds comprised 29 amplicon sequence variants (ASVs) affiliated with Acidimicrobiia, Alphaproteobacteria, Bacilli, Bacteroidia, Clostridia, Gammaproteobacteria, and unclassified_Proteobacteria. These vertically transmitted taxa are widely distributed in cotton plants. Through 16S rRNA gene-based function prediction analysis of the cotton microbiome, we preliminarily understood that there are potential differences in metabolic capabilities and phenotypic traits among microbiomes in different microhabitats.
Discussion:
In conclusion, this study demonstrated the crucial role of the microenvironment in influencing the cotton microbiome and offered insights into the structures and functions of the cotton seed microbiome, facilitating future crop yield enhancement through core seed microbiome regulation.
Introduction
The plant microbiome plays a pivotal role in plant health and forms a large biological system known as the holobiome (Vandenkoornhuyse et al., 2015), interacting with plants at various physiological, ecological, and evolutionary levels. This interaction commences at the seed stage, where microbiomes, particularly seed-endophytic microorganisms, significantly influence seedling growth vigor and the environmental adaptation of the plant (Frank et al., 2017). Seed-endophytic bacteria engage in intricate ecological relationships with plant hosts, including symbiosis, co-evolution, and benign parasitism. Importantly, plants selectively transmit beneficial taxa from one generation to the next (Ferreira et al., 2008; Hardoim et al., 2011; Truyens et al., 2014; Sánchez-López et al., 2018).
Situated in a warm temperate and inland region, Xinjiang is renowned for its ample heat and sunlight, conditions that are particularly conducive to cotton cultivation. Cotton (Gossypium spp.) holds a pivotal role in global agriculture, not merely as a strategic material integral to national economies but also as a significant contributor to economic benefits. Its wide-ranging applications span textiles, oil production, and the medical industry (Predeepa et al., 2022; Zhang et al., 2023). Leveraging the unique geographic advantages of Xinjiang, conducive to optimal cotton growth, this region accounts for over 80% of China’s cotton yield, thereby constituting nearly a quarter of global production (Ju et al., 2021). Cotton is vital for regional economic growth and the preservation of ecological balance (Adham et al., 2020). Nonetheless, the sustainability of cotton production is increasingly challenged by climate change, pest infestations, and various environmental stressors (Nelson, 2004). Consequently, rigorous research into the cotton seed microbiome is imperative for augmenting the crop’s resilience against such adversities.
The current study focuses on fundamental research into the cotton microbiome, particularly concerning rhizosphere and endophytic fungi (Nelson, 2018; Da Silveira et al., 2019; Hone, 2021; Bai, 2022; Hou, 2023), yet research on the cotton endophytic bacterial microbiome remains limited. Studies of plant microbiome have shown that research on rice microbiome reveals the origins of seed endophytes, the composition of core microbiome, and the interaction patterns between microorganisms and host plants (Walitang et al., 2019; Zhao, 2019; Zhang et al., 2022), providing an important reference framework for our study on the formation and function of cotton seed microbiome. Understanding the composition and mechanisms of the seed microbiome allows us to better infer and analyze potential characteristics of the cotton seed microbiome, particularly in response to various stresses. Recognizing the importance of priority effects, optimizing the seed microbiome methods holds significant potential for plant breeding and crop improvement (Mitter et al., 2017; Chesneau et al., 2020). Identifying the core plant microbiome serves as a valuable resource in enhancing defenses against environmental stresses and plant diseases (Gopal et al., 2013). Therefore, this study aims to fill the research gap in the cotton seed microbiome, which is crucial for improving cotton’s adaptability and productivity.
This study aimed to achieve the following specific objectives: (1) investigate the structure and functionality of microbiome associated with cotton, (2) identify the core taxa within cotton seed microbiome, and (3) examine the vertical diffusion of seed microorganisms and their distribution across various plant parts throughout the growth stages of cotton.
Materials and methods
Experimental setup and sampling
In Shihezi City, Xinjiang, located at 86°3’52.884” E and 44°18’36” N, we cultivated three cotton genotypes: R2069, XLZ63, and XLZ78. Complete cotton plants and corresponding rhizosphere samples were collected at the flowering stage, followed by seed collection at harvest. Each of the three distinct genotypes was planted in quadruplicate, with each plot measuring 18 m2 (3 m × 6 m) and managed with a fertilization regimen of 40 kg of nitrogen, 25 kg of phosphate (P2O5), and 20 kg of potash (K2O) per mu. Of the nitrogen fertilizer, 25% was used as base fertilizer, with the remaining 75% applied in stages. For phosphate and potash fertilizers, 70% and 50%, respectively, were used as base fertilizers, with the remainder used as top dressing. Approximately eight rounds of top dressing were conducted throughout the growing period, following the ‘one water, one fertilizer’ management strategy. To comprehensively assess the microbiomes associated with these cotton genotypes, we extracted total DNA from a broad range of samples, comprising 12 rhizosphere samples, 12 root samples, 12 stem samples, 12 leaf samples, 12 seed samples, and 12 seed-P (parental seed) samples, as illustrated in Figure 1A. Each sequenced tissue sample consisted of a composite from four individual plants within the same plot. Using sterile scissors, we cut the roots of the cotton plants and then centrifuged the cuttings at 12,000 rpm for 10 min to collect the rhizosphere, defined as the soil layer within 1–2 mm of the root. Each seed sample was standardized to 10 g. Root, stem, leaf, and seed samples were initially rinsed with tap water to remove the majority of soil and debris. Subsequently, these samples were immersed in 75% ethanol for 1 min for disinfection, followed by a 10-min soak in a 50 g/L sodium hypochlorite solution. Afterward, the samples were thoroughly rinsed three times with sterile water to remove all residual disinfectants. Finally, the surface moisture was gently blotted with sterile paper towels, and the samples were placed in a sterile environment to dry. Subsequently, these samples were individually stored at −80°C for further processing and subsequent analysis.
FIGURE 1
The measurements of organic matter (OM), alkali-hydrolyzable nitrogen (AN), available potassium (AK), available phosphorus (AP), electrical conductivity (EC), and pH in cotton rhizosphere soil were all conducted in accordance with the soil quality guidelines (Hill and Sparling, 2009).
Illumina and 16S rRNA gene sequencing
Total genomic DNA was extracted using a DNeasy PowerSoil Kit (MP Biomedicals, Germany). The extracted DNA was characterized using 1.0% agarose gel electrophoresis and quantified using a NanoDrop One spectrophotometer (Thermo Scientific, USA). The DNA samples were stored at −80°C for future investigations.
In terms of gene library construction, we used the universal primers 335F and 769R to amplify the V3–V4 regions of the 16S rRNA gene, which predominantly amplifies bacterial DNA in plant samples, thereby significantly reducing the likelihood of amplifying plant, archaeal, and other eukaryotic DNA (Dorn-In et al., 2015). Polymerase chain reaction (PCR) was performed using the Phanta Max Master Mix Kit P515 (Vazyme, China) under predefined thermocycling conditions. The procedure commenced with an initial denaturation step at 95°C for 3 min. The amplification comprised seven cycles, each including a denaturation step at 95°C for 45 s, followed by annealing at 65°C for 1 min with a gradual cooling rate of 2°C per cycle, and an extension step at 72°C for 90 s. This was followed by 30 consecutive cycles of stable-state amplification, comprising a 15-s denaturation step at 95°C, a 30-s annealing step at 50°C, and a 30-s extension step at 72°C. The amplification was completed with a final 5-min extension step at 72°C. The amplification products were validated using 1.0% agarose gel electrophoresis. High-throughput sequencing was conducted by BMKCloud Co. using Illumina HiSeq P250 technology, which yielded 250 bp paired-end reads.
In this study, raw data were initially subjected to demultiplexing and quality filtering using Trimmomatic (Bolger et al., 2014) (version 0.33), with quality filtering employing a 50 bp sliding window to trim the trailing bases when the average quality within the window fell below 20. Subsequently, primer sequences were identified and removed with Cutadapt (Martin, 2011) (version 1.9.1) using a maximum mismatch rate of 20% and a minimum overlap of 80%. Paired-end reads were then merged using USEARCH (Edgar, 2016) (version 10) and chimeric sequences were identified and removed with UCHIME (Edgar et al., 2011) (version 8.1). Afterward, the quality-controlled data were denoised using the DADA2 (Callahan et al., 2016; Bolyen et al., 2019) method within QIIME2 (Bolyen et al., 2019) (version 2020.6), setting the threshold to 0.005% of all sequencing reads to filter out ASVs. From a sample set that included 12 rhizospheres, 12 roots, 12 stems, 12 leaves, 12 seeds, and 12 seeds-P, 12,605,376 raw reads were obtained, with an average of 165,637 high-quality reads per sample. After quality trimming, each sequence was classified using the classify-sklearn tool within QIIME2 based on the Naive Bayes classifier (Wang et al., 2007), referencing the SILVA small subunit rRNA database (version 138) with a confidence threshold set at 0.7 (Quast et al., 2012). After excluding ASVs classified as mitochondrial, chloroplast, or archaeal, 123,828 bacterial ASVs were successfully identified from 8,261,802 non-chimeric reads, encompassing 50 bacterial phyla, 131 classes, and 2645 genus. α-diversity analyses were conducted at the ASV level based on Chao 1, Shannon, Simpson, and ACE diversity indices. In addition, beta diversity analyses were performed at the ASV level using Bray-Curtis dissimilarity. The asymptotic behavior of the sample dilution curves confirmed the representativeness of the sampling depth. All diversity analyses were completed on the BMKCloud microbiome diversity analysis platform (BMKCloud, 2022).
Species taxonomic composition and distribution of the seed microbiome
Taxonomic composition was assessed at the class and genus levels across six distinct microhabitats: rhizosphere, root, stem, leaf, seed, and seed-P. The top 20 most abundant classes and genera were also identified. The filtering criteria included a minimum count threshold of 1, presence in >10% of the samples, and consolidation of taxa with counts <10. This taxonomic analysis was performed on 50 bacterial phyla and 2645 genera. At the phylum level, taxa with ASV counts below 10 were categorized into the “other”. At the genus level, the first 50 dominant bacteria were identified, with the remainder classified as “other”, whereas only the top 20 positions were explicitly labeled. The microbiome distributions were visually represented using stacked histograms created using the ggplot2 tool in R (Wickham, 2016). To investigate the vertical distribution of microorganisms during the seed-to-seedling transition, a comparative analysis was conducted between microbiomes in different tissues of germinating seeds and seedlings. These microbiomes were subsequently compared to the bacterial and fungal populations present in the seeds before sowing. This study postulated that if there is a transfer of ASVs from seeds to various plant tissues (including roots, stems, leaves, seeds and seeds-P) of both seedlings and mature host plants, there should be a similarity in the microbiome composition across these tissues. Community dissimilarity was quantified using the vegan package in R programming language. The Jaccard dissimilarity index was used to measure the dissimilarity between the two communities by calculating the ratio of shared ASVs to the total number of ASVs, regardless of their individual abundances. To investigate the phenomenon of microbiome migration from seeds to host plants during the planting process, we initially identified ASVs present in both seeds and various cotton plant tissues (leaves, stems, roots, and progeny seeds). These specific ASVs were denoted as “seed source ASVs”. The relative proportion of seed-source ASVs was estimated by determining the ratio of the number of seed-source ASVs to the total number of ASVs in each tissue. Bar graphs were used to visually depict the composition of seed-derived ASVs, presenting their cumulative relative abundance.
Selection of core endophytes in cotton seeds
Common ASVs across all plant microhabitats were visualized using Venn diagrams with the InteractiVenn1 (Heberle et al., 2015). The term “common ASVs” refers to ASVs that are identified in all samples within each group chamber. The core ASVs present in cotton were determined using frequency analysis with a 70% threshold. The SourceTracker package (Knights et al., 2011) (version 1.0.1) in the R software (version 3.1.1) was used to analyze and identify the microbiome sources of cotton progeny seeds. This analysis employed a Bayesian algorithm that explores microbiome origins (source) in the target sample (sink), using the seed sample as the sink and the seed-P, along with each microhabitat samples, as the sources. Predictions were based on the distribution of community structures in the source and sink samples, enabling the estimation of the proportion of progeny seed composition originating from distinct cotton tissue parts and progeny seeds.
Network construction
We conducted a comprehensive analysis of the bacterial community networks associated with cotton within six distinct microhabitats. Networks were constructed at the ASV level using the SparCC algorithm, and significant interactions between ASVs (p < 0.05) were identified in the network from ASVs with relative abundance greater than 0.01%. These significant interactions indicate relationships of interdependence or influence observed within the microbial community. Revealing potential functional interactions within the microbial community is of great importance for understanding the structure and function of these communities. ASVs present in at least 25% of the data were evaluated using the construction criteria determined through random matrix theory (RMT). Network analysis, including correlation analysis and modularity identification, was performed using the Molecular Ecological Network Analysis (MENA) platform (Deng et al., 2012; Feng et al., 2022), accessible at http://ieg2.ou.edu/MENA/, as described by Deng and Feng. Network visualization was executed utilizing Gephi 0.9.2 software (Bastian et al., 2009), while network modularity assessment relied on an optimized greedy method. Furthermore, intramodular connectivity (Zi) and intermodular connectivity (Pi) were calculated for each network node to identify the potential keystone species. Nodes were classified into four categories based on their network structure: module centers, network centers, links, and peripheral nodes. Graphics in this study were generated using ggplot2 (version 3.3.2) (Wickham, 2016).
Functional gene prediction
To comprehensively elucidate the predominant microbiome biological processes in various cotton tissues, including the roots, stems, and leaves, we employed two distinct functional prediction methods. The first method involves a manual correlation of prokaryotic taxa, such as genera, with metabolic or ecological functions, such as nitrification, denitrification, or fermentation. This correlation relies on information derived from literature concerning cultured representative microorganisms. To facilitate this process, we utilized the FAPROTA database (version 1.2.6) as outlined by Louca et al. (2016). A Python script provided by FAPROTAX was used to convert the ASV table into a putative functional table. This conversion was based on the taxa detected in the samples and their corresponding functional annotations in the FAPROTAX database.
To achieve a more comprehensive prediction of microbiome phenotypes, we used BugBase (Ward et al., 2017). This tool can analyze 16S rRNA gene sequence data to predict a range of phenotypic characteristics, including redox conditions (aerobic or anaerobic), presence or absence of mobile elements, biofilm formation, Gram stain response (positive or negative), potential pathogenicity, and stress resistance.
Functional and phenotypic predictions obtained from the two techniques were further evaluated by comparison using the Wilcoxon rank-sum test within the Statistical Analysis of Metagenomic Profiles (STAMP v2.1.3) (STAMP, 2015) software. The main objective of this integrated approach was to provide a comprehensive perspective on the functional diversity and potential ecological roles of microbiomes within various cotton tissue types.
Statistical analysis and visualization
All statistical analyses were performed using R software (R Core Team, 2013) (version 2023.06.1 + 524) as the default platform, unless otherwise specified, with a predetermined statistical significance threshold of alpha = 0.05. Adjustments for multiple hypothesis testing were conducted using the false positive detection rate (FDR) when applicable. Sparsity corrections were conducted for the α-diversity analysis. The microbiome software package for R (Lahti and Shetty, 2019) (version 1.9.13) was used to obtain Shannon and Simpson diversity values using its alpha function. For microorganisms with relative abundances exceeding 0.5%, their distributions were visualized utilizing the ggplot2 package for R (Wickham, 2016) (version 3.2.1). Assessment of beta diversity among samples employed the weighted UniFrac algorithm, which was based on ASV abundance and composition. Finally, the Principal Coordinate Analysis (PCoA) tool in R (Gower, 1966) was employed for graphical analysis of sample similarities.
Results
Quality metrics of sequencing analysis
To methodically investigate cotton-associated microbiomes, we performed DNA amplification and sequencing targeting the V3–V4 region of the bacterial 16S rRNA gene. From a comprehensive dataset comprising 72 samples, a total of 12,605,376 raw reads were generated. These samples included 12 rhizosphere soil samples, 12 root samples, 12 stem samples, 12 leaf samples, 12 seed samples, and 12 seed-P samples. The average number of non-chimeric reads per sample was 114,747 (Supplementary Table 1). After rigorous quality control to eliminate chimeric and plant-derived reads, 8,261,802 non-chimeric reads were retained. These reads were subsequently classified into 123,828 bacterial ASVs and are detailed in Supplementary Table 2. The analysis of 16S rRNA ASVs included their distribution within a diverse spectrum of bacterial phyla and genera, with 50 phyla and 2,645 genera identified (Supplementary Table 2). To assess the comprehensiveness of the sequencing process, a dilution study was conducted. The results of this analysis confirmed the adequacy of our sequencing methodology for capturing the true microbiome diversity present in the samples (Supplementary Figure 1A).
Diversity and drivers of the cotton microbiome
We conducted a comprehensive investigation of cotton-associated microbiome diversity using α-diversity indices (abundance and richness, using the Chao1 and ACE indices). The diversity and evenness were assessed using the Shannon and Simpson indices for each sample. Microbiomes from different cotton genotypes were analyzed separately in the rhizospheres, roots, stems, leaves, seeds and seeds-P. The results indicated no significant differences in the alpha diversity of endophytic bacteria within seeds-P across different genotypes (Figures 1B–E and Supplementary Table 3). Additionally, the α-diversity indices of endophytic microbiomes in the rhizospheres, roots, stems, leaves, seeds and seeds-P of various cotton genotypes did not exhibit significant differences overall (Figures 1F–I and Supplementary Figure 2). This suggests that microbiome diversity in different ecological niches of cotton was minimally affected by the genotype. Conversely, microbiome α-diversity within each cotton site (rhizosphere, root, stem, leaf, seed and seed-P) exhibited significant variation across microhabitats (Figures 1G–M and Supplementary Table 3). These results suggest that microhabitats, rather than cotton genotypes, primarily dictate the microbiome diversity. To further confirm this observation, we performed a principal coordinate analysis (PCoA) to elucidate the underlying relationships among the factors (Figure 2). The PCoA results strongly indicated substantial differences in the distribution of the cotton microbiome across ecological niches (e.g., rhizospheres, roots, stems, leaves, seeds and seeds-P) (Figure 2A), with these distinctions not being attributed to cotton genotypes (Figure 2B). Microbiomes from various ecological niches exhibited convergence in rhizospheres, roots, stems, leaves, seeds and seeds-P, although significant overlap was observed in the internal microbiome diversity of roots, stems, and leaves (Figures 2C–H). We evaluated six physicochemical characteristics of rhizosphere soil samples from three cotton genotypes and observed no significant differences among these parameters, indicating that the soil fertility is comparable among the three genotypes (Supplementary Figure 1B).
FIGURE 2
Microhabitats have emerged as the predominant determinants of diversity within the cotton microbiome, revealing a distinct spatial distribution pattern in microbiome diversity associated with cotton.
Taxonomic composition and dynamics of the cotton microbiome
To comprehensively investigate the taxonomic characteristics and dynamics of the cotton microbiome under varying conditions, we conducted a systematic analysis across multiple microhabitats, including rhizospheres, roots, stems, leaves, seeds and seeds-P. Figure 3, Supplementary Figure 3, and Supplementary Table 4 illustrate different degrees of enrichment or depletion within the cotton microbiome at both the phylum and genus taxonomic levels. The analysis revealed Gammaproteobacteria as the dominant taxon across most microhabitats, constituting over 13% of the ASVs. Nevertheless, the abundance and distribution of Gammaproteobacteria exhibited significant variation among microhabitats (Supplementary Figure 3 and Supplementary Table 4). For example, Gammaproteobacteria accounted for 22.9% and 18.4% of ASVs in seeds-P and roots, respectively, whereas Alphaproteobacteria were the most abundant in rhizospheres, reaching 29.6% of ASVs. Certain bacterial groups (class), such as Clostridia, Bacilli, and Bacteroidia, showed a gradual increase from seeds-P to roots, stems, leaves, and seeds. Their abundance was relatively low in rhizosphere, but progressively increased within plant tissues, particularly in stems, leaves, and progeny seeds (Figures 3A–C and Supplementary Table 4). Particularly noteworthy is the change in microbiome composition from seed-P to seeds, with a marked decrease in the relative abundance of Gammaproteobacteria, and an increase in the relative abundance of Clostridia, Bacilli, and Bacteroidida.
FIGURE 3
In addition, Bacilli were identified as the predominant bacterial group during cotton development from root to seed, accompanied by a gradual decline in unclassified_Sphingomonadaceae and certain unclassified_Bacteria genera within the Alphaproteobacteria. This pattern indicated a significant influence of seed-derived microorganisms on the aboveground portion, whereas rhizosphere microbiomes were predominantly affected by the field environment.
Key taxa and co-occurrence networks of the bacteria
To comprehensively investigate the bacterial interactions within cotton microbiomes across distinct microhabitats, including rhizosphere, roots, stems, leaves, seeds, and seeds-P, we constructed six bacterial co-occurrence networks (Figure 4). Overall, these networks predominantly exhibited positive correlations among the bacterial taxa. A detailed analysis of network topological parameters revealed that microbiomes in stem (modularity index 0.438) and leaf (modularity index 0.448) microhabitats demonstrated more pronounced modular structures than those in roots (modularity index 0.427) and rhizosphere (modularity index 0.421) microhabitats (Figure 4B). Notably, the seeds-P exhibited the highest values for network modularity, average clustering coefficient, and average path length, indicative of structured, strongly correlated, and intricate microbiome networks (Figure 4). Conversely, daughter seeds displayed relatively lower values for these metrics, suggesting more fragmented and less-correlated microbiome networks. Among all assessed microhabitats, the ratio of positive to negative connections (i.e., edges) was lowest in leaf microhabitats (1.47) and highest in rhizosphere microhabitats (5.15), highlighting the varying degrees of ecologically significant effects on bacterial interactions across different microhabitats (Figure 4B).
FIGURE 4
Microbiomes often involve microorganisms that serve as central taxa, commonly referred to as central nodes, within their networks. These central nodes include modular centers and links, and their identification involves assessing the intra-module degree (Zi) and inter-module connectivity (Pi) of each ASV within the network. Across all networks, 279 connections, 35 module centers, and two network hubs were observed (Supplementary Figure 4 and Supplementary Table 5). The stem and leaf networks indicated a comparatively higher number of central nodes, with 138 and 110 nodes, respectively. In contrast, offspring seed and seed-P networks exhibited a notable difference in the central node count, featuring only 1 and 30 nodes, respectively (Figure 4B and Supplementary Table 5). Rhizosphere, stem, and leaf networks predominantly featured Bacilli as the major taxa. Furthermore, the stem and leaf networks identified crucial taxa (Supplementary Table 6), demonstrating the diverse ecological significance of different bacterial class within various cotton regions. Despite the limited number of central nodes in the seed network, these nodes exhibited the highest mesocosm centralities. This observation highlights their pivotal role in maintaining network structure stability, as indicated in Supplementary Table 5.
Endophytic microbiomes of cotton seed cores
Cotton microbiomes exhibit distinct ecological niche-dependent dynamics. However, a consistent set of core microbiomes within cotton seeds remains prevalent, demonstrating their relative stability and resilience to external environmental influences. To comprehensively explore the composition and characteristics of this endophytic core microbiome, we synthesized and analyzed extensive data from diverse cotton genotypes and microhabitats.
The study results revealed the presence of 23,546, 32,336, 25,432, 15,518, and 38,407 ASVs with frequencies exceeding 70% in samples obtained from the roots, stems, leaves, seeds, and seeds-P, respectively (Figure 5A). Particularly noteworthy was the identification of 1,839 shared ASVs in the Venn diagram among the root, stem, and leaf tissues, constituting 7.81%, 5.69%, and 7.23% of the ASVs within each tissue. In addition, 130 ASVs overlapped with the progeny seed roots, accounting for 0.84% of the total progeny seed ASVs. Meanwhile, 687 ASVs overlapped with the progeny seeds, representing 1.79% of the total progeny seed ASVs (Figure 5A and Supplementary Table 7). SourceTracker software and Bayesian algorithm-based analysis further elucidated the potential sources of microbiomes in the progeny seeds. Based on the microbiome structure of the source and sink samples, predictive analysis indicated that the largest proportion (11.34%) of ASVs in progeny seeds originated from seeds-P, whereas roots, stems, and leaves contributed 4.83%, 3.18%, and 8.32%, respectively. ASVs from the rhizosphere environment accounted for only 0.24% of the total ASVs, while the source of the remaining ASVs was unknown, comprising 72.09% (Figure 5B). Collectively, 334 core ASVs were successfully identified in five major cotton endophytic tissues (roots, stems, leaves, seeds and seeds-P), constituenting the core cotton endophytic microbiome. These core ASVs represented 0.87% of the relative abundance of progeny seeds. Further traceability analysis identified 29 core ASVs unique to progeny seeds distributed across seven microbiome classes, including Acidimicrobiia, Alphaproteobacteria, Bacilli, Bacteroidia, Clostridia, Gammaproteobacteria, and unclassified_Proteobacteria (Figures 5A, B and Supplementary Table 7). Among the core endophytic microbiome, Clostridia and Gammaproteobacteria were the most enriched taxa, each featuring 10 associated sequence variants (ASVs, Figure 5B). Remarkably, these common ASVs constituted 8.68% of the core cotton-endophytic microbiome. In particular, ASV_1070 of the Acidimicrobiia class exhibited substantial absolute abundance and a high occurrence frequency in the root microhabitat, surpassing other microhabitats (Figure 5C). Furthermore, three genera, Bacilli, Clostridia, and Gammaproteobacteria, exhibited significantly higher absolute abundance and frequency of occurrence in seeds-P than in other microhabitats, followed closely by seeds, while they were relatively scarce in rhizosphere microhabitats. These observations reinforced the hypothesis that the 29 ASVs, spanning seven different genera of microorganisms, were potentially transmitted through vertical pathways. This finding indicates selective transmission of only a few specific microorganisms during vertical transmission, suggesting a selective element in the transmission mechanism (Cordovez et al., 2019).
FIGURE 5
Predictive functional characterization of cotton bacterial communities
Comprehensive analyses using functional prediction tools such as FAPROTAX and Bugbase, which reference the SILVA rRNA and Greengenes databases, were employed to preliminarily explore the metabolic potential and phenotypic characteristics of microbiomes within various microhabitats of cotton (Gossypium spp.). These tools, based on 16S rRNA gene sequence analysis, aim to provide an initial understanding of the overall functional trends of microbiomes. Our findings indicate that the microbiome associated with roots exhibits heightened nitrate reduction capacity (Figure 6A), a crucial metabolic function for enhancing plant growth (Guilherme et al., 2019). Similarly, microbiomes in seeds-P demonstrated increased activities related to sulfur oxidation reduction and iron reduction (Figure 6C), potentially influencing their growth and survival in soil (Nakajima et al., 2019; Malviya et al., 2022). Moreover, microbiomes in microenvironments other than the rhizosphere showed higher activities in fermentation functions and anaerobicity (Figure 6B). A notable change was observed in the relative abundance of Gram-positive bacteria between seeds-P and seeds (Figure 6D), potentially reflecting an adaptation process that confers a survival advantage to Gram-positive microorganisms in seeds. Significant variations were also noted in various microbiome activities between seeds-P and seeds, including sulfate respiration, sulfide respiration, xylan degradation, and photoheterotrophy, potentially signifying shifts in intergenerational adaptation, impacting the viability and functional behavior of microorganisms with diverse roles in seeds (Moreira and Filho, 2016; Couturier et al., 2018; Graham et al., 2018).
FIGURE 6
The integrated analysis of these functional predictions provides preliminary insights into the functional mechanisms and ecological characteristics of microbiomes in various ecological niches of cotton plants, but should be viewed with caution once it is based only on predictions from short region of 16S rRNA gene sequences and may vary greatly among taxonomic groups and intra-specifically. These findings reveal the potential contributions of microbiomes to plant growth and health, while highlighting significant biological alterations that may occur across generations. Such changes could affect the overall stability and functionality of microbiomes, thereby impacting the ecological role of cotton in its environment. Further research is warranted to explore the accuracy of these functional predictions and to gain a deeper understanding of the potential long-term ecological implications for both plants and ecosystems.
Discussion
This study aimed to comprehensively analyze the diversity, structure, and functional properties of microbiomes inhabiting different microhabitats within cotton (Gossypium spp.). These microorganisms engage in diverse beneficial interactions with host plants, including mutualistic symbiosis, co-evolution, and benign parasitism (Hardoim et al., 2008; Kaga et al., 2009). In addition, host plants can transfer beneficial endophytic microorganisms vertically to their progenies (Mastretta et al., 2009; Compant et al., 2010). To identify crucial vertically transferred taxa within the cotton microbiome, we analyzed microorganism populations originating from five distinct microhabitats across the three cotton genotypes. High-throughput 16S rRNA gene sequencing and diversity assessments were employed to elucidate the collective influence of different microhabitats (e.g., rhizosphere, root, stem, leaf, seed and seed-P) and cotton genotypes on microbiome diversity and structure. Our findings revealed that microhabitats exerted a more substantial influence on microbiome diversity than genotypes. This is consistent with prior research, such as the study of rice microbiome, which identified microhabitats as the primary determinant rather than environmental factors or host genotypes (Zhang et al., 2022). Furthermore, this study demonstrated distinct patterns of microorganism population dynamics in different cotton tissues. Analysis of microbiome diversity and network interactions revealed a significant decrease in both microbiome diversity and network complexity from roots, stems, and leaves to seeds. This trend may suggest a selective effect by the plant host, preferring to vertically transmit beneficial microorganisms to its progeny (Jeong et al., 2021; Zhang et al., 2022). Utilizing a variety of bioinformatics approaches, we conducted a detailed analysis of the microbiomes associated with cotton plants, aiming to define their core microbiome and to thoroughly investigate the vertical transmission of endophytic bacteria in seeds. In this study, we particularly noted the significant impact of environmental changes on the microbiome structure between seeds-P and seeds. By comparing the microbiome composition of seed-P and seeds, we revealed the dynamics of vertical transmission of microbiomes within the cotton lifecycle. This finding suggests that the adaptability and evolution of microbiomes are not only related to the genetic traits of the plant itself but may also be significantly influenced by external environmental factors. Subsequently, we generated predictions regarding the functional characteristics of microbiomes in each microhabitat, thus contributing to our understanding of the potential interactions between plants and microorganism populations. The results of this study offer a comprehensive and validated body of knowledge, enhancing our understanding of the vertical transfer of microbiomes within seeds and shedding light on the co-evolutionary dynamics between hosts and microorganisms.
Alphaproteobacteria and Bacilli indicated noteworthy dynamics across diverse microenvironments within cotton. This study revealed that Alphaproteobacteria exhibited a widespread distribution in the rhizosphere environment with a high abundance, indicating their specific adaptation to external conditions. These microorganisms are involved in a wide array of biological activities such as nitrogen fixation, release of biotrophic factors, pathogen antagonism, and mineral bioactivation (Dedysh et al., 2004; Taulé et al., 2012; Castanheira et al., 2014; Wicaksono et al., 2023). Bacterial communities were notably enriched within seeds across various cotton microhabitats, which is consistent with previous studies of cotton seed microbiomes (Adams and Kloepper, 2002). The ability of Bacillus microorganisms to produce stress-tolerant spores adapted to extreme environments may elucidate their successful colonization within seeds, potentially contributing to seed resilience to adverse conditions and pathogenic threats. Additionally, this trait may induce vertical transmission between plants and microorganisms (Gagné-Bourque et al., 2015), suggesting that Bacillus could employ seeds as a medium for persistent ecological impacts throughout the plant lifecycle. At the class level, Clostridia, Bacteroidia, and Bacilli were identified as the abundant resident microorganisms in cotton. Clostridia and Bacteroidia are more crucial for cotton ecosystems than Bacilli. Clostridia are pivotal in anaerobic decomposition and maintenance of soil structure, potentially promoting the health and stability of plant root systems (Mowlick et al., 2013; Huang et al., 2019). Bacteroidia primarily contribute to nutrient uptake and plant health and may play a role in plant disease resistance and nitrogen fixation processes (Pan et al., 2023). These microbiome taxa potentially exert significant functional impacts on the ecological stability, growth, and adaptation to adversity in cotton seeds. Therefore, a comprehensive exploration of the distribution and functional attributes of these microorganisms across different growth stages and cotton tissues is essential for gaining deeper insights into their integrated roles within the cotton ecosystem. Within the Clostridia class, the relative abundance of Lachnospiraceae increased significantly from the rhizosphere to the seed microenvironment. Notably, Bacillus spp. accounted for 10.36% of the seed microbiome. Previous studies tracking bacterial titers in seedlings and employing scanning electron microscopy have demonstrated the presence of Bacillus at high titers in stems, leaves, and pods of oilseed rape plants after 1 month of initial treatment (Reva et al., 2004). The ability of Bacillus to propagate vertically between plants through seeds was further confirmed by the isolation of several strains originating from seedlings. Overall, our results substantiate the persistence of cotton seed endophytes within seeds and validate the theory of their uniform distribution throughout plant growth, including intra-tissue migration.
The core endophytic microbiome of seeds can engage in co-evolutionary processes with their host plants over the course of the plant life cycle and develop efficient dispersal mechanisms across multiple generations (Truyens et al., 2014; Frank et al., 2017; Shade et al., 2017; Berg and Raaijmakers, 2018; Nelson, 2018). Therefore, identifying core microbiome groups is of critical importance for comprehending stable and consistent constituents within complex microbiomes (Shade and Handelsman, 2012). Core taxa, discerned through network analysis, such as Bacilli and Cytophagales, can be pivotal ‘bridges’ connecting microbiomes with host plants. The presence and activity of these core taxa can significantly affect the overall stability and functionality of the microbiome networks (Berg, 2009; Schmidt et al., 2019; Ikeda et al., 2023). In this study, through in-depth analysis of microbiomes, we identified seven core bacterial classes residing within cotton seeds: Acidimicrobiia, Alphaproteobacteria, Bacilli, Bacteroidia, Clostridia, Gammaproteobacteria, and unclassified_Proteobacteria. The core bacterial classes identified in this study may possess critical biological functions related to soil nutrient cycling, plant growth promotion, and disease resistance. For instance, within the Alphaproteobacteria class, four genera: Bradyrhizobium, Azospirillales, Brevundimonas, and Allorhizobium_Neorhizobium_Pararhizobium_Rhizobium have been noted. The genus Bradyrhizobium, known for its significant role in nitrogen fixation, might play a vital part in enhancing soil nitrogen content surrounding cotton plants (Pitumpe Arachchige et al., 2020; Favero et al., 2022). The order Azospirillales is recognized for its plant growth-promoting capabilities and may play a crucial role in the development of cotton seeds (Zhao et al., 2022). Genera Brevundimonas and Allorhizobium_Neorhizobium_Pararhizobium_Rhizobium have demonstrated adaptability and ecological functions in extreme environments (Ding et al., 2021). In the Bacteroidia class, the genus Bacteroides is believed to play a significant role in the decomposition of organic matter, potentially facilitating the nutrient cycling in the soil surrounding cotton seeds, which is essential for cotton growth (Wang et al., 2016). Certain strains within the Bacilli class also exhibit significant nitrogen-fixing and phosphate-solubilizing activities (Gelfand and Ravcheev, 2016). These microbiome taxa not only contribute to the diversity of cotton seed ecosystems, but also possess the potential to promote cotton growth and enhance agricultural yields. Gammaproteobacteria, as identified in several studies, has been identified as a class of disease-resistant facultative bacteria with versatile functionalities, including microbiome growth inhibition, nitrogen fixation, phosphorus solubilization, iron carrier production, and phytohormone synthesis (Walterson and Stavrinides, 2015). Moreover, the genus Bacillus has been associated with the growth of various plant seeds, such as maize (Johnston-Monje and Raizada, 2011; Liu et al., 2012; Tkalec et al., 2022), rice (Mano et al., 2006; Zhang et al., 2019), and radish (Rezki et al., 2018), where it is the dominant genus in the endophytic microbiomes. Our observations revealed substantial enrichment of the genus Bacillus, which dominated the endophytic microbiomes of seeds during the transition from seeds-P to seeds. Comparable findings have been substantiated in other plant seeds, including Arabidopsis thaliana (Truyens et al., 2013), rice (Mano et al., 2006), and kidney beans (López-López et al., 2010). Cotton yellow wilt, often referred to as “cotton blight,” is a detrimental disease caused by the soil-borne pathogen Verticillium dahliae and has substantial adverse effects on cotton yield and quality (Wang et al., 2021). Notably, Bacillus strains belonging to the Bacillaceae family have demonstrated effectiveness in inhibiting the growth and propagation of Ralstonia solanacearum, providing a potential avenue for biocontrol to manage this challenging disease (Sussanti et al., 2023). Through functional prediction analysis, we discovered that bacterial communities within different microhabitats of cotton exhibit a variety of functional traits. For instance, bacterial communities in the rhizosphere microenvironment demonstrated significant nitrogen-fixing capabilities, whereas those closely associated with roots showed enhanced nitrate reduction capacity. Additionally, we observed noticeable differences in microbiome functions between seeds and their progenies, which may reflect the mutual selection between host and microorganisms, aiding the plant’s adaptation to diverse environmental conditions (Ferreira et al., 2008; Rudgers et al., 2009). These preliminary findings suggest a potential key role of the seed core microbiomes in establishing beneficial microbiome symbioses essential for the health of plant offspring (Truyens et al., 2014). To deepen our understanding of the interactions and co-evolutionary mechanisms between these endophytes and their host plants, future research should concentrate on the isolation and purification of these strains, coupled with comprehensive molecular-level investigations. This approach provides valuable insights into how these microorganisms establish and sustain beneficial symbiotic relationships within cotton seeds.
Conclusion
This study provides an extensive perspective on microbiome dynamics within cotton seeds throughout their growth and developmental cycles by employing amplicon-based microbiome analysis. This revealed the efficient migration of seed-derived bacteria into both aboveground (e.g., stems and leaves) and belowground (e.g., roots) plant structures. The identification of core microorganisms intimately associated with the interior of cotton seeds, pivotal network nodes, and core microbiome that underwent vertical transmission provides a theoretical foundation for microbiome regulation aimed at enhancing plant adaptation.
Future research should involve the isolation and functional characterization of these core microorganisms, along with an exploration of host-related factors, including gene expression, phytohormones, and metabolites that could affect their distribution during seed germination, seedling growth, and various stages of plant maturity. These investigations can enhance our understanding of how host factors regulate the temporal dynamics of microbiomes and how they influence plant physiological processes. In summary, this study establishes a robust theoretical and practical basis for customized microbiome applications in agriculture, industry, and environmental sciences.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA1028990.
Author contributions
CW: Data curation, Formal Analysis, Methodology, Writing – original draft, Software, Visualization. XZ: Data curation, Formal Analysis, Methodology, Writing – original draft. LD: Investigation, Methodology, Writing – original draft. YW: Data curation, Investigation, Writing – original draft. YR: Investigation, Writing – original draft. HY: Methodology, Writing – original draft. AW: Funding acquisition, Project administration, Writing – review and editing. YF: Investigation, Methodology, Writing – original draft. JY: Data curation, Investigation, Writing – original draft. RL: Project administration, Writing – review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the following grants: High-level Talent Research Startup Project of Shihezi University (No. RCZK202366), and the National College Student Innovation and Entrepreneurship Training Program (No. 202210759071).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1323342/full#supplementary-material
References
1
AdamsP. D.KloepperJ. W. (2002). Effect of host genotype on indigenous bacterial endophytes of cotton (Gossypium hirsutum L.).Plant Soil240181–189. 10.1023/A:1015840224564
2
AdhamB. M.ShavkatC. B.DilbarA. A.MakhbubaS. K.ShakhbozjonH. M. (2020). Cultivation of Cotton in Uzbekistan and Xinjiang (NW China) Using Drip Irrigation.IJITEE93547–3549. 10.35940/ijitee.B7235.019320
3
BaiJ. (2022). High-yield cultivation and pest control techniques of cotton in Xinjiang.Seed Sci. Technol.4022–24.
4
BastianM.HeymannS.JacomyM. (2009). Gephi: An open source software for exploring and manipulating networks.ICWSM3361–362. 10.1609/icwsm.v3i1.13937
5
BergG. (2009). Plant–microbe interactions promoting plant growth and health: perspectives for controlled use of microorganisms in agriculture.Appl. Microbiol. Biotechnol.8411–18. 10.1007/s00253-009-2092-7
6
BergG.RaaijmakersJ. M. (2018). Saving seed microbiomes.ISME J.121167–1170. 10.1038/s41396-017-0028-2
7
BMKCloud (2022). BMKCloud. Available online at: https://international.biocloud.net/zh/user/login?redirect=%2Fzh%2Fdashboard(accessed May 10, 2023).
8
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics30, 2114–2120. 10.1093/bioinformatics/btu170
9
BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC. C.Al-GhalithG. A.et al (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2.Nat. Biotechnol.37852–857. 10.1038/s41587-019-0209-9
10
CallahanB. J.McMurdieP. J.RosenM. J.HanA. W.JohnsonA. J. A.HolmesS. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data.Nat. Methods13581–583. 10.1038/nmeth.3869
11
CastanheiraN.DouradoA. C.AlvesP. I.Cortés-PalleroA. M.Delgado-RodríguezA. I.PrazeresÂet al (2014). Annual ryegrass-associated bacteria with potential for plant growth promotion.Microbiol. Res.169768–779. 10.1016/j.micres.2013.12.010
12
ChesneauG.Torres-CortesG.BriandM.DarrasseA.PreveauxA.MaraisC.et al (2020). Temporal dynamics of bacterial communities during seed development and maturation.FEMS Microbiol. Ecol.96fiaa190. 10.1093/femsec/fiaa190
13
CompantS.ClémentC.SessitschA. (2010). Plant growth-promoting bacteria in the rhizo- and endosphere of plants: Their role, colonization, mechanisms involved and prospects for utilization.Soil Biol. Biochem.42669–678. 10.1016/j.soilbio.2009.11.024
14
CordovezV.Dini-AndreoteF.CarriónV. J.RaaijmakersJ. M. (2019). Ecology and evolution of plant microbiomes.Annu. Rev. Microbiol.7369–88. 10.1146/annurev-micro-090817-062524
15
CouturierM.LadevèzeS.SulzenbacherG.CianoL.FanuelM.MoreauC.et al (2018). Lytic xylan oxidases from wood-decay fungi unlock biomass degradation.Nat. Chem. Biol.14306–310. 10.1038/nchembio.2558
16
Da SilveiraA. P. D.IórioR. D. P. F.MarcosF. C. C.FernandesA. O.De SouzaS. A. C. D.KuramaeE. E.et al (2019). Exploitation of new endophytic bacteria and their ability to promote sugarcane growth and nitrogen nutrition.Anton. Leeuw.112283–295. 10.1007/s10482-018-1157-y
17
DedyshS. N.RickeP.LiesackW. (2004). NifH and NifD phylogenies: an evolutionary basis for understanding nitrogen fixation capabilities of methanotrophic bacteria.Microbiology1501301–1313. 10.1099/mic.0.26585-0
18
DengY.JiangY.-H.YangY.HeZ.LuoF.ZhouJ. (2012). Molecular ecological network analyses.BMC Bioinformatics13:113. 10.1186/1471-2105-13-113
19
DingY.WeiR.WangL.YangC.WangH.LiH. (2021). Diversity and dynamics of microbial ecosystem on berry surface during the ripening of ecolly grape in Wuhai, China.Res. Square [Preprint]. 10.21203/rs.3.rs-656990/v1
20
Dorn-InS.BassittaR.SchwaigerK.BauerJ.HölzelC. S. (2015). Specific amplification of bacterial DNA by optimized so-called universal bacterial primers in samples rich of plant DNA.J. Microbiol. Methods11350–56. 10.1016/j.mimet.2015.04.001
21
EdgarR. C. (2016). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing.Bioinformatics10.1101/081257[Epub ahead of print].
22
EdgarR. C.HaasB. J.ClementeJ. C.QuinceC.KnightR. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics27, 2194–2200. 10.1093/bioinformatics/btr381
23
FaveroV. O.De CarvalhoR. H.LeiteA. B. C.Dos SantosD. M. T.De FreitasK. M.BoddeyR. M.et al (2022). Bradyrhizobium strains from Brazilian tropical soils promote increases in nodulation, growth and nitrogen fixation in mung bean.Appl. Soil Ecol.175:104461. 10.1016/j.apsoil.2022.104461
24
FengK.PengX.ZhangZ.GuS.HeQ.ShenW.et al (2022). iNAP: An integrated network analysis pipeline for microbiome studies.iMeta1:e13. 10.1002/imt2.13
25
FerreiraA.QuecineM. C.LacavaP. T.OdaS.AzevedoJ. L.AraðjoW. L. (2008). Diversity of endophytic bacteria from Eucalyptus species seeds and colonization of seedlings by Pantoea agglomerans.FEMS Microbiol. Lett.2878–14. 10.1111/j.1574-6968.2008.01258.x
26
FrankA.Saldierna GuzmánJ.ShayJ. (2017). Transmission of bacterial endophytes.Microorganisms5:70. 10.3390/microorganisms5040070
27
Gagné-BourqueF.MayerB. F.CharronJ.-B.ValiH.BertrandA.JabajiS. (2015). Accelerated growth rate and increased drought stress resilience of the model grass brachypodium distachyon colonized by Bacillus subtilis B26.PLoS One10:e0130456. 10.1371/journal.pone.0130456
28
GelfandM. S.RavcheevD. A. (2016). Nitrogen fixation and molecular oxygen: Comparative genomic reconstruction of transcription regulation in alphaproteobacteria.Front. Microbiol.7:1343. 10.3389/fmicb.2016.01343
29
GopalM.GuptaA.ThomasG. V. (2013). Bespoke microbiome therapy to manage plant diseases.Front. Microbiol.4:355. 10.3389/fmicb.2013.00355
30
GowerJ. C. (1966). Some distance properties of latent root and vector methods used in multivariate analysis.Biometrika53325–338. 10.1093/biomet/53.3-4.325
31
GrahamE. D.HeidelbergJ. F.TullyB. J. (2018). Potential for primary productivity in a globally-distributed bacterial phototroph.ISME J.121861–1866. 10.1038/s41396-018-0091-3
32
GuilhermeE. A.NascimentoC. S.LoboA. K. M.CarvalhoF. E. L.SilveiraJ. A. G. (2019). Nitrogen-utilization efficiency during early deficiency after a luxury consumption is improved by sustaining nitrate reductase activity and photosynthesis in cotton plants.Plant Soil443185–198. 10.1007/s11104-019-04214-7
33
HardoimP. R.AndreoteF. D.Reinhold-HurekB.SessitschA.Van OverbeekL. S.Van ElsasJ. D. (2011). Rice root-associated bacteria: insights into community structures across 10 cultivars: Insights into the bacterial community structure of rice cultivars.FEMS Microbiol. Ecol.77154–164. 10.1111/j.1574-6941.2011.01092.x
34
HardoimP. R.Van OverbeekL. S.ElsasJ. D. V. (2008). Properties of bacterial endophytes and their proposed role in plant growth.Trends Microbiol.16463–471. 10.1016/j.tim.2008.07.008
35
HeberleH.MeirellesG. V.Da SilvaF. R.TellesG. P.MinghimR. (2015). InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams.BMC Bioinformatics16:169. 10.1186/s12859-015-0611-3
36
HillR.SparlingG. (2009). “Soil quality monitoring,” in Land and soil monitoring: A guide for SoE and regional council reporting, Hamilton: Land Monitoring Forum, 27–88.
37
HoneH. (2021). Profiling, isolation and characterisation of beneficial microbes from the seed microbiomes of drought tolerant wheat.Sci. Rep.11:11916.
38
HouH. (2023). Study on the impact of meteorological disasters during cotton growth period in Xinjiang on cotton growth and countermeasures.Res. Agric. Disasters1322–24.
39
HuangX.LiuL.ZhaoJ.ZhangJ.CaiZ. (2019). The families Ruminococcaceae, Lachnospiraceae, and Clostridiaceae are the dominant bacterial groups during reductive soil disinfestation with incorporated plant residues.Appl. Soil Ecol.13565–72. 10.1016/j.apsoil.2018.11.011
40
IkedaS.OkazakiK.TakahashiH.TsurumaruH.MinamisawaK. (2023). Seasonal shifts in bacterial community structures in the lateral root of sugar beet grown in an andosol field in Japan.Microb. Environ.38:22071. 10.1264/jsme2.ME22071
41
JeongS.KimT.-M.ChoiB.KimY.KimE. (2021). Invasive Lactuca serriola seeds contain endophytic bacteria that contribute to drought tolerance.Sci. Rep.11:13307. 10.1038/s41598-021-92706-x
42
Johnston-MonjeD.RaizadaM. N. (2011). Conservation and diversity of seed associated endophytes in zea across boundaries of evolution, ethnography and ecology.PLoS One6:e20396. 10.1371/journal.pone.0020396
43
JuX.GuanJ.FanH.AnQ.WuR.ZhengJ. (2021). “Application of GEE in cotton monitoring of the 7th division of Xinjiang Production and Construction Corps,” in Proceedings of the 2021 9th International Conference on Agro-Geoinformatics (Agro-Geoinformatics), Shenzhen, 10.1109/Agro-Geoinformatics50104.2021.9530362
44
KagaH.ManoH.TanakaF.WatanabeA.KanekoS.MorisakiH. (2009). Rice seeds as sources of endophytic bacteria.Microb. Environ.24154–162. 10.1264/jsme2.ME09113
45
KnightsD.KuczynskiJ.CharlsonE. S.ZaneveldJ.MozerM. C.CollmanR. G.et al (2011). Bayesian community-wide culture-independent microbial source tracking.Nat. Methods8761–763. 10.1038/nmeth.1650
46
LahtiL.ShettyS. (2019). Tools For Microbiome Analysis in R.Groningen: University of Groningen.
47
LiuY.ZuoS.XuL.ZouY.SongW. (2012). Study on diversity of endophytic bacterial communities in seeds of hybrid maize and their parental lines.Arch. Microbiol.1941001–1012. 10.1007/s00203-012-0836-8
48
López-LópezA.RogelM. A.Ormeño-OrrilloE.Martínez-RomeroJ.Martínez-RomeroE. (2010). Phaseolus vulgaris seed-borne endophytic community with novel bacterial species such as Rhizobium endophyticum sp. nov.Syst. Appl. Microbiol.33322–327. 10.1016/j.syapm.2010.07.005
49
LoucaS.ParfreyL. W.DoebeliM. (2016). Decoupling function and taxonomy in the global ocean microbiome.Science3531272–1277. 10.1126/science.aaf4507
50
MartinM. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J.17:10. 10.14806/ej.17.1.200
51
MalviyaD.VarmaA.SinghU. B.SinghS.SinghH. V.SaxenaA. K. (2022). Sulfur-oxidizing bacteria from coal mine enhance sulfur nutrition in Pigeonpea (Cajanus cajan L.).Front. Environ. Sci.10:932402. 10.3389/fenvs.2022.932402
52
ManoH.TanakaF.WatanabeA.KagaH.OkunishiS.MorisakiH. (2006). Culturable Surface and Endophytic Bacterial Flora of the Maturing Seeds of Rice Plants (Oryza sativa) Cultivated in a Paddy Field.Microb. Environ.2186–100. 10.1264/jsme2.21.86
53
MastrettaC.TaghaviS.Van Der LelieD.MengoniA.GalardiF.GonnelliC.et al (2009). Endophytic bacteria from seeds of Nicotiana tabacum can reduce cadmium phytotoxicity.Int. J. Phytoremed.11251–267. 10.1080/15226510802432678
54
MitterB.PfaffenbichlerN.FlavellR.CompantS.AntonielliL.PetricA.et al (2017). A new approach to modify plant microbiomes and traits by introducing beneficial bacteria at flowering into progeny seeds.Front. Microbiol.8:11. 10.3389/fmicb.2017.00011
55
MoreiraL. R. S.FilhoE. X. F. (2016). Insights into the mechanism of enzymatic hydrolysis of xylan.Appl. Microbiol. Biotechnol.1005205–5214. 10.1007/s00253-016-7555-z
56
MowlickS.TakeharaT.KakuN.UekiK.UekiA. (2013). Proliferation of diversified clostridial species during biological soil disinfestation incorporated with plant biomass under various conditions.Appl. Microbiol. Biotechnol.978365–8379. 10.1007/s00253-012-4532-z
57
NakajimaT.KawanoY.OhtsuI.Maruyuama-NakashitaA.AllahhamA.SatoM.et al (2019). Effects of thiosulfate as a sulfur source on plant growth, metabolites accumulation and gene expression in arabidopsis and rice.Plant Cell Physiol.601683–1701. 10.1093/pcp/pcz082
58
NelsonE. B. (2004). Microbial dynamics and interactions in the spermosphere.Annu. Rev. Phytopathol.42271–309. 10.1146/annurev.phyto.42.121603.131041
59
NelsonE. B. (2018). The seed microbiome: Origins, interactions, and impacts.Plant Soil4227–34. 10.1007/s11104-017-3289-7
60
PanX.RaaijmakersJ. M.CarriónV. J. (2023). Importance of Bacteroidetes in host–microbe interactions and ecosystem functioning.Trends Microbiol.31959–971. 10.1016/j.tim.2023.03.018
61
Pitumpe ArachchigeP. S.RossoL. H. M.HanselF. D.RamundoB.TorresA. R.AsebedoR.et al (2020). Temporal biological nitrogen fixation pattern in soybean inoculated with Bradyrhizobium.Agrosyst. Geosci. Env.3:e20079. 10.1002/agg2.20079
62
PredeepaR.KumarR.AbrahamG. C.SubramanianT. S. (2022). Salt Tolerant Gossypium hirsutum L. Cultivars are Mostly Long-Staple Length Cultivars.Preprint10.21203/rs.3.rs-1210602/v1
63
QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools.Nucleic Acids Res.41D590–D596. 10.1093/nar/gks1219
64
R Core Team (2013). R: A Language and Environment for Statistical Computing.Vienna: R Core Team.
65
RevaO. N.DixeliusC.MeijerJ.PriestF. G. (2004). Taxonomic characterization and plant colonizing abilities of some bacteria related to Bacillus amyloliquefaciens and Bacillus subtilis.FEMS Microbiol. Ecol.48249–259. 10.1016/j.femsec.2004.02.003
66
RezkiS.CampionC.SimoneauP.JacquesM.-A.ShadeA.BarretM. (2018). Assembly of seed-associated microbial communities within and across successive plant generations.Plant Soil42267–79. 10.1007/s11104-017-3451-2
67
RudgersJ. A.AfkhamiM. E.RúaM. A.DavittA. J.HammerS.HuguetV. M. (2009). A fungus among us: broad patterns of endophyte distribution in the grasses.Ecology901531–1539. 10.1890/08-0116.1
68
Sánchez-LópezA. S.ThijsS.BeckersB.González-ChávezM. C.WeyensN.Carrillo-GonzálezR.et al (2018). Community structure and diversity of endophytic bacteria in seeds of three consecutive generations of Crotalaria pumila growing on metal mine residues.Plant Soil42251–66. 10.1007/s11104-017-3176-2
69
SchmidtJ. E.KentA. D.BrissonV. L.GaudinA. C. M. (2019). Agricultural management and plant selection interactively affect rhizosphere microbial community structure and nitrogen cycling.Microbiome7:146. 10.1186/s40168-019-0756-9
70
ShadeA.HandelsmanJ. (2012). Beyond the Venn diagram: the hunt for a core microbiome: The hunt for a core microbiome.Environ. Microbiol.144–12. 10.1111/j.1462-2920.2011.02585.x
71
ShadeA.JacquesM.-A.BarretM. (2017). Ecological patterns of seed microbiome diversity, transmission, and assembly.Curr. Opin. Microbiol.3715–22. 10.1016/j.mib.2017.03.010
72
STAMP (2015). Bioinformatics Software. Available online at: https://beikolab.cs.dal.ca/software/STAMP(accessed November 20, 2023).
73
SussantiD.KhairulU.RestiZ. (2023). The antagonistic potential of Bacillus spp. against Ralstonia solanacearum causing bacterial wilt disease in ginger (Zingiber officinale).IOP Conf. Ser. Earth Environ. Sci.1160:012044. 10.1088/1755-1315/1160/1/012044
74
TauléC.MarequeC.BarloccoC.HackembruchF.ReisV. M.SicardiM.et al (2012). The contribution of nitrogen fixation to sugarcane (Saccharum officinarum L.), and the identification and characterization of part of the associated diazotrophic bacterial community.Plant Soil35635–49. 10.1007/s11104-011-1023-4
75
TkalecV.MahnicA.GselmanP.RupnikM. (2022). Analysis of seed-associated bacteria and fungi on staple crops using the cultivation and metagenomic approaches.Folia Microbiol.67351–361. 10.1007/s12223-022-00958-5
76
TruyensS.WeyensN.CuypersA.VangronsveldJ. (2013). Changes in the population of seed bacteria of transgenerationally Cd-exposed Arabidopsis thaliana.Plant Biol. J.15971–981. 10.1111/j.1438-8677.2012.00711.x
77
TruyensS.WeyensN.CuypersA.VangronsveldJ. (2014). Bacterial seed endophytes: genera, vertical transmission and interaction with plants.Environ. Microbiol. Rep.740–50.
78
VandenkoornhuyseP.QuaiserA.DuhamelM.Le VanA.DufresneA. (2015). The importance of the microbiome of the plant holobiont.New Phytol.2061196–1206. 10.1111/nph.13312
79
WalitangD. I.KimC.JeonS.KangY.SaT. (2019). Conservation and transmission of seed bacterial endophytes across generations following crossbreeding and repeated inbreeding of rice at different geographic locations.MicrobiologyOpen8e00662. 10.1002/mbo3.662
80
WaltersonA. M.StavrinidesJ. (2015). Pantoea: insights into a highly versatile and diverse genus within the Enterobacteriaceae.FEMS Microbiol. Rev.39968–984. 10.1093/femsre/fuv027
81
WangF.LiX.LiY.HanJ.ChenY.ZengJ.et al (2021). Arabidopsis P4 ATPase-mediated cell detoxification confers resistance to Fusarium graminearum and Verticillium dahliae.Nat. Commun.126426. 10.1038/s41467-021-26727-5
82
WangH.QuY.LiD.AmbuchiJ. J.HeW.ZhouX.et al (2016). Cascade degradation of organic matters in brewery wastewater using a continuous stirred microbial electrochemical reactor and analysis of microbial communities.Sci. Rep.6:27023. 10.1038/srep27023
83
WangQ.GarrityG. M.TiedjeJ. M.ColeJ. R. (2007). Naïve bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.Appl. Environ. Microbiol.735261–5267. 10.1128/AEM.00062-07
84
WardT.LarsonJ.MeulemansJ.HillmannB.LynchJ.SidiropoulosD.et al (2017). BugBase predicts organism-level microbiome phenotypes.Bioinformatics10.1101/133462[Epub ahead of print].
85
WicaksonoW. A.MoraufC.MüllerH.AbdelfattahA.DonatC.BergG. (2023). The mature phyllosphere microbiome of grapevine is associated with resistance against Plasmopara viticola.Front. Microbiol.14:1149307. 10.3389/fmicb.2023.1149307
86
WickhamH. (2016). ggplot2: Elegant Graphics for Data Analysis.Cham: Springer International Publishing.
87
ZhangJ.ZhangC.YangJ.ZhangR.GaoJ.ZhaoX.et al (2019). Insights into endophytic bacterial community structures of seeds among various Oryza sativa L. rice genotypes.J. Plant Growth Regul.3893–102. 10.1007/s00344-018-9812-0
88
ZhangX.MaY.-N.WangX.LiaoK.HeS.ZhaoX.et al (2022). Dynamics of rice microbiomes reveal core vertically transmitted seed endophytes.Microbiome10:216. 10.1186/s40168-022-01422-9
89
ZhangZ.HuangJ.YaoY.PetersG.MacdonaldB.La RosaA. D.et al (2023). Environmental impacts of cotton and opportunities for improvement.Nat. Rev. Earth Environ.4703–715. 10.1038/s43017-023-00476-z
90
ZhaoX. (2019). Analysis of endophytic bacterial diversity in rice seeds and definition of core microbiome. Master’s thesis. Beijing: Chinese Academy of Agricultural Sciences.
91
ZhaoY.GuanD.LiuX.GaoG.-F.MengF.LiuB.et al (2022). Profound change in soil microbial assembly process and co-occurrence pattern in co-inoculation of Bradyrhizobium japonicum 5038 and Bacillus aryabhattai MB35-5 on soybean.Front. Microbiol.13:846359. 10.3389/fmicb.2022.846359
Summary
Keywords
high-throughput sequencing, cotton microbiome, vertical transfer, seed endophytic microbiome, functional traits
Citation
Wu C, Zhang X, Fan Y, Ye J, Dong L, Wang Y, Ren Y, Yong H, Liu R and Wang A (2024) Vertical transfer and functional characterization of cotton seed core microbiome. Front. Microbiol. 14:1323342. doi: 10.3389/fmicb.2023.1323342
Received
17 October 2023
Accepted
22 December 2023
Published
09 January 2024
Volume
14 - 2023
Edited by
Bernardo González, Adolfo Ibáñez University, Chile
Reviewed by
Md. Rashedul Islam, Agriculture and Agri-Food Canada (AAFC), Canada
Leonardo Magalhães Cruz, Federal University of Paraná, Brazil
Updates
Copyright
© 2024 Wu, Zhang, Fan, Ye, Dong, Wang, Ren, Yong, Liu and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Aiying Wang, way-sh@126.com
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.