Accurate Profiling of Microbial Communities from Massively Parallel Sequencing using Convex Optimization

We describe the Microbial Community Reconstruction ({\bf MCR}) Problem, which is fundamental for microbiome analysis. In this problem, the goal is to reconstruct the identity and frequency of species comprising a microbial community, using short sequence reads from Massively Parallel Sequencing (MPS) data obtained for specified genomic regions. We formulate the problem mathematically as a convex optimization problem and provide sufficient conditions for identifiability, namely the ability to reconstruct species identity and frequency correctly when the data size (number of reads) grows to infinity. We discuss different metrics for assessing the quality of the reconstructed solution, including a novel phylogenetically-aware metric based on the Mahalanobis distance, and give upper-bounds on the reconstruction error for a finite number of reads under different metrics. We propose a scalable divide-and-conquer algorithm for the problem using convex optimization, which enables us to handle large problems (with $\sim10^6$ species). We show using numerical simulations that for realistic scenarios, where the microbial communities are sparse, our algorithm gives solutions with high accuracy, both in terms of obtaining accurate frequency, and in terms of species phylogenetic resolution.Comment: To appear in SPIRE 1

Species-level functional profiling of metagenomes and metatranscriptomes.

Functional profiles of microbial communities are typically generated using comprehensive metagenomic or metatranscriptomic sequence read searches, which are time-consuming, prone to spurious mapping, and often limited to community-level quantification. We developed HUMAnN2, a tiered search strategy that enables fast, accurate, and species-resolved functional profiling of host-associated and environmental communities. HUMAnN2 identifies a community's known species, aligns reads to their pangenomes, performs translated search on unclassified reads, and finally quantifies gene families and pathways. Relative to pure translated search, HUMAnN2 is faster and produces more accurate gene family profiles. We applied HUMAnN2 to study clinal variation in marine metabolism, ecological contribution patterns among human microbiome pathways, variation in species' genomic versus transcriptional contributions, and strain profiling. Further, we introduce 'contributional diversity' to explain patterns of ecological assembly across different microbial community types

Transkingdom Networks: A Systems Biology Approach to Identify Causal Members of Host-Microbiota Interactions

Improvements in sequencing technologies and reduced experimental costs have resulted in a vast number of studies generating high-throughput data. Although the number of methods to analyze these "omics" data has also increased, computational complexity and lack of documentation hinder researchers from analyzing their high-throughput data to its true potential. In this chapter we detail our data-driven, transkingdom network (TransNet) analysis protocol to integrate and interrogate multi-omics data. This systems biology approach has allowed us to successfully identify important causal relationships between different taxonomic kingdoms (e.g. mammals and microbes) using diverse types of data

Metagenomic analysis of the saliva microbiome with merlin

In recent years, metagenomics has demonstrated to play an essential role on the study of the microorganisms that live in microbial communities, particularly those who inhabit the human body. Several bioinformatics tools and pipelines have been developed for the analysis of these data, but they usually only address one topic: to identify the taxonomic composition or to address the metabolic functional profile. This work aimed to implement a computational framework able to answer the two questions simultaneously. Merlin, a previously released software aiming at the reconstruction of genome-scale metabolic models for single organisms, was extended to deal with metagenomics data. It has an user-friendly and intuitive interface, being suitable for those with limited bioinformatics skills. The performance of the tool was evaluated with samples from the Human Microbiome Project, particularly from saliva. Overall, the results show the same patterns reported before: while the pathways needed for microbial life remain relatively stable, the community composition varies extensively among individuals

Patent Human Infections with the Whipworm, Trichuris trichiura, Are Not Associated with Alterations in the Faecal Microbiota

Background: The soil-transmitted helminth (STH), Trichuris trichiura colonises the human large intestine where it may modify inflammatory responses, an effect possibly mediated through alterations in the intestinal microbiota. We hypothesised that patent T. trichiura infections would be associated with altered faecal microbiota and that anthelmintic treatment would induce a microbiota resembling more closely that observed in uninfected individuals. Materials and Methods: School children in Ecuador were screened for STH infections and allocated to 3 groups: uninfected, T. trichiura only, and mixed infections with T. trichiura and Ascaris lumbricoides. A sample of uninfected children and those with T. trichiura infections only were given anthelmintic treatment. Bacterial community profiles in faecal samples were studied by 454 pyrosequencing of 16 S rRNA genes. Results: Microbiota analyses of faeces were done for 97 children: 30 were uninfected, 17 were infected with T. trichiura, and 50 with T. trichiura and A. lumbricoides. Post-treatment samples were analyzed for 14 children initially infected with T. trichiura alone and for 21 uninfected children. Treatment resulted in 100% cure of STH infections. Comparisons of the microbiota at different taxonomic levels showed no statistically significant differences in composition between uninfected children and those with T. trichiura infections. We observed a decreased proportional abundance of a few bacterial genera from the Clostridia class of Firmicutes and a reduced bacterial diversity among children with mixed infections compared to the other two groups, indicating a possible specific effect of A. lumbricoides infection. Anthelmintic treatment of children with T. trichiura did not alter faecal microbiota composition. Discussion: Our data indicate that patent human infections with T. trichiura may have no effect on faecal microbiota but that A. lumbricoides colonisation might be associated with a disturbed microbiota. Our results also catalogue the microbiota of rural Ecuadorians and indicate differences with individuals from more urban industrialised societies

Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences

Profiling phylogenetic marker genes, such as the 16S rRNA gene, is a key tool for studies of microbial communities but does not provide direct evidence of a community’s functional capabilities. Here we describe PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States), a computational approach to predict the functional composition of a metagenome using marker gene data and a database of reference genomes. PICRUSt uses an extended ancestral-state reconstruction algorithm to predict which gene families are present and then combines gene families to estimate the composite metagenome. Using 16S information, PICRUSt recaptures key findings from the Human Microbiome Project and accurately predicts the abundance of gene families in host-associated and environmental communities, with quantifiable uncertainty. Our results demonstrate that phylogeny and function are sufficiently linked that this ‘predictive metagenomic’ approach should provide useful insights into the thousands of uncultivated microbial communities for which only marker gene surveys are currently available

Genomic Expansion of Magnetotactic Bacteria Reveals an Early Common Origin of Magnetotaxis with Lineage-specific Evolution

The origin and evolution of magnetoreception, which in diverse prokaryotes and protozoa is known as magnetotaxis and enables these microorganisms to detect Earth’s magnetic field for orientation and navigation, is not well understood in evolutionary biology. The only known prokaryotes capable of sensing the geomagnetic field are magnetotactic bacteria (MTB), motile microorganisms that biomineralize intracellular, membrane-bounded magnetic single-domain crystals of either magnetite (Fe3O4) or greigite (Fe3S4) called magnetosomes. Magnetosomes are responsible for magnetotaxis in MTB. Here we report the first large-scale metagenomic survey of MTB from both northern and southern hemispheres combined with 28 genomes from uncultivated MTB. These genomes expand greatly the coverage of MTB in the Proteobacteria, Nitrospirae, and Omnitrophica phyla, and provide the first genomic evidence of MTB belonging to the Zetaproteobacteria and “Candidatus Lambdaproteobacteria” classes. The gene content and organization of magnetosome gene clusters, which are physically grouped genes that encode proteins for magnetosome biosynthesis and organization, are more conserved within phylogenetically similar groups than between different taxonomic lineages. Moreover, the phylogenies of core magnetosome proteins form monophyletic clades. Together, these results suggest a common ancient origin of iron-based (Fe3O4 and Fe3S4) magnetotaxis in the domain Bacteria that underwent lineage-specific evolution, shedding new light on the origin and evolution of biomineralization and magnetotaxis, and expanding significantly the phylogenomic representation of MTB

Microbial differences between dental plaque and historic dental calculus are related to oral biofilm maturation stage

Dental calculus, calcified oral plaque biofilm, contains microbial and host biomolecules that can be used to study historic microbiome communities and host responses. Dental calculus does not typically accumulate as much today as historically, and clinical oral microbiome research studies focus primarily on living dental plaque biofilm. However, plaque and calculus reflect different conditions of the oral biofilm, and the differences in microbial characteristics between the sample types have not yet been systematically explored. Here, we compare the microbial profiles of modern dental plaque, modern dental calculus, and historic dental calculus to establish expected differences between these substrates.- Background - Results -- Authentication of a preserved oral biofilm in calculus samples -- Dental calculus and plaque biofilm communities are distinct -- Health-associated communities of dental plaque and calculus are distinct -- Signatures of health and of disease are shared in modern and historic calculus samples -- Microbial community differences between health and disease in calculus are poorly resolved -- Absence of caries-specific microbial profiles in dental calculus -- Microbial co-exclusion patterns in plaque and calculus reflect biofilm maturity -- Microbial complexes in plaque and calculus -- Functional prediction in calculus is poorly predictive of health status -- Proteomic profiles of historic healthy site calculus -- Correlations between taxonomic, proteomic, and metabolomic profiles - Discussion - Conclusions - Materials and methods --Historic and modern calculus sample collection DNA extraction -- DNA library construction and high-throughput sequencing -- DNA sequence processing -- Genetic assessment of historic calculus sample preservation -- Genetic microbial taxonomic profiling -- Principal component analysis -- Assessment of differentially abundant taxa -- Sparse partial least squares-discriminant analysis -- Assessment of microbial co-exclusion patterns -- Gene functional categorization with SEED -- Proteomics -- Metabolomics -- Regularized canonical correlation analysi

The complete mitochondrial genome of the foodborne parasitic pathogen Cyclospora cayetanensis

Cyclospora cayetanensis is a human-specific coccidian parasite responsible for several food and water-related outbreaks around the world, including the most recent ones involving over 900 persons in 2013 and 2014 outbreaks in the USA. Multicopy organellar DNA such as mitochondrion genomes have been particularly informative for detection and genetic traceback analysis in other parasites. We sequenced the C. cayetanensis genomic DNA obtained from stool samples from patients infected with Cyclospora in Nepal using the Illumina MiSeq platform. By bioinformatically filtering out the metagenomic reads of non-coccidian origin sequences and concentrating the reads by targeted alignment, we were able to obtain contigs containing Eimeria-like mitochondrial, apicoplastic and some chromosomal genomic fragments. A mitochondrial genomic sequence was assembled and confirmed by cloning and sequencing targeted PCR products amplified from Cyclospora DNA using primers based on our draft assembly sequence. The results show that the C. cayetanensis mitochondrion genome is 6274 bp in length, with 33% GC content, and likely exists in concatemeric arrays as in Eimeria mitochondrial genomes. Phylogenetic analysis of the C. cayetanensis mitochondrial genome places this organism in a tight cluster with Eimeria species. The mitochondrial genome of C. cayetanensis contains three protein coding genes, cytochrome (cytb), cytochrome C oxidase subunit 1 (cox1), and cytochrome C oxidase subunit 3 (cox3), in addition to 14 large subunit (LSU) and nine small subunit (SSU) fragmented rRNA genes

Draft genome sequence of the cadmium-resistant strain JJU2, belonging to the family Hapalosiphonaceae of the cyanobacteria

Here, we report the genome of strain JJU2, a cyanobacterium of the family Hapalosiphonaceae known to be resistant to high cadmium levels, assembled from a nonaxenic, unialgal culture from Marinduque, Philippines. The draft genome is 7.1 Mb long with a GC content of 40.05% and contains 5,625 protein-coding gene