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.

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/.

Clonal heterogeneity and rates of specific chromosome gains are risk predictors in childhood high-hyperdiploid B-cell acute lymphoblastic leukemia.

B-cell acute lymphoblastic leukemia (B-ALL) is the commonest childhood cancer. High hyperdiploidy (HHD) identifies the most frequent cytogenetic subgroup in childhood B-ALL. Although hyperdiploidy represents an important prognostic factor in childhood B-ALL, the specific chromosome gains with prognostic value in HHD-B-ALL remain controversial, and the current knowledge about the hierarchy of chromosome gains, clonal heterogeneity and chromosomal instability in HHD-B-ALL remains very limited. We applied automated sequential-iFISH coupled with single-cell computational modeling to identify the specific chromosomal gains of the eight typically gained chromosomes in a large cohort of 72 primary diagnostic (DX, n = 62) and matched relapse (REL, n = 10) samples from HHD-B-ALL patients with either favorable or unfavorable clinical outcome in order to characterize the clonal heterogeneity, specific chromosome gains and clonal evolution. Our data show a high degree of clonal heterogeneity and a hierarchical order of chromosome gains in DX samples of HHD-B-ALL. The rates of specific chromosome gains and clonal heterogeneity found in DX samples differ between HHD-B-ALL patients with favorable or unfavorable clinical outcome. In fact, our comprehensive analyses at DX using a computationally defined risk predictor revealed low levels of trisomies +18+10 and low levels of clonal heterogeneity as robust relapse risk factors in minimal residual disease (MRD)-negative childhood HHD-B-ALL patients: relapse-free survival beyond 5 years: 22.1% versus 87.9%, P < 0.0001 and 33.3% versus 80%, P < 0.0001, respectively. Moreover, longitudinal analysis of matched DX-REL HHD-B-ALL samples revealed distinct patterns of clonal evolution at relapse. Our study offers a reliable prognostic sub-stratification of pediatric MRD-negative HHD-B-ALL patients.