Chia-Lang Hsu

Functional Analysis and Characterization of Differential Coexpression Networks

Differential coexpression analysis is emerging as a complement to conventional differential gene expression analysis. The identified differential coexpression links can be assembled into a differential coexpression network (DCEN) in response to environmental stresses or genetic changes. Differential coexpression analyses have been successfully used to identify condition-specific modules; however, the structural properties and biological significance of general DCENs have not been well investigated. Here, we analyzed two independent Saccharomyces cerevisiae DCENs constructed from large-scale time-course gene expression profiles in response to different situations. Topological analyses show that DCENs are tree-like networks possessing scale-free characteristics, but not small-world. Functional analyses indicate that differentially coexpressed gene pairs in DCEN tend to link different biological processes, achieving complementary or synergistic effects. Furthermore, the gene pairs lacking common transcription factors are sensitive to perturbation and hence lead to differential coexpression. Based on these observations, we integrated transcriptional regulatory information into DCEN and identified transcription factors that might cause differential coexpression by gain or loss of activation in response to different situations. Collectively, our results not only uncover the unique structural characteristics of DCEN but also provide new insights into interpretation of DCEN to reveal its biological significance and infer the underlying gene regulatory dynamics.

Co-expression analysis identifies long noncoding RNA SNHG1 as a novel predictor for event-free survival in neuroblastoma

Despite of the discovery of protein therapeutic targets and advancement in multimodal therapy, the survival chance of high-risk neuroblastoma (NB) patients is still less than 50%. MYCN amplification is a potent driver of NB, which exerts its oncogenic activity through either activating or inhibiting the transcription of target genes. Recently, long noncoding RNAs (lncRNAs) are reported to be altered in cancers including NB. However, lncRNAs that are altered by MYCN amplification and associated with outcome in high-risk NB patients are limitedly discovered. Herein, we examined the expression profiles of lncRNAs and protein-coding genes between MYCN amplified and MYCN non-amplified NB from microarray (n = 47) and RNA-seq datasets (n = 493). We identified 6 lncRNAs in common that were differentially expressed (adjusted P ≤ 0.05 and fold change ≥ 2) and subsequently validated by RT-qPCR. The co-expression analysis reveals lncRNA, SNHG1 and coding gene, TAF1D highly co-expressed in NB. Kaplan-Meier analysis shows that higher expression of SNHG1 is significantly associated with poor patient survival. Importantly, multivariate analysis confirms high expression of SNHG1 as an independent prognostic marker for event-free survival (EFS) (HR = 1.58, P = 2.36E-02). In conclusion, our study unveils that SNHG1 is up-regulated by MYCN amplification and could be a potential prognostic biomarker for high-risk NB intervention.

Integrating transcriptomics and proteomics to show that tanshinone IIA suppresses cell growth by blocking glucose metabolism in gastric cancer cells

BackgroundTanshinone IIA (TIIA) is a diterpene quinone extracted from the plant Danshen (Salvia miltiorrhiza) used in traditional Chinese herbal medicine. It has been reported to have anti-tumor potential against several kinds of cancer, including gastric cancer. In most solid tumors, a metabolic switch to glucose is a hallmark of cancer cells, which do this to provide nutrients for cell proliferation. However, the mechanism associated with glucose metabolism by which TIIA acts on gastric cancer cells remains to be elucidated.ResultsWe found that TIIA treatment is able to significantly inhibit cell growth and the proliferation of gastric cancer in a dose-dependent manner. Using next-generation sequencing-based RNA-seq transcriptomics and quantitative proteomics-isobaric tags for relative and absolute quantification (iTRAQ), we characterized the mechanism of TIIA regulation in gastric cancer cell line AGS. In total, 16,603 unique transcripts and 102 proteins were identified. After enrichment analysis, we found that TIIA regulated genes are involved in carbohydrate metabolism, the cell cycle, apoptosis, DNA damage and cytoskeleton reorganization. Our proteomics data revealed the downregulation of intracellular ATP levels, glucose-6-phosphate isomerase and L-lactate dehydrogenase B chains by TIIA, which might work with disorders of glucose metabolism and extracellular lactate levels to suppress cell proliferation. The up-regulation of p53 and down-regulation of AKT was shown in TIIA- treated cells, which indicates the transformation of oncogenes. Severe DNA damage, cell cycle arrest at the G2/M transition and apoptosis with cytoskeleton reorganization were detected in TIIA-treated gastric cancer cells.ConclusionsCombining transcriptomics and proteomics results, we propose that TIIA treatment could lead cell stresses, including nutrient deficiency and DNA damage, by inhibiting the glucose metabolism of cancer cells. This study provides an insight into how the TIIA regulatory metabolism in gastric cancer cells suppresses cell growth, and may help improve the development of cancer therapy.

Oncoprotein ZNF322A transcriptionally deregulates alpha-adducin, cyclin D1 and p53 to promote tumor growth and metastasis in lung cancer.

ZNF322A encoding a classical Cys2His2 zinc finger transcription factor was previously revealed as a potential oncogene in lung cancer patients. However, the oncogenic role of ZNF322A and its underlying mechanism in lung tumorigenesis remain elusive. Here we show ZNF322A protein overexpression in 123 Asian and 74 Caucasian lung cancer patients. Multivariate Cox regression analysis indicated that ZNF322A was an independent risk factor for a poor outcome in lung cancer, corroborating the Kaplan–Meier results that patients with ZNF322A protein overexpression had significantly poorer overall survival than other patients. Overexpression of ZNF322A promoted cell proliferation and soft agar growth by prolonging cell cycle in S phase in multiple lung cell lines, including the immortalized lung cell BEAS-2B. In addition, ZNF322A overexpression enhanced cell migration and invasion, whereas knockdown of ZNF322A reduced cell growth, invasion and metastasis abilities in vitro and in vivo. Quantitative proteomic analysis revealed potential ZNF322A-regulated downstream targets, including alpha-adducin (ADD1), cyclin D1 (CCND1), and p53. Using luciferase promoter activity assay combined with site-directed mutagenesis and sequential chromatin immunoprecipitation-PCR assay, we found that ZNF322A could form a complex with c-Jun and cooperatively activate ADD1 and CCND1 but repress p53 gene transcription by recruiting differential chromatin modifiers, such as histone deacetylase 3, in an AP-1 element dependent manner. Reconstitution experiments indicated that CCND1 and p53 were important to ZNF322A-mediated promotion of cell proliferation, whereas ADD1 was necessary for ZNF322A-mediated cell migration and invasion. Our results provide compelling evidence that ZNF322A overexpression transcriptionally dysregulates genes involved in cell growth and motility therefore contributes to lung tumorigenesis and poor prognosis.

Mucin glycosylating enzyme GALNT2 suppresses malignancy in gastric adenocarcinoma by reducing MET phosphorylation

Glycosylation affects malignancy in cancer. Here, we report that N- acetylgalactosaminyltransferase 2 (GALNT2), an enzyme that mediates the initial step of mucin type-O glycosylation, suppresses malignant phenotypes in gastric adenocarcinoma (GCA) by modifying MET (Hepatocyte growth factor receptor) activity. GALNT2 mRNA and protein were downregulated in GCAs, and this reduction was associated with more advanced disease stage and shorter recurrence-free survival. Suppressing GALNT2 expression in GCA cells increased cell growth, migration, and invasion in vitro, and tumor metastasis in vivo. GALNT2 knockdown enhanced phosphorylation of MET and decreased expression of the Tn antigen on MET. Inhibiting MET activity with PHA665752 decreased the malignant phenotypes caused by GALNT2 knockdown in GCA cells. Our results indicate that GALNT2 suppresses the malignant potential of GCA cells and provide novel insights into the significance of O-glycosylation in MET activity and GCA progression.

Unveiling MYCN regulatory networks in neuroblastoma via integrative analysis of heterogeneous genomics data

MYCN, an oncogenic transcription factor of the Myc family, is a major driver of neuroblastoma tumorigenesis. Due to the difficulty in drugging MYCN directly, revealing the molecules in MYCN regulatory networks will help to identify effective therapeutic targets for neuroblastoma therapy. Here we perform ChIP-sequencing and small RNA-sequencing of neuroblastoma cells to determine the MYCN-binding sites and MYCN-associated microRNAs, and integrate various types of genomic data to construct MYCN regulatory networks. The overall analysis indicated that MYCN-regulated genes were involved in a wide range of biological processes and could be used as signatures to identify poor-prognosis MYCN-non-amplified patients. Analysis of the MYCN binding sites showed that MYCN principally served as an activator. Using a computational approach, we identified 32 MYCN co-regulators, and some of these findings are supported by previous studies. Moreover, we investigated the interplay between MYCN transcriptional and microRNA post-transcriptional regulations and identified several microRNAs, such as miR-124-3p and miR-93-5p, which may significantly contribute to neuroblastoma pathogenesis. We also found MYCN and its regulated microRNAs acted together to repress the tumor suppressor genes. This work provides a comprehensive view of MYCN regulations for exploring therapeutic targets in neuroblastoma, as well as insights into the mechanism of neuroblastoma tumorigenesis.

Temporal Phosphoproteome Dynamics Induced by an ATP Synthase Inhibitor Citreoviridin

Citreoviridin, one of toxic mycotoxins derived from fungal species, can suppress lung cancer cell growth by inhibiting the activity of ectopic ATP synthase, but has limited effect on normal cells. However, the mechanism of citreoviridin triggering dynamic molecular responses in cancer cells remains unclear. Here, we performed temporal phosphoproteomics to elucidate the dynamic changes after citreoviridin treatment in cells and xenograft model. We identified a total of 829 phosphoproteins and demonstrated that citreoviridin treatment affects protein folding, cell cycle, and cytoskeleton function. Furthermore, response network constructed by mathematical modeling shows the relationship between the phosphorylated heat shock protein 90 β and mitogen-activated protein kinase signaling pathway. This work describes that citreoviridin suppresses cancer cell growth and mitogen-activated protein kinase/extracellular signal-regulated kinase signaling by site-specific dephosphorylation of HSP90AB1 on Serine 255 and provides perspectives in cancer therapeutic strategies.

GALNT6 expression enhances aggressive phenotypes of ovarian cancer cells by regulating EGFR activity

Ovarian cancer is the most lethal of the gynecologic malignancies. N-acetylgalactosaminyltransferase 6 (GALNT6), an enzyme that mediates the initial step of mucin type-O glycosylation, has been reported to regulate mammary carcinogenesis. However, the expression and role of GALNT6 in ovarian cancer are still unclear. Here we showed that high GALNT6 expression correlates with increased recurrence, lymph node metastasis, and chemoresistance in ovarian endometrioid and clear cell carcinomas; and higher GALNT6 levels are significantly associated with poorer patient survivals. GALNT6 knockdown with two independent siRNAs significantly suppressed viability, migration, and invasion of ovarian cancer cells. Using phospho-RTK array and Western blot analyses, we identified EGFR as a critical target of GALNT6. GALNT6 knockdown decreased phosphorylation of EGFR, whereas GALNT6 overexpression increased the phosphorylation. Lectin pull-down assays with Vicia villosa agglutinin (VVA) indicated that GALNT6 was able to modify O-glycans on EGFR. Moreover, the GALNT6-enhanced invasive behavior was significantly reversed by erlotinib, an EGFR inhibitor. Our results suggest that GALNT6 expression is associated with poor prognosis of ovarian cancer and enhances the aggressive behavior of ovarian cancer cells by regulating EGFR activity.

Discovering pathway cross-talks based on functional relations between pathways

BackgroundIn biological systems, pathways coordinate or interact with one another to achieve a complex biological process. Studying how they influence each other is essential for understanding the intricacies of a biological system. However, current methods rely on statistical tests to determine pathway relations, and may lose numerous biologically significant relations.ResultsThis study proposes a method that identifies the pathway relations by measuring the functional relations between pathways based on the Gene Ontology (GO) annotations. This approach identified 4,661 pathway relations among 166 pathways from Pathway Interaction Database (PID). Using 143 pathway interactions from PID as testing data, the function-based approach (FBA) is able to identify 93% of pathway interactions, better than the existing methods based on the shared components and protein-protein interactions. Many well-known pathway cross-talks are only identified by FBA. In addition, the false positive rate of FBA is significantly lower than others via pathway co-expression analysis.ConclusionsThis function-based approach appears to be more sensitive and able to infer more biologically significant and explainable pathway relations.

Identification of lncRNA functions in lung cancer based on associated protein-protein interaction modules

Long non-coding RNAs (lncRNAs) have been found to play important roles in various biological processes; however, many of their functions remain unclear. In this study, we present a novel approach to identify the lncRNA-associated protein-protein interaction (PPI) modules and ascertain their functions in human lung squamous cell carcinoma. We collected lncRNA and mRNA expression profiles of lung squamous cell carcinoma from The Cancer Genome Atlas. To identify the lncRNA-associated PPI modules, lncRNA-mRNA co-expression networks were first constructed based on the mutual ranks of expression correlations. Next, we examined whether the co-expressed mRNAs of a specific lncRNA were closely connected by PPIs. For this, a significantly connected mRNA set was considered to be the lncRNA-associated PPI module. Finally, the prospective functions of a lncRNA was inferred using Gene Ontology enrichment analysis on the associated module. We found that lncRNA-associated PPI modules were subtype-dependent and each subtype had unique molecular mechanisms. In addition, antisense lncRNAs and sense genes tended to be functionally associated. Our results might provide new directions for understanding lncRNA regulations in lung cancer. The analysis pipeline was implemented in a web tool, available at http://lncin.ym.edu.tw/.