WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ～2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)相似文献   

Association of chromosome 8q24 variants with prostate cancer risk in the Siberian region of Russia and meta-analysis

Compelling evidence demonstrates chromosome 8q24 as a prostate cancer susceptibility locus. In present work we studied whether the common variants of 8q24 region, rs6983267 and rs1447295, were associated with the sporadic prostate cancer risk in the Russian population. Polymorphisms were genotyped in 393 case and 384 control Russian Caucasian men from Siberia region. The A allele of rs1447295 was significantly associated with the risk of prostate cancer (OR[CI 95%] = 1.74 [1.26-2.4], p = 7.8 x 10(-4)). A common G-A haplotype for rs6983267 - rs1447295 also showed an association with prostate cancer risk in Russian population (OR[CI 95%] = 2.03 [1.1 - 3.75], p = 0.02). We performed a meta-analysis combining our results with previous studies to evaluate the association between studied SNPs and prostate cancer risk. Meta-analysis has strongly supported the association for these SNPs (p < 10(-6)). Accordingly our study confirms the association between chromosome 8q24 and prostate cancer risk.     

Association between Expression Quantitative Trait Loci and Metabolic Traits in Two Korean Populations

Most genome-wide association studies consider genes that are located closest to single nucleotide polymorphisms (SNPs) that are highly significant for those studies. However, the significance of the associations between SNPs and candidate genes has not been fully determined. An alternative approach that used SNPs in expression quantitative trait loci (eQTL) was reported previously for Crohn’s disease; it was shown that eQTL-based preselection for follow-up studies was a useful approach for identifying risk loci from the results of moderately sized GWAS. In this study, we propose an approach that uses eQTL SNPs to support the functional relationships between an SNP and a candidate gene in a genome-wide association study. The genome-wide SNP genotypes and 10 biochemical measures (fasting glucose levels, BUN, serum albumin levels, AST, ALT, gamma GTP, total cholesterol, HDL cholesterol, triglycerides, and LDL cholesterol) were obtained from the Korean Association Resource (KARE) consortium. The eQTL SNPs were isolated from the SNP dataset based on the RegulomeDB eQTL-SNP data from the ENCODE projects and two recent eQTL reports. A total of 25,658 eQTL SNPs were tested for their association with the 10 metabolic traits in 2 Korean populations (Ansung and Ansan). The proportion of phenotypic variance explained by eQTL and non-eQTL SNPs showed that eQTL SNPs were more likely to be associated with the metabolic traits genetically compared with non-eQTL SNPs. Finally, via a meta-analysis of the two Korean populations, we identified 14 eQTL SNPs that were significantly associated with metabolic traits. These results suggest that our approach can be expanded to other genome-wide association studies.     

Phosphodiesterase 8B gene variants are associated with serum TSH levels and thyroid function

Thyroid-stimulating hormone (TSH) controls thyroid growth and hormone secretion through binding to its G protein-coupled receptor (TSHR) and production of cyclic AMP (cAMP). Serum TSH is a sensitive indicator of thyroid function, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over a life span. By genotyping 362,129 SNPs in 4,300 Sardinians, we identified a strong association (p = 1.3 x 10(-11)) between alleles of rs4704397 and circulating TSH levels; each additional copy of the minor A allele was associated with an increase of 0.13 muIU/ml in TSH. The single-nucleotide polymorphism (SNP) is located in intron 1 of PDE8B, encoding a high-affinity cAMP-specific phosphodiesterase. The association was replicated in 4,158 individuals, including additional Sardinians and two genetically distant cohorts from Tuscany and the Old Order Amish (overall p value = 1.9 x 10(-20)). In addition to association of TSH levels with SNPs in PDE8B, our genome scan provided evidence for association with PDE10A and several biologically interesting candidates in a focused analysis of 24 genes. In particular, we found evidence for association of TSH levels with SNPs in the THRB (rs1505287, p = 7.3 x 10(-5)), GNAQ (rs10512065, p = 2.0 x 10(-4)), TG (rs2252696, p = 2.2 x 10(-3)), POU1F1 (rs1976324, p = 3.9 x 10(-3)), PDE4D (rs27178, p = 8.3 x 10(-3)), and TSHR (rs4903957, p = 8.6 x 10(-3)) loci. Overall, the results suggest a primary effect of PDE8B variants on cAMP levels in the thyroid. This would affect production of T4 and T3 and feedback to alter TSH release by the pituitary. PDE8B may thus provide a candidate target for the treatment of thyroid dysfunction.     

Eight common genetic variants associated with serum DHEAS levels suggest a key role in ageing mechanisms

Dehydroepiandrosterone sulphate (DHEAS) is the most abundant circulating steroid secreted by adrenal glands--yet its function is unknown. Its serum concentration declines significantly with increasing age, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. We conducted a meta-analysis of genome-wide association data with 14,846 individuals and identified eight independent common SNPs associated with serum DHEAS concentrations. Genes at or near the identified loci include ZKSCAN5 (rs11761528; p = 3.15 × 10(-36)), SULT2A1 (rs2637125; p = 2.61 × 10(-19)), ARPC1A (rs740160; p = 1.56 × 10(-16)), TRIM4 (rs17277546; p = 4.50 × 10(-11)), BMF (rs7181230; p = 5.44 × 10(-11)), HHEX (rs2497306; p = 4.64 × 10(-9)), BCL2L11 (rs6738028; p = 1.72 × 10(-8)), and CYP2C9 (rs2185570; p = 2.29 × 10(-8)). These genes are associated with type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins. Several SNPs were associated with changes in gene expression levels, and the related genes are connected to biological pathways linking DHEAS with ageing. This study provides much needed insight into the function of DHEAS.     

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each P ≤ 4 x 10(-6)). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls (P ≤ 0.001) and were highly significant in the combined dataset (P ≤ 6 x 10(-8)). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set P = 9 x 10(-9), replication set P = 7 x 10(-4), combined P = 2 x 10(-10)). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.     

Genome-wide scan on total serum IgE levels identifies FCER1A as novel susceptibility locus

High levels of serum IgE are considered markers of parasite and helminth exposure. In addition, they are associated with allergic disorders, play a key role in anti-tumoral defence, and are crucial mediators of autoimmune diseases. Total IgE is a strongly heritable trait. In a genome-wide association study (GWAS), we tested 353,569 SNPs for association with serum IgE levels in 1,530 individuals from the population-based KORA S3/F3 study. Replication was performed in four independent population-based study samples (total n = 9,769 individuals). Functional variants in the gene encoding the alpha chain of the high affinity receptor for IgE (FCER1A) on chromosome 1q23 (rs2251746 and rs2427837) were strongly associated with total IgE levels in all cohorts with P values of 1.85 x 10(-20) and 7.08 x 10(-19) in a combined analysis, and in a post-hoc analysis showed additional associations with allergic sensitization (P = 7.78 x 10(-4) and P = 1.95 x 10(-3)). The "top" SNP significantly influenced the cell surface expression of FCER1A on basophils, and genome-wide expression profiles indicated an interesting novel regulatory mechanism of FCER1A expression via GATA-2. Polymorphisms within the RAD50 gene on chromosome 5q31 were consistently associated with IgE levels (P values 6.28 x 10(-7)-4.46 x 10(-8)) and increased the risk for atopic eczema and asthma. Furthermore, STAT6 was confirmed as susceptibility locus modulating IgE levels. In this first GWAS on total IgE FCER1A was identified and replicated as new susceptibility locus at which common genetic variation influences serum IgE levels. In addition, variants within the RAD50 gene might represent additional factors within cytokine gene cluster on chromosome 5q31, emphasizing the need for further investigations in this intriguing region. Our data furthermore confirm association of STAT6 variation with serum IgE levels.     

A large-scale genetic association study confirms IL12B and leads to the identification of IL23R as psoriasis-risk genes

We performed a multitiered, case-control association study of psoriasis in three independent sample sets of white North American individuals (1,446 cases and 1,432 controls) with 25,215 genecentric single-nucleotide polymorphisms (SNPs) and found a highly significant association with an IL12B 3'-untranslated-region SNP (rs3212227), confirming the results of a small Japanese study. This SNP was significant in all three sample sets (odds ratio [OR](common) 0.64, combined P [Pcomb]=7.85x10(-10)). A Monte Carlo simulation to address multiple testing suggests that this association is not a type I error. The coding regions of IL12B were resequenced in 96 individuals with psoriasis, and 30 additional IL12B-region SNPs were genotyped. Haplotypes were estimated, and genotype-conditioned analyses identified a second risk allele (rs6887695) located approximately 60 kb upstream of the IL12B coding region that exhibited association with psoriasis after adjustment for rs3212227. Together, these two SNPs mark a common IL12B risk haplotype (OR(common) 1.40, Pcomb=8.11x10(-9)) and a less frequent protective haplotype (OR(common) 0.58, Pcomb=5.65x10(-12)), which were statistically significant in all three studies. Since IL12B encodes the common IL-12p40 subunit of IL-12 and IL-23, we individually genotyped 17 SNPs in the genes encoding the other chains of these cytokines (IL12A and IL23A) and their receptors (IL12RB1, IL12RB2, and IL23R). Haplotype analyses identified two IL23R missense SNPs that together mark a common psoriasis-associated haplotype in all three studies (OR(common) 1.44, Pcomb=3.13x10(-6)). Individuals homozygous for both the IL12B and the IL23R predisposing haplotypes have an increased risk of disease (OR(common) 1.66, Pcomb=1.33x10(-8)). These data, and the previous observation that administration of an antibody specific for the IL-12p40 subunit to patients with psoriasis is highly efficacious, suggest that these genes play a fundamental role in psoriasis pathogenesis.     

Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci

To identify genetic factors contributing to type 2 diabetes (T2D), we performed large-scale meta-analyses by using a custom ～50,000 SNP genotyping array (the ITMAT-Broad-CARe array) with ～2000 candidate genes in 39 multiethnic population-based studies, case-control studies, and clinical trials totaling 17,418 cases and 70,298 controls. First, meta-analysis of 25 studies comprising 14,073 cases and 57,489 controls of European descent confirmed eight established T2D loci at genome-wide significance. In silico follow-up analysis of putative association signals found in independent genome-wide association studies (including 8,130 cases and 38,987 controls) performed by the DIAGRAM consortium identified a T2D locus at genome-wide significance (GATAD2A/CILP2/PBX4; p = 5.7 × 10(-9)) and two loci exceeding study-wide significance (SREBF1, and TH/INS; p < 2.4 × 10(-6)). Second, meta-analyses of 1,986 cases and 7,695 controls from eight African-American studies identified study-wide-significant (p = 2.4 × 10(-7)) variants in HMGA2 and replicated variants in TCF7L2 (p = 5.1 × 10(-15)). Third, conditional analysis revealed multiple known and novel independent signals within five T2D-associated genes in samples of European ancestry and within HMGA2 in African-American samples. Fourth, a multiethnic meta-analysis of all 39 studies identified T2D-associated variants in BCL2 (p = 2.1 × 10(-8)). Finally, a composite genetic score of SNPs from new and established T2D signals was significantly associated with increased risk of diabetes in African-American, Hispanic, and Asian populations. In summary, large-scale meta-analysis involving a dense gene-centric approach has uncovered additional loci and variants that contribute to T2D risk and suggests substantial overlap of T2D association signals across multiple ethnic groups.     

Genome-Wide Association Study of Serum Creatinine Levels during Vancomycin Therapy

Vancomycin, a commonly used antibiotic, can be nephrotoxic. Known risk factors such as age, creatinine clearance, vancomycin dose / dosing interval, and concurrent nephrotoxic medications fail to accurately predict nephrotoxicity. To identify potential genomic risk factors, we performed a genome-wide association study (GWAS) of serum creatinine levels while on vancomycin in 489 European American individuals and validated findings in three independent cohorts totaling 439 European American individuals. In primary analyses, the chromosome 6q22.31 locus was associated with increased serum creatinine levels while on vancomycin therapy (most significant variant rs2789047, risk allele A, β = -0.06, p = 1.1 x 10^-7). SNPs in this region had consistent directions of effect in the validation cohorts, with a meta-p of 1.1 x 10^-7. Variation in this region on chromosome 6, which includes the genes TBC1D32/C6orf170 and GJA1 (encoding connexin43), may modulate risk of vancomycin-induced kidney injury.     

Genome-wide association study of bronchial asthma in the Volga-Ural region of Russia

Bronchial asthma is a chronic inflammatory respiratory disease that is caused by the complex interaction of environmental influences and genetic susceptibility. The first genome-wide association study of bronchial asthma discovered a significant association between SNPs within 17q12-21 genomic region and childhood bronchial asthma in individuals of European descent. Association with this genomic region was then replicated in a number of independent samples of European and Asian descent. Here we report results of the first genome-wide association study of bronchial asthma in the Volga-Ural region of Russia. The present study includes 358 unrelated patients with physician-diagnosed bronchial asthma and 369 disease-free control subjects of different ethnic origin (Russians, Tatars and Bashkirs). Genotyping of DNA samples was carried out using the Illumina Human610 quad array as a part of GABRIEL project (contract from the EC No LSHB-CT-2006-018996). After QC filtering procedures, a final set of 550915 SNPs genotyped in 330 cases and 348 controls was tested for association with bronchial asthma. Five markers on chromosome 17q12-21 showed statistically significant association with bronchial asthma (p < or = 4.79 x 10(-7)). SNP rs7216389 with the strongest evidence for association (p = 1.01 x 10(-7)) is located within the first intron of the GSDMB gene. Evidence for association was stronger with childhood-onset asthma (p = 1.97 x 10(-6) for SNP rs7216389) compared to late-onset asthma (p = 1.8 x 10(-4) for SNP rs7216389). Our replication study using three SNPs within GSDMB gene confirmed association with only childhood-onset asthma. In summary, these results suggest an important role for genetic variants within 17q12-q21 region in the development of bronchial asthma in the Volga-Ural region of Russia.     

A replication study of GWAS-derived lipid genes in Asian Indians: the chromosomal region 11q23.3 harbors loci contributing to triglycerides

Recent genome-wide association scans (GWAS) and meta-analysis studies on European populations have identified many genes previously implicated in lipid regulation. Validation of these loci on different global populations is important in determining their clinical relevance, particularly for development of novel drug targets for treating and preventing diabetic dyslipidemia and coronary artery disease (CAD). In an attempt to replicate GWAS findings on a non-European sample, we examined the role of six of these loci (CELSR2-PSRC1-SORT1 rs599839; CDKN2A-2B rs1333049; BUD13-ZNF259 rs964184; ZNF259 rs12286037; CETP rs3764261; APOE-C1-C4-C2 rs4420638) in our Asian Indian cohort from the Sikh Diabetes Study (SDS) comprising 3,781 individuals (2,902 from Punjab and 879 from the US). Two of the six SNPs examined showed convincing replication in these populations of Asian Indian origin. Our study confirmed a strong association of CETP rs3764261 with high-density lipoprotein cholesterol (HDL-C) (p?=?2.03×10(-26)). Our results also showed significant associations of two GWAS SNPs (rs964184 and rs12286037) from BUD13-ZNF259 near the APOA5-A4-C3-A1 genes with triglyceride (TG) levels in this Asian Indian cohort (rs964184: p?=?1.74×10(-17); rs12286037: p?=?1.58×10(-2)). We further explored 45 SNPs in a ～195 kb region within the chromosomal region 11q23.3 (encompassing the BUD13-ZNF259, APOA5-A4-C3-A1, and SIK3 genes) in 8,530 Asian Indians from the London Life Sciences Population (LOLIPOP) (UK) and SDS cohorts. Five more SNPs revealed significant associations with TG in both cohorts individually as well as in a joint meta-analysis. However, the strongest signal for TG remained with BUD13-ZNF259 (rs964184: p?=?1.06×10(-39)). Future targeted deep sequencing and functional studies should enhance our understanding of the clinical relevance of these genes in dyslipidemia and hypertriglyceridemia (HTG) and, consequently, diabetes and CAD.     

Associations between single nucleotide polymorphisms in iron-related genes and iron status in multiethnic populations

The existence of multiple inherited disorders of iron metabolism suggests genetic contributions to iron deficiency. We previously performed a genome-wide association study of iron-related single nucleotide polymorphisms (SNPs) using DNA from white men aged ≥ 25 y and women ≥ 50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF) ≤ 12 μg/L (cases) and controls (SF >100 μg/L in men, SF >50 μg/L in women). We report a follow-up study of white, African-American, Hispanic, and Asian HEIRS participants, analyzed for association between SNPs and eight iron-related outcomes. Three chromosomal regions showed association across multiple populations, including SNPs in the TF and TMPRSS6 genes, and on chromosome 18q21. A novel SNP rs1421312 in TMPRSS6 was associated with serum iron in whites (p = 3.7 × 10(-6)) and replicated in African Americans (p = 0.0012).Twenty SNPs in the TF gene region were associated with total iron-binding capacity in whites (p<4.4 × 10(-5)); six SNPs replicated in other ethnicities (p<0.01). SNP rs10904850 in the CUBN gene on 10p13 was associated with serum iron in African Americans (P = 1.0 × 10(-5)). These results confirm known associations with iron measures and give unique evidence of their role in different ethnicities, suggesting origins in a common founder.     

Genome-wide associated loci influencing interleukin (IL)-10, IL-1Ra, and IL-6 levels in African Americans

Interleukins (ILs) are key mediators of the immune response and inflammatory process. Plasma levels of IL-10, IL-1Ra, and IL-6 are associated with metabolic conditions, show large inter-individual variations, and are under strong genetic control. Therefore, elucidation of the genetic variants that influence levels of these ILs provides useful insights into mechanisms of immune response and pathogenesis of diseases. We conducted a genome-wide association study (GWAS) of IL-10, IL-1Ra, and IL-6 levels in 707 non-diabetic African Americans using 5,396,780 imputed and directly genotyped single nucleotide polymorphisms (SNPs) with adjustment for gender, age, and body mass index. IL-10 levels showed genome-wide significant associations (p < 5 × 10(-8)) with eight SNPs, the most significant of which was rs5743185 in the PMS1 gene (p = 2.30 × 10(-10)). We tested replication of SNPs that showed genome-wide significance in 425 non-diabetic individuals from West Africa, and successfully replicated rs17365948 in the YWHAZ gene (p = 0.02). IL-1Ra levels showed suggestive associations with two SNPs in the ASB3 gene (p = 2.55 × 10(-7)), ten SNPs in the IL-1 gene family (IL1F5, IL1F8, IL1F10, and IL1Ra, p = 1.04 × 10(-6) to 1.75 × 10(-6)), and 23 SNPs near the IL1A gene (p = 1.22 × 10(-6) to 1.63 × 10(-6)). We also successfully replicated rs4251961 (p = 0.009); this SNP was reported to be associated with IL-1Ra levels in a candidate gene study of Europeans. IL-6 levels showed genome-wide significant association with one SNP (RP11-314E23.1; chr6:133397598; p = 8.63 × 10(-9)). To our knowledge, this is the first GWAS on IL-10, IL-1Ra, and IL-6 levels. Follow-up of these findings may provide valuable insight into the pathobiology of IL actions and dysregulations in inflammation and human diseases.     

Evidence for a quantitative trait locus affecting low levels of apolipoprotein B and low density lipoprotein on chromosome 10 in Caucasian families

High plasma apolipoprotein B (apoB) and LDL cholesterol levels increase cardiovascular disease risk. These highly correlated measures may be partially controlled by common genetic polymorphisms. To identify chromosomal regions that contain genes causing low plasma levels of one or both parameters in Caucasian families ascertained for familial hypobetalipoproteinemia (FHBL), we conducted a whole-genome scan using 443 microsatellite markers typed in nine multigenerational families with at least two members with FHBL. Both variance components and regression-based linkage methods were used to identify regions of interest. Common linkage regions were identified for both measures on chromosomes 10q25.1-10q26.11 [maximum log of the odds (LOD) = 4.2 for LDL and 3.5 for apoB] and 6q24.3 (maximum LOD = 1.46 for LDL and 1.84 for apoB). There was also evidence for linkage to apoB on chromosome 13q13.2 (LOD = 1.97) and to LDL on chromosome 3p14.1 at 94 centimorgan (LOD = 1.52). Bivariate linkage analysis provided further evidence for loci contributing to both traits (6q24.3, LOD = 1.43; 10q25.1, LOD = 1.74). We evaluated single nucleotide polymorphisms (SNPs) in genes within our linkage regions to identify variants associated with apoB or LDL levels. The most significant finding was for rs2277205 in the 5' untranslated region of acyl-coenzyme A dehydrogenase short/branched chain and LDL (P = 10(-7)). Three additional SNPs were associated with apoB and/or LDL (P < 0.01). Although only the linkage signal on chromosome 10 reached genome-wide statistical significance, there are likely multiple chromosomal regions with variants that contribute to low levels of apoB and LDL and that may protect against coronary heart disease.     

Genome-wide association of BMI in African Americans

Recent genome-wide association studies (GWAS) have identified multiple novel loci associated with obesity in Europeans but results in other ethnicities are less convincing. Here, we report a two-stage GWAS of BMI in African Americans. The GWAS was performed using the Affymetrix 6.0 platform in 816 nondiabetic and 899 diabetic nephropathy subjects. 746,626 single-nucleotide polymorphisms (SNPs) were tested for association with BMI after adjustment for age, gender, disease status, and population structure. Sixty high scoring SNPs that showed nominal association in both GWAS cohorts were further replicated in 3,274 additional subjects in four replication cohorts and a meta-analysis was computed. Meta-analysis of 4,989 subjects revealed five SNPs (rs6794092, rs268972, rs2033195, rs815611, and rs6088887) at four loci showing consistent associations in both GWAS (P < 0.0001) and replication cohorts (P < 0.05) with combined P values range from 2.4 × 10(-6) to 5 × 10(-5). These loci are located near PP13439-TMEM212, CDH12, MFAP3-GALNT10, and FER1L4 and had effect sizes between 0.091 and 0.167 s.d. unit (or 0.67-1.24 kg/m(2)) of BMI for each copy of the effect allele. Our findings suggest the presence of novel loci potentially associated with adiposity in African Americans. Further replication and meta-analysis in African Americans and other populations will shed light on the role of these loci in different ethnic populations.     

GWAS of follicular lymphoma reveals allelic heterogeneity at 6p21.32 and suggests shared genetic susceptibility with diffuse large B-cell lymphoma

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL-associated locus on 6p21.32, rs2647012 (OR(combined) = 0.64, P(combined) = 2 × 10(-21)) located 962 bp away from rs10484561 (r(2)<0.1 in controls). After mutual adjustment, the associations at the two SNPs remained genome-wide significant (rs2647012:OR(adjusted) = 0.70, P(adjusted) = 4 × 10(-12); rs10484561:OR(adjusted) = 1.64, P(adjusted) = 5 × 10(-15)). Haplotype and coalescence analyses indicated that rs2647012 arose on an evolutionarily distinct haplotype from that of rs10484561 and tags a novel allele with an opposite (protective) effect on FL risk. Moreover, in a follow-up analysis of the top 6 FL-associated SNPs in 4,449 cases of other NHL subtypes, rs10484561 was associated with risk of diffuse large B-cell lymphoma (OR(combined) = 1.36, P(combined) = 1.4 × 10(-7)). Our results reveal the presence of allelic heterogeneity within the HLA class II region influencing FL susceptibility and indicate a possible shared genetic etiology with diffuse large B-cell lymphoma. These findings suggest that the HLA class II region plays a complex yet important role in NHL.     

Association of smoking behavior with an odorant receptor allele telomeric to the human major histocompatibility complex

Smoking behavior has been associated in two independent European cohorts with the most common Caucasian human leukocyte antigen (HLA) haplotype (A1-B8-DR3). We aimed to test whether polymorphic members of the two odorant receptor (OR) clusters within the extended HLA complex might be responsible for the observed association, by genotyping a cohort of Hungarian women in which the mentioned association had been found. One hundred and eighty HLA haplotypes from Centre d'Etude du Polymorphisme Humain families were analyzed in silico to identify single-nucleotide polymorphisms (SNPs) within OR genes that are in linkage disequilibrium with the A1-B8-DR3 haplotype, as well as with two other haplotypes indirectly linked to smoking behavior. A nonsynonymous SNP within the OR12D3 gene (rs3749971(T)) was found to be linked to the A1-B8-DR3 haplotype. This polymorphism leads to a (97)Thr --> Ile exchange that affects a putative ligand binding region of the OR12D3 protein. Smoking was found to be associated in the Hungarian cohort with the rs3749971(T) allele (p = 1.05 x 10(-2)), with higher significance than with A1-B8-DR3 (p = 2.38 x 10(-2)). Our results link smoking to a distinct OR allele, and demonstrate that the rs3749971(T) polymorphism is associated with the HLA haplotype-dependent differential recognition of cigarette smoke components, at least among Caucasian women.     

Regional heritability mapping method helps explain missing heritability of blood lipid traits in isolated populations

Single single-nucleotide polymorphism (SNP) genome-wide association studies (SSGWAS) may fail to identify loci with modest effects on a trait. The recently developed regional heritability mapping (RHM) method can potentially identify such loci. In this study, RHM was compared with the SSGWAS for blood lipid traits (high-density lipoprotein (HDL), low-density lipoprotein (LDL), plasma concentrations of total cholesterol (TC) and triglycerides (TG)). Data comprised 2246 adults from isolated populations genotyped using ∼300 000 SNP arrays. The results were compared with large meta-analyses of these traits for validation. Using RHM, two significant regions affecting HDL on chromosomes 15 and 16 and one affecting LDL on chromosome 19 were identified. These regions covered the most significant SNPs associated with HDL and LDL from the meta-analysis. The chromosome 19 region was identified in our data despite the fact that the most significant SNP in the meta-analysis (or any SNP tagging it) was not genotyped in our SNP array. The SSGWAS identified one SNP associated with HDL on chromosome 16 (the top meta-analysis SNP) and one on chromosome 10 (not reported by RHM or in the meta-analysis and hence possibly a false positive association). The results further confirm that RHM can have better power than SSGWAS in detecting causal regions including regions containing crucial ungenotyped variants. This study suggests that RHM can be a useful tool to explain some of the ‘missing heritability'' of complex trait variation.     

Use of a multiethnic approach to identify rheumatoid- arthritis-susceptibility loci, 1p36 and 17q12

We have previously shown that rheumatoid arthritis (RA) risk alleles overlap between different ethnic groups. Here, we utilize a multiethnic approach to show that we can effectively discover RA risk alleles. Thirteen putatively associated SNPs that had not yet exceeded genome-wide significance (p < 5 × 10(-8)) in our previous RA genome-wide association study (GWAS) were analyzed in independent sample sets consisting of 4,366 cases and 17,765 controls of European, African American, and East Asian ancestry. Additionally, we conducted an overall association test across all 65,833 samples (a GWAS meta-analysis plus the replication samples). Of the 13 SNPs investigated, four were significantly below the study-wide Bonferroni corrected p value threshold (p < 0.0038) in the replication samples. Two SNPs (rs3890745 at the 1p36 locus [p = 2.3 × 10(-12)] and rs2872507 at the 17q12 locus [p = 1.7 × 10(-9)]) surpassed genome-wide significance in all 16,659 RA cases and 49,174 controls combined. We used available GWAS data to fine map these two loci in Europeans and East Asians, and we found that the same allele conferred risk in both ethnic groups. A series of bioinformatic analyses identified TNFRSF14-MMEL1 at the 1p36 locus and IKZF3-ORMDL3-GSDMB at the 17q12 locus as the genes most likely associated with RA. These findings demonstrate empirically that a multiethnic approach is an effective strategy for discovering RA risk loci, and they suggest that combining GWASs across ethnic groups represents an efficient strategy for gaining statistical power.