Functional analysis of the apple fruit microbiome based on shotgun metagenomic sequencing of conventional and organic orchard samples

Fruits harbour abundant and diverse microbial communities that protect them from post-harvest pathogens. Identification of functional traits associated with a given microbiota can provide a better understanding of their potential influence. Here, we focused on the epiphytic microbiome of apple fruit. We suggest that shotgun metagenomic data can indicate specific functions carried out by different groups and provide information on their potential impact. Samples were collected from the surface of ‘Golden Delicious’ apples from four orchards that differ in their geographic location and management practice. Approximately 1 million metagenes were predicted based on a high-quality assembly. Functional profiling of the microbiome of fruits from orchards differing in their management practice revealed a functional shift in the microbiota. The organic orchard microbiome was enriched in pathways involved in plant defence activities; the conventional orchard microbiome was enriched in pathways related to the synthesis of antibiotics. The functional significance of the variations was explored using microbial network modelling algorithms to reveal the metabolic role of specific phylogenetic groups. The analysis identified several associations supported by other published studies. For example, the analysis revealed the nutritional dependencies of the Capnodiales group, including the Alternaria pathogen, on aromatic compounds.     

Milk microbiome signatures of subclinical mastitis-affected cattle analysed by shotgun sequencing

Aims: Metagenomic analysis of milk samples collected from Kankrej, Gir (Bos indicus) and crossbred (Bos taurus × B. indicus) cattle harbouring subclinical mastitis was carried out by next‐generation sequencing 454 GS‐FLX technology to elucidate the microbial community structure of cattle milk. Methods and Results: Milk samples from Kankrej, Gir and crossbred cattle were subjected to metagenomic profiling by pyrosequencing. The Metagenomic analysis produced 63·07, 11·09 and 7·87 million base pairs (Mb) of sequence data, assembled in 264 798, 56 114 and 36 762 sequences with an average read length of 238, 197 and 214 nucleotides in Kankrej, Gir and crossbred cattle, respectively. Phylogenetic and metabolic profiles by the web‐based tool MG‐RAST revealed that the members of Enterobacteriales were predominant in mastitic milk followed by Pseudomonadales, Bacillales and Lactobacillales. Around 56 different species with varying abundance were detected in the subclinically infected milk. Escherichia coli was found to be the most predominant species in Kankrej and Gir cattle followed by Pseudomonas aeruginosa, Pseudomonas mendocina, Shigella flexneri and Bacillus cereus. In crossbred cattle, Staphylococcus aureus followed by Klebsiella pneumoniae, Staphylococcus epidermidis and E. coli were detected in descending order. Metabolic profiling indicated fluoroquinolones, methicillin, copper, cobalt–zinc–cadmium as the groups of antibiotics and toxic compounds to which the organisms showed resistance. Sequences indicating potential of organisms exhibiting multidrug resistance against antibiotics and resistance to toxic compounds were also present. Interestingly, presence of bacteriophages against Staph. aureus, E. coli, Enterobacter and Yersinia species was also observed. Conclusions: The analysis identified potential infectious organisms in mastitis, resistance of organisms to antibiotics and chemical compounds and the natural resistance potential of dairy cows. Significance and Impact of the Study: The findings of this study may help in formulating strategies for the prevention and treatment of mastitis in dairy animals and consequently in reducing economic losses incurred because of it.     

Characterization of the viral microbiome in patients with severe lower respiratory tract infections, using metagenomic sequencing

The human respiratory tract is heavily exposed to microorganisms. Viral respiratory tract pathogens, like RSV, influenza and rhinoviruses cause major morbidity and mortality from respiratory tract disease. Furthermore, as viruses have limited means of transmission, viruses that cause pathogenicity in other tissues may be transmitted through the respiratory tract. It is therefore important to chart the human virome in this compartment. We have studied nasopharyngeal aspirate samples submitted to the Karolinska University Laboratory, Stockholm, Sweden from March 2004 to May 2005 for diagnosis of respiratory tract infections. We have used a metagenomic sequencing strategy to characterize viruses, as this provides the most unbiased view of the samples. Virus enrichment followed by 454 sequencing resulted in totally 703,790 reads and 110,931 of these were found to be of viral origin by using an automated classification pipeline. The snapshot of the respiratory tract virome of these 210 patients revealed 39 species and many more strains of viruses. Most of the viral sequences were classified into one of three major families; Paramyxoviridae, Picornaviridae or Orthomyxoviridae. The study also identified one novel type of Rhinovirus C, and identified a number of previously undescribed viral genetic fragments of unknown origin.     

Shallow shotgun sequencing of the microbiome recapitulates 16S amplicon results and provides functional insights

Prevailing 16S rRNA gene-amplicon methods for characterizing the bacterial microbiome of wildlife are economical, but result in coarse taxonomic classifications, are subject to primer and 16S copy number biases, and do not allow for direct estimation of microbiome functional potential. While deep shotgun metagenomic sequencing can overcome many of these limitations, it is prohibitively expensive for large sample sets. Here we evaluated the ability of shallow shotgun metagenomic sequencing to characterize taxonomic and functional patterns in the faecal microbiome of a model population of feral horses (Sable Island, Canada). Since 2007, this unmanaged population has been the subject of an individual-based, long-term ecological study. Using deep shotgun metagenomic sequencing, we determined the sequencing depth required to accurately characterize the horse microbiome. In comparing conventional vs. high-throughput shotgun metagenomic library preparation techniques, we validate the use of more cost-effective laboratory methods. Finally, we characterize similarities between 16S amplicon and shallow shotgun characterization of the microbiome, and demonstrate that the latter recapitulates biological patterns first described in a published amplicon data set. Unlike for amplicon data, we further demonstrate how shallow shotgun metagenomic data provide useful insights regarding microbiome functional potential which support previously hypothesized diet effects in this study system.     

Faecal microbiome sequences in relation to the egg-laying performance of hens using amplicon-based metagenomic association analysis

Exploring the composition and structure of the faecal microbial community improves the understanding of the role of the gut microbiota in the gastrointestinal function and the egg-laying performance of hens. Therefore, detection of hen–microbial interactions can explore a new breeding marker for the selection of egg production due to the important role of the gut microbiome in the host’s metabolism and health. Recently, the gut microbiota has been recognised as a regulator of host performance, which has led to investigations of the productive effects of changes in the faecal microbiome in various animals. In the present study, a metagenomics analysis was applied to characterise the composition and structural diversity of faecal microbial communities under two selections of egg-laying performance, high (H, n = 30) and low (L, n = 30), using 16S rRNA-based metagenomic association analysis. The most abundant bacterial compositions were estimated based on the operational classification units among samples and between the groups from metagenomic data sets. The results indicated that Firmicutes phylum has higher significant (P < 0.01) in the H group than in the L group. In addition, higher relative abundance phyla of Bacteroides and Fusobacteria were estimated in the H group than the L group, contrasting the phyla of Actinobacteria, Cyanobacteria and Proteobacteria were more relative abundance in the L group. The families (Lactobacillus, Bifidobacterium, Acinetobacter, Flavobacteriaceae, Lachnoclostridum and Rhodococcus) were more abundant in the H group based on the comparison between the H and L groups. Meanwhile, three types of phyla (Proteobacteria, Actinobacteria and Cyanobacteria) and six families (Acinetobacter, Avibacterium, Clostridium, Corynebacterium, Helicobacter and Peptoclostridium) were more abundant in the L group (P < 0.01). Overall, the selection of genotypes has enriched a relationship between the gut microbiota and the egg-laying performance. These findings suggest that the faecal microbiomes of chickens with high egg-laying performance have more diverse activities than those of chickens with low egg-laying performance, which may be related to the metabolism and health of the host and egg production variation.     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

A new era of the vaginal microbiome: advances using next-generation sequencing

Until recently, bacterial species that inhabit the human vagina have been primarily studied using organism-centric approaches. Understanding how these bacterial species interact with each other and the host vaginal epithelium is essential for a more complete understanding of vaginal health. Molecular approaches have already led to the identification of uncultivated bacterial taxa associated with bacterial vaginosis. Here, we review recent studies of the vaginal microbiome and discuss how culture-independent approaches, such as applications of next-generation sequencing, are advancing the field and shifting our understanding of how vaginal health is defined. This work may lead to improved diagnostic tools and treatments for women who suffer from, or are at risk for, vaginal imbalances, pregnancy complications, and sexually acquired infections. These approaches may also transform our understanding of how host genetic factors, physiological conditions (e.g., menopause), and environmental exposures (e.g., smoking, antibiotic usage) influence the vaginal microbiome.     

Neuropeptide Y effects on vasorelaxation and intestinal contraction in the Atlantic cod Gadus morhua

Neuropeptide Y (NPY) has prominent cardiovascular effects in mammals and sharks, but no such effect has previously been demonstrated in any teleost fish. In the Atlantic cod, we found that cod NPY (10(-10)-10(-6) M) relaxed celiac arteries precontracted with epinephrine, and weak contractions were elicited in intestinal ring preparations. A few NPY-immunoreactive nerve fibers were present along small gut arteries. The results suggest that cod NPY produces vasorelaxation both by a direct action on smooth muscle and by release of prostaglandins, but with no involvement of nitric oxide, leukotrienes, or endothelium-derived relaxing factors. An additional indirect effect involving another neurotransmitter may occur. Cod NPY (10(-7) M) and human NPY (10(-7) M) had identical effects on the vessels. Small differences only in the effects of porcine [Leu(31),Pro(34)]NPY, NPY-(13-36), and cod NPY suggest the presence of a Y(1) subfamily receptor, similar to the zebrafish Ya receptor. A physiological role for NPY in teleost vasculature is concluded, but surprisingly the effect, a vasodilation, is opposite to that in mammals and is mediated by prostaglandins.     

A metagenomic approach to characterization of the vaginal microbiome signature in pregnancy

While current major national research efforts (i.e., the NIH Human Microbiome Project) will enable comprehensive metagenomic characterization of the adult human microbiota, how and when these diverse microbial communities take up residence in the host and during reproductive life are unexplored at a population level. Because microbial abundance and diversity might differ in pregnancy, we sought to generate comparative metagenomic signatures across gestational age strata. DNA was isolated from the vagina (introitus, posterior fornix, midvagina) and the V5V3 region of bacterial 16S rRNA genes were sequenced (454FLX Titanium platform). Sixty-eight samples from 24 healthy gravidae (18 to 40 confirmed weeks) were compared with 301 non-pregnant controls (60 subjects). Generated sequence data were quality filtered, taxonomically binned, normalized, and organized by phylogeny and into operational taxonomic units (OTU); principal coordinates analysis (PCoA) of the resultant beta diversity measures were used for visualization and analysis in association with sample clinical metadata. Altogether, 1.4 gigabytes of data containing >2.5 million reads (averaging 6,837 sequences/sample of 493 nt in length) were generated for computational analyses. Although gravidae were not excluded by virtue of a posterior fornix pH >4.5 at the time of screening, unique vaginal microbiome signature encompassing several specific OTUs and higher-level clades was nevertheless observed and confirmed using a combination of phylogenetic, non-phylogenetic, supervised, and unsupervised approaches. Both overall diversity and richness were reduced in pregnancy, with dominance of Lactobacillus species (L. iners crispatus, jensenii and johnsonii, and the orders Lactobacillales (and Lactobacillaceae family), Clostridiales, Bacteroidales, and Actinomycetales. This intergroup comparison using rigorous standardized sampling protocols and analytical methodologies provides robust initial evidence that the vaginal microbial 16S rRNA gene catalogue uniquely differs in pregnancy, with variance of taxa across vaginal subsite and gestational age.     

Assessment of metagenomic assembly using simulated next generation sequencing data

Due to the complexity of the protocols and a limited knowledge of the nature of microbial communities, simulating metagenomic sequences plays an important role in testing the performance of existing tools and data analysis methods with metagenomic data. We developed metagenomic read simulators with platform-specific (Sanger, pyrosequencing, Illumina) base-error models, and simulated metagenomes of differing community complexities. We first evaluated the effect of rigorous quality control on Illumina data. Although quality filtering removed a large proportion of the data, it greatly improved the accuracy and contig lengths of resulting assemblies. We then compared the quality-trimmed Illumina assemblies to those from Sanger and pyrosequencing. For the simple community (10 genomes) all sequencing technologies assembled a similar amount and accurately represented the expected functional composition. For the more complex community (100 genomes) Illumina produced the best assemblies and more correctly resembled the expected functional composition. For the most complex community (400 genomes) there was very little assembly of reads from any sequencing technology. However, due to the longer read length the Sanger reads still represented the overall functional composition reasonably well. We further examined the effect of scaffolding of contigs using paired-end Illumina reads. It dramatically increased contig lengths of the simple community and yielded minor improvements to the more complex communities. Although the increase in contig length was accompanied by increased chimericity, it resulted in more complete genes and a better characterization of the functional repertoire. The metagenomic simulators developed for this research are freely available.     

Effects of cod bradykinin and its analogs on vascular and intestinal smooth muscle of the Atlantic cod, Gadus morhua.

The effects of [Arg(0),Trp(5),Leu(8)]-BK (cod [Arg(0)]BK) on vascular preparations from branches of the cod celiac artery and on longitudinal smooth muscle preparations from the cod intestine were investigated. Cod [Arg(0)]BK (3 x 10(-8) M) caused a relaxation of the celiac artery precontracted with adrenaline. The relaxation was abolished by the cyclooxygenase inhibitor indomethacin, suggesting that the effect is mediated through the release of prostaglandins, but there was no evidence for the involvement of leukotrienes or nitric oxide in the response. In the intestinal preparations, cod [Arg(0)]BK produced concentration-dependent contractions (pD(2) = 8.28 +/- 0.16). Experiments with N-terminally and C-terminally truncated analogs and with alanine-substituted analogs of cod [Arg(0)]BK demonstrate that the central amino acid Gly(4) and the C-terminal amino acids Leu(8) and Arg(9) are the most important in determining the conformation of the peptide that interacts with the receptor. The results indicate that the ligand binding properties of the cod BK receptor are considerably different from the receptor present in trout tissues and may resemble those of the mammalian B(2) receptor more closely.     

Metabolic reconstruction for metagenomic data and its application to the human microbiome

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes including metabolism and immune homeostasis in the human microbiome. Shotgun sequencing of such communities' metagenomes provides information complementary to organismal abundances from taxonomic markers, but the resulting data typically comprise short reads from hundreds of different organisms and are at best challenging to assemble comparably to single-organism genomes. Here, we describe an alternative approach to infer the functional and metabolic potential of a microbial community metagenome. We determined the gene families and pathways present or absent within a community, as well as their relative abundances, directly from short sequence reads. We validated this methodology using a collection of synthetic metagenomes, recovering the presence and abundance both of large pathways and of small functional modules with high accuracy. We subsequently applied this method, HUMAnN, to the microbial communities of 649 metagenomes drawn from seven primary body sites on 102 individuals as part of the Human Microbiome Project (HMP). This provided a means to compare functional diversity and organismal ecology in the human microbiome, and we determined a core of 24 ubiquitously present modules. Core pathways were often implemented by different enzyme families within different body sites, and 168 functional modules and 196 metabolic pathways varied in metagenomic abundance specifically to one or more niches within the microbiome. These included glycosaminoglycan degradation in the gut, as well as phosphate and amino acid transport linked to host phenotype (vaginal pH) in the posterior fornix. An implementation of our methodology is available at http://huttenhower.sph.harvard.edu/humann. This provides a means to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads, enabling the determination of community roles in the HMP cohort and in future metagenomic studies.     

Induction of Mx protein in Atlantic cod with poly I:C: immuno-cross reactive studies of antibodies to Atlantic salmon Mx with Atlantic cod

A polyclonal rabbit antiserum directed against the conserved region of the Atlantic salmon antiviral Mx1 protein was used to detect the putative Atlantic cod Mx protein using Western and dot blotting. A doublet band at about 75kDa and 65kDa was detected by Western blotting in kidney and spleen extracts of cod 3 and 4 days after i.p. injection with poly I:C but not in control fish injected with PBS. In blood leucocyte lysates, similar immunostaining could also be detected in Atlantic cod weakly after injection with PBS and more intensely after injection with poly I:C, suggesting some constitutive expression of Mx protein by leucocytes. Dot blot analysis showed that the Mx protein level was significantly higher in spleen, kidney, liver and gill of cod at least up to 4 days after injection with poly I:C when compared with the PBS-injected controls.     

A practical guide to amplicon and metagenomic analysis of microbiome data

Advances in high-throughput sequencing(HTS)have fostered rapid developments in the field of microbiome research,and massive microbiome datasets are now being generated.However,the diversity of software tools and the complexity of analysis pipelines make it difficult to access this field.Here,we systematically summarize the advantages and limitations of micro-biome methods.Then,we recommend specific pipelines for amplicon and metagenomic analyses,and describe commonly-used software and databases,to help researchers select the appropriate tools.Furthermore,we introduce statistical and visualization methods suit-able for microbiome analysis,including alpha-and beta-diversity,taxonomic composition,difference compar-isons,correlation,networks,machine learning,evolu-tion,source tracing,and common visualization styles to help researchers make informed choices.Finally,a step-by-step reproducible analysis guide is introduced.We hope this review will allow researchers to carry out data analysis more effectively and to quickly select the appropriate tools in order to efficiently mine the bio-logical significance behind the data.     

Susceptibility of juvenile Atlantic cod Gadus morhua to viral haemorrhagic septicaemia virus isolated from wild-caught Atlantic cod

A European strain of viral haemorrhagic septicaemia virus (VHSV) isolated from wild-caught cod Gadus morhua (H16/7/95) was shown to cause clinical disease and mortality in excess of 80% in juvenile Atlantic cod when administered by the intra-peritoneal (i.p.) route. No virus was recovered from cod cohabiting with experimentally infected fish at a ratio of 1:1, and no VHSV-associated mortality was demonstrated following immersion infection. External signs of disease in cod were the presence of exophthalmia and ascites. Virus was identified as VHSV by enzyme-linked immunosorbent assay (ELISA) and was recovered from both brain and organ pools (kidney, liver and spleen) of 100% of i.p. infected cod mortalities. Virus was also detected using an indirect immunofluorescence test on tissue imprints of kidney, liver, spleen and brain taken from moribund fish. The fact that cod were not susceptible to VHSV following waterborne exposure raises important questions surrounding the propagation, maintenance and impact of a naturally occurring reservoir of virus in the marine environment.     

Modeling the feasibility of whole genome shotgun sequencing using a pairwise end strategy

In pairwise end sequencing, sequences are determined from both ends of random subclones derived from a DNA target. Sufficiently similar overlapping end sequences are identified and grouped into contigs. When a clone's paired end sequences fall in different contigs, the contigs are connected together to form scaffolds. Increasingly, the goals of pairwise strategies are large and highly repetitive genomic targets. Here, we consider large-scale pairwise strategies that employ mixtures of subclone sizes. We explore the properties of scaffold formation within a hybrid theory/simulation mathematical model of a genomic target that contains many repeat families. Using this model, we evaluate problems that may arise, such as falsely linked end sequences (due either to random matches or to homologous repeats) and scaffolds that terminate without extending the full length of the target. We illustrate our model with an exploration of a strategy for sequencing the human genome. Our results show that, for a strategy that generates 10-fold sequence coverage derived from the ends of clones ranging in length from 2 to 150 kb, using an appropriate rule for detecting overlaps, we expect few false links while obtaining a single scaffold extending the length of each chromosome.     

RiboTagger: fast and unbiased 16S/18S profiling using whole community shotgun metagenomic or metatranscriptome surveys

Background

Taxonomic profiling of microbial communities is often performed using small subunit ribosomal RNA (SSU) amplicon sequencing (16S or 18S), while environmental shotgun sequencing is often focused on functional analysis. Large shotgun datasets contain a significant number of SSU sequences and these can be exploited to perform an unbiased SSU--based taxonomic analysis.

Results

Here we present a new program called RiboTagger that identifies and extracts taxonomically informative ribotags located in a specified variable region of the SSU gene in a high-throughput fashion.

Conclusions

RiboTagger permits fast recovery of SSU-RNA sequences from shotgun nucleic acid surveys of complex microbial communities. The program targets all three domains of life, exhibits high sensitivity and specificity and is substantially faster than comparable programs.