SpadaHC: a database to improve the classification of variants in hereditary cancer genes in the Spanish population.

Accurate classification of genetic variants is crucial for clinical decision-making in hereditary cancer. In Spain, genetic diagnostic laboratories have traditionally approached this task independently due to the lack of a dedicated resource. Here we present SpadaHC, a web-based database for sharing variants in hereditary cancer genes in the Spanish population. SpadaHC is implemented using a three-tier architecture consisting of a relational database, a web tool and a bioinformatics pipeline. Contributing laboratories can share variant classifications and variants from individuals in Variant Calling Format (VCF) format. The platform supports open and restricted access, flexible dataset submissions, automatic pseudo-anonymization, VCF quality control, variant normalization and liftover between genome builds. Users can flexibly explore and search data, receive automatic discrepancy notifications and access SpadaHC population frequencies based on many criteria. In February 2024, SpadaHC included 18 laboratory members, storing 1.17 million variants from 4306 patients and 16 343 laboratory classifications. In the first analysis of the shared data, we identified 84 genetic variants with clinically relevant discrepancies in their classifications and addressed them through a three-phase resolution strategy. This work highlights the importance of data sharing to promote consistency in variant classifications among laboratories, so patients and family members can benefit from more accurate clinical management. Database URL: https://spadahc.ciberisciii.es/.

Open-Source Bioinformatic Pipeline to Improve PMS2 Genetic Testing Using Short-Read NGS Data.

The molecular diagnosis of mismatch repair-deficient cancer syndromes is hampered by difficulties in sequencing the PMS2 gene, mainly owing to the PMS2CL pseudogene. Next-generation sequencing short reads cannot be mapped unambiguously by standard pipelines, compromising variant calling accuracy. This study aimed to provide a refined bioinformatic pipeline for PMS2 mutational analysis and explore PMS2 germline pathogenic variant prevalence in an unselected hereditary cancer (HC) cohort. PMS2 mutational analysis was optimized using two cohorts: 192 unselected HC patients for assessing the allelic ratio of paralogous sequence variants, and 13 samples enriched with PMS2 (likely) pathogenic variants screened previously by long-range genomic DNA PCR amplification. Reads were forced to align with the PMS2 reference sequence, except those corresponding to exon 11, where only those intersecting gene-specific invariant positions were considered. Afterward, the refined pipeline's accuracy was validated in a cohort of 40 patients and used to screen 5619 HC patients. Compared with our routine diagnostic pipeline, the PMS2_vaR pipeline showed increased technical sensitivity (0.853 to 0.956, respectively) in the validation cohort, identifying all previously PMS2 pathogenic variants found by long-range genomic DNA PCR amplification. Fifteen HC cohort samples carried a pathogenic PMS2 variant (15 of 5619; 0.285%), doubling the estimated prevalence in the general population. The refined open-source approach improved PMS2 mutational analysis accuracy, allowing its inclusion in the routine next-generation sequencing pipeline streamlining PMS2 screening.

Tumor analysis of MMR genes in Lynch-like syndrome: Challenges associated with results interpretation.

BACKGROUND: Up to 70% of suspected Lynch syndrome patients harboring MMR deficient tumors lack identifiable germline pathogenic variants in MMR genes, being referred to as Lynch-like syndrome (LLS). Previous studies have reported biallelic somatic MMR inactivation in a variable range of LLS-associated tumors. Moreover, translating tumor testing results into patient management remains controversial. Our aim is to assess the challenges associated with the implementation of tumoral MMR gene testing in routine workflows. METHODS: Here, we present the clinical characterization of 229 LLS patients. MMR gene testing was performed in 39 available tumors, and results were analyzed using two variant allele frequency (VAF) thresholds (≥5% and ≥10%). RESULTS AND DISCUSSION: More biallelic somatic events were identified at VAF ≥ 5% than ≥10% (35.9% vs. 25.6%), although the rate of nonconcordant results regarding immunohistochemical pattern increased (30.8% vs. 20.5%). Interpretation difficulties question the current utility of the identification of MMR somatic hits in the diagnostic algorithm of suspected LS cases.

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.

Variants of uncertain significance (VUS) in cancer predisposing genes: What are we learning from multigene panels?

One of the main factors influencing the clinical utility of genetic tests for cancer predisposition is the ability to provide actionable classifications (ie pathogenic or benign). However, a large fraction of the variants identified in cancer predisposing genes (CPGs) are of uncertain significance (VUS), and cannot be used for clinical purposes either to identify individuals at risk or to drive treatment. Here we analyze the current status of VUS identification in a subset of 24 CPGs included by the American College of Medical Genetics/Association for Molecular Pathology in the list of genes that should be considered for the return of incidental findings. To this purpose we retrieved published literature using different search strings according to the frequency of the condition and we extracted corresponding data from ClinVar. The total number of VUS has not decreased with time, due to widespread multigene panel testing, and the relative yield of VUS compared to pathogenic variants is higher in more recent studies, which tend to involve series not selected for the presence of specific high risk criteria. In addition, only few studies adopt gene specific interpretation criteria when these are available. Despite the large yield of VUS associated with multigene testing, the data obtained from such studies can be very useful for variant classification, especially for those variants that are more likely to be benign, since these are expected to be detected more frequently in a population that does not show gene specific manifestations. In addition, wider use of gene specific interpretation criteria should be promoted in order to optimize the interpretation process.

Screening of CNVs using NGS data improves mutation detection yield and decreases costs in genetic testing for hereditary cancer.

INTRODUCTION: Germline CNVs are important contributors to hereditary cancer. In genetic diagnostics, multiplex ligation-dependent probe amplification (MLPA) is commonly used to identify them. However, MLPA is time-consuming and expensive if applied to many genes, hence many routine laboratories test only a subset of genes of interest. METHODS AND RESULTS: We evaluated a next-generation sequencing (NGS)-based CNV detection tool (DECoN) as first-tier screening to decrease costs and turnaround time and expand CNV analysis to all genes of clinical interest in our diagnostics routine. We used DECoN in a retrospective cohort of 1860 patients where a limited number of genes were previously analysed by MLPA, and in a prospective cohort of 2041 patients, without MLPA analysis. In the retrospective cohort, 6 new CNVs were identified and confirmed by MLPA. In the prospective cohort, 19 CNVs were identified and confirmed by MLPA, 8 of these would have been lost in our previous MLPA-restricted detection strategy. Also, the number of genes tested by MLPA across all samples decreased by 93.0% in the prospective cohort. CONCLUSION: Including an in silico germline NGS CNV detection tool improved our genetic diagnostics strategy in hereditary cancer, both increasing the number of CNVs detected and reducing turnaround time and costs.

Solving the enigma of POLD1 p.V295M as a potential cause of increased cancer risk.

Germline variants that affect the proofreading activity of polymerases epsilon (POLE) and delta (POLD1) predispose to colorectal adenomas and carcinomas, among other cancers. All cancer-associated pathogenic variants reported to date consist of non-disruptive genetic changes affecting the sequence that codifies the exonuclease domain (ED). Generally, disruptive (frameshift, stop-gain) POLE and POLD1 variants and missense variants outside the ED do not predispose to cancer. However, this statement may not be true for some, very specific variants that would indirectly affect the proofreading activity of the corresponding polymerase. We evaluated, by using multiple approaches, the possibility that POLD1 c.883G>A; p.(Val295Met), -a variant located 9 amino acids upstream the ED and present in ~0.25% of hereditary cancer patients-, affects POLD1 proofreading activity. Our findings show cumulative evidence that support no alteration of the proofreading activity and lack of association with cancer. The variant is classified as likely benign according to the ACMG/AMP guidelines.

Breast and Prostate Cancer Risks for Male BRCA1 and BRCA2 Pathogenic Variant Carriers Using Polygenic Risk Scores.

BACKGROUND: Recent population-based female breast cancer and prostate cancer polygenic risk scores (PRS) have been developed. We assessed the associations of these PRS with breast and prostate cancer risks for male BRCA1 and BRCA2 pathogenic variant carriers. METHODS: 483 BRCA1 and 1318 BRCA2 European ancestry male carriers were available from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). A 147-single nucleotide polymorphism (SNP) prostate cancer PRS (PRSPC) and a 313-SNP breast cancer PRS were evaluated. There were 3 versions of the breast cancer PRS, optimized to predict overall (PRSBC), estrogen receptor (ER)-negative (PRSER-), or ER-positive (PRSER+) breast cancer risk. RESULTS: PRSER+ yielded the strongest association with breast cancer risk. The odds ratios (ORs) per PRSER+ standard deviation estimates were 1.40 (95% confidence interval [CI] =1.07 to 1.83) for BRCA1 and 1.33 (95% CI = 1.16 to 1.52) for BRCA2 carriers. PRSPC was associated with prostate cancer risk for BRCA1 (OR = 1.73, 95% CI = 1.28 to 2.33) and BRCA2 (OR = 1.60, 95% CI = 1.34 to 1.91) carriers. The estimated breast cancer odds ratios were larger after adjusting for female relative breast cancer family history. By age 85 years, for BRCA2 carriers, the breast cancer risk varied from 7.7% to 18.4% and prostate cancer risk from 34.1% to 87.6% between the 5th and 95th percentiles of the PRS distributions. CONCLUSIONS: Population-based prostate and female breast cancer PRS are associated with a wide range of absolute breast and prostate cancer risks for male BRCA1 and BRCA2 carriers. These findings warrant further investigation aimed at providing personalized cancer risks for male carriers and informing clinical management.

RNA assay identifies a previous misclassification of BARD1 c.1977A>G variant.

Case-control studies have shown an association of BARD1 with hereditary breast and/or ovarian cancer (HBOC) predisposition. BARD1 alternatively spliced isoforms are abundant and some are highly expressed in different cancer types. In addition, a number of BARD1 germline pathogenic variants have been reported among HBOC patients. In previous reports, BARD1 c.1977A>G variant has been classified as pathogenic since it produces a frameshift transcript lacking exons 2 to 9. In the present study, we sought to validate the mRNA splicing results previously published and to contribute with new evidence to refine the classification of this substitution according to ACMG/AMP guidelines. The presence of the variant was screened in patients and controls. RT-PCR was performed in order to compare the transcriptional profiles of two variant carriers and ten non-carrier controls. In addition, allele-specific expression was assessed. No differences in variant frequency were detected between patients and controls. The RNA assay confirmed the presence of the shorter transcript lacking exons 2-9, but it was detected both in carriers and non-carriers. Furthermore, allelic imbalance was discarded and no significant differences in the proportion of full-length and shorter transcript were detected between carriers and controls. The shorter transcript detected corresponds to BARD1 isoform η, constituted by exons 1, 10 and 11. Our results support that this transcript is a constitutive splicing product rather than an aberrant transcript caused by BARD1 c.1977A>G variant, and for this reason this variant should be considered as likely benign following ACMG/AMP guidelines.

Genetic Screening for TLR7 Variants in Young and Previously Healthy Men With Severe COVID-19.

Introduction: Loss-of-function TLR7 variants have been recently reported in a small number of males to underlie strong predisposition to severe COVID-19. We aimed to determine the presence of these rare variants in young men with severe COVID-19. Methods: We prospectively studied males between 18 and 50 years-old without predisposing comorbidities that required at least high-flow nasal oxygen to treat COVID-19. The coding region of TLR7 was sequenced to assess the presence of potentially deleterious variants. Results: TLR7 missense variants were identified in two out of 14 patients (14.3%). Overall, the median age was 38 (IQR 30-45) years. Both variants were not previously reported in population control databases and were predicted to be damaging by in silico predictors. In a 30-year-old patient a maternally inherited variant [c.644A>G; p.(Asn215Ser)] was identified, co-segregating in his 27-year-old brother who also contracted severe COVID-19. A second variant [c.2797T>C; p.(Trp933Arg)] was found in a 28-year-old patient, co-segregating in his 24-year-old brother who developed mild COVID-19. Functional testing of this variant revealed decreased type I and II interferon responses in peripheral mononuclear blood cells upon stimulation with the TLR7 agonist imiquimod, confirming a loss-of-function effect. Conclusions: This study supports a rationale for the genetic screening for TLR7 variants in young men with severe COVID-19 in the absence of other relevant risk factors. A diagnosis of TLR7 deficiency could not only inform on treatment options for the patient, but also enables pre-symptomatic testing of at-risk male relatives with the possibility of instituting early preventive and therapeutic interventions.

Paired Somatic-Germline Testing of 15 Polyposis and Colorectal Cancer-Predisposing Genes Highlights the Role of APC Mosaicism in de Novo Familial Adenomatous Polyposis.

Familial adenomatous polyposis (FAP) is an autosomal dominant syndrome responsible for 1% of colorectal cancers (CRCs). Up to 90% of classic FAPs are caused by inactivating mutations in APC, and mosaicism has been previously reported in 20% of de novo cases, usually linked to milder phenotypic manifestations. This study aimed to explore the prevalence of mosaicism in 11 unsolved cases of classic FAP and to evaluate the diagnostic yield of somatic testing. Paired samples of colorectal polyps, tumors, and/or mucosa were analyzed using a custom next-generation sequencing panel targeting 15 polyposis and CRC-predisposing genes. Whenever possible, the extension of mosaicism to blood or sperm was also examined. Of 11 patients with classic adenomatous polyposis, a mosaic pathogenic variant in APC was identified in 7 (64%). No other altered genes were identified. In two of seven patients (29%), mosaicism was found restricted to colonic tissues, whereas in five of seven patients (71%), it was extended to the blood. Germline affectation was confirmed in one patient. We report the first analysis at a somatic level of 15 genes associated with CRC susceptibility, which highlights the role of APC mosaicism in classic FAP etiology. The results further reinforce the importance of testing target tissues when blood test results are negative.

CNVfilteR: an R/Bioconductor package to identify false positives produced by germline NGS CNV detection tools.

SUMMARY: Germline copy-number variants (CNVs) are relevant mutations for multiple genetics fields, such as the study of hereditary diseases. However, available benchmarks show that all next-generation sequencing (NGS) CNV calling tools produce false positives. We developed CNVfilteR, an R package that uses the single-nucleotide variant calls usually obtained in germline NGS pipelines to identify those false positives. The package can detect both false deletions and false duplications. We evaluated CNVfilteR performance on callsets generated by 13 CNV calling tools on three whole-genome sequencing and 541 panel samples, showing a decrease of up to 44.8% in false positives and consistent F1-score increase. Using CNVfilteR to detect false-positive calls can improve the overall performance of existing CNV calling pipelines. AVAILABILITY AND IMPLEMENTATION: CNVfilteR is released under Artistic-2.0 License. Source code and documentation are freely available at Bioconductor (http://www.bioconductor.org/packages/CNVfilteR). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Challenge of Diagnosing Constitutional Mismatch Repair Deficiency Syndrome in Brain Malignancies from Young Individuals.

Biallelic germline mismatch repair (MMR) gene (MLH1, MSH2, MSH6, and PMS2) mutations are an extremely rare event that causes constitutional mismatch repair deficiency (CMMRD) syndrome. CMMRD is underdiagnosed and often debuts with pediatric malignant brain tumors. A high degree of clinical awareness of the CMMRD phenotype is needed to identify new cases. Immunohistochemical (IHC) assessment of MMR protein expression and analysis of microsatellite instability (MSI) are the first tools with which to initiate the study of this syndrome in solid malignancies. MMR IHC shows a hallmark pattern with absence of staining in both neoplastic and non-neoplastic cells for the biallelic mutated gene. However, MSI often fails in brain malignancies. The aim of this report is to draw attention to the peculiar IHC profile that characterizes CMMRD syndrome and to review the difficulties in reaching an accurate diagnosis by describing the case of two siblings with biallelic MSH6 germline mutations and brain tumors. Given the difficulties involved in early diagnosis of CMMRD we propose the use of the IHC of MMR proteins in all malignant brain tumors diagnosed in individuals younger than 25 years-old to facilitate the diagnosis of CMMRD and to select those neoplasms that will benefit from immunotherapy treatment.

BARD1 Pathogenic Variants are Associated with Triple-Negative Breast Cancer in a Spanish Hereditary Breast and Ovarian Cancer Cohort.

Only a small fraction of hereditary breast and/or ovarian cancer (HBOC) cases are caused by germline variants in the high-penetrance breast cancer 1 and 2 genes (BRCA1 and BRCA2). BRCA1-associated ring domain 1 (BARD1), nuclear partner of BRCA1, has been suggested as a potential HBOC risk gene, although its prevalence and penetrance are variable according to populations and type of tumor. We aimed to investigate the prevalence of BARD1 truncating variants in a cohort of patients with clinical suspicion of HBOC. A comprehensive BARD1 screening by multigene panel analysis was performed in 4015 unrelated patients according to our regional guidelines for genetic testing in hereditary cancer. In addition, 51,202 Genome Aggregation Database (gnomAD) non-Finnish, non-cancer European individuals were used as a control population. In our patient cohort, we identified 19 patients with heterozygous BARD1 truncating variants (0.47%), whereas the frequency observed in the gnomAD controls was 0.12%. We found a statistically significant association of truncating BARD1 variants with overall risk (odds ratio (OR) = 3.78; CI = 2.10-6.48; p = 1.16 × 10-5). This association remained significant in the hereditary breast cancer (HBC) group (OR = 4.18; CI = 2.10-7.70; p = 5.45 × 10-5). Furthermore, deleterious BARD1 variants were enriched among triple-negative BC patients (OR = 5.40; CI = 1.77-18.15; p = 0.001) compared to other BC subtypes. Our results support the role of BARD1 as a moderate penetrance BC predisposing gene and highlight a stronger association with triple-negative tumors.

Characteristics of Adrenocortical Carcinoma Associated With Lynch Syndrome.

CONTEXT: Lynch syndrome (LS) is the most common inherited colorectal and endometrial cancer syndrome, caused by germline mutations in DNA mismatch repair (MMR) genes. It is also characterized by an increased risk of other tumors with lower prevalence, such as adrenal cortical carcinoma (ACC), an endocrine tumor with an incidence of <2 cases/million individuals/year. Most ACC developed during childhood are associated with hereditary syndromes. In adults, this association is not as well established as in children. Previous studies showed a 3.2% prevalence of LS among patients with ACC. EVIDENCE ACQUISITION: The objective of this study is to determine the prevalence of ACC in a Spanish LS cohort and their molecular and histological characteristics. This retrospective study includes 634 patients from 220 LS families registered between 1999 and 2018. EVIDENCE SYNTHESIS: During the follow-up period, 3 patients were diagnosed with ACC (0.47%); all were carriers of a MSH2 germline mutation. The 3 ACC patients presented loss of expression of MSH2 and MSH6 proteins. One tumor analysis showed loss of heterozygosity of the MSH2 wildtype allele. Our findings support previous data that considered ACC as a LS spectrum tumor. CONCLUSION: MMR protein immunohistochemistry screening could be an efficient strategy to detect LS in patients with ACC.

ERCC3, a new ovarian cancer susceptibility gene?

BACKGROUND: Hereditary breast and ovarian cancer syndrome (HBOC) is an inherited disorder with an increased risk of breast cancer (BC) and ovarian cancers (OC). Mutations in BRCA1-BRCA2 explains less than a half of cases. In the last decade several genes with different penetrance have been associated with an increased risk of BC or OC. A recurrent heterozygous ERCC3 truncating mutation increases the risk for breast cancer in patients with Ashkenazi Jewish ancestry. Our study aimed to investigate the role of ERCC3 truncating variants in a cohort of patients with suspicion of HBOC. PATIENTS AND METHODS: ERCC3 screening by multigene-panel analysis in 1311 unrelated patients after our regional consensus for genetic testing in hereditary cancer was done. In addition, 453 Spanish cancer-free individuals and 51,343 GnomAD non-Finnish, non-cancer European individuals were used as control populations. RESULTS: We identified 13 patients with heterozygous ERCC3 truncating variants (0.99%). Five of them also carried a mutation in a high- /moderate-penetrance HBOC gene (BRCA1, BRCA2, CHEK2, and TP53) being Multilocus Inherited Neoplasia Alleles syndrome (MINAS) patients. The frequency in 453 Spanish controls was of 0.22%; similar to that observed in 51,343 non-Finnish European GnomAD population (0.24%). We found an almost statistically significant association of truncating ERCC3 variants with BC (odds ratio [OR] = 2.25, confidence interval [CI] = 0.6-5.93, P = 0.11), and we observed for the first time a significant association with OC (OR = 4.74, CI = 1-14.34, P = 0.028), that holds even after removing MINAS cases. CONCLUSIONS: To our knowledge, this is the largest HBOC series comprehensively analysed for ERCC3 mutations, and the first study identifying ERCC3 as a cancer risk for OC.

Improving Genetic Testing in Hereditary Cancer by RNA Analysis: Tools to Prioritize Splicing Studies and Challenges in Applying American College of Medical Genetics and Genomics Guidelines.

RNA analyses are a potent tool to identify spliceogenic effects of DNA variants, although they are time-consuming and cannot always be performed. We present splicing assays of 20 variants that represent a variety of mutation types in 10 hereditary cancer genes and attempt to incorporate these results into American College of Medical Genetics and Genomics (ACMG) classification guidelines. Sixteen single-nucleotide variants, 3 exon duplications, and 1 single-exon deletion were selected and prioritized by in silico algorithms. RNA was extracted from short-term lymphocyte cultures to perform RT-PCR and Sanger sequencing, and allele-specific expression was assessed whenever possible. Aberrant transcripts were detected in 14 variants (70%). Variant interpretation was difficult, especially comparing old classification standards to generic ACMG guidelines and a proposal was devised to weigh functional analyses at RNA level. According to the ACMG guidelines, only 12 variants were reclassified as pathogenic/likely pathogenic because the other two variants did not gather enough evidence. This study highlights the importance of RNA studies to improve variant classification. However, it also indicates the challenge of incorporating these results into generic ACMG guidelines and the need to refine these criteria gene specifically. Nevertheless, 60% of variants were reclassified, thus improving genetic counseling and surveillance for carriers of these variants.

Comprehensive analysis and ACMG-based classification of CHEK2 variants in hereditary cancer patients.

CHEK2 variants are associated with intermediate breast cancer risk, among other cancers. We aimed to comprehensively describe CHEK2 variants in a Spanish hereditary cancer (HC) cohort and adjust the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) guidelines for their classification. First, three CHEK2 frequent variants were screened in a retrospective Hereditary Breast and Ovarian Cancer cohort of 516 patients. After, the whole CHEK2 coding region was analyzed by next-generation sequencing in 1848 prospective patients with HC suspicion. We refined ACMG-AMP criteria and applied different combined rules to classify CHEK2 variants and define risk alleles. We identified 10 CHEK2 null variants, 6 missense variants with discordant interpretation in ClinVar database, and 35 additional variants of unknown significance. Twelve variants were classified as (likely)-pathogenic; two can also be considered "established risk-alleles" and one as "likely risk-allele." The prevalence of (likely)-pathogenic variants in the HC cohort was 0.8% (1.3% in breast cancer patients and 1.0% in hereditary nonpolyposis colorectal cancer patients). Here, we provide ACMG adjustment guidelines to classify CHEK2 variants. We hope that this study would be useful for variant classification of other genes with low effect variants.