Suggestion of GLYAT gene underlying variation of bone size and body lean mass as revealed by a bivariate genome-wide association study

Bone and muscle, two major tissue types of musculoskeletal system, have strong genetic determination. Abnormality in bone and/or muscle may cause musculoskeletal diseases such as osteoporosis and sarcopenia. Bone size phenotypes (BSPs), such as hip bone size (HBS), appendicular bone size (ABS), are genetically correlated with body lean mass (mainly muscle mass). However, the specific genes shared by these phenotypes are largely unknown. In this study, we aimed to identify the specific genes with pleiotropic effects on BSPs and appendicular lean mass (ALM). We performed a bivariate genome-wide association study (GWAS) by analyzing ~690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) followed by a replication study in 2,286 unrelated US Caucasians (558 males and 1,728 females). We identified 14 interesting single nucleotide polymorphisms (SNPs) that may contribute to variation of both BSPs and ALM, with p values <10^?6 in discovery stage. Among them, the association of three SNPs (rs2507838, rs7116722, and rs11826261) in/near GLYAT (glycine-N-acyltransferase) gene was replicated in US Caucasians, with p values ranging from 1.89 × 10^?3 to 3.71 × 10^?4 for ALM–ABS, from 5.14 × 10^?3 to 1.11 × 10^?2 for ALM–HBS, respectively. Meta-analyses yielded stronger association signals for rs2507838, rs7116722, and rs11826261, with pooled p values of 1.68 × 10^?8, 7.94 × 10^?8, 6.80 × 10^?8 for ALB–ABS and 1.22 × 10^?4, 9.85 × 10^?5, 3.96 × 10^?4 for ALM–HBS, respectively. Haplotype allele ATA based on these three SNPs was also associated with ALM–HBS and ALM–ABS in both discovery and replication samples. Interestingly, GLYAT was previously found to be essential to glucose metabolism and energy metabolism, suggesting the gene’s dual role in both bone development and muscle growth. Our findings, together with the prior biological evidence, suggest the importance of GLYAT gene in co-regulation of bone phenotypes and body lean mass.     

Genetic variants in the JAK1 gene confer higher risk of Behcet’s disease with ocular involvement in Han Chinese

Recent surveys have identified SLC22A4, SLC22A5, RUNX1, JAK1 as susceptibility genes for various immune-related diseases. An association study was performed in 738 Behcet’s patients with ocular involvement and 1,873 controls using the iPLEX system method. The first-stage study for 30 SNPs showed that SNPs rs2780815, rs310241, rs3790532 in JAK1 were associated with Behcet’s disease in Han Chinese (Pc_{(Bonferroni correction)} = 0.022–7.7 × 10^?3). The G allele and AA genotype of SNP rs2834643 in RUNX1 (Pc = 0.041–1.75 × 10^?3), but none of the other SNPs, were associated with Behcet’s disease. Haplotype analysis for the SLC22A4, SLC22A5 genes showed an increased tendency for AGTCTGCCGC frequency in patients compared with controls; however, the significance was lost after Bonferroni correction (P = 0.004, Pc > 0.05). Subsequently, we further replicated the significantly associated SNPs using another independent cohort. Replication and combining studies showed that three SNPs rs2780815, rs310241, rs3790532 in JAK1, but not SNP rs2834643 in RUNX1, were consistently associated with Behcet’s disease (replication: Pc = 0.012–9.60 × 10^?4; combining: Pc = 0.030–1.90 × 10^?4). SNPs rs2780815, rs310241, rs3790532 were estimated to confer a population attributable risk of 35.0, 28.0, 27.0 %, respectively. We found a strong association between HLA-B51 with Behcet’s disease in Chinese Han population (P = 1.35 × 10^?73; OR = 5.15; 95 % CI 4.28–6.19). GMDR analysis showed that no gene–gene interaction was detectable between JAK1 and HLA-B51. Logistic analysis indicated that the JAK1 gene was an independent risk factor for Behcet’s disease (P > 0.05). Real-time PCR analysis showed that no difference on the expression of JAK1 in PBMCs or LPS-stimulated PBMCs between individuals with the different rs1762780815 genotypes studied (P > 0.05). In conclusion, this study suggests that JAK1, but not SLC22A4, SLC22A5 and RUNX1, contributes to the genetic susceptibility to Behcet’s disease with ocular involvement.     

BET1L and TNRC6B associate with uterine fibroid risk among European Americans

Uterine fibroid (UFs) affect 77 % of women by menopause and account for $9.4 billion in healthcare costs each year. Although UFs are heritable, genetic risk is poorly understood. The first genome-wide association study (GWAS) of UFs was recently performed in a Japanese population, with reported genome-wide significance for single nucleotide polymorphisms (SNPs) across three chromosomal regions. We tested these SNPs for association with UFs in US cohorts. Women were enrolled in the Right from the Start (RFTS) cohort and the BioVU DNA repository. UF status in both cohorts was determined by pelvic imaging. We tested 65 candidate and haplotype-tagging SNPs for association with UFs presence using logistic regression in RFTS and the top three GWAS-associated SNPs in BioVU. We also combined association results from both cohorts using meta-analysis. 1,086 European American (EA) cases and 1,549 controls were examined. Two SNP associations replicated [blocked early in transport 1 homolog (BET1L) rs2280543, RFTS–BioVU meta-odds ratio (OR) = 0.67 95 % confidence interval (CI) 0.38–0.96, Q = 0.70, I = 0, p = 6.9 × 10^?3; trinucleotide repeat containing 6B (TNRC6B) rs12484776, RFTS–BioVU meta-OR = 1.21, 95 % CI 1.07–1.35, Q = 0.24, I = 28.37, p = 8.7 × 10^?3). Meta-analyses combining evidence from RFTS, BioVU, and prior GWAS showed little heterogeneity in effect sizes across studies, with meta-p values between 7.45 × 10^?8 and 3.89 × 10^?9, which were stronger than prior GWAS and supported associations observed for all previously identified loci. These data suggest common variants increase risk for UF in both EA and Japanese populations. However, further research is needed to assess the role of these genes across other racial groups.     

Identification of genetic variation that determines human trehalase activity and its association with type 2 diabetes

A prior linkage scan in Pima Indians identified a putative locus for type two diabetes (T2D) and body mass index (BMI) on chromosome 11q23-25. Association mapping across this region identified single nucleotide polymorphisms (SNPs) in the trehalase gene (TREH) that were associated with T2D. To assess the putative connection between trehalase activity and T2D, we performed a linkage study for trehalase activity in 570 Pima Indians who had measures of trehalase activity. Strong evidence of linkage of plasma trehalase activity (LOD = 7.0) was observed in the TREH locus. Four tag SNPs in TREH were genotyped in these subjects and plasma trehalase activity was highly associated with three SNPs: rs2276064, rs117619140 and rs558907 (p = 2.2 × 10^?11–1.4 × 10^?23), and the fourth SNP, rs10790256, was associated conditionally on these three (p = 2.9 × 10^?7). Together, the four tag SNPs explained 51 % of the variance in plasma trehalase activity and 79 % of the variance attributed to the linked locus. These four tag SNPs were further genotyped in 828 subjects used for association mapping of T2D, and rs558907 was associated with T2D (odds ratio (OR) 1.94, p = 0.002). To assess replication of the T2D association, all four tag SNPs were additionally genotyped in two non-overlapping samples of Native Americans. Rs558907 was reproducibly associated with T2D in 2,942 full-heritage Pima Indians (OR 1.27 p = 0.03) and 3,897 “mixed” heritage Native Americans (OR 1.21, p = 0.03), and the strongest evidence for association came from combining all samples (OR 1.27 p = 1.6 × 10^?4, n = 7,667). However, among 320 longitudinally studied subjects, measures of trehalase activity from a non-diabetic exam did not predict those who would eventually develop diabetes versus those who would remain non-diabetic (hazard ratio 0.94 per SD of trehalase activity, p = 0.29). We conclude that variants in TREH control trehalase activity, and although one of these variants is also reproducibly associated with T2D, it is likely that the effect of the SNP on risk of T2D occurs by a mechanism different than affecting trehalase activity. Alternatively, TREH variants may be tagging a nearby T2D locus.     

HLA-DPB1*04:01 allele is associated with non-obstructive azoospermia in Japanese patients

Azoospermia is defined by absence of sperm in the semen and can either be caused by obstruction of the seminal tract (obstructive azoospermia) or by defects in spermatogenesis (non-obstructive azoospermia, NOA). Previous studies reported that specific alleles and single nucleotide polymorphisms (SNPs) in the human leukocyte antigen (HLA) region were associated with NOA in East Asians. We attempt to expand upon previous findings by genotyping more HLA genes and to replicate SNP associations by focusing on Japanese NOA patients. HLA typing of six genes (HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1) was done on 355 NOA patients using SSO-Luminex assay while genotyping of two previously reported SNPs (rs498422 and rs3129878) was done on 443 patients and 544 fertile males using TaqMan assay. Association between the HLA alleles and SNP with NOA was assessed with Chi squared and logistic regression tests. We found that HLA-DPB1*04:01 [corrected p value, P _c 7.13 × 10^?6; odds ratio (OR) 2.52], DRB1*13:02 (P _c 4.93 × 10^?4, OR 1.97), DQB1*06:04 (P _c 8.94 × 10^?4, OR 1.91) and rs3129878 (p value 3.98 × 10^?4; OR 1.32) showed significant association with NOA, however, these loci are in linkage disequilibrium with each other. The conditional logistic regression tests showed that DPB1*04:01 is independently associated with NOA, confirming the involvement of the HLA region in the etiology of NOA in Japanese patients.     

A single nucleotide polymorphism in the FADS1/FADS2 gene is associated with plasma lipid profiles in two genetically similar Asian ethnic groups with distinctive differences in lifestyle

Recent genome-wide association studies (GWASs) showed that single nucleotide polymorphisms (SNPs) in FADS1/FADS2 were associated with plasma lipid concentrations in populations with European ancestry. We investigated the associations between the SNPs in FADS1/FADS2 and plasma concentrations of triglycerides, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) in two Asian groups, i.e., Japanese and Mongolians. The genotype of rs174547 (T/C), found to be associated with triglyceride and HDL-C concentrations in the GWAS, was determined in 21,004 Japanese and 1,203 Mongolian individuals. Genotype–phenotype association was assessed by using multiple linear regression models, assuming an additive model of inheritance. The copy number of the rs174547 C allele was significantly associated with increased triglyceride levels (P = 1.5 × 10^?6) and decreased HDL-C levels (P = 0.03) in the Japanese population. On the other hand, in the Mongolian population, the rs174547 C allele copy number was strongly associated with decreased LDL-C levels (P = 2.6 × 10^?6), but was not associated with triglyceride and HDL-C levels. The linkage disequilibrium pattern and haplotype structures of SNPs around the FADS1/FADS2 locus showed no marked dissimilarity between Japanese and Mongolian individuals. The present data indicate that the FADS1/FADS2 locus can be added to the growing list of loci involved in polygenic dyslipidemia in Asians. Furthermore, the variable effects of FADS1/FADS2 on plasma lipid profiles in Asians may result from differences in the dietary intake of polyunsaturated fatty acids, which serve as substrates for enzymes encoded by FADS1/FADS2.     

A meta-analysis of two genome-wide association studies to identify novel loci for maximum number of alcoholic drinks

Maximum number of alcoholic drinks consumed in a 24-h period (maxdrinks) is a heritable (>50 %) trait and is strongly correlated with vulnerability to excessive alcohol consumption and subsequent alcohol dependence (AD). Several genome-wide association studies (GWAS) have studied alcohol dependence, but few have concentrated on excessive alcohol consumption. We performed two GWAS using maxdrinks as an excessive alcohol consumption phenotype: one in 118 extended families (N = 2,322) selected from the Collaborative Study on the Genetics of Alcoholism (COGA), and the other in a case–control sample (N = 2,593) derived from the Study of Addiction: Genes and Environment (SAGE). The strongest association in the COGA families was detected with rs9523562 (p = 2.1 × 10^?6) located in an intergenic region on chromosome 13q31.1; the strongest association in the SAGE dataset was with rs67666182 (p = 7.1 × 10^?7), located in an intergenic region on chromosome 8. We also performed a meta-analysis with these two GWAS and demonstrated evidence of association in both datasets for the LMO1 (p = 7.2 × 10^?7) and PLCL1 genes (p = 4.1 × 10^?6) with maxdrinks. A variant in AUTS2 and variants in INADL, C15orf32 and HIP1 that were associated with measures of alcohol consumption in a meta-analysis of GWAS studies and a GWAS of alcohol consumption factor score also showed nominal association in the current meta-analysis. The present study has identified several loci that warrant further examination in independent samples. Among the top SNPs in each of the dataset (p ≤ 10^?4) far more showed the same direction of effect in the other dataset than would be expected by chance (p = 2 × 10^?3, 3 × 10^?6), suggesting that there are true signals among these top SNPs, even though no SNP reached genome-wide levels of significance.     

Large multiethnic Candidate Gene Study for C-reactive protein levels: identification of a novel association at CD36 in African Americans

C-reactive protein (CRP) is a heritable biomarker of systemic inflammation and a predictor of cardiovascular disease (CVD). Large-scale genetic association studies for CRP have largely focused on individuals of European descent. We sought to uncover novel genetic variants for CRP in a multiethnic sample using the ITMAT Broad-CARe (IBC) array, a custom 50,000 SNP gene-centric array having dense coverage of over 2,000 candidate CVD genes. We performed analyses on 7,570 African Americans (AA) from the Candidate gene Association Resource (CARe) study and race-combined meta-analyses that included 29,939 additional individuals of European descent from CARe, the Women’s Health Initiative (WHI) and KORA studies. We observed array-wide significance (p < 2.2 × 10^?6) for four loci in AA, three of which have been reported previously in individuals of European descent (IL6R, p = 2.0 × 10^?6; CRP, p = 4.2 × 10^?71; APOE, p = 1.6 × 10^?6). The fourth significant locus, CD36 (p = 1.6 × 10^?6), was observed at a functional variant (rs3211938) that is extremely rare in individuals of European descent. We replicated the CD36 finding (p = 1.8 × 10^?5) in an independent sample of 8,041 AA women from WHI; a meta-analysis combining the CARe and WHI AA results at rs3211938 reached genome-wide significance (p = 1.5 × 10^?10). In the race-combined meta-analyses, 13 loci reached significance, including ten (CRP, TOMM40/APOE/APOC1, HNF1A, LEPR, GCKR, IL6R, IL1RN, NLRP3, HNF4A and BAZ1B/BCL7B) previously associated with CRP, and one (ARNTL) previously reported to be nominally associated with CRP. Two novel loci were also detected (RPS6KB1, p = 2.0 × 10^?6; CD36, p = 1.4 × 10^?6). These results highlight both shared and unique genetic risk factors for CRP in AA compared to populations of European descent.     

Replication of genomewide association studies on age at menarche in the Korean population

Early menarche is associated with adverse health outcomes, including breast cancer, endometrial cancer, obesity, type 2 diabetes, and cardiovascular disease. Recently, a genomewide association study (GWAS) of age at menarche (AAM) in 104,533 individuals of European ancestry was reported by the ReproGen consortium. They identified 42 loci known and novel loci that were linked to age at menarche. Because age at menarche varies between ethnic groups, we decided to investigate if these results would be replicated in the Korean population. To this end, we examined the association of the SNPs reported in the ReproGen GWAS with AAM in 3,194 individuals from the Korean Genome and Epidemiology Study (KoGES) cohort. Genotype data for total 17 SNPs (6 genotyped SNPs and 11 imputed SNPs) were available for the association analysis using linear regression analysis for age at menarche with controlling current age, waist-to-hip ratio, and body mass index as the covariates. We found replication of the ReproGen study in two SNPs; one SNP (rs466639) in the retinoic acid receptor gamma gene (RXRG), showing a significant association with early menarche (beta = ?0.224 ± 0.065, p value = 5.2 × 10^?4, Bonferroni-corrected p value = 0.009), and the other (rs10899489), in GRB2 (growth factor receptor bound protein 2)-associated binding protein 2 (GAB2), linked to late menarche (beta = 0.140 ± 0.047, p value = 2.8 × 10^?3, Bonferroni-corrected p value = 0.049). This result possibly suggests that genetic factors governing AAM in the Korean population would be distinct from those in the Europeans, implying roles of modulating or interacting factors in determining AAM, including environmental factors such as nutritional status.     

Genome-wide association study of aspirin-exacerbated respiratory disease in a Korean population

Aspirin-exacerbated respiratory disease (AERD) is a nonallergic clinical syndrome characterized by a severe decline in forced expiratory volume in one second (FEV1) following the ingestion of non-steroidal anti-inflammatory drugs (NSAIDs) such as aspirin. The effects of genetic variants have not fully explained all of the observed individual differences to an aspirin challenge despite previous attempts to identify AERD-related genes. In the present study, we performed genome-wide association study (GWAS) and targeted association study in Korean asthmatics to identify new genetic factors associated with AERD. A total of 685 asthmatic patients without AERD and 117 subjects with AERD were used for the GWAS of the first stage, and 996 asthmatics without AERD and 142 subjects with AERD were used for a follow-up study. A total of 702 SNPs were genotyped using the GoldenGate assay with the VeraCode microbead. GWAS revealed the top-ranked variants in 3′ regions of the HLA-DPB1 gene. To investigate the detailed genetic effects of an associated region with the risk of AERD, a follow-up targeted association study with the 702 single nucleotide polymorphisms (SNPs) of 14 genes was performed on 802 Korean subjects. In a case–control analysis, HLA-DPB1 rs1042151 (Met105Val) shows the most significant association with the susceptibility of AERD (p = 5.11 × 10^?7; OR = 2.40). Moreover, rs1042151 also shows a gene dose for the percent decline of FEV1 after an aspirin challenge (p = 2.82 × 10^?7). Our findings show that the HLA-DPB1 gene polymorphism may be the most susceptible genetic factor for the risk of AERD in Korean asthmatics and confirm the importance of HLA-DPB1 in the genetic etiology of AERD.     

Genetic risk factors for myocardial infarction more clearly manifest for early age of first onset

Epidemiological genetics established that heritability in determining the risk of myocardial infarction (MI) is substantially greater when MI occurs early in life. However, the genetic architecture of early-onset and late-onset MI was not compared. We analyzed genotype frequencies of SNPs in/near 20 genes whose protein products are involved in the pathogenesis of atherosclerosis in two groups of Russian patients with MI: the first group included patients with age of first MI onset <60 years (N?=?230) and the second group with onset ≥60 years (N?=?174). The control group of corresponding ethnicity consisted of 193 unrelated volunteers without cardiovascular diseases (93 individuals were over 60 years). We found that in the group of patients with age of onset <60 years, SNPs FGB rs1800788*T, TGFB1 rs1982073*T/T, ENOS rs2070744*C and CRP rs1130864*T/T were associated with risk of MI, whereas in patients with age of onset ≥60 years, only TGFB1 rs1982073*T/T was associated with risk of MI. Using APSampler software, we found composite markers associated with MI only in patients with early onset: FGB rs1800788*T?+?TGFB1 rs1982073*T; FGB rs1800788*T?+?LPL rs328*C?+?IL4 rs2243250*C; FGB rs1800788*T?+?ENOS rs2070744*C (Fisher p values of 1.4?×?10^?6 to 2.2?×?10^?5; the permutation p values of 1.1?×?10^?5 to 3.0?×?10^?4; ORs?=?2.67–2.54). Alleles included in the combinations were associated with MI less significantly and with lower ORs than the combinations themselves. The result showed a substantially greater contribution of the genetic component in the development of MI if it occurs early in life, and demonstrated the usefulness of genetic testing for young people.     

Identification of BACH2 as a susceptibility gene for Graves’ disease in the Chinese Han population based on a three-stage genome-wide association study

The BACH2 gene regulates B cell differentiation and function and has been reported to be a shared susceptibility gene for several autoimmune diseases. Our previous genome-wide association study (GWAS) indicated that several single nucleotide polymorphisms (SNPs) in the BACH2 gene are associated with Graves’ disease (GD) in the Chinese Han population; however, the association did not achieve genome-wide significance levels. Recently, this association of BACH2 with GD was confirmed in Caucasians in the UK population, but fine mapping in this region has not yet been reported. Here, we provide a refined analysis of a 331-kb region in the BACH2 gene, which harbors 359 SNPs, using GWAS data from 1,442 GD patients and 1,468 controls. The SNPs rs2474619 and rs9344996 were implied as the independent variants associated with GD by forward and two-locus logistic regression analysis. We genotyped eight out of 10 tagSNPs with P < 1 × 10^?3 in 3,508 GD patients and 3,209 controls, the results also showed that rs2474619 was independently associated with GD in the combined population from GWAS and the second stage (P = 1.81 × 10^?5). The rs2474619 and rs9344996 were further genotyped in the third stage cohorts, and rs2474619 showed evidence of association with GD at genome-wide significance levels in the combined population (P = 3.28 × 10^?8, odds ratio = 1.13). The association of rs9344996 with GD can be explained by its linkage to rs2474619 in the combined population. Our study clearly demonstrated that BACH2 is a susceptibility gene for GD in the Chinese Han population and further supported rs2474619, in intron 2 of BACH2, is the best association signal with GD. However, the mechanism by which BACH2 confers increased risk of GD requires further study.     

Combined effect of five single nucleotide polymorphisms related to IL23/Th17 pathway in the risk of psoriasis

Psoriasis is a common immune-mediated inflammatory skin disease with strong genetic components, in which the IL23/Th17 pathway has been implicated. To explore the effective role in psoriasis, we genotyped five single nucleotide polymorphisms in genes related to IL23/Th17 pathway in 14,929 Han Chinese samples. A Bonferroni-corrected significant single-variant association was identified between rs1512970 within IL21 (odds ratio (OR)?=?1.07, 95 % confidence interval (CI)?=?1.02–1.13, P?=?4.94?×?10^?03). We failed to validate rs744166 (OR?=?1.06, 95 % CI?=?1.01–1.11, P?=?1.52?×?10^?02) and other three SNPs (P?=?2.48?×?10^?01?～?1.27?×?10^?02) to meet the single-variant association significance threshold. However, we found that their combined effect substantially contributed to the risk of psoriasis in our sample (P?=?3.91?×?10^?07) and the highest score group conferred the largest risk effect size (OR?=?1.22, 95 % CI?=?1.11–1.34, P?=?1.85?×?10^?05). Our results implicate the ethnic heterogeneity in the susceptibility of psoriasis and suggest common variants with weak effect in IL23/Th17 pathway, which do not show significance in conventional single-variant association study, may contribute to the risk of psoriasis. This study sheds light on the important role of IL23/Th17 pathway in the susceptibility of psoriasis.     

TNFAIP3 gene polymorphisms confer risk for Behcet’s disease in a Chinese Han population

The tumor necrosis factor alpha-inducible protein 3 (TNFAIP3) gene polymorphisms have recently been reported to be associated with the susceptibility to several immune-related diseases. This study was performed to evaluate the potential association of TNFAIP3 polymorphisms with Behcet’s disease (BD) in a Chinese Han population. Five single-nucleotide polymorphisms (SNPs), rs10499194, rs610604, rs7753873, rs5029928, and rs9494885 of TNFAIP3 were genotyped in 722 BD patients and 1,415 healthy controls using a PCR-restriction fragment length polymorphism assay. Allele and genotype frequencies were compared between patients and controls using the χ ² test. The results showed a significantly increased prevalence of the rs9494885 TC genotype and C allele in BD patients compared with controls (Bonferroni corrected p (p _c) = 1.83 × 10^?10, odds ratio (OR) [95 % CI] 2.03 [1.65–2.49]; p _c = 8.35 × 10^?10, OR [95 % CI] 1.81 [1.51–2.18], respectively).The frequency of the TT genotype and T allele of rs9494885 was markedly lower in BD patients than that in controls (p _c = 1.23 × 10^?10, OR [95 % CI] 0.50 [0.40–0.61]; p _c = 8.35 × 10^?10, OR [95 % CI] 0.55 [0.46–0.66], respectively). For rs10499194, a higher frequency of the CC genotype (p _c = 0.015, OR [95 % CI] 1.96 [1.30–2.97]) and C allele (p _c = 0.005, OR [95 % CI] 1.92 [1.28–2.90]), and a lower frequency of the TC genotype (p _c = 0.015, OR [95 % CI] 0.51 [0.34–0.77]) and T allele (p _c = 0.005, OR [95 % CI] 0.52 [0.35–2.97]) were found in BD patients. Concerning rs7753873, a higher frequency of the AC genotype (p _c = 0.015, OR [95 % CI] 1.49 [1.17–1.91]) and C allele (p _c = 0.025, OR [95 % CI] 1.39 [1.11–1.76]), and a lower frequency of the AA genotype (p _c = 0.03, OR [95 % CI] 0.68 [0.53–0.87]) and A allele (p _c = 0.025, OR [95 % CI] 0.72 [0.57–0.91]) were observed in BD patients. This study identified one strong risk SNP rs9494885 and two weak risk SNPs rs10499194 and rs7753873 of TNFAIP3 in Chinese Han BD patients.     

Genome-wide association study in Han Chinese identifies three novel loci for human height

Human height is a complex genetic trait with high heritability but discovery efforts in Asian populations are limited. We carried out a meta-analysis of genome-wide association studies (GWAS) for height in 6,534 subjects with in silico replication of 1,881 subjects in Han Chinese. We identified three novel loci reaching the genome-wide significance threshold (P < 5 × 10^?8), which mapped in or near ZNF638 (rs12612930, P = 2.02 × 10^?10), MAML2 (rs11021504, P = 7.81 × 10^?9), and C18orf12 (rs11082671, P = 1.87 × 10^?8). We also confirmed two loci previously reported in European populations including CS (rs3816804, P = 2.63 × 10^?9) and CYP19A1 (rs3751599, P = 4.80 × 10^?10). In addition, we provided evidence supporting 35 SNPs identified by previous GWAS (P < 0.05). Our study provides new insights into the genetic determination of biological regulation of human height.     

CUBN and NEBL common variants in the chromosome 10p13 linkage region are associated with multibacillary leprosy in Vietnam

Leprosy is caused by infection with Mycobacterium leprae and is classified clinically into paucibacillary (PB) or multibacillary (MB) subtypes based on the number of skin lesions and the bacillary index detected in skin smears. We previously identified a major PB susceptibility locus on chromosome region 10p13 in Vietnamese families by linkage analysis. In the current study, we conducted high-density association mapping of the 9.5 Mb linkage peak on chromosome region 10p13 covering 39 genes. Using leprosy per se and leprosy subtypes as phenotypes, we employed 294 nuclear families (303 leprosy cases, 63 % MB, 37 % PB) as a discovery sample and 192 nuclear families (192 cases, 55 % MB, 45 % PB) as a replication sample. Replicated significant association signals were revealed in the genes for cubilin (CUBN) and nebulette (NEBL). In the combined sample, the C allele (frequency 0.26) at CUBN SNP rs10904831 showed association [p = 1 × 10^?5; OR 0.52 (0.38–0.7)] with MB leprosy only. Likewise, allele T (frequency 0.42) at NEBL SNP rs11012461 showed association [p = 4.2 × 10^?5; OR 2.51 (1.6–4)] with MB leprosy only. These associations remained valid for the CUBN signal when taking into account the effective number of tests performed (type I error significance threshold = 2.4 × 10^?5). We used the results of our analyses to propose a new model for the genetic control of polarization of clinical leprosy.     

Replication of an African-American GWAS on blood pressure and hypertension in the Korean population

Hypertension is a complex disease that is caused by the interaction of multiple genetic and environmental risk factors, affecting 30% adult in industrialized countries. The identification of genetic factors that impact one’s predisposition to hypertension and its progression is an ongoing challenge. A genome-wide association study of African-Americans, who have one of the highest rates of hypertension in the world, was reported. We replicated the GWAS results in 8842 unrelated Koreans. Fifteen of the 30 reported SNPs were analyzed for their association with blood pressure and hypertension. Linear regression was used to analyze blood pressure as a quantitative trait in 7551 subjects, and a case-control study was performed using 1968 hypertensive cases and 4452 normotensive controls by logistic regression analysis. The quantitative trait study demonstrated a moderate association of 2 SNPs, rs9301196 (p=4.9×10^?3 for diastolic blood pressure) and rs2823756 (p=0.04 for systolic blood pressure), which were also associated with hypertension (p=0.042 and p=6.3×10^?3, respectively). Further, 3 SNPs were associated with hypertension (p=0.042 for rs7902529, p=0.027 for rs10135446, and p=0.01 for rs4613079) but not with blood pressure. Based on the moderate association signals and the low proportion of positive signals, this cross validation between African-Americans and Asians suggests that association studies of blood pressure traits require a larger number of subjects and a more refined design.     

Variants in BET1L and TNRC6B associate with increasing fibroid volume and fibroid type among European Americans

Uterine fibroids (UFs) affect 77 % of women by menopause and account for $9.4 billion in yearly healthcare costs. We recently replicated findings from the first UF genome-wide association study (GWAS), conducted in the Japanese. Here we tested these GWAS-discovered SNPs for association with UF characteristics to further assess whether risk varies by sub-phenotypes of UFs. Women were enrolled in Right from the Start (RFTS) and the BioVU DNA repository (BioVU). UF status was determined by pelvic imaging. We tested the top GWAS-associated SNPs for association with UF characteristics (RFTS: type, number, volume; BioVU: type) using covariate adjusted logistic and linear regression. We also combined association results of UF type using meta-analysis. 456 European American (EA) cases and 1,549 controls were examined. Trinucleotide repeat containing 6B (TNRC6B) rs12484776 associated with volume in RFTS (β = 0.40, 95 % CI 0.05–0.75, p = 0.024). RFTS analyses evaluating stratified quartiles of volume showed the strongest OR at rs12484776 for the largest volume (16.6–179.1 cc, odds ratio (OR) = 2.19, 95 % confidence interval (CI) 1.07–4.46, p = 0.031). Meta-analysis showed a strong association at blocked early in transport 1 homolog (BET1L) rs2280543 for intramural UFs (meta-OR = 0.51, standard error (SE) = 0.14, Q = 0.590, I = 0, p = 2.48 × 10^?6), which is stronger than the overall association with UF risk. This study is the first to evaluate these SNPs for association with UF characteristics and suggests these genes associate with increasing UF volume and protection from intramural UF in EAs.     

Linkage disequilibrium pattern and age-at-diagnosis are critical for replicating genetic associations across ethnic groups in leprosy

One of the persistent challenges of genetic association studies is the replication of genetic marker-disease associations across ethnic groups. Here, we conducted high-density association mapping of PARK2/PACRG SNPs with leprosy and identified 69 SNPs significantly associated with leprosy in 198 single-case Vietnamese leprosy families. A total of 56 associated SNPs localized to the overlapping promoter regions of PARK2/PACRG. For this region, multivariate analysis identified four SNPs belonging to two major SNP bins (rs1333955, rs7744433) and two single SNP bins (rs2023004, rs6936895) that capture the combined statistical evidence (P = 1.1 × 10^?5) for association among Vietnamese patients. Next, we enrolled a case–control sample of 364 leprosy cases and 370 controls from Northern India. We genotyped all subjects for 149 SNPs that capture >80 % of the genetic variation in the Vietnamese sample and found 24 SNPs significantly associated with leprosy. Multivariate analysis identified three SNPs (rs1333955, rs9356058 and rs2023004) that capture the association with leprosy (P < 10^?8). Hence, two SNPs (rs1333955 and rs2023004) were replicated by multivariate analysis between both ethnic groups. Marked differences in the linkage disequilibrium pattern explained some of the differences in univariate analysis between the two ethnic groups. In addition, the strength of association for two promoter region SNP bins was significantly stronger among young leprosy patients in the Vietnamese sample. The same trend was observed in the Indian sample, but due to the higher age-at-diagnosis of the patients the age effect was less pronounced.     

Identification of six novel susceptibility loci for dyslipidemia using longitudinal exome-wide association studies in a Japanese population

Recent genome-wide association studies have identified various dyslipidemia-related genetic variants. However, most studies were conducted in a cross-sectional manner. We thus performed longitudinal exome-wide association studies of dyslipidemia in a Japanese population. We used ~244,000 genetic variants and clinical data of 6022 Japanese individuals who had undergone annual health checkups for several years. After quality control, the association of dyslipidemia-related phenotypes with 24,691 single nucleotide polymorphisms (SNPs) was tested using the generalized estimating equation model. In total, 82 SNPs were significantly (P < 2.03 × 10^?6) associated with dyslipidemia phenotypes. Of these SNPs, four (rs74416240 of TCHP, rs925368 of GIT2, rs7969300 of ATXN2, and rs12231744 of NAA25) and two (rs34902660 of SLC17A3 and rs1042127 of CDSN) were identified as novel genetic determinants of hypo-HDL- and hyper-LDL-cholesterolemia, respectively. A replication study using the cross-sectional data of 8310 Japanese individuals showed the association of the six identified SNPs with dyslipidemia-related traits.