Search

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. The molecular characteristics of histologically normal appearing tissue adjacent to the tumor (NAT) from PTC patients are not well characterized. The aim of this study was to characterize the global gene expression profile of NAT and compare it with those of normal and tumor thyroid tissues. We performed total RNA sequencing with fresh frozen thyroid tissues from a cohort of three categories of samples including NAT, normal thyroid (N), and PTC tumor (T). Transcriptome analysis shows that NAT presents a unique gene expression profile, which was not associated with sex or the presence of lymphocytic thyroiditis. Among the differentially expressed genes (DEGs) of NAT vs N, 256 coding genes and 5 noncoding genes have been reported as cancer genes involved in cell proliferation, apoptosis, and/or tumorigenesis. Bioinformatics analysis with Ingenuity Pathway Analysis software revealed that “Cancer, Organismal Injury and Abnormalities, Cellular Response to Therapeutics, and Cellular Movement” were major dysregulated pathways in the NAT tissues. This study provides improved insight into the complexity of gene expression changes in the thyroid glands of patients with PTC.

BackgroundPapillary thyroid carcinoma (PTC) demonstrates high heritability and a low somatic mutation burden relative to other cancers. Therefore, the genetic risk predisposing to PTC is likely due to a combination of low penetrance variants. A recent genome-wide association study revealed the association of PTC with a missense variant, rs6793295, at 3q26 in a gene called Leucine Repeat Rich Containing 34 (LRRC34).MethodsWe report the mechanisms of PTC risk at 3q26 using a combination of overexpression, mass spectroscopy, knockdown, transcriptome profiling, migration assays and genetic analysis.ResultsWe observed differential binding of wild-type and missense LRRC34 to RANBP1. Overexpression of missense LRRC34 reduced RanGTP levels and increased apoptosis. We also identified a second linkage disequilibrium (LD) block upstream of LRRC34 containing regulatory variants with allele-specific expression. Transcriptome profiling of LRRC34 knockdown cells showed changes in genes involved with cellular movement. LRRC34 knockdown reduced the migration of thyroid cancer cell lines. Lastly, we assessed the relative contribution of PTC risk from each locus using haplotype analysis.ConclusionsOur study demonstrates two separate mechanisms, one in G protein signalling and the other in transcriptional control, dictating PTC risk at 3q26 using both biochemical and genetic techniques.

Background: Our genome sequencing analysis revealed a frameshift mutation in the shelterin gene TINF2 in a large family with individuals affected with papillary thyroid carcinoma (PTC) and melanoma. Here, we further characterized the mutation and screened for coding variants in the 6 shelterin genes in 24 families. Methods: Sanger sequencing was performed to screen for the TINF2 mutation in the key family. Quantitative reverse transcription-polymerase chain reaction (PCR) was used for TINF2 gene expression analysis. Exogenous expression and co-immunoprecipitation techniques were used for assessing TINF2 binding to TERF1. Relative telomere length (RTL) was quantified in DNAs from lymphocytes by using quantitative real-time PCR. Whole exome sequencing (WES) was performed in seven families with individuals affected with PTC and other cancer types. Screening for DNA variants in shelterin genes was performed by using whole genome sequencing data from 17 families and WES data from 7 further families. Results: The TINF2 mutation (TINF2 p.Trp198fs) showed complete co-segregation with PTC and melanoma in the key family. The mutation is not reported in databases and not identified in 23 other families we screened. The expression of TINF2 was borderline reduced in individuals with the mutation. The truncated TINF2 protein showed abolished binding to TERF1. The RTL in the individuals with the mutation was significantly longer when compared with those without the mutation from the same family as well as compared with 62 healthy controls. Among the 24 families, we identified 3 missense and 1 synonymous variant(s) in 2 shelterin genes (TINF2 and ACD). Conclusions: The rare frameshift mutation in the TINF2 gene and the associated longer telomere length suggest that dysregulated telomeres could be a mechanism predisposing to PTC and melanoma. DNA coding variants in shelterin genes are rare. Further studies are required to evaluate the roles of variants in shelterin genes in thyroid cancer and melanoma.

Background We studied a large family with 22 individuals affected with autosomal dominant hereditary spherocytosis (HS). Methods Genome‐wide linkage, whole‐genome sequencing (WGS), Sanger sequencing, RT‐PCR, and ToPO TA cloning analyses were performed. Results We revealed a heterozygous G>A transition in the 14q23 locus, at position +1 of the intron 8 donor splice site of the spectrin beta, erythrocytic (SPTB) gene. This splice variant (SPTB c.1064+1G>A) was confirmed by Sanger sequencing and showed complete co‐segregation with HS in the family. Further RT‐PCR reactions and sequencing analysis indicated that the variant leads to the exclusion of exon 8 and subsequent frameshift in exon 9 and a premature stop codon in SPTB. Translation of the altered allele would lead to a truncation with a loss of all spectrin repeat domains in SPTB protein. Conclusion This variant is novel and has not been found in any databases. We propose that this splice variant explains the spherocytosis phenotype observed in this large family., We studied a large family with 22 individuals affected with autosomal dominant hereditary spherocytosis (HS). Genome‐wide linkage and whole‐genome sequencing analyses revealed a heterozygous G>A transition in the 14q23 locus, at position +1 of the intron 8 donor splice site of the spectrin beta, erythrocytic (SPTB) gene.

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. Here, we conducted a genome-wide association study (GWAS) involving 1,554 cases/1,973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from the EPITHYR consortium. Our results confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 (rs334729) and 16q23.2 (rs16950982), which were associated with TSH levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. The frequency of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans and may explain part of the highest DTC incidence observed in Oceanians.Competing Interest StatementThe authors have declared no competing interest.Funding StatementINCA (grant number 9533) and ARC (grant number PGA120150202302), Ligue Nationale Contre le Cancer, Agence Nationale pour la Recherche (ANR), the Direction Generale de la Sante, the Agence Francaise de Securite Sanitaire de l alimentation, de lenvironnement et du travail (ANSES), CHILDTHYR EEC program, and the Fondation de France. JG and CX were the recipient of a PhD fellowship from Region Ile-de-France.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:All participants provided informed consent and each study was approved by their governing ethics committee.All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).Yes I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesThe data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions. (Less)

Papillary thyroid carcinoma (PTC) is the most common histotype of thyroid carcinoma. The heritability of PTC is high compared to other cancers, but its underlying causes are unknown. A recent genome-wide association study revealed the association of a variant at the 5q22 locus, rs73227498, with PTC predisposition. We report that rs17134155, a variant in high linkage disequilibrium with rs73227498, is located in an enhancer region downstream of coding transcripts of EPB41L4A. Rs17134155 showed significant enhancer activity in luciferase assays, and haplotypes containing the protective allele of this variant conferred a significantly lower risk of PTC. While the index SNP, rs73227498, acted as a significant cis-eQTL for expression of EPB41L4A, rs17134155 was a significant cis-sQTL for the alternative splicing of a non-coding transcript of EPB41L4A, called EPB41L4A-203. We also performed knockdown of EPB41L4A followed by microarray analysis. Some of the top differentially-expressed genes were represented among regulators of the WNT/β-catenin signaling pathway. Our results indicate that an enhancer region at 5q22 regulates the expression and splicing of EPB41L4A transcripts. We also provide evidence that EPB41L4A expression is involved in regulating growth and differentiation pathways, suggesting that decreased expression of EPB41L4A is a mechanism in the predisposition to PTC.

The first two genome wide association studies (GWAS) of papillary thyroid carcinoma (PTC) detected five variants associated with PTC. Two of them (rs944289 and rs116909374) are located at 14q13 making that locus an important target of research into the genetic predisposition to PTC. We aimed at uncovering other variants at 14q13 associated with PTC independently from the GWAS variants. We performed next generation sequencing of the 14q13 region and analyzed the allele frequencies of single nucleotide polymorphisms (SNPs) in n = 90 PTC cases vs. n = 379 EUR controls from the 1,000 Genome Project. The variants associated with PTC were validated in an Ohio cohort of n = 1,216 PTC cases and n = 1,416 controls. Next, we analyzed the association between SNPs and expression of nearby genes and clinical parameters. We showed that rs368187 was associated with PTC (OR = 1.31, p = 2.20 × 10-6 ). Rs1632250, Rs1863347 and rs1755787 showed association with classical PTC (cPTC) (n = 891; OR = 1.24, 2.22 × 10-3 , OR = 1.31, p = 2.15 × 10-4 and OR = 1.24, p = 2.06 × 10-3 , respectively) while variant rs28397092 showed association with follicular variant (n = 243; OR = 1.51, p = 1.36 × 10-3 ). Rs1863347 was associated with suppression of PTCSC3 in unaffected thyroid tissue (p = 0.026). Rs1632250, rs1863347 and rs1755787 showed association with multifocality (OR = 1.85, p = 0.001, OR = 1.98, p = 0.001 and OR = 1.76, p = 0.003 respectively) and N stage (OR = 1.79, p = 0.014, OR = 1.73, p = 0.023 and OR = 1.81, p = 0.013, respectively) in microPTC (n = 328) while rs368187 was associated with M stage (OR = 0.56, p = 0.034) in cPTC. Our results disclose multiple variants associated with PTC and clinical features in the 14q13 superlocus. We suggest that translational genotype/phenotype studies should take into account not only somatic mutations but also germline variants.

To identify and characterize the functional variants, regulatory gene networks, and potential binding targets of SMAD3 in the 15q22 thyroid cancer risk locus. We performed linkage disequilibrium (LD) and haplotype analyses to fine map the 15q22 locus. Luciferase reporter assays were applied to evaluate the regulatory effects of the candidate variants. Knockdown by small interfering RNA, microarray analysis, chromatin immunoprecipitation (ChIP) and quantitative real-time polymerase chain reaction assays were performed to reveal the regulatory gene network and identify its binding targets. We report a 25.6-kb haplotype within SMAD3 containing numerous single-nucleotide polymorphisms (SNPs) in high LD. SNPs rs17293632 and rs4562997 were identified as functional variants of SMAD3 by luciferase assays within the LD region. These variants regulate SMAD3 transcription in an allele-specific manner through enhancer elements in introns of SMAD3. Knockdown of SMAD3 in thyroid cancer cell lines revealed its regulatory gene network including two upregulated genes, SPRY4 and SPRY4-IT1. Sequence analysis and ChIP assays validated the actual binding of SMAD3 protein to multiple SMAD binding element sites in the region upstream of SPRY4. Our data provide a functional annotation of the 15q22 thyroid cancer risk locus.

Background: Familial non-medullary thyroid cancer (NMTC) accounts for a relatively small proportion of thyroid cancer cases, but it displays strong genetic predisposition. So far, only a few NMTC susceptible genes and low-penetrance variants contributing to NMTC have been described. This study aimed to identify rare germline variants that may predispose individuals to NMTC by sequencing a cohort of 17 NMTC families. Methods: Whole-genome sequencing and genome-wide linkage analysis were performed in 17 NMTC families. MendelScan and BasePlayer were applied to screen germline variants followed by customized filtering. The remaining candidate variants were subsequently validated by Sanger sequencing. A panel of 277 known cancer predisposition genes was also screened in these families. Results: A total of 41 rare coding candidate variants in 40 genes identified by whole-genome sequencing are reported, including 24 missense, five frameshift, five splice change, and seven nonsense variants. Sanger sequencing confirmed all 41 rare variants and proved their co-segregation with NMTC in the extended pedigrees. In silico functional analysis of the candidate genes using Ingenuity Pathway Analysis showed that cancer was the top category of "Diseases and Disorders." Additionally, a targeted search displayed six variants in known cancer predisposition genes, including one frameshift variant and five missense variants. Conclusions: The data identify rare germline variants that may play important roles in NMTC predisposition. It is proposed that in future research including functional characterization, these variants and genes be considered primary candidates for thyroid cancer predisposition.

Long noncoding RNAs (lncRNAs) regulate pathological processes, yet their potential roles in papillary thyroid carcinoma (PTC) are poorly understood.To profile transcriptionally dysregulated lncRNAs in PTC and identify lncRNAs associated with clinicopathological characteristics.We performed RNA sequencing of 12 paired PTC tumors and matched noncancerous tissues and correlated the expression of lncRNAs with clinical parameters. The 2 most significantly dysregulated lncRNAs were studied in an Ohio PTC cohort (n = 109) and in PTC data (n = 497) from The Cancer Genome Atlas.A combination of laboratory-based studies and computational analysis using clinical data and samples and a publically available database.Correlation between expression values and clinical parameters.We identified 218 lncRNAs showing differential expression in PTC (fold change ≥ 2.0, P.01). Significant correlation was observed between the expression of 2 lncRNAs (XLOC_051122 and XLOC_006074) and 1) lymph node metastasis (N stage) and 2) BRAF(V600E) mutation. Among patients with wild-type BRAF, the expression of these 2 lncRNAs showed significantly higher levels in the patients with lymph node metastasis. In silico analysis of these lncRNAs pinpointed cell movement and cellular growth and proliferation as targeted functions.Comprehensive expression screening identified 2 novel lncRNAs associated with risk factors of adverse prognosis in PTC patients. These lncRNAs may be novel players in PTC carcinogenesis.

Chromosomal aberrations and multiple genome-wide association studies (GWAS) have established a major hematopoietic quantitative trait locus in chromosome 6q23.3. The locus comprises an active enhancer region, in which some of the associated SNPs alter transcription factor binding. We now identify miR-3662 as a new functional driver contributing to the associated phenotypes. The GWAS SNPs are strongly associated with higher miR-3662 expression. Genome editing of rs66650371, a three-base-pair deletion, suggests a functional link between the SNP genotype and the abundance of miR-3662. Increasing miR-3662′s abundance increases colony formation in hematopoietic progenitor cells, particularly the erythroid lineage. In contrast, miR-3662 is not expressed in acute myeloid leukemia cells, and its overexpression has potent antileukemic effects in vitro and in vivo. Mechanistically, miR-3662 directly targets NF-κB–mediated transcription. Thus, miR-3662 is a new player of the hematopoietic 6q23.3 locus. Significance: The characterization of miR-3662 has identified a new actor in the prominent hematopoietic quantitative trait locus in chromosome 6q23.3. The mechanistic insights into miR-3662′s function may reveal novel or only partially known pathways for normal and malignant hematopoietic cell proliferation. Cancer Discov; 6(9); 1036–51. ©2016 AACR. This article is highlighted in the In This Issue feature, p. 932

We previously showed that a long noncoding RNA gene, PTCSC3, located close to the variant rs944289 that predisposes to papillary thyroid carcinoma (PTC) might target the S100A4 gene.The aim was to investigate the impact of PTCSC3 on S100A4 expression and its role in cancer development.S100A4 abundance was analyzed by quantitative PCR (qPCR) in unaffected and tumor tissue from n = 73 PTC patients. The expression of PTCSC3 and S100A4 was studied in BCPAP and TPC-1 cell lines with forced expression of PTCSC3 by qPCR. Expression of S100A4 target genes (VEGF and MMP-9) was studied in the BCPAP cell line with forced expression of PTCSC3 by qPCR, reverse transcriptase PCR, and Western blot. The impact of PTCSC3 on BCPAP motility and invasiveness was analyzed by the Transwell and Matrigel assays, respectively.This was a laboratory-based study using cells from clinical samples and thyroid cancer cell lines.We aimed to find evidence for a link between the expression of PTCSC3 and thyroid carcinogenesis.Expression data from PTC cell lines pinpointed S100A4 as the most significantly downregulated gene in the presence of PTCSC3. S100A4 was upregulated in tumor tissue (P = 9.33 × 10(-7)) while PTCSC3 was strongly downregulated (P = 2.2 × 10(-16)). S100A4 transcription was moderately correlated with PTCSC3 expression in unaffected thyroid tissue (r = 0.429, P = .0001), and strongly in unaffected tissue of patients with the risk allele of rs944289 (r = 0.685, P = 7.88 × 10(-5)). S100A4, VEGF, and MMP-9 were suppressed in the presence of PTCSC3 (P = .0051, P = .0090, and P =.0037, respectively). PTC cells expressing PTCSC3 showed reduction in motility and invasiveness (P = 4.52 × 10(-5) and P = 1.0 × 10(-4), respectively).PTCSC3 downregulates S100A4, leading to a reduction in cell motility and invasiveness. We propose that PTCSC3 impacts PTC predisposition and carcinogenesis through the S100A4 pathway.

The [A] allele of SNP rs965513 in 9q22 has been consistently shown to be highly associated with increased papillary thyroid cancer (PTC) risk with an odds ratio of ∼1.8 as determined by genome-wide association studies, yet the molecular mechanisms remain poorly understood. Previously, we noted that the expression of two genes in the region, forkhead box E1 (FOXE1) and PTC susceptibility candidate 2 (PTCSC2), is regulated by rs965513 in unaffected thyroid tissue, but the underlying mechanisms were not elucidated. Here, we fine-mapped the 9q22 region in PTC and controls and detected an ∼33-kb linkage disequilibrium block (containing the lead SNP rs965513) that significantly associates with PTC risk. Chromatin characteristics and regulatory element signatures in this block disclosed at least three regulatory elements functioning as enhancers. These enhancers harbor at least four SNPs (rs7864322, rs12352658, rs7847449, and rs10759944) that serve as functional variants. The variant genotypes are associated with differential enhancer activities and/or transcription factor binding activities. Using the chromosome conformation capture methodology, long-range looping interactions of these elements with the promoter region shared by FOXE1 and PTCSC2 in a human papillary thyroid carcinoma cell line (KTC-1) and unaffected thyroid tissue were found. Our results suggest that multiple variants coinherited with the lead SNP and located in long-range enhancers are involved in the transcriptional regulation of FOXE1 and PTCSC2 expression. These results explain the mechanism by which the risk allele of rs965513 predisposes to thyroid cancer.

Context: By genome-wide association studies, the risk allele [A] of SNP rs965513 predisposes strongly to papillary thyroid carcinoma (PTC). It is located in a gene-poor region of 9q22, some 60 kb from the FOXE1 gene. The underlying mechanisms remain to be discovered. Objective: Our objective was to identify novel transcripts in the 9q22 locus and correlate gene expression levels with the genotypes of rs965513. Design: We performed 3′ and 5′ rapid amplification of cDNA ends and RT-PCR to detect novel transcripts. One novel transcript was forcibly expressed in a cell line followed by gene expression array analysis. We genotyped rs965513 from PTC patients and measured gene expression levels by real-time RT-PCR in unaffected thyroid tissue and matched tumor. Setting: This was a laboratory-based study using cells from clinical tissue samples and a cancer cell line. Main Outcome Measures: We detected previously uncharacterized transcripts and evaluated the gene expression levels and the correlation with the risk allele of rs965513, age, gender, chronic lymphocyte thyroiditis (CLT), and TSH levels. Results: We found a novel long intergenic noncoding RNA gene and named it papillary thyroid cancer susceptibility candidate 2 (PTCSC2). Transcripts of PTCSC2 are down-regulated in PTC tumors. The risk allele [A] of rs965513 was significantly associated with low expression of unspliced PTCSC2, FOXE1, and TSHR in unaffected thyroid tissue. We also observed a significant association of age and CLT with PTCSC2 unspliced transcript levels. The correlation between the rs965513 genotype and the PTCSC2 unspliced transcript levels remained significant after adjusting for age, gender, and CLT. Forced expression of PTCSC2 in the BCPAP cell line affected the expression of a subset of noncoding and coding transcripts with enrichment of genes functionally involved in cell cycle and cancer. Conclusions: Our data suggest a role for PTCSC2, FOXE1, and TSHR in the predisposition to PTC.

Neuroblastoma rat sarcoma (RAS) viral oncogene homolog (NRAS), a small GTPase, is one of the most thoroughly studied oncogenes that controls cell growth, differentiation, and survival by facilitating signal transduction. Here, we identify four novel naturally occurring NRAS isoforms (isoforms 2-5) in addition to the canonical isoform (isoform 1). Expression analyses performed on a panel of several different human malignancies and matching normal tissue revealed distinct isoform expression patterns. Two of the novel isoforms were found in the nucleus and cytoplasm, whereas the others were exclusively cytoplasmic. The isoforms varied in their binding affinities to known downstream targets and differentially regulated the RAS signaling pathway. Strikingly, forced expression of isoform 5, which encodes only a 20-aa peptide, led to increased cell proliferation and to transformation by activation of known NRAS targets. These discoveries open new avenues in the study of NRAS.

Germline mutations in PMS2 are associated with Lynch syndrome (LS), the most common known cause of hereditary colorectal cancer. Mutation detection in PMS2 has been difficult due to the presence of several pseudogenes, but a custom-designed long-range PCR strategy now allows adequate mutation detection. Many mutations are unique. However, some mutations are observed repeatedly across individuals not known to be related due to the mutation being either recurrent, arising multiple times de novo at hot spots for mutations, or of founder origin, having occurred once in an ancestor. Previously, we observed 36 distinct mutations in a sample of 61 independently ascertained Caucasian probands of mixed European background with PMS2 mutations. Eleven of these mutations were detected in more than one individual not known to be related and of these, six were detected more than twice. These six mutations accounted for 31 (51%) ostensibly unrelated probands. Here, we performed genotyping and haplotype analysis in four mutations observed in multiple probands and found two (c.137G>T and exon 10 deletion) to be founder mutations and one (c.903G>T) a probable founder. One (c.1A>G) could not be evaluated for founder mutation status. We discuss possible explanations for the frequent occurrence of founder mutations in PMS2.

The great majority of thyroid cancers are of the non-medullary type. Here we report findings from a genome-wide association study of non-medullary thyroid cancer, including in total 3,001 patients and 287,550 controls from five study groups of European descent. Our results yield five novel loci (all with Pcombined, We thank the patients whose contribution made this work possible. The study was supported in part by grants from the National Institutes of Health: CA16058 (to the Ohio State University Comprehensive Cancer Center) and CA168505 (to M.D.R. and A.d.l.C.). This work was supported in part by The University of Texas MD Anderson Cancer Center start-up funds (to E.M.S.); an American Thyroid Association research grant (to E.M.S.), National Institutes of Health (NIH) grant U01 DE019765-01 (to E.M.S.); and Cancer Center Support Grant CA016672 to the University of Texas MD Anderson Cancer Center. This study makes use of data generated by the Genome of the Netherlands Project (GoNL). A full list of the investigators is available from http://www.nlgenome.nl. Funding for the GoNL project was provided by the Netherlands Organization for Scientific Research under award number 184021007, dated 9 July 2009 and made available as a Rainbow Project of the Biobanking and Biomolecular Research Infrastructure Netherlands (BBMRI-NL). The sequencing was carried out in collaboration with the Beijing Institute for Genomics (BGI). This study also makes use of data generated by the GTEx Consortium and we are grateful for early stage sharing of RNA-seq and genotypic data. The following acknowledgement refers to the GTEx data: The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health. Additional funds were provided by the NCI, NHGRI, NHLBI, NIDA, NIMH and NINDS. Donors were enrolled at Biospecimen Source Sites funded by NCI\SAIC-Frederick, Inc. (SAIC-F) subcontracts to the National Disease Research Interchange (10XS170), Roswell Park Cancer Institute (10XS171) and Science Care, Inc. (X10S172). The Laboratory, Data Analysis and Coordinating Center (LDACC) was funded through a contract (HHSN268201000029C) to The Broad Institute, Inc. Biorepository operations were funded through an SAIC-F subcontract to Van Andel Institute (10ST1035). Additional data repository and project management were provided by SAIC-F (HHSN261200800001E). The Brain Bank was supported by supplements to University of Miami grants DA006227 and DA033684 and to contract N01MH000028. Statistical Methods development grants were made to the University of Geneva (MH090941 and MH101814), the University of Chicago (MH090951, MH090937, MH101820 and MH101825), the University of North Carolina—Chapel Hill (MH090936 and MH101819), Harvard University (MH090948), Stanford University (MH101782), Washington University St Louis (MH101810) and the University of Pennsylvania (MH101822).

High BAALC expression levels are associated with poor outcome in cytogenetically normal acute myeloid leukemia (CN-AML) patients. Recently, miR-3151 was discovered in intron 1 of BAALC. To evaluate the prognostic significance of miR-3151 expression levels and to gain insight into the biologic and prognostic interplay between miR-3151 and its host, miR-3151 and BAALC expression were measured in pretreatment blood of 179 CN-AML patients. Gene-expression profiling and miRNA-expression profiling were performed using microarrays. High miR-3151 expression was associated with shorter disease-free and overall survival, whereas high BAALC expression predicted failure of complete remission and shorter overall survival. Patients exhibiting high expression of both miR-3151 and BAALC had worse outcome than patients expressing low levels of either gene or both genes. In gene-expression profiling, high miR-3151 expressers showed down-regulation of genes involved in transcriptional regulation, posttranslational modification, and cancer pathways. Two genes, FBXL20 and USP40, were validated as direct miR-3151 targets. The results of the present study show that high expression of miR-3151 is an independent prognosticator for poor outcome in CN-AML and affects different outcome end points than its host gene, BAALC. The combination of both markers identified a patient subset with the poorest outcome. This interplay between an intronic miR and its host may have important biologic implications.

Although the physiological basis of canonical or classical IκB kinase β (IKKβ)–nuclear factor κB (NF-κB) signaling pathway is well established, how alternative NF-κB signaling functions beyond its role in lymphoid development remains unclear. In particular, alternative NF-κB signaling has been linked with cellular metabolism, but this relationship is poorly understood. In this study, we show that mice deleted for the alternative NF-κB components IKKα or RelB have reduced mitochondrial content and function. Conversely, expressing alternative, but not classical, NF-κB pathway components in skeletal muscle stimulates mitochondrial biogenesis and specifies slow twitch fibers, suggesting that oxidative metabolism in muscle is selectively controlled by the alternative pathway. The alternative NF-κB pathway mediates this specificity by direct transcriptional activation of the mitochondrial regulator PPAR-γ coactivator 1β (PGC-1β) but not PGC-1α. Regulation of PGC-1β by IKKα/RelB also is mammalian target of rapamycin (mTOR) dependent, highlighting a cross talk between mTOR and NF-κB in muscle metabolism. Together, these data provide insight on PGC-1β regulation during skeletal myogenesis and reveal a unique function of alternative NF-κB signaling in promoting an oxidative metabolic phenotype.

Context: The family risk ratio for papillary thyroid carcinoma (PTC) is among the highest of all cancers. Collectively, familial cases (fPTC) and sporadic cases (sPTC) are not known to show molecular differences. However, one study reported that telomeres were markedly shorter and the telomerase reverse transcriptase (TERT) gene was amplified and up-regulated in germline DNA from patients with fPTC compared with sPTC. Objective: The aim of this study was to evaluate telomere length and TERT gene amplification and expression in blood samples of fPTC and sPTC patients in a genetically distinct population from the previous study. Design: In 42 fPTC and 65 sPTC patients, quantitative real-time PCR was employed to measure the relative telomere length (RTL) and TERT gene copy number and RNA level. To validate the results using alternative methods, we further studied a subset of the original cohort consisting of randomly chosen fPTC (n = 10) and sPTC (n = 14) patients and controls (n = 21) by assessing both telomere length by flow fluorescent in situ hybridization and TERT gene expression by quantitative real-time PCR. Results: RTL and TERT gene copy number did not differ between fPTC and sPTC (P = 0.957 and P = 0.998, respectively). The mean RTL and TERT gene expression were not significantly different among the groups of the validation series (P = 0.169 and P = 0.718, respectively). Conclusion: Our data show no difference between familial and sporadic PTC with respect to telomere length, TERT copy number, or expression in our cohort. Further investigations in additional cohorts of patients are desirable.

Mutations in the mismatch repair genes cause Lynch syndrome (LS), conferring high risk of colorectal, endometrial and some other cancers. After the same splice site mutation in the MLH1 gene (c.589-2A>G) had been observed in four ostensibly unrelated American families with typical LS cancers, its occurrence in comprehensive series of LS cases (Mayo Clinic, Germany and Italy) was determined. It occurred in 10 out of 995 LS mutation carriers (1.0%) diagnosed in the Mayo Clinic diagnostic laboratory. It did not occur among 1,803 cases tested for MLH1 mutations by the German HNPCC consortium, while it occurred in three probands and an additional five family members diagnosed in Italy. In the U.S., the splice site mutation occurs on a large (~4.8 Mb) shared haplotype that also harbors the variant c.2146G>A, which predicts a missense change in codon 716 referred to here as V716M. In Italy, it occurs on a different, shorter shared haplotype (~2.2 Mb) that does not carry V716M. The V716M variant was found to be present by itself in the U.S., German and Italian populations with individuals sharing a common haplotype of 280 kb, allowing us to calculate that the variant arose around 5,600 years ago (225 generations; 95% confidence interval 183–272). The splice site mutation in America arose or was introduced some 450 years ago (18 generations; 95% confidence interval 14–23); it accounts for 1.0% all LS in the Unites States and can be readily screened for.

Trimethylation of histone 3 lysine 27 (H3K27me3) is a critical epigenetic mark for the maintenance of gene silencing. Additional accumulation of DNA methylation in target loci is thought to cooperatively support this epigenetic silencing during tumorigenesis. However, molecular mechanisms underlying the complex interplay between the two marks remain to be explored. Here we show that activation of PI3K/AKT signaling can be a trigger of this epigenetic processing at many downstream target genes. We also find that DNA methylation can be acquired at the same loci in cancer cells, thereby reinforcing permanent repression in those losing the H3K27me3 mark. Because of a link between PI3K/AKT signaling and epigenetic alterations, we conducted epigenetic therapies in conjunction with the signaling-targeted treatment. These combined treatments synergistically relieve gene silencing and suppress cancer cell growth in vitro and in xenografts. The new finding has important implications for improving targeted cancer therapies in the future. Cancer Res; 71(5); 1752–62. ©2011 AACR.

Substantial evidence indicates that exposure to bisphenol A (BPA) during early development may increase breast cancer risk later in life. The changes may persist into puberty and adulthood, suggesting an epigenetic process being imposed in differentiated breast epithelial cells. The molecular mechanisms by which early memory of BPA exposure is imprinted in breast progenitor cells and then passed onto their epithelial progeny are not well understood. The aim of this study was to examine epigenetic changes in breast epithelial cells treated with low-dose BPA. We also investigated the effect of BPA on the ERα signaling pathway and global gene expression profiles. Compared to control cells, nuclear internalization of ERα was observed in epithelial cells preexposed to BPA. We identified 170 genes with similar expression changes in response to BPA. Functional analysis confirms that gene suppression was mediated in part through an ERα-dependent pathway. As a result of exposure to BPA or other estrogen-like chemicals, the expression of lysosomal-associated membrane protein 3 (LAMP3) became epigenetically silenced in breast epithelial cells. Furthermore, increased DNA methylation in the LAMP3 CpG island was this repressive mark preferentially occurred in ERα-positive breast tumors. These results suggest that the in vitro system developed in our laboratory is a valuable tool for exposure studies of BPA and other xenoestrogens in human cells. Individual and geographical differences may contribute to altered patterns of gene expression and DNA methylation in susceptible loci. Combination of our exposure model with epigenetic analysis and other biochemical assays can give insight into the heritable effect of low-dose BPA in human cells.

The genetic component of colorectal cancer (CRC) predisposition has been only partially explained. We recently suggested that a subtle decrease in the expression of one allele of the TGFBR1 gene was a heritable quantitative trait predisposing to CRC. Here, we refined the measurements of allele-specific expression (ASE) of TGFBR1 in a population-based series of CRC patients and controls. Five single-nucleotide polymorphisms (SNPs) in the 3′-untranslated region of the gene were genotyped and used for ASE determination by pyrosequencing. After eliminating non-informative samples and samples with RNA of insufficient quality 109 cases and 125 controls were studied. Allelic ratios ranged between 0.74 and 1.69 without evidence of bimodality or cutoff points for ‘ASE’ versus ‘non-ASE’. Treating ASE as a continuous variable, cases had non-significantly different values than controls (P = 0.081 when comparing means by permutation test). However, cases had significantly higher ASE values when comparing medians by permutation test (P = 0.0027) and when using Wilcoxon test (P = 0.0094). We conclude that with the present-day technology, ASE differences between individuals and between cases and controls are too subtle to be used to assess CRC risk. More advanced technology is expected to resolve this issue as well as the low informativity caused by the limited heterozygosity of transcribed SNPs.

Nuclear factor (NF)-kappaB is a positive regulator of tumour development and progression, but how it functions in normal cells leading to oncogenesis is not clear. As cellular senescence has proven to be an intrinsic tumour suppressor mechanism that cells must overcome to establish deregulated growth, we used primary fibroblasts to follow NF-kappaB function in cells transitioning from senescence to subsequent immortalization. Our findings show that RelA/p65(-/-) murine fibroblasts immortalize at considerably faster rates than RelA/p65(+/+) cells. The ability of RelA/p65(-/-) fibroblasts to escape senescence earlier is due to their genomic instability, characterized by high frequencies of DNA mutations, gene deletions and gross chromosomal translocations. This increase in genomic instability is closely related to a compromised DNA repair that occurs in both murine RelA/p65(-/-) fibroblasts and tissues. Significantly, these results can also be duplicated in human fibroblasts lacking NF-kappaB. Altogether, our findings present a fresh perspective on the role of NF-kappaB as a tumour suppressor, which acts in pre-neoplastic cells to maintain cellular senescence by promoting DNA repair and genomic stability.

Early exposure to xenoestrogens may predispose to breast cancer risk later in adult life. It is likely that long-lived, self-regenerating epithelial progenitor cells are more susceptible to these exposure injuries over time and transmit the injured memory through epigenetic mechanisms to their differentiated progeny. Here, we used progenitor-containing mammospheres as an in vitro exposure model to study this epigenetic effect. Expression profiling identified that, relative to control cells, 9.1% of microRNAs (82 of 898 loci) were altered in epithelial progeny derived from mammospheres exposed to a synthetic estrogen, diethylstilbestrol. Repressive chromatin marks, trimethyl Lys27 of histone H3 (H3K27me3) and dimethyl Lys9 of histone H3 (H3K9me2), were found at a down-regulated locus, miR-9-3, in epithelial cells preexposed to diethylstilbestrol. This was accompanied by recruitment of DNA methyltransferase 1 that caused an aberrant increase in DNA methylation of its promoter CpG island in mammosphere-derived epithelial cells on diethylstilbestrol preexposure. Functional analyses suggest that miR-9-3 plays a role in the p53-related apoptotic pathway. Epigenetic silencing of this gene, therefore, reduces this cellular function and promotes the proliferation of breast cancer cells. Promoter hypermethylation of this microRNA may be a hallmark for early breast cancer development, and restoration of its expression by epigenetic and microRNA-based therapies is another viable option for future treatment of this disease. [Cancer Res 2009;69(14):5936–45]

Papillary thyroid carcinoma (PTC) displays higher heritability than most other cancers. To search for genes predisposing to PTC, we performed a genome-wide linkage analysis in a large family with PTC and melanoma. Among several peaks the highest was at 8q24, with a maximum nonparametric linkage (NPL) score of 7.03. Linkage analysis was then broadened to comprise 25 additional PTC families that produced a maximum NPL score of 3.2, P = 0.007 at the 8q24 locus. Fine mapping with microsatellite markers was compatible with linkage to the 8q24 locus in 10 of the 26 families. In the large family, a ∼320 Kb haplotype was shared by individuals with PTC, melanoma, or benign thyroid disease, but not by unaffected individuals. A 12 Kb haplotype of 8 SNP markers within the larger haplotype was shared by 9 of the 10 families in which the 8q24 locus was compatible with linkage. The shared haplotype is located within 2 known overlapping protein-coding genes, thyroglobulin (TG) and Src-like adaptor (SLA). Resequencing of the coding and control regions of TG and SLA did not disclose putative mutations in PTC patients. Embedded in the TG-SLA region are three likely noncoding RNA genes, one of which (AK023948) harbors the 8-SNP haplotype. Resequencing of AK023948 and one of the other RNA genes did not reveal candidate mutations. Gene expression analysis indicated that AK023948 is significantly down-regulated in most PTC tumors. The putative noncoding RNA gene AK023948 is a candidate suseptibility gene for PTC. [Cancer Res 2009;69(2):625–31]

Posttranslationally modified histones and DNA hypermethylation frequently interplay to deregulate gene expression in cancer. We report that acute myeloid leukemia (AML) with an aberrant histone methyltransferase, the mixed lineage leukemia partial tandem duplication (MLL-PTD), exhibits increased global DNA methylation versus AML with MLL-wildtype (MLL-WT; P = .02). Among the differentially methylated genes, the SLC5A8 tumor suppressor gene (TSG) was more frequently hypermethylated (P = .003). In MLL-PTD+ cell lines having SLC5A8 promoter hypermethylation, incubation with decitabine activated SLC5A8 expression. Ectopic SLC5A8 expression enhanced histones H3 and H4 acetylation in response to the histone deacetylase inhibitor, valproate, consistent with the encoded protein—SMCT1—short-chain fatty acid transport function. In addition, enhanced cell death was observed in SMCT1-expressing MLL-PTD+ AML cells treated with valproate. Within the majority of MLL-PTD AML is a mechanism in which DNA hypermethylation silences a TSG that, together with MLL-PTD, can contribute further to aberrant chromatin remodeling and altered gene expression.

Large germline deletions within the mismatch repair gene MSH2 account for a significant proportion (up to 20%) of all deleterious mutations of this gene which are associated with Lynch syndrome. An exons 1 to 6 deletion of MSH2, originally reported in nine families, has been associated with a founding event within the United States, which genealogic studies had previously dated to 1727, and the number of present day carriers was estimated to be 18,981. Here, we report the development of a robust multiplex PCR which has assisted in the detection of 32 new families who carry the MSH2 American Founder Mutation (AFM). By offering testing to family members, 126 carriers of the AFM have been identified. Extensive genealogic studies have connected 27 of the 41 AFM families into seven extended pedigrees. These extended families have been traced back to around the 18th century without any evidence of further convergence between them. Characterization of the genomic sequence flanking the deletion and the identification of a common disease haplotype of between 0.6 and 2.3 Mb in all probands provides evidence for a common ancestor between these extended families. The DMLE+2.2 software predicts an age of ∼500 years (95% confidence interval, 425–625) for this mutation. Taken together, these data are suggestive of an earlier founding event than was first thought, which likely occurred in a European or a Native American population. The consequences of this finding would be that the AFM is significantly more frequent in the United States than was previously predicted. [Cancer Res 2008;68(7):2145–53]

Hypermethylation of CpG islands is a common epigenetic alteration associated with cancer. Global patterns of hypermethylation are tumor-type specific and nonrandom. The biological significance and the underlying mechanisms of tumor-specific aberrant promoter methylation remain unclear, but some evidence suggests that this specificity involves differential sequence susceptibilities, the targeting of DNA methylation activity to specific promoter sequences, or the selection of rare DNA methylation events during disease progression. Using restriction landmark genomic scanning on samples derived from tissue culture and in vivo models of T cell lymphomas, we found that MYC overexpression gave rise to a specific signature of CpG island hypermethylation. This signature reflected gene transcription profiles and was detected only in advanced stages of disease. The further inactivation of the Pten, p53, and E2f2 tumor suppressors in MYC-induced lymphomas resulted in distinct and diagnostic CpG island methylation signatures. Our data suggest that tumor-specific DNA methylation in lymphomas arises as a result of the selection of rare DNA methylation events during the course of tumor development. This selection appears to be driven by the genetic configuration of tumor cells, providing experimental evidence for a causal role of DNA hypermethylation in tumor progression and an explanation for the tremendous epigenetic heterogeneity observed in the evolution of human cancers. The ability to predict genome-wide epigenetic silencing based on relatively few genetic alterations will allow for a more complete classification of tumors and understanding of tumor cell biology., Author Summary Genetic and epigenetic alterations of the genome are common features of cancers. The relationship between these two types of alterations, however, remains unclear. One type of epigenetic modification—DNA methylation in promoter sequences of genes—is of particular interest, since tumor cells have different patterns of promoter methylation than normal cells. Previous studies on human tumor samples have suggested a link between genetic alterations and the induction of aberrant DNA methylation; however, this link has been difficult to rigorously assess because of the incredible genetic heterogeneity found in human cancer. In this study, a mouse model of T cell lymphoma was used to explore the relationship between genetic and epigenetic modifications experienced by tumor cells. By introducing defined genetic changes into preneoplastic T cells of mice, such as the overexpression of the MYC oncogene and the ablation of tumor suppressor genes, we could carefully evaluate how these genetic changes impacted promoter methylation profiles during development of lymphomas in vivo. We found that the introduction of different genetic insults resulted in unique and diagnostic profiles of promoter methylation. Understanding how these methylation signatures contribute to tumor progression could eventually have diagnostic, prognostic, and therapeutic value for human cancers.

Papillary thyroid carcinomas (PTCs) that invade into local structures are associated with a poor prognosis, but the mechanisms for PTC invasion are incompletely defined, limiting the development of new therapies. To characterize biological processes involved in PTC invasion, we analyzed the gene expression profiles of microscopically dissected intratumoral samples from central and invasive regions of seven widely invasive PTCs and normal thyroid tissue by oligonucleotide microarray and performed confirmatory expression and functional studies. In comparison with the central regions of primary PTCs, the invasive fronts overexpressed TGF β, NFκB and integrin pathway members, and regulators of small G proteins and CDC42. Moreover, reduced levels of mRNAs encoding proteins involved in cell–cell adhesion and communication were identified, consistent with epithelial-to-mesenchymal transition (EMT). To confirm that aggressive PTCs were characterized by EMT, 34 additional PTCs were examined for expression of vimentin, a hallmark of EMT. Overexpression of vimentin was associated with PTC invasion and nodal metastasis. Functional, in vitro studies demonstrated that vimentin was required both for the development and maintenance of a mesenchymal morphology and invasiveness in thyroid cancer cells. We conclude that EMT is common in PTC invasion and that vimentin regulates thyroid cancer EMT in vitro .

The B-Raf proto-oncogene serine/threonine kinase (BRAF) gene is the most frequently mutated gene in malignant melanoma (MM) and papillary thyroid cancer (PTC) and is causally involved in malignant cell transformation. Mutated BRAF is associated with an aggressive disease phenotype, thus making it a top candidate for targeted treatment strategies in MM and PTC. We show that BRAF mutations in both MM and PTC drive increased expression of oncomiR-3151, which is coactivated by the SP1/NF-κB complex. Knockdown of microRNA-3151 (miR-3151) with short hairpin RNAs reduces cell proliferation and increases apoptosis of MM and PTC cells. Using a targeted RNA sequencing approach, we mechanistically determined that miR-3151 directly targets TP53 and other members of the TP53 pathway. Reducing miR-3151's abundance increases TP53's mRNA and protein expression and favors its nuclear localization. Consequently, knockdown of miR-3151 also leads to caspase-3-dependent apoptosis. Simultaneous inhibition of aberrantly activated BRAF and knockdown of miR-3151 potentiates the effects of sole BRAF inhibition with the BRAF inhibitor vemurafenib and may provide a novel targeted therapeutic approach in BRAF-mutated MM and PTC patients. In conclusion, we identify miR-3151 as a previously unidentified player in MM and PTC pathogenesis, which is driven by BRAF-dependent and BRAF-independent mechanisms. Characterization of TP53 as a downstream effector of miR-3151 provides evidence for a causal link between BRAF mutations and TP53 inactivation.

The main nonmedullary form of thyroid cancer is papillary thyroid carcinoma (PTC) that accounts for 80–90% of all thyroid malignancies. Only 3–10% of PTC patients have a positive family history of PTC yet the familiality is one of the highest of all cancers as measured by case control studies. A handful of genes have been implicated accounting for a small fraction of this genetic predisposition. It was therefore of considerable interest that a mutation in the HABP2 gene was recently implicated in familial PTC. The present work was undertaken to examine the extent of HABP2 variant involvement in PTC. The HABP2 G534E variant (rs7080536) was genotyped in blood DNA from 179 PTC families (one affected individual per family), 1160 sporadic PTC cases and 1395 controls. RNA expression of HABP2 was tested by qPCR in RNA extracted from tumor and normal thyroid tissue from individuals that are homozygous wild-type or heterozygous for the variant. The variant was found to be present in 6.1% familial cases, 8.0% sporadic cases (2 individuals were homozygous for the variant) and 8.7% controls. The variant did not segregate with PTC in one large and 6 smaller families in which it occurred. In keeping with data from the literature and databases the expression of HABP2 was highest in the liver, much lower in 3 other tested tissues (breast, kidney, brain) but not found in thyroid. Given these results showing lack of any involvement we suggest that the putative role of variant HABP2 in PTC should be carefully scrutinized.

Background Traditionally, the CD56dimCD16+ subset of Natural Killer (NK) cells has been thought to mediate cellular cytotoxicity with modest cytokine secretion capacity. However, studies have suggested that this subset may exert a more diverse array of immunological functions. There exists a lack of well-developed functional models to describe the behavior of activated NK cells, and the interactions between signaling pathways that facilitate effector functions are not well understood. In the present study, a combination of genome-wide microarray analyses and systems-level bioinformatics approaches were utilized to elucidate the transcriptional landscape of NK cells activated via interactions with antibody-coated targets in the presence of interleukin-12 (IL-12). Methods We conducted differential gene expression analysis of CD56dimCD16+ NK cells following FcR stimulation in the presence or absence of IL-12. Next, we functionally characterized gene sets according to patterns of gene expression and validated representative genes using RT-PCR. IPA was utilized for biological pathway analysis, and an enriched network of interacting genes was generated using GeneMANIA. Furthermore, PAJEK and the HITS algorithm were employed to identify important genes in the network according to betweeness centrality, hub, and authority node metrics. Results Analyses revealed that CD56dimCD16+ NK cells co-stimulated via the Fc receptor (FcR) and IL-12R led to the expression of a unique set of genes, including genes encoding cytotoxicity receptors, apoptotic proteins, intracellular signaling molecules, and cytokines that may mediate enhanced cytotoxicity and interactions with other immune cells within inflammatory tissues. Network analyses identified a novel set of connected key players, BATF, IRF4, TBX21, and IFNG, within an integrated network composed of differentially expressed genes in NK cells stimulated by various conditions (immobilized IgG, IL-12, or the combination of IgG and IL-12). Conclusions These results are the first to address the global mechanisms by which NK cells mediate their biological functions when encountering antibody-coated targets within inflammatory sites. Moreover, this study has identified a set of high-priority targets for subsequent investigation into strategies to combat cancer by enhancing the anti-tumor activity of CD56dimCD16+ NK cells. Electronic supplementary material The online version of this article (doi:10.1186/s12920-015-0142-9) contains supplementary material, which is available to authorized users.

Papillary thyroid carcinoma (PTC) displays strong but so far largely uncharacterized heritability. Here we studied genetic predisposition in a family with six affected individuals. We genotyped all available family members and conducted whole exome sequencing of blood DNA from two affected individuals. Haplotype analysis and other genetic criteria narrowed our list of candidates to a germline variant in the serine/arginine repetitive matrix 2 gene (SRRM2). This heterozygous variant, c.1037C > T (Ser346Phe or S346F; rs149019598) cosegregated with PTC in the family. It was not found in 138 other PTC families. It was found in 7/1,170 sporadic PTC cases and in 0/1,404 controls (p = 0.004). The encoded protein SRRM2 (also called SRm300) is part of the RNA splicing machinery. To evaluate the possibility that the S346F missense mutation affects alternative splicing, we compared RNA-Seq data in leukocytes from three mutation carriers and three controls. Significant differences in alternative splicing were identified for 1,642 exons, of which a subset of 7 exons was verified experimentally. The results confirmed a higher ratio of inclusion of exons in mutation carriers. These data suggest that the S346F mutation in SRRM2 predisposes to PTC by affecting alternative splicing of unidentified downstream target genes.

In breast cancer and normal estrogen target tissues, estrogen receptor-alpha (ERalpha) signaling results in the establishment of spatiotemporal patterns of gene expression. Whereas primary target gene regulation by ERalpha involves recruitment of coregulatory proteins, coactivators, or corepressors, activation of these downstream promoters by receptor signaling may also involve partnership of ERalpha with other transcription factors. By using an integrated, genome-wide approach that involves ChIP-chip and computational modeling, we uncovered 13 ERalpha-responsive promoters containing both ERalpha and c-MYC binding elements located within close proximity (13-214 bp) to each other. Estrogen stimulation enhanced the c-MYC-ERalpha interaction and facilitated the association of ERalpha, c-MYC, and the coactivator TRRAP with these estrogen-responsive promoters, resulting in chromatin remodeling and increased transcription. These results suggest that ERalpha and c-MYC physically interact to stabilize the ERalpha-coactivator complex, thereby permitting other signal transduction pathways to fine-tune estrogen-mediated signaling networks.

Purpose: Repetitive ribosomal DNA (rDNA) genes are GC-rich clusters in the human genome. The aim of the study was to determine the methylation status of two rDNA subunits, the 18S and 28S genes, in ovarian tumors and to correlate methylation levels with clinicopathologic features in a cohort of ovarian cancer patients. Experimental Design: 18S and 28S rDNA methylation was examined by quantitative methylation-specific PCR in 74 late-stage ovarian cancers, 9 histologically uninvolved, and 11 normal ovarian surface epithelial samples. In addition, methylation and gene expression levels of 18S and 28S rDNAs in two ovarian cancer cell lines were examined by reverse transcription-PCR before and after treatment with the demethylating drug 5′-aza-2′-deoxycytidine. Results: The methylation level (amount of methylated rDNA/β-actin) of 18S and 28S rDNAs was significantly higher (P < 0.05) in tumors than in normal ovarian surface epithelial samples. Methylation of 18S and 28S rDNA was highly correlated (R2 = 0.842). Multivariate analysis by Cox regression found that rDNA hypermethylation [hazard ratio (HR), 0.25; P < 0.01], but not age (HR, 1.29; P = 0.291) and stage (HR, 1.09; P = 0.709), was independently associated with longer progression-free survival. In ovarian cancer cell lines, methylation levels of rDNA correlated with gene down-regulation and 5′-aza-2′-deoxycytidine treatment resulted in a moderate increase in 18S and 28S rDNA gene expressions. Conclusion: This is the first report of rDNA hypermethylation in ovarian tumors. Furthermore, rDNA methylation levels were higher in patients with long progression-free survival versus patients with short survival. Thus, rDNA methylation as a prognostic marker in ovarian cancer warrants further investigation.

Purpose Chronic lymphocytic leukemia (CLL) is a clinically heterogeneous disease for which natural history can be predicted based on the presence or absence of immunoglobulin (Ig) variable heavy chain (VH) gene mutations. Herein we report selective epigenetic silencing of the transcription factor TWIST2 (DERMO1) in Ig VH mutated CLL and describe a semiquantitative assay to study promoter methylation of this gene in primary tumor cells. Materials and Methods TWIST2 promoter methylation was identified by restriction landmark genome scanning. Southern blot (SB), bisulfite sequencing, and combined bisulfite restriction analysis (COBRA), and quantitative SB-COBRA was performed to study methylation of the TWIST2 promoter. Reverse transcription polymerase chain reaction assays were used to study TWIST2 expression in CLL cells. Results Following identification and confirmation of TWIST2 methylation in CLL patients, we demonstrated that expression of this transcription factor is related to the degree of promoter methylation. Expression of TWIST2 in a CLL cell line in which the promoter is methylated was increased following decitabine treatment. We next studied 53 patients by COBRA and demonstrated that 72% of patient samples with mutated Ig VH show TWIST2 methylation, while only 16% of patient samples with unmutated Ig VH were methylated (P < .001). In a subset of patients, methylation of TWIST2 correlated with mRNA expression. Conclusion TWIST2 is differentially methylated in CLL cells relative to Ig VH mutational status and can be quantitatively monitored by SB-COBRA. Based on the known role of TWIST2 in silencing p53 function in other malignancies, future studies should focus on the role of TWIST2 in CLL and related lymphoproliferative diseases.

The estrogen receptor alpha (ERalpha) regulates gene expression by either direct binding to estrogen response elements or indirect tethering to other transcription factors on promoter targets. To identify these promoter sequences, we conducted a genome-wide screening with a novel microarray technique called ChIP-on-chip. A set of 70 candidate ERalpha loci were identified and the corresponding promoter sequences were analyzed by statistical pattern recognition and comparative genomics approaches. We found mouse counterparts for 63 of these loci and classified 42 (67%) as direct ERalpha targets using classification and regression tree (CART) statistical model, which involves position weight matrix and human-mouse sequence similarity scores as model parameters. The remaining genes were considered to be indirect targets. To validate this computational prediction, we conducted an additional ChIP-on-chip assay that identified acetylated chromatin components in active ERalpha promoters. Of the 27 loci upregulated in an ERalpha-positive breast cancer cell line, 20 having mouse counterparts were correctly predicted by CART. This integrated approach, therefore, sets a paradigm in which the iterative process of model refinement and experimental verification will continue until an accurate prediction of promoter target sequences is derived.

Purpose We have previously conducted independent microarray expression analyses of the two most common types of nonmedullary thyroid carcinoma, namely papillary thyroid carcinoma (PTC) and follicular thyroid carcinoma (FTC). In this study, we sought to combine our data sets to shed light on the similarities and differences between these tumor types. Materials and Methods Microarray data from six PTCs, nine FTCs, and 13 normal thyroid samples were normalized to remove interlaboratory variability and then analyzed by unsupervised clustering, t test, and by comparison of absolute and change calls. Expression changes in four genes not previously implicated in thyroid carcinogenesis were verified by reverse transcriptase polymerase chain reaction on these same samples, together with eight additional FTC tumors. Results PTCs showed two distinct groups of genes that were either over- or underexpressed compared with normal thyroid, whereas the predominant changes in FTCs were of decreased expression. Five genes could collectively distinguish the two tumor types. PTCs showed overexpression of CITED1, claudin-10 (CLDN10), and insulin-like growth factor binding protein 6 (IGFBP6) but showed no change in expression of caveolin-1 (CAV1) or -2 (CAV2); conversely, FTCs did not express CLDN10 and had decreased expression of IGFBP6 and/or CAV1 and CAV2. Conclusion PTC and FTC show distinctive microarray expression profiles, suggesting that either they have different molecular origins or they diverge distinctly from a common origin. Furthermore, if verified in a larger series of tumors, these genes could, in combination with known tumor-specific chromosome translocations, form the basis of a valuable diagnostic tool.

CpG island methylation is an epigenetic alteration that contributes to tumorigenesis by transcriptional inactivation of genes. Little is known about the overall levels of CpG island methylation in chronic lymphocytic leukemia (CLL). To provide a baseline estimate of global aberrant methylation and identify target sequences for additional investigation, we performed Restriction Landmark Genomic Scanning on 10 CLL samples. Two methylation-sensitive landmark enzymes were used (NotI and AscI), allowing assessment of over 3000 CpG islands in each sample. Tumor-derived Restriction Landmark Genomic Scanning profiles were compared with profiles from CD19-selected B cells from normal volunteers and matched normal neutrophils from 4 CLL patients. We found 2.5–8.1% (mean 4.8%) of the CpG islands in CLL samples were aberrantly methylated compared with controls, and the methylation events had a nonrandom distribution (P < 0.0001). Furthermore, we identified 193 aberrantly methylated sequences, of which 93% have CpG island characteristics and 90% have homology to genes or expressed sequences. One such gene, the G protein-coupled metabotropic glutamate receptor 7 (GRM7), possibly inhibits cyclic AMP signaling in the induction of apoptosis. Bisulfite sequencing of GRM7 confirmed extensive CpG island methylation, and treatment with 5-aza-2′-deoxycytidine (decitabine) resulted in up-regulated expression of several genes in vitro with concurrent cellular depletion of DNMT1 protein. Our dual-enzyme global methylation study shows that CLL is characterized by widespread nonrandom CpG island methylation similar to other tumors and provides a panel of novel methylation targets that can be used in larger studies designed to assess impact on disease progression and survival.

Molecular mechanisms of leukemogenesis have been successfully unraveled by studying genes involved in simple rearrangements including balanced translocations and inversions. In contrast, little is known about genes altered in complex karyotypic abnormalities. We studied acute myeloid leukemia (AML) patients with complex karyotypes and abnormal chromosome 21. High-resolution bacterial artificial chromosome (BAC) array-based comparative genomic hybridization disclosed amplification predominantly in the 25- to 30-megabase (MB) region that harbors the APP gene (26.3 MB) and at position 38.7-39.1 MB that harbors the transcription factors ERG and ETS2 . Using oligonucleotide arrays, APP was by far the most overexpressed gene (mean fold change 19.74, P = 0.0003) compared to a control group of AML with normal cytogenetics; ERG and ETS2 also ranked among the most highly expressed chromosome 21 genes. Overexpression of APP and ETS2 correlated with genomic amplification, but high APP expression occurred even in a subset of AML patients with normal cytogenetics (10 of 64, 16%). APP encodes a glycoprotein of unknown function previously implicated in Alzheimer's disease, but not in AML. We hypothesize that APP and the transcription factors ERG and ETS2 are altered by yet unknown molecular mechanisms involved in leukemogenesis. Our results highlight the value of molecularly dissecting leukemic cells with complex karyotypes.

Objective The gene BAALC ( Brain And Acute Leukemia, Cytoplasmic ), a novel molecular marker involved in leukemia, is highly expressed in a subset of patients with acute leukemia and predictive of clinical outcome in patients with acute myeloid leukemia and normal karyotype. The role of BAALC in hematopoiesis and leukemogenesis is unknown. Material and methods We used real-time RT-PCR to show that BAALC is strongly expressed in CD34 + cells from the bone marrow and blood and only weakly expressed in total normal bone marrow and blood cells. Results Expression analyses of FACSorted cells revealed high BAALC transcript levels in CD34 + bone marrow cells including CD34 + /CD38 − , CD34 + /CD33 + , as well as CD34 + /CD19 + /CD10 + , CD34 + /CD7 + , and CD34 + /CD71 + /CD45 − cell fractions. Expression was significantly lower in all CD34 − fractions. In vitro differentiation of CD34 + bone marrow cells showed downregulation of BAALC and CD34 transcripts as early as day 4 in suspension cultures supplemented with lineage-specific cytokines (G-CSF, M-CSF, or EPO). In cultures with only lineage-unspecific cytokines (IL-3, SCF, GM-CSF), BAALC transcripts persisted up to day 20, while CD34 transcripts disappeared earlier. These observations suggest that expression of BAALC is stage specific. Conclusions BAALC expression is restricted to progenitor cells, and downregulation of BAALC occurs with cell differentiation. We postulate that BAALC represents a novel marker of an early progenitor cell common to the myeloid, lymphoid, and erythroid pathways.

TP53 does not fully comply with the Knudson model [Knudson, A. G., Jr. (1971) Proc. Natl. Acad. Sci. USA 68, 820–823] in that a reduction of constitutional expression of p53 may be sufficient for tumor predisposition . This finding suggests a gene-dosage effect for p53 function. To determine whether TP53 gene dosage affects the transcriptional regulation of target genes, we performed oligonucleotide-array gene expression analysis by using human cells with wild-type p53 (p53 +/+), or with one (p53 +/−), or both (p53 −/−) TP53 alleles disrupted by homologous recombination. We identified 35 genes whose expression is significantly correlated to the dosage of TP53 . These genes are involved in a variety of cellular processes including signal transduction, cell adhesion, and transcription regulation. Several of them are involved in neurogenesis and neural crest migration, developmental processes in which p53 is known to play a role. Motif search analysis revealed that of the genes highly expressed in p53 +/+ and +/− cells, several contain a putative p53 consensus binding site (bs), suggesting that they could be directly regulated by p53. Among those genes, we chose CSPG2 (which encodes versican) for further study because it contains a bona fide p53 bs in its first intron and its expression highly correlates with TP53 dosage. By using in vitro and in vivo assays, we showed CSPG2 to be directly transactivated by p53. In conclusion, we developed a strategy to demonstrate that many genes are affected by TP53 gene dosage for their expression. We report several candidate genes as potential downstream targets of p53 in nonstressed cells. Among them, CSPG2 is validated as being directly transactivated by p53. Our method provides a useful tool to elucidate additional mechanisms by which p53 exerts its functions.

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P, Author Summary Individuals with thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune thyroid diseases (AITD), which are common in the general population and associated with increased cardiovascular, metabolic and psychiatric morbidity and mortality. As the causative genes of TPOAbs and AITD remain largely unknown, we performed a genome-wide scan for TPOAbs in 18,297 individuals, with replication in 8,990 individuals. Significant associations were detected with variants at TPO, ATXN2, BACH2, MAGI3, and KALRN. Individuals carrying multiple risk variants also had a higher risk of increased thyroid-stimulating hormone levels (including subclinical and overt hypothyroidism), and a decreased risk of goiter. The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, and the MAGI3 variant was also associated with an increased risk of hypothyroidism. This first genome-wide scan for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. These results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which individuals are particularly at risk of developing clinical thyroid dysfunction.

Background: Over 15 genome-wide association studies (GWASs) have identified a narrow region in chromosome 6q23.3 as a major hematopoietic quantitative trait locus (QTL) for many hematologic traits. The associated phenotypes have yet to be fully explained by genes in the region, including MYB. We noted a newly annotated microRNA within the QTL, miR-3662, and sought to characterize its role in the QTL by exploring its function in normal and malignant hematopoiesis, specifically acute myeloid leukemia (AML). Results: To assess the association between the QTL-defining SNPs and miR-3662 abundance, we genotyped 10 GWAS SNPs in 200 blood samples from non-leukemic individuals. The risk alleles of all 10 SNPs associated with higher miR-3662 abundance. Homozygosity for rs66650371, a three base pair (bp) deletion, showed the strongest association (P To explore the role of miR-3662 in malignant hematopoiesis, we overexpressed miR-3662 in AML cell lines (Kasumi1, KG1a, MV4-11) and primary AML patient blasts (n = 12). Overexpression of miR-3662 reduced cell growth and viability in vitro, and significantly reduced tumor size in a murine xenograft model (scramble vs. miR-3662; n = 9 mice/group). We used targeted RNA sequencing to identify differentially expressed genes between forced miR-3662-expresssing- and control AML patient blasts using a panel of 361 candidate genes (TruSeq RNA platform, Illumina). IKBKB was identified as a promising target, and direct regulation of IKBKB's 3′-UTR by miR-3662 was confirmed via luciferase assay. MiR-3662 overexpression lowered IKBKB mRNA and protein expression in both HPCs and AML patient blasts. As IKBKB promotes activation and nuclear localization of NF-κB, overexpression of miR-3662 led to decreased NF-κB transactivation potential via luciferase assay. Confocal microscopy of HPCs and AML cells, as well as subcellular fraction immunoblots, revealed decreased NF-κB nuclear localization upon miR-3662 overexpression. Re-introduction of IKBKB partly rescued miR-3662's phenotype by increasing nuclear localization of NF-κB, decreasing growth of HPCs, and increasing viability of AML cells. Conclusion: The risk alleles of the 6q23.3 region associated with higher miR-3662 abundance, and miR-3662 abundance was functionally linked to rs66650371. We characterized miR-3662 as an agent in normal and malignant hematopoiesis, and demonstrated that miR-3662 reduces NF-κB activity by targeting IKBKB. Citation Format: Sophia E. Maharry, Christopher J. Walker, Sandya Liyanarachchi, Sujay Mehta, James S. Blachly, Denis C. Guttridge, Clara D. Bloomfield, Albert de la Chapelle, Ann-Kathrin Eisfeld. Dissection of the major hematopoietic quantitative trait locus in chromosome 6q23.3 identifies miR-3662 as a player in hematopoiesis and AML. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 1117.

A genome-wide association study of papillary thyroid carcinoma (PTC) pinpointed two independent SNPs (rs944289 and rs965513) located in regions containing no annotated genes (14q13.3 and 9q22.33, respectively). Here, we describe a unique, long, intergenic, noncoding RNA gene (lincRNA) named Papillary Thyroid Carcinoma Susceptibility Candidate 3 ( PTCSC3 ) located 3.2 kb downstream of rs944289 at 14q.13.3 and the expression of which is strictly thyroid specific. By quantitative PCR, PTCSC3 expression was strongly down-regulated ( P = 2.84 × 10 −14 ) in thyroid tumor tissue of 46 PTC patients and the risk allele (T) was associated with the strongest suppression (genotype [TT] ( n = 21) vs. [CT] ( n = 19), P = 0.004). In adjacent unaffected thyroid tissue, the genotype [TT] was associated with up-regulation of PTCSC3 ([TT] ( n = 21) vs. [CT] ( n = 19), P = 0.034). The SNP rs944289 was located in a binding site for the CCAAT/enhancer binding proteins (C/EBP) α and β. The risk allele destroyed the binding site in silico. Both C/EBPα and C/EBPβ activated the PTCSC3 promoter in reporter assays ( P = 0.0009 and P = 0.0014, respectively) and the risk allele reduced the activation compared with the nonrisk allele (C) ( P = 0.026 and P = 0.048, respectively). Restoration of PTCSC3 expression in PTC cell line cells (TPC-1 and BCPAP) inhibited cell growth ( P = 0.002 and P = 0.019, respectively) and affected the expression of genes involved in DNA replication, recombination and repair, cellular movement, tumor morphology, and cell death. Our data suggest that SNP rs944289 predisposes to PTC through a previously uncharacterized, long intergenic noncoding RNA gene ( PTCSC3 ) that has the characteristics of a tumor suppressor.

Glioblastoma multiforme (GBM), the most common and aggressive primary brain malignancy, is incurable despite the best combination of current cancer therapies. For the development of more effective therapies, discovery of novel candidate tumor drivers is urgently needed. Here, we report that peroxiredoxin 4 (PRDX4) is a putative tumor driver. PRDX4 levels were highly increased in a majority of human GBMs as well as in a mouse model of GBM. Reducing PRDX4 expression significantly decreased GBM cell growth and radiation resistance in vitro with increased levels of ROS, DNA damage, and apoptosis. In a syngenic orthotopic transplantation model, Prdx4 knockdown limited GBM infiltration and significantly prolonged mouse survival. These data suggest that PRDX4 can be a novel target for GBM therapies in the future.

International audience; BACKGROUND & AIMS: Expression of the netrin-1 dependence receptor UNC5C is reduced in many colorectal tumors; mice with the UNC5C mutations have increased progression of intestinal tumors. We investigated whether specific variants in UNC5C increase risk of colorectal cancer (CRC). METHODS: We analyzed the sequence of UNC5C in blood samples from 1801 patients with CRC and 4152 controls from 3 cohorts (France, United States, and Finland). Almost all cases from France and the United States had familial CRC; of the Finnish cases, 92 of 984 were familial. We analyzed whether CRC segregates with the UNC5C variant A628K in 3 families with histories of CRC. We also performed haplotype analysis to determine the origin of this variant. RESULTS: Of 817 patients with familial CRC, 14 had 1 of 4 different, unreported missense variants in UNC5C. The variants p.Asp353Asn (encodes D353N), p.Arg603Cys (encodes R603C), and p.Gln630Glu (encodes Q630E) did not occur significantly more often in cases than controls. The variant p.Ala628Lys (A628K) was detected in 3 families in the French cohort (odds ratio, 8.8; Wald's 95% confidence interval, 1.47-52.93; P = .03) and in 2 families in the US cohort (odds ratio, 1.9; P = .6) but was not detected in the Finnish cohort; UNC5C A628K segregated with CRC in families. Three families with A628K had a 109-kilobase identical haplotype that spanned most of UNC5C, indicating recent origin of this variant in white subjects (14 generations; 95% confidence interval, 6-36 generations). Transfection of HEK293T cells with UNC5C-A628K significantly reduced apoptosis compared with wild-type UNC5C, measured in an assay of active caspase-3. CONCLUSIONS: Inherited mutations in UNC5C prevent apoptosis and increase risk of CRC.