Methods for generating shotgun and mixed shotgun/paired-end libraries for the 454 DNA sequencer

Metatranscriptomic analysis of a high-sulfide aquatic spring reveals insights into sulfur cycling and unexpected aerobic metabolism

Zodletone spring is a sulfide-rich spring in southwestern Oklahoma characterized by shallow, microoxic, light-exposed spring water overlaying anoxic sediments. Previously, culture-independent 16S rRNA gene based diversity surveys have revealed that Zodletone spring source sediments harbor a highly diverse microbial community, with multiple lineages putatively involved in various sulfur-cycling processes. Here, we conducted a metatranscriptomic survey of microbial populations in Zodletone spring source sediments to characterize the relative prevalence and importance of putative phototrophic, chemolithotrophic, and heterotrophic microorganisms in the sulfur cycle, the identity of lineages actively involved in various sulfur cycling processes, and the interaction between sulfur cycling and other geochemical processes at the spring source. Sediment samples at the spring’s source were taken at three different times within a 24-h period for geochemical analyses and RNA sequencing. In depth mining of datasets for sulfur cycling transcripts revealed major sulfur cycling pathways and taxa involved, including an unexpected potential role of Actinobacteria in sulfide oxidation and thiosulfate transformation. Surprisingly, transcripts coding for the cyanobacterial Photosystem II D1 protein, methane monooxygenase, and terminal cytochrome oxidases were encountered, indicating that genes for oxygen production and aerobic modes of metabolism are actively being transcribed, despite below-detectable levels (<1 µM) of oxygen in source sediment. Results highlight transcripts involved in sulfur, methane, and oxygen cycles, propose that oxygenic photosynthesis could support aerobic methane and sulfide oxidation in anoxic sediments exposed to sunlight, and provide a viewpoint of microbial metabolic lifestyles under conditions similar to those seen during late Archaean and Proterozoic eons.

Dynamics of an experimental microbial invasion

The ecological dynamics underlying species invasions have been a major focus of research in macroorganisms for the last five decades. However, we still know little about the processes behind invasion by unicellular organisms. To expand our knowledge of microbial invasions, we studied the roles of propagule pressure, nutrient supply, and biotic resistance in the invasion success of a freshwater invasive alga, Prymnesium parvum, using microcosms containing natural freshwater microbial assemblages. Microcosms were subjected to a factorial design with two levels of nutrient-induced diversity and three levels of propagule pressure, and incubated for 7 d, during which P. parvum densities and microbial community composition were tracked. Successful invasion occurred in microcosms receiving high propagule pressure whereas nutrients or community diversity played no role in invasion success. Invaded communities experienced distinctive changes in composition compared with communities where the invasion was unsuccessful. Successfully invaded microbial communities had an increased abundance of fungi and ciliates, and decreased abundances of diatoms and cercozoans. Many of these changes mirrored the microbial community changes detected during a natural P. parvum bloom in the source system. This role of propagule pressure is particularly relevant for P. parvum in the reservoir-dominated southern United States because this species can form large, sustained blooms that can generate intense propagule pressures for downstream sites. Human impact and global climate change are currently causing widespread environmental changes in most southern US freshwater systems that may facilitate P. parvum establishment and, when coupled with strong propagule pressure, could put many more systems at risk for invasion.

Coral-zooxanthellae meta-transcriptomics reveals integrated response to pollutant stress

BACKGROUND

Corals represent symbiotic meta-organisms that require harmonization among the coral animal, photosynthetic zooxanthellae and associated microbes to survive environmental stresses. We investigated integrated-responses among coral and zooxanthellae in the scleractinian coral Acropora formosa in response to an emerging marine pollutant, the munitions constituent, 1,3,5-trinitro-1,3,5 triazine (RDX; 5 day exposures to 0 (control), 0.5, 0.9, 1.8, 3.7, and 7.2 mg/L, measured in seawater).

RESULTS

RDX accumulated readily in coral soft tissues with bioconcentration factors ranging from 1.1 to 1.5. Next-generation sequencing of a normalized meta-transcriptomic library developed for the eukaryotic components of the A. formosa coral holobiont was leveraged to conduct microarray-based global transcript expression analysis of integrated coral/zooxanthellae responses to the RDX exposure. Total differentially expressed transcripts (DET) increased with increasing RDX exposure concentrations as did the proportion of zooxanthellae DET relative to the coral animal. Transcriptional responses in the coral demonstrated higher sensitivity to RDX compared to zooxanthellae where increased expression of gene transcripts coding xenobiotic detoxification mechanisms (i.e. cytochrome P450 and UDP glucuronosyltransferase 2 family) were initiated at the lowest exposure concentration. Increased expression of these detoxification mechanisms was sustained at higher RDX concentrations as well as production of a physical barrier to exposure through a 40% increase in mucocyte density at the maximum RDX exposure. At and above the 1.8 mg/L exposure concentration, DET coding for genes involved in central energy metabolism, including photosynthesis, glycolysis and electron-transport functions, were decreased in zooxanthellae although preliminary data indicated that zooxanthellae densities were not affected. In contrast, significantly increased transcript expression for genes involved in cellular energy production including glycolysis and electron-transport pathways was observed in the coral animal.

CONCLUSIONS

Transcriptional network analysis for central energy metabolism demonstrated highly correlated responses to RDX among the coral animal and zooxanthellae indicative of potential compensatory responses to lost photosynthetic potential within the holobiont. These observations underscore the potential for complex integrated responses to RDX exposure among species comprising the coral holobiont and highlight the need to understand holobiont-species interactions to accurately assess pollutant impacts.

BAC-pool sequencing and analysis of large segments of A12 and D12 homoeologous chromosomes in upland cotton

Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.

Sequencing approach to analyze the role of quasispecies for classical swine fever

Classical swine fever virus (CSFV) is a positive-sense RNA virus with a high degree of genetic variability among isolates. High diversity is also found in virulence, with strains covering the complete spectrum from avirulent to highly virulent. The underlying genetic determinants are far from being understood. Since RNA polymerases of RNA viruses lack any proof-reading activity, different genome variations called haplotypes, occur during replication. A set of haplotypes is referred to as a viral quasispecies. Genetic variability can be a fitness advantage through facilitating of a more effective escape from the host immune response. In order to investigate the correlation of quasispecies composition and virulence in vivo, we analyzed next-generation sequencing data of CSFV isolates of varying virulence. Viral samples from pigs infected with the highly virulent isolates "Koslov" and "Brescia" showed higher quasispecies diversity and more nucleotide variability, compared to samples of pigs infected with low and moderately virulent isolates.

Gene Assembly from Chip-Synthesized Oligonucleotides

De novo synthesis of long double-stranded DNA constructs has a myriad of applications in biology and biological engineering. However, its widespread adoption has been hindered by high costs. Cost can be significantly reduced by using oligonucleotides synthesized on high-density DNA chips. However, most methods for using off-chip DNA for gene synthesis have failed to scale due to the high error rates, low yields, and high chemical complexity of the chip-synthesized oligonucleotides. We have recently demonstrated that some commercial DNA chip manufacturers have improved error rates, and that the issues of chemical complexity and low yields can be solved by using barcoded primers to accurately and efficiently amplify subpools of oligonucleotides. This article includes protocols for computationally designing the DNA chip, amplifying the oligonucleotide subpools, and assembling 500-800 basepair (bp) constructs.

Impact of host cell line adaptation on quasispecies composition and glycosylation of influenza A virus hemagglutinin

The genome of influenza A viruses is constantly changing (genetic drift) resulting in small, gradual changes in viral proteins. Alterations within antibody recognition sites of the viral membrane glycoproteins hemagglutinin (HA) and neuraminidase (NA) result in an antigenetic drift, which requires the seasonal update of human influenza virus vaccines. Generally, virus adaptation is necessary to obtain sufficiently high virus yields in cell culture-derived vaccine manufacturing. In this study detailed HA N-glycosylation pattern analysis was combined with in-depth pyrosequencing analysis of the virus genomic RNA. Forward and backward adaptation from Madin-Darby Canine Kidney (MDCK) cells to African green monkey kidney (Vero) cells was investigated for two closely related influenza A virus PR/8/34 (H1N1) strains: from the National Institute for Biological Standards and Control (NIBSC) or the Robert Koch Institute (RKI). Furthermore, stability of HA N-glycosylation patterns over ten consecutive passages and different harvest time points is demonstrated. Adaptation to Vero cells finally allowed efficient influenza A virus replication in Vero cells. In contrast, during back-adaptation the virus replicated well from the very beginning. HA N-glycosylation patterns were cell line dependent and stabilized fast within one (NIBSC-derived virus) or two (RKI-derived virus) successive passages during adaptation processes. However, during adaptation new virus variants were detected. These variants carried "rescue" mutations on the genomic level within the HA stem region, which result in amino acid substitutions. These substitutions finally allowed sufficient virus replication in the new host system. According to adaptation pressure the composition of the virus populations varied. In Vero cells a selection for "rescue" variants was characteristic. After back-adaptation to MDCK cells some variants persisted at indifferent frequencies, others slowly diminished and even dropped below the detection limit.

Highly pathogenic avian influenza virus subtype H5N1 escaping neutralization: more than HA variation

Influenza A viruses are one of the major threats in modern health care. Novel viruses arise due to antigenic drift and antigenic shift, leading to escape from the immune system and resulting in a serious problem for disease control. In order to investigate the escape process and to enable predictions of escape, we serially passaged influenza A H5N1 virus in vitro 100 times under immune pressure. The generated escape viruses were characterized phenotypically and in detail by full-genome deep sequencing. Mutations already found in natural isolates were detected, evidencing the in vivo relevance of the in vitro-induced amino acid substitutions. Additionally, several novel alterations were triggered. Altogether, the results imply that our in vitro system is suitable to study influenza A virus evolution and that it might even be possible to predict antigenic changes of influenza A viruses circulating in vaccinated populations.

Complete genome sequence of virulent duck enteritis virus (DEV) strain 2085 and comparison with genome sequences of virulent and attenuated DEV strains

We here report the complete genome sequence of the duck enteritis virus (DEV) wild-type strain 2085, an avian herpesvirus (GenBank ID: JF999965). The nucleotide sequence was derived from the 2085 genome cloned as an infectious bacterial artificial chromosome (BAC) clone. The DEV 2085 genome is 160,649-bp in length and encodes 78 predicted open reading frames (ORFs), a number identical to that identified for the attenuated DEV VAC strain (GenBank ID: EU082088.2). Comparison of the genome sequences DEV 2085 and VAC with partial sequences of the virulent CHv strain and the attenuated strain Clone-03 was carried out to identify nucleotide or amino acid polymorphisms that potentially contribute to DEV virulence. No amino acid changes were identified in 24 of the 78 ORFs, a result indicating high conservation in DEV independently of strain origin or virulence. In addition, 39 ORFs contain non-synonymous nucleotide substitutions, while 15 ORFs had nucleotide insertions or deletions, frame-shift mutations and/or non-synonymous nucleotide substitutions with an effect on ORF initiation or termination. In 7 of the 15 ORFs with high and 27 of the 39 ORFs with low variability, polymorphisms were exclusively found in DEV 2085, a finding that likely is a result of a different origin of 2085 (Europe) or VAC, Clone-03 and CHv (Eastern Asia). Five ORFs (UL2, UL12, US10, UL47 and UL41) with polymorphisms were identical between the virulent DEV 2085 and CHv but different from VAC or Clone-03. They, individually or in combination, may therefore represent DEV virulence factors. Our comparative analysis of four DEV sequences provides a comprehensive overview of DEV genome structure and identifies ORFs that are changed during serial virus passage.

A comprehensive deep sequencing strategy for full-length genomes of influenza A

Driven by the impact of influenza A viruses on human and animal health, much research is conducted on this pathogen. To support this research, we designed an all influenza A-embracing reverse transcription-PCR (RT-PCR) for the generation of DNA from influenza A virus negative strand RNA genome segments for full-length genome deep sequencing on a Genome Sequencer FLX instrument. For high reliability, the RT-PCRs are designed such that every genome segment is divided into two amplicons and for the most variable segments redundancy is included. Moreover, to minimize the risk of contamination of diagnostic real-time PCRs by sequencing amplicons, RT-PCR does not generate amplicons that are amenable to RT-qPCR detection. With the presented protocol we were able to generate virtually all amplicons (99.3% success rate) from isolates representing all so far known 16 hemagglutinin and 9 neuraminidase subtypes and from an additional 2009 pandemic influenza A H1N1 virus. Three isolates were sequenced to analyze the suitability of the DNA for sequencing. Moreover, we provide a short R script that disambiguates the sequences of the primers used. We show that using unambiguous primer sequences for read trimming prior to assembly with the genome sequencer assembler software results in higher quality of the final genome sequences. Using the disambiguated primer sequences, high quality full-length sequences for the three isolates used for sequencing trials could be established from the raw data in de novo assemblies.

Genotype and phenotypes of an intestine-adapted Escherichia coli K-12 mutant selected by animal passage for superior colonization

We previously isolated a spontaneous mutant of Escherichia coli K-12, strain MG1655, following passage through the streptomycin-treated mouse intestine, that has colonization traits superior to the wild-type parent strain (M. P. Leatham et al., Infect. Immun. 73:8039-8049, 2005). This intestine-adapted strain (E. coli MG1655*) grew faster on several different carbon sources than the wild type and was nonmotile due to deletion of the flhD gene. We now report the results of several high-throughput genomic analysis approaches to further characterize E. coli MG1655*. Whole-genome pyrosequencing did not reveal any changes on its genome, aside from the deletion at the flhDC locus, that could explain the colonization advantage of E. coli MG1655*. Microarray analysis revealed modest yet significant induction of catabolic gene systems across the genome in both E. coli MG1655* and an isogenic flhD mutant constructed in the laboratory. Catabolome analysis with Biolog GN2 microplates revealed an enhanced ability of both E. coli MG1655* and the isogenic flhD mutant to oxidize a variety of carbon sources. The results show that intestine-adapted E. coli MG1655* is more fit than the wild type for intestinal colonization, because loss of FlhD results in elevated expression of genes involved in carbon and energy metabolism, resulting in more efficient carbon source utilization and a higher intestinal population. Hence, mutations that enhance metabolic efficiency confer a colonization advantage.

Recovery of infectious virus from full-length cowpox virus (CPXV) DNA cloned as a bacterial artificial chromosome (BAC)

Transmission from pet rats and cats to humans as well as severe infection in felids and other animal species have recently drawn increasing attention to cowpox virus (CPXV). We report the cloning of the entire genome of cowpox virus strain Brighton Red (BR) as a bacterial artificial chromosome (BAC) in Escherichia coli and the recovery of infectious virus from cloned DNA. Generation of a full-length CPXV DNA clone was achieved by first introducing a mini-F vector, which allows maintenance of large circular DNA in E. coli, into the thymidine kinase locus of CPXV by homologous recombination. Circular replication intermediates were then electroporated into E. coli DH10B cells. Upon successful establishment of the infectious BR clone, we modified the full-length clone such that recombination-mediated excision of bacterial sequences can occur upon transfection in eukaryotic cells. This self-excision of the bacterial replicon is made possible by a sequence duplication within mini-F sequences and allows recovery of recombinant virus progeny without remaining marker or vector sequences. The in vitro growth properties of viruses derived from both BAC clones were determined and found to be virtually indistinguishable from those of parental, wild-type BR. Finally, the complete genomic sequence of the infectious clone was determined and the cloned viral genome was shown to be identical to that of the parental virus. In summary, the generated infectious clone will greatly facilitate studies on individual genes and pathogenesis of CPXV. Moreover, the vector potential of CPXV can now be more systematically explored using this newly generated tool.

Titration-free massively parallel pyrosequencing using trace amounts of starting material

Continuous efforts have been made to improve next-generation sequencing methods for increased robustness and for applications on low amounts of starting material. We applied double-stranded library protocols for the Roche 454 platform to avoid the yield-reducing steps associated with single-stranded library preparation, and applied a highly sensitive Taqman MGB-probe-based quantitative polymerase chain reaction (qPCR) method. The MGB-probe qPCR, which can detect as low as 100 copies, was used to quantify the amount of effective library, i.e. molecules that form functional clones in emulsion PCR. We also demonstrate that the distribution of library molecules on capture beads follows a Poisson distribution. Combining the qPCR and Poisson statistics, the labour-intensive and costly titration can be eliminated and trace amounts of starting material such as precious clinical samples, transcriptomes of small tissue samples and metagenomics on low biomass environments is applicable.

Increased throughput by parallelization of library preparation for massive sequencing

Background

Massively parallel sequencing systems continue to improve on data output, while leaving labor-intensive library preparations a potential bottleneck. Efforts are currently under way to relieve the crucial and time-consuming work to prepare DNA for high-throughput sequencing.

Methodology/Principal Findings

In this study, we demonstrate an automated parallel library preparation protocol using generic carboxylic acid-coated superparamagnetic beads and polyethylene glycol precipitation as a reproducible and flexible method for DNA fragment length separation. With this approach the library preparation for DNA sequencing can easily be adjusted to a desired fragment length. The automated protocol, here demonstrated using the GS FLX Titanium instrument, was compared to the standard manual library preparation, showing higher yield, throughput and great reproducibility. In addition, 12 libraries were prepared and uniquely tagged in parallel, and the distribution of sequence reads between these indexed samples could be improved using quantitative PCR-assisted pooling.

Conclusions/Significance

We present a novel automated procedure that makes it possible to prepare 36 indexed libraries per person and day, which can be increased to up to 96 libraries processed simultaneously. The yield, speed and robust performance of the protocol constitute a substantial improvement to present manual methods, without the need of extensive equipment investments. The described procedure enables a considerable efficiency increase for small to midsize sequencing centers.

Removal of deaminated cytosines and detection of in vivo methylation in ancient DNA

DNA sequences determined from ancient organisms have high error rates, primarily due to uracil bases created by cytosine deamination. We use synthetic oligonucleotides, as well as DNA extracted from mammoth and Neandertal remains, to show that treatment with uracil-DNA-glycosylase and endonuclease VIII removes uracil residues from ancient DNA and repairs most of the resulting abasic sites, leaving undamaged parts of the DNA fragments intact. Neandertal DNA sequences determined with this protocol have greatly increased accuracy. In addition, our results demonstrate that Neandertal DNA retains in vivo patterns of CpG methylation, potentially allowing future studies of gene inactivation and imprinting in ancient organisms.

Phylotyping and functional analysis of two ancient human microbiomes

BACKGROUND

The Human Microbiome Project (HMP) is one of the U.S. National Institutes of Health Roadmap for Medical Research. Primary interests of the HMP include the distinctiveness of different gut microbiomes, the factors influencing microbiome diversity, and the functional redundancies of the members of human microbiotas. In this present work, we contribute to these interests by characterizing two extinct human microbiotas.

METHODOLOGY/PRINCIPAL FINDINGS

We examine two paleofecal samples originating from cave deposits in Durango Mexico and dating to approximately 1300 years ago. Contamination control is a serious issue in ancient DNA research; we use a novel approach to control contamination. After we determined that each sample originated from a different human, we generated 45 thousand shotgun DNA sequencing reads. The phylotyping and functional analysis of these reads reveals a signature consistent with the modern gut ecology. Interestingly, inter-individual variability for phenotypes but not functional pathways was observed. The two ancient samples have more similar functional profiles to each other than to a recently published profile for modern humans. This similarity could not be explained by a chance sampling of the databases.

CONCLUSIONS/SIGNIFICANCE

We conduct a phylotyping and functional analysis of ancient human microbiomes, while providing novel methods to control for DNA contamination and novel hypotheses about past microbiome biogeography. We postulate that natural selection has more of an influence on microbiome functional profiles than it does on the species represented in the microbial ecology. We propose that human microbiomes were more geographically structured during pre-Columbian times than today.