Functional gene arrays-based analysis of fecal microbiomes in patients with liver cirrhosis

Background

Human gut microbiota plays an important role in the pathogenesis of cirrhosis complications. Although the phylogenetic diversity of intestinal microbiota in patients with liver cirrhosis has been examined in several studies, little is known about their functional composition and structure.

Results

To characterize the functional gene diversity of the gut microbiome in cirrhotic patients, we recruited a total of 42 individuals, 12 alcoholic cirrhosis patients, 18 hepatitis B virus (HBV)-related cirrhosis patients, and 12 normal controls. We determined the functional structure of these samples using a specific functional gene array, which is a combination of GeoChip for monitoring biogeochemical processes and HuMiChip specifically designed for analyzing human microbiomes. Our experimental data showed that the microbial community functional composition and structure were dramatically distinctive in the alcoholic cirrhosis. Various microbial functional genes involved in organic remediation, stress response, antibiotic resistance, metal resistance, and virulence were highly enriched in the alcoholic cirrhosis group compared to the control group and HBV-related cirrhosis group. Cirrhosis may have distinct influences on metabolic potential of fecal microbial communities. The abundance of functional genes relevant to nutrient metabolism, including amino acid metabolism, lipid metabolism, nucleotide metabolism, and isoprenoid biosynthesis, were significantly decreased in both alcoholic cirrhosis group and HBV-related cirrhosis group. Significant correlations were observed between functional gene abundances and Child-Pugh scores, such as those encoding aspartate-ammonia ligase, transaldolase, adenylosuccinate synthetase and IMP dehydrogenase.

Conclusions

Functional gene array was utilized to study the gut microbiome in alcoholic and HBV-related cirrhosis patients and controls in this study. Our array data indicated that the functional composition of fecal microbiomes was heavily influenced by cirrhosis, especially by alcoholic cirrhosis. This study provides new insights into the functional potentials and activity of gut microbiota in cirrhotic patients with different etiologies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-753) contains supplementary material, which is available to authorized users.     

Small RNAs from plants,bacteria and fungi within the order Hypocreales are ubiquitous in human plasma

Background

The human microbiome plays a significant role in maintaining normal physiology. Changes in its composition have been associated with bowel disease, metabolic disorders and atherosclerosis. Sequences of microbial origin have been observed within small RNA sequencing data obtained from blood samples. The aim of this study was to characterise the microbiome from which these sequences are derived.

Results

Abundant non-human small RNA sequences were identified in plasma and plasma exosomal samples. Assembly of these short sequences into longer contigs was the pivotal novel step in ascertaining their origin by BLAST searches. Most reads mapped to rRNA sequences. The taxonomic profiles of the microbes detected were very consistent between individuals but distinct from microbiomes reported at other sites. The majority of bacterial reads were from the phylum Proteobacteria, whilst for 5 of 6 individuals over 90% of the more abundant fungal reads were from the phylum Ascomycota; of these over 90% were from the order Hypocreales. Many contigs were from plants, presumably of dietary origin. In addition, extremely abundant small RNAs derived from human Y RNAs were detected.

Conclusions

A characteristic profile of a subset of the human microbiome can be obtained by sequencing small RNAs present in the blood. The source and functions of these molecules remain to be determined, but the specific profiles are likely to reflect health status. The potential to provide biomarkers of diet and for the diagnosis and prognosis of human disease is immense.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-933) contains supplementary material, which is available to authorized users.     

Metagenomic analysis of the Rhinopithecus bieti fecal microbiome reveals a broad diversity of bacterial and glycoside hydrolase profiles related to lignocellulose degradation

Background

The animal gastrointestinal tract contains a complex community of microbes, whose composition ultimately reflects the co-evolution of microorganisms with their animal host and the diet adopted by the host. Although the importance of gut microbiota of humans has been well demonstrated, there is a paucity of research regarding non-human primates (NHPs), especially herbivorous NHPs.

Results

In this study, an analysis of 97,942 pyrosequencing reads generated from Rhinopithecus bieti fecal DNA extracts was performed to help better understanding of the microbial diversity and functional capacity of the R. bieti gut microbiome. The taxonomic analysis of the metagenomic reads indicated that R. bieti fecal microbiomes were dominated by Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria phyla. The comparative analysis of taxonomic classification revealed that the metagenome of R. bieti was characterized by an overrepresentation of bacteria of phylum Fibrobacteres and Spirochaetes as compared with other animals. Primary functional categories were associated mainly with protein, carbohydrates, amino acids, DNA and RNA metabolism, cofactors, cell wall and capsule and membrane transport. Comparing glycoside hydrolase profiles of R. bieti with those of other animal revealed that the R. bieti microbiome was most closely related to cow rumen.

Conclusions

Metagenomic and functional analysis demonstrated that R. bieti possesses a broad diversity of bacteria and numerous glycoside hydrolases responsible for lignocellulosic biomass degradation which might reflect the adaptations associated with a diet rich in fibrous matter. These results would contribute to the limited body of NHPs metagenome studies and provide a unique genetic resource of plant cell wall degrading microbial enzymes. However, future studies on the metagenome sequencing of R. bieti regarding the effects of age, genetics, diet and environment on the composition and activity of the metagenomes are required.

Electronic supplementary material

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

Comparative metagenomics reveals host specific metavirulomes and horizontal gene transfer elements in the chicken cecum microbiome

Background

The complex microbiome of the ceca of chickens plays an important role in nutrient utilization, growth and well-being of these animals. Since we have a very limited understanding of the capabilities of most species present in the cecum, we investigated the role of the microbiome by comparative analyses of both the microbial community structure and functional gene content using random sample pyrosequencing. The overall goal of this study was to characterize the chicken cecal microbiome using a pathogen-free chicken and one that had been challenged with Campylobacter jejuni.

Methodology/Principal Findings

Comparative metagenomic pyrosequencing was used to generate 55,364,266 bases of random sampled pyrosequence data from two chicken cecal samples. SSU rDNA gene tags and environmental gene tags (EGTs) were identified using SEED subsystems-based annotations. The distribution of phylotypes and EGTs detected within each cecal sample were primarily from the Firmicutes, Bacteroidetes and Proteobacteria, consistent with previous SSU rDNA libraries of the chicken cecum. Carbohydrate metabolism and virulence genes are major components of the EGT content of both of these microbiomes. A comparison of the twelve major pathways in the SEED Virulence Subsystem (metavirulome) represented in the chicken cecum, mouse cecum and human fecal microbiomes showed that the metavirulomes differed between these microbiomes and the metavirulomes clustered by host environment. The chicken cecum microbiomes had the broadest range of EGTs within the SEED Conjugative Transposon Subsystem, however the mouse cecum microbiomes showed a greater abundance of EGTs in this subsystem. Gene assemblies (32 contigs) from one microbiome sample were predominately from the Bacteroidetes, and seven of these showed sequence similarity to transposases, whereas the remaining sequences were most similar to those from catabolic gene families.

Conclusion/Significance

This analysis has demonstrated that mobile DNA elements are a major functional component of cecal microbiomes, thus contributing to horizontal gene transfer and functional microbiome evolution. Moreover, the metavirulomes of these microbiomes appear to associate by host environment. These data have implications for defining core and variable microbiome content in a host species. Furthermore, this suggests that the evolution of host specific metavirulomes is a contributing factor in disease resistance to zoonotic pathogens.     

Characterization of the fecal microbiome from non-human wild primates reveals species specific microbial communities

Background

Host-associated microbes comprise an integral part of animal digestive systems and these interactions have a long evolutionary history. It has been hypothesized that the gastrointestinal microbiome of humans and other non-human primates may have played significant roles in host evolution by facilitating a range of dietary adaptations. We have undertaken a comparative sequencing survey of the gastrointestinal microbiomes of several non-human primate species, with the goal of better understanding how these microbiomes relate to the evolution of non-human primate diversity. Here we present a comparative analysis of gastrointestinal microbial communities from three different species of Old World wild monkeys.

Methodology/Principal Findings

We analyzed fecal samples from three different wild non-human primate species (black-and-white colobus [Colubus guereza], red colobus [Piliocolobus tephrosceles], and red-tailed guenon [Cercopithecus ascanius]). Three samples from each species were subjected to small subunit rRNA tag pyrosequencing. Firmicutes comprised the vast majority of the phyla in each sample. Other phyla represented were Bacterioidetes, Proteobacteria, Spirochaetes, Actinobacteria, Verrucomicrobia, Lentisphaerae, Tenericutes, Planctomycetes, Fibrobacateres, and TM7. Bray-Curtis similarity analysis of these microbiomes indicated that microbial community composition within the same primate species are more similar to each other than to those of different primate species. Comparison of fecal microbiota from non-human primates with microbiota of human stool samples obtained in previous studies revealed that the gut microbiota of these primates are distinct and reflect host phylogeny.

Conclusion/Significance

Our analysis provides evidence that the fecal microbiomes of wild primates co-vary with their hosts, and that this is manifested in higher intraspecies similarity among wild primate species, perhaps reflecting species specificity of the microbiome in addition to dietary influences. These results contribute to the limited body of primate microbiome studies and provide a framework for comparative microbiome analysis between human and non-human primates as well as a comparative evolutionary understanding of the human microbiome.     

The cervical microbiome over 7 years and a comparison of methodologies for its characterization

Background

The rapidly expanding field of microbiome studies offers investigators a large choice of methods for each step in the process of determining the microorganisms in a sample. The human cervicovaginal microbiome affects female reproductive health, susceptibility to and natural history of many sexually transmitted infections, including human papillomavirus (HPV). At present, long-term behavior of the cervical microbiome in early sexual life is poorly understood.

Methods

The V6 and V6–V9 regions of the 16S ribosomal RNA gene were amplified from DNA isolated from exfoliated cervical cells. Specimens from 10 women participating in the Natural History Study of HPV in Guanacaste, Costa Rica were sampled successively over a period of 5–7 years. We sequenced amplicons using 3 different platforms (Sanger, Roche 454, and Illumina HiSeq 2000) and analyzed sequences using pipelines based on 3 different classification algorithms (usearch, RDP Classifier, and pplacer).

Results

Usearch and pplacer provided consistent microbiome classifications for all sequencing methods, whereas RDP Classifier deviated significantly when characterizing Illumina reads. Comparing across sequencing platforms indicated 7%–41% of the reads were reclassified, while comparing across software pipelines reclassified up to 32% of the reads. Variability in classification was shown not to be due to a difference in read lengths. Six cervical microbiome community types were observed and are characterized by a predominance of either G. vaginalis or Lactobacillus spp. Over the 5–7 year period, subjects displayed fluctuation between community types. A PERMANOVA analysis on pairwise Kantorovich-Rubinstein distances between the microbiota of all samples yielded an F-test ratio of 2.86 (p<0.01), indicating a significant difference comparing within and between subjects’ microbiota.

Conclusions

Amplification and sequencing methods affected the characterization of the microbiome more than classification algorithms. Pplacer and usearch performed consistently with all sequencing methods. The analyses identified 6 community types consistent with those previously reported. The long-term behavior of the cervical microbiome indicated that fluctuations were subject dependent.     

Characteristic male urine microbiomes associate with asymptomatic sexually transmitted infection

Background

The microbiome of the male urogenital tract is poorly described but it has been suggested that bacterial colonization of the male urethra might impact risk of sexually transmitted infection (STI). Previous cultivation-dependent studies showed that a variety of non-pathogenic bacteria colonize the urethra but did not thoroughly characterize these microbiomes or establish links between the compositions of urethral microbiomes and STI.

Methodology/Findings

Here, we used 16S rRNA PCR and sequencing to identify bacteria in urine specimens collected from men who lacked symptoms of urethral inflammation but who differed in status for STI. All of the urine samples contained multiple bacterial genera and many contained taxa that colonize the human vagina. Uncultivated bacteria associated with female genital tract pathology were abundant in specimens from men who had STI.

Conclusions

Urine microbiomes from men with STI were dominated by fastidious, anaerobic and uncultivated bacteria. The same taxa were rare in STI negative individuals. Our findings suggest that the composition of male urine microbiomes is related to STI.     

Characterization of bacteria in biopsies of colon and stools by high throughput sequencing of the V2 region of bacterial 16S rRNA gene in human

Background

The characterization of the human intestinal microflora and their interactions with the host have been identified as key components in the study of intestinal disorders such as inflammatory bowel diseases. High-throughput sequencing has enabled culture-independent studies to deeply analyze bacteria in the gut. It is possible with this technology to systematically analyze links between microbes and the genetic constitution of the host, such as DNA polymorphisms and methylation, and gene expression.

Methods and Findings

In this study the V2 region of the bacterial 16S ribosomal RNA (rRNA) gene using 454 pyrosequencing from seven anatomic regions of human colon and two types of stool specimens were analyzed. The study examined the number of reads needed to ascertain differences between samples, the effect of DNA extraction procedures and PCR reproducibility, and differences between biopsies and stools in order to design a large scale systematic analysis of gut microbes. It was shown (1) that sequence coverage lower than 1,000 reads influenced quantitative and qualitative differences between samples measured by UniFrac distances. Distances between samples became stable after 1,000 reads. (2) Difference of extracted bacteria was observed between the two DNA extraction methods. In particular, Firmicutes Bacilli were not extracted well by one method. (3) Quantitative and qualitative difference in bacteria from ileum to rectum colon were not observed, but there was a significant positive trend between distances within colon and quantitative differences. Between sample type, biopsies or stools, quantitative and qualitative differences were observed.

Conclusions

Results of human colonic bacteria analyzed using high-throughput sequencing were highly dependent on the experimental design, especially the number of sequence reads, DNA extraction method, and sample type.     

Inflammation-associated enterotypes,host genotype,cage and inter-individual effects drive gut microbiota variation in common laboratory mice

Background

Murine models are a crucial component of gut microbiome research. Unfortunately, a multitude of genetic backgrounds and experimental setups, together with inter-individual variation, complicates cross-study comparisons and a global understanding of the mouse microbiota landscape. Here, we investigate the variability of the healthy mouse microbiota of five common lab mouse strains using 16S rDNA pyrosequencing.

Results

We find initial evidence for richness-driven, strain-independent murine enterotypes that show a striking resemblance to those in human, and which associate with calprotectin levels, a marker for intestinal inflammation. After enterotype stratification, we find that genetic, caging and inter-individual variation contribute on average 19%, 31.7% and 45.5%, respectively, to the variance in the murine gut microbiota composition. Genetic distance correlates positively to microbiota distance, so that genetically similar strains have more similar microbiota than genetically distant ones. Specific mouse strains are enriched for specific operational taxonomic units and taxonomic groups, while the ''cage effect'' can occur across mouse strain boundaries and is mainly driven by Helicobacter infections.

Conclusions

The detection of enterotypes suggests a common ecological cause, possibly low-grade inflammation that might drive differences among gut microbiota composition in mammals. Furthermore, the observed environmental and genetic effects have important consequences for experimental design in mouse microbiome research.     

MaxSSmap: a GPU program for mapping divergent short reads to genomes with the maximum scoring subsequence

Background

Programs based on hash tables and Burrows-Wheeler are very fast for mapping short reads to genomes but have low accuracy in the presence of mismatches and gaps. Such reads can be aligned accurately with the Smith-Waterman algorithm but it can take hours and days to map millions of reads even for bacteria genomes.

Results

We introduce a GPU program called MaxSSmap with the aim of achieving comparable accuracy to Smith-Waterman but with faster runtimes. Similar to most programs MaxSSmap identifies a local region of the genome followed by exact alignment. Instead of using hash tables or Burrows-Wheeler in the first part, MaxSSmap calculates maximum scoring subsequence score between the read and disjoint fragments of the genome in parallel on a GPU and selects the highest scoring fragment for exact alignment. We evaluate MaxSSmap’s accuracy and runtime when mapping simulated Illumina E.coli and human chromosome one reads of different lengths and 10% to 30% mismatches with gaps to the E.coli genome and human chromosome one. We also demonstrate applications on real data by mapping ancient horse DNA reads to modern genomes and unmapped paired reads from NA12878 in 1000 genomes.

Conclusions

We show that MaxSSmap attains comparable high accuracy and low error to fast Smith-Waterman programs yet has much lower runtimes. We show that MaxSSmap can map reads rejected by BWA and NextGenMap with high accuracy and low error much faster than if Smith-Waterman were used. On short read lengths of 36 and 51 both MaxSSmap and Smith-Waterman have lower accuracy compared to at higher lengths. On real data MaxSSmap produces many alignments with high score and mapping quality that are not given by NextGenMap and BWA. The MaxSSmap source code in CUDA and OpenCL is freely available from http://www.cs.njit.edu/usman/MaxSSmap.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-969) contains supplementary material, which is available to authorized users.     

CRISPR-Cas systems target a diverse collection of invasive mobile genetic elements in human microbiomes

Background

Bacteria and archaea develop immunity against invading genomes by incorporating pieces of the invaders'' sequences, called spacers, into a clustered regularly interspaced short palindromic repeats (CRISPR) locus between repeats, forming arrays of repeat-spacer units. When spacers are expressed, they direct CRISPR-associated (Cas) proteins to silence complementary invading DNA. In order to characterize the invaders of human microbiomes, we use spacers from CRISPR arrays that we had previously assembled from shotgun metagenomic datasets, and identify contigs that contain these spacers'' targets.

Results

We discover 95,000 contigs that are putative invasive mobile genetic elements, some targeted by hundreds of CRISPR spacers. We find that oral sites in healthy human populations have a much greater variety of mobile genetic elements than stool samples. Mobile genetic elements carry genes encoding diverse functions: only 7% of the mobile genetic elements are similar to known phages or plasmids, although a much greater proportion contain phage- or plasmid-related genes. A small number of contigs share similarity with known integrative and conjugative elements, providing the first examples of CRISPR defenses against this class of element. We provide detailed analyses of a few large mobile genetic elements of various types, and a relative abundance analysis of mobile genetic elements and putative hosts, exploring the dynamic activities of mobile genetic elements in human microbiomes. A joint analysis of mobile genetic elements and CRISPRs shows that protospacer-adjacent motifs drive their interaction network; however, some CRISPR-Cas systems target mobile genetic elements lacking motifs.

Conclusions

We identify a large collection of invasive mobile genetic elements in human microbiomes, an important resource for further study of the interaction between the CRISPR-Cas immune system and invaders.     

Multiplexed DNA sequence capture of mitochondrial genomes using PCR products

Background

To utilize the power of high-throughput sequencers, target enrichment methods have been developed. The majority of these require reagents and equipment that are only available from commercial vendors and are not suitable for the targets that are a few kilobases in length.

Methodology/Principal Findings

We describe a novel and economical method in which custom made long-range PCR products are used to capture complete human mitochondrial genomes from complex DNA mixtures. We use the method to capture 46 complete mitochondrial genomes in parallel and we sequence them on a single lane of an Illumina GA_II instrument.

Conclusions/Significance

This method is economical and simple and particularly suitable for targets that can be amplified by PCR and do not contain highly repetitive sequences such as mtDNA. It has applications in population genetics and forensics, as well as studies of ancient DNA.     

Genome size and DNA base composition of geophytes: the mirror of phenology and ecology?

Background and Aims

Genome size is known to affect various plant traits such as stomatal size, seed mass, and flower or shoot phenology. However, these associations are not well understood for species with very large genomes, which are laregly represented by geophytic plants. No detailed associations are known between DNA base composition and genome size or species ecology.

Methods

Genome sizes and GC contents were measured in 219 geophytes together with tentative morpho-anatomical and ecological traits.

Key Results

Increased genome size was associated with earliness of flowering and tendency to grow in humid conditions, and there was a positive correlation between an increase in stomatal size in species with extremely large genomes. Seed mass of geophytes was closely related to their ecology, but not to genomic parameters. Genomic DNA GC content showed a unimodal relationship with genome size but no relationship with species ecology.

Conclusions

Evolution of genome size in geophytes is closely related to their ecology and phenology and is also associated with remarkable changes in DNA base composition. Although geophytism together with producing larger cells appears to be an advantageous strategy for fast development of an organism in seasonal habitats, the drought sensitivity of large stomata may restrict the occurrence of geophytes with very large genomes to regions not subject to water stress.