The microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain of Escherichia coli harboring the pks island that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown.

Aspirin use has been associated with reduced ovarian cancer risk, yet the underlying biological mechanisms are not fully understood. To gain mechanistic insights, we assessed the association between prediagnosis low and regular-dose aspirin use and gene expression profiles in ovarian tumors.

Copy number variants (CNVs) have been implicated in many complex diseases. We examined whether inherited CNVs were associated with overall survival among women with invasive epithelial ovarian cancer. Germline DNA from 1,056 cases (494 deceased, average of 3.7 years follow-up) was interrogated with the Illumina 610 quad genome-wide array containing, after quality control exclusions, 581,903 single nucleotide polymorphisms (SNPs) and 17,917 CNV probes. Comprehensive analysis capitalized upon the strengths of three complementary approaches to CNV classification. First, to identify small CNVs, single markers were evaluated and, where associated with survival, consecutive markers were combined. Two chromosomal regions were associated with survival using this approach (14q31.3 rs2274736 p = 1.59 × 10(-6), p = 0.001; 22q13.31 rs2285164 p = 4.01 × 10(-5), p = 0.009), but were not significant after multiple testing correction. Second, to identify large CNVs, genome-wide segmentation was conducted to characterize chromosomal gains and losses, and association with survival was evaluated by segment. Four regions were associated with survival (1q21.3 loss p = 0.005, 5p14.1 loss p = 0.004, 9p23 loss p = 0.002, and 15q22.31 gain p = 0.002); however, again, after correcting for multiple testing, no regions were statistically significant, and none were in common with the single marker approach. Finally, to evaluate associations with general amounts of copy number changes across the genome, we estimated CNV burden based on genome-wide numbers of gains and losses; no associations with survival were observed (p > 0.40). Although CNVs that were not well-covered by the Illumina 610 quad array merit investigation, these data suggest no association between inherited CNVs and survival after ovarian cancer.

Cdc20 is an activator of the anaphase-promoting complex/cyclosome that initiates anaphase onset by ordering the destruction of cyclin B1 and securin in metaphase. To study the physiological significance of Cdc20 in higher eukaryotes, we generated hypomorphic mice that express small amounts of this essential cell cycle regulator. In this study, we show that these mice are healthy and not prone to cancer despite substantial aneuploidy. Cdc20 hypomorphism causes chromatin bridging and chromosome misalignment, revealing a requirement for Cdc20 in efficient sister chromosome separation and chromosome-microtubule attachment. We find that cyclin B1 is newly synthesized during mitosis via cytoplasmic polyadenylation element-binding protein-dependent translation, causing its rapid accumulation between prometaphase and metaphase of Cdc20 hypomorphic cells. Anaphase onset is significantly delayed in Cdc20 hypomorphic cells but not when translation is inhibited during mitosis. These data reveal that Cdc20 is particularly rate limiting for cyclin B1 destruction because of regulated de novo synthesis of this cyclin after prometaphase onset.

One difficult question facing researchers is how to prioritize SNPs detected from genetic association studies for functional studies. Often a list of the top M SNPs is determined based on solely the p-value from an association analysis, where M is determined by financial/time constraints. For many studies of complex diseases, multiple analyses have been completed and integrating these multiple sets of results may be difficult. One may also wish to incorporate biological knowledge, such as whether the SNP is in the exon of a gene or a regulatory region, into the selection of markers to follow-up. In this manuscript, we propose a Bayesian latent variable model (BLVM) for incorporating "features" about a SNP to estimate a latent "quality score", with SNPs prioritized based on the posterior probability distribution of the rankings of these quality scores. We illustrate the method using data from an ovarian cancer genome-wide association study (GWAS). In addition to the application of the BLVM to the ovarian GWAS, we applied the BLVM to simulated data which mimics the setting involving the prioritization of markers across multiple GWAS for related diseases/traits. The top ranked SNP by BLVM for the ovarian GWAS, ranked 2(nd) and 7(th) based on p-values from analyses of all invasive and invasive serous cases. The top SNP based on serous case analysis p-value (which ranked 197(th) for invasive case analysis), was ranked 8(th) based on the posterior probability of being in the top 5 markers (0.13). In summary, the application of the BLVM allows for the systematic integration of multiple SNP "features" for the prioritization of loci for fine-mapping or functional studies, taking into account the uncertainty in ranking.

The 10q26 locus in the second intron of FGFR2 is the locus most strongly associated with estrogen-receptor-positive breast cancer in genome-wide association studies. We conducted fine-scale mapping in case-control studies genotyped with a custom chip (iCOGS), comprising 41 studies (n = 89,050) of European ancestry, 9 Asian ancestry studies (n = 13,983), and 2 African ancestry studies (n = 2,028) from the Breast Cancer Association Consortium. We identified three statistically independent risk signals within the locus. Within risk signals 1 and 3, genetic analysis identified five and two variants, respectively, highly correlated with the most strongly associated SNPs. By using a combination of genetic fine mapping, data on DNase hypersensitivity, and electrophoretic mobility shift assays to study protein-DNA binding, we identified rs35054928, rs2981578, and rs45631563 as putative functional SNPs. Chromatin immunoprecipitation showed that FOXA1 preferentially bound to the risk-associated allele (C) of rs2981578 and was able to recruit ERα to this site in an allele-specific manner, whereas E2F1 preferentially bound the risk variant of rs35054928. The risk alleles were preferentially found in open chromatin and bound by Ser5 phosphorylated RNA polymerase II, suggesting that the risk alleles are associated with changes in transcription. Chromatin conformation capture demonstrated that the risk region was able to interact with the promoter of FGFR2, the likely target gene of this risk region. A role for FOXA1 in mediating breast cancer susceptibility at this locus is consistent with the finding that the FGFR2 risk locus primarily predisposes to estrogen-receptor-positive disease.

Intestinal fibrosis is a long-term complication in inflammatory bowel diseases (IBD) that frequently results in functional damage, bowel obstruction, and surgery. Interleukin (IL) 36 is a group of cytokines in the IL1 family with inflammatory effects. We studied the expression of IL36 and its receptor, interleukin 1 receptor like 2 (IL1RL2 or IL36R) in the development of intestinal fibrosis in human tissues and mice.

Height is a heritable and highly heterogeneous trait. Short stature affects 3% of the population and in most cases is genetic in origin. After excluding known causes, 67% of affected individuals remain without diagnosis. To identify novel candidate genes for short stature, we performed exome sequencing in 254 unrelated families with short stature of unknown cause and identified variants in 63 candidate genes in 92 (36%) independent families. Based on systematic characterization of variants and functional analysis including expression in chondrocytes, we classified 13 genes as strong candidates. Whereas variants in at least two families were detected for all 13 candidates, two genes had variants in 6 (UBR4) and 8 (LAMA5) families, respectively. To facilitate their characterization, we established a clustered network of 1025 known growth and short stature genes, which yielded 29 significantly enriched clusters, including skeletal system development, appendage development, metabolic processes, and ciliopathy. Eleven of the candidate genes mapped to 21 of these clusters, including CPZ, EDEM3, FBRS, IFT81, KCND1, PLXNA3, RASA3, SLC7A8, UBR4, USP45, and ZFHX3. Fifty additional growth-related candidates we identified await confirmation in other affected families. Our study identifies Mendelian forms of growth retardation as an important component of idiopathic short stature.

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

High-risk mutations in several genes predispose to both colorectal cancer (CRC) and endometrial cancer (EC). We therefore hypothesised that some lower-risk genetic variants might also predispose to both CRC and EC. Using CRC and EC genome-wide association series, totalling 13,265 cancer cases and 40,245 controls, we found that the protective allele [G] at one previously-identified CRC polymorphism, rs2736100 near TERT, was associated with EC risk (odds ratio (OR) = 1.08, P = 0.000167); this polymorphism influences the risk of several other cancers. A further CRC polymorphism near TERC also showed evidence of association with EC (OR = 0.92; P = 0.03). Overall, however, there was no good evidence that the set of CRC polymorphisms was associated with EC risk, and neither of two previously-reported EC polymorphisms was associated with CRC risk. A combined analysis revealed one genome-wide significant polymorphism, rs3184504, on chromosome 12q24 (OR = 1.10, P = 7.23 × 10(-9)) with shared effects on CRC and EC risk. This polymorphism, a missense variant in the gene SH2B3, is also associated with haematological and autoimmune disorders, suggesting that it influences cancer risk through the immune response. Another polymorphism, rs12970291 near gene TSHZ1, was associated with both CRC and EC (OR = 1.26, P = 4.82 × 10(-8)), with the alleles showing opposite effects on the risks of the two cancers.

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

Joubert syndrome (JBTS) is a severe recessive neurodevelopmental ciliopathy which can affect several organ systems. Mutations in known JBTS genes account for approximately half of the cases. By homozygosity mapping and whole-exome sequencing, we identified a novel locus, JBTS23, with a homozygous splice site mutation in KIAA0586 (alias TALPID3), a known lethal ciliopathy locus in model organisms. Truncating KIAA0586 mutations were identified in two additional patients with JBTS. One mutation, c.428delG (p.Arg143Lysfs*4), is unexpectedly common in the general population and may be a major contributor to JBTS. We demonstrate KIAA0586 protein localization at the basal body in human and mouse photoreceptors, as is common for JBTS proteins, and also in pericentriolar locations. We show that loss of TALPID3 (KIAA0586) function in animal models causes abnormal tissue polarity, centrosome length and orientation, and centriolar satellites. We propose that JBTS and other ciliopathies may in part result from cell polarity defects.

Rhizomelic chondrodysplasia punctata (RCDP) is a group of disorders with overlapping clinical features including rhizomelia, chondrodysplasia punctata, coronal clefts, cervical dysplasia, congenital cataracts, profound postnatal growth retardation, severe intellectual disability, and seizures. Mutations in PEX7, GNPAT, and AGPS, all involved in the plasmalogen-biosynthesis pathway, have been described in individuals with RCDP. Here, we report the identification of mutations in another gene in plasmalogen biosynthesis, fatty acyl-CoA reductase 1 (FAR1), in two families affected by severe intellectual disability, early-onset epilepsy, microcephaly, congenital cataracts, growth retardation, and spasticity. Exome analyses revealed a homozygous in-frame indel mutation (c.495_507delinsT [p.Glu165_Pro169delinsAsp]) in two siblings from a consanguineous family and compound-heterozygous mutations (c.[787C>T];[1094A>G], p.[Arg263(∗)];[Asp365Gly]) in a third unrelated individual. FAR1 reduces fatty acids to their respective fatty alcohols for the plasmalogen-biosynthesis pathway. To assess the pathogenicity of the identified mutations, we transfected human embryonic kidney 293 cells with plasmids encoding FAR1 with either wild-type or mutated constructs and extracted the lipids from the cells. We screened the lipids with gas chromatography and mass spectrometry and found that all three mutations abolished the reductase activity of FAR1, given that no fatty alcohols could be detected. We also observed reduced plasmalogens in red blood cells in one individual to a range similar to that seen in individuals with RCDP, further supporting abolished FAR1 activity. We thus expand the spectrum of clinical features associated with defects in plasmalogen biosynthesis to include FAR1 deficiency as a cause of syndromic severe intellectual disability with cataracts, epilepsy, and growth retardation but without rhizomelia.

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82-0.88]) and ER-negative (OR = 0.87 [0.82-0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91-0.95] and OR = 1.06 [1.03-1.09]) and ER-negative (OR = 0.95 [0.91-0.98] and OR = 1.08 [1.04-1.13]) disease. There was weaker evidence for iCHAV4, located 5' of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90-0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1-4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.

Previous studies have identified common germline variants nominally associated with breast cancer survival. These associations have not been widely replicated in further studies. The purpose of this study was to evaluate the association of previously reported SNPs with breast cancer-specific survival using data from a pooled analysis of eight breast cancer survival genome-wide association studies (GWAS) from the Breast Cancer Association Consortium.

A genome-wide association study identified LMO1, which encodes an LIM-domain-only transcriptional cofactor, as a neuroblastoma susceptibility gene that functions as an oncogene in high-risk neuroblastoma. Here we show that dβh promoter-mediated expression of LMO1 in zebrafish synergizes with MYCN to increase the proliferation of hyperplastic sympathoadrenal precursor cells, leading to a reduced latency and increased penetrance of neuroblastomagenesis. The transgenic expression of LMO1 also promoted hematogenous dissemination and distant metastasis, which was linked to neuroblastoma cell invasion and migration, and elevated expression levels of genes affecting tumor cell-extracellular matrix interaction, including loxl3, itga2b, itga3, and itga5. Our results provide in vivo validation of LMO1 as an important oncogene that promotes neuroblastoma initiation, progression, and widespread metastatic dissemination.

High throughput sequencing has greatly advanced disease gene identification, especially in heterogeneous entities. Despite falling costs this is still an expensive and laborious technique, particularly when studying large cohorts. To address this problem we applied Exome Pool-Seq as an economic and fast screening technology in neurodevelopmental disorders (NDDs). Sequencing of 96 individuals can be performed in eight pools of 12 samples on less than one Illumina sequencer lane. In a pilot study with 96 cases we identified 27 variants, likely or possibly affecting function. Twenty five of these were identified in 923 established NDD genes (based on SysID database, status November 2016) (ACTB, AHDC1, ANKRD11, ATP6V1B2, ATRX, CASK, CHD8, GNAS, IFIH1, KCNQ2, KMT2A, KRAS, MAOA, MED12, MED13L, RIT1, SETD5, SIN3A, TCF4, TRAPPC11, TUBA1A, WAC, ZBTB18, ZMYND11), two in 543 (SysID) candidate genes (ZNF292, BPTF), and additionally a de novo loss-of-function variant in LRRC7, not previously implicated in NDDs. Most of them were confirmed to be de novo, but we also identified X-linked or autosomal-dominantly or autosomal-recessively inherited variants. With a detection rate of 28%, Exome Pool-Seq achieves comparable results to individual exome analyses but reduces costs by >85%. Compared with other large scale approaches using Molecular Inversion Probes (MIP) or gene panels, it allows flexible re-analysis of data. Exome Pool-Seq is thus well suited for large-scale, cost-efficient and flexible screening in characterized but heterogeneous entities like NDDs.

Hematotoxicity is one of the major side effects of chemotherapy. The aim of this study was to examine the association between single nucleotide polymorphisms (SNPs) and hematotoxicity in breast cancer patients in a subset of patients of the SUCCESS prospective phase III chemotherapy study. All patients (n = 1678) received three cycles of 5-fluorouracil, epirubicin, and cyclophosphamide (FEC) followed by three cycles of docetaxel or docetaxel/gemcitabine, depending on randomization. Germline DNA was genotyped for 246 SNPs selected from a previous genome-wide association study (GWAS) in a panel of lymphoblastoid cell lines, with gemcitabine toxicity as the phenotype. All SNPs were tested for their value in predicting grade 3 or 4 neutropenic or leukopenic events (NLEs). Their prognostic value in relation to overall survival and disease-free survival was also tested. None of the SNPs was found to be predictive for NLEs during treatment with docetaxel/gemcitabine. Two SNPs in and close to the PIGB gene significantly improved the prediction of NLEs after FEC, in addition to the factors of age and body surface area. The top SNP (rs12050587) had an odds ratio of 1.38 per minor allele (95% confidence interval, 1.17 to 1.62). No associations were identified for predicting disease-free or overall survival. Genetic variance in the PIGB gene may play a role in determining interindividual differences in relation to hematotoxicity after FEC chemotherapy.

Cancer progresses through distinct stages, and mouse models recapitulating traits of this progression are frequently used to explore genetic, morphological, and pharmacological aspects of tumor development. To complement genomic investigations of this process, we here quantify phosphoproteomic changes in skin cancer development using the SILAC mouse technology coupled to high-resolution mass spectrometry. We distill protein expression signatures from our data that distinguish between skin cancer stages. A distinct phosphoproteome of the two stages of cancer progression is identified that correlates with perturbed cell growth and implicates cell adhesion as a major driver of malignancy. Importantly, integrated analysis of phosphoproteomic data and prediction of kinase activity revealed PAK4-PKC/SRC network to be highly deregulated in SCC but not in papilloma. This detailed molecular picture, both at the proteome and phosphoproteome level, will prove useful for the study of mechanisms of tumor progression.