MicroRNA‑mRNA integrated analysis based on a case of well‑differentiated thyroid cancer with both metastasis and metastatic recurrence

MiRNA expression deregulation correlates with the Oncotype DX® DCIS score

Background

Current clinical criteria do not discriminate well between women who will or those who will not develop ipsilateral invasive breast cancer (IBC), or a DCIS recurrence after a ductal carcinoma in situ (DCIS) diagnosis. The 12-gene Oncotype DX® DCIS assay (RT qPCR gene-based scoring system) was established and shown to predict the risk of subsequent ipsilateral IBC or DCIS recurrence. Recent studies have shown that microRNA (miRNA) expression deregulation can contribute to the development of IBC, but very few have evaluated miRNA deregulation in DCIS lesions. In this study, we sought to determine whether specific miRNA expression changes may correlate with Oncotype DX® DCIS scores.

Methods

For this study, we used archived formalin-fixed, paraffin-embedded (FFPE) specimens from 41 women diagnosed with DCIS between 2012 and 2018. The DCIS lesions were stratified into low (n = 26), intermediate (n = 10), and high (n = 5) risk score groups using the Oncotype DX® DCIS assay. Total RNA was extracted from DCIS lesions by macro-dissection of unstained FFPE sections, and next-generation small-RNA sequencing was performed. We evaluated the correlation between miRNA expression data and Oncotype score, as well as patient age. RT-qPCR validations were performed to validate the topmost differentially expressed miRNAs identified between the different risk score groups.

Results

MiRNA sequencing of 32 FFPE DCIS specimens from the three different risk group scores identified a correlation between expression deregulation of 17 miRNAs and Oncotype scores. Our analyses also revealed a correlation between the expression deregulation of 9 miRNAs and the patient’s age. Based on these results, a total of 15 miRNAs were selected for RT-qPCR validation. Of these, miR-190b (p = 0.043), miR-135a (p = 0.05), miR-205 (p = 0.00056), miR-30c (p = 0.011), and miR-744 (p = 0.038) showed a decreased expression in the intermediate/high Oncotype group when compared to the low-risk score group. A composite risk score was established using these 5 miRNAs and indicated a significant association between miRNA expression deregulation and the Oncotype DX® DCIS Score (p < 0.0021), between high/intermediate and low risk groups.

Conclusions

Our analyses identified a subset of 5 miRNAs able to discriminate between Oncotype DX® DCIS score subgroups. Together, our data suggest that miRNA expression analysis may add value to the predictive and prognostic evaluation of DCIS lesions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-022-01558-4.

Profiling and Integrated Analysis of Differentially Expressed MicroRNAs as Novel Biomarkers of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a heterogeneous disease that has multiple etiologies. It is the most common primary liver cancer, the sixth highest cause of cancer incidences, and the fourth highest cause of cancer-related deaths. The discovery of new biomarkers for the early detection, treatment, and prognosis of HCC would therefore be extremely useful. This study investigated differentially expressed ribonucleic acid (RNA) profiles by constructing a genome-wide profile of clinical samples. Differential expression analysis identified 1,280 differentially expressed messenger RNAs (dif-mRNAs), 99 differentially expressed microRNAs (dif-miRNAs), 181 differentially expressed long non-coding RNAs (dif-lncRNAs), and 31 differentially expressed circular RNAs (dif-circRNAs). Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) path analysis were then conducted on these differentially expressed RNAs, revealing that they were clearly related to cell division, foreign body metabolism, and ribosome assembly. A competing endogenous RNA (ceRNA) network was then constructed based on the regulatory dif-miRNA-dif-mRNA and dif-miRNA-dif-lncRNA relationships. These results were also verified using HCC data from the Cancer Genome Atlas (TCGA); seven dif-miRNAs were verified in clinical samples by real-time quantitative polymerase chain reaction (RT-qPCR). Kaplan-Meier survival analysis revealed that the expression levels of Hsa-miR-1269a, Hsa-miR-421, and Hsa-miR-190b were correlated with overall survival. (P &lt;0.05). Survival analysis of clinical samples showed that hsa-mir-1269a, hsa-mir-421 were associated with prognosis (p&lt;0.05).This study revealed the general expression characteristics of specific differentially expressed miRNAs using a ceRNA network constructed from HCC samples. Hsa-mir-1269a, hsa-mir-421 may be promising candidates.

Potential biomarkers of abnormal osseointegration of implants in type II diabetes mellitus

Background

Type II diabetes mellitus (T2DM) is an important risk factor for osseointegration of implants. The aim of this study was to explore key genes of T2DM affecting bone metabolism through bioinformatic analysis of published RNA sequencing data, identify potential biomarkers, and provide a reference for finding the molecular mechanism of abnormal osseointegration caused by T2DM.

Methods

We identified differentially expressed mRNAs and miRNAs from the Gene Expression Omnibus database using the R package ‘limma’ and analysed the predicted target genes using Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis and Gene Ontology analysis. At the same time, miRNA–mRNA interactions were explored using miRWalk 2.0.

Results

We constructed an miRNA-gene regulatory network and a protein–protein interaction network. The enrichment pathways of differentially expressed mRNAs included extracellular matrix receptor interactions, protein digestion and absorption, the PI3K-Akt signalling pathway, cytokine–cytokine receptor interactions, chemokine signalling pathways, and haematopoietic cell lineage functions. We analysed the expression of these differentially expressed mRNAs and miRNAs in T2DM rats and normal rats with bone implants and identified Smpd3, Itga10, and rno-mir-207 as possible key players in osseointegration in T2DM.

Conclusion

Smpd3, Itga10, and rno-mir-207 are possible biomarkers of osseointegration in T2DM. This study sheds light on the possible molecular mechanism of abnormal osseointegration caused by bone metabolism disorder in T2DM.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12903-021-01939-9.

Whole-Transcriptome Sequencing Identifies Key Differentially Expressed mRNAs, miRNAs, lncRNAs, and circRNAs Associated with CHOL

To systematically evaluate the whole-transcriptome sequencing data of cholangiocarcinoma (CHOL) to gain more insights into the transcriptomic landscape and molecular mechanism of this cancer, we performed whole-transcriptome sequencing based on the tumorous (C) and their corresponding non-tumorous adjacent to the tumors (CP) from eight CHOL patients. Subsequently, differential expression analysis was performed on the C and CP groups, followed by functional interaction prediction analysis to investigate gene-regulatory circuits in CHOL. In addition, The Cancer Genome Atlas (TCGA) for CHOL data was used to validate the results. In total, 2,895 differentially expressed messenger RNAs (dif-mRNAs), 56 differentially expressed microRNAs (dif-miRNAs), 151 differentially expressed long non-coding RNAs (dif-lncRNAs), and 110 differentially expressed circular RNAs (dif-circRNAs) were found in CHOL samples compared with controls. Enrichment analysis on those differentially expressed genes (DEGs) related to miRNA, lncRNA, and circRNA also identified the function of spliceosome. The downregulated hsa-miR-144-3p were significantly enriched in the competing endogenous RNA (ceRNA) complex network, which also included 7 upregulated and 13 downregulated circRNAs, 7 upregulated lncRNAs, and 90 upregulated and 40 downregulated mRNAs. Moreover, most of the DEGs and a few of the miRNAs (such as hsa-miR-144-3p) were successfully validated by TCGA data. The genes involved in RNA splicing and protein degradation processes and miR-144-3p may play fundamental roles in the pathogenesis of CHOL.

Dual Regulatory Mechanisms of Expression and Mutation Involving Metabolism-Related Genes FDFT1 and UQCR5 during CLM

Colorectal cancer (CRC) is the third most common cancer worldwide, and liver metastasis presents a major cause of CRC-associated death. Extensive genomic analysis has provided valuable insight into the pathogenesis and progression of CRC; however, a comprehensive proteogenomic characterization of CRC liver metastasis (CLM) has yet to be reported. Here, we analyzed the proteomes of 44 paired normal colorectal tissues and CRC tissues with or without liver metastasis, as well as analyzed genomics of CRC characterized previously by The Cancer Genome Atlas (TCGA) to conduct integrated proteogenomic analyses. We identified a total of 2,170 significantly deregulated proteins associated with CLM, 14.88% of which were involved in metabolic pathways. The mutated peptide number was found to have potential prognosis value, and somatic variants revealed two metabolism-related genes UQCR5 and FDFT1 that frequently mutated only in the liver metastatic cohort and displayed dysregulated protein abundance with biological function and clinical significance in CLM. Proteogenomic characterization and integrative and comparative genomic analysis provides functional context and prognostic value to annotate genomic abnormalities and affords a new paradigm for understanding human colon and rectal cancer liver metastasis.