Search

Inbreeding depression has been observed in many different species, but in humans a systematic analysis has been difficult so far. Here, analysing more than 1.3 million individuals, the authors show that a genomic inbreeding coefficient (FROH) is associated with disadvantageous outcomes in 32 out of 100 traits tested.

Relative to European Americans, African Americans have lower 25-hydroxyvitamin D (25OHD) and vitamin D binding protein (VDBP) concentrations, higher 1,25-dihydroxyvitamin D (1,25(OH)2D3) concentrations and bone mineral density (BMD), and paradoxically reduced burdens of calcified atherosclerotic plaque (subclinical atherosclerosis). To identify genetic factors contributing to vitamin D and BMD measures, association analysis of >14M variants was conducted in a maximum of 697 African American-Diabetes Heart Study participants with type 2 diabetes (T2D). The most significant association signals were detected for VDBP on chromosome 4; variants rs7041 (β = 0.44, SE = 0.019, P = 9.4x10-86) and rs4588 (β = 0.17, SE = 0.021, P = 3.5x10-08) in the group-specific component (vitamin D binding protein) gene (GC). These variants were found to be independently associated. In addition, rs7041 was also associated with bioavailable vitamin D (BAVD; β = 0.16, SE = 0.02, P = 3.3x10-19). Six rare variants were significantly associated with 25OHD, including a non-synonymous variant in HSPG2 (rs116788687; β = -1.07, SE = 0.17, P = 2.2x10-10) and an intronic variant in TNIK (rs143555701; β = -1.01, SE = 0.18, P = 9.0x10-10), both biologically related to bone development. Variants associated with 25OHD failed to replicate in African Americans from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of vitamin D metabolism and bone mineral density phenotypes in an African American population enriched for T2D could provide insight into ethnic specific differences in vitamin D metabolism and bone mineral density.

Abstract Whole genome sequences (WGS) enable discovery of rare variants which may contribute to missing heritability of coronary artery disease (CAD). To measure their contribution, we apply the GREML-LDMS-I approach to WGS of 4949 cases and 17,494 controls of European ancestry from the NHLBI TOPMed program. We estimate CAD heritability at 34.3% assuming a prevalence of 8.2%. Ultra-rare (minor allele frequency ≤ 0.1%) variants with low linkage disequilibrium (LD) score contribute ~50% of the heritability. We also investigate CAD heritability enrichment using a diverse set of functional annotations: i) constraint; ii) predicted protein-altering impact; iii) cis-regulatory elements from a cell-specific chromatin atlas of the human coronary; and iv) annotation principal components representing a wide range of functional processes. We observe marked enrichment of CAD heritability for most functional annotations. These results reveal the predominant role of ultra-rare variants in low LD on the heritability of CAD. Moreover, they highlight several functional processes including cell type-specific regulatory mechanisms as key drivers of CAD genetic risk.

Heavy alcohol consumption is an established risk factor for hypertension; the mechanism by which alcohol consumption impact blood pressure (BP) regulation remains unknown. We hypothesized that a genome-wide association study accounting for gene-alcohol consumption interaction for BP might identify additional BP loci and contribute to the understanding of alcohol-related BP regulation. We conducted a large two-stage investigation incorporating joint testing of main genetic effects and single nucleotide variant (SNV)-alcohol consumption interactions. In Stage 1, genome-wide discovery meta-analyses in ≈131K individuals across several ancestry groups yielded 3,514 SNVs (245 loci) with suggestive evidence of association (P < 1.0 x 10-5). In Stage 2, these SNVs were tested for independent external replication in ≈440K individuals across multiple ancestries. We identified and replicated (at Bonferroni correction threshold) five novel BP loci (380 SNVs in 21 genes) and 49 previously reported BP loci (2,159 SNVs in 109 genes) in European ancestry, and in multi-ancestry meta-analyses (P < 5.0 x 10-8). For African ancestry samples, we detected 18 potentially novel BP loci (P < 5.0 x 10-8) in Stage 1 that warrant further replication. Additionally, correlated meta-analysis identified eight novel BP loci (11 genes). Several genes in these loci (e.g., PINX1, GATA4, BLK, FTO and GABBR2) have been previously reported to be associated with alcohol consumption. These findings provide insights into the role of alcohol consumption in the genetic architecture of hypertension.

Changes in select adipose tissue volumes may differentially impact bone mineral density. This study was performed to assess cross-sectional and longitudinal relationships between computed tomography-determined visceral (VAT), subcutaneous (SAT), inter-muscular (IMAT), and pericardial adipose tissue (PAT) volumes with respective changes in thoracic vertebral and lumbar vertebral volumetric trabecular bone mineral density (vBMD) in African Americans with type 2 diabetes. Generalized linear models were fitted to test relationships between baseline and change in adipose volumes with change in vBMD in 300 African American-Diabetes Heart Study participants; adjustment was performed for age, sex, diabetes duration, study interval, smoking, hypertension, BMI, kidney function, and medications. Participants were 50% female with mean ± SD age 55.1±9.0 years, diabetes duration 10.2±7.2 years, and BMI 34.7±7.7 kg/m2. Over 5.3 ± 1.4 years, mean vBMD decreased in thoracic/lumbar spine, while mean adipose tissue volumes increased in SAT, IMAT, and PAT, but not VAT depots. In fully-adjusted models, changes in lumbar and thoracic vBMD were positively associated with change in SAT (β[SE] 0.045[0.011], p

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHRadjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHRadjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHRadjBMI and eight previously established loci at P < 5×10-8: seven for BMI, and one for WHRadjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHRadjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHRadjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (

BACKGROUND:Apolipoprotein L1 gene (APOL1) G1 and G2 renal-risk variants, common in populations with recent African ancestry, are strongly associated with non-diabetic nephropathy, end-stage kidney disease, and shorter allograft survival in deceased-donor kidneys (autosomal recessive inheritance). Circulating APOL1 protein is synthesized primarily in the liver and hydrodynamic gene delivery of APOL1 G1 and G2 risk variants has caused hepatic necrosis in a murine model. METHODS:To evaluate the impact of these variants in liver transplantation, this multicenter study investigated the association of APOL1 G1 and G2 alleles in deceased African American liver donors with allograft survival. Transplant recipients were followed for liver allograft survival using data from the Scientific Registry of Transplant Recipients. RESULTS:Of the 639 liver donors evaluated, 247 had no APOL1 risk allele, 300 had 1 risk allele, and 92 had 2 risk alleles. Graft failure assessed at 15 days, 6 months, 1 year and total was not significantly associated with donor APOL1 genotype (p-values = 0.25, 0.19, 0.67 and 0.89, respectively). CONCLUSIONS:In contrast to kidney transplantation, deceased-donor APOL1 G1 and G2 risk variants do not significantly impact outcomes in liver transplantation.

Type 2 diabetes (T2D) is the result of metabolic defects in insulin secretion and insulin sensitivity, yet most T2D loci identified to date influence insulin secretion. We hypothesized that T2D loci, particularly those affecting insulin sensitivity, can be identified through interaction with insulin secretion loci. To test this hypothesis, single nucleotide polymorphisms (SNPs) associated with acute insulin response to glucose (AIRg), a dynamic measure of first-phase insulin secretion, were identified in African Americans from the Insulin Resistance Atherosclerosis Family Study (IRASFS; n = 492 subjects). These SNPs were tested for interaction, individually and jointly as a genetic risk score (GRS), using genome-wide association study (GWAS) data from five cohorts (ARIC, CARDIA, JHS, MESA, WFSM; n = 2,725 cases, 4,167 controls) with T2D as the outcome. In single variant analyses, suggestively significant (Pinteraction

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.

Abstract Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

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)

Type 2 diabetes (T2D)-associated end-stage kidney disease (ESKD) is a complex disorder resulting from the combined influence of genetic and environmental factors. This study contains a comprehensive genetic analysis of putative nephropathy loci in 965 African American (AA) cases with T2D-ESKD and 1029 AA population-based controls extending prior findings. Analysis was based on 4,341 directly genotyped and imputed single nucleotide polymorphisms (SNPs) in 22 nephropathy candidate genes. After admixture adjustment and correction for multiple comparisons, 37 SNPs across eight loci were significantly associated (1.6E-05

Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is a negative regulator of T-cell activation associated with several autoimmune diseases, including systemic lupus erythematosus (SLE). Missense rs2476601 is associated with SLE in individuals with European ancestry. Since the rs2476601 risk allele frequency differs dramatically across ethnicities, we assessed robustness of PTPN22 association with SLE and its clinical sub-phenotypes across four ethnically diverse populations. Ten SNPs were genotyped in 8220 SLE cases and 7369 controls from in European-Americans (EA), African-Americans (AA), Asians (AS), and Hispanics (HS). We performed imputation-based association followed by conditional analysis to identify independent associations. Significantly associated SNPs were tested for association with SLE clinical sub-phenotypes, including autoantibody profiles. Multiple testing was accounted for by using false discovery rate. We successfully imputed and tested allelic association for 107 SNPs within the PTPN22 region and detected evidence of ethnic-specific associations from EA and HS. In EA, the strongest association was at rs2476601 (P = 4.7 × 10(-9), OR = 1.40 (95% CI = 1.25-1.56)). Independent association with rs1217414 was also observed in EA, and both SNPs are correlated with increased European ancestry. For HS imputed intronic SNP, rs3765598, predicted to be a cis-eQTL, was associated (P = 0.007, OR = 0.79 and 95% CI = 0.67-0.94). No significant associations were observed in AA or AS. Case-only analysis using lupus-related clinical criteria revealed differences between EA SLE patients positive for moderate to high titers of IgG anti-cardiolipin (aCL IgG >20) versus negative aCL IgG at rs2476601 (P = 0.012, OR = 1.65). Association was reinforced when these cases were compared to controls (P = 2.7 × 10(-5), OR = 2.11). Our results validate that rs2476601 is the most significantly associated SNP in individuals with European ancestry. Additionally, rs1217414 and rs3765598 may be associated with SLE. Further studies are required to confirm the involvement of rs2476601 with aCL IgG.

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22-24 (LOD=6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ~1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [P(meta) = 5.20×10(-14); odds ratio, 95% confidence interval = 0.82 (0.78-0.87)], and two missense variants, rs1990760 (Ala946Thr) [P(meta) = 3.08×10(-7); 0.88 (0.84-0.93)] and rs10930046 (Arg460His) [P(dom) = 1.16×10(-8); 0.70 (0.62-0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

We previously reported that the G allele of rs3853839 at 3'untranslated region (UTR) of Toll-like receptor 7 (TLR7) was associated with elevated transcript expression and increased risk for systemic lupus erythematosus (SLE) in 9,274 Eastern Asians [P = 6.5×10(-10), odds ratio (OR) (95%CI) = 1.27 (1.17-1.36)]. Here, we conducted trans-ancestral fine-mapping in 13,339 subjects including European Americans, African Americans, and Amerindian/Hispanics and confirmed rs3853839 as the only variant within the TLR7-TLR8 region exhibiting consistent and independent association with SLE (Pmeta = 7.5×10(-11), OR = 1.24 [1.18-1.34]). The risk G allele was associated with significantly increased levels of TLR7 mRNA and protein in peripheral blood mononuclear cells (PBMCs) and elevated luciferase activity of reporter gene in transfected cells. TLR7 3'UTR sequence bearing the non-risk C allele of rs3853839 matches a predicted binding site of microRNA-3148 (miR-3148), suggesting that this microRNA may regulate TLR7 expression. Indeed, miR-3148 levels were inversely correlated with TLR7 transcript levels in PBMCs from SLE patients and controls (R(2) = 0.255, P = 0.001). Overexpression of miR-3148 in HEK-293 cells led to significant dose-dependent decrease in luciferase activity for construct driven by TLR7 3'UTR segment bearing the C allele (P = 0.0003). Compared with the G-allele construct, the C-allele construct showed greater than two-fold reduction of luciferase activity in the presence of miR-3148. Reduced modulation by miR-3148 conferred slower degradation of the risk G-allele containing TLR7 transcripts, resulting in elevated levels of gene products. These data establish rs3853839 of TLR7 as a shared risk variant of SLE in 22,613 subjects of Asian, EA, AA, and Amerindian/Hispanic ancestries (Pmeta = 2.0×10(-19), OR = 1.25 [1.20-1.32]), which confers allelic effect on transcript turnover via differential binding to the epigenetic factor miR-3148.

Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR.Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula.The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4 × 10(-5)) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5 × 10(-4)) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5 × 10(-4)) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome.The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN.

We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P=1.71 × 10(-34) , OR=1.43[1.26-1.60]) and rs1234317-T (P=1.16 × 10(-28) , OR=1.38[1.24-1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5' region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5' risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait.

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.

Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10⁻⁸, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

Genetic case-control association studies often include data on clinical covariates, such as body mass index (BMI), smoking status, or age, that may modify the underlying genetic risk of case or control samples. For example, in type 2 diabetes, odds ratios for established variants estimated from low-BMI cases are larger than those estimated from high-BMI cases. An unanswered question is how to use this information to maximize statistical power in case-control studies that ascertain individuals on the basis of phenotype (case-control ascertainment) or phenotype and clinical covariates (case-control-covariate ascertainment). While current approaches improve power in studies with random ascertainment, they often lose power under case-control ascertainment and fail to capture available power increases under case-control-covariate ascertainment. We show that an informed conditioning approach, based on the liability threshold model with parameters informed by external epidemiological information, fully accounts for disease prevalence and non-random ascertainment of phenotype as well as covariates and provides a substantial increase in power while maintaining a properly controlled false-positive rate. Our method outperforms standard case-control association tests with or without covariates, tests of gene x covariate interaction, and previously proposed tests for dealing with covariates in ascertained data, with especially large improvements in the case of case-control-covariate ascertainment. We investigate empirical case-control studies of type 2 diabetes, prostate cancer, lung cancer, breast cancer, rheumatoid arthritis, age-related macular degeneration, and end-stage kidney disease over a total of 89,726 samples. In these datasets, informed conditioning outperforms logistic regression for 115 of the 157 known associated variants investigated (P-value = 1 × 10(-9)). The improvement varied across diseases with a 16% median increase in χ(2) test statistics and a commensurate increase in power. This suggests that applying our method to existing and future association studies of these diseases may identify novel disease loci.

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

Single nucleotide polymorphisms (SNPs) in MYH9 and APOL1 on chromosome 22 (c22) are powerfully associated with non-diabetic end-stage renal disease (ESRD) in African Americans (AAs). Many AAs diagnosed with type 2 diabetic nephropathy (T2DN) have non-diabetic kidney disease, potentially masking detection of DN genes. Therefore, genome-wide association analyses were performed using the Affymetrix SNP Array 6.0 in 966 AA with T2DN and 1,032 non-diabetic, non-nephropathy (NDNN) controls, with and without adjustment for c22 nephropathy risk variants. No associations were seen between FRMD3 SNPs and T2DN before adjusting for c22 variants. However, logistic regression analysis revealed seven FRMD3 SNPs significantly interacting with MYH9-a finding replicated in 640 additional AA T2DN cases and 683 NDNN controls. Contrasting all 1,592 T2DN cases with all 1,671 NDNN controls, FRMD3 SNPs appeared to interact with the MYH9 E1 haplotype (e.g., rs942280 interaction p-value = 9.3E⁻⁷ additive; odds ratio [OR] 0.67). FRMD3 alleles were associated with increased risk of T2DN only in subjects lacking two MYH9 E1 risk haplotypes (rs942280 OR = 1.28), not in MYH9 E1 risk allele homozygotes (rs942280 OR = 0.80; homogeneity p-value = 4.3E⁻⁴). Effects were weaker stratifying on APOL1. FRMD3 SNPS were associated with T2DN, not type 2 diabetes per se, comparing AAs with T2DN to those with diabetes lacking nephropathy. T2DN-associated FRMD3 SNPs were detectable in AAs only after accounting for MYH9, with differential effects for APOL1. These analyses reveal a role for FRMD3 in AA T2DN susceptibility and accounting for c22 nephropathy risk variants can assist in detecting DN susceptibility genes.

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Variants of genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) within the chromosome 1q32 locus linked to SLE, have been associated with multiple human diseases and may contribute to dysregulated complement activation predisposing to SLE. We assessed 60 SNPs covering the CFH-CFHRs region for association with SLE in 15,864 case-control subjects derived from four ethnic groups. Significant allelic associations with SLE were detected in European Americans (EA) and African Americans (AA), which could be attributed to an intronic CFH SNP (rs6677604, in intron 11, P(meta) = 6.6×10(-8), OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, P(meta) = 2.9×10(-7), OR = 1.17) rather than to previously identified disease-associated CFH exonic SNPs, including I62V, Y402H, A474A, and D936E. In addition, allelic association of rs6677604 with SLE was subsequently confirmed in Asians (AS). Haplotype analysis revealed that the underlying causal variant, tagged by rs6677604 and rs16840639, was localized to a ~146 kb block extending from intron 9 of CFH to downstream of CFHR1. Within this block, the deletion of CFHR3 and CFHR1 (CFHR3-1Δ), a likely causal variant measured using multiplex ligation-dependent probe amplification, was tagged by rs6677604 in EA and AS and rs16840639 in AA, respectively. Deduced from genotypic associations of tag SNPs in EA, AA, and AS, homozygous deletion of CFHR3-1Δ (P(meta) = 3.2×10(-7), OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (P(meta) = 3.5×10(-4), OR = 1.14). These results suggested that the CFHR3-1Δ deletion within the SLE-associated block, but not the previously described exonic SNPs of CFH, might contribute to the development of SLE in EA, AA, and AS, providing new insights into the role of complement regulators in the pathogenesis of SLE.

African Americans have been understudied in genome wide association studies of diabetes and related traits. In the current study, we examined the joint association of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) with fasting insulin and an index of insulin resistance (HOMA-IR) in the HyperGEN study, a family based study with proband ascertainment for hypertension. This analysis is restricted to 1,040 African Americans without diabetes. We generated allele specific CNV genotypes at 872,243 autosomal loci using Birdsuite, a freely available multi-stage program. Joint tests of association for SNPs and CNVs were performed using linear mixed models adjusting for covariates and familial relationships. Our results highlight SNPs associated with fasting insulin and HOMA-IR (rs6576507 and rs8026527, 3.7*10(-7)≤P≤1.1*10(-5)) near ATPase, class V, type 10A (ATP10A), and the L Type voltage dependent calcium channel (CACNA1D, rs1401492, P≤5.2*10(-6)). ATP10A belongs to a family of aminophospholipid-transporting ATPases and has been associated with type 2 diabetes in mice. CACNA1D has been linked to pancreatic beta cell generation in mice. The two most significant copy variable markers (rs10277702 and rs361367; P

Acetyl Coenzyme A carboxylase β (ACACB) is the rate-limiting enzyme in fatty acid oxidation, and continuous fatty acid oxidation in Acacb knock-out mice increases insulin sensitivity. Systematic human studies have not been performed to evaluate whether ACACB variants regulate gene expression and insulin sensitivity in skeletal muscle and adipose tissues. We sought to determine whether ACACB transcribed variants were associated with ACACB gene expression and insulin sensitivity in non-diabetic African American (AA) and European American (EA) adults.ACACB transcribed single nucleotide polymorphisms (SNPs) were genotyped in 105 EAs and 46 AAs whose body mass index (BMI), lipid profiles and ACACB gene expression in subcutaneous adipose and skeletal muscle had been measured. Allelic expression imbalance (AEI) was assessed in lymphoblast cell lines from heterozygous subjects in an additional EA sample (n = 95). Selected SNPs were further examined for association with insulin sensitivity in a cohort of 417 EAs and 153 AAs.ACACB transcribed SNP rs2075260 (A/G) was associated with adipose ACACB messenger RNA expression in EAs and AAs (p = 3.8×10(-5), dominant model in meta-analysis, Stouffer method), with the (A) allele representing lower gene expression in adipose and higher insulin sensitivity in EAs (p = 0.04). In EAs, adipose ACACB expression was negatively associated with age and sex-adjusted BMI (r = -0.35, p = 0.0002).Common variants within the ACACB locus appear to regulate adipose gene expression in humans. Body fat (represented by BMI) may further regulate adipose ACACB gene expression in the EA population.

It has been suggested that genetic susceptibility plays an important role in the pathogenesis of diabetic nephropathy. A large-scale genotyping analysis of gene-based single nucleotide polymorphisms (SNPs) in Japanese patients with type 2 diabetes identified the gene encoding acetyl-coenzyme A carboxylase beta (ACACB) as a candidate for a susceptibility to diabetic nephropathy; the landmark SNP was found in the intron 18 of ACACB (rs2268388: intron 18 +4139 C > T, p = 1.4x10(-6), odds ratio = 1.61, 95% confidence interval [CI]: 1.33-1.96). The association of this SNP with diabetic nephropathy was examined in 9 independent studies (4 from Japan including the original study, one Singaporean, one Korean, and two European) with type 2 diabetes. One case-control study involving European patients with type 1 diabetes was included. The frequency of the T allele for SNP rs2268388 was consistently higher among patients with type 2 diabetes and proteinuria. A meta-analysis revealed that rs2268388 was significantly associated with proteinuria in Japanese patients with type 2 diabetes (p = 5.35 x 10(-8), odds ratio = 1.61, 95% Cl: 1.35-1.91). Rs2268388 was also associated with type 2 diabetes-associated end-stage renal disease (ESRD) in European Americans (p = 6 x 10(-4), odds ratio = 1.61, 95% Cl: 1.22-2.13). Significant association was not detected between this SNP and nephropathy in those with type 1 diabetes. A subsequent in vitro functional analysis revealed that a 29-bp DNA fragment, including rs2268388, had significant enhancer activity in cultured human renal proximal tubular epithelial cells. Fragments corresponding to the disease susceptibility allele (T) had higher enhancer activity than those of the major allele. These results suggest that ACACB is a strong candidate for conferring susceptibility for proteinuria in patients with type 2 diabetes.

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes.Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: “Some College” (yes/no, for any education beyond high school) and “Graduated College” (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10−8) and suggestive (p < 1 × 10−6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals).Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue.Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Although the common genetic variants contributing to blood lipid levels have been studied, the contribution of rare variants is less understood. Here, the authors perform a rare coding and noncoding variant association study of blood lipid levels using whole genome sequencing data.

Hypercalcemia is a common and potentially serious electrolyte abnormality that is often observed in patients with chronic kidney disease (CKD). When malignancy is considered, parathyroid hormone-related protein (PTHrP) levels are often measured. PTHrP is produced by cancer cells and mimics the effects of parathyroid hormone (PTH) to elevate serum calcium concentrations. The amino and carboxy termini of PTHrP are of functional relevance. C-terminal PTHrP levels accumulate with CKD and can be elevated in normocalcemic CKD patients who lack malignancy. The existence of amino (N)-terminal and carboxy (C)-terminal PTHrP assays and how their concentrations are impacted by CKD are reviewed herein. The case of a patient on maintenance hemodialysis who developed prolonged hypercalcemia with elevated PTHrP concentrations is presented. The workup revealed suppressed intact PTH, low 25-hydroxyvitamin D, and 1,25-dihydroxyvitamin D levels. The initial PTHrP assay returned elevated. However, it was unappreciated that it was the C-terminal assay and the patient underwent an unnecessary search for malignancy. A subsequent N-terminal PTHrP assay returned within the normal range. Many commercial labs run the C-terminal PTHrP assay as their first-line test. This can lead to inaccurate differential diagnoses in hypercalcemic patients with CKD. We emphasize the need to specifically request N-terminal PTHrP assays in patients with advanced kidney disease when humoral hypercalcemia of malignancy is suspected.

apoM is a minor HDL apolipoprotein and carrier for sphingosine-1-phosphate (S1P). HDL apoM and S1P concentrations are inversely associated with atherosclerosis progression in rodents. We evaluated associations between plasma concentrations of S1P, plasma concentrations of apoM, and HDL apoM levels with prevalent subclinical atherosclerosis and mortality in the African American-Diabetes Heart Study participants (N = 545). Associations between plasma S1P, plasma apoM, and HDL apoM with subclinical atherosclerosis and mortality were assessed using multivariate parametric, nonparametric, and Cox proportional hazards models. At baseline, participants' median (25th percentile, 75th percentile) age was 55 (49, 62) years old and their coronary artery calcium (CAC) mass score was 26.5 (0.0, 346.5). Plasma S1P, plasma apoM, and HDL apoM were not associated with CAC. After 64 (57.6, 70.3) months of follow-up, 81 deaths were recorded. Higher concentrations of plasma S1P [odds ratio (OR) = 0.14, P = 0.01] and plasma apoM (OR = 0.10, P = 0.02), but not HDL apoM (P = 0.89), were associated with lower mortality after adjusting for age, sex, statin use, CAC, kidney function, and albuminuria. We conclude that plasma S1P and apoM concentrations are inversely and independently associated with mortality, but not CAC, in African Americans with type 2 diabetes after accounting for conventional risk factors.

APOL1 gene renal-risk variants are associated with nephropathy and CVD in African Americans; however, little is known about the circulating APOL1 variant proteins which reportedly bind to HDL. We examined whether APOL1 G1 and G2 renal-risk variant serum concentrations or lipoprotein distributions differed from nonrisk G0 APOL1 in African Americans without nephropathy. Serum APOL1 protein concentrations were similar regardless of APOL1 genotype. In addition, serum APOL1 protein was bound to protein complexes in two nonoverlapping peaks, herein referred to as APOL1 complex A (12.2 nm diameter) and complex B (20.0 nm diameter). Neither of these protein complexes associated with HDL or LDL. Proteomic analysis revealed that complex A was composed of APOA1, haptoglobin-related protein (HPR), and complement C3, whereas complex B contained APOA1, HPR, IgM, and fibronectin. Serum HPR was less abundant on complex B in individuals with G1 and G2 renal-risk variant genotypes, relative to G0 (P = 0.0002–0.037). These circulating complexes may play roles in HDL metabolism and susceptibility to CVD.

Two APOL1 gene variants, which likely evolved to protect individuals from African sleeping sickness, are strongly associated with nondiabetic kidney disease in individuals with recent African ancestry. Consistent with its role in trypanosome killing, the pro-death APOL1 protein is toxic to most cells, but its mechanism of cell death is poorly understood and little is known regarding its intracellular trafficking and secretion. Because the liver appears to be the main source of circulating APOL1, we examined its secretory behavior and mechanism of toxicity in hepatoma cells and primary human hepatocytes. APOL1 is poorly secreted in vitro, even in the presence of chemical chaper­ones; however, it is efficiently secreted in wild-type transgenic mice, suggesting that APOL1 secretion has specialized requirements that cultured cells fail to support. In hepatoma cells, inducible expression of APOL1 and its risk variants promoted cell death, with the G1 variant displaying the highest degree of toxicity. To explore the basis for APOL1-mediated cell toxicity, endoplasmic reticulum stress, pyroptosis, autophagy, and apoptosis were examined. Our results suggest that autophagy represents the predominant mechanism of APOL1-mediated cell death. Overall, these results increase our understanding of the basic biology and trafficking behavior of circulating APOL1 from the liver.