Evolutionarily conserved RNA secondary structures in coding and non-coding sequences at the 3' end of the hepatitis G virus/GB-virus C genome

Equine pegiviruses cause persistent infection of bone marrow and are not associated with hepatitis

Pegiviruses frequently cause persistent infection (as defined by >6 months), but unlike most other Flaviviridae members, no apparent clinical disease. Human pegivirus (HPgV, previously GBV-C) is detectable in 1–4% of healthy individuals and another 5–13% are seropositive. Some evidence for infection of bone marrow and spleen exists. Equine pegivirus 1 (EPgV-1) is not linked to disease, whereas another pegivirus, Theiler’s disease-associated virus (TDAV), was identified in an outbreak of acute serum hepatitis (Theiler’s disease) in horses. Although no subsequent reports link TDAV to disease, any association with hepatitis has not been formally examined. Here, we characterized EPgV-1 and TDAV tropism, sequence diversity, persistence and association with liver disease in horses. Among more than 20 tissue types, we consistently detected high viral loads only in serum, bone marrow and spleen, and viral RNA replication was consistently identified in bone marrow. PBMCs and lymph nodes, but not liver, were sporadically positive. To exclude potential effects of co-infecting agents in experimental infections, we constructed full-length consensus cDNA clones; this was enabled by determination of the complete viral genomes, including a novel TDAV 3’ terminus. Clone derived RNA transcripts were used for direct intrasplenic inoculation of healthy horses. This led to productive infection detectable from week 2–3 and persisting beyond the 28 weeks of study. We did not observe any clinical signs of illness or elevation of circulating liver enzymes. The polyprotein consensus sequences did not change, suggesting that both clones were fully functional. To our knowledge, this is the first successful extrahepatic viral RNA launch and the first robust reverse genetics system for a pegivirus. In conclusion, equine pegiviruses are bone marrow tropic, cause persistent infection in horses, and are not associated with hepatitis. Based on these findings, it may be appropriate to rename the group of TDAV and related viruses as EPgV-2.

Transmissible hepatitis in horses (Theiler’s disease) has been known for 100 years without knowledge of causative infectious agents. Recently, two novel equine pegiviruses (EPgV) were discovered. Whereas EPgV-1 was not associated to disease, the other was identified in an outbreak of acute serum hepatitis and therefore named Theiler’s disease-associated virus (TDAV). This finding was surprising since human and monkey pegiviruses typically cause long-term infection without associated clinical disease. Whereas no subsequent reports link TDAV to disease, the original association to hepatitis has not been formally examined. Here, we studied EPgV-1 and TDAV and found that their natural history of infection in horses were remarkably similar. Examination of various tissues identified the bone marrow as the primary site of replication for both viruses with no evidence of replication in the liver. To exclude potential effects of other infectious agents, we developed molecular full-length clones for EPgV-1 and TDAV and were able to initiate infection in horses using derived synthetic viral genetic material. This demonstrated long-term infection, but no association with hepatitis. These findings call into question the connection between TDAV, liver infection, and hepatitis in horses.

Evolutionary relationship of 5'-untranslated regions among Thai dengue-3 viruses, Bangkok isolates, during 24 year-evolution

OBJECTIVE

To study evolutionary relationship of the 5'untranslated regions (5'UTRs) in low passage dengue3 viruses (DEN3) isolated from hospitalized children with different clinical manifestations in Bangkok during 24 year-evolution (1977-2000) comparing to the DEN3 prototype (H87).

METHODS

The 5'UTRs of these Thai DEN3 and the H87 prototype were amplified by RT-PCR and sequenced. Their multiple sequence alignments were done by Codon Code Aligner v 4.0.4 software and their RNA secondary structures were predicted by MFOLD software. Replication of five Thai DEN3 candidates comparing to the H87 prototype were done in human (HepG2) and the mosquito (C6/36) cell lines.

RESULTS

Among these Thai DEN3, the completely identical sequences of their first 89 nucleotides, their high-order secondary structure of 5'UTRs and three SNPs including the predominant C90T, and two minor SNPs including A109G and A112G were found. The C90T of Thai DEN3, Bangkok isolates was shown predominantly before 1977. Five Thai DEN3 candidates with the predominant C90T were shown to replicate in human (HepG2) and the mosquito (C6/36) cell lines better than the H87 prototype. However, their highly conserved sequences as well as SNPs of the 5'UTR did not appear to correlate with their disease severity in human.

CONCLUSIONS

Our findings highlighted evolutionary relationship of the completely identical 89 nucleotide sequence, the high-order secondary structure and the predominant C90T of the 5'UTR of these Thai DEN3 during 24 year-evolution further suggesting to be their genetic markers and magic targets for future research on antiviral therapy as well as vaccine approaches of Thai DEN3.

Discovery and characterization of distinct simian pegiviruses in three wild African Old World monkey species

Within the Flaviviridae, the recently designated genus Pegivirus has expanded greatly due to new discoveries in bats, horses, and rodents. Here we report the discovery and characterization of three simian pegiviruses (SPgV) that resemble human pegivirus (HPgV) and infect red colobus monkeys (Procolobus tephrosceles), red-tailed guenons (Cercopithecus ascanius) and an olive baboon (Papio anubis). We have designated these viruses SPgVkrc, SPgVkrtg and SPgVkbab, reflecting their host species' common names, which include reference to their location of origin in Kibale National Park, Uganda. SPgVkrc and SPgVkrtg were detected in 47% (28/60) of red colobus and 42% (5/12) red-tailed guenons, respectively, while SPgVkbab infection was observed in 1 of 23 olive baboons tested. Infections were not associated with any apparent disease, despite the generally high viral loads observed for each variant. These viruses were monophyletic and equally divergent from HPgV and pegiviruses previously identified in chimpanzees (SPgVcpz). Overall, the high degree of conservation of genetic features among the novel SPgVs, HPgV and SPgVcpz suggests conservation of function among these closely related viruses. Our study describes the first primate pegiviruses detected in Old World monkeys, expanding the known genetic diversity and host range of pegiviruses and providing insight into the natural history of this genus.

Unmasking the information encoded as structural motifs of viral RNA genomes: a potential antiviral target

RNA viruses show enormous capacity to evolve and adapt to new cellular and molecular contexts, a consequence of mutations arising from errors made by viral RNA-dependent RNA polymerase during replication. Sequence variation must occur, however, without compromising functions essential for the completion of the viral cycle. RNA viruses are safeguarded in this respect by their genome carrying conserved information that does not code only for proteins but also for the formation of structurally conserved RNA domains that directly perform these critical functions. Functional RNA domains can interact with other regions of the viral genome and/or proteins to direct viral translation, replication and encapsidation. They are therefore potential targets for novel therapeutic strategies. This review summarises our knowledge of the functional RNA domains of human RNA viruses and examines the achievements made in the design of antiviral compounds that interfere with their folding and therefore their function.

Deep sequencing identifies two genotypes and high viral genetic diversity of human pegivirus (GB virus C) in rural Ugandan patients

Human pegivirus (HPgV), formerly 'GB virus C' or 'hepatitis G virus', is a member of the genus Flavivirus (Flaviviridae) that has garnered significant attention due to its inhibition of HIV, including slowing disease progression and prolonging survival in HIV-infected patients. Currently, there are six proposed HPgV genotypes that have roughly distinct geographical distributions. Genotypes 2 and 3 are the most comprehensively characterized, whereas those genotypes occurring on the African continent, where HPgV prevalence is highest, are less well studied. Using deep sequencing methods, we identified complete coding HPgV sequences in four of 28 patients (14.3%) in rural Uganda, east Africa. One of these sequences corresponds to genotype 1 and is the first complete genome of this genotype from east Africa. The remaining three sequences correspond to genotype 5, a genotype that was previously considered exclusively South African. All four positive samples were collected within a geographical area of less than 25 km(2), showing that multiple HPgV genotypes co-circulate in this area. Analysis of intra-host viral genetic diversity revealed that total single-nucleotide polymorphism frequency was approximately tenfold lower in HPgV than in hepatitis C virus. Finally, one patient was co-infected with HPgV and HIV, which, in combination with the high prevalence of HIV, suggests that this region would be a useful locale to study the interactions and co-evolution of these viruses.

Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA

BACKGROUND

Epstein-Barr virus (EBV) is a human herpesvirus implicated in cancer and autoimmune disorders. Little is known concerning the roles of RNA structure in this important human pathogen. This study provides the first comprehensive genome-wide survey of RNA and RNA structure in EBV.

RESULTS

Novel EBV RNAs and RNA structures were identified by computational modeling and RNA-Seq analyses of EBV. Scans of the genomic sequences of four EBV strains (EBV-1, EBV-2, GD1, and GD2) and of the closely related Macacine herpesvirus 4 using the RNAz program discovered 265 regions with high probability of forming conserved RNA structures. Secondary structure models are proposed for these regions based on a combination of free energy minimization and comparative sequence analysis. The analysis of RNA-Seq data uncovered the first observation of a stable intronic sequence RNA (sisRNA) in EBV. The abundance of this sisRNA rivals that of the well-known and highly expressed EBV-encoded non-coding RNAs (EBERs).

CONCLUSION

This work identifies regions of the EBV genome likely to generate functional RNAs and RNA structures, provides structural models for these regions, and discusses potential functions suggested by the modeled structures. Enhanced understanding of the EBV transcriptome will guide future experimental analyses of the discovered RNAs and RNA structures.

A comparison of genotype-phenotype maps for RNA and proteins

The relationship between the genotype (sequence) and the phenotype (structure) of macromolecules affects their ability to evolve new structures and functions. We here compare the genotype space organization of proteins and RNA molecules to identify differences that may affect this ability. To this end, we computationally study the genotype-phenotype relationship for short RNA and lattice proteins of a reduced monomer alphabet size, to make exhaustive analysis and direct comparison of their genotype spaces feasible. We find that many fewer protein molecules than RNA molecules fold, but they fold into many more structures than RNA. In consequence, protein phenotypes have smaller genotype networks whose member genotypes tend to be more similar than for RNA phenotypes. Neighborhoods in sequence space of a given radius around an RNA molecule contain more novel structures than for protein molecules. We compare this property to evidence from natural RNA and protein molecules, and conclude that RNA genotype space may be more conducive to the evolution of new structure phenotypes.

GB virus C (GBV-C) evolutionary patterns revealed by analyses of reference genomes, E2 and NS5B sequences amplified from viral strains circulating in the Lisbon area (Portugal)

GBV-C is a non-pathogenic virus that is largely dispersed in different human populations. The phylogenetic analysis of the 5'-untranslated region (5'UTR) of the GBV-C genome has led to the segregation of viral strains into six genotypes, but incongruent results are frequently obtained depending on the genome region analyzed. In this report, different phylogenetic approaches and multivariate statistics were combined to disclose evolutionary patterns that contribute to shape GBV-C evolution. The data here presented indicate: (i) that the phylogenetic noise was mostly determined by the size of the analyzed sequence, rather than by its position on the viral genome; (ii) that most genomic segments in the coding sequence seemed to evolve under a similar evolution model, which was different from that which best fits the 5'UTR, with overall large heterogeneity of rate change across the sequence; (iii) that due to saturation of transversions occurring in the 5'UTR at genetic distances <0.10, care should be taken in drawing conclusions about the tree topologies involving the deeper branches, especially when using distance-based methods; (iv) that a non-uniform distribution of InSi and dS occurs over the viral ORF highlighting regions of the viral genome with remarkably low levels of silent substitutions, and implying that the observed differences may contribute to the detected phylogenetic incongruences; and finally (v) that genetic recombination clearly impacts the GBV-C evolution extensively, this being shown for both reference genomes and NS5B GBV-C sequences amplified from Portuguese residents.

Predicted RNA secondary structures for the conserved regions in dengue virus

Dengue fever, dengue hemorrhagic fever and dengue shock syndrome are the prevalent mosquito borne viral infections worldwide. The dengue virus belongs to the genus flavivirus with conserved RNA domains peptidase_S7 and dexHc among its members. The secondary structures for RNA domains peptidase_S7 and DexHc are hence predicted and discussed with other known viral RNA structures to glean structural insights through comparison.

Computational modeling analyses of RNA secondary structures and phylogenetic inference of evolutionary conserved 5S rRNA in the prokaryotes

Bacteria are unicellular, ubiquitous microorganisms which grow on soil, acidic hot springs, radioactive wastes, etc. The genome of bacteria constitutes species specific conserved region. The 5S rRNA is one of the most conserved region determined in each bacteria and the size ranges between 110 and 148 bp. On this basis phylogenetic study of 37 bacterial strains was done which results in formation of seven clades and furthermore RNA secondary structure from each clade was made. The lowest free energy (delta G) of the 5S rRNA may divulge the most primitive bacteria and slow changes occurs throughout the evolution whereas the higher free energy indicates less stability during the evolution. The RNA secondary structure may provide new insights to understand bacteria evolution and stability.

Efficient algorithms for probing the RNA mutation landscape

The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Although secondary structure predictions from thermodynamics-based algorithms are not as accurate as methods employing comparative genomics, the former methods are the only available tools to investigate novel RNAs, such as the many RNAs of unknown function recently reported by the ENCODE consortium. In this paper, we generalize the McCaskill partition function algorithm to sum over the grand canonical ensemble of all secondary structures of all mutants of the given sequence. Specifically, our new program, RNAmutants, simultaneously computes for each integer k the minimum free energy structure MFE(k) and the partition function Z(k) over all secondary structures of all k-point mutants, even allowing the user to specify certain positions required not to mutate and certain positions required to base-pair or remain unpaired. This technically important extension allows us to study the resilience of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. We have successfully applied RNAmutants to investigate deleterious mutations (mutations that radically modify the secondary structure) in the Hepatitis C virus cis-acting replication element and to evaluate the evolutionary pressure applied on different regions of the HIV trans-activation response element. In particular, we show qualitative agreement between published Hepatitis C and HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work also predicts other deleterious mutations, which could be verified experimentally. Finally, we provide evidence that the 3′ UTR of the GB RNA virus C has been optimized to preserve evolutionarily conserved stem regions from a deleterious effect of pointwise mutations. We hope that there will be long-term potential applications of RNAmutants in de novo RNA design and drug design against RNA viruses. This work also suggests potential applications for large-scale exploration of the RNA sequence-structure network. Binary distributions are available at http://RNAmutants.csail.mit.edu/.

Evolution is a central concept in biology. This phenomenon can be observed at all levels of the organization of life—from single molecules to multicellular organisms. Here, we focus our attention on the implication of evolution at the level of nucleic acid sequences. In this context, RNA sequences presumably have been optimized by evolution to achieve specific functions. These functions are supported by a structure that can be determined using thermodynamics-based models and energy minimization techniques. In this work, we develop efficient algorithms for predicting energetically favorable mutations and study their impact on the stability of the structure. We use these techniques to reveal sequences under evolutionary pressure and design new methods to predict lethal mutations. Applications of our tool lead to a better understanding of the mutational process of some key regulatory elements of two important pathogenic RNA viruses—human immunodeficiency virus and hepatitis C virus.

Hepatitis C virus population dynamics during infection

Hepatitis C virus (HCV) behaves as an evolving viral quasispecies in its continuously changing environment. The study of HCV quasispecies population dynamics in experimental models and infected patients can provide useful information on factors involved in the HCV life cycle and pathogenicity. HCV quasispecies variability also has therapeutic implications, as the continuous generation and selection of fitter or truly resistant variants can allow the virus to escape control by antiviral drugs.

Conserved RNA secondary structures in Flaviviridae genomes

Presented here is a comprehensive computational survey of evolutionarily conserved secondary structure motifs in the genomic RNAs of the family Flaviviridae: This virus family consists of the three genera Flavivirus, Pestivirus and Hepacivirus and the group of GB virus C/hepatitis G virus with a currently uncertain taxonomic classification. Based on the control of replication and translation, two subgroups were considered separately: the genus Flavivirus, with its type I cap structure at the 5' untranslated region (UTR) and a highly structured 3' UTR, and the remaining three groups, which exhibit translation control by means of an internal ribosomal entry site (IRES) in the 5' UTR and a much shorter less-structured 3' UTR. The main findings of this survey are strong hints for the possibility of genome cyclization in hepatitis C virus and GB virus C/hepatitis G virus in addition to the flaviviruses; a surprisingly large number of conserved RNA motifs in the coding regions; and a lower level of detailed structural conservation in the IRES and 3' UTR motifs than reported in the literature. An electronic atlas organizes the information on the more than 150 conserved, and therefore putatively functional, RNA secondary structure elements.

Contributions of vesicular stomatitis virus to the understanding of RNA virus evolution

Vesicular stomatitis virus has been a preferred system to study evolution for several decades. New approaches to antiviral treatment, such as lethal mutagenesis, stem from investigations done with VSV. Recent work has shed new light in the way we view neutrality, a fundamental concept to understand the evolutionary history of RNA viruses.

Hepatitis C virus genetic variability: pathogenic and clinical implications

HCV is variable because of the properties of the viral RdRp, high levels of replication, and large population sizes. The Darwinian evolution of HCV has been characterized by the emergence of the HCV genotypes, including six main types and a large number of subtypes. The study of HCV genotype epidemiology provides useful information on the worldwide HCV epidemics. The HCV genotype is an important predictor of the response to IFN-alpha-based antiviral therapy, and genotype determination is currently used to tailor treatment indications. In addition, HCV circulates and behaves in infected individuals as mixtures of closely related but distinct viral populations referred to as quasispecies. This particular nature of the virus influences its transmission, the pathogenesis of liver disease and extra-hepatic manifestations, and the outcome during and after antiviral therapy or after transplantation for HCV-related end-stage liver disease. Further studies are needed to understand better the implications of HCV quasispecies diversity in the pathophysiology of HCV infection.

Reconstructing the origins of human hepatitis viruses

Infections with hepatitis B and C viruses (HBV, HCV) are widespread in human populations throughout the world, and are major causes of chronic liver disease and liver cancer. HBV, HCV and the related hepatitis G virus or GB virus C (referred to here as HGV/GBV-C) are capable of establishing persistent, frequently lifelong infections characterized by high levels of continuous replication. All three viruses show substantial genetic heterogeneity, which has allowed each to be classified into a number of distinct genotypes that have different geographical distributions and associations with different risk groups for infection. Information on their past transmission and epidemiology might be obtained by estimation of the time of divergence of the different genotypes of HCV, HBV and HGV/GBV-C using knowledge of their rates of sequence change. While information on the latter is limited to short observation periods and is therefore subject to considerable error and uncertainty, the relatively recent times of origin for genotype of each virus predicted by this method (HCV, 500-2000 years; HBV, 3000 years; HGV/GBV-C, 200 years) are quite incompatible with their epidemiological distributions in human populations. They also cannot easily be reconciled with the recent evidence for species-associated variants of HBV and HGV/GBV-C in a range of non-human primates. The apparent conservatism of viruses over long periods implied by their epidemiological distributions instead suggests that nucleotide sequence change may be subject to constraints peculiar to viruses with single-stranded genomes, or with overlapping reading frames that defy attempts to reconstruct evolution according to the principles of the 'molecular clock'. Large population sizes and intense selection pressures that optimize fitness may be additional factors that set virus evolution apart from that of their hosts.

The origin and evolution of hepatitis viruses in humans

The spread and origins of hepatitis C virus (HCV) in human populations have been the subject of extensive investigations, not least because of the importance this information would provide in predicting clinical outcomes and controlling spread of HCV in the future. However, in the absence of historical and archaeological records of infection, the evolution of HCV and other human hepatitis viruses can only be inferred indirectly from their epidemiology and by genetic analysis of contemporary virus populations. Some information on the history of the latter may be obtained by dating the time of divergence of various genotypes of HCV, hepatitis B virus (HBV) and the non-pathogenic hepatitis G virus (HGV)/GB virus-C (GBV-C). However, the relatively recent times predicted for the origin of these viruses fit poorly with their epidemiological distributions and the recent evidence for species-associated variants of HBV and HGV/GBV-C in a wide range of non-human primates. The apparent conservatism of viruses over long periods implied by these latter observations may be the result of constraints on sequence change peculiar to viruses with single-stranded genomes, or with overlapping reading frames. Large population sizes and intense selection pressures that optimize fitness may be the factors that set virus evolution apart from that of their hosts.