Small RNA Profiling Research Articles

Differential Expression of miRNAs Involved in Response to Candidatus Liberibacter asiaticus Infection in Mexican Lime at Early and Late Stages of Huanglongbing Disease

Huanglongbing (HLB) is one of the most destructive diseases threatening citriculture worldwide. This disease has been associated with α-proteobacteria species, namely Candidatus Liberibacter. Due to the unculturable nature of the causal agent, it has been difficult to mitigate the disease, and nowadays a cure is not available. MicroRNAs (miRNAs) are key regulators of gene expression, playing an essential role in abiotic and biotic stress in plants including antibacterial responses. However, knowledge derived from non-model systems including Candidatus Liberibacter asiaticus (CLas)-citrus pathosystem remains largely unknown. In this study, small RNA profiles from Mexican lime (Citrus aurantifolia) plants infected with CLas at asymptomatic and symptomatic stages were generated by sRNA-Seq, and miRNAs were obtained with ShortStack software. A total of 46 miRNAs, including 29 known miRNAs and 17 novel miRNAs, were identified in Mexican lime. Among them, six miRNAs were deregulated in the asymptomatic stage, highlighting the up regulation of two new miRNAs. Meanwhile, eight miRNAs were differentially expressed in the symptomatic stage of the disease. The target genes of miRNAs were related to protein modification, transcription factors, and enzyme-coding genes. Our results provide new insights into miRNA-mediated regulation in C. aurantifolia in response to CLas infection. This information will be useful to understand molecular mechanisms behind the defense and pathogenesis of HLB.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc.

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

Regenerative capacity of trophoblast stem cell-derived extracellular vesicles on mesenchymal stem cells

BackgroundHuman mesenchymal stem cells (MSCs) are therapeutic for clinical applications because of their excellent immunomodulatory and multiple lineage differentiation abilities at tissue injury sites. However, insufficient number of cells and lack of regenerative properties during in vitro expansion still limit the clinical applicability of MSC therapies. Here, we demonstrated a preconditioning strategy with trophoblast stem cell-derived extracellular vesicles (TSC-EVs) to boost the proliferation and regenerative capacity of MSCs.MethodsWe employed cell proliferation analyses such as CCK8 and BrdU assays to determine the proliferation-promoting role of TSC-EVs on MSCs. Osteogenic effects of TSC-EVs on MSCs were assessed by alkaline phosphatase (ALP) activity, calcium assays, and calvarial bone defect animal models. For skin regenerative effects, skin wound mice model was exploited to analyze wound-healing rate in this study, as well as immunofluorescence and histological staining evaluates. We also performed the small RNA profiling and RNA-sequencing analyzes to understand the cellular mechanism of TSC-EVs on MSCs.ResultsTSC-EVs significantly promoted MSC proliferation under xeno-free conditions and facilitated the therapeutic effects of MSCs, including osteogenesis, anti-senescence, and wound healing. Transcriptomic analysis also provided evidence that specific microRNAs in TSC-EVs and differentially expressed genes (DEGs) in TSC-EV-treated MSCs showed the possibility of TSC-EVs triggering the regenerative abilities of MSCs with cytokine interaction. Hence, we found that NGF/Akt signaling mediated the regenerative effects of TSC-EVs on MSCs as a particular cellular signaling pathway.ConclusionThe results of this study demonstrated the functional properties of TSC-EVs on MSCs for MSC-based therapeutic applications, suggesting that TSC-EVs may serve as a potential preconditioning source for MSC therapy in the clinical field of regenerative medicine.Graphical abstract

Gluten-free diet affects fecal small non-coding RNA profiles and microbiome composition in celiac disease supporting a host-gut microbiota crosstalk

ABSTRACT Current treatment for celiac disease (CD) is adhering to a gluten-free diet (GFD), although its long-term molecular effects are still undescribed. New molecular features detectable in stool may improve and facilitate noninvasive clinical management of CD. For this purpose, fecal small non-coding RNAs (sncRNAs) and gut microbiome profiles were concomitantly explored in CD subjects in relation to strict (or not) GFD adherence over time. In this observational study, we performed small RNA and shotgun metagenomic sequencing in stool from 63 treated CD (tCD) and 3 untreated subjects as well as 66 sex- and age-matched healthy controls. tCD included 51 individuals on strict GFD and with negative transglutaminase (TG) serology (tCD-TG-) and 12 symptomatic with not strict/short-time of GFD adherence and positive TG serology (tCD-TG+). Samples from additional 40 healthy adult individuals and a cohort of 19 untreated pediatric CD subjects and 19 sex/age matched controls were analyzed to further test the outcomes. Several miRNA and microbial profiles were altered in tCD subjects (adj. p < .05). Findings were validated in the external group of adult controls. In tCD-TG-, GFD duration correlated with five miRNA levels (p < .05): for miR-4533-3p and miR-2681-3p, the longer the diet adherence, the less the expression differed from controls. tCD-TG+ and untreated pediatric CD patients showed a similar miRNA dysregulation. Immune-response, trans-membrane transport and cell death pathways were enriched in targets of identified miRNAs. Bifidobacterium longum, Ruminococcus bicirculans, and Haemophilus parainfluenzae abundances shifted (adj. p < .05) with a progressive reduction of denitrification pathways with GFD length. Integrative analysis highlighted 121 miRNA-bacterial relationships (adj. p < .05). Specific molecular patterns in stool characterize CD subjects, reflecting either the long-term GFD effects or the gut inflammatory status, in case of a not strict/short-time adherence. Our findings suggest novel host-microbial interplays and could help the discovery of biomarkers for GFD monitoring over time.

Differential expression profiling of tRNA-Derived small RNAs and their potential roles in methamphetamine self-administered rats.

Transfer RNA-derived small RNAs (tsRNAs) are a novel class of short, non-coding RNAs that are closely associated with the pathogenesis of various diseases. Accumulating evidence has demonstrated their critical functional roles as regulatory factors in gene expression regulation, protein translation regulation, regulation of various cellular activities, immune mediation, and response to stress. However, the underlying mechanisms by which tRFs & tiRNAs affect methamphetamine-induced pathophysiological processes are largely unknown. In this study, we used a combination of small RNA sequencing, quantitative reverse transcription-polymerase chain reaction (qRT‒PCR), bioinformatics, and luciferase reporter assays to screen the expression profiles and identify the functional roles of tRFs and tiRNAs in the nucleus accumbens (NAc) of methamphetamine self-administration rat models. A total of 461 tRFs & tiRNAs were identified in the NAc of rats after 14days of methamphetamine self-administration training. Of those, 132 tRFs & tiRNAs were significantly differentially expressed: 59 were significantly upregulated, whereas 73 were significantly downregulated in the rats with methamphetamine self-administration. Decreased expression levels of tiRNA-1-34-Lys-CTT-1 and tRF-1-32-Gly-GCC-2-M2, as well as increased expression levels of tRF-1-16-Ala-TGC-4 in the METH group compared with the saline control were validated by using RT‒PCR. Then, bioinformatic analysis was performed to analyse the possible biological functions of tRFs & tiRNAs in methamphetamine-induced pathogenesis. Furthermore, tRF-1-32-Gly-GCC-2-M2 was identified to target BDNF using the luciferase reporter assay. An altered tsRNA expression pattern was proven, and tRF-1-32-Gly-GCC-2-M2 was shown to be involved in methamphetamine-induced pathophysiologic processes by targeting BDNF. The current study provides new insights for future investigations to explore the mechanisms and therapeutic methods for methamphetamine addiction.

Small RNA sequencing andidentification ofAndrographis paniculata miRNAs withpotential cross‑kingdom human gene targets.

Cross-species post-transcriptional regulatory potential of plant derived small non-coding microRNAs (miRNAs) has been well documented by plenteous studies. MicroRNAs are transferred to host cells via oral ingestion wherein they play a decisive role in regulation of host genes; thus, miRNAs have evolved as the nascent bioactive molecules imparting pharmacological values to traditionally used medicinal plants. The present study aims to investigate small RNA profiling in order to uncover the potential regulatory role of miRNAs derived from Andrographis paniculata, one of the most widely used herb by tribal communities for liver disorders and document the pharmacological properties of A. paniculata miRNAs. In this study, high-throughput sequencing method was used to generate raw data, ~ 60 million sequences were generated from A. paniculata leaves. Using computational tools and bioinformatics approach, analyses of 3,480,097 clean reads resulted in identification of 3440 known and 51 putative novel miRNAs regulating 1365 and 192 human genes respectively. Remarkably, the identified plausible novel miRNAs apa-miR-5, apa-miR-1, apa-miR-26, and apa-miR-30 are projected to target significant host genes including CDK6, IKBKB, TRAF3, CHD4, MECP2, and ADIPOQ. Subsequent annotations revealed probable involvement of the target genes in various pathways for instance p38-MAPK, AKT, AMPK, NF-Kβ, ERK, WNT signalling, MYD88 dependant cascade, and pathways in cancer. Various diseases such as human papilloma virus infection, Alzheimer's, Non-alcoholic Fatty Liver, Alcoholic liver diseases, HepatoCellular Carcinoma (HCC), and numerous other cancers were predominantly found to be linked with target genes. Our findings postulate novel interpretations regarding modulation of human transcripts by A. paniculata miRNAs and exhibit the regulation of human diseases by plant-derived miRNAs. Though our study elucidates miRNAs as novel therapeutic agents, however, experimental validations for assessment of therapeutic potential of these miRNAs are still warranted.

Sexual differentiation in dioecious cucurbits - A molecular perspective

Understanding the evolution and genetics of sex differences is significant in both pure and applied research. Sex determination in the family Cucurbitaceae is complex owing to the existence of multiple sexual phenotypes and evolutionary changes between different systems as well as variability in chromosome numbers and structures. Dioecy (or the separate sexed condition) is widespread in angiosperms and the recent availability of molecular markers, next-generation based genome sequencing and small RNA profiling had advanced understanding of the evolution and genetics of sexual dimorphism in several dioecious plant species. Many cucurbits including some economically important species are dioecious. In this review, we summarize knowledge about the molecular mechanisms associated with sexual differentiation in dioecious cucurbits including sexlinked regions and genes identified using molecular markers. We also discuss the roles of small regulatory molecules, phytohormones and transposon variation affecting sex phenotypes in dioecious cucurbits for commercial exploitation.

Small RNA Profiling by Next-Generation Sequencing Using High-Definition Adapters.

Next-generation sequencing (NGS) of small RNA (sRNA) cDNA libraries permits the identification and characterization of sRNA species de novo. However, the method through which these libraries are constructed can often introduce artifacts such as over- or underrepresentation of specific sequences or adapter oligonucleotides due to sequence biases held by the enzymes used. In this chapter we describe a protocol for sRNA library construction making use of high-definition (HD) adapters for the Illumina sequencing platform, which reduce ligation bias. This protocol leads to drastically reduced direct 5'/3' adapter ligation products and can be used for the synthesis of sRNA libraries from total RNA or sRNA of various plant, animal, and fungal samples. This protocol also includes a method for total RNA extraction from plant leaf and cultured cells or body fluids.

Role of Royal Jelly Treated Adipose‐Derived Stem Cell‐Extracellular Vesicles on Fibroblast Proliferation, Migration, and Collagen Production

Extracellular vesicles (EVs) secreted from adipose‐derived stem cells (ADSCs) are known to exhibit collagen synthesis and migration activity in fibroblasts. These EVs can be used as a regenerative medicine and for cosmetic medicinal materials. Royal jelly (RJ) is produced from the hypopharyngeal and mandibular glands of honeybees (nurse bee), which contains apicin and 10‐hydroxy‐2‐decenoic acid, compounds that promote fibroblast proliferation and collagen and elastin synthesis. Here, we investigated whether RJ could further enhance the physiological function of ADSC‐EVs on fibroblasts. The findings confirmed that lyophilized RJ and enzyme‐treated RJ enhanced the secretion of EVs by ADSCs, increasing the proliferation, migration, and collagen synthesis of fibroblasts compared to those induced by conventional ADSC‐EVs. Bioinformatics analysis was conducted using various software programs, based on the data obtained from comprehensive gene expression and small RNA profiling analyses through next‐generation sequencing; the results suggested that the RJ treatment of ADSC‐EVs significantly enhanced the expression of genes related to extracellular matrix composition. Our findings suggest that the ADSC‐EVs qualitatively altered using RJ could offer a more effective therapeutic material for regenerative and cosmetic medicine than conventional ADSCs‐EVs.

Expression profiles of small non-coding RNAs in breast cancer tumors characterize clinicopathological features and show prognostic and predictive potential

Precision medicine approaches are required for more effective therapies for cancer. As small non-coding RNAs (sncRNAs) have recently been suggested as intriguing candidates for cancer biomarkers and have shown potential also as novel therapeutic targets, we aimed at profiling the non-miRNA sncRNAs in a large sample set to evaluate their role in invasive breast cancer (BC). We used small RNA sequencing and 195 fresh-frozen invasive BC and 22 benign breast tissue samples to identify significant associations of small nucleolar RNAs, small nuclear RNAs, and miscellaneous RNAs with the clinicopathological features and patient outcome of BC. Ninety-six and five sncRNAs significantly distinguished (Padj < 0.01) invasive local BC from benign breast tissue and metastasized BC from invasive local BC, respectively. Furthermore, 69 sncRNAs significantly associated (Padj < 0.01) with the tumor grade, hormone receptor status, subtype, and/or tumor histology. Additionally, 42 sncRNAs were observed as candidates for prognostic markers and 29 for predictive markers for radiotherapy and/or tamoxifen response (P < 0.05). We discovered the clinical relevance of sncRNAs from each studied RNA type. By introducing new sncRNA biomarker candidates for invasive BC and validating the potential of previously described ones, we have guided the way for further research that is warranted for providing novel insights into BC biology.

Circulating Small Noncoding RNA Profiling as a Potential Biomarker of Atherosclerotic Plaque Composition in Type 1 Diabetes.

Cardiovascular disease (CVD) accounts for most deaths in patients with type 1 diabetes (T1D); however, the determinants of plaque composition are unknown. miRNAs regulate gene expression, participate in the development of atherosclerosis, and represent promising CVD biomarkers. This study analyzed the circulating miRNA expression profile in T1D with either carotid calcified (CCP) or fibrous plaque (CFP). Circulating small noncoding RNAs were sequenced and quantified using next-generation sequencing and bioinformatic analysis in an exploratory set of 26 subjects with T1D with CCP and in 25 with CFP. Then, in a validation set of 40 subjects with CCP, 40 with CFP, and 24 control subjects with T1D, selected miRNA expression was measured by digital droplet PCR. Putative gene targets enriched for pathways implicated in atherosclerosis/vascular calcification/diabetes were analyzed. The patients' main clinical characteristics were also recorded. miR-503-5p, let-7d-5p, miR-106b-3p, and miR-93-5p were significantly upregulated, while miR-10a-5p was downregulated in patients with CCP compared with CFP (all fold change >±1.5; P < 0.05). All candidate miRNAs showed a significant correlation with LDL-cholesterol, direct for the upregulated and inverse for the downregulated miRNA, in CCP. Many target genes of upregulated miRNAs in CCP participate in osteogenic differentiation, apoptosis, inflammation, cholesterol metabolism, and extracellular matrix organization. These findings characterize miRNAs and their signature in the regulatory network of carotid plaque phenotype in T1D, providing new insights into plaque pathophysiology and possibly novel biomarkers of plaque composition.

Small RNA profiling reveals the involvement of microRNA-mediated gene regulation in response to symbiosis in raspberry.

Dark septate endophytes (DSEs) can form reciprocal symbioses with most terrestrial plants, providing them with mineral nutrients in exchange for photosynthetic products. Although the mechanism of plant-DSEs is well understood at the transcriptional level, little is known about their post-transcriptional regulation, and microRNAs (miRNAs) for the symbiotic process of DSE infestation of raspberry have not been identified. In this study, we comprehensively identified the miRNAs of DSE-infested raspberry symbiosis using Illumina sequencing. A total of 361 known miRNAs and 95 novel miRNAs were identified in the roots. Similar to other dicotyledons, most of the identified raspberry miRNAs were 21 nt in length. Thirty-seven miRNAs were differentially expressed during colonization after inoculation with Phialocephala fortinii F5, suggesting a possible role for these miRNAs in the symbiotic process. Notably, two miRNAs (miR171h and miR396) previously reported to be responsive to symbiotic processes in alfalfa also had altered expression during raspberry symbiosis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggests that miRNAs are mainly involved in regulatory mechanisms, such as biological processes, cellular metabolic processes, biosynthesis of secondary metabolites, plant-pathogen interactions, and phytohormone signaling pathways. This study revealed the potential conservation of miRNA-mediated post-transcriptional regulation in symbiotic processes among plants and provides some novel miRNAs for understanding the regulatory mechanisms of DSE-raspberry symbiosis.

Small RNA sequencing and profiling of serum-derived exosomes from African swine fever virus-infected pigs.

African swine fever (ASF) virus (ASFV) is responsible for one of the most severe swine diseases worldwide, with a morbidity rate of up to 100%; no vaccines or antiviral medicines are available against the virus. Exosomal miRNAs from individual cells can regulate the immune response to infectious diseases. In this study, pigs were infected with an ASFV Pig/HN/07 strain that was classified as acute form, and exosomal miRNA expression in the serum of infected pigs was analyzed using small RNA sequencing (small RNA-seq). Twenty-seven differentially expressed (DE) miRNAs were identified in the ASFV-infected pigs compared to that in the uninfected controls. Of these, 10 were upregulated and 17 were downregulated in the infected pigs. All DE miRNAs were analyzed using gene ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and the DE miRNAs were found to be highly involved in T-cell receptor signaling, cGMP-PKG signaling, Toll-like receptor, MAPK signaling, and mTOR signaling pathways. Furthermore, the Cytoscape network analysis identified the network of interactions between DE miRNAs and target genes. Finally, the transcription levels of four miRNA genes (ssc-miR-24-3p, ssc-miR-130b-3p, ssc-let-7a, and ssc-let-7c) were examined using quantitative real-time PCR (qRT-PCR) and were found to be consistent with the small RNA-seq data. These DE miRNAs were associated with cellular genes involved in the pathways related to immune response, virus-host interactions, and several viral genes. Overall, our findings provide an important reference and improve our understanding of ASF pathogenesis and the immune or protective responses during an acute infection in the host.

An atlas of small RNAs from potato.

Small RNAs, including microRNAs (miRNAs), phased secondary small interfering RNAs (phasiRNA), and heterochromatic small interfering RNAs (hc-siRNA) are an essential component of gene regulation. To establish a broad potato small RNA atlas, we constructed an expression atlas of leaves, flowers, roots, and tubers of Desiree and Eva, which are commercially important potato (Solanum tuberosum) cultivars. All small RNAs identified were observed to be conserved between both cultivars, supporting the hypothesis that small RNAs have a low evolutionary rate and are mostly conserved between lineages. However, we also found that a few miRNAs showed differential accumulation between the two potato cultivars, and that hc-siRNAs have a tissue specific expression. We further identified dozens of reproductive and non-reproductive phasiRNAs originating from coding and noncoding regions that appeared to exhibit tissue-specific expression. Together, this study provides an extensive small RNA profiling of different potato tissues that might be used as a resource for future investigations.

The developmental miR-17–92 cluster and the Sfmbt2 miRNA cluster cannot rescue the abnormal embryonic development generated using obstructive epididymal environment-producing sperm in C57BL/6 J mice

BackgroundSperm, during epididymal transit, acquires microRNAs(miRNAs), which are crucial for embryonic development. However, whether sperm miRNAs influenced by an obstructive epididymal environment affect embryonic development remains unknown.MethodThe sham operation and vasectomy were performed in C57BL/6 J mice to create the control group (CON) and the obstructive epididymal environment group(OEE) group, respectively. The morphology of the testis and epididymis was observed using hematoxylin and eosin staining (HE staining) to establish the OEE mice model. The sperm quality test, intracytoplasmic sperm injection (ICSI), and epididymosomes fusion were employed to observe the effect of the obstructive epididymal environment on sperm and resultant embryonic development. The alteration of the sperm small RNA (sRNA) profile was analyzed by sRNA sequencing. RT-qPCR and DNA methylation were applied to observe the effect of obstructive epididymis on the expression of sperm miRNAs. The miRNAs microinjection was used to explore the impacts of sperm miRNAs on embryonic development.ResultsWe confirmed postoperative 8-week mice as the OEE mice model by examining the morphology of the testis and epididymis. In the OEE group, we observed that sperm quality degraded and the development potential of embryos was reduced, which can be saved by the normal epididymal environment. The sperm sRNA sequencing revealed that the expression of the developmental miR-17–92 cluster and the Sfmbt2 miRNA cluster was downregulated in the OEE group. The expression of these two miRNA clusters in epididymis was also downregulated and regulated by DNA methylation. However, the downregulation of either the miR-17–92 cluster or the Sfmbt2 miRNA cluster in normal zygotes did not impair embryonic development.ConclusionThe obstructive epididymal environment influences sperm quality and resultant embryonic development, as well as the abundance of the developmental miR-17–92 cluster and the Sfmbt2 miRNA cluster in sperm, but these miRNA clusters are not the cause of abnormal embryonic development. It implies that epididymis is important in early embryonic development and may play a potential role in sperm epigenome.

Assessing spermatozoal small ribonucleic acids and their relationship to blastocyst development in idiopathic infertile males

Clinical testing strategies for diagnosing male factor infertility are limited. A deeper analysis of spermatozoa-derived factors could potentially diagnose some cases of ‘unexplained infertility’. Spermatozoa carry a rich and dynamic profile of small RNAs, which have demonstrated potential developmental importance and association with fertility status. We used next-generation sequencing to correlate sperm small RNA profiles of normozoospermic males (n = 54) with differing blastocyst development rates, when using young donor oocytes. While ribosomal RNAs accounted for the highest number of sequencing reads, transfer RNA fragments of tRNAGly/GCC and tRNAVal-CAC were the most abundant sequences across all sperm samples. A total of 324 small RNAs were differentially expressed between samples with high (n = 18) and low (n = 14) blastocyst rates (p-adj < 0.05). Ninety three miRNAs were differentially expressed between these groups (p-adj < 0.05). Differentially expressed transfer RNA fragments included: 5'-tRF-Asp-GTC; 5'-tRF-Phe-GAA; and 3'-tRF-Ser-GCA. Differentially expressed miRNAs included: let-7f-2-5p; miR-4755-3p; and miR-92a-3p. This study provides the foundation on which to validate a clinical panel of fertility-related sperm small RNAs, as well as to pursue potential mechanisms through which they alter blastocyst development.

The Landscape of Small Nucleolar RNA Expression in Multiple Myeloma

Background Small nucleolar RNAs (snoRNAs) are non-coding RNAs involved in the maturation of ribosome RNAs and various other biological processes. Emerging evidence demonstrates that snoRNAs play an essential role in tumor initiation and progression. However, their role in multiple myeloma (MM) pathogenesis remains poorly characterized. To date, the profiling of small RNAs in MM has mainly focused on microRNAs and approaches are confined to PCR- or array-based techniques. Here, we used next generation sequencing (NGS) to explore the snoRNA expression patterns of MM patients and to identify the association between snoRNAs and mRNAs, molecular subgroups of MM and clinical parameters. Methods Purified CD138+ cells from 71 newly diagnosed MMs and 4 refractory/relapsed MMs were included. Three CD19+ B cell samples from healthy donors were used as control. Total RNA was isolated with the AllPrep DNA/RNA kit (Qiagen, Germany). The mRNA and snoRNA NGS libraries were prepared using the Illumina TruSeq stranded mRNA and the NEBNext multiplex small RNA library kit, respectively. Sequencing of mRNA and snoRNA was performed on an Illumina NextSeq and NovaSeq 6000 PE 100 S1 platform, respectively. Paired-end RNA-seq reads were mapped to the hg19 reference genome using STAR. Gene expressions were quantified using featureCounts. Differentially expressed mRNAs and snoRNAs were identified using DESeq2. The Kaplan-Meier method was used for survival analyses. Progression-free survival (PFS) time was measured from enrollment to relapse or death from any cause, whichever occurred first. Results Using NGS, we observed expression of 625 full-length snoRNAs. Compared to normal B cells high levels of snoRNA expression were detected in MM plasma cells. Looking for associations between snoRNA expression and molecular subgroups of MM, we found an upregulation of 50 snoRNAs in patients with a hyperdiploid karyotype. The pattern of upregulation supported a regulatory mechanism in this molecular subgroup that goes beyond a pure copy-number effect. Even in patients with multiple trisomies the overwhelming majority (82%) of upregulated snoRNAs were located on chromosome 15, including the SNORD115 family on 15q11. In patients with gain of 1q21, all snoRNAs located on 1q25 were upregulated and the highest expression level of those snoRNAs appeared in patients with amp1q, indicating a copy-number effect. SNORD78, which was recently shown to impact outcome in prostate cancer and lung cancer, was significantly overexpressed in gain1q patients and defined a subgroup with inferior PFS (p =0.019) and overall survival (p = 0.011). SCARNA22, which was upregulated in patients with a t(4;14), was another snoRNA with a negative impact on PFS (31.06 months vs unreached, p = 0.007). To further characterize MM with upregulation of prognostically relevant snoRNAs, we analyzed bulk RNAseq of CD138-enriched MM cells and performed a pathway enrichment analysis. Comparing SNORD78 high expression group with the low group, we found 537 upregulated genes, e.g.SH2D2A, and CDK1, and 229 downregulated genes. Differentially expressed genes were enriched for cell cycle and chromosome segregation. For SCARNA22, we found 963 up- and 219 down-regulated genes. These showed enrichment for cell adhesion (upregulated in SCARNA22high) and protein targeting to the membrane or the endoplasmic reticulum (downregulated in SCARNA22high). Conclusion In this study we provide first insights into the expression patterns of snoRNAs in MM by deep sequencing. We identified snoRNAs dysregulated in MM subgroups, with two of them being associated with poor outcome. This study also shows that expression of particular snoRNAs is associated with specific chromosome gains, which is a common element in MM. Disclosures No relevant conflicts of interest to declare.

Small Non-Coding RNA Profiles of Sorted Plasma Extracellular Vesicles: Technical Approach

Small Non-Coding RNA Profiles of Sorted Plasma Extracellular Vesicles: Technical Approach

Circulating Small Non-Coding RNA Profiles for Premature Cardiovascular Death in Hemodialysis Patients

Circulating Small Non-Coding RNA Profiles for Premature Cardiovascular Death in Hemodialysis Patients

SARS-CoV-2 RNAs are processed into 22-nt vsRNAs in Vero cells.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the global pandemic, resulting in great fatalities around the world. Although the antiviral roles of RNA interference (RNAi) have been well studied in plants, nematodes and insects, the antiviral roles of RNAi in mammalians are still debating as RNAi effect is suspected to be suppressed by interferon (IFN) signaling pathways in most cell types. To determine the role of RNAi in mammalian resistance to SARS-CoV-2, we studied the profiling of host small RNAs and SARS-CoV-2 virus-derived small RNAs (vsRNAs) in the early infection stages of Vero cells, an IFN-deficient cell line. We found that host microRNAs (miRNAs) were dysregulated upon SARS-CoV-2 infection, resulting in downregulation of microRNAs playing antiviral functions and upregulation of microRNAs facilitating viral proliferations. Moreover, vsRNA peaked at 22 nt at negative strand but not the positive strand of SARS-CoV-2 and formed successive Dicer-spliced pattern at both strands. Similar characteristics of vsRNAs were observed in IFN-deficient cell lines infected with Sindbis and Zika viruses. Together, these findings indicate that host cell may deploy RNAi pathway to combat SARS-CoV-2 infection in IFN-deficient cells, informing the alternative antiviral strategies to be developed for patients or tissues with IFN deficiency.