Genomic insights into versatile lifestyle of three new bacterial candidate phyla

Metagenomic explorations of the Earth’s biosphere enable the discovery of previously unknown bacterial lineages of phylogenetic and ecological significance. Here, we retrieved 11 metagenomic-assembled genomes (MAGs) affiliated to three new monophyletic bacterial lineages from the seawater of the Yap Trench. Phylogenomic analysis revealed that each lineage is a new bacterial candidate phylum, subsequently named Candidatus Qinglongiota, Candidatus Heilongiota, and Candidatus Canglongiota. Metabolic reconstruction of genomes from the three phyla suggested that they adopt a versatile lifestyle, with the potential to utilize various types of sugars, proteins, and/or short-chain fatty acids through anaerobic pathways. This was further confirmed by a global distribution map of the three phyla, indicating a preference for oxygen-limited or particle-attached niches, such as anoxic sedimentary environments. Of note, Candidatus Canglongiota genomes harbor genes for the complete Wood- Ljungdahl pathway and sulfate reduction that are similar to those identified in some sulfate-reducing bacteria. Evolutionary analysis indicated that gene gain and loss events, and horizontal gene transfer (HGT) play important roles in shaping the genomic and metabolic features of the three new phyla. This study presents the genomic insight into the ecology, metabolism, and evolution of three new phyla, which broadens the phylum-level diversity within the domain Bacteria.

     

Reductions in complexity of mitochondrial genomes in lichen‐forming fungi shed light on genome architecture of obligate symbioses

Symbioses among co‐evolving taxa are often marked by genome reductions such as a loss of protein‐coding genes in at least one of the partners as a means of reducing redundancy or intergenomic conflict. To explore this phenomenon in an iconic yet under‐studied group of obligate symbiotic organisms, mitochondrial genomes of 22 newly sequenced and annotated species of lichenized fungi were compared to 167 mitochondrial genomes of nonlichenized fungi. Our results demonstrate the first broad‐scale loss of atp9 from mitochondria of lichenized fungi. Despite key functions in mitochondrial energy production, we show that atp9 has been independently lost in three different lineages spanning 10 of the 22 studied species. A search for predicted, functional copies of atp9 among genomes of other symbionts involved in each lichen revealed the full‐length, presumably functional copies of atp9 in either the photosynthetic algal partner or in other symbiotic fungi in all 10 instances. Together, these data yield evidence of an obligate symbiotic relationship in which core genomic processes have been streamlined, likely due to co‐evolution.     

Gulf of Mexico blue hole harbors high levels of novel microbial lineages

Exploration of oxygen-depleted marine environments has consistently revealed novel microbial taxa and metabolic capabilities that expand our understanding of microbial evolution and ecology. Marine blue holes are shallow karst formations characterized by low oxygen and high organic matter content. They are logistically challenging to sample, and thus our understanding of their biogeochemistry and microbial ecology is limited. We present a metagenomic and geochemical characterization of Amberjack Hole on the Florida continental shelf (Gulf of Mexico). Dissolved oxygen became depleted at the hole’s rim (32 m water depth), remained low but detectable in an intermediate hypoxic zone (40–75 m), and then increased to a secondary peak before falling below detection in the bottom layer (80–110 m), concomitant with increases in nutrients, dissolved iron, and a series of sequentially more reduced sulfur species. Microbial communities in the bottom layer contained heretofore undocumented levels of the recently discovered phylum Woesearchaeota (up to 58% of the community), along with lineages in the bacterial Candidate Phyla Radiation (CPR). Thirty-one high-quality metagenome-assembled genomes (MAGs) showed extensive biochemical capabilities for sulfur and nitrogen cycling, as well as for resisting and respiring arsenic. One uncharacterized gene associated with a CPR lineage differentiated hypoxic from anoxic zone communities. Overall, microbial communities and geochemical profiles were stable across two sampling dates in the spring and fall of 2019. The blue hole habitat is a natural marine laboratory that provides opportunities for sampling taxa with under-characterized but potentially important roles in redox-stratified microbial processes.Subject terms: Metagenomics, Microbial ecology, Biogeochemistry     

候选门级辐射类群(CPR)细菌研究进展

候选门级辐射类群（candidate phyla radiation,CPR）是一类在自然界广泛分布的庞大细菌类群,成员众多,基于宏基因组测序结果已划分出75个门水平分支。CPR细菌是研究微生物进化与各类生境微生物协同作用的理想材料与系统,其在系统进化树中位于细菌域的最外侧,与古菌域外侧的DPANN超门毗邻,进化地位特殊。该类群具有独特的生物学特性,如细胞个体小,基因组精简,缺少三羧酸循环途径以及常见细菌电子传递链相关基因,通常需要与宿主共生、纯培养困难,目前仅在人类口腔中成功分离到一株CPR细菌（TM7菌株UB2523）。本文综述了CPR细菌的发现历程、细胞形态、代谢特征、在碳/氮/硫元素循环中的作用,以及不依赖于培养的组学技术在CPR细菌研究中的应用,可为深入探究CPR细菌的进化地位和驱动地球化学元素循环的方式,以及与宿主共生系统的互作机制提供参考。     

An interdependent metabolic patchwork in the nested symbiosis of mealybugs

Highly reduced genomes of 144-416 kilobases have been described from nutrient-provisioning bacterial symbionts of several insect lineages [1-5]. Some host insects have formed stable associations with pairs of bacterial symbionts that live in specialized cells and provide them with essential nutrients; genomic data from these systems have revealed remarkable levels of metabolic complementarity between the symbiont pairs [3, 4, 6, 7]. The mealybug Planococcus citri (Hemiptera: Pseudococcidae) contains dual bacterial symbionts existing with an unprecedented organization: an unnamed gammaproteobacteria, for which we propose the name Candidatus Moranella endobia, lives inside the betaproteobacteria Candidatus Tremblaya princeps [8]. Here we describe the complete genomes and metabolic contributions of these unusual nested symbionts. We show that whereas there is little overlap in retained genes involved in nutrient production between symbionts, several essential amino acid pathways in the mealybug assemblage require a patchwork of interspersed gene products from Tremblaya, Moranella, and possibly P.?citri. Furthermore, although Tremblaya has the smallest cellular genome yet described, it contains a genomic inversion present in both orientations in individual insects, starkly contrasting with the extreme structural stability typical of highly reduced bacterial genomes [4, 9, 10].     

Metabolic networks of Sodalis glossinidius: a systems biology approach to reductive evolution

Background

Genome reduction is a common evolutionary process affecting bacterial lineages that establish symbiotic or pathogenic associations with eukaryotic hosts. Such associations yield highly reduced genomes with greatly streamlined metabolic abilities shaped by the type of ecological association with the host. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, represents one of the few complete genomes available of a bacterium at the initial stages of this process. In the present study, genome reduction is studied from a systems biology perspective through the reconstruction and functional analysis of genome-scale metabolic networks of S. glossinidius.

Results

The functional profile of ancestral and extant metabolic networks sheds light on the evolutionary events underlying transition to a host-dependent lifestyle. Meanwhile, reductive evolution simulations on the extant metabolic network can predict possible future evolution of S. glossinidius in the context of genome reduction. Finally, knockout simulations in different metabolic systems reveal a gradual decrease in network robustness to different mutational events for bacterial endosymbionts at different stages of the symbiotic association.

Conclusions

Stoichiometric analysis reveals few gene inactivation events whose effects on the functionality of S. glossinidius metabolic systems are drastic enough to account for the ecological transition from a free-living to host-dependent lifestyle. The decrease in network robustness across different metabolic systems may be associated with the progressive integration in the more stable environment provided by the insect host. Finally, reductive evolution simulations reveal the strong influence that external conditions exert on the evolvability of metabolic systems.     

Comparison of the mesophilic cellulosome‐producing Clostridium cellulovorans genome with other cellulosome‐related clostridial genomes

Clostridium cellulovorans, an anaerobic and mesophilic bacterium, degrades native substrates in soft biomass such as corn fibre and rice straw efficiently by producing an extracellular enzyme complex called the cellulosome. Recently, we have reported the whole‐genome sequence of C. cellulovorans comprising 4220 predicted genes in 5.10 Mbp [Y. Tamaru et al., (2010) J. Bacteriol., 192: 901–902]. As a result, the genome size of C. cellulovorans was about 1 Mbp larger than that of other cellulosome‐producing clostridia, mesophilic C. cellulolyticum and thermophilic C. thermocellum. A total of 57 cellulosomal genes were found in the C. cellulovorans genome, and they coded for not only carbohydrate‐degrading enzymes but also a lipase, peptidases and proteinase inhibitors. Interestingly, two novel genes encoding scaffolding proteins were found in the genome. According to KEGG metabolic pathways and their comparison with 11 Clostridial genomes, gene expansion in the C. cellulovorans genome indicated mainly non‐cellulosomal genes encoding hemicellulases and pectin‐degrading enzymes. Thus, by examining genome sequences from multiple Clostridium species, comparative genomics offers new insight into genome evolution and the way natural selection moulds functional DNA sequence evolution. Our analysis, coupled with the genome sequence data, provides a roadmap for constructing enhanced cellulosome‐producing Clostridium strains for industrial applications such as biofuel production.     

Iterative subtractive binning of freshwater chronoseries metagenomes identifies over 400 novel species and their ecologic preferences

Recent advances in sequencing technology and bioinformatic pipelines have allowed unprecedented access to the genomes of yet-uncultivated microorganisms from diverse environments. However, the catalogue of freshwater genomes remains limited, and most genome recovery attempts in freshwater ecosystems have only targeted specific taxa. Here, we present a genome recovery pipeline incorporating iterative subtractive binning, and apply it to a time series of 100 metagenomic datasets from seven connected lakes and estuaries along the Chattahoochee River (Southeastern USA). Our set of metagenome-assembled genomes (MAGs) represents >400 yet-unnamed genomospecies, substantially increasing the number of high-quality MAGs from freshwater lakes. We propose names for two novel species: ‘Candidatus Elulimicrobium humile’ (‘Ca. Elulimicrobiota’, ‘Patescibacteria’) and ‘Candidatus Aquidulcis frankliniae’ (‘Chloroflexi’). Collectively, our MAGs represented about half of the total microbial community at any sampling point. To evaluate the prevalence of these genomospecies in the chronoseries, we introduce methodologies to estimate relative abundance and habitat preference that control for uneven genome quality and sample representation. We demonstrate high degrees of habitat-specialization and endemicity for most genomospecies in the Chattahoochee lakes. Wider ecological ranges characterized smaller genomes with higher coding densities, indicating an overall advantage of smaller, more compact genomes for cosmopolitan distributions.     

Polyploidy does not control all: Lineage‐specific average chromosome length constrains genome size evolution in ferns

Recent studies investigating the evolution of genome size diversity in ferns have shown that they have a distinctive genome profile compared with other land plants. Ferns are typically characterized by possessing medium‐sized genomes, although a few lineages have evolved very large genomes. Ferns are different from other vascular plant lineages as they are the only group to show evidence for a correlation between genome size and chromosome number. In this study, we aim to explore whether the evolution of fern genome sizes is not only shaped by chromosome number changes arising from polyploidy but also by constraints on the average amount of DNA per chromosome. We selected the genus Asplenium L. as a model genus to study the question because of the unique combination of a highly conserved base chromosome number and a high frequency of polyploidy. New genome size data for Asplenium taxa were combined with existing data and analyzed within a phylogenetic framework. Genome size varied substantially between diploid species, resulting in overlapping genome sizes among diploid and tetraploid spleenworts. The observed additive pattern indicates the absence of genome downsizing following polyploidy. The genome size of diploids varied non‐randomly and we found evidence for clade‐specific trends towards larger or smaller genomes. The 578‐fold range of fern genome sizes have arisen not only from repeated cycles of polyploidy but also through clade‐specific constraints governing accumulation and/or elimination of DNA.     

Single‐cell sequencing unveils the lifestyle and CRISPR‐based population history of Hydrotalea sp. in acid mine drainage

Acid mine drainage (AMD) is characterized by an acid and metal‐rich run‐off that originates from mining systems. Despite having been studied for many decades, much remains unknown about the microbial community dynamics in AMD sites, especially during their early development, when the acidity is moderate. Here, we describe draft genome assemblies from single cells retrieved from an early‐stage AMD sample. These cells belong to the genus Hydrotalea and are closely related to Hydrotalea flava. The phylogeny and average nucleotide identity analysis suggest that all single amplified genomes (SAGs) form two clades that may represent different strains. These cells have the genomic potential for denitrification, copper and other metal resistance. Two coexisting CRISPR‐Cas loci were recovered across SAGs, and we observed heterogeneity in the population with regard to the spacer sequences, together with the loss of trailer‐end spacers. Our results suggest that the genomes of Hydrotalea sp. strains studied here are adjusting to a quickly changing selective pressure at the microhabitat scale, and an important form of this selective pressure is infection by foreign DNA.     

Genome‐wide sequence information reveals recurrent hybridization among diploid wheat wild relatives

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome‐wide sequence information, independent data were obtained from genotyping‐by‐sequencing and a target‐enrichment experiment that returned 244 low‐copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent‐based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat‐group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat‐group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.     

A method for achieving complete microbial genomes and improving bins from metagenomics data

Metagenomics facilitates the study of the genetic information from uncultured microbes and complex microbial communities. Assembling complete genomes from metagenomics data is difficult because most samples have high organismal complexity and strain diversity. Some studies have attempted to extract complete bacterial, archaeal, and viral genomes and often focus on species with circular genomes so they can help confirm completeness with circularity. However, less than 100 circularized bacterial and archaeal genomes have been assembled and published from metagenomics data despite the thousands of datasets that are available. Circularized genomes are important for (1) building a reference collection as scaffolds for future assemblies, (2) providing complete gene content of a genome, (3) confirming little or no contamination of a genome, (4) studying the genomic context and synteny of genes, and (5) linking protein coding genes to ribosomal RNA genes to aid metabolic inference in 16S rRNA gene sequencing studies. We developed a semi-automated method called Jorg to help circularize small bacterial, archaeal, and viral genomes using iterative assembly, binning, and read mapping. In addition, this method exposes potential misassemblies from k-mer based assemblies. We chose species of the Candidate Phyla Radiation (CPR) to focus our initial efforts because they have small genomes and are only known to have one ribosomal RNA operon. In addition to 34 circular CPR genomes, we present one circular Margulisbacteria genome, one circular Chloroflexi genome, and two circular megaphage genomes from 19 public and published datasets. We demonstrate findings that would likely be difficult without circularizing genomes, including that ribosomal genes are likely not operonic in the majority of CPR, and that some CPR harbor diverged forms of RNase P RNA. Code and a tutorial for this method is available at https://github.com/lmlui/Jorg and is available on the DOE Systems Biology KnowledgeBase as a beta app.     

Ecological and genomic analyses of candidate phylum WPS-2 bacteria in an unvegetated soil

Members of the bacterial candidate phylum WPS-2 (or Eremiobacterota) are abundant in several dry, bare soil environments. In a bare soil deposited by an extinct iron–sulfur spring, we found that WPS-2 comprised up to 24% of the bacterial community and up to 10⁸ cells per g of soil based on 16S rRNA gene sequencing and quantification. A single genus-level cluster (Ca. Rubrimentiphilum) predominated in bare soils but was less abundant in adjacent forest. Nearly complete genomes of Ca. Rubrimentiphilum were recovered as single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs). Surprisingly, given the abundance of WPS-2 in bare soils, the genomes did not indicate any capacity for autotrophy, phototrophy, or trace gas metabolism. Instead, they suggest a predominantly aerobic organoheterotrophic lifestyle, perhaps based on scavenging amino acids, nucleotides, and complex oligopeptides, along with lithotrophic capacity on thiosulfate. Network analyses of the entire community showed that some species of Chloroflexi, Actinobacteria, and candidate phylum AD3 (or Dormibacterota) co-occurred with Ca. Rubrimentiphilum and may represent ecological or metabolic partners. We propose that Ca. Rubrimentiphilum act as efficient heterotrophic scavengers. Combined with previous studies, these data suggest that the phylum WPS-2 includes bacteria with diverse metabolic capabilities.     

Plasticity and epistasis strongly affect bacterial fitness after losing multiple metabolic genes

Many bacterial lineages lack seemingly essential metabolic genes. Previous work suggested selective benefits could drive the loss of biosynthetic functions from bacterial genomes when the corresponding metabolites are sufficiently available in the environment. However, the factors that govern this “genome streamlining” remain poorly understood. Here we determine the effect of plasticity and epistasis on the fitness of Escherichia coli genotypes from whose genome biosynthetic genes for one, two, or three different amino acids have been deleted. Competitive fitness experiments between auxotrophic mutants and prototrophic wild‐type cells in one of two carbon environments revealed that plasticity and epistasis strongly affected the mutants’ fitness individually and interactively. Positive and negative epistatic interactions were prevalent, yet on average cancelled each other out. Moreover, epistasis correlated negatively with the expected effects of combined auxotrophy‐causing mutations, thus producing a pattern of diminishing returns. Moreover, computationally analyzing 1,432 eubacterial metabolic networks revealed that most pairs of auxotrophies co‐occurred significantly more often than expected by chance, suggesting epistatic interactions and/or environmental factors favored these combinations. Our results demonstrate that both the genetic background and environmental conditions determine the adaptive value of a loss‐of‐biochemical‐function mutation and that fitness gains decelerate, as more biochemical functions are lost.     

Single‐cell genomics of uncultivated deep‐branching magnetotactic bacteria reveals a conserved set of magnetosome genes

While magnetosome biosynthesis within the magnetotactic Proteobacteria is increasingly well understood, much less is known about the genetic control within deep‐branching phyla, which have a unique ultrastructure and biosynthesize up to several hundreds of bullet‐shaped magnetite magnetosomes arranged in multiple bundles of chains, but have no cultured representatives. Recent metagenomic analysis identified magnetosome genes in the genus ‘Candidatus Magnetobacterium’ homologous to those in Proteobacteria. However, metagenomic analysis has been limited to highly abundant members of the community, and therefore only little is known about the magnetosome biosynthesis, ecophysiology and metabolic capacity in deep‐branching MTB. Here we report the analysis of single‐cell derived draft genomes of three deep‐branching uncultivated MTB. Single‐cell sorting followed by whole genome amplification generated draft genomes of Candidatus Magnetobacterium bavaricum and Candidatus Magnetoovum chiemensis CS‐04 of the Nitrospirae phylum. Furthermore, we present the first, nearly complete draft genome of a magnetotactic representative from the candidate phylum Omnitrophica, tentatively named Candidatus Omnitrophus magneticus SKK‐01. Besides key metabolic features consistent with a common chemolithoautotrophic lifestyle, we identified numerous, partly novel genes most likely involved in magnetosome biosynthesis of bullet‐shaped magnetosomes and their arrangement in multiple bundles of chains.     

Genome and species relationships in genus<Emphasis Type="Italic"> Avena</Emphasis> based on RAPD and AFLP molecular markers

Species and genome relationships among 11 diploid (A and C genomes), five tetraploid (AB and AC genomes) and two hexaploid (ACD genome) Avena taxa were investigated using amplified fragment length polymorphisms (AFLPs) and random amplified polymorphic DNA (RAPD) markers. The two primer pairs used for the AFLP reactions produced a total of 354 polymorphic bands, while 187 reproducible bands were generated using ten RAPD primers. Genetic similarities amongst the entries were estimated using the Jaccard and Dice algorithms, and cluster analyses were performed using UPGMA and neighbor joining methods. Principle coordinate analysis was also applied. The highest cophenetic correlation coefficient was obtained for the Jaccard algorithm and UPGMA clustering method (r=0.99 for AFLP and r=0.94 for RAPD). No major clustering differences were present between phenograms produced with AFLPs and RAPDs. Furthermore, data produced with AFLPs and RAPDs were highly correlated (r=0.92), indicating the reliability of our results. All A genome diploid taxa are clustered together according to their karyotype. The AB genome tetraploids were found to form a subcluster within the A_s genome diploids (AFLPs), indicating their near-autoploid origin. The AC genome tetraploids are clustered to the ACD genome hexaploids. Finally, the C genome diploids form an outer branch, indicating the major genomic divergence between the A and C genomes in Avena.Communicated by J.S. Heslop-Harrison