Optimization of miRNA-seq data preprocessing

Identifying miRNA as biomarker for breast cancer subtyping using association rule

- This paper presents a comprehensive study focused on breast cancer subtyping, utilizing a multifaceted approach that integrates feature selection, machine learning classifiers, and miRNA regulatory networks. The feature selection process begins with the CFS algorithm, followed by the Apriori algorithm for association rule generation, resulting in the identification of significant features tailored to Luminal A, Luminal B, HER-2 enriched, and Basal-like subtypes. The subsequent application of Random Forest (RF) and Support Vector Machine (SVM) classifiers yielded promising results, with the SVM model achieving an overall accuracy of 76.60 % and the RF model demonstrating robust performance at 80.85 %. Detailed accuracy metrics revealed strengths and areas for refinement, emphasizing the potential for optimizing subtype-specific recall. To explore the regulatory landscape in depth, an analysis of selected miRNAs was conducted using MIENTURNET, a tool for visualizing miRNA-target interactions. While FDR analysis raised concerns for HER-2 and Basal-like subtypes, Luminal A and Luminal B subtypes showcased significant miRNA-gene interactions. Functional enrichment analysis for Luminal A highlighted the role of Ovarian steroidogenesis, implicating specific miRNAs such as hsa-let-7c-5p and hsa-miR-125b-5p as potential diagnostic biomarkers and regulators of Luminal A breast cancer. Luminal B analysis uncovered associations with the MAPK signaling pathway, with miRNAs like hsa-miR-203a-3p and hsa-miR-19a-3p exhibiting potential diagnostic and therapeutic significance. In conclusion, this integrative approach combines machine learning techniques with miRNA analysis to provide a holistic understanding of breast cancer subtypes. The identified miRNAs and associated pathways offer insights into potential diagnostic biomarkers and therapeutic targets, contributing to the ongoing efforts to improve breast cancer diagnostics and personalized treatment strategies.

Dynamic changes in extracellular vesicle-associated miRNAs elicited by ultrasound in inflammatory bowel disease patients

Blood-based biomarkers that reliably indicate disease activity in the intestinal tract are an important unmet need in the management of patients with IBD. Extracellular vesicles (EVs) are cell-derived membranous microparticles, which reflect the cellular and functional state of their site of site of origin. As ultrasound waves may lead to molecular shifts of EV contents, we hypothesized that application of ultrasound waves on inflamed intestinal tissue in IBD may amplify the inflammation-specific molecular shifts in EVs like altered EV-miRNA expression, which in turn can be detected in the peripheral blood. 26 patients with IBD were included in the prospective clinical study. Serum samples were collected before and 30 min after diagnostic transabdominal ultrasound. Differential miRNA expression was analyzed by sequencing. Candidate inducible EV-miRNAs were functionally assessed in vitro by transfection of miRNA mimics and qPCR of predicted target genes. Serum EV-miRNA concentration at baseline correlated with disease severity, as determined by clinical activity scores and sonographic findings. Three miRNAs (miR-942-5p, mir-5588, mir-3195) were significantly induced by sonography. Among the significantly regulated EV-miRNAs, miR-942-5p was strongly induced in higher grade intestinal inflammation and correlated with clinical activity in Crohn's disease. Prediction of target regulation and transfection of miRNA mimics inferred a role of this EV-miRNA in regulating barrier function in inflammation. Induction of mir-5588 and mir-3195 did not correlate with inflammation grade. This proof-of-concept trial highlights the principle of induced molecular shifts in EVs from inflamed tissue through transabdominal ultrasound. These inducible EVs and their molecular cargo like miRNA could become novel biomarkers for intestinal inflammation in IBD.

Nomogram for soiling prediction in postsurgery hirschsprung children: a retrospective study

PURPOSE

The aim of this study was to develop a nomogram for predicting the probability of postoperative soiling in patients aged greater than 1 year operated for Hirschsprung disease (HSCR).

MATERIALS AND METHODS

The authors retrospectively analyzed HSCR patients with surgical therapy over 1 year of age from January 2000 and December 2019 at our department. Eligible patients were randomly categorized into the training and validation set at a ratio of 7:3. By integrating the least absolute shrinkage and selection operator [LASSO] and multivariable logistic regression analysis, crucial variables were determined for establishment of the nomogram. And, the performance of nomogram was evaluated by C-index, area under the receiver operating characteristic curve, calibration curves, and decision curve analysis. Meanwhile, a validation set was used to further assess the model.

RESULTS

This study enrolled 601 cases, and 97 patients suffered from soiling. Three risk factors, including surgical history, length of removed bowel, and surgical procedures were identified as predictive factors for soiling occurrence. The C-index was 0.871 (95% CI: 0.821-0.921) in the training set and 0.878 (95% CI: 0.811-0.945) in the validation set, respectively. And, the AUC was found to be 0.896 (95% CI: 0.855-0.929) in the training set and 0.866 (95% CI: 0.767-0.920) in the validation set. Additionally, the calibration curves displayed a favorable agreement between the nomogram model and actual observations. The decision curve analysis revealed that employing the nomogram to predict the risk of soiling occurrence would be advantageous if the threshold was between 1 and 73% in the training set and 3-69% in the validation set.

CONCLUSION

This study represents the first efforts to develop and validate a model capable of predicting the postoperative risk of soiling in patients aged greater than 1 year operated for HSCR. This model may assist clinicians in determining the individual risk of soiling subsequent to HSCR surgery, aiding in personalized patient care and management.

Nuclear localization of Argonaute 2 is affected by cell density and may relieve repression by microRNAs

Argonaute protein is associated with post-transcriptional control of cytoplasmic gene expression through miRNA-induced silencing complexes (miRISC). Specific cellular and environmental conditions can trigger AGO protein to accumulate in the nucleus. Localization of AGO is central to understanding miRNA action, yet the consequences of AGO being in the nucleus are undefined. We show nuclear enrichment of AGO2 in HCT116 cells grown in two-dimensional culture to high density, HCT116 cells grown in three-dimensional tumor spheroid culture, and human colon tumors. The shift in localization of AGO2 from cytoplasm to nucleus de-represses cytoplasmic AGO2-eCLIP targets that were candidates for canonical regulation by miRISC. Constitutive nuclear localization of AGO2 using an engineered nuclear localization signal increases cell migration. Critical RNAi factors also affect the localization of AGO2. Knocking out an enzyme essential for miRNA biogenesis, DROSHA, depletes mature miRNAs and restricts AGO2 localization to the cytoplasm, while knocking out the miRISC scaffolding protein, TNRC6, results in nuclear localization of AGO2. These data suggest that AGO2 localization and miRNA activity can be regulated depending on environmental conditions, expression of mature miRNAs, and expression of miRISC cofactors. Localization and expression of core miRISC protein machinery should be considered when investigating the roles of miRNAs.

Functional impact of multi-omic interactions in lung cancer

Lung tumors are a leading cause of cancer-related death worldwide. Lung cancers are highly heterogeneous on their phenotypes, both at the cellular and molecular levels. Efforts to better understand the biological origins and outcomes of lung cancer in terms of this enormous variability often require of high-throughput experimental techniques paired with advanced data analytics. Anticipated advancements in multi-omic methodologies hold potential to reveal a broader molecular perspective of these tumors. This study introduces a theoretical and computational framework for generating network models depicting regulatory constraints on biological functions in a semi-automated way. The approach successfully identifies enriched functions in analyzed omics data, focusing on Adenocarcinoma (LUAD) and Squamous cell carcinoma (LUSC, a type of NSCLC) in the lung. Valuable information about novel regulatory characteristics, supported by robust biological reasoning, is illustrated, for instance by considering the role of genes, miRNAs and CpG sites associated with NSCLC, both novel and previously reported. Utilizing multi-omic regulatory networks, we constructed robust models elucidating omics data interconnectedness, enabling systematic generation of mechanistic hypotheses. These findings offer insights into complex regulatory mechanisms underlying these cancer types, paving the way for further exploring their molecular complexity.

Systematic computational hunting for small RNAs derived from ncRNAs during dengue virus infection in endothelial HMEC-1 cells

In recent years, a population of small RNA fragments derived from non-coding RNAs (sfd-RNAs) has gained significant interest due to its functional and structural resemblance to miRNAs, adding another level of complexity to our comprehension of small-RNA-mediated gene regulation. Despite this, scientists need more tools to test the differential expression of sfd-RNAs since the current methods to detect miRNAs may not be directly applied to them. The primary reasons are the lack of accurate small RNA and ncRNA annotation, the multi-mapping read (MMR) placement, and the multicopy nature of ncRNAs in the human genome. To solve these issues, a methodology that allows the detection of differentially expressed sfd-RNAs, including canonical miRNAs, by using an integrated copy-number-corrected ncRNA annotation was implemented. This approach was coupled with sixteen different computational strategies composed of combinations of four aligners and four normalization methods to provide a rank-order of prediction for each differentially expressed sfd-RNA. By systematically addressing the three main problems, we could detect differentially expressed miRNAs and sfd-RNAs in dengue virus-infected human dermal microvascular endothelial cells. Although more biological evaluations are required, two molecular targets of the hsa-mir-103a and hsa-mir-494 (CDK5 and PI3/AKT) appear relevant for dengue virus (DENV) infections. Here, we performed a comprehensive annotation and differential expression analysis, which can be applied in other studies addressing the role of small fragment RNA populations derived from ncRNAs in virus infection.

Imprinted Dlk1-Gtl2 cluster miRNAs are potential epigenetic regulators of lamb fur quality

BACKGROUND

Tan and Hu sheep are well-known local breeds in China, producing lamb fur with unique ornamental and practical values highly appreciated by consumers worldwide. Fur quality is optimal at one month of age and gradually declines with time. Despite active research on its genetic mechanism using transcriptomic and whole genome bisulfite sequencing analysis, the main effective gene locus has not been found, and its regulatory mechanism is still unclear, which limits the breeding and improvement of fur traits.

RESULTS

Scapular skin samples from newborn (1-month old) and adult (24-month old) Tan sheep were utilized for small ribonucleic acid (RNA) sequencing Principal Component Analysis (PCA) showed that the newborn and adult groups were completely separated. Differential expression analysis of micro-RNAs (miRNAs) identified 32 up-regulated miRNAs and 48 down-regulated miRNAs in the newborn groups. All up-regulated miRNAs were located in the imprinted. Dlk1-Gtl2 locus on chromosome 18, whereas all down-regulated miRNAs were distributed across the sheep chromosomes, without a clear pattern of positional consistency. Further, by systematically analyzing the target genes and signaling pathways of all 32 up-regulated miRNAs, we found that the PI3K-AKT signaling pathway has the potential to be targeted and regulated by most of the miRNAs in the Dlk1-Gtl2 region. In addition, we also re-analyzed miRNA sequencing data from public databases on Hu lambs (full sibling Hu lambs with high- and low-quality fur characteristics). Again, it was found that most of the up-regulated miRNAs in lambs with high-quality fur were also located in the Dlk1-Gtl2 region, whereas this patter was not present for down-regulated miRNAs.

CONCLUSION

Sequencing of miRNAs in conjunction with public databases was employed to identify miRNAs within the imprinted Dlk1-Gtl2 region on chromosome 18, suggesting their potential roles as epigenetic regulators of fur traits. Small RNAs located at the Dlk1-Gtl2 locus were identified as having the potential to systematically regulate the PI3K-AKT signaling pathway, thereby indicating the relevance of the Dlk1-Gtl2/PI3K-AKT axis in the context of fur traits. Selection of parental specific expressed imprinted genes in the process of conserving and exploiting lamb fur traits should be emphasized.

Consequences of depleting TNRC6, AGO, and DROSHA proteins on expression of microRNAs

The potential for microRNAs (miRNAs) to regulate gene expression remains incompletely understood. DROSHA initiates the biogenesis of miRNAs while variants of Argonaute (AGO) and trinucleotide repeat containing six (TNRC6) family proteins form complexes with miRNAs to facilitate RNA recognition and gene regulation. Here we investigate the fate of miRNAs in the absence of these critical RNAi protein factors. Knockout of DROSHA expression reduces levels of some miRNAs annotated in miRBase but not others. The identity of miRNAs with reduced expression matches the identity of miRNAs previously identified by experimental approaches. The MirGeneDB resource offers the closest alignment with experimental results. In contrast, the loss of TNRC6 proteins had much smaller effects on miRNA levels. Knocking out AGO proteins, which directly contact the mature miRNA, decreased expression of the miRNAs most strongly associated with AGO2 as determined from enhanced crosslinking immunoprecipitation (AGO2-eCLIP). Evaluation of miRNA binding to endogenously expressed AGO proteins revealed that miRNA:AGO association was similar for AGO1, AGO2, AGO3, and AGO4. Our data emphasize the need to evaluate annotated miRNAs based on approximate cellular abundance, DROSHA-dependence, and physical association with AGO when forming hypotheses related to their function.

Nuclear Localization of Argonaute is affected by Cell Density and May Relieve Repression by microRNAs

Argonaute protein is associated with post-transcriptional control of cytoplasmic gene expression through miRNA-induced silencing complexes (miRISC). Specific cellular and environmental conditions can trigger AGO protein to accumulate in the nucleus. Localization of AGO is central to understanding miRNA action, yet the consequences of AGO being in the nucleus are undefined. We show nuclear enrichment of AGO2 in HCT116 cells grown in two-dimensional culture to high density, HCT116 cells grown in three-dimensional tumor spheroid culture, and human colon tumors. The shift in localization of AGO2 from cytoplasm to nucleus de-represses cytoplasmic AGO2-eCLIP targets that were candidates for canonical regulation by miRISC. Constitutive nuclear localization of AGO2 using an engineered nuclear localization signal increases cell migration. Critical RNAi factors also affect the localization of AGO2. Knocking out an enzyme essential for miRNA biogenesis, DROSHA, depletes mature miRNAs and restricts AGO2 localization to the cytoplasm, while knocking out the miRISC scaffolding protein, TNRC6, results in nuclear localization of AGO2. These data suggest that AGO2 localization and miRNA activity can be regulated depending on environmental conditions, expression of mature miRNAs, and expression of miRISC cofactors. Localization and expression of core miRISC protein machinery should be considered when investigating the roles of miRNAs.

Drug repurposing for Basal breast cancer subpopulations using modular network signatures

Breast cancer is characterized as being a heterogeneous pathology with a broad phenotype variability. Breast cancer subtypes have been developed in order to capture some of this heterogeneity. Each of these breast cancer subtypes, in turns retains varied characteristic features impacting diagnostic, prognostic and therapeutics. Basal breast tumors, in particular have been challenging in these regards. Basal breast cancer is often more aggressive, of rapid evolution and no tailor-made targeted therapies are available yet to treat it. Arguably, epigenetic variability is behind some of these intricacies. It is possible to further classify basal breast tumor in groups based on their non-coding transcriptome and methylome profiles. It is expected that these groups will have differences in survival as well as in sensitivity to certain classes of drugs. With this in mind, we implemented a computational learning approach to infer different subpopulations of basal breast cancer (from TCGA multi-omic data) based on their epigenetic signatures. Such epigenomic signatures were associated with different survival profiles; we then identified their associated gene co-expression network structure, extracted a signature based on modules within these networks, and use these signatures to find and prioritize drugs (in the LINCS dataset) that may be used to target these types of cancer. In this way we are introducing the analytical workflow for an epigenomic signature-based drug repurposing structure.

An exploratory approach to identify microRNAs as circulatory biomarker candidates for epilepsy-associated psychiatric comorbidities in an electrical post-status epilepticus model

Patients with epilepsy have a high risk of developing psychiatric comorbidities, and there is a particular need for early detection of these comorbidities. Here, in an exploratory, hypothesis-generating approach, we aimed to identify microRNAs as potential circulatory biomarkers for epilepsy-associated psychiatric comorbidities across different rat models of epilepsy. The identification of distress-associated biomarkers can also contribute to animal welfare assessment. MicroRNA expression profiles were analyzed in blood samples from the electrical post-status epilepticus (SE) model. Preselected microRNAs were correlated with behavioral and biochemical parameters in the electrical post-SE model, followed by quantitative real-time PCR validation in three additional well-described rat models of epilepsy. Six microRNAs (miR-376a, miR-429, miR-494, miR-697, miR-763, miR-1903) were identified showing a positive correlation with weight gain in the early post-insult phase as well as a negative correlation with social interaction, saccharin preference, and plasma BDNF. Real-time PCR validation confirmed miR-203, miR-429, and miR-712 as differentially expressed with miR-429 being upregulated across epilepsy models. While readouts from the electrical post-SE model suggest different microRNA candidates for psychiatric comorbidities, cross-model analysis argues against generalizability across models. Thus, further research is necessary to compare the predictive validity of rodent epilepsy models for detection and management of psychiatric comorbidities.

Rex: R-linked EXcel add-in for statistical analysis of medical and bioinformatics data

BACKGROUND

Microsoft Excel has substantial functionalities for data management and analyses, and has been the most popular software in this field. However, in spite of Excel's user-friendly interface and functionality for data management, it provides very few functions for in-depth statistical analyses, which has limited its wider application for this purpose.

OBJECTIVE

Here, we introduce Rex, an Excel add-in software implementing the powerful analytical and graphical functions of R within Excel.

METHODS

Rex was implemented using three types of programming software: R, JavaScript, and Microsoft VB.Net.

RESULTS

Rex provides a graphical user interface (GUI) through Excel, and statistical analysis can be conducted by pointing and clicking the menu without programming R. Rex covers a wide range of analyses from basic statistics to advanced analysis, including structural equation modeling, complex sampling design, and machine learning models, making it possible for researchers not skilled in using a command-line interface to conduct in-depth statistical analyses. Most Rex modules are available in a free version for non-commercial use, and it can be used for educational and public purposes.

CONCLUSION

In this article, we introduce the framework and features of Rex with illustrative examples of its implementation.

Distinct Cellular Origins and Differentiation Process Account for Distinct Oncogenic and Clinical Behaviors of Leiomyosarcomas

In leiomyosarcoma (LMS), a very aggressive disease, a relatively transcriptionally uniform subgroup of well-differentiated tumors has been described and is associated with poor survival. The question raised how differentiation and tumor progression, two apparently antagonist processes, coexist and allow tumor malignancy. We first identified the most transcriptionally homogeneous LMS subgroup in three independent cohorts, which we named 'hLMS'. The integration of multi-omics data and functional analysis suggests that hLMS originate from vascular smooth muscle cells and show that hLMS transcriptional program reflects both modulations of smooth muscle contraction activity controlled by MYOCD/SRF regulatory network and activation of the cell cycle activity controlled by E2F/RB1 pathway. We propose that the phenotypic plasticity of vascular smooth muscle cells coupled with MYOCD/SRF pathway amplification, essential for hLMS survival, concomitant with PTEN absence and RB1 alteration, could explain how hLMS balance this uncommon interplay between differentiation and aggressiveness.

Functional impact of multi-omic interactions in breast cancer subtypes

Multi-omic approaches are expected to deliver a broader molecular view of cancer. However, the promised mechanistic explanations have not quite settled yet. Here, we propose a theoretical and computational analysis framework to semi-automatically produce network models of the regulatory constraints influencing a biological function. This way, we identified functions significantly enriched on the analyzed omics and described associated features, for each of the four breast cancer molecular subtypes. For instance, we identified functions sustaining over-representation of invasion-related processes in the basal subtype and DNA modification processes in the normal tissue. We found limited overlap on the omics-associated functions between subtypes; however, a startling feature intersection within subtype functions also emerged. The examples presented highlight new, potentially regulatory features, with sound biological reasons to expect a connection with the functions. Multi-omic regulatory networks thus constitute reliable models of the way omics are connected, demonstrating a capability for systematic generation of mechanistic hypothesis.

Association between lipid peroxidation biomarkers and microRNA expression profiles

BACKGROUND

In-vitro and animal studies demonstrate that epigenetic regulation may play an important role in lipid peroxidation. No human study to date has directly evaluated microRNAs (miRNAs), as epigenetic modulators, in relation to systemic levels of lipid peroxidation.

OBJECTIVES

To evaluate associations between systemic levels of lipid peroxidation and miRNA expression profiles in women.

METHODS

Included in the analysis were 92 women aged 40-70 years, a subset of the Shanghai Women's Health Study (SWHS). Lipid peroxidation was assessed by urinary markers F₂-isoprostanes (F₂-IsoPs), the products of free radical-catalyzed peroxidation of arachidonic acid, and its major metabolite after β-oxidation, 2,3-dinor-5,6-dihydro-15-F_2t-IsoP (F₂-IsoP-M), with GC/NICI-MS assays. Expression levels of 798 miRNAs were quantified in peripheral plasma with NanoString nCounter assays. A multivariable linear regression model was used to examine the association between lipid peroxidation and miRNA expression.

RESULTS

After adjusting for potential confounders, 29 miRNAs and 213 miRNAs were associated with F₂-IsoPs and F₂-IsoP-M, respectively. When further controlling for multiple comparisons, none of these nominally significant associations with F₂-IsoPs was retained, whereas 71 of 213 miRNAs associated with F₂-IsoP-M remained. The predicted targets of the F₂-IsoP-M associated miRNAs were enriched for several lipid peroxidation-related processes such as PI3K/AKT, MAPK, FOXO and HIF-1 signaling pathways. Moreover, 10 miRNAs (miR-93-5p, miR-761, miR-301b-3p, miR-497-5p, miR-141-3p, miR-186-5p, miR-126-3p, miR-200b-3p, miR-520d-3p, and miR-363-3p) exhibited functional interactions with 50 unique mRNAs targets involved in the regulation of β-oxidation.

CONCLUSIONS

To our knowledge, this study, for the first time, provides human data suggesting that miRNA expression may be linked to lipid peroxidation products and their metabolism.

MicroRNA Content of Ewing Sarcoma Derived Extracellular Vesicles Leads to Biomarker Potential and Identification of a Previously Undocumented EWS-FLI1 Translocation

Objective: Ewing Sarcoma Family of Tumors (ESFT) are a highly aggressive pediatric bone and soft tissue malignancy with poor outcomes in the refractory and recurrent setting. Over 90% of Ewing Sarcoma (ES) tumors are driven by the pathognomonic EWS-ETS chimeric transcripts and their corresponding oncoproteins. It has been suggested that the EWS-ETS oncogenic action can mediate microRNA (miRNA) processing. Importantly, small extracellular vesicles (sEVs), including those frequently referred to as exosomes have been shown to be highly enriched with tumor-derived small RNAs such as miRNAs. We hypothesized that ESFT-specific sEVs are enriched with certain miRNAs which could be utilized toward an exo-miRNA biomarker signature specific to this disease. Methods: We performed miRNAseq to compare both the exo-derived and cell-derived miRNA content from 8 ESFT, 2 osteosarcoma, 2 non-cancerous cell lines, and pediatric plasma samples. Results: We found that sEVs derived from ESFT cells contained nearly 2-fold more number of unique individual miRNAs as compared to non-ESFT samples. Quantitative analysis of the differential enrichment of sEV miRNAs resulted in the identification of 62 sEV-miRNAs (exo-miRNAs) with significant (P < .05) enrichment variation between ESFT and non-ESFT sEV samples. To determine if we could utilize this miRNA signature to diagnose ESFT patients via a liquid biopsy, we analyzed the RNA content of total circulating sEVs isolated from 500 µL plasma from 5 pediatric ESFT patients, 2 pediatric osteosarcoma patients, 2 pediatric rhabdomyosarcoma patients, and 4 non-cancer pediatric controls. Pearson's clustering of 60 of the 62 candidate exo-miRNAs correctly identified 80% (4 of 5) of pathology confirmed ESFT patients. Importantly, RNAseq analysis of tumor tissue from the 1 outlier, revealed a previously uncharacterized EWS-FLI1 translocation.Conclusions: Taken together, these findings support the development and validation of an exo-miRNA-based liquid biopsy to aid in the diagnosis and monitoring of ESFT.

Confounding Factors Impacting microRNA Expression in Human Saliva: Methodological and Biological Considerations

There is growing interest in saliva microRNAs (miRNAs) as non-invasive biomarkers for human disease. Such an approach requires understanding how differences in experimental design affect miRNA expression. Variations in technical methodologies, coupled with inter-individual variability may reduce study reproducibility and generalizability. Another barrier facing salivary miRNA biomarker research is a lack of recognized “control miRNAs”. In one of the largest studies of human salivary miRNA to date (922 healthy individuals), we utilized 1225 saliva samples to quantify variability in miRNA expression resulting from aligner selection (Bowtie1 vs. Bowtie2), saliva collection method (expectorated vs. swabbed), RNA stabilizer (presence vs. absence), and individual biological factors (sex, age, body mass index, exercise, caloric intake). Differential expression analyses revealed that absence of RNA stabilizer introduced the greatest variability, followed by differences in methods of collection and aligner. Biological factors generally affected a smaller number of miRNAs. We also reported coefficients of variations for 643 miRNAs consistently present in saliva, highlighting several salivary miRNAs to serve as reference genes. Thus, the results of this analysis can be used by researchers to optimize parameters of salivary miRNA measurement, exclude miRNAs confounded by numerous biologic factors, and identify appropriate miRNA controls.

miR-381-3p Inhibits Intramuscular Fat Deposition through Targeting FABP3 by ceRNA Regulatory Network

Intramuscular fat (IMF) deposition is an important determinant of pork quality and a complex process facilitated by non-coding ceRNAs. In this study, 52 Berkshire × Anqing Sixwhite crossbred pigs were slaughtered to measure eight carcass and pork quality traits. Whole-transcriptome sequencing analysis was performed using longissimus dorsi samples of six low- and high-IMF samples; 34 ceRNA networks, based on 881, 394, 158 differentially expressed (DE) lncRNAs, miRNAs, and mRNAs, were constructed. Following weighted gene co-expression network analysis between the low and high IMF, only one ceRNA, lncRNA4789/miR-381-3p/FABP3, that showed similar DE trend in longissimus dorsi tissue was retained. Dual-luciferase reporter assays further indicated that FABP3 was a direct, functional target of miR-381-3p, where miR-381-3p overexpression inhibited the mRNA and protein expression of FABP3. In addition, overexpressed lncRNA4789 attenuated the effect of miR-381-3p on FABP3 by sponging miR-381-3p. Cell function verification experiment demonstrated that miR-381-3p suppressed IMF deposition by inhibiting preadipocyte cell differentiation and lipid droplet deposition via the suppression of FABP3 expression in the peroxisome proliferator-activated receptor signalling pathway, whereas lncRNA4789 rescued FABP3 expression by sponging miR-381-3p. Our study may aid in identifying novel molecular markers for its optimization in IMF which is of importance in breeding for improving pork quality.

Urinary extracellular vesicles and micro-RNA as markers of acute kidney injury after cardiac surgery

We hypothesised that measuring changes in urinary levels of EV and miR will predict the onset of acute kidney injury in cardiac surgery patients. The study was performed in the cohort of the REVAKI-2 trial. Urine samples were collected before and 24 h after the procedure from 94 cardiac surgery patients. Urinary particle concentrations and size distribution were assessed using NanoSight. EV derivation and levels were measured using flow cytometry. Samples from 10 selected patients were sequenced, and verification was performed with advanced TaqMan assays in samples from all patients. Urinary particle concentrations significantly increased in patients with AKI after surgery, with the percentage of EV positive for CD105 and β1-integrin also increasing. Pre-surgery podocalyxin-positive EV were significantly lower in patients with AKI. Their levels correlated with the severity of the injury. Pre-operative miR-125a-5p was expressed at lower levels in urine from patients with AKI when adjusted for urinary creatinine. Levels of miR-10a-5p were lower after surgery in AKI patients and its levels correlated with the severity of the injury. Pre-operative levels of podocalyxin EVs, urinary particle concentrations and miR-125a-5p had moderate AKI predictive value and, in a logistic model together with ICU lactate levels, offered good (AUC = 82%) AKI prediction.

miR-29a-5p Inhibits Prenatal Hair Placode Formation Through Targeting EDAR by ceRNA Regulatory Network

Hair placode formation is an important stage of hair follicle morphogenesis and it is a complex process facilitated by non-coding RNAs. In this study, we conducted whole transcriptome sequencing analysis of skin, heart, liver, lung, and kidney tissues of day 41 (E41) normal and hairless pig embryos, and respectively detected 15, 8, and 515 skin-specific differentially expressed (DE) lncRNAs, miRNAs, and mRNAs. Furthermore, 18 competing endogenous RNA (ceRNA) networks were constructed. Following weighted gene co-expression network analysis (WGCNA) of stages E39, E41, E45, E52, and E60, between normal and hairless pig embryos, only two ceRNAs (lncRNA2162.1/miR-29a-5p/BMPR1b and lncRNA627.1/miR-29a-5p/EDAR) that showed period-specific differential expression in E41 skin were retained. Dual-luciferase reporter assays further indicated that EDAR was a direct, functioning target of miR-29a-5p and that no binding site was found in BMPR1b. Moreover, miR-29a-5p overexpression inhibited the mRNA and protein expression of EDAR while no significant differential expression of BMPR1b was detected. In addition, over-expressed lncRNA627.1 reduces the expression of miR-29a-5p and increase EDAR expression while inhibits lncRNA627.1 resulted in a opposite expression trend. Cell proliferation result demonstrated that lower expression of EDAR and lncRNA627.1 inhibited hair placode precursor cells (HPPCs) proliferation in a manner similar to that shown by over-expressed miR-29a-5p. This study identified that miR-29a-5p inhibited HPPCs proliferation via the suppression of EDAR expression in the EDA/EDAR signaling pathway, while lncRNA627.1 rescues EDAR expression. Our study provides a basis for a better understanding of the mechanisms underlying the ceRNA complex, miR29a-5p/EDAR/lncRNA627.1, that could regulate hair placode formation, which may help decipher diseases affecting human hair.

The EGFR Signaling Modulates in Mesenchymal Stem Cells the Expression of miRNAs Involved in the Interaction with Breast Cancer Cells

We previously demonstrated that the epidermal growth factor receptor (EGFR) modulates in mesenchymal stem cells (MSCs) the expression of a number of genes coding for secreted proteins that promote breast cancer progression. However, the role of the EGFR in modulating in MSCs the expression of miRNAs potentially involved in the progression of breast cancer remains largely unexplored. Following small RNA-sequencing, we identified 36 miRNAs differentially expressed between MSCs untreated or treated with the EGFR ligand transforming growth factor α (TGFα), with a fold change (FC) < 0.56 or FC ≥ 1.90 (CI, 95%). KEGG analysis revealed a significant enrichment in signaling pathways involved in cancer development and progression. EGFR activation in MSCs downregulated the expression of different miRNAs, including miR-23c. EGFR signaling also reduced the secretion of miR-23c in conditioned medium from MSCs. Functional assays demonstrated that miR-23c acts as tumor suppressor in basal/claudin-low MDA-MB-231 and MDA-MB-468 cells, through the repression of IL-6R. MiR-23c downregulation promoted cell proliferation, migration and invasion of these breast cancer cell lines. Collectively, our data suggested that the EGFR signaling regulates in MSCs the expression of miRNAs that might be involved in breast cancer progression, providing novel information on the mechanisms that regulate the MSC-tumor cell cross-talk.

Differential microRNA profiles of intramuscular and secreted extracellular vesicles in human tissue-engineered muscle

Exercise affects the expression of microRNAs (miR/s) and muscle-derived extracellular vesicles (EVs). To evaluate sarcoplasmic and secreted miR expression in human skeletal muscle in response to exercise-mimetic contractile activity, we utilized a three-dimensional tissue-engineered model of human skeletal muscle ("myobundles"). Myobundles were subjected to three culture conditions: no electrical stimulation (CTL), chronic low frequency stimulation (CLFS), or intermittent high frequency stimulation (IHFS) for 7 days. RNA was isolated from myobundles and from extracellular vesicles (EVs) secreted by myobundles into culture media; miR abundance was analyzed by miRNA-sequencing. We used edgeR and a within-sample design to evaluate differential miR expression and Pearson correlation to evaluate correlations between myobundle and EV populations within treatments with statistical significance set at p < 0.05. Numerous miRs were differentially expressed between myobundles and EVs; 116 miRs were differentially expressed within CTL, 3 within CLFS, and 2 within IHFS. Additionally, 25 miRs were significantly correlated (18 in CTL, 5 in CLFS, 2 in IHFS) between myobundles and EVs. Electrical stimulation resulted in differential expression of 8 miRs in myobundles and only 1 miR in EVs. Several KEGG pathways, known to play a role in regulation of skeletal muscle, were enriched, with differentially overrepresented miRs between myobundle and EV populations identified using miEAA. Together, these results demonstrate that in vitro exercise-mimetic contractile activity of human engineered muscle affects both their expression of miRs and number of secreted EVs. These results also identify novel miRs of interest for future studies of the role of exercise in organ-organ interactions in vivo.

Electroacupuncture altered expression of microRNAs in Stat5 knockout obese mice

BACKGROUND

Increasing evidence shows that miRNAs contribute to the establishment and development of obesity by affecting many biological and pathological processes, such as adipocyte differentiation, hepatic lipid metabolism, insulin resistance, and neurological regulation of obesity. As a clinical intervention approach, acupuncture has been shown to be effective in the treatment of obesity and other metabolic diseases. Our previous whole genome study in central nervous system (CNS)-specific Stat5 knockout (NKO) obese mice found that electroacupuncture (EA) could reduce body weight and promote white browning.

OBJECTIVE

To clarify the effect of EA on miRNAs and understand how it regulates gene expression.

METHODS

Twelve-week-old male Stat5NKO mice with body weight 20% greater than that of Stat5fl/fl (control) mice were divided into a Stat5NKO (model) group and EA-treated Stat5NKO + EA group. A cohort of Stat5fl/fl mice of the same age were included as the control group. EA was administered under isoflurane anesthesia at unilateral ST36 and ST44 daily (left and right sides were treated every other day), 6 times per week for a total of 4 weeks. The miRNA profile was generated and miRNA regulatory networks were analyzed in the Stat5 nestin-cre mice before and after EA treatment. Autophagy-related proteins in adipocytes were detected after over-expression of miR27a.

RESULTS

EA altered abnormal miRNA expression, including miRNA27a expression, and reduced the autophagy-related proteins ATG5 and ATG12.

CONCLUSION

We found that EA could regulate miRNA27a-mediated autophagy-related proteins and promote white fat browning, which may contribute to weight loss. To our knowledge, this is the first report of miRNAs potentially driving the effect of EA on white fat browning through the autophagy process.

Discriminating Neoplastic from Nonneoplastic Tissues Using an miRNA-Based Deep Cancer Classifier

Next-generation sequencing has enabled the collection of large biological data sets, allowing novel molecular-based classification methods to be developed for increased understanding of disease. miRNAs are small regulatory RNA molecules that can be quantified using next-generation sequencing and are excellent classificatory markers. Herein, we adapt a deep cancer classifier (DCC) to differentiate neoplastic from nonneoplastic samples using comprehensive miRNA expression profiles from 1031 human breast and skin tissue samples. The classifier was fine-tuned and evaluated using 750 neoplastic and 281 nonneoplastic breast and skin tissue samples. Performance of the DCC was compared with two machine-learning classifiers: support vector machine and random forests. In addition, performance of feature extraction through the DCC was also compared with a developed feature selection algorithm, cancer specificity. The DCC had the highest performance of area under the receiver operating curve and high performance in both sensitivity and specificity, unlike machine-learning and feature selection models, which often performed well in one metric compared with the other. In particular, deep learning was shown to have noticeable advantages with highly heterogeneous data sets. In addition, our cancer specificity algorithm identified candidate biomarkers for differentiating neoplastic and nonneoplastic tissue samples (eg, miR-144 and miR-375 in breast cancer and miR-375 and miR-451 in skin cancer).

Identification of Exosomal MicroRNA Signature by Liquid Biopsy in Hereditary Hemorrhagic Telangiectasia Patients

Hereditary hemorrhagic telangiectasia (HHT) is a rare autosomal dominant vascular dysplasia characterized by epistaxis, mucocutaneous telangiectases, and arteriovenous malformations (AVM) in the visceral organs. The diagnosis of HHT is based on clinical Curaçao criteria, which show limited sensitivity in children and young patients. Here, we carried out a liquid biopsy by which we isolated total RNA from plasma exosome samples. A cohort of 15 HHT type 1 patients, 15 HHT type 2 patients, and 10 healthy relatives were analyzed. Upon gene expression data processing and normalization, a statistical analysis was performed to explore similarities in microRNA expression patterns among samples and detect differentially expressed microRNAs between HHT samples and the control group. We found a disease-associated molecular fingerprint of 35 miRNAs over-represented in HHT vs. controls, with eight being specific for HHT1 and 11 for HHT2; we also found 30 under-represented, including nine distinct for HHT1 and nine for HHT2. The analysis of the receiver operating characteristic (ROC) curves showed that eight miRNAs had good (AUC > 75%) or excellent (AUC > 90%) diagnosis value for HHT and even for type HHT1 and HHT2. In addition, we identified the cellular origin of these miRNAs among the cell types involved in the vascular malformations. Interestingly, we found that only some of them were incorporated into exosomes, which suggests a key functional role of these exosomal miRNAs in the pathophysiology of HHT.

Analysis of the Applicability of microRNAs in Peripheral Blood Leukocytes as Biomarkers of Sensitivity and Exposure to Fractionated Radiotherapy towards Breast Cancer

Biomarkers for predicting individual response to radiation and for dose verification are needed to improve radiotherapy. A biomarker should optimally show signal fidelity, meaning that its level is stable and proportional to the absorbed dose. miRNA levels in human blood serum were suggested as promising biomarkers. The aim of the present investigation was to test the miRNA biomarker in leukocytes of breast cancer patients undergoing external beam radiotherapy. Leukocytes were isolated from blood samples collected prior to exposure (control); on the day when a total dose of 2 Gy, 10 Gy, or 20 Gy was reached; and one month after therapy ended (46-50 Gy in total). RNA sequencing was performed and univariate analysis was used to analyse the effect of the radiation dose on the expression of single miRNAs. To check if combinations of miRNAs can predict absorbed dose, a multinomial logistic regression model was built using a training set from eight patients (representing 40 samples) and a validation set with samples from the remaining eight patients (15 samples). Finally, Broadside, an explorative interaction mining tool, was used to extract sets of interacting miRNAs. The most prominently increased miRNA was miR-744-5p, followed by miR-4461, miR-34a-5p, miR-6513-5p, miR-1246, and miR-454-3p. Decreased miRNAs were miR-3065-3p, miR-103a-2-5p, miR-30b-3p, and miR-5690. Generally, most miRNAs showed a relatively strong inter-individual variability and different temporal patterns over the course of radiotherapy. In conclusion, miR-744-5p shows promise as a stable miRNA marker, but most tested miRNAs displayed individual signal variability which, at least in this setting, may exclude them as sensitive biomarkers of radiation response.

Temporal Gene Expression Signature of Plasma Extracellular Vesicles-MicroRNAs from Post-Smolt Coho Salmon Challenged with Piscirickettsia salmonis

Piscirickettsiosis is the most important bacterial disease in the Chilean salmon industry, which has borne major economic losses due to failure to control it. Cells use extracellular vesicles (EVs) as an inter-cellular communicators to deliver several factors (e.g., microRNAs) that may regulate the responses of other cells. However, there is limited knowledge about the identification and characterization of EV-miRNAs in salmonids or the effect of infections on these. In this study, Illumina sequencing technology was used to identify Coho salmon plasma EV-miRNAs upon Piscirickettsia salmonis infection at four different time points. A total of 118 novels and 188 known EV-miRNAs, including key immune teleost miRNAs families (e.g., miR-146, miR-122), were identified. A total of 245 EV-miRNAs were detected as differentially expressed (FDR < 5%) in terms of control, with a clear down-regulation pattern throughout the disease. KEGG enrichment results of EV-miRNAs target genes showed that they were grouped mainly in cellular, stress, inflammation and immune responses. Therefore, it is hypothesized that P. salmonis could potentially benefit from unbalanced modulation response of Coho salmon EV-miRNAs in order to promote a hyper-inflammatory and compromised immune response through the suppression of different key immune host miRNAs during the course of the infection, as indicated by the results of this study.

Seasonal Variation in miR-328-3p and let-7d-3p Are Associated With Seasonal Allergies and Asthma Symptoms in Children

OBJECTIVE

MicroRNAs (miRs) are small non-coding RNA molecules of around 18-22 nucleotides that are key regulators of many biologic processes, particularly inflammation. The purpose of this study was to determine the association of circulating miRs from asthmatic children with seasonal variation in allergic inflammation and asthma symptoms.

METHODS

We used available small RNA sequencing on blood serum from 398 children with mild-to-moderate asthma from the Childhood Asthma Management Program. We used seasonal asthma symptom data at the study baseline and allergen affection status from baseline skin prick tests as primary outcomes. We identified differentially expressed (DE) miRs between pairs of seasons using DESeq2. Regression analysis was used to identify associations between allergy status to specific seasonal allergens and DE miRs in 4 seasons and between seasonal asthma symptom data and DE miRs. We performed pathway enrichment analysis for target genes of the DE miRs using DAVID.

RESULTS

After quality control, 398 samples underwent differential analysis between the 4 seasons. We found 52 unique miRs from a total of 81 DE miRs across seasons. Further investigation of the association between these miRs and sensitization to seasonal allergens using skin prick tests revealed that 26 unique miRs from a total of 38 miRs were significantly associated with a same-season allergen. Comparison between seasonal asthma symptom data revealed that 2 of these 26 miRs also had significant associations with asthma symptoms in the same seasons: miR-328-3p (P < 0.03) and let-7d-3p (P < 0.05). Enrichment analysis showed that the most enriched pathway clusters were Rap1, Ras, and MAPK signaling pathways.

CONCLUSION

Our results show seasonal variation in miR-328-3p and let-7d-3p are significantly associated with seasonal asthma symptoms and seasonal allergies. These indicate a potentially protective role for let-7d-3p and a deleterious role for miR-328-3p in asthmatics sensitized to mulberry. Further work will determine whether these miRs are drivers or results of the allergic response.

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.

Huntington's disease mice and human brain tissue exhibit increased G3BP1 granules and TDP43 mislocalization

Chronic cellular stress associated with neurodegenerative disease can result in the persistence of stress granule (SG) structures, membraneless organelles that form in response to cellular stress. In Huntington's disease (HD), chronic expression of mutant huntingtin generates various forms of cellular stress, including activation of the unfolded protein response and oxidative stress. However, it has yet to be determined whether SGs are a feature of HD neuropathology. We examined the miRNA composition of extracellular vesicles (EVs) present in the cerebrospinal fluid (CSF) of patients with HD and show that a subset of their target mRNAs were differentially expressed in the prefrontal cortex. Of these targets, SG components were enriched, including the SG-nucleating Ras GTPase-activating protein-binding protein 1 (G3BP1). We investigated localization and levels of G3BP1 and found a significant increase in the density of G3BP1-positive granules in the cortex and hippocampus of R6/2 transgenic mice and in the superior frontal cortex of the brains of patients with HD. Intriguingly, we also observed that the SG-associated TAR DNA-binding protein 43 (TDP43), a nuclear RNA/DNA binding protein, was mislocalized to the cytoplasm of G3BP1 granule-positive HD cortical neurons. These findings suggest that G3BP1 SG dynamics may play a role in the pathophysiology of HD.

HDL-enriched miR-30a-5p is associated with HDL-cholesterol levels and glucose metabolism in healthy men and women

Aim: To investigate the associations between high-density lipoprotein (HDL)-enriched miRNAs and the cardiometabolic profile of healthy men and women. Patients & methods: miRNAs were quantified using next-generation sequencing of miRNAs extracted from purified HDL and plasma from 17 healthy men and women couples. Results: Among the HDL-enriched miRNAs, miR-30a-5p correlated positively with HDL-cholesterol levels, whereas miR-144-5p and miR-30a-5p were negatively associated with fasting insulin levels and Homeostasis model assessment of insulin resistance index. Overall, miR-30a-5p, miR-150-5p and sex contributed to 45% of HDL-cholesterol variance. A model containing only miR-30a-5p, age and sex explained 41% of fasting glucose variance. Conclusion: HDL-enriched miRNAs, notably miR-30a-5p, are associated with cardiometabolic markers. These miRNAs could play a role in HDL's protective functions, particularly regarding glucose-insulin homeostasis.

EZcount: An all-in-one software for microRNA expression quantification from NGS sequencing data

MicroRNAs (miRNAs) are short endogenous molecules of RNA that influence cell regulation by suppressing genes. Their ubiquity throughout all branches of the tree of life has suggested their central role in many cellular functions. Nowadays, several personalized medicine applications rely on miRNAs as biomarkers for diagnoses, prognoses, and prediction of drug response. The increasing ease of sequencing miRNAs contrasts with the difficulty of accurately quantifying their concentration. The use of general purpose aligners is only a partial solution as they have limited possibilities to accurately solve ambiguous mapping due to the short length of these sequences. We developed EZcount, an all-in-one software that, with a single command, performs the entire quantification process: from raw fastq files to read counts. Experiments show that EZcount is more sensitive and accurate than methods based on sequence alignment, independently of the library preparation protocol and sequencing machine. The parallel architecture of EZcount makes it fast enough to process a sample in minutes using a standard workstation. EZcount runs on all of the most common operating systems (Linux, Windows and MacOS) and is freely available for download at https://gitlab.com/BioAlgo/miR-pipe. A detailed description of the datasets, the raw experimental results, and all the scripts used for testing are available as supplementary material.

An Information Theoretical Multilayer Network Approach to Breast Cancer Transcriptional Regulation

Breast cancer is a complex, highly heterogeneous disease at multiple levels ranging from its genetic origins and molecular processes to clinical manifestations. This heterogeneity has given rise to the so-called intrinsic or molecular breast cancer subtypes. Aside from classification, these subtypes have set a basis for differential prognosis and treatment. Multiple regulatory mechanisms-involving a variety of biomolecular entities-suffer from alterations leading to the diseased phenotypes. Information theoretical approaches have been found to be useful in the description of these complex regulatory programs. In this work, we identified the interactions occurring between three main mechanisms of regulation of the gene expression program: transcription factor regulation, regulation via noncoding RNA, and epigenetic regulation through DNA methylation. Using data from The Cancer Genome Atlas, we inferred probabilistic multilayer networks, identifying key regulatory circuits able to (partially) explain the alterations that lead from a healthy phenotype to different manifestations of breast cancer, as captured by its molecular subtype classification. We also found some general trends in the topology of the multi-omic regulatory networks: Tumor subtype networks present longer shortest paths than their normal tissue counterpart; epigenomic regulation has frequently focused on genes enriched for certain biological processes; CpG methylation and miRNA interactions are often part of a regulatory core of conserved interactions. The use of probabilistic measures to infer information regarding theoretical-derived multilayer networks based on multi-omic high-throughput data is hence presented as a useful methodological approach to capture some of the molecular heterogeneity behind regulatory phenomena in breast cancer, and potentially other diseases.

A comparative analysis of heart microRNAs in vertebrates brings novel insights into the evolution of genetic regulatory networks

BACKGROUND

During vertebrate evolution, the heart has undergone remarkable changes that lead to morphophysiological differences in the fully formed heart of these species, such as chamber septation, heart rate frequency, blood pressure, and cardiac output volume. Despite these differences, the heart developmental process is guided by a core gene set conserved across vertebrates. Nonetheless, the regulatory mechanisms controlling the expression of genes involved in heart development and maintenance are largely uncharted. MicroRNAs (miRNAs) have been described as important regulatory elements in several biological processes, including heart biology. These small RNA molecules are broadly conserved in sequence and genomic context in metazoans. Mutations may occur in miRNAs and/or genes that contribute to the establishment of distinct repertoires of miRNA-target interactions, thereby favoring the differential control of gene expression and, consequently, the origin of novel phenotypes. In fact, several studies showed that miRNAs are integrated into genetic regulatory networks (GRNs) governing specific developmental programs and diseases. However, studies integrating miRNAs in vertebrate heart GRNs under an evolutionary perspective are still scarce.

RESULTS

We comprehensively examined and compared the heart miRNome of 20 species representatives of the five major vertebrate groups. We found 54 miRNA families with conserved expression and a variable number of miRNA families with group-specific expression in fishes, amphibians, reptiles, birds, and mammals. We also detected that conserved miRNAs present higher expression levels and a higher number of targets, whereas the group-specific miRNAs present lower expression levels and few targets.

CONCLUSIONS

Both the conserved and group-specific miRNAs can be considered modulators orchestrating the core and peripheral genes of heart GRNs of vertebrates, which can be related to the morphophysiological differences and similarities existing in the heart of distinct vertebrate groups. We propose a hypothesis to explain evolutionary differences in the putative functional roles of miRNAs in the heart GRNs analyzed. Furthermore, we present new insights into the molecular mechanisms that could be helping modulate the diversity of morphophysiology in the heart organ of vertebrate species.

Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

Altered non-coding RNA expression profile in F1 progeny 1 year after parental irradiation is linked to adverse effects in zebrafish

Gamma radiation produces DNA instability and impaired phenotype. Previously, we observed negative effects on phenotype, DNA methylation, and gene expression profiles, in offspring of zebrafish exposed to gamma radiation during gametogenesis. We hypothesize that previously observed effects are accompanied with changes in the expression profile of non-coding RNAs, inherited by next generations. Non-coding RNA expression profile was analysed in F₁ offspring (5.5 h post-fertilization) by high-throughput sequencing 1 year after parental irradiation (8.7 mGy/h, 5.2 Gy total dose). Using our previous F₁-γ genome-wide gene expression data (GSE98539), hundreds of mRNAs were predicted as targets of differentially expressed (DE) miRNAs, involved in pathways such as insulin receptor, NFkB and PTEN signalling, linking to apoptosis and cancer. snRNAs belonging to the five major spliceosomal snRNAs were down-regulated in the F₁-γ group, Indicating transcriptional and post-transcriptional alterations. In addition, DEpiRNA clusters were associated to 9 transposable elements (TEs) (LTR, LINE, and TIR) (p = 0.0024), probable as a response to the activation of these TEs. Moreover, the expression of the lincRNAs malat-1, and several others was altered in the offspring F₁, in concordance with previously observed phenotypical alterations. In conclusion, our results demonstrate diverse gamma radiation-induced alterations in the ncRNA profiles of F₁ offspring observable 1 year after parental irradiation.

MicroRNA Profiles in Monocyte-Derived Macrophages Generated by Interleukin-27 and Human Serum: Identification of a Novel HIV-Inhibiting and Autophagy-Inducing MicroRNA

Interleukin-27 (IL-27) is a pleiotropic cytokine that influences the innate and adaptive immune systems. It inhibits viral infection and regulates the expression of microRNAs (miRNAs). We recently reported that macrophages differentiated from human primary monocytes in the presence of IL-27 and human AB serum resisted human immunodeficiency virus (HIV) infection and showed significant autophagy induction. In the current study, the miRNA profiles in these cells were investigated, especially focusing on the identification of novel miRNAs regulated by IL-27-treatment. The miRNA sequencing analysis detected 38 novel miRNAs. Real-time reverse transcription polymerase chain reaction (RT-PCR) analysis confirmed that IL-27 differentially regulated the expression of 16 of the 38 miRNAs. Overexpression of the synthesized miRNA mimics by transfection revealed that miRAB40 had potent HIV-inhibiting and autophagy-inducing properties. B18R, an interferon (IFN)-neutralization protein, partially suppressed both activities, indicating that the two functions were induced via IFN-dependent and -independent pathways. Although the target mRNA(s) of miRAB40 involving in the induction of both functions was unable to identify in this study, the discovery of miRAB40, a potential HIV-inhibiting and autophagy inducing miRNA, may provide novel insights into the miRNA (small none-coding RNA)-mediated regulation of HIV inhibition and autophagy induction as an innate immune response.

High Throughput Sequencing of MicroRNA in Rainbow Trout Plasma, Mucus, and Surrounding Water Following Acute Stress

Circulating plasma microRNAs (miRNAs) are well established as biomarkers of several diseases in humans and have recently been used as indicators of environmental exposures in fish. However, the role of plasma miRNAs in regulating acute stress responses in fish is largely unknown. Tissue and plasma miRNAs have recently been associated with excreted miRNAs; however, external miRNAs have never been measured in fish. The objective of this study was to identify the altered plasma miRNAs in response to acute stress in rainbow trout (Oncorhynchus mykiss), as well as altered miRNAs in fish epidermal mucus and the surrounding ambient water. Small RNA was extracted and sequenced from plasma, mucus, and water collected from rainbow trout pre- and 1 h-post a 3-min air stressor. Following small RNA-Seq and pathway analysis, we identified differentially expressed plasma miRNAs that targeted biosynthetic, degradation, and metabolic pathways. We successfully isolated miRNA from trout mucus and the surrounding water and detected differences in miRNA expression 1-h post air stress. The expressed miRNA profiles in mucus and water were different from the altered plasma miRNA profile, which indicated that the plasma miRNA response was not associated with or immediately reflected in external samples, which was further validated through qPCR. This research expands understanding of the role of plasma miRNA in the acute stress response of fish and is the first report of successful isolation and profiling of miRNA from fish mucus or samples of ambient water. Measurements of miRNA from plasma, mucus, or water can be further studied and have potential to be applied as non-lethal indicators of acute stress in fish.

The Identification of Plasma Exosomal miR-423-3p as a Potential Predictive Biomarker for Prostate Cancer Castration-Resistance Development by Plasma Exosomal miRNA Sequencing

Castration-resistant prostate cancer (CRPC) is the major cause of death from prostate cancer. Biomarkers to improve early detection and prediction of CRPC especially using non-invasive liquid biopsies could improve outcomes. Therefore, we investigated the plasma exosomal miRNAs associated with CRPC and their potential for development into non-invasive early detection biomarkers for resistance to treatment. RNA-sequencing, which generated approximately five million reads per patient, was performed to identify differentially expressed plasma exosomal miRNAs in 24 treatment-naive prostate cancer and 24 CRPC patients. RT-qPCR was used to confirm the differential expressions of six exosomal miRNAs, miR-423-3p, miR-320a, miR-99a-5p, miR-320d, miR-320b, and miR-150-5p (p = 7.3 × 10^-8, 0.0020, 0.018, 0.0028, 0.0013, and 0.0058, respectively) firstly in a validation cohort of 108 treatment-naive prostate cancer and 42 CRPC patients. The most significant differentially expressed miRNA, miR-423-3p, was shown to be associated with CRPC with area under the ROC curve (AUC) = 0.784. Combining miR-423-3p with prostate-specific antigen (PSA) enhanced the prediction of CRPC (AUC = 0.908). A separate research center validation with 30 treatment-naive and 30 CRPC patients also confirmed the differential expression of miR-423-3p (p = 0.016). Finally, plasma exosomal miR-423-3p expression in CRPC patients was compared to 36 non-CRPC patients under androgen depletion therapy, which showed significantly higher expression in CRPC than treated non-CRPC patients (p < 0.0001) with AUC = 0.879 to predict CRPC with no difference between treatment-naive and treated non-CRPC patients. Therefore, our findings demonstrate that a number of plasma exosomal miRNAs are associated with CRPC and miR-423-3p may serve as a biomarker for early detection/prediction of castration-resistance.

Network- and systems-based re-engineering of dendritic cells with non-coding RNAs for cancer immunotherapy

Dendritic cells (DCs) are professional antigen-presenting cells that induce and regulate adaptive immunity by presenting antigens to T cells. Due to their coordinative role in adaptive immune responses, DCs have been used as cell-based therapeutic vaccination against cancer. The capacity of DCs to induce a therapeutic immune response can be enhanced by re-wiring of cellular signalling pathways with microRNAs (miRNAs). Methods: Since the activation and maturation of DCs is controlled by an interconnected signalling network, we deploy an approach that combines RNA sequencing data and systems biology methods to delineate miRNA-based strategies that enhance DC-elicited immune responses. Results: Through RNA sequencing of IKKβ-matured DCs that are currently being tested in a clinical trial on therapeutic anti-cancer vaccination, we identified 44 differentially expressed miRNAs. According to a network analysis, most of these miRNAs regulate targets that are linked to immune pathways, such as cytokine and interleukin signalling. We employed a network topology-oriented scoring model to rank the miRNAs, analysed their impact on immunogenic potency of DCs, and identified dozens of promising miRNA candidates, with miR-15a and miR-16 as the top ones. The results of our analysis are presented in a database that constitutes a tool to identify DC-relevant miRNA-gene interactions with therapeutic potential (https://www.synmirapy.net/dc-optimization). Conclusions: Our approach enables the systematic analysis and identification of functional miRNA-gene interactions that can be experimentally tested for improving DC immunogenic potency.

mirPLS: a partial linear structure identifier method for cancer subtyping using microRNAs

MOTIVATION

MicroRNAs (miRNAs) are small non-coding RNAs that have been successfully identified to be differentially expressed in various cancers. However, some miRNAs were reported to be up-regulated in one subtype of a cancer but down-regulated in another, making overall associations between these miRNAs and the heterogeneous cancer non-linear. These non-linearly associated miRNAs, if identified, are thus informative for cancer subtyping.

RESULTS

Here, we propose mirPLS, a Partial Linear Structure identifier for miRNA data that simultaneously identifies miRNAs of linear or non-linear associations with cancer status when non-linearly associated miRNAs can then be used for subsequent cancer subtyping. Simulation studies showed that mirPLS can identify both non-linearly and linearly outcome-associated miRNAs more accurately than the comparison methods. Using the identified non-linearly associated miRNAs much improves the cancer subtyping accuracy. Applications to miRNA data of three different cancer types suggest that the cancer subtypes defined by the non-linearly associated miRNAs identified by mirPLS are consistently more predictive of patient survival and more biological meaningful.

AVAILABILITY AND IMPLEMENTATION

The R package mirPLS is available for downloading from https://github.com/pfruan/mirPLS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Circulating miRNA Increases the Diagnostic Accuracy of Chromogranin A in Metastatic Pancreatic Neuroendocrine Tumors

Simple Summary

Despite its varying sensitivity and decreased specificity, chromogranin A (CgA) is the most widely used biomarker for neuroendocrine tumors. The most common factor affecting its diagnostic accuracy is the use of proton pump inhibitors (PPIs). Our aim was to investigate circulating miRNA expression profiles in pancreatic neuroendocrine tumors (pNET) and pheochromocytomas/paragangliomas (PPGL) to find miRNAs which could be used as biomarkers along with CgA in these patients. MiRNA expression profiles were determined with next generation sequencing and validated by quantitative real time PCR in 74 samples obtained from patients and healthy volunteers treated with PPI. We observed a global downregulation of miRNAs in NET compared to controls. A set of miRNAs in combination with CgA resulted in the best discrimination of pNET irrespective of PPI treatment and a combination of miRNAs increased the diagnostic utility of CgA even in pNET patients with low CgA.

Abstract

Chromogranin A (CgA) is the most widely accepted biomarker for neuroendocrine tumors (NET) but its diagnostic accuracy is dependent on tumor type and the use of proton-pump inhibitors (PPI). We investigated the diagnostic value of circulating miRNAs along with CgA in pancreatic neuroendocrine tumors (pNET). 74 serum samples from patients with pNET (n = 25, nonfunctioning), pheochromocytoma/paraganglioma (PPGL, n = 20), healthy individuals with normal CgA (n = 29) including 10 samples from 5 healthy individuals with and without current PPI treatment were collected. MiRNA expression profiles were determined using next-generation sequencing, followed by validation with individual TaqMan assays. A global downregulation of miRNAs was observed in patients with NET compared to controls. MiRNA expression of 33 miRNAs was able to discriminate tumor samples from controls. No miRNA alone could be considered as an applicable biomarker for pNET or PPGL. However, using a logistic model, the combination of a set of miRNAs increased the discriminatory role of CgA irrespective of PPI treatment. In pNET patients with normal CgA level our regression model yielded high (89.4%) diagnostic accuracy (AUC: 0.904, sensitivity: 66.6%, specificity: 96.5%). A set of miRNAs increased the diagnostic utility of CgA in pNET even in patients with low CgA.

cfRNAs as biomarkers in oncology - still experimental or applied tool for personalized medicine already?

Currently, the challenges of contemporary oncology are focused mainly on the development of personalized medicine and precise treatment, which could be achieved through the use of molecular biomarkers. One of the biological molecules with great potential are circulating free RNAs (cfRNAs) which are present in various types of body fluids, such as blood, serum, plasma, and saliva. Also, different types of cfRNA particles can be distinguished depending on their length and function: microRNA (miRNA), PIWI-interacting RNA (piRNA), tRNA-derived RNA fragments (tRFs), circular RNA (circRNA), long non-coding RNA (lncRNA), and messenger RNA (mRNA). Moreover, cfRNAs occur in various forms: as a free molecule alone, in membrane vesicles, such as exosomes, or in complexes with proteins and lipids. One of the modern approaches for monitoring patient's condition is a "liquid biopsy" that provides a non-invasive and easily available source of circulating RNAs. Both the presence of specific cfRNA types as well as their concentration are dependent on many factors including cancer type or even reaction to treatment. Despite the possibility of using circulating free RNAs as biomarkers, there is still a lack of validated diagnostic panels, defined protocols for sampling, storing as well as detection methods. In this work we examine different types of cfRNAs, evaluate them as possible biomarkers, and analyze methods of their detection. We believe that further research on cfRNA and defining diagnostic panels could lead to better and faster cancer identification and improve treatment monitoring.

Neonatal Diet Impacts Circulatory miRNA Profile in a Porcine Model

microRNAs (miRNAs) are conserved non-coding small nucleotide molecules found in nearly all species and breastmilk. miRNAs present in breastmilk are very stable to freeze-thaw, RNase treatment, and low pH as they are protected inside exosomes. They are involved in regulating several physiologic and pathologic processes, including immunologic pathways, and we have demonstrated better immune response to vaccines in piglets fed with human milk (HM) in comparison to dairy-based formula (MF). To understand if neonatal diet impacts circulatory miRNA expression, serum miRNA expression was evaluated in piglets fed HM or MF while on their neonatal diet at postnatal day (PND) 21 and post-weaning to solid diet at PND 35 and 51. MF fed piglets showed increased expression of 14 miRNAs and decreased expression of 10 miRNAs, relative to HM fed piglets at PND 21. At PND 35, 9 miRNAs were downregulated in the MF compared to the HM group. At PND 51, 10 miRNAs were decreased and 17 were increased in the MF relative to HM suggesting the persistent effect of neonatal diet. miR-148 and miR-181 were decreased in MF compared to HM at PND 21. Let-7 was decreased at PND 35 while miR-199a and miR-199b were increased at PND 51 in MF compared to HM. Pathway analysis suggested that many of the miRNAs are involved in immune function. In conclusion, we observed differential expression of blood miRNAs at both PND 21 and PND 51. miRNA found in breastmilk were decreased in the serum of the MF group, suggesting that diet impacts circulating miRNA profiles at PND 21. The miRNAs continue to be altered at PND 51 suggesting a persistent effect of the neonatal diet. The sources of miRNAs in circulation need to be evaluated, as the piglets were fed the same solid diet leading up to PND 51 collections. In conclusion, the HM diet appears to have an immediate and persistent effect on the miRNA profile and likely regulates the pathways that impact the immune system and pose benefits to breastfed infants.

Small-Medium Extracellular Vesicles and Their miRNA Cargo in Retinal Health and Degeneration: Mediators of Homeostasis, and Vehicles for Targeted Gene Therapy

Photoreceptor cell death and inflammation are known to occur progressively in retinal degenerative diseases such as age-related macular degeneration (AMD). However, the molecular mechanisms underlying these biological processes are largely unknown. Extracellular vesicles (EV) are essential mediators of cell-to-cell communication with emerging roles in the modulation of immune responses. EVs, including exosomes, encapsulate and transfer microRNA (miRNA) to recipient cells and in this way can modulate the environment of recipient cells. Dysregulation of EVs however is correlated to a loss of cellular homeostasis and increased inflammation. In this work we investigated the role of isolated retinal small-medium sized EV (s-mEV) which includes exosomes in both the healthy and degenerating retina. Isolated s-mEV from normal retinas were characterized using dynamic light scattering, transmission electron microscopy and western blotting, and quantified across 5 days of photo-oxidative damage-induced degeneration using nanotracking analysis. Small RNAseq was used to characterize the miRNA cargo of retinal s-mEV isolated from healthy and damaged retinas. Finally, the effect of exosome inhibition on cell-to-cell miRNA transfer and immune modulation was conducted using systemic daily administration of exosome inhibitor GW4869 and in situ hybridization of s-mEV-abundant miRNA, miR-124-3p. Electroretinography and immunohistochemistry was performed to assess functional and morphological changes to the retina as a result of GW4869-induced exosome depletion. Results demonstrated an inverse correlation between s-mEV concentration and photoreceptor survivability, with a decrease in s-mEV numbers following degeneration. Small RNAseq revealed that s-mEVs contained uniquely enriched miRNAs in comparison to in whole retinal tissue, however, there was no differential change in the s-mEV miRNAnome following photo-oxidative damage. Exosome inhibition via the use of GW4869 was also found to exacerbate retinal degeneration, with reduced retinal function and increased levels of inflammation and cell death demonstrated following photo-oxidative damage in exosome-inhibited mice. Further, GW4869-treated mice displayed impaired translocation of photoreceptor-derived miR-124-3p to the inner retina during damage. Taken together, we propose that retinal s-mEV and their miRNA cargo play an essential role in maintaining retinal homeostasis through immune-modulation, and have the potential to be used in targeted gene therapy for retinal degenerative diseases.

mmannot: How to improve small-RNA annotation?

High-throughput sequencing makes it possible to provide the genome-wide distribution of small non coding RNAs in a single experiment, and contributed greatly to the identification and understanding of these RNAs in the last decade. Small non coding RNAs gather a wide collection of classes, such as microRNAs, tRNA-derived fragments, small nucleolar RNAs and small nuclear RNAs, to name a few. As usual in RNA-seq studies, the sequencing step is followed by a feature quantification step: when a genome is available, the reads are aligned to the genome, their genomic positions are compared to the already available annotations, and the corresponding features are quantified. However, problem arises when many reads map at several positions and while different strategies exist to circumvent this problem, all of them are biased. In this article, we present a new strategy that compares all the reads that map at several positions, and their annotations when available. In many cases, all the hits co-localize with the same feature annotation (a duplicated miRNA or a duplicated gene, for instance). When different annotations exist for a given read, we propose to merge existing features and provide the counts for the merged features. This new strategy has been implemented in a tool, mmannot, freely available at https://github.com/mzytnicki/mmannot.

Multi-Omic Regulation of the PAM50 Gene Signature in Breast Cancer Molecular Subtypes

Breast cancer is a disease that exhibits heterogeneity that goes from the genomic to the clinical levels. This heterogeneity is thought to be captured (at least partially) by the so-called breast cancer molecular subtypes. These molecular subtypes were initially defined based on the unsupervised clustering of gene expression and its correlate with histological, morphological, phenotypic and clinical features already known. Later, a 50-gene signature, PAM50, was defined in order to identify the biological subtype of a given sample within the clinical setting. The PAM50 signature was obtained by the use of unsupervised statistical methods, and therefore no limitation was set on the biological relevance (or lack of) of the selected genes beyond its predictive capacity. An open question that remains is what are the regulatory elements that drive the various expression behaviors of this set of genes in the different molecular subtypes. This question becomes more relevant as the measurement of more biological layers of regulation becomes accessible. In this work, we analyzed the gene expression regulation of the 50 genes in the PAM50 signature, in terms of (a) gene co-expression, (b) transcription factors, (c) micro-RNAs, and (d) methylation. Using data from the Cancer Genome Atlas (TCGA) for the Luminal A and B, Basal, and HER2-enriched molecular subtypes as well as normal tumor adjacent tissue, we identified predictors for gene expression through the use of an elastic net model. We compare and contrast the sets of identified regulators for the gene signature in each molecular subtype, and systematically compare them to current literature. We also identified a unique set of predictors for the expression of genes in the PAM50 signature associated with each of the molecular subtypes. Most selected predictors are exclusive for a PAM50 gene and predictors are not shared across subtypes. There are only 13 coding transcripts and 2 miRNAs selected for the four subtypes. MiR-21 and miR-10b connect almost all the PAM50 genes in all the subtypes and normal tissue, but do it in an exclusive manner, suggesting a cancer switch from miR-10b coordination in normal tissue to miR-21. The PAM50 gene sets of selected predictors that enrich for a function across subtypes, support that different regulatory molecular mechanisms are taking place. With this study we aim to a wider understanding of the regulatory mechanisms that differentiate the expression of the PAM50 signature, which in turn could perhaps help understand the molecular basis of the differences between the molecular subtypes.

Sensitive detection of EBV microRNAs across cancer spectrum reveals association with decreased survival in adult acute myelocytic leukemia

Epstein Barr virus (EBV) is the etiologic agent involved in numerous human cancers. After infecting the host, EBV establishes a latent infection, with low levels of messenger RNA (mRNA) and protein expression, evolved to evade immune recognition. Conversely, EBV microRNAs (miRNA) are expressed ubiquitously and abundantly within infected cells. Their role in tumor biology and clinical outcomes across the spectrum of cancer is not fully explained. Here, we applied our bioinformatics pipeline for quantitative EBV miRNA detection to examine sequencing data of 8,955 individual tumor samples across 27 tumor types representing the breadth of cancer. We uncover an association of intermediate levels of viral miRNA with decreased survival in adult acute myeloid leukemia (AML) patients (P = 0.00013). Prognostic modeling of this association suggests that increased EBV miRNA levels represent an independent risk factor for poor patient outcomes. Furthermore, we explore differences in expression between elevated and absent viral miRNA loads in adult AML tumors finding that EBV positivity was associated with proinflammatory signals. Together, given no associations were found for pediatric AML, our analyses suggests EBV positivity has the potential for being a prognostic biomarker and might represent a surrogate measure related to immune impairment in adult patients.