Genome Biology and Evolution of the Animal Kingdom

Estudi realitzat a partir d’una estada a la Institut J.W. Jenkinson Laboratory for Evolution and Development of the University of Oxford, Regne Unit, entre 2010 i 2012. He estat membre del laboratori del Professor Peter W.H. Holland com a becari post-doctoral Beatriu de Pinós des de setembre de 2010 al setembre de 2012. El nostre projecte de recerca se centra en l'anàlisi genòmic comparatiu del Regne Animal, tot explorant el contingut dels genomes a través de totes les branques de l'arbre dels animals. Totes les referències a les meves publicacions durant aquest post-doc es poden trobar a http://about.me/jordi_paps. Crec que el nombre i la qualitat dels resultats del meu post-doc, un total de 8 publicacions incloent dos articles a la prestigiosa revista Nature, són prova de l'èxit d'aquest post-doc. Prof Peter W. H. Holland (Departament de Zoologia de la Universitat d'Oxford) i jo som coautors de tres articles de genòmica comparativa, resultats directes d'aquest projecte: 1) comparació de families gèniques entre vertebrats invertebrats (Briefings in Functional Genomics), 2) el genoma de l'ostra (publicat a la revista Nature), i 3) els genomes de 6 platihelmints paràsits (acceptat també a Nature). A més, tenim altres 2 treballs en preparació. Un d'ells analitza l'evolució, expressió i funció dels gens Hox al a la tènia Hymenolepis. El perfil fi d'aquests gens clau del desenvolupament esclareix els canvis d'estil de vida dels organismes. A més, durant aquest últim post-doc he participat en diverses col•laboracions, incloent anàlisi de gens d'envelliment a cucs plans, un estudi sobre la filogènia del grup Gastrotricha, una revisió de l'evolució phylum Platyhelminthes, així com un capítol d'un llibre sobre l'evolució dels animals bilaterals. Finalment, gràcies a la beca Beatriu de Pinós, el Prof. Peter W.H. Holland m'ha convidat a formar part del seu equip com un investigador post-doctoral en el seu projecte ERC Advance actual sobre duplicacions genòmiques.

The genome of the fire ant Solenopsis invicta.

Ants have evolved very complex societies and are key ecosystem members. Some ants, such as the fire ant Solenopsis invicta, are also major pests. Here, we present a draft genome of S. invicta, assembled from Roche 454 and Illumina sequencing reads obtained from a focal haploid male and his brothers. We used comparative genomic methods to obtain insight into the unique features of the S. invicta genome. For example, we found that this genome harbors four adjacent copies of vitellogenin. A phylogenetic analysis revealed that an ancestral vitellogenin gene first underwent a duplication that was followed by possibly independent duplications of each of the daughter vitellogenins. The vitellogenin genes have undergone subfunctionalization with queen- and worker-specific expression, possibly reflecting differential selection acting on the queen and worker castes. Additionally, we identified more than 400 putative olfactory receptors of which at least 297 are intact. This represents the largest repertoire reported so far in insects. S. invicta also harbors an expansion of a specific family of lipid-processing genes, two putative orthologs to the transformer/feminizer sex differentiation gene, a functional DNA methylation system, and a single putative telomerase ortholog. EST data indicate that this S. invicta telomerase ortholog has at least four spliceforms that differ in their use of two sets of mutually exclusive exons. Some of these and other unique aspects of the fire ant genome are likely linked to the complex social behavior of this species.

Evolving olfactory systems on the fly.

The detection of odour stimuli in the environment is universally important for primal behaviours such as feeding, mating, kin interactions and escape responses. Given the ubiquity of many airborne chemical signals and the similar organisation of animal olfactory circuits, a fundamental question in our understanding of the sense of smell is how species-specific behavioural responses to odorants can evolve. Recent comparative genomic, developmental and physiological studies are shedding light on this problem by providing insights into the genetic mechanisms that underlie anatomical and functional evolution of the olfactory system. Here we synthesise these data, with a particular focus on insect olfaction, to address how new olfactory receptors and circuits might arise and diverge, offering glimpses into how odour-evoked behaviours could adapt to an ever-changing chemosensory world.

Homologous sex chromosomes in three deeply divergent anuran species.

Comparative genomic studies are revealing that, in sharp contrast with the strong stability found in birds and mammals, sex determination mechanisms are surprisingly labile in cold-blooded vertebrates, with frequent transitions between different pairs of sex chromosomes. It was recently suggested that, in context of this high turnover, some chromosome pairs might be more likely than others to be co-opted as sex chromosomes. Empirical support, however, is still very limited. Here we show that sex-linked markers from three highly divergent groups of anurans map to Xenopus tropicalis scaffold 1, a large part of which is homologous to the avian sex chromosome. Accordingly, the bird sex determination gene DMRT1, known to play a key role in sex differentiation across many animal lineages, is sex linked in all three groups. Our data provide strong support for the idea that some chromosome pairs are more likely than others to be co-opted as sex chromosomes because they harbor key genes from the sex determination pathway.

Evolutionary comparison provides evidence for pathogenicity of RMRP mutations.

Cartilage-hair hypoplasia (CHH) is a pleiotropic disease caused by recessive mutations in the RMRP gene that result in a wide spectrum of manifestations including short stature, sparse hair, metaphyseal dysplasia, anemia, immune deficiency, and increased incidence of cancer. Molecular diagnosis of CHH has implications for management, prognosis, follow-up, and genetic counseling of affected patients and their families. We report 20 novel mutations in 36 patients with CHH and describe the associated phenotypic spectrum. Given the high mutational heterogeneity (62 mutations reported to date), the high frequency of variations in the region (eight single nucleotide polymorphisms in and around RMRP), and the fact that RMRP is not translated into protein, prediction of mutation pathogenicity is difficult. We addressed this issue by a comparative genomic approach and aligned the genomic sequences of RMRP gene in the entire class of mammals. We found that putative pathogenic mutations are located in highly conserved nucleotides, whereas polymorphisms are located in non-conserved positions. We conclude that the abundance of variations in this small gene is remarkable and at odds with its high conservation through species; it is unclear whether these variations are caused by a high local mutation rate, a failure of repair mechanisms, or a relaxed selective pressure. The marked diversity of mutations in RMRP and the low homozygosity rate in our patient population indicate that CHH is more common than previously estimated, but may go unrecognized because of its variable clinical presentation. Thus, RMRP molecular testing may be indicated in individuals with isolated metaphyseal dysplasia, anemia, or immune dysregulation.

Control of gene expression in melanocytes and RPE by conserved distal regulatory elements

Résumé Durant le développement embryonnaire, les cellules pigmentaires des mammifères se développent à partir de deux origines différentes : les melanocytes se développent à partir de la crête neurale alors que les cellules de la rétine pigmentaire (RP) ont une origine neuronale. Un grand nombre de gènes sont impliqués dans la pigmentation dont les gènes de la famille tyrosinase à savoir Tyr, Tyrp1 et Dct. Certaines études ont suggéré que les gènes de la pigmentation sont régulés de manière différentielle dans les mélanocytes et dans la RP. Dans ce travail, les gènes de la famille tyrosinase ont été étudiés comme modèle de la régulation des gènes de la pigmentation par des éléments régulateurs agissant à distance. II a été montré que le promoteur du gène Tyrp1pouvait induire l'expression d'un transgène uniquement dans la RP alors que ce gène est aussi exprimé dans les mélanocytes comme le montre le phénotype des souris mutantes pour Tyrp1. Ce résultat suggère que les éléments régulateurs du promoteur sont suffisants pour l'expression dans la RP mais pas pour l'expression dans les mélanocytes. J'ai donc cherché à identifier la séquence qui régule l'expression dans les mélanocytes. Un chromosome artificiel bactérien (CAB) contenant le gène Tyrp1 s'est avéré suffisant pour induire l'expression dans les mélanocytes, comme démontré par la correction du phénotype mutant. La séquence de ce CAB contient plusieurs régions très conservées qui pourraient représenter de nouveaux éléments régulateurs. Par la suite, j'ai focalisé mon analyse sur une séquence située à -I5 kb qui s'est révélée être un amplificateur spécifique aux mélanocytes comme démontré par des expériences de cultures cellulaire et de transgenèse. De plus, une analyse poussée de cet élément a révélé que le facteur de transcription Sox 10 représentait un transactivateur de cet amplificateur. Comme pour Tyrp1, la régulation du gène tyrosinase est contrôlée par différents éléments régulateurs dans les mélanocytes et la RP. Il a été montré que le promoteur de tyrosinase n'était pas suffisant pour une forte expression dans les mélanocytes et la RP. De plus, l'analyse de la région située en amont a révélé la présence d'un amplificateur nécessaire à l'expression dans les mélanocytes à la position -15 kb. Cet amplificateur n'est toutefois pas actif dans la RP mais agit comme un répresseur dans ces cellules. Ces résultats indiquent que certains éléments nécessaires à l'expression dans les deux types de cellules pigmentaires sont absents de ces constructions. Comme pour Tyrp1, j'ai en premier lieu démontré qu'un CAB était capable de corriger le phénotype albinique, puis ai inséré un gène reporter (lacZ) dans le CAB par recombinaison homologue et ai finalement analysé l'expression du reporter en transgenèse. Ces souris ont montré une expression forte du lacZ dans les mélanocytes et la RP, ce qui indique que le CAB contient les séquences régulatrices nécessaires à l'expression correcte de tyrosinase. Afin de localiser plus précisément les éléments régulateurs, j'ai ensuite généré des délétions dans le CAB et analysé l'expression du lacZ en transgenèse. La comparaison de séquences génomiques provenant de différentes espèces a permis par la suite d'identifier des régions représentant de nouveaux éléments régulateurs potentiels. En utilisant cette approche, j'ai identifié une région qui se comporte comme un amplificateur dans la RP et qui est nécessaire à l'expression de tyrosinase dans ce tissu. De plus, j'ai identifié les facteurs de transcription Mitf et Sox10 comme transactivateurs de l'amplificateur spécifique aux mélanocytes situé à -15 kb. L'identification et la caractérisation des ces éléments régulateurs des gènes tyrosinase et Tyrp1confirme donc que la régulation différentielle des gènes dans les mélanocytes et la RP est liée à des éléments régulateurs séparés. Summary Pigment cells of mammals originate from two different lineages: melanocytes arise from the neural crest, whereas cells of the retinal pigment epithelium (RPE) originate from the optic cup of the developing forebrain. A large set of genes are involved in pigmentation, including the members of the tyrosinase gene family, namely tyrosinase, Tyrp1 and Dct. Previous studies have suggested that pigmentation genes are differentially regulated in melanocytes and RPE. In this work, the tyrosinase gene family was used as a model for studying the involvement of distal regulatory elements in pigment cell-specific gene expression. The promoter of the Tyrp1 gene has been shown to drive detectable transgene expression only to the RPE, even though the gene is also expressed in melanocytes as evident from Tyrp1-mutant mice. This indicates that the regulatory elements responsible for Tyrp1 gene expression in the RPE are not sufficient for expression in melanocytes. I thus searched for a putative melanocyte-specific regulatory sequence and demonstrate that a bacterial artificial chromosome (BAC) containing the Tyrp1 gene and surrounding sequences is able to target transgenic expression to melanocytes and to rescue the Tyrp1 b (brown) phenotype. This BAC contains several highly conserved non-coding sequences that might represent novel regulatory elements. I further focused on a sequence located at -15 kb which I identified as amelanocyte-specific enhancer as shown by cell culture and transgenic mice. In addition, further functional analysis identified the transcription factor Sox10 as being able to bind and transactivate this enhancer. As for Tyrp1, tyrosinase gene regulation is mediated by different cis-regulatory elements in melanocytes and RPE. It was shown that the tyrosinase promoter was not sufficient to confer strong and specific expression in melanocytes and RPE. Moreover, analysis of tyrosinase upstream sequence, revealed the presence of a specific enhancer at position -15 kb which was necessary to confer strong expression in melanocytes. This enhancer element however failed to act as an enhancer in the RPE, but rather repressed expression. This indicates that some regulatory elements required for tyrosinase expression in both RPE and melanocytes are still missing from these constructs. As for Tyrp1, I first demonstrated that a BAC containing the Tyr gene is able to rescue the Tyr c (albino) phenotype in mice, then I inserted a lacZ reporter gene in the BAC by homologous recombination, and finally analysed the pattern of lacZ expression in transgenic mice. These mice showed strong lacZ expression in both RPE and melanocytes, indicating that the BAC contains the regulatory sequences required for proper tyrosinase expression. In order to localize more precisely these regulatory elements, I have then generated several deletions in the BAC and analysed lacZ expression in transgenic mice. Multi-species comparative genomic analysis then allowed identifying conserved sequences that potentially represent novel regulatory elements. Using this experimental approach, I identified a region that behaves as a RPE-specific enhancer and that is required for tyrosinase expression in the retina] pigment epithelium. In addition, I identified the transcription factors Mitf and Sox l0 as being transactivators of the melanocyte-specific enhancer located at -l5 kb. The identification and characterization of these tyrosinase and Tyrp1 distal regulatory element supports the idea that separate regulatory sequences mediate differential gene expression in melanocytes and RPE.

Evolution of Bordetella pertussis post vaccination

Evolution of Bordetella pertussis post vaccination Whooping cough or pertussis is caused by the gram-negative bacterium Bordetella pertussis. It is a highly contiguous disease in the human respiratory tract. Characteristic of pertussis is a paroxysmal cough with whooping sound during gasps of breath after coughing episodes. It is potentially fatal to unvaccinated infants. The best approach to fight pertussis is to vaccinate. Vaccinations against pertussis have been available from the 1940s. Traditionally vaccines were whole-cell pertussis (wP) preparations as part of the combined diphtheria-tetanus-pertussis (DTP) vaccines. More recently acellular pertussis (aP) vaccines have replaced the wP vaccines in many countries. The aP vaccines are less reactogenic and can also be administered to school children and adults. There are several publications reporting variation in the i>B. pertussis virulence factors that are also aP vaccine antigens. This has occurred in the genes coding for pertussis toxin and pertactin about 15 to 30 years after the introduction of pertussis vaccines to immunisation programs. Resurgence of pertussis has also been reported in many countries with high vaccination coverage. In this study the evolution of B. pertussis was investigated in Finland, the United Kingdom, Poland, Serbia, China, Senegal and Kenya. These represent countries with a long history of high vaccination coverage with stable vaccines or changes in the vaccine formulation; countries which established high vaccination coverage late; and countries where vaccinations against pertussis were started late. With bacterial cytotoxicity and cytokine measurements, comparative genomic hybridisation, pulsed-field gel electrophoresis (PFGE), genotyping and serotyping it was found that changes in the vaccine composition can postpone the emergence of antigenic variants. It seems that the change in PFGE profiles and the loss of genetic material in the genome of B. pertussis are similar in most countries and the vaccine-induced immunity is selecting non-vaccine type strains. However, the differences in the formulation of the vaccines, the vaccination programs and in the coverage of pertussis vaccination have affected the speed and timing of these changes.

Upconverting diagnostics: Multiplex assays utilizing upconversion luminescence technology

Conventional diagnostics tests and technologies typically allow only a single analysis and result per test. The aim of this study was to propose robust and multiplex array-inwell test platforms based on oligonucleotide and protein arrays combining the advantages of simple instrumentation and upconverting phosphor (UCP) reporter technology. The UCPs are luminescent lanthanide-doped crystals that have a unique capability to convert infrared radiation into visible light. No autofluorescence is produced from the sample under infrared excitation enabling the development of highly sensitive assays. In this study, an oligonucleotide array-in-well hybridization assay was developed for the detection and genotyping of human adenoviruses. The study provided a verification of the advantages and potential of the UCP-based reporter technology in multiplex assays as well as anti-Stokes photoluminescence detection with a new anti- Stokes photoluminescence imager. The developed assay was technically improved and used to detect and genotype adenovirus types from clinical specimens. Based on the results of the epidemiological study, an outbreak of adenovirus type B03 was observed in the autumn of 2010. A quantitative array-in-well immunoassay was developed for three target analytes (prostate specific antigen, thyroid stimulating hormone, and luteinizing hormone). In this study, quantitative results were obtained for each analyte and the analytical sensitivities in buffer were in clinically relevant range. Another protein-based array-inwell assay was developed for multiplex serodiagnostics. The developed assay was able to detect parvovirus B19 IgG and adenovirus IgG antibodies simultaneously from serum samples according to reference assays. The study demonstrated that the UCPtechnology is a robust detection method for diverse multiplex imaging-based array-inwell assays.

Campylobacter dans différents environnements aquatiques : quantification et génotypage afin de mieux évaluer les risques potentiels d’infection pour l’être humain

Campylobacter est l’agent pathogène zoonotique responsable de la majorité des gastro-entérites d’origine bactérienne chez l’homme. Les produits de volaille représentent la principale source d’infection; toutefois, l’exposition peut également découler de contacts directs avec les animaux ou avec l’eau. Une forte variation saisonnière est présente dans les cas rapportés, qui n’est toujours pas élucidée : les eaux environnementales, sources d’infection connues, sont soupçonnées. Cette étude transversale a été réalisée dans la région Sud-Est du Québec (Canada) où Campylobacter fut quantifié et génotypé à partir de différentes sources d’eau (eaux de captage, récréatives et usées) et de cas cliniques afin d’évaluer les risques potentiels posé par l’eau environnementale. Différents essais PCR en temps réel furent appliqués à l’eau environnementale et comparés: 2 ont été sélectionnés pour leur spécificité et sensibilité de quantification. Les courbes standards ont été calibrées en utilisant la PCR digitale pour déterminer précisément les concentrations. Les isolats environnementaux et cliniques furent comparés génétiquement en utilisant le CGF (« comparative genomic fingerprinting »). Les eaux usées étaient plus contaminées que les eaux de captage et récréatives (3.9Log, 1.7Log et 1.0Log cellules/L en moyenne, respectivement). Six pour cent des isolats d’eaux environnementales étaient génétiquement similaires (100 % homologie) aux isolats cliniques. Les cas cliniques de campylobactériose d’été montraient des isolats avec davantage de similarités génétiques avec les isolats retrouvés dans l’eau environnementale comparativement aux autres saisons (p<0.01). Les faibles concentrations et similarités génétiques entre les isolats d’eau et cliniques suggèrent un risque de transmission possible, mais faible.

Origin of avian genome size and structure in non-avian dinosaurs

Avian genomes are small and streamlined compared with those of other amniotes by virtue of having fewer repetitive elements and less non-coding DNA(1,2). This condition has been suggested to represent a key adaptation for flight in birds, by reducing the metabolic costs associated with having large genome and cell sizes(3,4). However, the evolution of genome architecture in birds, or any other lineage, is difficult to study because genomic information is often absent for long-extinct relatives. Here we use a novel bayesian comparative method to show that bone-cell size correlates well with genome size in extant vertebrates, and hence use this relationship to estimate the genome sizes of 31 species of extinct dinosaur, including several species of extinct birds. Our results indicate that the small genomes typically associated with avian flight evolved in the saurischian dinosaur lineage between 230 and 250 million years ago, long before this lineage gave rise to the first birds. By comparison, ornithischian dinosaurs are inferred to have had much larger genomes, which were probably typical for ancestral Dinosauria. Using comparative genomic data, we estimate that genome-wide interspersed mobile elements, a class of repetitive DNA, comprised 5 - 12% of the total genome size in the saurischian dinosaur lineage, but was 7 - 19% of total genome size in ornithischian dinosaurs, suggesting that repetitive elements became less active in the saurischian lineage. These genomic characteristics should be added to the list of attributes previously considered avian but now thought to have arisen in non-avian dinosaurs, such as feathers(5), pulmonary innovations 6, and parental care and nesting

The leptospiral antigen Lp49 is a two-domain protein with putative protein binding function

Pathogenic Leptospira is the etiological agent of leptospirosis, a life-threatening disease that affects populations worldwide. Currently available vaccines have limited effectiveness and therapeutic interventions are complicated by the difficulty in making an early diagnosis of leptospirosis. The genome of Leptospira interrogans was recently sequenced and comparative genomic analysis contributed to the identification of surface antigens, potential candidates for development of new vaccines and serodiagnosis. Lp49 is a membrane-associated protein recognized by antibodies present in sera from early and convalescent phases of leptospirosis patients. Its crystal structure was determined by single-wavelength anomalous diffraction using selenomethionine-labelled crystals and refined at 2.0 angstrom resolution. Lp49 is composed of two domains and belongs to the all-beta-proteins class. The N-terminal domain folds in an immunoglobulin-like beta-sandwich structure, whereas the C-terminal domain presents a seven-bladed beta-propeller fold. Structural analysis of Lp49 indicates putative protein-protein binding sites, suggesting a role in Leptospira-host interaction. This is the first crystal structure of a leptospiral antigen described to date. (C) 2008 Elsevier Inc. All rights reserved.

alpha-Hydroxynitrile lyase protein from Xylella fastidiosa: Cloning, expression, and characterization

Xylella fastidiosa is a xylem-restricted plant pathogen that causes a range of diseases in several and important crops. Through comparative genomic sequence analysis many genes were identified and, among them, several potentially involved in plant-pathogen interaction. The experimental determination of the primary sequence of some markedly expressed proteins for X fastidiosa and the comparison with the nucleic acids sequence of genome identified one of them as being SCJ21.16 (XFa0032) gene product. The comparative analysis of this protein against SWISSPROT database, in special, resulted in similarity with a-hydroxynitrile lyase enzyme (HNL) from Arabidopsis thaliana, causing interest for being one of the most abundant proteins both in the whole cell extract as well as in the extracellular protein fraction. It is known that HNL enzyme are involved in a process termed ""cyanogenesis"", which catalyzes the dissociation of alpha-hydroxinitrile into carbonyle and HCN when plant tissue is damaged. Although the complete genome sequences of X.fastidiosa are available and the cyanogenesis process is well known, the biological role of this protein in this organism is not yet functionally characterized. In this study we presented the cloning, expression, characterization of recombinant HNL from X fastidiosa, and its probable function in the cellular metabolism. The successful cloning and heterologous expression in Escherichia coli resulted in a satisfactory amount of the recombinant HNL expressed in a soluble, and active form giving convenient access to pure enzyme for biochemical and structural studies. Finally, our results confirmed that the product of the gene XFa0032 can be positively assigned as FAD-independent HNLs. (C) 2009 Elsevier Ltd. All rights reserved.

alpha-Hydroxynitrile lyase protein from Xylella fastidiosa: Cloning, expression, and characterization

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Molecular and morphological data support the existence of a sexual cycle in species of the genus Paracoccidioides

The genus Paracoccidioides includes the thermodimorphic species Paracoccidioides brasiliensis and P. lutzii, both of which are etiologic agents of paracoccidioidomycosis, a systemic mycosis that affects humans in Latin America. Despite the common occurrence of a sexual stage among closely related fungi, this has not been observed with Paracoccidioides species, which have thus been considered asexual. Molecular evolutionary studies revealed recombination events within isolated populations of the genus Paracoccidioides, suggesting the possible existence of a sexual cycle. Comparative genomic analysis of all dimorphic fungi and Saccharomyces cerevisiae demonstrated the presence of conserved genes involved in sexual reproduction, including those encoding mating regulators such as MAT, pheromone receptors, pheromone-processing enzymes, and mating signaling regulators. The expression of sex-related genes in the yeast and mycelial phases of both Paracoccidioides species was also detected by realtime PCR, with nearly all of these genes being expressed preferentially in the filamentous form of the pathogens. In addition, the expression of sex-related genes was responsive to the putative presence of pheromone in the supernatants obtained from previous cocultures of strains of two different mating types. In vitro crossing of isolates of different mating types, discriminated by phylogenetic analysis of the α-box (MAT1-1) and the high-mobility-group (HMG) domain (MAT1-2), led to the identification of the formation of young ascocarps with constricted coiled hyphae related to the initial stage of mating. These genomic and morphological analyses strongly support the existence of a sexual cycle in species of the genus Paracoccidioides. © 2013, American Society for Microbiology.

Identification of duplicated and stress-inducible Aox2b gene co-expressed with Aox1 in species of the Medicago genus reveals a regulation linked to gene rearrangement in leguminous genomes

In flowering plants, alternative oxidase (Aox) is encoded by 3-5 genes distributed in 2 subfamilies (Aox1 and Aox2). In several species only Aox1 is reported as a stress-responsive gene, but in the leguminous Vigna unguiculata Aox2b is also induced by stress. In this work we investigated the Aox genes from two leguminous species of the Medicago genus (Medicago sativa and Medicago truncatula) which present one Aox1, one Aox2a and an Aox2b duplication (named here Aox2b1 and Aox2b2). Expression analyses by semi-quantitative RT-PCR in M. sativa revealed that Aox1, Aox2b1 and Aox2b2 transcripts increased during seed germination. Similar analyses in leaves and roots under different treatments (SA, PEG, H2O2 and cysteine) revealed that these genes are also induced by stress, but with peculiar spatio-temporal differences. Aox1 and Aox2b1 showed basal levels of expression under control conditions and were induced by stress in leaves and roots. Aox2b2 presented a dual behavior, i.e., it was expressed only under stress conditions in leaves, and showed basal expression levels in roots that were induced by stress. Moreover, Aox2a was expressed at higher levels in leaves and during seed germination than in roots and appeared to be not responsive to stress. The Aox expression profiles obtained from a M. truncatula microarray dataset also revealed a stress-induced co-expression of Aox1, Aox2b1 and Aox2b2 in leaves and roots. These results reinforce the stress-inducible co-expression of Aox1/Aox2b in some leguminous plants. Comparative genomic analysis indicates that this regulation is linked to Aox1/Aox2b proximity in the genome as a result of the gene rearrangement that occurred in some leguminous plants during evolution. The differential expression of Aox2b1/2b2 suggests that a second gene has been originated by recent gene duplication with neofunctionalization. © 2013 Elsevier GmbH. All rights reserved.