Mathematical modeling of Lynch syndrome carcinogenesis

Cancer is one of the leading causes of disease-related death worldwide. In recent years, large amounts of data on cancer genetics and molecular characteristics have become available and accumulated with increasing speed. However, the current understanding of cancer as a disease is still limited by the lack of suitable models that allow interpreting these data in proper ways. Thus, the highly interdisciplinary research field of mathematical oncology has evolved to use mathematics, modeling, and simulations to study cancer with the overall goal to improve clinical patient care. This dissertation aims at developing mathematical models and tools for different spatial scales of cancer development at the example of colorectal cancer in Lynch syndrome, the most common inherited colorectal cancer predisposition syndrome. We derive model-driven approaches for carcinogenesis at the DNA, cell, and crypt level, as well as data-driven methods for cancer-immune interactions at the DNA level and for the evaluation of diagnostic procedures at the Lynch syndrome population level. The developed models present an important step toward an improved understanding of hereditary cancer as a disease aiming at rapid implementation into clinical management guidelines and into the development of novel, innovative approaches for prevention and treatment

A simple approach for detecting HLA-A02 alleles in archival formalin-fixed paraffin-embedded tissue samples and an application example for studying cancer immunoediting

The HLA system represents a central component of the antigen presentation machinery. As every patient possesses a defined set of HLA molecules, only certain antigens can be presented on the cell surface. Thus, studying HLA type-dependent antigen presentation can improve the understanding of variation in susceptibility to various diseases, including infectious diseases and cancer. In archival formalin-fixed paraffin-embedded (FFPE) tissue, the HLA type is difficult to analyze because of fragmentation of DNA, hindering the application of commonly used assays that rely on long DNA stretches. Addressing these difficulties, we present a refined approach for characterizing presence or absence of HLA-A*02, the most common HLA-A allele in the Caucasian population, in archival samples. We validated our genotyping strategy in a cohort of 90 samples with HLA status obtained by an NGS-based method. 90% (n = 81) of the samples could be analyzed with the approach. For all of them, the presence or absence of HLA-A*02 alleles was correctly determined with the method, demonstrating 100% sensitivity and specificity (95% CI: 91.40%-100% and 91.19%-100%). Furthermore, we provide an example of application in an independent cohort of 73 FFPE microsatellite-unstable (MSI) colorectal cancer samples. As MSI cancer cells encompass a high number of mutations in coding microsatellites, leading to the generation of highly immunogenic frameshift peptide antigens, they are ideally suited for studying relations between the mutational landscape of tumor cells and interindividual differences in the immune system, including the HLA genotype. Overall, our method can help to promote studying HLA type-dependency during the pathogenesis of a wide range of diseases, making archival and historic tissue samples accessible for identifying HLA-A*02 alleles.Peer reviewe

Is HLA type a possible cancer risk modifier in Lynch syndrome?

Lynch syndrome (LS) is the most common inherited cancer syndrome. It is inherited via a monoallelic germline variant in one of the DNA mismatch repair (MMR) genes. LS carriers have a broad 30% to 80% risk of developing various malignancies, and more precise, individual risk estimations would be of high clinical value, allowing tailored cancer prevention and surveillance. Due to MMR deficiency, LS cancers are characterized by the accumulation of frameshift mutations leading to highly immunogenic frameshift peptides (FSPs). Thus, immune surveillance is proposed to inhibit the outgrowth of MMR-deficient cell clones. Recent studies have shown that immunoediting during the evolution of MMR-deficient cancers leads to a counter-selection of highly immunogenic antigens. The immunogenicity of FSPs is dependent on the antigen presentation. One crucial factor determining antigen presentation is the HLA genotype. Hence, a LS carrier's HLA genotype plays an important role in the presentation of FSP antigens to the immune system, and may influence the likelihood of progression from precancerous lesions to cancer. To address the challenge of clarifying this possibility including diverse populations with different HLA types, we have established the INDICATE initiative (Individual cancer risk by HLA type, ), an international network aiming at a systematic evaluation of the HLA genotype as a possible cancer risk modifier in LS. Here we summarize the current knowledge on the role of HLA type in cancer risk and outline future research directions to delineate possible association in the scenario of LS with genetically defined risk population and highly immunogenic tumors.Peer reviewe

The shared frameshift mutation landscape of microsatellite-unstable cancers suggests immunoediting during tumor evolution

The immune system can recognize and attack cancer cells, especially those with a high load of mutation-induced neoantigens. Such neoantigens are abundant in DNA mismatch repair (MMR)-deficient, microsatellite-unstable (MSI) cancers. MMR deficiency leads to insertion/deletion (indel) mutations at coding microsatellites (cMS) and to neoantigen-inducing translational frameshifts. Here, we develop a tool to quantify frameshift mutations in MSI colorectal and endometrial cancer. Our results show that frameshift mutation frequency is negatively correlated to the predicted immunogenicity of the resulting peptides, suggesting counterselection of cell clones with highly immunogenic frameshift peptides. This correlation is absent in tumors with Beta-2-microglobulin mutations, and HLA-A*02:01 status is related to cMS mutation patterns. Importantly, certain outlier mutations are common in MSI cancers despite being related to frameshift peptides with functionally confirmed immunogenicity, suggesting a possible driver role during MSI tumor evolution. Neoantigens resulting from shared mutations represent promising vaccine candidates for prevention of MSI cancers. DNA mismatch repair (MMR)-deficient cancers with microsatellite-instability are characterized by a high load of frameshift mutation-derived neoantigens. Here, by mapping the frameshift mutation landscape and predicting the immunogenicity of the resulting peptides, the authors show evidence of immunoediting in MMR-deficient colorectal and endometrial cancers.Peer reviewe

Fertility preservation for female patients with childhood, adolescent, and young adult cancer:recommendations from the PanCareLIFE Consortium and the International Late Effects of Childhood Cancer Guideline Harmonization Group

Female patients with childhood, adolescent, and young adult cancer are at increased risk for fertility impairment when treatment adversely affects the function of reproductive organs. Patients and their families desire biological children but substantial variations in clinical practice guidelines reduce consistent and timely implementation of effective interventions for fertility preservation across institutions. As part of the PanCareLIFE Consortium, and in collaboration with the International Late Effects of Childhood Cancer Guideline Harmonization Group, we reviewed the current literature and developed a clinical practice guideline for fertility preservation in female patients who were diagnosed with childhood, adolescent, and young adult cancer at age 25 years or younger, including guidance on risk assessment and available methods for fertility preservation. The Grading of Recommendations Assessment, Development and Evaluation methodology was used to grade the available evidence and to form the recommendations. This clinical practice guideline leverages existing evidence and international expertise to develop transparent recommendations that are easy to use to facilitate the care of female patients with childhood, adolescent, and young adult cancer who are at high risk for fertility impairment. A complete review of the existing evidence, including a quality assessment, transparent reporting of the guideline panel's decisions, and achievement of global interdisciplinary consensus, is an important result of this intensive collaboration.info:eu-repo/semantics/publishe

Possible modification of BRSK1 on the risk of alkylating chemotherapy-related reduced ovarian function

STUDY QUESTION: Do genetic variations in the DNA damage response pathway modify the adverse effect of alkylating agents on ovarian function in female childhood cancer survivors (CCS)? SUMMARY ANSWER: Female CCS carrying a common BR serine/threonine kinase 1 (BRSK1) gene variant appear to be at 2.5-fold increased odds of reduced ovarian function after treatment with high doses of alkylating chemotherapy. WHAT IS KNOWN ALREADY: Female CCS show large inter-individual variability in the impact of DNA-damaging alkylating chemotherapy, given as treatment of childhood cancer, on adult ovarian function. Genetic variants in DNA repair genes affecting ovarian function might explain this variability. STUDY DESIGN, SIZE, DURATION: CCS for the discovery cohort were identified from the Dutch Childhood Oncology Group (DCOG) LATER VEVO-study, a multi-centre retrospective cohort study evaluating fertility, ovarian reserve and risk of premature menopause among adult female 5-year survivors of childhood cancer. Female 5-year CCS, diagnosed with cancer and treated with chemotherapy before the age of 25 years, and aged 18 years or older at time of study were enrolled in the current study. Results from the discovery Dutch DCOG-LATER VEVO cohort (n = 285) were validated in the pan-European PanCareLIFE (n =465) and the USA-based St. Jude Lifetime Cohort (n = 391). PARTICIPANTS/MATERIALS, SETTING, METHODS: To evaluate ovarian function, anti-Miillerian hormone (AMH) levels were assessed in both the discovery cohort and the replication cohorts. Using additive genetic models in linear and logistic regression, five genetic variants involved in DNA damage response were analysed in relation to cyclophosphamide equivalent dose (CED) score and their impact on ovarian function. Results were then examined using fixed-effect meta-analysis. MAIN RESULTS AND THE ROLE OF CHANCE: Meta-analysis across the three independent cohorts showed a significant interaction effect (P= 3.0 x 10(-4)) between rs11668344 of BRSK 1 (allele frequency = 0.34) among CCS treated with high-dose alkylating agents (CED score >= 8000 mg/m(2)), resulting in a 2.5-fold increased odds of a reduced ovarian function (lowest AMH tertile) for CCS carrying one G allele compared to CCS without this allele (odds ratio genotype AA: 2.01 vs AG: 5.00). LIMITATIONS, REASONS FOR CAUTION: While low AMH levels can also identify poor responders in assisted reproductive technology, it needs to be emphasized that AMH remains a surrogate marker of ovarian function. WIDER IMPLICATIONS OF THE FINDINGS: Further research, validating our findings and identifying additional risk contributing genetic variants, may enable individualized counselling regarding treatment-related risks and necessity of fertility preservation procedures in girls with cancer

Colorectal cancer incidences in Lynch syndrome: a comparison of results from the prospective lynch syndrome database and the international mismatch repair consortium

Objective To compare colorectal cancer (CRC) incidences in carriers of pathogenic variants of the MMR genes in the PLSD and IMRC cohorts, of which only the former included mandatory colonoscopy surveillance for all participants. Methods CRC incidences were calculated in an intervention group comprising a cohort of confirmed carriers of pathogenic or likely pathogenic variants in mismatch repair genes (path_MMR) followed prospectively by the Prospective Lynch Syndrome Database (PLSD). All had colonoscopy surveillance, with polypectomy when polyps were identified. Comparison was made with a retrospective cohort reported by the International Mismatch Repair Consortium (IMRC). This comprised confirmed and inferred path_MMR carriers who were first- or second-degree relatives of Lynch syndrome probands. Results In the PLSD, 8,153 subjects had follow-up colonoscopy surveillance for a total of 67,604 years and 578 carriers had CRC diagnosed. Average cumulative incidences of CRC in path_MLH1 carriers at 70 years of age were 52% in males and 41% in females; for path_MSH2 50% and 39%; for path_MSH6 13% and 17% and for path_PMS2 11% and 8%. In contrast, in the IMRC cohort, corresponding cumulative incidences were 40% and 27%; 34% and 23%; 16% and 8% and 7% and 6%. Comparing just the European carriers in the two series gave similar findings. Numbers in the PLSD series did not allow comparisons of carriers from other continents separately. Cumulative incidences at 25 years were < 1% in all retrospective groups. Conclusions Prospectively observed CRC incidences (PLSD) in path_MLH1 and path_MSH2 carriers undergoing colonoscopy surveillance and polypectomy were higher than in the retrospective (IMRC) series, and were not reduced in path_MSH6 carriers. These findings were the opposite to those expected. CRC point incidence before 50 years of age was reduced in path_PMS2 carriers subjected to colonoscopy, but not significantly so

Colorectal cancer incidences in Lynch syndrome : a comparison of results from the prospective lynch syndrome database and the international mismatch repair consortium

Objective: To compare colorectal cancer (CRC) incidences in carriers of pathogenic variants of the MMR genes in the PLSD and IMRC cohorts, of which only the former included mandatory colonoscopy surveillance for all participants. Methods: CRC incidences were calculated in an intervention group comprising a cohort of confirmed carriers of pathogenic or likely pathogenic variants in mismatch repair genes (path_MMR) followed prospectively by the Prospective Lynch Syndrome Database (PLSD). All had colonoscopy surveillance, with polypectomy when polyps were identified. Comparison was made with a retrospective cohort reported by the International Mismatch Repair Consortium (IMRC). This comprised confirmed and inferred path_MMR carriers who were first- or second-degree relatives of Lynch syndrome probands. Results: In the PLSD, 8,153 subjects had follow-up colonoscopy surveillance for a total of 67,604 years and 578 carriers had CRC diagnosed. Average cumulative incidences of CRC in path_MLH1 carriers at 70 years of age were 52% in males and 41% in females; for path_MSH2 50% and 39%; for path_MSH6 13% and 17% and for path_PMS2 11% and 8%. In contrast, in the IMRC cohort, corresponding cumulative incidences were 40% and 27%; 34% and 23%; 16% and 8% and 7% and 6%. Comparing just the European carriers in the two series gave similar findings. Numbers in the PLSD series did not allow comparisons of carriers from other continents separately. Cumulative incidences at 25 years were < 1% in all retrospective groups. Conclusions: Prospectively observed CRC incidences (PLSD) in path_MLH1 and path_MSH2 carriers undergoing colonoscopy surveillance and polypectomy were higher than in the retrospective (IMRC) series, and were not reduced in path_MSH6 carriers. These findings were the opposite to those expected. CRC point incidence before 50 years of age was reduced in path_PMS2 carriers subjected to colonoscopy, but not significantly so.publishedVersionPeer reviewe

Dominantly inherited micro-satellite instable cancer – the four Lynch syndromes - an EHTG, PLSD position statement

The recognition of dominantly inherited micro-satellite instable (MSI) cancers caused by pathogenic variants in one of the four mismatch repair (MMR) genes MSH2, MLH1, MSH6 and PMS2 has modified our understanding of carcinogenesis. Inherited loss of function variants in each of these MMR genes cause four dominantly inherited cancer syndromes with different penetrance and expressivities: the four Lynch syndromes. No person has an “average sex “or a pathogenic variant in an “average Lynch syndrome gene” and results that are not stratified by gene and sex will be valid for no one. Carcinogenesis may be a linear process from increased cellular division to localized cancer to metastasis. In addition, in the Lynch syndromes (LS) we now recognize a dynamic balance between two stochastic processes: MSI producing abnormal cells, and the host’s adaptive immune system’s ability to remove them. The latter may explain why colonoscopy surveillance does not reduce the incidence of colorectal cancer in LS, while it may improve the prognosis. Most early onset colon, endometrial and ovarian cancers in LS are now cured and most cancer related deaths are after subsequent cancers in other organs. Aspirin reduces the incidence of colorectal and other cancers in LS. Immunotherapy increases the host immune system’s capability to destroy MSI cancers. Colonoscopy surveillance, aspirin prevention and immunotherapy represent major steps forward in personalized precision medicine to prevent and cure inherited MSI cancer

Mathematical modeling of multiple pathways in colorectal carcinogenesis using dynamical systems with Kronecker structure.

Like many other types of cancer, colorectal cancer (CRC) develops through multiple pathways of carcinogenesis. This is also true for colorectal carcinogenesis in Lynch syndrome (LS), the most common inherited CRC syndrome. However, a comprehensive understanding of the distribution of these pathways of carcinogenesis, which allows for tailored clinical treatment and even prevention, is still lacking. We suggest a linear dynamical system modeling the evolution of different pathways of colorectal carcinogenesis based on the involved driver mutations. The model consists of different components accounting for independent and dependent mutational processes. We define the driver gene mutation graphs and combine them using the Cartesian graph product. This leads to matrix components built by the Kronecker sum and product of the adjacency matrices of the gene mutation graphs enabling a thorough mathematical analysis and medical interpretation. Using the Kronecker structure, we developed a mathematical model which we applied exemplarily to the three pathways of colorectal carcinogenesis in LS. Beside a pathogenic germline variant in one of the DNA mismatch repair (MMR) genes, driver mutations in APC, CTNNB1, KRAS and TP53 are considered. We exemplarily incorporate mutational dependencies, such as increased point mutation rates after MMR deficiency, and based on recent experimental data, biallelic somatic CTNNB1 mutations as common drivers of LS-associated CRCs. With the model and parameter choice, we obtained simulation results that are in concordance with clinical observations. These include the evolution of MMR-deficient crypts as early precursors in LS carcinogenesis and the influence of variants in MMR genes thereon. The proportions of MMR-deficient and MMR-proficient APC-inactivated crypts as first measure for the distribution among the pathways in LS-associated colorectal carcinogenesis are compatible with clinical observations. The approach provides a modular framework for modeling multiple pathways of carcinogenesis yielding promising results in concordance with clinical observations in LS CRCs