Pyoderma gangrenosum (PG) and erythema nodosum (EN) are the most common cutaneous manifestations of inflammatory bowel disease (IBD) but little is known regarding their etiopathogenesis.
We performed a case control study comparing characteristics between IBD patients with a documented episode of PG (PG+) and/or EN (EN+) with those without PG (PG-) and EN (EN-). Data on clinical features were obtained by chart review. IBD related serology was determined using ELISA and genome-wide data generated using Illumina technology. Standard statistical tests for association were used.
We identified a total 92 cases of PG and 103 cases of EN with genetic and clinical characteristics, of which 64 PG and 55 EN were available for serological analyses. Fewer male subjects were identified in the PG(+) (OR 0.6, p=0.009) and EN(+) groups (OR 0.31, p=0<0.0001). Colonic disease, previous IBD related surgery, and non-cutaneous extra-intestinal manifestations were more common among both PG(+) and EN(+) patients compared to controls. PG(+) was associated with ANCA seropositivity (p=0.03) and higher ANCA level (p=0.02) in CD. Genetic associations with PG included known IBD loci (IL8RA (p=0.00003), and PRDM1 (0.03)) as well as with USP15 (4.8×10−6) and TIMP3 (5.6 ×10−7). Genetic associations with EN included known IBD susceptibility genes (PTGER4 (p=8.8×10−4), ITGAL (0.03)) as well as SOCS5 (9.64×10−6), CD207 3.14×10−6), ITGB3 (7.56×10−6) and rs6828740 (4q26)(p <5.0 × 10−8). Multivariable models using clinical, serologic, and genetic parameters predicted PG (AUC 0.8) and EN (AUC 0.97).
Cutaneous manifestations in IBD are associated with distinctive genetic characteristics as well as with the similar clinical characteristics including the development of other extra-intestinal manifestations suggesting shared and distinct etiologies.
Crohn's disease; extra-intestinal manifestations; inflammatory bowel disease; pyoderma gangrenosum; erythema nodosum; ulcerative colitis; genetics; serology
Background & Aims
Genetic susceptibility loci for Crohn’s disease (CD) are numerous, complex, and likely interact with undefined components of the environment. It has been a challenge to link the effects of particular loci to phenotypes of cells associated with pathogenesis of CD, such as Paneth cells. We investigated whether specific phenotypes of Paneth cells associated with particular genetic susceptibility loci can be used to define specific subtypes of CD.
We performed a retrospective analysis of 119 resection specimens collected from patients with CD at 2 separate medical centers. Paneth cell phenotypes were classified as normal or abnormal (with disordered, diminished, diffuse, or excluded granule phenotypes) based on lysozyme-positive secretory granule morphology. To uncover the molecular basis of the Paneth cell phenotypes, we developed methods to determine transcriptional profiles from whole-thickness and laser-capture microdissected, formalin-fixed, paraffin-embedded tissue sections.
The proportion of abnormal Paneth cells was associated with the number of CD-associated NOD2 risk alleles. The cumulative number of NOD2 and ATG16L1 risk alleles had an additive effect on the proportion of abnormal Paneth cells. Unsupervised clustering analysis of demographic and Paneth cell data divided patients into 2 principal subgroups, defined by high and low proportions of abnormal Paneth cells. The disordered and diffuse abnormal Paneth cell phenotypes were associated with an altered transcriptional signature of immune system activation. We observed an inverse correlation between abnormal Paneth cells and the presence of granuloma. Moreover, high proportions of abnormal Paneth cells were associated with shorter time to disease recurrence after surgery.
Histologic analysis of Paneth cell phenotypes can be used to divide patients with CD into subgroups with distinct pathognomonic and clinical features.
pathogenesis; prognostic factor; diagnosis; inflammatory bowel disease
The Puerto Rico population may be modeled as an admixed population with contributions from three continents: Sub-Saharan Africa, Ancient America, and Europe. Extending the study of the genetics of inflammatory bowel disease (IBD) to an admixed population such as Puerto Rico has the potential to shed light on IBD genes identified in studies of European populations, find new genes contributing to IBD susceptibility, and provide basic information on IBD for the care of US patients of Puerto Rican and Latino descent. In order to study the association between immune-related genes and Crohn’s disease (CD) and ulcerative colitis (UC) in Puerto Rico, we genotyped 1159 Puerto Rican cases, controls, and family members with the ImmunoChip. We also genotyped 832 subjects from the Human Genome Diversity Panel to provide data for estimation of global and local continental ancestry. Association of SNPs was tested by logistic regression corrected for global continental descent and family structure. We observed the association between Crohn’s disease and NOD2 (rs17313265, 0.28 in CD, 0.19 in controls, OR 1.5, p = 9×10−6) and IL23R (rs11209026, 0.026 in CD, 0.0.071 in controls, OR 0.4, p = 3.8×10−4). The haplotype structure of both regions resembled that reported for European populations and “local” continental ancestry of the IL23R gene was almost entirely of European descent. We also observed suggestive evidence for the association of the BAZ1A promoter SNP with CD (rs1200332, 0.45 in CD, 0.35 in controls, OR 1.5, p = 2×10−6). Our estimate of continental ancestry surrounding this SNP suggested an origin in Ancient America for this putative susceptibility region. Our observations underscored the great difference between global continental ancestry and local continental ancestry at the level of the individual gene, particularly for immune-related loci.
Estimates of the heritability of plasma fibrinogen concentration, an established predictor of cardiovascular disease (CVD), range from 34 to 50%. Genetic variants so far identified by genome-wide association (GWA) studies only explain a small proportion (< 2%) of its variation.
Methods and Results
We conducted a meta-analysis of 28 GWA studies, including more than 90,000 subjects of European ancestry, the first GWA meta-analysis of fibrinogen levels in 7 African Americans studies totaling 8,289 samples, and a GWA study in Hispanic-Americans totaling 1,366 samples. Evaluation for association of SNPs with clinical outcomes included a total of 40,695 cases and 85,582 controls for coronary artery disease (CAD), 4,752 cases and 24,030 controls for stroke, and 3,208 cases and 46,167 controls for venous thromboembolism (VTE). Overall, we identified 24 genome-wide significant (P<5×10−8) independent signals in 23 loci, including 15 novel associations, together accounting for 3.7% of plasma fibrinogen variation. Gene-set enrichment analysis highlighted key roles in fibrinogen regulation for the three structural fibrinogen genes and pathways related to inflammation, adipocytokines and thyrotrophin-releasing hormone signaling. Whereas lead SNPs in a few loci were significantly associated with CAD, the combined effect of all 24 fibrinogen-associated lead SNPs was not significant for CAD, stroke or VTE.
We identify 23 robustly associated fibrinogen loci, 15 of which are new. Clinical outcome analysis of these loci does not support a causal relationship between circulating levels of fibrinogen and CAD, stroke or VTE.
Fibrinogen; cardiovascular disease; genome-wide association study
The ability to identify patients with Crohn’s disease (CD) at highest risk of surgery would be invaluable in guiding therapy. Genome-wide association (GWA) studies have identified multiple IBD loci with unknown phenotypic consequences.
1) to identify associations between known and novel CD loci with early resective CD surgery2) to develop the best predictive model for time to surgery using a combination of phenotypic, serologic and genetic variables.
Genotyping was performed on 1115 subjects using Illumina-based Genome-wide technology. Univariate and multivariate analyses tested genetic associations with need for surgery within 5 years. Analyses were performed by testing known CD loci (n=71) and by performing a GWA study. Time to surgery was analyzed using Cox regression modeling. Clinical and serologic variables were included along with genotype to build predictive models for time to surgery.
Surgery occurred within 5 years in 239 subjects at a median time of 12 months. Three CD susceptibility loci were independently associated with surgery within 5 years (IL12B, IL23R, C11orf30). GWA identified novel putative loci associated with early surgery; 7q21 (CACNA2D1) and 9q34 (RXRA, COL5A1). The most predictive models of time to surgery included genetic and clinical risk factors. More than a 20% difference in frequency of progression to surgery was seen between the lowest and highest risk groups.
Progression to surgery is faster in CD patients with both genetic and clinical risk factors. IL12B is independently associated with need and time to early surgery in CD patients and justifies the investigation of novel and existing therapies that affect this pathway.
Keratoconus is a genetically heterogeneous corneal dystrophy. Previously, we performed two genome-wide linkage scans in a four generation autosomal dominant pedigree and repeatedly mapped a keratoconus locus to a genomic region located on chromosome 5q overlapping the gene encoding the inhibitor of calpains, calpastatin (CAST). To test whether variants in CAST gene are involved in genetic susceptibility to keratoconus we performed genetic testing of polymorphic markers in CAST gene in family and case-control panels of patients with keratoconus.
We genotyped SNPs (Single Nucleotide Polymorphisms) located in CAST gene in 262 patients in 40 Caucasian keratoconus families and in a Caucasian case-control panel with 304 cases and 518 controls. Generalized estimating equation models accounting for familial correlations implemented in GWAF program were used for association testing in families. Logistic regression models implemented in PLINK were performed to test associations in case-control samples.
Genetic testing of first set of seven SNPs in familial samples revealed two tentative nominally significant markers (rs4869307 p=0.03; rs27654: p=0.07). Additional genotyping of twelve tightly spaced SNPs identified CAST SNP rs4434401 to be associated with keratoconus in both familial and case-control panels with p values of 0.005 and 0.05, respectively; and with combined meta p value of familial and case-control cohorts of 0.002, or, after Bonferroni correction, 0.04.
Linkage analysis and genetic association support involvement of CAST gene in the genetic susceptibility to keratoconus. In-silico analysis of CAST expression suggests differential regulation of calpain/calpastatin system in cornea as a potential mechanism of functional defect.
Linkage analysis; genetic association; common variation; keratoconus; calpains; calpastatin
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 novel locus at 5q33 (GALNT10, rs7708584, p=3.4×10−11) and another at 7p15 when combined with 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 displayed directionally consistent effect estimates in our GWAS (binomial p=9.7×10−7), of which five 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.
Central corneal thickness (CCT) is a clinically important risk factor for primary open-angle glaucoma and keratoconus. Genetic factors controlling CCT in Latinos, the most populous minority population in the United States, are unclear. Here we describe the first genome-wide association study (GWAS) report of CCT in Latinos.
We performed a GWAS for CCT on 1768 Latinos recruited in the Los Angeles Latino Eye Study (LALES) using Illumina's HumanOmniExpress BeadChip (∼730K markers). To discover additional associated single-nucleotide polymorphisms (SNPs), we imputed SNPs based on the 1000 Genomes Project reference panels. All subjects were 40 years of age and older. We used linear regression with adjustment for age, sex, and principal components of genetic ancestry.
We replicated the involvement of several previously reported loci, such as RXRA-COL5A1, FOXO1, and ZNF469, for CCT in Latinos (P < 0.002). Moreover, we discovered novel SNPs, rs3118515, rs943423, rs3118594, and rs3132307, that reached GWAS significance (P < 5 × 10−8) in the uncharacterized LOC100506532 (gene type: miscRNA) for CCT in Latinos. By conditional analysis, we demonstrate that rs3118515 in this gene is responsible for the GWAS signal in the chromosome 9 RXRA-COL5A1 region in Latinos. Moreover, multiple sources of ENCODE evidence suggest that rs3118515 is in a regulatory region. Reverse-transcription PCR products indicated that transcripts of LOC100506532 surrounding rs3118515 were expressed in human corneas.
We discovered novel SNPs for CCT in Latinos and provided the first reported evidence of the corneal expression of LOC100506532. These results help to further increase our understanding of the genetic architecture of CCT.
In this first GWAS of central corneal thickness in Latinos, we discovered the involvement of novel SNPs, both directly genotyped and imputed, that reached genome-wide significance. This study provided the first reported evidence of the corneal expression of LOC100506532.
central corneal thickness; Latino; GWAS; rs3118515; LOC100506532
Single nucleotide polymorphisms (SNPs) located near or within the COL5A1 gene, at 9q34.2-q34.3 chromosomal region have been reported in association with central corneal thickness (CCT). Using family linkage analysis, we identified a keratoconus susceptibility locus at 9q34. These findings led us to perform an association study between COL5A1 variation and keratoconus susceptibility.
A Caucasian case–control cohort of 222 keratoconus patients and 3324 controls was selected as the discovery panel. An independent case–control panel of 304 cases and 518 controls and a family panel of 186 subjects were replicated for genotyping and association. Forty-four SNPs (21 for discovery and 23 for fine-mapping) spanning 300 kilobases in and around COL5A1 were genotyped and tested for genetic association. Logistic regression models implemented in PLINK were used to test for association in case controls. Generalized estimating equation models accounting for familial correlations implemented in genome-wide interaction analyses with family data were used for association testing in families.
Two CCT associated SNPs (rs1536482 and rs7044529 near and within COL5A1) were identified in the keratoconus discovery cohort (P values of 6.5 × 10−3 and 7.4 × 10−3). SNP rs1536482 was replicated in the second case–control sample (P = 0.02), and SNP rs7044529 was replicated in a keratoconus family panel (P = 0.03). Meta P values of rs1536482 and rs7044529 in the keratoconus cohorts were 1.5 × 10−4 (odds ratio [OR] = 1.30) and 2.9 × 10−3 (OR = 1.39). After Bonferroni correction, the association of SNP rs1536482 remained significant (P = 6.5 × 10−3).
SNPs in the COL5A1 region, which regulate normal variation in CCT, may play a role in the thinning associated with keratoconus.
Variants in the COL5A1 gene may contribute to genetic susceptibility to corneal thinning associated with keratoconus, in addition to their role in genetic regulation of normal variation in central corneal thickness.
keratoconus; association; COL5A1
Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10−6 were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10−8 for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10−8 for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10−8; RREB1: p = 5.7×10−8). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.
Central obesity is a marker of body fat distribution and is known to have a genetic underpinning. Few studies have reported genome-wide association study (GWAS) results among individuals of predominantly African ancestry (AA). We performed a collaborative meta-analysis in order to identify genetic loci associated with body fat distribution in AA individuals using waist circumference (WC) and waist to hip ratio (WHR) as measures of fat distribution, with and without adjustment for body mass index (BMI). We uncovered 2 genetic loci potentially associated with fat distribution: LHX2 in association with WC-adjusted-for-BMI and at RREB1 for WHR-adjusted-for-BMI. Six of fourteen previously reported loci for waist in EA populations were significant in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). These findings reinforce the concept that there are loci for body fat distribution that are independent of generalized adiposity.
Genotyping arrays are a cost effective approach when typing previously-identified genetic polymorphisms in large numbers of samples. One limitation of genotyping arrays with rare variants (e.g., minor allele frequency [MAF] <0.01) is the difficulty that automated clustering algorithms have to accurately detect and assign genotype calls. Combining intensity data from large numbers of samples may increase the ability to accurately call the genotypes of rare variants. Approximately 62,000 ethnically diverse samples from eleven Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium cohorts were genotyped with the Illumina HumanExome BeadChip across seven genotyping centers. The raw data files for the samples were assembled into a single project for joint calling. To assess the quality of the joint calling, concordance of genotypes in a subset of individuals having both exome chip and exome sequence data was analyzed. After exclusion of low performing SNPs on the exome chip and non-overlap of SNPs derived from sequence data, genotypes of 185,119 variants (11,356 were monomorphic) were compared in 530 individuals that had whole exome sequence data. A total of 98,113,070 pairs of genotypes were tested and 99.77% were concordant, 0.14% had missing data, and 0.09% were discordant. We report that joint calling allows the ability to accurately genotype rare variation using array technology when large sample sizes are available and best practices are followed. The cluster file from this experiment is available at www.chargeconsortium.com/main/exomechip.
Cerebral autosomal-dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a vascular dementia arising from abnormal arteriolar vascular smooth muscle cells. CADASIL results from mutations in Notch3 that alter the number of cysteine residues in the extracellular epidermal growth factor–like repeats, important for ligand binding. It is not known whether CADASIL mutations lead to loss or gain of Notch3 receptor function. To examine the functional consequences of CADASIL mutations, we engineered 4 CADASIL-like mutations into rat Notch3 and have shown that the presence of an unpaired cysteine does not impair cell-surface expression or ligand binding.
CADASIL; Notch3; vascular dementia; epidermal growth factor–like repeat; transmembrane receptor
Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10−8), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.
Men and women differ substantially regarding height, weight, and body fat. Interestingly, previous work detecting genetic effects for waist-to-hip ratio, to assess body fat distribution, has found that many of these showed sex-differences. However, systematic searches for sex-differences in genetic effects have not yet been conducted. Therefore, we undertook a genome-wide search for sexually dimorphic genetic effects for anthropometric traits including 133,723 individuals in a large meta-analysis and followed promising variants in further 137,052 individuals, including a total of 94 studies. We identified seven loci with significant sex-difference including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were significant in women, but not in men. Of interest is that sex-difference was only observed for waist phenotypes, but not for height or body-mass-index. We found no evidence for sex-differences with opposite effect direction for men and women. The PPARG locus is of specific interest due to its link to diabetes genetics and therapy. Our findings demonstrate the importance of investigating sex differences, which may lead to a better understanding of disease mechanisms with a potential relevance to treatment options.
Crohn’s disease (CD) and ulcerative colitis (UC), the two common forms of inflammatory bowel disease (IBD), affect over 2.5 million people of European ancestry with rising prevalence in other populations1. Genome-wide association studies (GWAS) and subsequent meta-analyses of CD and UC2,3 as separate phenotypes implicated previously unsuspected mechanisms, such as autophagy4, in pathogenesis and showed that some IBD loci are shared with other inflammatory diseases5. Here we expand knowledge of relevant pathways by undertaking a meta-analysis of CD and UC genome-wide association scans, with validation of significant findings in more than 75,000 cases and controls. We identify 71 new associations, for a total of 163 IBD loci that meet genome-wide significance thresholds. Most loci contribute to both phenotypes, and both directional and balancing selection effects are evident. Many IBD loci are also implicated in other immune-mediated disorders, most notably with ankylosing spondylitis and psoriasis. We also observe striking overlap between susceptibility loci for IBD and mycobacterial infection. Gene co-expression network analysis emphasizes this relationship, with pathways shared between host responses to mycobacteria and those predisposing to IBD.
There is evidence across several species for genetic control of phenotypic variation of complex traits1–4, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using 170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype)5–7, is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of 0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI8, possibly mediated by DNA methylation9,10. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.
Keratoconus is a condition in which the cornea progressively thins over time, and is a major cause for cornea transplantation. To identify keratoconus susceptibility regions, we performed a comprehensive genome-wide association study (GWAS) using a discovery and replication design. A discovery panel of 222 keratoconus Caucasian patients and 3324 Caucasian controls was genotyped using Illumina 370K beadchips. Further associated and fine-mapping single nucleotide polymorphisms (SNPs) (n= 4905) were genotyped in an independent replication case–control panel of 304 cases and 518 controls and a family panel of 307 subjects in 70 families. Logistic regression models implemented in PLINK were performed to test associations in case–control samples with and without principal component (PC) adjustments. Generalized estimation equation models accounting for familial correlations implemented in GWAF were used for association testing in families. No genome-wide associations were identified in the discovery GWAS panel. From the initial testing without adjustments for PCs, the top three SNPs located at 3p26 (rs6442925), 2q21.3 (rs4954218) and 19q13.3 (rs1428642) were identified with unadjusted P-values of 6.5 × 10−8, 2.4 × 10−7 and 3.1 × 10−7, respectively. After adjustments for PCs, rs1428642 became the most significant through the genome with a P-value of 1.4 × 10−6, while rs6442925 and rs4954218 were less significant (P= 1.9 × 10−5 and 2.6 × 10−4). SNP rs4954218 was confirmed in two independent replication panels with P-values of 0.004 and 0.009, respectively. Meta-analysis revealed a highest association at rs4954218 with adjusted P= 1.6 × 10−7 (unadjusted P= 1.2 × 10−9). These findings suggest SNP rs4954218, located near the RAB3GAP1 gene, previously reported to be associated with corneal malformation, is a potential susceptibility locus for keratoconus.
Genetic studies may help explain abnormalities of fat distribution in HIV-infected patients treated with antiretroviral therapy (ARV).
Subcutaneous adipose tissue (SAT) volume measured by magnetic resonance imaging (MRI) in leg, lower trunk, upper trunk, and arm was examined in 192 HIV-infected Caucasian men, ARV-treated from the Fat Redistribution and Metabolic Change in HIV infection (FRAM) study. Single nucleotide polymorphisms (SNPs) were assayed using the Illumina HumanCNV370-quad beadchip. Multivariate and univariate genome wide association analyses of the four SAT depots were implemented in PLINK software adjusted for age and ARV duration. Functional annotation analysis (FAA) using Ingenuity Systems Pathway Analysis tool (IPA) was carried out for markers with P<10-3 near known genes identified by multivariate analysis.
Loci (rs10504906, rs13267998, rs921231) in or near the anion exchanger solute carrier family 26, member 7 isoform a (SLC26A7) were strongly associated with upper trunk and arm SAT (9.8*10-7≤P<7.8*10-6). Loci (rs193139, rs7523050, rs1761621) in and near a gene rich region including G-protein-signaling modulator 2 (GPSM2) and syntaxin binding protein 3 (STXBP3) were significantly associated with lower body SAT depots (9.9*10-7≤P<9.5*10-6). GPSM2 is associated with cell division and cancer while STXBP3 is associated with glucose metabolism in adipoctyes. IPA identified atherosclerosis, mitochondrial function and T-Cell mediated apoptosis as processes related to SAT volume in HIV-infected individuals (P<5*10-3).
Our results are limited by the small sample size and replication is needed, however this genomic scan uncovered new genes associated with metabolism and inflammatory pathways that may affect SAT volume in ARV-treated HIV-infected patients.
HIV; HAART; GWAS; Subcutaneous Fat; SAT
Keratoconus is a bilateral noninflammatory progressive corneal disorder with complex genetic inheritance and a common cause for cornea transplantation in young adults. A genomewide linkage scan in keratoconus families identified a locus at 5q23.2, overlapping the gene coding for the lysyl oxidase (LOX). LOX encodes an enzyme responsible for collagen cross-linking in a variety of tissues including the cornea. Corneal collagen cross-linking with long-wave ultraviolet light and riboflavin is a promising new treatment for keratoconus. To determine whether LOX is a genetic determinant of the pathogenesis of keratoconus, we analyzed association results of LOX polymorphisms in two independent case-control samples and in keratoconus families.
Association results were analyzed of single-nucleotide polymorphisms (SNPs) in the LOX gene from a Genome-Wide Association Study (GWAS) investigation in two independent panels of patients with keratoconus and controls and in keratoconus families.
Evidence of association was found at SNPs rs10519694 and rs2956540 located in intron 4 of LOX in the GWAS discovery case-control panel with P values of 2.3 × 10−3 and 7 × 10−3, respectively. The same two SNPs were found to be associated with keratoconus by family-based association testing with P values of 2.7 × 10−3 and 7.7 × 10−4, respectively. Meta P values of 4.0 × 10−5 and 4.0 × 10−7 were calculated for SNPs rs10519694 and rs2956540 by analyzing case-control and family samples simultaneously. Sequencing of LOX exons in a subset of keratoconus patients identified two polymorphisms, rs1800449 and rs2288393, located in LOX transcripts I and II, associated with keratoconus in case-control and family samples with a meta P value of 0.02.
Results provided strong genetic evidence that LOX variants lead to increased susceptibility to developing of keratoconus.
Keratoconus is a progressive disorder of the cornea with complex genetic inheritance, which can lead to severe visual impairment or blindness. Here we describe the association between variation in the lysyl oxidase gene and keratoconus susceptibility in case-control and family studies.
Genome-wide association studies (GWAS) have identified multiple Crohn’s disease (CD) susceptibility loci, including association with non-coding intergenic single nucleotide polymorphisms (SNPs) at 10q21.
To fine-map the 10q21 locus, we genotyped 86 SNPs in 1632 CD cases and 961 controls and performed single marker and conditional analyses using logistic regression.
We observed association with CD risk spanning eleven SNPs (p< 0.001). The most significant association observed was at the nonsynonymous SNP rs7076156 (Ala62Thr) in ZNF365. The alanine allele was over-represented in CD (p= 5.23×10−7; OR= 1.39 [1.22–1.58]); allele frequency 76% CD, 69.7% controls). Conditional analysis on rs7076156 nullified all other significant associations, suggesting that this is the causative variant at this locus. Four isoforms of ZNF365 have previously been identified and rs7076156 is located in an exon unique to ZNF365 isoform D. We demonstrated, using RT-PCR, expression of ZNF365D in intestinal resections from both CD subjects and controls. We identified markedly reduced ZNF365D mean expression levels in EBV-transformed lymphoblastoid cell lines (ELCLs) from CD subjects homozygous for the risk allele (Ala). A whole-genome microarray expression study further suggested that the Ala62Thr change in ZNF365 isoform D is related to differential expression of the genes ARL4A, MKKS, RRAGD, SUMF2, TDR1 and ZNF148 in CD.
Collectively our data support the hypothesis that the nonsynonymous Ala62Thr SNP rs7076156 underlies the association between 10q21 and CD risk and suggests that this SNP acts by altering expression of genes under the control of ZNF365 isoform D.
GWAS; Crohn’s Disease; ZNF365; genetics; association study
Genotype imputation is a vital tool in genome-wide association studies (GWAS) and meta-analyses of multiple GWAS results. Imputation enables researchers to increase genomic coverage and to pool data generated using different genotyping platforms. HapMap samples are often employed as the reference panel. More recently, the 1000 Genomes Project resource is becoming the primary source for reference panels. Multiple GWAS and meta-analyses are targeting Latinos, the most populous, and fastest growing minority group in the US. However, genotype imputation resources for Latinos are rather limited compared to individuals of European ancestry at present, largely because of the lack of good reference data. One choice of reference panel for Latinos is one derived from the population of Mexican individuals in Los Angeles contained in the HapMap Phase 3 project and the 1000 Genomes Project. However, a detailed evaluation of the quality of the imputed genotypes derived from the public reference panels has not yet been reported. Using simulation studies, the Illumina OmniExpress GWAS data from the Los Angles Latino Eye Study and the MACH software package, we evaluated the accuracy of genotype imputation in Latinos. Our results show that the 1000 Genomes Project AMR + CEU + YRI reference panel provides the highest imputation accuracy for Latinos, and that also including Asian samples in the panel can reduce imputation accuracy. We also provide the imputation accuracy for each autosomal chromosome using the 1000 Genomes Project panel for Latinos. Our results serve as a guide to future imputation based analysis in Latinos.
genotype imputation; Latino; HapMap Project; 1000 Genomes Project
More than a thousand disease susceptibility loci have been identified via genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings generally remain to be defined. We utilize pooled next-generation sequencing to study 56 genes in regions associated to Crohn’s Disease in 350 cases and 350 controls. Follow up genotyping of 70 rare and low-frequency protein-altering variants (MAF ~ .001-.05) in nine independent case-control series (16054 CD patients, 12153 UC patients, 17575 healthy controls) identifies four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association to a novel, protective splice variant in CARD9 (p < 1e-16, OR ~ 0.29), as well as additional associations to coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by providing novel, rare, and likely functional variants that will empower functional experiments and predictive models.
Elevated serum urate levels can lead to gout and are associated with cardiovascular risk factors. We performed genome-wide association to search for genetic susceptibility loci for serum urate and gout, and investigated the causal nature of the associations of serum urate with gout and selected cardiovascular risk factors and coronary heart disease (CHD).
Methods and Results
Meta-analyses of genome-wide association studies (GWAS) were performed in 5 population-based cohorts of the CHARGE consortium for serum urate and gout in 28,283 white individuals. The effect of the most significant SNP at all genome-wide significant loci on serum urate was added to create a genetic urate score. Findings were replicated in the Women’s Genome Health Study (WGHS; n=22,054). SNPs at 8 genetic loci achieved genome-wide significance with serum urate levels (p-values 4×10−8 to 2×10−242; SLC22A11, GCKR, R3HDM2-INHBC region, RREB1, PDZK1, SLC2A9, ABCG2, SLC17A1). Only two loci [SLC2A9, ABCG2] showed genome-wide significant association with gout. The genetic urate score was strongly associated with serum urate and gout (odds ratio 12.4 per 100 umol/L; p-value=3×10−39), but not with blood pressure, glucose, eGFR, chronic kidney disease, or CHD. The lack of association between the genetic score and the latter phenotypes was also observed in WGHS.
The genetic urate score analysis suggested a causal relationship between serum urate and gout but did not provide evidence for one between serum urate and cardiovascular risk factors and CHD.
urate; gout; cardiovascular disease risk factors; genome-wide association study; Mendelian randomization
This is a meta-analysis of two genome-wide association studies that found evidence of association of keratoconus with polymorphisms in the promoter of the HGF gene. One polymorphism is associated with higher levels of serum HGF.
Keratoconus is a progressive disorder of the cornea that can lead to severe visual impairment or blindness. Although several genomic regions have been linked to rare familial forms of keratoconus, no genes have yet been definitively identified for common forms of the disease.
Two genome-wide association scans were undertaken in parallel. The first used pooled DNA from an Australian cohort, followed by typing of top-ranked single-nucleotide polymorphisms (SNPs) in individual DNA samples. The second was conducted in individually genotyped patients, and controls from the USA. Tag SNPs around the hepatocyte growth factor (HGF) gene were typed in three additional replication cohorts. Serum levels of HGF protein in normal individuals were assessed with ELISA and correlated with genotype.
The only SNP observed to be associated in both the pooled discovery and primary replication cohort was rs1014091, located upstream of the HGF gene. The nearby SNP rs3735520 was found to be associated in the individually typed discovery cohort (P = 6.1 × 10−7). Genotyping of tag SNPs around HGF revealed association at rs3735520 and rs17501108/rs1014091 in four of the five cohorts. Meta-analysis of all five datasets together yielded suggestive P values for rs3735520 (P = 9.9 × 10−7) and rs17501108 (P = 9.9 × 10−5). In addition, SNP rs3735520 was found to be associated with serum HGF level in normal individuals (P = 0.036).
Taken together, these results implicate genetic variation at the HGF locus with keratoconus susceptibility.
The metabolic syndrome (MetS) is defined as concomitant disorders of lipid and glucose metabolism, central obesity, and high blood pressure, with an increased risk of type 2 diabetes and cardiovascular disease. This study tests whether common genetic variants with pleiotropic effects account for some of the correlated architecture among five metabolic phenotypes that define MetS.
RESEARCH DESIGN AND METHODS
Seven studies of the STAMPEED consortium, comprising 22,161 participants of European ancestry, underwent genome-wide association analyses of metabolic traits using a panel of ∼2.5 million imputed single nucleotide polymorphisms (SNPs). Phenotypes were defined by the National Cholesterol Education Program (NCEP) criteria for MetS in pairwise combinations. Individuals exceeding the NCEP thresholds for both traits of a pair were considered affected.
Twenty-nine common variants were associated with MetS or a pair of traits. Variants in the genes LPL, CETP, APOA5 (and its cluster), GCKR (and its cluster), LIPC, TRIB1, LOC100128354/MTNR1B, ABCB11, and LOC100129150 were further tested for their association with individual qualitative and quantitative traits. None of the 16 top SNPs (one per gene) associated simultaneously with more than two individual traits. Of them 11 variants showed nominal associations with MetS per se. The effects of 16 top SNPs on the quantitative traits were relatively small, together explaining from ∼9% of the variance in triglycerides, 5.8% of high-density lipoprotein cholesterol, 3.6% of fasting glucose, and 1.4% of systolic blood pressure.
Qualitative and quantitative pleiotropic tests on pairs of traits indicate that a small portion of the covariation in these traits can be explained by the reported common genetic variants.