The evolution of genome compression and genomic novelty in RNA viruses

The application of single-molecule optical tweezers to study disease-related structural dynamics in RNA

RNA, a dynamic and flexible molecule with intricate three-dimensional structures, has myriad functions in disease development. Traditional methods, such as X-ray crystallography and nuclear magnetic resonance, face limitations in capturing real-time, single-molecule dynamics crucial for understanding RNA function. This review explores the transformative potential of single-molecule force spectroscopy using optical tweezers, showcasing its capability to directly probe time-dependent structural rearrangements of individual RNA molecules. Optical tweezers offer versatility in exploring diverse conditions, with the potential to provide insights into how environmental changes, ligands and RNA-binding proteins impact RNA behaviour. By enabling real-time observations of large-scale structural dynamics, optical tweezers emerge as an invaluable tool for advancing our comprehension of RNA structure and function. Here, we showcase their application in elucidating the dynamics of RNA elements in virology, such as the pseudoknot governing ribosomal frameshifting in SARS-CoV-2.

Hidden evolutionary constraints dictate the retention of coronavirus accessory genes

Coronaviruses exhibit many mechanisms of genetic innovation (1-5), including the acquisition of accessory genes that originate by capture of cellular genes or through duplication of existing viral genes (6,7). Accessory genes influence viral host range and cellular tropism, but little is known about how selection acts on these variable regions of virus genomes. We used experimental evolution of mouse hepatitis virus (MHV) encoding a cellular AKAP7 phosphodiesterase and an inactive native phosphodiesterase, NS2 (ref 8) to simulate the capture of a host gene and analyze its evolution. After courses of serial infection, the gene encoding inactive NS2, ORF2, unexpectedly remained intact, suggesting it is under cryptic constraint uncoupled from the function of NS2. In contrast, AKAP7 was retained under strong selection but rapidly lost under relaxed selection. Guided by the retention of ORF2 and similar patterns in related betacoronaviruses, we analyzed ORF8 of SARS-CoV-2, which arose via gene duplication6 and contains premature stop codons in several globally successful lineages. As with MHV ORF2, the coding-defective SARS-CoV-2 ORF8 gene remains largely intact, mirroring patterns observed during MHV experimental evolution, challenging assumptions on the dynamics of gene loss in virus genomes and extending these findings to viruses currently adapting to humans.

Insights into the genetic variability and evolutionary dynamics of tomato spotted wilt orthotospovirus in China

BACKGROUND

Viral diseases are posing threat to annual production and quality of tobacco in China. Recently, tomato spotted wilt orthotospovirus (TSWV) has been reported to infect three major crops including tobacco. Current study was aimed to investigate the population dynamics and molecular diversity of the TSWV. In the current study, to assess and identify the prevalence and evolutionary history of TSWV in tobacco crops in China, full-length genome sequences of TSWV isolates from tobacco, were identified and analyzed.

METHODS

After trimming and validation, sequences of new isolates were submitted to GenBank. We identified the full-length genomes of ten TSWV isolates, infecting tobacco plants from various regions of China. Besides these, six isolates were partially sequenced. Phylogenetic analysis was performed to assess the relativeness of newly identified sequences and corresponding sequences from GenBank. Recombination and population dynamics analysis was performed using RDP4, RAT, and statistical estimation. Reassortment analysis was performed using MegaX software.

RESULTS

Phylogenetic analysis of 41 newly identified sequences, depicted that the majority of the Chinese isolates have separate placement in the tree. RDP4 software predicted that RNA M of newly reported isolate YNKM-2 had a recombinant region spanning from 3111 to 3811 bp. The indication of parental sequences (YNKMXD and YNHHKY) from newly identified isolates, revealed the conservation of local TSWV population. Genetic diversity and population dynamics analysis also support the same trend. RNA M was highlighted to be more capable of mutating or evolving as revealed by data obtained from RDP4, RAT, population dynamics, and phylogenetic analyses. Reassortment analysis revealed that it might have happened in L segment of TSWV isolate YNKMXD (reported herein).

CONCLUSION

Taken together, this is the first detailed study revealing the pattern of TWSV genetic diversity, and population dynamics helping to better understand the ability of this pathogen to drastically reduce the tobacco production in China. Also, this is a valuable addition to the existing worldwide profile of TSWV, especially in China, where a few studies related to TSWV have been reported including only one complete genome of this virus isolated from tobacco plants.

Overlaps Between CDS Regions of Protein-Coding Genes in the Human Genome: A Case Study on the NR1D1-THRA Gene Pair

For several decades, it has been known that a substantial number of genes within human DNA exhibit overlap; however, the biological and evolutionary significance of these overlaps remain poorly understood. This study focused on investigating specific instances of overlap where the overlapping DNA region encompasses the coding DNA sequences (CDSs) of protein-coding genes. The results revealed that proteins encoded by overlapping CDSs exhibit greater disorder than those from nonoverlapping CDSs. Additionally, these DNA regions were identified as GC-rich. This could be partially attributed to the absence of stop codons from two distinct reading frames rather than one. Furthermore, these regions were found to harbour fewer single-nucleotide polymorphism (SNP) sites, possibly due to constraints arising from the overlapping state where mutations could affect two genes simultaneously.While elucidating these properties, the NR1D1-THRA gene pair emerged as an exceptional case with highly structured proteins and a distinctly conserved sequence across eutherian mammals. Both NR1D1 and THRA are nuclear receptors lacking a ligand-binding domain at their C-terminus, which is the region where these gene pairs overlap. The NR1D1 gene is involved in the regulation of circadian rhythm, while the THRA gene encodes a thyroid hormone receptor, and both play crucial roles in various physiological processes. This study suggests that, in addition to their well-established functions, the specifically overlapping CDS regions of these genes may encode protein segments with additional, yet undiscovered, biological roles.

Raman Fingerprints of SARS-CoV-2 Omicron Subvariants: Molecular Roots of Virological Characteristics and Evolutionary Directions

The latest RNA genomic mutation of SARS-CoV-2 virus, termed the Omicron variant, has generated a stream of highly contagious and antibody-resistant strains, which in turn led to classifying Omicron as a variant of concern. We systematically collected Raman spectra from six Omicron subvariants available in Japan (i.e., BA.1.18, BA.2, BA.4, BA.5, XE, and BA.2.75) and applied machine-learning algorithms to decrypt their structural characteristics at the molecular scale. Unique Raman fingerprints of sulfur-containing amino acid rotamers, RNA purines and pyrimidines, tyrosine phenol ring configurations, and secondary protein structures clearly differentiated the six Omicron subvariants. These spectral characteristics, which were linked to infectiousness, transmissibility, and propensity for immune evasion, revealed evolutionary motifs to be compared with the outputs of genomic studies. The availability of a Raman "metabolomic snapshot", which was then translated into a barcode to enable a prompt subvariant identification, opened the way to rationalize in real-time SARS-CoV-2 activity and variability. As a proof of concept, we applied the Raman barcode procedure to a nasal swab sample retrieved from a SARS-CoV-2 patient and identified its Omicron subvariant by coupling a commercially available magnetic bead technology with our newly developed Raman analyses.

The panzootic potential of SARS-CoV-2

Each year, SARS-CoV-2 is infecting an increasingly unprecedented number of species. In the present article, we combine mammalian phylogeny with the genetic characteristics of isolates found in mammals to elaborate on the host-range potential of SARS-CoV-2. Infections in nonhuman mammals mirror those of contemporary viral strains circulating in humans, although, in certain species, extensive viral circulation has led to unique genetic signatures. As in other recent studies, we found that the conservation of the ACE2 receptor cannot be considered the sole major determinant of susceptibility. However, we are able to identify major clades and families as candidates for increased surveillance. On the basis of our findings, we argue that the use of the term panzootic could be a more appropriate term than pandemic to describe the ongoing scenario. This term better captures the magnitude of the SARS-CoV-2 host range and would hopefully inspire inclusive policy actions, including systematic screenings, that could better support the management of this worldwide event.

A conserved Y-shaped RNA structure in the 3'UTR of chikungunya virus genome as a host-specialized element that modulates viral replication and evolution

RNA viral genomes compact information into functional RNA structures. Here, using chikungunya virus as a model, we investigated the structural requirements of conserved RNA elements in the 3' untranslated region (3'UTR) for viral replication in mosquito and mammalian cells. Using structural predictions and co-variation analysis, we identified a highly stable and conserved Y-shaped structure (SLY) at the end of the 3'UTR that is duplicated in the Asian lineage. Functional studies with mutant viruses showed that the SLY has host-specific functions during viral replication and evolution. The SLY positively modulates viral replication in mosquito cells but has the opposite effect in mammalian cells. Additional structural/functional analyses showed that maintaining the Y-shaped fold and specific nucleotides in the loop are critical for full SLY functionality and optimal viral replication in mosquito cells. Experimental adaptation of viruses with duplicated SLYs to mammalian cells resulted in the generation of heterogeneous viral populations comprising variants with diverse 3'UTRs, contrasting with the homogeneous populations from viruses without SLY copies. Altogether, our findings constitute the first evidence of an RNA secondary structure in the 3'UTR of chikungunya virus genome that plays host-dependent functions.

A Glimpse on the Evolution of RNA Viruses: Implications and Lessons from SARS-CoV-2

RNA viruses are characterised by extremely high genetic variability due to fast replication, large population size, low fidelity, and (usually) a lack of proofreading mechanisms of RNA polymerases leading to high mutation rates. Furthermore, viral recombination and reassortment may act as a significant evolutionary force among viruses contributing to greater genetic diversity than obtainable by mutation alone. The above-mentioned properties allow for the rapid evolution of RNA viruses, which may result in difficulties in viral eradication, changes in virulence and pathogenicity, and lead to events such as cross-species transmissions, which are matters of great interest in the light of current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemics. In this review, we aim to explore the molecular mechanisms of the variability of viral RNA genomes, emphasising the evolutionary trajectory of SARS-CoV-2 and its variants. Furthermore, the causes and consequences of coronavirus variation are explored, along with theories on the origin of human coronaviruses and features of emergent RNA viruses in general. Finally, we summarise the current knowledge on the circulating variants of concern and highlight the many unknowns regarding SARS-CoV-2 pathogenesis.

Analysis of insertions and extensions in the functional evolution of the ribonucleotide reductase family

Ribonucleotide reductases (RNRs) are used by all free-living organisms and many viruses to catalyze an essential step in the de novo biosynthesis of DNA precursors. RNRs are remarkably diverse by primary sequence and cofactor requirement, while sharing a conserved fold and radical-based mechanism for nucleotide reduction. In this work, we expand on our recent phylogenetic inference of the entire RNR family and describe the evolutionarily relatedness of insertions and extensions around the structurally homologous catalytic barrel. Using evo-velocity and sequence similarity network (SSN) analyses, we show that the N-terminal regulatory motif known as the ATP-cone domain was likely inherited from an ancestral RNR. By combining SSN analysis with AlphaFold2 predictions, we also show that the C-terminal extensions of class II RNRs can contain folded domains that share homology with an Fe-S cluster assembly protein. Finally, using sequence analysis and AlphaFold2, we show that the sequence motif of a catalytically essential insertion known as the finger loop is tightly coupled to the catalytic mechanism. Based on these results, we propose an evolutionary model for the diversification of the RNR family.

Expansion of viral genomes with viral protein genome linked copies

Virus protein-linked genome (VPg) proteins are required for replication. VPgs are duplicated in a subset of RNA viruses however their roles are not fully understood and the extent of viral genomes containing VPg copies has not been investigated in detail. Here, we generated a novel bioinformatics approach to identify VPg sequences in viral genomes using hidden Markov models (HMM) based on alignments of dicistrovirus VPg sequences. From metagenomic datasets of dicistrovirus genomes, we identified 717 dicistrovirus genomes containing VPgs ranging from a single copy to 8 tandem copies. The VPgs are classified into nine distinct types based on their sequence and length. The VPg types but not VPg numbers per viral genome followed specific virus clades, thus suggesting VPgs co-evolved with viral genomes. We also identified VPg duplications in aquamavirus and mosavirus genomes. This study greatly expands the number of viral genomes that contain VPg copies and indicates that duplicated viral sequences are more widespread than anticipated.

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

Droplet Microfluidics for High-Resolution Virology

Microfluidics has enabled a new era of cellular and molecular assays due to the small length scales, parallelization, and the modularity of various analysis and actuation functions. Droplet microfluidics, in particular, has been instrumental in providing new tools for biology with its ability to quickly and reproducibly generate drops that act as individual reactors. A notable beneficiary of this technology has been single-cell RNA sequencing, which has revealed new heterogeneities and interactions for the fundamental unit of life. However, viruses far surpass the diversity of cellular life, affect the dynamics of all ecosystems, and are a chronic source of global health crises. Despite their impact on the world, high-throughput and high-resolution viral profiling has been difficult, with conventional methods being limited to population-level averaging, large sample volumes, and few cultivable hosts. Consequently, most viruses have not been identified and studied. Droplet microfluidics holds the potential to address many of these limitations and offers new levels of sensitivity and throughput for virology. This Feature highlights recent efforts that have applied droplet microfluidics to the detection and study of viruses, including for diagnostics, virus-host interactions, and cell-independent virus assays. In combination with traditional virology methods, droplet microfluidics should prove a potent tool toward achieving a better understanding of the most abundant biological species on Earth.

A Novel Lineage of Cile-Like Viruses Discloses the Phylogenetic Continuum Across the Family Kitaviridae

An increasing number of plant species have been recognized or considered likely reservoirs of viruses transmitted by Brevipalpus mites. A tiny fraction of these viruses, primarily those causing severe economic burden to prominent crops, have been fully characterized. In this study, based on high-throughput sequencing, transmission electron microscopy analyses of virions in plant-infected tissues, viral transmission experiments, and the morphoanatomical identification of the involved Brevipalpus mites, we describe molecular and biological features of viruses representing three new tentative species of the family Kitaviridae. The genomes of Solanum violifolium ringspot virus (SvRSV, previously partially characterized), Ligustrum chlorotic spot virus (LigCSV), and Ligustrum leprosis virus (LigLV) have five open reading frames (ORFs) > 500 nts, two distributed in RNA1 and three in RNA2. RNA1 of these three viruses display the same genomic organization found in RNA1 of typical cileviruses, while their RNA2 are shorter, possessing only orthologs of genes p61, p32, and p24. LigCSV and LigLV are more closely related to each other than to SvRSV, but the identities between their genomic RNAs were lower than 70%. In gene-by-gene comparisons, ORFs from LigCSV and LigLV had the highest sequence identity values (nt sequences: 70–76% and deduced amino acid sequences: 74–83%). The next higher identity values were with ORFs from typical cileviruses, with values below 66%. Virions of LigLV (≈ 40 nm × 55 nm) and LigCSV (≈ 54 nm × 66 nm) appear almost spherical, contrasting with the bacilliform shape of SvRSV virions (≈ 47 nm × 101 nm). Mites collected from the virus-infected plants were identified as Brevipalpus papayensis, B. tucuman, and B. obovatus. Viruliferous B. papayensis mites successfully transmitted LigCSV to Arabidopsis thaliana. SvRSV, LigCSV, and LigLV seem to represent novel sub-lineages of kitaviruses that descent on parallel evolutionary branches from a common ancestor shared with the tentative cile-like virus hibiscus yellow blotch virus and typical cileviruses. Biological and molecular data, notably, the phylogenetic reconstruction based on the RdRp proteins in which strong support for monophyly of the family Kitaviridae is observed, mark an advance in the understanding of kitavirids.

Using networks to analyze and visualize the distribution of overlapping genes in virus genomes

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may increase the information content of compact genomes or influence the creation of new genes. Here we report a global comparative study of overlapping open reading frames (OvRFs) of 12,609 virus reference genomes in the NCBI database. We retrieved metadata associated with all annotated open reading frames (ORFs) in each genome record to calculate the number, length, and frameshift of OvRFs. Our results show that while the number of OvRFs increases with genome length, they tend to be shorter in longer genomes. The majority of overlaps involve +2 frameshifts, predominantly found in dsDNA viruses. Antisense overlaps in which one of the ORFs was encoded in the same frame on the opposite strand (−0) tend to be longer. Next, we develop a new graph-based representation of the distribution of overlaps among the ORFs of genomes in a given virus family. In the absence of an unambiguous partition of ORFs by homology at this taxonomic level, we used an alignment-free k-mer based approach to cluster protein coding sequences by similarity. We connect these clusters with two types of directed edges to indicate (1) that constituent ORFs are adjacent in one or more genomes, and (2) that these ORFs overlap. These adjacency graphs not only provide a natural visualization scheme, but also a novel statistical framework for analyzing the effects of gene- and genome-level attributes on the frequencies of overlaps.

Genetic Engineering Systems to Study Human Viral Pathogens from the Coronaviridae Family

The COVID-19 pandemic caused by the previously unknown SARS-CoV-2 Betacoronavirus made it extremely important to develop simple and safe cellular systems which allow manipulation of the viral genome and high-throughput screening of its potential inhibitors. In this review, we made an attempt at summarizing the currently existing data on genetic engineering systems used to study not only SARS-CoV-2, but also other viruses from the Coronaviridae family. In addition, the review covers the basic knowledge about the structure and the life cycle of coronaviruses.

Gene Overlapping as a Modulator of Begomovirus Evolution

In RNA viruses, which have high mutation—and fast evolutionary— rates, gene overlapping (i.e., genomic regions that encode more than one protein) is a major factor controlling mutational load and therefore the virus evolvability. Although DNA viruses use host high-fidelity polymerases for their replication, and therefore should have lower mutation rates, it has been shown that some of them have evolutionary rates comparable to those of RNA viruses. Notably, these viruses have large proportions of their genes with at least one overlapping instance. Hence, gene overlapping could be a modulator of virus evolution beyond the RNA world. To test this hypothesis, we use the genus Begomovirus of plant viruses as a model. Through comparative genomic approaches, we show that terminal gene overlapping decreases the rate of virus evolution, which is associated with lower frequency of both synonymous and nonsynonymous mutations. In contrast, terminal overlapping has little effect on the pace of virus evolution. Overall, our analyses support a role for gene overlapping in the evolution of begomoviruses and provide novel information on the factors that shape their genetic diversity.

Extending the Coding Potential of Viral Genomes with Overlapping Antisense ORFs: A Case for the De Novo Creation of the Gene Encoding the Antisense Protein ASP of HIV-1

Gene overprinting occurs when point mutations within a genomic region with an existing coding sequence create a new one in another reading frame. This process is quite frequent in viral genomes either to maximize the amount of information that they encode or in response to strong selective pressure. The most frequent scenario involves two different reading frames in the same DNA strand (sense overlap). Much less frequent are cases of overlapping genes that are encoded on opposite DNA strands (antisense overlap). One such example is the antisense ORF, asp in the minus strand of the HIV-1 genome overlapping the env gene. The asp gene is highly conserved in pandemic HIV-1 strains of group M, and it is absent in non-pandemic HIV-1 groups, HIV-2, and lentiviruses infecting non-human primates, suggesting that the ~190-amino acid protein that is expressed from this gene (ASP) may play a role in virus spread. While the function of ASP in the virus life cycle remains to be elucidated, mounting evidence from several research groups indicates that ASP is expressed in vivo. There are two alternative hypotheses that could be envisioned to explain the origin of the asp ORF. On one hand, asp may have originally been present in the ancestor of contemporary lentiviruses, and subsequently lost in all descendants except for most HIV-1 strains of group M due to selective advantage. Alternatively, the asp ORF may have originated very recently with the emergence of group M HIV-1 strains from SIVcpz. Here, we used a combination of computational and statistical approaches to study the genomic region of env in primate lentiviruses to shed light on the origin, structure, and sequence evolution of the asp ORF. The results emerging from our studies support the hypothesis of a recent de novo addition of the antisense ORF to the HIV-1 genome through a process that entailed progressive removal of existing internal stop codons from SIV strains to HIV-1 strains of group M, and fine tuning of the codon sequence in env that reduced the chances of new stop codons occurring in asp. Altogether, the study supports the notion that the HIV-1 asp gene encodes an accessory protein, providing a selective advantage to the virus.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

Biological factors in the synthetic construction of overlapping genes

BACKGROUND

Overlapping genes (OLGs) with long protein-coding overlapping sequences are disallowed by standard genome annotation programs, outside of viruses. Recently however they have been discovered in Archaea, diverse Bacteria, and Mammals. The biological factors underlying life's ability to create overlapping genes require more study, and may have important applications in understanding evolution and in biotechnology. A previous study claimed that protein domains from viruses were much better suited to forming overlaps than those from other cellular organisms - in this study we assessed this claim, in order to discover what might underlie taxonomic differences in the creation of gene overlaps.

RESULTS

After overlapping arbitrary Pfam domain pairs and evaluating them with Hidden Markov Models we find OLG construction to be much less constrained than expected. For instance, close to 10% of the constructed sequences cannot be distinguished from typical sequences in their protein family. Most are also indistinguishable from natural protein sequences regarding identity and secondary structure. Surprisingly, contrary to a previous study, virus domains were much less suitable for designing OLGs than bacterial or eukaryotic domains were. In general, the amount of amino acid change required to force a domain to overlap is approximately equal to the variation observed within a typical domain family. The resulting high similarity between natural sequences and those altered so as to overlap is mostly due to the combination of high redundancy in the genetic code and the evolutionary exchangeability of many amino acids.

CONCLUSIONS

Synthetic overlapping genes which closely resemble natural gene sequences, as measured by HMM profiles, are remarkably easy to construct, and most arbitrary domain pairs can be altered so as to overlap while retaining high similarity to the original sequences. Future work however will need to assess important factors not considered such as intragenic interactions which affect protein folding. While the analysis here is not sufficient to guarantee functional folding proteins, further analysis of constructed OLGs will improve our understanding of the origin of these remarkable genetic elements across life and opens up exciting possibilities for synthetic biology.

Top Three Strategies of ss(+)RNA Plant Viruses: Great Opportunists and Ecosystem Tuners with a Small Genome

ss(+)RNA viruses represent the dominant group of plant viruses. They owe their evolutionary superiority to the large number of mutations that occur during replication, courtesy of RNA-dependent RNA polymerase. Natural selection rewards successful viral subtypes, whose effective tuning of the ecosystem regulates the interactions between its participants. Thus, ss(+)RNA viruses act as shuttles for the functionally important genes of the participants in symbiotic relationships within the ecosystem, of which the most common ecological triad is “plant–virus–insect”. Due to their short life cycle and large number of offspring, RNA viruses act as skillful tuners of the ecosystem, which benefits both viruses and the system as a whole. A fundamental understanding of this aspect of the role played by viruses in the ecosystem makes it possible to apply this knowledge to the creation of DNA insecticides. In fact, since the genes that viruses are involved in transferring are functionally important for both insects and plants, silencing these genes (for example, in insects) can be used to regulate the pest population. RNA viruses are increasingly treated not as micropathogens but as necessary regulators of ecosystem balance.

Computational methods for inferring location and genealogy of overlapping genes in virus genomes: approaches and applications

Viruses may evolve to increase the amount of encoded genetic information by means of overlapping genes, which utilize several reading frames. Such overlapping genes may be especially impactful for genomes of small size, often serving a source of novel accessory proteins, some of which play a crucial role in viral pathogenicity or in promoting the systemic spread of virus. Diverse genome-based metrics were proposed to facilitate recognition of overlapping genes that otherwise may be overlooked during genome annotation. They can detect the atypical codon bias associated with the overlap (e.g. a statistically significant reduction in variability at synonymous sites) or other sequence-composition features peculiar to overlapping genes. In this review, I compare nine computational methods, discuss their strengths and limitations, and survey how they were applied to detect candidate overlapping genes in the genome of SARS-CoV-2, the etiological agent of COVID-19 pandemic.

Poorly Conserved P15 Proteins of Cileviruses Retain Elements of Common Ancestry and Putative Functionality: A Theoretical Assessment on the Evolution of Cilevirus Genomes

The genus Cilevirus groups enveloped single-stranded (+) RNA virus members of the family Kitaviridae, order Martellivirales. Proteins P15, scarcely conserved polypeptides encoded by cileviruses, have no apparent homologs in public databases. Accordingly, the open reading frames (ORFs) p15, located at the 5'-end of the viral RNA2 molecules, are considered orphan genes (ORFans). In this study, we have delved into ORFs p15 and the relatively poorly understood biochemical properties of the proteins P15 to posit their importance for viruses across the genus and theorize on their origin. We detected that the ORFs p15 are under purifying selection and that, in some viral strains, the use of synonymous codons is biased, which might be a sign of adaptation to their plant hosts. Despite the high amino acid sequence divergence, proteins P15 show the conserved motif [FY]-L-x(3)-[FL]-H-x-x-[LIV]-S-C-x-C-x(2)-C-x-G-x-C, which occurs exclusively in members of this protein family. Proteins P15 also show a common predicted 3D structure that resembles the helical scaffold of the protein ORF49 encoded by radinoviruses and the phosphoprotein C-terminal domain of mononegavirids. Based on the 3D structural similarities of P15, we suggest elements of common ancestry, conserved functionality, and relevant amino acid residues. We conclude by postulating a plausible evolutionary trajectory of ORFans p15 and the 5'-end of the RNA2 of cileviruses considering both protein fold superpositions and comparative genomic analyses with the closest kitaviruses, negeviruses, nege/kita-like viruses, and unrelated viruses that share the ecological niches of cileviruses.

Towards an integrative view of virus phenotypes

Understanding how phenotypes emerge from genotypes is a foundational goal in biology. As challenging as this task is when considering cellular life, it is further complicated in the case of viruses. During replication, a virus as a discrete entity (the virion) disappears and manifests itself as a metabolic amalgam between the virus and the host (the virocell). Identifying traits that unambiguously constitute a virus's phenotype is straightforward for the virion, less so for the virocell. Here, we present a framework for categorizing virus phenotypes that encompasses both virion and virocell stages and considers functional and performance traits of viruses in the context of fitness. Such an integrated view of virus phenotype is necessary for comprehensive interpretation of viral genome sequences and will advance our understanding of viral evolution and ecology.

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

Origin, Evolution and Stability of Overlapping Genes in Viruses: A Systematic Review

During their long evolutionary history viruses generated many proteins de novo by a mechanism called “overprinting”. Overprinting is a process in which critical nucleotide substitutions in a pre-existing gene can induce the expression of a novel protein by translation of an alternative open reading frame (ORF). Overlapping genes represent an intriguing example of adaptive conflict, because they simultaneously encode two proteins whose freedom to change is constrained by each other. However, overlapping genes are also a source of genetic novelties, as the constraints under which alternative ORFs evolve can give rise to proteins with unusual sequence properties, most importantly the potential for novel functions. Starting with the discovery of overlapping genes in phages infecting Escherichia coli, this review covers a range of studies dealing with detection of overlapping genes in small eukaryotic viruses (genomic length below 30 kb) and recognition of their critical role in the evolution of pathogenicity. Origin of overlapping genes, what factors favor their birth and retention, and how they manage their inherent adaptive conflict are extensively reviewed. Special attention is paid to the assembly of overlapping genes into ad hoc databases, suitable for future studies, and to the development of statistical methods for exploring viral genome sequences in search of undiscovered overlaps.

The HIV-1 Antisense Gene ASP: The New Kid on the Block

Viruses have developed incredibly creative ways of making a virtue out of necessity, including taking full advantage of their small genomes. Indeed, viruses often encode multiple proteins within the same genomic region by using two or more reading frames in both orientations through a process called overprinting. Complex retroviruses provide compelling examples of that. The human immunodeficiency virus type 1 (HIV-1) genome expresses sixteen proteins from nine genes that are encoded in the three positive-sense reading frames. In addition, the genome of some HIV-1 strains contains a tenth gene in one of the negative-sense reading frames. The so-called Antisense Protein (ASP) gene overlaps the HIV-1 Rev Response Element (RRE) and the envelope glycoprotein gene, and encodes a highly hydrophobic protein of ~190 amino acids. Despite being identified over thirty years ago, relatively few studies have investigated the role that ASP may play in the virus lifecycle, and its expression in vivo is still questioned. Here we review the current knowledge about ASP, and we discuss some of the many unanswered questions.

OLGenie: Estimating Natural Selection to Predict Functional Overlapping Genes

Purifying (negative) natural selection is a hallmark of functional biological sequences, and can be detected in protein-coding genes using the ratio of nonsynonymous to synonymous substitutions per site (d_N/d_S). However, when two genes overlap the same nucleotide sites in different frames, synonymous changes in one gene may be nonsynonymous in the other, perturbing d_N/d_S. Thus, scalable methods are needed to estimate functional constraint specifically for overlapping genes (OLGs). We propose OLGenie, which implements a modification of the Wei–Zhang method. Assessment with simulations and controls from viral genomes (58 OLGs and 176 non-OLGs) demonstrates low false-positive rates and good discriminatory ability in differentiating true OLGs from non-OLGs. We also apply OLGenie to the unresolved case of HIV-1’s putative antisense protein gene, showing significant purifying selection. OLGenie can be used to study known OLGs and to predict new OLGs in genome annotation. Software and example data are freely available at https://github.com/chasewnelson/OLGenie (last accessed April 10, 2020).

The Picornavirus Precursor 3CD Has Different Conformational Dynamics Compared to 3Cpro and 3Dpol in Functionally Relevant Regions

Viruses have evolved numerous strategies to maximize the use of their limited genetic material, including proteolytic cleavage of polyproteins to yield products with different functions. The poliovirus polyprotein 3CD is involved in important protein-protein, protein-RNA and protein-lipid interactions in viral replication and infection. It is a precursor to the 3C protease and 3D RNA-dependent RNA polymerase, but has different protease specificity, is not an active polymerase, and participates in other interactions differently than its processed products. These functional differences are poorly explained by the known X-ray crystal structures. It has been proposed that functional differences might be due to differences in conformational dynamics between 3C, 3D and 3CD. To address this possibility, we conducted nuclear magnetic resonance spectroscopy experiments, including multiple quantum relaxation dispersion, chemical exchange saturation transfer and methyl spin-spin relaxation, to probe conformational dynamics across multiple timescales. Indeed, these studies identified differences in conformational dynamics in functionally important regions, including enzyme active sites, and RNA and lipid binding sites. Expansion of the conformational ensemble available to 3CD may allow it to perform additional functions not observed in 3C and 3D alone despite having nearly identical lowest-energy structures.

Detection of SARS-CoV-2 from patient fecal samples by whole genome sequencing

BACKGROUND

SARS-CoV-2 has been detected not only in respiratory secretions, but also in stool collections. Here were sought to identify SARS-CoV-2 by enrichment next-generation sequencing (NGS) from fecal samples, and to utilize whole genome analysis to characterize SARS-CoV-2 mutational variations in COVID-19 patients.

RESULTS

Study participants underwent testing for SARS-CoV-2 from fecal samples by whole genome enrichment NGS (n = 14), and RT-PCR nasopharyngeal swab analysis (n = 12). The concordance of SARS-CoV-2 detection by enrichment NGS from stools with RT-PCR nasopharyngeal analysis was 100%. Unique variants were identified in four patients, with a total of 33 different mutations among those in which SARS-CoV-2 was detected by whole genome enrichment NGS.

CONCLUSION

These results highlight the potential viability of SARS-CoV-2 in feces, its ongoing mutational accumulation, and its possible role in fecal-oral transmission. This study also elucidates the advantages of SARS-CoV-2 enrichment NGS, which may be a key methodology to document complete viral eradication. Trial registration ClinicalTrials.gov, NCT04359836, Registered 24 April 2020, https://clinicaltrials.gov/ct2/show/NCT04359836?term=NCT04359836&draw=2&rank=1 ).

Tuning intrinsic disorder predictors for virus proteins

Many virus-encoded proteins have intrinsically disordered regions that lack a stable, folded three-dimensional structure. These disordered proteins often play important functional roles in virus replication, such as down-regulating host defense mechanisms. With the widespread availability of next-generation sequencing, the number of new virus genomes with predicted open reading frames is rapidly outpacing our capacity for directly characterizing protein structures through crystallography. Hence, computational methods for structural prediction play an important role. A large number of predictors focus on the problem of classifying residues into ordered and disordered regions, and these methods tend to be validated on a diverse training set of proteins from eukaryotes, prokaryotes, and viruses. In this study, we investigate whether some predictors outperform others in the context of virus proteins and compared our findings with data from non-viral proteins. We evaluate the prediction accuracy of 21 methods, many of which are only available as web applications, on a curated set of 126 proteins encoded by viruses. Furthermore, we apply a random forest classifier to these predictor outputs. Based on cross-validation experiments, this ensemble approach confers a substantial improvement in accuracy, e.g., a mean 36 per cent gain in Matthews correlation coefficient. Lastly, we apply the random forest predictor to severe acute respiratory syndrome coronavirus 2 ORF6, an accessory gene that encodes a short (61 AA) and moderately disordered protein that inhibits the host innate immune response. We show that disorder prediction methods perform differently for viral and non-viral proteins, and that an ensemble approach can yield more robust and accurate predictions.

Evolutionary history of cotranscriptional editing in the paramyxoviral phosphoprotein gene

The phosphoprotein gene of the paramyxoviruses encodes multiple protein products. The P, V, and W proteins are generated by transcriptional slippage. This process results in the insertion of non-templated guanosine nucleosides into the mRNA at a conserved edit site. The P protein is an essential component of the viral RNA polymerase and is encoded by a faithful copy of the gene in the majority of paramyxoviruses. However, in some cases, the non-essential V protein is encoded by default and guanosines must be inserted into the mRNA in order to encode P. The number of guanosines inserted into the P gene can be described by a probability distribution, which varies between viruses. In this article, we review the nature of these distributions, which can be inferred from mRNA sequencing data, and reconstruct the evolutionary history of cotranscriptional editing in the paramyxovirus family. Our model suggests that, throughout known history of the family, the system has switched from a P default to a V default mode four times; complete loss of the editing system has occurred twice, the canonical zinc finger domain of the V protein has been deleted or heavily mutated a further two times, and the W protein has independently evolved a novel function three times. Finally, we review the physical mechanisms of cotranscriptional editing via slippage of the viral RNA polymerase.

Genome Modularization Reveals Overlapped Gene Topology Is Necessary for Efficient Viral Reproduction

Sequence overlap between two genes is common across all genomes, with viruses having high proportions of these gene overlaps. Genome modularization and refactoring is the process of disrupting natural gene overlaps to separate coding sequences to enable their individual manipulation. The biological function and fitness effects of gene overlaps are not fully understood, and their effects on gene cluster and genome-level refactoring are unknown. The bacteriophage φX174 genome has ∼26% of nucleotides involved in encoding more than one gene. In this study we use an engineered φX174 phage containing a genome with all gene overlaps removed to show that gene overlap is critical to maintaining optimal viral fecundity. Through detailed phenotypic measurements we reveal that genome modularization in φX174 causes virion replication, stability, and attachment deficiencies. Quantitation of the complete phage proteome across an infection cycle reveals 30% of proteins display abnormal expression patterns. Taken together, we have for the first time comprehensively demonstrated that gene modularization severely perturbs the coordinated functioning of a bacteriophage replication cycle. This work highlights the biological importance of gene overlap in natural genomes and that reducing gene overlap disruption should be an integral part of future genome engineering projects.

Adaptation of turnip mosaic potyvirus to a specific niche reduces its genetic and environmental robustness

Robustness is the preservation of the phenotype in the face of genetic and environmental perturbations. It has been argued that robustness must be an essential fitness component of RNA viruses owed to their small and compacted genomes, high mutation rates and living in ever-changing environmental conditions. Given that genetic robustness might hamper possible beneficial mutations, it has been suggested that genetic robustness can only evolve as a side-effect of the evolution of robustness mechanisms specific to cope with environmental perturbations, a theory known as plastogenetic congruence. However, empirical evidences from different viral systems are contradictory. To test how adaptation to a particular environment affects both environmental and genetic robustness, we have used two strains of turnip mosaic potyvirus (TuMV) that differ in their degree of adaptation to Arabidopsis thaliana at a permissive temperature. We show that the highly adapted strain is strongly sensitive to the effect of random mutations and to changes in temperature conditions. In contrast, the non-adapted strain shows more robustness against both the accumulation of random mutations and drastic changes in temperature conditions. Together, these results are consistent with the predictions of the plastogenetic congruence theory, suggesting that genetic and environmental robustnesses may be two sides of the same coin for TuMV.

New insights into the evolutionary features of viral overlapping genes by discriminant analysis

Overlapping genes originate by a mechanism of overprinting, in which nucleotide substitutions in a pre-existing frame induce the expression of a de novo protein from an alternative frame. In this study, I assembled a dataset of 319 viral overlapping genes, which included 82 overlaps whose expression is experimentally known and the respective 237 homologs. Principal component analysis revealed that overlapping genes have a common pattern of nucleotide and amino acid composition. Discriminant analysis separated overlapping from non-overlapping genes with an accuracy of 97%. When applied to overlapping genes with known genealogy, it separated ancestral from de novo frames with an accuracy close to 100%. This high discriminant power was crucial to computationally design variants of de novo viral proteins known to possess selective anticancer toxicity (apoptin) or protection against neurodegeneration (X protein), as well as to detect two new potential overlapping genes in the genome of the new coronavirus SARS-CoV-2.

Properties and abundance of overlapping genes in viruses

Overlapping genes are commonplace in viruses and play an important role in their function and evolution. However, aside from studies on specific groups of viruses, relatively little is known about the extent and nature of gene overlap and its determinants in viruses as a whole. Here, we present an extensive characterisation of gene overlap in viruses through an analysis of reference genomes present in the NCBI virus genome database. We find that over half the instances of gene overlap are very small, covering <10 nt, and 84 per cent are <50 nt in length. Despite this, 53 per cent of all viruses still contained a gene overlap of 50 nt or larger. We also investigate several predictors of gene overlap such as genome structure (single- and double-stranded RNA and DNA), virus family, genome length, and genome segmentation. This revealed that gene overlap occurs more frequently in DNA viruses than in RNA viruses, and more frequently in single-stranded viruses than in double-stranded viruses. Genome segmentation is also associated with gene overlap, particularly in single-stranded DNA viruses. Notably, we observed a large range of overlap frequencies across families of all genome types, suggesting that it is a common evolutionary trait that provides flexible genome structures in all virus families.

A case for a negative-strand coding sequence in a group of positive-sense RNA viruses

Positive-sense single-stranded RNA viruses form the largest and most diverse group of eukaryote-infecting viruses. Their genomes comprise one or more segments of coding-sense RNA that function directly as messenger RNAs upon release into the cytoplasm of infected cells. Positive-sense RNA viruses are generally accepted to encode proteins solely on the positive strand. However, we previously identified a surprisingly long (∼1,000-codon) open reading frame (ORF) on the negative strand of some members of the family Narnaviridae which, together with RNA bacteriophages of the family Leviviridae, form a sister group to all other positive-sense RNA viruses. Here, we completed the genomes of three mosquito-associated narnaviruses, all of which have the long reverse-frame ORF. We systematically identified narnaviral sequences in public data sets from a wide range of sources, including arthropod, fungal, and plant transcriptomic data sets. Long reverse-frame ORFs are widespread in one clade of narnaviruses, where they frequently occupy >95 per cent of the genome. The reverse-frame ORFs correspond to a specific avoidance of CUA, UUA, and UCA codons (i.e. stop codon reverse complements) in the forward-frame RNA-dependent RNA polymerase ORF. However, absence of these codons cannot be explained by other factors such as inability to decode these codons or GC3 bias. Together with other analyses, we provide the strongest evidence yet of coding capacity on the negative strand of a positive-sense RNA virus. As these ORFs comprise some of the longest known overlapping genes, their study may be of broad relevance to understanding overlapping gene evolution and de novo origin of genes.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes

Background

Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown.

Results

Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes.

Conclusions

This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

An exploration of ambigrammatic sequences in narnaviruses

Narnaviruses have been described as positive-sense RNA viruses with a remarkably simple genome of ~3 kb, encoding only a highly conserved RNA-dependent RNA polymerase (RdRp). Many narnaviruses, however, are 'ambigrammatic' and harbour an additional uninterrupted open reading frame (ORF) covering almost the entire length of the reverse complement strand. No function has been described for this ORF, yet the absence of stops is conserved across diverse narnaviruses, and in every case the codons in the reverse ORF and the RdRp are aligned. The >3 kb ORF overlap on opposite strands, unprecedented among RNA viruses, motivates an exploration of the constraints imposed or alleviated by the codon alignment. Here, we show that only when the codon frames are aligned can all stop codons be eliminated from the reverse strand by synonymous single-nucleotide substitutions in the RdRp gene, suggesting a mechanism for de novo gene creation within a strongly conserved amino-acid sequence. It will be fascinating to explore what implications this coding strategy has for other aspects of narnavirus biology. Beyond narnaviruses, our rapidly expanding catalogue of viral diversity may yet reveal additional examples of this broadly-extensible principle for ambigrammatic-sequence development.

Harnessing host-virus evolution in antiviral therapy and immunotherapy

Pathogen resistance and development costs are major challenges in current approaches to antiviral therapy. The high error rate of RNA synthesis and reverse‐transcription confers genome plasticity, enabling the remarkable adaptability of RNA viruses to antiviral intervention. However, this property is coupled to fundamental constraints including limits on the size of information available to manipulate complex hosts into supporting viral replication. Accordingly, RNA viruses employ various means to extract maximum utility from their informationally limited genomes that, correspondingly, may be leveraged for effective host‐oriented therapies. Host‐oriented approaches are becoming increasingly feasible because of increased availability of bioactive compounds and recent advances in immunotherapy and precision medicine, particularly genome editing, targeted delivery methods and RNAi. In turn, one driving force behind these innovations is the increasingly detailed understanding of evolutionarily diverse host–virus interactions, which is the key concern of an emerging field, neo‐virology. This review examines biotechnological solutions to disease and other sustainability issues of our time that leverage the properties of RNA and DNA viruses as developed through co‐evolution with their hosts.

On the stability of sequences inserted into viral genomes

Viruses are widely used as vectors for heterologous gene expression in cultured cells or natural hosts, and therefore a large number of viruses with exogenous sequences inserted into their genomes have been engineered. Many of these engineered viruses are viable and express heterologous proteins at high levels, but the inserted sequences often prove to be unstable over time and are rapidly lost, limiting heterologous protein expression. Although virologists are aware that inserted sequences can be unstable, processes leading to insert instability are rarely considered from an evolutionary perspective. Here, we review experimental work on the stability of inserted sequences over a broad range of viruses, and we present some theoretical considerations concerning insert stability. Different virus genome organizations strongly impact insert stability, and factors such as the position of insertion can have a strong effect. In addition, we argue that insert stability not only depends on the characteristics of a particular genome, but that it will also depend on the host environment and the demography of a virus population. The interplay between all factors affecting stability is complex, which makes it challenging to develop a general model to predict the stability of genomic insertions. We highlight key questions and future directions, finding that insert stability is a surprisingly complex problem and that there is need for mechanism-based, predictive models. Combining theoretical models with experimental tests for stability under varying conditions can lead to improved engineering of viral modified genomes, which is a valuable tool for understanding genome evolution as well as for biotechnological applications, such as gene therapy.

A Simple Method to Detect Candidate Overlapping Genes in Viruses Using Single Genome Sequences

Overlapping genes in viruses maximize the coding capacity of their genomes and allow the generation of new genes without major increases in genome size. Despite their importance, the evolution and function of overlapping genes are often not well understood, in part due to difficulties in their detection. In addition, most bioinformatic approaches for the detection of overlapping genes require the comparison of multiple genome sequences that may not be available in metagenomic surveys of virus biodiversity. We introduce a simple new method for identifying candidate functional overlapping genes using single virus genome sequences. Our method uses randomization tests to estimate the expected length of open reading frames and then identifies overlapping open reading frames that significantly exceed this length and are thus predicted to be functional. We applied this method to 2548 reference RNA virus genomes and find that it has both high sensitivity and low false discovery for genes that overlap by at least 50 nucleotides. Notably, this analysis provided evidence for 29 previously undiscovered functional overlapping genes, some of which are coded in the antisense direction suggesting there are limitations in our current understanding of RNA virus replication.

Detection of antisense protein (ASP) RNA transcripts in individuals infected with human immunodeficiency virus type 1 (HIV-1)

The detection of antisense RNA is hampered by reverse transcription (RT) non-specific priming, due to the ability of RNA secondary structures to prime RT in the absence of specific primers. The detection of antisense RNA by conventional RT-PCR does not allow assessment of the polarity of the initial RNA template, causing the amplification of non-specific cDNAs. In this study we have developed a modified protocol for the detection of human immunodeficiency virus type 1 (HIV-1) antisense protein (ASP) RNA. Using this approach, we have identified ASP transcripts in CD4+ T cells isolated from five HIV-infected individuals, either untreated or under suppressive therapy. We show that ASP RNA can be detected in stimulated CD4+ T cells from both groups of patients, but not in unstimulated cells. We also show that in untreated patients, the patterns of expression of ASP and env are very similar, with the levels of ASP RNA being markedly lower than those of env. Treatment of cells from one viraemic patient with α-amanitin greatly reduces the rate of ASP RNA synthesis, suggesting that it is associated with RNA polymerase II, the central enzyme in the transcription of protein-coding genes. Our data represent the first nucleotide sequences obtained in patients for ASP, demonstrating that its transcription indeed occurs in those HIV-1 lineages in which the ASP open reading frame is present.

cRegions-a tool for detecting conserved cis-elements in multiple sequence alignment of diverged coding sequences

Identifying cis-acting elements and understanding regulatory mechanisms of a gene is crucial to fully understand the molecular biology of an organism. In general, it is difficult to identify previously uncharacterised cis-acting elements with an unknown consensus sequence. The task is especially problematic with viruses containing regions of limited or no similarity to other previously characterised sequences. Fortunately, the fast increase in the number of sequenced genomes allows us to detect some of these elusive cis-elements. In this work, we introduce a web-based tool called cRegions. It was developed to identify regions within a protein-coding sequence where the conservation in the amino acid sequence is caused by the conservation in the nucleotide sequence. The cRegion can be the first step in discovering novel cis-acting sequences from diverged protein-coding genes. The results can be used as a basis for future experimental analysis. We applied cRegions on the non-structural and structural polyproteins of alphaviruses as an example and successfully detected all known cis-acting elements. In this publication and in previous work, we have shown that cRegions is able to detect a wide variety of functional elements in DNA and RNA viruses. These functional elements include splice sites, stem-loops, overlapping reading frames, internal promoters, ribosome frameshifting signals and other embedded elements with yet unknown function. The cRegions web tool is available at http://bioinfo.ut.ee/cRegions/.

HIV-1 Reverse Transcriptase: A Metamorphic Protein with Three Stable States

There has been a steadily increasing appreciation of the fact that the relationship between protein sequence and structure is often sufficiently ambiguous to allow a single sequence to adopt alternative, stable folds. Living organisms have been able to utilize such metamorphic proteins in remarkable and unanticipated ways. HIV-1 reverse transcriptase is among the earliest such proteins identified and remains a unique example in which a functional heterodimer contains two, alternatively folded polymerase domains. Structural characterization of the p66 precursor protein combined with NMR spectroscopic and molecular modeling studies have provided insights into the factors underlying the metamorphic transition and the subunit-specific programmed unfolding step required to expose the protease cleavage site within the ribonuclease H domain, supporting the conversion of the p66/p66' precursor into the mature p66/p51 heterodimer.

Evolutionary Virology at 40

RNA viruses are diverse, abundant, and rapidly evolving. Genetic data have been generated from virus populations since the late 1970s and used to understand their evolution, emergence, and spread, culminating in the generation and analysis of many thousands of viral genome sequences. Despite this wealth of data, evolutionary genetics has played a surprisingly small role in our understanding of virus evolution. Instead, studies of RNA virus evolution have been dominated by two very different perspectives, the experimental and the comparative, that have largely been conducted independently and sometimes antagonistically. Here, we review the insights that these two approaches have provided over the last 40 years. We show that experimental approaches using in vitro and in vivo laboratory models are largely focused on short-term intrahost evolutionary mechanisms, and may not always be relevant to natural systems. In contrast, the comparative approach relies on the phylogenetic analysis of natural virus populations, usually considering data collected over multiple cycles of virus-host transmission, but is divorced from the causative evolutionary processes. To truly understand RNA virus evolution it is necessary to meld experimental and comparative approaches within a single evolutionary genetic framework, and to link viral evolution at the intrahost scale with that which occurs over both epidemiological and geological timescales. We suggest that the impetus for this new synthesis may come from methodological advances in next-generation sequencing and metagenomics.

A planarian nidovirus expands the limits of RNA genome size

RNA viruses are the only known RNA-protein (RNP) entities capable of autonomous replication (albeit within a permissive host environment). A 33.5 kilobase (kb) nidovirus has been considered close to the upper size limit for such entities; conversely, the minimal cellular DNA genome is in the 100–300 kb range. This large difference presents a daunting gap for the transition from primordial RNP to contemporary DNA-RNP-based life. Whether or not RNA viruses represent transitional steps towards DNA-based life, studies of larger RNA viruses advance our understanding of the size constraints on RNP entities and the role of genome size in virus adaptation. For example, emergence of the largest previously known RNA genomes (20–34 kb in positive-stranded nidoviruses, including coronaviruses) is associated with the acquisition of a proofreading exoribonuclease (ExoN) encoded in the open reading frame 1b (ORF1b) in a monophyletic subset of nidoviruses. However, apparent constraints on the size of ORF1b, which encodes this and other key replicative enzymes, have been hypothesized to limit further expansion of these viral RNA genomes. Here, we characterize a novel nidovirus (planarian secretory cell nidovirus; PSCNV) whose disproportionately large ORF1b-like region including unannotated domains, and overall 41.1-kb genome, substantially extend the presumed limits on RNA genome size. This genome encodes a predicted 13,556-aa polyprotein in an unconventional single ORF, yet retains canonical nidoviral genome organization and expression, as well as key replicative domains. These domains may include functionally relevant substitutions rarely or never before observed in highly conserved sites of RdRp, NiRAN, ExoN and 3CLpro. Our evolutionary analysis suggests that PSCNV diverged early from multi-ORF nidoviruses, and acquired additional genes, including those typical of large DNA viruses or hosts, e.g. Ankyrin and Fibronectin type II, which might modulate virus-host interactions. PSCNV's greatly expanded genome, proteomic complexity, and unique features–impressive in themselves–attest to the likelihood of still-larger RNA genomes awaiting discovery.

RNA viruses are the only known RNA-protein (RNP) entities capable of autonomous replication. The upper genome size for such entities was assumed to be <35 kb; conversely, the minimal cellular DNA genome is in the 100–300 kilobase (kb) range. This large difference presents a daunting gap for the proposed evolution of contemporary DNA-RNP-based life from primordial RNP entities. Here, we describe a nidovirus from planarians, named planarian secretory cell nidovirus (PSCNV), whose 41.1 kb genome is 23% larger than any riboviral genome yet discovered. This increase is nearly equivalent in size to the entire poliovirus genome, and it equips PSCNV with an unprecedented extra coding capacity to adapt. PSCNV has broken apparent constraints on the size of the genomic subregion that encodes core replication machinery in other nidoviruses, including coronaviruses, and has acquired genes not previously observed in RNA viruses. This virus challenges and advances our understanding of the limits to RNA genome size.

Overlapping genes and the proteins they encode differ significantly in their sequence composition from non-overlapping genes

Overlapping genes represent a fascinating evolutionary puzzle, since they encode two functionally unrelated proteins from the same DNA sequence. They originate by a mechanism of overprinting, in which point mutations in an existing frame allow the expression (the "birth") of a completely new protein from a second frame. In viruses, in which overlapping genes are abundant, these new proteins often play a critical role in infection, yet they are frequently overlooked during genome annotation. This results in erroneous interpretation of mutational studies and in a significant waste of resources. Therefore, overlapping genes need to be correctly detected, especially since they are now thought to be abundant also in eukaryotes. Developing better detection methods and conducting systematic evolutionary studies require a large, reliable benchmark dataset of known cases. We thus assembled a high-quality dataset of 80 viral overlapping genes whose expression is experimentally proven. Many of them were not present in databases. We found that overall, overlapping genes differ significantly from non-overlapping genes in their nucleotide and amino acid composition. In particular, the proteins they encode are enriched in high-degeneracy amino acids and depleted in low-degeneracy ones, which may alleviate the evolutionary constraints acting on overlapping genes. Principal component analysis revealed that the vast majority of overlapping genes follow a similar composition bias, despite their heterogeneity in length and function. Six proven mammalian overlapping genes also followed this bias. We propose that this apparently near-universal composition bias may either favour the birth of overlapping genes, or/and result from selection pressure acting on them.