Intrinsic gene changes determine the successful establishment of stable renal cancer cell lines from tumor tissue

wMKL: multi-omics data integration enables novel cancer subtype identification via weight-boosted multi-kernel learning

BACKGROUND

Cancer is a heterogeneous disease driven by complex molecular alterations. Cancer subtypes determined from multi-omics data can provide novel insight into personalised precision treatment. It is recognised that incorporating prior weight knowledge into multi-omics data integration can improve disease subtyping.

METHODS

We develop a weighted method, termed weight-boosted Multi-Kernel Learning (wMKL) which incorporates heterogeneous data types as well as flexible weight functions, to boost subtype identification. Given a series of weight functions, we propose an omnibus combination strategy to integrate different weight-related P-values to improve subtyping precision.

RESULTS

wMKL models each data type with multiple kernel choices, thus alleviating the sensitivity and robustness issue due to selecting kernel parameters. Furthermore, wMKL integrates different data types by learning weights of different kernels derived from each data type, recognising the heterogeneous contribution of different data types to the final subtyping performance. The proposed wMKL outperforms existing weighted and non-weighted methods. The utility and advantage of wMKL are illustrated through extensive simulations and applications to two TCGA datasets. Novel subtypes are identified followed by extensive downstream bioinformatics analysis to understand the molecular mechanisms differentiating different subtypes.

CONCLUSIONS

The proposed wMKL method provides a novel strategy for disease subtyping. The wMKL is freely available at https://github.com/biostatcao/wMKL .

Biological and genetic characterization of a newly established human external auditory canal carcinoma cell line, SCEACono2

Squamous cell carcinoma of the external auditory canal (EACSCC) is an extraordinarily rare and aggressive malignant disease. Establishment of EACSCC cell line with robust molecular characteristics is essential for the basic and translational research of EACSCC. Here, we show the newly established EACSCC cell line SCEACono2, derived from a patient with well-to-moderately differentiated EACSCC. We analyzed histologic and genetic features of SCEACono2 hiring multiple experiments, including next-generation sequencing (NGS). Immunocytochemical staining of SCEACono2 showed positivity of p53 and SCC1/2. Furthermore, SCEACono2 exhibited a unique characteristic that cytokeratin, vimentin as well as cancer stem cell markers (CD44, CD133, ALP and Oct3/4) were positive. SCEACono2 had an ability to form tumors at the temporal lesion xenograft nude mice model. NGS revealed that SCEACono2 harbored the somatic mutations of TP53 (p.G245S) and NOTCH1 (p.A465T). RNA-seq and downstream bioinformatics analysis revealed significant enrichment of genes involved in inflammation and cell adhesion in SCEACono2 compared to SCC-9 and HSC-4. STR profiling indicated no evidence of cross-contamination. In conclusion, SCEACono2 could serves as a promising and robust research resource of EACSCC in vitro and in vivo.

RNA Sequencing Reveals Alterations and Similarities in Cell Metabolism, Hypoxia and Immune Evasion in Primary Cell Cultures of Clear Cell Renal Cell Carcinoma

The treatment of advanced renal cell carcinoma remains a challenge. To develop novel therapeutic approaches, primary cell cultures as an in vitro model are considered more representative than commercial cell lines. In this study, we analyzed the gene expression of previously established primary cell cultures of clear cell renal cell carcinoma by bulk (3’m)RNA sequencing and compared it to the tissue of origin. The objectives were the identification of dysregulated pathways under cell culture conditions. Furthermore, we assessed the suitability of primary cell cultures for studying crucial biological pathways, including hypoxia, growth receptor signaling and immune evasion. RNA sequencing of primary cell cultures of renal cell carcinoma and a following Enrichr database analysis revealed multiple dysregulated pathways under cell culture conditions. 444 genes were significantly upregulated and 888 genes downregulated compared to the tissue of origin. The upregulated genes are crucial in DNA repair, cell cycle, hypoxia and metabolic shift towards aerobic glycolysis. A downregulation was observed for genes involved in pathways of immune cell differentiation and cell adhesion. We furthermore observed that 7275 genes have a similar mRNA expression in cell cultures and in tumor tissue, including genes involved in the immune checkpoint signaling or in pathways responsible for tyrosine kinase receptor resistance. Our findings confirm that primary cell cultures are a representative tool for specified experimental approaches. The results presented in this study give further valuable insights into the complex adaptation of patient-derived cells to a new microenvironment, hypoxia and other cell culture conditions, which are often neglected in daily research, and allow new translational and therapeutic approaches.

LAG3 (LAG-3, CD223) DNA methylation correlates with LAG3 expression by tumor and immune cells, immune cell infiltration, and overall survival in clear cell renal cell carcinoma

Background

Lymphocyte activating 3 (LAG3, LAG-3, CD223) is a promising target for immune checkpoint inhibition in clear cell renal cell carcinoma (KIRC). The aim of this study was to investigate the epigenetic regulation of LAG3 in KIRC by methylation.

Methods

We correlated quantitative LAG3 methylation levels with transcriptional activity, immune cell infiltration, and overall survival in a cohort of n=533 patients with KIRC and n=160 normal adjacent tissue (NAT) samples obtained from The Cancer Genome Atlas (TCGA). Furthermore, we analyzed LAG3 methylation in peripheral blood mononuclear cells (PBMCs) and KIRC cell lines. We validated correlations between LAG3 expression, immune cell infiltrates, survival, and methylation in an independent KIRC cohort (University Hospital Bonn (UHB) cohort, n=118) by means of immunohistochemistry and quantitative methylation-specific PCR.

Results

We found differential methylation profiles among PBMCs, NAT, KIRC cell lines, and KIRC tumor tissue. Methylation strongly correlated with LAG3 mRNA expression in KIRCs (TCGA cohort) and KIRC cell lines. In the UHB cohort, methylation correlated with LAG3-positive immune cells and tumor-intrinsic LAG3 protein expression. Furthermore, LAG3 methylation strongly correlated with signatures of distinct immune cell infiltrates, an interferon-y signature (TCGA cohort), and immunohistochemically quantified CD45⁺, CD8⁺, and CD4⁺ immune cell infiltrates (UHB cohort). LAG3 mRNA expression (TCGA cohort), methylation (both cohorts), and tumor cell-intrinsic protein expression (UHB cohort) was significantly associated with overall survival.

Conclusion

Our data suggest an epigenetic regulation of LAG3 expression in tumor and immune cells via DNA methylation. LAG3 expression and methylation is associated with a subset of KIRCs showing a distinct clinical course and immunogenicity. Our study provides rationale for further testing LAG3 DNA methylation as a predictive biomarker for response to LAG3 immune checkpoint inhibitors.

Establishment and characterization of a new human colon cancer cell line, PUMC-CRC1

Colorectal cancer (CRC) is one of the most common and fatal gastrointestinal cancers worldwide. Considering their diversity, the establishment of new continuous CRC cell lines with clear genetic backgrounds will provide useful tools for exploring molecular mechanisms, screening and evaluating antitumor drugs in CRC studies. Our de novo CRC cell line, PUMC-CRC1 (Peking Union Medical College Colorectal Cancer 1) was derived from a 47-year-old Chinese female patient diagnosed with moderately to poorly differentiated colon adenocarcinoma. Multiple experiments were used for full characterization. The new cell line was epithelial-like and was passaged for more than 40 times, with a population doubling time of 44 h in vitro, detected by cell counts. The cells exhibited complicated chromosomal abnormalities. The tumor formation rate in SCID mice was 100%. The xenograft tumor was adenocarcinoma with poor to moderate differentiation by Haematoxylin and Eosin staining (H&E) sections. Immunohistochemistry (IHC) analysis and next-generation sequencing (NGS) revealed microsatellite stable (MSS), APC (p.T1493fs) inactivation, KRAS (p.G12V) activation, and SMAD4 (p.V506A) mutation. Quality control of the cell line proved mycoplasma negative and identical STR profile with that of the original tissue, and no interspecific or intraspecific cross contamination was detected. In conclusion, PUMC-CRC1 was a newly established and well characterized human colon cancer cell line, which might be a good model for both in vitro and in vivo studies of the mechanism of colon cancer progression and the treatment strategies for MSS CRC.

FAM196B promotes proliferation and migration via regulating epithelial-mesenchymal transition in esophageal cancer

BACKGROUND

EC (esophageal cancer) is a common cancer among people in the world. The molecular mechanism of FAM196B (family with sequence similarity 196 member B) in EC is still unclear. This article aimed to clarify the role of FAM196B in EC.

METHODS

The expression of FAM196B in EC tissues was detected using qRT-PCR. The prognosis of FAM196B in EC patients was determined by log-rank kaplan-Meier survival analysis and Cox regression analysis. Furthermore, shRNA was used to knockdown the expression of FAM196B in EC cell lines. MTT, wound healing assays and western blot were used to determine the role of FAM196B in EC cells.

RESULTS

In our research, we found that the expression of FAM196B was up-regulated in EC tissues. The increased expression of FAM196B was significantly correlated with differentiation, lymph node metastasis, stage, and poor survival. The proliferation and migration of EC cells were inhibited after FAM196B-shRNA transfection in vitro and vivo. The western blot result showed that FAM196B could regulate EMT.

CONCLUSION

These results suggested that FAM196B severs as an oncogene and promotes cell proliferation and migration in EC. In addition, FAM196B may be a potential therapeutic target for EC patients.

Transmembrane protein ADAM29 facilitates cell proliferation, invasion and migration in clear cell renal cell carcinoma

Abnormal expression of ADAM29 has been frequently reported in several cancers, however, its role in clear cell renal cell carcinoma (ccRCC) has not evaluated in detail. Herein, we attempt to determine the biological role and the action mechanism of ADAM29 in ccRCC. Bioinformatics analysis based on the ccRCC RNA-Seq dataset from TCGA database revealed that ADAM29 was up-expressed in ccRCC tissues by comparison with normal tissues. And a significant increase of ADAM29 expression was also observed in 3 ccRCC cell lines (UT33A, Caki-1, and786-O) in comparison with normal cell line. Besides, high level of ADAM29 was found to be connected with the poor prognosis and could be considered as an independent prognosticator for patients with ccRCC. Furthermore, functional experiments in vitro demonstrated that ADAM29 promoted the growth, invasion and migration of ccRCC cells. Moreover, Western blot assays indicated that ADAM29 was positively correlated with the level of proliferation-related proteins Cyclin D1 and PCNA and motion-related proteins MMP9 and Snail. Our data indicate that ADAM29 acts as an oncogene that increases tumour cells proliferation, invasion and migration partly by regulating the expression of Cyclin D1/PCNA/MMP9/Snail, suggesting that ADAM29 may become a prognosticator and therapeutic candidate for ccRCC.

RNA sequencing reveals PNN and KCNQ1OT1 as predictive biomarkers of clinical outcome in stage III colorectal cancer patients treated with adjuvant chemotherapy

Five-year overall survival of stage III colorectal cancer (CRC) patients treated with standard adjuvant chemotherapy (ACHT) is highly variable. Genomic biomarkers and/or transcriptomic profiles identified lack of adequate validation. Aim of our study was to identify and validate molecular biomarkers predictive of ACHT response in stage III CRC patients by a transcriptomic approach. From a series of CRC patients who received ACHT, two stage III extreme cohorts (unfavorable vs. favorable prognosis) were selected. RNA-sequencing was performed from fresh frozen explants. Tumors were characterized for somatic mutations. Validation was performed in stage III CRC patients extracted from two GEO datasets. According to disease-free survival (DFS), 108 differentially expressed genes (104/4 up/downregulated in the unfavorable prognosis group) were identified. Among 104 upregulated genes, 42 belonged to olfactory signaling pathways, 62 were classified as pseudogenes (n = 17), uncharacterized noncoding RNA (n = 10), immune response genes (n = 4), microRNA (n = 1), cancer-related genes (n = 14) and cancer-unrelated genes (n = 16). Three out of four down-regulated genes were cancer-related. Mutational status (i.e., RAS, BRAF, PIK3CA) did not differ among the cohorts. In the validation cohort, multivariate analysis showed high PNN and KCNQ1OT1 expression predictive of shorter DFS in ACHT treated patients (p = 0.018 and p = 0.014, respectively); no difference was observed in untreated patients. This is the first study that identifies by a transcriptomic approach and validates PNN and KCNQ1OT1 as molecular biomarkers predictive of chemotherapy response in stage III CRC patients. After a further validation in an independent cohort, PNN and KCNQ1OT1 evaluation could be proposed to prospectively identify stage III CRC patients benefiting from ACHT.

Effect of FAM196B in human lung adenocarcinoma

Lung cancer is the leading cause of cancer-related mortality worldwide. The present study focused on the function of family with sequence similarity 196 member B (FAM196B) in lung adenocarcinoma (LUAD). We analyzed lung carcinoma patients' data from the Cancer Genome Atlas (TCGA), and verify the change of FAM196B expression in 30 LUAD tissues and matched adjacent non-tumor tissues (> 5 cm) by tissue microarray. To investigate the role of FAM196B in LUAD, we used lentivirus-mediated short hairpin RNA (shRNA) to knockdown FAM196B expression in human LUAD cell line A549 and H1299 and assessed it by RT-qPCR and western blot. Celigo Imaging Cytometry System, MTT assays and colony formation were conducted to evaluate cell proliferation. Cell apoptosis were assayed by Annexin V staining. We found that FAM196B was significantly upregulated (P=5.06E-06) in LUAD compared with adjacent non-tumor tissues. Cell proliferation was inhibited in FAM196B-silenced A549 and H1299 cells using Celigo Imaging Cytometry System, MTT assays and colony formation assays. Apoptosis rate was significantly increased in FAM196B-shRNA group than the control group. In conclusion, Knockdown of FAM196B can inhibit cell proliferation, cell colony formation capacity, and promote cell apoptosis in A549 and H1299 cell lines. FAM196B may be one novel potential targets for treating patients with LUAD.

Solute carrier family 34 member 2 overexpression contributes to tumor growth and poor patient survival in colorectal cancer

Solute carrier family 34 member 2 (SLC34A2) is a well-known sodium-dependent phosphate transporter that has recently been linked to cancer development. However, its specific oncogenic role remains controversial in numerous human malignancies, and is currently unknown in colorectal cancer (CRC). Therefore, in this study we firstly used Oncomine database to determine its expression in cancer tissues and found it is overexpressed in thyroid, ovarian and renal cancer, while it is opposite in lung, breast and pancreas cancer. Using qRT-PCR and western blot, we then demonstrated its overexpression in CRC tissues as compared with adjacent normal tissues (n = 20). In a retrospective cohort enrolling 190 CRC patients, we proved its expression was significantly correlated with N stage. Furthermore, high SLC34A2 expression is associated with higher postoperative metastasis rate and serves as an independent adverse factor affecting patient prognosis. In subgroup analysis, SLC34A2 expression could stratify the patient prognosis in stage II and III CRC, but failed in stage IV CRC. In cellular assays in vitro, knockdown of SLC34A2 dramatically inhibited the proliferation and colony formation, induced the apoptosis and arrests the cell cycle progression of HCT-116 CRC cells. In cellular assays in vivo, knockdown of SLC34A2 significantly inhibited the growth of xenografts, decreasing Ki-67 and proliferating cell nuclear antigen (PCNA) expression and increasing apoptosis rate. Taken together, our study indicates SLC34A2 plays a crucial promoting role in CRC development and therefore has great potential to be further developed as a reliable biomarker for CRC diagnosis and treatment.