IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth

Metagenomic Evidence for the Presence of Comammox Nitrospira-Like Bacteria in a Drinking Water System

Nitrification plays an important role in regulating the concentrations of inorganic nitrogen species in a range of environments, from drinking water and wastewater treatment plants to the oceans. Until recently, aerobic nitrification was considered to be a two-step process involving ammonia-oxidizing bacteria or archaea and nitrite-oxidizing bacteria. This process requires close cooperation between these two functional guilds for complete conversion of ammonia to nitrate, without the accumulation of nitrite or other intermediates, such as nitrous oxide, a potent greenhouse gas. The discovery of a single organism with the potential to oxidize both ammonia and nitrite adds a new dimension to the current understanding of aerobic nitrification, while presenting opportunities to rethink nitrogen management in engineered systems.

We report metagenomic evidence for the presence of a Nitrospira-like organism with the metabolic potential to perform the complete oxidation of ammonia to nitrate (i.e., it is a complete ammonia oxidizer [comammox]) in a drinking water system. This metagenome bin was discovered through shotgun DNA sequencing of samples from biologically active filters at the drinking water treatment plant in Ann Arbor, MI. Ribosomal proteins, 16S rRNA, and nxrA gene analyses confirmed that this genome is related to Nitrospira-like nitrite-oxidizing bacteria. The presence of the full suite of ammonia oxidation genes, including ammonia monooxygenase and hydroxylamine dehydrogenase, on a single ungapped scaffold within this metagenome bin suggests the presence of recently discovered comammox potential. Evaluations based on coverage and k-mer frequency distribution, use of two different genome-binning approaches, and nucleic acid and protein similarity analyses support the presence of this scaffold within the Nitrospira metagenome bin. The amoA gene found in this metagenome bin is divergent from those of canonical ammonia and methane oxidizers and clusters closely with the unusual amoA gene of comammox Nitrospira. This finding suggests that previously reported imbalances in abundances of nitrite- and ammonia-oxidizing bacteria/archaea may likely be explained by the capacity of Nitrospira-like organisms to completely oxidize ammonia. This finding might have significant implications for our understanding of microbially mediated nitrogen transformations in engineered and natural systems.

IMPORTANCE Nitrification plays an important role in regulating the concentrations of inorganic nitrogen species in a range of environments, from drinking water and wastewater treatment plants to the oceans. Until recently, aerobic nitrification was considered to be a two-step process involving ammonia-oxidizing bacteria or archaea and nitrite-oxidizing bacteria. This process requires close cooperation between these two functional guilds for complete conversion of ammonia to nitrate, without the accumulation of nitrite or other intermediates, such as nitrous oxide, a potent greenhouse gas. The discovery of a single organism with the potential to oxidize both ammonia and nitrite adds a new dimension to the current understanding of aerobic nitrification, while presenting opportunities to rethink nitrogen management in engineered systems.

Metagenomic investigation of the geologically unique Hellenic Volcanic Arc reveals a distinctive ecosystem with unexpected physiology

Hydrothermal vents represent a deep, hot, aphotic biosphere where chemosynthetic primary producers, fuelled by chemicals from Earth's subsurface, form the basis of life. In this study, we examined microbial mats from two distinct volcanic sites within the Hellenic Volcanic Arc (HVA). The HVA is geologically and ecologically unique, with reported emissions of CO2 -saturated fluids at temperatures up to 220°C and a notable absence of macrofauna. Metagenomic data reveals highly complex prokaryotic communities composed of chemolithoautotrophs, some methanotrophs, and to our surprise, heterotrophs capable of anaerobic degradation of aromatic hydrocarbons. Our data suggest that aromatic hydrocarbons may indeed be a significant source of carbon in these sites, and instigate additional research into the nature and origin of these compounds in the HVA. Novel physiology was assigned to several uncultured prokaryotic lineages; most notably, a SAR406 representative is attributed with a role in anaerobic hydrocarbon degradation. This dataset, the largest to date from submarine volcanic ecosystems, constitutes a significant resource of novel genes and pathways with potential biotechnological applications.

gbtools: Interactive Visualization of Metagenome Bins in R

Improvements in DNA sequencing technology have increased the amount and quality of sequences that can be obtained from metagenomic samples, making it practical to extract individual microbial genomes from metagenomic assemblies (“binning”). However, while many tools and methods exist for unsupervised binning with various statistical algorithms, there are few options for visualizing the results, even though visualization is vital to exploratory data analysis. We have developed gbtools, a software package that allows users to visualize metagenomic assemblies by plotting coverage (sequencing depth) and GC values of contigs, and also to annotate the plots with taxonomic information. Different sets of annotations, including taxonomic assignments from conserved marker genes or SSU rRNA genes, can be imported simultaneously; users can choose which annotations to plot. Bins can be manually defined from plots, or be imported from third-party binning tools and overlaid onto plots, such that results from different methods can be compared side-by-side. gbtools reports summary statistics of bins including marker gene completeness, and allows the user to add or subtract bins with each other. We illustrate some of the functions available in gbtools with two examples: the metagenome of Olavius algarvensis, a marine oligochaete worm that has up to five bacterial symbionts, and the metagenome of a synthetic mock community comprising 64 bacterial and archaeal strains. We show how instances of poor automated binning, sequencer GC% bias, and variation between samples can be quickly diagnosed by visualization, and demonstrate how the results from different binning tools can be combined and refined to yield manually curated bins with higher completeness. gbtools is open-source and written in R. The software package, documentation, and example data are available freely online at https://github.com/kbseah/genome-bin-tools.

The Road to Metagenomics: From Microbiology to DNA Sequencing Technologies and Bioinformatics

The study of microorganisms that pervade each and every part of this planet has encountered many challenges through time such as the discovery of unknown organisms and the understanding of how they interact with their environment. The aim of this review is to take the reader along the timeline and major milestones that led us to modern metagenomics. This new and thriving area is likely to be an important contributor to solve different problems. The transition from classical microbiology to modern metagenomics studies has required the development of new branches of knowledge and specialization. Here, we will review how the availability of high-throughput sequencing technologies has transformed microbiology and bioinformatics and how to tackle the inherent computational challenges that arise from the DNA sequencing revolution. New computational methods are constantly developed to collect, process, and extract useful biological information from a variety of samples and complex datasets, but metagenomics needs the integration of several of these computational methods. Despite the level of specialization needed in bioinformatics, it is important that life-scientists have a good understanding of it for a correct experimental design, which allows them to reveal the information in a metagenome.

A New N-Acyl Homoserine Lactone Synthase in an Uncultured Symbiont of the Red Sea Sponge Theonella swinhoei

Sponges harbor a remarkable diversity of microbial symbionts in which signal molecules can accumulate and enable cell-cell communication, such as quorum sensing (QS). Bacteria capable of QS were isolated from marine sponges; however, an extremely small fraction of the sponge microbiome is amenable to cultivation. We took advantage of community genome assembly and binning to investigate the uncultured majority of sponge symbionts. We identified a complete N-acyl-homoserine lactone (AHL)-QS system (designated TswIR) and seven partial luxI homologues in the microbiome of Theonella swinhoei. The TswIR system was novel and shown to be associated with an alphaproteobacterium of the order Rhodobacterales, here termed Rhodobacterales bacterium TS309. The tswI gene, when expressed in Escherichia coli, produced three AHLs, two of which were also identified in a T. swinhoei sponge extract. The taxonomic affiliation of the 16S rRNA of Rhodobacterales bacterium TS309 to a sponge-coral specific clade, its enrichment in sponge versus seawater and marine sediment samples, and the presence of sponge-specific features, such as ankyrin-like domains and tetratricopeptide repeats, indicate a likely symbiotic nature of this bacterium.

Analysis of metagenomic data reveals common features of halophilic viral communities across continents

Microbial communities from hypersaline ponds, dominated by halophilic archaea, are considered specific of such extreme conditions. The associated viral communities have accordingly been shown to display specific features, such as similar morphologies among different sites. However, little is known about the genetic diversity of these halophilic viral communities across the Earth. Here, we studied viral communities in hypersaline ponds sampled on the coast of Senegal (8-36% of salinity) using metagenomics approach, and compared them with hypersaline viromes from Australia and Spain. The specificity of hyperhalophilic viruses could first be demonstrated at a community scale, salinity being a strong discriminating factor between communities. For the major viral group detected in all samples (Caudovirales), only a limited number of halophilic Caudovirales clades were highlighted. These clades gather viruses from different continents and display consistent genetic composition, indicating that they represent related lineages with a worldwide distribution. Non-tailed hyperhalophilic viruses display a greater rate of gene transfer and recombination, with uncharacterized genes conserved across different kind of viruses and plasmids. Thus, hypersaline viral communities around the world appear to form a genetically consistent community that are likely to harbour new genes coding for enzymes specifically adapted to these environments.

A viability-linked metagenomic analysis of cleanroom environments: eukarya, prokaryotes, and viruses

BACKGROUND

Recent studies posit a reciprocal dependency between the microbiomes associated with humans and indoor environments. However, none of these metagenome surveys has considered the viability of constituent microorganisms when inferring impact on human health.

RESULTS

Reported here are the results of a viability-linked metagenomics assay, which (1) unveil a remarkably complex community profile for bacteria, fungi, and viruses and (2) bolster the detection of underrepresented taxa by eliminating biases resulting from extraneous DNA. This approach enabled, for the first time ever, the elucidation of viral genomes from a cleanroom environment. Upon comparing the viable biomes and distribution of phylotypes within a cleanroom and adjoining (uncontrolled) gowning enclosure, the rigorous cleaning and stringent control countermeasures of the former were observed to select for a greater presence of anaerobes and spore-forming microflora. Sequence abundance and correlation analyses suggest that the viable indoor microbiome is influenced by both the human microbiome and the surrounding ecosystem(s).

CONCLUSIONS

The findings of this investigation constitute the literature's first ever account of the indoor metagenome derived from DNA originating solely from the potential viable microbial population. Results presented in this study should prove valuable to the conceptualization and experimental design of future studies on indoor microbiomes aimed at inferring impact on human health.

Complete Genome Sequence of a New Member of the Marseilleviridae Recovered from the Brackish Submarine Spring in the Cassis Port-Miou Calanque, France

Marseilleviridae is a rapidly expanding family of Acanthamoeba-infecting large DNA viruses distributed worldwide. We report here the complete 349-kbp genome sequence of Port-Miou virus, which is surprisingly close to that of Lausannevirus (isolated from the Seine River upstream from Paris, France), despite the strong dissimilarities of their sampling locations.

Metagenomics of an Alkaline Hot Spring in Galicia (Spain): Microbial Diversity Analysis and Screening for Novel Lipolytic Enzymes

A fosmid library was constructed with the metagenomic DNA from the water of the Lobios hot spring (76°C, pH = 8.2) located in Ourense (Spain). Metagenomic sequencing of the fosmid library allowed the assembly of 9722 contigs ranging in size from 500 to 56,677 bp and spanning ~18 Mbp. 23,207 ORFs (Open Reading Frames) were predicted from the assembly. Biodiversity was explored by taxonomic classification and it revealed that bacteria were predominant, while the archaea were less abundant. The six most abundant bacterial phyla were Deinococcus-Thermus, Proteobacteria, Firmicutes, Acidobacteria, Aquificae, and Chloroflexi. Within the archaeal superkingdom, the phylum Thaumarchaeota was predominant with the dominant species “Candidatus Caldiarchaeum subterraneum.” Functional classification revealed the genes associated to one-carbon metabolism as the most abundant. Both taxonomic and functional classifications showed a mixture of different microbial metabolic patterns: aerobic and anaerobic, chemoorganotrophic and chemolithotrophic, autotrophic and heterotrophic. Remarkably, the presence of genes encoding enzymes with potential biotechnological interest, such as xylanases, galactosidases, proteases, and lipases, was also revealed in the metagenomic library. Functional screening of this library was subsequently done looking for genes encoding lipolytic enzymes. Six genes conferring lipolytic activity were identified and one was cloned and characterized. This gene was named LOB4Est and it was expressed in a yeast mesophilic host. LOB4Est codes for a novel esterase of family VIII, with sequence similarity to β-lactamases, but with unusual wide substrate specificity. When the enzyme was purified from the mesophilic host it showed half-life of 1 h and 43 min at 50°C, and maximal activity at 40°C and pH 7.5 with p-nitrophenyl-laurate as substrate. Interestingly, the enzyme retained more than 80% of maximal activity in a broad range of pH from 6.5 to 8.

Use of Substrate-Induced Gene Expression in Metagenomic Analysis of an Aromatic Hydrocarbon-Contaminated Soil

Metagenomics allows the study of genes related to xenobiotic degradation in a culture-independent manner, but many of these studies are limited by the lack of genomic context for metagenomic sequences. This study combined a phenotypic screen known as substrate-induced gene expression (SIGEX) with whole-metagenome shotgun sequencing. SIGEX is a high-throughput promoter-trap method that relies on transcriptional activation of a green fluorescent protein (GFP) reporter gene in response to an inducing compound and subsequent fluorescence-activated cell sorting to isolate individual inducible clones from a metagenomic DNA library. We describe a SIGEX procedure with improved library construction from fragmented metagenomic DNA and improved flow cytometry sorting procedures. We used SIGEX to interrogate an aromatic hydrocarbon (AH)-contaminated soil metagenome. The recovered clones contained sequences with various degrees of similarity to genes (or partial genes) involved in aromatic metabolism, for example, nahG (salicylate oxygenase) family genes and their respective upstream nahR regulators. To obtain a broader context for the recovered fragments, clones were mapped to contigs derived from de novo assembly of shotgun-sequenced metagenomic DNA which, in most cases, contained complete operons involved in aromatic metabolism, providing greater insight into the origin of the metagenomic fragments. A comparable set of contigs was generated using a significantly less computationally intensive procedure in which assembly of shotgun-sequenced metagenomic DNA was directed by the SIGEX-recovered sequences. This methodology may have broad applicability in identifying biologically relevant subsets of metagenomes (including both novel and known sequences) that can be targeted computationally by in silico assembly and prediction tools.

hybridSPAdes: an algorithm for hybrid assembly of short and long reads

MOTIVATION

Recent advances in single molecule real-time (SMRT) and nanopore sequencing technologies have enabled high-quality assemblies from long and inaccurate reads. However, these approaches require high coverage by long reads and remain expensive. On the other hand, the inexpensive short reads technologies produce accurate but fragmented assemblies. Thus, a hybrid approach that assembles long reads (with low coverage) and short reads has a potential to generate high-quality assemblies at reduced cost.

RESULTS

We describe hybridSPAdes algorithm for assembling short and long reads and benchmark it on a variety of bacterial assembly projects. Our results demonstrate that hybridSPAdes generates accurate assemblies (even in projects with relatively low coverage by long reads) thus reducing the overall cost of genome sequencing. We further present the first complete assembly of a genome from single cells using SMRT reads.

AVAILABILITY AND IMPLEMENTATION

hybridSPAdes is implemented in C++ as a part of SPAdes genome assembler and is publicly available at http://bioinf.spbau.ru/en/spades

CONTACT

d.antipov@spbu.ru

SUPPLEMENTARY INFORMATION

supplementary data are available at Bioinformatics online.

Draft Genome Sequence of the Moderately Heat-Tolerant Lactococcus lactis subsp. lactis bv. diacetylactis Strain GL2 from Algerian Dromedary Milk

Lactococcus lactis subsp. lactis bv. diacetylactis GL2 is a moderately thermotolerant lactic acid bacterium isolated from dromedary raw milk. Here, we present the draft genome sequence of this potential new dairy starter strain, which combines thermotolerance and the capacity to metabolize lactose, casein, and citrate.

Metagenome Analysis Exploiting High-Throughput Chromosome Conformation Capture (3C) Data

Microbial communities are complex and constitute important parts of our environment. Genomic analysis of these populations is a dynamic research area but remains limited by the difficulty in assembling full genomes of individual species. Recently, a new method for metagenome assembly/analysis based on chromosome conformation capture has emerged (meta3C). This approach quantifies the collisions experienced by DNA molecules to identify those sharing the same cellular compartments, allowing the characterization of genomes present within complex mixes of species. The exploitation of these chromosome 3D signatures holds promising perspectives for genome sequencing of discrete species in complex populations. It also has the potential to assign correctly extra-chromosomal elements, such as plasmids, mobile elements and phages, to their host cells.

The Utility of Genome Skimming for Phylogenomic Analyses as Demonstrated for Glycerid Relationships (Annelida, Glyceridae)

Glyceridae (Annelida) are a group of venomous annelids distributed worldwide from intertidal to abyssal depths. To trace the evolutionary history and complexity of glycerid venom cocktails, a solid backbone phylogeny of this group is essential. We therefore aimed to reconstruct the phylogenetic relationships of these annelids using Illumina sequencing technology. We constructed whole-genome shotgun libraries for 19 glycerid specimens and 1 outgroup species (Glycinde armigera). The chosen target genes comprise 13 mitochondrial proteins, 2 ribosomal mitochondrial genes, and 4 nuclear loci (18SrRNA, 28SrRNA, ITS1, and ITS2). Based on partitioned maximum likelihood as well as Bayesian analyses of the resulting supermatrix, we were finally able to resolve a robust glycerid phylogeny and identified three clades comprising the majority of taxa. Furthermore, we detected group II introns inside the cox1 gene of two analyzed glycerid specimens, with two different insertions in one of these species. Moreover, we generated reduced data sets comprising 10 million, 4 million, and 1 million reads from the original data sets to test the influence of the sequencing depth on assembling complete mitochondrial genomes from low coverage genome data. We estimated the coverage of mitochondrial genome sequences in each data set size by mapping the filtered Illumina reads against the respective mitochondrial contigs. By comparing the contig coverage calculated in all data set sizes, we got a hint for the scalability of our genome skimming approach. This allows estimating more precisely the number of reads that are at least necessary to reconstruct complete mitochondrial genomes in Glyceridae and probably non-model organisms in general.

Genomic and transcriptomic evidence for scavenging of diverse organic compounds by widespread deep-sea archaea

Microbial activity is one of the most important processes to mediate the flux of organic carbon from the ocean surface to the seafloor. However, little is known about the microorganisms that underpin this key step of the global carbon cycle in the deep oceans. Here we present genomic and transcriptomic evidence that five ubiquitous archaeal groups actively use proteins, carbohydrates, fatty acids and lipids as sources of carbon and energy at depths ranging from 800 to 4,950 m in hydrothermal vent plumes and pelagic background seawater across three different ocean basins. Genome-enabled metabolic reconstructions and gene expression patterns show that these marine archaea are motile heterotrophs with extensive mechanisms for scavenging organic matter. Our results shed light on the ecological and physiological properties of ubiquitous marine archaea and highlight their versatile metabolic strategies in deep oceans that might play a critical role in global carbon cycling.

misFinder: identify mis-assemblies in an unbiased manner using reference and paired-end reads

Background

Because of the short read length of high throughput sequencing data, assembly errors are introduced in genome assembly, which may have adverse impact to the downstream data analysis. Several tools have been developed to eliminate these errors by either 1) comparing the assembled sequences with some similar reference genome, or 2) analyzing paired-end reads aligned to the assembled sequences and determining inconsistent features alone mis-assembled sequences. However, the former approach cannot distinguish real structural variations between the target genome and the reference genome while the latter approach could have many false positive detections (correctly assembled sequence being considered as mis-assembled sequence).

Results

We present misFinder, a tool that aims to identify the assembly errors with high accuracy in an unbiased way and correct these errors at their mis-assembled positions to improve the assembly accuracy for downstream analysis. It combines the information of reference (or close related reference) genome and aligned paired-end reads to the assembled sequence. Assembly errors and correct assemblies corresponding to structural variations can be detected by comparing the genome reference and assembled sequence. Different types of assembly errors can then be distinguished from the mis-assembled sequence by analyzing the aligned paired-end reads using multiple features derived from coverage and consistence of insert distance to obtain high confident error calls.

Conclusions

We tested the performance of misFinder on both simulated and real paired-end reads data, and misFinder gave accurate error calls with only very few miscalls. And, we further compared misFinder with QUAST and REAPR. misFinder outperformed QUAST and REAPR by 1) identified more true positive mis-assemblies with very few false positives and false negatives, and 2) distinguished the correct assemblies corresponding to structural variations from mis-assembled sequence. misFinder can be freely downloaded from https://github.com/hitbio/misFinder.

Electronic supplementary material

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

Reconstructing 16S rRNA genes in metagenomic data

Metagenomic data, which contains sequenced DNA reads of uncultured microbial species from environmental samples, provide a unique opportunity to thoroughly analyze microbial species that have never been identified before. Reconstructing 16S ribosomal RNA, a phylogenetic marker gene, is usually required to analyze the composition of the metagenomic data. However, massive volume of dataset, high sequence similarity between related species, skewed microbial abundance and lack of reference genes make 16S rRNA reconstruction difficult. Generic de novo assembly tools are not optimized for assembling 16S rRNA genes. In this work, we introduce a targeted rRNA assembly tool, REAGO (REconstruct 16S ribosomal RNA Genes from metagenOmic data). It addresses the above challenges by combining secondary structure-aware homology search, zproperties of rRNA genes and de novo assembly. Our experimental results show that our tool can correctly recover more rRNA genes than several popular generic metagenomic assembly tools and specially designed rRNA construction tools.

AVAILABILITY AND IMPLEMENTATION

The source code of REAGO is freely available at https://github.com/chengyuan/reago.

De novo meta-assembly of ultra-deep sequencing data

We introduce a new divide and conquer approach to deal with the problem of de novo genome assembly in the presence of ultra-deep sequencing data (i.e. coverage of 1000x or higher). Our proposed meta-assembler Slicembler partitions the input data into optimal-sized ‘slices’ and uses a standard assembly tool (e.g. Velvet, SPAdes, IDBA_UD and Ray) to assemble each slice individually. Slicembler uses majority voting among the individual assemblies to identify long contigs that can be merged to the consensus assembly. To improve its efficiency, Slicembler uses a generalized suffix tree to identify these frequent contigs (or fraction thereof). Extensive experimental results on real ultra-deep sequencing data (8000x coverage) and simulated data show that Slicembler significantly improves the quality of the assembly compared with the performance of the base assembler. In fact, most of the times, Slicembler generates error-free assemblies. We also show that Slicembler is much more resistant against high sequencing error rate than the base assembler.

Availability and implementation: Slicembler can be accessed at http://slicembler.cs.ucr.edu/.

Contact:hamid.mirebrahim@email.ucr.edu

Misassembly detection using paired-end sequence reads and optical mapping data

Motivation: A crucial problem in genome assembly is the discovery and correction of misassembly errors in draft genomes. We develop a method called misSEQuel that enhances the quality of draft genomes by identifying misassembly errors and their breakpoints using paired-end sequence reads and optical mapping data. Our method also fulfills the critical need for open source computational methods for analyzing optical mapping data. We apply our method to various assemblies of the loblolly pine, Francisella tularensis, rice and budgerigar genomes. We generated and used stimulated optical mapping data for loblolly pine and F.tularensis and used real optical mapping data for rice and budgerigar.

Results: Our results demonstrate that we detect more than 54% of extensively misassembled contigs and more than 60% of locally misassembled contigs in assemblies of F.tularensis and between 31% and 100% of extensively misassembled contigs and between 57% and 73% of locally misassembled contigs in assemblies of loblolly pine. Using the real optical mapping data, we correctly identified 75% of extensively misassembled contigs and 100% of locally misassembled contigs in rice, and 77% of extensively misassembled contigs and 80% of locally misassembled contigs in budgerigar.

Availability and implementation:misSEQuel can be used as a post-processing step in combination with any genome assembler and is freely available at http://www.cs.colostate.edu/seq/.

Contact:muggli@cs.colostate.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Chemostat culture systems support diverse bacteriophage communities from human feces

BACKGROUND

Most human microbiota studies focus on bacteria inhabiting body surfaces, but these surfaces also are home to large populations of viruses. Many are bacteriophages, and their role in driving bacterial diversity is difficult to decipher without the use of in vitro ecosystems that can reproduce human microbial communities.

RESULTS

We used chemostat culture systems known to harbor diverse fecal bacteria to decipher whether these cultures also are home to phage communities. We found that there are vast viral communities inhabiting these ecosystems, with estimated concentrations similar to those found in human feces. The viral communities are composed entirely of bacteriophages and likely contain both temperate and lytic phages based on their similarities to other known phages. We examined the cultured phage communities at five separate time points over 24 days and found that they were highly individual-specific, suggesting that much of the subject-specificity found in human viromes also is captured by this culture-based system. A high proportion of the community membership is conserved over time, but the cultured communities maintain more similarity with other intra-subject cultures than they do to human feces. In four of the five subjects, estimated viral diversity between fecal and cultured communities was highly similar.

CONCLUSIONS

Because the diversity of phages in these cultured fecal communities have similarities to those found in humans, we believe these communities can serve as valuable ecosystems to help uncover the role of phages in human microbial communities.

Experimental Design and Bioinformatics Analysis for the Application of Metagenomics in Environmental Sciences and Biotechnology

Recent advances in DNA sequencing technologies have prompted the widespread application of metagenomics for the investigation of novel bioresources (e.g., industrial enzymes and bioactive molecules) and unknown biohazards (e.g., pathogens and antibiotic resistance genes) in natural and engineered microbial systems across multiple disciplines. This review discusses the rigorous experimental design and sample preparation in the context of applying metagenomics in environmental sciences and biotechnology. Moreover, this review summarizes the principles, methodologies, and state-of-the-art bioinformatics procedures, tools and database resources for metagenomics applications and discusses two popular strategies (analysis of unassembled reads versus assembled contigs/draft genomes) for quantitative or qualitative insights of microbial community structure and functions. Overall, this review aims to facilitate more extensive application of metagenomics in the investigation of uncultured microorganisms, novel enzymes, microbe-environment interactions, and biohazards in biotechnological applications where microbial communities are engineered for bioenergy production, wastewater treatment, and bioremediation.

Subtractive assembly for comparative metagenomics, and its application to type 2 diabetes metagenomes

Comparative metagenomics remains challenging due to the size and complexity of metagenomic datasets. Here we introduce subtractive assembly, a de novo assembly approach for comparative metagenomics that directly assembles only the differential reads that distinguish between two groups of metagenomes. Using simulated datasets, we show it improves both the efficiency of the assembly and the assembly quality of the differential genomes and genes. Further, its application to type 2 diabetes (T2D) metagenomic datasets reveals clear signatures of the T2D gut microbiome, revealing new phylogenetic and functional features of the gut microbial communities associated with T2D.

Identification and Resolution of Microdiversity through Metagenomic Sequencing of Parallel Consortia

To gain a predictive understanding of the interspecies interactions within microbial communities that govern community function, the genomic complement of every member population must be determined. Although metagenomic sequencing has enabled the de novo reconstruction of some microbial genomes from environmental communities, microdiversity confounds current genome reconstruction techniques. To overcome this issue, we performed short-read metagenomic sequencing on parallel consortia, defined as consortia cultivated under the same conditions from the same natural community with overlapping species composition. The differences in species abundance between the two consortia allowed reconstruction of near-complete (at an estimated >85% of gene complement) genome sequences for 17 of the 20 detected member species. Two Halomonas spp. indistinguishable by amplicon analysis were found to be present within the community. In addition, comparison of metagenomic reads against the consensus scaffolds revealed within-species variation for one of the Halomonas populations, one of the Rhodobacteraceae populations, and the Rhizobiales population. Genomic comparison of these representative instances of inter- and intraspecies microdiversity suggests differences in functional potential that may result in the expression of distinct roles in the community. In addition, isolation and complete genome sequence determination of six member species allowed an investigation into the sensitivity and specificity of genome reconstruction processes, demonstrating robustness across a wide range of sequence coverage (9× to 2,700×) within the metagenomic data set.

Anvi'o: an advanced analysis and visualization platform for 'omics data

Advances in high-throughput sequencing and ‘omics technologies are revolutionizing studies of naturally occurring microbial communities. Comprehensive investigations of microbial lifestyles require the ability to interactively organize and visualize genetic information and to incorporate subtle differences that enable greater resolution of complex data. Here we introduce anvi’o, an advanced analysis and visualization platform that offers automated and human-guided characterization of microbial genomes in metagenomic assemblies, with interactive interfaces that can link ‘omics data from multiple sources into a single, intuitive display. Its extensible visualization approach distills multiple dimensions of information about each contig, offering a dynamic and unified work environment for data exploration, manipulation, and reporting. Using anvi’o, we re-analyzed publicly available datasets and explored temporal genomic changes within naturally occurring microbial populations through de novo characterization of single nucleotide variations, and linked cultivar and single-cell genomes with metagenomic and metatranscriptomic data. Anvi’o is an open-source platform that empowers researchers without extensive bioinformatics skills to perform and communicate in-depth analyses on large ‘omics datasets.

Draft genome sequence of Halomonas meridiana R1t3 isolated from the surface microbiota of the Caribbean Elkhorn coral Acropora palmata

Members of the gammaproteobacterial genus Halomonas are common in marine environments. Halomonas and other members of the Oceanospirillales have recently been identified as prominent members of the surface microbiota of reef-building corals. Halomonas meridiana strain R1t3 was isolated from the surface mucus layer of the scleractinian coral Acropora palmata in 2005 from the Florida Keys. This strain was chosen for genome sequencing to provide insight into the role of commensal heterotrophic bacteria in the coral holobiont. The draft genome consists of 290 scaffolds, totaling 3.5 Mbp in length and contains 3397 protein-coding genes.

Sequencing of 15 622 gene-bearing BACs clarifies the gene-dense regions of the barley genome

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole-genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene-containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical-mapped gene-bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene-enriched BACs and are characterized by high recombination rates, there are also gene-dense regions with suppressed recombination. We made use of published map-anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D-genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley-Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map-based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene-dense but low recombination is particularly relevant.

Time-series metagenomic analysis reveals robustness of soil microbiome against chemical disturbance

Soil microbial communities have great potential for bioremediation of recalcitrant aromatic compounds. However, it is unclear which taxa and genes in the communities, and how they contribute to the bioremediation in the polluted soils. To get clues about this fundamental question here, time-course (up to 24 weeks) metagenomic analysis of microbial community in a closed soil microcosm artificially polluted with four aromatic compounds, including phenanthrene, was conducted to investigate the changes in the community structures and gene pools. The pollution led to drastic changes in the community structures and the gene sets for pollutant degradation. Complete degradation of phenanthrene was strongly suggested to occur by the syntrophic metabolism by Mycobacterium and the most proliferating genus, Burkholderia. The community structure at Week 24 (∼12 weeks after disappearance of the pollutants) returned to the structure similar to that before pollution. Our time-course metagenomic analysis of phage genes strongly suggested the involvement of the ‘kill-the-winner’ phenomenon (i.e. phage predation of Burkholderia cells) for the returning of the microbial community structure. The pollution resulted in a decrease in taxonomic diversity and a drastic increase in diversity of gene pools in the communities, showing the functional redundancy and robustness of the communities against chemical disturbance.

Metagenomic evidence for reciprocal particle exchange between the mainstem estuary and lateral bay sediments of the lower Columbia River

Lateral bays of the lower Columbia River estuary are areas of enhanced water retention that influence net ecosystem metabolism through activities of their diverse microbial communities. Metagenomic characterization of sediment microbiota from three disparate sites in two brackish lateral bays (Baker and Youngs) produced ∼100 Gbp of DNA sequence data analyzed subsequently for predicted SSU rRNA and peptide-coding genes. The metagenomes were dominated by Bacteria. A large component of Eukaryota was present in Youngs Bay samples, i.e., the inner bay sediment was enriched with the invasive New Zealand mudsnail, Potamopyrgus antipodarum, known for high ammonia production. The metagenome was also highly enriched with an archaeal ammonia oxidizer closely related to Nitrosoarchaeum limnia. Combined analysis of sequences and continuous, high-resolution time series of biogeochemical data from fixed and mobile platforms revealed the importance of large-scale reciprocal particle exchanges between the mainstem estuarine water column and lateral bay sediments. Deposition of marine diatom particles in sediments near Youngs Bay mouth was associated with a dramatic enrichment of Bacteroidetes (58% of total Bacteria) and corresponding genes involved in phytoplankton polysaccharide degradation. The Baker Bay sediment metagenome contained abundant Archaea, including diverse methanogens, as well as functional genes for methylotrophy and taxonomic markers for syntrophic bacteria, suggesting that active methane cycling occurs at this location. Our previous work showed enrichments of similar anaerobic taxa in particulate matter of the mainstem estuarine water column. In total, our results identify the lateral bays as both sources and sinks of biogenic particles significantly impacting microbial community composition and biogeochemical activities in the estuary.

Draft Genome Sequence of Erwinia billingiae OSU19-1, Isolated from a Pear Tree Canker

Plant-associated Erwinia include pathogenic and nonpathogenic species. We report the 5.6-Mb genome sequence of Erwinia billingiae OSU19-1, isolated from a canker on a pear tree inoculated with Erwinia amylovora. OSU19-1 and a closely related European isolate, E. billingiae Eb661^T, share many similarities including 40 kb of plasmid sequence.

Genome sequence of Anoxybacillus ayderensis AB04(T) isolated from the Ayder hot spring in Turkey

Species of Anoxybacillus are thermophiles and, therefore, their enzymes are suitable for many biotechnological applications. Anoxybacillus ayderensis AB04(T) (= NCIMB 13972(T) = NCCB 100050(T)) was isolated from the Ayder hot spring in Rize, Turkey, and is one of the earliest described Anoxybacillus type strains. The present work reports the cellular features of A. ayderensis AB04(T), together with a high-quality draft genome sequence and its annotation. The genome is 2,832,347 bp long (74 contigs) and contains 2,895 protein-coding sequences and 103 RNA genes including 14 rRNAs, 88 tRNAs, and 1 tmRNA. Based on the genome annotation of strain AB04(T), we identified genes encoding various glycoside hydrolases that are important for carbohydrate-related industries, which we compared with those of other, sequenced Anoxybacillus spp. Insights into under-explored industrially applicable enzymes and the possible applications of strain AB04(T) were also described.

Enhanced virome sequencing using targeted sequence capture

Metagenomic shotgun sequencing (MSS) is an important tool for characterizing viral populations. It is culture independent, requires no a priori knowledge of the viruses in the sample, and may provide useful genomic information. However, MSS can lack sensitivity and may yield insufficient data for detailed analysis. We have created a targeted sequence capture panel, ViroCap, designed to enrich nucleic acid from DNA and RNA viruses from 34 families that infect vertebrate hosts. A computational approach condensed ∼1 billion bp of viral reference sequence into <200 million bp of unique, representative sequence suitable for targeted sequence capture. We compared the effectiveness of detecting viruses in standard MSS versus MSS following targeted sequence capture. First, we analyzed two sets of samples, one derived from samples submitted to a diagnostic virology laboratory and one derived from samples collected in a study of fever in children. We detected 14 and 18 viruses in the two sets, comprising 19 genera from 10 families, with dramatic enhancement of genome representation following capture enrichment. The median fold-increases in percentage viral reads post-capture were 674 and 296. Median breadth of coverage increased from 2.1% to 83.2% post-capture in the first set and from 2.0% to 75.6% in the second set. Next, we analyzed samples containing a set of diverse anellovirus sequences and demonstrated that ViroCap could be used to detect viral sequences with up to 58% variation from the references used to select capture probes. ViroCap substantially enhances MSS for a comprehensive set of viruses and has utility for research and clinical applications.

Metatranscriptomic insights into polyphosphate metabolism in marine sediments

Microorganisms can influence inorganic phosphate (Pi) in pore waters, and thus the saturation state of phosphatic minerals, by accumulating and hydrolyzing intracellular polyphosphate (poly-P). Here we used comparative metatranscriptomics to explore microbial poly-P utilization in marine sediments. Sulfidic marine sediments from methane seeps near Barbados and from the Santa Barbara Basin (SBB) oxygen minimum zone were incubated under oxic and anoxic sulfidic conditions. Pi was sequestered under oxic conditions and liberated under anoxic conditions. Transcripts homologous to poly-P kinase type 2 (ppk2) were 6-22 × more abundant in metatranscriptomes from the anoxic incubations, suggesting that reversible poly-P degradation by Ppk2 may be an important metabolic response to anoxia by marine microorganisms. Overall, diverse taxa differentially expressed homologues of genes for poly-P degradation (ppk2 and exopolyphosphatase) under different incubation conditions. Sulfur-oxidizing microorganisms appeared to preferentially express genes for poly-P degradation under anoxic conditions, which may impact phosphorus cycling in a wide range of oxygen-depleted marine settings.

Draft Genome Sequence of an Alphaproteobacterium Associated with the Mediterranean Sponge Oscarella lobularis

While sequencing DNA purified from the homoscleromorph sponge Oscarella lobularis, we detected a large number of reads with strong similarity to available alphaproteobacteria gene sequences of family Rhodobacteraceae. Here, we present the genome sequence of this putative sponge symbiont that we propose to designate as "Candidatus Rhodobacter lobularis."

Environmental genes and genomes: understanding the differences and challenges in the approaches and software for their analyses

DNA-based taxonomic and functional profiling is widely used for the characterization of organismal communities across a rapidly increasing array of research areas that include the role of microbiomes in health and disease, biomonitoring, and estimation of both microbial and metazoan species richness. Two principal approaches are currently used to assign taxonomy to DNA sequences: DNA metabarcoding and metagenomics. When initially developed, each of these approaches mandated their own particular methods for data analysis; however, with the development of high-throughput sequencing (HTS) techniques they have begun to share many aspects in data set generation and processing. In this review we aim to define the current characteristics, goals and boundaries of each field, and describe the different software used for their analysis. We argue that an appreciation of the potential and limitations of each method can help underscore the improvements required by each field so as to better exploit the richness of current HTS-based data sets.

Discovery of phosphonic acid natural products by mining the genomes of 10,000 actinomycetes

Although natural products have been a particularly rich source of human medicines, activity-based screening results in a very high rate of rediscovery of known molecules. Based on the large number of natural product biosynthetic genes in microbial genomes, many have proposed "genome mining" as an alternative approach for discovery efforts; however, this idea has yet to be performed experimentally on a large scale. Here, we demonstrate the feasibility of large-scale, high-throughput genome mining by screening a collection of over 10,000 actinomycetes for the genetic potential to make phosphonic acids, a class of natural products with diverse and useful bioactivities. Genome sequencing identified a diverse collection of phosphonate biosynthetic gene clusters within 278 strains. These clusters were classified into 64 distinct groups, of which 55 are likely to direct the synthesis of unknown compounds. Characterization of strains within five of these groups resulted in the discovery of a new archetypical pathway for phosphonate biosynthesis, the first (to our knowledge) dedicated pathway for H-phosphinates, and 11 previously undescribed phosphonic acid natural products. Among these compounds are argolaphos, a broad-spectrum antibacterial phosphonopeptide composed of aminomethylphosphonate in peptide linkage to a rare amino acid N(5)-hydroxyarginine; valinophos, an N-acetyl l-Val ester of 2,3-dihydroxypropylphosphonate; and phosphonocystoximate, an unusual thiohydroximate-containing molecule representing a new chemotype of sulfur-containing phosphonate natural products. Analysis of the genome sequences from the remaining strains suggests that the majority of the phosphonate biosynthetic repertoire of Actinobacteria has been captured at the gene level. This dereplicated strain collection now provides a reservoir of numerous, as yet undiscovered, phosphonate natural products.

Utilizing de Bruijn graph of metagenome assembly for metatranscriptome analysis

Motivation: Metagenomics research has accelerated the studies of microbial organisms, providing insights into the composition and potential functionality of various microbial communities. Metatranscriptomics (studies of the transcripts from a mixture of microbial species) and other meta-omics approaches hold even greater promise for providing additional insights into functional and regulatory characteristics of the microbial communities. Current metatranscriptomics projects are often carried out without matched metagenomic datasets (of the same microbial communities). For the projects that produce both metatranscriptomic and metagenomic datasets, their analyses are often not integrated. Metagenome assemblies are far from perfect, partially explaining why metagenome assemblies are not used for the analysis of metatranscriptomic datasets.

Results: Here, we report a reads mapping algorithm for mapping of short reads onto a de Bruijn graph of assemblies. A hash table of junction k-mers (k-mers spanning branching structures in the de Bruijn graph) is used to facilitate fast mapping of reads to the graph. We developed an application of this mapping algorithm: a reference-based approach to metatranscriptome assembly using graphs of metagenome assembly as the reference. Our results show that this new approach (called TAG) helps to assemble substantially more transcripts that otherwise would have been missed or truncated because of the fragmented nature of the reference metagenome.

Availability and implementation: TAG was implemented in C++ and has been tested extensively on the Linux platform. It is available for download as open source at http://omics.informatics.indiana.edu/TAG.

Contact:yye@indiana.edu

Does human papillomavirus-negative condylomata exist?

Condylomata acuminata is caused by human papillomavirus (HPV). PCR with consensus primers will typically detect HPV in >96% of condylomata. Metagenomic sequencing has found that some "HPV-negative" condylomata do indeed contain HPV. We wished to perform a renewed evaluation of the "HPV-negative" condylomata using deeper metagenomics sequencing. Sequencing of whole genome amplified DNA from 40 apparently "HPV-negative" condylomata detected HPV in 37/40 specimens. We found 75 different HPV types, out of which 43 represented novel putative HPV types. Three types were cloned and established as HPV types 200, 201 and 202. Molluscum contagiosum virus was detected in 24 of the 40 samples. In summary, deep sequencing enables detection of HPV in almost all condylomata. "HPV-negative" condylomata might largely be explained by clinical misdiagnosis or the presence of viral variants, distantly related HPV types and/or low viral loads.

Extensive Genomic Diversity among Bovine-Adapted Staphylococcus aureus: Evidence for a Genomic Rearrangement within CC97

Staphylococcus aureus is an important pathogen associated with both human and veterinary disease and is a common cause of bovine mastitis. Genomic heterogeneity exists between S. aureus strains and has been implicated in the adaptation of specific strains to colonise particular mammalian hosts. Knowledge of the factors required for host specificity and virulence is important for understanding the pathogenesis and management of S. aureus mastitis. In this study, a panel of mastitis-associated S. aureus isolates (n = 126) was tested for resistance to antibiotics commonly used to treat mastitis. Over half of the isolates (52%) demonstrated resistance to penicillin and ampicillin but all were susceptible to the other antibiotics tested. S. aureus isolates were further examined for their clonal diversity by Multi-Locus Sequence Typing (MLST). In total, 18 different sequence types (STs) were identified and eBURST analysis demonstrated that the majority of isolates grouped into clonal complexes CC97, CC151 or sequence type (ST) 136. Analysis of the role of recombination events in determining S. aureus population structure determined that ST diversification through nucleotide substitutions were more likely to be due to recombination compared to point mutation, with regions of the genome possibly acting as recombination hotspots. DNA microarray analysis revealed a large number of differences amongst S. aureus STs in their variable genome content, including genes associated with capsule and biofilm formation and adhesion factors. Finally, evidence for a genomic arrangement was observed within isolates from CC97 with the ST71-like subgroup showing evidence of an IS431 insertion element having replaced approximately 30 kb of DNA including the ica operon and histidine biosynthesis genes, resulting in histidine auxotrophy. This genomic rearrangement may be responsible for the diversification of ST71 into an emerging bovine adapted subgroup.

Is the whole greater than the sum of its parts? De novo assembly strategies for bacterial genomes based on paired-end sequencing

BACKGROUND

Whole genome sequence construction is becoming increasingly feasible because of advances in next generation sequencing (NGS), including increasing throughput and read length. By simply overlapping paired-end reads, we can obtain longer reads with higher accuracy, which can facilitate the assembly process. However, the influences of different library sizes and assembly methods on paired-end sequencing-based de novo assembly remain poorly understood.

RESULTS

We used 250 bp Illumina Miseq paired-end reads of different library sizes generated from genomic DNA from Escherichia coli DH1 and Streptococcus parasanguinis FW213 to compare the assembly results of different library sizes and assembly approaches. Our data indicate that overlapping paired-end reads can increase read accuracy but sometimes cause insertion or deletions. Regarding genome assembly, merged reads only outcompete original paired-end reads when coverage depth is low, and larger libraries tend to yield better assembly results. These results imply that distance information is the most critical factor during assembly. Our results also indicate that when depth is sufficiently high, assembly from subsets can sometimes produce better results.

CONCLUSIONS

In summary, this study provides systematic evaluations of de novo assembly from paired end sequencing data. Among the assembly strategies, we find that overlapping paired-end reads is not always beneficial for bacteria genome assembly and should be avoided or used with caution especially for genomes containing high fraction of repetitive sequences. Because increasing numbers of projects aim at bacteria genome sequencing, our study provides valuable suggestions for the field of genomic sequence construction.

Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing

BACKGROUND

Natural and artificial selection following domestication has led to the existence of more than a hundred pig breeds, as well as incredible variation in phenotypic traits. Berkshire pigs are regarded as having superior meat quality compared to other breeds. As the meat production industry seeks selective breeding approaches to improve profitable traits such as meat quality, information about genetic determinants of these traits is in high demand. However, most of the studies have been performed using trained sensory panel analysis without investigating the underlying genetic factors. Here we investigate the relationship between genomic composition and this phenotypic trait by scanning for signatures of positive selection in whole-genome sequencing data.

RESULTS

We generated genomes of 10 Berkshire pigs at a total of 100.6 coverage depth, using the Illumina Hiseq2000 platform. Along with the genomes of 11 Landrace and 13 Yorkshire pigs, we identified genomic variants of 18.9 million SNVs and 3.4 million Indels in the mapped regions. We identified several associated genes related to lipid metabolism, intramuscular fatty acid deposition, and muscle fiber type which attribute to pork quality (TG, FABP1, AKIRIN2, GLP2R, TGFBR3, JPH3, ICAM2, and ERN1) by applying between population statistical tests (XP-EHH and XP-CLR). A statistical enrichment test was also conducted to detect breed specific genetic variation. In addition, de novo short sequence read assembly strategy identified several candidate genes (SLC25A14, IGF1, PI4KA, CACNA1A) as also contributing to lipid metabolism.

CONCLUSIONS

Results revealed several candidate genes involved in Berkshire meat quality; most of these genes are involved in lipid metabolism and intramuscular fat deposition. These results can provide a basis for future research on the genomic characteristics of Berkshire pigs.

Next generation sequencing under de novo genome assembly

The next generation sequencing (NGS) is an important process which assures inexpensive organization of vast size of raw sequence dataset over any traditional sequencing systems or methods. Various aspects of NGS such as template preparation, sequencing imaging and genome alignment and assembly outline the genome sequencing and alignment. Consequently, de Bruijn graph (dBG) is an important mathematical tool that graphically analyzes how the orientations are constructed in groups of nucleotides. Basically, dBG describes the formation of the genome segments in circular iterative fashions. Some pivotal dBG-based de novo algorithms and software packages such as T-IDBA, Oases, IDBA-tran, Euler, Velvet, ABySS, AllPaths, SOAPde novo and SOAPde novo2 are illustrated in this paper. Consequently, overlap layout consensus (OLC) graph-based algorithms also play vital role in NGS assembly. Some important OLC-based algorithms such as MIRA3, CABOG, Newbler, Edena, Mosaik and SHORTY are portrayed in this paper. It has been experimented that greedy graph-based algorithms and software packages are also vital for proper genome dataset assembly. A few algorithms named SSAKE, SHARCGS and VCAKE help to perform proper genome sequencing.

Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases

Studies have demonstrated that ~60%–80% of emerging infectious diseases (EIDs) in humans originated from wild life. Bats are natural reservoirs of a large variety of viruses, including many important zoonotic viruses that cause severe diseases in humans and domestic animals. However, the understanding of the viral population and the ecological diversity residing in bat populations is unclear, which complicates the determination of the origins of certain EIDs. Here, using bats as a typical wildlife reservoir model, virome analysis was conducted based on pharyngeal and anal swab samples of 4440 bat individuals of 40 major bat species throughout China. The purpose of this study was to survey the ecological and biological diversities of viruses residing in these bat species, to investigate the presence of potential bat-borne zoonotic viruses and to evaluate the impacts of these viruses on public health. The data obtained in this study revealed an overview of the viral community present in these bat samples. Many novel bat viruses were reported for the first time and some bat viruses closely related to known human or animal pathogens were identified. This genetic evidence provides new clues in the search for the origin or evolution pattern of certain viruses, such as coronaviruses and noroviruses. These data offer meaningful ecological information for predicting and tracing wildlife-originated EIDs.

Characterization of the genome and transcriptome of the blue tit Cyanistes caeruleus: polymorphisms, sex-biased expression and selection signals

Decoding genomic sequences and determining their variation within populations has potential to reveal adaptive processes and unravel the genetic basis of ecologically relevant trait variation within a species. The blue tit Cyanistes caeruleus--a long-time ecological model species--has been used to investigate fitness consequences of variation in mating and reproductive behaviour. However, very little is known about the underlying genetic changes due to natural and sexual selection in the genome of this songbird. As a step to bridge this gap, we assembled the first draft genome of a single blue tit, mapped the transcriptome of five females and five males to this reference, identified genomewide variants and performed sex-differential expression analysis in the gonads, brain and other tissues. In the gonads, we found a high number of sex-biased genes, and of those, a similar proportion were sex-limited (genes only expressed in one sex) in males and females. However, in the brain, the proportion of female-limited genes within the female-biased gene category (82%) was substantially higher than the proportion of male-limited genes within the male-biased category (6%). This suggests a predominant on-off switching mechanism for the female-limited genes. In addition, most male-biased genes were located on the Z-chromosome, indicating incomplete dosage compensation for the male-biased genes. We called more than 500,000 SNPs from the RNA-seq data. Heterozygote detection in the single reference individual was highly congruent between DNA-seq and RNA-seq calling. Using information from these polymorphisms, we identified potential selection signals in the genome. We list candidate genes which can be used for further sequencing and detailed selection studies, including genes potentially related to meiotic drive evolution. A public genome browser of the blue tit with the described information is available at http://public-genomes-ngs.molgen.mpg.de.

Metabolic diversity and ecological niches of Achromatium populations revealed with single-cell genomic sequencing

Large, sulfur-cycling, calcite-precipitating bacteria in the genus Achromatium represent a significant proportion of bacterial communities near sediment-water interfaces at sites throughout the world. Our understanding of their potentially crucial roles in calcium, carbon, sulfur, nitrogen, and iron cycling is limited because they have not been cultured or sequenced using environmental genomics approaches to date. We utilized single-cell genomic sequencing to obtain one incomplete and two nearly complete draft genomes for Achromatium collected at Warm Mineral Springs (WMS), FL. Based on 16S rRNA gene sequences, the three cells represent distinct and relatively distant Achromatium populations (91-92% identity). The draft genomes encode key genes involved in sulfur and hydrogen oxidation; oxygen, nitrogen and polysulfide respiration; carbon and nitrogen fixation; organic carbon assimilation and storage; chemotaxis; twitching motility; antibiotic resistance; and membrane transport. Known genes for iron and manganese energy metabolism were not detected. The presence of pyrophosphatase and vacuolar (V)-type ATPases, which are generally rare in bacterial genomes, suggests a role for these enzymes in calcium transport, proton pumping, and/or energy generation in the membranes of calcite-containing inclusions.

GMcloser: closing gaps in assemblies accurately with a likelihood-based selection of contig or long-read alignments

MOTIVATION

Genome assemblies generated with next-generation sequencing (NGS) reads usually contain a number of gaps. Several tools have recently been developed to close the gaps in these assemblies with NGS reads. Although these gap-closing tools efficiently close the gaps, they entail a high rate of misassembly at gap-closing sites.

RESULTS

We have found that the assembly error rates caused by these tools are 20-500-fold higher than the rate of errors introduced into contigs by de novo assemblers. We here describe GMcloser, a tool that accurately closes these gaps with a preassembled contig set or a long read set (i.e., error-corrected PacBio reads). GMcloser uses likelihood-based classifiers calculated from the alignment statistics between scaffolds, contigs and paired-end reads to correctly assign contigs or long reads to gap regions of scaffolds, thereby achieving accurate and efficient gap closure. We demonstrate with sequencing data from various organisms that the gap-closing accuracy of GMcloser is 3-100-fold higher than those of other available tools, with similar efficiency.

AVAILABILITY AND IMPLEMENTATION

GMcloser and an accompanying tool (GMvalue) for evaluating the assembly and correcting misassemblies except SNPs and short indels in the assembly are available at https://sourceforge.net/projects/gmcloser/.

CONTACT

shunichi.kosugi@riken.jp.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.