Implications of time-series gene expression profiles of replicative senescence

Although senescence has long been implicated in aging-associated pathologies, it is not clearly understood how senescent cells are linked to these diseases. To address this knowledge gap, we profiled cellular senescence phenotypes and mRNA expression patterns during replicative senescence in human diploid fibroblasts. We identified a sequential order of gain-of-senescence phenotypes: low levels of reactive oxygen species, cell mass/size increases with delayed cell growth, high levels of reactive oxygen species with increases in senescence-associated β-galactosidase activity (SA-β-gal), and high levels of SA-β-gal activity. Gene expression profiling revealed four distinct modules in which genes were prominently expressed at certain stages of senescence, allowing us to divide the process into four stages: early, middle, advanced, and very advanced. Interestingly, the gene expression modules governing each stage supported the development of the associated senescence phenotypes. Senescence-associated secretory phenotype-related genes also displayed a stage-specific expression pattern with three unique features during senescence: differential expression of interleukin isoforms, differential expression of interleukins and their receptors, and differential expression of matrix metalloproteinases and their inhibitory proteins. We validated these phenomena at the protein level using human diploid fibroblasts and aging Sprague-Dawley rat skin tissues. Finally, disease-association analysis of the modular genes also revealed stage-specific patterns. Taken together, our results reflect a detailed process of cellular senescence and provide diverse genome-wide information of cellular backgrounds for senescence.

Embryonic Stem Cell-Derived mmu-miR-291a-3p Inhibits Cellular Senescence in Human Dermal Fibroblasts Through the TGF-β Receptor 2 Pathway

Senescent cells accumulate in various tissues over time and contribute to tissue dysfunction and aging-associated phenotypes. Accumulating evidence suggests that cellular senescence can be inhibited through pharmacological intervention, as well as through treatment with soluble factors derived from embryonic stem cells (ESCs). In an attempt to investigate the anti-senescence factors secreted by ESCs, we analyzed mouse ESC-derived extracellular microRNAs in conditioned medium via microRNA array analysis. We selected mmu-miR-291a-3p as a putative anti-senescence factor via bioinformatics analysis. We validated its inhibitory effects on replicative, Adriamycin-induced, and ionizing radiation-induced senescence in human dermal fibroblasts. Treatment of senescent cells with mmu-miR-291a-3p decreased senescence-associated β-galactosidase activity, enhanced proliferative potential, and reduced mRNA and protein expression of TGF-β receptor 2, p53, and p21. mmu-miR-291a-3p in conditioned medium was enclosed in ESC-derived exosomes and exosomes purified from ESC conditioned medium inhibited cellular senescence. The inhibitory effects of mmu-miR-291a-3p were mediated through the TGF-β receptor 2 signaling pathway. Hsa-miR-371a-3p and hsa-miR-520e, the human homologs of mmu-miR-291a-3p, showed similar anti-senescence activity. Furthermore, mmu-miR-291a-3p accelerated the excisional skin wound healing process in aged mice. Our results indicate that the ESC-derived mmu-miR-291a-3p is a novel candidate agent that can be utilized for cell-free therapeutic intervention against aging and aging-related diseases.

Identification of a mitochondrial defect gene signature reveals NUPR1 as a key regulator of liver cancer progression

Many cancer cells require more glycolytic adenosine triphosphate production due to a mitochondrial respiratory defect. However, the roles of mitochondrial defects in cancer development and progression remain unclear. To address the role of transcriptomic regulation by mitochondrial defects in liver cancer cells, we performed gene expression profiling for three different cell models of mitochondrial defects: cells with chemical respiratory inhibition (rotenone, thenoyltrifluoroacetone, antimycin A, and oligomycin), cells with mitochondrial DNA depletion (Rho0), and liver cancer cells harboring mitochondrial defects (SNU354 and SNU423). By comparing gene expression in the three models, we identified 10 common mitochondrial defect-related genes that may be responsible for retrograde signaling from cancer cell mitochondria to the intracellular transcriptome. The concomitant expression of the 10 common mitochondrial defect genes is significantly associated with poor prognostic outcomes in liver cancers, suggesting their functional and clinical relevance. Among the common mitochondrial defect genes, we found that nuclear protein 1 (NUPR1) is one of the key transcription regulators. Knockdown of NUPR1 suppressed liver cancer cell invasion, which was mediated in a Ca(2+) signaling-dependent manner. In addition, by performing an NUPR1-centric network analysis and promoter binding assay, granulin was identified as a key downstream effector of NUPR1. We also report association of the NUPR1-granulin pathway with mitochondrial defect-derived glycolytic activation in human liver cancer.

CONCLUSION

Mitochondrial respiratory defects and subsequent retrograde signaling, particularly the NUPR1-granulin pathway, play pivotal roles in liver cancer progression.

Transcriptomic and histopathological analysis of cholangiolocellular differentiation trait in intrahepatic cholangiocarcinoma

BACKGROUND & AIMS

Intrahepatic cholangiocarcinoma (iCCA) is a heterogeneous entity with diverse aetiologies, morphologies and clinical outcomes. Recently, histopathological distinction of cholangiolocellular differentiation (CD) of iCCA has been suggested. However, its genome-wide molecular features and clinical significance remain unclear.

METHODS

Based on CD status, we stratified iCCAs into iCCA with CD (n=20) and iCCA without CD (n=102), and performed an integrative analysis using transcriptomic and clinicopathological profiles.

RESULTS

iCCA with CD revealed less aggressive histopathological features compared to iCCA without CD, and iCCA with CD showed favourable clinical outcomes of overall survival and time to recurrence than iCCA without CD (P<.05 for all). Transcriptomic profiling revealed that iCCA with CD resembled an inflammation-related subtype, while iCCA without CD resembled a proliferation subtype. In addition, we identified a CD signature that can predict prognostic outcomes of iCCA (CD_UP, n=486 and CD_DOWN, n=308). iCCAs were subgrouped into G1 (positivity for CRP and CDH2, 7%), G3 (positivity for S100P and TFF1, 32%) and G2 (the others, 61%). Prognostic outcomes for overall survival (P=.001) and time to recurrence (P=.017) were the most favourable in G1-iCCAs, intermediate in G2-iCCAs and the worst in G3-iCCAs. Similar result was confirmed in the iCCA set from GSE26566 (n=68).

CONCLUSIONS

CD signature was identified to predict the prognosis of iCCA. The combined evaluation of histology of CD and protein expression status of CRP, CDH2, TFF1 and S100P might help subtyping and predicting clinical outcomes of iCCA.

The Ubiquitin-like with PHD and Ring Finger Domains 1 (UHRF1)/DNA Methyltransferase 1 (DNMT1) Axis Is a Primary Regulator of Cell Senescence

As senescence develops, cells sequentially acquire diverse senescent phenotypes along with simultaneous multistage gene reprogramming. It remains unclear what acts as the key regulator of the collective changes in gene expression at initiation of senescent reprogramming. Here we analyzed time series gene expression profiles obtained in two different senescence models in human diploid fibroblasts: replicative senescence and H₂O₂-induced senescence. Our results demonstrate that suppression of DNA methyltransferase 1 (DNMT1)-mediated DNA methylation activity was an initial event prior to the display of senescent phenotypes. We identified seven DNMT1-interacting proteins, ubiquitin-like with PHD and ring finger domains 1 (UHRF1), EZH2, CHEK1, SUV39H1, CBX5, PARP1, and HELLS (also known as LSH (lymphoid-specific helicase) 1), as being commonly down-regulated at the same time point as DNMT1 in both senescence models. Knockdown experiments revealed that, among the DNMT1-interacting proteins, only UHRF1 knockdown suppressed DNMT1 transcription. However, UHRF1 overexpression alone did not induce DNMT1 expression, indicating that UHRF1 was essential but not sufficient for DNMT1 transcription. Although UHRF1 knockdown effectively induced senescence, this was significantly attenuated by DNMT1 overexpression, clearly implicating the UHRF1/DNMT1 axis in senescence. Bioinformatics analysis further identified WNT5A as a downstream effector of UHRF1/DNMT1-mediated senescence. Senescence-associated hypomethylation was found at base pairs -1569 to -1363 from the transcription start site of the WNT5A gene in senescent human diploid fibroblasts. As expected, WNT5A overexpression induced senescent phenotypes. Overall, our results indicate that decreased UHRF1 expression is a key initial event in the suppression of DNMT1-mediated DNA methylation and in the consequent induction of senescence via increasing WNT5A expression.

Dynamics of Genomic, Epigenomic, and Transcriptomic Aberrations during Stepwise Hepatocarcinogenesis

Hepatocellular carcinoma (HCC) undergoes a stepwise progression from liver cirrhosis to low-grade dysplastic nodule (LGDN), high-grade dysplastic nodule (HGDN), early HCC (eHCC), and progressed HCC (pHCC). Here, we profiled multilayered genomic, epigenomic, and transcriptomic aberrations in the stepwise hepatocarcinogenesis. Initial DNA methylation was observed in eHCC (e.g., DKK3, SALL3, and SOX1) while more extensive methylation was observed in pHCC. In addition, eHCCs showed an initial loss of DNA copy numbers of tumor suppressor genes in the 4q and 13q regions, thereby conferring survival benefits to cancer cells. Transcriptome analysis revealed that HGDNs expressed endoplasmic reticulum (ER) stress-related genes, while eHCC started to express oncogenes. Furthermore, integrative analysis indicated that expression of the serine peptidase inhibitor, Kazal type 1 (SPINK1), played a pivotal role in eHCC development. Significant demethylation of SPINK1 was observed in eHCC compared to HGDN. The study also demonstrated that ER stress may induce SPINK1 demethylation and expression in liver cancer cells. In conclusion, these results reveal the dynamics of multiomic aberrations during malignant conversion of liver cancer, thus providing novel pathobiological insights into hepatocarcinogenesis. SIGNIFICANCE: Multiomics profiling and integrative analyses of stepwise hepatocarcinogenesis reveal novel mechanistic and clinical insights into hepatocarcinogenesis.

Molecular classification of basal cell carcinoma of skin by gene expression profiling

Non-melanoma skin cancers (NMSC) including basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are more common kinds of skin cancer. Although these tumors share common pathological and clinical features, their similarity and heterogeneity at molecular levels are not fully elaborated yet. Here, by performing comparative analysis of gene expression profiling of BCC, SCC, and normal skin tissues, we could classify the BCC into three subtypes of classical, SCC-like, and normal-like BCCs. Functional enrichment and pathway analyses revealed the molecular characteristics of each subtype. The classical BCC showed the enriched expression and transcription signature with the activation of Wnt and Hedgehog signaling pathways, which were well known key features of BCC. By contrast, the SCC-like BCC was enriched with immune-response genes and oxidative stress-related genes. Network analysis revealed the PLAU/PLAUR as a key regulator of SCC-like BCC. The normal-like BCC showed prominent activation of metabolic processes particularly the fatty acid metabolism. The existence of these molecular subtypes could be validated in an independent dataset, which demonstrated the three subgroups of BCC with distinct functional enrichment. In conclusion, we suggest a novel molecular classification of BCC providing insights on the heterogeneous progression of BCC.

Perspectives of integrative cancer genomics in next generation sequencing era

The explosive development of genomics technologies including microarrays and next generation sequencing (NGS) has provided comprehensive maps of cancer genomes, including the expression of mRNAs and microRNAs, DNA copy numbers, sequence variations, and epigenetic changes. These genome-wide profiles of the genetic aberrations could reveal the candidates for diagnostic and/or prognostic biomarkers as well as mechanistic insights into tumor development and progression. Recent efforts to establish the huge cancer genome compendium and integrative omics analyses, so-called "integromics", have extended our understanding on the cancer genome, showing its daunting complexity and heterogeneity. However, the challenges of the structured integration, sharing, and interpretation of the big omics data still remain to be resolved. Here, we review several issues raised in cancer omics data analysis, including NGS, focusing particularly on the study design and analysis strategies. This might be helpful to understand the current trends and strategies of the rapidly evolving cancer genomics research.

Dissection of relationship between small heat shock proteins and mycobacterial diseases

Mycobacteria belong to a genus which has membership ranging from saprophytes to deadly pathogens that cause several infectious diseases affecting a large population of the world. Among them, tuberculosis and leprosy are the major granulomatous mycobacterial diseases. While there are successes and failures in the fight against these infections, mechanisms of pathogenesis continue to be a challenge to clinicians and biologists alike. Though it is known that both host and bacterial factors are important in the pathogenicity versus protection, all the triggers and responses are not known. Among various bacterial factors, small heat shock proteins (sSHPs) could be important targets for drug development, immunomodulation and serodiagnosis. sSHPs are the molecular chaperones that are believed to act as mantle for the mycobacteria against host's immune attack and facilitate the survival of pathogen in host body. Best studied small heat shock proteins in M. tuberculosis are sSHP16.3 and Acr2 while in M. leprae, it is 18 kD protein antigen. In this review, works on various aspects of small heat shock proteins which fall in 10 to 19 kD range have been summarized and some thoughts about future road-map have been put into.