Transcriptomic analysis of differential host gene expression upon uptake of symbionts: a case study with Symbiodinium and the major bioeroding sponge Cliona varians

Background: We have a limited understanding of genomic interactions that occur among partners for many symbioses. One of the most important symbioses in tropical reef habitats involves Symbiodinium. Most work examining Symbiodinium-host interactions involves cnidarian partners. To fully and broadly understand the conditions that permit Symbiodinium to procure intracellular residency, we must explore hosts from different taxa to help uncover universal cellular and genetic strategies for invading and persisting in host cells. Here, we present data from gene expression analyses involving the bioeroding sponge Cliona varians that harbors Clade G Symbiodinium. Results: Patterns of differential gene expression from distinct symbiont states (“normal”, “reinfected”, and “aposymbiotic”) of the sponge host are presented based on two comparative approaches (transcriptome sequencing and suppressive subtractive hybridization (SSH)). Transcriptomic profiles were different when reinfected tissue was compared to normal and aposymbiotic tissue. We characterized a set of 40 genes drawn from a pool of differentially expressed genes in “reinfected” tissue compared to “aposymbiotic” tissue via SSH. As proof of concept, we determined whether some of the differentially expressed genes identified above could be monitored in sponges grown under ecologically realistic field conditions. We allowed aposymbiotic sponge tissue to become re-populated by natural pools of Symbiodinium in shallow water flats in the Florida Keys, and we analyzed gene expression profiles for two genes found to be increased in expression in “reinfected” tissue in both the transcriptome and via SSH. These experiments highlighted the experimental tractability of C. varians to explore with precision the genetic events that occur upon establishment of the symbiosis. We briefly discuss lab- and field-based experimental approaches that promise to offer insights into the co-opted genetic networks that may modulate uptake and regulation of Symbiondinium populations in hospite. Conclusions: This work provides a sponge transcriptome, and a database of putative genes and genetic pathways that may be involved in Symbiodinium interactions. The relative patterns of gene expression observed in these experiments will need to be evaluated on a gene-by-gene basis in controlled and natural re-infection experiments. We argue that sponges offer particularly useful characteristics for discerning essential dimensions of the Symbiodinium niche. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-376) contains supplementary material, which is available to authorized users

A single-nucleotide polymorphism in the methylene tetrahydrofolate reductase (MTHFR) gene is associated with risk of radiation pneumonitis in lung cancer patients treated with thoracic radiation therapy

Background: To study the association between functional single nucleotide polymorphisms (SNPs) in candidate genes from oxidative stress pathways, and risk of radiation pneumonitis (RP) in patients treated with thoracic radiation therapy (RT) for locally advanced lung cancer (LC). Methods: We reviewed 136 patients treated with RT for LC between 2001 and 2007, and had prior genotyping of functional SNPs in oxidative stress genes including superoxide dismutase 2 (SOD2; rs4880) and methylenetetrahydrofolate reductase (MTHFR; rs1801131, rs1801133). RP events were retrospectively scored using the Common Terminology Criteria for Adverse Events, version 4.0. Cox proportional hazard regression was performed to identify clinical variables and genotypes associated with risk of grade ≥2 and grade ≥3 RP on univariate and multivariate analysis. P-values were corrected for multiple hypothesis testing. Results: With a median follow-up of 21.4 months, the incidence of ≥grade 2 RP was 29% and ≥grade 3 RP was 14%. On multivariate analysis, after adjusting for clinical factors such as concurrent chemotherapy, and consolidation docetaxel, and lung dosimetric parameters such as V20 and mean lung dose, MTHFR genotype (rs1801131; AA versus AC/CC) was significantly associated with risk of ≥grade 2 RP (Hazard ratio [HR]: 0.37; 95% confidence interval [CI]: 0.18-0.76; p=0.006, corrected p=0.018) and ≥grade 3 RP (HR: 0.21; 95% CI: 0.06-0.70; p=0.01; corrected p=0.03). SOD2 genotype was not associated with RP. Conclusions: Our study showed an association between MTHFR genotype and risk of clinically significant RP. Further study of MTHFR-related pathways may provide insight into the mechanisms behind RP

Acquired Resistance to KRAS (G12C) Inhibition in Cancer

BACKGROUND: Clinical trials of the KRAS inhibitors adagrasib and sotorasib have shown promising activity in cancers harboring KRAS glycine-to-cysteine amino acid substitutions at codon 12 (KRAS(G12C)). The mechanisms of acquired resistance to these therapies are currently unknown. METHODS: Among patients with KRAS(G12C) -mutant cancers treated with adagrasib monotherapy, we performed genomic and histologic analyses that compared pretreatment samples with those obtained after the development of resistance. Cell-based experiments were conducted to study mutations that confer resistance to KRAS(G12C) inhibitors. RESULTS: A total of 38 patients were included in this study: 27 with non-small-cell lung cancer, 10 with colorectal cancer, and 1 with appendiceal cancer. Putative mechanisms of resistance to adagrasib were detected in 17 patients (45% of the cohort), of whom 7 (18% of the cohort) had multiple coincident mechanisms. Acquired KRAS alterations included G12D/R/V/W, G13D, Q61H, R68S, H95D/Q/R, Y96C, and high-level amplification of the KRAS(G12C) allele. Acquired bypass mechanisms of resistance included MET amplification; activating mutations in NRAS, BRAF, MAP2K1, and RET; oncogenic fusions involving ALK, RET, BRAF, RAF1, and FGFR3; and loss-of-function mutations in NF1 and PTEN. In two of nine patients with lung adenocarcinoma for whom paired tissue-biopsy samples were available, histologic transformation to squamous-cell carcinoma was observed without identification of any other resistance mechanisms. Using an in vitro deep mutational scanning screen, we systematically defined the landscape of KRAS mutations that confer resistance to KRAS(G12C) inhibitors. CONCLUSIONS: Diverse genomic and histologic mechanisms impart resistance to covalent KRAS(G12C) inhibitors, and new therapeutic strategies are required to delay and overcome this drug resistance in patients with cancer. (Funded by Mirati Therapeutics and others; ClinicalTrials.gov number, NCT03785249.)

DW-MRI as a Biomarker to Compare Therapeutic Outcomes in Radiotherapy Regimens Incorporating Temozolomide or Gemcitabine in Glioblastoma

The effectiveness of the radiosensitizer gemcitabine (GEM) was evaluated in a mouse glioma along with the imaging biomarker diffusion-weighted magnetic resonance imaging (DW-MRI) for early detection of treatment effects. A genetically engineered murine GBM model [Ink4a-Arf−/− PtenloxP/loxP/Ntv-a RCAS/PDGF(+)/Cre(+)] was treated with gemcitabine (GEM), temozolomide (TMZ) +/− ionizing radiation (IR). Therapeutic efficacy was quantified by contrast-enhanced MRI and DW-MRI for growth rate and tumor cellularity, respectively. Mice treated with GEM, TMZ and radiation showed a significant reduction in growth rates as early as three days post-treatment initiation. Both combination treatments (GEM/IR and TMZ/IR) resulted in improved survival over single therapies. Tumor diffusion values increased prior to detectable changes in tumor volume growth rates following administration of therapies. Concomitant GEM/IR and TMZ/IR was active and well tolerated in this GBM model and similarly prolonged median survival of tumor bearing mice. DW-MRI provided early changes to radiosensitization treatment warranting evaluation of this imaging biomarker in clinical trials

Diffusion-weighted MRI as a biomarker of tumor radiation treatment response heterogeneity: A comparative study of whole-volume histogram analysis versus voxel-based functional diffusion map analysis

RATIONALE: Treatment of glioblastoma (GBM) remains challenging due in part to its histologic intratumoral heterogeneity that contributes to its overall poor treatment response. Our goal was to evaluate a voxel-based biomarker, the functional diffusion map (fDM), as an imaging biomarker to detect heterogeneity of tumor response in a radiation dose escalation protocol using a genetically engineered murine GBM model. EXPERIMENTAL DESIGN: Twenty-four genetically engineered murine GBM models [Ink4a-Arf-/-/Ptenloxp/loxp/Ntv-a RCAS/PDGF(+)/Cre(+)] were randomized in four treatment groups (n = 6 per group) consisting of daily doses of 0, 1, 2, and 4 Gy delivered for 5 days. Contrast-enhanced T1-weighted and diffusion-weighted magnetic resonance imaging (MRI) scans were acquired for tumor delineation and quantification of apparent diffusion coefficient (ADC) maps, respectively. MRI experiments were performed daily for a week and every 2 days thereafter. For each animal, the area under the curve (AUC) of the percentage change of the ADC (AUCADC) and that of the increase in fDM values (AUCfDM+) were determined within the first 5 days following therapy initiation. RESULTS: Animal survival increased with increasing radiation dose. Treatment induced a dose-dependent increase in tumor ADC values. The strongest correlation between survival and ADC measurements was observed using the AUCfDM+ metric (R2 = 0.88). CONCLUSION: This study showed that the efficacy of a voxel-based imaging biomarker (fDM) was able to detect spatially varying changes in tumors, which were determined to be a more sensitive predictor of overall response versus whole-volume tumor measurements (AUCADC). Finally, fDM provided for visualization of treatment-associated spatial heterogeneity within the tumor. © 2013 Neoplasia Press, Inc. All rights reserved

Evaluation of Concurrent Radiation, Temozolomide and ABT-888 Treatment Followed by Maintenance Therapy with Temozolomide and ABT-888 in a Genetically Engineered Glioblastoma Mouse Model

Despite the use of ionizing radiation (IR) and temozolomide (TMZ), outcome for glioblastoma (GBM) patients remains dismal. Poly (ADP-ribose) polymerase (PARP) is important in repair pathways for IR-induced DNA damage and TMZ-induced alkylation at N7-methylguanine and N3-methyldenine. However, optimized protocols for administration of PARP inhibitors have not been delineated. In this study, the PARP inhibitor ABT-888 was evaluated in combination with and compared to current standard-of-care in a genetically engineered mouse GBM model. Results demonstrated that concomitant TMZ/IR/ABT-888 with adjuvant TMZ/ABT-888 was more effective in inducing apoptosis and reducing proliferation with significant tumor growth delay and improved overall survival over concomitant TMZ/IR with adjuvant TMZ. Diffusion-weighted MRI, an early translatable response biomarker detected changes in tumors reflecting response at 1 day post TMZ/IR/ABT-888 treatment. This study provides strong scientific rationale for the development of an optimized dosing regimen for a PARP inhibitor with TMZ/IR for upfront treatment of GBM

MRI-Guided Stereotactic Biopsy of Murine GBM for Spatiotemporal Molecular Genomic Assessment

Brain tumor biopsies that are routinely performed in clinical settings significantly aid in diagnosis and staging. The aim of this study is to develop and evaluate a methodological image-guided approach that would allow for routine sampling of glioma tissue from orthotopic mouse brain tumor models. A magnetic resonance imaging-guided biopsy method is presented to allow for spatially precise stereotaxic sampling of a murine glioma coupled with genome-scale technology to provide unbiased characterization of intra- and intertumoral clonal heterogeneity. Longitudinal and multiregional sampling of intracranial tumors allows for successful collection of tumor biopsy samples, thus allowing for a pathway-enrichment analysis and a transcriptional profiling of RNA sequencing data. Spatiotemporal gene expression pattern variations revealing genomic heterogeneity were found