The Chara Genome: Secondary Complexity and Implications for Plant Terrestrialization.

Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.

A catalogue of recombination coldspots in interspecific tomato hybrids.

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae.

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Genome architecture and genetic diversity of allopolyploid okra (Abelmoschus esculentus).

The allopolyploid okra (Abelmoschus esculentus) unveiled telomeric repeats flanking distal gene-rich regions and short interstitial TTTAGGG telomeric repeats, possibly representing hallmarks of chromosomal speciation. Ribosomal RNA (rRNA) genes organize into 5S clusters, distinct from the 18S-5.8S-28S units, indicating an S-type rRNA gene arrangement. The assembly, in line with cytogenetic and cytometry observations, identifies 65 chromosomes and a 1.45 Gb genome size estimate in a haploid sibling. The lack of aberrant meiotic configurations implies limited to no recombination among sub-genomes. k-mer distribution analysis reveals 75% has a diploid nature and 15% heterozygosity. The configurations of Benchmarking Universal Single-Copy Ortholog (BUSCO), k-mer, and repeat clustering point to the presence of at least two sub-genomes one with 30 and the other with 35 chromosomes, indicating the allopolyploid nature of the okra genome. Over 130 000 putative genes, derived from mapped IsoSeq data and transcriptome data from public okra accessions, exhibit a low genetic diversity of one single nucleotide polymorphisms per 2.1 kbp. The genes are predominantly located at the distal chromosome ends, declining toward central scaffold domains. Long terminal repeat retrotransposons prevail in central domains, consistent with the observed pericentromeric heterochromatin and distal euchromatin. Disparities in paralogous gene counts suggest potential sub-genome differentiation implying possible sub-genome dominance. Amino acid query sequences of putative genes facilitated phenol biosynthesis pathway annotation. Comparison with manually curated reference KEGG pathways from related Malvaceae species reveals the genetic basis for putative enzyme coding genes that likely enable metabolic reactions involved in the biosynthesis of dietary and therapeutic compounds in okra.

The genome of Lactuca saligna, a wild relative of lettuce, provides insight into non-host resistance to the downy mildew Bremia lactucae.

Lactuca saligna L. is a wild relative of cultivated lettuce (Lactuca sativa L.), with which it is partially interfertile. Hybrid progeny suffer from hybrid incompatibility (HI), resulting in reduced fertility and distorted transmission ratios. Lactuca saligna displays broad-spectrum resistance against lettuce downy mildew caused by Bremia lactucae Regel and is considered a non-host species. This phenomenon of resistance in L. saligna is called non-host resistance (NHR). One possible mechanism behind this NHR is through the plant-pathogen interaction triggered by pathogen recognition receptors, including nucleotide-binding leucine-rich repeat (NLR) proteins and receptor-like kinases (RLKs). We report a chromosome-level genome assembly of L. saligna (accession CGN05327), leading to the identification of two large paracentric inversions (>50 Mb) between L. saligna and L. sativa. Genome-wide searches delineated the major resistance clusters as regions enriched in NLRs and RLKs. Three of the enriched regions co-locate with previously identified NHR intervals. RNA-seq analysis of Bremia-infected lettuce identified several differentially expressed RLKs in NHR regions. Three tandem wall-associated kinase-encoding genes (WAKs) in the NHR8 interval display particularly high expression changes at an early stage of infection. We propose RLKs as strong candidates for determinants of the NHR phenotype of L. saligna.

The genome of Chenopodium quinoa.

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

Corrigendum: The genome of Chenopodium quinoa.

This corrects the article DOI: 10.1038/nature21370.

Meiotic recombination profiling of interspecific hybrid F1 tomato pollen by linked read sequencing.

Genome wide screening of pooled pollen samples from a single interspecific F1 hybrid obtained from a cross between tomato, Solanum lycopersicum and its wild relative, Solanum pimpinellifolium using linked read sequencing of the haploid nuclei, allowed profiling of the crossover (CO) and gene conversion (GC) landscape. We observed a striking overlap between cold regions of CO in the male gametes and our previously established F6 recombinant inbred lines (RILs) population. COs were overrepresented in non-coding regions in the gene promoter and 5'UTR regions of genes. Poly-A/T and AT rich motifs were found enriched in 1 kb promoter regions flanking the CO sites. Non-crossover associated allelic and ectopic GCs were detected in most chromosomes, confirming that besides CO, GC represents also a source for genetic diversity and genome plasticity in tomato. Furthermore, we identified processed break junctions pointing at the involvement of both homology directed and non-homology directed repair pathways, suggesting a recombination machinery in tomato that is more complex than currently anticipated.

The role of strigolactones in P deficiency induced transcriptional changes in tomato roots.

BACKGROUND: Phosphorus (P) is an essential macronutrient for plant growth and development. Upon P shortage, plant responds with massive reprogramming of transcription, the Phosphate Starvation Response (PSR). In parallel, the production of strigolactones (SLs)-a class of plant hormones that regulates plant development and rhizosphere signaling molecules-increases. It is unclear, however, what the functional link is between these two processes. In this study, using tomato as a model, RNAseq was used to evaluate the time-resolved changes in gene expression in the roots upon P starvation and, using a tomato CAROTENOID CLEAVAGE DIOXYGENASES 8 (CCD8) RNAi line, what the role of SLs is in this. RESULTS: Gene ontology (GO)-term enrichment and KEGG analysis of the genes regulated by P starvation and P replenishment revealed that metabolism is an important component of the P starvation response that is aimed at P homeostasis, with large changes occurring in glyco-and galactolipid and carbohydrate metabolism, biosynthesis of secondary metabolites, including terpenoids and polyketides, glycan biosynthesis and metabolism, and amino acid metabolism. In the CCD8 RNAi line about 96% of the PSR genes was less affected than in wild-type (WT) tomato. For example, phospholipid biosynthesis was suppressed by P starvation, while the degradation of phospholipids and biosynthesis of substitute lipids such as sulfolipids and galactolipids were induced by P starvation. Around two thirds of the corresponding transcriptional changes depend on the presence of SLs. Other biosynthesis pathways are also reprogrammed under P starvation, such as phenylpropanoid and carotenoid biosynthesis, pantothenate and CoA, lysine and alkaloids, and this also partially depends on SLs. Additionally, some plant hormone biosynthetic pathways were affected by P starvation and also here, SLs are required for many of the changes (more than two thirds for Gibberellins and around one third for Abscisic acid) in the gene expression. CONCLUSIONS: Our analysis shows that SLs are not just the end product of the PSR in plants (the signals secreted by plants into the rhizosphere), but also play a major role in the regulation of the PSR (as plant hormone).

Comparative genomics of the nonlegume Parasponia reveals insights into evolution of nitrogen-fixing rhizobium symbioses.

Nodules harboring nitrogen-fixing rhizobia are a well-known trait of legumes, but nodules also occur in other plant lineages, with rhizobia or the actinomycete Frankia as microsymbiont. It is generally assumed that nodulation evolved independently multiple times. However, molecular-genetic support for this hypothesis is lacking, as the genetic changes underlying nodule evolution remain elusive. We conducted genetic and comparative genomics studies by using Parasponia species (Cannabaceae), the only nonlegumes that can establish nitrogen-fixing nodules with rhizobium. Intergeneric crosses between Parasponia andersonii and its nonnodulating relative Trema tomentosa demonstrated that nodule organogenesis, but not intracellular infection, is a dominant genetic trait. Comparative transcriptomics of P. andersonii and the legume Medicago truncatula revealed utilization of at least 290 orthologous symbiosis genes in nodules. Among these are key genes that, in legumes, are essential for nodulation, including NODULE INCEPTION (NIN) and RHIZOBIUM-DIRECTED POLAR GROWTH (RPG). Comparative analysis of genomes from three Parasponia species and related nonnodulating plant species show evidence of parallel loss in nonnodulating species of putative orthologs of NIN, RPG, and NOD FACTOR PERCEPTION Parallel loss of these symbiosis genes indicates that these nonnodulating lineages lost the potential to nodulate. Taken together, our results challenge the view that nodulation evolved in parallel and raises the possibility that nodulation originated ∼100 Mya in a common ancestor of all nodulating plant species, but was subsequently lost in many descendant lineages. This will have profound implications for translational approaches aimed at engineering nitrogen-fixing nodules in crop plants.

Novel routes towards bioplastics from plants: elucidation of the methylperillate biosynthesis pathway from Salvia dorisiana trichomes.

Plants produce a large variety of highly functionalized terpenoids. Functional groups such as partially unsaturated rings and carboxyl groups provide handles to use these compounds as feedstock for biobased commodity chemicals. For instance, methylperillate, a monoterpenoid found in Salvia dorisiana, may be used for this purpose, as it carries both an unsaturated ring and a methylated carboxyl group. The biosynthetic pathway of methylperillate in plants is still unclear. In this work, we identified glandular trichomes from S. dorisiana as the location of biosynthesis and storage of methylperillate. mRNA from purified trichomes was used to identify four genes that can encode the pathway from geranyl diphosphate towards methylperillate. This pathway includes a (-)-limonene synthase (SdLS), a limonene 7-hydroxylase (SdL7H, CYP71A76), and a perillyl alcohol dehydrogenase (SdPOHDH). We also identified a terpene acid methyltransferase, perillic acid O-methyltransferase (SdPAOMT), with homology to salicylic acid OMTs. Transient expression in Nicotiana benthamiana of these four genes, in combination with a geranyl diphosphate synthase to boost precursor formation, resulted in production of methylperillate. This demonstrates the potential of these enzymes for metabolic engineering of a feedstock for biobased commodity chemicals.

Host- and stage-dependent secretome of the arbuscular mycorrhizal fungus Rhizophagus irregularis.

Arbuscular mycorrhizal fungi form the most wide-spread endosymbiosis with plants. There is very little host specificity in this interaction, however host preferences as well as varying symbiotic efficiencies have been observed. We hypothesize that secreted proteins (SPs) may act as fungal effectors to control symbiotic efficiency in a host-dependent manner. Therefore, we studied whether arbuscular mycorrhizal (AM) fungi adjust their secretome in a host- and stage-dependent manner to contribute to their extremely wide host range. We investigated the expression of SP-encoding genes of Rhizophagus irregularis in three evolutionary distantly related plant species, Medicago truncatula, Nicotiana benthamiana and Allium schoenoprasum. In addition we used laser microdissection in combination with RNA-seq to study SP expression at different stages of the interaction in Medicago. Our data indicate that most expressed SPs show roughly equal expression levels in the interaction with all three host plants. In addition, a subset shows significant differential expression depending on the host plant. Furthermore, SP expression is controlled locally in the hyphal network in response to host-dependent cues. Overall, this study presents a comprehensive analysis of the R. irregularis secretome, which now offers a solid basis to direct functional studies on the role of fungal SPs in AM symbiosis.

Correcting palindromes in long reads after whole-genome amplification.

BACKGROUND: Next-generation sequencing requires sufficient DNA to be available. If limited, whole-genome amplification is applied to generate additional amounts of DNA. Such amplification often results in many chimeric DNA fragments, in particular artificial palindromic sequences, which limit the usefulness of long sequencing reads. RESULTS: Here, we present Pacasus, a tool for correcting such errors. Two datasets show that it markedly improves read mapping and de novo assembly, yielding results similar to these that would be obtained with non-amplified DNA. CONCLUSIONS: With Pacasus long-read technologies become available for sequencing targets with very small amounts of DNA, such as single cells or even single chromosomes.

Lost in diversity: the interactions between soil-borne fungi, biodiversity and plant productivity.

There is consensus that plant species richness enhances plant productivity within natural grasslands, but the underlying drivers remain debated. Recently, differential accumulation of soil-borne fungal pathogens across the plant diversity gradient has been proposed as a cause of this pattern. However, the below-ground environment has generally been treated as a 'black box' in biodiversity experiments, leaving these fungi unidentified. Using next generation sequencing and pathogenicity assays, we analysed the community composition of root-associated fungi from a biodiversity experiment to examine if evidence exists for host specificity and negative density dependence in the interplay between soil-borne fungi, plant diversity and productivity. Plant species were colonised by distinct (pathogenic) fungal communities and isolated fungal species showed negative, species-specific effects on plant growth. Moreover, 57% of the pathogenic fungal operational taxonomic units (OTUs) recorded in plant monocultures were not detected in eight plant species plots, suggesting a loss of pathogenic OTUs with plant diversity. Our work provides strong evidence for host specificity and negative density-dependent effects of root-associated fungi on plant species in grasslands. Our work substantiates the hypothesis that fungal root pathogens are an important driver of biodiversity-ecosystem functioning relationships.

Re-sequencing transgenic plants revealed rearrangements at T-DNA inserts, and integration of a short T-DNA fragment, but no increase of small mutations elsewhere.

KEY MESSAGE: Transformation resulted in deletions and translocations at T-DNA inserts, but not in genome-wide small mutations. A tiny T-DNA splinter was detected that probably would remain undetected by conventional techniques. We investigated to which extent Agrobacterium tumefaciens-mediated transformation is mutagenic, on top of inserting T-DNA. To prevent mutations due to in vitro propagation, we applied floral dip transformation of Arabidopsis thaliana. We re-sequenced the genomes of five primary transformants, and compared these to genomic sequences derived from a pool of four wild-type plants. By genome-wide comparisons, we identified ten small mutations in the genomes of the five transgenic plants, not correlated to the positions or number of T-DNA inserts. This mutation frequency is within the range of spontaneous mutations occurring during seed propagation in A. thaliana, as determined earlier. In addition, we detected small as well as large deletions specifically at the T-DNA insert sites. Furthermore, we detected partial T-DNA inserts, one of these a tiny 50-bp fragment originating from a central part of the T-DNA construct used, inserted into the plant genome without flanking other T-DNA. Because of its small size, we named this fragment a T-DNA splinter. As far as we know this is the first report of such a small T-DNA fragment insert in absence of any T-DNA border sequence. Finally, we found evidence for translocations from other chromosomes, flanking T-DNA inserts. In this study, we showed that next-generation sequencing (NGS) is a highly sensitive approach to detect T-DNA inserts in transgenic plants.

Gene and transposable element methylation in great tit (Parus major) brain and blood.

BACKGROUND: Studies on vertebrate DNA methylomes have revealed a regulatory role of tissue specific DNA methylation in relation to gene expression. However, it is not well known how tissue-specific methylation varies between different functional and structural components of genes and genomes. Using whole-genome bisulfite sequencing data we here describe both CpG and non-CpG methylation profiles of whole blood and brain tissue in relation to gene features, CpG-islands (CGIs), transposable elements (TE), and their functional roles in an ecological model species, the great tit (Parus major). RESULTS: We show that hypomethylation at the transcription start site (TSS) is enriched in genes with functional classes that relate directly to processes specific to each tissue type. We find that 6877 (~21 %) of the CGIs are differentially methylated between blood and brain, of which 1186 and 2055 are annotated to promoter and intragenic regions, respectively. We observe that CGI methylation in promoter regions is more conserved between tissues compared to CGI methylation in intra and inter-genic regions. Differentially methylated CGIs in promoter and intragenic regions are overrepresented in genomic loci linked to development, suggesting a distinct role for CGI methylation in regulating expression during development. Additionally, we find significant non-CpG methylation in brain but not in blood with a strong preference for methylation at CpA dinucleotide sites. Finally, CpG hypermethylation of TEs is significantly stronger in brain compared to blood, but does not correlate with TE activity. Surprisingly, TEs showed significant hypomethylation in non-CpG contexts which was negatively correlated with TE expression. CONCLUSION: The discovery that TSS methylation levels are directly linked to functional classes related to each tissue provides new insights in the regulatory role of DNA-methylation patterns. The dominant sequence motifs for brain non-CpG methylation, similar to those found in mammals, suggests that a conserved non-CpG regulatory mechanism was already present in the amniote ancestor. The negative correlation between brain non-CpG methylation and TE activity (not found for CpG methylation) suggests that non-CpG is the dominant regulatory form of methylation in TE silencing.

Living apart together: crosstalk between the core and supernumerary genomes in a fungal plant pathogen.

BACKGROUND: Eukaryotes display remarkable genome plasticity, which can include supernumerary chromosomes that differ markedly from the core chromosomes. Despite the widespread occurrence of supernumerary chromosomes in fungi, their origin, relation to the core genome and the reason for their divergent characteristics are still largely unknown. The complexity of genome assembly due to the presence of repetitive DNA partially accounts for this. RESULTS: Here we use single-molecule real-time (SMRT) sequencing to assemble the genome of a prominent fungal wheat pathogen, Fusarium poae, including at least one supernumerary chromosome. The core genome contains limited transposable elements (TEs) and no gene duplications, while the supernumerary genome holds up to 25 % TEs and multiple gene duplications. The core genome shows all hallmarks of repeat-induced point mutation (RIP), a defense mechanism against TEs, specific for fungi. The absence of RIP on the supernumerary genome accounts for the differences between the two (sub)genomes, and results in a functional crosstalk between them. The supernumerary genome is a reservoir for TEs that migrate to the core genome, and even large blocks of supernumerary sequence (>200 kb) have recently translocated to the core. Vice versa, the supernumerary genome acts as a refuge for genes that are duplicated from the core genome. CONCLUSIONS: For the first time, a mechanism was determined that explains the differences that exist between the core and supernumerary genome in fungi. Different biology rather than origin was shown to be responsible. A "living apart together" crosstalk exists between the core and supernumerary genome, accelerating chromosomal and organismal evolution.

A detailed analysis of the recombination landscape of the button mushroom Agaricus bisporus var. bisporus.

The button mushroom (Agaricus bisporus) is one of the world's most cultivated mushroom species, but in spite of its economic importance generation of new cultivars by outbreeding is exceptional. Previous genetic analyses of the white bisporus variety, including all cultivars and most wild isolates revealed that crossing over frequencies are low, which might explain the lack of introducing novel traits into existing cultivars. By generating two high quality whole genome sequence assemblies (one de novo and the other by improving the existing reference genome) of the first commercial white hybrid Horst U1, a detailed study of the crossover (CO) landscape was initiated. Using a set of 626 SNPs in a haploid offspring of 139 single spore isolates and whole genome sequencing on a limited number of homo- and heterokaryotic single spore isolates, we precisely mapped all COs showing that they are almost exclusively restricted to regions of about 100kb at the chromosome ends. Most basidia of A. bisporus var. bisporus produce two spores and pair preferentially via non-sister nuclei. Combined with the COs restricted to the chromosome ends, these spores retain most of the heterozygosity of the parent thus explaining how present-day white cultivars are genetically so close to the first hybrid marketed in 1980. To our knowledge this is the first example of an organism which displays such specific CO landscape.

Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing.

We explored genetic variation by sequencing a selection of 84 tomato accessions and related wild species representative of the Lycopersicon, Arcanum, Eriopersicon and Neolycopersicon groups, which has yielded a huge amount of precious data on sequence diversity in the tomato clade. Three new reference genomes were reconstructed to support our comparative genome analyses. Comparative sequence alignment revealed group-, species- and accession-specific polymorphisms, explaining characteristic fruit traits and growth habits in the various cultivars. Using gene models from the annotated Heinz 1706 reference genome, we observed differences in the ratio between non-synonymous and synonymous SNPs (dN/dS) in fruit diversification and plant growth genes compared to a random set of genes, indicating positive selection and differences in selection pressure between crop accessions and wild species. In wild species, the number of single-nucleotide polymorphisms (SNPs) exceeds 10 million, i.e. 20-fold higher than found in most of the crop accessions, indicating dramatic genetic erosion of crop and heirloom tomatoes. In addition, the highest levels of heterozygosity were found for allogamous self-incompatible wild species, while facultative and autogamous self-compatible species display a lower heterozygosity level. Using whole-genome SNP information for maximum-likelihood analysis, we achieved complete tree resolution, whereas maximum-likelihood trees based on SNPs from ten fruit and growth genes show incomplete resolution for the crop accessions, partly due to the effect of heterozygous SNPs. Finally, results suggest that phylogenetic relationships are correlated with habitat, indicating the occurrence of geographical races within these groups, which is of practical importance for Solanum genome evolution studies.

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.