Beth L. Cobb

Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM , PXK , KIAA1542 and other loci

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.

Identification of novel genetic susceptibility loci for Behçet's disease using a genome-wide association study

IntroductionBehçets disease is a chronic systemic inflammatory disease that remains incompletely understood. Herein, we perform the first genome-wide association study in Behçets disease.MethodsUsing DNA pooling technology and the Affymetrix 500K arrays, we identified possible candidate gene associations with Behçets disease in a cohort of 152 Behçets disease patients and 172 healthy ethnically matched controls. Genetic loci that were identified in the pooling study were genotyped in patients and controls using TaqMan genotyping technology.ResultsWe identified genetic associations between Behçets disease and single-nucleotide polymorphisms (SNPs) in KIAA1529, CPVL, LOC100129342, UBASH3B, and UBAC2 (odds ratio = 2.04, 2.26, 1.84, 1.71, and 1.61, respectively; P value = 4.2 × 10-5, 1.0 × 10-4, 3.0 × 10-4, 1.5 × 10-3, and 5.8 × 10-3, respectively). Among the associated SNPs, the Behçets disease-risk allele in rs2061634 leads to substitution of serine to cysteine at amino acid position 995 (S995C) in the KIAA1529 protein.ConclusionsUsing an unbiased whole-genome genetic association approach, we identified novel candidate genetic loci that are associated with increased susceptibility for Behçets disease. These findings will help to better understand the pathogenesis of Behçets disease and identify novel targets for therapeutic intervention.

Genes and Sjögren's syndrome.

The evidence for a strong genetic component conferring susceptibility to primary Sjögrens syndrome (SS) is mounting. Several associations with SS have been reported and provide evidence that the HLA region harbors important susceptibility loci and that multiple genes outside the HLA region play a role. Genetic discovery lags behind success observed in related autoimmune diseases. Identifying genetic factors that cause SS will allow more precise definition of pathogenic mechanisms leading to the overall SS phenotype and clinically heterogeneous subsets of patients. Critical opportunities are certain to follow for translation into improved diagnosis and therapies for SS and its spectrum diseases.

Phenome-wide association study (PheWAS) in EMR-linked pediatric cohorts, genetically links PLCL1 to speech language development and IL5-IL13 to Eosinophilic Esophagitis.

Objective: We report the first pediatric specific Phenome-Wide Association Study (PheWAS) using electronic medical records (EMRs). Given the early success of PheWAS in adult populations, we investigated the feasibility of this approach in pediatric cohorts in which associations between a previously known genetic variant and a wide range of clinical or physiological traits were evaluated. Although computationally intensive, this approach has potential to reveal disease mechanistic relationships between a variant and a network of phenotypes. Method: Data on 5049 samples of European ancestry were obtained from the EMRs of two large academic centers in five different genotyped cohorts. Recently, these samples have undergone whole genome imputation. After standard quality controls, removing missing data and outliers based on principal components analyses (PCA), 4268 samples were used for the PheWAS study. We scanned for associations between 2476 single-nucleotide polymorphisms (SNP) with available genotyping data from previously published GWAS studies and 539 EMR-derived phenotypes. The false discovery rate was calculated and, for any new PheWAS findings, a permutation approach (with up to 1,000,000 trials) was implemented. Results: This PheWAS found a variety of common variants (MAF > 10%) with prior GWAS associations in our pediatric cohorts including Juvenile Rheumatoid Arthritis (JRA), Asthma, Autism and Pervasive Developmental Disorder (PDD) and Type 1 Diabetes with a false discovery rate < 0.05 and power of study above 80%. In addition, several new PheWAS findings were identified including a cluster of association near the NDFIP1 gene for mental retardation (best SNP rs10057309, p = 4.33 × 10−7, OR = 1.70, 95%CI = 1.38 − 2.09); association near PLCL1 gene for developmental delays and speech disorder [best SNP rs1595825, p = 1.13 × 10−8, OR = 0.65(0.57 − 0.76)]; a cluster of associations in the IL5-IL13 region with Eosinophilic Esophagitis (EoE) [best at rs12653750, p = 3.03 × 10−9, OR = 1.73 95%CI = (1.44 − 2.07)], previously implicated in asthma, allergy, and eosinophilia; and association of variants in GCKR and JAZF1 with allergic rhinitis in our pediatric cohorts [best SNP rs780093, p = 2.18 × 10−5, OR = 1.39, 95%CI = (1.19 − 1.61)], previously demonstrated in metabolic disease and diabetes in adults. Conclusion: The PheWAS approach with re-mapping ICD-9 structured codes for our European-origin pediatric cohorts, as with the previous adult studies, finds many previously reported associations as well as presents the discovery of associations with potentially important clinical implications.

Two Functional Lupus-Associated BLK Promoter Variants Control Cell-Type- and Developmental-Stage-Specific Transcription

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

Genetic risk for malignant hyperthermia in non-anesthesia-induced myopathies.

Malignant hyperthermia (MH) is a pharmacogenetic, autosomal dominantly inherited disorder of skeletal muscle triggered by volatile anesthetics and infrequently by extreme exertion and heat exposure. MH has variable penetrance with an incidence ranging from 1 in 5000 to 1 in 50,000-100,000 anesthesias. Mutations in the ryanodine receptor gene, RYR1, are found in 50-70% of cases. We hypothesized that a portion of patients with drug-induced muscle diseases, unrelated to anesthesia, such as severe statin myopathy, have underlying genetic liability that may include RYR1 gene mutations. DNA samples were collected from 885 patients in 4 groups: severe statin myopathy (n=197), mild statin myopathy (n=163), statin-tolerant controls (n=133), and non-drug-induced myopathies of unknown etiology characterized by exercise-induced muscle pain and weakness (n=392). Samples were screened for 105 mutations and variants in 26 genes associated with 7 categories of muscle disease including 34 mutations and variants in the RYR1 gene. Disease-causing mutations or variants in RYR1 were present in 3 severe statin myopathy cases, 1 mild statin myopathy case, 8 patients with non-drug-induced myopathy, and none in controls. These results suggest that disease-causing mutations and certain variants in the RYR1 gene may contribute to underlying genetic risk for non-anesthesia-induced myopathies and should be included in genetic susceptibility screening in patients with severe statin myopathy and in patients with non-statin-induced myopathies of unknown etiology.

EMR-linked GWAS study: investigation of variation landscape of loci for body mass index in children

Common variations at the loci harboring the fat mass and obesity gene (FTO), MC4R, and TMEM18 are consistently reported as being associated with obesity and body mass index (BMI) especially in adult population. In order to confirm this effect in pediatric population five European ancestry cohorts from pediatric eMERGE-II network (CCHMC-BCH) were evaluated. Method: Data on 5049 samples of European ancestry were obtained from the Electronic Medical Records (EMRs) of two large academic centers in five different genotyped cohorts. For all available samples, gender, age, height, and weight were collected and BMI was calculated. To account for age and sex differences in BMI, BMI z-scores were generated using 2000 Centers of Disease Control and Prevention (CDC) growth charts. A Genome-wide association study (GWAS) was performed with BMI z-score. After removing missing data and outliers based on principal components (PC) analyses, 2860 samples were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and BMI was tested using linear regression adjusting for age, gender, and PC by cohort. The effects of SNPs were modeled assuming additive, recessive, and dominant effects of the minor allele. Meta-analysis was conducted using a weighted z-score approach. Results: The mean age of subjects was 9.8 years (range 2–19). The proportion of male subjects was 56%. In these cohorts, 14% of samples had a BMI ≥95 and 28 ≥ 85%. Meta analyses produced a signal at 16q12 genomic region with the best result of p = 1.43 × 10-7 [p(rec) = 7.34 × 10-8) for the SNP rs8050136 at the first intron of FTO gene (z = 5.26) and with no heterogeneity between cohorts (p = 0.77). Under a recessive model, another published SNP at this locus, rs1421085, generates the best result [z = 5.782, p(rec) = 8.21 × 10-9]. Imputation in this region using dense 1000-Genome and Hapmap CEU samples revealed 71 SNPs with p < 10-6, all at the first intron of FTO locus. When hetero-geneity was permitted between cohorts, signals were also obtained in other previously identified loci, including MC4R (rs12964056, p = 6.87 × 10-7, z = -4.98), cholecystokinin CCK (rs8192472, p = 1.33 × 10-6, z = -4.85), Interleukin 15 (rs2099884, p = 1.27 × 10-5, z = 4.34), low density lipoprotein receptor-related protein 1B [LRP1B (rs7583748, p = 0.00013, z = -3.81)] and near transmembrane protein 18 (TMEM18) (rs7561317, p = 0.001, z = -3.17). We also detected a novel locus at chromosome 3 at COL6A5 [best SNP = rs1542829, minor allele frequency (MAF) of 5% p = 4.35 × 10-9, z = 5.89]. Conclusion: An EMR linked cohort study demonstrates that the BMI-Z measurements can be successfully extracted and linked to genomic data with meaningful confirmatory results. We verified the high prevalence of childhood rate of overweight and obesity in our cohort (28%). In addition, our data indicate that genetic variants in the first intron of FTO, a known adult genetic risk factor for BMI, are also robustly associated with BMI in pediatric population.

Association of common variants in the human eyes shut ortholog (EYS) with statin-induced myopathy: evidence for additional functions of EYS.

Introduction: Of the nearly 38 million people in the USA who receive statin therapy, 0.1‐0.5% experience severe or life‐threatening myopathic side effects. Methods: We performed a genome‐wide association study (GWAS) in a group of patients with severe statin myopathy versus a statin‐tolerant group to identify genetic susceptibility loci. Results: Replication studies in independent groups of severe statin myopathy (n = 190) and statin‐tolerant controls (n = 130) resulted in the identification of three single‐nucleotide polymorphisms (SNPs), rs9342288, rs1337512, and rs3857532, in the eyes shut homolog (EYS) on chromosome 6 suggestive of an association with risk for severe statin myopathy (P = 0.0003–0.0008). Analysis of EYS cDNA demonstrated that EYS gene products are complex and expressed with relative abundance in the spinal cord as well as in the retina. Conclusion: Structural similarities of these EYS gene products to members of the Notch signaling pathway and to agrin suggest a possible functional role in the maintenance and regeneration of the structural integrity of skeletal muscle. Muscle Nerve, 2011

Genetic association between methyl-CpG binding protein 2 (MECP2) and primary Sjögren's syndrome

Primary Sjogrens syndrome (pSS) is a complex autoimmune disease that is incompletely understood. The principal manifestations of pSS are dry eyes and dry mouth due to lymphocytic infiltration of the salivary and lacrimal glands. Like lupus, women are affected by pSS over nine times more than with men. The association of human leucocyte antigen loci with pSS and other autoimmune diseases is well documented.1,–,3 A genetic association of the methyl-CpG binding protein 2 ( MECP2 ) gene with lupus has been reported and confirmed in multiple independent cohorts.4 5 Herein, we examine the genetic association between MECP2 and pSS in a large cohort of pSS patients and healthy controls. We studied a cohort of 460 European-derived independent pSS patients (423 women and 37 men) and 1828 ethnically matched normal healthy controls (1279 women and 549 men). All patients fulfiled the American–European Consensus Group classification criteria for pSS.6 Our study protocols were …

Practical considerations in genomic decision support: The eMERGE experience

Background: Genomic medicine has the potential to improve care by tailoring treatments to the individual. There is consensus in the literature that pharmacogenomics (PGx) may be an ideal starting point for real-world implementation, due to the presence of well-characterized drug-gene interactions. Clinical Decision Support (CDS) is an ideal avenue by which to implement PGx at the bedside. Previous literature has established theoretical models for PGx CDS implementation and discussed a number of anticipated real-world challenges. However, work detailing actual PGx CDS implementation experiences has been limited. Anticipated challenges include data storage and management, system integration, physician acceptance, and more. Methods: In this study, we analyzed the experiences of ten members of the Electronic Medical Records and Genomics (eMERGE) Network, and one affiliate, in their attempts to implement PGx CDS. We examined the resulting PGx CDS system characteristics and conducted a survey to understand the unanticipated implementation challenges sites encountered. Results: Ten sites have successfully implemented at least one PGx CDS rule in the clinical setting. The majority of sites elected to create an Omic Ancillary System (OAS) to manage genetic and genomic data. All sites were able to adapt their existing CDS tools for PGx knowledge. The most common and impactful delays were not PGx-specific issues. Instead, they were general IT implementation problems, with top challenges including team coordination/communication and staffing. The challenges encountered caused a median total delay in system go-live of approximately two months. Conclusions: These results suggest that barriers to PGx CDS implementations are generally surmountable. Moreover, PGx CDS implementation may not be any more difficult than other healthcare IT projects of similar scope, as the most significant delays encountered were not unique to genomic medicine. These are encouraging results for any institution considering implementing a PGx CDS tool, and for the advancement of genomic medicine.