Emma Ramsay

Germline mutations in RAD51D confer susceptibility to ovarian cancer

Recently, RAD51C mutations were identified in families with breast and ovarian cancer. This observation prompted us to investigate the role of RAD51D in cancer susceptibility. We identified eight inactivating RAD51D mutations in unrelated individuals from 911 breast-ovarian cancer families compared with one inactivating mutation identified in 1,060 controls (P = 0.01). The association found here was principally with ovarian cancer, with three mutations identified in the 59 pedigrees with three or more individuals with ovarian cancer (P = 0.0005). The relative risk of ovarian cancer for RAD51D mutation carriers was estimated to be 6.30 (95% CI 2.86–13.85, P = 4.8 × 10−6). By contrast, we estimated the relative risk of breast cancer to be 1.32 (95% CI 0.59–2.96, P = 0.50). These data indicate that RAD51D mutation testing may have clinical utility in individuals with ovarian cancer and their families. Moreover, we show that cells deficient in RAD51D are sensitive to treatment with a PARP inhibitor, suggesting a possible therapeutic approach for cancers arising in RAD51D mutation carriers.

Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case–control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10−5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10−4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10−9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

Germline RAD51C mutations confer susceptibility to ovarian cancer

To the Editor: In 2010, Meindl and colleagues proposed that germline RAD51C mutations confer high risk for breast and ovarian cancer, comparable to BRCA1 and BRCA2 mutations1,2. However, multiple follow-up studies have provided no supportive evidence that RAD51C mutations predispose to breast cancer3–12. Following the original report, we began investigating the role of other RAD51 paralogs in breast and ovarian cancer susceptibility. This led to our recent discovery that germline RAD51D mutations predispose to ovarian cancer13. We identified truncating RAD51D mutations in 8 of 911 familial breast-ovarian cancer pedigrees and 1 of 1,060 population controls. Our analysis of simultaneous association with both breast and ovarian cancer risk showed that RAD51D mutations confer a sixfold increased risk of ovarian cancer (relative risk (RR) = 6.30, 95% confidence interval (CI) = 2.86–13.85; P = 4.8 °— 10-6) but do not affect or cause only a small increase in breast cancer risk (RR = 1.32, 95% CI = 0.59– 2.96; P = 0.50). This result was supported by our analysis of 737 familial breast cancer pedigrees with no ovarian cancer, in which we detected no RAD51D mutations. These findings prompted us to reevaluate the role of RAD51C in cancer susceptibility. We sequenced the full coding region and intronexon boundaries of RAD51C in 1,132 probands from families with a history of ovarian cancer occurring with or without breast cancer, 272 individuals with ovarian cancer from a hospital-based unselected case series and 1,156 population-based controls (Supplementary Tables 1 and 2 and Supplementary Methods). We identified 12 mutations that result in premature protein truncation in cases compared to 1 such mutation in controls (P = 0.009) (Table 1 and Supplementary Fig. 1). Nine mutations were identified among the 1,132 familial cases, and there was a higher prevalence of mutations in families with multiple ovarian cancer cases: 4 mutations were detected in 311 families with 2 or more cases of ovarian cancer, and 2 mutations were detected in the 67 families with 3 or more cases of ovarian cancer. Three mutations were identified among the 272 individuals with ovarian cancer unselected for family history, suggesting that ~1% of ovarian cancer cases harbor germline RAD51C mutations. We also identified a total of 12 nonsynonymous RAD51C variants (Supplementary Table 3). Four variants were identified in cases and controls; only one, c.790G>A, encoding a p.Gly264Ser amino-acid change, showed any evidence of association with cancer (P = 0.02), consistent with other studies2,9,12. Of note, this variant is predicted to be benign by in silico analyses and has limited impact on RAD51C function2. The remaining eight nonsynonymous variants were each identified in a single individual; there was no significant difference in the overall frequency (P = 0.36), position or predicted functional effects of these variants between cases and controls (Supplementary Table 3). These data exemplify the inherent complexities of evaluating the clinical consequences of missense variants (outside simple Mendelian disorders) and underscore why non-truncating and truncating variants should be considered separately. Analyzing controls for specific rare variants detected in cases and concluding that their absence in controls is evidence of pathogenicity can result in over-interpretation of the data. Such findings confirm that the specific variant is rare but can seldom provide conclusive evidence of disease association. Full sequencing of the gene in both cases and controls is a more appropriate analysis, as it allows the spectrum of variants in cases and controls to be directly compared. Functional and conservation data can be useful in the evaluation of variants, but in vitro functional effects do not necessarily imply that the variant has clinical sequelae. Moreover, as we and others have shown (for example, in studies of the breast cancer susceptibility genes BRIP1 and ATM), such an assumption can result in incorrect attribution of pathogenicity14,15. Better information is provided when mutational and functional analyses are equally ascertained in both cases and controls. To estimate the risk associated with RAD51C mutations, we undertook modified segregation analysis, in which we simultaneously modeled the risks of ovarian and breast cancer and incorporated control data and information from the full pedigrees of mutation-positive and mutation-negative families (Supplementary Methods). The relative risk of ovarian cancer for RAD51C mutation carriers was estimated to be 5.88 (95% CI = 2.91–11.88; P = 7.65 × 10-7), which constitutes a >9% cumulative risk by age 80. In contrast, there was no evidence of an association with breast cancer (RR = 0.91, 95% CI = 0.45–1.86; P = 0.8). Thus, the cancer risk estimates for RAD51C mutations were similar to those estimated for RAD51D mutations13. These data are fully consistent with the results presented by Meindl et al. and provide a likely explanation for why Meindl et al. identified RAD51C mutations only in breast cancer cases that had relatives with ovarian cancer and not in 620 familial breast cancer pedigrees without ovarian cancer. As RAD51C

Mutations in the DNA methyltransferase gene, DNMT3A, cause an overgrowth syndrome with intellectual disability

Overgrowth disorders are a heterogeneous group of conditions characterized by increased growth parameters and other variable clinical features such as intellectual disability and facial dysmorphism. To identify new causes of human overgrowth, we performed exome sequencing in ten proband-parent trios and detected two de novo DNMT3A mutations. We identified 11 additional de novo mutations by sequencing DNMT3A in a further 142 individuals with overgrowth. The mutations alter residues in functional DNMT3A domains, and protein modeling suggests that they interfere with domain-domain interactions and histone binding. Similar mutations were not present in 1,000 UK population controls (13/152 cases versus 0/1,000 controls; P < 0.0001). Mutation carriers had a distinctive facial appearance, intellectual disability and greater height. DNMT3A encodes a DNA methyltransferase essential for establishing methylation during embryogenesis and is commonly somatically mutated in acute myeloid leukemia. Thus, DNMT3A joins an emerging group of epigenetic DNA- and histone-modifying genes associated with both developmental growth disorders and hematological malignancies.Overgrowth disorders are a heterogeneous group of conditions characterized by increased growth parameters and other variable clinical features such as intellectual disability and facial dysmorphism. To identify new causes of human overgrowth, we performed exome sequencing in ten proband-parent trios and detected two de novo DNMT3A mutations. We identified 11 additional de novo mutations by sequencing DNMT3A in a further 142 individuals with overgrowth. The mutations alter residues in functional DNMT3A domains, and protein modeling suggests that they interfere with domain-domain interactions and histone binding. Similar mutations were not present in 1,000 UK population controls (13/152 cases versus 0/1,000 controls; P < 0.0001). Mutation carriers had a distinctive facial appearance, intellectual disability and greater height. DNMT3A encodes a DNA methyltransferase essential for establishing methylation during embryogenesis and is commonly somatically mutated in acute myeloid leukemia. Thus, DNMT3A joins an emerging group of epigenetic DNA- and histone-modifying genes associated with both developmental growth disorders and hematological malignancies.

Germline mutations in the PAF1 complex gene CTR9 predispose to Wilms tumour

Wilms tumour is a childhood kidney cancer. Here we identify inactivating CTR9 mutations in 3 of 35 Wilms tumour families, through exome and Sanger sequencing. By contrast, no similar mutations are present in 1,000 population controls (P<0.0001). Each mutation segregates with Wilms tumour in the family and a second mutational event is present in available tumours. CTR9 is a key component of the polymerase-associated factor 1 complex which has multiple roles in RNA polymerase II regulation and is implicated in embryonic organogenesis and maintenance of embryonic stem cell pluripotency. These data establish CTR9 as a Wilms tumour predisposition gene and suggest it acts as a tumour suppressor gene.

CSN and CAVA: variant annotation tools for rapid, robust next-generation sequencing analysis in the clinical setting

BackgroundNext-generation sequencing (NGS) offers unprecedented opportunities to expand clinical genomics. It also presents challenges with respect to integration with data from other sequencing methods and historical data. Provision of consistent, clinically applicable variant annotation of NGS data has proved difficult, particularly of indels, an important variant class in clinical genomics. Annotation in relation to a reference genome sequence, the DNA strand of coding transcripts and potential alternative variant representations has not been well addressed. Here we present tools that address these challenges to provide rapid, standardized, clinically appropriate annotation of NGS data in line with existing clinical standards.MethodsWe developed a clinical sequencing nomenclature (CSN), a fixed variant annotation consistent with the principles of the Human Genome Variation Society (HGVS) guidelines, optimized for automated variant annotation of NGS data. To deliver high-throughput CSN annotation we created CAVA (Clinical Annotation of VAriants), a fast, lightweight tool designed for easy incorporation into NGS pipelines. CAVA allows transcript specification, appropriately accommodates the strand of a gene transcript and flags variants with alternative annotations to facilitate clinical interpretation and comparison with other datasets. We evaluated CAVA in exome data and a clinical BRCA1/BRCA2 gene testing pipeline.ResultsCAVA generated CSN calls for 10,313,034 variants in the ExAC database in 13.44 hours, and annotated the ICR1000 exome series in 6.5 hours. Evaluation of 731 different indels from a single individual revealed 92 % had alternative representations in left aligned and right aligned data. Annotation of left aligned data, as performed by many annotation tools, would thus give clinically discrepant annotation for the 339 (46 %) indels in genes transcribed from the forward DNA strand. By contrast, CAVA provides the correct clinical annotation for all indels. CAVA also flagged the 370 indels with alternative representations of a different functional class, which may profoundly influence clinical interpretation. CAVA annotation of 50 BRCA1/BRCA2 gene mutations from a clinical pipeline gave 100 % concordance with Sanger data; only 8/25 BRCA2 mutations were correctly clinically annotated by other tools.ConclusionsCAVA is a freely available tool that provides rapid, robust, high-throughput clinical annotation of NGS data, using a standardized clinical sequencing nomenclature.

Multi-stage genome-wide association study identifies new susceptibility locus for testicular germ cell tumour on chromosome 3q25

Recent genome-wide association studies (GWAS) and subsequent meta-analyses have identified over 25 SNPs at 18 loci, together accounting for >15% of the genetic susceptibility to testicular germ cell tumour (TGCT). To identify further common SNPs associated with TGCT, here we report a three-stage experiment, involving 4098 cases and 18 972 controls. Stage 1 comprised previously published GWAS analysis of 307 291 SNPs in 986 cases and 4946 controls. In Stage 2, we used previously published customised Illumina iSelect genotyping array (iCOGs) data across 694 SNPs in 1064 cases and 10 082 controls. Here, we report new genotyping of eight SNPs showing some evidence of association in combined analysis of Stage 1 and Stage 2 in an additional 2048 cases of TGCT and 3944 controls (Stage 3). Through fixed-effects meta-analysis across three stages, we identified a novel locus at 3q25.31 (rs1510272) demonstrating association with TGCT [per-allele odds ratio (OR) = 1.16, 95% confidence interval (CI) = 1.06-1.27; P = 1.2 × 10(-9)].

Mutations in the PP2A regulatory subunit B family genes PPP2R5B, PPP2R5C and PPP2R5D cause human overgrowth

Overgrowth syndromes comprise a group of heterogeneous disorders characterised by excessive growth parameters, often in association with intellectual disability. To identify new causes of human overgrowth, we have been undertaking trio-based exome sequencing studies in overgrowth patients and their unaffected parents. Prioritisation of functionally relevant genes with multiple unique de novo mutations revealed four mutations in protein phosphatase 2A (PP2A) regulatory subunit B family genes protein phosphatase 2, regulatory Subunit B’, beta (PPP2R5B); protein phosphatase 2, regulatory Subunit B’, gamma (PPP2R5C); and protein phosphatase 2, regulatory Subunit B’, delta (PPP2R5D). This observation in 3 related genes in 111 individuals with a similar phenotype is greatly in excess of the expected number, as determined from gene-specific de novo mutation rates (P = 1.43 × 10−10). Analysis of exome-sequencing data from a follow-up series of overgrowth probands identified a further pathogenic mutation, bringing the total number of affected individuals to 5. Heterozygotes shared similar phenotypic features including increased height, increased head circumference and intellectual disability. The mutations clustered within a region of nine amino acid residues in the aligned protein sequences (P = 1.6 × 10−5). We mapped the mutations onto the crystal structure of the PP2A holoenzyme complex to predict their molecular and functional consequences. These studies suggest that the mutations may affect substrate binding, thus perturbing the ability of PP2A to dephosphorylate particular protein substrates. PP2A is a major negative regulator of v-akt murine thymoma viral oncogene homolog 1 (AKT). Thus, our data further expand the list of genes encoding components of the phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K)/AKT signalling cascade that are disrupted in human overgrowth conditions.

Biallelic TRIP13 mutations predispose to Wilms tumor and chromosome missegregation

Through exome sequencing, we identified six individuals with biallelic loss-of-function mutations in TRIP13. All six developed Wilms tumor. Constitutional mosaic aneuploidies, microcephaly, developmental delay and seizures, which are features of mosaic variegated aneuploidy (MVA) syndrome, were more variably present. Through functional studies, we show that TRIP13-mutant patient cells have no detectable TRIP13 and have substantial impairment of the spindle assembly checkpoint (SAC), leading to a high rate of chromosome missegregation. Accurate segregation, as well as SAC proficiency, is rescued by restoring TRIP13 function. Individuals with biallelic TRIP13 or BUB1B mutations have a high risk of embryonal tumors, and here we show that their cells display severe SAC impairment. MVA due to biallelic CEP57 mutations, or of unknown cause, is not associated with embryonal tumors and cells from these individuals show minimal SAC deficiency. These data provide insights into the complex relationships between aneuploidy and carcinogenesis.

Mutations in Epigenetic Regulation Genes Are a Major Cause of Overgrowth with Intellectual Disability

To explore the genetic architecture of human overgrowth syndromes and human growth control, we performed experimental and bioinformatic analyses of 710 individuals with overgrowth (height and/or head circumference ≥+2 SD) and intellectual disability (OGID). We identified a causal mutation in 1 of 14 genes in 50% (353/710). This includes HIST1H1E, encoding histone H1.4, which has not been associated with a developmental disorder previously. The pathogenic HIST1H1E mutations are predicted to result in a product that is less effective in neutralizing negatively charged linker DNA because it has a reduced net charge, and in DNA binding and protein-protein interactions because key residues are truncated. Functional network analyses demonstrated that epigenetic regulation is a prominent biological process dysregulated in individuals with OGID. Mutations in six epigenetic regulation genes—NSD1, EZH2, DNMT3A, CHD8, HIST1H1E, and EED—accounted for 44% of individuals (311/710). There was significant overlap between the 14 genes involved in OGID and 611 genes in regions identified in GWASs to be associated with height (p = 6.84 × 10−8), suggesting that a common variation impacting function of genes involved in OGID influences height at a population level. Increased cellular growth is a hallmark of cancer and there was striking overlap between the genes involved in OGID and 260 somatically mutated cancer driver genes (p = 1.75 × 10−14). However, the mutation spectra of genes involved in OGID and cancer differ, suggesting complex genotype-phenotype relationships. These data reveal insights into the genetic control of human growth and demonstrate that exome sequencing in OGID has a high diagnostic yield.