A multiple-alignment based primer design algorithm for genetically highly variable DNA targets.

BACKGROUND: Primer design for highly variable DNA sequences is difficult, and experimental success requires attention to many interacting constraints. The advent of next-generation sequencing methods allows the investigation of rare variants otherwise hidden deep in large populations, but requires attention to population diversity and primer localization in relatively conserved regions, in addition to recognized constraints typically considered in primer design. RESULTS: Design constraints include degenerate sites to maximize population coverage, matching of melting temperatures, optimizing de novo sequence length, finding optimal bio-barcodes to allow efficient downstream analyses, and minimizing risk of dimerization. To facilitate primer design addressing these and other constraints, we created a novel computer program (PrimerDesign) that automates this complex procedure. We show its powers and limitations and give examples of successful designs for the analysis of HIV-1 populations. CONCLUSIONS: PrimerDesign is useful for researchers who want to design DNA primers and probes for analyzing highly variable DNA populations. It can be used to design primers for PCR, RT-PCR, Sanger sequencing, next-generation sequencing, and other experimental protocols targeting highly variable DNA samples.

Epitope-specific CD8+ T lymphocytes cross-recognize mutant simian immunodeficiency virus (SIV) sequences but fail to contain very early evolution and eventual fixation of epitope escape mutations during SIV infection.

Human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) evade containment by CD8(+) T lymphocytes through focused epitope mutations. However, because of limitations in the numbers of viral sequences that can be sampled, traditional sequencing technologies have not provided a true representation of the plasticity of these viruses or the intensity of CD8(+) T lymphocyte-mediated selection pressure. Moreover, the strategy by which CD8(+) T lymphocytes contain evolving viral quasispecies has not been characterized fully. In the present study we have employed ultradeep 454 pyrosequencing of virus and simultaneous staining of CD8(+) T lymphocytes with multiple tetramers in the SIV/rhesus monkey model to explore the coevolution of virus and the cellular immune response during primary infection. We demonstrated that cytotoxic T lymphocyte (CTL)-mediated selection pressure on the infecting virus was manifested by epitope mutations as early as 21 days following infection. We also showed that CD8(+) T lymphocytes cross-recognized wild-type and mutant epitopes and that these cross-reactive cell populations were present at a time when mutant forms of virus were present at frequencies of as low as 1 in 22,000 sequenced clones. Surprisingly, these cross-reactive cells became enriched in the epitope-specific CD8(+) T lymphocyte population as viruses with mutant epitope sequences largely replaced those with epitope sequences of the transmitted virus. These studies demonstrate that mutant epitope-specific CD8(+) T lymphocytes that are present at a time when viral mutant epitope sequences are detected at extremely low frequencies fail to contain the later accumulation and fixation of the mutant epitope sequences in the viral quasispecies.

Non-contiguous finished genome sequence of plant-growth promoting Serratia proteamaculans S4.

Serratia proteamaculans S4 (previously Serratia sp. S4), isolated from the rhizosphere of wild Equisetum sp., has the ability to stimulate plant growth and to suppress the growth of several soil-borne fungal pathogens of economically important crops. Here we present the non-contiguous, finished genome sequence of S. proteamaculans S4, which consists of a 5,324,944 bp circular chromosome and a 129,797 bp circular plasmid. The chromosome contains 5,008 predicted genes while the plasmid comprises 134 predicted genes. In total, 4,993 genes are assigned as protein-coding genes. The genome consists of 22 rRNA genes, 82 tRNA genes and 58 pseudogenes. This genome is a part of the project "Genomics of four rapeseed plant growth-promoting bacteria with antagonistic effect on plant pathogens" awarded through the 2010 DOE-JGI's Community Sequencing Program.

Complete genome sequence of Enterobacter sp. IIT-BT 08: A potential microbial strain for high rate hydrogen production.

Enterobacter sp. IIT-BT 08 belongs to Phylum: Proteobacteria, Class: Gammaproteobacteria, Order: Enterobacteriales, Family: Enterobacteriaceae. The organism was isolated from the leaves of a local plant near the Kharagpur railway station, Kharagpur, West Bengal, India. It has been extensively studied for fermentative hydrogen production because of its high hydrogen yield. For further enhancement of hydrogen production by strain development, complete genome sequence analysis was carried out. Sequence analysis revealed that the genome was linear, 4.67 Mbp long and had a GC content of 56.01%. The genome properties encode 4,393 protein-coding and 179 RNA genes. Additionally, a putative pathway of hydrogen production was suggested based on the presence of formate hydrogen lyase complex and other related genes identified in the genome. Thus, in the present study we describe the specific properties of the organism and the generation, annotation and analysis of its genome sequence as well as discuss the putative pathway of hydrogen production by this organism.

Improving genome assemblies by sequencing PCR products with PacBio.

Advances in sequencing technologies have dramatically reduced costs in producing high-quality draft genomes. However, there are still many contigs and possible misassembled regions in those draft genomes. Improving the quality of these genomes requires an efficient and economical means to close gaps and resequence some regions. Sequencing pooled gap region PCR products with Pacific Biosciences (PacBio) provides a significantly less expensive means for this need. We have developed a genome improvement pipeline with this strategy after decreasing a loading bias against larger PCR products in the PacBio process. Compared with Sanger technology, this approach is not only cost-effective but also can close gaps greater than 2.5 kb in a single round of reactions, and sequence through high GC regions as well as difficult secondary structures such as small hairpin loops.

Artificial polyploidy improves bacterial single cell genome recovery.

BACKGROUND: Single cell genomics (SCG) is a combination of methods whose goal is to decipher the complete genomic sequence from a single cell and has been applied mostly to organisms with smaller genomes, such as bacteria and archaea. Prior single cell studies showed that a significant portion of a genome could be obtained. However, breakages of genomic DNA and amplification bias have made it very challenging to acquire a complete genome with single cells. We investigated an artificial method to induce polyploidy in Bacillus subtilis ATCC 6633 by blocking cell division and have shown that we can significantly improve the performance of genomic sequencing from a single cell. METHODOLOGY/PRINCIPAL FINDINGS: We inhibited the bacterial cytoskeleton protein FtsZ in B.subtilis with an FtsZ-inhibiting compound, PC190723, resulting in larger undivided single cells with multiple copies of its genome. qPCR assays of these larger, sorted cells showed higher DNA content, have less amplification bias, and greater genomic recovery than untreated cells. SIGNIFICANCE: The method presented here shows the potential to obtain a nearly complete genome sequence from a single bacterial cell. With millions of uncultured bacterial species in nature, this method holds tremendous promise to provide insight into the genomic novelty of yet-to-be discovered species, and given the temporary effects of artificial polyploidy coupled with the ability to sort and distinguish differences in cell size and genomic DNA content, may allow recovery of specific organisms in addition to their genomes.

Temperature regulation of virulence factors in the pathogen Vibrio coralliilyticus.

Sea surface temperatures (SST) are rising because of global climate change. As a result, pathogenic Vibrio species that infect humans and marine organisms during warmer summer months are of growing concern. Coral reefs, in particular, are already experiencing unprecedented degradation worldwide due in part to infectious disease outbreaks and bleaching episodes that are exacerbated by increasing SST. For example, Vibrio coralliilyticus, a globally distributed bacterium associated with multiple coral diseases, infects corals at temperatures above 27 °C. The mechanisms underlying this temperature-dependent pathogenicity, however, are unknown. In this study, we identify potential virulence mechanisms using whole genome sequencing of V. coralliilyticus ATCC (American Type Culture Collection) BAA-450. Furthermore, we demonstrate direct temperature regulation of numerous virulence factors using proteomic analysis and bioassays. Virulence factors involved in motility, host degradation, secretion, antimicrobial resistance and transcriptional regulation are upregulated at the higher virulent temperature of 27 °C, concurrent with phenotypic changes in motility, antibiotic resistance, hemolysis, cytotoxicity and bioluminescence. These results provide evidence that temperature regulates multiple virulence mechanisms in V. coralliilyticus, independent of abundance. The ecological and biological significance of this temperature-dependent virulence response is reinforced by climate change models that predict tropical SST to consistently exceed 27 °C during the spring, summer and fall seasons. We propose V. coralliilyticus as a model Gram-negative bacterium to study temperature-dependent pathogenicity in Vibrio-related diseases.

Complete genome sequence of Dehalogenimonas lykanthroporepellens type strain (BL-DC-9(T)) and comparison to "Dehalococcoides" strains.

Dehalogenimonas lykanthroporepellens is the type species of the genus Dehalogenimonas, which belongs to a deeply branching lineage within the phylum Chloroflexi. This strictly anaerobic, mesophilic, non spore-forming, Gram-negative staining bacterium was first isolated from chlorinated solvent contaminated groundwater at a Superfund site located near Baton Rouge, Louisiana, USA. D. lykanthroporepellens was of interest for genome sequencing for two reasons: (a) an unusual ability to couple growth with reductive dechlorination of environmentally important polychlorinated aliphatic alkanes and (b) a phylogenetic position that is distant from previously sequenced bacteria. The 1,686,510 bp circular chromosome of strain BL-DC-9(T) contains 1,720 predicted protein coding genes, 47 tRNA genes, a single large subunit rRNA (23S-5S) locus, and a single, orphan, small subunit rRNA (16S) locus.

Genomics for key players in the N cycle from guinea pigs to the next frontier.

While sequencing methods were available in the late 1970s, it was not until the human genome project and a significant influx of funds for such research that this technology became high throughput. The fields of microbiology and microbial ecology, among many others, have been tremendously impacted over the years, to such an extent that the determination of complete microbial genome sequences is now commonplace. Given the lower costs of next-generation sequencing platforms, even small laboratories from around the world will be able to generate millions of base pairs of data, equivalent to entire genomes worth of sequence information. With this prospect just around the corner, it is timely to provide an overview of the genomics process: from sample preparation to some of the analytical methods used to gain functional knowledge from sequence information.

A rapid multiplex assay for nucleic acid-based diagnostics.

We have developed a rapid (under 4 hours), multiplex, nucleic acid assay, adapted to a microsphere array detection platform. We call this assay multiplex oligonucleotide ligation-PCR (MOL-PCR). Unlike other ligation-based assays that require multiple steps, our protocol consists of a single tube reaction, followed by hybridization to a Luminex microsphere array for detection. We demonstrate the ability of this assay to simultaneously detect diverse nucleic acid signatures (e.g., unique sequences, single nucleotide polymorphisms) in a single multiplex reaction. Detection probes consist of modular components that enable target detection, probe amplification, and subsequent capture onto microsphere arrays. To demonstrate the utility of our assay, we applied it to the detection of three biothreat agents, B. anthracis, Y. pestis, and F. tularensis. Combined with the ease and robustness of this assay, the results presented here show a strong potential of our assay for use in diagnostics and surveillance.

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

Mechanism of induction and suppression of antiviral immunity directed by virus-derived small RNAs in Drosophila.

The small RNA-directed viral immunity pathway in plants and invertebrates begins with the production by Dicer nuclease of virus-derived siRNAs (viRNAs), which guide specific antiviral silencing by Argonaute protein in an RNA-induced silencing complex (RISC). Molecular identity of the viral RNA precursor of viRNAs remains a matter of debate. Using Flock house virus (FHV) infection of Drosophila as a model, we show that replication of FHV positive-strand RNA genome produces an approximately 400 bp dsRNA from its 5' terminus that serves as the major Dicer-2 substrate. ViRNAs thus generated are loaded in Argonaute-2 and methylated at their 3' ends. Notably, FHV-encoded RNAi suppressor B2 protein interacts with both viral dsRNA and RNA replicase and inhibits production of the 5'-terminal viRNAs. Our findings, therefore, provide a model in which small RNA-directed viral immunity is induced during the initiation of viral progeny (+)RNA synthesis and suppressed by B2 inside the viral RNA replication complex.

The complete genome sequence of Bacillus thuringiensis Al Hakam.

Bacillus thuringiensis is an insect pathogen that is widely used as a biopesticide (E. Schnepf, N. Crickmore, J. Van Rie, D. Lereclus, J. Baum, J. Feitelson, D. R. Zeigler, and D. H. Dean, Microbiol. Mol. Biol. Rev. 62:775-806, 1998). Here we report the finished, annotated genome sequence of B. thuringiensis Al Hakam, which was collected in Iraq by the United Nations Special Commission (L. Radnedge, P. Agron, K. Hill, P. Jackson, L. Ticknor, P. Keim, and G. Andersen, Appl. Environ. Microbiol. 69:2755-2764, 2003).

Pathogenomic sequence analysis of Bacillus cereus and Bacillus thuringiensis isolates closely related to Bacillus anthracis.

Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis are closely related gram-positive, spore-forming bacteria of the B. cereus sensu lato group. While independently derived strains of B. anthracis reveal conspicuous sequence homogeneity, environmental isolates of B. cereus and B. thuringiensis exhibit extensive genetic diversity. Here we report the sequencing and comparative analysis of the genomes of two members of the B. cereus group, B. thuringiensis 97-27 subsp. konkukian serotype H34, isolated from a necrotic human wound, and B. cereus E33L, which was isolated from a swab of a zebra carcass in Namibia. These two strains, when analyzed by amplified fragment length polymorphism within a collection of over 300 of B. cereus, B. thuringiensis, and B. anthracis isolates, appear closely related to B. anthracis. The B. cereus E33L isolate appears to be the nearest relative to B. anthracis identified thus far. Whole-genome sequencing of B. thuringiensis 97-27and B. cereus E33L was undertaken to identify shared and unique genes among these isolates in comparison to the genomes of pathogenic strains B. anthracis Ames and B. cereus G9241 and nonpathogenic strains B. cereus ATCC 10987 and B. cereus ATCC 14579. Comparison of these genomes revealed differences in terms of virulence, metabolic competence, structural components, and regulatory mechanisms.