A phylogenomic and molecular markers based taxonomic framework for members of the order Entomoplasmatales: proposal for an emended order Mycoplasmatales containing the family Spiroplasmataceae and emended family Mycoplasmataceae comprised of six genera

Detection of Mycoplasma spp. from snakes from five different families

BACKGROUND

Mycoplasmas are an important cause of respiratory diseases in tortoises. In snakes, evidence of mycoplasma infections has been found almost exclusively in pythons. To better understand the occurrence of these bacteria in other snake species, samples submitted for routine testing for respiratory pathogens were also tested for mycoplasma by polymerase chain reaction (PCR). A total of 640 samples (mostly oral swabs) from snakes of 5 different families (Boidae n = 114, Colubridae n = 109, Elapidae n = 34, Pythonidae n = 301 and Viperidae n = 82) were included in the study. A genus-specific PCR (PCR1) developed for the detection of Mycoplasma [Mycoplasmopsis] agassizii and a pan-mycoplasma PCR (PCR2) were used. PCR products were sequenced for validation and phylogenetic analysis was performed. The sampled animals were from various owners and collections, all in human care at the time of sampling. Clinical background information was not provided.

RESULTS

Using PCR1, mycoplasmas were detected in 175 (175/640, 27%) samples (Boidae: 7/114, 6%; Colubridae: 3/109, 3%; Elapidae: 8/34, 24%; Pythonidae: 155/301, 51%; Viperidae: 2/82, 2%). A higher percentage of positive results were obtained using PCR2 (258/640, 40%; Boidae: 9/114, 8%; Colubridae: 25/109, 23%; Elapidae: 19/34, 56%; Pythonidae: 172/301, 57%; Viperidae: 33/82, 40%). The detected bacteria can be divided into at least 6 genetically diverse clusters representing different genera and species based on multiple sequence alignment and phylogenetic analysis.

CONCLUSIONS

These results show that diverse mycoplasmas are found in pythons and other snakes. Further investigations are necessary to evaluate the role of various mycoplasmas in respiratory diseases in snakes.

Discovery of deep-sea coral symbionts from a novel clade of marine bacteria with severely reduced genomes

Microbes perform critical functions in corals, yet most knowledge is derived from the photic zone. Here, we discover two mollicutes that dominate the microbiome of the deep-sea octocoral, Callogorgia delta, and likely reside in the mesoglea. These symbionts are abundant across the host's range, absent in the water, and appear to be rare in sediments. Unlike other mollicutes, they lack all known fermentative capabilities, including glycolysis, and can only generate energy from arginine provided by the coral host. Their genomes feature several mechanisms to interact with foreign DNA, including extensive CRISPR arrays and restriction-modification systems, which may indicate their role in symbiosis. We propose the novel family Oceanoplasmataceae which includes these symbionts and others associated with five marine invertebrate phyla. Its exceptionally broad host range suggests that the diversity of this enigmatic family remains largely undiscovered. Oceanoplasmataceae genomes are the most highly reduced among mollicutes, providing new insight into their reductive evolution and the roles of coral symbionts.

Delineating bacterial genera based on gene content analysis: a case study of the Mycoplasmatales-Entomoplasmatales clade within the class Mollicutes

Genome-based analysis allows for large-scale classification of diverse bacteria and has been widely adopted for delineating species. Unfortunately, for higher taxonomic ranks such as genus, establishing a generally accepted approach based on genome analysis is challenging. While core-genome phylogenies depict the evolutionary relationships among species, determining the correspondence between clades and genera may not be straightforward. For genotypic divergence, the percentage of conserved proteins and genome-wide average amino acid identity are commonly used, but often do not provide a clear threshold for classification. In this work, we investigated the utility of global comparisons and data visualization in identifying clusters of species based on their overall gene content and rationalized that such patterns can be integrated with phylogeny and other information such as phenotypes for improving taxonomy. As a proof of concept, we selected 177 representative genome sequences from the Mycoplasmatales-Entomoplasmatales clade within the class Mollicutes for a case study. We found that the clustering patterns corresponded to the current understanding of these organisms, namely the split into three above-genus groups: Hominis, Pneumoniae and Spiroplasma-Entomoplasmataceae-Mycoides. However, at the genus level, several important issues were found. For example, recent taxonomic revisions that split the Hominis group into three genera and Entomoplasmataceae into five genera are problematic, as those newly described or emended genera lack clear differentiations in gene content from one another. Moreover, several cases of misclassification were identified. These findings demonstrated the utility of this approach and its potential application to other bacteria.

Comparative genomic analysis identifies potential adaptive variation in Mycoplasma ovipneumoniae

Mycoplasma ovipneumoniae is associated with respiratory disease in wild and domestic Caprinae globally, with wide variation in disease outcomes within and between host species. To gain insight into phylogenetic structure and mechanisms of pathogenicity for this bacterial species, we compared M. ovipneumoniae genomes for 99 samples from 6 countries (Australia, Bosnia and Herzegovina, Brazil, China, France and USA) and 4 host species (domestic sheep, domestic goats, bighorn sheep and caribou). Core genome sequences of M. ovipneumoniae assemblies from domestic sheep and goats fell into two well-supported phylogenetic clades that are divergent enough to be considered different bacterial species, consistent with each of these two clades having an evolutionary origin in separate host species. Genome assemblies from bighorn sheep and caribou also fell within these two clades, indicating multiple spillover events, most commonly from domestic sheep. Pangenome analysis indicated a high percentage (91.4 %) of accessory genes (i.e. genes found only in a subset of assemblies) compared to core genes (i.e. genes found in all assemblies), potentially indicating a propensity for this pathogen to adapt to within-host conditions. In addition, many genes related to carbon metabolism, which is a virulence factor for Mycoplasmas, showed evidence for homologous recombination, a potential signature of adaptation. The presence or absence of annotated genes was very similar between sheep and goat clades, with only two annotated genes significantly clade-associated. However, three M. ovipneumoniae genome assemblies from asymptomatic caribou in Alaska formed a highly divergent subclade within the sheep clade that lacked 23 annotated genes compared to other assemblies, and many of these genes had functions related to carbon metabolism. Overall, our results suggest that adaptation of M. ovipneumoniae has involved evolution of carbon metabolism pathways and virulence mechanisms related to those pathways. The genes involved in these pathways, along with other genes identified as potentially involved in virulence in this study, are potential targets for future investigation into a possible genomic basis for the high variation observed in disease outcomes within and between wild and domestic host species.

Study of the influence of tributyrin-supplemented diets on the gut bacterial communities of rainbow trout (Oncorhynchus mykiss)

Dietary supplementation with triglyceride tributyrin (TBT), a butyrate precursor, has been associated with beneficial effects on fish health and improvements in the ability of carnivorous fish to tolerate higher levels of plant-based protein. In this study, we aimed to investigate the effects of a plant-based diet supplemented with TBT on the structural diversity and putative function of the digesta-associated bacterial communities of rainbow trout (Oncorhynchus mykiss). In addition to this, we also assessed the response of fish gut digestive enzyme activities and chyme metabolic profile in response to TBT supplementation. Our results indicated that TBT had no significant effects on the overall fish gut bacterial communities, digestive enzyme activities or metabolic profile when compared with non-supplemented controls. However, a more in-depth analysis into the most abundant taxa showed that diets at the highest TBT concentrations (0.2% and 0.4%) selectively inhibited members of the Enterobacterales order and reduced the relative abundance of a bacterial population related to Klebsiella pneumoniae, a potential fish pathogen. Furthermore, the predicted functional analysis of the bacterial communities indicated that increased levels of TBT were associated with depleted KEGG pathways related to pathogenesis. The specific effects of TBT on gut bacterial communities observed here are intriguing and encourage further studies to investigate the potential of this triglyceride to promote pathogen suppression in the fish gut environment, namely in the context of aquaculture.

AppIndels.com server: a web-based tool for the identification of known taxon-specific conserved signature indels in genome sequences. Validation of its usefulness by predicting the taxonomic affiliation of >700 unclassified strains of Bacillus species

Taxon-specific conserved signature indels (CSIs) in genes/proteins provide reliable molecular markers (synapomorphies) for unambiguous demarcation of taxa of different ranks in molecular terms and for genetic, biochemical and diagnostic studies. Because of their predictive abilities, the shared presence of known taxon-specific CSIs in genome sequences has proven useful for taxonomic purposes. However, the lack of a convenient method for identifying the presence of known CSIs in genome sequences has limited their utility for taxonomic and other studies. We describe here a web-based tool/server (AppIndels.com) that identifies the presence of known and validated CSIs in genome sequences and uses this information for predicting taxonomic affiliation. The utility of this server was tested by using a database of 585 validated CSIs, which included 350 CSIs specific for ≈45 Bacillales genera, with the remaining CSIs being specific for members of the orders Neisseriales, Legionellales and Chlorobiales, family Borreliaceae, and some Pseudomonadaceae species/genera. Using this server, genome sequences were analysed for 721 Bacillus strains of unknown taxonomic affiliation. Results obtained showed that 651 of these genomes contained significant numbers of CSIs specific for the following Bacillales genera/families: Alkalicoccus, 'Alkalihalobacillaceae', Alteribacter, Bacillus Cereus clade, Bacillus Subtilis clade, Caldalkalibacillus, Caldibacillus, Cytobacillus, Ferdinandcohnia, Gottfriedia, Heyndrickxia, Lederbergia, Litchfieldia, Margalitia, Mesobacillus, Metabacillus, Neobacillus, Niallia, Peribacillus, Priestia, Pseudalkalibacillus, Robertmurraya, Rossellomorea, Schinkia, Siminovitchia, Sporosarcina, Sutcliffiella, Weizmannia and Caryophanaceae. Validity of the taxon assignment made by the server was examined by reconstructing phylogenomic trees. In these trees, all Bacillus strains for which taxonomic predictions were made correctly branched with the indicated taxa. The unassigned strains likely correspond to taxa for which CSIs are lacking in our database. Results presented here show that the AppIndels server provides a useful new tool for predicting taxonomic affiliation based on shared presence of the taxon-specific CSIs. Some caveats in using this server are discussed.

Mariniplasma anaerobium gen. nov., sp. nov., a novel anaerobic marine mollicute, and proposal of three novel genera to reclassify members of Acholeplasma clusters II-IV

A novel strictly anaerobic chemoorganotrophic bacterium, designated Mahy22^T, was isolated from sulfidic bottom water of a shallow brackish meromictic lake in Japan. Cells of the strain were Gram-stain-negative, non-motile and coccoid in shape with diameters of about 600-800 nm. The temperature range for growth was 15-37 °C, with optimum growth at 30-32 °C. The pH range for growth was pH 6.2-8.9, with optimum growth at pH 7.2-7.4. The strain grew with NaCl concentrations of 5% or below (optimum, 2-3%). Growth of the strain was enhanced by the addition of thiosulfate. The major cellular fatty acids were C_16:0 and anteiso-C_15:0. Respiratory quinones were not detected. The complete genome sequence of strain Mahy22^T possessed a 1 885 846 bp circular chromosome and a 12 782 bp circular genetic element. The G+C content of the genome sequence was 30.1 mol%. Phylogenetic analysis based on the 16S rRNA gene revealed that the novel strain belonged to the family Acholeplasmataceae, class Mollicutes. The closest relative of strain Mahy22^T with a validly published name was Acholeplasma palmae J233^T with a 16S rRNA gene sequence similarity of 90.5%. Based on the results of polyphasic analysis, the name Mariniplasma anaerobium gen. nov., sp. nov. is proposed to accommodate strain Mahy22^T, along with reclassification of some Acholeplasma species into Alteracholeplasma gen. nov., Haploplasma gen. nov. and Paracholeplasma gen. nov.

A High-Throughput Short Sequence Typing Scheme for Serratia marcescens Pure Culture and Environmental DNA

Molecular typing methods are used to characterize the relatedness between bacterial isolates involved in infections. These approaches rely mostly on discrete loci or whole-genome sequencing (WGS) analyses of pure cultures. On the other hand, their application to environmental DNA profiling to evaluate epidemiological relatedness among patients and environments has received less attention. We developed a specific, high-throughput short sequence typing (HiSST) method for the opportunistic human pathogen Serratia marcescens. Genes displaying the highest polymorphism were retrieved from the core genome of 60 S. marcescens strains. Bioinformatics analyses showed that use of only three loci (within bssA, gabR, and dhaM) distinguished strains with a high level of efficiency. This HiSST scheme was applied to an epidemiological survey of S. marcescens in a neonatal intensive care unit (NICU). In a first case study, a strain responsible for an outbreak in the NICU was found in a sink drain of this unit, by using HiSST scheme and confirmed by WGS. The HiSST scheme was also applied to environmental DNA extracted from sink-environment samples. Diversity of S. marcescens was modest, with 11, 6, and 4 different sequence types (ST) of gabR, bssA, and dhaM loci among 19 sink drains, respectively. Epidemiological relationships among sinks were inferred on the basis of pairwise comparisons of ST profiles. Further research aimed at relating ST distribution patterns to environmental features encompassing sink location, utilization, and microbial diversity is needed to improve the surveillance and management of opportunistic pathogens. IMPORTANCE Serratia marcescens is an important opportunistic human pathogen, often multidrug resistant and involved in outbreaks of nosocomial infections in neonatal intensive care units. Here, we propose a quick and user-friendly method to select the best typing scheme for nosocomial outbreaks in relating environmental and clinical sources. This method, named high-throughput short sequence typing (HiSST), allows to distinguish strains and to explore the diversity profile of nonculturable S. marcescens. The application of HiSST profile analysis for environmental DNA offers new possibilities to track opportunistic pathogens, identify their origin, and relate their distribution pattern with environmental features encompassing sink location, utilization, and microbial diversity. Adaptation of the method to other opportunistic pathogens is expected to improve knowledge regarding their ecology, which is of significant interest for epidemiological risk assessment and elaborate outbreak mitigation strategies.

Rapid molecular evolution of Spiroplasma symbionts of Drosophila

Spiroplasma is a genus of Mollicutes whose members include plant pathogens, insect pathogens and endosymbionts of animals. Spiroplasma phenotypes have been repeatedly observed to be spontaneously lost in Drosophila cultures, and several studies have documented a high genomic turnover in Spiroplasma symbionts and plant pathogens. These observations suggest that Spiroplasma evolves quickly in comparison to other insect symbionts. Here, we systematically assess evolutionary rates and patterns of Spiroplasma poulsonii, a natural symbiont of Drosophila. We analysed genomic evolution of sHy within flies, and sMel within in vitro culture over several years. We observed that S. poulsonii substitution rates are among the highest reported for any bacteria, and around two orders of magnitude higher compared with other inherited arthropod endosymbionts. The absence of mismatch repair loci mutS and mutL is conserved across Spiroplasma, and likely contributes to elevated substitution rates. Further, the closely related strains sMel and sHy (>99.5 % sequence identity in shared loci) show extensive structural genomic differences, which potentially indicates a higher degree of host adaptation in sHy, a protective symbiont of Drosophila hydei. Finally, comparison across diverse Spiroplasma lineages confirms previous reports of dynamic evolution of toxins, and identifies loci similar to the male-killing toxin Spaid in several Spiroplasma lineages and other endosymbionts. Overall, our results highlight the peculiar nature of Spiroplasma genome evolution, which may explain unusual features of its evolutionary ecology.

Bacterial Vector-Borne Plant Diseases: Unanswered Questions and Future Directions

Vector-borne plant diseases have significant ecological and economic impacts, affecting farm profitability and forest composition throughout the world. Bacterial vector-borne pathogens have evolved sophisticated strategies to interact with their hemipteran insect vectors and plant hosts. These pathogens reside in plant vascular tissue, and their study represents an excellent opportunity to uncover novel biological mechanisms regulating intracellular pathogenesis and to contribute to the control of some of the world's most invasive emerging diseases. In this perspective, we highlight recent advances and major unanswered questions in the realm of bacterial vector-borne disease, focusing on liberibacters, phytoplasmas, spiroplasmas, and Xylella fastidiosa.

Necessity and rationale for the proposed name changes in the classification of Mollicutes species. Reply to: 'Recommended rejection of the names Malacoplasma gen. nov., Mesomycoplasma gen. nov., Metamycoplasma gen. nov., Metamycoplasmataceae fam. nov., Mycoplasmoidaceae fam. nov., Mycoplasmoidales ord. nov., Mycoplasmoides gen. nov., Mycoplasmopsis gen. nov. [Gupta, Sawnani, Adeolu, Alnajar and Oren 2018] and all proposed species comb. nov. placed therein', by M. Balish et al. (Int J Syst Evol Microbiol, 2019;69:3650-3653)

This response summarizes the highly disordered state of the Mollicutes taxonomy that existed until recently, where most Mollicutes taxa lacked proper circumscriptions and their names were not in accordance with the International Code of Nomenclature of Prokaryotes and illegitimate. We also summarize the comprehensive phylogenomic and comparative genomic studies forming the basis for the proposed changes in the classification of Mollicultes species. Our responses to the concerns raised by Balish et al., show that the proposed taxonomic changes do not violate any essential point of the Code. Instead the proposed name changes rectify numerous taxonomic anomalies that have long plagued the classification of Mollicutes species, leading to a better understanding of their evolutionary relationships and bringing their nomenclature in conformity with the Code.

Phylogenomic Analyses of Members of the Widespread Marine Heterotrophic Genus Pseudovibrio Suggest Distinct Evolutionary Trajectories and a Novel Genus, Polycladidibacter gen. nov

Bacteria belonging to the Pseudovibrio genus are widespread, metabolically versatile, and able to thrive as both free-living and host-associated organisms. Although more than 50 genomes are available, a comprehensive comparative genomics study to resolve taxonomic inconsistencies is currently missing. We analyzed all available genomes and used 552 core genes to perform a robust phylogenomic reconstruction. This in-depth analysis revealed the divergence of two monophyletic basal lineages of strains isolated from polyclad flatworm hosts, namely, Pseudovibrio hongkongensis and Pseudovibrio stylochi These strains have reduced genomes and lack sulfur-related metabolisms and major biosynthetic gene clusters, and their environmental distribution appears to be tightly associated with invertebrate hosts. We showed experimentally that the divergent strains are unable to utilize various sulfur compounds that, in contrast, can be utilized by the type strain Pseudovibrio denitrificans Our analyses suggest that the lineage leading to these two strains has been subject to relaxed purifying selection resulting in great gene loss. Overall genome relatedness indices (OGRI) indicate substantial differences between the divergent strains and the rest of the genus. While 16S rRNA gene analyses do not support the establishment of a different genus for the divergent strains, their substantial genomic, phylogenomic, and physiological differences strongly suggest a divergent evolutionary trajectory and the need for their reclassification. Therefore, we propose the novel genus Polycladidibacter gen. nov.IMPORTANCE The genus Pseudovibrio is commonly associated with marine invertebrates, which are essential for ocean health and marine nutrient cycling. Traditionally, the phylogeny of the genus has been based on 16S rRNA gene analysis. The use of the 16S rRNA gene or any other single marker gene for robust phylogenetic placement has recently been questioned. We used a large set of marker genes from all available Pseudovibrio genomes for in-depth phylogenomic analyses. We identified divergent monophyletic basal lineages within the Pseudovibrio genus, including two strains isolated from polyclad flatworms. These strains showed reduced sulfur metabolism and biosynthesis capacities. The phylogenomic analyses revealed distinct evolutionary trajectories and ecological adaptations that differentiate the divergent strains from the other Pseudovibrio members and suggest that they fall into a novel genus. Our data show the importance of widening the use of phylogenomics for better understanding bacterial physiology, phylogeny, and evolution.

Robust Demarcation of the Family Caryophanaceae (Planococcaceae) and Its Different Genera Including Three Novel Genera Based on Phylogenomics and Highly Specific Molecular Signatures

The family Caryophanaceae/Planococcaceae is a taxonomically heterogeneous assemblage of >100 species classified within 13 genera, many of which are polyphyletic. Exhibiting considerable phylogenetic overlap with other families, primarily Bacillaceae, the evolutionary history of this family, containing the potent mosquitocidal species Lysinibacillus sphaericus, remains incoherent. To develop a reliable phylogenetic and taxonomic framework for the family Caryophanaceae/Planococcaceae and its genera, we report comprehensive phylogenetic and comparative genomic analyses on 124 genome sequences from all available Caryophanaceae/Planococcaceae and representative Bacillaceae species. Phylogenetic trees were constructed based on multiple datasets of proteins including 819 core proteins for this group and 87 conserved Firmicutes proteins. Using the core proteins, pairwise average amino acid identity was also determined. In parallel, comparative analyses on protein sequences from these species have identified 92 unique molecular markers (synapomorphies) consisting of conserved signature indels that are specifically shared by either the entire family Caryophanaceae/Planococcaceae or different monophyletic clades present within this family, enabling their reliable demarcation in molecular terms. Based on multiple lines of investigations, 18 monophyletic clades can be reliably distinguished within the family Caryophanaceae/Planococcaceae based on their phylogenetic affinities and identified molecular signatures. Some of these clades are comprised of species from several polyphyletic genera within this family as well as other families. Based on our results, we are proposing the creation of three novel genera within the family Caryophanaceae/Planococcaceae, namely Metalysinibacillus gen. nov., Metasolibacillus gen. nov., and Metaplanococcus gen. nov., as well as the transfer of 25 misclassified species from the families Caryophanaceae/Planococcaceae and Bacillaceae into these three genera and in Planococcus, Solibacillus, Sporosarcina, and Ureibacillus genera. These amendments establish a coherent taxonomy and evolutionary history for the family Caryophanaceae/Planococcaceae, and the described molecular markers provide novel means for diagnostic, genetic, and biochemical studies. Lastly, we are also proposing a consolidation of the family Planococcaceae within the emended family Caryophanaceae.

Effect of Long-Term Farming Practices on Agricultural Soil Microbiome Members Represented by Metagenomically Assembled Genomes (MAGs) and Their Predicted Plant-Beneficial Genes

To follow the hypothesis that agricultural management practices affect structure and function of the soil microbiome regarding soil health and plant-beneficial traits, high-throughput (HT) metagenome analyses were performed on Chernozem soil samples from a long-term field experiment designated LTE-1 carried out at Bernburg-Strenzfeld (Saxony-Anhalt, Germany). Metagenomic DNA was extracted from soil samples representing the following treatments: (i) plough tillage with standard nitrogen fertilization and use of fungicides and growth regulators, (ii) plough tillage with reduced nitrogen fertilization (50%), (iii) cultivator tillage with standard nitrogen fertilization and use of fungicides and growth regulators, and (iv) cultivator tillage with reduced nitrogen fertilization (50%). Bulk soil (BS), as well as root-affected soil (RS), were considered for all treatments in replicates. HT-sequencing of metagenomic DNA yielded approx. 100 Giga bases (Gb) of sequence information. Taxonomic profiling of soil communities revealed the presence of 70 phyla, whereby Proteobacteria, Actinobacteria, Bacteroidetes, Planctomycetes, Acidobacteria, Thaumarchaeota, Firmicutes, Verrucomicrobia and Chloroflexi feature abundances of more than 1%. Functional microbiome profiling uncovered, i.a., numerous potential plant-beneficial, plant-growth-promoting and biocontrol traits predicted to be involved in nutrient provision, phytohormone synthesis, antagonism against pathogens and signal molecule synthesis relevant in microbe-plant interaction. Neither taxonomic nor functional microbiome profiling based on single-read analyses revealed pronounced differences regarding the farming practices applied. Soil metagenome sequences were assembled and taxonomically binned. The ten most reliable and abundant Metagenomically Assembled Genomes (MAGs) were taxonomically classified and metabolically reconstructed. Importance of the phylum Thaumarchaeota for the analyzed microbiome is corroborated by the fact that the four corresponding MAGs were predicted to oxidize ammonia (nitrification), thus contributing to the cycling of nitrogen, and in addition are most probably able to fix carbon dioxide. Moreover, Thaumarchaeota and several bacterial MAGs also possess genes with predicted functions in plant-growth-promotion. Abundances of certain MAGs (species resolution level) responded to the tillage practice, whereas the factors compartment (BS vs. RS) and nitrogen fertilization only marginally shaped MAG abundance profiles. Hence, soil management regimes promoting plant-beneficial microbiome members are very likely advantageous for the respective agrosystem, its health and carbon sequestration and accordingly may enhance plant productivity. Since Chernozem soils are highly fertile, corresponding microbiome data represent a valuable reference resource for agronomy in general.

Reconstructing the functions of endosymbiotic Mollicutes in fungus-growing ants

Mollicutes, a widespread class of bacteria associated with animals and plants, were recently identified as abundant abdominal endosymbionts in healthy workers of attine fungus-farming leaf-cutting ants. We obtained draft genomes of the two most common strains harbored by Panamanian fungus-growing ants. Reconstructions of their functional significance showed that they are independently acquired symbionts, most likely to decompose excess arginine consistent with the farmed fungal cultivars providing this nitrogen-rich amino-acid in variable quantities. Across the attine lineages, the relative abundances of the two Mollicutes strains are associated with the substrate types that foraging workers offer to fungus gardens. One of the symbionts is specific to the leaf-cutting ants and has special genomic machinery to catabolize citrate/glucose into acetate, which appears to deliver direct metabolic energy to the ant workers. Unlike other Mollicutes associated with insect hosts, both attine ant strains have complete phage-defense systems, underlining that they are actively maintained as mutualistic symbionts.