Matthew R. Henn

Naturally Occurring Dominant Resistance Mutations to Hepatitis C Virus Protease and Polymerase Inhibitors in Treatment-Naïve Patients

Resistance mutations to hepatitis C virus (HCV) nonstructural protein 3 (NS3) protease inhibitors in <1% of the viral quasispecies may still allow >1000‐fold viral load reductions upon treatment, consistent with their reported reduced replicative fitness in vitro. Recently, however, an R155K protease mutation was reported as the dominant quasispecies in a treatment‐naïve individual, raising concerns about possible full drug resistance. To investigate the prevalence of dominant resistance mutations against specifically targeted antiviral therapy for HCV (STAT‐C) in the population, we analyzed HCV genome sequences from 507 treatment‐naïve patients infected with HCV genotype 1 from the United States, Germany, and Switzerland. Phylogenetic sequence analysis and viral load data were used to identify the possible spread of replication‐competent, drug‐resistant viral strains in the population and to infer the consequences of these mutations upon viral replication in vivo. Mutations described to confer resistance to the protease inhibitors Telaprevir, BILN2061, ITMN‐191, SCH6 and Boceprevir; the NS5B polymerase inhibitor AG‐021541; and to the NS4A antagonist ACH‐806 were observed mostly as sporadic, unrelated cases, at frequencies between 0.3% and 2.8% in the population, including two patients with possible multidrug resistance. Collectively, however, 8.6% of the patients infected with genotype 1a and 1.4% of those infected with genotype 1b carried at least one dominant resistance mutation. Viral loads were high in the majority of these patients, suggesting that drug‐resistant viral strains might achieve replication levels comparable to nonresistant viruses in vivo. Conclusion: Naturally occurring dominant STAT‐C resistance mutations are common in treatment‐naïve patients infected with HCV genotype 1. Their influence on treatment outcome should further be characterized to evaluate possible benefits of drug resistance testing for individual tailoring of drug combinations when treatment options are limited due to previous nonresponse to peginterferon and ribavirin. (HEPATOLOGY 2008;48:1769–1778.)

Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection

Deep sequencing technologies have the potential to transform the study of highly variable viral pathogens by providing a rapid and cost-effective approach to sensitively characterize rapidly evolving viral quasispecies. Here, we report on a high-throughput whole HIV-1 genome deep sequencing platform that combines 454 pyrosequencing with novel assembly and variant detection algorithms. In one subject we combined these genetic data with detailed immunological analyses to comprehensively evaluate viral evolution and immune escape during the acute phase of HIV-1 infection. The majority of early, low frequency mutations represented viral adaptation to host CD8+ T cell responses, evidence of strong immune selection pressure occurring during the early decline from peak viremia. CD8+ T cell responses capable of recognizing these low frequency escape variants coincided with the selection and evolution of more effective secondary HLA-anchor escape mutations. Frequent, and in some cases rapid, reversion of transmitted mutations was also observed across the viral genome. When located within restricted CD8 epitopes these low frequency reverting mutations were sufficient to prime de novo responses to these epitopes, again illustrating the capacity of the immune response to recognize and respond to low frequency variants. More importantly, rapid viral escape from the most immunodominant CD8+ T cell responses coincided with plateauing of the initial viral load decline in this subject, suggestive of a potential link between maintenance of effective, dominant CD8 responses and the degree of early viremia reduction. We conclude that the early control of HIV-1 replication by immunodominant CD8+ T cell responses may be substantially influenced by rapid, low frequency viral adaptations not detected by conventional sequencing approaches, which warrants further investigation. These data support the critical need for vaccine-induced CD8+ T cell responses to target more highly constrained regions of the virus in order to ensure the maintenance of immunodominant CD8 responses and the sustained decline of early viremia.

Whole Genome Amplification and De novo Assembly of Single Bacterial Cells

Background Single-cell genome sequencing has the potential to allow the in-depth exploration of the vast genetic diversity found in uncultured microbes. We used the marine cyanobacterium Prochlorococcus as a model system for addressing important challenges facing high-throughput whole genome amplification (WGA) and complete genome sequencing of individual cells. Methodology/Principal Findings We describe a pipeline that enables single-cell WGA on hundreds of cells at a time while virtually eliminating non-target DNA from the reactions. We further developed a post-amplification normalization procedure that mitigates extreme variations in sequencing coverage associated with multiple displacement amplification (MDA), and demonstrated that the procedure increased sequencing efficiency and facilitated genome assembly. We report genome recovery as high as 99.6% with reference-guided assembly, and 95% with de novo assembly starting from a single cell. We also analyzed the impact of chimera formation during MDA on de novo assembly, and discuss strategies to minimize the presence of incorrectly joined regions in contigs. Conclusions/Significance The methods describe in this paper will be useful for sequencing genomes of individual cells from a variety of samples.

Genomic analysis of oceanic cyanobacterial myoviruses compared with T4-like myoviruses from diverse hosts and environments

T4-like myoviruses are ubiquitous, and their genes are among the most abundant documented in ocean systems. Here we compare 26 T4-like genomes, including 10 from non-cyanobacterial myoviruses, and 16 from marine cyanobacterial myoviruses (cyanophages) isolated on diverse Prochlorococcus or Synechococcus hosts. A core genome of 38 virion construction and DNA replication genes was observed in all 26 genomes, with 32 and 25 additional genes shared among the non-cyanophage and cyanophage subsets, respectively. These hierarchical cores are highly syntenic across the genomes, and sampled to saturation. The 25 cyanophage core genes include six previously described genes with putative functions (psbA, mazG, phoH, hsp20, hli03, cobS), a hypothetical protein with a potential phytanoyl-CoA dioxygenase domain, two virion structural genes, and 16 hypothetical genes. Beyond previously described cyanophage-encoded photosynthesis and phosphate stress genes, we observed core genes that may play a role in nitrogen metabolism during infection through modulation of 2-oxoglutarate. Patterns among non-core genes that may drive niche diversification revealed that phosphorus-related gene content reflects source waters rather than host strain used for isolation, and that carbon metabolism genes appear associated with putative mobile elements. As well, phages isolated on Synechococcus had higher genome-wide %G+C and often contained different gene subsets (e.g. petE, zwf, gnd, prnA, cpeT) than those isolated on Prochlorococcus. However, no clear diagnostic genes emerged to distinguish these phage groups, suggesting blurred boundaries possibly due to cross-infection. Finally, genome-wide comparisons of both diverse and closely related, co-isolated genomes provide a locus-to-locus variability metric that will prove valuable for interpreting metagenomic data sets.

Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives.

While most Ascomycetes tend to associate principally with plants, the dimorphic fungi Coccidioides immitis and Coccidioides posadasii are primary pathogens of immunocompetent mammals, including humans. Infection results from environmental exposure to Coccidiodies, which is believed to grow as a soil saprophyte in arid deserts. To investigate hypotheses about the life history and evolution of Coccidioides, the genomes of several Onygenales, including C. immitis and C. posadasii; a close, nonpathogenic relative, Uncinocarpus reesii; and a more diverged pathogenic fungus, Histoplasma capsulatum, were sequenced and compared with those of 13 more distantly related Ascomycetes. This analysis identified increases and decreases in gene family size associated with a host/substrate shift from plants to animals in the Onygenales. In addition, comparison among Onygenales genomes revealed evolutionary changes in Coccidioides that may underlie its infectious phenotype, the identification of which may facilitate improved treatment and prevention of coccidioidomycosis. Overall, the results suggest that Coccidioides species are not soil saprophytes, but that they have evolved to remain associated with their dead animal hosts in soil, and that Coccidioides metabolism genes, membrane-related proteins, and putatively antigenic compounds have evolved in response to interaction with an animal host.

The GAAS Metagenomic Tool and Its Estimations of Viral and Microbial Average Genome Size in Four Major Biomes

Metagenomic studies characterize both the composition and diversity of uncultured viral and microbial communities. BLAST-based comparisons have typically been used for such analyses; however, sampling biases, high percentages of unknown sequences, and the use of arbitrary thresholds to find significant similarities can decrease the accuracy and validity of estimates. Here, we present Genome relative Abundance and Average Size (GAAS), a complete software package that provides improved estimates of community composition and average genome length for metagenomes in both textual and graphical formats. GAAS implements a novel methodology to control for sampling bias via length normalization, to adjust for multiple BLAST similarities by similarity weighting, and to select significant similarities using relative alignment lengths. In benchmark tests, the GAAS method was robust to both high percentages of unknown sequences and to variations in metagenomic sequence read lengths. Re-analysis of the Sargasso Sea virome using GAAS indicated that standard methodologies for metagenomic analysis may dramatically underestimate the abundance and importance of organisms with small genomes in environmental systems. Using GAAS, we conducted a meta-analysis of microbial and viral average genome lengths in over 150 metagenomes from four biomes to determine whether genome lengths vary consistently between and within biomes, and between microbial and viral communities from the same environment. Significant differences between biomes and within aquatic sub-biomes (oceans, hypersaline systems, freshwater, and microbialites) suggested that average genome length is a fundamental property of environments driven by factors at the sub-biome level. The behavior of paired viral and microbial metagenomes from the same environment indicated that microbial and viral average genome sizes are independent of each other, but indicative of community responses to stressors and environmental conditions.

Enabling a Community to Dissect an Organism: Overview of the Neurospora Functional Genomics Project

A consortium of investigators is engaged in a functional genomics project centered on the filamentous fungus Neurospora, with an eye to opening up the functional genomic analysis of all the filamentous fungi. The overall goal of the four interdependent projects in this effort is to accomplish functional genomics, annotation, and expression analyses of Neurospora crassa, a filamentous fungus that is an established model for the assemblage of over 250,000 species of non yeast fungi. Building from the completely sequenced 43-Mb Neurospora genome, Project 1 is pursuing the systematic disruption of genes through targeted gene replacements, phenotypic analysis of mutant strains, and their distribution to the scientific community at large. Project 2, through a primary focus in Annotation and Bioinformatics, has developed a platform for electronically capturing community feedback and data about the existing annotation, while building and maintaining a database to capture and display information about phenotypes. Oligonucleotide-based microarrays created in Project 3 are being used to collect baseline expression data for the nearly 11,000 distinguishable transcripts in Neurospora under various conditions of growth and development, and eventually to begin to analyze the global effects of loss of novel genes in strains created by Project 1. cDNA libraries generated in Project 4 document the overall complexity of expressed sequences in Neurospora, including alternative splicing alternative promoters and antisense transcripts. In addition, these studies have driven the assembly of an SNP map presently populated by nearly 300 markers that will greatly accelerate the positional cloning of genes.

Dynamics of Dengue Disease Severity Determined by the Interplay Between Viral Genetics and Serotype-Specific Immunity

The complex interactions between serotype-specific immunity and viral genetics drive dengue disease outcomes in endemic populations. Dissecting Dengue Dynamics In some parts of the world, mosquito bites cause more than a maddening itch. These tiny insect vessels can host dengue viruses (DENVs), which the mosquitoes transmit to human beings. DENVs cause disease states that range in severity from asymptomatic infection to the potentially fatal dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS). Infection with one DENV serotype offers a person lifelong protection against that serotype and temporary immune protection from severe dengue disease caused by a different DENV serotype. But once this immune-protected period passes, the individual becomes at increased risk for more severe disease. “What line separates the rumble of the wheels from the howl of the wolves?”1 Now, OhAinle et al. dissect the complex interactions between DENV immunity and viral genetics that drive dengue disease outcomes in endemic Nicaraguan populations. To address this worldwide public health problem, the authors carried out two independent clinical studies of dengue epidemics in Managua, Nicaragua, from 2004/5 to 2008/9. In these populations, patients displayed a striking increase in severe dengue disease outcomes caused by DENV-2 infection across several epidemic seasons. The increase in severe cases in later seasons could be attributed, in part, to waning DENV-1 immunity. However, OhAinle et al. also determined that the increased risk of severe DENV-2–related disease correlated with substitution, in patients, of the Asian/American DENV-2 NI-1 clade with a new virus clade, NI-2B. The authors then molecularly characterized the various viral genomes, assessed the viruses’ fitness levels in vitro, and measured the extent of viremia in patient blood samples; experimental results supported the emergence of a fitter virus in later epidemic seasons. OhAinle et al. also found that the NI-1 clade displayed enhanced virulence specifically in children who had been immune to DENV-1; in contrast, DENV-3 immunity correlated with more severe disease resulting from NI-2B infections. These data highlight the complex interplay between viral genetics and population dynamics that increase the risk of severe dengue disease. This study offers insights into viral evolution and the effects of a patient’s immunological background on viral fitness and virulence. A detailed understanding of DENV epidemiology and pathogenesis should aid in the design of a much-needed dengue vaccine. 1From Invisible Cities by Italo Calvino The rapid spread of dengue is a worldwide public health problem. In two clinical studies of dengue in Managua, Nicaragua, we observed an abrupt increase in disease severity across several epidemic seasons of dengue virus serotype 2 (DENV-2) transmission. Waning DENV-1 immunity appeared to increase the risk of severe disease in subsequent DENV-2 infections after a period of cross-protection. The increase in severity coincided with replacement of the Asian/American DENV-2 NI-1 clade with a new virus clade, NI-2B. In vitro analyses of viral isolates from the two clades and analysis of viremia in patient blood samples support the emergence of a fitter virus in later, relative to earlier, epidemic seasons. In addition, the NI-1 clade of viruses was more virulent specifically in children who were immune to DENV-1, whereas DENV-3 immunity was associated with more severe disease among NI-2B infections. Our data demonstrate that the complex interaction between viral genetics and population dynamics of serotype-specific immunity contributes to the risk of severe dengue disease. Furthermore, this work provides insights into viral evolution and the interaction between viral and immunological determinants of viral fitness and virulence.

A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454

We present an automated, high throughput library construction process for 454 technology. Sample handling errors and cross-contamination are minimized via end-to-end barcoding of plasticware, along with molecular DNA barcoding of constructs. Automation-friendly magnetic bead-based size selection and cleanup steps have been devised, eliminating major bottlenecks and significant sources of error. Using this methodology, one technician can create 96 sequence-ready 454 libraries in 2 days, a dramatic improvement over the standard method.

Trends in Patterns of Dengue Transmission over 4 Years in a Pediatric Cohort Study in Nicaragua

BACKGROUND Dengue is the most prevalent mosquito-borne viral disease in humans and a major urban public health problem worldwide. METHODS A prospective cohort study of approximately 3800 children initially aged 2-9 years was established in Managua, Nicaragua, in 2004 to study the natural history of dengue transmission in an urban pediatric population. Blood samples from healthy subjects were collected annually prior to the dengue season, and identification of dengue cases occurred via enhanced passive surveillance at the study health center. RESULTS Over the first four years of the study, seroprevalence of anti-dengue virus (DENV) antibodies increased from 22%-40% in the 2-year-old cohort and 90%-95% in the 9-year-old cohort. The incidence of symptomatic dengue cases and the ratio of inapparent to symptomatic DENV infection varied substantially from year to year. The switch in dominant transmission from DENV-1 to DENV-2 was accompanied by an increase in disease severity but, paradoxically, a decrease in transmission. Phylogeographic analysis of full-length DENV-2 sequences revealed strong geographic clustering of dengue cases. CONCLUSIONS This large-scale cohort study of dengue in the Americas demonstrates year-to-year variation of dengue within a pediatric population, revealing expected patterns in transmission while highlighting the impact of interventions, climate, and viral evolution.