The ecoresponsive genome of Daphnia pulex

Daphnia stressor database: Taking advantage of a decade of Daphnia '-omics' data for gene annotation

Gene expression patterns help to measure and characterize the effect of environmental perturbations at the cellular and organism-level. Complicating interpretation is the presence of uncharacterized or "hypothetical" gene functions for a large percentage of genomes. This is particularly evident in Daphnia genomes, which contains many regions coding for "hypothetical proteins" and are significantly divergent from many of the available arthropod model species, but might be ecologically important. In the present study, we developed a gene expression database, the Daphnia stressor database (http://www.daphnia-stressordb.uni-hamburg.de/dsdbstart.php), built from 90 published studies on Daphnia gene expression. Using a comparative genomics approach, we used the database to annotate D. galeata transcripts. The extensive body of literature available for Daphnia species allowed to associate stressors with gene expression patterns. We believe that our stressor based annotation strategy allows for better understanding and interpretation of the functional role of the understudied hypothetical or uncharacterized Daphnia genes, thereby increasing our understanding of Daphnia's genetic and phenotypic variability.

Purification and partial biochemical characterization of recombinant lactate dehydrogenase 1 (LDH-1) of the white shrimp Litopenaeus vannamei

Lactate dehydrogenase (LDH) is a key enzyme to produce energy during hypoxia by anaerobic glycolysis. In the white shrimp Litopenaeus vannamei, two protein subunits (LDH-1 and LDH-2) were previously identified, deduced from two different transcripts that come from the same LDH gene by processing via mutually exclusive alternative splicing. LDH-1 contains exon five and LDH-2 contains exon six and the two proteins differ only in 15 amino acid residues. Both subunits were independently cloned and overexpressed in E. coli as a fusion protein containing a chitin binding domain. Previously, recombinant LDH-2 was successfully purified and characterized, but LDH-1 was insoluble and aggregated forming inclusion bodies. We report the production of soluble LDH-1 by testing different pHs in the buffers used to lyse the bacterial cells before the purification step and the characterization of the purified protein to show that the cDNA indeed codes for a functional and active protein. The recombinant native protein is a homotetramer of approximately 140 kDa composed by 36 kDa subunits and has higher affinity for pyruvate than for lactate. LDH-1 has an optimum pH of 7.5 and is stable between pH 8.0 and 9.0; pH data analysis showed two pKa values of 6.1 ± 0.15 and 8.8 ± 0.15 suggesting a histidine and asparagine, respectively, involved in the active site. The enzyme optimal temperature was 44 °C and it was stable between 20 and 60 °C. LDH-1 was slightly activated by NaCl, KCl and MgCl₂ and fully inhibited by ZnCl₂.

The Molecular Basis of Freshwater Adaptation in Prawns: Insights from Comparative Transcriptomics of Three Macrobrachium Species

Elucidating the molecular basis of adaptation to different environmental conditions is important because adaptive ability of a species can shape its distribution, influence speciation, and also drive a variety of evolutionary processes. For crustaceans, colonization of freshwater habitats has significantly impacted diversity, but the molecular basis of this process is poorly understood. In the current study, we examined three prawn species from the genus Macrobrachium (M. australiense, M. tolmerum, and M. novaehollandiae) to better understand the molecular basis of freshwater adaptation using a comparative transcriptomics approach. Each of these species naturally inhabit environments with different salinity levels; here, we exposed them to the same experimental salinity conditions (0‰ and 15‰), to compare expression patterns of candidate genes that previously have been shown to influence phenotypic traits associated with freshwater adaptation (e.g., genes associated with osmoregulation). Differential gene expression analysis revealed 876, 861, and 925 differentially expressed transcripts under the two salinities for M. australiense, M. tolmerum, and M. novaehollandiae, respectively. Of these, 16 were found to be unannotated novel transcripts and may be taxonomically restricted or orphan genes. Functional enrichment and molecular pathway mapping revealed 13 functionally enriched categories and 11 enriched molecular pathways that were common to the three Macrobrachium species. Pattern of selection analysis revealed 26 genes with signatures of positive selection among pairwise species comparisons. Overall, our results indicate that the same key genes and similar molecular pathways are likely to be involved with freshwater adaptation widely across this decapod group; with nonoverlapping sets of genes showing differential expression (mainly osmoregulatory genes) and signatures of positive selection (genes involved with different life history traits).

Ageing differently: Sex-dependent ageing rates in Daphnia magna

Ageing is defined as the gradual decline of normal physiological functions in a time-dependent manner. Significant progress has been made in characterising the regulatory processes involved in the mechanisms of ageing which would have been hindered without the use of model organisms. Use of alternative model organisms greatly diversifies our understanding of different factors underpinning the ageing process and the potential translation for human application. Unique characteristics make Daphnia an attractive model organism for research into mechanisms underlying ageing, such as transparent body, short generation time, well-characterised methylome, regenerative capabilities and available naturally occurring ecotypes. Most interestingly, genetically identical female and male Daphnia have evolved different average lifespans, providing a unique opportunity for understanding the underlying mechanisms of ageing and regulation of lifespan. Investigating sex differences in longevity could provide insight into principal mechanisms of ageing and lifespan regulation. In this study we provide evidence in support of establishing genetically identical female and male Daphnia as unique and valuable resources for research into mechanisms of ageing and begin to delineate the mechanisms involved in sex differences in lifespan. We identify significant differences between genders in physiological markers such as lifespan, growth rate, heart rate and swimming speed in addition to molecular markers such as lipid peroxidation product accumulation, thiol content decline and age-dependent decline in DNA damage repair efficiency. Overall, our data indicates that investigating sex differences in longevity in the clonal organism Daphnia under controlled laboratory conditions can provide insight into principal mechanisms of ageing and lifespan regulation.

The Evolution of Gene Expression Underlying Vision Loss in Cave Animals

Dissecting the evolutionary genetic processes underlying eye reduction and vision loss in obligate cave-dwelling organisms has been a long-standing challenge in evolutionary biology. Independent vision loss events in related subterranean organisms can provide critical insight into these processes as well as into the nature of convergent loss of complex traits. Advances in evolutionary developmental biology have illuminated the significant role of heritable gene expression variation in the evolution of new forms. Here, we analyze gene expression variation in adult eye tissue across the freshwater crayfish, representing four independent vision-loss events in caves. Species and individual expression patterns cluster by eye function rather than phylogeny, suggesting convergence in transcriptome evolution in independently blind animals. However, this clustering is not greater than what is observed in surface species with conserved eye function after accounting for phylogenetic expectations. Modeling expression evolution suggests that there is a common increase in evolutionary rates in the blind lineages, consistent with a relaxation of selective constraint maintaining optimal expression levels. This is evidence for a repeated loss of expression constraint in the transcriptomes of blind animals and that convergence occurs via a similar trajectory through genetic drift.

De novo European eel transcriptome provides insights into the evolutionary history of duplicated genes in teleost lineages

Paralogues pairs are more frequently observed in eels (Anguilla sp.) than in other teleosts. The paralogues often show low phylogenetic distances; however, they have been assigned to the third round of whole genome duplication (WGD), shared by all teleosts (3R), due to their conserved synteny. The apparent contradiction of low phylogenetic difference and 3R conserved synteny led us to study the duplicated gene complement of the freshwater eels. With this aim, we assembled de novo transcriptomes of two highly relevant freshwater eel species: The European (Anguilla anguilla) and the Japanese eel (Anguilla japonica). The duplicated gene complement was analysed in these transcriptomes, and in the genomes and transcriptomes of other Actinopterygii species. The study included an assessment of neutral genetic divergence (4dTv), synteny, and the phylogenetic origins and relationships of the duplicated gene complements. The analyses indicated a high accumulation of duplications (1217 paralogue pairs) among freshwater eel genes, which may have originated in a WGD event after the Elopomorpha lineage diverged from the remaining teleosts, and thus not at the 3R. However, very similar results were observed in the basal Osteoglossomorpha and Clupeocephala branches, indicating that the specific genomic regions of these paralogues may still have been under tetrasomic inheritance at the split of the teleost lineages. Therefore, two potential hypotheses may explain the results: i) The freshwater eel lineage experienced an additional WGD to 3R, and ii) Some duplicated genomic regions experienced lineage specific rediploidization after 3R in the ancestor to freshwater eels. The supporting/opposing evidence for both hypotheses is discussed.

Multiple stressors produce differential transcriptomic patterns in a stream-dwelling salamander

BACKGROUND

Global biodiversity is decreasing at an alarming rate and amphibians are at the forefront of this crisis. Understanding the factors that negatively impact amphibian populations and effectively monitoring their health are fundamental to addressing this epidemic. Plasma glucocorticoids are often used to assess stress in amphibians and other vertebrates, but these hormones can be extremely dynamic and impractical to quantify in small organisms. Transcriptomic responses to stress hormones in amphibians have been largely limited to laboratory models, and there have been few studies on vertebrates that have evaluated the impact of multiple stressors on patterns of gene expression. Here we examined the gene expression patterns in tail tissues of stream-dwelling salamanders (Eurycea tynerensis) chronically exposed to the stress hormone corticosterone under different temperature regimes.

RESULTS

We found unique transcriptional signatures for chronic corticosterone exposure that were independent of temperature variation. Several of the corticosterone responsive genes are known to be involved in immune system response (LY-6E), oxidative stress (GSTM2 and TRX), and tissue repair (A2M and FX). We also found many genes to be influenced by temperature (CIRBP, HSC71, HSP40, HSP90, HSP70, ZNF593). Furthermore, the expression patterns of some genes (GSTM2, LY-6E, UMOD, ZNF593, CIRBP, HSP90) show interactive effects of temperature and corticosterone exposure, compared to each treatment alone. Through a series of experiments we also showed that stressor induced patterns of expression were largely consistent across ages, life cycle modes, and tissue regeneration.

CONCLUSIONS

Outside of thermal stressors, the application of transcriptomes to monitor the health of non-human vertebrate systems has been vastly underinvestigated. Our study suggests that transcriptomic patterns harbor stressor specific signatures that can be highly informative for monitoring the diverse stressors of amphibian populations.

Lignocellulose degradation in isopods: new insights into the adaptation to terrestrial life

BACKGROUND

Isopods constitute a particular group of crustaceans that has successfully colonized all environments including marine, freshwater and terrestrial habitats. Their ability to use various food sources, especially plant biomass, might be one of the reasons of their successful spread. All isopods, which feed on plants and their by-products, must be capable of lignocellulose degradation. This complex composite is the main component of plants and is therefore an important nutrient source for many living organisms. Its degradation requires a large repertoire of highly specialized Carbohydrate-Active enZymes (called CAZymes) which are produced by the organism itself and in some cases, by its associated microbiota. The acquisition of highly diversified CAZymes could have helped isopods to adapt to their diet and to their environment, especially during land colonization.

RESULTS

To test this hypothesis, isopod host CAZomes (i.e. the entire CAZyme repertoire) were characterized in marine, freshwater and terrestrial species through a transcriptomic approach. Many CAZymes were identified in 64 isopod transcriptomes, comprising 27 de novo datasets. Our results show that marine, freshwater and terrestrial isopods exhibit different CAZomes, illustrating different strategies for lignocellulose degradation. The analysis of variations of the size of CAZy families shows these are expanded in terrestrial isopods while they are contracted in aquatic isopods; this pattern is probably resulting from the evolution of the host CAZomes during the terrestrial adaptation of isopods. We show that CAZyme gene duplications and horizontal transfers can be involved in adaptive divergence between isopod CAZomes.

CONCLUSIONS

Our characterization of the CAZomes in 64 isopods species provides new insights into the evolutionary processes that enabled isopods to conquer various environments, especially terrestrial ones.

Molecular cloning, expression pattern analysis, and in situ hybridization of a Transformer-2 gene in the oriental freshwater prawn, Macrobrachium nipponense (de Haan, 1849)

In this study, we isolated a full-length cDNA sequence from Macrobrachium nipponense and investigated its gene function. We named the gene Mntra-2a because of high similarities and close evolutionary divergence with arthropod tra-2. The full-length cDNA of Mntra-2a was 1293 bp, consisting of a 212 bp 5' UTR, a 268 bp 3' UTR, and an ORF of 813 bp encoding 270 amino acids. It contained an RNA recognition motif and a linker region. Real-time PCR analysis showed that Mntra-2a was highly expressed in the gonads of both males and females. Further in situ hybridization analysis showed that Mntra-2a was mainly located in oocytes and spermatocytes. During embryogenesis, Mntra-2a expression was higher in the cleavage and nauplius stages. During the ovarian reproductive cycle, Mntra-2a expression reached a peak at OvaryV and decreased to the lowest level at OvaryIV. These results indicated that Mntra-2a probably played important roles in embryonic development and early gonad development in M. nipponense. Our results provide basic information for further functional studies of tra-2 in M. nipponense.

Genome-Wide Identification and Expression Profiles of Myosin Genes in the Pacific White Shrimp, Litopenaeus vannamei

As the main structural protein of muscle fiber, myosin is essential for multiple cellular processes or functions, especially for muscle composition and development. Although the shrimp possess a well-developed muscular system, the knowledge about the myosin family in shrimp is far from understood. In this study, we performed comprehensive analysis on the myosin genes in the genome of the Pacific white shrimp, Litopenaeus vannamei. A total of 29 myosin genes were identified, which were classified into 14 subfamilies. Among them, Myo2 subfamily was significantly expanded in the penaeid shrimp genome. Most of the Myo2 subfamily genes were primarily expressed in abdominal muscle, which suggested that Myo2 subfamily genes might be responsible for the well-developed muscular system of the penaeid shrimp. In situ hybridization detection showed that the slow-type muscle myosin gene was mainly localized in pleopod muscle and superficial ventral muscle of the shrimp. This study provides valuable insights into the evolutionary and functional characterization of myosin genes in shrimps, which provides clues for us to understand the well-developed muscular system of shrimp.

Quantitative Proteomic Analysis to Understand the Mechanisms of Zinc Oxide Nanoparticle Toxicity to Daphnia pulex (Crustacea: Daphniidae): Comparing with Bulk Zinc Oxide and Zinc Salt

The widespread use of zinc oxide nanoparticles (ZnO NPs) has resulted in their release to the environment. There has been concern about the ecotoxicity of ZnO NPs, but little is known about their toxic mechanisms. In the present study, we conducted acute toxicity tests to show that ZnO NPs are more toxic to the freshwater crustacean Daphnia pulex compared to bulk ZnO or ZnSO₄·7H₂O. To provide an integrated and quantitative insights into the toxicity of ZnO NPs, we conducted isobaric tags for relative and absolute quantitation (iTRAQ) proteomic analysis, which detected 262, 331, and 360 differentially expressed proteins (DEPs) in D. pulex exposed to ZnO NPs, bulk ZnO, and ZnSO₄·7H₂O, respectively. Among the DEPs, 224 were shared among the three treatments. These proteins were related to energy metabolism, oxidative stress, and endoplasmic reticulum stress. The three forms of Zn all caused D. pulex to downregulate Chitinase expression, disrupt Ca²⁺ homeostasis, and reduce expression of digestive enzymes. Nevertheless, 29 proteins were expressed only in the ZnO NP treatment. In particular, histone (H3) and ribosomal proteins (L13) were obviously influenced under ZnO NP treatment. However, increased expression levels of h3 and l13 genes were not induced only in ZnO NP treatment, they were sensitive to Zn ions under the same exposure concentration. These results indicate that the three zinc substances have a similar mode of action and that released zinc ions are the main contributor to ZnO NP toxicity to D. pulex under a low concentration. Further investigation is needed to clarify whether a small proportion of DEPs or higher bioavailability cause ZnO NPs to be more toxic compared to bulk ZnO or ionic zinc.

The genome of the freshwater water flea Daphnia magna: A potential use for freshwater molecular ecotoxicology

The water flea Daphnia magna is a small planktonic cladoceran. D. magna has been used as a model species for ecotoxicology, as it is sensitive to environmental stressors and environmental changes. Since Daphnia is affected by culture environment and each population/strain has its own ecological and genetic characteristics, its population/strain-based genome information is useful for environmental genomic studies. In this study, we assembled and characterized the genome of D. magna. Using a high-density genetic map of D. magna xinb3, the draft genome was integrated to 10 linkage groups (LGs). The total length of the integrated genome was about 123 Mb with N50 = 10.1 Mb, and the number of scaffolds was 4193 including 10 LGs. A total of 15,721 genes were annotated after manual curation. Orthologous genes were characterized in the genome and compared with other genomes of Daphnia. In addition, we identified defense related genes such as cytochrome P450 (CYP) genes, glutathione S-transferase (GST) genes, and ATP-binding cassette (ABC) genes from the assembled D. magna genome for its potential use in molecular ecotoxicological studies in the freshwater environment. This genomic resource will be helpful to study for a better understanding on molecular mechanism in response to various pollutants.

Variation in transcriptional responses to copper exposure across Daphnia pulex lineages

Copper pollution is pervasive in aquatic habitats and is particularly harmful to invertebrates sensitive to environmental changes such as Daphnia pulex. Mechanisms of toxicity and tolerance to copper are not well understood. We used RNA-sequencing to investigate these mechanisms in three genetically distinct D. pulex clonal lineages with different histories of copper exposure. Upregulated genes after copper exposure were enriched with Gene Ontology (GO) categories involved in digestion, molting and growth, whereas downregulated genes after copper exposure were enriched in the metal-regulatory system, immune response and epigenetic modifications. The three D. pulex clones in our study show largely similar transcriptional patterns in response to copper, with only a total of twenty genes differentially expressed in a single clonal lineages. We also detected lower relative expression of some genes known to be important for copper tolerance, metallothionein and glutathione-S-transferase, in a sensitive lineage sampled from an uncontaminated habitat. Daphnia-specific genes (without orthologs outside the genus) and Daphnia-specific duplications (genes duplicated in the Daphnia lineage) were overrepresented in differentially expressed genes, highlighting an important role for newly emerged genes in tolerating environmental stressors. The results indicate that the D. pulex lineages tested in this study generally respond to copper stress using the same major pathways, but that the more resistant clone with previous copper exposure might be better able to regulate key genes. This finding highlights the important nuances in gene expression among clones, shaped by historical exposure and influencing copper tolerance.

Mini-review an insect-specific system for terrestrialization: Laccase-mediated cuticle formation

Insects are often regarded as the most successful group of animals in the terrestrial environment. Their success can be represented by their huge biomass and large impact on ecosystems. Among the factors suggested to be responsible for their success, we focus on the possibility that the cuticle might have affected the process of insects' evolution. The cuticle of insects, like that of other arthropods, is composed mainly of chitin and structural cuticle proteins. However, insects seem to have evolved a specific system for cuticle formation. Oxidation reaction of catecholamines catalyzed by a copper enzyme, laccase, is the key step in the metabolic pathway for hardening of the insect cuticle. Molecular phylogenetic analysis indicates that laccase functioning in cuticle sclerotization has evolved only in insects. In this review, we discuss a theory on how the insect-specific "laccase" function has been advantageous for establishing their current ecological position as terrestrial animals.

A natural antisense transcript of the fem-1 gene was found expressed in female gonads during the characterization, expression profile, and cellular localization of the fem-1 gene in Pacific white shrimp Penaeus vannamei

The fem-1 gene in Caenorhabditis elegans is involved in sex differentiation; it is specifically required for all aspects of male development. In this study, the full-length cDNA of the fem-1 (Pvfem-1) gene was isolated from the Pacific whiteleg shrimp Penaeus vannamei. The Pvfem-1 transcript is 3778 nt long and encodes a putative protein (PvFEM-1) of 638 amino acids that presented eight ankyrin repeats. The translated protein showed a significant (P < 0.05) structural similitude by superposition with C. elegans FEM-1 protein. Pvfem-1 expression was evaluated by qPCR and in situ hybridization (ISH) during embryogenesis, larval development, and gonads of both genders in subadult and adult life stages. Pvfem-1 was found expressed in brain, intestine, hepatopancreas, and in the gonads of both genders in subadults and adults when quantified by RT-qPCR. A significant finding was the discovery of a natural antisense transcript (NAT) of Pvfem-1 by ISH. It was present in the oocyte nucleus of subadult female shrimp gonads but was not seen within oocytes from adult females, although it was detected in follicular cells, suggesting a possible post-transcriptional regulation of Pvfem-1 in female gonad. Conversely, in males, no NAT was observed, and Pvfem-1 was found expressed in spermatogonia of both, subadult and adult shrimps indicating a function in male sexual differentiation and gametes generation. This study represents the first step for future functional analysis that is expected to contribute to clarifying the role of Pvfem-1 in sex differentiation and determination.

Expression of NK genes that are not part of the NK cluster in the onychophoran Euperipatoides rowelli (Peripatopsidae)

Background

NK genes are a group of homeobox transcription factors that are involved in various molecular pathways across bilaterians. They are typically divided into two subgroups, the NK cluster (NKC) and NK-linked genes (NKL). While the NKC genes have been studied in various bilaterians, corresponding data of many NKL genes are missing to date. To further investigate the ancestral roles of NK family genes, we analyzed the expression patterns of NKL genes in the onychophoran Euperipatoides rowelli.

Results

The NKL gene complement of E. rowelli comprises eight genes, including BarH, Bari, Emx, Hhex, Nedx, NK2.1, vax and NK2.2, of which only NK2.2 was studied previously. Our data for the remaining seven NKL genes revealed expression in different structures associated with the developing nervous system in embryos of E. rowelli. While NK2.1 and vax are expressed in distinct medial regions of the developing protocerebrum early in development, BarH, Bari, Emx, Hhex and Nedx are expressed in late developmental stages, after all major structures of the nervous system have been established. Furthermore, BarH and Nedx are expressed in distinct mesodermal domains in the developing limbs.

Conclusions

Comparison of our expression data to those of other bilaterians revealed similar patterns of NK2.1, vax, BarH and Emx in various aspects of neural development, such as the formation of anterior neurosecretory cells mediated by a conserved molecular mechanism including NK2.1 and vax, and the development of the central and peripheral nervous system involving BarH and Emx. A conserved role in neural development has also been reported from NK2.2, suggesting that the NKL genes might have been primarily involved in neural development in the last common ancestor of bilaterians or at least nephrozoans (all bilaterians excluding xenacoelomorphs). The lack of comparative data for many of the remaining NKL genes, including Bari, Hhex and Nedx currently hampers further evolutionary conclusions. Hence, future studies should focus on the expression of these genes in other bilaterians, which would provide a basis for comparative studies and might help to better understand the role of NK genes in the diversification of bilaterians.

Electronic supplementary material

The online version of this article (10.1186/s12861-019-0185-9) contains supplementary material, which is available to authorized users.

Carbonyl reductase sniffer from the model organism daphnia: Cloning, substrate determination and inhibitory sensitivity

Carbonyl reductases (CRs) represent a fundamental enzymatic defense mechanism against oxidative stress. While commonly two carbonyl reductases (CBR1 and CBR3) are found in mammalian genomes, invertebrate model organisms like Drosophila melanogaster express no CR but a functional homolog to human CBR1, termed sniffer. The importance of sniffer could be demonstrated in D. melanogaster where it protected against age-dependent neurodegeneration. Interestingly, the microcrustacean Daphnia harbors four copies of the CR gene (CR1, CR2, CR3, CR4) in addition to one sniffer gene. Due to this unique equipment Daphnia is an ideal model organism to investigate the function of sniffer. Recombinant sniffer from D. magna und D. pules were produced in E. coli, purified by Ni-affinity chromatography and tested with a variety of aliphatic and aromatic diketones, reactive aldehydes and precursors of advanced glycation end products (AGE). The highest catalytic activities were determined for sniffer from D. pulex with the aromatic dicarbonyls 9,10-phenanthrenequinone (k_cat/K_m = 2.6 s^-1 x μM^-1) and isatin (k_cat/K_m = 1.5 s^-1 x μM^-1). While sniffer from D. magna displayed preference for the same two substances, the respective catalytic activities were noticeably lower. Kinetic constants with aliphatic diketones were generally lower than those with aromatic dicarbonyls for both sniffer enzymes. The best aliphatic diketone as substrate for sniffer from D. magna and D. pulex was hexane-3,4-dione with k_cat/K_m = 0.23 s^-1 μM^-1 and k_cat/K_m = 0.35 s^-1 μM^-1, respectively. Poor or no detectable activity of the two sniffer enzymes was seen with the aliphatic diketones 2,5-hexanedione and 3,5-heptanedione, the aldehydes butanal, hexanal, decanal, crotonaldehyde, acrolein, trans-2-hexenal, and the AGE precursors glyoxal, methylglyoxal, furfural and glyceraldehyde, indicating no physiological function in the metabolism of short-chain aldehydes. Substrate inhibition for both sniffer enzymes was observed with the quinone substrates 1,4-naphthoquinone and 2-methyl-1,4-benzoquinone. From a variety of pesticides endosulfan turned out as an effective inhibitor of the sniffer enzymes (Ki = 9.2 μM for sniffer from D. magna, Ki = 12.0 μM for sniffer from D. pulex). In conclusion, the present results on sniffer from the protein superfamily of the short-chain dehydrogenases/reductases (SDR) in Daphnia ssp. complement earlier studies on carbonyl reductases in the same species and indicate that Daphnia is an interesting model to study the overall response to carbonyl stress.

How clonal are clones? A quest for loss of heterozygosity during asexual reproduction in Daphnia magna

Due to the lack of recombination, asexual organisms are predicted to accumulate mutations and show high levels of within‐individual allelic divergence (heterozygosity); however, empirical evidence for this prediction is largely missing. Instead, evidence of genome homogenization during asexual reproduction is accumulating. Ameiotic crossover recombination is a mechanism that could lead to long genomic stretches of loss of heterozygosity (LOH) and unmasking of mutations that have little or no effect in heterozygous state. Therefore, LOH might be an important force for inducing variation among asexual offspring and may contribute to the limited longevity of asexual lineages. To investigate the genetic consequences of asexuality, here we used high‐throughput sequencing of Daphnia magna for assessing the rate of LOH over a single generation of asexual reproduction. Comparing parthenogenetic daughters with their mothers at several thousand genetic markers generated by restriction site‐associated DNA (RAD) sequencing resulted in high LOH rate estimation that largely overlapped with our estimates for the error rate. To distinguish these two, we Sanger re‐sequenced the top 17 candidate RAD‐loci for LOH, and all of them proved to be false positives. Hence, even though we cannot exclude the possibility that short stretches of LOH occur in genomic regions not covered by our markers, we conclude that LOH does not occur frequently during asexual reproduction in D. magna and ameiotic crossovers are very rare or absent. This finding suggests that clonal lineages of D. magna will remain genetically homogeneous at least over time periods typically relevant for experimental work.

Genome-Wide Stress Responses to Copper and Arsenic in a Field Population of Daphnia

Over the past decade, significant advances have been made to unravel molecular mechanisms of stress response in different ecotoxicological model species. Within this study, we focus on population level transcriptomic responses of a natural population of Daphnia magna Straus, (1820), to heavy metals. We aim to characterize the population level transcriptomic responses, which include standing genetic variation, and improve our understanding on how populations respond to environmental stress at a molecular level. We studied population level responses to two heavy metals, copper and arsenic, and their binary mixture across time. Transcriptomic patterns identified significantly regulated gene families and genes at the population level including cuticle proteins and resilins. Furthermore, some of these differentially regulated gene families, such as cuticle proteins, were also significantly enriched for genetic variations including SNPs and MNPs. In general, genetic variation was observed in specific gene families, many of which are known to be involved in stress response. Overall, our results indicate that molecular stress responses can be identified within natural populations and that linking molecular mechanisms with genetic variation at the population level could contribute significantly to adverse outcome frameworks.

Genome-wide identification and expression of the entire 52 glutathione S-transferase (GST) subfamily genes in the Cu2+-exposed marine copepods Tigriopus japonicus and Paracyclopina nana

In this study, the entire glutathione S-transferases (GSTs), the major phase II detoxification enzyme, were identified in two marine copepod species Tigriopus japonicus and Paracyclopina nana. The genome-wide identification of GSTs in T. japonicus and P. nana resulted in 32 and 20 GSTs in total, respectively. Among the identified GSTs, two specific classes of GSTs, specifically sigma and delta/epsilon GSTs were the dominant form of cytosolic GSTs in T. japonicus, while delta/epsilon and mu classes were dominant cytosolic GSTs in P. nana. In addition, Membrane-Associated Proteins in Eicosanoid and Glutathione metabolism (MAPEG) family were found in relatively higher proportion compared to other classes. Moreover, sigma, delta/epsilon, and microsomal GSTs have shown to expand through tandem duplication. To validate the detoxification function of the identified GSTs, both copepods were exposed to copper (Cu²⁺) and the reactive oxygen species (ROS) level and GST activity were measured. With integration of phylogenetic analysis and xenobiotic-mediated GST mRNA expression patterns along with previous enzymatic activities, the functional divergence among species-specific GST genes was clearly observed. This study covers full identification of GST classes in two marine copepod species and their important role in marine environmental ecotoxicology.

Phototransduction Gene Expression and Evolution in Cave and Surface Crayfishes

In the absence of light in caves, animals have repeatedly evolved reduced eyes and visual systems. Whether the underlying genetic components remain intact in blind species remains unanswered across taxa. The freshwater crayfish have evolved to live in caves multiple times throughout their history; therefore, this system provides an opportunity to probe the genetic patterns and processes underlying repeated vision loss. Using transcriptomic data from the eyes of 14 species of cave and surface crayfishes, we identify the expression of 17 genes putatively related to visual phototransduction. We find a similarly complete repertoire of phototransduction gene families expressed in cave and surface species, but that the expression levels of those transcripts are consistently lower in cave species. We find statistical support for episodic positive selection, increased and decreased selection strength in caves, depending on the gene family. Analyses of gene expression evolution suggest convergent and possibly adaptive downregulation of these genes across eye-reduction events. Our results reveal a combination of evolutionary processes acting on the sequences and gene expression levels of vision-related genes underlying the loss of vision in caves.

Selection of suitable reference genes for gene expression analysis in gills and liver of fish under field pollution conditions

To understand the role of gene expression in adaptive variation, it is necessary to examine expression variation in an ecological context. Quantitative real-time PCR (qPCR) is considered the most accurate and reliable technique to measure gene expression and to validate the data obtained by RNA-seq; however, accurate normalization is crucial. In Chile, the freshwater silverside fish Basilichthys microlepidotus inhabits both polluted and nonpolluted areas, showing differential gene expression related to pollution. In this study, we infer the stability of six potential reference genes (tubulin alpha, hypoxanthine-guanine phosphoribosyltransferase, glyceraldehyde-3-phosphate dehydrogenase, beta-actin, 60S ribosomal protein L13, and 60S ribosomal protein L8) in the gills and liver of silverside individuals inhabiting polluted and nonpolluted areas. To validate the reference genes selected, the most and least stable reference genes were used to normalize two target transcripts, one for each organ. The RefFinder tool was used to analyze and identify the most stably expressed genes. The 60S ribosomal protein L8 gene was ranked as the most stable gene for both organs. Our results show that reference gene selection influences the detection of differences in the expression levels of target genes in different organs and, also highlighting candidate reference genes that could be used in field studies.

Effect of phthalates on development, reproduction, fat metabolism and lifespan in Daphnia magna

Phthalates are used as plasticizers to increase durability, resistivity and flexibility of plastic materials. The commonly used phthalate, diethylhexyl phthalate (DEHP) is used in different plastic materials like food packaging, toys and medical devices. DEHP has been linked to different toxicities in humans as well as in animals, and as a consequence other phthalates, including dibutyl phthalate (DBP) and diethyl phthalate (DEP) are being introduced. The increased use of phthalates has resulted in contamination of aquatic ecosystem and it directly threatens the aquatic life. In this study, we analyzed the effects of three phthalates DEHP, DEP and DBP using freshwater organism Daphnia magna. Although, exposure of the three phthalates at 1 and 10 μM did not result any lethality and hatching delay, the chronic exposure for 14 days resulted in reduction of body length. There was enhanced fat accumulation on exposure to all the phthalates, as indicated by oil red O staining. qRT-PCR analysis of genes involved in fat metabolism suggests that the increase in fat content could be due to inhibition of absorption and catabolism of fatty acids. Reproduction analysis showed that DBP and DEP did not alter fecundity but surprisingly, DEHP at 1 μM increased reproduction by 1.5 fold compared to control group. Phthalates also showed negative effect on lifespan as DEP at 10 μM and DBP at both 1 and 10 μM significantly reduced the lifespan. Our data indicates that along with the banned phthalate DEHP, the other substitute phthalates DEP and DBP could also have detrimental effect on aquatic organisms.

Evaluation of 4-nonylphenol and bisphenol A toxicity using multiple molecular biomarkers in the water flea Daphnia magna

Alkylphenols are well-known endocrine disruptors and may cause developmental and reproductive disorders in aquatic organisms. Daphnia magna is commonly used in ecotoxicological studies as a promising model species to investigate the effects of endocrine distruptors. In the present study, transcriptional modulation of eleven potential molecular indicators related to detoxification, antioxidant, development, and cellular stress was analyzed in D. magna exposed to different concentrations of bisphenol A (BPA) and 4-nonylphenol (4-NP) for 24 h and 48 h, using real-time qPCR. A hierarchical clustering analysis was applied to investigate relations among molecular markers depending on the compound, exposure duration, and concentration. Our findings suggested that GSH-related systems and stress proteins may be involved in cellular defense against BPA and 4-NP-mediated toxicity with different modes of action. Furthermore, these compounds may interrupt molting and reproduction in daphnids. In particular, D. magna GSH-related genes seem to be strongly affected by 4-NP exposure, indicating their potential as molecular biomarkers.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

Genome survey, high-resolution genetic linkage map construction, growth-related quantitative trait locus (QTL) identification and gene location in Scylla paramamosain

Scylla paramamosain is one of the most economically important crabs in China. In this study, the first genome survey sequencing of this crab was performed, and the results revealed that the estimated genome size was 1.21 Gb with high heterozygosity (1.3%). Then, RAD technology was used to construct a high-resolution linkage map for this species. A total of 24,444 single nucleotide polymorphism (SNP) makers were grouped into 47 linkage groups. The total length of the linkage groups was 3087.53 cM with a markers interval of 0.92 cM. With the aid of transcriptome and genome scaffold data, 4,271 markers were linked to genes, including several important growth-related genes such as transforming growth factor-beta regulator I, immune related-gene C-type lectin and ecdysone pathway gene broad-complex-like protein. Further, 442 markers, representing 279 QTLs, associated with 24 traits were identified, and of these markers, 78 were linked to genes. Some interesting genes, such as dedicator of cytokinesis protein 3, tenascin-X and DNA helicase MCM8, were believed to have important relationship with specific traits and merit further exploration. The results of this study will accelerate the genetic improvement and genome sequencing analysis of the mud crab.

IGFBP-rP1, a strongly conserved member of the androgenic hormone signalling pathway in Isopoda

The first protein which has been described to interact with the malacostracan Androgenic Gland Hormone (AGH) is a binding protein called IGFBP-rP1. It has been identified and studied in several species of decapods, in which its interaction with the masculinizing hormone and its expression patterns have been established in several ways. However, this protein remains uncharacterised to date in the other malacostracan orders, like Amphipoda and Isopoda, although they were historically the first ones in which the androgenic gland and the corresponding hormone were respectively described. In this article, we identified the IGFBP-rP1 of isopods and established its implication in the pathway of the AGH with a silencing approach in the model species Armadillidium vulgare. We also showed that this gene is expressed in all the tissues of males and females, with a similar pattern in animals infected with Wolbachia, a feminizing endosymbiont of several isopod species. The expression pattern did not differ during the development of uninfected and infected animals either. We finally studied the evolution of the IGFBP-rP1 in 68 isopod species, looking for conserved motifs and evidence of natural selection. Altogether, our results showed that this gene is constitutively expressed and strongly conserved in isopods, in which it likely constitutes a key element of the insulin/IGF signalling pathway. However, we also illustrated that IGFBP-rP1 is not sufficient on its own to explain the different developmental paths taken by the males and the females or feminized genetic males.

Penaeid shrimp genome provides insights into benthic adaptation and frequent molting

Crustacea, the subphylum of Arthropoda which dominates the aquatic environment, is of major importance in ecology and fisheries. Here we report the genome sequence of the Pacific white shrimp Litopenaeus vannamei, covering ~1.66 Gb (scaffold N50 605.56 Kb) with 25,596 protein-coding genes and a high proportion of simple sequence repeats (>23.93%). The expansion of genes related to vision and locomotion is probably central to its benthic adaptation. Frequent molting of the shrimp may be explained by an intensified ecdysone signal pathway through gene expansion and positive selection. As an important aquaculture organism, L. vannamei has been subjected to high selection pressure during the past 30 years of breeding, and this has had a considerable impact on its genome. Decoding the L. vannamei genome not only provides an insight into the genetic underpinnings of specific biological processes, but also provides valuable information for enhancing crustacean aquaculture.

The Pacific white shrimp Litopenaeus vannamei is an important aquaculture species and a promising model for crustacean biology. Here, the authors provide a reference genome assembly, and show that gene expansion is involved in the regulation of frequent molting as well as benthic adaptation of the shrimp.

Sensory Ecology of Predator-Induced Phenotypic Plasticity

Ecological communities are organized in trophic levels that share manifold interactions forming complex food webs. Infochemicals can further modify these interactions, e.g., by inducing defenses in prey. The micro-crustacean Daphnia is able to respond to predator-specific chemical cues indicating an increased predation risk. Daphnia shows plastic responses by adapting its morphology, behavior, and physiology, increasing organism, and population fitness. This stabilizes community structures. This review will describe the progress that has been made in understanding the high degree of plasticity observed in the model crustacean Daphnia. I summarize current knowledge on the processes of predator detection, ranging from the nature of biologically active chemical cues to the underlying neurophysiological mechanisms. With this, I aim to provide a comprehensive overview on the molecular mechanisms of ad hoc environmental phenotypic adaptation. In times of climate change and pollution understanding information transfer in aquatic systems is valuable as it will allow us to predict whether and how community structures are being affected.

Fitness and Genomic Consequences of Chronic Exposure to Low Levels of Copper and Nickel in Daphnia pulex Mutation Accumulation Lines

In at least some unicellular organisms, mutation rates are temporarily raised upon exposure to environmental stress, potentially contributing to the evolutionary response to stress. Whether this is true for multicellular organisms, however, has received little attention. This study investigated the effects of chronic mild stress, in the form of low-level copper and nickel exposure, on mutational processes in Daphnia pulex using a combination of mutation accumulation, whole genome sequencing and life-history assays. After over 100 generations of mutation accumulation, we found no effects of metal exposure on the rates of single nucleotide mutations and of loss of heterozygosity events, the two mutation classes that occurred in sufficient numbers to allow statistical analysis. Similarly, rates of decline in fitness, as measured by intrinsic rate of population increase and of body size at first reproduction, were negligibly affected by metal exposure. We can reject the possibility that Daphnia were insufficiently stressed to invoke genetic responses as we have previously shown rates of large-scale deletions and duplications are elevated under metal exposure in this experiment. Overall, the mutation accumulation lines did not significantly depart from initial values for phenotypic traits measured, indicating the lineage used was broadly mutationally robust. Taken together, these results indicate that the mutagenic effects of chronic low-level exposure to these metals are restricted to certain mutation classes and that fitness consequences are likely minor and therefore unlikely to be relevant in determining the evolutionary responses of populations exposed to these stressors.

Accelerated rates of large-scale mutations in the presence of copper and nickel

Mutation rate variation has been under intense investigation for decades. Despite these efforts, little is known about the extent to which environmental stressors accelerate mutation rates and influence the genetic load of populations. Moreover, most studies on stressors have focused on unicellular organisms and point mutations rather than large-scale deletions and duplications (copy number variations [CNVs]). We estimated mutation rates in Daphnia pulex exposed to low levels of environmental stressors as well as the effect of selection on de novo mutations. We conducted a mutation accumulation (MA) experiment in which selection was minimized, coupled with an experiment in which a population was propagated under competitive conditions in a benign environment. After an average of 103 generations of MA propagation, we sequenced 60 genomes and found significantly accelerated rates of deletions and duplications in MA lines exposed to ecologically relevant concentrations of metals. Whereas control lines had gene deletion and duplication rates comparable to other multicellular eukaryotes (1.8 × 10-6 per gene per generation), the presence of nickel and copper increased these rates fourfold. The realized mutation rate under selection was reduced to 0.4× that of control MA lines, providing evidence that CNVs contribute to mutational load. Our CNV breakpoint analysis revealed that nonhomologous recombination associated with regions of DNA fragility is the primary source of CNVs, plausibly linking metal-induced DNA strand breaks with higher CNV rates. Our findings suggest that environmental stress, in particular multiple stressors, can have profound effects on large-scale mutation rates and mutational load of multicellular organisms.

Symbiont type and environmental factors affect transcriptome-wide gene expression in the coral Montipora capitata

Reef-building corals may harbor genetically distinct lineages of endosymbiotic dinoflagellates in the genus Symbiodinium, which have been shown to affect important colony properties, including growth rates and resilience against environmental stress. However, the molecular processes underlying these differences are not well understood. In this study, we used whole transcriptome sequencing (RNA-seq) to assess gene expression differences between 27 samples of the coral Montipora capitata predominantly hosting two different Symbiodinium types in clades C and D. The samples were further characterized by their origin from two field sites on Hawai'i Island with contrasting environmental conditions. We found that transcriptome-wide gene expression profiles clearly separated by field site first, and symbiont clade second. With 273 differentially expressed genes (DEGs, 1.3% of all host transcripts), symbiont clade had a measurable effect on host gene expression, but the effect of field site proved almost an order of magnitude higher (1,957 DEGs, 9.6%). According to SNP analysis, we found moderate evidence for host genetic differentiation between field sites (F ST = 0.046) and among corals harboring alternative symbiont clades (F ST = 0.036), suggesting that site-related gene expression differences are likely due to a combination of local adaptation and acclimatization to environmental factors. The correlation between host gene expression and symbiont clade may be due to several factors, including host genotype or microhabitat selecting for alternative clades, host physiology responding to different symbionts, or direct modulation of host gene expression by Symbiodinium. However, the magnitude of these effects at the level of transcription was unexpectedly small considering the contribution of symbiont type to holobiont phenotype.

Complementary roles of photoperiod and temperature in environmental sex determination in Daphnia spp

Daphnia spp, a keystone genus in freshwater lentic habitats, are subject to environmental sex determination wherein environmental conditions dictate offspring sex and whether they reproduce asexually or sexually. The introduction of males into a population denotes the first step in the switch from asexual parthenogenetic reproduction to sexual reproduction. We tested the hypothesis that photoperiod and temperature co-regulate male sex determination and that these environmental stimuli would activate elements of the male sex determination signaling cascade. Results revealed that photoperiod was a critical cue in creating permissive conditions for male production. Further, under photoperiod-induced permissive conditions, male sex determination was temperature dependent. The two daphnid species evaluated, Daphnia pulex and D. magna, exhibited different temperature dependencies. D. pulex produced fewer males with increasing temperatures between 16-22°C, and D. magna exhibited the opposite trend. We found consistent expression patterns of key genes along the male sex determining signaling pathway in D. pulex independent of environmental stimuli. mRNA levels for the enzyme responsible for synthesis of the male sex determining hormone, methyl farnesoate, were elevated early in the reproductive cycle, followed by increased mRNA levels of the methyl farnesoate receptor subunits, Met and SRC. Environmental conditions that stimulated male offspring production significantly increased Met mRNA levels. Results indicate that male sex determination in daphnids is under the permissive control of photoperiod and the regulatory control of temperature. Further, these environmental cues may stimulate male sex determination by increasing levels of the Met subunit of the methyl farnesoate receptor.

The transcriptome of the marine calanoid copepod Temora longicornis under heat stress and recovery

Understanding the impacts of global change in zooplankton communities is crucial, as alterations in the zooplankton communities can affect entire marine ecosystems. Despite the economic and ecological importance of the calanoid copepod Temora longicornis in the Belgian part of the North Sea, molecular data is still very limited for this species. Using HiSeq Illumina sequencing, we sequenced the whole transcriptome of T. longicornis, after being exposed to realistic temperatures of 14 and 17 °C. After both an acute (1 day) and a more sustained (5 days) thermal exposure to 17 °C, we investigated gene expression differences with animals exposed to 14 °C, which may be critical for the thermal acclimation and resilience of this copepod species. We also studied the possibility of a short term stress recovery of a heat shock. A total of 179,569 transcripts were yielded, of which 44,985 putative ORF transcripts were identified. These transcripts were subsequently annotated into roughly 22,000 genes based on known sequences using Gene Ontology (GO) and KEGG databases. Temora only showed a mild response to both the temperature and the duration of the exposure. We found that the expression of 27 transcripts varied significantly with an increase in temperature of 3 °C, of which eight transcripts were differentially expressed after acute exposure only. Gene set enrichment analysis revealed that, overall, T. longicornis was more impacted by a sustained thermal exposure, rather than an immediate (acute) exposure, with two times as many enriched GO terms in the sustained treatment. We also identified several general stress responses independent of exposure time, such as modified protein synthesis, energy mobilisation, cuticle and chaperone proteins. Finally, we highlighted candidate genes of a possible recovery from heat exposure, identifying similar terms as those enriched in the heat treatments, i.e. related to for example energy metabolism, cuticle genes and extracellular matrix. The data presented in this study provides the first transcriptome available for T. longicornis which can be used for future genomic studies.

Daphnia galeata responds to the exposure to an ichthyosporean gut parasite by down-regulation of immunity and lipid metabolism

BACKGROUND

Regulatory circuits of infection in the emerging experimental model system, water flea Daphnia and their microparasites, remain largely unknown. Here we provide the first molecular insights into the response of Daphnia galeata to its highly virulent and common parasite Caullerya mesnili, an ichthyosporean that infects the gut epithelium. We generated a transcriptomic dataset using RNAseq from parasite-exposed (vs. control) Daphnia, at two time points (4 and 48 h) after parasite exposure.

RESULTS

We found a down-regulation of metabolism and immunity-related genes, at 48 h (but not 4 h) after parasite exposure. These genes are involved in lipid metabolism and fatty acid biosynthesis, as well as microbe recognition (e.g. c-type lectins) and pathogen attack (e.g. gut chitin).

CONCLUSIONS

General metabolic suppression implies host energy shift from reproduction to survival, which is in agreement with the known drastic reduction in Daphnia fecundity after Caullerya infection. The down-regulation of gut chitin indicates a possible interaction between the peritrophic matrix and the evading host immune system. Our study provides the first description of host transcriptional responses in this very promising host-parasite experimental system.

Comparative Genomics Reveals Thousands of Novel Chemosensory Genes and Massive Changes in Chemoreceptor Repertories across Chelicerates

Chemoreception is a widespread biological function that is essential for the survival, reproduction, and social communication of animals. Though the molecular mechanisms underlying chemoreception are relatively well known in insects, they are poorly studied in the other major arthropod lineages. Current availability of a number of chelicerate genomes constitutes a great opportunity to better characterize gene families involved in this important function in a lineage that emerged and colonized land independently of insects. At the same time, that offers new opportunities and challenges for the study of this interesting animal branch in many translational research areas. Here, we have performed a comprehensive comparative genomics study that explicitly considers the high fragmentation of available draft genomes and that for the first time included complete genome data that cover most of the chelicerate diversity. Our exhaustive searches exposed thousands of previously uncharacterized chemosensory sequences, most of them encoding members of the gustatory and ionotropic receptor families. The phylogenetic and gene turnover analyses of these sequences indicated that the whole-genome duplication events proposed for this subphylum would not explain the differences in the number of chemoreceptors observed across species. A constant and prolonged gene birth and death process, altered by episodic bursts of gene duplication yielding lineage-specific expansions, has contributed significantly to the extant chemosensory diversity in this group of animals. This study also provides valuable insights into the origin and functional diversification of other relevant chemosensory gene families different from receptors, such as odorant-binding proteins and other related molecules.

The genome of the water strider Gerris buenoi reveals expansions of gene repertoires associated with adaptations to life on the water

BACKGROUND

Having conquered water surfaces worldwide, the semi-aquatic bugs occupy ponds, streams, lakes, mangroves, and even open oceans. The diversity of this group has inspired a range of scientific studies from ecology and evolution to developmental genetics and hydrodynamics of fluid locomotion. However, the lack of a representative water strider genome hinders our ability to more thoroughly investigate the molecular mechanisms underlying the processes of adaptation and diversification within this group.

RESULTS

Here we report the sequencing and manual annotation of the Gerris buenoi (G. buenoi) genome; the first water strider genome to be sequenced thus far. The size of the G. buenoi genome is approximately 1,000 Mb, and this sequencing effort has recovered 20,949 predicted protein-coding genes. Manual annotation uncovered a number of local (tandem and proximal) gene duplications and expansions of gene families known for their importance in a variety of processes associated with morphological and physiological adaptations to a water surface lifestyle. These expansions may affect key processes associated with growth, vision, desiccation resistance, detoxification, olfaction and epigenetic regulation. Strikingly, the G. buenoi genome contains three insulin receptors, suggesting key changes in the rewiring and function of the insulin pathway. Other genomic changes affecting with opsin genes may be associated with wavelength sensitivity shifts in opsins, which is likely to be key in facilitating specific adaptations in vision for diverse water habitats.

CONCLUSIONS

Our findings suggest that local gene duplications might have played an important role during the evolution of water striders. Along with these findings, the sequencing of the G. buenoi genome now provides us the opportunity to pursue exciting research opportunities to further understand the genomic underpinnings of traits associated with the extreme body plan and life history of water striders.

De novo transcriptome assembly of brackish water flea Diaphanosoma celebensis based on short-term cadmium and benzo[a]pyrene exposure experiments

To develop a brackish water flea as a promising model for marine monitoring, Diaphanosoma celebensis were exposed to two pollutants, cadmium (Cd) and benzo[a]pyrene (BaP), which have different chemical characteristics and distinct modes of metabolic action on aquatic animals. Twenty-four hours after exposure to Cd (2 mg/L) or BaP (25 μg/L), whole body transcriptomes were analyzed. In total, 99.6 Mbp were assembled from nine libraries, resulting in 98,458 transcripts with an N50 of 1883 bp and an average contig length of 968 bp. Functional gene annotations were performed using Gene Ontology, Eukaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes pathway analyses. Cd significantly modulated endocrine and digestive enzyme system. Following BaP treatment, DNA repair and circadian rhythm related metabolisms were significantly modulated. Both the chemicals induced stress response and detoxification metabolism. This brackish water flea genomic information will be useful to monitor estuaries and coastal regions, as water fleas have been confirmed as promising sentinel models in freshwater ecosystems.

Proteomic Responses Under Cold Stress Reveal Unique Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei), one of the most widely cultured shrimp species in the world, often suffers from cold stress. To understand the molecular mechanism of cold tolerance in Pacific white shrimp, we conducted a proteomic analysis on two contrasting shrimp cultivars, namely, cold-tolerant Guihai2 (GH2) and cold-sensitive Guihai1 (GH1), under normal temperature (28°C), under cold stress (16°C), and during recovery to 28°C. In total, 3,349 proteins were identified, among which 2,736 proteins were quantified. Based on gene ontology annotations, differentially expressed proteins largely belonged to biological processes, cellular components, and molecular functions. KEGG pathway annotations indicated that the main changes were observed in the lysosome, ribosomes, and oxidative phosphorylation. Subcellular localization analysis showed a significant increase in proteins present in cytosol, extracellular regions, and mitochondria. Combining enrichment-based clustering analysis and qRT-PCR analysis, we found that glutathione S-transferase, zinc proteinase, m7GpppX diphosphatase, AP2 transcription complex, and zinc-finger transcription factors played a major role in the cold stress response in Pacific white shrimp. Moreover, structure proteins, including different types of lectin and DAPPUDRAFT, were indispensable for cold stress tolerance of the Pacific white shrimp. Results indicate the molecular mechanisms of the Pacific white shrimp in response to cold stress and provide new insight into breeding new cultivars with increased cold tolerance.

Evolutionary development and morphological modifications of the brain: an interview with Angelika Stollewerk

Angelika Stollewerk is a Reader at Queen Mary University of London, where her lab uses a diverse range of species to study the evolution of the arthropod nervous system. Angelika spoke to us about social spiders, the future of evo-devo, and open peer review.

Ecdysteroid and juvenile hormone biosynthesis, receptors and their signaling in the freshwater microcrustacean Daphnia

The two essential insect hormones, ecdysteroids and juvenile hormones, are possessed not only by insects, but also widely by arthropods, and regulate various developmental and physiological processes. In contrast to the abundant information about molecular endocrine mechanisms in insects, the knowledge of non-insect arthropod endocrinology is still limited. In this review, we summarize recent reports about the molecular basis of these two major insect hormones in the freshwater microcrustacean Daphnia, a keystone taxon in limnetic ecology and a bioindicator in environmental studies. Comprehensive comparisons of endocrine signaling pathways between insects and daphnids may shed light on the regulatory mechanisms of various biological phenomena and, moreover, evolutionary processes of arthropod species.

A New Standard for Crustacean Genomes: The Highly Contiguous, Annotated Genome Assembly of the Clam Shrimp Eulimnadia texana Reveals HOX Gene Order and Identifies the Sex Chromosome

Vernal pool clam shrimp (Eulimnadia texana) are a promising model system due to their ease of lab culture, short generation time, modest sized genome, a somewhat rare stable androdioecious sex determination system, and a requirement to reproduce via desiccated diapaused eggs. We generated a highly contiguous genome assembly using 46× of PacBio long read data and 216× of Illumina short reads, and annotated using Illumina RNAseq obtained from adult males or hermaphrodites. Of the 120 Mb genome 85% is contained in the largest eight contigs, the smallest of which is 4.6 Mb. The assembly contains 98% of transcripts predicted via RNAseq. This assembly is qualitatively different from scaffolded Illumina assemblies: It is produced from long reads that contain sequence data along their entire length, and is thus gap free. The contiguity of the assembly allows us to order the HOX genes within the genome, identifying two loci that contain HOX gene orthologs, and which approximately maintain the order observed in other arthropods. We identified a partial duplication of the Antennapedia complex adjacent to the few genes homologous to the Bithorax locus. Because the sex chromosome of an androdioecious species is of special interest, we used existing allozyme and microsatellite markers to identify the E. texana sex chromosome, and find that it comprises nearly half of the genome of this species. Linkage patterns indicate that recombination is extremely rare and perhaps absent in hermaphrodites, and as a result the location of the sex determining locus will be difficult to refine using recombination mapping.