To assemble or not to assemble: metagenomic profiling of microbially mediated biogeochemical pathways in complex communities

High-throughput profiling of microbial functional traits involved in various biogeochemical cycling pathways using shotgun metagenomic sequencing has been routinely applied in microbial ecology and environmental science. Multiple bioinformatics data processing approaches are available, including assembly-based (single-sample assembly and multi-sample assembly) and read-based (merged reads and raw data). However, it remains not clear how these different approaches may differ in data analyses and affect result interpretation. In this study, using two typical shotgun metagenome datasets recovered from geographically distant coastal sediments, the performance of different data processing approaches was comparatively investigated from both technical and biological/ecological perspectives. Microbially mediated biogeochemical cycling pathways, including nitrogen cycling, sulfur cycling and B12 biosynthesis, were analyzed. As a result, multi-sample assembly provided the most amount of usable information for targeted functional traits, at a high cost of computational resources and running time. Single-sample assembly and read-based analysis were comparable in obtaining usable information, but the former was much more time- and resource-consuming. Critically, different approaches introduced much stronger variations in microbial profiles than biological differences. However, community-level differences between the two sampling sites could be consistently observed despite the approaches being used. In choosing an appropriate approach, researchers shall balance the trade-offs between multiple factors, including the scientific question, the amount of usable information, computational resources and time cost. This study is expected to provide valuable technical insights and guidelines for the various approaches used for metagenomic data analysis.

Metagenomic profiling of placental tissue suggests DNA virus infection of the placenta is rare

It is widely recognized that pathogens can be transmitted across the placenta from mother to foetus. Recent re-evaluation of metagenomic studies indicates that the placenta has no unique microbiome of commensal bacteria. However, viral transmission across the placenta, including transmission of DNA viruses such as the human herpesviruses, is possible. A fuller understanding of which DNA virus sequence can be found in the placenta is required. We employed a metagenomic analysis to identify viral DNA sequences in placental metagenomes from full-term births (20 births), pre-term births (13 births), births from pregnancies associated with antenatal infections (12 births) or pre-term births with antenatal infections (three births). Our analysis found only a small number of DNA sequences corresponding to the genomes of human herpesviruses in four of the 48 metagenomes analysed. Therefore, our data suggest that DNA virus infection of the placenta is rare and support the concept that the placenta is largely free of pathogen infection.

Evaluation of the Microba Community Profiler for Taxonomic Profiling of Metagenomic Datasets From the Human Gut Microbiome

A fundamental goal of microbial ecology is to accurately determine the species composition in a given microbial ecosystem. In the context of the human microbiome, this is important for establishing links between microbial species and disease states. Here we benchmark the Microba Community Profiler (MCP) against other metagenomic classifiers using 140 moderate to complex in silico microbial communities and a standardized reference genome database. MCP generated accurate relative abundance estimates and made substantially fewer false positive predictions than other classifiers while retaining a high recall rate. We further demonstrated that the accuracy of species classification was substantially increased using the Microba Genome Database, which is more comprehensive than reference datasets used by other classifiers and illustrates the importance of including genomes of uncultured taxa in reference databases. Consequently, MCP classifies appreciably more reads than other classifiers when using their recommended reference databases. These results establish MCP as best-in-class with the ability to produce comprehensive and accurate species profiles of human gastrointestinal samples.

The Neuromodulator-Encoding sadA Gene Is Widely Distributed in the Human Skin Microbiome

Trace amines (TA) are endogenously produced in mammals, have a low concentration in the central nervous system (CNS), but trigger a variety of neurological effects and intervene in host cell communication. It emerged that neurotransmitters and TA are produced also by the microbiota. As it has been shown that TA contribute to wound healing, we examined the skin microbiome of probands using shotgun metagenomics. The phyla Actinobacteria, Proteobacteria, Firmicutes, and Bacteroidetes were predominant. Since SadA is a highly promiscuous TA-producing decarboxylase in Firmicutes, the skin microbiome was specifically examined for the presence of sadA-homologous genes. By mapping the reads of certain genes, we found that, although there were less reads mapping to sadA than to ubiquitous housekeeping genes (arcC and mutS), normalized reads counts were still >1000 times higher than those of rare control genes (icaA, icaB, and epiA). At protein sequence level SadA homologs were found in at least 7 phyla: Firmicutes, Actinobacteria, Proteobacteria, Bacteroidetes, Acidobacteria, Chloroflexi, and Cyanobacteria, and in 23 genera of the phylum Firmicutes. A high proportion of the genera that have a SadA homolog belong to the classical skin and intestinal microbiota. The distribution of sadA in so many different phyla illustrates the importance of horizontal gene transfer (HGT). We show that the sadA gene is widely distributed in the human skin microbiome. When comparing the sadA read counts in the probands, there was no correlation between age and gender, but an enormous difference in the sadA read counts in the microbiome of the individuals. Since sadA is involved in TA synthesis, it is likely that the TA content of the skin is correlated with the amount of TA producing bacteria in the microbiome. In this way, the microbiome-generated TA could influence signal transmission in the epithelial and nervous system.

RBUD: A New Functional Potential Analysis Approach for Whole Microbial Genome Shotgun Sequencing

Whole metagenome shotgun sequencing is a powerful approach to detect the functional potential of microbial communities. Currently, the read-based metagenomics profiling for established database (RBED) method is one of the two kinds of conventional methods for species and functional annotations. However, the databases, which are established based on test samples or specific reference genomes or protein sequences, limit the coverage of global microbial diversity. The other assembly-based metagenomics profiling for unestablished database (ABUD) method has a low utilization rate of reads, resulting in a lot of biological information loss. In this study, we proposed a new method, read-based metagenomics profiling for unestablished database (RBUD), based on Metagenome Database of Global Microorganisms (MDGM), to solve the above problems. To evaluate the accuracy and effectiveness of our method, the intestinal bacterial composition and function analyses were performed in both avian colibacillosis chicken cases and type 2 diabetes mellitus patients. Comparing to the existing methods, RBUD is superior in detecting proteins, percentage of reads mapping and ontological similarity of intestinal microbes. The results of RBUD are in better agreement with the classical functional studies on these two diseases. RBUD also has the advantages of fast analysis speed and is not limited by the sample size.