Functional metagenomics reveals novel salt tolerance loci from the human gut microbiome

Metagenomics is a powerful tool that allows for the culture-independent analysis of complex microbial communities. One of the most complex and dense microbial ecosystems known is that of the human distal colon, with cell densities reaching up to 10¹² per gram of faeces. With the majority of species as yet uncultured, there are an enormous number of novel genes awaiting discovery. In the current study, we conducted a functional screen of a metagenomic library of the human gut microbiota for potential salt-tolerant clones. Using transposon mutagenesis, three genes were identified from a single clone exhibiting high levels of identity to a species from the genus Collinsella (closest relative being Collinsella aerofaciens) (COLAER_01955, COLAER_01957 and COLAER_01981), a high G+C, Gram-positive member of the Actinobacteria commonly found in the human gut. The encoded proteins exhibit a strong similarity to GalE, MurB and MazG. Furthermore, pyrosequencing and bioinformatic analysis of two additional fosmid clones revealed the presence of an additional galE and mazG gene, with the highest level of genetic identity to Akkermansia muciniphila and Eggerthella sp. YY7918, respectively. Cloning and heterologous expression of the genes in the osmosensitive strain, Escherichia coli MKH13, resulted in increased salt tolerance of the transformed cells. It is hoped that the identification of atypical salt tolerance genes will help to further elucidate novel salt tolerance mechanisms, and will assist our increased understanding how resident bacteria cope with the osmolarity of the gastrointestinal tract.     

Metagenomic Identification of a Novel Salt Tolerance Gene from the Human Gut Microbiome Which Encodes a Membrane Protein with Homology to a brp/blh-Family β-Carotene 15,15′-Monooxygenase

The human gut microbiome consists of at least 3 million non-redundant genes, 150 times that of the core human genome. Herein, we report the identification and characterisation of a novel stress tolerance gene from the human gut metagenome. The locus, assigned brpA, encodes a membrane protein with homology to a brp/blh-family β-carotene monooxygenase. Cloning and heterologous expression of brpA in Escherichia coli confers a significant salt tolerance phenotype. Furthermore, when cultured in the presence of exogenous β-carotene, cell pellets adopt a red/orange pigmentation indicating the incorporation of carotenoids in the cell membrane.     

Identifying Keystone Species in the Human Gut Microbiome from Metagenomic Timeseries Using Sparse Linear Regression

Human associated microbial communities exert tremendous influence over human health and disease. With modern metagenomic sequencing methods it is now possible to follow the relative abundance of microbes in a community over time. These microbial communities exhibit rich ecological dynamics and an important goal of microbial ecology is to infer the ecological interactions between species directly from sequence data. Any algorithm for inferring ecological interactions must overcome three major obstacles: 1) a correlation between the abundances of two species does not imply that those species are interacting, 2) the sum constraint on the relative abundances obtained from metagenomic studies makes it difficult to infer the parameters in timeseries models, and 3) errors due to experimental uncertainty, or mis-assignment of sequencing reads into operational taxonomic units, bias inferences of species interactions due to a statistical problem called “errors-in-variables”. Here we introduce an approach, Learning Interactions from MIcrobial Time Series (LIMITS), that overcomes these obstacles. LIMITS uses sparse linear regression with boostrap aggregation to infer a discrete-time Lotka-Volterra model for microbial dynamics. We tested LIMITS on synthetic data and showed that it could reliably infer the topology of the inter-species ecological interactions. We then used LIMITS to characterize the species interactions in the gut microbiomes of two individuals and found that the interaction networks varied significantly between individuals. Furthermore, we found that the interaction networks of the two individuals are dominated by distinct “keystone species”, Bacteroides fragilis and Bacteroided stercosis, that have a disproportionate influence on the structure of the gut microbiome even though they are only found in moderate abundance. Based on our results, we hypothesize that the abundances of certain keystone species may be responsible for individuality in the human gut microbiome.     

Comparative Analysis of the Intestinal Bacterial and RNA Viral Communities from Sentinel Birds Placed on Selected Broiler Chicken Farms

There is a great deal of interest in characterizing the complex microbial communities in the poultry gut, and in understanding the effects of these dynamic communities on poultry performance, disease status, animal welfare, and microbes with human health significance. Investigations characterizing the poultry enteric virome have identified novel poultry viruses, but the roles these viruses play in disease and performance problems have yet to be fully characterized. The complex bacterial community present in the poultry gut influences gut development, immune status, and animal health, each of which can be an indicator of overall performance. The present metagenomic investigation was undertaken to provide insight into the colonization of specific pathogen free chickens by enteric microorganisms under field conditions and to compare the pre-contact intestinal microbiome with the altered microbiome following contact with poultry raised in the field. Analysis of the intestinal virome from contact birds (“sentinels”) placed on farms revealed colonization by members of the Picornaviridae, Picobirnaviridae, Reoviridae, and Astroviridae that were not present in pre-contact birds or present in proportionally lower numbers. Analysis of the sentinel gut bacterial community revealed an altered community in the post-contact birds, notably by members of the Lachnospiracea/Clostridium and Lactobacillus families and genera. Members of the avian enteric Reoviridae and Astroviridae have been well-characterized and have historically been implicated in poultry enteric disease; members of the Picobirnaviridae and Picornaviridae have only relatively recently been described in the poultry and avian gut, and their roles in the recognized disease syndromes and in poultry performance in general have not been determined. This metagenomic analysis has provided insight into the colonization of the poultry gut by enteric microbes circulating in commercial broiler flocks, and has identified enteric viruses and virus communities that warrant further study in order to understand their role(s) in avian gut health and disease.     

Gene-targeted metagenomic analysis of glucan-branching enzyme gene profiles among human and animal fecal microbiota

Glycoside hydrolases (GHs), the enzymes that breakdown complex carbohydrates, are a highly diversified class of key enzymes associated with the gut microbiota and its metabolic functions. To learn more about the diversity of GHs and their potential role in a variety of gut microbiomes, we used a combination of 16S, metagenomic and targeted amplicon sequencing data to study one of these enzyme families in detail. Specifically, we employed a functional gene-targeted metagenomic approach to the 1-4-α-glucan-branching enzyme (gBE) gene in the gut microbiomes of four host species (human, chicken, cow and pig). The characteristics of operational taxonomic units (OTUs) and operational glucan-branching units (OGBUs) were distinctive in each of hosts. Human and pig were most similar in OTUs profiles while maintaining distinct OGBU profiles. Interestingly, the phylogenetic profiles identified from 16S and gBE gene sequences differed, suggesting the presence of different gBE genes in the same OTU across different vertebrate hosts. Our data suggest that gene-targeted metagenomic analysis is useful for an in-depth understanding of the diversity of a particular gene of interest. Specific carbohydrate metabolic genes appear to be carried by distinct OTUs in different individual hosts and among different vertebrate species'' microbiomes, the characteristics of which differ according to host genetic background and/or diet.     

Sequence analysis of a salt tolerant metagenomic clone

Metagenome represent an unlimited resource for discovery of novel genes. Here we report, sequence analysis of a salt tolerant metagenomic clone (6B4) from a pond water metagenomic library. Clone 6B4 had an insert of 2254 bp with G+C composition of 64.06%. DNA sequence from 6B4 showed homology to DNA sequences from pro-teobacteria indicating origin of 6B4 metagenomic insert from a yet uncharacterized proteobacteria. Two encoded proteins from clone 6B4 showed match with ATP-depen-dent Clp protease adaptor protein (ClpS) and phasin, while two truncated encoded proteins showed match with poly-3-hydroxybutyrate synthase and permease. Clp complex is known to play a role in stress tolerance. Expression of ClpS from metagenomic clone is proposed to be responsible for salt tolerance of the metagenomic clone 6B4.     

Quality control of microbiota metagenomics by k-mer analysis

Background

The biological and clinical consequences of the tight interactions between host and microbiota are rapidly being unraveled by next generation sequencing technologies and sophisticated bioinformatics, also referred to as microbiota metagenomics. The recent success of metagenomics has created a demand to rapidly apply the technology to large case–control cohort studies and to studies of microbiota from various habitats, including habitats relatively poor in microbes. It is therefore of foremost importance to enable a robust and rapid quality assessment of metagenomic data from samples that challenge present technological limits (sample numbers and size). Here we demonstrate that the distribution of overlapping k-mers of metagenome sequence data predicts sequence quality as defined by gene distribution and efficiency of sequence mapping to a reference gene catalogue.

Results

We used serial dilutions of gut microbiota metagenomic datasets to generate well-defined high to low quality metagenomes. We also analyzed a collection of 52 microbiota-derived metagenomes. We demonstrate that k-mer distributions of metagenomic sequence data identify sequence contaminations, such as sequences derived from “empty” ligation products. Of note, k-mer distributions were also able to predict the frequency of sequences mapping to a reference gene catalogue not only for the well-defined serial dilution datasets, but also for 52 human gut microbiota derived metagenomic datasets.

Conclusions

We propose that k-mer analysis of raw metagenome sequence reads should be implemented as a first quality assessment prior to more extensive bioinformatics analysis, such as sequence filtering and gene mapping. With the rising demand for metagenomic analysis of microbiota it is crucial to provide tools for rapid and efficient decision making. This will eventually lead to a faster turn-around time, improved analytical quality including sample quality metrics and a significant cost reduction. Finally, improved quality assessment will have a major impact on the robustness of biological and clinical conclusions drawn from metagenomic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1406-7) contains supplementary material, which is available to authorized users.     

Frame-shifted APOBEC3A encodes two alternative proapoptotic proteins that target the mitochondrial network

The human APOBEC3A (A3A) cytidine deaminase is a powerful DNA mutator enzyme recognized as a major source of somatic mutations in tumor cell genomes. However, there is a discrepancy between APOBEC3A mRNA levels after interferon stimulation in myeloid cells and A3A detection at the protein level. To understand this difference, we investigated the expression of two novel alternative “A3Alt” proteins encoded in the +1-shifted reading frame of the APOBEC3A gene. A3Alt-L and its shorter isoform A3Alt-S appear to be transmembrane proteins targeted to the mitochondrial compartment that induce membrane depolarization and apoptosis. Thus, the APOBEC3A gene represents a new example wherein a single gene encodes two proapoptotic proteins, A3A cytidine deaminases that target the genome and A3Alt proteins that target mitochondria.     

Toll-like Receptor 11 (TLR11) Prevents Salmonella Penetration into the Murine Peyer Patches

Toll-like receptors (TLRs) are key molecular sensors used by the mammalian innate immune system to detect microorganisms. Although TLR functions in colonic immune homeostasis and tolerance to commensal bacteria have been intensively researched, the precise roles of different TLRs in response to pathogen infection in the gut remain elusive. Peyer patches are the major entrance of Salmonella infection and antigen transportation in intestine. Here, we report that, in contrast to TLR5 as a “carrier of Salmonella,” TLR11 works as a “blocker of Salmonella” to prevent highly invasive Salmonella from penetrating into the murine Peyer patches and spreading systemically. TLR11 plays an important role in mediating TNF-α induction and systemic inflammation in response to Salmonella infection. Remarkably, in mice lacking TLR11, apparent hemorrhages at Peyer patches are induced by highly invasive Salmonella, a phenotype resembling human Salmonella infection. Therefore, our results indicate a potentially important role for TLR11 in preventing murine intestinal infection and modulating antigen transportation in the gut and imply an important role for various TLRs in cooperation with tight control of pathogens penetrating into Peyer patches. The TLR11 knock-out mouse can serve as a good animal model to study Salmonella infection.     

MyTaxa: an advanced taxonomic classifier for genomic and metagenomic sequences

Determining the taxonomic affiliation of sequences assembled from metagenomes remains a major bottleneck that affects research across the fields of environmental, clinical and evolutionary microbiology. Here, we introduce MyTaxa, a homology-based bioinformatics framework to classify metagenomic and genomic sequences with unprecedented accuracy. The distinguishing aspect of MyTaxa is that it employs all genes present in an unknown sequence as classifiers, weighting each gene based on its (predetermined) classifying power at a given taxonomic level and frequency of horizontal gene transfer. MyTaxa also implements a novel classification scheme based on the genome-aggregate average amino acid identity concept to determine the degree of novelty of sequences representing uncharacterized taxa, i.e. whether they represent novel species, genera or phyla. Application of MyTaxa on in silico generated (mock) and real metagenomes of varied read length (100–2000 bp) revealed that it correctly classified at least 5% more sequences than any other tool. The analysis also showed that ∼10% of the assembled sequences from human gut metagenomes represent novel species with no sequenced representatives, several of which were highly abundant in situ such as members of the Prevotella genus. Thus, MyTaxa can find several important applications in microbial identification and diversity studies.     

Screening of Metagenomic and Genomic Libraries Reveals Three Classes of Bacterial Enzymes That Overcome the Toxicity of Acrylate

Acrylate is produced in significant quantities through the microbial cleavage of the highly abundant marine osmoprotectant dimethylsulfoniopropionate, an important process in the marine sulfur cycle. Acrylate can inhibit bacterial growth, likely through its conversion to the highly toxic molecule acrylyl-CoA. Previous work identified an acrylyl-CoA reductase, encoded by the gene acuI, as being important for conferring on bacteria the ability to grow in the presence of acrylate. However, some bacteria lack acuI, and, conversely, many bacteria that may not encounter acrylate in their regular environments do contain this gene. We therefore sought to identify new genes that might confer tolerance to acrylate. To do this, we used functional screening of metagenomic and genomic libraries to identify novel genes that corrected an E. coli mutant that was defective in acuI, and was therefore hyper-sensitive to acrylate. The metagenomic libraries yielded two types of genes that overcame this toxicity. The majority encoded enzymes resembling AcuI, but with significant sequence divergence among each other and previously ratified AcuI enzymes. One other metagenomic gene, arkA, had very close relatives in Bacillus and related bacteria, and is predicted to encode an enoyl-acyl carrier protein reductase, in the same family as FabK, which catalyses the final step in fatty-acid biosynthesis in some pathogenic Firmicute bacteria. A genomic library of Novosphingobium, a metabolically versatile alphaproteobacterium that lacks both acuI and arkA, yielded vutD and vutE, two genes that, together, conferred acrylate resistance. These encode sequential steps in the oxidative catabolism of valine in a pathway in which, significantly, methacrylyl-CoA is a toxic intermediate. These findings expand the range of bacteria for which the acuI gene encodes a functional acrylyl-CoA reductase, and also identify novel enzymes that can similarly function in conferring acrylate resistance, likely, again, through the removal of the toxic product acrylyl-CoA.     

The novel anti-CRISPR AcrIIA22 relieves DNA torsion in target plasmids and impairs SpyCas9 activity

To overcome CRISPR-Cas defense systems, many phages and mobile genetic elements (MGEs) encode CRISPR-Cas inhibitors called anti-CRISPRs (Acrs). Nearly all characterized Acrs directly bind Cas proteins to inactivate CRISPR immunity. Here, using functional metagenomic selection, we describe AcrIIA22, an unconventional Acr found in hypervariable genomic regions of clostridial bacteria and their prophages from human gut microbiomes. AcrIIA22 does not bind strongly to SpyCas9 but nonetheless potently inhibits its activity against plasmids. To gain insight into its mechanism, we obtained an X-ray crystal structure of AcrIIA22, which revealed homology to PC4-like nucleic acid–binding proteins. Based on mutational analyses and functional assays, we deduced that acrIIA22 encodes a DNA nickase that relieves torsional stress in supercoiled plasmids. This may render them less susceptible to SpyCas9, which uses free energy from negative supercoils to form stable R-loops. Modifying DNA topology may provide an additional route to CRISPR-Cas resistance in phages and MGEs.

Derived from phages of the gut microbiome, this study describes a CRISPR-Cas9 inhibitor with an unexpected mechanism; instead of binding Cas9 itself, this “anti-CRISPR” relaxes plasmid DNA and enables Cas9 evasion, suggesting that DNA topology is an underappreciated battleground in phage-bacterial conflicts.     

Selection of Orphan Rhs Toxin Expression in Evolved Salmonella enterica Serovar Typhimurium

Clonally derived bacterial populations exhibit significant genotypic and phenotypic diversity that contribute to fitness in rapidly changing environments. Here, we show that serial passage of Salmonella enterica serovar Typhimurium LT2 (StLT2) in broth, or within a mouse host, results in selection of an evolved population that inhibits the growth of ancestral cells by direct contact. Cells within each evolved population gain the ability to express and deploy a cryptic “orphan” toxin encoded within the rearrangement hotspot (rhs) locus. The Rhs orphan toxin is encoded by a gene fragment located downstream of the “main” rhs gene in the ancestral strain StLT2. The Rhs orphan coding sequence is linked to an immunity gene, which encodes an immunity protein that specifically blocks Rhs orphan toxin activity. Expression of the Rhs orphan immunity protein protects ancestral cells from the evolved lineages, indicating that orphan toxin activity is responsible for the observed growth inhibition. Because the Rhs orphan toxin is encoded by a fragmented reading frame, it lacks translation initiation and protein export signals. We provide evidence that evolved cells undergo recombination between the main rhs gene and the rhs orphan toxin gene fragment, yielding a fusion that enables expression and delivery of the orphan toxin. In this manner, rhs locus rearrangement provides a selective advantage to a subpopulation of cells. These observations suggest that rhs genes play important roles in intra-species competition and bacterial evolution.     

Phylogenetic Characterization of Two Novel Commensal Bacteria Involved with Innate Immune Homeostasis in Drosophila melanogaster

During a previous study on the molecular interaction between commensal bacteria and host gut immunity, two novel bacterial strains, A911^T and G707^T, were isolated from the gut of Drosophila melanogaster. In this study, these strains were characterized in a polyphasic taxonomic study using phenotypic, genetic, and chemotaxonomic analyses. We show that the strains represent novel species in the family Acetobacteraceae. Strain G707^T, a highly pathogenic organism, represents a new species in the genus Gluconobacter, “Gluconobacter morbifer” sp. nov. (type strain G707 = KCTC 22116^T = JCM 15512^T). Strain A911^T, dominantly present in the normal Drosphila gut community, represents a novel genus and species, designated “Commensalibacter intestini” gen. nov., sp. nov. (type strain A911 = KCTC 22117^T = JCM 15511^T).     

Metagenomics,paratransgenesis and the Anopheles microbiome: a portrait of the geographical distribution of the anopheline microbiota based on a meta-analysis of reported taxa

Anophelines harbour a diverse microbial consortium that may represent an extended gene pool for the host. The proposed effects of the insect microbiota span physiological, metabolic and immune processes. Here we synthesise how current metagenomic tools combined with classical culture-dependent techniques provide new insights in the elucidation of the role of the Anopheles-associated microbiota. Many proposed malaria control strategies have been based upon the immunomodulating effects that the bacterial components of the microbiota appear to exert and their ability to express anti-Plasmodium peptides. The number of identified bacterial taxa has increased in the current “omics” era and the available data are mostly scattered or in “tables” that are difficult to exploit. Published microbiota reports for multiple anopheline species were compiled in an Excel^® spreadsheet. We then filtered the microbiota data using a continent-oriented criterion and generated a visual correlation showing the exclusive and shared bacterial genera among four continents. The data suggested the existence of a core group of bacteria associated in a stable manner with their anopheline hosts. However, the lack of data from Neotropical vectors may reduce the possibility of defining the core microbiota and understanding the mosquito-bacteria interactive consortium.     

Enrichment and physiological characterization of a novel comammox Nitrospira indicates ammonium inhibition of complete nitrification

The recent discovery of bacteria within the genus Nitrospira capable of complete ammonia oxidation (comammox) demonstrated that the sequential oxidation of ammonia to nitrate via nitrite can also be performed within a single bacterial cell. Although comammox Nitrospira exhibit a wide distribution in natural and engineered ecosystems, information on their physiological properties is scarce due to the limited number of cultured representatives. Additionally, most available genomic information is derived from metagenomic sequencing and high-quality genomes of Nitrospira in general are limited. In this study, we obtained a high (90%) enrichment of a novel comammox species, tentatively named “Candidatus Nitrospira kreftii”, and performed a detailed genomic and physiological characterization. The complete genome of “Ca. N. kreftii” allowed reconstruction of its basic metabolic traits. Similar to Nitrospira inopinata, the enrichment culture exhibited a very high ammonia affinity (K_{m(app)_NH3} ≈ 0.040 ± 0.01 µM), but a higher nitrite affinity (K_{m(app)_NO2}- = 12.5 ± 4.0 µM), indicating an adaptation to highly oligotrophic environments. Furthermore, we observed partial inhibition of ammonia oxidation at ammonium concentrations as low as 25 µM. This inhibition of “Ca. N. kreftii” indicates that differences in ammonium tolerance rather than affinity could potentially be a niche determining factor for different comammox Nitrospira.Subject terms: Bacterial genomics, Environmental microbiology, Bacterial physiology     

The Motility Symbiont of the Termite Gut Flagellate Caduceia versatilis Is a Member of the “Synergistes” Group

The flagellate Caduceia versatilis in the gut of the termite Cryptotermes cavifrons reportedly propels itself not by its own flagella but solely by the flagella of ectosymbiotic bacteria. Previous microscopic observations have revealed that the motility symbionts are flagellated rods partially embedded in the host cell surface and that, together with a fusiform type of ectosymbiotic bacteria without flagella, they cover almost the entire surface. To identify these ectosymbionts, we conducted 16S rRNA clone analyses of bacteria physically associated with the Caduceia cells. Two phylotypes were found to predominate in the clone library and were phylogenetically affiliated with the “Synergistes” phylum and the order Bacteroidales in the Bacteroidetes phylum. Probes specifically targeting 16S rRNAs of the respective phylotypes were designed, and fluorescence in situ hybridization (FISH) was performed. As a result, the “Synergistes” phylotype was identified as the motility symbiont; the Bacteroidales phylotype was the fusiform ectobiont. The “Synergistes” phylotype was a member of a cluster comprising exclusively uncultured clones from the guts of various termite species. Interestingly, four other phylotypes in this cluster, including the one sharing 95% sequence identity with the motility symbiont, were identified as nonectosymbiotic, or free-living, gut bacteria by FISH. We thus suggest that the motility ectosymbiont has evolved from a free-living gut bacterium within this termite-specific cluster. Based on these molecular and previous morphological data, we here propose a novel genus and species, “Candidatus Tammella caduceiae,” for this unique motility ectosymbiont of Caducaia versatilis.