Gene and Allele-Specific Expression Underlying the Electric Signal Divergence in African Weakly Electric Fish.

In the African weakly electric fish genus Campylomormyrus, electric organ discharge signals are strikingly different in shape and duration among closely related species, contribute to prezygotic isolation, and may have triggered an adaptive radiation. We performed mRNA sequencing on electric organs and skeletal muscles (from which the electric organs derive) from 3 species with short (0.4 ms), medium (5 ms), and long (40 ms) electric organ discharges and 2 different cross-species hybrids. We identified 1,444 upregulated genes in electric organ shared by all 5 species/hybrid cohorts, rendering them candidate genes for electric organ-specific properties in Campylomormyrus. We further identified several candidate genes, including KCNJ2 and KLF5, and their upregulation may contribute to increased electric organ discharge duration. Hybrids between a short (Campylomormyrus compressirostris) and a long (Campylomormyrus rhynchophorus) discharging species exhibit electric organ discharges of intermediate duration and showed imbalanced expression of KCNJ2 alleles, pointing toward a cis-regulatory difference at this locus, relative to electric organ discharge duration. KLF5 is a transcription factor potentially balancing potassium channel gene expression, a crucial process for the formation of an electric organ discharge. Unraveling the genetic basis of the species-specific modulation of the electric organ discharge in Campylomormyrus is crucial for understanding the adaptive radiation of this emerging model taxon of ecological (perhaps even sympatric) speciation.

A new genome assembly of an African weakly electric fish (Campylomormyrus compressirostris, Mormyridae) indicates rapid gene family evolution in Osteoglossomorpha.

BACKGROUND: Teleost fishes comprise more than half of the vertebrate species. Within teleosts, most phylogenies consider the split between Osteoglossomorpha and Euteleosteomorpha/Otomorpha as basal, preceded only by the derivation of the most primitive group of teleosts, the Elopomorpha. While Osteoglossomorpha are generally species poor, the taxon contains the African weakly electric fish (Mormyroidei), which have radiated into numerous species. Within the mormyrids, the genus Campylomormyrus is mostly endemic to the Congo Basin. Campylomormyrus serves as a model to understand mechanisms of adaptive radiation and ecological speciation, especially with regard to its highly diverse species-specific electric organ discharges (EOD). Currently, there are few well-annotated genomes available for electric fish in general and mormyrids in particular. Our study aims at producing a high-quality genome assembly and to use this to examine genome evolution in relation to other teleosts. This will facilitate further understanding of the evolution of the osteoglossomorpha fish in general and of electric fish in particular. RESULTS: A high-quality weakly electric fish (C. compressirostris) genome was produced from a single individual with a genome size of 862 Mb, consisting of 1,497 contigs with an N50 of 1,399 kb and a GC-content of 43.69%. Gene predictions identified 34,492 protein-coding genes, which is a higher number than in the two other available Osteoglossomorpha genomes of Paramormyrops kingsleyae and Scleropages formosus. A Computational Analysis of gene Family Evolution (CAFE5) comparing 33 teleost fish genomes suggests an overall faster gene family turnover rate in Osteoglossomorpha than in Otomorpha and Euteleosteomorpha. Moreover, the ratios of expanded/contracted gene family numbers in Osteoglossomorpha are significantly higher than in the other two taxa, except for species that had undergone an additional genome duplication (Cyprinus carpio and Oncorhynchus mykiss). As potassium channel proteins are hypothesized to play a key role in EOD diversity among species, we put a special focus on them, and manually curated 16 Kv1 genes. We identified a tandem duplication in the KCNA7a gene in the genome of C. compressirostris. CONCLUSIONS: We present the fourth genome of an electric fish and the third well-annotated genome for Osteoglossomorpha, enabling us to compare gene family evolution among major teleost lineages. Osteoglossomorpha appear to exhibit rapid gene family evolution, with more gene family expansions than contractions. The curated Kv1 gene family showed seven gene clusters, which is more than in other analyzed fish genomes outside Osteoglossomorpha. The KCNA7a, encoding for a potassium channel central for EOD production and modulation, is tandemly duplicated which may related to the diverse EOD observed among Campylomormyrus species.

Intergenerational Effects of Early-Life Starvation on Life History, Consumption, and Transcriptome of a Holometabolous Insect.

AbstractIntergenerational effects, also known as parental effects in which the offspring phenotype is influenced by the parental phenotype, can occur in response to factors that occur not only in early but also in late parental life. However, little is known about how these parental life stage-specific environments interact with each other and with the offspring environment to influence offspring phenotypes, particularly in organisms that realize distinct niches across ontogeny. We examined the effects of parental larval starvation and adult reproductive environment on offspring traits under matching or mismatching offspring larval starvation conditions using the holometabolous, haplodiploid insect Athalia rosae (turnip sawfly). We show that parental larval starvation had trait-dependent intergenerational effects on both life history and consumption traits of offspring larvae, partly in interaction with offspring conditions, while there was no significant effect of parental adult reproductive environment. In addition, while offspring larval starvation led to numerous gene- and pathway-level expression differences, parental larval starvation impacted far fewer genes and only the ribosomal pathway. Our findings reveal that parental starvation evokes complex intergenerational effects on offspring life history traits, consumption patterns, and gene expression, although the effects are less pronounced than those of offspring starvation.

Triple RNA-Seq characterizes aphid gene expression in response to infection with unequally virulent strains of the endosymbiont Hamiltonella defensa.

BACKGROUND: Secondary endosymbionts of aphids provide benefits to their hosts, but also impose costs such as reduced lifespan and reproductive output. The aphid Aphis fabae is host to different strains of the secondary endosymbiont Hamiltonella defensa, which encode different putative toxins. These strains have very different phenotypes: They reach different densities in the host, and the costs and benefits (protection against parasitoid wasps) they confer to the host vary strongly. RESULTS: We used RNA-Seq to generate hypotheses on why four of these strains inflict such different costs to A. fabae. We found different H. defensa strains to cause strain-specific changes in aphid gene expression, but little effect of H. defensa on gene expression of the primary endosymbiont, Buchnera aphidicola. The highly costly and over-replicating H. defensa strain H85 was associated with strongly reduced aphid expression of hemocytin, a marker of hemocytes in Drosophila. The closely related strain H15 was associated with downregulation of ubiquitin-related modifier 1, which is related to nutrient-sensing and oxidative stress in other organisms. Strain H402 was associated with strong differential regulation of a set of hypothetical proteins, the majority of which were only differentially regulated in presence of H402. CONCLUSIONS: Overall, our results suggest that costs of different strains of H. defensa are likely caused by different mechanisms, and that these costs are imposed by interacting with the host rather than the host's obligatory endosymbiont B. aphidicola.

Temperature-dependent life history and transcriptomic responses in heat-tolerant versus heat-sensitive Brachionus rotifers.

Thermal stress response is an essential physiological trait that determines occurrence and temporal succession in nature, including response to climate change. We compared temperature-related demography in closely related heat-tolerant and heat-sensitive Brachionus rotifer species. We found significant differences in heat response, with the heat-sensitive species adopting a strategy of long survival and low population growth, while the heat-tolerant followed the opposite strategy. In both species, we examined the genetic basis of physiological variation by comparing gene expression across increasing temperatures. Comparative transcriptomic analyses identified shared and opposing responses to heat. Interestingly, expression of heat shock proteins (hsps) was strikingly different in the two species and mirrored differences in population growth rates, showing that hsp genes are likely a key component of a species' adaptation to different temperatures. Temperature induction caused opposing patterns of expression in further functional categories including energy, carbohydrate and lipid metabolism, and in genes related to ribosomal proteins. In the heat-sensitive species, elevated temperatures caused up-regulation of genes related to meiosis induction and post-translational histone modifications. This work demonstrates the sweeping reorganizations of biological functions that accompany temperature adaptation in these two species and reveals potential molecular mechanisms that might be activated for adaptation to global warming.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum.

BACKGROUND: Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts. RESULTS: We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes. CONCLUSIONS: These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

A new lineage of segmented RNA viruses infecting animals.

Metagenomic sequencing has revolutionised our knowledge of virus diversity, with new virus sequences being reported faster than ever before. However, virus discovery from metagenomic sequencing usually depends on detectable homology: without a sufficiently close relative, so-called 'dark' virus sequences remain unrecognisable. An alternative approach is to use virus-identification methods that do not depend on detecting homology, such as virus recognition by host antiviral immunity. For example, virus-derived small RNAs have previously been used to propose 'dark' virus sequences associated with the Drosophilidae (Diptera). Here, we combine published Drosophila data with a comprehensive search of transcriptomic sequences and selected meta-transcriptomic datasets to identify a completely new lineage of segmented positive-sense single-stranded RNA viruses that we provisionally refer to as the Quenyaviruses. Each of the five segments contains a single open reading frame, with most encoding proteins showing no detectable similarity to characterised viruses, and one sharing a small number of residues with the RNA-dependent RNA polymerases of single- and double-stranded RNA viruses. Using these sequences, we identify close relatives in approximately 20 arthropods, including insects, crustaceans, spiders, and a myriapod. Using a more conserved sequence from the putative polymerase, we further identify relatives in meta-transcriptomic datasets from gut, gill, and lung tissues of vertebrates, reflecting infections of vertebrates or of their associated parasites. Our data illustrate the utility of small RNAs to detect viruses with limited sequence conservation, and provide robust evidence for a new deeply divergent and phylogenetically distinct RNA virus lineage.

A Novel RNA Virus in the Parasitoid Wasp Lysiphlebus fabarum: Genomic Structure, Prevalence, and Transmission.

We report on a novel RNA virus infecting the wasp Lysiphlebus fabarum, a parasitoid of aphids. This virus, tentatively named "Lysiphlebus fabarum virus" (LysV), was discovered in transcriptome sequences of wasps from an experimental evolution study in which the parasitoids were allowed to adapt to aphid hosts (Aphis fabae) with or without resistance-conferring endosymbionts. Based on phylogenetic analyses of the viral RNA-dependent RNA polymerase (RdRp), LysV belongs to the Iflaviridae family in the order of the Picornavirales, with the closest known relatives all being parasitoid wasp-infecting viruses. We developed an endpoint PCR and a more sensitive qPCR assay to screen for LysV in field samples and laboratory lines. These screens verified the occurrence of LysV in wild parasitoids and identified the likely wild-source population for lab infections in Western Switzerland. Three viral haplotypes could be distinguished in wild populations, of which two were found in the laboratory. Both vertical and horizontal transmission of LysV were demonstrated experimentally, and repeated sampling of laboratory populations suggests that the virus can form persistent infections without obvious symptoms in infected wasps.

Mapping of Multiple Complementary Sex Determination Loci in a Parasitoid Wasp.

Sex determination has evolved in a variety of ways and can depend on environmental and genetic signals. A widespread form of genetic sex determination is haplodiploidy, where unfertilized, haploid eggs develop into males and fertilized diploid eggs into females. One of the molecular mechanisms underlying haplodiploidy in Hymenoptera, the large insect order comprising ants, bees, and wasps, is complementary sex determination (CSD). In species with CSD, heterozygosity at one or several loci induces female development. Here, we identify the genomic regions putatively underlying multilocus CSD in the parasitoid wasp Lysiphlebus fabarum using restriction-site associated DNA sequencing. By analyzing segregation patterns at polymorphic sites among 331 diploid males and females, we identify up to four CSD candidate regions, all on different chromosomes. None of the candidate regions feature evidence for homology with the csd gene from the honey bee, the only species in which CSD has been characterized, suggesting that CSD in L. fabarum is regulated via a novel molecular mechanism. Moreover, no homology is shared between the candidate loci, in contrast to the idea that multilocus CSD should emerge from duplications of an ancestral single-locus system. Taken together, our results suggest that the molecular mechanisms underlying CSD in Hymenoptera are not conserved between species, raising the question as to whether CSD may have evolved multiple times independently in the group.

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto.

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena).

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Rapid evolution of symbiont-mediated resistance compromises biological control of aphids by parasitoids.

There is growing interest in biological control as a sustainable and environmentally friendly way to control pest insects. Aphids are among the most detrimental agricultural pests worldwide, and parasitoid wasps are frequently employed for their control. The use of asexual parasitoids may improve the effectiveness of biological control because only females kill hosts and because asexual populations have a higher growth rate than sexuals. However, asexuals may have a reduced capacity to track evolutionary change in their host populations. We used a factorial experiment to compare the ability of sexual and asexual populations of the parasitoid Lysiphlebus fabarum to control caged populations of black bean aphids (Aphis fabae) of high and low clonal diversity. The aphids came from a natural population, and one-third of the aphid clones harbored Hamiltonella defensa, a heritable bacterial endosymbiont that increases resistance to parasitoids. We followed aphid and parasitoid population dynamics for 3 months but found no evidence that the reproductive mode of parasitoids affected their effectiveness as biocontrol agents, independent of host clonal diversity. Parasitoids failed to control aphids in most cases, because their introduction resulted in strong selection for clones protected by H. defensa. The increasingly resistant aphid populations escaped control by parasitoids, and we even observed parasitoid extinctions in many cages. The rapid evolution of symbiont-conferred resistance in turn imposed selection on parasitoids. In cages where asexual parasitoids persisted until the end of the experiment, they became dominated by a single genotype able to overcome the protection provided by H. defensa. Thus, there was evidence for parasitoid counteradaptation, but it was generally too slow for parasitoids to regain control over aphid populations. It appears that when pest aphids possess defensive symbionts, the presence of parasitoid genotypes able to overcome symbiont-conferred resistance is more important for biocontrol success than their reproductive mode.

Positive selection and comparative molecular evolution of reproductive proteins from New Zealand tree weta (Orthoptera, Hemideina).

Animal reproductive proteins, especially those in the seminal fluid, have been shown to have higher levels of divergence than non-reproductive proteins and are often evolving adaptively. Seminal fluid proteins have been implicated in the formation of reproductive barriers between diverging lineages, and hence represent interesting candidates underlying speciation. RNA-seq was used to generate the first male reproductive transcriptome for the New Zealand tree weta species Hemideina thoracica and H. crassidens. We identified 865 putative reproductive associated proteins across both species, encompassing a diverse range of functional classes. Candidate gene sequencing of nine genes across three Hemideina, and two Deinacrida species suggests that H. thoracica has the highest levels of intraspecific genetic diversity. Non-monophyly was observed in the majority of sequenced genes indicating that either gene flow may be occurring between the species, or that reciprocal monophyly at these loci has yet to be attained. Evidence for positive selection was found for one lectin-related reproductive protein, with an overall omega of 7.65 and one site in particular being under strong positive selection. This candidate gene represents the first step in the identification of proteins underlying the evolutionary basis of weta reproduction and speciation.

Parasitoid gene expression changes after adaptation to symbiont-protected hosts.

Reciprocal selection between aphids, their protective endosymbionts, and the parasitoid wasps that prey upon them offers an opportunity to study the basis of their coevolution. We investigated adaptation to symbiont-conferred defense by rearing the parasitoid wasp Lysiphlebus fabarum on aphids (Aphis fabae) possessing different defensive symbiont strains (Hamiltonella defensa). After ten generations of experimental evolution, wasps showed increased abilities to parasitize aphids possessing the H. defensa strain they evolved with, but not aphids possessing the other strain. We show that the two symbiont strains encode different toxins, potentially creating different targets for counter-adaptation. Phenotypic and behavioral comparisons suggest that neither life-history traits nor oviposition behavior differed among evolved parasitoid lineages. In contrast, comparative transcriptomics of adult female wasps identified a suite of differentially expressed genes among lineages, even when reared in a common, symbiont-free, aphid host. In concurrence with the specificity of each parasitoid lineages' infectivity, most differentially expressed parasitoid transcripts were also lineage-specific. These transcripts are enriched with putative venom toxins and contain highly expressed, potentially defensive viral particles. Together, these results suggest that wild populations of L. fabarum employ a complicated offensive arsenal with sufficient genetic variation for wasps to adapt rapidly and specifically to their hosts' microbial defenses.

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(2+) 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.

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.

Divergent transcriptional responses to low temperature among populations of alpine and lowland species of New Zealand stick insects (Micrarchus).

In widespread and genetically structured populations, temperature variation may lead to among-population differentiation of thermal biology. The New Zealand stick insect genus Micrarchus contains four species that inhabit different thermal environments, two of which are geographically widespread. RNA-Seq and quantitative PCR were used to investigate the transcriptional responses to cold shock among lowland and alpine species to identify cold-responsive transcripts that differ between the species and to determine whether there is intraspecific geographical variation in gene expression. We also used mitochondrial DNA, nuclear 28S ribosomal DNA and transcriptome-wide SNPs to determine phylogeographic structure and the potential for differences in genetic backgrounds to contribute to variation in gene expression. RNA-Seq identified 2160 unigenes that were differentially expressed as a result of low-temperature exposure across three populations from two species (M. hystriculeus and M. nov. sp. 2), with a majority (68% ± 20%) being population specific. This extensive geographical variation is consistent across years and is likely a result of background genetic differences among populations caused by genetic drift and possibly local adaptation. Responses to cold shock shared among alpine M. nov. sp. 2 populations included the enrichment of cuticular structure-associated transcripts, suggesting that cuticle modification may have accompanied colonization of low-temperature alpine environments and the development of a more cold-hardy phenotype.

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.

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

The endemic New Zealand alpine stick insect Micrarchus nov. sp. 2 regularly experiences sub-zero temperatures in the wild. 454-based RNA-Seq was used to generate a de novo transcriptome and differentiate between treatments to investigate the genetic basis of cold tolerance. Non cold-treated individuals were compared to those exposed to 0°C for 1 h followed by a 1 h recovery period at 20°C. We aligned 607,410 Roche 454 reads, generating a transcriptome of 5235 contigs. Differential expression analysis ranked candidate cold responsive genes for qPCR validation by P-value. The top nine up-regulated candidates, together with eight a priori targets identified from previous studies, had their relative expression quantified using qPCR. Three candidate cold responsive genes from the RNA-Seq data were verified as significantly up-regulated, annotated as: prolyl 4-hydroxylase subunit alpha-1 (P4HA1), staphylococcal nuclease domain-containing protein 1 (snd1) and cuticular protein analogous to peritrophins 3-D2 (Cpap3-d2). All three are novel candidate genes, illustrating the varied response to low temperature across insects.