Pathway activity inference for multiclass disease classification through a mathematical programming optimisation framework

Background

Applying machine learning methods on microarray gene expression profiles for disease classification problems is a popular method to derive biomarkers, i.e. sets of genes that can predict disease state or outcome. Traditional approaches where expression of genes were treated independently suffer from low prediction accuracy and difficulty of biological interpretation. Current research efforts focus on integrating information on protein interactions through biochemical pathway datasets with expression profiles to propose pathway-based classifiers that can enhance disease diagnosis and prognosis. As most of the pathway activity inference methods in literature are either unsupervised or applied on two-class datasets, there is good scope to address such limitations by proposing novel methodologies.

Results

A supervised multiclass pathway activity inference method using optimisation techniques is reported. For each pathway expression dataset, patterns of its constituent genes are summarised into one composite feature, termed pathway activity, and a novel mathematical programming model is proposed to infer this feature as a weighted linear summation of expression of its constituent genes. Gene weights are determined by the optimisation model, in a way that the resulting pathway activity has the optimal discriminative power with regards to disease phenotypes. Classification is then performed on the resulting low-dimensional pathway activity profile.

Conclusions

The model was evaluated through a variety of published gene expression profiles that cover different types of disease. We show that not only does it improve classification accuracy, but it can also perform well in multiclass disease datasets, a limitation of other approaches from the literature. Desirable features of the model include the ability to control the maximum number of genes that may participate in determining pathway activity, which may be pre-specified by the user. Overall, this work highlights the potential of building pathway-based multi-phenotype classifiers for accurate disease diagnosis and prognosis problems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0390-2) contains supplementary material, which is available to authorized users.     

Methodology for the inference of gene function from phenotype data

Background

Biomedical ontologies are increasingly instrumental in the advancement of biological research primarily through their use to efficiently consolidate large amounts of data into structured, accessible sets. However, ontology development and usage can be hampered by the segregation of knowledge by domain that occurs due to independent development and use of the ontologies. The ability to infer data associated with one ontology to data associated with another ontology would prove useful in expanding information content and scope. We here focus on relating two ontologies: the Gene Ontology (GO), which encodes canonical gene function, and the Mammalian Phenotype Ontology (MP), which describes non-canonical phenotypes, using statistical methods to suggest GO functional annotations from existing MP phenotype annotations. This work is in contrast to previous studies that have focused on inferring gene function from phenotype primarily through lexical or semantic similarity measures.

Results

We have designed and tested a set of algorithms that represents a novel methodology to define rules for predicting gene function by examining the emergent structure and relationships between the gene functions and phenotypes rather than inspecting the terms semantically. The algorithms inspect relationships among multiple phenotype terms to deduce if there are cases where they all arise from a single gene function.We apply this methodology to data about genes in the laboratory mouse that are formally represented in the Mouse Genome Informatics (MGI) resource. From the data, 7444 rule instances were generated from five generalized rules, resulting in 4818 unique GO functional predictions for 1796 genes.

Conclusions

We show that our method is capable of inferring high-quality functional annotations from curated phenotype data. As well as creating inferred annotations, our method has the potential to allow for the elucidation of unforeseen, biologically significant associations between gene function and phenotypes that would be overlooked by a semantics-based approach. Future work will include the implementation of the described algorithms for a variety of other model organism databases, taking full advantage of the abundance of available high quality curated data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0405-z) contains supplementary material, which is available to authorized users.     

PhenUMA: a tool for integrating the biomedical relationships among genes and diseases

Background

Several types of genetic interactions in humans can be directly or indirectly associated with the causal effects of mutations. These interactions are usually based on their co-associations to biological processes, coexistence in cellular locations, coexpression in cell lines, physical interactions and so on. In addition, pathological processes can present similar phenotypes that have mutations either in the same genomic location or in different genomic regions. Therefore, integrative resources for all of these complex interactions can help us prioritize the relationships between genes and diseases that are most deserving to be studied by researchers and physicians.

Results

PhenUMA is a web application that displays biological networks using information from biomedical and biomolecular data repositories. One of its most innovative features is to combine the benefits of semantic similarity methods with the information taken from databases of genetic diseases and biological interactions. More specifically, this tool is useful in studying novel pathological relationships between functionally related genes, merging diseases into clusters that share specific phenotypes or finding diseases related to reported phenotypes.

Conclusions

This framework builds, analyzes and visualizes networks based on both functional and phenotypic relationships. The integration of this information helps in the discovery of alternative pathological roles of genes, biological functions and diseases. PhenUMA represents an advancement toward the use of new technologies for genomics and personalized medicine.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0375-1) contains supplementary material, which is available to authorized users.     

A knowledge-based weighting framework to boost the power of genome-wide association studies

Background

We are moving to second-wave analysis of genome-wide association studies (GWAS), characterized by comprehensive bioinformatical and statistical evaluation of genetic associations. Existing biological knowledge is very valuable for GWAS, which may help improve their detection power particularly for disease susceptibility loci of moderate effect size. However, a challenging question is how to utilize available resources that are very heterogeneous to quantitatively evaluate the statistic significances.

Methodology/Principal Findings

We present a novel knowledge-based weighting framework to boost power of the GWAS and insightfully strengthen their explorative performance for follow-up replication and deep sequencing. Built upon diverse integrated biological knowledge, this framework directly models both the prior functional information and the association significances emerging from GWAS to optimally highlight single nucleotide polymorphisms (SNPs) for subsequent replication. In the theoretical calculation and computer simulation, it shows great potential to achieve extra over 15% power to identify an association signal of moderate strength or to use hundreds of whole-genome subjects fewer to approach similar power. In a case study on late-onset Alzheimer disease (LOAD) for a proof of principle, it highlighted some genes, which showed positive association with LOAD in previous independent studies, and two important LOAD related pathways. These genes and pathways could be originally ignored due to involved SNPs only having moderate association significance.

Conclusions/Significance

With user-friendly implementation in an open-source Java package, this powerful framework will provide an important complementary solution to identify more true susceptibility loci with modest or even small effect size in current GWAS for complex diseases.     

Evaluating the Good Ontology Design Guideline (GoodOD) with the Ontology Quality Requirements and Evaluation Method and Metrics (OQuaRE)

Objective

To (1) evaluate the GoodOD guideline for ontology development by applying the OQuaRE evaluation method and metrics to the ontology artefacts that were produced by students in a randomized controlled trial, and (2) informally compare the OQuaRE evaluation method with gold standard and competency questions based evaluation methods, respectively.

Background

In the last decades many methods for ontology construction and ontology evaluation have been proposed. However, none of them has become a standard and there is no empirical evidence of comparative evaluation of such methods. This paper brings together GoodOD and OQuaRE. GoodOD is a guideline for developing robust ontologies. It was previously evaluated in a randomized controlled trial employing metrics based on gold standard ontologies and competency questions as outcome parameters. OQuaRE is a method for ontology quality evaluation which adapts the SQuaRE standard for software product quality to ontologies and has been successfully used for evaluating the quality of ontologies.

Methods

In this paper, we evaluate the effect of training in ontology construction based on the GoodOD guideline within the OQuaRE quality evaluation framework and compare the results with those obtained for the previous studies based on the same data.

Results

Our results show a significant effect of the GoodOD training over developed ontologies by topics: (a) a highly significant effect was detected in three topics from the analysis of the ontologies of untrained and trained students; (b) both positive and negative training effects with respect to the gold standard were found for five topics.

Conclusion

The GoodOD guideline had a significant effect over the quality of the ontologies developed. Our results show that GoodOD ontologies can be effectively evaluated using OQuaRE and that OQuaRE is able to provide additional useful information about the quality of the GoodOD ontologies.     

High-Content,Image-Based Screening for Drug Targets in Yeast

Background

Drug discovery and development are predicated on elucidation of the potential mechanisms of action and cellular targets of candidate chemical compounds. Recent advances in high-content imaging techniques allow simultaneous analysis of a range of cellular events. In this study, we propose a novel strategy to identify drug targets by combining genetic screening and high-content imaging in yeast.

Methodology

In this approach, we infer the cellular functions affected by candidate drugs by comparing morphologic changes induced by the compounds with the phenotypes of yeast mutants.

Conclusions

Using this method and four well-characterized reagents, we successfully identified previously known target genes of the compounds as well as other genes involved with functionally related cellular pathways. This is the first demonstration of a genetic high-content assay that can be used to identify drug targets based on morphologic phenotypes of a reference mutant panel.     

Turning text into research networks: information retrieval and computational ontologies in the creation of scientific databases

Background

Web-based, free-text documents on science and technology have been increasing growing on the web. However, most of these documents are not immediately processable by computers slowing down the acquisition of useful information. Computational ontologies might represent a possible solution by enabling semantically machine readable data sets. But, the process of ontology creation, instantiation and maintenance is still based on manual methodologies and thus time and cost intensive.

Method

We focused on a large corpus containing information on researchers, research fields, and institutions. We based our strategy on traditional entity recognition, social computing and correlation. We devised a semi automatic approach for the recognition, correlation and extraction of named entities and relations from textual documents which are then used to create, instantiate, and maintain an ontology.

Results

We present a prototype demonstrating the applicability of the proposed strategy, along with a case study describing how direct and indirect relations can be extracted from academic and professional activities registered in a database of curriculum vitae in free-text format. We present evidence that this system can identify entities to assist in the process of knowledge extraction and representation to support ontology maintenance. We also demonstrate the extraction of relationships among ontology classes and their instances.

Conclusion

We have demonstrated that our system can be used for the conversion of research information in free text format into database with a semantic structure. Future studies should test this system using the growing number of free-text information available at the institutional and national levels.     

Derivation and validation of clinical phenotypes for COPD: a systematic review

Background

The traditional classification of COPD, which relies solely on spirometry, fails to account for the complexity and heterogeneity of the disease. Phenotyping is a method that attempts to derive a single or combination of disease attributes that are associated with clinically meaningful outcomes. Deriving phenotypes entails the use of cluster analyses, and helps individualize patient management by identifying groups of individuals with similar characteristics. We aimed to systematically review the literature for studies that had derived such phenotypes using unsupervised methods.

Methods

Two independent reviewers systematically searched multiple databases for studies that performed validated statistical analyses, free of definitive pre-determined hypotheses, to derive phenotypes among patients with COPD. Data were extracted independently.

Results

9156 citations were retrieved, of which, 8 studies were included. The number of subjects ranged from 213 to 1543. Most studies appeared to be biased: patients were more likely males, with severe disease, and recruited in tertiary care settings. Statistical methods used to derive phenotypes varied by study. The number of phenotypes identified ranged from 2 to 5. Two phenotypes, with poor longitudinal health outcomes, were common across multiple studies: young patients with severe respiratory disease, few cardiovascular co-morbidities, poor nutritional status and poor health status, and a phenotype of older patients with moderate respiratory disease, obesity, cardiovascular and metabolic co-morbidities.

Conclusions

The recognition that two phenotypes of COPD were often reported may have clinical implications for altering the course of the disease. This review also provided important information on limitations of phenotype studies in COPD and the need for improvement in future studies.

Electronic supplementary material

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

Deletions of chromosomal regulatory boundaries are associated with congenital disease

Background

Recent data from genome-wide chromosome conformation capture analysis indicate that the human genome is divided into conserved megabase-sized self-interacting regions called topological domains. These topological domains form the regulatory backbone of the genome and are separated by regulatory boundary elements or barriers. Copy-number variations can potentially alter the topological domain architecture by deleting or duplicating the barriers and thereby allowing enhancers from neighboring domains to ectopically activate genes causing misexpression and disease, a mutational mechanism that has recently been termed enhancer adoption.

Results

We use the Human Phenotype Ontology database to relate the phenotypes of 922 deletion cases recorded in the DECIPHER database to monogenic diseases associated with genes in or adjacent to the deletions. We identify combinations of tissue-specific enhancers and genes adjacent to the deletion and associated with phenotypes in the corresponding tissue, whereby the phenotype matched that observed in the deletion. We compare this computationally with a gene-dosage pathomechanism that attempts to explain the deletion phenotype based on haploinsufficiency of genes located within the deletions. Up to 11.8% of the deletions could be best explained by enhancer adoption or a combination of enhancer adoption and gene-dosage effects.

Conclusions

Our results suggest that enhancer adoption caused by deletions of regulatory boundaries may contribute to a substantial minority of copy-number variation phenotypes and should thus be taken into account in their medical interpretation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0423-1) contains supplementary material, which is available to authorized users.     

Helveticoside is a biologically active component of the seed extract of Descurainia sophia and induces reciprocal gene regulation in A549 human lung cancer cells

Background

Although the pharmacological activities of the seed extract of Descurainia sophia have been proven to be useful against cough, asthma, and edema, the biologically active components, particularly at the molecular level, remain elusive. Therefore, we aimed to identify the active component of an ethanol extract of D. sophia seeds (EEDS) by applying a systematic genomic approach.

Results

After treatment with EEDS, the dose-dependently expressed genes in A549 cells were used to query the Connectivity map to determine which small molecules could closely mimic EEDS in terms of whole gene expression. Gene ontology and pathway analyses were also performed to identify the functional involvement of the drug responsive genes. In addition, interaction network and enrichment map assays were implemented to measure the functional network structure of the drug-responsive genes. A Connectivity map analysis of differentially expressed genes resulted in the discovery of helveticoside as a candidate drug that induces a similar gene expression pattern to EEDS. We identified the presence of helveticoside in EEDS and determined that helveticoside was responsible for the dose-dependent gene expression induced by EEDS. Gene ontology and pathway analyses revealed that the metabolism and signaling processes in A549 cells were reciprocally regulated by helveticoside and inter-connected as functional modules. Additionally, in an ontological network analysis, diverse cancer type-related genes were found to be associated with the biological functions regulated by helveticoside.

Conclusions

Using bioinformatic analyses, we confirmed that helveticoside is a biologically active component of EEDS that induces reciprocal regulation of metabolism and signaling processes. Our approach may provide novel insights to the herbal research field for identifying biologically active components from extracts.

Electronic supplementary material

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

Prioritizing causal disease genes using unbiased genomic features

Background

Cardiovascular disease (CVD) is the leading cause of death in the developed world. Human genetic studies, including genome-wide sequencing and SNP-array approaches, promise to reveal disease genes and mechanisms representing new therapeutic targets. In practice, however, identification of the actual genes contributing to disease pathogenesis has lagged behind identification of associated loci, thus limiting the clinical benefits.

Results

To aid in localizing causal genes, we develop a machine learning approach, Objective Prioritization for Enhanced Novelty (OPEN), which quantitatively prioritizes gene-disease associations based on a diverse group of genomic features. This approach uses only unbiased predictive features and thus is not hampered by a preference towards previously well-characterized genes. We demonstrate success in identifying genetic determinants for CVD-related traits, including cholesterol levels, blood pressure, and conduction system and cardiomyopathy phenotypes. Using OPEN, we prioritize genes, including FLNC, for association with increased left ventricular diameter, which is a defining feature of a prevalent cardiovascular disorder, dilated cardiomyopathy or DCM. Using a zebrafish model, we experimentally validate FLNC and identify a novel FLNC splice-site mutation in a patient with severe DCM.

Conclusion

Our approach stands to assist interpretation of large-scale genetic studies without compromising their fundamentally unbiased nature.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0534-8) contains supplementary material, which is available to authorized users.     

Identification of genes associated with shell color in the black-lipped pearl oyster,Pinctada margaritifera

Background

Color polymorphism in the nacre of pteriomorphian bivalves is of great interest for the pearl culture industry. The nacreous layer of the Polynesian black-lipped pearl oyster Pinctada margaritifera exhibits a large array of color variation among individuals including reflections of blue, green, yellow and pink in all possible gradients. Although the heritability of nacre color variation patterns has been demonstrated by experimental crossing, little is known about the genes involved in these patterns. In this study, we identify a set of genes differentially expressed among extreme color phenotypes of P. margaritifera using a suppressive and subtractive hybridization (SSH) method comparing black phenotypes with full and half albino individuals.

Results

Out of the 358 and 346 expressed sequence tags (ESTs) obtained by conducting two SSH libraries respectively, the expression patterns of 37 genes were tested with a real-time quantitative PCR (RT-qPCR) approach by pooling five individuals of each phenotype. The expression of 11 genes was subsequently estimated for each individual in order to detect inter-individual variation. Our results suggest that the color of the nacre is partially under the influence of genes involved in the biomineralization of the calcitic layer. A few genes involved in the formation of the aragonite tablets of the nacre layer and in the biosynthesis chain of melanin also showed differential expression patterns. Finally, high variability in gene expression levels were observed within the black phenotypes.

Conclusions

Our results revealed that three main genetic processes were involved in color polymorphisms: the biomineralization of the nacreous and calcitic layers and the synthesis of pigments such as melanin, suggesting that color polymorphism takes place at different levels in the shell structure. The high variability of gene expression found within black phenotypes suggests that the present work should serve as a basis for future studies exploring more thoroughly the expression patterns of candidate genes within black phenotypes with different dominant iridescent colors.

Electronic supplementary material

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

ProKinO: an ontology for integrative analysis of protein kinases in cancer

Background

Protein kinases are a large and diverse family of enzymes that are genomically altered in many human cancers. Targeted cancer genome sequencing efforts have unveiled the mutational profiles of protein kinase genes from many different cancer types. While mutational data on protein kinases is currently catalogued in various databases, integration of mutation data with other forms of data on protein kinases such as sequence, structure, function and pathway is necessary to identify and characterize key cancer causing mutations. Integrative analysis of protein kinase data, however, is a challenge because of the disparate nature of protein kinase data sources and data formats.

Results

Here, we describe ProKinO, a protein kinase-specific ontology, which provides a controlled vocabulary of terms, their hierarchy, and relationships unifying sequence, structure, function, mutation and pathway information on protein kinases. The conceptual representation of such diverse forms of information in one place not only allows rapid discovery of significant information related to a specific protein kinase, but also enables large-scale integrative analysis of protein kinase data in ways not possible through other kinase-specific resources. We have performed several integrative analyses of ProKinO data and, as an example, found that a large number of somatic mutations (∼288 distinct mutations) associated with the haematopoietic neoplasm cancer type map to only 8 kinases in the human kinome. This is in contrast to glioma, where the mutations are spread over 82 distinct kinases. We also provide examples of how ontology-based data analysis can be used to generate testable hypotheses regarding cancer mutations.

Conclusion

We present an integrated framework for large-scale integrative analysis of protein kinase data. Navigation and analysis of ontology data can be performed using the ontology browser available at: http://vulcan.cs.uga.edu/prokino.     

An Integrated Approach to Uncover Driver Genes in Breast Cancer Methylation Genomes

Background

Cancer cells typically exhibit large-scale aberrant methylation of gene promoters. Some of the genes with promoter methylation alterations play “driver” roles in tumorigenesis, whereas others are only “passengers”.

Results

Based on the assumption that promoter methylation alteration of a driver gene may lead to expression alternation of a set of genes associated with cancer pathways, we developed a computational framework for integrating promoter methylation and gene expression data to identify driver methylation aberrations of cancer. Applying this approach to breast cancer data, we identified many novel cancer driver genes and found that some of the identified driver genes were subtype-specific for basal-like, luminal-A and HER2+ subtypes of breast cancer.

Conclusion

The proposed framework proved effective in identifying cancer driver genes from genome-wide gene methylation and expression data of cancer. These results may provide new molecular targets for potential targeted and selective epigenetic therapy.     

Impact of Elevated Plasma Serotonin on Global Gene Expression of Murine Megakaryocytes

Background

Serotonin (5-HT) is a biogenic amine that also acts as a mitogen and a developmental signal early in rodent embryogenesis. Genetic and pharmacological disruption of 5-HT signaling causes various diseases and disorders via mediating central nervous system, cardiovascular system, and serious abnormalities on a growing embryo. Today, neither the effective modulators on 5-HT signaling pathways nor the genes affected by 5-HT signal are well known yet.

Methodology/Principal Findings

In an attempt to identify the genes altered by 5-HT signaling pathways, we analyzed the global gene expression via the Illumina array platform using the mouse WG-6 v2.0 Expression BeadChip containing 45,281 probe sets representing 30,854 genes in megakaryocytes isolated from mice infused with 5-HT or saline. We identified 723 differentially expressed genes of which 706 were induced and 17 were repressed by elevated plasma 5-HT.

Conclusions/Significance

Hierarchical gene clustering analysis was utilized to represent relations between groups and clusters. Using gene ontology mining tools and canonical pathway analyses, we identified multiple biological pathways that are regulated by 5-HT: (i) cytoskeletal remodeling, (ii) G-protein signaling, (iii) vesicular transport, and (iv) apoptosis and survival. Our data encompass the first extensive genome-wide based profiling in the progenitors of platelets in response to 5-HT elevation in vivo.     

Messina: A Novel Analysis Tool to Identify Biologically Relevant Molecules in Disease

Background

Morphologically similar cancers display heterogeneous patterns of molecular aberrations and follow substantially different clinical courses. This diversity has become the basis for the definition of molecular phenotypes, with significant implications for therapy. Microarray or proteomic expression profiling is conventionally employed to identify disease-associated genes, however, traditional approaches for the analysis of profiling experiments may miss molecular aberrations which define biologically relevant subtypes.

Methodology/Principal Findings

Here we present Messina, a method that can identify those genes that only sometimes show aberrant expression in cancer. We demonstrate with simulated data that Messina is highly sensitive and specific when used to identify genes which are aberrantly expressed in only a proportion of cancers, and compare Messina to contemporary analysis techniques. We illustrate Messina by using it to detect the aberrant expression of a gene that may play an important role in pancreatic cancer.

Conclusions/Significance

Messina allows the detection of genes with profiles typical of markers of molecular subtype, and complements existing methods to assist the identification of such markers. Messina is applicable to any global expression profiling data, and to allow its easy application has been packaged into a freely-available stand-alone software package.     

Expression-based in silico screening of candidate therapeutic compounds for lung adenocarcinoma

Background

Lung adenocarcinom (AC) is the most common form of lung cancer. Currently, the number of medical options to deal with lung cancer is very limited. In this study, we aimed to investigate potential therapeutic compounds for lung adenocarcinoma based on integrative analysis.

Methodology/Principal Findings

The candidate therapeutic compounds were identified in a two-step process. First, a meta-analysis of two published microarray data was conducted to obtain a list of 343 differentially expressed genes specific to lung AC. In the next step, expression profiles of these genes were used to query the Connectivity-Map (C-MAP) database to identify a list of compounds whose treatment reverse expression direction in various cancer cells. Several compounds in the categories of HSP90 inhibitor, HDAC inhibitor, PPAR agonist, PI3K inhibitor, passed our screening to be the leading candidates. On top of the list, three HSP90 inhibitors, i.e. 17-AAG (also known as tanespimycin), monorden, and alvespimycin, showed significant negative enrichment scores. Cytotoxicity as well as effects on cell cycle regulation and apoptosis were evaluated experimentally in lung adenocarcinoma cell line (A549 or GLC-82) with or without treatment with 17-AAG. In vitro study demonstrated that 17-AAG alone or in combination with cisplatin (DDP) can significantly inhibit lung adenocarcinoma cell growth by inducing cell cycle arrest and apoptosis.

Conclusions/Significance

We have used an in silico screening to identify compounds for treating lung cancer. One such compound 17-AAG demonstrated its anti-lung AC activity by inhibiting cell growth and promoting apoptosis and cell cycle arrest.     

Integrated analyses of zebrafish miRNA and mRNA expression profiles identify miR-29b and miR-223 as potential regulators of optic nerve regeneration

Background

Unlike mammals, zebrafish have the ability to regenerate damaged parts of their central nervous system (CNS) and regain functionality of the affected area. A better understanding of the molecular mechanisms involved in zebrafish regeneration may therefore provide insight into how CNS repair might be induced in mammals. Although many studies have described differences in gene expression in zebrafish during CNS regeneration, the regulatory mechanisms underpinning the differential expression of these genes have not been examined.

Results

We used microarrays to analyse and integrate the mRNA and microRNA (miRNA) expression profiles of zebrafish retina after optic nerve crush to identify potential regulatory mechanisms that underpin central nerve regeneration. Bioinformatic analysis identified 3 miRNAs and 657 mRNAs that were differentially expressed after injury. We then combined inverse correlations between our miRNA expression and mRNA expression, and integrated these findings with target predictions from TargetScan Fish to identify putative miRNA-gene target pairs. We focused on two over-expressed miRNAs (miR-29b and miR-223), and functionally validated seven of their predicted gene targets using RT-qPCR and luciferase assays to confirm miRNA-mRNA binding. Gene ontology analysis placed the miRNA-regulated genes (eva1a, layna, nefmb, ina, si:ch211-51a6.2, smoc1, sb:cb252) in key biological processes that included cell survival/apoptosis, ECM-cytoskeleton signaling, and heparan sulfate proteoglycan binding,

Conclusion

Our results suggest a key role for miR-29b and miR-223 in zebrafish regeneration. The identification of miRNA regulation in a zebrafish injury model provides a framework for future studies in which to investigate not only the cellular processes required for CNS regeneration, but also how these mechanisms might be regulated to promote successful repair and return of function in the injured mammalian brain.

Electronic supplementary material

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