Nonsyndromic Cleft Lip and/or Palate (NSCLP) is regarded as a multifactorial condition in which clefting is an isolated phenotype, distinguished from the largely monogenic, syndromic forms which include clefts among a spectrum of phenotypes. Nonsyndromic clefting has been shown to arise through complex interactions between genetic and environmental factors. However, there is increasing evidence that the broad NSCLP classification may include a proportion of cases showing familial patterns of inheritance and contain highly penetrant deleterious variation in specific genes. Through exome sequencing of multi-case families ascertained in Bogota, Colombia, we identify 28 non-synonymous single nucleotide variants that are considered damaging by at least one predictive score. We discuss the functional impact of candidate variants identified. In one family we find a coding variant in the MSX1 gene which is predicted damaging by multiple scores. This variant is in exon 2, a highly conserved region of the gene. Previous sequencing has suggested that mutations in MSX1 may account for ~2% of NSCLP. Our analysis further supports evidence that a proportion of NSCLP cases arise through monogenic coding mutations, though further work is required to unravel the complex interplay of genetics and environment involved in facial clefting.

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

Genome wide association studies (GWAs) of breast cancer mortality have identified few potential associations. The concordance between these studies is unclear. In this study, we used a meta-analysis of two prognostic GWAs and a replication cohort to identify the strongest associations and to evaluate the loci suggested in previous studies. We attempt to identify those SNPs which could impact overall survival irrespective of the age of onset.

Many common diseases have a complex genetic basis in which large numbers of genetic variations combine with environmental factors to determine risk. However, quantifying such polygenic effects has been challenging. In order to address these difficulties we developed a global measure of the information content of an individual's genome relative to a reference population, which may be used to assess differences in global genome structure between cases and appropriate controls. Informally this measure, which we call relative genome information (RGI), quantifies the relative "disorder" of an individual's genome. In order to test its ability to predict disease risk we used RGI to compare single-nucleotide polymorphism genotypes from two independent samples of women with early-onset breast cancer with three independent sets of controls. We found that RGI was significantly elevated in both sets of breast cancer cases in comparison with all three sets of controls, with disease risk rising sharply with RGI. Furthermore, these differences are not due to associations with common variants at a small number of disease-associated loci, but rather are due to the combined associations of thousands of markers distributed throughout the genome. Our results indicate that the information content of an individual's genome may be used to measure the risk of a complex disease, and suggest that early-onset breast cancer has a strongly polygenic component.

Clonal proliferation in myeloproliferative neoplasms (MPN) is driven by somatic mutations in JAK2, CALR or MPL, but the contribution of inherited factors is poorly characterized. Using a three-stage genome-wide association study of 3,437 MPN cases and 10,083 controls, we identify two SNPs with genome-wide significance in JAK2(V617F)-negative MPN: rs12339666 (JAK2; meta-analysis P=1.27 × 10(-10)) and rs2201862 (MECOM; meta-analysis P=1.96 × 10(-9)). Two additional SNPs, rs2736100 (TERT) and rs9376092 (HBS1L/MYB), achieve genome-wide significance when including JAK2(V617F)-positive cases. rs9376092 has a stronger effect in JAK2(V617F)-negative cases with CALR and/or MPL mutations (Breslow-Day P=4.5 × 10(-7)), whereas in JAK2(V617F)-positive cases rs9376092 associates with essential thrombocythemia (ET) rather than polycythemia vera (allelic χ(2) P=7.3 × 10(-7)). Reduced MYB expression, previously linked to development of an ET-like disease in model systems, associates with rs9376092 in normal myeloid cells. These findings demonstrate that multiple germline variants predispose to MPN and link constitutional differences in MYB expression to disease phenotype.

Multiple genes have been implicated by association studies in altering inflammatory bowel disease (IBD) predisposition. Paediatric patients often manifest more extensive disease and a particularly severe disease course. It is likely that genetic predisposition plays a more substantial role in this group.

Recently, metric linkage disequilibrium (LD) maps that assign an LD unit (LDU) location for each marker have been developed (Maniatis et al. 2002). Here we present a multiple pairwise method for positional cloning by LD within a composite likelihood framework and investigate the operating characteristics of maps in physical units (kb) and LDU for two bodies of data (Daly et al. 2001; Jeffreys et al. 2001) on which current ideas of blocks are based. False-negative indications of a disease locus (type II error) were examined by selecting one single-nucleotide polymorphism (SNP) at a time as causal and taking its allelic count (0, 1, or 2, for the three genotypes) as a pseudophenotype, Y. By use of regression and correlation, association between every pseudophenotype and the allelic count of each SNP locus (X) was based on an adaptation of the Malecot model, which includes a parameter for location of the putative gene. By expressing locations in kb or LDU, greater power for localization was observed when the LDU map was fitted. The efficiency of the kb map, relative to the LDU map, to describe LD varied from a maximum of 0.87 to a minimum of 0.36, with a mean of 0.62. False-positive indications of a disease locus (type I error) were examined by simulating an unlinked causal SNP and the allele count was used as a pseudophenotype. The type I error was in good agreement with Wald's likelihood theorem for both metrics and all models that were tested. Unlike tests that select only the most significant marker, haplotype, or haploset, these methods are robust to large numbers of markers in a candidate region. Contrary to predictions from tagging SNPs that retain haplotype diversity, the sample with smaller size but greater SNP density gave less error. The locations of causal SNPs were estimated with the same precision in blocks and steps, suggesting that block definition may be less useful than anticipated for mapping a causal SNP. These results provide a guide to efficient positional cloning by SNPs and a benchmark against which the power of positional cloning by haplotype-based alternatives may be measured.

Missense variants in RNA-binding proteins (RBPs) underlie a spectrum of disease phenotypes, including amyotrophic lateral sclerosis, frontotemporal dementia, and inclusion body myopathy. Here, we present ten independent families with a severe, progressive muscular dystrophy, reminiscent of oculopharyngeal muscular dystrophy (OPMD) but of much earlier onset, caused by heterozygous frameshift variants in the RBP hnRNPA2/B1. All disease-causing frameshift mutations abolish the native stop codon and extend the reading frame, creating novel transcripts that escape nonsense-mediated decay and are translated to produce hnRNPA2/B1 protein with the same neomorphic C-terminal sequence. In contrast to previously reported disease-causing missense variants in HNRNPA2B1, these frameshift variants do not increase the propensity of hnRNPA2 protein to fibrillize. Rather, the frameshift variants have reduced affinity for the nuclear import receptor karyopherin β2, resulting in cytoplasmic accumulation of hnRNPA2 protein in cells and in animal models that recapitulate the human pathology. Thus, we expand the phenotypes associated with HNRNPA2B1 to include an early-onset form of OPMD caused by frameshift variants that alter its nucleocytoplasmic transport dynamics.

We have utilized a two-stage study design to search for SNPs associated with the survival of breast cancer patients treated with adjuvant chemotherapy. Our initial GWS data set consisted of 805 Finnish breast cancer cases (360 treated with adjuvant chemotherapy). The top 39 SNPs from this stage were analyzed in three independent data sets: iCOGS (n=6720 chemotherapy-treated cases), SUCCESS-A (n=3596), and POSH (n=518). Two SNPs were successfully validated: rs6500843 (any chemotherapy; per-allele HR 1.16, 95% C.I. 1.08-1.26, p=0.0001, p(adjusted)=0.0091), and rs11155012 (anthracycline therapy; per-allele HR 1.21, 95% C.I. 1.08-1.35, p=0.0010, p(adjusted)=0.0270). The SNP rs6500843 was found to specifically interact with adjuvant chemotherapy, independently of standard prognostic markers (p(interaction)=0.0009), with the rs6500843-GG genotype corresponding to the highest hazard among chemotherapy-treated cases (HR 1.47, 95% C.I. 1.20-1.80). Upon trans-eQTL analysis of public microarray data, the rs6500843 locus was found to associate with the expression of a group of genes involved in cell cycle control, notably AURKA, the expression of which also exhibited differential prognostic value between chemotherapy-treated and untreated cases in our analysis of microarray data. Based on previously published information, we propose that the eQTL genes may be connected to the rs6500843 locus via a RBFOX1-FOXM1 -mediated regulatory pathway.

No abstract available

Ovarian carcinosarcoma (OCS) is an aggressive and rare tumor type with limited treatment options. OCS is hypothesized to develop via the combination theory, with a single progenitor resulting in carcinomatous and sarcomatous components, or alternatively via the conversion theory, with the sarcomatous component developing from the carcinomatous component through epithelial-to-mesenchymal transition (EMT). In this study, we analyzed DNA variants from isolated carcinoma and sarcoma components to show that OCS from 18 women is monoclonal. RNA sequencing indicated that the carcinoma components were more mesenchymal when compared with pure epithelial ovarian carcinomas, supporting the conversion theory and suggesting that EMT is important in the formation of these tumors. Preclinical OCS models were used to test the efficacy of microtubule-targeting drugs, including eribulin, which has previously been shown to reverse EMT characteristics in breast cancers and induce differentiation in sarcomas. Vinorelbine and eribulin more effectively inhibited OCS growth than standard-of-care platinum-based chemotherapy, and treatment with eribulin reduced mesenchymal characteristics and N-MYC expression in OCS patient-derived xenografts. Eribulin treatment resulted in an accumulation of intracellular cholesterol in OCS cells, which triggered a downregulation of the mevalonate pathway and prevented further cholesterol biosynthesis. Finally, eribulin increased expression of genes related to immune activation and increased the intratumoral accumulation of CD8+ T cells, supporting exploration of immunotherapy combinations in the clinic. Together, these data indicate that EMT plays a key role in OCS tumorigenesis and support the conversion theory for OCS histogenesis. Targeting EMT using eribulin could help improve OCS patient outcomes.

Rhabdomyosarcomas (RMS) are predominantly paediatric sarcomas thought to originate from muscle precursor cells due to impaired myogenic differentiation. Despite intensive treatment, 5-year survival for patients with advanced disease remains low (< 30%), highlighting a need for novel therapies to improve outcomes. Differentiation therapeutics are agents that induce differentiation of cancer cells from malignant to benign. The histone methyltransferase, Enhancer of Zeste Homolog 2 (EZH2) suppresses normal skeletal muscle differentiation and is highly expressed in RMS tumours.

Identifying healthy carriers of germline pathogenic variants in high penetrance cancer susceptibility genes offers the potential for risk-reducing surgery. The NHS England National Genomic Test Directory offers germline and somatic testing to patients with certain cancers or rare and inherited diseases, or, in some cases, to their relatives. This review summarises current UK guidelines for risk-reducing surgical interventions available for individuals with no personal history of cancer, who are determined to carry germline pathogenic variants. An electronic literature search of NICE guidelines and PubMed citable articles was performed. NICE guidelines are available for bilateral mastectomy and are currently in development for risk-reducing bilateral salpingo-oophorectomy. Guidelines developed with affiliation to, or through relevant British Surgical Societies or international consensus, are available for risk-reducing hysterectomy, polypectomy, gastrectomy, and thyroidectomy. There is a disparity in the development and distribution of national guidelines for interventions amongst tumour types. Whilst we are focusing on UK guidelines, we anticipate they will be relevant much more generally and so of interest to a wider audience including where there are no national guidelines to refer to. We suggest that, as genetic testing becomes rapidly more accessible, guideline development for interventions should be more closely aligned to those for testing.

Linkage maps have been invaluable for the positional cloning of many genes involved in severe human diseases. Standard genetic linkage maps have been constructed for this purpose from the Centre d'Etude du Polymorphisme Humain and other panels, and have been widely used. Now that attention has shifted towards identifying genes predisposing to common disorders using linkage disequilibrium (LD) and maps of single nucleotide polymorphisms (SNPs), it is of interest to consider a standard LD map which is somewhat analogous to the corresponding map for linkage. We have constructed and evaluated a cosmopolitan LD map by combining samples from a small number of populations using published data from a 10-megabase region on chromosome 20. In support of a pilot study, which examined a number of small genomic regions with a lower density of markers, we have found that a cosmopolitan map, which serves all populations when appropriately scaled, recovers 91 to 95 per cent of the information within population-specific maps. Recombination hot spots appear to have a dominant role in shaping patterns of LD. The success of the cosmopolitan map might be attributed to the co-localisation of hot spots in all populations. Although there must be finer scale differences between populations due to other processes (mutation, drift, selection), the results suggest that a whole-genome standard LD map would indeed be a useful resource for disease gene mapping.

PLEXIND1 is upregulated in several cancers, including pancreatic ductal adenocarcinoma (PDAC). It is an established mediator of semaphorin signaling, and neuropilins are its known coreceptors. Herein, we report data to support the proposal that PLEXIND1 acts as a transforming growth factor beta (TGFβ) coreceptor, modulating cell growth through SMAD3 signaling. Our findings demonstrate that PLEXIND1 plays a pro-tumorigenic role in PDAC cells with oncogenic KRAS (KRASmut). We show in KRASmut PDAC cell lines (PANC-1, AsPC-1,4535) PLEXIND1 downregulation results in decreased cell viability (in vitro) and reduced tumor growth (in vivo). Conversely, PLEXIND1 acts as a tumor suppressor in the PDAC cell line (BxPC-3) with wild-type KRAS (KRASwt), as its reduced expression results in higher cell viability (in-vitro) and tumor growth (in vivo). Additionally, we demonstrate that PLEXIND1-mediated interactions can be selectively disrupted using a peptide based on its C-terminal sequence (a PDZ domain-binding motif), an outcome that may possess significant therapeutic implications. To our knowledge, this is the first report showing that (1) PLEXIND1 acts as a TGFβ coreceptor and mediates SMAD3 signaling, and (2) differential roles of PLEXIND1 in PDAC cell lines correlate with KRASmut and KRASwt status.

Focal adhesion kinase (FAK) localizes to focal adhesions and is overexpressed in many cancers. FAK can also translocate to the nucleus, where it binds to, and regulates, several transcription factors, including MBD2, p53 and IL-33, to control gene expression by unknown mechanisms. We have used ATAC-seq to reveal that FAK controls chromatin accessibility at a subset of regulated genes. Integration of ATAC-seq and RNA-seq data showed that FAK-dependent chromatin accessibility is linked to differential gene expression, including of the FAK-regulated cytokine and transcriptional regulator interleukin-33 (Il33), which controls anti-tumor immunity. Analysis of the accessibility peaks on the Il33 gene promoter/enhancer regions revealed sequences for several transcription factors, including ETS and AP-1 motifs, and we show that c-Jun, a component of AP-1, regulates Il33 gene expression by binding to its enhancer in a FAK kinase-dependent manner. This work provides the first demonstration that FAK controls transcription via chromatin accessibility, identifying a novel mechanism by which nuclear FAK regulates biologically important gene expression.

CXCR2 has been suggested to have both tumor-promoting and tumor-suppressive properties. Here we show that CXCR2 signaling is upregulated in human pancreatic cancer, predominantly in neutrophil/myeloid-derived suppressor cells, but rarely in tumor cells. Genetic ablation or inhibition of CXCR2 abrogated metastasis, but only inhibition slowed tumorigenesis. Depletion of neutrophils/myeloid-derived suppressor cells also suppressed metastasis suggesting a key role for CXCR2 in establishing and maintaining the metastatic niche. Importantly, loss or inhibition of CXCR2 improved T cell entry, and combined inhibition of CXCR2 and PD1 in mice with established disease significantly extended survival. We show that CXCR2 signaling in the myeloid compartment can promote pancreatic tumorigenesis and is required for pancreatic cancer metastasis, making it an excellent therapeutic target.

An understanding of linkage disequilibrium (LD) structures in the human genome underpins much of medical genetics and provides a basis for disease gene mapping and investigating biological mechanisms such as recombination and selection. Whole genome sequencing (WGS) provides the opportunity to determine LD structures at maximal resolution.

To identify genetic variants associated with breast cancer prognosis we conduct a meta-analysis of overall survival (OS) and disease-free survival (DFS) in 6042 patients from four cohorts. In young women, breast cancer is characterized by a higher incidence of adverse pathological features, unique gene expression profiles and worse survival, which may relate to germline variation. To explore this hypothesis, we also perform survival analysis in 2315 patients aged ≤ 40 years at diagnosis. Here, we identify two SNPs associated with early-onset DFS, rs715212 (P meta = 3.54 × 10-5) and rs10963755 (P meta = 3.91 × 10-4) in ADAMTSL1. The effect of these SNPs is independent of classical prognostic factors and there is no heterogeneity between cohorts. Most importantly, the association with rs715212 is noteworthy (FPRP <0.2) and approaches genome-wide significance in multivariable analysis (P multivariable = 5.37 × 10-8). Expression quantitative trait analysis provides tentative evidence that rs715212 may influence AREG expression (P eQTL = 0.035), although further functional studies are needed to confirm this association and determine a mechanism.