Microarray comparative genomic hybridization reveals genome-wide patterns of DNA gains and losses in post-Chernobyl thyroid cancer

OUP accepted manuscript

This study aimed to investigate the molecular mechanisms through which the intestinal microbiota and microRNAs (miRNAs) participate in colon cancer metastasis. Intestinal flora data, and the GSE29621 (messenger RNA/long non-coding RNA [mRNA/lncRNA]) and GSE29622 (miRNA) datasets, were downloaded from The Cancer Gene Atlas and Gene Expression Omnibus databases, respectively. Immune-related cells in M1 vs. M0 samples were analyzed using the Wilcoxon test. Furthermore, an lncRNA-miRNA-mRNA (competing endogenous RNA [ceRNA]) network was constructed, and survival analysis of RNAs in the network was performed. A total of 16 miRNA-genus co-expression pairs containing eight microbial genera and 15 miRNAs were screened; notably, Porphyromonas and Bifidobacterium spp. were found to be associated with most miRNAs, and has-miR-3943 was targeted by most microbial genera. Furthermore, five immune cell types, including activated natural killer cells, M1 macrophages, resting mast cells, activated mast cells and neutrophils, were differentially accumulated between the M1 and M0 groups. Enrichment analysis suggested that mRNAs related to colon cancer metastasis were mainly involved in pathways related to bacterial and immune responses. Survival analysis revealed that TMEM176A and PALM3 in the ceRNA network were significantly associated with the prognosis of patients with colon cancer. In conclusion, this study revealed a potential mechanism by which the intestinal microbiota influences the colon cancer microenvironment by targeting miRNAs.

Methylation of TMEM176A, a key ERK signaling regulator, is a novel synthetic lethality marker of ATM inhibitors in human lung cancer

Aim: The role of TMEM176A methylation in lung cancer and its therapeutic application remains unclear. Materials and methods: Nine lung cancer cell lines and 123 cases of cancer tissue samples were employed. Results: TMEM176A was methylated in 53.66% of primary lung cancer. Restoration of TMEM176A expression induced cell apoptosis and G2/M phase arrest, and inhibited colony formation, cell proliferation, migration and invasion. TMEM176A suppressed H1299 cell xenograft growth in mice. Methylation of TMEM176A activated ERK signaling and sensitized H1299 and H23 cells to AZD0156, an ATM inhibitor. Conclusion: The expression of TMEM176A is regulated by promoter region methylation. Methylation of TMEM176A is a potential lung cancer diagnostic marker and a novel synthetic lethal therapeutic marker for AZD0156.

Genomic landscape of metastatic papillary thyroid carcinoma and novel biomarkers for predicting distant metastasis

Papillary thyroid carcinoma (PTC) is the most common malignancy of the thyroid gland, with a relatively high cure rate. Distant metastasis (DM) of PTC is uncommon, but when it occurs, it significantly decreases the survival of PTC patients. The molecular mechanisms of DM in PTC have not been systematically studied. We performed whole exome sequencing and GeneseeqPrime (425 genes) panel sequencing of the primary tumor, plasma and matched white blood cell samples from 20 PTC with DM and 46 PTC without DM. We identified somatic mutations, gene fusions and copy number alterations and analyzed their relationships with DM of PTC. BRAF-V600E was identified in 73% of PTC, followed by RET fusions (14%) in a mutually exclusive manner (P < 0.0001). We found that gene fusions (RET, ALK or NTRK1) (P < 0.01) and chromosome 22q loss (P < 0.01) were independently associated with DM in both univariate and multivariate analyses. A nomogram model consisting of chromosome 22q loss, gene fusions and three clinical variables was built for predicting DM in PTC (C-index = 0.89). The plasma circulating tumor DNA (ctDNA) detection rate in PTC was only 38.9%; however, it was significantly associated with the metastatic status (P = 0.04), tumor size (P = 0.001) and invasiveness (P = 0.01). In conclusion, gene fusions and chromosome 22q loss were independently associated with DM in PTC and could serve as molecular biomarkers for predicting DM. The ctDNA detection rate was low in non-DM PTC but significantly higher in PTC with DM.

DNA Copy Number Variations as Markers of Mutagenic Impact

DNA copy number variation (CNV) occurs due to deletion or duplication of DNA segments resulting in a different number of copies of a specific DNA-stretch on homologous chromosomes. Implications of CNVs in evolution and development of different diseases have been demonstrated although contribution of environmental factors, such as mutagens, in the origin of CNVs, is poorly understood. In this review, we summarize current knowledge about mutagen-induced CNVs in human, animal and plant cells. Differences in CNV frequencies induced by radiation and chemical mutagens, distribution of CNVs in the genome, as well as adaptive effects in plants, are discussed. Currently available information concerning impact of mutagens in induction of CNVs in germ cells is presented. Moreover, the potential of CNVs as a new endpoint in mutagenicity test-systems is discussed.

Epigenetic silencing of TMEM176A activates ERK signaling in human hepatocellular carcinoma

Background

The role of TMEM176A in human hepatocellular carcinoma (HCC) is unknown. This study explored the epigenetic regulation and function of TMEM176A in human HCC.

Materials and methods

Twelve HCC cell lines and 126 cases of primary cancer were analyzed. Methylation-specific PCR, immunohistochemistry, flow cytometry, and xenograft mouse models were employed.

Results

TMEM176A was highly expressed in SNU387, SNU182, Huh1, and SNU475 cells; reduced expression was observed in HepG2 and PLC/PRF/5 cells; and no expression was found in SNU449, HBXF344, SMMC7721, Huh7, and LM3 cells. Unmethylation of the TMEM176A promoter was detected in SNU387, SNU182, Huh1, and SNU475 cells; partial methylation was observed in HepG2 and PLC/PRF/5 cells; and complete methylation was found in SNU449, HBXF344, SMMC7721, Huh7, and LM3 cells. Upon treatment with 5-Aza-2-deoxycytidine, re-expression of TMEM176A was detected in SNU449, HBXF344, SMMC7721, Huh7, and LM3 cells; increased expression of TMEM176A was observed in HepG2 and PLC/PRF/5 cells; and no expression changes were found in SNU387, SNU182, Huh1, and SNU475 cells. The TMEM176A promoter region was methylated in 75.4% (95/126) of primary human HCC. Reduced expression of TMEM176A was associated with promoter region methylation (P < 0.05). No association was found between TMEM176A promoter methylation and age, gender, HBV infection, liver cirrhosis, tumor size, lymph node metastasis, vessel cancerous embolus, number of lesions, and TNM stage (all P > 0.05). These results demonstrated that the expression of TMEM176A is regulated by promoter region methylation. Methylation of the TMEM176A promoter was significantly associated with tumor cell differentiation (P < 0.05) and was an independent prognostic factor for poor 3-year overall survival (OS, P < 0.05). TMEM176A expression induced cell apoptosis; inhibited cell proliferation, migration, and invasion; suppressed human HCC cell xenograft growth in mice; and inhibited ERK signaling in HCC cells.

Conclusion

The promoter region of TMEM176A is frequently methylated in human HCC, and the expression of TMEM176A is regulated by promoter region methylation. Methylation of the TMEM176A promoter may serve as a diagnostic and prognostic marker in HCC. TMEM176A suppresses HCC growth by inhibiting the ERK signaling pathway.

Epigenetic silencing of TMEM176A promotes esophageal squamous cell cancer development

The function of human transmembrane protein 176A (TMEM176A) in cancer remains unclear. To understand the function and mechanism of TMEM176A in human esophageal cancer development, 13 esophageal cancer cell lines and 267 cases of primary esophageal squamous cell cancer (ESCC) samples were analyzed by methylation specific PCR (MSP), flow cytometry, immunohistochemistry and transfection assays. TMEM176A was highly expressed in BIC1 cells and loss of TMEM176A expression was found in TE1, TE3, TE13, KYSE140, KYSE180, KYSE410, KYSE450, KYSE520, Segl, KYSE150, YES2 and COLO680N cells. Complete methylation was detected in TE1, TE3, TE13, KYSE140, KYSE180, KYSE410, KYSE450, KYSE520, Segl, KYSE150, YES2 and COLO680N cells, while unmethylation was detected in BIC1 cells. Restoration of TMEM176A expression was induced by 5-aza-2’-deoxycytidine treatment in methylated cell lines. TMEM176A was methylated in 66.7% (178/267) of primary esophageal cancer samples, and promoter region methylation was significantly associated with tumor differentiation (p<0.001) and loss off/reduced expression of TMEM176A (p<0.05). Methylation of TMEM176A was significantly associated with poor 5-year overall survival (p < 0.05). Cox proportional hazards model analysis suggest that TMEM176A methylation is an independent prognostic factor for poor 5-years OS. TMEM176A inhibited cell invasion and migration, and induced apoptosis in esophageal cancer cells. TMEM176A suppressed esophageal cancer cell growth both in vitro and in vivo. In conclusion, TMEM176A is frequently methylated in human ESCC and the expression of TMEM176A is regulated by promoter region methylation. TMEM176A methylation may serve as a diagnostic and prognostic marker in ESCC. TMEM176A is a potential tumor suppressor in human ESCC.

Methylation of TMEM176A is an independent prognostic marker and is involved in human colorectal cancer development

Colorectal cancer (CRC) is the third most common malignancy and the fourth most common cause of cancer related death worldwide. This study was designed to find tumor suppressors involved in CRC development by performing RNA-seq. Eight CRC cell lines and 130 cases of primary CRC samples were used. RNA-seq, methylation-specific PCR (MSP), flow cytometry, transwell assays, and a xenograft mouse model were used. Reduction of TMEM176A expression was confirmed in human CRC cells by RNA-seq. TMEM176A was expressed in LS180 and SW620 cells, loss of TMEM176A expression was observed in LOVO, HCT116, RKO, and DLD1 cells, and reduced TMEM176A expression was found in HT29 and SW480 cells. Unmethylation of the TMEM176A promoter was found in LS180 and SW620 cells, whereas complete methylation was found in LOVO, HCT116, RKO, and DLD1 cells, and partial methylation was found in HT29 and SW480 cells. Promoter region methylation correlated with loss of/reduced expression of TMEM176A. Re-expression of TMEM176A was induced by 5-aza-2'-deoxycytidine. TMEM176A was methylated in 50.77% of primary colorectal cancers. Methylation of TMEM176A was associated with tumor metastasis (P<0.05) and was an independent prognostic factor for 5-year overall survival (OS) according to Cox proportional hazards model analysis (P<0.05). TMEM176A induced apoptosis and inhibited cell migration and invasion in CRC cells. TMEM176A suppressed CRC cell growth both in vitro and in vivo. Our results suggest that expression of TMEM176A is regulated by promoter region methylation. TMEM176A methylation is an independent prognostic marker for 5-year OS in CRC, and may act as a tumor suppressor in CRC.

Genome-wide analysis of abdominal and pleural malignant mesothelioma with DNA arrays reveals both common and distinct regions of copy number alteration

Malignant mesothelioma (MM) is an aggressive tumor arising from mesothelial linings of the serosal cavities. Pleural space is the most common site, accounting for about 80% of cases, while peritoneum makes up the majority of the remaining 20%. While histologically similar, tumors from these sites are epidemiologically and clinically distinct and their attribution to asbestos exposure differs. We compared DNA array-based findings from 48 epithelioid peritoneal MMs and 41 epithelioid pleural MMs to identify similarities and differences in copy number alterations (CNAs). Losses in 3p (BAP1 gene), 9p (CDKN2A) and 22q (NF2) were seen in tumors from both tumor sites, although CDKN2A and NF2 losses were seen at a higher rate in pleural disease (p<0.01). Overall, regions of copy number gain were more common in peritoneal MM, whereas losses were more common in pleural MM, with regions of loss containing known tumor suppressor genes and regions of gain encompassing genes encoding receptor tyrosine kinase pathway members. Cases with known asbestos causation (n = 32 ) were compared with those linked to radiation exposure (n = 9 ). Deletions in 6q, 14q, 17p and 22q, and gain of 17q were seen in asbestos-associated but not radiation-related cases. As reported in post-radiation sarcoma, gains outnumbered losses in radiation-associated MM. The patterns of genomic imbalances suggest overlapping and distinct molecular pathways in MM of the pleura and peritoneum, and that differences in causation (i.e., asbestos vs. radiation) may account for some of these site-dependent differences.

Radiation signatures in childhood thyroid cancers after the Chernobyl accident: possible roles of radiation in carcinogenesis

After the Tokyo Electric Power Company Fukushima Daiichi nuclear power plant accident, cancer risk from low-dose radiation exposure has been deeply concerning. The linear no-threshold model is applied for the purpose of radiation protection, but it is a model based on the concept that ionizing radiation induces stochastic oncogenic alterations in the target cells. As the elucidation of the mechanism of radiation-induced carcinogenesis is indispensable to justify the concept, studies aimed at the determination of molecular changes associated with thyroid cancers among children who suffered effects from the Chernobyl nuclear accident will be overviewed. We intend to discuss whether any radiation signatures are associated with radiation-induced childhood thyroid cancers.

Genomic landscape of pancreatic neuroendocrine tumors

AIM: To investigate the prognostic role of genomic stability and copy number alterations (CNAs) pancreatic neuroendocrine tumors (PanNETs).

METHODS: A high-resolution array-based comparative genomic hybridization approach was utilized in order to investigate and quantify chromosomal aberrations in a panel of 37 primary PanNET and 11 metastatic samples. DNA samples were extracted from formalin-fixed and paraffin-embedded tumor specimen. Genomic findings were correlated with histopathological and immunohistochemical data. Moreover, the dataset was subjected to employing an unsupervised hierarchical clustering analysis approach utilizing Euclidean distance and average linkage and associations between genomically defined tumor groups and recurrent CNAs or clinicopathological features of the study group were assessed.

RESULTS: Numerous chromosomal aberrations were recurrently detected in both, primary tumor samples and metastases. Copy number gains were most frequently observed at 06p22.2-p22.1 (27.1%), 17p13.1 (20.8%), 07p21.3-p21.2 (18.8%), 09q34.11 (18.8%). Genomic losses were significantly less frequent and the only recurrent aberration affected 08q24.3 (6.3%). Moreover, we detected a high degree of genomic heterogeneity between primary tumors and metastatic lesions. Unsupervised hierarchical clustering of loci affected by CNAs in more than 3 primary tumor samples revealed two genetically distinct tumor groups as well as two chromosomal clusters of genomic imbalances indicating a small subset of tumors with common molecular features (13.5%). Aberrations affecting 6p22.2-22.1, 8q24.3, 9q34.11 and 17p13.1 (P = 0.011; 0.003; 0.003; 0.001), were significantly associated with a poorer survival prognosis.

CONCLUSION: This study suggests that several frequent CNAs in numerous candidate regions are involved in the pathogenesis and metastatic progression of PanNET.

The chernobyl tissue bank - a repository for biomaterial and data used in integrative and systems biology modeling the human response to radiation

The only unequivocal radiological effect of the Chernobyl accident on human health is the increase in thyroid cancer in those exposed in childhood or early adolescence. In response to the scientific interest in studying the molecular biology of thyroid cancer post Chernobyl, the Chernobyl Tissue Bank (CTB: www.chernobyltissuebank.com) was established in 1998. Thus far it is has collected biological samples from 3,861 individuals, and provided 27 research projects with 11,254 samples. The CTB was designed from its outset as a resource to promote the integration of research and clinical data to facilitate a systems biology approach to radiation related thyroid cancer. The project has therefore developed as a multidisciplinary collaboration between clinicians, dosimetrists, molecular biologists and bioinformaticians and serves as a paradigm for tissue banking in the omics era.

Discriminating gene expression signature of radiation-induced thyroid tumors after either external exposure or internal contamination

Both external radiation exposure and internal radionuclide contamination are well known risk factors in the development of thyroid epithelial tumors. The identification of specific molecular markers deregulated in radiation-induced thyroid tumors is important for the etiological diagnosis since neither histological features nor genetic alterations can discriminate between sporadic and radiation-induced tumors. Identification of highly discriminating markers in radiation-induced tumors is challenging as it relies on the ability to identify marker deregulation which is associated with a cellular stress that occurred many years before in the thyroid cells. The existence of such a signature is still controversial, as it was not found in several studies while a highly discriminating signature was found in both post-radiotherapy and post-Chernobyl series in other studies. Overall, published studies searching for radiation-induced thyroid tumor specificities, using transcriptomic, proteomic and comparative genomic hybridization approaches, and bearing in mind the analytical constraints required to analyze such small series of tumors, suggest that such a molecular signature could be found. In comparison with sporadic tumors, we highlight molecular similarities and specificities in tumors occurring after high-dose external radiation exposure, such as radiotherapy, and in post-Chernobyl tumors that occurred after internal 131I contamination. We discuss the relevance of signature extrapolation from series of tumors developing after high and low doses in the identification of tumors induced at very low doses of radiation.

A transcriptome signature distinguished sporadic from postradiotherapy radiation-induced sarcomas

Exposure to ionizing radiation is a known risk factor for cancer. However, up to now, rigorously defined scientific criteria that could establish case-by-case the radiation-induced (RI) origin of a tumour have been lacking. To identify genes that could constitute a RI signature, we compared the transcriptome of 12 sarcomas arising in the irradiation field of a primary tumour following radiotherapy with the transcriptome of 12 sporadic sarcomas. This learning/training set contained four leiomyosarcomas, four osteosarcomas and four angiosarcomas in each subgroup. We identified a signature of 135 genes discriminating RI from sporadic sarcomas. The robustness of this signature was tested by the blind case-by-case classification of an independent set of 36 sarcomas of various histologies. Thirty-one sarcomas were classified as RI or sporadic; it was not possible to propose an aetiology for the five others. After the code break, it was found that one sporadic sarcoma was misclassified as RI. Thus, the signature is robust with a sensitivity of 96%, a positive and a negative predictive value of 96 and 100%, respectively and a specificity of 62%. The functions of the genes of the signature suggest that RI sarcomas were subject to chronic oxidative stress probably due to mitochondrial dysfunction.

Genomic and gene expression signatures of radiation in medulloblastomas after low-dose irradiation in Ptch1 heterozygous mice

Accurate cancer risk assessment of low-dose radiation poses many challenges that are partly due to the inability to distinguish radiation-induced tumors from spontaneous ones. To elucidate characteristic features of radiation-induced tumors, we analyzed 163 medulloblastomas that developed either spontaneously or after X-ray irradiation at doses of 0.05-3 Gy in Ptch1 heterozygous mice. All spontaneous tumors showed loss of heterozygosity in broad regions on chromosome 13, with losses at all consecutive markers distal to Ptch1 locus (S-type). In contrast, all tumors that developed after 3 Gy irradiation exhibited interstitial losses around Ptch1 with distal markers retained (R-type). There was a clear dose-dependent increase in the proportion of R-type tumors within the intermediate dose range, indicating that the R-type change is a reliable radiation signature. Importantly, the incidence of R-type tumors increased significantly (P = 0.007) at a dose as low as 50 mGy. Integrated array-comparative genomic hybridization and expression microarray analyses demonstrated that expression levels of many genes around the Ptch1 locus faithfully reflected the signature-associated reduction in genomic copy number. Furthermore, 573 genes on other chromosomes were also expressed differently between S-type and R-type tumors. They include genes whose expression changes during early cerebellar development such as Plagl1 and Tgfb2, suggesting a recapitulation of gene subsets functioning at distinct developmental stages. These findings provide, for the first time, solid experimental evidence for a significant increase in cancer risk by low-dose radiation at diagnostic levels and imply that radiation-induced carcinogenesis accompanies both genomic and gene expression signatures.

Downregulation of ZIP kinase is associated with tumor invasion, metastasis and poor prognosis in gastric cancer

Deletion of 19p13 is one of the most frequent genetic changes in gastric carcinoma (GC), implying the existence of a tumor suppressor gene (TSG) that plays an important role in GC development. To identify the candidate TSG at 19p, array-comparative genomic hybridization (CGH) was applied to study DNA copy-number changes on chromosomes 3, 5p, 13, 16q and 19. The result showed that gains of 16q21, 19q13.1, 5p15.1 and 3q26.31, and losses of 3p21.32, 3p22.2, 19q13.33 and 19p13.3, were frequently detected by array-CGH. One candidate TSG, ZIP kinase (ZIPK), at 19p13.3 was further characterized by immunohistochemistry using a tissue microarray containing 172 primary GCs. Downregulation of ZIPK was detected in 111/162 informative GCs, which was significantly associated with invasion, metastasis and poorer prognosis of GC. To investigate the association of the downregulation of ZIPK with apoptosis, apoptosis assay (TUNEL) was used to compare the apoptotic index between GCs with normal expression and downregulation of ZIPK. TUNEL assay showed that the apoptotic index in GCs with normal ZIPK expression was significantly higher than that in GCs with downregulation of ZIPK (p < 0.001), indicating that ZIPK plays an important pro-apoptotic role in GC. Taken together, we demonstrated here that ZIPK is a tumor suppresser gene and plays an important role in GC development through its pro-apoptotic function. Downregulation of ZIPK can be used to evaluate tumor invasiveness, metastasis and to predict survival of GC.

Genome-wide studies in thyroid neoplasia

There is much interest in the application of genome biology to the field of thyroid neoplasia, despite the relatively low mortality rate associated with thyroid cancer in general. The principal reason for this interest is that the field of thyroid neoplasia stands to benefit from the application of genomic information to address a variety of pathologic and clinical issues. In addition to practical patient care issues, there is an excellent opportunity of expand the basic understanding of thyroid carcinogenesis. In this article, the most relevant genomic work on thyroid tumors performed to date is reviewed along with some general comments about the potential impact of genomic biology on thyroid pathology and the management of patients with thyroid nodules and cancer.

DNA copy-number instability in low-dose gamma-irradiated TK6 lymphoblastoid clones

Genomic instability that might occur early during low-dose, fractionated radiation exposures may be traceable in radiogenic compared to spontaneous cancers. Using a human 18K cDNA microarray-based comparative genome hybridization protocol, we measured changes in DNA copy number at over 14,000 loci in nine low-dose (137)Cs gamma-irradiated (acute exposure to 10 cGy/day x 21 days) and nine unirradiated TK6 clones and estimated locus-specific copy-number differences between them. Radiation induced copy-number hypervariability at thousands of loci across all chromosomes, with a sevenfold increase in low-level, randomly positioned DNA gains. Recurrent gains at 40 loci occurred among irradiated clones and were distributed nonrandomly across the genome, with the highest densities in 3q, 13q and 20q at sites that were hypodiploid without irradiation. Another nonrandomly distributed set of 94 loci exhibited relative recurrent gains from a hypodiploid state to a diploid state, suggesting hemizygous-to-homozygous transitions. Frequently recurring losses at 57 loci were concentrated on the single X-chromosome but were sparsely distributed at 0-2 loci per autosome. These results suggest induced mitotic homologous recombination as a possible mechanism of low-dose radiation-induced genomic instability. Genomic instability induced in TK6 cells resembled that seen in radiogenic tumors and suggests a way that radiation could induce genomic instability in preneoplastic cells.

A radiation-induced gene signature distinguishes post-Chernobyl from sporadic papillary thyroid cancers

We investigated selected gene targets to differentiate radiation-induced papillary thyroid cancers (PTCs) from other etiologies. Total RNA was isolated from 11 post-Chernobyl PTCs and 41 sporadic PTCs characterized by a more aggressive tumor type and lacking a radiation exposure history. RNA from 10 tumor samples from both groups was pooled and hybridized separately on a whole genome microarray for screening. Then 92 selected gene targets were examined quantitatively on each tumor sample using an RTQ-PCR-based low-density array (LDA). Screening for more than fivefold differences in gene expression between the groups by microarray detected 646 up-regulated and 677 down-regulated genes. Categorization of these genes revealed a significant (P < 0.0006) over-representation of the number of up-regulated genes coding for oxidoreductases, G-proteins and growth factors, while the number of genes coding for immunoglobulin appeared to be significantly down-regulated. With the LDA, seven genes (SFRP1, MMP1, ESM1, KRTAP2-1, COL13A1, BAALC and PAGE1) made a complete differentiation between the groups possible. Gene expression patterns known to be associated with a more aggressive tumor type in older patients appeared to be more pronounced in post-Chernobyl PTC, thus underlining the known aggressiveness of radiation-induced PTC. Seven genes were found that completely distinguished post-Chernobyl (PTC) from sporadic PTC.