Regulatory enhancer elements in solid tumours remain poorly characterized. Here we apply micro-scale chromatin profiling to survey the distal enhancer landscape of primary gastric adenocarcinoma (GC), a leading cause of global cancer mortality. Integrating 110 epigenomic profiles from primary GCs, normal gastric tissues and cell lines, we highlight 36,973 predicted enhancers and 3,759 predicted super-enhancers respectively. Cell-line-defined super-enhancers can be subclassified by their somatic alteration status into somatic gain, loss and unaltered categories, each displaying distinct epigenetic, transcriptional and pathway enrichments. Somatic gain super-enhancers are associated with complex chromatin interaction profiles, expression patterns correlated with patient outcome and dense co-occupancy of the transcription factors CDX2 and HNF4α. Somatic super-enhancers are also enriched in genetic risk SNPs associated with cancer predisposition. Our results reveal a genome-wide reprogramming of the GC enhancer and super-enhancer landscape during tumorigenesis, contributing to dysregulated local and regional cancer gene expression.
Gene expression is regulated by enhancers and super-enhancers, which can be identified by chromatin profiling. Here, the authors surveyed gastric cancer samples and cell lines to identify enhancer elements exhibiting somatic alterations.
Novel, tumor-specific drugs are urgently needed for a breakthrough in cancer therapy. Herein, we generated a first-in-class humanized antibody (PRL3-zumab) against PRL-3, an intracellular tumor-associated phosphatase upregulated in multiple human cancers, for unconventional cancer immunotherapies. We focused on gastric cancer (GC), wherein elevated PRL-3 mRNA levels significantly correlated with shortened overall survival of GC patients. PRL-3 protein was overexpressed in 85% of fresh-frozen clinical gastric tumor samples examined but not in patient-matched normal gastric tissues. Using human GC cell lines, we demonstrated that PRL3-zumab specifically blocked PRL-3+, but not PRL-3–, orthotopic gastric tumors. In this setting, PRL3-zumab had better therapeutic efficacy as a monotherapy, rather than simultaneous combination with 5-fluorouracil or 5-fluorouracil alone. PRL3-zumab could also prevent PRL-3+ tumor recurrence. Mechanistically, we found that intracellular PRL-3 antigens could be externalized to become “extracellular oncotargets” that serve as bait for PRL3-zumab binding to potentially bridge and recruit immunocytes into tumor microenvironments for killing effects on cancer cells. In summary, our results document a comprehensive cancer therapeutic approach to specific antibody-targeted therapy against the PRL-3 oncotarget as a case study for developing antibodies against other intracellular targets in drug discovery.
A humanized antibody against the phosphatase PRL-3 specifically inhibits PRL-3+ orthotopic gastric tumors in mice and protects against postsurgical tumor recurrence.
Metastasis is a significant health issue. The standard mode of care is combination of chemotherapy and targeted therapeutics but the 5-year survival rate remains low. New/better drug targets that can improve outcomes of patients with metastatic disease are needed. Metastasis is a complex process, with each step conferred by a set of genetic aberrations. Mapping the molecular changes associated with metastasis improves our understanding of the etiology of this disease and contributes to the pipeline of targeted therapeutics. Here, phosphoproteomics of a xenograft-derived in vitro model comprising 4 isogenic cell lines with increasing metastatic potential implicated Transient Receptor Potential Vanilloid subtype 4 in breast cancer metastasis. TRPV4 mRNA levels in breast, gastric and ovarian cancers correlated with poor clinical outcomes, suggesting a wide role of TRPV4 in human epithelial cancers. TRPV4 was shown to be required for breast cancer cell invasion and transendothelial migration but not growth/proliferation. Knockdown of Trpv4 significantly reduced the number of metastatic nodules in mouse xenografts leaving the size unaffected. Overexpression of TRPV4 promoted breast cancer cell softness, blebbing, and actin reorganization. The findings provide new insights into the role of TRPV4 in cancer extravasation putatively by reducing cell rigidity through controlling the cytoskeleton at the cell cortex.
Advances in high-throughput genomic technologies coupled with a growing number of genomic results potentially useful in clinical care have led to ground-breaking genomic medicine implementation programs in various nations. Many of these innovative programs capitalize on unique local capabilities arising from the structure of their health care systems or their cultural or political milieu, as well as from unusual burdens of disease or risk alleles. Many such programs are being conducted in relative isolation and might benefit from sharing of approaches and lessons learned in other nations. The National Human Genome Research Institute recently brought together 25 of these groups from around the world to describe and compare projects, examine the current state of implementation and desired near-term capabilities, and identify opportunities for collaboration to promote the responsible implementation of genomic medicine.
The wide variety of nascent programs in diverse settings demonstrates that implementation of genomic medicine is expanding globally in varied and highly innovative ways. Opportunities for collaboration abound in the areas of evidence generation, health information technology, education, workforce development, pharmacogenomics, and policy and regulatory issues. Several international organizations that are already facilitating effective research collaborations should engage to ensure implementation proceeds collaboratively without potentially wasteful duplication. Efforts to coalesce these groups around concrete but compelling signature projects, such as global eradication of genetically-mediated drug reactions or developing a truly global genomic variant data resource across a wide number of ethnicities, would accelerate appropriate implementation of genomics to improve clinical care world-wide.
medical genomics; implementation; global collaborations; practice standards; pharmacogenomics; personalized medicine; precision medicine
Cholangiocarcinoma (CCA) is the primary cancer of the bile duct system. The role of bile duct tissue microbiomes in CCA tumorigenesis is unestablished. To address this, sixty primary CCA tumors and matched normals, from both liver fluke (Opisthorchis viverrini) associated (OVa, n = 28) and non-O. viverrini associated (non-OVa, n = 32) cancers, were profiled using high-throughput 16S rRNA sequencing. A distinct, tissue-specific microbiome dominated by the bacterial families Dietziaceae, Pseudomonadaceae and Oxalobacteraceae was observed in bile duct tissues. Systemic perturbation of the microbiome was noted in tumor and paired normal samples (vs non-cancer normals) for several bacterial families with a significant increase in Stenotrophomonas species distinguishing tumors vs paired normals. Comparison of parasite associated (OVa) vs non-associated (non-OVa) groups identified enrichment for specific enteric bacteria (Bifidobacteriaceae, Enterobacteriaceae and Enterococcaceae). One of the enriched families, Bifidobacteriaceae, was found to be dominant in the O. viverrini microbiome, providing a mechanistic link to the parasite. Functional analysis and comparison of CCA microbiomes revealed higher potential for producing bile acids and ammonia in OVa tissues, linking the altered microbiota to carcinogenesis. These results define how the unique microbial communities resident in the bile duct, parasitic infections and the tissue microenvironment can influence each other, and contribute to cancer.
•Stenotrophomonas, implicated in bile duct infections, is enriched in tumor tissues of non-fluke related cholangiocarcinoma.•O. viverrini infection may alter composition of bile duct tissue microbiomes.•Enteric bacteria with metabolic outputs linked to carcinogenesis are enriched in O. viverrini associated tissue microbiomes.
The link between microbiota and cancer of the gastrointestinal (GI) tract has been extensively studied. However, the role of tissue microbiome in cholangiocarcinoma (CCA), cancer of the bile duct (an organ connected to the GI tract), is largely unknown. In this study, we detected intriguing compositional differences in the tissue microbiomes of liver fluke related and non-related CCA. Taken together, our data suggests a connection between parasitic infections, tissue microbiome alterations, tissue micro-environment changes and CCA development.
Microbiome; Cancer; Cholangiocarcinoma; Liver fluke
In gastric cancer (GC), the main subtypes (diffuse and intestinal types) differ in pathological characteristics, with diffuse GC exhibiting early disseminative and invasive behaviour. A distinctive feature of diffuse GC is loss of intercellular adhesion. Although widely attributed to mutations in the CDH1 gene encoding E-cadherin, a significant percentage of diffuse GC do not harbor CDH1 mutations. We found that the expression of the actin-modulating cytoskeletal protein, gelsolin, is significantly higher in diffuse-type compared to intestinal-type GCs, using immunohistochemical and microarray analysis. Furthermore, in GCs with wild-type CDH1, gelsolin expression correlated inversely with CDH1 gene expression. Downregulating gelsolin using siRNA in GC cells enhanced intercellular adhesion and E-cadherin expression, and reduced invasive capacity. Interestingly, hepatocyte growth factor (HGF) induced increased gelsolin expression, and gelsolin was essential for HGF-medicated cell scattering and E-cadherin transcriptional repression through Snail, Twist and Zeb2. The HGF-dependent effect on E-cadherin was found to be mediated by interactions between gelsolin and PI3K-Akt signaling. This study reveals for the first time a function of gelsolin in the HGF/cMet oncogenic pathway, which leads to E-cadherin repression and cell scattering in gastric cancer. Our study highlights gelsolin as an important pro-disseminative factor contributing to the aggressive phenotype of diffuse GC.
gelsolin; gastric cancer; E-Cadherin; hepatocyte growth factor (HGF); cancer invasion
Ion channels are well characterised in model organisms, principally because of the availability of functional genomic tools and datasets for these species. This contrasts the situation, for example, for parasites of humans and animals, whose genomic and biological uniqueness means that many genes and their products cannot be annotated. As ion channels are recognised as important drug targets in mammals, the accurate identification and classification of parasite channels could provide major prospects for defining unique targets for designing novel and specific anti-parasite therapies. Here, we established a reliable bioinformatic pipeline for the identification and classification of ion channels encoded in the genome of the cancer-causing liver fluke Opisthorchis viverrini, and extended its application to related flatworms affecting humans.
We built an ion channel identification + classification pipeline (called MuSICC), employing an optimised support vector machine (SVM) model and using the Kyoto Encyclopaedia of Genes and Genomes (KEGG) classification system. Ion channel proteins were first identified and grouped according to amino acid sequence similarity to classified ion channels and the presence and number of ion channel-like conserved and transmembrane domains. Predicted ion channels were then classified to sub-family using a SVM model, trained using ion channel features.
Following an evaluation of this pipeline (MuSICC), which demonstrated a classification sensitivity of 95.2 % and accuracy of 70.5 % for known ion channels, we applied it to effectively identify and classify ion channels in selected parasitic flatworms.
MuSICC provides a practical and effective tool for the identification and classification of ion channels of parasitic flatworms, and should be applicable to a broad range of organisms that are evolutionarily distant from taxa whose ion channels are functionally characterised.
Electronic supplementary material
The online version of this article (doi:10.1186/s13071-016-1428-2) contains supplementary material, which is available to authorized users.
Ion channels; Identification; Classification; Parasitic flatworms; Bioinformatic pipeline
Peritoneal dissemination is the most frequent, incurable metastasis occurring in patients with advanced gastric cancer (GC). However, molecular mechanisms driving peritoneal dissemination still remain poorly understood. Here, we aimed to provide novel insights into the molecular mechanisms that drive the peritoneal dissemination of GC. We performed combined expression analysis with in vivo-selected metastatic cell lines and samples from 200 GC patients to identify driver genes of peritoneal dissemination. The driver-gene functions associated with GC dissemination were examined using a mouse xenograft model. We identified a peritoneal dissemination-associated expression signature, whose profile correlated with those of genes related to development, focal adhesion, and the extracellular matrix. Among the genes comprising the expression signature, we identified that discoidin-domain receptor 2 (DDR2) as a potential regulator of peritoneal dissemination. The DDR2 was upregulated by the loss of DNA methylation and that DDR2 knockdown reduced peritoneal metastasis in a xenograft model. Dasatinib, an inhibitor of the DDR2 signaling pathway, effectively suppressed peritoneal dissemination. DDR2 was identified as a driver gene for GC dissemination from the combined expression signature and can potentially serve as a novel therapeutic target for inhibiting GC peritoneal dissemination.
Inherited cardiac conditions (ICCs) are characterised by marked genetic and allelic heterogeneity and require extensive sequencing for genetic characterisation. We iteratively optimised a targeted gene capture panel for ICCs that includes disease-causing, putatively pathogenic, research and phenocopy genes (n = 174 genes). We achieved high coverage of the target region on both MiSeq (>99.8 % at ≥20× read depth, n = 12) and NextSeq (>99.9 % at ≥20×, n = 48) platforms with 100 % sensitivity and precision for single nucleotide variants and indels across the protein-coding target on the MiSeq. In the final assay, 40 out of 43 established ICC genes informative in clinical practice achieved complete coverage (100 % at ≥20×). By comparison, whole exome sequencing (WES; ∼80×), deep WES (∼500×) and whole genome sequencing (WGS; ∼70×) had poorer performance (88.1, 99.2 and 99.3 % respectively at ≥20×) across the ICC target. The assay described here delivers highly accurate and affordable sequencing of ICC genes, complemented by accessible cloud-based computation and informatics. See Editorial in this issue (DOI: 10.1007/s12265-015-9667-8).
Electronic supplementary material
The online version of this article (doi:10.1007/s12265-016-9673-5) contains supplementary material, which is available to authorized users.
Inherited cardiac conditions; Targeted sequencing; Whole exome sequencing; Whole genome sequencing; Genetics; Diagnostics
Gene expression profiling has contributed greatly to cancer research. However, expression-driven biomarker discovery in metastatic gastric cancer (mGC) remains unclear. A gene expression profile predicting RAD001 response in refractory GC was explored in this study.
Total RNA isolated from 54 tumour specimens from patients with mGC, prior to RAD001 treatment, was analysed via the NanoString nCounter gene expression assay. This assay targeted 477 genes representing 10 different GC-related oncogenic signalling and molecular subtype-specific expression signatures. Gene expression profiles were correlated with patient clinicopathological variables.
NanoString data confirmed similar gene expression profiles previously identified by microarray analysis. Signature I with 3 GC subtypes (mesenchymal, metabolic and proliferative) showed approximately 90% concordance where the mesenchymal and proliferative subtypes were significantly associated with signet ring cell carcinoma and the WHO classified tubular adenocarcinoma GC, respectively (p=0.042). Single-gene-level correlations with patient clinicopathological variables showed strong associations between FHL1 expression (mesenchymal subtype) and signet ring cell carcinoma, and NEK2, OIP5, PRC1, TPX2 expression (proliferative subtype) with tubular adenocarcinoma (adjusted p<0.05). Increased BRCA2 (p=0.040) and MMP9 (p=0.045) expression was significantly associated with RAD001 good response and longer progression-free survival outcome (BRCA2, p=0.012, HR 0.370 95% CI (0.171 to 0.800); MMP9, p=0.010, HR 0.359 95% CI (0.166 to 0.779)). In contrast, increased BTC (p=0.035) expression was significantly associated with RAD001 poor response and poor progression-free survival (p=0.031, HR 2.336 95% CI (1.079 to 5.059) by univariate Cox regression analysis.
Microarray results are highly reproducible with NanoString nCounter gene expression profiling. Additionally, BRCA2 and MMP9 expression are potential predictive biomarkers for good response in RAD001-treated mGC.
gastric cancer; NanoString; RAD001; metastasis; gene expression
To date, most molecular investigations of schistosomatids have focused principally on blood flukes (schistosomes) of humans. Despite the clinical importance of cercarial dermatitis in humans caused by Trichobilharzia regenti and the serious neuropathologic disease that this parasite causes in its permissive avian hosts and accidental mammalian hosts, almost nothing is known about the molecular aspects of how this fluke invades its hosts, migrates in host tissues and how it interacts with its hosts’ immune system. Here, we explored selected aspects using a transcriptomic-bioinformatic approach. To do this, we sequenced, assembled and annotated the transcriptome representing two consecutive life stages (cercariae and schistosomula) of T. regenti involved in the first phases of infection of the avian host. We identified key biological and metabolic pathways specific to each of these two developmental stages and also undertook comparative analyses using data available for taxonomically related blood flukes of the genus Schistosoma. Detailed comparative analyses revealed the unique involvement of carbohydrate metabolism, translation and amino acid metabolism, and calcium in T. regenti cercariae during their invasion and in growth and development, as well as the roles of cell adhesion molecules, microaerobic metabolism (citrate cycle and oxidative phosphorylation), peptidases (cathepsins) and other histolytic and lysozomal proteins in schistosomula during their particular migration in neural tissues of the avian host. In conclusion, the present transcriptomic exploration provides new and significant insights into the molecular biology of T. regenti, which should underpin future genomic and proteomic investigations of T. regenti and, importantly, provides a useful starting point for a range of comparative studies of schistosomatids and other trematodes.
Despite the clinical importance of Trichobilharzia regenti in bird hosts and as a cause of cercarial dermatitis in humans, almost nothing is known about the molecular aspects of this fluke and its interactions with its hosts. Here, we sequenced, assembled and annotated the transcriptome representing two life stages (cercariae and schistosomula) of T. regenti involved in the first phases of infection of the bird host. We identified key biological and metabolic pathways specific to each of these two developmental stages and also undertook comparative analyses using data available for related flukes. Detailed analyses showed the unique involvement of carbohydrate metabolism, translation and amino acid metabolism, and calcium in T. regenti cercariae during invasion and in growth and development, as well as cell adhesion molecules, microaerobic metabolism (citrate cycle and oxidative phosphorylation), peptidases (cathepsins) and other histolytic and lysozomal proteins in schistosomula during migration in neural tissues. These molecular insights into T. regenti biology should support future genomic and proteomic investigations of T. regenti, and comparative studies of flatworms.
Despite the approval of numerous molecular targeted drugs, long-term antiproliferative efficacy is rarely achieved and therapy resistance remains a central obstacle of cancer care. Combined inhibition of multiple cancer-driving pathways promises to improve antiproliferative efficacy. HIF-1 is a driver of gastric cancer and considered to be an attractive target for therapy. We noted that gastric cancer cells are able to functionally compensate the stable loss of HIF-1α. Via transcriptomics we identified a group of upregulated genes in HIF-1α-deficient cells and hypothesized that these genes confer survival upon HIF-1α loss. Strikingly, simultaneous knock-down of HIF-1α and Annexin A1 (ANXA1), one of the identified genes, resulted in complete cessation of proliferation. Using stable isotope-resolved metabolomics, oxidative and reductive glutamine metabolism was found to be significantly impaired in HIF-1α/ANXA1-deficient cells, potentially explaining the proliferation defect. In summary, we present a conceptually novel application of stable gene inactivation enabling in-depth deconstruction of resistance mechanisms. In theory, this experimental approach is applicable to any cancer-driving gene or pathway and promises to identify various new targets for combination therapies.
cancer therapy; Annexin A1; cancer metabolism; HIF-1; induced essentiality
Trichinellosis is a globally important food-borne parasitic disease of humans caused by roundworms of the Trichinella complex. Extensive biological diversity is reflected in substantial ecological and genetic variability within and among Trichinella taxa, and major controversy surrounds the systematics of this complex. Here we report the sequencing and assembly of 16 draft genomes representing all 12 recognized Trichinella species and genotypes, define protein-coding gene sets and assess genetic differences among these taxa. Using thousands of shared single-copy orthologous gene sequences, we fully reconstruct, for the first time, a phylogeny and biogeography for the Trichinella complex, and show that encapsulated and non-encapsulated Trichinella taxa diverged from their most recent common ancestor ∼21 million years ago (mya), with taxon diversifications commencing ∼10−7 mya.
Trichinellosis is a globally important food-borne disease caused by roundworms of the Trichinella complex. Here the authors present genomic sequences representing all 12 recognized Trichinella species and genotypes, and reconstruct their phylogeny and biogeography.
Dacomitinib, an irreversible pan-HER inhibitor, had shown modest clinical activity in squamous cell carcinoma of head and neck (SCCHN) patients. Therefore, validated predictive biomarkers are required to identify patients most likely to benefit from this therapeutic option. To characterize the genetic landscape of cisplatin-treated SCCHN genomes and identify potential predictive biomarkers for dacomitinib sensitivity, we performed whole exome sequencing on 18 cisplatin-resistant metastatic SCCHN tumors and their matched germline DNA. Platinum-based chemotherapy elevated the mutation rates of SCCHN compared to chemotherapy-naïve SCCHNs. Cisplatin-treated SCCHN genomes uniquely exhibited a novel mutational signature characterized by C:G to A:T transversions at CCR sequence contexts that may have arisen due to error-prone translesional synthesis. Somatic mutations in REV3L, the gene encoding the catalytic subunit of DNA polymerase ζ involved in translesional synthesis, are significantly enriched in a subset of patients who derived extended clinical benefit to dacomitinib (P = 0.04). Functional assays showed that loss-of-function of REV3L dramatically enhanced the sensitivity of SCCHN cells to dacomitinib by the loss of both translesion synthesis and homologous recombination pathways. Our data suggest that the ‘platinum’ mutational signature and inactivation of REV3L may inform treatment options in patients of recurrent SCCHN.
Differences in gastric cancer (GC) clinical outcomes between patients in Asian and non-Asian countries has been historically attributed to variability in clinical management. However, recent international Phase III trials suggest that even with standardised treatments, GC outcomes differ by geography. Here, we investigated gene expression differences between Asian and non-Asian GCs, and if these molecular differences might influence clinical outcome.
We compared gene expression profiles of 1016 GCs from six Asian and three non-Asian GC cohorts, using a two-stage meta-analysis design and a novel biostatistical method (RUV-4) to adjust for technical variation between cohorts. We further validated our findings by computerised immunohistochemical analysis on two independent tissue microarray (TMA) cohorts from Asian and non-Asian localities (n=665).
Gene signatures differentially expressed between Asians and non-Asian GCs were related to immune function and inflammation. Non-Asian GCs were significantly enriched in signatures related to T-cell biology, including CTLA-4 signalling. Similarly, in the TMA cohorts, non-Asian GCs showed significantly higher expression of T-cell markers (CD3, CD45R0, CD8) and lower expression of the immunosuppressive T-regulatory cell marker FOXP3 compared to Asian GCs (p<0.05). Inflammatory cell markers CD66b and CD68 also exhibited significant cohort differences (p<0.05). Exploratory analyses revealed a significant relationship between tumour immunity factors, geographic locality-specific prognosis, and postchemotherapy outcomes.
Analyses of >1600 GCs suggest that Asian and non-Asian GCs exhibit distinct tumour immunity signatures related to T-cell function. These differences may influence geographical differences in clinical outcome, and the design of future trials particularly in immuno-oncology.
GASTRIC CANCER; GENE EXPRESSION; MOLECULAR PATHOLOGY; IMMUNOLOGY; MOLECULAR MECHANISMS
The blood fluke Schistosoma haematobium causes urogenital schistosomiasis, a neglected tropical disease (NTD) that affects more than 110 million people. Treating this disease by targeted or mass administration with a single chemical, praziquantel, carries the risk that drug resistance will develop in this pathogen. Therefore, there is an imperative to search for new drug targets in S. haematobium and other schistosomes. In this regard, protein kinases have potential, given their essential roles in biological processes and as targets for drugs already approved by the US Food and Drug Administration (FDA) for use in humans. In this context, we defined here the kinome of S. haematobium using a refined bioinformatic pipeline. We classified, curated and annotated predicted kinases, and assessed the developmental transcription profiles of kinase genes. Then, we prioritised a panel of kinases as potential drug targets and inferred chemicals that bind to them using an integrated bioinformatic pipeline. Most kinases of S. haematobium are very similar to those of its congener, S. mansoni, offering the prospect of designing chemicals that kill both species. Overall, this study provides a global insight into the kinome of S. haematobium and should assist the repurposing or discovery of drugs against schistosomiasis.
Schistosomiasis is a neglected tropical disease that affects more than 200 million people worldwide. The main disease-causing agents, Schistosoma japonicum, S. mansoni and S. haematobium, are blood flukes that have complex life cycles involving a snail intermediate host. In Asia, S. japonicum causes hepatointestinal disease (schistosomiasis japonica) and is challenging to control due to a broad distribution of its snail hosts and range of animal reservoir hosts. In China, extensive efforts have been underway to control this parasite, but genetic variability in S. japonicum populations could represent an obstacle to eliminating schistosomiasis japonica. Although a draft genome sequence is available for S. japonicum, there has been no previous study of molecular variation in this parasite on a genome-wide scale. In this study, we conducted the first deep genomic exploration of seven S. japonicum populations from mainland China, constructed phylogenies using mitochondrial and nuclear genomic data sets, and established considerable variation between some of the populations in genes inferred to be linked to key cellular processes and/or pathogen-host interactions. Based on the findings from this study, we propose that verifying intraspecific conservation in vaccine or drug target candidates is an important first step toward developing effective vaccines and chemotherapies against schistosomiasis.
Despite advances in therapeutics, outcomes for hepatocellular carcinoma (HCC) remain poor and there is an urgent need for efficacious systemic therapy. Unfortunately, drugs that are successful in preclinical studies often fail in the clinical setting, and we hypothesize that this is due to functional differences between primary tumors and commonly used preclinical models. In this study, we attempt to answer this question by comparing tumor morphology and gene expression profiles between primary tumors, xenografts and HCC cell lines.
Hep G2 cell lines and tumor cells from patient tumor explants were subcutaneously (ectopically) injected into the flank and orthotopically into liver parenchyma of Mus Musculus SCID mice. The mice were euthanized after two weeks. RNA was extracted from the tumors, and gene expression profiling was performed using the Gene Chip Human Genome U133 Plus 2.0. Principal component analyses (PCA) and construction of dendrograms were conducted using Partek genomics suite.
PCA showed that the commonly used HepG2 cell line model and its xenograft counterparts were vastly different from all fresh primary tumors. Expression profiles of primary tumors were also significantly divergent from their counterpart patient-derived xenograft (PDX) models, regardless of the site of implantation. Xenografts from the same primary tumors were more likely to cluster together regardless of site of implantation, although heat maps showed distinct differences in gene expression profiles between orthotopic and ectopic models.
The data presented here challenges the utility of routinely used preclinical models. Models using HepG2 were vastly different from primary tumors and PDXs, suggesting that this is not clinically representative. Surprisingly, site of implantation (orthotopic versus ectopic) resulted in limited impact on gene expression profiles, and in both scenarios xenografts differed significantly from the original primary tumors, challenging the long-held notion that orthotopic PDX model is the gold standard preclinical model for HCC.
Electronic supplementary material
The online version of this article (doi:10.1186/s12885-015-1814-8) contains supplementary material, which is available to authorized users.
Hepatocellular carcinoma; Ectopic; Orthotopic; Xenograft; HepG2 cell line
Resistance to chemotherapy is common in gastroesophageal cancer. Mechanisms of resistance are incompletely characterised and there are no predictive biomarkers in clinical practice for cytotoxic drugs. We used new cell line models to characterise novel chemotherapy resistance mechanisms and validated them in tumour specimens to identify new targets and biomarkers for gastroesophageal cancer.
Cell lines were selected for resistance to oxaliplatin, cisplatin and docetaxel and gene expression examined using Affymetrix Exon 1.0 ST arrays. Leads were validated by qRT-PCR and HPLC of tumour metabolites. Protein expression and pharmacological inhibition of lead target SPHK1 was evaluated in independent cell lines, and by immunohistochemistry in gastroesophageal cancer patients.
Genes with differential expression in drug resistant cell lines compared to the parental cell line they were derived from, were identified for each drug resistant cell line. Biological pathway analysis of these gene lists, identified over-represented pathways, and only 3 pathways - lysosome, sphingolipid metabolism and p53 signalling- were identified as over-represented in these lists for all three cytotoxic drugs investigated. The majority of genes differentially expressed in chemoresistant cell lines from these pathways, were involved in metabolism of glycosphingolipids and sphingolipids in lysosomal compartments suggesting that sphingolipids might be important mediators of cytotoxic drug resistance in gastroeosphageal cancers . On further investigation, we found that drug resistance (IC50) was correlated with increased sphingosine kinase 1(SPHK1) mRNA and also with decreased sphingosine-1-phosphate lysase 1(SGPL1) mRNA. SPHK1 and SGPL1 gene expression were inversely correlated. SPHK1:SGPL1 ratio correlated with increased cellular sphingosine-1-phosphate (S1P), and S1P correlated with drug resistance (IC50). High SPHK1 protein correlated with resistance to cisplatin (IC50) in an independent gastric cancer cell line panel and with survival of patients treated with chemotherapy prior to surgery but not in patients treated with surgery alone. Safingol a SPHK1 inhibitor, was cytotoxic as a single agent and acted synergistically with cisplatin in gastric cancer cell lines.
Agents that inhibit SPHK1 or S1P could overcome cytotoxic drug resistance in gastroesophageal cancer. There are several agents in early phase human trials including Safingol that could be combined with chemotherapy or used in patients progressing after chemotherapy.
Electronic supplementary material
The online version of this article (doi:10.1186/s12885-015-1718-7) contains supplementary material, which is available to authorized users.
Gastroesophageal cancer; Chemoresistance; Sphingosine-1-phosphate; Sphingosine kinase 1; Sphingosine-1-phopshate lyase
Carcinoma of the oral tongue (OTSCC) is the most common malignancy of the oral cavity, characterized by frequent recurrence and poor survival. The last three decades has witnessed a change in the OTSCC epidemiological profile, with increasing incidence in younger patients, females and never-smokers. Here, we sought to characterize the OTSCC genomic landscape and to determine factors that may delineate the genetic basis of this disease, inform prognosis and identify targets for therapeutic intervention.
Seventy-eight cases were subjected to whole-exome (n = 18) and targeted deep sequencing (n = 60).
While the most common mutation was in TP53, the OTSCC genetic landscape differed from previously described cohorts of patients with head and neck tumors: OTSCCs demonstrated frequent mutations in DST and RNF213, while alterations in CDKN2A and NOTCH1 were significantly less frequent. Despite a lack of previously reported NOTCH1 mutations, integrated analysis showed enrichments of alterations affecting Notch signaling in OTSCC. Importantly, these Notch pathway alterations were prognostic on multivariate analyses. A high proportion of OTSCCs also presented with alterations in drug targetable and chromatin remodeling genes. Patients harboring mutations in actionable pathways were more likely to succumb from recurrent disease compared with those who did not, suggesting that the former should be considered for treatment with targeted compounds in future trials.
Our study defines the Asian OTSCC mutational landscape, highlighting the key role of Notch signaling in oral tongue tumorigenesis. We also observed somatic mutations in multiple therapeutically relevant genes, which may represent candidate drug targets in this highly lethal tumor type.
Electronic supplementary material
The online version of this article (doi:10.1186/s13073-015-0219-2) contains supplementary material, which is available to authorized users.
Cancer is a genetic condition driven by a series of both inherited (germline) and tumour specific (somatic) mutations. The analysis of tumour to identify somatic alterations and generate a genetic profile has given rise to a personalised approach for treating oncology patients, guiding optimal therapy. Next generation sequencing (NGS) technology is now routinely utilised in oncology to screen for both genes known to be associated with improving treatment response, as well as, genes of uncertain significance to contribute to new findings in cancer treatment. In the development of genome wide mutation analysis of tumour tissue, germline DNA has also been routinely collected for comparative analysis in order to identity tumour specific mutations. As the identification of germline mutations is secondary to the main purpose of tumour testing for targeted treatment, the management of germline findings has been extensively debated, in particular, around if and when to communicate these findings to patients. The role of the genetic counsellor in this process has been pertinent to these discussions. Clinical practice, however, is shifting towards a preference of sequencing tumour tissue alone to characterise its molecular profile for reasons such as reducing cost and facilitating logistics of sample collection. This method raises the question of how the testing of tumour alone contributes to the identification of germline mutations. The laboratory at POLARIS, Singapore, has generated genomic data of 130 colon tumours using a panel of 84 genes, including 29 germline susceptibility genes associated with known inherited syndromes. The analysis of this data has given an insight into the proportion of mutations that are potentially germline, which in turn will impact how patients are consented for tumour profiling and the delivery of the results. In generating these findings, however, the analysis process has demonstrated a number of the complexities in interpreting the pathogenicity of the variants, as well as, determining the somatic or germline origin of the variant. The challenges in interpreting germline mutations from tumour next generation data, the importance of conveying a germline finding to the patient and the integration of the genetic counsellor in this process will be discussed.
Genetic counsellor; precision oncology; next generation sequencing (NGS); tumour profiling
Microsatellite instability (MSI) is a form of hypermutation that occurs in some tumors due to defects in cellular DNA mismatch repair. MSI is characterized by frequent somatic mutations (i.e., cancer-specific mutations) that change the length of simple repeats (e.g., AAAAA…., GATAGATAGATA...). Clinical MSI tests evaluate the lengths of a handful of simple repeat sites, while next-generation sequencing can assay many more sites and offers a much more complete view of their somatic mutation frequencies. Using somatic mutation data from the exomes of a 361-tumor training set, we developed classifiers to determine MSI status based on four machine-learning frameworks. All frameworks had high accuracy, and after choosing one we determined that it had >98% concordance with clinical tests in a separate 163-tumor test set. Furthermore, this classifier retained high concordance even when classifying tumors based on subsets of whole-exome data. We have released a CRAN R package, MSIseq, based on this classifier. MSIseq is faster and simpler to use than software that requires large files of aligned sequenced reads. MSIseq will be useful for genomic studies in which clinical MSI test results are unavailable and for detecting possible misclassifications by clinical tests.
Orexinergic neurotransmission is involved in mediating temperature responses to methamphetamine (Meth). In experiments in rats, SB-334867 (SB), an antagonist of orexin receptors (OX1R), at a dose of 10 mg/kg decreases late temperature responses (t>60 min) to an intermediate dose of Meth (5 mg/kg). A higher dose of SB (30 mg/kg) attenuates temperature responses to low dose (1 mg/kg) of Meth and to stress. In contrast, it significantly exaggerates early responses (t<60 min) to intermediate and high doses (5 and 10 mg/kg) of Meth. As pretreatment with SB also inhibits temperature response to the stress of injection, traditional statistical analysis of temperature responses is difficult.
We have developed a mathematical model that explains the complexity of temperature responses to Meth as the interplay between excitatory and inhibitory nodes. We have extended the developed model to include the stress of manipulations and the effects of SB. Stress is synergistic with Meth on the action on excitatory node. Orexin receptors mediate an activation of on both excitatory and inhibitory nodes by low doses of Meth, but not on the node activated by high doses (HD). Exaggeration of early responses to high doses of Meth involves disinhibition: low dose of SB decreases tonic inhibition of HD and lowers the activation threshold, while the higher dose suppresses the inhibitory component. Using a modeling approach to data assimilation appears efficient in separating individual components of complex response with statistical analysis unachievable by traditional data processing methods.
Aristolochic acid (AA) is a natural compound found in many plants of the Aristolochia genus, and these plants are widely used in traditional medicines for numerous conditions and for weight loss. Previous work has connected AA-mutagenesis to upper-tract urothelial cell carcinomas and hepatocellular carcinomas. We hypothesize that AA may also contribute to bladder cancer.
Here, we investigated the involvement of AA-mutagenesis in bladder cancer by sequencing bladder tumor genomes from two patients with known exposure to AA. After detecting strong mutational signatures of AA exposure in these tumors, we exome-sequenced and analyzed an additional 11 bladder tumors and analyzed publicly available somatic mutation data from a further 336 bladder tumors.
The somatic mutations in the bladder tumors from the two patients with known AA exposure showed overwhelming AA signatures. We also detected evidence of AA exposure in 1 out of 11 bladder tumors from Singapore and in 3 out of 99 bladder tumors from China. In addition, 1 out of 194 bladder tumors from North America showed a pattern of mutations that might have resulted from exposure to an unknown mutagen with a heretofore undescribed pattern of A > T mutations. Besides the signature of AA exposure, the bladder tumors also showed the CpG > TpG and activated-APOBEC signatures, which have been previously reported in bladder cancer.
This study demonstrates the utility of inferring mutagenic exposures from somatic mutation spectra. Moreover, AA exposure in bladder cancer appears to be more pervasive in the East, where traditional herbal medicine is more widely used. More broadly, our results suggest that AA exposure is more extensive than previously thought both in terms of populations at risk and in terms of types of cancers involved. This appears to be an important public health issue that should be addressed by further investigation and by primary prevention through regulation and education. In addition to opportunities for primary prevention, knowledge of AA exposure would provide opportunities for secondary prevention in the form of intensified screening of patients with known or suspected AA exposure.
Electronic supplementary material
The online version of this article (doi:10.1186/s13073-015-0161-3) contains supplementary material, which is available to authorized users.
Somatic mutations of TP53 are among the most common in cancer and germline mutations of TP53 (usually missense) can cause Li-Fraumeni syndrome (LFS). Recently, recurrent genomic rearrangements in intron 1 of TP53 have been described in osteosarcoma (OS), a highly malignant neoplasm of bone belonging to the spectrum of LFS tumors. Using whole-genome sequencing of OS, we found features of TP53 intron 1 rearrangements suggesting a unique mechanism correlated with transcription. Screening of 288 OS and 1,090 tumors of other types revealed evidence for TP53 rearrangements in 46 (16%) OS, while none were detected in other tumor types, indicating this rearrangement to be highly specific to OS. We revisited a four-generation LFS family where no TP53 mutation had been identified and found a 445 kb inversion spanning from the TP53 intron 1 towards the centromere. The inversion segregated with tumors in the LFS family. Cancers in this family had loss of heterozygosity, retaining the rearranged allele and resulting in TP53 expression loss. In conclusion, intron 1 rearrangements cause p53-driven malignancies by both germline and somatic mechanisms and provide an important mechanism of TP53 inactivation in LFS, which might in part explain the diagnostic gap of formerly classified “TP53 wild-type” LFS.
TP53; Li-Fraumeni syndrome; osteosarcoma; cancer genomics; structural variations