Genetic constraint at single amino acid resolution in protein domains improves missense variant prioritisation and gene discovery

: Background: One of the major hurdles in clinical genetics is interpreting the clinical consequences associated with germline missense variants in humans. Recent significant advances have leveraged natural variation observed in large-scale human populations to uncover genes or genomic regions that show a depletion of natural variation, indicative of selection pressure. We refer to this as “genetic constraint”. Although existing genetic constraint metrics have been demonstrated to be successful in prioritising genes or genomic regions associated with diseases, their spatial resolution is limited in distinguishing pathogenic variants from benign variants within genes. Methods: We aim to identify missense variants that are significantly depleted in the general human population. Given the size of currently available human populations with exome or genome sequencing data, it is not possible to directly detect depletion of individual missense variants, since the average expected number of observations of a variant at most positions is less than one. We instead focus on protein domains, grouping homologous variants with similar functional impacts to examine the depletion of natural variations within these comparable sets. To accomplish this, we develop the Homologous Missense Constraint (HMC) score. We utilise the Genome Aggregation Database (gnomAD) 125 K exome sequencing data and evaluate genetic constraint at quasi amino-acid resolution by combining signals across protein homologues. Results: We identify one million possible missense variants under strong negative selection within protein domains. Though our approach annotates only protein domains, it nonetheless allows us to assess 22% of the exome confidently. It precisely distinguishes pathogenic variants from benign variants for both early-onset and adult-onset disorders. It outperforms existing constraint metrics and pathogenicity meta-predictors in prioritising de novo mutations from probands with developmental disorders (DD). It is also methodologically independent of these, adding power to predict variant pathogenicity when used in combination. We demonstrate utility for gene discovery by identifying seven genes newly significantly associated with DD that could act through an altered-function mechanism. Conclusions: Grouping variants of comparable functional impacts is effective in evaluating their genetic constraint. HMC is a novel and accurate predictor of missense consequence for improved variant interpretation

Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci.

Multiple regulatory elements distant from their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we develop a capture Hi-C (cHi-C) approach to allow an agnostic characterization of these physical interactions on a genome-wide scale. Single-nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci allows us to identify key long-range chromatin interactions in cis and trans involving these loci

Phenotype and Clinical Outcomes of Titin Cardiomyopathy.

BACKGROUND: Improved understanding of dilated cardiomyopathy (DCM) due to titin truncation (TTNtv) may help guide patient stratification. OBJECTIVES: The purpose of this study was to establish relationships among TTNtv genotype, cardiac phenotype, and outcomes in DCM. METHODS: In this prospective, observational cohort study, DCM patients underwent clinical evaluation, late gadolinium enhancement cardiovascular magnetic resonance, TTN sequencing, and adjudicated follow-up blinded to genotype for the primary composite endpoint of cardiovascular death, and major arrhythmic and major heart failure events. RESULTS: Of 716 subjects recruited (mean age 53.5 ± 14.3 years; 469 men [65.5%]; 577 [80.6%] New York Heart Association function class I/II), 83 (11.6%) had TTNtv. Patients with TTNtv were younger at enrollment (49.0 years vs. 54.1 years; p = 0.002) and had lower indexed left ventricular mass (5.1 g/m2 reduction; padjusted = 0.03) compared with patients without TTNtv. There was no difference in biventricular ejection fraction between TTNtv+/- groups. Overall, 78 of 604 patients (12.9%) met the primary endpoint (median follow-up 3.9 years; interquartile range: 2.0 to 5.8 years), including 9 of 71 patients with TTNtv (12.7%) and 69 of 533 (12.9%) without. There was no difference in the composite primary outcome of cardiovascular death, heart failure, or arrhythmic events, for patients with or without TTNtv (hazard ratio adjusted for primary endpoint: 0.92 [95% confidence interval: 0.45 to 1.87]; p = 0.82). CONCLUSIONS: In this large, prospective, genotype-phenotype study of ambulatory DCM patients, we show that prognostic factors for all-cause DCM also predict outcome in TTNtv DCM, and that TTNtv DCM does not appear to be associated with worse medium-term prognosis

Characterising the loss-of-function impact of 5' untranslated region variants in 15,708 individuals

Upstream open reading frames (uORFs) are tissue-specific cis-regulators of protein translation. Isolated reports have shown that variants that create or disrupt uORFs can cause disease. Here, in a systematic genome-wide study using 15,708 whole genome sequences, we show that variants that create new upstream start codons, and variants disrupting stop sites of existing uORFs, are under strong negative selection. This selection signal is significantly stronger for variants arising upstream of genes intolerant to loss-of-function variants. Furthermore, variants creating uORFs that overlap the coding sequence show signals of selection equivalent to coding missense variants. Finally, we identify specific genes where modification of uORFs likely represents an important disease mechanism, and report a novel uORF frameshift variant upstream of NF2 in neurofibromatosis. Our results highlight uORF-perturbing variants as an under-recognised functional class that contribute to penetrant human disease, and demonstrate the power of large-scale population sequencing data in studying non-coding variant classes. Upstream open reading frames (uORFs), located in 5' untranslated regions, are regulators of downstream protein translation. Here, Whiffin et al. use the genomes of 15,708 individuals in the Genome Aggregation Database (gnomAD) to systematically assess the deleteriousness of variants creating or disrupting uORFs.Peer reviewe

Phenotype and Clinical Outcomes of Titin Cardiomyopathy

BACKGROUND: Improved understanding of dilated cardiomyopathy (DCM) due to titin truncation (TTNtv) may help guide patient stratification. OBJECTIVES: The purpose of this study was to establish relationships among TTNtv genotype, cardiac phenotype, and outcomes in DCM. METHODS: In this prospective, observational cohort study, DCM patients underwent clinical evaluation, late gadolinium enhancement cardiovascular magnetic resonance, TTN sequencing, and adjudicated follow-up blinded to genotype for the primary composite endpoint of cardiovascular death, and major arrhythmic and major heart failure events. RESULTS: Of 716 subjects recruited (mean age 53.5 ± 14.3 years; 469 men [65.5%]; 577 [80.6%] New York Heart Association function class I/II), 83 (11.6%) had TTNtv. Patients with TTNtv were younger at enrollment (49.0 years vs. 54.1 years; p = 0.002) and had lower indexed left ventricular mass (5.1 g/m2 reduction; padjusted = 0.03) compared with patients without TTNtv. There was no difference in biventricular ejection fraction between TTNtv+/- groups. Overall, 78 of 604 patients (12.9%) met the primary endpoint (median follow-up 3.9 years; interquartile range: 2.0 to 5.8 years), including 9 of 71 patients with TTNtv (12.7%) and 69 of 533 (12.9%) without. There was no difference in the composite primary outcome of cardiovascular death, heart failure, or arrhythmic events, for patients with or without TTNtv (hazard ratio adjusted for primary endpoint: 0.92 [95% confidence interval: 0.45 to 1.87]; p = 0.82). CONCLUSIONS: In this large, prospective, genotype-phenotype study of ambulatory DCM patients, we show that prognostic factors for all-cause DCM also predict outcome in TTNtv DCM, and that TTNtv DCM does not appear to be associated with worse medium-term prognosis

Phenotypic expression and outcomes in individuals with rare genetic variants of hypertrophic cardiomyopathy

BACKGROUND: Hypertrophic cardiomyopathy (HCM) is caused by rare variants in sarcomere-encoding genes, but little is known about the clinical significance of these variants in the general population. OBJECTIVES: The goal of this study was to compare lifetime outcomes and cardiovascular phenotypes according to the presence of rare variants in sarcomere-encoding genes among middle-aged adults. METHODS: This study analyzed whole exome sequencing and cardiac magnetic resonance imaging in UK Biobank participants stratified according to sarcomere-encoding variant status. RESULTS: The prevalence of rare variants (allele frequency <0.00004) in HCM-associated sarcomere-encoding genes in 200,584 participants was 2.9% (n = 5,712; 1 in 35), and the prevalence of variants pathogenic or likely pathogenic for HCM (SARC-HCM-P/LP) was 0.25% (n = 493; 1 in 407). SARC-HCM-P/LP variants were associated with an increased risk of death or major adverse cardiac events compared with controls (hazard ratio: 1.69; 95% confidence interval [CI]: 1.38-2.07; P < 0.001), mainly due to heart failure endpoints (hazard ratio: 4.23; 95% CI: 3.07-5.83; P < 0.001). In 21,322 participants with both cardiac magnetic resonance imaging and whole exome sequencing, SARC-HCM-P/LP variants were associated with an asymmetric increase in left ventricular maximum wall thickness (10.9 ± 2.7 mm vs 9.4 ± 1.6 mm; P < 0.001), but hypertrophy (≥13 mm) was only present in 18.4% (n = 9 of 49; 95% CI: 9%-32%). SARC-HCM-P/LP variants were still associated with heart failure after adjustment for wall thickness (hazard ratio: 6.74; 95% CI: 2.43-18.7; P < 0.001). CONCLUSIONS: In this population of middle-aged adults, SARC-HCM-P/LP variants have low aggregate penetrance for overt HCM but are associated with an increased risk of adverse cardiovascular outcomes and an attenuated cardiomyopathic phenotype. Although absolute event rates are low, identification of these variants may enhance risk stratification beyond familial disease

Quantifying evidence toward pathogenicity for rare phenotypes: The case of succinate dehydrogenase genes, SDHB and SDHD.

PURPOSE: The weight of the evidence to attach to observation of a novel rare missense variant in SDHB or SDHD in individuals with the rare neuroendocrine tumors, pheochromocytomas and paragangliomas (PCC/PGL), is uncertain. METHODS: We compared the frequency of SDHB and SDHD very rare missense variants (VRMVs) in 6328 and 5847 cases of PCC/PGL, respectively, with that of population controls to generate a pan-gene VRMV likelihood ratio (LR). Via windowing analysis, we measured regional enrichments of VRMVs to calculate the domain-specific VRMV-LR (DS-VRMV-LR). We also calculated subphenotypic LRs for variant pathogenicity for various clinical, histologic, and molecular features. RESULTS: We estimated the pan-gene VRMV-LR to be 76.2 (54.8-105.9) for SDHB and 14.8 (8.7-25.0) for SDHD. Clustering analysis revealed an SDHB enriched region (ɑɑ 177-260, P = .001) for which the DS-VRMV-LR was 127.2 (64.9-249.4) and an SDHD enriched region (ɑɑ 70-114, P = .000003) for which the DS-VRMV-LR was 33.9 (14.8-77.8). Subphenotypic LRs exceeded 6 for invasive disease (SDHB), head-and-neck disease (SDHD), multiple tumors (SDHD), family history of PCC/PGL, loss of SDHB staining on immunohistochemistry, and succinate-to-fumarate ratio >97 (SDHB, SDHD). CONCLUSION: Using methodology generalizable to other gene-phenotype dyads, the LRs relating to rarity and phenotypic specificity for a single observation in PCC/PGL of a SDHB/SDHD VRMV can afford substantial evidence toward pathogenicity

A retrospective observational study of the relationship between single nucleotide polymorphisms associated with the risk of developing Colorectal cancer and survival

Background: There is variability in clinical outcome for patients with apparently the same stage colorectal cancer (CRC). Single nucleotide polymorphisms (SNPs) mapping to chromosomes 1q41, 3q26.2, 6p21, 8q23.3, 8q24.21, 10p14, 11q13, 11q23.1, 12q13.13, 14q22, 14q22.2, 15q13.3, 16q22.1, 18q21.1, 19q13.11, 20p12, 20p12.3, 20q13.33 and Xp22 have robustly been shown to be associated with the risk of developing CRC. Since germline variation can also influence patient outcome the relationship between these SNPs and patient survivorship from CRC was examined. Methods: All enrolled into the National Study of Colorectal Cancer Genetics (NSCCG) were genotyped for 1q41, 3q26.2, 6p21, 8q23.3, 8q24.21, 10p14, 11q13, 11q23.1, 12q13.13, 14q22, 14q22.2, 15q13.3, 16q22.1, 18q21.1, 19q13.11, 20p12, 20p12.3, 20q13.33 and xp22 SNPs. Linking this information to the National Cancer Data Repository allowed patient genotype to be related to survival. Results: The linked dataset consisted of 4,327 individuals. 14q22.22 genotype defined by the SNP rs4444235 showed a significant association with overall survival. Specifically, the C allele was associated with poorer observed survival (per allele hazard ratio 1.13, 95% confidence interval 1.05-1.22, P = 0.0015). Conclusion: The CRC susceptibility SNP rs4444235 also appears to exert an influence in modulating patient survival and warrants further evaluation as a potential prognostic marker

Moderate excess alcohol consumption and adverse cardiac remodelling in dilated cardiomyopathy

Objective: The effect of moderate excess alcohol consumption is widely debated and has not been well defined in dilated cardiomyopathy (DCM). There is need for a greater evidence base to help advise patients. We sought to evaluate the effect of moderate excess alcohol consumption on cardiovascular structure, function and outcomes in DCM. Methods: Prospective longitudinal observational cohort study. Patients with DCM (n=604) were evaluated for a history of moderate excess alcohol consumption (UK government guidelines; &gt;14 units/week for women, &gt;21 units/week for men) at cohort enrollment, had cardiovascular magnetic resonance and were followed up for the composite endpoint of cardiovascular death, heart failure and arrhythmic events. Patients meeting criteria for alcoholic cardiomyopathy were not recruited. Results: DCM patients with a history of moderate excess alcohol consumption (n=98, 16%) had lower biventricular function and increased chamber dilatation of the left ventricle, right ventricle and left atrium, as well as increased left ventricular hypertrophy compared to patients without moderate alcohol consumption. They were more likely to be male (alcohol excess group– n =92, 94% vs n =306, 61%, p=&lt;0.001). After adjustment for biological sex, moderate excess alcohol was not associated with adverse cardiac structure. There was no difference in mid-wall myocardial fibrosis between groups. Prior moderate excess alcohol consumption did not affect prognosis (HR 1.29, 0.73 to 2.26, p=0.38) during median follow up of 3.9 years. Conclusion: Dilated cardiomyopathy patients with moderate excess alcohol consumption have adverse cardiac structure and function at presentation but this is largely due to biological sex. Alcohol may contribute to sex-specific phenotypic differences in DCM. These findings help to inform lifestyle discussions for patients with dilated cardiomyopathy