Functional consequences of a rare missense BARD1 c.403G>A germline mutation identified in a triple-negative breast cancer patient

Incorporation of next-generation sequencing into clinical practice continues to expand the list of variants of unknown significance (VUS), making it challenging to appropriately interpret the clinical significance of such rare mutations in terms of pathogenicity and treatment options for the patient [1‐3]. Analysis of the whole-genome sequencing data from a triple-negative breast cancer (TNBC) patient’s germline DNA uncovered a c.403G>A (p.Asp135Asn) mutation in BARD1 (Additional file 1). The BRCA1-associated ring domain 1 (BARD1) protein is a binding partner of BRCA1 and is essential for DNA damage repair [4], which is reported to be associated with breast cancer susceptibility [2, 3]. In contrast to BRCA1 and BRCA2, for which numerous pathogenic mutations and benign variants have been identified [5], validated, and used for clinical decision-making, the clinical significance of a specific BARD1 mutation remains unclear and is classified as VUS [6]. The newly identified locus is located near the ring domain of BARD1 and highly conserved across different species (Fig. 1a). The germline variant was validated by orthogonal Sanger sequencing (Fig. 1b), and it was inherited from the father of the patient (Fig. 1c). This c.403G>A mutation was only reported in one case from 138,632 whole-exome and whole-genome sequences in the Genome Aggregation Database (gnomAD [7]) with no reports on its clinical relevance. We predicted the pathogenicity of the rare missense c.403G>A mutation using four computational tools with conflicting prediction results: tolerated by SIFT [8] and benign by PolyPhen 2 [9], while damaged by Mutation Assessor [10] and SNPs&GO [11].

Importantly, computational analysis of the genomics profiling from the patient’s tumor samples showed that the patient’s tumor exhibited extensive copy number changes, loss of heterozygosity, and large-scale structural variations across the entire genome (Fig. 1d), indicating a typical genomic mutational signature resulting from deficiency in homologous recombination DNA damage repair [12]. Moreover, the mRNA expression profile of the tumor clustered closely with the basal-like and immune-suppressed (BLIS) subgroup in a large cohort of 465 TNBCs from Fudan University Shanghai Cancer Center (FUSCC) (Fig. 1e) [13], a study partially inspired by the patient. In addition, the mutation profile of BARD1 in the FUSCC TNBC cohort is shown in Fig. 1f, indicating that the germline mutation rate in BARD1 is low and the mutation loci were heterogeneous. Taken together, we hypothesized that the c.403G>A mutation in BARD1 might be damaging by impairing the homologous recombination capacity of the cells.

To explore the in vitro functions of the c.403G>A (p.Asp135Asn) mutation in BARD1, reconstitution with wild-type BARD1 (BARD1^WTres) and mutant BARD1 (BARD1^D135Nres) was conducted in two BARD1 knockout breast cancer cells (T47D and MDA-MA-468) by the CRISPR-Cas9 system. As shown in Fig. 2a, the BARD1 protein was undetectable in the knockout cells, whereas comparable expression of the BARD1 protein was detected in rescued cells (BARD1^WTres and BARD1^D135Nres) by Western blot using anti-BARD1 antibodies.

Impairment of DNA damage repair was reflected by increased sensitivity of tumor cells to DNA damage. In the γ-irradiation (IR) induced DNA damage experiment, BARD1 knockout T47D and MDA-MB-468 cells showed significantly reduced clonogenic survival compared with control cells after irradiation (Fig. 2b, P< 0.01). As expected, ectopic expression of BARD1^WTres resulted in full rescue of the reduced clonogenic survival, whereas ectopic expression of BARD1^D135Nres showed partial rescue and the clonogenic survival was still significantly reduced compared with control cells (Fig. 2b). Similar results were observed in the DNA damage experiment with DNA inter-strand crosslinking agent mitomycin C (MMC). BARD1-deficient cells expressing BARD1^D135Nres were more sensitive to MMC than cells expressing wild-type BARD1^WTres as determined by the CCK8 cell proliferation assay (Fig. 2b). Moreover, BARD1-deficient cells expressing BARD1^D135Nres were also more sensitive to an investigational PARP inhibitor HS-10160 than cells expressing BARD1^WTres by the CCK8 cell proliferation assays (Fig. 2b).

Mechanically, BARD1 was thought to play an important role in DNA damage repair through direct interaction with BRCA1 by forming a BARD1-BRCA1 heterodimer [14]. However, immunoprecipitation (IP) assays showed that the c.403G>A mutation in BARD1 had little effect on the interaction between BARD1 and BRCA1 (Fig. 3a). Recently, Zhao et al. reported that the BARD1 region (123-261) is indispensable for the interaction between BARD1 and RAD51, and RAD51-mediated homologous DNA pairing [15]. This raises the possibility that the p.D135N germline mutation may affect the function of BRAD1 in DNA damage repair with RAD51-mediated homologous recombination pairing. On the other hand, we found that the interaction between BARD1^D135Nres and RAD51 was significantly weaker compared to that between BARD1^WTres and RAD51 in T47D and HEK293T cells after irradiation (Fig. 3b). Taken together, these data suggested that the c.403G>A mutation in BARD1 increased sensitivity to PARP inhibitor treatment in vitro, which may be caused by the reduction of its DNA damage-induced association with RAD51.

To investigate the functions of the c.403G>A mutation in BARD1 in vivo, two individual mouse lines (PM#5 and PM#8) with a point mutation c.718G>A in mice (p.Asp127Asn) corresponding to the c.G403A (p.Asp135Asn) mutation in human BARD1 were obtained using the CRISPR/Cas9 system. Mice with the c.718G>A point mutation (PM#5 and PM#8) were significantly more susceptible to irradiation than the wild-type mice (Fig. 4a, b). With a sub-lethal dose (7 Gy) (i.e., 10% of the LD50 for wild-type C57BL6J male mice) of total body irradiation (TBI), the time to reach 50% lethality for wild-type control mice was nearly 4 days or longer compared to the homozygotes point mutation lines, and the P values of log-rank tests all showed a statistically significant difference (P = 0.045 for PM#5 and 0.025 for PM#8).

Furthermore, one line with a deletion of four amino acids (FGDA) indicated as #121SRASL"FGDA"ERKKNSIKMW#138 (referred to as Del4aa) and another line with a deletion of six amino acids (SLFGDA) indicated as #121RK"SLFNDA"SRA"SLFGDA"ERKKNSIKMW#138 (referred to as Del6aa) were also obtained and undergone with the same irradiation experiments. Similarly, the deletion lines of Del4aa and Del6aa displayed a statistically significantly lower survival compared to littermate wild-type control mice (P = 0.0006 for Del4aa and 0.012 for Del4aa; Fig. 4c, d). This result provided strong evidence for a pivotal role of this short and conserved sequence to maintain BARD1 functions in mammals.

Our results added evidence that the inherited c.403G>A mutation in the highly conserved functional domain of BARD1 appears to suggest a favorable response of a triple-negative breast cancer patient to a PARP inhibitor, thus benefiting patients beyond carriers of BRCA1 or BRCA2 germline mutations. Incorporation of next-generation sequencing into clinical practice continues to expand the list of VUS in DNA damage repair genes, posing particular challenges in the clinical decision for carriers of the VUS about personalized drug therapy and genetic counseling. Our integrated approach by combining the identification of germline mutation with somatic and transcriptomic profiling of cancer patients, followed by functional assays at the molecular, cellular, and animal levels, will accelerate the interpretation of VUS for precision medicine.