Conserved Secondary Structures in Viral mRNAs

Ligand cross-feeding resolves bacterial vitamin B12 auxotrophies

Cobalamin (vitamin B12, herein referred to as B12) is an essential cofactor for most marine prokaryotes and eukaryotes1,2. Synthesized by a limited number of prokaryotes, its scarcity affects microbial interactions and community dynamics2-4. Here we show that two bacterial B12 auxotrophs can salvage different B12 building blocks and cooperate to synthesize B12. A Colwellia sp. synthesizes and releases the activated lower ligand α-ribazole, which is used by another B12 auxotroph, a Roseovarius sp., to produce the corrin ring and synthesize B12. Release of B12 by Roseovarius sp. happens only in co-culture with Colwellia sp. and only coincidently with the induction of a prophage encoded in Roseovarius sp. Subsequent growth of Colwellia sp. in these conditions may be due to the provision of B12 by lysed cells of Roseovarius sp. Further evidence is required to support a causative role for prophage induction in the release of B12. These complex microbial interactions of ligand cross-feeding and joint B12 biosynthesis seem to be widespread in marine pelagic ecosystems. In the western and northern tropical Atlantic Ocean, bacteria predicted to be capable of salvaging cobinamide and synthesizing only the activated lower ligand outnumber B12 producers. These findings add new players to our understanding of B12 supply to auxotrophic microorganisms in the ocean and possibly in other ecosystems.

VirClust-A Tool for Hierarchical Clustering, Core Protein Detection and Annotation of (Prokaryotic) Viruses

Recent years have seen major changes in the classification criteria and taxonomy of viruses. The current classification scheme, also called "megataxonomy of viruses", recognizes six different viral realms, defined based on the presence of viral hallmark genes (VHGs). Within the realms, viruses are classified into hierarchical taxons, ideally defined by the phylogeny of their shared genes. To enable the detection of shared genes, viruses have first to be clustered, and there is currently a need for tools to assist with virus clustering and classification. Here, VirClust is presented. It is a novel, reference-free tool capable of performing: (i) protein clustering, based on BLASTp and Hidden Markov Models (HMMs) similarities; (ii) hierarchical clustering of viruses based on intergenomic distances calculated from their shared protein content; (iii) identification of core proteins and (iv) annotation of viral proteins. VirClust has flexible parameters both for protein clustering and for splitting the viral genome tree into smaller genome clusters, corresponding to different taxonomic levels. Benchmarking on a phage dataset showed that the genome trees produced by VirClust match the current ICTV classification at family, sub-family and genus levels. VirClust is freely available, as a web-service and stand-alone tool.

IFN-Induced PARPs—Sensors of Foreign Nucleic Acids?

Cells have developed different strategies to cope with viral infections. Key to initiating a defense response against viruses is the ability to distinguish foreign molecules from their own. One central mechanism is the perception of foreign nucleic acids by host proteins which, in turn, initiate an efficient immune response. Nucleic acid sensing pattern recognition receptors have evolved, each targeting specific features to discriminate viral from host RNA. These are complemented by several RNA-binding proteins that assist in sensing of foreign RNAs. There is increasing evidence that the interferon-inducible ADP-ribosyltransferases (ARTs; PARP9—PARP15) contribute to immune defense and attenuation of viruses. However, their activation, subsequent targets, and precise mechanisms of interference with viruses and their propagation are still largely unknown. Best known for its antiviral activities and its role as RNA sensor is PARP13. In addition, PARP9 has been recently described as sensor for viral RNA. Here we will discuss recent findings suggesting that some PARPs function in antiviral innate immunity. We expand on these findings and integrate this information into a concept that outlines how the different PARPs might function as sensors of foreign RNA. We speculate about possible consequences of RNA binding with regard to the catalytic activities of PARPs, substrate specificity and signaling, which together result in antiviral activities.

Structural and Functional RNA Motifs of SARS-CoV-2 and Influenza A Virus as a Target of Viral Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic, whereas the influenza A virus (IAV) causes seasonal epidemics and occasional pandemics. Both viruses lead to widespread infection and death. SARS-CoV-2 and the influenza virus are RNA viruses. The SARS-CoV-2 genome is an approximately 30 kb, positive sense, 5' capped single-stranded RNA molecule. The influenza A virus genome possesses eight single-stranded negative-sense segments. The RNA secondary structure in the untranslated and coding regions is crucial in the viral replication cycle. The secondary structure within the RNA of SARS-CoV-2 and the influenza virus has been intensively studied. Because the whole of the SARS-CoV-2 and influenza virus replication cycles are dependent on RNA with no DNA intermediate, the RNA is a natural and promising target for the development of inhibitors. There are a lot of RNA-targeting strategies for regulating pathogenic RNA, such as small interfering RNA for RNA interference, antisense oligonucleotides, catalytic nucleic acids, and small molecules. In this review, we summarized the knowledge about the inhibition of SARS-CoV-2 and influenza A virus propagation by targeting their RNA secondary structure.

RNA G-quadruplex forming regions from SARS-2, SARS-1 and MERS coronoviruses

COVID-19 (Corona Virus Disease 2019), SARS (Severe Acute Respiratory Syndrome) and MERS (Middle East Respiratory Syndrome) are infectious diseases each caused by coronavirus outbreaks. Small molecules and other therapeutics are rapidly being developed to treat these diseases, but the threat of new variants and outbreaks argue for the identification of additional viral targets. Here we identify regions in each of the three coronavirus genomes that are able to form G-quadruplex (G4) structures. G4s are structures formed by DNA or RNA with a core of two or more stacked planes of guanosine tetrads. In recent years, numerous DNA and RNA G4s have emerged as promising pharmacological targets for the treatment of cancer and viral infection. We use a combination of bioinformatics and biophysical approaches to identify conserved RNA G4 regions from the ORF1A and S sequences of SARS-CoV, SARS-CoV-2 and MERS-CoV. Although a general depletion of G4-forming regions is observed in coronaviridae, the preservation of these selected G4 sequences support a significance in viral replication. Targeting these RNA structures may represent a new antiviral strategy against these viruses distinct from current approaches that target viral proteins.

New Microviridae isolated from Sulfitobacter reveals two cosmopolitan subfamilies of single-stranded DNA phages infecting marine and terrestrial Alphaproteobacteria

The Microviridae family represents one of the major clades of single-stranded DNA (ssDNA) phages. Their cultivated members are lytic and infect Proteobacteria, Bacteroidetes, and Chlamydiae. Prophages have been predicted in the genomes from Bacteroidales, Hyphomicrobiales, and Enterobacteriaceae and cluster within the 'Alpavirinae', 'Amoyvirinae', and Gokushovirinae. We have isolated 'Ascunsovirus oldenburgi' ICBM5, a novel phage distantly related to known Microviridae. It infects Sulfitobacter dubius SH24-1b and uses both a lytic and a carrier-state life strategy. Using ICBM5 proteins as a query, we uncovered in publicly available resources sixty-five new Microviridae prophages and episomes in bacterial genomes and retrieved forty-seven environmental viral genomes (EVGs) from various viromes. Genome clustering based on protein content and phylogenetic analysis showed that ICBM5, together with Rhizobium phages, new prophages, episomes, and EVGs cluster within two new phylogenetic clades, here tentatively assigned the rank of subfamily and named 'Tainavirinae' and 'Occultatumvirinae'. They both infect Rhodobacterales. Occultatumviruses also infect Hyphomicrobiales, including nitrogen-fixing endosymbionts from cosmopolitan legumes. A biogeographical assessment showed that tainaviruses and occultatumviruses are spread worldwide, in terrestrial and marine environments. The new phage isolated here sheds light onto new and diverse branches of the Microviridae tree, suggesting that much of the ssDNA phage diversity remains in the dark.

Conserved RNA secondary structure in Cherry virus A 5'-UTR associated with translation regulation

Background

A variety of cis-acting RNA elements with structures in the 5′- or 3′-untranslated region (UTR) of viral genomes play key roles in viral translation. Cherry virus A (CVA) is a member of the genus Capillovirus in the family Betaflexiviridae. It has a positive single-stranded RNA genome of ~ 7400 nucleotides (nt). The length of the CVA 5′-UTR is ~ 100 nt; however, the function of this long UTR has not yet been reported.

Methods

Molecular and phylogenetic analyses were performed on 75 CVA sequences, which could be divided into four groups, and the RNA secondary structure was predicted in four CVA 5′-UTR types. These four CVA 5′-UTR types were then inserted upstream of the firefly luciferase reporter gene FLuc (FLuc), and in vitro translation of the corresponding transcripts was evaluated using wheat germ extract (WGE). Then, in-line structure probing was performed to reveal the conserved RNA structures in CVA-5′UTR.

Results

The four CVA 5′-UTR types appeared to have a conserved RNA structure, and the FLuc construct containing these four CVA 5′-UTR types increased the translation of FLuc by 2–3 folds, suggesting weak translation enhancement activity. Mutations in CVA 5′-UTR suppressed translation, suggesting that the conserved RNA structure was important for function.

Conclusion

The conserved RNA secondary structure was identified by structural evolution analysis of different CVA isolates and was found to regulate translation.

A comparative analysis depicting the disease characteristics and phylogenetic signature of human cytomegalovirus infection in Human Immunodeficiency Virus 1 seropositive patients with end-organ retinitis and gastro-enteric diseases

During advanced HIV infection, Human Cytomegalovirus (HCMV) has been proven to produce devitalizing end-organ diseases (EOD). The interactive co-existence of HIV and HCMV has been reported by many researchers and has been suggested to be linked with a more aggressive disease state. This study has been designed to bring forward an assessment of the clinical risk factors capable of defining the conditions of HCMV induced retinitis and gastro-enteric diseases among HIV1 seropositive patients. We also intended to analyse the phylogenetic variation if any, among the infecting virus types inducing the two separate clinical conditions. The patients were arranged in three different groups; (Group 1 with 26 individuals and group 2 and group 3 with 25 individuals each) based on their current status of HIV and HCMV infections. Serum ELISA, qualitative and quantitative detection of HCMV DNA, Real time mRNA expression study, sequencing, and phylogenetic analysis were performed. All statistical analyses and graphs were exercised using relevant software. We found that in HIV patients with HCMV induced end-organ diseases the components of the CXCL9, 10, 11-CXCR3 chemokine pathway is highly expressed with significant differences existing among patients with retinitis and gastrointestinal disease. We found that the gL gene sequences from the retinitis (HR) group clustered almost separately from that of the gastroenteritis (HG) group in the phylogenetic tree. It may be suggested that a form of natural selection pressure is working on the clinical HCMV strains creating a slight divergence in their phylogenetic lineage thereby helping them adapt to the particular tissue microenvironment they are colonizing.

Secondary Structure of Subgenomic RNA M of SARS-CoV-2

SARS-CoV-2 belongs to the Coronavirinae family. Like other coronaviruses, SARS-CoV-2 is enveloped and possesses a positive-sense, single-stranded RNA genome of ~30 kb. Genomic RNA is used as the template for replication and transcription. During these processes, positive-sense genomic RNA (gRNA) and subgenomic RNAs (sgRNAs) are created. Several studies presented the importance of the genomic RNA secondary structure in SARS-CoV-2 replication. However, the structure of sgRNAs has remained largely unsolved so far. In this study, we probed the sgRNA M model of SARS-CoV-2 in vitro. The presented model molecule includes 5'UTR and a coding sequence of gene M. This is the first experimentally informed secondary structure model of sgRNA M, which presents features likely to be important in sgRNA M function. The knowledge of sgRNA M structure provides insights to better understand virus biology and could be used for designing new therapeutics.

Conformational dynamics and energetics of viral RNA recognition by lab-evolved proteins

The conserved and structured elements in viral RNA genomes interact with proteins to regulate various events in the viral life cycle and have become key targets for developing novel therapeutic approaches. We probe physical interactions between lab-evolved proteins and a viral RNA element from the HIV-1 genome. Specifically, we study the role of an arginine-rich loop in recognition of designed proteins by the viral RNA element. We report free energy calculations to quantitatively estimate the protein/RNA binding energetics, focusing on the mutations of arginine residues involved in recognition of the major groove of RNA by proteins.

A Mathematical Analysis of HDV Genotypes: From Molecules to Cells

Hepatitis D virus (HDV) is classified according to eight genotypes. The various genotypes are included in the HDVdb database, where each HDV sequence is specified by its genotype. In this contribution, a mathematical analysis is performed on RNA sequences in HDVdb. The RNA folding predicted structures of the Genbank HDV genome sequences in HDVdb are classified according to their coarse-grain tree-graph representation. The analysis allows discarding in a simple and efficient way the vast majority of the sequences that exhibit a rod-like structure, which is important for the virus replication, to attempt to discover other biological functions by structure consideration. After the filtering, there remain only a small number of sequences that can be checked for their additional stem-loops besides the main one that is known to be responsible for virus replication. It is found that a few sequences contain an additional stem-loop that is responsible for RNA editing or other possible functions. These few sequences are grouped into two main classes, one that is well-known experimentally belonging to genotype 3 for patients from South America associated with RNA editing, and the other that is not known at present belonging to genotype 7 for patients from Cameroon. The possibility that another function besides virus replication reminiscent of the editing mechanism in HDV genotype 3 exists in HDV genotype 7 has not been explored before and is predicted by eigenvalue analysis. Finally, when comparing native and shuffled sequences, it is shown that HDV sequences belonging to all genotypes are accentuated in their mutational robustness and thermodynamic stability as compared to other viruses that were subjected to such an analysis.

Computational search of hybrid human/SARS-CoV-2 dsRNA reveals unique viral sequences that diverge from those of other coronavirus strains

The role of the RNAi/Dicer/Ago system in degrading RNA viruses has been elusive in mammals in the past, which has prompted authors to think that interferon (IFN) synthesis is essential in this clade, relegating the RNAi defense strategy against viral infection as an accessory function. However, recent publications highlight the existence of abundant viral small interference and micro RNAs (VsiRNAs and VmiRNAs) in both cell-line and whole organism based experiments, indicating a contribution of these molecules in host responses and/or viral replication. We explore the theoretical possibility that RNAi triggered by SARS-CoV-2 might degrade some host transcripts in the opposite direction, although this hypothesis seems counterintuitive. The SARS-CoV-2 genome was therefore computationally searched for exact intrapairing within the viral RNA and exact hybrid pairing with the human transcriptome over a minimum of 20 bases in length. Minimal segments of 20-base lengths of SARS-CoV-2 RNA were found based on the theoretical matching with existing complementary strands in the human host transcriptome. Few human genes potentially annealing with SARS-CoV-2 RNA, including mitochondrial deubiquitinase USP30, the subunit of ubiquitin protein ligase complex FBXO21 and two long noncoding RNAs, were retrieved. The hypothesis that viral-originated RNAi might mediate degradation of host transcriptome messages was corroborated by published high throughput sequencing of RNA from infected tissues and cultured cells, clinical observation and phylogenetic comparative analysis, indicating a strong specificity of these SARS-CoV-2 hybrid pairing sequences for human genomes.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

Non-optimal codon usage preferences of coronaviruses determine their promiscuity for infecting multiple hosts

Circulating animal coronaviruses occasionally infect humans. The SARS-CoV-2 is responsible for the current worldwide outbreak of COVID-19 that has resulted in 2 112 844 deaths as of late January 2021. We compared genetic code preferences in 496 viruses, including 34 coronaviruses and 242 corresponding hosts, to uncover patterns that distinguish single- and 'promiscuous' multiple-host-infecting viruses. Based on a codon usage preference score, promiscuous viruses were shown to significantly employ nonoptimal codons, namely codons that involve 'wobble' binding to anticodons, as compared to single-host viruses. The codon adaptation index (CAI) and the effective number of codons (ENC) were calculated for all viruses and hosts. Promiscuous viruses were less adapted hosts vs single-host viruses (P-value = 4.392e-11). All coronaviruses exploit nonoptimal codons to infect multiple hosts. We found that nonoptimal codon preferences at the beginning of viral coding sequences enhance the translational efficiency of viral proteins within the host. Finally, coronaviruses lack endogenous RNA degradation motifs to a significant degree, thereby increasing viral mRNA burden and infection load. To conclude, we found that promiscuously infecting coronaviruses prefer nonoptimal codon usage to remove degradation motifs from their RNAs and to dramatically increase their viral RNA production rates.

Dynamic Molecular Epidemiology Reveals Lineage-Associated Single-Nucleotide Variants That Alter RNA Structure in Chikungunya Virus

Chikungunya virus (CHIKV) is an emerging Alphavirus which causes millions of human infections every year. Outbreaks have been reported in Africa and Asia since the early 1950s, from three CHIKV lineages: West African, East Central South African, and Asian Urban. As new outbreaks occurred in the Americas, individual strains from the known lineages have evolved, creating new monophyletic groups that generated novel geographic-based lineages. Building on a recently updated phylogeny of CHIKV, we report here the availability of an interactive CHIKV phylodynamics dataset, which is based on more than 900 publicly available CHIKV genomes. We provide an interactive view of CHIKV molecular epidemiology built on Nextstrain, a web-based visualization framework for real-time tracking of pathogen evolution. CHIKV molecular epidemiology reveals single nucleotide variants that change the stability and fold of locally stable RNA structures. We propose alternative RNA structure formation in different CHIKV lineages by predicting more than a dozen RNA elements that are subject to perturbation of the structure ensemble upon variation of a single nucleotide.

Insight into the distinctive paradigm of Human Cytomegalovirus associated intrahepatic and extrahepatic cholestasis in neonates

Human Cytomegalovirus has been implicated as a probable cause for the development of hepatic cholestasis among neonates. Our study tried to ascertain the exact demographic, biochemical and immunological markers to differentially diagnose patients with HCMV associated intrahepatic and extrahepatic cholestasis and also decipher the phylogenetic variability among the viral strains infecting the two groups. A total of 110 neonates collected over a span of 2 years were selected for the study classified into four different groups based on the presence of hepatic cholestasis and active HCMV infection. Our analysis predicted that total Cholesterol, GGT, ALP and TNFα were the only significant biological markers with exact cut-off scores, capable of distinguishing between HCMV associated intrahepatic and extrahepatic cholestasis. We confirmed that in patients belonging to both of these groups, the inflammasome is activated and the extent of this activation is more or less same except for the initial activators NLRP3 and AIM2 respectively. When we performed two separate phylogenetic analyses with HCMV gM and gN gene sequences, we found that in both cases the sequences from the IHC and EHC groups formed almost separate phylogenetic clusters. Our study has shown that the HCMV clinical strains infecting at intrahepatic and extrahepatic sites are phylogenetically segregated as distinct clusters. These two separate groups show different physiological as well as immunological modulations while infecting a similar host.

A cross-sectional study exploring disease characteristics and phylogenetic nature of human cytomegalovirus among infected neonates with congenital nephrotic syndrome

BACKGROUND

Congenital nephrotic syndrome (CNS) is a rare but serious condition which affects neonates and is caused by monogenic defects of glomerular structural proteins or congenital viral infections. Several reports have established a causal relationship between human cytomegalovirus (HCMV) intrauterine infection and CNS, but thorough study assessing parameters has not yet been done.

METHODS

This study aimed to ascertain significant demographic, biochemical, serological, inflammatory and etiological parameters with 12 months follow-up to clinically identify and monitor neonates with HCMV-associated CNS and sought to decipher the phylogenetic nature of infecting strains. Differences between four patient groups (neonates < 4 weeks old) with or without CNS and HCMV infection were compared by unpaired t testing and one-way analysis of variance (ANOVA). Linear regression was performed to assess statistical significance among individual groups. Maximum-likelihood-based phylogenetic analysis was performed with HCMV gH gene sequences to compare clinically isolated and referenced NCBI strains. This was further supported by analysis of effective number of codons (ENc), codon adaptation index (CAI) and mRNA structural variation.

RESULTS

Patients with HCMV-associated CNS were found to have significant variations in many studied parameters compared with controls. The majority of clinical strains formed a separate phylogenetic cluster defining them as somewhat distinct from standard reference strains, which was supported by the other analyses.

CONCLUSION

This study defined parameters for monitoring cases of HCMV-associated CNS, which suggest the possible existence of a selection force acting and rendering these HCMV strains able to infect selective host tissues and cause specific disease types.