Comprehensive Evaluation of Protein Coding Mononucleotide Microsatellites in Microsatellite-Unstable Colorectal Cancer

What can we learn about acid-base transporters in cancer from studying somatic mutations in their genes?

Acidosis is a chemical signature of the tumour microenvironment that challenges intracellular pH homeostasis. The orchestrated activity of acid-base transporters of the solute-linked carrier (SLC) family is critical for removing the end-products of fermentative metabolism (lactate/H+) and maintaining a favourably alkaline cytoplasm. Given the critical role of pH homeostasis in enabling cellular activities, mutations in relevant SLC genes may impact the oncogenic process, emerging as negatively or positively selected, or as driver or passenger mutations. To address this, we performed a pan-cancer analysis of The Cancer Genome Atlas simple nucleotide variation data for acid/base-transporting SLCs (ABT-SLCs). Somatic mutation patterns of monocarboxylate transporters (MCTs) were consistent with their proposed essentiality in facilitating lactate/H+ efflux. Among all cancers, tumours of uterine corpus endometrial cancer carried more ABT-SLC somatic mutations than expected from median tumour mutation burden. Among these, somatic mutations in SLC4A3 had features consistent with meaningful consequences on cellular fitness. Definitive evidence for ABT-SLCs as 'cancer essential' or 'driver genes' will have to consider microenvironmental context in genomic sequencing because bulk approaches are insensitive to pH heterogeneity within tumours. Moreover, genomic analyses must be validated with phenotypic outcomes (i.e. SLC-carried flux) to appreciate the opportunities for targeting acid-base transport in cancers.

Profiling the molecular and clinical landscape of glioblastoma utilizing the Oncology Research Information Exchange Network brain cancer database

Background

Glioblastoma exhibits aggressive growth and poor outcomes despite treatment, and its marked variability renders therapeutic design and prognostication challenging. The Oncology Research Information Exchange Network (ORIEN) database contains complementary clinical, genomic, and transcriptomic profiling of 206 glioblastoma patients, providing opportunities to identify novel associations between molecular features and clinical outcomes.

Methods

Survival analyses were performed using the Logrank test, and clinical features were evaluated using Wilcoxon and chi-squared tests with q-values derived via Benjamini-Hochberg correction. Mutational analyses utilized sample-level enrichments from whole exome sequencing data, and statistical tests were performed using the one-sided Fisher Exact test with Benjamini-Hochberg correction. Transcriptomic analyses utilized a student's t-test with Benjamini-Hochberg correction. Expression fold changes were processed with Ingenuity Pathway Analysis to determine pathway-level alterations between groups.

Results

Key findings include an association of MUC17, SYNE1, and TENM1 mutations with prolonged overall survival (OS); decreased OS associated with higher epithelial growth factor receptor (EGFR) mRNA expression, but not with EGFR amplification or mutation; a 14-transcript signature associated with OS > 2 years; and 2 transcripts associated with OS < 1 year.

Conclusions

Herein, we report the first clinical, genomic, and transcriptomic analysis of ORIEN glioblastoma cases, incorporating sample reclassification under updated 2021 diagnostic criteria. These findings create multiple avenues for further investigation and reinforce the value of multi-institutional consortia such as ORIEN in deepening our knowledge of intractable diseases such as glioblastoma.

Single-Cell Mononucleotide Microsatellite Analysis Reveals Differential Insertion-Deletion Dynamics in Mouse T Cells

Microsatellite sequences are particularly prone to slippage during DNA replication, forming insertion-deletion loops that, if left unrepaired, result in de novo mutations (expansions or contractions of the repeat array). Mismatch repair (MMR) is a critical DNA repair mechanism that corrects these insertion-deletion loops, thereby maintaining microsatellite stability. MMR deficiency gives rise to the molecular phenotype known as microsatellite instability (MSI). By sequencing MMR-proficient and -deficient (Mlh1+/+ and Mlh1−/−) single-cell exomes from mouse T cells, we reveal here several previously unrecognized features of in vivo MSI. Specifically, mutational dynamics of insertions and deletions were different on multiple levels. Factors that associated with propensity of mononucleotide microsatellites to insertions versus deletions were: microsatellite length, nucleotide composition of the mononucleotide tract, gene length and transcriptional status, as well replication timing. Here, we show on a single-cell level that deletions — the predominant MSI type in MMR-deficient cells — are preferentially associated with longer A/T tracts, long or transcribed genes and later-replicating genes.

Predicting stability of DNA bulge at mononucleotide microsatellite

Mononucleotide microsatellites are clinically and forensically crucial DNA sequences due to their high mutability and abundance in the human genome. As a mutagenic intermediate of an indel in a microsatellite and a consequence of probe hybridization after such mutagenesis, a bulge with structural degeneracy sliding within a microsatellite is formed. Stability of such dynamic bulges, however, is still poorly understood despite their critical role in cancer genomics and neurological disease studies. In this paper, we have built a model that predicts the thermodynamics of a sliding bulge at a microsatellite. We first identified 40 common bulge states that can be assembled into any sliding bulges, and then characterized them with toehold exchange energy measurement and the partition function. Our model, which is the first to predict the free energy of sliding bulges with more than three repeats, can infer the stability penalty of a sliding bulge of any sequence and length with a median prediction error of 0.22 kcal/mol. Patterns from the prediction clearly explain landscapes of microsatellites observed in the literature, such as higher mutation rates of longer microsatellites and C/G microsatellites.

Functional characterization of the sodium/hydrogen exchanger 8 and its role in proliferation of colonic epithelial cells

Intestinal NaCl, HCO₃^-, and fluid absorption are strongly dependent on apical Na⁺/H⁺ exchange. The intestine expresses three presumably apical sodium-hydrogen exchanger (NHE) isoforms: NHE2, NHE3, and NHE8. We addressed the role of NHE8 [solute carrier 9A8 (SLC9A8)] and its interplay with NHE2 (SLC9A2) in luminal proton extrusion during acute and chronic enterocyte acidosis and studied the differential effects of NHE8 and NHE2 on enterocyte proliferation. In contrast to NHE3, which was upregulated in differentiated versus undifferentiated colonoids, the expression of NHE2 and NHE8 remained constant during differentiation of colonoids and Caco2Bbe cells. Heterogeneously expressed Flag-tagged rat (r)Nhe8 and human (h)NHE8 translocated to the apical membrane of Caco2Bbe cells. rNhe8 and hNHE8, when expressed in NHE-deficient PS120 fibroblasts showed higher sensitivity to HOE642 compared to NHE2. Lentiviral shRNA knockdown of endogenous NHE2 in Caco2Bbe cells (C2Bbe/shNHE2) resulted in a decreased steady-state intracellular pH (pH_i), an increased NHE8 mRNA expression, and augmented NHE8-mediated apical NHE activity. Lentiviral shRNA knockdown of endogenous NHE8 in Caco2Bbe cells (C2Bbe/shNHE8) resulted in a decreased steady-state pH_i as well, accompanied by decreased NHE2 mRNA expression and activity, which together contributed to reduced apical NHE activity in the NHE8-knockdown cells. Chronic acidosis increased NHE8 but not NHE2 mRNA expression. Alterations in NHE2 and NHE8 expression/activity affected proliferation, with C2Bbe/shNHE2 cells having lower and C2Bbe/shNHE8 having higher proliferative capacity, accompanied by amplified ERK1/2 signaling pathway and increased EGFR expression in the latter cell line. Thus, both Na⁺/H⁺ exchangers have distinct functions during cellular homeostasis by triggering different signaling pathways to regulate cellular proliferation and pH_i control.

No evidence of EMAST in whole genome sequencing data from 248 colorectal cancers

Microsatellite instability (MSI) is caused by defective DNA mismatch repair (MMR), and manifests as accumulation of small insertions and deletions (indels) in short tandem repeats of the genome. Another form of repeat instability, elevated microsatellite alterations at selected tetranucleotide repeats (EMAST), has been suggested to occur in 50% to 60% of colorectal cancer (CRC), of which approximately one quarter are accounted for by MSI. Unlike for MSI, the criteria for defining EMAST is not consensual. EMAST CRCs have been suggested to form a distinct subset of CRCs that has been linked to a higher tumor stage, chronic inflammation, and poor prognosis. EMAST CRCs not exhibiting MSI have been proposed to show instability of di- and trinucleotide repeats in addition to tetranucleotide repeats, but lack instability of mononucleotide repeats. However, previous studies on EMAST have been based on targeted analysis of small sets of marker repeats, often in relatively few samples. To gain insight into tetranucleotide instability on a genome-wide level, we utilized whole genome sequencing data from 227 microsatellite stable (MSS) CRCs, 18 MSI CRCs, 3 POLE-mutated CRCs, and their corresponding normal samples. As expected, we observed tetranucleotide instability in all MSI CRCs, accompanied by instability of mono-, di-, and trinucleotide repeats. Among MSS CRCs, some tumors displayed more microsatellite mutations than others as a continuum, and no distinct subset of tumors with the previously proposed molecular characters of EMAST could be observed. Our results suggest that tetranucleotide repeat mutations in non-MSI CRCs represent stochastic mutation events rather than define a distinct CRC subclass.

Whole genome sequencing of metastatic colorectal cancer reveals prior treatment effects and specific metastasis features

In contrast to primary colorectal cancer (CRC) little is known about the genomic landscape of metastasized CRC. Here we present whole genome sequencing data of metastases of 429 CRC patients participating in the pan-cancer CPCT-02 study (NCT01855477). Unsupervised clustering using mutational signature patterns highlights three major patient groups characterized by signatures known from primary CRC, signatures associated with received prior treatments, and metastasis-specific signatures. Compared to primary CRC, we identify additional putative (non-coding) driver genes and increased frequencies in driver gene mutations. In addition, we identify specific genes preferentially affected by microsatellite instability. CRC-specific 1kb-10Mb deletions, enriched for common fragile sites, and LINC00672 mutations are associated with response to treatment in general, whereas FBXW7 mutations predict poor response specifically to EGFR-targeted treatment. In conclusion, the genomic landscape of mCRC shows defined changes compared to primary CRC, is affected by prior treatments and contains features with potential clinical relevance.

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.

Unique genomic and neoepitope landscapes across tumors: a study across time, tissues, and space within a single lynch syndrome patient

Lynch syndrome (LS) arises in patients with pathogenic germline variants in DNA mismatch repair genes. LS is the most common inherited cancer predisposition condition and confers an elevated lifetime risk of multiple cancers notably colorectal and endometrial carcinomas. A distinguishing feature of LS associated tumors is accumulation of variants targeting microsatellite repeats and the potential for high tumor specific neoepitope levels. Recurrent somatic variants targeting a small subset of genes have been identified in tumors with microsatellite instability. Notably these include frameshifts that can activate immune responses and provide vaccine targets to affect the lifetime cancer risk associated with LS. However the presence and persistence of targeted neoepitopes across multiple tumors in single LS patients has not been rigorously studied. Here we profiled the genomic landscapes of five distinct treatment naïve tumors, a papillary transitional cell renal cell carcinoma, a duodenal carcinoma, two metachronous colorectal carcinomas, and multi-regional sampling in a triple-negative breast tumor, arising in a LS patient over 10 years. Our analyses suggest each tumor evolves a unique complement of variants and that vaccines based on potential neoepitopes from one tissue may not be effective across all tumors that can arise during the lifetime of LS patients.

Stabilize and connect: the role of LARP7 in nuclear non-coding RNA metabolism

La and La-related proteins (LARPs) are characterized by a common RNA interaction platform termed the La module. This structural hallmark allows LARPs to pervade various aspects of RNA biology. The metazoan LARP7 protein binds to the 7SK RNA as part of a 7SK small nuclear ribonucleoprotein (7SK snRNP), which inhibits the transcriptional activity of RNA polymerase II (Pol II). Additionally, recent findings revealed unanticipated roles of LARP7 in the assembly of other RNPs, as well as in the modification, processing and cellular transport of RNA molecules. Reduced levels of functional LARP7 have been linked to cancer and Alazami syndrome, two seemingly unrelated human diseases characterized either by hyperproliferation or growth retardation. Here, we review the intricate regulatory networks centered on LARP7 and assess how malfunction of these networks may relate to the etiology of LARP7-linked diseases.

Comprehensive analysis of indels in whole-genome microsatellite regions and microsatellite instability across 21 cancer types

Microsatellites are repeats of 1- to 6-bp units, and approximately 10 million microsatellites have been identified across the human genome. Microsatellites are vulnerable to DNA mismatch errors and have thus been used to detect cancers with mismatch repair deficiency. To reveal the mutational landscape of microsatellite repeat regions at the genome level, we analyzed approximately 20.1 billion microsatellites in 2717 whole genomes of pan-cancer samples across 21 tissue types. First, we developed a new insertion and deletion caller (MIMcall) that takes into consideration the error patterns of different types of microsatellites. Among the 2717 pan-cancer samples, our analysis identified 31 samples, including colorectal, uterus, and stomach cancers, with a higher proportion of mutated microsatellite (≥0.03), which we defined as microsatellite instability (MSI) cancers of genome-wide level. Next, we found 20 highly mutated microsatellites that can be used to detect MSI cancers with high sensitivity. Third, we found that replication timing and DNA shape were significantly associated with mutation rates of microsatellites. Last, analysis of mutations in mismatch repair genes showed that somatic SNVs and short indels had larger functional impacts than germline mutations and structural variations. Our analysis provides a comprehensive picture of mutations in the microsatellite regions and reveals possible causes of mutations, as well as provides a useful marker set for MSI detection.

Comprehensive evaluation of coding region point mutations in microsatellite-unstable colorectal cancer

Microsatellite instability (MSI) leads to accumulation of an excessive number of mutations in the genome, mostly small insertions and deletions. MSI colorectal cancers (CRCs), however, also contain more point mutations than microsatellite‐stable (MSS) tumors, yet they have not been as comprehensively studied. To identify candidate driver genes affected by point mutations in MSI CRC, we ranked genes based on mutation significance while correcting for replication timing and gene expression utilizing an algorithm, MutSigCV. Somatic point mutation data from the exome kit‐targeted area from 24 exome‐sequenced sporadic MSI CRCs and respective normals, and 12 whole‐genome‐sequenced sporadic MSI CRCs and respective normals were utilized. The top 73 genes were validated in 93 additional MSI CRCs. The MutSigCV ranking identified several well‐established MSI CRC driver genes and provided additional evidence for previously proposed CRC candidate genes as well as shortlisted genes that have to our knowledge not been linked to CRC before. Two genes, SMARCB1 and STK38L, were also functionally scrutinized, providing evidence of a tumorigenic role, for SMARCB1 mutations in particular.

Exome and immune cell score analyses reveal great variation within synchronous primary colorectal cancers

Background

Approximately 4% of colorectal cancer (CRC) patients have at least two simultaneous cancers in the colon. Due to the shared environment, these synchronous CRCs (SCRCs) provide a unique setting to study colorectal carcinogenesis. Understanding whether these tumours are genetically similar or distinct is essential when designing therapeutic approaches.

Methods

We performed exome sequencing of 47 primary cancers and corresponding normal samples from 23 patients. Additionally, we carried out a comprehensive mutational signature analysis to assess whether tumours had undergone similar mutational processes and the first immune cell score analysis (IS) of SCRC to analyse the interplay between immune cell invasion and mutation profile in both lesions of an individual.

Results

The tumour pairs shared only few mutations, favouring different mutations in known CRC genes and signalling pathways and displayed variation in their signature content. Two tumour pairs had discordant mismatch repair statuses. In majority of the pairs, IS varied between primaries. Differences were not explained by any clinicopathological variable or mutation burden.

Conclusions

The study shows major diversity within SCRCs. Rather than rely on data from one tumour, our study highlights the need to evaluate both tumours of a synchronous pair for optimised targeted therapy.

Prevalence of elevated microsatellite alterations at selected tetranucleotide repeats in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) prognosis remains poor even after complete resection owing to no valuable biomarkers for recurrence and chemosensitivity. Tumors not expressing MSH3 show elevated microsatellite alterations at selected tetranucleotide repeats (EMAST). EMAST reportedly occurs in several tumors. In colorectal cancer (CRC), EMAST was reportedly correlated with 5-fluorouracil (5-FU) sensitivity. However, EMAST prevalence in PDAC and its significance as a prognostic biomarker are unknown. This study aimed to investigate EMAST prevalence in PDAC and the associations between EMAST and pathological factors, EMAST and prognosis, and EMAST and MSH3 expression via immunohistochemistry (IHC). We assessed 40 PDAC patients undergoing surgery. Genomic DNA was extracted from tumors and normal tissues. EMAST and microsatellite instability-high (MSI-H) were analyzed using five polymorphic tetranucleotide markers and five mononucleotide markers, respectively. Tumor sections were stained for MSH3, and staining intensity was evaluated via the Histoscore (H-score). Eighteen of 40 (45%) PDAC patients were EMAST-positive; however, none were MSI-H-positive. Clinicopathological characteristics including overall survival (OS) and recurrence-free survival (RFS) were not significantly different between EMAST-positive and EMAST-negative patients (P = 0.45, 0.98 respectively). IHC was performed to evaluate MSH3 protein expression levels for the PDAC tissue specimens. H-scores of EMAST-positive patients ranged from 0 to 300 (median, 40) and those of EMAST-negative patients ranged from 0 to 300 (median, 170). MSH3 protein was not significantly downregulated in EMAST-positive patients (P = 0.07). This study is a preliminary study and the number of cases investigated was small, and thus, study of a larger cohort will reveal the clinical implication of EMAST.

NHE8 Deficiency Promotes Colitis-Associated Cancer in Mice via Expansion of Lgr5-Expressing Cells

BACKGROUND & AIMS

Lgr5 overexpression has been detected in colorectal cancers (CRCs), including some cases of colitis-associated CRCs. In colitis-associated CRCs, chronic inflammation is a contributing factor in carcinogenesis. We recently reported that intestinal Na⁺/H⁺ exchanger isoform 8 (NHE8) plays an important role in intestinal mucosal protection and that loss of NHE8 expression results in an ulcerative colitis-like condition. Therefore, we hypothesized that NHE8 may be involved in the development of intestinal tumors.

METHODS

We assessed NHE8 expression in human CRCs by immunohistochemistry and studied tumor burden in NHE8 knockout (KO) mice using an azoxymethane/dextran sodium sulfate colon cancer model. We also evaluated cell proliferation in HT29^NHE8KO cells and assessed tumor growth in NOD scid gamma (NSG) mice xenografted with HT29^NHE8KO cells. To verify if a relationship exists between Lgr5 and NHE8 expression, we analyzed Lgr5 expression in NHE8KO mice by polymerase chain reaction and in situ hybridization. Lgr5 expression and cell proliferation in the absence of NHE8 were confirmed in colonic organoid cultures. The expression of β-catenin and c-Myc also were analyzed to evaluate Wnt/β-catenin activation.

RESULTS

NHE8 was undetectable in human CRC tissues. Although only 9% of NHE8 wild-type mice showed tumorigenesis in the azoxymethane/dextran sodium sulfate colon cancer model, almost 10 times more NHE8KO mice (89%) developed tumors. In the absence of NHE8, a higher colony formation unit was discovered in HT29^NHE8KO cells. In NSG mice, larger tumors developed at the site where HT29^NHE8KO cells were injected compared with HT29^{NHE8 wild type} cells. Furthermore, NHE8 deficiency resulted in increased Lgr5 expression in the colon, in HT29-derived tumors, and in colonoids. The absence of NHE8 also increased Wnt/β-catenin activation.

CONCLUSIONS

NHE8 might be an intrinsic factor that regulates Wnt/β-catenin in the intestine.

Identification of Positively and Negatively Selected Driver Gene Mutations Associated With Colorectal Cancer With Microsatellite Instability

BACKGROUND & AIMS

Recent studies have shown that cancers arise as a result of the positive selection of driver somatic events in tumor DNA, with negative selection playing only a minor role, if any. However, these investigations were concerned with alterations at nonrepetitive sequences and did not take into account mutations in repetitive sequences that have very high pathophysiological relevance in the tumors showing microsatellite instability (MSI) resulting from mismatch repair deficiency investigated in the present study.

METHODS

We performed whole-exome sequencing of 47 MSI colorectal cancers (CRCs) and confirmed results in an independent cohort of 53 MSI CRCs. We used a probabilistic model of mutational events within microsatellites, while adapting pre-existing models to analyze nonrepetitive DNA sequences. Negatively selected coding alterations in MSI CRCs were investigated for their functional and clinical impact in CRC cell lines and in a third cohort of 164 MSI CRC patients.

RESULTS

Both positive and negative selection of somatic mutations in DNA repeats was observed, leading us to identify the expected true driver genes associated with the MSI-driven tumorigenic process. Several coding negatively selected MSI-related mutational events (n = 5) were shown to have deleterious effects on tumor cells. In the tumors in which deleterious MSI mutations were observed despite the negative selection, they were associated with worse survival in MSI CRC patients (hazard ratio, 3; 95% CI, 1.1-7.9; P = .03), suggesting their anticancer impact should be offset by other as yet unknown oncogenic processes that contribute to a poor prognosis.

CONCLUSIONS

The present results identify the positive and negative driver somatic mutations acting in MSI-driven tumorigenesis, suggesting that genomic instability in MSI CRC plays a dual role in achieving tumor cell transformation. Exome sequencing data have been deposited in the European genome-phenome archive (accession: EGAS00001002477).