Spread of Chikungunya Virus East/Central/South African Genotype in Northeast Brazil.

We investigated an outbreak of exanthematous illness in Maceió by using molecular surveillance; 76% of samples tested positive for chikungunya virus. Genetic analysis of 23 newly generated genomes identified the East/Central/South African genotype, suggesting that this lineage has persisted since mid-2014 in Brazil and may spread in the Americas and beyond.

An ensemble strategy that significantly improves de novo assembly of microbial genomes from metagenomic next-generation sequencing data.

Next-generation sequencing (NGS) approaches rapidly produce millions to billions of short reads, which allow pathogen detection and discovery in human clinical, animal and environmental samples. A major limitation of sequence homology-based identification for highly divergent microorganisms is the short length of reads generated by most highly parallel sequencing technologies. Short reads require a high level of sequence similarities to annotated genes to confidently predict gene function or homology. Such recognition of highly divergent homologues can be improved by reference-free (de novo) assembly of short overlapping sequence reads into larger contigs. We describe an ensemble strategy that integrates the sequential use of various de Bruijn graph and overlap-layout-consensus assemblers with a novel partitioned sub-assembly approach. We also proposed new quality metrics that are suitable for evaluating metagenome de novo assembly. We demonstrate that this new ensemble strategy tested using in silico spike-in, clinical and environmental NGS datasets achieved significantly better contigs than current approaches.

Unifying the spatial epidemiology and molecular evolution of emerging epidemics.

We introduce a conceptual bridge between the previously unlinked fields of phylogenetics and mathematical spatial ecology, which enables the spatial parameters of an emerging epidemic to be directly estimated from sampled pathogen genome sequences. By using phylogenetic history to correct for spatial autocorrelation, we illustrate how a fundamental spatial variable, the diffusion coefficient, can be estimated using robust nonparametric statistics, and how heterogeneity in dispersal can be readily quantified. We apply this framework to the spread of the West Nile virus across North America, an important recent instance of spatial invasion by an emerging infectious disease. We demonstrate that the dispersal of West Nile virus is greater and far more variable than previously measured, such that its dissemination was critically determined by rare, long-range movements that are unlikely to be discerned during field observations. Our results indicate that, by ignoring this heterogeneity, previous models of the epidemic have substantially overestimated its basic reproductive number. More generally, our approach demonstrates that easily obtainable genetic data can be used to measure the spatial dynamics of natural populations that are otherwise difficult or costly to quantify.

The perils of pathogen discovery: origin of a novel parvovirus-like hybrid genome traced to nucleic acid extraction spin columns.

Next-generation sequencing was used for discovery and de novo assembly of a novel, highly divergent DNA virus at the interface between the Parvoviridae and Circoviridae. The virus, provisionally named parvovirus-like hybrid virus (PHV), is nearly identical by sequence to another DNA virus, NIH-CQV, previously detected in Chinese patients with seronegative (non-A-E) hepatitis. Although we initially detected PHV in a wide range of clinical samples, with all strains sharing ∼99% nucleotide and amino acid identity with each other and with NIH-CQV, the exact origin of the virus was eventually traced to contaminated silica-binding spin columns used for nucleic acid extraction. Definitive confirmation of the origin of PHV, and presumably NIH-CQV, was obtained by in-depth analyses of water eluted through contaminated spin columns. Analysis of environmental metagenome libraries detected PHV sequences in coastal marine waters of North America, suggesting that a potential association between PHV and diatoms (algae) that generate the silica matrix used in the spin columns may have resulted in inadvertent viral contamination during manufacture. The confirmation of PHV/NIH-CQV as laboratory reagent contaminants and not bona fide infectious agents of humans underscores the rigorous approach needed to establish the validity of new viral genomes discovered by next-generation sequencing.

Infectivity in chimpanzees (Pan troglodytes) of plasma collected before HCV RNA detectability by FDA-licensed assays: implications for transfusion safety and HCV infection outcomes.

Serial plasma aliquots (50 mL) obtained from 10 commercial donors who converted from hepatitis C virus (HCV) RNA negative to positive were transfused into 2 chimpanzees to assess infectivity during early HCV infection. Plasma, obtained 4 days before HCV RNA detectability by licensed assays, transmitted HCV infection to chimpanzee X355. The infectious PCR-negative plasma was subsequently shown to be positive in 2 of 23 replicates using a sensitive transcription-mediated amplification (TMA) assay, and estimated to contain 1.2 HCV RNA copies/mL (60 copies/50 mL transfused). Plasma units obtained up to 8 weeks earlier were not infectious in a second susceptible chimp, even when from donors with low-level, intermittent HCV RNA detection. Chimp x355 developed acute viremia with subsequent seroconversion, but cleared both virus and Ab in 17 weeks. When rechallenged 38 months later with 6000 RNA copies/mL from the same donor, X355 was transiently reinfected and again rapidly lost all HCV markers. We conclude that: (1) transfusions can transmit HCV infection before RNA detection, but the interval of test-negative infectivity is very brief; (2) early "blips" of HCV RNA appear noninfectious and can be ignored when calculating residual transfusion risk; and (3) markers of HCV infection can be lost rapidly after exposure to low-dose inocula.

The fecal viral flora of wild rodents.

The frequent interactions of rodents with humans make them a common source of zoonotic infections. To obtain an initial unbiased measure of the viral diversity in the enteric tract of wild rodents we sequenced partially purified, randomly amplified viral RNA and DNA in the feces of 105 wild rodents (mouse, vole, and rat) collected in California and Virginia. We identified in decreasing frequency sequences related to the mammalian viruses families Circoviridae, Picobirnaviridae, Picornaviridae, Astroviridae, Parvoviridae, Papillomaviridae, Adenoviridae, and Coronaviridae. Seventeen small circular DNA genomes containing one or two replicase genes distantly related to the Circoviridae representing several potentially new viral families were characterized. In the Picornaviridae family two new candidate genera as well as a close genetic relative of the human pathogen Aichi virus were characterized. Fragments of the first mouse sapelovirus and picobirnaviruses were identified and the first murine astrovirus genome was characterized. A mouse papillomavirus genome and fragments of a novel adenovirus and adenovirus-associated virus were also sequenced. The next largest fraction of the rodent fecal virome was related to insect viruses of the Densoviridae, Iridoviridae, Polydnaviridae, Dicistroviriade, Bromoviridae, and Virgaviridae families followed by plant virus-related sequences in the Nanoviridae, Geminiviridae, Phycodnaviridae, Secoviridae, Partitiviridae, Tymoviridae, Alphaflexiviridae, and Tombusviridae families reflecting the largely insect and plant rodent diet. Phylogenetic analyses of full and partial viral genomes therefore revealed many previously unreported viral species, genera, and families. The close genetic similarities noted between some rodent and human viruses might reflect past zoonoses. This study increases our understanding of the viral diversity in wild rodents and highlights the large number of still uncharacterized viruses in mammals.

A divergent clade of circular single-stranded DNA viruses from pig feces.

Using metagenomics and molecular cloning methods, we characterized five novel small, circular viral genomes from pig feces that are distantly related to chimpanzee and porcine stool-associated circular viruses, (ChiSCV and PoSCV1). Phylogenetic analysis placed these viruses into a highly divergent clade of this rapidly growing new viral family. This new clade of viruses, provisionally named porcine stool-associated circular virus 2 and 3 (PoSCV2 and PoSCV3), encodes a stem-loop structure (presumably the origin of DNA replication) in the small intergenic region and a replication initiator protein commonly found in other biological systems that replicate their genomes via the rolling-circle mechanism. Furthermore, these viruses also exhibit three additional overlapping open reading frames in the large intergenic region between the capsid and replication initiator protein genes.

Differences in HIV-specific T cell responses between HIV-exposed and -unexposed HIV-seronegative individuals.

HIV-1-specific T lymphocyte responses in individuals exposed to HIV-1 but who remain persistently seronegative (HESNs) have been reported in some but not all previous studies. This study was designed to resolve unequivocally the question of whether HESNs make HIV-1-specific T cell responses. We performed a blind investigation to measure HIV-1-specific T cell responses in both HIV-1-serodiscordant couples and HIV-1-unexposed seronegative controls (HUSNs). We found low-frequency HIV-1-specific T cells in both HESNs and HUSNs but show that the response rates were higher over time in the former (P = 0.01). Furthermore, the magnitudes of the HIV-1-specific T cell responses were significantly higher among responding HESNs than among HUSNs over time (P = 0.002). In both groups, responses were mediated by CD4 T cells. The responses were mapped to single peptides, which often corresponded to epitopes restricted by multiple HLA-DR types that have previously been detected in HIV-1-infected patients. HIV-1-specific T cell responses in HUSNs and some HESNs likely represent cross-reactivity to self or foreign non-HIV-1 antigens. The significantly greater T cell responses in HESNs, including in two who were homozygous for CCR5Δ32, demonstrates that HIV-1-specific T cell responses can be induced or augmented by exposure to HIV-1 without infection.

Viral nucleic acids in live-attenuated vaccines: detection of minority variants and an adventitious virus.

Metagenomics and a panmicrobial microarray were used to examine eight live-attenuated viral vaccines. Viral nucleic acids in trivalent oral poliovirus (OPV), rubella, measles, yellow fever, varicella-zoster, multivalent measles/mumps/rubella, and two rotavirus live vaccines were partially purified, randomly amplified, and pyrosequenced. Over half a million sequence reads were generated covering from 20 to 99% of the attenuated viral genomes at depths reaching up to 8,000 reads per nucleotides. Mutations and minority variants, relative to vaccine strains, not known to affect attenuation were detected in OPV, mumps virus, and varicella-zoster virus. The anticipated detection of endogenous retroviral sequences from the producer avian and primate cells was confirmed. Avian leukosis virus (ALV), previously shown to be noninfectious for humans, was present as RNA in viral particles, while simian retrovirus (SRV) was present as genetically defective DNA. Rotarix, an orally administered rotavirus vaccine, contained porcine circovirus-1 (PCV1), a highly prevalent nonpathogenic pig virus, which has not been shown to be infectious in humans. Hybridization of vaccine nucleic acids to a panmicrobial microarray confirmed the presence of endogenous retroviral and PCV1 nucleic acids. Deep sequencing and microarrays can therefore detect attenuated virus sequence changes, minority variants, and adventitious viruses and help maintain the current safety record of live-attenuated viral vaccines.

Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection.

The precise identification of the HIV-1 envelope glycoprotein (Env) responsible for productive clinical infection could be instrumental in elucidating the molecular basis of HIV-1 transmission and in designing effective vaccines. Here, we developed a mathematical model of random viral evolution and, together with phylogenetic tree construction, used it to analyze 3,449 complete env sequences derived by single genome amplification from 102 subjects with acute HIV-1 (clade B) infection. Viral env genes evolving from individual transmitted or founder viruses generally exhibited a Poisson distribution of mutations and star-like phylogeny, which coalesced to an inferred consensus sequence at or near the estimated time of virus transmission. Overall, 78 of 102 subjects had evidence of productive clinical infection by a single virus, and 24 others had evidence of productive clinical infection by a minimum of two to five viruses. Phenotypic analysis of transmitted or early founder Envs revealed a consistent pattern of CCR5 dependence, masking of coreceptor binding regions, and equivalent or modestly enhanced resistance to the fusion inhibitor T1249 and broadly neutralizing antibodies compared with Envs from chronically infected subjects. Low multiplicity infection and limited viral evolution preceding peak viremia suggest a finite window of potential vulnerability of HIV-1 to vaccine-elicited immune responses, although phenotypic properties of transmitted Envs pose a formidable defense.

Novel circular DNA viruses in stool samples of wild-living chimpanzees.

Viral particles in stool samples from wild-living chimpanzees were analysed using random PCR amplification and sequencing. Sequences encoding proteins distantly related to the replicase protein of single-stranded circular DNA viruses were identified. Inverse PCR was used to amplify and sequence multiple small circular DNA viral genomes. The viral genomes were related in size and genome organization to vertebrate circoviruses and plant geminiviruses but with a different location for the stem-loop structure involved in rolling circle DNA replication. The replicase genes of these viruses were most closely related to those of the much smaller (approximately 1 kb) plant nanovirus circular DNA chromosomes. Because the viruses have characteristics of both animal and plant viruses, we named them chimpanzee stool-associated circular viruses (ChiSCV). Further metagenomic studies of animal samples will greatly increase our knowledge of viral diversity and evolution.

Cardioviruses are genetically diverse and cause common enteric infections in South Asian children.

Cardioviruses cause enteric infections in mice and rats which when disseminated have been associated with myocarditis, type 1 diabetes, encephalitis, and multiple sclerosis-like symptoms. Cardioviruses have also been detected at lower frequencies in other mammals. The Cardiovirus genus within the Picornaviridae family is currently made up of two viral species, Theilovirus and Encephalomyocarditis virus. Until recently, only a single strain of cardioviruses (Vilyuisk virus within the Theilovirus species) associated with a geographically restricted and prevalent encephalitis-like condition had been reported to occur in humans. A second theilovirus-related cardiovirus (Saffold virus [SAFV]) was reported in 2007 and subsequently found in respiratory secretions from children with respiratory problems and in stools of both healthy and diarrheic children. Using viral metagenomics, we identified RNA fragments related to SAFV in the stools of Pakistani and Afghani children with nonpolio acute flaccid paralysis (AFP). We sequenced three near-full-length genomes, showing the presence of divergent strains of SAFV and preliminary evidence of a distant recombination event between the ancestors of the Theiler-like viruses of rats and those of human SAFV. Further VP1 sequencing showed the presence of five new SAFV genotypes, doubling the reported genetic diversity of human and animal theiloviruses combined. Both AFP patients and healthy children in Pakistan were found to be excreting SAFV at high frequencies of 9 and 12%, respectively. Further studies are needed to examine the roles of these highly common and diverse SAFV genotypes in nonpolio AFP and other human diseases.

Evidence for persistent low-level viremia in individuals who control human immunodeficiency virus in the absence of antiretroviral therapy.

A subset of antiretroviral-untreated, human immunodeficiency virus (HIV)-infected individuals are able to maintain undetectable plasma HIV RNA levels in the absence of antiretroviral therapy. These "elite" controllers are of high interest as they may provide novel insights regarding host mechanisms of virus control. The degree to which these individuals have residual plasma viremia has not been well defined. We performed a longitudinal study of 46 elite controllers, defined as HIV-seropositive, antiretroviral-untreated individuals with plasma HIV RNA levels of <50 to 75 copies/ml. The median duration of HIV diagnosis was 13 years, the median baseline CD4(+) T-cell count was 753 cells/mm(3), and the median duration of follow-up was 16 months. Plasma and cellular HIV RNA levels were measured using the transcription-mediated amplification (TMA) assay (estimated limit of detection of <3.5 copies RNA/ml). A total of 1,117 TMA assays were performed (median of five time points/subject and four replicates/time point). All but one subject had detectable plasma HIV RNA on at least one time point, and 15 (33%) subjects had detectable RNA at all time points. The majority of controllers also had detectable cell-associated RNA and proviral DNA. A mixed-effect linear model showed no strong evidence of change in plasma RNA levels over time. In conclusion, the vast majority (98%) of elite controllers had measurable plasma HIV RNA, often at levels higher than that observed in antiretroviral-treated patients. This confirms the failure to eradicate the virus, even in these unique individuals who are able to reduce plasma viremia to very low levels without antiretroviral therapy.

A novel picornavirus associated with gastroenteritis.

A novel picornavirus genome was sequenced, showing 42.6%, 35.2%, and 44.6% of deduced amino acid identities corresponding to the P1, P2, and P3 regions, respectively, of the Aichi virus. Divergent strains of this new virus, which we named salivirus, were detected in 18 stool samples from Nigeria, Tunisia, Nepal, and the United States. A statistical association was seen between virus shedding and unexplained cases of gastroenteritis in Nepal (P = 0.0056). Viruses with approximately 90% nucleotide similarity, named klassevirus, were also recently reported in three cases of unexplained diarrhea from the United States and Australia and in sewage from Spain, reflecting a global distribution and supporting a pathogenic role for this new group of picornaviruses.

Rapid identification of known and new RNA viruses from animal tissues.

Viral surveillance programs or diagnostic labs occasionally obtain infectious samples that fail to be typed by available cell culture, serological, or nucleic acid tests. Five such samples, originating from insect pools, skunk brain, human feces and sewer effluent, collected between 1955 and 1980, resulted in pathology when inoculated into suckling mice. In this study, sequence-independent amplification of partially purified viral nucleic acids and small scale shotgun sequencing was used on mouse brain and muscle tissues. A single viral agent was identified in each sample. For each virus, between 16% to 57% of the viral genome was acquired by sequencing only 42-108 plasmid inserts. Viruses derived from human feces or sewer effluent belonged to the Picornaviridae family and showed between 80% to 91% amino acid identities to known picornaviruses. The complete polyprotein sequence of one virus showed strong similarity to a simian picornavirus sequence in the provisional Sapelovirus genus. Insects and skunk derived viral sequences exhibited amino acid identities ranging from 25% to 98% to the segmented genomes of viruses within the Reoviridae family. Two isolates were highly divergent: one is potentially a new species within the orthoreovirus genus, and the other is a new species within the orbivirus genus. We demonstrate that a simple, inexpensive, and rapid metagenomics approach is effective for identifying known and highly divergent new viruses in homogenized tissues of acutely infected mice.

Modeling sequence evolution in acute HIV-1 infection.

We describe a mathematical model and Monte Carlo (MC) simulation of viral evolution during acute infection. We consider both synchronous and asynchronous processes of viral infection of new target cells. The model enables an assessment of the expected sequence diversity in new HIV-1 infections originating from a single transmitted viral strain, estimation of the most recent common ancestor (MRCA) of the transmitted viral lineage, and estimation of the time to coalesce back to the MRCA. We also calculate the probability of the MRCA being the transmitted virus or an evolved variant. Excluding insertions and deletions, we assume HIV-1 evolves by base substitution without selection pressure during the earliest phase of HIV-1 infection prior to the immune response. Unlike phylogenetic methods that follow a lineage backwards to coalescence, we compare the observed data to a model of the diversification of a viral population forward in time. To illustrate the application of these methods, we provide detailed comparisons of the model and simulations results to 306 envelope sequences obtained from eight newly infected subjects at a single time point. The data from 68 patients were in good agreement with model predictions, and hence compatible with a single-strain infection evolving under no selection pressure. The diversity of the samples from the other two patients was too great to be explained by the model, suggesting multiple HIV-1-strains were transmitted. The model can also be applied to longitudinal patient data to estimate within-host viral evolutionary parameters.

Multiple independent origins of a protease inhibitor resistance mutation in salvage therapy patients.

BACKGROUND: Combination anti-viral therapies have reduced treatment failure rates by requiring multiple specific mutations to be selected on the same viral genome to impart high-level drug resistance. To determine if the common protease inhibitor resistance mutation L90M is only selected once or repeatedly on different HIV genetic backbones during the course of failed anti-viral therapies we analyzed a linked region of the viral genome during the evolution of multi-drug resistance. RESULTS: Using L90M allele specific PCR we amplified and sequenced gag-pro regions linked to very early L90M containing HIV variants prior to their emergence and detection as dominant viruses in 15 failed salvage therapy patients. The early minority L90M linked sequences were then compared to those of the later L90M viruses that came to dominate the plasma quasispecies. Using Bayesian evolutionary analysis sampling trees the emergence of L90M containing viruses was seen to take place on multiple occasion in 5 patients, only once for 2 patients and an undetermined number of time for the remaining 8 patients. CONCLUSION: These results indicate that early L90M mutants can frequently be displaced by viruses carrying independently selected L90M mutations rather than by descendents of the earlier mutants.

Small Circular Rep-Encoding Single-Stranded DNA Genomes in Peruvian Diarrhea Virome.

Metagenomic analysis of diarrhea samples revealed the presence of numerous human enteric viruses and small circular Rep-encoding single-stranded DNA (CRESS-DNA) genomes. One such genome was related to smacoviruses, while eight others were related to genomes reported in the feces of different mammals. The tropism of these CRESS-DNA viruses remains unknown.

Genetic Characterization and Classification of Human and Animal Sapoviruses.

Sapoviruses (SaVs) are enteric caliciviruses that have been detected in multiple mammalian species, including humans, pigs, mink, dogs, sea lions, chimpanzees, and rats. They show a high level of diversity. A SaV genome commonly encodes seven nonstructural proteins (NSs), including the RNA polymerase protein NS7, and two structural proteins (VP1 and VP2). We classified human and animal SaVs into 15 genogroups (G) based on available VP1 sequences, including three newly characterized genomes from this study. We sequenced the full length genomes of one new genogroup V (GV), one GVII and one GVIII porcine SaV using long range RT-PCR including newly designed forward primers located in the conserved motifs of the putative NS3, and also 5' RACE methods. We also determined the 5'- and 3'-ends of sea lion GV SaV and canine GXIII SaV. Although the complete genomic sequences of GIX-GXII, and GXV SaVs are unavailable, common features of SaV genomes include: 1) "GTG" at the 5'-end of the genome, and a short (9~14 nt) 5'-untranslated region; and 2) the first five amino acids (M [A/V] S [K/R] P) of the putative NS1 and the five amino acids (FEMEG) surrounding the putative cleavage site between NS7 and VP1 were conserved among the chimpanzee, two of five genogroups of pig (GV and GVIII), sea lion, canine, and human SaVs. In contrast, these two amino acid motifs were clearly different in three genogroups of porcine (GIII, GVI and GVII), and bat SaVs. Our results suggest that several animal SaVs have genetic similarities to human SaVs. However, the ability of SaVs to be transmitted between humans and animals is uncertain.