ASGARD: an open-access database of annotated transcriptomes for emerging model arthropod species

Proteomic analysis of Sarcoptes scabiei reveals that proteins differentially expressed between eggs and female adult stages are involved predominantly in genetic information processing, metabolism and/or host-parasite interactions

Presently, there is a dearth of proteomic data for parasitic mites and their relationship with the host animals. Here, using a high throughput LC-MS/MS-based approach, we undertook the first comprehensive, large-scale proteomic investigation of egg and adult female stages of the scabies mite, Sarcoptes scabiei-one of the most important parasitic mites of humans and other animals worldwide. In total, 1,761 S. scabiei proteins were identified and quantified with high confidence. Bioinformatic analyses revealed differentially expressed proteins to be involved predominantly in biological pathways or processes including genetic information processing, energy (oxidative phosphorylation), nucleotide, amino acid, carbohydrate and/or lipid metabolism, and some adaptive processes. Selected, constitutively and highly expressed proteins, such as peptidases, scabies mite inactivated protease paralogues (SMIPPs) and muscle proteins (myosin and troponin), are proposed to be involved in key biological processes within S. scabiei, host-parasite interactions and/or the pathogenesis of scabies. These proteomic data will enable future molecular, biochemical and physiological investigations of early developmental stages of S. scabiei and the discovery of novel interventions, targeting the egg stage, given its non-susceptibility to acaricides currently approved for the treatment of scabies in humans.

CerealESTDb: A Comprehensive Resource for Abiotic Stress-Responsive Annotated ESTs With Predicted Genes, Gene Ontology, and Metabolic Pathways in Major Cereal Crops

Cereals are the most important food crops and are considered key contributors to global food security. Loss due to abiotic stresses in cereal crops is limiting potential productivity in a significant manner. The primary reasons for abiotic stresses are abrupt temperature, variable rainfall, and declining nutrient status of the soil. Varietal development is the key to sustaining productivity under influence of multiple abiotic stresses and must be studied in context with genomics and molecular breeding. Recently, advances in a plethora of Next Generation Sequencing (NGS) based methods have accelerated the enormous genomic data generation associated with stress-induced transcripts such as microarray, RNAseq, Expressed Sequenced Tag (ESTs), etc. Many databases related to microarray and RNA-seq based transcripts have been developed and profusely utilized. However, an abundant amount of transcripts related to abiotic stresses in various cereal crops arising from EST technology are available but still remain underutilized in absence of a consolidated database. In this study, an attempt has been made with a primary goal to integrate, analyse, and characterise the available resources of ESTs responsive to abiotic stresses in major cereals. The developed CerealESTdb presents a customisable search in two different ways in the form of searchable content for easy access and potential use. This database comprises ESTs from four major cereal crops, namely rice (Oryza sativa L.), wheat (Triticum aestivum L.), sorghum (Sorghum bicolour L.), and maize (Zea mays L.), under a set of abiotic stresses. The current statistics of this cohesive database consists of 55,826 assembled EST sequences, 51,791 predicted genes models, and their 254,609 gene ontology terms including extensive information on 1,746 associated metabolic pathways. We anticipate that developed CerealESTdb will be helpful in deciphering the knowledge of complex biological phenomena under abiotic stresses to accelerate the molecular breeding programs towards the development of crop cultivars resilient to abiotic stresses. The CerealESTdb is publically available with the URL http://cabgrid.res.in/CerealESTDb.

InsectBase 2.0: a comprehensive gene resource for insects

Insects are the largest group of animals on the planet and have a huge impact on human life by providing resources, transmitting diseases, and damaging agricultural crop production. Recently, a large amount of insect genome and gene data has been generated. A comprehensive database is highly desirable for managing, sharing, and mining these resources. Here, we present an updated database, InsectBase 2.0 (http://v2.insect-genome.com/), covering 815 insect genomes, 25 805 transcriptomes and >16 million genes, including 15 045 111 coding sequences, 3 436 022 3'UTRs, 4 345 664 5'UTRs, 112 162 miRNAs and 1 293 430 lncRNAs. In addition, we used an in-house standard pipeline to annotate 1 434 653 genes belonging to 164 gene families; 215 986 potential horizontally transferred genes; and 419 KEGG pathways. Web services such as BLAST, JBrowse2 and Synteny Viewer are provided for searching and visualization. InsectBase 2.0 serves as a valuable platform for entomologists and researchers in the related communities of animal evolution and invertebrate comparative genomics.

A conserved role for arrow in posterior axis patterning across Arthropoda

Segmentation is a key characteristic of Arthropoda that is linked to the evolutionary success of this lineage. It has previously been shown in both vertebrates and short germ insects that posterior segmentation requires canonical Wnt (cWnt) signaling, which maintains the expression of Caudal and the posterior growth zone; disruption of cWnt signaling incurs posterior truncations in these lineages due to the loss of the tail bud. However, comparable datasets for Wnt signaling are limited outside of holometabolous insects, due to incomparable phenotypic spectra and inefficacy of gene misexpression methods in certain model species. We applied RNA interference (RNAi) against the Wnt co-receptor arrow (arr), a key member of the cWnt signaling pathway in holometabolous insects and vertebrates, to examine posterior axis elongation of the cobweb spider Parasteatoda tepidariorum (short germ embryogenesis; one Wnt8 homolog), the cricket Gryllus bimaculatus (intermediate germ; one Wnt8 homolog), and the milkweed bug Oncopeltus fasciatus (short germ; two Wnt8 homologs). Knockdown of arr in insects resulted in posterior truncations affecting the gnathos through the abdomen in O. fasciatus, whereas posterior truncations only affected the T3 segment through the abdomen in G. bimaculatus. Spider embryos with disrupted arr expression exhibited defects along the entire axis, including segmentation defects throughout the germband. RNA-Seq-based differential gene expression analysis of severe Ptep-arr loss-of-function phenotypes at two developmental stages was used to confirm that knockdown of Ptep-arr results in systemic disruption of the Wnt pathway. Intriguingly, we found that knockdown of arr did not abrogate Wnt8 expression in any of the three species, with cad expression additionally retained in severe loss-of-function phenotypes in the cricket and the spider. Together with data from a holometabolous insect, our results suggest that cWnt signaling is not required for maintenance of Wnt8 expression across Arthropoda. These outcomes underscore the diagnostic power of differential gene expression analyses in characterizing catastrophic phenotypes in emerging model species.

The common house spider Parasteatoda tepidariorum

The common house spider Parasteatoda tepidariorum, belonging to the Chelicerata in the phylum Arthropoda, has emerged as an experimental system for studying mechanisms of development from an evolutionary standpoint. In this article, we review the distinct characteristics of P. tepidariorum, the major research questions relevant to this organism, and the available key methods and resources. P. tepidariorum has a relatively short lifecycle and, once mated, periodically lays eggs. The morphogenetic field of the P. tepidariorum embryo is cellular from an early stage and exhibits stepwise symmetry-breaking events and stripe-forming processes that are associated with body axes formation and segmentation, respectively, before reaching the arthropod phylotypic stage. Self-regulatory capabilities of the embryonic field are a prominent feature in P. tepidariorum. The mechanisms and logic underlying the evolvability of heritable patterning systems at the phylum level could be one of the major avenues of research investigated using this animal. The sequenced genome reveals whole genome duplication (WGD) within chelicerates, which offers an invertebrate platform for investigating the potential roles of WGD in animal diversification and evolution. The development and evolution of lineage-specific organs, including the book lungs and the union of spinnerets and silk glands, are attractive subjects of study. Studies using P. tepidariorum can benefit from the use of parental RNA interference, microinjection applications (including cell labeling and embryonic RNA interference), multicolor fluorescence in situ hybridization, and laser ablation as well as rich genomic and transcriptomic resources. These techniques enable functional gene discoveries and the uncovering of cellular and molecular insights.

The amphipod crustacean Parhyale hawaiensis: An emerging comparative model of arthropod development, evolution, and regeneration

Recent advances in genetic manipulation and genome sequencing have paved the way for a new generation of research organisms. The amphipod crustacean Parhyale hawaiensis is one such system. Parhyale are easy to rear and offer large broods of embryos amenable to injection, dissection, and live imaging. Foundational work has described Parhyale embryonic development, while advancements in genetic manipulation using CRISPR-Cas9 and other techniques, combined with genome and transcriptome sequencing, have enabled its use in studies of arthropod development, evolution, and regeneration. This study introduces Parhyale development and life history, a catalog of techniques and resources for Parhyale research, and two case studies illustrating its power as a comparative research system. This article is categorized under: Comparative Development and Evolution > Evolutionary Novelties Adult Stem Cells, Tissue Renewal, and Regeneration > Regeneration Comparative Development and Evolution > Model Systems Comparative Development and Evolution > Body Plan Evolution.

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.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

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

Non-insect crustacean models in developmental genetics including an encomium to Parhyale hawaiensis

The impressive diversity of body plans, lifestyles and segmental specializations exhibited by crustaceans (barnacles, copepods, shrimps, crabs, lobsters and their kin) provides great material to address longstanding questions in evolutionary developmental biology. Recent advances in forward and reverse genetics and in imaging approaches applied in the amphipod Parhyale hawaiensis and other emerging crustacean model species have made it possible to probe the molecular and cellular basis of crustacean diversity. A number of biological and technical qualities like the slow tempo and holoblastic cleavage mode, the stereotypy of many cellular processes, the functional and morphological diversity of limbs along the body axis, and the availability of various experimental manipulations, have made Parhyale a powerful system to study normal development and regeneration.

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.

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.

Ablation of a single cell from eight-cell embryos of the amphipod crustacean Parhyale hawaiensis

The amphipod Parhyale hawaiensis is a small crustacean found in intertidal marine habitats worldwide. Over the past decade, Parhyale has emerged as a promising model organism for laboratory studies of development, providing a useful outgroup comparison to the well studied arthropod model organism Drosophila melanogaster. In contrast to the syncytial cleavages of Drosophila, the early cleavages of Parhyale are holoblastic. Fate mapping using tracer dyes injected into early blastomeres have shown that all three germ layers and the germ line are established by the eight-cell stage. At this stage, three blastomeres are fated to give rise to the ectoderm, three are fated to give rise to the mesoderm, and the remaining two blastomeres are the precursors of the endoderm and germ line respectively. However, blastomere ablation experiments have shown that Parhyale embryos also possess significant regulatory capabilities, such that the fates of blastomeres ablated at the eight-cell stage can be taken over by the descendants of some of the remaining blastomeres. Blastomere ablation has previously been described by one of two methods: injection and subsequent activation of phototoxic dyes or manual ablation. However, photoablation kills blastomeres but does not remove the dead cell body from the embryo. Complete physical removal of specific blastomeres may therefore be a preferred method of ablation for some applications. Here we present a protocol for manual removal of single blastomeres from the eight-cell stage of Parhyale embryos, illustrating the instruments and manual procedures necessary for complete removal of the cell body while keeping the remaining blastomeres alive and intact. This protocol can be applied to any Parhyale cell at the eight-cell stage, or to blastomeres of other early cleavage stages. In addition, in principle this protocol could be applicable to early cleavage stage embryos of other holoblastically cleaving marine invertebrates.

Subdivision of arthropod cap-n-collar expression domains is restricted to Mandibulata

BACKGROUND

The monophyly of Mandibulata - the division of arthropods uniting pancrustaceans and myriapods - is consistent with several morphological characters, such as the presence of sensory appendages called antennae and the eponymous biting appendage, the mandible. Functional studies have demonstrated that the patterning of the mandible requires the activity of the Hox gene Deformed and the transcription factor cap-n-collar (cnc) in at least two holometabolous insects: the fruit fly Drosophila melanogaster and the beetle Tribolium castaneum. Expression patterns of cnc from two non-holometabolous insects and a millipede have suggested conservation of the labral and mandibular domains within Mandibulata. However, the activity of cnc is unknown in crustaceans and chelicerates, precluding understanding of a complete scenario for the evolution of patterning of this appendage within arthropods. To redress these lacunae, here we investigate the gene expression of the ortholog of cnc in Parhyale hawaiensis, a malacostracan crustacean, and two chelicerates: the harvestman Phalangium opilio, and the scorpion Centruroides sculpturatus.

RESULTS

In the crustacean P. hawaiensis, the segmental expression of Ph-cnc is the same as that reported previously in hexapods and myriapods, with two distinct head domains in the labrum and the mandibular segment. In contrast, Po-cnc and Cs-cnc expression is not enriched in the labrum of either chelicerate, but instead is expressed at comparable levels in all appendages. In further contrast to mandibulate orthologs, the expression domain of Po-cnc posterior to the labrum is not confined within the expression domain of Po-Dfd.

CONCLUSIONS

Expression data from two chelicerate outgroup taxa suggest that the signature two-domain head expression pattern of cnc evolved at the base of Mandibulata. The observation of the archetypal labral and mandibular segment domains in a crustacean exemplar supports the synapomorphic nature of mandibulate cnc expression. The broader expression of Po-cnc with respect to Po-Dfd in chelicerates further suggests that the regulation of cnc by Dfd was also acquired at the base of Mandibulata. To test this hypothesis, future studies examining panarthropod cnc evolution should investigate expression of the cnc ortholog in arthropod outgroups, such as Onychophora and Tardigrada.

Identification of a putative germ plasm in the amphipod Parhyale hawaiensis

Background

Specification of the germ line is an essential event during the embryonic development of sexually reproducing animals, as germ line cells are uniquely capable of giving rise to the next generation. Animal germ cells arise through either inheritance of a specialized, maternally supplied cytoplasm called 'germ plasm’ or though inductive signaling by somatic cells. Our understanding of germ cell determination is based largely on a small number of model organisms. To better understand the evolution of germ cell specification, we are investigating this process in the amphipod crustacean Parhyale hawaiensis. Experimental evidence from previous studies demonstrated that Parhyale germ cells are specified through inheritance of a maternally supplied cytoplasmic determinant; however, this determinant has not been identified.

Results

Here we show that the one-cell stage Parhyale embryo has a distinct cytoplasmic region that can be identified by morphology as well as the localization of germ line-associated RNAs. Removal of this cytoplasmic region results in a loss of embryonic germ cells, supporting the hypothesis that it is required for specification of the germ line. Surprisingly, we found that removal of this distinct cytoplasm also results in aberrant somatic cell behaviors, as embryos fail to gastrulate.

Conclusions

Parhyale hawaiensis embryos have a specialized cytoplasm that is required for specification of the germ line. Our data provide the first functional evidence of a putative germ plasm in a crustacean and provide the basis for comparative functional analysis of germ plasm formation within non-insect arthropods.

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

Most genomic resources available for insects represent the Holometabola, which are insects that undergo complete metamorphosis like beetles and flies. In contrast, the Hemimetabola (direct developing insects), representing the basal branches of the insect tree, have very few genomic resources. We have therefore created a large and publicly available transcriptome for the hemimetabolous insect Gryllus bimaculatus (cricket), a well-developed laboratory model organism whose potential for functional genetic experiments is currently limited by the absence of genomic resources. cDNA was prepared using mRNA obtained from adult ovaries containing all stages of oogenesis, and from embryo samples on each day of embryogenesis. Using 454 Titanium pyrosequencing, we sequenced over four million raw reads, and assembled them into 21,512 isotigs (predicted transcripts) and 120,805 singletons with an average coverage per base pair of 51.3. We annotated the transcriptome manually for over 400 conserved genes involved in embryonic patterning, gametogenesis, and signaling pathways. BLAST comparison of the transcriptome against the NCBI non-redundant protein database (nr) identified significant similarity to nr sequences for 55.5% of transcriptome sequences, and suggested that the transcriptome may contain 19,874 unique transcripts. For predicted transcripts without significant similarity to known sequences, we assessed their similarity to other orthopteran sequences, and determined that these transcripts contain recognizable protein domains, largely of unknown function. We created a searchable, web-based database to allow public access to all raw, assembled and annotated data. This database is to our knowledge the largest de novo assembled and annotated transcriptome resource available for any hemimetabolous insect. We therefore anticipate that these data will contribute significantly to more effective and higher-throughput deployment of molecular analysis tools in Gryllus.

Evidence against a germ plasm in the milkweed bug Oncopeltus fasciatus, a hemimetabolous insect

Primordial germ cell (PGC) formation in holometabolous insects like Drosophila melanogaster relies on maternally synthesised germ cell determinants that are asymmetrically localised to the oocyte posterior cortex. Embryonic nuclei that inherit this "germ plasm" acquire PGC fate. In contrast, historical studies of basally branching insects (Hemimetabola) suggest that a maternal requirement for germ line genes in PGC specification may be a derived character confined principally to Holometabola. However, there have been remarkably few investigations of germ line gene expression and function in hemimetabolous insects. Here we characterise PGC formation in the milkweed bug Oncopeltus fasciatus, a member of the sister group to Holometabola, thus providing an important evolutionary comparison to members of this clade. We examine the transcript distribution of orthologues of 19 Drosophila germ cell and/or germ plasm marker genes, and show that none of them localise asymmetrically within Oncopeltus oocytes or early embryos. Using multiple molecular and cytological criteria, we provide evidence that PGCs form after cellularisation at the site of gastrulation. Functional studies of vasa and tudor reveal that these genes are not required for germ cell formation, but that vasa is required in adult males for spermatogenesis. Taken together, our results provide evidence that Oncopeltus germ cells may form in the absence of germ plasm, consistent with the hypothesis that germ plasm is a derived strategy of germ cell specification in insects.