De novo transcriptome analysis of wing development-related signaling pathways in Locusta migratoria manilensis and Ostrinia furnacalis (Guenée)

Genetics of flight in spongy moths (Lymantria dispar ssp.): functionally integrated profiling of a complex invasive trait

BACKGROUND

Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi-species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability-female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight-relevant genetic divergence between ESM and ASM: reduced representation genome-wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight-associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution.

RESULTS

Candidate gene sets were typically non-overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight-relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits.

CONCLUSIONS

Our study identified a suite of flight-associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

A Combinatorial Single-Molecule Real-Time and Illumina Sequencing Analysis of Postembryonic Gene Expression in the Asian Citrus Psyllid Diaphorina citri

Huanglongbing (HLB) is a systemic plant disease caused by 'Candidatus Liberibacter asiaticus (CLas)' and transmitted by Diaphorina citri. D. citri acquires the CLas bacteria in the nymph stage and transmits it in the adult stage, indicating that molting from the nymph to adult stages is crucial for HLB transmission. However, the available D. citri reference genomes are incomplete, and gene function studies have been limited to date. In the current research, PacBio single-molecule real-time (SMRT) and Illumina sequencing were performed to investigate the transcriptome of D. citri nymphs and adults. In total, 10,641 full-length, non-redundant transcripts (FLNRTs), 594 alternative splicing (AS) events, 4522 simple sequence repeats (SSRs), 1086 long-coding RNAs (lncRNAs), 281 transcription factors (TFs), and 4459 APA sites were identified. Furthermore, 3746 differentially expressed genes (DEGs) between nymphs and adults were identified, among which 30 DEGs involved in the Hippo signaling pathway were found. Reverse transcription-quantitative PCR (RT-qPCR) further validated the expression levels of 12 DEGs and showed a positive correlation with transcriptome data. Finally, the spatiotemporal expression pattern of genes involved in the Hippo signaling pathway exhibited high expression in the D. citri testis, ovary, and egg. Silencing of the D. citri transcriptional co-activator (DcYki) gene significantly increased D. citri mortality and decreased the cumulative molting. Our results provide useful information and a reliable data resource for gene function research of D. citri.

Notch Signaling in Insect Development: A Simple Pathway with Diverse Functions

Notch signaling is an evolutionarily conserved pathway which functions between adjacent cells to establish their distinct identities. Despite operating in a simple mechanism, Notch signaling plays remarkably diverse roles in development to regulate cell fate determination, organ growth and tissue patterning. While initially discovered and characterized in the model insect Drosophila melanogaster, recent studies across various insect species have revealed the broad involvement of Notch signaling in shaping insect tissues. This review focuses on providing a comprehensive picture regarding the roles of the Notch pathway in insect development. The roles of Notch in the formation and patterning of the insect embryo, wing, leg, ovary and several specific structures, as well as in physiological responses, are summarized. These results are discussed within the developmental context, aiming to deepen our understanding of the diversified functions of the Notch signaling pathway in different insect species.

Cuticle Protein LmACP19 Is Required for the Stability of Epidermal Cells in Wing Development and Morphogenesis of Locusta migratoria

Insect wing consists of a double layer of epidermal cells that produce and secrete the dorsal and ventral cuticular components. It is important for the stability of epidermal cells during wing development and morphogenesis, but its specific gene expression and physiological function during this process remain unclear. In our previous work, a wing cuticle protein gene LmACP19 was identified in Locusta migratoria based on transcriptomic data. Here, we report on its roles in wing development and morphogenesis. LmACP19 encodes a chitin-binding protein belonging to RR-2 subfamily of CPR family, which is highly homologous to CP19-like proteins in other insect species. RT-qPCR analysis revealed that LmACP19 is highly expressed in wing pads of fifth-instar nymphs, and its encoded protein is located in two layers of epidermal cells but not in the cuticle. Suppression of LmACP19 by RNA interference led to abnormal wing pad and wing morphogenesis with curved, unclosed, and wrinkled phenotypes during nymph-to-nymph and nymph-to-adult transition, respectively. Furthermore, deficiency of LmACP19 affected arrangement of epidermal cells, resulting in apoptosis. Our results indicate that LmACP19 is indispensable for wing development and normal morphological structure by maintaining the stability of epidermal cells during L. migratoria molting.

The transcription factor of the Hippo signaling pathway, LmSd, regulates wing development in Locusta migratoria

Scalloped (Sd) is transcription factor that regulates cell proliferation and organ growth in the Hippo pathway. In the present research, LmSd was identified and characterized, and found to encode an N-terminal TEA domain and a C-terminal YBD domain. qRT-PCR showed that the LmSd transcription level was highest in the fifth-instar nymphs and very little was expressed in embryos. Tissue-specific analyses showed that LmSd was highly expressed in the wing. Immunohistochemistry indicated that LmSd was highly abundant in the head, prothorax, and legs during embryonic development. LmSd dsRNA injection resulted in significantly down-regulated transcription and protein expression levels compared with dsGFP injection. Gene silencing of LmSd resulted in deformed wings that were curved, wrinkled, and failed to fully expand. Approximately 40% of the nymphs had wing pads that were not able to close normally during molting from fifth-instar nymphs into adults. After silencing of LmSd, the transcription levels of cell division genes were suppressed and the expression levels of apoptosis genes were significantly up-regulated. Our results reveal that LmSd plays an important role in wing formation and development by controlling cell proliferation and inhibiting apoptosis.

Panaxadiol saponins treatment caused the subtle variations in the global transcriptional state of Asiatic corn borer, Ostrinia furnacalis

Background

Methods

Results

Conclusion

Comparative genomics sheds light on the predatory lifestyle of accipitrids and owls

Raptors are carnivorous birds including accipitrids (Accipitridae, Accipitriformes) and owls (Strigiformes), which are diurnal and nocturnal, respectively. To examine the evolutionary basis of adaptations to different light cycles and hunting behavior between accipitrids and owls, we de novo assembled besra (Accipiter virgatus, Accipitridae, Accipitriformes) and oriental scops owl (Otus sunia, Strigidae, Strigiformes) draft genomes. Comparative genomics demonstrated four PSGs (positively selected genes) (XRCC5, PRIMPOL, MDM2, and SIRT1) related to the response to ultraviolet (UV) radiation in accipitrids, and one PSG (ALCAM) associated with retina development in owls, which was consistent with their respective diurnal/nocturnal predatory lifestyles. We identified five accipitrid-specific and two owl-specific missense mutations and most of which were predicted to affect the protein function by PolyPhen-2. Genome comparison showed the diversification of raptor olfactory receptor repertoires, which may reflect an important role of olfaction in their predatory lifestyle. Comparison of TAS2R gene (i.e. linked to tasting bitterness) number in birds with different dietary lifestyles suggested that dietary toxins were a major selective force shaping the diversity of TAS2R repertoires. Fewer TAS2R genes in raptors reflected their carnivorous diet, since animal tissues are less likely to contain toxins than plant material. Our data and findings provide valuable genomic resources for studying the genetic mechanisms of raptors' environmental adaptation, particularly olfaction, nocturnality and response to UV radiation.

A head transcriptome provides insights into odorant binding proteins of the bamboo grasshopper

The bamboo grasshopper Ceracris kiangsu is a famous bamboo pest in China. The identification of genes involved in olfactory behavior of C. kiangsu is necessary for better understanding the molecular basis and expression profiles of behavior ecology. However, necessary genomic and transcriptomic data are lacking in the species, limiting control efficiency. The primary objective of this study was to find and describe odorant binding proteins in the head of the bamboo grasshopper. We performed the paired-end sequencing on an Illumina Hiseq2000 following the vendor's recommended protocol. Functional annotation was performed by comparison with public databases. OBP genes were first identified using BLASTN and BLASTX results from our C. kiangsu datebase, which was established from the date of transcriptome sequencing. The gene-specific primers were used to conduct RT-PCR to detect the tissue distribution of OBPs using a SYBR Premix ExTaq kit following the manufacturer's instructions with a real-time thermal cycler. We obtained more than 133 million clean reads derived from the C. Kiangsu heads using the next-generation sequencing, which were assembled into 260,822 unique sequences (average 814 bp). We have detected eight putative odorant binding protein genes (OBPs) of C. kiangsu for the first time, and analyzed the expression profiles of the OBPs in different tissues (head, antenna, mouthpart, body and leg). Our results reveal that the eight OBPs display a clear divergence, strongly indicating that they possessed diverse functions, and thus provides comprehensive sequence analysis for elucidating the molecular basis of OBPs in C. kiangsu. In addition, we find that the relative expression levels of OBP1, OBP2 and OBP8 are significantly higher in the antennae as compared to the other OBP genes, suggesting that these three OBP genes play crucial roles in the locust's odorant discrimination. In general, this is the first study to characterize the complete head transcriptome of C. kiangsu using high-throughput sequencing. The study opens a window for functional characterization of the OBPs of C. kiangsu, with potential for new or refined applications of semiochemicals for control of this notorious pest.

De Novo Assembly and Characterization of the Xenocatantops brachycerus Transcriptome

Grasshoppers are common pests but also have high nutritional and commercial potential. Xenocatantops brachycerus Willemse (Orthoptera: Acrididae) is an economically important grasshopper species that is reared in China. Using the IlluminaHiSeqTM 4000 platform, three transcriptomes of the adult male, adult female, and nymph of X. brachycerus were sequenced. A total of 133,194,848 clean reads were obtained and de novo assembled into 43,187 unigenes with an average length of 964 bp (N50 of 1799 bp); of these, 24,717 (57.23%) unigenes matched known proteins. Based on these annotations, many putative transcripts related to X. brachycerus growth, development, environmental adaptability, and metabolism of nutritional components and bioactive components were identified. In addition, the expression profiles of all three transcriptome datasets were analyzed, and many differentially expressed genes were detected using RSEM and PossionDis. Unigenes. Unigenes with functions associated with growth and development exhibited higher transcript levels at the nymph stage, and unigenes associated with environmental adaptability showed increased transcription in the adults. These comprehensive X. brachycerus transcriptomic data will provide a useful molecular resource for gene prediction, molecular marker development, and studies on signaling pathways in this species and will serve as a reference for the efficient use of other grasshoppers.

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

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

A new A-P compartment boundary and organizer in holometabolous insect wings

Decades of research on the highly modified wings of Drosophila melanogaster has suggested that insect wings are divided into two Anterior-Posterior (A-P) compartments separated by an axis of symmetry. This axis of symmetry is created by a developmental organizer that establishes symmetrical patterns of gene expression that in turn pattern the A-P axis of the wing. Butterflies possess more typical insect wings and butterfly wing colour patterns provide many landmarks for studies of wing structure and development. Using eyespot colour pattern variation in Vanessa butterflies, here we show an additional A-P axis of symmetry running between wing sectors 3 and 4. Boundaries of Drosophila mitotic clones suggest the existence of a previously undetected Far-Posterior (F-P) compartment boundary that coincides with this additional A-P axis. A similar compartment boundary is evident in butterfly mosaic gynandromorphs. We suggest that this additional compartment boundary and its associated developmental organizer create an axis of wing colour pattern symmetry and a gene expression-based combinatorial code, permitting each insect wing compartment to acquire a unique identity and allowing for the individuation of butterfly eyespots.

spalt is functionally conserved in Locusta and Drosophila to promote wing growth

Locusta has strong fly wings to ensure its long distance migration, but the molecular mechanism that regulates the Locusta wing development is poorly understood. To address the developmental mechanism of the Locusta flying wing, we cloned the Dpp target gene spalt (sal) and analyzed its function in wing growth in the Locusta. The Locusta wing size is apparently reduced with vein defects when sal is interfered by injection of dsRNA, indicating that sal is required for locust wing growth and vein formation. This function is conserved during the Drosophila wing development. To better understand sal’s function in wing growth, we then used Drosophila wing disc as a model for further study. We found that sal promotes cell proliferation in the whole wing disc via positive regulation of a microRNA bantam. Our results firstly unravel sal’s function in the Locusta wing growth and confirm a highly conserved function of sal in Locusta and Drosophila.

Transcriptome Comparison Analysis of Ostrinia furnacalis in Four Developmental Stages

The Asian corn borer, Ostrinia furnacalis, is one of the most destructive pests of maize and causes huge losses in maize yield each year. In order to characterize the different developmental stages, a high-throughput sequencing platform was employed to perform de novo transcriptome assembly and gene expression analysis for the egg, larva, pupa and adult stages. Approximately 185 million reads were obtained, trimmed, and assembled into 42,638 unigenes with an average length of 801.94 bp and an N50 length of 1,152 bp. These unigene sequences were annotated and classified by performing Gene Ontology (GO), Cluster of Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional classifications. Comparison of the gene expression profiles of the two transitional stages revealed dramatic differences. Some differentially expressed genes are associated with digestion, cuticularization olfactory recognition and wing formation as well as growth and development. In total, 12 putative insect development-related genes were identified. Real-time quantitative PCR (RT-qPCR) results and sequencing based on relative expression levels of randomly selected genes confirmed these expression patterns. These data represent the most comprehensive transcriptomic resource currently available for O. furnacalis and will facilitate the study of developmental pathways, cuticularization, wing formation and olfactory recognition.

De novo transcriptome sequencing and analysis of Coccinella septempunctata L. in non-diapause, diapause and diapause-terminated states to identify diapause-associated genes

BACKGROUND

The most common ladybird beetle, Coccinella septempunctata L., is an excellent predator of crop pests such as aphids and white flies, and it shows a wide range of adaptability, a large appetite and a high reproductive ability. Diapause research plays an important role in the artificial propagation and shelf-life extension of insect products. Although this lady beetle's regulatory, physiological and biochemical characteristics in the diapause period are well understood, the molecular mechanism of diapause remains unknown. Therefore, we collected female adults in three different states, i.e., non-diapause, diapause and diapause termination, for transcriptome sequencing.

RESULTS

After transcriptome sequencing using the Illumina HiSeq 2500 platform with pretreatment, a total of 417.6 million clean reads from nine samples were filtered using the program FASTX (version 0.0). Additionally, 106,262 contigs were assembled into 82,820 unigenes with an average length of 921 bp and an N50 of 1,241 bp. All of the unigenes were annotated through BLASTX alignment against the Nr or UniProt database, and 37,872 unigenes were matched. We performed further analysis of these unigenes using the Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Through pairwise comparisons of the non-diapause (ND), diapause (D), and diapause-terminated (DT) groups, 3,501 and 1,427 differentially expressed genes (DEGs) were identified between D and ND and between DT and D, respectively. Moreover, 443 of the DEGs were specifically expressed during the diapause period (i.e., DEGs that were expressed at the highest or lowest levels during diapause compared with the other stages). GO function and KEGG pathway enrichment were performed on all DEGs and showed that RNA-directed DNA polymerase activity and fatty acid metabolism were significantly affected. Furthermore, eight specific expressed genes were selected for validation using qRT-PCR. Among these eight genes, seven genes were up-regulated, and one gene was down-regulated; the change trends of the eight genes were the same between the qRT-PCR and RNA-seq analysis results.

CONCLUSIONS

In this study, a new method for collecting and identifying diapause insects was described. We generated a vast quantity of transcriptome data from C. septempunctata L., providing a resource for gene function research. The diapause-associated genes that we identified establish a foundation for future studies on the molecular mechanisms of diapause.

Comparative transcriptomic analysis of immune responses of the migratory locust, Locusta migratoria, to challenge by the fungal insect pathogen, Metarhizium acridum

Background

The migratory locust, Locusta migratoria manilensis, is an immensely destructive agricultural pest that forms a devastating and voracious gregarious phase. The fungal insect pathogen, Metarhizium acridum, is a specialized locust pathogen that has been used as a potent mycoinsecticide for locust control. Little, however, is known about locust immune tissue, i.e. fat body and hemocyte, responses to challenge by this fungus.

Methods

RNA-seq (RNA sequencing) technology were applied to comparatively examine the different roles of locust fat body and hemocytes, the two major contributors to the insect immune response, in defense against M. acridum. According to the sequence identity to homologies of other species explored immune response genes, immune related unigenes were screened in all transcriptome wide range from locust and the differential expressed genes were identified in these two tissues, respectively.

Results

Analysis of differentially expressed locust genes revealed 4660 and 138 up-regulated, and 1647 and 23 down-regulated transcripts in the fat body and hemocytes, respectively after inoculation with M. acridum spores. GO (Gene Ontology) enrichment analysis showed membrane biogenesis related proteins and effector proteins significantly differentially expressed in hemocytes, while the expression of energy metabolism and development related transcripts were enriched in the fat body after fungal infection. A total of 470 immune related unigenes were identified, including members of the three major insect immune pathways, i.e. Toll, Imd (immune deficiency) and JAK/STAT (janus kinase/signal transduction and activator of transcription). Of these, 58 and three were differentially expressed in the insect fat body or hemocytes after infection, respectively. Of differential expressed transcripts post challenge, 43 were found in both the fat body and hemocytes, including the LmLys4 lysozyme, representing a microbial cell wall targeting enzyme.

Conclusions

These data indicate that locust fat body and hemocytes adopt different strategies in response to M. acridum infection. Fat body gene expression after M. acridum challenge appears to function mainly through activation of innate immune related genes, energy metabolism and development related genes. Hemocyte responses attempt to limit fungal infection primarily through regulation of membrane related genes and activation of cellular immune responses and release of humoral immune factors.

Electronic supplementary material

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

Kiwi genome provides insights into evolution of a nocturnal lifestyle

BACKGROUND

Kiwi, comprising five species from the genus Apteryx, are endangered, ground-dwelling bird species endemic to New Zealand. They are the smallest and only nocturnal representatives of the ratites. The timing of kiwi adaptation to a nocturnal niche and the genomic innovations, which shaped sensory systems and morphology to allow this adaptation, are not yet fully understood.

RESULTS

We sequenced and assembled the brown kiwi genome to 150-fold coverage and annotated the genome using kiwi transcript data and non-redundant protein information from multiple bird species. We identified evolutionary sequence changes that underlie adaptation to nocturnality and estimated the onset time of these adaptations. Several opsin genes involved in color vision are inactivated in the kiwi. We date this inactivation to the Oligocene epoch, likely after the arrival of the ancestor of modern kiwi in New Zealand. Genome comparisons between kiwi and representatives of ratites, Galloanserae, and Neoaves, including nocturnal and song birds, show diversification of kiwi's odorant receptors repertoire, which may reflect an increased reliance on olfaction rather than sight during foraging. Further, there is an enrichment of genes influencing mitochondrial function and energy expenditure among genes that are rapidly evolving specifically on the kiwi branch, which may also be linked to its nocturnal lifestyle.

CONCLUSIONS

The genomic changes in kiwi vision and olfaction are consistent with changes that are hypothesized to occur during adaptation to nocturnal lifestyle in mammals. The kiwi genome provides a valuable genomic resource for future genome-wide comparative analyses to other extinct and extant diurnal ratites.