Comprehensive gene expression analysis of 5'-end of mRNA identified novel intronic transcripts associated with hepatocellular carcinoma

AXL transcriptionally up-regulates TMEM14A expression to mediate cell proliferation in non-small-cell lung cancer cells

In non-small cell lung cancer (NSCLC), the receptor tyrosine kinase AXL has been identified as a potent activator of tumor progression and resistance to therapies. However, the molecular mechanisms behind AXL-mediated oncogenesis remain elusive. Current study thus aimed to uncover potential downstream genes regulated by AXL in NSCLC. Through transcriptomic RNA sequencing of AXL-silenced NSCLC cells, TMEM14A was identified as a significantly up-regulated gene. Clinical evaluations using GEPIA2 revealed that TMEM14A mRNA expression was notably higher in lung adenocarcinoma (LUAD) tumor tissues compared to normal tissues. Further, significantly increased TMEM14A levels were associated with poorer overall survival in LUAD patients. Experimentally, silencing TMEM14A in NSCLC cells led to reduced cellular proliferation and ATP levels, highlighting a key role of TMEM14A in NSCLC progression. Moreover, our promoter analysis demonstrated that AXL-mediated regulation of TMEM14A transcription could involve binding of transcription factors STAT and NF-κB to 5'-promoter of TMEM14A. Collectively, current study unveils TMEM14A as a novel downstream target of AXL, suggesting its potential as a therapeutic target to counteract resistance in future NSCLC patients undergoing AXL-targeted therapies.

Genomic deciphering of sex determination and unique immune system of a potential model species rare minnow (Gobiocypris rarus)

Gobiocypris rarus is sensitive to environmental pollution, especially to heavy metal and grass carp reovirus (GCRV). Hence, it has potential utility as a biological monitor. Genetic deciphering of its unique immune system will advance our understanding of its unique adaptive strategies, which provide cues for its better application. A de novo genome of rare minnow was obtained, and its sex determination mechanism is ZZ/ZW. We identified several specific mutation genes and specific lost genes of rare minnow, and these might be related to the sensitivity of rare minnow to environmental stimuli. We also analyzed the gene expression level of different organs/tissues and found that several IFIT genes may play key roles in GCRV resistance. In addition, knockout of the gene PCDH10L indicates that PCDH10L affects Pb²⁺-induced mortality in rare minnow. Rare minnow is ready for genetic manipulation and shows potential as an emerging experimental model.

Ginsenoside Rh2 Suppresses Breast Cancer Cell Proliferation by Epigenetically Regulating the Long Noncoding RNA C3orf67-AS1

Ginsenoside Rh2, a major bioactive ingredient abundant in red ginseng, has an antiproliferative effect on various cancer cells. In this study, we report a novel long noncoding RNA, C3orf67-AS1, which was identified as being hypermethylated at a CpG site of the promoter by Rh2 in MCF-7 cancer cells. Rh2-induced hypermethylation was responsible for the lower gene expression; the expression was recovered following treatment with a methyltransferase inhibitor, 5-aza-2^'-deoxycytidine. When C3orf67-AS1 was downregulated by a siRNA, the cell growth rate was decreased, demonstrating the RNA's oncogenic activity. Accordingly, breast cancer patients showed a lower methylation and higher expression level of C3orf67-AS1. Within 800 kb flanking C3orf67-AS1 on the chromosome, eight genes were found, and four genes including C3orf67 (the sense strand gene of C3orf67-AS1) were downregulated by Rh2. In particular, C3orf67 was downregulated when C3orf67-AS1 was suppressed by a siRNA; however, the expression of C3orf67-AS1 was not affected by C3orf67. Taken together, this study identifies a novel noncoding RNA, C3orf67-AS1, of which the expression could be suppressed by Rh2 via promoter methylation, thereby mediating the anti-proliferative effect of the ginsenoside.

Classification of early and late stage liver hepatocellular carcinoma patients from their genomics and epigenomics profiles

Background

Liver Hepatocellular Carcinoma (LIHC) is one of the major cancers worldwide, responsible for millions of premature deaths every year. Prediction of clinical staging is vital to implement optimal therapeutic strategy and prognostic prediction in cancer patients. However, to date, no method has been developed for predicting the stage of LIHC from the genomic profile of samples.

Methods

The Cancer Genome Atlas (TCGA) dataset of 173 early stage (stage-I), 177 late stage (stage-II, Stage-III and stage-IV) and 50 adjacent normal tissue samples for 60,483 RNA transcripts and 485,577 methylation CpG sites, was extensively analyzed to identify the key transcriptomic expression and methylation-based features using different feature selection techniques. Further, different classification models were developed based on selected key features to categorize different classes of samples implementing different machine learning algorithms.

Results

In the current study, in silico models have been developed for classifying LIHC patients in the early vs. late stage and cancerous vs. normal samples using RNA expression and DNA methylation data. TCGA datasets were extensively analyzed to identify differentially expressed RNA transcripts and methylated CpG sites that can discriminate early vs. late stages and cancer vs. normal samples of LIHC with high precision. Naive Bayes model developed using 51 features that combine 21 CpG methylation sites and 30 RNA transcripts achieved maximum MCC (Matthew’s correlation coefficient) 0.58 with an accuracy of 78.87% on the validation dataset in discrimination of early and late stage. Additionally, the prediction models developed based on 5 RNA transcripts and 5 CpG sites classify LIHC and normal samples with an accuracy of 96–98% and AUC (Area Under the Receiver Operating Characteristic curve) 0.99. Besides, multiclass models also developed for classifying samples in the normal, early and late stage of cancer and achieved an accuracy of 76.54% and AUC of 0.86.

Conclusion

Our study reveals stage prediction of LIHC samples with high accuracy based on the genomics and epigenomics profiling is a challenging task in comparison to the classification of cancerous and normal samples. Comprehensive analysis, differentially expressed RNA transcripts, methylated CpG sites in LIHC samples and prediction models are available from CancerLSP (http://webs.iiitd.edu.in/raghava/cancerlsp/).

TMEM Proteins in Cancer: A Review

A transmembrane protein (TMEM) is a type of protein that spans biological membranes. Many of them extend through the lipid bilayer of the plasma membrane but others are located to the membrane of organelles. The TMEM family gathers proteins of mostly unknown functions. Many studies showed that TMEM expression can be down- or up-regulated in tumor tissues compared to adjacent healthy tissues. Indeed, some TMEMs such as TMEM48 or TMEM97 are defined as potential prognostic biomarkers for lung cancer. Furthermore, experimental evidence suggests that TMEM proteins can be described as tumor suppressors or oncogenes. TMEMs, such as TMEM45A and TMEM205, have also been implicated in tumor progression and invasion but also in chemoresistance. Thus, a better characterization of these proteins could help to better understand their implication in cancer and to allow the development of improved therapy strategies in the future. This review gives an overview of the implication of TMEM proteins in cancer.

Knockdown of TMEM14A expression by RNAi inhibits the proliferation and invasion of human ovarian cancer cells

Transmembrane protein 14A (TMEM14A) is a member of TMEMs. Alterations in TMEMs expression have been identified in several types of cancer, but the expression and function of TMEM14A in ovarian cancer is still unclear. Here, analysis on the expression data of the Cancer Genome Atlas (TCGA) ovarian serous cystadenocarcinoma (OV) dataset demonstrated the overexpression of TMEM14A in ovarian cancer tissues compared with normal tissues, which was consistent with our real-time PCR analysis on ovarian cancer and normal tissues collected from 30 patients. In addition, TMEM14A knockdown in two ovarian cancer cell lines, A2780 and HO-8910, reduced cell proliferation, causes cell cycle arrest and suppressed cell invasion. Moreover, silencing of TMEM14A notably repressed G1/S cell cycle transition and cell invasion via down-regulating the expression of cell cycle related proteins (Cyclin D1, Cyclin E and PCNA) and metastasis-related proteins (MMP-2 and MMP-9), respectively. TMEM14A knockdown significantly reduced the phosphorylation status of Smad2 and Smad3, downstream effectors of TGF-β signalling. In summary, these results indicate that TMEM14A has a pro-tumorigenic effect in ovarian cancer cells, suggesting an important role of this protein in ovarian cancer oncogenesis and metastasis.

Expression of an estrogen-regulated variant transcript of the peroxisomal branched chain fatty acid oxidase ACOX2 in breast carcinomas

Background

Alternate transcripts from a single gene locus greatly enhance the combinatorial flexibility of the human transcriptome. Different patterns of exon usage have been observed when comparing normal tissue to cancers, suggesting that variant transcripts may play a role in the tumor phenotype.

Methods

Ribonucleic acid-sequencing (RNA-seq) data from breast cancer samples was used to identify an intronic start variant transcript of Acyl-CoA oxidase 2, ACOX2 (ACOX2-i9). Difference in expression between Estrogen Receptor (ER) positive and ER negative patients was assessed by the Wilcoxon rank sum test, and the findings validated in The Cancer Genome Atlas (TCGA) breast cancer dataset (BRCA). ACOX2-i9 expression was also assessed in cell lines using both quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) and Western blot analysis. Knock down by short hairpin RNA (shRNA) and colony formation assays were used to determine whether ACOX2-i9 expression would influence cellular fitness. The effect of ACOX2-i9 expression on patient survival was assessed by the Kaplan-Meier survival function, and association to clinical parameters was analyzed using a Fisher exact test.

Results

The expression and translation of ACOX2-i9 into a 25 kDa protein was demonstrated in HepG2 cells as well as in several breast cancer cell lines. shRNA knock down of the ACOX2-i9 variant resulted in decreased cell viability of T47D and MDA-MB 436 cells. Moreover, expression of ACOX2-i9 was shown to be estrogen regulated, being induced by propyl pyrazoletriol and inhibited by tamoxifen and fulvestrant in ER+ T47D and Mcf-7 cells, but not in the ER- MDA-MB 436 cell line. This variant transcript showed expression predominantly in ER-positive breast tumors as assessed in our initial set of 53 breast cancers and further validated in 87 tumor/normal pairs from the TCGA breast cancer dataset, and expression was associated with better outcome in ER positive patients.

Conclusions

ACOX2-i9 is specifically enriched in ER+ breast cancers where expression of the variant is associated with improved outcome. These data identify variant ACOX2 as a potential novel therapeutic biomarker in ER+ breast tumors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1510-8) contains supplementary material, which is available to authorized users.

Race-associated biological differences among Luminal A breast tumors

African-American (AA) women have higher breast cancer-specific mortality rates. A higher prevalence of the worse outcome Basal-like breast cancer subtype contributes to this, but AA women also have higher mortality even within the more favorable outcome Luminal A breast cancers. These differences may reflect treatment or health care access issues, inherent biological differences, or both. To identify potential biological differences by race among Luminal A breast cancers, gene expression data from 108 CAU and 57 AA breast tumors were analyzed. Race-associated genes were evaluated for associations with survival. Finally, expression of race- and survival-associated genes was evaluated in normal tissue of AA and CAU women. Six genes (ACOX2, MUC1, CRYBB2, PSPH, SQLE, TYMS) were differentially expressed by race among Luminal A breast cancers and were associated with survival (HR <0.8, HR >1.25). For all six genes, tumors in AA had higher expression of poor prognosis genes (CRYBB2, PSPH, SQLE, TYMS) and lower expression of good prognosis genes (ACOX2, MUC1). A score based on all six genes predicted survival in a large independent dataset (HR = 1.9 top vs. bottom quartile, 95% CI: 1.4-2.5). For four genes, normal tissue of AA and CAU women showed similar expression (ACOX2, MUC1, SQLE, TYMS); however, the poor outcome-associated genes CRYBB2 and PSPH were more highly expressed in AA versus CAU women's normal tissue. This analysis identified gene expression differences that may contribute to mortality disparities and suggests that among Luminal A breast tumors there are biological differences between AA and CAU patients. Some of these differences (CRYBB2 and PSPH) may exist from the earliest stages of tumor development, or may even precede malignancy.

Identification of novel targets for antiangiogenic therapy by comparing the gene expressions of tumor and normal endothelial cells

Targeting tumor angiogenesis is an established strategy for cancer therapy. Because angiogenesis is not limited to pathological conditions such as cancer, molecular markers that can distinguish between physiological and pathological angiogenesis are required to develop more effective and safer approaches for cancer treatment. To identify such molecules, we determined the gene expression profiles of murine tumor endothelial cells (mTEC) and murine normal endothelial cells using DNA microarray analysis followed by quantitative reverse transcription–polymerase chain reaction analysis. We identified 131 genes that were differentially upregulated in mTEC. Functional analysis using siRNA-mediated gene silencing revealed five novel tumor endothelial cell markers that were involved in the proliferation or migration of mTEC. The expression of DEF6 and TMEM176B was upregulated in tumor vessels of human renal cell carcinoma specimens, suggesting that they are potential targets for antiangiogenic intervention for renal cell carcinoma. Comparative gene expression analysis revealed molecular differences between tumor endothelial cells and normal endothelial cells and identified novel tumor endothelial cell markers that may be exploited to target tumor angiogenesis for cancer treatment.

Abnormal accumulation of human transmembrane (TMEM)-176A and 176B proteins is associated with cancer pathology

Transmembrane (TMEM)-176A and 176B proteins belong to the MS4A family of proteins whose function in the immune system remains unclear. TMEM176A transcripts were previously shown to be elevated in liver cancer or kidney tissue with proteinuria, while marked changes in TMEM176B transcripts have been found in tolerated tissue allografts and neoplastic fibroblasts. To study the functional relationship between human TMEM176A and 176B and their putative link to cancer, we used polymerase chain reaction and biochemical assays. Here, we show that TMEM176A and 176B are widely expressed in all human tissues examined. Co-immunoprecipitation of heterologously expressed TMEM176A and 176B revealed direct physical interaction. To determine the relevance of such interaction to cancer pathology, we analyzed biopsied tissue samples from a variety of normal and cancer tissues. Our data reveal that human TMEM176A and 176B protein levels are significantly elevated in lymphoma, but not in normal tissues. The protein levels of TMEM176A are also significantly increased in lung carcinoma. Finally, analysis of the protein expression ratio of TMEM176A over 176B showed significant differences between normal and cancer tissues of the breast, lymph, skin, and liver, which indicates that both TMEM proteins could be potential useful markers for certain human cancers.

Molecular characteristics of non-cancerous liver tissue in non-B non-C hepatocellular carcinoma

Although chronic infection with hepatitis B virus (HBV) and/or hepatitis C virus (HCV) are the most important risk factors for the development of hepatocellular carcinoma (HCC) worldwide, the proportion of HCC patients negative for the hepatitis B surface antigen and hepatitis C antibody, so-called "non-B non-C HCC", is rapidly increasing, especially in Japan. The background liver diseases of non-B non-C HCC patients can be multifactorial, including occult HBV infection and non-alcoholic steatohepatitis. It is reasonable to investigate the non-cancerous liver tissues to identify the potential molecular mechanisms responsible for the processes of hepatocarcinogenesis of non-B non-C HCC. However, to date, only a few studies have focused on this research concept based on the idea of "field cancerization". This review highlights the potential importance of the molecular analysis of non-cancerous liver tissues to clarify the molecular characteristics in patients with non-B non-C HCC. A better understanding of the molecular mechanisms underlying the individual predisposition to non-B non-C HCC will lead to improvements in the prevention, early diagnosis and treatment of this neoplastic disease.

Progression-free survival in ovarian cancer is reflected in epigenetic DNA methylation profiles

OBJECTIVE

Many patients with ovarian cancer disease relapse within 6 months after adjuvant chemotherapy, with a limited prognosis. Epigenetic modifications have been shown to play an important role in tumor development and formation. Therefore, global analysis of DNA methylation patterns might reveal specific CpG sites that correlate with progression-free interval (PFI) after therapy.

METHODS

Twenty samples of advanced ovarian cancer with a predominantly serous papillary histological subtype were subjected to DNA methylation profiling. Illumina HumanMethylation27 BeadChip technology was used for simultaneous analysis of 27,578 CpG sites in >14,000 genes.

RESULTS

Differential DNA methylation of various cytosines correlated with PFI. However, this becomes only significant by classification according to PFI with a cutoff of >28 months. Longer survival was associated with hypomethylation at specific CpG sites (e.g. GREB1, TGIF and TOB1) and hypermethylation in other genes (e.g. TMCO5, PTPRN and GUCY2C). Gene ontology analysis revealed that differentially methylated genes were significantly overrepresented in the categories telomere organization, mesoderm development and immune regulation.

CONCLUSION

Epigenetic modifications at specific CpG sites correlate with PFI in ovarian cancer. Therefore, such analysis might be of prognostic value.

Identification of a secretory protein c19orf10 activated in hepatocellular carcinoma

The identification of genes involved in tumor growth is crucial for the development of inventive anticancer treatments. Here, we have cloned a 17-kDa secretory protein encoded by c19orf10 from hepatocellular carcinoma (HCC) serial analysis of gene expression libraries. Gene expression analysis indicated that c19orf10 was overexpressed in approximately two-thirds of HCC tissues compared to the adjacent noncancerous liver tissues, and its expression was significantly positively correlated with that of alpha-fetoprotein (AFP). Overexpression of c19orf10 enhanced cell proliferation of AFP-negative HLE cells, whereas knockdown of c19orf10 inhibited cell proliferation of AFP-positive Hep3B and HuH7 cells along with G1 cell cycle arrest. Supplementation of recombinant c19orf10 protein in culture media enhanced cell proliferation in HLE cells, and this effect was abolished by the addition of antibodies developed against c19orf10. Intriguingly, c19orf10 could regulate cell proliferation through the activation of Akt/mitogen-activated protein kinase pathways. Taken together, these data suggest that c19orf10 might be one of the growth factors and potential molecular targets activated in HCC.