An integrative framework to prioritize genes in more than 500 loci associated with body mass index.

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss.

A large meta-analysis identifies genes associated with anterior uveitis.

Anterior Uveitis (AU) is the inflammation of the anterior part of the eye, the iris and ciliary body and is strongly associated with HLA-B*27. We report AU exome sequencing results from eight independent cohorts consisting of 3,850 cases and 916,549 controls. We identify common genome-wide significant loci in HLA-B (OR = 3.37, p = 1.03e-196) and ERAP1 (OR = 0.86, p = 1.1e-08), and find IPMK (OR = 9.4, p = 4.42e-09) and IDO2 (OR = 3.61, p = 6.16e-08) as genome-wide significant genes based on the burden of rare coding variants. Dividing the cohort into HLA-B*27 positive and negative individuals, we find ERAP1 haplotype is strongly protective only for B*27-positive AU (OR = 0.73, p = 5.2e-10). Investigation of B*27-negative AU identifies a common signal near HLA-DPB1 (rs3117230, OR = 1.26, p = 2.7e-08), risk genes IPMK and IDO2, and several additional candidate risk genes, including ADGFR5, STXBP2, and ACHE. Taken together, we decipher the genetics underlying B*27-positive and -negative AU and identify rare and common genetic signals for both subtypes of disease.

Life is pain: Fibromyalgia as a nexus of multiple liability distributions.

Fibromyalgia is a complex disease of unclear etiology that is complicated by difficulties in diagnosis, treatment, and clinical heterogeneity. To clarify this etiology, healthcare-based data are leveraged to assess the influences on fibromyalgia in several domains. Prevalence is less than 1% of females in our population register data, and about 1/10th that in males. Fibromyalgia often presents with co-occurring conditions including back pain, rheumatoid arthritis, and anxiety. More comorbidities are identified with hospital-associated biobank data, falling into three broad categories of pain-related, autoimmune, and psychiatric disorders. Selecting representative phenotypes with published genome-wide association results for polygenic scoring, we confirm genetic predispositions to psychiatric, pain sensitivity, and autoimmune conditions show associations with fibromyalgia, although these may differ by ancestry group. We conduct a genome-wide association analysis of fibromyalgia in biobank samples, which did not result in any genome-wide significant loci; further studies with increased sample size are necessary to identify specific genetic effects on fibromyalgia. Overall, fibromyalgia appears to have strong clinical and likely genetic links to several disease categories, and could usefully be understood as a composite manifestation of these etiological sources.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways.

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Insights From a Large-Scale Whole-Genome Sequencing Study of Systolic Blood Pressure, Diastolic Blood Pressure, and Hypertension.

BACKGROUND: The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure. METHODS: We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants. RESULTS: Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings (P<5×10-8). Among them, a rare intergenic variant at novel locus, LOC100506274, was associated with lower systolic blood pressure in stage-1 (beta [SE]=-32.6 [6.0]; P=4.99×10-8) but not stage-2 analysis (P=0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=-0.36 [0.07]; P=4.18×10-7) and attained genome-wide significance in stage-2 (beta [SE]=-0.29 [0.03]; P=7.28×10-23). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings (P<1×10-6 and P<1×10-4, respectively). DISCUSSION: We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

Population-scale analysis of common and rare genetic variation associated with hearing loss in adults.

To better understand the genetics of hearing loss, we performed a genome-wide association meta-analysis with 125,749 cases and 469,497 controls across five cohorts. We identified 53/c loci affecting hearing loss risk, including common coding variants in COL9A3 and TMPRSS3. Through exome sequencing of 108,415 cases and 329,581 controls, we observed rare coding associations with 11 Mendelian hearing loss genes, including additive effects in known hearing loss genes GJB2 (Gly12fs; odds ratio [OR] = 1.21, P = 4.2 × 10-11) and SLC26A5 (gene burden; OR = 1.96, P = 2.8 × 10-17). We also identified hearing loss associations with rare coding variants in FSCN2 (OR = 1.14, P = 1.9 × 10-15) and KLHDC7B (OR = 2.14, P = 5.2 × 10-30). Our results suggest a shared etiology between Mendelian and common hearing loss in adults. This work illustrates the potential of large-scale exome sequencing to elucidate the genetic architecture of common disorders where both common and rare variation contribute to risk.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed.

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Whole-Genome Sequencing Association Analyses of Stroke and Its Subtypes in Ancestrally Diverse Populations From Trans-Omics for Precision Medicine Project.

BACKGROUND AND PURPOSE: Stroke is the leading cause of death and long-term disability worldwide. Previous genome-wide association studies identified 51 loci associated with stroke (mostly ischemic) and its subtypes among predominantly European populations. Using whole-genome sequencing in ancestrally diverse populations from the Trans-Omics for Precision Medicine (TOPMed) Program, we aimed to identify novel variants, especially low-frequency or ancestry-specific variants, associated with all stroke, ischemic stroke and its subtypes (large artery, cardioembolic, and small vessel), and hemorrhagic stroke and its subtypes (intracerebral and subarachnoid). METHODS: Whole-genome sequencing data were available for 6833 stroke cases and 27 116 controls, including 22 315 European, 7877 Black, 2616 Hispanic/Latino, 850 Asian, 54 Native American, and 237 other ancestry participants. In TOPMed, we performed single variant association analysis examining 40 million common variants and aggregated association analysis focusing on rare variants. We also combined TOPMed European populations with over 28 000 additional European participants from the UK BioBank genome-wide array data through meta-analysis. RESULTS: In the single variant association analysis in TOPMed, we identified one novel locus 13q33 for large artery at whole-genome-wide significance (P<5.00×10-9) and 4 novel loci at genome-wide significance (P<5.00×10-8), all of which need confirmation in independent studies. Lead variants in all 5 loci are low-frequency but are more common in non-European populations. An aggregation of synonymous rare variants within the gene C6orf26 demonstrated suggestive evidence of association for hemorrhagic stroke (P<3.11×10-6). By meta-analyzing European ancestry samples in TOPMed and UK BioBank, we replicated several previously reported stroke loci including PITX2, HDAC9, ZFHX3, and LRCH1. CONCLUSIONS: We represent the first association analysis for stroke and its subtypes using whole-genome sequencing data from ancestrally diverse populations. While our findings suggest the potential benefits of combining whole-genome sequencing data with populations of diverse genetic backgrounds to identify possible low-frequency or ancestry-specific variants, they also highlight the need to increase genome coverage and sample sizes.

ERAP1, ERAP2, and Two Copies of HLA-Aw19 Alleles Increase the Risk for Birdshot Chorioretinopathy in HLA-A29 Carriers.

Purpose: Birdshot chorioretinopathy (BSCR) is strongly associated with HLA-A29. This study was designed to elucidate the genetic modifiers of BSCR in HLA-A29 carriers. Methods: We sequenced the largest BSCR cohort to date, including 286 cases and 108 HLA-A29-positive controls to determine genome-wide common and rare variant associations. We further typed the HLA alleles of cases and 45,386 HLA-A29 controls of European ancestry to identify HLA alleles that associate with BSCR risk. Results: Carrying a second allele that belongs to the HLA-Aw19 broad antigen family (including HLA-A29, -A30, -A31, and -A33) increases the risk for BSCR (odds ratio [OR] = 4.44; P = 2.2e-03). This result was validated by comparing allele frequencies to large HLA-A29-controlled cohorts (n = 45,386; OR > 2.5; P < 1.3e-06). We also confirm that ERAP1 and ERAP2 haplotypes modulate disease risk. A meta-analysis with an independent dataset confirmed that ERAP1 and ERAP2 haplotypes modulate the risk for disease at a genome-wide significant level: ERAP1-rs27432 (OR = 2.46; 95% confidence interval [CI], 1.85-3.26; P = 4.07e-10), an expression quantitative trait locus (eQTL) decreasing ERAP1 expression; and ERAP2-rs10044354 (OR = 1.95; 95% CI, 1.55-2.44; P = 6.2e-09), an eQTL increasing ERAP2 expression. Furthermore, ERAP2-rs2248374 that disrupts ERAP2 expression is protective (OR = 0.56; 95% CI, 0.45-0.70; P = 2.39e-07). BSCR risk is additively increased when combining ERAP1/ERAP2 risk genotypes with two copies of HLA-Aw19 alleles (OR = 13.53; 95% CI, 3.79-54.77; P = 1.17e-05). Conclusions: The genetic factors increasing BSCR risk demonstrate a pattern of increased processing, as well as increased presentation of ERAP2-specific peptides. This suggests a mechanism in which exceeding a peptide presentation threshold activates the immune response in choroids of A29 carriers.

Use of the PsycheMERGE Network to Investigate the Association Between Depression Polygenic Scores and White Blood Cell Count.

Importance: Although depression is a common psychiatric disorder, its underlying biological basis remains poorly understood. Pairing depression polygenic scores with the results of clinical laboratory tests can reveal biological processes involved in depression etiology and in the physiological changes resulting from depression. Objective: To characterize the association between depression polygenic scores and an inflammatory biomarker, ie, white blood cell count. Design, Setting, and Participants: This genetic association study was conducted from May 19, 2019, to June 5, 2021, using electronic health record data from 382 452 patients across 4 health care systems. Analyses were conducted separately in each health care system and meta-analyzed across all systems. Primary analyses were conducted in Vanderbilt University Medical Center's biobank. Replication analyses were conducted across 3 other PsycheMERGE sites: Icahn School of Medicine at Mount Sinai, Mass General Brigham, and the Million Veteran Program. All patients with available genetic data and recorded white blood cell count measurements were included in the analyses. Primary analyses were conducted in individuals of European descent and then repeated in a population of individuals of African descent. Exposures: Depression polygenic scores. Main Outcomes and Measures: White blood cell count. Results: Across the 4 PsycheMERGE sites, there were 382 452 total participants of European ancestry (18.7% female; median age, 57.9 years) and 12 383 participants of African ancestry (61.1% female; median age, 39.0 [range, birth-90.0 years]). A laboratory-wide association scan revealed a robust association between depression polygenic scores and white blood cell count (ß, 0.03; SE, 0.004; P = 1.07 × 10-17), which was replicated in a meta-analysis across the 4 health care systems (ß, 0.03; SE, 0.002; P = 1.03 × 10-136). Mediation analyses suggested a bidirectional association, with white blood cell count accounting for 2.5% of the association of depression polygenic score with depression diagnosis (95% CI, 2.2%-20.8%; P = 2.84 × 10-70) and depression diagnosis accounting for 9.8% of the association of depression polygenic score with white blood cell count (95% CI, 8.4%-11.1%; P = 1.78 × 10-44). Mendelian randomization provided additional support for an association between increased white blood count and depression risk, but depression modeled as the exposure showed no evidence of an influence on white blood cell counts. Conclusions and Relevance: This genetic association study found that increased depression polygenic scores were associated with increased white blood cell count, and suggests that this association may be bidirectional. These findings highlight the potential importance of the immune system in the etiology of depression and may motivate future development of clinical biomarkers and targeted treatment options for depression.

The Genetic Architecture of Depression in Individuals of East Asian Ancestry: A Genome-Wide Association Study.

Importance: Most previous genome-wide association studies (GWAS) of depression have used data from individuals of European descent. This limits the understanding of the underlying biology of depression and raises questions about the transferability of findings between populations. Objective: To investigate the genetics of depression among individuals of East Asian and European descent living in different geographic locations, and with different outcome definitions for depression. Design, Setting, and Participants: Genome-wide association analyses followed by meta-analysis, which included data from 9 cohort and case-control data sets comprising individuals with depression and control individuals of East Asian descent. This study was conducted between January 2019 and May 2021. Exposures: Associations of genetic variants with depression risk were assessed using generalized linear mixed models and logistic regression. The results were combined across studies using fixed-effects meta-analyses. These were subsequently also meta-analyzed with the largest published GWAS for depression among individuals of European descent. Additional meta-analyses were carried out separately by outcome definition (clinical depression vs symptom-based depression) and region (East Asian countries vs Western countries) for East Asian ancestry cohorts. Main Outcomes and Measures: Depression status was defined based on health records and self-report questionnaires. Results: There were a total of 194 548 study participants (approximate mean age, 51.3 years; 62.8% women). Participants included 15 771 individuals with depression and 178 777 control individuals of East Asian descent. Five novel associations were identified, including 1 in the meta-analysis for broad depression among those of East Asian descent: rs4656484 (ß = -0.018, SE = 0.003, P = 4.43x10-8) at 1q24.1. Another locus at 7p21.2 was associated in a meta-analysis restricted to geographically East Asian studies (ß = 0.028, SE = 0.005, P = 6.48x10-9 for rs10240457). The lead variants of these 2 novel loci were not associated with depression risk in European ancestry cohorts (ß = -0.003, SE = 0.005, P = .53 for rs4656484 and ß = -0.005, SE = 0.004, P = .28 for rs10240457). Only 11% of depression loci previously identified in individuals of European descent reached nominal significance levels in the individuals of East Asian descent. The transancestry genetic correlation between cohorts of East Asian and European descent for clinical depression was r = 0.413 (SE = 0.159). Clinical depression risk was negatively genetically correlated with body mass index in individuals of East Asian descent (r = -0.212, SE = 0.084), contrary to findings for individuals of European descent. Conclusions and Relevance: These results support caution against generalizing findings about depression risk factors across populations and highlight the need to increase the ancestral and geographic diversity of samples with consistent phenotyping.

Myopathic Cardiac Genotypes Increase Risk for Myocarditis.

Impairments in certain cardiac genes confer risk for myocarditis in children. To determine the extent of this association, we performed genomic sequencing in predominantly adult patients with acute myocarditis and matched control subjects. Putatively deleterious variants in a broad set of cardiac genes were found in 19 of 117 acute myocarditis cases vs 34 of 468 control subjects (P = 0.003). Thirteen genes classically associated with cardiomyopathy or neuromuscular disorders with cardiac involvement were implicated, including >1 associated damaging variant in DYSF, DSP, and TTN. Phenotypes of subjects who have acute myocarditis with or without deleterious variants were similar, indicating that genetic testing is necessary to differentiate them.

Burden of Cardiomyopathic Genetic Variation in Lethal Pediatric Myocarditis.

BACKGROUND: Acute myocarditis (AM) is a well-known cause of sudden death and heart failure, often caused by prevalent viruses. We previously showed that some pediatric AM correlates with putatively damaging variants in genes related to cardiomyocyte structure and function. We sought to evaluate whether deleterious cardiomyopathic variants were enriched among fatal pediatric AM cases in New York City compared with ancestry-matched controls. METHODS: Twenty-four children (aged 3 weeks to 20 years) with death due to AM were identified through autopsy records; histologies were reviewed to confirm that all cases met Dallas criteria for AM and targeted panel sequencing of 57 cardiomyopathic genes was performed. Controls without cardiovascular disease were identified from a pediatric database and matched by genetic ancestry to cases using principal components from exome sequencing. Rates of putative deleterious variations (DV) were compared between cases and controls. Where available, AM tissues underwent viral analysis by polymerase chain reaction. RESULTS: DV were identified in 4 of 24 AM cases (16.7%), compared with 2 of 96 age and ancestry-matched controls (2.1%, P=0.014). Viral causes were proven for 6 of 8 AM cases (75%), including the one DV+ case where tissue was available for testing. DV+ cases were more likely to be female, have no evidence of chronic inflammation, and associate with sudden cardiac death than DV- cases. CONCLUSIONS: Deleterious variants in genes related to cardiomyocyte integrity are more common in children with fatal AM than controls, likely conferring susceptibility. Additionally, genetically mediated AM may progress more rapidly and be more severe.

Rare Coding Variants Associated With Electrocardiographic Intervals Identify Monogenic Arrhythmia Susceptibility Genes: A Multi-Ancestry Analysis.

BACKGROUND: Alterations in electrocardiographic (ECG) intervals are well-known markers for arrhythmia and sudden cardiac death (SCD) risk. While the genetics of arrhythmia syndromes have been studied, relations between electrocardiographic intervals and rare genetic variation at a population level are poorly understood. METHODS: Using a discovery sample of 29 000 individuals with whole-genome sequencing from Trans-Omics in Precision Medicine and replication in nearly 100 000 with whole-exome sequencing from the UK Biobank and MyCode, we examined associations between low-frequency and rare coding variants with 5 routinely measured electrocardiographic traits (RR, P-wave, PR, and QRS intervals and corrected QT interval). RESULTS: We found that rare variants associated with population-based electrocardiographic intervals identify established monogenic SCD genes (KCNQ1, KCNH2, and SCN5A), a controversial monogenic SCD gene (KCNE1), and novel genes (PAM and MFGE8) involved in cardiac conduction. Loss-of-function and pathogenic SCN5A variants, carried by 0.1% of individuals, were associated with a nearly 6-fold increased odds of the first-degree atrioventricular block (P=8.4×10-5). Similar variants in KCNQ1 and KCNH2 (0.2% of individuals) were associated with a 23-fold increased odds of marked corrected QT interval prolongation (P=4×10-25), a marker of SCD risk. Incomplete penetrance of such deleterious variation was common as over 70% of carriers had normal electrocardiographic intervals. CONCLUSIONS: Our findings indicate that large-scale high-depth sequence data and electrocardiographic analysis identifies monogenic arrhythmia susceptibility genes and rare variants with large effects. Known pathogenic variation in conventional arrhythmia and SCD genes exhibited incomplete penetrance and accounted for only a small fraction of marked electrocardiographic interval prolongation.

Toward a fine-scale population health monitoring system.

Understanding population health disparities is an essential component of equitable precision health efforts. Epidemiology research often relies on definitions of race and ethnicity, but these population labels may not adequately capture disease burdens and environmental factors impacting specific sub-populations. Here, we propose a framework for repurposing data from electronic health records (EHRs) in concert with genomic data to explore the demographic ties that can impact disease burdens. Using data from a diverse biobank in New York City, we identified 17 communities sharing recent genetic ancestry. We observed 1,177 health outcomes that were statistically associated with a specific group and demonstrated significant differences in the segregation of genetic variants contributing to Mendelian diseases. We also demonstrated that fine-scale population structure can impact the prediction of complex disease risk within groups. This work reinforces the utility of linking genomic data to EHRs and provides a framework toward fine-scale monitoring of population health.

Implementing genomic screening in diverse populations.

BACKGROUND: Population-based genomic screening has the predicted ability to reduce morbidity and mortality associated with medically actionable conditions. However, much research is needed to develop standards for genomic screening and to understand the perspectives of people offered this new testing modality. This is particularly true for non-European ancestry populations who are vastly underrepresented in genomic medicine research. Therefore, we implemented a pilot genomic screening program in the BioMe Biobank in New York City, where the majority of participants are of non-European ancestry. METHODS: We initiated genomic screening for well-established genes associated with hereditary breast and ovarian cancer syndrome (HBOC), Lynch syndrome (LS), and familial hypercholesterolemia (FH). We evaluated and included an additional gene (TTR) associated with hereditary transthyretin amyloidosis (hATTR), which has a common founder variant in African ancestry populations. We evaluated the characteristics of 74 participants who received results associated with these conditions. We also assessed the preferences of 7461 newly enrolled BioMe participants to receive genomic results. RESULTS: In the pilot genomic screening program, 74 consented participants received results related to HBOC (N = 26), LS (N = 6), FH (N = 8), and hATTR (N = 34). Thirty-three of 34 (97.1%) participants who received a result related to hATTR were self-reported African American/African (AA) or Hispanic/Latinx (HL), compared to 14 of 40 (35.0%) participants who received a result related to HBOC, LS, or FH. Among the 7461 participants enrolled after the BioMe protocol modification to allow the return of genomic results, 93.4% indicated that they would want to receive results. Younger participants, women, and HL participants were more likely to opt to receive results. CONCLUSIONS: The addition of TTR to a pilot genomic screening program meant that we returned results to a higher proportion of AA and HL participants, in comparison with genes traditionally included in genomic screening programs in the USA. We found that the majority of participants in a multi-ethnic biobank are interested in receiving genomic results for medically actionable conditions. These findings increase knowledge about the perspectives of diverse research participants on receiving genomic results and inform the broader implementation of genomic medicine in underrepresented patient populations.

Inflamed Ulcerative Colitis Regions Associated With MRGPRX2-Mediated Mast Cell Degranulation and Cell Activation Modules, Defining a New Therapeutic Target.

BACKGROUND & AIMS: Recent literature has implicated a key role for mast cells in murine models of colonic inflammation, but their role in human ulcerative colitis (UC) is not well established. A major advance has been the identification of mrgprb2 (human orthologue, MRGPX2) as mediating IgE-independent mast cell activation. We sought to define mechanisms of mast cell activation and MRGPRX2 in human UC. METHODS: Colon tissues were collected from patients with UC for bulk RNA sequencing and lamina propria cells were isolated for MRGPRX2 activation studies and single-cell RNA sequencing. Genetic association of all protein-altering G-protein coupled receptor single-nucleotide polymorphism was performed in an Ashkenazi Jewish UC case-control cohort. Variants of MRGPRX2 were transfected into Chinese hamster ovary (CHO) and human mast cell (HMC) 1.1 cells to detect genotype-dependent effects on ß-arrestin recruitment, IP-1 accumulation, and phosphorylated extracellular signal-regulated kinase. RESULTS: Mast cell-specific mediators and adrenomedullin (proteolytic precursor of PAMP-12, an MRGPRX2 agonist) are up-regulated in inflamed compared to uninflamed UC. MRGPRX2 stimulation induces carboxypeptidase secretion from inflamed UC. Of all protein-altering GPCR alleles, a unique variant of MRGPRX2, Asn62Ser, was most associated with and was bioinformatically predicted to alter arrestin recruitment. We validated that the UC protective serine allele enhances ß-arrestin recruitment, decreases IP-1, and increases phosphorylated extracellular signal-regulated kinase with MRGPRX2 agonists. Single-cell RNA sequencing defines that adrenomedullin is expressed by activated fibroblasts and epithelial cells and that interferon gamma is a key upstream regulator of mast cell gene expression. CONCLUSION: Inflamed UC regions are distinguished by MRGPRX2-mediated activation of mast cells, with decreased activation observed with a UC-protective genetic variant. These results define cell modules of UC activation and a new therapeutic target.

A genotype-first approach to exploring Mendelian cardiovascular traits with clear external manifestations.

PURPOSE: The purpose of this study is to use a genotype-first approach to explore highly penetrant, autosomal dominant cardiovascular diseases with external features, the RASopathies and Marfan syndrome (MFS), using biobank data. METHODS: This study uses exome sequencing and corresponding phenotypic data from Mount Sinai's BioMe (n = 32,344) and the United Kingdom Biobank (UKBB; n = 49,960). Variant curation identified pathogenic/likely pathogenic (P/LP) variants in RASopathy genes and FBN1. RESULTS: Twenty-one subjects harbored P/LP RASopathy variants; three (14%) were diagnosed, and another 46% had ≥1 classic Noonan syndrome (NS) feature. Major NS features (short stature [9.5% p = 7e-5] and heart anomalies [19%, p < 1e-5]) were less frequent than expected. Prevalence of hypothyroidism/autoimmune disorders was enriched compared with biobank populations (p = 0.007). For subjects with FBN1 P/LP variants, 14/41 (34%) had a MFS diagnosis or highly suggestive features. Five of 15 participants (33%) with echocardiographic data had aortic dilation, fewer than expected (p = 8e-6). Ectopia lentis affected only 15% (p < 1e-5). CONCLUSIONS: Substantial fractions of individuals harboring P/LP variants with partial or full phenotypic matches to a RASopathy or MFS remain undiagnosed, some not meeting diagnostic criteria. Routine population genotyping would enable multidisciplinary care and avoid life-threatening events.

Common and Rare Variant Prediction and Penetrance of IBD in a Large, Multi-ethnic, Health System-based Biobank Cohort.

BACKGROUND AND AIMS: Polygenic risk scores (PRS) may soon be used to predict inflammatory bowel disease (IBD) risk in prevention efforts. We leveraged exome-sequence and single nucleotide polymorphism (SNP) array data from 29,358 individuals in the multiethnic, randomly ascertained health system-based BioMe biobank to define effects of common and rare IBD variants on disease prediction and pathophysiology. METHODS: PRS were calculated from European, African American, and Ashkenazi Jewish (AJ) reference case-control studies, and a meta-GWAS run using all three association datasets. PRS were then combined using regression to assess which combination of scores best predicted IBD status in European, AJ, Hispanic, and African American cohorts in BioMe. Additionally, rare variants were assessed in genes associated with very early-onset IBD (VEO-IBD), by estimating genetic penetrance in each BioMe population. RESULTS: Combining risk scores based on association data from distinct ancestral populations improved IBD prediction for every population in BioMe and significantly improved prediction among European ancestry UK Biobank individuals. Lower predictive power for non-Europeans was observed, reflecting in part substantially lower African IBD case-control reference sizes. We replicated associations for two VEO-IBD genes, ADAM17 and LRBA, with high dominant model penetrance in BioMe. Autosomal recessive LRBA risk alleles are associated with severe, early-onset autoimmunity; we show that heterozygous carriage of an African-predominant LRBA protein-altering allele is associated with significantly decreased LRBA and CTLA-4 expression with T-cell activation. CONCLUSIONS: Greater genetic diversity in African populations improves prediction across populations, and generalizes some VEO-IBD genes. Increasing African American IBD case-collections should be prioritized to reduce health disparities and enhance pathophysiological insight.