Telomere-lengthening germline variants predispose to a syndromic papillary thyroid cancer subtype.

Papillary thyroid cancer (PTC) is the most common endocrine malignancy. 10% to 15% of individuals show familial clustering with three or more affected members, but the factors underlying this risk are unknown. In a group of recently studied individuals with POT1 pathogenic variants and ultra-long telomere length, PTC was the second most common solid tumor. We tested whether variants in POT1 and four other telomere-maintenance genes associated with familial cancer underlie PTC susceptibility. Among 470 individuals, we identified pathogenic or likely pathogenic variants in three genes encoding telomere-binding proteins: POT1, TINF2, and ACD. They were found in 4.5% and 1.5% of familial and unselected cases, respectively. Individuals harboring these variants had ultra-long telomere length, and 15 of 18 (83%) developed other cancers, of which melanoma, lymphoma, and sarcoma were most common. Among individuals with PTC and melanoma, 22% carried a deleterious germline variant, suggesting that a long telomere syndrome might be clinically recognizable. Successive generations had longer telomere length than their parents and, at times, developed more cancers at younger ages. Tumor sequencing identified a single oncogenic driver, BRAF p.Val600Glu, in 10 of 10 tumors studied, but no telomere-maintenance mechanism, including at the TERT promoter. These data identify a syndromic subset of PTCs with locus heterogeneity and telomere lengthening as a convergent mechanism. They suggest these germline variants lower the threshold to cancer by obviating the need for an acquired telomere-maintenance mechanism in addition to sustaining the longevity of oncogenic mutations.

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.

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.

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.

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.

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.

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.

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.

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.

NRAS isoforms differentially affect downstream pathways, cell growth, and cell transformation.

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.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease.

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<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis 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. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

The polymorphism rs944289 predisposes to papillary thyroid carcinoma through a large intergenic noncoding RNA gene of tumor suppressor type.

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.

Heritable polymorphism predisposes to high BAALC expression in acute myeloid leukemia.

Overexpression of the brain and acute leukemia, cytoplasmic (BAALC) gene is implicated in myeloid leukemogenesis and associated with poor outcome in both acute myeloid leukemia (AML) and acute lymphoblastic leukemia patients. Additionally, high BAALC expression occurs in glioblastoma, melanoma, and childhood gastrointestinal stroma tumors, suggesting an oncogenic role for BAALC. However, the mechanisms underlying the deregulated expression are unknown. We hypothesized that a common heritable genetic feature located in cis might account for overexpression of BAALC in an allele-specific manner. By sequencing the genomic region of BAALC we identified nine informative single nucleotide polymorphisms (SNPs) and tested them for a possible association with BAALC expression levels. We show that BAALC overexpression occurs in the presence of the T allele of SNP rs62527607[GT], which creates a binding site for the activating RUNX1 transcription factor in the BAALC promoter region. The mechanism is demonstrated experimentally in vitro using luciferase reporter assays and electrophoretic mobility shift assay (EMSA) analysis. The association of high BAALC expression with the T allele and its correlations with RUNX1 expresser status are shown in vivo in a test set (n = 253) and validation set (n = 105) of samples from cytogenetically normal AML patients from different populations. Thus, we identify a heritable genomic feature predisposing to overexpression of an oncogene, thereby possibly leading to enhanced AML leukemogenesis. Our findings further suggest that genomic variants might become useful tools in the practice of personalized medicine.

miR-3151 interplays with its host gene BAALC and independently affects outcome of patients with cytogenetically normal acute myeloid leukemia.

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.

The evaluation and analysis of irritable bowel syndrome-related short videos on social media (TikTok).

BACKGROUND: TikTok is one of the fastest growing social media platforms. Irritable bowel syndrome (IBS) has recently become a trending topic of interest among TikTok users. AIM: To better understand the quality and accuracy of information presented in the most popular IBS-relevant videos on TikTok. METHODS: We reviewed videos with the tag 'IBS'. We excluded those not relevant to IBS or lasting <10 s or >10 min. Baseline characteristics about the videos were collected. Two independent reviewers assessed each video using DISCERN and Patient Education Materials and Assessment Tool (PEMAT) tools, two validated instruments to assess the quality of patient education materials. RESULTS: Of 100 videos, 33% were uploaded by participants with a defined medical background. The median DISCERN score of videos uploaded by participants with a medical background was 2.43 (2.00-3.10); from participants with a non-medical background, it was 1.37 (1.23-1.70) (p < 0.01). The median PEMAT Understandability scores of videos uploaded by participants with or without a medical background were 92.86 (86.61-95.00) and 80.95 (75.76-89.58), respectively (p < 0.01). The median PEMAT Actionability scores of videos uploaded by participants with or without a medical background were 100.00 (66.67-100.00) and 0.00, respectively (0.00-45.83; p < 0.01). CONCLUSION: Videos posted by medical professionals are easier to understand and to act on, and are more reliable and unbiased, and more likely to recommend shared decision making about treatment.

Presumed Pathogenic Germ Line and Somatic Variants in African American Thyroid Cancer.

Background: African American (AA) thyroid cancer patients have worse prognoses than European Americans (EA), which has been attributed to both health care disparities and possible genetic differences. We investigated the impact of both germ line and somatic variants on clinical outcome in a cohort of AA nonmedullary thyroid cancer (NMTC) patients who had received therapeutic intervention from cancer centers. Methods: Whole-exome sequencing was performed on DNA from available blood/normal tissues (N = 37) and paired tumor samples (N = 32) collected from 37 and 29 AA NMTC patients, respectively. Variants with Combined Annotation Depletion Dependent (CADD) score of ≥20 and VarSome Clinical classification of likely pathogenic or pathogenic were classified as presumed pathogenic germ line or somatic variants (PPGVs/PPSVs). PPGVs/PPSVs in cancer-related genes and PPGVs in cardiovascular risk genes were further investigated, and PPGVs/PPSVs associated with African (AFR) ancestry were identified. Results: Among 17 PPGVs identified in 16 cancer predisposition or known cancer-related genes, only WRN was previously known to associate with NMTC predisposition. Among PPSVs, BRAFV600E was most the prevalent and detected in 12 of the 29 (41%) tumors. Examining PPGVs/PPSVs among three patients who died from NMTC, one patient who died from papillary thyroid carcinoma/anaplastic thyroid carcinoma (PTC/ATC) led us to speculate that the PPGV ERCC4R799W may have increased the risk of PPSV TP53R273H acquisition. Among PPGVs identified in 18 cardiovascular risk genes, PPGVs in SC5NA, GYG1, CBS, CFTR, and SI are known to have causal and pathogenic implications in cardiovascular disease. Conclusion: In this cohort, most AA-NMTC patients exhibit favorable outcomes after therapeutic intervention given at cancer centers, suggesting that health care disparity is the major contributor for worse prognoses among AA-NMTC patients. Nevertheless, the clinical impact of PPGVs that might facilitate the acquisition of TP53 tumor mutations, and/or PPGVs that predispose individuals to adverse cardiovascular events, which could be exacerbated by therapy-induced cardiotoxicity, needs to be further explored. Integrated analysis of PPGV/PPSV profiles among NMTC patients with different stages of disease may help to identify NMTC patients who require close monitoring or proactive intervention.

PDPR Gene Variants Predisposing to Papillary Thyroid Cancer.

Background: Papillary thyroid cancer (PTC) is the predominant subtype of thyroid cancer (THCA), and it can cluster in families with an autosomal dominant (AD) inheritance pattern. The aim of this study was to identify novel genes and mechanisms underlying PTC susceptibility. Methods: Our previous investigation of 17 AD PTC families led us to conduct a deeper analysis on one family (Family Q) with whole-genome sequencing data from 3 PTC-affected individuals. In addition, 323 sporadic THCA cases from Avatar data and 12 familial adenomatous polyposis (FAP) individuals with secondary THCA were screened for pyruvate dehydrogenase phosphatase regulatory (PDPR) variants. CRISPR-Cas9 was used to create PDPR-deficient THCA (TPC1) and transformed normal thyroid cell lines (N-Thyori3-1) to study the metabolic consequences of PDPR loss. Results: We found truncating PDPR splice donor variants (NM_017990.4:c.361 + 1G>C) in all affected PTC Family Q members, and another PDPR splice donor variant (NM_017990.4:c.443 + 1G>C) in a sporadic PTC case. In addition, an ultra-rare missense variant was found in an FAP-PTC patient. The PDPR-deficient cells presented with elevated phosphorylation of pyruvate dehydrogenase and altered glucose metabolism, implying that PDPR plays an essential part in regulating glucose metabolism in thyroid cells. Conclusions: Our finding of novel truncating germline variants in PDPR in Family Q and additional cohorts suggests a role for PDPR loss in PTC predisposition. Also, somatic and RNA sequencing from the thyroid carcinoma (Firehouse Legacy) data showed that PDPR gene expression is much lower in THCA tumor tissue compared with matching normal tissue. Thus, PDPR appears to have a loss of function effect on THCA tumorigenesis.

Pheochromocytoma in rats with multiple endocrine neoplasia (MENX) shares gene expression patterns with human pheochromocytoma.

Pheochromocytomas are rare neoplasias of neural crest origin arising from chromaffin cells of the adrenal medulla and sympathetic ganglia (extra-adrenal pheochromocytoma). Pheochromocytoma that develop in rats homozygous for a loss-of-function mutation in p27Kip1 (MENX syndrome) show a clear progression from hyperplasia to tumor, offering the possibility to gain insight into tumor pathobiology. We compared the gene-expression signatures of both adrenomedullary hyperplasia and pheochromocytoma with normal rat adrenal medulla. Hyperplasia and tumor show very similar transcriptome profiles, indicating early determination of the tumorigenic signature. Overrepresentation of developmentally regulated neural genes was a feature of the rat lesions. Quantitative RT-PCR validated the up-regulation of 11 genes, including some involved in neural development: Cdkn2a, Cdkn2c, Neurod1, Gal, Bmp7, and Phox2a. Overexpression of these genes precedes histological changes in affected adrenal glands. Their presence at early stages of tumorigenesis indicates they are not acquired during progression and may be a result of the lack of functional p27Kip1. Adrenal and extra-adrenal pheochromocytoma development clearly follows diverged molecular pathways in MENX rats. To correlate these findings to human pheochromocytoma, we studied nine genes overexpressed in the rat lesions in 46 sporadic and familial human pheochromocytomas. The expression of GAL, DGKH, BMP7, PHOX2A, L1CAM, TCTE1, EBF3, SOX4, and HASH1 was up-regulated, although with different frequencies. Immunohistochemical staining detected high L1CAM expression selectively in 27 human pheochromocytomas but not in 140 nonchromaffin neuroendocrine tumors. These studies reveal clues to the molecular pathways involved in rat and human pheochromocytoma and identify previously unexplored biomarkers for clinical use.

Amnionless, essential for mouse gastrulation, is mutated in recessive hereditary megaloblastic anemia.

The amnionless gene, Amn, on mouse chromosome 12 encodes a type I transmembrane protein that is expressed in the extraembryonic visceral layer during gastrulation. Mice homozygous with respect to the amn mutation generated by a transgene insertion have no amnion. The embryos are severely compromised, surviving to the tenth day of gestation but seem to lack the mesodermal layers that normally produce the trunk. The Amn protein has one transmembrane domain separating a larger, N-terminal extracellular region and a smaller, C-terminal cytoplasmic region. The extracellular region harbors a cysteine-rich domain resembling those occurring in Chordin, found in Xenopus laevis embryos, and Sog, found in Drosophila melanogaster. As these cysteine-rich domains bind bone morphogenetic proteins (Bmps), it has been speculated that the cysteine-rich domain in Amn also binds Bmps. We show that homozygous mutations affecting exons 1-4 of human AMN lead to selective malabsorption of vitamin B12 (a phenotype associated with megaloblastic anemia 1, MGA1; OMIM 261100; refs. 5,6) in otherwise normal individuals, suggesting that the 5' end of AMN is dispensable for embryonic development but necessary for absorption of vitamin B12. When the 5' end of AMN is truncated by mutations, translation is initiated from alternative downstream start codons.

The RCAN1.4 Metastasis Suppressor Is Hypermethylated at Intron 1 in Thyroid Cancer.

Background: Regulator of calcineurin 1.4 (RCAN1.4) is a functionally downregulated metastasis progression suppressor (MPS) in thyroid cancer; however, the mechanisms for RCAN1.4 loss in thyroid cancer have not yet been reported. The RCAN1.4 promoter and gene contain several cytosine-guanine (CG)-rich regions, some of which are reported to be hypermethylated in nonthyroid tissues. We, therefore, hypothesized that RCAN1.4 downregulation in thyroid cancer was in part due to hypermethylation. Methods: Studies were performed in 5 thyroid cancer cell lines (TPC1, FTC133, BCPAP, C643, and 8505C) with different genetic drivers, and in 18 paired normal and thyroid cancer human thyroid cancer tissues. Basal RCAN1.4 messenger RNA (mRNA) and protein levels were assessed in all of the cell lines. Cell lines with lowest RCAN1.4 expression levels were treated with the DNA methyl transferase inhibitor, decitabine. Normal/tumor tissue pairs were analyzed for methylation of three CG-rich regions both by capture of methylated DNA by MBD2 protein and by methylation-specific polymerase chain reaction (MSPCR). Results: In all assessed cell lines, RCAN1.4 mRNA and protein levels increased after decitabine treatment. In silico analysis of the RCAN1.4 gene identified 3 CG-rich regions as possible methylation targets: 1 in the proximal promoter and 2 in intron 1. Hypermethylation of the intron 1 CG-rich regions was identified by both the capture method and MSPCR. In contrast, hypermethylation of the CG-rich region of the proximal promoter was not identified. Gene expression confirmed that hypermethylation in thyroid cancer samples in intron 1 of RCAN1.4 was associated with lower levels of RCAN1.4 mRNA. Finally, the cancer samples demonstrated increased NFE2L3 expression, a downstream marker of functional RCAN1.4 loss. Conclusions: The MPS gene, RCAN1.4, is downregulated in thyroid cancer cells and human thyroid cancer in part by hypermethylation of CG-rich regions in intron 1.