Multiethnic genome-wide and HLA association study of total serum IgE level

Genetic and Immunological Pathogenesis of Atopic Dermatitis

Type 2 immune-mediated diseases give a clear answer to the issue of nature (genetics) versus nurture (environment). Both genetics and environment play vital complementary roles in the development of atopic dermatitis (AD). As a key component of the atopic march, AD demonstrates the interactive nature of genetic and environmental contributions to atopy. From sequence variants in the epithelial barrier gene encoding FLG to the hygiene hypothesis, AD combines a broad array of contributions into a single syndrome. This review will focus on the genetic contribution to AD and where genetics facilitates the elicitation or enhancement of AD pathogenesis.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

Meta-regression of genome-wide association studies to estimate age-varying genetic effects

Fixed-effect meta-analysis has been used to summarize genetic effects on a phenotype across multiple Genome-Wide Association Studies (GWAS) assuming a common underlying genetic effect. Genetic effects may vary with age (or other characteristics), and not allowing for this in a GWAS might lead to bias. Meta-regression models between study heterogeneity and allows effect modification of the genetic effects to be explored. The aim of this study was to explore the use of meta-analysis and meta-regression for estimating age-varying genetic effects on phenotypes. With simulations we compared the performance of meta-regression to fixed-effect and random -effects meta-analyses in estimating (i) main genetic effects and (ii) age-varying genetic effects (SNP by age interactions) from multiple GWAS studies under a range of scenarios. We applied meta-regression on publicly available summary data to estimate the main and age-varying genetic effects of the FTO SNP rs9939609 on Body Mass Index (BMI). Fixed-effect and random-effects meta-analyses accurately estimated genetic effects when these did not change with age. Meta-regression accurately estimated both main genetic effects and age-varying genetic effects. When the number of studies or the age-diversity between studies was low, meta-regression had limited power. In the applied example, each additional minor allele (A) of rs9939609 was inversely associated with BMI at ages 0 to 3, and positively associated at ages 5.5 to 13. Our findings challenge the assumption that genetic effects are consistent across all ages and provide a method for exploring this. GWAS consortia should be encouraged to use meta-regression to explore age-varying genetic effects.

Epigenome-wide DNA methylation association study of circulating IgE levels identifies novel targets for asthma

BACKGROUND

Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases.

METHODS

We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables.

FINDINGS

We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P < 7.94E-04) that can be explored as potential therapeutic targets.

INTERPRETATION

By integrating eQTM and MR analyses in general and clinical asthma populations, our findings provide a deeper understanding of the multidimensional inter-relations of DNA methylation, gene expression, and IgE levels.

FUNDING

US NIH/NHLBI grants: P01HL132825, K99HL159234. N01-HC-25195 and HHSN268201500001I.

A Review of the Epigenetic Clock: Emerging Biomarkers for Asthma and Allergic Disease

DNA methylation (DNAm) is a dynamic, age-dependent epigenetic modification that can be used to study interactions between genetic and environmental factors. Environmental exposures during critical periods of growth and development may alter DNAm patterns, leading to increased susceptibility to diseases such as asthma and allergies. One method to study the role of DNAm is the epigenetic clock-an algorithm that uses DNAm levels at select age-informative Cytosine-phosphate-Guanine (CpG) dinucleotides to predict epigenetic age (EA). The difference between EA and calendar age (CA) is termed epigenetic age acceleration (EAA) and reveals information about the biological capacity of an individual. Associations between EAA and disease susceptibility have been demonstrated for a variety of age-related conditions and, more recently, phenotypes such as asthma and allergic diseases, which often begin in childhood and progress throughout the lifespan. In this review, we explore different epigenetic clocks and how they have been applied, particularly as related to childhood asthma. We delve into how in utero and early life exposures (e.g., smoking, air pollution, maternal BMI) result in methylation changes. Furthermore, we explore the potential for EAA to be used as a biomarker for asthma and allergic diseases and identify areas for further study.

Sequence variant affects GCSAML splicing, mast cell specific proteins, and risk of urticaria

Urticaria is a skin disorder characterized by outbreaks of raised pruritic wheals. In order to identify sequence variants associated with urticaria, we performed a meta-analysis of genome-wide association studies for urticaria with a total of 40,694 cases and 1,230,001 controls from Iceland, the UK, Finland, and Japan. We also performed transcriptome- and proteome-wide analyses in Iceland and the UK. We found nine sequence variants at nine loci associating with urticaria. The variants are at genes participating in type 2 immune responses and/or mast cell biology (CBLB, FCER1A, GCSAML, STAT6, TPSD1, ZFPM1), the innate immunity (C4), and NF-κB signaling. The most significant association was observed for the splice-donor variant rs56043070[A] (hg38: chr1:247556467) in GCSAML (MAF = 6.6%, OR = 1.24 (95%CI: 1.20-1.28), P-value = 3.6 × 10^-44). We assessed the effects of the variants on transcripts, and levels of proteins relevant to urticaria pathophysiology. Our results emphasize the role of type 2 immune response and mast cell activation in the pathogenesis of urticaria. Our findings may point to an IgE-independent urticaria pathway that could help address unmet clinical need.

Race/ethnicity-stratified fine-mapping of the MHC locus reveals genetic variants associated with late-onset asthma

Introduction: Asthma is a chronic disease of the airways that impairs normal breathing. The etiology of asthma is complex and involves multiple factors, including the environment and genetics, especially the distinct genetic architecture associated with ancestry. Compared to early-onset asthma, little is known about genetic predisposition to late-onset asthma. We investigated the race/ethnicity-specific relationship among genetic variants within the major histocompatibility complex (MHC) region and late-onset asthma in a North Carolina-based multiracial cohort of adults. Methods: We stratified all analyses by self-reported race (i.e., White and Black) and adjusted all regression models for age, sex, and ancestry. We conducted association tests within the MHC region and performed fine-mapping analyses conditioned on the race/ethnicity-specific lead variant using whole-genome sequencing (WGS) data. We applied computational methods to infer human leukocyte antigen (HLA) alleles and residues at amino acid positions. We replicated findings in the UK Biobank. Results: The lead signals, rs9265901 on the 5' end of HLA-B, rs55888430 on HLA-DOB, and rs117953947 on HCG17, were significantly associated with late-onset asthma in all, White, and Black participants, respectively (OR = 1.73, 95%CI: 1.31 to 2.14, p = 3.62 × 10^-5; OR = 3.05, 95%CI: 1.86 to 4.98, p = 8.85 × 10^-6; OR = 19.5, 95%CI: 4.37 to 87.2, p = 9.97 × 10^-5, respectively). For the HLA analysis, HLA-B*40:02 and HLA-DRB1*04:05, HLA-B*40:02, HLA-C*04:01, and HLA-DRB1*04:05, and HLA-DRB1*03:01 and HLA-DQB1 were significantly associated with late-onset asthma in all, White, and Black participants. Conclusion: Multiple genetic variants within the MHC region were significantly associated with late-onset asthma, and the associations were significantly different by race/ethnicity group.

Human germline heterozygous gain-of-function STAT6 variants cause severe allergic disease

STAT6 (signal transducer and activator of transcription 6) is a transcription factor that plays a central role in the pathophysiology of allergic inflammation. We have identified 16 patients from 10 families spanning three continents with a profound phenotype of early-life onset allergic immune dysregulation, widespread treatment-resistant atopic dermatitis, hypereosinophilia with esosinophilic gastrointestinal disease, asthma, elevated serum IgE, IgE-mediated food allergies, and anaphylaxis. The cases were either sporadic (seven kindreds) or followed an autosomal dominant inheritance pattern (three kindreds). All patients carried monoallelic rare variants in STAT6 and functional studies established their gain-of-function (GOF) phenotype with sustained STAT6 phosphorylation, increased STAT6 target gene expression, and TH2 skewing. Precision treatment with the anti-IL-4Rα antibody, dupilumab, was highly effective improving both clinical manifestations and immunological biomarkers. This study identifies heterozygous GOF variants in STAT6 as a novel autosomal dominant allergic disorder. We anticipate that our discovery of multiple kindreds with germline STAT6 GOF variants will facilitate the recognition of more affected individuals and the full definition of this new primary atopic disorder.

Novel insights into atopic dermatitis

Recent research into the pathophysiology and treatment of atopic dermatitis (AD) has shown notable progress. An increasing number of aspects of the immune system are being implicated in AD, including the epithelial barrier, TH2 cytokines, and mast cells. Major advances in therapeutics were made in biologic cytokine and receptor antagonists and among Janus kinase inhibitors. We focus on these areas and address new insights into AD epidemiology, biomarkers, endotypes, prevention, and comorbidities. Going forward, we expect future mechanistic insights and therapeutic advances to broaden physicians' ability to diagnose and manage AD patients, and perhaps to find a cure for this chronic condition.

Targeted capture enrichment and sequencing identifies HLA variants associated with the severity of COVID-19

BACKGROUND

Coronavirus disease 2019 (COVID-19) is currently a global pandemic. The pathogenesis of severe COVID-19 has been widely investigated, but it is still unclear. Human leukocyte antigen (HLA) plays a central role in immune response, and its variants might be related to COVID-19 progression and severity.

OBJECTIVE

To investigate the hypothesis that individual HLA variations could alter the course of COVID-19 and might be associated with the severity of COVID-19.

METHODS

In this study, we conducted an HLA targeted capture enrichment and sequencing of severe COVID-19 patients matched to mild cases. A total of 16 COVID-19 patients, confirmed by SARS-CoV-2 viral RNA polymerase-chain-reaction (PCR) test and chest computed tomography (CT) scan, were enrolled in this study. The HLA targeted capture enrichment and sequencing were conducted. HLA typing was performed by comparing contigs with IPD-IMGT/HLA Database.

RESULTS

In this study, 139 four-digit resolution HLA alleles were acquired. The results showed that HLA-DRB3*01:01 allele was significantly associated with the severity of COVID-19 (odds ratio [OR] = 27.64, 95% confidence interval [CI] = 1.35-560.50, P = 0.0064). And HLA-K*01:01 might be a potential risk factor for COVID-19 severity (OR = 0.11, 95% CI = 0.017-0.66, P = 0.019), but HLA-K*01:02 might be a protective factor (OR = 7.50, 95% CI = 1.48-37.92, P = 0.019).

CONCLUSION

Three non-classical HLA alleles, including HLA-DRB3*01:01, HLA-K*01:01, HLA-K*01:02 were identified to be associated with the severity of COVID-19 by comparing mild and severe patients. The current findings would be helpful for exploring the influence of HLA gene polymorphisms on the development and severity of COVID-19.

Recent progress in the genetic and epigenetic underpinnings of atopy

In the past 2 years, there continue to be advances in our understanding of the genetic and epigenetic underpinnings of atopy pertaining to disease risk and disease severity. The joint role of genetics and the environment has been emphasized in multiple studies. Combining genetics with family history, biomarkers, and comorbidities is further refining our ability to predict the development of individual atopic diseases as well as the advancement of the atopic march. Polygenic risk scores will be an important next step for the field moving toward clinical translation of the genetic findings thus far. A systems biology approach, as illustrated by studies of the microbiome and epigenome, will be necessary to fully understand disease development and to develop increasingly targeted therapeutics.

Genetic Determinants in Airways Obstructive Diseases: The Case of Asthma Chronic Obstructive Pulmonary Disease Overlap

Genome-wide association studies (GWAS) of asthma and chronic obstructive pulmonary disease (COPD) with ever-increasing sample sizes have found multiple genetic loci associated with either disease. However, there are few intersecting loci between asthma and COPD. GWAS specifically focused on asthma-COPD overlap (ACO) have been limited by smaller sample sizes and the lack of a consistent definition of ACO that has also hampered clinical and epidemiologic studies. Other genomic techniques, such as gene expression profiling, are feasible with smaller sample sizes. Genetic analyses of objective measures of airway reactivity and allergy/T2 inflammation biomarkers in COPD studies may be another strategy to overcome limitations in ACO definitions.

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank

BACKGROUND

The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma.

METHODS

We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping.

RESULTS

Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11).

CONCLUSIONS

At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.