Transcriptomic characterization of lung pericytes in systemic sclerosis-associated pulmonary fibrosis.

Systemic sclerosis (SSc) is a chronic disease characterized by fibrosis and vascular abnormalities in the skin and internal organs, including the lung. SSc-associated pulmonary fibrosis (SSc-PF) is the leading cause of death in SSc patients. Pericytes are key regulators of vascular integrity and endothelial function. The role that pericytes play in SSc-PF remains unclear. We compared the transcriptome of pericytes from SSc-PF lungs (SScL) to pericytes from normal lungs (NORML). We identified 1,179 differentially expressed genes in SScL pericytes. Pathways enriched in SScL pericytes included prostaglandin, PI3K-AKT, calcium, and vascular remodeling signaling. Decreased cyclic AMP production and altered phosphorylation of AKT in response to prostaglandin E2 in SScL pericytes demonstrate the functional consequence of changes in the prostaglandin pathway that may contribute to fibrosis. The transcriptomic signature of SSc lung pericytes suggests that they promote vascular dysfunction and contribute to the loss of protection against lung inflammation and fibrosis.

First Characterization of the Transcriptome of Lung Fibroblasts of SSc Patients and Healthy Donors of African Ancestry.

Systemic sclerosis (SSc) is a connective tissue disorder that results in fibrosis of the skin and visceral organs. SSc-associated pulmonary fibrosis (SSc-PF) is the leading cause of death amongst SSc patients. Racial disparity is noted in SSc as African Americans (AA) have a higher frequency and severity of disease than European Americans (EA). Using RNAseq, we determined differentially expressed genes (DEGs; q < 0.1, log2FC > |0.6|) in primary pulmonary fibroblasts from SSc lungs (SScL) and normal lungs (NL) of AA and EA patients to characterize the unique transcriptomic signatures of AA-NL and AA-SScL fibroblasts using systems-level analysis. We identified 69 DEGs in "AA-NL vs. EA-NL" and 384 DEGs in "AA-SScL vs. EA-SScL" analyses, and a comparison of disease mechanisms revealed that only 7.5% of DEGs were commonly deregulated in AA and EA patients. Surprisingly, we also identified an SSc-like signature in AA-NL fibroblasts. Our data highlight differences in disease mechanisms between AA and EA SScL fibroblasts and suggest that AA-NL fibroblasts are in a "pre-fibrosis" state, poised to respond to potential fibrotic triggers. The DEGs and pathways identified in our study provide a wealth of novel targets to better understand disease mechanisms leading to racial disparity in SSc-PF and develop more effective and personalized therapies.

Investigating the effects of chronic low-dose radiation exposure in the liver of a hypothermic zebrafish model.

Mankind's quest for a manned mission to Mars is placing increased emphasis on the development of innovative radio-protective countermeasures for long-term space travel. Hibernation confers radio-protective effects in hibernating animals, and this has led to the investigation of synthetic torpor to mitigate the deleterious effects of chronic low-dose-rate radiation exposure. Here we describe an induced torpor model we developed using the zebrafish. We explored the effects of radiation exposure on this model with a focus on the liver. Transcriptomic and behavioural analyses were performed. Radiation exposure resulted in transcriptomic perturbations in lipid metabolism and absorption, wound healing, immune response, and fibrogenic pathways. Induced torpor reduced metabolism and increased pro-survival, anti-apoptotic, and DNA repair pathways. Coupled with radiation exposure, induced torpor led to a stress response but also revealed maintenance of DNA repair mechanisms, pro-survival and anti-apoptotic signals. To further characterise our model of induced torpor, the zebrafish model was compared with hepatic transcriptomic data from hibernating grizzly bears (Ursus arctos horribilis) and active controls revealing conserved responses in gene expression associated with anti-apoptotic processes, DNA damage repair, cell survival, proliferation, and antioxidant response. Similarly, the radiation group was compared with space-flown mice revealing shared changes in lipid metabolism.

Reduced Cathepsin L expression and secretion into the extracellular milieu contribute to lung fibrosis in systemic sclerosis.

OBJECTIVES: Lung fibrosis is the leading cause of death in SSc, with no cure currently available. Antifibrotic Endostatin (ES) production does not reach therapeutic levels in SSc patients, suggesting a deficit in its release from Collagen XVIII by the main cleavage enzyme, Cathepsin L (CTSL). Thus, elucidating a potential deficit in CTSL expression and activity unravels an underlying molecular cause for SSc-driven lung fibrosis. METHODS: Fibrosis was induced experimentally using TGF-ß in vitro, in primary human lung fibroblasts (pLFs), and ex vivo, in human lung tissues. ES and CTSL expression was quantified using ELISA, RT-qPCR, immunoblotting or immunofluorescence. Recombinant NC1-FLAG peptide was used to assess CTSL cleavage activity. CTSL expression was also compared between SSc vs normal (NL)-derived pLFs and lung tissues. RESULTS: ES levels were significantly reduced in media conditioned by TGF-ß-induced pLFs. TGF-ß-stimulated pLFs significantly reduced expression and secretion of CTSL into the extracellular matrix (ECM). CTSL was also sequestered in its inactive form into extracellular vesicles, further reducing its availability in the ECM. Media conditioned by TGF-ß-induced pLFs showed reduced cleavage of NC1-Flag and reduced release of the antifibrotic ES fragment. SSc-derived pLFs and lung tissues expressed significantly lower levels of CTSL compared with NL. CONCLUSIONS: Our findings identify CTSL as a protein protective against lung fibrosis via its activation of antifibrotic ES, and whose expression in SSc pLFs and lung tissues is suppressed. Identifying strategies to boost CTSL endogenous levels in SSc patients could serve as a viable therapeutic strategy.

miRNA and lncRNA Expression Networks Modulate Cell Cycle and DNA Repair Inhibition in Senescent Prostate Cells.

Cellular senescence is a state of permanent growth arrest that arises once cells reach the limit of their proliferative capacity. It creates an inflammatory microenvironment favouring the initiation and progression of various age-related diseases, including prostate cancer. Non-coding RNAs (ncRNAs) have emerged as important regulators of cellular gene expression. Nonetheless, very little is known about the interplay of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) and how deregulation of ncRNA networks promotes cellular senescence. To investigate this, human prostate epithelial cells were cultured through different passages until senescent, and their RNA was extracted and sequenced using RNA sequencing (RNAseq) and microRNA sequencing (miRNA-seq) miRNAseq. Differential expression (DE) gene analysis was performed to compare senescent and proliferating cells with Limma, miRNA-target interactions with multiMiR, lncRNA-target interactions using TCGA data and network evaluation with miRmapper. We found that miR-335-3p, miR-543 and the lncRNAs H19 and SMIM10L2A all play central roles in the regulation of cell cycle and DNA repair processes. Expression of most genes belonging to these pathways were down-regulated by senescence. Using the concept of network centrality, we determined the top 10 miRNAs and lncRNAs, with miR-335-3p and H19 identified as the biggest hubs for miRNAs and lncRNA respectively. These ncRNAs regulate key genes belonging to pathways involved in cell senescence and prostate cancer demonstrating their central role in these processes and opening the possibility for their use as biomarkers or therapeutic targets to mitigate against prostate ageing and carcinogenesis.

Antifibrotic factor KLF4 is repressed by the miR-10/TFAP2A/TBX5 axis in dermal fibroblasts: insights from twins discordant for systemic sclerosis.

OBJECTIVES: Systemic sclerosis (SSc) is a complex disease of unknown aetiology in which inflammation and fibrosis lead to multiple organ damage. There is currently no effective therapy that can halt the progression of fibrosis or reverse it, thus studies that provide novel insights into disease pathogenesis and identify novel potential therapeutic targets are critically needed. METHODS: We used global gene expression and genome-wide DNA methylation analyses of dermal fibroblasts (dFBs) from a unique cohort of twins discordant for SSc to identify molecular features of this pathology. We validated the findings using in vitro, ex vivo and in vivo models. RESULTS: Our results revealed distinct differentially expressed and methylated genes, including several transcription factors involved in stem cell differentiation and developmental programmes (KLF4, TBX5, TFAP2A and homeobox genes) and the microRNAs miR-10a and miR-10b which target several of these deregulated genes. We show that KLF4 expression is reduced in SSc dFBs and its expression is repressed by TBX5 and TFAP2A. We also show that KLF4 is antifibrotic, and its conditional knockout in fibroblasts promotes a fibrotic phenotype. CONCLUSIONS: Our data support a role for epigenetic dysregulation in mediating SSc susceptibility in dFBs, illustrating the intricate interplay between CpG methylation, miRNAs and transcription factors in SSc pathogenesis, and highlighting the potential for future use of epigenetic modifiers as therapies.

PDGF Promotes Dermal Fibroblast Activation via a Novel Mechanism Mediated by Signaling Through MCHR1.

Systemic sclerosis (SSc) is an autoimmune disease characterized by vasculopathy and excessive fibrosis of the skin and internal organs. To this day, no effective treatments to prevent the progression of fibrosis exist, and SSc patients have disabilities and reduced life expectancy. The need to better understand pathways that drive SSc and to find therapeutic targets is urgent. RNA sequencing data from SSc dermal fibroblasts suggested that melanin-concentrating hormone receptor 1 (MCHR1), one of the G protein-coupled receptors regulating emotion and energy metabolism, is abnormally deregulated in SSc. Platelet-derived growth factor (PDGF)-BB stimulation upregulated MCHR1 mRNA and protein levels in normal human dermal fibroblasts (NHDF), and MCHR1 silencing prevented the PDGF-BB-induced expression of the profibrotic factors transforming growth factor beta 1 (TGFß1) and connective tissue growth factor (CTGF). PDGF-BB bound MCHR1 in membrane fractions of NHDF, and the binding was confirmed using surface plasmon resonance (SPR). MCHR1 inhibition blocked PDGF-BB modulation of intracellular cyclic adenosine monophosphate (cAMP). MCHR1 silencing in NHDF reduced PDGF-BB signaling. In summary, MCHR1 promoted the fibrotic response in NHDF through modulation of TGFß1 and CTGF production, intracellular cAMP levels, and PDGF-BB-induced signaling pathways, suggesting that MCHR1 plays an important role in mediating the response to PDGF-BB and in the pathogenesis of SSc. Inhibition of MCHR1 should be considered as a novel therapeutic strategy in SSc-associated fibrosis.

Identification of Impacted Pathways and Transcriptomic Markers as Potential Mediators of Pulmonary Fibrosis in Transgenic Mice Expressing Human IGFBP5.

Pulmonary fibrosis is a serious disease characterized by extracellular matrix (ECM) component overproduction and remodeling. Insulin-like growth factor-binding protein 5 (IGFBP5) is a conserved member of the IGFBP family of proteins that is overexpressed in fibrotic tissues and promotes fibrosis. We used RNA sequencing (RNAseq) to identify differentially expressed genes (DEGs) between primary lung fibroblasts (pFBs) of homozygous (HOMO) transgenic mice expressing human IGFBP5 (hIGFBP5) and wild type mice (WT). The results of the differential expression analysis showed 2819 DEGs in hIGFBP5 pFBs. Functional enrichment analysis confirmed the pro-fibrotic character of IGFBP5 and revealed its impact on fundamental signaling pathways, including cytokine-cytokine receptor interaction, focal adhesion, AGE-RAGE signaling, calcium signaling, and neuroactive ligand-receptor interactions, to name a few. Noticeably, 7% of the DEGs in hIGFBP5-expressing pFBs are receptors and integrins. Furthermore, hub gene analysis revealed 12 hub genes including Fpr1, Bdkrb2, Mchr1, Nmur1, Cnr2, P2ry14, and Ptger3. Validation assays were performed to complement the RNAseq data. They confirmed significant differences in the levels of the corresponding proteins in cultured pFBs. Our study provides new insights into the molecular mechanism(s) of IGFBP5-associated pulmonary fibrosis through possible receptor interactions that drive fibrosis and tissue remodeling.

Induced Torpor as a Countermeasure for Low Dose Radiation Exposure in a Zebrafish Model.

The development of the Artemis programme with the goal of returning to the moon is spurring technology advances that will eventually take humans to Mars and herald a new era of interplanetary space travel. However, long-term space travel poses unique challenges including exposure to ionising radiation from galactic cosmic rays and potential solar particle events, exposure to microgravity and specific nutritional challenges arising from earth independent exploration. Ionising radiation is one of the major obstacles facing future space travel as it can generate oxidative stress and directly damage cellular structures such as DNA, in turn causing genomic instability, telomere shortening, extracellular-matrix remodelling and persistent inflammation. In the gastrointestinal tract (GIT) this can lead to leaky gut syndrome, perforations and motility issues, which impact GIT functionality and affect nutritional status. While current countermeasures such as shielding from the spacecraft can attenuate harmful biological effects, they produce harmful secondary particles that contribute to radiation exposure. We hypothesised that induction of a torpor-like state would confer a radioprotective effect given the evidence that hibernation extends survival times in irradiated squirrels compared to active controls. To test this hypothesis, a torpor-like state was induced in zebrafish using melatonin treatment and reduced temperature, and radiation exposure was administered twice over the course of 10 days. The protective effects of induced-torpor were assessed via RNA sequencing and qPCR of mRNA extracted from the GIT. Pathway and network analysis were performed on the transcriptomic data to characterise the genomic signatures in radiation, torpor and torpor + radiation groups. Phenotypic analyses revealed that melatonin and reduced temperature successfully induced a torpor-like state in zebrafish as shown by decreased metabolism and activity levels. Genomic analyses indicated that low dose radiation caused DNA damage and oxidative stress triggering a stress response, including steroidal signalling and changes to metabolism, and cell cycle arrest. Torpor attenuated the stress response through an increase in pro-survival signals, reduced oxidative stress via the oxygen effect and detection and removal of misfolded proteins. This proof-of-concept model provides compelling initial evidence for utilizing an induced torpor-like state as a potential countermeasure for radiation exposure.

Systemic sclerosis biomarkers detection in the secretome of TGFß1-activated primary human lung fibroblasts.

TGFß1 is a profibrotic mediator that contributes to a broad spectrum of pathologies, including systemic sclerosis-associated pulmonary fibrosis (SSc-PF). However, the secretome of TGFß1-stimulated primary human normal lung (NL) fibroblasts has not been well characterized. Using fluorescent 2-dimensional gel electrophoresis (2D-PAGE) and differential gel electrophoresis (DIGE) followed by Mass Spectrometry, we identified 37 differentially secreted proteins in the conditioned media of TGFß1-activated NL fibroblasts and generated a protein-protein association network of the TGFß1 secretome using STRING. Functional enrichment revealed that several biological processes and pathways characteristic of PF were enriched. Additionally, by comparing the TGFß1 secretome of NL fibroblasts to proteomic biomarkers from biological fluids of systemic sclerosis (SSc) patients, we identified 11 overlapping proteins. Together our data validate the TGFß1-induced secretome of NL fibroblasts as a valid in vitro model that reflects SSc biomarkers and identify potential therapeutic targets for SSc-PF. SIGNIFICANCE: All proteins secreted by fibroblasts into the extracellular space, representing the secretome, promote cell-to-cell communication as well as tissue homeostasis, immune mechanisms, developmental regulation, proteolysis, development of the extracellular matrix (ECM) and cell adhesion. Therefore, it is crucial to understand how TGFß1, a well-known profibrotic cytokine, modulates the secretome of pulmonary fibroblasts, and how the TGFß1-induced secretome resembles biomarkers in SSc. Using functional enrichment analysis, key pathways and hub proteins can be identified and studied as potential therapeutic targets for pulmonary fibrosis.

Prominence of IL6, IGF, TLR, and Bioenergetics Pathway Perturbation in Lung Tissues of Scleroderma Patients With Pulmonary Fibrosis.

Scleroderma-associated pulmonary fibrosis (SSc-PF) and idiopathic pulmonary fibrosis (IPF) are two of many chronic fibroproliferative diseases that are responsible for nearly 45% of all deaths in developed countries. While sharing several pathobiological characteristics, they also have very distinct features. Currently no effective anti-fibrotic treatments exist that can halt the progression of PF or reverse it. Our goal is to uncover potential gene targets for the development of anti-fibrotic therapies efficacious in both diseases, and those specific to SSc-PF, by identifying universal pathways and molecules driving fibrosis in SSc-PF and IPF tissues as well as those unique to SSc-PF. Using DNA microarray data, a meta-analysis of the differentially expressed (DE) genes in SSc-PF and IPF lung tissues (diseased vs. normal) was performed followed by a full systems level analysis of the common and unique transcriptomic signatures obtained. Protein-protein interaction networks were generated to identify hub proteins and explore the data using the centrality principle. Our results suggest that therapeutic strategies targeting IL6 trans-signaling, IGFBP2, IGFL2, and the coagulation cascade may be efficacious in both SSc-PF and IPF. Further, our data suggest that the expression of matrikine-producing collagens is also perturbed in PF. Lastly, an overall perturbation of bioenergetics, specifically between glycolysis and fatty acid metabolism, was uncovered in SSc-PF. Our findings provide insights into potential targets for the development of anti-fibrotic therapies that could be effective in both IPF and SSc-PF.

The effects of rosiglitazone on the neonatal rat cardiomyocyte transcriptome: a temporal analysis.

Aim: The objective was to determine via high-throughput RNA sequencing the temporal effects of rosiglitazone (Avandia®) on the neonatal rat ventricular myocyte transcriptome. Materials & methods: Neonatal rat ventricular myocytes (NRVMs) were exposed to rosiglitazone in vitro. Meta analyses utilized temporal comparisons of 0.5 h control versus 0.5 h treatment, 0.5 h treatment versus 24 h treatment and 24 h treatment versus 48 h treatment. Results: Time dependent responses were observed. At 0.5 h, the PI3K-AKT signaling pathway was impacted. At 24 h endoplasmic reticulum activity and protein degradation were altered. At 48 h, oxytocin signaling was perturbed. Conclusion: The effects of rosiglitazone occured early and increased in magnitude over time. A protective molecular response was triggered at 24 h and maintained until 48 h. In parallel, a response that can cause cardiac damage was activated. Our findings suggest that rosiglitazone has deleterious effects.

Genome-Wide Analysis of Low Dose Bisphenol-A (BPA) Exposure in Human Prostate Cells.

Endocrine disrupting compounds (EDCs) have the potential to cause adverse effects on wild-life and human health. Two important EDCs are the synthetic estrogen 17α-ethynylestradiol (EE2) and bisphenol-A (BPA) both of which are xenoestrogens (XEs) as they bind the estrogen receptor and dis-rupt estrogen physiology in mammals and other vertebrates. In the recent years the influence of XEs on oncogenes, specifically in relation to breast and prostate cancer has been the subject of considerable study. METHODOLOGY: In this study, healthy primary human prostate epithelial cells (PrECs) were exposed to environmentally relevant concentrations of BPA (5nM and 25nM BPA) and interrogated using a whole genome microarray. RESULTS: Exposure to 5 and 25nM BPA resulted in 7,182 and 7,650 differentially expressed (DE) genes, respectively in treated PrECs. Exposure to EE2 had the greatest effect on the PrEC transcriptome (8,891 DE genes). CONCLUSION: We dissected and investigated the nature of the non-estrogenic gene signature associated with BPA with a focus on transcripts relevant to epigenetic modifications. The expression of transcripts encoding nuclear hormone receptors as well as histone and DNA methylation, modifying enzymes were significantly perturbed by exposure to BPA.

Transcriptomic analysis of short-term 17α-ethynylestradiol exposure in two Californian sentinel fish species sardine (Sardinops sagax) and mackerel (Scomber japonicus).

Endocrine disrupting chemicals (EDCs) are substances which disrupt normal functioning of the endocrine system by interfering with hormone regulated physiological pathways. Aquatic environments provide the ultimate reservoir for many EDCs as they enter rivers and the ocean via effluent discharges and accumulate in sediments. One EDC widely dispersed in municipal wastewater effluent discharges is 17α-ethynylestradiol (EE2), which is one of the most widely prescribed medicines. EE2 is a bio-active estrogen employed in the majority of oral contraceptive pill formulations. As evidence of the health risks posed by EDCs mount, there is an urgent need to improve diagnostic tools for monitoring the effects of pollutants. As the cost of high throughput sequencing (HTS) diminishes, transcriptional profiling of an organism in response to EDC perturbation presents a cost-effective way of screening a wide range of endocrine responses. Coastal pelagic filter feeding fish species analyzed using HTS provide an excellent tool for EDC risk assessment in the marine environment. Unfortunately, there are limited genome sequence data and annotation for many of these species including Pacific sardine (Sardinops sagax) and chub mackerel (Scomber japonicus), which limits the utility of molecular tools such as HTS to interrogate the effects of endocrine disruption. In this study, we carried out RNA sequencing (RNAseq) of liver RNA harvested from wild sardine and mackerel exposed for 5 h under laboratory conditions to a concentration of 12.5 pM EE2 in the tank water. We developed an analytical framework for transcriptomic analyses of species with limited genomic information. EE2 exposure altered expression patterns of key genes involved in important metabolic and physiological processes. The systems approach presented here provides a powerful tool for obtaining a comprehensive picture of endocrine disruption in aquatic organisms.

Systems analysis of the liver transcriptome in adult male zebrafish exposed to the non-ionic surfactant nonylphenol.

Nonylphenol (NP) arises from the environmental degradation of nonylphenol ethoxylates. It is a ubiquitous environmental contaminant and has been detected at levels up to 167 nM in rivers in the United States. NP is an endocrine disruptor (ED) that can act as an agonist for estrogen receptors. The Adverse Outcome Pathway (AOP) framework defines an adverse outcome as the causal result of a series of molecular initiating events (MIEs) and key events (KEs) that lead to altered phenotypes. This study examined the liver transcriptome after a 21 day exposure to NP and 17ß-estradiol (E2) by exploiting the zebrafish (Danio rerio) as a systems toxicology model. The goal of this study was to tease out non-estrogenic genomic signatures associated with NP exposure using DNA microarray and RNA sequencing. Our experimental design included E2 as a positive and potent estrogenic control in order to effectively compare and contrast the 2 compounds. This approach allowed us to identify hepatic transcriptomic perturbations that could serve as MIEs for adverse health outcomes in response to NP. Our results revealed that exposure to NP was associated with differential expression (DE) of genes associated with the development of steatosis, disruption of metabolism, altered immune response, and metabolism of reactive oxygen species, further highlighting NP as a chemical of emerging concern (CEC).

Interplay Between MicroRNAs and Targeted Genes in Cellular Homeostasis of Adult Zebrafish (Danio rerio).

BACKGROUND: Cellular homeostasis is regulated by the intricate interplay between a plethora of signaling pathways and "energetic sensors" in organs. In order to maintain energy balance, induction or repression of metabolic pathways must be regulated and act in concert with the energetic demands of the cell at a given point in time. A new class of small noncoding RNAs, the microRNAs (miRNAs), has added yet further complexity to the control of metabolic homeostasis. OBJECTIVE: Understanding the damages induced by toxins in the liver and the intestine as well as the interplay between the miRNome and transcriptome first requires baseline characterization in these tissues in healthy animals under cellular homeostasis. METHODS: The liver (main site for detoxification) and the gut (primary exposure routes for contaminant exposure) were dissected out (wildtype fish), total and small RNA extracted, mRNA and miRNA libraries constructed and subjected to high throughput sequencing. Differential Expression (DE) analysis was performed comparing liver with gut and an "miRNA matrix" that integrates the miRNA-seq and mRNA-seq data was constructed. RESULTS: Both the miRNome and transcriptome of the liver and gut tissues were characterized and putative novel miRNAs were identified. Exploration of the "miRNA matrix" regulatory network revealed that miRNAs uniquely expressed in the liver or gut tissue regulated fundamental cellular processes important for both organs, and that commonly expressed miRNAs in both tissues regulated biological processes that were specific to either the liver or the gut. CONCLUSION: The result of our analyses revealed new insights into microRNA function in these tissues.

De Novo Hepatic Transcriptome Assembly and Systems Level Analysis of Three Species of Dietary Fish, Sardinops sagax, Scomber japonicus, and Pleuronichthys verticalis.

The monitoring of marine species as sentinels for ecosystem health has long been a valuable tool worldwide, providing insight into how both anthropogenic pollution and naturally occurring phenomena (i.e., harmful algal blooms) may lead to human and animal dietary concerns. The marine environments contain many contaminants of anthropogenic origin that have sufficient similarities to steroid and thyroid hormones, to potentially disrupt normal endocrine physiology in humans, fish, and other animals. An appropriate understanding of the effects of these endocrine disrupting chemicals (EDCs) on forage fish (e.g., sardine, anchovy, mackerel) can lead to significant insight into how these contaminants may affect local ecosystems in addition to their potential impacts on human health. With advancements in molecular tools (e.g., high-throughput sequencing, HTS), a genomics approach offers a robust toolkit to discover putative genetic biomarkers in fish exposed to these chemicals. However, the lack of available sequence information for non-model species has limited the development of these genomic toolkits. Using HTS and de novo assembly technology, the present study aimed to establish, for the first time for Sardinops sagax (Pacific sardine), Scomber japonicas (Pacific chub mackerel) and Pleuronichthys verticalis (hornyhead turbot), a de novo global transcriptome database of the liver, the primary organ involved in detoxification. The assembled transcriptomes provide a foundation for further downstream validation, comparative genomic analysis and biomarker development for future applications in ecotoxicogenomic studies, as well as environmental evaluation (e.g., climate change) and public health safety (e.g., dietary screening).

miRmapper: A Tool for Interpretation of miRNA⁻mRNA Interaction Networks.

It is estimated that 30% of all genes in the mammalian cells are regulated by microRNA (miRNAs). The most relevant miRNAs in a cellular context are not necessarily those with the greatest change in expression levels between healthy and diseased tissue. Differentially expressed (DE) miRNAs that modulate a large number of messenger RNA (mRNA) transcripts ultimately have a greater influence in determining phenotypic outcomes and are more important in a global biological context than miRNAs that modulate just a few mRNA transcripts. Here, we describe the development of a tool, "miRmapper", which identifies the most dominant miRNAs in a miRNA⁻mRNA network and recognizes similarities between miRNAs based on commonly regulated mRNAs. Using a list of miRNA⁻target gene interactions and a list of DE transcripts, miRmapper provides several outputs: (1) an adjacency matrix that is used to calculate miRNA similarity utilizing the Jaccard distance; (2) a dendrogram and (3) an identity heatmap displaying miRNA clusters based on their effect on mRNA expression; (4) a miRNA impact table and (5) a barplot that provides a visual illustration of this impact. We tested this tool using nonmetastatic and metastatic bladder cancer cell lines and demonstrated that the most relevant miRNAs in a cellular context are not necessarily those with the greatest fold change. Additionally, by exploiting the Jaccard distance, we unraveled novel cooperative interactions between miRNAs from independent families in regulating common target mRNAs; i.e., five of the top 10 miRNAs act in synergy.

Systems Analysis of the Liver Transcriptome in Adult Male Zebrafish Exposed to the Plasticizer (2-Ethylhexyl) Phthalate (DEHP).

The organic compound diethylhexyl phthalate (DEHP) represents a high production volume chemical found in cosmetics, personal care products, laundry detergents, and household items. DEHP, along with other phthalates causes endocrine disruption in males. Exposure to endocrine disrupting chemicals has been linked to the development of several adverse health outcomes with apical end points including Non-Alcoholic Fatty Liver Disease (NAFLD). This study examined the adult male zebrafish (Danio rerio) transcriptome after exposure to environmental levels of DEHP and 17α-ethinylestradiol (EE2) using both DNA microarray and RNA-sequencing technologies. Our results show that exposure to DEHP is associated with differentially expressed (DE) transcripts associated with the disruption of metabolic processes in the liver, including perturbation of five biological pathways: 'FOXA2 and FOXA3 transcription factor networks', 'Metabolic pathways', 'metabolism of amino acids and derivatives', 'metabolism of lipids and lipoproteins', and 'fatty acid, triacylglycerol, and ketone body metabolism'. DE transcripts unique to DEHP exposure, not observed with EE2 (i.e. non-estrogenic effects) exhibited a signature related to the regulation of transcription and translation, and ruffle assembly and organization. Collectively our results indicate that exposure to low DEHP levels modulates the expression of liver genes related to fatty acid metabolism and the development of NAFLD.

The Plasticizer Bisphenol A Perturbs the Hepatic Epigenome: A Systems Level Analysis of the miRNome.

Ubiquitous exposure to bisphenol A (BPA), an endocrine disruptor (ED), has raised concerns for both human and ecosystem health. Epigenetic factors, including microRNAs (miRNAs), are key regulators of gene expression during cancer. The effect of BPA exposure on the zebrafish epigenome remains poorly characterized. Zebrafish represents an excellent model to study cancer as the organism develops a disease that resembles human cancer. Using zebrafish as a systems toxicology model, we hypothesized that chronic BPA-exposure impacts the miRNome in adult zebrafish and establishes an epigenome more susceptible to cancer development. After a 3 week exposure to 100 nM BPA, RNA from the liver was extracted to perform high throughput mRNA and miRNA sequencing. Differential expression (DE) analyses comparing BPA-exposed to control specimens were performed using established bioinformatics pipelines. In the BPA-exposed liver, 6188 mRNAs and 15 miRNAs were differently expressed (q ≤ 0.1). By analyzing human orthologs of the DE zebrafish genes, signatures associated with non-alcoholic fatty liver disease (NAFLD), oxidative phosphorylation, mitochondrial dysfunction and cell cycle were uncovered. Chronic exposure to BPA has a significant impact on the liver miRNome and transcriptome in adult zebrafish with the potential to cause adverse health outcomes including cancer.