Search

Background The impact of raltegravir-resistant HIV-1 minority variants (MVs) on raltegravir treatment failure is unknown. Illumina sequencing offers greater throughput than 454, but sequence analysis tools for viral sequencing are needed. We evaluated Illumina and 454 for the detection of HIV-1 raltegravir-resistant MVs. Methods A5262 was a single-arm study of raltegravir and darunavir/ritonavir in treatment-naive patients. Pre-treatment plasma was obtained from 5 participants with raltegravir resistance at the time of virologic failure. A control library was created by pooling integrase clones at predefined proportions. Multiplexed sequencing was performed with Illumina and 454 platforms at comparable costs. Illumina sequence analysis was performed with the novel snp-assess tool and 454 sequencing was analyzed with V-Phaser. Results Illumina sequencing resulted in significantly higher sequence coverage and a 0.095% limit of detection. Illumina accurately detected all MVs in the control library at ≥0.5% and 7/10 MVs expected at 0.1%. 454 sequencing failed to detect any MVs at 0.1% with 5 false positive calls. For MVs detected in the patient samples by both 454 and Illumina, the correlation in the detected variant frequencies was high (R2 = 0.92, P

Dengue virus is the most prevalent mosquito-borne virus worldwide. In this study, we used pyrosequencing to analyze the whole viral genome of two mouse-adapted strains, D2S10 and D2S20, that induce a dengue hemorrhagic fever/dengue shock syndrome (DHF/DSS)-like lethal disease in mice lacking the type I and/or type II interferon receptors. Previous experiments with D2S10 indicated that N124D and K128E mutations in the envelope protein were responsible for the severe disease induced in mice compared to its parental strain PL046. Here we demonstrate that D2S20 is more virulent than D2S10 and captured the presence of five key amino acid mutations – T70I, N83D, and K122I in envelope (E), and A62T in nonstructural protein 2A (NS2A) and G605V in nonstructural protein 5 (NS5) – that may account for this. These findings set the foundation for further dissection of the viral determinants responsible for dengue disease manifestations in mouse models.

Hepatitis C virus (HCV) is an enveloped, positive-strand RNA virus belonging to the Hepacivirus genus in the Flaviviridae family. Seven confirmed genotypes (1–7) are generally distinguished by phylogenetic methods and pair-wise distance calculations [1, 2]. On the basis of full genome nucleotide sequences, HCV genotypes diverged from each other by a pair-wise distance of >30%. Individual genotypes can be further divided into more closely related subtypes that diverged by a pair-wise distance of 15%–30%. All viral genotypes retain their repertoire of colinear structural and nonstructural genes, thereby facilitating preliminary genotype classification on the basis of partial genome sequences of short fragments (approximately 300–400 nucleotides) in the structural core/E1 region and the nonstructural NS5B region [3]. However, full genome sequences remain indispensible for the detection of genome recombination events. While only recombinant HCV genotype 2 k/1b viruses have been found to actively circulate in the population so far [4–7], recombinants of genotypes 2/5, 2b/1b, 2b/1a, and 2i/6p have been detected in single isolates from humans [8–11]. Furthermore, the frequency of HCV intergenotype and intragenotype recombination may be underestimated because of the lack of robust detection methods [12, 13]. Although genotyping based on core/E1 or NS5B sequences has resulted in the provisional classification of a large number of subtype variants within each genotype, at least 1 but preferably ≥2 full genome sequences are required to confirm subtype designation [2, 12]. Accurate genotype and subtype classification is clinically important because major genotypes differ considerably in their response rates to treatment with pegylated interferon and ribavirin and with directly acting antiviral drugs (DAAs) that are designed mainly against genotype 1 isolates. The limited efficacy of DAAs against other genotypes has been shown for genotype 3 [14], and treatment response rates for drug regimens containing boceprevir (an NS3 protease inhibitor) or BMS-790052 (an NS5A inhibitor) were found to vary even between subtypes 1a and 1b [15, 16]. Different baseline frequencies of drug resistance mutations may account for these differences [17], and such variation might become more significant if interferon-free regimens containing 1 or more DAAs could indeed displace the current standard of care [18]. Until clinically tolerable DAAs with high barriers to resistance against the complete spectrum of HCV genotypes become available [18, 19], accurate subtype determination and, possibly, drug resistance profiling of the individual's viral population may guide the optimal choice of drugs. Full genome deep sequencing, as performed in this study, allows the identification of drug resistance mutations across the genome that may exist as dominant or minor variants in treatment-naive patients, thereby informing the design of therapeutic regimens.

Marine viruses impose a heavy mortality on their host bacteria, whereas at the same time the degree of viral resistance in marine bacteria appears to be high. Antagonistic coevolution—the reciprocal evolutionary change of interacting species—might reconcile these observations, if it leads to rapid and dynamic levels of viral resistance. Here we demonstrate the potential for extensive antagonistic coevolution between the ecologically important marine cyanobacterium Synechococcus and a lytic virus. In a 6-mo-long replicated chemostat experiment, Synechococcus sp. WH7803 and the virus (RIM8) underwent multiple coevolutionary cycles, leading to the rapid diversification of both host and virus. Over the course of the experiment, we detected between 4 and 13 newly evolved viral phenotypes (differing in host range) and between 4 and 11 newly evolved Synechococcus phenotypes (differing in viral resistance) in each chemostat. Genomic analysis of isolates identified several candidate genes in both the host and virus that might influence their interactions. Notably, none of the viral candidates were tail fiber genes, thought to be the primary determinants of host range in tailed bacteriophages, highlighting the difficulty in generalizing results from bacteriophage infecting γ-Proteobacteria. Finally, we show that pairwise virus–host coevolution may have broader community consequences; coevolution in the chemostat altered the sensitivity of Synechoccocus to a diverse suite of viruses, as well as the virus’ ability to infect additional Synechococcus strains. Our results indicate that rapid coevolution may contribute to the generation and maintenance of Synechococcus and virus diversity and thereby influence viral-mediated mortality of these key marine bacteria.

Most RNA viruses exist in their hosts as a heterogeneous population of related variants. Due to error prone replication, mutants are constantly generated which may differ in individual fitness from the population as a whole. Here we characterize three WNV isolates that contain, along with full-length genomes, mutants with large internal deletions to structural and nonstructural protein-coding regions. The isolates were all obtained from lorikeets that died from WNV at the Rio Grande Zoo in Albuquerque, NM between 2005 and 2007. The deletions are approximately 2 kb, in frame, and result in the elimination of the complete envelope, and portions of the prM and NS-1 proteins. In Vero cell culture, these internally deleted WNV genomes function as defective interfering particles, reducing the production of full-length virus when introduced at high multiplicities of infection. In mosquitoes, the shortened WNV genomes reduced infection and dissemination rates, and virus titers overall, and were not detected in legs or salivary secretions at 14 or 21 days post-infection. In mice, inoculation with internally deleted genomes did not attenuate pathogenesis relative to full-length or infectious clone derived virus, and shortened genomes were not detected in mice at the time of death. These observations provide evidence that large deletions may occur within flavivirus populations more frequently than has generally been appreciated and suggest that they impact population phenotype minimally. Additionally, our findings suggest that highly similar mutants may frequently occur in particular vertebrate hosts.

The Coccolithoviridae are a group of viruses which infect the marine coccolithophorid microalga Emiliania huxleyi . The Emiliania huxleyi viruses (known as EhVs) described herein have 160- to 180-nm diameter icosahedral structures, have genomes of approximately 400 kbp, and consist of more than 450 predicted coding sequences (CDSs). Here, we describe the genomic features of four newly sequenced coccolithoviruses (EhV-88, EhV-201, EhV-207, and EhV-208) together with their draft genome sequences and their annotations, highlighting the homology and heterogeneity of these genomes to the EhV-86 model reference genome.

Emiliania huxleyi virus 202 (EhV-202) is a member of the Coccolithoviridae , a group of viruses that infect the marine coccolithophorid Emiliania huxleyi . EhV-202 has a 160- to 180-nm-diameter icosahedral structure and a genome of approximately 407 kbp, consisting of 485 coding sequences (CDSs). Here we describe the genomic features of EhV-202, together with a draft genome sequence and its annotation, highlighting the homology and heterogeneity of this genome in comparison with the EhV-86 reference genome.

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.

While human leukocyte antigen B57 (HLA-B57) is associated with the spontaneous clearance of hepatitis C virus (HCV), the mechanisms behind this control remain unclear. Immunodominant CD8 + T cell responses against the B57-restricted epitopes comprised of residues 2629 to 2637 of nonstructural protein 5B (NS5B 2629–2637 ) (KSKKTPMGF) and E2 541–549 (NTRPPLGNW) were recently shown to be crucial in the control of HCV infection. Here, we investigated whether the selection of deleterious cytotoxic T lymphocyte (CTL) escape mutations in the NS5B KSKKTPMGF epitope might impair viral replication and contribute to the B57-mediated control of HCV. Common CTL escape mutations in this epitope were identified from a cohort of 374 HCV genotype 1a-infected subjects, and their impact on HCV replication assessed using a transient HCV replicon system. We demonstrate that while escape mutations at residue 2633 (position 5) of the epitope had little or no impact on HCV replication in vitro , mutations at residue 2629 (position 1) substantially impaired replication. Notably, the deleterious mutations at position 2629 were tightly linked in vivo to upstream mutations at residue 2626, which functioned to restore the replicative defects imparted by the deleterious escape mutations. These data suggest that the selection of costly escape mutations within the immunodominant NS5B KSKKTPMGF epitope may contribute in part to the control of HCV replication in B57-positive individuals and that persistence of HCV in B57-positive individuals may involve the development of specific secondary compensatory mutations. These findings are reminiscent of the selection of deleterious CTL escape and compensatory mutations by HLA-B57 in HIV-1 infection and, thus, may suggest a common mechanism by which alleles like HLA-B57 mediate protection against these highly variable pathogens.

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., Author Summary The ability of HIV-1 and other highly variable pathogens to rapidly mutate to escape vaccine-induced immune responses represents a major hurdle to the development of effective vaccines to these highly persistent pathogens. Application of next-generation or deep sequencing technologies to the study of host pathogens could significantly improve our understanding of the mechanisms by which these pathogens subvert host immunity, and aid in the development of novel vaccines and therapeutics. Here, we developed a 454 deep sequencing approach to enable the sensitive detection of low-frequency viral variants across the entire HIV-1 genome. When applied to the acute phase of HIV-1 infection we observed that the majority of early, low frequency mutations represented viral adaptations to host cellular immune responses, evidence of strong host immunity developing during the early decline of peak viral load. Rapid viral escape from the most dominant immune responses however correlated with loss of this initial viral control, suggestive of the importance of mounting immune responses against more conserved regions of the virus. These data provide a greater understanding of the early evolutionary events subverting the ability of host immune responses to control early HIV-1 replication, yielding important insight into the design of more effective vaccine strategies.

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts., National Institute of Allergy and Infectious Diseases (U.S.), National Institutes of Health. Department of Health and Human Services (contract HHSN266200400001C), National Institutes of Health. Department of Health and Human Services(contract HHSN2722009000018C), Brazil. National Council for Scientific and Technological Development

University of Massachusetts Medical School. Worcester, MA, USA. University of Alabama at Birmingham. Birmingham, Al, USA. Biomed, Universidad de Carabobo. Maracay, Venezuela. Biomed, Universidad de Carabobo. Maracay, Venezuela. University of Massachusetts Medical School. Worcester, MA, USA. The Broad Institute of MIT and Harvard. Cambridge, MA, USA. University of Massachusetts Medical School. Worcester, MA, USA. Universidad Central de Venezuela. Caracas, Venezuela / Banco Municipal de Sangre del Distrito Capital. Caracas, Venezuela. University of Massachusetts Medical School. Worcester, MA, USA. Faculdade de Medicina de S?o Jos? do Rio Preto. Departamento de Doen?as Infecciosas e Parasit?rias. Laborat?rio de Pesquisas em Virologia. S?o Jos? do Rio Preto, SP, Brasil. Faculdade de Medicina de S?o Jos? do Rio Preto. Departamento de Doen?as Infecciosas e Parasit?rias. Laborat?rio de Pesquisas em Virologia. S?o Jos? do Rio Preto, SP, Brasil. Minist?rio da Sa?de. Secretaria de Vigil?ncia em Sa?de. Instituto Evandro Chagas. Instituto Nacional de Ci?ncia e Tecnologia para Febres Hemorr?gicas. Ananindeua, PA, Brasil / Universidade Federal de Minas Gerais. Instituto de Ci?ncias Biol?gicas. Departamento de Bioqu?mica e Imunologia. Instituto Nacional de Ci?ncia e Tecnologia em Dengue. Belo Horizonte, MG, Brasil. Minist?rio da Sa?de. Secretaria de Vigil?ncia em Sa?de. Instituto Evandro Chagas. Instituto Nacional de Ci?ncia e Tecnologia para Febres Hemorr?gicas. Ananindeua, PA, Brasil / Universidade Federal de Minas Gerais. Instituto de Ci?ncias Biol?gicas. Departamento de Bioqu?mica e Imunologia. Instituto Nacional de Ci?ncia e Tecnologia em Dengue. Belo Horizonte, MG, Brasil. Centers for Disease Control and Prevention. SanJuan, Puerto Rico, USA. Centers for Disease Control and Prevention. SanJuan, Puerto Rico, USA. Universidad Industrial de Santander. Santander, Colombia. Massachusetts Institute of Technology. Division of Heath Science and Technology. Cambridge, MA, USA. University of Alabama at Birmingham. Birmingham, Al, USA. The Broad Institute of MIT and Harvard. Cambridge, MA, USA. The Broad Institute of MIT and Harvard. Cambridge, MA, USA. University of Massachusetts Medical School. Worcester, MA, USA. Dengue virus currently causes 50?100 million infections annually. Comprehensive knowledge about the evolution of Dengue in response to selection pressure is currently unavailable, but would greatly enhance vaccine design efforts. In the current study, we sequenced 187 new dengue virus serotype 3 (DENV-3) genotype III whole genomes isolated from Asia and the Americas. We analyzed them together with previously-sequenced isolates to gain a more detailed understanding of the evolutionary adaptations existing in this prevalent American serotype. In order to analyze the phylogenetic dynamics of DENV-3 during outbreak periods; we incorporated datasets of 48 and 11 sequences spanning two major outbreaks in Venezuela during 2001 and 2007?2008, respectively. Our phylogenetic analysis of newly sequenced viruses shows that subsets of genomes cluster primarily by geographic location, and secondarily by time of virus isolation. DENV-3 genotype III sequences from Asia are significantly divergent from those from the Americas due to their geographical separation and subsequent speciation. We measured amino acid variation for the E protein by calculating the Shannon entropy at each position between Asian and American genomes. We found a cluster of seven amino acid substitutions having high variability within E protein domain III, which has previously been implicated in serotype-specific neutralization escape mutants. No novel mutations were found in the E protein of sequences isolated during either Venezuelan outbreak. Shannon entropy analysis of the NS5 polymerase mature protein revealed that a G374E mutation, in a region that contributes to interferon resistance in other flaviviruses by interfering with JAK-STAT signaling was present in both the Asian and American sequences from the 2007?2008 Venezuelan outbreak, but was absent in the sequences from the 2001 Venezuelan outbreak. In addition to E, several NS5 amino acid changes were unique to the 2007?2008 epidemic in Venezuela and may give additional insight into the adaptive response of DENV-3 at the population level.

Dengue is one of the most important infectious diseases of humans and has spread throughout much of the tropical and subtropical world. Despite this widespread dispersal, the determinants of dengue transmission in endemic populations are not well understood, although essential for virus control. To address this issue we performed a phylogeographic analysis of 751 complete genome sequences of dengue 1 virus (DENV-1) sampled from both rural (Dong Thap) and urban (Ho Chi Minh City) populations in southern Viet Nam during the period 2003–2008. We show that DENV-1 in Viet Nam exhibits strong spatial clustering, with likely importation from Cambodia on multiple occasions. Notably, multiple lineages of DENV-1 co-circulated in Ho Chi Minh City. That these lineages emerged at approximately the same time and dispersed over similar spatial regions suggests that they are of broadly equivalent fitness. We also observed an important relationship between the density of the human host population and the dispersion rate of dengue, such that DENV-1 tends to move from urban to rural populations, and that densely populated regions within Ho Chi Minh City act as major transmission foci. Despite these fluid dynamics, the dispersion rates of DENV-1 are relatively low, particularly in Ho Chi Minh City where the virus moves less than an average of 20 km/year. These low rates suggest a major role for mosquito-mediated dispersal, such that DENV-1 does not need to move great distances to infect a new host when there are abundant susceptibles, and imply that control measures should be directed toward the most densely populated urban environments., Author Summary Although dengue is a major cause of morbidity in many tropical and subtropical regions of the world, little is known about how the causative virus (dengue virus, DENV) spreads through endemic populations. To address this issue we undertook a phylogeny-based analysis of 751 complete genome sequences of DENV-1 sampled from patients in southern Vietnam during 2003–2008. We show that multiple viral lineages co-circulate within the urban area of Ho Chi Minh City (HCMC), and spread at approximately equivalent rates through overlapping geographical areas, suggesting that they are of equivalent fitness. We also observed that DENV-1 within HCMC tended to disperse from more to less densely populated regions, and that this city was the source population for DENV-1 in the rural area of Dong Thap. Despite the high prevalence of DENV-1 in southern Vietnam, viral dispersion rates were relatively low, especially in HCMC where they averaged less then 20 km/year. Such a low rate is consistent with predominantly mosquito-borne spatial dispersal of DENV-1 in this urban setting containing a large number of susceptibles. Together, these results suggest that dengue control measures such as insecticide spraying should be directed toward the most densely populated regions of localities where the virus is endemic.

In Cambodia, dengue virus (DENV) was first isolated in 1963 and has become endemic with peak epidemic during raining season. Since 2000, the Dengue National Control Program has reported from 10,000 to 40,000 cases per year with fatality rates ranging from 0.7 to 1.7. All four dengue serotypes are found circulating in Cambodia with alternative predominance of serotypes DENV-2 and DENV-3. The DENV-1 represents from 5% to 20% of all circulating viruses, depending upon the year. In this work, 79 clinical strains of DENV-1 were isolated between 2000 and 2009 and their genome fully sequenced. Four distinct lineages with different dynamics were identified. The main evolutionary drive was negative selective pressure but each lineage was characterized by the presence of specific mutations acquired through evolution. Coexistence, extinction and replacement of lineages occurred over the 10-year period. Lineages 1, 2 and 3 were all detected since 2000–2002 and disappeared in 2003, 2004–2005 and 2007, respectively. Lineages 1 and 2 displayed different dynamics. Lineage 1 was very diverse whereas lineage 2 was very homogeneous. Lineage 4 which derived from lineage 3 in 2003 remained the only one at the end of the sampling period in 2008–2009 owing to a selective sweep. The lineages dynamic of DENV-1 viruses and consequences for molecular epidemiology are discussed.

The Coccolithoviridae is a recently discovered group of viruses that infect the marine coccolithophorid Emiliania huxleyi. Emiliania huxleyi virus 84 (EhV-84) has a 160 -180 nm diameter icosahedral structure and a genome of approximately 400 kbp. Here we describe the structural and genomic features of this virus, together with a near complete draft genome sequence (~99%) and its annotation. This is the fourth genome sequence of a member of the coccolithovirus family.

Traditional study designs do not identify acute asymptomatic or pre-symptomatic dengue virus (DENV) infections, thus limiting our understanding of immunologic and viral factors that modulate infection outcome. In the 2006 and 2007 dengue seasons, we conducted a pilot index cluster study in Managua, Nicaragua, in which 442 persons living within 50 meters of 22 index cases identified through an ongoing pediatric cohort study were evaluated for DENV infection. Post-enrollment and pre-enrollment DENV infections were confirmed in 12 (2.7%) and 19 (4.3%) contacts, respectively. Five (42%) post-enrollment infections were asymptomatic, and DENV-2 was identified in 9 (75%) infections. Phylogenetic analysis with full-length DENV genomic sequence from contacts, index cases, and cohort dengue cases indicated focal transmission and infection outside the local area. We demonstrate the feasibility of identification of acute asymptomatic and pre-symptomatic cases in urban Latin America, the first report of such a study in the Americas, and identify age and concomitant immunity to DENV of contacts as a key factor in index cluster study design.

Background Comprehensive sequence characterization across the MHC is important for successful organ transplantation and genetic association studies. To this end, we have developed an automated sample preparation, molecular barcoding and multiplexing protocol for the amplification and sequence-determination of class I HLA loci. We have coupled this process to a novel HLA calling algorithm to determine the most likely pair of alleles at each locus. Results We have benchmarked our protocol with 270 HapMap individuals from four worldwide populations with 96.4% accuracy at 4-digit resolution. A variation of this initial protocol, more suitable for large sample sizes, in which molecular barcodes are added during PCR rather than library construction, was tested on 95 HapMap individuals with 98.6% accuracy at 4-digit resolution. Conclusions Next-generation sequencing on the 454 FLX Titanium platform is a reliable, efficient, and scalable technology for HLA typing.

We have sequenced the genomes of 18 isolates of the closely related human pathogenic fungi Coccidioides immitis and Coccidioides posadasii to more clearly elucidate population genomic structure, bringing the total number of sequenced genomes for each species to 10. Our data confirm earlier microsatellite-based findings that these species are genetically differentiated, but our population genomics approach reveals that hybridization and genetic introgression have recently occurred between the two species. The directionality of introgression is primarily from C. posadasii to C. immitis, and we find more than 800 genes exhibiting strong evidence of introgression in one or more sequenced isolates. We performed PCRbased sequencing of one region exhibiting introgression in 40 C. immitis isolates to confirm and better define the extent of gene flow between the species. We find more coding sequence than expected by chance in the introgressed regions, suggesting that natural selection may play a role in the observed genetic exchange. We find notable heterogeneity in repetitive sequence composition among the sequenced genomes and present the first detailed genome-wide profile of a repeat-induced point mutation (RIP) process distinctly different from what has been observed in Neurospora. We identify promiscuous HLA-I and HLA-II epitopes in both proteomes and discuss the possible implications of introgression and population genomic data for public health and vaccine candidate prioritization. This study highlights the importance of population genomic data for detecting subtle but potentially important phenomena such as introgression. [Supplemental material is available online at http://www.genome.org. Genome assemblies from this study have been submitted to NCBI (http://www.ncbi.nlm.nih.gov/) under Genome Project IDs 10735 and 12821. SNPs have been submitted to dbSNP (http://www.ncbi.nlm.nih.gov/projects/SNP/) under submitter handle BROAD_GENOMEBIO and submission batch CI_0001 (ss numbers are available in the Supplemental material). Sequencing data have been submitted to the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra) under accession nos. SRX012717–SRX012723.]

Background Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage. Methodology/Principal Findings To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling. Conclusions/Significance These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.

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., Author Summary Metagenomics uses DNA or RNA sequences isolated directly from the environment to determine what viruses or microorganisms exist in natural communities and what metabolic activities they encode. Typically, metagenomic sequences are compared to annotated sequences in public databases using the BLAST search tool. Our methods, implemented in the Genome relative Abundance and Average Size (GAAS) software, improve the way BLAST searches are processed to estimate the taxonomic composition of communities and their average genome length. GAAS provides a more accurate picture of community composition by correcting for a systematic sampling bias towards larger genomes, and is useful in situations where organisms with small genomes are abundant, such as disease outbreaks caused by small RNA viruses. Microbial average genome length relates to environmental complexity and the distribution of genome lengths describes community diversity. A study of the average genome length of viruses and microorganisms in four different biomes using GAAS on 169 metagenomes showed significantly different average genome sizes between biomes, and large variability within biomes as well. This also revealed that microbial and viral average genome sizes in the same environment are independent of each other, which reflects the different ways that microorganisms and viruses respond to stress and environmental conditions.

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.

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.

We report the discovery and validation of a set of single nucleotide polymorphisms (SNPs) between the reference Neurospora crassa strain Oak Ridge and the Mauriceville strain (FGSC 2555), of sufficient density to allow fine mapping of most loci. Sequencing of Mauriceville cDNAs and alignment to the completed genomic sequence of the Oak Ridge strain identified 19,087 putative SNPs. Of these, a subset was validated by cleaved amplified polymorphic sequence (CAPS), a simple and robust PCR-based assay that reliably distinguishes between SNP alleles. Experimental confirmation resulted in the development of 250 CAPS markers distributed evenly over the genome. To demonstrate the applicability of this map, we used bulked segregant analysis followed by interval mapping to locate the csp-1 mutation to a narrow region on LGI. Subsequently, we refined mapping resolution to 74 kbp by developing additional markers, resequenced the candidate gene, NCU02713.3, in the mutant background, and phenocopied the mutation by gene replacement in the WT strain. Together, these techniques demonstrate a generally applicable and straightforward approach for the isolation of novel genes from existing mutants. Data on both putative and validated SNPs are deposited in a customized public database at the Broad Institute, which encourages augmentation by community users.

Infection with genotype 4 of the Hepatitis C virus is common in Africa and the Mediterranean area, but has also been found at increasing frequencies in injection drug users in Europe and North America. Full length viral sequences to characterize viral diversity and structure have recently become available mostly for subtype 4a, and studies in Egypt and Saudi Arabia, where high proportions of subtype 4a infected patients exist, have begun to establish optimized treatment regimens. However knowledge about other subtype variants of genotype 4 present in less developed African states is lacking. In this study the full coding region from so far poorly characterized variants of HCV genotype 4 was amplified and sequenced using a long range PCR technique. Sequences were analyzed with respect to phylogenetic relationship, possible recombination and prominent sequence characteristics compared to other known HCV strains. We present for the first time two full-length sequences from the HCV genotype 4k, in addition to five strains from HCV genotypes 4d and 4f. Reference sequences for accurate HCV genotyping are required for optimized treatment, and a better knowledge of the global viral sequence diversity is needed to guide vaccines or new drugs effective in the world wide epidemic.

Acidobacteria constitute a still poorly studied phylum that is well represented in soils. Recent studies suggest that members of this phylum may be also abundant in deep-sea plankton, but their relative abundance and ecological role in this ecosystem are completely unknown. A recent screening of three metagenomic deep-sea libraries of bathypelagic plankton from the South Atlantic (1000 m depth), the Adriatic (1000 m depth) and the Ionian (3000 m depth) seas in the Mediterranean revealed an unexpected relative proportion of acidobacterial fosmids, which affiliated to the Solibacterales (Group 3), to the Group 11 and, most frequently, to the Group 6 of this diverse phylum. Here, we present the comparative analysis of 11 acidobacterial genome fragments containing the rrn operon from these Mediterranean libraries. A highly conserved syntenic region spanning up to 30 kb and containing up to 25 open reading frames was shared by Group 6 Acidobacteria. Synteny was also partially conserved in distantly related acidobacterial genome fragments derived from a metagenomic soil library, indicating a remarkable conservation of this genomic region within these Acidobacteria. A search for Acidobacteria-specific hits in directly comparable, available fosmid-end sequences from soil and marine metagenomic libraries showed a significant increase of their relative proportion in plankton libraries as a function of increasing depth reaching, at high depth, levels nearly comparable to those of soil. Thus, our results suggest that Acidobacteria are abundant and represent a significant proportion of the microbial community in the deep-sea ecosystem.

Dermatophytes are the most common cause of fungal infections worldwide, impacting millions of individuals annually. In the United States alone, this translates into an economic impact on the health care system estimated to exceed $400 million a year for treatment alone (8, 46). A recent survey conducted in 16 European countries has shown that more than one-third (35 to 40%) of the 90,000 participants were suffering with a fungal foot disease, mainly caused by dermatophytes (4). In a recent study, 22 to 50% of children in a U.S. day care center exhibited symptoms of dermatophytic scalp infections (2). Despite the prominence of dermatophyte infections and their resulting socioeconomic consequences, the research and medical communities lack a sophisticated understanding of these organisms' biology. Consequently, effective preventatives and therapeutics are lacking. These deficiencies are in large part due to the lack of genetic tools available for the study of these fungi and their host specificities. Several research efforts are now poised to produce genomic and molecular resources that will enable the molecular characterization of dermatophytes. These resources will soon be available for use by the larger research community to address many questions about the biology and pathogenesis of dermatophytes. This review will elaborate on the current status of these resources and their importance to the study of the dermatophytes.

Airborne microorganisms have been studied for centuries, but the majority of this research has relied on cultivation-dependent surveys that may not capture all of the microbial diversity in the atmosphere. As a result, our understanding of airborne microbial ecology is limited despite the relevance of airborne microbes to human health, various ecosystem functions, and environmental quality. Cultivation-independent surveys of small-subunit rRNA genes were conducted in order to identify the types of airborne bacteria and fungi found at a single site (Boulder, CO) and the temporal variability in the microbial assemblages over an 8-day period. We found that the air samples were dominated by ascomycete fungi of the Hypocreales order and a diverse array of bacteria, including members of the proteobacterial and Cytophaga-Flavobacterium-Bacteroides groups that are commonly found in comparable culture-independent surveys of airborne bacteria. Bacterium/fungus ratios varied by 2 orders of magnitude over the sampling period, and we observed large shifts in the phylogenetic diversity of bacteria present in the air samples collected on different dates, shifts that were not likely to be related to local meteorological conditions. We observed more phylogenetic similarity between bacteria collected from geographically distant sites than between bacteria collected from the same site on different days. These results suggest that outdoor air may harbor similar types of bacteria regardless of location and that the short-term temporal variability in airborne bacterial assemblages can be very large.

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 acccomplish 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 nonyeast 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.

The last decade has witnessed a revolution in the genomics of the fungal kingdom. Since the sequencing of the first fungus in 1996, the number of available fungal genome sequences has increased by an order of magnitude. Over 40 complete fungal genomes have been publicly released with an equal number currently being sequenced—representing the widest sampling of genomes from any eukaryotic kingdom. Moreover, many of these sequenced species form clusters of related organisms designed to enable comparative studies. These data provide an unparalleled opportunity to study the biology and evolution of this medically, industrially, and environmentally important kingdom. In addition, fungi also serve as model organisms for all eukaryotes. The available fungal genomic resource, coupled with the experimental tractability of the fungi, is accelerating research into the fundamental aspects of eukaryotic biology. We provide here an overview of available fungal genomes and highlight some of the biological insights that have been derived through their analysis. We also discuss insights into the fundamental cellular biology shared between fungi and other eukaryotic organisms.

The Coccolithoviridae are a recently discovered group of viruses that infect the marine coccolithophorid Emiliania huxleyi . Emiliania huxleyi virus 203 (EhV-203) has a 160- to 180-nm-diameter icosahedral structure and a genome of approximately 400 kbp, consisting of 464 coding sequences (CDSs). Here we describe the genomic features of EhV-203 together with a draft genome sequence and its annotation, highlighting the homology and heterogeneity of this genome in comparison with the EhV-86 reference genome.

Summary: PriSM is a set of algorithms designed to select and match degenerate primer pairs for the amplification of viral genomes. The design of panels of hundreds of primer pairs takes just hours using this program, compared with days using a manual approach. PriSM allows for rapid in silico optimization of primers for downstream applications such as sequencing. As a validation, PriSM was used to create an amplification primer panel for human immunodeficiency virus (HIV) Clade B. Availability: The program is freely available for use at: www.broadinstitute.org/perl/seq/specialprojects/primerDesign.cgi Contact: nlennon@broadinstitute.org Supplementary information: Supplementary data are available at Bioinformatics online.

To quantify and characterize N and C isotopic fractionation effects due to fungal transformation of organic substrates in forest ecosystems, we performed a field study in California and a meta-analysis of three additional studies conducted by others across the Northern Hemisphere. Basidiomycete fungal biomass was consistently enriched for the heavier isotope for C relative to substrate and either enriched or depleted for N relative to atmospheric N. Extent and pattern of fractionation was very variable, but the distinction between ectomycorrhizal and saprotrophic basidiomycetes was strongly supported, particularly when dual isotope analyses were performed. This differentiation, which we call the "EM-SAP Divide" holds for studies within a restricted ecosystem, but becomes less distinct over larger geographical regions, removing the rationale for using direct isotopic values from single specimens as diagnostic of ecophysiological role. For C, the EM-SAP Divide seems to reflect substrate effects, potentially due to differential access to recently synthesized versus recycled organic compounds, rather than distinct physiological pathways. Once substrate and ecophysiological role effects are removed, our meta-analysis suggests the existence of more than one mechanism causing C fractionations in fungi which is found equally in ectomycorrhizal and saprotrophic fungi. Similarly, a multimodal distribution of δ

Stable isotope analysis is a major tool used in ecosystem studies to establish pathways and rates of C exchange between various ecosystem components. Little is known about isotopic effects of many such components, especially microbes. Here we report on the discovery of an unexpected pattern of C isotopic discrimination by basidiomycete fungi with far-reaching consequences for our understanding of isotopic processing in ecosystems where these microbes mediate material transfers across trophic levels. We measured fractionation effects on three ecologically relevant basidiomycete species under controlled laboratory conditions. Sucrose derived from C3and C4plants is fractionated differentially by these microbes in a taxon-specific manner. The differentiation between mycorrhizal and saprotrophic fungi observed in the field by others is not explained by intrinsic discrimination patterns. Fractionation occurs during sugar uptake and is sensitive to the nonrandom distribution of stable isotopes in the sucrose molecule. The balance between respiratory physiology and fermentative physiology modulates the degree of fractionation. These discoveries disprove the assumption that fungal C processing does not significantly alter the distribution of stable C isotopes and provide the basis for a reevaluation of ecosystem models based on isotopic evidence that involve C transfer across microbial interfaces. We provide a mechanism to account for the observed differential discrimination effects.

Massively parallel sequencing offers the possibility of revolutionizing the study of viral populations by providing ultra deep sequencing (tens to hundreds of thousand fold coverage) of complete viral genomes. However, differentiation of true low frequency variants from sequencing errors remains challenging. We developed a software package, V-Phaser 2, for inferring intrahost diversity within viral populations. This program adds three major new methodologies to the state of the art: a technique to efficiently utilize paired end read data for calling phased variants, a new strategy to represent and infer length polymorphisms, and an in line filter for erroneous calls arising from systematic sequencing artifacts. We have also heavily optimized memory and run time performance. This combination of algorithmic and technical advances allows V-Phaser 2 to fully utilize extremely deep paired end sequencing data (such as generated by Illumina sequencers) to accurately infer low frequency intrahost variants in viral populations in reasonable time on a standard desktop computer. V-Phaser 2 was validated and compared to both QuRe and the original V-Phaser on three datasets obtained from two viral populations: a mixture of eight known strains of West Nile Virus (WNV) sequenced on both 454 Titanium and Illumina MiSeq and a mixture of twenty-four known strains of WNV sequenced only on 454 Titanium. V-Phaser 2 outperformed the other two programs in both sensitivity and specificity while using more than five fold less time and memory. We developed V-Phaser 2, a publicly available software tool (V-Phaser 2 can be accessed via: http://www.broadinstitute.org/scientific-community/science/projects/viral-genomics/v-phaser-2 and is freely available for academic use) that enables the efficient analysis of ultra-deep sequencing data produced by common next generation sequencing platforms for viral populations.

Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing methods and obtain the high sequence coverage needed to detect rare genetic variants (97% sensitivity and >97% specificity on control read sets. On data derived from a patient after four years of HIV-1 infection, V-Phaser detected 2,015 variants across the ∼10 kb genome, including 603 rare variants (, Author Summary New sequencing technologies provide unprecedented resolution to study pathogen populations, such as the single stranded RNA viruses HIV, dengue (DENV), and West Nile (WNV), and how they evolve within infected individuals in response to immune, therapeutic, and vaccine pressures. While these new technologies provide high volumes of data, these data contain process errors. To detect biological variants, especially those occurring at low frequencies in the population, these technologies require a method to differentiate biological variants from process errors with high sensitivity and specificity. To address this challenge, we introduce the V-Phaser algorithm, which distinguished the covariation of biological variants from that of process errors. We validate the method by measuring how frequently it correctly identifies variants and errors on actual read sets with known variation. Further, using data derived from a patient following four years of HIV-1 infection, we show that V-Phaser can detect biological variants at frequencies comparable to approaches that require prior knowledge. V-Phaser is available for download at: http://www.broadinstitute.org/scientific-community/software.

Viral diseases such as HIV/AIDS and Dengue have an enormous impact on human health worldwide. Despite this, application of new sequencing technologies to viral genomics has lagged. We are using genome sequence data to study how populations of single stranded RNA viruses, including HIV, Dengue, West Nile and Hepatitis C, evolve within infected individuals in response to host immune, therapeutic and vaccine pressures. To support this, we have developed high-throughput sequencing, assembly and population profiling pipelines based on 454 and Illumina technology that are tuned to the specific needs of viral sequencing. These strategies can capture full genome sequences and can profile sequence diversity at each residue in the genome with unprecedented sensitivity. Our analytical pipeline for 454 and Illumina data (i) generates complete genome assemblies from short sequencing reads derived from populations with high rates of variation; and (ii) detects and quantifies variants in these populations with high sensitivity, while differentiating true variants from process errors. Initial results with our viral sequencing and analysis pipeline are extremely promising. We detected rare variants to below 1% frequency, revolutionizing our ability to accurately assess the earliest events in viral evolution. We have demonstrated effective assessment of genome-wide diversity during acute HIV infection, enabling rapid, affordable, and highly sensitive identification of the earliest cellular immune responses to HIV. This has allowed us to detect earlier evolutionary events, demonstrating, for example, that HIV cytotoxic T-lymphocyte (CTL) escape can occur much faster than previously known. In addition, we have shown that the extent of intra-host diversity in Flaviviruses such as DENV and WNV is different between these closely related viruses with the latter exhibiting greater genetic diversity. The results presented here demonstrate the power of scalable, next generation sequencing-based methodologies as a genome-wide and unbiased global approach to profiling genomic diversity in intra-host populations of single stranded RNA viruses.

A better description of the extent and structure of genetic diversity in dengue virus (DENV) in endemic settings is central to its eventual control. To this end we determined the complete coding region sequence of 187 DENV-2 genomes and 68 E genes from viruses sampled from Vietnamese patients between 1995 and 2009. Strikingly, an episode of genotype replacement was observed, with Asian 1 lineage viruses entirely displacing the previously dominant Asian/American lineage viruses. This genotype replacement event also seems to have occurred within DENV-2 in Thailand and Cambodia, suggestive of a major difference in viral fitness. To determine the cause of this major evolutionary event we compared both the infectivity of the Asian 1 and Asian/American genotypes in mosquitoes and their viraemia levels in humans. Although there was little difference in infectivity in mosquitoes, we observed significantly higher plasma viraemia levels in paediatric patients infected with Asian 1 lineage viruses relative to Asian/American viruses, a phenotype that is predicted to result in a higher probability of human-to-mosquito transmission. These results provide a mechanistic basis to a marked change in the genetic structure of DENV-2 and more broadly underscore that an understanding of DENV evolutionary dynamics can inform the development of vaccines and anti-viral drugs., Author Summary Dengue virus (DENV) is a mosquito transmitted RNA virus. One of the most characteristic patterns of DENV evolution is that viral lineages, including whole genotypes, are born and die-out on a regular basis. However, the precise cause of such lineage turnover is unclear. To address this issue we explored the genome-scale genetic diversity and evolution of dengue serotype 2 virus (DENV-2) in Viet Nam, a country with a high burden of dengue disease. We observed that DENV viruses assigned to the Asian 1 lineage were likely introduced into southern Viet Nam in the late 1990's (perhaps from Thailand) and subsequently displaced the resident Asian/American lineage serotype 2 viruses. A similar scenario seems to have occurred in Thailand and Cambodia, such that there appears to have been a region-wide proliferation of Asian 1 lineage viruses. Investigation of Vietnamese patients experiencing DENV-2 infection revealed that Asian 1 viruses also attain higher virus levels in the blood than viruses of the Asian/American lineage. This difference in virus titre is likely to have a profound impact on viral fitness by increasing the probability of mosquito transmission, and therein providing Asian 1 lineage viruses with a selective advantage.

An automated method for constructing libraries for 454 sequencing significantly reduces the cost and time required., 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.