Corrections for mRNA extraction and sample normalization errors find increased mRNA levels may compensate for cancer haplo-insufficiency

High Matrix Metalloproteinase 28 Expression is Associated with Poor Prognosis in Pancreatic Adenocarcinoma

Purpose

Pancreatic adenocarcinoma (PAAD) is a devastating disease with high mortality and morbidity. Matrix metalloproteinase 28 (MMP28) has been associated with carcinogenesis of many human cancers. However, little is known about the potential prognostic value and underlying regulatory mechanisms of MMP28 in PAAD.

Methods

The relationship between MMP28 expression level and various clinicopathological parameters was analyzed in TCGA-PAAD cohorts. MMP28-correlated genes in the TCGA-PAAD cohort were identified and enrichment analysis according to the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes was conducted using LinkedOmics. Protein–protein interaction and transcription factors-miRNA co-regulatory networks were constructed with the use of NetworkAnalyst. Then, the distribution of immune cells related to MMP28 expression in blood was analyzed using the Human Protein Atlas, and the tumor microenvironment of PAAD was analyzed by the TIMER 2.0 database. To investigate the biological function of MMP28 in PAAD, siRNA was constructed to knock down the MMP28 gene in vitro.

Results

High MMP28 expression is associated with poor overall survival and disease-free survival in PAAD patients. The expression of MMP28 in PAAD is most significantly correlated with KRT19, IL1RN, and ANXA2 genes. Network analysis revealed that MIR-181 family, TAFs, and CDC6 are potential regulators of MMP28. Furthermore, naive CD4⁺ T cell, naive CD8⁺ T cell, and mucosal-associated invariant T cell enrichment in blood were correlated with MMP28 expression. Furthermore, high MMP28 expression was correlated with a decrease in B cell, naive CD4⁺ T cell, naive CD8⁺ T cell, and endothelial cell presence in the tumor microenvironment in PAAD. Finally, genetic knockdown of MMP28 could restrain the proliferation, migration, and invasion of PAAD cells.

Conclusion

Our findings indicate that high MMP28 expression in PAAD is associated with cancer progression, invasion, and metastasis. Hence, MMP28 might serve as an independent prognostic biomarker and a prospective therapeutic target for PAAD.

Identification and characterization of genes with absolute mRNA abundances changes in tumor cells with varied transcriptome sizes

BACKGROUND

The amount of RNA per cell, namely the transcriptome size, may vary under many biological conditions including tumor. If the transcriptome size of two cells is different, direct comparison of the expression measurements on the same amount of total RNA for two samples can only identify genes with changes in the relative mRNA abundances, i.e., cellular mRNA concentration, rather than genes with changes in the absolute mRNA abundances.

RESULTS

Our recently proposed RankCompV2 algorithm identify differentially expressed genes (DEGs) through comparing the relative expression orderings (REOs) of disease samples with that of normal samples. We reasoned that both the mRNA concentration and the absolute abundances of these DEGs must have changes in disease samples. In simulation experiments, this method showed excellent performance for identifying DEGs between normal and disease samples with different transcriptome sizes. Through analyzing data for ten cancer types, we found that a significantly higher proportion of the DEGs with absolute mRNA abundance changes overlapped or directly interacted with known cancer driver genes and anti-cancer drug targets than that of the DEGs only with mRNA concentration changes alone identified by the traditional methods. The DEGs with increased absolute mRNA abundances were enriched in DNA damage-related pathways, while DEGs with decreased absolute mRNA abundances were enriched in immune and metabolism associated pathways.

CONCLUSIONS

Both the mRNA concentration and the absolute abundances of the DEGs identified through REOs comparison change in disease samples in comparison with normal samples. In cancers these genes might play more important upstream roles in carcinogenesis.