Cross-species genomics identifies DLG2 as a tumor suppressor in osteosarcoma

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.

SIGNIFICANCE

The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Decoding the Impact of Tumor Microenvironment in Osteosarcoma Progression and Metastasis

Osteosarcoma (OS) is a heterogeneous, highly metastatic bone malignancy in children and adolescents. Despite advancements in multimodal treatment strategies, the prognosis for patients with metastatic or recurrent disease has not improved significantly in the last four decades. OS is a highly heterogeneous tumor; its genetic background and the mechanism of oncogenesis are not well defined. Unfortunately, no effective molecular targeted therapy is currently available for this disease. Understanding osteosarcoma's tumor microenvironment (TME) has recently gained much interest among scientists hoping to provide valuable insights into tumor heterogeneity, progression, metastasis, and the identification of novel therapeutic avenues. Here, we review the current understanding of the TME of OS, including different cellular and noncellular components, their crosstalk with OS tumor cells, and their involvement in tumor progression and metastasis. We also highlight past/current clinical trials targeting the TME of OS for effective therapies and potential future therapeutic strategies with negligible adverse effects.

Large Scale Comparative Deconvolution Analysis of the Canine and Human Osteosarcoma Tumor Microenvironment Uncovers Conserved Clinically Relevant Subtypes

Osteosarcoma is a relatively rare but aggressive cancer of the bones with a shortage of effective biomarkers. Although less common in humans, Osteosarcomas are fairly common in adult pet dogs and have been shown to share many similarities with their human analogs. In this work, we analyze bulk transcriptomic data of 213 primary and 100 metastatic Osteosarcoma samples from 210 pet dogs enrolled in nation-wide clinical trials to uncover three Tumor Microenvironment (TME)-based subtypes: Immune Enriched (IE), Immune Enriched Dense Extra-Cellular Matrix-like (IE-ECM) and Immune Desert (ID) with distinct cell type compositions, oncogenic pathway activity and chromosomal instability. Furthermore, leveraging bulk transcriptomic data of canine primary tumors and their matched metastases from different sites, we characterize how the Osteosarcoma TME evolves from primary to metastatic disease in a standard of care clinical setting and assess its overall impact on clinical outcomes of canines. Most importantly, we find that TME-based subtypes of canine Osteosarcomas are conserved in humans and predictive of progression free survival outcomes of human patients, independently of known prognostic biomarkers such as presence of metastatic disease at diagnosis and percent necrosis following chemotherapy. In summary, these results demonstrate the power of using canines to model the human Osteosarcoma TME and discover novel biomarkers for clinical translation.

Single-molecule approaches for DNA damage detection and repair: A focus on Repair Assisted Damage Detection (RADD)

The human genome is continually exposed to various stressors, which can result in DNA damage, mutations, and diseases. Among the different types of DNA damage, single-strand lesions are commonly induced by external stressors and metabolic processes. Accurate detection and quantification of DNA damage are crucial for understanding repair mechanisms, assessing environmental impacts, and evaluating response to therapy. However, traditional techniques have limitations in sensitivity and the ability to detect multiple types of damage. In recent years, single-molecule fluorescence approaches have emerged as powerful tools for precisely localizing and quantifying DNA damage. Repair Assisted Damage Detection (RADD) is a single-molecule technique that employs specific repair enzymes to excise damaged bases and incorporates fluorescently labeled nucleotides to visualize the damage. This technique provides valuable insights into repair efficiency and sequence-specific damage. In this review, we discuss the principles and applications of RADD assays, highlighting their potential for enhancing our understanding of DNA damage and repair processes.

Use of tumor suppressor genes of naked mole rats for human cancer treatment

Cancer not only has a significant prevalence in the human population, but is also leading cause of death in animals. Despite a long history, our battle against cancer continues. Cross-species comparative genomics offers insight into shared genes and pathways by analyzing genomic data across species, enhancing our understanding of cancer mechanisms, evolutionary processes, and possible therapeutic targets. However, no previous study has demonstrated the inhibitory effects of tumor suppressor genes from one species on tumor cells from another. The naked mole rat is the only mammal yet to be found with cancer that is attributed to its tumor suppressor genes. In this study, we constructed phylogenetic trees and assessed the anti-tumor activity of two suppressor genes, programmed cell death molecule 5 (PDCD5) and dickkopf 3 (DKK3), from rats, mice, and humans. DKK3 robustly inhibited the proliferation of breast cancer cells within and across species due to its highly conserved protein sequence. However, the cross-species inhibitory effect of PDCD5 on breast cancer cells was inconsistent due to significant sequence variations. Intriguingly, PDCD5 from the naked mole rat demonstrated potent anti-tumor activity against breast cancer cells from mice, rats, and humans, surpassing that of PDCD5 from parental species. Our results demonstrate that the suppressor genes from the naked mole rat have a strong inhibitory effect on cancer cells, indicating that the powerful anti-cancer functions of the naked mole rat may be useful for human tumor treatment.

Distinct mechanisms of PTEN inactivation in dogs and humans highlight convergent molecular events that drive cell division in the pathogenesis of osteosarcoma

A hallmark of osteosarcoma in both human and canine tumors is somatic fragmentation and rearrangement of chromosome structure which leads to recurrent increases and decreases in DNA copy number. The PTEN gene has been implicated as an important tumor suppressor in osteosarcoma via forward genetic screens. Here, we analyzed copy number changes, promoter methylation and transcriptomes to better understand the role of PTEN in canine and human osteosarcoma. Reduction in PTEN copy number was observed in 23 of 95 (25%) of the canine tumors examined leading to corresponding decreases in PTEN transcript levels from RNA-Seq samples. Unexpectedly, canine tumors with an intact PTEN locus had higher levels of PTEN transcripts than human tumors. This variation in transcript abundance was used to evaluate the role of PTEN in osteosarcoma biology. Decreased PTEN copy number and transcript level was observed in - and likely an important driver of - increases in cell cycle transcripts in four independent canine transcriptional datasets. In human osteosarcoma, homozygous copy number loss was not observed, instead increased methylation of the PTEN promoter was associated with increased cell cycle transcripts. Somatic modification of PTEN, either by homozygous deletion in dogs or by promoter methylation in humans, is clinically relevant to osteosarcoma, because the cell cycle related transcripts are associated with patient outcomes. The PTEN gene is part of a syntenic rearrangement unique to the canine genome, making it susceptible to somatic loss of both copies of distal chromosome 26 which also includes the FAS death receptor. SIGNIFICANCE STATEMENT: PTEN function is abrogated by different mechanisms in canine and human osteosarcoma tumors leading to uncontrolled cell cycling. Somatic loss of this canine specific syntenic region may help explain why the canine genome appears to be uniquely susceptible to osteosarcoma. Syntenic arrangement, in the context of copy number change, may lead to synergistic interactions that in turn modify species specific cancer risk. Comparative models of tumorigenesis may utilize different driver mechanisms.

Deciphering the Signaling Mechanisms of Osteosarcoma Tumorigenesis

Osteosarcoma (OS) is the predominant primary bone tumor in the pediatric and adolescent populations. It has high metastatic potential, with the lungs being the most common site of metastasis. In contrast to many other sarcomas, OS lacks conserved translocations or genetic mutations; instead, it has heterogeneous abnormalities, including somatic DNA copy number alteration, ploidy, chromosomal amplification, and chromosomal loss and gain. Unfortunately, clinical outcomes have not significantly improved in over 30 years. Currently, no effective molecularly targeted therapies are available for this disease. Several genomic studies showed inactivation in the tumor suppressor genes, including p53, RB, and ATRX, and hyperactivation of the tumor promoter genes, including MYC and MDM2, in OS. Alterations in the major signaling pathways, including the PI3K/AKT/mTOR, JAK/STAT, Wnt/β-catenin, NOTCH, Hedgehog/Gli, TGF-β, RTKs, RANK/RANKL, and NF-κB signaling pathways, have been identified in OS development and metastasis. Although OS treatment is currently based on surgical excision and systematic multiagent therapies, several potential targeted therapies are in development. This review focuses on the major signaling pathways of OS, and we propose a biological rationale to consider novel and targeted therapies in the future.

Over-expression of ADAR1 in mice does not initiate or accelerate cancer formation in vivo

Adenosine to inosine editing (A-to-I) in regions of double stranded RNA (dsRNA) is mediated by adenosine deaminase acting on RNA 1 (ADAR1) or ADAR2. ADAR1 and A-to-I editing levels are increased in many human cancers. Inhibition of ADAR1 has emerged as a high priority oncology target, however, whether ADAR1 overexpression enables cancer initiation or progression has not been directly tested. We established a series of in vivo models to allow overexpression of full-length ADAR1, or its individual isoforms, to test if increased ADAR1 expression was oncogenic. Widespread over-expression of ADAR1 or the p110 or p150 isoforms individually as sole lesions was well tolerated and did not result in cancer initiation. Therefore, ADAR1 overexpression alone is not sufficient to initiate cancer. We demonstrate that endogenous ADAR1 and A-to-I editing increased upon immortalization in murine cells, consistent with the observations from human cancers. We tested if ADAR1 over-expression could co-operate with cancer initiated by loss of tumour suppressors using a model of osteosarcoma. We did not see a disease potentiating or modifying effect of overexpressing ADAR1 or its isoforms in the models assessed. We conclude that increased ADAR1 expression and A-to-I editing in cancers is most likely a consequence of tumor formation.

Cellular dynamics of distinct skeletal cells and the development of osteosarcoma

Bone contributes to the maintenance of vital biological activities. At the cellular level, multiple types of skeletal cells, including skeletal stem and progenitor cells (SSPCs), osteoblasts, chondrocytes, marrow stromal cells, and adipocytes, orchestrate skeletal events such as development, aging, regeneration, and tumorigenesis. Osteosarcoma (OS) is a primary malignant tumor and the main form of bone cancer. Although it has been proposed that the cellular origins of OS are in osteogenesis-related skeletal lineage cells with cancer suppressor gene mutations, its origins have not yet been fully elucidated because of a poor understanding of whole skeletal cell diversity and dynamics. Over the past decade, the advent and development of single-cell RNA sequencing analyses and mouse lineage-tracing approaches have revealed the diversity of skeletal stem and its lineage cells. Skeletal stem cells (SSCs) in the bone marrow endoskeletal region have now been found to efficiently generate OS and to be robust cells of origin under p53 deletion conditions. The identification of SSCs may lead to a more limited redefinition of bone marrow mesenchymal stem/stromal cells (BM-MSCs), and this population has been thought to contain cells from which OS originates. In this mini-review, we discuss the cellular diversity and dynamics of multiple skeletal cell types and the origin of OS in the native in vivo environment in mice. We also discuss future challenges in the study of skeletal cells and OS.

Bridging clinic and wildlife care with AI-powered pan-species computational pathology

Cancers occur across species. Understanding what is consistent and varies across species can provide new insights into cancer initiation and evolution, with significant implications for animal welfare and wildlife conservation. We build a pan-species cancer digital pathology atlas (panspecies.ai) and conduct a pan-species study of computational comparative pathology using a supervised convolutional neural network algorithm trained on human samples. The artificial intelligence algorithm achieves high accuracy in measuring immune response through single-cell classification for two transmissible cancers (canine transmissible venereal tumour, 0.94; Tasmanian devil facial tumour disease, 0.88). In 18 other vertebrate species (mammalia = 11, reptilia = 4, aves = 2, and amphibia = 1), accuracy (range 0.57-0.94) is influenced by cell morphological similarity preserved across different taxonomic groups, tumour sites, and variations in the immune compartment. Furthermore, a spatial immune score based on artificial intelligence and spatial statistics is associated with prognosis in canine melanoma and prostate tumours. A metric, named morphospace overlap, is developed to guide veterinary pathologists towards rational deployment of this technology on new samples. This study provides the foundation and guidelines for transferring artificial intelligence technologies to veterinary pathology based on understanding of morphological conservation, which could vastly accelerate developments in veterinary medicine and comparative oncology.

Estimates of Genomic Heritability and the Marker-Derived Gene for Re(Production) Traits in Xinggao Sheep

Xinggao sheep are a breed of Chinese domestic sheep that are adapted to the extremely cold climatic features of the Hinggan League in China. The economically vital reproductive trait of ewes (litter size, LS) and productive traits of lambs (birth weight, BWT; weaning weight, WWT; and average daily gain, ADG) are expressed in females and later in life after most of the selection decisions have been made. This study estimated the genetic parameters for four traits to explore the genetic mechanisms underlying the variation, and we performed genome-wide association study (GWAS) tests on a small sample size to identify novel marker trait associations (MTAs) associated with prolificacy and growth. We detected two suggestive significant single-nucleotide polymorphisms (SNPs) associated with LS and eight significant SNPs for BWT, WWT, and ADG. These candidate loci and genes also provide valuable information for further fine-mapping of QTLs and improvement of reproductive and productive traits in sheep.

EVI2B Is a Prognostic Biomarker and Is Correlated with Monocyte and Macrophage Infiltration in Osteosarcoma Based on an Integrative Analysis

Osteosarcoma (OS) is the most common malignant bone tumor. However, treatment strategies have not changed over the past 30 years. The relationship between OS and the immune microenvironment may provide a basis for the establishment of novel therapeutic targets. In this study, a large-scale gene expression dataset (GSE42352) was used to identify key genes in OS. A Target-OS dataset from the Cancer Genome Atlas was used as a validation set. Ecotropic viral integration site 2B (EVI2B) was significantly upregulated in OS tumor samples. Differentially expressed genes (DEGs) were identified between samples with high and low EVI2B expression in both the test and validation cohorts. The top three functions of DEGs determined by a gene set enrichment analysis (GSEA) were chemokine signaling, cytokine-cytokine receptor interaction, and Human T-cell leukemia virus 1 infection. A prognostic prediction model including EVI2B, DOCK2, and CD33 was constructed by a Cox regression analysis. This model indicated that EVI2B is an independent protective prognostic marker in OS. An analysis of immune infiltration further showed that high EVI2B expression levels were correlated with high levels of macrophage infiltration. Protein expression data derived from the Human Protein Atlas suggested EVI2B to be highly expressed in monocytes. Finally, we validated the elevated expression of EVI2B in OS cell lines and OS tissue samples; these results were consistent with those of the analyses of the GSE42352 and Target-OS datasets. Our integrative bioinformatics analysis and experimental results provide clear evidence for the prognostic value of EVI2B in OS and its close relationship with monocyte and macrophage infiltration.

The association of gene polymorphisms with milk production and mastitis resistance phenotypic traits in dairy cattle

The aim of this study was to evaluate the association between gene polymorphisms (SNPs) and mastitis indicators and their relationship with milk production profitability in dairy herd.A functional analysis was also performed of five genes containing the studied SNPs and those located close by. DNA was isolated from the hair bulb of 320 dairy cows kept in three herds and SNP-microarray analysis was performed. The data on 299 cows was subjected to final statistical analysis using AI-REML method with one-trait repeatability test-day animal model and pedigree information using the DMU4 package. Five from 35 SNPs significantly associated with mastitis indicators or production traits and located within a gene or no more than 500,000 nucleotides from the gene were selected for the functional and economic analysis. A questionnaire was also developed to collect associated economic data of 219 cows from three herds, such as the value of milk production and direct costs incurred over three years; this allowed the gross margin, direct profitability index and direct costs incurred to produce one liter of milk to be determined, among others. None of the five studied SNPs were related to protein content. The rs110785912(T/A), found near CXCR4, and rs136813430(T/C), located in the TLR4 gene exon, were associated with lnSCC, while rs110455063(C/G), located near IGFI, was associated with milk yield, fat and total solid contents. rs109421300(T/C), associated with fat/protein content ratio, as well as fat and total solid content, is located in the DGAT1 gene intron. rs41587003(A/C), located in the DLG2 gene intron, was associated with lactose content. The economic analysis revealed differences between the variants of the three tested SNPs. The T/C variant of the rs136813430(T/C) SNP was characterized by the highest gross margin, the highest direct profitability index and the lowest costs incurred to produce 1 liter of milk. The T/A variant of rs110785912(T/A) was related to low lnSCC and was characterized by the highest direct profitability index. In turn, the C/C variant of the rs41587003(T/C) was related to the lowest level of lactose and the highest costs of milk production. It appears that rs136813430(T/C) may be the most promising of the tested SNPs for increasing the profitability of milk production. To our knowledge, it is the first effort to assess directly a correlation between the DNA polymorphism and economic output of a dairy enterprise.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Non-coding RNAs as potential biomarkers in osteosarcoma

Osteosarcoma (OS) is a primary solid malignant tumor that occurs most frequently in the metaphysis of long bones. More likely to happen to children and adolescents. OS has high mortality and disability rate. However, the etiology and pathogenesis of OS have not been fully understood till now. Due to the lack of effective biomarkers, OS cannot be precisely detected in the early stage. With the application of next-generation and high-throughput sequencing, more and more abnormally expressed non-coding RNAs(ncRNAs) have been identified in OS. Growing evidences have suggested the ncRNAs, such as microRNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), have played an important role in the tumorigenesis and progression of OS. Thus, they can be served as novel biomarkers for diagnosis, treatment and prognosis. This review summarized the application of ncRNA as biomarkers in OS in detail, and discussed the limitation and future improvement of the potential biomarkers.

Comparative Evaluation of Tumor-Infiltrating Lymphocytes in Companion Animals: Immuno-Oncology as a Relevant Translational Model for Cancer Therapy

Despite the important role of preclinical experiments to characterize tumor biology and molecular pathways, there are ongoing challenges to model the tumor microenvironment, specifically the dynamic interactions between tumor cells and immune infiltrates. Comprehensive models of host-tumor immune interactions will enhance the development of emerging treatment strategies, such as immunotherapies. Although in vitro and murine models are important for the early modelling of cancer and treatment-response mechanisms, comparative research studies involving veterinary oncology may bridge the translational pathway to human studies. The natural progression of several malignancies in animals exhibits similar pathogenesis to human cancers, and previous studies have shown a relevant and evaluable immune system. Veterinary oncologists working alongside oncologists and cancer researchers have the potential to advance discovery. Understanding the host-tumor-immune interactions can accelerate drug and biomarker discovery in a clinically relevant setting. This review presents discoveries in comparative immuno-oncology and implications to cancer therapy.

Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma

Metastatic cutaneous squamous cell carcinoma (CSCC) is a highly morbid disease requiring radical surgery and adjuvant therapy, which is associated with a poor prognosis. Yet, compared to other advanced malignancies, relatively little is known of the genomic landscape of metastatic CSCC. We have previously reported the mutational signatures and mutational patterns of CCCTC-binding factor (CTCF) regions in metastatic CSCC. However, many other genomic components (indel signatures, non-coding drivers, and structural variants) of metastatic CSCC have not been reported. To this end, we performed whole genome sequencing on lymph node metastases and blood DNA from 25 CSCC patients with regional metastases of the head and neck. We designed a multifaceted computational analysis at the whole genome level to provide a more comprehensive perspective of the genomic landscape of metastatic CSCC. In the non-coding genome, 3′ untranslated region (3′UTR) regions of EVC (48% of specimens), PPP1R1A (48% of specimens), and ABCA4 (20% of specimens) along with the tumor-suppressing long non-coding RNA (lncRNA) LINC01003 (64% of specimens) were significantly functionally altered (Q-value < 0.05) and represent potential non-coding biomarkers of CSCC. Recurrent copy number loss in the tumor suppressor gene PTPRD was observed. Gene amplification was much less frequent, and few genes were recurrently amplified. Single nucleotide variants driver analyses from three tools confirmed TP53 and CDKN2A as recurrently mutated genes but also identified C9 as a potential novel driver in this disease. Furthermore, indel signature analysis highlighted the dominance of ID signature 13 (ID13) followed by ID8 and ID9. ID9 has previously been shown to have no association with skin melanoma, unlike ID13 and ID8, suggesting a novel pattern of indel variation in metastatic CSCC. The enrichment analysis of various genetically altered candidates shows enrichment of “TGF-beta regulation of extracellular matrix” and “cell cycle G1 to S check points.” These enriched terms are associated with genetic instability, cell proliferation, and migration as mechanisms of genomic drivers of metastatic CSCC.

Origin and Therapies of Osteosarcoma

Simple Summary

Osteosarcoma is the most common malignant bone tumor in children, with a 5-year survival rate ranging from 70% to 20% depending on the aggressiveness of the disease. The current treatments have not evolved over the past four decades due in part to the genetic complexity of the disease and its heterogeneity. This review will summarize the current knowledge of OS origin, diagnosis and therapies.

Abstract

Osteosarcoma (OS) is the most frequent primary bone tumor, mainly affecting children and young adults. Despite therapeutic advances, the 5-year survival rate is 70% but drastically decreases to 20–30% for poor responders to therapies or for patients with metastasis. No real evolution of the survival rates has been observed for four decades, explained by poor knowledge of the origin, difficulties related to diagnosis and the lack of targeted therapies for this pediatric tumor. This review will describe a non-exhaustive overview of osteosarcoma disease from a clinical and biological point of view, describing the origin, diagnosis and therapies.

Inflammation suppresses DLG2 expression decreasing inflammasome formation

Purpose

Loss of expression of DLG2 has been identified in a number of cancers to contribute to the disease by resulting in increased tumor cell proliferation and poor survival. In light of the previous evidence that DLG2 alters the cell cycle and affects proliferation, combined with indications that DLG2 is involved in NLRP3 inflammasome axis we speculated that DLG2 has an immune function. So far, there is no data that clearly elucidates this role, and this study was designed to investigate DLG2 in inflammatory colon disease and in colon cancer as well as its impact on inflammasome induction.

Methods

The DLG2 expression levels were established in publicly available inflammation, colon cancer and mouse model datasets. The overexpression and silencing of DLG2 in colon cancer cells were used to determine the effect of DLG2 expression on the activation of the inflammasome and subsequent cytokine release.

Results

The expression of DLG2 is repressed in inflammatory colon diseases IBD and Ulcerative colitis as well as colorectal cancer tissue compared to healthy individuals. We subsequently show that induction with inflammatory agents in cell and animal models results in a biphasic alteration of DLG2 with an initial increase followed by an ensuing decrease. DLG2 overexpression leads to a significant increase in expression of IL1B, IκBζ and BAX, components that result in inflammasome formation. DLG2 silencing in THP1 cells resulted in increased release of IL-6 into the microenvironment which once used to treat bystander COLO205 cells resulted in an increase in STAT3 phosphorylation and an increase proliferating cells and more cells in the G2/M phase. Restoration of DLG2 to the colon resulted in reduced AKT and S6 signaling.

Conclusion

DLG2 expression is altered in response to inflammation in the gut as well as colon cancer, resulting in altered ability to form inflammasomes.

Trial registration

NCT03072641.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00432-022-04029-7.

Canine and murine models of osteosarcoma

Osteosarcoma (OS) is the most common malignant bone tumor in children. Despite efforts to develop and implement new therapies, patient outcomes have not measurably improved since the 1980s. Metastasis continues to be the main source of patient mortality, with 30% of cases developing metastatic disease within 5 years of diagnosis. Research models are critical in the advancement of cancer research and include a variety of species. For example, xenograft and patient-derived xenograft (PDX) mouse models provide opportunities to study human tumor cells in vivo while transgenic models have offered significant insight into the molecular mechanisms underlying OS development. A growing recognition of naturally occurring cancers in companion species has led to new insights into how veterinary patients can contribute to studies of cancer biology and drug development. The study of canine cases, including the use of diagnostic tissue archives and clinical trials, offers a potential mechanism to further canine and human cancer research. Advancement in the field of OS research requires continued development and appropriate use of animal models. In this review, animal models of OS are described with a focus on the mouse and tumor-bearing pet dog as parallel and complementary models of human OS.

Nuclear translocation of p85β promotes tumorigenesis of PIK3CA helical domain mutant cancer

PI3Ks consist of p110 catalytic subunits and p85 regulatory subunits. PIK3CA, encoding p110α, is frequently mutated in human cancers. Most PIK3CA mutations are clustered in the helical domain or the kinase domain. Here, we report that p85β disassociates from p110α helical domain mutant protein and translocates into the nucleus through a nuclear localization sequence (NLS). Nuclear p85β recruits deubiquitinase USP7 to stabilize EZH1 and EZH2 and enhances H3K27 trimethylation. Knockout of p85β or p85β NLS mutant reduces the growth of tumors harboring a PIK3CA helical domain mutation. Our studies illuminate a novel mechanism by which PIK3CA helical domain mutations exert their oncogenic function. Finally, a combination of Alpelisib, a p110α-specific inhibitor, and an EZH inhibitor, Tazemetostat, induces regression of xenograft tumors harboring a PIK3CA helical domain mutation, but not tumors with either a WT PIK3CA or a PIK3CA kinase domain mutation, suggesting that the drug combination could be an effective therapeutic approach for PIK3CA helical domain mutant tumors.

The mechanisms behind the oncogenic role of the PIK3CA helical domain mutant is poorly understood. Here, the authors show that its oncogenic function depends on the release of p85β from mutated p110α, its translocation to the nucleus and the consequent increased activity of EZH proteins.

Comprehensive Analysis of DNA 5-Methylcytosine and N6-Adenine Methylation by Nanopore Sequencing in Hepatocellular Carcinoma

DNA methylation is a widespread epigenetic signal in human genome. With Nanopore technology, differential methylation modifications including 5-methylcytosine (5mC) and 6-methyladenine (6mA) can be identified. 5mC is the most important modification in mammals, although 6mA may also function in growth and development as well as in pathogenesis. While the role of 5mC at CpG islands in promoter regions associated with transcriptional regulation has been well studied, but the relationship between 6mA and transcription is still unclear. Thus, we collected two pairs of tumor tissues and adjacent normal tissues from hepatocellular carcinoma (HCC) surgical samples for Nanopore sequencing and transcriptome sequencing. It was found that 2,373 genes had both 5mC and 6mA, along with up- and down-regulated methylation sites. These genes were regarded as unstable methylation genes. Compared with 6mA, 5mC had more inclined distribution of unstable methylation sites. Chi-square test showed that the levels of 5mC were consistent with both up- and down-regulated genes, but 6mA was not significant. Moreover, the top three unstable methylation genes, TBC1D3H, CSMD1, and ROBO2, were all related to cancer. Transcriptome and survival analyses revealed four potential tumor suppressor genes including KCNIP4, CACNA1C, PACRG, and ST6GALNAC3. In this study, we firstly proposed to combine 5mC and 6mA methylation sites to explore functional genes, and further research found top of these unstable methylation genes might be functional and some of them could serve as potential tumor suppressor genes. Our study provided a new solution for epigenetic regulation research and therapy of HCC.

Silencing of AKIP1 Suppresses the Proliferation, Migration, and Epithelial-Mesenchymal Transition Process of Glioma Cells by Upregulating DLG2

Gliomas, the most prevalent brain tumors, account for nearly one-third of the all brain and central nervous system (CNS) tumors diagnosed in the USA. The purpose of this study was to discuss the important role of A kinase-interacting protein 1 (AKIP1) in glioma and reveal the potential mechanism. After prediction by CCLE, the expression of AKIP1 was determined by qRT-PCR and western blot. The impacts of AKIP1 knockdown on the proliferation, migration, and invasion were then measured by MTT, colony formation assay, wound healing, and transwell assays. Western blot was used to assess the protein levels of migration and epithelial-mesenchymal transition- (EMT-) related factors. Subsequently, the expression of Disks Large Homolog 2 (DLG2) was predicted by bioinformatics analyses, and the interaction between AKIP1 and DLG2 was confirmed by IP assay, qRT-PCR, and western blot. Finally, DLG2 was further downregulated in glioma cells to detect the association between AKIP1 and DLG2 in the cellular functions of glioma. It was demonstrated that AKIP1 exhibited a high level in glioma cells, and interference of AKIP1 led to reductions in the proliferation, migration, invasion, and EMT of glioma cells. DLG2, which was lowly expressed in glioma cells, demonstrated a negative link to AKIP2. Inhibition of both AKIP2 and DLG2 counteracted the inhibited cellular behaviors on account of AKIP1 interference. To be concluded, this study presented evidence that AKIP1 silencing suppressed the progression of glioma via targeting DLG2, which could provide novel insights to impede the development of glioma.

Immunohistochemical Characterisation of GLUT1, MMP3 and NRF2 in Osteosarcoma

Osteosarcoma (OSA) is an aggressive bone malignancy. Unlike many other malignancies, OSA outcomes have not improved in recent decades. One challenge to the development of better diagnostic and therapeutic methods for OSA has been the lack of well characterized experimental model systems. Spontaneous OSA in dogs provides a good model for the disease seen in people and also remains an important veterinary clinical challenge. We recently used RNA sequencing and qRT-PCR to provide a detailed molecular characterization of OSA relative to non-malignant bone in dogs. We identified differential mRNA expression of the solute carrier family 2 member 1 (SLC2A1/GLUT1), matrix metallopeptidase 3 (MMP3) and nuclear factor erythroid 2–related factor 2 (NFE2L2/NRF2) genes in canine OSA tissue in comparison to paired non-tumor tissue. Our present work characterizes protein expression of GLUT1, MMP3 and NRF2 using immunohistochemistry. As these proteins affect key processes such as Wnt activation, heme biosynthesis, glucose transport, understanding their expression and the enriched pathways and gene ontologies enables us to further understand the potential molecular pathways and mechanisms involved in OSA. This study further supports spontaneous OSA in dogs as a model system to inform the development of new methods to diagnose and treat OSA in both dogs and people.

11q Deletion or ALK Activity Curbs DLG2 Expression to Maintain an Undifferentiated State in Neuroblastoma

High-risk neuroblastomas typically display an undifferentiated or poorly differentiated morphology. It is therefore vital to understand molecular mechanisms that block the differentiation process. We identify an important role for oncogenic ALK-ERK1/2-SP1 signaling in the maintenance of undifferentiated neural crest-derived progenitors through the repression of DLG2, a candidate tumor suppressor gene in neuroblastoma. DLG2 is expressed in the murine "bridge signature" that represents the transcriptional transition state when neural crest cells or Schwann cell precursors differentiate to chromaffin cells of the adrenal gland. We show that the restoration of DLG2 expression spontaneously drives neuroblastoma cell differentiation, highlighting the importance of DLG2 in this process. These findings are supported by genetic analyses of high-risk 11q deletion neuroblastomas, which identified genetic lesions in the DLG2 gene. Our data also suggest that further exploration of other bridge genes may help elucidate the mechanisms underlying the differentiation of NC-derived progenitors and their contribution to neuroblastomas.

Translational oncotargets for immunotherapy: From pet dogs to humans

Preclinical studies in rodent models have been a pivotal role in human clinical research, but many of them fail in the translational process. Spontaneous tumors in pet dogs have the potential to bridge the gap between preclinical models and human clinical trials. Their natural occurrence in an immunocompetent system overcome the limitations of preclinical rodent models. Due to its reasonable cellular, molecular, and genetic homology to humans, the pet dog represents a valuable model to accelerate the translation of preclinical studies to clinical trials in humans, actually with benefits for both species. Moreover, their unique genetic features of breeding and breed-related mutations have contributed to assess and optimize therapeutics in individuals with different genetic backgrounds. This review aims to outline four main immunotherapy approaches - cancer vaccines, adaptive T-cell transfer, antibodies, and cytokines -, under research in veterinary medicine and how they can serve the clinical application crosstalk with humans.

Hydrogen deuterium exchange mass spectrometry identifies the dominant paratope in CD20 antigen binding to the NCD1.2 monoclonal antibody

A comparative canine–human therapeutics model is being developed in B-cell lymphoma through the generation of a hybridoma cell that produces a murine monoclonal antibody specific for canine CD20. The hybridoma cell produces two light chains, light chain-3, and light chain-7. However, the contribution of either light chain to the authentic full-length hybridoma derived IgG is undefined. Mass spectrometry was used to identify only one of the two light chains, light chain-7, as predominating in the full-length IgG. Gene synthesis created a recombinant murine–canine chimeric monoclonal antibody expressing light chain-7 that reconstituted the IgG binding to CD20. Using light chain-7 as a reference sequence, hydrogen deuterium exchange mass spectrometry was used to identify the dominant CDR region implicated in CD20 antigen binding. Early in the deuteration reaction, the CD20 antigen suppressed deuteration at CDR3 (V_H). In later time points, deuterium suppression occurred at CDR2 (V_H) and CDR2 (V_L), with the maintenance of the CDR3 (V_H) interaction. These data suggest that CDR3 (V_H) functions as the dominant antigen docking motif and that antibody aggregation is induced at later time points after antigen binding. These approaches define a methodology for fine mapping of CDR contacts using nested enzymatic reactions and hydrogen deuterium exchange mass spectrometry. These data support the further development of an engineered, synthetic canine–murine monoclonal antibody, focused on CDR3 (V_H), for use as a canine lymphoma therapeutic that mimics the human–murine chimeric anti-CD20 antibody Rituximab.

The Landscape of the Heritable Cancer Genome

Genome-wide association studies (GWAS) have found hundreds of single-nucleotide polymorphisms (SNP) associated with increased risk of cancer. However, the amount of heritable risk explained by SNPs is limited, leaving most of the cancer heritability unexplained. Tumor sequencing projects have shown that causal mutations are enriched in genic regions. We hypothesized that SNPs located in protein coding genes and nearby regulatory regions could explain a significant proportion of the heritable risk of cancer. To perform gene-level heritability analysis, we developed a new method, called Bayesian Gene Heritability Analysis (BAGHERA), to estimate the heritability explained by all genotyped SNPs and by those located in genic regions using GWAS summary statistics. BAGHERA was specifically designed for low heritability traits such as cancer and provides robust heritability estimates under different genetic architectures. BAGHERA-based analysis of 38 cancers reported in the UK Biobank showed that SNPs explain at least 10% of the heritable risk for 14 of them, including late onset malignancies. We then identified 1,146 genes, called cancer heritability genes (CHG), explaining a significant proportion of cancer heritability. CHGs were involved in hallmark processes controlling the transformation from normal to cancerous cells. Importantly, 60 of them also harbored somatic driver mutations, and 27 are tumor suppressors. Our results suggest that germline and somatic mutation information could be exploited to identify subgroups of individuals at higher risk of cancer in the broader population and could prove useful to establish strategies for early detection and cancer surveillance. SIGNIFICANCE: This study describes a new statistical method to identify genes associated with cancer heritability in the broader population, creating a map of the heritable cancer genome with gene-level resolution.See related commentary by Bader, p. 2586.

The loss of DLG2 isoform 7/8, but not isoform 2, is critical in advanced staged neuroblastoma

BACKGROUND

Neuroblastoma is a childhood neural crest tumor showing large clinical and genetic heterogeneity, one form displaying 11q-deletion is very aggressive. It has been shown that 11q-deletion results in decreased expression of DLG2, a gene residing in the deleted region. DLG2 has a number of different isoforms with the main difference is the presence or absence of a L27 domain. The L27 domain containing DLG proteins can form complexes with CASK/MPP and LIN7 protein family members, which will control cell polarity and signaling.

METHODS

We evaluated the DLG gene family and the LIN7 gene family for their expression in differently INSS staged neuroblastoma from publically available data and primary tumors, we included two distinct DLG1 and DLG2 N-terminal transcript isoforms encoding L27 domains for their expression. Functionality of DLG2 isoforms and of LIN7A were evaluated in the 11q-deleted neuroblastoma cell line SKNAS.

RESULTS

In neuroblastoma only two DLG2 isoforms were expressed: isoform 2 and isoform 7/8. Using the array data we could determine that higher expression of DLG members that contain L27 domains correlated to better survival and prognosis. Whilst DLG1 showed a decrease in both isoforms with increased INSS stage, only the full length L27 containing DLG2 transcripts DLG2-isoform 7/8 showed a decrease in expression in high stage neuroblastoma. We could show that the protein encoded by DLG2-isoform 7 could bind to LIN7A, and increased DLG2-isoform 7 gene expression increased the expression of LIN7A, this reduced neuroblastoma cell proliferation and viability, with increased BAX/BCL2 ratio indicating increased apoptosis.

CONCLUSION

We have provided evidence that gene expression of the L27 domain containing DLG2-isoform 7/8 but not L27 domain lacking DLG2-isoform 2 is disrupted in neuroblastoma, in particular in the aggressive subsets of tumors. The presence of the complete L27 domain allows for the binding to LIN7A, which will control cell polarity and signaling, thus affecting cancer cell viability.

Transformed Canine and Murine Mesenchymal Stem Cells as a Model for Sarcoma with Complex Genomics

Simple Summary

Sarcomas are rare cancers of mesenchymal origin, the majority of which are characterized by many copy number alterations, amplifications, or deletions. Because of these complex genomics, it is notoriously difficult to identify driver events of malignant transformation. In this study, we show that murine and canine mesenchymal stem cells (MSCs) can be used to model spontaneous malignant transformation towards sarcomas with complex genomics. We show that these MSCs have an abnormal karyotype, many structural variants, and point mutations at whole genome sequencing analysis, and form sarcomas after injection into mice. Our cross-species analysis reveals that p53 loss is an early event in sarcomagenesis, and it was shown that MSCs with a knock-out in Trp53 transform earlier compared to wild-type MSCs. Our study points to the importance of p53 loss in the transformation process towards sarcomas with complex genomics.

Abstract

Sarcomas are rare mesenchymal tumors with a broad histological spectrum, but they can be divided into two groups based on molecular pathology: sarcomas with simple or complex genomics. Tumors with complex genomics can have aneuploidy and copy number gains and losses, which hampers the detection of early, initiating events in tumorigenesis. Often, no benign precursors are known, which is why good models are essential. The mesenchymal stem cell (MSC) is the presumed cell of origin of sarcoma. In this study, MSCs of murine and canine origin are used as a model to identify driver events for sarcomas with complex genomic alterations as they transform spontaneously after long-term culture. All transformed murine but not canine MSCs formed sarcomas after subcutaneous injection in mice. Using whole genome sequencing, spontaneously transformed murine and canine MSCs displayed a complex karyotype with aneuploidy, point mutations, structural variants, inter-chromosomal translocations, and copy number gains and losses. Cross-species analysis revealed that point mutations in Tp53/Trp53 are common in transformed murine and canine MSCs. Murine MSCs with a cre-recombinase induced deletion of exon 2–10 of Trp53 transformed earlier compared to wild-type murine MSCs, confirming the contribution of loss of p53 to spontaneous transformation. Our comparative approach using transformed murine and canine MSCs points to a crucial role for p53 loss in the formation of sarcomas with complex genomics.

Genomics and transcriptomics in veterinary oncology

The sequencing of the canine genome, combined with additional genomic technologies, has created opportunities for research linking veterinary genomics with naturally occurring cancer in dogs. Also, as numerous canine cancers have features in common with human cancers, comparative studies can be performed to evaluate the use of cancers in dogs as models for human cancer. There have been several reviews of veterinary genomics but, to the best of our knowledge, there has been no comprehensive review of the literature of canine cancer genomics. PubMed and CAB Abstracts databases were searched to retrieve relevant literature using the search terms ‘veterinary’, ‘cancer’ or ‘oncology’, and ‘genomics’ or ‘transcriptomics’. Results were manually assessed and grouped based on the techniques used, the cancer type investigated and genomic lesions targeted. The search resulted in the retrieval of 44 genomic and transcriptomic studies, with the most common technique employed being comparative genomic hybridization. Across both fields, the most commonly studied cancer type was canine osteosarcoma. Genomic and transcriptomic aberrations in canine cancer often reflected those reported in the corresponding human cancers. Analysis of the literature indicated that employing genomic and transcriptomic technologies has been instrumental in developing the understanding of the origin, development and pathogenesis of several canine cancers. However, their use in canine oncology is at an early phase, and there appears to be comparatively little understanding of certain canine cancer types in contrast to their human forms. Aberrations detected in all tumors were tabulated, and the results for osteosarcoma, lymphoma and leukemia, mast cell tumor, transmissible venereal tumor and urothelial carcinoma discussed in detail.

Patient Derived Xenografts for Genome-Driven Therapy of Osteosarcoma

Osteosarcoma (OS) is a rare malignant primary tumor of mesenchymal origin affecting bone. It is characterized by a complex genotype, mainly due to the high frequency of chromothripsis, which leads to multiple somatic copy number alterations and structural rearrangements. Any effort to design genome-driven therapies must therefore consider such high inter- and intra-tumor heterogeneity. Therefore, many laboratories and international networks are developing and sharing OS patient-derived xenografts (OS PDX) to broaden the availability of models that reproduce OS complex clinical heterogeneity. OS PDXs, and new cell lines derived from PDXs, faithfully preserve tumor heterogeneity, genetic, and epigenetic features and are thus valuable tools for predicting drug responses. Here, we review recent achievements concerning OS PDXs, summarizing the methods used to obtain ectopic and orthotopic xenografts and to fully characterize these models. The availability of OS PDXs across the many international PDX platforms and their possible use in PDX clinical trials are also described. We recommend the coupling of next-generation sequencing (NGS) data analysis with functional studies in OS PDXs, as well as the setup of OS PDX clinical trials and co-clinical trials, to enhance the predictive power of experimental evidence and to accelerate the clinical translation of effective genome-guided therapies for this aggressive disease.

Unraveling the chaotic genomic landscape of primary and metastatic canine appendicular osteosarcoma with current sequencing technologies and bioinformatic approaches

Osteosarcoma is a rare disease in children but is one of the most common cancers in adult large breed dogs. The mutational landscape of both the primary and pulmonary metastatic tumor in two dogs with appendicular osteosarcoma (OSA) was comprehensively evaluated using an automated whole genome sequencing, exome and RNA-seq pipeline that was adapted for this study for use in dogs. Chromosomal lesions were the most common type of mutation. The mutational landscape varied substantially between dogs but the lesions within the same patient were similar. Copy number neutral loss of heterozygosity in mutant TP53 was the most significant driver mutation and involved a large region in the middle of chromosome 5. Canine and human OSA is characterized by loss of cell cycle checkpoint integrity and DNA damage response pathways. Mutational profiling of individual patients with canine OSA would be recommended prior to targeted therapy, given the heterogeneity seen in our study and previous studies.

Gene signatures with predictive and prognostic survival values in human osteosarcoma

Osteosarcoma is a common malignancy seen mainly in children and adolescents. The disease is characterized by poor overall prognosis and lower survival due to a lack of predictive markers. Many gene signatures with diagnostic, prognostic, and predictive values were evaluated to achieve better clinical outcomes. Two public data series, GSE21257 and UCSC Xena, were used to identify the minimum number of robust genes needed for a predictive signature to guide prognosis of patients with osteosarcoma. The lasso regression algorithm was used to analyze sequencing data from TCGA-TARGET, and methods such as Cox regression analysis, risk factor scoring, receiving operating curve, KMplot prognosis analysis, and nomogram were used to characterize the prognostic predictive power of the identified genes. Their utility was assessed using the GEO osteosarcoma dataset. Finally, the functional enrichment analysis of the identified genes was performed. A total of twenty-gene signatures were found to have a good prognostic value for predicting patient survival. Gene ontology analysis showed that the key genes related to osteosarcoma were categorized as peptide–antigen binding, clathrin-coated endocytic vesicle membrane, peptide binding, and MHC class II protein complex. The osteosarcoma related genes in these modules were significantly enriched in the processes of antigen processing and presentation, phagocytosis, cell adhesion molecules, Staphylococcus aureus infection. Twenty gene signatures were identified related to osteosarcoma, which would be helpful for predicting prognosis of patients with OS. Further, these signatures can be used to determine the subtypes of osteosarcoma.

Comparative Immunology and Immunotherapy of Canine Osteosarcoma

Approximately 800 people are diagnosed with osteosarcoma (OSA) per year in the USA. Although 70% of patients with localized OSA are cured with multiagent chemotherapy and surgical resection, the prognosis for patients with metastatic or relapsed disease is guarded. The small number of patients diagnosed annually contributes to an incomplete understanding of disease pathogenesis, and challenges in performing appropriately powered clinical trials and detecting correlative biomarkers of response. While mouse models of OSA are becoming increasingly sophisticated, they generally fail to accurately recapitulate tumor heterogeneity, tumor microenvironment (TME), systemic immune dysfunction, and the clinical features of tumor recurrence, metastases, and chemoresistance, which influence outcome. Pet dogs spontaneously develop OSA with an incidence that is 30-50 times higher than humans. Canine OSA parallels the human disease in its clinical presentation, biological behavior, genetic complexity, and therapeutic management. However, despite therapy, most dogs die from metastatic disease within 1 year of diagnosis. Since OSA occurs in immune-competent dogs, immune factors that sculpt tumor immunogenicity and influence responses to immune modulation are in effect. In both species, immune modulation has shown beneficial effects on patient outcome and work is now underway to identify the most effective immunotherapies, combination of immunotherapies, and correlative biomarkers that will further improve clinical response. In this chapter, the immune landscape of canine OSA and the immunotherapeutic strategies used to modulate antitumor immunity in dogs with the disease will be reviewed. From this immunological viewpoint, the value of employing dogs with spontaneous OSA to accelerate and inform the translation of immunotherapies into the human clinic will be underscored.

Improving human cancer therapy through the evaluation of pet dogs

Comparative oncology clinical trials play an important and growing role in cancer research and drug development efforts. These trials, typically conducted in companion (pet) dogs, allow assessment of novel anticancer agents and combination therapies in a veterinary clinical setting that supports serial biologic sample collections and exploration of dose, schedule and corresponding pharmacokinetic/pharmacodynamic relationships. Further, an intact immune system and natural co-evolution of tumour and microenvironment support exploration of novel immunotherapeutic strategies. Substantial improvements in our collective understanding of the molecular landscape of canine cancers have occurred in the past 10 years, facilitating translational research and supporting the inclusion of comparative studies in drug development. The value of the approach is demonstrated in various clinical trial settings, including single-agent or combination response rates, inhibition of metastatic progression and randomized comparison of multiple agents in a head-to-head fashion. Such comparative oncology studies have been purposefully included in the developmental plan for several US FDA-approved and up-and-coming anticancer drugs. Challenges for this field include keeping pace with technology and data dissemination/harmonization, improving annotation of the canine genome and immune system, and generation of canine-specific validated reagents to support integration of correlative biology within clinical trial efforts.

Circ0106714 inhibits tumorigenesis of colorectal cancer by sponging miR-942-5p and releasing DLG2 via Hippo-YAP signaling

This study aimed to investigate the role of circ0106714-miR-942-5p-discs large homolog 2 (DLG2), a novel interactome, in colorectal cancer (CRC). Circ0106714 was found to be the most significantly downregulated circular RNA in CRC using a bioinformatics method, and we researched whether the ability of circ0106714 to sponge miR-942-5p and release DLG2 could affect CRC development via Hippo-YES-associated protein (YAP) signaling. We first employed qRT-PCR and immunoblotting to detect messenger RNA (mRNA) and protein expression, respectively. Live imaging of mice tumor xenografts was then conducted to study the effect of circ0106714 on tumor progression in vivo. Reporter gene assays were subsequently conducted to verify the predicted targeting relationship between circ0106714, miR-942-5p, and DLG2 mRNA in SW480 and HCT116 cell lines. As well as using flow cytometry for both apoptosis and cell cycle profile analyses, CCK-8 and clone foci formation assays were performed to assess cell survival. Wound healing assay and transwell invasion assay were later carried out to evaluate the migration and invasion of the cell lines. Findings revealed that circ0106714 and DLG2 were significantly downregulated, while miR-942-5p was significantly upregulated in human CRC tissues and cell lines. However, circ0106714 upregulation significantly suppressed tumor progression in vivo and inhibited the malignancy phenotypes of tumor cells in vitro by targeting miR-942-5p. Also discovered in this research was that miR-942-5p could directly target DLG2 mRNA, thus enhancing the malignancy phenotypes of CRC cells. We even found that DLG2 overexpression resulted in enhanced phosphorylation of YAP, a critical downstream effector of DLG2. This downstream effector was demonstrated to have a tumor-suppressive capacity in CRC cell lines. In sum, circ0106714 could suppress CRC by sponging miR-942-5p and releasing DLG2, thus promoting YAP phosphorylation.

Somatic structural variation targets neurodevelopmental genes and identifies SHANK2 as a tumor suppressor in neuroblastoma

Neuroblastoma is a malignancy of the developing sympathetic nervous system that accounts for 12% of childhood cancer deaths. Like many childhood cancers, neuroblastoma shows a relative paucity of somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) compared to adult cancers. Here, we assessed the contribution of somatic structural variation (SV) in neuroblastoma using a combination of whole-genome sequencing (WGS) of tumor-normal pairs (n = 135) and single-nucleotide polymorphism (SNP) genotyping of primary tumors (n = 914). Our study design allowed for orthogonal validation and replication across platforms. SV frequency, type, and localization varied significantly among high-risk tumors. MYCN nonamplified high-risk tumors harbored an increased SV burden overall, including a significant excess of tandem duplication events across the genome. Genes disrupted by SV breakpoints were enriched in neuronal lineages and associated with phenotypes such as autism spectrum disorder (ASD). The postsynaptic adapter protein-coding gene, SHANK2, located on Chromosome 11q13, was disrupted by SVs in 14% of MYCN nonamplified high-risk tumors based on WGS and 10% in the SNP array cohort. Expression of SHANK2 was low across human-derived neuroblastoma cell lines and high-risk neuroblastoma tumors. Forced expression of SHANK2 in neuroblastoma cells resulted in significant growth inhibition (P = 2.6 × 10-2 to 3.4 × 10-5) and accelerated neuronal differentiation following treatment with all-trans retinoic acid (P = 3.1 × 10-13 to 2.4 × 10-30). These data further define the complex landscape of somatic structural variation in neuroblastoma and suggest that events leading to deregulation of neurodevelopmental processes, such as inactivation of SHANK2, are key mediators of tumorigenesis in this childhood cancer.

Spontaneous and Induced Animal Models for Cancer Research

Considering the complexity of the current framework in oncology, the relevance of animal models in biomedical research is critical in light of the capacity to produce valuable data with clinical translation. The laboratory mouse is the most common animal model used in cancer research due to its high adaptation to different environments, genetic variability, and physiological similarities with humans. Beginning with spontaneous mutations arising in mice colonies that allow for pursuing studies of specific pathological conditions, this area of in vivo research has significantly evolved, now capable of generating humanized mice models encompassing the human immune system in biological correlation with human tumor xenografts. Moreover, the era of genetic engineering, especially of the hijacking CRISPR/Cas9 technique, offers powerful tools in designing and developing various mouse strains. Within this article, we will cover the principal mouse models used in oncology research, beginning with behavioral science of animals vs. humans, and continuing on with genetically engineered mice, microsurgical-induced cancer models, and avatar mouse models for personalized cancer therapy. Moreover, the area of spontaneous large animal models for cancer research will be briefly presented.

S100A4 mRNA-protein relationship uncovered by measurement noise reduction

Intrinsic biological fluctuation and/or measurement error can obscure the association of gene expression patterns between RNA and protein levels. Appropriate normalization of reverse-transcription quantitative PCR (RT-qPCR) data can reduce technical noise in transcript measurement, thus uncovering such relationships. The accuracy of gene expression measurement is often challenged in the context of cancer due to the genetic instability and "splicing weakness" involved. Here, we sequenced the poly(A) cancer transcriptome of canine osteosarcoma using mRNA-Seq. Expressed sequences were resolved at the level of two consecutive exons to enable the design of exon-border spanning RT-qPCR assays and ranked for stability based on the coefficient of variation (CV). Using the same template type for RT-qPCR validation, i.e. poly(A) RNA, avoided skewing of stability assessment by circular RNAs (circRNAs) and/or rRNA deregulation. The strength of the relationship between mRNA expression of the tumour marker S100A4 and its proportion score of quantitative immunohistochemistry (qIHC) was introduced as an experimental readout to fine-tune the normalization choice. Together with the essential logit transformation of qIHC scores, this approach reduced the noise of measurement as demonstrated by uncovering a highly significant, strong association between mRNA and protein expressions of S100A4 (Spearman's coefficient ρ = 0.72 (p = 0.006)). KEY MESSAGES: • RNA-seq identifies stable pairs of consecutive exons in a heterogeneous tumour. • Poly(A) RNA templates for RT-qPCR avoid bias from circRNA and rRNA deregulation. • HNRNPL is stably expressed across various cancer tissues and osteosarcoma. • Logit transformed qIHC score better associates with mRNA amount. • Quantification of minor S100A4 mRNA species requires poly(A) RNA templates and dPCR.

Low DLG2 gene expression, a link between 11q-deleted and MYCN-amplified neuroblastoma, causes forced cell cycle progression, and predicts poor patient survival

Background

Neuroblastoma (NB) is a childhood neural crest tumor. There are two groups of aggressive NBs, one with MYCN amplification, and another with 11q chromosomal deletion; these chromosomal aberrations are generally mutually exclusive. The DLG2 gene resides in the 11q-deleted region, thus makes it an interesting NB candidate tumor suppressor gene.

Methods

We evaluated the association of DLG2 gene expression in NB with patient outcomes, stage and MYCN status, using online microarray data combining independent NB patient data sets. Functional studies were also conducted using NB cell models and the fruit fly.

Results

Using the array data we concluded that higher DLG2 expression was positively correlated to patient survival. We could also see that expression of DLG2 was inversely correlated with MYCN status and tumor stage. Cell proliferation was lowered in both 11q-normal and 11q-deleted NB cells after DLG2 over expression, and increased in 11q-normal NB cells after DLG2 silencing. Higher level of DLG2 increased the percentage of cells in the G2/M phase and decreased the percentage of cells in the G1 phase. We detected increased protein levels of Cyclin A and Cyclin B in fruit fly models either over expressing dMyc or with RNAi-silenced dmDLG, indicating that both events resulted in enhanced cell cycling. Induced MYCN expression in NB cells lowered DLG2 gene expression, which was confirmed in the fly; when dMyc was over expressed, the dmDLG protein level was lowered, indicating a link between Myc over expression and low dmDLG level.

Conclusion

We conclude that low DLG2 expression level forces cell cycle progression, and that it predicts poor NB patient survival. The low DLG2 expression level could be caused by either MYCN-amplification or 11q-deletion.

Graphical abstract

Predictive modeling for cancer drug discovery using canine models

INTRODUCTION

Rodent models of cancer lack many features associated with the disease in humans. Because dogs closely share an environment with humans, as well as comparable pathophysiology of cancer, they represent a powerful model with which to study novel approaches to cancer treatment.

AREAS COVERED

The authors summarize the weaknesses of rodent models of cancer and the ongoing need for better animal models with which to study potential therapeutic approaches. The homology of cancer in dogs and humans is described, along with examples specific to several common cancer types.

EXPERT OPINION

Laboratory mice and rats will continue to play a central role in cancer research; however, because of a variety of limitations, pet dogs with spontaneous cancer offer unique opportunities for research and should be included in the preclinical development of therapeutic compounds. Environmental homology between dogs and humans, along with biological and molecular similarities present circumstances that strengthen the translational rigor of studies conducted using canine patients. Progress will depend on a sufficient number of dogs to be diagnosed with cancer and available for use in studies; and essential to this will be the availability of enhanced resources for diagnosis of cancer in canine patients and reliable coordination between research scientists, veterinarians, and physicians.

Phenotypic characterization and genome-wide association studies of digital cushion thickness in Holstein cows

The bovine digital cushion is a compression pad between the distal phalanx and sole and has been associated with claw horn disruption lesions. Digital cushion thickness (DCT) is estimated to be moderately heritable. Therefore, the objectives of our study were to examine influences of management and environment on DCT and to identify genetic markers and candidate genes associated with DCT. In a cohort of 502 Holsteins from 5 farms in New York State, DCT and body condition score (BCS) were collected twice, at <137 d prepartum and from 86 to 127 d in milk, corresponding to periods when the digital cushion is thickest and thinnest, respectively, as determined by previous research. Cows underwent sonographic examination of the digital cushion evaluated at the typical sole ulcer site for the right front and hind foot. Linear mixed models were conducted on DCT with the fixed effects of time point, digit, wither height, sacral height, BCS group, and multiple farm system variables separately and included random effects to control for the random subset of cows per farm, repeated measures, and multiple measurements from each cow. The phenotypic results indicated that DCT varied by sample time point, sacral height, parity, digit, BCS group, and wither height. For the genotypic study, 447 DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA). Quality assessment of markers and samples provided a final data set of 431 samples and 579,449 markers. Genome-wide association studies were conducted for DCT testing inheritance models and genetic variation of digit, foot, time point, and average thickness. One marker passed the Bonferroni correction threshold and 26 passed false discovery rate from 4 genome-wide association studies with covariates of sequencing batch plate, parity group, BCS, wither height, and sacral height. Ten candidate genes were identified, with 2 genes on Bos taurus autosomes 24 and 29 involved in biological functions related to the digital cushion: MC4R and DLG2 were related to fat deposition and bone growth, respectively. The genetic markers discovered in this study have the opportunity to be used in breeding programs using genomic selection to select against claw horn disruption lesions and lameness due to associations between the markers and DCT. Further studies on the biologically plausible candidate genes may identify causative genetic variants and how they relate to DCT through gene regulation, expression, structure, or copy number variation.

What's new in bone forming tumours of the skeleton?

Bone tumours are difficult to diagnose and treat, as they are rare and over 60 different subtypes are recognised. The emergence of next-generation sequencing has partly elucidated the molecular mechanisms behind these tumours, including the group of bone forming tumours (osteoma, osteoid osteoma, osteoblastoma and osteosarcoma). Increased knowledge on the molecular mechanism could help to identify novel diagnostic markers and/or treatment options. Osteoid osteoma and osteoblastoma are bone forming tumours without malignant potential that have overlapping morphology. They were recently shown to carry FOS and-to a lesser extent-FOSB rearrangements suggesting that these tumours are closely related. The presence of these rearrangements could help discriminate these entities from other lesions with woven bone deposition. Osteosarcoma is a malignant bone forming tumour for which different histological subtypes are recognised. High-grade osteosarcoma is the prototype of a complex karyotype tumour, and extensive research exploring its molecular background has identified phenomena like chromothripsis and kataegis and some recurrent alterations. Due to lack of specificity, this has not led to a valuable novel diagnostic marker so far. Nevertheless, these studies have also pointed towards potential targetable drivers of which the therapeutic merit remains to be further explored.

Effects of Epigallocatechin Gallate (EGCG) on Urinary Bladder Urothelial Carcinoma-Next-Generation Sequencing and Bioinformatics Approaches

Background and objectives: Bladder urothelial carcinoma is the most common type of genitourinary cancer. Patients with bladder cancer may have limited treatment efficacy related to drug toxicity, resistance or adverse effects, and novel therapeutic strategies to enhance treatment efficacy or increase sensitivity to drugs are of high clinical importance. Epigallocatechin gallate (EGCG) is a polyphenolic compound found in green tea leaves, and a potential anti-cancer agent in various cancer types through modulating and regulating multiple signaling pathways. The current study aimed to explore the role and novel therapeutic targets of EGCG on bladder urothelial carcinoma. Materials and Methods: The BFTC-905 cells, human urinary bladder transitional cell carcinoma (TCC) cell line, were treated with EGCG or water for 24 hours, and the expression profiles of mRNAs and microRNAs were analyzed using next generation sequencing (NGS). The enriched biological functions were determined using different bioinformatics databases. Results: A total of 108 differentially expressed genes in EGCG-treated bladder TCC cells were identified, which were mainly involved in nicotinamide adenine dinucleotide (NAD) biogenesis, inflammatory response and oxidation-reduction metabolism. Moreover, several microRNA-mRNA interactions that potentially participated in the response of bladder TCC to EGCG treatment, including miR-185-3p- ARRB1 (arrestin beta 1), miR-3116- MGAT5B (alpha-1,6-mannosylglycoprotein 6-beta-N-acetylglucosaminyltransferase B), miR-31-5p-TNS1 (tensin 1), miR-642a-5p-TNS1, miR-1226-3p- DLG2 (discs large homolog 2), miR-484-DLG2, and miR-22-3p- PPM1K (protein phosphatase 1K). Conclusions: The current findings provide insights into novel therapeutic targets and underlying mechanisms of action of EGCG treatment in bladder cancer.

Osteosarcoma in the Post Genome Era: Preclinical Models and Approaches to Identify Tractable Therapeutic Targets

PURPOSE OF REVIEW

Osteosarcoma (OS) is the most common cancer of bone, yet is classified as a rare cancer. Treatment and outcomes for OS have not substantively changed in several decades. While the decoding of the OS genome greatly advanced the understanding of the mutational landscape of OS, immediately actionable therapeutic targets were not apparent. Here we describe recent preclinical models that can be leveraged to identify, test, and prioritize therapeutic candidates.

RECENT FINDINGS

The generation of multiple high fidelity murine models of OS, the spontaneous disease that arises in pet dogs, and the establishment of a diverse collection of patient-derived OS xenografts provide a robust preclinical platform for OS. These models enable evidence to be accumulated across multiple stages of preclinical evaluation. Chemical and genetic screening has identified therapeutic targets, often demonstrating cross species activity. Clinical trials in both PDX models and in canine OS have effectively tested new therapies for prioritization. Improving clinical outcomes in OS has proven elusive. The integrated target discovery and testing possible through a cross species platform provides validation of a putative target and may enable the rigorous evaluation of new therapies in models where endpoints can be rapidly assessed.

Cross-Species Gene Expression Analysis Reveals Gene Modules Implicated in Human Osteosarcoma

Background: Osteosarcoma (OS) is one of the malignant bone tumors occurring in both human and canine, and in both of them, it is characterized by a high rate of metastasis and poor prognosis. Cross-species analysis reveals previously neglected molecular or signaling pathways involved in the progression of diseases, and dogs are genetically comparable to humans and live in similar environments. Therefore, the aim of this study was to find out OS hub genes through a cross-species analysis.

Materials and Methods: All the human and canine OS gene expression data obtained by the Affymetrix platform were collected. After quality assessment and normalization, co-expression network was performed using weighted gene co-expression network analysis (WGCNA). Species-specific modules and consensus modules were identified. Protein–protein interaction (PPI) networks analysis was performed based on consensus gene modules. Then, consensus modules were functionally annotated and correlated with clinical traits. Hub nodes were identified by a subnetwork analysis of PPI network and WGCNA module membership. Modules of interest and hub nodes were validated in an external data set.

Results: Three modules for the human network, seven modules for the canine network, and four consensus modules were identified. The consensus module 3 (C3) showed a significant correlation with the metastatic status in the training data set and a significant correlation with metastasis-free survival in the external data set. Cluster of differentiation 86 (CD86) was identified as the hub gene of C3, showing a significant correlation with metastasis-free survival.

Conclusion: Genes in C3 play an important role in OS metastasis, whereas CD86 might be a potential molecular biomarker for OS metastasis.

Canine osteosarcoma genome sequencing identifies recurrent mutations in DMD and the histone methyltransferase gene SETD2

Osteosarcoma (OS) is a rare, metastatic, human adolescent cancer that also occurs in pet dogs. To define the genomic underpinnings of canine OS, we performed multi-platform analysis of OS tumors from 59 dogs, including whole genome sequencing (n = 24) and whole exome sequencing (WES; n = 13) of primary tumors and matched normal tissue, WES (n = 10) of matched primary/metastatic/normal samples and RNA sequencing (n = 54) of primary tumors. We found that canine OS recapitulates features of human OS including low point mutation burden (median 1.98 per Mb) with a trend towards higher burden in metastases, high structural complexity, frequent TP53 (71%), PI3K pathway (37%), and MAPK pathway mutations (17%), and low expression of immune-associated genes. We also identified novel features of canine OS including putatively inactivating somatic SETD2 (42%) and DMD (50%) aberrations. These findings set the stage for understanding OS development in dogs and humans, and establish genomic contexts for future comparative analyses.