Dealing with Errors in QCA

This paper discusses five strategies to deal with five types of errors in Qualitative Comparative Analysis (QCA): condition errors, systematic errors, random errors, calibration errors, and deviant case errors. These strategies are the comparative inspection of complex, intermediary, and parsimonious solutions; the use of an adjustment factor, the use of probabilistic criteria, the test of the robustness of calibration parameters, and the use of a frequency threshold for observed combinations of conditions. The strategies are systematically reviewed, assessed, and evaluated as regards their applicability, advantages, limitations, and complementarities.

Simultaneous B(0)- and B(1)+-map acquisition for fast localized shim, frequency, and RF power determination in the heart at 3 T.

High-field (>or=3 T) cardiac MRI is challenged by inhomogeneities of both the static magnetic field (B(0)) and the transmit radiofrequency field (B(1)+). The inhomogeneous B fields not only demand improved shimming methods but also impede the correct determination of the zero-order terms, i.e., the local resonance frequency f(0) and the radiofrequency power to generate the intended local B(1)+ field. In this work, dual echo time B(0)-map and dual flip angle B(1)+-map acquisition methods are combined to acquire multislice B(0)- and B(1)+-maps simultaneously covering the entire heart in a single breath hold of 18 heartbeats. A previously proposed excitation pulse shape dependent slice profile correction is tested and applied to reduce systematic errors of the multislice B(1)+-map. Localized higher-order shim correction values including the zero-order terms for frequency f(0) and radiofrequency power can be determined based on the acquired B(0)- and B(1)+-maps. This method has been tested in 7 healthy adult human subjects at 3 T and improved the B(0) field homogeneity (standard deviation) from 60 Hz to 35 Hz and the average B(1)+ field from 77% to 100% of the desired B(1)+ field when compared to more commonly used preparation methods.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

On the use of spatially discrete data to compute energy and mass balance

In many practical applications the state of field soils is monitored by recording the evolution of temperature and soil moisture at discrete depths. We theoretically investigate the systematic errors that arise when mass and energy balances are computed directly from these measurements. We show that, even with no measurement or model errors, large residuals might result when finite difference approximations are used to compute fluxes and storage term. To calculate the limits set by the use of spatially discrete measurements on the accuracy of balance closure, we derive an analytical solution to estimate the residual on the basis of the two key parameters: the penetration depth and the distance between the measurements. When the thickness of the control layer for which the balance is computed is comparable to the penetration depth of the forcing (which depends on the thermal diffusivity and on the forcing period) large residuals arise. The residual is also very sensitive to the distance between the measurements, which requires accurately controlling the position of the sensors in field experiments. We also demonstrate that, for the same experimental setup, mass residuals are sensitively larger than the energy residuals due to the nonlinearity of the moisture transport equation. Our analysis suggests that a careful assessment of the systematic mass error introduced by the use of spatially discrete data is required before using fluxes and residuals computed directly from field measurements.

High precision anatomy for MEG.

Precise MEG estimates of neuronal current flow are undermined by uncertain knowledge of the head location with respect to the MEG sensors. This is either due to head movements within the scanning session or systematic errors in co-registration to anatomy. Here we show how such errors can be minimized using subject-specific head-casts produced using 3D printing technology. The casts fit the scalp of the subject internally and the inside of the MEG dewar externally, reducing within session and between session head movements. Systematic errors in matching to MRI coordinate system are also reduced through the use of MRI-visible fiducial markers placed on the same cast. Bootstrap estimates of absolute co-registration error were of the order of 1mm. Estimates of relative co-registration error were <1.5mm between sessions. We corroborated these scalp based estimates by looking at the MEG data recorded over a 6month period. We found that the between session sensor variability of the subject's evoked response was of the order of the within session noise, showing no appreciable noise due to between-session movement. Simulations suggest that the between-session sensor level amplitude SNR improved by a factor of 5 over conventional strategies. We show that at this level of coregistration accuracy there is strong evidence for anatomical models based on the individual rather than canonical anatomy; but that this advantage disappears for errors of greater than 5mm. This work paves the way for source reconstruction methods which can exploit very high SNR signals and accurate anatomical models; and also significantly increases the sensitivity of longitudinal studies with MEG.

Systematic difference between blood pressure readings caused by cuff type.

In this study we determine whether blood pressure readings using a cuff of fixed size systematically differed from readings made with a triple-bladder cuff (Tricuff) that automatically adjusts bladder width to arm circumference and assessed subsequent clinical and epidemiological effects. Blood pressure was measured with a standard cuff or a Tricuff in 454 patients visiting an outpatient clinic in the Seychelles (Indian Ocean). Overall means of within-individual standard cuff-Tricuff differences in systolic and diastolic blood pressures were examined in relation to arm circumference and sex. The standard cuff-Tricuff difference in systolic and diastolic blood pressures increased monotonically with circumference (from 4.7 +/- 0.8/3.2 +/- 0.7 mm Hg for arm circumference of 30 to 31 cm to 10.0 +/- 1.1/8.0 +/- 0.9 mm Hg for arm circumference &gt; or = 36 cm) and was larger in women than men. Multivariate linear regression indicated independent effects of arm circumference and sex. Forty percent of subjects with a diastolic blood pressure of &gt; or = 95 mm Hg measured with a standard cuff had values less than 95 mm Hg measured with a Tricuff. Extrapolation to the entire population of the Seychelles decreased the prevalence of blood pressure greater than or equal to 160/95 mm Hg by 11.5% and 24.0% in men and women, respectively, aged 35 to 64 years. The age-adjusted effect of body mass index on systolic and diastolic blood pressures decreased twofold using blood pressure readings made with a Tricuff instead of a standard cuff.(ABSTRACT TRUNCATED AT 250 WORDS)

Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

Performance of matrix-assisted laser desorption ionization-time of flight mass spectrometry for identification of bacterial strains routinely isolated in a clinical microbiology laboratory.

Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) has recently been introduced in diagnostic microbiology laboratories for the identification of bacterial and yeast strains isolated from clinical samples. In the present study, we prospectively compared MALDI-TOF MS to the conventional phenotypic method for the identification of routine isolates. Colonies were analyzed by MALDI-TOF MS either by direct deposition on the target plate or after a formic acid-acetonitrile extraction step if no valid result was initially obtained. Among 1,371 isolates identified by conventional methods, 1,278 (93.2%) were putatively identified to the species level by MALDI-TOF MS and 73 (5.3%) were identified to the genus level, but no reliable identification was obtained for 20 (1.5%). Among the 1,278 isolates identified to the species level by MALDI-TOF MS, 63 (4.9%) discordant results were initially identified. Most discordant results (42/63) were due to systematic database-related taxonomical differences, 14 were explained by poor discrimination of the MALDI-TOF MS spectra obtained, and 7 were due to errors in the initial conventional identification. An extraction step was required to obtain a valid MALDI-TOF MS identification for 25.6% of the 1,278 valid isolates. In conclusion, our results show that MALDI-TOF MS is a fast and reliable technique which has the potential to replace conventional phenotypic identification for most bacterial strains routinely isolated in clinical microbiology laboratories.

Death caused by cardioinhibitory reflex cardiac arrest--a systematic review of cases.

Forensic pathologists often refer to the cardioinhibitory reflex cardiac arrest (CiRCA) following short neck trauma as a mechanism of death. We sought via a systematic review of the literature to identify circumstances under which carotid bifurcation stimulation could lead to death. Two independent reviewers selected case studies or reports from Medline, ISI Web of Knowledge, and Embase. Circumstances and contributory factors were extracted for each case. From the available data, authors independently assessed whether CiRCA was highly probable (no alternative explanation possible), probable (alternative explanation possible), or unlikely (alternative explanation highly probable). A narrative approach was used to define circumstances in which CiRCA remained possible. From the 48 published cases evoking CiRCA as a possible cause of death between 1881 and 2009, 28 were most likely to result of other mechanism of death (i.e., cerebral hypoxia due to carotid compression, mechanical asphyxia, myocardial infarction). CiRCA remained possible for 20 cases (including five based on anecdotal evidence only) with only one case with no alternative explanation other than CiRCA. Our findings support the presumption that reflex cardiac arrhythmia due to carotid bifurcation stimulation cannot provoke death alone. Actual state of knowledge suggests CiRCA might be contributory to death in the presence of drug abuse and/or cardiac pathology, often associated with physical and/or mental excitation.

Safety and occlusion rates of surgical treatment of unruptured intracranial aneurysms: a systematic review and meta-analysis of the literature from 1990 to 2011.

BACKGROUND AND PURPOSE: Surgical clipping of unruptured intracranial aneurysms (UIAs) has recently been challenged by the emergence of endovascular treatment. We performed an updated systematic review and meta-analysis on the surgical treatment of UIAs, in an attempt to determine the aneurysm occlusion rates and safety of surgery in the modern era. METHODS: A detailed protocol was developed prior to conducting the review according to the Cochrane Collaboration guidelines. Electronic databases spanning January 1990-April 2011 were searched, complemented by hand searching. Heterogeneity was assessed using I(2), and publication bias with funnel plots. Surgical mortality and morbidity were analysed with weighted random effect models. RESULTS: 60 studies with 9845 patients harbouring 10 845 aneurysms were included. Mortality occurred in 157 patients (1.7%; 99% CI 0.9% to 3.0%; I(2)=82%). Unfavourable outcomes, including death, occurred in 692 patients (6.7%; 99% CI 4.9% to 9.0%; I(2)=85%). Morbidity rates were significantly greater in higher quality studies, and with large or posterior circulation aneurysms. Reported morbidity rates decreased over time. Studies were generally of poor quality; funnel plots showed heterogeneous results and publication bias, and data on aneurysm occlusion rates were scant. CONCLUSIONS: In studies published between 1990 and 2011, clipping of UIAs was associated with 1.7% mortality and 6.7% overall morbidity. The reputed durability of clipping has not been rigorously documented. Due to the quality of the included studies, the available literature cannot properly guide clinical decisions.

PFO closure vs. medical therapy in cryptogenic stroke or transient ischemic attack: A systematic review and meta-analysis.

BACKGROUND/OBJECTIVES: This study aims to assess whether patent foramen ovale (PFO) closure is superior to medical therapy in preventing recurrence of cryptogenic ischemic stroke or transient ischemic attack (TIA). METHODS: We searched PubMed for randomized trials which compared PFO closure with medical therapy in cryptogenic stroke/TIA using the items: "stroke or cerebrovascular accident or TIA" and "patent foramen ovale or paradoxical embolism" and "trial or study". RESULTS: Among 650 potentially eligible articles, 3 were included including 2303 patients. There was no statistically significant difference between PFO-closure and medical therapy in ischemic stroke recurrence (1.91% vs. 2.94% respectively, OR: 0.64, 95%CI: 0.37-1.10), TIA (2.08% vs. 2.42% respectively, OR: 0.87, 95%CI: 0.50-1.51) and death (0.60% vs. 0.86% respectively, OR: 0.71, 95%CI: 0.28-1.82). In subgroup analysis, there was significant reduction of ischemic strokes in the AMPLATZER PFO Occluder arm vs. medical therapy (1.4% vs. 3.04% respectively, OR: 0.46, 95%CI: 0.21-0.98, relative-risk-reduction: 53.2%, absolute-risk-reduction: 1.6%, number-needed-to-treat: 61.8) but not in the STARFlex device (2.7% vs. 2.8% with medical therapy, OR: 0.93, 95%CI: 0.45-2.11). Compared to medical therapy, the number of patients with new-onset atrial fibrillation (AF) was similar in the AMPLATZER PFO Occluder arm (0.72% vs. 1.28% respectively, OR: 1.81, 95%CI: 0.60-5.42) but higher in the STARFlex device (0.64% vs. 5.14% respectively, OR: 8.30, 95%CI: 2.47-27.84). CONCLUSIONS: This meta-analysis does not support PFO closure for secondary prevention with unselected devices in cryptogenic stroke/TIA. In subgroup analysis, selected closure devices may be superior to medical therapy without increasing the risk of new-onset AF, however. This observation should be confirmed in further trials using inclusion criteria for patients with high likelihood of PFO-related stroke recurrence.

Systematic identification of genomic markers of drug sensitivity in cancer cells.

Clinical responses to anticancer therapies are often restricted to a subset of patients. In some cases, mutated cancer genes are potent biomarkers for responses to targeted agents. Here, to uncover new biomarkers of sensitivity and resistance to cancer therapeutics, we screened a panel of several hundred cancer cell lines--which represent much of the tissue-type and genetic diversity of human cancers--with 130 drugs under clinical and preclinical investigation. In aggregate, we found that mutated cancer genes were associated with cellular response to most currently available cancer drugs. Classic oncogene addiction paradigms were modified by additional tissue-specific or expression biomarkers, and some frequently mutated genes were associated with sensitivity to a broad range of therapeutic agents. Unexpected relationships were revealed, including the marked sensitivity of Ewing's sarcoma cells harbouring the EWS (also known as EWSR1)-FLI1 gene translocation to poly(ADP-ribose) polymerase (PARP) inhibitors. By linking drug activity to the functional complexity of cancer genomes, systematic pharmacogenomic profiling in cancer cell lines provides a powerful biomarker discovery platform to guide rational cancer therapeutic strategies.

Alcohol Brief Intervention in Primary Care: A Systematic Review

BACKGROUND: The purpose of the present review was to evaluate the evidence of the effectiveness of brief interventions aimed at reducing chronic alcohol use and harm related to alcohol consumption, conducted among individuals actively attending primary care but who were not seeking help for alcohol problems. METHODS: Randomised trials reporting at-least one outcome related to alcohol consumption and conducted in outpatients who were actively attending primary care centre or provider were selected using Cochrane Central Register of Controlled Trials, MEDLINE, PsycINFO, ISI Web of Science, ETOH database, and bibliographies of the retrieved references and previous reviews. Selection and data abstraction were performed independently and in duplicate. We assessed validity of the studies and performed a meta-analysis for studies reporting alcohol consumption at 6 or 12 months follow up. RESULTS: We included 24 reports, reporting results of 19 trials and including 5,639 individuals. Seventeen trials reported a measure of alcohol consumption, eight reporting a significant effect of intervention. The meta-analysis showed a mean pooled difference of -41 (95% CI: −54; −28) g of pure ethanol per week in favour of brief intervention group. Evidences for other outcomes (laboratory values, health related quality of life, morbidity and mortality, health care utilisation) were inconclusive. CONCLUSION: Our systematic review indicated that brief intervention might be effective for both men and women in reducing alcohol consumption compared to a controlled intervention, in a primary health care population. The meta-analysis confirmed the reduction in alcohol consumption at 6 and 12 month. Further research should precise the components of effectiveness of brief intervention and the evidence of effects on morbidity, mortality, and quality of life related outcomes.

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.