Inherited myeloproliferative neoplasm risk affects haematopoietic stem cells.

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P < 5.0 × 10-8), 7 of which have not been previously reported. We find that there is a shared genetic architecture between MPN risk and several haematopoietic traits from distinct lineages; that there is an enrichment for MPN risk variants within accessible chromatin of HSCs; and that increased MPN risk is associated with longer telomere length in leukocytes and other clonal haematopoietic states-collectively suggesting that MPN risk is associated with the function and self-renewal of HSCs. We use gene mapping to identify modulators of HSC biology linked to MPN risk, and show through targeted variant-to-function assays that CHEK2 and GFI1B have roles in altering the function of HSCs to confer disease risk. Overall, our results reveal a previously unappreciated mechanism for inherited MPN risk through the modulation of HSC function.

Genome-wide Modeling of Polygenic Risk Score in Colorectal Cancer Risk.

Accurate colorectal cancer (CRC) risk prediction models are critical for identifying individuals at low and high risk of developing CRC, as they can then be offered targeted screening and interventions to address their risks of developing disease (if they are in a high-risk group) and avoid unnecessary screening and interventions (if they are in a low-risk group). As it is likely that thousands of genetic variants contribute to CRC risk, it is clinically important to investigate whether these genetic variants can be used jointly for CRC risk prediction. In this paper, we derived and compared different approaches to generating predictive polygenic risk scores (PRS) from genome-wide association studies (GWASs) including 55,105 CRC-affected case subjects and 65,079 control subjects of European ancestry. We built the PRS in three ways, using (1) 140 previously identified and validated CRC loci; (2) SNP selection based on linkage disequilibrium (LD) clumping followed by machine-learning approaches; and (3) LDpred, a Bayesian approach for genome-wide risk prediction. We tested the PRS in an independent cohort of 101,987 individuals with 1,699 CRC-affected case subjects. The discriminatory accuracy, calculated by the age- and sex-adjusted area under the receiver operating characteristics curve (AUC), was highest for the LDpred-derived PRS (AUC = 0.654) including nearly 1.2 M genetic variants (the proportion of causal genetic variants for CRC assumed to be 0.003), whereas the PRS of the 140 known variants identified from GWASs had the lowest AUC (AUC = 0.629). Based on the LDpred-derived PRS, we are able to identify 30% of individuals without a family history as having risk for CRC similar to those with a family history of CRC, whereas the PRS based on known GWAS variants identified only top 10% as having a similar relative risk. About 90% of these individuals have no family history and would have been considered average risk under current screening guidelines, but might benefit from earlier screening. The developed PRS offers a way for risk-stratified CRC screening and other targeted interventions.

Assessing thyroid cancer risk using polygenic risk scores.

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.

Mutational landscape and clinical outcome of patients with de novo acute myeloid leukemia and rearrangements involving 11q23/KMT2A.

Balanced rearrangements involving the KMT2A gene, located at 11q23, are among the most frequent chromosome aberrations in acute myeloid leukemia (AML). Because of numerous fusion partners, the mutational landscape and prognostic impact of specific 11q23/KMT2A rearrangements are not fully understood. We analyzed clinical features of 172 adults with AML and recurrent 11q23/KMT2A rearrangements, 141 of whom had outcome data available. We compared outcomes of these patients with outcomes of 1,097 patients without an 11q23/KMT2A rearrangement categorized according to the 2017 European LeukemiaNet (ELN) classification. Using targeted next-generation sequencing, we investigated the mutational status of 81 leukemia/cancer-associated genes in 96 patients with 11q23/KMT2A rearrangements with material for molecular studies available. Patients with 11q23/KMT2A rearrangements had a low number of additional gene mutations (median, 1; range 0 to 6), which involved the RAS pathway (KRAS, NRAS, and PTPN11) in 32% of patients. KRAS mutations occurred more often in patients with t(6;11)(q27;q23)/KMT2A-AFDN compared with patients with the other 11q23/KMT2A subsets. Specific gene mutations were too infrequent in patients with specific 11q23/KMT2A rearrangements to assess their associations with outcomes. We demonstrate that younger (age <60 y) patients with t(9;11)(p22;q23)/KMT2A-MLLT3 had better outcomes than patients with other 11q23/KMT2A rearrangements and those without 11q23/KMT2A rearrangements classified in the 2017 ELN intermediate-risk group. Conversely, outcomes of older patients (age ≥60 y) with t(9;11)(p22;q23) were poor and comparable to those of the ELN adverse-risk group patients. Our study shows that patients with an 11q23/KMT2A rearrangement have distinct mutational patterns and outcomes depending on the fusion partner.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium.

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We 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 and 16q23.2, which were associated with thyroid-stimulating hormone 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. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Circulating Levels of Insulin-like Growth Factor 1 and Insulin-like Growth Factor Binding Protein 3 Associate With Risk of Colorectal Cancer Based on Serologic and Mendelian Randomization Analyses.

BACKGROUND & AIMS: Human studies examining associations between circulating levels of insulin-like growth factor 1 (IGF1) and insulin-like growth factor binding protein 3 (IGFBP3) and colorectal cancer risk have reported inconsistent results. We conducted complementary serologic and Mendelian randomization (MR) analyses to determine whether alterations in circulating levels of IGF1 or IGFBP3 are associated with colorectal cancer development. METHODS: Serum levels of IGF1 were measured in blood samples collected from 397,380 participants from the UK Biobank, from 2006 through 2010. Incident cancer cases and cancer cases recorded first in death certificates were identified through linkage to national cancer and death registries. Complete follow-up was available through March 31, 2016. For the MR analyses, we identified genetic variants associated with circulating levels of IGF1 and IGFBP3. The association of these genetic variants with colorectal cancer was examined with 2-sample MR methods using genome-wide association study consortia data (52,865 cases with colorectal cancer and 46,287 individuals without [controls]) RESULTS: After a median follow-up period of 7.1 years, 2665 cases of colorectal cancer were recorded. In a multivariable-adjusted model, circulating level of IGF1 associated with colorectal cancer risk (hazard ratio per 1 standard deviation increment of IGF1, 1.11; 95% confidence interval [CI] 1.05-1.17). Similar associations were found by sex, follow-up time, and tumor subsite. In the MR analyses, a 1 standard deviation increment in IGF1 level, predicted based on genetic factors, was associated with a higher risk of colorectal cancer risk (odds ratio 1.08; 95% CI 1.03-1.12; P = 3.3 × 10-4). Level of IGFBP3, predicted based on genetic factors, was associated with colorectal cancer risk (odds ratio per 1 standard deviation increment, 1.12; 95% CI 1.06-1.18; P = 4.2 × 10-5). Colorectal cancer risk was associated with only 1 variant in the IGFBP3 gene region (rs11977526), which also associated with anthropometric traits and circulating level of IGF2. CONCLUSIONS: In an analysis of blood samples from almost 400,000 participants in the UK Biobank, we found an association between circulating level of IGF1 and colorectal cancer. Using genetic data from 52,865 cases with colorectal cancer and 46,287 controls, a higher level of IGF1, determined by genetic factors, was associated with colorectal cancer. Further studies are needed to determine how this signaling pathway might contribute to colorectal carcinogenesis.

Double somatic mismatch repair gene pathogenic variants as common as Lynch syndrome among endometrial cancer patients.

OBJECTIVE: Lynch syndrome is the most common cause of inherited endometrial cancer, attributable to germline pathogenic variants (PV) in mismatch repair (MMR) genes. Tumor microsatellite instability (MSI-high) and MMR IHC abnormalities are characteristics of Lynch syndrome. Double somatic MMR gene PV also cause MSI-high endometrial cancers. The aim of this study was to determine the relative frequency of Lynch syndrome and double somatic MMR PV. METHODS: 341 endometrial cancer patients enrolled in the Ohio Colorectal Cancer Prevention Initiative at The Ohio State University Comprehensive Cancer Center from 1/1/13-12/31/16. All tumors underwent immunohistochemical (IHC) staining for the four MMR proteins, MSI testing, and MLH1 methylation testing if the tumor was MMR-deficient (dMMR). Germline genetic testing for Lynch syndrome was undertaken for all cases with dMMR tumors lacking MLH1 methylation. Tumor sequencing followed if a germline MMR gene PV was not identified. RESULTS: Twenty-seven percent (91/341) of tumors were either MSI-high or had abnormal IHC indicating dMMR. As expected, most dMMR tumors had MLH1 methylation; (69, 75.8% of the dMMR cases; 20.2% of total). Among the 22 (6.5%) cases with dMMR not explained by methylation, 10 (2.9% of total) were found to have Lynch syndrome (6 MSH6, 3 MSH2, 1 PMS2). Double somatic MMR PV accounted for the remaining 12 dMMR cases (3.5% of total). CONCLUSIONS: Since double somatic MMR gene PV are as common as Lynch syndrome among endometrial cancer patients, paired tumor and germline testing for patients with non-methylated dMMR tumor may be the most efficient approach for LS screening.

Variants in LRRC34 reveal distinct mechanisms for predisposition to papillary thyroid carcinoma.

BACKGROUND: Papillary 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). METHODS: We report the mechanisms of PTC risk at 3q26 using a combination of overexpression, mass spectroscopy, knockdown, transcriptome profiling, migration assays and genetic analysis. RESULTS: We 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. CONCLUSIONS: Our 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.

Adiposity, metabolites, and colorectal cancer risk: Mendelian randomization study.

BACKGROUND: Higher adiposity increases the risk of colorectal cancer (CRC), but whether this relationship varies by anatomical sub-site or by sex is unclear. Further, the metabolic alterations mediating the effects of adiposity on CRC are not fully understood. METHODS: We examined sex- and site-specific associations of adiposity with CRC risk and whether adiposity-associated metabolites explain the associations of adiposity with CRC. Genetic variants from genome-wide association studies of body mass index (BMI) and waist-to-hip ratio (WHR, unadjusted for BMI; N = 806,810), and 123 metabolites from targeted nuclear magnetic resonance metabolomics (N = 24,925), were used as instruments. Sex-combined and sex-specific Mendelian randomization (MR) was conducted for BMI and WHR with CRC risk (58,221 cases and 67,694 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium, Colorectal Cancer Transdisciplinary Study, and Colon Cancer Family Registry). Sex-combined MR was conducted for BMI and WHR with metabolites, for metabolites with CRC, and for BMI and WHR with CRC adjusted for metabolite classes in multivariable models. RESULTS: In sex-specific MR analyses, higher BMI (per 4.2 kg/m2) was associated with 1.23 (95% confidence interval (CI) = 1.08, 1.38) times higher CRC odds among men (inverse-variance-weighted (IVW) model); among women, higher BMI (per 5.2 kg/m2) was associated with 1.09 (95% CI = 0.97, 1.22) times higher CRC odds. WHR (per 0.07 higher) was more strongly associated with CRC risk among women (IVW OR = 1.25, 95% CI = 1.08, 1.43) than men (IVW OR = 1.05, 95% CI = 0.81, 1.36). BMI or WHR was associated with 104/123 metabolites at false discovery rate-corrected P ≤ 0.05; several metabolites were associated with CRC, but not in directions that were consistent with the mediation of positive adiposity-CRC relations. In multivariable MR analyses, associations of BMI and WHR with CRC were not attenuated following adjustment for representative metabolite classes, e.g., the univariable IVW OR for BMI with CRC was 1.12 (95% CI = 1.00, 1.26), and this became 1.11 (95% CI = 0.99, 1.26) when adjusting for cholesterol in low-density lipoprotein particles. CONCLUSIONS: Our results suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. Adiposity was associated with numerous metabolic alterations, but none of these explained associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.

Clinical characteristics of patients with colorectal cancer with double somatic mismatch repair mutations compared with Lynch syndrome.

BACKGROUND: Patients 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. METHODS: We 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. RESULTS: Of 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. CONCLUSIONS: Individuals 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.

MYH9 binds to lncRNA gene PTCSC2 and regulates FOXE1 in the 9q22 thyroid cancer risk locus.

A locus on chromosome 9q22 harbors a SNP (rs965513) firmly associated with risk of papillary thyroid carcinoma (PTC). The locus also comprises the forkhead box E1 (FOXE1) gene, which is implicated in thyroid development, and a long noncoding RNA (lncRNA) gene, papillary thyroid cancer susceptibility candidate 2 (PTCSC2). How these might interact is not known. Here we report that PTCSC2 binds myosin-9 (MYH9). In a bidirectional promoter shared by FOXE1 and PTCSC2, MYH9 inhibits the promoter activity in both directions. This inhibition can be reversed by PTCSC2, which acts as a suppressor. RNA knockdown of FOXE1 in primary thyroid cells profoundly interferes with the p53 pathway. We propose that the interaction between the lncRNA, its binding protein MYH9, and the coding gene FOXE1 underlies the predisposition to PTC triggered by rs965513.

Identification of NRAS isoform 2 overexpression as a mechanism facilitating BRAF inhibitor resistance in malignant melanoma.

Activating mutations in BRAF are found in 50% of melanomas and although treatment with BRAF inhibitors (BRAFi) is effective, resistance often develops. We now show that recently discovered NRAS isoform 2 is up-regulated in the setting of BRAF inhibitor resistance in melanoma, in both cell lines and patient tumor tissues. When isoform 2 was overexpressed in BRAF mutant melanoma cell lines, melanoma cell proliferation and in vivo tumor growth were significantly increased in the presence of BRAFi treatment. shRNA-mediated knockdown of isoform 2 in BRAFi resistant cells restored sensitivity to BRAFi compared with controls. Signaling analysis indicated decreased mitogen-activated protein kinase (MAPK) pathway signaling and increased phosphoinositol-3-kinase (PI3K) pathway signaling in isoform 2 overexpressing cells compared with isoform 1 overexpressing cells. Immunoprecipitation of isoform 2 validated a binding affinity of this isoform to both PI3K and BRAF/RAF1. The addition of an AKT inhibitor to BRAFi treatment resulted in a partial restoration of BRAFi sensitivity in cells expressing high levels of isoform 2. NRAS isoform 2 may contribute to resistance to BRAFi by facilitating PI3K pathway activation.

Fine mapping of 14q13 reveals novel variants associated with different histological subtypes of papillary thyroid carcinoma.

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.

PD-1 Blockade in Tumors with Mismatch-Repair Deficiency.

BACKGROUND: Somatic mutations have the potential to encode "non-self" immunogenic antigens. We hypothesized that tumors with a large number of somatic mutations due to mismatch-repair defects may be susceptible to immune checkpoint blockade. METHODS: We conducted a phase 2 study to evaluate the clinical activity of pembrolizumab, an anti-programmed death 1 immune checkpoint inhibitor, in 41 patients with progressive metastatic carcinoma with or without mismatch-repair deficiency. Pembrolizumab was administered intravenously at a dose of 10 mg per kilogram of body weight every 14 days in patients with mismatch repair-deficient colorectal cancers, patients with mismatch repair-proficient colorectal cancers, and patients with mismatch repair-deficient cancers that were not colorectal. The coprimary end points were the immune-related objective response rate and the 20-week immune-related progression-free survival rate. RESULTS: The immune-related objective response rate and immune-related progression-free survival rate were 40% (4 of 10 patients) and 78% (7 of 9 patients), respectively, for mismatch repair-deficient colorectal cancers and 0% (0 of 18 patients) and 11% (2 of 18 patients) for mismatch repair-proficient colorectal cancers. The median progression-free survival and overall survival were not reached in the cohort with mismatch repair-deficient colorectal cancer but were 2.2 and 5.0 months, respectively, in the cohort with mismatch repair-proficient colorectal cancer (hazard ratio for disease progression or death, 0.10 [P<0.001], and hazard ratio for death, 0.22 [P=0.05]). Patients with mismatch repair-deficient noncolorectal cancer had responses similar to those of patients with mismatch repair-deficient colorectal cancer (immune-related objective response rate, 71% [5 of 7 patients]; immune-related progression-free survival rate, 67% [4 of 6 patients]). Whole-exome sequencing revealed a mean of 1782 somatic mutations per tumor in mismatch repair-deficient tumors, as compared with 73 in mismatch repair-proficient tumors (P=0.007), and high somatic mutation loads were associated with prolonged progression-free survival (P=0.02). CONCLUSIONS: This study showed that mismatch-repair status predicted clinical benefit of immune checkpoint blockade with pembrolizumab. (Funded by Johns Hopkins University and others; ClinicalTrials.gov number, NCT01876511.).

The role of SMAD3 in the genetic predisposition to papillary thyroid carcinoma.

PURPOSE: To identify and characterize the functional variants, regulatory gene networks, and potential binding targets of SMAD3 in the 15q22 thyroid cancer risk locus. METHODS: 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. RESULTS: 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. CONCLUSION: Our data provide a functional annotation of the 15q22 thyroid cancer risk locus.

Multiple functional variants in long-range enhancer elements contribute to the risk of SNP rs965513 in thyroid cancer.

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.

MicroRNA-3151 inactivates TP53 in BRAF-mutated human malignancies.

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.

Immunoglobulin transcript sequence and somatic hypermutation computation from unselected RNA-seq reads in chronic lymphocytic leukemia.

Immunoglobulins (Ig) are produced by B lymphocytes as secreted antibodies or as part of the B-cell receptor. There is tremendous diversity of potential Ig transcripts (>1 × 10(12)) as a result of hundreds of germ-line gene segments, random nucleotide incorporation during joining of gene segments into a complete transcript, and the process of somatic hypermutation at individual nucleotides. This recombination and mutation process takes place in the maturing B cell and is responsible for the diversity of potential epitope recognition. Cancers arising from mature B cells are characterized by clonal production of Ig heavy (IGH@) and light chain transcripts, although whether the sequence has undergone somatic hypermutation is dependent on the maturation stage at which the neoplastic clone arose. Chronic lymphocytic leukemia (CLL) is the most common leukemia in adults and arises from a mature B cell with either mutated or unmutated IGH@ transcripts, the latter having worse prognosis and the assessment of which is routinely performed in the clinic. Currently, IGHV mutation status is assessed by Sanger sequencing and comparing the transcript to known germ-line genes. In this paper, we demonstrate that complete IGH@ V-D-J sequences can be computed from unselected RNA-seq reads with results equal or superior to the clinical procedure: in the only discordant case, the clinical transcript was out-of-frame. Therefore, a single RNA-seq assay can simultaneously yield gene expression profile, SNP and mutation information, as well as IGHV mutation status, and may one day be performed as a general test to capture multidimensional clinically relevant data in CLL.

Frequent PIK3CA Mutations in Colorectal and Endometrial Tumors With 2 or More Somatic Mutations in Mismatch Repair Genes.

BACKGROUND & AIMS: Some colorectal and endometrial tumors with microsatellite instability not attributable to MLH1 hypermethylation or germline mutations contain 2 or more somatic mutations in genes encoding mismatch repair (MMR) proteins. We sought to define the molecular phenotype of this newly recognized tumor subtype. METHODS: From 2 prospective studies of the efficacy of screening for Lynch syndrome, we identified patients with colorectal and endometrial tumors who had 2 or more somatic (but not germline) mutations in genes encoding MMR proteins (double somatic). We determined the frequencies of tumor mutations in PIK3CA, BRAF, KRAS, NRAS, and PTEN by targeted next-generation sequencing and used logistic-regression models to compare them with those from patients with Lynch syndrome, MLH1-hypermethylated, or microsatellite-stable tumors. We validated our findings using independent data sets from The Cancer Genome Atlas. RESULTS: Among colorectal cancer cases, we found that 14 of 21 (67%) patients with double somatic tumors also had PIK3CA mutations, compared with 4 of 18 (22%) tumors from patients with Lynch syndrome, 2 of 10 (20%) tumors with MLH1 hypermethylation, and 12 of 78 (15%) tumors with microsatellite stability (P < .0001 for patients with double somatic tumors vs other subgroups). Mutations in PIK3CA were detected in all 13 patients with double somatic endometrial cancers (P = .04 compared with other subgroups). We did not detect BRAF mutations in patients with double somatic colorectal tumors or Lynch syndrome. We found highly similar results in a validation cohort from The Cancer Genome Atlas (113 patients with colorectal tumors, 178 endometrial tumors); 100% of double somatic cases had a somatic mutation in PIK3CA (P < .0001 compared with other subgroups). CONCLUSIONS: Most patients with colorectal or endometrial tumors with 2 or more somatic (but not germline) mutations in MMR proteins also have mutations in PIK3CA; mutations in PIK3CA are detected at substantially higher frequencies in these double somatic tumors than in other microsatellite-instability subgroups. PIK3CA mutation status might be used to identify a specific group of colorectal tumors, and to select treatment or determine prognosis.