OmicsVolcano: software for intuitive visualization and interactive exploration of high-throughput biological data

Cell cycle traverse rate predicts long-term outcomes in a multi-institutional cohort of patients with triple-negative breast cancer

BACKGROUND

Ki67 index (KI) and mitotic index (MI) are proliferation markers with established prognostic value in breast carcinomas. While KI is evaluated immunohistochemically and reported as a percentage, MI is determined visually and reflects total mitotic cells in 10 high-power fields. Our objective was to integrate KI and MI into a novel metric; the cell cycle traverse rate (CCTR). Given the lack of prognostic and predictive biomarkers in TNBC, we sought to assess the potential of CCTR as a risk-stratification tool for chemotherapy-treated TNBC patients from two independent cohorts: the Nottingham group (n = 124) and the Norway group (n = 71).

METHODS

We evaluated the ability of CCTR to predict survival after adjuvant chemotherapy for TNBC patients (n = 195) in two independent cohorts. Using immunohistochemistry and RNA sequencing, we determined the differences in immunohistochemical biomarkers, gene ontologies, molecular pathways and immune cell fractions based on CCTR.

RESULTS

TNBC shows a significantly lower median CCTR compared to luminal A (p < 0.01), luminal B (p < 0.01), and HER2+ samples (p < 0.01). CCTR outperformed both KI and MI in effectively risk-stratifying TNBC patients suggesting that combining KI and MI into a single metric, namely CCTR, could serve as a superior prognostic marker for Breast Cancer Specific Survival (BCSS) (p = 0.041). CCTR-high group exhibited enriched expression of various oncogenic signatures, including angiogenesis, epithelial-to-mesenchymal transition (EMT), Hedgehog signaling, hypoxia, Notch signaling, PI3K-AKT-mTOR signaling, TGFβ signaling, p53 signaling, and TNFα signaling via NFκB. These findings suggest the potential involvement of these pathways in the aggressiveness and clinical outcomes of TNBC patients.

CONCLUSIONS

Collectively, these findings suggest that CCTR offers superior predictive information compared to KI and MI alone with respect to long-term outcomes from adjuvant chemotherapy in patients with TNBC that may guide treatment decision making.

Quantitative subcellular reconstruction reveals a lipid mediated inter-organelle biogenesis network

The structures and functions of organelles in cells depend on each other but have not been systematically explored. We established stable knockout cell lines of peroxisomal, Golgi and endoplasmic reticulum genes identified in a whole-genome CRISPR knockout screen for inducers of mitochondrial biogenesis stress, showing that defects in peroxisome, Golgi and endoplasmic reticulum metabolism disrupt mitochondrial structure and function. Our quantitative total-organelle profiling approach for focussed ion beam scanning electron microscopy revealed in unprecedented detail that specific organelle dysfunctions precipitate multi-organelle biogenesis defects, impair mitochondrial morphology and reduce respiration. Multi-omics profiling showed a unified proteome response and global shifts in lipid and glycoprotein homeostasis that are elicited when organelle biogenesis is compromised, and that the resulting mitochondrial dysfunction can be rescued with precursors for ether-glycerophospholipid metabolic pathways. This work defines metabolic and morphological interactions between organelles and how their perturbation can cause disease.

Temporal landscape of mitochondrial proteostasis governed by the UPRmt

Breakdown of mitochondrial proteostasis activates quality control pathways including the mitochondrial unfolded protein response (UPRmt) and PINK1/Parkin mitophagy. However, beyond the up-regulation of chaperones and proteases, we have a limited understanding of how the UPRmt remodels and restores damaged mitochondrial proteomes. Here, we have developed a functional proteomics framework, termed MitoPQ (Mitochondrial Proteostasis Quantification), to dissect the UPRmt's role in maintaining proteostasis during stress. We find essential roles for the UPRmt in both protecting and repairing proteostasis, with oxidative phosphorylation metabolism being a central target of the UPRmt. Transcriptome analyses together with MitoPQ reveal that UPRmt transcription factors drive independent signaling arms that act in concert to maintain proteostasis. Unidirectional interplay between the UPRmt and PINK1/Parkin mitophagy was found to promote oxidative phosphorylation recovery when the UPRmt failed. Collectively, this study defines the network of proteostasis mediated by the UPRmt and highlights the value of functional proteomics in decoding stressed proteomes.

Uniformly shaped harmonization combines human transcriptomic data from different platforms while retaining their biological properties and differential gene expression patterns

Introduction: Co-normalization of RNA profiles obtained using different experimental platforms and protocols opens avenue for comprehensive comparison of relevant features like differentially expressed genes associated with disease. Currently, most of bioinformatic tools enable normalization in a flexible format that depends on the individual datasets under analysis. Thus, the output data of such normalizations will be poorly compatible with each other. Recently we proposed a new approach to gene expression data normalization termed Shambhala which returns harmonized data in a uniform shape, where every expression profile is transformed into a pre-defined universal format. We previously showed that following shambhalization of human RNA profiles, overall tissue-specific clustering features are strongly retained while platform-specific clustering is dramatically reduced. Methods: Here, we tested Shambhala performance in retention of fold-change gene expression features and other functional characteristics of gene clusters such as pathway activation levels and predicted cancer drug activity scores. Results: Using 6,793 cancer and 11,135 normal tissue gene expression profiles from the literature and experimental datasets, we applied twelve performance criteria for different versions of Shambhala and other methods of transcriptomic harmonization with flexible output data format. Such criteria dealt with the biological type classifiers, hierarchical clustering, correlation/regression properties, stability of drug efficiency scores, and data quality for using machine learning classifiers. Discussion: Shambhala-2 harmonizer demonstrated the best results with the close to 1 correlation and linear regression coefficients for the comparison of training vs validation datasets and more than two times lesser instability for calculation of drug efficiency scores compared to other methods.

Identification of Proximity Interactors of Mammalian Nucleoid Proteins by BioID

Mitochondrial nucleoids are compact nucleoprotein complexes, in which mtDNA is located, replicated, and transcribed. Several proteomic approaches have been previously employed to identify nucleoid proteins; however, a consensus list of nucleoid-associated proteins has not been generated. Here we describe a proximity-biotinylation assay, BioID, which allows identification of proximity interactors of mitochondrial nucleoid proteins. It uses a promiscuous biotin ligase fused to a protein of interest which covalently attaches biotin to lysine residues of its proximal neighbors. Biotinylated proteins can be further enriched by a biotin-affinity purification and identified by mass-spectrometry. BioID can identify transient and weak interactions and can be used to identify changes in the interactions upon different cellular treatments, for different protein isoforms or for pathogenic variants.

Deleterious variants in CRLS1 lead to cardiolipin deficiency and cause an autosomal recessive multi-system mitochondrial disease

Mitochondrial diseases are a group of inherited diseases with highly varied and complex clinical presentations. Here, we report four individuals, including two siblings, affected by a progressive mitochondrial encephalopathy with biallelic variants in the cardiolipin biosynthesis gene CRLS1. Three affected individuals had a similar infantile presentation comprising progressive encephalopathy, bull's eye maculopathy, auditory neuropathy, diabetes insipidus, autonomic instability, cardiac defects and early death. The fourth affected individual presented with chronic encephalopathy with neurodevelopmental regression, congenital nystagmus with decreased vision, sensorineural hearing loss, failure to thrive and acquired microcephaly. Using patient-derived fibroblasts, we characterized cardiolipin synthase 1 (CRLS1) dysfunction that impaired mitochondrial morphology and biogenesis, providing functional evidence that the CRLS1 variants cause mitochondrial disease. Lipid profiling in fibroblasts from two patients further confirmed the functional defect demonstrating reduced cardiolipin levels, altered acyl-chain composition and significantly increased levels of phosphatidylglycerol, the substrate of CRLS1. Proteomic profiling of patient cells and mouse Crls1 knockout cell lines identified both endoplasmic reticular and mitochondrial stress responses, and key features that distinguish between varying degrees of cardiolipin insufficiency. These findings support that deleterious variants in CRLS1 cause an autosomal recessive mitochondrial disease, presenting as a severe encephalopathy with multi-systemic involvement. Furthermore, we identify key signatures in cardiolipin and proteome profiles across various degrees of cardiolipin loss, facilitating the use of omics technologies to guide future diagnosis of mitochondrial diseases.

Mitochondrial mistranslation modulated by metabolic stress causes cardiovascular disease and reduced lifespan

Changes in the rate and fidelity of mitochondrial protein synthesis impact the metabolic and physiological roles of mitochondria. Here we explored how environmental stress in the form of a high-fat diet modulates mitochondrial translation and affects lifespan in mutant mice with error-prone (Mrps12^{ep / ep} ) or hyper-accurate (Mrps12^{ha / ha} ) mitochondrial ribosomes. Intriguingly, although both mutations are metabolically beneficial in reducing body weight, decreasing circulating insulin and increasing glucose tolerance during a high-fat diet, they manifest divergent (either deleterious or beneficial) outcomes in a tissue-specific manner. In two distinct organs that are commonly affected by the metabolic disease, the heart and the liver, Mrps12^{ep / ep} mice were protected against heart defects but sensitive towards lipid accumulation in the liver, activating genes involved in steroid and amino acid metabolism. In contrast, enhanced translational accuracy in Mrps12^{ha / ha} mice protected the liver from a high-fat diet through activation of liver proliferation programs, but enhanced the development of severe hypertrophic cardiomyopathy and led to reduced lifespan. These findings reflect the complex transcriptional and cell signalling responses that differ between post-mitotic (heart) and highly proliferative (liver) tissues. We show trade-offs between the rate and fidelity of mitochondrial protein synthesis dictate tissue-specific outcomes due to commonly encountered stressful environmental conditions or aging.