Search

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

AbstractChronic pain affects more than 50 million Americans. Treatments remain inadequate, in large part, because the pathophysiological mechanisms underlying the development of chronic pain remain poorly understood. Pain biomarkers could potentially identify and measure biological pathways and phenotypical expressions that are altered by pain, provide insight into biological treatment targets, and help identify at-risk patients who might benefit from early intervention. Biomarkers are used to diagnose, track, and treat other diseases, but no validated clinical biomarkers exist yet for chronic pain. To address this problem, the National Institutes of Health Common Fund launched the Acute to Chronic Pain Signatures (A2CPS) program to evaluate candidate biomarkers, develop them into biosignatures, and discover novel biomarkers for chronification of pain after surgery. This article discusses candidate biomarkers identified by A2CPS for evaluation, including genomic, proteomic, metabolomic, lipidomic, neuroimaging, psychophysical, psychological, and behavioral measures. Acute to Chronic Pain Signatures will provide the most comprehensive investigation of biomarkers for the transition to chronic postsurgical pain undertaken to date. Data and analytic resources generatedby A2CPS will be shared with the scientific community in hopes that other investigators will extract valuable insights beyond A2CPS's initial findings. This article will review the identified biomarkers and rationale for including them, the current state of the science on biomarkers of the transition from acute to chronic pain, gaps in the literature, and how A2CPS will address these gaps.

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of Asian-Ancestry (AsA) disproportionately experience more severe SLE compared to individuals of European-Ancestry (EA), including increased renal involvement and tissue damage. However, the mechanisms underlying elevated severity in the AsA population remain unclear. Here, we utilized available gene expression data and genotype data based on all non-HLA SNP associations in EA and AsA SLE patients detected using the Immunochip genotyping array. We identified 2778 ancestry-specific and 327 trans-ancestry SLE-risk polymorphisms. Genetic associations were examined using connectivity mapping and gene signatures based on predicted biological pathways and were used to interrogate gene expression datasets. SLE-associated pathways in AsA patients included elevated oxidative stress, altered metabolism and mitochondrial dysfunction, whereas SLE-associated pathways in EA patients included a robust interferon response (type I and II) related to enhanced cytosolic nucleic acid sensing and signaling. An independent dataset derived from summary genome-wide association data in an AsA cohort was interrogated and identified similar molecular pathways. Finally, gene expression data from AsA SLE patients corroborated the molecular pathways predicted by SNP associations. Identifying ancestry-related molecular pathways predicted by genetic SLE risk may help to disentangle the population differences in clinical severity that impact AsA and EA individuals with SLE.

IntroductionAfrican Americans are at increased risk for type 2 diabetes.ObjectivesThis work aimed to examine metabolomic signature of glucose homeostasis in African Americans.MethodsWe used an untargeted liquid chromatography-mass spectrometry metabolomic approach to comprehensively profile 727 plasma metabolites among 571 African Americans from the Insulin Resistance Atherosclerosis Family Study (IRAS-FS) and investigate the associations between these metabolites and both the dynamic (SI, insulin sensitivity; AIR, acute insulin response; DI, disposition index; and SG, glucose effectiveness) and basal (HOMA-IR and HOMA-B) measures of glucose homeostasis using univariate and regularized regression models. We also compared the results with our previous findings in the IRAS-FS Mexican Americans.ResultsWe confirmed increased plasma metabolite levels of branched-chain amino acids and their metabolic derivatives, 2-aminoadipate, 2-hydroxybutyrate, glutamate, arginine and its metabolic derivatives, carbohydrate metabolites, and medium- and long-chain fatty acids were associated with insulin resistance, while increased plasma metabolite levels in the glycine, serine and threonine metabolic pathway were associated with insulin sensitivity. We also observed a differential ancestral effect of glutamate on glucose homeostasis with significantly stronger effects observed in African Americans than those previously observed in Mexican Americans.ConclusionWe extended the observations that metabolites are useful biomarkers in the identification of prediabetes in individuals at risk of type 2 diabetes in African Americans. We revealed, for the first time, differential ancestral effect of certain metabolites (i.e., glutamate) on glucose homeostasis traits. Our study highlights the need for additional comprehensive metabolomic studies in well-characterized multiethnic cohorts.

BackgroundCoronary artery calcification (CAC) and carotid artery intima-media thickness (cIMT) are measures of subclinical atherosclerosis in asymptomatic individuals and strong risk factors for cardiovascular disease. Type 2 diabetes (T2D) is an independent cardiovascular disease risk factor that accelerates atherosclerosis.MethodsWe performed meta-analyses of genome-wide association studies in up to 2500 T2D individuals of European ancestry (EA) and 1590 T2D individuals of African ancestry with or without exclusion of prevalent cardiovascular disease, for CAC measured by cardiac computed tomography, and 3608 individuals of EA and 838 individuals of African ancestry with T2D for cIMT measured by ultrasonography within the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium.ResultsWe replicated 2 loci (rs9369640 and rs9349379 near PHACTR1 and rs10757278 near CDKN2B) for CAC and one locus for cIMT (rs7412 and rs445925 near APOE-APOC1) that were previously reported in the general EA populations. We identified one novel CAC locus (rs8000449 near CSNK1A1L/LINC00547/POSTN at 13q13.3) at P=2.0×10-8 in EA. No additional loci were identified with the meta-analyses of EA and African ancestry. The expression quantitative trait loci analysis with nearby expressed genes derived from arterial wall and metabolic tissues from the Genotype-Tissue Expression project pinpoints POSTN, encoding a matricellular protein involved in bone formation and bone matrix organization, as the potential candidate gene at this locus. In addition, we found significant associations (P

Decline in insulin sensitivity due to dysfunction of adipose tissue (AT) is one of the earliest pathogenic events in type 2 diabetes. We hypothesize that differential DNA methylation (DNAm) controls insulin sensitivity and obesity by modulating transcript expression in AT. Integrating AT DNAm profiles with transcript profile data measured in a cohort of 230 African Americans (AAs) from the African American Genetics of Metabolism and Expression cohort, we performed cis-expression quantitative trait methylation (cis-eQTM) analysis to identify epigenetic regulatory loci for glucometabolic trait-associated transcripts. We identified significantly associated cytosine-guanine dinucleotide regions for 82 transcripts (false discovery rate [FDR]-P < 0.05). The strongest eQTM locus was observed for the proopiomelanocortin (POMC; ρ = -0.632, P = 4.70 × 10-27) gene. Epigenome-wide association studies (EWAS) further identified 155, 46, and 168 cytosine-guanine dinucleotide regions associated (FDR-P < 0.05) with the Matsuda index, SI, and BMI, respectively. Intersection of EWAS, transcript level to trait association, and eQTM results, followed by causal inference test identified significant eQTM loci for 23 genes that were also associated with Matsuda index, SI, and/or BMI in EWAS. These associated genes include FERMT3, ITGAM, ITGAX, and POMC In summary, applying an integrative multiomics approach, our study provides evidence for DNAm-mediated regulation of gene expression at both previously identified and novel loci for many key AT transcripts influencing insulin resistance and obesity.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

ObjectivePopulations of Mexican American ancestry are at an increased risk for nonalcoholic fatty liver disease. The objective of this study was to determine whether loci in known and novel genes were associated with variation in aspartate aminotransferase (AST) (n = 3,644), alanine aminotransferase (ALT) (n = 3,595), and gamma-glutamyl transferase (GGT) (n = 1,577) levels by conducting the first genome-wide association study (GWAS) of liver enzymes, which commonly measure liver function, in individuals of Mexican American ancestry.MethodsLevels of AST, ALT, and GGT were determined by enzymatic colorimetric assays. A multi-cohort GWAS of individuals of Mexican American ancestry was performed. Single-nucleotide polymorphisms (SNP) were tested for association with liver outcomes by multivariable linear regression using an additive genetic model. Association analyses were conducted separately in each cohort, followed by a nonparametric meta-analysis.ResultsIn the PNPLA3 gene, rs4823173 (P = 3.44 × 10-10 ), rs2896019 (P = 7.29 × 10-9 ), and rs2281135 (P = 8.73 × 10-9 ) were significantly associated with AST levels. Although not genome-wide significant, these same SNPs were the top hits for ALT (P = 7.12 × 10-8 , P = 1.98 × 10-7 , and P = 1.81 × 10-7 , respectively). The strong correlation (r2  = 1.0) for these SNPs indicated a single hit in the PNPLA3 gene. No genome-wide significant associations were found for GGT.ConclusionsPNPLA3, a locus previously identified with ALT, AST, and nonalcoholic fatty liver disease in European and Japanese GWAS, is also associated with liver enzymes in populations of Mexican American ancestry.

PurposeTo describe the study protocol and baseline characteristics of the African Descent and Glaucoma Evaluation Study (ADAGES) III.DesignCross-sectional, case-control study.ParticipantsThree thousand two hundred sixty-six glaucoma patients and control participants without glaucoma of African or European descent were recruited from 5 study centers in different regions of the United States.MethodsIndividuals of African descent (AD) and European descent (ED) with primary open-angle glaucoma (POAG) and control participants completed a detailed demographic and medical history interview. Standardized height, weight, and blood pressure measurements were obtained. Saliva and blood samples to provide serum, plasma, DNA, and RNA were collected for standardized processing. Visual fields, stereoscopic disc photographs, and details of the ophthalmic examination were obtained and transferred to the University of California, San Diego, Data Coordinating Center for standardized processing and quality review.Main outcome measuresParticipant gender, age, race, body mass index, blood pressure, history of smoking and alcohol use in POAG patients and control participants were described. Ophthalmic measures included intraocular pressure, visual field mean deviation, central corneal thickness, glaucoma medication use, or past glaucoma surgery. Ocular conditions, including diabetic retinopathy, age-related macular degeneration, and past cataract surgery, were recorded.ResultsThe 3266 ADAGES III study participants in this report include 2146 AD POAG patients, 695 ED POAG patients, 198 AD control participants, and 227 ED control participants. The AD POAG patients and control participants were significantly younger (both, 67.4 years) than ED POAG patients and control participants (73.4 and 70.2 years, respectively). After adjusting for age, AD POAG patients had different phenotypic characteristics compared with ED POAG patients, including higher intraocular pressure, worse visual acuity and visual field mean deviation, and thinner corneas (all P < 0.001). Family history of glaucoma did not differ between AD and ED POAG patients.ConclusionsWith its large sample size, extensive specimen collection, and deep phenotyping of AD and ED glaucoma patients and control participants from different regions in the United States, the ADAGES III genomics study will address gaps in our knowledge of the genetics of POAG in this high-risk population.

PurposeTo find genetic contributions to glaucoma in African Americans.DesignCross-sectional, case-control study.ParticipantsOne thousand eight hundred seventy-five primary open-angle glaucoma (POAG) patients and 1709 controls, self-identified as being of African descent (AD), from the African Descent and Glaucoma Evaluation Study (ADAGES) III and Wake Forest School of Medicine.MethodsMegaChip genotypes were imputed to Thousand Genomes data. Association of single nucleotide polymorphisms (SNPs) with POAG and advanced POAG was tested by linear mixed model correcting for relatedness and population stratification. Genetic risk scores were tested by receiver operator characteristic curves (ROC-AUCs).Main outcome measuresPrimary open-angle glaucoma defined by visual field loss without other nonocular conditions (n = 1875). Advanced POAG was defined by age-based mean deviation of visual field (n = 946).ResultsEighteen million two hundred eighty-one thousand nine hundred twenty SNPs met imputation quality of r2 > 0.7 and minor allele frequency > 0.005. Association of a novel locus, EN04, was observed for advanced POAG (rs185815146 β, 0.36; standard error, 0.065; P < 3×10-8). For POAG, an AD signal was observed at the 9p21 European descent (ED) POAG signal (rs79721419; P < 6.5×10-5) independent of the previously observed 9p21 ED signal (rs2383204; P < 2.3×10-5) by conditional analyses. An association with POAG in FNDC3B (rs111698934; P < 3.9×10-5) was observed, not in linkage disequilibrium (LD) with the previously reported ED SNP. Additional previously identified loci associated with POAG in persons of AD were: 8q22, AFAP1, and TMC01. An AUC of 0.62 was observed with an unweighted genetic risk score comprising 11 SNPs in candidate genes. Two additional risk scores were studied by using a penalized matrix decomposition with cross-validation; risk scores of 50 and 400 SNPs were identified with ROC of AUC = 0.74 and AUC = 0.94, respectively.ConclusionsA novel association with advanced POAG in the EN04 locus was identified putatively in persons of AD. In addition to this finding, this genome-wide association study in POAG patients of AD contributes to POAG genetics by identification of novel signals in prior loci (9p21), as well as advancing the fine mapping of regions because of shorter average LD (FNDC3B). Although not useful without confirmation and clinical trials, the use of genetic risk scores demonstrated that considerable AD-specific genetic information remains in these data.

Carotid artery intima media thickness (cIMT) and carotid plaque are measures of subclinical atherosclerosis associated with ischemic stroke and coronary heart disease (CHD). Here, we undertake meta-analyses of genome-wide association studies (GWAS) in 71,128 individuals for cIMT, and 48,434 individuals for carotid plaque traits. We identify eight novel susceptibility loci for cIMT, one independent association at the previously-identified PINX1 locus, and one novel locus for carotid plaque. Colocalization analysis with nearby vascular expression quantitative loci (cis-eQTLs) derived from arterial wall and metabolic tissues obtained from patients with CHD identifies candidate genes at two potentially additional loci, ADAMTS9 and LOXL4. LD score regression reveals significant genetic correlations between cIMT and plaque traits, and both cIMT and plaque with CHD, any stroke subtype and ischemic stroke. Our study provides insights into genes and tissue-specific regulatory mechanisms linking atherosclerosis both to its functional genomic origins and its clinical consequences in humans.

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

ObjectiveThis study aimed to explore the genetic mechanisms of regional fat deposition, which is a strong risk factor for metabolic diseases beyond total adiposity.MethodsA genome-wide association study of 7,757,139 single-nucleotide polymorphisms (SNPs) in 983 Mexican Americans (nmale  = 403; nfemale  = 580) from the Insulin Resistance Atherosclerosis Family Study was performed. Association analyses were performed with and without sex stratification for subcutaneous adipose tissue, visceral adipose tissue (VAT), and visceral-subcutaneous ratio (VSR) obtained from computed tomography.ResultsThe strongest signal identified was SNP rs2185405 (minor allele frequencies [MAF] = 40%; PVAT  = 1.98 × 10-8 ) with VAT. It is an intronic variant of the GLIS family zinc finger 3 gene (GLIS3). In addition, SNP rs12657394 (MAF = 19%) was associated with VAT in males (Pmale  = 2.39×10-8 ; Pfemale  = 2.5 × 10-3 ). It is located intronically in the serum response factor binding protein 1 gene (SRFBP1). On average, male carriers of the variant had 24.6 cm2 increased VAT compared with noncarriers. Subsequently, genome-wide SNP-sex interaction analysis was performed. SNP rs10913233 (MAF = 14%; Pint  = 3.07 × 10-8 ) in PAPPA2 and rs10923724 (MAF = 38%; Pint  = 2.89 × 10-8 ) upstream of TBX15 were strongly associated with the interaction effect for VSR.ConclusionsSix loci were identified with genome-wide significant associations with fat deposition and interactive effects. These results provided genetic evidence for a differential basis of fat deposition between genders.

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P

Type 2 diabetes (T2D) affects Latinos at twice the rate seen in populations of European descent. We recently identified a risk haplotype spanning SLC16A11 that explains ∼20% of the increased T2D prevalence in Mexico. Here, through genetic fine-mapping, we define a set of tightly linked variants likely to contain the causal allele(s). We show that variants on the T2D-associated haplotype have two distinct effects: (1) decreasing SLC16A11 expression in liver and (2) disrupting a key interaction with basigin, thereby reducing cell-surface localization. Both independent mechanisms reduce SLC16A11 function and suggest SLC16A11 is the causal gene at this locus. To gain insight into how SLC16A11 disruption impacts T2D risk, we demonstrate that SLC16A11 is a proton-coupled monocarboxylate transporter and that genetic perturbation of SLC16A11 induces changes in fatty acid and lipid metabolism that are associated with increased T2D risk. Our findings suggest that increasing SLC16A11 function could be therapeutically beneficial for T2D. VIDEO ABSTRACT.

Insertions and deletions (INDELs) represent a significant fraction of interindividual variation in the human genome yet their contribution to phenotypes is poorly understood. To confirm the quality of imputed INDELs and investigate their roles in mediating cardiometabolic phenotypes, genome-wide association and linkage analyses were performed for 15 phenotypes with 1,273,952 imputed INDELs in 1,024 Mexican-origin Americans. Imputation quality was validated using whole exome sequencing with an average kappa of 0.93 in common INDELs (minor allele frequencies [MAFs] ≥ 5%). Association analysis revealed one genome-wide significant association signal for the cholesterylester transfer protein gene (CETP) with high-density lipoprotein levels (rs36229491, P = 3.06 × 10-12 ); linkage analysis identified two peaks with logarithm of the odds (LOD) > 5 (rs60560566, LOD = 5.36 with insulin sensitivity (SI ) and rs5825825, LOD = 5.11 with adiponectin levels). Suggestive overlapping signals between linkage and association were observed: rs59849892 in the WSC domain containing 2 gene (WSCD2) was associated and nominally linked with SI (P = 1.17 × 10-7 , LOD = 1.99). This gene has been implicated in glucose metabolism in human islet cell expression studies. In addition, rs201606363 was linked and nominally associated with low-density lipoprotein (P = 4.73 × 10-4 , LOD = 3.67), apolipoprotein B (P = 1.39 × 10-3 , LOD = 4.64), and total cholesterol (P = 1.35 × 10-2 , LOD = 3.80) levels. rs201606363 is an intronic variant of the UBE2F-SCLY (where UBE2F is ubiquitin-conjugating enzyme E2F and SCLY is selenocysteine lyase) fusion gene that may regulate cholesterol through selenium metabolism. In conclusion, these results confirm the feasibility of imputing INDELs from array-based single nucleotide polymorphism (SNP) genotypes. Analysis of these variants using association and linkage replicated previously identified SNP signals and identified multiple novel INDEL signals. These results support the inclusion of INDELs into genetic studies to more fully interrogate the spectrum of genetic variation.

Family-based methods are a potentially powerful tool to identify trait-defining genetic variants in extended families, particularly when used to complement conventional association analysis. We utilized two-point linkage analysis and single variant association analysis to evaluate whole exome sequencing (WES) data from 1205 Hispanic Americans (78 families) from the Insulin Resistance Atherosclerosis Family Study. WES identified 211,612 variants above the minor allele frequency threshold of ≥0.005. These variants were tested for linkage and/or association with 50 cardiometabolic traits after quality control checks. Two-point linkage analysis yielded 10,580,600 logarithm of the odds (LOD) scores with 1148 LOD scores ≥3, 183 LOD scores ≥4, and 29 LOD scores ≥5. The maximal novel LOD score was 5.50 for rs2289043:T>C, in UNC5C with subcutaneous adipose tissue volume. Association analysis identified 13 variants attaining genome-wide significance (P T in APOA5 and triglyceride levels (P  =  3.67 × 10-10 ). Overall, there was a 5.2-fold increase in the number of informative variants detected by WES compared to exome chip analysis in this population, nearly 30% of which were novel variants relative to the Database of Single Nucleotide Polymorphisms (dbSNP) build 138. Thus, integration of results from two-point linkage and single-variant association analysis from WES data enabled identification of novel signals potentially contributing to cardiometabolic traits.

Linkage studies of complex genetic diseases have been largely replaced by genome-wide association studies, due in part to limited success in complex trait discovery. However, recent interest in rare and low-frequency variants motivates re-examination of family-based methods. In this study, we investigated the performance of two-point linkage analysis for over 1.6 million single-nucleotide polymorphisms (SNPs) combined with single variant association analysis to identify high impact variants, which are both strongly linked and associated with cardiometabolic traits in up to 1414 Hispanics from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of all 50 phenotypes yielded 83 557 000 LOD (logarithm of the odds) scores, with 9214 LOD scores ⩾3.0, 845 ⩾4.0 and 89 ⩾5.0, with a maximal LOD score of 6.49 (rs12956744 in the LAMA1 gene for tumor necrosis factor-α (TNFα) receptor 2). Twenty-seven variants were associated with P4, including variants in the NFIB gene under a linkage peak with TNFα receptor 2 levels on chromosome 9. Linkage regions of interest included a broad peak (31 Mb) on chromosome 1q with acute insulin response (max LOD=5.37). This region was previously documented with type 2 diabetes in family-based studies, providing support for the validity of these results. Overall, we have demonstrated the utility of two-point linkage and association in comprehensive genome-wide array-based SNP genotypes.

Variation in body fat distribution contributes to the metabolic sequelae of obesity. The genetic determinants of body fat distribution are poorly understood. The goal of this study was to gain new insights into the underlying genetics of body fat distribution by conducting sample-size-weighted fixed-effects genome-wide association meta-analyses in up to 9,594 women and 8,738 men of European, African, Hispanic and Chinese ancestry, with and without sex stratification, for six traits associated with ectopic fat (hereinafter referred to as ectopic-fat traits). In total, we identified seven new loci associated with ectopic-fat traits (ATXN1, UBE2E2, EBF1, RREB1, GSDMB, GRAMD3 and ENSA; P < 5 × 10-8; false discovery rate < 1%). Functional analysis of these genes showed that loss of function of either Atxn1 or Ube2e2 in primary mouse adipose progenitor cells impaired adipocyte differentiation, suggesting physiological roles for ATXN1 and UBE2E2 in adipogenesis. Future studies are necessary to further explore the mechanisms by which these genes affect adipocyte biology and how their perturbations contribute to systemic metabolic disease.

BackgroundThe presence of population structure in a sample may confound the search for important genetic loci associated with disease. Our four samples in the Family Investigation of Nephropathy and Diabetes (FIND), European Americans, Mexican Americans, African Americans, and American Indians are part of a genome- wide association study in which population structure might be particularly important. We therefore decided to study in detail one component of this, individual genetic ancestry (IGA). From SNPs present on the Affymetrix 6.0 Human SNP array, we identified 3 sets of ancestry informative markers (AIMs), each maximized for the information in one the three contrasts among ancestral populations: Europeans (HAPMAP, CEU), Africans (HAPMAP, YRI and LWK), and Native Americans (full heritage Pima Indians). We estimate IGA and present an algorithm for their standard errors, compare IGA to principal components, emphasize the importance of balancing information in the ancestry informative markers (AIMs), and test the association of IGA with diabetic nephropathy in the combined sample.ResultsA fixed parental allele maximum likelihood algorithm was applied to the FIND to estimate IGA in four samples: 869 American Indians; 1385 African Americans; 1451 Mexican Americans; and 826 European Americans. When the information in the AIMs is unbalanced, the estimates are incorrect with large error. Individual genetic admixture is highly correlated with principle components for capturing population structure. It takes ~700 SNPs to reduce the average standard error of individual admixture below 0.01. When the samples are combined, the resulting population structure creates associations between IGA and diabetic nephropathy.ConclusionsThe identified set of AIMs, which include American Indian parental allele frequencies, may be particularly useful for estimating genetic admixture in populations from the Americas. Failure to balance information in maximum likelihood, poly-ancestry models creates biased estimates of individual admixture with large error. This also occurs when estimating IGA using the Bayesian clustering method as implemented in the program STRUCTURE. Odds ratios for the associations of IGA with disease are consistent with what is known about the incidence and prevalence of diabetic nephropathy in these populations.

Type 2 diabetes (T2D) is a heterogeneous disorder with contributions from peripheral insulin resistance and β-cell dysfunction. For minimization of phenotypic heterogeneity, quantitative intermediate phenotypes characterizing basal glucose homeostasis (insulin resistance and HOMA of insulin resistance [HOMAIR] and of β-cell function [HOMAB]) have shown promise in relatively large samples. We investigated the utility of dynamic measures of glucose homeostasis (insulin sensitivity [SI] and acute insulin response [AIRg]) evaluating T2D-susceptibility variants (n = 57) in Hispanic Americans from the GUARDIAN Consortium (n = 2,560). Basal and dynamic measures were genetically correlated (HOMAB-AIRg: ρG = 0.28-0.73; HOMAIR-SI: ρG = -0.73 to -0.83) with increased heritability for the dynamic measure AIRg Significant association of variants with dynamic measures (P < 8.77 × 10(-4)) was observed. A pattern of superior performance of AIRg was observed for well-established loci including MTNR1B (P = 9.46 × 10(-12)), KCNQ1 (P = 1.35 × 10(-4)), and TCF7L2 (P = 5.10 × 10(-4)) with study-wise statistical significance. Notably, significant association of MTNR1B with AIRg (P < 1.38 × 10(-9)) was observed in a population one-fourteenth the size of the initial discovery cohort. These observations suggest that basal and dynamic measures provide different views and levels of sensitivity to discrete elements of glucose homeostasis. Although more costly to obtain, dynamic measures yield significant results that could be considered physiologically "closer" to causal pathways and provide insight into the discrete mechanisms of action.

Knowledge of the genetic basis of the type 2 diabetes (T2D)-related quantitative traits fasting glucose (FG) and insulin (FI) in African ancestry (AA) individuals has been limited. In non-diabetic subjects of AA (n = 20,209) and European ancestry (EA; n = 57,292), we performed trans-ethnic (AA+EA) fine-mapping of 54 established EA FG or FI loci with detailed functional annotation, assessed their relevance in AA individuals, and sought previously undescribed loci through trans-ethnic (AA+EA) meta-analysis. We narrowed credible sets of variants driving association signals for 22/54 EA-associated loci; 18/22 credible sets overlapped with active islet-specific enhancers or transcription factor (TF) binding sites, and 21/22 contained at least one TF motif. Of the 54 EA-associated loci, 23 were shared between EA and AA. Replication with an additional 10,096 AA individuals identified two previously undescribed FI loci, chrX FAM133A (rs213676) and chr5 PELO (rs6450057). Trans-ethnic analyses with regulatory annotation illuminate the genetic architecture of glycemic traits and suggest gene regulation as a target to advance precision medicine for T2D. Our approach to utilize state-of-the-art functional annotation and implement trans-ethnic association analysis for discovery and fine-mapping offers a framework for further follow-up and characterization of GWAS signals of complex trait loci.

Leptin is abnormally elevated in the plasma of patients with systemic lupus erythematosus (SLE), where it is thought to promote and/or sustain pro-inflammatory responses. Whether this association could reflect an increased genetic susceptibility to develop SLE is not known, and studies of genetic associations with leptin-related polymorphisms in SLE patients have been so far inconclusive. Here we genotyped DNA samples from 15,706 SLE patients and healthy matched controls from four different ancestral groups, to correlate polymorphisms of genes of the leptin pathway to risk for SLE. It was found that although several SNPs showed weak associations, those associations did not remain significant after correction for multiple testing. These data do not support associations between defined leptin-related polymorphisms and increased susceptibility to develop SLE.

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

Background and purposeThe majority of genome-wide association studies (GWAS) of stroke have focused on European-ancestry populations; however, none has been conducted in African Americans, despite the disproportionately high burden of stroke in this population. The Consortium of Minority Population Genome-Wide Association Studies of Stroke (COMPASS) was established to identify stroke susceptibility loci in minority populations.MethodsUsing METAL, we conducted meta-analyses of GWAS in 14 746 African Americans (1365 ischemic and 1592 total stroke cases) from COMPASS, and tested genetic variants with P

Genetic variants at chromosomal region 11q23.3, near the gene ETS1, have been associated with systemic lupus erythematosus (SLE), or lupus, in independent cohorts of Asian ancestry. Several recent studies have implicated ETS1 as a critical driver of immune cell function and differentiation, and mice deficient in ETS1 develop an SLE-like autoimmunity. We performed a fine-mapping study of 14,551 subjects from multi-ancestral cohorts by starting with genotyped variants and imputing to all common variants spanning ETS1. By constructing genetic models via frequentist and Bayesian association methods, we identified 16 variants that are statistically likely to be causal. We functionally assessed each of these variants on the basis of their likelihood of affecting transcription factor binding, miRNA binding, or chromatin state. Of the four variants that we experimentally examined, only rs6590330 differentially binds lysate from B cells. Using mass spectrometry, we found more binding of the transcription factor signal transducer and activator of transcription 1 (STAT1) to DNA near the risk allele of rs6590330 than near the non-risk allele. Immunoblot analysis and chromatin immunoprecipitation of pSTAT1 in B cells heterozygous for rs6590330 confirmed that the risk allele increased binding to the active form of STAT1. Analysis with expression quantitative trait loci indicated that the risk allele of rs6590330 is associated with decreased ETS1 expression in Han Chinese, but not other ancestral cohorts. We propose a model in which the risk allele of rs6590330 is associated with decreased ETS1 expression and increases SLE risk by enhancing the binding of pSTAT1.

Insulin sensitivity, insulin secretion, insulin clearance, and glucose effectiveness exhibit strong genetic components, although few studies have examined their genetic architecture or influence on type 2 diabetes (T2D) risk. We hypothesized that loci affecting variation in these quantitative traits influence T2D. We completed a multicohort genome-wide association study to search for loci influencing T2D-related quantitative traits in 4,176 Mexican Americans. Quantitative traits were measured by the frequently sampled intravenous glucose tolerance test (four cohorts) or euglycemic clamp (three cohorts), and random-effects models were used to test the association between loci and quantitative traits, adjusting for age, sex, and admixture proportions (Discovery). Analysis revealed a significant (P < 5.00 × 10(-8)) association at 11q14.3 (MTNR1B) with acute insulin response. Loci with P < 0.0001 among the quantitative traits were examined for translation to T2D risk in 6,463 T2D case and 9,232 control subjects of Mexican ancestry (Translation). Nonparametric meta-analysis of the Discovery and Translation cohorts identified significant associations at 6p24 (SLC35B3/TFAP2A) with glucose effectiveness/T2D, 11p15 (KCNQ1) with disposition index/T2D, and 6p22 (CDKAL1) and 11q14 (MTNR1B) with acute insulin response/T2D. These results suggest that T2D and insulin secretion and sensitivity have both shared and distinct genetic factors, potentially delineating genomic components of these quantitative traits that drive the risk for T2D.

We previously identified a low-frequency (1.1 %) coding variant (G45R; rs200573126) in the adiponectin gene (ADIPOQ) which was the basis for a multipoint microsatellite linkage signal (LOD = 8.2) for plasma adiponectin levels in Hispanic families. We have empirically evaluated the ability of data from targeted common variants, exome chip genotyping, and genome-wide association study data to detect linkage and association to adiponectin protein levels at this locus. Simple two-point linkage and association analyses were performed in 88 Hispanic families (1,150 individuals) using 10,958 SNPs on chromosome 3. Approaches were compared for their ability to map the functional variant, G45R, which was strongly linked (two-point LOD = 20.98) and powerfully associated (p value = 8.1 × 10(-50)). Over 450 SNPs within a broad 61 Mb interval around rs200573126 showed nominal evidence of linkage (LOD > 3) but only four other SNPs in this region were associated with p values

Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5-TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5-TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10(-49); OR = 1.38-1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3 (P-valuesEU = 10(-27)-10(-32), OR = 1.7-1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credible set of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5-TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5-TNPO3.

Obesity is growing epidemic affecting 35% of adults in the United States. Previous genome-wide association studies (GWAS) have identified numerous loci associated with obesity. However, the majority of studies have been completed in Caucasians focusing on total body measures of adiposity. Here we report the results from genome-wide and exome chip association studies focusing on total body measures of adiposity including body mass index (BMI), percent body fat (PBF) and measures of fat deposition including waist circumference (WAIST), waist-hip ratio (WHR), subcutaneous adipose tissue (SAT), and visceral adipose tissue (VAT) in Hispanic Americans (nmax = 1263) from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Five SNPs from two novel loci attained genome-wide significance (P

Apolipoprotein L1 gene (APOL1) G1 and G2 coding variants are strongly associated with chronic kidney disease (CKD) in African Americans (AAs). Here APOL1 association was tested with baseline estimated glomerular filtration rate (eGFR), urine albumin:creatinine ratio (UACR), and prevalent cardiovascular disease (CVD) in 2571 AAs from the Systolic Blood Pressure Intervention Trial (SPRINT), a trial assessing effects of systolic blood pressure reduction on renal and CVD outcomes. Logistic regression models that adjusted for potentially important confounders tested for association between APOL1 risk variants and baseline clinical CVD (myocardial infarction, coronary, or carotid artery revascularization) and CKD (eGFR under 60 ml/min per 1.73 m(2) and/or UACR over 30 mg/g). AA SPRINT participants were 45.3% female with a mean (median) age of 64.3 (63) years, mean arterial pressure 100.7 (100) mm Hg, eGFR 76.3 (77.1) ml/min per 1.73 m(2), and UACR 49.9 (9.2) mg/g, and 8.2% had clinical CVD. APOL1 (recessive inheritance) was positively associated with CKD (odds ratio 1.37, 95% confidence interval 1.08-1.73) and log UACR estimated slope (β) 0.33) and negatively associated with eGFR (β -3.58), all significant. APOL1 risk variants were not significantly associated with prevalent CVD (1.02, 0.82-1.27). Thus, SPRINT data show that APOL1 risk variants are associated with mild CKD but not with prevalent CVD in AAs with a UACR under 1000 mg/g.

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15 × 10(-94)

Linkage analysis of complex traits has had limited success in identifying trait-influencing loci. Recently, coding variants have been implicated as the basis for some biomedical associations. We tested whether coding variants are the basis for linkage peaks of complex traits in 42 African-American (n = 596) and 90 Hispanic (n = 1,414) families in the Insulin Resistance Atherosclerosis Family Study (IRASFS) using Illumina HumanExome Beadchips. A total of 92,157 variants in African Americans (34%) and 81,559 (31%) in Hispanics were polymorphic and tested using two-point linkage and association analyses with 37 cardiometabolic phenotypes. In African Americans 77 LOD scores greater than 3 were observed. The highest LOD score was 4.91 with the APOE SNP rs7412 (MAF = 0.13) with plasma apolipoprotein B (ApoB). This SNP was associated with ApoB (P-value = 4 × 10(-19)) and accounted for 16.2% of the variance in African Americans. In Hispanic families, 104 LOD scores were greater than 3. The strongest evidence of linkage (LOD = 4.29) was with rs5882 (MAF = 0.46) in CETP with HDL. CETP variants were strongly associated with HDL (0.00049 < P-value

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

ObjectiveThe GUARDIAN (Genetics UndeRlying DIAbetes in HispaNics) consortium is described, along with heritability estimates and genetic and environmental correlations of insulin sensitivity and metabolic clearance rate of insulin (MCRI).MethodsGUARDIAN is comprised of seven cohorts, consisting of 4,336 Mexican-American individuals in 1,346 pedigrees. Insulin sensitivity (SI ), MCRI, and acute insulin response (AIRg) were measured by frequently sampled intravenous glucose tolerance test in four cohorts. Insulin sensitivity (M, M/I) and MCRI were measured by hyperinsulinemic-euglycemic clamp in three cohorts. Heritability and genetic and environmental correlations were estimated within the family cohorts (totaling 3,925 individuals) using variance components.ResultsAcross studies, age, and gender-adjusted heritability of insulin sensitivity (SI , M, M/I) ranged from 0.23 to 0.48 and of MCRI from 0.35 to 0.73. The ranges for the genetic correlations were 0.91 to 0.93 between SI and MCRI; and -0.57 to -0.59 for AIRg and MCRI (all P

ObjectiveLupus nephritis (LN) is a severe manifestation of systemic lupus erythematosus (SLE) that exhibits familial aggregation and may progress to end-stage renal disease (ESRD). LN is more prevalent among African Americans than among European Americans. This study was undertaken to investigate the hypothesis that the apolipoprotein L1 gene (APOL1) nephropathy risk alleles G1/G2, common in African Americans and rare in European Americans, contribute to the ethnic disparity in risk.MethodsAPOL1 G1 and G2 nephropathy alleles were genotyped in 855 African American SLE patients with LN-ESRD (cases) and 534 African American SLE patients without nephropathy (controls) and tested for association under a recessive genetic model, by logistic regression.ResultsNinety percent of the SLE patients were female. The mean ± SD age at SLE diagnosis was significantly lower in LN-ESRD cases than in SLE non-nephropathy controls (27.3 ± 10.9 years versus 39.5 ± 12.2 years). The mean ± SD time from SLE diagnosis to development of LN-ESRD in cases was 7.3 ± 7.2 years. The G1/G2 risk alleles were strongly associated with SLE-ESRD, with 25% of cases and 12% of controls having 2 nephropathy alleles (odds ratio [OR] 2.57, recessive model P = 1.49 × 10(-9)), and after adjustment for age, sex, and ancestry admixture (OR 2.72, P = 6.23 × 10(-6)). The age-, sex-, and admixture-adjusted population attributable risk for ESRD among patients with G1/G2 polymorphisms was 0.26, compared to 0.003 among European American patients. The mean time from SLE diagnosis to ESRD development was ∼2 years earlier among individuals with APOL1 risk genotypes (P = 0.01).ConclusionAPOL1 G1/G2 alleles strongly impact the risk of LN-ESRD in African Americans, as well as the time to progression to ESRD. The high frequency of these alleles in African Americans with near absence in European Americans explains an important proportion of the increased risk of LN-ESRD in African Americans.

Genome-wide association studies (GWAS) have identified 36 loci associated with body mass index (BMI), predominantly in populations of European ancestry. We conducted a meta-analysis to examine the association of >3.2 million SNPs with BMI in 39,144 men and women of African ancestry and followed up the most significant associations in an additional 32,268 individuals of African ancestry. We identified one new locus at 5q33 (GALNT10, rs7708584, P = 3.4 × 10(-11)) and another at 7p15 when we included data from the GIANT consortium (MIR148A-NFE2L3, rs10261878, P = 1.2 × 10(-10)). We also found suggestive evidence of an association at a third locus at 6q16 in the African-ancestry sample (KLHL32, rs974417, P = 6.9 × 10(-8)). Thirty-two of the 36 previously established BMI variants showed directionally consistent effect estimates in our GWAS (binomial P = 9.7 × 10(-7)), five of which reached genome-wide significance. These findings provide strong support for shared BMI loci across populations, as well as for the utility of studying ancestrally diverse populations.

ObjectiveLittle is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate-related genes biological candidates for disease susceptibility. We analyzed variation in reactive intermediate genes for association with SLE in 2 populations with African ancestry.MethodsA total of 244 single-nucleotide polymorphisms (SNP) from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls). Single-marker, haplotype, and 2-locus interaction tests were computed for these populations.ResultsThe glutathione reductase gene GSR (rs2253409; p = 0.0014, OR 1.26, 95% CI 1.09-1.44) was the most significant single SNP association in AA. In the Gullah, the NADH dehydrogenase NDUFS4 (rs381575; p = 0.0065, OR 2.10, 95% CI 1.23-3.59) and NO synthase gene NOS1 (rs561712; p = 0.0072, OR 0.62, 95% CI 0.44-0.88) were most strongly associated with SLE. When both populations were analyzed together, GSR remained the most significant effect (rs2253409; p = 0.00072, OR 1.26, 95% CI 1.10-1.44). Haplotype and 2-locus interaction analyses also uncovered different loci in each population.ConclusionThese results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.

ObjectivesThe Xq28 region containing IRAK1 and MECP2 has been identified as a risk locus for systemic lupus erythematosus (SLE) in previous genetic association studies. However, due to the strong linkage disequilibrium between IRAK1 and MECP2, it remains unclear which gene is affected by the underlying causal variant(s) conferring risk of SLE.MethodsWe fine-mapped ≥136 SNPs in a ∼227 kb region on Xq28, containing IRAK1, MECP2 and seven adjacent genes (L1CAM, AVPR2, ARHGAP4, NAA10, RENBP, HCFC1 and TMEM187), for association with SLE in 15 783 case-control subjects derived from four different ancestral groups.ResultsMultiple SNPs showed strong association with SLE in European Americans, Asians and Hispanics at p

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

ObjectiveSystemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway.MethodsWe analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively.ResultsWe found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression.ConclusionOur results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

ObjectiveAmerican Indian-Europeans, Asians, and African Americans have an excess morbidity from systemic lupus erythematosus (SLE) and a higher prevalence of lupus nephritis than do Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and sociodemographic characteristics and clinical features in a large cohort of American Indian-European SLE patients.MethodsA total of 2,116 SLE patients of American Indian-European origin and 4,001 SLE patients of European descent for whom we had clinical data were included in the study. Genotyping of 253 continental ancestry-informative markers was performed on the Illumina platform. Structure and Admixture software were used to determine genetic ancestry proportions of each individual. Logistic regression was used to test the association between genetic ancestry and sociodemographic and clinical characteristics. Odds ratios (ORs) were calculated with 95% confidence intervals (95% CIs).ResultsThe average American Indian genetic ancestry of 2,116 SLE patients was 40.7%. American Indian genetic ancestry conferred increased risks of renal involvement (P < 0.0001, OR 3.50 [95% CI 2.63- 4.63]) and early age at onset (P < 0.0001). American Indian ancestry protected against photosensitivity (P < 0.0001, OR 0.58 [95% CI 0.44-0.76]), oral ulcers (P < 0.0001, OR 0.55 [95% CI 0.42-0.72]), and serositis (P < 0.0001, OR 0.56 [95% CI 0.41-0.75]) after adjustment for age, sex, and age at onset. However, age and sex had stronger effects than genetic ancestry on malar rash, discoid rash, arthritis, and neurologic involvement.ConclusionIn general, American Indian genetic ancestry correlates with lower sociodemographic status and increases the risk of developing renal involvement and SLE at an earlier age.

ObjectiveSystemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. We undertook to study the role of TRAF6 as a candidate gene for SLE, since it plays a major role in several signaling pathways that are important for immunity and organ development.MethodsFifteen single-nucleotide polymorphisms (SNPs) across TRAF6 were evaluated in 7,490 SLE patients and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios (ORs) with 95% confidence intervals (95% CIs) were calculated.ResultsEvidence of associations was detected in multiple SNPs. The best overall P values were obtained for SNPs rs5030437 and rs4755453 (P = 7.85 × 10(-5) and P = 4.73 × 10(-5) , respectively) without significant heterogeneity among populations (P = 0.67 and P = 0.50, respectively, in Q statistic). In addition, SNP rs540386, which was previously reported to be associated with rheumatoid arthritis (RA), was found to be in linkage disequilibrium with these 2 SNPs (r(2) = 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis P = 9.15 × 10(-4) , OR 0.89 [95% CI 0.83-0.95]). The presence of thrombocytopenia improved the overall results in different populations (meta-analysis P = 1.99 × 10(-6) , OR 0.57 [95% CI 0.45-0.72], for rs5030470). Finally, evidence of family-based association in 34 African American pedigrees with the presence of thrombocytopenia was detected in 1 available SNP (rs5030437) with a Z score magnitude of 2.28 (P = 0.02) under a dominant model.ConclusionOur data indicate the presence of association of TRAF6 with SLE, consistent with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non-Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10(-06)) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn's disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non-MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.

Diabetic nephropathy (DN) is a leading cause of mortality and morbidity in patients with type 1 and type 2 diabetes. The multicenter FIND consortium aims to identify genes for DN and its associated quantitative traits, e.g. the urine albumin:creatinine ratio (ACR). Herein, the results of whole-genome linkage analysis and a sparse association scan for ACR and a dichotomous DN phenotype are reported in diabetic individuals.A genomewide scan comprising more than 5,500 autosomal single nucleotide polymorphism markers (average spacing of 0.6 cM) was performed on 1,235 nuclear and extended pedigrees (3,972 diabetic participants) ascertained for DN from African-American (AA), American-Indian (AI), European-American (EA) and Mexican-American (MA) populations.Strong evidence for linkage to DN was detected on chromosome 6p (p = 8.0 × 10(-5), LOD = 3.09) in EA families as well as suggestive evidence for linkage to chromosome 7p in AI families. Regions on chromosomes 3p in AA, 7q in EA, 16q in AA and 22q in MA displayed suggestive evidence of linkage for urine ACR. The linkage peak on chromosome 22q overlaps the MYH9/APOL1 gene region, previously implicated in AA diabetic and nondiabetic nephropathies.These results strengthen the evidence for previously identified genomic regions and implicate several novel loci potentially involved in the pathogenesis of DN.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti-dsDNA autoantibody production, a SLE-related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti-dsDNA autoantibody positive (anti-dsDNA +, n = 811) and anti-dsDNA autoantibody negative (anti-dsDNA -, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti-dsDNA + SLE. Far fewer and weaker associations were observed for anti-dsDNA - SLE. For example, rs7574865 in STAT4 had an OR for anti-dsDNA + SLE of 1.77 (95% CI 1.57-1.99, p = 2.0E-20) compared to an OR for anti-dsDNA - SLE of 1.26 (95% CI 1.12-1.41, p = 2.4E-04), with p(heterogeneity)