Search

OBJECTIVE: To determine the feasibility of universal genetic testing of women with newly diagnosed breast cancer, to estimate the incidence of pathogenic gene variants and their impact on patient management, and to evaluate patient and clinician acceptance of universal testing. DESIGN, SETTING, PARTICIPANTS: Prospective study of women with invasive or high grade in situ breast cancer and unknown germline status discussed at the Parkville Breast Service (Melbourne) multidisciplinary team meeting. Women were recruited to the pilot (12 June 2020 - 22 March 2021) and expansion phases (17 October 2021 - 8 November 2022) of the Mutational Assessment of newly diagnosed breast cancer using Germline and tumour genomICs (MAGIC) study. MAIN OUTCOME MEASURES: Germline testing by DNA sequencing, filtered for nineteen hereditary breast and ovarian cancer genes that could be classified as actionable; only pathogenic variants were reported. Surveys before and after genetic testing assessed pilot phase participants' perceptions of genetic testing, and psychological distress and cancer-specific worry. A separate survey assessed clinicians' views on universal testing. RESULTS: Pathogenic germline variants were identified in 31 of 474 expanded study phase participants (6.5%), including 28 of 429 women with invasive breast cancer (6.5%). Eighteen of the 31 did not meet current genetic testing eligibility guidelines (probability of a germline pathogenic variant ≥ 10%, based on CanRisk, or Manchester score ≥ 15). Clinical management was changed for 24 of 31 women after identification of a pathogenic variant. Including 68 further women who underwent genetic testing outside the study, 44 of 542 women carried pathogenic variants (8.1%). Acceptance of universal testing was high among both patients (90 of 103, 87%) and clinicians; no decision regret or adverse impact on psychological distress or cancer-specific worry were reported. CONCLUSION: Universal genetic testing following the diagnosis

BACKGROUND: PALB2 is the most important contributor to familial breast cancer after BRCA1 and BRCA2. Large genomic rearrangements (LGRs) in BRCA1 and BRCA2 are routinely assessed in clinical testing and are a significant contributor to the yield of actionable findings. In contrast, the contribution of LGRs in PALB2 has not been systematically studied. METHODS: We performed targeted sequencing and real-time qPCR validation to identify LGRs in PALB2 in 5770 unrelated patients with familial breast cancer and 5741 cancer-free control women from the same Australian population. RESULTS: Seven large deletions ranging in size from 0.96 kbp to 18.07 kbp involving PALB2 were identified in seven cases, while no LGRs were identified in any of the controls. Six LGRs were considered pathogenic as they included one or more exons of PALB2 and disrupted the WD40 domain at the C terminal end of the PALB2 protein while one LGR only involved a partial region of intron 10 and was considered a variant of unknown significance. Altogether, pathogenic LGRs identified in this study accounted for 10.3% (6 of 58) of the pathogenic PALB2 variants detected among the 5770 families with familial breast cancer. CONCLUSIONS: Our data show that a clinically important proportion of PALB2 pathogenic mutations in Australian patients with familial breast cancer are LGRs. Such observations have provided strong support for inclusion of PALB2 LGRs in routine clinical genetic testing.

BACKGROUND: Breast cancers (BCs) that arise in individuals heterozygous for a germline pathogenic variant in a susceptibility gene, such as BRCA1 and BRCA2, PALB2, and RAD51C, have been shown to exhibit biallelic loss in the respective genes and be associated with triple-negative breast cancer (TNBC) and distinctive somatic mutational signatures. Tumor sequencing thus presents an orthogonal approach to assess the role of candidate genes in BC development. METHODS: Exome sequencing was performed on paired normal-breast tumor DNA from 124 carriers of germline loss-of-function (LoF) or missense variant carriers in 15 known and candidate BC predisposition genes identified in the BEACCON case-control study. Biallelic inactivation and association with tumor genome features including mutational signatures and homologous recombination deficiency (HRD) score were investigated. RESULTS: BARD1-carrying TNBC (4 of 5) displayed biallelic loss and associated high HRD scores and mutational signature 3, as did a RAD51D-carrying TNBC and ovarian cancer. Biallelic loss was less frequent in BRIP1 BCs (4 of 13) and had low HRD scores. In contrast to other established BC genes, BCs from carriers of CHEK2 LoF (6 of 17) or missense (2 of 20) variant had low rates of biallelic loss. Exploratory analysis of BC from carriers of LoF variants in candidate genes such as BLM, FANCM, PARP2, and RAD50 found little evidence of biallelic inactivation. CONCLUSIONS: BARD1 and RAD51D behave as classic BRCA-like predisposition genes with biallelic inactivation, but this was not observed for any of the candidate genes. However, as demonstrated for CHEK2, the absence of biallelic inactivation does not provide definitive evidence against the gene's involvement in BC predisposition.

Clonal analysis of tumour sequencing data enables the evaluation of the relationship of histologically distinct synchronous lesions, such as co-existing benign areas, and temporally distinct tumours, such as primary-recurrence comparisons. In this review, we summarise statistical approaches that are commonly employed to define tumour clonal relatedness using data from bulk DNA technologies. We discuss approaches using total copy number, allele-specific copy number and mutation data, and the relative genomic resolution required for analysis and summarise some of the current tools for inferring clonal relationships. We argue that the impact of the biological context is critical in selecting any particular approach, such as the relative genomic complexity of the lesions being compared, and we recommend considering this context before employing any method to a new dataset.

While protein-truncating variants in RAD51C have been shown to predispose to triple-negative (TN) breast cancer (BC) and ovarian cancer, little is known about the pathogenicity of missense (MS) variants. The frequency of rare RAD51C MS variants was assessed in the BEACCON study of 5734 familial BC cases and 14,382 population controls, and findings were integrated with tumour sequencing data from 21 cases carrying a candidate variant. Collectively, a significant enrichment of rare MS variants was detected in cases (MAF < 0.001, OR 1.57, 95% CI 1.00-2.44, p = 0.05), particularly for variants with a REVEL score >0.5 (OR 3.95, 95% CI 1.40-12.01, p = 0.006). Sequencing of 21 tumours from 20 heterozygous and 1 homozygous carriers of nine candidate MS variants identified four cases with biallelic inactivation through loss of the wild-type allele, while six lost the variant allele and ten that remained heterozygous. Biallelic loss of the wild-type alleles corresponded strongly with ER- and TN breast tumours, high homologous recombination deficiency scores and mutational signature 3. Using this approach, the p.Gly264Ser variant, which was previously suspected to be pathogenic based on small case-control analyses and loss of activity in in vitro functional assays, was shown to be benign with similar prevalence in cases and controls and seven out of eight tumours showing no biallelic inactivation or characteristic mutational signature. Conversely, evaluation of case-control findings and tumour sequencing data identified p.Ile144Thr, p.Arg212His, p.Gln143Arg and p.Gly114Arg as variants warranting further investigation.

Pheochromocytomas (PC) and paragangliomas (PG) are rare neuroendocrine tumors associated with autonomic nerves. Here we use single-nuclei RNA-seq and bulk-tissue gene-expression data to characterize the cellular composition of PCPG and normal adrenal tissues, refine tumor gene-expression subtypes and make clinical and genotypic associations. We confirm seven PCPG gene-expression subtypes with significant genotype and clinical associations. Tumors with mutations in VHL, SDH-encoding genes (SDHx) or MAML3-fusions are characterized by hypoxia-inducible factor signaling and neoangiogenesis. PCPG have few infiltrating lymphocytes but abundant macrophages. While neoplastic cells transcriptionally resemble mature chromaffin cells, early chromaffin and neuroblast markers are also features of some PCPG subtypes. The gene-expression profile of metastatic SDHx-related PCPG indicates these tumors have elevated cellular proliferation and a lower number of non-neoplastic Schwann-cell-like cells, while GPR139 is a potential theranostic target. Our findings therefore clarify the diverse transcriptional programs and cellular composition of PCPG and identify biomarkers of potential clinical significance.

Gene expression by RNA polymerase II (RNAPII) is tightly controlled by cyclin-dependent kinases (CDKs) at discrete checkpoints during the transcription cycle. The pausing checkpoint following transcription initiation is primarily controlled by CDK9. We discovered that CDK9-mediated, RNAPII-driven transcription is functionally opposed by a protein phosphatase 2A (PP2A) complex that is recruited to transcription sites by the Integrator complex subunit INTS6. PP2A dynamically antagonizes phosphorylation of key CDK9 substrates including DSIF and RNAPII-CTD. Loss of INTS6 results in resistance to tumor cell death mediated by CDK9 inhibition, decreased turnover of CDK9 phospho-substrates, and amplification of acute oncogenic transcriptional responses. Pharmacological PP2A activation synergizes with CDK9 inhibition to kill both leukemic and solid tumor cells, providing therapeutic benefit in vivo. These data demonstrate that fine control of gene expression relies on the balance between kinase and phosphatase activity throughout the transcription cycle, a process dysregulated in cancer that can be exploited therapeutically.Copyright © 2021 Elsevier Inc.

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.Copyright © 2021, The Author(s).

Gene expression by RNA polymerase II (RNAPII) is tightly controlled by cyclin-dependent kinases (CDKs) at discrete checkpoints during the transcription cycle. The pausing checkpoint following transcription initiation is primarily controlled by CDK9. We discovered that CDK9-mediated, RNAPII-driven transcription is functionally opposed by a protein phosphatase 2A (PP2A) complex that is recruited to transcription sites by the Integrator complex subunit INTS6. PP2A dynamically antagonizes phosphorylation of key CDK9 substrates including DSIF and RNAPII-CTD. Loss of INTS6 results in resistance to tumor cell death mediated by CDK9 inhibition, decreased turnover of CDK9 phospho-substrates, and amplification of acute oncogenic transcriptional responses. Pharmacological PP2A activation synergizes with CDK9 inhibition to kill both leukemic and solid tumor cells, providing therapeutic benefit in vivo. These data demonstrate that fine control of gene expression relies on the balance between kinase and phosphatase activity throughout the transcription cycle, a process dysregulated in cancer that can be exploited therapeutically.Copyright © 2021 Elsevier Inc.

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.Copyright © 2021, The Author(s).

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.

Breast cancer (BC) has a significant heritable component but the genetic contribution remains unresolved in the majority of high-risk BC families. This study aims to investigate the monogenic causes underlying the familial aggregation of BC beyond BRCA1 and BRCA2, including the identification of new predisposing genes. A total of 11,511 non-BRCA familial BC cases and population-matched cancer-free female controls in the BEACCON study were investigated in two sequencing phases: 1303 candidate genes in up to 3892 cases and controls, followed by validation of 145 shortlisted genes in an additional 7619 subjects. The coding regions and exon-intron boundaries of all candidate genes and 14 previously proposed BC genes were sequenced using custom designed sequencing panels. Pedigree and pathology data were analysed to identify genotype-specific associations. The contribution of ATM, PALB2 and CHEK2 to BC predisposition was confirmed, but not RAD50 and NBN. An overall excess of loss-of-function (LoF) (OR 1.27, p = 9.05 × 10-9) and missense (OR 1.27, p = 3.96 × 10-73) variants was observed in the cases for the 145 candidate genes. Leading candidates harbored LoF variants with observed ORs of 2-4 and individually accounted for no more than 0.79% of the cases. New genes proposed by this study include NTHL1, WRN, PARP2, CTH and CDK9. The new candidate BC predisposition genes identified in BEACCON indicate that much of the remaining genetic causes of high-risk BC families are due to genes in which pathogenic variants are both very rare and convey only low to moderate risk.

Acute myeloid leukemia (AML) is a heterogeneous disease with poor prognosis and limited treatment strategies. Determining the role of cell-extrinsic regulators of leukemic cells is vital to gain clinical insights into the biology of AML. Iron is a key extrinsic regulator of cancer, but its systemic regulation remains poorly explored in AML. To address this question, we studied iron metabolism in patients with AML at diagnosis and explored the mechanisms involved using the syngeneic MLL-AF9-induced AML mouse model. We found that AML is a disorder with a unique iron profile, not associated with inflammation or transfusion, characterized by high ferritin, low transferrin, high transferrin saturation (TSAT), and high hepcidin. The increased TSAT in particular, contrasts with observations in other cancer types and in anemia of inflammation. Using the MLL-AF9 mouse model of AML, we demonstrated that the AML-induced loss of erythroblasts is responsible for iron redistribution and increased TSAT. We also show that AML progression is delayed in mouse models of systemic iron overload and that elevated TSAT at diagnosis is independently associated with increased overall survival in AML. We suggest that TSAT may be a relevant prognostic marker in AML.

Multimodal single-cell RNA sequencing enables the precise mapping of transcriptional and phenotypic features of cellular differentiation states but does not allow for simultaneous integration of critical posttranslational modification data. Here, we describe SUrface-protein Glycan And RNA-seq (SUGAR-seq), a method that enables detection and analysis of N-linked glycosylation, extracellular epitopes, and the transcriptome at the single-cell level. Integrated SUGAR-seq and glycoproteome analysis identified tumor-infiltrating T cells with unique surface glycan properties that report their epigenetic and functional state.

OBJECTIVE: Mucinous ovarian carcinoma (MOC) is an uncommon ovarian cancer histotype that responds poorly to conventional chemotherapy regimens. Although long overall survival outcomes can occur with early detection and optimal surgical resection, recurrent and advanced disease are associated with extremely poor survival. There are no current guidelines specifically for the systemic management of recurrent MOC. We analyzed data from a large cohort of women with MOC to evaluate the potential for clinical utility from a range of systemic agents. METHODS: We analyzed gene copy number (n = 191) and DNA sequencing data (n = 184) from primary MOC to evaluate signatures of mismatch repair deficiency and homologous recombination deficiency, and other genetic events. Immunohistochemistry data were collated for ER, CK7, CK20, CDX2, HER2, PAX8 and p16 (n = 117-166). RESULTS: Molecular aberrations noted in MOC that suggest a match with current targeted therapies include amplification of ERBB2 (26.7%) and BRAF mutation (9%). Observed genetic events that suggest potential efficacy for agents currently in clinical trials include: KRAS/NRAS mutations (66%), TP53 missense mutation (49%), RNF43 mutation (11%), ARID1A mutation (10%), and PIK3CA/PTEN mutation (9%). Therapies exploiting homologous recombination deficiency (HRD) may not be effective in MOC, as only 1/191 had a high HRD score. Mismatch repair deficiency was similarly rare (1/184). CONCLUSIONS: Although genetically diverse, MOC has several potential therapeutic targets. Importantly, the lack of response to platinum-based therapy observed clinically corresponds to the lack of a genomic signature associated with HRD, and MOC are thus also unlikely to respond to PARP inhibition.

High-grade serous ovarian carcinoma (HGSOC) has a significant hereditary component, approximately half of which cannot be explained by known genes. To discover genes, we analyse germline exome sequencing data from 516 BRCA1/2-negative women with HGSOC, focusing on genes enriched with rare, protein-coding loss-of-function (LoF) variants. Overall, there is a significant enrichment of rare protein-coding LoF variants in the cases (p < 0.0001, chi-squared test). Only thirty-four (6.6%) have a pathogenic variant in a known or proposed predisposition gene. Few genes have LoF mutations in more than four individuals and the majority are detected in one individual only. Forty-three highly-ranked genes are identified with three or more LoF variants that are enriched by three-fold or more compared to GnomAD. These genes represent diverse functional pathways with relatively few involved in DNA repair, suggesting that much of the remaining heritability is explained by previously under-explored genes and pathways.

Intraductal papillomas (IDP) are challenging breast findings because of their variable risk of progression to malignancy. The molecular events driving IDP development and genomic features of malignant progression are poorly understood. In this study, genome-wide CNA and/or targeted mutation analysis was performed on 44 cases of IDP, of which 20 cases had coexisting ductal carcinoma in situ (DCIS), papillary DCIS or invasive ductal carcinoma (IDC). CNA were rare in pure IDP, but 69% carried an activating PIK3CA mutation. Among the synchronous IDP cases, 55% (11/20) were clonally related to the synchronous DCIS and/or IDC, only one of which had papillary histology. In contrast to pure IDP, PIK3CA mutations were absent from clonal cases. CNAs in any of chromosomes 1, 16 or 11 were significantly enriched in clonal IDP lesions compared to pure and non-clonal IDP. The observation that 55% of IDP are clonal to DCIS/IDC indicates that IDP can be a direct precursor for breast carcinoma, not limited to the papillary type. The absence of PIK3CA mutations and presence of CNAs in IDP could be used clinically to identify patients at high risk of progression to carcinoma.

Mucinous ovarian carcinoma (MOC) is a unique subtype of ovarian cancer with an uncertain etiology, including whether it genuinely arises at the ovary or is metastatic disease from other organs. In addition, the molecular drivers of invasive progression, high-grade and metastatic disease are poorly defined. We perform genetic analysis of MOC across all histological grades, including benign and borderline mucinous ovarian tumors, and compare these to tumors from other potential extra-ovarian sites of origin. Here we show that MOC is distinct from tumors from other sites and supports a progressive model of evolution from borderline precursors to high-grade invasive MOC. Key drivers of progression identified are TP53 mutation and copy number aberrations, including a notable amplicon on 9p13. High copy number aberration burden is associated with worse prognosis in MOC. Our data conclusively demonstrate that MOC arise from benign and borderline precursors at the ovary and are not extra-ovarian metastases.Copyright © 2019, The Author(s).

Mucinous ovarian carcinoma (MOC) is a unique subtype of ovarian cancer with an uncertain etiology, including whether it genuinely arises at the ovary or is metastatic disease from other organs. In addition, the molecular drivers of invasive progression, high-grade and metastatic disease are poorly defined. We perform genetic analysis of MOC across all histological grades, including benign and borderline mucinous ovarian tumors, and compare these to tumors from other potential extra-ovarian sites of origin. Here we show that MOC is distinct from tumors from other sites and supports a progressive model of evolution from borderline precursors to high-grade invasive MOC. Key drivers of progression identified are TP53 mutation and copy number aberrations, including a notable amplicon on 9p13. High copy number aberration burden is associated with worse prognosis in MOC. Our data conclusively demonstrate that MOC arise from benign and borderline precursors at the ovary and are not extra-ovarian metastases.Copyright © 2019, The Author(s).

The BCL6 Corepressor (BCOR) is a component of a variant Polycomb repressive complex 1 (PRC1) that is essential for normal development. Recurrent mutations in the BCOR gene have been identified in acute myeloid leukaemia and myelodysplastic syndrome among other cancers; however, its function remains poorly understood. Here we examine the role of BCOR in haematopoiesis in vivo using a conditional mouse model that mimics the mutations observed in haematological malignancies. Inactivation of Bcor in haematopoietic stem cells (HSCs) results in expansion of myeloid progenitors and co-operates with oncogenic KrasG12D in the initiation of an aggressive and fully transplantable acute leukaemia. Gene expression analysis and chromatin immunoprecipitation sequencing reveals differential regulation of a subset of PRC1-target genes including HSC-associated transcription factors such as Hoxa7/9. This study provides mechanistic understanding of how BCOR regulates cell fate decisions and how loss of function contributes to the development of leukaemia.Copyright © 2019, The Author(s).

The BCL6 Corepressor (BCOR) is a component of a variant Polycomb repressive complex 1 (PRC1) that is essential for normal development. Recurrent mutations in the BCOR gene have been identified in acute myeloid leukaemia and myelodysplastic syndrome among other cancers; however, its function remains poorly understood. Here we examine the role of BCOR in haematopoiesis in vivo using a conditional mouse model that mimics the mutations observed in haematological malignancies. Inactivation of Bcor in haematopoietic stem cells (HSCs) results in expansion of myeloid progenitors and co-operates with oncogenic KrasG12D in the initiation of an aggressive and fully transplantable acute leukaemia. Gene expression analysis and chromatin immunoprecipitation sequencing reveals differential regulation of a subset of PRC1-target genes including HSC-associated transcription factors such as Hoxa7/9. This study provides mechanistic understanding of how BCOR regulates cell fate decisions and how loss of function contributes to the development of leukaemia.

Mucinous ovarian carcinoma (MOC) is a unique subtype of ovarian cancer with an uncertain etiology, including whether it genuinely arises at the ovary or is metastatic disease from other organs. In addition, the molecular drivers of invasive progression, high-grade and metastatic disease are poorly defined. We perform genetic analysis of MOC across all histological grades, including benign and borderline mucinous ovarian tumors, and compare these to tumors from other potential extra-ovarian sites of origin. Here we show that MOC is distinct from tumors from other sites and supports a progressive model of evolution from borderline precursors to high-grade invasive MOC. Key drivers of progression identified are TP53 mutation and copy number aberrations, including a notable amplicon on 9p13. High copy number aberration burden is associated with worse prognosis in MOC. Our data conclusively demonstrate that MOC arise from benign and borderline precursors at the ovary and are not extra-ovarian metastases.

BACKGROUND: Genome-wide association studies (GWASs) have identified numerous single-nucleotide polymorphisms (SNPs) associated with small increases in breast cancer risk. Studies to date suggest that some SNPs alter the expression of the associated genes, which potentially mediates risk modification. On this basis, we hypothesised that some of these genes may be enriched for rare coding variants associated with a higher breast cancer risk. METHODS: The coding regions and exon-intron boundaries of 56 genes that have either been proposed by GWASs to be the regulatory targets of the SNPs and/or located < 500 kb from the risk SNPs were sequenced in index cases from 1043 familial breast cancer families that previously had negative test results for BRCA1 and BRCA2 mutations and 944 population-matched cancer-free control participants from an Australian population. Rare (minor allele frequency ≤ 0.001 in the Exome Aggregation Consortium and Exome Variant Server databases) loss-of-function (LoF) and missense variants were studied. RESULTS: LoF variants were rare in both the cases and control participants across all the candidate genes, with only 38 different LoF variants observed in a total of 39 carriers. For the majority of genes (n = 36), no LoF variants were detected in either the case or control cohorts. No individual gene showed a significant excess of LoF or missense variants in the cases compared with control participants. Among all candidate genes as a group, the total number of carriers with LoF variants was higher in the cases than in the control participants (26 cases and 13 control participants), as was the total number of carriers with missense variants (406 versus 353), but neither reached statistical significance (p = 0.077 and p = 0.512, respectively). The genes contributing most of the excess of LoF variants in the cases included TET2, NRIP1, RAD51B and SNX32 (12 cases versus 2 control participants), whereas ZNF283 and CASP8 contributed largely to the exces