Developmental gene discovery in a hemimetabolous insect: de novo assembly and annotation of a transcriptome for the cricket Gryllus bimaculatus

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Insights into the genomic evolution of insects from cricket genomes

Most of our knowledge of insect genomes comes from Holometabolous species, which undergo complete metamorphosis and have genomes typically under 2 Gb with little signs of DNA methylation. In contrast, Hemimetabolous insects undergo the presumed ancestral process of incomplete metamorphosis, and have larger genomes with high levels of DNA methylation. Hemimetabolous species from the Orthopteran order (grasshoppers and crickets) have some of the largest known insect genomes. What drives the evolution of these unusual insect genome sizes, remains unknown. Here we report the sequencing, assembly and annotation of the 1.66-Gb genome of the Mediterranean field cricket Gryllus bimaculatus, and the annotation of the 1.60-Gb genome of the Hawaiian cricket Laupala kohalensis. We compare these two cricket genomes with those of 14 additional insects and find evidence that hemimetabolous genomes expanded due to transposable element activity. Based on the ratio of observed to expected CpG sites, we find higher conservation and stronger purifying selection of methylated genes than non-methylated genes. Finally, our analysis suggests an expansion of the pickpocket class V gene family in crickets, which we speculate might play a role in the evolution of cricket courtship, including their characteristic chirping.

Ylla, Extavour et al. use genomic data from crickets to investigate the evolution of large genome sizes and DNA methylation events in insects. Their findings indicate that transposable element activity drove genome expansion in hemimetabolous insects, such as crickets and grasshoppers, and that DNA methylation is predominant in conserved genes.

Adaptation of codon and amino acid use for translational functions in highly expressed cricket genes

BACKGROUND

For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems.

RESULTS

Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization.

CONCLUSIONS

Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

Striking parallels between dorsoventral patterning in Drosophila and Gryllus reveal a complex evolutionary history behind a model gene regulatory network

Dorsoventral pattering relies on Toll and BMP signalling in all insects studied so far, with variations in the relative contributions of both pathways. Drosophila and the beetle Tribolium share extensive dependence on Toll, while representatives of more distantly related lineages like the wasp Nasonia and bug Oncopeltus rely more strongly on BMP signalling. Here, we show that in the cricket Gryllus bimaculatus, an evolutionarily distant outgroup, Toll has, like in Drosophila, a direct patterning role for the ventral half of the embryo. In addition, Toll polarises BMP signalling, although this does not involve the conserved BMP inhibitor Sog/Chordin. Finally, Toll activation relies on ovarian patterning mechanisms with striking similarity to Drosophila. Our data suggest two surprising hypotheses: (1) that Toll's patterning function in Gryllus and Drosophila is the result of convergent evolution or (2) a Drosophila-like system arose early in insect evolution and was extensively altered in multiple independent lineages.

Extensive Linkage and Genetic Coupling of Song and Preference Loci Underlying Rapid Speciation in Laupala Crickets

In nature, closely related species commonly display divergent mating behaviors, suggesting a central role for such traits in the origin of species. Elucidating the genetic basis of divergence in these traits is necessary to understand the evolutionary process leading to reproductive barriers and speciation. The rapidly speciating Hawaiian crickets of the genus Laupala provides an ideal system for dissecting the genetic basis of mating behavior divergence. In Laupala, closely related species differ markedly in male song pulse rate and female preference for pulse rate. These behaviors play an important role in determining mating patterns. Previous studies identified a genetic architecture consisting of numerous small to moderate effect loci causing interspecific differences in pulse rate and preference, including colocalizing pulse rate and preference QTL on linkage group one (LG1). To further interrogate these QTL, we conduct a fine mapping study using high-density SNP linkage maps. With improved statistical power and map resolution, we provide robust evidence for genetic coupling between song and preference, along with two additional pulse rate QTL on LG1, revealing a more resolved picture of the genetic architecture underlying mating behavior divergence. Our sequence-based genetic map, along with dramatically narrowed QTL confidence intervals, allowed us to annotate genes within the QTL regions and identify several exciting candidate genes underlying variation in pulse rate and preference divergence. Such knowledge suggests potential molecular mechanisms underlying the evolution of behavioral barriers.

Transcriptome analysis of life stages of the house cricket, Acheta domesticus, to improve insect crop production

To develop genetic resources for the improvement of insects as food, we sequenced transcripts from embryos, one-day hatchlings, three nymphal stages, and male and female adults of the house cricket, Acheta domesticus. A draft transcriptome was assembled from more than 138 million sequences combined from all life stages and sexes. The draft transcriptome assembly contained 45,866 contigs, and more than half were similar to sequences at NCBI (e value < e⁻³). The highest sequence identity was found in sequences from the termites Cryptotermes secundus and Zootermopsis nevadensis. Sequences with identity to Gregarina niphandrodes suggest that these crickets carry the parasite. Among all life stages, there were 5,042 genes with differential expression between life stages (significant at p < 0.05). An enrichment analysis of gene ontology terms from each life stage or sex highlighted genes that were important to biological processes in cricket development. We further characterized genes that may be important in future studies of genetically modified crickets for improved food production, including those involved in RNA interference, and those encoding prolixicin and hexamerins. The data represent an important first step in our efforts to provide genetically improved crickets for human consumption and livestock feed.

Characterization of plexinA and two distinct semaphorin1a transcripts in the developing and adult cricket Gryllus bimaculatus

Guidance cues act during development to guide growth cones to their proper targets in both the central and peripheral nervous systems. Experiments in many species indicate that guidance molecules also play important roles after development, though less is understood about their functions in the adult. The Semaphorin family of guidance cues, signaling through Plexin receptors, influences the development of both axons and dendrites in invertebrates. Semaphorin functions have been extensively explored in Drosophila melanogaster and some other Dipteran species, but little is known about their function in hemimetabolous insects. Here, we characterize sema1a and plexA in the cricket Gryllus bimaculatus. In fact, we found two distinct predicted Sema1a proteins in this species, Sema1a.1 and Sema1a.2, which shared only 48% identity at the amino acid level. We include a phylogenetic analysis that predicted that many other insect species, both holometabolous and hemimetabolous, express two Sema1a proteins as well. Finally, we used in situ hybridization to show that sema1a.1 and sema1a.2 expression patterns were spatially distinct in the embryo, and both roughly overlap with plexA. All three transcripts were also expressed in the adult brain, mainly in the mushroom bodies, though sema1a.2 was expressed most robustly. sema1a.2 was also expressed strongly in the adult thoracic ganglia while sema1a.1 was only weakly expressed and plexA was undetectable.

Hox genes limit germ cell formation in the short germ insect Gryllus bimaculatus

Hox genes are conserved transcription factor-encoding genes that specify the identity of body regions in bilaterally symmetrical animals. In the cricket Gryllus bimaculatus, a member of the hemimetabolous insect group Orthoptera, the induction of a subset of mesodermal cells to form the primordial germ cells (PGCs) is restricted to the second through the fourth abdominal segments (A2 to A4). In numerous insect species, the Hox genes Sex-combs reduced (Scr), Antennapedia (Antp), Ultrabithorax (Ubx), and abdominal-A (abd-A) jointly regulate the identities of middle and posterior body segments, suggesting that these genes may restrict PGC formation to specific abdominal segments in G. bimaculatus Here we show that reducing transcript levels of some or all of these Hox genes results in supernumerary and/or ectopic PGCs, either individually or in segment-specific combinations, suggesting that the role of these Hox genes is to limit PGC development with respect to their number, segmental location, or both. These data provide evidence of a role for this ancient group of genes in PGC development.

Genome-wide analysis of developmental stage-specific transcriptome in Bradysia odoriphaga

Bradysia odoriphaga is a serious pest of the Chinese chive; however, detailed information regarding the developmental stage-specific gene expression patterns of B. odoriphaga is not yet available. In this study, RNA sequencing (RNA-seq) was performed to determine the gene expression patterns of developmental stages including the eggs, second instar larvae, fourth instar larvae, pupae, and adults of B. odoriphaga. Analysis of 15 samples revealed an average of 89.56% of the clean reads could be mapped onto the assembled UniGene database. Cluster tree analysis showed that the expression patterns were stage-specific and that samples of the second and fourth instar larvae clustered in one group, while those of eggs, pupae, and adults clustered in another group. Differential expression unigenes (DEUs) for sequential developmental stages were between 3314 and 10,632. A total of 1910-7756 DEUs of sequential developmental stages were assigned into 45-56 gene ontology categories and 1165-3845 DEUs were mapped into Kyoto Encyclopedia of Genes and Genomes pathways. The expression of DEUs related to growth and development showed that hormone receptors highly expressed in the pupal stage, while chitinases were highly expressed in the larval stage. The results of quantitative real time polymerase chain reaction (qRT-PCR) and RNA-seq expression agreed well for 12 growth- and development-related unigenes. This study identified DEUs for sequential developmental stages of B. odoriphaga. Gene Ontology classifications and KEGG pathway identification of DEUs not only provide information useful for understanding insect growth and development but also for exploring novel approaches to control B. odoriphaga.

The Genetics of Mating Song Evolution Underlying Rapid Speciation: Linking Quantitative Variation to Candidate Genes for Behavioral Isolation

Differences in mating behaviors evolve early during speciation, eventually contributing to reproductive barriers between species. Knowledge of the genetic and genomic basis of these behaviors is therefore integral to a causal understanding of speciation. Acoustic behaviors are often part of the mating ritual in animal species. The temporal rhythms of mating songs are notably species-specific in many vertebrates and arthropods and often underlie assortative mating. Despite discoveries of mutations that disrupt the temporal rhythm of these songs, we know surprisingly little about genes affecting naturally occurring variation in the temporal pattern of singing behavior. In the rapidly speciating Hawaiian cricket genus Laupala, the striking species variation in song rhythms constitutes a behavioral barrier to reproduction between species. Here, we mapped the largest-effect locus underlying interspecific variation in song rhythm between two Laupala species to a narrow genomic region, wherein we find no known candidate genes affecting song temporal rhythm in Drosophila Whole-genome sequencing, gene prediction, and functional annotation of this region reveal an exciting and promising candidate gene, the putative cyclic nucleotide-gated ion channel-like gene, for natural variation in mating behavior. Identification and molecular characterization of the candidate gene reveals a nonsynonymous mutation in a conserved binding domain, suggesting that ion channels are important targets of selection on rhythmic signaling during establishment of behavioral isolation and rapid speciation.

The Cricket Gryllus bimaculatus: Techniques for Quantitative and Functional Genetic Analyses of Cricket Biology

All extant species are an outcome of nature's "experiments" during evolution, and hence multiple species need to be studied and compared to gain a thorough understanding of evolutionary processes. The field of evolutionary developmental biology (evo-devo) aspires to expand the number of species studied, because most functional genetic studies in animals have been limited to a small number of "traditional" model organisms, many of which belong to the same phylum (Chordata). The phylum Arthropoda, and particularly its component class Insecta, possesses many important characteristics that are considered favorable and attractive for evo-devo research, including an astonishing diversity of extant species and a wide disparity in body plans. The development of the most thoroughly investigated insect genetic model system to date, the fruit fly Drosophila melanogaster (a holometabolous insect), appears highly derived with respect to other insects and indeed with respect to most arthropods. In comparison, crickets (a basally branching hemimetabolous insect lineage compared to the Holometabola) are thought to embody many developmental features that make them more representative of insects. Here we focus on crickets as emerging models to study problems in a wide range of biological areas and summarize the currently available molecular, genomic, forward and reverse genetic, imaging and computational tool kit that has been established or adapted for cricket research. With an emphasis on the cricket species Gryllus bimaculatus, we highlight recent efforts made by the scientific community in establishing this species as a laboratory model for cellular biology and developmental genetics. This broad toolkit has the potential to accelerate many traditional areas of cricket research, including studies of adaptation, evolution, neuroethology, physiology, endocrinology, regeneration, and reproductive behavior. It may also help to establish newer areas, for example, the use of crickets as animal infection model systems and human food sources.

Satellite DNAs are conserved and differentially transcribed among Gryllus cricket species

Satellite DNA (satDNA) is an abundant class of non-coding repetitive DNA that is preferentially found as tandemly repeated arrays in gene-poor heterochromatin but is also present in gene-rich euchromatin. Here, we used DNA- and RNA-seq from Gryllus assimilis to address the content and transcriptional patterns of satDNAs. We also mapped RNA-seq libraries for other Gryllus species against the satDNAs found in G. assimilis and G. bimaculatus genomes to investigate their evolutionary conservation and transcriptional profiles in Gryllus. Through DNA-seq read clustering analysis using RepeatExplorer, dotplots analysis and fluorescence in situ hybridization mapping, we found that ∼4% of the G. assimilis genome is represented by 11 well-defined A + T-rich satDNA families. These are mainly located in heterochromatic areas, with some repeats able to form high-order repeat structures. By in silico transcriptional analysis we identified satDNAs that are conserved in Gryllus but differentially transcribed. The data regarding satDNA presence in G. assimilis genome were discussed in an evolutionary context, with transcriptional data enabling comparisons between sexes and across tissues when possible. We discuss hypotheses for the conservation and transcription of satDNAs in Gryllus, which might result from their role in sexual differentiation at the chromatin level, heterochromatin formation and centromeric function.

De novo assembly of a transcriptome for the cricket Gryllus bimaculatus prothoracic ganglion: An invertebrate model for investigating adult central nervous system compensatory plasticity

The auditory system of the cricket, Gryllus bimaculatus, demonstrates an unusual amount of anatomical plasticity in response to injury, even in adults. Unilateral removal of the ear causes deafferented auditory neurons in the prothoracic ganglion to sprout dendrites across the midline, a boundary they typically respect, and become synaptically connected to the auditory afferents of the contralateral ear. The molecular basis of this sprouting and novel synaptogenesis in the adult is not understood. We hypothesize that well-conserved developmental guidance cues may recapitulate their guidance functions in the adult in order to facilitate this compensatory growth. As a first step in testing this hypothesis, we have generated a de novo assembly of a prothoracic ganglion transcriptome derived from control and deafferented adult individuals. We have mined this transcriptome for orthologues of guidance molecules from four well-conserved signaling families: Slit, Netrin, Ephrin, and Semaphorin. Here we report that transcripts encoding putative orthologues of most of the candidate developmental ligands and receptors from these signaling families were present in the assembly, indicating expression in the adult G. bimaculatus prothoracic ganglion.

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.

Evolutionary coincidence of adaptive changes in exuperantia and the emergence of bicoid in Cyclorrhapha (Diptera)

The great radiation in the infraorder Cyclorrhapha involved several morphological and molecular changes, including important changes in anterior egg development. During Drosophila oogenesis, exuperantia (exu) is critical for localizing bicoid (bcd) messenger RNA (mRNA) to the anterior region of the oocyte. Because it is phylogenetically older than bcd, which is exclusive to Cyclorrhapha, we hypothesize that exu has undergone adaptive changes to enable this new function. Although exu has been well studied in Drosophila, there is no functional or transcriptional information about it in any other Diptera. Here, we investigate exu in the South American fruit fly Anastrepha fraterculus, a Cyclorrhapha of great agricultural importance that have lost bcd, aiming to understand the evolution of exu in this infraorder. We assessed its pattern of gene expression in A. fraterculus by analyzing transcriptomes from cephalic and reproductive tissues. A combination of next-generation data with classical sequencing procedures enabled identification of the structure of exu and its alternative transcripts in this species. In addition to the sex-specific isoforms described for Drosophila, we found that not only exu is expressed in heads, but this is mediated by two transcripts with a specific 5′UTR exon—likely a result from usage of a third promoter. Furthermore, we tested the hypothesis that exu is evolving under positive selection in Cyclorrhapha after divergence from lower Diptera. We found evidence of positive selection at two important exu domains, EXO-like and SAM-like, both involved with mRNA binding during bcd mRNA localization in Drosophila, which could reflect its cooptation for the new function of bcd mRNA localization in Cyclorrhapha.

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 profiling of ontogeny in the acridid grasshopper Chorthippus biguttulus

Acridid grasshoppers (Orthoptera:Acrididae) are widely used model organisms for developmental, evolutionary, and neurobiological research. Although there has been recent influx of orthopteran transcriptomic resources, many use pooled ontogenetic stages obscuring information about changes in gene expression during development. Here we developed a de novo transcriptome spanning 7 stages in the life cycle of the acridid grasshopper Chorthippus biguttulus. Samples from different stages encompassing embryonic development through adults were used for transcriptomic profiling, revealing patterns of differential gene expression that highlight processes in the different life stages. These patterns were validated with semi-quantitative RT-PCR. Embryonic development showed a strongly differentiated expression pattern compared to all of the other stages and genes upregulated in this stage were involved in signaling, cellular differentiation, and organ development. Our study is one of the first to examine gene expression during post-embryonic development in a hemimetabolous insect and we found that only the fourth and fifth instars had clusters of genes upregulated during these stages. These genes are involved in various processes ranging from synthesis of biogenic amines to chitin binding. These observations indicate that post-embryonic ontogeny is not a continuous process and that some instars are differentiated. Finally, genes upregulated in the imago were generally involved in aging and immunity. Our study highlights the importance of looking at ontogeny as a whole and indicates promising directions for future research in orthopteran development.

Characterization and analysis of a de novo transcriptome from the pygmy grasshopper Tetrix japonica

The pygmy grasshopper Tetrix japonica is a common insect distributed throughout the world, and it has the potential for use in studies of body colour polymorphism, genomics and the biology of Tetrigoidea (Insecta: Orthoptera). However, limited biological information is available for this insect. Here, we conducted a de novo transcriptome study of adult and larval T. japonica to provide a better understanding of its gene expression and develop genomic resources for future work. We sequenced and explored the characteristics of the de novo transcriptome of T. japonica using Illumina HiSeq 2000 platform. A total of 107 608 206 paired-end clean reads were assembled into 61 141 unigenes using the trinity software; the mean unigene size was 771 bp, and the N50 length was 1238 bp. A total of 29 225 unigenes were functionally annotated to the NCBI nonredundant protein sequences (Nr), NCBI nonredundant nucleotide sequences (Nt), a manually annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of putative genes that are potentially involved in pigment pathways, juvenile hormone (JH) metabolism and signalling pathways were identified in the T. japonica transcriptome. Additionally, 165 769 and 156 796 putative single nucleotide polymorphisms occurred in the adult and larvae transcriptomes, respectively, and a total of 3162 simple sequence repeats were detected in this assembly. This comprehensive transcriptomic data for T. japonica will provide a usable resource for gene predictions, signalling pathway investigations and molecular marker development for this species and other pygmy grasshoppers.

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.

Deep Sequencing-Based Transcriptome Analysis Reveals the Regulatory Mechanism of Bemisia tabaci (Hemiptera: Aleyrodidae) Nymph Parasitized by Encarsia sophia (Hymenoptera: Aphelinidae)

The whitefly Bemisia tabaci is a genetically diverse complex with multiple cryptic species, and some are the most destructive invasive pests of many ornamentals and crops worldwide. Encarsia sophia is an autoparasitoid wasp that demonstrated high efficiency as bio-control agent of whiteflies. However, the immune mechanism of B. tabaci parasitization by E. sophia is unknown. In order to investigate immune response of B. tabaci to E. Sophia parasitization, the transcriptome of E. sophia parasitized B. tabaci nymph was sequenced by Illumina sequencing. De novo assembly generated 393,063 unigenes with average length of 616 bp, in which 46,406 unigenes (15.8% of all unigenes) were successfully mapped. Parasitization by E. sophia had significant effects on the transcriptome profile of B. tabaci nymph. A total of 1482 genes were significantly differentially expressed, of which 852 genes were up-regulated and 630 genes were down-regulated. These genes were mainly involved in immune response, development, metabolism and host signaling pathways. At least 52 genes were found to be involved in the host immune response, 33 genes were involved in the development process, and 29 genes were involved in host metabolism. Taken together, the assembled and annotated transcriptome sequences provided a valuable genomic resource for further understanding the molecular mechanism of immune response of B. tabaci parasitization by E. sophia.

A genes eye view of ontogeny: de novo assembly and profiling of the Gryllus rubens transcriptome

Crickets (Orthoptera:Gryllidae) are widely used model organisms for developmental, evolutionary, neurobiological and behavioural research. Here, we developed a de novo transcriptome from pooled RNA-seq Illumina data spanning seven stages in the life cycle of Gryllus rubens. Approximately 705 Mbp of data was assembled and filtered to form 27 312 transcripts. We were able to annotate 52% of our transcripts using BLAST and assign at least one gene ontology term to 41%. Pooled samples from three different ontogenetic stages were used for transcriptomic profiling revealing patterns of differential gene expression that highlight processes in the different life stages. Embryonic and early instar development was enriched for ecdysteroid metabolism, cytochrome P450s and glutathione production. Late instar development was enriched for regulation of gene expression and many of the genes highly expressed during this stage were involved in conserved developmental signalling pathways suggesting that these developmental pathways are active beyond embryonic development. Adults were enriched for fat transport (mostly relating to egg production) and production of octopamine, an important neurohormone. We also identified genes involved in conserved developmental pathways (Hedgehog, Hippo, Wnt, JAK/STAT, TGF-beta, Notch, and MEK/ERK). This is the first transcriptome spanning ontogeny in Gryllus rubens and a valuable resource for future work on development and evolution in Orthoptera.

The insect central complex as model for heterochronic brain development-background, concepts, and tools

The adult insect brain is composed of neuropils present in most taxa. However, the relative size, shape, and developmental timing differ between species. This diversity of adult insect brain morphology has been extensively described while the genetic mechanisms of brain development are studied predominantly in Drosophila melanogaster. However, it has remained enigmatic what cellular and genetic mechanisms underlie the evolution of neuropil diversity or heterochronic development. In this perspective paper, we propose a novel approach to study these questions. We suggest using genome editing to mark homologous neural cells in the fly D. melanogaster, the beetle Tribolium castaneum, and the Mediterranean field cricket Gryllus bimaculatus to investigate developmental differences leading to brain diversification. One interesting aspect is the heterochrony observed in central complex development. Ancestrally, the central complex is formed during embryogenesis (as in Gryllus) but in Drosophila, it arises during late larval and metamorphic stages. In Tribolium, it forms partially during embryogenesis. Finally, we present tools for brain research in Tribolium including 3D reconstruction and immunohistochemistry data of first instar brains and the generation of transgenic brain imaging lines. Further, we characterize reporter lines labeling the mushroom bodies and reflecting the expression of the neuroblast marker gene Tc-asense, respectively.

The transcriptional repressor Blimp-1 acts downstream of BMP signaling to generate primordial germ cells in the cricket Gryllus bimaculatus

Segregation of the germ line from the soma is an essential event for transmission of genetic information across generations in all sexually reproducing animals. Although some well-studied systems such as Drosophila and Xenopus use maternally inherited germ determinants to specify germ cells, most animals, including mice, appear to utilize zygotic inductive cell signals to specify germ cells during later embryogenesis. Such inductive germ cell specification is thought to be an ancestral trait of Bilateria, but major questions remain as to the nature of an ancestral mechanism to induce germ cells, and how that mechanism evolved. We previously reported that BMP signaling-based germ cell induction is conserved in both the mouse Mus musculus and the cricket Gryllus bimaculatus, which is an emerging model organism for functional studies of induction-based germ cell formation. In order to gain further insight into the functional evolution of germ cell specification, here we examined the Gryllus ortholog of the transcription factor Blimp-1 (also known as Prdm1), which is a widely conserved bilaterian gene known to play a crucial role in the specification of germ cells in mice. Our functional analyses of the Gryllus Blimp-1 ortholog revealed that it is essential for Gryllus primordial germ cell development, and is regulated by upstream input from the BMP signaling pathway. This functional conservation of the epistatic relationship between BMP signaling and Blimp-1 in inductive germ cell specification between mouse and cricket supports the hypothesis that this molecular mechanism regulated primordial germ cell specification in a last common bilaterian ancestor.

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.

Short germ insects utilize both the ancestral and derived mode of Polycomb group-mediated epigenetic silencing of Hox genes

In insect species that undergo long germ segmentation, such as Drosophila, all segments are specified simultaneously at the early blastoderm stage. As embryogenesis progresses, the expression boundaries of Hox genes are established by repression of gap genes, which is subsequently replaced by Polycomb group (PcG) silencing. At present, however, it is not known whether patterning occurs this way in a more ancestral (short germ) mode of embryogenesis, where segments are added gradually during posterior elongation. In this study, two members of the PcG family, Enhancer of zeste (E(z)) and Suppressor of zeste 12 (Su(z)12), were analyzed in the short germ cricket, Gryllus bimaculatus. Results suggest that although stepwise negative regulation by gap and PcG genes is present in anterior members of the Hox cluster, it does not account for regulation of two posterior Hox genes, abdominal-A (abd-A) and Abdominal-B (Abd-B). Instead, abd-A and Abd-B are predominantly regulated by PcG genes, which is the mode present in vertebrates. These findings suggest that an intriguing transition of the PcG-mediated silencing of Hox genes may have occurred during animal evolution. The ancestral bilaterian state may have resembled the current vertebrate mode of regulation, where PcG-mediated silencing of Hox genes occurs before their expression is initiated and is responsible for the establishment of individual expression domains. Then, during insect evolution, the repression by transcription factors may have been acquired in anterior Hox genes of short germ insects, while PcG silencing was maintained in posterior Hox genes.

Embryonic development of the cricket Gryllus bimaculatus

Extensive research into Drosophila melanogaster embryogenesis has improved our understanding of insect developmental mechanisms. However, Drosophila development is thought to be highly divergent from that of the ancestral insect and arthropod in many respects. We therefore need alternative models for arthopod development that are likely to be more representative of basally-branching clades. The cricket Gryllus bimaculatus is such a model, and currently has the most sophisticated functional genetic toolkit of any hemimetabolous insect. The existing cricket embryonic staging system is fragmentary, and it is based on morphological landmarks that are not easily visible on a live, undissected egg. To address this problem, here we present a complementary set of "egg stages" that serve as a guide for identifying the developmental progress of a cricket embryo from fertilization to hatching, based solely on the external appearance of the egg. These stages were characterized using a combination of brightfield timelapse microscopy, timed brightfield micrographs, confocal microscopy, and measurements of egg dimensions. These egg stages are particularly useful in experiments that involve egg injection (including RNA interference, targeted genome modification, and transgenesis), as injection can alter the speed of development, even in control treatments. We also use 3D reconstructions of fixed embryo preparations to provide a comprehensive description of the morphogenesis and anatomy of the cricket embryo during embryonic rudiment assembly, germ band formation, elongation, segmentation, and appendage formation. Finally, we aggregate and schematize a variety of published developmental gene expression patterns. This work will facilitate further studies on G. bimaculatus development, and serve as a useful point of reference for other studies of wild type and experimentally manipulated insect development in fields from evo-devo to disease vector and pest management.

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.

A comprehensive transcriptome provides candidate genes for sex determination/differentiation and SSR/SNP markers in yellow catfish

Sex dimorphic growth pattern has significant theory and application implications in fish. Recently, a Y- and X-specific allele marker-assisted sex control technique has been developed for mass production of all-male population in yellow catfish (Pelteobagrus fulvidraco), but the genetic information for sex determination and sex control breeding has remained unclear. Here, we attempted to provide the first insight into a comprehensive transcriptome covering multiple tissues from XX females, XY males, and YY super-males of yellow catfish by using 454 GS-FLX platform, for a better assembly and gene coverage. A total of 1,202,933 high quality reads (about 540 Mbp) were obtained and assembled into 28,297 contigs and 141,951 singletons. BLASTX searches against the NCBI non-redundant protein database (nr) led a total of 52,564 unique sequences including 18,748 contigs and 33,816 singletons to match 25,669 known or predicted unique proteins. All of them with annotated function were categorized by gene ontology (GO) analysis, and 712 were assigned to reproduction and reproductive process. Some potential genes relevant to reproductive system including steroid hormone biosynthesis and GnRH (gonadotropin-releasing hormone) signaling pathway were further identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis; and at least 21 sex determination and differentiation-related genes, such as Dmrt1, Sox9a/b, Cyp19b, WT1, and AMH were identified and characterized. Additionally, a total of 82,794 simple sequence repeats (SSRs), 26,450 single nucleotide polymorphisms (SNPs), and 4,145 insertions and deletions (INDELs) were revealed from the transcriptome data. Therefore, the current transcriptome resources highlight further studies on sex-control breeding in yellow catfish and will benefit future studies on reproduction and sex determination in teleost fish.

De novo sequencing, assembly and characterization of antennal transcriptome of Anomala corpulenta Motschulsky (Coleoptera: Rutelidae)

Background

Anomala corpulenta is an important insect pest and can cause enormous economic losses in agriculture, horticulture and forestry. It is widely distributed in China, and both larvae and adults can cause serious damage. It is difficult to control this pest because the larvae live underground. Any new control strategy should exploit alternatives to heavily and frequently used chemical insecticides. However, little genetic research has been carried out on A. corpulenta due to the lack of genomic resources. Genomic resources could be produced by next generation sequencing technologies with low cost and in a short time. In this study, we performed de novo sequencing, assembly and characterization of the antennal transcriptome of A. corpulenta.

Results

Illumina sequencing technology was used to sequence the antennal transcriptome of A. corpulenta. Approximately 76.7 million total raw reads and about 68.9 million total clean reads were obtained, and then 35,656 unigenes were assembled. Of these unigenes, 21,463 of them could be annotated in the NCBI nr database, and, among the annotated unigenes, 11,154 and 6,625 unigenes could be assigned to GO and COG, respectively. Additionally, 16,350 unigenes could be annotated in the Swiss-Prot database, and 14,499 unigenes could map onto 258 pathways in the KEGG Pathway database. We also found 24 unigenes related to OBPs, 6 to CSPs, and in total 167 unigenes related to chemodetection. We analyzed 4 OBPs and 3CSPs sequences and their RT-qPCR results agreed well with their FPKM values.

Conclusion

We produced the first large-scale antennal transcriptome of A. corpulenta, which is a species that has little genomic information in public databases. The identified chemodetection unigenes can promote the molecular mechanistic study of behavior in A. corpulenta. These findings provide a general sequence resource for molecular genetics research on A. corpulenta.

A comprehensive reference transcriptome resource for the common house spider Parasteatoda tepidariorum

Parasteatoda tepidariorum is an increasingly popular model for the study of spider development and the evolution of development more broadly. However, fully understanding the regulation and evolution of P. tepidariorum development in comparison to other animals requires a genomic perspective. Although research on P. tepidariorum has provided major new insights, gene analysis to date has been limited to candidate gene approaches. Furthermore, the few available EST collections are based on embryonic transcripts, which have not been systematically annotated and are unlikely to contain transcripts specific to post-embryonic stages of development. We therefore generated cDNA from pooled embryos representing all described embryonic stages, as well as post-embryonic stages including nymphs, larvae and adults, and using Illumina HiSeq technology obtained a total of 625,076,514 100-bp paired end reads. We combined these data with 24,360 ESTs available in GenBank, and 1,040,006 reads newly generated from 454 pyrosequencing of a mixed-stage embryo cDNA library. The combined sequence data were assembled using a custom de novo assembly strategy designed to optimize assembly product length, number of predicted transcripts, and proportion of raw reads incorporated into the assembly. The de novo assembly generated 446,427 contigs with an N50 of 1,875 bp. These sequences obtained 62,799 unique BLAST hits against the NCBI non-redundant protein data base, including putative orthologs to 8,917 Drosophila melanogaster genes based on best reciprocal BLAST hit identity compared with the D. melanogaster proteome. Finally, we explored the utility of the transcriptome for RNA-Seq studies, and showed that this resource can be used as a mapping scaffold to detect differential gene expression in different cDNA libraries. This resource will therefore provide a platform for future genomic, gene expression and functional approaches using P. tepidariorum.

BMP signaling is required for the generation of primordial germ cells in an insect

Two modes of germ cell formation are known in animals. Specification through maternally inherited germ plasm occurs in many well-characterized model organisms, but most animals lack germ plasm by morphological and functional criteria. The only known alternative mechanism is induction, experimentally described only in mice, which specify germ cells through bone morphogenetic protein (BMP) signal-mediated induction of a subpopulation of mesodermal cells. Until this report, no experimental evidence of an inductive germ cell signal for specification has been available outside of vertebrates. Here we provide functional genetic experimental evidence consistent with a role for BMP signaling in germ cell formation in a basally branching insect. We show that primordial germ cells of the cricket Gryllus bimaculatus transduce BMP signals and require BMP pathway activity for their formation. Moreover, increased BMP activity leads to ectopic and supernumerary germ cells. Given the commonality of BMP signaling in mouse and cricket germ cell induction, we suggest that BMP-based germ cell formation may be a shared ancestral mechanism in animals.