The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus

Identification and function of microRNAs in hemipteran pests: A review

Hemiptera is one of the most significant orders of insect pests, including whiteflies, true bugs, aphids, planthoppers, psyllids, and so forth, which have led to substantial economic losses in agricultural industries and have significantly affected food yields through their ability to suck the phloem sap of crops and transmit numerous bacterial and viral pathogens. Therefore, explorations of pest-specific, eco-friendly and easy-to-adopt technologies for hemipteran pest control are urgently needed. To the best of our knowledge, microRNAs (miRNAs), which are endogenous non-coding small RNAs approximately 22 nucleotides in length, are involved in regulating gene expression via the direct recognition and binding of the 3'-untranslated region (3'-UTR) of target messenger RNAs (mRNAs) or by acting as a center of a competitive endogenous RNA (ceRNA) network at the post-transcriptional level. This review systematically outlines the characterization and functional investigation of the miRNA biogenesis pathway in hemipteran pests, such as whiteflies, true bugs, aphids and planthoppers. In addition, we explored the results of small RNA sequencing and functional observations of miRNAs in these pests, and the results suggest that the numerous miRNAs obtained and annotated via high-throughput sequencing technology and bioinformatic analyses contribute to molting development, fitness, wing polyphenism, symbiont interactions and insecticide resistance in hemipteran pests. Finally, we summarize current advances and propose a framework for future research to extend the current data and address potential limitations in the investigation and application of hemipteran miRNAs.

De novo transcriptome assembly, annotation and SSR mining data of Hellula undalis (Fabr.) (Lepidoptera: Pyralidae), the cabbage webworm

BACKGROUND

The cabbage webworm, Hellula undalis (Fabricius) (Lepidoptera: Pyralidae), is a significant pest of brassicas and other cruciferous plants in warm regions worldwide. Transcriptome analysis is valuable for investigation of molecular mechanisms underlying the insect development and reproduction. De novo assembly is particularly useful for acquiring complete transcriptome information of insect species when there is no reference genome available. In case of Hellula undalis, only 17 nucleotide records are currently available throughout NCBI nucleotide database. Genes associated with metabolic processes, general development, reproduction, defense and functional genomics were not previously predicted in the Hellula undalis at the genomic level.

METHODS & RESULTS

To address this issue, we constructed Hellula undalis transcriptome using Illumina NovaSeq6000 technology. Approximately 48 million 150 bp paired-end reads were obtained from sequencing. A total of 30,451 contigs were generated by de novo assembly of sample and were compared with the sequences in the NCBI non-redundant protein database (Nr). In total, 71 % of contigs were matched to known proteins in public databases including Nr, Gene Ontology (GO), and Cluster Orthologous Gene Database (COG), and then, contigs were mapped to 123 via functional annotation against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, we compared the ortholog gene family of the Hullula undalis, transcriptome to Spodoptera frugiperda, spodotera litura and spodoptera littoralis and found that 391 orthologous gene families are specific to Hullula undalis. A total of 1,913 potential SSRs was discovered in Hullula undalis contigs.

CONCLUSIONS

This study is the first transcriptome data for Hullula undalis. Additionally, it serves as a valuable resource for identifying target genes and developing effective and environmentally friendly strategies for pest control.

Distinct developmental mechanisms influence sexual dimorphisms in the milkweed bug Oncopeltus fasciatus

Sexual dimorphism is common in animals. The most complete model of sex determination comes from Drosophila melanogaster, where the relative dosage of autosomes and X chromosomes leads indirectly to sex-specific transcripts of doublesex (dsx). Female Dsx interacts with a mediator complex protein encoded by intersex (ix) to activate female development. In males, the transcription factor encoded by fruitless (fru) promotes male-specific behaviour. The genetics of sex determination have been examined in a small number of other insects, yet several questions remain about the plesiomorphic state. Is dsx required for female and male development? Is fru conserved in male behaviour or morphology? Are other components such as ix functionally conserved? To address these questions, we report expression and functional tests of dsx, ix and fru in the hemipteran Oncopeltus fasciatus, characterizing three sexual dimorphisms. dsx prevents ix phenotypes in all sexes and dimorphic traits in the milkweed bug. ix and fru are expressed across the body, in females and males. fru and ix also affect the genitalia of both sexes, but have effects limited to different dimorphic structures in different sexes. These results reveal roles for ix and fru distinct from other insects, and demonstrate distinct development mechanisms in different sexually dimorphic structures.

Dynamics of maternal gene expression in Rhodnius prolixus

The study of developmental processes in Rhodnius prolixus has recently advanced with the sequencing of the genome. In this work, we analyze the maternal gene expression driving oogenesis and early embryogenesis in R. prolixus. We examined the transcriptional profile of mRNAs to establish the genes expressed across the ovary, unfertilized eggs and different embryonic stages of R. prolixus until the formation of the germ band anlage (0, 12, 24, and 48 h post egg laying). We identified 81 putative maternal and ovary-related genes and validated their expression by qRT-PCR. We validate the function of the ortholog gene Bicaudal-D (Rp-BicD) by in situ hybridization and parental RNAi. Consistent with a role in oogenesis and early development of R. prolixus, we show that lack of Rp-BicD does not significantly affect oogenesis but impairs the formation of the blastoderm. Based on our findings, we propose three times of action for maternal genes during oogenesis and embryogenesis in R. prolixus.

Identification and activity of monoamine oxidase in the orb-weaving spider Larinioides cornutus

Monoamine oxidase (MAO) is a mitochondrial membrane-bound enzyme that catalyzes the oxidative deamination of monoamines in a wide array of organisms. While the enzyme monoamine oxidase has been studied extensively in its role in moderating behavior in mammals, there is a paucity of research investigating this role in invertebrates, where the latter utilizes this enzyme in a major pathway to degrade monoamines. There is especially a dismal lack of information on how MAO influences activity in invertebrates, particularly in account of the circadian cycle. Previous studies revealed MAO degrades serotonin and norepinephrine in arachnids, but did not investigate other critically important compounds like octopamine. Larinioides cornutus is a species of orb-weaving spider that exhibits diel fluctuations in behavior, specifically levels of aggression. The monoamines octopamine and serotonin have been shown to influence aggressive behaviors in L. cornutus, thus this species was used to investigate if MAO is a potential site of regulation throughout the day. Not only did gene expression of MAO orthologs and MAO activity fluctuate at different times of day, but the enzymatic activity was substrate-specific producing a higher level of degradation of octopamine as compared to serotonin in vitro. This study further supports evidence that MAO has an active role in monoamine inactivation in invertebrates and provides a first look at how MAO ultimately may be regulating behavior in an invertebrate.

Fog signaling has diverse roles in epithelial morphogenesis in insects

The Drosophila Fog pathway represents one of the best-understood signaling cascades controlling epithelial morphogenesis. During gastrulation, Fog induces apical cell constrictions that drive the invagination of mesoderm and posterior gut primordia. The cellular mechanisms underlying primordia internalization vary greatly among insects and recent work has suggested that Fog signaling is specific to the fast mode of gastrulation found in some flies. On the contrary, here we show in the beetle Tribolium, whose development is broadly representative for insects, that Fog has multiple morphogenetic functions. It modulates mesoderm internalization and controls a massive posterior infolding involved in gut and extraembryonic development. In addition, Fog signaling affects blastoderm cellularization, primordial germ cell positioning, and cuboidal-to-squamous cell shape transitions in the extraembryonic serosa. Comparative analyses with two other distantly related insect species reveals that Fog's role during cellularization is widely conserved and therefore might represent the ancestral function of the pathway.

Shifting roles of Drosophila pair-rule gene orthologs: segmental expression and function in the milkweed bug Oncopeltus fasciatus

The discovery of pair-rule genes (PRGs) in Drosophila revealed the existence of an underlying two-segment-wide prepattern directing embryogenesis. The milkweed bug Oncopeltus fasciatus, a hemimetabolous insect, is a more representative arthropod: most of its segments form sequentially after gastrulation. Here, we report the expression and function of orthologs of the complete set of nine Drosophila PRGs in Oncopeltus Seven Of-PRG-orthologs are expressed in stripes in the primordia of every segment, rather than every other segment; Of-runt is PR-like and several orthologs are also expressed in the segment addition zone. RNAi-mediated knockdown of Of-odd-skipped, paired and sloppy-paired impacted all segments, with no indication of PR-like register. We confirm that Of-E75A is expressed in PR-like stripes, although it is not expressed in this way in Drosophila, demonstrating the existence of an underlying PR-like prepattern in Oncopeltus These findings reveal that a switch occurred in regulatory circuits, leading to segment formation: while several holometabolous insects are 'Drosophila-like', using PRG orthologs for PR patterning, most Of-PRGs are expressed segmentally in Oncopeltus, a more basally branching insect. Thus, an evolutionarily stable phenotype - segment formation - is directed by alternate regulatory pathways in diverse species.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Comparative transcriptome analysis explores maternal to zygotic transition during Eriocheir sinensis early embryogenesis

The maternal genome directs almost all aspects of early animal development. As development proceeds, the elimination of maternal gene products and zygotic genome activation (ZGA) occur during the maternal to zygotic transition (MZT). To study the molecular mechanisms regulating this developmental event in Eriocheir sinensis, RNA-Seq technology was applied to generate comprehensive information on transcriptome dynamics during early embryonic stages. In total, 32,088 annotated unigenes were obtained from the transcriptomes of fertilized eggs and embryos at the cleavage (2-4 cell) and blastula stage. A total of 566 maternal genes and 1165 zygotic genes were isolated, among which 103 and 266 genes were predicted conserved maternal transcripts (COMATs) and conserved zygotic transcripts (COZYTs), respectively. The COMATs performed housekeeping gene functions and may be essential for initiating early embryogenesis of the Bilateria. Furthermore, 87, 76 and 117 differentially expressed genes associated with the MZT, morphogenesis and immunity were identified when compared the three transcriptomic datasets. We also unmask that the MZT takes place around the cleavage stage, when the genes involved in the clearance of maternal gene products and the ZGA were significantly up-regulated. Taken together, these datasets provide a valuable resource for understanding the mechanisms of early developmental events in E. sinensis, and facilitate further studies on molecular mechanisms of asynchronous development in crabs.

Transcriptomic Analysis of Pseudoscorpion Venom Reveals a Unique Cocktail Dominated by Enzymes and Protease Inhibitors

Transcriptomic and genomic analyses have illuminated the diversity of venoms in three of the four venomous arachnid orders (scorpions, spiders, and ticks). To date, no venom gland transcriptome analysis has been available for pseudoscorpions, the fourth venomous arachnid lineage. To redress this gap, we sequenced an mRNA library generated from the venom glands of the species Synsphyronus apimelus (Garypidae). High-throughput sequencing by the Illumina protocol, followed by de novo assembly, resulted in a total of 238,331 transcripts. From those, we annotated 131 transcripts, which code for putative peptides/proteins with similar sequences to previously reported venom components available from different arachnid species in protein databases. Transcripts putatively coding for enzymes showed the richest diversity, followed by other venom components such as peptidase inhibitors, cysteine-rich peptides, and thyroglobulin 1-like peptides. Only 11 transcripts were found that code for putatively low molecular mass spider toxins. This study constitutes the first report of the diversity of components within pseudoscorpion venom.

A proteomic insight into vitellogenesis during tick ovary maturation

Ticks are arthropod ectoparasites of importance for public and veterinary health. The understanding of tick oogenesis and embryogenesis could contribute to the development of novel control methods. However, to date, studies on the temporal dynamics of proteins during ovary development were not reported. In the present study we followed protein profile during ovary maturation. Proteomic analysis of ovary extracts was performed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) using shotgun strategy, in addition to dimethyl labelling-based protein quantification. A total of 3,756 proteins were identified, which were functionally annotated into 30 categories. Circa 80% of the annotated proteins belong to categories related to basal metabolism, such as protein synthesis and modification machineries, nuclear regulation, cytoskeleton, proteasome machinery, transcriptional machinery, energetic metabolism, extracellular matrix/cell adhesion, immunity, oxidation/detoxification metabolism, signal transduction, and storage. The abundance of selected proteins involved in yolk uptake and degradation, as well as vitellin accumulation during ovary maturation, was assessed using dimethyl-labelling quantification. In conclusion, proteins identified in this study provide a framework for future studies to elucidate tick development and validate candidate targets for novel control methods.

By land, air, and sea: hemipteran diversity through the genomic lens

Thanks to a recent spate of sequencing projects, the Hemiptera are the first hemimetabolous insect order to achieve a critical mass of species with sequenced genomes, establishing the basis for comparative genomics of the bugs. However, as the most speciose hemimetabolous order, there is still a vast swathe of the hemipteran phylogeny that awaits genomic representation across subterranean, terrestrial, and aquatic habitats, and with lineage-specific and developmentally plastic cases of both wing polyphenisms and flightlessness. In this review, we highlight opportunities for taxonomic sampling beyond obvious pest species candidates, motivated by intriguing biological features of certain groups as well as the rich research tradition of ecological, physiological, developmental, and particularly cytogenetic investigation that spans the diversity of the Hemiptera.

Identification of immune-related genes in gill cells of Japanese eels (Anguilla japonica) in adaptation to water salinity changes

The changes in ambient salinity influence ion and water homeostasis, hormones secretion, and immune response in fish gills. The physiological functions of hormones and ion transporters in the regulation of gill-osmoregulation have been widely studied, however the modulation of immune response under salinity changes is not determined. Using transcriptome sequencing, we obtained a comprehensive profile of osmo-responsive genes in gill cells of Japanese eel (Anguilla japonica). Herein, we applied bioinformatics analysis to identify the immune-related genes that were significantly higher expressed in gill pavement cells (PVCs) and mitochondrial-rich cells (MRCs) in freshwater (FW) than seawater (SW) adapted fish. We validated the data using the real-time qPCR, which showed a high correlation between the RNA-seq and real-time qPCR data. In addition, the immunohistochemistry results confirmed the changes of the expression of selected immune-related genes, including C-reactive protein (CRP) in PVCs, toll-like receptor 2 (TLR2) in MRCs and interleukin-1 receptor type 2 (IL-1R2) in both PVCs and MRCs. Collectively our results demonstrated that those immune-related genes respond to salinity changes, and might trigger related special signaling pathways and network. This study provides new insights into the impacts of ambient salinity changes on adaptive immune response in fish gill cells.

Transcriptome sequencing and differential gene expression analysis of the schistosome-transmitting snail Oncomelania hupensis inhabiting hilly and marshland regions

The freshwater snail Oncomelania hupensis is the unique intermediate host of the blood fluke Schistosoma japonicum, which is the major cause of schistosomiasis. The snail inhabits two contrasting environments: the hilly and marshland regions. The hilly snails are smaller in size and have the typical smooth shell, whereas the marshland snails are larger and possess the ribbed shell. To reveal the differences in gene expression between the hilly and marshland snails, a total of six snails, three per environment, were individually examined by RNA sequencing technology. All paired-end reads were assembled into contigs from which 34,760 unigenes were predicted. Based on single nucleotide polymorphisms, principal component analysis and neighbor-joining clustering revealed two distinct clusters of hilly and marshland snails. Analysis of expression changes between environments showed that upregulated genes relating to immunity and development were enriched in hilly snails, while those associated with reproduction were over-represented in marshland snails. Eight differentially expressed genes between the two types of snails were validated by qRT-PCR. Our study identified candidate genes that could be targets for future functional studies, and provided a link between expression profiling and ecological adaptation of the snail that may have implications for schistosomiasis control.

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.

Temporal transcriptomic profiling of the ant-feeding assassin bug Acanthaspis cincticrus reveals a biased expression of genes associated with predation in nymphs

Acanthaspis cincticrus (Stål) is an assassin bug with a specialized camouflaging behavior to ambush ants in the nymphal stages. In this study, we comprehensively sequenced all the life stages of A. cincticrus, including the eggs, five nymph instars, female and male adults using Illumina HiSeq technology. We obtained 176 million clean sequence reads. The assembled 84,055 unigenes were annotated and classified functionally based on protein databases. Among the unigenes, 29.03% were annotated by one or more databases, suggesting their well-conserved functions. Comparison of the gene expression profiles in the egg, nymph and adult stages revealed certain bias. Functional enrichment analysis of significantly differentially expressed genes (SDEGs) showed positive correlation with specific physiological processes within each stage, including venom, aggression, olfactory recognition as well as growth and development. Relative expression of ten SDEGs involved in predation process was validated using quantitative real-time PCR (qRT-PCR).

Comparative transcriptomics of Entelegyne spiders (Araneae, Entelegynae), with emphasis on molecular evolution of orphan genes

Next-generation sequencing technology is rapidly transforming the landscape of evolutionary biology, and has become a cost-effective and efficient means of collecting exome information for non-model organisms. Due to their taxonomic diversity, production of interesting venom and silk proteins, and the relative scarcity of existing genomic resources, spiders in particular are excellent targets for next-generation sequencing (NGS) methods. In this study, the transcriptomes of six entelegyne spider species from three genera (Cicurina travisae, C. vibora, Habronattus signatus, H. ustulatus, Nesticus bishopi, and N. cooperi) were sequenced and de novo assembled. Each assembly was assessed for quality and completeness and functionally annotated using gene ontology information. Approximately 100 transcripts with evidence of homology to venom proteins were discovered. After identifying more than 3,000 putatively orthologous genes across all six taxa, we used comparative analyses to identify 24 instances of positively selected genes. In addition, between ~ 550 and 1,100 unique orphan genes were found in each genus. These unique, uncharacterized genes exhibited elevated rates of amino acid substitution, potentially consistent with lineage-specific adaptive evolution. The data generated for this study represent a valuable resource for future phylogenetic and molecular evolutionary research, and our results provide new insight into the forces driving genome evolution in taxa that span the root of entelegyne spider phylogeny.

Transcriptome Characterization of Dendrolimus punctatus and Expression Profiles at Different Developmental Stages

The pine moth Dendrolimus punctatus (Walker) is a common insect pest that confers serious damage to conifer forests in south of China. Extensive physiology and ecology studies on D. punctatus have been carried out, but the lack of genetic information has limited our understanding of the molecular mechanisms behind its development and resistance. Using RNA-seq approach, we characterized the transcriptome of this pine moth and investigated its developmental expression profiles during egg, larval, pupal, and adult stages. A total of 107.6 million raw reads were generated that were assembled into 70,664 unigenes. More than 30% unigenes were annotated by searching for homology in protein databases. To better understand the process of metamorphosis, we pairwise compared four developmental phases and obtained 17,624 differential expression genes. Functional enrichment analysis of differentially expressed genes showed positive correlation with specific physiological activities of each stage, and these results were confirmed by qRT-PCR experiments. This study provides a valuable genomic resource of D. punctatus covering all its developmental stages, and will promote future studies on biological processes at the molecular level.

Transcriptomic alterations in Daphnia magna embryos from mothers exposed to hypoxia

Hypoxia occurs when dissolved oxygen (DO) falls below 2.8mgL(-1) in aquatic environments. It can cause trans-generational effects not only in fish, but also in the water fleas Daphnia. In this study, transcriptome sequencing analysis was employed to identify transcriptomic alterations induced by hypoxia in embryos of Daphnia magna, with an aim to investigate the mechanism underlying the trans-generational effects caused by hypoxia in Daphnia. The embryos (F1) were collected from adults (F0) that were previously exposed to hypoxia (or normoxia) for their whole life. De novo transcriptome assembly identified 18270 transcripts that were matched to the UniProtKB/Swiss-Prot database and resulted in 7419 genes. Comparative transcriptome analysis showed 124 differentially expressed genes, including 70 up- and 54 down-regulated genes under hypoxia. Gene ontology analysis further highlighted three clusters of genes which revealed acclimatory changes of haemoglobin, suppression in vitellogenin gene family and histone modifications. Specifically, the expressions of histone H2B, H3, H4 and histone deacetylase 4 (HDAC4) were deregulated. This study suggested that trans-generational effects of hypoxia on Daphnia may be mediated through epigenetic regulations of histone modifications.

Phylogenomics and Divergence Dating of Fungus-Farming Ants (Hymenoptera: Formicidae) of the Genera Sericomyrmex and Apterostigma

Fungus-farming ("attine") ants are model systems for studies of symbiosis, coevolution, and advanced eusociality. A New World clade of nearly 300 species in 15 genera, all attine ants cultivate fungal symbionts for food. In order to better understand the evolution of ant agriculture, we sequenced, assembled, and analyzed transcriptomes of four different attine ant species in two genera: three species in the higher-attine genus Sericomyrmex and a single lower-attine ant species, Apterostigma megacephala, representing the first genomic data for either genus. These data were combined with published genomes of nine other ant species and the honey bee Apis mellifera for phylogenomic and divergence-dating analyses. The resulting phylogeny confirms relationships inferred in previous studies of fungus-farming ants. Divergence-dating analyses recovered slightly older dates than most prior analyses, estimating that attine ants originated 53.6-66.7 million of years ago, and recovered a very long branch subtending a very recent, rapid radiation of the genus Sericomyrmex. This result is further confirmed by a separate analysis of the three Sericomyrmex species, which reveals that 92.71% of orthologs have 99% - 100% pairwise-identical nucleotide sequences. We searched the transcriptomes for genes of interest, most importantly argininosuccinate synthase and argininosuccinate lyase, which are functional in other ants but which are known to have been lost in seven previously studied attine ant species. Loss of the ability to produce the amino acid arginine has been hypothesized to contribute to the obligate dependence of attine ants upon their cultivated fungi, but the point in fungus-farming ant evolution at which these losses occurred has remained unknown. We did not find these genes in any of the sequenced transcriptomes. Although expected for Sericomyrmex species, the absence of arginine anabolic genes in the lower-attine ant Apterostigma megacephala strongly suggests that the loss coincided with the origin of attine ants.

De Novo Assembly and Developmental Transcriptome Analysis of the Small White Butterfly Pieris rapae

The small white butterfly Pieris rapae is one of the most destructive pests of Brassicaceae. Yet little is understood about its genes involved in development. To facilitate research on P. rapae, we sequenced the transcriptome of P. rapae during six developmental stages, including the egg, three larval stages, the pupa, and the adult. In total, 240 million high-quality reads were obtained. De novo assembly generated 96,069 unigenes with an average length of 1353 nt. Of these, 31,629 unigenes had homologs as determined by a blastx search against the NR database with a cut-off e-value of 10⁻⁵. Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted to functionally annotate those genes. Then, 849 genes involved in seven canonical development signaling pathway were identified, including dozens of key genes such as Hippo, Notch, and JAK2. A total of 21,883 differentially expressed (cut-off of 2-fold) unigenes were detected across the developmental stages, most of which were found between the egg and first larval stages. Interestingly, only 34 differentially expressed unigenes, most of which are cuticle protein related genes, were detected with a cut-off of 2¹⁰-fold. Furthermore, we identified 32 heat shock protein (Hsp) genes that were expressed with complete open reading frames. Based on phylogenetic trees of the Hsp genes, we found that Hsp genes with close evolutionary relationships had similar expression pattern. Additionally, partial pattern recognition receptors genes were found to be developmental regulated. This study provides comprehensive sequence resources for P. rapae and numerous differential expressed genes, and these findings will lay the foundation for future functional genomics studies on this species.

De novo assembly of a tadpole shrimp (Triops newberryi) transcriptome and preliminary differential gene expression analysis

Next-generation sequencing techniques, such as RNA sequencing, have provided a wealth of genomic information for nonmodel species. Transcriptomic information can be used to quantify the patterns of gene expression, which can identify how environmental differences invoke organismal stress responses and provide a gauge in predicting species adaptability. In our study, we used RNA sequencing to characterize the first transcriptome from a naupliar tadpole shrimp (Triops newberryi) to identify the genes expressed during the early life history stages and which could be important for future genomic studies. RNA was extracted from naupliar T. newberryi that were reared in a laboratory-controlled setting and in two different water types, a native and a non-native condition. A total of six replicates, three per condition, were sequenced with the Illumina Hi-Seq 2000 achieving 365 M 50-nt reads. High-quality reads were produced and de novo assembly was used to construct a T. newberryi transcriptome that was approximately 24.8 M base pairs. More than 10 000 peptides were predicted from the assembly, and genes were sorted into gene ontology categories. The use of different water conditions allowed for a preliminary differential gene expression analysis in order to compare the changes in gene expression between conditions. There were 299 differentially expressed genes between water conditions that might serve as a focal point for future genomic studies of Triops acclimation to different environments. The Triops transcriptome could serve as vital genomic information for additional studies on Branchiopod crustaceans.

The benefits of expanding studies of trait exaggeration to hemimetabolous insects and beyond morphology

Trait exaggeration, well known to naturalists and evolutionary biologists, has recently become a prominent research subject in the modern field of Evolutionary Developmental Biology. A large number of traits that can be considered as cases of exaggeration exist in nature. Yet, the field has almost exclusively focused on the study of growth-related exaggerated traits in a selection of holometabolous insects. The absence of the hemimetabola from studies of exaggeration leaves a significant gap in our understanding of the development and evolution of such traits. Here we argue that efforts to understand the mechanisms of trait exaggeration would benefit from expanding the study subjects to include other kinds of exaggeration and other model species.

Insights into Sexual Precocity of Female Oriental River Prawn Macrobrachium nipponense through Transcriptome Analysis

Background

The oriental river prawn (Macrobrachium nipponense) is the most prevalent aquaculture species in China. The sexual precocity in this species has received considerable attention in recent years because more and more individuals matured at a small size, which devalues the commercial production. In this study, we developed deep-coverage transcriptomic sequencing data for the ovaries of sexually precocious and normal sexually mature M. nipponense using next-generation RNA sequencing technology and attempted to provide the first insight into the molecular regulatory mechanism of sexual precocity in this species.

Results

A total of 63,336 unigenes were produced from the ovarian cDNA libraries of sexually precocious and normal sexually mature M. nipponense using Illumina HiSeq 2500 platform. Through BLASTX searches against the NR, STRING, Pfam, Swissprot and KEGG databases, 15,134 unigenes were annotated, accounting for 23.89% of the total unigenes. 5,195 and 3,227 matched unigenes were categorized by GO and COG analysis respectively. 15,908 unigenes were consequently mapped into 332 KEGG pathways, and many reproduction-related pathways and genes were identified. Moreover, 26,008 SSRs were identified from 18,133 unigenes. 80,529 and 80,516 SNPs were yielded from ovarian libraries of sexually precocious and normal sexually mature prawn, respectively, and 29,851 potential SNPs between these two groups were also predicted. After comparing the ovarian libraries of sexually precocious and normal sexually mature prawn, 549 differentially expressed genes (DEGs) and 9 key DEGs that may be related to sexual precocity of M. nipponense were identified. 20 DEGs were selected for validation by quantitative real-time PCR (QPCR) and 19 DEGs show consistent expression between QPCR and RNAseq-based differential expression analysis datasets.

Conclusion

This is the first report on the large-scale RNA sequencing of ovaries of sexually precocious and normal sexually mature M. nipponense. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of M. nipponense. The large number of candidate SNPs and SSRs detected in this study could be used as genetic markers for population genetics and functional genomics in this species. More importantly, many DEGs, especially nine key DEGs between sexually precocious and normal sexually mature prawns were identified, which will dramatically improve understanding of molecular regulatory mechanism of sexual precocity of this species.

A de novo transcriptome of the Malpighian tubules in non-blood-fed and blood-fed Asian tiger mosquitoes Aedes albopictus: insights into diuresis, detoxification, and blood meal processing

Background. In adult female mosquitoes, the renal (Malpighian) tubules play an important role in the post-prandial diuresis, which removes excess ions and water from the hemolymph of mosquitoes following a blood meal. After the post-prandial diuresis, the roles that Malpighian tubules play in the processing of blood meals are not well described.

Methods. We used a combination of next-generation sequencing (paired-end RNA sequencing) and physiological/biochemical assays in adult female Asian tiger mosquitoes (Aedes albopictus) to generate molecular and functional insights into the Malpighian tubules and how they may contribute to blood meal processing (3–24 h after blood ingestion).

Results/Discussion. Using RNA sequencing, we sequenced and assembled the first de novo transcriptome of Malpighian tubules from non-blood-fed (NBF) and blood-fed (BF) mosquitoes. We identified a total of 8,232 non-redundant transcripts. The Malpighian tubules of NBF mosquitoes were characterized by the expression of transcripts associated with active transepithelial fluid secretion/diuresis (e.g., ion transporters, water channels, V-type H⁺-ATPase subunits), xenobiotic detoxification (e.g., cytochrome P450 monoxygenases, glutathione S-transferases, ATP-binding cassette transporters), and purine metabolism (e.g., xanthine dehydrogenase). We also detected the expression of transcripts encoding sodium calcium exchangers, G protein coupled-receptors, and septate junctional proteins not previously described in mosquito Malpighian tubules. Within 24 h after a blood meal, transcripts associated with active transepithelial fluid secretion/diuresis exhibited a general downregulation, whereas those associated with xenobiotic detoxification and purine catabolism exhibited a general upregulation, suggesting a reinvestment of the Malpighian tubules’ molecular resources from diuresis to detoxification. Physiological and biochemical assays were conducted in mosquitoes and isolated Malpighian tubules, respectively, to confirm that the transcriptomic changes were associated with functional consequences. In particular, in vivo diuresis assays demonstrated that adult female mosquitoes have a reduced diuretic capacity within 24 h after a blood meal. Moreover, biochemical assays in isolated Malpighian tubules showed an increase in glutathione S-transferase activity and the accumulation of uric acid (an end product of purine catabolism) within 24 h after a blood meal. Our data provide new insights into the molecular physiology of Malpighian tubules in culicine mosquitoes and reveal potentially important molecular targets for the development of chemical and/or gene-silencing insecticides that would disrupt renal function in mosquitoes.

Transcriptomic analysis reveals specific osmoregulatory adaptive responses in gill mitochondria-rich cells and pavement cells of the Japanese eel

BACKGROUND

Homeostasis of ions and water is important for the maintenance of cellular functions. The regulation of the homeostasis is particularly important in euryhaline fish that migrate between freshwater (FW) and seawater (SW) environments. The fish gill, the major tissue that forms an interface separating the extracellular fluids and external water environment, has an effective transport system to maintain and regulate a constant body osmolality. In fish gills, the two major epithelial cells, pavement cells (PVCs) and mitochondria-rich cells (MRCs), are known to play key and complementary roles in ion transport at the interface. Discovering the robust mechanisms underlying the two cell types' response to osmotic stress would benefit our understanding of the fundamental mechanism allowing PVCs and MRCs to handle osmotic stress. Owing to the limited genomic data available on estuarine species, existing knowledge in this area is slim. In this study, transcriptome analyses were conducted using PVCs and MRCs isolated from Japanese eels adapted to FW or SW environments to provide a genome-wide molecular study to unravel the fundamental processes at work.

RESULTS

The study identified more than 12,000 transcripts in the gill cells. Interestingly, remarkable differential expressed genes (DEGs) were identified in PVCs (970 transcripts) instead of MRCs (400 transcripts) in gills of fish adapted to FW or SW. Since PVCs cover more than 90 % of the gill epithelial surface, the greater change in gene expression patterns in PVCs in response to external osmolality is anticipated. In the integrity pathway analysis, 19 common biological functions were identified in PVCs and MRCs. In the enriched signaling pathways analysis, most pathways differed between PVCs and MRCs; 14 enriched pathways were identified in PVCs and 12 in MRCs. The results suggest that the osmoregulatory responses in PVCs and MRCs are cell-type specific, which supports the complementary functions of the cells in osmoregulation.

CONCLUSIONS

This is the first study to provide transcriptomic analysis of PVCs and MRCs in gills of eels adapted to FW or SW environments. It describes the cell-type specific transcriptomic network in different tonicity. The findings consolidate the known osmoregulatory pathways and provide molecular insight in osmoregulation. The presented data will be useful for researchers to select their targets for further studies.

De Novo Assembly and Transcriptome Characterization of Canine Retina Using High-Throughput Sequencing

We performed transcriptome sequencing of canine retinal tissue by 454 GS-FLX and Ion Torrent PGM platforms. RNA-Seq analysis by CLC Genomics Workbench mapped expression of 10,360 genes. Gene ontology analysis of retinal transcriptome revealed abundance of transcripts known to be involved in vision associated processes. The de novo assembly of the sequences using CAP3 generated 29,683 contigs with mean length of 560.9 and N50 of 619 bases. Further analysis of contigs predicted 3,827 full-length cDNAs and 29,481 (99%) open reading frames (ORFs). In addition, 3,782 contigs were assigned to 316 KEGG pathways which included melanogenesis, phototransduction, and retinol metabolism with 33, 15, and 11 contigs, respectively. Among the identified microsatellites, dinucleotide repeats were 68.84%, followed by trinucleotides, tetranucleotides, pentanucleotides, and hexanucleotides in proportions of 25.76, 9.40, 2.52, and 0.96%, respectively. This study will serve as a valuable resource for understanding the biology and function of canine retina.

Origin and diversification of wings: Insights from a neopteran insect

Winged insects underwent an unparalleled evolutionary radiation, but mechanisms underlying the origin and diversification of wings in basal insects are sparsely known compared with more derived holometabolous insects. In the neopteran species Oncopeltus fasciatus, we manipulated wing specification genes and used RNA-seq to obtain both functional and genomic perspectives. Combined with previous studies, our results suggest the following key steps in wing origin and diversification. First, a set of dorsally derived outgrowths evolved along a number of body segments including the first thoracic segment (T1). Homeotic genes were subsequently co-opted to suppress growth of some dorsal flaps in the thorax and abdomen. In T1 this suppression was accomplished by Sex combs reduced, that when experimentally removed, results in an ectopic T1 flap similar to prothoracic winglets present in fossil hemipteroids and other early insects. Global gene-expression differences in ectopic T1 vs. T2/T3 wings suggest that the transition from flaps to wings required ventrally originating cells, homologous with those in ancestral arthropod gill flaps/epipods, to migrate dorsally and fuse with the dorsal flap tissue thereby bringing new functional gene networks; these presumably enabled the T2/T3 wing's increased size and functionality. Third, "fused" wings became both the wing blade and surrounding regions of the dorsal thorax cuticle, providing tissue for subsequent modifications including wing folding and the fit of folded wings. Finally, Ultrabithorax was co-opted to uncouple the morphology of T2 and T3 wings and to act as a general modifier of hindwings, which in turn governed the subsequent diversification of lineage-specific wing forms.

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.

Transcriptome Changes during the Life Cycle of the Red Sponge, Mycale phyllophila (Porifera, Demospongiae, Poecilosclerida)

Sponges are an ancient metazoan group with broad ecological, evolutionary, and biotechnological importance. As in other marine invertebrates with a biphasic life cycle, the developing sponge undergoes a significant morphological, physiological, and ecological transformation during settlement and metamorphosis. In this study, we compare new transcriptome datasets for three life cycle stages of the red sponge (Mycale phyllophila) to test whether gene expression (as in the model poriferan, Amphimedon queenslandica) also varies more after settlement and metamorphosis. In contrast to A. queenslandica, we find that the transcriptome of M. phyllophila changes more during the earlier pre-competent larva/post-larva transition that spans these defining events. We also find that this transition is marked by a greater frequency of significantly up-regulated Gene Ontology terms including those for morphogenesis, differentiation, and development and that the transcriptomes of its pre-competent larvae and adult are distinct. The life cycle transcriptome variation between M. phyllophila and A. queenslandica may be due to their long separate evolutionary histories and corresponding differences in developmental rates and timing. This study now calls for new transcriptome datasets of M. phyllophila and other sponges, which will allow for tests of the generality of our life cycle expression differences and for the greater exploitation of poriferans in both basic and applied research.

Codon and Amino Acid Usage Are Shaped by Selection Across Divergent Model Organisms of the Pancrustacea

In protein-coding genes, synonymous codon usage and amino acid composition correlate to expression in some eukaryotes, and may result from translational selection. Here, we studied large-scale RNA-seq data from three divergent arthropod models, including cricket (Gryllus bimaculatus), milkweed bug (Oncopeltus fasciatus), and the amphipod crustacean Parhyale hawaiensis, and tested for optimization of codon and amino acid usage relative to expression level. We report strong signals of AT3 optimal codons (those favored in highly expressed genes) in G. bimaculatus and O. fasciatus, whereas weaker signs of GC3 optimal codons were found in P. hawaiensis, suggesting selection on codon usage in all three organisms. Further, in G. bimaculatus and O. fasciatus, high expression was associated with lowered frequency of amino acids with large size/complexity (S/C) scores in favor of those with intermediate S/C values; thus, selection may favor smaller amino acids while retaining those of moderate size for protein stability or conformation. In P. hawaiensis, highly transcribed genes had elevated frequency of amino acids with large and small S/C scores, suggesting a complex dynamic in this crustacean. In all species, the highly transcribed genes appeared to favor short proteins, high optimal codon usage, specific amino acids, and were preferentially involved in cell-cycling and protein synthesis. Together, based on examination of 1,680,067, 1,667,783, and 1,326,896 codon sites in G. bimaculatus, O. fasciatus, and P. hawaiensis, respectively, we conclude that translational selection shapes codon and amino acid usage in these three Pancrustacean arthropods.

Comparative analysis of zygotic developmental genes in Rhodnius prolixus genome shows conserved features on the tracheal developmental pathway

Most of the in-depth studies on insect developmental genetic have been carried out in the fruit fly Drosophila melanogaster, an holometabolous insect, so much more still remains to be studied in hemimetabolous insects. Having Rhodnius prolixus sequenced genome available, we search for orthologue genes of zygotic signaling pathways, segmentation, and tracheogenesis in the R. prolixus genome and in three species of Triatoma genus transcriptomes, concluding that there is a high level of gene conservation. We also study the function of two genes required for tracheal system development in D. melanogaster - R. prolixus orthologues: trachealess (Rp-trh) and empty spiracles (Rp-ems). From that we see that Rp-trh is required for early tracheal development since Rp-trh RNAi shows that the primary tracheal branches fail to form. On the other hand, Rp-ems is implied in the proper formation of the posterior tracheal branches, in a similar way to D. melanogaster. These results represent the initial characterization of the genes involved in the tracheal development of an hemimetabolous insect building a bridge between the current genomic era and V. Wigglesworth's classical studies on insects' respiratory system physiology.

Generation of a Transcriptome in a Model Lepidopteran Pest, Heliothis virescens, Using Multiple Sequencing Strategies for Profiling Midgut Gene Expression

Heliothine pests such as the tobacco budworm, Heliothis virescens (F.), pose a significant threat to production of a variety of crops and ornamental plants and are models for developmental and physiological studies. The efforts to develop new control measures for H. virescens, as well as its use as a relevant biological model, are hampered by a lack of molecular resources. The present work demonstrates the utility of next-generation sequencing technologies for rapid molecular resource generation from this species for which lacks a sequenced genome. In order to amass a de novo transcriptome for this moth, transcript sequences generated from Illumina, Roche 454, and Sanger sequencing platforms were merged into a single de novo transcriptome assembly. This pooling strategy allowed a thorough sampling of transcripts produced under diverse environmental conditions, developmental stages, tissues, and infections with entomopathogens used for biological control, to provide the most complete transcriptome to date for this species. Over 138 million reads from the three platforms were assembled into the final set of 63,648 contigs. Of these, 29,978 had significant BLAST scores indicating orthologous relationships to transcripts of other insect species, with the top-hit species being the monarch butterfly (Danaus plexippus) and silkworm (Bombyx mori). Among identified H. virescens orthologs were immune effectors, signal transduction pathways, olfactory receptors, hormone biosynthetic pathways, peptide hormones and their receptors, digestive enzymes, and insecticide resistance enzymes. As an example, we demonstrate the utility of this transcriptomic resource to study gene expression profiling of larval midguts and detect transcripts of putative Bacillus thuringiensis (Bt) Cry toxin receptors. The substantial molecular resources described in this study will facilitate development of H. virescens as a relevant biological model for functional genomics and for new biological experimentation needed to develop efficient control efforts for this and related Noctuid pest moths.

Dynamic BMP signaling polarized by Toll patterns the dorsoventral axis in a hemimetabolous insect

Toll-dependent patterning of the dorsoventral axis in Drosophila represents one of the best understood gene regulatory networks. However, its evolutionary origin has remained elusive. Outside the insects Toll is not known for a patterning function, but rather for a role in pathogen defense. Here, we show that in the milkweed bug Oncopeltus fasciatus, whose lineage split from Drosophila's more than 350 million years ago, Toll is only required to polarize a dynamic BMP signaling network. A theoretical model reveals that this network has self-regulatory properties and that shallow Toll signaling gradients are sufficient to initiate axis formation. Such gradients can account for the experimentally observed twinning of insect embryos upon egg fragmentation and might have evolved from a state of uniform Toll activity associated with protecting insect eggs against pathogens.

Transcriptomic responses of corpuscle of Stannius gland of Japanese eels (Anguilla japonica) to changes in water salinity

Physiological studies of a unique endocrine gland in fish, named corpuscles of Stannius (CS), described a Ca2(+)-regulatory function for this gland mediated by stanniocalcin-1, a hypocalcemic polypeptide hormone. However, to date, the endocrine functions of the glands have not been completely elucidated.We hypothesized that other unidentified active principles in the glands are involved in the regulation of plasma ion (Na(+), Ca2(+)) and/or blood pressure. In this study, transcriptome sequencing of CS glands was performed using Japanese eels (Anguilla japonica) adapted to freshwater (FW) or seawater (SW) to reveal the presence and differential expression of genes encoding proteins related to the ion-osmoregulatory and pressor functions. We acquired a total of 14.1 Mb and 12.1 Mb quality-trimmed reads from the CS glands collected from FW and SW adapted eels, respectively. The de novo assembly resulted in 9254 annotated genes. Among them, 475 genes were differentially expressed with 357 up- and 118 down-regulated in the SW group. Gene ontology analysis further demonstrated the presence of natriuresis and pressor related genes. In summary, ours is the first study using high-throughput sequencing to identify gene targets that could explain the physiological importance of the CS glands.

We can't all be supermodels: the value of comparative transcriptomics to the study of non-model insects

Insects are the most diverse group of organisms on the planet. Variation in gene expression lies at the heart of this biodiversity and recent advances in sequencing technology have spawned a revolution in researchers' ability to survey tissue-specific transcriptional complexity across a wide range of insect taxa. Increasingly, studies are using a comparative approach (across species, sexes and life stages) that examines the transcriptional basis of phenotypic diversity within an evolutionary context. In the present review, we summarize much of this research, focusing in particular on three critical aspects of insect biology: morphological development and plasticity; physiological response to the environment; and sexual dimorphism. A common feature that is emerging from these investigations concerns the dynamic nature of transcriptome evolution as indicated by rapid changes in the overall pattern of gene expression, the differential expression of numerous genes with unknown function, and the incorporation of novel, lineage-specific genes into the transcriptional profile.

Tissue-specific transcriptome assemblies of the marine medaka Oryzias melastigma and comparative analysis with the freshwater medaka Oryzias latipes

Background

The marine medaka Oryzias melastigma has been demonstrated as a novel model for marine ecotoxicological studies. However, the lack of genome and transcriptome reference has largely restricted the use of O. melastigma in the assessment of in vivo molecular responses to environmental stresses and the analysis of biological toxicity in the marine environment. Although O. melastigma is believed to be phylogenetically closely related to Oryzias latipes, the divergence between these two species is still largely unknown. Using Illumina high-throughput RNA sequencing followed by de novo assembly and comprehensive gene annotation, we provided transcriptomic resources for the brain, liver, ovary and testis of O. melastigma. We also investigated the possible extent of divergence between O. melastigma and O. latipes at the transcriptome level.

Results

More than 14,000 transcripts across brain, liver, ovary and testis in marine medaka were annotated, of which 5880 transcripts were orthologous between O. melastigma and O. latipes. Tissue-enriched genes were identified in O. melastigma, and Gene Ontology analysis demonstrated the functional specificity of the annotated genes in respective tissue. Lastly, the identification of marine medaka-enriched transcripts suggested the necessity of generating transcriptome dataset of O. melastigma.

Conclusions

Orthologous transcripts between O. melastigma and O. latipes, tissue-enriched genes and O. melastigma-enriched transcripts were identified. Genome-wide expression studies of marine medaka require an assembled transcriptome, and this sequencing effort has generated a valuable resource of coding DNA for a non-model species. This transcriptome resource will aid future studies assessing in vivo molecular responses to environmental stresses and those analyzing biological toxicity in the marine environment.

Electronic supplementary material

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

Sugarcane giant borer transcriptome analysis and identification of genes related to digestion

Sugarcane is a widely cultivated plant that serves primarily as a source of sugar and ethanol. Its annual yield can be significantly reduced by the action of several insect pests including the sugarcane giant borer (Telchin licus licus), a lepidopteran that presents a long life cycle and which efforts to control it using pesticides have been inefficient. Although its economical relevance, only a few DNA sequences are available for this species in the GenBank. Pyrosequencing technology was used to investigate the transcriptome of several developmental stages of the insect. To maximize transcript diversity, a pool of total RNA was extracted from whole body insects and used to construct a normalized cDNA database. Sequencing produced over 650,000 reads, which were de novo assembled to generate a reference library of 23,824 contigs. After quality score and annotation, 43% of the contigs had at least one BLAST hit against the NCBI non-redundant database, and 40% showed similarities with the lepidopteran Bombyx mori. In a further analysis, we conducted a comparison with Manduca sexta midgut sequences to identify transcripts of genes involved in digestion. Of these transcripts, many presented an expansion or depletion in gene number, compared to B. mori genome. From the sugarcane giant borer (SGB) transcriptome, a number of aminopeptidase N (APN) cDNAs were characterized based on homology to those reported as Cry toxin receptors. This is the first report that provides a large-scale EST database for the species. Transcriptome analysis will certainly be useful to identify novel developmental genes, to better understand the insect's biology and to guide the development of new strategies for insect-pest control.

De novo assembly and transcriptome analysis of the Mediterranean fruit fly Ceratitis capitata early embryos

The agricultural pest Ceratitis capitata, also known as the Mediterranean fruit fly or Medfly, belongs to the Tephritidae family, which includes a large number of other damaging pest species. The Medfly has been the first non-drosophilid fly species which has been genetically transformed paving the way for designing genetic-based pest control strategies. Furthermore, it is an experimentally tractable model, in which transient and transgene-mediated RNAi have been successfully used. We applied Illumina sequencing to total RNA preparations of 8–10 hours old embryos of C. capitata, This developmental window corresponds to the blastoderm cellularization stage. In summary, we assembled 42,614 transcripts which cluster in 26,319 unique transcripts of which 11,045 correspond to protein coding genes; we identified several hundreds of long ncRNAs; we found an enrichment of transcripts encoding RNA binding proteins among the highly expressed transcripts, such as CcTRA-2, known to be necessary to establish and, most likely, to maintain female sex of C. capitata. Our study is the first de novo assembly performed for Ceratitis capitata based on Illumina NGS technology during embryogenesis and it adds novel data to the previously published C. capitata EST databases. We expect that it will be useful for a variety of applications such as gene cloning and phylogenetic analyses, as well as to advance genetic research and biotechnological applications in the Medfly and other related Tephritidae.

T cell transcripts and T cell activities in the gills of the teleost fish sea bass (Dicentrarchus labrax)

The gills of fish are a mucosal tissue that contains T cells involved in the recognition of non-self and pathogens, and in this work we describe some features of gill-associated T cells of European sea bass, a marine model species. A whole transcriptome was obtained by deep sequencing of RNA from unstimulated gills that has been analyzed for the presence of T cell-related transcripts. Of the putative expressed sequences identified in the transcriptome, around 30 were related to main functions related to T cells including Th1/Th2/Th17/Treg cell subpopulations, thus suggesting their possible presence in the branchial epithelium. The number of T cells in the gills of sea bass, measured with the specific T cell mAb DLT15 range from 10% to 20%, and IHC analysis shows their abundance and distribution in the epithelium. Leukocytes from gills are able to proliferate in the presence of lectins ConA and PHA, as measured by flow cytometry using CFSE fluorescence incorporation, and during proliferation the number of T cells counted by immunofluorescence increased. In lectin-proliferating cells the expression of T cell-related genes TRβ, TRγ, CD4, CD8α, CD45 and IL-10 increased dramatically. Our data represent a first analysis on T cell genes and on basic T cell activities of fish gills, and suggest the presence of functionally active subpopulations of T lymphocytes in this tissue.

Comparative transcriptome analysis of two oysters, Crassostrea gigas and Crassostrea hongkongensis provides insights into adaptation to hypo-osmotic conditions

Environmental salinity creates a key barrier to limit the distribution of most aquatic organisms. Adaptation to osmotic fluctuation is believed to be a factor facilitating species diversification. Adaptive evolution often involves beneficial mutations at more than one locus. Bivalves hold great interest, with numerous species living in waters, as osmoconformers, who maintain the osmotic pressure balance mostly by free amino acids. In this study, 107,076,589 reads from two groups of Crassostrea hongkongensis were produced and the assembled into 130,629 contigs. Transcripts putatively involved in stress-response, innate immunity and cell processes were identified according to Gene ontology and KEGG pathway analyses. Comparing with the transcriptome of C. gigas to characterize the diversity of transcripts between species with osmotic divergence, we identified 182,806 high-quality single nucleotide polymorphisms (SNPs) for C. hongkongensis, and 196,779 SNPs for C. gigas. Comparison of 11,602 pairs of putative orthologs allowed for identification of 14 protein-coding genes that experienced strong positive selection (Ka/Ks>1). In addition, 45 genes that may show signs of moderate positive selection (1 ≥ Ka/Ks>0.5) were also identified. Based on Ks ratios and divergence time between the two species published previously, we estimated a neutral transcriptome-wide substitution mutation rate of 1.39 × 10(-9) per site per year. Several genes were differentially expressed across the control and treated groups of each species. This is the first time to sequence the transcriptome of C. hongkongensis and provide the most comprehensive transcriptomic resource available for it. The increasing amount of transcriptome data on Crassostrea provides an excellent resource for phylogenetic analysis. A large number of SNPs identified in this work are expected to provide valuable resources for future marker and genotyping assay development. The analysis of natural selection provides an innovative view on the adaptation within species and sets the basis for future genetic and evolutionary studies.

RNA sequencing and de novo assembly of the digestive gland transcriptome in Mytilus galloprovincialis fed with toxinogenic and non-toxic strains of Alexandrium minutum

Background

The Mediterranean mussel Mytilus galloprovincialis is marine bivalve with a relevant commercial importance as well as a key sentinel organism for the biomonitoring of environmental pollution. Here we report the RNA sequencing of the mussel digestive gland, performed with the aim: a) to produce a high quality de novo transcriptome assembly, thus improving the genetic and molecular knowledge of this organism b) to provide an initial assessment of the response to paralytic shellfish poisoning (PSP) on a molecular level, in order to identify possible molecular markers of toxin accumulation.

Results

The comprehensive de novo assembly and annotation of the transcriptome yielded a collection of 12,079 non-redundant consensus sequences with an average length of 958 bp, with a high percentage of full-length transcripts. The whole-transcriptome gene expression study indicated that the accumulation of paralytic toxins produced by the dinoflagellate Alexandrium minutum over a time span of 5 days scarcely affected gene expression, but the results need further validation with a greater number of biological samples and naturally contaminated specimens.

Conclusion

The digestive gland reference transcriptome we produced significantly improves the data collected from previous sequencing efforts and provides a basic resource for expanding functional genomics investigations in M. galloprovincialis. Although not conclusive, the results of the RNA-seq gene expression analysis support the classification of mussels as bivalves refractory to paralytic shellfish poisoning and point out that the identification molecular biomarkers of PSP in the digestive gland of this organism is problematic.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-722) contains supplementary material, which is available to authorized users.