Nicola Whiffin

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.

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.

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.

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.

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.

Using high-resolution variant frequencies to empower clinical genome interpretation

PurposeWhole-exome and whole-genome sequencing have transformed the discovery of genetic variants that cause human Mendelian disease, but discriminating pathogenic from benign variants remains a daunting challenge. Rarity is recognized as a necessary, although not sufficient, criterion for pathogenicity, but frequency cutoffs used in Mendelian analysis are often arbitrary and overly lenient. Recent very large reference datasets, such as the Exome Aggregation Consortium (ExAC), provide an unprecedented opportunity to obtain robust frequency estimates even for very rare variants.MethodsWe present a statistical framework for the frequency-based filtering of candidate disease-causing variants, accounting for disease prevalence, genetic and allelic heterogeneity, inheritance mode, penetrance, and sampling variance in reference datasets.ResultsUsing the example of cardiomyopathy, we show that our approach reduces by two-thirds the number of candidate variants under consideration in the average exome, without removing true pathogenic variants (false-positive rate<0.001).ConclusionWe outline a statistically robust framework for assessing whether a variant is “too common” to be causative for a Mendelian disorder of interest. We present precomputed allele frequency cutoffs for all variants in the ExAC dataset.

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.

The TERT variant rs2736100 is associated with colorectal cancer risk

Background:Polymorphic variation at the 5p15.33 (TERT–CLPTM1L) locus is associated with the risk of many cancers but a relationship with colorectal cancer (CRC) risk has yet to be defined.Methods:We used data from six genome-wide association studies (GWAS) of CRC, linkage disequilibrium mapping and imputation, to examine the relationship between 73 single-nucleotide polymorphisms at 5p15.33 and CRC risk in detail.Results:rs2736100, which localises to intron 2 of TERT, provided the strongest evidence of an association with CRC (P=2.28 × 10−4). The association was also shown in an independent series of 10 047 CRC cases and 6918 controls (P=0.02). A meta-analysis of all seven studies (totalling 16 039 cases, 16 430 controls) provided increased evidence of association (P=2.49 × 10−5; per allele odds ratio=1.07). The association of rs2736100 on CRC risk was shown to be independent of 15 low-penetrance variants previously identified.Conclusion:The rs2736100 association demonstrates an influence of variation at 5p15.33 on CRC risk and further evidence that the 5p15.33 (TERT–CLPTM1L) locus has pleiotropic effects (reflecting generic or lineage-specific effects) on cancer risk.

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.

Adaptation and validation of the ACMG/AMP variant classification framework for MYH7 -associated inherited cardiomyopathies: recommendations by ClinGen’s Inherited Cardiomyopathy Expert Panel

PurposeIntegrating genomic sequencing in clinical care requires standardization of variant interpretation practices. The Clinical Genome Resource has established expert panels to adapt the American College of Medical Genetics and Genomics/Association for Molecular Pathology classification framework for specific genes and diseases. The Cardiomyopathy Expert Panel selected MYH7, a key contributor to inherited cardiomyopathies, as a pilot gene to develop a broadly applicable approach.MethodsExpert revisions were tested with 60 variants using a structured double review by pairs of clinical and diagnostic laboratory experts. Final consensus rules were established via iterative discussions.ResultsAdjustments represented disease-/gene-informed specifications (12) or strength adjustments of existing rules (5). Nine rules were deemed not applicable. Key specifications included quantitative frameworks for minor allele frequency thresholds, the use of segregation data, and a semiquantitative approach to counting multiple independent variant occurrences where fully controlled case-control studies are lacking. Initial inter-expert classification concordance was 93%. Internal data from participating diagnostic laboratories changed the classification of 20% of the variants (n = 12), highlighting the critical importance of data sharing.ConclusionThese adapted rules provide increased specificity for use in MYH7-associated disorders in combination with expert review and clinical judgment and serve as a stepping stone for genes and disorders with similar genetic and clinical characteristics.