De novo transcriptome of Ischnura elegans provides insights into sensory biology, colour and vision genes

Increasing complexity of opsin expression across stomatopod development

Stomatopods are well studied for their unique visual systems, which can consist of up to 16 different photoreceptor types and 33 opsin proteins expressed in the adults of some species. The light-sensing abilities of larval stomatopods are comparatively less well understood with limited information about the opsin repertoire of these early-life stages. Early work has suggested that larval stomatopods may not possess the extensive light detection abilities found in their adult counterparts. However, recent studies have shown that these larvae may have more complex photosensory systems than previously thought. To examine this idea at the molecular level, we characterized the expression of putative light-absorbing opsins across developmental stages, from embryo to adult, in the stomatopod species Pullosquilla thomassini using transcriptomic methods with a special focus on ecological and physiological transition periods. Opsin expression during the transition from the larval to the adult stage was further characterized in the species Gonodactylaceus falcatus. Opsin transcripts from short, middle, and long wavelength-sensitive clades were found in both species, and analysis of spectral tuning sites suggested differences in absorbance within these clades. This is the first study to document the changes in opsin repertoire across development in stomatopods, providing novel evidence for light detection across the visual spectrum in larvae.

Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, non-model insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression.

Genome assembly, sex-biased gene expression and dosage compensation in the damselfly Ischnura elegans

The evolution of sex chromosomes, and patterns of sex-biased gene expression and dosage compensation, are poorly known among early winged insects such as odonates. We assembled and annotated the genome of Ischnura elegans (blue-tailed damselfly), which, like other odonates, has a male-hemigametic sex-determining system (X0 males, XX females). By identifying X-linked genes in I. elegans and their orthologs in other insect genomes, we found homologies between the X chromosome in odonates and chromosomes of other orders, including the X chromosome in Coleoptera. Next, we showed balanced expression of X-linked genes between sexes in adult I. elegans, i.e. evidence of dosage compensation. Finally, among the genes in the sex-determining pathway only fruitless was found to be X-linked, while only doublesex showed sex-biased expression. This study reveals partly conserved sex chromosome synteny and independent evolution of dosage compensation among insect orders separated by several hundred million years of evolutionary history.

An integrated transcriptomic and proteomic approach to identify the main Torymus sinensis venom components

During oviposition, ectoparasitoid wasps not only inject their eggs but also a complex mixture of proteins and peptides (venom) in order to regulate the host physiology to benefit their progeny. Although several endoparasitoid venom proteins have been identified, little is known about the components of ectoparasitoid venom. To characterize the protein composition of Torymus sinensis Kamijo (Hymenoptera: Torymidae) venom, we used an integrated transcriptomic and proteomic approach and identified 143 venom proteins. Moreover, focusing on venom gland transcriptome, we selected additional 52 transcripts encoding putative venom proteins. As in other parasitoid venoms, hydrolases, including proteases, phosphatases, esterases, and nucleases, constitute the most abundant families in T. sinensis venom, followed by protease inhibitors. These proteins are potentially involved in the complex parasitic syndrome, with different effects on the immune system, physiological processes and development of the host, and contribute to provide nutrients to the parasitoid progeny. Although additional in vivo studies are needed, initial findings offer important information about venom factors and their putative host effects, which are essential to ensure the success of parasitism.

Pigmentation and color pattern diversity in Odonata

The order Odonata (dragonflies and damselflies) comprises diurnal insects with well-developed vision, showing diverse colors in adult wings and bodies. It is one of the most ancestral winged insect groups. Because Odonata species use visual cues to recognize each other, color patterns have been investigated from ecological and evolutionary viewpoints. Here we review the recent progress on molecular mechanisms of pigmentation, especially focused on light-blue coloration. Results from histology and pigment analysis showed that ommochrome pigments on the proximal layer and pteridine pigments on the distal layer of the epidermis are essential for light-blue coloration. We also summarize genes involved in the biosynthesis of three major insect pigments conserved across insects and discuss that gene-functional analysis deserves future studies.

Exceptional diversity of opsin expression patterns in Neogonodactylus oerstedii (Stomatopoda) retinas

Stomatopod crustaceans possess some of the most complex animal visual systems, including at least 16 spectrally distinct types of photoreceptive units (e.g., assemblages of photoreceptor cells). Here we fully characterize the set of opsin genes expressed in retinal tissues and determine expression patterns of each in the stomatopod Neogonodactylus oerstedii Using a combination of transcriptome and RACE sequencing, we identified 33 opsin transcripts expressed in each N. oerstedii eye, which are predicted to form 20 long-wavelength-sensitive, 10 middle-wavelength-sensitive, and three UV-sensitive visual pigments. Observed expression patterns of these 33 transcripts were highly unusual in five respects: 1) All long-wavelength and short/middle-wavelength photoreceptive units expressed multiple opsins, while UV photoreceptor cells expressed single opsins; 2) most of the long-wavelength photoreceptive units expressed at least one middle-wavelength-sensitive opsin transcript; 3) the photoreceptors involved in spatial, motion, and polarization vision expressed more transcripts than those involved in color vision; 4) there is a unique opsin transcript that is expressed in all eight of the photoreceptive units devoted to color vision; and 5) expression patterns in the peripheral hemispheres of the eyes suggest visual specializations not previously recognized in stomatopods. Elucidating the expression patterns of all opsin transcripts expressed in the N. oerstedii retina reveals the potential for previously undocumented functional diversity in the already complex stomatopod eye and is a first step toward understanding the functional significance of the unusual abundance of opsins found in many arthropod species' visual systems.

Retracted: The evolutionary history of colour polymorphism in Ischnura damselflies

The above article from Journal of Evolutionary Biology, published online on 24 May 2018 in Wiley Online Library (http://wileyonlinelibrary.com), has been retracted on the request of the authors and with the agreement of the Journal's Editor in Chief Wolf Blanckenhorn and John Wiley & Sons, following disagreement on potential corrections to the article after publication. The decision to retract followed significant issues with the methods and analyses of the manuscript that were originally not uncovered during peer-review, but which were subsequently brought to the Journal's attention following publication of the Article on Early View. [Correction added on 2 July 2021, after first online publication: retraction statement has been modified.].

Linking crop traits to transcriptome differences in a progeny population of tetraploid potato

BACKGROUND

Potato is the third most consumed crop in the world. Breeding for traits such as yield, product quality and pathogen resistance are main priorities. Identifying molecular signatures of these and other important traits is important in future breeding efforts. In this study, a progeny population from a cross between a breeding line, SW93-1015, and a cultivar, Désirée, was studied by trait analysis and RNA-seq in order to develop understanding of segregating traits at the molecular level and identify transcripts with expressional correlation to these traits. Transcript markers with predictive value for field performance applicable under controlled environments would be of great value for plant breeding.

RESULTS

A total of 34 progeny lines from SW93-1015 and Désirée were phenotyped for 17 different traits in a field in Nordic climate conditions and controlled climate settings. A master transcriptome was constructed with all 34 progeny lines and the parents through a de novo assembly of RNA-seq reads. Gene expression data obtained in a controlled environment from the 34 lines was correlated to traits by different similarity indices, including Pearson and Spearman, as well as DUO, which calculates the co-occurrence between high and low values for gene expression and trait. Our study linked transcripts to traits such as yield, growth rate, high laying tubers, late and tuber blight, tuber greening and early flowering. We found several transcripts associated to late blight resistance and transcripts encoding receptors were associated to Dickeya solani susceptibility. Transcript levels of a UBX-domain protein was negatively associated to yield and a GLABRA2 expression modulator was negatively associated to growth rate.

CONCLUSION

In our study, we identify 100's of transcripts, putatively linked based on expression with 17 traits of potato, representing both well-known and novel associations. This approach can be used to link the transcriptome to traits. We explore the possibility of associating the level of transcript expression from controlled, optimal environments to traits in a progeny population with different methods introducing the application of DUO for the first time on transcriptome data. We verify the expression pattern for five of the putative transcript markers in another progeny population.

Charting oat (Avena sativa) embryo and endosperm transcription factor expression reveals differential expression of potential importance for seed development

Uniquely, oat, among cereals, accumulates an appreciable amount of oil in the endosperm together with starch. Oat is also recognized for its soluble fibers in the form of β-glucans. Despite high and increasing interest in oat yield and quality, the genetic and molecular understanding of oat grain development is still very limited. Transcription factors (TFs) are important regulatory components for plant development, product quality and yield. This study aimed to develop a workflow to determine seed tissue specificity of transcripts encoding transcription factors to reveal differential expression of potential importance for storage compound deposition and quality characters in oat. We created a workflow through the de novo assembly of sequenced seed endosperm and embryo, and publicly available oat seed RNAseq dataset, later followed by TF identification. RNAseq data were assembled into 33,878 transcripts with approximately 90% completeness. A total of 3875 putative TF encoding transcripts were identified from the oat hybrid assemblies. Members of the B3, bHLH, bZIP, C3H, ERF, NAC, MYB and WRKY families were the most abundant TF transcripts. A total of 514 transcripts which were differentially expressed between embryo and endosperm were identified with a threshold of 16-fold expression difference. Among those, 36 TF transcript homologs, belonging to 7 TF families, could be identified through similarity search in wheat embryo and endosperm EST libraries of NCBI Unigene database, and almost all the closest homologs were specifically expressed in seed when explored in WheatExp database. We verified our findings by cloning, sequencing and finally confirming differential expression of two TF encoding transcripts in oat seed embryo and endosperm. The developed workflow for identifying tissue-specific transcription factors allows further functional characterization of specific genes to increase our understanding of grain filling and quality.

Pterin-pigmented nanospheres create the colours of the polymorphic damselfly Ischnura elegans

Animal colours commonly act as signals for mates or predators. In many damselfly species, both sexes go through a developmental colour change as adults, and females often show colour polymorphism, which may have a function in mate choice, avoidance of mating harassment and camouflage. In the blue-tailed damselfly, Ischnura elegans, young males are bright green and turn blue as they reach maturity. Females are red ( rufescens) or violet ( violacea) as immatures and, when mature, either mimic the blue colour of the males ( androchrome), or acquire an inconspicuous olive-green ( infuscans) or olive-brown ( obsoleta). The genetic basis of these differences is still unknown. Here, we quantify the colour development of all morphs of I. elegans and investigate colour formation by combining anatomical data and reflectance spectra with optical finite-difference time-domain simulations. While the coloration primarily arises from a disordered assembly of nanospheres in the epidermis, morph-dependent changes result from adjustments in the composition of pterin pigments within the nanospheres, and from associated shifts in optical density. Other pigments fine-tune hue and brilliance by absorbing stray light. These mechanisms produce an impressive palette of colours and offer guidance for genetic studies on the evolution of colour polymorphism and visual communication.

Laboratory Rearing System for Ischnura senegalensis (Insecta: Odonata) Enables Detailed Description of Larval Development and Morphogenesis in Dragonfly

In an attempt to establish an experimental dragonfly model, we developed a laboratory rearing system for the blue-tailed damselfly, Ischnura senegalensis. Adoption of multi-well plastic plates as rearing containers enabled mass-rearing of isolated larvae without cannibalism and convenient microscopic monitoring of individual larvae. Feeding Artemia brine shrimps to younger larvae and Tubifex worms for older larvae resulted in low mortality, synchronized ecdysis, and normal development of the larvae. We continuously monitored the development of 118 larvae every day, of which 49 individuals (41.5%) reached adulthood. The adult insects were fed with Drosophila flies in wet plastic cages, attained reproductive maturity in a week, copulated, laid fertilized eggs, and produced progeny. The final larval instar varied from 9^th to 12^th, with the 11^th instar (56.5%) and the 12^th instar (24.2%) constituting the majority. From the 1^st instar to the penultimate instar, the duration of each instar was relatively short, mainly ranging from three to 11 days. Afterwards, the duration of each instar was prolonged, reaching 7-25 days for the penultimate instar and 14-28 days for the final instar. Some larvae of final, penultimate and younger instars were subjected to continuous and close morphological examinations, which enabled developmental staging of larvae based on size, shape, and angle of compound eyes and other morphological traits. This laboratory rearing system may facilitate the understanding of physiological, biochemical, and molecular mechanisms underlying metamorphosis, hormonal control, morphogenesis, body color polymorphism, and other biological features of dragonflies.

Ommochromes in invertebrates: biochemistry and cell biology

Ommochromes are widely occurring coloured molecules of invertebrates, arising from tryptophan catabolism through the so-called Tryptophan → Ommochrome pathway. They are mainly known to mediate compound eye vision, as well as reversible and irreversible colour patterning. Ommochromes might also be involved in cell homeostasis by detoxifying free tryptophan and buffering oxidative stress. These biological functions are directly linked to their unique chromophore, the phenoxazine/phenothiazine system. The most recent reviews on ommochrome biochemistry were published more than 30 years ago, since when new results on the enzymes of the ommochrome pathway, on ommochrome photochemistry as well as on their antiradical capacities have been obtained. Ommochromasomes are the organelles where ommochromes are synthesised and stored. Hence, they play an important role in mediating ommochrome functions. Ommochromasomes are part of the lysosome-related organelles (LROs) family, which includes other pigmented organelles such as vertebrate melanosomes. Ommochromasomes are unique because they are the only LRO for which a recycling process during reversible colour change has been described. Herein, we provide an update on ommochrome biochemistry, photoreactivity and antiradical capacities to explain their diversity and behaviour both in vivo and in vitro. We also highlight new biochemical techniques, such as quantum chemistry, metabolomics and crystallography, which could lead to major advances in their chemical and functional characterisation. We then focus on ommochromasome structure and formation by drawing parallels with the well-characterised melanosomes of vertebrates. The biochemical, genetic, cellular and microscopic tools that have been applied to melanosomes should provide important information on the ommochromasome life cycle. We propose LRO-based models for ommochromasome biogenesis and recycling that could be tested in the future. Using the context of insect compound eyes, we finally emphasise the importance of an integrated approach in understanding the biological functions of ommochromes.

Transcriptomic analysis of male and female Schistosoma mekongi adult worms

Background

Schistosoma mekongi is one of five major causative agents of human schistosomiasis and is endemic to communities along the Mekong River in southern Lao People’s Democratic Republic (Laos) and northern Cambodia. Sporadic cases of schistosomiasis have been reported in travelers and immigrants who have visited endemic areas. Schistosoma mekongi biology and molecular biology is poorly understood, and few S. mekongi gene and transcript sequences are available in public databases.

Results

Transcriptome sequencing (RNA-Seq) of male and female S. mekongi adult worms (a total of three biological replicates for each sex) were analyzed and the results demonstrated that approximately 304.9 and 363.3 million high-quality clean reads with quality Q30 (> 90%) were obtained from male and female adult worms, respectively. A total of 119,604 contigs were assembled with an average length of 1273 nt and an N50 of 2017 nt. From the contigs, 20,798 annotated protein sequences and 48,256 annotated transcript sequences were obtained using BLASTP and BLASTX searches against the UniProt Trematoda database. A total of 4658 and 3509 transcripts were predominantly expressed in male and female worms, respectively. Male-biased transcripts were mostly involved in structural organization while female-biased transcripts were typically involved in cell differentiation and egg production. Interestingly, pathway enrichment analysis suggested that genes involved in the phosphatidylinositol signaling pathway may play important roles in the cellular processes and reproductive systems of S. mekongi worms.

Conclusions

We present comparative transcriptomic analyses of male and female S. mekongi adult worms, which provide a global view of the S. mekongi transcriptome as well as insights into differentially-expressed genes associated with each sex. This work provides valuable information and sequence resources for future studies of gene function and for ongoing whole genome sequencing efforts in S. mekongi.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3086-z) contains supplementary material, which is available to authorized users.

Molecular and ecological signatures of an expanding hybrid zone

Many species are currently changing their distributions and subsequently form sympatric zones with hybridization between formerly allopatric species as one possible consequence. The damselfly Ischnura elegans has recently expanded south into the range of its ecologically and morphologically similar sister species Ischnura graellsii. Molecular work shows ongoing introgression between these species, but the extent to which this species mixing is modulated by ecological niche use is not known. Here, we (1) conduct a detailed population genetic analysis based on molecular markers and (2) model the ecological niche use of both species in allopatric and sympatric regions. Population genetic analyses showed chronic introgression between I. elegans and I. graellsii across a wide part of Spain, and admixture analysis corroborated this, showing that the majority of I. elegans from the sympatric zone could not be assigned to either the I. elegans or I. graellsii species cluster. Niche modeling demonstrated that I. elegans has modified its environmental niche following hybridization and genetic introgression with I. graellsii, making niche space of introgressed I. elegans populations more similar to I. graellsii. Taken together, this corroborates the view that adaptive introgression has moved genes from I. graellsii into I. elegans and that this process is enabling Spanish I. elegans to occupy a novel niche, further facilitating its expansion. Our results add to the growing evidence that hybridization can play an important and creative role in the adaptive evolution of animals.

Transcriptional profiling of long non-coding RNAs in mantle of Crassostrea gigas and their association with shell pigmentation

Long non-coding RNAs (lncRNAs) play crucial roles in diverse biological processes and have drawn extensive attention in the past few years. However, lncRNAs remain poorly understood about expression and roles in Crassostrea gigas, a potential model organism for marine molluscan studies. Here, we systematically identified lncRNAs in the mantles of C. gigas from four full-sib families characterized by white, black, golden, and partially pigmented shell. Using poly(A)-independent and strand-specific RNA-seq, a total of 441,205,852 clean reads and 12,243 lncRNA transcripts were obtained. LncRNA transcripts were relatively short with few exons and low levels of expression in comparison to protein coding mRNA transcripts. A total of 427 lncRNAs and 349 mRNAs were identified to differentially express among six pairwise groups, mainly involving in biomineralization and pigmentation through functional enrichment. Furthermore, a total of 6 mRNAs and their cis-acting lncRNAs were predicted to involve in synthesis of melanin, carotenoid, tetrapyrrole, or ommochrome. Of them, chorion peroxidase and its cis-acting lincRNA TCONS_00951105 are implicated in playing an essential role in the melanin synthetic pathway. Our studies provided the first systematic characterization of lncRNAs catalog expressed in oyster mantle, which may facilitate understanding the molecular regulation of shell colour diversity and provide new insights into future selective breeding of C. gigas for aquaculture.

Transcriptome profiling with focus on potential key genes for wing development and evolution in Megaloprepus caerulatus, the damselfly species with the world's largest wings

The evolution, development and coloration of insect wings remains a puzzling subject in evolutionary research. In basal flying insects such as Odonata, genomic research regarding bauplan evolution is still rare. Here we focus on the world's largest odonate species-the "forest giant" Megaloprepus caerulatus, to explore its potential for looking deeper into the development and evolution of wings. A recently discovered cryptic species complex in this genus previously considered monotypic is characterized by morphological differences in wing shape and color patterns. As a first step toward understanding wing pattern divergence and pathways involved in adaptation and speciation at the genomic level, we present a transcriptome profiling of M. caerulatus using RNA-Seq and compare these data with two other odonate species. The de novo transcriptome assembly consists of 61,560 high quality transcripts and is approximately 93% complete. For almost 75% of the identified transcripts a possible function could be assigned: 48,104 transcripts had a hit to an InterPro protein family or domain, and 28,653 were mapped to a Gene Ontology term. In particular, we focused on genes related to wing development and coloration. The comparison with two other species revealed larva-specific genes and a conserved 'core' set of over 8,000 genes forming orthologous clusters with Ischnura elegans and Ladona fulva. This transcriptome may provide a first point of reference for future research in odonates addressing questions surrounding the evolution of wing development, wing coloration and their role in speciation.

Comparative transcriptomics reveal developmental turning points during embryogenesis of a hemimetabolous insect, the damselfly Ischnura elegans

Identifying transcriptional changes during embryogenesis is of crucial importance for unravelling evolutionary, molecular and cellular mechanisms that underpin patterning and morphogenesis. However, comparative studies focusing on early/embryonic stages during insect development are limited to a few taxa. Drosophila melanogaster is the paradigm for insect development, whereas comparative transcriptomic studies of embryonic stages of hemimetabolous insects are completely lacking. We reconstructed the first comparative transcriptome covering the daily embryonic developmental progression of the blue-tailed damselfly Ischnura elegans (Odonata), an ancient hemimetabolous representative. We identified a "core" set of 6,794 transcripts - shared by all embryonic stages - which are mainly involved in anatomical structure development and cellular nitrogen compound metabolic processes. We further used weighted gene co-expression network analysis to identify transcriptional changes during Odonata embryogenesis. Based on these analyses distinct clusters of transcriptional active sequences could be revealed, indicating that embryos at different development stages have their own transcriptomic profile according to the developmental events and leading to sequential reprogramming of metabolic and developmental genes. Interestingly, a major change in transcriptionally active sequences is correlated with katatrepsis (revolution) during mid-embryogenesis, a 180° rotation of the embryo within the egg and specific to hemimetabolous insects.

Genomic Features of the Damselfly Calopteryx splendens Representing a Sister Clade to Most Insect Orders

Insects comprise the most diverse and successful animal group with over one million described species that are found in almost every terrestrial and limnic habitat, with many being used as important models in genetics, ecology, and evolutionary research. Genome sequencing projects have greatly expanded the sampling of species from many insect orders, but genomic resources for species of certain insect lineages have remained relatively limited to date. To address this paucity, we sequenced the genome of the banded demoiselle, Calopteryx splendens, a damselfly (Odonata: Zygoptera) belonging to Palaeoptera, the clade containing the first winged insects. The 1.6 Gbp C. splendens draft genome assembly is one of the largest insect genomes sequenced to date and encodes a predicted set of 22,523 protein-coding genes. Comparative genomic analyses with other sequenced insects identified a relatively small repertoire of C. splendens detoxification genes, which could explain its previously noted sensitivity to habitat pollution. Intriguingly, this repertoire includes a cytochrome P450 gene not previously described in any insect genome. The C. splendens immune gene repertoire appears relatively complete and features several genes encoding novel multi-domain peptidoglycan recognition proteins. Analysis of chemosensory genes revealed the presence of both gustatory and ionotropic receptors, as well as the insect odorant receptor coreceptor gene (OrCo) and at least four partner odorant receptors (ORs). This represents the oldest known instance of a complete OrCo/OR system in insects, and provides the molecular underpinning for odonate olfaction. The C. splendens genome improves the sampling of insect lineages that diverged before the radiation of Holometabola and offers new opportunities for molecular-level evolutionary, ecological, and behavioral studies.

De novo transcriptome of the mayfly Cloeon viridulum and transcriptional signatures of Prometabola

Mayflies (Ephemeroptera) display many primitive characters and a unique type of metamorphosis (Prometabola). However, information on the genomes and transcriptomes of this insect group is limited. The RNA sequencing study presented here generated the first de novo transcriptome assembly of Cloeon viridulum (Ephemeroptera: Baetidae), and compared gene expression signatures among the young larva (YL), mature larva (ML), subimago (SI), and imago (IM) stages of this mayfly. The transcriptome, based on 88 Gb of sequence data, comprised a set of 81,185 high quality transcripts. The number of differentially expressed genes (DEGs) in YL vs. ML, ML vs. SI, and SI vs. IM, was 4,825, 1,584, and 1,278, respectively, according to the reads per kilobase of transcript per million mapped reads analysis, assuming a false discovery rate <0.05 and a fold change >2. Gene enrichment analysis revealed that these DEGs were enriched in the "chitin metabolic process", "germ cell development", "steroid hormone biosynthesis", and "cutin, suberine, and wax biosynthesis" pathways. Finally, the expression pattern of a selected group of candidate signature genes for Prometabola, including vestigial, methoprene-tolerant, wingless, and broad-complex were confirmed by quantitative real time-PCR analysis. The Q-PCR analysis of larval, subimaginal, and imaginal stages of C. viridulum suggests that the development of mayflies more closely resembles hemimetamorphosis than holometamorphosis.

Transcriptome of Pterospermum kingtungense provides implications on the mechanism underlying its rapid vegetative growth and limestone adaption

Pterospermum kingtungense C.Y.Wu ex Hsue is a typical tree species living in the relatively adverse limestone habitat. Due to its excellent wood quality and big size, it is an important timber resource which caused its endangered. We firstly provide the data resources by reporting an annotated transcriptome assembly. 203 million unique Illumina RNA-seq reads were produced with totally 50,333 transcripts, among which 48,778 transcripts were annotated. By a global comparison of homology between P. kingtungense and cacao, we identified 9,507 single copy orthologues and 990 P. kingtungense specific genes. GO enrichment analyses indicate that P. kingtungense specific genes are enriched in defense response, implying potential adaptation to limestone environment. As to cell compartment, the genes are enriched in thylakoid component. Consistently, KEGG enrichment indicates that genes are enriched in photosynthesis. In addition, we identified two genes under positive selection in P. kingtungense species. These results suggest that P. kingtungense have strong photosynthetic capacity, which related to vegetation growth. Our work provides the genomic resources of a limestone specific tree with economic importance to local society and suggests possible mechanism on its characteristics on the limestone adaption and excellent wood properties, which will be important for its conservation and sustainable utilization.

Transcriptome profiling in the damselfly Ischnura elegans identifies genes with sex-biased expression

Background

Sexual dimorphism occurs widely across the animal kingdom and has profound effects on evolutionary trajectories. Here, we investigate sex-specific gene expression in Ischnura elegans (Odonata: dragonflies and damselflies), a species with pronounced sexual differences including a female-limited colour polymorphism with two female-like gynochrome morphs and one male-mimicking, androchrome morph. Whole-organism transcriptome profiling and sex-biased gene expression analysis was conducted on adults of both sexes (pooling all females as well as separating the three morphs) to gain insights into genes and pathways potentially associated with sexual development and sexual conflict.

Results

The de novo transcriptome assembly was of high quality and completeness (54 k transcripts; 99.6% CEGMA score; 55% annotated). We identified transcripts of several relevant pathways, including transcripts involved in sex determination, hormone biosynthesis, pigmentation and innate immune signalling. A total of 1,683 genes were differentially expressed (DE) between males and all females (1,173 were female-biased; 510 male-biased). The DE genes were associated with sex-specific physiological and reproductive processes, olfaction, pigmentation (ommochrome and melanin), hormone (ecdysone) biosynthesis and innate immunity signalling pathways. Comparisons between males and each female morph category showed that the gynochromes differed more from males than the androchrome morph.

Conclusions

This is the first study to characterize sex-biased gene expression in odonates, one of the most ancient extant insect orders. Comparison between I. elegans sexes revealed expression differences in several genes related to sexual differences in behaviour and development as well as morphology. The differential expression of several olfactory genes suggests interesting sexual components in the detection of odours, pheromones and environmental volatiles. Up-regulation of pigmentation pathways in females indicates a prominent role of ommochrome pigments in the formation of the genetically controlled female colour polymorphism. Finally, the female-biased expression of several immunity genes suggests a stronger immune response in females, possibly related to the high levels of male mating harassment and recurrent matings in this species, both of which have been shown to injure females and expose them to sexually transmitted diseases and toxins contained in seminal fluids.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3334-6) contains supplementary material, which is available to authorized users.

Odonata (dragonflies and damselflies) as a bridge between ecology and evolutionary genomics

Odonata (dragonflies and damselflies) present an unparalleled insect model to integrate evolutionary genomics with ecology for the study of insect evolution. Key features of Odonata include their ancient phylogenetic position, extensive phenotypic and ecological diversity, several unique evolutionary innovations, ease of study in the wild and usefulness as bioindicators for freshwater ecosystems worldwide. In this review, we synthesize studies on the evolution, ecology and physiology of odonates, highlighting those areas where the integration of ecology with genomics would yield significant insights into the evolutionary processes that would not be gained easily by working on other animal groups. We argue that the unique features of this group combined with their complex life cycle, flight behaviour, diversity in ecological niches and their sensitivity to anthropogenic change make odonates a promising and fruitful taxon for genomics focused research. Future areas of research that deserve increased attention are also briefly outlined.

Color vision and color formation in dragonflies

Dragonflies including damselflies are colorful and large-eyed insects, which show remarkable sexual dimorphism, color transition, and color polymorphism. Recent comprehensive visual transcriptomics has unveiled an extraordinary diversity of opsin genes within the lineage of dragonflies. These opsin genes are differentially expressed between aquatic larvae and terrestrial adults, as well as between dorsal and ventral regions of adult compound eyes. Recent topics of color formation in dragonflies are also outlined. Non-iridescent blue color is caused by coherent light scattering from the quasiordered nanostructures, whereas iridescent color is produced by multilayer structures. Wrinkles or wax crystals sometimes enhances multilayer structural colors. Sex-specific and stage-specific color differences in red dragonflies is attributed to redox states of ommochrome pigments.

The inheritance of female colour polymorphism in Ischnura genei (Zygoptera: Coenagrionidae), with observations on melanism under laboratory conditions

Current research on female colour polymorphism in Ischnura damselflies suggests that a balanced fitness trade-off between morphotypes contributes to the maintenance of polymorphism inside populations. The genetic inheritance system constitutes a key factor to understand morph fluctuation and fitness. Ischnura genei, an endemic species of some Mediterranean islands, has three female colour morphs, including one androchrome (male-coloured) and two gynochromes. In this study, we reared two generations of I. genei under laboratory conditions and tested male behavioural responses to female colour morphs in the field. We recorded ontogenetic colour changes and studied morph frequency in three populations from Sardinia (Italy). Morph frequencies of laboratory crosses can be explained by a model based on an autosomal locus with three alleles and sex-restricted expression, except for one crossing of 42 families with unexpected offspring. The allelic dominance relationship was androchrome > infuscans > aurantiaca. Old individuals reared in the laboratory exhibited different levels of melanism in variable extent depending on sex and morph. Results of model presentations indicate a male preference for gynochrome females and the lack of recognition of androchromes as potential mates. Aurantiaca females were the most frequent morph in the field (63–87%). Further studies in other populations and islands are needed to understand the maintenance of this polymorphism.

Short read sequencing assembly revealed the complete mitochondrial genome of Ischnura elegans Vander Linden, 1820 (Odonata: Zygoptera)

Damselflies of the genus Ischnura emerge as organisms with high potential in ecological, evolutionary and developmental research at the base of flying insects. Ischnura elegans and Ischnura hastata are for example one of the few odonate species where a complete life cycle over generations can be reared under laboratory conditions. We here report the complete mitochondrial genome of Ischnura elegans as a valuable genomic resource for future eco-evo-devo studies at the base of flying insects. The genome has a total length of 15,962 bp and displays all typical features of Odonata (dragonflies and damselflies) mitochondrial genomes in gene content and order as well as A + T content. Start and stop codons of all protein-coding genes are consistent. Most interestingly, we found four intergenic spacer regions and a long A + T rich (control) region of 1196 bp, which is almost double the size of the close relative Ischnura pumilio. We assume that the adequate insert size and iterative mapping may be more efficient in assembling this duplicated and repetitive region.

Cloning and expression profile of ionotropic receptors in the parasitoid wasp Microplitis mediator (Hymenoptera: Braconidae)

Ionotropic receptors (IRs) mainly detect the acids and amines having great importance in many insect species, representing an ancient olfactory receptor family in insects. In the present work, we performed RNAseq of Microplitis mediator antennae and identified seventeen IRs. Full-length MmedIRs were cloned and sequenced. Phylogenetic analysis of the Hymenoptera IRs revealed that ten MmedIR genes encoded "antennal IRs" and seven encoded "divergent IRs". Among the IR25a orthologous groups, two genes, MmedIR25a.1 and MmedIR25a.2, were found in M. mediator. Gene structure analysis of MmedIR25a revealed a tandem duplication of IR25a in M. mediator. The tissue distribution and development specific expression of the MmedIR genes suggested that these genes showed a broad expression profile. Quantitative gene expression analysis showed that most of the genes are highly enriched in adult antennae, indicating the candidate chemosensory function of this family in parasitic wasps. Using immunocytochemistry, we confirmed that one co-receptor, MmedIR8a, was expressed in the olfactory sensory neurons. Our data will supply fundamental information for functional analysis of the IRs in parasitoid wasp chemoreception.

Carbon dioxide detection in adult Odonata

The present paper shows, by means of single-cell recordings, responses of antennal sensory neurons of the damselfly Ischnura elegans when stimulated by air streams at different CO2 concentrations. Unlike most insects, but similarly to termites, centipedes and ticks, Odonata possess sensory neurons strongly inhibited by CO2, with the magnitude of the off-response depending upon the CO2 concentration. The Odonata antennal sensory neurons responding to CO2 are also sensitive to airborne odors; in particular, the impulse frequency is increased by isoamylamine and decreased by heptanoic and pentanoic acid. Further behavioral investigations are necessary to assign a biological role to carbon dioxide detection in Odonata.

Assembly, Assessment, and Availability of De novo Generated Eukaryotic Transcriptomes

De novo assembly of a complete transcriptome without the need for a guiding reference genome is attractive, particularly where the cost and complexity of generating a eukaryote genome is prohibitive. The transcriptome should not however be seen as just a quick and cheap alternative to building a complete genome. Transcriptomics allows the understanding and comparison of spatial and temporal samples within an organism, and allows surveying of multiple individuals or closely related species. De novo assembly in theory allows the building of a complete transcriptome without any prior knowledge of the genome. It also allows the discovery of alternate splice forms of coding RNAs and also non-coding RNAs, which are often missed by proteomic approaches, or are incompletely annotated in genome studies. The limitations of the method are that the generation of a truly complete assembly is unlikely, and so we require some methods for the assessment of the quality and appropriateness of a generated transcriptome. Whilst no single consensus pipeline or tool is agreed as optimal, various algorithms, and easy to use software do exist making transcriptome generation a more common approach. With this expansion of data, questions still exist relating to how do we make these datasets fully discoverable, comparable and most useful to understand complex biological systems?

Morphological Characters and Transcriptome Profiles Associated with Black Skin and Red Skin in Crimson Snapper (Lutjanus erythropterus)

In this study, morphology observation and illumina sequencing were performed on two different coloration skins of crimson snapper (Lutjanus erythropterus), the black zone and the red zone. Three types of chromatophores, melanophores, iridophores and xanthophores, were organized in the skins. The main differences between the two colorations were in the amount and distribution of the three chromatophores. After comparing the two transcriptomes, 9200 unigenes with significantly different expressions (ratio change ≥ 2 and q-value ≤ 0.05) were found, of which 5972 were up-regulated in black skin and 3228 were up-regulated in red skin. Through the function annotation, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the differentially transcribed genes, we excavated a number of uncharacterized candidate pigment genes as well as found the conserved genes affecting pigmentation in crimson snapper. The patterns of expression of 14 pigment genes were confirmed by the Quantitative real-time PCR analysis between the two color skins. Overall, this study shows a global survey of the morphological characters and transcriptome analysis of the different coloration skins in crimson snapper, and provides valuable cellular and genetic information to uncover the mechanism of the formation of pigment patterns in snappers.

De novo assembly and sex-specific transcriptome profiling in the sand fly Phlebotomus perniciosus (Diptera, Phlebotominae), a major Old World vector of Leishmania infantum

Background

The phlebotomine sand fly Phlebotomus perniciosus (Diptera: Psychodidae, Phlebotominae) is a major Old World vector of the protozoan Leishmania infantum, the etiological agent of visceral and cutaneous leishmaniases in humans and dogs, a worldwide re-emerging diseases of great public health concern, affecting 101 countries. Despite the growing interest in the study of this sand fly species in the last years, the development of genomic resources has been limited so far. To increase the available sequence data for P. perniciosus and to start studying the molecular basis of the sexual differentiation in sand flies, we performed whole transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females and de novo transcriptome assembly.

Results

We assembled 55,393 high quality transcripts, of which 29,292 were unique, starting from adult whole body male and female pools. 11,736 transcripts had at least one functional annotation, including full-length low abundance salivary transcripts, 981 transcripts were classified as putative long non-coding RNAs and 244 transcripts encoded for putative novel proteins specific of the Phlebotominae sub-family. Differential expression analysis identified 8590 transcripts significantly biased between sexes. Among them, some show relaxation of selective constraints when compared to their orthologs of the New World sand fly species Lutzomyia longipalpis.

Conclusions

In this paper, we present a comprehensive transcriptome resource for the sand fly species P. perniciosus built from short-read RNA-seq and we provide insights into sex-specific gene expression at adult stage. Our analysis represents a first step towards the identification of sex-specific genes and pathways and a foundation for forthcoming investigations into this important vector species, including the study of the evolution of sex-biased genes and of the sexual differentiation in phlebotomine sand flies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2088-x) contains supplementary material, which is available to authorized users.

NMDA receptor activation upstream of methyl farnesoate signaling for short day-induced male offspring production in the water flea, Daphnia pulex

BACKGROUND

The cladoceran crustacean Daphnia pulex produces female offspring by parthenogenesis under favorable conditions, but in response to various unfavorable external stimuli, it produces male offspring (environmental sex determination: ESD). We recently established an innovative system for ESD studies using D. pulex WTN6 strain, in which the sex of the offspring can be controlled simply by changes in the photoperiod: the long-day and short-day conditions can induce female and male offspring, respectively. Taking advantage of this system, we demonstrated that de novo methyl farnesoate (MF) synthesis is necessary for male offspring production. These results indicate the key role of innate MF signaling as a conductor between external environmental stimuli and the endogenous male developmental pathway. Despite these findings, the molecular mechanisms underlying up- and downstream signaling of MF have not yet been well elucidated in D. pulex.

RESULTS

To elucidate up- and downstream events of MF signaling during sex determination processes, we compared the transcriptomes of daphnids reared under the long-day (female) condition with short-day (male) and MF-treated (male) conditions. We found that genes involved in ionotropic glutamate receptors, known to mediate the vast majority of excitatory neurotransmitting processes in various organisms, were significantly activated in daphnids by the short-day condition but not by MF treatment. Administration of specific agonists and antagonists, especially for the N-methyl-D-aspartic acid (NMDA) receptor, strongly increased or decreased, respectively, the proportion of male-producing mothers. Moreover, we also identified genes responsible for male production (e.g., protein kinase C pathway-related genes). Such genes were generally shared between the short-day reared and MF-treated daphnids.

CONCLUSIONS

We identified several candidate genes regulating ESD which strongly suggests that these genes may be essential factors for male offspring production as an upstream regulator of MF signaling in D. pulex. This study provides new insight into the fundamental mechanisms underlying how living organisms alter their phenotypes in response to various external environments.