Variant Location Is a Novel Risk Factor for Individuals With Arrhythmogenic Cardiomyopathy Due to a Desmoplakin (DSP) Truncating Variant.

BACKGROUND: Truncating variants in desmoplakin (DSPtv) are an important cause of arrhythmogenic cardiomyopathy; however the genetic architecture and genotype-specific risk factors are incompletely understood. We evaluated phenotype, risk factors for ventricular arrhythmias, and underlying genetics of DSPtv cardiomyopathy. METHODS: Individuals with DSPtv and any cardiac phenotype, and their gene-positive family members were included from multiple international centers. Clinical data and family history information were collected. Event-free survival from ventricular arrhythmia was assessed. Variant location was compared between cases and controls, and literature review of reported DSPtv performed. RESULTS: There were 98 probands and 72 family members (mean age at diagnosis 43±8 years, 59% women) with a DSPtv, of which 146 were considered clinically affected. Ventricular arrhythmia (sudden cardiac arrest, sustained ventricular tachycardia, appropriate implantable cardioverter defibrillator therapy) occurred in 56 (33%) individuals. DSPtv location and proband status were independent risk factors for ventricular arrhythmia. Further, gene region was important with variants in cases (cohort n=98; Clinvar n=167) more likely to occur in the regions resulting in nonsense mediated decay of both major DSP isoforms, compared with n=124 genome aggregation database control variants (148 [83.6%] versus 29 [16.4%]; P<0.0001). CONCLUSIONS: In the largest series of individuals with DSPtv, we demonstrate that variant location is a novel risk factor for ventricular arrhythmia, can inform variant interpretation, and provide critical insights to allow for precision-based clinical management.

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.

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.

Genetic Variants Associated With Cancer Therapy-Induced Cardiomyopathy.

BACKGROUND: Cancer therapy-induced cardiomyopathy (CCM) is associated with cumulative drug exposures and preexisting cardiovascular disorders. These parameters incompletely account for substantial interindividual susceptibility to CCM. We hypothesized that rare variants in cardiomyopathy genes contribute to CCM. METHODS: We studied 213 patients with CCM from 3 cohorts: retrospectively recruited adults with diverse cancers (n=99), prospectively phenotyped adults with breast cancer (n=73), and prospectively phenotyped children with acute myeloid leukemia (n=41). Cardiomyopathy genes, including 9 prespecified genes, were sequenced. The prevalence of rare variants was compared between CCM cohorts and The Cancer Genome Atlas participants (n=2053), healthy volunteers (n=445), and an ancestry-matched reference population. Clinical characteristics and outcomes were assessed and stratified by genotypes. A prevalent CCM genotype was modeled in anthracycline-treated mice. RESULTS: CCM was diagnosed 0.4 to 9 years after chemotherapy; 90% of these patients received anthracyclines. Adult patients with CCM had cardiovascular risk factors similar to the US population. Among 9 prioritized genes, patients with CCM had more rare protein-altering variants than comparative cohorts ( P≤1.98e-04). Titin-truncating variants (TTNtvs) predominated, occurring in 7.5% of patients with CCM versus 1.1% of The Cancer Genome Atlas participants ( P=7.36e-08), 0.7% of healthy volunteers ( P=3.42e-06), and 0.6% of the reference population ( P=5.87e-14). Adult patients who had CCM with TTNtvs experienced more heart failure and atrial fibrillation ( P=0.003) and impaired myocardial recovery ( P=0.03) than those without. Consistent with human data, anthracycline-treated TTNtv mice and isolated TTNtv cardiomyocytes showed sustained contractile dysfunction unlike wild-type ( P=0.0004 and P<0.002, respectively). CONCLUSIONS: Unrecognized rare variants in cardiomyopathy-associated genes, particularly TTNtvs, increased the risk for CCM in children and adults, and adverse cardiac events in adults. Genotype, along with cumulative chemotherapy dosage and traditional cardiovascular risk factors, improves the identification of patients who have cancer at highest risk for CCM. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov . Unique identifiers: NCT01173341; AAML1031; NCT01371981.

Recurrent Coding Sequence Variation Explains Only A Small Fraction of the Genetic Architecture of Colorectal Cancer.

Whilst common genetic variation in many non-coding genomic regulatory regions are known to impart risk of colorectal cancer (CRC), much of the heritability of CRC remains unexplained. To examine the role of recurrent coding sequence variation in CRC aetiology, we genotyped 12,638 CRCs cases and 29,045 controls from six European populations. Single-variant analysis identified a coding variant (rs3184504) in SH2B3 (12q24) associated with CRC risk (OR = 1.08, P = 3.9 × 10(-7)), and novel damaging coding variants in 3 genes previously tagged by GWAS efforts; rs16888728 (8q24) in UTP23 (OR = 1.15, P = 1.4 × 10(-7)); rs6580742 and rs12303082 (12q13) in FAM186A (OR = 1.11, P = 1.2 × 10(-7) and OR = 1.09, P = 7.4 × 10(-8)); rs1129406 (12q13) in ATF1 (OR = 1.11, P = 8.3 × 10(-9)), all reaching exome-wide significance levels. Gene based tests identified associations between CRC and PCDHGA genes (P < 2.90 × 10(-6)). We found an excess of rare, damaging variants in base-excision (P = 2.4 × 10(-4)) and DNA mismatch repair genes (P = 6.1 × 10(-4)) consistent with a recessive mode of inheritance. This study comprehensively explores the contribution of coding sequence variation to CRC risk, identifying associations with coding variation in 4 genes and PCDHG gene cluster and several candidate recessive alleles. However, these findings suggest that recurrent, low-frequency coding variants account for a minority of the unexplained heritability of CRC.

A new GWAS and meta-analysis with 1000Genomes imputation identifies novel risk variants for colorectal cancer.

Genome-wide association studies (GWAS) of colorectal cancer (CRC) have identified 23 susceptibility loci thus far. Analyses of previously conducted GWAS indicate additional risk loci are yet to be discovered. To identify novel CRC susceptibility loci, we conducted a new GWAS and performed a meta-analysis with five published GWAS (totalling 7,577 cases and 9,979 controls of European ancestry), imputing genotypes utilising the 1000 Genomes Project. The combined analysis identified new, significant associations with CRC at 1p36.2 marked by rs72647484 (minor allele frequency [MAF] = 0.09) near CDC42 and WNT4 (P = 1.21 × 10(-8), odds ratio [OR] = 1.21 ) and at 16q24.1 marked by rs16941835 (MAF = 0.21, P = 5.06 × 10(-8); OR = 1.15) within the long non-coding RNA (lncRNA) RP11-58A18.1 and ~500 kb from the nearest coding gene FOXL1. Additionally we identified a promising association at 10p13 with rs10904849 intronic to CUBN (MAF = 0.32, P = 7.01 × 10(-8); OR = 1.14). These findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to CRC. Additionally, our analysis further demonstrates that imputation can be used to exploit GWAS data to identify novel disease-causing variants.

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.

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.

Putative cis-regulatory drivers in colorectal cancer.

The cis-regulatory effects responsible for cancer development have not been as extensively studied as the perturbations of the protein coding genome in tumorigenesis. To better characterize colorectal cancer (CRC) development we conducted an RNA-sequencing experiment of 103 matched tumour and normal colon mucosa samples from Danish CRC patients, 90 of which were germline-genotyped. By investigating allele-specific expression (ASE) we show that the germline genotypes remain important determinants of allelic gene expression in tumours. Using the changes in ASE in matched pairs of samples we discover 71 genes with excess of somatic cis-regulatory effects in CRC, suggesting a cancer driver role. We correlate genotypes and gene expression to identify expression quantitative trait loci (eQTLs) and find 1,693 and 948 eQTLs in normal samples and tumours, respectively. We estimate that 36% of the tumour eQTLs are exclusive to CRC and show that this specificity is partially driven by increased expression of specific transcription factors and changes in methylation patterns. We show that tumour-specific eQTLs are more enriched for low CRC genome-wide association study (GWAS) P values than shared eQTLs, which suggests that some of the GWAS variants are tumour specific regulatory variants. Importantly, tumour-specific eQTL genes also accumulate more somatic mutations when compared to the shared eQTL genes, raising the possibility that they constitute germline-derived cancer regulatory drivers. Collectively the integration of genome and the transcriptome reveals a substantial number of putative somatic and germline cis-regulatory cancer changes that may have a role in tumorigenesis.

Architecture of inherited susceptibility to colorectal cancer: a voyage of discovery.

This review looks back at five decades of research into genetic susceptibility to colorectal cancer (CRC) and the insights these studies have provided. Initial evidence of a genetic basis of CRC stems from epidemiological studies in the 1950s and is further provided by the existence of multiple dominant predisposition syndromes. Genetic linkage and positional cloning studies identified the first high-penetrance genes for CRC in the 1980s and 1990s. More recent genome-wide association studies have identified common low-penetrance susceptibility loci and provide support for a polygenic model of disease susceptibility. These observations suggest a high proportion of CRC may arise in a group of susceptible individuals as a consequence of the combined effects of common low-penetrance risk alleles and rare variants conferring moderate CRC risks. Despite these advances, however, currently identified loci explain only a small fraction of the estimated heritability to CRC. It is hoped that a new generation of sequencing projects will help explain this missing heritability.

Identification of susceptibility loci for colorectal cancer in a genome-wide meta-analysis.

To identify common variants influencing colorectal cancer (CRC) risk, we performed a meta-analysis of five genome-wide association studies, comprising 5626 cases and 7817 controls of European descent. We conducted replication of top ranked single nucleotide polymorphisms (SNPs) in additional series totalling 14 037 cases and 15 937 controls, identifying a new CRC risk locus at 10q24.2 [rs1035209; odds ratio (OR) = 1.13, P = 4.54 × 10(-11)]. We also performed meta-analysis of our studies, with previously published data, of several recently purported CRC risk loci. We failed to find convincing evidence for a previously reported genome-wide association at rs11903757 (2q32.3). Of the three additional loci for which evidence of an association in Europeans has been previously described we failed to show an association between rs59336 (12q24.21) and CRC risk. However, for the other two SNPs, our analyses demonstrated new, formally significant associations with CRC. These are rs3217810 intronic in CCND2 (12p13.32; OR = 1.19, P = 2.16 × 10(-10)) and rs10911251 near LAMC1 (1q25.3; OR = 1.09, P = 1.75 × 10(-8)). Additionally, we found some evidence to support a relationship between, rs647161, rs2423297 and rs10774214 and CRC risk originally identified in East Asians in our European datasets. Our findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to CRC.

Deciphering the genetic architecture of low-penetrance susceptibility to colorectal cancer.

Recent genome-wide association studies (GWASs) have identified common variants at 16 autosomal regions influencing the risk of developing colorectal cancer (CRC). To decipher the genetic basis of the association signals at these loci, we performed a meta-analysis of data from five GWASs, totalling 5626 cases and 7817 controls, using imputation to recover un-typed genotypes. To enhance our ability to discover low-frequency risk variants, in addition to using 1000 Genomes Project data as a reference panel, we made use of high-coverage sequencing data on 253 individuals, 199 with early-onset familial CRC. For 13 of the regions, it was possible to refine the association signal identifying a smaller region of interest likely to harbour the functional variant. Our analysis did not provide evidence that any of the associations at the 16 loci being a consequence of synthetic associations rather than linkage disequilibrium with a common risk variant.

Spatiotemporal organization of Aurora-B by APC/CCdh1 after mitosis coordinates cell spreading through FHOD1.

Spatiotemporal regulation of mitotic kinase activity underlies the extensive rearrangement of cellular components required for cell division. One highly dynamic mitotic kinase is Aurora-B (AurB), which has multiple roles defined by the changing localisation of the chromosome passenger complex (CPC) as cells progress through mitosis, including regulation of cytokinesis and abscission. Like other mitotic kinases, AurB is a target of the anaphase-promoting complex (APC/C) ubiquitin ligase during mitotic exit, but it is not known if APC/C-mediated destruction plays any specific role in controlling AurB activity. We have examined the contribution of the Cdh1 coactivator-associated APC/C(Cdh1) to the organization of AurB activity as cells exit mitosis and re-enter interphase. We report that APC/C(Cdh1)-dependent proteolysis restricts a cell-cortex-associated pool of active AurB in space and time. In early G1 phase this pool of AurB is found at protrusions associated with cell spreading. AurB retention at the cortex depends on a formin, FHOD1, critically required to organize the cytoskeleton after division. We identify AurB phosphorylation sites in FHOD1 and show that phosphomutant FHOD1 is impaired in post-mitotic assembly of oriented actin cables. We propose that Cdh1 contributes to spatiotemporal organization of AurB activity and that organization of FHOD1 activity by AurB contributes to daughter cell spreading after mitosis.

Common variation near CDKN1A, POLD3 and SHROOM2 influences colorectal cancer risk.

We performed a meta-analysis of five genome-wide association studies to identify common variants influencing colorectal cancer (CRC) risk comprising 8,682 cases and 9,649 controls. Replication analysis was performed in case-control sets totaling 21,096 cases and 19,555 controls. We identified three new CRC risk loci at 6p21 (rs1321311, near CDKN1A; P = 1.14 × 10(-10)), 11q13.4 (rs3824999, intronic to POLD3; P = 3.65 × 10(-10)) and Xp22.2 (rs5934683, near SHROOM2; P = 7.30 × 10(-10)) This brings the number of independent loci associated with CRC risk to 20 and provides further insight into the genetic architecture of inherited susceptibility to CRC.

Relationship between 16 susceptibility loci and colorectal cancer phenotype in 3146 patients.

Recent genome-wide association studies have identified single-nucleotide polymorphisms at 16 genetic loci associated with colorectal cancer risk: rs6691170 (1q41), rs10936599 (3q26.2), rs16892766 (8q23.3), rs6983267 (8q24.21), rs10795668 (10p14), rs3802842 (11q23.1), rs11169552 (12q13.13), rs4444235, rs1957636 (14q22.2), rs4779584 (15q13.3), rs9929218 (16q22.1), rs4939827 (18q21.1), rs10411210 (19q13.11), rs961253 and rs4813802 (20p12.3) and rs4925386 (20q13.33). In the present study, we examined whether these variants are preferentially associated with tumour subtype-tumour site, stage, degree of differentiation and microsatellite instability status-in 3146 patients. Several loci showed statistically significant associations with specific phenotypes notably rs6691170 and rs3802842 associated with microsatellite stable rectal disease; rs4779584, rs961253 and rs4813802 associated with microsatellite stable colonic disease and rs4444235 and rs4925386 with microsatellite instability colonic disease. These findings are consistent with pathogenic variants in loci differentially impacting on distinct morphogenetic pathways consistent with aetiologically different risk factors in the development of colorectal cancer.

Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer.

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.

MLH1-93G > A is a risk factor for MSI colorectal cancer.

The -93G > A (rs1800734) polymorphism within the core promoter region of the MutL homolog 1 (MLH1) gene has recently been proposed as a low penetrance variant for colorectal cancer (CRC). We evaluated the significance of rs1800734 on CRC risk by genotyping 10 409 CRC cases and 6965 controls. The per allele odds ratio (OR) for all CRC-associated MLH1-93G > A was 1.06 (P = 0.037). Using a subset of 3132 cases with known microsatellite instability (MSI) status, the risk was shown to be confined to microsatellite instability-high (MSI-H) CRC; OR = 1.39 (P = 1.45 × 10(-4)). A meta-analysis of our study and four smaller published studies (totalling 801 cases, 10 890 controls) provided for increased evidence of relationship between MLH1-93G > A and MSI-H CRC risk (P = 3.43 × 10(-12)). The impact of MLH1-93G > A on CRC risk was shown to be independent of the 14 low penetrance loci for CRC identified by recent genome-wide association studies. These data provide further evidence that MLH1-93G > A is a low-penetrance variant for CRC and support the proposition that MLH1-93G > A acts as marker for a somatic event defining a specific CRC subtype.