Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

Identification of HMX1 target genes: a predictive promoter model approach.

PURPOSE: A homozygous mutation in the H6 family homeobox 1 (HMX1) gene is responsible for a new oculoauricular defect leading to eye and auricular developmental abnormalities as well as early retinal degeneration (MIM 612109). However, the HMX1 pathway remains poorly understood, and in the first approach to better understand the pathway's function, we sought to identify the target genes. METHODS: We developed a predictive promoter model (PPM) approach using a comparative transcriptomic analysis in the retina at P15 of a mouse model lacking functional Hmx1 (dmbo mouse) and its respective wild-type. This PPM was based on the hypothesis that HMX1 binding site (HMX1-BS) clusters should be more represented in promoters of HMX1 target genes. The most differentially expressed genes in the microarray experiment that contained HMX1-BS clusters were used to generate the PPM, which was then statistically validated. Finally, we developed two genome-wide target prediction methods: one that focused on conserving PPM features in human and mouse and one that was based on the co-occurrence of HMX1-BS pairs fitting the PPM, in human or in mouse, independently. RESULTS: The PPM construction revealed that sarcoglycan, gamma (35kDa dystrophin-associated glycoprotein) (Sgcg), teashirt zinc finger homeobox 2 (Tshz2), and solute carrier family 6 (neurotransmitter transporter, glycine) (Slc6a9) genes represented Hmx1 targets in the mouse retina at P15. Moreover, the genome-wide target prediction revealed that mouse genes belonging to the retinal axon guidance pathway were targeted by Hmx1. Expression of these three genes was experimentally validated using a quantitative reverse transcription PCR approach. The inhibitory activity of Hmx1 on Sgcg, as well as protein tyrosine phosphatase, receptor type, O (Ptpro) and Sema3f, two targets identified by the PPM, were validated with luciferase assay. CONCLUSIONS: Gene expression analysis between wild-type and dmbo mice allowed us to develop a PPM that identified the first target genes of Hmx1.

Amplification of 16S rRNA gene sequences to differentiate two highly related bradyrhizobia species

A 16S rRNA gene PCR-based assay was developed aiming at a fast molecular diagnostic method to differentiate the two phylogenetically closely related species Bradyrhizobium japonicum and B. elkanii, isolated from soybean nodules, in order to identify those more competitive and comprising greater nitrogen fixation ability for use in the formulation of commercial inoculants. The assay used was able to discriminate ten reference strains belonging to these two Bradyrhizobium species, as well as to efficiently identify 37 strains isolated from fields cultivated with soybean.

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model.

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

Gene expression profiling reveals consistent differences between clinical samples of human leukaemias and their model cell lines.

Microarray gene expression profiles of fresh clinical samples of chronic myeloid leukaemia in chronic phase, acute promyelocytic leukaemia and acute monocytic leukaemia were compared with profiles from cell lines representing the corresponding types of leukaemia (K562, NB4, HL60). In a hierarchical clustering analysis, all clinical samples clustered separately from the cell lines, regardless of leukaemic subtype. Gene ontology analysis showed that cell lines chiefly overexpressed genes related to macromolecular metabolism, whereas in clinical samples genes related to the immune response were abundantly expressed. These findings must be taken into consideration when conclusions from cell line-based studies are extrapolated to patients.

Gene transfer of a soluble IL-1 type 2 receptor-Ig fusion protein improves cardiac allograft survival in rats.

OBJECTIVE: Interleukin-1 (IL-1) mediates ischemia-reperfusion injury and graft inflammation after heart transplantation. IL-1 affects target cells through two distinct types of transmembrane receptors, type-1 receptor (IL-1R1), which transduces the signal, and the non-signaling type-2 receptor (IL-1R2), which acts as a ligand sink that subtracts IL-1beta from IL-1R1. We analyzed the efficacy of adenovirus (Ad)-mediated gene transfer of a soluble IL-1R2-Ig fusion protein in delaying cardiac allograft rejection and the mechanisms underlying the protective effect. METHODS: IL-1 inhibition by IL-1R2-Ig was tested using an in vitro functional assay whereby endothelial cells preincubated with AdIL-1R2-Ig or control virus were stimulated with recombinant IL-1beta or tumor necrosis factor-alpha (TNF-alpha), and urokinase-type plasminogen activator (u-PA) induction was measured by zymography. AdIL-1R2-Ig was delivered to F344 rat donor hearts ex vivo, which were placed in the abdominal position in LEW hosts. Intragraft inflammatory cell infiltrates and proinflammatory cytokine expression were analyzed by immunohistochemistry and real-time reverse transcriptase-polymerase chain reaction (RT-PCR), respectively. RESULTS: IL-1R2-Ig specifically inhibited IL-1beta-induced u-PA responses in vitro. IL-1R2-Ig gene transfer reduced intragraft monocytes/macrophages and CD4(+) cell infiltrates (p<0.05), TNF-alpha and transforming growth factor-beta (TGF-beta) expression (p<0.05), and prolonged graft survival (15.6+/-5.7 vs 10.3+/-2.5 days with control vector and 10.1+/-2.1 days with buffer alone; p<0.01). AdIL-1R2-Ig combined with a subtherapeutic regimen of cyclosporin A (CsA) was superior to CsA alone (19.4+/-3.0 vs 15.9+/-1.8 days; p<0.05). CONCLUSIONS: Soluble IL-1 type-2 receptor gene transfer attenuates cardiac allograft rejection in a rat model. IL-1 inhibition may be useful as an adjuvant therapy in heart transplantation.

Abnormal social behaviors and altered gene expression rates in a mouse model for Potocki-Lupski syndrome.

The Potocki-Lupski syndrome (PTLS) is associated with a microduplication of 17p11.2. Clinical features include multiple congenital and neurobehavioral abnormalities and autistic features. We have generated a PTLS mouse model, Dp(11)17/+, that recapitulates some of the physical and neurobehavioral phenotypes present in patients. Here, we investigated the social behavior and gene expression pattern of this mouse model in a pure C57BL/6-Tyr(c-Brd) genetic background. Dp(11)17/+ male mice displayed normal home-cage behavior but increased anxiety and increased dominant behavior in specific tests. A subtle impairment in the preference for a social target versus an inanimate target and abnormal preference for social novelty (the preference to explore an unfamiliar mouse versus a familiar one) was also observed. Our results indicate that these animals could provide a valuable model to identify the specific gene(s) that confer abnormal social behaviors and that map within this delimited genomic deletion interval. In a first attempt to identify candidate genes and for elucidating the mechanisms of regulation of these important phenotypes, we directly assessed the relative transcription of genes within and around this genomic interval. In this mouse model, we found that candidates genes include not only most of the duplicated genes, but also normal-copy genes that flank the engineered interval; both categories of genes showed altered expression levels in the hippocampus of Dp(11)17/+ mice.

An operon of three transcriptional regulators controls horizontal gene transfer of the integrative and conjugative element ICEclc in Pseudomonas knackmussii B13.

The integrative and conjugative element ICEclc is a mobile genetic element in Pseudomonas knackmussii B13, and an experimental model for a widely distributed group of elements in Proteobacteria. ICEclc is transferred from specialized transfer competent cells, which arise at a frequency of 3-5% in a population at stationary phase. Very little is known about the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene operon encoded by ICEclc, which exerts global control on transfer initiation. The operon consists of three consecutive regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that MfsR autoregulates expression of the operon, whereas TciR is a global activator of ICEclc gene expression, but no clear role was yet found for MarR. Deletion of mfsR increases expression of tciR and marR, causing the proportion of transfer competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two orders of magnitude loss in population viability, individual cell growth arrest and loss of ICEclc. This indicates that autoregulation is an important feature maintaining ICE transfer but avoiding fitness loss. Bioinformatic analysis showed that the mfsR-marR-tciR operon is unique for ICEclc and a few highly related ICE, whereas tciR orthologues occur more widely in a large variety of suspected ICE among Proteobacteria.

Calreticulin levels determine onset of early muscle denervation by fast motoneurons of ALS model mice.

Although the precise signaling mechanisms underlying the vulnerability of some sub-populations of motoneurons in ALS remain unclear, critical factors such as metallo-proteinase 9 expression, neuronal activity and endoplasmic reticulum stress have been shown to be involved. In the context of SOD1(G93A) ALS mouse model, we previously showed that a two-fold decrease in calreticulin (CRT) is occurring in the vulnerable fast motoneurons. Here, we asked to which extent the decrease in CRT levels was causative to muscle denervation and/or motoneuron degeneration. Toward this goal, a hemizygous deletion of the crt gene in SOD1(G93A) mice was generated since the complete ablation of crt is embryonic lethal. We observed that SOD1(G93A);crt(+/-) mice display increased and earlier muscle weakness and muscle denervation compared to SOD1(G93A) mice. While CRT reduction in motoneurons leads to a strong upregulation of two factors important in motoneuron dysfunction, ER stress and mTOR activation, it does not aggravate motoneuron death. Our results underline a prevalent role for CRT levels in the early phase of muscle denervation and support CRT regulation as a potential therapeutic approach.

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.

Synthesis of the two isomers of the potential sex pheromone of Thaumetopoea pityocampa (Lepidoptera, Notodontidae) and related model compounds

The synthesis of the major component of the sex pheromone secretion of the processionary moth, Tkawnztopoeja pltyocampa (Denis and Schiff.) (Lepidoptera, Notodontidae), (Z)-13-hexadecen-ll-ynyl acetate (1), the corresponding (E)-isomer (2) and the four structurally related model compounds (Z⁄E,Z,Z)-5,9,13-hexadecatrienyl acetate (3), (Z⁄E,Z,Z)-3,7,ll-hexadecatrienyl acetate (4), (Z⁄E,E,Z)-7,9,13-hexadecatrienyl acetate (5) and (Z)-7-hexadecen-5-ynyl acetate (6) is described.

The Pseudomonas aeruginosa antimetabolite L -2-amino-4-methoxy-trans-3-butenoic acid (AMB) is made from glutamate and two alanine residues via a thiotemplate-linked tripeptide precursor.

The Pseudomonas aeruginosa toxin L-2-amino-4-methoxy-trans-3-butenoic acid (AMB) is a non-proteinogenic amino acid which is toxic for prokaryotes and eukaryotes. Production of AMB requires a five-gene cluster encoding a putative LysE-type transporter (AmbA), two non-ribosomal peptide synthetases (AmbB and AmbE), and two iron(II)/α-ketoglutarate-dependent oxygenases (AmbC and AmbD). Bioinformatics analysis predicts one thiolation (T) domain for AmbB and two T domains (T1 and T2) for AmbE, suggesting that AMB is generated by a processing step from a precursor tripeptide assembled on a thiotemplate. Using a combination of ATP-PPi exchange assays, aminoacylation assays, and mass spectrometry-based analysis of enzyme-bound substrates and pathway intermediates, the AmbB substrate was identified to be L-alanine (L-Ala), while the T1 and T2 domains of AmbE were loaded with L-glutamate (L-Glu) and L-Ala, respectively. Loading of L-Ala at T2 of AmbE occurred only in the presence of AmbB, indicative of a trans loading mechanism. In vitro assays performed with AmbB and AmbE revealed the dipeptide L-Glu-L-Ala at T1 and the tripeptide L-Ala-L-Glu-L-Ala attached at T2. When AmbC and AmbD were included in the assay, these peptides were no longer detected. Instead, an L-Ala-AMB-L-Ala tripeptide was found at T2. These data are in agreement with a biosynthetic model in which L-Glu is converted into AMB by the action of AmbC, AmbD, and tailoring domains of AmbE. The importance of the flanking L-Ala residues in the precursor tripeptide is discussed.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Learning-induced gene expression in the heads of two Nasonia species that differ in long-term memory formation.

BACKGROUND: Cellular processes underlying memory formation are evolutionary conserved, but natural variation in memory dynamics between animal species or populations is common. The genetic basis of this fascinating phenomenon is poorly understood. Closely related species of Nasonia parasitic wasps differ in long-term memory (LTM) formation: N. vitripennis will form transcription-dependent LTM after a single conditioning trial, whereas the closely-related species N. giraulti will not. Genes that were differentially expressed (DE) after conditioning in N. vitripennis, but not in N. giraulti, were identified as candidate genes that may regulate LTM formation. RESULTS: RNA was collected from heads of both species before and immediately, 4 or 24 hours after conditioning, with 3 replicates per time point. It was sequenced strand-specifically, which allows distinguishing sense from antisense transcripts and improves the quality of expression analyses. We determined conditioning-induced DE compared to naïve controls for both species. These expression patterns were then analysed with GO enrichment analyses for each species and time point, which demonstrated an enrichment of signalling-related genes immediately after conditioning in N. vitripennis only. Analyses of known LTM genes and genes with an opposing expression pattern between the two species revealed additional candidate genes for the difference in LTM formation. These include genes from various signalling cascades, including several members of the Ras and PI3 kinase signalling pathways, and glutamate receptors. Interestingly, several other known LTM genes were exclusively differentially expressed in N. giraulti, which may indicate an LTM-inhibitory mechanism. Among the DE transcripts were also antisense transcripts. Furthermore, antisense transcripts aligning to a number of known memory genes were detected, which may have a role in regulating these genes. CONCLUSION: This study is the first to describe and compare expression patterns of both protein-coding and antisense transcripts, at different time points after conditioning, of two closely related animal species that differ in LTM formation. Several candidate genes that may regulate differences in LTM have been identified. This transcriptome analysis is a valuable resource for future in-depth studies to elucidate the role of candidate genes and antisense transcription in natural variation in LTM formation.