Women with epithelial ovarian cancer (EOC) are usually treated with platinum/taxane therapy after cytoreductive surgery but there is considerable inter-individual variation in response. To identify germline single-nucleotide polymorphisms (SNPs) that contribute to variations in individual responses to chemotherapy, we carried out a multi-phase genome-wide association study (GWAS) in 1,244 women diagnosed with serous EOC who were treated with the same first-line chemotherapy, carboplatin and paclitaxel. We identified two SNPs (rs7874043 and rs72700653) in TTC39B (best P=7x10-5, HR=1.90, for rs7874043) associated with progression-free survival (PFS). Functional analyses show that both SNPs lie in a putative regulatory element (PRE) that physically interacts with the promoters of PSIP1, CCDC171 and an alternative promoter of TTC39B. The C allele of rs7874043 is associated with poor PFS and showed increased binding of the Sp1 transcription factor, which is critical for chromatin interactions with PSIP1. Silencing of PSIP1 significantly impaired DNA damage-induced Rad51 nuclear foci and reduced cell viability in ovarian cancer lines. PSIP1 (PC4 and SFRS1 Interacting Protein 1) is known to protect cells from stress-induced apoptosis, and high expression is associated with poor PFS in EOC patients. We therefore suggest that the minor allele of rs7874043 confers poor PFS by increasing PSIP1 expression.

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.

Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.

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.

Clinical genetic testing is commercially available for rs61764370, an inherited variant residing in a KRAS 3' UTR microRNA binding site, based on suggested associations with increased ovarian and breast cancer risk as well as with survival time. However, prior studies, emphasizing particular subgroups, were relatively small. Therefore, we comprehensively evaluated ovarian and breast cancer risks as well as clinical outcome associated with rs61764370.

Survival after a diagnosis of breast cancer varies considerably between patients, and some of this variation may be because of germline genetic variation. We aimed to identify genetic markers associated with breast cancer-specific survival.

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

To identify factors that are responsible for the phenotypic differences between transient chondrocytes within human osteophytes prone to endochondral ossification and permanent chondrocytes within articular cartilage persisting for decades.

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.

Spermatogenesis is an efficient and complex system of continuous cell differentiation. Previous studies investigating the transcriptomes of different cell populations in the testis relied either on sorting cells, cell depletion, or juvenile animals where not all stages of spermatogenesis have been completed. We present single-cell RNA sequencing (scRNA-Seq) data of 2,500 cells from the testes of two 8-week-old C57Bl/6J mice. Our dataset includes all spermatogenic stages from preleptotene to condensing spermatids as well as individual spermatogonia, Sertoli and Leydig cells. The data capture the full continuity of the meiotic and postmeiotic stages of spermatogenesis, and is thus ideally suited for marker discovery, network inference and similar analyses for which temporal ordering of differentiation processes can be exploited. Furthermore, it can serve as a reference for future studies involving single-cell RNA-Seq in mice where spermatogenesis is perturbed.

We evaluated genetics of hyperuricemia and gout, their interaction with kidney function and medication intake in chronic kidney disease (CKD) patients. Genome-wide association studies (GWAS) of urate and gout were performed in 4941 CKD patients in the German Chronic Kidney Disease (GCKD) study. Effect estimates of 26 known urate-associated population-based single nucleotide polymorphisms (SNPs) were examined. Interactions of urate-associated variants with urate-altering medications and clinical characteristics of gout were evaluated. Genome-wide significant associations with serum urate and gout were identified for known loci at SLC2A9 and ABCG2, but not for novel loci. Effects of the 26 known SNPs were of similar magnitude in CKD patients compared to population-based individuals, except for SNPs at ABCG2 that showed greater effects in CKD. Gene-medication interactions were not significant when accounting for multiple testing. Associations with gout in specific joints were significant for SLC2A9 rs12498742 in wrists and midfoot joints. Known genetic variants in SLC2A9 and ABCG2 were associated with urate and gout in a CKD cohort, with effect sizes for ABCG2 significantly greater in CKD compared to the general population. CKD patients are at high risk of gout due to reduced kidney function, diuretics intake and genetic predisposition, making treatment to target challenging.

We studied the association of the serum levels of the microRNA family members miR-320a/-b/-c with clinico-pathological data to assess their applicability as diagnostic biomarker in prostate cancer (PCa) patients. The levels of miR-320a/-b/-c in 3 groups were evaluated by qRT-PCR (145 patients with PCa, 31 patients with benign prostatic hyperplasia (BPH) and 19 healthy controls). The levels of the three family members of miR-320 were directly correlated within each group (P < 0.001), but they differed significantly among the three groups (P < 0.001). The serum levels of the miR-320 family members were significantly increased in older patients compared to younger patients (≤ 66 years vs. > 66 years, P ≤ 0.001). In addition, the levels of all three miR-320 family members were significantly different in patients with low tumor stage compared with those with high tumor stage (miR-320a: P = 0.034; miR-320b: P = 0.006; miR-320c: P = 0.007) and in patients with low serum PSA compared with those with high serum PSA (≤ 4 ng vs. > 4 ng; miR-320a: P = 0.003; miR-320b: P = 0.003; miR-320c: P = 0.006). The levels of these miRNAs were inversely correlated with serum PSA levels. Detection in the serum samples of PCa patients with or without PSA relapse revealed higher levels of miR-320a/-b/-c in the group without PSA relapse before/after radical prostatectomy than in that with PCa relapse. In summary, the differences among the PCa/BPH/healthy control groups with respect to miR-320a/-b/-c levels in conjunction with higher levels in patients without a PSA relapse than in those with a relapse suggest the diagnostic potential of these miRNA-320 family members in PCa patients.

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

The receptor activator of NF-κB (RANK) pathway is involved in bone health as well as breast cancer (BC) pathogenesis and progression. Whereas the therapeutic implication of this pathway is established for the treatment of osteoporosis and bone metastases, the application in adjuvant BC is currently investigated. As genetic variants in this pathway have been described to influence bone health, aim of this study was the prognostic relevance of genetic variants in RANK and RANKL. Single nucleotide polymorphisms in RANK(L) (rs1054016/rs1805034/rs35211496) were genotyped and analyzed with regard to bone metastasis-free survival (BMFS), disease-free survival, and overall survival for a retrospective cohort of 1251 patients. Cox proportional hazard models were built to examine the prognostic influence in addition to commonly established prognostic factors. The SNP rs1054016 seems to influence BMFS. Patients with two minor alleles had a more favorable prognosis than patients with at least one common allele (HR 0.37 (95% CI: 0.17, 0.84)), whereas other outcome parameters remained unaffected. rs1805034 and rs35211496 had no prognostic relevance. The effect of rs1054016(RANKL) adds to the evidence that the RANK pathway plays a role in BC pathogenesis and progression with respect to BMFS, emphasizing the connection between BC and bone health.

Invasive lobular breast cancer (ILC) accounts for 10-15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09-1.18), P = 6.0 × 10(-10); P-het for ILC vs IDC ER+ tumors = 1.8 × 10(-4)). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes.

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.

Obsessive-Compulsive Disorder (OCD) is a common and chronic disorder in which a person has uncontrollable, reoccurring thoughts and behaviours. It is a complex genetic condition and, in case of early onset (EO), the patients manifest a more severe phenotype, and an increased heritability. Large (>500 kb) copy number variations (CNVs) previously associated with autism and schizophrenia have been reported in OCD. Recently, rare CNVs smaller than 500 kb overlapping risk loci for other neurodevelopmental conditions have also been reported in OCD, stressing the importance of examining CNVs of any size range. The aim of this study was to further investigate the role of rare and small CNVs in the aetiology of EO-OCD.

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.