Medical Records-Based Genetic Studies of the Complement System.

BACKGROUND: Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts. METHODS: We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network. RESULTS: In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; ß=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; ß=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; ß=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (ß=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (ß=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation. CONCLUSIONS: We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype.

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome-wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single-nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA-B herpes zoster (shingles) association and discovered a novel zoster-associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

Identification of Four Novel Loci in Asthma in European American and African American Populations.

RATIONALE: Despite significant advances in knowledge of the genetic architecture of asthma, specific contributors to the variability in the burden between populations remain uncovered. OBJECTIVES: To identify additional genetic susceptibility factors of asthma in European American and African American populations. METHODS: A phenotyping algorithm mining electronic medical records was developed and validated to recruit cases with asthma and control subjects from the Electronic Medical Records and Genomics network. Genome-wide association analyses were performed in pediatric and adult asthma cases and control subjects with European American and African American ancestry followed by metaanalysis. Nominally significant results were reanalyzed conditioning on allergy status. MEASUREMENTS AND MAIN RESULTS: The validation of the algorithm yielded an average of 95.8% positive predictive values for both cases and control subjects. The algorithm accrued 21,644 subjects (65.83% European American and 34.17% African American). We identified four novel population-specific associations with asthma after metaanalyses: loci 6p21.31, 9p21.2, and 10q21.3 in the European American population, and the PTGES gene in African Americans. TEK at 9p21.2, which encodes TIE2, has been shown to be involved in remodeling the airway wall in asthma, and the association remained significant after conditioning by allergy. PTGES, which encodes the prostaglandin E synthase, has also been linked to asthma, where deficient prostaglandin E2 synthesis has been associated with airway remodeling. CONCLUSIONS: This study adds to understanding of the genetic architecture of asthma in European Americans and African Americans and reinforces the need to study populations of diverse ethnic backgrounds to identify shared and unique genetic predictors of asthma.

Genetic regulation of serum IgA levels and susceptibility to common immune, infectious, kidney, and cardio-metabolic traits.

Immunoglobulin A (IgA) mediates mucosal responses to food antigens and the intestinal microbiome and is involved in susceptibility to mucosal pathogens, celiac disease, inflammatory bowel disease, and IgA nephropathy. We performed a genome-wide association study of serum IgA levels in 41,263 individuals of diverse ancestries and identified 20 genome-wide significant loci, including 9 known and 11 novel loci. Co-localization analyses with expression QTLs prioritized candidate genes for 14 of 20 significant loci. Most loci encoded genes that produced immune defects and IgA abnormalities when genetically manipulated in mice. We also observed positive genetic correlations of serum IgA levels with IgA nephropathy, type 2 diabetes, and body mass index, and negative correlations with celiac disease, inflammatory bowel disease, and several infections. Mendelian randomization supported elevated serum IgA as a causal factor in IgA nephropathy. African ancestry was consistently associated with higher serum IgA levels and greater frequency of IgA-increasing alleles compared to other ancestries. Our findings provide novel insights into the genetic regulation of IgA levels and its potential role in human disease.

Quantifying the phenome-wide disease burden of obesity using electronic health records and genomics.

OBJECTIVE: High BMI is associated with many comorbidities and mortality. This study aimed to elucidate the overall clinical risk of obesity using a genome- and phenome-wide approach. METHODS: This study performed a phenome-wide association study of BMI using a clinical cohort of 736,726 adults. This was followed by genetic association studies using two separate cohorts: one consisting of 65,174 adults in the Electronic Medical Records and Genomics (eMERGE) Network and another with 405,432 participants in the UK Biobank. RESULTS: Class 3 obesity was associated with 433 phenotypes, representing 59.3% of all billing codes in individuals with severe obesity. A genome-wide polygenic risk score for BMI, accounting for 7.5% of variance in BMI, was associated with 296 clinical diseases, including strong associations with type 2 diabetes, sleep apnea, hypertension, and chronic liver disease. In all three cohorts, 199 phenotypes were associated with class 3 obesity and polygenic risk for obesity, including novel associations such as increased risk of renal failure, venous insufficiency, and gastroesophageal reflux. CONCLUSIONS: This combined genomic and phenomic systematic approach demonstrated that obesity has a strong genetic predisposition and is associated with a considerable burden of disease across all disease classes.

A Mendelian Randomization Approach Using 3-HMG-Coenzyme-A Reductase Gene Variation to Evaluate the Association of Statin-Induced Low-Density Lipoprotein Cholesterol Lowering With Noncardiovascular Disease Phenotypes.

Importance: Observational studies suggest that statins, which inhibit 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase, may be associated with beneficial effects in many noncardiovascular diseases. Objective: To construct a weighted HMG-CoA reductase (HMGCR) gene genetic risk score (GRS) using variants in the HMGCR gene affecting low-density lipoprotein cholesterol as an instrumental variable for mendelian randomization analyses to test associations with candidate noncardiovascular phenotypes previously associated with statin use in observational studies. Design, Setting, and Participants: This cohort study included 53 385 unrelated adults of European ancestry with genome-wide genotypes available from BioVU (a practice-based biobank, used for discovery) and 30 444 unrelated adults with European ancestry available in the Electronic Medical Records and Genomics (eMERGE; a research consortium that conducts genetic research using electronic medical records, used for replication). The study was conducted from February 6, 2015, through April 31, 2019; data analysis was performed from August 26, 2019, through December 22, 2020. Interventions: An HMGCR GRS was calculated. Main Outcomes and Measures: The association between the HMGCR GRS and the presence or absence of 22 noncardiovascular phenotypes previously associated with statin use in clinical studies. Results: Of the 53 385 individuals in BioVU, 29 958 (56.1%) were women; mean (SD) age was 59.9 (15.6) years. The finding between the HMGCR GRS and the noncardiovascular phenotypes of interest in this cohort was significant only for type 2 diabetes. An HMGCR GRS equivalent to a 10-mg/dL decrease in the low-density lipoprotein cholesterol level was associated with an increased risk of type 2 diabetes (odds ratio [OR], 1.09; 95% CI, 1.04-1.15; P = 5.58 × 10-4). The HMGCR GRS was not associated with other phenotypes; the closest were increased risk of Parkinson disease (OR, 1.30; 95% CI, 1.07-1.58; P = .007) and kidney failure (OR, 1.18; 95% CI, 1.05-1.34; P = .008). Of the 30 444 individuals in eMERGE, 16 736 (55.0%) were women; mean (SD) age was 68.7 (15.4) years. The association between the HMGCR GRS and type 2 diabetes was replicated in this cohort (OR, 1.09; 95% CI, 1.01-1.17; P = .02); however, the HMGCR GRS was not associated with Parkinson disease (OR, 0.93; 95% CI, 0.75-1.16; P = .53) and kidney failure (OR, 1.18; 95% CI, 0.98-1.41; P = .08) in the eMERGE cohort. Conclusions and Relevance: A mendelian randomization approach using variants in the HMGCR gene replicated the association between statin use and increased type 2 diabetes risk but provided no strong evidence for pleiotropic effects of statin-induced decrease of the low-density lipoprotein cholesterol level on other diseases.

Patients with Proliferative Lupus Nephritis Have Autoantibodies That React to Moesin and Demonstrate Increased Glomerular Moesin Expression.

Kidney involvement in systemic lupus erythematosus (SLE)-termed lupus nephritis (LN)-is a severe manifestation of SLE that can lead to end-stage kidney disease (ESKD). LN is characterized by immune complex deposition and inflammation in the glomerulus. We tested the hypothesis that autoantibodies targeting podocyte and glomerular cell proteins contribute to the development of immune complex formation in LN. We used Western blotting with SLE sera from patients with and without LN to identify target antigens in human glomerular and cultured human-derived podocyte membrane proteins. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), we identified the proteins in the gel regions corresponding to reactive bands observed with sera from LN patients. We identified 102 proteins that were present in both the podocyte and glomerular samples. We identified 10 high-probability candidates, including moesin, using bioinformatic analysis. Confirmation of moesin as a target antigen was conducted using immunohistochemical analysis (IHC) of kidney biopsy tissue and enzyme-linked immunosorbent assay (ELISA) to detect circulating antibodies. By IHC, biopsies from patients with proliferative lupus nephritis (PLN, class III/IV) demonstrated significantly increased glomerular expression of moesin (p < 0.01). By ELISA, patients with proliferative LN demonstrated significantly increased antibodies against moesin (p < 0.01). This suggests that moesin is a target glomerular antigen in lupus nephritis.

Associations of Genetically Predicted Lp(a) (Lipoprotein [a]) Levels With Cardiovascular Traits in Individuals of European and African Ancestry.

BACKGROUND: Lp(a) (lipoprotein [a]) levels are higher in individuals of African ancestry (AA) than in individuals of European ancestry (EA). We examined associations of genetically predicted Lp(a) levels with (1) atherosclerotic cardiovascular disease subtypes: coronary heart disease, cerebrovascular disease, peripheral artery disease, and abdominal aortic aneurysm and (2) nonatherosclerotic cardiovascular disease phenotypes, stratified by ancestry. METHODS: We performed (1) Mendelian randomization analyses for previously reported cardiovascular associations and (2) Mendelian randomization-phenome-wide association analyses for novel associations. Analyses were stratified by ancestry in electronic Medical Records and Genomics, United Kingdom Biobank, and Million Veteran Program cohorts separately and in a combined cohort of 804 507 EA and 103 580 AA participants. RESULTS: In Mendelian randomization analyses using the combined cohort, a 1-SD genetic increase in Lp(a) level was associated with atherosclerotic cardiovascular disease subtypes in EA-odds ratio and 95% CI for coronary heart disease 1.28 (1.16-1.41); cerebrovascular disease 1.14 (1.07-1.21); peripheral artery disease 1.22 (1.11-1.34); abdominal aortic aneurysm 1.28 (1.17-1.40); in AA, the effect estimate was lower than in EA and nonsignificant for coronary heart disease 1.11 (0.99-1.24) and cerebrovascular disease 1.06 (0.99-1.14) but similar for peripheral artery disease 1.16 (1.01-1.33) and abdominal aortic aneurysm 1.34 (1.11-1.62). In EA, a 1-SD genetic increase in Lp(a) level was associated with aortic valve disorders 1.34 (1.10-1.62), mitral valve disorders 1.18 (1.09-1.27), congestive heart failure 1.12 (1.05-1.19), and chronic kidney disease 1.07 (1.01-1.14). In AA, no significant associations were noted for aortic valve disorders 1.08 (0.94-1.25), mitral valve disorders 1.02 (0.89-1.16), congestive heart failure 1.02 (0.95-1.10), or chronic kidney disease 1.05 (0.99-1.12). Mendelian randomization-phenome-wide association analyses identified novel associations in EA with arterial thromboembolic disease, nonaortic aneurysmal disease, atrial fibrillation, cardiac conduction disorders, and hypertension. CONCLUSIONS: Many cardiovascular associations of genetically increased Lp(a) that were significant in EA were not significant in AA. Lp(a) was associated with atherosclerotic cardiovascular disease in four major arterial beds in EA but only with peripheral artery disease and abdominal aortic aneurysm in AA. Additionally, novel cardiovascular associations were detected in EA.

Ethical Considerations Related to Return of Results from Genomic Medicine Projects: The eMERGE Network (Phase III) Experience.

We examined the Institutional Review Board (IRB) process at 9 academic institutions in the electronic Medical Records and Genomics (eMERGE) Network, for proposed electronic health record-based genomic medicine studies, to identify common questions and concerns. Sequencing of 109 disease related genes and genotyping of 14 actionable variants is being performed in ~28,100 participants from the 9 sites. Pathogenic/likely pathogenic variants in actionable genes are being returned to study participants. We examined each site's research protocols, informed-consent materials, and interactions with IRB staff. Research staff at each site completed questionnaires regarding their IRB interactions. The time to prepare protocols for IRB submission, number of revisions and time to approval ranged from 10-261 days, 0-11, and 11-90 days, respectively. IRB recommendations related to the readability of informed consent materials, specifying the full range of potential risks, providing options for receiving limited results or withdrawal, sharing of information with family members, and establishing the mechanisms to answer participant questions. IRBs reviewing studies that involve the return of results from genomic sequencing have a diverse array of concerns, and anticipating these concerns can help investigators to more effectively engage IRBs.

Suboptimal Clinical Documentation in Young Children with Severe Obesity at Tertiary Care Centers.

Background and Objectives. The prevalence of severe obesity in children has doubled in the past decade. The objective of this study is to identify the clinical documentation of obesity in young children with a BMI ≥ 99th percentile at two large tertiary care pediatric hospitals. Methods. We used a standardized algorithm utilizing data from electronic health records to identify children with severe early onset obesity (BMI ≥ 99th percentile at age <6 years). We extracted descriptive terms and ICD-9 codes to evaluate documentation of obesity at Boston Children's Hospital and Cincinnati Children's Hospital and Medical Center between 2007 and 2014. Results. A total of 9887 visit records of 2588 children with severe early onset obesity were identified. Based on predefined criteria for documentation of obesity, 21.5% of children (13.5% of visits) had positive documentation, which varied by institution. Documentation in children first seen under 2 years of age was lower than in older children (15% versus 26%). Documentation was significantly higher in girls (29% versus 17%, p < 0.001), African American children (27% versus 19% in whites, p < 0.001), and the obesity focused specialty clinics (70% versus 15% in primary care and 9% in other subspecialty clinics, p < 0.001). Conclusions. There is significant opportunity for improvement in documentation of obesity in young children, even years after the 2007 AAP guidelines for management of obesity.

Imputation and quality control steps for combining multiple genome-wide datasets.

The electronic MEdical Records and GEnomics (eMERGE) network brings together DNA biobanks linked to electronic health records (EHRs) from multiple institutions. Approximately 51,000 DNA samples from distinct individuals have been genotyped using genome-wide SNP arrays across the nine sites of the network. The eMERGE Coordinating Center and the Genomics Workgroup developed a pipeline to impute and merge genomic data across the different SNP arrays to maximize sample size and power to detect associations with a variety of clinical endpoints. The 1000 Genomes cosmopolitan reference panel was used for imputation. Imputation results were evaluated using the following metrics: accuracy of imputation, allelic R (2) (estimated correlation between the imputed and true genotypes), and the relationship between allelic R (2) and minor allele frequency. Computation time and memory resources required by two different software packages (BEAGLE and IMPUTE2) were also evaluated. A number of challenges were encountered due to the complexity of using two different imputation software packages, multiple ancestral populations, and many different genotyping platforms. We present lessons learned and describe the pipeline implemented here to impute and merge genomic data sets. The eMERGE imputed dataset will serve as a valuable resource for discovery, leveraging the clinical data that can be mined from the EHR.