The ViennaRNA web services

Complete mitochondrial genomes of Dactylogyrus crucifer and Dactylogyrus zandti reveal distinct patterns of codon usage within Dactylogyrus

Monogeneans of the genus Dactylogyrus Diesing, 1850, the largest genus in the family Dactylogyridae, mostly parasitize the gills of cyprinoid hosts; however, only 3 Dactylogyrus' mitochondrial genomes (mitogenomes) are studied so far. The aim of this research is to extend our understanding of the mitogenomes of Dactylogyrus. We sequenced the mitogenomes of D. crucifer and D. zandti isolated from Rutilus rutilus and Abramis brama orientalis in northwest China, and then we compared these mitogenomes with other monogeneans. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of D. crucifer and D. zandti and characterized the mitogenomes to understand the gene structure, gene identity, the secondary structures of the 22 tRNA genes, and relative synonymous codon usage. We used the analytic Bayesian Information and Maximum Likelihood methods to determine their associated phylogenetic trees. The mitogenomes of D. crucifer and D. zandti were 14,403 and 18,584 bp, respectively. Organization and positioning of these genes were in accordance with Dactylogyrus lamellatus and Dactylogyrus tuba. The nucleotide composition of Dactylogyridae was different from other families of Monogenea, and the A+T count of genus Dactylogyrus (54 - 58.4 %) was lower than other genus species of the family Dactylogyridea (63.9 - 78.4 %) in protein-coding genes. Dactylogyrus members displayed a codon usage bias. The relative synonymous codon used by Dactylogyrus was not conserved and was lower than other monogeneans. The codon use patterns of closely-related species isolated from closely-related hosts were identical. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrus isolated from host subfamily Leuciscinae formed a sister-group. Our results contributed significantly to an increased database of mitogenomes, more than 50 %, for Dactylogyrus that may help future studies of mitochondrial genes and codon uses for the analysis of monogenean phylogenetics.

Morphology, complete mitochondrial genome, and molecular phylogeny of Rhabdias macrocephalum n. sp. (Nematoda: Rhabdiasidae) from Diploderma splendidum (Reptilia: Agamidae)

Species of the genus Rhabdias Stiles & Hassall, 1905 are common parasitic nematodes occurring in the lungs of amphibians and reptiles worldwide. In the present study, Rhabdias macrocephalum n. sp. is described using integrated morphological methods (light and scanning electron microscopy) and molecular approaches (sequencing of the nuclear 28S and ITS regions, and mitochondrial cox1, cox2, and 12S genes) based on specimens collected from the green striped tree dragon Diploderma splendidum (Barbour & Dunn) (Reptilia: Agamidae) in China. The complete mitochondrial genome of R. macrocephalum n. sp. was sequenced and annotated: it is 14,819 bp in length, including 12 protein coding genes (missing atp8), 22 tRNA genes, 2 rRNA genes and three non-coding regions. The gene arrangement of R. macrocephalum n. sp. is different from all of the currently available mitogenomes of nematodes and represents a novel type of mitochondrial gene arrangement reported in Nematoda. Molecular phylogenetic results based on the ITS + 28S data support the monophyly of Entomelas, Pneumonema, Serpentirhabdias, and Rhabdias, and showed R. macrocephalum n. sp. forming a most basal lineage in Rhabdias.

Genetic compensation between ribosomal protein paralogs mediated by a cognate circular RNA

Inter-regulation between related genes, such as ribosomal protein (RP) paralogs, has been observed to be important for genetic compensation and paralog-specific functions. However, how paralogs communicate to modulate their expression levels is unknown. Here, we report a circular RNA involved in the inter-regulation between RP paralogs RpL22 and RpL22-like during Drosophila spermatogenesis. Both paralogs are mutually regulated by the circular stable intronic sequence RNA (sisRNA) circRpL22(NE,3S) produced from the RpL22 locus. RpL22 represses itself and RpL22-like. Interestingly, circRpL22 binds to RpL22 to repress RpL22-like, but not RpL22, suggesting that circRpL22 modulates RpL22's function. circRpL22 is in turn controlled by RpL22-like, which regulates RpL22 binding to circRpL22 to indirectly modulate RpL22. This circRpL22-centric inter-regulatory circuit enables the loss of RpL22-like to be genetically compensated by RpL22 upregulation to ensure robust male germline development. Thus, our study identifies sisRNA as a possible mechanism of genetic crosstalk between paralogous genes.

Genome-wide profiling of Hfq-bound RNAs reveals the iron-responsive small RNA RusT in Caulobacter crescentus

The alphaproteobacterium Caulobacter crescentus thrives in oligotrophic environments and is able to optimally exploit minimal resources by entertaining an intricate network of gene expression control mechanisms. Numerous transcriptional activators and repressors have been reported to contribute to these processes, but only few studies have focused on regulation at the post-transcriptional level in C. crescentus. Small RNAs (sRNAs) are a prominent class of regulators of bacterial gene expression, and most sRNAs characterized today engage in direct base-pairing interactions to modulate the translation and/or stability of target mRNAs. In many cases, the ubiquitous RNA chaperone, Hfq, contributes to the establishment of RNA-RNA interactions. Although the deletion of the hfq gene is associated with a severe loss of fitness in C. crescentus, the RNA ligands of the chaperone have remained largely unexplored. Here we report on the identification of coding and non-coding transcripts associated with Hfq in C. crescentus and demonstrate Hfq-dependent post-transcriptional regulation in this organism. We show that the Hfq-bound sRNA RusT is transcriptionally controlled by the NtrYX two-component system and induced in response to iron starvation. By combining RusT pulse expression with whole-genome transcriptome analysis, we determine 16 candidate target transcripts that are deregulated, many of which encode outer membrane transporters. We hence suggest RusT to support remodeling of the C. crescentus cell surface when iron supplies are limited.IMPORTANCEThe conserved RNA-binding protein Hfq contributes significantly to the adaptation of bacteria to different environmental conditions. Hfq not only stabilizes associated sRNAs but also promotes inter-molecular base-pairing interactions with target transcripts. Hfq plays a pivotal role for growth and survival, controlling central metabolism and cell wall synthesis in the oligotroph Caulobacter crescentus. However, direct evidence for Hfq-dependent post-transcriptional regulation and potential oligotrophy in C. crescentus has been lacking. Here, we identified sRNAs and mRNAs associated with Hfq in vivo, and demonstrated the requirement of Hfq for sRNA-mediated regulation, particularly of outer membrane transporters in C. crescentus.

Molecular phylogeny of the family Rhabdiasidae (Nematoda: Rhabditida), with morphology, genetic characterization and mitochondrial genomes of Rhabdias kafunata and R. bufonis

BACKGROUND

The family Rhabdiasidae (Nematoda: Rhabditida) is a globally distributed group of nematode parasites, with over 110 species parasitic mainly in amphibians and reptiles. However, the systematic position of the family Rhabdiasidae in the order Rhabditida remains unsolved, and the evolutionary relationships among its genera are still unclear. Moreover, the present knowledge of the mitochondrial genomes of rhabdiasids remains limited.

METHODS

Two rhabdiasid species: Rhabdias kafunata Sata, Takeuchi & Nakano, 2020 and R. bufonis (Schrank, 1788) collected from the Asiatic toad Bufo gargarizans Cantor (Amphibia: Anura) in China, were identified based on morphology (light and scanning electron microscopy) and molecular characterization (sequencing of the nuclear 28S and ITS regions and mitochondrial cox1 and 12S genes). The complete mitochondrial genomes of R. kafunata and R. bufonis were also sequenced and annotated for the first time. Moreover, phylogenetic analyses based on the amino acid sequences of 12 protein-coding genes (PCGs) of the mitochondrial genomes were performed to clarify the systematic position of the family Rhabdiasidae in the order Rhabditida using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic analyses based on the 28S + ITS sequences, were also inferred to assess the evolutionary relationships among the genera within Rhabdiasidae.

RESULTS

The detailed morphology of the cephalic structures, vulva and eggs in R. kafunata and R. bufonis was revealed using scanning electron microscopy (SEM) for the first time. The characterization of 28S and ITS regions of R. kafunata was reported for the first time. The mitogenomes of R. kafunata and R. bufonis are 15,437 bp and 15,128 bp long, respectively, and both contain 36 genes, including 12 PCGs (missing atp8). Comparative mitogenomics revealed that the gene arrangement of R. kafunata and R. bufonis is different from all of the currently available mitogenomes of nematodes. Phylogenetic analyses based on the ITS + 28S data showed Neoentomelas and Kurilonema as sister lineages, and supported the monophyly of Entomelas, Pneumonema, Serpentirhabdias and Rhabdias. Mitochondrial phylogenomic results supported Rhabdiasidae as a member of the superfamily Rhabditoidea in the suborder Rhabditina, and its occurrance as sister to the family Rhabditidae.

CONCLUSIONS

The complete mitochondrial genome of R. kafunata and R. bufonis were reported for the first time, and two new gene arrangements of mitogenomes in Nematoda were revealed. Mitogenomic phylogenetic results indicated that the family Rhabdiasidae is a member of Rhabditoidea in Rhabditina, and is closely related to Rhabditidae. Molecular phylogenies based on the ITS + 28S sequence data supported the validity of Kurilonema, and showed that Kurilonema is sister to Neoentomelas. The present phylogenetic results also indicated that the ancestors of rhabdiasids seem to have initially infected reptiles, then spreading to amphibians.

Discovery of novel replication proteins for large plasmids in cyanobacteria and their potential applications in genetic engineering

Numerous cyanobacteria capable of oxygenic photosynthesis possess multiple large plasmids exceeding 100 kbp in size. These plasmids are believed to have distinct replication and distribution mechanisms, as they coexist within cells without causing incompatibilities between plasmids. However, information on plasmid replication proteins (Rep) in cyanobacteria is limited. Synechocystis sp. PCC 6803 hosts four large plasmids, pSYSM, pSYSX, pSYSA, and pSYSG, but Rep proteins for these plasmids, except for CyRepA1 on pSYSA, are unknown. Using Autonomous Replication sequencing (AR-seq), we identified two potential Rep genes in Synechocystis 6803, slr6031 and slr6090, both located on pSYSX. The corresponding Rep candidates, Slr6031 and Slr6090, share structural similarities with Rep-associated proteins of other bacteria and homologs were also identified in various cyanobacteria. We observed autonomous replication activity for Slr6031 and Slr6090 in Synechococcus elongatus PCC 7942 by fusing their genes with a construct expressing GFP and introducing them via transformation. The slr6031/slr6090-containing plasmids exhibited lower copy numbers and instability in Synechococcus 7942 cells compared to the expression vector pYS. While recombination occurred in the case of slr6090, the engineered plasmid with slr6031 coexisted with plasmids encoding CyRepA1 or Slr6090 in Synechococcus 7942 cells, indicating the compatibility of Slr6031 and Slr6090 with CyRepA1. Based on these results, we designated Slr6031 and Slr6090 as CyRepX1 (Cyanobacterial Rep-related protein encoded on pSYSX) and CyRepX2, respectively, demonstrating that pSYSX is a plasmid with "two Reps in one plasmid." Furthermore, we determined the copy number and stability of plasmids with cyanobacterial Reps in Synechococcus 7942 and Synechocystis 6803 to elucidate their potential applications. The novel properties of CyRepX1 and 2, as revealed by this study, hold promise for the development of innovative genetic engineering tools in cyanobacteria.

miRNAs in Heart Development and Disease

Cardiovascular diseases (CVD) are a group of disorders that affect the heart and blood vessels. They include conditions such as myocardial infarction, coronary artery disease, heart failure, arrhythmia, and congenital heart defects. CVDs are the leading cause of death worldwide. Therefore, new medical interventions that aim to prevent, treat, or manage CVDs are of prime importance. MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression at the posttranscriptional level and play important roles in various biological processes, including cardiac development, function, and disease. Moreover, miRNAs can also act as biomarkers and therapeutic targets. In order to identify and characterize miRNAs and their target genes, scientists take advantage of computational tools such as bioinformatic algorithms, which can also assist in analyzing miRNA expression profiles, functions, and interactions in different cardiac conditions. Indeed, the combination of miRNA research and bioinformatic algorithms has opened new avenues for understanding and treating CVDs. In this review, we summarize the current knowledge on the roles of miRNAs in cardiac development and CVDs, discuss the challenges and opportunities, and provide some examples of recent bioinformatics for miRNA research in cardiovascular biology and medicine.

Asgard archaeal selenoproteome reveals a roadmap for the archaea-to-eukaryote transition of selenocysteine incorporation machinery

Selenocysteine (Sec) is encoded by the UGA codon that normally functions as a stop signal and is specifically incorporated into selenoproteins via a unique recoding mechanism. The translational recoding of UGA as Sec is directed by an unusual RNA structure, the SECIS element. Although archaea and eukaryotes adopt similar Sec encoding machinery, the SECIS elements have no similarities to each other with regard to sequence and structure. We analyzed >400 Asgard archaeal genomes to examine the occurrence of both Sec encoding system and selenoproteins in this archaeal superphylum, the closest prokaryotic relatives of eukaryotes. A comprehensive map of Sec utilization trait has been generated, providing the most detailed understanding of the use of this nonstandard amino acid in Asgard archaea so far. By characterizing the selenoproteomes of all organisms, several selenoprotein-rich phyla and species were identified. Most Asgard archaeal selenoprotein genes possess eukaryotic SECIS-like structures with varying degrees of diversity. Moreover, euryarchaeal SECIS elements might originate from Asgard archaeal SECIS elements via lateral gene transfer, indicating a complex and dynamic scenario of the evolution of SECIS element within archaea. Finally, a roadmap for the transition of eukaryotic SECIS elements from archaea was proposed, and selenophosphate synthetase may serve as a potential intermediate for the generation of ancestral eukaryotic SECIS element. Our results offer new insights into a deeper understanding of the evolution of Sec insertion machinery.

Characterization of the complete mitochondrial genomes of the zoonotic parasites Bolbosoma nipponicum and Corynosoma villosum (Acanthocephala: Polymorphida) and the molecular phylogeny of the order Polymorphida

Acanthocephalans of the order Polymorphida mainly parasitic in birds and mammals, are of veterinary, medical and economic importance. However, the evolutionary relationships of its 3 families (Centrorhynchidae, Polymorphidae and Plagiorhynchidae) remain under debate. Additionally, some species of Polymorphida (i.e. Bolbosoma spp. and Corynosoma spp.) are recognized as zoonotic parasites, associated with human acanthocephaliasis, but the mitochondrial genomes for representatives of Bolbosoma and Corynosoma have not been reported so far. In the present study, the complete mitochondrial genomes B. nipponicum and C. villosum (Acanthocephala: Polymorphidae) are reported for the first time, which are 14 296 and 14 241 bp in length, respectively, and both contain 36 genes [including 12 PCGs, 22 tRNA genes and 2 rRNA genes] and 2 non-coding regions (NCR1 and NCR2). The gene arrangement of some tRNAs in the mitogenomes of B. nipponicum and C. villosum differs from that found in all other acanthocephalans, except Polymorphus minutus. Phylogenetic results based on concatenated amino acid (AA) sequences of the 12 protein-coding genes (PCGs) strongly supported that the family Polymorphidae is a sister to the Centrorhynchidae rather than the Plagiorhynchidae, and also confirmed the sister relationship of the genera Bolbosoma and Corynosoma in the Polymorphidae based on the mitogenomic data for the first time. Our present findings further clarified the phylogenetic relationships of the 3 families Plagiorhynchidae, Centrorhynchidae and Polymorphidae, enriched the mitogenome data of the phylum Acanthocephala (especially the order Polymorphida), and provided the resource of genetic data for diagnosing these 2 pathogenic parasites of human acanthocephaliasis.

Targeting non-coding RNAs: Perspectives and challenges of in-silico approaches

The growing information currently available on the central role of non-coding RNAs (ncRNAs) including microRNAs (miRNAS) and long non-coding RNAs (lncRNAs) for chronic and degenerative human diseases makes them attractive therapeutic targets. RNAs carry out different functional roles in human biology and are deeply deregulated in several diseases. So far, different attempts to therapeutically target the 3D RNA structures with small molecules have been reported. In this scenario, the development of computational tools suitable for describing RNA structures and their potential interactions with small molecules is gaining more and more interest. Here, we describe the most suitable strategies to study ncRNAs through computational tools. We focus on methods capable of predicting 2D and 3D ncRNA structures. Furthermore, we describe computational tools to identify, design and optimize small molecule ncRNA binders. This review aims to outline the state of the art and perspectives of computational methods for ncRNAs over the past decade.

GRIK1-AS1 deficiency accelerates endometriosis progression by boosting DNMT1-dependent SFRP1 promoter methylation in endometrial stromal cells

BACKGROUND

Endometriosis, a gynecological disease that affects up to 10% of women, is a major cause of pain and infertility. Deregulation of the epigenome is accountable for the onset and progression of endometriosis, although its exact mechanism is unknown. The purpose of the current study is to examine the role of the long non-coding RNA (lncRNA) GRIK1-AS1 in the epigenetic regulation of endometrial stromal cell proliferation and the development of endometriosis.

METHODS

Endometriosis datasets were screened to identify GRIKI-AS1 as dramatically declining in endometriosis. Gain or loss of function endometrial stromal cell (ESC) models were established. The anti-proliferation phenotype was investigated using in vitro and in vivo experiments. Epigenetic regulatory network analyses were conducted to suggest the intrinsic molecular mechanism.

RESULTS

With bioinformatic and clinical data, we observed that GRIK1-AS1 and SFRP1 were expressed at low levels in endometriosis. Overexpressed GRIK1-AS1 inhibited ESC proliferation, while SFRP1 knockdown rescued the antiproliferative ability of GRIK1-AS1. Specifically, methylation-dependent expression inhibition of SFRP1 was revealed in ESCs. Mechanistically, GRIK1-AS1 hampers the occupancy of DNMT1 in SRFP1 promoter, leading to hypomethylation of SFRP1 and upregulated SFRP1 expression, thereby potentially suppressing Wnt signaling and its adverse proliferative effect. Therapeutically, lentivirus-mediated upregulation of GRIK1-AS1 inhibited endometriosis disease progression in vivo.

CONCLUSIONS

Our study is a proof-of-concept demonstration for GRIKI-AS1-associated endometriosis pathogenesis and highlights a potential intervention target.

G-quadruplex forming sequences in the genes coding for cytochrome P450 enzymes and their potential roles in drug metabolism

The majority of drugs are metabolized by cytochrome P450 (CYP) enzymes, primarily belonging to the CYP1, CYP2 and CYP3 families. Genetic variations are the main cause of inter-individual differences in drug response, which constitutes a major concern in pharmacotherapy. G-quadruplexes (G4s), are non-canonical DNA and RNA secondary structures formed by guanine-rich sequences. G4s have been implicated in cancer and gene regulation. In this study, we investigated putative G4-forming sequences (PQSs) in the CYP genes. Our findings reveal a high density of PQSs in the full genes of CYP family 2. Moreover, we observe an increased density of PQSs in the promoters of CYP family 1 genes compared to non-CYP450 genes. Importantly, stable PQSs were also identified in all studied CYP genes. Subsequently, we assessed the impact of the most frequently reported genetic mutations in the selected genes and the possible effect of these mutations on G4 formation as well as on the thermodynamic stability of predicted G4s. We found that 4 SNPs overlap G4 sequences and lead to mutated DNA and RNA G4 forming sequences in their context. Notably, the mutation in the CYP2C9 gene, which is associated with impaired (S)-warfarin metabolism in patients, alters a G4 sequence. We then demonstrated that at least 10 of the 13 chosen cytochrome P450 G4 candidates form G-quadruplex structures in vitro, using a combination of spectroscopic methods. In conclusion, our findings indicate the potential role of G-quadruplexes in the regulation of cytochrome genes, and emphasize the importance of G-quadruplexes in drug metabolism.

First complete mitochondrial genome of the tribe Coccini (Hemiptera, Coccomorpha, Coccidae) and its phylogenetic implications

Soft scale insects (Hemiptera, Coccidae) are important pests of various agricultural and horticultural crops and ornamental plants. They have negative impacts on agriculture and forestry. The tribe Coccini represents one of the most ancient evolutionary lineages of soft scale insects. However, no complete Coccini mitochondrial genome (mitogenome) is available in public databases. Here, we described the complete mitogenome of Coccushesperidum L., 1758. The 15,566 bp mitogenome of C.hesperidum had a high A+T content (83.4%) and contained a typical set of 37 genes, with 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs) and two ribosomal RNA genes (rRNAs). Only seven tRNAs had the typical clover-leaf secondary structure and the remaining tRNAs lacked the DHU arm, TψC arm or both. Moreover, a comparative analysis of all reported scale insect mitogenomes from GenBank database was performed. The mitogenomes of scale insects showed high similarities in base composition and A+T content. Additionally, our phylogenetic analysis confirmed the monophyly of Coccomorpha and revealed that the archaeococcoids were the most basal lineage within Coccomorpha, while Ericeruspela and Didesmococcuskoreanus, belonging to Coccidae, were often mixed with Aclerdidae, making Coccidae a paraphyletic group. These findings expand the mitogenome database of scale insects and provide new insights on mitogenome evolution for future studies across different insect groups.

Evidence that an Unnamed Isometric Virus Associated with Potato Rugose Disease in Peru Is a New Species of Genus Torradovirus

A previously uncharacterized torradovirus species infecting potatoes was detected by high-throughput sequencing from field samples from Peru and in customs intercepts in potato tubers that originated from South America in the United States of America and the Netherlands. This new potato torradovirus showed high nucleotide sequence identity to an unidentified isometric virus (SB26/29), which was associated with a disease named potato rugose stunting in southern Peru characterized over two decades ago. Thus, this virus is tentatively named potato rugose stunting virus (PotRSV). The genome of PotRSV isolates sequenced in this study were composed of two polyadenylated RNA segments. RNA1 ranges from 7,086 to 7,089 nt and RNA2 from 5,228 to 5,230 nt. RNA1 encodes a polyprotein containing the replication block (helicase-protease-polymerase), whereas RNA2 encodes a polyprotein cleaved into a movement protein and the three capsid proteins (CPs). Pairwise comparison among PotRSV isolates revealed amino acid identity values greater than 86% in the protease-polymerase (Pro-Pol) region and greater than 82% for the combined CPs. The closest torradovirus species, squash chlorotic leaf spot virus, shares amino acid identities of ∼58 and ∼41% in the Pro-Pol and the combined CPs, respectively. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Analysis of G-quadruplex forming sequences in podocytes-marker genes and their potential roles in inherited glomerular diseases

Nephrotic Syndrome is the most widespread pediatric kidney disorder. Genetic alterations in podocyte genes are thought to be responsible for the disease. G-quadruplexes are non-conventional guanine-rich DNA and RNA structures, which are commonly found in regulatory regions. This study examined the potential G-quadruplexes forming sequences in the promoters and gene bodies of podocyte-marker genes. High G-quadruplexes density was found in the vascular endothelial growth facto, cluster of differentiation-151, integrin subunit beta-4, metalloendopeptidase, Wilms tumor-1, integrin subunit beta-3, synaptopodin, and nephrin promoters. Vascular endothelial growth facto, cluster of differentiation-151 and integrin subunit beta-4 had the highest G-quadruplexes density in their gene bodies and promoters. Additionally, highly stable G-quadruplexes forming sequences were identified within all podocyte-marker genes. Furthermore, it is hypothesized that Wilms tumor-1 is capable of controlling the transcription of podocalyxin by binding to two possible G-quadruplexes forming motifs. We next analyzed the most frequently reported genetic mutations in the selected genes for their effect on DNA G-quadruplexes formation, and the thermodynamic stability of predicted RNA G-quadruplexes, using RNAfold. Importantly, the missense mutation c.121_122del in the nephrin gene reported in patients with NS type 1 affected DNA G-quadruplexes formation in this region as well as the thermodynamic stability of the corresponding RNA G-quadruplexes. Overall, we report the potential regulatory roles of G-quadruplexes in the etiology of nephrotic syndrome and their possible use as drug targets to treat kidney diseases.

Genomic characterization of dengue virus serotype 2 during dengue outbreak and endemics in Hangzhou, Zhejiang (2017-2019)

Introduction

Dengue fever (DF) is a mosquito-borne viral disease caused by the dengue virus (DENV). In recent years, Hangzhou has undergone a DF epidemic, particularly in 2017, with an outbreak of 1,128 patients. The study aimed to investigate the genetic diversity and molecular evolution among the DF clinical isolates during and after the outbreak to aid in mapping its spread.

Methods

To understand the genetic diversity, 74 DENV-2 strains were isolated from DF epidemic cases between 2017 and 2019. Combining whole genome sequencing (WGS) technology, additional phylogenetic, haplotype, amino acid (AA) substitution, and recombination analyses were performed.

Results

The results revealed that strains from 2017 were closely related to those from Singapore, Malaysia, and Thailand, indicating an imported international transmission. Local strains from 2018 were clustered with those recovered from 2019 and were closely associated with Guangzhou isolates, suggesting a within-country transmission after the significant outbreak in 2017. Compared to DENV-2 virus P14337 (Thailand/0168/1979), a total of 20 AA substitutions were detected. Notably, V431I, T2881I, and K3291T mutations only occurred in indigenous cases from 2017, and A1402T, V1457I, Q2777E, R3189K, and Q3310R mutations were exclusively found in imported cases from 2018 to 2019. The recombination analysis indicated that a total of 14 recombination events were observed.

Conclusion

This study may improve our understanding of DENV transmission in Hangzhou and provide further insight into DENV-2 transmission and the local vaccine choice.

Genome-wide identification and characterization of Puccinia striiformis-responsive lncRNAs in Triticum aestivum

Wheat stripe rust (yellow rust) caused by Puccinia striiformis f. sp. tritici (Pst) is a serious biotic stress factor limiting wheat production worldwide. Emerging evidence demonstrates that long non-coding RNAs (lncRNAs) participate in various developmental processes in plants via post-transcription regulation. In this study, RNA sequencing (RNA-seq) was performed on a pair of near-isogenic lines-rust resistance line FLW29 and rust susceptible line PBW343-which differed only in the rust susceptibility trait. A total of 6,807 lncRNA transcripts were identified using bioinformatics analyses, among which 10 lncRNAs were found to be differentially expressed between resistance and susceptible lines. In order to find the target genes of the identified lncRNAs, their interactions with wheat microRNA (miRNAs) were predicted. A total of 199 lncRNAs showed interactions with 65 miRNAs, which further target 757 distinct mRNA transcripts. Moreover, detailed functional annotations of the target genes were used to identify the candidate genes, pathways, domains, families, and transcription factors that may be related to stripe rust resistance response in wheat plants. The NAC domain protein, disease resistance proteins RPP13 and RPM1, At1g58400, monodehydroascorbate reductase, NBS-LRR-like protein, rust resistance kinase Lr10-like, LRR receptor, serine/threonine-protein kinase, and cysteine proteinase are among the identified targets that are crucial for wheat stripe rust resistance. Semiquantitative PCR analysis of some of the differentially expressed lncRNAs revealed variations in expression profiles of two lncRNAs between the Pst-resistant and Pst-susceptible genotypes at least under one condition. Additionally, simple sequence repeats (SSRs) were also identified from wheat lncRNA sequences, which may be very useful for conducting targeted gene mapping studies of stripe rust resistance in wheat. These findings improved our understanding of the molecular mechanism responsible for the stripe rust disease that can be further utilized to develop wheat varieties with durable resistance to this disease.

Dynamic Transcriptional Landscape of Mycobacterium smegmatis under Cold Stress

Bacterial adaptation to cold stress requires wide transcriptional reprogramming. However, the knowledge of molecular mechanisms underlying the cold stress response of mycobacteria is limited. We conducted comparative transcriptomic analysis of Mycobacterium smegmatis subjected to cold shock. The growth of M. smegmatis cultivated at 37 °C was arrested just after exposure to cold (acclimation phase) but later (by 24 h) was resumed at a much slower rate (adaptation phase). Transcriptomic analyses revealed distinct gene expression patterns corresponding to the two phases. During the acclimation phase, differential expression was observed for genes associated with cell wall remodeling, starvation response, and osmotic pressure stress, in parallel with global changes in the expression of transcription factors and the downregulation of ribosomal genes, suggesting an energy-saving strategy to support survival. At the adaptation phase, the expression profiles were recovered, indicating restoration of the processes repressed earlier. Comparison of transcriptional responses in M. smegmatis with those in other bacteria revealed unique adaptation strategies developed by mycobacteria. Our findings shed light on the molecular mechanisms underlying M. smegmatis survival under cold stress. Further research should clarify whether the discovered transcriptional mechanisms exist in other mycobacterial species, including pathogenic Mycobacterium tuberculosis, which could be important for transmission control.

Mitogenomic phylogenies suggest the resurrection of the subfamily Porrocaecinae and provide insights into the systematics of the superfamily Ascaridoidea (Nematoda: Ascaridomorpha), with the description of a new species of Porrocaecum

BACKGROUND

The family Toxocaridae is a group of zooparasitic nematodes of veterinary, medical and economic significance. However, the evolutionary relationship of Porrocaecum and Toxocara, both genera currently classified in Toxocaridae, and the monophyly of the Toxocaridae remain under debate. Moreover, the validity of the subgenus Laymanicaecum in the genus Porrocaecum is open to question. Due to the scarcity of an available genetic database, molecular identification of Porrocaecum nematodes is still in its infancy.

METHODS

A number of Porrocaecum nematodes collected from the Eurasian marsh harrier Circus aeruginosus (Linnaeus) (Falconiformes: Accipitridae) in the Czech Republic were identified using integrated morphological methods (light and scanning electron microscopy) and molecular techniques (sequencing and analyzing the nuclear 18S, 28S and ITS regions). The complete mitochondrial genomes of the collected nematode specimens and of Porrocaecum (Laymanicaecum) reticulatum (Linstow, 1899) were sequenced and annotated for the first time. Phylogenetic analyses of ascaridoid nematodes based on the amino acid sequences of 12 protein-coding genes of mitochondrial genomes were performed using maximum likelihood and Bayesian inference.

RESULTS

A new species of Porrocaecum, named P. moraveci n. sp., is described based on the morphological and genetic evidence. The mitogenomes of P. moraveci n. sp. and P. reticulatum both contain 36 genes and are 14,517 and 14,210 bp in length, respectively. Comparative mitogenomics revealed that P. moraveci n. sp. represents the first known species with three non-coding regions and that P. reticulatum has the lowest overall A + T content in the mitogenomes of ascaridoid nematodes tested to date. Phylogenetic analyses showed the representatives of Toxocara clustered together with species of the family Ascarididae rather than with Porrocaecum and that P. moraveci n. sp. is a sister to P. reticulatum.

CONCLUSIONS

The characterization of the complete mitochondrial genomes of P. moraveci n. sp. and P. reticulatum is reported for the first time. Mitogenomic phylogeny analyses indicated that the family Toxocaridae is non-monophyletic and that the genera Porrocaecum and Toxocara do not have an affinity. The validity of the subgenus Laymanicaecum in Porrocaecum was also rejected. Our results suggest that: (i) Toxocaridae should be degraded to a subfamily of the Ascarididae that includes only the genus Toxocara; and (ii) the subfamily Porrocaecinae should be resurrected to include only the genus Porrocaecum. The present study enriches the database of ascaridoid mitogenomes and provides a new insight into the systematics of the superfamily Ascaridoidea.

Variable Region Sequences Influence 16S rRNA Performance

16S rRNA gene sequences are commonly analyzed for taxonomic and phylogenetic studies because they contain variable regions that can help distinguish different genera. However, intra-genus distinction using variable region homology is often impossible due to the high overall sequence identities among closely related species, even though some residues may be conserved within respective species. Using a computational method that included the allelic diversity within individual genomes, we discovered that certain Escherichia and Shigella species can be distinguished by a multi-allelic 16S rRNA variable region single nucleotide polymorphism (SNP). To evaluate the performance of 16S rRNAs with altered variable regions, we developed an in vivo system that measures the acceptance and distribution of variant 16S rRNAs into a large pool of natural versions supporting normal translation and growth. We found that 16S rRNAs containing evolutionarily disparate variable regions were underpopulated both in ribosomes and in active translation pools, even for an SNP. Overall, this study revealed that variable region sequences can substantially influence the performance of 16S rRNAs and that this biological constraint can be leveraged to justify refining taxonomic assignments of variable region sequence data. IMPORTANCE This study reevaluates the notion that 16S rRNA gene variable region sequences are uninformative for intra-genus classification and that single nucleotide variations within them have no consequence to strains that bear them. We demonstrated that the performance of 16S rRNAs in Escherichia coli can be negatively impacted by sequence changes in variable regions, even for single nucleotide changes that are native to closely related Escherichia and Shigella species; thus, biological performance is likely constraining the evolution of variable regions in bacteria. Further, the native nucleotide variations we tested occur in all strains of their respective species and across their multiple 16S rRNA gene copies, suggesting that these species evolved beyond what would be discerned from a consensus sequence comparison. Therefore, this work also reveals that the multiple 16S rRNA gene alleles found in most bacteria can provide more informative phylogenetic and taxonomic detail than a single reference allele.

The complete mitochondrial genomes of the Leuciscus baicalensis and Rutilus rutilus: a detailed genomic comparison among closely related species of the Leuciscinae subfamily

Cyprinidae is the largest family in the order of freshwater fish Cypriniformes. Increased subfamily members of Cyprinidae have been suggested to be re-classified for decades. In this study, we sequenced the mitochondrial genomes (mitogenomes) of Leuciscus baicalensis and Rutilus rutilus collected from northwest China and compared with other closely related species to determine their associated family or subfamily. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of Leuciscus baicalensis and Rutilus rutilus and characterized the mitogenomes by the gene structure, gene order, and the secondary structures of the 22 tRNA genes. We compared mitogenome features of Leuciscinae with other subfamilies in Cyprinidae. We used the analytic Bayesian Information and Maximum Likelihood methods to determine phylogenetic trees of 13 PCGs. The mitogenomes of Leuciscus baicalensis and Rutilus rutilus were 16,607 bp and 16,606 bp, respectively. Organization and location of these genes were consistent with already studied Leuciscinae fishes. Synonymous codon usage was conservative in Leuciscinae as compared with other subfamilies in Cyprinidae. Phylogenetic analysis indicated that Leuciscinae was a monophyletic group, and genus Leuciscus was a paraphyletic group. Our approach, for the first time, of studying comparative mitochondrial genomics and phylogenetics together provided a supportive platform to the analysis of population genetics and phylogeny for Leuciscinae. Our results indicated a promising potential of comparative mitochondrial genomics in the manifestation of phylogenetic relationships between fishes, leading us to a suggestion that mitogenomes should be routinely considered in clarifying phylogenetics of family and subfamily members of fish.

Thirteen dubious ways to detect conserved structural RNAs

Covariation induced by compensatory base substitutions in RNA alignments is a great way to deduce conserved RNA structure, in principle. In practice, success depends on many factors, importantly the quality and depth of the alignment and the choice of covariation statistic. Measuring covariation between pairs of aligned positions is easy. However, using covariation to infer evolutionarily conserved RNA structure is complicated by other extraneous sources of covariation such as that resulting from homologous sequences having evolved from a common ancestor. In order to provide evidence of evolutionarily conserved RNA structure, a method to distinguish covariation due to sources other than RNA structure is necessary. Moreover, there are several sorts of artifactually generated covariation signals that can further confound the analysis. Additionally, some covariation signal is difficult to detect due to incomplete comparative data. Here, we investigate and critically discuss the practice of inferring conserved RNA structure by comparative sequence analysis. We provide new methods on how to approach and decide which of the numerous long non-coding RNAs (lncRNAs) have biologically relevant structures.

Actively transcribed rDNA and distal junction (DJ) sequence are involved in association of NORs with nucleoli

Although they are organelles without a limiting membrane, nucleoli have an exclusive structure, built upon the rDNA-rich acrocentric short arms of five human chromosomes (nucleolar organizer regions or NORs). This has raised the question: what are the structural features of a chromosome required for its inclusion in a nucleolus? Previous work has suggested that sequences adjacent to the tandemly repeated rDNA repeat units (DJ, distal junction sequence) may be involved, and we have extended such studies by addressing several issues related to the requirements for the association of NORs with nucleoli. We exploited both a set of somatic cell hybrids containing individual human acrocentric chromosomes and a set of Human Artificial Chromosomes (HACs) carrying different parts of a NOR, including an rDNA unit or DJ or PJ (proximal junction) sequence. Association of NORs with nucleoli was increased when constituent rDNA was transcribed and may be also affected by the status of heterochromatin blocks formed next to the rDNA arrays. Furthermore, our data suggest that a relatively small size DJ region, highly conserved in evolution, is also involved, along with the rDNA repeats, in the localization of p-arms of acrocentric chromosomes in nucleoli. Thus, we infer a cooperative action of rDNA sequence-stimulated by its activity-and sequences distal to rDNA contributing to incorporation into nucleoli. Analysis of NOR sequences also identified LncRNA_038958 in the DJ, a candidate transcript with the region of the suggested promoter that is located close to the DJ/rDNA boundary and contains CTCF binding sites. This LncRNA may affect RNA Polymerase I and/or nucleolar activity. Our findings provide the basis for future studies to determine which RNAs and proteins interact critically with NOR sequences to organize the higher-order structure of nucleoli and their function in normal cells and pathological states.

Phylomitogenomic Analyses Provided Further Evidence for the Resurrection of the Family Pseudoacanthocephalidae (Acanthocephala: Echinorhynchida)

The phylum Acanthocephala is an important monophyletic group of parasites, with adults parasitic in the digestive tracts of all major vertebrate groups. Acanthocephalans are of veterinary, medical, and economic importance due to their ability to cause disease in domestic animals, wildlife, and humans. However, the current genetic data for acanthocephalans are sparse, both in terms of the proportion of taxa surveyed and the number of genes sequenced. Consequently, the basic molecular phylogenetic framework for the phylum is still incomplete. In the present study, we reported the first complete mitochondrial genome from a representative of the family Pseudoacanthocephalidae Petrochenko, 1956. The mitogenome of Pseudoacanthocephalus bufonis (Shipley, 1903) is 14,056 bp in length, contains 36 genes (12 protein-coding genes (PCGs) (lacking atp8), 22 tRNA genes, and 2 rRNA genes (rrnL and rrnS)) and two non-coding regions (NCR1 and NCR2), and displayed the highest GC-skew in the order Echinorhynchida. Phylogenetic results of maximum likelihood (ML) and Bayesian inference (BI) using the amino acid sequences of 12 protein-coding genes in different models provided further evidence for the resurrection of the family Pseudoacanthocephalidae and also supported that the order Echinorhynchida is paraphyletic. A monophyletic clade comprising P. bufonis and Cavisoma magnum suggests a close affinity between Pseudoacanthocephalidae and Cavisomatidae. Our phylogenetic analyses also showed that Polymorphidae has a closer relationship with Centrorhynchidae than Plagiorhynchidae in the monophyletic order Polymorphida.

A hybrid approach to assess the structural impact of long noncoding RNA mutations uncovers key NEAT1 interactions in colorectal cancer

Long noncoding RNAs (lncRNAs) are emerging players in cancer and they entail potential as prognostic biomarkers or therapeutic targets. Earlier studies have identified somatic mutations in lncRNAs that are associated with tumor relapse after therapy, but the underlying mechanisms behind these associations remain unknown. Given the relevance of secondary structure for the function of some lncRNAs, some of these mutations may have a functional impact through structural disturbance. Here, we examined the potential structural and functional impact of a novel A > G point mutation in NEAT1 that has been recurrently observed in tumors of colorectal cancer patients experiencing relapse after treatment. Here, we used the nextPARS structural probing approach to provide first empirical evidence that this mutation alters NEAT1 structure. We further evaluated the potential effects of this structural alteration using computational tools and found that this mutation likely alters the binding propensities of several NEAT1-interacting miRNAs. Differential expression analysis on these miRNA networks shows upregulation of Vimentin, consistent with previous findings. We propose a hybrid pipeline that can be used to explore the potential functional effects of lncRNA somatic mutations.

Regulation of RNase E during the UV stress response in the cyanobacterium Synechocystis sp. PCC 6803

Endoribonucleases govern the maturation and degradation of RNA and are indispensable in the posttranscriptional regulation of gene expression. A key endoribonuclease in Gram-negative bacteria is RNase E. To ensure an appropriate supply of RNase E, some bacteria, such as Escherichia coli, feedback-regulate RNase E expression via the rne 5'-untranslated region (5' UTR) in cis. However, the mechanisms involved in the control of RNase E in other bacteria largely remain unknown. Cyanobacteria rely on solar light as an energy source for photosynthesis, despite the inherent ultraviolet (UV) irradiation. In this study, we first investigated globally the changes in gene expression in the cyanobacterium Synechocystis sp. PCC 6803 after a brief exposure to UV. Among the 407 responding genes 2 h after UV exposure was a prominent upregulation of rne mRNA level. Moreover, the enzymatic activity of RNase E rapidly increased as well, although the protein stability decreased. This unique response was underpinned by the increased accumulation of full-length rne mRNA caused by the stabilization of its 5' UTR and suppression of premature transcriptional termination, but not by an increased transcription rate. Mapping of RNA 3' ends and in vitro cleavage assays revealed that RNase E cleaves within a stretch of six consecutive uridine residues within the rne 5' UTR, indicating autoregulation. These observations suggest that RNase E in cyanobacteria contributes to reshaping the transcriptome during the UV stress response and that its required activity level is secured at the RNA level despite the enhanced turnover of the protein.

Regulation of pSYSA defense plasmid copy number in Synechocystis through RNase E and a highly transcribed asRNA

Synthetic biology approaches toward the development of cyanobacterial producer strains require the availability of appropriate sets of plasmid vectors. A factor for the industrial usefulness of such strains is their robustness against pathogens, such as bacteriophages infecting cyanobacteria. Therefore, it is of great interest to understand the native plasmid replication systems and the CRISPR-Cas based defense mechanisms already present in cyanobacteria. In the model cyanobacterium Synechocystis sp. PCC 6803, four large and three smaller plasmids exist. The ~100 kb plasmid pSYSA is specialized in defense functions by encoding all three CRISPR-Cas systems and several toxin-antitoxin systems. The expression of genes located on pSYSA depends on the plasmid copy number in the cell. The pSYSA copy number is positively correlated with the expression level of the endoribonuclease E. As molecular basis for this correlation we identified the RNase E-mediated cleavage within the pSYSA-encoded ssr7036 transcript. Together with a cis-encoded abundant antisense RNA (asRNA1), this mechanism resembles the control of ColE1-type plasmid replication by two overlapping RNAs, RNA I and II. In the ColE1 mechanism, two non-coding RNAs interact, supported by the small protein Rop, which is encoded separately. In contrast, in pSYSA the similar-sized protein Ssr7036 is encoded within one of the interacting RNAs and it is this mRNA that likely primes pSYSA replication. Essential for plasmid replication is furthermore the downstream encoded protein Slr7037 featuring primase and helicase domains. Deletion of slr7037 led to the integration of pSYSA into the chromosome or the other large plasmid pSYSX. Moreover, the presence of slr7037 was required for successful replication of a pSYSA-derived vector in another model cyanobacterium, Synechococcus elongatus PCC 7942. Therefore, we annotated the protein encoded by slr7037 as Cyanobacterial Rep protein A1 (CyRepA1). Our findings open new perspectives on the development of shuttle vectors for genetic engineering of cyanobacteria and of modulating the activity of the entire CRISPR-Cas apparatus in Synechocystis sp. PCC 6803.

Potential of siRNA in COVID-19 therapy: Emphasis on in silico design and nanoparticles based delivery

Small interfering RNA (siRNA)-mediated mRNA degradation approach have imparted its eminence against several difficult-to-treat genetic disorders and other allied diseases. Viral outbreaks and resulting pandemics have repeatedly threatened public health and questioned human preparedness at the forefront of drug design and biomedical readiness. During the recent pandemic caused by the SARS-CoV-2, mRNA-based vaccination strategies have paved the way for a new era of RNA therapeutics. RNA Interference (RNAi) based approach using small interfering RNA may complement clinical management of the COVID-19. RNA Interference approach will primarily work by restricting the synthesis of the proteins required for viral replication, thereby hampering viral cellular entry and trafficking by targeting host as well as protein factors. Despite promising benefits, the stability of small interfering RNA in the physiological environment is of grave concern as well as site-directed targeted delivery and evasion of the immune system require immediate attention. In this regard, nanotechnology offers viable solutions for these challenges. The review highlights the potential of small interfering RNAs targeted toward specific regions of the viral genome and the features of nanoformulations necessary for the entrapment and delivery of small interfering RNAs. In silico design of small interfering RNA for different variants of SARS-CoV-2 has been discussed. Various nanoparticles as promising carriers of small interfering RNAs along with their salient properties, including surface functionalization, are summarized. This review will help tackle the real-world challenges encountered by the in vivo delivery of small interfering RNAs, ensuring a safe, stable, and readily available drug candidate for efficient management of SARS-CoV-2 in the future.

The identification of conserved sequence features of co-translationally decayed mRNAs and upstream open reading frames in angiosperm transcriptomes

RNA turnover is essential in maintaining messenger RNA (mRNA) homeostasis during various developmental stages and stress responses. Co-translational mRNA decay (CTRD), a process in which mRNAs are degraded while still associated with translating ribosomes, has recently been discovered to function in yeast and three angiosperm transcriptomes. However, it is still unclear how prevalent CTRD across the plant lineage. Moreover, the sequence features of co-translationally decayed mRNAs have not been well-studied. Here, utilizing a collection of publicly available degradome sequencing datasets for another seven angiosperm transcriptomes, we have confirmed that CTRD is functioning in at least 10 angiosperms and likely throughout the plant lineage. Additionally, we have identified sequence features shared by the co-translationally decayed mRNAs in these species, implying a possible conserved triggering mechanism for this pathway. Given that degradome sequencing datasets can also be used to identify actively translating upstream open reading frames (uORFs), which are quite understudied in plants, we have identified numerous actively translating uORFs in the same 10 angiosperms. These findings reveal that actively translating uORFs are prevalent in plant transcriptomes, some of which are conserved across this lineage. We have also observed conserved sequence features in the regions flanking these uORFs' stop codons that might contribute to ribosome stalling at these sequences. Finally, we discovered that there were very few overlaps between the mRNAs harboring actively translating uORFs and those sorted into the co-translational decay pathway in the majority of the studied angiosperms, suggesting that these two processes might be nearly mutually exclusive in those species. In total, our findings provide the identification of CTRD and actively translating uORFs across a broad collection of plants and provide novel insights into the important sequence features associated with these collections of mRNAs and regulatory elements, respectively.

Thermodynamic and structural characterization of an EBV infected B-cell lymphoma transcriptome

Epstein-Barr virus (EBV) is a widely prevalent human herpes virus infecting over 95% of all adults and is associated with a variety of B-cell cancers and induction of multiple sclerosis. EBV accomplishes this in part by expression of coding and noncoding RNAs and alteration of the host cell transcriptome. To better understand the structures which are forming in the viral and host transcriptomes of infected cells, the RNA structure probing technique Structure-seq2 was applied to the BJAB-B1 cell line (an EBV infected B-cell lymphoma). This resulted in reactivity profiles and secondary structural analyses for over 10000 human mRNAs and lncRNAs, along with 19 lytic and latent EBV transcripts. We report in-depth structural analyses for the human MYC mRNA and the human lncRNA CYTOR. Additionally, we provide a new model for the EBV noncoding RNA EBER2 and provide the first reported model for the EBV tandem terminal repeat RNA. In-depth thermodynamic and structural analyses were carried out with the motif discovery tool ScanFold and RNAfold prediction tool; subsequent covariation analyses were performed on resulting models finding various levels of support. ScanFold results for all analyzed transcripts are made available for viewing and download on the user-friendly RNAStructuromeDB.

Oropouche virus as an emerging cause of acute febrile illness in Colombia

Arbovirus infections are frequent causes of acute febrile illness (AFI) in tropical countries. We conducted health facility-based AFI surveillance at four sites in Colombia (Cucuta, Cali, Villavicencio, Leticia) during 2019-2022. Demographic, clinical and risk factor data were collected from persons with AFI that consented to participate in the study (n = 2,967). Serologic specimens were obtained and tested for multiple pathogens by RT-PCR and rapid test (Antigen/IgM), with 20.7% identified as dengue positive from combined testing. Oropouche virus (OROV) was initially detected in serum by metagenomic next-generation sequencing (mNGS) and virus target capture in a patient from Cúcuta. Three additional infections from Leticia were confirmed by conventional PCR, sequenced, and isolated in tissue culture. Phylogenetic analysis determined there have been at least two independent OROV introductions into Colombia. To assess OROV spread, a RT-qPCR dual-target assay was developed which identified 87/791 (10.9%) viremic cases in AFI specimens from Cali (3/53), Cucuta (3/19), Villavicencio (38/566), and Leticia (43/153). In parallel, an automated anti-nucleocapsid antibody assay detected IgM in 27/503 (5.4%) and IgG in 92/568 (16.2%) patients screened, for which 24/68 (35.3%) of PCR positives had antibodies. Dengue was found primarily in people aged <18 years and linked to several clinical manifestations (weakness, skin rash and petechiae), whereas Oropouche cases were associated with the location, climate phase, and odynophagia symptom. Our results confirm OROV as an emerging pathogen and recommend increased surveillance to determine its burden as a cause of AFI in Colombia.

5′ Untranslated mRNA Regions Allow Bypass of Host Cell Translation Inhibition by Legionella pneumophila

Legionella pneumophila grows within membrane-bound vacuoles in alveolar macrophages during human disease. Pathogen manipulation of the host cell is driven by bacterial proteins translocated through a type IV secretion system (T4SS).

Self-amplifying mRNA vaccines: Mode of action, design, development and optimization

The mRNA-based vaccines are quality-by-design (QbD) immunotherapies that provide safe, tunable, scalable, streamlined and potent treatment possibilities against different types of diseases. The self-amplifying mRNA (saRNA) vaccines, as a highly advantageous class of mRNA vaccines, are inspired by the intracellular self-multiplication nature of some positive-sense RNA viruses. Such vaccine platforms provide a relatively increased expression level of vaccine antigen(s) together with self-adjuvanticity properties. Lined with the QbD saRNA vaccines, essential optimizations improve the stability, safety, and immunogenicity of the vaccine constructs. Here, we elaborate on the concepts and mode-of-action of mRNA and saRNA vaccines, articulate the potential limitations or technical bottlenecks, and explain possible solutions or optimization methods in the process of their design and development.

Genome-Wide Screening for Pathogenic Proteins and microRNAs Associated with Parasite-Host Interactions in Trypanosoma brucei

Tsetse flies are a type of blood-sucking insect living in diverse locations in sub-Saharan Africa. These insects can transmit the unicellular parasite Trypanosoma brucei (T. brucei) which causes African trypanosomiasis in mammals. There remain huge unmet needs for prevention, early detection, and effective treatments for this disease. Currently, few studies have investigated the molecular mechanisms of parasite-host interactions underlying African trypanosomiasis, mainly due to a lack of understanding of the T. brucei genome. In this study, we dissected the genomic and transcriptomic profiles of T. brucei by annotating the genome and analyzing the gene expression. We found about 5% of T. brucei proteins in the human proteome, while more than 80% of T. brucei protein in other trypanosomes. Sequence alignment analysis showed that 142 protein homologs were shared among T. brucei and mammalian genomes. We identified several novel proteins with pathogenic potential supported by their molecular functions in T. brucei, including 24 RNA-binding proteins and six variant surface glycoproteins. In addition, 26 novel microRNAs were characterized, among which five miRNAs were not found in the mammalian genomes. Topology analysis of the miRNA-gene network revealed three genes (RPS27A, UBA52 and GAPDH) involved in the regulation of critical pathways related to the development of African trypanosomiasis. In conclusion, our work opens a new door to understanding the parasite-host interaction mechanisms by resolving the genome and transcriptome of T. brucei.

Transient RNA structures cause aberrant influenza virus replication and innate immune activation

During infection, the influenza A virus RNA polymerase produces both full-length and aberrant RNA molecules, such as defective viral genomes (DVGs) and mini viral RNAs (mvRNAs). Subsequent innate immune activation involves the binding of host pathogen receptor retinoic acid-inducible gene I (RIG-I) to viral RNAs. However, it is not clear what factors determine which influenza A virus RNAs are RIG-I agonists. Here, we provide evidence that RNA structures, called template loops (t-loops), stall the viral RNA polymerase and contribute to innate immune activation by mvRNAs during influenza A virus infection. Impairment of replication by t-loops depends on the formation of an RNA duplex near the template entry and exit channels of the RNA polymerase, and this effect is enhanced by mutation of the template exit path from the RNA polymerase active site. Overall, these findings are suggestive of a mechanism involving polymerase stalling that links aberrant viral replication to the activation of the innate immune response.

Immunomodulatory LncRNA on antisense strand of ICAM-1 augments SARS-CoV-2 infection-associated airway mucoinflammatory phenotype

Noncoding RNAs are important regulators of mucoinflammatory response, but little is known about the contribution of airway long noncoding RNAs (lncRNAs) in COVID-19. RNA-seq analysis showed a more than 4-fold increased expression of IL-6, ICAM-1, CXCL-8, and SCGB1A1 inflammatory factors; MUC5AC and MUC5B mucins; and SPDEF, FOXA3, and FOXJ1 transcription factors in COVID-19 patient nasal samples compared with uninfected controls. A lncRNA on antisense strand to ICAM-1 or LASI was induced 2-fold in COVID-19 patients, and its expression was directly correlated with viral loads. A SARS-CoV-2-infected 3D-airway model largely recapitulated these clinical findings. RNA microscopy and molecular modeling indicated a possible interaction between viral RNA and LASI lncRNA. Notably, blocking LASI lncRNA reduced the SARS-CoV-2 replication and suppressed MUC5AC mucin levels and associated inflammation, and select LASI-dependent miRNAs (e.g., let-7b-5p and miR-200a-5p) were implicated. Thus, LASI lncRNA represents an essential facilitator of SARS-CoV-2 infection and associated airway mucoinflammatory response.

Small RNAs in the Antarctic bacterium Pseudomonas extremaustralis responsive to oxygen availability and oxidative stress

Bacterial small non-coding RNAs (sRNAs) play key roles as genetic regulators, mediating in the adaptability to changing environmental conditions and stress responses. In this work, we analysed putative sRNAs identified by RNA-seq experiments in different aeration conditions in the extremophile bacterium P. extremaustralis. These analyses allowed the identification of 177 putative sRNAs under aerobiosis (A), microaerobiosis (M) and microaerobiosis after H₂ O₂ exposure (m-OS). The size and transcription profile of eight sRNAs with differential expression were verified by Northern blot. sRNA40, with unknown function but conserved in other Pseudomonas species, was selected to perform overexpression experiments followed by RNA-seq analysis. The overexpression of sRNA40 in P. extremaustralis resulted in significant expression changes of 19 genes with 14 differentially upregulated and five downregulated. Among the upregulated genes, eight transcripts corresponded to components of secretion systems, such as gspH, gspK, and gspM, belonging to the Type II secretion system, and rspO and rspP from Type III secretion system. Our results showed a novel sRNA which expression was triggered by low oxygen levels, and whose overexpression was associated with upregulation of selected components of protein secretion systems.

Implication of genetic variants in primary microRNA processing sites in the risk of multiple sclerosis

Background

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system with a well-established genetic contribution to susceptibility. Over 200 genetic regions have been linked to the inherited risk of developing MS, but the disease-causing variants and their functional effects at the molecular level are still largely unresolved. We hypothesised that MS-associated single-nucleotide polymorphisms (SNPs) affect the recognition and enzymatic cleavage of primary microRNAs (pri-miRNAs).

Methods

Our study focused on 11 pri-miRNAs (9 primate-specific) that are encoded in genetic risk loci for MS. The levels of mature miRNAs and potential isoforms (isomiRs) produced from those pri-miRNAs were measured in B cells obtained from the peripheral blood of 63 MS patients and 28 healthy controls. We tested for associations between SNP genotypes and miRNA expression in cis using quantitative trait locus (cis-miR-eQTL) analyses. Genetic effects on miRNA stem-loop processing efficiency were verified using luciferase reporter assays. Potential direct miRNA target genes were identified by transcriptome profiling and computational binding site assessment.

Findings

Mature miRNAs and isomiRs from hsa-mir-26a-2, hsa-mir-199a-1, hsa-mir-4304, hsa-mir-4423, hsa-mir-4464 and hsa-mir-4492 could be detected in all B-cell samples. When MS patient subgroups were compared with healthy controls, a significant differential expression was observed for miRNAs from the 5’ and 3’ strands of hsa-mir-26a-2 and hsa-mir-199a-1. The cis-miR-eQTL analyses and reporter assays pointed to a slightly more efficient Drosha-mediated processing of hsa-mir-199a-1 when the MS risk allele T of SNP rs1005039 is present. On the other hand, the MS risk allele A of SNP rs817478, which substitutes the first C in a CNNC sequence motif, was found to cause a markedly lower efficiency in the processing of hsa-mir-4423. Overexpression of hsa-mir-199a-1 inhibited the expression of 60 protein-coding genes, including IRAK2, MIF, TNFRSF12A and TRAF1. The only target gene identified for hsa-mir-4423 was TMEM47.

Interpretation

We found that MS-associated SNPs in sequence determinants of pri-miRNA processing can affect the expression of mature miRNAs. Our findings complement the existing literature on the dysregulation of miRNAs in MS. Further studies on the maturation and function of miRNAs in different cell types and tissues may help to gain a more detailed functional understanding of the genetic basis of MS.

Funding

This study was funded by the Rostock University Medical Center (FORUN program, grant: 889002), Sanofi Genzyme (grant: GZ-2016-11560) and Merck Serono GmbH (Darmstadt, Germany, an affiliate of Merck KGaA, CrossRef Funder ID: 10.13039/100009945, grant: 4501860307). NB was supported by the Stiftung der Deutschen Wirtschaft (sdw) and the FAZIT foundation. EP was supported by the Landesgraduiertenförderung Mecklenburg-Vorpommern.

Cytosolic microRNA-inducible nuclear translocation of Cas9 protein for disease-specific genome modification

MicroRNA-dependent mRNA decay plays an important role in gene silencing by facilitating posttranscriptional and translational repression. Inspired by this intrinsic nature of microRNA-mediated mRNA cleavage, here, we describe a microRNA-targeting mRNA as a switch platform called mRNA bridge mimetics to regulate the translocation of proteins. We applied the mRNA bridge mimetics platform to Cas9 protein to confer it the ability to translocate into the nucleus via cleavage of the nuclear export signal. This system performed programmed gene editing in vitro and in vivo. Combinatorial treatment with cisplatin and miR-21-EZH2 axis-targeting CRISPR Self Check-In improved sensitivity to chemotherapeutic drugs in vivo. Using the endogenous microRNA-mediated mRNA decay mechanism, our platform is able to remodel a cell's natural biology to allow the entry of precise drugs into the nucleus, devoid of non-specific translocation. The mRNA bridge mimetics strategy is promising for applications in which the reaction must be controlled via intracellular stimuli and modulates Cas9 proteins to ensure safe genome modification in diseased conditions.

Genome-Wide TSS Distribution in Three Related Clostridia with Normalized Capp-Switch Sequencing

Transcription initiation is a tightly regulated process that is crucial for many aspects of prokaryotic physiology. High-throughput transcription start site (TSS) mapping can shed light on global and local regulation of transcription initiation, which in turn may help us understand and predict microbial behavior. In this study, we used Capp-Switch sequencing to determine the TSS positions in the genomes of three model solventogenic clostridia: Clostridium acetobutylicum ATCC 824, C. beijerinckii DSM 6423, and C. beijerinckii NCIMB 8052. We first refined the approach by implementing a normalization pipeline accounting for gene expression, yielding a total of 12,114 mapped TSSs across the species. We further compared the distributions of these sites in the three strains. Results indicated similar distribution patterns at the genome scale, but also some sharp differences, such as for the butyryl-CoA synthesis operon, particularly when comparing C. acetobutylicum to the C. beijerinckii strains. Lastly, we found that promoter structure is generally poorly conserved between C. acetobutylicum and C. beijerinckii. A few conserved promoters across species are discussed, showing interesting examples of how TSS determination and comparison can improve our understanding of gene expression regulation at the transcript level. IMPORTANCE Solventogenic clostridia have been employed in industry for more than a century, initially being used in the acetone-butanol-ethanol (ABE) fermentation process for acetone and butanol production. Interest in these bacteria has recently increased in the context of green chemistry and sustainable development. However, our current understanding of their genomes and physiology limits their optimal use as industrial solvent production platforms. The gene regulatory mechanisms of solventogenesis are still only partly understood, impeding efforts to increase rates and yields. Genome-wide mapping of transcription start sites (TSSs) for three model solventogenic Clostridium strains is an important step toward understanding mechanisms of gene regulation in these industrially important bacteria.

A flavivirus-inducible gene expression system that modulates broad-spectrum antiviral activity against dengue and Zika viruses

Incidence of dengue virus (DENV) and Zika virus (ZIKV), two mosquito-borne flaviviruses, is increasing in large parts of the world. Vaccination and medication for these diseases are unsatisfactory. Here, we developed a novel antiviral approach, using a virus-inducible gene expression system, to block virus replication and transmission. Constructs containing the smallest replication units of dengue virus serotype 2 (DENV2) with negative-stranded DENV2 artificial genomes and genes of interest were established in an Aedes aegypti cell line, resulting in expression of target genes after DENV2 infection. Green fluorescent protein (GFP) assays confirmed the system was virus-inducible. When we used one of two apoptosis-related genes, A. aegypti michelob_x (AaMx) and inhibitor of apoptosis (IAP)-antagonist michelob_x-like protein (AaIMP) instead of GFP, the production of viral RNA and proteins were inhibited for all five viruses tested (DENV1-4 and ZIKV), and effector caspase activity was induced. The system thus inhibited the production of infectious virus particles in vitro, and in mosquitoes it did so after DENV2 infection. This is a novel broad-spectrum antiviral approach using a flavivirus-inducible gene-expression system, which could lead to new avenues for mosquito-borne disease control.

SPENCER: a comprehensive database for small peptides encoded by noncoding RNAs in cancer patients

As an increasing number of noncoding RNAs (ncRNAs) have been suggested to encode short bioactive peptides in cancer, the exploration of ncRNA-encoded small peptides (ncPEPs) is emerging as a fascinating field in cancer research. To assist in studies on the regulatory mechanisms of ncPEPs, we describe here a database called SPENCER (http://spencer.renlab.org). Currently, SPENCER has collected a total of 2806 mass spectrometry (MS) data points from 55 studies, covering 1007 tumor samples and 719 normal samples. Using an MS-based proteomics analysis pipeline, SPENCER identified 29 526 ncPEPs across 15 different cancer types. Specifically, 22 060 of these ncPEPs were experimentally validated in other studies. By comparing tumor and normal samples, the identified ncPEPs were divided into four expression groups: tumor-specific, upregulated in cancer, downregulated in cancer, and others. Additionally, since ncPEPs are potential targets for neoantigen-based cancer immunotherapy, SPENCER also predicted the immunogenicity of all the identified ncPEPs by assessing their MHC-I binding affinity, stability, and TCR recognition probability. As a result, 4497 ncPEPs curated in SPENCER were predicted to be immunogenic. Overall, SPENCER will be a useful resource for investigating cancer-associated ncPEPs and may boost further research in cancer.

OUP accepted manuscript

DNA-encoded library (DEL) technology is a powerful tool for small molecule identification in drug discovery, yet the reported DEL selection strategies were applied primarily on protein targets in either purified form or in cellular context. To expand the application of this technology, we employed DEL selection on an RNA target HIV-1 TAR (trans-acting responsive region), but found that the majority of signals were resulted from false positive DNA–RNA binding. We thus developed an optimized selection strategy utilizing RNA patches and competitive elution to minimize unwanted DNA binding, followed by k-mer analysis and motif search to differentiate false positive signal. This optimized strategy resulted in a very clean background in a DEL selection against Escherichia coli FMN Riboswitch, and the enriched compounds were determined with double digit nanomolar binding affinity, as well as similar potency in functional FMN competition assay. These results demonstrated the feasibility of small molecule identification against RNA targets using DEL selection. The developed experimental and computational strategy provided a promising opportunity for RNA ligand screening and expanded the application of DEL selection to a much wider context in drug discovery.

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as PBC.

Mutate-and-chemical-shift-fingerprint (MCSF) to characterize excited states in RNA using NMR spectroscopy

It is important to understand the dynamics and higher energy structures of RNA, called excited states, to achieve better understanding of RNA function. R1ρ relaxation dispersion NMR spectroscopy (RD) determines chemical shift differences between the most stable, ground state and the short-lived, low-populated excited states. We describe a procedure for deducing the excited state structure from these chemical shift differences using the mutate-and-chemical-shift-fingerprint (MCSF) method, which requires ~2-6 weeks and moderate understanding of NMR and RNA structure. We recently applied the MCSF methodology to elucidate the excited state of microRNA 34a targeting the SIRT1 mRNA and use this example to demonstrate the analysis. The protocol comprises the following steps: (i) determination of the secondary structure of the excited state from RD chemical shift data, (ii) design of trapped excited state RNA, (iii) validation of the excited state structure by NMR, and (iv) MCSF analysis comparing the chemical shifts of the trapped excited state with the RD-derived chemical shift differences. MCSF enables observation of the short-lived RNA structures, which can be functionally and structurally characterized by entrapment.

RNAPhaSep: a resource of RNAs undergoing phase separation

Liquid-liquid phase separation (LLPS) partitions cellular contents, underlies the formation of membraneless organelles and plays essential biological roles. To date, most of the research on LLPS has focused on proteins, especially RNA-binding proteins. However, accumulating evidence has demonstrated that RNAs can also function as ‘scaffolds’ and play essential roles in seeding or nucleating the formation of granules. To better utilize the knowledge dispersed in published literature, we here introduce RNAPhaSep (http://www.rnaphasep.cn), a manually curated database of RNAs undergoing LLPS. It contains 1113 entries with experimentally validated RNA self-assembly or RNA and protein co-involved phase separation events. RNAPhaSep contains various types of information, including RNA information, protein information, phase separation experiment information and integrated annotation from multiple databases. RNAPhaSep provides a valuable resource for exploring the relationship between RNA properties and phase behaviour, and may further enhance our comprehensive understanding of LLPS in cellular functions and human diseases.

Distinct Evolutionary Profiles and Functions of microRNA156 and microRNA529 in Land Plants

MicroRNA156 (miR156) and miR529 have high sequence similarity and recognize overlapping sites in the same target genes, SQUAMOSA promoter binding protein-like (SPL or SBP box) genes, making it difficult to accurately distinguish their roles in regulatory networks that affect numerous biological functions. Here, we collected data about miR156 and miR529 family members from representative land plants and performed sequence comparisons, phylogenetic analysis, small RNA sequencing, and parallel analysis of RNA ends (PARE) analysis to dissect their evolutionary and functional differences. Although miR156 and miR529 are highly similar, there are differences in their mismatch-sensitive regions, which are essential for target recognition. In land plants, miR156 precursors are conserved mainly within the hairpin region, whereas miR529 precursors are conserved outside the hairpin region, including both the 5’ and 3’ arms. Phylogenetic analysis showed that MIR156 and MIR529 evolved independently, through divergent evolutionary patterns. The two genes also exhibit different expression patterns, with MIR529 preferentially expressed in reproductive tissues and MIR156 in other tissues. PARE analysis revealed that miR156 and miR529 possess specific targets in addition to common targets in maize, pointing to functional differences between them. Based on our findings, we developed a method for the rapid identification of miR529 and miR156 family members and uncovered the evolutionary divergence of these families, providing insights into their different regulatory roles in plant growth and development.

The temperature-regulated DEAD-box RNA helicase CrhR interactome: Autoregulation and photosynthesis-related transcripts

RNA helicases play crucial functions in RNA biology. In plants, RNA helicases are encoded by large gene families, performing roles in abiotic stress responses, development, the post-transcriptional regulation of gene expression as well as house-keeping functions. Several of these RNA helicases are targeted to the organelles, mitochondria and chloroplasts. Cyanobacteria are the direct evolutionary ancestors of plant chloroplasts. The cyanobacterium Synechocystis 6803 encodes a single DEAD-box RNA helicase, CrhR, that is induced by a range of abiotic stresses, including low temperature. Though the ΔcrhR mutant exhibits a severe cold-sensitive phenotype, the physiological function(s) performed by CrhR have not been described. To identify transcripts interacting with CrhR, we performed RNA co-immunoprecipitation with extracts from a Synechocystis crhR deletion mutant expressing the FLAG-tagged native CrhR or a K57A mutated version with an anticipated enhanced RNA binding. The composition of the interactome was strikingly biased towards photosynthesis-associated and redox-controlled transcripts. A transcript highly enriched in all experiments was the crhR mRNA, suggesting an auto-regulatory molecular mechanism. The identified interactome explains the described physiological role of CrhR in response to the redox poise of the photosynthetic electron transport chain and characterizes CrhR as an enzyme with a diverse range of transcripts as molecular targets.