PHGDH-related microcephalic dwarfism in two fetuses: Expanding the phenotypical spectrum of L-serine biosynthesis defect

Defects in L-serine biosynthesis are a group of autosomal recessive diseases resulting in a wide phenotypic spectrum ranging from viable to lethal presentations and caused by variants in the three genes encoding the L-serine biosynthesis enzymes, PHGDH, PSAT1, and PSPH. Neu-Laxova syndrome (NLS) is the fetal form of this group, characterized by multiple congenital anomalies including severe intrauterine growth retardation, cutaneous lesions extending from ichthyosis to severe restrictive dermopathy with ectropion and eclabion, edema, microcephaly, central nervous system abnormalities, and flexion contractures. Here we report on two unrelated fetuses with an attenuated phenotype of NLS, that initially evoked Taybi-Linder syndrome. They carry biallelic pathogenic variants in the PHGDH gene. These observations expand the phenotypic continuum of L-serine biosynthesis defects, and illustrate the phenotypic overlap between NLS and microcephalic primordial dwarfism.

Deficiency of the minor spliceosome component U4atac snRNA secondarily results in ciliary defects in human and zebrafish

In the human genome, about 750 genes contain one intron excised by the minor spliceosome. This spliceosome comprises its own set of snRNAs, among which U4atac. Its noncoding gene, RNU4ATAC, has been found mutated in Taybi-Linder (TALS/microcephalic osteodysplastic primordial dwarfism type 1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes. These rare developmental disorders, whose physiopathological mechanisms remain unsolved, associate ante- and post-natal growth retardation, microcephaly, skeletal dysplasia, intellectual disability, retinal dystrophy, and immunodeficiency. Here, we report bi-allelic RNU4ATAC mutations in five patients presenting with traits suggestive of the Joubert syndrome (JBTS), a well-characterized ciliopathy. These patients also present with traits typical of TALS/RFMN/LWS, thus widening the clinical spectrum of RNU4ATAC-associated disorders and indicating ciliary dysfunction as a mechanism downstream of minor splicing defects. Intriguingly, all five patients carry the n.16G>A mutation, in the Stem II domain, either at the homozygous or compound heterozygous state. A gene ontology term enrichment analysis on minor intron-containing genes reveals that the cilium assembly process is over-represented, with no less than 86 cilium-related genes containing at least one minor intron, among which there are 23 ciliopathy-related genes. The link between RNU4ATAC mutations and ciliopathy traits is supported by alterations of primary cilium function in TALS and JBTS-like patient fibroblasts, as well as by u4atac zebrafish model, which exhibits ciliopathy-related phenotypes and ciliary defects. These phenotypes could be rescued by WT but not by pathogenic variants-carrying human U4atac. Altogether, our data indicate that alteration of cilium biogenesis is part of the physiopathological mechanisms of TALS/RFMN/LWS, secondarily to defects of minor intron splicing.

Mutations in the non-coding RNU4ATAC gene affect the homeostasis and function of the Integrator complex

Various genetic diseases associated with microcephaly and developmental defects are due to pathogenic variants in the U4atac small nuclear RNA (snRNA), a component of the minor spliceosome essential for the removal of U12-type introns from eukaryotic mRNAs. While it has been shown that a few RNU4ATAC mutations result in impaired binding of essential protein components, the molecular defects of the vast majority of variants are still unknown. Here, we used lymphoblastoid cells derived from RNU4ATAC compound heterozygous (g.108_126del;g.111G>A) twin patients with MOPD1 phenotypes to analyze the molecular consequences of the mutations on small nuclear ribonucleoproteins (snRNPs) formation and on splicing. We found that the U4atac108_126del mutant is unstable and that the U4atac111G>A mutant as well as the minor di- and tri-snRNPs are present at reduced levels. Our results also reveal the existence of 3'-extended snRNA transcripts in patients' cells. Moreover, we show that the mutant cells have alterations in splicing of INTS7 and INTS10 minor introns, contain lower levels of the INTS7 and INTS10 proteins and display changes in the assembly of Integrator subunits. Altogether, our results show that compound heterozygous g.108_126del;g.111G>A mutations induce splicing defects and affect the homeostasis and function of the Integrator complex.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

Diagnostic chest X-rays and breast cancer risk among women with a hereditary predisposition to breast cancer unexplained by a BRCA1 or BRCA2 mutation

Background

Diagnostic ionizing radiation is a risk factor for breast cancer (BC). BC risk increases with increased dose to the chest and decreases with increased age at exposure, with possible effect modification related to familial or genetic predisposition. While chest X-rays increase the BC risk of BRCA1/2 mutation carriers compared to non-carriers, little is known for women with a hereditary predisposition to BC but who tested negative for a BRCA1 or BRCA2 (BRCA1/2) mutation.

Methods

We evaluated the effect of chest X-rays from diagnostic medical procedures in a dataset composed of 1552 BC cases identified through French family cancer clinics and 1363 unrelated controls. Participants reported their history of X-ray exposures in a detailed questionnaire and were tested for 113 DNA repair genes. Logistic regression and multinomial logistic regression models were used to assess the association with BC.

Results

Chest X-ray exposure doubled BC risk. A 3% increased BC risk per additional exposure was observed. Being 20 years old or younger at first exposure or being exposed before first full-term pregnancy did not seem to modify this risk. Birth after 1960 or carrying a rare likely deleterious coding variant in a DNA repair gene other than BRCA1/2 modified the effect of chest X-ray exposure.

Conclusion

Ever/never chest X-ray exposure increases BC risk 2-fold regardless of age at first exposure and, by up to 5-fold when carrying 3 or more rare variants in a DNA repair gene.

Further studies are needed to evaluate other DNA repair genes or variants to identify those which could modify radiation sensitivity. Identification of subpopulations that are more or less susceptible to ionizing radiation is important and potentially clinically relevant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01456-1.

5' Region Large Genomic Rearrangements in the BRCA1 Gene in French Families: Identification of a Tandem Triplication and Nine Distinct Deletions with Five Recurrent Breakpoints

Simple Summary

Large genomic rearrangements in BRCA1 consisting of deletions/duplications of one or several exons are complex events, often occurring in the 5′ region. We characterized 10 events in 20 families: one large triplication classified as benign and nine large deletions classified as pathogenic. The breakpoint localization will certainly help to further understand the chromatin structure in regions sensitive to rearrangement.

Abstract

Background: Large genomic rearrangements (LGR) in BRCA1 consisting of deletions/duplications of one or several exons have been found throughout the gene with a large proportion occurring in the 5′ region from the promoter to exon 2. The aim of this study was to better characterize those LGR in French high-risk breast/ovarian cancer families. Methods: DNA from 20 families with one apparent duplication and nine deletions was analyzed with a dedicated comparative genomic hybridization (CGH) array, high-resolution BRCA1 Genomic Morse Codes analysis and Sanger sequencing. Results: The apparent duplication was in fact a tandem triplication of exons 1 and 2 and part of intron 2 of BRCA1, fully characterized here for the first time. We calculated a causality score with the multifactorial model from data obtained from six families, classifying this variant as benign. Among the nine deletions detected in this region, eight have never been identified. The breakpoints fell in six recurrent regions and could confirm some specific conformation of the chromatin. Conclusions: Taken together, our results firmly establish that the BRCA1 5′ region is a frequent site of different LGRs and highlight the importance of the segmental duplication and Alu sequences, particularly the very high homologous region, in the mechanism of a recombination event. This also confirmed that those events are not systematically deleterious.

Gene- and pathway-level analyses of iCOGS variants highlight novel signaling pathways underlying familial breast cancer susceptibility

Single‐nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome‐wide association studies involving mostly unselected population‐based case‐control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC‐prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP‐level analysis using index cases and controls, we performed pedigree‐based association tests to capture transmission information in the sibships. We also performed gene‐ and pathway‐level analyses to maximize the power to detect associations with lower‐frequency SNPs or those with modest effect sizes. While SNP‐level analyses identified 18 loci, gene‐level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the “index case‐control” analysis, we built pathway‐derived polygenic risk scores (PRS) and assessed their performance in the population‐based CECILE study and in a data set composed of GENESIS‐affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP‐level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models.

What's new?

Genetic studies have identified more than 180 single‐nucleotide polymorphisms (SNPs) associated with breast cancer susceptibility, but these studies are reaching their limits. Here, the authors evaluated SNPs in the iCOGS genotyping array using a multilevel approach, including single variant, gene, and pathway analyses. They measured the contribution of the SNPs to breast cancer in patients who have a sister with breast cancer but do not carry a BRCA1/2 mutation. They showed that a pathway‐derived polygenic risk score performed poorly in the general population, and that the best predictive model must include family history.

A de novo frameshift pathogenic variant in TBR1 identified in autism without intellectual disability

Background

In order to be able to provide accurate genetic counseling to patients with Autism Spectrum Disorder (ASD), it is crucial to identify correlations between heterogeneous phenotypes and genetic alterations. Among the hundreds of de novo pathogenic variants reported in ASD, single-nucleotide variations and small insertions/deletions were reported in TBR1. This gene encodes a transcription factor that plays a key role in brain development. Pathogenic variants in TBR1 are often associated with severe forms of ASD, including intellectual disability and language impairment.

Methods

Adults diagnosed with ASD but without intellectual disability (diagnosis of Asperger syndrome, according to the DSM-IV) took part in a genetic consultation encompassing metabolic assessments, a molecular karyotype and the screening of a panel of 268 genes involved in intellectual disability, ASD and epilepsy. In addition, the patient reported here went through a neuropsychological assessment, structural magnetic resonance imaging and magnetic resonance spectroscopy measurements.

Results

Here, we report the case of a young adult male who presents with a typical form of ASD. Importantly, this patient presents with no intellectual disability or language impairment, despite a de novo heterozygous frameshift pathogenic variant in TBR1, leading to an early premature termination codon (c.26del, p.(Pro9Leufs*12)).

Conclusion

Based on this case report, we discuss the role of TBR1 in general brain development, language development, intellectual disability and other symptoms of ASD. Providing a detailed clinical description of the individuals with such pathogenic variants should help to understand the genotype-phenotype relationships in ASD.

New insights into minor splicing-a transcriptomic analysis of cells derived from TALS patients

Minor intron splicing plays a central role in human embryonic development and survival. Indeed, biallelic mutations in RNU4ATAC, transcribed into the minor spliceosomal U4atac snRNA, are responsible for three rare autosomal recessive multimalformation disorders named Taybi-Linder (TALS/MOPD1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes, which associate numerous overlapping signs of varying severity. Although RNA-seq experiments have been conducted on a few RFMN patient cells, none have been performed in TALS, and more generally no in-depth transcriptomic analysis of the ∼700 human genes containing a minor (U12-type) intron had been published as yet. We thus sequenced RNA from cells derived from five skin, three amniotic fluid, and one blood biosamples obtained from seven unrelated TALS cases and from age- and sex-matched controls. This allowed us to describe for the first time the mRNA expression and splicing profile of genes containing U12-type introns, in the context of a functional minor spliceosome. Concerning RNU4ATAC-mutated patients, we show that as expected, they display distinct U12-type intron splicing profiles compared to controls, but that rather unexpectedly mRNA expression levels are mostly unchanged. Furthermore, although U12-type intron missplicing concerns most of the expressed U12 genes, the level of U12-type intron retention is surprisingly low in fibroblasts and amniocytes, and much more pronounced in blood cells. Interestingly, we found several occurrences of introns that can be spliced using either U2, U12, or a combination of both types of splice site consensus sequences, with a shift towards splicing using preferentially U2 sites in TALS patients' cells compared to controls.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5′ and 3′ consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Enhancer Domains in Gastrointestinal Stromal Tumor Regulate KIT Expression and Are Targetable by BET Bromodomain Inhibition

Gastrointestinal stromal tumor (GIST) is a mesenchymal neoplasm characterized by activating mutations in the related receptor tyrosine kinases KIT and PDGFRA. GIST relies on expression of these unamplified receptor tyrosine kinase (RTK) genes through a large enhancer domain, resulting in high expression levels of the oncogene required for tumor growth. Although kinase inhibition is an effective therapy for many patients with GIST, disease progression from kinase-resistant mutations is common and no other effective classes of systemic therapy exist. In this study, we identify regulatory regions of the KIT enhancer essential for KIT gene expression and GIST cell viability. Given the dependence of GIST upon enhancer-driven expression of RTKs, we hypothesized that the enhancer domains could be therapeutically targeted by a BET bromodomain inhibitor (BBI). Treatment of GIST cells with BBIs led to cell-cycle arrest, apoptosis, and cell death, with unique sensitivity in GIST cells arising from attenuation of the KIT enhancer domain and reduced KIT gene expression. BBI treatment in KIT-dependent GIST cells produced genome-wide changes in the H3K27ac enhancer landscape and gene expression program, which was also seen with direct KIT inhibition using a tyrosine kinase inhibitor (TKI). Combination treatment with BBI and TKI led to superior cytotoxic effects in vitro and in vivo, with BBI preventing tumor growth in TKI-resistant xenografts. Resistance to select BBI in GIST was attributable to drug efflux pumps. These results define a therapeutic vulnerability and clinical strategy for targeting oncogenic kinase dependency in GIST. SIGNIFICANCE: Expression and activity of mutant KIT is essential for driving the majority of GIST neoplasms, which can be therapeutically targeted using BET bromodomain inhibitors.

Familial breast cancer and DNA repair genes: Insights into known and novel susceptibility genes from the GENESIS study, and implications for multigene panel testing

Pathogenic variants in BRCA1 and BRCA2 only explain the underlying genetic cause of about 10% of hereditary breast and ovarian cancer families. Because of cost‐effectiveness, multigene panel testing is often performed even if the clinical utility of testing most of the genes remains questionable. The purpose of our study was to assess the contribution of rare, deleterious‐predicted variants in DNA repair genes in familial breast cancer (BC) in a well‐characterized and homogeneous population. We analyzed 113 DNA repair genes selected from either an exome sequencing or a candidate gene approach in the GENESIS study, which includes familial BC cases with no BRCA1 or BRCA2 mutation and having a sister with BC (N = 1,207), and general population controls (N = 1,199). Sequencing data were filtered for rare loss‐of‐function variants (LoF) and likely deleterious missense variants (MV). We confirmed associations between LoF and MV in PALB2, ATM and CHEK2 and BC occurrence. We also identified for the first time associations between FANCI, MAST1, POLH and RTEL1 and BC susceptibility. Unlike other associated genes, carriers of an ATM LoF had a significantly higher risk of developing BC than carriers of an ATM MV (OR_LoF = 17.4 vs. OR_MV = 1.6; p_Het = 0.002). Hence, our approach allowed us to specify BC relative risks associated with deleterious‐predicted variants in PALB2, ATM and CHEK2 and to add MAST1, POLH, RTEL1 and FANCI to the list of DNA repair genes possibly involved in BC susceptibility. We also highlight that different types of variants within the same gene can lead to different risk estimates.

What's new?

Pathogenic variants in BRCA1 and BRCA2 only explain the genetic cause of about 10% of hereditary breast and ovarian cancer families, and the clinical usefulness of testing other genes following the recent introduction of cost‐effective multigene panel sequencing in diagnostics laboratories remains questionable. This large case‐control study describes genetic variation in 113 DNA repair genes and specifies breast cancer relative risks associated with rare deleterious‐predicted variants in PALB2, ATM, and CHEK2. Importantly, different types of variants within the same gene can lead to different risk estimates. The results may help improve risk prediction models and define gene‐specific consensus management guidelines.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

Prediction of Breast and Prostate Cancer Risks in Male BRCA1 and BRCA2 Mutation Carriers Using Polygenic Risk Scores

Purpose BRCA1/2 mutations increase the risk of breast and prostate cancer in men. Common genetic variants modify cancer risks for female carriers of BRCA1/2 mutations. We investigated-for the first time to our knowledge-associations of common genetic variants with breast and prostate cancer risks for male carriers of BRCA1/ 2 mutations and implications for cancer risk prediction. Materials and Methods We genotyped 1,802 male carriers of BRCA1/2 mutations from the Consortium of Investigators of Modifiers of BRCA1/2 by using the custom Illumina OncoArray. We investigated the combined effects of established breast and prostate cancer susceptibility variants on cancer risks for male carriers of BRCA1/2 mutations by constructing weighted polygenic risk scores (PRSs) using published effect estimates as weights. Results In male carriers of BRCA1/2 mutations, PRS that was based on 88 female breast cancer susceptibility variants was associated with breast cancer risk (odds ratio per standard deviation of PRS, 1.36; 95% CI, 1.19 to 1.56; P = 8.6 × 10-6). Similarly, PRS that was based on 103 prostate cancer susceptibility variants was associated with prostate cancer risk (odds ratio per SD of PRS, 1.56; 95% CI, 1.35 to 1.81; P = 3.2 × 10-9). Large differences in absolute cancer risks were observed at the extremes of the PRS distribution. For example, prostate cancer risk by age 80 years at the 5th and 95th percentiles of the PRS varies from 7% to 26% for carriers of BRCA1 mutations and from 19% to 61% for carriers of BRCA2 mutations, respectively. Conclusion PRSs may provide informative cancer risk stratification for male carriers of BRCA1/2 mutations that might enable these men and their physicians to make informed decisions on the type and timing of breast and prostate cancer risk management.

Evaluation of Polygenic Risk Scores for Breast and Ovarian Cancer Risk Prediction in BRCA1 and BRCA2 Mutation Carriers

Background

Genome-wide association studies (GWAS) have identified 94 common single-nucleotide polymorphisms (SNPs) associated with breast cancer (BC) risk and 18 associated with ovarian cancer (OC) risk. Several of these are also associated with risk of BC or OC for women who carry a pathogenic mutation in the high-risk BC and OC genes BRCA1 or BRCA2. The combined effects of these variants on BC or OC risk for BRCA1 and BRCA2 mutation carriers have not yet been assessed while their clinical management could benefit from improved personalized risk estimates.

Methods

We constructed polygenic risk scores (PRS) using BC and OC susceptibility SNPs identified through population-based GWAS: for BC (overall, estrogen receptor [ER]-positive, and ER-negative) and for OC. Using data from 15 252 female BRCA1 and 8211 BRCA2 carriers, the association of each PRS with BC or OC risk was evaluated using a weighted cohort approach, with time to diagnosis as the outcome and estimation of the hazard ratios (HRs) per standard deviation increase in the PRS.

Results

The PRS for ER-negative BC displayed the strongest association with BC risk in BRCA1 carriers (HR = 1.27, 95% confidence interval [CI] = 1.23 to 1.31, P =  8.2×10 -53 ). In BRCA2 carriers, the strongest association with BC risk was seen for the overall BC PRS (HR = 1.22, 95% CI = 1.17 to 1.28, P =  7.2×10 -20 ). The OC PRS was strongly associated with OC risk for both BRCA1 and BRCA2 carriers. These translate to differences in absolute risks (more than 10% in each case) between the top and bottom deciles of the PRS distribution; for example, the OC risk was 6% by age 80 years for BRCA2 carriers at the 10th percentile of the OC PRS compared with 19% risk for those at the 90th percentile of PRS.

Conclusions

BC and OC PRS are predictive of cancer risk in BRCA1 and BRCA2 carriers. Incorporation of the PRS into risk prediction models has promise to better inform decisions on cancer risk management.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

Association of breast cancer risk in BRCA1 and BRCA2 mutation carriers with genetic variants showing differential allelic expression: identification of a modifier of breast cancer risk at locus 11q22.3

PURPOSE

Cis-acting regulatory SNPs resulting in differential allelic expression (DAE) may, in part, explain the underlying phenotypic variation associated with many complex diseases. To investigate whether common variants associated with DAE were involved in breast cancer susceptibility among BRCA1 and BRCA2 mutation carriers, a list of 175 genes was developed based of their involvement in cancer-related pathways.

METHODS

Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of ~320 SNPs located in the vicinity of these genes with breast and ovarian cancer risks in 15,252 BRCA1 and 8211 BRCA2 mutation carriers ascertained from 54 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2.

RESULTS

We identified a region on 11q22.3 that is significantly associated with breast cancer risk in BRCA1 mutation carriers (most significant SNP rs228595 p = 7 × 10-6). This association was absent in BRCA2 carriers (p = 0.57). The 11q22.3 region notably encompasses genes such as ACAT1, NPAT, and ATM. Expression quantitative trait loci associations were observed in both normal breast and tumors across this region, namely for ACAT1, ATM, and other genes. In silico analysis revealed some overlap between top risk-associated SNPs and relevant biological features in mammary cell data, which suggests potential functional significance.

CONCLUSION

We identified 11q22.3 as a new modifier locus in BRCA1 carriers. Replication in larger studies using estrogen receptor (ER)-negative or triple-negative (i.e., ER-, progesterone receptor-, and HER2-negative) cases could therefore be helpful to confirm the association of this locus with breast cancer risk.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Fine-Scale Mapping at 9p22.2 Identifies Candidate Causal Variants That Modify Ovarian Cancer Risk in BRCA1 and BRCA2 Mutation Carriers

Population-based genome wide association studies have identified a locus at 9p22.2 associated with ovarian cancer risk, which also modifies ovarian cancer risk in BRCA1 and BRCA2 mutation carriers. We conducted fine-scale mapping at 9p22.2 to identify potential causal variants in BRCA1 and BRCA2 mutation carriers. Genotype data were available for 15,252 (2,462 ovarian cancer cases) BRCA1 and 8,211 (631 ovarian cancer cases) BRCA2 mutation carriers. Following genotype imputation, ovarian cancer associations were assessed for 4,873 and 5,020 SNPs in BRCA1 and BRCA 2 mutation carriers respectively, within a retrospective cohort analytical framework. In BRCA1 mutation carriers one set of eight correlated candidate causal variants for ovarian cancer risk modification was identified (top SNP rs10124837, HR: 0.73, 95%CI: 0.68 to 0.79, p-value 2× 10-16). These variants were located up to 20 kb upstream of BNC2. In BRCA2 mutation carriers one region, up to 45 kb upstream of BNC2, and containing 100 correlated SNPs was identified as candidate causal (top SNP rs62543585, HR: 0.69, 95%CI: 0.59 to 0.80, p-value 1.0 × 10-6). The candidate causal in BRCA1 mutation carriers did not include the strongest associated variant at this locus in the general population. In sum, we identified a set of candidate causal variants in a region that encompasses the BNC2 transcription start site. The ovarian cancer association at 9p22.2 may be mediated by different variants in BRCA1 mutation carriers and in the general population. Thus, potentially different mechanisms may underlie ovarian cancer risk for mutation carriers and the general population.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND

Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk.

METHOD

We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation.

RESULTS

Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05.

CONCLUSION

This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.

Refining the phenotypical and mutational spectrum of Taybi-Linder syndrome

Taybi-Linder syndrome (TALS, OMIM 210710) is a rare autosomal recessive disorder belonging to the group of microcephalic osteodysplastic primordial dwarfisms (MOPD). This syndrome is characterized by short stature, skeletal anomalies, severe microcephaly with brain malformations and facial dysmorphism, and is caused by mutations in RNU4ATAC. RNU4ATAC is transcribed into a non-coding small nuclear RNA which is a critical component of the minor spliceosome. We report here four foetuses and four unrelated patients with RNU4ATAC mutations. We provide antenatal descriptions of this rare syndrome including unusual features found in two twin foetuses with compound heterozygosity for two rare mutations who presented with mild intrauterine growth retardation and atypical dysmorphic facial features. We also carried out a literature review of the patients described up to now with RNU4ATAC mutations, affected either with TALS or Roifman syndrome, a recently described allelic disorder.

No clinical utility of KRAS variant rs61764370 for ovarian or breast cancer

OBJECTIVE

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.

METHODS

Centralized genotyping and analysis were performed for 140,012 women enrolled in the Ovarian Cancer Association Consortium (15,357 ovarian cancer patients; 30,816 controls), the Breast Cancer Association Consortium (33,530 breast cancer patients; 37,640 controls), and the Consortium of Modifiers of BRCA1 and BRCA2 (14,765 BRCA1 and 7904 BRCA2 mutation carriers).

RESULTS

We found no association with risk of ovarian cancer (OR=0.99, 95% CI 0.94-1.04, p=0.74) or breast cancer (OR=0.98, 95% CI 0.94-1.01, p=0.19) and results were consistent among mutation carriers (BRCA1, ovarian cancer HR=1.09, 95% CI 0.97-1.23, p=0.14, breast cancer HR=1.04, 95% CI 0.97-1.12, p=0.27; BRCA2, ovarian cancer HR=0.89, 95% CI 0.71-1.13, p=0.34, breast cancer HR=1.06, 95% CI 0.94-1.19, p=0.35). Null results were also obtained for associations with overall survival following ovarian cancer (HR=0.94, 95% CI 0.83-1.07, p=0.38), breast cancer (HR=0.96, 95% CI 0.87-1.06, p=0.38), and all other previously-reported associations.

CONCLUSIONS

rs61764370 is not associated with risk of ovarian or breast cancer nor with clinical outcome for patients with these cancers. Therefore, genotyping this variant has no clinical utility related to the prediction or management of these cancers.

Identification of four novel susceptibility loci for oestrogen receptor negative breast cancer

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.

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.

Abstract P2-10-01: The BRCA-1 polymorphism (major homozygous rs5820483) is associated with higher expression of phosphorylated -IGF-1 receptor

Background: Data from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) indicates that common variation in the BRCA1 locus modifies the penetrance of mutations. The rs5820483 variation is one such polymorphism in which patients homozygous for the major allele have a greater risk of breast cancer. The BRCA1 gene contains 24 exons that generate a full-length (FL) protein. Major alternatively spliced isoforms include skipping of exon 11 (Δ11) and a truncated form (IRIS), each of which appear to have distinct roles in cellular processes. rs5820483 is the only variant that maps to intron 10, flanking exon 11. Experimental models suggest that deregulation of the FL/D11 isoform ratio can have profound functional consequences. We have shown that cells from patients with the major rs5820483 allele express more FL relative to Δ11 compared to cells from patients with the minor allele. Breast cancer exhibit increased IGF-1 and IGF-1 activity in patients with BRCA1 mutations relative to breast cancer in women without BRCA1 mutations. IGF-1 produced in the mammary stroma mediates estrogen dependent proliferation. Here we test whether the rs5820483 allele affects IGF-1 activity in breast tissue of BRCA1 mutation carriers.

Design: We assessed the zygosity of the rs5820483 alleles in 28 BRCA1 mutation carriers. Phosphorylated IGF-1 receptor (p-IGF-1R rabbit polyclonal Ab39398, Abcam, MA) was measured by immunofluorescence and Ki-67 was measured by immunohistochemistry in sections from histologically normal breast tissue from these cases were compared to heterozygous controls. Image analysis was used to assess the intensity of the p-IGF-1R immunofluorescence at the epithelial cell membrane with appropriate controls. Ki67 mmunohistochemistry was also performed on these specimens.

Results:The rs5820483 allele was heterozygous in 10 specimens, homozygous for the major allele in 11 and homozygous for the minor allele in 7. The intensity of p-IGF-1R was higher in major homozygous cases (140±54 SD) than in either the minor homozygous (99±24) or heterozygous cases (95±44). The frequency of Ki-67+ cells was higher in the major homozygous case (3.1±2.9 SD) than in either the minor homozygous (1.5±1.9 SD) or heterozygous cases (2.2±2.5 SD). However, neither p-IGF-IR immunoreactivity nor frequency of Ki-67+ cells was statistically different between groups.

Conclusion: Breast cancer risk in BRCA1 mutation carriers is modified by common genetic variation at the corresponding locus. We have identified a variation at rs5820483 that affects the isoform ratio. Our preliminary analysis suggests that IGF-1 activity increases in those women homozygous for the major allele, concomitant with increased Ki-67. Corroboration of this analysis in a larger series is ongoing. A statistically significant difference might have fundamental implications for cancer prevention in those carriers.

Citation Format: Singh B, Bochaca II, Ruiz de Garibay G, Damiola F, Mazoyer S, Antoniou A, Chenevix-Trench G, Pujana MA, Barcellos-Hoff MH, Kleinberg D. The BRCA-1 polymorphism (major homozygous rs5820483) is associated with higher expression of phosphorylated -IGF-1 receptor. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P2-10-01.

BRCA2 Polymorphic Stop Codon K3326X and the Risk of Breast, Prostate, and Ovarian Cancers

BACKGROUND

The K3326X variant in BRCA2 (BRCA2*c.9976A>T; p.Lys3326*; rs11571833) has been found to be associated with small increased risks of breast cancer. However, it is not clear to what extent linkage disequilibrium with fully pathogenic mutations might account for this association. There is scant information about the effect of K3326X in other hormone-related cancers.

METHODS

Using weighted logistic regression, we analyzed data from the large iCOGS study including 76 637 cancer case patients and 83 796 control patients to estimate odds ratios (ORw) and 95% confidence intervals (CIs) for K3326X variant carriers in relation to breast, ovarian, and prostate cancer risks, with weights defined as probability of not having a pathogenic BRCA2 variant. Using Cox proportional hazards modeling, we also examined the associations of K3326X with breast and ovarian cancer risks among 7183 BRCA1 variant carriers. All statistical tests were two-sided.

RESULTS

The K3326X variant was associated with breast (ORw = 1.28, 95% CI = 1.17 to 1.40, P = 5.9x10(-) (6)) and invasive ovarian cancer (ORw = 1.26, 95% CI = 1.10 to 1.43, P = 3.8x10(-3)). These associations were stronger for serous ovarian cancer and for estrogen receptor-negative breast cancer (ORw = 1.46, 95% CI = 1.2 to 1.70, P = 3.4x10(-5) and ORw = 1.50, 95% CI = 1.28 to 1.76, P = 4.1x10(-5), respectively). For BRCA1 mutation carriers, there was a statistically significant inverse association of the K3326X variant with risk of ovarian cancer (HR = 0.43, 95% CI = 0.22 to 0.84, P = .013) but no association with breast cancer. No association with prostate cancer was observed.

CONCLUSIONS

Our study provides evidence that the K3326X variant is associated with risk of developing breast and ovarian cancers independent of other pathogenic variants in BRCA2. Further studies are needed to determine the biological mechanism of action responsible for these associations.

Mutation screening of MIR146A/B and BRCA1/2 3'-UTRs in the GENESIS study

Although a wide number of breast cancer susceptibility alleles associated with various levels of risk have been identified to date, about 50% of the heritability is still missing. Although the major BRCA1 and BRCA2 genes are being extensively screened for truncating and missense variants in breast and/or ovarian cancer families, potential regulatory variants affecting their expression remain largely unexplored. In an attempt to identify such variants, we focused our attention on gene regulation mediated by microRNAs (miRs). We screened two genes, MIR146A and MIR146B, producing miR-146a and miR-146b-5p, respectively, that regulate BRCA1, and the 3'- untranslated regions (3'-UTRs) of BRCA1 and BRCA2 in the GENESIS French national case/control study (BRCA1- and BRCA2-negative breast cancer cases with at least one sister with breast cancer and matched controls). We identified one rare variant in MIR146A, four in MIR146B, five in BRCA1 3'-UTR and one in BRCA2 3'-UTR in 716 index cases and 619 controls. Among these 11 rare variants, 7 were identified each in 1 index case. None of the three relevant MIR146A/MIR146B variants affected the pre-miR sequences. The potential causality of the four relevant BRCA1/BRCA2 3'-UTRs variants was evaluated with luciferase reporter assays and co-segregation studies, as well as with bioinformatics analyses to predict miRs-binding sites, RNA secondary structures and RNA accessibility. This is the first study to report the screening of miR genes and of BRCA2 3'-UTR in a large series of familial breast cancer cases. None of the variant identified in this study gave convincing evidence of potential pathogenicity.

GENESIS: a French national resource to study the missing heritability of breast cancer

Background

Less than 20 % of familial breast cancer patients who undergo genetic testing for BRCA1 and BRCA2 carry a pathogenic mutation in one of these two genes. The GENESIS (GENE SISter) study was designed to identify new breast cancer susceptibility genes in women attending cancer genetics clinics and with no BRCA1/2 mutation.

Methods

The study involved the French national network of family cancer clinics. It was based on enrichment in genetic factors of the recruited population through case selection relying on familial criteria, but also on the consideration of environmental factors and endophenotypes like mammary density or tumor characteristics to assess potential genetic heterogeneity. One of the initial aims of GENESIS was to recruit affected sibpairs. Siblings were eligible when index cases and at least one affected sister were diagnosed with infiltrating mammary or ductal adenocarcinoma, with no BRCA1/2 mutation. In addition, unrelated controls and unaffected sisters were recruited. The enrolment of patients, their relatives and their controls, the collection of the clinical, epidemiological, familial and biological data were centralized by a coordinating center.

Results

Inclusion of participants started in February 2007 and ended in December 2013. A total of 1721 index cases, 826 affected sisters, 599 unaffected sisters and 1419 controls were included. 98 % of participants completed the epidemiological questionnaire, 97 % provided a blood sample, and 76 % were able to provide mammograms. Index cases were on average 59 years old at inclusion, were born in 1950, and were 49.7 years of age at breast cancer diagnosis. The mean age at diagnosis of affected sisters was slightly higher (51.4 years). The representativeness of the control group was verified.

Conclusions

The size of the study, the availability of biological specimens and the clinical data collection together with the detailed and complete epidemiological questionnaire make this a unique national resource for investigation of the missing heritability of breast cancer, by taking into account environmental and life style factors and stratifying data on endophenotypes to decrease genetic heterogeneity.

Targeted Sequencing of the Mitochondrial Genome of Women at High Risk of Breast Cancer without Detectable Mutations in BRCA1/2

Breast Cancer is a complex multifactorial disease for which high-penetrance mutations have been identified. Approaches used to date have identified genomic features explaining about 50% of breast cancer heritability. A number of low- to medium penetrance alleles (per-allele odds ratio < 1.5 and 4.0, respectively) have been identified, suggesting that the remaining heritability is likely to be explained by the cumulative effect of such alleles and/or by rare high-penetrance alleles. Relatively few studies have specifically explored the mitochondrial genome for variants potentially implicated in breast cancer risk. For these reasons, we propose an exploration of the variability of the mitochondrial genome in individuals diagnosed with breast cancer, having a positive breast cancer family history but testing negative for BRCA1/2 pathogenic mutations. We sequenced the mitochondrial genome of 436 index breast cancer cases from the GENESIS study. As expected, no pathogenic genomic pattern common to the 436 women included in our study was observed. The mitochondrial genes MT-ATP6 and MT-CYB were observed to carry the highest number of variants in the study. The proteins encoded by these genes are involved in the structure of the mitochondrial respiration chain, and variants in these genes may impact reactive oxygen species production contributing to carcinogenesis. More functional and epidemiological studies are needed to further investigate to what extent variants identified may influence familial breast cancer risk.

Occurrence of a non deleterious gene conversion event in the BRCA1 gene

The duplication in the primate lineage of a portion of the breast and ovarian cancer susceptibility gene BRCA1 has created a BRCA1 pseudogene 45 kb away. Non-allelic homologous recombination (NAHR) between BRCA1 and BRCA1P1 has generated recurrent deleterious germ-line 37-kb deletions encompassing the first two exons of BRCA1, accounting for several breast and ovarian cancer families in various populations. In principle, NAHR intermediates resolution could also lead through a non-crossover configuration to interlocus gene conversion (IGC), but none had been described as yet. Here, we report for the first time an IGC event identified in a breast and ovarian cancer family involving exactly the same segment as that involved in the 37-kb deletions. Close examination of the consequences of this IGC event showed that it does not impact BRCA1 expression. Detailed analysis of the regions of homology between BRCA1 and its pseudogene revealed the specificity of the segment where recombination systematically occurs.

FANCM c.5791C>T nonsense mutation (rs144567652) induces exon skipping, affects DNA repair activity and is a familial breast cancer risk factor

Numerous genetic factors that influence breast cancer risk are known. However, approximately two-thirds of the overall familial risk remain unexplained. To determine whether some of the missing heritability is due to rare variants conferring high to moderate risk, we tested for an association between the c.5791C>T nonsense mutation (p.Arg1931*; rs144567652) in exon 22 of FANCM gene and breast cancer. An analysis of genotyping data from 8635 familial breast cancer cases and 6625 controls from different countries yielded an association between the c.5791C>T mutation and breast cancer risk [odds ratio (OR) = 3.93 (95% confidence interval (CI) = 1.28-12.11; P = 0.017)]. Moreover, we performed two meta-analyses of studies from countries with carriers in both cases and controls and of all available data. These analyses showed breast cancer associations with OR = 3.67 (95% CI = 1.04-12.87; P = 0.043) and OR = 3.33 (95% CI = 1.09-13.62; P = 0.032), respectively. Based on information theory-based prediction, we established that the mutation caused an out-of-frame deletion of exon 22, due to the creation of a binding site for the pre-mRNA processing protein hnRNP A1. Furthermore, genetic complementation analyses showed that the mutation influenced the DNA repair activity of the FANCM protein. In summary, we provide evidence for the first time showing that the common p.Arg1931* loss-of-function variant in FANCM is a risk factor for familial breast cancer.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

INTRODUCTION

Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

METHODS

We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

RESULTS

We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

CONCLUSIONS

This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Association of type and location of BRCA1 and BRCA2 mutations with risk of breast and ovarian cancer

IMPORTANCE

Limited information about the relationship between specific mutations in BRCA1 or BRCA2 (BRCA1/2) and cancer risk exists.

OBJECTIVE

To identify mutation-specific cancer risks for carriers of BRCA1/2.

DESIGN, SETTING, AND PARTICIPANTS

Observational study of women who were ascertained between 1937 and 2011 (median, 1999) and found to carry disease-associated BRCA1 or BRCA2 mutations. The international sample comprised 19,581 carriers of BRCA1 mutations and 11,900 carriers of BRCA2 mutations from 55 centers in 33 countries on 6 continents. We estimated hazard ratios for breast and ovarian cancer based on mutation type, function, and nucleotide position. We also estimated RHR, the ratio of breast vs ovarian cancer hazard ratios. A value of RHR greater than 1 indicated elevated breast cancer risk; a value of RHR less than 1 indicated elevated ovarian cancer risk.

EXPOSURES

Mutations of BRCA1 or BRCA2.

MAIN OUTCOMES AND MEASURES

Breast and ovarian cancer risks.

RESULTS

Among BRCA1 mutation carriers, 9052 women (46%) were diagnosed with breast cancer, 2317 (12%) with ovarian cancer, 1041 (5%) with breast and ovarian cancer, and 7171 (37%) without cancer. Among BRCA2 mutation carriers, 6180 women (52%) were diagnosed with breast cancer, 682 (6%) with ovarian cancer, 272 (2%) with breast and ovarian cancer, and 4766 (40%) without cancer. In BRCA1, we identified 3 breast cancer cluster regions (BCCRs) located at c.179 to c.505 (BCCR1; RHR = 1.46; 95% CI, 1.22-1.74; P = 2 × 10(-6)), c.4328 to c.4945 (BCCR2; RHR = 1.34; 95% CI, 1.01-1.78; P = .04), and c. 5261 to c.5563 (BCCR2', RHR = 1.38; 95% CI, 1.22-1.55; P = 6 × 10(-9)). We also identified an ovarian cancer cluster region (OCCR) from c.1380 to c.4062 (approximately exon 11) with RHR = 0.62 (95% CI, 0.56-0.70; P = 9 × 10(-17)). In BRCA2, we observed multiple BCCRs spanning c.1 to c.596 (BCCR1; RHR = 1.71; 95% CI, 1.06-2.78; P = .03), c.772 to c.1806 (BCCR1'; RHR = 1.63; 95% CI, 1.10-2.40; P = .01), and c.7394 to c.8904 (BCCR2; RHR = 2.31; 95% CI, 1.69-3.16; P = .00002). We also identified 3 OCCRs: the first (OCCR1) spanned c.3249 to c.5681 that was adjacent to c.5946delT (6174delT; RHR = 0.51; 95% CI, 0.44-0.60; P = 6 × 10(-17)). The second OCCR spanned c.6645 to c.7471 (OCCR2; RHR = 0.57; 95% CI, 0.41-0.80; P = .001). Mutations conferring nonsense-mediated decay were associated with differential breast or ovarian cancer risks and an earlier age of breast cancer diagnosis for both BRCA1 and BRCA2 mutation carriers.

CONCLUSIONS AND RELEVANCE

Breast and ovarian cancer risks varied by type and location of BRCA1/2 mutations. With appropriate validation, these data may have implications for risk assessment and cancer prevention decision making for carriers of BRCA1 and BRCA2 mutations.

Assessing associations between the AURKA-HMMR-TPX2-TUBG1 functional module and breast cancer risk in BRCA1/2 mutation carriers

While interplay between BRCA1 and AURKA-RHAMM-TPX2-TUBG1 regulates mammary epithelial polarization, common genetic variation in HMMR (gene product RHAMM) may be associated with risk of breast cancer in BRCA1 mutation carriers. Following on these observations, we further assessed the link between the AURKA-HMMR-TPX2-TUBG1 functional module and risk of breast cancer in BRCA1 or BRCA2 mutation carriers. Forty-one single nucleotide polymorphisms (SNPs) were genotyped in 15,252 BRCA1 and 8,211 BRCA2 mutation carriers and subsequently analyzed using a retrospective likelihood approach. The association of HMMR rs299290 with breast cancer risk in BRCA1 mutation carriers was confirmed: per-allele hazard ratio (HR) = 1.10, 95% confidence interval (CI) 1.04-1.15, p = 1.9 x 10(-4) (false discovery rate (FDR)-adjusted p = 0.043). Variation in CSTF1, located next to AURKA, was also found to be associated with breast cancer risk in BRCA2 mutation carriers: rs2426618 per-allele HR = 1.10, 95% CI 1.03-1.16, p = 0.005 (FDR-adjusted p = 0.045). Assessment of pairwise interactions provided suggestions (FDR-adjusted pinteraction values > 0.05) for deviations from the multiplicative model for rs299290 and CSTF1 rs6064391, and rs299290 and TUBG1 rs11649877 in both BRCA1 and BRCA2 mutation carriers. Following these suggestions, the expression of HMMR and AURKA or TUBG1 in sporadic breast tumors was found to potentially interact, influencing patients' survival. Together, the results of this study support the hypothesis of a causative link between altered function of AURKA-HMMR-TPX2-TUBG1 and breast carcinogenesis in BRCA1/2 mutation carriers.

Identification of six new susceptibility loci for invasive epithelial ovarian cancer

Genome-wide association studies (GWAS) have identified 12 epithelial ovarian cancer (EOC) susceptibility alleles. The pattern of association at these loci is consistent in BRCA1 and BRCA2 mutation carriers who are at high risk of EOC. After imputation to 1000 Genomes Project data, we assessed associations of 11 million genetic variants with EOC risk from 15,437 cases unselected for family history and 30,845 controls and from 15,252 BRCA1 mutation carriers and 8,211 BRCA2 mutation carriers (3,096 with ovarian cancer), and we combined the results in a meta-analysis. This new study design yielded increased statistical power, leading to the discovery of six new EOC susceptibility loci. Variants at 1p36 (nearest gene, WNT4), 4q26 (SYNPO2), 9q34.2 (ABO) and 17q11.2 (ATAD5) were associated with EOC risk, and at 1p34.3 (RSPO1) and 6p22.1 (GPX6) variants were specifically associated with the serous EOC subtype, all with P < 5 × 10(-8). Incorporating these variants into risk assessment tools will improve clinical risk predictions for BRCA1 and BRCA2 mutation carriers.

Candidate genetic modifiers for breast and ovarian cancer risk in BRCA1 and BRCA2 mutation carriers

BACKGROUND

BRCA1 and BRCA2 mutation carriers are at substantially increased risk for developing breast and ovarian cancer. The incomplete penetrance coupled with the variable age at diagnosis in carriers of the same mutation suggests the existence of genetic and nongenetic modifying factors. In this study, we evaluated the putative role of variants in many candidate modifier genes.

METHODS

Genotyping data from 15,252 BRCA1 and 8,211 BRCA2 mutation carriers, for known variants (n = 3,248) located within or around 445 candidate genes, were available through the iCOGS custom-designed array. Breast and ovarian cancer association analysis was performed within a retrospective cohort approach.

RESULTS

The observed P values of association ranged between 0.005 and 1.000. None of the variants was significantly associated with breast or ovarian cancer risk in either BRCA1 or BRCA2 mutation carriers, after multiple testing adjustments.

CONCLUSION

There is little evidence that any of the evaluated candidate variants act as modifiers of breast and/or ovarian cancer risk in BRCA1 or BRCA2 mutation carriers.

IMPACT

Genome-wide association studies have been more successful at identifying genetic modifiers of BRCA1/2 penetrance than candidate gene studies.

Associations of common breast cancer susceptibility alleles with risk of breast cancer subtypes in BRCA1 and BRCA2 mutation carriers

INTRODUCTION

More than 70 common alleles are known to be involved in breast cancer (BC) susceptibility, and several exhibit significant heterogeneity in their associations with different BC subtypes. Although there are differences in the association patterns between BRCA1 and BRCA2 mutation carriers and the general population for several loci, no study has comprehensively evaluated the associations of all known BC susceptibility alleles with risk of BC subtypes in BRCA1 and BRCA2 carriers.

METHODS

We used data from 15,252 BRCA1 and 8,211 BRCA2 carriers to analyze the associations between approximately 200,000 genetic variants on the iCOGS array and risk of BC subtypes defined by estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and triple-negative- (TN) status; morphologic subtypes; histological grade; and nodal involvement.

RESULTS

The estimated BC hazard ratios (HRs) for the 74 known BC alleles in BRCA1 carriers exhibited moderate correlations with the corresponding odds ratios from the general population. However, their associations with ER-positive BC in BRCA1 carriers were more consistent with the ER-positive associations in the general population (intraclass correlation (ICC) = 0.61, 95% confidence interval (CI): 0.45 to 0.74), and the same was true when considering ER-negative associations in both groups (ICC = 0.59, 95% CI: 0.42 to 0.72). Similarly, there was strong correlation between the ER-positive associations for BRCA1 and BRCA2 carriers (ICC = 0.67, 95% CI: 0.52 to 0.78), whereas ER-positive associations in any one of the groups were generally inconsistent with ER-negative associations in any of the others. After stratifying by ER status in mutation carriers, additional significant associations were observed. Several previously unreported variants exhibited associations at P <10(-6) in the analyses by PR status, HER2 status, TN phenotype, morphologic subtypes, histological grade and nodal involvement.

CONCLUSIONS

Differences in associations of common BC susceptibility alleles between BRCA1 and BRCA2 carriers and the general population are explained to a large extent by differences in the prevalence of ER-positive and ER-negative tumors. Estimates of the risks associated with these variants based on population-based studies are likely to be applicable to mutation carriers after taking ER status into account, which has implications for risk prediction.

Estimation of the RNU2 macrosatellite mutation rate by BRCA1 mutation tracing

Large tandem repeat sequences have been poorly investigated as severe technical limitations and their frequent absence from the genome reference hinder their analysis. Extensive allelotyping of this class of variation has not been possible until now and their mutational dynamics are still poorly known. In order to estimate the mutation rate of a macrosatellite, we analysed in detail the RNU2 locus, which displays at least 50 different alleles containing 5-82 copies of a 6.1 kb repeat unit. Mining data from the 1000 Genomes Project allowed us to precisely estimate copy numbers of the RNU2 repeat unit using read depth of coverage. This further revealed significantly different mean values in various recent modern human populations, favoring a scenario of fast evolution of this locus. Its proximity to a disease gene with numerous founder mutations, BRCA1, within the same linkage disequilibrium block, offered the unique opportunity to trace RNU2 arrays over a large timescale. Analysis of the transmission of RNU2 arrays associated with one ‘private’ mutation in an extended kindred and four founder mutations in multiple kindreds gave an estimation by maximum likelihood of 5 × 10⁻³ mutations per generation, which is close to that of microsatellites.

DNA glycosylases involved in base excision repair may be associated with cancer risk in BRCA1 and BRCA2 mutation carriers

Single Nucleotide Polymorphisms (SNPs) in genes involved in the DNA Base Excision Repair (BER) pathway could be associated with cancer risk in carriers of mutations in the high-penetrance susceptibility genes BRCA1 and BRCA2, given the relation of synthetic lethality that exists between one of the components of the BER pathway, PARP1 (poly ADP ribose polymerase), and both BRCA1 and BRCA2. In the present study, we have performed a comprehensive analysis of 18 genes involved in BER using a tagging SNP approach in a large series of BRCA1 and BRCA2 mutation carriers. 144 SNPs were analyzed in a two stage study involving 23,463 carriers from the CIMBA consortium (the Consortium of Investigators of Modifiers of BRCA1 and BRCA2). Eleven SNPs showed evidence of association with breast and/or ovarian cancer at p<0.05 in the combined analysis. Four of the five genes for which strongest evidence of association was observed were DNA glycosylases. The strongest evidence was for rs1466785 in the NEIL2 (endonuclease VIII-like 2) gene (HR: 1.09, 95% CI (1.03-1.16), p = 2.7 × 10(-3)) for association with breast cancer risk in BRCA2 mutation carriers, and rs2304277 in the OGG1 (8-guanine DNA glycosylase) gene, with ovarian cancer risk in BRCA1 mutation carriers (HR: 1.12 95%CI: 1.03-1.21, p = 4.8 × 10(-3)). DNA glycosylases involved in the first steps of the BER pathway may be associated with cancer risk in BRCA1/2 mutation carriers and should be more comprehensively studied.

Direct visualization of the highly polymorphic RNU2 locus in proximity to the BRCA1 gene

Although the breast cancer susceptibility gene BRCA1 is one of the most extensively characterized genetic loci, much less is known about its upstream variable number tandem repeat element, the RNU2 locus. RNU2 encodes the U2 small nuclear RNA, an essential splicing element, but this locus is missing from the human genome assembly due to the inherent difficulty in the assembly of repetitive sequences. To fill the gap between RNU2 and BRCA1, we have reconstructed the physical map of this region by re-examining genomic clone sequences of public databases, which allowed us to precisely localize the RNU2 array 124 kb telomeric to BRCA1. We measured by performing FISH analyses on combed DNA for the first time the exact number of repeats carried by each of the two alleles in 41 individuals and found a range of 6-82 copies and a level of heterozygosity of 98%. The precise localisation of the RNU2 locus in the genome reference assembly and the implementation of a new technical tool to study it will make the detailed exploration of this locus possible. This recently neglected macrosatellite could be valuable for evaluating the potential role of structural variations in disease due to its location next to a major cancer susceptibility gene.

BRCA1-Dependent Translational Regulation in Breast Cancer Cells

BRCA1 (Breast Cancer 1) has been implicated in a number of cellular processes, including transcription regulation, DNA damage repair and protein ubiquitination. We previously demonstrated that BRCA1 interacts with PABP1 (Poly(A)-Binding Protein 1) and that BRCA1 modulates protein synthesis through this interaction. To identify the mRNAs that are translationally regulated by BRCA1, we used a microarray analysis of polysome-bound mRNAs in BRCA1-depleted and non-depleted MCF7 cells. Our findings show that BRCA1 modifies the translational efficiency of approximately 7% of the mRNAs expressed in these cells. Further analysis revealed that several processes contributing to cell surveillance such as cell cycle arrest, cell death, cellular growth and proliferation, DNA repair and gene expression, are largely enriched for the mRNAs whose translation is impacted by BRCA1. The BRCA1-dependent translation of these species of mRNAs therefore uncovers a novel mechanism through which BRCA1 exerts its onco-suppressive role. In addition, the BRCA1-dependent translation of mRNAs participating in unexpected functions such as cellular movement, nucleic acid metabolism or protein trafficking is indicative of novel functions for BRCA1. Finally, this study contributes to the identification of several markers associated with BRCA1 deficiency and to the discovery of new potential anti-neoplastic therapeutic targets.

Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer

TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ∼480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10(-7)), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10(-8)) and BRCA1 mutation carrier (P = 1.1 × 10(-5)) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10(-14)), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10(-15)) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10(-12)) and BRCA1 mutation carrier (P = 1.6 × 10(-14)) breast and invasive ovarian (P = 1.3 × 10(-11)) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant.

Identification of a BRCA2-specific modifier locus at 6p24 related to breast cancer risk

Common genetic variants contribute to the observed variation in breast cancer risk for BRCA2 mutation carriers; those known to date have all been found through population-based genome-wide association studies (GWAS). To comprehensively identify breast cancer risk modifying loci for BRCA2 mutation carriers, we conducted a deep replication of an ongoing GWAS discovery study. Using the ranked P-values of the breast cancer associations with the imputed genotype of 1.4 M SNPs, 19,029 SNPs were selected and designed for inclusion on a custom Illumina array that included a total of 211,155 SNPs as part of a multi-consortial project. DNA samples from 3,881 breast cancer affected and 4,330 unaffected BRCA2 mutation carriers from 47 studies belonging to the Consortium of Investigators of Modifiers of BRCA1/2 were genotyped and available for analysis. We replicated previously reported breast cancer susceptibility alleles in these BRCA2 mutation carriers and for several regions (including FGFR2, MAP3K1, CDKN2A/B, and PTHLH) identified SNPs that have stronger evidence of association than those previously published. We also identified a novel susceptibility allele at 6p24 that was inversely associated with risk in BRCA2 mutation carriers (rs9348512; per allele HR = 0.85, 95% CI 0.80-0.90, P = 3.9 × 10(-8)). This SNP was not associated with breast cancer risk either in the general population or in BRCA1 mutation carriers. The locus lies within a region containing TFAP2A, which encodes a transcriptional activation protein that interacts with several tumor suppressor genes. This report identifies the first breast cancer risk locus specific to a BRCA2 mutation background. This comprehensive update of novel and previously reported breast cancer susceptibility loci contributes to the establishment of a panel of SNPs that modify breast cancer risk in BRCA2 mutation carriers. This panel may have clinical utility for women with BRCA2 mutations weighing options for medical prevention of breast cancer.

Genome-wide association study in BRCA1 mutation carriers identifies novel loci associated with breast and ovarian cancer risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7 × 10(-8), HR = 1.14, 95% CI: 1.09-1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4 × 10(-8), HR = 1.27, 95% CI: 1.17-1.38) and 4q32.3 (rs4691139, P = 3.4 × 10(-8), HR = 1.20, 95% CI: 1.17-1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10(-4)). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%-50% compared to 81%-100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.