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

Abstract

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.

Keywords: Asthma; Copy number variants; Exome sequencing; Fine-mapping; Genetic association; Human leukocyte antigen; UK Biobank.

Fig. 1

Number of UK Biobank participants meeting study selection criteria. *Asthma cases defined as having self-reported doctor-diagnosed asthma or an ICD10 code for asthma in hospital inpatient records. **Controls defined as having no self-reported doctor-diagnosed asthma, no ICD10 code for asthma, and no self-reported asthma from nurse’s interview

Fig. 2

Association of HLA variants with asthma in UK Biobank. Locuszoom plots of the HLA region showing association of HLA region variants with asthma (top panel) and association of HLA region variants with asthma conditional on the SNP with the top posterior inclusion probability in credible set 1 (bottom panel). Variants in credible set 1 are highlighted in red and variants in credible set 2 are highlighted in blue