Exome-wide evidence of compound heterozygous effects across common phenotypes in the UK Biobank

Exome-sequencing association studies have successfully linked rare protein-coding variation to risk of thousands of diseases. However, the relationship between rare deleterious compound heterozygous (CH) variation and their phenotypic impact has not been fully investigated. Here, we leverage advances in statistical phasing to accurately phase rare variants (MAF ~ 0.001%) in exome sequencing data from 175,587 UK Biobank (UKBB) participants, which we then systematically annotate to identify putatively deleterious CH coding variation. We show that 6.5% of individuals carry such damaging variants in the CH state, with 90% of variants occurring at MAF < 0.34%. Using a logistic mixed model framework, systematically accounting for relatedness, polygenic risk, nearby common variants, and rare variant burden, we investigate recessive effects in common complex diseases. We find six exome-wide significant (P < 1.68 × 10 - 7 ) and 17 nominally significant (P < 5.25 × 10 - 5 ) gene-trait associations. Among these, only four would have been identified without accounting for CH variation in the gene. We further incorporate age-at-diagnosis information from primary care electronic health records, to show that genetic phase influences lifetime risk of disease across 20 gene-trait combinations (FDR < 5%). Using a permutation approach, we find evidence for genetic phase contributing to disease susceptibility for a collection of gene-trait pairs, including FLG-asthma (P = 0.00205 ) and USH2A-visual impairment (P = 0.0084 ). Taken together, we demonstrate the utility of phasing large-scale genetic sequencing cohorts for robust identification of the phenome-wide consequences of compound heterozygosity.

Filaggrin loss-of-function mutations, atopic dermatitis and risk of actinic keratosis: results from two cross-sectional studies

BACKGROUND

Common loss-of-function mutations in filaggrin gene (FLG) represent a strong genetic risk factor for atopic dermatitis (AD). Homozygous mutation carriers typically display ichthyosis vulgaris (IV) and many have concomitant AD. Previously, homozygous, but not heterozygous, filaggrin gene mutations have been associated with squamous cell carcinomas.

OBJECTIVE

The first objective was to examine the association between FLG mutations and actinic keratosis (AK). The second objective was to investigate the occurrence of AK in patients with IV and AD, respectively.

METHODS

FLG mutation status in patients with AK was compared with controls from the general population. Furthermore, based on nationwide data from Danish registers, we compared the risk of AK in patients with IV, AD and psoriasis, respectively.

RESULTS

The prevalence of homozygous FLG mutations was significantly higher in the AK group (n = 4, 0.8%) in comparison with the control group (n = 18, 0.2%), whereas the prevalence of heterozygous FLG mutations was lower. In hospital registry data, patients with AD exhibited an increased risk of AK than did psoriasis controls (adjusted OR 1.46; [95% CI 1.12-1.90]), whereas no difference in risk was observed between patients with IV and AD.

CONCLUSIONS

This study indicates an increased susceptibility to AK in individuals with homozygous, but not heterozygous, FLG mutations and in patients with AD compared to psoriasis. Whether a reduction or absence of epidermal filaggrin could contribute to the susceptibility to AK in patients with IV and AD is unknown and additional research is needed to further explore this relationship.

Predicting the combined effect of multiple genetic variants

Background

Many genetic variants have been identified in the human genome. The functional effects of a single variant have been intensively studied. However, the joint effects of multiple variants in the same genes have been largely ignored due to their complexity or lack of data. This paper uses HMMvar, a hidden Markov model based approach, to investigate the combined effect of multiple variants from the 1000 Genomes Project. Two tumor suppressor genes, TP53 and phosphatase and tensin homolog (PTEN), are also studied for the joint effect of compensatory indel variants.

Results

Results show that there are cases where the joint effect of having multiple variants in the same genes is significantly different from that of a single variant. The deleterious effect of a single indel variant can be alleviated by their compensatory indels in TP53 and PTEN. Compound mutations in two genes, β-MHC and MyBP-C, leading to severer cardiovascular disease compared to single mutations, are also validated.

Conclusions

This paper extends the functionality of HMMvar, a tool for assigning a quantitative score to a variant, to measure not only the deleterious effect of a single variant but also the joint effect of multiple variants. HMMvar is the first tool that can predict the functional effects of both single and general multiple variations on proteins. The precomputed scores for multiple variants from the 1000 Genomes Project and the HMMvar package are available at https://bioinformatics.cs.vt.edu/zhanglab/HMMvar/