key: cord-0937346-tdhu4fg7 authors: Grehl, Claudius; Schultheiß, Christoph; Hoffmann, Katrin; Binder, Mascha; Altmann, Thomas; Grosse, Ivo; Kuhlmann, Markus title: Detection of SARS-CoV-2 Derived Small RNAs and Changes in Circulating Small RNAs Associated with COVID-19 date: 2021-08-11 journal: Viruses DOI: 10.3390/v13081593 sha: 7d21b18ff75b581fece0a752eb5c41e202a2bc0d doc_id: 937346 cord_uid: tdhu4fg7 Cleavage of double-stranded RNA is described as an evolutionary conserved host defense mechanism against viral infection. Small RNAs are the product and triggers of post transcriptional gene silencing events. Up until now, the relevance of this mechanism for SARS-CoV-2-directed immune responses remains elusive. Herein, we used high throughput sequencing to profile the plasma of active and convalescent COVID-19 patients for the presence of small circulating RNAs. The existence of SARS-CoV-2 derived small RNAs in plasma samples of mild and severe COVID-19 cases is described. Clusters of high siRNA abundance were discovered, homologous to the nsp2 3′-end and nsp4 virus sequence. Four virus-derived small RNA sequences have the size of human miRNAs, and a target search revealed candidate genes associated with ageusia and long COVID symptoms. These virus-derived small RNAs were detectable also after recovery from the disease. The additional analysis of circulating human miRNAs revealed differentially abundant miRNAs, discriminating mild from severe cases. A total of 29 miRNAs were reduced or absent in severe cases. Several of these are associated with JAK-STAT response and cytokine storm. A novel disease has been spreading in a pandemic manner since the end of 2019. This disease was named COVID-19 (coronavirus disease 2019), as it was found to be associated with the zoonotic Severe Acute Respiratory Syndrome Corona Virus type 2 (SARS-CoV-2). SARS-CoV-2 is an enveloped, positive-sense, single-stranded RNA betacoronavirus of the family Coronaviridae. In humans, coronaviruses cause infections of the respiratory tract and can lead to acute respiratory distress syndrome (ARDS), which may result in widespread inflammation of many organs and death [1] . The case fatality ratio of SARS-CoV-2 is approximately 100 times higher than the common influenza infection (case fatality ratio: 0.01-0.04 deaths per 100 positive cases, [2] ). It ranges between 1 and 5 per 100 cases, depending on local medicinal standards (MERS: 34.4, SARS; 9.5 [3] ) and vulnerability, e.g., Trisomy 21 [4] . This high rate is mainly driven by the lack of effective treatment of severe cases. Cleavage of double-stranded RNA is a well-studied host mechanism, involved in viral defense [5, 6] in the plant and animal kingdoms. This initial defense reaction is triggered by the presence of double-stranded RNA (dsRNA). Presence of dsRNA is a typical indication Do we see circulating small RNAs corresponding to parts of viral genome, and can they be used as marker for ongoing or cleared infections? • Do we see circulating small RNAs corresponding to parts of viral genome with distinct features? • Do we see circ-miRNAs that suggest RNA-based regulation/fine-tuning of the host immune response that differs in mild and severe cases, with potential as biomarker [14, 15] ? Blood was collected during March and August 2020 at the University Hospital Halle (Saale) from patients (Caucasian) in close proximity to Halle (Saale) in Germany. Plasma was isolated by centrifugation of whole blood from EDTA vacutainers for 15 min at 2000× g. Samples were stored at −80 • C before use. Blood collection was performed under institutional review board approvals, numbers 2020-039 and 11/17, in accordance with the Declaration of Helsinki. Informed written consent was obtained from all participating patients or legal representatives. In total, 100 µL of plasma were used for the heat inactivated analysis, according to [19] and processed by GenXPro (Frankfurt am Main, Germany). Levels of antibodies directed against the spike (S1) and the nucleocapsid protein (NCP) of SARS-CoV-2 in plasma were determined using the Anti-SARS-CoV-2-ELISA IgA/IgG and Anti-SARS-CoV-2-NCP-ELISA ELISA kits from Euroimmun AG (Lübeck, Germany), according to manufacturer's instructions. Plasma RNA was extracted using an internal purification protocol of GenXPro GmbH, based on silica columns. The small RNA libraries for Illumina sequencing were prepared using the TrueQuant small RNA kit (GenXPro GmbH, Frankfurt am Main, Germany) as described by [20] . The libraries were sequenced using an illumina NextSeq500 instrument (GenXPro GmbH, Frankfurt am Main, Germany). For detection of SARS-CoV-2 associated small RNAs, the raw reads of the individual samples, as well as the grouped samples (mild/severe) were mapped against the SARS-CoV-2 genome sequence (NCBI Genome Reference Sequence NC_045512) in forward and reverse orientation. CLC-Workbench 11 software was used with default settings, allowing one mismatch for mapping of reads against the virus genome. To screen for SARS-CoV-2-associated small RNAs, plasma was collected from eight PCR-positive COVID-19 patients within the first three weeks after symptom onset, as well as from one healthy donor and one patient with bacterial pneumonia as controls. Five of the COVID-19 patients had severe courses and required extracorporeal membrane oxygenation (ECMO) in an intensive care unit (ICU). Three of the five severe cases succumbed to the disease. The remaining three had mild to moderate disease symptoms. Of these, three were sampled during active disease phase within the first 17 days after symptom onset, and one in the convalescence phase 41 days post infection. Two samples during active infection were available for one of the mild cases, with a sample each in the active and in the convalescence phase for one patient. The sampling, clinical and laboratory findings of the cohort and NGS read output are summarized in Table 1, Supplementary Table S1 , Tab Samples and Supplementary Figure S2A ,B. Comprehensive immunological data are described for five of these patients [21] . For the identification of small SARS-CoV-2 derived RNAs, we analyzed plasma of patients with PCR-confirmed presence of SARS-CoV-2 (mild and recovered, n = 5; severe, n = 5) and negative controls (n = 2) without a present or previous SARS-CoV-2 infection. RNA was extracted, the fraction of small RNAs enriched and subjected to library preparation. The read number after adapter trimming ranged from 0.8 to 6.9 million (Table 1 ) with a median of 1.4 million reads per sample. For detection of SARS-CoV-2 associated small RNAs, these raw reads of the individual samples, as well as the grouped samples (severe/mild), were mapped against the SARS-CoV-2 genome sequence in sense and antisense orientation, allowing one mismatch (Supplementary Table S1 for numbers and sequences of mapped reads in sense and antisense). The summarized mapping results are displayed in Table 2 . In mild/recovered and severe samples substantial small virus-derived RNAs are detected. The results show n = 85 (sense: 68, antisense: 17) in mild cases and n = 171 (sense: 118, antisense: 53) in severe cases in both orientations, while the sense orientation is higher in abundance. The size of the mapped RNAs ranged from 15 to 19 nts. The healthy individual did not reveal small RNAs complementary to the SARS-CoV-2 genome. The pneumonia disease control had three reads mapping to less SARS-CoV-2 specific regions. One of these reads was assigned to the polyA structure in the 3 -region of the virus genome, and the second was mapped to the 3 -end of the SARS-CoV-2 genome which is homologous to e.g., Vibrio and Acinetobacter genomes. Acinetobacter baumannii is a known cause of both communityand hospital-acquired infections, including pneumonia [22] . The 15 nt long antisense read shares complete homology to the human ERCC excision repair 6 like-transcript, and might represent a degradation product. The read counts were normalized to RPM (reads per million) for quantification. The abundance of reads per sample is displayed in Figure 1 . The highest number of reads (35 RPM) was detected in the sample of a mild case 10 days after first symptoms (pt68), where the lowest number of read output was also noted. In addition, a high abundance of SARS-CoV-2 reads was found in the severe cases in pt1 and pt44 (5 and 74 days after infection,~10 RPM). In samples of recovered patients (pt7_4 and pt50), virus-associated RNAs were detectable (27 and 41 days after infection, respectively). A low abundance of reads was found across the virus sequence. However, particular clusters of reads were detectable in the 5 -region of the virus genome associated with the nsp2 encoding sequence. It is interesting to note that in all patients where the SARS-CoV-2 specific antibodies were not [yet] present (Supplementary Figure S2B) , the abundance of virus dependent small RNAs was higher (pt1, pt44 and pt68, Supplementary Figure S2C ). Two particular clusters of reads are found to be associated with the SARS-CoV-2 genome sequence. Both are located in the nsp2 encoding region (Figure 2A ), in a forward orientation within ORF1. Reads associated with cluster 1 (CCTTCAATGGGGGATGTCC) were detected in all PCR positive tested samples. A higher abundance of cluster 2 reads (ACCCTAAGAGGGGTGTAT), located upstream of cluster 1, was only detected in severe cases. The normalized abundance of cluster 1 matching reads is shown in Figure S2B . Furthermore, the highest number of reads was found in the sample of a mild case (pt68), 10 days after infection. No reads of this cluster were found in the pneumonia patient and the healthy individual. In samples pt7_4 and pt50, representing recovered cases, reads from this cluster were detectable. A detailed graphical description of cluster 1 and 2 are given in Supplementary Figure S2D . As human microRNAs can have a size between 19 and 24 nts [23] , we were focusing on the 19 nt long small RNA associated with SARS-CoV-2. Here, four sequences were identified, two of them with a higher abundance. Relative positions in the virus genome are displayed in Supplementary Figure S2E . The first sequence with a high abundance was CCTTCAATGGGGGATGTCC (at position 1058 in the virus genome, 10 reads), which was also located in cluster 1. A sequence specific target gene search by (psRNATarget [24] ), adjusted for the human transcriptome, revealed six genes with an expectation score < 4 (Supplementary Table S1 ). Among these six genes, GFRA1 (GDNF family receptor alpha-1; MIM 601496) was found, which has been reported as a receptor required for taste perception [25] and function of the olfactory system [26] . Using the sequence CATATTCAGTGGATGGTTG (pos. 9616, 11 reads) of the higher abundant 19 nucleotide-sized RNA as input for the target search, we found ARL6IP6, TRPM3, HTR2A and FGF2 among the seven top potential targets (Supplementary Table S1 ). ARL6IP6 (MIM 616495) is an ADP-Ribosylation-Like Factor 6 Interacting Protein 6, associated with ischemic stroke [27] . HTR2A (MIM 18135) is a Serotonin 5-HT-2A Receptor with effects on sleep quality [28] . Reduction of TRPM3 (MIM 608961) is described to be associated with chronic fatigue syndrome/Myalgic encephalomyelitis patients [29] , and FGF2 (MIM 134920) has been proposed to play a role in intussusceptive angiogenesis during COVID-19 [30] . The third sequence with low abundance was CTCATTCAAGGAGGACTTA (pos. 24945), with five potential target genes predicted (Supplementary Table S1 ). The fourth sequence was CAATGCTCATGGATTGTTG (pos. 29496) in antisense orientation at the 5 -end of ORF9 encoding the nucleocapsid. For this sequence, eight potential target RNAs were identified (Supplementary Table S1 ). Among these, PSG1 and SLFN13 were found: PSG1 (MIM 176390) is a pregnancy-specific beta-1-glycoprotein 1, regulating P1GF and VREGF-A. The upregulation of both genes has been described as a COVID-19 marker [31] . SLFN13 (MIM 614957) is a member of the SCHLAFEN family, and is suggested to mediate antiviral responses to both influenza A and B RNA virus infections [32] . Circulating small RNAs in the blood are described as potent biomarkers for several diseases and general inflammatory reactions [12, [14] [15] [16] [17] . For a detailed analysis of circulating miRNAs in the sample set, the reads were mapped against the publicly available miRBase [33] , revealing information for long non-coding RNAs, transfer-RNAs, structural RNAs and miRNAs. Principal component analysis (PCA) indicated a common separation of all SARS-CoV-2 infection samples from the healthy individual (HD, Supplementary Figure S3A ). The sample from the pneumonia patient (ICU2, negative control) is among the group of SARS-CoV-2 infection samples. As the patient was affected by a bacterial infection, this might indicate a common pro-inflammatory pattern of circ-miRNA. The patients which succumbed to the disease (pt1, pt25 and pt44) have the most diverse circ-miRNA pattern, and the samples are located apart from the core cluster. Differentially abundant circ-miRNAs were identified by comparison of the set of mild cases (n = 3) with severe cases (n = 5, p < 0.05, DESeq [34] , Supplementary Figure S3B) . Twenty of the thirty-one differentially abundant miRNAs showed a significant (p < 0.05) lower abundance (log2fold change) for severe cases. Eleven circ-miRNAs had positive log2fold-change values, indicating higher abundance in severe cases. In total, we found 1119 Homo sapiens mircroRNAs (hsa-miRNAs), 17,871 human Ensembl transcript ids (ncRNAs) and 392 tRNAs. According to our analysis, 210 ensemble transcripts (including a high number of non-coding RNAs) were differentially abundant, while only four tRNAs (tRNA-His-GTG-1-4, tRNA-Lys-TTT, tRNA-Arg-TCT and tRNA-Ser-GCT:~3 log2fold decreased) showed a significant differential abundance (p < 0.05). The complete list of differentially abundant RNAs is given in Supplementary Table S2 . Differentially abundant circ-miRNAs sorted by log-fold change and association are shown in Table 3 . The available literature reference for their identified context is given in Supplementary Table S2 . The diagrams in Figure 3 describe the detailed read abundance for the top six differentially abundant circ-miRNAs: has-miRNA320a-3hashsa-miRNA629-5p, hsa-miRNA29a-3p, has-miRNA-324-3p, has-miRNA-185-5hasnd hsa-miRNA4516, whereby the first five of the six clearly distinguish mild from severe cases in our cohort. hsa-miRNA320 (Figure 3 ) was reported to be unregulated, in association with heart failure in patients with diabetes mellitus [35] [36] [37] . The entire subset of miRNA320 (b,c,d) was found to be up-regulated in severe cases here. hsa-miRNA29a-3p was reported as a biomarker for the diagnosis of tuberculosis [38] . has-miRNA625-3p is used as biomarker for Malignant Pleural Mesothelioma [39] (up-regulated) and Nontuberculous Mycobacterial Pulmonary Disease [40] (up-regulated). Its counterparhashsa-miRNA625-5p, was reported to suppress inflammatory responses by targeting AKT2 (MIM 164731) in human bronchial epithelial cells [41] . This miRNA was found to be strongly down-regulated (−1.64 log2fold change, p = 0.0468) in all severe cases. hsa-miRNA4516 was reported to be involved in STAT3 (MIM 102582) regulation and also associated with psoriasis [42] . hsa-miRNA19a (Figure 3 ) was absent in the severe samples analyzed. hsa-miR19 is involved in many steps of the inflammatory response and coupled to cytokine response [43] . Due to its absence in severe samples and its low read number, hsa-miR19 did not show up among the DESeq generated table of significant DE miRNAs. In addition, hsa-miR-122-3p and hsa-miR-122b-5p were found among the differential abundant circ-miRs (−4 log2fold change, p = 0.02), which are already proposed as a potential target of human miRNAs in SARS-CoV-2 for RNA-based drug discovery [44] . The mapping of the plasma-derived small RNAs reads against the SARS-CoV-2 genome resulted in the identification of SARS-CoV-2 derived sequences. These sequences are uniquely present in samples of SARS-CoV-2 positive patients, and are considered as SARS-CoV-2 derived small RNAs. The size of the identified reads ranged from 15 to 19 nts (Supplementary Table S1 ). These reads are relatively short in length compared to the expected size of 21 +/− 1 nts [45, 46] for DICER products. Up to now, the length of 21 +/− 1 nts has been considered a hallmark of DICER products. This difference in size might be an explanation of why other studies failed to observe differences in virus-derived small RNAs [47] . However, based on their size, we were not able to identify the mechanistic origin of these molecules in the cell-free plasma. One option might be that partial degradation, cleavage or misprocessing is generating artificial byproducts of virus replication. Another origin might be processing or generation of smaller products by DICER and miRNA processing, followed by an export process. This is supported by the presence of human miRNAs (ranging from 16 to 27 nts) with a different size than the canonical 21mers miRNAs as products of the human DICER [23] in the serum. Some of the identified sequences are overlapping with in silico predicted regions encoding small RNAs [48] , based on folding and cleavage mechanisms. As shown in Supplementary Figure S2 , in all our tested patients we saw a negative relation between the presence of SARS-CoV-2 specific antibodies and the abundance of SARS-CoV-2 derived small RNAs. In the samples (pt1, pt67-1 and pt68) taken at an early timepoint after infection, a substantial/high abundance of small RNAs was detectable. Although the number of analyzed sample was low, the presence of virus-derived small RNAs suggests dicing of viral dsRNA as an early defense mechanism (4-10 days after infection) for virus suppression. Cleavage of dsRNA might represent a first hurdle, inactivating the intruding viruses, while the antibody defense is not yet established. The identified virus-derived sequences of 19 nts length have the potential to mimic human miRNAs and suppress target mRNA function, based on their complementarity. The encoding of miRNA-like small RNAs in a viral genome [7] was described for the bovine leukemia virus [49] , the Ebola virus [50] and lately for SARS-CoV-2 [13] . Among the predicted target mRNAs, GFRA1 (GDNF family receptor alpha-1), HTR2A, TRPM3 and FGF2 are found. Post transcriptional gene suppression mediated by the virus-derived small RNAs of these transcripts might explain some of the COVID-19 associated symptoms. These transient symptoms, such as loss of taste and smell, insomnia, chronic fatigue syndrome and elevated risk for thrombosis [51] might be associated with the transient suppression of the identified transcripts. This finding opens a novel perspective on the molecular regulation and importance of small virus-derived RNAs as a regulator for virus related processes in the host. In contrast to [13] , no virus-derived small RNAs larger than 20 were detected in our serum samples. The comparison of the SARS-CoV-2 mapped reads with the human genome validated the virus specificity. Only the low-complexity reads with mappings to the unstructured polyA virus region showed an additional mapping to human sequences. The presence of these sequences in the samples of control patients without a SARS-CoV-2 infection is also indicative of their origin from the human genome, and therefore an artificial random alignment to the virus genome. Therefore, these reads are considered misaligned and are excluded from the SARS-CoV-2 derived small RNAs. Considering the ongoing evolution of the SARS-CoV-2, reflected by its entropy and mutation rate [52] , it is interesting to note that a high level of entropy (0.5) is detectable at position 1058-1060 in close proximity to the identified small RNA cluster. However, up to now, the mutations of novel detected variants have not affected the region covered by the small RNA cluster. Due to the small number of samples, we were only able to show the presence of the virus-derived small RNAs, and it was not possible to show a quantitative correlation with patient gender or duration of infection. A cluster of reads with a higher abundance was found homologous to the 3 -end of the nsp2 encoding region. These SARS-CoV-2 specific reads might be suitable as biomarkers for an early, ongoing and cleared infection. The detection of a virus-derived microRNA as a biomarker for severe COVID-19 cases was described [13] . The size and abundance of the clustered small RNAs are relatively high, and comparable to some human miRNAs. As such, they can be detected and quantified by stem-loop qPCR [53] . As identified in samples pt7_4 and pt50, they are substantially detectable in patients that have recovered from COVID-19 21-and 48days post-infection. In these previously SARS-CoV-2 PCR positive patients, no RNA virus particles were detectable with standard PCR methods anymore when the samples for the siRNA analyses were taken. Although the recurrence of RNAemia is still under debate [54] [55] [56] , it is generally accepted that viral RNAs can be detectable, independent of COVID-19, up to three weeks after infection [57] . Therefore, the presence of the detected virus-derived small RNAs is in in agreement with reported RNA presence. The abundance of the small RNAs from the cluster region ( Figure 2B ) is comparable to other miRNAs, proposed as serum-derived biomarkers [58] , e.g., hsa-miR7-5p, hsa-miR370-3p and hsa-miR221-5p. MicroRNAs are small RNA molecules with a suppressive function on their target mRNAs. Their expression is regulated by developmental as well as external factors. Their extracellular presence in serum has been described previously, although the detailed function of the most circulating miRNAs is still elusive. However, their presence and abundance are used as biomarkers, associated with disease responses [12, [14] [15] [16] [17] 59, 60] . Here we analyzed the abundance of the circ-miRNAs and compared mild versus severe COVD-19 cases in cell free plasma. As shown in Table 3 and Supplementary Figure S4 (modified after [43] ), the amount of several circ-miRNAs associated with inflammation responses are changed in severe COVID-19 cases in contrast to the mild ones. This refers to mir-4516, miR-320 and miR-145b-5p, for example. This strong response at the miRNA level contributing to the immune response regulation agrees with the clinical description of patients requiring ICU admission [61, 62] . This involves higher concentrations of GCSF, IP10, MCP1, MIP1A and TNFα, suggesting that the cytokine storm was associated with disease severity. Based on an in silico analysis [63] , miR362-5p (3.6 log2fold change, p = 0.009, Table 2 ) was proposed for targeting the RNA encoding the virus receptor ACE2 and let7F2, targeting the RNA encoding trans-membrane-protein TMPRSS2. Both proteins were found to be involved in the SARS-CoV-2 infection processes. While miR362-5p abundance is increased in ICU cases, let-7F is reduced (Supplementary Figure S5) . No overlap of differential abundant miRNAs from blood tissue is found [64] . Furthermore, several symptoms of COVID-19 might be correlated with the differential abundance of some circulating miRNAS, as the subfamily of hsa-miRNA320 was found to be increased in severe cases ( Figure 3 ). The complete miRNA320 family is used as biomarker for lung-associated diseases [65] , and is down-regulated upon adenovirus infection [66] . A low abundance of hsa-miRNA320b was described as a predictive biomarker for survival rate of COPD patients [67] . Therefore, the SARS-CoV-2 induced deregulation of the miRNA320 family is involved in the infection response, reflecting the molecular status of the lung tissue. hsa-miRNA625-3p (Figure 3) , a circ-miRNA used as biomarker for Malignant Pleural Mesothelioma [39] (up-regulated) and Nontuberculous Mycobacterial Pulmonary Disease [40] (up-regulated) was found to be strongly down-regulated (−1.64 log2fold change, p = 0.0468) in all severe cases. Its counterpart, hsa-miRNA625-5p, was reported to suppress the inflammatory responses by targeting AKT2 in human bronchial epithelial cells [41] . Overexpression of hsa-miR-502-3p (here: 1.7 log2fold) was described to be involved in the regulation of the STAT3 pathway [68] . hsa-miR-363-3p (here: −2.5 log2fold) was described as a biomarker for an acute coronary syndrome [69] . In addition, miRNA hsa-miR-126-3p was found to be reduced in severe cases. This finding was recently confirmed by qPCR analysis, comparing healthy persons and COVID-19 patients [70] . hsa-miR-126-3p (here: −2.1 log2fold) has been shown to control the DNA damage response by repressing ATM protein kinase activity in endothelial cells, and is also a possible contributor to endothelial and lung tissue dysfunction [71] . Our data give a first insight into a novel field of SARS-CoV-2 virus defense. We are aware that the cohort size was very limited. Within our data sample set, we are dealing with a sex bias and have only limited information about characteristics and other predispositions that may affect the levels of circulating microRNAs; therefore, the presented results should be seen as a first insight, and further conclusions might require additional experimental pursuits to prove. Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/v13081593/s1, Figure S1 : Illustration for Introduction, Figure S2 : Illustration of Sampling, antibody status and SARS-CoV-2 derived small RNAs in patients (pt), Figure S3 : Analysis of circhsa-miRNAs, Figure S4 : Selected hsa-miRNAs associated with inflammation response, Figure S5 : Abundance of selected has-miRNAs, Table S1 : SARS-CoV-2 derived reads, Table S2 : circ-miRNAs. Institutional Review Board Statement: The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board approvals numbers 2020-039 and 11/17. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study. Data Availability Statement: The original reads sequences are available for download at: https:// doi.ipk-gatersleben.de/DOI/42ac1f14-ab8e-440f-bd7b-0faa465cfbce/39046ce0-3b02-4428-9083-f923 60d0b2cb/2/1847940088. DOI of the dataset released in this manuscript were constructed using the e!DAL system [72] . COVID-19: Immunology and treatment options Investigating the impact of influenza on excess mortality in all ages in Italy during recent seasons (2013/14-2016/17 seasons) A novel coronavirus emerging in China-Key questions for impact assessment COVID-19 mortality risk in Down syndrome: Results from a cohort study of 8 million adults RNA interference: Antiviral defense and genetic tool Origins and evolution of eukaryotic RNA interference Virus-encoded microRNAs Coronavirus nonstructural protein 15 mediates evasion of dsRNA sensors and limits apoptosis in macrophages RNAi in the regulation of mammalian viral infections Using siRNA in prophylactic and therapeutic regimens against SARS coronavirus in Rhesus macaque Small interfering RNA effectively inhibits the expression of SARS coronavirus membrane gene at two novel targeting sites Comparing the microRNA spectrum between serum and plasma A virus-derived microRNA-like small RNA serves as a serum biomarker to prioritize the COVID-19 patients at high risk of developing severe disease Serum miR-92 a-3 p as a new potential biomarker for diagnosis of Kawasaki disease with coronary artery lesions Circulating microRNAs as diagnostic biomarkers for cardiovascular diseases Extracellular microRNAs in human circulation are associated with miRISC complexes that are accessible to anti-AGO2 antibody and can bind target mimic oligonucleotides Large differences in small RNA composition between human biofluids The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak Methods of inactivation of SARS-CoV-2 for downstream biological assays Non-invasive approach for evaluation of pulmonary hypertension using extracellular vesicle-associated small non-coding RNA Next-generation sequencing of T and B cell receptor repertoires from COVID-19 patients showed signatures associated with severity of disease Diagnosis of severe community-acquired pneumonia caused by Acinetobacter baumannii through next-generation sequencing: A case report The sequence structures of human microRNA molecules and their implications A plant small RNA target analysis server (2017 release) Taste perception: From the tongue to the testis Critical role of GFR1 in the development and function of the main olfactory system Genome-wide association analysis of ischemic stroke in young adults. G3 Genes Genomes Genet Polymorphisms in the 5-HTR2A gene related to obstructive sleep apnea syndrome Novel identification and characterisation of transient receptor potential melastatin 3 ion channels on natural killer cells and B lymphocytes: Effects on cell signalling in chronic fatigue syndrome/myalgic encephalomyelitis patients Intussusceptive angiogenesis in Covid-19: Hypothesis on the significance and focus on the possible role of FGF2 Placental growth factor level in plasma predicts COVID-19 severity and in-hospital mortality Schlafen 14 (SLFN14) is a novel antiviral factor involved in the control of viral replication From microRNA sequences to function Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 Nuclear miR-320 mediates diabetes-induced cardiac dysfunction by activating transcription of fatty acid metabolic genes to cause lipotoxicity in the heart Downregulation of miR-320 alleviates endoplasmic reticulum stress and inflammatory response in 3T3-L1 adipocytes MiR-320 regulates cardiomyocyte apoptosis induced by ischemia-reperfusion injury by targeting AKIP1 MicroRNA hsa-miR-29a-3p is a plasma biomarker for the differential diagnosis and monitoring of tuberculosis Increased Circulating miR-625-3p: A potential biomarker for patients with malignant pleural mesothelioma miRNA expression profiles and potential as biomarkers in nontuberculous mycobacterial pulmonary disease miR-625-5p suppresses inflammatory responses by targeting AKT2 in human bronchial epithelial cells Furin expression in patients with psoriasis-A patient cohort endangered to SARS-COV2? Front A guide to miRNAs in inflammation and innate immune responses miRCOVID-19: Potential targets of human miRNAs in SARS-CoV-2 for RNA-based drug discovery Methods for detection and study of virus-derived small RNAs produced from the intramolecular base-pairing region of the picornavirus genome Influenza A virus-derived siRNAs increase in the absence of NS1 yet fail to inhibit virus replication No evidence for viral small RNA production and antiviral function of Argonaute 2 in human cells Novel SARS-CoV-2 encoded small RNAs in the passage to humans RNA virus microRNA that mimics a B-cell oncomiR Autonomic dysfunction in "long COVID": Rationale, physiology and management strategies Nextstrain: Real-time tracking of pathogen evolution Protocol: A highly sensitive RT-PCR method for detection and quantification of microRNAs Recurrence of SARS-CoV-2 viral RNA in recovered COVID-19 patients: A narrative review SARS-CoV-2 antibody responses correlate with resolution of RNAemia but are short-lived in patients with mild illness Viral RNA level, serum antibody responses, and transmission risk in recovered COVID-19 patients with recurrent positive SARS-CoV-2 RNA test results: A population-based observational cohort study Detection of novel coronavirus (SARS-CoV-2) RNA in peripheral blood specimens Serum miRNA as a possible biomarker in the diagnosis of bipolar II disorder Serum microRNA expression profile as a biomarker for the diagnosis of pertussis miRNAs can be generally associated with human pathologies as exemplified for miR-144 Clinical features of patients infected with 2019 novel coronavirus in Altered cytokine levels and immune responses in patients with SARS-CoV-2 infection and related conditions Identification of novel hypothalamic microRNAs as promising therapeutics for SARS-CoV-2 by regulating ACE2 and TMPRSS2 expression: An in silico analysis The noncoding and coding transcriptional landscape of the peripheral immune response in patients with COVID-19 Correlation of plasma exosomal microRNAs with the efficacy of immunotherapy in EGFR/ALK wild-type advanced non-small cell lung cancer Fluctuating expression of microRNAs in adenovirus infected cells Low miR-150-5p and miR-320b expression predicts reduced survival of COPD patients Overexpressed circ-RPL15 predicts poor survival and promotes the progression of gastric cancer via regulating miR-502-3p/OLFM4/STAT3 pathway miR-941 as a promising biomarker for acute coronary syndrome Circulating cardiovascular microRNAs in critically ill COVID-19 patients Downregulation of microRNA-126 augments DNA damage response in cigarette smokers and patients with chronic obstructive pulmonary disease ! DAL-A framework to store, share and publish research data We thank Sandra Driesslein (Leibniz Institute of Plant Genetics and Crop Plant Research, IPK Gatersleben, Germany) for her excellent technical assistance and Helmut Bäumlein and Rhonda Meyer for their critical reading and constructive comments on the manuscript. We thank Christoph Wosiek for excellent technical assistance. We thank all patients and donors for donating their blood for this project. The authors have no conflict of interests.