Tissue-specific transcriptomics in the field cricket Teleogryllus oceanicus

Genome-wide identification and expression analysis of the mating-responsive genes in the male accessory glands of Spodoptera litura (Lepidoptera: Noctuidae)

BACKGROUND

Mating elicits significant changes in gene expression and leads to subsequent physiological and behavioural modifications in insects. The reproductive success of both sexes is contributed immensely by the male accessory gland (MAG) proteins that are transferred along with sperms to the female reproductive tract during mating where they facilitate several processes that modify the post-mating behaviour. The mating-responsive genes in the MAGs have been identified and reported in many insects but have not been well-characterized in the important agricultural pest Spodoptera litura. Here, we present RNA sequencing analysis to identify mating-responsive genes from the accessory glands of virgin males and males interrupted during mating.

RESULTS

Overall, 91,744 unigenes were generated after clustering the assembled transcript sequences of both samples, while the total number of transcripts annotated was 48,708 based on sequence homology against the non-redundant (NR) database. Comparative transcriptomics analysis revealed 16,969 genes that were differentially expressed between the two groups, including 9814 up-regulated and 7155 down-regulated genes. Among the top 80 genes that were selected for heat map analysis, several prominent genes including odorant binding protein, cytochrome P450, heat shock proteins, juvenile hormone binding protein, carboxypeptidases and serine protease were differentially expressed.

CONCLUSIONS

The identified genes are known or predicted to promote several processes that modify the female post-mating behaviour. Future studies with the individual MAG protein or in combination will be required to recognize the precise mechanisms by which these proteins alter female physiology and reproductive behaviour. Thus, our study provides essential data to address fundamental questions about reproduction within and among insects and also paves way for further exploration of the functions of these proteins in female insects.

Genomics and genome editing techniques of crickets, an emerging model insect for biology and food science

Most tools available for manipulating gene function in insects have been developed for holometabolous species. In contrast, functional genetics tools for the Hemimetabola are highly underdeveloped. This is a barrier both to understanding ancestral insect biology, and to optimizing contemporary study and manipulation of particular large hemimetabolous orders of crucial economic and agricultural importance like the Orthoptera. For orthopteran insects, including crickets, the rapid spread of next-generation sequencing technology has made transcriptome data available for a wide variety of species over the past decade. Furthermore, whole genome sequences of orthopteran insects with relatively large genome sizes are now available. With these new genome assemblies and the development of genome editing technologies such as the CRISPR-Cas9 system, it has become possible to create gene knock-out and knock-in strains in orthopteran insects. As a result, orthopteran species should become increasingly feasible for laboratory study not only in research fields that have traditionally used insects, but also in agricultural fields that use them as food and feed. In this review, we summarize these recent advances and their relevance to such applications.

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.

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.

Socially cued seminal fluid gene expression mediates responses in ejaculate quality to sperm competition risk

There is considerable evidence that males will increase the number of sperm ejaculated in response to sperm competition risk. However, whether they have the capacity to adjust seminal fluid components of the ejaculate has received less attention. Male crickets (Teleogryllus oceanicus) have been shown to adjust the viability of sperm in their ejaculate in response to sperm competition risk. Here we show that socially mediated plasticity in sperm viability is probably due, at least in part, to male adjustments in the protein composition of the seminal fluid. Seven seminal fluid protein genes were found to have an increased expression in males exposed to rival calls. Increased expression of these genes was correlated with increased sperm viability in whole ejaculates, and gene knockdown confirmed that at least one of these proteins promotes sperm viability. Our results lend support for recent theoretical models that predict complex responses in male allocation to seminal fluid composition in response to sperm competition risk.

Genes Integral to the Reproductive Function of Male Reproductive Tissues Drive Heterogeneity in Evolutionary Rates in Japanese Quail

Early comparative genomics studies originally uncovered a nonintuitive pattern; genes involved in reproduction appeared to evolve more rapidly than other classes of genes. Currently, the emerging consensus is that genes encoding reproductive proteins evolve under variable selective pressures, producing more heterogeneous divergence patterns than previously appreciated. Here, we investigate a facet of that heterogeneity and explore the factors that drive male reproductive tissue-based heterogeneity in evolutionary rates. In Japanese quail (Coturnix japonica), genes with enriched expression in the testes evolve much more rapidly than those enriched in the foam gland (FG), a novel gland that secretes an airy foam that males transfer to females during mating. We compared molecular evolutionary patterns among (1) genes with induced expression in breeding vs. wintering conditions for both tissues and (2) genes that encode foam proteins (FPs) vs. those with varying degrees of expression specificity in the FG. We report two major findings. First, genes upregulated in breeding condition testes evolve exceptionally rapidly, while those induced in breeding condition FGs evolve slowly. These differences hold even after correcting for hormonally-dependent gene expression and chromosomal location. Second, genes encoding FPs are extremely conserved in terms of gene identity and sequence. Together, these finding suggest that genes involved in the reproductive function of each tissue drive the marked rate of heterogeneity.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

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.

Australian black field crickets show changes in neural gene expression associated with socially-induced morphological, life-history, and behavioral plasticity

BACKGROUND

Ecological and evolutionary model organisms have provided extensive insight into the ecological triggers, adaptive benefits, and evolution of life-history driven developmental plasticity. Despite this, we still have a poor understanding of the underlying genetic changes that occur during shifts towards different developmental trajectories. The goal of this study is to determine whether we can identify underlying gene expression patterns that can describe the different life-history trajectories individuals follow in response to social cues of competition. To do this, we use the Australian black field cricket (Teleogryllus commodus), a species with sex-specific developmental trajectories moderated by the density and quality of calls heard during immaturity. In this study, we manipulated the social information males and females could hear by rearing individuals in either calling or silent treatments. We next used RNA-Seq to develop a reference transcriptome to study changes in brain gene expression at two points prior to sexual maturation.

RESULTS

We show accelerated development in both sexes when exposed to calling; changes were also seen in growth, lifespan, and reproductive effort. Functional relationships between genes and phenotypes were apparent from ontological enrichment analysis. We demonstrate that increased investment towards traits such as growth and reproductive effort were often associated with the expression of a greater number of genes with similar effect, thus providing a suite of candidate genes for future research in this and other invertebrate organisms.

CONCLUSIONS

Our results provide interesting insight into the genomic underpinnings of developmental plasticity and highlight the potential of a genomic exploration of other evolutionary theories such as condition dependence and sex-specific developmental strategies.

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.

Starvation resistance and tissue-specific gene expression of stress-related genes in a naturally inbred ant population

Starvation is one of the most common and severe stressors in nature. Not only does it lead to death if not alleviated, it also forces the starved individual to allocate resources only to the most essential processes. This creates energetic trade-offs which can lead to many secondary challenges for the individual. These energetic trade-offs could be exacerbated in inbred individuals, which have been suggested to have a less efficient metabolism. Here, we studied the effect of inbreeding on starvation resistance in a natural population of Formica exsecta ants, with a focus on survival and tissue-specific expression of stress, metabolism and immunity-related genes. Starvation led to large tissue-specific changes in gene expression, but inbreeding had little effect on most of the genes studied. Our results illustrate the importance of studying stress responses in different tissues instead of entire organisms.

Sexual selection and population divergence I: The influence of socially flexible cuticular hydrocarbon expression in male field crickets (Teleogryllus oceanicus)

Debates about how coevolution of sexual traits and preferences might promote evolutionary diversification have permeated speciation research for over a century. Recent work demonstrates that the expression of such traits can be sensitive to variation in the social environment. Here, we examined social flexibility in a sexually selected male trait-cuticular hydrocarbon (CHC) profiles-in the field cricket Teleogryllus oceanicus and tested whether population genetic divergence predicts the extent or direction of social flexibility in allopatric populations. We manipulated male crickets' social environments during rearing and then characterized CHC profiles. CHC signatures varied considerably across populations and also in response to the social environment, but our prediction that increased social flexibility would be selected in more recently founded populations exposed to fluctuating demographic environments was unsupported. Furthermore, models examining the influence of drift and selection failed to support a role of sexual selection in driving population divergence in CHC profiles. Variation in social environments might alter the dynamics of sexual selection, but our results align with theoretical predictions that the role social flexibility plays in modulating evolutionary divergence depends critically on whether responses to variation in the social environment are homogeneous across populations, or whether gene by social environment interactions occur.

Proteome analysis of male accessory gland secretions in oriental fruit flies reveals juvenile hormone-binding protein, suggesting impact on female reproduction

In insects, the accessory gland proteins (Acps) secreted by male accessory glands (MAGs) account for the majority of seminal fluids proteins. Mixed with sperm, they are transferred to the female at mating and so impact reproduction. In this project, we identified 2,927 proteins in the MAG secretions of the oriental fruit fly Bactrocera dorsalis, an important agricultural pest worldwide, using LC-MS analysis, and all sequences containing open reading frames were analyzed using signalP. In total, 90 Acps were identified. About one third (26) of these 90 Acps had a specific functional description, while the other two thirds (64) had no functional description including dozens of new classes of proteins. Hence, several of these novel Acps were abundant in the MAG secretions, and we confirmed their MAG-specific expression by qPCR. Finally and interestingly, one of these novel proteins was functionally predicted as juvenile hormone-binding protein, suggesting the impact of Acps with reproductive events in the female. Our results will aid in the development of an experimental method to identify Acps in insects, and in turn this information with new Acps in B. dorsalis will pave the way of further exploration their function in reproduction and potential development as new insecticide targets.

What's in the Gift? Towards a Molecular Dissection of Nuptial Feeding in a Cricket

Nuptial gifts produced by males and transferred to females during copulation are common in insects. Yet, their precise composition and subsequent physiological effects on the female recipient remain unresolved. Male decorated crickets Gryllodes sigillatus transfer a spermatophore to the female during copulation that is composed of an edible gift, the spermatophylax, and the ampulla that contains the ejaculate. After transfer of the spermatophore, the female detaches the spermatophylax and starts to eat it while sperm from the ampulla are evacuated into the female reproductive tract. When the female has finished consuming the spermatophylax, she detaches the ampulla and terminates sperm transfer. Hence, one simple function of the spermatophylax is to ensure complete sperm transfer by distracting the female from prematurely removing the ampulla. However, the majority of orally active components of the spermatophylax itself and their subsequent effects on female behavior have not been identified. Here, we report the first analysis of the proteome of the G. sigillatus spermatophylax and the transcriptome of the male accessory glands that make these proteins. The accessory gland transcriptome was assembled into 17,691 transcripts whilst about 30 proteins were detected within the mature spermatophylax itself. Of these 30 proteins, 18 were encoded by accessory gland encoded messages. Most spermatophylax proteins show no similarity to proteins with known biological functions and are therefore largely novel. A spermatophylax protein shows similarity to protease inhibitors suggesting that it may protect the biologically active components from digestion within the gut of the female recipient. Another protein shares similarity with previously characterized insect polypeptide growth factors suggesting that it may play a role in altering female reproductive physiology concurrent with fertilization. Characterization of the spermatophylax proteome provides the first step in identifying the genes encoding these proteins in males and in understanding their biological functions in the female recipient.

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.

eulerAPE: drawing area-proportional 3-Venn diagrams using ellipses

Venn diagrams with three curves are used extensively in various medical and scientific disciplines to visualize relationships between data sets and facilitate data analysis. The area of the regions formed by the overlapping curves is often directly proportional to the cardinality of the depicted set relation or any other related quantitative data. Drawing these diagrams manually is difficult and current automatic drawing methods do not always produce appropriate diagrams. Most methods depict the data sets as circles, as they perceptually pop out as complete distinct objects due to their smoothness and regularity. However, circles cannot draw accurate diagrams for most 3-set data and so the generated diagrams often have misleading region areas. Other methods use polygons to draw accurate diagrams. However, polygons are non-smooth and non-symmetric, so the curves are not easily distinguishable and the diagrams are difficult to comprehend. Ellipses are more flexible than circles and are similarly smooth, but none of the current automatic drawing methods use ellipses. We present eulerAPE as the first method and software that uses ellipses for automatically drawing accurate area-proportional Venn diagrams for 3-set data. We describe the drawing method adopted by eulerAPE and we discuss our evaluation of the effectiveness of eulerAPE and ellipses for drawing random 3-set data. We compare eulerAPE and various other methods that are currently available and we discuss differences between their generated diagrams in terms of accuracy and ease of understanding for real world data.

Genomics: moving behavioural ecology beyond the phenotypic gambit

Researchers studying the adaptive significance of behaviour typically assume that genetic mechanisms will not inhibit evolutionary trajectories, an assumption commonly known as the 'phenotypic gambit'. Although the phenotypic gambit continues to be a useful heuristic for behavioural ecology, here we discuss how genomic methods provide new tools and conceptual approaches that are relevant to behavioural ecology. We first describe how the concept of a genetic toolkit for behaviour can allow behavioural ecologists to synthesize both genomic and ecological information when assessing behavioural adaptation. Then we show how gene expression profiles can be viewed as complex phenotypic measurements, used to (1) predict behaviour, (2) evaluate phenotypic plasticity and (3) devise methods to manipulate behaviour in order to test adaptive hypotheses. We propose that advances in genomics and bioinformatics may allow researchers to overcome some of the logistical obstacles that motivated the inception of the phenotypic gambit. Behavioural ecology and genomics are mutually informative, providing potential synergy that could lead to powerful advances in the field of animal behaviour.

Insights into the development and evolution of exaggerated traits using de novo transcriptomes of two species of horned scarab beetles

Scarab beetles exhibit an astonishing variety of rigid exo-skeletal outgrowths, known as “horns”. These traits are often sexually dimorphic and vary dramatically across species in size, shape, location, and allometry with body size. In many species, the horn exhibits disproportionate growth resulting in an exaggerated allometric relationship with body size, as compared to other traits, such as wings, that grow proportionately with body size. Depending on the species, the smallest males either do not produce a horn at all, or they produce a disproportionately small horn for their body size. While the diversity of horn shapes and their behavioural ecology have been reasonably well studied, we know far less about the proximate mechanisms that regulate horn growth. Thus, using 454 pyrosequencing, we generated transcriptome profiles, during horn growth and development, in two different scarab beetle species: the Asian rhinoceros beetle, Trypoxylus dichotomus, and the dung beetle, Onthophagus nigriventris. We obtained over half a million reads for each species that were assembled into over 6,000 and 16,000 contigs respectively. We combined these data with previously published studies to look for signatures of molecular evolution. We found a small subset of genes with horn-biased expression showing evidence for recent positive selection, as is expected with sexual selection on horn size. We also found evidence of relaxed selection present in genes that demonstrated biased expression between horned and horn-less morphs, consistent with the theory of developmental decoupling of phenotypically plastic traits.