Genomic heterogeneity as a barrier to precision oncology in urothelial cancer.

Precision oncology relies on the accurate molecular characterization of individual patients with cancer at the time of treatment initiation. However, tumor molecular profiles are not static, and cancers continually evolve because of ongoing mutagenesis and clonal selection. Here, we performed genomic analyses of primary tumors, metastases, and plasma collected from individual patients to define the concordance of actionable genomic alterations and to identify drivers of metastatic disease progression. We observed a high degree of discordance of actionable genomic alterations, with 23% discordant between primary and metastatic disease sites. Among chromatin-modifying genes, ARID1A mutations, when discordant, were exclusive to the metastatic tumor samples. Our findings indicate that the high degree of lesion-to-lesion genomic heterogeneity may be a barrier to precision oncology approaches for bladder cancer and that circulating tumor DNA profiling may be preferred to tumor sequencing for a subset of patients.

Identification of a Synthetic Lethal Relationship between Nucleotide Excision Repair Deficiency and Irofulven Sensitivity in Urothelial Cancer.

PURPOSE: Cisplatin-based chemotherapy is a first-line treatment for muscle-invasive and metastatic urothelial cancer. Approximately 10% of bladder urothelial tumors have a somatic missense mutation in the nucleotide excision repair (NER) gene, ERCC2, which confers increased sensitivity to cisplatin-based chemotherapy. However, a significant subset of patients is ineligible to receive cisplatin-based therapy due to medical contraindications, and no NER-targeted approaches are available for platinum-ineligible or platinum-refractory ERCC2-mutant cases. EXPERIMENTAL DESIGN: We used a series of NER-proficient and NER-deficient preclinical tumor models to test sensitivity to irofulven, an abandoned anticancer agent. In addition, we used available clinical and sequencing data from multiple urothelial tumor cohorts to develop and validate a composite mutational signature of ERCC2 deficiency and cisplatin sensitivity. RESULTS: We identified a novel synthetic lethal relationship between tumor NER deficiency and sensitivity to irofulven. Irofulven specifically targets cells with inactivation of the transcription-coupled NER (TC-NER) pathway and leads to robust responses in vitro and in vivo, including in models with acquired cisplatin resistance, while having minimal effect on cells with intact NER. We also found that a composite mutational signature of ERCC2 deficiency was strongly associated with cisplatin response in patients and was also associated with cisplatin and irofulven sensitivity in preclinical models. CONCLUSIONS: Tumor NER deficiency confers sensitivity to irofulven, a previously abandoned anticancer agent, with minimal activity in NER-proficient cells. A composite mutational signature of NER deficiency may be useful in identifying patients likely to respond to NER-targeting agents, including cisplatin and irofulven.See related commentary by Jiang and Greenberg, p. 1833.

AKT1 E17K Inhibits Cancer Cell Migration by Abrogating ß-Catenin Signaling.

Mutational activation of the PI3K/AKT pathway is among the most common pro-oncogenic events in human cancers. The clinical utility of PI3K and AKT inhibitors has, however, been modest to date. Here, we used CRISPR-mediated gene editing to study the biological consequences of AKT1 E17K mutation by developing an AKT1 E17K-mutant isogenic system in a TP53-null background. AKT1 E17K expression under the control of its endogenous promoter enhanced cell growth and colony formation, but had a paradoxical inhibitory effect on cell migration and invasion. The mechanistic basis by which activated AKT1 inhibited cell migration and invasion was increased E-cadherin expression mediated by suppression of ZEB1 transcription via altered ß-catenin subcellular localization. This phenotypic effect was AKT1-specific, as AKT2 activation had the opposite effect, a reduction in E-cadherin expression. Consistent with the opposing effects of AKT1 and AKT2 activation on E-cadherin expression, a pro-migratory effect of AKT1 activation was not observed in breast cancer cells with PTEN loss or expression of an activating PIK3CA mutation, alterations which induce the activation of both AKT isoforms. The results suggest that the use of AKT inhibitors in patients with breast cancer could paradoxically accelerate metastatic progression in some genetic contexts and may explain the frequent coselection for CDH1 mutations in AKT1-mutated breast tumors. IMPLICATIONS: AKT1 E17K mutation in breast cancer impairs migration/invasiveness via sequestration of ß-catenin to the cell membrane leading to decreased ZEB1 transcription, resulting in increased E-cadherin expression and a reversal of epithelial-mesenchymal transition.

Activation of Nrf2 in the liver is associated with stress resistance mediated by suppression of the growth hormone-regulated STAT5b transcription factor.

The transcription factor Nrf2 (encoded by Nfe2l2) induces expression of numerous detoxifying and antioxidant genes in response to oxidative stress. The cytoplasmic protein Keap1 interacts with and represses Nrf2 function. Computational approaches were developed to identify factors that modulate Nrf2 in a mouse liver gene expression compendium. Forty-eight Nrf2 biomarker genes were identified using profiles from the livers of mice in which Nrf2 was activated genetically in Keap1-null mice or chemically by a potent activator of Nrf2 signaling. The rank-based Running Fisher statistical test was used to determine the correlation between the Nrf2 biomarker genes and a test set of 81 profiles with known Nrf2 activation status demonstrating a balanced accuracy of 96%. For a large number of factors examined in the compendium, we found consistent relationships between activation of Nrf2 and feminization of the liver transcriptome through suppression of the male-specific growth hormone (GH)-regulated transcription factor STAT5b. The livers of female mice exhibited higher Nrf2 activation than male mice in untreated or chemical-treated conditions. In male mice, Nrf2 was activated by treatment with ethinyl estradiol, whereas in female mice, Nrf2 was suppressed by treatment with testosterone. Nrf2 was activated in 5 models of disrupted GH signaling containing mutations in Pit1, Prop1, Ghrh, Ghrhr, and Ghr. Out of 59 chemical treatments that activated Nrf2, 36 exhibited STAT5b suppression in the male liver. The Nrf2-STAT5b coupling was absent in in vitro comparisons of chemical treatments. Treatment of male and female mice with 11 chemicals that induce oxidative stress led to activation of Nrf2 to greater extents in females than males. The enhanced basal and inducible levels of Nrf2 activation in females relative to males provides a molecular explanation for the greater resistance often seen in females vs. males to age-dependent diseases and chemical-induced toxicity.

Accumulation of medium-chain, saturated fatty acyl moieties in seed oils of transgenic Camelina sativa.

With its high seed oil content, the mustard family plant Camelina sativa has gained attention as a potential biofuel source. As a bioenergy crop, camelina has many advantages. It grows on marginal land with low demand for water and fertilizer, has a relatively short life cycle, and is stress tolerant. As most other crop seed oils, camelina seed triacylglycerols (TAGs) consist of mostly long, unsaturated fatty acyl moieties, which is not desirable for biofuel processing. In our efforts to produce shorter, saturated chain fatty acyl moieties in camelina seed oil for conversion to jet fuel, a 12:0-acyl-carrier thioesterase gene, UcFATB1, from California bay (Umbellularia californica Nutt.) was expressed in camelina seeds. Up to 40% of short chain laurate (C12:0) and myristate (C14:0) were present in TAGs of the seed oil of the transgenics. The total oil content and germination rate of the transgenic seeds were not affected. Analysis of positions of these two fatty acyl moieties in TAGs indicated that they were present at the sn-1 and sn-3 positions, but not sn-2, on the TAGs. Suppression of the camelina KASII genes by RNAi constructs led to higher accumulation of palmitate (C16:0), from 7.5% up to 28.5%, and further reduction of longer, unsaturated fatty acids in seed TAGs. Co-transformation of camelina with both constructs resulted in enhanced accumulation of all three medium-chain, saturated fatty acids in camelina seed oils. Our results show that a California bay gene can be successfully used to modify the oil composition in camelina seed and present a new biological alternative for jet fuel production.

Disruption of STAT5b-Regulated Sexual Dimorphism of the Liver Transcriptome by Diverse Factors Is a Common Event.

Signal transducer and activator of transcription 5b (STAT5b) is a growth hormone (GH)-activated transcription factor and a master regulator of sexually dimorphic gene expression in the liver. Disruption of the GH hypothalamo-pituitary-liver axis controlling STAT5b activation can lead to metabolic dysregulation, steatosis, and liver cancer. Computational approaches were developed to identify factors that disrupt STAT5b function in a mouse liver gene expression compendium. A biomarker comprised of 144 STAT5b-dependent genes was derived using comparisons between wild-type male and wild-type female mice and between STAT5b-null and wild-type mice. Correlations between the STAT5b biomarker gene set and a test set comprised of expression datasets (biosets) with known effects on STAT5b function were evaluated using a rank-based test (the Running Fisher algorithm). Using a similarity p-value ≤ 10(-4), the test achieved a balanced accuracy of 99% and 97% for detection of STAT5b activation or STAT5b suppression, respectively. The STAT5b biomarker gene set was then used to identify factors that activate (masculinize) or suppress (feminize) STAT5b function in an annotated mouse liver and primary hepatocyte gene expression compendium of ~1,850 datasets. Disruption of GH-regulated STAT5b is a common phenomenon in liver in vivo, with 5% and 29% of the male datasets, and 11% and 13% of the female datasets, associated with masculinization or feminization, respectively. As expected, liver STAT5b activation/masculinization occurred at puberty and suppression/feminization occurred during aging and in mutant mice with defects in GH signaling. A total of 70 genes were identified that have effects on STAT5b activation in genetic models in which the gene was inactivated or overexpressed. Other factors that affected liver STAT5b function were shown to include fasting, caloric restriction and infections. Together, these findings identify diverse factors that perturb the hypothalamo-pituitary-liver GH axis and disrupt GH-dependent STAT5b activation in mouse liver.

A photorespiratory bypass increases plant growth and seed yield in biofuel crop Camelina sativa.

BACKGROUND: Camelina sativa is an oilseed crop with great potential for biofuel production on marginal land. The seed oil from camelina has been converted to jet fuel and improved fuel efficiency in commercial and military test flights. Hydrogenation-derived renewable diesel from camelina is environmentally superior to that from canola due to lower agricultural inputs, and the seed meal is FDA approved for animal consumption. However, relatively low yield makes its farming less profitable. Our study is aimed at increasing camelina seed yield by reducing carbon loss from photorespiration via a photorespiratory bypass. Genes encoding three enzymes of the Escherichia coli glycolate catabolic pathway were introduced: glycolate dehydrogenase (GDH), glyoxylate carboxyligase (GCL) and tartronic semialdehyde reductase (TSR). These enzymes compete for the photorespiratory substrate, glycolate, convert it to glycerate within the chloroplasts, and reduce photorespiration. As a by-product of the reaction, CO2 is released in the chloroplast, which increases photosynthesis. Camelina plants were transformed with either partial bypass (GDH), or full bypass (GDH, GCL and TSR) genes. Transgenic plants were evaluated for physiological and metabolic traits. RESULTS: Expressing the photorespiratory bypass genes in camelina reduced photorespiration and increased photosynthesis in both partial and full bypass expressing lines. Expression of partial bypass increased seed yield by 50-57 %, while expression of full bypass increased seed yield by 57-73 %, with no loss in seed quality. The transgenic plants also showed increased vegetative biomass and faster development; they flowered, set seed and reached seed maturity about 1 week earlier than WT. At the transcriptional level, transgenic plants showed differential expression in categories such as respiration, amino acid biosynthesis and fatty acid metabolism. The increased growth of the bypass transgenics compared to WT was only observed in ambient or low CO2 conditions, but not in elevated CO2 conditions. CONCLUSIONS: The photorespiratory bypass is an effective approach to increase photosynthetic productivity in camelina. By reducing photorespiratory losses and increasing photosynthetic CO2 fixation rates, transgenic plants show dramatic increases in seed yield. Because photorespiration causes losses in productivity of most C3 plants, the bypass approach may have significant impact on increasing agricultural productivity for C3 crops.

Screening a mouse liver gene expression compendium identifies modulators of the aryl hydrocarbon receptor (AhR).

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that mediates the biological and toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), dioxin-like compounds (DLC) as well as some drugs and endogenous tryptophan metabolites. Short-term activation of AhR can lead to hepatocellular steatosis, and chronic activation can lead to liver cancer in mice and rats. Analytical approaches were developed to identify biosets in a genomic database in which AhR activity was altered. A set of 63 genes was identified (the AhR gene expression biomarker) that was dependent on AhR for regulation after exposure to TCDD or benzo[a]pyrene and includes the known AhR targets Cyp1a1 and Cyp1b1. A fold-change rank-based test (Running Fisher's test; p-value ≤ 10(-4)) was used to evaluate the similarity between the AhR biomarker and a test set of 37 and 41 biosets positive or negative, respectively for AhR activation. The test resulted in a balanced accuracy of 95%. The rank-based test was used to identify factors that activate or suppress AhR in an annotated mouse liver/mouse primary hepatocyte gene expression database of ∼ 1850 comparisons. In addition to the expected activation of AhR by TCDD and DLC, AhR was activated by AP20189 and phenformin. AhR was suppressed by phenobarbital and 1,4-Bis[2-(3,5-dichloropyridyloxy)] benzene (TCPOBOP) in a constitutive activated receptor (CAR)-dependent manner and pregnenolone-16α-carbonitrile in a pregnane X receptor (PXR)-dependent manner. Inactivation of individual genes in nullizygous models led to AhR activation (Pxr, Ghrhr, Taf10) or suppression (Ahr, Ilst6st, Hnf1a). This study describes a novel screening strategy for identifying factors in mouse liver that perturb AhR in a gene expression compendium.

Identification of chemical modulators of the constitutive activated receptor (CAR) in a gene expression compendium.

The nuclear receptor family member constitutive activated receptor (CAR) is activated by structurally diverse drugs and environmentally-relevant chemicals leading to transcriptional regulation of genes involved in xenobiotic metabolism and transport. Chronic activation of CAR increases liver cancer incidence in rodents, whereas suppression of CAR can lead to steatosis and insulin insensitivity. Here, analytical methods were developed to screen for chemical treatments in a gene expression compendium that lead to alteration of CAR activity. A gene expression biomarker signature of 83 CAR-dependent genes was identified using microarray profiles from the livers of wild-type and CAR-null mice after exposure to three structurally-diverse CAR activators (CITCO, phenobarbital, TCPOBOP). A rank-based algorithm (Running Fisher's algorithm (p-value ≤ 10(-4))) was used to evaluate the similarity between the CAR biomarker signature and a test set of 28 and 32 comparisons positive or negative, respectively, for CAR activation; the test resulted in a balanced accuracy of 97%. The biomarker signature was used to identify chemicals that activate or suppress CAR in an annotated mouse liver/primary hepatocyte gene expression database of ~1850 comparisons. CAR was activated by 1) activators of the aryl hydrocarbon receptor (AhR) in wild-type but not AhR-null mice, 2) pregnane X receptor (PXR) activators in wild-type and to lesser extents in PXR-null mice, and 3) activators of PPARα in wild-type and PPARα-null mice. CAR was consistently activated by five conazole fungicides and four perfluorinated compounds. Comparison of effects in wild-type and CAR-null mice showed that the fungicide propiconazole increased liver weight and hepatocyte proliferation in a CAR-dependent manner, whereas the perfluorinated compound perfluorooctanoic acid (PFOA) increased these endpoints in a CAR-independent manner. A number of compounds suppressed CAR coincident with increases in markers of inflammation including acetaminophen, concanavalin A, lipopolysaccharide, and 300 nm silica particles. In conclusion, we have shown that a CAR biomarker signature coupled with a rank-based similarity method accurately predicts CAR activation. This analytical approach, when applied to a gene expression compendium, increased the universe of known chemicals that directly or indirectly activate CAR, highlighting the promiscuous nature of CAR activation and signaling through activation of other xenobiotic-activated receptors.

Methods for RNA profiling of gravi-responding plant tissues.

Plant transcriptional responses to gravity stimulation by reorientation are among the fastest measured in any tissue or species. Upon reorientation, changes in abundance of specific mRNAs can be measured within seconds or minutes, for plastid or nuclear encoded genes, respectively. Identifying fast gravity-induced transcripts has been made possible by the development of high-throughput technology for qualitative and quantitative RNA analysis. RNA profiling has undergone further rapid development due to its enormous potential in basic sciences and medical applications. We describe here the current and most widely used methods to profile the changes in an entire transcriptome by high-throughput sequencing of RNA fractions (RNAseq) and single gene transcript analysis using real-time quantitative reverse-transcriptase polymerase chain reaction (qRT-PCR).

Identification of modulators of the nuclear receptor peroxisome proliferator-activated receptor α (PPARα) in a mouse liver gene expression compendium.

The nuclear receptor family member peroxisome proliferator-activated receptor α (PPARα) is activated by therapeutic hypolipidemic drugs and environmentally-relevant chemicals to regulate genes involved in lipid transport and catabolism. Chronic activation of PPARα in rodents increases liver cancer incidence, whereas suppression of PPARα activity leads to hepatocellular steatosis. Analytical approaches were developed to identify biosets (i.e., gene expression differences between two conditions) in a genomic database in which PPARα activity was altered. A gene expression signature of 131 PPARα-dependent genes was built using microarray profiles from the livers of wild-type and PPARα-null mice after exposure to three structurally diverse PPARα activators (WY-14,643, fenofibrate and perfluorohexane sulfonate). A fold-change rank-based test (Running Fisher's test (p-value ≤ 10(-4))) was used to evaluate the similarity between the PPARα signature and a test set of 48 and 31 biosets positive or negative, respectively for PPARα activation; the test resulted in a balanced accuracy of 98%. The signature was then used to identify factors that activate or suppress PPARα in an annotated mouse liver/primary hepatocyte gene expression compendium of ~1850 biosets. In addition to the expected activation of PPARα by fibrate drugs, di(2-ethylhexyl) phthalate, and perfluorinated compounds, PPARα was activated by benzofuran, galactosamine, and TCDD and suppressed by hepatotoxins acetaminophen, lipopolysaccharide, silicon dioxide nanoparticles, and trovafloxacin. Additional factors that activate (fasting, caloric restriction) or suppress (infections) PPARα were also identified. This study 1) developed methods useful for future screening of environmental chemicals, 2) identified chemicals that activate or suppress PPARα, and 3) identified factors including diets and infections that modulate PPARα activity and would be hypothesized to affect chemical-induced PPARα activity.