Genome-wide association study identifies a novel susceptibility gene for serum TSH levels in Chinese populations

Novel insights of disulfidptosis-mediated immune microenvironment regulation in atherosclerosis based on bioinformatics analyses

Atherosclerosis (AS) is the leading cause of coronary heart disease, which is the primary cause of death worldwide. Recent studies have identified disulfidptosis as a new type of cell death that may be involved in onset and development of many diseases. However, the role of disulfidptosis in AS is not clear. In this study, bioinformatics analysis and experiments in vivo and in vitro were performed to evaluate the potential relationship between disulfidptosis and AS. AS-related sequencing data were obtained from the Gene Expression Omnibus (GEO). Bioinformatics techniques were used to evaluate differentially expressed genes (DEGs) associated with disulfidptosis-related AS. Hub genes were screened using least absolute shrinkage and selection operator (LASSO) and random forests (RF) methods. In addition, we established a foam cell model in vitro and an AS mouse model in vivo to verify the expressions of hub genes. In addition, we constructed a diagnostic nomogram with hub genes to predict progression of AS. Finally, the consensus clustering method was used to establish two different subtypes, and associations between subtypes and immunity were explored. As the results, 9 disulfidptosis-related AS DEGs were identified from GSE28829 and GSE43292 datasets. Evaluation of DEGs using LASSO and RF methods resulted in identification of 4 hub genes (CAPZB, DSTN, MYL6, PDLIM1), which were analyzed for diagnostic value using ROC curve analysis and verified in vitro and in vivo. Furthermore, a nomogram including hub genes was established that accurately predicted the occurrence of AS. The consensus clustering algorithm was used to separate patients with early atherosclerotic plaques and patients with advanced atherosclerotic plaques into two disulfidptosis subtypes. Cluster B displayed higher levels of infiltrating immune cells, which indicated that patients in cluster B may have a positive immune response for progression of AS. In summary, disulfidptosis-related genes including CAPZB, DSTN, MYL6, and PDLIM1 may be diagnostic markers and therapeutic targets for AS. In addition, these genes are closely related to immune cells, which may inform immunotherapy for AS.

Genome-wide association studies of thyroid-related hormones, dysfunction, and autoimmunity among 85,421 Chinese pregnancies

Maintaining normal thyroid function is crucial in pregnancy, yet thyroid dysfunction and the presence of thyroid peroxidase antibodies (TPOAb) affect 0.5% to 18% of pregnant women. Here, we conducted a genome-wide association study (GWAS) of eight thyroid traits, including two thyroid-related hormones, four thyroid dysfunctions, and two thyroid autoimmunity measurements among 85,421 Chinese pregnant women to investigate the genetic basis of thyroid function during pregnancy. Our study identified 176 genetic loci, including 125 previously unknown genome-wide associations. Joint epidemiological and Mendelian randomization analyses revealed significant associations between the gestational thyroid phenotypes and gestational complications, birth outcomes, and later-age health outcomes. Specifically, genetically elevated thyroid-stimulating hormone (TSH) levels during pregnancy correlated with lower glycemic levels, reduced blood pressure, and longer gestational duration. Additionally, TPOAb and thyroid functions during pregnancy share genetic correlations with later-age thyroid and cardiac disorders. These findings provide insights into the genetic determinants of thyroid traits during pregnancy, which may lead to new therapeutics, early pre-diagnosis and preventive strategies starting from early adulthood.

Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle

Locating the genetic variation of important livestock and poultry economic traits is essential for genetic improvement in breeding programs. Identifying the candidate genes for the productive ability of Huaxi cattle was one crucial element for practical breeding. Based on the genotype and phenotype data of 1,478 individuals and the RNA-seq data of 120 individuals contained in 1,478 individuals, we implemented genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and Fisher's combined test (FCT) to identify the candidate genes for the carcass trait, the weight of longissimus dorsi muscle (LDM). The results indicated that GWAS, TWAS, and FCT identified seven candidate genes for LDM altogether: PENK was located by GWAS and FCT, PPAT was located by TWAS and FCT, and XKR4, MTMR3, FGFRL1, DHRS4, and LAP3 were only located by one of the methods. After functional analysis of these candidate genes and referring to the reported studies, we found that they were mainly functional in the progress of the development of the body and the growth of muscle cells. Combining advanced breeding techniques such as gene editing with our study will significantly accelerate the genetic improvement for the future breeding of Huaxi cattle.

Suggestive evidence of the genetic association of TMOD1 and PTCSC2 polymorphisms with thyroid carcinoma in the Chinese Han population

BACKGROUND

The purpose of this study was to survey the associations of six single nucleotide polymorphisms (SNPs) in the TMOD1 and PTCSC2 genes with thyroid carcinoma (TC).

METHOD

Peripheral blood samples were obtained from 510 patients with TC and 509 normal controls. Six SNPs were genotyped by the Agena MassARRAY platform. Logistic regression was used to evaluate the association between SNPs and TC susceptibility by calculating odds ratios (ORs) and 95% confidence intervals (CIs). SNP-SNP interactions were analyzed by multifactor dimensionality reduction (MDR).

RESULTS

Our study showed that rs925489 (OR = 1.45, p = 0.011) and rs965513 (OR = 1.40, p = 0.021) were significantly associated with an increased risk of TC. Rs10982622 decreased TC risk (OR = 0.74, p = 0.025). Further stratification analysis showed that rs10982622 reduced the susceptibility to TC in patients aged ≤ 45 years (OR = 0.69, p = 0.019) and in females (OR = 0.61, p = 0.014). Rs925489 increased TC risk in people aged > 45 years (OR = 1.54, p = 0.044) and in males (OR = 2.34, p = 0.003). In addition, rs965513 was related to an increased risk of TC in males (OR = 2.14, p = 0.007). Additionally, haplotypes in the block (rs925489|rs965513) significantly increased TC risk (p < 0.05). The best predictive model for TC was the combination of rs1052270, rs10982622, rs1475545, rs16924016, and rs925489.

CONCLUSION

TMOD1 and PTCSC2 polymorphisms were separately correlated with a remarkable decrease and increase in TC risk based on the analysis.

Genetic Susceptibility to Joint Occurrence of Polycystic Ovary Syndrome and Hashimoto's Thyroiditis: How Far Is Our Understanding?

Polycystic ovary syndrome (PCOS) and Hashimoto’s thyroiditis (HT) are endocrine disorders that commonly occur among young women. A higher prevalence of HT in women with PCOS, relative to healthy individuals, is observed consistently. Combined occurrence of both diseases is associated with a higher risk of severe metabolic and reproductive complications. Genetic factors strongly impact the pathogenesis of both PCOS and HT and several susceptibility loci associated with a higher risk of both disorders have been identified. Furthermore, some candidate gene polymorphisms are thought to be functionally relevant; however, few genetic variants are proposed to be causally associated with the incidence of both disorders together.

Genetic Determinants of Lung Cancer Prognosis in Never Smokers: A Pooled Analysis in the International Lung Cancer Consortium

BACKGROUND

Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO).

METHODS

Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression-based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887).

RESULTS

One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41-2.07), P = 3.60 × 10-8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10-7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P < 10-6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38-0.66), P = 6.92 × 10-7] and Japanese populations [HR = 0.79 (0.67-0.94), P = 0.007].

CONCLUSIONS

Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants.

IMPACT

We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.

Automated feature extraction from population wearable device data identified novel loci associated with sleep and circadian rhythms

Wearable devices have been increasingly used in research to provide continuous physical activity monitoring, but how to effectively extract features remains challenging for researchers. To analyze the generated actigraphy data in large-scale population studies, we developed computationally efficient methods to derive sleep and activity features through a Hidden Markov Model-based sleep/wake identification algorithm, and circadian rhythm features through a Penalized Multi-band Learning approach adapted from machine learning. Unsupervised feature extraction is useful when labeled data are unavailable, especially in large-scale population studies. We applied these two methods to the UK Biobank wearable device data and used the derived sleep and circadian features as phenotypes in genome-wide association studies. We identified 53 genetic loci with p<5×10-8 including genes known to be associated with sleep disorders and circadian rhythms as well as novel loci associated with Body Mass Index, mental diseases and neurological disorders, which suggest shared genetic factors of sleep and circadian rhythms with physical and mental health. Further cross-tissue enrichment analysis highlights the important role of the central nervous system and the shared genetic architecture with metabolism-related traits and the metabolic system. Our study demonstrates the effectiveness of our unsupervised methods for wearable device data when additional training data cannot be easily acquired, and our study further expands the application of wearable devices in population studies and genetic studies to provide novel biological insights.

The Genetic Basis of Thyroid Function: Novel Findings and New Approaches

CONTEXT

Genetic factors are major determinants of thyroid function. Over the last two decades, multiple genetic variants have been associated with variations in normal range thyroid function tests. Most recently, a large-scale genome-wide association study (GWAS) doubled the number of known variants associated with normal range thyrotropin (TSH) and free thyroxine (FT4) levels.

EVIDENCE ACQUISITION

This review summarizes the results of genetic association studies on normal range thyroid function and explores how these genetic variants can be used in future studies to improve our understanding of thyroid hormone regulation and disease.

EVIDENCE SYNTHESIS

Serum TSH and FT4 levels are determined by multiple genetic variants on virtually all levels of the hypothalamus-pituitary-thyroid (HPT) axis. Functional follow-up studies on top of GWAS hits has the potential to discover new key players in thyroid hormone regulation, as exemplified by the identification of the thyroid hormone transporter SLC17A4 and the metabolizing enzyme AADAT. Translational studies may use these genetic variants to investigate causal associations between thyroid function and various outcomes in Mendelian Randomization (MR) studies, to identify individuals with an increased risk of thyroid dysfunction, and to predict the individual HPT axis setpoint.

CONCLUSIONS

Recent genetic studies have greatly improved our understanding of the genetic basis of thyroid function, and have revealed novel pathways involved in its regulation. In addition, these findings have paved the way for various lines of research that can improve our understanding of thyroid hormone regulation and thyroid diseases, as well as the potential use of these markers in future clinical practice.

Variants in oxidative stress-related genes affect the chemosensitivity through Nrf2-mediated signaling pathway in biliary tract cancer

BACKGROUND

Oxidative stress and their effectors play critical roles in carcinogenesis and chemoresistance. However, the role of oxidative stress-related genes variants in biliary tract cancer (BTC) chemoresistance remains unknown. In this work, we aim to investigate oxidative stress-dependent molecular mechanisms underlying chemoresistance, and find potential biomarkers to predict chemotherapy response for BTC.

METHODS

Sixty-six SNPs in 21 oxidative stress-related genes were genotyped and analyzed in 367 BTC patients. Immunoblot, immunohistochemical, immunofluorescent, quantitative PCR, chromatin immunoprecipitation analysis and study of animal xenograft models were performed to discover oxidative stress-related susceptibility genes underlying chemoresistance mechanism of BTC.

FINDINGS

We found that 3 functional polymorphisms (CAT_rs769217, GPX4_rs4807542, and GSR_rs3779647), which were shown to affect their respective gene expression levels, modified the effect of chemotherapy on overall survival (OS). We then demonstrated that knockdown of GPX4, CAT, or GSR induced chemoresistance through elevation of ROS level and activation of Nrf2-ABCG2 pathway in BTC cell lines. Moreover, the association between Nrf2 expression and BTC prognosis is only found in patients who received chemotherapy. Knockdown of Nrf2 enhanced chemosensitivity or even eliminated postoperative recurrence in BTC xenograft mouse models. Importantly, upon chemotherapy treatment patients harboring high oxidative stress-related score received higher survival benefit from adjuvant chemotherapy compared with patients with low oxidative stress-related score.

INTERPRETATION

The result of our study suggests, for the first time, that the oxidative stress-related score calculated by combining variations in CAT, GPX4, and GSR or Nrf2 expression could be used for predicting the chemosensitivity of BTC patients. FUND: This work was supported by the National Science Foundation of China, Foundation of Shanghai Shen Kang Hospital Development Center, and Shanghai Outstanding Academic Leaders Plan.

Mutation screening of the TSHR gene in 220 Chinese patients with congenital hypothyroidism

BACKGROUND

Defects in the human thyroid stimulating hormone receptor (TSHR) gene are reported to be one of the causes of congenital hypothyroidism (CH). We aimed to identify mutations in Chinese patients with CH and analyze the relationships between TSHR phenotypes and clinical phenotypes.

METHODS

220 patients with primary CH were screened for TSHR mutations by performing next-generation sequencing. All the exons and exon-intron boundaries of TSHR were analyzed. The function of 8 mutants in TSHR were further investigated in vitro.

RESULTS

Among 220 patients with CH, 15 distinct TSHR mutations were identified in 13 patients (5.91%, 13/220, including our previous reported 110 patients, carried with 10 mutations in 8 patients). We found five distinct mutations in the additional cohort of 110 CH patients and identified 7 mutations (including a novel mutation, p.S567R) were loss-of-function mutations.

CONCLUSION

Our study indicated that the prevalence of TSHR mutations was 5.91% among studied Chinese patients with CH. One novel TSHR variant was found and four genetic alterations revealed important role of the Ile216, Ala275, Asn372, Ser567 residues in signaling.

Predicting adult height from DNA variants in a European-Asian admixed population

Accurate genomic profiling for adult height is of high practical relevance in forensics genetics. Adult height is a classical reference trait in the field of human complex trait genetics characterized by highly polygenic nature and relatively high heritability. A meta-analysis of genome-wide association studies by the Genetic Investigation of Anthropocentric Traits (GIANT) consortium has identified 697 DNA variants associated with adult height in Europeans; however, whether these variants will still be informative in non-Europeans is still in question. The present study investigated the predictive power of these 697 height-associated SNPs in 687 Uyghurs of European-Asian admixed origin. Among all GIANT SNPs, 11% showed nominally significant association (6.78 × 10^-4 < p < 0.05) with adult height in the Uyghur population and among the significant SNPs 77% of allele effects were in the same direction as those in Europeans reported in the GIANT study. Fitting linear and logistic models using a polygenic score consisting of all GIANT SNPs resulted in an 80-20 cross-validated mean R² of 10.08% (95% CI 3.16-18.40%) for quantitative height prediction and a mean AUC value of 0.65 (95% CI 0.57-0.72%) for qualitative "above average" prediction. Fine-tuning the SNP set using their association p values considerably improved the prediction results (number of SNPs = 62, R² = 15.59%, 95% CI 6.80-25.71%; AUC = 0.70, 95% CI 62-0.77) in the Uyghurs. Overall, our findings demonstrate substantial differences between the European and Asian populations in the genetics of adult height, emphasizing the importance of population heterogeneity underlying the genetic architecture of adult height.

Functional Urate-Associated Genetic Variants Influence Expression of lincRNAs LINC01229 and MAFTRR

Genetic variation in the genomic regulatory landscape likely plays a crucial role in the pathology of disease. Non-coding variants associated with disease can influence the expression of long intergenic non-coding RNAs (lincRNAs), which in turn function in the control of protein-coding gene expression. Here, we investigate the function of two independent serum urate-associated signals (SUA1 and SUA2) in close proximity to lincRNAs and an enhancer that reside ∼60 kb and ∼300 kb upstream of MAF, respectively. Variants within SUA1 are expression quantitative trait loci (eQTL) for LINC01229 and MAFTRR, both co-expressed with MAF. We have also identified that variants within SUA1 are trans-eQTL for genes that are active in kidney- and serum urate-relevant pathways. Serum urate-associated variants rs4077450 and rs4077451 within SUA2 lie within an enhancer that recruits the transcription factor HNF4α and forms long range interactions with LINC01229 and MAFTRR. The urate-raising alleles of rs4077450 and rs4077451 increase enhancer activity and associate with increased expression of LINC01229. We show that the SUA2 enhancer region drives expression in the zebrafish pronephros, recapitulating endogenous MAF expression. Depletion of MAFTRR and LINC01229 in HEK293 cells in turn lead to increased MAF expression. Collectively, our results are consistent with serum urate variants mediating long-range transcriptional regulation of the lincRNAs LINC01229 and MAFTRR and urate relevant genes (e.g., SLC5A8 and EHHADH) in trans.

Correlation of TSHR and CTLA-4 Single Nucleotide Polymorphisms with Graves Disease

This study was designed to explore the association between Graves disease (GD) and thyroid-stimulating hormone receptor (TSHR) and cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) single nucleotide polymorphisms (SNPs). We studied a total of 1217 subjects from a Han population in northern Anhui province in China. Six SNPs within TSHR (rs179247, rs12101261, rs2284722, rs4903964, rs2300525, and rs17111394) and four SNPs within CTLA-4 (rs10197319, rs231726, rs231804, and rs1024161) were genotyped via a Taqman probe technique using a Fluidigm EP1 platform. The TSHR alleles rs179247-G, rs12101261-C, and rs4903964-G were negatively correlated with GD, whereas the rs2284722-A and rs17111394-C alleles were positively correlated with GD. Analyzing TSHR SNPs at rs179247, rs2284722, rs12101261, and rs4903964 yielded 8 different haplotypes. There were positive correlations between GD risk and the haplotypes AGTA and AATA (OR = 1.27, 95%CI = 1.07-1.50, P = 0.005; OR = 1.45, 95%CI = 1.21-1.75, P < 0.001, respectively). There were negative correlations between GD risk and the haplotype GGCG (OR = 0.56, 95%CI = 0.46-0.67, P < 0.001). With respect to haplotypes based on SNPs at the TSHR rs2300525 and rs17111394 loci, the CC haplotype was positively correlated with GD risk (OR = 1.32, 95%CI = 1.08-1.60, P = 0.006). Analyzing CTLA-4 SNPs at rs231804, rs1024161, and rs231726 yielded four haplotypes, of which AAA was positively correlated with GD risk (OR = 1.21, 95%CI = 1.02-1.43, P = 0.029). Polymorphisms at rs179247, rs12101261, rs2284722, rs4903964, and rs17111394 were associated with GD susceptibility. Haplotypes of both TSHR and CTLA-4 were additionally related to GD risk.

Genome-Wide Association Study Reveals Distinct Genetic Susceptibility of Thyroid Nodules From Thyroid Cancer

CONTEXT

Thyroid nodules are very common, and 7% to 15% of them are diagnosed as thyroid cancer. However, the inherited genetic risk factors for thyroid nodules and their associations with thyroid cancer remain unknown.

OBJECTIVE

To identify the genetic variants associated with susceptibility to thyroid nodules in comparison with thyroid cancer.

DESIGN AND SETTING

We performed a three-stage genome-wide association study for thyroid nodules. The discovery stage involved a genome-wide scan of 811 subjects with thyroid nodules and 691 subjects with a normal thyroid from a population-based cohort. Replication studies were conducted in an additional 1981 cases and 3100 controls from the participants of a health checkup. We also performed expression quantitative trait loci analysis of public data.

RESULTS

The most robust association was observed in TRPM3 (rs4745021) in the joint analysis (OR, 1.26; P = 6.12 × 10-8) and meta-analysis (OR, 1.28; P = 2.11 × 10-8). Signals at MBIP/NKX2-1 were replicated but did not reach genome-wide significance in the joint analysis (rs2415317, P = 4.62 × 10-5; rs944289, P = 8.68 × 10-5). The expression quantitative trait loci analysis showed that TRPM3 expression was associated with the rs4745021 genotype in thyroid tissues.

CONCLUSIONS

To the best of our knowledge, we have performed the first genome-wide association study of thyroid nodules and identified a susceptibility locus associated with thyroid nodules, suggesting that thyroid nodules have a genetic predisposition distinct from that of thyroid cancer.

The genetic characteristics of congenital hypothyroidism in China by comprehensive screening of 21 candidate genes

OBJECTIVE

Congenital hypothyroidism (CH), the most common neonatal metabolic disorder, is characterized by impaired neurodevelopment. Although several candidate genes have been associated with CH, comprehensive screening of causative genes has been limited.

DESIGN AND METHODS

One hundred ten patients with primary CH were recruited in this study. All exons and exon-intron boundaries of 21 candidate genes for CH were analyzed by next-generation sequencing. And the inheritance pattern of causative genes was analyzed by the study of family pedigrees.

RESULTS

Our results showed that 57 patients (51.82%) carried biallelic mutations (containing compound heterozygous mutations and homozygous mutations) in six genes (DUOX2, DUOXA2, DUOXA1, TG, TPO and TSHR) involved in thyroid hormone synthesis. Autosomal recessive inheritance of CH caused by mutations in DUOX2, DUOXA2, TG and TPO was confirmed by analysis of 22 family pedigrees. Notably, eight mutations in four genes (FOXE1, NKX2-1, PAX8 and HHEX) that lead to thyroid dysgenesis were identified in eight probands. These mutations were heterozygous in all cases and hypothyroidism was not observed in parents of these probands.

CONCLUSIONS

Most cases of congenital hypothyroidism in China were caused by thyroid dyshormonogenesis rather than thyroid dysgenesis. This study identified previously reported causative genes for 57/110 Chinese patients and revealed DUOX2 was the most frequently mutated gene in these patients. Our study expanded the mutation spectrum of CH in Chinese patients, which was significantly different from Western countries.

Genome-Wide Association Studies of Autoimmune Thyroid Diseases, Thyroid Function, and Thyroid Cancer

Thyroid diseases, including autoimmune thyroid diseases and thyroid cancer, are known to have high heritability. Family and twin studies have indicated that genetics plays a major role in the development of thyroid diseases. Thyroid function, represented by thyroid stimulating hormone (TSH) and free thyroxine (T4), is also known to be partly genetically determined. Before the era of genome-wide association studies (GWAS), the ability to identify genes responsible for susceptibility to thyroid disease was limited. Over the past decade, GWAS have been used to identify genes involved in many complex diseases, including various phenotypes of the thyroid gland. In GWAS of autoimmune thyroid diseases, many susceptibility loci associated with autoimmunity (human leukocyte antigen [HLA], protein tyrosine phosphatase, non-receptor type 22 [PTPN22], cytotoxic T-lymphocyte associated protein 4 [CTLA4], and interleukin 2 receptor subunit alpha [IL2RA]) or thyroid-specific genes (thyroid stimulating hormone receptor [TSHR] and forkhead box E1 [FOXE1]) have been identified. Regarding thyroid function, many susceptibility loci for levels of TSH and free T4 have been identified through genome-wide analyses. In GWAS of differentiated thyroid cancer, associations at FOXE1, MAP3K12 binding inhibitory protein 1 (MBIP)-NK2 homeobox 1 (NKX2-1), disrupted in renal carcinoma 3 (DIRC3), neuregulin 1 (NRG1), and pecanex-like 2 (PCNXL2) have been commonly identified in people of European and Korean ancestry, and many other susceptibility loci have been found in specific populations. Through GWAS of various thyroid-related phenotypes, many susceptibility loci have been found, providing insights into the pathogenesis of thyroid diseases and disease co-clustering within families and individuals.

Tumor-suppressive function of UNC5D in papillary thyroid cancer

Background

Studies have shown an association of the UNC5D gene with kidney and bladder cancer and neuroblastoma. We investigated whether UNC5D acts as a tumor suppressor in papillary thyroid carcinoma (PTC).

Methods

Primary PTC tumors and matched normal thyroid tissues were obtained from 112 patients to detect UNC5D mRNA by real-time PCR. Genomic DNA sequencing was performed to detect BRAF mutation in PTC tumors. The association between UNC5D expression and clinicopathological data from PTC patients was reviewed retrospectively. PTC-derived cancer cell lines TPC-1 and K1 with stable transfection of UNC5D were used to investigate the functions of UNC5D. Flow cytometry, CCK-8, Transwell assay and scratch tests were used to examine cell cycle distribution, proliferation and migration.

Results

The expression of UNC5D was significantly decreased in PTC compared with adjacent normal thyroid tissues. Lower UNC5D expression was significantly associated with aggressive tumor behaviors, such as lymph node metastasis and BRAF mutation. Overexpression of UNC5D significantly suppressed malignant cell behaviors, including cell proliferation and migration, as well as tumor growth in vivo.

Conclusions

These findings suggest a potential tumor suppressor role of UNC5D in PTC progression; and provide insight into potential clinical relevance for the prognosis of PTC.

Genetics of thyroid function

Recent studies show that subtle variations in thyroid function, including subclinical thyroid dysfunction, and even variation in thyroid function within the normal range, are associated with morbidity and mortality. It is estimated that 40-65% of the inter-individual variation in serum TSH and FT4 levels is determined by genetic factors. To identify these factors, various linkage and candidate gene studies have been performed in the past, which have identified only a few genes. In the last decade, genome-wide association studies identified many new genes, while recent whole-genome sequencing efforts have also been proven to be effective. In the current review, we provide a systematic overview of these studies, including strengths and limitations. We discuss new techniques which will further clarify the genetic basis of thyroid function in the near future, as well as the potential use of these genetic markers in personalizing the management of thyroid disease patients.

New insights into thyroid hormone action

Thyroid hormones (TH) are endocrine messengers essential for normal development and function of virtually every vertebrate. The hypothalamic-pituitary-thyroid axis is exquisitely modulated to maintain nearly constant TH (T4 and T3) levels in circulation. However peripheral tissues and the CNS control the intracellular availability of TH, suggesting that circulating concentrations of TH are not fully representative of what each cell type sees. Indeed, recent work in the field has identified that TH transporters, deiodinases and thyroid hormone receptor coregulators can strongly control tissue-specific sensitivity to a set amount of TH. Furthermore, the mechanism by which the thyroid hormone receptors regulate target gene expression can vary by gene, tissue and cellular context. This review will highlight novel insights into the machinery that controls the cellular response to TH, which include unique signaling cascades. These findings shed new light into the pathophysiology of human diseases caused by abnormal TH signaling.

Carriers of a VEGFA enhancer polymorphism selectively binding CHOP/DDIT3 are predisposed to increased circulating levels of thyroid-stimulating hormone

BACKGROUND

Levels of serum thyroid-stimulating hormone (TSH) indicate thyroid function, because thyroid hormone negatively controls TSH release. Genetic variants in the vascular endothelial growth factor A (VEGFA) gene are associated with TSH levels. The aim of this study was to characterise the association of VEGFA variants with TSH in a Danish cohort and to identify and characterise functional variants.

METHODS

We performed an association study of the VEGFA locus for circulating TSH levels in 8445 Danish individuals. Lead variants were tested for allele-specific effects in vitro using luciferase reporter and gel-shift assays.

RESULTS

Four SNPs in VEGFA were associated with circulating TSH (rs9472138, rs881858, rs943080 and rs4711751). For rs881858, the presence of each G-allele was associated with a corresponding decrease in TSH levels of 2.3% (p=8.4×10^-9) and an increase in circulating free T4 levels (p=0.0014). The SNP rs881858 is located in a binding site for CHOP (C/EBP homology protein) and c/EBPβ (ccaat enhancer binding protein β). Reporter-gene analysis showed increased basal enhancer activity of the rs881858 A-allele versus the G-allele (34.5±9.9% (average±SEM), p=0.0012), while co-expression of CHOP effectively suppressed the rs881858 A-allele activity. The A-allele showed stronger binding to CHOP in gel-shift assays.

CONCLUSIONS

VEGF is an important angiogenic signal required for tissue expansion. We show that VEGFA variation giving allele-specific response to transcription factors with overlapping binding sites associate closely with circulating TSH levels. Because CHOP is induced by several types of intracellular stress, this indicates that cellular stress could be involved in the normal or pathophysiological response of the thyroid to TSH.

TRIAL REGISTRATION NUMBER

NCT00289237, NCT00316667; Results.

A single nucleotide polymorphism associated with isolated cleft lip and palate, thyroid cancer and hypothyroidism alters the activity of an oral epithelium and thyroid enhancer near FOXE1

Three common diseases, isolated cleft lip and cleft palate (CLP), hypothyroidism and thyroid cancer all map to the FOXE1 locus, but causative variants have yet to be identified. In patients with CLP, the frequency of coding mutations in FOXE1 fails to account for the risk attributable to this locus, suggesting that the common risk alleles reside in nearby regulatory elements. Using a combination of zebrafish and mouse transgenesis, we screened 15 conserved non-coding sequences for enhancer activity, identifying three that regulate expression in a tissue specific pattern consistent with endogenous foxe1 expression. These three, located -82.4, -67.7 and +22.6 kb from the FOXE1 start codon, are all active in the oral epithelium or branchial arches. The -67.7 and +22.6 kb elements are also active in the developing heart, and the -67.7 kb element uniquely directs expression in the developing thyroid. Within the -67.7 kb element is the SNP rs7850258 that is associated with all three diseases. Quantitative reporter assays in oral epithelial and thyroid cell lines show that the rs7850258 allele (G) associated with CLP and hypothyroidism has significantly greater enhancer activity than the allele associated with thyroid cancer (A). Moreover, consistent with predicted transcription factor binding differences, the -67.7 kb element containing rs7850258 allele G is significantly more responsive to both MYC and ARNT than allele A. By demonstrating that this common non-coding variant alters FOXE1 expression, we have identified at least in part the functional basis for the genetic risk of these seemingly disparate disorders.

Genetic determination of the hypothalamic-pituitary-thyroid axis: where do we stand?

For a long time it has been known that both hypo- and hyperthyroidism are associated with an increased risk of morbidity and mortality. In recent years, it has also become clear that minor variations in thyroid function, including subclinical dysfunction and variation in thyroid function within the reference range, can have important effects on clinical endpoints, such as bone mineral density, depression, metabolic syndrome, and cardiovascular mortality. Serum thyroid parameters show substantial interindividual variability, whereas the intraindividual variability lies within a narrow range. This suggests that every individual has a unique hypothalamus-pituitary-thyroid axis setpoint that is mainly determined by genetic factors, and this heritability has been estimated to be 40-60%. Various mutations in thyroid hormone pathway genes have been identified in persons with thyroid dysfunction or altered thyroid function tests. Because these causes are rare, many candidate gene and linkage studies have been performed over the years to identify more common variants (polymorphisms) associated with thyroid (dys)function, but only a limited number of consistent associations have been found. However, in the past 5 years, advances in genetic research have led to the identification of a large number of new candidate genes. In this review, we provide an overview of the current knowledge about the polygenic basis of thyroid (dys)function. This includes new candidate genes identified by genome-wide approaches, what insights these genes provide into the genetic basis of thyroid (dys)function, and which new techniques will help to further decipher the genetic basis of thyroid (dys)function in the near future.

Association between rs12045440 Polymorphism in the CAPZB Intron and Serum TSH Concentrations in Chinese Thyroid Tumor Patients

The aim of this study was to investigate the possible influence of different genotypes of the lead single nucleotide polymorphisms (SNPs) rs10917468 and rs12045440 in the CAPZB gene on the thyroid function in papillary thyroid carcinoma (PTC) and benign thyroid neoplasm (BN) patients. In the study, a significant association was detected between rs12045440 and serum TSH concentrations in thyroid tumor patients (p = 0.001). After the adjustment of relevant covariates, the difference between the mean serum TSH levels in different genotypes of rs12045440 was still significant in the BN group (p = 0.003) but was not significant in the PTC cases (p = 0.115). No significant association of rs10917468 with TSH levels was found. The SNP rs12045440 was associated with the serum TSH concentrations in Chinese thyroid tumor patients, especially in benign thyroid tumor cases.