RNAseq analysis of treatment-dependent signaling changes during inflammation in a mouse cutaneous wound healing model.

BACKGROUND: Despite proven therapeutic effects in inflammatory conditions, the specific mechanisms of phytochemical therapies are not well understood. The transcriptome effects of Traumeel (Tr14), a multicomponent natural product, and diclofenac, a non-selective cyclooxygenase (COX) inhibitor, were compared in a mouse cutaneous wound healing model to identify both known and novel pathways for the anti-inflammatory effect of plant-derived natural products. METHODS: Skin samples from abraded mice were analyzed by single-molecule, amplification-free RNAseq transcript profiling at 7 points between 12 and 192 h after injury. Immediately after injury, the wounds were treated with either diclofenac, Tr14, or placebo control (n = 7 per group/time). RNAseq levels were compared between treatment and control at each time point using a systems biology approach. RESULTS: At early time points (12-36 h), both control and Tr14-treated wounds showed marked increase in the inducible COX2 enzyme mRNA, while diclofenac-treated wounds did not. Tr14, in contrast, modulated lipoxygenase transcripts, especially ALOX12/15, and phospholipases involved in arachidonate metabolism. Notably, Tr14 modulated a group of cell-type specific markers, including the T cell receptor, that could be explained by an overarching effect on the type of cells that were recruited into the wound tissue. CONCLUSIONS: Tr14 and diclofenac had very different effects on the COX/LOX synthetic pathway after cutaneous wounding. Tr14 allowed normal autoinduction of COX2 mRNA, but suppressed mRNA levels for key enzymes in the leukotriene synthetic pathway. Tr14 appeared to have a broad 'phytocellular' effect on the wound transcriptome by altering the balance of cell types present in the wound.

The Landscape of long noncoding RNA classification.

Advances in the depth and quality of transcriptome sequencing have revealed many new classes of long noncoding RNAs (lncRNAs). lncRNA classification has mushroomed to accommodate these new findings, even though the real dimensions and complexity of the noncoding transcriptome remain unknown. Although evidence of functionality of specific lncRNAs continues to accumulate, conflicting, confusing, and overlapping terminology has fostered ambiguity and lack of clarity in the field in general. The lack of fundamental conceptual unambiguous classification framework results in a number of challenges in the annotation and interpretation of noncoding transcriptome data. It also might undermine integration of the new genomic methods and datasets in an effort to unravel the function of lncRNA. Here, we review existing lncRNA classifications, nomenclature, and terminology. Then, we describe the conceptual guidelines that have emerged for their classification and functional annotation based on expanding and more comprehensive use of large systems biology-based datasets.

Functional annotation of the vlinc class of non-coding RNAs using systems biology approach.

Functionality of the non-coding transcripts encoded by the human genome is the coveted goal of the modern genomics research. While commonly relied on the classical methods of forward genetics, integration of different genomics datasets in a global Systems Biology fashion presents a more productive avenue of achieving this very complex aim. Here we report application of a Systems Biology-based approach to dissect functionality of a newly identified vast class of very long intergenic non-coding (vlinc) RNAs. Using highly quantitative FANTOM5 CAGE dataset, we show that these RNAs could be grouped into 1542 novel human genes based on analysis of insulators that we show here indeed function as genomic barrier elements. We show that vlinc RNAs genes likely function in cisto activate nearby genes. This effect while most pronounced in closely spaced vlinc RNA-gene pairs can be detected over relatively large genomic distances. Furthermore, we identified 101 vlinc RNA genes likely involved in early embryogenesis based on patterns of their expression and regulation. We also found another 109 such genes potentially involved in cellular functions also happening at early stages of development such as proliferation, migration and apoptosis. Overall, we show that Systems Biology-based methods have great promise for functional annotation of non-coding RNAs.

The Bromodomain protein BRD4 controls HOTAIR, a long noncoding RNA essential for glioblastoma proliferation.

Bromodomain and extraterminal (BET) domain proteins have emerged as promising therapeutic targets in glioblastoma and many other cancers. Small molecule inhibitors of BET bromodomain proteins reduce expression of several oncogenes required for Glioblastoma Multiforme (GBM) progression. However, the mechanism through which BET protein inhibition reduces GBM growth is not completely understood. Long noncoding RNAs (lncRNAs) are important epigenetic regulators with critical roles in cancer initiation and malignant progression, but mechanistic insight into their expression and regulation by BET bromodomain inhibitors remains elusive. In this study, we used Helicos single molecule sequencing to comprehensively profile lncRNAs differentially expressed in GBM, and we identified a subset of GBM-specific lncRNAs whose expression is regulated by BET proteins. Treatment of GBM cells with the BET bromdomain inhibitor I-BET151 reduced levels of the tumor-promoting lncRNA HOX transcript antisense RNA (HOTAIR) and restored the expression of several other GBM down-regulated lncRNAs. Conversely, overexpression of HOTAIR in conjunction with I-BET151 treatment abrogates the antiproliferative activity of the BET bromodomain inhibitor. Moreover, chromatin immunoprecipitation analysis demonstrated binding of Bromodomain Containing 4 (BRD4) to the HOTAIR promoter, suggesting that BET proteins can directly regulate lncRNA expression. Our data unravel a previously unappreciated mechanism through which BET proteins control tumor growth of glioblastoma cells and suggest that modulation of lncRNA networks may, in part, mediate the antiproliferative effects of many epigenetic inhibitors currently in clinical trials for cancer and other diseases.

Complex interaction between dengue virus replication and expression of miRNA-133a.

BACKGROUND: Dengue virus (DENV) is the most common vector-borne viral infection worldwide with approximately 390 million cases and 25,000 reported deaths each year. MicroRNAs (miRNAs) are small non-coding RNA molecules responsible for the regulation of gene expression by repressing mRNA translation or inducing mRNA degradation. Although miRNAs possess antiviral activity against many mammalian-infecting viruses, their involvement in DENV replication is poorly understood. METHODS: Here, we explored the relationship between DENV and cellular microRNAs using bioinformatics tools. We overexpressed miRNA-133a in Vero cells to test its role in DENV replication and analyzed its expression using RT-qPCR. Furthermore, the expression of polypyrimidine tract binding protein (PTB), a protein involved in DENV replication, was analyzed by western blot. In addition, we profiled miRNA-133a expression in Vero cells challenged with DENV-2, using Taqman miRNA. RESULTS: Bioinformatic analysis revealed that the 3' untranslated region (3'UTR) of the DENV genome of all four DENV serotypes is targeted by several cellular miRNAs, including miRNA-133a. We found that overexpression of synthetic miRNA-133a suppressed DENV replication. Additionally, we observed that PTB transcription , a miRNA-133a target, is down-regulated during DENV infection. Based in our results we propose that 3'UTR of DENV down-regulates endogenous expression of miRNA-133a in Vero cells during the first hours of infection. CONCLUSIONS: miRNA-133a regulates DENV replication possibly through the modulation of a host factor such as PTB. Further investigations are needed to verify whether miRNA-133a has an anti-DENV effect in vivo.

Regulation of chromatin structure by long noncoding RNAs: focus on natural antisense transcripts.

In the decade following the publication of the Human Genome, noncoding RNAs (ncRNAs) have reshaped our understanding of the broad landscape of genome regulation. During this period, natural antisense transcripts (NATs), which are transcribed from the opposite strand of either protein or non-protein coding genes, have vaulted to prominence. Recent findings have shown that NATs can exert their regulatory functions by acting as epigenetic regulators of gene expression and chromatin remodeling. Here, we review recent work on the mechanisms of epigenetic modifications by NATs and their emerging role as master regulators of chromatin states. Unlike other long ncRNAs, antisense RNAs usually regulate their counterpart sense mRNA in cis by bridging epigenetic effectors and regulatory complexes at specific genomic loci. Understanding the broad range of effects of NATs will shed light on the complex mechanisms that regulate chromatin remodeling and gene expression in development and disease.

Dark matter RNA illuminates the puzzle of genome-wide association studies.

In the past decade, numerous studies have made connections between sequence variants in human genomes and predisposition to complex diseases. However, most of these variants lie outside of the charted regions of the human genome whose function we understand; that is, the sequences that encode proteins. Consequently, the general concept of a mechanism that translates these variants into predisposition to diseases has been lacking, potentially calling into question the validity of these studies. Here we make a connection between the growing class of apparently functional RNAs that do not encode proteins and whose function we do not yet understand (the so-called 'dark matter' RNAs) and the disease-associated variants. We review advances made in a different genomic mapping effort - unbiased profiling of all RNA transcribed from the human genome - and provide arguments that the disease-associated variants exert their effects via perturbation of regulatory properties of non-coding RNAs existing in mammalian cells.

On the importance of small changes in RNA expression.

The analysis of the differential expression of genes has been the key goal of many molecular biology methods for decades and will remain with us for decades to come. It constitutes a fundamental resource at our disposal for determining the relationship between products of transcription, biology and disease. The completed genome sequencing of many common species allowed microarrays and RNA sequencing (RNAseq) to become major tools in Systems Biology. However, we estimate that at least half of all experiments ignore transcripts that change less than some subjectively chosen threshold, typically around 2-3 fold. Here we show that a majority of the informative RNAs and differentially expressed transcripts can exhibit fold changes less than 2. We use highly quantitative single-molecule sequencing of total cellular RNA derived from a time course of inflammatory response, a process critical to a large number of diseases. Furthermore, we show that enrichment of biologically-relevant functions occurs even at very low fold changes in RNA levels. In addition, we show that most of the common statistical methods can reliably detect transcripts with low fold change when as few as 3 biological replicates are sequenced using single-molecule based RNAseq. In conclusion, given the prevalence of expression profiling in current research, the loss of data in half of all expression studies results in a significant, yet needless drain on the discovery process.

Intronic RNAs constitute the major fraction of the non-coding RNA in mammalian cells.

BACKGROUND: The function of RNA from the non-coding (the so called "dark matter") regions of the genome has been a subject of considerable recent debate. Perhaps the most controversy is regarding the function of RNAs found in introns of annotated transcripts, where most of the reads that map outside of exons are usually found. However, it has been reported that the levels of RNA in introns are minor relative to those of the corresponding exons, and that changes in the levels of intronic RNAs correlate tightly with that of adjacent exons. This would suggest that RNAs produced from the vast expanse of intronic space are just pieces of pre-mRNAs or excised introns en route to degradation. RESULTS: We present data that challenges the notion that intronic RNAs are mere by-standers in the cell. By performing a highly quantitative RNAseq analysis of transcriptome changes during an inflammation time course, we show that intronic RNAs have a number of features that would be expected from functional, standalone RNA species. We show that there are thousands of introns in the mouse genome that generate RNAs whose overall abundance, which changes throughout the inflammation timecourse, and other properties suggest that they function in yet unknown ways. CONCLUSIONS: So far, the focus of non-coding RNA discovery has shied away from intronic regions as those were believed to simply encode parts of pre-mRNAs. Results presented here suggest a very different situation--the sequences encoded in the introns appear to harbor a yet unexplored reservoir of novel, functional RNAs. As such, they should not be ignored in surveys of functional transcripts or other genomic studies.

Insights from the HuR-interacting transcriptome: ncRNAs, ubiquitin pathways, and patterns of secondary structure dependent RNA interactions.

The HuR protein regulates the expression of thousands of cellular transcripts by modulating mRNA splicing, trafficking, translation, and stability. Although it serves as a model of RNA-protein interactions, many features of HuR's interactions with RNAs remain unknown. In this report, we deployed the cryogenic RNA immunoprecipitation technique to analyze HuR-interacting RNAs with the Affymetrix all-exon microarray platform. We revealed several thousand novel HuR-interacting RNAs, including hundreds of non-coding RNAs such as natural antisense transcripts from stress responsive loci. To gain insight into the mechanisms of specificity and sensitivity of HuR's interaction with its target RNAs, we searched HuR-interacting RNAs for composite patterns of primary sequence and secondary structure. We provide evidence that secondary structures of 66-75 nucleotides enhance HuR's recognition of its specific RNA targets composed of short primary sequence patterns. We validated thousands of these RNAs by analysis of overlap with recently published findings, including HuR's interaction with RNAs in the pathways of RNA splicing and stability. Finally, we observed a striking enrichment for members of ubiquitin ligase pathways among the HuR-interacting mRNAs, suggesting a new role for HuR in the regulation of protein degradation to mirror its known function in protein translation.

Primary hepatocyte culture supports hepatitis C virus replication: a model for infection-associated hepatocarcinogenesis.

UNLABELLED: Analysis of progressive changes in hepatic gene expression that underlie hepatocarcinogenesis following hepatitis C virus (HCV) infection require examination of long-term cultures of normally differentiating primary human hepatocytes. We report a culture system of primary hepatocytes that support productive replication of infectious HCV. Hepatic functions were analyzed by reverse-transcription polymerase chain reaction amplification of total cell RNA from cultures maintained in serum-free defined medium for up to 190 days. Sustained hepatic function was assessed by expression of albumin, alpha-fetoprotein, cytochrome P4502E1, cytokeratin-18, type-1 collagen, transforming growth factor-beta 1, matrix metalloproteinase-2 (MMP-2), MMP-13, and interferon alpha-receptors 1 and 2. Normally differentiated human primary hepatocytes supported productive replication of infectious clones of HCV genotypes 1a, 1b, and 2a; virus infection was inhibited by antibodies against CD81 virus entry factor. Virus released into the culture media of HCV-infected primary hepatocytes repeatedly passage to naïve hepatocytes. Replication of the three HCV genotypes shows interferon sensitivity observed in natural infections. CONCLUSION: Sustained cultures of physiologic host cells for the propagation of infectious HCV strains should accelerate studies of host response to HCV infection and progressive liver disease.

The majority of total nuclear-encoded non-ribosomal RNA in a human cell is 'dark matter' un-annotated RNA.

BACKGROUND: Discovery that the transcriptional output of the human genome is far more complex than predicted by the current set of protein-coding annotations and that most RNAs produced do not appear to encode proteins has transformed our understanding of genome complexity and suggests new paradigms of genome regulation. However, the fraction of all cellular RNA whose function we do not understand and the fraction of the genome that is utilized to produce that RNA remain controversial. This is not simply a bookkeeping issue because the degree to which this un-annotated transcription is present has important implications with respect to its biologic function and to the general architecture of genome regulation. For example, efforts to elucidate how non-coding RNAs (ncRNAs) regulate genome function will be compromised if that class of RNAs is dismissed as simply 'transcriptional noise'. RESULTS: We show that the relative mass of RNA whose function and/or structure we do not understand (the so called 'dark matter' RNAs), as a proportion of all non-ribosomal, non-mitochondrial human RNA (mt-RNA), can be greater than that of protein-encoding transcripts. This observation is obscured in studies that focus only on polyA-selected RNA, a method that enriches for protein coding RNAs and at the same time discards the vast majority of RNA prior to analysis. We further show the presence of a large number of very long, abundantly-transcribed regions (100's of kb) in intergenic space and further show that expression of these regions is associated with neoplastic transformation. These overlap some regions found previously in normal human embryonic tissues and raises an interesting hypothesis as to the function of these ncRNAs in both early development and neoplastic transformation. CONCLUSIONS: We conclude that 'dark matter' RNA can constitute the majority of non-ribosomal, non-mitochondrial-RNA and a significant fraction arises from numerous very long, intergenic transcribed regions that could be involved in neoplastic transformation.

Propagation of kinetic uncertainties through a canonical topology of the TLR4 signaling network in different regions of biochemical reaction space.

BACKGROUND: Signal transduction networks represent the information processing systems that dictate which dynamical regimes of biochemical activity can be accessible to a cell under certain circumstances. One of the major concerns in molecular systems biology is centered on the elucidation of the robustness properties and information processing capabilities of signal transduction networks. Achieving this goal requires the establishment of causal relations between the design principle of biochemical reaction systems and their emergent dynamical behaviors. METHODS: In this study, efforts were focused in the construction of a relatively well informed, deterministic, non-linear dynamic model, accounting for reaction mechanisms grounded on standard mass action and Hill saturation kinetics, of the canonical reaction topology underlying Toll-like receptor 4 (TLR4)-mediated signaling events. This signaling mechanism has been shown to be deployed in macrophages during a relatively short time window in response to lipopolysaccharide (LPS) stimulation, which leads to a rapidly mounted innate immune response. An extensive computational exploration of the biochemical reaction space inhabited by this signal transduction network was performed via local and global perturbation strategies. Importantly, a broad spectrum of biologically plausible dynamical regimes accessible to the network in widely scattered regions of parameter space was reconstructed computationally. Additionally, experimentally reported transcriptional readouts of target pro-inflammatory genes, which are actively modulated by the network in response to LPS stimulation, were also simulated. This was done with the main goal of carrying out an unbiased statistical assessment of the intrinsic robustness properties of this canonical reaction topology. RESULTS: Our simulation results provide convincing numerical evidence supporting the idea that a canonical reaction mechanism of the TLR4 signaling network is capable of performing information processing in a robust manner, a functional property that is independent of the signaling task required to be executed. Nevertheless, it was found that the robust performance of the network is not solely determined by its design principle (topology), but this may be heavily dependent on the network's current position in biochemical reaction space. Ultimately, our results enabled us the identification of key rate limiting steps which most effectively control the performance of the system under diverse dynamical regimes. CONCLUSIONS: Overall, our in silico study suggests that biologically relevant and non-intuitive aspects on the general behavior of a complex biomolecular network can be elucidated only when taking into account a wide spectrum of dynamical regimes attainable by the system. Most importantly, this strategy provides the means for a suitable assessment of the inherent variational constraints imposed by the structure of the system when systematically probing its parameter space.

Cardiovascular inflammation and lesion cell apoptosis: a novel connection via the interferon-inducible immunoproteasome.

OBJECTIVE: Increasing evidence suggests that chronic inflammation contributes to atherogenesis, and that acute inflammatory events cause plaque rupture, thrombosis, and myocardial infarction. The present studies examined how inflammatory factors, such as interferon-gamma (IFNgamma), cause increased sensitivity to apoptosis in vascular lesion cells. METHODS AND RESULTS: Cells from the fibrous cap of human atherosclerotic lesions were sensitized by interferon-gamma (IFNgamma) to Fas-induced apoptosis, in a Bcl-X(L) reversible manner. Microarray profiling identified 72 INFgamma-induced transcripts with potential relevance to apoptosis. Half could be excluded because they were induced by IRF-1 overexpression, which did not sensitize to apoptosis. IFNgamma treatment strongly reduced Mcl-1, phospho-Bcl-2 (ser70), and phospho-Bcl-X(L) (ser62) protein levels. Candidate transcripts were modulated by siRNA, overexpression, or inhibitors to assess the effect on IFNgamma-induced Fas sensitivity. Surprisingly, siRNA knockdown of PSMB8 (LMP7), an "immunoproteasome" component, reversed IFNgamma-induced sensitivity to Fas ligation and prevented Fas/IFNgamma-induced degradation of Mcl-1, but did not protect p-Bcl-2 or p-Bcl-X(L). Proteasome inhibition markedly increased Mcl-1, p-Bcl-2, and p-Bcl-X(L) levels after IFNgamma treatment. CONCLUSIONS: Although critical for antigen presentation, the immunoproteasome appears to be a key link between inflammatory factors and the control of vascular cell apoptosis and may thus be an important factor in plaque rupture and myocardial infarction.

Biomarker discovery in attention deficit hyperactivity disorder: RNA sequencing of whole blood in discordant twin and case-controlled cohorts.

BACKGROUND: A variety of DNA-based methods have been applied to identify genetic markers of attention deficit hyperactivity disorder (ADHD), but the connection to RNA-based gene expression has not been fully exploited. METHODS: Using well defined cohorts of discordant, monozygotic twins from the Michigan State University Twin Registry, and case-controlled ADHD cases in adolescents, the present studies utilized advanced single molecule RNA sequencing to identify expressed changes in whole blood RNA in ADHD. Multiple analytical strategies were employed to narrow differentially expressed RNA targets to a small set of potential biomarkers of ADHD. RESULTS: RNA markers common to both the discordant twin study and case-controlled subjects further narrowed the putative targets, some of which had been previously associated with ADHD at the DNA level. The potential role of several differentially expressed genes, including ABCB5, RGS2, GAK, GIT1 and 3 members of the galactose metabolism pathway (GALE, GALT, GALK1) are substantiated by prior associations to ADHD and by established mechanistic connections to molecular pathways relevant to ADHD and behavioral control. CONCLUSIONS: The convergence of DNA, RNA, and metabolic data suggests these may be promising targets for diagnostics and therapeutics in ADHD.

Noncoding RNAs: couplers of analog and digital information in nervous system function?

The mammalian nervous system expresses numerous noncoding RNAs (ncRNAs). We propose that ncRNAs are capable of coupling the digital information universe of nucleic acids with the analog universe of cellular protein interactions. ncRNAs could contribute to the success of the organism's information processing in several ways. First, ncRNAs would allow for efficient coupling of energy with information, wherein less energy is required to represent and process more information, condensed in analog and digital form, into smaller spatial and temporal domains, ideal for the environments found in neural tissues. Second, ncRNAs would permit the rapid acquisition of information from the environment, along with the rapid flexible processing and elimination of that information when it is no longer necessary. Third, ncRNAs would facilitate accelerated evolution of an organism's information content and functional computational systems. This emerging panorama might open new dimensions of information processing in the nervous system.

Identification of novel GLI1 target genes and regulatory circuits in human cancer cells.

Hedgehog (HH) signaling is involved in many physiological processes, and pathway deregulation can result in a wide range of malignancies. Glioma-associated oncogene 1 (GLI1) is a transcription factor and a terminal effector of the HH cascade. Despite its crucial role in tumorigenesis, our understanding of the GLI1 cellular targets is quite limited. In this study, we identified multiple new GLI1 target genes using a combination of different genomic surveys and then subjected them to in-depth validation in human cancer cell lines. We were able to validate >90% of the new targets, which were enriched in functions involved in neurogenesis and regulation of transcription, in at least one type of follow-up experiment. Strikingly, we found that RNA editing of GLI1 can modulate effects on the targets. Furthermore, one of the top targets, FOXS1, a gene encoding a transcription factor previously implicated in nervous system development, was shown to act in a negative feedback loop limiting the cellular effects of GLI1 in medulloblastoma and rhabdomyosarcoma cells. Moreover, FOXS1 is both highly expressed and positively correlated with GLI1 in medulloblastoma samples of the Sonic HH subgroup, further arguing for the existence of FOXS1/GLI1 interplay in human tumors. Consistently, high FOXS1 expression predicts longer relapse-free survival in breast cancer. Overall, our findings open multiple new avenues in HH signaling pathway research and have potential for translational implications.

Meta-analysis of microarray results: challenges, opportunities, and recommendations for standardization.

Microarray profiling of gene expression is a powerful tool for discovery, but the ability to manage and compare the resulting data can be problematic. Biological, experimental, and technical variations between studies of the same phenotype/phenomena create substantial differences in results. The application of conventional meta-analysis to raw microarray data is complicated by differences in the type of microarray used, gene nomenclatures, species, and analytical methods. An alternative approach to combining multiple microarray studies is to compare the published gene lists which result from the investigators' analyses of the raw data, as implemented in Lists of Lists Annotated (LOLA: www.lola.gwu.edu) and L2L (depts.washington.edu/l2l/). The present review considers both the potential value and the limitations of databasing and enabling the comparison of results from different microarray studies. Further, a major impediment to cross-study comparisons is the absence of a standard for reporting microarray study results. We propose a reporting standard: standard microarray results template (SMART), which will facilitate the integration of microarray studies.

Nuclear Factor 90(NF90) targeted to TAR RNA inhibits transcriptional activation of HIV-1.

BACKGROUND: Examination of host cell-based inhibitors of HIV-1 transcription may be important for attenuating viral replication. We describe properties of a cellular double-stranded RNA binding protein with intrinsic affinity for HIV-1 TAR RNA that interferes with Tat/TAR interaction and inhibits viral gene expression. RESULTS: Utilizing TAR affinity fractionation, North-Western blotting, and mobility-shift assays, we show that the C-terminal variant of nuclear factor 90 (NF90ctv) with strong affinity for the TAR RNA, competes with Tat/TAR interaction in vitro. Analysis of the effect of NF90ctv-TAR RNA interaction in vivo showed significant inhibition of Tat-transactivation of HIV-1 LTR in cells expressing NF90ctv, as well as changes in histone H3 lysine-4 and lysine-9 methylation of HIV chromatin that are consistent with the epigenetic changes in transcriptionally repressed gene. CONCLUSION: Structural integrity of the TAR element is crucial in HIV-1 gene expression. Our results show that perturbation Tat/TAR RNA interaction by the dsRNA binding protein is sufficient to inhibit transcriptional activation of HIV-1.

Deep Sequencing Transcriptome Analysis of Murine Wound Healing: Effects of a Multicomponent, Multitarget Natural Product Therapy-Tr14.

Wound healing involves an orchestrated response that engages multiple processes, such as hemostasis, cellular migration, extracellular matrix synthesis, and in particular, inflammation. Using a murine model of cutaneous wound repair, the transcriptome was mapped from 12 h to 8 days post-injury, and in response to a multicomponent, multi-target natural product, Tr14. Using single-molecule RNA sequencing (RNA-seq), there were clear temporal changes in known transcripts related to wound healing pathways, and additional novel transcripts of both coding and non-coding genes. Tr14 treatment modulated >100 transcripts related to key wound repair pathways, such as response to wounding, wound contraction, and cytokine response. The results provide the most precise and comprehensive characterization to date of the transcriptome's response to skin damage, repair, and multicomponent natural product therapy. By understanding the wound repair process, and the effects of natural products, it should be possible to intervene more effectively in diseases involving aberrant repair.