Person:

Kurki, Mitja

Summary Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. We describe the aggregation and analysis of high-quality exome (protein-coding region) sequence data for 60,706 individuals of diverse ethnicities generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of truncating variants with 72% having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human “knockout” variants in protein-coding genes.

Background: Intracranial aneurysms (IAs), abdominal aortic aneurysms (AAAs), and thoracic aortic aneurysms (TAAs) all have a familial predisposition. Given that aneurysm types are known to co‐occur, we hypothesized that there may be shared genetic risk factors for IAs, AAAs, and TAAs. Methods and Results: We performed a mega‐analysis of 1000 Genomes Project‐imputed genome‐wide association study (GWAS) data of 4 previously published aneurysm cohorts: 2 IA cohorts (in total 1516 cases, 4305 controls), 1 AAA cohort (818 cases, 3004 controls), and 1 TAA cohort (760 cases, 2212 controls), and observed associations of 4 known IA, AAA, and/or TAA risk loci (9p21, 18q11, 15q21, and 2q33) with consistent effect directions in all 4 cohorts. We calculated polygenic scores based on IA‐, AAA‐, and TAA‐associated SNPs and tested these scores for association to case‐control status in the other aneurysm cohorts; this revealed no shared polygenic effects. Similarly, linkage disequilibrium–score regression analyses did not show significant correlations between any pair of aneurysm subtypes. Last, we evaluated the evidence for 14 previously published aneurysm risk single‐nucleotide polymorphisms through collaboration in extended aneurysm cohorts, with a total of 6548 cases and 16 843 controls (IA) and 4391 cases and 37 904 controls (AAA), and found nominally significant associations for IA risk locus 18q11 near RBBP8 to AAA (odds ratio [OR]=1.11; P=4.1×10−5) and for TAA risk locus 15q21 near FBN1 to AAA (OR=1.07; P=1.1×10−3). Conclusions: Although there was no evidence for polygenic overlap between IAs, AAAs, and TAAs, we found nominally significant effects of two established risk loci for IAs and TAAs in AAAs. These two loci will require further replication.

Disruptive and damaging ultra-rare variants (URVs) in highly constrained (HC) genes are enriched in individuals with neurodevelopmental disorders. In the general population, this class of variants was associated with a decrease in years of education (YOE; −3.1 months; P-value=3.3×10−8). This effect was stronger among high brain-expressed genes and explained more YOE variance than pathogenic copy number variation, but less than common variants. Disruptive and damaging URVs in HC genes influence the determinants of YOE in the general population.

Protein-truncating variants protective against human disease provide in vivo validation of therapeutic targets. Here we used targeted sequencing to conduct a search for protein-truncating variants conferring protection against inflammatory bowel disease exploiting knowledge of common variants associated with the same disease. Through replication genotyping and imputation we found that a predicted protein-truncating variant (rs36095412, p.R179X, genotyped in 11,148 ulcerative colitis patients and 295,446 controls, MAF=up to 0.78%) in RNF186, a single-exon ring finger E3 ligase with strong colonic expression, protects against ulcerative colitis (overall P=6.89 × 10−7, odds ratio=0.30). We further demonstrate that the truncated protein exhibits reduced expression and altered subcellular localization, suggesting the protective mechanism may reside in the loss of an interaction or function via mislocalization and/or loss of an essential transmembrane domain.

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10–100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10−16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Identifying causal factors for Mendelian and common diseases is an ongoing challenge in medical genetics 1 . Population bottleneck events, such as in the Finnish population history, enrich some homozygous variants to higher frequencies which facilitates the identification of variants that cause disease with recessive inheritance 2,3 . Here, we study homozygous and heterozygous effects of 44,370 coding variants on 2,444 disease phenotypes using nationwide electronic health record (EHR) data of 176,899 Finns. We find associations for homozygous genotypes across a broad spectrum of phenotypes including known associations to retinal dystrophy, and novel associations to adult-onset cataract and female infertility. 13/20 of recessive disease associations would have been missed by the additive model traditionally used in genome-wide association studies (GWAS). We further use these results to find many known Mendelian variants whose inheritance cannot be adequately described by a traditional definition of dominant or recessive. In particular, we find significant heterozygous effects of variants known to cause disease with recessive inheritance, as well as for reported benign variants. Our results demonstrate how biobanks, particularly in founder populations, can broaden our understanding of the sometimes more complex dosage effects of Mendelian variants on disease.