Introduction

Invasive lobular cancer (ILC) is the most common special histological subtype of breast cancer representing 10–15% of cases overall. ILC is associated with a higher risk of inheritance than invasive ductal cancers (IDC) which are also termed 'no special type' and represent 70% of cases overall [1, 2]. Pathogenic germline variants (PGVs) in one of the known high-risk breast cancer susceptibility genes e.g. BRCA1, BRCA2, PALB2 confer a high risk of developing breast cancer and are enriched in familial cases. Notably, lobular breast cancer (LBC) is less common as a proportion of breast cancers in women with BRCA1 PGVs compared to women without BRCA1 PGVs [3,4,5,6] As the main pathway to breast cancer is from basal progenitor cells leading to the predominant ductal triple negative breast cancer, this route excludes the predominantly estrogen receptor positive lobular cancer. However, individuals with BRCA2 PGVs have a similar proportion of LBC to women without BRCA2 PGVs. Familial breast cancer not associated with variants in BRCA1 or BRCA2 is more likely to be lobular than those with a BRCA1 or BRCA2 variant [7, 8]. In our first pathology update of the Manchester Scoring system, 11% of familial breast cancers that tested negative for BRCA1/2 were lobular, but lobular cancer was present in only 1.6% of index BRCA1 cases [8].

A recent increase in the use of multi-gene panels to screen for breast cancer-associated gene variants beyond BRCA1 and BRCA2 has started to identify PGVs in additional genes. Several studies have reported the association for the diffuse gastric cancer predisposition gene, CDH1, in lobular, rather than ductal breast cancer [9,10,11,12]. Petridis et al. screened six genes for an association with LBC and found that variants in BRCA2, CHEK2, and PALB2 in addition to CDH1 were all enriched in women with lobular compared to ductal cancer [5].

Only PGVs in the CDH1 gene have been convincingly shown to cause a specific high risk of LBC with lifetime risks of about 40% in women [13, 14]. However, CDH1 germline PGVs are very rare with a population incidence of 1 in 10–100,000 [15]. Indeed, a study of 1434 women with lobular breast cancer only found five (0.35%) with a CDH1 PGV [5]. As such there are still many individuals with LBC and a family history of breast cancer in whom no PGV has been identified.

Here, we present an analysis of the frequency of detectable pathogenic variants including those in an extended panel of nine non-BRCA1/2 genes in a cohort of 302 women with LBC.

Methods

Patient materials

Women were eligible for this study if they had a histologically confirmed diagnosis of LBC (lobular invasive lobular carcinoma). A total of 302 women affected with LBC were included in the study. Two-hundred and sixty four of these women were seen at the Manchester Centre for Genomic Medicine (MCGM) and 259 (98.1%) of these had a family history of breast or ovarian cancer or a second primary breast/ovarian cancer. The first families ascertained through MCGM were from 1990, although testing for BRCA1/2 only started in 1996. A further 38 women took part in the population-based study, Predicting the Risk Of Cancer At Screening (PROCAS), held in Greater Manchester [16]. Additionally, 1567 women without a breast cancer diagnosis at entry (aged 46–73 years) who were also recruited to the PROCAS study were included as controls. To obtain a population average dataset we included 124 women who subsequently developed breast cancer to provide an ~ 8% population risk of breast cancer at a median last age of follow up of 69 years. Clinical or research consent was given for extended testing of breast cancer associated genes (approval from the North Manchester Research Ethics Committee (reference 09/H1008/81(PROCAS) and 08/H1006/77)).

Genetic screening

For women that were seen at the MCGM, DNA was extracted from lymphocytes and all these samples were initially sequenced for PGVs in BRCA1/2 by a combination of Next Generation Sequencing and MLPA. A subset of 134 of the patients that did not harbour a PGV in either of those genes underwent extended testing. All women recruited through PROCAS provided a saliva sample for DNA extraction. They underwent panel testing as part of the Breast Cancer Risk after Diagnostic Gene Sequencing (BRIDGES) study [17]. The BRIDGES study performed sequencing of 33 genes (AKT1, ATM, BARD1, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, EPCAM, FAM175A, FANCC, FANCM, GEN1, MEN1, MRE11A, MSH2, MSH6, MUTYH, NBN, NF1, PALB2, PIK3CA, PMS2, PPM1D, PTEN, RAD50, RAD51C, RAD51D, RECQL, RINT1, STK11, TP53 and XRCC2). Forty-six women underwent clinical panel testing (the panel included: ATM, BAP1, BRCA1, BRCA2, CDH1, CHEK2, MLH1, NBN, PALB2, PTEN, RAD50, RAD51D, RECQL, TP53, and for 20 of these women BARD1, MSH6, MSH2 and PMS2 were also included) and 50 women underwent targeted exome sequencing through the Beijing Genomics Institute (BGI), Shenzhen, Guangdong, China. A minimum of nine genes were screened in each patient that underwent extended testing. These genes were: ATM, CDH1, CHEK2, NBN, PALB2, PTEN, RAD50, RAD51D, and TP53 (Supplemental Fig. 1).

Samples that underwent panel testing though BRIDGES did not undergo copy number analysis. The samples that underwent targeted exome sequencing only had copy number testing for BRCA1/2.

Variants were classified according to the ACMG guidelines [18]. All identified pathogenic variants were confirmed through the clinical diagnostics laboratory. Only variants that were classified as ‘likely pathogenic’ or ‘pathogenic’ are reported here.

Tumour pathology information was obtained for each case through hospital records, and cancer registries as previously described [19]. The probability of a BRCA1/2 PGV was determined using the Manchester score (MS) for each affected individual [20]. Genes were considered actionable for breast cancer if they had published data confirming at least a twofold relative risk of breast cancer [21]. Associations between LBC and PGVs in the screened genes were calculated by Fisher’s exact test using GraphPad Prism 8.

Results

Pathogenic variants

A total of 302 women affected with invasive LBC were initially screened for PGVs in BRCA1/2. Within this group 22 (7.28%) PGVs were identified in BRCA1/2 (5 in BRCA1 (1.66%) and 17 in BRCA2 (5.63%)). Of the 280 women who tested negative for BRCA1/2 PGVs, 134 (47.9%) had sufficient DNA available for extended panel testing of at least ATM, CDH1, CHEK2, NBN, PALB2, PTEN, RAD50, RAD51D, and TP53. This resulted in the detection of an additional 13 PGVs of which 10 (7.46%) could be considered to be actionable with at least a confirmed two fold increased breast cancer risk [21] (4 in ATM, 1 in CDH1, 1 in CHEK2, 2 in PALB2, and 2 in TP53) (Table 1). Additional PGVs genes with a less clear association with breast cancer were found in MSH6, RAD50, and NBN. Thus, in total, PGVs were detected in 35 individuals (11.59%). MSH6 was not assayed on every clinical gene panel and was therefore screened in 108 of the 134 cases. Twenty-five women with lobular carcinoma in situ were also screened for PGVs, but none of them carried a PGV.

Table 1 Pathogenic variants identified

In the control group of 1567 women, a total of 36 PGVs (2.30%) were detected. A first degree family history of breast cancer was present in 11.9% of PROCAS controls who provided a saliva DNA (and were analysed in BRIDGES) compared to 11.3% of PROCAS controls who did not provide a DNA sample. This indicates that there is no bias towards family history in providing DNA samples. There were eleven PGVs in BRCA1/2 (2 in BRCA1 and 9 in BRCA2), as well as six in ATM, six in CHEK2, three in PALB2, three in NBN, one in MSH6, and six in RAD50.

Odds ratios (ORs) for each gene are presented in Table 2. The ORs for BRCA2 (OR = 10.33 (95%CI 4.58–23.95; P < 0.0001)), BRCA1 (OR = 13.17 (95%CI 2.83–66.38; P = 0.0017)), and ATM (OR = 8.01 95%CI 2.52–29.92; P = 0.0053) confirmed elevated LBC risks with the lower 95% confidence intervals (CI) above twofold. We were unable to confirm this increase in LBC risk for PALB2 and CHEK2. CDH1 was non assessable due to absence of a reliable population incidence and the identification of only one case, previously reported [22]. All patients were also tested for CHEK2 c.1100delC. We only identified one woman with LBC and six controls with this PGV (OR = 0.86; 95%CI 0.08–5.23; P  ≥ 0.9999).

Table 2 Association of pathogenic variants in women with lobular breast cancer

ATM variant segregation

One of the individuals with an ATM PGV did not have a family history of breast cancer. For only one family, an additional family member was available to perform segregation analysis. This family member was affected with ductal breast cancer rather than lobular breast cancer and also carried the same ATM PGV. The other two families also had mixed lobular and ductal pathologies.

Receptor status

Of the 302 women, receptor status was known for 220 (72.82%). The majority were estrogen receptor (ER) positive (214 cases, 97.27%). Of the ER positive cases, HER2 receptor status was also known in 137 cases and, of these, 128 (93.43%) were HER2-negative (Table 3).

Table 3 Distribution of pathogenic variants according to Manchester score and receptor status

In three of the only six ER negative cases, a PGV was present (BRCA1(1); BRCA2(1); MSH6(1)). The majority of PGVs were identified in women with ER positive (14 PGVs: BRCA1(2); BRCA2(7); TP53(1); ATM(2); CDH1(1); CHEK2(1)) or ER positive/HER2- negative receptor status (11 PGVs: BRCA1(1); BRCA2(6); PALB2(1); ATM(2); RAD50(1)). One PGV was found in a woman with ER positive/HER2 positive tumour receptor status (BRCA2) and six PGVs (BRCA1(1); BRCA2(2); TP53(1); PALB2(1); NBN(1)) were found in women whose receptor status was unknown (Table 3).

Manchester score

In order to assess the probability of PGVs in BRCA1/2, MS was determined for all affected women. The majority of women had a MS < 20 (232/302) and a detection rate of 4.74%. The detection rate of BRCA1/2 PGVs in women with a MS > 20 was 15.71% (Table 3).

Only 43 of the 302 (14.23%) women did not have a family history of breast cancer and four of these carried an actionable PGV (1 in BRCA1, 1 in BRCA2 and 1 in TP53 and 1 in ATM). The BRCA1 and BRCA2 PGVs were identified in women diagnosed with a first primary lobular cancer who then developed a second primary breast or ovarian cancer, respectively, whereas the TP53 PGV was identified in a woman with a first primary lobular breast cancer at age 21 who then developed bilateral disease.

Discussion

In this cohort of 302 women diagnosed with LBC, we identified a PGV in 11.59%. Interestingly, the detection rate of PGVs beyond BRCA1/2 was 7.46%, which is similar to the BRCA1/2 detection rate (7.28%). Although the overall rate was low, the comparable rates of detection between BRCA1/2 and non-BRCA1/2 variants implies that extended testing may be particularly beneficial in women with LBC and a family history of breast cancer.

We identified four PGVs in ATM, which is equal to the number of PGVs detected in BRCA1 even though ATM was screened in less than half the individuals where BRCA1 was tested. The OR of 8.01 and two fold increased risk at the lower end of the confidence interval supports an association of ATM PGVs with LBC, which is consistent with Lu et al. [23] This group identified six PGVs in ATM in 369 (1.63%) patients with LBC (OR = 3.50; 95%CI 1.10–9.73), although this is a lower frequency than in our study (4/134 (2.99%)).

In line with previous observations, the four ATM PGVs identified in this study were found in women with ER + (/HER-) tumour characteristics [24, 25], which is also the most common tumour type seen in LBC [26].

In this study, we found as many PGVs in ATM as in PALB2 and CHEK2 combined. The recent large study by Petridis et al. [5] found that BRCA2, PALB2 and CHEK2 PGVs were the most prevalent in women affected with LBC [5]. However, ATM was not investigated in the Petridis study. A recent large American-based study of germline genetic testing criteria in 3907 women with breast cancer, identified 43 ATM variants, but did not distinguish between women with ductal and lobular cancers [27].

CDH1 is known to predispose to diffuse gastric cancer and has more recently been associated with LBC [5, 12]. In our cohort, we only identified one CDH1 PGV in a LBC family and there was no history of gastric cancer in this family. Another recent study reported an association with MSH6 (PGV identified in 7 of 590 patients) [23], although this was questioned as being due to potential sequencing errors [28]. We were not able to evaluate this fully as MSH6 variants were screened in only 80.6% of the samples tested and our study was not powered to refute this association.

We identified two cases with a TP53 PGV. This is in contrast to Petridis et al. who did not identify any TP53 PGVs in an unselected series of 1434 lobular cancers [5]. A lack of association with lobular cancer and germline TP53 was also suggested by Ditchi et al. in 2019 [6] but this is based on one case in only 57 carriers and the 95% CI do not exclude a 10% rate. One of our women was diagnosed with very early onset, bilateral breast cancer (the contralateral tumour was reported as invasive ductal carcinoma histologically). Very early onset ductal breast cancer is associated with TP53 PGVs and therefore potentially explains the TP53 PGV in this lobular case [29]. Given the extremely high odds ratios for TP53 in very early onset bilateral and familial cases an increased risk of lobular cancer is nonetheless still possible.

Interestingly, an MS of ≥ 20 is usually associated with a probability of BRCA1/2 involvement of around 35%; we only identified BRCA1/2 PGVs in 11/75 (14.67%) in women with a MS ≥ 20 despite the MS already including a downward adjustment for LBC [8]. This almost certainly indicates that even the current reduction in score for lobular breast cancer of -2 is insufficient to reflect BRCA1 risk. The absence of an association of LBC with BRCA1 has already been noted by ourselves and others [6], with only 2/342 (0.58%) LBCs in BRCA1 breast cancer patients in one study [6]. The higher odds ratio for BRCA1 in our study is therefore likely to be explained by the strong family history of breast and/or ovarian cancer in relatives rather than a true association with BRCA1. Although the MS is specifically used for BRCA1/2, it does contain an element of assessing high risk for inherited breast cancer by scoring all breast cancers in the lineage and increasing the score in younger cases. Inclusion of the remaining genes only increased this to 20% with an identified PGV. The relative absence of family histories of ovarian cancers in the series (only 43/302 and 8/134 of those with panel testing) would suggest that the absence of PGVs in so many with high MS is likely due to missing heritability. This indicates there is still a greater proportion of the heritability unexplained for lobular cancer.

Although we have confirmed associations with well described breast cancer genes and lobular cancer, most of the ORs in our study and the literature are in keeping with the ORs for overall breast cancer risk, which does not suggest that these genes have a specific link to lobular cancers. In particular the OR of only 0.86 for c.1100delC in CHEK2 is less than would be expected for ER + ductal cancer and does not support CHEK2 being particularly associated with LBC. However, the ORs for ATM of 8.01 fold in this study and 3.50 fold from a previous study [23] if confirmed in further work would imply a more specific link to lobular cancers. The ORs for BRCA1 in contrast were as high as expected for overall breast cancer risk; however this may reflect the rather low carrier frequency in the control population of only two of 1567 individuals, as well as the strong family histories of breast and/or ovarian cancer.

One of the limitations of this study was that we were unable to account for any copy number variations (CNVs) in samples tested through the BRIDGES study. Therefore, the results for some of the genes may be underestimated. For example, we have estimated that CNVs in BRCA1 account for 20% of the PGVs in non-Jewish families [30]. We also did not have detailed pathology of the breast cancers in the families of those with familial breast cancer, although a number of families did have confirmed lobular family history, including the family with the CDH1 variant. Also, as only a subset of the full cohort has been tested for the extended gene panel, there may be still undetected PGVs. In contrast to many other studies, we have been able to use local control samples rather than infer odds ratios from frequencies on gnomAD [31]. The relatively high frequency of BRCA1/2 in PROCAS controls of 0.56% does not appear to be due to bias as 1st degree breast cancer family history was not more frequent in subjects who provided a DNA sample. It is also very similar to a study of 50,726 adult biobank volunteers in the USA which found a 0.5% rate [32]. Another possible limitation is that there was no central pathology review. However, both of our TP53-positive cases were reviewed as part of their treatment at a major cancer centre (the Christie) and our study reflects actual practice, as pathologies are not usually reviewed in order to undertake genetic testing.

In summary, we have identified an association between ATM and an increased risk of LBC, but it is likely that further familial lobular cancer genes remain to be discovered as only 20% of patients with LBC and a family history of breast cancer with MS ≥ 20 were explained by currently known breast cancer predisposition genes.