Analysis of the ImmunoChip single nucleotide polymorphism (SNP) array in 2816 individuals, comprising the most common subtypes (oligoarticular and RF negative polyarticular) of juvenile idiopathic arthritis (JIA) and 13056 controls strengthens the evidence for association to three known JIA-risk loci (HLA, PTPN22 and PTPN2) and has identified fourteen risk loci reaching genome-wide significance (p < 5 × 10-8) for the first time. Eleven additional novel regions showed suggestive evidence for association with JIA (p < 1 × 10-6). Dense-mapping of loci along with bioinformatic analysis has refined the association to one gene for eight regions, highlighting crucial pathways, including the IL-2 pathway, in JIA disease pathogenesis. The entire ImmunoChip loci, HLA region and the top 27 loci (p < 1 × 10-6) explain an estimated 18%, 13% and 6% risk of JIA, respectively. Analysis of the ImmunoChip dataset, the largest cohort of JIA cases investigated to date, provides new insight in understanding the genetic basis for this childhood autoimmune disease.
Complex genetic disorders are a result of a combination of genetic and non-genetic factors, all potentially interacting. Machine learning methods hold the potential to identify multi-locus and environmental associations thought to drive complex genetic traits. Decision trees, a popular machine learning technique, offer a computationally low complexity algorithm capable of detecting associated sets of SNPs of arbitrary size, including modern genome-wide SNP scans. However, interpretation of the importance of an individual SNP within these trees can present challenges.
We present a new decision tree algorithm denoted as Bagged Alternating Decision Trees (BADTrees) that is based on identifying common structural elements in a bootstrapped set of ADTrees. The algorithm is order nk2, where n is the number of SNPs considered and k is the number of SNPs in the tree constructed. Our simulation study suggests that BADTrees have higher power and lower type I error rates than ADTrees alone and comparable power with lower type I error rates compared to logistic regression. We illustrate the application of these data using simulated data as well as from the Lupus Large Association Study 1 (7822 SNPs in 3548 individuals). Our results suggest that BADTrees holds promise as a low computational order algorithm for detecting complex combinations of SNP and environmental factors associated with disease.
Machine Learning; Genetic Association; Gene-Gene Interaction; Multi-locus Models
African Americans (AAs) are predisposed to non-diabetic (non-DM) end-stage renal disease (ESRD) and studies have shown a genetic component to this risk. Rare mutations in ACTN4 (α-actinin-4) an actin binding protein expressed in podocytes cause familial focal segmental glomerulosclerosis.
We assessed the contribution of coding variants in ACTN4 to non-DM ESRD risk in AAs. Nineteen exons, 2800 bases of the promoter and 392 bases of the 3’ untranslated region of ACTN4 were sequenced in 96 AA non-DM ESRD cases and 96 non-nephropathy controls (384 chromosomes). Sixty-seven single nucleotide polymorphisms (SNPs) including 51 novel SNPs were identified. The SNPs comprised 33 intronic, 21 promoter, 12 exonic, and 1 3’ variant. Sixty-two of the SNPs were genotyped in 296 AA non-DM ESRD cases and 358 non-nephropathy controls.
One SNP, rs10404257, was associated with non-DM ESRD (p<1.0E-4, odds ratio (OR)=0.76, confidence interval (CI)=0.59–0.98; additive model). Forty-seven SNPs had minor allele frequencies less than 5%. These SNPs were segregated into risk and protective SNPs and each category was collapsed into a single marker, designated by the presence or absence of any rare allele. The presence of any rare allele at a risk SNP was significantly associated with non-DM ESRD (p = 0.001, dominant model). The SNPs with the strongest evidence for association (n = 20) were genotyped in an independent set of 467 non-DM ESRD cases and 279 controls. Although, rs10404257 was not associated in this replication sample, when the samples were combined rs10404257 was modestly associated (p=0.032, OR=0.78, CI=0.63–0.98; dominant model). SNPs were tested for interaction with markers in the APOL1 gene, previously associated with non-DM ESRD in AAs and rs10404257 was modestly associated (p = 0.0261, additive model).
This detailed evaluation of ACTN4 variation revealed limited evidence of association with non-DM ESRD in AAs.
ACTN4; non-diabetic ESRD; FSGS; kidney; hypertensive nephrosclerosis; African Americans
Genome-wide association studies have identified novel genetic factors that contribute to intracranial aneurysm (IA) susceptibility. We sought to confirm previously reported loci, to identify novel risk factors, and to evaluate the contribution of these factors to familial and sporadic IA.
We utilized 2 complementary samples, one recruited on the basis of a dense family history of IA (discovery sample 1: 388 IA cases and 397 controls) and the other without regard to family history (discovery sample 2: 1095 IA cases and 1286 controls). Imputation was used to generate a common set of single nucleotide polymorphisms (SNP) across samples, and a logistic regression model was used to test for association in each sample. Results from each sample were then combined in a meta-analysis.
There was only modest overlap in the association results obtained in the 2 samples. In neither sample did results reach genome-wide significance. However, the meta-analysis yielded genome-wide significance for SNP on chromosome 9p (CDKN2BAS; rs6475606; P=3.6×10−8) and provided further evidence to support the previously reported association of IA with SNP in SOX17 on chromosome 8q (rs1072737; P=8.7×10−5). Analyses suggest that the effect of smoking acts multiplicatively with the SNP genotype, and smoking has a greater effect on risk than SNP genotype.
In addition to replicating several previously reported loci, we provide further evidence that the association on chromosome 9p is attributable to variants in CDKN2BAS (also known as ANRIL, an antisense noncoding RNA).
genome-wide association study; intracranial aneurysm
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.
Little 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. This study analyzed variation in reactive intermediate genes for association with SLE in two populations with African ancestry.
A total of 244 SNPs 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) and. Single-marker, haplotype, and two-locus interaction tests were computed for these populations.
The glutathione reductase gene GSR (rs2253409, P=0.0014, OR [95% CI]=1.26 [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 [95%CI]=2.10 [1.23–3.59]) and nitric oxide synthase gene NOS1 (rs561712, P=0.0072, OR [95%CI]=0.62 [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 [95%CI]=1.26 [1.10–1.44]). Haplotype and two-locus interaction analyses also uncovered different loci in each population.
These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
We have conducted a GWAS in a Caucasian cohort of juvenile idiopathic arthritis (JIA) patients and have previously published findings limited to autoimmune loci shared with other diseases. The goal of this study was to identify novel JIA predisposing loci using genome-wide approaches.
The Discovery cohort consisted of Caucasian JIA cases (814) and local controls (658) genotyped on the Affymetrix SNP 6.0 Array along with 2400 out-of-study controls. A replication study consisted of 10 SNPs genotyped in 1744 cases and 7010 controls from the US and Europe.
Analysis within the Discovery cohort provided evidence of associations at 3q13 within C3orf1 and near CD80 (rs4688011, OR=1.37, P=1.88×10−6), and 10q21 near the gene JMJD1C [rs6479891, odds ratio (OR) =1.59, P=6.1×10−8; rs12411988, OR=1.57, P=1.16×10−7 and rs10995450, OR = 1.31, P=6.74×10−5]. Meta-analysis continued to provide evidence for association for these 4 SNPs (rs4688011, P=3.6×10−7, rs6479891, P=4.33×10−5; rs12411988, P=2.71×10−5; and rs10995450, 5.39×10−5;). Gene expression data from 68 JIA cases and 23 local controls showed cis eQTL associations for C3orf1 SNP rs4688011 (P=0.024 or P=0.034, depending on probe set) and the JMJD1C SNPs rs6479891 and rs12411988 (P=0.01 and P=0.008, respectively). A variance component liability model estimated that common SNP variation accounts for ~1/3 of JIA susceptibility.
Genetic association results and correlated gene expression findings provide evidence of association at 3q13 and 10q21 for JIA and offer novel genes as plausible candidates in disease pathology.
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.
Systemic lupus erythematosus (SLE/lupus) is a complex disease in which the body's immune cells cause inflammation in one or more systems to cause the associated morbidity. Hormones, the environment and genes are all causal contributors to SLE and over the past several years the genetic component of SLE has been firmly established. Several genes which are regulators of the immune system are associated with disease risk. We have established one of these, the tumour-necrosis family superfamily member 4 (TNFSF4) gene, as a lupus susceptibility gene in Northern Europeans. A major obstacle in pinpointing the marker(s) at TNFSF4 which best explain the risk of SLE has been the strong correlation (linkage disequilibrium, LD) between adjacent markers across the TNFSF4 region in this population. To address this, we have typed polymorphisms in several populations in addition to the European groups. The mixed ancestry of these populations gives a different LD pattern than that found in Europeans, presenting a method of pinpointing the section of the TNFSF4 region which results in SLE susceptibility. The Non-European populations have allowed identification of a polymorphism likely to regulate expression of TNFSF4 to increase susceptibility to SLE.
Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (λS = ~30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 × 10−7 < Poverall < 1.6 × 10−23; odds ratio 0.82–1.62)in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 × 10−5) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at ≥9 other loci (P < 2 × 10−7). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.
Non-alcoholic fatty liver disease (NAFLD) is a highly prevalent condition, particularly among Hispanic Americans. A genetic variant in PNPLA3 (rs738409) has been identified as a strong predictor of hepatic fat content.
To confirm the association of this variant with NAFLD in two minority cohorts, Hispanic Americans and African Americans, in whom liver density was quantified by computed tomography (CT).
This analysis was conducted in the Insulin Resistance Atherosclerosis (IRAS) Family Study. Participants were recruited from the general community and included 843 Hispanic American and 371 African American adults aged 18–81 years. A single variant in PNPLA3 (rs738409) was genotyped. Liver density was calculated in Hounsfield Units from abdominal CT scans.
Single nucleotide polymorphism (SNP) rs738409 was strongly associated with reduced liver density (i.e. NAFLD) in Hispanic Americans (1.18 × 10−9) and in African Americans (P = 4.99 × 10−6). The association followed an additive genetic model with the G allele conferring risk. The allele was two times more common in Hispanic Americans than in African Americans (40 vs 19%), consistent with the greater prevalence of NAFLD in Hispanic Americans (24 vs 9%). The SNP explained 4.4 and 5.6% of the variance of the adjusted liver density outcome in Hispanic Americans and African Americans, respectively.
We confirmed the association of a PNPLA3 variant with NAFLD in Hispanic Americans and African Americans, suggesting that PNPLA3 contributes to the variation in NAFLD across multiple ethnicities. This study adds to the growing evidence that some of the ethnic variation in NAFLD is genetic.
African Americans; computed tomography; genetic epidemiology; hepatic steatosis; Hispanic Americans; non-alcoholic fatty liver disease; PNPLA3
Despite intensive anti-hypertensive therapy there was a high incidence of renal end-points in participants of the African American Study of Kidney Disease and Hypertension (AASK) cohort. To better understand this, coding variants in the apolipoprotein L1 (APOL1) and the non-muscle myosin heavy chain 9 (MYH9) genes were evaluated for an association with hypertension-attributed nephropathy and clinical outcomes in a case-control study. Clinical data and DNA were available for 675 AASK participant cases and 618 African American non-nephropathy control individuals. APOL1 G1 and G2, and MYH9 E1 variants along with 44 ancestry informative markers were genotyped with allele frequency differences between cases and controls analyzed by logistic regression multivariable models adjusting for ancestry, age, and gender. In recessive models, APOL1 risk variants were significantly associated with kidney disease in all cases compared to controls with an odds ratio of 2.57. In AASK cases with more advanced disease, such as a baseline urine protein to creatinine ratio over 0.6 g/g or a serum creatinine over 3 mg/dL during follow-up, the association was strengthened with odds ratios of 6.29 and 4.61, respectively. APOL1 risk variants were consistently associated with renal disease progression across medication classes and blood pressure targets. Thus, kidney disease in AASK participants was strongly associated with APOL1 renal risk variants.
Patients with type 2 diabetes (T2D) are at elevated risk for cardiovascular disease (CVD) events and mortality. Recent studies have assessed the impact of genetic variants affecting high-density lipoprotein cholesterol (HDL) concentrations on CVD risk in the general population. This study examined the utility of HDL-associated single nucleotide polymorphisms (SNPs) for CVD risk prediction in European Americans with T2D enrolled in the Diabetes Heart Study (DHS).
Genetic risk scores (GRS) of HDL-associated SNPs were constructed and evaluated for potential associations with mortality and with coronary artery calcified atherosclerotic plaque (CAC), a measure of subclinical CVD strongly associated with CVD events and mortality. Two sets of SNPs were used to construct GRS; while all SNPs were selected primarily for their impacts on HDL, one set of SNPs had pleiotropic effects on other lipid parameters, while the other set lacked effects on low-density lipoprotein cholesterol (LDL) or triglyceride concentrations.
The GRS were specifically associated with HDL concentrations (4.90 × 10-7 < p < 0.02) in models adjusted for age, sex, and body mass index (BMI), but were not associated with LDL or triglycerides. Cox proportional hazards regression analysis suggested the HDL-associated GRS had no impact on risk of CVD-mortality (0.48 < p < 0.99) in models adjusted for other known CVD risk factors. However, associations between several of the GRS and CAC were observed (3.85 × 10-4 < p < 0.03) in models adjusted for other known CVD risk factors.
The GRS analyzed in this study provide a tool for assessment of HDL-associated SNPs and their impact on CVD risk in T2D. The observed associations between several of the GRS and CAC suggest a potential role for HDL-associated SNPs on subclinical CVD risk in patients with T2D.
High-density lipoprotein cholesterol; Type 2 diabetes; Coronary artery calcified plaque; Mortality; Genetic risk score
Recently, a genome-wide association scan was completed in the IRAS (Insulin Resistance Atherosclerosis Study) Family Study (IRASFS) Hispanic-American cohort. Multiple single-nucleotide polymorphisms (SNPs) in the G-protein signaling 6 (RGS6) gene were found to be associated with adiposity phenotypes. RGS6 has shown downstream antagonistic interplay with opioid receptors, targets of fatty/sugary food agonists. The possibility that RGS6 promotes tolerance and tachyphylaxis among the opioid receptor is a plausible pathway for overconsuming fat/sugar-laden food. Therefore, we hypothesized that RGS6 variants are associated with intake of fatty/sugary foods. In 932 Hispanics from San Antonio and San Luis Valley, CO, the following dietary intake variables were assessed using the Block Brief 2000 food frequency questionnaire: total calories, total fat, % calories from fat, % calories from saturated fat, protein, % calories from protein, carbohydrates, % calories from carbohydrates, and daily frequency of servings of fats/oils/sweets. We tested for association between 23 SNPs in RGS6 and dietary intake using a variance components measured genotype approach. All models were adjusted for gender, recruitment site, admixture, BMI, and age. Using an additive genetic model, rs1402064 was associated with higher intake of fats/oils/sweets, total calories, total fat and saturated fat (P = 0.0007, 0.026, 0.023, and 0.024). SNPs rs847330 and rs847354 were associated with higher intake of fats/oils/sweets (P = 0.002 and 0.018), total fat (P = 0.040 and 0.048) and saturated fat (P = 0.044 and 0.041). Finally, rs769148 was associated with higher intake of fats/oils/sweets (P = 0.002). RGS6 is a new candidate gene for adiposity traits that may be associated with a behavioral tendency toward fat-laden food intake.
To investigate and refine the relationships among systemic lupus erythematosus (SLE) and related autoantibodies, interferon-α (IFN-α), and various ancestral backgrounds.
We investigated quantitatively defined genetic ancestry through principal component analysis in place of self-reported ancestry.
African ancestry was found to be associated with presence of anti-RNP antibody (p = 0.0026), and anti-RNP was correlated with high levels of IFN-α (p = 2.8 × 10−5).
Our data support a model in which African ancestry increases the likelihood of SLE-associated autoantibody formation, which subsequently results in higher levels of serum IFN-α.
SYSTEMIC LUPUS ERYTHEMATOSUS; AUTOANTIBODIES; INTERFERONS; GENETICS
Genetic variants in myocardial sodium and potassium channel genes are associated with prolonged QT interval and increased risk of sudden death. It is unclear whether these genetic variants remain relevant in subjects with underlying conditions such as diabetes that are associated with prolonged QT interval.
We tested single nucleotide polymorphisms (SNPs) in five candidate genes for association with QT interval in a family-based study of subjects with type 2 diabetes mellitus (T2DM). Thirty-six previously reported SNPs were genotyped in KCNQ1, HERG, SCN5A, KCNE1, and KCNE2 in 901 European Americans from 366 families. The heart rate-corrected (QTc) durations were determined using the Marquette 12SL program. Associations between the QTc interval and the genotypes were evaluated using SOLAR adjusting for age, gender, T2DM status, and body mass index.
Within KCNQ1 there was weak evidence for association between the minor allele of IVS12+14T>C and increased QTc (p=0.02). The minor allele of rs2236609 in KCNE1 trended toward significance with longer QTc (p=0.06), while the minor allele of rs1805123 in HERG trended toward significance with shorter QTc (p=0.07). However, no statistically significant associations were observed between the remaining SNPs and QTc variation.
We found weak evidence of association between three previously-reported SNPs and QTc interval duration. While it appears as though genetic variants in previously identified candidate genes may be associated with QT duration in subjects with diabetes, the clinical implications of these associations in diabetic subjects at high risk for sudden death remains to be determined.
QT interval; diabetes; association study; genetics; ion channels
Systemic lupus erythematosus (SLE) is a sexually dimorphic autoimmune disease which is more common in women, but affected men often experience a more severe disease. The genetic basis of sexual dimorphism in SLE is not clearly defined. A study was undertaken to examine sex-specific genetic effects among SLE susceptibility loci.
A total of 18 autosomal genetic susceptibility loci for SLE were genotyped in a large set of patients with SLE and controls of European descent, consisting of 5932 female and 1495 male samples. Sex-specific genetic association analyses were performed. The sex–gene interaction was further validated using parametric and nonparametric methods. Aggregate differences in sex-specific genetic risk were examined by calculating a cumulative genetic risk score for SLE in each individual and comparing the average genetic risk between male and female patients.
A significantly higher cumulative genetic risk for SLE was observed in men than in women. (P = 4.52×10−8) A significant sex–gene interaction was seen primarily in the human leucocyte antigen (HLA) region but also in IRF5, whereby men with SLE possess a significantly higher frequency of risk alleles than women. The genetic effect observed in KIAA1542 is specific to women with SLE and does not seem to have a role in men.
The data indicate that men require a higher cumulative genetic load than women to develop SLE. These observations suggest that sex bias in autoimmunity could be influenced by autosomal genetic susceptibility loci.
Polymorphisms in the non-muscle myosin IIA gene (MYH9) are associated with focal segmental glomerulosclerosis (FSGS) and non-diabetic end-stage renal disease (ESRD) in African Americans and FSGS in European Americans. We tested for association of single nucleotide polymorphisms (SNPs) in MYH9 with T2DM–ESRD in European Americans; additionally, three APOL1 gene variants were evaluated.
Fifteen MYH9 SNPs and two APOL1 SNPs plus a 6-bp deletion were genotyped in 1963 European Americans, 536 cases with T2DM–ESRD and 1427 non-nephropathy controls (467 with T2DM and 960 without diabetes).
Comparing T2DM–ESRD cases with the 467 T2DM non-nephropathy controls, single variant associations trending toward significance were detected with SNPs rs4821480, rs2032487 and rs4281481 comprising part of the major MYH9 E1 risk haplotype [P-values 0.053–0.055 recessive, odds ratio (OR) 6.08–6.14]. Comparing T2DM–ESRD cases to all 1427 non-nephropathy controls, we confirmed evidence of association in these three SNPs as well as in the fourth E1 SNP (rs3752462) (P-values 0.017–0.035, OR 1.41–3.72). APOL1 G1/G2 nephropathy risk variants were rare in individuals of European American heritage, present in 0.28% of chromosomes in T2DM–ESRD cases and 0.32% of controls.
MYH9 SNPs rs4821480, rs2032487, rs4281481 and rs3752462 are associated with T2DM–ESRD susceptibility in European Americans. The APOL1 risk variants are not present at appreciable frequency in this cohort with T2DM–ESRD. Therefore, polymorphisms in MYH9 appear to influence nephropathy risk in this sample.
APOL1; diabetic nephropathy; end-stage renal disease; MYH9; type 2 diabetes mellitus
African American; APOL1; end-stage renal disease; FSGS; hypertensive nephrosclerosis
Familial aggregation of non-diabetic end stage renal disease (ESRD) is found in African Americans and variants in the apolipoprotein L1 gene (APOL1) contribute to this risk. To detect genetic associations with milder forms of nephropathy in high-risk families, analyses were performed using generalized estimating equations to assess relationships between kidney disease phenotypes and APOL1 variants in 786 relatives of 470 families. Adjusting for familial correlations, 23.1, 46.7, and 30.2 percent of genotyped relatives possessed two, one, or no APOL1 risk variants, respectively. Relatives with two compared to one or no risk variants had statistically indistinguishable median systolic blood pressure, urine albumin to creatinine ratio, estimated GFR (MDRD equation) and serum cystatin C levels. After adjusting for age, gender, age at ESRD in families, and African ancestry, significant associations were detected between APOL1 with overt proteinuria and estimated GFR (CKD-EPI equation), with a trend toward significance for quantitative albuminuria. Thus, relatives of African Americans with non-diabetic ESRD are enriched for APOL1 risk variants. After adjustment, two APOL1 risk variants weakly predict mild forms of kidney disease. Second hits appear necessary for the initiation of APOL1-associated nephropathy.
African American; APOL1; end-stage renal disease; FSGS; kidney; screening
Cardiac manifestations of neonatal lupus, comprising atrioventricular conduction defects and cardiomyopathy, occur in fetuses exposed to anti-Ro/SSA antibodies, and carry substantial mortality. There is strong evidence of a genetic contribution to the risk. This study was undertaken to evaluate single-nucleotide polymorphisms (SNPs) for associations with cardiac neonatal lupus.
Children of European ancestry with cardiac neonatal lupus (n = 116) were genotyped using the Illumina 370K SNP platform and merged with 3,351 controls. Odds ratios (ORs) and 95% confidence intervals (95% CIs) for association with cardiac neonatal lupus were determined.
The 17 most significant associations with cardiac neonatal lupus were found in the HLA region. The region near the MICB gene showed the strongest variant (rs3099844; Pdom = 4.52 × 10−10, OR 3.34 [95% CI 2.29–4.89]), followed by a missense variant within C6orf10 (rs7775397; Pdom = 1.35 × 10−9, OR 3.30), which lies between NOTCH4 and BTNL2, and several SNPs near the tumor necrosis factor α gene, including rs2857595 (Padd = 1.96 × 10−9, OR 2.37), rs2230365 (Padd = 1.00 × 10−3, OR 0.46), and rs3128982 (Padd = 6.40 × 10−6, OR 1.86). Outside the HLA region, an association was detected at 21q22, upstream of the transcription regulator ets-related isoform 1 (rs743446; P = 5.45 × 10−6, OR 2.40). HLA notwithstanding, no individual locus previously implicated in autoimmune diseases achieved genome-wide significance.
These results suggest that variation near genes related to inflammatory and apoptotic responses may promote cardiac injury initiated by passively acquired autoantibodies.
Neonatal lupus (NL) occurs in fetuses exposed to maternal anti-SSA/Ro and/or SSB/La antibodies, although the mothers may not manifest any clinical disease. A focus on transmission of risk factors for NL from maternal grandparents to mothers may yield dividends towards understanding the aggregation of autoantibodies and genetic factors in families.
51 mothers of children with cardiac and/or cutaneous NL, 48 maternal grandmothers and 35 maternal grandfathers in the Research Registry for Neonatal Lupus were interrogated for clinical symptoms (questionnaire) and laboratory assessments including anti-SSA/Ro and SSB/La antibodies (ELISA) and genotype at rs1800629 (TNFα-308) and rs7775397 (C6orf10) (allelic discrimination). The transmission disequilibrium test (TDT) was computed to test for nonrandom transmission from maternal grandparents to the NL mothers.
The phenotypic feature held in common in NL mothers was the autoantibody, and not the clinical profile; 7 had lupus, 14 had Sjogren's syndrome, 7 had both, and 23 were asymptomatic. NL mothers were significantly enriched in variant allelic frequencies at both TNFα-308 and C6orf10. The grandparents of NL children carried minimal burden for autoimmune disease or abnormal antibody production and were not enriched in the genetic risk factors. However, the TDT analysis showed significant excess transmission of the risk alleles at both TNFα-308 (P=3.93×10−4, OR=6.67) and C6orf10 (P=3.74×10−5, OR=35) to NL mothers.
NL mothers are enriched for the TNFα-308 and C6orf10 variant alleles, which are preferentially inherited from the asymptomatic maternal grandparents. These findings support the hypothesis that the development of NL and genetic etiology are multigenerational.
The 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.
We fine-mapped ≥136 SNPs in a ~227kb region on Xq28, containing IRAK1, MECP2 and 7 adjacent genes (L1CAM, AVPR2, ARHGAP4, NAA10, RENBP, HCFC1 and TMEM187), for association with SLE in 15,783 case-control subjects derived from 4 different ancestral groups.
Multiple SNPs showed strong association with SLE in European Americans, Asians and Hispanics at P<5×10−8 with consistent association in subjects with African ancestry. Of these, 6 SNPs located in the TMEM187-IRAK1-MECP2 region captured the underlying causal variant(s) residing in a common risk haplotype shared by all 4 ancestral groups. Among them, rs1059702 best explained the Xq28 association signals in conditional testings and exhibited the strongest P value in trans-ancestral meta-analysis (Pmeta=1.3×10−27, OR=1.43), and thus was considered to be the most-likely causal variant. The risk allele of rs1059702 results in the amino acid substitution S196F in IRAK1 and had previously been shown to increase NF-κB activity in vitro. We also found that the homozygous risk genotype of rs1059702 was associated with lower mRNA levels of MECP2, but not IRAK1, in SLE patients (P=0.0012) and healthy controls (P=0.0064).
These data suggest contributions of both IRAK1 and MECP2 to SLE susceptibility.
Systemic Lupus Erythematosus; Gene Polymorphism; Xq28; IRAK1; MECP2
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 (R2 = 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.
Systemic lupus erythematosus (SLE) is a debilitating autoimmune disease contributed to by excessive innate immune activation involving toll-like receptors (TLRs, particularly TLR7/8/9) and type I interferon (IFN) signaling pathways. TLR7 responds against RNA–containing nuclear antigens and activates IFN-α pathway, playing a pivotal role in the development of SLE. While a genomic duplication of Tlr7 promotes lupus-like disease in the Y-linked autoimmune accelerator (Yaa) murine model, the lack of common copy number variations at TLR7 in humans led us to identify a functional single nucleotide polymorphism (SNP), rs3853839 at 3′ UTR of the TLR7 gene, associated with SLE susceptibility in Eastern Asians. In this study, we fine-mapped the TLR7-TLR8 region and confirmed rs3853839 exhibiting the strongest association with SLE in European Americans, African Americans, and Amerindian/Hispanics. Individuals carrying the risk G allele of rs3853839 exhibited increased TLR7 expression at the both mRNA and protein level and decreased transcript degradation. MicroRNA-3148 (miR-3148) downregulated the expression of non-risk allele (C) containing transcripts preferentially, suggesting a likely mechanism for increased TLR7 levels in risk-allele carriers. This trans-ancestral mapping provides evidence for the global association with SLE risk at rs3853839, which resides in a microRNA–gene regulatory site affecting TLR7 expression.
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 [Pmeta = 5.20×10−14; odds ratio, 95% confidence interval = 0.82 (0.78–0.87)], and two missense variants, rs1990760 (Ala946Thr) [Pmeta = 3.08×10−7; 0.88 (0.84–0.93)] and rs10930046 (Arg460His) [Pdom = 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.
African-Americans (AA) are at increased risk of systemic lupus erythematosus (SLE), but the genetic basis of this risk increase is largely unknown. We used admixture mapping to localize disease-causing genetic variants that differ in frequency across populations. This approach is advantageous for localizing susceptibility genes in recently admixed populations like AA. Our genome-wide admixture scan identified seven admixture signals, and we followed the best signal at 2q22–24 with fine-mapping, imputation-based association analysis and experimental validation. We identified two independent coding variants and a non-coding variant within the IFIH1 gene associated with SLE. 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.
Haptoglobin (HP) is an acute phase protein that binds to freely circulating hemoglobin. HP exists as two distinct forms, HP1 and HP2. The longer HP2 form has been associated with cardiovascular (CVD) events and mortality in individuals with type 2 diabetes (T2DM).
This study examined the association of HP genotypes with subclinical CVD, T2DM risk, and associated risk factors in a T2DM-enriched sample. Haptoglobin genotypes were determined in 1208 European Americans (EA) from 473 Diabetes Heart Study (DHS) families via PCR. Three promoter SNPs (rs5467, rs5470, and rs5471) were also genotyped.
Analyses revealed association between HP2-2 duplication and increased carotid intima-media thickness (IMT; p = 0.001). No association between HP and measures of calcified arterial plaque were observed, but the HP polymorphism was associated with triglyceride concentrations (p = 0.005) and CVD mortality (p = 0.04). We found that the HP2-2 genotype was associated with increased T2DM risk with an odds ratio (OR) of 1.49 (95% CI 1.18-1.86, p = 6.59x10-4). Promoter SNPs were not associated with any traits.
This study suggests association between the HP duplication and IMT, triglycerides, CVD mortality, and T2DM in an EA population enriched for T2DM. Lack of association with atherosclerotic calcified plaque likely reflect differences in the pathogenesis of these CVD phenotypes. HP variation may contribute to the heritable risk for CVD complications in T2DM.
Haptoglobin; Genetic polymorphism; Cardiovascular disease; Type 2 diabetes