The Type 1 Diabetes Genetics Consortium (T1DGC) has assembled and genotyped a large collection of multiplex families for the purpose of mapping genomic regions linked to type 1 diabetes. In the current study, we tested for evidence of loci associated with type 1 diabetes utilizing genome-wide linkage scan data and family-based association methods.

Admixture mapping recently identified MYH9 as a susceptibility gene for idiopathic focal segmental glomerulosclerosis (FSGS), HIV-associated nephropathy (HIVAN) and end-stage kidney disease attributed to hypertension (H-ESKD) in African Americans (AA). MYH9 encodes the heavy chain of non-muscle myosin IIA, a cellular motor involved in motility. A haplotype and its tagging SNPs spanning introns 12-23 were most strongly associated with kidney disease (OR 2-7; P < 10(-8), recessive). To narrow the region of association and identify potential causal variation, we performed a dense-mapping study using 79 MYH9 SNPs in AA populations with FSGS, HIVAN and H-ESKD (typed for a subset of 46 SNPs), for a total of 2496 cases and controls. The strongest associations were for correlated SNPs rs5750250, rs2413396 and rs5750248 in introns 13, 14 and 15, a region of 5.6 kb. Rs5750250 showed OR 5.0, 8.0 and 2.8; P = 2 x 10(-17), 2 x 10(-10) and 3 x 10(-22), respectively, for FSGS, HIVAN and H-ESKD; OR 5.7; P = 9 x 10(-27) for combined FSGS and HIVAN, recessive. An independent association was observed for rs11912763 in intron 33. Neither the highly associated SNPs nor the results of resequencing MYH9 in 40 HIVAN or FSGS cases and controls revealed non-synonymous changes that could account for the disease associations. Rs2413396 and one of the highly associated SNPs in intron 23, rs4821480, are predicted splicing motif modifiers. Rs5750250 combined with rs11912763 had receiver operator characteristic (ROC) C statistics of 0.80, 0.73 and 0.65 for HIVAN, FSGS and H-ESKD, respectively, allowing prediction of genetic risk by typing two SNPs.

Blood specimen collection at an early study visit is often included in observational studies or clinical trials for analysis of secondary outcome biomarkers. A common protocol is to store buffy coat specimens for future DNA isolation and these may remain in frozen storage for many years. It is uncertain if the DNA remains suitable for modern genome wide association (GWA) genotyping.

A genome-wide association scan of the Genetics of Kidneys in Diabetes (GoKinD) collections identified four novel susceptibility loci, located on chromosomes 7p14.3, 9q21.32, 11p15.4, and 13q33.3 associated with type 1 diabetic nephropathy. A recent evaluation of these loci in Japanese patients with type 2 diabetes supported an association at the 13q33.3 locus. To follow up these findings, we determined whether single-nucleotide polymorphisms (SNPs) at these same four loci were associated with diabetic nephropathy in the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes collection. A total of 6 SNPs across these loci were genotyped in 646 normoalbuminuric controls and in 743 nephropathy patients of European ancestry. A significant association was identified at the 13q33.3 locus (rs9521445: P = 4.4 × 10(-3)). At this same locus, rs1411766 was also significantly associated with type 2 diabetic nephropathy (P = 0.03). Meta-analysis of these data with those of the Japanese and GoKinD collections significantly improved the strength of the association (P = 9.7 × 10(-9)). In addition, there was a significant association at the 11p15.4 locus (rs451041: P = 0.02). Thus, associations identified in the GoKinD collections on chromosomes 11p15.4 (near the CARS gene) and 13q33.3 (within an intergenic region between MYO16 and IRS2) are susceptibility loci of kidney disease common to both type 1 and 2 diabetes.

To discover new rheumatoid arthritis (RA) risk loci, we systematically examined 370 SNPs from 179 independent loci with P < 0.001 in a published meta-analysis of RA genome-wide association studies (GWAS) of 3,393 cases and 12,462 controls. We used Gene Relationships Across Implicated Loci (GRAIL), a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three were convincingly validated: CD2-CD58 (rs11586238, P = 1 x 10(-6) replication, P = 1 x 10(-9) overall), CD28 (rs1980422, P = 5 x 10(-6) replication, P = 1 x 10(-9) overall) and PRDM1 (rs548234, P = 1 x 10(-5) replication, P = 2 x 10(-8) overall). An additional four were replicated (P < 0.0023): TAGAP (rs394581, P = 0.0002 replication, P = 4 x 10(-7) overall), PTPRC (rs10919563, P = 0.0003 replication, P = 7 x 10(-7) overall), TRAF6-RAG1 (rs540386, P = 0.0008 replication, P = 4 x 10(-6) overall) and FCGR2A (rs12746613, P = 0.0022 replication, P = 2 x 10(-5) overall). Many of these loci are also associated to other immunologic diseases.

Case-control genetic studies of complex human diseases can be confounded by population stratification. This issue can be addressed using panels of ancestry informative markers (AIMs) that can provide substantial population substructure information. Previously, we described a panel of 128 SNP AIMs that were designed as a tool for ascertaining the origins of subjects from Europe, Sub-Saharan Africa, Americas, and East Asia.

Accounting for population genetic substructure is important in reducing type 1 errors in genetic studies of complex disease. As efforts to understand complex genetic disease are expanded to different continental populations the understanding of genetic substructure within these continents will be useful in design and execution of association tests. In this study, population differentiation (Fst) and Principal Components Analyses (PCA) are examined using >200 K genotypes from multiple populations of East Asian ancestry. The population groups included those from the Human Genome Diversity Panel [Cambodian, Yi, Daur, Mongolian, Lahu, Dai, Hezhen, Miaozu, Naxi, Oroqen, She, Tu, Tujia, Naxi, Xibo, and Yakut], HapMap [ Han Chinese (CHB) and Japanese (JPT)], and East Asian or East Asian American subjects of Vietnamese, Korean, Filipino and Chinese ancestry. Paired Fst (Wei and Cockerham) showed close relationships between CHB and several large East Asian population groups (CHB/Korean, 0.0019; CHB/JPT, 00651; CHB/Vietnamese, 0.0065) with larger separation with Filipino (CHB/Filipino, 0.014). Low levels of differentiation were also observed between Dai and Vietnamese (0.0045) and between Vietnamese and Cambodian (0.0062). Similarly, small Fst's were observed among different presumed Han Chinese populations originating in different regions of mainland of China and Taiwan (Fst's <0.0025 with CHB). For PCA, the first two PC's showed a pattern of relationships that closely followed the geographic distribution of the different East Asian populations. PCA showed substructure both between different East Asian groups and within the Han Chinese population. These studies have also identified a subset of East Asian substructure ancestry informative markers (EASTASAIMS) that may be useful for future complex genetic disease association studies in reducing type 1 errors and in identifying homogeneous groups that may increase the power of such studies.

There is increasing evidence from genome-wide association studies for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children, yet the effects of protein-coding variants on ALL risk have not been systematically evaluated. Here we show a missense variant in CDKN2A associated with the development of ALL at genome-wide significance (rs3731249, P=9.4 × 10(-23), odds ratio=2.23). Functional studies indicate that this hypomorphic variant results in reduced tumour suppressor function of p16(INK4A), increases the susceptibility to leukaemic transformation of haematopoietic progenitor cells, and is preferentially retained in ALL tumour cells. Resequencing the CDKN2A-CDKN2B locus in 2,407 childhood ALL cases reveals 19 additional putative functional germline variants. These results provide direct functional evidence for the influence of inherited genetic variation on ALL risk, highlighting the important and complex roles of CDKN2A-CDKN2B tumour suppressors in leukaemogenesis.

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of The Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value <6.61 × 10(-4)), including a rare missense variant in ABCA6 (rs77542162, p.Cys1359Arg, frequency 0.034), which is predicted to be deleterious. The frequency of this ABCA6 variant is 3.65-fold increased in the Dutch and its effect (βLDL-C=0.135, βTC=0.140) is estimated to be very similar to those observed for single variants in well-known lipid genes, such as LDLR.

Patients with type 2 diabetes are at increased risk of age-related cognitive decline and dementia. Neuroimaging measures such as white matter lesion volume, brain volume, and fractional anisotropy may reflect the pathogenesis of these cognitive declines, and genetic factors may contribute to variability in these measures. This study examined multiple neuroimaging measures in 465 participants from 238 families with extensive genotype data in the type 2 diabetes enriched Diabetes Heart Study-Mind cohort. Heritability of these phenotypes and their association with candidate single-nucleotide polymorphisms (SNPs), and SNP data from genome- and exome-wide arrays were explored. All neuroimaging measures analyzed were significantly heritable (ĥ(2) = 0.55-0.99 in unadjusted models). Seventeen candidate SNPs (from 16 genes/regions) associated with neuroimaging phenotypes in prior studies showed no significant evidence of association. A missense variant (rs150706952, A432V) in PLEKHG4B from the exome-wide array was significantly associated with white matter mean diffusivity (p = 3.66 × 10(-7)) and gray matter mean diffusivity (p = 2.14 × 10(-7)). This analysis suggests genetic factors contribute to variation in neuroimaging measures in a population enriched for metabolic disease and other associated comorbidities.

Resistant hypertension (RHTN), defined by lack of blood pressure (BP) control despite treatment with at least 3 antihypertensive drugs, increases cardiovascular risk compared with controlled hypertension. Yet, there are few data on genetic variants associated with RHTN.

Copy number variants (CNVs) have been proposed as a possible source of 'missing heritability' in complex human diseases. Two studies of type 1 diabetes (T1D) found null associations with common copy number polymorphisms, but CNVs of low frequency and high penetrance could still play a role. We used the Log-R-ratio intensity data from a dense single nucleotide polymorphism (SNP) array, ImmunoChip, to detect rare CNV deletions (rDELs) and duplications (rDUPs) in 6808 T1D cases, 9954 controls and 2206 families with T1D-affected offspring. Initial analyses detected CNV associations. However, these were shown to be false-positive findings, failing replication with polymerase chain reaction. We developed a pipeline of quality control (QC) tests that were calibrated using systematic testing of sensitivity and specificity. The case-control odds ratios (OR) of CNV burden on T1D risk resulting from this QC pipeline converged on unity, suggesting no global frequency difference in rDELs or rDUPs. There was evidence that deletions could impact T1D risk for a small minority of cases, with enrichment for rDELs longer than 400 kb (OR = 1.57, P = 0.005). There were also 18 de novo rDELs detected in affected offspring but none for unaffected siblings (P = 0.03). No specific CNV regions showed robust evidence for association with T1D, although frequencies were lower than expected (most less than 0.1%), substantially reducing statistical power, which was examined in detail. We present an R-package, plumbCNV, which provides an automated approach for QC and detection of rare CNVs that can facilitate equivalent analyses of large-scale SNP array datasets.

Advanced glycation end-products (AGEs) are a diverse group of molecules produced by the non-enzymatic addition of glucose to proteins, lipids, and nucleic acids. AGE levels have been associated with hyperglycemia and diabetic complications, especially in animal models, but less clearly in human studies. We measured total serum AGEs using an enzyme linked immunosorbant assay (ELISA) in 506 subjects from 246 families in the Diabetes Heart Study (DHS)/DHS MIND Study (n=399 type 2 diabetes (T2D)-affected). Single nucleotide polymorphisms (SNPs) in several candidate genes, including known AGE receptors, were tested for their influence on circulating AGE levels. The genetic analysis was expanded to include an exploratory genome-wide association study (GWAS) and exome chip analysis of AGEs (≈440,000 SNPs). AGEs were found to be highly heritable (h(2)=0.628, p=8.96 × 10(-10)). While no SNPs from candidate genes were significantly associated after Bonferroni correction, rs1035798 in the gene AGER was the most significantly associated (p=0.007). Additionally, rs7198427, in MT1A, showed a nominally significant p-value (p=0.0099). No SNPs from the GWAS or exome studies were identified after correction for multiple comparisons; however, rs17054480 in the PALLD2 gene on chromosome 4 showed the strongest association (p=7.77 × 10(-7)). Five SNPs at two loci (ISCA2/NPC2 and FBXO33) had p-values of less than 2.0 × 10(-5) and three additional SNPs (rs716326 in MACROD2, and rs6795197 and rs6765857 in ZBTB38) showed a nominal association with p-values of less than 1.0 × 10(-5).These findings provide a foundation for further investigation into the genetic component of circulating AGEs.

Despite the substantial burden of hypertension in US minority populations, few genetic studies of blood pressure have been conducted in Hispanics and African Americans, and it is unclear whether many of the established loci identified in European-descent populations contribute to blood pressure variation in non-European descent populations. Using the Metabochip array, we sought to characterize the genetic architecture of previously identified blood pressure loci, and identify novel cardiometabolic variants related to systolic and diastolic blood pressure in a multi-ethnic US population including Hispanics (n = 19,706) and African Americans (n = 18,744). Several known blood pressure loci replicated in African Americans and Hispanics. Fourteen variants in three loci (KCNK3, FGF5, ATXN2-SH2B3) were significantly associated with blood pressure in Hispanics. The most significant diastolic blood pressure variant identified in our analysis, rs2586886/KCNK3 (P = 5.2 x 10-9), also replicated in independent Hispanic and European-descent samples. African American and trans-ethnic meta-analysis data identified novel variants in the FGF5, ULK4 and HOXA-EVX1 loci, which have not been previously associated with blood pressure traits. Our identification and independent replication of variants in KCNK3, a gene implicated in primary hyperaldosteronism, as well as a variant in HOTTIP (HOXA-EVX1) suggest that further work to clarify the roles of these genes may be warranted. Overall, our findings suggest that loci identified in European descent populations also contribute to blood pressure variation in diverse populations including Hispanics and African Americans-populations that are understudied for hypertension genetic risk factors.

Massively parallel whole-genome sequencing (WGS) data have ushered in a new era in human genetics. These data are now being used to understand the role of rare variants in complex traits and to advance the goals of precision medicine. The technological and computing advances that have enabled us to generate WGS data on thousands of individuals have also outpaced our ability to perform analyses in scientifically and statistically rigorous and thoughtful ways. The past several years have witnessed the application of whole-exome sequencing (WES) to complex traits and diseases. From our analysis of NHLBI Exome Sequencing Project (ESP) data, not only have a number of important disease and complex trait association findings emerged, but our collective experience offers some valuable lessons for WGS initiatives. These include caveats associated with generating automated pipelines for quality control and analysis of rare variants; the importance of studying minority populations; sample size requirements and efficient study designs for identifying rare-variant associations; and the significance of incidental findings in population-based genetic research. With the ESP as an example, we offer guidance and a framework on how to conduct a large-scale association study in the era of WGS.

Using the Immunochip custom SNP array, which was designed for dense genotyping of 186 loci identified through genome-wide association studies (GWAS), we analyzed 11,475 individuals with rheumatoid arthritis (cases) of European ancestry and 15,870 controls for 129,464 markers. We combined these data in a meta-analysis with GWAS data from additional independent cases (n = 2,363) and controls (n = 17,872). We identified 14 new susceptibility loci, 9 of which were associated with rheumatoid arthritis overall and five of which were specifically associated with disease that was positive for anticitrullinated peptide antibodies, bringing the number of confirmed rheumatoid arthritis risk loci in individuals of European ancestry to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at 6 loci and identified association to low-frequency variants at 4 loci. Bioinformatic analyses generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Acetyl coenzyme A carboxylase B gene (ACACB) single nucleotide polymorphism (SNP) rs2268388 is reproducibly associated with type 2 diabetes (T2DM)-associated nephropathy (DN). ACACB knock-out mice are also protected from obesity. This study assessed relationships between rs2268388, body mass index (BMI) and gene expression in multiple populations, with and without T2DM. Among subjects without T2DM, rs2268388 DN risk allele (T) associated with higher BMI in Pima Indian children (n = 2021; p-additive = 0.029) and African Americans (AAs) (n = 177; p-additive = 0.05), with a trend in European Americans (EAs) (n = 512; p-additive = 0.09), but not Germans (n = 858; p-additive = 0.765). Association with BMI was seen in a meta-analysis including all non-T2DM subjects (n = 3568; p-additive = 0.02). Among subjects with T2DM, rs2268388 was not associated with BMI in Japanese (n = 2912) or EAs (n = 1149); however, the T allele associated with higher BMI in the subset with BMI≥30 kg/m(2) (n = 568 EAs; p-additive = 0.049, n = 196 Japanese; p-additive = 0.049). Association with BMI was strengthened in a T2DM meta-analysis that included an additional 756 AAs (p-additive = 0.080) and 48 Hong Kong Chinese (p-additive = 0.81) with BMI≥30 kg/m(2) (n = 1575; p-additive = 0.0033). The effect of rs2268388 on gene expression revealed that the T risk allele associated with higher ACACB messenger levels in adipose tissue (41 EAs and 20 AAs with BMI>30 kg/m(2); p-additive = 0.018) and ACACB protein levels in the liver tissue (mixed model p-additive = 0.03, in 25 EA bariatric surgery patients with BMI>30 kg/m(2) for 75 exams). The T allele also associated with higher hepatic triglyceride levels. These data support a role for ACACB in obesity and potential roles for altered lipid metabolism in susceptibility to DN.

Genome-wide association studies (GWAS) have identified common variants of modest-effect size at hundreds of loci for common autoimmune diseases; however, a substantial fraction of heritability remains unexplained, to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, most studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach fails to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon-sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 UK residents of white European origin, comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci have a negligible role in common autoimmune disease susceptibility. These results do not support the rare-variant synthetic genome-wide-association hypothesis (in which unobserved rare causal variants lead to association detected at common tag variants). Many known autoimmune disease risk loci contain multiple, independently associated, common and low-frequency variants, and so genes at these loci are a priori stronger candidates for harbouring rare coding-region variants than other genes. Our data indicate that the missing heritability for common autoimmune diseases may not be attributable to the rare coding-region variant portion of the allelic spectrum, but perhaps, as others have proposed, may be a result of many common-variant loci of weak effect.

Narrow arterioles in the retina have been shown to predict hypertension as well as other vascular diseases, likely through an increase in the peripheral resistance of the microcirculatory flow. In this study, we performed a genome-wide association study in 18,722 unrelated individuals of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium and the Blue Mountain Eye Study, to identify genetic determinants associated with variations in retinal arteriolar caliber. Retinal vascular calibers were measured on digitized retinal photographs using a standardized protocol. One variant (rs2194025 on chromosome 5q14 near the myocyte enhancer factor 2C MEF2C gene) was associated with retinal arteriolar caliber in the meta-analysis of the discovery cohorts at genome-wide significance of P-value <5×10(-8). This variant was replicated in an additional 3,939 individuals of European ancestry from the Australian Twins Study and Multi-Ethnic Study of Atherosclerosis (rs2194025, P-value = 2.11×10(-12) in combined meta-analysis of discovery and replication cohorts). In independent studies of modest sample sizes, no significant association was found between this variant and clinical outcomes including coronary artery disease, stroke, myocardial infarction or hypertension. In conclusion, we found one novel loci which underlie genetic variation in microvasculature which may be relevant to vascular disease. The relevance of these findings to clinical outcomes remains to be determined.