Sandya Liyanarachchi

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.

Background: Germline pathogenic variants in CHEK2 are associated with a moderate increase in the lifetime risk for breast cancer. Increased risk for other cancers, including non-medullary thyroid cancer (NMTC), has also been suggested. To date, data implicating CHEK2 variants in NMTC predisposition primarily derive from studies within Poland, driven by a splice site variant (c.444+1G>A) that is uncommon in other populations. In contrast, the predominant CHEK2 variants in non-Polish populations are c.1100del and c.470T>C/p.I157T, representing 61.1% and 63.8%, respectively, of all CHEK2 pathogenic variants in two large US-based commercial laboratory datasets. To further delineate the impact of common CHEK2 variants on thyroid cancer, we aimed to investigate the association of three CHEK2 founder variants (c.444+1G>A, c.1100del, and c.470T>C/p.Ile157Thr) on NMTC susceptibility in three groups of unselected NMTC patients. Methods: The presence of three CHEK2 founder variants were assessed within three groups: (1) 1544 NMTC patients (and 1593 controls) from previously published GWAS analyses, (2) 789 NMTC patients with germline exome sequencing (ORIEN Avatar), and (3) 499 NMTC patients with germline sequence data available in The Cancer Genome Atlas (TCGA). A case-control study design was utilized with odds ratios calculated by comparison of all three groups with the OSU GWAS control group. Results: The predominant Polish variant (c.444+1G>A) was present in only one case. The proportion of patients with c.1100del was 0.92% in the GWAS group, 1.65% in the ORIEN Avatar group, and 0.80% in the TCGA group. The odds ratios (with 95% confidence intervals) for NMTC associated with c.1100del were 1.71 (0.73-4.29), 2.64 (0.95-7.63), and 2.5 (0.63-8.46), respectively. The proportion of patients with c.470T>C/p.I157T was 0.91% in the GWAS group, 0.76% in the ORIEN Avatar group, and 0.80% in the TCGA group, respectively. The odds ratios (with 95% confidence intervals) for NMTC associated with c.470T>C/p.I157T were 1.75 (0.74-4.39), 1.52 (0.42-4.96), and 2.31 (0.58-7.90), respectively. Conclusions: Our analyses of unselected patients with NMTC suggest that CHEK2 variants c.1100del and c.470T>C/p.I157T have only a modest impact on thyroid cancer risk. These results provide important information for providers regarding the relatively low magnitude of thyroid cancer risk associated with these CHEK2 variants.

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, BRAF ^V600E 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 ERCC4 ^R799W may have increased the risk of PPSV TP53 ^R273H 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.

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.

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

Familial non-medullary thyroid cancer (FNMTC) is clinically defined as two or more first-degree relatives with NMTC and appears to follow an autosomal dominant inheritance pattern. Approximately 5–7% of NMTC is hereditary and affects multiple generations with a young age of onset. The primary aim of this study was to determine the age-specific penetrance of NMTC in individuals from a large family with FNMTC with a previously identified private mutation at 4q32, with a secondary aim to determine the penetrance for benign thyroid disease in this family. We present a large family with NMTC in which we had previously described a culpable mutation. Participants provided their personal medical history and family history. The germline 4q32 A > C mutation was detected in 34 of 68 tested individuals. Age-specific penetrance of thyroid cancer and benign thyroid disease was determined using the inverted Kaplan–Meier method of segregation analysis. Individuals who tested positive for the 4q32 mutation have a 68.9% (95% CI 46.5–88.7) risk of developing thyroid cancer by age 70 and a 65.3% (95% CI 46.0–83.8) risk of developing benign thyroid disease by age 70. The 4q32 A > C mutation significantly increases the risk to develop thyroid cancer but not benign thyroid disease in members of this family. The female:male sex ratio of 1.33 that we observed in affected mutation carriers differs greatly from the ratio of approximately 3:1 observed in PTC, supporting a central role of the mutation. Early thyroid surveillance with annual ultrasound is recommended to individuals testing positive for this private familial mutation.

Background: Papillary thyroid cancer (PTC) is reported to be highly heritable in epidemiological studies. Genome-wide association studies (GWAS) have uncovered several variants associated with PTC predisposition. It remains unknown whether these variants might contribute to better clinical stratification of PTC patients. Methods: In order to assess the usefulness of germline genetic analyses in the management of PTC patients, the genotypes of five variants (rs965513, rs944289, rs116909374, rs2439302, and rs966423) were determined in 1216 PTC patients and 1416 controls. Additionally, the expression of seven genes located close to GWAS variants ( PTCSC3 , MBIP , NKX2-1 , FOXE1 , DIRC3 , PTCSC2 , and NRG1 ) were measured in 73 PTC paired tumor/normal tissues, respectively. Next, the association was analyzed between the genotypes of the germline variants and the levels of gene expression with clinical/pathological features such as age, sex, TNM staging, multifocality status, extrathyroidal expansion, and MACIS score. Results: The risk allele of rs965513 was associated with larger tumor size ( p = 0.025) and extrathyroidal expansion (odd ratio [OR] = 1.29, p = 0.045). The variant rs2439302 showed association with lymph node metastasis (OR = 1.24, p = 0.016), and multifocality status of the tumor (OR = 1.24, p = 0.012). The expression of MBIP was associated with T stage ( p = 0.010). MBIP and PTCSC3 displayed lower expression in PTC tissue in males than in females ( p = 0.025 and p = 0.036, respectively). NKX2-1 displayed lower expression in patients with N1 stage ( p = 0.040). Conclusions: The studied germline risk alleles predisposing to PTC were associated with a more aggressive course of the disease reflected by larger tumor diameter, higher multifocality rate, and more advanced N stage at the time of diagnosis. These results show that germline variants not only predispose to PTC but also might impact its clinical course. However, these associations were only moderate, and further large multi-ethnic studies are required to evaluate the usefulness of these germline variants in the clinical stratification of PTC patients.

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

Inherited malabsorption of cobalamin (Cbl) causes hematological and neurological abnormalities that can be fatal. Three genes have been implicated in Cbl malabsorption; yet, only about 10% of ~400-500 reported cases have been molecularly studied to date. Recessive mutations in CUBN or AMN cause Imerslund-Gräsbeck Syndrome (IGS), while recessive mutations in GIF cause Intrinsic Factor Deficiency (IFD). IGS and IFD differ in that IGS usually presents with proteinuria, which is not observed in IFD. The genetic heterogeneity and numerous differential diagnoses make clinical assessment difficult.

Background: Imerslund-Gräsbeck syndrome (IGS) was described just over 50 years ago by Olga Imerslund and Ralph Gräsbeck and colleagues. IGS is caused by specific malabsorption of cobalamin (Cbl) due to bi-allelic mutations in either the cubilin gene (CUBN) or the human amnionless homolog (AMN). Mutations in the two genes are commonly seen in founder populations or in societies with a high degree of consanguineous marriages. One particular mutation in AMN, c.208-2A>G, causing an out-of-frame loss of exon 4 in the mRNA, is responsible for some 15% of IGS cases globally. We present evidence that this founder mutation causes a substantial percentage of cases among diverse ethnicities and that the mutation is as old as human civilization. Methods: Partial genotyping indicated a founder event but its presence in diverse peoples of Arabic, Turkish, Jewish, and Hispanic ancestry suggested that the mutation might be recurrent. We therefore studied the flanking sequence spanning 3.5 Mb to elucidate the origin of the haplotype and estimate the age of the mutation using a Bayesian inference method based on observed linkage disequilibrium. Results: The mutation's distribution, the size of the shared haplotype, and estimates of growth rate and carrier frequency indicated that the mutation was a single prehistoric event. Dating back to the ancient Middle East around 11,600 BC, the mutation predates the advent of writing, farming, and the monotheistic religions of the region. Conclusions: This mutation causes over 50% of the IGS cases among Arabic, Turkish, and Sephardic Jewish families, making it a primary target for genetic screening among diverse IGS cases originating from the Middle East. Thus, rare founder mutations may cause a substantial number of cases, even among diverse ethnicities not usually thought to be related.