RESUMO
Plants have powerful defense mechanisms and extensive immune receptor repertoires, yet crop monocultures are prone to epidemic diseases. Rice (Oryza sativa) is susceptible to many diseases, such as rice blast caused by Magnaporthe oryzae. Varietal resistance of rice to blast relies on intracellular nucleotide binding, leucine-rich repeat (NLR) receptors that recognize specific pathogen molecules and trigger immune responses. In the Yuanyang terraces in southwest China, rice landraces rarely show severe losses to disease whereas commercial inbred lines show pronounced field susceptibility. Here, we investigate within-landrace NLR sequence diversity of nine rice landraces and eleven modern varieties using complexity reduction techniques. We find that NLRs display high sequence diversity in landraces, consistent with balancing selection, and that balancing selection at NLRs is more pervasive in landraces than modern varieties. Notably, modern varieties lack many ancient NLR haplotypes that are retained in some landraces. Our study emphasizes the value of standing genetic variation that is maintained in farmer landraces as a resource to make modern crops and agroecosystems less prone to disease. The conservation of landraces is, therefore, crucial for ensuring food security in the face of dynamic biotic and abiotic threats.
RESUMO
Magnaporthe AVRs and ToxB-like (MAX) effectors constitute a family of secreted virulence proteins in the fungus Pyricularia oryzae (syn. Magnaporthe oryzae), which causes blast disease on numerous cereals and grasses. In spite of high sequence divergence, MAX effectors share a common fold characterized by a ß-sandwich core stabilized by a conserved disulfide bond. In this study, we investigated the structural landscape and diversity within the MAX effector repertoire of P. oryzae. Combining experimental protein structure determination and in silico structure modeling we validated the presence of the conserved MAX effector core domain in 77 out of 94 groups of orthologs (OG) identified in a previous population genomic study. Four novel MAX effector structures determined by NMR were in remarkably good agreement with AlphaFold2 (AF2) predictions. Based on the comparison of the AF2-generated 3D models we propose a classification of the MAX effectors superfamily in 20 structural groups that vary in the canonical MAX fold, disulfide bond patterns, and additional secondary structures in N- and C-terminal extensions. About one-third of the MAX family members remain singletons, without strong structural relationship to other MAX effectors. Analysis of the surface properties of the AF2 MAX models also highlights the high variability within the MAX family at the structural level, potentially reflecting the wide diversity of their virulence functions and host targets.
Assuntos
Ascomicetos , Proteínas Fúngicas , Doenças das Plantas , Proteínas Fúngicas/química , Proteínas Fúngicas/metabolismo , Proteínas Fúngicas/genética , Ascomicetos/genética , Ascomicetos/patogenicidade , Ascomicetos/metabolismo , Doenças das Plantas/microbiologia , Modelos Moleculares , Conformação Proteica , Virulência , Fatores de Virulência/genética , Fatores de Virulência/química , Fatores de Virulência/metabolismo , Sequência de AminoácidosRESUMO
Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.
Assuntos
Ascomicetos , Hordeum , Doenças das Plantas , Hordeum/microbiologia , Doenças das Plantas/microbiologia , França , Ascomicetos/genética , Interações Hospedeiro-Patógeno/genética , Filogenia , VernalizaçãoRESUMO
Pyricularia oryzae (syn. Magnaporthe oryzae), is a filamentous ascomycete that causes a major disease called blast on cereal crops, as well as on a wide variety of wild and cultivated grasses. Blast diseases have a tremendous impact worldwide particularly on rice and on wheat, where the disease emerged in South America in the 1980s, before spreading to Asia and Africa. Its economic importance, coupled with its amenability to molecular and genetic manipulation, have inspired extensive research efforts aiming at understanding its biology and evolution. In the past 40 years, this plant-pathogenic fungus has emerged as a major model in molecular plant-microbe interactions. In this review, we focus on the clarification of the taxonomy and genetic structure of the species and its host range determinants. We also discuss recent molecular studies deciphering its lifecycle. TAXONOMY: Kingdom: Fungi, phylum: Ascomycota, sub-phylum: Pezizomycotina, class: Sordariomycetes, order: Magnaporthales, family: Pyriculariaceae, genus: Pyricularia. HOST RANGE: P. oryzae has the ability to infect a wide range of Poaceae. It is structured into different host-specialized lineages that are each associated with a few host plant genera. The fungus is best known to cause tremendous damage to rice crops, but it can also attack other economically important crops such as wheat, maize, barley, and finger millet. DISEASE SYMPTOMS: P. oryzae can cause necrotic lesions or bleaching on all aerial parts of its host plants, including leaf blades, sheaths, and inflorescences (panicles, spikes, and seeds). Characteristic symptoms on leaves are diamond-shaped silver lesions that often have a brown margin and whose appearance is influenced by numerous factors such as the plant genotype and environmental conditions. USEFUL WEBSITES Resources URL Genomic data repositories http://genome.jouy.inra.fr/gemo/ Genomic data repositories http://openriceblast.org/ Genomic data repositories http://openwheatblast.net/ Genome browser for fungi (including P. oryzae) http://fungi.ensembl.org/index.html Comparative genomics database https://mycocosm.jgi.doe.gov/mycocosm/home T-DNA mutant database http://atmt.snu.kr/ T-DNA mutant database http://www.phi-base.org/ SNP and expression data https://fungidb.org/fungidb/app/.
Assuntos
Ascomicetos , Hordeum , Ascomicetos/genética , Produtos Agrícolas , TriticumRESUMO
Fungal pathogens cause devastating disease in crops. Understanding the evolutionary origin of pathogens is essential to the prediction of future disease emergence and the potential of pathogens to disperse. The fungus Pyrenophora teres f. teres causes net form net blotch (NFNB), an economically significant disease of barley. In this study, we have used 104 P. teres f. teres genomes from four continents to explore the population structure and demographic history of the fungal pathogen. We showed that P. teres f. teres is structured into populations that tend to be geographically restricted to different regions. Using Multiple Sequentially Markovian Coalescent and machine learning approaches we demonstrated that the demographic history of the pathogen correlates with the history of barley, highlighting the importance of human migration and trade in spreading the pathogen. Exploring signatures of natural selection, we identified several population-specific selective sweeps that colocalized with genomic regions enriched in putative virulence genes, and loci previously identified as determinants of virulence specificities by quantitative trait locus analyses. This reflects rapid adaptation to local hosts and environmental conditions of P. teres f. teres as it spread with barley. Our research highlights how human activities can contribute to the spread of pathogens that significantly impact the productivity of field crops.
Assuntos
Ascomicetos , Hordeum , Humanos , Hordeum/genética , Hordeum/microbiologia , Domesticação , Ascomicetos/genética , Doenças das Plantas/genética , Doenças das Plantas/microbiologia , Locos de Características Quantitativas/genéticaRESUMO
Neurospora crassa clade A is a model system for genetics, biochemistry, molecular biology, and experimental evolution. Here, we present the draft genome sequences of four isolates of N. crassa clade B. These data represent a valuable resource to investigate the population biology and evolutionary history of N. crassa sensu lato.
RESUMO
The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.
Assuntos
Genômica , Sordariales , Humanos , Filogenia , Genômica/métodos , Genoma , Sordariales/genética , Sequência de Bases , Evolução MolecularRESUMO
Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.
Assuntos
Aminoácidos , Ascomicetos , Virulência/genética , Sequência de Aminoácidos , Ascomicetos/genéticaRESUMO
Traditional agrosystems, where humans, crops and microbes have coevolved over long periods, can serve as models to understand the ecoevolutionary determinants of disease dynamics and help the engineering of durably resistant agrosystems. Here, we investigated the genetic and phenotypic relationship between rice (Oryza sativa) landraces and their rice blast pathogen (Pyricularia oryzae) in the traditional Yuanyang terraces of flooded rice paddies in China, where rice landraces have been grown and bred over centuries without significant disease outbreaks. Analyses of genetic subdivision revealed that indica rice plants clustered according to landrace names. Three new diverse lineages of rice blast specific to the Yuanyang terraces coexisted with lineages previously detected at the worldwide scale. Population subdivision in the pathogen population did not mirror pattern of population subdivision in the host. Measuring the pathogenicity of rice blast isolates on landraces revealed generalist life history traits. Our results suggest that the implementation of disease control strategies based on the emergence or maintenance of a generalist lifestyle in pathogens may sustainably reduce the burden of disease in crops.
Assuntos
Variação Genética , Oryza , Humanos , Oryza/genética , Melhoramento Vegetal , Produtos Agrícolas , China , Doenças das Plantas/genética , Doenças das Plantas/microbiologiaRESUMO
Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.
Assuntos
Ascomicetos , Sordariales , Genes Fúngicos Tipo Acasalamento/genética , Reprodução/genética , Ascomicetos/genética , Sordariales/genética , Recombinação Genética/genética , EsporosRESUMO
Invasions by fungal plant pathogens pose a significant threat to the health of agricultural ecosystems. Despite limited standing genetic variation, many invasive fungal species can adapt and spread rapidly, resulting in significant losses to crop yields. Here, we report on the population genomics of Colletotrichum truncatum, a polyphagous pathogen that can infect more than 460 plant species, and an invasive pathogen of soybean in Brazil. We study the whole-genome sequences of 18 isolates representing 10 fields from two major regions of soybean production. We show that Brazilian C. truncatum is subdivided into three phylogenetically distinct lineages that exchange genetic variation through hybridization. Introgression affects 2%-30% of the nucleotides of genomes and varies widely between the lineages. We find that introgressed regions comprise secreted protein-encoding genes, suggesting possible co-evolutionary targets for selection in those regions. We highlight the inherent vulnerability of genetically uniform crops in the agro-ecological environment, particularly when faced with pathogens that can take full advantage of the opportunities offered by an increasingly globalized world. Finally, we discuss "the means, motive and opportunity" of fungal pathogens and how they can become invasive species of crops. We call for more population genomic studies because such analyses can help identify geographical areas and pathogens that pose a risk, thereby helping to inform control strategies to better protect crops in the future.
Assuntos
Ecossistema , Introgressão Genética , Doenças das Plantas/genética , Doenças das Plantas/microbiologia , Evolução Biológica , Glycine max/genética , Glycine max/microbiologiaRESUMO
[This corrects the article DOI: 10.1371/journal.ppat.1010687.].
RESUMO
Many species of fungal plant pathogens coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. The rice blast fungus Pyricularia oryzae is a widespread model plant pathogen displaying population subdivision. However, most studies of natural variation in P. oryzae have been limited in genomic or geographic resolution, and host adaptation is the only factor that has been investigated extensively as a contributor to population subdivision. In an effort to complement previous studies, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus covering a broad geographical range. Using single-nucleotide polymorphism genotyping data for 886 isolates sampled from 152 sites in 51 countries, we showed that population subdivision of P. oryzae in one recombining and three clonal lineages with broad distributions persisted with deeper sampling. We also extended previous findings by showing further population subdivision of the recombining lineage into one international and three Asian clusters, and by providing evidence that the three clonal lineages of P. oryzae were found in areas with different prevailing environmental conditions, indicating niche separation. Pathogenicity tests and bioinformatic analyses using an extended set of isolates and rice varieties indicated that partial specialization to rice subgroups contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong additional barriers to gene flow between clonal lineages. Our results demonstrate that the spread of a fungal pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct, widely distributed, lineages.
Assuntos
Magnaporthe , Oryza , Ascomicetos , Variação Genética , Magnaporthe/genética , Oryza/microbiologia , Doenças das Plantas/microbiologiaRESUMO
Effectors are small and very diverse proteins secreted by fungi and translocated in plant cells during infection. Among them, MAX effectors (for Magnaporthe Avrs and ToxB) were identified as a family of effectors that share an identical fold topology despite having highly divergent sequences. They are mostly secreted by ascomycetes from the Magnaporthe genus, a fungus that causes the rice blast, a plant disease leading to huge crop losses. As rice is the first source of calories in many countries, especially in Asia and Africa, this constitutes a threat for world food security. Hence, a better understanding of these effectors, including structural and functional characterization, constitutes a strategic milestone in the fight against phytopathogen fungi and may give clues for the development of resistant varieties of rice. We report here the near complete 1H, 15 N and 13C NMR resonance assignment of three new putative MAX effectors (MAX47, MAX60 and MAX67). Secondary structure determination using TALOS-N and CSI.3 demonstrates a high content of ß-strands in all the three proteins, in agreement with the canonic ß-sandwich structure of MAX effectors. This preliminary study provides foundations for further structural characterization, that will help in turn to improve sequence predictions of other MAX effectors through data mining.
Assuntos
Ascomicetos , Magnaporthe , Oryza , Ascomicetos/metabolismo , Proteínas Fúngicas/química , Magnaporthe/metabolismo , Ressonância Magnética Nuclear Biomolecular , Oryza/metabolismo , Oryza/microbiologiaRESUMO
Hybridization and adaptation to new hosts are important mechanisms of fungal disease emergence. Evaluating the risk of emergence of hybrids with enhanced virulence is then key to develop sustainable crop disease management. We evaluated this risk in Venturia inaequalis, the fungus responsible for the common and serious scab disease on Rosaceae hosts, including apple, pyracantha, and loquat. Field isolates from these three hosts and progenies obtained from five crosses between formae speciales isolates collected from pyracantha (f. sp. pyracantha) and apple (f. sp. pomi) were tested for their pathogenicity on the three hosts. We confirmed a strict host specificity between isolates from apple and pyracantha and showed that most isolates were able to cause disease on loquat. None of the 251 progeny obtained from five crosses between V. inaequalis f. sp. pyracantha and V. inaequalis f. sp. pomi could infect apple. If confirmed on more crosses, the inability of the hybrids to infect apple could lead to a novel biocontrol strategy based on a sexual hijacking of V. inaequalis f. sp. pomi by a massive introduction of V. inaequalis f. sp. pyracantha in apple orchards. This strategy, analogous to the sterile insect approach, could lead to the collapse of the population size of V. inaequalis and dramatically reduce the use of chemicals in orchards.
Assuntos
Ascomicetos , Malus , Ascomicetos/genética , Fungos do Gênero Venturia , Malus/microbiologia , Doenças das Plantas/microbiologia , Doenças das Plantas/prevenção & controleRESUMO
Many fungal plant pathogens encompass multiple populations specialized on different plant species. Understanding the factors underlying pathogen adaptation to their hosts is a major challenge of evolutionary microbiology, and it should help to prevent the emergence of new specialized pathogens on novel hosts. Previous studies have shown that French populations of the gray mold pathogen Botrytis cinerea parasitizing tomato and grapevine are differentiated from each other, and have higher aggressiveness on their host of origin than on other hosts, indicating some degree of host specialization in this polyphagous pathogen. Here, we aimed at identifying the genomic features underlying the specialization of B. cinerea populations to tomato and grapevine. Based on whole genome sequences of 32 isolates, we confirmed the subdivision of B. cinerea pathogens into two genetic clusters on grapevine and another, single cluster on tomato. Levels of genetic variation in the different clusters were similar, suggesting that the tomato-specific cluster has not recently emerged following a bottleneck. Using genome scans for selective sweeps and divergent selection, tests of positive selection based on polymorphism and divergence at synonymous and nonsynonymous sites, and analyses of presence and absence variation, we identified several candidate genes that represent possible determinants of host specialization in the tomato-associated population. This work deepens our understanding of the genomic changes underlying the specialization of fungal pathogen populations.
Assuntos
Botrytis , Solanum lycopersicum , Botrytis/genética , França , Genética Populacional , Solanum lycopersicum/microbiologia , Metagenômica , Doenças das Plantas/microbiologiaRESUMO
Population genetics has been a key discipline in phytopathology for many years. The recent rise in cost-effective, high-throughput DNA sequencing technologies, allows sequencing of dozens, if not hundreds of specimens, turning population genetics into population genomics and opening up new, exciting opportunities as described in this Focus Issue. Without the limitations of genetic markers and the availability of whole or near whole-genome data, population genomics can give new insights into the biology, evolution and adaptation, and dissemination patterns of plant-associated microbes.
Assuntos
Metagenômica , Doenças das Plantas , Genética Populacional , Genômica , FilogeniaRESUMO
Blast disease is a notorious fungal disease leading to dramatic yield losses on major food crops such as rice and wheat. The causal agent, Pyricularia oryzae, encompasses different lineages, each having a different host range. Host shifts are suspected to have occurred in this species from Setaria spp. to rice and from Lolium spp. to wheat. The emergence of blast disease on maize in Iran was observed for the first time in the north of the country in 2012. We later identified blast disease in two additional regions of Iran: Gilan in 2013 and Golestan in 2016. Epidemics on the weed barnyard grass (Echinochloa spp.) were also observed in the same maize fields. Here, we showed that P. oryzae is the causal agent of this disease on both hosts. Pathogenicity assays in the greenhouse revealed that strains from maize can infect barnyard grass and conversely. However, genotyping with simple sequence repeat markers and comparative genomics showed that strains causing field epidemics on maize and on barnyard grass are different, although they belong to the same previously undescribed clade of P. oryzae. Phylogenetic analyses including these strains and a maize strain collected in Gabon in 1985 revealed two independent host-range expansion events from barnyard grass to maize. Comparative genomics between maize and barnyard grass strains revealed the presence or absence of five candidate genes associated with host specificity on maize, with the deletion of a small genomic region possibly responsible for adaptation to maize. This recent emergence of P. oryzae on maize provides a case study to understand host range expansion. Epidemics on maize raise concerns about potential yield losses on this crop in Iran and potential geographic expansion of the disease.
Assuntos
Echinochloa , Magnaporthe , Oryza , Ascomicetos , Irã (Geográfico) , Filogenia , Doenças das Plantas , Zea maysRESUMO
The filamentous fungus Podospora anserina has been used as a model organism for more than 100 years and has proved to be an invaluable resource in numerous areas of research. Throughout this period, P. anserina has been embroiled in a number of taxonomic controversies regarding the proper name under which it should be called. The most recent taxonomic treatment proposed to change the name of this important species to Triangularia anserina. The results of past name changes of this species indicate that the broader research community is unlikely to accept this change, which will lead to nomenclatural instability and confusion in literature. Here, we review the phylogeny of the species closely related to P. anserina and provide evidence that currently available marker information is insufficient to resolve the relationships amongst many of the lineages. We argue that it is not only premature to propose a new name for P. anserina based on current data, but also that every effort should be made to retain P. anserina as the current name to ensure stability and to minimise confusion in scientific literature. Therefore, we synonymise Triangularia with Podospora and suggest that either the type species of Podospora be moved to P. anserina from P. fimiseda or that all species within the Podosporaceae be placed in the genus Podospora.
RESUMO
The Honghe Hani rice terraces system (HHRTS) is a traditional rice cultivation system where Hani people cultivate remarkably diverse rice varieties. Recent introductions of modern rice varieties to the HHRTS have significantly increased the severity of rice diseases within the terraces. Here, we determine the impacts of these recent introductions on the composition of the rice-associated microbial communities. We confirm that the HHRTS contains a range of both traditional HHRTS landraces and introduced modern rice varieties and find differences between the microbial communities of these two groups. However, this introduction of modern rice varieties has not strongly impacted the overall diversity of the HHRTS rice microbial community. Furthermore, we find that the rice varieties (i.e. groups of closely related genotypes) have significantly structured the rice microbial community composition (accounting for 15%-22% of the variance) and that the core microbial community of HHRTS rice plants represents less than 3.3% of all the microbial taxa identified. Collectively, our study suggests a highly diverse HHRTS rice holobiont (host with its associated microbes) where the diversity of rice hosts mirrors the diversity of their microbial communities. Further studies will be needed to better determine how such changes might impact the sustainability of the HHRTS.