CoDP: predicting the impact of unclassified genetic variants in MSH6 by the combination of different properties of the protein

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Familial Intraductal Papillary Mucinous Neoplasm Associated With the Germline MSH6 Missense Variant and Progression of Pancreatic cancer

OBJECTIVES

Intraductal papillary mucinous neoplasm (IPMN) in individuals with at least one first-degree relative with IPMN is defined as familial IPMN. However, few studies have reported on familial IPMN, its clinical characteristics, or the associated genetic factors.

MATERIALS AND METHODS

We report the case of a 58-year-old woman with multifocal IPMN and a mural nodule in the pancreatic body. The patient underwent a distal pancreatectomy and developed pancreatic head cancer 1 year and 6 months postoperatively. The patient had a family history of multifocal IPMN in her father. Therefore, a genetic predisposition to IPMN and pancreatic cancer was suspected. The patient was analyzed for germline variants, and the resected IPMN was subjected to immunohistochemical and somatic variant analyses.

RESULTS

Next-generation sequencing revealed a heterozygous germline missense variant in exon 5 of MSH6 (c.3197A>G; Tyr1066Cys). The pathogenicity of this variant of uncertain significance was suspected based on multiple in silico analyses, and the same MSH6 variant was identified in the patient's father's colonic adenoma. The mural nodule in the pancreatic body was pathologically diagnosed as a high-grade IPMN with ossification and somatic KRAS and PIK3CA variants.

CONCLUSIONS

This case revealed a possible genetic factor for familial IPMN development and presented interesting clinicopathological findings.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

A loss-of-function variant in ZCWPW1 causes human male infertility with sperm head defect and high DNA fragmentation

BACKGROUND

Male infertility is a global health issue. The more causative genes related to human male infertility should be further explored. The essential role of Zcwpw1 in male mouse fertility has been established and the role of ZCWPW1 in human reproduction needs further investigation to verify.

METHODS

An infertile man with oligoasthenoteratozoospermia phenotype and his parents were recruited from West China Second University Hospital, Sichuan University. A total of 200 healthy Han Chinese volunteers without any evidence of infertility were recruited as normal controls, while an additional 150 infertile individuals were included to assess the prevalence of ZCWPW1 variants in a sporadic male sterile population. The causative gene variant was identified by Whole-exome sequencing and Sanger sequencing. The phenotype of the oligoasthenoteratozoospermia was determined by Papanicolaou staining, immunofluorescence staining and electron microscope. In-vitro experiments, western blot and in-silicon analysis were applied to assess the pathogenicity of the identified variant. Additionally, we examined the influence of the variant on the DNA fragmentation and DNA repair capability by Sperm Chromatin Dispersion and Neutral Comet Assay.

RESULTS

The proband exhibits a phenotype of oligoasthenoteratozoospermia, his spermatozoa show head defects by semen examination, Papanicolaou staining and electron microscope assays. Whole-exome sequencing and Sanger sequencing found the proband carries a homozygous ZCWPW1 variant (c.1064C > T, p. P355L). Immunofluorescence analysis shows a significant decrease in ZCWPW1 expression in the proband's sperm. By exogenous expression with ZCWPW1 mutant plasmid in vitro, the obvious declined expression of ZCWPW1 with the mutation is validated in HEK293T. After being treated by hydroxyurea, MUT-ZCWPW1 transfected cells and empty vector transfected cells have a higher level of γ-H2AX, increased tail DNA and reduced H3K9ac level than WT-ZCWPW1 transfected cells. Furthermore, the Sperm Chromatin Dispersion assay revealed the proband's spermatozoa have high DNA fragmentation.

CONCLUSIONS

It is the first report that a novel homozygous missense mutation in ZCWPW1 caused human male infertility with sperm head defects and high DNA fragmentation. This finding enriches the gene variant spectrum and etiology of oligoasthenoteratozoospermia.

Computational approaches for predicting variant impact: An overview from resources, principles to applications

One objective of human genetics is to unveil the variants that contribute to human diseases. With the rapid development and wide use of next-generation sequencing (NGS), massive genomic sequence data have been created, making personal genetic information available. Conventional experimental evidence is critical in establishing the relationship between sequence variants and phenotype but with low efficiency. Due to the lack of comprehensive databases and resources which present clinical and experimental evidence on genotype-phenotype relationship, as well as accumulating variants found from NGS, different computational tools that can predict the impact of the variants on phenotype have been greatly developed to bridge the gap. In this review, we present a brief introduction and discussion about the computational approaches for variant impact prediction. Following an innovative manner, we mainly focus on approaches for non-synonymous variants (nsSNVs) impact prediction and categorize them into six classes. Their underlying rationale and constraints, together with the concerns and remedies raised from comparative studies are discussed. We also present how the predictive approaches employed in different research. Although diverse constraints exist, the computational predictive approaches are indispensable in exploring genotype-phenotype relationship.

Computational study of the impact of nucleotide variations on highly conserved proteins: In the case of actin

Sequencing of individual human genomes enables studying relationship among nucleotide variations, amino acid substitutions, effect on protein structures and diseases. Many studies have found general tendencies, for instance, that pathogenic variations tend to be found in the buried regions of the protein structures, that benign variations tend to be found on the surface of the proteins, and that variations on evolutionary conserved residues tend to be pathogenic. These tendencies were deduced from globular proteins with standard evolutionary changes in amino acid sequences. In this study, we investigated the variation distribution on actin, one of the highly conserved proteins. Many nucleotide variations and three-dimensional structures of actin have been registered in databases. By combining those data, we found that variations buried inside the protein were rather benign and variations on the surface of the protein were pathogenic. This idiosyncratic distribution of the variation impact is likely ascribed to the extensive use of the surface of the protein for protein-protein interactions in actin.

Characterization of the Genomic and Immunological Diversity of Malignant Brain Tumors Through Multi-Sector Analysis

Despite some success in secondary brain metastases, targeted or immune-based therapies have shown limited efficacy against primary brain malignancies such as glioblastoma (GBM). While the intratumoral heterogeneity of GBM is implicated in treatment resistance, it remains unclear whether this diversity is observed within brain metastases and to what extent cancer-cell intrinsic heterogeneity sculpts the local immune microenvironment. Here, we profiled the immunogenomic state of 93 spatially distinct regions from 30 malignant brain tumors through whole exome, RNA, and TCR-sequencing. Our analyses identified differences between primary and secondary malignancies with gliomas displaying more spatial heterogeneity at the genomic and neoantigen level. Additionally, this spatial diversity was recapitulated in the distribution of T cell clones where some gliomas harbored highly expanded but spatially restricted clonotypes. This study defines the immunogenomic landscape across a cohort of malignant brain tumors and contains implications for the design of targeted and immune-based therapies against intracranial malignancies.

VIPdb, a genetic Variant Impact Predictor Database

Genome sequencing identifies vast number of genetic variants. Predicting these variants' molecular and clinical effects is one of the preeminent challenges in human genetics. Accurate prediction of the impact of genetic variants improves our understanding of how genetic information is conveyed to molecular and cellular functions, and is an essential step towards precision medicine. Over one hundred tools/resources have been developed specifically for this purpose. We summarize these tools as well as their characteristics, in the genetic Variant Impact Predictor Database (VIPdb). This database will help researchers and clinicians explore appropriate tools, and inform the development of improved methods. VIPdb can be browsed and downloaded at https://genomeinterpretation.org/vipdb.

A portrait of germline mutation in Brazilian at-risk for hereditary breast cancer

PURPOSE

Knowledge about the germline mutational spectrum among Brazilian with hereditary breast and ovarian cancer (HBOC) is limited. Only five studies have performed comprehensive BRCA sequencing, corresponding to 1041 individuals among a Brazilian population of over 207 million people. Herein we aimed to determine the clinical and molecular characteristics of Brazilian patients who underwent oncogenetic counseling and genetic testing of a panel of high-risk and moderate-risk genes from 2009 to 2017.

METHODS

Massively parallel sequencing was applied in 157 individuals (132 breast cancer-affected and 25 breast cancer-unaffected individuals) selected according NCCN criteria for hereditary breast cancer. Analysis of mutation segregation in family members was performed by capillary bidirectional sequencing, clinical response after treament and survival analysis was estimated by Kaplan-Meier.

RESULTS

Nineteen germline variants were identified,15 pathogenic and 4 VUS (Variants of Uncertain Significance) in 27 individuals (27/157; 17% P < 0.0001) distributed among 7 genes. Sixty-eight percent of patients (13/19) harbor mutation in BRCA genes and 32% (6/19) in moderate risk genes. This is the first study reporting ATR deleterious germline mutation in association with hereditary breast cancer. Cancer-affected patients with moderate- risk mutation present a more aggressive phenotype, with bilateral cancer (25% vs. 13%, P = 0.0305), high-grade tumors (79.2% vs. 46.3%, P = 0.0001) and triple-negative (50% vs. 22.4%, P < 0.0001). However, no difference in the 5 years overall survival was observed between BRCA and moderate risk groups.

CONCLUSIONS

This work highlights the benefits of large-scale sequencing for oncogenetic counseling and extends our understanding about the genetics of hereditary breast cancer in the multi-ethnic Brazilian population.

Must Peutz-Jeghers syndrome patients have the LKB1/STK11 gene mutation? A case report and review of the literature

Peutz-Jeghers syndrome (PJS) is an autosomal dominant inherited disease, which is characterized by mucocutaneous pigmentation and multiple gastrointestinal hamartoma polyps. The germline mutation of LKB1/STK11 gene on chromosome 19p13.3 is considered to be the hereditary cause of PJS. However, must a patient with PJS have the LKB1/STK11 gene mutation? We here report a case of a male patient who had typical manifestations of PJS and a definite family history, but did not have LKB1/STK11 gene mutation. By means of high-throughput sequencing technology, only mutations in APC gene (c.6662T > C: p.Met2221Thr) and MSH6 gene (c.3488A > T: p.Glu1163Val) were detected. The missense mutations in APC and MSH6 gene may lead to abnormalities in structure and function of their expression products, and may result in the occurrence of PJS. This study suggests that some other genetic disorders may cause PJS besides LKB1/STK11 gene mutation.

TumorNext-Lynch-MMR: a comprehensive next generation sequencing assay for the detection of germline and somatic mutations in genes associated with mismatch repair deficiency and Lynch syndrome

The current algorithm for Lynch syndrome diagnosis is highly complex with multiple steps which can result in an extended time to diagnosis while depleting precious tumor specimens. Here we describe the analytical validation of a custom probe-based NGS tumor panel, TumorNext-Lynch-MMR, which generates a comprehensive genetic profile of both germline and somatic mutations that can accelerate and streamline the time to diagnosis and preserve specimen. TumorNext-Lynch-MMR can detect single nucleotide variants, small insertions and deletions in 39 genes that are frequently mutated in Lynch syndrome and colorectal cancer. Moreover, the panel provides microsatellite instability status and detects loss of heterozygosity in the five Lynch genes; MSH2, MSH6, MLH1, PMS2 and EPCAM. Clinical cases are described that highlight the assays ability to differentiate between somatic and germline mutations, precisely classify variants and resolve discordant cases.

Suspected Lynch syndrome associated MSH6 variants: A functional assay to determine their pathogenicity

Lynch syndrome (LS) is a hereditary cancer predisposition caused by inactivating mutations in DNA mismatch repair (MMR) genes. Mutations in the MSH6 DNA MMR gene account for approximately 18% of LS cases. Many LS-associated sequence variants are nonsense and frameshift mutations that clearly abrogate MMR activity. However, missense mutations whose functional implications are unclear are also frequently seen in suspected-LS patients. To conclusively diagnose LS and enroll patients in appropriate surveillance programs to reduce morbidity as well as mortality, the functional consequences of these variants of uncertain clinical significance (VUS) must be defined. We present an oligonucleotide-directed mutagenesis screen for the identification of pathogenic MSH6 VUS. In the screen, the MSH6 variant of interest is introduced into mouse embryonic stem cells by site-directed mutagenesis. Subsequent selection for MMR-deficient cells using the DNA damaging agent 6-thioguanine (6TG) allows the identification of MMR abrogating VUS because solely MMR-deficient cells survive 6TG exposure. We demonstrate the efficacy of the genetic screen, investigate the phenotype of 26 MSH6 VUS and compare our screening results to clinical data from suspected-LS patients carrying these variant alleles.

The colorectal and endometrial cancer predisposition Lynch syndrome (LS) is caused by an inherited heterozygous defect in one of four DNA mismatch repair (MMR) genes. Deleterious mutations (e.g., protein-deleting or -truncating) in DNA MMR genes unambiguously allow for the clinical diagnosis LS and hence enable appropriate surveillance measures to be taken to reduce cancer risk and ensure early detection of tumors. However, currently about one-third of detected MMR gene variants are subtle with less clear functional consequences: missense mutations affecting a single amino acid may be innocuous, hence not causing LS, or partially or fully destroy protein function. As long as uncertainty exists about their pathogenicity, such mutations are labeled ‘variants of uncertain (clinical) significance’ (VUS). VUS hamper genetic counseling and therefore the need for functional testing of VUS is widely recognized. To functionally annotate MMR gene VUS, we have developed a high content cellular assay in which the VUS is introduced in a cell culture by oligonucleotide-directed gene modification. Should the VUS be deleterious for MMR, the modified cells survive exposure to the guanine analog 6-thioguanine (6TG) and 6TG-resistant colonies appear. Should the mutation not affect MMR, no colonies appear. Here we present the adaptation and application of this protocol to the functional annotation of variants of the MMR gene MSH6. Implementation of our assay in clinical genetics laboratories will provide clinicians with information for proper counseling of mutation carriers and treatment of their of tumors.

Functional analysis of rare variants in mismatch repair proteins augments results from computation-based predictive methods

The cancer-predisposing Lynch Syndrome (LS) arises from germline mutations in DNA mismatch repair (MMR) genes, predominantly MLH1, MSH2, MSH6, and PMS2. A major challenge for clinical diagnosis of LS is the frequent identification of variants of uncertain significance (VUS) in these genes, as it is often difficult to determine variant pathogenicity, particularly for missense variants. Generic programs such as SIFT and PolyPhen-2, and MMR gene-specific programs such as PON-MMR and MAPP-MMR, are often used to predict deleterious or neutral effects of VUS in MMR genes. We evaluated the performance of multiple predictive programs in the context of functional biologic data for 15 VUS in MLH1, MSH2, and PMS2. Using cell line models, we characterized VUS predicted to range from neutral to pathogenic on mRNA and protein expression, basal cellular viability, viability following treatment with a panel of DNA-damaging agents, and functionality in DNA damage response (DDR) signaling, benchmarking to wild-type MMR proteins. Our results suggest that the MMR gene-specific classifiers do not always align with the experimental phenotypes related to DDR. Our study highlights the importance of complementary experimental and computational assessment to develop future predictors for the assessment of VUS.

Integrative analysis of hereditary nonpolyposis colorectal cancer: the contribution of allele-specific expression and other assays to diagnostic algorithms

The identification of germline variants predisposing to hereditary nonpolyposis colorectal cancer (HNPCC) is crucial for clinical management of carriers, but several probands remain negative for such variants or bear variants of uncertain significance (VUS). Here we describe the results of integrative molecular analyses in 132 HNPCC patients providing evidences for improved genetic testing of HNPCC with traditional or next generation methods. Patients were screened for: germline allele-specific expression (ASE), nucleotide variants, rearrangements and promoter methylation of mismatch repair (MMR) genes; germline EPCAM rearrangements; tumor microsatellite instability (MSI) and immunohistochemical (IHC) MMR protein expression. Probands negative for pathogenic variants of MMR genes were screened for germline APC and MUTYH sequence variants. Most germline defects identified were sequence variants and rearrangements of MMR genes. Remarkably, altered germline ASE of MMR genes was detected in 8/22 (36.5%) probands analyzed, including 3 cases negative at other screenings. Moreover, ASE provided evidence for the pathogenic role and guided the characterization of a VUS shared by 2 additional probands. No germline MMR gene promoter methylation was observed and only one EPCAM rearrangement was detected. In several cases, tumor IHC and MSI diverged from germline screening results. Notably, APC or biallelic MUTYH germline defects were identified in 2/19 probands negative for pathogenic variants of MMR genes. Our results show that ASE complements gDNA-based analyses in the identification of MMR defects and in the characterization of VUS affecting gene expression, increasing the number of germline alterations detected. An appreciable fraction of probands negative for MMR gene variants harbors APC or MUTYH variants. These results indicate that germline ASE analysis and screening for APC and MUTYH defects should be included in HNPCC diagnostic algorithms.

Molecular and clinical characteristics of MSH6 germline variants detected in colorectal cancer patients

The MSH6 gene is one of the mismatch repair genes involved in Lynch syndrome and its mutations account for 10-20% of Lynch syndrome. Although previous studies suggested that the difference of the geographical region affects the clinical phenotype of Lynch syndrome, there has been no report on the detailed features of Japanese Lynch syndrome patients carrying an MSH6 mutation. The aim of the present study was to investigate the clinical and molecular features of MSH6 mutation carriers in Japan. Surgically resected 1720 colorectal carcinoma specimens were screened by microsatellite instability (MSI) testing and the MSI-high cases were subjected to a germline mutation analysis of the mismatch repair genes MLH1, MSH2 and MSH6. We investigated the clinical and molecular features of the MSH6 variants, such as the family cancer history, pathological findings, immunohistochemistry, methylation status of the MLH1 promoter and BRAF mutation in the colorectal tumor. Furthermore, the impact of the missense variants on MSH6 protein was predicted by using in silico tools. We identified nine novel pathogenic mutations and eight unclassified missense variants. Among the eight missense variants, three were suspected pathogenic by in silico analysis. We also found that most colorectal cancers in the MSH6 mutation carrier were diagnosed after the age of 50 and were localized distally. Furthermore, the mean age at diagnosis of endometrial cancer in Japanese MSH6 mutation carriers (49.2 years) was earlier than previous reports from Western countries (56.5 years). These results may improve the surveillance program for Japanese MSH6 mutation carriers.