Computational methods to create and analyze a digital gene expression atlas of embryo development from microscopy images

Abstract The creation of atlases, or digital models where information from different subjects can be combined, is a field of increasing interest in biomedical imaging. When a single image does not contain enough information to appropriately describe the organism under study, it is then necessary to acquire images of several individuals, each of them containing complementary data with respect to the rest of the components in the cohort. This approach allows creating digital prototypes, ranging from anatomical atlases of human patients and organs, obtained for instance from Magnetic Resonance Imaging, to gene expression cartographies of embryo development, typically achieved from Light Microscopy. Within such context, in this PhD Thesis we propose, develop and validate new dedicated image processing methodologies that, based on image registration techniques, bring information from multiple individuals into alignment within a single digital atlas model. We also elaborate a dedicated software visualization platform to explore the resulting wealth of multi-dimensional data and novel analysis algo-rithms to automatically mine the generated resource in search of bio¬logical insights. In particular, this work focuses on gene expression data from developing zebrafish embryos imaged at the cellular resolution level with Two-Photon Laser Scanning Microscopy. Disposing of quantitative measurements relating multiple gene expressions to cell position and their evolution in time is a fundamental prerequisite to understand embryogenesis multi-scale processes. However, the number of gene expressions that can be simultaneously stained in one acquisition is limited due to optical and labeling constraints. These limitations motivate the implementation of atlasing strategies that can recreate a virtual gene expression multiplex. The developed computational tools have been tested in two different scenarios. The first one is the early zebrafish embryogenesis where the resulting atlas constitutes a link between the phenotype and the genotype at the cellular level. The second one is the late zebrafish brain where the resulting atlas allows studies relating gene expression to brain regionalization and neurogenesis. The proposed computational frameworks have been adapted to the requirements of both scenarios, such as the integration of partial views of the embryo into a whole embryo model with cellular resolution or the registration of anatom¬ical traits with deformable transformation models non-dependent on any specific labeling. The software implementation of the atlas generation tool (Match-IT) and the visualization platform (Atlas-IT) together with the gene expression atlas resources developed in this Thesis are to be made freely available to the scientific community. Lastly, a novel proof-of-concept experiment integrates for the first time 3D gene expression atlas resources with cell lineages extracted from live embryos, opening up the door to correlate genetic and cellular spatio-temporal dynamics. La creación de atlas, o modelos digitales, donde la información de distintos sujetos puede ser combinada, es un campo de creciente interés en imagen biomédica. Cuando una sola imagen no contiene suficientes datos como para describir apropiadamente el organismo objeto de estudio, se hace necesario adquirir imágenes de varios individuos, cada una de las cuales contiene información complementaria respecto al resto de componentes del grupo. De este modo, es posible crear prototipos digitales, que pueden ir desde atlas anatómicos de órganos y pacientes humanos, adquiridos por ejemplo mediante Resonancia Magnética, hasta cartografías de la expresión genética del desarrollo de embrionario, típicamente adquiridas mediante Microscopía Optica. Dentro de este contexto, en esta Tesis Doctoral se introducen, desarrollan y validan nuevos métodos de procesado de imagen que, basándose en técnicas de registro de imagen, son capaces de alinear imágenes y datos provenientes de múltiples individuos en un solo atlas digital. Además, se ha elaborado una plataforma de visualization específicamente diseñada para explorar la gran cantidad de datos, caracterizados por su multi-dimensionalidad, que resulta de estos métodos. Asimismo, se han propuesto novedosos algoritmos de análisis y minería de datos que permiten inspeccionar automáticamente los atlas generados en busca de conclusiones biológicas significativas. En particular, este trabajo se centra en datos de expresión genética del desarrollo embrionario del pez cebra, adquiridos mediante Microscopía dos fotones con resolución celular. Disponer de medidas cuantitativas que relacionen estas expresiones genéticas con las posiciones celulares y su evolución en el tiempo es un prerrequisito fundamental para comprender los procesos multi-escala característicos de la morfogénesis. Sin embargo, el número de expresiones genéticos que pueden ser simultáneamente etiquetados en una sola adquisición es reducido debido a limitaciones tanto ópticas como del etiquetado. Estas limitaciones requieren la implementación de estrategias de creación de atlas que puedan recrear un multiplexado virtual de expresiones genéticas. Las herramientas computacionales desarrolladas han sido validadas en dos escenarios distintos. El primer escenario es el desarrollo embrionario temprano del pez cebra, donde el atlas resultante permite constituir un vínculo, a nivel celular, entre el fenotipo y el genotipo de este organismo modelo. El segundo escenario corresponde a estadios tardíos del desarrollo del cerebro del pez cebra, donde el atlas resultante permite relacionar expresiones genéticas con la regionalización del cerebro y la formación de neuronas. La plataforma computacional desarrollada ha sido adaptada a los requisitos y retos planteados en ambos escenarios, como la integración, a resolución celular, de vistas parciales dentro de un modelo consistente en un embrión completo, o el alineamiento entre estructuras de referencia anatómica equivalentes, logrado mediante el uso de modelos de transformación deformables que no requieren ningún marcador específico. Está previsto poner a disposición de la comunidad científica tanto la herramienta de generación de atlas (Match-IT), como su plataforma de visualización (Atlas-IT), así como las bases de datos de expresión genética creadas a partir de estas herramientas. Por último, dentro de la presente Tesis Doctoral, se ha incluido una prueba conceptual innovadora que permite integrar los mencionados atlas de expresión genética tridimensionales dentro del linaje celular extraído de una adquisición in vivo de un embrión. Esta prueba conceptual abre la puerta a la posibilidad de correlar, por primera vez, las dinámicas espacio-temporales de genes y células.

Comparison of dynamic effects of high-speed traffic load on ballasted track using a simplified two-dimensional and full three-dimensional model.

The vertical dynamic actions transmitted by railway vehicles to the ballasted track infrastructure is evaluated taking into account models with different degree of detail. In particular, we have studied this matter from a two-dimensional (2D) finite element model to a fully coupled three-dimensional (3D) multi-body finite element model. The vehicle and track are coupled via a non-linear Hertz contact mechanism. The method of Lagrange multipliers is used for the contact constraint enforcement between wheel and rail. Distributed elevation irregularities are generated based on power spectral density (PSD) distributions which are taken into account for the interaction. The numerical simulations are performed in the time domain, using a direct integration method for solving the transient problem due to the contact nonlinearities. The results obtained include contact forces, forces transmitted to the infrastructure (sleeper) by railpads and envelopes of relevant results for several track irregularities and speed ranges. The main contribution of this work is to identify and discuss coincidences and differences between discrete 2D models and continuum 3D models, as wheel as assessing the validity of evaluating the dynamic loading on the track with simplified 2D models

Molecular basis of salt tolerance in Physcomitrella Patens model plant: potassium homeostasis and physiological roles of chx transporters

El suelo salino impone un estrés abiótico importante que causa graves problemas en la agricultura ya que la mayoría de los cultivos se ven afectados por la salinidad debido a efectos osmóticos y tóxicos. Por ello, la contaminación y la escasez de agua dulce, la salinización progresiva de tierras y el aumento exponencial de la población humana representan un grave problema que amenaza la seguridad alimentaria mundial para las generaciones futuras. Por lo tanto, aumentar la tolerancia a la salinidad de los cultivos es un objetivo estratégico e ineludible para garantizar el suministro de alimentos en el futuro. Mantener una óptima homeostasis de K+ en plantas que sufren estrés salino es un objetivo importante en el proceso de obtención de plantas tolerantes a la salinidad. Aunque el modelo de la homeostasis de K+ en las plantas está razonablemente bien descrito en términos de entrada de K+, muy poco se sabe acerca de los genes implicados en la salida de K+ o de su liberación desde la vacuola. En este trabajo se pretende aclarar algunos de los mecanismos implicados en la homeostasis de K+ en plantas. Para ello se eligió la briofita Physcomitrella patens, una planta no vascular de estructura simple y de fase haploide dominante que, entre muchas otras cualidades, hacen que sea un modelo ideal. Lo más importante es que no sólo P. patens es muy tolerante a altas concentraciones de Na+, sino que también su posición filogenética en la evolución de las plantas abre la posibilidad de estudiar los cambios claves que, durante el curso de la evolución, se produjeron en las diversas familias de los transportadores de K+. Se han propuesto varios transportadores de cationes como candidatos que podrían tener un papel en la salida de K+ o su liberación desde la vacuola, especialmente miembros de la familia CPA2 que contienen las familias de transportadores KEA y CHX. En este estudio se intenta aumentar nuestra comprensión de las funciones de los transportadores de CHX en las células de las plantas usando P. patens, como ya se ha dicho. En esta especie, se han identificado cuatro genes CHX, PpCHX1-4. Dos de estos genes, PpCHX1 y PpCHX2, se expresan aproximadamente al mismo nivel que el gen PpACT5, y los otros dos genes muestran una expresión muy baja. La expresión de PpCHX1 y PpCHX2 en mutantes de Escherichia coli defectivos en el transporte de K+ restauraron el crecimiento de esta cepa en medios con bajo contenido de K+, lo que viii sugiere que la entrada de K+ es energizada por un mecanismo de simporte con H+. Por otra parte, estos transportadores suprimieron el defecto asociado a la mutación kha1 en Saccharomyces cerevisiae, lo que sugiere que podrían mediar un antiporte en K+/H+. La proteína PpCHX1-GFP expresada transitoriamente en protoplastos de P. patens co-localizó con un marcador de Golgi. En experimentos similares, la proteína PpCHX2-GFP localizó aparentemente en la membrana plasmática y tonoplasto. Se construyeron las líneas mutantes simples de P. patens ΔPpchx1 y ΔPpchx2, y también el mutante doble ΔPpchx2 ΔPphak1. Los mutantes simples crecieron normalmente en todas las condiciones ensayadas y mostraron flujos de entrada normales de K+ y Rb+; la mutación ΔPpchx2 no aumentó el defecto de las plantas ΔPphak1. En experimentos a largo plazo, las plantas ΔPpchx2 mostraron una retención de Rb+ ligeramente superior que las plantas silvestres, lo que sugiere que PpCHX2 promueve la transferencia de Rb+ desde la vacuola al citosol o desde el citosol al medio externo, actuando en paralelo con otros transportadores. Sugerimos que transportadores de K+ de varias familias están involucrados en la homeostasis de pH de orgánulos ya sea mediante antiporte K+/H+ o simporte K+-H+.ix ABSTRACT Soil salinity is a major abiotic stress causing serious problems in agriculture as most crops are affected by it. Moreover, the contamination and shortage of freshwater, progressive land salinization and exponential increase of human population aggravates the problem implying that world food security may not be ensured for the next generations. Thus, a strategic and an unavoidable goal would be increasing salinity tolerance of plant crops to secure future food supply. Maintaining an optimum K+ homeostasis in plants under salinity stress is an important trait to pursue in the process of engineering salt tolerant plants. Although the model of K+ homeostasis in plants is reasonably well described in terms of K+ influx, very little is known about the genes implicated in K+ efflux or release from the vacuole. In this work, we aim to clarify some of the mechanisms involved in K+ homeostasis in plants. For this purpose, we chose the bryophyte plant Physcomitrella patens, a nonvascular plant of simple structure and dominant haploid phase that, among many other characteristics, makes it an ideal model. Most importantly, not only P. patens is very tolerant to high concentrations of Na+, but also its phylogenetic position in land plant evolution opens the possibility to study the key changes that occurred in K+ transporter families during the course of evolution. Several cation transporter candidates have been proposed to have a role in K+ efflux or release from the vacuole especially members of the CPA2 family which contains the KEA and CHX transporter families. We intended in this study to increase our understanding of the functions of CHX transporters in plant cells using P. patens, in which four CHX genes have been identified, PpCHX1-4. Two of these genes, PpCHX1 and PpCHX2, are expressed at approximately the same level as the PpACT5 gene, but the other two genes show an extremely low expression. PpCHX1 and PpCHX2 restored growth of Escherichia coli mutants on low K+-containing media, suggesting they mediated K+ uptake that may be energized by symport with H+. In contrast, these genes suppressed the defect associated to the kha1 mutation in Saccharomyces cerevisiae, which suggest that they might mediate K+/H+ antiport. PpCHX1-GFP protein transiently expressed in P. patens protoplasts co-localized with a Golgi marker. In similar experiments, the PpCHX2-GFP protein appeared to localize to tonoplast and plasma x membrane. We constructed the ΔPpchx1 and ΔPpchx2 single mutant lines, and the ΔPpchx2 ΔPphak1 double mutant. Single mutant plants grew normally under all the conditions tested and exhibited normal K+ and Rb+ influxes; the ΔPpchx2 mutation did not increase the defect of ΔPphak1 plants. In long-term experiments, ΔPpchx2 plants showed a slightly higher Rb+ retention than wild type plants, which suggests that PpCHX2 mediates the transfer of Rb+ from either the vacuole to the cytosol or from the cytosol to the external medium in parallel with other transporters. We suggest that K+ transporters of several families are involved in the pH homeostasis of organelles by mediating either K+/H+ antiport or K+-H+ symport.

Competitive exclusion in a two-especies chemotasis model

We consider a mathematical model for the spatio-temporal evolution of two biological species in a competitive situation. Besides diffusing, both species move toward higher concentrations of a chemical substance which is produced by themselves. The resulting system consists of two parabolic equations with Lotka–Volterra-type kinetic terms and chemotactic cross-diffusion, along with an elliptic equation describing the behavior of the chemical. We study the question in how far the phenomenon of competitive exclusion occurs in such a context. We identify parameter regimes for which indeed one of the species dies out asymptotically, whereas the other reaches its carrying capacity in the large time limit.

On a two species chemotaxis model with slow chemical diffucion.

In this paper we consider a system of three parabolic equations modeling the behavior of two biological species moving attracted by a chemical factor. The chemical substance verifies a parabolic equation with slow diffusion. The system contains second order terms in the first two equations modeling the chemotactic effects. We apply an iterative method to obtain the global existence of solutions using that the total mass of the biological species is conserved. The stability of the homogeneous steady states is studied by using an energy method. A final example is presented to illustrate the theoretical results.

Hepatocyte gene therapy in a large animal: A neonatal bovine model of citrullinemia

The development of gene-replacement therapy for inborn errors of metabolism has been hindered by the limited number of suitable large-animal models of these diseases and by inadequate methods of assessing the efficacy of treatment. Such methods should provide sensitive detection of expression in vivo and should be unaffected by concurrent pharmacologic and dietary regimens. We present the results of studies in a neonatal bovine model of citrullinemia, an inborn error of urea-cycle metabolism characterized by deficiency of argininosuccinate synthetase and consequent life-threatening hyperammonemia. Measurements of the flux of nitrogen from orally administered 15NH4 to [15N]urea were used to determine urea-cycle activity in vivo. In control animals, these isotopic measurements proved to be unaffected by pharmacologic treatments. Systemic administration of a first-generation E1-deleted adenoviral vector expressing human argininosuccinate synthetase resulted in transduction of hepatocytes and partial correction of the enzyme defect. The isotopic method showed significant restoration of urea synthesis. Moreover, the calves showed clinical improvement and normalization of plasma glutamine levels after treatment. The results show the clinical efficacy of treating a large-animal model of an inborn error of hepatocyte metabolism in conjunction with a method for sensitively measuring correction in vivo. These studies will be applicable to human trials of the treatment of this disorder and other related urea-cycle disorders.

A third member of the synapsin gene family

Synapsins are a family of neuron-specific synaptic vesicle-associated phosphoproteins that have been implicated in synaptogenesis and in the modulation of neurotransmitter release. In mammals, distinct genes for synapsins I and II have been identified, each of which gives rise to two alternatively spliced isoforms. We have now cloned and characterized a third member of the synapsin gene family, synapsin III, from human DNA. Synapsin III gives rise to at least one protein isoform, designated synapsin IIIa, in several mammalian species. Synapsin IIIa is associated with synaptic vesicles, and its expression appears to be neuron-specific. The primary structure of synapsin IIIa conforms to the domain model previously described for the synapsin family, with domains A, C, and E exhibiting the highest degree of conservation. Synapsin IIIa contains a novel domain, termed domain J, located between domains C and E. The similarities among synapsins I, II, and III in domain organization, neuron-specific expression, and subcellular localization suggest a possible role for synapsin III in the regulation of neurotransmitter release and synaptogenesis. The human synapsin III gene is located on chromosome 22q12–13, which has been identified as a possible schizophrenia susceptibility locus. On the basis of this localization and the well established neurobiological roles of the synapsins, synapsin III represents a candidate gene for schizophrenia.

Substitution of a mutant α2a-adrenergic receptor via “hit and run” gene targeting reveals the role of this subtype in sedative, analgesic, and anesthetic-sparing responses in vivo

Norepinephrine contributes to antinociceptive, sedative, and sympatholytic responses in vivo, and α2 adrenergic receptor (α2AR) agonists are used clinically to mimic these effects. Lack of subtype-specific agonists has prevented elucidation of the role that each α2AR subtype (α2A, α2B, and α2C) plays in these central effects. Here we demonstrate that α2AR agonist-elicited sedative, anesthetic-sparing, and analgesic responses are lost in a mouse line expressing a subtly mutated α2AAR, D79N α2AAR, created by two-step homologous recombination. These functional changes are accompanied by failure of the D79N α2AAR to inhibit voltage-gated Ca2+ currents and spontaneous neuronal firing, a measure of K+ current activation. These results provide definitive evidence that the α2AAR subtype is the primary mediator of clinically important central actions of α2AR agonists and suggest that the D79N α2AAR mouse may serve as a model for exploring other possible α2AAR functions in vivo.

Phylogenomic analysis of natural products biosynthetic gene clusters allows discovery of arseno-organic metabolites in model streptomycetes

We are indebted with Marnix Medema, Paul Straight and Sean Rovito, for useful discussions and critical reading of the manuscript, as well as with Alicia Chagolla and Yolanda Rodriguez of the MS Service of Unidad Irapuato, Cinvestav, and Araceli Fernandez for technical support in high-performance computing. This work was funded by Conacyt Mexico (grants No. 179290 and 177568) and FINNOVA Mexico (grant No. 214716) to FBG. PCM was funded by Conacyt scholarship (No. 28830) and a Cinvestav posdoctoral fellowship. JF and JFK acknowledge funding from the College of Physical Sciences, University of Aberdeen, UK.

Two sterol regulatory element-like sequences mediate up-regulation of caveolin gene transcription in response to low density lipoprotein free cholesterol

Caveolae form the terminus for a major pathway of intracellular free cholesterol (FC) transport. Caveolin mRNA levels in confluent human skin fibroblasts were up-regulated following increased uptake of low density lipoprotein (LDL) FC. The increase induced by FC was not associated with detectable change in mRNA stability, indicating that caveolin mRNA levels were mediated at the level of gene transcription. A total of 924 bp of 5′ flanking region of the caveolin gene were cloned and sequenced. The promoter sequence included three G+C-rich potential sterol regulatory elements (SREs), a CAAT sequence and a Sp1 consensus sequence. Deletional mutagenesis of individual SRE-like sequences indicated that of these two (at −646 and −395 bp) were essential for the increased transcription rates mediated by LDL-FC, whereas the third was inconsequential. Gel shift analysis of protein binding from nuclear extracts to these caveolin promoter DNA sequences, together with DNase I footprinting, confirmed nucleoprotein binding to the SRE-like elements as part of the transcriptional response to LDL-FC. A supershift obtained with antibody to SRE-binding protein 1 (SPEBP-1) indicated that this protein binds at −395 bp. There was no reaction at −395 bp with anti-Sp1 antibody nor with either antibody at −646 bp. The cysteine protease inhibitor N-acetyl-leu-leu-norleucinal (ALLN), which inhibits SREBP catabolism, superinhibited caveolin mRNA levels regardless of LDL-FC. This finding suggests that SREBP inhibits caveolin gene transcription in contrast to its stimulating effect on other promoters. The findings of this study are consistent with the postulated role for caveolin as a regulator of cellular FC homeostasis in quiescent peripheral cells, and the coordinate regulation by SREBP of FC influx and efflux.

The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

A factor IX-deficient mouse model for hemophilia B gene therapy

We have generated a mouse where the clotting factor IX (FIX) gene has been disrupted by homologous recombination. The FIX nullizygous (−/−) mouse was devoid of factor IX antigen in plasma. Consistent with the bleeding disorder, the factor IX coagulant activities for wild-type (+/+), heterozygous (+/−), and homozygous (−/−) mice were 92%, 53%, and <5%, respectively, in activated partial thromboplastin time assays. Plasma factor IX activity in the deficient mice (−/−) was restored by introducing wild-type murine FIX gene via adenoviral vectors. Thus, these factor IX-deficient mice provide a useful animal model for gene therapy studies of hemophilia B.

Gene replacement by homologous recombination in the multicellular green alga Volvox carteri

With only two different cell types, the haploid green alga Volvox represents the simplest multicellular model system. To facilitate genetic investigations in this organism, the occurrence of homologous recombination events was investigated with the intent of developing methods for gene replacement and gene disruption. First, homologous recombination between two plasmids was demonstrated by using overlapping nonfunctional fragments of a recombinant arylsulfatase gene (tubulin promoter/arylsulfatase gene). After bombardment of Volvox reproductive cells with DNA-coated gold microprojectiles, transformants expressing arylsulfatase constitutively were recovered, indicating the presence of the machinery for homologous recombination in Volvox. Second, a well characterized loss-of-function mutation in the nuclear nitrate reductase gene (nitA) with a single G → A nucleotide exchange in a 5′-splice site was chosen as a target for gene replacement. Gene replacement by homologous recombination was observed with a reasonably high frequency only if the replacement vector containing parts of the functional nitrate reductase gene contained only a few nucleotide exchanges. The ratio of homologous to random integration events ranged between 1:10 and 1:50, i.e., homologous recombination occurs frequently enough in Volvox to apply the powerful tool of gene disruption for functional studies of novel genes.

Synergistic inhibition of tumor growth in a murine mammary adenocarcinoma model by combinational gene therapy using IL-12, pro-IL-18, and IL-1β converting enzyme cDNA

We report here that a cancer gene therapy protocol using a combination of IL-12, pro-IL-18, and IL-1β converting enzyme (ICE) cDNA expression vectors simultaneously delivered via gene gun can significantly augment antitumor effects, evidently by generating increased levels of bioactive IL-18 and consequently IFN-γ. First, we compared the levels of IFN-γ secreted by mouse splenocytes stimulated with tumor cells transfected with various test genes, including IL-12 alone; pro-IL-18 alone; pro-IL-18 and ICE; IL-12 and pro-IL-18; and IL-12, pro-IL-18, and ICE. Among these treatments, the combination of IL-12, pro-IL-18, and ICE cDNA resulted in the highest level of IFN-γ production from splenocytes in vitro, and similar results were obtained when these same treatments were delivered to the skin of a mouse by gene gun and IFN-γ levels were measured at the skin transfection site in vivo. Furthermore, the triple gene combinatorial gene therapy protocol was the most effective among all tested groups at suppressing the growth of TS/A (murine mammary adenocarcinoma) tumors previously implanted intradermally at the skin site receiving DNA transfer by gene gun on days 6, 8, 10, and 12 after tumor implantation. Fifty percent of mice treated with the combined three-gene protocol underwent complete tumor regression. In vivo depletion experiments showed that this antitumor effect was CD8+ T cell-mediated and partially IFN-γ-dependent. These results suggest that a combinatorial gene therapy protocol using a mixture of IL-12, pro-IL-18, and ICE cDNAs can confer potent antitumor activities against established TS/A tumors via cytotoxic CD8+ T cells and IFN-γ-dependent pathways.