Reproductive switching analysis of Daphnia similoides between sexual female and parthenogenetic female by transcriptome comparison

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

Pyroclastic Dust from Arequipa-Peru Decorated with Iron Oxide Nanoparticles and Their Ecotoxicological Properties in Water Flea D. magna

A novel magnetic composite made of Peruvian pyroclastic dust material decorated with maghemite nanoparticles was synthesized and characterized using a variety of analytic techniques. The 13 nm maghemite nanoparticles were grown on the pyroclastic dust using the conventional coprecipitation chemical route. A short-term acute assay was developed to study the ecotoxicological behavior of the water flea, Daphnia magna. A 24 h-lethal concentration (LC50) value equal to 123.6 mg L-1 was determined only for the magnetic composite. While the pyroclastic dust material did not exhibit a lethal concentration, it caused morphologically significant changes (p < 0.05) for heart and tail parameters at high concentrations. Morphologies exposed to the magnetic composite above the 24 h-LC50 revealed less tolerance and significant changes in the body, heart, antenna, and eye. Hence, it affects biomarker growth and swimming. The reproduction rate was not affected by the raw pyroclastic dust material. However, the number of individuals showed a decrease with increasing composite concentrations. The present study indicates the LC50 value, which can be used as a reference concentration for in-situ water cleaning with this material without damaging or changing the Daphnia magna ecosystem.

Ecotoxicological Properties of Pure and Phosphorus-Containing Graphene Oxide Bidimensional Sheets in Daphnia magna

In this work, the synthesis and structural, thermal, vibrational, morphological, and electronic characterization of 2D-like pure graphene oxide (GO) and phosphorus-containing graphene oxide (GOP) sheets were investigated. The average thicknesses of GO and GOP were 0.8 μm and 3.1 μm, respectively. The electron energy-loss spectroscopy spectra were used to analyze the differences in the C-K and O-K energy edge bands between GO and GOP. In addition, colloidal stability was studied using dynamic light scattering and zeta potential physicochemical techniques, determining that as the concentration increases, the hydrodynamic diameter and electrostatic stability of GO and GOP increase. The colloidal stability was quite important to ensure the interaction between the suspended solid phase and the biomarker. The 2D-like materials were used to determine their ecotoxicological properties, such as the medium lethal concentration, a crucial parameter for understanding ecotoxicity. Acute ecotoxicity experiments (24 h) were conducted in triplicate to obtain robust statistics, with corresponding mean lethal concentration (LC50) of 11.4 mg L-1 and 9.8 mg L-1 for GO and GOP, respectively. The morphological parameters of GO and GOP were compared with a negative control. However, only the case of GO was analyzed, since the Daphnia magna (D. magna) set exposed to GOP died before completing the time required for morphological analysis. The results indicate that the GOP sample is more toxic than the GO, both during and after exposure. Furthermore, the morphological parameters with the greatest statistically significant changes (p<0.05) were associated with the heart and body, while the eye and tail showed less significant changes.

Acute toxicity, oxidative stress, and apoptosis due to short-term triclosan exposure and multi- and transgenerational effects on in vivo endpoints, antioxidant defense, and DNA damage response in the freshwater water flea Daphnia magna

In this study, we measured the acute toxicity of triclosan (TCS) in neonate and adult Daphnia magna water fleas. The median lethal concentrations were 184.689 and 349.511 μg/L, respectively. Oxidative stress induced by TCS was analyzed based on changes in reactive oxygen species (ROS) content and antioxidant enzymatic activities in D. magna. Based on these endpoints, TCS concentrations of 50 and 100 μg/L induced oxidative stress. However, several apoptosis-mediated proteins showed TCS-induced oxidative-stress damage in response to 25 μg/L, indicating that apoptotic proteins were the most sensitive mediators. We also evaluated the multi- and transgenerational effects of TCS on D. magna over three generations in terms of various in vivo endpoints, DNA damage responses, and biochemical reactions. The transgenerational group exposed to TCS exhibited greater negative impacts on antioxidant responses, DNA fragmentation status, and biological endpoints compared with the multigenerational exposure group, leading to decreased reproductive rates and higher ROS content. The transcriptional expression levels of glutathione S-transferase genes in the transgenerational exposure group were upregulated compared to those in the multigenerational group but were fully recovered in F2 offspring. Our findings provide an in-depth understanding of the adaptive effects of multigenerational exposure to TCS.

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.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Effects of toxic Microcystis aeruginosa on the expression of Hox genes in Daphnia similoides sinensis

Lake eutrophication and cyanobacterial blooms have become worldwide environmental issues. Under cyanobacterial blooms (especially Microcystis), Daphnia spp. can transfer beneficial information to their offspring in order to improve adaptability. Hox genes are important regulatory factors of transcription in metazoans, and are involved in the growth and development of organisms. However, the mechanisms of Microcystis on the expression of Hox genes in Daphnia are unclear. In this study, the effects of Microcystis aeruginosa on Hox gene expression in the mothers and offspring (F1) of two Daphnia similoides sinensis clones were investigated using a mixed diet of M. aeruginosa and Scenedesmus obliquus. Compared with the 100%S food treatment, the survival rates at the end of the experiment of clone 1‐F1 in the food treatments containing M. aeruginosa were significantly lower, but it was significantly higher for clone 2‐F1 in the 20%M + 80%S food treatment. Moreover, the survival rates at the end of the experiment of clone 1‐F1 in the food treatments containing M. aeruginosa were significantly higher than those of their mother. Based on previous transcriptome data, 14 Hox genes of D. similoides sinensis were identified, including Abd‐B, CDX‐1, Dll, HOX‐1, HOX‐2, HOXA1, HOXA2, HOXB3, HOXB3‐2, HOXB7, HOXC4, HOXC7, HOXC8, and HOXD10. The expressions of Abd‐B, HOX‐2, HOXA1, HOXC7, and HOXD10 of clone 2‐mothers in the 40%M + 60%S food treatment were 2.9–22.5 times as high as in the 100%S food treatment, whereas the expressions of CDX‐1, HOX‐1, HOXB3, and HOXD10 of clone 1‐mothers were 4.8–13.1 times at same food level. The expression of HOXA2, HOXC7, HOXC8, and HOXD10 of clone 1‐F1 in the 40%M + 60%S food treatment was 8.2–21.1 times as high as in the 100%S food treatment. However, compared with the 100%S food treatment, the expressions of CDX‐1 in the mothers and F1 of clone 2 and HOXB7 in the mothers of clone 1 in the food treatments containing M. aeruginosa were significantly lower (p < .05). Our results suggest that the offspring (F1) produced by D. similoides sinensis mother pre‐exposed to toxic M. aeruginosa had stronger adaptability to M. aeruginosa than their mothers. Moreover, Hox gene expressions of D. similoides sinensis had obvious differences between clones under stress of toxic M. aeruginosa.

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

Acute effects of UVB radiation on the survival, growth, development, and reproduction of Daphniopsis tibetana Sars (Crustacea: Cladocera)

Daphniopsis tibetana Sars lives in elevation, usually with strong solar UV radiation. We speculate that UV may have an effect on the ecology and evolutionary biology of this species. However, the regulatory effect and mechanism of UV on D. tibetana have not been studied previously. Here, our results showed that UVB could act as a positive factor in the relative body lengths, reproductive parameters, and population growth parameters of D. tibetana when UVB radiation is 20-170 mJ cm^-2, compared with the control group. Strikingly, these parameters were highest at 120 mJ cm^-2. To explore the mechanism underlying the UVB irradiation effects, we conducted a transcriptome analysis using the Trinity platform. The results indicated that differentially regulated genes were mostly enriched in lipid transport and lipid localization by Gene Ontology (GO) enrichment analysis of 146 differentially expressed genes (83 upregulated and 63 downregulated). This is the first study of UVB radiation of D. tibetana to reveal genes that may have crucial roles in survival, growth, and reproduction and could be candidates for future functional studies. Additionally, the study could supply a substantial resource for investigating and elucidating lipids that could play important roles in a physiological context.

Transcriptional and Translational Relationship in Environmental Stress: RNAseq and ITRAQ Proteomic Analysis Between Sexually Reproducing and Parthenogenetic Females in Moina micrura

Moina micrura is a kind of small-bodied water flea within the family Moinidae. Similar to Daphnia, M. micrura could also switch its reproduction mode from parthenogenetic female (PF) to sexual female (SF) to adapt to the external environment. To uncover the mechanisms of reproductive switching in M. micrura, we used both RNA-Seq and iTRAQ analyses to investigate the differentially expressed genes (DEGs) and their protein products between SF and PF in M. micrura. A total of 1665 DEGs (702 up-regulated, 963 down-regulated) and 600 differentially expressed proteins (DEPs) (102 up-regulated, 498 down-regulated) were detected in SF. Correlation analyses indicated that 31 genes were expressed significantly differentially at both transcriptomic and proteomic levels, including 15 up-regulated genes and 16 down-regulated genes in SF. Meanwhile, our data also showed that 528 DEPs have discordant expression at transcript level, implying post-transcriptional (including translational) regulation. These top up-regulated genes and their protein products in SF were mainly grouped into the globin-related family, vitellogenin-related family, cuticle-related family, Hsp-related family and methyltransferases-related family, which were all involved in the reproductive switching in Daphnia. In contrast, a cluster of orthologous groups revealed that up-regulated genes and their protein products in PF were strongly associated with the metabolic process, which may be responsible for rapid population proliferation in M. micrura.

Differential gene transcription across the life cycle in Daphnia magna using a new all genome custom-made microarray

Background

Unravelling the link between genes and environment across the life cycle is a challenging goal that requires model organisms with well-characterized life-cycles, ecological interactions in nature, tractability in the laboratory, and available genomic tools. Very few well-studied invertebrate model species meet these requirements, being the waterflea Daphnia magna one of them. Here we report a full genome transcription profiling of D. magna during its life-cycle. The study was performed using a new microarray platform designed from the complete set of gene models representing the whole transcribed genome of D. magna.

Results

Up to 93% of the existing 41,317 D. magna gene models showed differential transcription patterns across the developmental stages of D. magna, 59% of which were functionally annotated. Embryos showed the highest number of unique transcribed genes, mainly related to DNA, RNA, and ribosome biogenesis, likely related to cellular proliferation and morphogenesis of the several body organs. Adult females showed an enrichment of transcripts for genes involved in reproductive processes. These female-specific transcripts were essentially absent in males, whose transcriptome was enriched in specific genes of male sexual differentiation genes, like doublesex.

Conclusion

Our results define major characteristics of transcriptional programs involved in the life-cycle, differentiate males and females, and show that large scale gene-transcription data collected in whole animals can be used to identify genes involved in specific biological and biochemical processes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4725-7) contains supplementary material, which is available to authorized users.

Effects of Microcystis aeruginosa on the expression of nuclear receptor genes in Daphnia similoides sinensis

Nuclear receptor (NR) genes form a conserved superfamily, which is involved in organism metabolism, reproduction, development, homeostasis, and resource allocation. Microcystis aeruginosa can inhibit the growth and reproduction of Daphnia. However, whether M. aeruginosa can affect the expression of Daphnia NR genes is unknown. In total, 18 NRs were identified in this study based on previous Daphnia similoides sinensis transcriptome data. In treatments containing M. aeruginosa, the gene expression of the NR1 subfamily (E75a, E75b, HR3, HR96, NHR-1, HR97a, HR97g, and NHR97) and the NR2 subfamily (RXR, TLL, PNR, and SVP) were down-regulated 59% and 79%, respectively. In treatments containing M. aeruginosa, although the expression of 78% of the genes showed a similar trend in clones 1 and 2, the expression of 42% of the genes in clone 3 showed the opposite trend compared to clones 1 and 2, suggesting that the adaptability and molecular mechanism differ in individuals with different Microcystis tolerance genotypes.

RNA-seq approach to analysis of gene expression profiles in dark green islands and light green tissues of Cucumber mosaic virus-infected Nicotiana tabacum

Dark green islands (DGIs) surrounded by light green tissues (LGTs) are common leaf symptoms of plants that are systemically infected by various viruses that induce leaf mosaic in infected plants. The inoculation of Cucumber mosaic virus (CMV) in Nicotiana tabacum produced a commonly occurring sequence of classic patterns of DGIs and LGTs. Previous studies confirmed that there are significant differences between DGIs and LGTs in terms of physiology, biochemistry and molecular biology, but the mechanisms by which DGIs form remain unclear. To investigate the global gene expression changes that occur in these special tissues, individual differential gene expression tag libraries were constructed from three total RNA samples isolated from DGIs, LGTs and control plants (CK) and were sequenced using an Illumina HiSeqTM 2000. An analysis of differentially expressed genes (DEGs) and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. These analyses revealed the differences between DGIs, LGTs and CK. GO enrichment and KEGG pathway analyses suggested that several pathways related to photosynthesis and chlorophyll metabolism were enriched in DGIs compared to LGTs and CK. Several pathways related to apoptosis were significantly up-regulated in LGTs compared to DGIs. Additionally, we identified sets of DEGs that may be related to the formation or development of DGIs and LGTs. Our systematic analyses provide comprehensive transcriptomic information regarding DGIs and LGTs in CMV-infected N. tabacum. These data will help characterize the detailed mechanisms of DGI and LGT formation.