Symbioses between fungi and bacteria: from mechanisms to impacts on biodiversity.

Symbiotic interactions between fungi and bacteria range from positive to negative. They are ubiquitous in free-living as well as host-associated microbial communities worldwide. Yet, the impact of fungal-bacterial symbioses on the organization and dynamics of microbial communities is uncertain. There are two reasons for this uncertainty: (1) knowledge gaps in the understanding of the genetic mechanisms underpinning fungal-bacterial symbioses and (2) prevailing interpretations of ecological theory that favor antagonistic interactions as drivers stabilizing biological communities despite the existence of models emphasizing contributions of positive interactions. This review synthesizes information on fungal-bacterial symbioses common in the free-living microbial communities of the soil as well as in host-associated polymicrobial biofilms. The interdomain partnerships are considered in the context of the relevant community ecology models, which are discussed critically.

Convergent reductive evolution and host adaptation in Mycoavidus bacterial endosymbionts of Mortierellaceae fungi.

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

Sequencing the Genomes of the First Terrestrial Fungal Lineages: What Have We Learned?

The first genome sequenced of a eukaryotic organism was for Saccharomyces cerevisiae, as reported in 1996, but it was more than 10 years before any of the zygomycete fungi, which are the early-diverging terrestrial fungi currently placed in the phyla Mucoromycota and Zoopagomycota, were sequenced. The genome for Rhizopus delemar was completed in 2008; currently, more than 1000 zygomycete genomes have been sequenced. Genomic data from these early-diverging terrestrial fungi revealed deep phylogenetic separation of the two major clades-primarily plant-associated saprotrophic and mycorrhizal Mucoromycota versus the primarily mycoparasitic or animal-associated parasites and commensals in the Zoopagomycota. Genomic studies provide many valuable insights into how these fungi evolved in response to the challenges of living on land, including adaptations to sensing light and gravity, development of hyphal growth, and co-existence with the first terrestrial plants. Genome sequence data have facilitated studies of genome architecture, including a history of genome duplications and horizontal gene transfer events, distribution and organization of mating type loci, rDNA genes and transposable elements, methylation processes, and genes useful for various industrial applications. Pathogenicity genes and specialized secondary metabolites have also been detected in soil saprobes and pathogenic fungi. Novel endosymbiotic bacteria and viruses have been discovered during several zygomycete genome projects. Overall, genomic information has helped to resolve a plethora of research questions, from the placement of zygomycetes on the evolutionary tree of life and in natural ecosystems, to the applied biotechnological and medical questions.

Molecular Dialogues between Early Divergent Fungi and Bacteria in an Antagonism versus a Mutualism.

Fungal-bacterial symbioses range from antagonisms to mutualisms and remain one of the least understood interdomain interactions despite their ubiquity as well as ecological and medical importance. To build a predictive conceptual framework for understanding interactions between fungi and bacteria in different types of symbioses, we surveyed fungal and bacterial transcriptional responses in the mutualism between Rhizopus microsporus (Rm) (ATCC 52813, host) and its Mycetohabitans (formerly Burkholderia) endobacteria versus the antagonism between a nonhost Rm (ATCC 11559) and Mycetohabitans isolated from the host, at two time points, before and after partner physical contact. We found that bacteria and fungi sensed each other before contact and altered gene expression patterns accordingly. Mycetohabitans did not discriminate between the host and nonhost and engaged a common set of genes encoding known as well as novel symbiosis factors. In contrast, responses of the host versus nonhost to endobacteria were dramatically different, converging on the altered expression of genes involved in cell wall biosynthesis and reactive oxygen species (ROS) metabolism. On the basis of the observed patterns, we formulated a set of hypotheses describing fungal-bacterial interactions and tested some of them. By conducting ROS measurements, we confirmed that nonhost fungi increased production of ROS in response to endobacteria, whereas host fungi quenched their ROS output, suggesting that ROS metabolism contributes to the nonhost resistance to bacterial infection and the host ability to form a mutualism. Overall, our study offers a testable framework of predictions describing interactions of early divergent Mucoromycotina fungi with bacteria.IMPORTANCE Animals and plants interact with microbes by engaging specific surveillance systems, regulatory networks, and response modules that allow for accommodation of mutualists and defense against antagonists. Antimicrobial defense responses are mediated in both animals and plants by innate immunity systems that owe their functional similarities to convergent evolution. Like animals and plants, fungi interact with bacteria. However, the principles governing these relations are only now being discovered. In a study system of host and nonhost fungi interacting with a bacterium isolated from the host, we found that bacteria used a common gene repertoire to engage both partners. In contrast, fungal responses to bacteria differed dramatically between the host and nonhost. These findings suggest that as in animals and plants, the genetic makeup of the fungus determines whether bacterial partners are perceived as mutualists or antagonists and what specific regulatory networks and response modules are initiated during each encounter.

A TAL effector-like protein of an endofungal bacterium increases the stress tolerance and alters the transcriptome of the host.

Symbioses of bacteria with fungi have only recently been described and are poorly understood. In the symbiosis of Mycetohabitans (formerly Burkholderia) rhizoxinica with the fungus Rhizopus microsporus, bacterial type III (T3) secretion is known to be essential. Proteins resembling T3-secreted transcription activator-like (TAL) effectors of plant pathogenic bacteria are encoded in the three sequenced Mycetohabitans spp. genomes. TAL effectors nuclear-localize in plants, where they bind and activate genes important in disease. The Burkholderia TAL-like (Btl) proteins bind DNA but lack the N- and C-terminal regions, in which TAL effectors harbor their T3 and nuclear localization signals, and activation domain. We characterized a Btl protein, Btl19-13, and found that, despite the structural differences, it can be T3-secreted and can nuclear-localize. A btl19-13 gene knockout did not prevent the bacterium from infecting the fungus, but the fungus became less tolerant to cell membrane stress. Btl19-13 did not alter transcription in a plant-based reporter assay, but 15 R. microsporus genes were differentially expressed in comparisons both of the fungus infected with the wild-type bacterium vs. the mutant and with the mutant vs. a complemented strain. Southern blotting revealed btl genes in 14 diverse Mycetohabitans isolates. However, banding patterns and available sequences suggest variation, and the btl19-13 phenotype could not be rescued by a btl gene from a different strain. Our findings support the conclusion that Btl proteins are effectors that act on host DNA and play important but varied or possibly host genotype-specific roles in the M. rhizoxinica-R. microsporus symbiosis.

Narnaviruses: novel players in fungal-bacterial symbioses.

Rhizopus microsporus is an early-diverging fungal species with importance in ecology, agriculture, food production, and public health. Pathogenic strains of R. microsporus harbor an intracellular bacterial symbiont, Mycetohabitans (formerly named Burkholderia). This vertically transmitted bacterial symbiont is responsible for the production of toxins crucial to the pathogenicity of Rhizopus and remarkably also for fungal reproduction. Here we show that R. microsporus can live not only in symbiosis with bacteria but also with two viral members of the genus Narnavirus. Our experiments revealed that both viruses replicated similarly in the growth conditions we tested. Viral copies were affected by the developmental stage of the fungus, the substrate, and the presence or absence of Mycetohabitans. Absolute quantification of narnaviruses in isolated asexual sporangiospores and sexual zygospores indicates their vertical transmission. By curing R. microsporus of its viral and bacterial symbionts and reinfecting bacteria to reestablish symbiosis, we demonstrate that these viruses affect fungal biology. Narnaviruses decrease asexual reproduction, but together with Mycetohabitans, are required for sexual reproductive success. This fungal-bacterial-viral system represents an outstanding model to investigate three-way microbial symbioses and their evolution.

Recombination contributes to population diversification in the polyploid intestinal symbiont Epulopiscium sp. type B.

Epulopiscium sp. type B (Lachnospiraceae) is an exceptionally large, highly polyploid, intestinal symbiont of the coral reef dwelling surgeonfish Naso tonganus. These obligate anaerobes do not form mature endospores and reproduce solely through the production of multiple intracellular offspring. This likely makes them dependent on immediate transfer to a receptive host for dispersal. During reproduction, only a small proportion of Epulopiscium mother-cell DNA is inherited. To explore the impact of this unusual viviparous lifestyle on symbiont population dynamics, we investigated Epulopiscium sp. type B and their fish hosts collected over the course of two decades, at island and reef habitats near Lizard Island, Australia. Using multi-locus sequence analysis, we found that recombination plays an important role in maintaining diversity of these symbionts and yet populations exhibit linkage disequilibrium (LD). Symbiont populations showed spatial but not temporal partitioning. Surgeonfish are long-lived and capable of traveling long distances, yet the population structures of Epulopiscium suggest that adult fish tend to not roam beyond a limited locale. Codiversification analyses and traits of this partnership suggest that while symbionts are obligately dependent on their host, the host has a facultative association with Epulopiscium. We suggest that congression of unlinked markers contributes to LD estimates in this and other recombinant populations of bacteria. The findings here inform our understanding of evolutionary processes within intestinal Lachnospiraceae populations.

Distribution and population structure of endobacteria in arbuscular mycorrhizal fungi at North Atlantic dunes.

Arbuscular mycorrhizal fungi (AMF, Glomeromycotina), in addition to forming symbioses with the majority of land plants, harbor vertically transmitted endosymbiotic bacteria 'Candidatus Glomeribacter gigasporarum' (CaGg) and 'Candidatus Moeniiplasma glomeromycotorum' (CaMg). CaGg is a nonessential mutualist of AMF, whereas the lifestyle of CaMg is unknown. To start unraveling the interactions between AMF and their endosymbionts in nature, we examined diversity and distribution of AMF-associated endobacteria in North Atlantic dunes at Cape Cod. Of nearly 500 foredune AMF isolates successfully genotyped during a systematic study, 94% were classified as Gigasporaceae. Two percent of all AMF spores harbored CaGg, and 88% contained CaMg. CaGg was found only in the Gigasporaceae, whereas CaMg was present in Gigasporaceae, Acaulosporaceae, and Diversisporaceae. Incidence of CaGg across AMF was not affected by any of the environmental parameters measured, whereas distribution of CaMg in one of the fungal hosts was impacted by plant density. CaMg populations associated with AMF individuals displayed high levels of genetic diversity but no evidence of gene flow, suggesting that host physical proximity is not sufficient to facilitate horizontal transmission of CaMg. Finally, in addition to a novel lineage of CaGg, we discovered that AMF likely harbor Burkholderia-related bacteria with close phylogenetic affinity to free-living Burkholderia and endobacteria of other Mucoromycota fungi.

Biology of Fungi and Their Bacterial Endosymbionts.

Heritable symbioses, in which endosymbiotic bacteria (EB) are transmitted vertically between host generations, are an important source of evolutionary novelties. A primary example of such symbioses is the eukaryotic cell with its EB-derived organelles. Recent discoveries suggest that endosymbiosis-related innovations can be also found in associations formed by early divergent fungi in the phylum Mucoromycota with heritable EB from two classes, Betaproteobacteria and Mollicutes. These symbioses exemplify novel types of host-symbiont interactions. Studies of these partnerships fuel theoretical models describing mechanisms that stabilize heritable symbioses, control the rate of molecular evolution, and enable the establishment of mutualisms. Lastly, by altering host phenotypes and metabolism, these associations represent an important instrument for probing the basic biology of the Mucoromycota hosts, which remain one of the least explored filamentous fungi.

Bacterial endosymbionts influence host sexuality and reveal reproductive genes of early divergent fungi.

Many heritable mutualisms, in which beneficial symbionts are transmitted vertically between host generations, originate as antagonisms with parasite dispersal constrained by the host. Only after the parasite gains control over its transmission is the symbiosis expected to transition from antagonism to mutualism. Here, we explore this prediction in the mutualism between the fungus Rhizopus microsporus (Rm, Mucoromycotina) and a beta-proteobacterium Burkholderia, which controls host asexual reproduction. We show that reproductive addiction of Rm to endobacteria extends to mating, and is mediated by the symbiont gaining transcriptional control of the fungal ras2 gene, which encodes a GTPase central to fungal reproductive development. We also discover candidate G-protein-coupled receptors for the perception of trisporic acids, mating pheromones unique to Mucoromycotina. Our results demonstrate that regulating host asexual proliferation and modifying its sexual reproduction are sufficient for the symbiont's control of its own transmission, needed for antagonism-to-mutualism transition in heritable symbioses. These properties establish the Rm-Burkholderia symbiosis as a powerful system for identifying reproductive genes in Mucoromycotina.

'Candidatus Moeniiplasma glomeromycotorum', an endobacterium of arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal fungi (AMF, subphylum Glomeromycotina) are symbionts of most terrestrial plants. They commonly harbour endobacteria of a largely unknown biology, referred to as MRE (Mollicutes/mycoplasma-related endobacteria). Here, we propose to accommodate MRE in the novel genus 'Candidatus Moeniiplasma.' Phylogeny reconstructions based on the 16S rRNA gene sequences cluster 'Ca.Moeniiplasma' with representatives of the class Mollicutes, whereas phylogenies derived from amino acid sequences of 19 genes indicate that it is a discrete lineage sharing ancestry with the members of the family Mycoplasmataceae. Cells of 'Ca.Moeniiplasma' reside directly in the host cytoplasm and have not yet been cultivated. They are coccoid, ~500 nm in diameter, with an electron-dense layer outside the plasma membrane. However, the draft genomes of 'Ca.Moeniiplasma' suggest that this structure is not a Gram-positive cell wall. The evolution of 'Ca.Moeniiplasma' appears to be driven by an ultrarapid rate of mutation accumulation related to the loss of DNA repair mechanisms. Moreover, molecular evolution patterns suggest that, in addition to vertical transmission, 'Ca.Moeniiplasma' is able to transmit horizontally among distinct Glomeromycotina host lineages and exchange genes. On the basis of these unique lifestyle features, the new species 'Candidatus Moeniiplasma glomeromycotorum' is proposed.

Lipid metabolic changes in an early divergent fungus govern the establishment of a mutualistic symbiosis with endobacteria.

The recent accumulation of newly discovered fungal-bacterial mutualisms challenges the paradigm that fungi and bacteria are natural antagonists. To understand the mechanisms that govern the establishment and maintenance over evolutionary time of mutualisms between fungi and bacteria, we studied a symbiosis of the fungus Rhizopus microsporus (Mucoromycotina) and its Burkholderia endobacteria. We found that nonhost R. microsporus, as well as other mucoralean fungi, interact antagonistically with endobacteria derived from the host and are not invaded by them. Comparison of gene expression profiles of host and nonhost fungi during interaction with endobacteria revealed dramatic changes in expression of lipid metabolic genes in the host. Analysis of the host lipidome confirmed that symbiosis establishment was accompanied by specific changes in the fungal lipid profile. Diacylglycerol kinase (DGK) activity was important for these lipid metabolic changes, as its inhibition altered the fungal lipid profile and caused a shift in the host-bacterial interaction into an antagonism. We conclude that adjustments in host lipid metabolism during symbiosis establishment, mediated by DGKs, are required for the mutualistic outcome of the Rhizopus-Burkholderia symbiosis. In addition, the neutral and phospholipid profiles of R. microsporus provide important insights into lipid metabolism in an understudied group of oleaginous Mucoromycotina. Lastly, our study revealed that the DGKs involved in the symbiosis form a previously uncharacterized clade of DGK domain proteins.

Defying Muller's Ratchet: Ancient Heritable Endobacteria Escape Extinction through Retention of Recombination and Genome Plasticity.

Heritable endobacteria, which are transmitted from one host generation to the next, are subjected to evolutionary forces that are different from those experienced by free-living bacteria. In particular, they suffer consequences of Muller's ratchet, a mechanism that leads to extinction of small asexual populations due to fixation of slightly deleterious mutations combined with the random loss of the most-fit genotypes, which cannot be recreated without recombination. Mycoplasma-related endobacteria (MRE) are heritable symbionts of fungi from two ancient lineages, Glomeromycota (arbuscular mycorrhizal fungi) and Mucoromycotina Previous studies revealed that MRE maintain unusually diverse populations inside their hosts and may have been associated with fungi already in the early Paleozoic. Here we show that MRE are vulnerable to genomic degeneration and propose that they defy Muller's ratchet thanks to retention of recombination and genome plasticity. We suggest that other endobacteria may be capable of raising similar defenses against Muller's ratchet.

Nondegenerative Evolution in Ancient Heritable Bacterial Endosymbionts of Fungi.

Bacterial endosymbionts are critical to the existence of many eukaryotes. Among them, vertically transmitted endobacteria are uniquely typified by reduced genomes and molecular evolution rate acceleration relative to free-living taxa. These patterns are attributable to genetic drift-dominated degenerative processes associated with reproductive dependence on the host. The degenerative evolution scenario is well supported in endobacteria with strict vertical transmission, such as essential mutualists of insects. In contrast, heritable endosymbionts that are nonessential to their hosts and engage occasionally in horizontal transmission are expected to display deviations from the degenerative evolution model. To explore evolution patterns in such nonessential endobacteria, we focused on Candidatus Glomeribacter gigasporarum ancient heritable mutualists of fungi. Using a collection of genomes, we estimated in Glomeribacter mutation rate at 2.03 × 10(-9) substitutions per site per year and effective population size at 1.44 × 10(8) Both fall within the range of values observed in free-living bacteria. To assess the ability of Glomeribacter to purge slightly deleterious mutations, we examined genome-wide dN/dS values and distribution patterns. We found that these dN/dS profiles cluster Glomeribacter with free-living bacteria as well as with other nonessential endosymbionts, while distinguishing it from essential heritable mutualists of insects. Finally, our evolutionary simulations revealed that the molecular evolution rate acceleration in Glomeribacter is caused by limited recombination in a largely clonal population rather than by increased fixation of slightly deleterious mutations. Based on these patterns, we propose that genome evolution in Glomeribacter is nondegenerative and exemplifies a departure from the model of degenerative evolution in heritable endosymbionts.

The role of mobile genetic elements in evolutionary longevity of heritable endobacteria.

The movement of mobile genetic elements (MGEs), including bacteriophages, insertion sequence (IS) elements, and integrative and conjugative elements (ICEs) can have profound effects on bacterial evolution by introducing novel genes, or disrupting the existing ones. Obligate endobacteria are a distinctive group of bacteria that reside within the intracellular compartments of their eukaryotic hosts. Many obligate endobacteria are reproductively dependent on their hosts. Vertical transmission, in addition to degenerative genome contraction and loss of MGEs, makes heritable endobacteria vulnerable to Muller's ratchet, a process that jeopardizes evolutionary longevity of small populations. Mycoplasma-related endobacteria (MRE) are ancient heritable endosymbionts of arbuscular mycorrhizal fungi. Their genomes harbour numerous MGEs. To explore the significance of MGEs in the evolution of MRE and other obligate endobacteria, we analyze the impact of transmission mode, recombination, and evolutionary age on the maintenance of MGEs. Furthermore, we discuss the ability of MGEs to act as sites of gene conversion and recombination in endobacterial genomes. We propose that MGEs are important instruments of genome shuffling, contributing to population heterogeneity and evolutionary longevity in heritable obligate endobacteria.

Minimal genomes of mycoplasma-related endobacteria are plastic and contain host-derived genes for sustained life within Glomeromycota.

Arbuscular mycorrhizal fungi (AMF, Glomeromycota) colonize roots of the majority of terrestrial plants. They provide essential minerals to their plant hosts and receive photosynthates in return. All major lineages of AMF harbor endobacteria classified as Mollicutes, and known as mycoplasma-related endobacteria (MRE). Except for their substantial intrahost genetic diversity and ability to transmit vertically, virtually nothing is known about the life history of these endobacteria. To understand MRE biology, we sequenced metagenomes of three MRE populations, each associated with divergent AMF hosts. We found that each AMF species harbored a genetically distinct group of MRE. Despite vertical transmission, all MRE populations showed extensive chromosomal rearrangements, which we attributed to genetic recombination, activity of mobile elements, and a history of plectroviral invasion. The MRE genomes are characterized by a highly reduced gene content, indicating metabolic dependence on the fungal host, with the mechanism of energy production remaining unclear. Several MRE genes encode proteins with domains involved in protein-protein interactions with eukaryotic hosts. In addition, the MRE genomes harbor genes horizontally acquired from AMF. Some of these genes encode small ubiquitin-like modifier (SUMO) proteases specific to the SUMOylation systems of eukaryotes, which MRE likely use to manipulate their fungal host. The extent of MRE genome plasticity and reduction, along with the large number of horizontally acquired host genes, suggests a high degree of adaptation to the fungal host. These features, together with the ubiquity of the MRE-Glomeromycota associations, emphasize the significance of MRE in the biology of Glomeromycota.

Molecular evolution patterns reveal life history features of mycoplasma-related endobacteria associated with arbuscular mycorrhizal fungi.

The mycoplasma-related endobacteria (MRE), representing a recently discovered lineage of Mollicutes, are widely distributed across arbuscular mycorrhizal fungi (AMF, Glomeromycota). AMF colonize roots of most terrestrial plants and improve plant mineral nutrient uptake in return for plant-assimilated carbon. The role of MRE in the biology of their fungal hosts is unknown. To start characterizing this association, we assessed partitioning of MRE genetic diversity within AMF individuals and across the AMF phylogeographic range. We further used molecular evolution patterns to make inferences about MRE codivergence with AMF, their lifestyle and antiquity of the Glomeromycota-MRE association. While we did not detect differentiation between MRE derived from different continents, high levels of diversity were apparent in MRE populations within AMF host individuals. MRE exhibited significant codiversification with AMF over ecological time and the absence of codivergence over evolutionary time. Moreover, genetic recombination was evident in MRE. These patterns indicate that, while MRE transmission is predominantly vertical, their complex intrahost populations are likely generated by horizontal transmission and recombination. Based on predictions of evolutionary theory, we interpreted these observations as a suggestion that MRE may be antagonists of AMF. Finally, we detected a marginally significant signature of codivergence of MRE with Glomeromycota and the Endogone lineage of Mucoromycotina, implying that the symbiosis between MRE and fungi may predate the divergence between these two groups of fungi.

Detection of a novel intracellular microbiome hosted in arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal fungi (AMF) are important members of the plant microbiome. They are obligate biotrophs that colonize the roots of most land plants and enhance host nutrient acquisition. Many AMF themselves harbor endobacteria in their hyphae and spores. Two types of endobacteria are known in Glomeromycota: rod-shaped Gram-negative Candidatus Glomeribacter gigasporarum, CaGg, limited in distribution to members of the Gigasporaceae family, and coccoid Mollicutes-related endobacteria, Mre, widely distributed across different lineages of AMF. The goal of the present study is to investigate the patterns of distribution and coexistence of the two endosymbionts, CaGg and Mre, in spore samples of several strains of Gigaspora margarita. Based on previous observations, we hypothesized that some AMF could host populations of both endobacteria. To test this hypothesis, we performed an extensive investigation of both endosymbionts in G. margarita spores sampled from Cameroonian soils as well as in the Japanese G. margarita MAFF520054 isolate using different approaches (molecular phylotyping, electron microscopy, fluorescence in situ hybridization and quantitative real-time PCR). We found that a single AMF host can harbour both types of endobacteria, with Mre population being more abundant, variable and prone to recombination than the CaGg one. Both endosymbionts seem to retain their genetic and lifestyle peculiarities regardless of whether they colonize the host alone or together. These findings show for the first time that fungi support an intracellular bacterial microbiome, in which distinct types of endobacteria coexist in a single cell.

Ribosomal RNA gene diversity, effective population size, and evolutionary longevity in asexual glomeromycota.

Arbuscular mycorrhizal fungi (phylum Glomeromycota) are among the oldest and most successful symbionts of land plants. With no evidence of sexual reproduction, their evolutionary success is inconsistent with the prediction that asexual taxa are vulnerable to extinction due to accumulation of deleterious mutations. To explore why Glomeromycota defy this prediction, we studied ribosomal RNA (rRNA) gene evolution in the Claroideoglomus lineage and estimated effective population size, N(e) , in C. etunicatum. We found that rRNA genes of these fungi exhibit unusual and complex patterns of molecular evolution. In C. etunicatum, these patterns can be collectively explained by an unexpectedly large N(e) combined with imperfect genome-wide and population-level rRNA gene repeat homogenization. The mutations accumulated in rRNA gene sequences indicate that natural selection is effective at purging deleterious mutations in the Claroideoglomus lineage, which is also consistent with the large N(e) of C. etunicatum. We propose that in the near absence of recombination, asexual reproduction involving massively multinucleate spores typical for Glomeromycota is responsible for the improved efficacy of selection relative to drift. We postulate that large effective population sizes contribute to the evolutionary longevity of Glomeromycota.