ATM germline pathogenic variants affect outcomes in children with ataxia-telangiectasia and hematological malignancies.

Ataxia-telangiectasia (A-T) is an autosomal-recessive disorder caused by pathogenic variants (PVs) of the ATM gene. Children with A-T are predisposed to hematological malignancies. We aimed to investigate their characteristics and outcomes in order to generate data-based treatment recommendations. In this multinational, observational study we report 202 patients aged ≤25 years with A-T and hematological malignancies from 25 countries. Ninety-one patients (45%) presented with mature B-cell lymphomas, 82 (41%) with acute lymphoblastic leukemia/lymphoma, 21(10%) with Hodgkin lymphoma and eight (4%) with other hematological malignancies. Four-year overall survival and event-free survival (EFS) were 50.8% (95% CI 43.6-59.1) and 47.9% (95% CI 40.8-56.2), respectively. Cure rates have not significantly improved over the last four decades (p=.76). The major cause of treatment failure was treatment-related mortality (TRM) with a four-year cumulative incidence of 25.9% (95% CI 19.5-32.4). Germline ATM PVs were categorized as null or hypomorphic and patients with available genetic data (n=110) were classified as having absent (n=81) or residual (n=29) ATM kinase activity. Four-year EFS was 39.4% (95% CI 29-53.3) vs 78.7% (95% CI 63.7-97.2), (p<.001), and TRM rates were 37.6% (95% CI 26.4-48.7) vs 4.0% (95% CI 0-11.8), (p=.017), for those with absent and residual ATM kinase activity, respectively. Absence of ATM kinase activity was independently associated with decreased EFS (HR=0.362, 95% CI 0.16-0.82; p=.009) and increased TRM (HR=14.11, 95% CI 1.36-146.31; p=.029). Patients with A-T and leukemia/lymphoma may benefit from de-escalated therapy for patients with absent ATM kinase activity and near-standard therapy regimens for those with residual kinase activity.

The BRCA2 R2645G variant increases DNA binding and induces hyper-recombination.

BRCA2 tumor suppressor protein ensures genome integrity by mediating DNA repair via homologous recombination (HR). This function is executed in part by its canonical DNA binding domain located at the C-terminus (BRCA2CTD), the only folded domain of the protein. Most germline pathogenic missense variants are located in this highly conserved region which binds to single-stranded DNA (ssDNA) and to the acidic protein DSS1. These interactions are essential for the HR function of BRCA2. Here, we report that the variant R2645G, identified in breast cancer and located at the DSS1 interface, unexpectedly increases the ssDNA binding activity of BRCA2CTDin vitro. Human cells expressing this variant display a hyper-recombination phenotype, chromosomal instability in the form of chromatid gaps when exposed to DNA damage, and increased PARP inhibitor sensitivity. In mouse embryonic stem cells (mES), this variant alters viability and confers sensitivity to cisplatin and Mitomycin C. These results suggest that BRCA2 interaction with ssDNA needs to be tightly regulated to limit HR and prevent chromosomal instability and we propose that this control mechanism involves DSS1. Given that several missense variants located within this region have been identified in breast cancer patients, these findings might have clinical implications for carriers.

Mosaic BRCA1 promoter methylation contribution in hereditary breast/ovarian cancer pedigrees.

PURPOSE: Mosaic BRCA1 promoter methylation (BRCA1meth) increases the risk of early-onset breast cancer, triple-negative breast cancer and ovarian cancer. As mosaic BRCA1meth are believed to occur de novo, their role in family breast/ovarian cancer has not been assessed. PATIENTS: Blood-derived DNA from 20 unrelated affected cases from families with aggregation of breast/ovarian cancer, but with no germline pathogenic variants in BRCA1/2, PALB2 or RAD51C/D, were screened by methylation-sensitive high-resolution melting. CpG analysis was performed by pyrosequencing on blood and buccal swab. Two probands carried a pathogenic variant in a moderate-penetrance gene (ATM and BARD1), and 8 of 18 others (44%) carried BRCA1meth (vs none of the 20 age-matched controls). Involvement of BRCA1 in tumourigenesis in methylated probands was demonstrated in most tested cases by detection of a loss of heterozygosity and a homologous recombination deficiency signature. Among the eight methylated probands, two had relatives with breast cancer with detectable BRCA1meth in blood, including one with high methylation levels in two non-tumour tissues. CONCLUSIONS: The high prevalence of mosaic BRCA1meth in patients with breast/ovarian cancer with affected relatives, as well as this first description of a family aggregation of mosaic BRCA1meth, shows how this de novo event can contribute to hereditary breast/ovarian cancer pedigrees.

Validation of the BOADICEA model in a prospective cohort of BRCA1/2 pathogenic variant carriers.

BACKGROUND: No validation has been conducted for the BOADICEA multifactorial breast cancer risk prediction model specifically in BRCA1/2 pathogenic variant (PV) carriers to date. Here, we evaluated the performance of BOADICEA in predicting 5-year breast cancer risks in a prospective cohort of BRCA1/2 PV carriers ascertained through clinical genetic centres. METHODS: We evaluated the model calibration and discriminatory ability in the prospective TRANsIBCCS cohort study comprising 1614 BRCA1 and 1365 BRCA2 PV carriers (209 incident cases). Study participants had lifestyle, reproductive, hormonal, anthropometric risk factor information, a polygenic risk score based on 313 SNPs and family history information. RESULTS: The full multifactorial model considering family history together with all other risk factors was well calibrated overall (E/O=1.07, 95% CI: 0.92 to 1.24) and in quintiles of predicted risk. Discrimination was maximised when all risk factors were considered (Harrell's C-index=0.70, 95% CI: 0.67 to 0.74; area under the curve=0.79, 95% CI: 0.76 to 0.82). The model performance was similar when evaluated separately in BRCA1 or BRCA2 PV carriers. The full model identified 5.8%, 12.9% and 24.0% of BRCA1/2 PV carriers with 5-year breast cancer risks of <1.65%, <3% and <5%, respectively, risk thresholds commonly used for different management and risk-reduction options. CONCLUSION: BOADICEA may be used to aid personalised cancer risk management and decision-making for BRCA1 and BRCA2 PV carriers. It is implemented in the free-access CanRisk tool (https://www.canrisk.org/).

Tumor mutational burden assessment and standardized bioinformatics approach using custom NGS panels in clinical routine.

BACKGROUND: High tumor mutational burden (TMB) was reported to predict the efficacy of immune checkpoint inhibitors (ICIs). Pembrolizumab, an anti-PD-1, received FDA-approval for the treatment of unresectable/metastatic tumors with high TMB as determined by the FoundationOne®CDx test. It remains to be determined how TMB can also be calculated using other tests. RESULTS: FFPE/frozen tumor samples from various origins were sequenced in the frame of the Institut Curie (IC) Molecular Tumor Board using an in-house next-generation sequencing (NGS) panel. A TMB calculation method was developed at IC (IC algorithm) and compared to the FoundationOne® (FO) algorithm. Using IC algorithm, an optimal 10% variant allele frequency (VAF) cut-off was established for TMB evaluation on FFPE samples, compared to 5% on frozen samples. The median TMB score for MSS/POLE WT tumors was 8.8 mut/Mb versus 45 mut/Mb for MSI/POLE-mutated tumors. When focusing on MSS/POLE WT tumor samples, the highest median TMB scores were observed in lymphoma, lung, endometrial, and cervical cancers. After biological manual curation of these cases, 21% of them could be reclassified as MSI/POLE tumors and considered as "true TMB high." Higher TMB values were obtained using FO algorithm on FFPE samples compared to IC algorithm (40 mut/Mb [10-3927] versus 8.2 mut/Mb [2.5-897], p < 0.001). CONCLUSIONS: We herein propose a TMB calculation method and a bioinformatics tool that is customizable to different NGS panels and sample types. We were not able to retrieve TMB values from FO algorithm using our own algorithm and NGS panel.

MSH3: a confirmed predisposing gene for adenomatous polyposis.

BACKGROUND: The MSH3 gene is part of the DNA mismatch repair system, but has never been shown to be involved in Lynch syndrome. A first report of four patients from two families, bearing biallelic MSH3 germline variants, with a phenotype of attenuated colorectal adenomatous polyposis raised the question of its involvement in hereditary cancer predisposition. The patients' tumours exhibited elevated microsatellite alterations at selected tetranucleotide repeats (EMAST), a hallmark of MSH3 deficiency. METHODS: We report five new unrelated patients with MSH3-associated polyposis. We describe their personal and familial history and study the EMAST phenotype in various normal and tumour samples, which are relevant findings based on the rarity of this polyposis subtype so far. RESULTS: All patients had attenuated colorectal adenomatous polyposis, with duodenal polyposis in two cases. Both women had breast carcinomas. EMAST phenotype was present at various levels in different samples of the five patients, confirming the MSH3 deficiency, with a gradient of instability in polyps depending on their degree of dysplasia. The negative EMAST phenotype ruled out the diagnosis of germline MSH3 deficiency for two patients: one homozygous for a benign variant and one with a monoallelic large deletion. CONCLUSION: This report lends further credence to biallelic MSH3 germline pathogenic variants being involved in colorectal and duodenal adenomatous polyposis. Large-scale studies may help clarify the tumour spectrum and associated risks. Ascertainment of EMAST may help with the interpretation of variants of unknown significance. We recommend adding MSH3 to dedicated diagnostic gene panels.

Germline HPF1 retrogene insertion in RB1 gene involved in cancer predisposition.

About half of the human genome is composed of repeated sequences derived from mobile elements, mainly retrotransposons, generally without pathogenic effect. Familial forms of retinoblastoma are caused by germline pathogenic variants in RB1 gene. Here, we describe a family with retinoblastoma affecting a father and his son. No pathogenic variant was identified after DNA analysis of RB1 gene coding sequence and exon-intron junctions. However, RB1 mRNA analysis showed a chimeric transcript with insertion of 114 nucleotides from HPF1 gene inside RB1 gene. This chimeric transcript led to an insertion of 38 amino acids in functional domain of retinoblastoma protein. Subsequent DNA analysis in RB1 intron 17 revealed the presence of a full-length HPF1 retrogene insertion in opposite orientation. Functional assay shows that this insertion has a deleterious impact on retinoblastoma protein function. This is the first report of a full-length retrogene insertion involved in human Mendelian disease leading to a chimeric transcript and a non-functional chimeric protein. Some retrogene insertions may be missed by standard diagnostic genetic testing, so contribution of retrogene insertions to human disease may be underestimated. The increasing use of whole genome sequencing in diagnostic settings will help to get a more comprehensive view of retrogenes.

APC germline pathogenic variants and epithelial ovarian cancer: causal or coincidental findings?

APC germline pathogenic variants result in predisposition to familial adenomatous polyposis and extraintestinal tumours such as desmoid fibromatosis, medulloblastomas and thyroid cancers. They have also been recently involved in ovarian microcystic stromal tumours. APC inactivation has been described at the tumour level in epithelial ovarian cancers (EOCs). Here, we report the identification of APC germline pathogenic variants in two patients diagnosed with premenopausal EOC in early 30s, with no other pathogenic variant detected in the known ovarian cancer predisposing genes. Subsequent tumour analysis showed neither a second hit of APC inactivation nor ß-catenin activation. Both tumours did not have a homologous recombination (HR) deficiency, pointing towards the implication of other genes than those involved in HR. APC may contribute to the carcinogenesis of EOC in a multifactorial context. Further studies are required to clarify the role of APC in predisposition to EOC.

Adaptive nanopore sequencing to determine pathogenicity of BRCA1 exonic duplication.

BRCA1 and BRCA2 are tumour suppressor genes that have been characterised as predisposition genes for the development of hereditary breast and ovarian cancers among other malignancies. The molecular diagnosis of this predisposition syndrome is based on the detection of inactivating variants of any type in those genes. But in the case of structural variants, functional consequences can be difficult to assess using standard molecular methods, as the precise resolution of their sequence is often impossible with short-read next generation sequencing techniques. It has been recently demonstrated that Oxford Nanopore long-read sequencing technology can accurately and rapidly provide genetic diagnoses of Mendelian diseases, including those linked to pathogenic structural variants. Here, we report the accurate resolution of a germline duplication event of exons 18-20 of BRCA1 using Nanopore sequencing with adaptive sampling target enrichment. This allowed us to classify this variant as pathogenic within a short timeframe of 10 days. This study provides a proof-of-concept that nanopore adaptive sampling is a highly efficient technique for the investigation of structural variants of tumour suppressor genes in a clinical context.

First report of medulloblastoma in a patient with MUTYH-associated polyposis.

AIMS: The mutY DNA glycosylase encoded by the MUTYH gene prevents G:C → T:A transversions through the base excision repair DNA repair system. Germline biallelic pathogenic variants in MUTYH cause an adenomatous polyposis called MUTYH-associated polyposis (MAP), an autosomal recessive disease (OMIM: 608456), with an increased risk of colorectal cancer. Digestive lesions in this context show an excess of G:C → T:A transversions, individualising a specific mutational signature associated with MUTYH deficiency called signature SBS36. Predisposition to other tumours in patients with germline biallelic pathogenic variants in MUTYH is suspected but remains unclear. We report the first case of medulloblastoma in a patient with MAP, carrying the homozygous pathogenic variant c.1227_1228dup, p.(Glu410Glyfs*43) in MUTYH. METHODS: Whole exome sequencing was performed on the medulloblastoma to enlighten single nucleotide variants of interest, microsatellite status and mutational signature. The objective was to determine the involvement of MUTYH deficiency in the oncogenesis of this medulloblastoma. RESULTS: The medulloblastoma has the mutational signature SBS36 and driver pathogenic variants in CTNNB1, PTCH1 and KDM6A corresponding to G:C → T:A transversions, suggesting a role of MUTYH deficiency in oncogenesis. CONCLUSIONS: Therefore, medulloblastoma could be a rare manifestation associated with germline biallelic pathogenic variants in MUTYH.

Hereditary cancer predispositions: Comparison of multigene panel sequencing on fresh-frozen breast/ovarian tumor versus blood.

In breast or ovarian cancer (BC/OC) patients with evocative personal and/or family history, multigene panel sequencing is performed on blood to diagnose hereditary predispositions. Additionally, BRCA1/BRCA2 testing can be performed on tumor sample for therapeutic purpose. The accuracy of multigene panel tumor analysis on BC/OC to detect predisposing germline pathogenic variants (gPV) has not been precisely assessed. By comparing sequencing data from blood and fresh-frozen tumor we show that tumor genomic instability causes pitfalls to consider when performing tumor testing to detect gPV. Even if loss of heterozygosity increases germline signal in most cases, somatic copy number variants (CNV) can mask germline CNV and collapse point gPV variant allele frequency (VAF). Moreover, VAF does not allow an accurate distinction between germline and somatic pathogenic variants.

Combined-modality treatment for locally advanced cervical cancer in a woman with Bloom-like syndrome: A case report and review of the literature.

We report the case of a 46-year-old woman with Bloom-like syndrome affected with locally advanced cervical cancer. She was treated with induction chemotherapy and radical radiation therapy concurrent with chemotherapy (carboplatin and paclitaxel). She was able to complete treatment, but grade III toxicities were observed. The limited relevant literature is presented. We conclude that the management of patients with DNA repair deficiency is challenging for the team in charge because of the potentially high sensitivity to treatment and the lack of clear recommendations in the literature. The main objective remains to deliver the optimal treatment while reducing toxicities.

Testicular Sertoli cell tumour and potentially testicular Leydig cell tumour are features of DICER1 syndrome.

DICER1 syndrome is a rare paediatric autosomal dominant inherited disorder predisposing to various benign and malignant tumours. It is caused by a germline pathogenic variant in DICER1, and the second hit for tumour development is usually a missense hotspot pathogenic variant in the DICER1 ribonuclease IIIb domain. While DICER1 predisposing variants account for about 60% of ovarian Sertoli-Leydig cell tumours, no DICER1-related testicular stromal tumours have been described. Here we report the first two cases of testicular stromal tumours in children carrying a DICER1 germline pathogenic variant: a case of Sertoli cell tumour and a case of Leydig cell tumour diagnosed at 2 and 12 years of age, respectively. A somatic DICER1 hotspot pathogenic variant was detected in the Sertoli cell tumour. This report extends the spectrum of DICER1-related tumours to include testicular Sertoli cell tumour and potentially testicular Leydig cell tumour. Diagnosis of a testicular Sertoli cell tumour should prompt DICER1 genetic testing so that patients with a DICER1 germline pathogenic variant can benefit from established surveillance guidelines. DICER1 genetic evaluation may be considered for testicular Leydig cell tumour. Our findings suggest that miRNA dysregulation underlies the aetiology of some testicular stromal tumours.

First estimates of diffuse gastric cancer risks for carriers of CTNNA1 germline pathogenic variants.

BACKGROUND: Pathogenic variants (PV) of CTNNA1 are found in families fulfilling criteria for hereditary diffuse gastric cancer (HDGC) but no risk estimates were available until now. The aim of this study is to evaluate diffuse gastric cancer (DGC) risks for carriers of germline CTNNA1 PV. METHODS: Data from published CTNNA1 families were updated and new families were identified through international collaborations. The cumulative risk of DGC by age for PV carriers was estimated with the genotype restricted likelihood (GRL) method, taking into account non-genotyped individuals and conditioning on all observed phenotypes and genotypes of the index case to obtain unbiased estimates. A non-parametric (NP) and the Weibull functions were used to model the shape of penetrance function with the GRL. Kaplan-Meier incidence curve and standardised incidence ratios were also computed. A 'leave-one-out' strategy was used to evaluate estimate uncertainty. RESULTS: Thirteen families with 46 carriers of PV were included. The cumulative risks of DGC at 80 years for carriers of CTNNA1 PV are 49% and 57%, respectively with the Weibull GRL and NP GRL methods. Risk ratios to population incidence reach particularly high values at early ages and decrease with age. At 40 years, they are equal to 65 and 833, respectively with the Weibull GRL and NP GRL. CONCLUSION: This is the largest series of CTNNA1 families that provides the first risk estimates of GC. These data will help to improve management and surveillance for these patients and support inclusion of CTNNA1 in germline testing panels.

Knowledge, acceptability and personal attitude toward pre-implantation 1 genetic testing (PGT) and pre-natal diagnosis (PND) for females carrying BRCA pathogenic variant according to fertility preservation experience.

PURPOSE: Preimplantation genetic testing (PGT-M) and prenatal diagnosis (PND) followed by medical termination of pregnancy when the fetus is affected are two procedures developed to avoid the transmission of a severe hereditary disease which can be proposed to females that carried BRCA pathogenic variants. These females can also be offered fertility preservation (FP) when diagnosed with cancer or even before a malignancy occurs. The aim of the study was to evaluate the acceptability and personal attitude of women carrying a BRCA mutation toward techniques that can prevent BRCA transmission to their progeny. METHODS: Female mutated for BRCA1 or BRCA2 were invited to complete an online survey of 49 queries anonymously between June and August 2022. RESULTS: A total of 87 participants responded to the online survey. Overall, 86.2% of women considered that PGT-M should be proposed to all BRCA mutation carriers regardless of the severity of the family history, and 47.1% considered or would consider PGT-M for themselves. For PND, these percentages were lower reaching 66.7% and 29.9%, respectively. Females with personal history of breast cancer or FP achievement were more prone to undergo PND for themselves despite the overall acceptability of this procedure. Among the subgroup who had undergone FP (n = 58), there was no significant difference in acceptance of principle and personal attitude toward PGT-M and PND compared to the group without FP. CONCLUSION: BRCA pathogenic variants female carriers do need information about reproductive issues, even if they are not prone to undergo PGT-M nor PND for themselves. TRIAL REGISTRATION NUMBER: N/A.

Novel germline MET pathogenic variants in French patients with papillary renal cell carcinomas type I.

Hereditary papillary renal cell carcinoma (HPRC) is a rare inherited renal cancer syndrome characterized by bilateral and multifocal papillary type 1 renal tumors (PRCC1). Activating germline pathogenic variants of the MET gene were identified in HPRC families. We reviewed the medical and molecular records of a large French series of 158 patients screened for MET oncogenic variants. MET pathogenic and likely pathogenic variants rate was 12.4% with 40.6% among patients with familial PRCC1 and 5% among patients with sporadic PRCC1. The phenotype in cases with MET pathogenic and likely pathogenic variants was characteristic: PRCC1 tumors were mainly bilateral (84.3%) and multifocal (87.5%). Histologically, six out of seven patients with MET pathogenic variant harbored biphasic squamoid alveolar PRCC. Genetic screening identified one novel pathogenic variant MET c.3389T>C, p.(Leu1130Ser) and three novel likely pathogenic variants: MET c.3257A>T, p.(His1086Leu); MET c.3305T>C, p.(Ile1102Thr) and MET c.3373T>G, p.(Cys1125Gly). Functional assay confirmed their oncogenic effect as they induced an abnormal focus formation. The genotype-phenotype correlation between MET pathogenic variants and PRCC1 presentation should encourage to widen the screening, especially toward nonfamilial PRCC1. This precise phenotype also constitutes a strong argument for the classification of novel missense variants within the tyrosine kinase domain when functional assays are not accessible.

Feasibility of personalized screening and prevention recommendations in the general population through breast cancer risk assessment: results from a dedicated risk clinic.

PURPOSE: A personalized approach to prevention and early detection based on known risk factors should contribute to early diagnosis and treatment of breast cancer. We initiated a risk assessment clinic for all women wishing to undergo an individual breast cancer risk assessment. METHODS: Women underwent a complete breast cancer assessment including a questionnaire, mammogram with evaluation of breast density, collection of saliva sample, consultation with a radiologist, and a breast cancer specialist. Women aged 40 or older, with 0 or 1 first-degree relative with breast cancer diagnosed after the age of 40 were eligible for risk assessment using MammoRisk, a machine learning-based tool that provides an individual 5-year estimated risk of developing breast cancer based on the patient's clinical data and breast density, with or without polygenic risk scores (PRSs). DNA was extracted from saliva samples for genotyping of 76 single-nucleotide polymorphisms. The individual risk was communicated to the patient, with individualized screening and prevention recommendations. RESULTS: A total of 290 women underwent breast cancer assessment, among which 196 women (68%) were eligible for risk assessment using MammoRisk (median age 52, range 40-72). When PRS was added to MammoRisk, 40% (n = 78) of patients were assigned a different risk category, with 28% (n = 55) of patients changing from intermediate to moderate or high risk. CONCLUSION: Individual risk assessment is feasible in the general population. Screening recommendations could be given based on individual risk. The use of PRS changed the risk score and screening recommendations in 40% of women.

Highly Sensitive Detection Method of DICER1 Tumor Hotspot Mutations by Drop-off Droplet Digital PCR.

BACKGROUND: DICER1 syndrome is an autosomal dominant inherited syndrome predisposing to various benign and malignant tumors, mainly occurring in children and young adults, requiring broad surveillance starting at birth with repeated irradiating imaging exams and sedations for young patients. It is caused by monoallelic germline pathogenic variants in the DICER1 gene. More than 90% of tumors bear an additional somatic DICER1 missense hotspot mutation, as a second hit, involving 1 of 6 codons clustered in exons 24 and 25. We designed and in vitro validated a drop-off droplet digital PCR (ddPCR) system to scan all DICER1 hotspot codons, allowing for a liquid biopsy test, an alternative to sedation and radiation exposure. METHODS: Three drop-off ddPCR assays were designed, with 2 TaqMan probes per assay, 1 complementary to the wild-type sequence of the region containing hotspots and another 1 used as a reference. Eight tumor-derived DNAs and 5 synthetic oligonucleotides bearing DICER1 hotspot mutations were tested. RESULTS: All tested mutations were detected, with a limit of detection ranging from 0.07% to 0.31% for codons p. E1705, p. D1709, and p. D1713 in exon 24 and from 0.06% to 0.15% for codons p. G1809, p. D1810, and p. E1813 in exon 25. CONCLUSIONS: The high sensitivity of this method is compatible with its use for plasma circulating tumor DNA (ctDNA) analysis for early tumor detection in DICER1 syndrome patients. It may reduce the need for radiation exposure and sedation in surveillance protocols and may also improve patient prognosis. Clinical trials are needed to evaluate ctDNA analysis in these patients.

Genetic clinicians' confidence in BOADICEA comprehensive breast cancer risk estimates and counselees' psychosocial outcomes: A prospective study.

Counseling for familial breast cancer focuses on communicating the gene test result (GENE) to counselees, but risk prediction models have become more complex by including non-genetic risk factors (NGRF) and polygenic risk scores (PRS). We examined genetic clinicians' confidence in counseling and counselees' psychosocial outcomes, using the BOADICEA risk prediction tool with different categories of risk factors as input. A prospective observational study in Dutch, French and German genetic clinics was performed including 22 clinicians, and 406 of 460 (88.3%) eligible cancer-unaffected women at high breast cancer risk assessed at pre-test and 350 (76.1%) at post-test. We performed multilevel analyses accounting for the clinician, and counselees' characteristics. Overall, risk estimates category by GENE versus GENE+ NGRF, or GENE+NGRF+PRS differed in 11% and 25% of counselees, respectively. In multilevel analyses, clinicians felt less confident in counseling when the full model provided lower breast cancer risks than GENE (i.e., in 8% of cases). Older counselees expressed higher breast cancer risk perception and worries about the hereditary predisposition when the full model provided higher breast cancer risks than GENE only. Genetic clinicians appear confident with breast cancer risk comprehensive models, which seem only to affect perceptions of older counselees.

Boosting GWAS using biological networks: A study on susceptibility to familial breast cancer.

Genome-wide association studies (GWAS) explore the genetic causes of complex diseases. However, classical approaches ignore the biological context of the genetic variants and genes under study. To address this shortcoming, one can use biological networks, which model functional relationships, to search for functionally related susceptibility loci. Many such network methods exist, each arising from different mathematical frameworks, pre-processing steps, and assumptions about the network properties of the susceptibility mechanism. Unsurprisingly, this results in disparate solutions. To explore how to exploit these heterogeneous approaches, we selected six network methods and applied them to GENESIS, a nationwide French study on familial breast cancer. First, we verified that network methods recovered more interpretable results than a standard GWAS. We addressed the heterogeneity of their solutions by studying their overlap, computing what we called the consensus. The key gene in this consensus solution was COPS5, a gene related to multiple cancer hallmarks. Another issue we observed was that network methods were unstable, selecting very different genes on different subsamples of GENESIS. Therefore, we proposed a stable consensus solution formed by the 68 genes most consistently selected across multiple subsamples. This solution was also enriched in genes known to be associated with breast cancer susceptibility (BLM, CASP8, CASP10, DNAJC1, FGFR2, MRPS30, and SLC4A7, P-value = 3 × 10-4). The most connected gene was CUL3, a regulator of several genes linked to cancer progression. Lastly, we evaluated the biases of each method and the impact of their parameters on the outcome. In general, network methods preferred highly connected genes, even after random rewirings that stripped the connections of any biological meaning. In conclusion, we present the advantages of network-guided GWAS, characterize their shortcomings, and provide strategies to address them. To compute the consensus networks, implementations of all six methods are available at https://github.com/hclimente/gwas-tools.