Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis

Comparative Analysis of Transposable Elements Reveals the Diversity of Transposable Elements in Decapoda and Their Effects on Genomic Evolution

Transposable elements (TEs) are mobile genetic elements that exist in the host genome and exert considerable influence on the evolution of the host genome. Since crustaceans, including decapoda, are considered ideal models for studying the relationship between adaptive evolution and TEs, TEs were identified and classified in the genomes of eight decapoda species and one diplostraca species (as the outgroup) using two strategies, namely homology-based annotation and de novo annotation. The statistics and classification of TEs showed that their proportion in the genome and their taxonomic composition in decapoda were different. Moreover, correlation analysis and transcriptome data demonstrated that there were more PIF-Harbinger TEs in the genomes of Eriocheir sinensis and Scylla paramamosain, and the expression patterns of PIF-Harbingers were significantly altered under air exposure stress conditions. These results signaled that PIF-Harbingers expanded in the genome of E. sinensis and S. paramamosain and might be related to their air exposure tolerance levels. Meanwhile, sequence alignment revealed that some Jockey-like sequences (JLSs) with high similarity to specific regions of the White spot syndrome virus (WSSV) genome existed in all eight decapod species. At the same time, phylogenetic comparison exposed that the phylogenetic tree constructed by JLSs was not in agreement with that of the species tree, and the distribution of each branch was significantly different. The abovementioned results signaled that these WSSV-specific JLSs might transfer horizontally and contribute to the emergence of WSSV. This study accumulated data for expanding research on TEs in decapod species and also provided new insights and future direction for the breeding of stress-resistant and disease-resistant crab breeds.

Chromosome-scale genome assembly of the freshwater cladoceran crustacean Chydorus sphaericus: A resource for discovery of genes responsive to ecological challenges

The genus Chydorus Leach 1816 (family Chydoridae) is a diverse and ecologically important taxon within freshwater ecosystems. Despite having been widely used in ecological, evolutionary and eco-toxicological studies, no high-quality genomic resource is available for any member of the genus. Here, we present a high-quality chromosome-level assembly of the C. sphaericus genome by combining 7.40 Gb (∼ 50 × coverage) PacBio reads, 19.28 Gb (∼ 135 × coverage) Illumina paired-end reads, and 34.04 Gb Hi-C reads. Our genome assembly is approximately 151 Mb, with contig and scaffold N50 lengths of 1.09 Mb and 13.70 Mb, respectively. The assembly captured 94.9% of the complete eukaryotic BUSCO. Repetitive elements accounted for 17.6% of the genome, and 13,549 protein-coding genes were predicted (based on transcriptome sequencing data, ab-initio or homology-based prediction), of which 96.4% have been functionally annotated in the NCBI-NR database. We identified 303 gene families specific to C. sphaericus, mainly families enriched in functions related to immune response, visual senses and detoxification. Interestingly, we also found 53 significantly expanded gene families in C. sphaericus, mostly with functions related to detoxification. This high-quality assembly genome will act as a reference genome for C. sphaericus and benefit studies on functional and comparative genomics of Chydorus and other crustaceans.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

Exploring mitogenome evolution in Branchiopoda (Crustacea) lineages reveals gene order rearrangements in Cladocera

The class Branchiopoda, whose origin dates back to Cambrian, includes ~ 1200 species which mainly occupy freshwater habitats. The phylogeny and systematics of the class have been debated for long time, until recent phylogenomic analyses allowed to better clarify the relationships among major clades. Based on these data, the clade Anostraca (fairy and brine shrimps) is sister to all other branchiopods, and the Notostraca (tadpole shrimps) results as sister group to Diplostraca, which includes Laevicaudata + Spinicaudata (clam shrimps) and Cladoceromorpha (water fleas + Cyclestherida). In the present analysis, thanks to an increased taxon sampling, a complex picture emerges. Most of the analyzed mitogenomes show the Pancrustacea gene order while in several other taxa they are found rearranged. These rearrangements, though, occur unevenly among taxa, most of them being found in Cladocera, and their taxonomic distribution does not agree with the phylogeny. Our data also seems to suggest the possibility of potentially homoplastic, alternative gene order within Daphniidae.