Search

Supplementary Table 1. Linkage disequilibrium patterns of the tag-SNPs used for genotyping the 6q23.3 locus. Supplementary Table 2. Cytogenetic and molecular information on AML patients used for functional studies (n=12). Supplementary Table 3. French-American-British (FAB) classification35 and molecular information on AML patients used to determine miR-3662's abundance (n=8). Supplementary Table 4. Canonical pathway analysis of the miR-3662-associated gene expression signature. Supplementary Figure 1. Transcription factor binding according to the transcription factor chip data from ENCODE. Supplementary Figure 2. Electrophoretic mobility shift assay comparing the binding affinity of the alleles of rs66650371 and rs9483788. Supplementary Figure 3. Endogenous miR-3662 expression levels of hematopoietic progenitor (HP) cells during differentiation, total bone marrow aspirate of three non-leukemic donors (total BM 1-3), different populations of differentiated peripheral blood cells, and three AML cell lines. Supplementary Figure 4. Top panel, Macroscopic pictures of the spleens of three mice of the scramble and miR-3662-infected groups (organs harvested post-mortem). All mice had a massive splenomegaly compared to the un-injected, sacrificed control mouse. Bottom panel, images of spleen histologies (40x enlargement). Slides were stained for CD45 to proof MV4-11 origin of the leukemia. Supplementary Figure 5. Endogenous abundance of miR-3662 and IKBKB in patient samples and cell lines. Supplementary Figure 6. Comparison of the relative miR-3662 abundance of AML patient blasts and AML cell lines before (black) and after (red) forced miR-3662 expression with the lentiviral expression construct.

Supplementary Methods, Supplementary References Supplementary Table S1: Inherited regions of homozygosity in >1% of 1,798 non-leukemic individuals Supplementary Table S2: Acquired uniparental disomies detected in 425 cytogenetically normal AML patients Supplementary Table S3. Associations between recurrent UPDs with pretreatment patient characteristics for patients with cytogenetically normal acute myeloid leukemia Supplementary Table S5. Variant allele fraction of gene mutations that co-occurred with UPDs Supplementary Table S6. Allelic ratio of FLT3 internal tandem duplications that co-occurred with UPDs Supplementary Table S7: Copy number gains and losses in 425 patients with cytogenetically normal acute myeloid leukemia

Purpose:Uniparental disomy (UPD) is a way cancer cells duplicate a mutated gene, causing loss of heterozygosity (LOH). Patients with cytogenetically normal acute myeloid leukemia (CN-AML) do not have microscopically detectable chromosome abnormalities, but can harbor UPDs. We examined the prognostic significance of UPDs and frequency of LOH in patients with CN-AML.Experimental Design: We examined the frequency and prognostic significance of UPDs in a set of 425 adult patients with de novo CN-AML who were previously sequenced for 81 genes typically mutated in cancer. Associations of UPDs with outcome were analyzed in the 315 patients with CN-AML younger than 60 years.Results:We detected 127 UPDs in 109 patients. Most UPDs were large and typically encompassed all or most of the affected chromosome arm. The most common UPDs occurred on chromosome arms 13q (7.5% of patients), 6p (2.8%), and 11p (2.8%). Many UPDs significantly cooccurred with mutations in genes they encompassed, including 13q UPD with FLT3-internal tandem duplication (FLT3-ITD; P < 0.001), and 11p UPD with WT1 mutations (P = 0.02). Among patients younger than 60 years, UPD of 11p was associated with longer overall survival (OS) and 13q UPD with shorter disease-free survival (DFS) and OS. In multivariable models that accounted for known prognostic markers, including FLT3-ITD and WT1 mutations, UPD of 13q maintained association with shorter DFS, and UPD of 11p maintained association with longer OS.Conclusions:LOH mediated by UPD is a recurrent feature of CN-AML. Detection of UPDs of 13q and 11p might be useful for genetic risk stratification of patients with CN-AML.

Mutations detected in 425 CN-AML patients

Purpose: Aberrant DNA methylation, now recognized as a contributing factor to neoplasia, often shows definitive gene/sequence preferences unique to specific cancer types. Correspondingly, distinct combinations of methylated loci can function as biomarkers for numerous clinical correlates of ovarian and other cancers.Experimental Design: We used a microarray approach to identify methylated loci prognostic for reduced progression-free survival (PFS) in advanced ovarian cancer patients. Two data set classification algorithms, Significance Analysis of Microarray and Prediction Analysis of Microarray, successfully identified 220 candidate PFS-discriminatory methylated loci. Of those, 112 were found capable of predicting PFS with 95% accuracy, by Prediction Analysis of Microarray, using an independent set of 40 advanced ovarian tumors (from 20 short-PFS and 20 long-PFS patients, respectively). Additionally, we showed the use of these predictive loci using two bioinformatics machine-learning algorithms, Support Vector Machine and Multilayer Perceptron.Conclusion: In this report, we show that highly prognostic DNA methylation biomarkers can be successfully identified and characterized, using previously unused, rigorous classifying algorithms. Such ovarian cancer biomarkers represent a promising approach for the assessment and management of this devastating disease.

Supplementary Table 1 from Origins and Prevalence of the American Founder Mutation of MSH2

Supplementary Table 2 from A Susceptibility Locus for Papillary Thyroid Carcinoma on Chromosome 8q24

Supplementary Table 2 from Xenoestrogen-Induced Epigenetic Repression of microRNA-9-3 in Breast Epithelial Cells

Supplementary Figure 1 from A Susceptibility Locus for Papillary Thyroid Carcinoma on Chromosome 8q24

Supplementary Figures 1-6, Tables 1-15, Methods from Epigenetic Silencing Mediated through Activated PI3K/AKT Signaling in Breast Cancer

Supplementary Figure 2 from A Susceptibility Locus for Papillary Thyroid Carcinoma on Chromosome 8q24

Supplementary Table 1 from A Susceptibility Locus for Papillary Thyroid Carcinoma on Chromosome 8q24

Supplementary Figure 1 from Origins and Prevalence of the American Founder Mutation of MSH2

Supplementary Table 1 from Xenoestrogen-Induced Epigenetic Repression of microRNA-9-3 in Breast Epithelial Cells

The interplay between histone modifications and promoter hypermethylation provides a causative explanation for epigenetic gene silencing in cancer. Less is known about the upstream initiators that direct this process. Here, we report that the Cystatin M (CST6) tumor suppressor gene is concurrently down-regulated with other loci in breast epithelial cells cocultured with cancer-associated fibroblasts (CAF). Promoter hypermethylation of CST6 is associated with aberrant AKT1 activation in epithelial cells, as well as the disabled INNP4B regulator resulting from the suppression by CAFs. Repressive chromatin, marked by trimethyl-H3K27 and dimethyl-H3K9, and de novo DNA methylation is established at the promoter. The findings suggest that microenvironmental stimuli are triggers in this epigenetic cascade, leading to the long-term silencing of CST6 in breast tumors. Our present findings implicate a causal mechanism defining how tumor stromal fibroblasts support neoplastic progression by manipulating the epigenome of mammary epithelial cells. The result also highlights the importance of direct cell-cell contact between epithelial cells and the surrounding fibroblasts that confer this epigenetic perturbation. Because this two-way interaction is anticipated, the described coculture system can be used to determine the effect of epithelial factors on fibroblasts in future studies. [Cancer Res 2008;68(24):10257–66]

Supplementary Table 2 from Origins and Prevalence of the American Founder Mutation of MSH2

Supplementary Table 3 from A Susceptibility Locus for Papillary Thyroid Carcinoma on Chromosome 8q24

Supplementary Methods, Figures 1-5, Tables 1-5 from Breast Cancer–Associated Fibroblasts Confer AKT1-Mediated Epigenetic Silencing of Cystatin M in Epithelial Cells

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]

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.

Progressive metastasis is the proximate cause of cancer-related mortality for the majority of patients with solid tumors. Identifying both drivers, and gatekeepers of late stage progression (metastasis progression suppressors; MPS) is crucial to improve biomarkers and treatments. We previously identified Regulator of Calcineurin 1.4 (RCAN1.4) as a MPS gene in thyroid cancer that functions in part through overexpression of the cap-n-collar transcription factor NFE2L3 (NRF3). However, the mechanisms for RCAN1.4 loss in thyroid cancer, and whether or not this is reversible, had not be addressed. Because the RCAN 1.4 promoter contains CG rich regions reported to be hypermethyalated in other tissues, we hypothesized that RCAN 1.4 downregulation in thyroid cancer was similarly regulated. Studies were performed in both cell lines and de-identified human papillary thyroid cancer (PTC) tumor and normal tissue pairs from our Endocrine Neoplasia Repository tumor bank. Five thyroid cancer cell lines (TPC1, FTC133, BCPAP, C643, and 8505C) with different genetic drivers were selected based on low expression of RCAN 1.4. To determine if lines display methylation-regulated expression of RCAN 1.4, cells were treated with 1-10uM decitabine, a DNA methyl transferase inhibitor. In all cell lines, RCAN1.4 mRNA and protein levels increased with decitabine treatment, consistent with hypermethylation. Analysis of the literature and in silico analysis of the RCAN 1.4 gene identified CG rich regions as putative methylation targets in the proximal promoter and in intron 1. Using two different methods; 1) isolation of MseI-digestion DNA fragments using CpG MethylQuest beads followed by qPCR and 2) methylation-specific PCR (MSP) of bisulfite-modified DNA, hypermethylation of the intron 1 CG rich region was identified. There was no evidence of hypermethylation of the CG rich region of the proximal promoter. Methylation of the two regions also was assessed on 18 pairs PTC and adjacent normal tissues. Similar to the cell lines, hypermethylation was identified at intron 1 of RCAN 1.4 in PTC versus normal tissues. QRT-PCR results confirmed lower levels of RCAN1-4 mRNA in the PTC vs normal tissue. These samples also demonstrated increased NRF3 expression versus normal tissue. Taken together, these data demonstrate that RCAN 1.4 is downregulated in thyroid cancer cells and human PTCs by hypermethylation of CG rich regions in intron 1 and that this effect is reversible using demethylating agents. Presentation: Saturday, June 11, 2022 11:45 a.m. - 12:00 p.m.

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.

Purpose:Uniparental disomy (UPD) is a way cancer cells duplicate a mutated gene, causing loss of heterozygosity (LOH). Patients with cytogenetically normal acute myeloid leukemia (CN-AML) do not have microscopically detectable chromosome abnormalities, but can harbor UPDs. We examined the prognostic significance of UPDs and frequency of LOH in patients with CN-AML.Experimental Design: We examined the frequency and prognostic significance of UPDs in a set of 425 adult patients with de novo CN-AML who were previously sequenced for 81 genes typically mutated in cancer. Associations of UPDs with outcome were analyzed in the 315 patients with CN-AML younger than 60 years.Results:We detected 127 UPDs in 109 patients. Most UPDs were large and typically encompassed all or most of the affected chromosome arm. The most common UPDs occurred on chromosome arms 13q (7.5% of patients), 6p (2.8%), and 11p (2.8%). Many UPDs significantly cooccurred with mutations in genes they encompassed, including 13q UPD with FLT3-internal tandem duplication (FLT3-ITD; P < 0.001), and 11p UPD with WT1 mutations (P = 0.02). Among patients younger than 60 years, UPD of 11p was associated with longer overall survival (OS) and 13q UPD with shorter disease-free survival (DFS) and OS. In multivariable models that accounted for known prognostic markers, including FLT3-ITD and WT1 mutations, UPD of 13q maintained association with shorter DFS, and UPD of 11p maintained association with longer OS.Conclusions:LOH mediated by UPD is a recurrent feature of CN-AML. Detection of UPDs of 13q and 11p might be useful for genetic risk stratification of patients with CN-AML.

Background: Thyroid cancer patients with radioiodine-refractory (RAI-R) disease, resulting from insufficient RAI delivery and/or RAI resistance, have increased mortality and limited treatment options. To date, studies have largely focused on tumor mutations associated with different stages of disease, which could provide prognostic value for RAI-R disease. It was hypothesized that germline variants contributing to intrinsic differences in iodine metabolism, tumor microenvironment, and/or immune surveillance are associated with RAI-R disease. Methods: Whole-genome genotyping data analysis was performed on 1145 Caucasian (CAU) patients, 244 of whom were RAI-R, and 55 African American (AA) patients, nine of whom were RAI-R. Germline-variant association studies were conducted using candidate genes involved in iodine metabolism or DNA-damage repair, as well as genome-wide association analysis. Initial data indicated several notable variants in a small number of patients (n = 7), who were later determined to be AA patients of >80% African ancestry (n = 37). This led to the study focusing on germline single nucleotide polymorphisms uniquely associated with RAI-R AA patients. Sanger sequencing was performed to validate risk alleles and identify the incidence of the common somatic mutations BRAF(V600E), NRAS(Q61R), and HRAS(Q61R) in AA patients whose primary tumor samples were available (28/55). Results: TG, BRCA1, and NSMCE2 haplotypes were identified as being uniquely associated with RAI-R AA patients of >80% African ancestry. All patients with the TG haplotype (n = 4) had a biochemical incomplete response to RAI therapy. Patients with the NSMCE2 haplotype (n = 4) were diagnosed at a young age (13, 17, 17, and 26 years old) with distant metastatic disease at initial diagnosis. The BRCA1 haplotype co-occurred in three out of four patients with the NSMCE2 haplotype. The incidence of BRAF(V600E) appears lower in papillary thyroid carcinomas from AA patients of >80% African ancestry (3/14; 21%) than in AA patients of

BackgroundPatients with colorectal cancer (CRC) with mismatch repair-deficient (dMMR) tumours without MLH1 methylation or germline MMR pathogenic variants (PV) were previously thought to have Lynch syndrome (LS). It is now appreciated that they can have double somatic (DS) MMR PVs. We explored the clinical characteristics between patients with DS tumours and LS in two population-based cohorts.MethodsWe included patients with CRC from Ohio 2013–2016 and Iceland 2000–2009. All had microsatellite instability testing and/or immunohistochemistry (IHC) of MMR proteins, and MLH1 methylation testing when indicated. Germline next-generation sequencing was performed for all with dMMR tumours; tumour sequencing followed for patients with unexplained dMMR. Clinical characteristics of DS patients and patients with LS were compared.ResultsOf the 232 and 51 patients with non-methylated dMMR tumours in the Ohio and Iceland cohorts, respectively, 57.8% (n=134) and 45.1% (n=23) had LS, 32.8% (n=76) and 31.4% (n=16) had DS PVs, 6% (n=14) and 9.8% (n=5) were unexplained and 4.3% (n=10) and 13.7% (n=7) had incorrect IHC. Age of diagnosis for DS patients was older than patients with LS (p=3.73×10−4) in the two cohorts. Patients with LS were more likely to meet Amsterdam II criteria (OR=15.81, p=8.47×10−6) and have multiple LS-associated tumours (OR=6.67, p=3.31×10−5). Absence of MLH1/PMS2 was predictive of DS PVs; isolated MSH6 and PMS2 absence was predictive of LS in both cohorts.ConclusionsIndividuals with LS are 15× more likely to meet Amsterdam II criteria and >5× more likely to have multiple cancers as compared with those with DS tumours. Furthermore, isolated loss of MSH6 or PMS2 protein predicts LS.

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: 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.

Genome-wide association studies (GWASs) have identified at least 10 single-nucleotide polymorphisms (SNPs) associated with papillary thyroid cancer (PTC) risk. Most of these SNPs are common variants with small to moderate effect sizes. Here we assessed the combined genetic effects of these variants on PTC risk by using summarized GWAS results to build polygenic risk score (PRS) models in three PTC study groups from Ohio (1,544 patients and 1,593 controls), Iceland (723 patients and 129,556 controls), and the United Kingdom (534 patients and 407,945 controls). A PRS based on the 10 established PTC SNPs showed a stronger predictive power compared with the clinical factors model, with a minimum increase of area under the receiver-operating curve of 5.4 percentage points ( P ≤ 1.0 × 10 −9 ). Adding an extended PRS based on 592,475 common variants did not significantly improve the prediction power compared with the 10-SNP model, suggesting that most of the remaining undiscovered genetic risk in thyroid cancer is due to rare, moderate- to high-penetrance variants rather than to common low-penetrance variants. Based on the 10-SNP PRS, individuals in the top decile group of PRSs have a close to sevenfold greater risk (95% CI, 5.4–8.8) compared with the bottom decile group. In conclusion, PRSs based on a small number of common germline variants emphasize the importance of heritable low-penetrance markers in PTC.