The Optimal Confidence Region for a Random Parameter

Under a two-level hierarchical model, suppose that the distribution of the random parameter is known or can be estimated well. Data are generated via a fixed, but unobservable realization of this parameter. In this paper, we derive the smallest confidence region of the random parameter under a joint Bayesian/frequentist paradigm. On average this optimal region can be much smaller than the corresponding Bayesian highest posterior density region. The new estimation procedure is appealing when one deals with data generated under a highly parallel structure, for example, data from a trial with a large number of clinical centers involved or genome-wide gene-expession data for estimating individual gene- or center-specific parameters simultaneously. The new proposal is illustrated with a typical microarray data set and its performance is examined via a small simulation study.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Integrative analysis of RUNX1 downstream pathways and target genes

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications.

The sequence of addition of IL-3 and IgE-dependent/IgE-independent stimuli regulates RANKL expression in human basophils

Background: Receptor Activator of Nuclear Factor kappaB Ligand (RANKL), a member of the TNF superfamily, contributes to the imbalance of bone resorption and immunoregulation in rheumatoid arthritis. In mice, collagen induced arthritis was exacerbated by IL-3 and anti-IgER antibodies, two mediators activating basophils that are known as effector cells of allergy. Interestingly, our unpublished microarray data revealed that IL-3 induces RANKL mRNA in human basophils. Here we further investigate under which conditions human basophils express surface and/or soluble RANKL. Methods: One part of purified human basophils was co-stimulated with IL-3 and either IgE-dependent or IgE-independent stimuli. The other part of purified basophils was first primed with IL-3 and subsequently triggered with IgE-dependent or IgE-independent stimuli. Expression of surface and soluble RANKL were detected by flow cytometry, ELISA and real-time PCR. Results: By flow cytometry we show that IL-3 induces de novo expression of surface RANKL on human basophils in a time and dose dependent manner. Co-stimulation of basophils with IL-3 and an IgE-dependent stimulus reduces IL-3-induced expression of surface RANKL in a dose dependent manner while IgE-independent stimuli have no effect. In contrast, both IgE-dependent and IgE-independent stimuli enhance expression of surface and soluble RANKL in basophils that were first primed with IL-3 and then triggered. Real-time PCR analysis shows that surface hRANKL1 and soluble hRANKL3 are induced by IL-3 and reduced by co-stimulation with IL-3 and an IgE-dependent stimulus and thus confirms our flow cytometry data. Conclusion: RANKL expression in human basophils is not only dependent on IL-3 and IgE-dependent/IgE-independent stimuli but also on the sequence of their addition to cell culture. Based on our data, we suggest that basophils might have previously unidentified functions in bone resorption or immunoregulation via RANKL.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

Loss of tumor suppressor mir-203 mediates overexpression of LIM and SH3 Protein 1 (LASP1) in high-risk prostate cancer thereby increasing cell proliferation and migration

Several studies have linked overexpression of the LIM and SH3 domain protein 1 (LASP1) to progression of breast, colon, liver, and bladder cancer. However, its expression pattern and role in human prostate cancer (PCa) remained largely undefined. Analysis of published microarray data revealed a significant overexpression of LASP1 in PCa metastases compared to parental primary tumors and normal prostate epithelial cells. Subsequent gene-set enrichment analysis comparing LASP1-high and -low PCa identified an association of LASP1 with genes involved in locomotory behavior and chemokine signaling. These bioinformatic predictions were confirmed in vitro as the inducible short hairpin RNA-mediated LASP1 knockdown impaired migration and proliferation in LNCaP prostate cancer cells. By immunohistochemical staining and semi-quantitative image analysis of whole tissue sections we found an enhanced expression of LASP1 in primary PCa and lymph node metastases over benign prostatic hyperplasia. Strong cytosolic and nuclear LASP1 immunoreactivity correlated with PSA progression. Conversely, qRT-PCR analyses for mir-203, which is a known translational suppressor of LASP1 in matched RNA samples revealed an inverse correlation of LASP1 protein and mir-203 expression. Collectively, our results suggest that loss of mir-203 expression and thus uncontrolled LASP1 overexpression might drive progression of PCa.

Transcriptional analysis of left-sided colitis, pancolitis, and ulcerative colitis-associated dysplasia.

BACKGROUND It is unknown why patients with extensive ulcerative colitis (UC) have a higher risk of colorectal cancer compared with patients with left-sided UC. This study characterizes the inflammatory processes in left-sided UC, pancolitis, and UC-associated dysplasia at the transcriptional level to identify potential biomarkers and transcripts of importance for the carcinogenic behavior of chronic inflammation. METHODS The Affymetrix GeneChip Human Genome U133 Plus 2.0 was applied on colonic biopsies from UC patients with left-sided UC, pancolitis, dysplasia, and controls. Reverse transcription polymerase chain reaction and immunohistochemistry were performed for validating selected transcripts in the initial cohort and in 2 independent cohorts of patients with UC. Microarray data were analyzed by principal component analysis, and reverse transcription polymerase chain reaction and immunohistochemistry data by the Wilcoxon's rank-sum test. RESULTS The principal component analysis results revealed separate clusters for left-sided UC, pancolitis, dysplasia, and controls. Close clustering of dysplastic and pancolitic samples indicated similarities in gene expression. Indeed, 101 and 656 parallel upregulated and downregulated transcripts, respectively, were identified in specimens from dysplasia and pancolitis. Validation of selected transcripts hereof identified insulin receptor alpha (INSRA) and MAP kinase interacting serine/threonine kinase 2 (MKNK2) with an enhanced expression in dysplasia compared with left-sided UC and controls, whereas laminin γ2 (LAMC2) was found with a lower expression in dysplasia compared with the remaining 3 groups. CONCLUSIONS This study demonstrates pancolitis and left-sided UC as distinct inflammatory processes at the transcriptional level, and identifies INSRA, MKNK2, and LAMC2 as potential critical transcripts in the inflammation-driven preneoplastic process of UC.

Metastatic lung cancer in a new mouse model inheriting p53 and latent K-ras mutations

Lung cancer is a devastating disease with very poor prognosis. The design of better treatments for patients would be greatly aided by mouse models that closely resemble the human disease. The most common type of human lung cancer is adenocarcinoma with frequent metastasis. Unfortunately, current models for this tumor are inadequate due to the absence of metastasis. Based on the molecular findings in human lung cancer and metastatic potential of osteosarcomas in mutant p53 mouse models, I hypothesized that mice with both K-ras and p53 missense mutations might develop metastatic lung adenocarcinomas. Therefore, I incorporated both K-rasLA1 and p53RI72HΔg alleles into mouse lung cells to establish a more faithful model for human lung adenocarcinoma and for translational and mechanistic studies. Mice with both mutations ( K-rasLA1/+ p53R172HΔg/+) developed advanced lung adenocarcinomas with similar histopathology to human tumors. These lung adenocarcinomas were highly aggressive and metastasized to multiple intrathoracic and extrathoracic sites in a pattern similar to that seen in lung cancer patients. This mouse model also showed gender differences in cancer related death and developed pleural mesotheliomas in 23.2% of them. In a preclinical study, the new drug Erlotinib (Tarceva) decreased the number and size of lung lesions in this model. These data demonstrate that this mouse model most closely mimics human metastatic lung adenocarcinoma and provides an invaluable system for translational studies. ^ To screen for important genes for metastasis, gene expression profiles of primary lung adenocarcinomas and metastases were analyzed. Microarray data showed that these two groups were segregated in gene expression and had 79 highly differentially expressed genes (more than 2.5 fold changes and p<0.001). Microarray data of Bub1b, Vimentin and CCAM1 were validated in tumors by quantitative real-time PCR (QPCR). Bub1b , a mitotic checkpoint gene, was overexpressed in metastases and this correlated with more chromosomal abnormalities in metastatic cells. Vimentin, a marker of epithelial-mesenchymal transition (EMT), was also highly expressed in metastases. Interestingly, Twist, a key EMT inducer, was also highly upregulated in metastases by QPCR, and this significantly correlated with the overexpression of Vimentin in the same tumors. These data suggest EMT occurs in lung adenocarcinomas and is a key mechanism for the development of metastasis in K-ras LA1/+ p53R172HΔg/+ mice. Thus, this mouse model provides a unique system to further probe the molecular basis of metastatic lung cancer.^

A mechanism of insulin-like growth factor binding protein 2-induced cell mobility in glioma

Insulin-like growth factor binding protein 2 (IGFBP2) is a protein known to be overexpressed in a majority of glioblastoma multiforme (GBM) tumors. While it is known the IGFBP2 is involved in promoting GBM tumor cell invasion, no mechanism exists for how the protein is involved in signal transduction pathways leading to enhanced cell invasion. ^ We follow up on preliminary microarray data on IGFBP2-overexpressing GBM cells and protein sequence analysis of IGFBP2 in generating the hypothesis that IGFBP2 interacts with integnn α5 in regulating cell mobility. Microarray data showing upregulation of integrin α5 by IGFBP2 is validated and evidence of protein-protein interaction between IGFBP2 and integrin α5 is found. The exact binding domain on IGFBP2 responsible for its interaction with integrin α5 is also determined, confirming our initial findings and reaffirming that the IGFBP2/integrin α5 interaction is specific. Disruption of this interaction resulted in attenuation of IGFBP2-enhanced cell mobility. Further, we found that cell mobility is only enhanced when IGFBP2 and integrin α5 are both overexpressed and able to interact with each other. ^ We also determined fibronectin to be a critical player in the activation of the IGFBP2/integrin α5 pathway. The activation of this pathway appears to be progressive and initiates once GBM cells have sufficiently established anchorage. ^

Involvement of lipocalin 2 in leukemia and breast cancer

Human lipocalin 2 is described as the neutrophil gelatinase-associated lipocalin (NGAL). The lipocalin 2 gene encodes a small, secreted glycoprotein that possesses a variety of functions, of which the best characterized function is organic iron binding activity. Elevated NGAL expression has been observed in many human cancers including breast, colorectal, pancreatic and ovarian cancers. I focused on the characterization of NGAL function in chronic myelogenous leukemia (CML) and breast cancer. Using the leukemic xenograft mouse model, we demonstrated that over-expression of NGAL in K562 cells, a leukemic cell line, led to a higher apoptotic rate and an atrophy phenotype in the spleen of inoculated mice compared to K562 cells alone. These results indicate that NGAL plays a primary role in suppressing hematopoiesis by inducing apoptosis within normal hematopoietic cells. In the breast cancer project, we analyzed two microarray data sets of breast cancer cell lines ( n = 54) and primary breast cancer samples (n = 318), and demonstrated that high NGAL expression is significantly correlated with several tumor characteristics, including negative estrogen receptor (ER) status, positive HER2 status, high tumor grade, and lymph node metastasis. Ectopic NGAL expression in non-aggressive (ZR75.1 and MCF7) cells led to aggressive tumor phenotypes in vitro and in vivo. Conversely, knockdown of NGAL expression in various breast cancer cell lines by shRNA lentiviral infection significantly decreased migration, invasion, and metastasis activities of tumor cells both in vitro and in vivo . It has been previously reported that transgenic mice with a mutation in the region of trans-membrane domain (V664E) of HER2 develop mammary tumors that progress to lung metastasis. However, we observed that genetic deletion of the 24p3 gene, a mouse homolog of NGAL, in HER2 transgenic mice by breeding with 24p3-null mice resulted in a significant delay of mammary tumor formation and reduction of lung metastasis. Strikingly, we also found that treatment with affinity purified 24p3 antibodies in the 4T1 breast cancer mice strongly reduced lung metastasis. Our studies provide evidence that NGAL plays a critical role in breast cancer development and progression, and thus NGAL has potential as a new therapeutic target in breast cancer.^

Regulation of two P450 isoforms: CYP2D6 and CYP4F11

Cytochromes P450 catalyze a monooxygenase reaction in which molecular oxygen is split and one oxygen atom is incorporated into the substrate. As a whole, P450 researchers have focused most of their attention on substrate metabolism and relatively little on how these enzymes are regulated. This study will focus on the regulation of two P450 isoforms known as, CYP2D6 and CYP4F11. ^ The human CYP2D gene locus contains two pseudogenes and one functional gene known as CYP2D6. This locus is highly polymorphic and produces several alternatively spliced transcripts from the pseudogene CYP2D7. My objective was to understand the role of SV5-in (splice variant 5), one of several alternative splice variants transcribed from the CYP2D7 pseudogene. My results indicate that SV5-in mRNA causes an increase in CYP2D6 protein levels and suggest that there is a role for SV5-in in regulation of CYP2D6 expression. ^ Second, CYP4F11 is a recently discovered and uncharacterized isoform, derived from the CYP4F subfamily. It metabolizes several clinically relevant drugs (i.e.—erythromycin and benzphetamine) and some endogenous inflammatory mediators (i.e.—LTB4). After evaluation of microarray data, I observed an increase in CYP4F11 mRNA levels from wild-type HCT116 cells compared to p53-null cells. Our objectives were to explore and understand this connection between p53 and CYP4F11. Microarray data were confirmed by Q-PCR, after which this effect was again observed at the protein level via Western blot and again at the promoter level via luciferase assay and chromatin immunoprecipitation. Our results indicate that p53 protein regulates expression of CYP4F11 mRNA and protein through CYP4F11 promoter binding (note that p53 binding to CYP4F11 DNA was not shown to be direct). These results signify a whole new level of regulation of drug metabolizing enzymes by p53. ^ An understanding of CYP4F11 regulation by p53 could help us understand another pathway leading to apoptosis or cell growth arrest. This can aid future drug studies and discover new drug metabolism pathways under the control of a tumor suppressor protein. An understanding of the CYP2D6 regulation pathway could illuminate the role of non-coding RNAs in the P450 field and potentially explain several inter-individual drug response variations observed in clinical medicine that are not yet completely explained by genotyping analysis. ^

Inflammatory Breast Cancer: The Immune Perspective

Inflammatory breast cancer (IBC) is the most insidious form of locally advanced disease. Although rare and less than 2% of all breast cancer, IBC is responsible for up to 10% of all breast cancer deaths. Despite the name, very little is known about the role of inflammation or immune mediators in IBC. Therefore, we analyzed blood samples from IBC patients and non-IBC patients, as well as healthy donor controls to establish an IBC-specific profile of peripheral blood leukocyte phenotype and function of T cells and dendritic cells and serum inflammatory cytokines. Emerging evidence suggests that host factors in the microenviromement may interact with underlying IBC genetics to promote the aggressive nature of the tumor. An integral part of the metastatic process involves epithelial to mesenchymal transition (EMT) where primary breast cancer cells gain motility and stem cell-like features that allow distant seeding. Interestingly, the IBC consortium microarray data found no clear evidence for EMT in IBC tumor tissues. It is becoming increasingly evident that inflammatory factors can induce EMT. However, it is unknown if EMT-inducing soluble factors secreted by activated immune cells in the IBC microenvironment canπ account for the absence of EMT in studies of the tumor cells themselves. We hypothesized that soluble factors from immune cells are capable of inducing EMT in IBC. We tested the ability of immune conditioned media to induce EMT in IBC cells. We found that soluble factors from activated immune cells are able to induce the expression of EMT-related factors in IBC cells along with increased migration and invasion. Specifically, the pro-inflammatory cytokines TNF-α, IL-6 and TGF-β were able to induce EMT and blocking these factors in conditioned media abated the induction of EMT. Surprisingly, unique to IBC cells, this process was related to increased levels of E-cadherin expression and adhesion, reminiscent of the characteristic tightly packed tumor emboli seen in IBC samples. This data offers insight into the unique pathology of IBC by suggesting that tumor immune interactions in the tumor microenvironment contribute to the aggressive nature of IBC implying that immune induced inflammation can be a novel therapeutic target. Specifically, we showed that soluble factors secreted by activated immune cells are capable of inducing EMT in IBC cells and may mediate the persistent E-cadherin expression observed in IBC. This data suggests that immune mediated inflammation may contribute to the highly aggressive nature of IBC and represents a potential therapeutic target that warrants further investigation.

Molecular and biochemical analyses of wounding in mesocarp of peach (Prunus persica L. Batsch) ripe fruits

The physiological and molecular responses of ripe fruit to wounding were evaluated in two peach (Prunus persica) varieties ('Glohaven', GH, melting and 'BigTop', BT, slow melting nectarine) by comparing mesocarp samples from wedges (as in minimal processing) and whole fruit as the control. Slight differences between the two varieties were detected in terms of ethylene production, whereas total phenol and flavonoid concentrations, and PPO and POD enzyme activities showed a general increase in wounded GH but not in BT. This was associated with the better appearance of the BT wedges at the end of the experimental period (72 h). Microarray (genome-wide ?PEACH3.0) analysis revealed that a total number of 2218 genes were differentially expressed (p < 0.01, log2 fold change expression ratio >1 or <-1) in GH 24 h after wounding compared to the control. This number was much lower (1208) in BT. According to the enrichment analysis, cell wall, plasma membrane, response to stress, secondary metabolic processes, oxygen binding were the GO categories over-represented among the GH up-regulated genes, whereas plasma membrane and response to endogenous stimulus were the categories over-represented among the down-regulated genes. Only 32 genes showed a common expression trend in the two varieties 24 h after wounding, whereas a total of 512 genes (with highly represented transcription factors), displayed opposite behavior. Quantitative RT-PCR analysis confirmed the microarray data for 18 out of a total of 20 genes selected. Specific WRKY, AP2/ERF and HSP20 genes were markedly up-regulated in wounded GH, indicating the activation of regulatory and signaling mechanisms probably related to different hormone categories. Compared to BT, the expression of specific genes involved in phenylpropanoid and triterpenoid biosynthetic pathways showed a more pronounced induction in GH, highlighting the difference between the two peach varieties in terms of molecular responses to wounding in the mesocarp tissue.

Biogeochemistry and microbiology of groundwater during acetate and bicarbonate amendment to an alluvial aquifer

These data are from a field experiment conducted in a shallow alluvial aquifer along the Colorado River in Rifle, Colorado, USA. In this experiment, bicarbonate-promoted uranium desorption and acetate amendment were combined and compared to an acetate amendment-only experiment in the same experimental plot. Data include names and location data for boreholes, geochemical data for all the boreholes between June 1, 2010 and January 1, 2011, microarray data provided as signal to noise ratio (SNR) for individual microarray probes, microarray data provided as signal to noise ratio (SNR) by Genus.

Transcriptional response to cardiac injury in the zebrafish: systematic identification of genes with highly concordant activity across in vivo models

BACKGROUND Zebrafish is a clinically-relevant model of heart regeneration. Unlike mammals, it has a remarkable heart repair capacity after injury, and promises novel translational applications. Amputation and cryoinjury models are key research tools for understanding injury response and regeneration in vivo. An understanding of the transcriptional responses following injury is needed to identify key players of heart tissue repair, as well as potential targets for boosting this property in humans. RESULTS We investigated amputation and cryoinjury in vivo models of heart damage in the zebrafish through unbiased, integrative analyses of independent molecular datasets. To detect genes with potential biological roles, we derived computational prediction models with microarray data from heart amputation experiments. We focused on a top-ranked set of genes highly activated in the early post-injury stage, whose activity was further verified in independent microarray datasets. Next, we performed independent validations of expression responses with qPCR in a cryoinjury model. Across in vivo models, the top candidates showed highly concordant responses at 1 and 3 days post-injury, which highlights the predictive power of our analysis strategies and the possible biological relevance of these genes. Top candidates are significantly involved in cell fate specification and differentiation, and include heart failure markers such as periostin, as well as potential new targets for heart regeneration. For example, ptgis and ca2 were overexpressed, while usp2a, a regulator of the p53 pathway, was down-regulated in our in vivo models. Interestingly, a high activity of ptgis and ca2 has been previously observed in failing hearts from rats and humans. CONCLUSIONS We identified genes with potential critical roles in the response to cardiac damage in the zebrafish. Their transcriptional activities are reproducible in different in vivo models of cardiac injury.