Recent advancements in experimental and computational methods for RNA secondary structure detection have revealed the crucial role of RNA structural elements in diverse molecular processes within living cells. It has been demonstrated that the secondary structure of the entire viral genome is often responsible for performing crucial functions in the viral life cycle and also influences virus evolution. To investigate the role of viral RNA secondary structure, alongside experimental techniques, the use of bioinformatics tools is important for analyzing various secondary structure patterns, including hairpin loops, internal loops, multifurcations, external loops, bulges, stems, and pseudoknots. Here, we have introduced a Python package for analyzing RNA secondary structure elements in viral genomes, which includes the recognition of common secondary structure patterns, the generation of descriptive statistics for these structural elements, and the provision of their basic properties. We applied the developed package to analyze the secondary structures of complete viral genomes collected from the literature, aiming to gain insights into viral function and evolution. Both the package and the collection of secondary structures of viral genomes are available at http://github.com/KazanovLab/RNAsselem.

The structural transitions RNAs undergo during trafficking are not well-understood. Here, we used the well-developed yeast Ty1 retrotransposon to provide the first structural model of genome (g) RNA in the nucleus from a retrovirus-like transposon. Through a detailed comparison of nuclear Ty1 gRNA structure with those established in the cytoplasm, virus-like particles (VLPs), and synthesized in vitro, we detected Ty1 gRNA structural alternations that occur during retrotransposition. Full-length Ty1 gRNA serves as the mRNA for Gag and Gag-Pol proteins and as the genome that is reverse transcribed within VLPs. We show that about 60% of base pairs predicted for the nuclear Ty1 gRNA appear in the cytoplasm, and active translation does not account for such structural differences. Most of the shared base pairs are represented by short-range interactions, while the long-distance pairings seem unique for each compartment. Highly structured motifs tend to be preserved after nuclear export of Ty1 gRNA. In addition, our study highlights the important role of Ty1 Gag in mediating critical RNA:RNA interactions required for retrotransposition.

Under the oxidative stress condition, the small RNA (sRNA) Oxys that acts as essential post-transcriptional regulators of gene expression is produced and plays a regulatory function with the assistance of the RNA chaperone Hfq protein. Interestingly, experimental studies found that the N48A mutation of Hfq protein could enhance the binding affinity with OxyS while resulting in defection of gene regulation. But, how the Hfq protein interacts with sRNA Oxys and the origin of the stronger affinity of N48A mutation are both unclear. In this paper, molecular dynamics (MD) simulations were performed on the complex structure of Hfq and OxyS to explore their binding mechanism. The molecular mechanics generalized Born surface area (MM/GBSA) and interaction entropy (IE) method were combined to calculate the binding free energy between Hfq and OxyS sRNA, and the computational result is in excellent correlation with the experimental result. Per-residue decomposition of the binding free energy revealed that the enhanced binding ability of the N48A mutation mainly comes from the increased van der Waals interactions (vdW). This research explores the binding mechanism between Oxys and chaperone protein Hfq, and revealed the origin of the strong binding affinity of N48A mutation. The results provided important insights on the mechanism of gene expression regulation affected by protein mutations.

Starting from December 2019, coronavirus disease 2019 (COVID-19) has emerged as a once-in-a-century pandemic with deadly consequences, which urgently calls for new treatments, cures and supporting apparatuses. Remdesivir was reported by World Health Organization (WHO) as the most promising drug currently available for the treatment of COVID-19. Here, we use molecular dynamics simulations and free energy perturbation methods to study the inhibition mechanism of remdesivir to its target SARS-CoV-2 virus RNA-dependent RNA polymerase (RdRp). In the absence of a crystal structure of the SARS-CoV-2 RdRp, we first construct the homology model of this polymerase based on a previously available structure of SARS-CoV NSP12 RdRp (with a sequence identify of 95.8%). We then build the putative binding mode by aligning the remdesivir + RdRp complex to the ATP bound poliovirus RdRp. The putative binding structure is further optimized with molecular dynamics simulations and demonstrated to be stable, indicating a reasonable binding mode for remdesivir. The relative binding free energy of remdesivir is calculated to be -8.28 ± 0.65 kcal/mol, much stronger than the natural substrate ATP (-4.14 ± 0.89 kcal/mol) which is needed for the polymerization. The ~800-fold improvement in the K_d from remdesivir over ATP indicates an effective replacement of APT in blocking of the RdRp binding pocket. Key residues D618, S549 and R555 are found to be the contributors to the binding affinity of remdesivir. These findings demonstrate that remdesivir can potentially act as a SARS-CoV-2 RNA-chain terminator, effectively stopping its RNA reproduction, with key residues also identified for future lead optimization and/or drug resistance studies.

Rice grassy stunt virus (RGSV) a member of Tenuivirus family, is very potent and destructive which effects rice crop in many countries, particularly China. Non coding RNAs have important functions in development and epigenetic regulation of gene expression in numerous organsisms. There is three type of small non coding RNAs have been found in eukaryotes, which are small interferring RNAs (siRNAs), microRNAs (miRNAs) and piwi interacting RNAs (piRNAs). Small RNAs (sRNAs) origination is from the infecting virus which is known as virus-derived small interfering RNAs (vsiRNAs), has responsibility for RNA silencing in plants. Virus-induced gene silencing (VIGS) is mainly dependent on RNA silencing (RNAi). Interestingly, RNA silencing happens in plants during viral infections. RNAi technique showed significant results in Nephotettix cincticeps. RNAi technique demonstrated the gene silencing of planthopper Nilaparvata lugens. The proteins P5, pcf4, Dnj, psn5, and pn6 act as potential movement proteins and serve as silencing suppressors for RGSV. VsiRNAs originate from dsRNA molecules which require Dicer-like (DCL) proteins, RNA dependent RNA polymerase (RdRP) proteins, and Argonaute (AGO) proteins. RdRP uses ssRNA for perfect RNA amplification process and can also be used for DCL dependent secondary vsiRNA formation. VSRs interfere with the movement of signals during silencing mechanism. Moreover, intercellular movement of viruses is facilitated by virus-encoded movement protein. RNAi is found in many eukaryotes which are related to transcriptional or post-transcriptional regulation by gene suppression. Transcription is bidirectional in ssDNA viruses which are originated from dsRNA molecules. In this review, we highlighted the biology of Rice grassy stunt virus and its insect vector and its silencing suppressors. This work will be helpful for plant virologists to understand the whole biogenesis mechanism for rice viruses especially RGSV.

As a major infectious disease in chickens, Newcastle disease causes considerable economic losses in the poultry industry, especially in developing countries where there is limited access to effective vaccination. Therefore, enhancing resistance to the virus in commercial chickens through breeding is a promising way to promote poultry production. In this study, we investigated gene expression changes at 2 and 6 dpi after infection at day21 with a lentogenic Newcastle disease virus in a commercial egg-laying chicken hybrid using RNA sequencing analysis. By comparing NDV challenged and nonchallenged groups, 526 differentially expressed genes (DEGs) (FDR < 0.05) were identified at 2 dpi, and only 36 at 6 dpi. For the DEGs at 2 dpi, IPA analysis predicted inhibition of multiple signaling pathways in response to NDV that regulate immune cell development and activity, neurogenesis and angiogenesis. Upregulation of Interferon Induced Protein with Tetratricopeptide Repeats 5 (IFIT5) in response to NDV was consistent between the current and most previous studies. Sprouty RTK Signaling Antagonist 1 (SPRY1), a DEG in the current study is located in a significant QTL associated with virus load at 6 dpi in the same population. These identified pathways and DEGs provide potential targets to further study breeding strategy to enhance NDV resistance in chickens.

A method for sequencing a protein from a codon sequence is proposed. An unfolded protein molecule is threaded through a nano-sized pore in an electrolytic cell carboxyl end first and held with a voltage such that only the first residue is exposed in the trans chamber of the cell. A tRNA molecule in trans with matching anticodon for the residue binds itself to the latter in the presence of suitable catalysts. It is then cleaved and transferred to an extended electrolytic cell with N pores, 40 ≤ N ≤ 61, in N individual cis chambers and a single trans chamber. Each pore holds an RNA molecule ending in a unique codon that is held exposed in the trans chamber. In the presence of suitable catalysts the anticodon in the transferred tRNA binds with the codon of a matching RNA molecule. By reversing the voltages in the extended cell every RNA molecule except the one to which the transferred tRNA is bound retracts into its cis chamber, this identifies the residue unambiguously. The detected residue in the first cell is cleaved and the process repeated. Unlike in other nanopore-based methods, it suffices to detect the occurrence of a current blockade without having to measure the pore current precisely. A simplified but more time-consuming version that uses only the first cell is also described. In either case no a priori information about the protein is needed so de novo sequencing is possible. A feasibility analysis of the proposed scheme is presented.

We report a prebiotically relevant solution to the N1-ribosylation of pyrimidine nucleobases, a well-known challenge in the RNA World hypothesis. It is found that the presence of metal cations and clay mineral enables the previously unachievable direct ribosylation of uracil, providing by far the highest yield. Spectroscopy and chromatography analyses confirmed the formation of ribosylated uracil. The method can also be extended to the ribosylation of 2-pyrimidinone. These findings are also compatible with the metal-doped-clay model developed by our lab for the unified route of the selection of ribose and subsequent syntheses of nucleotide and RNA.

Abstract: Since ancient times aging has also been regarded as a disease, and humankind has always strived to extend the natural lifespan. Analyzing the genes involved in aging and disease allows for finding important indicators and biological markers for pathologies and possible therapeutic targets. An example of the use of omics technologies is the research regarding aging and the rare and fatal premature aging syndrome progeria (Hutchinson-Gilford progeria syndrome, HGPS). In our study, we focused on the in silico analysis of differentially expressed genes (DEGs) in progeria and aging, using a publicly available RNA Seq dataset (GEO dataset GSE113957) and a variety of bioinformatics tools. We identified several genes that appear to be involved both in natural aging and progeria (KRT8, KRT18, ACKR4, CCL2, UCP2, ADAMTS15, ACTN4P1, WNT16, IGFBP2). Further analyzing these genes and the pathways involved confirmed their possible roles in aging, suggesting the need for further in vitro and in vivo research. The graphical abstract illustrates the analysis workflow we used and will introduce in the following as an example to demonstrate the power of omics and bioinformatics.

Lung cancer is the most common cause of cancer deaths worldwide, and lung adenocarcinoma (LUAD) is the most common histological subtype. However, the prognostic and predictive outcomes differ because of the heterogeneity of programmed cell death. The purpose of this work is to investigate and develop a cuproptosis-associated lncRNA-based LUAD prediction marker. We firstly performed bioinformatic analysis of the Cuprotosis database and The Cancer Genome Atlas (TCGA) database to obtain 19 cuprotosis-related gene datasets and transcriptional data for LUAD. Univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analysis were utilized to construct cuproptosis-associated lncRNA modes. LUAD patients were thus classified into high-risk and low-risk categories based on prognostic risk values, with a median of It acted as a boundary. Risk models were evaluated and validated using Kaplan-Meier analysis, principal component analysis (PCA), gene set enrichment analysis (GSEA) and nomograms. Utilizing the TCGA-LUAD dataset, we identified seven predicted cuproptosis-associated lncRNAs in tumor microenvironment to create the risk model. 95.54% (214/224) of high-risk category tumor samples included cuproptosis-associated gene alterations, compared to 85.65% (203/237) of low-risk category tumor samples, with TP53 accounting for the bulk of occurrences. According to these findings, risk value was superior to other clinical variables and tumor mutation burden as a predictor of 1-, 3-, and 5-year overall survival (OS). The predictive validity of the cuproptosis-associated lncRNA-based risk model for LUAD is high, and this may have implications for how lung cancer patients are treated individually.

Aluminum (Al) toxicity is the most common factor limiting the growth of alfalfa in acidic soil conditions. Melatonin (MT), a significant pleiotropic molecule present in both plants and animals, has shown promise in mitigating Al toxicity in various plant species. This study aims to elucidate the underlying mechanism by which melatonin alleviates Al toxicity in alfalfa through a combined physiological and transcriptomic analysis. The results reveal that 5 μM MT significantly alleviates Al-induced growth inhibition and preserves the morphology of root tip cells. Moreover, MT application reduces Al accumulation in alfalfa roots and shoots by 28.6% and 27.6%, respectively. Additionally, MT plays a crucial role in scavenging Al-induced excess H2O2 by enhancing the activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), consequently reducing malondialdehyde (MDA) levels. The RNA-seq results reveal that MT application significantly upregulates the expression of xyloglucan endotransglucosylase/hydrolase (XTH) and carbon metabolism-related genes, including those involved in the glycolysis process, as well as sucrose and starch metabolism, suggesting that MT application may mitigate Al toxicity by facilitating the binding of Al to the cell walls, thereby reducing intracellular Al accumulation, and by improving respiration and the content of sucrose and trehalose. Taken together, our study demonstrates that MT alleviates Al toxicity in alfalfa by restoring redox homeostasis, enhancing Al binding to the cell wall, and modulating carbon metabolism. This research advances our understanding of the mechanisms underlying MT's effectiveness in mitigating Al toxicity, offering potential strategies for reducing Al accumulation in alfalfa.

MicroRNAs (miRNAs) are a class of small non-coding RNAs, which function as regulators of gene expression and contribute in numerous physiological processes. However, little is known referring to miRNAs function in insect chemosensation. In the current study, nine small RNA libraries were constructed and sequenced from the antennae of nymphs, adult males and females of Apolygus lucorum. In total, 399 miRNAs were identified including 275 known and 124 novel miRNAs. Known miRNAs were classified into 71 families, amongst which, 23 families were insect-specific. Expression profile analysis showed that miR-7-5p_1 was the most abundant miRNAs in the antennae of A. lucorum. Altogether, 69708 target genes related to biogenesis, membrane and binding activities were predicted for 399 miRNAs. Particularly, 15 miRNAs were found to target 16 olfactory genes. These miRNAs could be involved in regulation of olfactory-associated genes ex-pression. Comparing the antennae of nymphs, adult males and females, 94 miRNAs were found to be differentially expressed. The expression levels of some differentially expressed miRNAs measured by qPCR were consistent with sequencing results. This study provides a global miRNAs transcriptome in the antennae of A. lucorum and valuable information for further investigation on miRNA-mRNA interactions, especially the functions of miRNAs in regulating chemosensation.

Since its initial characterization in 2016, the interferon stimulated gene Shiftless (SHFL) has proven to be a critical piece of the innate immune response to viral infection. SHFL expression stringently restricts the replication of multiple DNA, RNA, and retroviruses with an extraordinary diversity of mechanisms that differ from one virus to the next. These inhibitory strategies include the negative regulation of viral RNA stability, translation, and even the manipulation of RNA granule formation during viral infection. Even more surprisingly, SHFL is the first human protein found to directly inhibit the activity of the -1 programmed ribosomal frameshift, a translation recoding strategy utilized across nearly all domains of life and a several human viruses. Recent literature has shown that SHFL expression also significantly impacts viral pathogenesis in mouse models, highlighting its in-vivo efficacy. To help reconcile the many mechanisms by which SHFL restricts viral replication, we provide here a comprehensive review of this complex ISG, its influence over viral RNA fate, and the implications of its functions on the virus-host arms race for control of the cell.

Repetitive nature of the ribosomal DNA (rDNA) gene makes the sequencing of hundreds copies of mammalian 45S rDNA (about 45-kb per copy) extremely difficult and its assembly is often excluded. Increasing evidence shows that 45S rDNA variations (copy number or single nucleotide), structural ribosomal RNA (rRNA) transcript variants, and non-structural rRNA transcripts (sense and anti-sense long noncoding rRNAs that include promoter rRNAs, and rRNA-derived fragments) play essential roles in mammalian development and diseases. Complete pictures of the hundreds copies of 45S rDNA and their rRNA transcripts require further innovation in sequencing techniques that include bioinformatics. The advancements in mammalian rDNA and rRNA sequencings and the discoveries of novel functions of the rDNA variants and rRNA transcripts are discussed here.

Acute myeloid leukemias are characterized by genetic alterations that lead to a complete or partial block at various stages of myeloid differentiation, leading to continual proliferation of myeloid progenitor cells. These alterations can include a well-documented chromosomal translocation that leads to the synthesis of the PML-retinoic acid receptor fusion protein that leads to a block in differentiation and the ability to induce differentiation with all-trans retinoic acid and other small molecules. In essence, in large part, the disease is due to stem cell-like cells that are blocked in their ability to undergo terminal differentiation to more mature, non-proliferative cell types. We discovered a long non-coding RNA (lncRNA), lncRNAp53int1, encoded from the first intron of the p53 tumor suppressor gene. Interestingly, lncRNAp53int1 appears to play an important regulatory role in maintaining the undifferentiated and proliferative state in a manner independent from p53 gene expression. lncRNAp53int1 is expressed at high levels in proliferating, undifferentiated promyelocytes and its levels fall significantly during their terminal differentiation to non-proliferating granulocytes or macrophages. Furthermore, forced overexpression of lncRNAp53int1 blocks terminal differentiation. The observation that the expression of this lncRNA is repressed during terminal differentiation of human leukemic cells and its expression blocks differentiation of leukemic myeloid progenitor cells, is likely to have clinical implications with regards to differentiation therapy for myeloid leukemia. Ongoing research aims to identify the molecular partners that it associates with and the signaling pathways that this lncRNA participates in with the long-term goal of identifying molecules that inhibit lncRNAp53int1 activity and serve as potential therapeutic drugs. The goal of this mini-review is to disseminate these novel and interesting findings to readers and to describe where we think this research is headed and the impact it can have on myeloid leukemias.

Intramuscular fat content (IFC) and meat color are vital indicators affecting pork quality. A significant positive correlation between IFC and redness meat color (a* value) indicates that these two traits are likely to be regulated by shared molecular pathways. To identify candidate genes, hub genes, and signaling pathways that regulate these two traits, we measured the IFC and a* value in 147 hybrid pigs and selected individuls with extreme phenotypes for transcriptome analysis. The results revealed 485 and 394 overlapping differentially expressed genes (DEGs) using DESeq2, limma, and edgeR packages affecting IFC and a* value, respectively. Weighted gene co-expression network analysis (WGCNA) identified four modules significantly correlated with IFC and a* value. Moreover, we integrated functional enrichment analysis results based on DEGs, GSEA, and WGCNA conditions to identify candidate genes, and identified 47 and 53 candidate genes affecting IFC and a* value, respectively. The protein-protein interaction (PPI) network analysis of candidate genes showed that 6 and 14 hub genes affecting IFC and a* value, respectively. These genes mainly participated in various pathways related to lipid metabolism and redox reactions. Notably, four crucial hub genes (MYC, SOX9, CEBPB, and PPAGRC1A) were shared for these two traits. After functional annotation of these four hub genes, we hypothesized that the SOX9/CEBPB/PPARGC1A axis could co-regulate lipid metabolism and myoglobin redox response. Further research on these hub genes, especially the SOX9/CEBPB/PPARGC1A axis, will help to understand the molecular mechanism of the co-regulation of IFC and a* value, which will provide the theoretical basis for improving pork quality.

The chitin synthase B gene is a key enzyme in the chitin synthesis of insect peritrophic matrix (PM), affecting insects' feeding behavior. The chitin synthase B gene has been cloned from the transcription library of Heortia vitessoides Moore. RT-qPCR showed that HvChsb was highly expressed in the larval stage of H. vitessoides, especially on the first day of the prepupal stage, and in the midgut of larvae and the abdomen of adults. After starvation treatment, HvChsb was found to be significantly inhibited with time. After 48 h of starvation, the feeding experiment showed that HvChsb increased with the prolongation of the refeeding time. The experimental data showed that feeding affected the expression of HvChsb. HvChsb was effectively silenced by RNA interference, so its function was lost, significantly decreasing the survival rate of H. vitessoides. The survival rate from larval-to-pupal stages was only 43.33%, accompanied by abnormal phenotypes. It can be seen that HvChsb plays a key role in the average growth and development of H. vitessoides.

B chromosomes (B) are supernumerary elements found in many taxonomic groups. Most B chromosomes are rich in heterochromatin and composed of abundant repetitive sequences, especially transposable elements (TEs). Bs origin is generally linked to the A chromosome complement (A). The first report of a B chromosome in African cichlids was on Astatotilapia latifasciata, which can harbor 0, 1 or 2 B chromosomes. Classical cytogenetics studies found high TE content on the species B chromosome. In this study, we aim to understand TE composition and expression on A. latifasciata genome and its relation to the B chromosome. We use bioinformatics analysis to explore TEs genome organization and also their composition on the B chromosome. Bioinformatics findings were validated by fluorescent in situ hybridization (FISH) and real-time PCR (qPCR). A. latifasciata has a TE content similar to other cichlid fishes and several expanded elements on its B chromosome. With RNA sequencing data (RNA-seq) we showed that all major TE classes are transcribed in brain, muscle and male/female gonads. The evaluation of TE expression between B- and B+ individuals showed that few elements have differential expression among groups and expanded B elements were not highly transcribed. Putative silencing mechanisms may the acting on the B chromosome of A. latifasciata to prevent adverse consequences of repeat transcription and mobilization in the genome.

Emerging paradigms in interferon (IFN) biology suggest a dynamic INF induced interactome that extends through broader Interferon Stimulated Gene (ISG)- induction, which implicates interferon- ISG coordinated cross-talk with mRNA processing, post-translational modification and metabolic processes that underlie pathological (viral, autoimmune and tumor biology) and physiological (stem cell regenerative pathways) processes. INF immune responses can also be triggered by endogenous host-derived molecules that are generated in response to cellular stress or hemostasis imbalance to establish tissue repair and regeneration in first place, however, overactivation or lack of countermeasures can result in host tissue damage. The proteases are integral to viral and tumor pathology, and importantly serine proteases TMPRSS2 and trypsin have been identified as important molecular determinants underlying COVID-19 pathology, and emergence of coronaviruses cultured in vitro, respectively. We propose that pathogen associated proteases can act as novel stress-inducers to facilitate viral- competent immunomodulation. We term it as Protease Induced Transcriptomic/ epi-Transcriptomic Reshaping (PITTR) of host cells to counter cellular stress. We present a novel experimental model and our preliminary findings of trypsin- primed Caco-2 cells (CPT) that result in translational halt comparable to cells grown under serum-starvation conditions (CSS). CPT at escalating trypsin concentration (CPT- EC) induce upregulation of selective proteins that majorly map to ribosomal, RNA transport, and spliceosome ribonucleoproteins (RNPs). The inclusion of proinflammatory IL1-b to CPT (CPT- IL) resulted in global overexpression of proteins comparable to Caco-2 cells cultured in growth-factor rich serum conditions (CFBS), indicating a likely de-repression of trypsin- induced translational halt. Caco-2 cells display abortive interferon proteome under differential trypsin conditions (CPT, CPT-EC and CPT-IL), which is marked by complete lack of INF generation despite induction of intermediate ISGs, suggestive of protease (trypsin)- dependent regulation of INF response. Viruses regulate the proteome of stress granules (SGs) that are induced to cope transient translational halt as a central adaptive response to pathogen induced cellular stress. The integral components of SGs include non-translating mRNA, ribonucleoproteins (RNPs) and RNA binding proteins (RBPs), which together form biological condensates through a biophysical process involving weak electrostatic interactions through intrinsically disordered regions in RBPs resulting in liquid- liquid phase separation. We compared the CPT- EC proteome to the Mammalian Stress Granules Proteome (MSGP) database to explore potential RBPs that could possibly regulate INF response (and could act as potential anti-viral targets). Notably, differentially upregulated RNPs and potential RBPs from ISG family including ADAR and PRKRA, and RNA helicases implicated in viral pathogenesis were found to be upregulated in the CPT- EC proteome further strengthening the role of proteases (trypsin) in regulating INF pathways independent of the pathogen. We propose that the supplementation of viable SARS-CoV-2 viral loads to trypsin- primed host cells could recapitulate an infectious disease model, which may closely phenocopy pathogen- driven inflammation and signaling events. Based on the global downregulation of seven SERPINS (serine protease inhibitors) linked to thromboinflammation in our LCMS profiling data, we support the candidature of serine protease inhibitors for protease mediated viral pathologies. COVID-19 is increasingly linked to coagulopathy and resemblance to Neutrophil Extracellular Trap (NET) related thromboinflammatory features; SERPIN A1AT (alpha 1 antitrypsin) being a potent neutrophil- elastase inhibitor and a negative regulator of coagulation complement pathway may be a promising candidate for establishing hemostasis rebalancing in COVID-19 pathology.

Advances in next-generation sequencing have facilitated the discovery of a multitude of long non-coding RNAs (lncRNAs) with pleiotropic functions in cellular processes, disease and viral pathogenesis. It came as no surprise when viruses were also revealed to transcribe their own lncRNAs. Among them, gammaherpesviruses, one of the three subfamilies of the Herpesviridae, code their largest number. These structurally and functionally intricate non-coding (nc) transcripts modulate cellular and viral gene expression to maintain viral latency or prompt lytic reactivation. The lncRNAs allow the virus to escape cytosolic surveillance, sequester and re-localize essential cellular factors and modulate the cell cycle and proliferation. Some viral lncRNAs act as “messenger molecules”, transferring information about viral infection to neighboring cells. This broad range of lncRNA functions is achieved through lncRNA structure-mediated interactions with effector molecules of viral and host origin, including other RNAs, proteins and DNAs. In this review, we discuss examples of gammaherpesvirus-encoded lncRNAs, emphasize their unique structural attributes, and link them to viral life cycle, pathogenesis and disease progression. We will address their potential as novel targets for drug discovery and propose future directions to explore lncRNA structure and function relationship.

The Ro60 autoantigen is an RNA-binding protein, and a core component of nucleocytoplasmic ribonucleoprotein (RNP) complexes. Ro60 is essential in RNA metabolism, cell stress response pathways, and cellular homeostasis. It stabilises and mediates quality control and cellular dis-tribution of YRNAs, retroelement transcripts, and misfolded RNAs. Ro60 transcriptional dysregulation or loss of function can result in the generation and release of RNA fragments from YRNAs and other small RNAs. Small RNA fragments can instigate an inflammatory cascade through endosomal TLRs and cytoplasmic RNA sensors, which typically sense patho-gen-associated molecular patterns (PAMPs) and mount the first line of defence against invading pathogens. However, recognition of host-originating RNA moieties from Ro60 RNP complexes can activate inflammatory response pathways and compromise self-tolerance. Autoreactive B cells may produce antibodies targeting extracellular Ro60 RNP complexes. Ro60 autoantibodies serve as diagnostic markers for various autoimmune diseases including Sjögren's disease (SjD) and systemic lupus erythematosus (SLE), and they may also act as predictive markers for an-ti-drug antibody response among rheumatic patients. Understanding Ro60’s structure, function, and role in self-tolerance can enhance our understanding of the underlying molecular mecha-nisms of autoimmune conditions.

The COVID-19 pandemic stimulated attention to the medicinal applications of messenger RNA (mRNA). mRNA is expected to be applied not only to vaccines but also to regenerative medicines. The purity of mRNA is important for its medicinal application. However, the current mRNA synthesis techniques have problems, e.g., contamination of undesired 5’-uncapped mRNA and double-stranded RNAs. Recently, our group developed a completely capped mRNA synthesis technology that contributed to the progress of mRNA research. The introduction of chemically modified nucleosides, e.g., N1-methylpseudouridine and 5-methylcytidine, have been demonstrated by Dr. Karikó and Dr. Weissman which opened the practical application of mRNA for vaccines and regenerative medicines. Dr. Yamanaka reported the production of Induced Pluripotent Stem Cells (iPSCs) by the introduction of 4 types of genes using a retrovirus vector. iPSCs are widely used for research about regenerative medicines and the preparation of disease models to screen new drug candidates. Of the Yamanaka factors, Klf4 and c-Myc are oncogenes, and there is a tumor development risk if these are integrated into genomic DNA. Therefore, regenerative medicines using mRNA, which has no risk of genome insertion, have attracted attention. In this review, we summarized the synthesis of mRNA and its application for regenerative medicine.

Nicotinamide is an important functional compound and, in the form of nicotinamide adenine dinucleotide (NAD), is used as a cofactor by protein-based enzymes to catalyze redox reactions. In the context of the RNA world hypothesis, it is therefore reasonable to assume that ancestral ribozymes could have used cofactors such as NAD or its simpler analog nicotinamide riboside (NAR) to catalyze redox reactions. The only described example of such an engineered ribozyme uses a nicotinamide moiety bound to the ribozyme through non-covalent interactions. Covalent attachment of NAR to RNA could be advantageous, but the demonstration of such scenarios to date has suffered from the chemical instability of both NAR and its reduced form NARH, making their use in oligonucleotide synthesis less straightforward. Here, we review the literature describing the chemical properties of the oxidized and reduced species of NAR, their synthesis, and previous attempts to incorporate either species into RNA. We discuss how to overcome the stability problem and succeed in generating RNA structures incorporating NAR.

Abstract: RNase Y is an endoribonuclease that governs global mRNA turnover/processing in Ba-cillus subtilis and likely many other bacteria. The enzyme is anchored at the cell membrane, a pseudo-compartmentalization that appears coherent with its task to initiate decay of mRNAs translated predominantly at the cell periphery. However, the reason and consequence for RNase Y attachment to the membrane are largely unknown. We have studied a strain expressing wild-type levels of a cytoplasmic form of RNase Y from the chromosomal locus. This strain exhibits a slow growth phenotype similar to an RNase Y null mutant. Genome-wide data show a significant effect on the expression of hundreds of genes. Certain RNA substrates clearly require RNase Y membrane attachment but others do not. We found no correlation between mRNA stabilization and the cellular location or function of encoded proteins. The Y-complex, an RNase Y specificity factor, appears to also recognize the cytoplasmic form of the enzyme which is able to restore wild-type levels of the corresponding transcripts. We propose that membrane attachment of RNase Y is required for the functional interaction with many coding and non-coding RNAs, for limiting the cleavage of certain substrates and potentially to avoid unfavorable competition between RNase Y and other ribonucleases like RNase J that have a similar evolutionarily conserved cleavage speci-ficity.

Comparative metatranscriptomics of bacterial and yeast communities of two milk kefir beverages (MKAA1 and MKAA2) obtained by fermentation with two different frozen stocks of the kefir grain CIDCA AGK1, and differing in rheological features and organic acid amounts production was carried out to figure out the relation between differences in physicochemical and rheological properties of kefir and the microbial active profile associated with each product. The dominance of lactic acid bacteria and yeast and a marginal amount of acetic acid bacteria marked the microbiome. The bacterial families Lactobacillaceae and Streptococcaceae account for almost all bacterial gene transcripts, with Lactobacillus helveticus, L. kefiranofaciens, L. gallinarum and Lactococcus lactis being most frequent in the microbiome of MKAA1 beverage, and L. kefiranofaciens, Lc. lactis and Leuconostoc mesenteroides being most prevalent in the MKAA2. Dipodascaceae and Saccharomycetaceae were the leading yeast families represented by Yarrowia lipolytica, Saccharomyces unisporus, and Kluyveromyces marxianus. The MKAA1 and MKAA2 shared >75% KEGG Ortologs (KO) in bacteria and yeast libraries. Considerable decreases in total expressed genes (KEGG Ortolog) assigned to Lactobacillus helveticus and L. gallinarum might be related to the variations in rheological features of the beverages, probably by compromising the interrelations with L. kefiranofaciens which might explain the variations in rheological features of the beverages.

Equine herpesvirus-1 (EHV-1) causes respiratory diseases, abortion, and encephalomyelitis in horses. EHV-1 immediate-early (IE) protein, essential for viral replication, is transactivated by the binding of a multiprotein complex including the open reading frame 12 (ORF12) and some host factors to the IE promoter region. Promoter-associated non-coding RNAs (pancRNAs), which are transcribed from bidirectional promoters, regulate the transcription of neighboring genes in mammals and pathogens. In this study, we identified a novel pancRNA transcribed from across the areas of the 5’-untranslated region and promoter of EHV-1 IE and named it IE pancRNA. IE pancRNA and mRNA were simultaneously expressed in EHV-1-infected RN33B-A68B2M cells. This pancRNA was also transcribed in RK13 and E. Derm cells, which are highly susceptible to EHV-1 infection. Furthermore, IE pancRNA upregulated IE gene expression in the presence of ORF12, and stable expression of IE pancRNA increased the number of EHV-1-infected RN33B-A68B2M cells. These results suggest that IE pancRNAs facilitate EHV-1 proliferation by promoting IE gene expression.

Inflammatory activation within the brain was linked to a decrease in cognitive abilities, however the molecular mechanisms implicated in the development of inflammatory-related cognitive dysfunction and its prevention are poorly understood. This study compared responses of hippocampal transcriptomes 3 months after the striatal infusion of lipopolysaccharide, alone (LPS; 30 µg) resulting in memory loss, or with dexamethasone (DEX; 5 mg/kg intraperitoneal) pretreatment, which abolished the long-term LPS-induced memory impairment. After LPS alone, a significant elevation in the expression of immunity/inflammatory-linked genes, including chemokine (Cxcl13), cytokines (Il1b, Tnfsf13b), and major histocompatibility complex (MHC) class II members (Cd74, RT1-Ba, RT1-Bb, RT1-Da, RT1-Db1) was revealed. DEX pretreatment did not change expression of these genes, but significantly affected expression of genes encoding ion channels, primarily calcium and potassium channels, regulators of glutamate (Slc1a2, Grm5, Grin2a) and GABA (Gabrr2, Gabrb2) neurotransmission which enriched in such GO biological processes as “Regulation of transmembrane transport”, “Cognition”, “Learning”, “Neurogenesis”, and “Nervous system development”. Taken together, the data suggest that: (1) pretreatment with DEX did not markedly affect LPS-induced prolonged inflammatory response; (2) DEX pretreatment can affect processes associated with glutamatergic signaling and nervous system development, possibly involved by that in the recovery of memory impairment induced by LPS

Kiwifruit is a perennial horticultural crop species with high nutritional and economic value. However, various pathogens stress brings a serious threat to kiwifruit. Pentatricopeptide repeat proteins (PPRs) are characterized by tandem repetitions of 35 amino acid motifs, and have been found to function in plant RNA editing. However, the roles of PPRs in plant development and disease resistance remain unclear. In this study, we performed a genome-wide identification and characterization of PPR gene family in two kiwifruit species (Actinidia chinensis (Ach), and Actinidia eriantha (Ace)) with markedly different disease resistance. A total of 497 and 499 PPRs were identified in Ach and Ace, respectively. All the kiwifruit PPRs could be phylogenetically divided into four subfamilies. There were about 40.68% PPRs predicted to be localized to mitochondria or chloroplasts. Synteny analysis showed that the expansion of kiwifruit PPRs mainly came from segmental duplication. Based on RNA-seq data from the fruit during 12 periods of development and maturity, weighted correlation network analysis suggested that two PPRs Actinidia20495.t1 and Actinidia15159.t1 were involved in the fruit development and maturation. In addition, we observed different gene expression of PPRs and chloroplast RNA editing profiles between resistant and susceptible kiwifruits after pathogen infection, indicating the roles of PPRs in stress response by modulating the editing extend of mRNA. The differentially expressed upstream transcription factors of PPRs were further identified, they may regulate resistance adaption by modulating the PPRs expression. These results provide a solid foundation for further analyses of the functions and molecular evolution of PPRs, in particular, for clarifying the resistance mechanisms in kiwifruits and breeding new cultivars with high resistance.

Background: It is assumed that the prevalence of hepatitis D in HBsAg-positive individuals reaches 4.5-13% in the world, and on average about 3% in Europe. Data from several European countries, including Slovakia, are missing or are from an older period. Methods: We analyzed all available data on hepatitis D from Slovakia, including reports from the Slovak Public Health Authority and the results of one prospective study, and three smaller surveys. The determination of anti-HDV IgG and IgM antibodies and/or HDV RNA was used to detect hepatitis D. Results: In the years 2005-2022, no confirmed case of acute or chronic HDV infection was reported in Slovakia. The presented survey includes a total of 343 patients, of which 126 were asymptomatic HBsAg carriers, 33 acute hepatitis B, and 184 chronic hepatitis B cases. In a recent prospective study of 206 HBsAg-positive patients who were completely serologically and virologically examined for hepatitis B and D, only 1 anti-HDV IgG positive and no anti-HDV IgM or HDV RNA positive cases were detected. In other smaller surveys, 2 anti-HDV IgG positive patients were found without the possibility of HDV RNA confirmation. In total, only 3 of 329 HBsAg-positive patients (0.91%) tested positive for anti-HDV IgG antibodies and none of 220 tested positive for HDV RNA. Conclusion: The available data show that Slovakia is one of the countries with a very low prevalence of HDV infection reaching less than 1% in HBsAg-positive patients. Routine testing for hepatitis D is lacking in Slovakia, and therefore it is necessary to implement testing of all HBsAg-positive individuals according to international recommendations.

The formation and development of potato tissues and organs is a complex process regulated by a variety of genes and environmental factors. However, the regulatory mechanisms underlying the growth and development are still unclear. In this study, we used autotetraploid potato JC14 as experimental subject to analyze the transcriptome of root, stem and leaf at seedling, tuber formation and tuber expansion stages to explore the spatio-temporal expression pattern of genes and genetic development characteristics. The results identified thousands of differentially expressed genes and KEGG pathway enrichment analysis showed that these genes were mainly involved in defense response and carbohydrate metabolism pathways. A total of 12 co-expressed Gene modules were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and 4 modules were screened out with the highest correlation with potato stem developmental traits. Core genes in the network were further investigated and functionally annotated by computing the connectivity of genes within the module. The results unveiled number of hub genes in stems at different developmental stages, including carbohydrate metabolism related genes, the defense response related genes, and transcription factors. These findings provide important leads for further understanding of the molecular regulation and genetic mechanisms of potato tissue development.

Multisubunit RNA polymerases (RNAP) carry out transcription in all domains of life; during vi-rus infection, RNAPs are targeted by transcription factors encoded by either the cell or the virus, resulting in the global repression of transcription with distinct outcomes for different host-virus combinations. These repressors serve as versatile molecular probes to study RNAP mechanisms, as well as they aid the exploration of druggable sites for the development of new antibiotics. Here, we review the mechanisms and structural basis of RNAP inhibition by the viral repressor RIP and the crenarchaeal negative regulator TFS4, which follow distinct strategies. RIP operates by occluding the DNA-binding channel and mimicking the initiation factor TFB/TFIIB. RIP binds tightly to the clamp and locks it into one fixed position, thereby preventing conformational oscil-lations that are critical for RNAP function as it progresses through the transcription cycle. TFS4 engages with RNAP in a similar manner to transcript cleavage factors such as TFS/TFIIS through the NTP-entry channel; TFS4 interferes with the trigger loop and bridge helix within the active site by occlusion and allosteric mechanisms, respectively. The conformational changes of RNAP described above are universally conserved and are also seen in inactive dimers of eukaryotic RNAPI and several inhibited RNAP complexes of both bacterial and eukaryotic RNA polymer-ases, including inactive states that precede transcription termination. A comparison of target sites and inhibitory mechanisms reveals that proteinaceous repressors and RNAP-specific antibiotics use surprisingly common ways to inhibit RNAP function.

The major threats to the sustainable supply of forest tree products are adverse climate, pests and diseases. Climate change, exemplified by increased drought, poses a unique threat to global forest health. This is attributed to the unpredictable behavior of forest pathosystems, which can favor fungal pathogens over the host under persistent drought stress conditions in the future. Currently, the effects of drought on tree resistance against pathogens are hypothetical, thus research is needed to identify these correlations. Norway spruce (Picea abies) is one of the most economically important tree species in Europe, and is considered highly vulnerable to changes in climate. Dedicated experiments to investigate how disturbances will affect the Norway spruce - Heterobasidion sp. pathosystem are important, in order to develop different strategies to limit the spread of H. annosum s.l. under the predicted climate change. Here, we report a transcriptional study to compare Norway spruce gene expressions to evaluate the effects of water availability and the infection of Heterobasidion parviporum. We performed inoculation studies of three-year-old saplings in a greenhouse (purchased from a nursery). Norway spruce saplings were treated in either high (+) or low (-) water groups: high water group received double the water amount than the low water group. RNA was extracted and sequenced. Similarly, we quantified gene expression levels of candidate genes in biotic stress and jasmonic acid (JA) signaling pathways using qRT-PCR, through which we discovered a unique preferential defense response of H. parviporum-infected Norway spruce under drought stress at the molecular level. Disturbances related to water availability, especially low water conditions can have negative effects on the tree host and benefit the infection ability of the pathogens in the host. From our RNA-seq analysis, 114 differentially expressed gene regions were identified between high (+) and low (-) water groups under pathogen attack. None of these gene pathways were identified to be differentially expressed from both non-treated and mock-control treatments between high (+) and low (-) water groups. Finally, only four genes were found to be associated with drought in all treatments.

Many positive sense RNA viruses, especially those infecting plants, are known to experience stringent, stochastic population bottlenecks inside the cells they invade, but exactly how and why these populations become bottlenecked are unclear. A model proposed ten years ago advocates that such bottlenecks are evolutionarily favored because they cause the isolation of individual viral variants in separate cells. Such isolation in turn allows the viral variants to manifest the phenotypic differences they encode. Recently published observations lend mechanistic support to this model, and prompt us to refine the model with novel molecular details. The refined model, designated Bottleneck, Isolate, Amplify, Select (BIAS), postulates that these viruses impose population bottlenecks on themselves by encoding bottleneck-enforcing proteins (BNEPs) that function in a concentration-dependent manner. In cells simultaneously invaded by numerous virions of the same virus, BNEPs reach the bottleneck-ready concentration sufficiently early to arrest nearly all internalized viral genomes. As a result, very few (as few as one) viral genomes stochastically escape to initiate reproduction. Repetition of this process in successively infected cells isolate viral genomes with different mutations in separate cells. This isolation prevents mutant viruses encoding defective viral proteins from hitchhiking on sister genome-encoded products, leading to the swift purging of such mutants. Importantly, genome isolation also ensures viral genomes harboring beneficial mutations accrue the cognate benefit exclusively to themselves, leading to the fixation of such beneficial mutations. Further interrogation of the BIAS hypothesis promises to deepen our understanding of virus evolution, and inspire new solutions to virus disease mitigation.

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from liver, muscle, hypothalamus, pituitary and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

The pandemic respiratory disease COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), emerged in Wuhan in December 2019 and then spread throughout the world; Italy was the most affected European country. Despite the close pet-human contact, little is known about the predisposition of pets to SARS-CoV-2. Among these, felines are the most susceptible. In this study, a domestic cat with clear symptoms of pneumonia, confirmed by Rx imaging, was found to be infected by SARS-CoV-2 using quantitative RT–qPCR from a nasal swab. This is the first Italian study reporting on the request of the scientific community to focus attention on the possible role of pets as a SARS-CoV-2 reservoir. An important question remains unanswered: did the cat die from SARS-CoV-2 infection?

This study aimed to rank cancers by the strength of relationship between comprehensive mRNA expression of the most harmful or protective genes and patient survival. Using TCGA dataset including RNA-SEQ and clinical data, we investigated not only gene specific prognostic availability, but also comprehensive prognostic availability of prognostic genes filtered by cox coefficient, and ranked cancers by specially designed prognostic indicator. Through Kaplan-Meier plots, we checked that cancers vary in the strength of influence of prognostic genes, and they follow as the rank. Developing treatment with method to reduce or increase expression of biomarkers for specific cancer which ranked bottom, it would be not efficient in high probability. The results of this study can be a scientific evidence for that.

Viroids are known the smallest plant pathogens, and although their genome sequences do not encode proteins, they can cause disease in economically important crops. In order to control viroid diseases and mitigate their damage, genetic resources used for breeding of the viroid-resistant crop have been searched, but the practical resistant trait has not been found in almost all viroid-crop combinations, as well as the tolerant trait. Due to the difficulty in exploiting naturally occurring resistance or tolerance, various effective strategies have been devised to control viroid diseases using non-transforming or transforming techniques. Meanwhile, extensive findings related to viroid resistance and tolerance may lead to confer resistance or tolerance to viroid infection by combining with the recently emerged new technologies (e.g., spray-induced gene silencing and genome-editing technologies), which are believed to be more environmentally viable and acceptable to the general public than previously reported approaches. In particular, some genome-modified crops produced by the latter technology are already on the market. In this review, we comprehensively summarize the current status about investigation of naturally occurring genetic traits for viroid resistance and tolerance, accumulating knowledge about host factors involved in viroid pathogenicity, and various basic technologies developed to try to possible viroid disease control strategies. Furthermore, we discuss prospects and challenges for the achievement of more effective, practical, and sustainable disease control of viroid.

In "genome editing," a type of genetic engineering, DNA is added, taken away, or changed in an organism's genome. Biomedicine, biotechnology, and synthetic biology, just to name a few, have all gained a lot from how often this method is used. Before the editing process can start, a Double-Strand Breaks (DSB) must be made in the DNA at a specific gene. Researchers have made nucleases, which are sometimes called "molecular scissors," to fix this DSB. Transcription activator-like effector nucleases (TALENs), Zinc-finger nucleases (ZFNs), and homing endonucleases are all examples of protein targets that have been studied and changed. The last part gave an overview of these study projects and explained how to make Cas12a for specific uses. The applications of Cas12a have been extensively studied in recent years, but the protein still holds a lot of promise as a treatment and screening tool. In order to give a quick review of CRISPR-Cas12a and its uses, we will briefly talk about the structure and function of the different parts of the reaction pathway that lead up to the catalysis of the target DNA. Cas12a uses a multistep process to make sure that it is selecting the right DNA. This is a good trait for a device that changes the DNA because it makes it less likely that something bad will happen. Even though data shows that a new CRISPR RNAs (crRNA) molecule can stop Cas12a from randomly destroying ssDNA, it may still hurt the host cell while trying to change the genome. The action of Cas12a catalysis, which is a powerful tool for changing the genome, needs to be changed so that it can be more easily controlled and managed. Using what we know about how Cas12a works, we have made mutants that are less active on ssDNA and more active on dsDNA. So far, only Cas9 and Cas12a, which are both part of the CRISPR family, have been used to change the genome. Because of how similar and different these two endonucleases are, CRISPR can now be used for multiple scientific purposes. Cas12a is better than Cas9 because it makes double-strand breaks (DSBs) that favor Homology-Directed Repair (HDR) over Non-homologous end joining (NHEJ) and can change more than one copy of the genome at the same time. Both Cas9 and Cas12a are being changed to make them better at detecting Protospacer Adjacent Motif (PAM) than they are in their natural state. This will make it possible to start focusing on more genes.

Changes in epitranscriptome with N6‐methyladenine (m6A) modification could be involved in the development of multiple diseases, which might be a prevalent modification of messenger RNAs (mRNAs) in eukaryotes. The m6A modification might be performed through the action of methyltransferases, demethylases, and methylation-binding proteins. Importantly, the m6A methylation may be associated with various neurological disorders including Alzheimer's disease (AD), Parkinson's disease (PD), depression, aging related diseases, and/or aging itself. In addition, the m6A methylation might functionally regulate the eukaryotic transcriptome by influencing splicing, export, subcellular localization, translation, stability, and decay of mRNAs. Neurodegenerative diseases may possess a wide variety of phenotypes depending on the neurons that degenerate at occasions. Interestingly, increasing number of evidences have indicated that m6A modification could modulate the expression of autophagy-related genes and promote autophagy in neuronal cells. Oxidative stresses such as reactive oxygen species (ROS) could stimulate the m6A RNA methylation, which may be also related to the regulation of autophagy and/or the development of neurodegenerative diseases. Both m6A modification and autophagy could also play critical roles in regulating health condition of neurons. Therefore, a comprehensive understanding of the m6A and autophagy relationship in human diseases may benefit in developing therapeutic strategies in the future. This paper reviews advances in the understanding of the regulatory mechanisms of m6A modification in the occurrence and development of neurodegenerative diseases and/or aging, discussing the possible therapeutic procedures related to mechanisms of m6A RNA methylation and autophagy.

Non-coding RNA regions have long been considered the "dark matter" of the genome, with their significance in gene regulation remaining enigmatic. However, recent advances in genomics and molecular biology have illuminated the pivotal role that non-coding RNA regions play in the fine-tuned orchestration of hemoglobin expression. Hemoglobin, the oxygen-carrying protein vital for life, is tightly regulated, and any disruption in its production can lead to severe medical conditions, such as hemoglobinopathies. This review delves into the intricate world of non-coding RNA regions, exploring how they act as master conductors in controlling hemoglobin gene expression. We will discuss various classes of non-coding RNAs, including long non-coding RNAs, microRNAs, and circular RNAs, and their specific mechanisms in fine-tuning hemoglobin levels. Additionally, we will highlight emerging therapeutic approaches that target non-coding RNA regions to ameliorate hemoglobin-related disorders. Understanding the critical role of non-coding RNA regions in hemoglobin regulation opens new avenues for precision medicine and gene therapies, offering hope for those affected by hemoglobinopathies and related conditions. This comprehensive review sheds light on the significance of these non-coding RNA regions, demonstrating that the "silent" players in the genome are, in fact, the conductors of harmony in the symphony of hemoglobin regulation.

In the face of rising global demand and unsustainable production methods, cultivated crustacean meat (CCM) is proposed as an alternative means to produce delicious lobster, shrimp and crab products. Cultivated meat requires starting stem cells that may vary in terms of potency, and propensity to proliferate or differentiate into myogenic (muscle-related) tissues. We suggest that regenerating crustacean limbs harbor a range of stem cells suitable for cultivated meat production and present a humane tissue source because they can be accessed non-lethally. To investigate stem cell activity in this model, we have conducted RNA-Seq analysis across six stages of claw regeneration in the crayfish Cherax quadricarinatus (four pre-molt and two post-molt stages), along with histology and real time quantitative PCR (qPCR). Differential expression analysis revealed expression patterns of significantly up- or down-regulated genes in 14 functional categories. Relative expression of select stem cell target genes (including endocrine, myogenic, cell cycle and pluripotency factors) highlighted two growth peaks during pre-molt regeneration and limited stem cell activity after the molt. Histology and qPCR results support the overall finding that pre-molt limb regeneration stages harbor more stem cell activity, likely making them a more prolific cell source for CCM development.

Epitranscriptomic RNA modifications play a crucial role in the posttranscriptional regulation of gene expression. N6-methyladenosine (m6A) is the most prevalent internal modification of eukaryotic RNA and plays a pivotal role in RNA fate. RNA m6A modification is regulated by a group of cellular proteins, methyltransferases (writers) and demethylases (erasers), which add and remove the methyl group from adenosine, respectively. m6A modification is recognized by a group of cellular RNA-binding proteins (readers) that specifically bind to m6A-modified RNA, mediating effects on RNA stability, splicing, transport, and translation. The functional significance of m6A modification of viral and cellular RNA is an active area of virology research. In this review, we summarize and analyze current literature on m6A modification of HIV-1 RNA, the multifaceted functions of m6A in regulating HIV-1 replication, and the role of viral RNA m6A modification in evading the innate immune responses to infection. Furthermore, we briefly discuss future directions and therapeutic implications of mechanistic studies of HIV-1 epitranscriptomic modifications.

Appropriate vaccine selection is crucial in the control of foot-and-mouth disease (FMD). Vaccination can prevent clinical disease and reduces viral shedding, but there is a lack of cross-protection between the seven serotypes and their sublineages, making the selection of an adequately protective vaccine difficult. Since the exact composition of their vaccines is not consistently disclosed by all manufacturers, incompatibility of the strains used for vaccination with regionally circulating strains can cause vaccination campaigns to fail. Here, we present a deep sequencing approach for polyvalent inactivated FMD vaccines that can identify all component strains by their genome sequences. The genomes of all strains of a commercial pentavalent FMD vaccine were de-novo assembled and the vaccine composition determined semi-quantitatively. The genome assembly required high stringency parameters to prevent misassemblies caused by conserved regions of the genome shared by related strains. In contrast, reference-guided assembly is only recommended in cases where the number of strains is previously known and appropriate reference sequences are available. The presented approach can be applied not only to any inactivated whole-virus FMD vaccine, but also to vaccine quality testing in general and allows for better decision-making for vaccines with unknown composition.

Nuclear export of viral mRNAs, is an essential step in the HIV replication cycle. This role is played by a small regulatory protein of HIV-1 called Rev.The N-terminal region of Rev contains an arginine-rich sequence. The arginine-rich motif (ARM) is located between amino acids 38-50 and forms an alpha-helical secondary structure. Expression of the structural proteins of human immunodeficiency virus type 1 requires the direct interaction of multiple copies of the viral Rev protein with its highly structured RNA target sequence, the Rev Response element (RRE). The major viral proteins are not produced if this transport of RNA is stopped. Therefore, knowledge of Rev structure is essential for understanding of its cooperative binding to the RRE, for understanding the mechanism of HIV infection and for the development of antiviral drugs that interfere with Rev’s essential functions and for acknowledgment of good candidate drugs for treatment of AIDS. To understand how REV interact with RRE element of HIV-RNA and its formation of oligomeric complex it is better to characterize the domain wise structure of REV with regard in function of each domain. Due to lack of structural data on Rev no single compound is reported as inhibitor of REV expect antiviral drugs. Identification of a high-affinity RNA-binding site for the human immunodeficiency virus type 1 Rev Protein is much more important. The ARM is a highly specific sequence which allows for the multimerization of Rev and also binding of REV with RNA. Here we are first time exploring the structural characteristics of REV protein both in free form and in complex with RNA at domain function level especially explore the role of ARM motif in REV HIV-1 protein as RNA binding sites by molecular dynamics (MD) simulation and homology modeling studies. Results indicate that the arginine-rich motif (ARM) is crucial in stability of this complex. The residues ARG38, 39, 41, 43, 44, 48, 50, and ASN40 are most interacting with nucleobases of RRE in Crystal structure of Rev and Rev-response-element RNA complex. Our study plays a major role in elaboration of binding of RNA with REV and pave the way for further investigation for therapeutically agent for HIV.

Cells form a bipolar spindle during mitosis to ensure accurate chromosome segregation. Proper spindle architecture is established by a set of kinesin motors and microtubule-associated proteins. In most eukaryotes, kinesin-5 motors are essential for this process, and genetic or chemical inhibition of their activity leads to the emergence of monopolar spindles and cell death. However, these deficiencies can be rescued by simultaneous inactivation of kinesin-14 motors, as they counteract kinesin-5. We conducted detailed genetic analyses in fission yeast to understand the mechanisms driving spindle assembly in the absence of kinesin-5. Here we show that deletion of the nrp1 gene, which encodes a putative RNA-binding protein with unknown function, can rescue temperature sensitivity caused by cut7-22, a fission yeast kinesin-5 mutant. Interestingly, kinesin-14/Klp2 levels on the spindles in the cut7 mutants were significantly reduced by the nrp1 deletion, although the total levels of Klp2 and the stability of spindle microtubules remained unaffected. Moreover, RNA-binding motifs of Nrp1 are essential for its cytoplasmic localization and function. We have also found that a portion of Nrp1 is spatially and functionally sequestered by chaperone-based protein aggregates upon mild heat stress and limits cell division at high temperatures. We propose that Nrp1 might be involved in post-transcriptional regulation through its RNA-binding ability to promote the loading of Klp2 on the spindle microtubules.

Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. Cucumerinum (Foc) of cucumber is a destructive soil-borne disease in cucumber (Cucumis sativus.L). There is limited knowledge on molecular mechanisms to FW resistance-mediated defense responses in cucumber. In this study, metabolome and transcriptome profiling were carried out in two FW, resistant (NR) and susceptible (NS), near isogenic lines (NILs) before and after Foc inoculation. The NILs showed consistent and stable FW resistances in multiple greenhouse and laboratory screening tests. Widely targeted metabolomic analysis identified differentially accumulated metabolites (DMAs) with significantly more accumulation in NR in response to Foc infection including many phenolic acid and flavonoid compounds from the flavonoid biosynthesis pathway. Transcriptome analysis identified differentially expressed genes (DEGs) between the NILs upon Foc inoculation including genes for secondary metabolite biosynthesis in or transcription factor genes regulating the flavonoid biosynthesis pathway. Joint analysis of the metabolomic and transcriptomic data identified DAMs and DEGs closely associated with biosynthesis of phenolic acid and flavonoid DAMs. The association of these compounds with NR-conferred FW resistance was exemplified by in vivo assays in which two phenolic acid compounds, bis (2-ethylhexyl) phthalate and diisooctyl phthalate, and the flavonoid compound gallocatechin 3-O-gallate were shown to have significant inhibitory effect on Foc growth. The antifungal effect of three compounds seems to be a novel finding. We conclude that phenolic acids and flavonoids play an important role in NR mediated FW resistance in cucumber.

RNA editing involves the insertion, deletion or substitution of single nucleotides within a RNA molecule, without altering the DNA sequence. Adenosine to inosine (A-to-I) editing consists of an RNA modification where single adenosines along the RNA sequence are converted into inosines. Such a biochemical transformation is catalyzed by enzymes belonging to the family of adenosine deaminases acting on RNA (ADARs) and occurs either co- or post-transcriptionally. Initially, the A-to-I RNA editing phenomenon was discovered and studied in messenger RNAs (mRNAs), where it can influence RNA splicing and cause the recoding of codon sequences. The employment of more powerful, high-throughput detection methods has recently revealed that A-to-I editing widely occurs in non-coding RNAs, including microRNAs (miRNAs). MiRNAs are a class of small regulatory non-coding RNAs (ncRNAs) acting as translation inhibitors, known to exert relevant roles in controlling cell cycle, proliferation, and cancer development. Indeed, a growing number of recent researches have evidenced the importance of miRNA editing in cancer biology by exploiting various detection and validation methods. Herein, we briefly overview early and currently available A-to-I miRNA editing detection and validation methods and discuss the significance of A-to-I miRNA editing in human cancer.

Blackcurrant leaf chlorosis associated virus (BCLCaV) was detected recently by next-generation sequencing (NGS) and proposed as a new and distinct species in the genus Idaeovirus. Genomic components of BCLCaV that were detected and confirmed include: 1) RNA-1 that is monocistronic and encodes the replicase complex; 2) a bicistronic RNA-2 that encodes a movement protein (MP) and the coat protein (CP) of the virus, with open reading frames (ORF) that overlap by a single adenine (A) nucleotide (nt) representing the third position of an opal stop codon of the MP ORF2a and the first position of the start codon of the CP ORF2b; 3) a subgenomic form of RNA-2 (RNA-3) that contains ORF2b; and 4) a concatenated form of RNA-2 that consists of a complementary and inverted RNA-3 conjoined to the full-length RNA-2. Analysis of NGS-derived paired-end reads revealed the existence of bridge reads encompassing the 3’-terminus and 5’-terminus of RNA-2 or RNA-3 of BCLCaV. The full RNA-2 or RNA-3 could be amplified using outward facing or abutting primers; also RNA-2/RNA-3 could be detected even after three consecutive RNase R enzyme treatments with denaturation at 95 ^oC preceding each digestion. Evidence was obtained indicating that there are circular forms of BCLCaV RNA-2 and RNA-3.

The oral microbiome, studied by sampling the saliva or oral rinse, has been long thought to have diagnostic capacity for head and neck cancers (HNC). However, previous reports on the HNC oral microbiome provide inconsistent results. The aim of this study is to consolidate these datasets and determine the oral microbial composition between HNC patients to healthy and premalignant individuals. We analyzed 16 published head and neck cancer (HNC) short-read 16S rRNA sequencing datasets, specifically targeting the V3V4, V4 and V4V5 regions. These datasets included saliva and oral rinse samples from donors with HNC, as well as from healthy and premalignant donors. Differences in diversities and microbial abundance were determined. HNC saliva displayed lower alpha diversity than healthy donors. In contrast, the opposite trend was observed for oral rinse samples. Beta diversity were largely similar across different patient types. Similar oral phyla were detected for all samples, but proportions were largely dependent on sample type (i.e. saliva or oral rinse) and primer set utilised for 16S rRNA sequencing. Neisseria, Leptotrichia and Megasphaera were elevated in healthy saliva, while Mycoplasma was elevated in HNC saliva. Oral rinse and saliva displayed similar enrichment for Fusobacterium, while Veillonella, Alloprevotella, and Campylobacter had conflicting results. The sparse partial least squares discriminant analysis model performed effectively in discriminating HNC from healthy or premalignant patients using V3V4 saliva (AUC = 0.888), and V3V4 oral rinse (AUC = 0.928), while poor discriminative capacity was observed for V4 saliva (AUC = 0.688). In conclusion, our meta-analysis highlighted the limitations of 16S rRNA sequencing, particularly due to variations across study batches, primer sets (i.e. V3V4, V4), and sample types. Hence, caution should be exercised when interpreting 16S rRNA sequencing results across studies, especially when different primer sets and sample types are used.

Legume-rhizobial symbiosis leads to to the formation of root nodules, where rhizobia reside and develop into bacteroids to reduce nitrogen into ammonium for plant growth, which leaves an opening question as how plant immunity is attenuated in nodules in the presence of large number of foreign bacteria. In Medicago truncatula, mutation in NAD1 (Nodules with Activated Defense 1) only produces necrotic nodules with overactivated immunity, indicating NAD1 is an indispensable component required for suppressing immunity in nodules. In this study, a dual RNA-seq transcriptomic technology was performed to extensively analyze gene expressions in nad1-1 mutant nodules. We indenifited 89 DEGs in symbiotic nitrogen fixation and 89 DEGs in immunity in Medicago truncatula at 10 dpi. Simultaneously, we indenifited 27 DEGs in fix and nif genes in Sinorhizobium meliloti. Then we identified 56 DEGs related to stress, including ROS, NO, and NCR, most of which were upregulated in Sinorhizobium meliloti. Our analyses of nitrogen fixation-defective plant nad1-1 mutants with overactivated defense indicate that host use plant immunity to control massive bacterial colonization during early and late symbiotic stages. Our findings provide insight into the function of NAD1 in improvement or inhibition of plant immune response to maintain symbiosis during rhizobial endosymbiosis.

Circular RNAs (circRNAs) served as covalently closed single-stranded RNAs have been proposed to influence plant development and stress resistance. Grapevine is the most economically valuable fruit crops cultivated worldwide and threaten by various abiotic stresses. Herein, we reported that a circRNA (Vv-circPTCD1) processed from the second exon of a pentatricopeptide repeat family gene PTCD1 was preferentially expressed in leaves and responded to salt and drought but not heat stress in grapevine. Additionally, the second exon sequence of PTCD1 was highly conserved but the biogenesis of Vv-circPTCD1 is species-dependent in plants. It was further found that the overexpressed Vv-circPTCD1 can slightly decreased abundance of the cognate host gene and the neighboring genes were barely affected in grapevine callus. Furthermore, we also successfully overexpressed the Vv-circPTCD1, and found that the Vv-circPTCD1 deteriorated the growth during heat, salt, and drought stresses in Arabidopsis. However, the biological effects on grapevine callus were not always consistent with that of Arabidopsis. Interestingly, we found that the transgenic plants of linear counterpart sequence also conferred the same phenotypes as that of circRNA during the three stress conditions, no matter what species. Those results imply that although the sequences are conserved, the biogenesis and functions of Vv-circPTCD1 are species-dependent. Our results indicate that the plant circRNA function investigation should be conducted in homologous species, which support valuable reference for further plant circRNA studies.

Barth Syndrome, a rare X-linked disorder affecting 1:300,000 live births, results from defects in Tafazzin, an acyltransferase that remodels cardiolipin and is essential for mitochondrial respiration. Barth Syndrome patients develop cardiomyopathy, muscular hypotonia and cyclic neutropenia during childhood, rarely surviving to middle age. At present, no effective therapy exists and downstream transcriptional effects of Tafazzin dysfunction are incompletely understood. To identify novel, cell-specific, pathological pathways that mediate heart dysfunction, we performed single-nucleus RNA-sequencing (snRNA-seq) on wild-type (WT) and Tafazzin-knockout (Taz-KO) mouse hearts. We determined differentially expressed genes (DEGs) and inferred predicted cell-cell communication networks from these data. Surprisingly, DEGs were distributed heterogeneously across the cell types, with fibroblasts, cardiomyocytes, endothelial cells, macrophages, adipocytes and pericytes exhibiting the greatest number of DEGs between genotypes. One differentially expressed gene was detected for the lymphatic endothelial and mesothelial cell types, while no significant DEGs were found in the lymphocytes. A Gene Ontology (GO) analysis of these DEGs showed cell-specific effects on biological processes such as fatty acid metabolism in adipocytes and cardiomyocytes, increased translation in cardiomyocytes, endothelial cells and fibroblasts, in addition to other cell-specific processes. Analysis of ligand-receptor pair expression, to infer intercellular communication patterns, revealed the strongest dysregulated communication involved adipocytes and cardiomyocytes. For the knockout hearts, there was a strong loss of ligand-receptor pair expression involving adipocytes, and cardiomyocyte expression of ligand-receptor pairs underwent reorganization. These findings suggest that adipocyte and cardiomyocyte mitochondria may be most sensitive to mitochondrial Tafazzin deficiency and that rescuing adipocyte mitochondrial dysfunction, in addition to cardiomyocyte mitochondrial dysfunction, may provide therapeutic benefit in Barth Syndrome patients.

DDX3 is an ATP-dependent RNA helicase involved in multiple cellular activities, including RNA metabolism and innate immunity. DDX3 is known to assist the replication of some viruses while restricting some others through direct interaction with the viral proteins. However, the role of DDX3 in the replication of the hepatitis E virus (HEV) is unknown. In this study, DDX3 is shown to interact with the HEV capsid protein and provide an indispensable role in HEV replication. The DDX3 C-terminal domain was demonstrated to interact with the capsid protein, which was previously demonstrated to inhibit the production of type I interferons. Knockdown of DDX3 compromised the capsid protein-mediated blockage of interferon induction. Notably, DDX3 silencing led to a significant reduction in HEV replication. Also, the ATPase activity of DDX3 is required for the HEV replication as an ATPase-null mutant DDX3 failed to rescue the viral replication in the DDX3-silenced cells. These results demonstrate a pro-viral role of DDX3 in HEV replication, providing further insights into the virus-cell interactions.

Remdesivir (RDV) has garnered much hope for its moderate anti-COVID-19 effects, but its limited amelioration of survival in hospitalized patient causes a huge controversy over the applicability of RDV to COVID-19 treatment. Developing strategies to improve its antivirus efficacy is urgently required. As anticipated, RDV exhibits similar behavior with other nucleotide analogs to disrupt the metabolism of natural endogenous ribonucleotides (RNs) and deoxyribonucleotides (dRNs). Alterations in endogenous RNs and dRNs play a critical role in virus replication as well as other key cellular functions. Thus elucidation of the disturbances of RDV on RNs and dRNs could help to understand its exact mechanism of action. Here, the metabolic profiling determined by liquid chromatography–mass spectrometry method showed a general increase in the abundance of nucleotides and a more than 2-fold increase for specific nucleotides. However, the variation of pyrimidine ribonucleotides was relative slight or even contrary, resulting in obvious imbalance between purine and pyrimidine ribonucleotides, which implied the obstacle of RDV to pyrimidine synthesis and could further block the transcription and replication of viral RNA. Additionally, the extreme disequilibrium between cytidine triphosphate (CTP) and cytidine monophosphate might result from the inhibition of CTP synthase and provide a metabolic target for the treatment of COVID-19 infection. Since nucleotides metabolism pathways are vulnerable to nucleotide analogues and are liable to be the regulation targets, it is promising to enhance the efficacy of RDV through co-administration with CTP synthase inhibitors or de novo pyrimidine synthesis inhibitors to exacerbate the imbalance of nucleotide pools.

Gestational Diabetes Mellitus (GDMs) is worldly wide pregnancy complication. Gestational Diabetes can significantly impact the fetus development. However, the effects of high glucose on embryological development post-fertilization are yet to be researched. Danio rerio embryos are great model to study the embryonic development. In this study, the effects on embryological (morphological and genetic) development were examined in the presence of a high glucose environment that mimic the developing fetus in pregnant women with GMD. Fertilized Zebrafish embryos were treated with normal media and high glucose for 5 days from 3 hours post-fertilization(hpf) to 96hpf respectively as control and experimental groups. Morphological changes are recorded with Microscope images. Hatch rate and heart rate are compared between groups in set up time points. RNA-Seq is performed to examine the gene changes in experimental group. Glucose delayed the zebrafish embryo development by slowing the hatch rate about 24 hours. The brain, heart, and tail started showing morphology smaller in the glucose group compared to the control group at 24 hpf. Heart rate is faster in the glucose group compared to the control group on day 2 and 3 with a statistically significant difference. Among zebrafish whole genome, the significantly changed genes are 556 genes up-regulated and 1118 genes down-regulated respectively in high glucose group. The metabolic and Wnt pathways are altered under high glucose condition. These conditions contribute to significant physiological differences that may provide insight on functionality post-embryological development.

The objective of this study was to identify key transcription factors involved in lipid metabolism and immune response related to the differentially expressed genes (DEG) from the liver samples of 35 pig model for metabolic diseases fed diets containing either 1.5 or 3.0% soybean oil (SOY1.5 or SOY3.0). A total of 281 DEG between SOY1.5 and SOY3.0 diets (log2fold-change ≥ 1 or ≤ −1; FDR-corrected p-value < 0.1) were identified, in which 129 were down-regulated and 152 were up-regulated in SOY1.5 group. The functional annotation analysis detected transcription factors linked to lipid homeostasis and immune response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key transcription factors related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive the food technology and science.

The ongoing COVID-19 pandemic dictated new priorities in biomedicine research. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is a single-stranded positive-sense RNA virus. In this pilot study, we optimized our padlock assay to visualize genomic/subgenomic regions using formalin-fixed paraffin-embedded placental samples obtained from a confirmed case of COVID-19. SARS-CoV-2 RNA was localized in trophoblastic cells. We also checked the presence of the virion by immunolocalization of its glycoprotein spike. In addition, we imaged mitochondria of placental villi keeping in mind that the mitochondrion has been suggested as a potential residence of the SARS-CoV-2 genome. Indeed, we observed a substantial overlapping of SARS-CoV-2 RNA and mitochondria in trophoblastic cells. This intriguing linkage correlated with an aberrant mitochondrial network. Overall, to our knowledge, this is the first study that provides the evidence of a co-localization of the SARS-CoV-2 genome and mitochondria in SARS-CoV-2 infected tissue. These findings also support the notion that SARS-CoV-2 infection could reprogram mitochondrial activity in highly specialized maternal/fetal interface.

Background: Jerusalem artichoke (Helianthus tuberosus L.) is tolerant to salinity stress and has high economic value. The salt tolerance mechanisms of Jerusalem artichoke are still unclear. Especially in the early stage of Jerusalem artichoke exposure to salt stress, the plant physiology, biochemistry and gene transcription are likely to undergo large changes. Elucidating these changes may be of great significance to understanding the salt tolerance mechanisms of it. Results: We obtained high-quality transcriptome from leaves and roots of Jerusalem artichoke exposed to salinity (300 mM NaCl) for 0 h, 6 h, 12 h, 24 h and 48 h, with 150,129 unigenes and 9023 DEGs (Differentially Expressed Genes). The RNA-seq data were clustered into time-dependent groups (nine clusters each in leaves and roots); gene functions were distributed evenly among the groups convergence. KEGG enrichment analysis showed the genes related to plant hormone signal transduction were enriched in almost all treatment comparisons. Under salt stress, genes belongs to PYL (abscisic acid receptor PYR / PYL family), PP2C (Type 2C protein phosphatases), GH3 (Gretchen Hagen3), ETR (ethylene receptor), EIN2/3 (ethylene-insensitive protein 2/3), JAZ (Genes such as jasmonate ZIM-domain gene) and MYC2 (Transcription factor MYC2) had extremely similar expression patterns. The results of qPCR of 12 randomly selected genes confirmed the accuracy of RNA-seq. Conclusions: Under the impact of high salinity (300mM) environment, Jerusalem artichoke in the seedling stage was difficult to survive for a long time, and the phenotype was severe in the short term. Based on the expression of genes on the time scale, we found that the distribution of gene functions in time is relatively even. Upregulation of the phytohormone signal transduction had a crucial role in the response of Jerusalem artichoke seedlings to salt stress, the genes of abscisic acid, auxin, ethylene, and jasmonic acid had the most obvious change pattern.

In recent years, the role of non-coding RNAs (ncRNAs) in regulating cell physiology has begun to be better understood. Recent discoveries in viral molecular biology have revealed that such cellular functions are disturbed during viral infections mainly due to host cell ncRNAs, cellular factors, and virus-derived ncRNAs. Apart from the interplay between those molecules, other interactions derive from the specific folding of RNA virus genomes. These fulfill canonical regulation functions such as replication, translation, and viral packaging. In some cases, folds serve as precursors of small viral RNAs whose biogenesis is not yet clearly understood. Since ncRNAs and RNA viral genomes modulate complex molecular and cellular processes in viral infections, a new taxonomy is being proposed here overarching three main categories, considering the current information about ncRNA interactions in some well-known viral infections. The first category shows examples of host ncRNAs associated with the trigger of the immune response under viral infections. The second category describes interactions between the virus and host ncRNAs. The last category shows how the shape of the RNA viral genome is essential in processing RNAs derived from viruses. Finally, we introduce evidence of how these three categories can also work as a framework in order to organize known interactions of ncRNAs and cellular factors under DENV infection. This new taxonomy of interactions provides a comprehensive framework for organizing the ncRNA regulatory roles in the context of viral interactions and an RNA world.

Cesiumtrifluoroacetate (CsTFA) is a key chemical for RNA-based stable isotope analyses to link the structure and function of microbial communities. We report a protocol to easily synthesize CsTFA from Cesiumcarbonate (Cs₂CO₃₎ and Trifluoroacetate (TFA) and show that self-synthesized CsTFA behaves similar to commercial CsTFA in the separation of isotopically labelled and unlabelled E. coli RNA.

The rise of data science in biology stimulates interdisciplinary collaborations to address fundamental questions. Here, we report the outcome of the first SINFONIA symposium focused on revealing the mechanisms governing plant reproductive development across biological scales. The intricate and dynamic target networks of known regulators of flower development remain poorly understood. To analyze development from the genome to the final floral organ morphology, high-resolution data that capture spatiotemporal regulatory activities are necessary and require advanced computational methods for analysis and modeling. Moreover, frameworks to share data, practices and approaches that facilitate the combination of varied expertise to advance the field are called for. Training young researchers in interdisciplinary approaches and science communication offers the opportunity to establish a collaborative mindset to shape future research.

Single cell RNA-seq (scRNA-seq) profiles conceal temporal and spatial tissue developmental information. De novo reconstruction of single cell temporal trajectory has been fairly addressed, but reverse engineering single cell 3D spatial tissue localization is hitherto landmark based, and de novo spatial reconstruction is a compelling computational open problem. Here we show that a new algorithm - named D-CE - for coalescent embedding of single cell transcriptomic networks can address this open problem. We rely merely on the spatial information encoded in the expression patterns of developmental signal transcription factor (DST) genes, and we find that D-CE of cell-cell association DST-transcriptomic networks reliably reconstructs the Geo-seq or single cell samples’ 3D spatial tissue distribution. Comparison to the novoSpaRC and CSOmap (recent and only available de novo 3D spatial reconstruction methods) on 16 datasets and 681 reconstructions, reveals a significantly distinctive superior performance of D-CE.

Although rRT-PCR is the gold standard method for SARS-CoV-2 detection, some factors, such as amplification inhibitors presence, lead to false-negative results. Here we describe differences between rRT-PCR results for SARS-CoV-2 infection in normal and diluted samples, simulating the need for dilution due to amplification inhibitors presence. Viral RNA extraction of nasopharyngeal swabs samples from 20 patients previously detected as 'Negative' and 21 patients detected as 'Positive' for SARS-CoV-2 was realized with the EasyExtract DNA-RNA (Interprise®). rRT-PCR was realized with OneStep/COVID-19 (IBMP) kit with normal and diluted (80µl of H₂O RNAse free) samples, totaling 82 tests. The results indicate that there is an average variation (ɑ < 0.05) delaying Cq between the amplification results of internal control (IC), N Gene (NG), and ORF-1ab (OF) of 1.811 Cq, 3.840 Cq, and 3.842 Cq, respectively. The extraction kit does not completely purify the inhibitor compounds, therefore non-amplification by inhibitors may occur. In this study, we obtained a 19.04% false-negative diagnosis after sample dilution, and this process reduces the efficiency of rRT-PCR to 29.80% for detecting SARS-CoV-2. Knowing the rRT-PCR standards of diluted samples can help in the identification of false-negative cases, and consequently avoid a wrong diagnosis.

Long noncoding RNA (lncRNA) are implicated in various genetic diseases and cancer, attributed to their critical role in gene regulation. RNA sequencing is used to capture their transcripts from certain cell types or conditions. For some studies, lncRNA interactions with other biomolecules have also been captured, which can give clues to their mechanisms of action. Complementary \textit{in silico} methods have been proposed to predict non-coding nature of transcripts and to analyze available RNA interaction data. Here we provide a critical review of such methods and identify associated challenges. Broadly, these can be categorized as reference-based and reference-free or \textit{ab initio}, with the former category of methods requiring a comprehensive annotated reference. The \textit{ab initio} methods can make use of machine learning classifiers that are trained on features extracted from sequences, making them suitable to predict novel transcripts, especially in non-model species. Machine learning approaches such as Logistic Regression, Support Vector Machines, Random Forest, and Deep Learning are commonly used. Initial approaches relied on basic sequential features to train the model, whereas the use of secondary structural features appears to be a promising approach for functional annotation. However, adding secondary features will result in model complexities, thus demanding an algorithm that can handle it and furthermore, considerably increasing the utilization of computation resources. Computational strategies combining identification and functional annotation which can be easily customized are currently lacking. These can be of immense value to accelerate research in this class of RNAs.

Emergence of zoonotic-human pathogens is proven to be a lethal threat to public health, and RNA virus including influenza viruses, severe acute respiratory syndrome coronavirus, middle east respiratory syndrome coronavirus, Wuhan coronavirus (COVID-19), plays a pivotal role. As those viruses as airborne microorganisms spread mainly by tiny airborne particles, it is important to de-active those airborne particles before their entry into human bodies. In this study, we investigated the effect of far infrared (FIR) radiation on inhibition of airborne microorganisms. The result confirmed that double stand DNA from airborne microorganisms were stable under mild FIR radiation. However, single strand RNA from them was found to be sensitive to FIR radiation, indicating that RNA virus in airborne particles is instable under FIR radiation. Based on this observation, two models on usage of FIR radiation to prevent RNA virus transmission and cure RNA virus infection were proposed, implying that FIR radiation might be a cheap, convenient, and efficient method in clinic to treat RNA virus.

Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is one of many severe diseases that threaten bread wheat (Triticum aestivum L.) yield and quality worldwide. The discovery and deployment of powdery mildew resistance genes (Pm) can prevent this disease epidemic in wheat. In a previous study, we transferred the powdery mildew resistance gene Pm57 from Aegilops searsii into common wheat and cytogenetically mapped the gene in a chromosome region with the fraction length (FL) 0.75-0.87, which represents 12% of 2Ss#1 segment on the long arm of chromosome 2Ss#1. In this study, we performed RNA-Seq on three infected and mock-infected wheat-Ae. searsii 2Ss#1 introgression lines with Bgt-isolates inoculation at 0, 12, 24, and 48 hours after inoculation. Then we designed 79 molecular markers based on transcriptome sequences and physically mapped them to Ae. searsii chromosome 2Ss#1- in seven intervals. We used these markers to identify 46 wheat-Ae. searsii 2Ss#1 recombinants induced by ph1b, a deletion mutant of pairing homoelogous (Ph) genes. Analysis of the 46 ph1b-induced 2Ss#1L recombinants with different Bgt-responses using 28 2Ss#1L-specific molecular markers in the interval FL0.72-0.87 where Pm57 is located, and the flanking intervals, we physically mapped Pm57 gene on the long arm of 2Ss#1 in a 5.13 Mb genomic region, which was flanked by markers X67593 (773.72 Mb) and X62492 (778.85 Mb). By comparative synteny analysis of the corresponding region on chromosome 2B in Chinese spring (T. aestivum L.) with other model species we identified ten genes that are putative plant defense-related (R) genes which includes six coiled-coil nucleotide-binding site-leucine-rich repeat (CNL), three nucleotide-binding site-leucine-rich repeat (NL) and a leucine-rich receptor-like repeat (RLP) encoding proteins. This study will lay a foundation for further cloning of Pm57, and benefit the understanding of interactions between resistance genes of wheat and powdery mildew pathogens.

Although there is a several array of diagnostic and therapeutic choices for pancreatic cancer in recent years, a crucial medical approach for the refractory disease is still needed. Oligonucleotide therapeutics, such as those based on antisense RNAs, RNA interference, aptamers and decoys, are promising agents against pancreatic cancer because they identify a specific nucleotide sequence or protein and interfere with gene expression as molecular-targeted agents. Within just the past quarter-century, the diversity and feasibility of these drugs as diagnostic or therapeutic tools have dramatically increased. Actually, there have been several clinical and preclinical studies of oligonucleotides for patients with pancreatic cancer so far. To support the discovery of effective diagnostic or therapeutic options by using oligonucleotide-based strategies in the absence of satisfactory therapies for long-term survival and the rising trend of diseases, we summarize the current clinical trials of oligonucleotide therapeutics for pancreatic cancer patients with underlying preclinical or scientific data and focus on the possibility of oligonucleotides to target pancreatic cancer in clinical implications.

Our recent work has focused on the application of infinite group theory and related algebraic geometric tools in the context of transcription factors and microRNAs. We were able to differentiate between “healthy" nucleotide sequences and disrupted sequences that may be associated with various diseases. In this paper, we extend our efforts to the study of messenger RNA metabolism, showcasing the power of our approach. We investigate (i) mRNA translation in prokaryotes and eukaryotes, (ii) polyadenylation in eukaryotes, which is crucial for nuclear export, translation, stability, and splicing of mRNA, (iii) miRNAs involved in RNA silencing and post-transcriptional regulation of gene expression, and (iv) the identification of disrupted sequences that could lead to potential illnesses. To achieve this, we employ (a) infinite (finitely generated) groups fp, with generators representing the r+1 distinct nucleotides and a relation between them [e.g., the consensus sequence in the mRNA translation (i), the poly(A) tail in item (ii), and the miRNA seed in item (iii)], (b) aperiodicity theory, which connects “healthy groups" fp to free groups Fr of rank r and their profinite completion F^r, and (c) the representation theory of groups fp over the space-time-spin group SL2(C), highlighting the role of surfaces with isolated singularities in the character variety. Our approach could potentially contribute to the understanding of the molecular mechanisms underlying various diseases and help develop new diagnostic or therapeutic strategies.

Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder characterized by unexplained left ventricular hypertrophy, with or without left ventricular outflow tract (LVOT) obstruction. Single nuclei RNA-sequencing (snRNA-seq) of both obstructive and nonobstructive HCM patient samples have revealed alterations in communication between various cell types but a direct and integrated comparison between the two HCM phenotypes has not been reported. We performed a bioinformatic analysis of HCM snRNA-seq datasets from obstructive and nonobstructive patient samples to identify differentially expressed genes and distinctive patterns of intercellular communication. Differential gene expression analysis revealed 37 differentially expressed genes, predominantly in cardiomyocytes but also in other cell types, relevant to aging, muscle contraction, cell motility and the extracellular matrix. Intercellular communication was generally reduced in HCM, affecting the extracellular matrix, growth factor binding, integrin binding, PDGF binding and SMAD binding, but with increases in adenylate cyclase binding, calcium channel inhibitor activity, and serine-threonine kinase activity in nonobstructive HCM. Increases in neuron to leukocyte and dendritic cell communication, in fibroblast to leukocyte and dendritic cell communication and in endothelial cell communication to other cell types, largely through changes in expression of integrin-b1 and its cognate ligands, were also noted. These findings indicate both common and distinct physiological mechanisms affecting the pathogenesis of obstructive and nonobstructive HCM and provide opportunities for personalized management of different HCM phenotypes.

The Pashmina and Barbari are two famous goat breeds found in the wide areas of the Indo-Pak region. Pashmina is famous for its long hair-fiber (Cashmere) production while Barbari is not-selected for this trait. So, the mRNA expression profiling in the skin samples of both breeds would be an attractive and judicious approach for detecting putative genes involved in this valued trait. Here, we performed differential gene expression analysis on publicly available RNA-Seq data from both breeds. Out of 44,617,994 filtered reads of Pashmina and 55,995,999 of Barbari which are 76.48% and 73.69% mapped to the ARS1 reference transcriptome assembly respectively. A pairwise comparison of both breeds resulted in 47,159 normalized expressed transcripts while 8,414 transcripts are differentially expressed above the significant threshold. Among these, 4,788 are upregulated in Pashmina while 3,626 transcripts are upregulated in Barbari. Fifty-nine transcripts harbor 57 genes including 32 LOC genes and 24 are annotated genes which were selected on the basis of TMM counts > 500. Genes with ectopic expressions other than uncharacterized and LOC symbol genes are Keratins (KRT) and Keratin Associated Proteins (KRTAPs), CystatinA&6, TCHH, SPRR4, PPIA, SLC25A4, S100A11, DMKN, LOR, ANXA2, PRR9 and SFN. All of these genes are likely to be involved in keratinocyte differentiation, sulfur matrix proteins, dermal papilla cells, hair follicles proliferation, hair curvature, wool fiber diameter, hair transition, hair shaft differentiation and its keratinization. These differentially expressed reported genes are critically valuable for enhancing the quality and quantity of the pashmina fiber and overall breed improvement. This study will also provide important information on hair follicle differentiation for further enrichment analyses and introducing this valued trait to other goat breeds as well.

The regulatory mechanism of PCSK7 gene is still unknown, although its encoded protein PC7 is the most ancient and highly conserved of all proprotein convertases and exhibits enzymatic and non-enzymatic functions in liver triglyceride regulation. Bioinformatics algorithms were used to predict regulatory microRNAs (miRNAs) of PCSK7 expression. This led to the identification of 4 miRNAs, namely miR-125a-5p, miR-143-3p, miR-409-3p, and miR-320a-3p, with potential binding sites on the 3’-untranslated region (3′-UTR) of human PCSK7 mRNA. The expression patterns of these miRNAs and PCSK7 mRNA were assessed in three different cell lines with quantitative polymerase chain reaction (qPCR), which revealed reciprocal expression patterns between the expression levels of the four selected miRNAs and PCSK7. Next, the interactions and effects of these miRNAs on PCSK7 expression levels were investigated via cell-based expression analysis, dual-luciferase assay, and Western blot analysis. The data revealed that PCSK7 mRNA levels decreased in cells transfected with vectors overexpressing miR-125a-5p, miR-143-3p, and miR-409-3p, but not miR-320a-3p. The dual-luciferase assay demonstrated that the above 3 miRNAs could directly interact with putative target sites in PCSK7 3′-UTR and regulate its expression, whereas miR-320-3p exhibited no interaction. Western blot analysis further revealed that overexpression of miR-125a-5p in Huh7 cells inhibits the expression and ability of PC7 to cleave human transferrin receptor 1. Our results support a regulatory role of these miRNAs in PCSK7 expression and function and opens the way to assess their roles in the regulation of PC7 activity in vivo in the development of hepatic steatosis.

High molecular weight glutenin subunits (HMW-GS) and Low molecular weight glutenin subunits (LMW-GS) in mature grains play important roles in the formation of glutenin macropolymer and gluten quality. To characterize the expressed glutenin genes of the bread wheat variety Xinmai 26 during seed development, in this study, we measured the dough rheological properties of mature grains through farinograph and gluten testing system, and revealed its strong gluten quality. The compositions of HMW-GS and LMW-GS were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS). Further, a total of 18 full-length transcripts were obtained by third-generation RNA sequencing, including 5 transcripts of HMW-GS genes and 13 transcripts of LMW-GS genes (8 intact genes and 5 pseudogenes). The deduced protein structures of transcript sequences exhibit the typical HMW-GS and LMW-GS structural characteristics. Moreover, a specific functional marker was developed to make better use of the extra cysteine residue of 1Dx5 subunit. This study provides an efficient method to accurately identify glutenin genes in bread wheat through matching full-length transcripts to their spectrum of glutenin, which will be helpful to understand the contributions of glutenin genes to wheat quality.

UPF3B is a constituent of the classical nonsense-mediated mRNA decay (NMD) pathway that degrades both the aberrant transcripts and a set of physiological transcripts. In higher eukaryotes, UPF3B have significant biochemical functions in diverse cellular processes including NMD and translation. UPF3B plays a crucial role in neuronal development and differentiation. Next-generation sequencing technologies identified several loss-of-function mutations in the UPF3B gene that results in neuro-developmental disorders in humans. To uncover the mechanistic role of UPF3B in neuronal functions, we have generated the UPF3B-knockout mammalian cell line model system using CRISPR-Cas9 gene editing method. RNA-Sequencing Analysis of cellular transcriptome from UPF3B-KO cells identified specific genes involved in cell growth and neuronal functions. Altered expression of genes related to the axon guidance pathway delineated the UPF3B function to regulate the neuron-specific genes. Functional enrichment analysis identified the genes involved in the disorders related to mental health and intellectual disability. Our study has the potential to identify the direct players of intellectual disability and will have broader implications.

Ticks transmit a wide variety of pathogens including bacteria, parasites and viruses. Over the last decade, numerous novel viruses have been described in arthropods, including ticks, and their characterization has provided new insights into RNA virus diversity and evolution. However, little is known about their ability to infect vertebrates. As very few studies have described the diversity of viruses present in ticks from the Caribbean, we implemented an RNA-sequencing approach on Amblyomma variegatum and Rhipicephalus microplus ticks collected from cattle in Guadeloupe and Martinique. Among the viral communities infecting Caribbean ticks, we selected four viruses belonging to the Chuviridae, Phenuiviridae and Flaviviridae families for further characterization and designing antibody screening tests. While viral prevalence in individual tick samples revealed high infection rates, suggesting a high level of exposure of Caribbean cattle to these viruses, no seropositive animals were detected. These results suggest that the Chuviridae- and Phenuiviridae-related viruses identified in the present study are more likely tick endosymbionts, raising the question of the epidemiological significance of their occurrence in ticks, especially regarding their possible impact on tick biology and vector capacity. The characterization of these viruses might open the door to new ways of preventing and controlling tick-borne diseases.

Injuries to the anterior cruciate ligament (ACL) often result in post-traumatic osteoarthritis (PTOA). To better understand the molecular mechanisms behind PTOA development following ACL injury, we profiled ACL injury-induced gene expression changes in knee joints of three mouse strains with varying susceptibility to OA: STR/ort (highly susceptible), C57BL/6 (moderately susceptible) and super-healer MRL/MpJ (not susceptible). Right knee joints of the mice were injured using a non-invasive tibial compression injury model that closely mimics ACL rupture in humans and global gene expression was quantified before and at 1-day, 1-week, and 2-weeks post-injury using RNA-seq. Following injury, STR/ort displayed severe cartilage degeneration while MRL/MpJ had little cartilage damage. Gene expression analysis suggested that prolonged inflammation and elevated catabolic activity in STR/ort injured joints, compared to the other two strains may be responsible for the severe PTOA phenotype observed in this strain. MRL/MpJ had the lowest expression values for several inflammatory cytokines and catabolic enzymes activated in response to ACL injury. Furthermore, we identified several genes highly expressed in MRL/MpJ compared to the other two strains including B4galnt2 and Tpsab1 which may contribute to enhanced healing in the MRL/MpJ. Overall, this study has increased our knowledge of early molecular changes associated with PTOA development.

The effects of hypothermia on neonatal encephalopathy may vary regionally, spatially, and cytopathologically in the gyrencephalic neocortex with manifestations potentially influenced by seizures that alter the severity and distribution of neuropathology. We developed a neonatal piglet survival model of hypoxic-ischemic (HI) encephalopathy and hypothermia with continuous encephalographic (cEEG) monitoring for seizures to study injury in neocortex. Neonatal piglets were randomized to naïve, HI-normothermia (NT), overnight hypothermia (HT) initiated 2 hours after HI, sham-NT, or sham-HT treatments. Some piglets within sham and HI groups received cEEG monitoring during recovery. Survival was 2-7 days (piglets were unmedicated and those with poor recovery and unresolving seizures were euthanized early); there was no differences in survival among groups (p&gt;0.078). Neuropathology was assessed by hematoxylin and eosin staining and immunohistochemistry for RNA Binding FOX-1 Homolog 3 (Rbfox3/NeuN). Normal and ischemic-necrotic neurons were counted (layers II-VI collectively) in somatosensory, motor, and prefrontal cortices, identified by cytoarchitecture and connectomics, and in inferior parietal cortex by layer. Seizure burden was determined. HI-NT piglets had reduced normal/total neuron ratio and increased ischemic-necrotic/total neuron ratio relative to naïve, sham-NT, and sham-HT piglets in anterior and posterior motor and somatosensory cortices. Frontal cortex was vulnerable in the prefrontal lateral bank after HI-NT with ischemic-necrosis. HI-HT piglets had higher normal/total neuron ratios and lower ischemic-necrotic/total neuron ratios than HI-NT piglets in anterior/posterior motor and somatosensory cortices. Total normal neuron density in layer III of inferior parietal cortex was reduced in HI-NT piglets compared to sham piglets and was protected by HT. Laminar analysis of Rbfox3 in somatosensory cortex revealed three types of neurons: Rbfox3-positive/normal, Rbfox3-positive/ischemic-necrotic, and Rbfox3-depleted. HI piglets had increased Rbfox3-depleted/total neuron ratio in layers II and III compared to sham-NT. Cortical neuron Rbfox3-depletion was partly rescued by HT. Seizure burden was more severe in HI piglets compared to sham piglets, but HI-NT and HI-HT piglets were similar. We conclude that 1) the neonatal piglet neocortex has a suprasylvian spatial vulnerability to HI and seizures; 2) HT protects against neuropathology in functionally different regions of the neonatal gyrencephalic neocortex; 3) neurons in neonatal neocortex have a limited cytopathology repertoire seen by H&amp;E staining; 4) higher seizure burden correlates with more ischemic-necrotic neurons in somatosensory cortex; 5) seizure presence associates with damage spread to inferior parietal cortex; 6) seizure presence is insensitive to HT, and 7) Rbfox3 immunophenotyping identifies a novel neuronal RNA splicing protein nuclear-depletion pathology that appears reversible by HT. This work demonstrates that HT protection of neocortex in neonatal HI is topographic and laminar, seizure unmitigating, and restores depleted neuronal RNA splicing factor.

Exosomal microRNAs (miRNAs) from cancer cells play a key role in mediating the oral squamous cell carcinoma (OSCC) microenvironment. The objective aimed to investigate how long non-coding RNA (lncRNA) MEG3 affects OSCC angiogenesis through exosomal miR-421. Global miRNA microarray analysis and quantitative real-time PCR (qRT-PCR) were performed to determine the level of miRNAs in OSCC cells-derived exosomes. Cell migration, invasion, tube formation, immunohistochemistry and hemoglobin concentration were used to study the effects of exosomal miR-421 in angiogenesis. Western blotting was used to determine the expression level of HS2ST1 and VEGFR2-related downstream proteins. MiRNA array and qRT-PCR identified upregulation of miR-421 in OSCC cell-derived exosomes. Furthermore, exosome-miR-421 can be taken up by human umbilical vein endothelial cells (HUVEC) and then targets HS2ST1 through VEGF-mediated ERK and AKT phosphorylation, thereby promoting HUVEC migration, invasion, and tube formation. Additionally, enforced expression of lncRNA MEG3 in OSCC cells reduced exosomal miR-421 levels and then increased HS2ST1 expression, thereby reducing the VEGF/VEGFR2 pathway in HUVEC. Our results demonstrate a novel mechanism by which lncRNA MEG3 can act as a tumor suppressor and regulate endothelial angiogenesis through the exosomal miR-421/HS2ST1 axis, which provides a potential therapeutic strategy for OSCC angiogenesis.

Type 2 Diabetes (T2D) is a major public health concern. We have previously shown that the biochemical functions of the gut microbiome contribute to the postprandial glucose response and built a machine learned (ML) model that can personalize dietary choices using a stool metatranscriptomic test. We integrated this data-driven personalized nutrition process into the Viome Precision Nutrition Program (VPNP), which includes sample collection kits (capillary blood and stool), an AI/ML nutritional recommendation engine, and an app for digital delivery of dietary recommendations. To test the efficacy of VPNP, we conducted a decentralized, blinded, randomized, placebo-controlled trial in a prediabetic and diabetic adult population in the USA, with HbA1c as the primary endpoint. For this pilot trial, we recruited 27 participants and randomized them into two arms. The placebo arm included sample collection, digital delivery of the USDA dietary recommendations, and physical delivery of supplements filled with inert material. The VPNP (interventional) arm included digital delivery of data-driven and personalized dietary recommendations, and physical delivery of precision supplements. The nutritional intervention lasted ~90 days. The ANCOVA analysis method was employed to control for potential confounders measured at baseline. VPNP intervention (n=12) reduced HbA1c more than the placebo (n=15) in a clinically significant manner (HbA1c difference of 0.42%, p = 0.028). We demonstrate that a precision nutrition program based on stool and blood metatranscriptomic data and AI/ML analyses can significantly reduce HbA1c relative to USDA recommended nutrition in a prediabetic and diabetic population. These results suggest that integrating data-driven precision nutrition into the healthcare system could reduce the burden of T2D.

Cytosine (C) to uracil (U) RNA editing is one of the most important post-transcriptional processes, however exploring C-to-U editing events efficiently within the crop mitochondrial genome remains a challenge. An improving predictive RNA editor for crop mitochondrial genomes, iPReditor-CMG, was proposed, which was based on SVM, three common crop mitochondrial genomes and self-sequenced tobacco mitochondrial ATPase. After multi-combination feature extracting, high-dimension feature screening and multi-test independent predicting, the results showed that the average accuracy of intraspecific prediction was 0.85, and the highest value even up to 0.91, which outperformed the previous reference models. While the prediction accuracies were 0.78 between dicotyledons and no more than 0.56 between dicotyledons and monocotyledons, implying a possible similarity in C-to-U editing mechanisms among close relatives. The best model was finally identified with an independent test accuracy of 0.91 and an area under the curve of 0.88, and further suggested that five unreported feature sequences TGACA, ACAAC, GTAGA, CCGTT and TAACA were closely associated with the editing phenomenon. Multiple evaluation findings supported that the iPReditor-CMG could be effectively applied to predict crop mitochondrial editing sites, which may contribute to insight into their recognition mechanisms and even other post-transcriptional events in crop mitochondria.

Alphaviruses, flaviviruses, measles viruses and rhabdoviruses are enveloped single-stranded RNA viruses, which have been engineered as expression vector systems for recombinant protein expression and vaccine development. Due to the presence of non-structural genes encoding the replicase complex, a 200,000-fold amplification of viral RNA occurs in the cytoplasm of infected cells providing extreme transgene expression levels, which is why they are named self-replicating RNA viruses. Expression of surface proteins of pathogens causing infectious disease and tumor antigens provide the basis for vaccine development against infectious diseases and cancer. The self-replicating RNA viral vectors can be administered as replicon RNA, recombinant viral particles, or layered DNA/RNA replicons. Self-replicating RNA viral vectors have been applied for vaccine development against influenza virus, HIV, hepatitis B virus, human papilloma virus, Ebola virus and recently coronaviruses, especially SARS-CoV-2 the causative agent of the COVID-19 pandemic. Measles virus and rhabdovirus vector-based SARS-CoV-2 vaccine candidates have been subjected to clinical trials. Moreover, RNA vaccine candidates based on self-amplifying alphaviruses have also been evaluated in clinical settings. Various cancers such as brain, breast, lung, ovarian, prostate cancer and melanoma have also been targeted for vaccine development. Robust immune responses and protection have been demonstrated in animal models. Clinical trials have shown good safety and target-specific immune responses. Ervebo, the VSV-based vaccine against Ebola virus disease has been approved for human use.

Translation of the genetic code into proteins is the main process across all life and ribosomes are ancient cellular machines uniquely enabling this information transformation. We provide a brief overview of the recent advances in linking the ribosomal structure and evolution. Based on these insights into ribosomal organisation across time, we propose that early replication and protein biosynthesis functions were inseparable and in fact were performed by the same ancient RNA molecule, the riboreplisome. Riboreplisome hypothesis helps to address issues of non-Darwinian evolution and complicated starting point that are characteristic to the RNA world, protein world and RNA:protein mixed co-development theories. We suggest that the riboreplisome is the missing link and a molecular machine connecting chemical and biological evolution paths, by being capable of basic genetic and feature selection functions in a cell- or cell-free setting. The riboreplisome hypothesis allows ease of sequential, genetically uninterrupted emergence and sophistication of the genetic code and its decoding machinery, and provides plausible explanations to the origins of the three main RNA types involved in the decoding: the ribosomal, transfer and messenger RNA. Furthermore, riboreplisome can help explaining the co-evolution of aminoacylation machinery, the driving force behind selective gene transcription and expression, and the cell-like compartmentalisation. While we may never find the original riboreplisome again, we might continue to discover different molecular remnants of its prior existence across the existing biological RNA, which, once identified or resurrected, can be useful in synthetic biology applications.

Nucleobase-containing molecules are compounds essential in biology due to the fundamental role of nucleic acids, and in particular G-quadruplex DNA and RNA in life. Moreover, some molecules different from nucleic acids isolated from different vegetal sources or microorganisms show nucleobase moieties in their structure. Nucleoamino acids and peptidyl nucleosides belong to this molecular class. Closely related to the above, nucleopeptides also known as nucleobase-bearing peptides are chimeric derivatives of synthetic origin and more rarely isolated from plants. Herein the self-assembly properties of a vast number of structures, belonging to the nucleic acid and nucleo amino acid/nucleopeptide family, is explored in the light of the recent scientific literature. Moreover, several technologically relevant properties such as the hydrogelation ability of some of the nucleobase-containing derivatives are reviewed in order to open the door to future experimental investigations of newly-inspired nucleobase-driven hydrogels. Nucleobase-containing molecules, such as mononucleosides, DNA, RNA, quadruplex (G4)-forming oligonucleotides and nucleopeptides, are paramount in gel and hydrogel formation owing to their distinctive molecular attributes and ability to self-assemble in biomolecular nanosystems with the most diverse applications in different fields of biomedicine and nanotechnology. In fact, these molecules and their gels present numerous advantages, underscoring their significance and applicability in both material science and biomedicine. Their versatility, capability for molecular recognition, responsiveness to stimuli, biocompatibility, and biodegradability collectively contribute to their prominence in modern nanotechnology and biomedicine. In this review, we emphasize the critical role of nucleobase-containing molecules of different nature in pioneering novel materials with multifaceted applications highlighting their potential in therapy, diagnostics, and new nanomaterials fabrication as required in order to address the numerous current biomedical and nanotechnological challenges.

Studies on animal virome have mainly concentrated on chordates and medically significant invertebrates, often overlooking sylvatic mosquitoes, constituting the major part of mosquito species diversity. Despite their potential role in arbovirus transmission, the viromes of sylvatic mosquitoes remain largely unexplored. These mosquitoes may also harbor insect-specific viruses (ISVs), affecting arboviral transmission dynamics. The Cerrado biome, known for rapid deforestation and its status as a biodiversity hotspot, offers an ideal setting for investigating mosquito viromes due to potential zoonotic spillover risks from land use changes. This study aimed to characterize the viromes of sylvatic mosquitoes collected from various locations within Minas Gerais state, Brazil. Total RNA was extracted from mosquito pools of Psorophora albipes, Sabethes albiprivus, Sa. chloropterus, Psorophora ferox, and Coquillettidia venezuelensis species, followed by high-throughput sequencing (HTS). Bioinformatic analysis included quality control, contig assembly, and viral detection. Sequencing data analysis revealed 11 near-complete viral genomes (new viruses are indicated with asterisks) across seven viral families and one unassigned genus. These included: Xinmoviridae (Ferox mosquito mononega-like virus* and Albipes mosquito Gordis-like virus*), Phasmaviridae (Sabethes albiprivus phasmavirus*), Lispiviridae (Pedras lispivirus isolate MG), Iflaviridae (Sabethes albiprivus iflavivirus*), Virgaviridae (Buriti virga-like virus isolate MG and Sabethes albiprivus virgavirus 1*), Flaviviridae (Psorophora ferox flavivirus*), Mesoniviridae (Alphamesonivirus cavallyense isolate MG), and the genus Negevirus (Biggie virus isolate MG virus and Coquillettidia venezuelensis negevirus*). Phylogenetic analyses revealed evolutionary relationships with varying degrees of similarity to known sequences in distinct host species. In some cases, new viruses and potential new genera were discovered, emphasizing the importance of studying sylvatic mosquito virome to understand viral diversity and evolution in ecologically sensitive regions. Moreover, the presence of ISVs and potential novel arboviruses underscores the need for ongoing surveillance and control strategies to mitigate the risk of emerging infectious diseases.

Intratumor heterogeneity (ITH) is associated with therapeutic resistance and poor prognosis in cancer patients, and attributed to genetic, epigenetic, and microenvironmental factors. We developed a new computational platform, GATHER, for geostatistical modeling of single cell RNA-seq data to synthesize high-resolution and continuous gene expression landscapes of a given tumor sample. Such landscapes allow GATHER to map the enriched regions of pathways of interest in the tumor space and identify genes that have spatial differential expressions at locations representing specific phenotypic contexts using measures based on optimal transport. GATHER provides new applications of spatial entropy measures for quantification and objective characterization of ITH. It includes new tools for insightful visualization of spatial transcriptomic phenomena. We illustrate the capabilities of GATHER using real data from breast cancer tumor to study hallmarks of cancer in the phenotypic contexts defined by cancer associated fibroblasts.

In this study, we resequenced the mitogenomes of P. megacephalus from different regions using a “primer bridging” approach with Sanger sequencing technologies. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between the distribution of species and their genetic distance. However, our research revealed the discovery of the “missing” ATP8 gene in P. megacephalus as well as three other previously published Polypedates, as a result of manual data curation. Using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and determined that COIII, ND2, ND4, and ND6 were reduced to 0.64±0.24, 0.55±0.34, 0.44±0.21 and 0.65±0.17 under the low-temperature stress (8℃), respectively, as compared with controls (P<0.05). Remarkably, the transcript of long non-coding RNA (LncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8℃). Antisense ND6 gene expression showed a downward trend but this was not significant. These results reveal that modulations of protein-coding mitogenomes genes and lncRNAs of P. megacephalus play a crucial part in the molecular response to cold stress.

With the rapid development of gene therapy technology in recent years, its abuse as a method of sports doping in athletic has become a concern. However, there is still room for improvement in gene doping testing methods, and a robust animal model needs to be developed. Therefore, the purposes of this study were to establish a model of gene doping using recombinant adeno-associated virus vector-9 including the human erythropoietin gene (rAAV9-hEPO), and to establish a relevant testing method. First, it was attempted to establish the model using rAAV9-hEPO on mice. The results showed a significant increase in erythrocyte volume accompanied by an increase in spleen weight, confirming the validity of the model. Next, we attempted to detect proof of gene doping by targeting DNA and RNA. Direct proof of gene doping was detected using a TaqMan-qPCR assay with certain primers/probes. In addition, some indirect proof was identified in RNAs through the combination of a TB Green qPCR assay with RNA sequencing. Taken together, these results could provide the foundation for an effective test for gene doping in human athletes in the future.

Sugarcane mosaic virus (SCMV) is a (genus, Potyvirus; family, Potyviridae) is a widespread, deleterious, and the most damaging pathogen of sugarcane ((Saccharum officinarum L. and Saccharum spp.) which causes a substantial barrier to producing high sugarcane earning. Sugarcane mosaic disease (SCMD) is caused by, single or compound infection of SCMV, disseminated by the several aphid vectors in a non-persistent manner. SCMV has flexuous filamentous particle of 700-750 nm long, which encapsidated a positive-sense, single-stranded RNA molecule of 9575 nucleotides. RNA interference (RNAi)-mediated antiviral innate immunity is an evolutionary conserved, key biological process in eukaryotes and has evolved as an antiviral defence system to interfere with viral genomes for controlling infections in plants. The current study aims to analyze sugarcane (Saccharum officinarum L. and Saccharum spp.) locus-derived microRNAs (sof-miRNAs/ssp-miRNAs) with predicted potential for targeting the SCMV +ssRNA-encoded mRNAs, using a predictive approach that involves five algorithms. The ultimate goal of this research is to mobilize the in silico- predicted endogenous sof-miRNAs/ssp-miRNAs to experimentally trigger the catalytic RNAi pathway and generate sugarcane cultivars to evaluate the potential antiviral resistance surveillance ability and capacity for SCMV. Experimentally validated mature sugarcane (S. officinarum, 2n = 8X = 80) and (S. spp., 2n = 100-120) sof-miRNA/ssp-miRNA sequences (n = 28) were downloaded from the miRBase database and aligned with the SCMV genome (KY548506). Among the 28 targeting mature locus-derived sof-miRNAs/ssp-miRNAs evaluated, one sugarcane miRNA homolog, sof-miR159c, was identified to have a predicted miRNA binding site, at nucleotide position 3847 of the SCMV genome targeting CI ORF. To verify the accuracy of the target prediction accuracy and to determine, whether the sugarcane sof-miRNA/ssp-miRNA could bind the predicted SCMV mRNA target(s), we constructed an integrated Circos plot. A genome-wide in-silico-predicted miRNA-mediated target gene regulatory network has implicated to validate interactions necessary to warrant in vivo analysis. The current work provides valuable computational evidence for the generation of SCMV-resistant sugarcane cultivars.

Trichothecene mycotoxin T-2 toxin is an important environmental pollutant and poses a global threat to human and animal health. T-2 toxin could induce nephrotoxicity, however, the precise molecular mechanism remains unclear. In this study, the mice were intraperitoneally administrated with a single dose of 2 mg/kg T-2 toxin. The kidney function and ultrastructural observation were also assessed and found T-2 toxin exposure could cause kidney damage. Transmission electron microscopy showed a significant swelling and vacuolization of the mitochondria of renal cell. A total of 1122, 58, and 391 differential expression genes (DEGs) that were predominantly expressed in kidney tissues at 1 d, 3 d, and 7 d after T-2 toxin exposure, respectively. Early transcriptomic changes at 1 d after exposure were down-regulation of DEGs involved in the cell cycle, p53 signaling pathway, and cellular senescence, while up-regulation of DEGs referred to the ribosome pathway. The temporal variations in gene expression pattern in T-2 toxin exposure in mice kidney were presented and cellular metabolism was disturbed at the exposure recovery period with 7 d. In conclusion, this study, for the first time, provided a comprehensive comparative transcriptomic analysis of gene regulation in T-2 toxin exposure-induced nephrotoxicity at different temporal periods and explored the nephrotoxicity mechanism of T-2 toxin at the mRNA level.

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and secondary meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperm class are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom, and determined the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

Mycoviruses are widespread across the major taxa of fungi including yeasts. Mycoviruses in the yeast fungus Geotrichum candidum are not well studied with three G. candidum-associated viral species characterized to date, all of which belong to the Totiviridae genus Totivirus. In this study, we report the molecular characteristics of another two totiviruses co-infecting isolate Gc6 of G. candidum. The two totiviruses were tentatively named Geotrichum candidum totivirus 2 (GcTV2-Gc6) and GcTV4-Gc6). Both viruses have typical genome organization of totiviruses comprising two ORFs encoding capsid protein (CP) and RNA-dependent RNA polymerase (RdRp) at the N and C termini, respectively. Genomes of GcTV2-Gc6 and GcTV4-Gc6 are 4592 and 4530 bp long, respectively. Both viruses contain the -1 frameshifting elements and their proteins could be expressed as a single fusion protein. GcTV2-Gc6 is closely related to a totivirus isolated from the same host whereas GcTV4-Gc6 is related to insect-associated totiviruses. The phylogenetic analysis indicated that GcTV2-Gc6 and GcTV4-Gc6 belong to two different sister clades, I-A and I-B, respectively. It is interesting that all viruses identified from G. candidum belong to the genus Totivirus. It is possible that the RNA interference (RNAi) mechanism is not active to suppress totivirus accumulation in isolate Gc6.

Olfaction plays a crucial role in many important behavioral activities of insects, such as finding sexual partners, locating hosts, and selecting oviposition sites. Ionotropic receptors (IRs) play a central role in detecting chemosensory information from the environment and guiding insect behaviors and are potential target genes for pest control. Empoasca onukii Matsuda is a major pest of the tea plant Camellia sinensis (L.) O. Ktze, and seriously influences tea yields and quality. In this study, the ionotropic receptor gene EnouIR25a in E. onukii was cloned, and the expression pattern of EnouIR25a was detected in various tissues. Behavioral responses of E. onukii to volatile compounds emitted by tea plants were determined by the olfactometer bioassay and field trials. To further explore the function of EnouIR25a in olfactory recognition of compounds, RNAi (RNA interference) of EnouIR25a was carried out by ingestion of in vitro synthesized dsRNAs. The CDS length of EnouIR25a was 1266 bp and it encoded a 48.87 kD protein. EnouIR25a was enriched in the antennae of E. onukii. E. onukii were more significantly attracted by 1-phenylethanol at concentration of 100 µl/ml. Feeding with dsRNA-IR25a significantly downregulated the expression level of EonuIR25a, after 3 h of treatment, which disturbed the behavioral responses of E. onukii to 1-phenylethanol at concentration of 100µl/ml. The response rate of E. onukii to 1-phenylethanol was significantly decreased after dsRNA-IR25a treatment for 12 h. In summary, the ionotropic receptor gene EnouIR25a was highly expressed in the antennae of E. onukii and was involved in olfactory recognition of the tea plant volatile 1-phenylethanol. The present study may help us to use the ionotropic receptor gene as a target for the behavioral manipulation of E. onukii in the future.

The innate immune system is the first line of host defense sensing viral infection. Manganese (Mn) has recently been found to be involved in the activation of innate immune DNA sensing cGAS-STING pathway and subsequent anti-DNA virus function. However, whether Mn2+ mediates host defense against RNA viruses is still unclear. In the current study, we demonstrated that Mn2+ exhibited antiviral effects against various animal and human viruses including RNA viruses such as PRRSVs and the VSV as well as a DNA virus such as the HSV1 in a dose-dependent manner. Moreover, cGAS and STING were both investigated in the Mn2+ mediated antiviral roles using the knockout cells made by CRISPR-Cas9 approach. Unexpectedly, the results revealed that neither cGAS knockout nor STING knockout had any effect on Mn2+ mediated antiviral functions. Nevertheless, we verified that Mn2+ promoted the activation of cGAS-STING signaling pathway. These findings suggest that Mn2+ has broad spectrum antiviral activities in a cGAS-STING pathway independent manner. This study also provides significant insights into redundant mechanisms participating in the Mn2+ antiviral functions and also indicates new target for Mn2+ antiviral therapeutics.

The central dogma of molecular biology dictates that, with only a few exceptions, information proceeds from DNA to protein through an RNA intermediate. Examining the enigmatic steps from prebiotic to biological chemistry, we take another road suggesting that primordial peptides acted as template for the self-assembly of the first nucleic acids polymers. Arguing in favour of a sort of archaic “reverse translation” from proteins to RNA, our basic premise is a Hadean Hearth where key biomolecules such as amino acids, polypeptides, purines, pyrimidines, nucleosides and nucleotides were available under different prebiotically plausible conditions, including meteorites delivery, shallow ponds and hydrothermal vents scenarios. Supporting a protein-first scenario alternative to the RNA world hypothesis, we propose the primeval occurrence of short peptides termed “selective amino acid- and nucleotide-matching oligopeptides” (henceforward SANMAOs) that noncovalently bind at the same time the polymerized amino acids and the single nucleotides dispersed in the prebiotic milieu. We describe the chemical features of this hypothetical oligopeptide, its biological plausibility and its virtues from an evolutionary perspective. We provide a theoretical example of SANMAO’s selective pairing between amino acids and nucleosides, simulating a poly-Glycine peptide that acts as a template to build a purinic chain corresponding to the glycine’s extant triplet codon GGG. Further, we discuss how SANMAO might have endorsed the formation of low-fidelity RNA’s polymerized strains, well before the appearance of the accurate genetic material’s transmission ensured by the current translation apparatus.

The pandemic of COVID-19 ravaged most countries and made the healthcare system go for a toss. The impact of the disease is different in each patient and it progresses differently. Based on the severity, the COVID-19 infection is stratified into three main categories- mild, moderate, and severe. In this study, we performed a transcriptomic study of different stages and studied the progression of the disease. The study was based on an Indian population of 28 COVID-19 patients, which were classified into different groups. Our analysis has shown that as the disease progresses, the genes involved in the degranulation of the neutrophils and galactose metabolism increase. Furthermore, we identified the hub proteins in each stage. TB is one of the comorbidities of COVID-19 and a comparative study was done to identify the preserved module of genes in both. Enrichment analysis showed that the members of this module are significantly involved in translation and ribosome synthesis.

Single-cell RNA sequencing data facilitates investigation of cell heterogeneity and subpopulations as well as differentially abundant states however modern single-cell RNA sequencing datasets are growing in size and complexity requiring advances in the bioinformatic methods that analyze them. Many methods exist for each step of analysis including read alignment, normalization, quality control, batch effect correction, imputation and dimensionality reduction. With so many options to choose from at each step of the analysis, benchmarking and a synthesis of the literature on the methods available is necessary to inform biological researchers on the most optimal workflow for their data. Here, recent key methods of analysis are highlighted with a focus on methods that facilitate identification of cell subpopulations and differentially abundant cell states. With a constantly expanding toolset for each step in single-cell RNA sequencing dataset analysis, biological researchers should stay informed to utilize the most applicable methods for their own analyses.

In this paper we establish a topological framework of τ-structures to quantify the evolutionary transitions between two RNA sequence-structure pairs. τ-structures developed here consist of a pair of RNA secondary structures together with a non-crossing partial matching between the two backbones. The loop complex of a τ-structure captures the intersections of loops in both secondary structures. We compute the loop homology of τ-structures. We show that only the zeroth, first and second homology groups are free. In particular, we prove that the rank of the second homology group equals the number γ of certain arc-components in a τ-structure, and the rank of the first homology is given by γ−χ+1, where χ is the Euler characteristic of the loop complex.

It was thought until the 1990s that the eukaryotic translation machinery was unable to translate a circular RNA. However internal ribosome entry sites (IRESs) and m⁶A-induced ribosome engagement sites (MIRESs) were discovered, promoting 5’end-independent translation initiation. Today a new family of so-called “non-coding” circular RNAs (circRNAs) has emerged, revealing the pivotal role of 5’end-independent translation. CircRNAs have a strong impact on translational control via their sponge function, and form a new mRNA family as they are translated into proteins with pathophysiological roles. While there is no more doubt about translation of covalently closed circRNA, the linearity of canonical mRNA is only theoretical: it has been shown for more than thirty years that polysomes exhibit a circular form and mRNA functional circularization has been demonstrated in the 1990s by the interaction of initiation factor eIF4G with poly(A) binding protein. More recently, additional mechanisms of 3’-5’ interaction have been reported, including m⁶A modification. Functional circularization enhances translation via ribosome recycling and acceleration of the translation initiation rate. This update of covalently and non-covalently circular mRNA translation landscape shows that RNA circular shape might be the rule for translation with an important impact on disease development and biotechnological applications.

Thyroid cancer is the most prevalent malignancy of the endocrine system and the ninth most common cancer globally. Despite the advances in the management of thyroid cancer, there are critical issues with the diagnosis and treatment of thyroid cancer that result in the poor overall survival of undifferentiated and metastatic thyroid cancer patients. Recent studies have revealed the role of different non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) that are dysregulated during thyroid cancer development or the acquisition of resistance to therapeutics, and may play key roles in treatment failure and poor prognosis of the thyroid cancer patients. Here, we systematically review the emerging roles and molecular mechanisms of ncRNAs that regulate thyroid tumorigenesis and drug response. We then propose the potential clinical implications of ncRNAs as novel diagnostic and prognostic biomarkers for thyroid cancer.

Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) has claimed nearly 180,000 lives and continues to spread. There are currently no approved medications or vaccines for this new coronavirus. Studies have shown that the positive RNA genome of SARS-CoV-2 contains unique features, including a 12-base sequence inserted between the two subunits of viral receptor protein Spike. This inserted sequence facilitates the cleavage of Spike by the cellular proteases Furin and TMPRSS2, leading to the fusion of virus and host cell membranes. Current studies are mostly focused on the SARS-CoV-2 Spike protein and its interacting cellular proteins ACE2, Furin, and TMPRSS2. RNA structural studies are limited and little is known about the potential impact of the 12-base sequence insert on the secondary structure of SARS-CoV-2 genomic RNA and/or its transcripts. Here, by using local and global RNA secondary structure predictions, we show that the novel 12-base insert of SARS-CoV-2 genome likely induces a major RNA secondary structure change.

The outbreak of SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) has caused severe damage to the world. With the support of classic evolutionary theories and population genetics principles, many studies on the origin of SARS-CoV-2 have revealed encouraging results but meanwhile are still under debate. We are concerned with the validity of applying classic evolutionary theories and formula to the evolution of RNA viruses. We have raised several factors like the RNA replication feature and the RNA modification systems of the hosts, which might jeopardize the validity of the application of classic methods to analyze the SARS-CoV-2 data.