Molecular subtypes of high-grade serous ovarian cancer across racial groups and gene expression platforms.

BACKGROUND: High-grade serous carcinoma (HGSC) gene expression subtypes are associated with differential survival. We characterized HGSC gene expression in Black individuals and considered whether gene expression differences by self-identified race may contribute to poorer HGSC survival among Black versus White individuals. METHODS: We included newly generated RNA-Seq data from Black and White individuals, and array-based genotyping data from four existing studies of White and Japanese individuals. We used K-means clustering, a method with no predefined number of clusters or dataset-specific features, to assign subtypes. Cluster- and dataset-specific gene expression patterns were summarized by moderated t-scores. We compared cluster-specific gene expression patterns across datasets by calculating the correlation between the summarized vectors of moderated t-scores. Following mapping to The Cancer Genome Atlas (TCGA)-derived HGSC subtypes, we used Cox proportional hazards models to estimate subtype-specific survival by dataset. RESULTS: Cluster-specific gene expression was similar across gene expression platforms and racial groups. Comparing the Black population to the White and Japanese populations, the immunoreactive subtype was more common (39% versus 23%-28%) and the differentiated subtype less common (7% versus 22%-31%). Patterns of subtype-specific survival were similar between the Black and White populations with RNA-Seq data; compared to mesenchymal cases, the risk of death was similar for proliferative and differentiated cases and suggestively lower for immunoreactive cases (Black population HR=0.79 [0.55, 1.13], White population HR=0.86 [0.62, 1.19]). CONCLUSIONS: While the prevalence of HGSC subtypes varied by race, subtype-specific survival was similar. IMPACT: HGSC subtypes can be consistently assigned across platforms and self-identified racial groups.

Molecular subtypes of high-grade serous ovarian cancer across racial groups and gene expression platforms.

Introduction: High-grade serous carcinoma (HGSC) gene expression subtypes are associated with differential survival. We characterized HGSC gene expression in Black individuals and considered whether gene expression differences by race may contribute to poorer HGSC survival among Black versus non-Hispanic White individuals. Methods: We included newly generated RNA-Seq data from Black and White individuals, and array-based genotyping data from four existing studies of White and Japanese individuals. We assigned subtypes using K-means clustering. Cluster- and dataset-specific gene expression patterns were summarized by moderated t-scores. We compared cluster-specific gene expression patterns across datasets by calculating the correlation between the summarized vectors of moderated t-scores. Following mapping to The Cancer Genome Atlas (TCGA)-derived HGSC subtypes, we used Cox proportional hazards models to estimate subtype-specific survival by dataset. Results: Cluster-specific gene expression was similar across gene expression platforms. Comparing the Black study population to the White and Japanese study populations, the immunoreactive subtype was more common (39% versus 23%-28%) and the differentiated subtype less common (7% versus 22%-31%). Patterns of subtype-specific survival were similar between the Black and White populations with RNA-Seq data; compared to mesenchymal cases, the risk of death was similar for proliferative and differentiated cases and suggestively lower for immunoreactive cases (Black population HR=0.79 [0.55, 1.13], White population HR=0.86 [0.62, 1.19]). Conclusions: A single, platform-agnostic pipeline can be used to assign HGSC gene expression subtypes. While the observed prevalence of HGSC subtypes varied by race, subtype-specific survival was similar.

Transgenic Expression of Nrf2 Induces a Pro-Reductive Stress and Adaptive Cardiac Remodeling in the Mouse.

Nuclear factor, erythroid 2 like 2 (Nfe2l2 or Nrf2), is a transcription factor that protects cells by maintaining a homeostatic redox state during stress. The constitutive expression of Nrf2 (CaNrf2-TG) was previously shown to be pathological to the heart over time. We tested a hypothesis that the cardiac-specific expression of full length Nrf2 (mNrf2-TG) would moderately increase the basal antioxidant defense, triggering a pro-reductive environment leading to adaptive cardiac remodeling. Transgenic and non-transgenic (NTG) mice at 7−8 months of age were used to analyze the myocardial transcriptome, structure, and function. Next generation sequencing (NGS) for RNA profiling and qPCR-based validation of the NGS data, myocardial redox levels, and imaging (echocardiography) were performed. Transcriptomic analysis revealed that out of 14,665 identified mRNAs, 680 were differently expressed (DEG) in TG hearts. Of 680 DEGs, 429 were upregulated and 251 were downregulated significantly (FC > 2.0, p < 0.05). Gene set enrichment analysis revealed that the top altered pathways were (a) Nrf2 signaling, (b) glutathione metabolism and (c) ROS scavenging. A comparative analysis of the glutathione redox state in the hearts demonstrated significant differences between pro-reductive vs. hyper-reductive conditions (233 ± 36.7 and 380 ± 68.7 vs. 139 ± 8.6 µM/mg protein in mNrf2-TG and CaNrf2-TG vs. NTG). Genes involved in fetal development, hypertrophy, cytoskeletal rearrangement, histone deacetylases (HDACs), and GATA transcription factors were moderately increased in mNrf2-TG compared to CaNrf2-TG. Non-invasive echocardiography analysis revealed an increase in systolic function (ejection fraction) in mNrf2-TG, suggesting an adaptation, as opposed to pathological remodeling in CaNrf2-TG mice experiencing a hyper-reductive stress, leading to reduced survival (40% at 60 weeks). The effects of excess Nrf2-driven antioxidant transcriptome revealed a pro-reductive condition in the myocardium leading to an adaptive cardiac remodeling. While pre-conditioning the myocardial redox with excess antioxidants (i.e., pro-reductive state) could be beneficial against oxidative stress, a chronic pro-reductive environment in the myocardium might transition the adaptation to pathological remodeling.

Identification of Nrf2-responsive microRNA networks as putative mediators of myocardial reductive stress.

Although recent advances in the treatment of acute coronary heart disease have reduced mortality rates, few therapeutic strategies exist to mitigate the progressive loss of cardiac function that manifests as heart failure. Nuclear factor, erythroid 2 like 2 (Nfe2l2, Nrf2) is a transcriptional regulator that is known to confer transient myocardial cytoprotection following acute ischemic insult; however, its sustained activation paradoxically causes a reductive environment characterized by excessive antioxidant activity. We previously identified a subset of 16 microRNAs (miRNA) significantly diminished in Nrf2-ablated (Nrf2-/-) mouse hearts, leading to the hypothesis that increasing levels of Nrf2 activation augments miRNA induction and post-transcriptional dysregulation. Here, we report the identification of distinct miRNA signatures (i.e. "reductomiRs") associated with Nrf2 overexpression in a cardiac-specific and constitutively active Nrf2 transgenic (caNrf2-Tg) mice expressing low (TgL) and high (TgH) levels. We also found several Nrf2 dose-responsive miRNAs harboring proximal antioxidant response elements (AREs), implicating these "reductomiRs" as putative meditators of Nrf2-dependent post-transcriptional regulation. Analysis of mRNA-sequencing identified a complex network of miRNAs and effector mRNAs encoding known pathological hallmarks of cardiac stress-response. Altogether, these data support Nrf2 as a putative regulator of cardiac miRNA expression and provide novel candidates for future mechanistic investigation to understand the relationship between myocardial reductive stress and cardiac pathophysiology.

Histone Deacetylase Inhibition Has Targeted Clinical Benefit in ARID1A-Mutated Advanced Urothelial Carcinoma.

Histone deacetylase (HDAC) inhibition has sporadic clinical efficacy in urothelial carcinoma; the genomic basis for clinical response is not known. In two separate phase I clinical trials testing pharmacokinetic aspects of HDAC inhibitors in advanced solid tumors, we identified one patient with advanced urothelial carcinoma who had a complete response to belinostat, and one patient with advanced urothelial carcinoma who had a partial response to panobinostat. The archived tumors of the responders were genomically characterized in comparison to others with urothelial carcinoma on the trials. Urothelial carcinoma cell lines treated with panobinostat and belinostat were studied to elucidate the mechanisms of benefit. Notably, the urothelial carcinoma tumors that responded to HDAC inhibition had ARID1A mutations. ARID1A mutations were also noted in the tumors of three patients who had stable disease as their best response to HDAC inhibition. Corroborating the basis of sensitivity, transcriptional profiling of platinum-resistant ARID1A-mutated HT1197 cells treated with panobinostat reveals negative enrichment for both cyto-proliferative (MYC and E2F targets) and DNA repair gene sets, and positive enrichment for TP53 and inflammatory gene sets. Our study identifies ARID1A loss as a basis for clinical response to pan HDAC inhibition and offers avenues for potential rational therapeutic combinations with HDAC inhibitors in advanced urothelial carcinoma.

Publisher Correction: Combating subclonal evolution of resistant cancer phenotypes.

The originally published version of this Article contained an error in Figure 4. In panel a, grey boxes surrounding the subclones associated with patients #2 and #4 obscured adjacent portions of the heatmap. This error has now been corrected in both the PDF and HTML versions of the Article.

Studying Selenoprotein mRNA Translation Using RNA-Seq and Ribosome Profiling.

Deep sequencing of ribosome protected mRNA footprints, also called ribosome profiling or Ribo-Seq, is a relatively new methodology well suited to address questions regarding the mechanisms and efficiency of protein expression. Specifically, the ability of this technique to quantify ribosome abundance with codon resolution enables experiments aimed at studying many aspects of translation, including gene-specific translational efficiency, translation of regulatory upstream short open reading frames, sites of ribosome pausing, and most importantly for selenoproteins, the efficiency by which UGA codons are redefined to encode selenocysteine. Here, we describe a streamlined protocol that was developed in our lab to process mammalian tissue to produce the requisite matched ribosome profiling and RNA-Seq libraries for deep sequencing.

Combating subclonal evolution of resistant cancer phenotypes.

Metastatic breast cancer remains challenging to treat, and most patients ultimately progress on therapy. This acquired drug resistance is largely due to drug-refractory sub-populations (subclones) within heterogeneous tumors. Here, we track the genetic and phenotypic subclonal evolution of four breast cancers through years of treatment to better understand how breast cancers become drug-resistant. Recurrently appearing post-chemotherapy mutations are rare. However, bulk and single-cell RNA sequencing reveal acquisition of malignant phenotypes after treatment, including enhanced mesenchymal and growth factor signaling, which may promote drug resistance, and decreased antigen presentation and TNF-α signaling, which may enable immune system avoidance. Some of these phenotypes pre-exist in pre-treatment subclones that become dominant after chemotherapy, indicating selection for resistance phenotypes. Post-chemotherapy cancer cells are effectively treated with drugs targeting acquired phenotypes. These findings highlight cancer's ability to evolve phenotypically and suggest a phenotype-targeted treatment strategy that adapts to cancer as it evolves.

Leucine Differentially Regulates Gene-Specific Translation in Mouse Skeletal Muscle.

Background: Amino acids, especially leucine, are particularly effective in promoting protein synthesis. Leucine is known to increase the rate of protein synthesis in skeletal muscle through the mechanistic target of rapamycin complex 1-dependent, as well as -independent, signaling pathways. However, the overall translation program is poorly defined, and it is unknown how the activation of these pathways differentially controls the translation of specific mRNAs.Objective: Ribosome profiling and RNA sequencing were used to precisely define the translational program activated by an acute oral dose of leucine.Methods: Adult male C57BL/6 mice were deprived of food overnight before the delivery of an acute dose of l-leucine (9.4 mg) (n = 6) or vehicle (n = 5) and tissues collected 30 min later. Ribosome footprints and total RNA were isolated and subjected to deep sequencing. Changes in gene-specific mRNA abundance and ribosome occupancy were determined between the leucine-treated and control groups by aligning sequence reads to Reference Sequence database mRNAs and applying statistical features of the Bioconductor package edgeR.Results: Our data revealed mRNA features that confer translational control of skeletal muscle mRNAs in response to an acute dose of leucine. The subset of skeletal muscle mRNAs that are activated consists largely of terminal oligopyrimidine mRNAs (false discovery rate: <0.05), whereas those with reduced translation had 5' untranslated regions with increased length. Only the small nuclear RNAs, which are required for ribosome biogenesis, were significantly altered in RNA abundance. The inferred functional translational program activated by dietary leucine includes increased protein synthesis capacity and energy metabolism, upregulation of sarcomere-binding proteins, modulation of circadian rhythm, and suppression of select immune components.Conclusions: These results clarify the translation program acutely stimulated by leucine in mouse skeletal muscle and establish new methodologies for use in future studies of skeletal muscle disease or aging and further examination of downstream effects of leucine on gene expression.

A comprehensive multi-omics approach uncovers adaptations for growth and survival of Pseudomonas aeruginosa on n-alkanes.

BACKGROUND: Examination of complex biological systems has long been achieved through methodical investigation of the system's individual components. While informative, this strategy often leads to inappropriate conclusions about the system as a whole. With the advent of high-throughput "omic" technologies, however, researchers can now simultaneously analyze an entire system at the level of molecule (DNA, RNA, protein, metabolite) and process (transcription, translation, enzyme catalysis). This strategy reduces the likelihood of improper conclusions, provides a framework for elucidation of genotype-phenotype relationships, and brings finer resolution to comparative genomic experiments. Here, we apply a multi-omic approach to analyze the gene expression profiles of two closely related Pseudomonas aeruginosa strains grown in n-alkanes or glycerol. RESULTS: The environmental P. aeruginosa isolate ATCC 33988 consumed medium-length (C10-C16) n-alkanes more rapidly than the laboratory strain PAO1, despite high genome sequence identity (average nucleotide identity >99%). Our data shows that ATCC 33988 induces a characteristic set of genes at the transcriptional, translational and post-translational levels during growth on alkanes, many of which differ from those expressed by PAO1. Of particular interest was the lack of expression from the rhl operon of the quorum sensing (QS) system, resulting in no measurable rhamnolipid production by ATCC 33988. Further examination showed that ATCC 33988 lacked the entire lasI/lasR arm of the QS response. Instead of promoting expression of QS genes, ATCC 33988 up-regulates a small subset of its genome, including operons responsible for specific alkaline proteases and sphingosine metabolism. CONCLUSION: This work represents the first time results from RNA-seq, microarray, ribosome footprinting, proteomics, and small molecule LC-MS experiments have been integrated to compare gene expression in bacteria. Together, these data provide insights as to why strain ATCC 33988 is better adapted for growth and survival on n-alkanes.

Combining quantitative genetic footprinting and trait enrichment analysis to identify fitness determinants of a bacterial pathogen.

Strains of Extraintestinal Pathogenic Escherichia c oli (ExPEC) exhibit an array of virulence strategies and are a major cause of urinary tract infections, sepsis and meningitis. Efforts to understand ExPEC pathogenesis are challenged by the high degree of genetic and phenotypic variation that exists among isolates. Determining which virulence traits are widespread and which are strain-specific will greatly benefit the design of more effective therapies. Towards this goal, we utilized a quantitative genetic footprinting technique known as transposon insertion sequencing (Tn-seq) in conjunction with comparative pathogenomics to functionally dissect the genetic repertoire of a reference ExPEC isolate. Using Tn-seq and high-throughput zebrafish infection models, we tracked changes in the abundance of ExPEC variants within saturated transposon mutant libraries following selection within distinct host niches. Nine hundred and seventy bacterial genes (18% of the genome) were found to promote pathogen fitness in either a niche-dependent or independent manner. To identify genes with the highest therapeutic and diagnostic potential, a novel Trait Enrichment Analysis (TEA) algorithm was developed to ascertain the phylogenetic distribution of candidate genes. TEA revealed that a significant portion of the 970 genes identified by Tn-seq have homologues more often contained within the genomes of ExPEC and other known pathogens, which, as suggested by the first axiom of molecular Koch's postulates, is considered to be a key feature of true virulence determinants. Three of these Tn-seq-derived pathogen-associated genes--a transcriptional repressor, a putative metalloendopeptidase toxin and a hypothetical DNA binding protein--were deleted and shown to independently affect ExPEC fitness in zebrafish and mouse models of infection. Together, the approaches and observations reported herein provide a resource for future pathogenomics-based research and highlight the diversity of factors required by a single ExPEC isolate to survive within varying host environments.

Zinc-dependent regulation of the Adh1 antisense transcript in fission yeast.

In yeast, Adh1 (alcohol dehydrogenase 1) is an abundant zinc-binding protein that is required for the conversion of acetaldehyde to ethanol. Through transcriptome profiling of the Schizosaccharomyces pombe genome, we identified a natural antisense transcript at the adh1 locus that is induced in response to zinc limitation. This antisense transcript (adh1AS) shows a reciprocal expression pattern to that of the adh1 mRNA partner. In this study, we show that increased expression of the adh1AS transcript in zinc-limited cells is necessary for the repression of adh1 gene expression and that the increased level of the adh1AS transcript in zinc-limited cells is a result of two mechanisms. At the transcriptional level, the adh1AS transcript is expressed at a high level in zinc-limited cells. In addition to this transcriptional control, adh1AS transcripts preferentially accumulate in zinc-limited cells when the adh1AS transcript is expressed from a constitutive promoter. This secondary mechanism requires the simultaneous expression of adh1. Our studies reveal how multiple mechanisms can synergistically control the ratio of sense to antisense transcripts and highlight a novel mechanism by which adh1 gene expression can be controlled by cellular zinc availability.

Antibody detection of translocations in Ewing sarcoma.

The detection of chromosomal translocations has important implications in the diagnosis, prognosis and treatment of patients with cancer. Current approaches to translocation detection have significant shortcomings, including limited sensitivity and/or specificity, and difficulty in application to formalin-fixed paraffin-embedded (FFPE) clinical samples. We developed a new approach called antibody detection of translocations (ADOT) that avoids the shortcomings of current techniques. ADOT combines a transcriptional microarray-based approach with a novel antibody-based detection method. ADOT allows for the accurate and sensitive identification of translocations and provides exon-level information about the fusion transcript. ADOT can detect translocations in poor-quality unprocessed total ribonucleic acid (RNA). Furthermore, the technique is readily generalizable to detect any potential fusion transcript, including previously undescribed fusions. We demonstrate the feasibility of ADOT by examples in which both known and unknown Ewing sarcoma translocations are identified from cell lines, tumour xenografts and FFPE primary tumours. These results demonstrate that ADOT may be an effective approach for translocation analysis in clinical specimens with significant RNA degradation and may offer a novel diagnostic tool for translocation-based cancers.

Polymorphisms in the NPY2R gene show significant associations with BMI that are additive to FTO, MC4R, and NPFFR2 gene effects.

Neuropeptide Y (NPY) is an appetite hormone that acts centrally to control feeding behavior. The 5' and exon 2 regions of NPY2R, one of five NPY receptor genes, have been weakly and inconsistently implicated with obesity. With the ATG start site of the gene at the beginning of exon 2, single-nucleotide polymorphisms (SNPs) across intron 1 may show stronger associations with obesity than expected. Two 5' SNPs, three intron 1 SNPs, and one synonymous exon 2 SNP were genotyped on 2,985 white Utah subjects. Previously associated FTO, NPY, NPY1R, MC4R, PPARGC1A, OR7D4, and four NPFFR2 SNPs were also genotyped and related to BMI. One NPY2R 5' SNP (rs12649641, P = 0.008), an exon 2 SNP (rs2880415, P = 0.009), and an intron 1 SNP (rs17376826, P = 7 × 10(-6)) were each significantly associated with BMI. All three SNPs, plus FTO (rs9939609, P = 1.5 × 10(-6)) and two NPFFR2 SNPs (rs4129733, P = 3.7 × 10(-13) and rs11940196, 4.2 × 10(-10)) remained significant in a multiple regression additive model. Diplotypes using the estimated haplotypes of NPY2R, NPFFR2, and MC4R were significantly associated with BMI (P = 1.0 × 10(-10), 3.2 × 10(-8), and 1.1 × 10(-4), respectively). Haplotypes of NPY2R, NPFFR2, and MC4R, plus the FTO SNP, explained 9.6% of the BMI variance. SNP effect sizes per allele for the four genes ranged from 0.8 to 3.5 kg/m(2). We conclude that haplotypes containing the rs17376826 SNP in intron 1 of NPY2R have strong associations with BMI, some NPFFR2 haplotypes are strongly protective against or increase risk of obesity, and both NPY2R and NPFFR2 play important roles in obesity predisposition independent of FTO and MC4R.

Using VAAST to identify an X-linked disorder resulting in lethality in male infants due to N-terminal acetyltransferase deficiency.

We have identified two families with a previously undescribed lethal X-linked disorder of infancy; the disorder comprises a distinct combination of an aged appearance, craniofacial anomalies, hypotonia, global developmental delays, cryptorchidism, and cardiac arrhythmias. Using X chromosome exon sequencing and a recently developed probabilistic algorithm aimed at discovering disease-causing variants, we identified in one family a c.109T>C (p.Ser37Pro) variant in NAA10, a gene encoding the catalytic subunit of the major human N-terminal acetyltransferase (NAT). A parallel effort on a second unrelated family converged on the same variant. The absence of this variant in controls, the amino acid conservation of this region of the protein, the predicted disruptive change, and the co-occurrence in two unrelated families with the same rare disorder suggest that this is the pathogenic mutation. We confirmed this by demonstrating a significantly impaired biochemical activity of the mutant hNaa10p, and from this we conclude that a reduction in acetylation by hNaa10p causes this disease. Here we provide evidence of a human genetic disorder resulting from direct impairment of N-terminal acetylation, one of the most common protein modifications in humans.

Next generation tools for genomic data generation, distribution, and visualization.

BACKGROUND: With the rapidly falling cost and availability of high throughput sequencing and microarray technologies, the bottleneck for effectively using genomic analysis in the laboratory and clinic is shifting to one of effectively managing, analyzing, and sharing genomic data. RESULTS: Here we present three open-source, platform independent, software tools for generating, analyzing, distributing, and visualizing genomic data. These include a next generation sequencing/microarray LIMS and analysis project center (GNomEx); an application for annotating and programmatically distributing genomic data using the community vetted DAS/2 data exchange protocol (GenoPub); and a standalone Java Swing application (GWrap) that makes cutting edge command line analysis tools available to those who prefer graphical user interfaces. Both GNomEx and GenoPub use the rich client Flex/Flash web browser interface to interact with Java classes and a relational database on a remote server. Both employ a public-private user-group security model enabling controlled distribution of patient and unpublished data alongside public resources. As such, they function as genomic data repositories that can be accessed manually or programmatically through DAS/2-enabled client applications such as the Integrated Genome Browser. CONCLUSIONS: These tools have gained wide use in our core facilities, research laboratories and clinics and are freely available for non-profit use. See http://sourceforge.net/projects/gnomex/, http://sourceforge.net/projects/genoviz/, and http://sourceforge.net/projects/useq.

Dynamic transcriptome of Schizosaccharomyces pombe shown by RNA-DNA hybrid mapping.

We have determined the high-resolution strand-specific transcriptome of the fission yeast S. pombe under multiple growth conditions using a novel RNA-DNA hybridization mapping (HybMap) technique. HybMap uses an antibody against an RNA-DNA hybrid to detect RNA molecules hybridized to a high-density DNA oligonucleotide tiling microarray. HybMap showed exceptional dynamic range and reproducibility, and allowed us to identify strand-specific coding, noncoding and structural RNAs, as well as previously unknown RNAs conserved in distant yeast species. Notably, we found that virtually the entire euchromatic genome (including intergenics) is transcribed, with heterochromatin dampening intergenic transcription. We identified features including large numbers of condition-specific noncoding RNAs, extensive antisense transcription, new properties of antisense transcripts and induced divergent transcription. Furthermore, our HybMap data informed the efficiency and locations of RNA splicing genome-wide. Finally, we observed strand-specific transcription islands around tRNAs at heterochromatin boundaries inside centromeres. Here, we discuss these new features in terms of organism fitness and transcriptome evolution.

Gains in sensitivity with a device that mixes microarray hybridization solution in a 25-microm-thick chamber.

A microarray hybridization system that allows mixing in volumes comparable to those used by glass coverslips is presented. This system is composed of a disposable flexible lid that binds to 1 in. x 3 in. glass slides via an adhesive gasket, forming a uniform 25-microm-thick hybridization chamber. This chamber rests on a base unit for temperature control. The lid contains two air-driven bladders that continuously mix the hybridization fluid. Mixing enhances sensitivity from a typical microarray experiment 2-3-fold. Mixing is particularly effective at high spotted probe and low labeled target concentrations and overcoming local target depletion that occurs when homologous probes are spotted in close proximity. Mixing appears to be compatible with most hybridization conditions; however, mix versus no-mix control experiments should be performed. Also covered are a number of microfluidic issues related to manufacturing, filling, mixing, and packaging.