Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

Defects of ciliogenesis have been implicated in a wide range of human phenotypes and play a crucial role in signal transduction and cell-cycle coordination. We used homozygosity mapping in two families with autosomal-recessive short-rib polydactyly syndrome Majewski type to identify mutations in NEK1 as an underlying cause of this lethal osteochondrodysplasia. NEK1 encodes a serine/threonine kinase with proposed function in DNA double-strand repair, neuronal development, and coordination of cell-cycle-associated ciliogenesis. We found that absence of functional full-length NEK1 severely reduces cilia number and alters ciliar morphology in vivo. We further substantiate a proposed digenic diallelic inheritance of ciliopathies by the identification of heterozygous mutations in NEK1 and DYNC2H1 in an additional family. Notably, these findings not only increase the broad spectrum of ciliar disorders, but suggest a correlation between the degree of defective microtubule or centriole elongation and organization and the severity of the resulting phenotype.

An increasing number of genes involved in chromatin structure and epigenetic regulation has been implicated in a variety of developmental disorders, often including intellectual disability. By trio exome sequencing and subsequent mutational screening we now identified two de novo frameshift mutations and one de novo missense mutation in CTCF in individuals with intellectual disability, microcephaly, and growth retardation. Furthermore, an individual with a larger deletion including CTCF was identified. CTCF (CCCTC-binding factor) is one of the most important chromatin organizers in vertebrates and is involved in various chromatin regulation processes such as higher order of chromatin organization, enhancer function, and maintenance of three-dimensional chromatin structure. Transcriptome analyses in all three individuals with point mutations revealed deregulation of genes involved in signal transduction and emphasized the role of CTCF in enhancer-driven expression of genes. Our findings indicate that haploinsufficiency of CTCF affects genomic interaction of enhancers and their regulated gene promoters that drive developmental processes and cognition.

RNA polymerase II interacts with various other complexes and factors to ensure correct initiation, elongation, and termination of mRNA transcription. One of these proteins is SR-related CTD-associated factor 4 (SCAF4), which is important for correct usage of polyA sites for mRNA termination. Using exome sequencing and international matchmaking, we identified nine likely pathogenic germline variants in SCAF4 including two splice-site and seven truncating variants, all residing in the N-terminal two thirds of the protein. Eight of these variants occurred de novo, and one was inherited. Affected individuals demonstrated a variable neurodevelopmental disorder characterized by mild intellectual disability, seizures, behavioral abnormalities, and various skeletal and structural anomalies. Paired-end RNA sequencing on blood lymphocytes of SCAF4-deficient individuals revealed a broad deregulation of more than 9,000 genes and significant differential splicing of more than 2,900 genes, indicating an important role of SCAF4 in mRNA processing. Knockdown of the SCAF4 ortholog CG4266 in the model organism Drosophila melanogaster resulted in impaired locomotor function, learning, and short-term memory. Furthermore, we observed an increased number of active zones in larval neuromuscular junctions, representing large glutamatergic synapses. These observations indicate a role of CG4266 in nervous system development and function and support the implication of SCAF4 in neurodevelopmental phenotypes. In summary, our data show that heterozygous, likely gene-disrupting variants in SCAF4 are causative for a variable neurodevelopmental disorder associated with impaired mRNA processing.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

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.

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.

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.

A recent meta-analysis of multiple genome-wide association and follow-up endometrial cancer case-control datasets identified a novel genetic risk locus for this disease at chromosome 14q32.33. To prioritize the functional SNP(s) and target gene(s) at this locus, we employed an in silico fine-mapping approach using genotyped and imputed SNP data for 6,608 endometrial cancer cases and 37,925 controls of European ancestry. Association and functional analyses provide evidence that the best candidate causal SNP is rs2494737. Multiple experimental analyses show that SNP rs2494737 maps to a silencer element located within AKT1, a member of the PI3K/AKT/MTOR intracellular signaling pathway activated in endometrial tumors. The rs2494737 risk A allele creates a YY1 transcription factor-binding site and abrogates the silencer activity in luciferase assays, an effect mimicked by transfection of YY1 siRNA. Our findings suggest YY1 is a positive regulator of AKT1, mediating the stimulatory effects of rs2494737 increasing endometrial cancer risk. Identification of an endometrial cancer risk allele within a member of the PI3K/AKT signaling pathway, more commonly activated in tumors by somatic alterations, raises the possibility that well tolerated inhibitors targeting this pathway could be candidates for evaluation as chemopreventive agents in individuals at high risk of developing endometrial cancer.

Heterozygous copy-number variants and SNPs of CNTNAP2 and NRXN1, two distantly related members of the neurexin superfamily, have been repeatedly associated with a wide spectrum of neuropsychiatric disorders, such as developmental language disorders, autism spectrum disorders, epilepsy, and schizophrenia. We now identified homozygous and compound-heterozygous deletions and mutations via molecular karyotyping and mutational screening in CNTNAP2 and NRXN1 in four patients with severe mental retardation (MR) and variable features, such as autistic behavior, epilepsy, and breathing anomalies, phenotypically overlapping with Pitt-Hopkins syndrome. With a frequency of at least 1% in our cohort of 179 patients, recessive defects in CNTNAP2 appear to significantly contribute to severe MR. Whereas the established synaptic role of NRXN1 suggests that synaptic defects contribute to the associated neuropsychiatric disorders and to severe MR as reported here, evidence for a synaptic role of the CNTNAP2-encoded protein CASPR2 has so far been lacking. Using Drosophila as a model, we now show that, as known for fly Nrx-I, the CASPR2 ortholog Nrx-IV might also localize to synapses. Overexpression of either protein can reorganize synaptic morphology and induce increased density of active zones, the synaptic domains of neurotransmitter release. Moreover, both Nrx-I and Nrx-IV determine the level of the presynaptic active-zone protein bruchpilot, indicating a possible common molecular mechanism in Nrx-I and Nrx-IV mutant conditions. We therefore propose that an analogous shared synaptic mechanism contributes to the similar clinical phenotypes resulting from defects in human NRXN1 and CNTNAP2.

Variants affecting the function of different subunits of the BAF chromatin-remodelling complex lead to various neurodevelopmental syndromes, including Coffin-Siris syndrome. Furthermore, variants in proteins containing PHD fingers, motifs recognizing specific histone tail modifications, have been associated with several neurological and developmental-delay disorders. Here, we report eight heterozygous de novo variants (one frameshift, two splice site, and five missense) in the gene encoding the BAF complex subunit double plant homeodomain finger 2 (DPF2). Affected individuals share common clinical features described in individuals with Coffin-Siris syndrome, including coarse facial features, global developmental delay, intellectual disability, speech impairment, and hypoplasia of fingernails and toenails. All variants occur within the highly conserved PHD1 and PHD2 motifs. Moreover, missense variants are situated close to zinc binding sites and are predicted to disrupt these sites. Pull-down assays of recombinant proteins and histone peptides revealed that a subset of the identified missense variants abolish or impaire DPF2 binding to unmodified and modified H3 histone tails. These results suggest an impairment of PHD finger structural integrity and cohesion and most likely an aberrant recognition of histone modifications. Furthermore, the overexpression of these variants in HEK293 and COS7 cell lines was associated with the formation of nuclear aggregates and the recruitment of both wild-type DPF2 and BRG1 to these aggregates. Expression analysis of truncating variants found in the affected individuals indicated that the aberrant transcripts escape nonsense-mediated decay. Altogether, we provide compelling evidence that de novo variants in DPF2 cause Coffin-Siris syndrome and propose a dominant-negative mechanism of pathogenicity.

Intellectual disability (ID) is a clinically and genetically heterogeneous common condition that remains etiologically unresolved in the majority of cases. Although several hundred diseased genes have been identified in X-linked, autosomal-recessive, or syndromic types of ID, the establishment of an etiological basis remains a difficult task in unspecific, sporadic cases. Just recently, de novo mutations in SYNGAP1, STXBP1, MEF2C, and GRIN2B were reported as relatively common causes of ID in such individuals. On the basis of a patient with severe ID and a 2.5 Mb microdeletion including ARID1B in chromosomal region 6q25, we performed mutational analysis in 887 unselected patients with unexplained ID. In this cohort, we found eight (0.9%) additional de novo nonsense or frameshift mutations predicted to cause haploinsufficiency. Our findings indicate that haploinsufficiency of ARID1B, a member of the SWI/SNF-A chromatin-remodeling complex, is a common cause of ID, and they add to the growing evidence that chromatin-remodeling defects are an important contributor to neurodevelopmental disorders.

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1.

The growth of an individual is deeply influenced by the regulation of cell growth and division, both of which also contribute to a wide variety of pathological conditions, including cancer, diabetes, and inflammation. To identify a major regulator of human growth, we performed positional cloning in an autosomal recessive type of profound short stature, anauxetic dysplasia. Homozygosity mapping led to the identification of novel mutations in the RMRP gene, which was previously known to cause two milder types of short stature with susceptibility to cancer, cartilage hair hypoplasia, and metaphyseal dysplasia without hypotrichosis. We show that different RMRP gene mutations lead to decreased cell growth by impairing ribosomal assembly and by altering cyclin-dependent cell cycle regulation. Clinical heterogeneity is explained by a correlation between the level and type of functional impairment in vitro and the severity of short stature or predisposition to cancer. Whereas the cartilage hair hypoplasia founder mutation affects both pathways intermediately, anauxetic dysplasia mutations do not affect B-cyclin messenger RNA (mRNA) levels but do severely incapacitate ribosomal assembly via defective endonucleolytic cleavage. Anauxetic dysplasia mutations thus lead to poor processing of ribosomal RNA while allowing normal mRNA processing and, therefore, genetically separate the different functions of RNase MRP.

Generalized pustular psoriasis (GPP) is a severe multi-systemic inflammatory disease characterized by neutrophilic pustulosis and triggered by pro-inflammatory IL-36 cytokines in skin. While 19%-41% of affected individuals harbor bi-allelic mutations in IL36RN, the genetic cause is not known in most cases. To identify and characterize new pathways involved in the pathogenesis of GPP, we performed whole-exome sequencing in 31 individuals with GPP and demonstrated effects of mutations in MPO encoding the neutrophilic enzyme myeloperoxidase (MPO). We discovered eight MPO mutations resulting in MPO -deficiency in neutrophils and monocytes. MPO mutations, primarily those resulting in complete MPO deficiency, cumulatively associated with GPP (p = 1.85E-08; OR = 6.47). The number of mutant MPO alleles significantly differed between 82 affected individuals and >4,900 control subjects (p = 1.04E-09); this effect was stronger when including IL36RN mutations (1.48E-13) and correlated with a younger age of onset (p = 0.0018). The activity of four proteases, previously implicated as activating enzymes of IL-36 precursors, correlated with MPO deficiency. Phorbol-myristate-acetate-induced formation of neutrophil extracellular traps (NETs) was reduced in affected cells (p = 0.015), and phagocytosis assays in MPO-deficient mice and human cells revealed altered neutrophil function and impaired clearance of neutrophils by monocytes (efferocytosis) allowing prolonged neutrophil persistence in inflammatory skin. MPO mutations contribute significantly to GPP's pathogenesis. We implicate MPO as an inflammatory modulator in humans that regulates protease activity and NET formation and modifies efferocytosis. Our findings indicate possible implications for the application of MPO inhibitors in cardiovascular diseases. MPO and affected pathways represent attractive targets for inducing resolution of inflammation in neutrophil-mediated skin diseases.

The link of chromatin remodeling to both neurodevelopment and cancer has recently been highlighted by the identification of mutations affecting BAF chromatin-remodeling components, such as ARID1B, in individuals with intellectual disability and cancer. However, the underlying molecular mechanism(s) remains unknown. Here, we show that ARID1B is a repressor of Wnt/β-catenin signaling. Through whole-transcriptome analysis, we find that in individuals with intellectual disability and ARID1B loss-of-function mutations, Wnt/β-catenin target genes are upregulated. Using cellular models of low and high Wnt/β-catenin activity, we demonstrate that knockdown of ARID1B activates Wnt/β-catenin target genes and Wnt/β-catenin-dependent transcriptional reporters in a β-catenin-dependent manner. Reciprocally, forced expression of ARID1B inhibits Wnt/β-catenin signaling downstream of the β-catenin destruction complex. Both endogenous and exogenous ARID1B associate with β-catenin and repress Wnt/β-catenin-mediated transcription through the BAF core subunit BRG1. Accordingly, mutations in ARID1B leading to partial or complete deletion of its BRG1-binding domain, as is often observed in intellectual disability and cancers, compromise association with β-catenin, and the resultant ARID1B mutant proteins fail to suppress Wnt/β-catenin signaling. Finally, knockdown of ARID1B in mouse neuroblastoma cells leads to neurite outgrowth through β-catenin. The data suggest that aberrations in chromatin-remodeling factors, such as ARID1B, might contribute to neurodevelopmental abnormalities and cancer through deregulation of developmental and oncogenic pathways, such as the Wnt/β-catenin signaling pathway.

Genome-wide association studies (GWASs) have revealed SNP rs889312 on 5q11.2 to be associated with breast cancer risk in women of European ancestry. In an attempt to identify the biologically relevant variants, we analyzed 909 genetic variants across 5q11.2 in 103,991 breast cancer individuals and control individuals from 52 studies in the Breast Cancer Association Consortium. Multiple logistic regression analyses identified three independent risk signals: the strongest associations were with 15 correlated variants (iCHAV1), where the minor allele of the best candidate, rs62355902, associated with significantly increased risks of both estrogen-receptor-positive (ER(+): odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.21-1.27, ptrend = 5.7 × 10(-44)) and estrogen-receptor-negative (ER(-): OR = 1.10, 95% CI = 1.05-1.15, ptrend = 3.0 × 10(-4)) tumors. After adjustment for rs62355902, we found evidence of association of a further 173 variants (iCHAV2) containing three subsets with a range of effects (the strongest was rs113317823 [pcond = 1.61 × 10(-5)]) and five variants composing iCHAV3 (lead rs11949391; ER(+): OR = 0.90, 95% CI = 0.87-0.93, pcond = 1.4 × 10(-4)). Twenty-six percent of the prioritized candidate variants coincided with four putative regulatory elements that interact with the MAP3K1 promoter through chromatin looping and affect MAP3K1 promoter activity. Functional analysis indicated that the cancer risk alleles of four candidates (rs74345699 and rs62355900 [iCHAV1], rs16886397 [iCHAV2a], and rs17432750 [iCHAV3]) increased MAP3K1 transcriptional activity. Chromatin immunoprecipitation analysis revealed diminished GATA3 binding to the minor (cancer-protective) allele of rs17432750, indicating a mechanism for its action. We propose that the cancer risk alleles act to increase MAP3K1 expression in vivo and might promote breast cancer cell survival.

The high-grade serous ovarian cancer (HGSOC) risk locus at chromosome 1p34.3 resides within a frequently amplified genomic region signifying the presence of an oncogene. Here, we integrate in silico variant-to-function analysis with functional studies to characterize the oncogenic potential of candidate genes in the 1p34.3 locus. Fine mapping of genome-wide association statistics identified candidate causal SNPs local to H3K27ac-demarcated enhancer regions that exhibit allele-specific binding for CTCF in HGSOC and normal fallopian tube secretory epithelium cells (FTSECs). SNP risk associations colocalized with eQTL for six genes (DNALI1, GNL2, RSPO1, SNIP1, MEAF6, and LINC01137) that are more highly expressed in carriers of the risk allele, and three (DNALI1, GNL2, and RSPO1) were upregulated in HGSOC compared to normal ovarian surface epithelium cells and/or FTSECs. Increased expression of GNL2 and MEAF6 was associated with shorter survival in HGSOC with 1p34.3 amplifications. Despite its activation of β-catenin signaling, RSPO1 overexpression exerted no effects on proliferation or colony formation in our study of ovarian cancer and FTSECs. Instead, GNL2, MEAF6, and SNIP1 silencing impaired in vitro ovarian cancer cell growth. Additionally, GNL2 silencing diminished xenograft tumor formation, whereas overexpression stimulated proliferation and colony formation in FTSECs. GNL2 influences 60S ribosomal subunit maturation and global protein synthesis in ovarian cancer and FTSECs, providing a potential mechanism of how GNL2 upregulation might promote ovarian cancer development and mediate genetic susceptibility of HGSOC.