The ecoresponsive genome of Daphnia pulex

High-quality genome of Diaphanosoma dubium provides insights into molecular basis of its broad ecological adaptation

Diaphanosoma dubium Manuilova, 1964, is a widespread planktonic water flea in Asian freshwater. Although sharing similar ecological roles with species of Daphnia, studies on D. dubium and its congeners are still few and lacking a genome for the further studies. Here, we assembled a high quality and chromosome level genome of D. dubium by combining long reads sequencing and Hi-C technologies. The total length of assembled genome was 101.8 Mb, with 98.92 Mb (97.2%) anchored into 22 chromosomes. Through comparative genomic analysis, we found the genes, involved in anti-ROS, detoxification, protein digestion, germ cells regulation and protection, underwent expansion in D. dubium. These genes and their expansion helpfully explain its widespread geographical distribution and dominance in eutrophic waters. This study provides insight into the adaptive evolution of D. dubium at genomic perspectives, and the present high quality genomic resource will be a footstone for future omics studies of the species and its congeners.

Genome Assembly of a Relict Arabian Species of Daphnia O. F. Müller (Crustacea: Cladocera) Adapted to the Desert Life

The water flea Daphnia O.F. Müller 1776 (Crustacea: Cladocera) is an important model of recent evolutionary biology. Here, we report a complete genome of Daphnia (Ctenodaphnia) arabica (Crustacea: Cladocera), recently described species endemic to deserts of the United Arab Emirates. In this study, genome analysis of D. arabica was carried out to investigate its genomic differences, complexity as well as its historical origins within the subgenus Daphnia (Ctenodaphnia). Hybrid genome assembly of D. arabica resulted in ~116 Mb of the assembled genome, with an N50 of ~1.13 Mb (BUSCO score of 99.2%). From the assembled genome, in total protein coding, 5374 tRNA and 643 rRNA genes were annotated. We found that the D. arabica complete genome differed from those of other Daphnia species deposited in the NCBI database but was close to that of D. cf. similoides. However, its divergence time estimate sets D. arabica in the Mesozoic, and our demographic analysis showed a great reduction in its genetic diversity compared to other Daphnia species. Interestingly, the population expansion in its diversity occurred during the megadrought climate around 100 Ka ago, reflecting the adaptive feature of the species to arid and drought-affected environments. Moreover, the PFAM comparative analysis highlights the presence of the important domain SOSS complex subunit C in D. arabica, which is missing in all other studied species of Daphnia. This complex consists of a few subunits (A, B, C) working together to maintain the genome stability (i.e., promoting the reparation of DNA under stress). We propose that this domain could play a role in maintaining the fitness and survival of this species in the desert environment. The present study will pave the way for future research to identify the genes that were gained or lost in this species and identify which of these were key factors to its adaptation to the harsh desert environment.

Assembly of a Draft Genome for the Mouse Ectoparasite Myocoptes musculinus

Myocoptes musculinus is a common ectoparasite of wild mice and is occasionally found on research mice. Infestations of research mice are often subclinical but can cause severe dermatitis. Perhaps more importantly, infestations can cause immunologic reactions that may alter research outcomes, and most animal research facilities strive to prevent or eliminate mites from their mouse colonies. M. musculinus infestations are currently detected by using microscopic evaluation of the fur and skin and PCR assays of pelt swabs targeting the rRNA genes of this mite. In our facility, we encountered multiple, false-positive 18S rRNA PCR results from a closed mouse colony. We could not identify the source of the false positives even after performing PCR analysis of other Myocoptes gene targets using assays developed from the few other target genomic sequences available for M. musculinus or Myocoptes japonensis in public databases. This situation highlighted the limited genetic resources available for development of diagnostic tests specific for this ectoparasite. To expand the available genetic resources, we generated a metagenome of M. musculinus derived by sequencing from fur plucks of an infected mouse. We also determined the completeness of this metagenome and compared it with those of related mites.

Pilot testing and optimization of a larval fathead minnow high throughput transcriptomics assay

Concentrations at which global gene expression profiles in cells or animals exposed to a test substance start to differ significantly from those of controls have been proposed as an alternative point of departure for use in screening level hazard assessment. The present study describes pilot testing of a high throughput compatible transcriptomics assay with larval fathead minnows. One day post hatch fathead minnows were exposed to eleven different concentrations of three metals, three selective serotonin reuptake inhibitors, and four neonicotinoid-like compounds for 24 h and concentration response modeling was applied to whole body gene expression data. Transcriptomics-based points of departure (tPODs) were consistently lower than effect concentrations reported in apical endpoint studies in fish. However, larval fathead minnow-based tPODs were not always lower than concentrations reported to elicit apical toxicity in other aquatic organisms like crustaceans or insects. Random in silico subsampling of data from the pilot assays was used to evaluate various assay design and acceptance considerations such as transcriptome coverage, number of replicate individuals to sequence per treatment, and minimum number of differentially expressed genes to produce a reliable tPOD estimate. Results showed a strong association between the total number of genes for which a concentration response relationship could be derived and the overall variability in the resulting tPOD estimates. We conclude that, for our current assay design and analysis pipeline, tPODs based on fewer than 15 differentially expressed genes are likely to be unreliable for screening and that interindividual variability in gene expression profiles appears to be a more significant driver of tPOD variability than sample size alone. Results represent initial steps toward developing high throughput transcriptomics assays for use in ecological hazard screening.

Deep Diversity: Extensive Variation in the Components of Complex Visual Systems across Animals

Understanding the molecular underpinnings of the evolution of complex (multi-part) systems is a fundamental topic in biology. One unanswered question is to what the extent do similar or different genes and regulatory interactions underlie similar complex systems across species? Animal eyes and phototransduction (light detection) are outstanding systems to investigate this question because some of the genetics underlying these traits are well characterized in model organisms. However, comparative studies using non-model organisms are also necessary to understand the diversity and evolution of these traits. Here, we compare the characteristics of photoreceptor cells, opsins, and phototransduction cascades in diverse taxa, with a particular focus on cnidarians. In contrast to the common theme of deep homology, whereby similar traits develop mainly using homologous genes, comparisons of visual systems, especially in non-model organisms, are beginning to highlight a "deep diversity" of underlying components, illustrating how variation can underlie similar complex systems across taxa. Although using candidate genes from model organisms across diversity was a good starting point to understand the evolution of complex systems, unbiased genome-wide comparisons and subsequent functional validation will be necessary to uncover unique genes that comprise the complex systems of non-model groups to better understand biodiversity and its evolution.

Short- and long-term single and combined effects of microplastics and chromium on the freshwater water flea Daphnia magna

In this study, we investigated the individual and combined effects of microplastics (MPs) and chromium (Cr) on the freshwater water flea Daphnia magna by measuring mortality, bioaccumulation, antioxidative response, multixenobiotic resistance activity, and sestrin-related mitochondrial biogenesis in short-term assays and in vivo endpoints including reproduction and adult survival rate in long-term assays. Exposure to MPs, Cr, and their combination caused significant deleterious effects and acute toxicity in D. magna. Alterations in oxidative stress occurred in the groups treated with MPs and Cr alone and together. However, upon co-exposure to MPs, the Cr concentration, measured by inductively coupled plasma optical emission spectroscopy, decreased, suggesting that MPs and Cr interact with each other. Based on enzymatic activities, we noted a decrease in MP egestion via inhibition of P-glycoprotein activity in the MP-exposed groups, and multidrug resistance-associated protein activity increased in some of the MP-exposed animals depending on Cr concentration. On the other hand, MP exposure seemed to lead to mitochondrial transcription dysfunction induced by Cr via sestrin-related mitochondrial biogenesis. Overall, these results indicate that co-exposure to MPs and Cr causes acute toxicity in D. magna but lacks the chronic toxicity (21 days) and mitochondrial dysfunction caused by Cr exposure alone.

The diversity of invertebrate visual opsins spanning Protostomia, Deuterostomia, and Cnidaria

Across eumetazoans, the ability to perceive and respond to visual stimuli is largely mediated by opsins, a family of proteins belonging to the G protein-coupled receptor (GPCR) superclass. Lineage-specific gains and losses led to a striking diversity in the numbers, types, and spectral sensitivities conferred by visual opsin gene expression. Here, we review the diversity of visual opsins and differences in opsin gene expression from well-studied protostome, invertebrate deuterostome, and cnidarian groups. We discuss the functional significance of opsin expression differences and spectral tuning among lineages. In some cases, opsin evolution has been linked to the detection of relevant visual signals, including sexually selected color traits and host plant features. In other instances, variation in opsins has not been directly linked to functional or ecological differences. Overall, the array of opsin expression patterns and sensitivities across invertebrate lineages highlight the diversity of opsins in the eumetazoan ancestor and the labile nature of opsins over evolutionary time.

Degradation of Azithromycin from aqueous solution using Chlorine-ferrous- oxidation: ANN-GA modeling and Daphnia magna biotoxicity test assessment

Azithromycin (AZM), an antibacterial considered one of the most consumed drugs, especially during the period against the Covid 19 pandemic, and it is one of the persistent contaminants that can be released into aquatic ecosystems. The purpose of this study is to determine the efficacy of a Fenton-like process (chlorine/iron) for the degradation of AZM in an aqueous medium by determining the impact of several factors (the initial concentration of (FeSO₄, NaClO, pollutant), and the initial pH) on the degradation rate. The Response Surface Methodology (RSM) based on the Box-Wilson design as well as the Artificial Neural Network (ANN) modeling combined with a genetic algorithm (GA) approaches were used to determine the optimal levels of the selected variables and the optimal rate of degradation. The quadratic model of multi-linear regression developed indicated that the optimal conditions were a concentration of chlorine of 600 μM, the concentration of AZM is 32.8 mg/L, the mass of the catalyst FeSO₄ is 3.5 mg and a pH of 2.5, these optimal values gave a predicted and experimental yield of 64.05% and 70% respectively, the lack of fit test in RSM modeling (F₀ = 3.31 which is inferior to Fcritic (0.05, 10.4) = 5.96) indicates that the true regression function is not linear therefore, the ANN-GA modeling as non-linear regression indicated that the optimal conditions were a concentration of chlorine of 256 μM, the concentration of AZM is 5 mg/L, the mass of the catalyst FeSO₄ is 9.5 mg and a pH of 2.8, these optimal values gave a predicted and experimental yield of 79.69% and close to 80% respectively, Furthermore, biotoxicity tests were conducted to confirm the performance of our process using bio-indicators called daphnia (Daphnia magna), which demonstrated the efficacy of the like-Fenton process after 4 h of degradation.

Crustacean conundrums: a review of opsin diversity and evolution

Knowledge of crustacean vision is lacking compared to the more well-studied vertebrates and insects. While crustacean visual systems are typically conserved morphologically, the molecular components (i.e. opsins) remain understudied. This review aims to characterize opsin diversity across crustacean lineages for an integrated view of visual system evolution. Using publicly available data from 95 species, we identified opsin sequences and classified them by clade. Our analysis produced 485 putative visual opsins and 141 non-visual opsins. The visual opsins were separated into six clades: long wavelength sensitive (LWS), middle wavelength sensitive (MWS) 1 and 2, short wavelength or ultraviolet sensitive (SWS/UVS) and a clade of thecostracan opsins, with multiple LWS and MWS opsin copies observed. The SWS/UVS opsins were relatively conserved in most species. The crustacean classes Cephalocarida, Remipedia and Hexanauplia exhibited reduced visual opsin diversity compared to others, with the malacostracan decapods having the highest opsin diversity. Non-visual opsins were identified from all investigated classes except Cephalocarida. Additionally, a novel clade of non-visual crustacean-specific, R-type opsins (Rc) was discovered. This review aims to provide a framework for future research on crustacean vision, with an emphasis on the need for more work in spectral characterization and molecular analysis. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

Daphnia as a Sentinel Species for Environmental Health Protection: A Perspective on Biomonitoring and Bioremediation of Chemical Pollution

Despite available technology and the knowledge that chemical pollution damages human and ecosystem health, chemical pollution remains rampant, ineffectively monitored, rarely prevented, and only occasionally mitigated. We present a framework that helps address current major challenges in the monitoring and assessment of chemical pollution by broadening the use of the sentinel species Daphnia as a diagnostic agent of water pollution. And where prevention has failed, we propose the application of Daphnia as a bioremediation agent to help reduce hazards from chemical mixtures in the environment. By applying "omics" technologies to Daphnia exposed to real-world ambient chemical mixtures, we show improvements at detecting bioactive components of chemical mixtures, determining the potential effects of untested chemicals within mixtures, and identifying targets of toxicity. We also show that using Daphnia strains that naturally adapted to chemical pollution as removal agents of ambient chemical mixtures can sustainably improve environmental health protection. Expanding the use of Daphnia beyond its current applications in regulatory toxicology has the potential to improve both the assessment and the remediation of environmental pollution.

Toxicity by descent: A comparative approach for chemical hazard assessment

Toxicology is traditionally divided between human and eco-toxicology. In the shared pursuit of environmental health, this separation does not account for discoveries made in the comparative studies of animal genomes. Here, we provide evidence on the feasibility of understanding the health impact of chemicals on all animals, including ecological keystone species and humans, based on a significant number of conserved genes and their functional associations to health-related outcomes across much of animal diversity. We test four conditions to understand the value of comparative genomics data to inform mechanism-based human and environmental hazard assessment: (1) genes that are most fundamental for health evolved early during animal evolution; (2) the molecular functions of pathways are better conserved among distantly related species than the individual genes that are members of these pathways; (3) the most conserved pathways among animals are those that cause adverse health outcomes when disrupted; (4) gene sets that serve as molecular signatures of biological processes or disease-states are largely enriched by evolutionarily conserved genes across the animal phylogeny. The concept of homology is applied in a comparative analysis of gene families and pathways among invertebrate and vertebrate species compared with humans. Results show that over 70% of gene families associated with disease are shared among the greatest variety of animal species through evolution. Pathway conservation between invertebrates and humans is based on the degree of conservation within vertebrates and the number of interacting genes within the human network. Human gene sets that already serve as biomarkers are enriched by evolutionarily conserved genes across the animal phylogeny. By implementing a comparative method for chemical hazard assessment, human and eco-toxicology converge towards a more holistic and mechanistic understanding of toxicity disrupting biological processes that are important for health and shared among animals (including humans).

The interplay between prior selection, mild intermittent exposure, and acute severe exposure in phenotypic and transcriptional response to hypoxia

Hypoxia has profound and diverse effects on aerobic organisms, disrupting oxidative phosphorylation and activating several protective pathways. Predictions have been made that exposure to mild intermittent hypoxia may be protective against more severe exposure and may extend lifespan. Here we report the lifespan effects of chronic, mild, intermittent hypoxia, and short-term survival in acute severe hypoxia in four clones of Daphnia magna originating from either permanent or intermittent habitats. We test the hypothesis that acclimation to chronic mild intermittent hypoxia can extend lifespan through activation of antioxidant and stress-tolerance pathways and increase survival in acute severe hypoxia through activation of oxygen transport and storage proteins and adjustment to carbohydrate metabolism. Unexpectedly, we show that chronic hypoxia extended the lifespan in the two clones originating from intermittent habitats but had the opposite effect in the two clones from permanent habitats, which also showed lower tolerance to acute hypoxia. Exposure to chronic hypoxia did not protect against acute hypoxia; to the contrary, Daphnia from the chronic hypoxia treatment had lower acute hypoxia tolerance than normoxic controls. Few transcripts changed their abundance in response to the chronic hypoxia treatment in any of the clones. After 12 h of acute hypoxia treatment, the transcriptional response was more pronounced, with numerous protein-coding genes with functionality in oxygen transport, mitochondrial and respiratory metabolism, and gluconeogenesis, showing upregulation. While clones from intermittent habitats showed somewhat stronger differential expression in response to acute hypoxia than those from permanent habitats, contrary to predictions, there were no significant hypoxia-by-habitat of origin or chronic-by-acute treatment interactions. GO enrichment analysis revealed a possible hypoxia tolerance role by accelerating the molting cycle and regulating neuron survival through upregulation of cuticular proteins and neurotrophins, respectively.

Endocrine-disrupting potential and toxicological effect of para-phenylphenol on Daphnia magna

Several phenol derivatives are suspected endocrine disruptors and have received attention in risk assessment studies for several decades owing to the structural similarity between estrogens and phenolic compounds. We assessed the endocrine disrupting effect of the phenolic compound para-phenylphenol (PPP) through acute tests and evaluating chronic endpoints in an invertebrate model, Daphnia magna. Exposure of D. magna to PPP induced substantial adverse effects, namely, reduced fecundity, slowed growth rate, delayed first brood, and a reduction in neonate size. Furthermore, we investigated the mRNA expression of relevant genes to elucidate the mechanism of endocrine disruption by PPP. Exposure of D. magna to PPP induced the substantial downregulation of genes and markers related to reproduction and development, such as EcR-A, EcR-B, Jhe, and Vtg. Consequently, we demonstrated that PPP has an endocrine disrupting effect on reproduction and development in D. magna.

Morphological and molecular evolution of hadal amphipod’s eggs provides insights into embryogenesis under high hydrostatic pressure

Hadal zones are unique habitats characterized by high hydrostatic pressure (HHP) and scarce food supplies. The ability of eggs of species dwelling in hadal zones to develop into normal embryo under high hydrostatic pressure is an important evolutionary and developmental trait. However, the mechanisms underlying the development of eggs of hadal-dwelling species remain unknown due to the difficulty of sampling ovigerous females. Here, morphological and transcriptome analyses of eggs of the “supergiant” amphipod Alicella gigantea collected from the New Britain Trench were conducted. The morphology of A. gigantea eggs, including size, was assessed and the ultrastructure of the eggshell was investigated by scanning electron microscopy. Transcriptome sequencing and molecular adaptive evolution analysis of A. gigantea eggs showed that, as compared with shallow-water Gammarus species, genes exhibiting accelerated evolution and the positively selected genes were mostly related to pathways associated with “mitosis” and “chitin-based embryonic cuticle biosynthetic process”, suggesting that “normal mitosis maintenance” and “cuticle development and protection” are the two main adaptation strategies for survival of eggs in hadal environments. In addition, the concentration of trimethylamine oxide (TMAO), an important osmotic regulator, was significantly higher in the eggs of hadal amphipods as compared to those of shallow-water species, which might promote the eggs’ adaptation abilities. Morphological identification, evolutionary analysis, and the trimethylamine oxide concentration of A. gigantea eggs will facilitate a comprehensive overview of the piezophilic adaptation of embryos in hadal environments and provide a strategy to analyze embryogenesis under high hydrostatic pressure.

Thermal stress and mutation accumulation increase heat shock protein expression in Daphnia

Understanding the short- and long-term consequences of climate change is a major challenge in biology. For aquatic organisms, temperature changes and drought can lead to thermal stress and habitat loss, both of which can ultimately lead to higher mutation rates. Here, we examine the effect of high temperature and mutation accumulation on gene expression at two loci from the heat shock protein (HSP) gene family, HSP60 and HSP90. HSPs have been posited to serve as 'mutational capacitors' given their role as molecular chaperones involved in protein folding and degradation, thus buffering against a wide range of cellular stress and destabilization. We assayed changes in HSP expression across 5 genotypes of Daphnia magna, a sentinel species in ecology and environmental biology, with and without acute exposure to thermal stress and accumulated mutations. Across genotypes, HSP expression increased ~ 6× in response to heat and ~ 4× with mutation accumulation, individually. Both factors simultaneously (lineages with high mutation loads exposed to high heat) increased gene expression ~ 23×-much more than that predicted by an additive model. Our results corroborate suggestions that HSPs can buffer against not only the effects of heat, but also mutations-a combination of factors both likely to increase in a warming world.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10682-022-10209-1.

Using the Daphnia magna Transcriptome to Distinguish Water Source: Wetland and Stormwater Case Studies

A major challenge in ecotoxicology is accurately and sufficiently measuring chemical exposures and biological effects given the presence of complex and dynamic contaminant mixtures in surface waters. It is impractical to quantify all chemicals in such matrices over space and time, and even if it were practical, concomitant biological effects would not be elucidated. Our study examined the performance of the Daphnia magna transcriptome to detect distinct responses across three water sources in Minnesota: laboratory (well) waters, wetland waters, and storm waters. Pyriproxyfen was included as a gene expression and male neonate production positive control to examine whether gene expression resulting from exposure to this well-studied juvenoid hormone analog can be detected in complex matrices. Laboratory-reared (<24 h) D. magna were exposed to a water source and/or pyriproxyfen for 16 days to monitor phenotypic changes or 96 h to examine gene expression responses using Illumina HiSeq 2500 (10 million reads per library, 50-bp paired end [2 × 50]). The results indicated that a unique gene expression profile was produced for each water source. At 119 ng/L pyriproxyfen (~25% effect concentration) for male neonate production, as expected, the Doublesex1 gene was up-regulated. In descending order, gene expression patterns were most discernable with respect to pyriproxyfen exposure status, season of stormwater sample collection, and wetland quality, as indicated by the index of biological integrity. However, the biological implications of the affected genes were not broadly clear given limited genome resources for invertebrates. Our study provides support for the utility of short-term whole-organism transcriptomic testing in D. magna to discern sample type, but highlights the need for further work on invertebrate genomics. Environ Toxicol Chem 2022;41:2107-2123. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Identification and Evolution Analysis of the Complete Methyl Farnesoate Biosynthesis and Related Pathway Genes in the Mud Crab, Scylla paramamosain

The sesquiterpenoid hormone methyl farnesoate (MF) plays a vital role during crustacean development, which is mainly evidenced by its varied titers during different developmental stages. However, the biosynthesis pathways of MF remain obscure to some extent. In this study, we identified the complete MF biosynthesis and related pathway genes in Scylla paramamosain, including three involved in acetyl-CoA metabolism, eight in the mevalonate pathway, five in the sesquiterpenoids synthesis pathway, and five in the methionine cycle pathway. Bioinformatics, genomic structure, and phylogenetic analysis indicated that the JH biosynthesis genes might have experienced evolution after species differentiation. The mRNA tissue distribution analysis revealed that almost all genes involving in or relating to MF syntheses were highly expressed in the mandibular organ (MO), among which juvenile hormone acid methyltransferase was exclusively expressed in the MO, suggesting that most of these genes might mainly function in MF biosynthesis and that the methionine cycle pathway genes might play a crucial regulatory role during MF synthesis. In addition, the phylogenetic and tissue distribution analysis of the cytochrome P450 CYP15-like gene suggested that the epoxidized JHs might exist in crustaceans, but are mainly synthesized in hepatopancreas rather than the MO. Finally, we also found that betaine-homocysteine S-methyltransferase genes were lost in insects while methionine synthase was probably lost in most insects except Folsomia candida, indicating a regulatory discrepancy in the methionine cycle between crustaceans and insects. This study might increase our understanding of synthetic metabolism tailored for sesquiterpenoid hormones in S. paramamosain and other closely related species.

Transcriptome analysis of different life-history stages and screening of male-biased genes in Daphnia sinensis

Background

In the life history of Daphnia, the reproductive mode of parthenogenesis and sexual reproduction alternate in aquatic ecosystem, which are often affected by environmental and genetic factors. Recently, the sex-biased genes are of great significance for clarifying the origin and evolution of reproductive transformation and the molecular regulation mechanism of sex determination in Daphnia. Although some genes on reproductive transition of Daphnia had been researched, molecular mechanism on the maintenance of sexually dimorphic phenotypes of Daphnia are still not well known, including differentially expressed genes in different life-history stages.

Results

In this study, four life-history stages of Daphnia sinensis, juvenile female (JF), parthenogenetic female (PF), sexual female (SF) and male (M), were performed for transcriptome, and male-biased genes were screened. A total of 110437 transcripts were obtained and assembled into 22996 unigenes. In the four life-history stages (JF, PF, SF and M), the number of unique unigenes is respectively 2863, 445, 437 and 586, and the number of common unigenes is 9708. The differentially expressed genes (DEGs) between male and other three female stages (M vs JF, M vs PF and M vs SF) were 4570, 4358 and 2855, respectively. GO gene enrichment analysis showed that the up-regulated genes in male were mainly enriched in hydrolase activity and peptidase activity. Thirty-six genes in male were significantly higher expression than in the three female stages, including one Doublesex (Dsx) gene, one laminin gene, five trypsin genes and one serine protease genes, and one chitin synthase gene and two chitinase genes.

Conclusions

Our results showed that thirty-six candidate genes may be as the male-biased genes involving in the maintenance of sexually dimorphic phenotypes. This work will provide a reference for further exploring the functional genes related to sex differentiation in Daphnia species. Moreover, according to previous investigations, we thought that the expression level of functional genes may be related to the life-history stages of organisms, and may be also affected by different Daphnia species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08824-x.

Daphnia as a versatile model system in ecology and evolution

Water fleas of the genus Daphnia have been a model system for hundreds of years and is among the best studied ecological model organisms to date. Daphnia are planktonic crustaceans with a cyclic parthenogenetic life-cycle. They have a nearly worldwide distribution, inhabiting standing fresh- and brackish water bodies, from small temporary pools to large lakes. Their predominantly asexual reproduction allows for the study of phenotypes excluding genetic variation, enabling us to separate genetic from non-genetic effects. Daphnia are often used in studies related to ecotoxicology, predator-induced defence, host–parasite interactions, phenotypic plasticity and, increasingly, in evolutionary genomics. The most commonly studied species are Daphnia magna and D. pulex, for which a rapidly increasing number of genetic and genomic tools are available. Here, I review current research topics, where the Daphnia model system plays a critical role.

The Gluopsins: Opsins without the Retinal Binding Lysine

Opsins allow us to see. They are G-protein-coupled receptors and bind as ligand retinal, which is bound covalently to a lysine in the seventh transmembrane domain. This makes opsins light-sensitive. The lysine is so conserved that it is used to define a sequence as an opsin and thus phylogenetic opsin reconstructions discard any sequence without it. However, recently, opsins were found that function not only as photoreceptors but also as chemoreceptors. For chemoreception, the lysine is not needed. Therefore, we wondered: Do opsins exists that have lost this lysine during evolution? To find such opsins, we built an automatic pipeline for reconstructing a large-scale opsin phylogeny. The pipeline compiles and aligns sequences from public sources, reconstructs the phylogeny, prunes rogue sequences, and visualizes the resulting tree. Our final opsin phylogeny is the largest to date with 4956 opsins. Among them is a clade of 33 opsins that have the lysine replaced by glutamic acid. Thus, we call them gluopsins. The gluopsins are mainly dragonfly and butterfly opsins, closely related to the RGR-opsins and the retinochromes. Like those, they have a derived NPxxY motif. However, what their particular function is, remains to be seen.

Nanoplastic Toxicity: Insights and Challenges from Experimental Model Systems

Nanoplastic particles (NPs) can be produced or derived from the degradation of several daily used products and can therefore be found in the air, water, and food. Every day, these microscopic particles are confronted by different routes of exposure. Recent investigations have shown the internalization of these particles, differing in size and modification, in vivo in aquatic organisms and terrestrial organisms, as well as in vitro in different human cell lines. During the last years, the number of studies investigating the effects of NPs using widely different model systems and experimental approaches is exponentially growing, thus providing information about NPs, especially about polystyrene particle toxicity on health. To facilitate the grasping of the most relevant information, an overview is provided on the toxic effects of NPs coming from studies in cellular systems and in vivo in model organisms and on aspects which can be of particular relevance for particle toxicity (e.g., particle internalization mechanisms and structural modifications). Major achievements and gaps in the field as well as the point of view on how more systematic studies and exploitation of in vivo model organisms may improve the knowledge on important aspects of NPs are also pointed out.

Photosynthetic activity in both algae and cyanobacteria changes in response to cues of predation

A plethora of adaptive responses to predation has been described in microscopic aquatic producers. Although the energetic costs of these responses are expected, with their consequences going far beyond an individual, their underlying molecular and metabolic mechanisms are not fully known. One, so far hardly considered, is if and how the photosynthetic efficiency of phytoplankton might change in response to the predation cues. Our main aim was to identify such responses in phytoplankton and to detect if they are taxon-specific. We exposed seven algae and seven cyanobacteria species to the chemical cues of an efficient consumer, Daphnia magna, which was fed either a green alga, Acutodesmus obliquus, or a cyanobacterium, Synechococcus elongatus (kairomone and alarm cues), or was not fed (kairomone alone). In most algal and cyanobacterial species studied, the quantum yield of photosystem II increased in response to predator fed cyanobacterium, whereas in most of these species the yield did not change in response to predator fed alga. Also, cyanobacteria tended not to respond to a non-feeding predator. The modal qualitative responses of the electron transport rate were similar to those of the quantum yield. To our best knowledge, the results presented here are the broadest scan of photosystem II responses in the predation context so far.

Phenotypic and transcriptional response of Daphnia pulicaria to the combined effects of temperature and predation

Daphnia, an ecologically important zooplankton species in lakes, shows both genetic adaptation and phenotypic plasticity in response to temperature and fish predation, but little is known about the molecular basis of these responses and their potential interactions. We performed a factorial experiment exposing laboratory-propagated Daphnia pulicaria clones from two lakes in the Sierra Nevada mountains of California to normal or high temperature (15°C or 25°C) in the presence or absence of fish kairomones, then measured changes in life history and gene expression. Exposure to kairomones increased upper thermal tolerance limits for physiological activity in both clones. Cloned individuals matured at a younger age in response to higher temperature and kairomones, while size at maturity, fecundity and population intrinsic growth were only affected by temperature. At the molecular level, both clones expressed more genes differently in response to temperature than predation, but specific genes involved in metabolic, cellular, and genetic processes responded differently between the two clones. Although gene expression differed more between clones from different lakes than experimental treatments, similar phenotypic responses to predation risk and warming arose from these clone-specific patterns. Our results suggest that phenotypic plasticity responses to temperature and kairomones interact synergistically, with exposure to fish predators increasing the tolerance of Daphnia pulicaria to stressful temperatures, and that similar phenotypic responses to temperature and predator cues can be produced by divergent patterns of gene regulation.

A thorough annotation of the krill transcriptome offers new insights for the study of physiological processes

The krill species Euphausia superba plays a critical role in the food chain of the Antarctic ecosystem. Significant changes in climate conditions observed in the Antarctic Peninsula region in the last decades have already altered the distribution of krill and its reproductive dynamics. A deeper understanding of the adaptation capabilities of this species is urgently needed. The availability of a large body of RNA-seq assays allowed us to extend the current knowledge of the krill transcriptome. Our study covered the entire developmental process providing information of central relevance for ecological studies. Here we identified a series of genes involved in different steps of the krill moulting cycle, in the reproductive process and in sexual maturation in accordance with what was already described in previous works. Furthermore, the new transcriptome highlighted the presence of differentially expressed genes previously unknown, playing important roles in cuticle development as well as in energy storage during the krill life cycle. The discovery of new opsin sequences, specifically rhabdomeric opsins, one onychopsin, and one non-visual arthropsin, expands our knowledge of the krill opsin repertoire. We have collected all these results into the KrillDB2 database, a resource combining the latest annotation of the krill transcriptome with a series of analyses targeting genes relevant to krill physiology. KrillDB2 provides in a single resource a comprehensive catalog of krill genes; an atlas of their expression profiles over all RNA-seq datasets publicly available; a study of differential expression across multiple conditions. Finally, it provides initial indications about the expression of microRNA precursors, whose contribution to krill physiology has never been reported before.

The reproductive effects of the cancer chemotherapy agent, Carmofur, on Daphnia magna are mediated by its metabolite, 5-Fluorouracil

Carmofur is an antineoplastic agent that inhibits ceramidase, a key enzyme in the sphingolipid pathway. Previous research suggests carmofur represses reproductive maturity in Daphnia magna. The purpose of this experiment was to confirm carmofur's effects on fecundity and reproductive maturity over two generations. A chronic toxicity test found reproductive maturity was delayed from 9 to 19 days by 0.80 μM carmofur with a 99.7% drop in reproduction, probably caused by delayed ovarian development. Second generation effects were even greater with 0% reproductive success at 0.40 μM. To our surprise, carmofur was not measured in the media by HPLC 24 h after exposure. Previous research indicated that carmofur is unstable in water and hydrolyzed into 5-fluorouracil (5-FU). Therefore, the chronic toxicity study was repeated with 5-FU and similar effects on reproductive maturity were observed at similar concentrations despite very different acute toxicities (48 h carmofur LC50 = 1.93 μM; 5-FU LC50 = 207 μM). 5-FU delayed reproductive maturity from 9 to 21 days with a 71.12% drop in reproduction at 0.80 μM and greater effects in the 2^nd generation similar to carmofur. 5-FU was found stable in aquatic media and HPLC confirmed 5-FU was hydrolyzed from carmofur within 24 h. In conclusion, carmofur and 5-FU reduce fecundity because they delay reproductive maturity and ovarian development in Daphnia magna. We conclude that the reproductive effects observed after carmofur treatment are primarily mediated by its breakdown product, 5-FU. This further underscores the importance of measuring chemical concentrations and evaluating chemical metabolism and decomposition when determining toxicity, especially of chemotherapeutic agents.Clinical trials registration Not applicable.

The genome of a mangrove plant, Avicennia marina, provides insights into adaptation to coastal intertidal habitats

Whole-genome duplication, gene family and lineage-specific genes analysis based on high-quality genome reveal the adaptation mechanisms of Avicennia marina to coastal intertidal habitats. Mangrove plants grow in a complex habitat of coastal intertidal zones with high salinity, hypoxia, etc. Therefore, it is an interesting question how mangroves adapt to the unique intertidal environment. Here, we present a chromosome-level genome of the Avicennia marina, a typical true mangrove with a size of 480.43 Mb, contig N50 of 11.33 Mb and 30,956 annotated protein-coding genes. We identified 621 Avicennia-specific genes that are mainly related to flavonoid and lignin biosynthesis, auxin homeostasis and response to abiotic stimulus. We found that A. marina underwent a novel specific whole-genome duplication, which is in line with a brief era of global warming that occurred during the paleocene-eocene maximum. Comparative genomic and transcriptomic analyses outline the distinct evolution and sophisticated regulations of A. marina adaptation to the intertidal environments, including expansion of photosynthesis and oxidative phosphorylation gene families, unique genes and pathways for antibacterial, detoxifying antioxidant and reactive oxygen species scavenging. In addition, we also analyzed salt gland secretion-related genes, and those involved in the red bark-related flavonoid biosynthesis, while significant expansions of key genes such as NHX, 4CL, CHS and CHI. High-quality genomes in future investigations will facilitate the understand of evolution of mangrove and improve breeding.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Genome-wide identification of 194 G protein-coupled receptor (GPCR) genes from the water flea Daphnia magna

In crustaceans, G protein-coupled receptors (GPCRs) are the largest transmembrane receptor family and function by mediating various environmental stimuli in cells. Understanding GPCR signaling is crucial to better understanding of crustacean endocrinology. GPCRs evolved from early eukaryotes, and genome-wide identification of GPCRs in a particular taxon can provide insight into evolutionary tendencies and adaptive strategies of GPCR response to environmental stimuli. Here, we identified 194 full-length GPCR genes in the water flea Daphnia magna that can be divided into five distinct classes (A, B, C, F, and other). A strong orthologous relationship for amine, neuropeptide, and opsin receptors was found in the phylogenetic comparison of D. magna GPCRs to those of humans and two well-known insects (Drosophila melanogaster and Solenopsis invicta). Our results based on phylogenetic relationships suggest that most GPCRs subfamilies have undergone sporadic evolutionary processes for adaptation to environmental pressures. Despite the dynamics of GPCR evolution, some GPCRs are highly conserved between species. This study provides a better understanding of the evolution of GPCRs and expands our knowledge of the potential physiological mechanisms in D. magna in response to various environmental stimuli.

Genome of a giant isopod, Bathynomus jamesi, provides insights into body size evolution and adaptation to deep-sea environment

BACKGROUND

The deep-sea may be regarded as a hostile living environment, due to low temperature, high hydrostatic pressure, and limited food and light. Isopods, a species-rich group of crustaceans, are widely distributed across different environments including the deep sea and as such are a useful model for studying adaptation, migration, and speciation. Similar to other deep-sea organisms, giant isopods have larger body size than their shallow water relatives and have large stomachs and fat bodies presumably to store organic reserves. In order to shed light on the genetic basis of these large crustaceans adapting to the oligotrophic environment of deep-sea, the high-quality genome of a deep-sea giant isopod Bathynomus jamesi was sequenced and assembled.

RESULTS

B. jamesi has a large genome of 5.89 Gb, representing the largest sequenced crustacean genome to date. Its large genome size is mainly attributable to the remarkable proliferation of transposable elements (84%), which may enable high genome plasticity for adaptive evolution. Unlike its relatives with small body size, B. jamesi has expanded gene families related to pathways of thyroid and insulin hormone signaling that potentially contribute to its large body size. Transcriptomic analysis showed that some expanded gene families related to glycolysis and vesicular transport were specifically expressed in its digestive organs. In addition, comparative genomics and gene expression analyses in six tissues suggested that B. jamesi has inefficient lipid degradation, low basal metabolic rate, and bulk food storage, suggesting giant isopods adopt a more efficient mechanism of nutrient absorption, storage, and utilization to provide sustained energy supply for their large body size.

CONCLUSIONS

Taken together, the giant isopod genome may provide a valuable resource for understanding body size evolution and adaptation mechanisms of macrobenthic organisms to deep-sea environments.

A deep learning approach for orphan gene identification in moso bamboo (Phyllostachys edulis) based on the CNN + Transformer model

Background

Orphan gene play an important role in the environmental stresses of many species and their identification is a critical step to understand biological functions. Moso bamboo has high ecological, economic and cultural value. Studies have shown that the growth of moso bamboo is influenced by various stresses. Several traditional methods are time-consuming and inefficient. Hence, the development of efficient and high-accuracy computational methods for predicting orphan genes is of great significance.

Results

In this paper, we propose a novel deep learning model (CNN + Transformer) for identifying orphan genes in moso bamboo. It uses a convolutional neural network in combination with a transformer neural network to capture k-mer amino acids and features between k-mer amino acids in protein sequences. The experimental results show that the average balance accuracy value of CNN + Transformer on moso bamboo dataset can reach 0.875, and the average Matthews Correlation Coefficient (MCC) value can reach 0.471. For the same testing set, the Balance Accuracy (BA), Geometric Mean (GM), Bookmaker Informedness (BM), and MCC values of the recurrent neural network, long short-term memory, gated recurrent unit, and transformer models are all lower than those of CNN + Transformer, which indicated that the model has the extensive ability for OG identification in moso bamboo.

Conclusions

CNN + Transformer model is feasible and obtains the credible predictive results. It may also provide valuable references for other related research. As our knowledge, this is the first model to adopt the deep learning techniques for identifying orphan genes in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04702-1.

Invasions of an obligate asexual daphnid species support the nearly neutral theory

To verify the "nearly neutral theory (NNT)," the ratio of nonsynonymous to synonymous substitutions (dN/dS) was compared among populations of different species. To determine the validity of NNT, however, populations that are genetically isolated from each other but share the same selection agents and differ in size should be compared. Genetically different lineages of obligate asexual Daphnia pulex invading Japan from North America are an ideal example as they satisfy these prerequisites. Therefore, we analyzed the whole-genome sequences of 18 genotypes, including those of the two independently invaded D. pulex lineages (JPN1 and JPN2) and compared the dN/dS ratio between the lineages. The base substitution rate of each genotype demonstrated that the JPN1 lineage having a larger distribution range diverged earlier and thus was older than the JPN2 lineage. Comparisons of the genotypes within lineages revealed that changes in dN/dS occurred after the divergence and were larger in the younger lineage, JPN2. These results imply that the JPN1 lineage has been more effectively subjected to purification selections, while slightly deteriorating mutations are less purged in JPN2 with smaller population size. Altogether, the lineage-specific difference in the dN/dS ratio for the obligate asexual D. pulex was well explained by the NNT.

Transcriptional responses of Daphnia magna exposed to Akaki river water

Pollution of the aquatic environment is a global problem, with industrial waste, farming effluents, sewage, and wastewater as the main contributors. Many pollutants are biologically active at low concentrations, resulting in sublethal effects, which makes it a highly complex situation and difficult to assess. In many places, such as the Akaki river in Ethiopia, the pollution situation has resulted in streams with minimal presence of invertebrates or vertebrates. As it is difficult to perform a complete chemical analysis of the waters, the present study focused on using gene expression analysis as a biological end point to determine the effects of Akaki river contaminants. The present study was conducted using the small planktonic crustacean Daphnia magna with toxicogenomic molecular markers. Daphnia magna neonates were exposed to Akaki water samples collected from two different sites on the river and analyzed for mortality and expression of genes involved in different biological pathways. Despite the poor quality of Akaki river water, 48 h acute toxicity tests showed no mortality. Interestingly, analysis of sublethal toxicogenomic responses showed that exposure to Akaki water altered the expression of 25 out of 37 genes involved in metal regulation, immune response, oxidative stress, respiration, reproduction, and development. The toxicogenomic data gives insight into the mechanisms involved in causing potential adverse effects to aquatic biota harboring the Akaki river system.

Historical exposure to chemicals reduces tolerance to novel chemical stress in Daphnia (waterflea)

Until the last few decades, anthropogenic chemicals used in most production processes have not been comprehensively assessed for their risk and impact on wildlife and humans. They are transported globally and usually end up in the environment as unintentional pollutants, causing long‐term adverse effects. Modern toxicology practices typically use acute toxicity tests of unrealistic concentrations of chemicals to determine their safe use, missing pathological effects arising from long‐term exposures to environmentally relevant concentrations. Here, we study the transgenerational effect of environmentally relevant concentrations of five chemicals on the priority list of international regulatory frameworks on the keystone species Daphnia magna. We expose Daphnia genotypes resurrected from the sedimentary archive of a lake with a known history of chemical pollution to the five chemicals to understand how historical exposure to chemicals influences adaptive responses to novel chemical stress. We measure within‐ and transgenerational plasticity in fitness‐linked life history traits following exposure of “experienced” and “naive” genotypes to novel chemical stress. As the revived Daphnia originate from the same genetic pool sampled at different times in the past, we are able to quantify the long‐term evolutionary impact of chemical pollution by studying genome‐wide diversity and identifying functional pathways affected by historical chemical stress. Our results suggest that historical exposure to chemical stress causes reduced genome‐wide diversity, leading to lower cross‐generational tolerance to novel chemical stress. Lower tolerance is underpinned by reduced gene diversity at detoxification, catabolism and endocrine genes in experienced genotypes. We show that these genes sit within pathways that are conserved and potential chemical targets in other species, including humans.

Chromosome-level genome of the globe skimmer dragonfly (Pantala flavescens)

BACKGROUND

The globe skimmer dragonfly (Pantala flavescens) is a notable Odonata insect distributed in nature fields and farmlands worldwide, and it is commonly recognized as a natural enemy because it preys on agricultural pests and health pests. As one of the sister groups of winged insects, odonatan species are key to understanding the evolution of insect wings.

FINDINGS

We present a high-quality reference genome of P. flavescens, which is the first chromosome-level genome in the Palaeoptera (Odonata and Ephemeroptera). The assembled genome size was 662 Mb, with a contig N50 of 16.2 Mb. Via Hi-C scaffolding, 648 Mb (97.9%) of contig sequences were clustered, ordered, and assembled into 12 large scaffolds, each corresponding to a natural chromosome. The X chromosome was identified by sequence coverage depth. The repetitive sequences and gene density of the X chromosome are similar to those of autosomal sequences, but the X chromosome shows a much lower degree of heterozygosity. Our analysis shows that the effective population size experienced 3 declining events, which may have been caused by climate change and environmental pollution.

CONCLUSIONS

The genome of P. flavescens provides more information on the biology and evolution of insects and will help for the use of this species in pest control.

Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.

A wide-field micro-computed tomography detector: micron resolution at half-centimetre scale

Ideal three-dimensional imaging of complex samples made up of micron-scale structures extending over mm to cm, such as biological tissues, requires both wide field of view and high resolution. For existing optics and detectors used for micro-CT (computed tomography) imaging, sub-micron pixel resolution can only be achieved for fields of view of <2 mm. This article presents a unique detector system with a 6 mm field-of-view image circle and 0.5 µm pixel size that can be used in micro-CT units utilizing both synchrotron and commercial X-ray sources. A resolution-test pattern with linear microstructures and whole adult Daphnia magna were imaged at beamline 8.3.2 of the Berkeley Advanced Light Source. Volumes of 10000 × 10000 × 7096 isotropic 0.5 µm voxels were reconstructed over a 5.0 mm × 3.5 mm field of view. Measurements in the projection domain confirmed a 0.90 µm measured spatial resolution that is largely Nyquist-limited. This unprecedented combination of field of view and resolution dramatically reduces the need for sectional scans and computational stitching for large samples, ultimately offering the means to elucidate changes in tissue and cellular morphology in the context of larger, whole, intact model organisms and specimens. This system is also anticipated to benefit micro-CT imaging in materials science, microelectronics, agricultural science and biomedical engineering.

Emerging trends in nanoparticle toxicity and the significance of using Daphnia as a model organism

Nanoparticle production is on the rise due to its many uses in the burgeoning nanotechnology industry. Although nanoparticles have growing applications, there is great concern over their environmental impact due to their inevitable release into the environment. With uncertainty of environmental concentration and risk to aquatic organisms, the microcrustacean Daphnia spp. has emerged as an important freshwater model organism for risk assessment of nanoparticles because of its biological properties, including parthenogenetic reproduction; small size and short generation time; wide range of endpoints for ecotoxicological studies; known genome, useful for providing mechanistic information; and high sensitivity to environmental contaminants and other stressors. In this review, we (1) highlight the advantages of using Daphnia as an experimental model organism for nanotoxicity studies, (2) summarize the impacts of nanoparticle physicochemical characteristics on toxicity in relation to Daphnia, and (3) summarize the effects of nanoparticles (including nanoplastics) on Daphnia as well as mechanisms of toxicity, and (4) highlight research uncertainties and recommend future directions necessary to develop a deeper understanding of the fate and toxicity of nanoparticles and for the development of safer and more sustainable nanotechnology.

iLoci: robust evaluation of genome content and organization for provisional and mature genome assemblies

We introduce a new framework for genome analyses based on parsing an annotated genome assembly into distinct interval loci (iLoci), available as open-source software as part of the AEGeAn Toolkit (https://github.com/BrendelGroup/AEGeAn). We demonstrate that iLoci provide an alternative coordinate system that is robust to changes in assembly and annotation versions and facilitates granular quality control of genome data. We discuss how statistics computed on iLoci reflect various characteristics of genome content and organization and illustrate how these statistics can be used to establish a baseline for assessment of the completeness and accuracy of the data. We also introduce a well-defined measure of relative genome compactness and compute other iLocus statistics that reveal genome-wide characteristics of gene arrangements in the whole genome context. Given the fast pace of assembly/annotation updates, our AEGeAn Toolkit fills a niche in computational genomics based on deriving persistent and species-specific genome statistics. Gene structure model-centric iLoci provide a precisely defined coordinate system that can be used to store assembly/annotation updates that reflect either stable or changed assessments. Large-scale application of the approach revealed species- and clade-specific genome organization in precisely defined computational terms, promising intriguing forays into the forces of shaping genome structure as more and more genome assemblies are being deposited.

Development of transgenic Daphnia magna for visualizing homology-directed repair of DNA

In the crustacean Daphnia magna, studying homology-directed repair (HDR) is important to understand genome maintenance during parthenogenesis, effects of environmental toxicants on the genome, and improvement of HDR-mediated genome editing. Here we developed a transgenic D. magna that expresses green fluorescence protein (GFP) upon HDR occurrence. We utilized the previously established reporter plasmid named DR-GFP that has a mutated eGFP gene (SceGFP) and the tandemly located donor GFP gene fragment (iGFP). Upon double-strand break (DSB) introduction on SceGFP, the iGFP gene fragment acts as the HDR template and restores functional eGFP expression. We customized this reporter plasmid to allow bicistronic expression of the mCherry gene under the control of the D. magna EF1α-1 promoter/enhancer. By CRISPR/Cas-mediated knock-in of this plasmid via non-homologous joining, we generated the transgenic D. magna that expresses mCherry ubiquitously, suggesting that the DR-GFP reporter gene is expressed in most cells. Introducing DSB on the SceGFP resulted in eGFP expression and this HDR event could be detected by fluorescence, genomic PCR, and quantitative reverse-transcription PCR, suggesting this line could be used for evaluating HDR. The established reporter line might expand our understanding of the HDR mechanism and also improve the HDR-based gene-editing system in this species.

Evaluating Multiple Stressor Effects on Benthic–Pelagic Freshwater Communities in Systems of Different Complexities: Challenges in Upscaling

Upscaling of ecological effects from indoor microcosms to outdoor mesocosms bridging the gap between controlled laboratory conditions and highly complex natural environments poses several challenges: typical standard water types used in laboratory experiments are not feasible in large outdoor experiments. Additionally, moving from the micro- to meso-scale, biodiversity is enhanced. We performed an indoor microcosm experiment to determine the effects of agricultural run-off (ARO) on a defined benthic–pelagic community comprising primary producers and primary consumers, exposed to ambient summer temperature and +3.5 °C. Treatments were replicated in two water types (standard Volvic and Munich well water). We then scaled up to outdoor mesocosms using an ARO concentration gradient and +3 °C warming above ambient temperature, using Munich well water. We included the same benthic macroorganisms but more complex periphyton and plankton communities. All the functional groups were affected by stressors in the microcosms, and a shift from macrophyte to phytoplankton dominance was observed. While effects were present, they were less pronounced in the mesocosms, where a higher biodiversity may have modified the responses of the system to the stressors. The stressor effects observed in controlled experiments may thus be masked in more complex outdoor experiments, but should not be interpreted as “no effects”.

Rethinking opsins

Opsins, the protein moieties of animal visual photo-pigments, have emerged as moonlighting proteins with diverse, light-dependent and -independent physiological functions. This raises the need to revise some basic assumptions concerning opsin expression, structure, classification, and evolution.

Metabolic Phenotype of Daphnia Under Hypoxia: Macroevolution, Microevolution, and Phenotypic Plasticity

Daphnia is a freshwater crustacean that is able to upregulate hemoglobin (Hb) in response to hypoxia, imparting a red color. We combine multiple field surveys across season with a lab experiment to evaluate changes in the metabolic phenotype of Daphnia in relation to environmental hypoxia. Looking at the zooplankton community, we found that D. pulicaria was restricted to lakes with a hypoxic hypolimnion. Comparing D. pulicaria with different amounts of Hb, red animals showed higher mRNA levels for several Hb genes, whereas most glycolytic genes showed red/pale differences of less than 50%. We also observed seasonal changes in the metabolic phenotype that differed between red and pale animals. Hb was upregulated early in the season in hypoxic lakes, and a relationship between Hb and lactate dehydrogenase only emerged later in the season in a temporal pattern that was lake specific. To evaluate whether these differences were due to specific lake environments or microevolutionary differences, we tested the induction of genes under controlled hypoxia in isofemale lines from each of four lakes. We found a strong response to 18 h hypoxia exposure in both Hb and lactate dehydrogenase mRNA, although the magnitude of the acute response was greater than the steady state differences in mRNA levels between pale and red Daphnia. The baseline expression of Hb and lactate dehydrogenase also varied between isofemale lines with different lake origins. These results, in combination with comparison of glycogen measurements, suggests that Hb functions primarily to facilitate oxygen delivery, mitigating systemic hypoxia, rather than an oxygen store. The combination of lab and field studies suggest that the metabolic phenotype of the animal is influenced by both microevolutionary differences (within and between lakes) as well as the spatial and temporal environmental heterogeneity of the lakes. The differences between Daphnia species, and the unexpected lack of hypoxia sensitivity of select glycolytic genes provide evidence of macroevolutionary differences in metabolic strategies to cope with hypoxia.

Juvenile hormone synthesis and signaling disruption triggering male offspring induction and population decline in cladocerans (water flea): Review and adverse outcome pathway development

Juvenile hormone (JH) are a family of multifunctional hormones regulating larval development, molting, metamorphosis, reproduction, and phenotypic plasticity in arthropods. Based on its importance in arthropod life histories, many insect growth regulators (IGRs) mimicking JH have been designed to control harmful insects in agriculture and aquaculture. These JH analogs (JHAs) may also pose hazards to nontarget species by causing unexpected endocrine-disrupting (ED) effects such as molting and metamorphosis defects, larval lethality, and disruption of the sexual identity. This critical review summarizes the current knowledge of the JH-mediated effects in the freshwater cladoceran crustaceans such as Daphnia species on JHA-triggered endocrine disruptive outputs to establish a systematic understanding of JHA effects. Based on the current knowledge, adverse outcome pathways (AOPs) addressing the JHA-mediated ED effects in cladoceran leading to male offspring production and subsequent population decline were developed. The weight of evidence (WoE) of AOPs was assessed according to established guidelines. The review and AOP development aim to present the current scientific understanding of the JH pathway and provide a robust reference for the development of tiered testing strategies and new risk assessment approaches for JHAs in future ecotoxicological research and regulatory processes.

Studying phenotypic variation and DNA methylation across development, ecology and evolution in the clonal marbled crayfish: a paradigm for investigating epigenotype-phenotype relationships in macro-invertebrates

Animals can produce different phenotypes from the same genome during development, environmental adaptation and evolution, which is mediated by epigenetic mechanisms including DNA methylation. The obligatory parthenogenetic marbled crayfish, Procambarus virginalis, whose genome and methylome are fully established, proved very suitable to study this issue in detail. Comparison between developmental stages and DNA methylation revealed low expression of Dnmt methylation and Tet demethylation enzymes from the spawned oocyte to the 256 cell embryo and considerably increased expression thereafter. The global 5-methylcytosine level was 2.78% at mid-embryonic development and decreased slightly to 2.41% in 2-year-old adults. Genetically identical clutch-mates raised in the same uniform laboratory setting showed broad variation in morphological, behavioural and life history traits and differences in DNA methylation. The invasion of diverse habitats in tropical to cold-temperate biomes in the last 20 years by the marbled crayfish was associated with the expression of significantly different phenotypic traits and DNA methylation patterns, despite extremely low genetic variation on the whole genome scale, suggesting the establishment of epigenetic ecotypes. The evolution of marbled crayfish from its parent species Procambarus fallax by autotriploidy a few decades ago was accompanied by a significant increase in body size, fertility and life span, a 20% reduction of global DNA methylation and alteration of methylation in hundreds of genes, suggesting that epigenetic mechanisms were involved in speciation and fitness enhancement. The combined analysis of phenotypic traits and DNA methylation across multiple biological contexts in the laboratory and field in marbled crayfish may serve as a blueprint for uncovering the role of epigenetic mechanisms in shaping of phenotypes in macro-invertebrates.

The freshwater water flea Daphnia magna NIES strain genome as a resource for CRISPR/Cas9 gene targeting: The glutathione S-transferase omega 2 gene

The water flea Daphnia magna is a small freshwater planktonic animal in the Cladocera. In this study, we assembled the genome of the D. magna NIES strain, which is widely used for gene targeting but has no reported genome. We used the long-read sequenced data of the Oxford nanopore sequencing tool for assembly. Using 3,231 genetic markers, the draft genome of the D. magna NIES strain was built into ten linkage groups (LGs) with 483 unanchored contigs, comprising a genome size of 173.47 Mb. The N50 value of the genome was 12.54 Mb and the benchmarking universal single-copy ortholog value was 98.8%. Repeat elements in the D. magna NIES genome were 40.8%, which was larger than other Daphnia spp. In the D. magna NIES genome, 15,684 genes were functionally annotated. To assess the genome of the D. magna NIES strain for CRISPR/Cas9 gene targeting, we selected glutathione S-transferase omega 2 (GST-O2), which is an important gene for the biotransformation of arsenic in aquatic organisms, and targeted it with an efficient make-up (25.0%) of mutant lines. In addition, we measured reactive oxygen species and antioxidant enzymatic activity between wild type and a mutant of the GST-O2 targeted D. magna NIES strain in response to arsenic. In this study, we present the genome of the D. magna NIES strain using GST-O2 as an example of gene targeting, which will contribute to the construction of deletion mutants by CRISPR/Cas9 technology.

Transcriptomic analysis of brine shrimp Artemia franciscana across a wide range of salinities

Brine shrimp Artemia franciscana, a commercially important species, can thrive in a wide range of salinities and is commonly found in hypersaline lakes and solar salterns. Transcriptome analysis can enhance the understanding of the adaptative mechanisms of brine shrimp in aquaculture. RNA sequencing (RNAseq) data was generated from A. franciscana adults that were salt-adapted for 2-4 weeks at five salinities: 35, 50, 100, 150, and 230 psu. Long-read isoform sequencing (IsoSeq) data was used to construct a high-quality transcriptome assembly. Also, the gene expression patterns in A. franciscana adults were examined. Notably, the transcriptional response of A. franciscana's acclimation to intermediate salinities (50-150 psu) displayed frequently and differentially U-shaped or inverted U-shaped expression patterns. In addition, the types of genes showing two nonmonotonic expression patterns were distinct from each other. The coordinated shifts in gene expression suggest different homeostatic strategies of A. franciscana at specific salinities; such strategies may enhance population fitness at extreme salinities. Our study should promote a scientific concept for the gene expression patterns of A. franciscana along a broad salinity gradient, and a variety of salinity and prey should be monitored for testing the gene expression pattern of this important aquaculture species.

Foster thy young: enhanced prediction of orphan genes in assembled genomes

Proteins encoded by newly-emerged genes ('orphan genes') share no sequence similarity with proteins in any other species. They provide organisms with a reservoir of genetic elements to quickly respond to changing selection pressures. Here, we systematically assess the ability of five gene prediction pipelines to accurately predict genes in genomes according to phylostratal origin. BRAKER and MAKER are existing, popular ab initio tools that infer gene structures by machine learning. Direct Inference is an evidence-based pipeline we developed to predict gene structures from alignments of RNA-Seq data. The BIND pipeline integrates ab initio predictions of BRAKER and Direct inference; MIND combines Direct Inference and MAKER predictions. We use highly-curated Arabidopsis and yeast annotations as gold-standard benchmarks, and cross-validate in rice. Each pipeline under-predicts orphan genes (as few as 11 percent, under one prediction scenario). Increasing RNA-Seq diversity greatly improves prediction efficacy. The combined methods (BIND and MIND) yield best predictions overall, BIND identifying 68% of annotated orphan genes, 99% of ancient genes, and give the highest sensitivity score regardless dataset in Arabidopsis. We provide a light weight, flexible, reproducible, and well-documented solution to improve gene prediction.

Major Evolutionary Transitions and the Roles of Facilitation and Information in Ecosystem Transformations

A small number of extraordinary “Major Evolutionary Transitions” (METs) have attracted attention among biologists. They comprise novel forms of individuality and information, and are defined in relation to organismal complexity, irrespective of broader ecosystem-level effects. This divorce between evolutionary and ecological consequences qualifies unicellular eukaryotes, for example, as a MET although they alone failed to significantly alter ecosystems. Additionally, this definition excludes revolutionary innovations not fitting into either MET type (e.g., photosynthesis). We recombine evolution with ecology to explore how and why entire ecosystems were newly created or radically altered – as Major System Transitions (MSTs). In doing so, we highlight important morphological adaptations that spread through populations because of their immediate, direct-fitness advantages for individuals. These are Major Competitive Transitions, or MCTs. We argue that often multiple METs and MCTs must be present to produce MSTs. For example, sexually-reproducing, multicellular eukaryotes (METs) with anisogamy and exoskeletons (MCTs) significantly altered ecosystems during the Cambrian. Therefore, we introduce the concepts of Facilitating Evolutionary Transitions (FETs) and Catalysts as key events or agents that are insufficient themselves to set a MST into motion, but are essential parts of synergies that do. We further elucidate the role of information in MSTs as transitions across five levels: (I) Encoded; (II) Epigenomic; (III) Learned; (IV) Inscribed; and (V) Dark Information. The latter is ‘authored’ by abiotic entities rather than biological organisms. Level IV has arguably allowed humans to produce a MST, and V perhaps makes us a FET for a future transition that melds biotic and abiotic life into one entity. Understanding the interactive processes involved in past major transitions will illuminate both current events and the surprising possibilities that abiotically-created information may produce.