Michèle M. Sale

Robust relationship inference in genome-wide association studies

MOTIVATION Genome-wide association studies (GWASs) have been widely used to map loci contributing to variation in complex traits and risk of diseases in humans. Accurate specification of familial relationships is crucial for family-based GWAS, as well as in population-based GWAS with unknown (or unrecognized) family structure. The family structure in a GWAS should be routinely investigated using the SNP data prior to the analysis of population structure or phenotype. Existing algorithms for relationship inference have a major weakness of estimating allele frequencies at each SNP from the entire sample, under a strong assumption of homogeneous population structure. This assumption is often untenable. RESULTS Here, we present a rapid algorithm for relationship inference using high-throughput genotype data typical of GWAS that allows the presence of unknown population substructure. The relationship of any pair of individuals can be precisely inferred by robust estimation of their kinship coefficient, independent of sample composition or population structure (sample invariance). We present simulation experiments to demonstrate that the algorithm has sufficient power to provide reliable inference on millions of unrelated pairs and thousands of relative pairs (up to 3rd-degree relationships). Application of our robust algorithm to HapMap and GWAS datasets demonstrates that it performs properly even under extreme population stratification, while algorithms assuming a homogeneous population give systematically biased results. Our extremely efficient implementation performs relationship inference on millions of pairs of individuals in a matter of minutes, dozens of times faster than the most efficient existing algorithm known to us. AVAILABILITY Our robust relationship inference algorithm is implemented in a freely available software package, KING, available for download at http://people.virginia.edu/∼wc9c/KING.

Genetic risk factors for ischaemic stroke and its subtypes (the METASTROKE collaboration): a meta-analysis of genome-wide association studies.

Summary Background Various genome-wide association studies (GWAS) have been done in ischaemic stroke, identifying a few loci associated with the disease, but sample sizes have been 3500 cases or less. We established the METASTROKE collaboration with the aim of validating associations from previous GWAS and identifying novel genetic associations through meta-analysis of GWAS datasets for ischaemic stroke and its subtypes. Methods We meta-analysed data from 15 ischaemic stroke cohorts with a total of 12 389 individuals with ischaemic stroke and 62 004 controls, all of European ancestry. For the associations reaching genome-wide significance in METASTROKE, we did a further analysis, conditioning on the lead single nucleotide polymorphism in every associated region. Replication of novel suggestive signals was done in 13 347 cases and 29 083 controls. Findings We verified previous associations for cardioembolic stroke near PITX2 (p=2·8×10−16) and ZFHX3 (p=2·28×10−8), and for large-vessel stroke at a 9p21 locus (p=3·32×10−5) and HDAC9 (p=2·03×10−12). Additionally, we verified that all associations were subtype specific. Conditional analysis in the three regions for which the associations reached genome-wide significance (PITX2, ZFHX3, and HDAC9) indicated that all the signal in each region could be attributed to one risk haplotype. We also identified 12 potentially novel loci at p<5×10−6. However, we were unable to replicate any of these novel associations in the replication cohort. Interpretation Our results show that, although genetic variants can be detected in patients with ischaemic stroke when compared with controls, all associations we were able to confirm are specific to a stroke subtype. This finding has two implications. First, to maximise success of genetic studies in ischaemic stroke, detailed stroke subtyping is required. Second, different genetic pathophysiological mechanisms seem to be associated with different stroke subtypes. Funding Wellcome Trust, UK Medical Research Council (MRC), Australian National and Medical Health Research Council, National Institutes of Health (NIH) including National Heart, Lung and Blood Institute (NHLBI), the National Institute on Aging (NIA), the National Human Genome Research Institute (NHGRI), and the National Institute of Neurological Disorders and Stroke (NINDS).

Linkage of Genetic Markers on Human Chromosomes 20 and 12 to NIDDM in Caucasian Sib Pairs With a History of Diabetic Nephropathy

The potential contribution of maturity-onset diabetes of the young (MODY) genes to NIDDM susceptibility in African-American and Caucasian NIDDM-affected sibling pairs with a history of adult-onset diabetic nephropathy has been evaluated. Evidence for linkage to NIDDM was found with polymorphic loci that map to the long arms of human chromosomes 20 and 12 in regions containing the MODY1 and MODY3 genes. Nonparametric analysis of chromosome 20 inheritance data collected with the MODYl-linked marker D20S197 provides evidence forlinkage to NIDDM with a P value of 0.005 in Caucasian sib pairs using affected sibpair (ASP) analyses. Nonparametric analysis of chromosome 12 inheritance data collected with the MODY3-linked markers D12S349 and D12S86 provides evidence for linkage to NIDDM with P values of 0.04 and 0.006, respectively, in Caucasian sib pairs using similar analyses. No evidence for linkage of MODY1 and MODY3 markers to NIDDM in African-American sib pairs was observed. In addition, no evidence for linkage to MODY2 (glucokinase-associated MODY) was observed with either study population. Results of multipoint maximum logarithm of odds (LOD) score analysis were consistent with the ASP results. A maximum LOD score of 1.48 was calculated for linkage to MODYl-linked loci and 1.45 to MODY3-linked loci in Caucasian sib pairs. Tabulation of allele sharing in affected sib pairs with D20S197 and D12S349 suggests that affected sibling pairs may inherit susceptibility genes simultaneously from chromosome 20 and chromosome 12. The results suggest that genes contributing to NIDDM in the general Caucasian population are located in the regions containing the MODY1 and MODY3 genes.

Association Analysis in African Americans of European-Derived Type 2 Diabetes Single Nucleotide Polymorphisms From Whole-Genome Association Studies

OBJECTIVE— Several whole-genome association studies have reported identification of type 2 diabetes susceptibility genes in various European-derived study populations. Little investigation of these loci has been reported in other ethnic groups, specifically African Americans. Striking differences exist between these populations, suggesting they may not share identical genetic risk factors. Our objective was to examine the influence of type 2 diabetes genes identified in whole-genome association studies in a large African American case-control population. RESEARCH DESIGN AND METHODS— Single nucleotide polymorphisms (SNPs) in 12 loci (e.g., TCF7L2, IDE/KIF11/HHEX, SLC30A8, CDKAL1, PKN2, IGF2BP2, FLJ39370, and EXT2/ALX4) associated with type 2 diabetes in European-derived populations were genotyped in 993 African American type 2 diabetic and 1,054 African American control subjects. Additionally, 68 ancestry-informative markers were genotyped to account for the impact of admixture on association results. RESULTS— Little evidence of association was observed between SNPs, with the exception of those in TCF7L2, and type 2 diabetes in African Americans. One TCF7L2 SNP (rs7903146) showed compelling evidence of association with type 2 diabetes (admixture-adjusted additive P [Pa] = 1.59 × 10−6). Only the intragenic SNP on 11p12 (rs9300039, dominant P [Pd] = 0.029) was also associated with type 2 diabetes after admixture adjustments. Interestingly, four of the SNPs are monomorphic in the Yoruba population of the HAPMAP project, with only the risk allele from the populations of European descent present. CONCLUSIONS— Results suggest that these variants do not significantly contribute to interindividual susceptibility to type 2 diabetes in African Americans. Consequently, genes contributing to type 2 diabetes in African Americans may, in part, be different from those in European-derived study populations. High frequency of risk alleles in several of these genes may, however, contribute to the increased prevalence of type 2 diabetes in African Americans.

A GENOME WIDE ASSOCIATION STUDY FOR DIABETIC NEPHROPATHY GENES IN AFRICAN AMERICANS

A genome-wide association study was performed using the Affymetrix 6.0 chip to identify genes associated with diabetic nephropathy in African Americans. Association analysis was performed adjusting for admixture in 965 type 2 diabetic African American patients with end-stage renal disease (ESRD) and in 1029 African Americans without type 2 diabetes or kidney disease as controls. The top 724 single nucleotide polymorphisms (SNPs) with evidence of association to diabetic nephropathy were then genotyped in a replication sample of an additional 709 type 2 diabetes-ESRD patients and 690 controls. SNPs with evidence of association in both the original and replication studies were tested in additional African American cohorts consisting of 1246 patients with type 2 diabetes without kidney disease and 1216 with non-diabetic ESRD to differentiate candidate loci for type 2 diabetes-ESRD, type 2 diabetes, and/or all-cause ESRD. Twenty-five SNPs were significantly associated with type 2 diabetes-ESRD in the genome-wide association and initial replication. Although genome-wide significance with type 2 diabetes was not found for any of these 25 SNPs, several genes, including RPS12, LIMK2, and SFI1 are strong candidates for diabetic nephropathy. A combined analysis of all 2890 patients with ESRD showed significant association SNPs in LIMK2 and SFI1 suggesting that they also contribute to all-cause ESRD. Thus, our results suggest that multiple loci underlie susceptibility to kidney disease in African Americans with type 2 diabetes and some may also contribute to all-cause ESRD.

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10−6 were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10−8 for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10−8 for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10−8; RREB1: p = 5.7×10−8). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.

Genetic dissection of the human leukocyte antigen region by use of haplotypes of Tasmanians with multiple sclerosis.

Association of multiple sclerosis (MS) with the human leukocyte antigen (HLA) class II haplotype DRB1*1501-DQB1*0602 is the most consistently replicated finding of genetic studies of the disease. However, the high level of linkage disequilibrium (LD) in the HLA region has hindered the identification of other loci that single-marker tests for association are unlikely to resolve. In order to address this issue, we generated haplotypes spanning 14.754 Mb (5 cM) across the entire HLA region. The haplotypes, which were inferred by genotyping relatives of 152 patients with MS and 105 unaffected control subjects of Tasmanian ancestry, define a genomic segment from D6S276 to D6S291, including 13 microsatellite markers integrated with allele-typing data for DRB1 and DQB1. Association to the DRB1*1501-DQB1*0602 haplotype was replicated. In addition, we found that the class I/extended class I region, defined by a genomic segment of approximately 400 kb between MOGCA and D6S265, harbors genes that independently increase risk of, or provide protection from, MS. Log-linear modeling analysis of constituent haplotypes that represent genomic regions containing class I (MOGCA-D6S265), class III (TNFa-TNFd-D6S273), and class II (DRB1-DQB1) genes indicated that having class I and class II susceptibility variants on the same haplotype provides an additive effect on risk. Moreover, we found no evidence for a disease locus in the class III region defined by a 150-kb genomic segment containing the TNF locus and 14 other genes. A global overview of LD performed using GOLD identified two discrete blocks of LD in the HLA region that correspond well with previous findings. We propose that the analysis of haplotypes, by use of the types of approaches outlined in the present article, should make it possible to more accurately define the contribution of the HLA to MS.

Mutations at the BLK locus linked to maturity onset diabetes of the young and β-cell dysfunction

Maturity-onset diabetes of the young (MODY) is a subtype of diabetes defined by an autosomal pattern of inheritance and a young age at onset, often before age 25. MODY is genetically heterogeneous, with 8 distinct MODY genes identified to date and more believed to exist. We resequenced 732 kb of genomic sequence at 8p23 in 6 MODY families unlinked to known MODY genes that showed evidence of linkage at that location. Of the 410 sequence differences that we identified, 5 had a frequency <1% in the general population and segregated with diabetes in 3 of the families, including the 2 showing the strongest support for linkage at this location. The 5 mutations were all placed within 100 kb corresponding to the BLK gene. One resulted in an Ala71Thr substitution; the other 4 were noncoding and determined decreased in vitro promoter activity in reporter gene experiments. We found that BLK—a nonreceptor tyrosine-kinase of the src family of proto-oncogenes—is expressed in β-cells where it enhances insulin synthesis and secretion in response to glucose by up-regulating transcription factors Pdx1 and Nkx6.1. These actions are greatly attenuated by the Ala71Thr mutation. These findings point to BLK as a previously unrecognized modulator of β-cell function, the deficit of which may lead to the development of diabetes.

Variants in Intron 13 of the ELMO1 Gene are Associated with Diabetic Nephropathy in African Americans

Variants in the engulfment and cell motility 1 (ELMO1) gene are associated with nephropathy due to type 2 diabetes mellitus (T2DM) in a Japanese cohort. We comprehensively evaluated this gene in African American (AA) T2DM patients with end‐stage renal disease (ESRD). Three hundred and nine HapMap tagging SNPs and 9 reportedly associated SNPs were genotyped in 577 AA T2DM‐ESRD patients and 596 AA non‐diabetic controls, plus 43 non‐diabetic European American controls and 45 Yoruba Nigerian samples for admixture adjustment. Replication analyses were conducted in 558 AA with T2DM‐ESRD and 564 controls without diabetes. Extension analyses included 328 AA with T2DM lacking nephropathy and 326 with non‐diabetic ESRD. The original and replication analyses confirmed association with four SNPs in intron 13 (permutation p‐values for combined analyses = 0.001–0.003), one in intron 1 (P = 0.004) and one in intron 5 (P = 0.002) with T2DM‐associated ESRD. In a subsequent combined analysis of all 1,135 T2DM‐ESRD cases and 1,160 controls, an additional 7 intron 13 SNPs produced evidence of association (P = 3.5 × 10−5– P = 0.05). No associations were seen with these SNPs in those with T2DM lacking nephropathy or with ESRD due to non‐diabetic causes. Variants in intron 13 of the ELMO1 gene appear to confer risk for diabetic nephropathy in AA.

Variants of the Transcription Factor 7-Like 2 (TCF7L2) Gene Are Associated With Type 2 Diabetes in an African-American Population Enriched for Nephropathy

OBJECTIVE—Recently, variants in the TCF7L2 gene have been reported to be associated with type 2 diabetes across multiple Europid populations, but only one small sample of African-American type 2 diabetic patients has been examined. Our objective was to investigate the importance of TCF7L2 in a larger African-American case-control population. RESEARCH DESIGN AND METHODS—We investigated single nucleotide polymorphisms (SNPs) in six known type 2 diabetes genes in 577 African-American case subjects with type 2 diabetes enriched for nephropathy and 596 African-American control subjects. Additionally, we genotyped 70 ancestry-informative markers (AIMs) to apply adjustments for differences in ancestral proportions. RESULTS—The most significant associations were observed with TCF7L2 intron 3 SNPs rs7903146 (additive P = 4.10 × 10−6, odds ratio [OR] 1.51; admixture-adjusted Pa = 3.77 × 10−6) and rs7901695 (P = 0.001, OR 1.30; Pa = 0.003). The 2-SNP haplotype containing these SNPs was also associated with type 2 diabetes (P = 3 × 10−5). Modest associations were also seen with TCF7L2 intron 4 SNPs rs7895340, rs11196205, and rs12255372 (0.01 < P < 0.05; 0.03 < Pa < 0.08), as well as with ATP-sensitive inwardly rectifying potassium channel subunit Kir6.2 (KCNJ11) and hepatocyte nuclear factor 4-α (HNF4A) SNPs (0.01 < P < 0.05; 0.01 < Pa < 0.41). No significant associations were detected with genotyped calpain 10 (CAPN10), peroxisome proliferator–activated receptor γ (PPARG), and transcription factor 1 (TCF1) SNPs. CONCLUSIONS—This study indicates that variants in the TCF7L2 gene significantly contribute to diabetes susceptibility in African-American populations.