Pathogenicity and selective constraint on variation near splice sites

Figure 1.. Signals of purifying selection around splice sites. (A) Selective constraint across splicing region in 13,750 unaffected parents of DDD probands and more than 60,000 aggregated exomes from ExAC. Mutability-adjusted proportion of singletons (MAPS) with 95% confidence intervals (CIs) shown for Ensembl's Variant Effect Predictor (VEP) annotated exonic sites, extended splice acceptor and splice donor regions, the last base of the exon, split by reference nucleotide, and grouped sites in the polypyrimidine tract (PolyPy) region, split by changes from a pyrimidine to a purine (PyPu) versus all other changes. (B) Proportion of variants with 95% CI in 13,750 unaffected parents of DDD probands that fall within genes with high probability of loss-of-function intolerance (pLI > 0.9) across VEP annotated exonic sites, extended splice acceptor and splice donor regions, the last base of the exon, split by reference nucleotide, and grouped sites in the PolyPy region, split by changes from a pyrimidine to a purine (PyPu) versus all other changes. Lower panel shows splice acceptor and splice donor consensus sequences, based on our exons of interest.

Figure 2.. De novo mutations (DNMs) around splice sites. Enrichment of DNMs across the splicing region in 7833 DDD probands. (A) Numbers of observed and expected DNMs across the splicing region in known dominant and recessive DD genes, as well as in non-DD–associated genes, with FDR-corrected Poisson P-values. Splice acceptor and splice donor consensus sequences are shown below, as in Figure 1. (B) Aggregation of observed and expected numbers of DNMs in the PolyPy region, with changes from a pyrimidine to a purine (PyPu) and all other changes shown separately for known dominant and recessive DD genes, as well as non-DD–associated genes. (C) Positive predictive values (PPVs) for DNMs in dominant DD-associated genes in positions across the splicing region, as well as VEP annotated stop gained and missense changes, calculated from observed and expected numbers of DNMs. (D) Enrichment (observed/expected) of DNMs by gene probability of pLI split into sextiles for donor+5, pyrimidine to purine PolyPy, and synonymous sites. pLI scores encompassed by each sextile: 1 = 5.36 × 10−91–0.000000605, 2 = 0.000000609–0.000558185, 3 = 0.000559475–0.027905143, 4 = 0.027908298–0.377456159, 5 =…

Figure 3.. Clinical classifications of noncanonical near-splice DNMs. Relationship between clinical classifications of 38 splice region DNMs in undiagnosed DDD probands and PPVs calculated using observed and expected numbers of DNMs in 7833 probands.

Figure 4.. Selective constraint and pathogenicity scores. MAPS, with 95% CI, calculated for pathogenicity score brackets (least to most severe) in 13,750 unaffected parents from the DDD project, with Spearman's rank correlation coefficient.

Figure 5.. Pathogenicity scores for observed near-splice site DNMs. Cumulative percentage of DNMs in known dominant DD genes with decreasing pathogenicity score bracket, shown with canonical splice site positions included (left) and excluded (right). (AUC*) area under curve.