The genome of the marine monogonont rotifer Brachionus plicatilis: Genome-wide expression profiles of 28 cytochrome P450 genes in response to chlorpyrifos and 2-ethyl-phenanthrene

Single and combined effects of increased temperature and methylmercury on different stages of the marine rotifer Brachionus plicatilis

Rapid, anthropogenic activity-induced global warming is a severe problem that not only raises water temperatures but also shifts aquatic environments by increasing the bioavailability of heavy metals (HMs), with potentially complicated effects on aquatic organisms, including small aquatic invertebrates. For this paper, we investigated the combined effects of temperature (23 and 28 °C) and methylmercury (MeHg) by measuring physiological changes, bioaccumulation, oxidative stress, antioxidants, and the mitogen-activated protein kinase signaling pathway in the marine rotifer Brachionus plicatilis. High temperature and MeHg adversely affected the survival rate, lifespan, and population of rotifers, and bioaccumulation, oxidative stress, and biochemical reactions depended on the developmental stage, with neonates showing higher susceptibility than adults. These findings demonstrate that increased temperature enhances potentially toxic effects from MeHg, and susceptibility differs with the developmental stage. This study provides a comprehensive understanding of the combined effects of elevated temperature and MeHg on rotifers. ENVIRONMENTAL IMPLICATION: Methylmercury (MeHg) is a widespread and harmful heavy metal that can induce lethal effects on aquatic organisms in even trace amounts. The toxicity of metals can vary depending on various environmental conditions. In particular, rising temperatures are considered a major factor affecting bioavailability and toxicity by changing the sensitivity of organisms. However, there are few studies on the combinational effects of high temperatures and MeHg on aquatic animals, especially invertebrates. Our research would contribute to understanding the actual responses of aquatic organisms to complex aquatic environments.

Application of Moina mongolica and Daphniopsis tibetana in marine ecotoxicology studies in China: A review

Moina mongolica and Daphniopsis tibetana are typical saline Cladocera in China that are characterized by a wide salinity range, rapid reproduction, and high-density culture. In this paper, the characteristics and life history of M. mongolica and D. tibetana are reviewed. The application of these two species of Cladocera to ecotoxicology in recent years is also summarized from the aspects of environmental factors and environmental pollutants, including ultraviolet B radiation, temperature, salinity, alkalinity, pH, heavy metals, harmful red tide, pesticides, and persistent organic pollutants. Additionally, the toxicity sensitivity of saline Cladocera in different reproductive statuses and inter-generational embryos is discussed. Finally, the need to enhance knowledge of the molecular genomics, population dynamics, and strategies for protection of saline Cladocera, along with the need for establishment of estuarine and marine environmental monitoring standards are discussed. Overall, this review highlights the potential for using these Cladocera species as indicator organisms for estuarine and marine ecotoxicology studies.

Changes in transcriptome regulations of a marine rotifer Brachionus plicatilis under methylmercury stress

Mercury (Hg), a heavy metal pollutant worldwide, can be transformed into methylmercury (MeHg) by various aquatic microorganisms in water, thus accumulating along the aquatic food chain and posing a particular challenge to human health. Zooplankton plays a crucial role in aquatic ecosystems and serves as a major component of the food chain. To evaluate the effects of MeHg on the rotifer Brachionus plicatilis and reveal the underlying mechanism of these effects, we exposed B. plicatilis to MeHg by either direct immersion or by feeding with MeHg-poisoned Chlorella pyrenoidesa, respectively, and conducted a transcriptomic analysis. The results showed that B. plicatilis directly exposed to MeHg by immersion showed significant enrichment of the glutathione metabolism pathway for detoxification of MeHg. In addition, the exposure to MeHg by feeding induced a significant enrichment of lysosome and notch signaling pathways of rotifers, supporting the hypothesis that MeHg can induce autophagy dysfunction in cells and disturb the nervous system of rotifers. In two different routes of MeHg exposure, the pathway of cytochrome P450 in rotifers showed significant enrichment for resisting MeHg toxicity. Our results suggest further studies on the potential mechanism and biological responses of MeHg toxicity in other links of the aquatic food chain.

Metabolic pathways of methylmercury in rotifer Brachionus plicatilis

Methylmercury (MeHg) readily accumulates in aquatic organisms while transferring and amplifying in the aquatic food chains. This study firstly explores the in vivo accumulation sites and metabolic regulation of MeHg in the rotifer Brachionus plicatilis by aggregation-induced emission fluorogen (AIEgen) and metabolomics. Fluorescent image analysis by AIEgen showed that MeHg in B. plicatilis mainly occured in the ciliary corona, esophagus, mastax, stomach and intestine in the direct absorption group. In the other group, where B. plicatilis were indirectly supplied with MeHg via food intake, the accumulation of MeHg in the rotifer occurred in the ciliary corona, various digestive organs, and the pedal gland. However, the MeHg accumulated in the rotifer is difficult to metabolize outside the body. Metabolomics analysis showed that the significant enrichment of ABC transporters was induced by the direct exposure of rotifers to dissolved MeHg. In contrast, exposure of rotifers to MeHg via food intake appeared to influence carbon, galactose, alanine, aspartate and glutamate metabolisms. Besides, the disturbed biological pathways such as histidine metabolism, beta-alanine metabolism and pantothenate and CoA biosynthesis in rotifers may be associated with L-aspartic acid upregulation in the feeding group. The significant enrichment of ABC transporters and carbon metabolism in rotifers may be related to the accumulation of MeHg in the intestine of rotifers. In both pathways of MeHg exposure, the arginine biosynthesis and metabolism of rotifers were disturbed, which may support the hypothesis that rotifers produce more energy to resist MeHg toxicity. This study provides new insight into the accumulation and toxicity mechanisms of MeHg on marine invertebrates from the macro and micro perspectives.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

A genome-wide screening of the 70 kDa heat shock protein (HSP70) genes in the rotifer Brachionus plicatilis sensu stricto with a characterization of two heat-inducible HSP70 genes

The 70 kDa heat shock proteins (HSP70s) and the constitutive members of the HSP70 family (heat shock cognates; HSC70s) play essential roles in various biological processes. The number of hsp70/hsc70 in the database is rapidly increasing because of their importance and the automatic annotation of newly sequenced genomes. However, accumulating evidence indicates that neither hsp70 nor hsc70 forms a monophyletic gene family, raising the need to reconsider the annotation strategy based on the traditional concept of the inducible HSP70 and constitutive HSC70s. The main aim of this study is to establish a systematic scheme to annotate hsp70-like genes taking the latest phylogenetic insights into account. We cloned two hsp70s from the rotifer Brachionus plicatilis sensu stricto (s.s.), an emerging model in evolutionary genetics, and conducted a genome-wide screening of B. plicatilis s.s. hsp70s using the two sequences as queries. A total of 15 hsp70-like genes were found, and 7 of them encoded distant members of the HSP70 family, the function of which largely remains unknown. Eight canonical hsp70s were annotated according to a recently proposed nomenclature based on the molecular evolution: e.g., HSP70cA1/B1 for the cytosolic lineage, HSP70er1 for the endoplasmic reticulum lineage, and HSP70m1 for the mitochondrial lineage. The two cloned hsp70s, HSP70cB1 and HSP70cB2, ubiquitously increased their mRNA levels up to 7.5 times after heat treatment as assessed by semi-quantitative PCR, real-time PCR, and in situ hybridization. This systematic screening incorporating a reasonable update to the annotation strategy would provide a useful example for future HSP70 studies, especially those in non-traditional model organisms.

Highly efficient CRISPR-mediated gene editing in a rotifer

Rotifers have been studied in the laboratory and field for over 100 years in investigations of microevolution, ecological dynamics, and ecotoxicology. In recent years, rotifers have emerged as a model system for modern studies of the molecular mechanisms of genome evolution, development, DNA repair, aging, life history strategy, and desiccation tolerance. However, a lack of gene editing tools and transgenic strains has limited the ability to link genotype to phenotype and dissect molecular mechanisms. To facilitate genetic manipulation and the creation of reporter lines in rotifers, we developed a protocol for highly efficient, transgenerational, CRISPR-mediated gene editing in the monogonont rotifer Brachionus manjavacas by microinjection of Cas9 protein and synthetic single-guide RNA into the vitellaria of young amictic (asexual) females. To demonstrate the efficacy of the method, we created knockout mutants of the developmental gene vasa and the DNA mismatch repair gene mlh3. More than half of mothers survived injection and produced offspring. Genotyping these offspring and successive generations revealed that most carried at least 1 CRISPR-induced mutation, with many apparently mutated at both alleles. In addition, we achieved precise CRISPR-mediated knock-in of a stop codon cassette in the mlh3 locus, with half of injected mothers producing F2 offspring with an insertion of the cassette. Thus, this protocol produces knockout and knock-in CRISPR/Cas9 editing with high efficiency, to further advance rotifers as a model system for biological discovery.

Effects of temperature and combinational exposures on lipid metabolism in aquatic invertebrates

Studies of changes in fatty acids in response to environmental temperature changes have been conducted in many species, particularly mammals. However, few studies have considered aquatic invertebrates, even though they are particularly vulnerable to changes in environmental temperature. In this review, we summarize the process by which animals synthesize common fatty acids and point out differences between the fatty acid profiles of vertebrates and those of aquatic invertebrates. Unlike vertebrates, some aquatic invertebrates can directly synthesize polyunsaturated fatty acids (PUFAs), which can be used to respond to temperature changes. Various studies have shown that aquatic invertebrates increase the degree of saturation in their fatty acids through an increase in saturated fatty acid production or a decrease in PUFAs as the temperature increases. In addition, we summarize recent studies that have examined the complex effects of temperature and combinational stressors to determine whether the degree of saturation in aquatic invertebrates is influenced by other factors. The combined effects of carbon dioxide partial pressure, food quality, starvation, salinity, and chemical exposures have been confirmed, and fatty acid profile changes in response to high temperature were greater than those from combinational stressors.

Metabolism deficiency and oxidative stress induced by plastic particles in the rotifer Brachionus plicatilis: Common and distinct phenotypic and transcriptomic responses to nano- and microplastics

Growing experimental data on the adverse effects of microplastic pollution on marine biota indicate that the size of the plastic particles is a key determinant of toxicity. Here, we investigated size-dependent toxicity at different levels of biological organizations in the marine rotifer Brachionus plicatilis, from bioaccumulation as an initiating event to adverse in-vivo outcomes, with ecotoxicogenomic approach to elucidate the size-dependent toxicity of microplastics. Nanoplastics strongly retarded the reproduction and population growth of B. plicatilis, while microplastics were associated with moderate effects. This size dependency could be attributed to the selective induction of oxidative stress by nanoplastic exposure in addition to a metabolic deficiency, which was a common toxicity mechanism with both nano- and microplastic exposure as predicted by transcriptomic analysis. Our findings suggested that metabolic deficiency is a shared toxicity mechanism of nano- and microplastics, while oxidative stress might be responsible for the stronger toxicity of nanoplastics.

Oxidative stress-mediated synergistic deleterious effects of nano- and microplastics in the hypoxia-conditioned marine rotifer Brachionus plicatilis

While pollution due to nano- and micro-sized plastics (NMPs) and hypoxic conditions both occur in coastal areas, the deleterious potential of co-exposure to hypoxia and NMPs (hypoxia and micro-sized plastics, HMPs; hypoxia and nano-sized plastics, HNPs) is largely unclear. Here, we provide evidence for multigenerational effects of HMP and HNP in the marine rotifer Brachionus plicatilis by investigating changes in its life traits, antioxidant system, and hypoxia-inducible factor (HIF) pathway using an orthogonal experimental design, with nanoscale and microscale particles measuring 0.05 μm and 6.0 μm in diameter, respectively, and hypoxic conditions of 0.5 mg/L for six generations. Combined exposure to NMPs and hypoxia caused a significant decrease in fecundity and overproduction of reactive oxygen species (ROS). The HIF pathway and circadian clock genes were also significantly upregulated in response to HMP and HNP exposure. In particular, synergistic deleterious effects of HNP were evident, suggesting that size-dependent toxicity can be a major driver of the effects of hypoxia and NMP co-exposure. After several generations of exposure, ROS levels returned to basal levels and transcriptomic resilience was observed, although rotifer reproduction remained suppressed. These findings help eluciating the underlying molecular mechanisms involved in responses to plastic pollution in hypoxic conditions.

Comparative genome analysis of the monogonont marine rotifer Brachionus manjavacas Australian strain: Potential application for ecotoxicology and environmental genomics

This is the first study to analyze the whole-genome sequence of B. manjavacas Australian (Aus.) strain through combination of Oxford Nanopore long-read seq, resulting in a total length of 108.1 Mb and 75 contigs. Genome-wide detoxification related gene families in B. manjavacas Aus. strain were comparatively analyzed with those previously identified in other Brachionus spp., including B. manjavacas German (Ger.) strain. Most of the subfamilies in detoxification related families (CYPs, GSTs, and ABCs) were highly conserved and confirmed orthologous relationship with Brachionus spp., and with accumulation of genome data, clear differences between genomic repertoires were demonstrated the marine and the freshwater species. Furthermore, strain-specific genetic variations were present between the Aus. and Ger. strains of B. manjavacas. This whole-genome analysis provides in-depth review on the genomic structural differences for detoxification-related gene families and further provides useful information for comparative ecotoxicological studies and evolution of detoxification mechanisms in Brachionus spp.

G protein-coupled receptors (GPCRs) in rotifers and cladocerans: Potential applications in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology

The G protein-coupled receptor (GPCR) superfamily plays a fundamental role in both sensory functions and the regulation of homeostasis, and is highly conserved across the eukaryote taxa. Its functional diversity is related to a conserved seven-transmembrane core and invariant set of intracellular signaling mechanisms. The interplay between these properties is key to the evolutionary success of GPCR. As this superfamily originated from a common ancestor, GPCR genes have evolved via lineage-specific duplications through the process of adaptation. Here we summarized information on GPCR gene families in rotifers and cladocerans based on their evolutionary position in aquatic invertebrates and their potential application in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology. Phylogenetic analyses were conducted to examine the evolutionary significance of GPCR gene families and to provide structural insight on their role in aquatic invertebrates. In particular, most GPCR gene families have undergone sporadic evolutionary processes, but some GPCRs are highly conserved across species despite the dynamics of GPCR evolution. Overall, this review provides a better understanding of GPCR evolution in aquatic invertebrates and expand our knowledge of the potential application of these receptors in various fields.

The genome of the freshwater monogonont rotifer Brachionus rubens: Identification of phase I, II, and III detoxification genes

Monogonont rotifers are common species in aquatic environments and make model species for ecotoxicology studies. Whole genomes of several species of the genus Brachionus have been assembled, but no information on the freshwater rotifer Brachionus rubens has been reported. In this study, the whole-genome sequence of B. rubens was successfully assembled using NextDenovo. The total length of the genome was 132.7 Mb (N50 = 2.51 Mb), including 122 contigs. The GC contents accounted for 29.96% of the genome. Aquatic organisms are always exposed to various external stresses, and a comprehensive genomic analysis is needed to better understand the adverse effects on organisms. This paper focuses on the ecotoxicological aspect and conducted genome analysis of representative gene families involved in detoxification mechanisms against environmental stressors. Specifically, we identified cytochrome P450 genes (CYPs) of phase I, glutathione S-transferase genes (GSTs) of phase II, and ATP-binding cassette transporter genes (ABCs) of phase III in the genome of B. rubens. Gene duplications were found in CYP, GST, and ABC genes, as is the case for other Brachionus rotifers. Our results suggest that these detoxification-related gene families have evolved in a species-specific and/or lineage-specific manner. This paper improves our understanding of how the freshwater Brachionus rotifers respond to environmental stressors in a molecular ecotoxicology context.

The Genome of the Marine Rotifer Brachionus manjavacas: Genome-Wide Identification of 310 G Protein-Coupled Receptor (GPCR) Genes

The marine rotifer Brachionus manjavacas is widely used in ecological, ecotoxicological, and ecophysiological studies. The reference genome of B. manjavacas is a good starting point to uncover the potential molecular mechanisms of responses to various environmental stressors. In this study, we assembled the whole-genome sequence (114.1 Mb total, N50 = 6.36 Mb) of B. manjavacas, consisting of 61 contigs with 18,527 annotated genes. To elucidate the potential ligand-receptor signaling pathways in marine Brachionus rotifers in response to environmental signals, we identified 310 G protein-coupled receptor (GPCR) genes in the B. manjavacas genome after comparing them with three other species, including the minute rotifer Proales similis, Drosophila melanogaster, and humans (Homo sapiens). The 310 full-length GPCR genes were categorized into five distinct classes: A (262), B (26), C (7), F (2), and other (13). Most GPCR gene families showed sporadic evolutionary processes, but some classes were highly conserved between species as shown in the minute rotifer P. similis. Overall, these results provide potential clues for in silico analysis of GPCR-based signaling pathways in the marine rotifer B. manjavacas and will expand our knowledge of ligand-receptor signaling pathways in response to various environmental signals in rotifers.

First Insights into the Repertoire of Secretory Lectins in Rotifers

Due to their high biodiversity and adaptation to a mutable and challenging environment, aquatic lophotrochozoan animals are regarded as a virtually unlimited source of bioactive molecules. Among these, lectins, i.e., proteins with remarkable carbohydrate-recognition properties involved in immunity, reproduction, self/nonself recognition and several other biological processes, are particularly attractive targets for biotechnological research. To date, lectin research in the Lophotrochozoa has been restricted to the most widespread phyla, which are the usual targets of comparative immunology studies, such as Mollusca and Annelida. Here we provide the first overview of the repertoire of the secretory lectin-like molecules encoded by the genomes of six target rotifer species: Brachionus calyciflorus, Brachionus plicatilis, Proales similis (class Monogononta), Adineta ricciae, Didymodactylos carnosus and Rotaria sordida (class Bdelloidea). Overall, while rotifer secretory lectins display a high molecular diversity and belong to nine different structural classes, their total number is significantly lower than for other groups of lophotrochozoans, with no evidence of lineage-specific expansion events. Considering the high evolutionary divergence between rotifers and the other major sister phyla, their widespread distribution in aquatic environments and the ease of their collection and rearing in laboratory conditions, these organisms may represent interesting targets for glycobiological studies, which may allow the identification of novel carbohydrate-binding proteins with peculiar biological properties.

Comparative analysis reveals within-population genome size variation in a rotifer is driven by large genomic elements with highly abundant satellite DNA repeat elements

BACKGROUND

Eukaryotic genomes are known to display an enormous variation in size, but the evolutionary causes of this phenomenon are still poorly understood. To obtain mechanistic insights into such variation, previous studies have often employed comparative genomics approaches involving closely related species or geographically isolated populations within a species. Genome comparisons among individuals of the same population remained so far understudied-despite their great potential in providing a microevolutionary perspective to genome size evolution. The rotifer Brachionus asplanchnoidis represents one of the most extreme cases of within-population genome size variation among eukaryotes, displaying almost twofold variation within a geographic population.

RESULTS

Here, we used a whole-genome sequencing approach to identify the underlying DNA sequence differences by assembling a high-quality reference genome draft for one individual of the population and aligning short reads of 15 individuals from the same geographic population including the reference individual. We identified several large, contiguous copy number variable regions (CNVs), up to megabases in size, which exhibited striking coverage differences among individuals, and whose coverage overall scaled with genome size. CNVs were of remarkably low complexity, being mainly composed of tandemly repeated satellite DNA with only a few interspersed genes or other sequences, and were characterized by a significantly elevated GC-content. CNV patterns in offspring of two parents with divergent genome size and CNV patterns in several individuals from an inbred line differing in genome size demonstrated inheritance and accumulation of CNVs across generations.

CONCLUSIONS

By identifying the exact genomic elements that cause within-population genome size variation, our study paves the way for studying genome size evolution in contemporary populations rather than inferring patterns and processes a posteriori from species comparisons.

Genomics and transcriptomics of epizoic Seisonidea (Rotifera, syn. Syndermata) reveal strain formation and gradual gene loss with growing ties to the host

Background

Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensulato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea.

Results

Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria).

Conclusion

Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07857-y.

Species-specific effects of iron on temperate and tropical marine rotifers in reproduction, lipid and ROS metabolisms

Two euryhaline rotifers, the temperate species Brachionus plicatilis and tropical species Brachionus rotundiformis, were used to investigate the effects of iron (FeSO₄·7H₂O), an essential trace metal, on reproductive patterns and lifetables, including the metabolism of lipid and reactive oxygen species (ROS). B. plicatilis was more sensitive to iron with regard to sexual reproduction. While iron had no significant effect on the population growth at 0-48 μg/mL, it caused a decrease in the resting egg production. B. plicatilis exposed to 6 and 12 μg/mL of iron showed an increase in the intracellular ROS levels and a decrease in the neutral lipid content in sexual organs, accompanied by downregulation of antioxidant components CuZnSOD and two cytochromes (CYP clan 2&3). These patterns suggested that iron-induced oxidative stress was not neutralized by its antioxidant defense system, thus negatively affecting the fecundity of fertilized mictic females. However, B. rotundiformis showed a dose-dependent increase in population growth with extended lifespan and positive sexual reproduction in response to 0-24 μg/mL iron. Furthermore, compared to Fe-exposed B. plicatilis, B. rotundiformis showed better antioxidant mechanism, whereas genes involved in lipid synthesis (citrate lyase, mitochondrial CYP) and reproduction (vasa, sirtuin-2) were significantly upregulated compared to the control, implying that B. rotundiformis was likely to have higher resilience in response to iron-induced oxidative stress. These findings suggest that iron is likely to cause interspecific interactions in the B. plicatilis species complex, whereas the tropical species B. rotundiformis may have evolved an effective defense mechanism against iron-induced stress.

The genome of the minute marine rotifer Proales similis: Genome-wide identification of 401 G protein-coupled receptor (GPCR) genes

The minute marine rotifer Proales similis is a potential model species for ecotoxicological and ecophysiological studies. Therefore, the provision of whole-genome data for P. similis is an easy way to deepen understanding of the molecular mechanisms involved in response to various environmental stressors. In this research, we assembled the whole-genome sequence (32.7 Mb total, N50 = 2.42 Mb) of P. similis, consisting of 15 contigs with 10,785 annotated genes. To understand the ligand-receptor signaling pathway in rotifers in response to environmental cues, we identified 401 G protein-coupled receptor (GPCR) genes in the P. similis genome and compared them with those from other species. The 401 full-length GPCR genes were classified into five distinct classes: A (363), B (18), C (7), F (2), and other (11). Most GPCR gene families have undergone sporadic evolutionary processes. However, some classes were highly conserved between species. Overall, this result provides new information about GPCR-based signaling pathways and the evolution of GPCRs in the minute rotifer P. similis, and it expands our knowledge of ligand-receptor signaling pathways in response to various environmental cues.

The genome of the euryhaline rotifer Brachionus paranguensis: Potential use in molecular ecotoxicology

Brachionus spp. rotifers have been proposed as model organisms for ecotoxicological studies. We analyzed the whole-genome sequence of B. paranguensis through NextDenovo, resulting in a total length of 106.2 Mb and 71 contigs. The N50 and the GC content were 4.13 Mb and 28%, respectively. A total of 18,501 genes were predicted within the genome of B. paranguensis. Prominent detoxification-related gene families of phase I and II detoxifications have been investigated. In parallel with other Brachionus rotifers, high gene expansion was observed in CYP clan 3 and GST sigma class in B. paranguensis. Moreover, species-specific expansion of sulfotransferase (SULTs) and gain of UDP-glucuronosyltransferases (UGTs) through horizontal gene transfer has been specifically found within B. plicatilis complex. This whole-genome analysis of B. paranguensis provides a basis for molecular ecotoxicological studies and provides useful information for comparative studies of the evolution of detoxification mechanisms in Brachionus spp.

The genome of the freshwater monogonont rotifer Brachionus angularis: Identification of phase I, II, and III detoxification genes and their roles in molecular ecotoxicology

Brachionus spp. rotifers, which are widely distributed in aquatic environments, have been proposed as model organisms for ecotoxicological studies. Although the genomes of several rotifers belonging to the genus Brachionus have been assembled, the genome for the freshwater rotifer Brachionus angularis remains unknown. In this study, we analyzed the whole-genome sequence of B. angularis, which revealed a total length of 56.5 Mb and 21 contigs. The N50 and the GC content were 5.42 Mb and 23.66%, respectively. A total of 13,952 genes were predicted. Of them, we identified the main detoxification-related gene families, including those for cytochrome P450, glutathione S-transferase (GST), and the ATP-binding cassette transporter. In comparison with other Brachionus rotifers, massive species-specific expansion in GST sigma genes was found in B. angularis. This whole-genome analysis of B. angularis provides a basis for molecular ecotoxicological studies and provides useful biological tools for comparative studies of the evolution of detoxification mechanisms in Brachionus spp.

Application of heavy-ion-beam irradiation to breeding large rotifer

In larviculture facilities, rotifers are generally used as an initial food source, while a proper size of live feeds to connect rotifer and Artemia associated with fish larval growth is needed. The improper management of feed size and density induces mass mortality and abnormal development of fish larvae. To improve the survival and growth of target larvae, this study applied carbon and argon heavy-ion-beam irradiation in mutation breeding to select rotifer mutants with larger lorica sizes. The optimal irradiation conditions of heavy-ion beam were determined with lethality, reproductivity, mutant frequency, and morphometric characteristics. Among 56 large mutants, TYC78, TYC176, and TYA41 also showed active population growth. In conclusion, (1) heavy-ion-beam irradiation was defined as an efficient tool for mutagenesis of rotifers and (2) the aforementioned 3 lines that have larger lorica length and active population growth may be used as a countermeasure of live feed size gap during fish larviculcure.

Evolutionary dynamics of transposable elements in bdelloid rotifers

Transposable elements (TEs) are selfish genomic parasites whose ability to spread autonomously is facilitated by sexual reproduction in their hosts. If hosts become obligately asexual, TE frequencies and dynamics are predicted to change dramatically, but the long-term outcome is unclear. Here, we test current theory using whole-genome sequence data from eight species of bdelloid rotifers, a class of invertebrates in which males are thus far unknown. Contrary to expectations, we find a variety of active TEs in bdelloid genomes, at an overall frequency within the range seen in sexual species. We find no evidence that TEs are spread by cryptic recombination or restrained by unusual DNA repair mechanisms. Instead, we find that that TE content evolves relatively slowly in bdelloids and that gene families involved in RNAi-mediated TE suppression have undergone significant expansion, which might mitigate the deleterious effects of active TEs and compensate for the consequences of long-term asexuality.

Genome-wide identification and structural analysis of heat shock protein gene families in the marine rotifer Brachionus spp.: Potential application in molecular ecotoxicology

Heat shock proteins (Hsp) are class of conserved and ubiquitous stress proteins present in all living organisms from primitive to higher level. Various studies have demonstrated multiple cellular functions of Hsp in living organisms as an important biomarker in response to abiotic and biotic stressors including temperature, salinity, pH, hypoxia, environmental pollutants, and pathogens. However, full understanding on the mechanism and pathway involved in the induction of Hsp still remains challenging, especially in aquatic invertebrates. In this study, the entire Hsp family and subfamily members in the marine rotifers Brachionus spp., one of the cosmopolitan ecotoxicological model organisms, have been genome-widely identified. In Brachionus spp. Hsp family was comprised of Hsp10, small hsp (sHsp), Hsp40, Hsp60, Hsp70/105, and Hsp90, with highest number of genes found within Hsp40 DnaJ homolog subfamily C members. Also, the differences in the orientation of the conserved motifs within Hsp family may have induced differences in transcriptional gene modulation in response to thermal stress in Brachionus koreanus. Overall, Hsp family-specific domains were highly conserved in all three Brachionus spp., relative to Homo sapiens and across other animal taxa and these findings will be helpful for future ecotoxicological studies focusing on Hsps.

Genome-wide identification and transcriptional modulation of histone variants and modification related genes in the low pH-exposed marine rotifer Brachionus koreanus

Histone modification is considered to be a major epigenetic control mechanism. These modifications (e.g. acetylation, phosphorylation, and methylation) may affect the interaction of histones with DNA and/or regulate DNA-based processes (e.g., recombination, repair, replication, and transcription) and chromatin remodeling complexes. Despite their significance in metazoan life and evolution, few studies have been conducted to identify genes undergoing epigenetic control modification in aquatic invertebrates. In this study, we identified whole core histones (70 total genes) and post-translational modification (PTM) histone genes (63 total genes) in the marine rotifer Brachionus koreanus through whole-genome analysis, and annotated them according to the human nomenclature. Notably, upon comparative analysis of cis-regulatory motif sequences, we found that B. koreanus core histone protein structures were similar to those of mammals. Furthermore, to examine the effect of parental low pH stress on the offspring's epigenetic regulation, we investigated the expression of PTM genes in two generations of B. koreanus exposed to low pH conditions. Given that the B. koreanus genome does not possess DNA methyltransferase 1 and 3 genes, we concluded that histone genes could be involved as an important epigenetic mechanism in B. koreanus. Therefore, the histone-associated genes identified in this study could be useful for ecotoxicological studies and facilitate the application of chromatin immunoprecipitation sequencing using high-throughput DNA sequencing based on the genome-wide identification of transcription factor binding sites in rotifers.

Genome-wide identification of DNA double-strand break repair genes and transcriptional modulation in response to benzo[α]pyrene in the monogonont rotifer Brachionus spp

The DNA repair system has evolved from the common ancestor of all life forms and its function is highly conserved within eukaryotes. In this study, to reveal the role of DNA double-strand break repair (DSB) genes in response to benzo[α]pyrene (B[α]P), we first identified DSB genes in relation to homologous recombination and non-homologous end joining events in four Brachionus rotifer spp.: B. calyciflorus, B. koreanus, B. plicatilis, and B. rotundiformis. In all the Brachionus spp., 39 orthologous genes to human DSB repair genes were identified. Furthermore, three genes in B. koreanus, two genes in B. plicatilis, and one gene in B. calyciflorus and B. rotundiformis were present as duplicated genes, indicating that these genes were diversified over speciation in the genus Brachionus. Moreover, we compared DSB repair genes on the gene structures in four monogonont Brachionus rotifers and the bdelloid rotifer Adineta vaga, which possesses highly efficient DNA repair ability. The transcriptional responses of four monogonont Brachionus rotifers in response to B[α]P exposure showed how B[α]P exposure led to DSBs and subsequently recruited DNA DSB repair pathways in the rotifer B. koreanus. Taken together, this study provides a better understanding of the potential role of DSB repair genes in the monogonont rotifer Brachionus spp. in response to B[α]P.

The genome of the marine monogonont rotifer Brachionus rotundiformis and insight into species-specific detoxification components in Brachionus spp

The monogonont rotifer Brachionus spp. have been widely used for ecotoxicological studies because of their advantages as one of the most suitable laboratory experimental species. In the present study, we obtained and assembled the whole genome sequence of the rotifer Brachionus rotundiformis, consisting of 13,612 annotated genes with 213 scaffolds and 58 Mb in total length. Focusing on ecotoxicological aspects, we conducted a comparative genome analysis on the gene families involved in detoxification, including four to six sulfotransferase gene families, seven uridine 5'-diphospho-glucuronosyltransferase gene families, and 58, 61, or 70 ATP-binding cassette genes in the genus Brachionus including Brachionus koreanus and Brachionus plicatilis. Our results suggest that these gene families have undergone a species- and/or lineage-specific evolution in response to the surrounding environmental pressure. Our genome resource for B. rotundiformis would be highly useful for future ecotoxicological studies and also provides a better understanding on the view of evolutionary mechanism of detoxification in the genus Brachionus spp.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Effects of temperature changes on life parameters, oxidative stress, and antioxidant defense system in the monogonont marine rotifer Brachionus plicatilis

Global warming is a big concern for all organisms and many efforts have been made to reveal the potential effects of temperature elevation on aquatic organisms. However, limited studies on molecular mechanistic approaches on physiological effects due to temperature changes are available. Here, we investigated the effects of temperature changes on life parameters (e.g., population growth [total number of rotifers], and lifespan), oxidative stress levels and antioxidant activities (e.g., glutathione S-transferase [GST], catalase [CAT], superoxide dismutase [SOD]) with expression levels in the monogonont marine rotifer Brachionus plicatilis. The changes in temperatures led to significant reduction (P < 0.05) in lifespan, possibly due to significant decrease (P < 0.05) in antioxidant activities, reducing the potential to cope with significant elevation in the temperature-induced oxidative stress in B. plicatilis. To further assess the actual induction and clearance of reactive oxygen species (ROS), N-acetyl-L-cysteine was used to examine whether the temperature-induced oxidative stress could be successfully scavenged. Furthermore, expression patterns of the antioxidant-related genes (GSTs, SODs, and CATs) were down- or upregulated (P < 0.05) in response to different temperatures in B. plicatilis. Overall, these findings indicate that ROS-mediated oxidative stress led to cellular damage and antioxidant defense system, resulting in deleterious effects on life parameters in rotifer.

The genome of the marine rotifer Brachionus koreanus sheds light on the antioxidative defense system in response to 2-ethyl-phenanthrene and piperonyl butoxide

BRACHIONUS: spp. (Rotifera: Monogononta) have been introduced as ecotoxicological model-organisms that are widely distributed in aquatic environments. Among the Brachionus spp., the monogonont rotifer Brachionus koreanus has been widely used for ecology, ecotoxicology, and evolution, thus, providing the whole genome data of B. koreanus is important for further understandings of in-depth molecular mechanisms. In this study, the completed assembly and characterization of the B. koreanus genome resulted in a total length of 85.7 Mb with 14,975 annotated genes. The final number of scaffolds was 567 with an N50 value and a GC content of 1.86 Mb and 24.35 %, respectively. Based on the fully constructed genome database, a total of 24 CYPs, 23 GSTs, two SODs, and a single CAT genes were identified and analyzed antioxidant activities (CAT, SOD, and GST), and transcriptional regulation of the entire CYPs, GSTs, SODs, and CAT in response to 2-ethyl-phenanthrene (2-ethyl-PHE) and piperonyl butoxide (PBO), to demonstrate the usefulness of the whole genome library of B. koreanus in response xenobiotic-induced oxidative stress. The assembled B. koreanus genome will provide a better understanding on the molecular ecotoxicology in the view of molecular mechanisms underlying toxicological responses, particularly on xenobiotic detoxification processes in the rotifer B. koreanus.