Copy-Number Variation of the Glucose Transporter Gene SLC2A3 and Congenital Heart Defects in the 22q11.2 Deletion Syndrome

The 22q11.2 deletion syndrome (22q11DS; velocardiofacial/DiGeorge syndrome; VCFS/DGS) is the most common microdeletion syndrome and the phenotypic presentation is highly variable. Approximately 65% of individuals with 22q11DS have a congenital heart defect (CHD), mostly of the conotruncal type, and/or an aortic arch defect. The etiology of this phenotypic variability is not currently known. We hypothesized that copy-number variants (CNVs) outside the 22q11.2 deleted region might increase the risk of being born with a CHD in this sensitized population. Genotyping with Affymetrix SNP Array 6.0 was performed on two groups of subjects with 22q11DS separated by time of ascertainment and processing. CNV analysis was completed on a total of 949 subjects (cohort 1, n = 562; cohort 2, n = 387), 603 with CHDs (cohort 1, n = 363; cohort 2, n = 240) and 346 with normal cardiac anatomy (cohort 1, n = 199; cohort 2, n = 147). Our analysis revealed that a duplication of SLC2A3 was the most frequent CNV identified in the first cohort. It was present in 18 subjects with CHDs and 1 subject without (p = 3.12 × 10⁻³, two-tailed Fisher’s exact test). In the second cohort, the SLC2A3 duplication was also significantly enriched in subjects with CHDs (p = 3.30 × 10⁻², two-tailed Fisher’s exact test). The SLC2A3 duplication was the most frequent CNV detected and the only significant finding in our combined analysis (p = 2.68 × 10⁻⁴, two-tailed Fisher’s exact test), indicating that the SLC2A3 duplication might serve as a genetic modifier of CHDs and/or aortic arch anomalies in individuals with 22q11DS.     

Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

Genome-wide Association Study Identifies TNFSF15 and POU2AF1 as Susceptibility Loci for Primary Biliary Cirrhosis in the Japanese Population

For the identification of susceptibility loci for primary biliary cirrhosis (PBC), a genome-wide association study (GWAS) was performed in 963 Japanese individuals (487 PBC cases and 476 healthy controls) and in a subsequent replication study that included 1,402 other Japanese individuals (787 cases and 615 controls). In addition to the most significant susceptibility region, human leukocyte antigen (HLA), we identified two significant susceptibility loci, TNFSF15 (rs4979462) and POU2AF1 (rs4938534) (combined odds ratio [OR] = 1.56, p = 2.84 × 10⁻¹⁴ for rs4979462, and combined OR = 1.39, p = 2.38 × 10⁻⁸ for rs4938534). Among 21 non-HLA susceptibility loci for PBC identified in GWASs of individuals of European descent, three loci (IL7R, IKZF3, and CD80) showed significant associations (combined p = 3.66 × 10⁻⁸, 3.66 × 10⁻⁹, and 3.04 × 10⁻⁹, respectively) and STAT4 and NFKB1 loci showed suggestive association with PBC (combined p = 1.11 × 10⁻⁶ and 1.42 × 10⁻⁷, respectively) in the Japanese population. These observations indicated the existence of ethnic differences in genetic susceptibility loci to PBC and the importance of TNF signaling and B cell differentiation for the development of PBC in individuals of European descent and Japanese individuals.     

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10⁻⁴⁰, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1^∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10⁻⁴³) and HLA-DQα1 47 (p = 4.02 × 10⁻⁴⁶), 56, and 76 (both p = 1.84 × 10⁻⁴⁵) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10⁻⁶, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10⁻⁶, OR = 1.20), and REL (rs115674477, p = 1.10 × 10⁻⁵, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.     

Genetic Variants at 6p21.1 and 7p15.3 Are Associated with Risk of Multiple Cancers in Han Chinese

Cancer susceptibility loci identified in reported genome-wide association studies (GWAS) are often tumor-specific; however, evidence of pleiotropy of some genes/loci has also been observed and biologically plausible. We hypothesized that there are important regions in the genome harboring genetic variants associated with risk of multiple types of cancer. In the current study, we attempted to map genetic variants that have consistent effects on risk of multiple cancers using our existing genome-wide scan data of lung cancer, noncardia gastric cancer, and esophageal squamous-cell carcinoma with overall 5,368 cases and 4,006 controls (GWAS stage), followed by a further evaluation in additional 9,001 cases with one of these cancer types and 11,436 controls (replication stage). Five variants satisfying the criteria of pleiotropy with p values from 1.10 × 10⁻⁸ to 8.96 × 10⁻⁶ for genome-wide scans of three cancer types were further evaluated in the replication stage. We found consistent associations of rs2494938 at 6p21.1 and rs2285947 at 7p15.3 with these three cancers in both GWAS and replication stages. In combined samples of GWAS and replication stages, the minor alleles of rs2494938 and rs2285947 were significantly associated with an increased risk of the cancers (odds ratio [OR] = 1.15, 95% confidence interval [CI], 1.10–1.19 and OR = 1.17, 95% CI, 1.12–1.21), with the p values being 1.20 × 10⁻¹² and 1.26 × 10⁻¹⁶, respectively, which are at a genome-wide significance level. Our findings highlight the potential importance of variants at 6p21.1 and 7p15.3 in the susceptibility to multiple cancers.     

HLA-DPB1 and HLA Class I Confer Risk of and Protection from Narcolepsy

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1^∗01:02-DQB1^∗06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1^∗01:03-DPB1^∗04:02 (DP0402; odds ratio [OR] = 0.51 [0.38–0.67], p = 1.01 × 10⁻⁶) and HLA-DPA1^∗01:03-DPB1^∗04:01 (DP0401; OR = 0.61 [0.47–0.80], p = 2.07 × 10⁻⁴) and predisposing effects of HLA-DPB1^∗05:01 in Asians (OR = 1.76 [1.34–2.31], p = 4.71 × 10⁻⁰⁵). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38–0.55] p = 8.99 × 10⁻¹⁷) and DP0501 (OR = 1.38 [1.18–1.61], p = 7.11 × 10⁻⁵). HLA-class-II-independent associations with HLA-A^∗11:01 (OR = 1.32 [1.13–1.54], p = 4.92 × 10⁻⁴), HLA-B^∗35:03 (OR = 1.96 [1.41–2.70], p = 5.14 × 10⁻⁵), and HLA-B^∗51:01 (OR = 1.49 [1.25–1.78], p = 1.09 × 10⁻⁵) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.     

Fine Mapping Seronegative and Seropositive Rheumatoid Arthritis to Shared and Distinct HLA Alleles by Adjusting for the Effects of Heterogeneity

Despite progress in defining human leukocyte antigen (HLA) alleles for anti-citrullinated-protein-autoantibody-positive (ACPA⁺) rheumatoid arthritis (RA), identifying HLA alleles for ACPA-negative (ACPA⁻) RA has been challenging because of clinical heterogeneity within clinical cohorts. We imputed 8,961 classical HLA alleles, amino acids, and SNPs from Immunochip data in a discovery set of 2,406 ACPA⁻ RA case and 13,930 control individuals. We developed a statistical approach to identify and adjust for clinical heterogeneity within ACPA⁻ RA and observed independent associations for serine and leucine at position 11 in HLA-DRβ1 (p = 1.4 × 10⁻¹³, odds ratio [OR] = 1.30) and for aspartate at position 9 in HLA-B (p = 2.7 × 10⁻¹², OR = 1.39) within the peptide binding grooves. These amino acid positions induced associations at HLA-DRB1^∗03 (encoding serine at 11) and HLA-B^∗08 (encoding aspartate at 9). We validated these findings in an independent set of 427 ACPA⁻ case subjects, carefully phenotyped with a highly sensitive ACPA assay, and 1,691 control subjects (HLA-DRβ1 Ser11+Leu11: p = 5.8 × 10⁻⁴, OR = 1.28; HLA-B Asp9: p = 2.6 × 10⁻³, OR = 1.34). Although both amino acid sites drove risk of ACPA⁺ and ACPA⁻ disease, the effects of individual residues at HLA-DRβ1 position 11 were distinct (p < 2.9 × 10⁻¹⁰⁷). We also identified an association with ACPA⁺ RA at HLA-A position 77 (p = 2.7 × 10⁻⁸, OR = 0.85) in 7,279 ACPA⁺ RA case and 15,870 control subjects. These results contribute to mounting evidence that ACPA⁺ and ACPA⁻ RA are genetically distinct and potentially have separate autoantigens contributing to pathogenesis. We expect that our approach might have broad applications in analyzing clinical conditions with heterogeneity at both major histocompatibility complex (MHC) and non-MHC regions.     

Identification of Multiple Genetic Susceptibility Loci in Takayasu Arteritis

Takayasu arteritis is a rare inflammatory disease of large arteries. The etiology of Takayasu arteritis remains poorly understood, but genetic contribution to the disease pathogenesis is supported by the genetic association with HLA-B^∗52. We genotyped ∼200,000 genetic variants in two ethnically divergent Takayasu arteritis cohorts from Turkey and North America by using a custom-designed genotyping platform (Immunochip). Additional genetic variants and the classical HLA alleles were imputed and analyzed. We identified and confirmed two independent susceptibility loci within the HLA region (r² < 0.2): HLA-B/MICA (rs12524487, OR = 3.29, p = 5.57 × 10⁻¹⁶) and HLA-DQB1/HLA-DRB1 (rs113452171, OR = 2.34, p = 3.74 × 10⁻⁹; and rs189754752, OR = 2.47, p = 4.22 × 10⁻⁹). In addition, we identified and confirmed a genetic association between Takayasu arteritis and the FCGR2A/FCGR3A locus on chromosome 1 (rs10919543, OR = 1.81, p = 5.89 × 10⁻¹²). The risk allele in this locus results in increased mRNA expression of FCGR2A. We also established the genetic association between IL12B and Takayasu arteritis (rs56167332, OR = 1.54, p = 2.18 × 10⁻⁸).     

APLP2 Regulates Refractive Error and Myopia Development in Mice and Humans

Myopia is the most common vision disorder and the leading cause of visual impairment worldwide. However, gene variants identified to date explain less than 10% of the variance in refractive error, leaving the majority of heritability unexplained (“missing heritability”). Previously, we reported that expression of APLP2 was strongly associated with myopia in a primate model. Here, we found that low-frequency variants near the 5’-end of APLP2 were associated with refractive error in a prospective UK birth cohort (n = 3,819 children; top SNP rs188663068, p = 5.0 × 10⁻⁴) and a CREAM consortium panel (n = 45,756 adults; top SNP rs7127037, p = 6.6 × 10⁻³). These variants showed evidence of differential effect on childhood longitudinal refractive error trajectories depending on time spent reading (gene x time spent reading x age interaction, p = 4.0 × 10⁻³). Furthermore, Aplp2 knockout mice developed high degrees of hyperopia (+11.5 ± 2.2 D, p < 1.0 × 10⁻⁴) compared to both heterozygous (-0.8 ± 2.0 D, p < 1.0 × 10⁻⁴) and wild-type (+0.3 ± 2.2 D, p < 1.0 × 10⁻⁴) littermates and exhibited a dose-dependent reduction in susceptibility to environmentally induced myopia (F(2, 33) = 191.0, p < 1.0 × 10⁻⁴). This phenotype was associated with reduced contrast sensitivity (F(12, 120) = 3.6, p = 1.5 × 10⁻⁴) and changes in the electrophysiological properties of retinal amacrine cells, which expressed Aplp2. This work identifies APLP2 as one of the “missing” myopia genes, demonstrating the importance of a low-frequency gene variant in the development of human myopia. It also demonstrates an important role for APLP2 in refractive development in mice and humans, suggesting a high level of evolutionary conservation of the signaling pathways underlying refractive eye development.     

Genome-Wide Association Study of Celiac Disease in North America Confirms FRMD4B as New Celiac Locus

We performed a genome-wide association study (GWAS) of 1550 North American celiac disease cases and 3084 controls. Twelve SNPs, distributed across four regions (3p21.31, 4q27, 6q15, 6q25), were significantly associated with disease (p-value <1.0×10⁻⁷), and a further seven SNPs, across four additional regions (1q24.3, 10p15.1, 6q22.31, 17q21.32) had suggestive evidence (1.0×10⁻⁷ < p-value < 1.0×10⁻⁶). This study replicated a previous suggestive association within FRMD4B (3p14.1), confirming it as a celiac disease locus. All four regions with significant associations and two regions with suggestive results (1q24.3, 10p15.1) were known disease loci. The 6q22.31 and 10p11.23 regions were not replicated. A total of 410 SNPs distributed across the eight significant and suggestive regions were tested for association with dermatitis herpetiformis and microscopic colitis. Preliminary, suggestive statistical evidence for association with the two traits was found at chromosomes 3p21.31, 6q15, 6q25, 1q24.3 and 10p11.23, with future studies being required to validate the reported associations.     

A Higher Mutational Burden in Females Supports a “Female Protective Model” in Neurodevelopmental Disorders

Increased male prevalence has been repeatedly reported in several neurodevelopmental disorders (NDs), leading to the concept of a “female protective model.” We investigated the molecular basis of this sex-based difference in liability and demonstrated an excess of deleterious autosomal copy-number variants (CNVs) in females compared to males (odds ratio [OR] = 1.46, p = 8 × 10⁻¹⁰) in a cohort of 15,585 probands ascertained for NDs. In an independent autism spectrum disorder (ASD) cohort of 762 families, we found a 3-fold increase in deleterious autosomal CNVs (p = 7 × 10⁻⁴) and an excess of private deleterious single-nucleotide variants (SNVs) in female compared to male probands (OR = 1.34, p = 0.03). We also showed that the deleteriousness of autosomal SNVs was significantly higher in female probands (p = 0.0006). A similar bias was observed in parents of probands ascertained for NDs. Deleterious CNVs (>400 kb) were maternally inherited more often (up to 64%, p = 10⁻¹⁵) than small CNVs < 400 kb (OR = 1.45, p = 0.0003). In the ASD cohort, increased maternal transmission was also observed for deleterious CNVs and SNVs. Although ASD females showed higher mutational burden and lower cognition, the excess mutational burden remained, even after adjustment for those cognitive differences. These results strongly suggest that females have an increased etiological burden unlinked to rare deleterious variants on the X chromosome. Carefully phenotyped and genotyped cohorts will be required for identifying the symptoms, which show gender-specific liability to mutational burden.     

A Genome-Wide Association Study of Optic Disc Parameters

The optic nerve head is involved in many ophthalmic disorders, including common diseases such as myopia and open-angle glaucoma. Two of the most important parameters are the size of the optic disc area and the vertical cup-disc ratio (VCDR). Both are highly heritable but genetically largely undetermined. We performed a meta-analysis of genome-wide association (GWA) data to identify genetic variants associated with optic disc area and VCDR. The gene discovery included 7,360 unrelated individuals from the population-based Rotterdam Study I and Rotterdam Study II cohorts. These cohorts revealed two genome-wide significant loci for optic disc area, rs1192415 on chromosome 1p22 (p = 6.72×10⁻¹⁹) within 117 kb of the CDC7 gene and rs1900004 on chromosome 10q21.3-q22.1 (p = 2.67×10⁻³³) within 10 kb of the ATOH7 gene. They revealed two genome-wide significant loci for VCDR, rs1063192 on chromosome 9p21 (p = 6.15×10⁻¹¹) in the CDKN2B gene and rs10483727 on chromosome 14q22.3-q23 (p = 2.93×10⁻¹⁰) within 40 kbp of the SIX1 gene. Findings were replicated in two independent Dutch cohorts (Rotterdam Study III and Erasmus Rucphen Family study; N = 3,612), and the TwinsUK cohort (N = 843). Meta-analysis with the replication cohorts confirmed the four loci and revealed a third locus at 16q12.1 associated with optic disc area, and four other loci at 11q13, 13q13, 17q23 (borderline significant), and 22q12.1 for VCDR. ATOH7 was also associated with VCDR independent of optic disc area. Three of the loci were marginally associated with open-angle glaucoma. The protein pathways in which the loci of optic disc area are involved overlap with those identified for VCDR, suggesting a common genetic origin.     

Runs of Homozygosity: Association with Coronary Artery Disease and Gene Expression in Monocytes and Macrophages

Runs of homozygosity (ROHs) are recognized signature of recessive inheritance. Contributions of ROHs to the genetic architecture of coronary artery disease and regulation of gene expression in cells relevant to atherosclerosis are not known. Our combined analysis of 24,320 individuals from 11 populations of white European ethnicity showed an association between coronary artery disease and both the count and the size of ROHs. Individuals with coronary artery disease had approximately 0.63 (95% CI: 0.4–0.8) excess of ROHs when compared to coronary-artery-disease-free control subjects (p = 1.49 × 10⁻⁹). The average total length of ROHs was approximately 1,046.92 (95% CI: 634.4–1,459.5) kb greater in individuals with coronary artery disease than control subjects (p = 6.61 × 10⁻⁷). None of the identified individual ROHs was associated with coronary artery disease after correction for multiple testing. However, in aggregate burden analysis, ROHs favoring increased risk of coronary artery disease were much more common than those showing the opposite direction of association with coronary artery disease (p = 2.69 × 10⁻³³). Individual ROHs showed significant associations with monocyte and macrophage expression of genes in their close proximity—subjects with several individual ROHs showed significant differences in the expression of 44 mRNAs in monocytes and 17 mRNAs in macrophages when compared to subjects without those ROHs. This study provides evidence for an excess of homozygosity in coronary artery disease in outbred populations and suggest the potential biological relevance of ROHs in cells of importance to the pathogenesis of atherosclerosis.     

Maternal Modifiers and Parent-of-Origin Bias of the Autism-Associated 16p11.2 CNV

Recurrent deletions and duplications at chromosomal region 16p11.2 are a major genetic contributor to autism but also associate with a wider range of pediatric diagnoses, including intellectual disability, coordination disorder, and language disorder. In order to investigate the potential genetic basis for phenotype variability, we assessed the parent of origin of the 16p11.2 copy-number variant (CNV) and the presence of additional CNVs in 126 families for which detailed phenotype data were available. Among de novo cases, we found a strong maternal bias for the origin of deletions (59/66, 89.4% of cases, p = 2.38 × 10⁻¹¹), the strongest such effect so far observed for a CNV associated with a microdeletion syndrome. In contrast to de novo events, we observed no transmission bias for inherited 16p11.2 CNVs, consistent with a female meiotic hotspot of unequal crossover driving this maternal bias. We analyzed this 16p11.2 CNV cohort for the presence of secondary CNVs and found a significant maternal transmission bias for secondary deletions (32 maternal versus 14 paternal, p = 1.14 × 10⁻²). Of the secondary deletions that disrupted a gene, 82% were either maternally inherited or de novo (p = 4.3 × 10⁻³). Nine probands carry secondary CNVs that disrupt genes associated with autism and/or intellectual disability risk variants. Our findings demonstrate a strong bias toward maternal origin of 16p11.2 de novo deletions as well as a maternal transmission bias for secondary deletions that contribute to the clinical outcome on a background sensitized by the 16p11.2 CNV.     

Fine Mapping on Chromosome 10q22-q23 Implicates Neuregulin 3 in Schizophrenia

Linkage studies have implicated 10q22-q23 as a schizophrenia (SZ) susceptibility locus in Ashkenazi Jewish (AJ) and Han Chinese from Taiwan populations. To further explore our previous linkage signal in the AJ population (NPL score: 4.27, empirical p = 2 × 10⁻⁵), we performed a peakwide association fine mapping study by using 1414 SNPs across ~12.5 Mb in 10q22-q23. We genotyped 1515 AJ individuals, including 285 parent-child trios, 173 unrelated cases, and 487 unrelated controls. We analyzed the binary diagnostic phenotype of SZ and 9 heritable quantitative traits derived from a principal components factor analysis of 73 items from our consensus diagnostic ratings and direct assessment interviews. Although no marker withstood multiple test correction for association with the binary SZ phenotype, we found strong evidence of association by using the “delusion” factor as the quantitative trait at three SNPs (rs10883866, rs10748842, and rs6584400) located in a 13 kb interval in intron 1 of Neuregulin 3 (NRG3). Our best p value from family-based association analysis was 7.26 × 10⁻⁷. We replicated this association in the collection of 173 unrelated AJ cases (p = 1.55 × 10⁻²), with a combined p value of 2.30 × 10⁻⁷. After performing 10,000 permutations of each of the phenotypes, we estimated the empirical study-wide significance across all 9 factors (90,000 permutations) to be p = 2.7 × 10⁻³. NRG3 is primarily expressed in the central nervous system and is one of three paralogs of NRG1, a gene strongly implicated in SZ. These biological properties together with our linkage and association results strongly support NRG3 as a gene involved in SZ.     

A Genome-Wide Association Study Uncovers a Genetic Locus Associated with Thoracic-to-Hip Ratio in Koreans

The thoracic-to-hip circumference ratio (THR) is an anthropometric marker recently described as a predictor of type 2 diabetes. In this study, we performed a genome-wide association study (GWAS) followed by confirmatory analyses to identify genetic markers associated with THR. A total of 7,240 Korean subjects (4,988 for the discovery stage and 2,252 for the confirmatory analyses) were recruited for this study, and genome-wide single nucleotide polymorphism (SNP) genotyping of the initial 4,988 individuals was performed using Affymetrix Human SNP array 5.0. Linear regression analysis was then performed to adjust for the effects of age, sex, and current diabetes medication status on the THR of the study subjects. In the initial discovery stage, there was a statistically nominal association between minor alleles of SNP markers on chromosomes 4, 8, 10, and 12, and THR changes (p < 5.0 × 10⁻⁶). The subsequent confirmatory analyses of these markers, however, only detected a significant association between two SNPs in the HECTD4 gene and decreased THRs. Notably, this association was detected in male (rs11066280: p = 1.14 × 10⁻²; rs2074356: p = 1.10 × 10⁻²), but not in female subjects. Meanwhile, the combined results from the two analyses (initial and confirmatory) indicated that minor alleles of these two intronic variants exhibited a significant genome-wide association with decreased THR in the male subjects (n = 3,155; rs11066280: effect size = −0.008624, p = 6.19 × 10⁻⁹; rs2074356: effect size = −0.008762, p = 1.89 × 10⁻⁸). Furthermore, minor alleles of these two SNPs exhibited protective effects on patients’ risks for developing type 2 diabetes. In conclusion, we have identified two genetic variations in HECTD4 that are associated with THR, particularly in men.