Germline EPHB2 Receptor Variants in Familial Colorectal Cancer

Familial clustering of colorectal cancer occurs in 15–20% of cases, however recognized cancer syndromes explain only a small fraction of this disease. Thus, the genetic basis for the majority of hereditary colorectal cancer remains unknown. EPHB2 has recently been implicated as a candidate tumor suppressor gene in colorectal cancer. The aim of this study was to evaluate the contribution of EPHB2 to hereditary colorectal cancer. We screened for germline EPHB2 sequence variants in 116 population-based familial colorectal cancer cases by DNA sequencing. We then estimated the population frequencies and characterized the biological activities of the EPHB2 variants identified. Three novel nonsynonymous missense alterations were detected. Two of these variants (A438T and G787R) result in significant residue changes, while the third leads to a conservative substitution in the carboxy-terminal SAM domain (V945I). The former two variants were found once in the 116 cases, while the V945I variant was present in 2 cases. Genotyping of additional patients with colorectal cancer and control subjects revealed that A438T and G787R represent rare EPHB2 alleles. In vitro functional studies show that the G787R substitution, located in the kinase domain, causes impaired receptor kinase activity and is therefore pathogenic, whereas the A438T variant retains its receptor function and likely represents a neutral polymorphism. Tumor tissue from the G787R variant case manifested loss of heterozygosity, with loss of the wild-type allele, supporting a tumor suppressor role for EPHB2 in rare colorectal cancer cases. Rare germline EPHB2 variants may contribute to a small fraction of hereditary colorectal cancer

A large-scale meta-analysis to refine colorectal cancer risk estimates associated with MUTYH variants

Background: defective DNA repair has a causal role in hereditary colorectal cancer (CRC). Defects in the base excision repair gene MUTYH are responsible for MUTYH-associated polyposis and CRC predisposition as an autosomal recessive trait. Numerous reports have suggested MUTYH mono-allelic variants to be low penetrance risk alleles. We report a large collaborative meta-analysis to assess and refine CRC risk estimates associated with bi-allelic and mono-allelic MUTYH variants and investigate age and sex influence on risk. Methods: MUTYH genotype data were included from 20 565 cases and 15 524 controls. Three logistic regression models were tested: a crude model; adjusted for age and sex; adjusted for age, sex and study. Results: all three models produced very similar results. MUTYH bi-allelic carriers demonstrated a 28-fold increase in risk (95% confidence interval (CI): 6.95-115). Significant bi-allelic effects were also observed for G396D and Y179C/G396D compound heterozygotes and a marginal mono-allelic effect for variant Y179C (odds ratio (OR)=1.34; 95% CI: 1.00-1.80). A pooled meta-analysis of all published and unpublished datasets submitted showed bi-allelic effects for MUTYH, G396D and Y179C (OR=10.8, 95% CI: 5.02-23.2; OR=6.47, 95% CI: 2.33-18.0; OR=3.35, 95% CI: 1.14-9.89) and marginal mono-allelic effect for variants MUTYH (OR=1.16, 95% CI: 1.00-1.34) and Y179C alone (OR=1.34, 95% CI: 1.01-1.77). Conclusions: overall, this large study refines estimates of disease risk associated with mono-allelic and bi-allelic MUTYH carriers

Genetic variant predictors of gene expression provide new insight into risk of colorectal cancer

Genome-wide association studies have reported 56 independently associated colorectal cancer (CRC) risk variants, most of which are non-coding and believed to exert their effects by modulating gene expression. The computational method PrediXcan uses cis-regulatory variant predictors to impute expression and perform gene-level association tests in GWAS without directly measured transcriptomes. In this study, we used reference datasets from colon (n = 169) and whole blood (n = 922) transcriptomes to test CRC association with genetically determined expression levels in a genome-wide analysis of 12,186 cases and 14,718 controls. Three novel associations were discovered from colon transverse models at FDR ≤ 0.2 and further evaluated in an independent replication including 32,825 cases and 39,933 controls. After adjusting for multiple comparisons, we found statistically significant associations using colon transcriptome models with TRIM4 (discovery P = 2.2 × 10- 4, replication P = 0.01), and PYGL (discovery P = 2.3 × 10- 4, replication P = 6.7 × 10- 4). Interestingly, both genes encode proteins that influence redox homeostasis and are related to cellular metabolic reprogramming in tumors, implicating a novel CRC pathway linked to cell growth and proliferation. Defining CRC risk regions as one megabase up- and downstream of one of the 56 independent risk variants, we defined 44 non-overlapping CRC-risk regions. Among these risk regions, we identified genes associated with CRC (P < 0.05) in 34/44 CRC-risk regions. Importantly, CRC association was found for two genes in the previously reported 2q25 locus, CXCR1 and CXCR2, which are potential cancer therapeutic targets. These findings provide strong candidate genes to prioritize for subsequent laboratory follow-up of GWAS loci. This study is the first to implement PrediXcan in a large colorectal cancer study and findings highlight the utility of integrating transcriptome data in GWAS for discovery of, and biological insight into, risk loci

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.Swedish Research Council et al.Manuscrip

Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148

Genome wide association studies (GWAS) have mapped multiple independent cancer susceptibility loci to chr5p15.33. Here, we show that fine-mapping of pancreatic and testicular cancer GWAS within one of these loci (Region 2 in CLPTM1L) focuses the signal to nine highly correlated SNPs. Of these, rs36115365-C associated with increased pancreatic and testicular but decreased lung cancer and melanoma risk, and exhibited preferred protein-binding and enhanced regulatory activity. Transcriptional gene silencing of this regulatory element repressed TERT expression in an allele-specific manner. Proteomic analysis identifies allele-preferred binding of Zinc finger protein 148 (ZNF148) to rs36115365-C, further supported by binding of purified recombinant ZNF148. Knockdown of ZNF148 results in reduced TERT expression, telomerase activity and telomere length. Our results indicate that the association with chr5p15.33-Region 2 may be explained by rs36115365, a variant influencing TERT expression via ZNF148 in a manner consistent with elevated TERT in carriers of the C allele

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

BACKGROUND: Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders. METHODS: We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors. RESULTS: Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged. CONCLUSIONS: Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders

Exploration of Shared Genetic Architecture Between Subcortical Brain Volumes and Anorexia Nervosa

In MRI scans of patients with anorexia nervosa (AN), reductions in brain volume are often apparent. However, it is unknown whether such brain abnormalities are influenced by genetic determinants that partially overlap with those underlying AN. Here, we used a battery of methods (LD score regression, genetic risk scores, sign test, SNP effect concordance analysis, and Mendelian randomization) to investigate the genetic covariation between subcortical brain volumes and risk for AN based on summary measures retrieved from genome-wide association studies of regional brain volumes (ENIGMA consortium, n = 13,170) and genetic risk for AN (PGC-ED consortium, n = 14,477). Genetic correlations ranged from − 0.10 to 0.23 (all p > 0.05). There were some signs of an inverse concordance between greater thalamus volume and risk for AN (permuted p = 0.009, 95% CI: [0.005, 0.017]). A genetic variant in the vicinity of ZW10, a gene involved in cell division, and neurotransmitter and immune system relevant genes, in particular DRD2, was significantly associated with AN only after conditioning on its association with caudate volume (pFDR = 0.025). Another genetic variant linked to LRRC4C, important in axonal and synaptic development, reached significance after conditioning on hippocampal volume (pFDR = 0.021). In this comprehensive set of analyses and based on the largest available sample sizes to date, there was weak evidence for associations between risk for AN and risk for abnormal subcortical brain volumes at a global level (that is, common variant genetic architecture), but suggestive evidence for effects of single genetic markers. Highly powered multimodal brain- and disorder-related genome-wide studies are needed to further dissect the shared genetic influences on brain structure and risk for AN

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale(1-3). Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter(4); identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation(5,6); analyses timings and patterns of tumour evolution(7); describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity(8,9); and evaluates a range of more-specialized features of cancer genomes(8,10-18).Peer reviewe

Functional informed genome-wide interaction analysis of body mass index, diabetes and colorectal cancer risk

BACKGROUND: Body mass index (BMI) and diabetes are established risk factors for colorectal cancer (CRC), likely through perturbations in metabolic traits (e.g. insulin resistance and glucose homeostasis). Identification of interactions between variation in genes and these metabolic risk factors may identify novel biologic insights into CRC etiology. METHODS: To improve statistical power and interpretation for gene-environment interaction (G × E) testing, we tested genetic variants that regulate expression of a gene together for interaction with BMI (kg/m2 ) and diabetes on CRC risk among 26 017 cases and 20 692 controls. Each variant was weighted based on PrediXcan analysis of gene expression data from colon tissue generated in the Genotype-Tissue Expression Project for all genes with heritability ≥1%. We used a mixed-effects model to jointly measure the G × E interaction in a gene by partitioning the interactions into the predicted gene expression levels (fixed effects), and residual G × E effects (random effects). G × BMI analyses were stratified by sex as BMI-CRC associations differ by sex. We used false discovery rates to account for multiple comparisons and reported all results with FDR <0.2. RESULTS: Among 4839 genes tested, genetically predicted expressions of FOXA1 (P = 3.15 × 10-5 ), PSMC5 (P = 4.51 × 10-4 ) and CD33 (P = 2.71 × 10-4 ) modified the association of BMI on CRC risk for men; KIAA0753 (P = 2.29 × 10-5 ) and SCN1B (P = 2.76 × 10-4 ) modified the association of BMI on CRC risk for women; and PTPN2 modified the association between diabetes and CRC risk in both sexes (P = 2.31 × 10-5 ). CONCLUSIONS: Aggregating G × E interactions and incorporating functional information, we discovered novel genes that may interact with BMI and diabetes on CRC risk