Diversification of mitochondrial genome of Daphnia galeata (Cladocera, Crustacea): Comparison with phylogenetic consideration of the complete sequences of clones isolated from five lakes in Japan

In silico prospecting of the mtDNA of Macrobrachium amazonicum from transcriptome data

BACKGROUND

Macrobrachium amazonicum is a freshwater prawn widely distributed in South America that is undergoing speciation, so the denomination "M. amazonicum complex" is used for it. The mitochondrial cytochrome c oxidase subunit I (COI) gene has been used to elucidate this speciation, but heteroplasmies and pseudogenes have been recorded, making separation difficult. Obtaining genes from cDNA (RNA) rather than genomic DNA is an effective tool to mitigate those two types of occurrences. The aim of this study was to assemble in silico the mitochondrial DNA (mtDNA) of the Amazonian coastal population of M. amazonicum inhabiting the state of Pará.

RESULTS

Sequences were obtained from the prawn's transcriptome using the de novo approach. Six libraries of cDNA from the androgen gland, hepatopancreas, and muscle tissue were used. The mtDNA of M. amazonicum was 14,960 bp in length. It contained 13 protein-coding genes, 21 complete transfer RNAs, and the 12S and 16S subunits of ribosomal RNA. All regions were found on the light strand except tRNAGln, which was on the heavy strand. The control region (D-loop) was not recovered, making for a gap of 793 bp. The cladogram showed the formation of the well-defined Macrobrachium clade, with high support value in the established branches (91-100). The three-dimensional spatial conformation of the mtDNA-encoded proteins showed that most of them were mainly composed of major α-helices that typically shows in those proteins inserted in the membrane (mitochondrial).

CONCLUSIONS

It was possible to assemble a large part of the mitochondrial genome of M. amazonicum in silico using data from other genomes deposited in GenBank and to validate it through the similarities between its COI and 16S genes and those from animals of the same region deposited in GenBank. Depositing the M. amazonicum mtDNA sequences in GenBank may help solve the taxonomic problems recorded for the species, in addition to providing complete sequences of candidate coding genes for use as biomarkers in ecological studies.

Climate has contributed to population diversification of Daphnia galeata across Eurasia

Climate is a fundamental abiotic factor that plays a key role in driving the evolution, distribution and population diversification of species. However, there have been few investigations of genomic signatures of adaptation to local climatic conditions in cladocerans. Here, we have provided the first high-quality chromosome-level genome assembly (~143 Mb, scaffold N50 12.6 Mb) of the waterflea, Daphnia galeata, and investigated genomic variation in 22 populations from Central Europe and Eastern China. Our ecological-niche models suggested that the historic distribution of D. galeata in Eurasia was significantly affected by Quaternary climate fluctuations. We detected pronounced genomic and morphometric divergences between European and Chinese D. galeata populations. Such divergences could be partly explained by genomic signatures of thermal adaptation to distinct climate regimes: a set of candidate single-nucleotide polymorphisms (SNPs) potentially associated with climate were detected. These SNPs were in genes significantly enriched in the Gene ontology terms "determination of adult lifespan" and "translation repressor activity", and especially, mthl5 and SOD1 involved in the IIS pathway, and EIF4EBP2 involved in the target of the rapamycin signalling pathway. Our study indicates that certain alleles might be associated with particular temperature regimes, playing a functional role in shaping the population structure of D. galeata at a large geographical scale. These results highlight the potential role of molecular variation in the response to climate variation, in the context of global climate change.

The complete mitochondrial genome of water flea Ceriodaphnia dubia (Crustacea: Cladocera) NIES strain

Water flea Ceriodaphnia dubia has been widely used for risk assessments of chemicals and environmental contamination. In this study, the complete mitochondrial genome (mitogenome) of this species NIES strain was determined using short-read high throughput and long-read sequencing technologies. The mitogenome of C. dubia was 15,170 bp in length and consisted of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), and 22 transfer RNAs (tRNAs). The gene order was identical to the pattern conserved across crustaceans. The complete mitogenome of the NIES strain will serve as genetical reference in ecological risk assessments in Japan, as well as resources for future phylogenetical studies using cladocerans.

Analysis of the complete mitogenome of Daphnia galeata from the Han River, South Korea: structure comparison and control region evolution

Daphnia galeata is an important plankton in aquatic ecosystems. As a widely distributed species, D. galeata has been found throughout the Holarctic region. Understanding the genetic diversity and evolution of D. galeata requires the accumulation of genetic information from different locations. Even though the mitochondrial genome (mitogenome) sequence of D. galeata has already been reported, little is known about the evolution of its mitochondrial control region. In this study, D. galeata samples were collected from the Han River on the Korean Peninsula and its partial nd2 gene was sequenced for haplotype network analysis. This analysis showed that four clades of D. galeata were present in the Holarctic region. Moreover, the D. galeata examined in this study belonged to clade D and was specific to South Korea. The mitogenome of D. galeata from the Han River showed similar gene content and structure compared to sequences reported from Japan. Furthermore, the structure of control region of the Han River was similar to those of Japanese clones and differed substantially from European clone. Finally, a phylogenetic analysis based on the amino acid sequences of 13 protein-coding genes (PCGs) indicated that D. galeata from the Han River formed a cluster with clones collected from Lakes Kasumigaura, Shirakaba, and Kizaki in Japan. The differences in control region structure and stem and loop structure reflect the different evolutionary directions of the mitogenomes from Asian and European clones. These findings improve our understanding of the mitogenome structure and genetic diversity of D. galeata.

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.

Daphnia japonica sp. nov. (Crustacea: Cladocera) an eastern Palearctic montane species with mitochondrial discordance

The Daphnia longispina complex (Crustacea: Cladocera) contains several keystone freshwater species such as D. longispina O.F. Müller (D. rosea Sars is a junior synonym), D. galeata Sars, D. cucullata Sars, and D. dentifera Forbes. The complex is common throughout the Holarctic, but there are several geographic regions where local forms have been assigned to European species names based on a superficial morphological resemblance. Here we examine the species status of a form that was previously assigned to D. rosea from a montane bog pond on Honshu, Japan. We used two nuclear non-coding loci (nDNA), mitochondrial sequences (the ND2 protein-coding region) and morphology for evidence. The mitochondrial gene evidence supported the existence of a divergent lineage that is more closely related to D. galeata than to D. dentifera. However, morphology and the nuclear DNA data indicated a lineage that is most closely related to D. dentifera. As our evidence supported the existence of a cohesive divergent lineage, we described a new species, Daphnia japonica sp. nov. Recognition of local and subalpine diversity in this group is critical as ongoing anthropogenic disturbance has been associated with introductions, local extirpations, and hybridization.

Complete mitochondrial genomes of three fairy shrimps from snowmelt pools in Japan

Background

Fairy shrimps belong to order Anostraca, class Branchiopoda, subphylum Crustacea, and phylum Arthropoda. Three fairy shrimp species (Eubranchipus uchidai, E. asanumai, and E. hatanakai) that inhabit snowmelt pools are currently known in Japan. Whole mitochondrial genomes are useful genetic information for conducting phylogenetic analyses. Mitochondrial genome sequences for Branchiopoda members are gradually being collated.

Results

Six whole mitochondrial genomes from the three Eubranchipus species are presented here. Eubranchipus species share the anostracan pattern of gene arrangement in their mitochondrial genomes. The mitochondrial genomes of the Eubranchipus species have a higher GC content than those of other anostracans. Accelerated substitution rates in the lineage of Eubranchipus species were observed.

Conclusion

This study is the first to obtain whole mitochondrial genomes for Far Eastern Eubranchipus species. We show that the nucleotide sequences of cytochrome oxidase subunit I and the 16S ribosomal RNA of E. asanumai presented in a previous study were nuclear mitochondrial DNA segments. Higher GC contents and accelerated substitution rates are specific characteristics of the mitochondrial genomes of Far Eastern Eubranchipus. The results will be useful for further investigations of the evolution of Anostraca as well as Branchiopoda.

A new species group from the Daphnia curvirostris species complex (Cladocera: Anomopoda) from the eastern Palaearctic: taxonomy, phylogeny and phylogeography

The eastern Palaearctic is a centre of diversity for freshwater cladocerans (Crustacea), but little is known about the evolution and taxonomy of this diversity. Daphnia curvirostris is a Holarctic species complex that has most of its diversity in the eastern Palaearctic. We examined the phylogeography, rates of evolution and taxonomic status for each clade of the D. curvirostris complex using morphological and genetic evidence from four genes. The cybertaxonomical and morphological evidence supported an eastern Palaearctic clade, with at least four species (described here as the Daphnia korovchinskyi sp. nov. group) having diagnostic morphological characters. We also detected convergent morphological characters in the D. curvirostris complex that provided information about species boundaries. Two of the new species (Daphnia koreana sp. nov. and Daphnia ishidai sp. nov.) are known from single ponds and are threatened by human activity. Divergence time estimates suggested an ancient origin (12–28 Mya) for the D. korovchinskyi group, but these estimates are complicated by the small number of calibration points.

Mitochondrial genome diversity and evolution in Branchiopoda (Crustacea)

BACKGROUND

The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes.

RESULTS

Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values.

CONCLUSIONS

Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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