A Genome Wide Association Study of Mathematical Ability Reveals an Association at Chromosome 3q29, a Locus Associated with Autism and Learning Difficulties: A Preliminary Study

Mathematical ability is heritable, but few studies have directly investigated its molecular genetic basis. Here we aimed to identify specific genetic contributions to variation in mathematical ability. We carried out a genome wide association scan using pooled DNA in two groups of U.K. samples, based on end of secondary/high school national academic exam achievement: high (n = 419) versus low (n = 183) mathematical ability while controlling for their verbal ability. Significant differences in allele frequencies between these groups were searched for in 906,600 SNPs using the Affymetrix GeneChip Human Mapping version 6.0 array. After meeting a threshold of p<1.5×10⁻⁵, 12 SNPs from the pooled association analysis were individually genotyped in 542 of the participants and analyzed to validate the initial associations (lowest p-value 1.14 ×10⁻⁶). In this analysis, one of the SNPs (rs789859) showed significant association after Bonferroni correction, and four (rs10873824, rs4144887, rs12130910 rs2809115) were nominally significant (lowest p-value 3.278 × 10⁻⁴). Three of the SNPs of interest are located within, or near to, known genes (FAM43A, SFT2D1, C14orf64). The SNP that showed the strongest association, rs789859, is located in a region on chromosome 3q29 that has been previously linked to learning difficulties and autism. rs789859 lies 1.3 kbp downstream of LSG1, and 700 bp upstream of FAM43A, mapping within the potential promoter/regulatory region of the latter. To our knowledge, this is only the second study to investigate the association of genetic variants with mathematical ability, and it highlights a number of interesting markers for future study.     

Variants in TF and HFE Explain ∼40% of Genetic Variation in Serum-Transferrin Levels

Only a small proportion of genetic variation in complex traits has been explained by SNPs from genome-wide association studies (GWASs). We report the results from two GWASs for serum markers of iron status (serum iron, serum transferrin, transferrin saturation with iron, and serum ferritin), which are important in iron overload (e.g., hemochromatosis) and deficiency (e.g., anemia) conditions. We performed two GWASs on samples of Australians of European descent. In the first GWAS, 411 adolescent twins and their siblings were genotyped with 100K SNPs. rs1830084, 10.8 kb 3′ of TF, was significantly associated with serum transferrin (p total association test = 1.0 × 10⁻⁹; p within-family test = 2.2 × 10⁻⁵). In the second GWAS on an independent sample of 459 female monozygotic (MZ) twin pairs genotyped with 300K SNPs, we found rs3811647 (within intron 11 of TF, HapMap CEU r² with rs1830084 = 0.86) was significantly associated with serum transferrin (p = 3.0 × 10⁻¹⁵). In the second GWAS, we found two additional and independent SNPs on TF (rs1799852 and rs2280673) and confirmed the known C282Y mutation in HFE to be independently associated with serum transferrin. The three variants in TF (rs3811647, rs1799852 and rs2280673) plus the HFE C282Y mutation explained ~40% of genetic variation in serum transferrin (p = 7.8 × 10⁻²⁵). These findings are potentially important for our understanding of iron metabolism and of regulation of hepatic protein secretion, and also strongly support the hypothesis that the genetic architecture of some endophenotypes may be simpler than that of disease.     

Genetic Analysis Identifies DDR2 as a Novel Gene Affecting Bone Mineral Density and Osteoporotic Fractures in Chinese Population

DDR2 gene, playing an essential role in regulating osteoblast differentiation and chondrocyte maturation, may influence bone mineral density (BMD) and osteoporosis, but the genetic variations actually leading to the association remain to be elucidated. Therefore, the aim of this study was to investigate whether the genetic variants in DDR2 are associated with BMD and fracture risk. This study was performed in three samples from two ethnicities, including 1,300 Chinese Han subjects, 700 Chinese Han subjects (350 with osteoporotic hip fractures and 350 healthy controls) and 2,286 US white subjects. Twenty-eight SNPs in DDR2 were genotyped and tested for associations with hip BMD and fractures. We identified 3 SNPs in DDR2 significantly associated with hip BMD in the Chinese population after multiple testing adjustments, which were rs7521233 (P = 1.06×10⁻⁴, β: −0.018 for allele C), rs7553831 (P = 1.30×10⁻⁴, β: −0.018 for allele T), and rs6697469 (P = 1.59×10⁻³, β: −0.015 for allele C), separately. These three SNPs were in high linkage disequilibrium. Haplotype analyses detected two significantly associated haplotypes, including one haplotype in block 2 (P = 9.54×10⁻⁴, β: −0.016) where these three SNPs located. SNP rs6697469 was also associated with hip fractures (P = 0.043, OR: 1.42) in the Chinese population. The effect on fracture risk was consistent with its association with lower BMD. However, in the white population, we didn’t observe significant associations with hip BMD. eQTL analyses revealed that SNPs associated with BMD also affected DDR2 mRNA expression levels in Chinese. Our findings, together with the prior biological evidence, suggest that DDR2 could be a new candidate for osteoporosis in Chinese population. Our results also reveal an ethnic difference, which highlights the need for further genetic studies in each ethnic group.     

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.     

Genome-Wide Association Study for Serum Complement C3 and C4 Levels in Healthy Chinese Subjects

Complement C3 and C4 play key roles in the main physiological activities of complement system, and their deficiencies or over-expression are associated with many clinical infectious or immunity diseases. A two-stage genome-wide association study (GWAS) was performed for serum levels of C3 and C4. The first stage was conducted in 1,999 healthy Chinese men, and the second stage was performed in an additional 1,496 subjects. We identified two SNPs, rs3753394 in CFH gene and rs3745567 in C3 gene, that are significantly associated with serum C3 levels at a genome-wide significance level (P = 7.33×10⁻¹¹ and P = 1.83×10⁻⁹, respectively). For C4, one large genomic region on chromosome 6p21.3 is significantly associated with serum C4 levels. Two SNPs (rs1052693 and rs11575839) were located in the MHC class I area that include HLA-A, HLA-C, and HLA-B genes. Two SNPs (rs2075799 and rs2857009) were located 5′ and 3′ of C4 gene. The other four SNPs, rs2071278, rs3763317, rs9276606, and rs241428, were located in the MHC class II region that includes HLA-DRA, HLA-DRB, and HLA-DQB genes. The combined P-values for those eight SNPs ranged from 3.19×10⁻²² to 5.62×10⁻⁹⁷. HBsAg-positive subjects have significantly lower C3 and C4 protein concentrations compared with HBsAg-negative subjects (P<0.05). Our study is the first GWAS report which shows genetic components influence the levels of complement C3 and C4. Our significant findings provide novel insights of their related autoimmune, infectious diseases, and molecular mechanisms.     

Endometriosis Is Associated with Rare Copy Number Variants

Endometriosis is a complex gynecological condition that affects 6–10% of women in their reproductive years and is defined by the presence of endometrial glands and stroma outside the uterus. Twin, family, and genome-wide association (GWA) studies have confirmed a genetic role, yet only a small part of the genetic risk can be explained by SNP variation. Copy number variants (CNVs) account for a greater portion of human genetic variation than SNPs and include more recent mutations of large effect. CNVs, likely to be prominent in conditions with decreased reproductive fitness, have not previously been examined as a genetic contributor to endometriosis. Here we employ a high-density genotyping microarray in a genome-wide survey of CNVs in a case-control population that includes 2,126 surgically confirmed endometriosis cases and 17,974 population controls of European ancestry. We apply stringent quality filters to reduce the false positive rate common to many CNV-detection algorithms from 77.7% to 7.3% without noticeable reduction in the true positive rate. We detected no differences in the CNV landscape between cases and controls on the global level which showed an average of 1.92 CNVs per individual with an average size of 142.3 kb. On the local level we identify 22 CNV-regions at the nominal significance threshold (P<0.05), which is greater than the 8.15 CNV-regions expected based on permutation analysis (P<0.001). Three CNV''s passed a genome-wide P-value threshold of 9.3×10⁻⁴; a deletion at SGCZ on 8p22 (P = 7.3×10⁻⁴, OR = 8.5, Cl = 2.3–31.7), a deletion in MALRD1 on 10p12.31 (P = 5.6×10⁻⁴, OR = 14.1, Cl = 2.7–90.9), and a deletion at 11q14.1 (P = 5.7×10⁻⁴, OR = 33.8, Cl = 3.3–1651). Two SNPs within the 22 CNVRs show significant genotypic association with endometriosis after adjusting for multiple testing; rs758316 in DPP6 on 7q36.2 (P = 0.0045) and rs4837864 in ASTN2 on 9q33.1 (P = 0.0002). Together, the CNV-loci are detected in 6.9% of affected women compared to 2.1% in the general population.     

Systematic Evaluation of Genetic Variants for Polycystic Ovary Syndrome in a Chinese Population

To date, eleven genome-wide significant (GWS) loci (P < 5×10⁻⁸) for polycystic ovary syndrome (PCOS) have been identified through genome-wide association studies (GWAS). Some of the risk loci have been selected for replications and validated in multiple ethnicities, however, few previous studies investigated all loci. Scanning all the GWAS variants would demonstrate a more informative profile of variance they explained. Thus, we analyzed all the 17 single nucleotide polymorphisms (SNPs) mapping to the 11 GWAS loci in an independent sample set of 800 Chinese subjects with PCOS and 1110 healthy controls systematically. Variants of rs3802457 in C9orf3 locus (P = 5.99×10⁻⁴) and rs13405728 in LHCGR locus (P = 3.73×10⁻⁴) were significantly associated with PCOS after the strict Bonferroni correction in our data set. The further haplotype analysis indicated that in the block of C9orf3 gene (rs4385527 and rs3802457), GA haplotype played a protective role in PCOS (8.7 vs 5.0, P = 9.85×10⁻⁶, OR = 0.548, 95%CI = 0.418–0.717), while GG haplotype was found suffering from an extraordinarily increased risk of PCOS (73.6% vs79.2%, P = 3.41×10⁻⁵, OR = 1.394, 95%CI = 1.191–1.632). Moreover, the directions of effects for all SNPs were consistent with previous GWAS reports (P = 1.53×10⁻⁵). Polygenic score analysis demonstrated that these 17 SNPs have a significant capacity on predicting case-control status in our samples (P = 7.17×10⁻⁹), meanwhile all these gathered 17 SNPs explained about 2.40% of variance. Our findings supported that C9orf3 and LHCGR loci variants were vital susceptibility of PCOS.     

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.     

Susceptibility Loci Associated with Specific and Shared Subtypes of Lymphoid Malignancies

The genetics of lymphoma susceptibility reflect the marked heterogeneity of diseases that comprise this broad phenotype. However, multiple subtypes of lymphoma are observed in some families, suggesting shared pathways of genetic predisposition to these pathologically distinct entities. Using a two-stage GWAS, we tested 530,583 SNPs in 944 cases of lymphoma, including 282 familial cases, and 4,044 public shared controls, followed by genotyping of 50 SNPs in 1,245 cases and 2,596 controls. A novel region on 11q12.1 showed association with combined lymphoma (LYM) subtypes. SNPs in this region included rs12289961 near LPXN, (P_LYM = 3.89×10⁻⁸, OR = 1.29) and rs948562 (P_LYM = 5.85×10⁻⁷, OR = 1.29). A SNP in a novel non-HLA region on 6p23 (rs707824, P_NHL = 5.72×10⁻⁷) was suggestive of an association conferring susceptibility to lymphoma. Four SNPs, all in a previously reported HLA region, 6p21.32, showed genome-wide significant associations with follicular lymphoma. The most significant association with follicular lymphoma was for rs4530903 (P_FL = 2.69×10⁻¹², OR = 1.93). Three novel SNPs near the HLA locus, rs9268853, rs2647046, and rs2621416, demonstrated additional variation contributing toward genetic susceptibility to FL associated with this region. Genes implicated by GWAS were also found to be cis-eQTLs in lymphoblastoid cell lines; candidate genes in these regions have been implicated in hematopoiesis and immune function. These results, showing novel susceptibility regions and allelic heterogeneity, point to the existence of pathways of susceptibility to both shared as well as specific subtypes of lymphoid malignancy.     

Genome-Wide Association Study of Late-Onset Myasthenia Gravis: Confirmation of TNFRSF11A and Identification of ZBTB10 and Three Distinct HLA Associations

To investigate the genetics of late-onset myasthenia gravis (LOMG), we conducted a genome-wide association study imputation of >6 million single nucleotide polymorphisms (SNPs) in 532 LOMG cases (anti–acetylcholine receptor [AChR] antibody positive; onset age ≥50 years) and 2,128 controls matched for sex and population substructure. The data confirm reported TNFRSF11A associations (rs4574025, P = 3.9 × 10⁻⁷, odds ratio [OR] 1.42) and identify a novel candidate gene, ZBTB10, achieving genome-wide significance (rs6998967, P = 8.9 × 10⁻¹⁰, OR 0.53). Several other SNPs showed suggestive significance including rs2476601 (P = 6.5 × 10⁻⁶, OR 1.62) encoding the PTPN22 R620W variant noted in early-onset myasthenia gravis (EOMG) and other autoimmune diseases. In contrast, EOMG-associated SNPs in TNIP1 showed no association in LOMG, nor did other loci suggested for EOMG. Many SNPs within the major histocompatibility complex (MHC) region showed strong associations in LOMG, but with smaller effect sizes than in EOMG (highest OR ~2 versus ~6 in EOMG). Moreover, the strongest associations were in opposite directions from EOMG, including an OR of 0.54 for DQA1*05:01 in LOMG (P = 5.9 × 10⁻¹²) versus 2.82 in EOMG (P = 3.86 × 10⁻⁴⁵). Association and conditioning studies for the MHC region showed three distinct and largely independent association peaks for LOMG corresponding to (a) MHC class II (highest attenuation when conditioning on DQA1), (b) HLA-A and (c) MHC class III SNPs. Conditioning studies of human leukocyte antigen (HLA) amino acid residues also suggest potential functional correlates. Together, these findings emphasize the value of subgrouping myasthenia gravis patients for clinical and basic investigations and imply distinct predisposing mechanisms in LOMG.     

Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits

The obesity epidemic is responsible for a substantial economic burden in developed countries and is a major risk factor for type 2 diabetes and cardiovascular disease. The disease is the result not only of several environmental risk factors, but also of genetic predisposition. To take advantage of recent advances in gene-mapping technology, we executed a genome-wide association scan to identify genetic variants associated with obesity-related quantitative traits in the genetically isolated population of Sardinia. Initial analysis suggested that several SNPs in the FTO and PFKP genes were associated with increased BMI, hip circumference, and weight. Within the FTO gene, rs9930506 showed the strongest association with BMI (p = 8.6 ×10⁻⁷), hip circumference (p = 3.4 × 10⁻⁸), and weight (p = 9.1 × 10⁻⁷). In Sardinia, homozygotes for the rare “G” allele of this SNP (minor allele frequency = 0.46) were 1.3 BMI units heavier than homozygotes for the common “A” allele. Within the PFKP gene, rs6602024 showed very strong association with BMI (p = 4.9 × 10⁻⁶). Homozygotes for the rare “A” allele of this SNP (minor allele frequency = 0.12) were 1.8 BMI units heavier than homozygotes for the common “G” allele. To replicate our findings, we genotyped these two SNPs in the GenNet study. In European Americans (N = 1,496) and in Hispanic Americans (N = 839), we replicated significant association between rs9930506 in the FTO gene and BMI (p-value for meta-analysis of European American and Hispanic American follow-up samples, p = 0.001), weight (p = 0.001), and hip circumference (p = 0.0005). We did not replicate association between rs6602024 and obesity-related traits in the GenNet sample, although we found that in European Americans, Hispanic Americans, and African Americans, homozygotes for the rare “A” allele were, on average, 1.0–3.0 BMI units heavier than homozygotes for the more common “G” allele. In summary, we have completed a whole genome–association scan for three obesity-related quantitative traits and report that common genetic variants in the FTO gene are associated with substantial changes in BMI, hip circumference, and body weight. These changes could have a significant impact on the risk of obesity-related morbidity in the general population.     

Genome-wide Association Study Reveals Multiple Nasopharyngeal Carcinoma-Associated Loci within the HLA Region at Chromosome 6p21.3

Nasopharyngeal carcinoma (NPC) is a multifactorial malignancy closely associated with genetic factors and Epstein-Barr virus infection. To identify the common genetic variants linked to NPC susceptibility, we conducted a genome-wide association study (GWAS) in 277 NPC patients and 285 healthy controls within the Taiwanese population, analyzing 480,365 single-nucleotide polymorphisms (SNPs). Twelve statistically significant SNPs were identified and mapped to chromosome 6p21.3. Associations were replicated in two independent sets of case-control samples. Two of the most significant SNPs (rs2517713 and rs2975042; p_combined = 3.9 × 10⁻²⁰ and 1.6 × 10⁻¹⁹, respectively) were located in the HLA-A gene. Moreover, we detected significant associations between NPC and two genes: specifically, gamma aminobutyric acid b receptor 1 (GABBR1) (rs29232; p_combined = 8.97 × 10⁻¹⁷) and HLA-F (rs3129055 and rs9258122; p_combined = 7.36 × 10⁻¹¹ and 3.33 × 10⁻¹⁰, respectively). Notably, the association of rs29232 remained significant (residual p < 5 × 10⁻⁴) after adjustment for age, gender, and HLA-related SNPs. Furthermore, higher GABA_B receptor 1 expression levels can be found in the tumor cells in comparison to the adjacent epithelial cells (p < 0.001) in NPC biopsies, implying a biological role of GABBR1 in NPC carcinogenesis. To our knowledge, it is the first GWAS report of NPC showing that multiple loci (HLA-A, HLA-F, and GABBR1) within chromosome 6p21.3 are associated with NPC. Although some of these relationships may be attributed to linkage disequilibrium between the loci, the findings clearly provide a fresh direction for the study of NPC development.     

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.     

Identification of New Genetic Risk Variants for Type 2 Diabetes

Although more than 20 genetic susceptibility loci have been reported for type 2 diabetes (T2D), most reported variants have small to moderate effects and account for only a small proportion of the heritability of T2D, suggesting that the majority of inter-person genetic variation in this disease remains to be determined. We conducted a multistage, genome-wide association study (GWAS) within the Asian Consortium of Diabetes to search for T2D susceptibility markers. From 590,887 SNPs genotyped in 1,019 T2D cases and 1,710 controls selected from Chinese women in Shanghai, we selected the top 2,100 SNPs that were not in linkage disequilibrium (r²<0.2) with known T2D loci for in silico replication in three T2D GWAS conducted among European Americans, Koreans, and Singapore Chinese. The 5 most promising SNPs were genotyped in an independent set of 1,645 cases and 1,649 controls from Shanghai, and 4 of them were further genotyped in 1,487 cases and 3,316 controls from 2 additional Chinese studies. Consistent associations across all studies were found for rs1359790 (13q31.1), rs10906115 (10p13), and rs1436955 (15q22.2) with P-values (per allele OR, 95%CI) of 6.49×10⁻⁹ (1.15, 1.10–1.20), 1.45×10⁻⁸ (1.13, 1.08–1.18), and 7.14×10⁻⁷ (1.13, 1.08–1.19), respectively, in combined analyses of 9,794 cases and 14,615 controls. Our study provides strong evidence for a novel T2D susceptibility locus at 13q31.1 and the presence of new independent risk variants near regions (10p13 and 15q22.2) reported by previous GWAS.     

Genetic Variants Associated with Lipid Profiles in Chinese Patients with Type 2 Diabetes

Dyslipidemia is a strong risk factor for cardiovascular disease among patients with type 2 diabetes (T2D). The aim of this study was to identify lipid-related genetic variants in T2D patients of Han Chinese ancestry. Among 4,908 Chinese T2D patients who were not taking lipid-lowering medications, single nucleotide polymorphisms (SNPs) in seven genes previously found to be associated with lipid traits in genome-wide association studies conducted in populations of European ancestry (ABCA1, GCKR, BAZ1B, TOMM40, DOCK7, HNF1A, and HNF4A) were genotyped. After adjusting for multiple covariates, SNPs in ABCA1, GCKR, BAZ1B, TOMM40, and HNF1A were identified as significantly associated with triglyceride levels in T2D patients (P < 0.05). The associations between the SNPs in ABCA1 (rs3890182), GCKR (rs780094), and BAZ1B (rs2240466) remained significant even after correction for multiple testing (P = 8.85×10⁻³, 7.88×10⁻⁷, and 2.03×10⁻⁶, respectively). BAZ1B (rs2240466) also was associated with the total cholesterol level (P = 4.75×10⁻²). In addition, SNP rs157580 in TOMM40 was associated with the low-density lipoprotein cholesterol level (P = 6.94×10⁻³). Our findings confirm that lipid-related genetic loci are associated with lipid profiles in Chinese patients with type 2 diabetes.     

Genome-wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10⁻⁵ were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10⁻⁸, OR = 0.69, 95% CI [0.60–0.79]; rs6457617, P = 1.14×10⁻⁷ and rs9275245, P = 1.39×10⁻⁷. Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10⁻⁵, OR = 1.36 [1.18–1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10⁻⁵) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10⁻¹⁰, OR:1.25), TNIP1 (P = 4.68×10⁻⁹, OR:1.31), and RHOB loci (P = 3.17×10⁻⁶, OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis.     

Identification of a Functional Genetic Variant at 16q12.1 for Breast Cancer Risk: Results from the Asia Breast Cancer Consortium

Genetic factors play an important role in the etiology of breast cancer. We carried out a multi-stage genome-wide association (GWA) study in over 28,000 cases and controls recruited from 12 studies conducted in Asian and European American women to identify genetic susceptibility loci for breast cancer. After analyzing 684,457 SNPs in 2,073 cases and 2,084 controls in Chinese women, we evaluated 53 SNPs for fast-track replication in an independent set of 4,425 cases and 1,915 controls of Chinese origin. Four replicated SNPs were further investigated in an independent set of 6,173 cases and 6,340 controls from seven other studies conducted in Asian women. SNP rs4784227 was consistently associated with breast cancer risk across all studies with adjusted odds ratios (95% confidence intervals) of 1.25 (1.20−1.31) per allele (P = 3.2×10⁻²⁵) in the pooled analysis of samples from all Asian samples. This SNP was also associated with breast cancer risk among European Americans (per allele OR  = 1.19, 95% CI  = 1.09−1.31, P = 1.3×10⁻⁴, 2,797 cases and 2,662 controls). SNP rs4784227 is located at 16q12.1, a region identified previously for breast cancer risk among Europeans. The association of this SNP with breast cancer risk remained highly statistically significant in Asians after adjusting for previously-reported SNPs in this region. In vitro experiments using both luciferase reporter and electrophoretic mobility shift assays demonstrated functional significance of this SNP. These results provide strong evidence implicating rs4784227 as a functional causal variant for breast cancer in the locus 16q12.1 and demonstrate the utility of conducting genetic association studies in populations with different genetic architectures.     

Genome-Wide Association Study Identifies Major Loci for Carcass Weight on BTA14 in Hanwoo (Korean Cattle)

This genome-wide association study (GWAS) was conducted to identify major loci that are significantly associated with carcass weight, and their effects, in order to provide increased understanding of the genetic architecture of carcass weight in Hanwoo. This genome-wide association study identified one major chromosome region ranging from 23 Mb to 25 Mb on chromosome 14 as being associated with carcass weight in Hanwoo. Significant Bonferroni-corrected genome-wide associations (P<1.52×10⁻⁶) were detected for 6 Single Nucleotide Polymorphic (SNP) loci for carcass weight on chromosome 14. The most significant SNP was BTB-01280026 (P = 4.02×10⁻¹¹), located in the 25 Mb region on Bos taurus autosome 14 (BTA14). The other 5 significant SNPs were Hapmap27934-BTC-065223 (P = 4.04×10⁻¹¹) in 25.2 Mb, BTB-01143580 (P = 6.35×10⁻¹¹) in 24.3 Mb, Hapmap30932-BTC-011225 (P = 5.92×10⁻¹⁰) in 24.8 Mb, Hapmap27112-BTC-063342 (P = 5.18×10⁻⁹) in 25.4 Mb, and Hapmap24414-BTC-073009 (P = 7.38×10⁻⁸) in 25.4 Mb, all on BTA 14. One SNP (BTB-01143580; P = 6.35×10⁻¹¹) lies independently from the other 5 SNPs. The 5 SNPs that lie together showed a large Linkage disequilibrium (LD) block (block size of 553 kb) with LD coefficients ranging from 0.53 to 0.89 within the block. The most significant SNPs accounted for 6.73% to 10.55% of additive genetic variance, which is quite a large proportion of the total additive genetic variance. The most significant SNP (BTB-01280026; P = 4.02×10⁻¹¹) had 16.96 kg of allele substitution effect, and the second most significant SNP (Hapmap27934-BTC-065223; P = 4.04×10⁻¹¹) had 18.06 kg of effect on carcass weight, which correspond to 44% and 47%, respectively, of the phenotypic standard deviation for carcass weight in Hanwoo cattle. Our results demonstrated that carcass weight was affected by a major Quantitative Trait Locus (QTL) with a large effect and by many SNPs with small effects that are normally distributed.     

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.