1. |
- Andersson, E, et al.
(author)
-
Near full-length HIV-1 genome sequencing in newly diagnosed individuals in Sweden
- 2019
-
In: VIRUS EVOLUTION. - : Oxford University Press (OUP). - 2057-1577. ; 5, s. S2-S2
-
Conference paper (other academic/artistic)abstract
- The Swedish HIV-1 epidemic is characterized by a high diversity in HIV subtypes and recombinants as a result of migration. To study the time from infection through viral diversification, transmission patterns, and drug resistance in minor quasispecies, a robust protocol for pan-genotypic near full-length HIV-1 genome (HIV-NFLG) next-generation sequencing (NGS) is key. Our group has established two protocols for HIV-NFLG on the Illumina platform that we aim to compare and, if necessary, modify to find a method with optimized coverage, depth, and subtype inclusivity. Zanini et al. (https://doi.org/10.7554/eLife.11282.001) have developed a method with one-step RT-PCR with six overlapping primer sets, followed by NGS and quality filtering and assembly with in-house methods. Aralaguppe et al. (https://doi.org/10.1016/j.jviromet.2016.07.010) have designed amplification in two fragments, followed by multiplexed NGS and quality control and assembly with Iterative Virus Assembler and VICUNA. Both methods have high coverage per nucleotide and low error rates in amplification and sequencing and can reliably identify SNPs at 1 per cent of the viral population with linkage within the quasispecies. Subtype inclusivity remains a challenge even though both methods show success in amplifying and sequencing subtypes B, C, and the common recombinants 01_AE and 02_AG. Therefore, we aim to evaluate and optimize our NFLG NGS methods on a panel of patient samples that more completely reflects HIV-1 diversity in Sweden. Patient samples from fifty treatment-naïve viremic individuals representing the genotypic HIV-1 panorama in Sweden, including CRFs, are being amplified and sequenced by both protocols. Coverage of the genome, error rate, and possible depth of quasispecies analysis is being evaluated. We will compare number of reads, coverage across the HIV genome, and representation of minor single nucleotide variants as well as subtype inclusivity and impact of plasma RNA levels. To do this we will use an in-house bioinformatic pipeline. The NFLG sequences will also be analyzed with phylogenetic tools for determination of subtypes including CRFs and URFs.
|
|
2. |
- Bennedbaek, M., et al.
(author)
-
Phylogenetic analysis of HIV-1 shows frequent cross-country transmission and local population expansions
- 2021
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 7:2
-
Journal article (peer-reviewed)abstract
- Understanding of pandemics depends on the characterization of pathogen collections from well-defined and demographically diverse cohorts. Since its emergence in Congo almost a century ago, Human Immunodeficiency Virus Type 1 (HIV-1) has geographically spread and genetically diversified into distinct viral subtypes. Phylogenetic analysis can be used to reconstruct the ancestry of the virus to better understand the origin and distribution of subtypes. We sequenced two 3.6-kb amplicons of HIV-1 genomes from 3,197 participants in a clinical trial with consistent and uniform sampling at sites across 35 countries and analyzed our data with another 2,632 genomes that comprehensively reflect the HIV-1 genetic diversity. We used maximum likelihood phylogenetic analysis coupled with geographical information to infer the state of ancestors. The majority of our sequenced genomes (n=2,501) were either pure subtypes (A-D, F, and G) or CRF01_AE. The diversity and distribution of subtypes across geographical regions differed; USA showed the most homogenous subtype population, whereas African samples were most diverse. We delineated transmission of the four most prevalent subtypes in our dataset (A, B, C, and CRF01_AE), and our results suggest both continuous and frequent transmission of HIV-1 over country borders, as well as single transmission events being the seed of endemic population expansions. Overall, we show that coupling of genetic and geographical information of HIV-1 can be used to understand the origin and spread of pandemic pathogens.
|
|
3. |
- Boswell, M. T., et al.
(author)
-
Intrahost evolution of the HIV-2 capsid correlates with progression to AIDS
- 2022
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 8:2
-
Journal article (peer-reviewed)abstract
- HIV-2 infection will progress to AIDS in most patients without treatment, albeit at approximately half the rate of HIV-1 infection. HIV-2 capsid (p26) amino acid polymorphisms are associated with lower viral loads and enhanced processing of T cell epitopes, which may lead to protective Gag-specific T cell responses common in slower progressors. Lower virus evolutionary rates, and positive selection on conserved residues in HIV-2 env have been associated with slower progression to AIDS. In this study we analysed 369 heterochronous HIV-2 p26 sequences from 12 participants with a median age of 30 years at enrolment. CD4% change over time was used to stratify participants into relative faster and slower progressor groups. We analysed p26 sequence diversity evolution, measured site-specific selection pressures and evolutionary rates, and determined if these evolutionary parameters were associated with progression status. Faster progressors had lower CD4% and faster CD4% decline rates. Median pairwise sequence diversity was higher in faster progressors (5.7x10-3 versus 1.4x10-3 base substitutions per site, P<0.001). p26 evolved under negative selection in both groups (dN/dS=0.12). Median virus evolutionary rates were higher in faster than slower progressors – synonymous rates: 4.6x10-3 vs. 2.3x10-3; and nonsynonymous rates: 6.9x10-4 vs. 2.7x10-4 substitutions/site/year, respectively. Virus evolutionary rates correlated negatively with CD4% change rates (ρ = -0.8, P=0.02), but not CD4% level. The signature amino acid at p26 positions 6, 12 and 119 differed between faster (6A, 12I, 119A) and slower (6G, 12V, 119P) progressors. These amino acid positions clustered near to the TRIM5α/p26 hexamer interface surface. p26 evolutionary rates were associated with progression to AIDS and were mostly driven by synonymous substitutions. Nonsynonymous evolutionary rates were an order of magnitude lower than synonymous rates, with limited amino acid sequence evolution over time within hosts. These results indicate HIV-2 p26 may be an attractive therapeutic target.
|
|
4. |
- Chang, Wei-Shan, et al.
(author)
-
Novel hepatitis D-like agents in vertebrates and invertebrates
- 2019
-
In: Virus Evolution. - : OXFORD UNIV PRESS. - 2057-1577. ; 5:2
-
Journal article (peer-reviewed)abstract
- Hepatitis delta virus (HDV) is the smallest known RNA virus, encoding a single protein. Until recently, HDV had only been identified in humans, where it is strongly associated with co-infection with hepatitis B virus (HBV). However, the recent discovery of HDV-like viruses in metagenomic samples from birds and snakes suggests that this virus has a far longer evolutionary history. Herein, using additional meta-transcriptomic data, we show that highly divergent HDV-like viruses are also present in fish, amphibians, and invertebrates, with PCR and Sanger sequencing confirming the presence of the invertebrate HDV-like viruses. Notably, the novel viruses identified here share genomic features characteristic of HDV, such as a circular genome of only approximately 1.7 kb in length, and self-complementary, unbranched rod-like structures. Coiled-coil domains, leucine zippers, conserved residues with essential biological functions, and isoelectronic points similar to those in the human hepatitis delta virus antigens (HDAgs) were also identified in the putative non-human viruses. Importantly, none of these novel HDV-like viruses were associated with hepadnavirus infection, supporting the idea that the HDV-HBV association may be specific to humans. Collectively, these data not only broaden our understanding of the diversity and host range of HDV, but also shed light on its origin and evolutionary history.
|
|
5. |
- de Klerk, A, et al.
(author)
-
Conserved recombination patterns across coronavirus subgenera
- 2022
-
In: Virus evolution. - : Oxford University Press (OUP). - 2057-1577. ; 8:2, s. veac054-
-
Journal article (peer-reviewed)abstract
- Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination break points at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination break points across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination break-point hotspot locations. We find that while the locations of recombination break points are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination break points most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination break-point distributions in coronavirus genomes sampled from nature.
|
|
6. |
|
|
7. |
- Esbjornsson, J., et al.
(author)
-
HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries
- 2016
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 2:1
-
Journal article (peer-reviewed)abstract
- Increased knowledge about HIV-1 transmission dynamics in different transmission groups and geographical regions is fundamental for assessing and designing prevention efforts against HIV-1 spread. Since the first reported cases of HIV infection during the early 1980s, the HIV-1 epidemic in the Nordic countries has been dominated by HIV-1 subtype B and MSM transmission. HIV-1 pol sequences and clinical data of 51 per cent of all newly diagnosed HIV-1 infections in Sweden, Denmark, and Finland in the period 2000-2012 (N = 3,802) were analysed together with a large reference sequence dataset (N = 4,537) by trend analysis and phylogenetics. Analysis of the eight dominating subtypes and CRFs in the Nordic countries (A, B, C, D, G, CRF01_AE, CRF02_AG, and CRF06_cpx) showed that the subtype B proportion decreased while the CRF proportion increased over the study period. A majority (57 per cent) of the Nordic sequences formed transmission clusters, with evidence of mixing both geographically and between transmission groups. Detailed analyses showed multiple occasions of transmissions from MSM to heterosexuals and that active transmission clusters more often involved single than multiple Nordic countries. The strongest geographical link was between Denmark and Sweden. Finally, Denmark had a larger proportion of heterosexual domestic spread of HIV-1 subtype B (75 per cent) compared with Sweden (49 per cent) and Finland (57 per cent). We describe different HIV-1 transmission patterns between countries and transmission groups in a large geographical region. Our results may have implications for public health interventions in targeting HIV-1 transmission networks and identifying where to introduce such interventions.
|
|
8. |
- Gaunt, Michael W., et al.
(author)
-
Widespread interspecific phylogenetic tree incongruence between mosquito-borne and insect-specific flaviviruses at hotspots originally identified in Zika virus
- 2022
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 8:1
-
Journal article (peer-reviewed)abstract
- Intraspecies (homologous) phylogenetic incongruence, or 'tree conflict' between different loci within the same genome of mosquito-borne flaviviruses (MBFV), was first identified in dengue virus (DENV) and subsequently in Japanese encephalitis virus (JEV), St Louis encephalitis virus, and Zika virus (ZIKV). Recently, the first evidence of phylogenetic incongruence between interspecific members of the MBFV was reported in ZIKV and its close relative, Spondweni virus. Uniquely, these hybrid proteomes were derived from four incongruent trees involving an Aedes-associated DENV node (1 tree) and three different Culex-associated flavivirus nodes (3 trees). This analysis has now been extended across a wider spectrum of viruses within the MBFV lineage targeting the breakpoints between phylogenetic incongruent loci originally identified in ZIKV. Interspecies phylogenetic incongruence at these breakpoints was identified in 10 of 50 viruses within the MBFV lineage, representing emergent Aedes and Culex-associated viruses including JEV, West Nile virus, yellow fever virus, and insect-specific viruses. Thus, interspecies phylogenetic incongruence is widespread amongst the flaviviruses and is robustly associated with the specific breakpoints that coincide with the interspecific phylogenetic incongruence previously identified, inferring they are 'hotspots'. The incongruence amongst the emergent MBFV group was restricted to viruses within their respective associated epidemiological boundaries. This MBFV group was RY-coded at the third codon position ('wobble codon') to remove transition saturation. The resulting 'wobble codon' trees presented a single topology for the entire genome that lacked any robust evidence of phylogenetic incongruence between loci. Phylogenetic interspecific incongruence was therefore observed for exactly the same loci between amino acid and the RY-coded 'wobble codon' alignments and this incongruence represented either a major part, or the entire genomes. Maximum likelihood codon analysis revealed positive selection for the incongruent lineages. Positive selection could result in the same locus producing two opposing trees. These analyses for the clinically important MBFV suggest that robust interspecific phylogenetic incongruence resulted from amino acid selection. Convergent or parallel evolutions are evolutionary processes that would explain the observation, whilst interspecific recombination is unlikely.
|
|
9. |
|
|
10. |
- Henningsson, R., et al.
(author)
-
Disseqt-distribution-based modeling of sequence space time dynamics
- 2019
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 5:2
-
Journal article (peer-reviewed)abstract
- Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenningsson/DISSEQT.jl, accessed 23 June 2019) and Synapse (https://www.synapse.org/#!Synapse: syn11425758, accessed 23 June 2019), covering the entire workflow from read alignment to visualization of results. Our pipeline is centered around robust dimension and model reduction algorithms for analysis of genotypic data with additional capabilities for including phenotypic features to explore dynamic genotype-phenotype maps. We illustrate its utility and capacity with examples from evolving RNA virus populations, which present one of the highest degrees of genetic heterogeneity within a given population found in nature. Using our pipeline, we empirically reconstruct the evolutionary trajectories of evolving populations in sequence space and genotype-phenotype fitness landscapes. We show that while sequence space is vastly multidimensional, the relevant genetic space of evolving microbial populations is of intrinsically low dimension. In addition, evolutionary trajectories of these populations can be faithfully monitored to identify the key minority genotypes contributing most to evolution. Finally, we show that empirical fitness landscapes, when reconstructed to include minority variants, can predict phenotype from genotype with high accuracy.
|
|
11. |
- Hostager, R, et al.
(author)
-
Hepatitis C virus genotype 1 and 2 recombinant genomes and the phylogeographic history of the 2k/1b lineage
- 2019
-
In: Virus evolution. - : Oxford University Press (OUP). - 2057-1577. ; 5:2, s. vez041-
-
Journal article (peer-reviewed)abstract
- Recombination is an important driver of genetic diversity, though it is relatively uncommon in hepatitis C virus (HCV). Recent investigation of sequence data acquired from HCV clinical trials produced twenty-one full-genome recombinant viruses belonging to three putative inter-subtype forms 2b/1a, 2b/1b, and 2k/1b. The 2k/1b chimera is the only known HCV circulating recombinant form (CRF), provoking interest in its genetic structure and origin. Discovered in Russia in 1999, 2k/1b cases have since been detected throughout the former Soviet Union, Western Europe, and North America. Although 2k/1b prevalence is highest in the Caucasus mountain region (i.e., Armenia, Azerbaijan, and Georgia), the origin and migration patterns of CRF 2k/1b have remained obscure due to a paucity of available sequences. We assembled an alignment which spans the entire coding region of the HCV genome containing all available 2k/1b sequences (>500 nucleotides; n = 109) sampled in ninteen countries from public databases (102 individuals), additional newly sequenced genomic regions (from 48 of these 102 individuals), unpublished isolates with newly sequenced regions (5 additional individuals), and novel complete genomes (2 additional individuals) generated in this study. Analysis of this expanded dataset reconfirmed the monophyletic origin of 2k/1b with a recombination breakpoint at position 3,187 (95% confidence interval: 3,172–3,202; HCV GT1a reference strain H77). Phylogeography is a valuable tool used to reveal viral migration dynamics. Inference of the timed history of spread in a Bayesian framework identified Russia as the ancestral source of the CRF 2k/1b clade. Further, we found evidence for migration routes leading out of Russia to other former Soviet Republics or countries under the Soviet sphere of influence. These findings suggest an interplay between geopolitics and the historical spread of CRF 2k/1b.
|
|
12. |
- Low, Matthew, et al.
(author)
-
Origins, diversity, and adaptive evolution of DWV in the honey bees of the Azores: the impact of the invasive mite Varroa destructor
- 2024
-
In: Virus evolution. - 2057-1577. ; 10
-
Journal article (peer-reviewed)abstract
- Deformed wing virus (DWV) is a honey bee virus, whose emergence from relative obscurity is driven by the recent host-switch, adaptation, and global dispersal of the ectoparasitic mite Varroa destructor (a highly efficient vector of DWV) to reproduction on honey bees (Apis mellifera). Our study examines how varroa affects the continuing evolution of DWV, using the Azores archipelago, where varroa is present on only three out of the eight Islands, as a natural experimental system for comparing different evolutionary conditions and trajectories. We combined qPCR of 494 honey bee colonies sampled across the archipelago with amplicon deep sequencing to reveal how the DWV genetic landscape is altered by varroa. Two of the varroa-free Islands were also free of DWV, while a further two Islands were intriguingly dominated by the rare DWV-C major variant. The other four Islands, including the three varroa-infested Islands, were dominated by the common DWV-A and DWV-B variants. The varroa-infested Islands had, as expected, an elevated DWV prevalence relative to the uninfested Islands, but not elevated DWV loads, due the relatively high prevalence and loads of DWV-C on the varroa-free Islands. This establishes the Azores as a stable refuge for DWV-C and provides the most convincing evidence to date that at least some major strains of DWV may be capable of not just surviving, but actually thriving in honey bees in the absence of varroa-mediated transmission. We did not detect any change in DWV genetic diversity associated with island varroa status but did find a positive association of DWV diversity with virus load, irrespective of island varroa status.
|
|
13. |
- Martin, DP, et al.
(author)
-
RDP5: a computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets
- 2021
-
In: Virus evolution. - : Oxford University Press (OUP). - 2057-1577. ; 7:1, s. veaa087-
-
Journal article (peer-reviewed)abstract
- For the past 20 years, the recombination detection program (RDP) project has focused on the development of a fast, flexible, and easy to use Windows-based recombination analysis tool. Whereas previous versions of this tool have relied on considerable user-mediated verification of detected recombination events, the latest iteration, RDP5, is automated enough that it can be integrated within analysis pipelines and run without any user input. The main innovation enabling this degree of automation is the implementation of statistical tests to identify recombination signals that could be attributable to evolutionary processes other than recombination. The additional analysis time required for these tests has been offset by algorithmic improvements throughout the program such that, relative to RDP4, RDP5 will still run up to five times faster and be capable of analyzing alignments containing twice as many sequences (up to 5000) that are five times longer (up to 50 million sites). For users wanting to remove signals of recombination from their datasets before using them for downstream phylogenetics-based molecular evolution analyses, RDP5 can disassemble detected recombinant sequences into their constituent parts and output a variety of different recombination-free datasets in an array of different alignment formats. For users that are interested in exploring the recombination history of their datasets, all the manual verification, data management and data visualization components of RDP5 have been extensively updated to minimize the amount of time needed by users to individually verify and refine the program’s interpretation of each of the individual recombination events that it detects.
|
|
14. |
- Martin, DP, et al.
(author)
-
RDP5: a computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets
- 2021
-
In: Virus evolution. - : Oxford University Press (OUP). - 2057-1577. ; 7:1, s. veaa087-
-
Journal article (peer-reviewed)abstract
- For the past 20 years, the recombination detection program (RDP) project has focused on the development of a fast, flexible, and easy to use Windows-based recombination analysis tool. Whereas previous versions of this tool have relied on considerable user-mediated verification of detected recombination events, the latest iteration, RDP5, is automated enough that it can be integrated within analysis pipelines and run without any user input. The main innovation enabling this degree of automation is the implementation of statistical tests to identify recombination signals that could be attributable to evolutionary processes other than recombination. The additional analysis time required for these tests has been offset by algorithmic improvements throughout the program such that, relative to RDP4, RDP5 will still run up to five times faster and be capable of analyzing alignments containing twice as many sequences (up to 5000) that are five times longer (up to 50 million sites). For users wanting to remove signals of recombination from their datasets before using them for downstream phylogenetics-based molecular evolution analyses, RDP5 can disassemble detected recombinant sequences into their constituent parts and output a variety of different recombination-free datasets in an array of different alignment formats. For users that are interested in exploring the recombination history of their datasets, all the manual verification, data management and data visualization components of RDP5 have been extensively updated to minimize the amount of time needed by users to individually verify and refine the program’s interpretation of each of the individual recombination events that it detects.
|
|
15. |
- Nduva, George M., et al.
(author)
-
Quantifying rates of HIV-1 flow between risk groups and geographic locations in Kenya : A country-wide phylogenetic study
- 2022
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 8:1
-
Journal article (peer-reviewed)abstract
- In Kenya, HIV-1 key populations including men having sex with men (MSM), people who inject drugs (PWID) and female sex workers (FSW) are thought to significantly contribute to HIV-1 transmission in the wider, mostly heterosexual (HET) HIV-1 transmission network. However, clear data on HIV-1 transmission dynamics within and between these groups are limited. We aimed to empirically quantify rates of HIV-1 flow between key populations and the HET population, as well as between different geographic regions to determine HIV-1 'hotspots' and their contribution to HIV-1 transmission in Kenya. We used maximum-likelihood phylogenetic and Bayesian inference to analyse 4058 HIV-1 pol sequences (representing 0.3 per cent of the epidemic in Kenya) sampled 1986-2019 from individuals of different risk groups and regions in Kenya. We found 89 per cent within-risk group transmission and 11 per cent mixing between risk groups, cyclic HIV-1 exchange between adjoining geographic provinces and strong evidence of HIV-1 dissemination from (i) West-to-East (i.e. higher-to-lower HIV-1 prevalence regions), and (ii) heterosexual-to-key populations. Low HIV-1 prevalence regions and key populations are sinks rather than major sources of HIV-1 transmission in Kenya. Targeting key populations in Kenya needs to occur concurrently with strengthening interventions in the general epidemic.
|
|
16. |
- Ortiz-Baez, Ayda Susana, et al.
(author)
-
Meta-transcriptomics reveals potential virus transfer between Aedes communis mosquitoes and their parasitic water mites
- 2022
-
In: Virus Evolution. - : Oxford University Press. - 2057-1577. ; 8:2
-
Journal article (peer-reviewed)abstract
- Arthropods harbor a largely undocumented diversity of RNA viruses. Some arthropods, like mosquitoes, can transmit viruses to vertebrates but are themselves parasitized by other arthropod species, such as mites. Very little is known about the viruses of these ectoparasites and how they move through the host-parasite relationship. To address this, we determined the virome of both mosquitoes and the mites that feed on them. The mosquito Aedes communis is an abundant and widely distributed species in Sweden, in northern Europe. These dipterans are commonly parasitized by water mite larvae (Trombidiformes: Mideopsidae) that are hypothesized to impose negative selection pressures on the mosquito by reducing fitness. In turn, viruses are dual-host agents in the mosquito-mite interaction. We determined the RNA virus diversity of mite-free and mite-detached mosquitoes, as well as their parasitic mites, using meta-transcriptomic sequencing. Our results revealed an extensive RNA virus diversity in both mites and mosquitoes, including thirty-seven putative novel RNA viruses that cover a wide taxonomic range. Notably, a high proportion of viruses (20/37) were shared between mites and mosquitoes, while a limited number of viruses were present in a single host. Comparisons of virus composition and abundance suggest potential virus transfer between mosquitoes and mites during their symbiotic interaction. These findings shed light on virome diversity and ecology in the context of arthropod host-parasite-virus relationships.
|
|
17. |
- Pettersson, Mats E., et al.
(author)
-
ZBED6 binding motifs correlate with endogenous retroviruses and Syncytin genes
- 2021
-
In: Virus Evolution. - : Oxford University Press. - 2057-1577. ; 7:1
-
Journal article (peer-reviewed)abstract
- Retroviruses have infiltrated vertebrate germlines for millions of years as inherited endogenous retroviruses (ERVs). Mammalian genomes host large numbers of ERVs and transposable elements (TEs), including retrotransposons and DNA transposons, that contribute to genomic innovation and evolution as coopted genes and regulators of diverse functions. To explore features distinguishing coopted ERVs and TEs from other integrations, we focus on the potential role of ZBED6 and repeated ERV domestication as repurposed Syncytin genes. The placental mammal-specific ZBED6 is a DNA transposon-derived transcription regulator and we demonstrate that its binding motifs are associated with distinct Syncytins and that ZBED6 binding motifs are 2- to 3-fold more frequent in ERVs than in flanking DNA. Our observations suggest that ZBED6 could contribute an extended regulatory role of genomic expression, utilizing ERVs as platforms for genomic innovation and evolution.
|
|
18. |
- Porter, Ashleigh F., et al.
(author)
-
Novel hepaci- and pegi-like viruses in native Australian wildlife and non-human primates
- 2020
-
In: Virus Evolution. - : OXFORD UNIV PRESS. - 2057-1577. ; 6:2
-
Journal article (peer-reviewed)abstract
- The Flaviviridae family of positive-sense RNA viruses contains important pathogens of humans and other animals, including Zika virus, dengue virus, and hepatitis C virus. The Flaviviridae are currently divided into four genera-Hepacivirus, Pegivirus, Pestivirus, and Flavivirus-each with a diverse host range. Members of the genus Hepacivirus are associated with an array of animal species, including humans, non-human primates, other mammalian species, as well as birds and fish, while the closely related pegiviruses have been identified in a variety of mammalian taxa, also including humans. Using a combination of total RNA and whole-genome sequencing we identified four novel hepaci-like viruses and one novel variant of a known hepacivirus in five species of Australian wildlife. The hosts infected comprised native Australian marsupials and birds, as well as a native gecko (Gehyra lauta). From these data we identified a distinct marsupial clade of hepaci-like viruses that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots (Perameles nasuta). Distinct lineages of hepaci-like viruses associated with geckos and birds were also identified. By mining the SRA database we similarly identified three new hepaci-like viruses from avian and primate hosts, as well as two novel pegi-like viruses associated with primates. The phylogenetic history of the hepaci-and pegi-like viruses as a whole, combined with co-phylogenetic analysis, provided support for virus-host co-divergence over the course of vertebrate evolution, although with frequent cross-species virus transmission. Overall, our work highlights the diversity of the Hepacivirus and Pegivirus genera as well as the uncertain phylogenetic distinction between.
|
|
19. |
|
|
20. |
- Wallace, Megan A., et al.
(author)
-
The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe
- 2021
-
In: Virus Evolution. - : Oxford University Press (OUP). - 2057-1577. ; 7:1
-
Journal article (peer-reviewed)abstract
- Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.
|
|
21. |
|
|
22. |
|
|
23. |
- Wille, Michelle, et al.
(author)
-
Evolutionary features of a prolific subtype of avian influenza A virus in European waterfowl
- 2022
-
In: Virus Evolution. - : Oxford University Press. - 2057-1577. ; 8:2
-
Journal article (peer-reviewed)abstract
- Avian influenza A virus (AIV) is ubiquitous in waterfowl and is detected annually at high prevalence in waterfowl during the Northern Hemisphere autumn. Some AIV subtypes are globally common in waterfowl, such as H3N8, H4N6, and H6N2, and are detected in the same populations at a high frequency, annually. In order to investigate genetic features associated to the long-term maintenance of common subtypes in migratory ducks, we sequenced 248 H4 viruses isolated across 8 years (2002-9) from mallards (Anas platyrhynchos) sampled in southeast Sweden. Phylogenetic analyses showed that both H4 and N6 sequences fell into three distinct lineages, structured by year of isolation. Specifically, across the 8 years of the study, we observed lineage replacement, whereby a different HA lineage circulated in the population each year. Analysis of deduced amino acid sequences of the HA lineages illustrated key differences in regions of the globular head of hemagglutinin that overlap with established antigenic sites in homologous hemagglutinin H3, suggesting the possibility of antigenic differences among these HA lineages. Beyond HA, lineage replacement was common to all segments, such that novel genome constellations were detected across years. A dominant genome constellation would rapidly amplify in the duck population, followed by unlinking of gene segments as a result of reassortment within 2-3 weeks following introduction. These data help reveal the evolutionary dynamics exhibited by AIV on both annual and decadal scales in an important reservoir host.
|
|
24. |
- Wille, Michelle, et al.
(author)
-
Where do all the subtypes go? : Temporal dynamics of H8-H12 influenza A viruses in waterfowl
- 2018
-
In: Virus Evolution. - : Oxford University Press. - 2057-1577. ; 4:2
-
Journal article (peer-reviewed)abstract
- Influenza A virus (IAV) is ubiquitous in waterfowl. In the northern hemisphere IAV prevalence is highest during the autumn and coincides with a peak in viral subtype diversity. Although haemagglutinin subtypes H1-H12 are associated with waterfowl hosts, subtypes H8-H12 are detected very infrequently. To better understand the role of waterfowl in the maintenance of these rare subtypes, we sequenced H8-H12 viruses isolated from Mallards (Anas platyrhynchos) from 2002 to 2009. These rare viruses exhibited varying ecological and phylodynamic features. The Eurasian clades of H8 and H12 phylogenies were dominated by waterfowl sequences; mostly viruses sequenced in this study. H11, once believed to be a subtype that infected charadriiformes (shorebirds), exhibited patterns more typical of common virus subtypes. Finally, subtypes H9 and H10, which have maintained lineages in poultry, showed markedly different patterns: H10 was associated with all possible NA subtypes and this drove HA lineage diversity within years. Rare viruses belonging to subtypes H8-H12 were highly reassorted, indicating that these rare subtypes are part of the broader IAV pool. Our results suggest that waterfowl play a role in the maintenance of these rare subtypes, but we recommend additional sampling of non-traditional hosts to better understand the reservoirs of these rare viruses.
|
|
25. |
|
|
26. |
- Zanini, F, et al.
(author)
-
In vivo mutation rates and the landscape of fitness costs of HIV-1
- 2017
-
In: Virus evolution. - : Oxford University Press (OUP). - 2057-1577. ; 3:1, s. vex003-
-
Journal article (peer-reviewed)abstract
- Mutation rates and fitness costs of deleterious mutations are difficult to measure in vivo but essential for a quantitative understanding of evolution. Using whole genome deep sequencing data from longitudinal samples during untreated HIV-1 infection, we estimated mutation rates and fitness costs in HIV-1 from the dynamics of genetic variation. At approximately neutral sites, mutations accumulate with a rate of 1.2 × 10−5 per site per day, in agreement with the rate measured in cell cultures. We estimated the rate from G to A to be the largest, followed by the other transitions C to T, T to C, and A to G, while transversions are less frequent. At other sites, mutations tend to reduce virus replication. We estimated the fitness cost of mutations at every site in the HIV-1 genome using a model of mutation selection balance. About half of all non-synonymous mutations have large fitness costs (>10 percent), while most synonymous mutations have costs <1 percent. The cost of synonymous mutations is especially low in most of pol where we could not detect measurable costs for the majority of synonymous mutations. In contrast, we find high costs for synonymous mutations in important RNA structures and regulatory regions. The intra-patient fitness cost estimates are consistent across multiple patients, indicating that the deleterious part of the fitness landscape is universal and explains a large fraction of global HIV-1 group M diversity.
|
|
27. |
- Zhao, LL, et al.
(author)
-
Phylogenetic estimation of the viral fitness landscape of HIV-1 set-point viral load
- 2022
-
In: Virus evolution. - : Oxford University Press (OUP). - 2057-1577. ; 8:1, s. veac022-
-
Journal article (peer-reviewed)abstract
- Set-point viral load (SPVL), a common measure of human immunodeficiency virus (HIV)-1 virulence, is partially determined by viral genotype. Epidemiological evidence suggests that this viral property has been under stabilising selection, with a typical optimum for the virus between 104 and 105 copies of viral RNA per ml. Here we aimed to detect transmission fitness differences between viruses from individuals with different SPVLs directly from phylogenetic trees inferred from whole-genome sequences. We used the local branching index (LBI) as a proxy for transmission fitness. We found that LBI is more sensitive to differences in infectiousness than to differences in the duration of the infectious state. By analysing subtype-B samples from the Bridging the Evolution and Epidemiology of HIV in Europe project, we inferred a significant positive relationship between SPVL and LBI up to approximately 105 copies/ml, with some evidence for a peak around this value of SPVL. This is evidence of selection against low values of SPVL in HIV-1 subtype-B strains, likely related to lower infectiousness, and perhaps a peak in the transmission fitness in the expected range of SPVL. The less prominent signatures of selection against higher SPVL could be explained by an inherent limit of the method or the deployment of antiretroviral therapy.
|
|