Small sample sizes in high-throughput miRNA screens: A common pitfall for the identification of miRNA biomarkers

Circulating extracellular vesicles and small non-coding RNAs cargo in idiopathic inflammatory myopathies reveal differences across myositis subsets

OBJECTIVE

To investigate the epigenetic footprint of idiopathic inflammatory myopathies (IIM) through characterization of circulating extracellular vesicles (EVs) and the expression of EV-derived small non-coding RNAs (sncRNAs).

METHODS

In this cross-sectional study, EVs were isolated by size-exclusion chromatography from plasma of patients with IIM and age- and sex-matched healthy donors (HD). EV-derived sncRNAs were sequenced and quantified using Next-Generation Sequencing (NGS). Following quality control and normalization, filtered count reads were used for differential microRNA (miRNA) and piwi-interacting RNA (piRNA) expression analyses. Putative gene targets enriched for pathways implicated in IIM were analyzed. Patients' clinical and laboratory characteristics at the time of sampling were recorded.

RESULTS

Forty-seven IIM patients and 45 HD were enrolled. MiR-486-5p (p < 0.01), miR-122-5p, miR-192-5p, and miR-32-5p were significantly upregulated (p < 0.05 for all), while miR-142-3p (p < 0.001), miR-141-3p (p < 0.01), let-7a-5p (p < 0.05) and miR-3613-5p (p < 0.05) downregulated in EVs from IIM patients versus HD. MiR-486-5p was associated with raised muscle enzymes levels. Several target genes of up/downregulated miRNAs in IIM participate in inflammation, necroptosis, interferon and immune signaling. Six piRNAs were significantly dysregulated in IIM EVs versus HD (p < 0.05). Within IIM, miR-335-5p was selectively upregulated and miR-27a-5p downregulated in dermatomyositis (n = 21, p < 0.01). Finally, plasma EV levels were significantly increased in cancer-associated myositis (CAM, n = 12) versus non-CAM IIM (n = 35, p = 0.02) and HD (p < 0.01). EVs cargo in CAM was significantly enriched of let-7f-5p and depleted of miR-143-3p.

CONCLUSION

Through an unbiased screening of EV-derived sncRNAs, we characterize miRNAs and piRNAs in the EVs cargo as potential biomarkers and modifiers of diverse IIM phenotypes.

Identification of Differentially Expressed Genes and microRNAs in the Gray and White Feather Follicles of Shitou Geese

The Shitou goose, a highly recognized indigenous breed with gray plumage originating from Chaozhou Raoping in Guangdong Province, China, is renowned for being the largest goose species in the country. Notably, during the pure breeding process of Shitou geese, approximately 2% of the offspring in each generation unexpectedly exhibited white plumage. To better understand the mechanisms underlying white plumage color formation in Shitou geese, we conducted a comparative transcriptome analysis between white and gray feather follicles, aiming to identify key genes and microRNAs that potentially regulate white plumage coloration in this unique goose breed. Our results revealed a number of pigmentation genes, encompassing TYR, TYRP1, EDNRB2, MLANA, SOX10, SLC45A2, GPR143, TRPM1, OCA2, ASIP, KIT, and SLC24A5, which were significantly down-regulated in the white feather follicles of Shitou geese. Among these genes, EDNRB2 and KIT emerged as the most promising candidate genes for white plumage coloration in Shitou geese. Additionally, our analysis also uncovered 46 differentially expressed miRNAs. Of these, miR-144-y may play crucial roles in the regulation of feather pigmentation. Furthermore, the expression of novel-m0086-5p, miR-489-y, miR-223-x, miR-7565-z, and miR-3535-z exhibits a significant negative correlation with the expression of pigmentation genes including TYRP1, EDNRB2, MLANA, SOX10, TRPM1, and KIT, suggesting these miRNAs may indirectly regulate the expression of these genes, thereby influencing feather color. Our findings provide valuable insights into the genetic mechanisms underlying white plumage coloration in Shitou geese and contribute to the broader understanding of avian genetics and coloration research.

Non-Coding RNA as Biomarkers and Their Role in the Pathogenesis of Gastric Cancer-A Narrative Review

Non-coding RNAs (ncRNAs) represent a broad family of molecules that regulate gene expression, including microRNAs, long non-coding RNAs and circular RNAs, amongst others. Dysregulated expression of ncRNAs alters gene expression, which is implicated in the pathogenesis of several malignancies and inflammatory diseases. Gastric cancer is the fifth most frequently diagnosed cancer and the fourth most common cause of cancer-related death. Studies have found that altered expression of ncRNAs may contribute to tumourigenesis through regulating proliferation, apoptosis, drug resistance and metastasis. This review describes the potential use of ncRNAs as diagnostic and prognostic biomarkers. Moreover, we discuss the involvement of ncRNAs in the pathogenesis of gastric cancer, including their interactions with the members of major signalling pathways.

A combined circulating microRNA panel predicts the risk of vascular calcification in community-dwelling older adults with age strata differences

BACKGROUND

Older adults have a higher risk of developing vascular calcification (VC). Circulating miRNAs can be potential risk indicators. However, prior studies used single miRNA mostly, whereas miRNA panels were rarely evaluated. We aimed to examine whether a miRNA panel outperformed each miRNA alone, and analyzed whether advanced age affected VC risk predictive performance offered by the miRNA panel.

METHODS

We prospectively enrolled older adults (age ≥65 years) during their annual health checkup in 2017, and examined their VC severity followed by analyzing sera for VC regulatory miRNAs (miR-125b-5p, miR-125b-3p, and miR-378a-3p). We used multiple regression analyses to determine associations between each miRNA or a 3-combind panel and VC risk, followed by area under the receiver-operating-characteristics curve (AUROC) analysis. Participants were further divided to those of 65-75 and ≥75 years for comparison.

RESULTS

From 199 older adults screened, 169 (median age, 73.3 years) with available calcification assessment were analyzed, among whom 74.6 % having VC. Those with VC had significantly lower circulating miR-125b-5p, miR-125b-3p, and miR-378a-3p levels than those without. Regression analyses showed that the 3-combined miRNA panel exhibited significant associations with VC risk, with significantly higher AUROC than those of models based on individual miRNA. Importantly, in those ≥75 years, the miRNA-predicted risk of VC was more prominent than that in the 65-75 years group.

CONCLUSION

A miRNA panel for VC risk prediction might outperform individual miRNA alone in older adults, and advanced age modified the association between circulating miRNAs and the risk of VC.

No Detectable Differences in microRNA Plasma Levels between Diabetic Hypertensive Patients with and without Incident Subclinical Atrial Fibrillation

Background: Long-term rhythm monitoring (LTRM) can detect undiagnosed atrial fibrillation (AF) in patients at risk of AF and stroke. Circulating microRNAs (miRNAs), which have been shown to play a role in atrial electrical and structural remodelling, could help to select patients who would benefit most from LTRM. The aim of this study was to investigate whether patients with diabetes mellitus (DM) and hypertension and screen-detected subclinical AF (SCAF) using an insertable cardiac monitor (ICM) have significantly different plasma baseline levels of five selected miRNAs playing a role in the modulation of atrial electrical and structural remodelling (miR-21-5p, miR-29b-3p, miR-150-5p, miR-328-3p, and miR-432-5p) compared to those without SCAF. Methods: This study was performed at the outpatient clinic of a secondary academic teaching hospital between December 2013 and November 2015. Eligible patients were ≥65 years of age with DM and hypertension but without known heart diseases. All patients received an ICM. On the day of ICM implantation, blood samples for the measurement of plasma levels of the five miRNAs were drawn. In this post hoc analysis, we investigated their expression by reverse transcription-quantitative polymerase chain reaction. MiRNA plasma levels in patients with and without newly detected SCAF were compared. Results: We included 82 consecutive patients (median age of 71.3 years (IQR 67.4-75.1)), who were followed for a median of 588 days (IQR: 453-712 days). Seventeen patients (20.7%) had ICM-detected SCAF. Plasma levels of miR-328-3p, miR-29b-3p, miR-21-5p, miR-432-5p, and miR-150-5p were slightly but not significantly different in patients with incident SCAF compared with patients without. Conclusions: In patients with hypertension and DM, newly detected SCAF was not significantly associated with changes in expression levels of miR-21-5p, miR-29b-3p, miR-150-5p, miR-328-3p, and miR-432-5p.

The Diagnostic Value of microRNA Expression Analysis in Detecting Intraductal Papillomas in Patients with Pathological Nipple Discharge

Patients with pathological nipple discharge (PND) often undergo local surgical procedures because standard radiologic imaging fails to identify the underlying cause. MicroRNA (MiRNA) expression analysis of nipple fluid holds potential for distinguishing between breast diseases. This study aimed to compare miRNA expression levels between nipple fluids from patients with PND to identify possible relevant miRNAs that could differentiate between intraductal papillomas and no abnormalities in the breast tissue. Nipple fluid samples from patients with PND without radiological and pathological suspicion for malignancy who underwent a ductoscopy procedure were analyzed. We used univariate and multivariate regression analyses to identify nipple fluid miRNAs differing between pathologically confirmed papillomas and breast tissue without abnormalities. A total of 27 nipple fluid samples from patients with PND were included for miRNA expression analysis. Out of the 22 miRNAs examined, only miR-145-5p was significantly differentially expressed (upregulated) in nipple fluid from patients with an intraductal papilloma compared to patients showing no breast abnormalities (OR 4.76, p = 0.046), with a diagnostic accuracy of 92%. miR-145-5p expression in nipple fluid differs for intraductal papillomas and breast tissue without abnormalities and, therefore, has potential as a diagnostic marker to signal presence of papillomas in PND patients. However, further refinement and validation in clinical trials are necessary to establish its clinical applicability.

MicroRNA signature from extracellular vesicles of HCV/HIV co-infected individuals differs from HCV mono-infected

Hepatitis C virus (HCV) coinfection with human immunodeficiency virus (HIV) has a detrimental impact on disease progression. Increasing evidence points to extracellular vesicles (EVs) as important players of the host-viral cross-talk. The microRNAs (miRNAs), as essential components of EVs cargo, are key regulators of normal cellular processes and also promote viral replication, viral pathogenesis, and disease progression. We aimed to characterize the plasma-derived EVs miRNA signature of chronic HCV infected and HIV coinfected patients to unravel the molecular mechanisms of coinfection. EVs were purified and characterized from 50 plasma samples (21 HCV mono- and 29 HCV/HIV co-infected). EV-derived small RNAs were isolated and analyzed by massive sequencing. Known and de novo miRNAs were identified with miRDeep2. Significant differentially expressed (SDE) miRNA identification was performed with generalized linear models and their putative dysregulated biological pathways were evaluated. Study groups were similar for most clinical and epidemiological characteristics. No differences were observed in EVs size or concentration between groups. Therefore, HCV/HIV co-infection condition did not affect the concentration or size of EVs but produced a disturbance in plasma-derived EVs miRNA cargo. Thus, a total of 149 miRNAs were identified (143 known and 6 de novo) leading to 37 SDE miRNAs of which 15 were upregulated and 22 downregulated in HCV/HIV co-infected patients. SDE miRNAs regulate genes involved in inflammation, fibrosis, and cancer, modulating different biological pathways related to HCV and HIV pathogenesis. These findings may help to develop new generation biomarkers and treatment strategies, in addition to elucidate the mechanisms underlying virus-host interaction. KEY MESSAGES: HCV and HCV/HIV displayed similar plasma-EV size and concentration. EVs- derived miRNA profile was characterized by NGS. 37 SDE miRNAs between HCV and HCV/HIV were observed. SDE miRNAs regulate genes involved in inflammation, fibrosis and cancer.

The cancer glycocode as a family of diagnostic biomarkers, exemplified by tumor-associated gangliosides

One unexploited family of cancer biomarkers comprise glycoproteins, carbohydrates, and glycolipids (the Tumor Glycocode).A class of glycolipid cancer biomarkers, the tumor-marker gangliosides (TMGs) are presented here as potential diagnostics for detecting cancer, especially at early stages, as the biological function of TMGs makes them etiological. We propose that a quantitative matrix of the Cancer Biomarker Glycocode and artificial intelligence-driven algorithms will expand the menu of validated cancer biomarkers as a step to resolve some of the challenges in cancer diagnosis, and yield a combination that can identify a specific cancer, in a tissue-agnostic manner especially at early stages, to enable early intervention. Diagnosis is critical to reducing cancer mortality but many cancers lack efficient and effective diagnostic tests, especially for early stage disease. Ideal diagnostic biomarkers are etiological, samples are preferably obtained via non-invasive methods (e.g. liquid biopsy of blood or urine), and are quantitated using assays that yield high diagnostic sensitivity and specificity for efficient diagnosis, prognosis, or predicting response to therapy. Validated biomarkers with these features are rare. While the advent of proteomics and genomics has led to the identification of a multitude of proteins and nucleic acid sequences as cancer biomarkers, relatively few have been approved for clinical use. The use of multiplex arrays and artificial intelligence-driven algorithms offer the option of combining data of known biomarkers; however, for most, the sensitivity and the specificity are below acceptable criteria, and clinical validation has proven difficult. One strategic solution to this problem is to expand the biomarker families beyond those currently exploited. One unexploited family of cancer biomarkers comprise glycoproteins, carbohydrates, and glycolipids (the Tumor Glycocode). Here, we focus on a family of glycolipid cancer biomarkers, the tumor-marker gangliosides (TMGs). We discuss the diagnostic potential of TMGs for detecting cancer, especially at early stages. We include prior studies from the literature to summarize findings for ganglioside quantification, expression, detection, and biological function and its role in various cancers. We highlight the examples of TMGs exhibiting ideal properties of cancer diagnostic biomarkers, and the application of GD2 and GD3 for diagnosis of early stage cancers with high sensitivity and specificity. We propose that a quantitative matrix of the Cancer Biomarker Glycocode and artificial intelligence-driven algorithms will expand the menu of validated cancer biomarkers as a step to resolve some of the challenges in cancer diagnosis, and yield a combination that can identify a specific cancer, in a tissue-agnostic manner especially at early stages, to enable early intervention.

Blood transcriptomic signature in type-2 biomarker-low severe asthma and asthma control

BACKGROUND

Patients with type-2 (T2) cytokine-low severe asthma often have persistent symptoms despite suppression of T2 inflammation with corticosteroids.

OBJECTIVES

We sought to analyze whole blood transcriptome from 738 samples in T2-biomarker-high/-low patients with severe asthma to relate transcriptomic signatures to T2 biomarkers and asthma symptom scores.

METHODS

Bulk RNA-seq data were generated for blood samples (baseline, week 24, week 48) from 301 participants recruited to a randomized clinical trial of corticosteroid optimization in severe asthma. Unsupervised clustering, differential gene expression analysis, and pathway analysis were performed. Patients were grouped by T2-biomarker status and symptoms. Associations between clinical characteristics and differentially expressed genes (DEGs) associated with biomarker and symptom levels were investigated.

RESULTS

Unsupervised clustering identified 2 clusters; cluster 2 patients were blood eosinophil-low/symptom-high and more likely to be receiving oral corticosteroids (OCSs). Differential gene expression analysis of these clusters, with and without stratification for OCSs, identified 2960 and 4162 DEGs, respectively. Six hundred twenty-seven of 2960 genes remained after adjusting for OCSs by subtracting OCS signature genes. Pathway analysis identified dolichyl-diphosphooligosaccharide biosynthesis and assembly of RNA polymerase I complex as significantly enriched pathways. No stable DEGs were associated with high symptoms in T2-biomarker-low patients, but numerous associated with elevated T2 biomarkers, including 15 that were upregulated at all time points irrespective of symptom level.

CONCLUSIONS

OCSs have a considerable effect on whole blood transcriptome. Differential gene expression analysis demonstrates a clear T2-biomarker transcriptomic signature, but no signature was found in association with T2-biomarker-low patients, including those with a high symptom burden.

Identification of hsa_circ_0001445 of a novel circRNA-miRNA-mRNA regulatory network as potential biomarker for coronary heart disease

Objects

To evaluate the hsa_circ_0001445 level in peripheral blood leukocytes of patients with coronary heart disease (CHD) and its related clinical factors, and predict its circRNA-miRNA-mRNA regulatory network in CHD pathogenesis via bioinformatics analysis.

Methods

Peripheral blood leukocytes were isolated from the whole blood samples of 94 CHD patients (aged 65.96 ± 9.78 years old) and 126 healthy controls (aged 60.75 ± 8.81 years old). qRT-PCR was used to quantify the expression level of circRNA and subsequently analyze its association with CHD clinical parameters. Via bioinformatics algorithm and GEO datasets, differential miRNA expression was evaluated using the Limma package. A miRNA-mRNA regulatory network was predicted by cyTargetLinker. ClusterProfiler was employed to perform functional enrichment analysis of the circRNA network to investigate its role in CHD pathogenesis.

Results

The expression of hsa_circ_0001445 in peripheral blood leukocytes of CHD patients was downregulated compared with that of healthy controls. Positive correlations were evident between hsa_circ_0001445 expression level and the levels of hemoglobin, triglycerides, high- and low-density lipoprotein cholesterol. A significant negative correlation was also found between hsa_circ_0001445 expression level and age and the neutrophil level. Low expression of hsa_circ_0001445 exhibited a discriminatory ability between CHD patients and healthy controls with a sensitivity of 67.5% and a specificity of 76.6% (p &lt; 0.05). By bioinformatics analysis, 405 gene ontology terms were identified. The Kyoto Encyclopedia of Genes and Genomes terms focused principally on the PI3K-Akt signaling pathway. hsa_circ_0001445 was associated with the expression of three miRNAs that may regulate 18 genes involved in KEGG processes: hsa-miR-507, hsa-miR-375–3p, and hsa-miR-942–5p.

Conclusion

The hsa_circ_0001445 level in peripheral blood leukocytes may serve as a biomarker for CHD diagnosis. Our work on circRNA-miRNA-mRNA networks suggests a potential role for hsa_circ_0001445 in CHD development.

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

OBJECTIVE

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).

RESEARCH DESIGN AND METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

Development of circulating microRNA-based biomarkers for medical decision-making: a friendly reminder of what should NOT be done

Circulating cell-free microRNAs (miRNAs) represent a major reservoir for biomarker discovery. Unfortunately, their implementation in clinical practice is limited due to a profound lack of reproducibility. The great technical variability linked to major pre-analytical and analytical caveats makes the interpretation of circulating cell-free miRNA data challenging and leads to inconsistent findings. Additional efforts directed to standardization are fundamental. Several well-established protocols are currently used by independent groups worldwide. Nonetheless, there are some specific aspects in specimen collection and processing, sample handling, miRNA quantification, and data analysis that should be considered to ensure reproducibility of results. Here, we have addressed this challenge using an alternative approach. We have highlighted and discussed common pitfalls that negatively impact the robustness of circulating miRNA quantification and their application for clinical decision-making. Furthermore, we provide a checklist usable by investigators to facilitate and ensure the control of the whole miRNA quantification and analytical process. We expect that these recommendations improve the reproducibility of findings, and ultimately, facilitate the incorporation of circulating miRNA profiles into clinical practice as the next generation of disease biomarkers.

Potential of Circulating miRNAs as Molecular Markers in Mood Disorders and Associated Suicidal Behavior

Mood disorders are the most prevalent psychiatric disorders associated with significant disability, morbidity, and mortality. The risk of suicide is associated with severe or mixed depressive episodes in patients with mood disorders. However, the risk of suicide increases with the severity of depressive episodes and is often presented with higher incidences in bipolar disorder (BD) patients than in patients with major depression (MDD). Biomarker study in neuropsychiatric disorders is critical for developing better treatment plans by facilitating more accurate diagnosis. At the same time, biomarker discovery also provides more objectivity to develop state-of-the-art personalized medicine with increased accuracy through clinical interventions. Recently, colinear changes in miRNA expression between brain and systemic circulation have added great interest in examining their potential as molecular markers in mental disorders, including MDD, BD, and suicidality. A present understanding of circulating miRNAs in body fluids implicates their role in managing neuropsychiatric conditions. Most notably, their use as prognostic and diagnostic markers and their potential role in treatment response have significantly advanced our knowledge base. The present review discusses circulatory miRNAs and their underlying possibilities to be used as a screening tool for assessing major psychiatric conditions, including MDD, BD, and suicidal behavior.

A study of differential microRNA expression profile in migraine: the microMIG exploratory study

BACKGROUND

Several studies have described potential microRNA (miRNA) biomarkers associated with migraine, but studies are scarcely reproducible primarily due to the heterogeneous variability of participants. Increasing evidence shows that disease-related intrinsic factors together with lifestyle (environmental factors), influence epigenetic mechanisms and in turn, diseases. Hence, the main objective of this exploratory study was to find differentially expressed miRNAs (DE miRNA) in peripheral blood mononuclear cells (PBMC) of patients with migraine compared to healthy controls in a well-controlled homogeneous cohort of non-menopausal women.

METHODS

Patients diagnosed with migraine according to the International Classification of Headache Disorders (ICHD-3) and healthy controls without familial history of headache disorders were recruited. All participants completed a very thorough questionnaire and structured-interview in order to control for environmental factors. RNA was extracted from PBMC and a microarray system (GeneChip miRNA 4.1 Array chip, Affymetrix) was used to determine the miRNA profiles between study groups. Principal components analysis and hierarchical clustering analysis were performed to study samples distribution and random forest (RF) algorithms were computed for the classification task. To evaluate the stability of the results and the prediction error rate, a bootstrap (.632 + rule) was run through all the procedure. Finally, a functional enrichment analysis of selected targets was computed through protein-protein interaction networks.

RESULTS

After RF classification, three DE miRNA distinguished study groups in a very homogeneous female cohort, controlled by factors such as demographics (age and BMI), life-habits (physical activity, caffeine and alcohol consumptions), comorbidities and clinical features associated to the disease: miR-342-3p, miR-532-3p and miR-758-5p. Sixty-eight target genes were predicted which were linked mainly to enriched ion channels and signaling pathways, neurotransmitter and hormone homeostasis, infectious diseases and circadian entrainment.

CONCLUSIONS

A 3-miRNA (miR-342-3p, miR-532-3p and miR-758-5p) novel signature has been found differentially expressed between controls and patients with migraine. Enrichment analysis showed that these pathways are closely associated with known migraine pathophysiology, which could lead to the first reliable epigenetic biomarker set. Further studies should be performed to validate these findings in a larger and more heterogeneous sample.

Missing data in multi-omics integration: Recent advances through artificial intelligence

Biological systems function through complex interactions between various 'omics (biomolecules), and a more complete understanding of these systems is only possible through an integrated, multi-omic perspective. This has presented the need for the development of integration approaches that are able to capture the complex, often non-linear, interactions that define these biological systems and are adapted to the challenges of combining the heterogenous data across 'omic views. A principal challenge to multi-omic integration is missing data because all biomolecules are not measured in all samples. Due to either cost, instrument sensitivity, or other experimental factors, data for a biological sample may be missing for one or more 'omic techologies. Recent methodological developments in artificial intelligence and statistical learning have greatly facilitated the analyses of multi-omics data, however many of these techniques assume access to completely observed data. A subset of these methods incorporate mechanisms for handling partially observed samples, and these methods are the focus of this review. We describe recently developed approaches, noting their primary use cases and highlighting each method's approach to handling missing data. We additionally provide an overview of the more traditional missing data workflows and their limitations; and we discuss potential avenues for further developments as well as how the missing data issue and its current solutions may generalize beyond the multi-omics context.

Evaluation of altered miRNA expression pattern to predict COVID-19 severity

Outbreak of COVID-19 pandemic in December 2019 affected millions of people globally. After substantial research, several biomarkers for COVID-19 have been validated however no specific and reliable biomarker for the prognosis of patients with COVID-19 infection exists. Present study was designed to identify specific biomarkers to predict COVID-19 severity and tool for formulating treatment. A small cohort of subjects (n = 43) were enrolled and categorized in four study groups; Dead (n = 16), Severe (n = 10) and Moderate (n = 7) patients and healthy controls (n = 10). Small RNA sequencing was done on Illumina platform after isolation of microRNA from peripheral blood. Differential expression (DE) of miRNA (patients groups compared to control) revealed 118 down-regulated and 103 up-regulated known miRNAs with fold change (FC) expression ≥2 folds and p ≤ 0.05. DE miRNAs were then subjected to functional enrichment and network analysis. Bioinformatic analysis resulted in 31 miRNAs (24 Down-regulated; 7 up-regulated) significantly associated with COVID-19 having AUC>0.8 obtained from ROC curve. Seventeen out of 31 DE miRNAs have been linked to COVID-19 in previous studies. Three miRNAs, hsa-miR-147b-5p and hsa-miR-107 (down-regulated) and hsa-miR-1299 (up-regulated) showed significant unique DE in Dead patients. Another set of 4 miRNAs, hsa-miR-224-5p (down-regulated) and hsa-miR-4659b-3p, hsa-miR-495-3p and hsa-miR-335-3p were differentially up-regulated uniquely in Severe patients. Members of three miRNA families, hsa-miR-20, hsa-miR-32 and hsa-miR-548 were significantly down-regulated in all patients group in comparison to healthy controls. Thus a distinct miRNA expression profile was observed in Dead, Severe and Moderate COVID-19 patients. Present study suggests a panel of miRNAs which identified in COVID-19 patients and could be utilized as potential diagnostic biomarkers for predicting COVID-19 severity.

Key methodological challenges in detecting circulating miRNAs in different biofluids

The technological advancement in diagnostic techniques has immensely improved the capability of predicting disease progression. Yet, there is a great interest in developing newer biomarkers that can enhance disease risk prediction thereby minimising the associated morbidity and mortality. Circulating miRNAs, a non-coding RNA molecule, are critical regulators in the pathophysiology of various complex multifactorial diseases. In recent years, circulating miRNAs have been enormously studied and are considered as an emerging biomarker due to their easy accessibility, stability, and detection by sequence-specific amplification methods. However, there is a distinct lack of consensus regarding the preanalytical factors such as preferred sample selection, methodological aspects, etc that may independently or together influence the detection of circulating miRNAs resulting in erroneous expression profiles. Therefore, the present review makes an attempt to highlight the various pre-analytical and analytical factors that can potentially influence the circulating miRNA levels. Literature on circulating miRNA's stability, processing and quantitation in different biofluids along with the effect of various controllable and uncontrollable factors influencing circulating miRNA expression have been summarised in the current review.

microRNAs as biomarkers of risk of major adverse cardiovascular events in atrial fibrillation

Atrial fibrillation is a complex and multifactorial disease. Although prophylactic anticoagulation has great benefits in avoiding comorbidities, adverse cardiovascular events still occur and thus in recent decades, many resources have been invested in the identification of useful markers in the prevention of the risk of MACE in these patients. As such, microRNAs, that are small non-coding RNAs whose function is to regulate gene expression post-transcriptionally, have a relevant role in the development of MACE. miRNAs, have been investigated for many years as potential non-invasive biomarkers of several diseases. Different studies have shown their utility in the diagnosis and prognosis of cardiovascular diseases. In particular, some studies have associated the presence of certain miRNAs in plasma with the development of MACE in AF. Despite these results, there are still many efforts to be done to allow the clinical use of miRNAs. The lack of standardization concerning the methodology in purifying and detecting miRNAs, still provides contradictory results. miRNAs also have a functional impact in MACE in AF through the dysregulation of immunothrombosis. Indeed, miRNAs may be a link between MACE and inflammation, through the regulation of neutrophil extracellular traps that are a key element in the establishment and evolution of thrombotic events. The use of miRNAs as therapy against thromboinflammatory processes should also be a future approach to avoid the occurrence of MACE in atrial fibrillation.

Dental implant material related changes in molecular signatures in peri-implantitis - A systematic review and integrative analysis of omics in-vitro studies

OBJECTIVE

Since peri-implantitis differs clinically and histopathologically from periodontitis, implant wear debris is considered to play a role in the destructive processes. This work aims to systematically review if titanium particles affect oral-related cells through changes in molecular signatures (e.g., transcriptome, proteome, epigenome), thereby promoting peri-implantitis.

METHODS

Leveraging three literature databases (Medline, Embase, Cochrane) a systematic search based on a priori defined PICOs was conducted: '-omics' studies examining titanium exposure in oral-related cells. After risk of bias assessments, lists of differentially expressed genes, proteins, and results of functional enrichment analyses were compiled. The significance of overlapping genes across multiple studies was assessed via Monte Carlo simulation and their ranking was verified using rank aggregation.

RESULTS

Out of 2104 screened articles we found 12 eligible publications. A significant overlap of gene expression in oral-related cells exposed to titanium particles was found in four studies. Furthermore, changes in biological processes like immune/inflammatory or stress response as well as toll-like receptor (TLR) and mitogen-activated protein kinase (MAPK) signaling pathways were linked to titanium in transcriptome and proteome studies. Epigenetic changes caused by titanium were detected but inconsistent.

CONCLUSION

An influence of titanium implant wear debris on the development and progression of peri-implantitis is plausible but needs to be proven in further studies. Limitations arise from small sample sizes of included studies and insufficient publication of re-analyzable data.

The potential of three whole blood microRNAs to predict outcome and monitor treatment response in sarcoid-bearing equids

MicroRNAs (miRNAs) have been proposed as biomarkers for equine sarcoid (ES) disease. In this study, the suitability of three whole blood miRNAs to diagnose ES and to predict and monitor the outcome of therapy was explored. Using reverse transcription-quantitative polymerase chain reaction (RT-qPCR), expression levels of eca-miR-127, eca-miR-379, and eca-miR-432 in whole blood of ES-affected equids before and at least one year after therapy were compared to those of unaffected control equids. Associations of age, sex, species, diagnosis, and therapy outcome with miRNA expression levels were examined using general linear models. In total, 48 ES-affected equids and 47 control equids were recruited. From the affected animals, 31 responded favorably to treatment, and 17 demonstrated a failure of therapy. None of the tested miRNAs were influenced by age. Male equids showed increased expression of eca-miR-127 compared to females and horses showed higher expression levels of eca-miR-379 and eca-miR-432 than donkeys. Eca-miR-127 was confirmed as a diagnostic discriminator between ES-affected and control equids. No difference in miRNA profiles before therapy was found when comparing ES-affected equids with success vs. failure of therapy. Eca-miR-379 and eca-miR-432 decreased over time in horses where therapy was successful, but not in those cases where it failed. Biological variables influence equine whole blood miRNA expression, which may complicate biomarker validation. While none of the tested miRNAs could predict the response to therapy in ES-affected equids and eca-miR-127 showed poor diagnostic accuracy for ES, eca-miR-379 and eca-miR-432 miRNAs might allow refinement of monitoring of success of ES therapy.

The effects of maternal cigarette smoking on cadmium and lead levels, miRNA expression and biochemical parameters across the feto-placental unit

Several miRNAs have been previously identified to be associated with cigarette smoke and/or the toxic metals cadmium (Cd) and lead (Pb). The aim of this study was to investigate the associations of maternal cigarette smoking with cadmium (Cd) and lead (Pb) levels, candidate miRNA expression and biochemical parameters across the feto-placental unit. miRNAs were isolated according to protocols provided by manufacturer from 72 healthy postpartum women using Qiagens' kits based on phenol/guanidine samples lysis and silica-membrane purification of total RNA. Candidate miRNAs (miR-1537, miR-190b, miR-16, miR-21, and miR-146a) were quantified by real-time PCR. Biochemical parameters were analyzed in plasma samples by standardized and harmonized enzymatic methods using appropriate calibrators, while CRP was determined by immunoturbidimetric method. Concentration of Cd and Pb in whole blood and placenta samples were measured by inductively coupled plasma mass spectroscopy. Cd levels in smokers were higher in all of the analyzed compartments of the feto-placental unit, Pb in maternal blood and placenta than non-smokers. Smokers also had a higher expression of miR-16 in maternal and miR-146a in cord plasma, and lower expression of miR-21 in the placenta in comparison to non-smokers. Urate concentrations in the maternal plasma of smokers were lower than this value in non-smokers. The study has demonstrated that maternal smoking was associated with toxic metals (Cd and Pb) levels, urate concentration and alteration of miRNA expression. Given that the effects of maternal smoking on miRNA expression are inadequate, all compartments of the feto-placental unit should be analyzed to obtain a complete picture. This paper is the first to report on the results of expression of cellular and circulating miRNAs simultaneously in maternal and fetal compartments and in the placenta.

Analysis of serum circulating MicroRNAs level in Malaysian patients with gestational diabetes mellitus

Gestational diabetes mellitus (GDM) is a severe global issue that requires immediate attention. MicroRNA expression abnormalities are possibly disease-specific and may contribute to GDM pathological processes. To date, there is limited data on miRNA profiling in GDM, especially that involves a longitudinal study. Here, we performed miRNA expression profiling in the entire duration of pregnancy (during pregnancy until parturition and postpartum) using a miRNA- polymerase chain reaction array (miRNA-PCRArray) and in-silico analysis to identify unique miRNAs expression and their anticipated target genes in Malay maternal serum. MiRNA expression levels and their unique potential as biomarkers were explored in this work. In GDM patients, the expression levels of hsa-miR-193a, hsa-miR-21, hsa-miR-23a, and hsa-miR-361 were significantly increased, but miR-130a was significantly downregulated. The area under the curve (AUC) and receiver operating characteristic (ROC) curve study demonstrated that hsa-miR-193a (AUC = 0.89060 ± 04,470, P = 0.0001), hsa-miR-21 (AUC = 0.89500 ± 04,411, P = 0.0001), and miR-130a (AUC = 0.6939 ± 0.05845, P = 0.0025) had potential biomarker features in GDM. In-silico analysis also revealed that KLF (Kruppel-Like family of transcription factor), ZNF25 (Zinc finger protein 25), AFF4 (ALF transcription elongation factor 4), C1orf143 (long intergenic non-protein coding RNA 2869), SRSF2 (serine and arginine rich splicing factor 2), and ZNF655 (Zinc finger protein 655) were prominent genes targeted by the common nodes of miR23a, miR130, miR193a, miR21, and miR361.Our findings suggest that circulating microRNAs in the first trimester has the potential for GDM screening in the Malay population.

Associations of plasma miRNAs with waist circumference and insulin resistance among women with polycystic ovary syndrome - Pilot study

BACKGROUND

Insulin resistance (IR) and central obesity are common in polycystic ovary syndrome (PCOS), but pathomechanisms for IR in PCOS are not established. Circulating microRNAs (miRNAs) are non-invasive biomarkers of epigenetic regulation that may contribute to the pathogenesis of IR and central adiposity in PCOS.

METHODS

We conducted a pilot study to examine associations of circulating miRNAs with IR and central adiposity among women with PCOS (n = 11) using high-throughput miRNA sequencing. We fit generalized linear models examining associations of waist circumference and HOMA-IR with plasma miRNAs. We used false discovery rate (FDR)-adjusted cutoff p < 0.1 to correct for multiple testing. We used miRDB's Gene Ontology (GO) tool to identify predicted pathways for top hits.

RESULTS

Mean age and BMI of participants were 27.9 years and 32.5 kg/m2, respectively. Lower levels of miR-1294 were associated with higher waist circumference (β = -0.10, FDR = 0.095). While no miRNAs were associated with HOMA-IR at our FDR cut off <0.1, 11 miRNAs were associated with waist circumference and 14 miRNAs with HOMA-IR at unadjusted p < 0.01, including members of the highly conserved miR-17/92 cluster and miR-1294 (β = -0.10, p < 0.001). The GO analysis of miR-1294 identified 54 overrepresented pathways, including "negative regulation of insulin receptor signaling" (FDR = 0.019), and 6 underrepresented pathways.

CONCLUSIONS

Plasma miR-1294 along with members of the miR-17/92 cluster and miRNAs involved in insulin signaling may be associated with central obesity and insulin resistance in PCOS. Larger studies among women with and without PCOS are needed to validate these findings.

MiRNA expression profiling in adenocarcinoma and squamous cell lung carcinoma reveals both common and specific deregulated microRNAs

The current study investigated the expression signatures of miRNAs in lung adenocarcinoma (LUAD) and squamous cell lung carcinoma (LUSC). miRNA profiling was performed using microarray in 12 LUAD and 12 LUSC samples and adjacent normal tissues. In LUAD, 107 miRNAs were significantly deregulated, whereas 235 miRNAs were deregulated in LUSC. Twenty-six miRNAs were common between the 2 cancer subtypes and 8 were prioritized for validation, in addition to 6 subtype-specific miRNAs. The RT-qPCR validation samples included 50 LUAD, 50 LUSC, and adjacent normal tissues. Eight miRNAs were validated in LUAD: 3 upregulated - miR-7-5p, miR-375-5p, miR-6785-3p, and 5 downregulated - miR-101-3p, miR-139-5p, miR-140-3p, miR-144-3p, miR-195-5p. Ten miRNAs were validated in the LUSC group: 3 upregulated - miR-7-5p, miR-21-3p, miR-650, and 7 downregulated - miR-95-5p, miR-140-3p, miR-144-3p, miR-195-5p, miR-375, miR-744-3p, and miR-4689-3p. Reactome pathway analysis revealed that the target genes of the deregulated miRNAs in LUAD were significantly enriched in cell cycle, membrane trafficking, gene expression processes, and EGFR signaling, while in LUSC, they were enriched in the immune system, transcriptional regulation by TP53, and FGFR signaling. This study identified distinct miRNA profiles in LUSC and LUAD, which are common and specific miRNAs that could be further investigated as biomarkers for diagnosis and prognosis.

Circulating microRNA signatures associated with disease severity and outcome in COVID-19 patients

Background

SARS-CoV-2 induces a spectrum of clinical conditions ranging from asymptomatic infection to life threatening severe disease. Host microRNAs have been involved in the cytokine storm driven by SARS-CoV-2 infection and proposed as candidate biomarkers for COVID-19.

Methods

To discover signatures of circulating miRNAs associated with COVID-19, disease severity and mortality, small RNA-sequencing was performed on serum samples collected from 89 COVID-19 patients (34 severe, 29 moderate, 26 mild) at hospital admission and from 45 healthy controls (HC). To search for possible sources of miRNAs, investigation of differentially expressed (DE) miRNAs in relevant human cell types in vitro.

Results

COVID-19 patients showed upregulation of miRNAs associated with lung disease, vascular damage and inflammation and downregulation of miRNAs that inhibit pro-inflammatory cytokines and chemokines, angiogenesis, and stress response. Compared with mild/moderate disease, patients with severe COVID-19 had a miRNA signature indicating a profound impairment of innate and adaptive immune responses, inflammation, lung fibrosis and heart failure. A subset of the DE miRNAs predicted mortality. In particular, a combination of high serum miR-22-3p and miR-21-5p, which target antiviral response genes, and low miR-224-5p and miR-155-5p, targeting pro-inflammatory factors, discriminated severe from mild/moderate COVID-19 (AUROC 0.88, 95% CI 0.80-0.95, p<0.0001), while high leukocyte count and low levels of miR-1-3p, miR-23b-3p, miR-141-3p, miR-155-5p and miR-4433b-5p predicted mortality with high sensitivity and specificity (AUROC 0.95, 95% CI 0.89-1.00, p<0.0001). In vitro experiments showed that some of the DE miRNAs were modulated directly by SARS-CoV-2 infection in permissive lung epithelial cells.

Conclusions

We discovered circulating miRNAs associated with COVID-19 severity and mortality. The identified DE miRNAs provided clues on COVID-19 pathogenesis, highlighting signatures of impaired interferon and antiviral responses, inflammation, organ damage and cardiovascular failure as associated with severe disease and death.

miR-338-5p Levels and Cigarette Smoking are Associated With Neuropathic Pain Severity in Individuals With Spinal Cord Injury: Preliminary Findings From a Genome-Wide microRNA Expression Profiling Screen

OBJECTIVE

To identify microRNA biomarkers and clinical factors associated with neuropathic pain after spinal cord injury.

DESIGN

Cross-sectional, secondary analysis of baseline data collected from ongoing clinical studies. Using a genome-wide microRNA screening approach, we studied differential microRNA expression in serum from 43 adults with spinal cord injury enrolled in ongoing clinical studies. Least squares regression was used to identify associations between microRNA expression, clinical factors, and neuropathic pain severity.

SETTING

Community-dwelling individuals with spinal cord injury.

PARTICIPANTS

Participants (N=43) were at least 18 years old with spinal cord injury, with 28 reporting neuropathic pain and 15 reporting no neuropathic pain.

INTERVENTIONS

Not applicable.

MAIN OUTCOME MEASURES

Pain presence, type, and intensity were assessed with the International Spinal Cord Injury Pain Basic Data Set. Serum microRNA normalized deep sequencing counts were quantified from blood samples. Participant demographic factors, injury characteristics, medication use, and health habits were collected via questionnaire.

RESULTS

miR-338-5p expression and history of cigarette smoking were associated with and explained 37% of the variance in neuropathic pain severity (R2=0.37, F2,18=5.31, P=.02) independent of other clinical factors. No association was identified between miR-338-5p levels and nociceptive pain severity.

CONCLUSIONS

Our findings suggest that miR-338-5p and cigarette smoking may both play a role in the development or maintenance of neuropathic pain after spinal cord injury. While additional work is needed to confirm these findings, validated target analysis suggests a neuroprotective role of miR-338-5p in modulating neuroinflammation and neuronal apoptosis and that its downregulation may result in maladaptive neuroplastic mechanisms contributing to neuropathic pain after spinal cord injury.

Circulating microRNAs differentiate fast-progressing from slow-progressing and non-progressing knee osteoarthritis in the Osteoarthritis Initiative cohort

Introduction:

The objective of this study is to identify circulating microRNAs that distinguish fast-progressing radiographic knee osteoarthritis (OA) in the Osteoarthritis Initiative cohort by applying microRNA-sequencing.

Methods:

Participants with Kellgren–Lawrence (KL) grade 0/1 at baseline were included (N = 106). Fast-progressors were defined by an increase to KL 3/4 by 4-year follow-up (N = 20), whereas slow-progressors showed an increase to KL 2/3/4 only at 8-year follow-up (N = 35). Non-progressors remained at KL 0/1 by 8-year follow-up (N = 51). MicroRNA-sequencing was performed on plasma collected at baseline and 4-year follow-up from the same participants. Negative binomial models were fitted to identify differentially expressed (DE) microRNAs. Penalized logistic regression (PLR) analyses were performed to select combinations of DE microRNAs that distinguished fast-progressors. Area under the receiver operating characteristic curves (AUC) were constructed to evaluate predictive ability.

Results:

DE analyses revealed 48 microRNAs at baseline and 2 microRNAs at 4-year follow-up [false discovery rate (FDR) < 0.05] comparing fast-progressors with both slow-progressors and non-progressors. Among these were hsa-miR-320b, hsa-miR-320c, hsa-miR-320d, and hsa-miR-320e, which were predicted to target gene families, including members of the 14-3-3 gene family, involved in signal transduction. PLR models included miR-320 members as top predictors of fast-progressors and yielded AUC ranging from 82.6 to 91.9, representing good accuracy.

Conclusion:

The miR-320 family is associated with fast-progressing radiographic knee OA and merits further investigation as potential biomarkers and mechanistic drivers of knee OA.

Comparative Study of Classification Algorithms for Various DNA Microarray Data

Microarrays are applications of electrical engineering and technology in biology that allow simultaneous measurement of expression of numerous genes, and they can be used to analyze specific diseases. This study undertakes classification analyses of various microarrays to compare the performances of classification algorithms over different data traits. The datasets were classified into test and control groups based on five utilized machine learning methods, including MultiLayer Perceptron (MLP), Support Vector Machine (SVM), Decision Tree (DT), Random Forest (RF), and k-Nearest Neighbors (KNN), and the resulting accuracies were compared. k-fold cross-validation was used in evaluating the performance and the result was analyzed by comparing the performances of the five machine learning methods. Through the experiments, it was observed that the two tree-based methods, DT and RF, showed similar trends in results and the remaining three methods, MLP, SVM, and DT, showed similar trends. DT and RF generally showed worse performance than other methods except for one dataset. This suggests that, for the effective classification of microarray data, selecting a classification algorithm that is suitable for data traits is crucial to ensure optimum performance.

Consensus guidelines for the validation of qRT-PCR assays in clinical research by the CardioRNA consortium

Despite promising findings, quantitative PCR (qPCR)-based tests for RNA quantification have experienced serious limitations in their clinical application. The noticeable lack of technical standardization remains a huge obstacle in the translation of qPCR-based tests. The incorporation of qPCR-based tests into the clinic will benefit from guidelines for clinical research assay validation. This will ultimately impact the clinical management of the patient, including diagnosis, prognosis, prediction, monitoring of the therapeutic response, and evaluation of toxicity. However, clear assay validation protocols for biomarker investigation in clinical trials using molecular assays are currently lacking. Here, we will focus on the necessary steps, including sample acquisition, processing and storage, RNA purification, target selection, assay design, and experimental design, that need to be taken toward the appropriate validation of qRT-PCR assays in clinical research. These recommendations can fill the gap between research use only (RUO) and in vitro diagnostics (IVD). Our contribution provides a tool for basic and clinical research for the development of validated assays in the intermediate steps of biomarker research. These guidelines are based on the current understanding and consensus within the EU-CardioRNA COST Action consortium (www.cardiorna.eu). Their applicability encompasses all clinical areas.

SARS-CoV-2 causes a significant stress response mediated by small RNAs in the blood of COVID-19 patients

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a serious impact on the world. In this study, small RNAs from the blood of COVID-19 patients with moderate or severe symptoms were extracted for high-throughput sequencing and analysis. Interestingly, the levels of a special group of tRNA-derived small RNAs (tsRNAs) were found to be dramatically upregulated after SARS-CoV-2 infection, particularly in coronavirus disease 2019 (COVID-19) patients with severe symptoms. In particular, the 3′CCA tsRNAs from tRNA-Gly were highly consistent with the inflammation indicator C-reactive protein (CRP). In addition, we found that the majority of significantly changed microRNAs (miRNAs) were associated with endoplasmic reticulum (ER)/unfolded protein response (UPR) sensors, which may lead to the induction of proinflammatory cytokine and immune responses. This study found that SARS-CoV-2 infection caused significant changes in the levels of stress-associated small RNAs in patient blood and their potential functions. Our research revealed that the cells of COVID-19 patients undergo tremendous stress and respond, which can be reflected or regulated by small non-coding RNA (sncRNAs), thus providing potential thought for therapeutic intervention in COVID-19 by modulating small RNA levels or activities.

Perspectives and challenges of epigenetic determinants of childhood obesity: A systematic review

The tremendous increase in childhood obesity prevalence over the last few decades cannot merely be explained by genetics and evolutionary changes in the genome, implying that gene-environment interactions, such as epigenetic modifications, likely play a major role. This systematic review aims to summarize the evidence of the association between epigenetics and childhood obesity. A literature search was performed via PubMed and Scopus engines using a combination of terms related to epigenetics and pediatric obesity. Articles studying the association between epigenetic mechanisms (including DNA methylation and hydroxymethylation, non-coding RNAs, and chromatin and histones modification) and obesity and/or overweight (or any related anthropometric parameters) in children (0-18 years) were included. The risk of bias was assessed with a modified Newcastle-Ottawa scale for non-randomized studies. One hundred twenty-one studies explored epigenetic changes related to childhood obesity. DNA methylation was the most widely investigated mechanism (N = 101 studies), followed by non-coding RNAs (N = 19 studies) with evidence suggestive of an association with childhood obesity for DNA methylation of specific genes and microRNAs (miRNAs). One study, focusing on histones modification, was identified. Heterogeneity of findings may have hindered more insights into the epigenetic changes related to childhood obesity. Gaps and challenges that future research should face are herein described.

Diagnostic and prognostic potential of eight whole blood microRNAs for equine sarcoid disease

MicroRNAs have been proposed as biomarkers for equine sarcoids, the most prevalent equine skin tumors globally. This study served to validate the diagnostic and prognostic potential of whole blood microRNAs identified in a previous study for long-term equine sarcoid diagnosis and outcome prediction. Based on findings of a clinical examination at the age of 3 years and a follow-up following a further 5–12 years, 32 Franches-Montagnes and 45 Swiss Warmblood horses were assigned to four groups: horses with regression (n = 19), progression (n = 9), new occurrences of sarcoid lesions (n = 19) and tumor-free control horses (n = 30). The expression levels for eight microRNAs (eca-miR-127, eca-miR-432, eca-miR-24, eca-miR-125a-5p, eca-miR-134, eca-miR-379, eca-miR-381, eca-miR-382) were analyzed through reverse transcription quantitative polymerase chain reaction in whole blood samples collected on initial examination. Associations of sex, breed, diagnosis, and prognosis with microRNA expression levels were examined using multivariable analysis of variance. Sex and breed influenced the expression level of five and two microRNAs, respectively. Eca-miR-127 allowed discrimination between sarcoid-affected and tumor-free horses. No variation in microRNA expression was found when comparing horses with sarcoid regression and progression. Expression levels of eca-miR-125a-5p and eca-miR-432 varied in male horses that developed sarcoids throughout the study period in comparison to male control horses. While none of the investigated miRNAs was validated for predicting the prognosis of sarcoid regression / progression within young horses with this condition, two miRNAs demonstrated potential to predict if young male (though not female) tumor-free horse can develop sarcoids within the following years. Sex- and breed- biased miRNAs exist within the equine species and have an impact on biomarker discovery.

hsa-MiR-19a-3p and hsa-MiR-19b-3p Are Associated with Spinal Cord Injury-Induced Neuropathic Pain: Findings from a Genome-Wide MicroRNA Expression Profiling Screen

Neuropathic pain in spinal cord injury (SCI) is associated with inflammation in both the peripheral and central nervous system (CNS), which may contribute to the initiation and maintenance of persistent pain. An understanding of factors contributing to neuroinflammation may lead to new therapeutic targets for neuropathic pain. Moreover, novel circulating biomarkers of neuropathic pain may facilitate earlier and more effective treatment. MicroRNAs (miRNAs) are short, non-coding single-stranded RNA that have emerged as important biomarkers and molecular mediators in physiological and pathological conditions. Using a genome-wide miRNA screening approach, we studied differential miRNA expression in plasma from 68 healthy, community-dwelling adults with and without SCI enrolled in ongoing clinical studies. We detected 2367 distinct miRNAs. Of these, 383 miRNAs were differentially expressed in acute SCI or chronic SCI versus no SCI and 71 were differentially expressed in chronic neuropathic pain versus no neuropathic pain. We selected homo sapiens (hsa)-miR-19a-3p and hsa-miR-19b-3p for additional analysis based on p-value, fold change, and their known role as regulators of neuropathic pain and neuroinflammation. Both hsa-miR-19a-3p and hsa-miR-19b-3p levels were significantly higher in those with chronic SCI and severe neuropathic pain versus those with chronic SCI and no neuropathic pain. In confirmatory studies, both hsa-miR-19a-3p and hsa-miR-19b-3p have moderate to strong discriminative ability to distinguish between those with and without pain. After adjusting for opioid use, hsa-miR-19b-3p levels were positively associated with pain interference with mood. Because hsa-miR-19 levels have been shown to change in response to exercise, folic acid, and resveratrol, these studies suggest that miRNAs are potential targets of therapeutic interventions.

Non-Coding RNAs in Human Breast Milk: A Systematic Review

Breast milk is the primary source of nutrition and hydration for the newborn infant but also plays an important role in the child’s first immune defense. Additionally, several breast milk factors have been implicated in immune-related health outcomes later in life, including immunoglobulins, cytokines, chemokines, growth factors and, more recently, non-coding RNA (ncRNA) species. In this systematic review, we provide a comprehensive summary of the current literature on endogenous ncRNAs found in human breast milk. Thirty (30) relevant studies were identified and, whilst the majority studies focused on microRNAs (miRNAs), there is evidence that breast milk contains high quantities of RNA which also include long-coding RNAs, circular RNAs, as well as other short RNAs and fragmented tRNA and rRNAs. Among studies investigating miRNAs, miR-148a-3p, miR-30a/d-5p, miR-22-3p, miR-146b-5p, miR-200a/c-3p, and the 5p end of the let-7 miRNAs were commonly reported among the top 10 miRNAs in the cell, lipid, and skim milk fractions of breast milk. Methodological difference and small sample sizes limit the possibility of conclusively identifying which maternal and infant characteristics affect the miRNA profile. The highly expressed miRNAs were generally reported to be similar across lactational stage, milk fraction, maternal and infant characteristics, or infant growth and health. All the same, individual studies identify potential differences in miRNA expression levels which should be confirmed by future studies. Stability, uptake, and physiological functions of miRNAs were also considered in several studies. Breast milk miRNAs are relatively resistant to a range of harsh conditions and uptake experiments suggest that extracellular vesicles containing miRNAs and circular RNAs can be taken up by intestinal epithelial cells. Although the evidence regarding the functional effect of breast milk miRNAs is limited, the predicted functions range from metabolic and biosynthetic processes to signaling pathways, cellular adhesion, communication, growth, and differentiation. Finally, this systematic review highlights some of the methodological challenges and knowledge gaps which can help direct future research in this field. In particular, it is important to further investigate the bioavailability of miRNAs in different milk fractions, and to characterize other ncRNAs which are largely unstudied.

Postpartum circulating microRNA enhances prediction of future type 2 diabetes in women with previous gestational diabetes

AIMS/HYPOTHESIS

Type 2 diabetes mellitus is a major cause of morbidity and death worldwide. Women with gestational diabetes mellitus (GDM) have greater than a sevenfold higher risk of developing type 2 diabetes in later life. Accurate methods for postpartum type 2 diabetes risk stratification are lacking. Circulating microRNAs (miRNAs) are well recognised as biomarkers/mediators of metabolic disease. We aimed to determine whether postpartum circulating miRNAs can predict the development of type 2 diabetes in women with previous GDM.

METHODS

In an observational study, plasma samples were collected at 12 weeks postpartum from 103 women following GDM pregnancy. Utilising a discovery approach, we measured 754 miRNAs in plasma from type 2 diabetes non-progressors (n = 11) and type 2 diabetes progressors (n = 10) using TaqMan-based real-time PCR on an OpenArray platform. Machine learning algorithms involving penalised logistic regression followed by bootstrapping were implemented.

RESULTS

Fifteen miRNAs were selected based on their importance in discriminating type 2 diabetes progressors from non-progressors in our discovery cohort. The levels of miRNA miR-369-3p remained significantly different (p < 0.05) between progressors and non-progressors in the validation sample set (n = 82; 71 non-progressors, 11 progressors) after adjusting for age and correcting for multiple comparisons. In a clinical model of prediction of type 2 diabetes that included six traditional risk factors (age, BMI, pregnancy fasting glucose, postpartum fasting glucose, cholesterol and triacylglycerols), the addition of the circulating miR-369-3p measured at 12 weeks postpartum improved the prediction of future type 2 diabetes from traditional AUC 0.83 (95% CI 0.68, 0.97) to an AUC 0.92 (95% CI 0.84, 1.00).

CONCLUSIONS

This is the first demonstration of miRNA-based type 2 diabetes prediction in women with previous GDM. Improved prediction will facilitate early lifestyle/drug intervention for type 2 diabetes prevention.

Application of miRNA-seq in neuropsychiatry: A methodological perspective

MiRNAs are emerging as key molecules to study neuropsychiatric diseases. However, despite the large number of methodologies and software for miRNA-seq analyses, there is little supporting literature for researchers in this area. This review focuses on evaluating how miRNA-seq has been used to study neuropsychiatric diseases to date, analyzing both the main findings discovered and the bioinformatics workflows and tools used from a methodological perspective. The objective of this review is two-fold: first, to evaluate current miRNA-seq procedures used in neuropsychiatry; and second, to offer comprehensive information that can serve as a guide to new researchers in bioinformatics. After conducting a systematic search (from 2016 to June 30, 2020) of articles using miRNA-seq in neuropsychiatry, we have seen that it has already been used for different types of studies in three main categories: diagnosis, prognosis, and mechanism. We carefully analyzed the bioinformatics workflows of each study, observing a high degree of variability with respect to the tools and methods used and several methodological complexities that are identified and discussed in this review.

Peripheral blood microRNAs and the COVID-19 patient: methodological considerations, technical challenges and practice points

The COVID-19 emergency pandemic resulting from infection with SARS-CoV-2 represents a major threat to public health worldwide. There is an urgent clinical demand for easily accessible tools to address weaknesses and gaps in the management of COVID-19 patients. In this context, transcriptomic profiling of liquid biopsies, especially microRNAs (miRNAs), has recently emerged as a robust source of potential clinical indicators for medical decision-making. Nevertheless, the analysis of the circulating miRNA signature and its translation to clinical practice requires strict control of a wide array of methodological details. In this review, we indicate the main methodological aspects that should be addressed when evaluating the circulating miRNA profiles in COVID-19 patients, from preanalytical and analytical variables to the experimental design, impact of confounding, analysis of the data and interpretation of the findings, among others. Additionally, we provide practice points to ensure the rigour and reproducibility of miRNA-based biomarker investigations of this condition.Abbreviations: ACE: angiotensin-converting enzyme; ARDS: acute respiratory distress syndrome; COVID-19: coronavirus disease 2019; ERDN: early Detection Research Network; LMWH: low molecular weight heparin; miRNA: microRNA; ncRNA: noncoding RNA; SARS-CoV-2: severe acute respiratory syndrome coronavirus-2; SOP: standard operating procedure.

Plasma microRNA expression levels in HIV-1-positive patients receiving antiretroviral therapy

MicroRNAs (miRNAs/miRs) may serve as therapeutic agents or targets in diseases in which the expression of proteins plays an important role. The aim of the present study was to compare the expression levels of specific miRNAs, as well as their correlation with markers of response to antiretroviral (ARV) therapy, in patients with human immunodeficiency virus type 1 (HIV-1) infection with and without resistance to highly active antiretroviral therapy (HAART). Methods: miRNA assays were performed on plasma samples obtained from 20 HIV-1-positive patients. A total of ten patients were divided into two groups: HAART-responsive and HAART-resistant (n=5 per group). Commercial arrays were subsequently used to identify 84 miRNAs. A total of three differentially expressed miRNAs were selected and analyzed by quantitative PCR (qPCR). Five other patients were subsequently added to each group for a new relative expression analysis. The absolute expression level of the two miRNAs was obtained and compared using the Student’s t test. Receiver operating characteristic (ROC) curves were used to identify patients with antiretroviral therapy (ART) resistance.

Results: The array analysis revealed that miR-15b-5p, miR-16-5p, miR-20a-5p, miR-26a-5p, miR-126-3p and miR-150-5p were down-regulated in patients with HAART-resistance comparing with HAART-responsive. The expression levels of miR-16-5p, miR-26a-5p and miR-150-5p were confirmed using qPCR. The area under the ROC curve was 1.0 for the three miRNAs.

Conclusions: The lower expression levels of miR-16-5p and miR-26a-5p in patients with HAART-resistance suggested that these may serve as potential biomarkers for the identification of HAART-responsive patients.

Improvement, identification, and target prediction for miRNAs in the porcine genome by using massive, public high-throughput sequencing data

Despite the broad variety of available microRNA (miRNA) research tools and methods, their application to the identification, annotation, and target prediction of miRNAs in nonmodel organisms is still limited. In this study, we collected nearly all public sRNA-seq data to improve the annotation for known miRNAs and identify novel miRNAs that have not been annotated in pigs (Sus scrofa). We newly annotated 210 mature sequences in known miRNAs and found that 43 of the known miRNA precursors were problematic due to redundant/missing annotations or incorrect sequences. We also predicted 811 novel miRNAs with high confidence, which was twice the current number of known miRNAs for pigs in miRBase. In addition, we proposed a correlation-based strategy to predict target genes for miRNAs by using a large amount of sRNA-seq and RNA-seq data. We found that the correlation-based strategy provided additional evidence of expression compared with traditional target prediction methods. The correlation-based strategy also identified the regulatory pairs that were controlled by nonbinding sites with a particular pattern, which provided abundant complementarity for studying the mechanism of miRNAs that regulate gene expression. In summary, our study improved the annotation of known miRNAs, identified a large number of novel miRNAs, and predicted target genes for all pig miRNAs by using massive public data. This large data-based strategy is also applicable for other nonmodel organisms with incomplete annotation information.

MiR-223-3p and miR-122-5p as circulating biomarkers for plaque instability

BACKGROUND

In this study, we discovered and validated candidate microRNA (miRNA) biomarkers for coronary artery disease (CAD).

METHOD

Candidate tissue-derived miRNAs from atherosclerotic plaque material in patients with stable coronary artery disease (SCAD) (n=14) and unstable coronary artery disease (UCAD) (n=25) were discovered by qPCR-based arrays. We validated differentially expressed miRNAs, along with seven promising CAD-associated miRNAs from the literature, in the serum of two large cohorts (n=395 and n=1000) of patients with SCAD and UCAD and subclinical atherosclerosis (SubA) and controls, respectively.

RESULT

From plaque materials (discovery phase), miR-125b-5p and miR-193b-3p were most upregulated in SCAD, whereas miR-223-3p and miR-142-3p were most upregulated in patients with UCAD. Subsequent validation in serum from patients with UCAD, SCAD, SubA and controls demonstrated significant upregulation of miR-223-3p, miR-133a-3p, miR-146-3p and miR-155-5p. The ischaemia-related miR-499-5p was also highly upregulated in patients with UCAD compared with the other groups (SCAD OR 20.63 (95% CI 11.16 to 38.15), SubA OR 96.10 (95% CI 40.13 to 230.14) and controls OR 15.73 (95% CI 7.80 to 31.72)). However, no significant difference in miR-499-5p expression was observed across SCAD, SubA and controls. MiR-122-5p was the only miRNA to be significantly upregulated in the serum of both patients with UCAD and SCAD.

CONCLUSION

In conclusion, miR-122-5p and miR-223-3p might be markers of plaque instability.

State of the Field in Multi-Omics Research: From Computational Needs to Data Mining and Sharing

Multi-omics, variously called integrated omics, pan-omics, and trans-omics, aims to combine two or more omics data sets to aid in data analysis, visualization and interpretation to determine the mechanism of a biological process. Multi-omics efforts have taken center stage in biomedical research leading to the development of new insights into biological events and processes. However, the mushrooming of a myriad of tools, datasets, and approaches tends to inundate the literature and overwhelm researchers new to the field. The aims of this review are to provide an overview of the current state of the field, inform on available reliable resources, discuss the application of statistics and machine/deep learning in multi-omics analyses, discuss findable, accessible, interoperable, reusable (FAIR) research, and point to best practices in benchmarking. Thus, we provide guidance to interested users of the domain by addressing challenges of the underlying biology, giving an overview of the available toolset, addressing common pitfalls, and acknowledging current methods' limitations. We conclude with practical advice and recommendations on software engineering and reproducibility practices to share a comprehensive awareness with new researchers in multi-omics for end-to-end workflow.

Sequencing identifies a distinct signature of circulating microRNAs in early radiographic knee osteoarthritis

OBJECTIVE

MicroRNAs act locally and systemically to impact osteoarthritis (OA) pathophysiology, but comprehensive profiling of the circulating miRNome in early vs late stages of OA has yet to be conducted. Sequencing has emerged as the preferred method for microRNA profiling since it offers high sensitivity and specificity. Our objective was to sequence the miRNome in plasma from 91 patients with early [Kellgren-Lawrence (KL) grade 0 or 1 (n = 41)] or late [KL grade 3 or 4 (n = 50)] symptomatic radiographic knee OA to identify unique microRNA signatures in each disease state.

DESIGN

MicroRNA libraries were prepared using the QIAseq miRNA Library Kit and sequenced on the Illumina NextSeq 550. Counts were produced for microRNAs captured in miRBase and for novel microRNAs. Statistical, bioinformatics, and computational biology approaches were used to refine and interpret the final list of microRNAs.

RESULTS

From 215 differentially expressed microRNAs (FDR < 0.01), 97 microRNAs showed an increase or decrease in expression in ≥85% of samples in the early OA group as compared to the median expression in the late OA group. Increasing this threshold to ≥95%, seven microRNAs were identified: hsa-miR-335-3p, hsa-miR-199a-5p, hsa-miR-671-3p, hsa-miR-1260b, hsa-miR-191-3p, hsa-miR-335-5p, and hsa-miR-543. Four novel microRNAs were present in ≥50% of early OA samples and had 27 predicted gene targets in common with the prioritized set of predicted gene targets from the 97 microRNAs, suggesting common underlying mechanisms.

CONCLUSION

Sequencing of well-characterized patient cohorts produced unbiased profiling of the circulating miRNome and identified a unique panel of 11 microRNAs in early radiographic knee OA.

Expression differences of miR-142-5p between treatment-naïve chronic myeloid leukemia patients responding and non-responding to imatinib therapy suggest a link to oncogenic ABL2, SRI, cKIT and MCL1 signaling pathways critical for development of therapy resistance

Background

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm characterized by constitutive activity of the tyrosine kinase BCR-ABL1. Although the introduction of tyrosine kinase inhibitors (TKIs) has substantially improved patients’ prognosis, drug resistance remains one of the major challenges in CML therapy. MicroRNAs (miRNAs), a class of short non-coding RNAs acting as post-transcriptional regulators, are implicated in CML progression and drug resistance. The aim of the present study was to analyze the miRNA expression profiles of 45 treatment-naïve CML patients in chronic phase (28 peripheral blood and 17 bone marrow samples) with respect to future response to imatinib therapy.

Methods

TaqMan low density arrays were used to analyze the miRNA expression pattern of the patient samples. For selected microRNAs, reporter gene assays were performed to study their ability to regulate CML associated target genes.

Results

Significant lower expression levels of miR-142-5p were identified in both, peripheral blood and bone marrow samples of future non-responders suggesting a potential tumor suppressor role of this miRNA. This was supported by reporter gene assays that identified the survival, proliferation and invasion promoting CML related genes ABL2, cKIT, MCL1 and SRI as targets of miR-142-5p and miR-365a-3p, the latter identified as potential biomarker in peripheral blood samples.

Conclusion

MiR-142-5p and to a certain extend also miR-365a-3p were able to discriminate treatment-naïve CML patients not responding to imatinib in the course of their treatment from patients, who responded to therapy. However, further large-scale studies should clarify if the identified miRNAs have the potential as predictive biomarkers for TKI resistance.

Circulating microRNA in Heart Failure - Practical Guidebook to Clinical Application

Heart failure (HF) remains a major cause of death and disability worldwide. Currently, B-type natriuretic peptide and N-terminal pro-brain natriuretic peptide are diagnostic biomarkers used in HF. Although very sensitive, they are not specific enough and do not allow the prediction or early diagnosis of HF. Many ongoing studies focus on determining the underlying cause and understanding the mechanisms of HF on the cellular level. MicroRNAs (miRNAs) are non-coding RNAs which control the majority of cellular processes and therefore are considered to have a potential clinical application in HF. In this review, we aim to provide synthesized information about miRNAs associated with ejection fraction, HF etiology, diagnosis, and prognosis, as well as outline therapeutic application of miRNAs in HF. Further, we discuss methodological challenges associated with the analysis of miRNAs and provide recommendations for defining a study population, collecting blood samples, and selecting detection methods to study miRNAs in a reliable and reproducible way. This review is intended to be an accessible tool for clinicians interested in the field of miRNAs and HF.

Placental microRNAs: Responders to environmental chemicals and mediators of pathophysiology of the human placenta

MicroRNAs (miRNAs) are epigenetic modifiers that play an important role in the regulation of the expression of genes across the genome. miRNAs are expressed in the placenta as well as other organs, and are involved in several biological processes including the regulation of trophoblast differentiation, migration, invasion, proliferation, apoptosis, angiogenesis and cellular metabolism. Related to their role in disease process, miRNAs have been shown to be differentially expressed between normal placentas and placentas obtained from women with pregnancy/health complications such as preeclampsia, gestational diabetes mellitus, and obesity. This dysregulation indicates that miRNAs in the placenta likely play important roles in the pathogenesis of diseases during pregnancy. Furthermore, miRNAs in the placenta are susceptible to altered expression in relation to exposure to environmental toxicants. With relevance to the placenta, the dysregulation of miRNAs in both placenta and blood has been associated with maternal exposures to several toxicants. In this review, we provide a summary of miRNAs that have been assessed in the context of human pregnancy-related diseases and in relation to exposure to environmental toxicants in the placenta. Where data are available, miRNAs are discussed in their context as biomarkers of exposure and/or disease, with comparisons made across-tissue types, and conservation across studies detailed.

MiRNA-200c, MiRNA-139 and ln RNA H19; new predictors of treatment response in H-pylori- induced gastric ulcer or progression to gastric cancer

Recent evidence indicates that the pathogenesis of gastric ulcer and progression to gastric cancer could be attributed to altered inflammatory/immunological response and associated differential non-coding RNAs expression signatures. However, co-expression profiling of lncRNA-miRNAs in GU/GC patients are scarcely focused on. Therefore, in the present study the expression of H19 and related miRNAs including miR-139, and miR-200 were assayed in the plasma samples of treatment responsive GU vs nonresponsive GC patients. This study is a case-control study carried out on 130 subjects recruited from the Gastrointestinal Endoscopy Unit in Al-Kasr Al-Aini Hospital, in Egypt. All recruited patients were diagnosed with H-pylori infection, 50 of them were gastric cancer patients (GC), with previous H-pylori induced gastric ulcer but were treatment non-respondent. Real-time PCR was performed to evaluate the expression level of serum non-coding RNA; miRNA-200c, miR-139, Ln RNA H19 in patients with peptic ulcer treatment non-respondent, who progressed to GC vs non-progressed gastric ulcer patients (GU) (n = 50), and compared to early diagnosed H-pylori-gastric ulcer patients (n = 30). The association between these miRNAs and the FGF-18/FGF-R signaling indicators of H-pylori-GC pathogenesis were then investigated. RESULTS: showed that the H19 level was significantly elevated while miR-139 and miR-200c expression were significantly down-regulated in GC patients, compared to GU participants (P < 0.01). The herein investigated ncRNAs are correlated to the disease duration with Ln H19 being significantly correlated with all inflammatory markers; TNF-α, INF-γ, TAC, MMP-9, and FGF18/FGFR2. A significant correlation was also observed between miRNA 200c and each of miRNA 139 and FGFR2. Moreover, ROC analysis revealed that miRNA 200c showed the highest AUC (0.906) and 81.2% sensitivity and 100% specificity. Moreover, the combined analysis of miRNA 200c/miRNA 139 revealed superior AUC (0.96) and 93% sensitivity and 100% specificity, than each separately. As for discriminative accuracy between stages III to IV of gastric cancer, LncRNA H19 showed the highest diagnostic accuracy (95.5%), specificity (100%), and sensitivity (90.9%). The current study demonstrated that the combination of serum miRNA 200c/miRNA 139 expression levels (down-regulation) could provide a new potential prognostic panel for GU predictive response and potential sequelae. In conclusion, LncRNA H19 and related miRNAs, miRNA 200c/miRNA 139, could serve as a potential diagnostic biomarker for early gastric cancer diagnosis.

Exploiting Circulating MicroRNAs as Biomarkers in Psychiatric Disorders

Non-invasive peripheral biomarkers play a significant role in both disease diagnosis and progression. In the past few years, microRNA (miRNA) expression changes in circulating peripheral tissues have been found to be correlative with changes in neuronal tissues from patients with neuropsychiatric disorders. This is a notable quality of a biomolecule to be considered as a biomarker for both prognosis and diagnosis of disease. miRNAs, members of the small non-coding RNA family, have recently gained significant attention due to their ability to epigenetically influence almost every aspect of brain functioning. Empirical evidence suggests that miRNA-associated changes in the brain are often translated into behavioral changes. Current clinical understanding further implicates their role in the management of major psychiatric conditions, including major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SZ). This review aims to critically evaluate the potential advantages and disadvantages of miRNAs as diagnostic/prognostic biomarkers in psychiatric disorders as well as in treatment response.