Transcriptomic profiling of induced steatosis in human and mouse precision-cut liver slices.

There is a high need for predictive human ex vivo models for non-alcoholic fatty liver disease (NAFLD). About a decade ago, precision-cut liver slices (PCLSs) have been established as an ex vivo assay for humans and other organisms. In the present study, we use transcriptomics by RNASeq to profile a new human and mouse PCLSs based assay for steatosis in NAFLD. Steatosis as quantified by an increase of triglycerides after 48 h in culture, is induced by incremental supplementation of sugars (glucose and fructose), insulin, and fatty acids (palmitate, oleate). We mirrored the experimental design for human vs. mouse liver organ derived PCLSs and profiled each organ at eight different nutrient conditions after 24 h and 48 h time in culture. Thus, the provided data allows a comprehensive analysis of the donor, species, time, and nutrient factor specific regulation of gene expression in steatosis, despite the heterogeneity of the human tissue samples. Exemplified this is demonstrated by ranking homologous gene pairs by convergent or divergent expression pattern across nutrient conditions.

Transcriptome analysis of AAV-induced retinopathy models expressing human VEGF, TNF-α, and IL-6 in murine eyes.

Retinopathies are multifactorial diseases with complex pathologies that eventually lead to vision loss. Animal models facilitate the understanding of the pathophysiology and identification of novel treatment options. However, each animal model reflects only specific disease aspects and understanding of the specific molecular changes in most disease models is limited. Here, we conducted transcriptome analysis of murine ocular tissue transduced with recombinant Adeno-associated viruses (AAVs) expressing either human VEGF-A, TNF-α, or IL-6. VEGF expression led to a distinct regulation of extracellular matrix (ECM)-associated genes. In contrast, both TNF-α and IL-6 led to more comparable gene expression changes in interleukin signaling, and the complement cascade, with TNF-α-induced changes being more pronounced. Furthermore, integration of single cell RNA-Sequencing data suggested an increase of endothelial cell-specific marker genes by VEGF, while TNF-α expression increased the expression T-cell markers. Both TNF-α and IL-6 expression led to an increase in macrophage markers. Finally, transcriptomic changes in AAV-VEGF treated mice largely overlapped with gene expression changes observed in the oxygen-induced retinopathy model, especially regarding ECM components and endothelial cell-specific gene expression. Altogether, our study represents a valuable investigation of gene expression changes induced by VEGF, TNF-α, and IL-6 and will aid researchers in selecting appropriate animal models for retinopathies based on their agreement with the human pathophysiology.

Transcriptional comparison of adult human primary Retinal Pigment Epithelium, human pluripotent stem cell-derived Retinal Pigment Epithelium, and ARPE19 cells.

The therapeutic potential of pluripotent stem cells is great as they promise to usher in a new era of medicine where cells or organs may be prescribed to replace dysfunctional tissue. At the forefront are efforts in the eye to develop this technology as it lends itself to in vivo monitoring and sophisticated non-invasive imaging modalities. In the retina, retinal pigment epithelium (RPE) is the most promising replacement cell as it has a single layer, is relatively simple to transplant, and is associated with several eye diseases. However, after transplantation, the cells may transform and cause complications. This transformation may be partially due to incomplete maturation. With the goal of learning how to mature RPE, we compared induced pluripotent stem cell-derived RPE (iPSC-RPE) cells with adult human primary RPE (ahRPE) cells and the immortalized human ARPE-19 line. We cultured ARPE-19, iPSC-RPE, and ahRPE cells for one month, and evaluated morphology, RPE marker staining, and transepithelial electrical resistance (TEER) as quality control indicators. We then isolated RNA for bulk RNA-sequencing and DNA for genotyping. We genotyped ahRPE lines for the top age-related macular degeneration (AMD) and proliferative vitreoretinopathy (PVR) risk allele polymorphisms. Transcriptome data verified that both adult and iPSC-RPE exhibit similar RPE gene expression signatures, significantly higher than ARPE-19. In addition, in iPSC-RPE, genes relating to stem cell maintenance, retina development, and muscle contraction were significantly upregulated compared to ahRPE. We compared ahRPE to iPSC-RPE in a model of epithelial-mesenchymal transition (EMT) and observed an increased sensitivity of iPSC-RPE to producing contractile aggregates in vitro which resembles incident reports upon transplantation. P38 inhibition was capable of inhibiting iPSC-RPE-derived aggregates. In summary, we find that the transcriptomic signature of iPSC-RPE conveys an immature RPE state which may be ameliorated by targeting "immature" gene regulatory networks.

Molecular characterization and cell type composition deconvolution of fibrosis in NAFLD.

Non-alcoholic fatty liver disease (NAFLD) is the most common cause of liver disease worldwide. In adults with NAFLD, fibrosis can develop and progress to liver cirrhosis and liver failure. However, the underlying molecular mechanisms of fibrosis progression are not fully understood. Using total RNA-Seq, we investigated the molecular mechanisms of NAFLD and fibrosis. We sequenced liver tissue from 143 adults across the full spectrum of fibrosis stage including those with stage 4 fibrosis (cirrhosis). We identified gene expression clusters that strongly correlate with fibrosis stage including four genes that have been found consistently across previously published transcriptomic studies on NASH i.e. COL1A2, EFEMP2, FBLN5 and THBS2. Using cell type deconvolution, we estimated the loss of hepatocytes versus gain of hepatic stellate cells, macrophages and cholangiocytes with advancing fibrosis stage. Hepatocyte-specific functional analysis indicated increase of pro-apoptotic pathways and markers of bipotent hepatocyte/cholangiocyte precursors. Regression modelling was used to derive predictors of fibrosis stage. This study elucidated molecular and cell composition changes associated with increasing fibrosis stage in NAFLD and defined informative gene signatures for the disease.

High-throughput identification of synthetic riboswitches by barcode-free amplicon-sequencing in human cells.

Synthetic riboswitches mediating ligand-dependent RNA cleavage or splicing-modulation represent elegant tools to control gene expression in various applications, including next-generation gene therapy. However, due to the limited understanding of context-dependent structure-function relationships, the identification of functional riboswitches requires large-scale-screening of aptamer-effector-domain designs, which is hampered by the lack of suitable cellular high-throughput methods. Here we describe a fast and broadly applicable method to functionally screen complex riboswitch libraries (~1.8 × 104 constructs) by cDNA-amplicon-sequencing in transiently transfected and stimulated human cells. The self-barcoding nature of each construct enables quantification of differential mRNA levels without additional pre-selection or cDNA-manipulation steps. We apply this method to engineer tetracycline- and guanine-responsive ON- and OFF-switches based on hammerhead, hepatitis-delta-virus and Twister ribozymes as well as U1-snRNP polyadenylation-dependent RNA devices. In summary, our method enables fast and efficient high-throughput riboswitch identification, thereby overcoming a major hurdle in the development cascade for therapeutically applicable gene switches.

Transcriptomic characterisation of the optimised rat model of Walker 256 breast cancer cell-induced bone pain.

In patients with breast cancer, metastases of cancer cells to the axial skeleton may cause excruciating pain, particularly in the advanced stages. The current drug treatments available to alleviate this debilitating pain condition often lack efficacy and/or produce undesirable side effects. Preclinical animal models of cancer-induced bone pain are key to studying the mechanisms that cause this pain and for the success of drug discovery programs. In a previous study conducted in our laboratory, we validated and characterised the rat model of Walker 256 cell-induced bone pain, which displayed several key resemblances to the human pain condition. However, gene level changes that occur in the pathophysiology of cancer-induced bone pain in this preclinical model are unknown. Hence, in this study, we performed the transcriptomic characterisation of the Walker 256 cell line cultured in vitro to predict the molecular genetic profile of this cell line. We also performed transcriptomic characterisation of the Walker 256 cell-induced bone pain model in rats using the lumbar spinal cord and lumbar dorsal root ganglia tissues. Here we show that the Walker 256 cell line resembles the basal-B molecular subtype of human breast cancer cell lines. We also identify several genes that may underpin the progression of pain hypersensitivities in this condition, however, this needs further confirmatory studies. These transcriptomic insights have the potential to direct future studies aimed at identifying various mechanisms underpinning pain hypersensitivities in this model that may also assist in discovery of novel pain therapeutics for breast cancer-induced bone pain.

Molecular phenotyping of multiple mouse strains under metabolic challenge uncovers a role for Elovl2 in glucose-induced insulin secretion.

OBJECTIVE: In type 2 diabetes (T2D), pancreatic ß cells become progressively dysfunctional, leading to a decline in insulin secretion over time. In this study, we aimed to identify key genes involved in pancreatic beta cell dysfunction by analyzing multiple mouse strains in parallel under metabolic stress. METHODS: Male mice from six commonly used non-diabetic mouse strains were fed a high fat or regular chow diet for three months. Pancreatic islets were extracted and phenotypic measurements were recorded at 2 days, 10 days, 30 days, and 90 days to assess diabetes progression. RNA-Seq was performed on islet tissue at each time-point and integrated with the phenotypic data in a network-based analysis. RESULTS: A module of co-expressed genes was selected for further investigation as it showed the strongest correlation to insulin secretion and oral glucose tolerance phenotypes. One of the predicted network hub genes was Elovl2, encoding Elongase of very long chain fatty acids 2. Elovl2 silencing decreased glucose-stimulated insulin secretion in mouse and human ß cell lines. CONCLUSION: Our results suggest a role for Elovl2 in ensuring normal insulin secretory responses to glucose. Moreover, the large comprehensive dataset and integrative network-based approach provides a new resource to dissect the molecular etiology of ß cell failure under metabolic stress.

Genome-wide transcriptional profiling of skin and dorsal root ganglia after ultraviolet-B-induced inflammation.

Ultraviolet-B (UVB)-induced inflammation produces a dose-dependent mechanical and thermal hyperalgesia in both humans and rats, most likely via inflammatory mediators acting at the site of injury. Previous work has shown that the gene expression of cytokines and chemokines is positively correlated between species and that these factors can contribute to UVB-induced pain. In order to investigate other potential pain mediators in this model we used RNA-seq to perform genome-wide transcriptional profiling in both human and rat skin at the peak of hyperalgesia. In addition we have also measured transcriptional changes in the L4 and L5 DRG of the rat model. Our data show that UVB irradiation produces a large number of transcriptional changes in the skin: 2186 and 3888 genes are significantly dysregulated in human and rat skin, respectively. The most highly up-regulated genes in human skin feature those encoding cytokines (IL6 and IL24), chemokines (CCL3, CCL20, CXCL1, CXCL2, CXCL3 and CXCL5), the prostanoid synthesising enzyme COX-2 and members of the keratin gene family. Overall there was a strong positive and significant correlation in gene expression between the human and rat (R = 0.8022). In contrast to the skin, only 39 genes were significantly dysregulated in the rat L4 and L5 DRGs, the majority of which had small fold change values. Amongst the most up-regulated genes in DRG were REG3B, CCL2 and VGF. Overall, our data shows that numerous genes were up-regulated in UVB irradiated skin at the peak of hyperalgesia in both human and rats. Many of the top up-regulated genes were cytokines and chemokines, highlighting again their potential as pain mediators. However many other genes were also up-regulated and might play a role in UVB-induced hyperalgesia. In addition, the strong gene expression correlation between species re-emphasises the value of the UVB model as translational tool to study inflammatory pain.

A comparison of RNA-seq and exon arrays for whole genome transcription profiling of the L5 spinal nerve transection model of neuropathic pain in the rat.

BACKGROUND: The past decade has seen an abundance of transcriptional profiling studies of preclinical models of persistent pain, predominantly employing microarray technology. In this study we directly compare exon microarrays to RNA-seq and investigate the ability of both platforms to detect differentially expressed genes following nerve injury using the L5 spinal nerve transection model of neuropathic pain. We also investigate the effects of increasing RNA-seq sequencing depth. Finally we take advantage of the "agnostic" approach of RNA-seq to discover areas of expression outside of annotated exons that show marked changes in expression following nerve injury. RESULTS: RNA-seq and microarrays largely agree in terms of the genes called as differentially expressed. However, RNA-seq is able to interrogate a much larger proportion of the genome. It can also detect a greater number of differentially expressed genes than microarrays, across a wider range of fold changes and is able to assign a larger range of expression values to the genes it measures. The number of differentially expressed genes detected increases with sequencing depth. RNA-seq also allows the discovery of a number of genes displaying unusual and interesting patterns of non-exonic expression following nerve injury, an effect that cannot be detected using microarrays. CONCLUSION: We recommend the use of RNA-seq for future high-throughput transcriptomic experiments in pain studies. RNA-seq allowed the identification of a larger number of putative candidate pain genes than microarrays and can also detect a wider range of expression values in a neuropathic pain model. In addition, RNA-seq can interrogate the whole genome regardless of prior annotations, being able to detect transcription from areas of the genome not currently annotated as exons. Some of these areas are differentially expressed following nerve injury, and may represent novel genes or isoforms. We also recommend the use of a high sequencing depth in order to detect differential expression for genes with low levels of expression.

Phenocopy--a strategy to qualify chemical compounds during hit-to-lead and/or lead optimization.

A phenocopy is defined as an environmentally induced phenotype of one individual which is identical to the genotype-determined phenotype of another individual. The phenocopy phenomenon has been translated to the drug discovery process as phenotypes produced by the treatment of biological systems with new chemical entities (NCE) may resemble environmentally induced phenotypic modifications. Various new chemical entities exerting inhibition of the kinase activity of Transforming Growth Factor ß Receptor I (TGF-ßR1) were qualified by high-throughput RNA expression profiling. This chemical genomics approach resulted in a precise time-dependent insight to the TGF-ß biology and allowed furthermore a comprehensive analysis of each NCE's off-target effects. The evaluation of off-target effects by the phenocopy approach allows a more accurate and integrated view on optimized compounds, supplementing classical biological evaluation parameters such as potency and selectivity. It has therefore the potential to become a novel method for ranking compounds during various drug discovery phases.

Genetics and beyond--the transcriptome of human monocytes and disease susceptibility.

BACKGROUND: Variability of gene expression in human may link gene sequence variability and phenotypes; however, non-genetic variations, alone or in combination with genetics, may also influence expression traits and have a critical role in physiological and disease processes. METHODOLOGY/PRINCIPAL FINDINGS: To get better insight into the overall variability of gene expression, we assessed the transcriptome of circulating monocytes, a key cell involved in immunity-related diseases and atherosclerosis, in 1,490 unrelated individuals and investigated its association with >675,000 SNPs and 10 common cardiovascular risk factors. Out of 12,808 expressed genes, 2,745 expression quantitative trait loci were detected (P<5.78x10(-12)), most of them (90%) being cis-modulated. Extensive analyses showed that associations identified by genome-wide association studies of lipids, body mass index or blood pressure were rarely compatible with a mediation by monocyte expression level at the locus. At a study-wide level (P<3.9x10(-7)), 1,662 expression traits (13.0%) were significantly associated with at least one risk factor. Genome-wide interaction analyses suggested that genetic variability and risk factors mostly acted additively on gene expression. Because of the structure of correlation among expression traits, the variability of risk factors could be characterized by a limited set of independent gene expressions which may have biological and clinical relevance. For example expression traits associated with cigarette smoking were more strongly associated with carotid atherosclerosis than smoking itself. CONCLUSIONS/SIGNIFICANCE: This study demonstrates that the monocyte transcriptome is a potent integrator of genetic and non-genetic influences of relevance for disease pathophysiology and risk assessment.

Into the unknown: expression profiling without genome sequence information in CHO by next generation sequencing.

The arrival of next-generation sequencing (NGS) technologies has led to novel opportunities for expression profiling and genome analysis by utilizing vast amounts of short read sequence data. Here, we demonstrate that expression profiling in organisms lacking any genome or transcriptome sequence information is feasible by combining Illumina's mRNA-seq technology with a novel bioinformatics pipeline that integrates assembled and annotated Chinese hamster ovary (CHO) sequences with information derived from related organisms. We applied this pipeline to the analysis of CHO cells which were chosen as a model system owing to its relevance in the production of therapeutic proteins. Specifically, we analysed CHO cells undergoing butyrate treatment which is known to affect cell cycle regulation and to increase the specific productivity of recombinant proteins. By this means, we identified sequences for >13,000 CHO genes which added sequence information of approximately 5000 novel genes to the CHO model. More than 6000 transcript sequences are predicted to be complete, as they covered >95% of the corresponding mouse orthologs. Detailed analysis of selected biological functions such as DNA replication and cell cycle control, demonstrated the potential of NGS expression profiling in organisms without extended genome sequence to improve both data quantity and quality.

CHO gene expression profiling in biopharmaceutical process analysis and design.

Increase in both productivity and product yields in biopharmaceutical process development with recombinant protein producing mammalian cells can be mainly attributed to the advancements in cell line development, media, and process optimization. Only recently, genome-scale technologies enable a system-level analysis to elucidate the complex biomolecular basis of protein production in mammalian cells promising an increased process understanding and the deduction of knowledge-based approaches for further process optimization. Here, the use of gene expression profiling for the analysis of a low titer (LT) and high titer (HT) fed batch process using the same IgG producing CHO cell line was investigated. We found that gene expression (i) significantly differed in HT versus LT process conditions due to differences in applied chemically defined, serum-free media, (ii) changed over the time course of the fed batch processes, and that (iii) both metabolic pathways and 14 biological functions such as cellular growth or cell death were affected. Furthermore, detailed analysis of metabolism in a standard process format revealed the potential use of transcriptomics for rational media design as is shown for the case of lipid metabolism where the product titer could be increased by about 20% based on a lipid modified basal medium. The results demonstrate that gene expression profiling can be an important tool for mammalian biopharmaceutical process analysis and optimization.

miR-155 inhibition sensitizes CD4+ Th cells for TREG mediated suppression.

BACKGROUND: In humans and mice naturally occurring CD4(+)CD25(+) regulatory T cells (nTregs) are a thymus-derived subset of T cells, crucial for the maintenance of peripheral tolerance by controlling not only potentially autoreactive T cells but virtually all cells of the adaptive and innate immune system. Recent work using Dicer-deficient mice irrevocably demonstrated the importance of miRNAs for nTreg cell-mediated tolerance. PRINCIPAL FINDINGS: DNA-Microarray analyses of human as well as murine conventional CD4(+) Th cells and nTregs revealed a strong up-regulation of mature miR-155 (microRNA-155) upon activation in both populations. Studying miR-155 expression in FoxP3-deficient scurfy mice and performing FoxP3 ChIP-Seq experiments using activated human T lymphocytes, we show that the expression and maturation of miR-155 seem to be not necessarily regulated by FoxP3. In order to address the functional relevance of elevated miR-155 levels, we transfected miR-155 inhibitors or mature miR-155 RNAs into freshly-isolated human and mouse primary CD4(+) Th cells and nTregs and investigated the resulting phenotype in nTreg suppression assays. Whereas miR-155 inhibition in conventional CD4(+) Th cells strengthened nTreg cell-mediated suppression, overexpression of mature miR-155 rendered these cells unresponsive to nTreg cell-mediated suppression. CONCLUSION: Investigation of FoxP3 downstream targets, certainly of bound and regulated miRNAs revealed the associated function between the master regulator FoxP3 and miRNAs as regulators itself. miR-155 is shown to be crucially involved in nTreg cell mediated tolerance by regulating the susceptibility of conventional human as well as murine CD4(+) Th cells to nTreg cell-mediated suppression.