categoryCompare, an analytical tool based on feature annotations

Polycomb deficiency drives a FOXP2-high aggressive state targetable by epigenetic inhibitors

Inhibitors of the Polycomb Repressive Complex 2 (PRC2) histone methyltransferase EZH2 are approved for certain cancers, but realizing their wider utility relies upon understanding PRC2 biology in each cancer system. Using a genetic model to delete Ezh2 in KRAS-driven lung adenocarcinomas, we observed that Ezh2 haplo-insufficient tumors were less lethal and lower grade than Ezh2 fully-insufficient tumors, which were poorly differentiated and metastatic. Using three-dimensional cultures and in vivo experiments, we determined that EZH2-deficient tumors were vulnerable to H3K27 demethylase or BET inhibitors. PRC2 loss/inhibition led to de-repression of FOXP2, a transcription factor that promotes migration and stemness, and FOXP2 could be suppressed by BET inhibition. Poorly differentiated human lung cancers were enriched for an H3K27me3-low state, representing a subtype that may benefit from BET inhibition as a single therapy or combined with additional EZH2 inhibition. These data highlight diverse roles of PRC2 in KRAS-driven lung adenocarcinomas, and demonstrate the utility of three-dimensional cultures for exploring epigenetic drug sensitivities for cancer.

Shoc2 controls ERK1/2-driven neural crest development by balancing components of the extracellular matrix

The extracellular signal-regulated kinase (ERK1/2) pathway is essential in embryonic development. The scaffold protein Shoc2 is a critical modulator of ERK1/2 signals, and mutations in the shoc2 gene lead to the human developmental disease known as Noonan-like syndrome with loose anagen hair (NSLH). The loss of Shoc2 and the shoc2 NSLH-causing mutations affect the tissues of neural crest (NC) origin. In this study, we utilized the zebrafish model to dissect the role of Shoc2-ERK1/2 signals in the development of NC. These studies established that the loss of Shoc2 significantly altered the expression of transcription factors regulating the specification and differentiation of NC cells. Using comparative transcriptome analysis of NC-derived cells from shoc2 CRISPR/Cas9 mutant larvae, we found that Shoc2-mediated signals regulate gene programs at several levels, including expression of genes coding for the proteins of extracellular matrix (ECM) and ECM regulators. Together, our results demonstrate that Shoc2 is an essential regulator of NC development. This study also indicates that disbalance in the turnover of the ECM may lead to the abnormalities found in NSLH patients.

Deletion of pancreas-specific miR-216a reduces beta-cell mass and inhibits pancreatic cancer progression in mice

miRNAs have crucial functions in many biological processes and are candidate biomarkers of disease. Here, we show that miR-216a is a conserved, pancreas-specific miRNA with important roles in pancreatic islet and acinar cells. Deletion of miR-216a in mice leads to a reduction in islet size, β-cell mass, and insulin levels. Single-cell RNA sequencing reveals a subpopulation of β-cells with upregulated acinar cell markers under a high-fat diet. miR-216a is induced by TGF-β signaling, and inhibition of miR-216a increases apoptosis and decreases cell proliferation in pancreatic cells. Deletion of miR-216a in the pancreatic cancer-prone mouse line KrasG12D;Ptf1aCreER reduces the propensity of pancreatic cancer precursor lesions. Notably, circulating miR-216a levels are elevated in both mice and humans with pancreatic cancer. Collectively, our study gives insights into how β-cell mass and acinar cell growth are modulated by a pancreas-specific miRNA and also suggests miR-216a as a potential biomarker for diagnosis of pancreatic diseases.

Untargeted Lipidomics of Non-Small Cell Lung Carcinoma Demonstrates Differentially Abundant Lipid Classes in Cancer vs. Non-Cancer Tissue

Lung cancer remains the leading cause of cancer death worldwide and non-small cell lung carcinoma (NSCLC) represents 85% of newly diagnosed lung cancers. In this study, we utilized our untargeted assignment tool Small Molecule Isotope Resolved Formula Enumerator (SMIRFE) and ultra-high-resolution Fourier transform mass spectrometry to examine lipid profile differences between paired cancerous and non-cancerous lung tissue samples from 86 patients with suspected stage I or IIA primary NSCLC. Correlation and co-occurrence analysis revealed significant lipid profile differences between cancer and non-cancer samples. Further analysis of machine-learned lipid categories for the differentially abundant molecular formulas identified a high abundance sterol, high abundance and high m/z sphingolipid, and low abundance glycerophospholipid metabolic phenotype across the NSCLC samples. At the class level, higher abundances of sterol esters and lower abundances of cardiolipins were observed suggesting altered stearoyl-CoA desaturase 1 (SCD1) or acetyl-CoA acetyltransferase (ACAT1) activity and altered human cardiolipin synthase 1 or lysocardiolipin acyltransferase activity respectively, the latter of which is known to confer apoptotic resistance. The presence of a shared metabolic phenotype across a variety of genetically distinct NSCLC subtypes suggests that this phenotype is necessary for NSCLC development and may result from multiple distinct genetic lesions. Thus, targeting the shared affected pathways may be beneficial for a variety of genetically distinct NSCLC subtypes.

Differential gene expression analysis of palbociclib-resistant TNBC via RNA-seq

Purpose

The management of triple-negative breast cancer (TNBC) remains a significant clinical challenge due to the lack of effective targeted therapies. Inhibitors of the cyclin-dependent kinases 4 and 6 (CDK4/6) are emerging as promising therapeutic agents against TNBC; however, cells can rapidly acquire resistance through multiple mechanisms that are yet to be identified. Therefore, determining the mechanisms underlying resistance to CDK4/6 inhibition is crucial to develop combination therapies that can extend the efficacy of the CDK4/6 inhibitors or delay resistance. This study aims to identify differentially expressed genes (DEG) associated with acquired resistance to palbociclib in ER− breast cancer cells.

Methods

We performed next-generation transcriptomic sequencing (RNA-seq) and pathway analysis in ER− MDA-MB-231 palbociclib-sensitive (231/pS) and palbociclib-resistant (231/pR) cells.

Results

We identified 2247 up-regulated and 1427 down-regulated transcripts in 231/pR compared to 231/pS cells. DEGs were subjected to functional analysis using Gene Ontology (GO) and the KEGG database which identified many transduction pathways associated with breast cancer, including the PI3K/AKT, PTEN and mTOR pathways. Additionally, Ingenuity Pathway Analysis (IPA) revealed that resistance to palbociclib is closely associated with altered cholesterol and fatty acid biosynthesis suggesting that resistance to palbociclib may be dependent on lipid metabolic reprograming.

Conclusion

This study provides evidence that lipid metabolism is altered in TNBC with acquired resistance to palbociclib. Further studies are needed to determine if the observed lipid metabolic rewiring can be exploited to overcome therapy resistance in TNBC.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-021-06127-5.

Loss of CLN3, the gene mutated in juvenile neuronal ceroid lipofuscinosis, leads to metabolic impairment and autophagy induction in retinal pigment epithelium

Juvenile neuronal ceroid lipofuscinosis (JNCL, aka. juvenile Batten disease or CLN3 disease) is a lysosomal storage disease characterized by progressive blindness, seizures, cognitive and motor failures, and premature death. JNCL is caused by mutations in the Ceroid Lipofuscinosis, Neuronal 3 (CLN3) gene, whose function is unclear. Although traditionally considered a neurodegenerative disease, CLN3 disease displays eye-specific effects: Vision loss not only is often one of the earliest symptoms of JNCL, but also has been reported in non-syndromic CLN3 disease. Here we described the roles of CLN3 protein in maintaining healthy retinal pigment epithelium (RPE) and normal vision. Using electroretinogram, fundoscopy and microscopy, we showed impaired visual function, retinal autofluorescent lesions, and RPE disintegration and metaplasia/hyperplasia in a Cln3 ~ 1 kb-deletion mouse model [1] on C57BL/6J background. Utilizing a combination of biochemical analyses, RNA-Seq, Seahorse XF bioenergetic analysis, and Stable Isotope Resolved Metabolomics (SIRM), we further demonstrated that loss of CLN3 increased autophagic flux, suppressed mTORC1 and Akt activities, enhanced AMPK activity, and up-regulated gene expression of the autophagy-lysosomal system in RPE-1 cells, suggesting autophagy induction. This CLN3 deficiency induced autophagy induction coincided with decreased mitochondrial oxygen consumption, glycolysis, the tricarboxylic acid (TCA) cycle, and ATP production. We also reported for the first time that loss of CLN3 led to glycogen accumulation despite of impaired glycogen synthesis. Our comprehensive analyses shed light on how loss of CLN3 affect autophagy and metabolism. This work suggests possible links among metabolic impairment, autophagy induction and lysosomal storage, as well as between RPE atrophy/degeneration and vision loss in JNCL.

Transcriptomic Profiling Identifies Differentially Expressed Genes in Palbociclib-Resistant ER+ MCF7 Breast Cancer Cells

Acquired resistance to cyclin-dependent kinases 4 and 6 (CDK4/6) inhibition in estrogen receptor-positive (ER+) breast cancer remains a significant clinical challenge. Efforts to uncover the mechanisms underlying resistance are needed to establish clinically actionable targets effective against resistant tumors. In this study, we sought to identify differentially expressed genes (DEGs) associated with acquired resistance to palbociclib in ER+ breast cancer. We performed next-generation transcriptomic RNA sequencing (RNA-seq) and pathway analysis in ER+ MCF7 palbociclib-sensitive (MCF7/pS) and MCF7 palbociclib-resistant (MCF7/pR) cells. We identified 2183 up-regulated and 1548 down-regulated transcripts in MCF7/pR compared to MCF7/pS cells. Functional analysis of the DEGs using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database identified several pathways associated with breast cancer, including ‘cell cycle’, ‘DNA replication’, ‘DNA repair’ and ‘autophagy’. Additionally, Ingenuity Pathway Analysis (IPA) revealed that resistance to palbociclib is closely associated with deregulation of several key canonical and metabolic pathways. Further studies are needed to determine the utility of these DEGs and pathways as therapeutics targets against ER+ palbociclib-resistant breast cancer.

Dataset for dose and time-dependent transcriptional response to ionizing radiation exposure

Exposure to ionizing radiation associated with highly energetic and charged heavy particles is an inherent risk astronauts face in long duration space missions. We have previously considered the transcriptional effects that three levels of radiation (0.3 Gy, 1.5 Gy, and 3.0 Gy) have at an immediate time point (1 hr) post-exposure [1]. Our analysis of these results suggest effects on transcript levels that could be modulated at lower radiation doses [2]. In addition, a time dependent effect is likely to be present. Therefore, in order to develop a lab-on-a-chip approach for detection of radiation exposure in terms of both radiation level and time since exposure, we developed a time- and dose-course study to determine appropriate sensitive and specific transcript biomarkers that are detectable in blood samples. The data described herein was developed from a study measuring exposure to 0.15 Gy, 0.30 Gy, and 1.5 Gy of radiation at 1 hr, 2 hr, and 6 hr post-exposure using Affymetrix® GeneChip® PrimeView™ microarrays. This report includes raw gene expression data files from the resulting microarray experiments representing typical radiation exposure levels an astronaut may experience as part of a long duration space mission. The data described here is available in NCBI's Gene Expression Omnibus (GEO), accession GSE63952.

Comparative gene expression analysis in melanocytes driven by tumor cell-derived exosomes

Abundant with organelle-like membranous structures, the tumor microenvironment is composed of cancer cells that secrete exosomes. Studies have shown that these secreted exosomes transport RNA and active molecules to other cells to reshape the tumor microenvironment and promote tumor growth. In fact, we found that exosomes derived from melanoma cells drive pre-malignant transition in primary melanocytes. However, there is little available in the scientific literature on how exosomes modulate melanocytes in the microenvironment to optimize conditions for tumor progression and metastasis. We therefore focused this current study on identifying these conditions genetically. Through RNA sequencing, we analyzed gene expression levels of melanocytes driven by exosomes derived from melanoma and lung cancer cells compared with those without exosome controls. Significant differences were found in gene expression patterns of melanocytes driven by exosomes derived from melanoma and lung cancer cells. In the melanocytes responding to exosomes derived from melanoma cells, genes of lipopolysaccharide and regulation of leukocyte chemotaxis were predominant. In the melanocytes responding to exosomes derived from lung cancer cells, genes of DNA replication and mitotic nuclear division played an important role. These results provide further mechanistic understanding of tumor progression promoted by tumor-derived exosomes. This will also help identify potential therapeutic targets for melanoma progression.

Advances in gene ontology utilization improve statistical power of annotation enrichment

Gene-annotation enrichment is a common method for utilizing ontology-based annotations in gene and gene-product centric knowledgebases. Effective utilization of these annotations requires inferring semantic linkages by tracing paths through edges in the ontological graph, referred to as relations. However, some relations are semantically problematic with respect to scope, necessitating their omission or modification lest erroneous term mappings occur. To address these issues, we created the Gene Ontology Categorization Suite, or GOcats-a novel tool that organizes the Gene Ontology into subgraphs representing user-defined concepts, while ensuring that all appropriate relations are congruent with respect to scoping semantics. Here, we demonstrate the improvements in annotation enrichment by re-interpreting edges that would otherwise be omitted by traditional ancestor path-tracing methods. Specifically, we show that GOcats' unique handling of relations improves enrichment over conventional methods in the analysis of two different gene-expression datasets: a breast cancer microarray dataset and several horse cartilage development RNAseq datasets. With the breast cancer microarray dataset, we observed significant improvement (one-sided binomial test p-value = 1.86E-25) in 182 of 217 significantly enriched GO terms identified from the conventional path traversal method when GOcats' path traversal was used. We also found new significantly enriched terms using GOcats, whose biological relevancy has been experimentally demonstrated elsewhere. Likewise, on the horse RNAseq datasets, we observed a significant improvement in GO term enrichment when using GOcat's path traversal: one-sided binomial test p-values range from 1.32E-03 to 2.58E-44.

HNRNPA2/B1 is upregulated in endocrine-resistant LCC9 breast cancer cells and alters the miRNA transcriptome when overexpressed in MCF-7 cells

MicroRNAs are dysregulated in breast cancer. Heterogeneous Nuclear Ribonucleoprotein A2/B1 (HNRNPA2/B1) is a reader of the N(6)-methyladenosine (m6A) mark in primary-miRNAs (pri-miRNAs) and promotes DROSHA processing to precursor-miRNAs (pre-miRNAs). We examined the expression of writers, readers, and erasers of m6A and report that HNRNPA2/B1 expression is higher in tamoxifen-resistant LCC9 breast cancer cells as compared to parental, tamoxifen-sensitive MCF-7 cells. To examine how increased expression of HNRNPA2/B1 affects miRNA expression, HNRNPA2/B1 was transiently overexpressed (~5.4-fold) in MCF-7 cells for whole genome miRNA profiling (miRNA-seq). 148 and 88 miRNAs were up- and down-regulated, respectively, 48 h after transfection and 177 and 172 up- and down-regulated, respectively, 72 h after transfection. MetaCore Enrichment analysis identified progesterone receptor action and transforming growth factor β (TGFβ) signaling via miRNA in breast cancer as pathways downstream of the upregulated miRNAs and TGFβ signaling via SMADs and Notch signaling as pathways of the downregulated miRNAs. GO biological processes for mRNA targets of HNRNPA2/B1-regulated miRNAs included response to estradiol and cell-substrate adhesion. qPCR confirmed HNRNPA2B1 downregulation of miR-29a-3p, miR-29b-3p, and miR-222 and upregulation of miR-1266-5p, miR-1268a, miR-671-3p. Transient overexpression of HNRNPA2/B1 reduced MCF-7 sensitivity to 4-hydroxytamoxifen and fulvestrant, suggesting a role for HNRNPA2/B1 in endocrine-resistance.

Detection of Differentially Expressed Cleavage Site Intervals Within 3' Untranslated Regions Using CSI-UTR Reveals Regulated Interaction Motifs

The length of untranslated regions at the 3' end of transcripts (3'UTRs) is regulated by alternate polyadenylation (APA). 3'UTRs contain regions that harbor binding motifs for regulatory molecules. However, the mechanisms that coordinate the 3'UTR length of specific groups of transcripts are not well-understood. We therefore developed a method, CSI-UTR, that models 3'UTR structure as tandem segments between functional alternative-polyadenylation sites (termed cleavage site intervals-CSIs). This approach facilitated (1) profiling of 3'UTR isoform expression changes and (2) statistical enrichment of putative regulatory motifs. CSI-UTR analysis is UTR-annotation independent and can interrogate legacy data generated from standard RNA-Seq libraries. CSI-UTR identified a set of CSIs in human and rodent transcriptomes. Analysis of RNA-Seq datasets from neural tissue identified differential expression events within 3'UTRs not detected by standard gene-based differential expression analyses. Further, in many instances 3'UTR and CDS from the same gene were regulated differently. This modulation of motifs for RNA-interacting molecules with potential condition-dependent and tissue-specific RNA binding partners near the polyA signal and CSI junction may play a mechanistic role in the specificity of alternative polyadenylation. Source code, CSI BED files and example datasets are available at: https://github.com/UofLBioinformatics/CSI-UTR.

Coupling of PARP1-mediated chromatin structural changes to transcriptional RNA polymerase II elongation and cotranscriptional splicing

Background

Recently, we showed that PARP1 is involved in cotranscriptional splicing, possibly by bridging chromatin to RNA and recruiting splicing factors. It also can influence alternative splicing decisions through the regulation of RNAPII elongation. In this study, we investigated the effect of PARP1-mediated chromatin changes on RNAPII movement, during transcription and alternative splicing.

Results

We show that RNAPII pauses at PARP1–chromatin structures within the gene body. Knockdown of PARP1 abolishes this RNAPII pausing, suggesting that PARP1 may regulate RNAPII elongation. Additionally, PARP1 alters nucleosome deposition and histone post-translational modifications at specific exon–intron boundaries, thereby affecting RNAPII movement. Lastly, genome-wide analyses confirmed that PARP1 influences changes in RNAPII elongation by either reducing or increasing the rate of RNAPII elongation depending on the chromatin context.

Conclusions

These studies suggest a context-specific effect of PARP1–chromatin binding on RNA polymerase movement and provide a platform to delineate PARP1’s role in RNA biogenesis and processing.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0261-1) contains supplementary material, which is available to authorized users.

Transcriptomic response of breast cancer cells to anacardic acid

Anacardic acid (AnAc), a potential dietary agent for preventing and treating breast cancer, inhibited the proliferation of estrogen receptor α (ERα) positive MCF-7 and MDA-MB-231 triple negative breast cancer cells. To characterize potential regulators of AnAc action, MCF-7 and MDA-MB-231 cells were treated for 6 h with purified AnAc 24:1n5 congener followed by next generation transcriptomic sequencing (RNA-seq) and network analysis. We reported that AnAc-differentially regulated miRNA transcriptomes in each cell line and now identify AnAc-regulated changes in mRNA and lncRNA transcript expression. In MCF-7 cells, 80 AnAc-responsive genes were identified, including lncRNA MIR22HG. More AnAc-responsive genes (886) were identified in MDA-MB-231 cells. Only six genes were commonly altered by AnAc in both cell lines: SCD, INSIG1, and TGM2 were decreased and PDK4, GPR176, and ZBT20 were increased. Modeling of AnAc-induced gene changes suggests that AnAc inhibits monounsaturated fatty acid biosynthesis in both cell lines and increases endoplasmic reticulum stress in MDA-MB-231 cells. Since modeling of downregulated genes implicated NFκB in MCF-7, we confirmed that AnAc inhibited TNFα-induced NFκB reporter activity in MCF-7 cells. These data identify new targets and pathways that may account for AnAc’s anti-proliferative and pro-apoptotic activity.

Pas de deux: An Intricate Dance of Anther Smut and Its Host

The successful interaction between pathogen/parasite and host requires a delicate balance between fitness of the former and survival of the latter. To optimize fitness a parasite/pathogen must effectively create an environment conducive to reproductive success, while simultaneously avoiding or minimizing detrimental host defense response. The association between Microbotryum lychnidis-dioicae and its host Silene latifolia serves as an excellent model to examine such interactions. This fungus is part of a species complex that infects species of the Caryophyllaceae, replacing pollen with the fungal spores. In the current study, transcriptome analyses of the fungus and its host were conducted during discrete stages of bud development so as to identify changes in fungal gene expression that lead to spore development and to identify changes associated with infection in the host plant. In contrast to early biotrophic phase stages of infection for the fungus, the latter stages involve tissue necrosis and in the case of infected female flowers, further changes in the developmental program in which the ovary aborts and a pseudoanther is produced. Transcriptome analysis via Illumina RNA sequencing revealed enrichment of fungal genes encoding small secreted proteins, with hallmarks of effectors and genes found to be relatively unique to the Microbotryum species complex. Host gene expression analyses also identified interesting sets of genes up-regulated, including those involving stress response, host defense response, and several agamous-like MADS-box genes (AGL61 and AGL80), predicted to interact and be involved in male gametophyte development.

Regulation of miR-29b-1/a transcription and identification of target mRNAs in CHO-K1 cells

miR-29b and miR-29a transcript levels were reported to increase in exponentially growing CHO-K1 cells. Here, we examine the regulation of miR-29b-1/a in CHO-K1 cells. We observed that 4-hydroxytamoxifen (4-OHT) increased pri-miR-29b-1 and pri-miR-29a transcription in CHO-K1 cells by activating endogenous estrogen receptor α (ERα). DICER, an established, bona fide target of miR-29b-1/a, was shown to be regulated by 4-OHT in CHO-K1 cells. We showed that miR-29b-1 and miR-29a serve a repressive role in cell proliferation, migration, invasion, and colony formation in CHO-K1 cells. To identify other targets of miR-29b-1 and miR-29a, RNA sequencing was performed by transfecting cells with anti-miR-29a, which inhibits both miR-29a and miR-29b-1, pre-miR-29b-1, and/or pre-miR-29a. In silico network analysis in MetaCore™ identified common and unique putative gene targets of miR-29b-1 and miR-29a. Pathway analysis of identified putative miR-29 targets were related to cell adhesion, cytoskeletal remodeling, and development. Further inquiry revealed regulation of pathways mediating responses to growth factor stimulus and cell cycle regulation.

Aberrant coordination geometries discovered in the most abundant metalloproteins

Metalloproteins bind and utilize metal ions for a variety of biological purposes. Due to the ubiquity of metalloprotein involvement throughout these processes across all domains of life, how proteins coordinate metal ions for different biochemical functions is of great relevance to understanding the implementation of these biological processes. Toward these ends, we have improved our methodology for structurally and functionally characterizing metal binding sites in metalloproteins. Our new ligand detection method is statistically much more robust, producing estimated false positive and false negative rates of ∼0.11% and ∼1.2%, respectively. Additional improvements expand both the range of metal ions and their coordination number that can be effectively analyzed. Also, the inclusion of additional quality control filters has significantly improved structure-function Spearman correlations as demonstrated by rho values greater than 0.90 for several metal coordination analyses and even one rho value above 0.95. Also, improvements in bond-length distributions have revealed bond-length modes specific to chemical functional groups involved in multidentation. Using these improved methods, we analyzed all single metal ion binding sites with Zn, Mg, Ca, Fe, and Na ions in the wwPDB, producing statistically rigorous results supporting the existence of both a significant number of unexpected compressed angles and subsequent aberrant metal ion coordination geometries (CGs) within structurally known metalloproteins. By recognizing these aberrant CGs in our clustering analyses, high correlations are achieved between structural and functional descriptions of metal ion coordination. Moreover, distinct biochemical functions are associated with aberrant CGs versus nonaberrant CGs. Proteins 2017; 85:885-907. © 2016 Wiley Periodicals, Inc.

Nuclear respiratory factor-1 and bioenergetics in tamoxifen-resistant breast cancer cells

Acquired tamoxifen (TAM) resistance is a significant clinical problem in treating patients with estrogen receptor α (ERα)+ breast cancer. We reported that ERα increases nuclear respiratory factor-1 (NRF-1), which regulates nuclear-encoded mitochondrial gene transcription, in MCF-7 breast cancer cells and NRF-1 knockdown stimulates apoptosis. Whether NRF-1 and target gene expression is altered in endocrine resistant breast cancer cells is unknown. We measured NRF-1and metabolic features in a cell model of progressive TAM-resistance. NRF-1 and its target mitochondrial transcription factor A (TFAM) were higher in TAM-resistant LCC2 and LCC9 cells than TAM-sensitive MCF-7 cells. Using extracellular flux assays we observed that LCC1, LCC2, and LCC9 cells showed similar oxygen consumption rate (OCR), but lower mitochondrial reserve capacity which was correlated with lower Succinate Dehydrogenase Complex, Subunit B in LCC1 and LCC2 cells. Complex III activity was lower in LCC9 than MCF-7 cells. LCC1, LCC2, and LCC9 cells had higher basal extracellular acidification (ECAR), indicating higher aerobic glycolysis, relative to MCF-7 cells. Mitochondrial bioenergetic responses to estradiol and 4-hydroxytamoxifen were reduced in the endocrine-resistant cells compared to MCF-7 cells. These results suggest the acquisition of altered metabolic phenotypes in response to long term antiestrogen treatment may increase vulnerability to metabolic stress.

Data set for transcriptional response to depletion of the Shoc2 scaffolding protein

The Suppressor of Clear, Caenorhabditis elegans Homolog (SHOC2) is a scaffold protein that positively modulates activity of the RAS/ERK1/2 MAP kinase signaling cascade. We set out to understand the ERK1/2 pathway transcriptional response transduced through the SHOC2 scaffolding module. This data article describes raw gene expression within triplicates of kidney fibroblast-like Cos1 cell line expressing non-targeting shRNA (Cos-NT) and triplicates of Cos1 cells depleted of SHOC2 using shRNA (Cos-LV1) upon activation of ERK1/2 pathway by the Epidermal Growth Factor Receptor (EGFR). The data referred here is available in NCBI׳s Gene Expression Omnibus (GEO), accession GEO: GSE67063 as well as NCBI׳s Sequence Read Archive (SRA), accession SRA: SRP056324. A complete analysis of the results can be found in “Shoc2-tranduced ERK1/2 motility signals – Novel insights from functional genomics”(Jeoung et al., 2016) [1].

Involvement of PARP1 in the regulation of alternative splicing

Specialized chromatin structures such as nucleosomes with specific histone modifications decorate exons in eukaryotic genomes, suggesting a functional connection between chromatin organization and the regulation of pre-mRNA splicing. Through profiling the functional location of Poly (ADP) ribose polymerase, we observed that it is associated with the nucleosomes at exon/intron boundaries of specific genes, suggestive of a role for this enzyme in alternative splicing. Poly (ADP) ribose polymerase has previously been implicated in the PARylation of splicing factors as well as regulation of the histone modification H3K4me3, a mark critical for co-transcriptional splicing. In light of these studies, we hypothesized that interaction of the chromatin-modifying factor, Poly (ADP) ribose polymerase with nucleosomal structures at exon–intron boundaries, might regulate pre-mRNA splicing. Using genome-wide approaches validated by gene-specific assays, we show that depletion of PARP1 or inhibition of its PARylation activity results in changes in alternative splicing of a specific subset of genes. Furthermore, we observed that PARP1 bound to RNA, splicing factors and chromatin, suggesting that Poly (ADP) ribose polymerase serves as a gene regulatory hub to facilitate co-transcriptional splicing. These studies add another function to the multi-functional protein, Poly (ADP) ribose polymerase, and provide a platform for further investigation of this protein’s function in organizing chromatin during gene regulatory processes.

Transcriptional profile of immediate response to ionizing radiation exposure

Astronauts participating in long duration space missions are likely to be exposed to ionizing radiation associated with highly energetic and charged heavy particles. Previously proposed gene biomarkers for radiation exposure include phosphorylated H2A Histone Family, Member X (γH2AX), Tumor Protein 53 (TP53), and Cyclin-Dependent Kinase Inhibitor 1A (CDKN1A). However, transcripts of these genes may not be the most suitable biomarkers for radiation exposure due to a lack of sensitivity or specificity. As part of a larger effort to develop lab-on-a-chip methods for detecting radiation exposure events using blood samples, we designed a dose–course microarray study in order to determine coding and non-coding RNA transcripts undergoing differential expression immediately following radiation exposure. The main goal was to elicit a small set of sensitive and specific radiation exposure biomarkers at low, medium, and high levels of ionizing radiation exposure. Four separate levels of radiation were considered: 0 Gray (Gy) control; 0.3 Gy; 1.5 Gy; and 3.0 Gy with four replicates at each radiation level. This report includes raw gene expression data files from the resulting microarray experiments from all three radiation levels ranging from a lower, typical exposure than an astronaut might see (0.3 Gy) to high, potentially lethal, levels of radiation (3.0 Gy). The data described here is available in NCBI's Gene Expression Omnibus (GEO), accession GSE64375.

Genome-Wide Profiling of PARP1 Reveals an Interplay with Gene Regulatory Regions and DNA Methylation

Poly (ADP-ribose) polymerase-1 (PARP1) is a nuclear enzyme involved in DNA repair, chromatin remodeling and gene expression. PARP1 interactions with chromatin architectural multi-protein complexes (i.e. nucleosomes) alter chromatin structure resulting in changes in gene expression. Chromatin structure impacts gene regulatory processes including transcription, splicing, DNA repair, replication and recombination. It is important to delineate whether PARP1 randomly associates with nucleosomes or is present at specific nucleosome regions throughout the cell genome. We performed genome-wide association studies in breast cancer cell lines to address these questions. Our studies show that PARP1 associates with epigenetic regulatory elements genome-wide, such as active histone marks, CTCF and DNase hypersensitive sites. Additionally, the binding of PARP1 to chromatin genome-wide is mutually exclusive with DNA methylation pattern suggesting a functional interplay between PARP1 and DNA methylation. Indeed, inhibition of PARylation results in genome-wide changes in DNA methylation patterns. Our results suggest that PARP1 controls the fidelity of gene transcription and marks actively transcribed gene regions by selectively binding to transcriptionally active chromatin. These studies provide a platform for developing our understanding of PARP1’s role in gene regulation.

Inorganic Arsenic-induced cellular transformation is coupled with genome wide changes in chromatin structure, transcriptome and splicing patterns

Background

Arsenic (As) exposure is a significant worldwide environmental health concern. Low dose, chronic arsenic exposure has been associated with a higher than normal risk of skin, lung, and bladder cancer, as well as cardiovascular disease and diabetes. While arsenic-induced biological changes play a role in disease pathology, little is known about the dynamic cellular changes resulting from arsenic exposure and withdrawal.

Results

In these studies, we sought to understand the molecular mechanisms behind the biological changes induced by arsenic exposure. A comprehensive global approach was employed to determine genome-wide changes to chromatin structure, transcriptome patterns and splicing patterns in response to chronic low dose arsenic and its subsequent withdrawal. Our results show that cells exposed to chronic low doses of sodium arsenite have distinct temporal and coordinated chromatin, gene expression, and miRNA changes consistent with differentiation and activation of multiple biochemical pathways. Most of these temporal patterns in gene expression are reversed when arsenic is withdrawn. However, some gene expression patterns remained altered, plausibly as a result of an adaptive response by cells. Additionally, the correlation of changes to gene expression and chromatin structure solidify the role of chromatin structure in gene regulatory changes due to arsenite exposure. Lastly, we show that arsenite exposure influences gene regulation both at the initiation of transcription as well as at the level of splicing.

Conclusions

Our results show that adaptation of cells to iAs-mediated EMT is coupled to changes in chromatin structure effecting differential transcriptional and splicing patterns of genes. These studies provide new insights into the mechanism of iAs-mediated pathology, which includes epigenetic chromatin changes coupled with changes to the transcriptome and splicing patterns of key genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1295-9) contains supplementary material, which is available to authorized users.