Skip to main content
  • Research article
  • Open access
  • Published:

Exploring the link between MORF4L1 and risk of breast cancer

Abstract

Introduction

Proteins encoded by Fanconi anemia (FA) and/or breast cancer (BrCa) susceptibility genes cooperate in a common DNA damage repair signaling pathway. To gain deeper insight into this pathway and its influence on cancer risk, we searched for novel components through protein physical interaction screens.

Methods

Protein physical interactions were screened using the yeast two-hybrid system. Co-affinity purifications and endogenous co-immunoprecipitation assays were performed to corroborate interactions. Biochemical and functional assays in human, mouse and Caenorhabditis elegans models were carried out to characterize pathway components. Thirteen FANCD2-monoubiquitinylation-positive FA cell lines excluded for genetic defects in the downstream pathway components and 300 familial BrCa patients negative for BRCA1/2 mutations were analyzed for genetic mutations. Common genetic variants were genotyped in 9,573 BRCA1/2 mutation carriers for associations with BrCa risk.

Results

A previously identified co-purifying protein with PALB2 was identified, MRG15 (MORF4L1 gene). Results in human, mouse and C. elegans models delineate molecular and functional relationships with BRCA2, PALB2, RAD51 and RPA1 that suggest a role for MRG15 in the repair of DNA double-strand breaks. Mrg15-deficient murine embryonic fibroblasts showed moderate sensitivity to γ-irradiation relative to controls and reduced formation of Rad51 nuclear foci. Examination of mutants of MRG15 and BRCA2 C. elegans orthologs revealed phenocopy by accumulation of RPA-1 (human RPA1) nuclear foci and aberrant chromosomal compactions in meiotic cells. However, no alterations or mutations were identified for MRG15/MORF4L1 in unclassified FA patients and BrCa familial cases. Finally, no significant associations between common MORF4L1 variants and BrCa risk for BRCA1 or BRCA2 mutation carriers were identified: rs7164529, Ptrend = 0.45 and 0.05, P2df = 0.51 and 0.14, respectively; and rs10519219, Ptrend = 0.92 and 0.72, P2df = 0.76 and 0.07, respectively.

Conclusions

While the present study expands on the role of MRG15 in the control of genomic stability, weak associations cannot be ruled out for potential low-penetrance variants at MORF4L1 and BrCa risk among BRCA2 mutation carriers.

Introduction

Genes that when mutated cause Fanconi anemia (FA) and/or influence breast cancer (BrCa) susceptibility functionally converge on a homology-directed DNA damage repair process [1]. That is, 15 FA genes (FANCs) and genes with high-penetrance, moderate-penetrance or low-penetrance mutations for BrCa encode for proteins cooperating in a defined FA/BrCa signaling pathway [2–6]. Remarkably, germline bi-allelic and mono-allelic loss-of-function mutations in four of these genes cause FA and BrCa, respectively: FANCD1/BRCA2 [7, 8], FANCJ/BRIP1 [9–12], FANCN/PALB2 [13–15], and the recently identified FA-like/BrCa mutated gene FANCO/RAD51C [3, 4]. These observations partially endorse perturbation of the DNA damage response as fundamental in leading to breast carcinogenesis. In addition to the main effects on susceptibility, variation in RAD51 - a gene encoding for a component of this pathway and paralog of RAD51C - modifies BrCa risk among BRCA2 but not BRCA1 mutation carriers [16]. Notably, RAD51 interacts with BRCA1 and BRCA2 [17, 18] to regulate double-strand breaks repair by homologous recombination [19].

While genes with low-penetrance and/or modifier alleles can be linked to diverse biological processes, the FA/BrCa pathway is still incomplete [2, 20]. To gain deeper insight into the molecular and functional FA/BrCa wiring diagram and the fundamental biological process(es) influencing cancer risk, we screened for novel protein physical interactions of known pathway components. Consistent with previous results on protein complex memberships [21, 22], we identified a physical interaction between PALB2 and MRG15. Results from the analysis of MRG15/MORF4L1 in unclassified FA patients and familial BrCa cases did not reveal pathological alterations; nonetheless, a weak modifier effect among carriers of BRCA2 mutations cannot be ruled out.

Materials and methods

Yeast two-hybrid design and screens

Following indications of increased sensitivity in the yeast two-hybrid (Y2H) system [23, 24], we designed multiple baits of each FA/BrCa pathway protein according to family domains defined by Pfam [25] and intrinsically disordered regions predicted by PONDR [26], as well as full-length ORFs. Proteome-scale Y2H screens were carried out using the mating strategy [27] and two different cDNA libraries as sources of prey, of human fetal brain or spleen (ProQuest; Invitrogen, Carlsbad, CA, USA). Bait fragments were obtained by RT-PCR using cDNAs derived from healthy lymphocytes, with the primers indicated in Additional file 1 and were subsequently cloned into the Gateway pDONR201 (Invitrogen) vector. Baits were 5'-sequenced so that they were confirmed, they did not show changes relative to publicly available sequence information and they were in-frame. Fragments were then transferred to the pPC97 yeast expression vector (Invitrogen) to be fused with the DNA-binding domain of Gal4. Constructs were transformed into the AH109 (Clontech, Palo Alto, CA, USA) yeast strain for screens (Y187 mate strain) using selective medium lacking histidine and supplemented with 10 mM 3-amino-triazole (Sigma-Aldrich, Taufkirchen, Germany) to test the interaction-dependent transactivation of the HIS3 reporter. Baits had previously been examined for self-activation at 3-amino-triazole concentrations in the range 10 to 80 mM. Under these conditions, >107 transformants were screened for each bait. Positive colonies were grown in selective medium for three cycles (10 to 15 days) to avoid unspecific cDNA contaminants, prior to PCR amplification and sequence identification of prey [28].

Microarray data analysis

The similarity of expression profiles was evaluated by calculating Pearson correlation coefficients using normalized (gcRMA) expression levels from the Human GeneAtlas U133A dataset [29] [Gene Expression Omnibus:GSE1133]. Comparisons were made for all possible microarray probe pairs.

Co-immunoprecipitation and co-affinity purification assays

For co-affinity purification (co-AP) assays, plasmids (1.5 μg) were transfected into HEK293/HeLa cells in six-well format using Lipofectamine 2000 (Invitrogen). Cells were then cultured for 48 hours and lysates prepared in buffer containing 50 mM Tris-HCl (pH 7.5), 100 to 150 mM NaCl, 0.5% Nonidet P-40, 1 mM ethylenediamine tetraacetic acid, and protease inhibitor cocktail (Roche Molecular Biochemicals, Indianapolis, IN, USA). Lysates were clarified twice by centrifugation at 13,000 × g before purification of protein complexes using sepharose beads (GE Healthcare, Piscataway, NJ, USA) for 1 hour at 4°C. Purified complexes and control lysate samples were resolved in Tris-glycine SDS-PAGE gels, then transferred to Invitrolon PVDF membranes (Invitrogen) or IMMOBILON PVDF (Millipore Corporation, Billerica, MA, USA), and target proteins were identified by detection of horseradish peroxidase-labeled antibody complexes with chemiluminescence using the ECL/ECL-Plus Western Blotting Detection Kit (GE Healthcare) or the Pierce ECL Western Blotting Substrate (Thermo Fisher Scientific, Waltham, MA, USA) following standard protocols. In some cases, samples were resolved in NuPAGE Novex 4 to 12% Bis-Tris or 3 to 8% Tris-Acetate Gels (Invitrogen). GST/GST-importin co-APs were performed as previously described [30].

For endogenous co-immunoprecipitation (co-IP) assays, cell cultures were washed with PBS and lysed at 0.5 × 107 to 1 × 107 cells/ml in NETN buffers (20 mM Tris pH 7.5, 1 mM ethylenediamine tetraacetic acid and 0.5% NP-40) containing 100 to 350 mM NaCl plus protease inhibitor cocktail (Roche Molecular Biochemicals). In some assays, supplementary phosphatase (10 to 50 mM NaF) or proteasome (MG132; Sigma-Aldrich) inhibitors were added to the solutions. Lysates were pre-cleared with protein-A sepharose beads (GE Healthcare), incubated with antibodies (2.5 to 5 μg) for 2 hours to overnight at 4°C with rotation, and then with protein-A beads for 1 hour at 4°C with rotation. Beads were collected by centrifugation and washed four times with lysis buffer prior to gel analysis.

Survival and iRNA-based assays

For evaluation of survival, 3 × 105 cells were seeded in duplicate in 60-mm dishes and left to recover for 24 hours. Cultures were then exposed to mitomycin-C or γ-radiation at the indicated doses. Next, 72 hours after the treatment, cells were rinsed with PBS, harvested by trypsinization and counted. Survival is reported as the percentage relative to untreated controls. Each siRNA (Additional file 2) was transfected for two successive rounds (24 hours apart) at a final concentration of 20 nM using Lipofectamine RNAiMAX reagent (Invitrogen) according to the manufacturer's instructions. After 4 days, cultures were treated with mitomycin-C or γ-radiation. Stealth siRNA Lo GC (12935-200; Invitrogen) was used as a negative control.

Immunofluorescence microscopy and antibodies

Cells were grown on glass cover slips and fixed using standard paraformaldehyde solution. Pre-extraction with PBS containing 0.5% Triton X-100 for 5 minutes at room temperature was used in some experiments. Staining was performed overnight at 4°C using appropriate primary antibody dilutions. Samples were then washed three times with 0.02% Tween 20 in PBS, incubated for 30 minutes at room temperature with Alexa fluor-conjugated secondary antibodies (Molecular Probes, Invitrogen), washed three times with 0.02% Tween 20 in PBS, and mounted on 4,6-diamidino-2-phenylindole-containing VECTASHIELD solution (Vector Laboratories, Peterborough, UK). Images were obtained using a Leica CTR-6000 microscope (Leica, Buffalo Grove, IL, USA).

Purified negative control IgGs of different species were purchased from Santa Cruz Biotechnology, Inc. (Santa Cruz, CA, USA). Anti-tag antibodies used were anti-HA (12CA5 and Y11; Santa Cruz Biotechnology), anti-HIS (H15; Santa Cruz Biotechnology) and anti-MYC (9E10; Sigma-Aldrich). Other antibodies used were anti-ACTN (ACTN05 C4; Abcam, Cambridge, UK), anti-Actb (8226; Abcam), anti-ATR (09-070; Millipore), anti-BRCA2 (Ab-1; Calbiochem-EMD Biosciences, San Diego, CA, USA), anti-CHEK2 (H300; Santa Cruz Biotechnology), anti-CHUK (ab54626; Abcam), anti-FANCD2 (ab2187; Abcam), anti-phospho-Ser139-H2AX (JBW301; Millipore), anti-KPNA1 (ab6035 and ab55387; Abcam), anti-MRG15 (N2-14; Novus Biologicals, Littleton, CO, USA; 1-235 ab37602; Abcam; and 15C [31–34]), anti-NFKB1 (H119; Santa Cruz Biotechnology), anti-p84 (ab487; Abcam), anti-PALB2 (675-725; Novus Biologicals), anti-PPHLN1 (ab69569; Abcam), anti-RAD51 (H92; Santa Cruz Biotechnology), anti-RPA1 (C88375; LifeSpan BioSciences, Seattle, WA, USA), anti-TOP3A (N20; Santa Cruz Biotechnology), anti-TRF2 (36; BD Transduction Laboratories, Mississauga, ON, USA), anti-TSNAX (3179C2a; Santa Cruz Biotechnology), and anti-USP1 (AP130a; Abgent, San Diego, CA, USA). Secondary horseradish peroxidase-linked antibodies were purchased from GE Healthcare and Abcam.

Caenorhabditis elegansstudies

Worms were cultured according to standard protocols, maintained on NGM agar seeded with Escherichia Coli OP50 [35]. The Bristol N2 strain was used as the wild-type strain. Strains carrying mutations studied here were provided by the Caenorhabditis Genetics Center (University of Minnesota, Minneapolis, MN, USA): DW104 brc-2(tm1086) III/hT2[bli-4(e937) let-?(q782) qIs48](I;III); VC1873: rad-51(ok2218) IV/nT1[qIs51](IV;V); and XA6226 mrg-1(qa6200)/qC1 dpy-19(e1259) glp-1(q339)[qIs26]. Gonads from gravid adults were dissected out with fine-gauge needles to perform a standard immunofluorescence. Primary antibodies were rat anti-RPA-1 (1:500) and rabbit anti-RAD-51 (1:100). Secondary antibodies were anti-rat Alexa 488 and anti-rabbit Alexa 568 (Invitrogen). Gonads were mounted with ProLong® Gold antifade reagent with 4,6-diamidino-2-phenylindole (Invitrogen). The cell-permeable SYTO 12 Green-Fluorescent Nucleic Acid Stain (Invitrogen) was used to label apoptotic cell death.

Study samples, genotyping and statistical analysis

All participants were enrolled under Institutional Review Boards or ethics committee approval at each participating center, and gave written informed consent. Research was conducted in accordance with the Declaration of Helsinki.

The MORF4L1 genomic sequence was obtained from the University of California at Santa Cruz Genome Browser version hg18 and intronic primers were designed using the web-based program Primer3 [36]. Extracts from 13 unclassified FANCD2 monoubiquitinylation-proficient FA cell lines, without mutations in FANCJ, FANCD1, FANCN, FANCO, or FANCP, and including six cases with deficient RAD51 nuclear foci formation, were examined by immunoblotting using the anti-MRG15 15C antibody [31–34]. These samples were also sequenced on all annotated MORF4L1 exons and exon-intron boundaries using primers shown in Additional file 3.

BRCA1 and BRCA2 mutation carriers were enrolled through 18 centers participating in the CIMBA and following previously detailed criteria [37, 38]. The following individual and clinical data were collected: year of birth, mutation description, ethnicity, country of residence, age at last follow-up, age at diagnosis of BrCa or at ovarian cancer diagnosis, age at bilateral prophylactic mastectomy, and age at bilateral prophylactic oophorectomy.

Genotyping was performed at the corresponding centers using 5' to 3' nuclease-based assays (TaqMan; Applied Biosystems, Foster City, CA, USA), except for an iPLEX assay carried out at the Queensland Institute of Medical Research (Brisbane, Australia) and containing EMBRACE, FCCC, GEORGETOWN, HEBCS, HEBON, ILUH, kConFab, Mayo Clinic, PBCS, SWE-BRCA and UPENN carriers. Results of these assays were centralized and analyzed for quality control as previously described [37]. Based on these criteria, one study was excluded from the analysis.

Hazard ratio (HR) estimates were obtained using Cox regression models under both standard regression analysis and under a weighted cohort approach to allow for the retrospective study design and the nonrandom sampling of affected and unaffected mutation carriers [39]. Analyses were stratified by birth cohort (<1940, 1940 to 1949, 1950 to 1959 and ≥1960), ethnicity and study center. A robust variance estimate was used to account for familial correlations. Time to diagnosis of BrCa from birth was modeled by censoring at the first of the following events: bilateral prophylactic mastectomy, BrCa diagnosis, ovarian cancer diagnosis, death and last date known to be alive. Subjects were considered affected if they were censored at BrCa diagnosis and unaffected otherwise. The weighted cohort approach involves assigning weights separately to affected and unaffected individuals such that the weighted observed incidences in the sample agree with established estimates for mutation carriers [39]. This approach has been shown to adjust for the bias in the HR estimates resulting from the ascertainment criteria used, which leads to an oversampling of affected women. Weights were assigned separately for carriers of mutations in BRCA1 and BRCA2 and by age interval (<25, 25 to 29, 30 to 34, 35 to 39, 40 to 44, 45 to 49, 50 to 54, 55 to 59, 60 to 64, 65 to 69 and ≥70). P values were derived from the robust score test.

Results

Protein physical interactions

The Y2H system was used to identify physical interactions for components of the FA/BrCa signaling pathway. In an initial phase, we screened for interactors of 12 proteins, which included the products of the FANCJ and FANCN genes (BRIP1 and PALB2, respectively) [9–11, 15], CHEK2 as linked to BrCa risk [40], and known molecular and/or functional interactors of FA/BrCa proteins (ATR, BLM, ERCC1, ERCC4, H2AFX, RAD51, TOP3A, TOPBP1 and USP1; see Additional file 1). To increase interactome coverage, we used specific protein domains or defined regions as baits, in addition to full-length ORFs, and screened >107 transformants belonging to two different cDNA sources (see Materials and methods). Multiple baits were thus screened for each protein based on Pfam-based family domain similarities [25] and on predicted intrinsically disordered regions using the PONDR algorithm [26]. Intrinsically disordered regions are defined as lacking a fixed tertiary structure and appear to be more common in nuclear proteins and involved in the cell cycle, transcription and signaling regulation processes [41, 42]. A total of 33 baits were screened for the 12 target proteins (Additional file 1).

Two previously demonstrated and six novel, potential physical interactions were identified through the Y2H screens (Additional file 4). Consistent with the physical interaction between their products, analysis of transcriptomic data identified significant expression correlations across normal human samples for most gene pairs (Additional file 5). The known interactions were BLM-MLH1 [43] and ERCC4-ERCC1 [44], through a predicted disordered region and a family domain, respectively (Additional file 6). The potential physical interactions included a previously described protein complex membership between PALB2 and MRG15 (also known as the MORF4-like 1 gene product) [21, 22]. To corroborate the Y2H results, we performed co-AP and co-IP assays, which suggested reliability for four of the interactions: CHEK2-NFKB1, PALB2-MRG15, TOP3A-TSNAX and USP1-KPNA1 (Additional file 7). TOP3A was originally co-purified with, among others, BLM, FANCA and replication proteins [45]. TSNAX (also known as translin (TSN)-associated factor X) was previously found to interact physically with MORF4 family associated protein 1-like 1 [46], and USP1 and KPNA1 were co-purified [47]. With the exception of MRG15 (see below), however, protein depletion assays did not show cellular sensitivity to γ-irradiation or mitomycin-C for any of the potential pathway components (siRNAs detailed in Additional file 2).

MRG15 is a chromo domain-containing protein present in histone acetyltransferase and deacetylase complexes [34], and the MRG15 ortholog in Drosophila melanogaster has been co-purified in histone chaperone complexes with a known BRCA2 interactor in humans, EMSY [48]. Consistent with a potential role in DNA damage repair, EAF3, the MORF family ortholog in Saccharomyces cerevisiae, was shown to interact genetically with radiation-sensitive (RAD) genes [49]. As previously shown [21, 22], MRGX, a close homolog of MRG15, also co-purified with PALB2 (Additional file 8). Consistent with the interaction domains delineated by the Y2H results, a MRG15 mutant lacking the C-terminal leucine zipper domain but not the N-terminal chromo domain was unable to interact with PALB2 (Additional file 8). Similarly, the helix-loop-helix region in MRGX was necessary for co-purification with PALB2 (Additional file 8). Together, these results support the identification of a physical interaction between PALB2 and MRG15, and probably MRGX.

MRG15 and DNA damage repair

According to the putative role of MRG15 in the repair of DNA double-strand breaks, murine embryonic fibroblasts (MEFs) derived from littermate embryos with the Morf4l1-/- genotype showed greater sensitivity (measured as cellular survival) to γ-irradiation than wild-type controls (Figure 1). The level of radiation sensitivity was moderate when compared with Atm-deficient MEFs (Figure 1). Milder sensitivity to mitomycin-C of cell cultures depleted of MRG15, relative to BRCA2 and PALB2, was also previously described [21]. In our study, however, deficiency of Mrg15 and depletion of MRG15 in MEFs and in HeLa and MCF10A cells, respectively, did not lead to a statistically significant increase in mitomycin-C-induced cell death or to G2/M phase cell cycle arrest and FANCD2 monoubiquitinylation (Additional file 9 shows results for HeLa cells). The observed milder effect and the use of different cell types may explain the discrepancy regarding mitomycin-C sensitivity when MRG15/Mrg15 is fully or partially depleted.

Figure 1
figure 1

Mrg15 deficiency confers sensitivity to γ-radiation. Mrg15-null murine embryonic fibroblasts (MEFs) show intermediate sensitivity to γ-radiation relative to controls (WT, wild-type; Morf4l1-/-, Mrg15-deficient; and Atm-/-, Atm-deficient). *Significant differences between WT and Morf4l1-/- MEFs (one-tailed t test, P < 0.01).

Contrary to the results for MRG15/Mrg15, radiation sensitivity phenotypes were not observed with assays for MRGX - also consistent with the previous study [21] - and for the potential novel interactor of TOP3A, TSNAX (data not shown). In agreement with the known role of TOP3A in telomere maintenance [50], however, an EmGFP-tagged TSNAX protein co-localized in specific nuclear structures with the telomere-binding protein TRF2 (Additional file 10). The major partner of TSNAX, TSN, was initially identified as a protein that binds to breakpoint junctions [51] and with high affinity to repeat sequences [52]. Although there is no evidence linking TSN to processes where recombination is necessary, there is some suggestion of a role in the DNA damage response [53]. Intriguingly, telomere shortening has been linked to FA pathology [54–56], and some FANC products were demonstrated to participate in telomere maintenance [57–59]. These observations lead to speculation that interactions between TSN-TSNAX-TOP3A may play a role in DNA damage repair and telomere maintenance by signaling through the FA/BrCa pathway.

In previous work, MRG15 appeared necessary for the association of BRCA2/PALB2/RAD51 with chromatin and the formation of nuclear foci following γ-irradiation [21]. In keeping with these observations, Morf4l1-/- MEFs showed lower numbers of Rad51 nuclear foci after γ-irradiation - discovered across time points and using clones or unselected cell cultures (Figure 2a shows results for clones). On the other hand, Morf4l1-/- MEFs showed lower expression levels of Brca1 and Brca2, but results were variable for Rad51 (Figure 2b) - Palb2 levels could not be assessed because the antibodies tested did not cross-react in mouse cell extracts. The result for Brca2 appeared to disagree with a previous study using human cell models [22]; however, another study showed reduction of BRCA2 through transient depletion of MRG15 but not MRGX [21]. This relationship for MRG15 could therefore be reminiscent of the role of PALB2 in stabilizing BRCA2 [60]. Together, these data suggest the involvement of MRG15 in the repair of DNA double-strand breaks through relationships with BRCA2, PALB2 and RAD51.

Figure 2
figure 2

Mrg15 deficiency impairs Rad51 foci formation and reduces Brca1 and Brca2 levels. (a) Left panel: number of cells with Rad51 nuclear foci (>4 foci per nuclei) in wild-type and Morf4l1-/- murine embryonic fibroblast (MEF) clones after (16 hours) treatment with 10 Gy. **Significant difference (two-tailed t test, P < 0.001). Right panel: representative images of Rad51 and pS139-H2ax immunodetection in cultures counted above for foci. DAPI, 4,6-diamidino-2-phenylindole; IR, γ-irradiated. (b) Levels of Brca1, Brca2 and Rad51, and control Actb, in whole cell extracts of Morf4l1-/- MEFs and wild-type counterparts (three cell clones of each genotype are shown).

Caenorhabditis elegansmutants of MRG15 and BRCA2 orthologs

The BRCA2 and RAD51 C. elegans orthologs (named BRC-2 and RAD-51, respectively) interact physically and regulate homologous recombination, so that brc-2 mutants fail to locate RAD-51 to sites of double-strand breaks present in meiosis or induced by DNA damage agents [61]. The hallmarks of brc-2 mutants in the germline are therefore lack of RAD-51 foci formation in parallel with an accumulation of RPA-1 at presumptive double-strand breaks, chromosomal abnormalities at diakinesis and, consequently, an increase in apoptotic corpses [61, 62]. C. elegans has an ortholog for the MORF human protein family (named MRG-1), which, like its mammalian counterparts, associates with chromatin and is required for embryo survival and cell proliferation [63, 64]. On the strength of this evidence, the functional link between BRC-2/BRCA2 and MRG-1/MRG15 was further investigated by assessing the phenocopy between brc-2 and mrg-1 mutants (tm1086 and qa6200, respectively).

Similar to brc-2 mutants, disruption of mrg-1 was linked to a remarkable increase in the number of RPA-1 foci in meiotic cells relative to wild-type animals (Figure 3a). While a wild-type animal presented, on average, three or four RPA-1 foci per nucleus, mrg-1 mutants commonly exhibited nuclei with more than 10 foci (Figure 3b). Two different patterns for RPA-1 staining were observed among mrg-1 mutant germ cell nuclei: one consisted of discrete foci similar to those observed in brc-2 mutants (Figure 3a, arrow), while the other showed more intense and diffuse staining (Figure 3a, arrowhead). Although RAD-51 staining was mainly nuclear in mrg-1 mutants - contrary to brc-2 mutants [61] - it was rather diffuse and often intense when compared with the usual pattern of discrete foci only observed in wild-type animals (Figure 3a and Additional file 11). Finally, mrg-1 mutants frequently showed aberrant chromosomal compaction (Figure 3a, asterisk) and, as expected, an increase in cell death revealed by SYTO-12 staining (Figure 3c). Together, these data further endorse the involvement of MRG-1/MRG15 in the control of genomic stability and suggest that perturbation of its function may activate the nonhomologous end-joining DNA damage repair process, as proposed for alteration of BRC-2 [61].

Figure 3
figure 3

Phenotypic study of Caenorhabditis elegans brc-2 and mrg-1 mutants. (a) Representative images of meiotic cells at the distal part, near the gonad bend. RAD-51 foci are bright and nuclear in wild-type (WT) animals whereas RAD-51 foci appear less intense and weakly diffuse in the cytoplasm, reduced but often dispersed and intense in the nuclei, or absent in brc-2, mrg-1 and rad-51 mutants, respectively. There are more RPA-1 nuclear foci in each of the three mutants than in WT animals. 4,6-Diamidino-2-phenylindole (DAPI) panels are merged with the red channel (for WT and brc-2 mutant) and with the green channel (for rad-51 mutant). *Abnormal chromosomal compaction. (b) Quantitation of RAD-51 and RPA-1 foci per nuclei in several germ cell lines of WT animals and brc-2 and mrg-1 mutant animals. Number of cells scored (n) and standard deviation of the mean indicated. **Significant differences relative to WT (Mann-Whitney U test, P < 0.001). (c) SYTO-12 staining in synchronized adult worms. Left top panel: an animal heterozygous for the brc-2 mutation (according to green fluorescent protein expression at the pharynx) shows WT SYTO-12 staining (that is, one to two labeled cells at the gonad bend). Right top and left bottom panels: an increase in SYTO-12-positive cells in the germline of brc-2 and mrg-1 mutants, respectively. Right bottom panel: magnification of the highlighted area in the left panel.

MORF4L1, Fanconi anemia and breast cancer risk

Having identified molecular and functional relationships for MRG15 in the repair of DNA double-strand breaks, we next evaluated the existence of alterations or mutations of MRG15/MORF4L1 in FA and BrCa patients. Immunoblotting of MRG15 using extracts of 13 FANCD2-monoubiquitinylation-positive FA cell lines - excluded for genetic defects in the downstream genes FANCD1/BRCA2, FANCJ/BRIP1, FANCN/PALB2, FANCO/RAD51C and FANCP/SLX4, and thus unclassifiable in terms of subtype - failed to show gross reduction of protein expression. This negative result included the analysis of six patient-derived FA cell lines defective for RAD51 foci (Additional file 12). Sequencing of MORF4L1 in these lines detected a few base substitutions and single base deletions deeper in the introns, and only annotated common variants in the exons (data not shown). Parallel to FA, we hypothesized that germline mutations or common variants in MORF4L1 may confer moderate/low risk of BrCa and/or modify cancer risk among BRCA1 and/or BRCA2 mutation carriers. Direct sequencing of MORF4L1 exons and flanking sequences in 300 patients with strong familial aggregation of BrCa but without detected mutations in BRCA1 or BRCA2, and belonging to two populations (United Kingdom, Institute of Cancer Research; Spain, Catalan Institute of Oncology), did not reveal pathogenic changes either. This negative result is consistent with a recent report in a similar setting by another group [65]. Nevertheless, given the extremely low frequency of high/moderate-penetrance mutations of other components of the FA/BrCa pathway [3, 12, 14] and the possible involvement in other cancer types [66], further investigation of MORF4L1 may be warranted.

The public results of the genome-wide association study conducted by the CGEMS initiative [67] suggest that common variation at the linkage disequilibrium block containing MORF4L1 is associated with BrCa risk (P2df < 0.01) (Figure 4a). Based on this observation, we genotyped two SNPs in a series of 9,573 BRCA1/2 mutation carriers collected through 18 centers participating in CIMBA: rs7164529 and rs10519219, with D' = 1 and r2 = 0.08. After quality control and Hardy-Weinberg equilibrium checks, Cox regression analysis revealed no significant associations between the SNPs and BrCa risk for BRCA1 or BRCA2 mutation carriers (rs7164529, Ptrend = 0.45 and 0.05, P2df = 0.51 and 0.14, respectively; rs10519219, Ptrend = 0.92 and 0.72, P2df = 0.76 and 0.07, respectively; Table 1). There was some suggestion of association with increased BrCa risk for BRCA2 mutation carriers under the recessive model for rs10519219 (P = 0.033) (Figure 4b and Additional file 13). Under the multiplicative model, there was no evidence of heterogeneity in the HRs of rs7164529 between studies (P = 0.66 and 0.21 for BRCA1 and BRCA2 mutation carriers, respectively) but some suggestion for rs10519219 among BRCA2 mutation carriers (P = 0.041). If an effect exists, the HR estimates for BRCA2 mutation carriers due to minor genotypes of rs7164529 or rs10519219 are in the opposite direction to those obtained in the general population (Table 1). Studying cancer susceptibility in mouse models has revealed opposite allele effects across different genetic backgrounds [68]. In this context, having a potential serial model of function between BRCA2 and MRG15, the effect of MORF4L1 alleles on BrCa risk might differ depending on the genetic/functional status of BRCA2/BRCA2: that is, wild-type in the general population versus altered or absent in BRCA2 mutation carriers. On the other hand, common predisposition alleles differentially associate with BrCa risk among BRCA1 and BRCA2 mutation carriers [16, 37, 69], which suggests differences in the influence of a given biological process on carcinogenesis between the two types of carriers.

Figure 4
figure 4

Variation at the MORF4L1 locus and breast cancer risk. (a) SNPs with previous suggestive evidence of association with breast cancer (BrCa) risk in the general population (P2df < 0.01) [67], genes and the linkage disequilibrium structure around MORF4L1 in HapMap Caucasians (data release 27). (b) Hazard ratio (HR) estimates of association of rs7164529 (top panels) and rs10519219 (bottom panels) with BrCa risk among BRCA1 (left panels) and BRCA2 (right panels) mutation carriers. Graphs show HRs and 95% confidence intervals of heterozygotes and minor allele homozygotes for all participating centers except for rs10519219 and relatively small groups (less than five individuals with the minor genotype). Size of the rectangle is proportional to the corresponding study precision.

Table 1 Association between variants at the MORF4L1 locus and breast cancer risk

We performed a number of sensitivity analyses to investigate the robustness of our results. Inclusion of prophylactic oophorectomy as a time-dependent covariate did not influence risk estimations (Pregression coefficients > 0.10). Some suggestion of association was revealed when prevalent cases, defined as those diagnosed >5 years before recruitment, were excluded from the analyses: rs7164529 per-allele model, BRCA2 n = 2,803, HR = 1.09, 95% confidence interval = 1.00 to 1.18, P = 0.048; and rs10519219 recessive model, BRCA2 n = 2,633, HR = 1.78, 95% confidence interval = 1.12 to 2.87, P = 0.027. Finally, data were also analyzed using a weighted cohort approach [39] to allow for the retrospective study design and, in particular, the nonrandom sampling of affected and unaffected mutation carriers. This yielded similar results to those shown in Table 1 for the per-allele and two-degrees-of-freedom models (rs7164529, BRCA1 weighted HR (wHR) = 1.04 to 1.08, BRCA2 wHR = 1.03 to 1.12; and rs10519219, BRCA1 wHR = 0.98 to 1.08, BRCA2 wHR = 0.95 to 1.59), but the rs10519219 association under the recessive model was no longer statistically significant (BRCA2 wHR = 1.62, 95% confidence interval = 0.97 to 2.70, P = 0.062) (Additional file 13). No evidence of heterogeneity was observed in any case for the wHRs (P > 0.30).

Discussion

Given the evidence across biological levels and species models, we hypothesized that perturbation of MRG15 function through genetic mutations or common alleles might be at the root of some cases of FA and/or BrCa. The results of our study, in addition to a recent publication on BrCa [65], indicate that in all probability the germline mutations in MORF4L1, if any, are not at the root of FA or BrCa. Next, analysis of common genetic variation at the MORF4L1 locus in BRCA1 and BRCA2 mutation carriers has not identified significant associations under the principal models. However, weak associations for risk among the latter group under the additive (rs7164529) and recessive (rs10519219) models might exist. Notably, in addition to the molecular and functional data presented, while MRG15 was demonstrated to co-purify with both BRCA1 and BRCA2, it only appeared to be necessary for the recruitment of BRCA2 (and PALB2/RAD51), but not of BRCA1, at sites of DNA damage [21]. Taken together, these observations suggest that the potential link between MORF4L1 and risk of BrCa warrants further assessment in larger sets of BRCA2 mutations and in additional case-control studies.

Conclusions

Studies in human, mouse and C. elegans models link MRG15 to the repair of DNA double-strand breaks, possibly through molecular and/or functional interactions with BRCA2, PALB2, RAD51 and RPA1. No pathogenic alterations of MRG15 or MORF4L1 have been observed in FA patients unclassified in terms of subtype or in familial BrCa cases negative for mutations in BRCA1 or BRCA2. Finally, no significant association with BrCa risk among BRCA1 and BRCA2 mutation carriers has been revealed for two common genetic variants at the MORF4L1 locus. Given a potentially weak and specific effect among BRCA2 mutation carriers, however, analyses in a larger series may be warranted.

Abbreviations

BrCa:

breast cancer

CGEMS:

Cancer and Genetics Markers of Susceptibility

CIMBA:

Consortium of Investigators of Modifiers of BRCA1/2

co-AP:

co-affinity purification

co-IP:

co-immunoprecipitation

df:

degrees of freedom

EMBRACE:

Epidemiological Study of BRCA1 and BRCA2 Mutation Carriers

FA:

Fanconi anemia

FCCC:

Fox Chase Cancer Center

GEORGETOWN:

Georgetown University

HEBCS:

Helsinki Breast Cancer Study

HEBON:

Hereditary Breast and Ovarian Cancer Research Group Netherlands

HR:

hazard ratio

ILUH:

Iceland Landspitali - University Hospital

iRNA:

interfering RNA

kConFab:

Kathleen Cuningham Foundation Consortium for Research into Familial Breast Cancer

MEF:

murine embryonic fibroblast

MORF:

mortality factor

ORF:

open reading frame

PBCS:

Pisa Breast Cancer Study

PBS:

phosphate-buffered saline

PCR:

polymerase chain reaction

RPA:

replication protein

RT:

reverse transcription

siRNA:

small interfering RNA

SNP:

single nucleotide polymorphism

SWE-BRCA:

Swedish Breast Cancer

TSN:

translin

UPENN:

University of Pennsylvania

wHR:

weighted hazard ratio

Y2H:

yeast two-hybrid.

References

  1. Wang W: Emergence of a DNA-damage response network consisting of Fanconi anaemia and BRCA proteins. Nat Rev Genet. 2007, 8: 735-748.

    Article  CAS  PubMed  Google Scholar 

  2. D'Andrea AD: Susceptibility pathways in Fanconi's anemia and breast cancer. N Engl J Med. 2010, 362: 1909-1919. 10.1056/NEJMra0809889.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Meindl A, Hellebrand H, Wiek C, Erven V, Wappenschmidt B, Niederacher D, Freund M, Lichtner P, Hartmann L, Schaal H, Ramser J, Honisch E, Kubisch C, Wichmann HE, Kast K, Deissler H, Engel C, Muller-Myhsok B, Neveling K, Kiechle M, Mathew CG, Schindler D, Schmutzler RK, Hanenberg H: Germline mutations in breast and ovarian cancer pedigrees establish RAD51C as a human cancer susceptibility gene. Nat Genet. 2010, 42: 410-414. 10.1038/ng.569.

    Article  CAS  PubMed  Google Scholar 

  4. Vaz F, Hanenberg H, Schuster B, Barker K, Wiek C, Erven V, Neveling K, Endt D, Kesterton I, Autore F, Fraternali F, Freund M, Hartmann L, Grimwade D, Roberts RG, Schaal H, Mohammed S, Rahman N, Schindler D, Mathew CG: Mutation of the RAD51C gene in a Fanconi anemia-like disorder. Nat Genet. 2010, 42: 406-409. 10.1038/ng.570.

    Article  CAS  PubMed  Google Scholar 

  5. Crossan GP, van der Weyden L, Rosado IV, Langevin F, Gaillard PH, McIntyre RE, Gallagher F, Kettunen MI, Lewis DY, Brindle K, Arends MJ, Adams DJ, Patel KJ: Disruption of mouse Slx4, a regulator of structure-specific nucleases, phenocopies Fanconi anemia. Nat Genet. 2011, 43: 147-152. 10.1038/ng.752.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Kim Y, Lach FP, Desetty R, Hanenberg H, Auerbach AD, Smogorzewska A: Mutations of the SLX4 gene in Fanconi anemia. Nat Genet. 2011, 43: 142-146. 10.1038/ng.750.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Howlett NG, Taniguchi T, Olson S, Cox B, Waisfisz Q, De Die-Smulders C, Persky N, Grompe M, Joenje H, Pals G, Ikeda H, Fox EA, D'Andrea AD: Biallelic inactivation of BRCA2 in Fanconi anemia. Science. 2002, 297: 606-609. 10.1126/science.1073834.

    Article  CAS  PubMed  Google Scholar 

  8. Wooster R, Bignell G, Lancaster J, Swift S, Seal S, Mangion J, Collins N, Gregory S, Gumbs C, Micklem G: Identification of the breast cancer susceptibility gene BRCA2. Nature. 1995, 378: 789-792. 10.1038/378789a0.

    Article  CAS  PubMed  Google Scholar 

  9. Levran O, Attwooll C, Henry RT, Milton KL, Neveling K, Rio P, Batish SD, Kalb R, Velleuer E, Barral S, Ott J, Petrini J, Schindler D, Hanenberg H, Auerbach AD: The BRCA1-interacting helicase BRIP1 is deficient in Fanconi anemia. Nat Genet. 2005, 37: 931-933. 10.1038/ng1624.

    Article  CAS  PubMed  Google Scholar 

  10. Levitus M, Waisfisz Q, Godthelp BC, de Vries Y, Hussain S, Wiegant WW, Elghalbzouri-Maghrani E, Steltenpool J, Rooimans MA, Pals G, Arwert F, Mathew CG, Zdzienicka MZ, Hiom K, De Winter JP, Joenje H: The DNA helicase BRIP1 is defective in Fanconi anemia complementation group J. Nat Genet. 2005, 37: 934-935. 10.1038/ng1625.

    Article  CAS  PubMed  Google Scholar 

  11. Litman R, Peng M, Jin Z, Zhang F, Zhang J, Powell S, Andreassen PR, Cantor SB: BACH1 is critical for homologous recombination and appears to be the Fanconi anemia gene product FANCJ. Cancer Cell. 2005, 8: 255-265. 10.1016/j.ccr.2005.08.004.

    Article  CAS  PubMed  Google Scholar 

  12. Seal S, Thompson D, Renwick A, Elliott A, Kelly P, Barfoot R, Chagtai T, Jayatilake H, Ahmed M, Spanova K, North B, McGuffog L, Evans DG, Eccles D, Easton DF, Stratton MR, Rahman N: Truncating mutations in the Fanconi anemia J gene BRIP1 are low-penetrance breast cancer susceptibility alleles. Nat Genet. 2006, 38: 1239-1241. 10.1038/ng1902.

    Article  CAS  PubMed  Google Scholar 

  13. Erkko H, Xia B, Nikkila J, Schleutker J, Syrjakoski K, Mannermaa A, Kallioniemi A, Pylkas K, Karppinen SM, Rapakko K, Miron A, Sheng Q, Li G, Mattila H, Bell DW, Haber DA, Grip M, Reiman M, Jukkola-Vuorinen A, Mustonen A, Kere J, Aaltonen LA, Kosma VM, Kataja V, Soini Y, Drapkin RI, Livingston DM, Winqvist R: A recurrent mutation in PALB2 in Finnish cancer families. Nature. 2007, 446: 316-319. 10.1038/nature05609.

    Article  CAS  PubMed  Google Scholar 

  14. Rahman N, Seal S, Thompson D, Kelly P, Renwick A, Elliott A, Reid S, Spanova K, Barfoot R, Chagtai T, Jayatilake H, McGuffog L, Hanks S, Evans DG, Eccles D, Easton DF, Stratton MR: PALB2, which encodes a BRCA2-interacting protein, is a breast cancer susceptibility gene. Nat Genet. 2007, 39: 165-167. 10.1038/ng1959.

    Article  CAS  PubMed  Google Scholar 

  15. Reid S, Schindler D, Hanenberg H, Barker K, Hanks S, Kalb R, Neveling K, Kelly P, Seal S, Freund M, Wurm M, Batish SD, Lach FP, Yetgin S, Neitzel H, Ariffin H, Tischkowitz M, Mathew CG, Auerbach AD, Rahman N: Biallelic mutations in PALB2 cause Fanconi anemia subtype FA-N and predispose to childhood cancer. Nat Genet. 2007, 39: 162-164. 10.1038/ng1947.

    Article  CAS  PubMed  Google Scholar 

  16. Antoniou AC, Sinilnikova OM, Simard J, Léoné M, Dumont M, Neuhausen SL, Struewing JP, Stoppa-Lyonnet D, Barjhoux L, Hughes DJ, Coupier I, Belotti M, Lasset C, Bonadona V, Bignon YJ, Genetic Modifiers of Cancer Risk in BRCA1/2 Mutation Carriers Study (GEMO), Rebbeck TR, Wagner T, Lynch HT, Domchek SM, Nathanson KL, Garber JE, Weitzel J, Narod SA, Tomlinson G, Olopade OI, Godwin A, Isaacs C, Jakubowska A, Lubinski J, et al: RAD51 135G→C modifies breast cancer risk among BRCA2 mutation carriers: results from a combined analysis of 19 studies. Am J Hum Genet. 2007, 81: 1186-1200. 10.1086/522611.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Scully R, Chen J, Plug A, Xiao Y, Weaver D, Feunteun J, Ashley T, Livingston DM: Association of BRCA1 with Rad51 in mitotic and meiotic cells. Cell. 1997, 88: 265-275. 10.1016/S0092-8674(00)81847-4.

    Article  CAS  PubMed  Google Scholar 

  18. Sharan SK, Morimatsu M, Albrecht U, Lim DS, Regel E, Dinh C, Sands A, Eichele G, Hasty P, Bradley A: Embryonic lethality and radiation hypersensitivity mediated by Rad51 in mice lacking Brca2. Nature. 1997, 386: 804-810. 10.1038/386804a0.

    Article  CAS  PubMed  Google Scholar 

  19. Thorslund T, West SC: BRCA2: a universal recombinase regulator. Oncogene. 2007, 26: 7720-7730. 10.1038/sj.onc.1210870.

    Article  CAS  PubMed  Google Scholar 

  20. Alpi AF, Patel KJ: Monoubiquitylation in the Fanconi anemia DNA damage response pathway. DNA Repair (Amst). 2009, 8: 430-435. 10.1016/j.dnarep.2009.01.019.

    Article  CAS  Google Scholar 

  21. Hayakawa T, Zhang F, Hayakawa N, Ohtani Y, Shinmyozu K, Nakayama J, Andreassen PR: MRG15 binds directly to PALB2 and stimulates homology-directed repair of chromosomal breaks. J Cell Sci. 2010, 123: 1124-1130. 10.1242/jcs.060178.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Sy SM, Huen MS, Chen J: MRG15 is a novel PALB2 interacting factor involved in homologous recombination. J Biol Chem. 2009, 284: 21127-21131. 10.1074/jbc.C109.023937.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Pujana MA, Han JD, Starita LM, Stevens KN, Tewari M, Ahn JS, Rennert G, Moreno V, Kirchhoff T, Gold B, Assmann V, Elshamy WM, Rual JF, Levine D, Rozek LS, Gelman RS, Gunsalus KC, Greenberg RA, Sobhian B, Bertin N, Venkatesan K, Ayivi-Guedehoussou N, Sole X, Hernandez P, Lazaro C, Nathanson KL, Weber BL, Cusick ME, Hill DE, Offit K, et al: Network modeling links breast cancer susceptibility and centrosome dysfunction. Nat Genet. 2007, 39: 1338-1349. 10.1038/ng.2007.2.

    Article  CAS  PubMed  Google Scholar 

  24. Boxem M, Maliga Z, Klitgord N, Li N, Lemmens I, Mana M, de Lichtervelde L, Mul JD, van de Peut D, Devos M, Simonis N, Yildirim MA, Cokol M, Kao HL, de Smet AS, Wang H, Schlaitz AL, Hao T, Milstein S, Fan C, Tipsword M, Drew K, Galli M, Rhrissorrakrai K, Drechsel D, Koller D, Roth FP, Iakoucheva LM, Dunker AK, Bonneau R, et al: A protein domain-based interactome network for C. elegans early embryogenesis. Cell. 2008, 134: 534-545. 10.1016/j.cell.2008.07.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz HR, Ceric G, Forslund K, Eddy SR, Sonnhammer EL, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2008, 36: D281-D288. 10.1093/nar/gkn226.

    Article  CAS  PubMed  Google Scholar 

  26. Romero P, Obradovic Z, Dunker AK: Natively disordered proteins: functions and predictions. Appl Bioinformatics. 2004, 3: 105-113. 10.2165/00822942-200403020-00005.

    Article  CAS  PubMed  Google Scholar 

  27. Walhout AJ, Vidal M: High-throughput yeast two-hybrid assays for large-scale protein interaction mapping. Methods. 2001, 24: 297-306. 10.1006/meth.2001.1190.

    Article  CAS  PubMed  Google Scholar 

  28. Vidalain PO, Boxem M, Ge H, Li S, Vidal M: Increasing specificity in high-throughput yeast two-hybrid experiments. Methods. 2004, 32: 363-370. 10.1016/j.ymeth.2003.10.001.

    Article  CAS  PubMed  Google Scholar 

  29. Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Depping R, Steinhoff A, Schindler SG, Friedrich B, Fagerlund R, Metzen E, Hartmann E, Kohler M: Nuclear translocation of hypoxia-inducible factors (HIFs): involvement of the classical importin a/b pathway. Biochim Biophys Acta. 2008, 1783: 394-404. 10.1016/j.bbamcr.2007.12.006.

    Article  CAS  PubMed  Google Scholar 

  31. Chen M, Takano-Maruyama M, Pereira-Smith OM, Gaufo GO, Tominaga K: MRG15, a component of HAT and HDAC complexes, is essential for proliferation and differentiation of neural precursor cells. J Neurosci Res. 2009, 87: 1522-1531. 10.1002/jnr.21976.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Garcia SN, Kirtane BM, Podlutsky AJ, Pereira-Smith OM, Tominaga K: Mrg15 null and heterozygous mouse embryonic fibroblasts exhibit DNA-repair defects post exposure to gamma ionizing radiation. FEBS Lett. 2007, 581: 5275-5281. 10.1016/j.febslet.2007.10.017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Tominaga K, Kirtane B, Jackson JG, Ikeno Y, Ikeda T, Hawks C, Smith JR, Matzuk MM, Pereira-Smith OM: MRG15 regulates embryonic development and cell proliferation. Mol Cell Biol. 2005, 25: 2924-2937. 10.1128/MCB.25.8.2924-2937.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Pardo PS, Leung JK, Lucchesi JC, Pereira-Smith OM: MRG15, a novel chromodomain protein, is present in two distinct multiprotein complexes involved in transcriptional activation. J Biol Chem. 2002, 277: 50860-50866. 10.1074/jbc.M203839200.

    Article  CAS  PubMed  Google Scholar 

  35. Stiernagle T: Maintenance of C. elegans. WormBook. 2006, 1: 1-11.

    Google Scholar 

  36. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  37. Antoniou AC, Sinilnikova OM, McGuffog L, Healey S, Nevanlinna H, Heikkinen T, Simard J, Spurdle AB, Beesley J, Chen X, Neuhausen SL, Ding YC, Couch FJ, Wang X, Fredericksen Z, Peterlongo P, Peissel B, Bonanni B, Viel A, Bernard L, Radice P, Szabo CI, Foretova L, Zikan M, Claes K, Greene MH, Mai PL, Rennert G, Lejbkowicz F, Andrulis IL, et al: Common variants in LSP1, 2q35 and 8q24 and breast cancer risk for BRCA1 and BRCA2 mutation carriers. Hum Mol Genet. 2009, 18: 4442-4456. 10.1093/hmg/ddp372.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. CIMBA. [http://www.srl.cam.ac.uk/consortia/cimba/eligibility/eligibility.html]

  39. Antoniou AC, Goldgar DE, Andrieu N, Chang-Claude J, Brohet R, Rookus MA, Easton DF: A weighted cohort approach for analysing factors modifying disease risks in carriers of high-risk susceptibility genes. Genet Epidemiol. 2005, 29: 1-11. 10.1002/gepi.20074.

    Article  PubMed  Google Scholar 

  40. Meijers-Heijboer H, van den Ouweland A, Klijn J, Wasielewski M, de Snoo A, Oldenburg R, Hollestelle A, Houben M, Crepin E, van Veghel-Plandsoen M, Elstrodt F, van Duijn C, Bartels C, Meijers C, Schutte M, McGuffog L, Thompson D, Easton D, Sodha N, Seal S, Barfoot R, Mangion J, Chang-Claude J, Eccles D, Eeles R, Evans DG, Houlston R, Murday V, Narod S, Peretz T, et al: Low-penetrance susceptibility to breast cancer due to CHEK2(*)1100delC in noncarriers of BRCA1 or BRCA2 mutations. Nat Genet. 2002, 31: 55-59. 10.1038/ng879.

    Article  CAS  PubMed  Google Scholar 

  41. Ward JJ, Sodhi JS, McGuffin LJ, Buxton BF, Jones DT: Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. J Mol Biol. 2004, 337: 635-645. 10.1016/j.jmb.2004.02.002.

    Article  CAS  PubMed  Google Scholar 

  42. Lobley A, Swindells MB, Orengo CA, Jones DT: Inferring function using patterns of native disorder in proteins. PLoS Comput Biol. 2007, 3: e162-10.1371/journal.pcbi.0030162.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Pedrazzi G, Perrera C, Blaser H, Kuster P, Marra G, Davies SL, Ryu GH, Freire R, Hickson ID, Jiricny J, Stagljar I: Direct association of Bloom's syndrome gene product with the human mismatch repair protein MLH1. Nucleic Acids Res. 2001, 29: 4378-4386. 10.1093/nar/29.21.4378.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Park CH, Sancar A: Formation of a ternary complex by human XPA, ERCC1, and ERCC4(XPF) excision repair proteins. Proc Natl Acad Sci USA. 1994, 91: 5017-5021. 10.1073/pnas.91.11.5017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Meetei AR, Sechi S, Wallisch M, Yang D, Young MK, Joenje H, Hoatlin ME, Wang W: A multiprotein nuclear complex connects Fanconi anemia and Bloom syndrome. Mol Cell Biol. 2003, 23: 3417-3426. 10.1128/MCB.23.10.3417-3426.2003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Rual JF, Venkatesan K, Hao T, Hirozane-Kishikawa T, Dricot A, Li N, Berriz GF, Gibbons FD, Dreze M, Ayivi-Guedehoussou N, Klitgord N, Simon C, Boxem M, Milstein S, Rosenberg J, Goldberg DS, Zhang LV, Wong SL, Franklin G, Li S, Albala JS, Lim J, Fraughton C, Llamosas E, Cevik S, Bex C, Lamesch P, Sikorski RS, Vandenhaute J, Zoghbi HY, et al: Towards a proteome-scale map of the human protein-protein interaction network. Nature. 2005, 437: 1173-1178. 10.1038/nature04209.

    Article  CAS  PubMed  Google Scholar 

  47. Sowa ME, Bennett EJ, Gygi SP, Harper JW: Defining the human deubiquitinating enzyme interaction landscape. Cell. 2009, 138: 389-403. 10.1016/j.cell.2009.04.042.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Moshkin YM, Kan TW, Goodfellow H, Bezstarosti K, Maeda RK, Pilyugin M, Karch F, Bray SJ, Demmers JA, Verrijzer CP: Histone chaperones ASF1 and NAP1 differentially modulate removal of active histone marks by LID-RPD3 complexes during NOTCH silencing. Mol Cell. 2009, 35: 782-793. 10.1016/j.molcel.2009.07.020.

    Article  CAS  PubMed  Google Scholar 

  49. Collins SR, Miller KM, Maas NL, Roguev A, Fillingham J, Chu CS, Schuldiner M, Gebbia M, Recht J, Shales M, Ding H, Xu H, Han J, Ingvarsdottir K, Cheng B, Andrews B, Boone C, Berger SL, Hieter P, Zhang Z, Brown GW, Ingles CJ, Emili A, Allis CD, Toczyski DP, Weissman JS, Greenblatt JF, Krogan NJ: Functional dissection of protein complexes involved in yeast chromosome biology using a genetic interaction map. Nature. 2007, 446: 806-810. 10.1038/nature05649.

    Article  CAS  PubMed  Google Scholar 

  50. Temime-Smaali N, Guittat L, Wenner T, Bayart E, Douarre C, Gomez D, Giraud-Panis MJ, Londono-Vallejo A, Gilson E, Amor-Gueret M, Riou JF: Topoisomerase IIIα is required for normal proliferation and telomere stability in alternative lengthening of telomeres. Embo J. 2008, 27: 1513-1524.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Aoki K, Suzuki K, Sugano T, Tasaka T, Nakahara K, Kuge O, Omori A, Kasai M: A novel gene, translin, encodes a recombination hotspot binding protein associated with chromosomal translocations. Nat Genet. 1995, 10: 167-174. 10.1038/ng0695-167.

    Article  CAS  PubMed  Google Scholar 

  52. Jacob E, Pucshansky L, Zeruya E, Baran N, Manor H: The human protein translin specifically binds single-stranded microsatellite repeats, d(GT)n, and G-strand telomeric repeats, d(TTAGGG)n: a study of the binding parameters. J Mol Biol. 2004, 344: 939-950. 10.1016/j.jmb.2004.09.095.

    Article  CAS  PubMed  Google Scholar 

  53. Jaendling A, McFarlane RJ: Biological roles of translin and translin-associated factor-X: RNA metabolism comes to the fore. Biochem J. 2010, 429: 225-234. 10.1042/BJ20100273.

    Article  CAS  PubMed  Google Scholar 

  54. Leteurtre F, Li X, Guardiola P, Le Roux G, Sergere JC, Richard P, Carosella ED, Gluckman E: Accelerated telomere shortening and telomerase activation in Fanconi's anaemia. Br J Haematol. 1999, 105: 883-893. 10.1046/j.1365-2141.1999.01445.x.

    Article  CAS  PubMed  Google Scholar 

  55. Ball SE, Gibson FM, Rizzo S, Tooze JA, Marsh JC, Gordon-Smith EC: Progressive telomere shortening in aplastic anemia. Blood. 1998, 91: 3582-3592.

    Article  CAS  PubMed  Google Scholar 

  56. Callén E, Samper E, Ramírez MJ, Creus A, Marcos R, Ortega JJ, Olivé T, Badell I, Blasco MA, Surrallés J: Breaks at telomeres and TRF2-independent end fusions in Fanconi anemia. Hum Mol Genet. 2002, 11: 439-444. 10.1093/hmg/11.4.439.

    Article  PubMed  Google Scholar 

  57. Fan Q, Zhang F, Barrett B, Ren K, Andreassen PR: A role for monoubiquitinated FANCD2 at telomeres in ALT cells. Nucleic Acids Res. 2009, 37: 1740-1754. 10.1093/nar/gkn995.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Franco S, van de Vrugt HJ, Fernández P, Aracil M, Arwert F, Blasco MA: Telomere dynamics in Fancg-deficient mouse and human cells. Blood. 2004, 104: 3927-3935. 10.1182/blood-2003-10-3626.

    Article  CAS  PubMed  Google Scholar 

  59. Rhee DB, Wang Y, Mizesko M, Zhou F, Haneline L, Liu Y: FANCC suppresses short telomere-initiated telomere sister chromatid exchange. Hum Mol Genet. 2010, 19: 879-887. 10.1093/hmg/ddp556.

    Article  CAS  PubMed  Google Scholar 

  60. Xia B, Sheng Q, Nakanishi K, Ohashi A, Wu J, Christ N, Liu X, Jasin M, Couch FJ, Livingston DM: Control of BRCA2 cellular and clinical functions by a nuclear partner, PALB2. Mol Cell. 2006, 22: 719-729. 10.1016/j.molcel.2006.05.022.

    Article  CAS  PubMed  Google Scholar 

  61. Martin JS, Winkelmann N, Petalcorin MI, McIlwraith MJ, Boulton SJ: RAD-51-dependent and -independent roles of a Caenorhabditis elegans BRCA2-related protein during DNA double-strand break repair. Mol Cell Biol. 2005, 25: 3127-3139. 10.1128/MCB.25.8.3127-3139.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Ko E, Lee J, Lee H: Essential role of brc-2 in chromosome integrity of germ cells in C. elegans. Mol Cells. 2008, 26: 590-594.

    CAS  PubMed  Google Scholar 

  63. Fujita M, Takasaki T, Nakajima N, Kawano T, Shimura Y, Sakamoto H: MRG-1, a mortality factor-related chromodomain protein, is required maternally for primordial germ cells to initiate mitotic proliferation in C. elegans. Mech Dev. 2002, 114: 61-69. 10.1016/S0925-4773(02)00058-8.

    Article  CAS  PubMed  Google Scholar 

  64. Takasaki T, Liu Z, Habara Y, Nishiwaki K, Nakayama J, Inoue K, Sakamoto H, Strome S: MRG-1, an autosome-associated protein, silences X-linked genes and protects germline immortality in Caenorhabditis elegans. Development. 2007, 134: 757-767. 10.1242/dev.02771.

    Article  CAS  PubMed  Google Scholar 

  65. Rio Frio T, Haanpaa M, Pouchet C, Pylkas K, Vuorela M, Tischkowitz M, Winqvist R, Foulkes WD: Mutation analysis of the gene encoding the PALB2-binding protein MRG15 in BRCA1/2-negative breast cancer families. J Hum Genet. 2010, 55: 842-843. 10.1038/jhg.2010.112.

    Article  PubMed  CAS  Google Scholar 

  66. Jones S, Hruban RH, Kamiyama M, Borges M, Zhang X, Parsons DW, Lin JC, Palmisano E, Brune K, Jaffee EM, Iacobuzio-Donahue CA, Maitra A, Parmigiani G, Kern SE, Velculescu VE, Kinzler KW, Vogelstein B, Eshleman JR, Goggins M, Klein AP: Exomic sequencing identifies PALB2 as a pancreatic cancer susceptibility gene. Science. 2009, 324: 217-10.1126/science.1171202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Hunter DJ, Kraft P, Jacobs KB, Cox DG, Yeager M, Hankinson SE, Wacholder S, Wang Z, Welch R, Hutchinson A, Wang J, Yu K, Chatterjee N, Orr N, Willett WC, Colditz GA, Ziegler RG, Berg CD, Buys SS, McCarty CA, Feigelson HS, Calle EE, Thun MJ, Hayes RB, Tucker M, Gerhard DS, Fraumeni JF, Hoover RN, Thomas G, Chanock SJ: A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer. Nat Genet. 2007, 39: 870-874. 10.1038/ng2075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Jorgenson TC, Williams BR, Wendland A, Bilger A, Sandgren EP, Drinkwater NR: Identification of susceptibility loci in a mouse model of KRASG12D-driven pancreatic cancer. Cancer Res. 2010, 70: 8398-8406. 10.1158/0008-5472.CAN-09-3980.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Antoniou AC, Spurdle AB, Sinilnikova OM, Healey S, Pooley KA, Schmutzler RK, Versmold B, Engel C, Meindl A, Arnold N, Hofmann W, Sutter C, Niederacher D, Deissler H, Caldes T, Kampjarvi K, Nevanlinna H, Simard J, Beesley J, Chen X, Neuhausen SL, Rebbeck TR, Wagner T, Lynch HT, Isaacs C, Weitzel J, Ganz PA, Daly MB, Tomlinson G, Olopade OI, et al: Common breast cancer-predisposition alleles are associated with breast cancer risk in BRCA1 and BRCA2 mutation carriers. Am J Hum Genet. 2008, 82: 937-948. 10.1016/j.ajhg.2008.02.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank Dr Andre Nussenzweig, Dr Laura Tusell and Dr Anton Gartner for providing Atm-deficient MEFs, anti-TRF2 and anti-RAD-51/RPA-1 antibodies, respectively. They also wish to thank all study participants, clinicians and centers for their valuable contribution, and the CGEMS initiative for making their genome-wide association study results available. The CIMBA data management is supported by Cancer Research - UK. Funding support is further detailed in Additional file 14.

Author information

Authors and Affiliations

Authors

Consortia

Corresponding authors

Correspondence to Julián Cerón, Kaoru Tominaga, Jordi Surrallés or Miguel Angel Pujana.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

The project was conceived and the experiments and data analyses were coordinated by JS and MAP. The Y2H design and screens were performed by GM, CAM, LG-B and MAP. The co-AP/co-IP assays, biochemical and/or cell biology studies of FA/BrCa pathway components were performed by GM, CAM, LG-B, HA, FKP, RD and MAP. The studies of MORF4L1/MRG15 in MEFs and the co-AP assays were performed by ET, OMP-S and KT. The studies of mitomycin-C and γ-radiation sensitivity, and FANCD2 monoubiquitinylation were performed by MB, MJR, MC, GH and JS. Statistical analyses were performed by NB, DC and MAP with the support of LM and ACA. The studies in C. elegans were performed by MP and JC. MORF4L1 sequencing was carried out by SS, AR and NR in the United Kingdom, and CL, IB, JBr, JF-R and MAP in Spain. The study of cell lines from FA patients was performed by JK, KN and DS. The study of CIMBA carriers was coordinated and executed by DFE, LM, ACA and GC-T. iPLEX genotyping was performed by XC and JBee. Classification of BRCA1/2 mutations was performed by SH and OMS. DNA samples and clinical data of carriers were contributed by: DFE, SP, MC, CTO, DF, RP, DGE, FL, RE, LI, CC, RD, K-RO, JC, FD, SH, CB, PJM and MP (EMBRACE); PP, SM, BP, DF, GR, MB, AV, BP, LO, ALP, AS, LB and PR (CONSIT TEAM); SH, AS, XC, JB and GC-T (kConFab); MAR, SV, MAT-L, MPV, CJA, DB, MGEMA, TAO, MJB, HEJM-H and FBLH (HEBON); DEG, SB, EMJ, AM, JLH and MBD (BCFR); KH, AB, JR, GB-B, HE and MS-A (SWE-BRCA); BK, YL, RM and EF (SMC); SMD, KLN and TRB (UPENN); OTJ (ILUH); FJC, XW and ZF (Mayo); TC (HCSC); TH and HN (HEBCS); UH and DT (DKFZ); MAC (PBCS); AKG (FCCC); ENI, RJ, OMS, DS-L, SM, CV-P, LC, AP, Y-JB, NU, J-PP, PV, SFF, M-AC-R and IM (BFBOCC and GEMO Study Collaborators); and CL, IB and J Brunet (ICO). AO, JBen, JBu and VM helped with data analysis and interpretation, and contributed with reagents. The manuscript was written by MAP. All authors read and approved the final manuscript.

Griselda Martrat, Christopher A Maxwell, Emiko Tominaga, Montserrat Porta-de-la-Riva, Núria Bonifaci, Laia Gómez-Baldó, Massimo Bogliolo contributed equally to this work.

Electronic supplementary material

13058_2010_9505_MOESM1_ESM.xls

Additional file 1: Y2H baits for 12 proteins in the FA/BrCa signaling pathway. Supplementary Table 1 containing details of the design of Y2H baits for 12 proteins in the FA/BrCa signaling pathway. (XLS 28 KB)

13058_2010_9505_MOESM2_ESM.xls

Additional file 2: siRNAs used in the present study. Supplementary Table 2 containing details of the siRNAs used in the present study. (XLS 22 KB)

13058_2010_9505_MOESM3_ESM.xls

Additional file 3: Primers for sequencing of MORF4L1. Supplementary Table 3 containing details of primers used for sequencing of MORF4L1. (XLS 20 KB)

13058_2010_9505_MOESM4_ESM.xls

Additional file 4: FA/BrCa signaling pathway components. Supplementary Table 4 containing details of known and potential FA/BrCa signaling pathway components identified through Y2H screens. (XLS 32 KB)

13058_2010_9505_MOESM5_ESM.PDF

Additional file 5: Gene co-expression. Supplementary Figure 1 containing results of the gene co-expression analysis. (PDF 638 KB)

13058_2010_9505_MOESM6_ESM.PDF

Additional file 6: Four bait designs and Y2H results. Supplementary Figure 2 containing details of four bait designs and the Y2H results. (PDF 618 KB)

13058_2010_9505_MOESM7_ESM.PDF

Additional file 7: Co-AP and co-IP assays. Supplementary Figure 3 containing results of the co-AP and co-IP assays. (PDF 3 MB)

13058_2010_9505_MOESM8_ESM.PDF

Additional file 8: Co-AP assays involving MRG15 and MRGX. Supplementary Figure 4 containing results of co-AP assays involving MRG15 and MRGX. (PDF 1 MB)

13058_2010_9505_MOESM9_ESM.PDF

Additional file 9: siRNA-mediated depletion of MRG15 and FANCD2 monoubiquitinylation. Supplementary Figure 5 containing results of siRNA-mediated depletion of MRG15 and FANCD2 monoubiquitinylation. (PDF 1 MB)

13058_2010_9505_MOESM10_ESM.PDF

Additional file 10: TRF2 and TSNAX co-localization. Supplementary Figure 6 containing results of TRF2 and TSNAX co-localization. (PDF 2 MB)

13058_2010_9505_MOESM11_ESM.eps

Additional file 11: Immunodetection of RAD-51 and RPA-1. Supplementary Figure 7 containing results for immunodetection of RAD-51 and RPA-1 in wild-type animals and in brc-2 and mrg-1 C. elegans mutant animals. (EPS 980 KB)

13058_2010_9505_MOESM12_ESM.PDF

Additional file 12: MRG15 in extracts of unclassified FA cell lines. Supplementary Figure 8 containing results for the analysis of MRG15 in extracts of unclassified FA cell lines. (PDF 1 MB)

13058_2010_9505_MOESM13_ESM.xls

Additional file 13: BrCa risk estimates for rs7164529 and rs10519219. Supplementary Table 5 containing BrCa risk estimates (HR and wHR) for rs7164529 (additive model) and rs10519219 (recessive model) among BRCA2 mutation carriers across participating centers. (XLS 24 KB)

Additional file 14: Funding support. Supplementary document containing details of funding support. (DOC 50 KB)

Authors’ original submitted files for images

Rights and permissions

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Cite this article

Martrat, G., Maxwell, C.A., Tominaga, E. et al. Exploring the link between MORF4L1 and risk of breast cancer. Breast Cancer Res 13, R40 (2011). https://doi.org/10.1186/bcr2862

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/bcr2862

Keywords