Integrated mobile genetic elements in Thaumarchaeota

To explore the diversity of mobile genetic elements (MGE) associated with archaea of the phylum Thaumarchaeota, we exploited the property of most MGE to integrate into the genomes of their hosts. Integrated MGE (iMGE) were identified in 20 thaumarchaeal genomes amounting to 2 Mbp of mobile thaumarchaeal DNA. These iMGE group into five major classes: (i) proviruses, (ii) casposons, (iii) insertion sequence-like transposons, (iv) integrative-conjugative elements and (v) cryptic integrated elements. The majority of the iMGE belong to the latter category and might represent novel families of viruses or plasmids. The identified proviruses are related to tailed viruses of the order Caudovirales and to tailless icosahedral viruses with the double jelly-roll capsid proteins. The thaumarchaeal iMGE are all connected within a gene sharing network, highlighting pervasive gene exchange between MGE occupying the same ecological niche. The thaumarchaeal mobilome carries multiple auxiliary metabolic genes, including multicopper oxidases and ammonia monooxygenase subunit C (AmoC), and stress response genes, such as those for universal stress response proteins (UspA). Thus, iMGE might make important contributions to the fitness and adaptation of their hosts. We identified several iMGE carrying type I-B CRISPR-Cas systems and spacers matching other thaumarchaeal iMGE, suggesting antagonistic interactions between coexisting MGE and symbiotic relationships with the ir archaeal hosts.     

A curated data resource of 214K metagenomes for characterization of the global antimicrobial resistome

The growing threat of antimicrobial resistance (AMR) calls for new epidemiological surveillance methods, as well as a deeper understanding of how antimicrobial resistance genes (ARGs) have been transmitted around the world. The large pool of sequencing data available in public repositories provides an excellent resource for monitoring the temporal and spatial dissemination of AMR in different ecological settings. However, only a limited number of research groups globally have the computational resources to analyze such data. We retrieved 442 Tbp of sequencing reads from 214,095 metagenomic samples from the European Nucleotide Archive (ENA) and aligned them using a uniform approach against ARGs and 16S/18S rRNA genes. Here, we present the results of this extensive computational analysis and share the counts of reads aligned. Over 6.76∙10⁸ read fragments were assigned to ARGs and 3.21∙10⁹ to rRNA genes, where we observed distinct differences in both the abundance of ARGs and the link between microbiome and resistome compositions across various sampling types. This collection is another step towards establishing global surveillance of AMR and can serve as a resource for further research into the environmental spread and dynamic changes of ARGs.

The growing threat of antimicrobial resistance (AMR) calls for new epidemiological surveillance methods and a deeper understanding of how resistance genes are transmitted around the world. This study presents a large-scale remapping of sequencing reads of publicly available metagenomic datasets that can be used to monitor the global prevalence of AMR genes.     

Nature of factors stimulating cobalamin production by Achromobacter cobalamini

The object of this work was to study the nature of factors contained in molasses and maize extract and stimulating cobalaminogenesis in Achromobacter cobalamini. The activity of substrate fractions was analyzed to show that the stimulating substance was precipitated on the cation exchanger and eluted from it with HCl. The factor was found to be an organic nitrogen base readily soluble in water and ethanol but insoluble in ether, chloroform and methanol. It was stable upon heating in concentrated HCl. Betaine in the composition of molasses and choline in the composition of maize extract had similar properties. Their addition to the growth medium produced the same effect as that of molasses and maize extract. It is concluded therefore that cobalaminogenesis is stimulated in A. cobalamini by betaine in molasses and by choline in maize extract.     

Phylogenetic congruence and ecological coherence in terrestrial Thaumarchaeota

Thaumarchaeota form a ubiquitously distributed archaeal phylum, comprising both the ammonia-oxidising archaea (AOA) and other archaeal groups in which ammonia oxidation has not been demonstrated (including Group 1.1c and Group 1.3). The ecology of AOA in terrestrial environments has been extensively studied using either a functional gene, encoding ammonia monooxygenase subunit A (amoA) or 16S ribosomal RNA (rRNA) genes, which show phylogenetic coherence with respect to soil pH. To test phylogenetic congruence between these two markers and to determine ecological coherence in all Thaumarchaeota, we performed high-throughput sequencing of 16S rRNA and amoA genes in 46 UK soils presenting 29 available contextual soil characteristics. Adaptation to pH and organic matter content reflected strong ecological coherence at various levels of taxonomic resolution for Thaumarchaeota (AOA and non-AOA), whereas nitrogen, total mineralisable nitrogen and zinc concentration were also important factors associated with AOA thaumarchaeotal community distribution. Other significant associations with environmental factors were also detected for amoA and 16S rRNA genes, reflecting different diversity characteristics between these two markers. Nonetheless, there was significant statistical congruence between the markers at fine phylogenetic resolution, supporting the hypothesis of low horizontal gene transfer between Thaumarchaeota. Group 1.1c Thaumarchaeota were also widely distributed, with two clusters predominating, particularly in environments with higher moisture content and organic matter, whereas a similar ecological pattern was observed for Group 1.3 Thaumarchaeota. The ecological and phylogenetic congruence identified is fundamental to understand better the life strategies, evolutionary history and ecosystem function of the Thaumarchaeota.     

Metabolic engineering of cobalamin (vitamin B12) production in Bacillus megaterium

Cobalamin (vitamin B₁₂) production in Bacillus megaterium has served as a model system for the systematic evaluation of single and multiple directed molecular and genetic optimization strategies. Plasmid and genome-based overexpression of genes involved in vitamin B₁₂ biosynthesis, including cbiX, sirA, modified hemA, the operons hemAXCDBL and cbiXJCDETLFGAcysG^AcbiYbtuR, and the regulatory gene fnr, significantly increased cobalamin production. To reduce flux along the heme branch of the tetrapyrrole pathway, an antisense RNA strategy involving silencing of the hemZ gene encoding coproporphyrinogen III oxidase was successfully employed. Feedback inhibition of the initial enzyme of the tetrapyrrole biosynthesis, HemA, by heme was overcome by stabilized enzyme overproduction. Similarly, the removal of the B₁₂ riboswitch upstream of the cbiXJCDETLFGAcysG^AcbiYbtuR operon and the recombinant production of three different vitamin B₁₂ binding proteins (glutamate mutase GlmS, ribonucleotide triphosphate reductase RtpR and methionine synthase MetH) partly abolished B₁₂-dependent feedback inhibition. All these strategies increased cobalamin production in B. megaterium. Finally, combinations of these strategies enhanced the overall intracellular vitamin B₁₂ concentrations but also reduced the volumetric cellular amounts by placing the organism under metabolic stress.     

Iron transporters in marine prokaryotic genomes and metagenomes

In the pelagic environment, iron is a scarce but essential micronutrient. The iron acquisition capabilities of selected marine bacteria have been investigated, but the recent proliferation of marine prokaryotic genomes and metagenomes offers a more comprehensive picture of microbial iron uptake pathways in the ocean. Searching these data sets, we were able to identify uptake mechanisms for Fe(3+), Fe(2+) and iron chelates (e.g. siderophore and haem iron complexes). Transport of iron chelates is accomplished by TonB-dependent transporters (TBDTs). After clustering the TBDTs from marine prokaryotic genomes, we identified TBDT clusters for the transport of hydroxamate and catecholate siderophore iron complexes and haem using gene neighbourhood analysis and co-clustering of TBDTs of known function. The genomes also contained two classes of siderophore biosynthesis genes: NRPS (non-ribosomal peptide synthase) genes and NIS (NRPS Independent Siderophore) genes. The most common iron transporters, in both the genomes and metagenomes, were Fe(3+) ABC transporters. Iron uptake-related TBDTs and siderophore biosynthesis genes were less common in pelagic marine metagenomes relative to the genomic data set, in part because Pelagibacter ubique and Prochlorococcus species, which almost entirely lacked these Fe uptake systems, dominate the metagenomes. Our results are largely consistent with current knowledge of iron speciation in the ocean, but suggest that in certain niches the ability to acquire siderophores and/or haem iron chelates is beneficial.     

Global phylogeography of chitinase genes in aquatic metagenomes

Phylogeny-based analysis of chitinase and 16S rRNA genes from metagenomic data suggests that salinity is a major driver for the distribution of both chitinolytic and total bacterial communities in aquatic systems. Additionally, more acidic chitinase proteins were observed with increasing salinity. Congruent habitat separation was further observed for both genes according to latitude and proximity to the coastline. However, comparison of chitinase and 16S rRNA genes extracted from different geographic locations showed little congruence in distribution. There was no indication that dispersal limited the global distribution of either gene.     

Advances in recovery of novel biocatalysts from metagenomes

Metagenomics has accelerated the process of discovery of novel biocatalysts by enabling scientists to tap directly into the entire diversity of enzymes held within natural microbial populations. Their characterization has revealed a great deal of valuable information about enzymatic activity in terms of factors which influence their stability and activity under a wide range of conditions. Many of the biocatalysts have particular properties making them suitable for biotechnological applications. A diverse array of strategies has been developed to optimize each step of the process of generating and screening metagenomic libraries for novel biocatalysts. This review covers the diversity of metagenome-derived enzymes characterized to date, and the strategies currently being developed to optimize discovery of novel metagenomic biocatalysts.     

Sequencing effort dictates gene discovery in marine microbial metagenomes

Massive metagenomic sequencing combined with gene prediction methods were previously used to compile the gene catalogue of the ocean and host-associated microbes. Global expeditions conducted over the past 15 years have sampled the ocean to build a catalogue of genes from pelagic microbes. Here we undertook a large sequencing effort of a perturbed Red Sea plankton community to uncover that the rate of gene discovery increases continuously with sequencing effort, with no indication that the retrieved 2.83 million non-redundant (complete) genes predicted from the experiment represented a nearly complete inventory of the genes present in the sampled community (i.e., no evidence of saturation). The underlying reason is the Pareto-like distribution of the abundance of genes in the plankton community, resulting in a very long tail of millions of genes present at remarkably low abundances, which can only be retrieved through massive sequencing. Microbial metagenomic projects retrieve a variable number of unique genes per Tera base-pair (Tbp), with a median value of 14.7 million unique genes per Tbp sequenced across projects. The increase in the rate of gene discovery in microbial metagenomes with sequencing effort implies that there is ample room for new gene discovery in further ocean and holobiont sequencing studies.     

Role of coral reefs in global ocean production

Coral reefs cover some 600 thousand square kilometres of the earth's surface (0.17% of the ocean surface). First order estimates show coral reefs to contribute about 0.05% of the estimated net CO₂ fixation rate of the global oceans. Gross CO₂ fixation is relatively high (of the order 700×10¹² g C year^-1), but most of this material is recycled within the reefs. Excess (net) production of organic material (E) is much smaller, of the order 20×10¹² g C year^-1. We estimate that 3×10¹² g C year^-1 (15% ofE) is buried in reef structure, 2×10¹² g C year^-1 (10% ofE) is available for sustained human harvest, and the remaining 75% ofE is available for export from coral reefs to adjacent areas. Comparison of estimates for net production by reefs and their surrounding oceans indicates that the excess production by coral reefs is similar to new production in the photic zone of oligotrophic oceans. Consequently, estimates for global ocean production should as a first approximation include reefal areas with the surrounding ocean when assigning average net production rates. While there are significant uncertainties in these numbers, it can be concluded that organic production by reefs plays a relatively minor role in the global scale of fluxes and storage of elements. In comparison, the companion process of biologically-mediated inorganic carbon precipitation represents a major role for reefs. While reef production does respond on local scales to variation in ocean climate, neither the absolute rates nor the amount accumulated into organic pools appear to be either sensitive indicators or accurate recorders of climatic change in most reef systems. Similarly, the productivity of most reefs should be little affected by currently predicted environmental changes resulting from the Greenhouse effect.Authorship sequence does not imply seniority, it is deliberately alphabetical     

Dynamics of Autotrophic Marine Planktonic Thaumarchaeota in the East China Sea

The ubiquitous and abundant distribution of ammonia-oxidizing Thaumarchaeota in marine environments is now well documented, and their crucial role in the global nitrogen cycle has been highlighted. However, the potential contribution of Thaumarchaeota in the carbon cycle remains poorly understood. Here we present for the first time a seasonal investigation on the shelf region (bathymetry≤200 m) of the East China Sea (ECS) involving analysis of both thaumarchaeal 16S rRNA and autotrophy-related genes (acetyl-CoA carboxylase gene, accA). Quantitative PCR results clearly showed a higher abundance of thaumarchaeal 16S and accA genes in late-autumn (November) than summer (August), whereas the diversity and community structure of autotrophic Thaumarchaeota showed no statistically significant difference between different seasons as revealed by thaumarchaeal accA gene clone libraries. Phylogenetic analysis indicated that shallow ecotypes dominated the autotrophic Thaumarchaeota in the ECS shelf (86.3% of total sequences), while a novel non-marine thaumarchaeal accA lineage was identified in the Changjiang estuary in summer (when freshwater plumes become larger) but not in autumn, implying that Changjiang freshwater discharge played a certain role in transporting terrestrial microorganisms to the ECS. Multivariate statistical analysis indicated that the biogeography of the autotrophic Thaumarchaeota in the shelf water of the ECS was influenced by complex hydrographic conditions. However, an in silico comparative analysis suggested that the diversity and abundance of the autotrophic Thaumarchaeota might be biased by the ‘universal’ thaumarchaeal accA gene primers Cren529F/Cren981R since this primer set is likely to miss some members within particular phylogenetic groups. Collectively, this study improved our understanding of the biogeographic patterns of the autotrophic Thaumarchaeota in temperate coastal waters, and suggested that new accA primers with improved coverage and sensitivity across phylogenetic groups are needed to gain a more thorough understanding of the role of the autotrophic Thaumarchaeota in the global carbon cycle.     

Strain/species identification in metagenomes using genome-specific markers

Shotgun metagenome sequencing has become a fast, cheap and high-throughput technology for characterizing microbial communities in complex environments and human body sites. However, accurate identification of microorganisms at the strain/species level remains extremely challenging. We present a novel k-mer-based approach, termed GSMer, that identifies genome-specific markers (GSMs) from currently sequenced microbial genomes, which were then used for strain/species-level identification in metagenomes. Using 5390 sequenced microbial genomes, 8 770 321 50-mer strain-specific and 11 736 360 species-specific GSMs were identified for 4088 strains and 2005 species (4933 strains), respectively. The GSMs were first evaluated against mock community metagenomes, recently sequenced genomes and real metagenomes from different body sites, suggesting that the identified GSMs were specific to their targeting genomes. Sensitivity evaluation against synthetic metagenomes with different coverage suggested that 50 GSMs per strain were sufficient to identify most microbial strains with ≥0.25× coverage, and 10% of selected GSMs in a database should be detected for confident positive callings. Application of GSMs identified 45 and 74 microbial strains/species significantly associated with type 2 diabetes patients and obese/lean individuals from corresponding gastrointestinal tract metagenomes, respectively. Our result agreed with previous studies but provided strain-level information. The approach can be directly applied to identify microbial strains/species from raw metagenomes, without the effort of complex data pre-processing.     

Variation in tropical reef symbiont metagenomes defined by secondary metabolism

The complex evolution of secondary metabolism is important in biology, drug development, and synthetic biology. To examine this problem at a fine scale, we compared the genomes and chemistry of 24 strains of uncultivated cyanobacteria, Prochloron didemni, that live symbiotically with tropical ascidians and that produce natural products isolated from the animals. Although several animal species were obtained along a >5500 km transect of the Pacific Ocean, P. didemni strains are >97% identical across much of their genomes, with only a few exceptions concentrated in secondary metabolism. Secondary metabolic gene clusters were sporadically present or absent in identical genomic locations with no consistent pattern of co-occurrence. Discrete mutations were observed, leading to new chemicals that we isolated from animals. Functional cassettes encoding diverse chemicals are exchanged among a single population of symbiotic P. didemni that spans the tropical Pacific, providing the host animals with a varying arsenal of secondary metabolites.     

Molecular adaptation of ammonia monooxygenase during independent pH specialization in Thaumarchaeota

Microbes are abundant in nature and often highly adapted to local conditions. While great progress has been made in understanding the ecological factors driving their distribution in complex environments, the underpinning molecular‐evolutionary mechanisms are rarely dissected. Therefore, we scrutinized the coupling of environmental and molecular adaptation in Thaumarchaeota, an abundant archaeal phylum with a key role in ammonia oxidation. These microbes are adapted to a diverse spectrum of environmental conditions, with pH being a key factor shaping their contemporary distribution and evolutionary diversification. We integrated high‐throughput sequencing data spanning a broad representation of ammonia‐oxidizing terrestrial lineages with codon modelling analyses, testing the hypothesis that ammonia monooxygenase subunit A (AmoA) – a highly conserved membrane protein crucial for ammonia oxidation and classical marker in microbial ecology – underwent adaptation during specialization to extreme pH environments. While purifying selection has been an important factor limiting AmoA evolution, we identified episodic shifts in selective pressure at the base of two phylogenetically distant lineages that independently adapted to acidic conditions and subsequently gained lasting ecological success. This involved nonconvergent selective mechanisms (positive selection vs. selection acting on variants fixed during an episode of relaxed selection) leading to unique sets of amino acid substitutions that remained fixed across the radiation of both acidophilic lineages, highlighting persistent adaptive value in acidic environments. Our data demonstrates distinct trajectories of AmoA evolution despite convergent phenotypic adaptation, suggesting that microbial environmental specialization can be associated with diverse signals of molecular adaptation, even for marker genes employed routinely by microbial ecologists.     

Discovering functional novelty in metagenomes: examples from light-mediated processes

The emerging coverage of diverse habitats by metagenomic shotgun data opens new avenues of discovering functional novelty using computational tools. Here, we apply three different concepts for predicting novel functions within light-mediated microbial pathways in five diverse environments. Using phylogenetic approaches, we discovered two novel deep-branching subfamilies of photolyases (involved in light-mediated repair) distributed abundantly in high-UV environments. Using neighborhood approaches, we were able to assign seven novel functional partners in luciferase synthesis, nitrogen metabolism, and quorum sensing to BLUF domain-containing proteins (involved in light sensing). Finally, by domain analysis, for RcaE proteins (involved in chromatic adaptation), we predict 16 novel domain architectures that indicate novel functionalities in habitats with little or no light. Quantification of protein abundance in the various environments supports our findings that bacteria utilize light for sensing, repair, and adaptation far more widely than previously thought. While the discoveries illustrate the opportunities in function discovery, we also discuss the immense conceptual and practical challenges that come along with this new type of data.