1
|
Andrews KR, Besser TE, Stalder T, Top EM, Baker KN, Fagnan MW, New DD, Schneider GM, Gal A, Andrews-Dickert R, Hunter SS, Beckmen KB, Christensen L, Justice-Allen A, Konetchy D, Lehman CP, Manlove K, Miyasaki H, Nordeen T, Roug A, Cassirer EF. Comparative genomic analysis identifies potential adaptive variation in Mycoplasma ovipneumoniae. Microb Genom 2024; 10:001279. [PMID: 39213169 PMCID: PMC11364169 DOI: 10.1099/mgen.0.001279] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2024] [Accepted: 07/16/2024] [Indexed: 09/04/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Kimberly R. Andrews
- Institute for Interdisciplinary Data Sciences, University of Idaho, Moscow, ID, USA
| | - Thomas E. Besser
- Department of Veterinary Microbiology and Pathology, Washington State University College of Veterinary Medicine, Pullman, WA, USA
| | - Thibault Stalder
- Department of Biological Sciences, University of Idaho, Moscow, ID, USA
| | - Eva M. Top
- Department of Biological Sciences, University of Idaho, Moscow, ID, USA
| | - Katherine N. Baker
- Institute for Interdisciplinary Data Sciences, University of Idaho, Moscow, ID, USA
| | - Matthew W. Fagnan
- Institute for Interdisciplinary Data Sciences, University of Idaho, Moscow, ID, USA
| | - Daniel D. New
- Institute for Interdisciplinary Data Sciences, University of Idaho, Moscow, ID, USA
| | - G. Maria Schneider
- Institute for Interdisciplinary Data Sciences, University of Idaho, Moscow, ID, USA
| | - Alexandra Gal
- Department of Biological Sciences, University of Idaho, Moscow, ID, USA
| | - Rebecca Andrews-Dickert
- Department of Physiology and Pharmacology, College of Osteopathic Medicine, Sam Houston State University, Conroe, TX, USA
| | - Samuel S. Hunter
- Institute for Interdisciplinary Data Sciences, University of Idaho, Moscow, ID, USA
| | | | - Lauren Christensen
- Department of Animal, Veterinary and Food Sciences, University of Idaho, Moscow ID, USA
| | | | - Denise Konetchy
- Department of Animal, Veterinary and Food Sciences, University of Idaho, Moscow ID, USA
| | | | - Kezia Manlove
- Department of Wildland Resources and Ecology Center, Utah State University, Logan, UT, USA
| | | | - Todd Nordeen
- Nebraska Game and Parks Commission, Alliance, NE, USA
| | - Annette Roug
- Utah Division of Wildlife Resources, Salt Lake City, UT, USA
| | | |
Collapse
|
2
|
Pascoal F, Tomasino MP, Piredda R, Quero GM, Torgo L, Poulain J, Galand PE, Fuhrman JA, Mitchell A, Tinta T, Turk Dermastia T, Fernandez-Guerra A, Vezzi A, Logares R, Malfatti F, Endo H, Dąbrowska AM, De Pascale F, Sánchez P, Henry N, Fosso B, Wilson B, Toshchakov S, Ferrant GK, Grigorov I, Vieira FRJ, Costa R, Pesant S, Magalhães C. Inter-comparison of marine microbiome sampling protocols. ISME COMMUNICATIONS 2023; 3:84. [PMID: 37598259 PMCID: PMC10439934 DOI: 10.1038/s43705-023-00278-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 06/19/2023] [Accepted: 06/23/2023] [Indexed: 08/21/2023]
Abstract
Research on marine microbial communities is growing, but studies are hard to compare because of variation in seawater sampling protocols. To help researchers in the inter-comparison of studies that use different seawater sampling methodologies, as well as to help them design future sampling campaigns, we developed the EuroMarine Open Science Exploration initiative (EMOSE). Within the EMOSE framework, we sampled thousands of liters of seawater from a single station in the NW Mediterranean Sea (Service d'Observation du Laboratoire Arago [SOLA], Banyuls-sur-Mer), during one single day. The resulting dataset includes multiple seawater processing approaches, encompassing different material-type kinds of filters (cartridge membrane and flat membrane), three different size fractionations (>0.22 µm, 0.22-3 µm, 3-20 µm and >20 µm), and a number of different seawater volumes ranging from 1 L up to 1000 L. We show that the volume of seawater that is filtered does not have a significant effect on prokaryotic and protist diversity, independently of the sequencing strategy. However, there was a clear difference in alpha and beta diversity between size fractions and between these and "whole water" (with no pre-fractionation). Overall, we recommend care when merging data from datasets that use filters of different pore size, but we consider that the type of filter and volume should not act as confounding variables for the tested sequencing strategies. To the best of our knowledge, this is the first time a publicly available dataset effectively allows for the clarification of the impact of marine microbiome methodological options across a wide range of protocols, including large-scale variations in sampled volume.
Collapse
Affiliation(s)
- Francisco Pascoal
- Interdisciplinary Centre of Marine and Environmental Research, University of Porto, Terminal de Cruzeiros do Porto de Leixões, Av. General Norton de Matos s/n, 4450-208, Porto, Portugal
- Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, rua do Campo Alegre s/n, 4169- 007, Porto, Portugal
| | - Maria Paola Tomasino
- Interdisciplinary Centre of Marine and Environmental Research, University of Porto, Terminal de Cruzeiros do Porto de Leixões, Av. General Norton de Matos s/n, 4450-208, Porto, Portugal
| | - Roberta Piredda
- Integrative Marine Ecology Department, Stazione Zoologica Anton Dohrn, Naples, Italy
| | - Grazia Marina Quero
- Institute for Biological Resources and Marine Biotechnologies, National Research Council (IRBIM-CNR), Largo Fiera della Pesca 2, 60125, Ancona, Italy
| | - Luís Torgo
- Faculty of Computer Science, Dalhousie University, Halifax, NS, Canada
| | - Julie Poulain
- Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, 2 Rue Gaston Crémieux, 91057, Evry, France
| | - Pierre E Galand
- Sorbonne Université, CNRS, Laboratoire d'Écogéochimie des Environnements Benthiques (LECOB), Observatoire Océanologique de Banyuls, Banyuls-sur-Mer, France
| | - Jed A Fuhrman
- Marine & Environmental Biology, Department of Biological Sciences, University of Southern California (USC), Los Angeles, CA, USA
| | - Alex Mitchell
- EMBL's European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridgeshire, CB10 1SD, UK
| | - Tinkara Tinta
- National Institute of Biology, Marine Biology Station Piran, Piran, Slovenia
| | | | - Antonio Fernandez-Guerra
- Lundbeck Foundation GeoGenetics Centre, GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
| | - Alessandro Vezzi
- Department of Biology, University of Padua, Via U. Bassi 58/B, 35131, Padua, Italy
| | - Ramiro Logares
- Institute of Marine Sciences (ICM), CSIC. Passeig Marítim de la Barceloneta, 37-49, ES08003, Barcelona, Spain
| | | | - Hisashi Endo
- Bioinformatics Center, Institute for Chemical Research, Kyoto University, Gokasho, Uji, Japan
| | - Anna Maria Dąbrowska
- Department of Marine Ecology, Institute of Oceanology Polish Academy of Sciences, Sopot, Poland
| | - Fabio De Pascale
- Department of Biology, University of Padua, Via U. Bassi 58/B, 35131, Padua, Italy
| | - Pablo Sánchez
- Institute of Marine Sciences (ICM), CSIC. Passeig Marítim de la Barceloneta, 37-49, ES08003, Barcelona, Spain
| | - Nicolas Henry
- Sorbonne Université, CNRS, Station Biologique de Roscoff, AD2M ECOMAP, UMR 7144, Roscoff, France
- CNRS, FR2424, ABiMS, Station Biologique de Roscoff, Sorbonne Université, Roscoff, France
| | - Bruno Fosso
- Department of Biosciences, Biotechnologies and Environment, University of Bari, 70126, Bari, Italy
| | - Bryan Wilson
- Department of Biology, John Krebs Field Station, University of Oxford, Wytham, OX2 8QJ, UK
| | | | | | - Ivo Grigorov
- Technical University of Denmark, National Institute of Aquatic Resources, Kgs. Lyngby, Denmark
| | | | - Rodrigo Costa
- Department of Bioengineering, Instituto Superior Técnico, University of Lisbon, Av. Rovisco Pais, 1049-001, Lisbon, Portugal
- Institute for Bioengineering and Biosciences (iBB) and i4HB-Institute for Health and Bioeconomy, Instituto Superior Técnico, University of Lisbon, Av. Rovisco Pais, 1049-001, Lisbon, Portugal
| | - Stéphane Pesant
- EMBL's European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridgeshire, CB10 1SD, UK.
| | - Catarina Magalhães
- Interdisciplinary Centre of Marine and Environmental Research, University of Porto, Terminal de Cruzeiros do Porto de Leixões, Av. General Norton de Matos s/n, 4450-208, Porto, Portugal.
- Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, rua do Campo Alegre s/n, 4169- 007, Porto, Portugal.
| |
Collapse
|
3
|
Liao J, Shenhav L, Urban JA, Serrano M, Zhu B, Buck GA, Korem T. Microdiversity of the vaginal microbiome is associated with preterm birth. Nat Commun 2023; 14:4997. [PMID: 37591872 PMCID: PMC10435516 DOI: 10.1038/s41467-023-40719-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Accepted: 08/07/2023] [Indexed: 08/19/2023] Open
Abstract
Preterm birth (PTB) is the leading cause of neonatal morbidity and mortality. The vaginal microbiome has been associated with PTB, yet the mechanisms underlying this association are not fully understood. Understanding microbial genetic adaptations to selective pressures, especially those related to the host, may yield insights into these associations. Here, we analyze metagenomic data from 705 vaginal samples collected during pregnancy from 40 women who delivered preterm spontaneously and 135 term controls from the Multi-Omic Microbiome Study-Pregnancy Initiative. We find that the vaginal microbiome of pregnancies that ended preterm exhibited unique genetic profiles. It was more genetically diverse at the species level, a result which we validate in an additional cohort, and harbored a higher richness and diversity of antimicrobial resistance genes, likely promoted by transduction. Interestingly, we find that Gardnerella species drove this higher genetic diversity, particularly during the first half of the pregnancy. We further present evidence that Gardnerella spp. underwent more frequent recombination and stronger purifying selection in genes involved in lipid metabolism. Overall, our population genetics analyses reveal associations between the vaginal microbiome and PTB and suggest that evolutionary processes acting on vaginal microbes may play a role in adverse pregnancy outcomes such as PTB.
Collapse
Affiliation(s)
- Jingqiu Liao
- Program for Mathematical Genomics, Department of Systems Biology, Columbia University Irving Medical Center, New York, NY, USA.
- Department of Civil and Environmental Engineering, Virginia Tech, Blacksburg, VA, USA.
| | - Liat Shenhav
- Center for Studies in Physics and Biology, Rockefeller University, New York, NY, USA
| | - Julia A Urban
- Program for Mathematical Genomics, Department of Systems Biology, Columbia University Irving Medical Center, New York, NY, USA
| | - Myrna Serrano
- Department of Microbiology and Immunology, School of Medicine, Virginia Commonwealth University, Richmond, VA, USA
- Center for Microbiome Engineering and Data Analysis, Virginia Commonwealth University, Richmond, VA, USA
| | - Bin Zhu
- Department of Microbiology and Immunology, School of Medicine, Virginia Commonwealth University, Richmond, VA, USA
- Center for Microbiome Engineering and Data Analysis, Virginia Commonwealth University, Richmond, VA, USA
| | - Gregory A Buck
- Department of Microbiology and Immunology, School of Medicine, Virginia Commonwealth University, Richmond, VA, USA
- Center for Microbiome Engineering and Data Analysis, Virginia Commonwealth University, Richmond, VA, USA
- Department of Computer Science, School of Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Tal Korem
- Program for Mathematical Genomics, Department of Systems Biology, Columbia University Irving Medical Center, New York, NY, USA.
- Department of Obstetrics and Gynecology, Columbia University Irving Medical Center, New York, NY, USA.
- CIFAR Azrieli Global Scholars program, CIFAR, Toronto, ON, Canada.
| |
Collapse
|
4
|
Liao J, Shenhav L, Urban JA, Serrano M, Zhu B, Buck GA, Korem T. Microdiversity of the Vaginal Microbiome is Associated with Preterm Birth. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.13.523991. [PMID: 36711990 PMCID: PMC9882146 DOI: 10.1101/2023.01.13.523991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Preterm birth (PTB) is the leading cause of neonatal morbidity and mortality. The vaginal microbiome has been associated with PTB, yet the mechanisms underlying this association are not fully understood. Understanding microbial genetic adaptations to selective pressures, especially those related to the host, may yield new insights into these associations. To this end, we analyzed metagenomic data from 705 vaginal samples collected longitudinally during pregnancy from 40 women who delivered preterm spontaneously and 135 term controls from the Multi-Omic Microbiome Study-Pregnancy Initiative (MOMS-PI). We find that the vaginal microbiome of pregnancies that ended preterm exhibits unique genetic profiles. It is more genetically diverse at the species level, a result which we validate in an additional cohort, and harbors a higher richness and diversity of antimicrobial resistance genes, likely promoted by transduction. Interestingly, we find that Gardnerella species, a group of central vaginal pathobionts, are driving this higher genetic diversity, particularly during the first half of the pregnancy. We further present evidence that Gardnerella spp. undergoes more frequent recombination and stronger purifying selection in genes involved in lipid metabolism. Overall, our results reveal novel associations between the vaginal microbiome and PTB using population genetics analyses, and suggest that evolutionary processes acting on the vaginal microbiome may play a vital role in adverse pregnancy outcomes such as preterm birth.
Collapse
Affiliation(s)
- Jingqiu Liao
- Program for Mathematical Genomics, Department of Systems Biology, Columbia University Irving Medical Center, New York, NY, USA
- Department of Civil and Environmental Engineering, Virginia Tech, Blacksburg, VA, USA
| | - Liat Shenhav
- Center for Studies in Physics and Biology, Rockefeller University, New York, NY, USA
| | - Julia A. Urban
- Program for Mathematical Genomics, Department of Systems Biology, Columbia University Irving Medical Center, New York, NY, USA
| | - Myrna Serrano
- Department of Microbiology and Immunology, School of Medicine, Virginia Commonwealth University, Richmond, VA, USA
- Center for Microbiome Engineering and Data Analysis, Virginia Commonwealth University, Richmond, VA, USA
| | - Bin Zhu
- Department of Microbiology and Immunology, School of Medicine, Virginia Commonwealth University, Richmond, VA, USA
- Center for Microbiome Engineering and Data Analysis, Virginia Commonwealth University, Richmond, VA, USA
| | - Gregory A. Buck
- Department of Microbiology and Immunology, School of Medicine, Virginia Commonwealth University, Richmond, VA, USA
- Center for Microbiome Engineering and Data Analysis, Virginia Commonwealth University, Richmond, VA, USA
- Department of Computer Science, School of Engineering, Virginia Commonwealth University, Richmond, VA, USA
| | - Tal Korem
- Program for Mathematical Genomics, Department of Systems Biology, Columbia University Irving Medical Center, New York, NY, USA
- Department of Obstetrics and Gynecology, Columbia University Irving Medical Center, New York, NY, USA
- CIFAR Azrieli Global Scholars program, CIFAR, Toronto, Canada
| |
Collapse
|
5
|
Roda-Garcia JJ, Haro-Moreno JM, López-Pérez M. Evolutionary pathways for deep-sea adaptation in marine planktonic Actinobacteriota. Front Microbiol 2023; 14:1159270. [PMID: 37234526 PMCID: PMC10205998 DOI: 10.3389/fmicb.2023.1159270] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2023] [Accepted: 04/24/2023] [Indexed: 05/28/2023] Open
Abstract
The deep ocean, one of the largest ecosystems on earth, is dominated by microorganisms that are keystones in the regulation of biogeochemical cycles. However, the evolutionary pathways underlying the specific adaptations required (e.g., high pressure and low temperature) by this unique niche remain understudied. Here, we analyzed the first representatives belonging to the order Acidimicrobiales, a group of marine planktonic Actinobacteriota, that specifically inhabits the aphotic zone of the oceanic water column (>200 m). Compared with their epipelagic counterparts, deep-sea representatives showed the same evolution in genome architecture with higher GC content, longer intergenic spaces as well as higher nitrogen (N-ARSC) and lower carbon (C-ARSC) content in encoded amino acid residue side chains consistent with the higher nitrogen concentration and lower carbon concentration in deep waters compared to the photic zone. Metagenomic recruitment showed distribution patterns that allowed the description of different ecogenomic units within the three deep water-associated genera defined by our phylogenomic analyses (UBA3125, S20-B6 and UBA9410). The entire genus UBA3125 was found exclusively associated with oxygen minimum zones linked to the acquisition of genes involved in denitrification. Genomospecies of genus S20-B6 recruited in samples from both mesopelagic (200-1,000 m) and bathypelagic (1000-4,000 m) zones, including polar regions. Diversity in the genus UBA9410 was higher, with genomospecies widely distributed in temperate zones, others in polar regions, and the only genomospecies associated with abyssal zones (>4,000 m). At the functional level, groups beyond the epipelagic zone have a more complex transcriptional regulation including in their genomes a unique WhiB paralog. In addition, they showed higher metabolic potential for organic carbon and carbohydrate degradation as well as the ability to accumulate glycogen as a source of carbon and energy. This could compensate for energy metabolism in the absence of rhodopsins, which is only present in genomes associated with the photic zone. The abundance in deep samples of cytochrome P450 monooxygenases associated with the genomes of this order suggests an important role in remineralization of recalcitrant compounds throughout the water column.
Collapse
|
6
|
Roda-Garcia JJ, Haro-Moreno JM, Rodriguez-Valera F, Almagro-Moreno S, López-Pérez M. Single-amplified genomes reveal most streamlined free-living marine bacteria. Environ Microbiol 2023. [PMID: 36755376 DOI: 10.1111/1462-2920.16348] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 02/04/2023] [Indexed: 02/10/2023]
Abstract
Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.
Collapse
Affiliation(s)
- Juan J Roda-Garcia
- Evolutionary Genomics Group, Departamento Producción Vegetal y Microbiología, Universidad Miguel Hernández, Alicante, Spain
| | - Jose M Haro-Moreno
- Evolutionary Genomics Group, Departamento Producción Vegetal y Microbiología, Universidad Miguel Hernández, Alicante, Spain
| | - Francisco Rodriguez-Valera
- Evolutionary Genomics Group, Departamento Producción Vegetal y Microbiología, Universidad Miguel Hernández, Alicante, Spain
| | - Salvador Almagro-Moreno
- Burnett School of Biomedical Sciences, University of Central Florida, Orlando, Florida, USA.,National Center for Integrated Coastal Research, University of Central Florida, Orlando, Florida, USA
| | - Mario López-Pérez
- Evolutionary Genomics Group, Departamento Producción Vegetal y Microbiología, Universidad Miguel Hernández, Alicante, Spain
| |
Collapse
|
7
|
Cyanobacterial Blooms Are Not a Result of Positive Selection by Freshwater Eutrophication. Microbiol Spectr 2022; 10:e0319422. [PMID: 36445094 PMCID: PMC9769789 DOI: 10.1128/spectrum.03194-22] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
Long-standing cyanobacterial harmful algal blooms (CyanoHABs) are known to result from synergistic interaction between elevated nutrients and superior ecophysiology of cyanobacteria. However, it remains to be determined whether CyanoHABs are a result of positive selection by eutrophic waters. To address this, we conducted molecular evolutionary analyses on the genomes of 9 bloom-forming cyanobacteria, combined with pangenomics and metatranscriptomics. The results showed no positive selection by water eutrophication. Instead, all homologous genes in the species are under strong purifying selection based on the ratio of divergence at nonsynonymous and synonymous sites (dN/dS) and phylogeny. The dN/dS < 0.85 (median = 0.3) for all homologous genes are similar between the genes in the pathways driving CyanoHABs and housekeeping functions. Phylogenetic support for non-positive selection comes from the mixed clustering of strains: strains of the same species from diverse geographic origins form the same clusters, while strains from the same origins form different clusters. Further support lies in the codon adaptation index (CAI) and single nucleotide polymorphism (SNP). The CAI ranged from 0.42 to 0.9 (mean = 0.75), which indicates high-level codon usage bias; the pathways for CyanoHABs and housekeeping functions showed a similar CAI. Interestingly, CAI was negatively correlated with gene expression in 3 metatranscriptomes. The numbers of SNPs were concentrated around 5 to 50. As the SNP number increases, the gene expression level decreases. These negative correlations agree with the population-level dN/dS and phylogeny in supporting purifying selection in bloom-forming cyanobacteria. In summary, superior ecophysiology appears to be acquired prior to water eutrophication. IMPORTANCE CyanoHABs are global environmental hazards, and their mechanisms of action are being intensively investigated. On an ecological scale, CyanoHABs are consequences of synergistic interactions between biological functions and elevated nutrients in eutrophic waters. On an evolutionary scale, one important question is how bloom-forming cyanobacteria acquire these superior biological functions. There are several possibilities, including adaptive evolution and horizontal gene transfer. Here, we explored the possibility of positive selection. We reasoned that there are two possible periods for cyanobacteria to acquire these functions: before the onset of water eutrophication or during water eutrophication. Either way, there should be molecular signatures in protein sequences for positive selection. Interestingly, we found no positive selection by water eutrophication, but strong purifying selection instead on nearly all the genes, suggesting these superior functions aiding CyanoHABs are acquired prior to water eutrophication.
Collapse
|
8
|
Okazaki Y, Nakano SI, Toyoda A, Tamaki H. Long-Read-Resolved, Ecosystem-Wide Exploration of Nucleotide and Structural Microdiversity of Lake Bacterioplankton Genomes. mSystems 2022; 7:e0043322. [PMID: 35938717 PMCID: PMC9426551 DOI: 10.1128/msystems.00433-22] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Accepted: 07/06/2022] [Indexed: 12/24/2022] Open
Abstract
Reconstruction of metagenome-assembled genomes (MAGs) has become a fundamental approach in microbial ecology. However, a MAG is hardly complete and overlooks genomic microdiversity because metagenomic assembly fails to resolve microvariants among closely related genotypes. Aiming at understanding the universal factors that drive or constrain prokaryotic genome diversification, we performed an ecosystem-wide high-resolution metagenomic exploration of microdiversity by combining spatiotemporal (2 depths × 12 months) sampling from a pelagic freshwater system, high-quality MAG reconstruction using long- and short-read metagenomic sequences, and profiling of single nucleotide variants (SNVs) and structural variants (SVs) through mapping of short and long reads to the MAGs, respectively. We reconstructed 575 MAGs, including 29 circular assemblies, providing high-quality reference genomes of freshwater bacterioplankton. Read mapping against these MAGs identified 100 to 101,781 SNVs/Mb and 0 to 305 insertions, 0 to 467 deletions, 0 to 41 duplications, and 0 to 6 inversions for each MAG. Nonsynonymous SNVs were accumulated in genes potentially involved in cell surface structural modification to evade phage recognition. Most (80.2%) deletions overlapped with a gene coding region, and genes of prokaryotic defense systems were most frequently (>8% of the genes) overlapped with a deletion. Some such deletions exhibited a monthly shift in their allele frequency, suggesting a rapid turnover of genotypes in response to phage predation. MAGs with extremely low microdiversity were either rare or opportunistic bloomers, suggesting that population persistency is key to their genomic diversification. The results concluded that prokaryotic genomic diversification is driven primarily by viral load and constrained by a population bottleneck. IMPORTANCE Identifying intraspecies genomic diversity (microdiversity) is crucial to understanding microbial ecology and evolution. However, microdiversity among environmental assemblages is not well investigated, because most microbes are difficult to culture. In this study, we performed cultivation-independent exploration of bacterial genomic microdiversity in a lake ecosystem using a combination of short- and long-read metagenomic analyses. The results revealed the broad spectrum of genomic microdiversity among the diverse bacterial species in the ecosystem, which has been overlooked by conventional approaches. Our ecosystem-wide exploration further allowed comparative analysis among the genomes and genes and revealed factors behind microbial genomic diversification, namely, that diversification is driven primarily by resistance against viral infection and constrained by the population size.
Collapse
Affiliation(s)
- Yusuke Okazaki
- Institute for Chemical Research, Kyoto University, Uji, Kyoto, Japan
- Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology, Tsukuba, Ibaraki, Japan
| | - Shin-ichi Nakano
- Center for Ecological Research, Kyoto University, Otsu, Shiga, Japan
| | - Atsushi Toyoda
- Advanced Genomics Center, National Institute of Genetics, Mishima City, Shizuoka, Japan
| | - Hideyuki Tamaki
- Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology, Tsukuba, Ibaraki, Japan
| |
Collapse
|
9
|
Gao Y, Wu M. Microbial genomic trait evolution is dominated by frequent and rare pulsed evolution. SCIENCE ADVANCES 2022; 8:eabn1916. [PMID: 35857501 PMCID: PMC9286504 DOI: 10.1126/sciadv.abn1916] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Accepted: 06/02/2022] [Indexed: 06/12/2023]
Abstract
On the macroevolutionary time scale, does trait evolution proceed gradually or by rapid bursts (pulses) separated by prolonged periods of stasis or slow evolution? Although studies have shown that pulsed evolution is prevalent in animals, our knowledge about the tempo and mode of evolution across the tree of life is very limited. This long-standing debate calls for a test in bacteria and archaea, the most ancient and diverse forms of life with unique population genetic properties. Using a likelihood-based framework, we show that pulsed evolution is not only present but also prevalent and predominant in microbial genomic trait evolution. We detected two distinct types of pulsed evolution (small frequent and large rare jumps) that are predicted by the punctuated equilibrium and quantum evolution theories. Our findings suggest that major bacterial lineages could have originated in quick bursts and that pulsed evolution is a common theme across the tree of life.
Collapse
|
10
|
Pascoal F, Costa R, Assmy P, Duarte P, Magalhães C. Exploration of the Types of Rarity in the Arctic Ocean from the Perspective of Multiple Methodologies. MICROBIAL ECOLOGY 2022; 84:59-72. [PMID: 34405249 DOI: 10.1007/s00248-021-01821-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2021] [Accepted: 07/12/2021] [Indexed: 06/13/2023]
Abstract
The Arctic Ocean is facing rapid environmental changes with cascading effects on the entire Arctic marine ecosystem. However, we have a limited understanding of the consequences such changes have on bacteria and archaea (prokaryotes) at the base of the marine food web. In this study, we show how the prokaryotic rare biosphere behaves over a range of highly heterogeneous environmental conditions using 16S rRNA gene reads from amplicon and metagenome sequencing data from seawater samples collected during the Norwegian young sea ICE expedition between late winter and early summer. The prokaryotic rare biosphere was analyzed using different approaches: amplicon sequence variants and operational taxonomic units from the 16S rRNA gene amplicons and operational taxonomic units from the 16S rRNA genes of the metagenomes. We found that prokaryotic rare biosphere communities are specific to certain water masses, and that the majority of the rare taxa identified were always rare and disappeared in at least one sample under changing conditions, suggesting their high sensitivity to environmental heterogeneity. In addition, our methodological comparison revealed a good performance of 16S rRNA gene amplicon sequencing in describing rare biosphere patterns, while the metagenome-derived data were better to capture a significant diversity of so-far uncultivated rare taxa. Our analysis on the dynamics of the rare prokaryotic biosphere, by combining different methodological approaches, improves the description of the types of rarity predicted from Community Assembly theory in the Arctic Ocean.
Collapse
Affiliation(s)
- Francisco Pascoal
- Interdisciplinary Centre of Marine and Environmental Research (CIIMAR), University of Porto, Porto, Portugal.
- Faculty of Sciences, University of Porto, Porto, Portugal.
| | - Rodrigo Costa
- Department of Bioengineering, Institute for Bioengineering and Biosciences (iBB), Instituto Superior Técnico, University of Lisbon, Lisbon, Portugal
- Centre of Marine Sciences (CCMAR), University of Algarve, Faro, Portugal
- U.S. Department of Energy Joint Genome Institute, Berkeley, CA, USA
- Environmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Philipp Assmy
- Norwegian Polar Institute, Fram Centre, N-9296, Tromsø, Norway
| | - Pedro Duarte
- Norwegian Polar Institute, Fram Centre, N-9296, Tromsø, Norway
| | - Catarina Magalhães
- Interdisciplinary Centre of Marine and Environmental Research (CIIMAR), University of Porto, Porto, Portugal
- Faculty of Sciences, University of Porto, Porto, Portugal
- School of Science & Engineering, University of Waikato, Hamilton, New Zealand
- Ocean Frontier Institute, Dalhousie University, Halifax, NS, Canada
| |
Collapse
|
11
|
Martinez-Gutierrez CA, Aylward FO. Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales. PLoS Genet 2022; 18:e1010220. [PMID: 35605022 PMCID: PMC9166353 DOI: 10.1371/journal.pgen.1010220] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 06/03/2022] [Accepted: 04/26/2022] [Indexed: 12/30/2022] Open
Abstract
The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel's lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.
Collapse
Affiliation(s)
| | - Frank O. Aylward
- Department of Biological Sciences, Virginia Tech, Blacksburg, Virginia, United States of America
- Center for Emerging, Zoonotic, and Arthropod-borne Pathogens, Virginia Tech, Blacksburg, Virginia, United States of America
| |
Collapse
|
12
|
Gao Y, Wu M. Modeling Pulsed Evolution and Time-independent Variation Improves the Confidence Level of Ancestral and Hidden State Predictions. Syst Biol 2022; 71:1225-1232. [PMID: 35212761 DOI: 10.1093/sysbio/syac016] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Revised: 02/14/2022] [Accepted: 02/22/2022] [Indexed: 11/14/2022] Open
Abstract
Ancestral state reconstruction is not only a fundamental tool for studying trait evolution, but also very useful for predicting the unknown trait values (hidden states) of extant species. A well-known problem in ancestral and hidden state predictions is that the uncertainty associated with predictions can be so large that predictions themselves are of little use. Therefore, for meaningful interpretation of predicted traits and hypothesis testing, it is prudent to accurately assess the uncertainty of the predictions. Commonly used constant-rate Brownian motion (BM) model fails to capture the complexity of tempo and mode of trait evolution in nature, making predictions under the BM model vulnerable to lack-of-fit errors from model misspecification. Using empirical data (mammalian body size and bacterial genome size), we show that the distribution of residual Z-scores under the BM model is neither homoscedastic nor normal as expected. Consequently, the 95% confidence intervals (CIs) of predicted traits are so unreliable that the actual coverage probability ranges from 33% (strongly permissive) to 100% (strongly conservative). Alternative methods such as BayesTraits and StableTraits that allow variable rates in evolution improve the predictions but are computationally expensive. Here we develop RasperGade, a method of ancestral and hidden state prediction that uses the Levy process to explicitly model gradual evolution, pulsed evolution and time-independent variation. Using the same empirical data, we show that RasperGade outperforms both BayesTraits and StableTraits in providing reliable confidence estimates and is orders-of-magnitude faster. Our results suggest that, when predicting the ancestral and hidden states of continuous traits, the rate variation should always be assessed and the quality of confidence estimates should always be examined.
Collapse
Affiliation(s)
- Yingnan Gao
- Department of Biology, University of Virginia, Charlottesville, VA 22904, USA
| | - Martin Wu
- Department of Biology, University of Virginia, Charlottesville, VA 22904, USA
| |
Collapse
|
13
|
Baker LJ, Reich HG, Kitchen SA, Grace Klinges J, Koch HR, Baums IB, Muller EM, Thurber RV. The coral symbiont Candidatus Aquarickettsia is variably abundant in threatened Caribbean acroporids and transmitted horizontally. THE ISME JOURNAL 2022; 16:400-411. [PMID: 34363004 PMCID: PMC8776821 DOI: 10.1038/s41396-021-01077-8] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Revised: 06/28/2021] [Accepted: 07/22/2021] [Indexed: 02/07/2023]
Abstract
The symbiont "Candidatus Aquarickettsia rohweri" infects a diversity of aquatic hosts. In the threatened Caribbean coral, Acropora cervicornis, Aquarickettsia proliferates in response to increased nutrient exposure, resulting in suppressed growth and increased disease susceptibility and mortality of coral. This study evaluated the extent, as well as the ecology and evolution of Aquarickettsia infecting threatened corals, Ac. cervicornis, and Ac. palmata and their hybrid ("Ac. prolifera"). Aquarickettsia was found in all acroporids, with coral host and geographic location impacting the infection magnitude. Phylogenomic and genome-wide single-nucleotide variant analysis of Aquarickettsia found phylogenetic clustering by geographic region, not by coral taxon. Analysis of Aquarickettsia fixation indices suggests multiple sequential infections of the same coral colony are unlikely. Furthermore, relative to other Rickettsiales species, Aquarickettsia is undergoing positive selection, with Florida populations experiencing greater positive selection relative to other Caribbean locations. This may be due in part to Aquarickettsia proliferating in response to greater nutrient stress in Florida, as indicated by greater in situ replication rates in these corals. Aquarickettsia was not found to significantly codiversify with either the coral animal or the coral's algal symbiont (Symbiodinium "fitti"). Quantitative PCR analysis showed that gametes, larvae, recruits, and juveniles from susceptible, captive-reared coral genets were not infected with Aquarickettsia. Thus, horizontal transmission of Aquarickettsia via coral mucocytes or an unidentified host is more likely. The prevalence of Aquarickettsia in Ac. cervicornis and its high abundance in the Florida coral population suggests that coral disease mitigation efforts focus on preventing early infection via horizontal transmission.
Collapse
Affiliation(s)
- Lydia J Baker
- Department of Microbiology, Oregon State University, Corvallis, OR, USA.
| | - Hannah G Reich
- Department of Biology, The Pennsylvania State University, University Park, PA, USA
| | - Sheila A Kitchen
- Division of Biology and Biological Engineering, California Institute of Science and Technology, Pasadena, CA, USA
| | - J Grace Klinges
- Department of Microbiology, Oregon State University, Corvallis, OR, USA
| | - Hanna R Koch
- Coral Restoration Program, Mote Marine Laboratory, Summerland Key, FL, USA
| | - Iliana B Baums
- Department of Biology, The Pennsylvania State University, University Park, PA, USA
| | - Erinn M Muller
- Coral Restoration Program, Mote Marine Laboratory, Summerland Key, FL, USA
| | | |
Collapse
|
14
|
Environmental stress leads to genome streamlining in a widely distributed species of soil bacteria. THE ISME JOURNAL 2022; 16:423-434. [PMID: 34408268 PMCID: PMC8776746 DOI: 10.1038/s41396-021-01082-x] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Revised: 07/14/2021] [Accepted: 07/28/2021] [Indexed: 02/07/2023]
Abstract
Bacteria have highly flexible pangenomes, which are thought to facilitate evolutionary responses to environmental change, but the impacts of environmental stress on pangenome evolution remain unclear. Using a landscape pangenomics approach, I demonstrate that environmental stress leads to consistent, continuous reduction in genome content along four environmental stress gradients (acidity, aridity, heat, salinity) in naturally occurring populations of Bradyrhizobium diazoefficiens (widespread soil-dwelling plant mutualists). Using gene-level network and duplication functional traits to predict accessory gene distributions across environments, genes predicted to be superfluous are more likely lost in high stress, while genes with multi-functional roles are more likely retained. Genes with higher probabilities of being lost with stress contain significantly higher proportions of codons under strong purifying and positive selection. Gene loss is widespread across the entire genome, with high gene-retention hotspots in close spatial proximity to core genes, suggesting Bradyrhizobium has evolved to cluster essential-function genes (accessory genes with multifunctional roles and core genes) in discrete genomic regions, which may stabilise viability during genomic decay. In conclusion, pangenome evolution through genome streamlining are important evolutionary responses to environmental change. This raises questions about impacts of genome streamlining on the adaptive capacity of bacterial populations facing rapid environmental change.
Collapse
|
15
|
Wei SA, Xu R, Ji YY, Ding ZW, Zou YZ. Deduction and exploration of the evolution and function of vertebrate GFPT family. Genes Genomics 2022; 44:175-185. [PMID: 35038160 DOI: 10.1007/s13258-021-01188-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Accepted: 11/06/2021] [Indexed: 11/27/2022]
Abstract
BACKGROUND Glutamine-fructose-6-phosphate aminotransferase (GFPT) is a key factor in the hexosamine metabolism pathway. It regulates the downstream factor O-GlcNAc to change cell function and plays an important role in the metabolism and immune process of tissues and organs. However, the evolutionary relationship of GFPT family proteins in vertebrates has not been elucidated. OBJECTIVE To deduce and explore the evolution and function of vertebrate GFPT family. METHODS 18 GFPT sequences were obtained from Homo sapiens (H. sapiens), Trachypithecus francoisi (T. francoisi), Mus musculus (M. musculus), Rattus norvegicus (R. norvegicus), Gallus gallus (G. gallus), Zootoca vivipara (Z. vivipara), Xenopus tropicalis (X. tropicalis), Danio rerio (D. rerio), Rhincodon typus (R. typus), Plasmodium relictum from National Center for Biotechnology Information (NCBI). The physical and chemical characteristics and molecular evolution of GFPT family proteins and nucleic acid sequences were analyzed by ClustalX2, Gene Doc, MEGA-X, SMART, Datamonkey, R etc. RESULTS: Based on the neighbor-joining (NJ) phylogenetic tree and evolution fingerprints, GFPT family members of vertebrates can be divided into two groups: the GFPT1 group and the GFPT2 group. Seven positive selection sites were identified by IFEL and integrated methods mixed effects model of evolution (MEME) and fixed effects likelihood (REL). Finally, we predicted 28 phosphorylation sites and 18 ubiquitousness sites in the human GFPT1 sequence, 10 phosphorylation sites, and five ubiquitousness sites in GFPT2. Gene ontology (GO) analyzes the protein molecules and KEGG signaling pathways of vertebrates interacting with GFPT family proteins. CONCLUSIONS Our work confirmed that higher animals GFPT family may have differentiated GFPT1 and GFPT2, which meets their own functional needs. This knowledge answers the question what the origin and evolution of GFPT family in vertebrates and provided the basis for disease treatment and function research of GFPT protein.
Collapse
Affiliation(s)
- Si-Ang Wei
- Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital, Fudan University, 180 Feng Lin Road, Shanghai, 200032, People's Republic of China
| | - Ran Xu
- Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital, Fudan University, 180 Feng Lin Road, Shanghai, 200032, People's Republic of China
| | - Yu-Yao Ji
- Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital, Fudan University, 180 Feng Lin Road, Shanghai, 200032, People's Republic of China
| | - Zhi-Wen Ding
- Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital, Fudan University, 180 Feng Lin Road, Shanghai, 200032, People's Republic of China.
| | - Yun-Zeng Zou
- Shanghai Institute of Cardiovascular Diseases, Zhongshan Hospital, Fudan University, 180 Feng Lin Road, Shanghai, 200032, People's Republic of China.
| |
Collapse
|
16
|
Shoemaker WR, Chen D, Garud NR. Comparative Population Genetics in the Human Gut Microbiome. Genome Biol Evol 2022; 14:evab116. [PMID: 34028530 PMCID: PMC8743038 DOI: 10.1093/gbe/evab116] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/22/2021] [Indexed: 11/13/2022] Open
Abstract
Genetic variation in the human gut microbiome is responsible for conferring a number of crucial phenotypes like the ability to digest food and metabolize drugs. Yet, our understanding of how this variation arises and is maintained remains relatively poor. Thus, the microbiome remains a largely untapped resource, as the large number of coexisting species in the microbiome presents a unique opportunity to compare and contrast evolutionary processes across species to identify universal trends and deviations. Here we outline features of the human gut microbiome that, while not unique in isolation, as an assemblage make it a system with unparalleled potential for comparative population genomics studies. We consciously take a broad view of comparative population genetics, emphasizing how sampling a large number of species allows researchers to identify universal evolutionary dynamics in addition to new genes, which can then be leveraged to identify exceptional species that deviate from general patterns. To highlight the potential power of comparative population genetics in the microbiome, we reanalyze patterns of purifying selection across ∼40 prevalent species in the human gut microbiome to identify intriguing trends which highlight functional categories in the microbiome that may be under more or less constraint.
Collapse
Affiliation(s)
- William R Shoemaker
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, California, USA
| | - Daisy Chen
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, California, USA
| | - Nandita R Garud
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, California, USA
- Department of Human Genetics, University of California, Los Angeles, California, USA
| |
Collapse
|
17
|
Phylogenomics of SAR116 Clade Reveals Two Subclades with Different Evolutionary Trajectories and an Important Role in the Ocean Sulfur Cycle. mSystems 2021; 6:e0094421. [PMID: 34609172 PMCID: PMC8547437 DOI: 10.1128/msystems.00944-21] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The SAR116 clade within the class Alphaproteobacteria represents one of the most abundant groups of heterotrophic bacteria inhabiting the surface of the ocean. The small number of cultured representatives of SAR116 (only two to date) is a major bottleneck that has prevented an in-depth study at the genomic level to understand the relationship between genome diversity and its role in the marine environment. In this study, we use all publicly available genomes to provide a genomic overview of the phylogeny, metabolism, and biogeography within the SAR116 clade. This increased genomic diversity has led to the discovery of two subclades that, despite coexisting in the same environment, display different properties in their genomic makeup. One represents a novel subclade for which no pure cultures have been isolated and is composed mainly of single-amplified genomes (SAGs). Genomes within this subclade showed convergent evolutionary trajectories with more streamlined features, such as low GC content (ca. 30%), short intergenic spacers (<22 bp), and strong purifying selection (low ratio of nonsynonymous to synonymous polymorphisms [dN/dS]). Besides, they were more abundant in metagenomic databases recruiting at the deep chlorophyll maximum. Less abundant and restricted to the upper photic layers of the global ocean, the other subclade of SAR116, enriched in metagenome-assembled genomes (MAGs), included the only two pure cultures. Genomic analysis suggested that both clades have a significant role in the sulfur cycle with differences in the way both clades can metabolize dimethylsulfoniopropionate (DMSP). IMPORTANCE The SAR116 clade of Alphaproteobacteria is a ubiquitous group of heterotrophic bacteria inhabiting the surface of the ocean, but the information about their ecology and population genomic diversity is scarce due to the difficulty of getting pure culture isolates. The combination of single-cell genomics and metagenomics has become an alternative approach to study these kinds of microbes. Our results expand the understanding of the genomic diversity, distribution, and lifestyles within this clade and provide evidence of different evolutionary trajectories in the genomic makeup of the two subclades that could serve to illustrate how evolutionary pressure can drive different adaptations to the same environment. Therefore, the SAR116 clade represents an ideal model organism for the study of the evolutionary streamlining of genomes in microbes that have relatively close relatedness to each other.
Collapse
|
18
|
Evolutionary stasis of a deep subsurface microbial lineage. THE ISME JOURNAL 2021; 15:2830-2842. [PMID: 33824425 PMCID: PMC8443664 DOI: 10.1038/s41396-021-00965-3] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 02/26/2021] [Accepted: 03/15/2021] [Indexed: 02/01/2023]
Abstract
Sulfate-reducing bacteria Candidatus Desulforudis audaxviator (CDA) were originally discovered in deep fracture fluids accessed via South African gold mines and have since been found in geographically widespread deep subsurface locations. In order to constrain models for subsurface microbial evolution, we compared CDA genomes from Africa, North America and Eurasia using single cell genomics. Unexpectedly, 126 partial single amplified genomes from the three continents, a complete genome from of an isolate from Eurasia, and metagenome-assembled genomes from Africa and Eurasia shared >99.2% average nucleotide identity, low frequency of SNP's, and near-perfectly conserved prophages and CRISPRs. Our analyses reject sample cross-contamination, recent natural dispersal, and unusually strong purifying selection as likely explanations for these unexpected results. We therefore conclude that the analyzed CDA populations underwent only minimal evolution since their physical separation, potentially as far back as the breakup of Pangea between 165 and 55 Ma ago. High-fidelity DNA replication and repair mechanisms are the most plausible explanation for the highly conserved genome of CDA. CDA presents a stark contrast to the current model organisms in microbial evolutionary studies, which often develop adaptive traits over far shorter periods of time.
Collapse
|
19
|
Hoffert M, Anderson RE, Reveillaud J, Murphy LG, Stepanauskas R, Huber JA. Genomic Variation Influences Methanothermococcus Fitness in Marine Hydrothermal Systems. Front Microbiol 2021; 12:714920. [PMID: 34489903 PMCID: PMC8417812 DOI: 10.3389/fmicb.2021.714920] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 07/31/2021] [Indexed: 11/13/2022] Open
Abstract
Hydrogenotrophic methanogens are ubiquitous chemoautotrophic archaea inhabiting globally distributed deep-sea hydrothermal vent ecosystems and associated subseafloor niches within the rocky subseafloor, yet little is known about how they adapt and diversify in these habitats. To determine genomic variation and selection pressure within methanogenic populations at vents, we examined five Methanothermococcus single cell amplified genomes (SAGs) in conjunction with 15 metagenomes and 10 metatranscriptomes from venting fluids at two geochemically distinct hydrothermal vent fields on the Mid-Cayman Rise in the Caribbean Sea. We observed that some Methanothermococcus lineages and their transcripts were more abundant than others in individual vent sites, indicating differential fitness among lineages. The relative abundances of lineages represented by SAGs in each of the samples matched phylogenetic relationships based on single-copy universal genes, and genes related to nitrogen fixation and the CRISPR/Cas immune system were among those differentiating the clades. Lineages possessing these genes were less abundant than those missing that genomic region. Overall, patterns in nucleotide variation indicated that the population dynamics of Methanothermococcus were not governed by clonal expansions or selective sweeps, at least in the habitats and sampling times included in this study. Together, our results show that although specific lineages of Methanothermococcus co-exist in these habitats, some outcompete others, and possession of accessory metabolic functions does not necessarily provide a fitness advantage in these habitats in all conditions. This work highlights the power of combining single-cell, metagenomic, and metatranscriptomic datasets to determine how evolution shapes microbial abundance and diversity in hydrothermal vent ecosystems.
Collapse
Affiliation(s)
- Michael Hoffert
- Biology Department, Carleton College, Northfield, MN, United States.,Finch Therapeutics Group, Somerville, MA, United States
| | - Rika E Anderson
- Biology Department, Carleton College, Northfield, MN, United States
| | - Julie Reveillaud
- Maladies Infectieuses et Vecteurs: Ecologie, Génétique, Evolution et Contrôle, University of Montpellier, Institut National de la Recherche Agronomique, Centre National de la Recherche Scientifique, Institut de Recherche Pour le Développement, Montpellier, France
| | | | | | - Julie A Huber
- Marine Chemistry and Geochemistry, Woods Hole Oceanographic Institution, Woods Hole, MA, United States
| |
Collapse
|
20
|
Banerjee R, Chaudhari NM, Lahiri A, Gautam A, Bhowmik D, Dutta C, Chattopadhyay S, Huson DH, Paul S. Interplay of Various Evolutionary Modes in Genome Diversification and Adaptive Evolution of the Family Sulfolobaceae. Front Microbiol 2021; 12:639995. [PMID: 34248865 PMCID: PMC8267890 DOI: 10.3389/fmicb.2021.639995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Accepted: 05/06/2021] [Indexed: 11/21/2022] Open
Abstract
Sulfolobaceae family, comprising diverse thermoacidophilic and aerobic sulfur-metabolizing Archaea from various geographical locations, offers an ideal opportunity to infer the evolutionary dynamics across the members of this family. Comparative pan-genomics coupled with evolutionary analyses has revealed asymmetric genome evolution within the Sulfolobaceae family. The trend of genome streamlining followed by periods of differential gene gains resulted in an overall genome expansion in some species of this family, whereas there was reduction in others. Among the core genes, both Sulfolobus islandicus and Saccharolobus solfataricus showed a considerable fraction of positively selected genes and also higher frequencies of gene acquisition. In contrast, Sulfolobus acidocaldarius genomes experienced substantial amount of gene loss and strong purifying selection as manifested by relatively lower genome size and higher genome conservation. Central carbohydrate metabolism and sulfur metabolism coevolved with the genome diversification pattern of this archaeal family. The autotrophic CO2 fixation with three significant positively selected enzymes from S. islandicus and S. solfataricus was found to be more imperative than heterotrophic CO2 fixation for Sulfolobaceae. Overall, our analysis provides an insight into the interplay of various genomic adaptation strategies including gene gain-loss, mutation, and selection influencing genome diversification of Sulfolobaceae at various taxonomic levels and geographical locations.
Collapse
Affiliation(s)
- Rachana Banerjee
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
| | - Narendrakumar M. Chaudhari
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
| | - Abhishake Lahiri
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| | - Anupam Gautam
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Kolkata, India
| | - Debaleena Bhowmik
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| | - Chitra Dutta
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
| | - Sujay Chattopadhyay
- JIS Institute of Advanced Studies and Research, JIS University, Kolkata, India
| | - Daniel H. Huson
- Institute for Bioinformatics and Medical Informatics, University of Tübingen, Tübingen, Germany
- Cluster of Excellence: Controlling Microbes to Fight Infection, Tübingen, Germany
| | - Sandip Paul
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| |
Collapse
|
21
|
Graham ED, Tully BJ. Marine Dadabacteria exhibit genome streamlining and phototrophy-driven niche partitioning. THE ISME JOURNAL 2021; 15:1248-1256. [PMID: 33230264 PMCID: PMC8115339 DOI: 10.1038/s41396-020-00834-5] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Revised: 10/27/2020] [Accepted: 11/05/2020] [Indexed: 01/29/2023]
Abstract
The remineralization of organic material via heterotrophy in the marine environment is performed by a diverse and varied group of microorganisms that can specialize in the type of organic material degraded and the niche they occupy. The marine Dadabacteria are cosmopolitan in the marine environment and belong to a candidate phylum for which there has not been a comprehensive assessment of the available genomic data to date. Here in, we assess the functional potential of the marine pelagic Dadabacteria in comparison to members of the phylum that originate from terrestrial, hydrothermal, and subsurface environments. Our analysis reveals that the marine pelagic Dadabacteria have streamlined genomes, corresponding to smaller genome sizes and lower nitrogen content of their DNA and predicted proteome, relative to their phylogenetic counterparts. Collectively, the Dadabacteria have the potential to degrade microbial dissolved organic matter, specifically peptidoglycan and phospholipids. The marine Dadabacteria belong to two clades with apparent distinct ecological niches in global metagenomic data: a clade with the potential for photoheterotrophy through the use of proteorhodopsin, present predominantly in surface waters up to 100 m depth; and a clade lacking the potential for photoheterotrophy that is more abundant in the deep photic zone.
Collapse
Affiliation(s)
- Elaina D Graham
- Department of Biological Sciences, University of Southern California, Los Angeles, CA, USA
| | - Benjamin J Tully
- Department of Biological Sciences, University of Southern California, Los Angeles, CA, USA.
- Center for Dark Energy Biosphere Investigations, University of Southern California, Los Angeles, CA, USA.
| |
Collapse
|
22
|
Moniruzzaman M, Weinheimer AR, Martinez-Gutierrez CA, Aylward FO. Widespread endogenization of giant viruses shapes genomes of green algae. Nature 2020; 588:141-145. [PMID: 33208937 DOI: 10.1038/s41586-020-2924-2] [Citation(s) in RCA: 71] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 09/01/2020] [Indexed: 12/21/2022]
Abstract
Endogenous viral elements (EVEs)-viruses that have integrated their genomes into those of their hosts-are prevalent in eukaryotes and have an important role in genome evolution1,2. The vast majority of EVEs that have been identified to date are small genomic regions comprising a few genes2, but recent evidence suggests that some large double-stranded DNA viruses may also endogenize into the genome of the host1. Nucleocytoplasmic large DNA viruses (NCLDVs) have recently become of great interest owing to their large genomes and complex evolutionary origins3-6, but it is not yet known whether they are a prominent component of eukaryotic EVEs. Here we report the widespread endogenization of NCLDVs in diverse green algae; these giant EVEs reached sizes greater than 1 million base pairs and contained as many as around 10% of the total open reading frames in some genomes, substantially increasing the scale of known viral genes in eukaryotic genomes. These endogenized elements often shared genes with host genomic loci and contained numerous spliceosomal introns and large duplications, suggesting tight assimilation into host genomes. NCLDVs contain large and mosaic genomes with genes derived from multiple sources, and their endogenization represents an underappreciated conduit of new genetic material into eukaryotic lineages that can substantially impact genome composition.
Collapse
Affiliation(s)
| | | | | | - Frank O Aylward
- Department of Biological Sciences, Virginia Tech, Blacksburg, VA, USA.
| |
Collapse
|