Identification of cold-responsive genes in a New Zealand alpine stick insect using RNA-Seq

Comparative transcriptome analysis of false codling moth, Thaumatotibia leucotreta in response to high and low-temperature treatments

The false codling moth (FCM), Thaumatotibia leucotreta, is a major quarantine pest native to Africa. Physical postharvest phytosanitary measures such as cold and heat treatments are championed to control its spread to new regions. However, the molecular changes that T. leucotreta undergoes as it attempts to adjust to its surroundings during the treatments and withstand the extreme temperatures remain largely unknown. The current study employs RNA-seq using the next-generation Illumina HiSeq platform to produce transcriptome profiles for differential gene expression analysis of T. leucotreta larvae under thermal stress. The transcriptome assembly analysis revealed 226,067 transcripts, clustering into 127,018 unigenes. In comparison to the 25 °C treated group, 874, 91, 159, and 754 individual differentially expressed genes (DEGs) co-regulated at -10, 0, 40, and 50 °C, respectively were discovered. Annotation of the DEGs by gene ontology (GO) revealed several genes, previously implicated in low and high-temperature stresses, including heat shock proteins, cytochrome P450, cuticle proteins, odorant binding proteins, and immune system genes. Kyoto Encyclopedia of Genes and Genomics (KEGG) classification analysis revealed that substantive DEGs were those involved in metabolic pathways such as thiamine, purine, folate, and glycerolipid metabolism pathways. The RT-qPCR validation of several significantly up- and down-regulated DEGs showed congruence between RNA-seq and qPCR data. This baseline study lays a foundation for future research into the molecular mechanisms underlying T. leucotreta's cold/heat tolerance by providing a thorough differential gene expression analysis that has identified multiple genes that may be associated with the insect's ability to withstand cold and heat.

Candidate genes for domestication and resistance to cold climate according to whole genome sequencing data of Russian cattle and sheep breeds

It is known that different species of animals, when living in the same environmental conditions, can form similar phenotypes. The study of the convergent evolution of several species under the influence of the same environmental factor makes it possible to identify common mechanisms of genetic adaptation. Local cattle and sheep breeds have been formed over thousands of years under the influence of domestication, as well as selection aimed at adaptation to the local environment and meeting human needs. Previously, we identified a number of candidate genes in genome regions potentially selected during domestication and adaptation to the climatic conditions of Russia, in local breeds of cattle and sheep using whole genome genotyping data. However, these data are of low resolution and do not reveal most nucleotide substitutions. The aim of the work was to create, using the whole genome sequencing data, a list of genes associated with domestication, selection and adaptation in Russian cattle and sheep breeds, as well as to identify candidate genes and metabolic pathways for selection for cold adaptation. We used our original data on the search for signatures of selection in the genomes of Russian cattle (Yakut, Kholmogory, Buryat, Wagyu) and sheep (Baikal, Tuva) breeds. We used the HapFLK, DCMS, FST and PBS methods to identify DNA regions with signatures of selection. The number of candidate genes in potentially selective regions was 946 in cattle and 151 in sheep. We showed that the studied Russian cattle and sheep breeds have at least 10 genes in common, apparently involved in the processes of adaptation/selection, including adaptation to a cold climate, including the ASTN2, PM20D1, TMEM176A, and GLIS1 genes. Based on the intersection with the list of selected genes in at least two Arctic/Antarctic mammal species, 20 and 8 genes, have been identified in cattle and sheep, respectively, that are potentially involved in cold adaptation. Among them, the most promising for further research are the ASPH, NCKAP5L, SERPINF1, and SND1 genes. Gene ontology analysis indicated the existence of possible common biochemical pathways for adaptation to cold in domestic and wild mammals associated with cytoskeleton disassembly and apoptosis.

Comparative transcriptome analysis of Liriomyza trifolii (Burgess) and Liriomyza sativae (Blanchard) (Diptera: Agromyzidae) in response to rapid cold hardening

The ability of insets to react efficiently to fluctuation in temperature is crucial for them to survive in variable surroundings. Rapid cold hardening (RCH) is a process that increase cold tolerance in most insect species. The molecular mechanisms of RCH remain largely unknown, and whether it is associated with transcriptional changes is unclear. In this study, we compared the transcriptomes of Liriomyza trifolii and L. sativae exposed to RCH to investigate the transcript abundance due to RCH in both species. RNA-seq revealed 93,166 assembled unigenes, and 34,303 of these were annotated in the L. trifolii and L. sativae transcriptome libraries. After a 4-h treatment at 1°C (RCH) compared with control, 268 and 606 unigenes were differentially expressed in L. trifolii and L. sativae, respectively. When comparing pupae exposed to 2h cold shock directly with pupae went through 4h acclimation prior to 2h cold shock, 60 and 399 unigenes were differentially expressed in L trifolii and L sativae, respectively. Genes that were commonly expressed in both L. trifolii and L. sativae, included cytochrome P450, cuticular protein, glucose dehydrogenase, solute carrier family 22 and cationic amino acid transporter. Additionally, several pathways including galactose metabolism and peroxisome were significantly enriched during RCH. Our results show that the transcriptional response is correlated with RCH in the pupal stage of the two Liriomyza species, but more transcriptional changes were identified in L sativae than in L. trifolii.

Differential Transcriptome Analysis Reveals Genes Related to Low- and High-Temperature Stress in the Fall Armyworm, Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is regarded as one of the world’s most harmful plant pests. This research examines the molecular response processes of FAW to low temperature (4°C) and high temperature (40°C) when gene expression is compared to controls (25°C). A total of 211,967 unigenes were collected, at least 14,338 of which were annotated with gene descriptions, gene ontology terms, and metabolic pathways. There were 50 Gene Ontology (GO) functional sub-groups and 21 EggNOG words as a result. Differentially expresses genes (DEGs) with log₂FC ≥ 2 were identified and compared at various temperatures. In comparison to the 25°C treated group, we discovered 199 and 1,248 individual DEGs co-regulated at 4 and 40°C, respectively. Comparing transcriptome profiles for differential gene expression revealed a number of DEGs, including cytochrome P450, odorant binding proteins (OBPs), and immune system genes previously implicated in cold and high temperature stresses. The enrichment pathways were identified using Kyoto Encyclopedia of Genes and Genomics (KEGG) analysis, and heatmaps of similar unigenes from both treatment groups (T4 and T40) were plotted. We used quantitative reverse transcription PCR (RT-qPCR) to confirm the RNA-seq data on 10 up- and down-regulated DEGs. These findings provide a foundation for future understanding of FAW adaptation mechanisms and the underlying basis underlying the response to low and high temperatures.

Genome-wide gene expression profiles of the pea aphid (Acyrthosiphon pisum) under cold temperatures provide insights into body color variation

Cold temperatures are one of the factors influencing color polymorphisms in Acyrthosiphon pisum, resulting in a change from a red to greenish color. Here we characterized gene expression profiles of A. pisum under different low temperatures (1°C, 4°C, 8°C, and 14°C) and durations (3, 6, 12, and 24 h). The number of differentially expressed genes (DEGs) increased as temperatures decreased and time increased, but only a small number of significant DEGs were identified. Genes involved in pigment metabolism were downregulated. An interaction network analysis for 506 common DEGs in comparisons among aphids exposed to 1°C for four durations indicated that a cytochrome P450 gene (CYP, LOC112935894) significantly downregulated may interact with a carotenoid metabolism gene (LOC100574964), similar to other genes encoding CYP, lycopene dehydrogenase and fatty acid synthase. We proposed that the body color shift in A. pisum responding to low temperatures may be regulated by CYPs.

Comparison of gene expression in the red imported fire ant (Solenopsis invicta) under different temperature conditions

The red imported fire ant (RIFA), Solenopsis invicta Buren is native to South America and is known as a global problematic invasive species. This study focused on the molecular response of RIFA by comparing gene expression profiles after exposing ants to low (10 °C) and high (40 °C) temperature stress and comparing them to untreated controls (30 °C). A total of 99,085 unigenes (the clustered non-redundant transcripts that are filtered from the longest assembled contigs) were obtained, of which 19,154 were annotated with gene descriptions, gene ontology terms, and metabolic pathways. 86 gene ontology (GO) functional sub-groups and 23 EggNOG terms resulted. Differentially expressed genes (DEGs) with log₂FC ≥ 10 were screened and were compared at different temperatures. We found 203, 48, and 66 specific DEGs co-regulated at 10, 20, and 40 °C. Comparing transcriptome profiles for differential gene expression resulted in various DE genes, including cytochrome P450, NADH dehydrogenase subunit 1, cuticle protein and heat shock protein (HSP), which have previously been reported to be involved in cold and high temperature resistance. GO analysis revealed that antioxidant activity is up-regulated under high temperature stress. We verified the RNA-seq data by qPCR on 20 up- and down-regulated DEGs. These findings provide a basis for future understanding of the adaptation mechanisms of RIFA and the molecular mechanisms underlying the response to low and high temperatures.

Expression analysis of genes related to cold tolerance in Dendroctonus valens

Pine beetles are well known in North America for their widespread devastation of pine forests. However, Dendroctonus valens LeConte is an important invasive forest pest in China also. Adults and larvae of this bark beetle mainly winter at the trunks and roots of Pinus tabuliformis and Pinus sylvestris; larvae, in particular, result in pine weakness or even death. Since the species was introduced from the United States to Shanxi in 1998, its distribution has spread northward. In 2017, it invaded a large area at the junction of Liaoning, Inner Mongolia and Hebei provinces, showing strong cold tolerance. To identify genes relevant to cold tolerance and the process of overwintering, we sequenced the transcriptomes of wintering and non-wintering adult and larval D. valens using the Illumina HiSeq platform. Differential expression analysis methods for other non-model organisms were used to compare transcript abundances in adults and larvae at two time periods, followed by the identification of functions and metabolic pathways related to genes associated with cold tolerance. We detected 4,387 and 6,091 differentially expressed genes (DEGs) between sampling dates in larvae and adults, respectively, and 1,140 common DEGs, including genes encoding protein phosphatase, very long-chain fatty acids protein, cytochrome P450, and putative leucine-rich repeat-containing proteins. In a Gene Ontology (GO) enrichment analysis, 1,140 genes were assigned to 44 terms, with significant enrichment for cellulase activity, hydrolase activity, and carbohydrate metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) classification and enrichment analyses showed that the lysosomal and purine metabolism pathways involved the most DEGs, the highly enriched terms included autophagy-animal, pentose and glucuronate interconversions and lysosomal processes. We identified 140 candidate genes associated with cold tolerance, including genes with established roles in this trait (e.g., genes encoding trehalose transporter, fructose-1,6-bisphosphatase, and trehalase). Our comparative transcriptome analysis of adult and larval D. valens in different conditions provides basic data for the discovery of key genes and molecular mechanisms underlying cold tolerance.

Multi-level analysis of reproduction in an Antarctic midge identifies female and male accessory gland products that are altered by larval stress and impact progeny viability

The Antarctic midge, Belgica antarctica, is a wingless, non-biting midge endemic to Antarctica. Larval development requires at least 2 years, but adults live only 2 weeks. The nonfeeding adults mate in swarms and females die shortly after oviposition. Eggs are suspended in a gel of unknown composition that is expressed from the female accessory gland. This project characterizes molecular mechanisms underlying reproduction in this midge by examining differential gene expression in whole males, females, and larvae, as well as in male and female accessory glands. Functional studies were used to assess the role of the gel encasing the eggs, as well as the impact of stress on reproductive biology. RNA-seq analyses revealed sex- and development-specific gene sets along with those associated with the accessory glands. Proteomic analyses were used to define the composition of the egg-containing gel, which is generated during multiple developmental stages and derived from both the accessory gland and other female organs. Functional studies indicate the gel provides a larval food source as well as a buffer for thermal and dehydration stress. All of these function are critical to juvenile survival. Larval dehydration stress directly reduces production of storage proteins and key accessory gland components, a feature that impacts adult reproductive success. Modeling reveals that bouts of dehydration may have a significant impact on population growth. This work lays a foundation for further examination of reproduction in midges and provides new information related to general reproduction in dipterans. A key aspect of this work is that reproduction and stress dynamics, currently understudied in polar organisms, are likely to prove critical in determining how climate change will alter their survivability.

Expansion of LINEs and species-specific DNA repeats drives genome expansion in Asian Gypsy Moths

Two subspecies of Asian gypsy moth (AGM), Lymantria dispar asiatica and L. dispar japonica, pose a serious alien invasive threat to North American forests. Despite decades of research on the ecology and biology of this pest, limited AGM-specific genomic resources are currently available. Here, we report on the genome sequences and functional content of these AGM subspecies. The genomes of L.d. asiatica and L.d. japonica are the largest lepidopteran genomes sequenced to date, totaling 921 and 999 megabases, respectively. Large genome size in these subspecies is driven by the accumulation of specific classes of repeats. Genome-wide metabolic pathway reconstructions suggest strong genomic signatures of energy-related pathways in both subspecies, dominated by metabolic functions related to thermogenesis. The genome sequences reported here will provide tools for probing the molecular mechanisms underlying phenotypic traits that are thought to enhance AGM invasiveness.

Comparative Analysis of Transcriptome Responses to Cold Stress in Galeruca daurica (Coleoptera: Chrysomelidae)

Galeruca daurica (Joannis) has become a new insect pest in the Inner Mongolia grasslands since 2009, and its larvae and eggs have strong cold tolerance. To get a deeper insight into its molecular mechanisms of cold stress responses, we performed de novo transcriptome assembly for G. daurica by RNA-Seq and compared the transcriptomes of its larvae exposed to five different temperature treatments (-10, -5, 0, 5, and 25°C for 1 h and then recovered at 25°C for 1 h), respectively. Compared with the control (25°C), the numbers of differentially expressed genes (DEGs) decreased from 1,821 to 882, with the temperature declining from 5 to -10°C. Moreover, we obtained 323 coregulated DEGs under different low temperatures. Under four low temperatures (-10, -5, 0, and 5°C), a large number of genes were commonly upregulated during recovery from cold stresses, including those related to cuticle protein, followed by cytochrome P450, clock protein, fatty acid synthase, and fatty acyl-CoA reductase; meanwhile, lots of genes encoding cuticle protein, RNA replication protein, RNA-directed DNA polymerase, and glucose dehydrogenase were commonly downregulated. Our findings provide important clues for further investigations of key genes and molecular mechanisms involved in the adaptation of G. daurica to harsh environments.

New Zealand Tree and Giant Wētā (Orthoptera) Transcriptomics Reveal Divergent Selection Patterns in Metabolic Loci

Exposure to low temperatures requires an organism to overcome physiological challenges. New Zealand wētā belonging to the genera Hemideina and Deinacrida are found across a wide range of thermal environments and therefore subject to varying selective pressures. Here we assess the selection pressures across the wētā phylogeny, with a particular emphasis on identifying genes under positive or diversifying selection. We used RNA-seq to generate transcriptomes for all 18 Deinacrida and Hemideina species. A total of 755 orthologous genes were identified using a bidirectional best-hit approach, with the resulting gene set encompassing a diverse range of functional classes. Analysis of ortholog ratios of synonymous to nonsynonymous amino acid changes found 83 genes that are under positive selection for at least one codon. A wide variety of Gene Ontology terms, enzymes, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways are represented among these genes. In particular, enzymes involved in oxidative phosphorylation, melanin synthesis, and free-radical scavenging are represented, consistent with physiological and metabolic changes that are associated with adaptation to alpine environments. Structural alignment of the transcripts with the most codons under positive selection revealed that the majority of sites are surface residues, and therefore have the potential to influence the thermostability of the enzyme, with the exception of prophenoloxidase where two residues near the active site are under selection. These proteins provide interesting candidates for further analysis of protein evolution.

Comparative transcriptome analysis of Glyphodes pyloalis Walker (Lepidoptera: Pyralidae) reveals novel insights into heat stress tolerance in insects

BACKGROUND

Heat tolerance is a key parameter that affects insect distribution and abundance. Glyphodes pyloalis Walker (Lepidoptera: Pyralidae) is a devastating pest of mulberry in the main mulberry-growing regions and can cause tremendous losses to sericulture by directly feeding on mulberry leaves and transmitting viruses to Bombyx mori. Moreover, G. pyloalis shows a prominent capacity for adaptation to daily and seasonal temperature fluctuations and can survive several hours under high temperature. To date, the molecular mechanism underlying the outstanding adaptability of this pest to high temperature remains unclear.

RESULTS

In this study, we performed comparative transcriptome analyses on G. pyloalis exposed to 25 and 40 °C for 4 h. We obtained 34,034 unigenes and identified 1275 and 1222 genes significantly upregulated or downregulated, respectively, by heat stress. Data from the transcriptome analyses indicated that some processes involved in heat tolerance are conserved, such as high expression of heat shock protein (HSP) genes and partial repression of metabolism progress. In addition, vitamin digestion and absorption pathways and detoxification pathways identified here provided new insights for the investigation of the molecular mechanisms of heat stress tolerance. Furthermore, transcriptome analysis indicated that immune and phosphatidylinositol signaling system have a close relationship with heat tolerance. In addition, the expression patterns of ten randomly selected genes, such as HSP and cytochrome P450, were consistent with the transcriptome results obtained through quantitative real-time PCR.

CONCLUSIONS

Comparisons among transcriptome results revealed the upregulation of HSPs and genes involved in redox homeostasis, detoxication, and immune progress. However, many metabolism progresses, such as glycolysis/gluconeogenesis and fatty acid biosynthesis, were partially repressed. The results reflected that the heat tolerance of G. pyloalis is a fairly complicated process and related to a broad range of physiological regulations. Our study can improve understanding on the mechanisms of insect thermal tolerance.

Differential transcriptome analysis reveals genes related to cold tolerance in seabuckthorn carpenter moth, Eogystia hippophaecolus

Seabuckthorn carpenter moth, Eogystia hippophaecolus (Lepidoptera: Cossidae), is an important pest of sea buckthorn (Hippophae rhamnoides), which is a shrub that has significant ecological and economic value in China. E. hippophaecolus is highly cold tolerant, but limited studies have been conducted to elucidate the molecular mechanisms underlying its cold resistance. Here we sequenced the E. hippophaecolus transcriptome using RNA-Seq technology and performed de novo assembly from the short paired-end reads. We investigated the larval response to cold stress by comparing gene expression profiles between treatments. We obtained 118,034 unigenes, of which 22,161 were annotated with gene descriptions, conserved domains, gene ontology terms, and metabolic pathways. These resulted in 57 GO terms and 193 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. By comparing transcriptome profiles for differential gene expression, we identified many differentially expressed proteins and genes, including heat shock proteins and cuticular proteins which have previously been reported to be involved in cold resistance of insects. This study provides a global transcriptome analysis and an assessment of differential gene expression in E. hippophaecolus under cold stress. We found seven differential expressed genes in common between developmental stages, which were verified with qPCR. Our findings facilitate future genomic studies aimed at improving our understanding of the molecular mechanisms underlying the response of insects to low temperatures.

Cold acclimation wholly reorganizes the Drosophila melanogaster transcriptome and metabolome

Cold tolerance is a key determinant of insect distribution and abundance, and thermal acclimation can strongly influence organismal stress tolerance phenotypes, particularly in small ectotherms like Drosophila. However, there is limited understanding of the molecular and biochemical mechanisms that confer such impressive plasticity. Here, we use high-throughput mRNA sequencing (RNA-seq) and liquid chromatography – mass spectrometry (LC-MS) to compare the transcriptomes and metabolomes of D. melanogaster acclimated as adults to warm (rearing) (21.5 °C) or cold conditions (6 °C). Cold acclimation improved cold tolerance and led to extensive biological reorganization: almost one third of the transcriptome and nearly half of the metabolome were differentially regulated. There was overlap in the metabolic pathways identified via transcriptomics and metabolomics, with proline and glutathione metabolism being the most strongly-supported metabolic pathways associated with increased cold tolerance. We discuss several new targets in the study of insect cold tolerance (e.g. dopamine signaling and Na⁺-driven transport), but many previously identified candidate genes and pathways (e.g. heat shock proteins, Ca²⁺ signaling, and ROS detoxification) were also identified in the present study, and our results are thus consistent with and extend the current understanding of the mechanisms of insect chilling tolerance.

De Novo Transcriptome Analysis of the Common New Zealand Stick Insect Clitarchus hookeri (Phasmatodea) Reveals Genes Involved in Olfaction, Digestion and Sexual Reproduction

Phasmatodea, more commonly known as stick insects, have been poorly studied at the molecular level for several key traits, such as components of the sensory system and regulators of reproduction and development, impeding a deeper understanding of their functional biology. Here, we employ de novo transcriptome analysis to identify genes with primary functions related to female odour reception, digestion, and male sexual traits in the New Zealand common stick insect Clitarchus hookeri (White). The female olfactory gene repertoire revealed ten odorant binding proteins with three recently duplicated, 12 chemosensory proteins, 16 odorant receptors, and 17 ionotropic receptors. The majority of these olfactory genes were over-expressed in female antennae and have the inferred function of odorant reception. Others that were predominantly expressed in male terminalia (n = 3) and female midgut (n = 1) suggest they have a role in sexual reproduction and digestion, respectively. Over-represented transcripts in the midgut were enriched with digestive enzyme gene families. Clitarchus hookeri is likely to harbour nine members of an endogenous cellulase family (glycoside hydrolase family 9), two of which appear to be specific to the C. hookeri lineage. All of these cellulase sequences fall into four main phasmid clades and show gene duplication events occurred early in the diversification of Phasmatodea. In addition, C. hookeri genome is likely to express γ-proteobacteria pectinase transcripts that have recently been shown to be the result of horizontal transfer. We also predicted 711 male terminalia-enriched transcripts that are candidate accessory gland proteins, 28 of which were annotated to have molecular functions of peptidase activity and peptidase inhibitor activity, two groups being widely reported to regulate female reproduction through proteolytic cascades. Our study has yielded new insights into the genetic basis of odour detection, nutrient digestion, and male sexual traits in stick insects. The C. hookeri reference transcriptome, together with identified gene families, provides a comprehensive resource for studying the evolution of sensory perception, digestive systems, and reproductive success in phasmids.

Sticky Genomes: Using NGS Evidence to Test Hybrid Speciation Hypotheses

Hypotheses of hybrid origin are common. Here we use next generation sequencing to test a hybrid hypothesis for a non-model insect with a large genome. We compared a putative hybrid triploid stick insect species (Acanthoxyla geisovii) with its putative paternal diploid taxon (Clitarchus hookeri), a relationship that provides clear predictions for the relative genetic diversity within each genome. The parental taxon is expected to have comparatively low allelic diversity that is nested within the diversity of the hybrid daughter genome. The scale of genome sequencing required was conveniently achieved by extracting mRNA and sequencing cDNA to examine expressed allelic diversity. This allowed us to test hybrid-progenitor relationships among non-model organisms with large genomes and different ploidy levels. Examination of thousands of independent loci avoids potential problems produced by the silencing of parts of one or other of the parental genomes, a phenomenon sometimes associated with the process of stabilisation of a hybrid genome. Transcript assembles were assessed for evidence of paralogs and/or alternative splice variants before proceeding. Comparison of transcript assemblies was not an appropriate measure of genetic variability, but by mapping reads back to clusters derived from each species we determined levels of allelic diversity. We found greater cDNA sequence diversity among alleles in the putative hybrid species (Acanthoxyla geisovii) than the non-hybrid. The allelic diversity within the putative paternal species (Clitachus hookeri) nested within the hybrid-daughter genome, supports the current view of a hybrid-progenitor relationship for these stick insect species. Next generation sequencing technology provides opportunities for testing evolutionary hypotheses with non-model organisms, including, as here, genomes that are large due to polyploidy.

PROTEOMICS ANALYSIS OF OVEREXPRESSED PLASMA PROTEINS IN RESPONSE TO COLD ACCLIMATION IN Ostrinia furnacalis

Many insects in temperate regions overwinter in diapause. In these insects, one of the metabolic adaptations to cold stress is the synthesis of responsive proteins. Using proteomic analysis, an investigation aimed to a better understanding of the molecular adaptation mechanisms to cold stress was carried out in Ostrinia furnacalis larva. Proteins were extracted from the larval hemolymph collected from both control and overwintering larva. By polyethylene glycol precipitation, approximately 560 protein spots were separated and visualized on two-dimensional (2D) gels after silver staining. Eighteen protein spots were found to be upregulated in overwinter larval plasma in different patterns. As an initial work, 13 of these proteins were identified using MALDI TOF/TOF MS. The differentially overexpressed proteins include heat shock 70 kDa cognate protein, small heat shock protein (sHSP), putative aliphatic nitrilase, arginine kinase, phosphoglyceromutase, triosephosphateisomerase, and glutathione transferase. Alterations in the levels of these proteins were further confirmed by qPCR. This study is the first analysis of differentially expressed plasma proteins in O. furnacalis diapause larvae under extremely low temperature conditions and gives new insights into the acclimation mechanisms responsive to cold stress. Our results also support the idea that energy metabolism, alanine and proline metabolism, and antioxidative reaction act in the cold acclimation of O. furnacalis diapause larvae.

Parallel molecular routes to cold adaptation in eight genera of New Zealand stick insects

The acquisition of physiological strategies to tolerate novel thermal conditions allows organisms to exploit new environments. As a result, thermal tolerance is a key determinant of the global distribution of biodiversity, yet the constraints on its evolution are not well understood. Here we investigate parallel evolution of cold tolerance in New Zealand stick insects, an endemic radiation containing three montane-occurring species. Using a phylogeny constructed from 274 orthologous genes, we show that stick insects have independently colonized montane environments at least twice. We compare supercooling point and survival of internal ice formation among ten species from eight genera, and identify both freeze tolerance and freeze avoidance in separate montane lineages. Freeze tolerance is also verified in both lowland and montane populations of a single, geographically widespread, species. Transcriptome sequencing following cold shock identifies a set of structural cuticular genes that are both differentially regulated and under positive sequence selection in each species. However, while cuticular proteins in general are associated with cold shock across the phylogeny, the specific genes at play differ among species. Thus, while processes related to cuticular structure are consistently associated with adaptation for cold, this may not be the consequence of shared ancestral genetic constraints.

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.

The accuracy, feasibility and challenges of sequencing short tandem repeats using next-generation sequencing platforms

To date we have little knowledge of how accurate next-generation sequencing (NGS) technologies are in sequencing repetitive sequences beyond known limitations to accurately sequence homopolymers. Only a handful of previous reports have evaluated the potential of NGS for sequencing short tandem repeats (microsatellites) and no empirical study has compared and evaluated the performance of more than one NGS platform with the same dataset. Here we examined yeast microsatellite variants from both long-read (454-sequencing) and short-read (Illumina) NGS platforms and compared these to data derived through Sanger sequencing. In addition, we investigated any locus-specific biases and differences that might have resulted from variability in microsatellite repeat number, repeat motif or type of mutation. Out of 112 insertion/deletion variants identified among 45 microsatellite amplicons in our study, we found 87.5% agreement between the 454-platform and Sanger sequencing in frequency of variant detection after Benjamini-Hochberg correction for multiple tests. For a subset of 21 microsatellite amplicons derived from Illumina sequencing, the results of short-read platform were highly consistent with the other two platforms, with 100% agreement with 454-sequencing and 93.6% agreement with the Sanger method after Benjamini-Hochberg correction. We found that the microsatellite attributes copy number, repeat motif and type of mutation did not have a significant effect on differences seen between the sequencing platforms. We show that both long-read and short-read NGS platforms can be used to sequence short tandem repeats accurately, which makes it feasible to consider the use of these platforms in high-throughput genotyping. It appears the major requirement for achieving both high accuracy and rare variant detection in microsatellite genotyping is sufficient read depth coverage. This might be a challenge because each platform generates a consistent pattern of non-uniform sequence coverage, which, as our study suggests, may affect some types of tandem repeats more than others.

Cloning and expression analysis of four heat shock protein genes in Ericerus pela (Homoptera: Coccidae)

To explore the function of small heat shock protein genes (shsps) and hsp70 in Ericerus pela, we cloned the full-length cDNA sequences of hsp21.5, hsp21.7, hsp70, and hsc70 and the genomic sequence of hsc70. Open reading frames of the four hsps were 570, 564, 1,908, and 1,962 base pairs (bp), respectively, which encode proteins with calculated molecular mass of 21.5, 21.7, 69.8, and 71.6 kDa. Amino acid sequence analysis revealed the presence of the conserved Hsp motifs in all four proteins. The genomic DNA of hsc70 had four introns. ep-hsp21.5 was orthologous and ep-hsp21.7 was species specific. Expression of all four transcripts during heat or cold stress and development was examined by quantitative real-time polymerase chain reaction. All four hsps were upregulated during heat or cold stress in female adults, indicating a correlation between the four hsps and heat or cold-stress tolerance in female adults. ep-hsp21.7 and ep-hsp70 were upregulated during heat stress in male larvae, implying a correlation between the two hsps and heat-stress tolerance in male larvae. The four ep-hsps were also upregulated during the developmental process in males, and ep-hsp21.5, ep-hsp70, and ep-hsc70 were upregulated in females, which indicates their possible role in the developmental regulation of E. pela.

De novo transcriptome of the desert beetle Microdera punctipennis (Coleoptera: Tenebrionidae) using Illumina RNA-seq technology

Insects in Tenebrionidae have unique stress adaptations that allow them to survive temperature extremes. We report here a gene expression profiling of Microdera punctipennis, a beetle in desert region, to gain a global view of its environmental adaptations. A total of 48,158,004 reads were obtained by transcriptome sequencing, and the de novo assembly yielded 56,348 unigenes with an average length of 666 bp. Based on similarity searches with a cut-off E-value of 10⁻⁵ against two protein sequence databases, 41,109 of the unigenes (about 72.96 %) were matched to known proteins. An in-depth analysis of the data revealed a large number of genes were associated with environmental stress, including genes that encode heat shock proteins, antifreeze proteins, and enzymes such as chitinase, trehalose, and trehalose-6-phosphate synthase. This study generated a substantial number of M. punctipennis transcript sequences that can be used to discover novel genes associated with stress adaptation. These sequences are a valuable resource for future studies of the desert beetle and other insects in Tenebrionidae. Transcriptome analysis based on Illumina paired-end sequencing is a powerful approach for gene discovery and molecular marker development for non-model species.

Positive selection in glycolysis among Australasian stick insects

BACKGROUND

The glycolytic pathway is central to cellular energy production. Selection on individual enzymes within glycolysis, particularly phosphoglucose isomerase (Pgi), has been associated with metabolic performance in numerous organisms. Nonetheless, how whole energy-producing pathways evolve to allow organisms to thrive in different environments and adopt new lifestyles remains little explored. The Lanceocercata radiation of Australasian stick insects includes transitions from tropical to temperate climates, lowland to alpine habitats, and winged to wingless forms. This permits a broad investigation to determine which steps within glycolysis and what sites within enzymes are the targets of positive selection. To address these questions we obtained transcript sequences from seven core glycolysis enzymes, including two Pgi paralogues, from 29 Lanceocercata species.

RESULTS

Using maximum likelihood methods a signature of positive selection was inferred in two core glycolysis enzymes. Pgi and Glyceraldehyde 3-phosphate dehydrogenase (Gaphd) genes both encode enzymes linking glycolysis to the pentose phosphate pathway. Positive selection among Pgi paralogues and orthologues predominately targets amino acids with residues exposed to the protein's surface, where changes in physical properties may alter enzyme performance.

CONCLUSION

Our results suggest that, for Lancerocercata stick insects, adaptation to new stressful lifestyles requires a balance between maintaining cellular energy production, efficiently exploiting different energy storage pools and compensating for stress-induced oxidative damage.

The importance of tissue specificity for RNA-seq: highlighting the errors of composite structure extractions

Background

A composite biological structure, such as an insect head or abdomen, contains many internal structures with distinct functions. Composite structures are often used in RNA-seq studies, though it is unclear how expression of the same gene in different tissues and structures within the same structure affects the measurement (or even utility) of the resulting patterns of gene expression. Here we determine how complex composite tissue structure affects measures of gene expression using RNA-seq.

Results

We focus on two structures in the honey bee (the sting gland and digestive tract) both contained within one larger structure, the whole abdomen. For each of the three structures, we used RNA-seq to identify differentially expressed genes between two developmental stages, nurse bees and foragers. Based on RNA-seq for each structure-specific extraction, we found that RNA-seq with composite structures leads to many false negatives (genes strongly differentially expressed in particular structures which are not found to be differentially expressed within the composite structure). We also found a significant number of genes with one pattern of differential expression in the tissue-specific extraction, and the opposite in the composite extraction, suggesting multiple signals from such genes within the composite structure. We found these patterns for different classes of genes including transcription factors.

Conclusions

Many RNA-seq studies currently use composite extractions, and even whole insect extractions, when tissue and structure specific extractions are possible. This is due to the logistical difficultly of micro-dissection and unawareness of the potential errors associated with composite extractions. The present study suggests that RNA-seq studies of composite structures are prone to false negatives and difficult to interpret positive signals for genes with variable patterns of local expression. In general, our results suggest that RNA-seq on large composite structures should be avoided unless it is possible to demonstrate that the effects shown here do not exist for the genes of interest.