Dataset Information

Extensive C->U transition biases in the genomes of a wide range of mammalian RNA viruses; potential associations with transcriptional mutations, damage- or host-mediated editing of viral RNA.

ABSTRACT: The rapid evolution of RNA viruses has been long considered to result from a combination of high copying error frequencies during RNA replication, short generation times and the consequent extensive fixation of neutral or adaptive changes over short periods. While both the identities and sites of mutations are typically modelled as being random, recent investigations of sequence diversity of SARS coronavirus 2 (SARS-CoV-2) have identified a preponderance of C->U transitions, proposed to be driven by an APOBEC-like RNA editing process. The current study investigated whether this phenomenon could be observed in datasets of other RNA viruses. Using a 5% divergence filter to infer directionality, 18 from 36 datasets of aligned coding region sequences from a diverse range of mammalian RNA viruses (including Picornaviridae, Flaviviridae, Matonaviridae, Caliciviridae and Coronaviridae) showed a >2-fold base composition normalised excess of C->U transitions compared to U->C (range 2.1x-7.5x), with a consistently observed favoured 5' U upstream context. The presence of genome scale RNA secondary structure (GORS) was the only other genomic or structural parameter significantly associated with C->U/U->C transition asymmetries by multivariable analysis (ANOVA), potentially reflecting RNA structure dependence of sites targeted for C->U mutations. Using the association index metric, C->U changes were specifically over-represented at phylogenetically uninformative sites, potentially paralleling extensive homoplasy of this transition reported in SARS-CoV-2. Although mechanisms remain to be functionally characterised, excess C->U substitutions accounted for 11-14% of standing sequence variability of structured viruses and may therefore represent a potent driver of their sequence diversification and longer-term evolution.

SUBMITTER: Simmonds P

PROVIDER: S-EPMC8195396 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Understanding the microbiome of ticks in Australia is of considerable interest given the ongoing debate over whether Lyme disease and its causative agent, the bacterium Borrelia burgdorferisensu lato, are present in Australia. The diversity of bacteria infecting Australian ticks has been studied using both culture- and metagenomics-based techniques. However, little is known about the virome of Australian ticks, including whether this includes viruses with the potential to infect mammals. We used a meta-transcriptomics approach to reveal the diversity and evolution of viruses from Australian ticks collected from two locations on the central east coast of Australia, including metropolitan Sydney. From this we identified 19 novel RNA viruses belonging to 12 families, as well as 1 previously described RNA virus. The majority of these viruses were related to arthropod-associated viruses, suggesting that they do not utilize mammalian hosts. However, two novel viruses discovered in ticks feeding on bandicoot marsupials clustered closely within the mammal-associated hepacivirus and pestivirus groups (family Flaviviridae). Another bandicoot tick yielded a novel coltivirus (family Reoviridae), a group of largely tick-associated viruses containing the known human pathogen Colorado tick fever virus and its relative, Eyach virus. Importantly, our transcriptomic data provided no evidence for the presence of B. burgdorferisensu lato in any tick sample, providing further evidence against the presence of Lyme disease in Australia. In sum, this study reveals that Australian ticks harbor a diverse virome, including some viruses that merit additional screening in the context of emerging infectious disease.IMPORTANCE Each year a growing number of individuals along the east coast of Australia experience debilitating disease following tick bites. As there is no evidence for the presence of the causative agent of Lyme disease, Borrelia burgdorferisensu lato, in Australian ticks, the etiological basis of this disease syndrome remains controversial. To characterize the viruses associated with Australian ticks, particularly those that might be associated with mammalian infection, we performed unbiased RNA sequencing on 146 ticks collected across two locations along the coast of New South Wales, Australia. This revealed 19 novel RNA viruses from a diverse set of families. Notably, three of these viruses clustered with known mammalian viruses, including a novel coltivirus that was related to the human pathogen Colorado tick fever virus.

Project description:CpG and UpA dinucleotides are under-represented in vertebrate genomes, whereas most invertebrates only show a bias against UpA. RNA viruses are thought to have evolved genomes that resemble the dinucleotide composition of their hosts, possibly to avoid restriction by the zinc-finger antiviral protein (ZAP). By performing a comprehensive analysis of RNA viruses, we show that, whereas UpA dinucleotides are similarly under-represented irrespective of viral genome composition or host, important differences are observed for CpG. The tendency for vertebrate-infecting viruses to have stronger CpG bias than invertebrate-infecting viruses is not universal. Rather, it is mainly driven by single-stranded (ss) RNA(+) viruses. Conversely, ssRNA(-) viruses have a dinucleotide composition that is unrelated to the host clade. Also, these viruses, especially those in the order Bunyavirales, are extremely CpG-depleted. By focusing on specific viral families, we also show that, even for vertebrate ssRNA(+) viruses, ZAP is unlikely to be a driver of CpG depletion. Consistently, CpG dinucleotides tend to be preferentially depleted in A/U-rich contexts in both vertebrate- and invertebrate-infecting viruses. Finally, within the same viral genomes, individual viral open reading frames (ORFs) can display different CpG content. Analysis of SARS-CoV-2 revealed a remarkable depletion of CpG dinucleotides in ORF1ab and S, but not in N and M. Thus, these results do not support the view that an adaptive shift for CpG depletion in the SARS-CoV-2 lineage occurred as an innate immunity evasion strategy. Our data provide a better understanding of viral evolution and inform approaches based on the modulation of CpG to generate attenuated viruses.IMPORTANCEAkin to a molecular signature, dinucleotide composition can be exploited by the zinc-finger antiviral protein (ZAP) to restrict CpG-rich (and UpA-rich) RNA viruses. ZAP evolved in tetrapods, and it is not encoded by invertebrates and fish. Because a systematic analysis is missing, we analyzed the genomes of RNA viruses that infect vertebrates or invertebrates. We show that vertebrate single-stranded (ss) RNA(+) viruses and, to a lesser extent, double-stranded RNA viruses tend to have stronger CpG bias than invertebrate viruses. Conversely, ssRNA(-) viruses have similar dinucleotide composition whether they infect vertebrates or invertebrates. Analysis of ssRNA(+) viruses that infect mammals, reptiles, and fish indicated that ZAP is unlikely to be a major driver of CpG depletion. We also show that, compared to other coronaviruses, the genome of SARS-CoV-2 is not homogeneously CpG-depleted. Our study provides new insights into virus evolution and strategies for recoding RNA virus genomes.

Project description:BACKGROUND:The family Mesoniviridae (order Nidovirales) comprises of a group of positive-sense, single-stranded RNA ([+]ssRNA) viruses isolated from mosquitoes. FINDINGS:Thirteen novel insect-specific virus isolates were obtained from mosquitoes collected in Indonesia, Thailand and the USA. By electron microscopy, the virions appeared as spherical particles with a diameter of ~50 nm. Their 20,129 nt to 20,777 nt genomes consist of positive-sense, single-stranded RNA with a poly-A tail. Four isolates from Houston, Texas, and one isolate from Java, Indonesia, were identified as variants of the species Alphamesonivirus-1 which also includes Nam Dinh virus (NDiV) from Vietnam and Cavally virus (CavV) from Côte d'Ivoire. The eight other isolates were identified as variants of three new mesoniviruses, based on genome organization and pairwise evolutionary distances: Karang Sari virus (KSaV) from Java, Bontag Baru virus (BBaV) from Java and Kalimantan, and Kamphaeng Phet virus (KPhV) from Thailand. In comparison with NDiV, the three new mesoniviruses each contained a long insertion (180 - 588 nt) of unknown function in the 5' region of ORF1a, which accounted for much of the difference in genome size. The insertions contained various short imperfect repeats and may have arisen by recombination or sequence duplication. CONCLUSIONS:In summary, based on their genome organizations and phylogenetic relationships, thirteen new viruses were identified as members of the family Mesoniviridae, order Nidovirales. Species demarcation criteria employed previously for mesoniviruses would place five of these isolates in the same species as NDiV and CavV (Alphamesonivirus-1) and the other eight isolates would represent three new mesonivirus species (Alphamesonivirus-5, Alphamesonivirus-6 and Alphamesonivirus-7). The observed spatiotemporal distribution over widespread geographic regions and broad species host range in mosquitoes suggests that mesoniviruses may be common in mosquito populations worldwide.

Dataset Information

Extensive C->U transition biases in the genomes of a wide range of mammalian RNA viruses; potential associations with transcriptional mutations, damage- or host-mediated editing of viral RNA.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets