Genomic epidemiology reveals multiple introductions of Zika virus into the United States.

Zika virus (ZIKV) is causing an unprecedented epidemic linked to severe congenital abnormalities. In July 2016, mosquito-borne ZIKV transmission was reported in the continental United States; since then, hundreds of locally acquired infections have been reported in Florida. To gain insights into the timing, source, and likely route(s) of ZIKV introduction, we tracked the virus from its first detection in Florida by sequencing ZIKV genomes from infected patients and Aedes aegypti mosquitoes. We show that at least 4 introductions, but potentially as many as 40, contributed to the outbreak in Florida and that local transmission is likely to have started in the spring of 2016-several months before its initial detection. By analysing surveillance and genetic data, we show that ZIKV moved among transmission zones in Miami. Our analyses show that most introductions were linked to the Caribbean, a finding corroborated by the high incidence rates and traffic volumes from the region into the Miami area. Our study provides an understanding of how ZIKV initiates transmission in new regions.

Genomic Epidemiology of Early SARS-CoV-2 Transmission Dynamics, Gujarat, India.

Limited genomic sampling in many high-incidence countries has impeded studies of severe respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic epidemiology. Consequently, critical questions remain about the generation and global distribution of virus genetic diversity. We investigated SARS-CoV-2 transmission dynamics in Gujarat, India, during the state's first epidemic wave to shed light on spread of the virus in one of the regions hardest hit by the pandemic. By integrating case data and 434 whole-genome sequences sampled across 20 districts, we reconstructed the epidemic dynamics and spatial spread of SARS-CoV-2 in Gujarat. Our findings indicate global and regional connectivity and population density were major drivers of the Gujarat outbreak. We detected >100 virus lineage introductions, most of which appear to be associated with international travel. Within Gujarat, virus dissemination occurred predominantly from densely populated regions to geographically proximate locations that had low population density, suggesting that urban centers contributed disproportionately to virus spread.

Genomic Surveillance of Yellow Fever Virus Epizootic in São Paulo, Brazil, 2016 - 2018.

São Paulo, a densely inhabited state in southeast Brazil that contains the fourth most populated city in the world, recently experienced its largest yellow fever virus (YFV) outbreak in decades. YFV does not normally circulate extensively in São Paulo, so most people were unvaccinated when the outbreak began. Surveillance in non-human primates (NHPs) is important for determining the magnitude and geographic extent of an epizootic, thereby helping to evaluate the risk of YFV spillover to humans. Data from infected NHPs can give more accurate insights into YFV spread than when using data from human cases alone. To contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo, we generated and analysed virus genomic data and epizootic case data from NHPs in São Paulo. We report the occurrence of three spatiotemporally distinct phases of the outbreak in São Paulo prior to February 2018. We generated 51 new virus genomes from YFV positive cases identified in 23 different municipalities in São Paulo, mostly sampled from NHPs between October 2016 and January 2018. Although we observe substantial heterogeneity in lineage dispersal velocities between phylogenetic branches, continuous phylogeographic analyses of generated YFV genomes suggest that YFV lineages spread in São Paulo at a mean rate of approximately 1km per day during all phases of the outbreak. Viral lineages from the first epizootic phase in northern São Paulo subsequently dispersed towards the south of the state to cause the second and third epizootic phases there. This alters our understanding of how YFV was introduced into the densely populated south of São Paulo state. Our results shed light on the sylvatic transmission of YFV in highly fragmented forested regions in São Paulo state and highlight the importance of continued surveillance of zoonotic pathogens in sentinel species.

Epidemiologic History and Genetic Diversity Origins of Chikungunya and Dengue Viruses, Paraguay.

Paraguay has been severely affected by emergent Zika and chikungunya viruses, and dengue virus is endemic. To learn more about the origins of genetic diversity and epidemiologic history of these viruses in Paraguay, we deployed portable sequencing technologies to strengthen genomic surveillance and determine the evolutionary and epidemic history of arthropod-borne viruses (arboviruses). Samples stored at the Paraguay National Central Laboratory were sequenced and subjected to phylogenetic analysis. Among 33 virus genomes generated, we identified 2 genotypes of chikungunya and 2 serotypes of dengue virus that circulated in Paraguay during 2014-2018; the main source of these virus lineages was estimated to be Brazil. The evolutionary history inferred by our analyses precisely matched the available travel history of the patients. The genomic surveillance approach used was valuable for describing the epidemiologic history of arboviruses and can be used to determine the origins and evolution of future arbovirus outbreaks.

Molecular epidemiology reveals the role of war in the spread of HIV in Ukraine.

Ukraine has one of the largest HIV epidemics in Europe, historically driven by people who inject drugs (PWID). The epidemic showed signs of stabilization in 2012, but the recent war in eastern Ukraine may be reigniting virus spread. We investigated the movement of HIV-infected people within Ukraine before and during the conflict. We analyzed HIV-1 subtype-A pol nucleotide sequences sampled during 2012-2015 from 427 patients of 24 regional AIDS centers and used phylogeographic analysis to reconstruct virus movement among different locations in Ukraine. We then tested for correlations between reported PWID behaviors and reconstructed patterns of virus spread. Our analyses suggest that Donetsk and Lugansk, two cities not controlled by the Ukrainian government in eastern Ukraine, were significant exporters of the virus to the rest of the country. Additional analyses showed that viral dissemination within the country changed after 2013. Spearman correlation analysis showed that incoming virus flow was correlated with the number of HIV-infected internally displaced people. Additionally, there was a correlation between more intensive virus movement and locations with a higher proportion of PWID practicing risky sexual behaviors. Our findings suggest that effective prevention responses should involve internally displaced people and people who frequently travel to war-affected regions. Scale-up of harm reduction services for PWID will be an important factor in preventing new local HIV outbreaks in Ukraine.

Yellow Fever Virus Reemergence and Spread in Southeast Brazil, 2016-2019.

The recent reemergence of yellow fever virus (YFV) in Brazil has raised serious concerns due to the rapid dissemination of the virus in the southeastern region. To better understand YFV genetic diversity and dynamics during the recent outbreak in southeastern Brazil, we generated 18 complete and nearly complete genomes from the peak of the epidemic curve from nonhuman primates (NHPs) and human infected cases across the Espírito Santo and Rio de Janeiro states. Genomic sequencing of 18 YFV genomes revealed the estimated timing, source, and likely routes of yellow fever virus transmission and dispersion during one of the largest outbreaks ever registered in Brazil. We showed that during the recent epidemic, YFV was reintroduced from Minas Gerais to the Espírito Santo and Rio de Janeiro states multiple times between 2016 and 2019. The analysis of data from portable sequencing could identify the corridor of spread of YFV. These findings reinforce the idea that continued genomic surveillance strategies can provide information on virus genetic diversity and transmission dynamics that might assist in understanding arbovirus epidemics.IMPORTANCE Arbovirus infections in Brazil, including yellow fever, dengue, zika, and chikungunya, result in considerable morbidity and mortality and are pressing public health concerns. However, our understanding of these outbreaks is hampered by the limited availability of genomic data. In this study, we investigated the genetic diversity and spatial distribution of YFV during the current outbreak by analyzing genomic data from areas in southeastern Brazil not covered by other previous studies. To gain insights into the routes of YFV introduction and dispersion, we tracked the virus by sequencing YFV genomes sampled from nonhuman primates and infected patients from the southeastern region. Our study provides an understanding of how YFV initiates transmission in new Brazilian regions and illustrates that genomics in the field can augment traditional approaches to infectious disease surveillance and control.

Transmission of hepatitis C virus in HIV-positive and PrEP-using MSM in England.

We sought to characterize risk factors and patterns of HCV transmission amongst men who have sex with men (MSM). MSM with recently acquired HCV (AHCV) were prospectively recruited ('clinic cohort') between January and September 2017. Clinical data and risk behaviours were identified and blood obtained for HCV whole genome sequencing. Phylogenetic analyses were performed, using sequences from this cohort and two other AHCV cohorts, to identify transmission clusters. Sixteen (40.0%) men in the clinic cohort were HIV-negative MSM. HIV-negative MSM were younger than HIV-positive MSM; most (81.3%) had taken HIV PrEP in the preceding year. Eighteen men (45.0%) reported injection drug use; most (34, 85.0%) reported noninjection drug use in the last year. Most in both groups reported condomless anal sex, fisting and sex in a group environment. Few (7, 17.5%) men thought partners may have had HCV. There were 52 sequences in the HCV genotype 1a phylogeny, 18 from the clinic cohort and 34 from other AHCV cohorts; 47 (90.4%) clustered with ≥1 other sequence. There were 7 clusters of 2-27 sequences; 6 clusters contained HIV-negative and HIV-positive MSM and 1 cluster only HIV-positive MSM. Four of these clusters were part of larger clusters first described in 2007. PrEP-using MSM are at risk of HCV, sharing similar risk factors to HIV-positive MSM. Phylogenetics highlights that PrEP-using and HIV-positive MSM are involved in the same HCV transmission networks. Few men demonstrated HCV awareness and risk reduction strategies should be expanded.

Genomic detection of a virus lineage replacement event of dengue virus serotype 2 in Brazil, 2019.

BACKGROUND Despite efforts to mitigate the impact of dengue virus (DENV) epidemics, the virus remains a public health problem in tropical and subtropical regions around the world. Most DENV cases in the Americas between January and July 2019 were reported in Brazil. São Paulo State in the southeast of Brazil has reported nearly half of all DENV infections in the country. OBJECTIVES To understand the origin and dynamics of the 2019 DENV outbreak. METHODS Here using portable nanopore sequencing we generated20 new DENV genome sequences from viremic patients with suspected dengue infection residing in two of the most-affected municipalities of São Paulo State, Araraquara and São José do Rio Preto. We conducted a comprehensive phylogenetic analysis with 1,630 global DENV strains to better understand the evolutionary history of the DENV lineages that currently circulate in the region. FINDINGS The new outbreak strains were classified as DENV2 genotype III (American/Asian genotype). Our analysis shows that the 2019 outbreak is the result of a novel DENV lineage that was recently introduced to Brazil from the Caribbean region. Dating phylogeographic analysis suggests that DENV2-III BR-4 was introduced to Brazil in or around early 2014, possibly from the Caribbean region. MAIN CONCLUSIONS Our study describes the early detection of a newly introduced and rapidly-expanding DENV2 virus lineage in Brazil.

Characterising West Nile virus epidemiology in Israel using a transmission suitability index.

BackgroundClimate is a major factor in the epidemiology of West Nile virus (WNV), a pathogen increasingly pervasive worldwide. Cases increased during 2018 in Israel, the United States and Europe.AimWe set to retrospectively understand the spatial and temporal determinants of WNV transmission in Israel, as a case study for the possible effects of climate on virus spread.MethodsWe employed a suitability index to WNV, parameterising it with prior knowledge pertaining to a bird reservoir and Culex species, using local time series of temperature and humidity as inputs. The predicted suitability index was compared with confirmed WNV cases in Israel (2016-2018).ResultsThe suitability index was highly associated with WNV cases in Israel, with correlation coefficients of 0.91 (p value = 4 × 10- 5), 0.68 (p = 0.016) and 0.9 (p = 2 × 10- 4) in 2016, 2017 and 2018, respectively. The fluctuations in the number of WNV cases between the years were explained by higher area under the index curve. A new WNV seasonal mode was identified in the south-east of Israel, along the Great Rift Valley, characterised by two yearly peaks (spring and autumn), distinct from the already known single summer peak in the rest of Israel.ConclusionsBy producing a detailed geotemporal estimate of transmission potential and its determinants in Israel, our study promotes a better understanding of WNV epidemiology and has the potential to inform future public health responses. The proposed approach further provides opportunities for retrospective and prospective mechanistic modelling of WNV epidemiology and its associated climatic drivers.

Early Genomic Detection of Cosmopolitan Genotype of Dengue Virus Serotype 2, Angola, 2018.

We used portable genome sequencing to investigate reported dengue virus transmission in Angola. Our results show that autochthonous transmission of dengue serotype 2 (cosmopolitan genotype) occurred in January 2018.

Metagenomic Next-Generation Sequencing of the 2014 Ebola Virus Disease Outbreak in the Democratic Republic of the Congo.

We applied metagenomic next-generation sequencing (mNGS) to detect Zaire Ebola virus (EBOV) and other potential pathogens from whole-blood samples from 70 patients with suspected Ebola hemorrhagic fever during a 2014 outbreak in Boende, Democratic Republic of the Congo (DRC) and correlated these findings with clinical symptoms. Twenty of 31 patients (64.5%) tested in Kinshasa, DRC, were EBOV positive by quantitative reverse transcriptase PCR (qRT-PCR). Despite partial degradation of sample RNA during shipping and handling, mNGS followed by EBOV-specific capture probe enrichment in a U.S. genomics laboratory identified EBOV reads in 22 of 70 samples (31.4%) versus in 21 of 70 (30.0%) EBOV-positive samples by repeat qRT-PCR (overall concordance = 87.1%). Reads from Plasmodium falciparum (malaria) were detected in 21 patients, of which at least 9 (42.9%) were coinfected with EBOV. Other positive viral detections included hepatitis B virus (n = 2), human pegivirus 1 (n = 2), Epstein-Barr virus (n = 9), and Orungo virus (n = 1), a virus in the Reoviridae family. The patient with Orungo virus infection presented with an acute febrile illness and died rapidly from massive hemorrhage and dehydration. Although the patient's blood sample was negative by EBOV qRT-PCR testing, identification of viral reads by mNGS confirmed the presence of EBOV coinfection. In total, 9 new EBOV genomes (3 complete genomes, and an additional 6 ≥50% complete) were assembled. Relaxed molecular clock phylogenetic analysis demonstrated a molecular evolutionary rate for the Boende strain 4 to 10× slower than that of other Ebola lineages. These results demonstrate the utility of mNGS in broad-based pathogen detection and outbreak surveillance.

Birth of a W sex chromosome by horizontal transfer of Wolbachia bacterial symbiont genome.

Sex determination is a fundamental developmental pathway governing male and female differentiation, with profound implications for morphology, reproductive strategies, and behavior. In animals, sex differences between males and females are generally determined by genetic factors carried by sex chromosomes. Sex chromosomes are remarkably variable in origin and can differ even between closely related species, indicating that transitions occur frequently and independently in different groups of organisms. The evolutionary causes underlying sex chromosome turnover are poorly understood, however. Here we provide evidence indicating that Wolbachia bacterial endosymbionts triggered the evolution of new sex chromosomes in the common pillbug Armadillidium vulgare We identified a 3-Mb insert of a feminizing Wolbachia genome that was recently transferred into the pillbug nuclear genome. The Wolbachia insert shows perfect linkage to the female sex, occurs in a male genetic background (i.e., lacking the ancestral W female sex chromosome), and is hemizygous. Our results support the conclusion that the Wolbachia insert is now acting as a female sex-determining region in pillbugs, and that the chromosome carrying the insert is a new W sex chromosome. Thus, bacteria-to-animal horizontal genome transfer represents a remarkable mechanism underpinning the birth of sex chromosomes. We conclude that sex ratio distorters, such as Wolbachia endosymbionts, can be powerful agents of evolutionary transitions in sex determination systems in animals.

Molecular Evolution, Diversity, and Adaptation of Influenza A(H7N9) Viruses in China.

The substantial increase in prevalence and emergence of antigenically divergent or highly pathogenic influenza A(H7N9) viruses during 2016-17 raises concerns about the epizootic potential of these viruses. We investigated the evolution and adaptation of H7N9 viruses by analyzing available data and newly generated virus sequences isolated in Guangdong Province, China, during 2015-2017. Phylogenetic analyses showed that circulating H7N9 viruses belong to distinct lineages with differing spatial distributions. Hemagglutination inhibition assays performed on serum samples from patients infected with these viruses identified 3 antigenic clusters for 16 strains of different virus lineages. We used ancestral sequence reconstruction to identify parallel amino acid changes on multiple separate lineages. We inferred that mutations in hemagglutinin occur primarily at sites involved in receptor recognition or antigenicity. Our results indicate that highly pathogenic strains likely emerged from viruses circulating in eastern Guangdong Province during March 2016 and are associated with a high rate of adaptive molecular evolution.

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.

Distinct Zika Virus Lineage in Salvador, Bahia, Brazil.

Sequencing of isolates from patients in Bahia, Brazil, where most Zika virus cases in Brazil have been reported, resulted in 11 whole and partial Zika virus genomes. Phylogenetic analyses revealed a well-supported Bahia-specific Zika virus lineage, which indicates sustained Zika virus circulation in Salvador, Bahia's capital city, since mid-2014.

Genetic diversity and phenotypic associations of feline caliciviruses from cats in Switzerland.

Feline calicivirus (FCV) is a common viral pathogen in domestic cats worldwide. The variable regions of the capsid (VP1) gene of FCV have one of the highest recorded rates of molecular evolution. Understanding the genetic diversity and phylogeny of FCV is a prerequisite to exploring the epidemiology and pathogenesis of this virus and to the development of efficacious vaccine strategies. In this study, we undertook a nationwide molecular characterization of FCV using for the first time nearly complete capsid (VP1) gene sequences. Sequences from 66 FCV samples were used to investigate the correlation between viral phylogeny and several traits, including geographic origin, signalment, husbandry, FCV vaccination and co-infections. Codon-based nucleotide alignment showed that individual nucleotides and their corresponding amino acid sites were either invariant or highly variable. Using a threshold of 20 % genetic distance in variable region E, FCV samples were grouped into 52 strains, 10 of which comprised two to three samples. Significant associations between FCV phylogeny and host characteristics were found, specifically the pedigree status of the cats, and two well-supported lineages were identified in which the current FCV strain definition was confounded. No correlation between viral genetic distances and geographic distances was evident. The greater resolution of the FCV phylogeny in this study compared to previous studies can be attributed to our use of more conserved regions of the capsid (VP1) gene; nonetheless, our results were still hampered by sequence saturation. The study highlights the need for whole-genome sequences for FCV phylogeny studies.

The genome of the nucleopolyhedrosis-causing virus from Tipula oleracea sheds new light on the Nudiviridae family.

UNLABELLED: A large double-stranded DNA (dsDNA) virus that produces occlusion bodies, typical of baculoviruses, has been described to infect crane fly larvae of the genus Tipula (Diptera, Tipulidae). Because of a lack of genomic data, this virus has remained unclassified. Electron microscopy of an archival virus isolated from Tipula oleracea, T. oleracea nudivirus (ToNV), showed irregularly shaped occlusion bodies measuring from 2 to 5 µm in length and 2 µm in middiameter, filled with rod-shape virions containing single nucleocapsids within a bilayer envelope. Whole-genome amplification and Roche 454 sequencing revealed a complete circular genome sequence of 145.7 kb, containing five direct repeat regions. We predicted 131 open reading frames, including a homolog of the polyhedrin gene encoding the major occlusion body protein of T. paludosa nucleopolyhedrovirus (NPV). BLAST searches demonstrated that ToNV had 21 of the 37 baculovirus core genes but shared 52 genes with nudiviruses (NVs). Phylogenomic analyses indicated that ToNV clearly belongs to the Nudiviridae family but should probably be assigned to a new genus. Among nudiviruses, ToNV was most closely related to the Penaeus monodon NV and Heliothis zea NV clade but distantly related to Drosophila innubia NV, the other nudivirus infecting a Diptera. Lastly, ToNV was found to be most closely related to the nuvidirus ancestor of bracoviruses. This was also reflected in terms of gene content, as ToNV was the only known exogenous virus harboring homologs of the Cc50C22.6 and 27b (Cc50C22.7) genes found in the nudiviral genomic cluster involved in bracovirus particle production. IMPORTANCE: The Nudiviridae is a family of arthropod dsDNA viruses from which striking cases of endogenization have been reported (i.e., symbiotic bracoviruses deriving from a nudivirus and the endogenous nudivirus of the brown planthopper). Although related to baculoviruses, relatively little is known about the genomic diversity of exogenous nudiviruses. Here, we characterized, morphologically and genetically, an archival sample of the Tipula oleracea nudivirus (ToNV), which has the particularity of forming occlusion bodies. Comparative genomic and phylogenomic analyses showed ToNV to be to date the closest known relative of the exogenous ancestor of bracoviruses and that ToNV should be assigned to a new genus. Moreover, we revised the homology relationships of nudiviral genes and identified a new set of 32 core genes for the Nudiviridae, of which 21 were also baculovirus core genes. These findings provide important insights into the evolutionary history of large arthropod dsDNA viruses.

Functional annotation of Cotesia congregata bracovirus: identification of viral genes expressed in parasitized host immune tissues.

UNLABELLED: Bracoviruses (BVs) from the Polydnaviridae family are symbiotic viruses used as biological weapons by parasitoid wasps to manipulate lepidopteran host physiology and induce parasitism success. BV particles are produced by wasp ovaries and injected along with the eggs into the caterpillar host body, where viral gene expression is necessary for wasp development. Recent sequencing of the proviral genome of Cotesia congregata BV (CcBV) identified 222 predicted virulence genes present on 35 proviral segments integrated into the wasp genome. To date, the expressions of only a few selected candidate virulence genes have been studied in the caterpillar host, and we lacked a global vision of viral gene expression. In this study, a large-scale transcriptomic analysis by 454 sequencing of two immune tissues (fat body and hemocytes) of parasitized Manduca sexta caterpillar hosts allowed the detection of expression of 88 CcBV genes expressed 24 h after the onset of parasitism. We linked the expression profiles of these genes to several factors, showing that different regulatory mechanisms control viral gene expression in the host. These factors include the presence of signal peptides in encoded proteins, diversification of promoter regions, and, more surprisingly, gene position on the proviral genome. Indeed, most genes for which expression could be detected are localized in particular proviral regions globally producing higher numbers of circles. Moreover, this polydnavirus (PDV) transcriptomic analysis also reveals that a majority of CcBV genes possess at least one intron and an arthropod transcription start site, consistent with an insect origin of these virulence genes. IMPORTANCE: Bracoviruses (BVs) are symbiotic polydnaviruses used by parasitoid wasps to manipulate lepidopteran host physiology, ensuring wasp offspring survival. To date, the expressions of only a few selected candidate BV virulence genes have been studied in caterpillar hosts. We performed a large-scale analysis of BV gene expression in two immune tissues of Manduca sexta caterpillars parasitized by Cotesia congregata wasps. Genes for which expression could be detected corresponded to genes localized in particular regions of the viral genome globally producing higher numbers of circles. Our study thus brings an original global vision of viral gene expression and paves the way to the determination of the regulatory mechanisms enabling the expression of BV genes in targeted organisms, such as major insect pests. In addition, we identify sequence features suggesting that most BV virulence genes were acquired from insect genomes.

Genomic diversity in European Spodoptera exigua multiple nucleopolyhedrovirus isolates.

Key virus traits such as virulence and transmission strategies rely on genetic variation that results in functional changes in the interactions between hosts and viruses. Here, comparative genomic analyses of seven isolates of Spodoptera exigua multiple nucleopolyhedrovirus (SeMNPV) with differing phenotypes were employed to pinpoint candidate genes that may be involved in host-virus interactions. These isolates obtained after vertical or horizontal transmission of infection in insects differed in virulence. Apart from one genome containing a piggyBac transposon, all European SeMNPV isolates had a similar genome size and content. Complete genome analyses of single nucleotide polymorphisms and insertions/deletions identified mutations in 48 ORFs that could result in functional changes. Among these, 13 ORFs could be correlated with particular phenotypic characteristics of SeMNPV isolates. Mutations were found in all gene functional classes and most of the changes we highlighted could potentially be associated with differences in transmission. The regulation of DNA replication (helicase, lef-7) and transcription (lef-9, p47) might be important for the establishment of sublethal infection prior to and following vertical transmission. Virus-host cell interactions also appear instrumental in the modulation of viral transmission as significant mutations were detected in virion proteins involved in primary (AC150) or secondary infections (ME35) and in apoptosis inhibition (IAP2, AC134). Baculovirus populations naturally harbour high genomic variation located in genes involved at different levels of the complex interactions between virus and host during the course of an infection. The comparative analyses performed here suggest that the differences in baculovirus virulence and transmission phenotypes involve multiple molecular pathways.