The genome of the crustacean Parhyale hawaiensis, a model for animal development, regeneration, immunity and lignocellulose digestion

Ago4-piRNA complex is a key component of genomic immune system against transposon expression in Penaeus monodon

Argonaute proteins are key constituents of small RNA-guided regulatory pathways. In crustaceans, members of the AGO subfamily of Argonaute proteins that play vital roles in immune defense are well studied, while proteins of the PIWI subfamily are less established. PmAgo4 of the black tiger shrimp, Penaeus monodon, though phylogenetically clustered with the AGO subfamily, has distinctive roles of the PIWI subfamily in safeguarding the genome from transposon invasion and controlling germ cell development. This study explored a molecular mechanism by which PmAgo4 regulates transposon expression in the shrimp germline. PmAgo4-associated small RNAs were co-immunoprecipitated from shrimp testis lysate using a PmAgo4-specific polyclonal antibody. RNA-seq revealed a majority of 26-27 nt long small RNAs in the PmAgo4-IP fraction suggesting that PmAgo4 is predominantly associated with piRNAs. Mapping of these piRNAs on nucleotide sequences of two gypsy and a mariner-like transposons of P. monodon suggested that most piRNAs were originated from the antisense strand of transposons. Suppression of PmAgo4 expression by a specific dsRNA elevated the expression levels of the three transposons while decreasing the levels of transposon-related piRNAs. Taken together, these results imply that PmAgo4 exerts its suppressive function on transposons by controlling the biogenesis of transposon-related piRNAs and thus, provides a defense mechanism against transposon invasion in shrimp germline cells.

New genes helped acorn barnacles adapt to a sessile lifestyle

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

Genomic Basis of Freshwater Adaptation in the Palaemonid Prawn Genus Macrobrachium: Convergent Evolution Following Multiple Independent Colonization Events

Adaptation to different salinity environments can enhance morphological and genomic divergence between related aquatic taxa. Species of prawns in the genus Macrobrachium naturally inhabit different osmotic niches and possess distinctive lifecycle traits associated with salinity tolerance. This study was conducted to investigate the patterns of adaptive genomic divergence during freshwater colonization in 34 Macrobrachium species collected from four continents; Australia, Asia, North and South America. Genotyping-by-sequencing (GBS) technique identified 5018 loci containing 82,636 single nucleotide polymorphisms (SNPs) that were used to reconstruct a phylogenomic tree. An additional phylogeny was reconstructed based on 43 candidate genes, previously identified as being potentially associated with freshwater adaptation. Comparison of the two phylogenetic trees revealed contrasting topologies. The GBS tree indicated multiple independent continent-specific invasions into freshwater by Macrobrachium lineages following common marine ancestry, as species with abbreviated larval development (ALD), i.e., species having a full freshwater life history, appeared reciprocally monophyletic within each continent. In contrast, the candidate gene tree showed convergent evolution for all ALD species worldwide, forming a single, well-supported clade. This latter pattern is likely the result of common evolutionary pressures selecting key mutations favored in continental freshwater habitats Results suggest that following multiple independent invasions into continental freshwaters at different evolutionary timescales, Macrobrachium taxa experienced adaptive genomic divergence, and in particular, convergence in the same genomic regions with parallel shifts in specific conserved phenotypic traits, such as evolution of larger eggs with abbreviated larval developmental.

Scavenger receptor B2, a type III membrane pattern recognition receptor, senses LPS and activates the IMD pathway in crustaceans

The immune deficiency (IMD) pathway is critical for elevating host immunity in both insects and crustaceans. The IMD pathway activation in insects is mediated by peptidoglycan recognition proteins, which do not exist in crustaceans, suggesting a previously unidentified mechanism involved in crustacean IMD pathway activation. In this study, we identified a Marsupenaeus japonicus B class type III scavenger receptor, SRB2, as a receptor for activation of the IMD pathway. SRB2 is up-regulated upon bacterial challenge, while its depletion exacerbates bacterial proliferation and shrimp mortality via abolishing the expression of antimicrobial peptides. The extracellular domain of SRB2 recognizes bacterial lipopolysaccharide (LPS), while its C-terminal intracellular region containing a cryptic RHIM-like motif interacts with IMD, and activates the pathway by promoting nuclear translocation of RELISH. Overexpressing shrimp SRB2 in Drosophila melanogaster S2 cells potentiates LPS-induced IMD pathway activation and diptericin expression. These results unveil a previously unrecognized SRB2-IMD axis responsible for antimicrobial peptide induction and restriction of bacterial infection in crustaceans and provide evidence of biological diversity of IMD signaling in animals. A better understanding of the innate immunity of crustaceans will permit the optimization of prevention and treatment strategies against the arising shrimp diseases.

Behavioral circatidal rhythms require Bmal1 in Parhyale hawaiensis

Organisms living in the intertidal zone are exposed to a particularly challenging environment. In addition to daily changes in light intensity and seasonal changes in photoperiod and weather patterns, they experience dramatic oscillations in environmental conditions due to the tides. To anticipate tides, and thus optimize their behavior and physiology, animals occupying intertidal ecological niches have acquired circatidal clocks. Although the existence of these clocks has long been known, their underlying molecular components have proven difficult to identify, in large part because of the lack of an intertidal model organism amenable to genetic manipulation. In particular, the relationship between the circatidal and circadian molecular clocks, and the possibility of shared genetic components, has been a long-standing question. Here, we introduce the genetically tractable crustacean Parhyale hawaiensis as a system for the study of circatidal rhythms. First, we show that P. hawaiensis exhibits robust 12.4-h rhythms of locomotion that can be entrained to an artificial tidal regimen and are temperature compensated. Using CRISPR-Cas9 genome editing, we then demonstrate that the core circadian clock gene Bmal1 is required for circatidal rhythms. Our results thus demonstrate that Bmal1 is a molecular link between circatidal and circadian clocks and establish P. hawaiensis as a powerful system to study the molecular mechanisms underlying circatidal rhythms and their entrainment.

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

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

Imaging Parhyale hawaiensis embryogenesis with frequency domain photoacoustic microscopy: A novel tool in developmental biology

We present the application of a low-cost frequency domain photoacoustic (FDPA) microscope for the label-free imaging of live developing embryos of the crustacean model organism Parhyale hawaiensis. By modulating the intensity of a continuous wave laser source at 9.5 MHz, we achieve the excitation of monochromatic PA waves, which are detected to provide amplitude and phase recordings. The data are subsequently processed to generate accurate maximum amplitude projection and surface reconstructions, delineating the morphological features of the embryos with high resolution and contrast. The findings of this study pave the way for the broader adoption of inexpensive PA diagnostic techniques in developmental biology, shedding light on various fundamental processes in established and emerging model organisms.

The amphipod Parhyale hawaiensis as a promising model in ecotoxicology

Near-shore marine/estuarine environments play an important role in the functioning of the marine ecosystem and are extremely vulnerable to the presence of chemical pollution. The ability to investigate the effects of pollution is limited by a lack of model organisms for which sufficient ecotoxicological information is available, and this is particularly true for tropical regions. The circumtropical marine amphipod Parhyale hawaiensis has become an important model organism in various disciplines, and here we summarize the scientific literature regarding the emergence of this model within ecotoxicology. P. hawaiensis is easily cultured in the laboratory and standardized ecotoxicity protocols have been developed and refined (e.g., miniaturized), and effects of toxicants on acute toxicity (Cd, Cu, Zn, Ag, ammonia, dyes, pesticides, environmental samples), genotoxicity as comet assay/micronuclei, and gene expression (Ag ion and Ag nanoparticles) and regeneration (pesticides) have been published. Methods for determination of internal concentrations of metals (Cu and Ag) and organic substances (synthetic dye) in hemolymph were successfully developed providing sources for the establishment of toxicokinetics models in aquatic amphipods. Protocols to evaluate reproduction and growth, for testing immune responses and DNA damage in germ cells are under way. The sensitivity of P. hawaiensis, measured as 50% lethal concentration (LC50), is in the same range as other amphipods. The combination of feasibility to culture P. hawaiensis in laboratory, the recent protocols for ecotoxicity evaluation and the rapidly expanding knowledge on its biology make it especially attractive as a model organism and promising tool for risk assessment evaluations in tropical environments.

"Omics" data unveil early molecular response underlying limb regeneration in the Chinese mitten crab, Eriocheir sinensis

Limb regeneration is a fascinating and medically interesting trait that has been well preserved in arthropod lineages, particularly in crustaceans. However, the molecular mechanisms underlying arthropod limb regeneration remain largely elusive. The Chinese mitten crab Eriocheir sinensis shows strong regenerative capacity, a trait that has likely allowed it to become a worldwide invasive species. Here, we report a chromosome-level genome of E. sinensis as well as large-scale transcriptome data during the limb regeneration process. Our results reveal that arthropod-specific genes involved in signal transduction, immune response, histone methylation, and cuticle development all play fundamental roles during the regeneration process. Particularly, Innexin2-mediated signal transduction likely facilitates the early stage of the regeneration process, while an effective crustacean-specific prophenoloxidase system (ProPo-AS) plays crucial roles in the initial immune response. Collectively, our findings uncover novel genetic pathways pertaining to arthropod limb regeneration and provide valuable resources for studies on regeneration from a comparative perspective.

Convergent evolution of barnacles and molluscs sheds lights in origin and diversification of calcareous shell and sessile lifestyle

The calcareous shell and sessile lifestyle are the representative phenotypes of many molluscs, which happen to be present in barnacles, a group of unique crustaceans. The origin of these phenotypes is unclear, but it may be embodied in the convergent genetics of such distant groups (interphylum). Herein, we perform comprehensive comparative genomics analysis in barnacles and molluscs, and reveal a genome-wide strong convergent molecular evolution between them, including coexpansion of biomineralization and organic matrix genes for shell formation, and origination of lineage-specific orphan genes for settlement. Notably, the expanded biomineralization gene encoding alkaline phosphatase evolves a novel, highly conserved motif that may trigger the origin of barnacle shell formation. Unlike molluscs, barnacles adopt novel organic matrices and cement proteins for shell formation and settlement, respectively, and their calcareous shells have potentially originated from the cuticle system of crustaceans. Therefore, our study corroborates the idea that selection pressures driving convergent evolution may strongly act in organisms inhabiting similar environments regardless of phylogenetic distance. The convergence signatures shed light on the origin of the shell and sessile lifestyle of barnacles and molluscs. In addition, notable non-convergence signatures are also present and may contribute to morphological and functional specificities.

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

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

The Evolution of Glycoside Hydrolase Family 1 in Insects Related to Their Adaptation to Plant Utilization

Insects closely interact with plants with multiple genes involved in their interactions. β-glucosidase, constituted mainly by glycoside hydrolase family 1 (GH1), is a crucial enzyme in insects to digest plant cell walls and defend against natural enemies with sequestered plant metabolites. To gain more insights into the role of this enzyme in plant-insect interactions, we analyzed the evolutionary history of the GH1 gene family with publicly available insect genomes. We found that GH1 is widely present in insects, while the gene numbers are significantly higher in insect herbivores directly feeding on plant cell walls than in other insects. After reconciling the insect GH1 gene tree with a species tree, we found that the patterns of duplication and loss of GH1 genes differ among insect orders, which may be associated with the evolution of their ecology. Furthermore, the majority of insects' GH1 genes were tandem-duplicated and subsequently went through neofunctionalization. This study shows the evolutionary history of an important gene family GH1 in insects and facilitates our understanding of the evolution of insect-plant interactions.

Crustacean leg regeneration restores complex microanatomy and cell diversity

Animals can regenerate complex organs, yet this process frequently results in imprecise replicas of the original structure. In the crustacean Parhyale, embryonic and regenerating legs differ in gene expression dynamics but produce apparently similar mature structures. We examine the fidelity of Parhyale leg regeneration using complementary approaches to investigate microanatomy, sensory function, cellular composition, and cell molecular profiles. We find that regeneration precisely replicates the complex microanatomy and spatial distribution of external sensory organs and restores their sensory function. Single-nuclei sequencing shows that regenerated and uninjured legs are indistinguishable in terms of cell-type composition and transcriptional profiles. This remarkable fidelity highlights the ability of organisms to achieve identical outcomes via distinct processes.

The organizing role of Wnt signaling pathway during arthropod posterior growth

Wnt signaling pathways are recognized for having major roles in tissue patterning and cell proliferation. In the last years, remarkable progress has been made in elucidating the molecular and cellular mechanisms that underlie sequential segmentation and axial elongation in various arthropods, and the canonical Wnt pathway has emerged as an essential factor in these processes. Here we review, with a comparative perspective, the current evidence concerning the participation of this pathway during posterior growth, its degree of conservation among the different subphyla within Arthropoda and its relationship with the rest of the gene regulatory network involved. Furthermore, we discuss how this signaling pathway could regulate segmentation to establish this repetitive pattern and, at the same time, probably modulate different cellular processes precisely coupled to axial elongation. Based on the information collected, we suggest that this pathway plays an organizing role in the formation of the body segments through the regulation of the dynamic expression of segmentation genes, via controlling the caudal gene, at the posterior region of the embryo/larva, that is necessary for the correct sequential formation of body segments in most arthropods and possibly in their common segmented ancestor. On the other hand, there is insufficient evidence to link this pathway to axial elongation by controlling its main cellular processes, such as convergent extension and cell proliferation. However, conclusions are premature until more studies incorporating diverse arthropods are carried out.

Distinct gene expression dynamics in developing and regenerating crustacean limbs

Regenerating animals have the ability to reproduce body parts that were originally made in the embryo and subsequently lost due to injury. Understanding whether regeneration mirrors development is an open question in most regenerative species. Here, we take a transcriptomics approach to examine whether leg regeneration shows similar temporal patterns of gene expression as leg development in the embryo, in the crustacean Parhyale hawaiensis. We find that leg development in the embryo shows stereotypic temporal patterns of gene expression. In contrast, the dynamics of gene expression during leg regeneration show a higher degree of variation related to the physiology of individual animals. A major driver of this variation is the molting cycle. We dissect the transcriptional signals of individual physiology and regeneration to obtain clearer temporal signals marking distinct phases of leg regeneration. Comparing the transcriptional dynamics of development and regeneration we find that, although the two processes use similar sets of genes, the temporal patterns in which these genes are deployed are different and cannot be systematically aligned.

Identification and classification of cis-regulatory elements in the amphipod crustacean Parhyale hawaiensis

Emerging research organisms enable the study of biology that cannot be addressed using classical 'model' organisms. New data resources can accelerate research in such animals. Here, we present new functional genomic resources for the amphipod crustacean Parhyale hawaiensis, facilitating the exploration of gene regulatory evolution using this emerging research organism. We use Omni-ATAC-seq to identify accessible chromatin genome-wide across a broad time course of Parhyale embryonic development. This time course encompasses many major morphological events, including segmentation, body regionalization, gut morphogenesis and limb development. In addition, we use short- and long-read RNA-seq to generate an improved Parhyale genome annotation, enabling deeper classification of identified regulatory elements. We discover differential accessibility, predict nucleosome positioning, infer transcription factor binding, cluster peaks based on accessibility dynamics, classify biological functions and correlate gene expression with accessibility. Using a Minos transposase reporter system, we demonstrate the potential to identify novel regulatory elements using this approach. This work provides a platform for the identification of novel developmental regulatory elements in Parhyale, and offers a framework for performing such experiments in other emerging research organisms.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01302-6.

The crustacean model Parhyale hawaiensis

Arthropods are the most abundant and diverse animals on earth. Among them, pancrustaceans are an ancient and morphologically diverse group, comprising a wide range of aquatic and semi-aquatic crustaceans as well as the insects, which emerged from crustacean ancestors to colonize most terrestrial habitats. Within insects, Drosophila stands out as one of the most powerful animal models, making major contributions to our understanding of development, physiology and behavior. Given these attributes, crustaceans provide a fertile ground for exploring biological diversity through comparative studies. However, beyond insects, few crustaceans are developed sufficiently as experimental models to enable such studies. The marine amphipod Parhyale hawaiensis is currently the best established crustacean system, offering year-round accessibility to developmental stages, transgenic tools, genomic resources, and established genetics and imaging approaches. The Parhyale research community is small but diverse, investigating the evolution of development, regeneration, aspects of sensory biology, chronobiology, bioprocessing and ecotoxicology.

Measuring concentrations of a dye in the hemolymph of a marine amphipod: Development of a protocol for exposure assessment

The increasing pollution of aquatic environments due to old and emerging contaminants requires the development of integrative methods for exposure assessment. Internal concentrations are a reliable way to estimate total exposure of contaminants originated from different routes (water, sediment, and food). We developed a protocol to evaluate the concentration of a dye, C.I. Disperse Red 1, in the hemolymph of Parhyale hawaiensis, a marine amphipod. LOD and LOQ were satisfactory to detect the dye in all hemolymph samples. The concentration detected in the hemolymph varied related to exposure time and dye concentration (0.003 to 0.086 μg mL^-1). Polynomial regression model was the best fit. The protocol was reliable to detect and quantify dye exposure in marine amphipods and can be considered for future assessments of estuarine and marine regions under the influence of dye processing plants. The method possibly can be easily adapted to other amphipods and other azo dyes.

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

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

The NGS Magic Pudding: A Nanopore-Led Long-Read Genome Assembly for the Commercial Australian Freshwater Crayfish, Cherax destructor

Cherax destructor, the yabby, is an iconic Australian freshwater crayfish species, which, similar to other major invertebrate groups, is grossly under-represented in genomic databases. The yabby is also the principal commercial freshwater crustacean species in Australia subject to explotation via inland fisheries and aquaculture. To address the genomics knowledge gap for this species and explore cost effective and efficient methods for genome assembly, we generated 106.8 gb of Nanopore reads and performed a long-read only assembly of the Cherax destructor genome. On a mini-server configured with an ultra-fast swap space, the de novo assembly took 131 h (∼5.5 days). Genome polishing with 126.3 gb of PCR-Free Illumina reads generated an assembled genome size of 3.3 gb (74.6% BUSCO completeness) with a contig N50 of 80,900 bp, making it the most contiguous for freshwater crayfish genome assemblies. We found an unusually large number of cellulase genes within the yabby genome which is relevant to understanding the nutritional biology, commercial feed development, and ecological role of this species and crayfish more generally. These resources will be useful for genomic research on freshwater crayfish and our methods for rapid and super-efficient genome assembly will have wide application.

The embryonic transcriptome of Parhyale hawaiensis reveals different dynamics of microRNAs and mRNAs during the maternal-zygotic transition

Parhyale hawaiensis has emerged as the crustacean model of choice due to its tractability, ease of imaging, sequenced genome, and development of CRISPR/Cas9 genome editing tools. However, transcriptomic datasets spanning embryonic development are lacking, and there is almost no annotation of non-protein-coding RNAs, including microRNAs. We have sequenced microRNAs, together with mRNAs and long non-coding RNAs, in Parhyale using paired size-selected RNA-seq libraries at seven time-points covering important transitions in embryonic development. Focussing on microRNAs, we annotate 175 loci in Parhyale, 88 of which have no known homologs. We use these data to annotate the microRNAome of 37 crustacean genomes, and suggest a core crustacean microRNA set of around 61 sequence families. We examine the dynamic expression of microRNAs and mRNAs during the maternal-zygotic transition. Our data suggest that zygotic genome activation occurs in two waves in Parhyale with microRNAs transcribed almost exclusively in the second wave. Contrary to findings in other arthropods, we do not predict a general role for microRNAs in clearing maternal transcripts. These data significantly expand the available transcriptomics resources for Parhyale, and facilitate its use as a model organism for the study of small RNAs in processes ranging from embryonic development to regeneration.

Widespread retention of ohnologs in key developmental gene families following whole-genome duplication in arachnopulmonates

Whole-genome duplications (WGDs) have occurred multiple times during animal evolution, including in lineages leading to vertebrates, teleosts, horseshoe crabs, and arachnopulmonates. These dramatic events initially produce a wealth of new genetic material, generally followed by extensive gene loss. It appears, however, that developmental genes such as homeobox genes, signaling pathway components and microRNAs are frequently retained as duplicates (so-called ohnologs) following WGD. These not only provide the best evidence for WGD, but an opportunity to study its evolutionary consequences. Although these genes are well studied in the context of vertebrate WGD, similar comparisons across the extant arachnopulmonate orders are patchy. We sequenced embryonic transcriptomes from two spider species and two amblypygid species and surveyed three important gene families, Hox, Wnt, and frizzled, across these and 12 existing transcriptomic and genomic resources for chelicerates. We report extensive retention of putative ohnologs, further supporting the ancestral arachnopulmonate WGD. We also found evidence of consistent evolutionary trajectories in Hox and Wnt gene repertoires across three of the six arachnopulmonate orders, with interorder variation in the retention of specific paralogs. We identified variation between major clades in spiders and are better able to reconstruct the chronology of gene duplications and losses in spiders, amblypygids, and scorpions. These insights shed light on the evolution of the developmental toolkit in arachnopulmonates, highlight the importance of the comparative approach within lineages, and provide substantial new transcriptomic data for future study.

Marine amphipods (Parhyale hawaiensis) as an alternative feed for the lined seahorse (Hippocampus erectus, Perri 1810): nutritional value and feeding trial

Finding new alternatives to traditional live preys such as Artemia and rotifers, which do not always promote optimal fish growth and survival, is required for the successful aquaculture of highly specialized predatory species, including seahorses. The present study assessed the nutritional value of an interesting marine amphipod (Parhyale hawaiensis), and evaluates through a feeding trial its potential use as a natural prey for 10-months lined seahorse, Hippocampus erectus. P. hawaiensis showed high levels of valuable lipids (20.4–26.7% on dry matter basis) and polyunsaturated fatty acids (PUFAs) ( 26.4–41% of total FAs), including the long-chain PUFAs (LC-PUFAs) arachidonic acid (ARA) (2.9–7.7%), eicosapentaenoic acid (EPA) (4.3–6.5%) and docosahexaenoic acid (DHA) (2.1–6.2%). A comparison between wild-captured and cultured amphipods revealed a significant improvement of the amphipod FA profile in terms of DHA%, total omega-3 (n3) FAs and n3/n6 ratio when employing both a conventional amphipod culture based on a commercial shrimp diet, and, to a lesser extent, a large (3,500 L) biofloc system. Seahorses fed with frozen/wild amphipods, either singly or in combination with Artemia enriched with Super Selco® (INVE Aquaculture, Belgium) for 57 days, substantially improved seahorse growth and FA profiles in terms of ARA, EPA and DHA%, including indices associated to marine sources, such as Σn3 and n3/n6, compared to a diet based solely on enriched Artemia. These results support the use of marine amphipods as an alternative food organism for juvenile H. erectus and suggest a potential use for general marine aquaculture.

Mitochondrial Genomes, Phylogenetic Associations, and SNP Recovery for the Key Invasive Ponto-Caspian Amphipods in Europe

The Ponto-Caspian region is the main donor of invasive amphipods to freshwater ecosystems, with at least 13 species successfully established in European inland waters. Dikerogammarus spp. and Pontogammarus robustoides are among the most successful, due to their strong invasive impact on local biota. However, genomic knowledge about these invaders is scarce, while phylogeography and population genetics have been based on short fragments of mitochondrial markers or nuclear microsatellites. In this study, we provide: (i) a reconstruction of six mitogenomes for four invasive gammarids (D. villosus, D. haemobaphes, D. bispinosus, and P. robustoides); (ii) a comparison between the structure of the newly obtained mitogenomes and those from the literature; (iii) SNP calling rates for individual D. villosus and D. haemobaphes from different invasion sites across Europe; and (iv) the first time-calibrated full mitogenome phylogeny reconstruction of several Ponto-Caspian taxa. We found that, in comparison to other gammarids, the mitogenomes of Ponto-Caspian species show a translocation between the tRNA-E and tRNA-R positions. Phylogenetic reconstruction using the mitogenomes identified that Ponto-Caspian gammarids form a well-supported group that originated in the Miocene. Our study supports paraphyly in the family Gammaridae. These provided mitogenomes will serve as vital genetic resources for the development of new markers for PCR-based identification methods and demographic studies.

Evolutionary transition to XY sex chromosomes associated with Y-linked duplication of a male hormone gene in a terrestrial isopod

Sex chromosomes are highly variable in some taxonomic groups, but the evolutionary mechanisms underlying this diversity are not well understood. In terrestrial isopod crustaceans, evolutionary turnovers in sex chromosomes are frequent, possibly caused by Wolbachia, a vertically-transmitted endosymbiont causing male-to-female sex reversal. Here, we use surgical manipulations and genetic crosses, plus genome sequencing, to examine sex chromosomes in the terrestrial isopod Trachelipus rathkei. Although an earlier cytogenetics study suggested a ZZ/ZW sex chromosome system in this species, we surprisingly find multiple lines of evidence that in our study population, sex is determined by an XX/XY system. Consistent with a recent evolutionary origin for this XX/XY system, the putative male-specific region of the genome is small. The genome shows evidence of Y-linked duplications of the gene encoding the androgenic gland hormone, a major component of male sexual differentiation in isopods. Our analyses also uncover sequences horizontally acquired from past Wolbachia infections, consistent with the hypothesis that Wolbachia may have interfered with the evolution of sex determination in T. rathkei. Overall, these results provide evidence for the co-occurrence of multiple sex chromosome systems within T. rathkei, further highlighting the relevance of terrestrial isopods as models for the study of sex chromosome evolution.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

DNA methylation atlas and machinery in the developing and regenerating annelid Platynereis dumerilii

BACKGROUND

Methylation of cytosines in DNA (5mC methylation) is a major epigenetic modification that modulates gene expression and constitutes the basis for mechanisms regulating multiple aspects of embryonic development and cell reprogramming in vertebrates. In mammals, 5mC methylation of promoter regions is linked to transcriptional repression. Transcription regulation by 5mC methylation notably involves the nucleosome remodeling and deacetylase complex (NuRD complex) which bridges DNA methylation and histone modifications. However, less is known about regulatory mechanisms involving 5mC methylation and their function in non-vertebrate animals. In this paper, we study 5mC methylation in the marine annelid worm Platynereis dumerilii, an emerging evolutionary and developmental biology model capable of regenerating the posterior part of its body post-amputation.

RESULTS

Using in silico and experimental approaches, we show that P. dumerilii displays a high level of DNA methylation comparable to that of mammalian somatic cells. 5mC methylation in P. dumerilii is dynamic along the life cycle of the animal and markedly decreases at the transition between larval to post-larval stages. We identify a full repertoire of mainly single-copy genes encoding the machinery associated with 5mC methylation or members of the NuRD complex in P. dumerilii and show that this repertoire is close to the one inferred for the last common ancestor of bilaterians. These genes are dynamically expressed during P. dumerilii development and regeneration. Treatment with the DNA hypomethylating agent Decitabine impairs P. dumerilii larval development and regeneration and has long-term effects on post-regenerative growth.

CONCLUSIONS

Our data reveal high levels of 5mC methylation in the annelid P. dumerilii, highlighting that this feature is not specific to vertebrates in the bilaterian clade. Analysis of DNA methylation levels and machinery gene expression during development and regeneration, as well as the use of a chemical inhibitor of DNA methylation, suggest an involvement of 5mC methylation in P. dumerilii development and regeneration. We also present data indicating that P. dumerilii constitutes a promising model to study biological roles and mechanisms of DNA methylation in non-vertebrate bilaterians and to provide new knowledge about evolution of the functions of this key epigenetic modification in bilaterian animals.

A chromosome-level assembly of the black tiger shrimp (Penaeus monodon) genome facilitates the identification of growth-associated genes

To salvage marine ecosystems from fishery overexploitation, sustainable and efficient aquaculture must be emphasized. The knowledge obtained from available genome sequence of marine organisms has accelerated marine aquaculture in many cases. The black tiger shrimp (Penaeus monodon) is one of the most prominent cultured penaeid shrimps (Crustacean) with an average annual global production of half a million tons in the last decade. However, its currently available genome assemblies lack the contiguity and completeness required for accurate genome annotation due to the highly repetitive nature of the genome and technical difficulty in extracting high-quality, high-molecular weight DNA. Here, we report the first chromosome-level whole-genome assembly of P. monodon. The combination of long-read Pacific Biosciences (PacBio) and long-range Chicago and Hi-C technologies enabled a successful assembly of this first high-quality genome sequence. The final assembly covered 2.39 Gb (92.3% of the estimated genome size) and contained 44 pseudomolecules, corresponding to the haploid chromosome number. Repetitive elements occupied a substantial portion of the assembly (62.5%), the highest of the figures reported among crustacean species. The availability of this high-quality genome assembly enabled the identification of genes associated with rapid growth in the black tiger shrimp through the comparison of hepatopancreas transcriptome of slow-growing and fast-growing shrimps. The results highlighted several growth-associated genes. Our high-quality genome assembly provides an invaluable resource for genetic improvement and breeding penaeid shrimp in aquaculture. The availability of P. monodon genome enables analyses of ecological impact, environment adaptation and evolution, as well as the role of the genome to protect the ecological resources by promoting sustainable shrimp farming.

The American lobster genome reveals insights on longevity, neural, and immune adaptations

The American lobster, Homarus americanus, is integral to marine ecosystems and supports an important commercial fishery. This iconic species also serves as a valuable model for deciphering neural networks controlling rhythmic motor patterns and olfaction. Here, we report a high-quality draft assembly of the H. americanus genome with 25,284 predicted gene models. Analysis of the neural gene complement revealed extraordinary development of the chemosensory machinery, including a profound diversification of ligand-gated ion channels and secretory molecules. The discovery of a novel class of chimeric receptors coupling pattern recognition and neurotransmitter binding suggests a deep integration between the neural and immune systems. A robust repertoire of genes involved in innate immunity, genome stability, cell survival, chemical defense, and cuticle formation represents a diversity of defense mechanisms essential to thrive in the benthic marine environment. Together, these unique evolutionary adaptations contribute to the longevity and ecological success of this long-lived benthic predator.

Comparative genomics of the coconut crab and other decapod crustaceans: exploring the molecular basis of terrestrial adaptation

Background

The complex life cycle of the coconut crab, Birgus latro, begins when an obligate terrestrial adult female visits the intertidal to hatch zoea larvae into the surf. After drifting for several weeks in the ocean, the post-larval glaucothoes settle in the shallow subtidal zone, undergo metamorphosis, and the early juveniles then subsequently make their way to land where they undergo further physiological changes that prevent them from ever entering the sea again. Here, we sequenced, assembled and analyzed the coconut crab genome to shed light on its adaptation to terrestrial life. For comparison, we also assembled the genomes of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus, and the short-tailed marine-living red king crab, Paralithodes camtschaticus. Our selection of the latter two organisms furthermore allowed us to explore parallel evolution of the crab-like form in anomurans.

Results

All three assembled genomes are large, repeat-rich and AT-rich. Functional analysis reveals that the coconut crab has undergone proliferation of genes involved in the visual, respiratory, olfactory and cytoskeletal systems. Given that the coconut crab has atypical mitochondrial DNA compared to other anomurans, we argue that an abundance of kif22 and other significantly proliferated genes annotated with mitochondrial and microtubule functions, point to unique mechanisms involved in providing cellular energy via nuclear protein-coding genes supplementing mitochondrial and microtubule function. We furthermore detected in the coconut crab a significantly proliferated HOX gene, caudal, that has been associated with posterior development in Drosophila, but we could not definitively associate this gene with carcinization in the Anomura since it is also significantly proliferated in the ornate spiny lobster. However, a cuticle-associated coatomer gene, gammacop, that is significantly proliferated in the coconut crab, may play a role in hardening of the adult coconut crab abdomen in order to mitigate desiccation in terrestrial environments.

Conclusion

The abundance of genomic features in the three assembled genomes serve as a source of hypotheses for future studies of anomuran environmental adaptations such as shell-utilization, perception of visual and olfactory cues in terrestrial environments, and cuticle sclerotization. We hypothesize that the coconut crab exhibits gene proliferation in lieu of alternative splicing as a terrestrial adaptation mechanism and propose life-stage transcriptomic assays to test this hypothesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07636-9.

The Chinese mitten crab genome provides insights into adaptive plasticity and developmental regulation

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general.

First annotated draft genomes of nonmarine ostracods (Ostracoda, Crustacea) with different reproductive modes

Ostracods are one of the oldest crustacean groups with an excellent fossil record and high importance for phylogenetic analyses but genome resources for this class are still lacking. We have successfully assembled and annotated the first reference genomes for three species of nonmarine ostracods; two with obligate sexual reproduction (Cyprideis torosa and Notodromas monacha) and the putative ancient asexual Darwinula stevensoni. This kind of genomic research has so far been impeded by the small size of most ostracods and the absence of genetic resources such as linkage maps or BAC libraries that were available for other crustaceans. For genome assembly, we used an Illumina-based sequencing technology, resulting in assemblies of similar sizes for the three species (335-382 Mb) and with scaffold numbers and their N50 (19-56 kb) in the same orders of magnitude. Gene annotations were guided by transcriptome data from each species. The three assemblies are relatively complete with BUSCO scores of 92-96. The number of predicted genes (13,771-17,776) is in the same range as Branchiopoda genomes but lower than in most malacostracan genomes. These three reference genomes from nonmarine ostracods provide the urgently needed basis to further develop ostracods as models for evolutionary and ecological research.

Chitin Synthesis and Degradation in Crustaceans: A Genomic View and Application

Chitin is among the most important components of the crustacean cuticular exoskeleton and intestinal peritrophic matrix. With the progress of genomics and sequencing technology, a large number of gene sequences related to chitin metabolism have been deposited in the GenBank database in recent years. Here, we summarized the genes and pathways associated with the biosynthesis and degradation of chitins in crustaceans based on genomic analyses. We found that chitin biosynthesis genes typically occur in single or two copies, whereas chitin degradation genes are all multiple copies. Moreover, the chitinase genes are significantly expanded in most crustacean genomes. The gene structure and expression pattern of these genes are similar to those of insects, albeit with some specific characteristics. Additionally, the potential applications of the chitin metabolism genes in molting regulation and immune defense, as well as industrial chitin degradation and production, are also summarized in this review.

A Survey of Transposon Landscapes in the Putative Ancient Asexual Ostracod Darwinula stevensoni

How asexual reproduction shapes transposable element (TE) content and diversity in eukaryotic genomes remains debated. We performed an initial survey of TE load and diversity in the putative ancient asexual ostracod Darwinula stevensoni. We examined long contiguous stretches of DNA in clones from a genomic fosmid library, totaling about 2.5 Mb, and supplemented these data with results on TE abundance and diversity from an Illumina draft genome. In contrast to other TE studies in putatively ancient asexuals, which revealed relatively low TE content, we found that at least 19% of the fosmid dataset and 26% of the genome assembly corresponded to known transposons. We observed a high diversity of transposon families, including LINE, gypsy, PLE, mariner/Tc, hAT, CMC, Sola2, Ginger, Merlin, Harbinger, MITEs and helitrons, with the prevalence of DNA transposons. The predominantly low levels of sequence diversity indicate that many TEs are or have recently been active. In the fosmid data, no correlation was found between telomeric repeats and non-LTR retrotransposons, which are present near telomeres in other taxa. Most TEs in the fosmid data were located outside of introns and almost none were found in exons. We also report an N-terminal Myb/SANT-like DNA-binding domain in site-specific R4/Dong non-LTR retrotransposons. Although initial results on transposable loads need to be verified with high quality draft genomes, this study provides important first insights into TE dynamics in putative ancient asexual ostracods.

A chromosome-level genome of the mud crab (Scylla paramamosain estampador) provides insights into the evolution of chemical and light perception in this crustacean

Mud crabs, found throughout the Indo-Pacific region, are coastal species that are important fisheries resources in many tropical and subtropical Asian countries. Here, we present a chromosome-level genome assembly of a mud crab (Scylla paramamosain). The genome is 1.55 Gb (contig N50 191 kb) in length and encodes 17,821 proteins. The heterozygosity of the assembled genome was estimated to be 0.47%. Effective population size analysis suggested that an initial large population size of this species was maintained until 200 thousand years ago. The contraction of cuticle protein and opsin genes compared with Litopenaeus vannamei is assumed to be correlated with shell hardness and light perception ability, respectively. Furthermore, the analysis of three chemoreceptor gene families, the odorant receptor (OR), gustatory receptor (GR) and ionotropic receptor (IR) families, suggested that the mud crab has no OR genes and shows a contraction of GR genes and expansion of IR genes. The numbers of the three gene families were similar to those in three other decapods but different from those in two nondecapods and insects. In addition, IRs were more diversified in decapods than in nondecapod crustaceans, and most of the expanded IRs in the mud crab genome were clustered with the antennal IR clades. These findings suggested that IRs might exhibit more diverse functions in decapods than in nondecapods, which may compensate for the smaller number of GR genes. Decoding the S. paramamosain genome not only provides insight into the genetic changes underpinning ecological traits but also provides valuable information for improving the breeding and aquaculture of this species.

Sex Determination and Differentiation in Decapod and Cladoceran Crustaceans: An Overview of Endocrine Regulation

Mechanisms underlying sex determination and differentiation in animals are known to encompass a diverse array of molecular clues. Recent innovations in high-throughput sequencing and mass spectrometry technologies have been widely applied in non-model organisms without reference genomes. Crustaceans are no exception. They are particularly diverse among the Arthropoda and contain a wide variety of commercially important fishery species such as shrimps, lobsters and crabs (Order Decapoda), and keystone species of aquatic ecosystems such as water fleas (Order Branchiopoda). In terms of decapod sex determination and differentiation, previous approaches have attempted to elucidate their molecular components, to establish mono-sex breeding technology. Here, we overview reports describing the physiological functions of sex hormones regulating masculinization and feminization, and gene discovery by transcriptomics in decapod species. Moreover, this review summarizes the recent progresses of studies on the juvenile hormone-driven sex determination system of the branchiopod genus Daphnia, and then compares sex determination and endocrine systems between decapods and branchiopods. This review provides not only substantial insights for aquaculture research, but also the opportunity to re-organize the current and future trends of this field.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Chromosome-scale genome assembly of the sea louse Caligus rogercresseyi by SMRT sequencing and Hi-C analysis

Caligus rogercresseyi, commonly known as sea louse, is an ectoparasite copepod that impacts the salmon aquaculture in Chile, causing losses of hundreds of million dollars per year. In this study, we report a chromosome-scale assembly of the sea louse (C. rogercresseyi) genome based on single-molecule real-time sequencing (SMRT) and proximity ligation (Hi-C) analysis. Coding RNAs and non-coding RNAs, and specifically long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) were identified through whole transcriptome sequencing from different life stages. A total of 23,686 protein-coding genes and 12,558 non-coding RNAs were annotated. In addition, 6,308 lncRNAs and 5,774 miRNAs were found to be transcriptionally active from larvae to adult stages. Taken together, this genomic resource for C. rogercresseyi represents a valuable tool to develop sustainable control strategies in the salmon aquaculture industry.

Phylogenomics of the Hyalella amphipod species-flock of the Andean Altiplano

Species diversification in ancient lakes has enabled essential insights into evolutionary theory as they embody an evolutionary microcosm compared to continental terrestrial habitats. We have studied the high-altitude amphipods of the Andes Altiplano using mitogenomic, nuclear ribosomal and single-copy nuclear gene sequences obtained from 36 Hyalella genomic libraries, focusing on species of the Lake Titicaca and other water bodies of the Altiplano northern plateau. Results show that early Miocene South American lineages have recently (late Pliocene or early Pleistocene) diversified in the Andes with a striking morphological convergence among lineages. This pattern is consistent with the ecological opportunities (access to unoccupied resources, initial relaxed selection on ecologically-significant traits and low competition) offered by the lacustrine habitats established after the Andean uplift.

Global Analysis of Gene Expression Profiles Provides Novel Insights into the Development and Evolution of the Large Crustacean Eriocheir sinensis

Chinese mitten crab (Eriocheir sinensis) is an important aquaculture species in Crustacea. Functional analysis, although essential, has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, transcriptome sequencing was conducted on 59 samples representing diverse developmental stages (fertilized eggs, zoea, megalopa, three sub-stages of larvae, juvenile crabs, and adult crabs) and different tissues (eyestalk, hepatopancreas, and muscle from juvenile crabs, and eyestalk, hepatopancreas, muscle, heart, stomach, gill, thoracic ganglia, intestine, ovary, and testis from adult crabs) of E. sinensis. A comprehensive reference transcriptome was assembled, including 19,023 protein-coding genes. Hierarchical clustering based on 128 differentially expressed cuticle-related genes revealed two distinct expression patterns during the early larval developmental stages, demonstrating the distinct roles of these genes in “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting. Phylogenetic analysis of 1406 one-to-one orthologous gene families identified from seven arthropod species and Caenorhabditis elegans strongly supported the hypothesis that Malacostraca and Branchiopoda do not form a monophyletic group. Furthermore, Branchiopoda is more phylogenetically closely related to Hexapoda, and the clade of Hexapoda and Branchiopoda and the clade of Malacostraca belong to the Pancrustacea. This study offers a high-quality transcriptome resource for E. sinensis and demonstrates the evolutionary relationships of major arthropod groups. The differentially expressed genes identified in this study facilitate further investigation of the cuticle-related gene expression networks which are likely associated with “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting.

A transcription factor glial cell missing (Gcm) in the freshwater crayfish Pacifastacus leniusculus

The transcription factor glial cell missing, Gcm, is known to be an important protein in the determination of glial cell fate as well as embryonic plasmatocyte differentiation in Drosophila melanogaster. So far, no function for Gcm in crustaceans has been reported. In this study, we show the cDNA sequence of a Gcm homologue in the freshwater crayfish Pacifastacus leniusculus. The P. leniusculus Gcm transcript is expressed exclusively in brain and nervous tissue, and by in situ hybridization we show that the expression is restricted to a small number of large cells with morphology similar to neurosecretory cells. Furthermore, we show that the expression of Gcm coincides with the expression of a Repo homologue, that is induced in expression by Gcm in Drosophila. Moreover, the Gcm transcript is increased shortly and transiently after injection of cystamine, a substance that inhibits transglutaminase and also strongly affects the movement behavior of crayfish. This finding of Gcm transcripts in a subpopulation of brain cells in very low numbers may enable more detailed studies about Gcm in adult crustaceans.

Knockout of crustacean leg patterning genes suggests that insect wings and body walls evolved from ancient leg segments

The origin of insect wings has long been debated. Central to this debate is whether wings are a novel structure on the body wall resulting from gene co-option, or evolved from an exite (outgrowth; for example, a gill) on the leg of an ancestral crustacean. Here, we report the phenotypes for the knockout of five leg patterning genes in the crustacean Parhyale hawaiensis and compare these with their previously published phenotypes in Drosophila and other insects. This leads to an alignment of insect and crustacean legs that suggests that two leg segments that were present in the common ancestor of insects and crustaceans were incorporated into the insect body wall, moving the proximal exite of the leg dorsally, up onto the back, to later form insect wings. Our results suggest that insect wings are not novel structures, but instead evolved from existing, ancestral structures.

Comparative Genomics of Strictly Vertically Transmitted, Feminizing Microsporidia Endosymbionts of Amphipod Crustaceans

Microsporidia are obligate intracellular eukaryotic parasites of vertebrates and invertebrates. Microsporidia are usually pathogenic and undergo horizontal transmission or a mix of horizontal and vertical transmission. However, cases of nonpathogenic microsporidia, strictly vertically transmitted from mother to offspring, have been reported in amphipod crustaceans. Some of them further evolved the ability to feminize their nontransmitting male hosts into transmitting females. However, our understanding of the evolution of feminization in microsporidia is hindered by a lack of genomic resources. We report the sequencing and analysis of three strictly vertically transmitted microsporidia species for which feminization induction has been demonstrated (Nosema granulosis) or is strongly suspected (Dictyocoela muelleri and Dictyocoela roeselum), along with a draft genome assembly of their host Gammarus roeselii. Contrary to horizontally transmitted microsporidia that form environmental spores that can be purified, feminizing microsporidia cannot be easily isolated from their host cells. Therefore, we cosequenced symbiont and host genomic DNA and devised a computational strategy to obtain genome assemblies for the different partners. Genomic comparison with feminizing Wolbachia bacterial endosymbionts of isopod crustaceans indicated independent evolution of feminization in microsporidia and Wolbachia at the molecular genetic level. Feminization thus represents a remarkable evolutionary convergence of eukaryotic and prokaryotic microorganisms. Furthermore, a comparative genomics analysis of microsporidia allowed us to identify several candidate genes for feminization, involving functions such as DNA binding and membrane fusion. The genomic resources we generated contribute to establish Gammarus roeselii and its microsporidia symbionts as a new model to study the evolution of symbiont-mediated feminization.

A genome-wide screening for RNAi pathway proteins in Acari

BACKGROUND

RNA interference (RNAi) is a highly conserved, sequence-specific gene silencing mechanism present in Eukaryotes. Three RNAi pathways are known, namely micro-RNA (miRNA), piwi-interacting RNA (piRNA) and short interfering RNA (siRNA). However, little knowledge exists about the proteins involved in these pathways in Acari. Moreover, variable successes has been obtained in gene knockdown via siRNA pathway in their functional genomics and management. We hypothesized that the clue may be in the variability of the composition and the efficacy of siRNA machinery among Acari.

RESULTS

Both comparative genomic analyses and domain annotation suggest that all the analyzed species have homologs of putative core proteins that mediate cleaving of targeted genes via the three RNAi pathways. We identified putative homologs of Caenorhabditis elegans RNA-dependent RNA polymerase (RdRP) protein in all species though no secondary Argonaute homologs that operate with this protein in siRNA amplification mechanism were found, suggesting that the siRNA amplification mechanism present in Acari may be distinct from that described in C. elegans. Moreover, the genomes of these species do not encode homologs of C. elegans systemic RNAi defective-1 (Sid-1) protein that mediate silencing of the mRNA target throughout the treated organisms suggesting that the phenomena of systemic RNAi that has been reported in some Acari species probably occur through a different mechanism. However, homologs of putative RNAi spreading defective-3 (Rsd-3) protein and scavenger receptors namely Eater and SR-CI that mediate endocytosis cellular update of dsRNA in C. elegans and Drosophila melanogaster were found in Acari genomes. This result suggests that cellular dsRNA uptake in Acari is endocytosis-dependent. Detailed phylogenetic analyses of core RNAi pathway proteins in the studied species revealed that their evolution is compatible with the proposed monophyletic evolution of this group.

CONCLUSIONS

Our analyses have revealed the potential activity of all three pathways in Acari. Still, much experimental work remains to be done to confirm the mechanisms behind these pathways in particular those that govern systemic/parental RNAi and siRNA amplification in Acari. Disclosure of these mechanisms will facilitate the development of new and specific management tools for the harmful species and enrichment of the beneficial species.

Two sets of candidate crustacean wing homologues and their implication for the origin of insect wings

The origin of insect wings is a biological mystery that has fascinated scientists for centuries. Identification of tissues homologous to insect wings from lineages outside of Insecta will provide pivotal information to resolve this conundrum. Here, through expression and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) functional analyses in Parhyale, we show that a gene network similar to the insect wing gene network (preWGN) operates both in the crustacean terga and in the proximal leg segments, suggesting that the evolution of a preWGN precedes the emergence of insect wings, and that from an evo-devo perspective, both of these tissues qualify as potential crustacean wing homologues. Combining these results with recent wing origin studies in insects, we discuss the possibility that both tissues are crustacean wing homologues, which supports a dual evolutionary origin of insect wings (that is, novelty through a merger of two distinct tissues). These outcomes have a crucial impact on the course of the intellectual battle between the two historically competing wing origin hypotheses.

Pharmacological profiling of stretch activated channels in proprioceptive neurons

Proprioception in mammals and invertebrates occurs through stretch activated ion channels (SACs) localized in sensory endings. In mammals, the primary organs for proprioception are the intrafusal muscle spindles embedded within extrafusal muscle. In invertebrates there are varied types of sensory organs, from chordotonal organs spanning joints to muscle receptor organs (MRO) which are analogous to the mammalian muscle spindles that monitor stretch of muscle fibers. A subset of SACs are the PIEZO channels. They are comprised of a distinct type of protein sequence and are similar among species, from mammals to invertebrates. We screened several new agents (YODA 1, JEDI 2, OB 1 and DOOKU) which have been identified to act on SACs of the PIEZO 1 subtype. JEDI 2 increased activity in the crayfish MRO but not the crab chordotonal organs. The SACs of the crustacean proprioceptors have not been satisfactorily pharmacologically classified, nor has their molecular makeup been identified. We screened these pharmacological agents on model sensory organs in crustaceans to learn more about their subtype classification and compare genomic profiles of related species.

------Widespread conservation and lineage-specific diversification of genome-wide DNA methylation patterns across arthropods

Cytosine methylation is an ancient epigenetic modification yet its function and extent within genomes is highly variable across eukaryotes. In mammals, methylation controls transposable elements and regulates the promoters of genes. In insects, DNA methylation is generally restricted to a small subset of transcribed genes, with both intergenic regions and transposable elements (TEs) depleted of methylation. The evolutionary origin and the function of these methylation patterns are poorly understood. Here we characterise the evolution of DNA methylation across the arthropod phylum. While the common ancestor of the arthropods had low levels of TE methylation and did not methylate promoters, both of these functions have evolved independently in centipedes and mealybugs. In contrast, methylation of the exons of a subset of transcribed genes is ancestral and widely conserved across the phylum, but has been lost in specific lineages. A similar set of genes is methylated in all species that retained exon-enriched methylation. We show that these genes have characteristic patterns of expression correlating to broad transcription initiation sites and well-positioned nucleosomes, providing new insights into potential mechanisms driving methylation patterns over hundreds of millions of years.