Microtubule associated protein 9 inhibits liver tumorigenesis by suppressing ERCC3

Integrative multiomic analysis identifies distinct molecular subtypes of NAFLD in a Chinese population

Nonalcoholic fatty liver disease (NAFLD) has become a common health care burden worldwide. The high heterogeneity of NAFLD remains elusive and impairs outcomes of clinical diagnosis and pharmacotherapy. Several NAFLD classifications have been proposed on the basis of clinical, genetic, alcoholic, or serum metabolic analyses. Yet, accurately predicting the progression of NAFLD to cirrhosis or hepatocellular carcinoma (HCC) in patients remains a challenge. Here, on the basis of a Chinese cohort of patients, we classified NAFLD into three distinct molecular subtypes (NAFLD-mSI, NAFLD-mSII, and NAFLD-mSIII) using integrative multiomics including whole-genome sequencing (WGS), proteomics, phosphoproteomics, lipidomics, and metabolomics across a broad range of liver, blood, and urine specimens. We found that NAFLD-mSI had higher expression of CYP1A2 and CYP3A4, which alleviate hepatic steatosis through mediating free fatty acid/bile acid-mTOR-FXR/PPARα signaling. NAFLD-mSII displayed an elevated risk of liver cirrhosis along with increased hepatic infiltration of M1 and M2 macrophages because of lipid-triggered hepatic CCL2 and CRP production. NAFLD-mSIII exhibited a potential risk for HCC development by increased transcription of CEBPB- and ERCC3-regulated oncogenes because of activation of the EGF-EGFR/CHKA/PI3K-PDK1-AKT cascade. Next, we validated the existence of these three NAFLD molecular subtypes in an external cohort comprising 92 patients with NAFLD across three different Chinese hospitals. These findings may aid in understanding the molecular features underlying NAFLD heterogeneity, thereby facilitating clinical diagnosis and treatment strategies with the aim of preventing the development of liver cirrhosis and HCC.

Etiologically Significant microRNAs in Hepatitis B Virus-Induced Hepatocellular Carcinoma

Hepatitis B virus (HBV) infection has been causally linked to hepatocellular carcinoma (HCC) in more than 50% cases. MicroRNAs (miRNAs) play cross-cutting mechanistic roles in the complex interplay between viral pathogenesis, host survival, and clinical outcomes. The present study set out to identify etiologically significant human miRNAs associated with HBV infection in liver-related pathologies leading to HCC. In diverse tissue types, we assembled 573 miRNAs differentially expressed in HBV-associated liver pathologies, HBV infection, fibrosis, cirrhosis, acute on chronic liver failure, and HCC. Importantly, 43 human differentially expressed miRNAs (hDEmiRs) were regulated in serum/plasma and liver tissue of patients with HBV-positive conditions. However, only two hDEmiRs, hsa-miR-21-5p and hsa-miR-143-3p, were regulated across all disease conditions. To shortlist the functional miRNAs in HBV-induced HCC pathogenesis, a reverse bioinformatics analysis was performed using eight GEO datasets and the TCGA database containing the list of differentially regulated mRNAs in HCC. A comparative study using these data with the identified targets of hDEmiRs, a set of unidirectionally regulated hDEmiRs with the potential to modulate mRNAs in HCC, were found. Moreover, our study identified five miRNAs; hsa-miR-98-5p, hsa-miR-193b-3p, hsa-miR-142-5p, hsa-miR-522-5p, and hsa-miR-370-3p targeting PIGC, KNTC1, CSTF2, SLC41A2, and RAB17, respectively, in HCC. These hDEmiRs and their targets could be pivotal in HBV infection and subsequent liver pathologies modulating HCC clinical progression. HBV infection is the largest contributor to HCC, and the present study comprises the first of its kind compendium of hDEmiRs related to HBV-related pathologies.

Molecular characterization of MAP9 in the photoreceptor sensory cilia as a modifier in canine RPGRIP1-associated cone-rod dystrophy

Photoreceptors possess a highly specialized primary cilium containing expanded ciliary membrane discs called the outer segment. The photoreceptor cilium is essential for the maintenance of the outer segment, and pathogenic variants in more than 50 cilia-related genes have been identified as causing non-syndromic inherited retinal diseases in patients. The retinitis pigmentosa GTPase regulator interacting protein 1 (RPGRIP1) is a structural protein localized to the photoreceptor cilium and biallelic RPGRIP1 variants have been associated with non-syndromic human inherited retinal diseases. In a canine cone-rod dystrophy model, a naturally occurring 44-bp exonic insertion in RPGRIP1 (RPGRIP1ins44/ins44) is the primary disease locus while an additional homozygous variant in MAP9 (microtubule associated protein 9) (MAP9aff/aff) acts as a modifier associated with early disease onset. MAP9 was originally identified as a microtubule-binding protein stabilizing microtubule structure during both mitosis and interphase in human cell lines. However, the roles of MAP9 in primary cilia, including photoreceptor neurosensory cilia, have not been well understood. Hence, we characterized the pathogenic phenotypes associated with homozygous MAP9 variant, and investigated the molecular function of MAP9 in primary cilia using the RPGRIP1-associated oligogenic canine cone-rod dystrophy model as well as cultured cells. Both functionally and structurally, the RPGRIP1ins44/ins44 MAP9aff/aff retina exhibited progressive cone photoreceptor degeneration starting earlier than the retina affected by RPGRIP1ins44/ins44 alone. Based on immunostaining of canine retinal sections and cultured cells, we found that MAP9 is prominently localized in the basal body of primary cilia and played an important role in maintaining the structure of ciliary microtubule axoneme. These findings suggest that the affected MAP9, together with mutant RPGRIP1, is deprived of critical roles in cilia organization and maintenance resulting in altered cilia structure and function giving rise to early onset and accelerated disease progression in the RPGRIP1ins44/ins44 MAP9aff/aff double homozygote cone-rod dystrophy canine model.

In-depth quantitative proteomics analysis revealed C1GALT1 depletion in ECC-1 cells mimics an aggressive endometrial cancer phenotype observed in cancer patients with low C1GALT1 expression

BACKGROUND

Endometrial cancer (EC) is the most common cancer of the female reproductive organs. Despite the good overall prognosis of most low-grade ECs, FIGO I and FIGO II patients might experience tumor recurrence and worse prognosis. The study of alterations related to EC pathogenesis might help to get insights into underlying mechanisms involved in EC development and progression.

METHODS

Core tumoral samples were used to investigate the role of C1GALT1 in EC by immunohistochemistry (IHC). ECC-1 cells were used as endometrioid EC model to investigate the effect of C1GALT1 depletion using C1GALT1 specific shRNAs. SILAC quantitative proteomics analyses and cell-based assays, PCR, qPCR, WB, dot-blot and IHC analyses were used to identify, quantify and validate dysregulation of proteins.

RESULTS

Low C1GALT1 protein expression levels associate to a more aggressive phenotype of EC. Out of 5208 proteins identified and quantified by LC-MS/MS, 100 proteins showed dysregulation (log2fold-change ≥ 0.58 or ≤-0.58) in the cell protein extracts and 144 in the secretome of C1GALT1 depleted ECC-1 cells. Nine dysregulated proteins were validated. Bioinformatics analyses pointed out to an increase in pathways associated with an aggressive phenotype. This finding was corroborated by loss-of-function cell-based assays demonstrating higher proliferation, invasion, migration, colony formation and angiogenesis capacity in C1GALT1 depleted cells. These effects were associated to the overexpression of ANXA1, as demonstrated by ANXA1 transient silencing cell-based assays, and thus, correlating C1GALT and ANXA1 protein expression and biological effects. Finally, the negative protein expression correlation found by proteomics between C1GALT1 and LGALS3 was confirmed by IHC.

CONCLUSION

C1GALT1 stably depleted ECC-1 cells mimic an EC aggressive phenotype observed in patients and might be useful for the identification and validation of EC markers of progression.

Comprehensive analysis of ERCC3 prognosis value and ceRNA network in AML

BACKGROUND

Acute myeloid leukemia (AML) is a hematological malignancy with high molecular and clinical heterogeneity, and is the most common type of acute leukemia in adults. Due to limited treatment options, AML is prone to relapse and has a poor prognosis. Excision repair cross-complementing 3 (ERCC3) is an important member of nucleotide excision repair (NER) that is overexpressed in types of solid cancers and potentially regarded as a prognostic factor. However, its role in AML remains unclear. The purpose of this study was to explore ERCC3 expression and functions in AML.

METHODS

The Cancer Genome Atlas (TCGA) and GEO (Gene Expression Omnibus) were used to test the accuracy of ERCC3 expression levels for AML diagnosis. Using online databases and R packages, we also explored the signaling pathway, epigenetic regulation, infiltration of immune cells, clinical prognostic value, and ceRNA network in AML.

RESULTS

Our results revealed that ERCC3 expression was increased in AML and that high ERCC3 expression had good value for disease-free survival and overall survival in AML patients who underwent allogeneic hematopoietic stem cell transplantation (allo-HSCT). We found that ERCC3 and co-expressed genes were mainly involved in chemical carcinogenesis/reactive oxygen species, ubiquitin-mediated protein degradation and oxidative phosphorylation. In addition, almost all the m⁶A-related coding genes (except GF2BP1) were positively associated with ERCC3 expression. We also constructed a ceRNA regulatory network containing ERCC3 in AML and identified 6 pairs of ceRNA networks, indicating that ERCC3 expression is regulated by a noncoding RNA system.

CONCLUSION

This study demonstrated that ERCC3 was overexpressed in AML and that high ERCC3 expression can be considered a biomarker conducive to allo-HSCT in AML patients.

Identification and preliminary validation of synovial tissue-specific genes and their-mediated biological mechanisms in rheumatoid arthritis

BACKGROUND

Rheumatoid arthritis (RA) is a chronic inflammatory autoimmune disease. It is well known that the formation of positive feedback between synovial hyperplasia and inflammatory infiltration is intimately associated with the occurrence and development of RA. However, the exact mechanisms still remain unknown, making the early diagnosis and therapy of RA difficult. This study was designed to identify prospective diagnostic and therapeutic biomarkers, as well as their-mediated biological mechanisms in RA.

METHODS

Three microarray datasets (GSE36700, GSE77298 and GSE153015) and two RNA-sequencing datasets (GSE89408 and GSE112656) of synovial tissues, as well as three other microarray datasets (GSE101193, GSE134087 and GSE94519) of peripheral blood were downloaded for integrated analysis. The differently expressed genes (DEGs) were identified by "limma" package of R software. Then, weight gene co-expression analysis and gene set enrichment analysis were performed to investigate synovial tissue-specific genes and their-mediated biological mechanisms in RA. The expression of candidate genes and their diagnostic value for RA were verified by quantitative real-time PCR and receiver operating characteristic (ROC) curve, respectively. Relevant biological mechanisms were explored through cell proliferation and colony formation assay. The suggestive anti-RA compounds were discovered by CMap analysis.

RESULTS

We identified a total of 266 DEGs, which were mainly enriched in cellular proliferation and migration, infection and inflammatory immune signaling pathways. Bioinformatics analysis and molecular validation revealed 5 synovial tissue-specific genes, which exhibited excellent diagnostic value for RA. The infiltration level of immune cells in RA synovial tissue was significantly higher than that in control individuals. Moreover, preliminary molecular experiments suggested that these characteristic genes may be responsible for the high proliferation potential of RA fibroblast-like synoviocytes (FLSs). Finally, 8 small molecular compounds with anti-RA potential were obtained.

CONCLUSIONS

We have proposed 5 potential diagnostic and therapeutic biomarkers (CDK1, TTK, HMMR, DLGAP5, and SKA3) in synovial tissues that may contribute to the pathogenesis of RA. These findings may shed light on the early diagnosis and therapy of RA.

Proteomic profiling reveals neuronal ion channel dysregulation and cellular responses to DNA damage-induced cell cycle arrest and senescence in human neuroblastoma SH-SY5Y cells exposed to cypermethrin

Cypermethrin (CYP), a synthetic pyrethroid of class II, is widely used as a pesticide worldwide. The primary target of cypermethrin is a voltage-gated sodium channel. The neurotoxicity of CYP has been extensively studied in terms of affecting neuronal development, increasing cellular oxidative stress, and apoptosis. However, little is known about how it affects the expression of channel proteins involved in synaptic transmission, as well as the effects of cypermethrin on DNA damage and cell cycle processes. We found that the ligand and voltage-gated calcium channels and proteins involved in synaptic transmission including NMDA 1 receptor subunit, alpha 1A-voltage-dependent calcium channel, synaptotagmin-17, and synaptojanin-2 were downregulated in CYP-treated cells. After 48h of CYP exposure, cell viability was reduced with flattened and enlarged morphology. The levels of 23 proteins regulating cell cycle processes were altered in CYP-treated cells, according to a proteomic study. The cell cycle analysis showed elevated G₀/G₁ cell cycle arrest and DNA fragmentation at the sub-G₀ stage after CYP exposure. CYP treatment also increased senescence-associated β-galactosidase positive cells, DNA damage, and apoptotic markers. Taken together, the current study showed that cypermethrin exposure caused DNA damage and hastened cellular senescence and apoptosis via disrupting cell cycle regulation. In addition, despite its primary target sodium channel, CYP might cause synaptic dysfunction via the downregulation of synaptic proteins and dysregulation of synapse-associated ion channels.

Genome-wide DNA methylation profiling of stomach cancer in the ethnic population of Mizoram, North East India

Stomach cancer is the fifth most common cancer in terms of prevalence and incidence and the fourth leading cause of mortality in men and women worldwide. It is well-established that aberrant DNA methylation in cells can lead to carcinogenesis. The primary objective of our study was to investigate the aberrant DNA methylation status of genes associated with stomach cancer with a particular reference to the ethnic population of Mizoram, North East India. The site-level analysis identified 2883 CpG sites differentially methylated, representing ~922 genes. Out of which 476 Differentially Methylated Positions (DMPs) were promoter-associated, 452 DMPs were hypermethylated, and 24 were hypomethylated. The region-level analysis identified 462 Differentially Methylated Regions (DMRs) corresponding to ~320 genes, of which ~281 genes were hypermethylated and ~ 40 genes were hypomethylated. TCGA analysis showed that some of the genes had been previously implicated in other cancers including stomach cancer. Five hypermethylated genes were selected as candidate genes for further investigations and they have shown to be novel and could serve as candidate hypermethylation biomarkers for stomach cancer in this particular ethnic group.

MetastaSite: Predicting metastasis to different sites using deep learning with gene expression data

Deep learning has massive potential in predicting phenotype from different omics profiles. However, deep neural networks are viewed as black boxes, providing predictions without explanation. Therefore, the requirements for these models to become interpretable are increasing, especially in the medical field. Here we propose a computational framework that takes the gene expression profile of any primary cancer sample and predicts whether patients' samples are primary (localized) or metastasized to the brain, bone, lung, or liver based on deep learning architecture. Specifically, we first constructed an AutoEncoder framework to learn the non-linear relationship between genes, and then DeepLIFT was applied to calculate genes' importance scores. Next, to mine the top essential genes that can distinguish the primary and metastasized tumors, we iteratively added ten top-ranked genes based upon their importance score to train a DNN model. Then we trained a final multi-class DNN that uses the output from the previous part as an input and predicts whether samples are primary or metastasized to the brain, bone, lung, or liver. The prediction performances ranged from AUC of 0.93-0.82. We further designed the model's workflow to provide a second functionality beyond metastasis site prediction, i.e., to identify the biological functions that the DL model uses to perform the prediction. To our knowledge, this is the first multi-class DNN model developed for the generic prediction of metastasis to various sites.

The origin of bladder cancer from mucosal field effects

Whole-organ mapping was used to study molecular changes in the evolution of bladder cancer from field effects. We identified more than 100 dysregulated pathways, involving immunity, differentiation, and transformation, as initiators of carcinogenesis. Dysregulation of interleukins signified the involvement of inflammation in the incipient phases of the process. An aberrant methylation/expression of multiple HOX genes signified dysregulation of the differentiation program. We identified three types of mutations based on their geographic distribution. The most common were mutations restricted to individual mucosal samples that targeted uroprogenitor cells. Two types of mutations were associated with clonal expansion and involved large areas of mucosa. The α mutations occurred at low frequencies while the β mutations increased in frequency with disease progression. Modeling revealed that bladder carcinogenesis spans 10-15 years and can be divided into dormant and progressive phases. The progressive phase lasted 1-2 years and was driven by β mutations.

MAP9 Exhibits Protumor Activities and Immune Escape toward Bladder Cancer by Mediating TGF-β1 Pathway

To investigate more potential targets for the treatment of human bladder cancer, quantitative reverse transcription polymerase chain reaction (qRT-PCR) and high-content screening (HCS) analysis were performed, and microtubule-associated protein 9 (MAP9), which had the strongest proliferation inhibition from 809 downregulated genes, has been selected. MAP9 is responsible for bipolar spindle assembly and is involved in the progression of many types of tumors; however, its role in bladder cancer (BC) remains unknown. Expressive levels of MAP9 in BC tissues were determined through immunohistochemistry, and the clinical significance of MAP9 in BC was analyzed. Short hairpin ribonucleic acid- (ShRNA-) MAP9 was used to construct stable MAP9 knockdown BC cell lines. The proliferative abilities of MAP9 were measured through assays in vivo and in vitro, and the migrated and invasive abilities of MAP9 were analyzed via in vitro experiments. Quantitative reverse transcription PCR, western blotting, coimmunoprecipitation (Co-IP), and rescue assays were used to identify downstream targets of MAP9. MAP9 expression increased in the tumor tissues, and its increased level was negatively correlated with prognosis. Further, the loss of MAP9 caused decreased BC cell proliferation via inducing the growth 1/synthesis (G1/S) cell cycle arrest in vitro and slowed tumor growth in vivo. In addition, MAP9 silencing attenuated BC cell migration and invasion. Moreover, we found that the growth 1/synthesis (G1/S) cell cycle-related genes and the epithelial mesenchymal transition (EMT) marker levels decreased after silencing MAP9. Finally, we found that the transforming growth factor beta 1 (TGF-β1) pathway is activated as a mediator for MAP9 to regulate genes related to the G1/S cell cycle and EMT. MAP9 promotes BC progression and immune escape activity through the TGF-β1 pathway and is a potential novel target for therapies of BC.

Overexpression of ERCC3 is associated with poor prognosis in patients with pancreatic cancer

Pancreatic cancer is associated with poor prognosis due to limited therapeutic options. Excision repair cross-complementing 3 (ERCC3) is an important member of nucleotide excision repair (NER) that is overexpressed in some cancers and may be regarded as a poor prognostic factor. Yet, its role in pancreatic cancer remains unclear. This study aimed to investigate the expression and functions of ERCC3 in pancreatic cancer patients and its relation with clinicopathological features. Our data suggested that the protein expression level of ERCC3 was higher in tumor tissues than in adjacent tissues. In addition, the expression of ERCC3 has shown to be associated with the tumor extent (p=0.035). Besides, analysis of the dataset in The Cancer Genome Atlas (TCGA) revealed that high expression of ERCC3 was associated with poor overall survival in pancreatic cancer patients (p=0.0136). In Cox regression analysis, ERCC3 was an independent prognostic factor for overall survival in pancreatic cancer (p<0.001). Furthermore, our in vitro data further suggested that the overexpression of ERCC3 significantly promoted pancreatic cancer (BxPC-3, CFPAC-1, and PANC-1 cells) proliferation, invasion, and migration. Taken together, this study suggested that high expression of ERCC3 might be a poor prognostic factor in human pancreatic cancer and might be used as a promising therapeutic target for pancreatic cancer treatment.