Pitfalls in variant annotation for hereditary cancer diagnostics: The example of Illumina® VariantStudio®.

Next Generation Sequencing (NGS), and specifically targeted panel sequencing is the state-of-the-art in clinical genetic diagnosis of Mendelian diseases. However, the bioinformatics analysis and interpretation of the generated data can be challenging. A spotlight on the default transcript selection of a user-friendly, commercially available software that is widely used by genetics professionals, i.e. Illumina® VariantStudio®, is presented. For the sake of comparison, we employed Ensembl VEP, an open-source command-line tool, as it provides flexibility regarding transcript selection. The analysis of NGS data deriving from sequencing of 857 germline DNA samples of cancer patients indicated a concordance of 82.82% between the two software programs. Significantly, using the default transcript configuration of VariantStudio®, we failed to annotate correctly 11.45% of the identified loss-of-function variants. Our results underline the importance of cautious software and transcript selection and the need for reliable, white-box data analysis, along with bioinformatics expertise in clinical diagnostics.

CanVaS: Documenting the genetic variation spectrum of Greek cancer patients.

National genetic variation registries vastly increase the level of detail for the relevant population, while directly affecting patient management. Herein, we report CanVaS, a Cancer Variation reSource aiming to document the genetic variation of cancer patients in Greece. CanVaS comprises germline genetic data from 7,363 Greek individuals with a personal and/or family history of malignancy. The data set incorporates approximately 24,000 functionally annotated rare variants in 97 established or suspected cancer susceptibility genes. For each variant, allele frequency for the Greek population, interpretation for clinical significance, anonymized family and segregation information, as well as phenotypic traits of the carriers, are included. Moreover, information on the geographic distribution of the variants across the country is provided, enabling the study of Greek population isolates. Direct comparisons between Greek (sub)populations with relevant genetic resources are supported, allowing fine-grain localized adjustment of guidelines and clinical decision-making. Most importantly, anonymized data are available for download, while the Leiden Open Variation Database schema is adopted, enabling integration/interconnection with central resources. CanVaS could become a stepping-stone for a countrywide effort to characterize the cancer genetic variation landscape, concurrently supporting national and international cancer research. The database can be accessed at: http://ithaka.rrp.demokritos.gr/CanVaS.

One in three highly selected Greek patients with breast cancer carries a loss-of-function variant in a cancer susceptibility gene.

BACKGROUND: Gene panel testing has become the norm for assessing breast cancer (BC) susceptibility, but actual cancer risks conferred by genes included in panels are not established. Contrarily, deciphering the missing hereditability on BC, through identification of novel candidates, remains a challenge. We aimed to investigate the mutation prevalence and spectra in a highly selected cohort of Greek patients with BC, questioning an extensive number of genes, implicated in cancer predisposition and DNA repair, while calculating gene-specific BC risks that can ultimately lead to important associations. METHODS: To further discern BC susceptibility, a comprehensive 94-cancer gene panel was implemented in a cohort of 1382 Greek patients with BC, highly selected for strong family history and/or very young age (<35 years) at diagnosis, followed by BC risk calculation, based on a case-control analysis. RESULTS: Herein, 31.5% of patients tested carried pathogenic variants (PVs) in 28 known, suspected or candidate BC predisposition genes. In total, 24.8% of the patients carried BRCA1/2 loss-of-function variants. An additional 6.7% carried PVs in additional genes, the vast majority of which can be offered meaningful clinical changes. Significant association to BC predisposition was observed for ATM, PALB2, TP53, RAD51C and CHEK2 PVs. Primarily, compared with controls, RAD51C PVs and CHEK2 damaging missense variants were associated with high (ORs 6.19 (Exome Aggregation Consortium (ExAC)) and 12.6 (Fabulous Ladies Over Seventy (FLOSSIES)), p<0.01) and moderate BC risk (ORs 3.79 (ExAC) and 5.9 (FLOSSIES), p<0.01), respectively. CONCLUSION: Studying a large and unique cohort of highly selected patients with BC, deriving from a population with founder effects, provides important insight on distinct associations, pivotal for patient management.

BRCA1 and BRCA2 germline testing in Cretan isolates reveals novel and strong founder effects.

Germline BRCA1 and BRCA2 loss-of-function variants have been linked to increased breast and ovarian cancer risk, with more than 5,000 distinct pathogenic variants being reported worldwide. Among individuals of Greek descent, the BRCA1/2 variant spectrum is heterogeneous, but characterized by strong founder effects. As patients from certain geographical regions of Greece (like Crete) were underrepresented in previous studies, we hypothesized that isolated Cretans, a southern Greece islanders' population with distinct demographic, cultural and genetic features, could harbor founder BRCA1/2 mutations. A total of 304 breast or/and ovarian cancer patients of Cretan descent, fulfilling NCCN criteria for genetic testing, were tested by NGS or Sanger sequencing, followed by MLPA. Haplotype analysis was subsequently performed to investigate potential founder effects of recurrent alleles. Overall, 16.5% (50/304) of the tested patients carried 22 different pathogenic variants; 48% in BRCA1, 52% in BRCA2. Three variants, namely two in BRCA2 (Δexons 12 and 13 and c.7806-2A>T) and one in BRCA1 (c.5492del), constituting approximately half (48%) of all detected pathogenic variants, were shown to have a founder effect, with all carriers sharing common haplotypes. Remarkably, these variants were confined to Cretans and have not been identified in other regions of Greece. The high prevalence of specific BRCA1/2 pathogenic variants among Cretans, provides the possibility of cost- and time-efficient screening of the Cretan population. Integrating this knowledge in local public health services may have a significant impact on cancer prevention, and may serve as a starting point for the implementation of testing on a population level.

Functional characterization of CHEK2 variants in a Saccharomyces cerevisiae system.

Genetic testing for cancer predisposition leads to the identification of a number of variants with uncertain significance. To some extent, variants of BRCA1/2 have been classified, in contrast to variants of other genes. CHEK2 is a typical example, in which a large number of variants of unknown clinical significance were identified and still remained unclassified. Herein, the CHEK2 variant assessment was performed through an in vivo, yeast-based, functional assay. In total, 120 germline CHEK2 missense variants, distributed along the protein sequence, and two large in-frame deletions were tested, originating from genetic test results in breast cancer families, or selected from the ClinVar database. Of these, 32 missense and two in-frame deletions behaved as non-functional, 73 as functional, and 15 as semi-functional, after comparing growth rates of each strain with positive and negative controls. The majority of non-functional variants were localized in the CHK2 kinase and forkhead-associated domains. In vivo results from the non-functional variants were in agreement with in silico predictions, and, where available, with strong breast cancer family history, to a great extent. The results of the largest, to date, yeast-based assay, evaluating CHEK2 variants, can complement and assist in the classification of rare CHEK2 variants with unclear clinical significance.

PALB2 c.2257C>T truncating variant is a Greek founder and is associated with high breast cancer risk.

PALB2 loss-of-function variants play an important role in breast, pancreatic and possibly, ovarian and gastric cancer susceptibility. Their frequency can be influenced by founder effects, already described in some populations. Herein, we have assessed the possible founder effect of PALB2 c.2257C>T (p.Arg753*) truncating variant among Greek breast cancer patients, while investigating possible correlations with cancer diagnoses. Following a lead deriving from a background study of highly selected Greek breast cancer patients, a total of 2496 breast and 697 ovarian cancer patients were directly genotyped for the PALB2 c.2257C>T truncating variant. Consequently, haplotype analysis was conducted on identified carriers, using seven microsatellite markers. The prevalence of the PALB2 variant was 0.24% (6/2496) and 0.14% (1/697) among breast and ovarian cases, respectively. Family history seems to be an important factor for the variant identification, although not reaching statistical significance. Microsatellite analysis on 12 carriers revealed a common shared haplotype, spanning a chromosomal region of ~1.2 Mb; the variant was possibly introduced in the Greek population ~1600 years ago. The variant confers high breast cancer risk, as illustrated by comparison with publicly available control groups. Genetic testing for PALB2, especially for the Greek founder c.2257C>T truncating variant, should be seriously considered in Greek breast cancer cases, since such findings could assist appropriate clinical management for the patients and their families.

Germline deleterious mutations in genes other than BRCA2 are infrequent in male breast cancer.

PURPOSE: Male breast cancer (MBC) is a rare cancer entity, with mutations in BRCA1 and BRCA2 genes accounting for ~ 10% of patients. Multiple-gene sequencing has already entered clinical practice for female breast cancer, whereas the performance of panel testing in MBC has not been studied extensively. Therefore, the aim of this study was to evaluate the clinical utility of panel testing for MBC, by the largest gene panel used so far, through investigation of patients deriving from a population with known founder effects. METHODS: Genomic DNA from one hundred and two Greek MBC patients, unselected for age and family history, was used to prepare libraries which capture the entire coding regions of 94 cancer genes. RESULTS: Loss-of-function (LoF) mutations were found in 12.7% of the cases, distributed in six genes: BRCA2, ATM, BRCA1, CHEK2, PMS2, and FANCL. BRCA2 mutations were the most frequent, followed by ATM mutations, accounting for 6.9 and 2%, respectively, while mutations in other genes were detected in single cases. Age at diagnosis or family history was not predictive of mutation status. Beyond mutations in established breast cancer predisposing genes, LoF mutations in PMS2 and FANCL among MBC patients are reported here for the first time. CONCLUSIONS: Our findings, using the largest gene panel for MBC patients so far, indicate that BRCA testing should be the primary concern for MBC patients. Until sufficient evidence arises from larger studies, multiple-gene panels may be of limited benefit for MBC and their families, at least for MBC patients of specific descent.

Contribution of RAD51D germline mutations in breast and ovarian cancer in Greece.

RAD51D gene's protein product is known to be involved in the DNA repair mechanism by homologous recombination. RAD51D germline mutations have been recently associated with ovarian and breast cancer (OC and BC, respectively) predisposition. Our aim was to evaluate the frequency of hereditary RAD51D mutations in Greek patients. To address this, we have screened for RAD51D germline mutations 609 BRCA1- and BRCA2-negative patients diagnosed with OC, unselected for age or family history, and 569 BC patients diagnosed under 55 years and with an additional relative with BC or OC. We identified four pathogenic mutations in four unrelated individuals with family history of BC and/or OC. Three of the RAD51D carriers had developed BC, while the other one was an OC patient, thus accounting for a mutation frequency of 0.16% in the OC cohort and 0.53% in the BC cohort. One of the detected mutations is novel (c.738 + 1G > A), whereas the rest had been detected previously (p.Gln151Ter, p.Arg186Ter, and p.Arg300Ter). It is noteworthy that the 4 carrier families had 13 BC cases and only 4 OC cases. Our data support that RAD51D should be implemented into the comprehensive multigene panel, as mutation carriers may benefit from the administration of PARP inhibitors.

Genetic Testing of Breast Cancer Patients with Very Early-Onset Breast Cancer (≤30 Years) Yields a High Rate of Germline Pathogenic Variants, Mainly in the BRCA1, TP53, and BRCA2 Genes.

Early-onset breast cancer constitutes a major criterion for genetic testing referral. Nevertheless, studies focusing on breast cancer patients (≤30 years) are limited. We investigated the contribution and spectrum of known breast-cancer-associated genes in 267 Greek women with breast cancer ≤30 years while monitoring their clinicopathological characteristics and outcomes. In this cohort, a significant proportion (39.7%) carried germline pathogenic variants (PVs) distributed in 8 genes. The majority, namely 36.7%, involved BRCA1, TP53, and BRCA2. PVs in BRCA1 were the most prevalent (28.1%), followed by TP53 (4.5%) and BRCA2 (4.1%) PVs. The contribution of PVs in CHEK2, ATM, PALB2, PTEN, and RAD51C was limited to 3%. In the patient group ≤26 years, TP53 PVs were significantly higher compared to the group 26-30 years (p = 0.0023). A total of 74.8% of TP53 carriers did not report a family history of cancer. Carriers of PVs receiving neoadjuvant chemotherapy showed an improved event-free survival (p < 0.0001) compared to non-carriers. Overall, many women with early-onset breast cancer carry clinically actionable variants, mainly in the BRCA1/2 and TP53 genes. The inclusion of timely testing of TP53 in these patients provides essential information for appropriate clinical management. This is important for countries where reimbursement involves the cost of genetic analysis of BRCA1/2 only.

POLQ inhibition elicits an immune response in homologous recombination-deficient pancreatic adenocarcinoma via cGAS/STING signaling.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy that harbors mutations in homologous recombination-repair (HR-repair) proteins in 20%-25% of cases. Defects in HR impart a specific vulnerability to poly ADP ribose polymerase inhibitors and platinum-containing chemotherapy in tumor cells. However, not all patients who receive these therapies respond, and many who initially respond ultimately develop resistance. Inactivation of the HR pathway is associated with the overexpression of polymerase theta (Polθ, or POLQ). This key enzyme regulates the microhomology-mediated end-joining (MMEJ) pathway of double-strand break (DSB) repair. Using human and murine HR-deficient PDAC models, we found that POLQ knockdown is synthetically lethal in combination with mutations in HR genes such as BRCA1 and BRCA2 and the DNA damage repair gene ATM. Further, POLQ knockdown enhances cytosolic micronuclei formation and activates signaling of cyclic GMP-AMP synthase-stimulator of interferon genes (cGAS-STING), leading to enhanced infiltration of activated CD8+ T cells in BRCA2-deficient PDAC tumors in vivo. Overall, POLQ, a key mediator in the MMEJ pathway, is critical for DSB repair in BRCA2-deficient PDAC. Its inhibition represents a synthetic lethal approach to blocking tumor growth while concurrently activating the cGAS-STING signaling pathway to enhance tumor immune infiltration, highlighting what we believe to be a new role for POLQ in the tumor immune environment.

Single-cell RNA sequencing reveals the effects of chemotherapy on human pancreatic adenocarcinoma and its tumor microenvironment.

The tumor microenvironment (TME) in pancreatic ductal adenocarcinoma (PDAC) is a complex ecosystem that drives tumor progression; however, in-depth single cell characterization of the PDAC TME and its role in response to therapy is lacking. Here, we perform single-cell RNA sequencing on freshly collected human PDAC samples either before or after chemotherapy. Overall, we find a heterogeneous mixture of basal and classical cancer cell subtypes, along with distinct cancer-associated fibroblast and macrophage subpopulations. Strikingly, classical and basal-like cancer cells exhibit similar transcriptional responses to chemotherapy and do not demonstrate a shift towards a basal-like transcriptional program among treated samples. We observe decreased ligand-receptor interactions in treated samples, particularly between TIGIT on CD8 + T cells and its receptor on cancer cells, and identify TIGIT as the major inhibitory checkpoint molecule of CD8 + T cells. Our results suggest that chemotherapy profoundly impacts the PDAC TME and may promote resistance to immunotherapy.

MARGINAL: An Automatic Classification of Variants in BRCA1 and BRCA2 Genes Using a Machine Learning Model.

Implementation of next-generation sequencing (NGS) for the genetic analysis of hereditary diseases has resulted in a vast number of genetic variants identified daily, leading to inadequate variant interpretation and, consequently, a lack of useful clinical information for treatment decisions. Herein, we present MARGINAL 1.0.0, a machine learning (ML)-based software for the interpretation of rare BRCA1 and BRCA2 germline variants. MARGINAL software classifies variants into three categories, namely, (likely) pathogenic, of uncertain significance and (likely) benign, implementing the criteria established by the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP). We first annotated BRCA1 and BRCA2 variants using various sources. Then, we automatically implemented the ACMG-AMP criteria, and we finally constructed the ML model for variant classification. To maximize accuracy, we compared the performance of eight different ML algorithms in a classification scheme based on a serial combination of two classifiers. The model showed high predictive abilities with maximum accuracy of 92% and 98%, recall of 92% and 98% and specificity of 90% and 98% for the first and second classifiers, respectively. Our results indicate that using a gene and disease-specific ML automated software for clinical variant evaluation can minimize conflicting interpretations.

CHEK2 Pathogenic Variants in Greek Breast Cancer Patients: Evidence for Strong Associations with Estrogen Receptor Positivity, Overuse of Risk-Reducing Procedures and Population Founder Effects.

CHEK2 germline pathogenic variants predispose to breast cancer and possibly to other malignancies, with their spectrum and frequency being variable among populations. Τhe majority of CHEK2-associated breast tumors are hormone receptor positive; however, relevant clinical outcomes are not well defined. Herein, we illustrate the histopathological characteristics and clinical outcomes of 52 Greek breast cancer patients who are CHEK2 carriers. Genetic analysis was performed by Sanger/massively parallel sequencing, followed by MLPA. Subsequent haplotype analysis investigated possible founder effects. Blood relatives were offered cascade testing. CHEK2 variant spectrum was characterized by variability, while influenced by founder effects. The majority of carriers, i.e., 60.8%, were diagnosed with breast cancer before the age of 45. Notably, 91.5% of breast tumors were hormone receptor positive. Hormone therapy and mastectomy at diagnosis seem to have a positive trend on overall survival, after a median follow-up of 9.5 years. Remarkably, 41.9% of patients underwent risk-reducing surgery, one third of which involved salpingo-oophorectomy. Nearly half of families responded to cascade testing. Our data highlight the need for guideline-adherent choices, based on the evidence that CHEK2 carriers are at moderate risk for breast cancer and no risk for ovarian cancer, while underscore the possible role of chemoprevention with tamoxifen.

Genetic analysis and clinical description of Greek patients with Peutz-Jeghers syndrome: Creation of a National Registry.

Peutz-Jeghers syndrome (PJS) is a rare autosomal dominant disorder caused by germline mutations in the STK11 tumor suppressor gene. PJS patients face a cumulative cancer risk as high as 93% for all sites combined. The present study reports the spectrum of STK11 mutations in eight families with clinical diagnosis of PJS, summarizes the clinical characteristics of sixteen mutation carriers and launches a National Registry for PJS in Greece. STK11 loss-of-function (LoF) mutations were detected in 87.5% of index patients. Carriers presented with their first manifestation at a median age of 24.9 years, while early-onset breast cancer was the most frequent malignancy observed, highlighting the need for breast surveillance. Out of the deleterious STK11 mutations identified, two were novel: c.375_376delGT and c.676_679dupAACG, with 57.2% of these potentially occurring de novo. Using all available clinical and genetic data, the National Registry for Greek PJS was established in an attempt to better characterize the syndrome and raise awareness among patients and clinicians (available at https://www.peutzjeghersgreece.org). This is the first comprehensive genetic analysis and clinical characterization of Greek PJS patients, where a high incidence of breast cancer was observed and the first attempt to centralize all data in a National Registry.