Integrated genomic analysis of mitochondrial RNA processing in human cancers

Translatomics to explore dynamic differences in immunocytes in the tumor microenvironment

Protein is an essential component of all living organisms and is primarily responsible for life activities; furthermore, its synthesis depends on a highly complex and accurate translation system. For proteins, the regulation at the translation level exceeds the sum of that during transcription, mRNA degradation, and protein degradation. Therefore, it is necessary to study regulation at the translation level. Imbalance in the translation process may change the cellular landscape, which not only leads to the occurrence, maintenance, progression, invasion, and metastasis of cancer but also affects the function of immune cells and changes the tumor microenvironment. Detailed analysis of transcriptional and protein atlases is needed to better understand how gene translation occurs. However, a more rigorous direct correlation between mRNA and protein levels is needed, which somewhat limits further studies. Translatomics is a technique for capturing and sequencing ribosome-related mRNAs that can effectively identify translation changes caused by ribosome stagnation and local translation abnormalities during cancer occurrence to further understand the changes in the translation landscape of cancer cells themselves and immune cells in the tumor microenvironment, which can provide new strategies and directions for tumor treatment.

High expression of novel biomarker TBRG4 promotes the progression and invasion of oral squamous cell carcinoma

BACKGROUND

Transforming growth factor β regulator 4 (TBRG4) is a potential prognostic indicator in various cancers, especially squamous cell carcinomas, and is associated with disease amelioration and poor outcomes. The study aimed to assess the expression pattern of TBRG4 in patients with operable oral squamous cell carcinoma (OSCC) to understand its role in tumour progression using indicators of disease outcome like tumour stage, grade, nodal metastasis, and pattern of invasion.

METHODS

TBRG4 expression was assessed by analyzing 51 cancer and adjacent non-cancerous tissues of OSCC patients using quantitative real-time PCR, and Western blot. TBRG4 expression was also analysed in The Cancer Genome Atlas Head-Neck Squamous Cell Carcinoma (TCGA-HNSC) dataset using the UALCAN tool (http://ualcan.path.uab.edu/). The relationship between TBRG4 expression and the patient's prognosis was analysed using Kaplan-Meier plotter.

RESULTS

Both mRNA and protein levels of TBRG4 were significantly increased in OSCC tissues. The TBGR4 expression was significantly associated with advanced stages (III and IV) and the worst pattern of invasion (WPOI-4 and 5). High TBRG4 expression was also significantly associated with reduced overall survival (p = 0.011). In addition, the analysis of TBRG4 gene expression and clinical data from TCGA, identified that TBRG4 was highly expressed in HPV negative OSCC patients and positively correlated with worst overall survival.

CONCLUSION

The present study suggests that the high expression of TBRG4 might serve as a novel prognostic biomarker for HPV-negative OSCC, which can be validated by future additional investigations in larger cohorts along with functional studies.

Dissecting the alternation landscape of mitochondrial metabolism-related genes in lung adenocarcinoma and their latent mechanisms

Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer with high incidence and unsatisfactory prognosis. The majority of LUAD patients eventually succumb to local and/or distinct metastatic recurrence. Genomic research of LUAD has broadened our understanding of this disease's biology and improved target therapies. However, the alternation landscape and characteristics of mitochondrial metabolism-related genes (MMRGs) in LUAD progression remain poorly understood. We performed a comprehensive analysis to identify the function and mechanism of MMRGs in LUAD based on the TCGA and GEO databases, which might offer therapeutic values for clinical researchers. Then, we figured out three hub prognosis-associated MMRGs (also termed as PMMRGs: ACOT11, ALDH2, and TXNRD1) that were engaged in the evolution of LUAD. To investigate the correlation between clinicopathological characteristics and MMRGs, we divided LUAD samples into two clusters (C1 and C2) based on key MMRGs. In addition, important pathways and the immune infiltration landscape affected by LUAD clusters were also delineated. Further, we nominated potential regulatory mechanisms underlying the MMRGs in LUAD development and progression. In conclusion, our integrative analysis enables a more comprehensive understanding of the mutation landscape of MMRGs in LUAD and provides an opportunity for more precise treatment.

Omics-based approaches for the systematic profiling of mitochondrial biology

Mitochondria are essential for cellular functions such as metabolism and apoptosis. They dynamically adapt to the changing environmental demands by adjusting their protein, nucleic acid, metabolite, and lipid contents. In addition, the mitochondrial components are modulated on different levels in response to changes, including abundance, activity, and interaction. A wide range of omics-based approaches has been developed to be able to explore mitochondrial adaptation and how mitochondrial function is compromised in disease contexts. Here, we provide an overview of the omics methods that allow us to systematically investigate the different aspects of mitochondrial biology. In addition, we show examples of how these methods have provided new biological insights. The emerging use of these toolboxes provides a more comprehensive understanding of the processes underlying mitochondrial function.

The emergent role of mitochondrial RNA modifications in metabolic alterations

Mitochondrial epitranscriptomics refers to the modifications occurring in all the different RNA types of mitochondria. Although the number of mitochondrial RNA modifications is less than those in cytoplasm, substantial evidence indicates that they play a critical role in accurate protein synthesis. Recent evidence supported those modifications in mitochondrial RNAs also have crucial implications in mitochondrial-related diseases. In the light of current knowledge about the involvement, the association between mitochondrial RNA modifications and diseases arises from studies focusing on mutations in both mitochondrial and nuclear DNA genes encoding enzymes involved in such modifications. Here, we review the current evidence available for mitochondrial RNA modifications and their role in metabolic disorders, and we also explore the possibility of using them as promising targets for prevention and early detection. Finally, we discuss future directions of mitochondrial epitranscriptomics in these metabolic alterations, and how these RNA modifications may offer a new diagnostic and theragnostic avenue for preventive purposes. This article is categorized under: RNA Processing > RNA Editing and Modification.

Targeted Mitochondrial Epigenetics: A New Direction in Alzheimer’s Disease Treatment

Mitochondrial epigenetic alterations are closely related to Alzheimer’s disease (AD), which is described in this review. Reports of the alteration of mitochondrial DNA (mtDNA) methylation in AD demonstrate that the disruption of the dynamic balance of mtDNA methylation and demethylation leads to damage to the mitochondrial electron transport chain and the obstruction of mitochondrial biogenesis, which is the most studied mitochondrial epigenetic change. Mitochondrial noncoding RNA modifications and the post-translational modification of mitochondrial nucleoproteins have been observed in neurodegenerative diseases and related diseases that increase the risk of AD. Although there are still relatively few mitochondrial noncoding RNA modifications and mitochondrial nuclear protein post-translational modifications reported in AD, we have reason to believe that these mitochondrial epigenetic modifications also play an important role in the AD process. This review provides a new research direction for the AD mechanism, starting from mitochondrial epigenetics. Further, this review summarizes therapeutic approaches to targeted mitochondrial epigenetics, which is the first systematic summary of therapeutic approaches in the field, including folic acid supplementation, mitochondrial-targeting antioxidants, and targeted ubiquitin-specific proteases, providing a reference for therapeutic targets for AD.

Identification of human mitochondrial RNA cleavage sites and candidate RNA processing factors

BACKGROUND

The human mitochondrial genome is transcribed as long strands of RNA containing multiple genes, which require post-transcriptional cleavage and processing to release functional gene products that play vital roles in cellular energy production. Despite knowledge implicating mitochondrial post-transcriptional processes in pathologies such as cancer, cardiovascular disease and diabetes, very little is known about the way their function varies on a human population level and what drives changes in these processes to ultimately influence disease risk. Here, we develop a method to detect and quantify mitochondrial RNA cleavage events from standard RNA sequencing data and apply this approach to human whole blood data from > 1000 samples across independent cohorts.

RESULTS

We detect 54 putative mitochondrial RNA cleavage sites that not only map to known gene boundaries, short RNA ends and RNA modification sites, but also occur at internal gene positions, suggesting novel mitochondrial RNA cleavage junctions. Inferred RNA cleavage rates correlate with mitochondrial-encoded gene expression across individuals, suggesting an impact on downstream processes. Furthermore, by comparing inferred cleavage rates to nuclear genetic variation and gene expression, we implicate multiple genes in modulating mitochondrial RNA cleavage (e.g. MRPP3, TBRG4 and FASTKD5), including a potentially novel role for RPS19 in influencing cleavage rates at a site near to the MTATP6-COX3 junction that we validate using shRNA knock down data.

CONCLUSIONS

We identify novel cleavage junctions associated with mitochondrial RNA processing, as well as genes newly implicated in these processes, and detect the potential impact of variation in cleavage rates on downstream phenotypes and disease processes. These results highlight the complexity of the mitochondrial transcriptome and point to novel mechanisms through which nuclear-encoded genes can potentially influence key mitochondrial processes.

Mitochondrial Genetic and Epigenetic Regulations in Cancer: Therapeutic Potential

Mitochondria are dynamic organelles managing crucial processes of cellular metabolism and bioenergetics. Enabling rapid cellular adaptation to altered endogenous and exogenous environments, mitochondria play an important role in many pathophysiological states, including cancer. Being under the control of mitochondrial and nuclear DNA (mtDNA and nDNA), mitochondria adjust their activity and biogenesis to cell demands. In cancer, numerous mutations in mtDNA have been detected, which do not inactivate mitochondrial functions but rather alter energy metabolism to support cancer cell growth. Increasing evidence suggests that mtDNA mutations, mtDNA epigenetics and miRNA regulations dynamically modify signalling pathways in an altered microenvironment, resulting in cancer initiation and progression and aberrant therapy response. In this review, we discuss mitochondria as organelles importantly involved in tumorigenesis and anti-cancer therapy response. Tumour treatment unresponsiveness still represents a serious drawback in current drug therapies. Therefore, studying aspects related to genetic and epigenetic control of mitochondria can open a new field for understanding cancer therapy response. The urgency of finding new therapeutic regimens with better treatment outcomes underlines the targeting of mitochondria as a suitable candidate with new therapeutic potential. Understanding the role of mitochondria and their regulation in cancer development, progression and treatment is essential for the development of new safe and effective mitochondria-based therapeutic regimens.

Dysregulated mitochondrial and cytosolic tRNA m1A methylation in Alzheimer's disease

RNA modifications affect many aspects of RNA metabolism and are involved in the regulation of many different biological processes. Mono-methylation of adenosine in the N1 position, N1-methyladensoine (m1A), is a reversible modification that is known to target rRNAs and tRNAs. m1A has been shown to increase tRNA structural stability and induce correct tRNA folding. Recent studies have begun to associate the dysregulation of epitranscriptomic control with age-related disorders such as Alzheimer's disease. Here, we applied the newly developed m1A-quant-seq approach to map the brain abundant m1A RNA modification in the cortex of an Alzheimer's disease mouse model, 5XFAD. We observed hypomethylation in both mitochondrial and cytosolic tRNAs in 5XFAD mice compared with wild type. Furthermore, the main enzymes responsible for the addition of m1A in mitochondrial (TRMT10C, HSD17B10) and cytosolic tRNAs (TRMT61A) displayed decreased expression in 5XFAD compared with wild-type mice. Knockdown of these enzymes results in a more severe phenotype in a Drosophila tau model, and differential m1A methylation is correlated with differences in mature mitochondrial tRNA expression. Collectively, this work suggests that hypo m1A modification in tRNAs may play a role in Alzheimer's disease pathogenesis.

A Novel Mitochondrial-Related Nuclear Gene Signature Predicts Overall Survival of Lung Adenocarcinoma Patients

Background: Lung cancer is the leading cause of cancer-related death worldwide, of which lung adenocarcinoma (LUAD) is one of the main histological subtypes. Mitochondria are vital for maintaining the physiological function, and their dysfunction has been found to be correlated with tumorigenesis and disease progression. Although, some mitochondrial-related genes have been found to correlate with the clinical outcomes of multiple tumors solely. The integrated relationship between nuclear mitochondrial genes (NMGs) and the prognosis of LUAD remains unclear.

Methods: The list of NMGs, gene expression data, and related clinical information of LUAD were downloaded from public databases. Bioinformatics methods were used and obtained 18 prognostic related NMGs to construct a risk signature.

Results: There were 18 NMGs (NDUFS2, ATP8A2, SCO1, COX14, COA6, RRM2B, TFAM, DARS2, GARS, YARS2, EFG1, GFM1, MRPL3, MRPL44, ISCU, CABC1, HSPD1, and ETHE1) identified by LASSO regression analysis. The mRNA expression of these 18 genes was positively correlated with their relative linear copy number alteration (CNA). Meanwhile, the established risk signature could effectively distinguish high- and low-risk patients, and its predictive capacity was validated in three independent gene expression omnibus (GEO) cohorts. Notably, a significantly lower prevalence of actionable EGFR alterations was presented in patients with high-risk NMGs signature but accompanied with a more inflame immune tumor microenvironment. Additionally, multicomponent Cox regression analysis showed that the model was stable when risk score, tumor stage, and lymph node stage were considered, and the 1-, 3-, and 5-year AUC were 0.74, 0.75, and 0.70, respectively.

Conclusion: Together, this study established a signature based on NMGs that is a prognostic biomarker for LUAD patients and has the potential to be widely applied in future clinical settings.

The role of RNA m5C modification in cancer metastasis

Epigenetic modification plays a crucial regulatory role in the biological processes of eukaryotic cells. The recent characterization of DNA and RNA methylation is still ongoing. Tumor metastasis has long been an unconquerable feature in the fight against cancer. As an inevitable component of the epigenetic regulatory network, 5-methylcytosine is associated with multifarious cellular processes and systemic diseases, including cell migration and cancer metastasis. Recently, gratifying progress has been achieved in determining the molecular interactions between m⁵C writers (DNMTs and NSUNs), demethylases (TETs), readers (YTHDF2, ALYREF and YBX1) and RNAs. However, the underlying mechanism of RNA m⁵C methylation in cell mobility and metastasis remains unclear. The functions of m⁵C writers and readers are believed to regulate gene expression at the post-transcription level and are involved in cellular metabolism and movement. In this review, we emphatically summarize the recent updates on m⁵C components and related regulatory networks. The content will be focused on writers and readers of the RNA m⁵C modification and potential mechanisms in diseases. We will discuss relevant upstream and downstream interacting molecules and their associations with cell migration and metastasis.

Quantifying RNA modifications by mass spectrometry: a novel source of biomarkers in oncology

Despite significant progress in targeted therapies, cancer recurrence remains a major cause of mortality worldwide. Identification of accurate biomarkers, through molecular profiling in healthy and cancer patient samples, will improve diagnosis and promote personalized medicine. While genetic and epigenetic alterations of DNA are currently exploited as cancer biomarkers, their robustness is limited by tumor heterogeneity. Recently, cancer-associated changes in RNA marks have emerged as a promising source of diagnostic and prognostic biomarkers. RNA epigenetics (also known as epitranscriptomics) is an emerging field in which at least 150 chemical modifications in all types of RNA (mRNA, tRNA, lncRNA, rRNA, and microRNA) have been detected. These modifications fine-tune gene expression in both physiological and pathological processes. A growing number of studies have established links between specific modified nucleoside levels in solid/liquid biopsies, and cancer onset and progression. In this review, we highlight the potential role of epitranscriptomic markers in refining cancer diagnosis and/or prognosis. RNA modification patterns may contain important information for establishing an initial diagnosis, monitoring disease evolution, and predicting response to treatment. Furthermore, recent developments in mass spectrometry allow reliable quantification of RNA marks in solid biopsies and biological fluids. We discuss the great potential of mass spectrometry for identifying epitranscriptomic biomarker signatures in cancer diagnosis. While there are various methods to quantify modified nucleosides, most are unable to detect and quantify more than one type of RNA modification at a time. Mass spectrometry analyses, especially GC-MS/MS and LC-MS/MS, overcome this limitation and simultaneously detect modified nucleosides by multiple reaction monitoring. Indeed, several groups are currently validating mass spectrometry methods that quantify several nucleosides at one time in liquid biopsies. The challenge now is to exploit these powerful analytical tools to establish epitranscriptomic signatures that should open new perspectives in personalized medicine. This review summarizes the growing clinical field of analysis of RNA modifications and discusses pre-analytical and analytical approaches, focusing in particular on the development of new mass spectrometry tools and their clinical applications.

Potential biomarkers and targets of mitochondrial dynamics

Mitochondrial dysfunction contributes to the imbalance of cellular homeostasis and the development of diseases, which is regulated by mitochondria-associated factors. The present review aims to explore the process of the mitochondrial quality control system as a new source of the potential diagnostic biomarkers and/or therapeutic targets for diseases, including mitophagy, mitochondrial dynamics, interactions between mitochondria and other organelles (lipid droplets, endoplasmic reticulum, endosomes, and lysosomes), as well as the regulation and posttranscriptional modifications of mitochondrial DNA/RNA (mtDNA/mtRNA). The direct and indirect influencing factors were especially illustrated in understanding the interactions among regulators of mitochondrial dynamics. In addition, mtDNA/mtRNAs and proteomic profiles of mitochondria in various lung diseases were also discussed as an example. Thus, alternations of mitochondria-associated regulators can be a new category of biomarkers and targets for disease diagnosis and therapy.

The Uprising of Mitochondrial DNA Biomarker in Cancer

Cancer is a heterogeneous group of diseases, the progression of which demands an accumulation of genetic mutations and epigenetic alterations of the human nuclear genome or possibly in the mitochondrial genome as well. Despite modern diagnostic and therapeutic approaches to battle cancer, there are still serious concerns about the increase in death from cancer globally. Recently, a growing number of researchers have extensively focused on the burgeoning area of biomarkers development research, especially in noninvasive early cancer detection. Intergenomic cross talk has triggered researchers to expand their studies from nuclear genome-based cancer researches, shifting into the mitochondria-mediated associations with carcinogenesis. Thus, it leads to the discoveries of established and potential mitochondrial biomarkers with high specificity and sensitivity. The research field of mitochondrial DNA (mtDNA) biomarkers has the great potential to confer vast benefits for cancer therapeutics and patients in the future. This review seeks to summarize the comprehensive insights of nuclear genome cancer biomarkers and their usage in clinical practices, the intergenomic cross talk researches that linked mitochondrial dysfunction to carcinogenesis, and the current progress of mitochondrial cancer biomarker studies and development.

Synaptic processes and immune-related pathways implicated in Tourette syndrome

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

Mitochondrial tRNA methylation in Alzheimer's disease and progressive supranuclear palsy

BACKGROUND

Methylation of mitochondrial tRNAs (mt-tRNA) at the 9th position ("p9 site") is known to impact translational efficiency and downstream mitochondrial function; however, direct assessment of mt-RNA methylation is challenging. Recent RNA sequence-based methods have been developed to reliably identify post-transcriptional methylation. Though p9 methylation has been studied in healthy human populations and in the context of cancer, it has not yet been analyzed in neurodegenerative disease, where mitochondrial dysfunction is a prominent and early hallmark of disease progression.

METHODS

Mitochondrial p9 methylation was inferred from multi-allelic calls in RNA-seq data. Gene-based association studies were performed in FUMA. Correlations between nuclear gene expression and p9 methylation were tested using Spearman's rho. Fisher's Exact test was used in PANTHER and IPA to test for overrepresentation and enrichment of biological processes and pathways in the top nuclear genes correlated with p9 methylation.

RESULTS

Variable methylation was observed at 11 p9 sites in post-mortem cerebellar tissue of elderly subjects who were either healthy or diagnosed with Alzheimer's disease (AD), progressive supranuclear palsy (PSP) or pathological aging (PA). Similarities in degree of methylation were observed between AD and PSP. Certain nuclear encoded genes were identified as significantly associated with p9 methylation. Expression of 5300 nuclear encoded genes was significantly correlated with p9 methylation, with AD and PSP subjects exhibiting similar expression profiles. Overrepresentation and enrichment testing using the top transcripts revealed enrichment for a number of molecular processes, terms and pathways including many of which that were mitochondrial-related.

CONCLUSION

With mitochondrial dysfunction being an established hallmark of neurodegenerative disease pathophysiology, this work sheds light on the potential molecular underpinnings of this dysfunction. Here we show overlap in cerebellar pathophysiology between common tauopathies such as Alzheimer's disease and progressive supranuclear palsy. Whether p9 hypermethylation is a cause or consequence of pathology remains an area of focus.

Analysis of mitochondrial m1A/G RNA modification reveals links to nuclear genetic variants and associated disease processes

RNA modifications affect the stability and function of RNA species, regulating important downstream processes. Modification levels are often dynamic, varying between tissues and individuals, although it is not always clear what modulates this or what impact it has on biological systems. Here, we quantify variation in m1A/G RNA modification levels at functionally important positions in the human mitochondrial genome across 11,552 samples from 39 tissue/cell types and find that modification levels are associated with mitochondrial transcript processing. We identify links between mitochondrial RNA modification levels and genetic variants in the nuclear genome, including a missense mutation in LONP1, and find that genetic variants within MRPP3 and TRMT61B are associated with RNA modification levels across a large number of tissues. Genetic variants linked to RNA modification levels are associated with multiple disease/disease-related phenotypes, including blood pressure, breast cancer and psoriasis, suggesting a role for mitochondrial RNA modification in complex disease.

Ali et al. analyze publicly available RNA-seq data from different tissues to quantify variation in m1A/G methylation levels in mitochondrial RNAs. They show a link between mitochondrial m1A/G modification levels and nuclear genetic variants, many of which are associated with disease.

Mitoepigenetics and Its Emerging Roles in Cancer

In human beings, there is a ∼16,569 bp circular mitochondrial DNA (mtDNA) encoding 22 tRNAs, 12S and 16S rRNAs, 13 polypeptides that constitute the central core of ETC/OxPhos complexes, and some non-coding RNAs. Recently, mtDNA has been shown to have some covalent modifications such as methylation or hydroxylmethylation, which play pivotal epigenetic roles in mtDNA replication and transcription. Post-translational modifications of proteins in mitochondrial nucleoids such as mitochondrial transcription factor A (TFAM) also emerge as essential epigenetic modulations in mtDNA replication and transcription. Post-transcriptional modifications of mitochondrial RNAs (mtRNAs) including mt-rRNAs, mt-tRNAs and mt-mRNAs are important epigenetic modulations. Besides, mtDNA or nuclear DNA (n-DNA)-derived non-coding RNAs also play important roles in the regulation of translation and function of mitochondrial genes. These evidences introduce a novel concept of mitoepigenetics that refers to the study of modulations in the mitochondria that alter heritable phenotype in mitochondria itself without changing the mtDNA sequence. Since mitochondrial dysfunction contributes to carcinogenesis and tumor development, mitoepigenetics is also essential for cancer. Understanding the mode of actions of mitoepigenetics in cancers may shade light on the clinical diagnosis and prevention of these diseases. In this review, we summarize the present study about modifications in mtDNA, mtRNA and nucleoids and modulations of mtDNA/nDNA-derived non-coding RNAs that affect mtDNA translation/function, and overview recent studies of mitoepigenetic alterations in cancer.

Nuclear genetic regulation of the human mitochondrial transcriptome

Mitochondria play important roles in cellular processes and disease, yet little is known about how the transcriptional regime of the mitochondrial genome varies across individuals and tissues. By analyzing >11,000 RNA-sequencing libraries across 36 tissue/cell types, we find considerable variation in mitochondrial-encoded gene expression along the mitochondrial transcriptome, across tissues and between individuals, highlighting the importance of cell-type specific and post-transcriptional processes in shaping mitochondrial-encoded RNA levels. Using whole-genome genetic data we identify 64 nuclear loci associated with expression levels of 14 genes encoded in the mitochondrial genome, including missense variants within genes involved in mitochondrial function (TBRG4, MTPAP and LONP1), implicating genetic mechanisms that act in trans across the two genomes. We replicate ~21% of associations with independent tissue-matched datasets and find genetic variants linked to these nuclear loci that are associated with cardio-metabolic phenotypes and Vitiligo, supporting a potential role for variable mitochondrial-encoded gene expression in complex disease.

Mitochondrial RNA Expression and Single Nucleotide Variants in Association with Clinical Parameters in Primary Breast Cancers

The human mitochondrial DNA (mtDNA) encodes 37 genes, including thirteen proteins essential for the respiratory chain, and RNAs functioning in the mitochondrial translation apparatus. The total number of mtDNA molecules per cell (mtDNA content) is variable between tissue types and also between tumors and their normal counterparts. For breast cancer, tumors tend to be depleted in their mtDNA content compared to adjacent normal mammary tissue. Various studies have shown that primary breast tumors harbor somatic mtDNA variants. A decrease in mtDNA content or the presence of somatic variants could indicate a reduced mitochondrial function within breast cancer. In this explorative study we aimed to further understand genomic changes and expression of the mitochondrial genome within breast cancer, by analyzing RNA sequencing data of primary breast tumor specimens of 344 cases. We demonstrate that somatic variants detected at the mtRNA level are representative for somatic variants in the mtDNA. Also, the number of somatic variants within the mitochondrial transcriptome is not associated with mutational processes impacting the nuclear genome, but is positively associated with age at diagnosis. Finally, we observe that mitochondrial expression is related to ER status. We conclude that there is a large heterogeneity in somatic mutations of the mitochondrial genome within primary breast tumors, and differences in mitochondrial expression among breast cancer subtypes. The exact impact on metabolic differences and clinical relevance deserves further study.

Dissecting the expression landscape of mitochondrial genes in lung squamous cell carcinoma and lung adenocarcinoma

Lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD) are the two major subtypes of lung cancer. To explore mitochondrial respiratory gene expression profiles in LUSC and LUAD, RNA sequencing data from The Cancer Genome Atlas was used for comprehensive analyses to establish the molecular characteristics of LUSC and LUAD. To elucidate expression profiles, subtypes were defined using unsupervised clustering of mitochondrial gene expression data. Differences in nuclear gene expression levels, signaling pathways and tumor microenvironments between subtypes were investigated. The analysis revealed that mitochondrial respiratory genes were generally expressed at lower levels in tumor tissues compared with matched control tissues. The expression of mitochondrially encoded NADH dehydrogenase 5 or 6 was associated with tumor progression in LUAD and LUSC. Patients were clustered into three subgroups based on the expression profile of 13 mitochondrial protein-encoding genes, and patients in Cluster 3 exhibited poor survival rates compared with patients from Cluster 1. Furthermore, this association was also observed in another independent data set. Further analyses of the expression of nuclear-encoded genes in the three clusters revealed the enrichment of several cancer-associated signaling pathways in Cluster 3, particularly the apoptotic signaling pathway, suggesting a potential association between the decreased expression of mitochondrial DNA genes and increased tumor aggressiveness. Furthermore, the analyses of immune cell compositions in the tumor microenvironment detected a significant increase in the proportion of CD4+ T cells and a decrease in the proportion of macrophages in LUAD compared with LUSC (P=0.0000104 and P=0.0000105, respectively). In conclusion, the present study revealed an association between the expression patterns of mitochondrial-encoded genes and lung cancer, which may contribute to novel therapeutic strategies for patients with LUSC and LUAD.