Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Monophyletic origin of multiple clonal lineages in an asexual fish (Poecilia formosa).

Despite the advantage of avoiding the costs of sexual reproduction, asexual vertebrates are very rare and often considered evolutionarily disadvantaged when compared to sexual species. Asexual species, however, may have advantages when colonizing (new) habitats or competing with sexual counterparts. They are also evolutionary older than expected, leaving the question whether asexual vertebrates are not only rare because of their 'inferior' mode of reproduction but also because of other reasons. A paradigmatic model system is the unisexual Amazon molly, Poecilia formosa, that arose by hybridization of the Atlantic molly, Poecilia mexicana, as the maternal ancestor, and the sailfin molly, Poecilia latipinna, as the paternal ancestor. Our extensive crossing experiments failed to resynthesize asexually reproducing (gynogenetic) hybrids confirming results of previous studies. However, by producing diploid eggs, female F(1) -hybrids showed apparent preadaptation to gynogenesis. In a range-wide analysis of mitochondrial sequences, we examined the origin of P. formosa. Our analyses point to very few or even a single origin(s) of its lineage, which is estimated to be approximately 120,000 years old. A monophyletic origin was supported from nuclear microsatellite data. Furthermore, a considerable degree of genetic variation, apparent by high levels of clonal microsatellite diversity, was found. Our molecular phylogenetic evidence and the failure to resynthesize the gynogenetic P. formosa together with the old age of the species indicate that some unisexual vertebrates might be rare not because they suffer the long-term consequences of clonal reproduction but because they are only very rarely formed as a result of complex genetic preconditions necessary to produce viable and fertile clonal genomes and phenotypes ('rare formation hypothesis').

A gene-rich, transcriptionally active environment and the pre-deposition of repressive marks are predictive of susceptibility to KRAB/KAP1-mediated silencing.

AbstractBACKGROUND: KRAB-ZFPs (Krüppel-associated box domain-zinc finger proteins) are vertebrate-restricted transcriptional repressors encoded in the hundreds by the mouse and human genomes. They act via an essential cofactor, KAP1, which recruits effectors responsible for the formation of facultative heterochromatin. We have recently shown that KRAB/KAP1 can mediate long-range transcriptional repression through heterochromatin spreading, but also demonstrated that this process is at times countered by endogenous influences.METHOD: To investigate this issue further we used an ectopic KRAB-based repressor. This system allowed us to tether KRAB/KAP1 to hundreds of euchromatic sites within genes, and to record its impact on gene expression. We then correlated this KRAB/KAP1-mediated transcriptional effect to pre-existing genomic and chromatin structures to identify specific characteristics making a gene susceptible to repression.RESULTS: We found that genes that were susceptible to KRAB/KAP1-mediated silencing carried higher levels of repressive histone marks both at the promoter and over the transcribed region than genes that were insensitive. In parallel, we found a high enrichment in euchromatic marks within both the close and more distant environment of these genes.CONCLUSION: Together, these data indicate that high levels of gene activity in the genomic environment and the pre-deposition of repressive histone marks within a gene increase its susceptibility to KRAB/KAP1-mediated repression.

DNA polymorphism of schistosomes and their snail hosts

Analysis of the genomes of schistosomes and one of their intermediate hosts, Biomphalaria glabrata, using Random Amplified Polymorphic DNA (RAPD) demonstrated that intraspecific genetic polymorphism in the parasite is limited but in the snail is highly pronounced. This suggests an important role for the snail in the determination of the epidemiology of the disease. In addition to their intraspecific stability, schistosome derived RAPDs exhibit a high level of interspecific polymorphism and are thus ideal for the construction of phylogenetic trees. For the detection of intraspecific polymorphisms extensive variation in the mitochondrial DNA is being exploited for the development of a PCR based test for Schistosoma mansoni. Gene level polymorphisms are being analyzed by Low Stringency Single Specific Primer PCR.

Desarrollo de una colección de Líneas Casi Isogénicas (NIL) para el estudio de caracteres de interés agronómico en especies del género Fragaria

Fragaria vesca posee varias características que la hacen interesante como especie modelo en la familia Rosaceae. Además, su naturaleza diploide permite sortear la complejidad genética de la fresa cultivada. Considerando que los genomas de las especies diploides y octoploides de Fragaria mantienen una relación de sintenia y colinearidad, estamos desarrollando una colección de Líneas Casi Isogénicas (NIL) en Fragaria diploide, usando Fragaria bucharica, una variedad asiática, como donante de introgresiones y Fragaria vesca var. Reine des Vallées, una variedad francesa comúnmente cultivada en España para usos industriales, como parental recurrente. Para obtener introgresiones de F. bucharica en un fondo genético homogéneo de F. vesca, se desarrolló un retrocruzamiento y se obtuvo una población BC1 que fue analizada para la presencia de alelos de F. bucharica a lo largo de los 7 grupos de ligamiento (LG) del mapa de referencia de Fragaria. Los individuos con bajo número de introgresiones de gran tamaño, que en conjunto abarcaron todo el genoma de F. bucharica, fueron seleccionados y retrocruzados con el parental recurrente. La descendencia fue analizada para los loci segregantes y los individuos BC2 con 1 ó 2 introgresiones fueron seleccionados como líneas donadoras de introgresiones de las NIL. Hasta el momento se han descrito tres fenotipos dominantes diferentes bajo condiciones de invernadero en varias NIL heterozigóticas: Las plantas con introgresiones en LG2 presentan estolones. Las plantas con introgresiones en LG1 y LG3 presentan floración temprana. Las plantas con introgresiones en LG6 presentan floración estacional. Actualmente se está trabajando en la selección y caracterización de introgresiones de pequeño tamaño mediante la autofecundación de las líneas seleccionadas. La caracterización fenotípica de los recombinantes seleccionados permitirá localizar y estimar los patrones de herencia de los genes implicados en los caracteres descritos, además de la identificación de nuevos caracteres recesivos que aparecerán en homozigosis.

Origen de la multicel·lularitat : una aproximació genòmica i funcional

En el present estudi s'analitza l'origen i evolució de 2 molècules claus pera entendre la multicel·lularitat dels animals: les molècules d'adhesió integrines i els factors de transcripció T-box. S’utilitzen els genomes recentment publicats de protists unicel•lulars parents propers dels animals. S’analitza l’origen i evolució d’aquests gens mitjançant anàlisi filogènic, determinació de motius funcionals i també tècniques de biologia molecular. A més, es documenta un cas de transferència gènica horitzontal des d'un eucariota cap a un procariota, fenomen poc habitual. Les principals conclusions són que tant l’adhesoma d'integrina com els gens T-box tenen un origen molt anterior als animals, en un context unicel•lular, i que després foren cooptats pel llinatge multicel•lular dels animals.

Bioinformática: identificar genes en una interfaz gráfica vía web para la comparación de genomas

Este trabajo desarrolla el proceso de diseño e implementación de una interfaz web que permite la exploración en detalle de las relaciones entre genomas completos. La interfaz permite la comparación simultánea de nueve genomas, representando en cada gráfica las relaciones entre cada par de genomas junto los genes identificados de cada uno de ellos. Es capaz de trabajar con genomas del dominio Eukaryota y se adapta a la capacidad de cómputo de la máquina cliente. La información representada son MUMs (Maximal Unique Matching, secuencia máxima y única encontrada en ambos genomas) y SuperMUMs (agrupación de MUMs mediante Approximate String Matching). Los datos son previamente calculados y accesibles desde un servidor web.

Analysis of Respiratory Syncytial Virus in Clinical Samples by Reverse Transcriptase-Polymerase Chain Reaction Restriction Mapping

The aim of this study was to develop a polymerase chain reaction (PCR) for the detection of respiratory syncytial virus (RSV) genomes. The primers were designed from published sequences and selected from conserved regions of the genome encoding for the N protein of subgroups A and B of RSV. PCR was applied to 20 specimens from children admitted to the respiratory ward of "William Soler" Pediatric Hospital in Havana City with a clinical diagnosis of bronchiolitis. The PCR was compared with viral isolation and with an indirect immunofluorescence technique that employs monoclonal antibodies of subgroups A and B. Of 20 nasopharyngeal exudates, 10 were found positive by the three assayed methods. In only two cases, samples that yielded positive RNA-PCR were found negative by indirect immunofluorescence and cell culture. Considering viral isolation as the "gold standard" technique, RNA-PCR had 100% sensitivity and 80% specificity. RNA-PCR is a specific and sensitive technique for the detection of the RSV genome. Technical advantages are discussed

An activation-independent role of transcription factors in insulator function.

Chromatin insulators are defined as transcriptionally neutral elements that prevent negative or positive influence from extending across chromatin to a promoter. Here we show that yeast subtelomeric anti-silencing regions behave as boundaries to telomere-driven silencing and also allow discontinuous propagation of silent chromatin. These two facets of insulator activity, boundary and silencing discontinuity, can be recapitulated by tethering various transcription activation domains to tandem sites on DNA. Importantly, we show that these insulator activities do not involve direct transcriptional activation of the reporter promoter. These findings predict that certain promoters behave as insulators and partition genomes in functionally independent domains.

Parasite Genome Projects and the Trypanosoma cruzi Genome Initiative

Since the start of the human genome project, a great number of genome projects on other "model" organism have been initiated, some of them already completed. Several initiatives have also been started on parasite genomes, mainly through support from WHO/TDR, involving North-South and South-South collaborations, and great hopes are vested in that these initiatives will lead to new tools for disease control and prevention, as well as to the establishment of genomic research technology in developing countries. The Trypanosoma cruzi genome project, using the clone CL-Brener as starting point, has made considerable progress through the concerted action of more than 20 laboratories, most of them in the South. A brief overview of the current state of the project is given

arrayMap 2014: an updated cancer genome resource.

Somatic copy number aberrations (CNA) represent a mutation type encountered in the majority of cancer genomes. Here, we present the 2014 edition of arrayMap (http://www.arraymap.org), a publicly accessible collection of pre-processed oncogenomic array data sets and CNA profiles, representing a vast range of human malignancies. Since the initial release, we have enhanced this resource both in content and especially with regard to data mining support. The 2014 release of arrayMap contains more than 64,000 genomic array data sets, representing about 250 tumor diagnoses. Data sets included in arrayMap have been assembled from public repositories as well as additional resources, and integrated by applying custom processing pipelines. Online tools have been upgraded for a more flexible array data visualization, including options for processing user provided, non-public data sets. Data integration has been improved by mapping to multiple editions of the human reference genome, with the majority of the data now being available for the UCSC hg18 as well as GRCh37 versions. The large amount of tumor CNA data in arrayMap can be freely downloaded by users to promote data mining projects, and to explore special events such as chromothripsis-like genome patterns.

Genomics of "chlamydia"-related bacteria

AbstractThe Chlamydiales order is an important bacterial phylum that comprises some of the most successful human pathogens such as Chlamydia trachomatis, the leading infectious cause of blindness worldwide. Since some years, several new bacteria related to Chlamydia have been discovered in clinical or environmental samples and might represent emerging pathogens. The genome sequencing of classical Chlamydia has brought invaluable information on these obligate intracellular bacteria otherwise difficult to study due to the lack of tools to perform basic genetic manipulation. The recent emergence of high-throughput sequencing technologies yielding millions of reads in a short time lowered the costs of genome sequencing and thus represented a unique opportunity to study Chlamydia-re\ated bacteria. Based on the sequencing and the analysis of Chlamydiales genomes, this thesis provides significant insights into the genetic determinants of the intracellular lifestyle, the pathogenicity, the metabolism and the evolution of Chlamydia-related bacteria. A first approach showed the efficacy of rapid sequencing coupled to proteomics to identify immunogenic proteins. This method, particularly useful for an emerging pathogen such as Parachlamydia acanthamoebae, enabled us to discover good candidates for the development of diagnostic tools that would permit to evaluate at larger scale the role of this bacterium in disease. Second, the complete genome of Waddlia chondrophila, a potential agent of miscarriage, encodes numerous virulence factors to manipulate its host cell and resist to environmental stresses. The reconstruction of metabolic pathways showed that the bacterium possesses extensive capabilities compared to related organisms. However, it is still incapable of synthesizing some essential components and thus has to import them from its host. Third, the genome comparison of Protochlamydia naegleriophila to its closest known relative Protochlamydia amoebophila revealed a particular evolutionary dynamic with the occurrence of an unexpected genome rearrangement. Fourth, a phylogenetic analysis of P. acanthamoebae and Legionella drancourtii identified several genes probably exchanged by horizontal gene transfer with other intracellular bacteria that might occur within their amoebal host. These genes often encode mechanisms for resistance to metal or toxic compounds. As a whole, the analysis of the different genomes enabled us to highlight a large diversity in size, GC percentage, repeat content as well as plasmid organization. The abundant genomic data obtained during this thesis have a wide impact since they provide the necessary bases for detailed investigations on countless aspects of the biology and the evolution of Chlamydia-related bacteria, whether in wet lab or by bioinformatical analyses.RésuméL'ordre des Chlamydiales est un important phylum bactérien qui comprend de nombreuses espèces pathogènes pour l'homme et les animaux, dont Chlamydia trachomatis, responsable du trachome, la cause majeure de cécité d'origine infectieuse à travers le monde. Durant ces dernières décennies, de nombreuses bactéries apparentées aux Chlamydia ont été découvertes dans des échantillons environnementaux ou cliniques mais leur éventuel rôle pathogène dans le développement de maladies reste peu connu. Ces bactéries sont des intracellulaires obligatoires car elles ont besoin d'une cellule hôte pour se multiplier, ce qui rend leur étude particulièrement difficile. Le développement de nouvelles technologies permettant de séquencer le génome d'un organisme rapidement et à moindre coût ainsi que l'essor des méthodes d'analyse s'y rapportant représentent une opportunité exceptionnelle d'étudier ces organismes. Dans ce contexte, cette thèse démontre l'utilité de la génomique pour développer de nouveaux outils diagnostiques ainsi que pour étudier le métabolisme de ces bactéries, leurs facteurs de virulence et leur évolution.Ainsi, une première approche a illustré l'utilité d'un séquençage rapide pour obtenir les informations nécessaires à l'identification de protéines qui sont reconnues par des anticorps humains ou animaux. Cette méthode, particulièrement utile pour un pathogène émergent tel que Parachlamydia acanthamoebae, a permis de découvrir de bons candidats pour le développement d'un outil diagnostique qui permettrait d'évaluer à plus large échelle le rôle de cette bactérie notamment dans la pneumonie. L'analyse du contenu génique de Waddlia chondrophila, un autre germe qui pourrait être impliqué dans les avortements et tes fausses-couches, a en outre mis en évidence la présence de nombreux facteurs connus qui lui permettent de manipuler son hôte. Cette bactérie possède de plus grandes capacités métaboliques que les autres Chlamydia, mais elle est incapable de synthétiser certains composants et doit donc les importer de son hôte pour subvenir à ses besoins. La comparaison du génome de Protochlamydia naegleriophila à son plus proche parent, Protochlamydia amoebophila, a dévoilé une évolution dynamique particulière avec l'occurrence d'un réarrangement majeur inattendu après la séparation de ces deux espèces. En outre, ces études ont montré l'occurrence de plusieurs transferts de gène avec d'autres organismes plus éloignés, notamment d'autres intracellulaires d'amibes, souvent pour l'acquisition de mécanismes de résistances à des composés toxiques. Les données génomiques acquises durant ce travail posent les fondements nécessaires a de nombreuses analyses qui permettront progressivement de mieux comprendre de nombreux aspects de ces bactéries fascinantes.

Short term evolution of a highly transmissible methicillin-resistant Staphylococcus aureus clone (ST228) in a tertiary care hospital.

Staphylococcus aureus is recognized as one of the major human pathogens and is by far one of the most common nosocomial organisms. The genetic basis for the emergence of highly epidemic strains remains mysterious. Studying the microevolution of the different clones of S. aureus is essential for identifying the forces driving pathogen emergence and spread. The aim of the present study was to determine the genetic changes characterizing a lineage belonging to the South German clone (ST228) that spread over ten years in a tertiary care hospital in Switzerland. For this reason, we compared the whole genome of eight isolates recovered between 2001 and 2008 at the Lausanne hospital. The genetic comparison of these isolates revealed that their genomes are extremely closely related. Yet, a few more important genetic changes, such as the replacement of a plasmid, the loss of large fragments of DNA, or the insertion of transposases, were observed. These transfers of mobile genetic elements shaped the evolution of the ST228 lineage that spread within the Lausanne hospital. Nevertheless, although the strains analyzed differed in their dynamics, we have not been able to link a particular genetic element with spreading success. Finally, the present study showed that new sequencing technologies improve considerably the quality and quantity of information obtained for a single strain; but this information is still difficult to interpret and important investments are required for the technology to become accessible for routine investigations.

Segmental copy number variation shapes tissue transcriptomes

Abstract : Copy number variation (CNV) of DNA segments has recently gained considerable interest as a source of genetic variation likely to play a role in phenotypic diversity and evolution. Much effort has been put into the identification and mapping of regions that vary in copy number among seemingly normal individuals, both in humans and in a number of model organisms, using both bioinformatic and hybridization-based methods. Synteny studies suggest the existence of CNV hotspots in mammalian genomes, often in connection with regions of segmental duplication. CNV alleles can be in equilibrium within a population, but can also arise de novo between generations, illustrating the highly dynamic nature of these regions. A small number of studies have assessed the effect of CNV on single loci, however, at the genome-wide scale, the functional impact of CNV remains poorly studied. We have explored the influence of CNV on gene expression, first using the Williams-Beuren syndrome (WBS) associated deletion as a model, and second at the genome-wide scale in inbred mouse strains. We found that the WBS deletion influences the expression levels not only of the hemizygous genes, but also affects the euploid genes mapping nearby. Consistently, on a genome wide scale we observe that CNV genes are expressed at more variable levels than genes that do not vary in copy number. Likewise, CNVs influence the relative expression levels of genes that map to the flank of the genome rearrangements, thus globally influencing tissue transcriptomes. Further studies are warranted to complete cataloguing and fine mapping of CNV regions, as well as to elucidate the different mechanisms by which CNVs influence gene expression. Résumé : La variation en nombre de copies (copy number variation ou CNV) de segments d'ADN suscite un intérêt en tant que variation génétique susceptible de jouer un r81e dans la diversité phénotypique et l'évolution. Les régions variables en nombre de copies parmi des individus apparemment normaux ont été cartographiées et cataloguées au moyen de puces à ADN et d'analyse bioinformatique. L'étude de la synténie entre plusieurs espèces de mammifères laisse supposer l'existence de régions à haut taux de variation, souvent liées à des duplications segmentaires. Les allèles CNV peuvent être en équilibre au sein d'une population ou peuvent apparaître de novo. Ces faits illustrent la nature hautement dynamique de ces régions. Quelques études se sont penchées sur l'effet de la variation en nombre de copies de loci isolés, cependant l'impact de ce phénomène n'a pas été étudié à l'échelle génomique. Nous avons examiné l'influence des CNV sur l'expression des gènes. Dans un premier temps nous avons utilisé la délétion associée au syndrome de Williams-Beuren (WBS), puis, dans un second temps, nous avons poursuivi notre étude à l'échelle du génome, dans des lignées consanguines de souris. Nous avons établi que la délétion WBS influence l'expression non seulement des gènes hémizygotes, mais également celle des gènes euploïdes voisins. A l'échelle génomique, nous observons des phénomènes concordants. En effet, l'expression des gènes variant en nombre de copies est plus variable que celles des gènes ne variant pas. De plus, à l'instar de la délétion WBS, les CNV influencent l'expression des gènes adjacents, exerçant ainsi un impact global sur les profils d'expression dans les tissus. Résumé pour un large public : De nombreuses maladies ont pour cause un défaut génétique. Parmi les types de mutations, on compte la disparition (délétion) d'une partie de notre génome ou sa duplication. Bien que l'on connaisse les anomalies associées à certaines maladies, les mécanismes moléculaires par lesquels ces réarrangements de notre matériel génétique induisent les maladies sont encore méconnus. C'est pourquoi nous nous sommes intéressés à la régulation des gènes dans les régions susceptibles à délétion ou duplication. Dans ce travail, nous avons démontré que les délétions et les duplications influencent la régulation des gènes situés à proximité, et que ces changements interviennent dans plusieurs organes.

Peroxisomal Delta(3),Delta(2)-enoyl CoA isomerases and evolution of cytosolic paralogues in embryophytes

Delta(3),Delta(2)-enoyl CoA isomerase (ECI) is an enzyme that participates in the degradation of unsaturated fatty acids through the beta-oxidation cycle. Three genes encoding Delta(3),Delta(2)-enoyl CoA isomerases and named AtECI1, AtECI2 and AtECI3 have been identified in Arabidopsis thaliana. When expressed heterologously in Saccharomyces cerevisiae, all three ECI proteins were targeted to the peroxisomes and enabled the yeast Deltaeci1 mutant to degrade 10Z-heptadecenoic acid, demonstrating Delta(3),Delta(2)-enoyl CoA isomerase activity in vivo. Fusion proteins between yellow fluorescent protein and AtECI1 or AtECI2 were targeted to the peroxisomes in onion epidermal cells and Arabidopsis root cells, but a similar fusion protein with AtECI3 remained in the cytosol for both tissues. AtECI3 targeting to peroxisomes in S. cerevisiae was dependent on yeast PEX5, while expression of Arabidopsis PEX5 in yeast failed to target AtECI3 to peroxisomes. AtECI2 and AtECI3 are tandem duplicated genes and show a high level of amino acid conservation, except at the C-terminus; AtECI2 ends with the well conserved peroxisome targeting signal 1 (PTS1) terminal tripeptide PKL, while AtECI3 possesses a divergent HNL terminal tripeptide. Evolutionary analysis of ECI genes in plants revealed several independent duplication events, with duplications occurring in rice and Medicago truncatula, generating homologues with divergent C-termini and no recognizable PTS1. All plant ECI genes analyzed, including AtECI3, are under negative purifying selection, implying functionality of the cytosolic AtECI3. Analysis of the mammalian and fungal genomes failed to identify cytosolic variants of the Delta(3),Delta(2)-enoyl CoA isomerase, indicating that evolution of cytosolic Delta(3),Delta(2)-enoyl CoA isomerases is restricted to the plant kingdom