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

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.

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.

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.

Physical linkage and mate preference generate linkage disequilibrium for behavioral isolation in two parapatric crickets

Behavioral isolation is a potent barrier to gene flow and a source of striking diversity in the animal kingdom. However, it remains unclear if the linkage disequilibrium (LD) between sex‐specific traits required for behavioral isolation results mostly from physical linkage between signal and preference loci or from directional mate preferences. Here, we test this in the field crickets Gryllus rubens and G. texensis. These closely related species diverged with gene flow and have strongly differentiated songs and preference functions for the mate calling song rhythm. We map quantitative trait loci for signal and preference traits (pQTL) as well as for gene expression associated with these traits (eQTL). We find strong, positive genetic covariance between song traits and between song and preference. Our results show that this is in part explained by incomplete physical linkage: although both linked pQTL and eQTL couple male and female traits, major effect loci for different traits were never on the same chromosome. We suggest that the finely tuned, highly divergent preference functions are likely an additional source of LD between male and female traits in this system. Furthermore, pleiotropy of gene expression presents an underappreciated mechanism to link sexually dimorphic phenotypes.

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.

Transcriptome Analysis of Cinnamomum chago: A Revelation of Candidate Genes for Abiotic Stress Response and Terpenoid and Fatty Acid Biosyntheses

Cinnamomum chago, an endangered species endemic to Yunnan province, possesses large economic and phylogenetic values in Lauraceae. However, the genomic information of this species remains relatively unexplored. In this study, we used RNAseq technology to characterize and annotate the C. chago transcriptome and identify candidate genes involved in special metabolic pathways and gene-associated simple sequence repeats (SSRs) and single-nucleotide polymorphism (SNP). A total of 129,097 unigenes, with a mean length of 667 bp and an N50 length of 1,062 bp, were assembled. Among these genes, 56,887 (44.07%) unigenes were successfully annotated using at least one database. Furthermore, 47 and 46 candidate genes were identified in terpenoid biosynthesis and fatty acid biosynthesis, respectively. A total of 22 candidate genes participated in at least one abiotic stress response of C. chago. Additionally, a total of 25,654 SSRs and 640 SNPs were also identified. Based on these potential loci, 55 novel expressed sequence tag (EST)-SSR primers were successfully developed. This work provides comprehensive transcriptomic data that can be used to establish a valuable information platform for gene prediction, signaling pathway investigation, and molecular marker development for C. chago and other related species. Such a platform can facilitate further studies on germplasm conservation and utilization of Lauraceae species.

De novo Assembly of the Burying Beetle Nicrophorus orbicollis (Coleoptera: Silphidae) Transcriptome Across Developmental Stages with Identification of Key Immune Transcripts

Burying beetles (Nicrophorus spp.) are among the relatively few insects that provide parental care while not belonging to the eusocial insects such as ants or bees. This behavior incurs energy costs as evidenced by immune deficits and shorter life-spans in reproducing beetles. In the absence of an assembled transcriptome, relatively little is known concerning the molecular biology of these beetles. This work details the assembly and analysis of the Nicrophorus orbicollis transcriptome at multiple developmental stages. RNA-Seq reads were obtained by next-generation sequencing and the transcriptome was assembled using the Trinity assembler. Validation of the assembly was performed by functional characterization using Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Differential expression analysis highlights developmental stage-specific expression patterns, and immunity-related transcripts are discussed. The data presented provides a valuable molecular resource to aid further investigation into immunocompetence throughout this organism's sexual development.

Demography and selection shape transcriptomic divergence in field crickets

Gene flow, demography, and selection can result in similar patterns of genomic variation and disentangling their effects is key to understanding speciation. Here, we assess transcriptomic variation to unravel the evolutionary history of Gryllus rubens and Gryllus texensis, cryptic field cricket species with highly divergent mating behavior. We infer their demographic history and screen their transcriptomes for footprints of selection in the context of the inferred demography. We find strong support for a long history of bidirectional gene flow, which ceased during the late Pleistocene, and a bottleneck in G. rubens consistent with a peripatric origin of this species. Importantly, the demographic history has likely strongly shaped patterns of genetic differentiation (empirical F_ST distribution). Concordantly, F_ST -based selection detection uncovers a large number of outliers, likely comprising many false positives, echoing recent theoretical insights. Alternative genetic signatures of positive selection, informed by the demographic history of the sibling species, highlighted a smaller set of loci; many of these are candidates for controlling variation in mating behavior. Our results underscore the importance of demography in shaping overall patterns of genetic divergence and highlight that examining both demography and selection facilitates a more complete understanding of genetic divergence during speciation.

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.