Canonical Wnt signaling is involved in gastric carcinogenesis. The aim of this study was to identify the link between Wnt signaling and aurora kinase A (AURKA), a target for the treatment of gastrointestinal cancers. Publicly available microarray data were used to identify phenotype-specific protein-protein interaction (PPI) subnetworks. The in silico analysis revealed a gastric cancer-specific PPI subnetwork consisting of 2745 proteins and 50,935 interactions. We focused on the link of AURKA to a Wnt-specific interaction module consisting of 92 proteins. There was a direct association of AURKA with Rac GTPase-activating protein 1 (RACGAP1), as well as with CTNBB1 (β-catenin) and CDKN1A as second-order interactors. Differential expression analysis revealed a significant downregulation of both AURKA and RACGAP1 in gastric cancer compared to noncancer controls. Biopsies from a prospective cohort of 56 patients with gastric cancer (32 intestinal type, 24 diffuse type) and 20 noncancer controls were used for validation of the identified targets. The RT-PCR data confirmed a strong correlation of AURKA and RACGAP1 gene expression both in the tumor, the tumor-adjacent and the tumor-distant mucosa. RACGAP1 in the tumor was also associated with CTNBB1 expression, and inversely associated with CDKN1A gene expression. Immunohistochemistry confirmed expression of the RACGAP1 protein in gastric cancer and the tumor-adjacent mucosa. RACGAP1 expression was not associated with tumor stage, grading, Lauren type, Helicobacter pylori infection, or age. In conclusion, AURKA is directly associated with the expression of RACGAP1, a modulator of the canonical Wnt signaling pathway.

Hepatocellular carcinoma (HCC) is a leading cause of global cancer mortality. However, little is known about the precise molecular mechanisms involved in tumor formation and pathogenesis. The primary goal of this study was to elucidate genome-wide molecular networks involved in development of HCC with multiple etiologies by exploring high quality microarray data. We undertook a comparative network analysis across 264 human microarray profiles monitoring transcript changes in healthy liver, liver cirrhosis, and HCC with viral and alcoholic etiologies. Gene co-expression profiling was used to derive a consensus gene relevance network of HCC progression that consisted of 798 genes and 2,012 links. The HCC interactome was further confirmed to be phenotype-specific and non-random. Additionally, we confirmed that co-expressed genes are more likely to share biological function, but not sub-cellular localization. Analysis of individual HCC genes revealed that they are topologically central in a human protein-protein interaction network. We used quantitative RT-PCR in a cohort of normal liver tissue (n = 8), hepatitis C virus (HCV)-induced chronic liver disease (n = 9), and HCC (n = 7) to validate co-expressions of several well-connected genes, namely ASPM, CDKN3, NEK2, RACGAP1, and TOP2A. We show that HCC is a heterogeneous disorder, underpinned by complex cross talk between immune response, cell cycle, and mRNA translation pathways. Our work provides a systems-wide resource for deeper understanding of molecular mechanisms in HCC progression and may be used further to define novel targets for efficient treatment or diagnosis of this disease.

Functional loss of TDP-43, an RNA-binding protein genetically and pathologically linked to ALS and FTD, leads to inclusion of cryptic exons in hundreds of transcripts during disease. Cryptic exons can promote degradation of affected transcripts, deleteriously altering cellular function through loss-of-function mechanisms. However, the possibility of de novo protein synthesis from cryptic exon transcripts has not been explored. Here, we show that mRNA transcripts harboring cryptic exons generate de novo proteins both in TDP-43 deficient cellular models and in disease. Using coordinated transcriptomic and proteomic studies of TDP-43 depleted iPSC-derived neurons, we identified numerous peptides that mapped to cryptic exons. Cryptic exons identified in iPSC models were highly predictive of cryptic exons expressed in brains of patients with TDP-43 proteinopathy, including cryptic transcripts that generated de novo proteins. We discovered that inclusion of cryptic peptide sequences in proteins altered their interactions with other proteins, thereby likely altering their function. Finally, we showed that these de novo peptides were present in CSF from patients with ALS. The demonstration of cryptic exon translation suggests new mechanisms for ALS pathophysiology downstream of TDP-43 dysfunction and may provide a strategy for novel biomarker development.

The neuromuscular junction (NMJ) is the peripheral synapse formed between a motor neuron axon terminal and a muscle fibre. NMJs are thought to be the primary site of peripheral pathology in many neuromuscular diseases, but innervation/denervation status is often assessed qualitatively with poor systematic criteria across studies, and separately from 3D morphological structure. Here, we describe the development of 'NMJ-Analyser', to comprehensively screen the morphology of NMJs and their corresponding innervation status automatically. NMJ-Analyser generates 29 biologically relevant features to quantitatively define healthy and aberrant neuromuscular synapses and applies machine learning to diagnose NMJ degeneration. We validated this framework in longitudinal analyses of wildtype mice, as well as in four different neuromuscular disease models: three for amyotrophic lateral sclerosis (ALS) and one for peripheral neuropathy. We showed that structural changes at the NMJ initially occur in the nerve terminal of mutant TDP43 and FUS ALS models. Using a machine learning algorithm, healthy and aberrant neuromuscular synapses are identified with 95% accuracy, with 88% sensitivity and 97% specificity. Our results validate NMJ-Analyser as a robust platform for systematic and structural screening of NMJs, and pave the way for transferrable, and cross-comparison and high-throughput studies in neuromuscular diseases.

Neoadjuvant therapy followed by surgery is the standard of care for locally advanced esophageal adenocarcinoma (EAC). Unfortunately, response to neoadjuvant chemotherapy (NAC) is poor (20-37%), as is the overall survival benefit at five years (9%). The EAC genome is complex and heterogeneous between patients, and it is not yet understood whether specific mutational patterns may result in chemotherapy sensitivity or resistance. To identify associations between genomic events and response to NAC in EAC, a comparative genomic analysis was performed in 65 patients with extensive clinical and pathological annotation using whole-genome sequencing (WGS). We defined response using Mandard Tumor Regression Grade (TRG), with responders classified as TRG1-2 (n = 27) and non-responders classified as TRG4-5 (n =38). We report a higher non-synonymous mutation burden in responders (median 2.08/Mb vs. 1.70/Mb, p = 0.036) and elevated copy number variation in non-responders (282 vs. 136/patient, p < 0.001). We identified copy number variants unique to each group in our cohort, with cell cycle (CDKN2A, CCND1), c-Myc (MYC), RTK/PIK3 (KRAS, EGFR) and gastrointestinal differentiation (GATA6) pathway genes being specifically altered in non-responders. Of note, NAV3 mutations were exclusively present in the non-responder group with a frequency of 22%. Thus, lower mutation burden, higher chromosomal instability and specific copy number alterations are associated with resistance to NAC.

Esophageal adenocarcinoma (EAC) incidence is increasing while 5-year survival rates remain less than 15%. A lack of experimental models has hampered progress. We have generated clinically annotated EAC organoid cultures that recapitulate the morphology, genomic, and transcriptomic landscape of the primary tumor including point mutations, copy number alterations, and mutational signatures. Karyotyping of organoid cultures has confirmed polyclonality reflecting the clonal architecture of the primary tumor. Furthermore, subclones underwent clonal selection associated with driver gene status. Medium throughput drug sensitivity testing demonstrates the potential of targeting receptor tyrosine kinases and downstream mediators. EAC organoid cultures provide a pre-clinical tool for studies of clonal evolution and precision therapeutics.

Cancers occurring at the gastroesophageal junction (GEJ) are classified as predominantly esophageal or gastric, which is often difficult to decipher. We hypothesized that the transcriptomic profile might reveal molecular subgroups which could help to define the tumor origin and behavior beyond anatomical location. The gene expression profiles of 107 treatment-naïve, intestinal type, gastroesophageal adenocarcinomas were assessed by the Illumina-HTv4.0 beadchip. Differential gene expression (limma), unsupervised subgroup assignment (mclust) and pathway analysis (gage) were undertaken in R statistical computing and results were related to demographic and clinical parameters. Unsupervised assignment of the gene expression profiles revealed three distinct molecular subgroups, which were not associated with anatomical location, tumor stage or grade (p > 0.05). Group 1 was enriched for pathways involved in cell turnover, Group 2 was enriched for metabolic processes and Group 3 for immune-response pathways. Patients in group 1 showed the worst overall survival (p = 0.019). Key genes for the three subtypes were confirmed by immunohistochemistry. The newly defined intrinsic subtypes were analyzed in four independent datasets of gastric and esophageal adenocarcinomas with transcriptomic data available (RNAseq data: OCCAMS cohort, n = 158; gene expression arrays: Belfast, n = 63; Singapore, n = 191; Asian Cancer Research Group, n = 300). The subgroups were represented in the independent cohorts and pooled analysis confirmed the prognostic effect of the new subtypes. In conclusion, adenocarcinomas at the GEJ comprise three distinct molecular phenotypes which do not reflect anatomical location but rather inform our understanding of the key pathways expressed.

Variants of UNC13A, a critical gene for synapse function, increase the risk of amyotrophic lateral sclerosis and frontotemporal dementia1-3, two related neurodegenerative diseases defined by mislocalization of the RNA-binding protein TDP-434,5. Here we show that TDP-43 depletion induces robust inclusion of a cryptic exon in UNC13A, resulting in nonsense-mediated decay and loss of UNC13A protein. Two common intronic UNC13A polymorphisms strongly associated with amyotrophic lateral sclerosis and frontotemporal dementia risk overlap with TDP-43 binding sites. These polymorphisms potentiate cryptic exon inclusion, both in cultured cells and in brains and spinal cords from patients with these conditions. Our findings, which demonstrate a genetic link between loss of nuclear TDP-43 function and disease, reveal the mechanism by which UNC13A variants exacerbate the effects of decreased TDP-43 function. They further provide a promising therapeutic target for TDP-43 proteinopathies.

The epithelial to mesenchymal transition (EMT) is a key cellular process underlying cancer progression, with multiple intermediate states whose molecular hallmarks remain poorly characterised. To fill this gap, we present a method to robustly evaluate EMT transformation in individual tumours based on transcriptomic signals. We apply this approach to explore EMT trajectories in 7180 tumours of epithelial origin and identify three macro-states with prognostic and therapeutic value, attributable to epithelial, hybrid E/M and mesenchymal phenotypes. We show that the hybrid state is relatively stable and linked with increased aneuploidy. We further employ spatial transcriptomics and single cell datasets to explore the spatial heterogeneity of EMT transformation and distinct interaction patterns with cytotoxic, NK cells and fibroblasts in the tumour microenvironment. Additionally, we provide a catalogue of genomic events underlying distinct evolutionary constraints on EMT transformation. This study sheds light on the aetiology of distinct stages along the EMT trajectory, and highlights broader genomic and environmental hallmarks shaping the mesenchymal transformation of primary tumours.

Head and neck squamous cell carcinoma (HNSCC) is a molecularly and spatially heterogeneous disease frequently characterized by impairment of immunosurveillance mechanisms. Despite recent success with immunotherapy treatment, disease progression still occurs quickly after treatment in the majority of cases, suggesting the need to improve patient selection strategies. In the quest for biomarkers that may help inform response to checkpoint blockade, we characterized the tumor microenvironment (TME) of 162 HNSCC primary tumors of diverse etiologic and spatial origin, through gene expression and IHC profiling of relevant immune proteins, T-cell receptor (TCR) repertoire analysis, and whole-exome sequencing. We identified five HNSCC TME categories based on immune/stromal composition: (i) cytotoxic, (ii) plasma cell rich, (iii) dendritic cell rich, (iv) macrophage rich, and (v) immune-excluded. Remarkably, the cytotoxic and plasma cell rich subgroups exhibited a phenotype similar to tertiary lymphoid structures (TLS), which have been previously linked to immunotherapy response. We also found an increased richness of the TCR repertoire in these two subgroups and in never smokers. Mutational patterns evidencing APOBEC activity were enriched in the plasma cell high subgroup. Furthermore, specific signal propagation patterns within the Ras/ERK and PI3K/AKT pathways associated with distinct immune phenotypes. While traditionally CD8/CD3 T-cell infiltration and immune checkpoint expression (e.g., PD-L1) have been used in the patient selection process for checkpoint blockade treatment, we suggest that additional biomarkers, such as TCR productive clonality, smoking history, and TLS index, may have the ability to pull out potential responders to benefit from immunotherapeutic agents.

Esophageal adenocarcinoma (EAC) has a poor outcome, and targeted therapy trials have thus far been disappointing owing to a lack of robust stratification methods. Whole-genome sequencing (WGS) analysis of 129 cases demonstrated that this is a heterogeneous cancer dominated by copy number alterations with frequent large-scale rearrangements. Co-amplification of receptor tyrosine kinases (RTKs) and/or downstream mitogenic activation is almost ubiquitous; thus tailored combination RTK inhibitor (RTKi) therapy might be required, as we demonstrate in vitro. However, mutational signatures showed three distinct molecular subtypes with potential therapeutic relevance, which we verified in an independent cohort (n = 87): (i) enrichment for BRCA signature with prevalent defects in the homologous recombination pathway; (ii) dominant T>G mutational pattern associated with a high mutational load and neoantigen burden; and (iii) C>A/T mutational pattern with evidence of an aging imprint. These subtypes could be ascertained using a clinically applicable sequencing strategy (low coverage) as a basis for therapy selection.

Lapatinib, a dual EGFR and HER2 inhibitor has shown disappointing results in clinical trials of metastatic oesophago-gastric adenocarcinomas (OGAs), and in vitro studies suggest that MET, IGFR, and HER3 confer resistance. This trial applied Lapatinib in the curative neoadjuvant setting and investigated the feasibility and utility of additional endoscopy and biopsy for assessment of resistance mechanisms ex vivo and in vivo.

Current strategies to guide selection of neoadjuvant therapy in oesophageal adenocarcinoma (OAC) are inadequate. We assessed the ability of a DNA damage immune response (DDIR) assay to predict response following neoadjuvant chemotherapy in OAC.

Transcriptomics is a developing field with new methods of analysis being produced which may hold advantages in price, accuracy, or information output. QuantSeq is a form of 3' sequencing produced by Lexogen which aims to obtain similar gene-expression information to RNA-seq with significantly fewer reads, and therefore at a lower cost. QuantSeq is also able to provide information on differential polyadenylation. We applied both QuantSeq at low read depth and total RNA-seq to the same two sets of mouse spinal cord RNAs, each comprised by four controls and four mutants related to the neurodegenerative disease amyotrophic lateral sclerosis. We found substantial differences in which genes were found to be significantly differentially expressed by the two methods. Some of this difference likely due to the difference in number of reads between our QuantSeq and RNA-seq data. Other sources of difference can be explained by the differences in the way the two methods handle genes with different primary transcript lengths and how likely each method is to find a gene to be differentially expressed at different levels of overall gene expression. This work highlights how different methods aiming to assess expression difference can lead to different results.

The poor outcomes in esophageal adenocarcinoma (EAC) prompted us to interrogate the pattern and timing of metastatic spread. Whole-genome sequencing and phylogenetic analysis of 388 samples across 18 individuals with EAC showed, in 90% of patients, that multiple subclones from the primary tumor spread very rapidly from the primary site to form multiple metastases, including lymph nodes and distant tissues-a mode of dissemination that we term 'clonal diaspora'. Metastatic subclones at autopsy were present in tissue and blood samples from earlier time points. These findings have implications for our understanding and clinical evaluation of EAC.

Many lines of evidence have highlighted the role played by heterogeneous nuclear ribonucleoproteins in amyotrophic lateral sclerosis. In this study, we have aimed to identify transcripts co-regulated by TAR DNA-binding protein 43 kDa and highly conserved heterogeneous nuclear ribonucleoproteins which have been previously shown to regulate TAR DNA-binding protein 43 kDa toxicity (deleted in azoospermia-associated protein 1, heterogeneous nuclear ribonucleoprotein -Q, -D, -K and -U). Using the transcriptome analyses, we have uncovered that Nitric Oxide Synthase 1 Adaptor Protein mRNA is a direct TAR DNA-binding protein 43 kDa target, and in flies, its modulation alone can rescue TAR DNA-binding protein 43 kDa pathology. In primary mouse cortical neurons, we show that TAR DNA-binding protein 43 kDa mediated downregulation of Nitric Oxide Synthase 1 Adaptor Protein expression strongly affects the NMDA-receptor signalling pathway. In human patients, the downregulation of Nitric Oxide Synthase 1 Adaptor Protein mRNA strongly correlates with TAR DNA-binding protein 43 kDa proteinopathy as measured by cryptic Stathmin-2 and Unc-13 homolog A cryptic exon inclusion. Overall, our results demonstrate that Nitric Oxide Synthase 1 Adaptor Protein may represent a novel disease-relevant gene, potentially suitable for the development of new therapeutic strategies.

A variety of mutational processes drive cancer development, but their dynamics across the entire disease spectrum from pre-cancerous to advanced neoplasia are poorly understood. We explore the mutagenic processes shaping oesophageal adenocarcinoma tumorigenesis in 997 instances comprising distinct stages of this malignancy, from Barrett Oesophagus to primary tumours and advanced metastatic disease. The mutational landscape is dominated by the C[T > C/G]T substitution enriched signatures SBS17a/b, which are linked with TP53 mutations, increased proliferation, genomic instability and disease progression. The APOBEC mutagenesis signature is a weak but persistent signal amplified in primary tumours. We also identify prevalent alterations in DNA damage repair pathways, with homologous recombination, base and nucleotide excision repair and translesion synthesis mutated in up to 50% of the cohort, and surprisingly uncoupled from transcriptional activity. Among these, the presence of base excision repair deficiencies show remarkably poor prognosis in the cohort. In this work, we provide insights on the mutational aetiology and changes enabling the transition from pre-neoplastic to advanced oesophageal adenocarcinoma.

Tumor mass dormancy is the key intermediate step between immune surveillance and cancer progression, yet due to its transitory nature it has been difficult to capture and characterize. Little is understood of its prevalence across cancer types and of the mutational background that may favor such a state. While this balance is finely tuned internally by the equilibrium between cell proliferation and cell death, the main external factors contributing to tumor mass dormancy are immunological and angiogenic. To understand the genomic and cellular context in which tumor mass dormancy may develop, we comprehensively profiled signals of immune and angiogenic dormancy in 9,631 cancers from the Cancer Genome Atlas and linked them to tumor mutagenesis. We find evidence for immunological and angiogenic dormancy-like signals in 16.5% of bulk sequenced tumors, with a frequency of up to 33% in certain tissues. Mutations in the CASP8 and HRAS oncogenes were positively selected in dormant tumors, suggesting an evolutionary pressure for controlling cell growth/apoptosis signals. By surveying the mutational damage patterns left in the genome by known cancer risk factors, we found that aging-induced mutations were relatively depleted in these tumors, while patterns of smoking and defective base excision repair were linked with increased tumor mass dormancy. Furthermore, we identified a link between APOBEC mutagenesis and dormancy, which comes in conjunction with immune exhaustion and may partly depend on the expression of the angiogenesis regulator PLG as well as interferon and chemokine signals. Tumor mass dormancy also appeared to be impaired in hypoxic conditions in the majority of cancers. The microenvironment of dormant cancers was enriched in cytotoxic and regulatory T cells, as expected, but also in macrophages and showed a reduction in inflammatory Th17 signals. Finally, tumor mass dormancy was linked with improved patient survival outcomes. Our analysis sheds light onto the complex interplay between dormancy, exhaustion, APOBEC activity and hypoxia, and sets directions for future mechanistic explorations.

Nuclear depletion and cytoplasmic aggregation of the RNA-binding protein TDP-43 is the hallmark of ALS, occurring in over 97% of cases. A key consequence of TDP-43 nuclear loss is the de-repression of cryptic exons. Whilst TDP-43 regulated cryptic splicing is increasingly well catalogued, cryptic alternative polyadenylation (APA) events, which define the 3' end of last exons, have been largely overlooked, especially when not associated with novel upstream splice junctions. We developed a novel bioinformatic approach to reliably identify distinct APA event types: alternative last exons (ALE), 3'UTR extensions (3'Ext) and intronic polyadenylation (IPA) events. We identified novel neuronal cryptic APA sites induced by TDP-43 loss of function by systematically applying our pipeline to a compendium of publicly available and in house datasets. We find that TDP-43 binding sites and target motifs are enriched at these cryptic events and that TDP-43 can have both repressive and enhancing action on APA. Importantly, all categories of cryptic APA can also be identified in ALS and FTD post mortem brain regions with TDP-43 proteinopathy underlining their potential disease relevance. RNA-seq and Ribo-seq analyses indicate that distinct cryptic APA categories have different downstream effects on transcript and translation. Intriguingly, cryptic 3'Exts occur in multiple transcription factors, such as ELK1, SIX3, and TLX1, and lead to an increase in wild-type protein levels and function. Finally, we show that an increase in RNA stability leading to a higher cytoplasmic localisation underlies these observations. In summary, we demonstrate that TDP-43 nuclear depletion induces a novel category of cryptic RNA processing events and we expand the palette of TDP-43 loss consequences by showing this can also lead to an increase in normal protein translation.

Barrett's esophagus (BE) involves a metaplastic replacement of native esophageal squamous epithelium (Sq) by columnar-intestinalized mucosa, and it is the main risk factor for Barrett-related adenocarcinoma (BAc). Ultra-conserved regions (UCRs) are a class non-coding sequences that are conserved in humans, mice and rats. More than 90% of UCRs are transcribed (T-UCRs) in normal tissues, and are altered at transcriptional level in tumorigenesis. To identify the T-UCR profiles that are dysregulated in Barrett's mucosa transformation, microarray analysis was performed on a discovery set of 51 macro-dissected samples obtained from 14 long-segment BE patients. Results were validated in an independent series of esophageal biopsy/surgery specimens and in two murine models of Barrett's esophagus (i.e. esophagogastric-duodenal anastomosis). Progression from normal to BE to adenocarcinoma was each associated with specific and mutually exclusive T-UCR signatures that included up-regulation of uc.58-, uc.202-, uc.207-, and uc.223- and down-regulation of uc.214+. A 9 T-UCR signature characterized BE versus Sq (with the down-regulation of uc.161-, uc.165-, and uc.327-, and the up-regulation of uc.153-, uc.158-, uc.206-, uc.274-, uc.472-, and uc.473-). Analogous BE-specific T-UCR profiles were shared by human and murine lesions. This study is the first demonstration of a role for T-UCRs in the transformation of Barrett's mucosa.