Wenhao Li
Faxiang Wang
Shouwen Jiang
Binbin Pan
Qi Liu
Qianghua Xu- 1Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, College of Marine Sciences, Shanghai Ocean University, Shanghai, China
- 2Shanghai Engineering Research Center of Hadal Science and Technology, College of Marine Sciences, Shanghai Ocean University, Shanghai, China
- 3Key Laboratory of Aquaculture Resources and Utilization, Ministry of Education, College of Fisheries and Life Sciences, Shanghai Ocean University, Shanghai, China
- 4National Distant-water Fisheries Engineering Research Center, Shanghai Ocean University, Shanghai, China
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.
Introduction
Embryogenesis is a highly vulnerable developmental stage that effects the entire life cycle of the organism (Scholz et al., 2008) and ensures production of the target phenotype (Arambourou et al., 2017). Zygote-specific protein-coding genes are expressed during embryogenesis and early organogenesis. Embryogenesis is a complex biological process that integrates mitosis, epiboly, and morphogenesis, including actin organization, microtubule assembly, DNA repair, and cell division (Ishii et al., 2004; Lan et al., 2018). During embryogenesis, abiotic stressors, such as temperature, oxygen concentration, environmental chemical compounds, and high hydrostatic pressure (HHP), can induce failure of embryonic development (Young et al., 1997; Villalobos et al., 2006; Mestre et al., 2008). Embryogenesis is also influenced by various endogenous and exogenous factors that determine the distribution and diversification of the species (Scholz et al., 2008; Ángeles-González et al., 2020).
As the most successful group of animals on the planet, arthropods account for the highest biomass and species diversity (Bar-On et al., 2018). Amphipods are among the most common arthropods in aquatic environments, particularly in hadal zones with the depths exceeding 8,000 m (Blankenship et al., 2006). The hadal zone (depth, 6,000–11000 m) is characterized by low temperatures, lack of light, limited food supply, and HHP (Jamieson et al., 2010). As the most defining environmental factor, HHP increases by 1 atm (or about 0.101 MPa) for every 10 m in water depth (Yancey et al., 2014; Yancey, 2020). Generally, HHP can break protein structures, cause DNA damage, and reduce cell membrane fluidity (Brooks, 2014; Lan et al., 2018). To survive in a hadal zone, chemical chaperones, especially organic osmolytes, facilitate adaptation to HHP. For example, trimethylamine oxide (TMAO) is an important osmotic regulator produced by oxidation of trimethylamine (TMA) via flavin-dependent monooxygenases (FMOs) (Canyelles et al., 2018). The concentration of TMAO is reported to gradually increase with depth, suggesting an important function for organisms occupying hadal zones (Yancey et al., 2014; Qin et al., 2021; Liu et al., 2022). The ability of amphipods to thrive in such extreme habitats and embryogenesis of hadal-dwelling amphipods have attracted much interest.
HHP, a dominant feature of hadal zones, is also reported to delay embryonic development (Mestre et al., 2008), increase the risk of malformations (Mestre et al., 2008), and induce chromosome doubling (Zhu et al., 2017) and aberrations (Glover et al., 2020). Previous studies have reported HHP-induced mitogynogenesis of cultured fish species (Yamamoto, 1999; Lin et al., 2016; Zhu et al., 2017; Glover et al., 2020). However, the pressure tolerance of eggs of other aquatic species remains unknown. The development of the amphipod embryo occurs via ectogenesis, which is defined as growth of an organism in an environment outside of the body. The amphipod’ eggs developing normally in hadal zones under HHP is an important evolutionary developmental trait.
However, investigations of embryonic development of hadal amphipods are limited due to low fecundity (Ritchie et al., 2017), small brood size (Grassle & Sanders, 1973), and fasting during gestation (Perrone et al., 2002; Blankenship-Williams & Levin, 2009). In 2001, a hadal lysianassoid female (Uristes chastaini) that was collected with a baited trap while brooding allowed for the first recording of a “live egg” (Blankenship et al., 2006; Blankenship-Williams & Levin, 2009). In 2018, a tiny brood of hatching Alicella gigantea (Crustacea, Amphipoda, Lysianassoidea) was collected from the New Britain Trench (7.02S, 149.16E) at a depth of 8,824 m with the use of an autonomous deep-ocean lander. Detailed information about this lander vehicle and the sample collection process is described by Chan et al. (2020; 2022; 2021). To the best of our knowledge, this was the first collection of A. gigantea eggs from a hadal zone.
A. gigantea is a common inhabitant of hadal zones. In spite of the limited food supply in hadal zones, this species exhibits incredible gigantism, with the adult body length reaching 240–340 mm (Harrison et al., 1983; Barnard & Ingram, 1986; Jamieson et al., 2013). Numerous studies of this “supergiant” amphipod have been conducted, including a mitochondrial genome study of the phylogeny and adaptation mechanisms (Li et al., 2019) and a trace element studies suggested the use of amphipods as potential indicators of trace element bioavailability in hadal zones (Zhu et al., 2022). Transcriptome and evolutionary analyses conducted by Li et al. (2021) identified seven positively selected genes (PSGs) of the “supergiant” A. gigantea involved in inositol phosphate metabolism, insulin signaling, and glycogenesis signaling, which are ultimately related to growth and proliferation. However, the disjunctive geographical distribution (Jamieson et al., 2013; Ritchie et al., 2015), low observation frequency (Jamieson et al., 2013), and limited genomic resources have hindered elucidation of the adaptation of embryos to the environmental extremes of hadal zones (Jamieson et al., 2013).
Owing to the gigantic genome size of 34.79 pg (34.02 Gb) (Ritchie et al., 2017), expected repeat-rich genomic regions or genome duplication (Ritchie et al., 2017; Wang et al., 2021), assembly of the A. gigantea genome remains challenging. Therefore, transcript assembly of the A. gigantea egg is an effective option for genome-wide comparative studies of the molecular evolution of embryogenesis of hadal amphipods. Here, the transcriptome of the A. gigantea egg was sequenced. Further analysis of the morphology, positive selection, gene evolution rate, and concentration determinations of the TAMO of the A. gigantea eggs will facilitate a comprehensive overview of the piezophilic adaptation of embryos in hadal zones and provide a system to study embryogenesis under HHP.
Materials and methods
Sample collection and species identification
During an expedition of the New Britain Trench (8,824 m, 7.02◦S, 149.16◦E) in the West Pacific Ocean conducted in 2016, our group collected 10 A. gigantea individuals with attached eggs with a full ocean depth lander vehicle launched from the “Zhang Jian” research vessel that were immediately stored at −80 °C immediately. Detailed information about the lander vehicle and sampling process is described by Chan et al. (2020; 2022; 2021) and Li et al. (2021). The specimens were stored in bubble chamber full of dry ice and delivered to Shanghai Ocean University.
Post-spawning eggs adhered to the ventral brood pouch of the brooding female amphipods were carefully separated. The species of the eggs and females were identified based on conventional morphological traits and mitochondrial molecular markers. The A. gigantea specimens were identified by mitochondrial cytochrome oxidase subunit Ⅰ (CO I) and 16S rRNA gene markers, which were described by Chan et al. (2020). The COⅠ was amplified using specific primers: COI-F: 5′-GGTCAACAAATCATAAAGATATTGG-3′ and COI-R: 5′-TAAACTTCAGGGTGACCAAAAAATCA-3′. The 16S rRNA gene was amplified using specific primers 16SF: 5′-GACGACAAGACCCTAAAAGC-3′ and 16SR: 5′-CGC TGTTATCCCTAAAGTA-3′.
Egg measurement and scanning electron microscopy
Eggs were separated from the female amphipods and measured under a ZEISS SteREO Discovery. V20 stereo microscope (Carl Zeiss AG, Oberkochen, Germany). Since the eggs of A. gigantea are oval, measurements were made along the major and minor axes. All measurements were repeated at least three times and the mean values were calculated.
The A. gigantea eggs were imaged with a scanning electron microscope. Briefly, the frozen eggs were fixed in 2.5% cold glutaraldehyde solution (4°C) for at least 12 h, as described previously (Tommasini et al., 1989), then dehydrated with a graded series of pre-cooled ethanol (30%, 50%, 70%, 80%, and 90%) for 20 min at 4°C at each concentration. Afterward, the eggs were stored at room temperature and dehydrated with 100% ethanol for 20 min, which was repeated five times. To prevent collapse of the eggshell, the dehydrated eggs were dried with the critical point drying method using a Leica EM CPD300 critical point dryer (Leica Microsystems GmbH, Wetzlar, Germany). The eggs were subsequently coated with gold film using a sputter coater (E-1045, Hitachi High-Tech Science Corporation, Tokyo, Japan) and observed with a scanning electron microscope (SU8010; Hitachi High-Tech Science Corporation).
RNA extraction, transcriptome sequencing and assembly, and gene functional annotation
Total RNA was extracted from the eggs using TRIzol™ Reagent (Thermo Fisher Scientific, Waltham, MA, United States) in accordance with the manufacturer’s instructions. In this study, three biological replicates of each individual were sequenced. A cDNA library was constructed using a VAHTS Universal V6 RNA-seq Library Prep Kit for Illumina (Vazyme Biotech Co., Ltd., Nanjing, China) and paired-end sequenced using a NovaSeq 6,000 System (Illumina, Inc., San Diego, CA, United States) with a read length of 150 bp in accordance with the manufacturer’s instructions.
To obtain high-confidence assembled transcripts, low-quality bases and adapter sequences were trimmed from the raw sequencing data using the Trimmomatic tool (v0.33; http://www.usadellab.org/cms/index.php?page=trimmomatic) with the following parameters: HEADCROP, 3; AVGQUAL, 30; TRAILING, 20; and MINLEN, 36. De novo assembly of the resulting clean reads was performed with Trinity software (v2.4.0; https://github.com/trinityrnaseq/trinityrnaseq/wiki) generated at k-mer = 25 (Haas et al., 2013). The assembled transcripts (contig cut-off > 300 bp) were used for subsequent analysis. The longest isoforms of the genes were selected as unigenes by a custom Perl script to eliminate redundant isoforms. TransDecoder software (https://github.com/TransDecoder/TransDecoder/releases) was used to predict the open reading frames of the non-redundant transcripts. The DIAMOND BLASTX sequence aligner (https://github.com/bbuchfink/diamond) was used to annotate the assembled transcripts of A. gigantea eggs in the National Center for Biotechnology Information (NCBI) non-redundant (NR) database (https://www.ncbi.nlm.nih.gov/refseq/about/nonredundantproteins/) with default parameters (e value = 1e-10). Only the best hit was maintained. All contaminated reads were removed and only transcripts identified as Arthropoda taxa were annotated with InterProScan (https://www.ebi.ac.uk/interpro/), Swiss-Prot (https://www.expasy.org/resources/uniprotkb-swiss-prot), NR, KOG (https://mycocosm.jgi.doe.gov/help/kogbrowser.jsf), CDD (http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml), eggNOG (http://eggnog5.embl.de/#/app/home), Kaas (https://www.genome.jp/kegg/kaas/), and blast2go (v2.5; https://www.blast2go.com/). The gene ontology (GO) annotation results were combined with the results obtained with InterProScan, eggNOG, and blast2go. The Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation results were integrated with the results obtained with Kaas and eggNOG.
Phylogenetic analysis
For construction of a phylogenetic tree, the embryo transcriptome or genome data of 11 Arthropoda species were obtained from multiple sources. The embryo transcriptome data of Gammarus chevreuxi (ERR1097462–ERR1097470) were downloaded from the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra), cleaned, and assembled using the same methodology as for the A. gigantea embryo. In addition, eight genome sequences and gene annotation files (i.e., Hyalella azteca) (Poynton et al., 2018), Penaeus vannamei (Zhang et al., 2019), Armadillidium nasatum, Armadillidium vulgare, Cryptotermes secundus, Daphnia pulex (Colbourne et al., 2011), Eurytemora affinis, and Trinorchestia longiramus) were download from the NCBI website (https://www.ncbi.nlm.nih.gov/genome/), while the genome sequence of Parhyale hawaiensis (Kao et al., 2016) was download from the Janelia Research Campus website (https://research.janelia.org/Pavlopoulos/). Orthologous protein sequences of each gene were grouped using the OrthoMcl genome-scale algorithm (v2.0.9; https://orthomcl.org/orthomcl/app). In total, 610 groups of single-copy orthologs of 11 species were identified with the use of a manual Perl script. MAFFT multiple sequence alignment software (v7.407; https://mafft.cbrc.jp/alignment/software/) was used to align the sequences, which were concentrated into a supergene with a length of 778,175 aa. Gblocks software (0.91b; http://phylogeny.lirmm.fr/phylo_cgi/one_task.cgi?task_type=gblocks) was used to identify the conserved regions. Finally, 80,824 positions (10%) of 1,742 selected blocks were retained. The topology of the phylogenetic tree was inferred with RaxML software (v8.1.24; https://cme.h-its.org/exelixis/web/software/raxml/) using the best model LG + I + G + F in accordance with the Akaike information criterion and Bayesian information criterion, as determined with ProtTest software (v3.4; https://github.com/ddarriba/prottest3) with C. secundus as an outgroup (bootstrap value, 1,000). The graphical viewer Figtree (v1.4.4; http://tree.bio.ed.ac.uk/software/figtree/) was used to visualize the cladogram tree.
Detection of genes exhibiting accelerated evolution and positively selected genes
Single-copy orthologues were selected among A. gigantea, G. chevreuxi, P. hawaiensis, T. longiramus, and H. azteca using OrthoMcl software. Considering the distant divergence time, other species were not searched for single-copy orthologues. PRANK software (https://www.ebi.ac.uk/research/goldman/software/prank/software) was used to align the codon sequences based on a codon substitution matrix. The conserved region was obtained using Gblocks software (0.91b) with the parameter “−t = c”. Signatures of genes exhibiting accelerated evolution and PSGs along a specific branch were individually detected by branch or branch-site models implemented in the CodeML program of the PAML (phylogenetic analysis by maximum likelihood) package (version 4.9 g) (Yang, 2007). In order to identify genes exhibiting accelerated evolution and PSGs in hadal amphipods, the A. gigantea lineage was set as the “foreground” and other shallow-water amphipods as the “background”. The nonsynonymous-synonymous substitution rate ratio dN/dS (ω) was used to detect the selective pressure in whole gene level and with codon frequencies calculated average nucleotide frequencies at the three codon positions (F3x4). Then, Bayes Empirical Bayes (BEB) analysis was used to identify amino acid sites under positive selection. The genes with p-value of the likelihood ratio test less than 0.05 were defined as PSGs. GO and KEGG enrichment analyses were implemented with a hypergeometric test using a custom Python script. Plotting and statistical analysis were performed with R software (https://www.r-project.org/) and in-house Python scripts. The GOChord plot, in which the genes are linked via ribbons to their assigned terms, was generated using GOplot R package (Walter et al., 2015). Orthologous genes of A. gigantea and D. melanogaster were detected based on reciprocal-best-BLAST-hits (RBH). Protein–protein interaction networks (PPIs) were constructed by STRINIG database (https://cn.string-db.org/; Organisms: Drosophila melanogaster). The visualization of PPI networks was implemented using igraph R package (Csardi and Nepusz, 2006).
Determination of trimethylamine oxide and trimethylamine concentrations in eggs
The concentrations of TMAO and TMA in the eggs of A. gigantea, Cherax quadricarinatus, Macrobrachium rosenbergii, and Procambarus clarkii were determined. Measurements of each egg were repeated nine times. Briefly, 20 mg of each egg were added to 1,200 μL of a 4:1 solution of CH3OH and water in a centrifuge tube. Then, the eggs were crushed with a high-throughput tissue crusher. Homogenization of the egg tissue revealed that the A. gigantea eggshell was quite flexible. The egg suspension was ultrasonicated three times in an ice bath, then chilled to −20°C for 20 min and centrifuged at 13,000 rpm and 4°C for 10 min. The supernatant was collected, diluted by 20 fold, and 400-ml aliquots were stored at 4°C until analysis of the TMAO and TMA concentrations by ultra-performance liquid chromatography–tandem mass spectrometry (Waters Corporation, Milford, MA, United States), as described by Liu et al. (2022).
Results
Morphology of A. gigantea eggs
The post-spawning, but not the free-spawning, “supergiant” hadal amphipods bearing eggs (Figure 1A) were classified as A. gigantea based on external morphological traits (Barnard & Ingram, 1986). The eggs collected from the ventral pouch were identified as A. gigantea based on conventional morphological traits and mitochondrial molecular markers. Due to the amphipoda egg-collecting behavior (Sheader, 1996) and frozen conditions, the developmental stages of the collected eggs could not be accurately determined. Since no rudimentary embryo was detected microscopically, it was concluded that the sampled eggs were in the early embryogenesis.
FIGURE 1. Morphological characteristics of the hadal A. gigantea’s eggs. (A). Representative image of A. gigantea with eggs (lateral view and vertical view). Scale bar = 10 mm; (B). The image of one A. gigantea egg. Horizontal and vertical line segments indicate the minor axis and the major axis, respectively. Scale bar = 1 mm; (C). The length of the minor and the major axis of the measured eggs. (D) and (E). The ultrastructure of eggshell by scanning electron microscopy (SEM). Scale bar = 100 μm in D and scale bar = 50 μm in E; (F). The density of the hole diameter in the endochorion of the A. gigantea’s eggs.
A. gigantea eggs are oval (Figures 1A,B) with the minor axis ranging from 3.30 to 8.54 mm and major axis ranging from 6.95 to 14.88 mm (Figure 1C). A. gigantea eggs can also be considered “supergiant” as compared with the eggs of other amphipod species (Serra, 2003). Scanning electron microscopy revealed that the eggshell of A. gigantea was composed of a fibrillar exochorion (Figure 1D) and a porous endochorion layer (Figures 1D,E), similar to the eggshell structure of the fruit flies Ceratitis capitata (Margaritis, 1985) and Rhagoletis cerasi (Mouzaki & Margaritis, 1991). The mean diameter of the pores of the endochorion was >10 μm, with the largest diameter up to 20 μm (Figure 1F), similar to the perforations of the alveolar layer of the eggshell of Triops cancriformis, a crustacean surviving in temporary freshwater pools (Tommasini et al., 1989), but much larger than that of three species of neotropical Eulimnadia (Crustacea: Branchiopoda: Spinicaudata) (Rabet et al., 2012).
Sequencing, assembly statistics, and gene annotation
The transcriptome sequencing results of the A. gigantea eggs are summarized in Supplementary Table S1. Sequencing with the Illumina platform returned 179.4 M reads and 26.9 GB of raw data. After filtering, there were 168.9 M (94.15%) clean reads and 24.8 GB of clean data. The trinity-based de novo transcriptome assembly consisted of 184,706 contigs. The longest isoform was extracted as a unigene, which consisted of 121,238 contigs. The contig N50 consisted of 628 bp with a total length (nt) of 73.69 M and GC content of 40.30%. Among these, 11,403 contigs were matched with GO annotations. The transcripts with “best-hit” homologs in the NCBI NR database were mostly mapped to Malacostraca, the largest of the six classes of crustaceans. Among these transcripts, 23,273 were mapped to H. azteca (Amphipoda), 3,999 to Penaeus vannamei (Decapoda), 643 to Armadillidium vulgare (Isopoda) and 166 to C. secundus (Blattodea).
Phylogenetic relationships
A detailed phylogenetic tree constructed from published embryo transcriptomes and whole genome data combined with the transcriptome data of the A. gigantea eggs clearly showed that G. chevreuxi, a shallow-water amphipod, was most closely related to A. gigantea, followed by P. hawaiensis, T. longiramus, and H. azteca (Figure 2). All four of these species are members of the order Amphipoda and, thus, used as reference for subsequent analysis.
FIGURE 2. Phylogenetic tree of the 11 arthropods.
Genes exhibiting accelerated evolution and positively selected genes of the hadal lineage
Genetic mutations that conferred an advantage to embryogenesis in the hadal environment were maintained throughout the evolutionary process. Rapidly evolving genes and PSGs were identified from 5-way orthologs genes of A. gigantea and shallow-water arthropods, respectively. The lineage of A. gigantea was set as the “foreground” phylogeny, and other nodes as the “background” phylogeny. In total, 161 genes exhibiting accelerated evolution were identified in the eggs of hadal amphipods’ (Supplementary Table S2). Enriched pathways associated with genes exhibiting accelerated evolution of A. gigantea eggs are shown in Figure 3A. Genes exhibiting accelerated evolution were mainly related to the terms “regulation of neuron apoptotic process”, “regulation of neuron death”, “oocyte construction”, “regulation of oocyte development”, “centrosome duplication”, “spindle”, “microtubule-based process”, “microtubule cytoskeleton”, and “cytoskeleton” (Figure 3A). Then, a network of genes exhibiting accelerated evolution was generated, which included genes associated with embryonic development, such as dpp (Figure 3B).
FIGURE 3. Accelerated evolution of the embryogenesis genes in the eggs of A. gigantea. (A). The bar plot shows the GO enrichment results of the accelerated evolution genes; (B). A network of the accelerated evolution genes in the A. gigantea’s eggs.
Overall, 272 PSGs were identified in the eggs of hadal amphipods (Supplementary Table S3). Most of the PSGs (182/272) possessed zero, one or two positively selected sites (Figure 4A). The distribution of the positively selected sites of the PSGs of amphipods in extreme hadal environments may provide a reference for subsequent studies of biological adaptability. Moreover, GO analysis showed that most of the PSGs were associated with significantly enriched pathways (p < 0.05) involved in cuticle development, DNA repair, cell division, and microtubule regulation (Figures 4B,C). These results demonstrate that HHP in the hadal environment could damage DNA, disintegrate microtubules, and disrupt cell division, suggesting that A. gigantea eggs have adapted to HHP in these aspects.
FIGURE 4. Positively selected genes expressed in the eggs of A. gigantea. (A). The distribution of PSGs by positively selected site number; (B). Bubble plot for the enriched KEGG terms; (C). Chord diagram of GO terms and related genes; (D). A network of positively selected genes in the A. gigantea’s eggs.
The over-represented GO terms associated with the PSGs also included “chitin-based embryonic cuticle biosynthetic process” and “cuticle development” (Figures 4B,C). The cuticle covering the entire surface of A. gigantea has multiple functions, such as forming a rigid exoskeleton for muscle attachment, temporary food storage, and formation of a major environmental barrier to prevent physical injury (Vincent and Wegst, 2004). The cuticle is composed of a composite material consisting of crystalline nanofibres embedded in a matrix of proteins, polyphenols, and water, with small amounts of lipids. Seven genes (i.e., Ci, AKT, UAP, ihog, arrb1, EcR, and grh) involved in cuticle development were positively selected (Figure 4C). Thus, the relationships among these PSGs were investigated (Figure 4D). The core genes involved with significantly enriched GO terms might interact with surrounding genes. However, further studies are needed to elucidate the specific relationships among these genes and the underlying mechanisms.
Determinations of trimethylamine oxide concentrations
TMAO plays important roles in the regulation of osmotic pressure in mammals, hadal zone dwellers, and even microbes (Yancey et al., 2014; Qin et al., 2021; Liu et al., 2022). The TMAO content was much higher in the eggs of A. gigantea than crustaceans inhabiting shallow waters, while there was no significant difference in TMA contents (Figure 5). These results are similar to those of previous studies that TMAO concentrations are significantly higher in hadal amphipods with similar levels of TMA among different tissues (Liu et al., 2022). TMAO is produced by the oxidation of TMA under the action of the flavin monooxygenases (FMOs) (Canyelles et al., 2018), while TMA is mainly formed by the consumption of carnitine choline by gut microorganisms (Wang et al., 2011). Based on the finding that A. gigantea eggs possess relatively high concentrations of TMAO (Figure 5), TMAO might be partially maternally provided, since A. gigantea might not generate TMA in the very early stage of embryogenesis.
FIGURE 5. The concentrations of the TMAO and TMA in the eggs of A. gigantea. (A) and (B), The content of the TMAO and TMA in the eggs of the A. gigantea, P. clarkii, C. quadricarinatus, and M. rosenbergii.
Moreover, FMO1, one of the main flavin monooxygenase genes to transform TMA into TMAO (Canyelles et al., 2018), was identified as a PSG in A. gigantea eggs (Figure 4D and Supplementary Table S3). In the embryo, FMO1 might catalyze maternally acquired TMA to TMAO to enhance survival of A. gigantea eggs in hadal environments.
Discussion
Evolutionary adaptation of adult organisms dwelling in deep-sea environments and hadal zones has been extensively studied in various species, including teleosts (Lan et al., 2018; Wang et al., 2019), bivalve molluscs (Sun et al., 2017), and amphipods (Lan et al., 2017; Li et al., 2021), but rarely in embryos because of the lack of accessibility (Grassle & Sanders, 1973; Perrone et al., 2002; Blankenship-Williams & Levin, 2009; Ritchie et al., 2017). Sampling of gestational amphipods in hadal zones is especially difficult because of limited responses to bait traps and avoided cannibalism, defined as interspecific predation (Fox 1975; Perrone et al., 2002; Blankenship-Williams & Levin, 2009). Some aquatic species are reported to care for eggs in the post-spawning period in deep-sea environments (Seibel et al., 2000; Seibel et al., 2005).
According to previous studies, larval and juvenile amphipods and fish species tend to occupy shallower waters than adult populations (Blankenship et al., 2006; Gerringer, 2019; Weston et al., 2021). Reduced hydrostatic pressure is beneficial to developing instars (Blankenship et al., 2006). However, it is unclear whether brooding amphipods occupy shallower waters in response to pressure or dwell at the same depth as female adults that are not reproductively active (Gerringer, 2019). In this study, brooding amphipod and eggs were collected in the hadal zone, indicating that A. gigantea eggs are tolerant to HHP. However, HHP can rupture the egg membrane, resulting in abnormal mass and extrusion of the cytoplasm (Mestre et al., 2008). Therefore, further studies are warranted to explore the mechanisms underlying the adaptation of A. gigantea eggs to hadal environments.
Reproductive traits of hadal species
As compared to shallow-water crustacean species, A. gigantea eggs are extremely large, consistent with the characteristics of the eggs of other deep-sea and hadal species (Thorson, 1950; Auel, 2004). Higher energetic investment, such as the production of vitellogenin, is needed to support the development of eggs of deep-sea species (Auel, 2004). However, the larger size of female Arctic deep-sea copepods is not correlated to the larger egg size (Auel, 2004). And no correlation between female size and egg size was found among individuals, such as Astacus leptodactylus (Harlioğlu & Türkgülü, 2000). In fact, the egg size of gammarid amphipods is reportedly inversely proportional to temperature during oogenesis. For example, the volume of winter eggs is 60% greater than that of summer eggs (Sheader, 1996). Therefore, the huge egg size of hadal amphipods is possibly an adaptive reproductive trait. The reproductive strategy of K-strategy is to produce few young with high energy content (McLeod et al., 1981). The harsh hadal environment is characterized by fierce competition for resources, less predation pressure, and poor food supply, which results in the K reproductive strategy in A. gigantea. As compared to the high reproduction rate of the r reproductive strategy, the K reproductive strategy has evolved as more efficient due to smaller number of eggs, larger egg size, and longer lifespan (Rushton, 1987). In other words, although the overall number is comparatively decreased, egg viability is improved, which may be an adaptation strategy of hadal species to extreme environments. Furthermore, larger eggs might be a ubiquitous characteristic of animals dwelling in harsh environments, such as Antarctic fish surviving in chronic cold (Shandikov & Faleeva, 1992). Therefore, the morphological characteristics of eggs might underlie successful adaptation of amphipods in hadal zones.
Effects of high hydrostatic pressure on tubulin
HHP has different effects on the assembly of tubulin, actin, myosin II, and cytokeratin in mammalian cells. For example, microtubule is sensitive to high-pressure and excessive hydrostatic pressure can affect the structure of microtubules and cause depolymerization (Zhu et al., 2007; Gao et al., 2018; Yang et al., 2021). HHP is also reported to affect the expression of tubulin (Tahir et al., 1988), regulate development of the cytoskeleton, and induce cytoskeletal deformations (Stramer & Martin, 2005). The sensitivity to high pressure can also change the structure and function of kinesin-microtubule complex. Notably, the GO terms of genes exhibiting accelerated evolution and PSGs were enriched in pathways associated with microtubules and the cytoskeleton (Figures 3, 4), indicating that these genes of A. gigantea eggs may have undergone adaptive changes to cope with the HHP in hadal environments.
Genes exhibiting accelerated evolution related to neuronal cell death
The formation and maintenance of neuronal circuits in the central nervous system is a key process during development (Stephan et al., 2012). Hydrostatic pressure effects on the central nervous system, which might trigger high pressure neurological syndrome (HPNS) (Brauer 1984). HPNS like symptoms such as spasms, convulsions, and tremors have been observed in the shallow-water invertebrates under high pressures, and some genes involved in neuronal DNA damage also showed up-regulations in those shallow-water invertebrates (Morris et al., 2015). Moreover, the neurological impairments, such as convulsions, were also observed in fish under HHP (Sebert 2002). Notably, the roles of adaptive evolution of the nervous system in high hydrostatic pressure deserve better attentions. The neuronal cell death, apoptosis and necrosis can be triggered by excitotoxic damage (Fricker et al., 2018). In vitro, pressure would lead to apoptosis by the induction of oxidative stress (Liu et al., 2007) in B35 and PC12 neuronal cell lines (Agar et al., 2000) and in retinal ganglion cells as well (Agar et al., 2006). In this study, genes exhibiting accelerated evolution, including Psidin/Naa25, Naam, Torsin and TER94, were significantly enriched in “regulation of neuron apoptotic process” and “regulation of neuron death” (Figure 3A). In Drosophila, Psidin/Naa25 is essential for neuron survival and axon targeting (Stephan et al., 2012). Naa25-deficient cells would increase cell apoptosis (Xu et al., 2021). Moreover, the mRNA expression of Drosophila nicotinamidase (D-NAAM) was found to be up-regulated by the low oxidative stress, which might protect the neuronal cells from oxidative stress-induced cell death (Balan et al., 2008). The inhibition of Torsin and overexpression of TER94 were also reported to induce apoptosis in the cultured cells and the Drosophila eyes, respectively (Higashiyama et al., 2002; Chen et al., 2010).
Evolution of cell-cycle regulators contributes to adaptation to hadal environments
The primary features of embryogenesis are vigorous cell division and proliferation. However, mitosis of eukaryotic cells is inhibited at 30 MPa (Seibel et al., 2000; Gao et al., 2018). In echinoids, high pressure affects the success of early embryonic cleavage (Young et al., 1997). On the other hand, as a mutagen, HHP can induce point mutations (Wang et al., 2015). Evolutionary analysis identified PSGs associated with the terms “mitotic cell cycle”, “supervision of G1 DNA damage checkpoint”, “regulation of G1/S transition”, and “initiation of DNA duplication” (Figure 4). These biological processes are regulated by cyclins and various transcription factors, including Ci, FBX O 31, PSME4, and AKT (Figure 4).
As a transcription factor, Ci binds to the promoter of cyclin E (a G1/S cyclin) to upregulate expression and induce DNA replication and transition to the S phase in cells arrested in the G1 phase (Duman-Scheel et al., 2002). Overexpression of Ci is also reported to regulate transcription of cyclin D (Duman-Scheel et al., 2002). The dual-function in cell growth and proliferation suggests that Ci might be a candidate gene to explain the gigantism of A. gigantea.
FBXO31 plays an essential role in G1 arrest and subsequent DNA damage. In mitosis, FBXO31 directly interacts and degrades cyclin D1 via recognition of the phosphorylated Thr286 residue and promotes ubiquitination of Cyclin D1 (Santra et al., 2009). PSME4A/PA200, a 200-kDa nuclear protein, activates the proteasome by promoting degradation of acetylated histones responding to DNA double-strand breaks (Ustrell et al., 2002; Mandemaker et al., 2018). PSME4A/PA200 is also reported to facilitate DNA repair (Ustrell et al., 2002), replication (MacLellan et al., 2015), and genome stability (Blickwedehl et al., 2008). AKT1 was associated with GO terms associated with pleiotropic effects, including “regulation of cell cycle”, “carbohydrate homeostasis”, and “chitin-based embryonic cuticle biosynthetic process” (Lal et al., 2009). In Drosophila melanogaster, AKT phosphorylates CDKN1B to control progression of the cell cycle at the G1 stage (Lal et al., 2009).
Therefore, PSGs (e.g., Ci, FBX O 31, PSME4, and AKT) expressed in the eggs of A. gigantea might regulate mitosis and counteract HHP-induced mutagenesis. So, as an overall evolutionary tendency, cell division is important to the survival of A. gigantea eggs in hadal environments.
Embryonic cuticle development in hadal environments
Multilayered eggshells are produced by follicle cells during oogenesis (Noguerón et al., 2000). In many arthropods, a cuticular egg envelope (blastodermal cuticle in crustaceans; serosal cuticle in insects) is synthesized during early embryogenesis (Lu et al., 2022). During the course of embryogenesis, the blastodermal cuticle is gradually substitutes maternal multilayered eggshell to protect the embryos of arthropods against harsh environmental factors (Vargas et al., 2021).
In the present study, some genes related to embryonic cuticle development, namely cuticle-producing and -organizing factors, were positively selected (Figure 4). Chitin fibers contribute to the stiffness and mechanical resistance of the cuticle (Mrak et al., 2017). UDP-N-acetylhexosamine (UDP-HexNAc) pyrophosphorylase (UAP) is involved in the biosynthesis of chitin, an essential component of the larval cuticle (Stramer & Martin, 2005). Ecdysteroids play an important role in the deposition of chitin in the serosal cuticle. Upregulated expression of ecdysone, the active form of ecdysteroid, triggers the formation of the serosal cuticle (Lagueux et al., 1979; Goltsev et al., 2009).
In addition, the ecdysone receptor (EcR) is also important for development of the embryonic cuticle. EcR is required for maximal expression of the cuticle protein-coding gene EMWCP10 (Wang et al., 2010). Tissue-autonomous EcR is required for concurrent organ morphogenesis in the Drosophila embryo. EcR is a major regulator of tissue development and growth in the marine salmonid ectoparasite Lepeophtheirus salmonis (Copepoda, Caligidae).
The transcription factor grh is also involved in cuticle formation. Drosophila grh regulates epidermal cuticle formation and plays essential roles in wound healing (Mace et al., 2005; Stramer & Martin, 2005). Positive selection of these genes indicates that the embryonic cuticle development of A. gigantea is involved in adaptation to hadal environments.
Trimethylamine oxide in the eggs of hadal amphipods
Proteins tend to lose folding ability under high pressure. Thus, species that dwell in hadal zones need pressure-counteractant to improve protein stability and functional adaptability. TMAO is reportedly important for adaptation of various species to hadal environments (Yancey, 2020). Stress-induced perturbations to the tertiary structure of proteins transcend the protective capacity of TMAO, the heat shock protein could be induced in spiny dogfish shark (Villalobos & Renfro, 2007; MacLellan et al., 2015). As a chemical chaperone, TMAO may be the first line of defense against stress-induced damage and replace functions of other molecular chaperones.
Previous studies have reported that TMAO concentrations are significantly higher in the tissues of hadal species than shallow-water species (Gillett et al., 1997; Samerotte et al., 2007; Liu et al., 2022). Likewise, the results of the present study revealed that TMAO contents were significantly higher in the eggs of hadal amphipods than those of shallow-water crustaceans, while TMA concentrations were similar (Figure 5). TMAO is an important molecule for the survival of hadal species. The high content of TMAO in A. gigantea eggs might be important to ensure normal egg development under HHP conditions. FMO1 is highly expressed in the human embryo (Koukouritaki et al., 2002) and may be involved in maternal interactions (Ishida et al., 2007). Thus, the positively selected FMO1 gene and TMAO play an important function in the eggs of amphipods in hadal environments.
Developmental regulation of amphipod eggs in hadal environments
PSGs combined with genes exhibiting rapid evolution were divided into two groups based on associations with mitosis in the early cleavage stage and cuticle development (Figure 6). Genes involved in the early cleavage stage (i.e., Ci, Fbxo31, and AKT) are mainly related to the GO terms “centrosome duplication”, “regulation of oocyte development”, “spindle”, “microtubule cytoskeleton”, and “regulation of G1/S transition of mitotic cell cycle” (Figure 6). Positive selection of these genes ensures normal cleavage and development of A. gigantea eggs in the early stage.
FIGURE 6. A schematic representation of the adaptive molecular evolution of the A. gigantea embryo.
The PSGs ihog, UAP, ECR, and GRH are mainly related to the molting cycle, chitin-based embryonic cuticle biological processes, and cuticle development and secretion. These PSGs might facilitate adaptation of the A. gigantea cuticle to hadal environments (Figure 6). In Drosophila, dpp is important to embryonic development, formation of the dorsal-ventral axis, and establishment of the early embryonic shape (Auman & Chipman, 2017), and might also regulate related genes involved in the formation of the Drosophila eggshell and appendages (Osterfield et al., 2017). The rapid evolution of dpp indicates that the A. gigantea eggshell has undergone adaptive evolution to cope with harsh hadal environments. Meanwhile, the large amount of TMAO in A. gigantea eggs is believed to maintain the stability and normal function of proteins involved with egg development.
By comparing the results of the eggs with the adults transcriptome analysis (Li et al., 2021), we found the PSGs obtained from the eggs and the adults were both involved in “cell cycle regulation”. However, in the adults whole body, PSGs were closely related to meiotic or mitotic process, such as “meiotic cell cycle”, “female meiotic nuclear division”, “spindle pole”, “mitotic prometaphase” and “condensed chromosome” (Li et al., 2021). Meiosis is a special type of cell division for gametogenesis in gonad. In the eggs, PSGs were enriched in mitotic process, such as “Mitotic DNA damage checkpoint” and “Regulation of G1/S transition of mitotic cell cycle”. The early embryo is characterized by intensive mitotic cell division. Furthermore, the embryonic development of arthropods has sufficient nutrient supply from yolk components (Silveira et al., 2006). This might be the reason why “response to starvation” was enriched in whole body but not enriched in PSGs in the eggs. Since cuticle formation happens during embryonic development of Arthropoda (Ostrowski et al., 2002), the eggs are excellent specimen to investigate the genetic control of cuticle development. The “cuticle development” and “chitin-based embryonic cuticle biosynthetic process” were enriched in the eggs, which certainly can’t be unfolded by the adults study. Collectively, these findings deepen the current understanding of embryogenetic adaptations of hadal amphipods to HHP and provide valuable genetic resources for adaptation to hadal environments.
Data availability statement
All sequencing data associated with this project were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive database (BioProject Accession Number: PRJNA841239).
Ethics statement
Experimental protocols involved dead animals in this study.
Author contributions
QX conceived the experiments and led the whole project. WL and FW analyzed the data. SJ extracted the RNA and performed RNA sequencing. QL measured the TMAO and TMA concentrations of the eggs. BP designed the lander vehicle for sample collection. WL and QX wrote the manuscript. And QX contributed to edits to the manuscript.
Acknowledgments
We would like to thank Shanghai Rainbowfish Ocean Technology Co., Ltd. For sample collection. We also thank Professor Weicheng Cui’s and Jiasong Fang’s research group members and other people for sample collection. This work was supported in part by the Funding Project of the National Key Research and Development Program of China (2018YFC0310600, 2018YFD0900601), the National Natural Science Foundation of China (Grant Nos. 31772826), and the major scientific innovation project from Shanghai Committee of Education (2017–01–07–00–10-E00060).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.987409/full#supplementary-material
References
Agar, A., Li, S., Agarwal, N., Coroneo, M. T., and Hill, M. A. (2006). Retinal ganglion cell line apoptosis induced by hydrostatic pressure. Brain Res. 1086 (1), 191–200. doi:10.1016/j.brainres.2006.02.061
Agar, A., Yip, S. S., Hill, M. A., and Coroneo, M. T. (2000). Pressure related apoptosis in neuronal cell lines. J. Neurosci. Res. 60 (4), 495–503. doi:10.1002/(sici)1097-4547(20000515)60:4<495:aid-jnr8>3.0.co;2-s
Ángeles-González, L. E., Lima, F. D., Caamal-Monsreal, C., Díaz, F., and Rosas, C. (2020). Exploring the effects of warming seas by using the optimal and pejus temperatures of the embryo of three Octopoda species in the Gulf of Mexico. J. Therm. Biol. 94, 102753. doi:10.1016/j.jtherbio.2020.102753
Arambourou, H., Chaumot, A., Vulliet, E., Daniele, G., Delorme, N., Abbaci, K., et al. (2017). Phenotypic defects in newborn Gammarus fossarum (Amphipoda) following embryonic exposure to fenoxycarb. Ecotoxicol. Environ. Saf. 144, 193–199. doi:10.1016/j.ecoenv.2017.06.017
Auel, H. (2004). Egg size and reproductive adaptations among Arctic deep-sea copepods (Calanoida, Paraeuchaeta). Helgol Mar. Res. 58 (3), 147–153. doi:10.1007/s10152-004-0179-1
Auman, T., and Chipman, A. D. (2017). The evolution of gene regulatory networks that define arthropod body plans. Integr. Comp. Biol. 57 (3), 523–532. doi:10.1093/icb/icx035
Balan, V., Miller, G. S., Kaplun, L., Balan, K., Chong, Z. Z., and Li, F. (2008). Life span extension and neuronal cell protection by Drosophila nicotinamidase. J. Biol. Chem. 283 (41), 27810–27819. doi:10.1074/jbc.M804681200
Bar-On, Y. M., Phillips, R., and Milo, R. (2018). The biomass distribution on Earth. Proc. Natl. Acad. Sci. U. S. A. 115 (25), 6506–6511. doi:10.1073/pnas.1711842115
Barnard, J. L., and Ingram, C. L. (1986). The supergiant amphipod Alicella gigantea chevreux from the north pacific gyre. J. Crustacean Biol. 6 (4), 825–839. doi:10.2307/1548395
Blankenship, L. E., Yayanos, A. A., Cadien, D. B., and Levin, L. A. (2006). Vertical zonation patterns of scavenging amphipods from the Hadal zone of the Tonga and Kermadec Trenches. Deep Sea Res. Part I Oceanogr. Res. Pap. 53 (1), 48–61. doi:10.1016/j.dsr.2005.09.006
Blankenship-Williams, L. E., and Levin, L. A. (2009). Living deep: A synopsis of hadal trench ecology. Mar. Technol. Soc. J. 43 (5), 137–143. doi:10.4031/mtsj.43.5.23
Blickwedehl, J., Agarwal, M., Seong, C., Pandita, R. K., Melendy, T., Sung, P., et al. (2008). Role for proteasome activator PA200 and postglutamyl proteasome activity in genomic stability. Proc. Natl. Acad. Sci. U. S. A. 105 (42), 16165–16170. doi:10.1073/pnas.0803145105
Brauer, R. W. (1984). Hydrostatic pressure effects on the central nervous system: Perspectives and outlook. Philos. Trans. R. Soc. Lond B Biol. Sci. 304 (1118), 17–30. doi:10.1098/rstb.1984.0005
Brooks, N. J. (2014). Pressure effects on lipids and bio-membrane assemblies. IUCrJ 1 (6), 470–477. doi:10.1107/S2052252514019551
Canyelles, M., Tondo, M., Cedó, L., Farràs, M., Escolà-Gil, J. C., and Blanco-Vaca, F. (2018). Trimethylamine N-oxide: A link among diet, gut microbiota, gene regulation of liver and intestine cholesterol homeostasis and hdl function. Int. J. Mol. Sci. 19 (10), 3228. doi:10.3390/ijms19103228
Chan, J., Geng, D., Pan, B., Zhang, Q., and Xu, Q. (2022). Gut microbial divergence between three hadal amphipod species from the isolated hadal trenches. Microb. Ecol. 84 (2), 627–637. doi:10.1007/s00248-021-01851-3
Chan, J., Geng, D., Pan, B., Zhang, Q., and Xu, Q. (2021). Metagenomic insights into the structure and function of intestinal microbiota of the hadal amphipods. Front. Microbiol. 12, 668989. doi:10.3389/fmicb.2021.668989
Chan, J., Pan, B., Geng, D., Zhang, Q., Zhang, S., Guo, J., et al. (2020). Genetic diversity and population structure analysis of three deep-sea amphipod species from geographically isolated hadal trenches in the Pacific ocean. Biochem. Genet. 58 (1), 157–170. doi:10.1007/s10528-019-09935-z
Chen, X. P., Hu, X. H., Wu, S. H., Zhang, Y. W., Xiao, B., and Shang, H. F. (2010). RNA interference-mediated inhibition of wild-type Torsin A expression increases apoptosis caused by oxidative stress in cultured cells. Neurochem. Res. 35 (8), 1214–1223. doi:10.1007/s11064-010-0177-4
Colbourne, J. K., Pfrender, M. E., Gilbert, D., Thomas, W. K., Tucker, A., Oakley, T. H., et al. (2011). The ecoresponsive genome of Daphnia pulex. Science 331 (6017), 555–561. doi:10.1126/science.1197761
Csardi, G., and Nepusz, T. (2006). The igraph software package for complex network research. InterJournal, complex Syst. 1695 (5), 1–9.
Duman-Scheel, M., Weng, L., Xin, S., and Du, W. (2002). Hedgehog regulates cell growth and proliferation by inducing Cyclin D and Cyclin E. Nature 417 (6886), 299–304. doi:10.1038/417299a
Fox, L. R. (1975). Cannibalism in natural populations. Annu. Rev. Ecol. Syst. 6 (1), 87–106. doi:10.1146/annurev.es.06.110175.000511
Fricker, M., Tolkovsky, A. M., Borutaite, V., Coleman, M., and Brown, G. C. (2018). Neuronal cell death. Physiol. Rev. 98 (2), 813–880. doi:10.1152/physrev.00011.2017
Gao, M., Berghaus, M., Mobitz, S., Schuabb, V., Erwin, N., Herzog, M., et al. (2018). On the origin of microtubules' high-pressure sensitivity. Biophysical J. 114 (5), 1080–1090. doi:10.1016/j.bpj.2018.01.021
Gerringer, M. E. (2019). On the success of the hadal snailfishes. Integr. Org. Biol. 1 (1), obz004. doi:10.1093/iob/obz004
Gillett, M. B., Suko, J. R., Santoso, F. O., and Yancey, P. H. (1997). Elevated levels of trimethylamine oxide in muscles of deep-sea gadiform teleosts: A high-pressure adaptation. J. Exp. Zoology 279 (4), 386–391. doi:10.1002/(sici)1097-010x(19971101)279:4<386:aid-jez8>3.0.co;2-k
Glover, K. A., Harvey, A. C., Hansen, T. J., Fjelldal, P. G., Besnier, F. N., Bos, J. B., et al. (2020). Chromosome aberrations in pressure-induced triploid Atlantic salmon. BMC Genet. 21 (1), 59. doi:10.1186/s12863-020-00864-0
Goltsev, Y., Rezende, G. L., Vranizan, K., Lanzaro, G., Valle, D., and Levine, M. (2009). Developmental and evolutionary basis for drought tolerance of the Anopheles gambiae embryo. Dev. Biol. 330 (2), 462–470. doi:10.1016/j.ydbio.2009.02.038
Grassle, J. F., and Sanders, H. L. (1973). Life histories and the role of disturbance. Deep Sea Res. Oceanogr. Abstr. 20 (7), 643–659.
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8 (8), 1494–1512. doi:10.1038/nprot.2013.084
Harlioğlu, M. M., and Türkgülü, I. (2000). The relationship between egg size and female size infreshwater crayfish, Astacus leptodactylu. Aquac. Int. 8 (1), 95–98.
Harrison, C. S., Hida, T. S., and Seki, M. P. (1983). Hawaiin seabird feeding ecology. Wildl. Monogr. 85, 1–71.
Higashiyama, H., Hirose, F., Yamaguchi, M., Inoue, Y., Fujikake, N., Matsukage, A., et al. (2002). Identification of ter94, Drosophila VCP, as a modulator of polyglutamine-induced neurodegeneration. Cell Death Differ. 9 (3), 264–273. doi:10.1038/sj/cdd/440095510.1038/sj.cdd.4400955
Ishida, M., Ohashi, S., Kizaki, Y., Naito, J., Horiguchi, K., and Harigaya, T. (2007). Expression profiling of mouse placental lactogen II and its correlative genes using a cDNA microarray analysis in the developmental mouse placenta. J. reproduction Dev. 53 (1), 69–76. doi:10.1262/jrd.18002
Ishii, A., Sato, T., Wachi, M., Nagai, K., and Kato, C. (2004). Effects of high hydrostatic pressure on bacterial cytoskeleton FtsZ polymers in vivo and in vitro. Microbiol. Read. 150 (6), 1965–1972. doi:10.1099/mic.0.26962-0
Jamieson, A. J., Fujii, T., Mayor, D. J., Solan, M., and Priede, I. G. (2010). Hadal trenches: The ecology of the deepest places on earth. Trends Ecol. Evol. 25 (3), 190–197. doi:10.1016/j.tree.2009.09.009
Jamieson, A. J., Lacey, N. C., Lörz, A. N., Rowden, A. A., and Piertney, S. B. (2013). The supergiant amphipod Alicella gigantea (Crustacea: Alicellidae) from hadal depths in the kermadec trench, SW Pacific Ocean. Deep Sea Res. Part II Top. Stud. Oceanogr. 92, 107–113. doi:10.1016/j.dsr2.2012.12.002
Kao, D., Lai, A. G., Stamataki, E., Rosic, S., Konstantinides, N., Jarvis, E., et al. (2016). The genome of the crustacean Parhyale hawaiensis, a model for animal development, regeneration, immunity and lignocellulose digestion. Elife 5, e20062. doi:10.7554/eLife.20062
Koukouritaki, S. B., Simpson, P., Yeung, C. K., Rettie, A. E., and Hines, R. N. (2002). Human hepatic flavin-containing monooxygenases 1 (FMO1) and 3 (FMO3) developmental expression. Pediatr. Res. 51 (2), 236–243. doi:10.1203/00006450-200202000-00018
Lagueux, M., Hetru, C., Goltzene, F., Kappler, C., and Hoffmann, J. (1979). Ecdysone titre and metabolism in relation to cuticulogenesis in embryos of Locusta migratoria. J. Insect Physiology 25 (9), 709–723. doi:10.1016/0022-1910(79)90123-9
Lal, M. A., Bae, D., Camilli, T. C., Patierno, S. R., and Ceryak, S. (2009). AKT1 mediates bypass of the G1/S checkpoint after genotoxic stress in normal human cells. Cell Cycle 8 (10), 1589–1602. doi:10.4161/cc.8.10.8547
Lan, Y., Sun, J., Tian, R., Bartlett, D. H., Li, R., Wong, Y. H., et al. (2017). Molecular adaptation in the world's deepest-living animal: Insights from transcriptome sequencing of the hadal amphipod Hirondellea gigas. Mol. Ecol. 26 (14), 3732–3743. doi:10.1111/mec.14149
Lan, Y., Sun, J., Xu, T., Chen, C., Tian, R., Qiu, J. W., et al. (2018). De novo transcriptome assembly and positive selection analysis of an individual deep-sea fish. BMC Genomics 19 (1), 394. doi:10.1186/s12864-018-4720-z
Li, J. Y., Song, Z. L., Yan, G. Y., and He, L. S. (2019). The complete mitochondrial genome of the largest amphipod, Alicella gigantea: Insight into its phylogenetic relationships and deep sea adaptive characters. Int. J. Biol. Macromol. 141, 570–577. doi:10.1016/j.ijbiomac.2019.09.050
Li, W., Wang, F., Jiang, S., Pan, B., Chan, J., and Xu, Q. (2021). The adaptive evolution and gigantism mechanisms of the hadal “supergiant” amphipod Alicella gigantea. Front. Mar. Sci. 8, 743663. doi:10.3389/fmars.2021.743663
Lin, Z., Zhu, X., Zhang, T., You, F., Wu, Z., and Cao, Y. (2016). Effects of hydrostatic pressure on microtubule organization and nucleus changes in gynogenetically activated eggs of olive flounder (Paralichthys olivaceus). Theriogenology 85 (9), 1610–1624. doi:10.1016/j.theriogenology.2016.01.020
Liu, Q., Jiang, S., Li, W., Pan, B., and Xu, Q. (2022). Trimethylamine N-oxide (TMAO) and Trimethylamine (TMA) determinations of two hadal amphipods. J. Mar. Sci. Eng. 10 (4), 454. doi:10.3390/jmse10040454
Liu, Q., Ju, W. K., Crowston, J. G., Xie, F., Perry, G., Smith, M. A., et al. (2007). Oxidative stress is an early event in hydrostatic pressure induced retinal ganglion cell damage. Invest. Ophthalmol. Vis. Sci. 48 (10), 4580–4589. doi:10.1167/iovs.07-0170
Lu, J. B., Guo, J. S., Chen, X., Cheng, C., Luo, X. M., Zhang, X. Y., et al. (2022). Chitin synthase 1 and five cuticle protein genes are involved in serosal cuticle formation during early embryogenesis to enhance eggshells in Nilaparvata lugens. Insect Sci. 29 (2), 363–378. doi:10.1111/1744-7917.12937organizing factors, were positively selected
Mace, K. A., Pearson, J. C., and McGinnis, W. (2005). An epidermal barrier wound repair pathway in Drosophila is mediated by grainy head. Science 308 (5720), 381–385. doi:10.1126/science.1107573
MacLellan, R. J., Tunnah, L., Barnett, D., Wright, P. A., MacCormack, T., and Currie, S. (2015). Chaperone roles for TMAO and HSP70 during hyposmotic stress in the spiny dogfish shark (Squalus acanthias). J. Comp. Physiology B 185 (7), 729–740. doi:10.1007/s00360-015-0916-6
Mandemaker, I. K., Geijer, M. E., Kik, I., Bezstarosti, K., Rijkers, E., Raams, A., et al. (2018). DNA damage-induced replication stress results in PA200-proteasome-mediated degradation of acetylated histones. EMBO Rep. 19 (10), e45566. doi:10.15252/embr.201745566
Margaritis, L. H. (1985). Comparative study of the eggshell of the fruit flies Dacus oleae and Ceratitis capitata (Diptera: Trypetidae). Can. J. Zoology 63 (9), 2194–2206. doi:10.1139/z85-324
McLeod, M., Hornbach, D. J., Guttman, S., Way, C., and Burky, A. J. (1981). Environmental heterogeneity, genetic polymorphism, and reproductive strategies. Am. Nat. 118 (1), 129–134. doi:10.1086/283808
Mestre, N. C., Thatje, S., and Tyler, P. A. (2008). The ocean is not deep enough: Pressure tolerances during early ontogeny of the blue mussel Mytilus edulis. Proc. R. Soc. B Biol. Sci. 276 (1657), 717–726. doi:10.1098/rspb.2008.1376
Morris, J. P., Thatje, S., Ravaux, J., Shillito, B., Fernando, D., and Hauton, C. (2015). Acute combined pressure and temperature exposures on a shallow-water crustacean: Novel insights into the stress response and high pressure neurological syndrome. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 181, 9–17. doi:10.1016/j.cbpa.2014.10.028
Mouzaki, D. G., and Margaritis, L. H. (1991). The eggshell of the cherry fly Rhagoletis cerasi. Tissue Cell 23 (5), 745–754. doi:10.1016/0040-8166(91)90027-q
Mrak, P., Bogataj, U., Strus, J., and Znidarsic, N. (2017). Cuticle morphogenesis in crustacean embryonic and postembryonic stages. Arthropod Struct. Dev. 46 (1), 77–95. doi:10.1016/j.asd.2016.11.001
Noguerón, M. I., Mauzy-Melitz, D., and Waring, G. L. (2000). Drosophila dec-1 eggshell proteins are differentially distributed via a multistep extracellular processing and localization pathway. Dev. Biol. 225 (2), 459–470. doi:10.1006/dbio.2000.9805
Osterfield, M., Berg, C. A., and Shvartsman, S. Y. (2017). Epithelial patterning, morphogenesis, and evolution: Drosophila eggshell as a model. Dev. Cell 41 (4), 337–348. doi:10.1016/j.devcel.2017.02.018
Ostrowski, S., Dierick, H. A., and Bejsovec, A. (2002). Genetic control of cuticle formation during embryonic development of Drosophila melanogaster. Genetics 161 (1), 171–182. doi:10.1093/genetics/161.1.171
Perrone, F. M., Dell'Anno, A., Danovaro, R., Croce, N. D., and Thurston, M. H. (2002). Population biology of hirondellea sp. nov. (Amphipoda: Gammaridea: Lysianassoidea) from the atacama trench (south-east Pacific ocean). J. Mar. Biol. Assoc. U. K. 82 (3), 419–425. doi:10.1017/s0025315402005672
Poynton, H. C., Hasenbein, S., Benoit, J. B., Sepulveda, M. S., Poelchau, M. F., Hughes, D. S. T., et al. (2018). The toxicogenome of Hyalella azteca: A model for sediment ecotoxicology and evolutionary toxicology. Environ. Sci. Technol. 52 (10), 6009–6022. doi:10.1021/acs.est.8b00837
Qin, Q. L., Wang, Z. B., Su, H. N., Chen, X. L., Miao, J., Wang, X. J., et al. (2021). Oxidation of trimethylamine to trimethylamine N-oxide facilitates high hydrostatic pressure tolerance in a generalist bacterial lineage. Sci. Adv. 7 (13), eabf9941. doi:10.1126/sciadv.abf9941
Rabet, N., Godinho, L. B., Montero, D., and Lacau, S. (2012). Exploration of the egg shell structure of three neotropical eulimnadiaspecies: A new insight into genus taxonomy (Crustacea: Branchiopoda: Spinicaudata). Stud. Neotropical Fauna Environ. 47 (3), 221–226. doi:10.1080/01650521.2012.711967
Ritchie, H., Jamieson, A. J., and Piertney, S. B. (2017). Genome size variation in deep-sea amphipods. R. Soc. open Sci. 4 (9), 170862. doi:10.1098/rsos.17086210.6084/m9
Ritchie, H., Jamieson, A. J., and Piertney, S. B. (2015). Phylogenetic relationships among hadal amphipods of the Superfamily Lysianassoidea: Implications for taxonomy and biogeography. Deep Sea Res. Part I Oceanogr. Res. Pap. 105, 119–131. doi:10.1016/j.dsr.2015.08.014
Rushton, J. P. (1987). Toward a theory of human multiple birthing: Sociobiology and r/K reproductive strategies. Acta Genet. medicae gemellologiae twin Res. 36 (3), 289–296. doi:10.1017/s0001566000006048
Samerotte, A. L., Drazen, J. C., Brand, G. L., Seibel, B. A., and Yancey, P. H. (2007). Correlation of trimethylamine oxide and habitat depth within and among species of teleost fish: An analysis of causation. Physiological Biochem. zoology 80 (2), 197–208. doi:10.1086/510566
Santra, M. K., Wajapeyee, N., and Green, M. R. (2009). F-box protein FBXO31 mediates cyclin D1 degradation to induce G1 arrest after DNA damage. Nature 459 (7247), 722–725. doi:10.1038/nature08011
Scholz, S., Fischer, S., Gundel, U., Kuster, E., Luckenbach, T., and Voelker, D. (2008). The zebrafish embryo model in environmental risk assessment--applications beyond acute toxicity testing. Environ. Sci. Pollut. Res. Int. 15 (5), 394–404. doi:10.1007/s11356-008-0018-z
Sebert, P. (2002). Fish at high pressure: A hundred year history. Comp. Biochem. Physiology Part A Mol. Integr. Physiology 131 (3), 575–585. doi:10.1016/s1095-6433(01)00509-8
Seibel, B. A., Hochberg, F. G., and Carlini, D. B. (2000). Life history of gonatus onyx (cephalopoda: Teuthoidea): deep-sea spawning and post-spawning egg care. Mar. Biol. 137 (3), 519–526. doi:10.1007/s002270000359
Seibel, B. A., Robison, B. H., and Haddock, S. H. (2005). Post-spawning egg care by a squid. Nature 438 (7070), 929. doi:10.1038/438929a
Serra, G. G. (2003). Some observations on the life history of the freshwater amphipod Echinogammarus longisetosus Pinkster, 1973 (Gammaridae) from Catalonia (Spain, N Iberian peninsula). Animal Biodivers. Conservation 26 (1), 31–39.
Shandikov, G. A., and Faleeva, T. I. (1992). Features of gametogenesis and sexual cycles of six notothenioid fishes from East Antarctica. Polar Biol. 11 (8), 615–621. doi:10.1007/bf00237956
Sheader, M. (1996). Factors influencing egg size in the gammarid amphipod Gammarus insensibilis. Mar. Biol. 124 (4), 519–526. doi:10.1007/bf00351033
Silveira, A. B., Castro-Santos, J., Senna, R., Logullo, C., Fialho, E., and Silva-Neto, M. A. (2006). Tick vitellin is dephosphorylated by a protein tyrosine phosphatase during egg development: Effect of dephosphorylation on VT proteolysis. Insect Biochem. Mol. Biol. 36 (3), 200–209. doi:10.1016/j.ibmb.2006.01.015
Stephan, D., Sanchez-Soriano, N., Loschek, L. F., Gerhards, R., Gutmann, S., Storchova, Z., et al. (2012). Drosophila Psidin regulates olfactory neuron number and axon targeting through two distinct molecular mechanisms. J. Neurosci. 32 (46), 16080–16094. doi:10.1523/JNEUROSCI.3116-12.2012
Stramer, B., and Martin, P. (2005). Cell biology: Master regulators of sealing and healing. Curr. Biol. 15 (11), R425–R427. doi:10.1016/j.cub.2005.05.034
Sun, J., Zhang, Y., Xu, T., Zhang, Y., Mu, H., Zhang, Y., et al. (2017). Adaptation to deep-sea chemosynthetic environments as revealed by mussel genomes. Nat. Ecol. Evol. 1 (5), 121. doi:10.1038/s41559-017-0121
Tahir, S. K., Symington, A. L., and Zimmerman, A. M. (1988). Pressure sensitivity of tubulin expression in Tetrahymena. Cell Biol. Int. Rep. 12 (12), 1005–1009. doi:10.1016/0309-1651(88)90025-2
Thorson, G. (1950). Reproductive and larval ecology of marine bottom invertebrates. Biol. Rev. 25 (1), 1–45. doi:10.1111/j.1469-185x.1950.tb00585.x
Tommasini, S., Scanabissisabelli, F., and Trentini, M. (1989). Scanning electron microscope study of eggshell development in Triops cancriformis (Bosc)(Crustacea, Notostraca). Vie Milieu/Life Environ. 39 (1), 29–32.
Ustrell, V., Hoffman, L., Pratt, G., and Rechsteiner, M. (2002). PA200, a nuclear proteasome activator involved in DNA repair. EMBO J. 21 (13), 3516–3525. doi:10.1093/emboj/cdf333
Vargas, H. C. M., Panfilio, K. A., Roelofs, D., and Rezende, G. L. (2021). Increase in egg resistance to desiccation in springtails correlates with blastodermal cuticle formation: Eco-evolutionary implications for insect terrestrialization. J. Exp. Zoology. Part B, Mol. Dev. Evol. 336 (8), 606–619. doi:10.1002/jez.b.22979
Villalobos, A. R., and Renfro, J. L. (2007). Trimethylamine oxide suppresses stress-induced alteration of organic anion transport in choroid plexus. J. Exp. Biol. 210 (3), 541–552. doi:10.1242/jeb.02681
Villalobos, F. B., Tyler, P. A., and Young, C. M. (2006). Temperature and pressure tolerance of embryos and larvae of the atlantic seastars Asterias rubens and Marthasterias glacialis (echinodermata: Asteroidea): Potential for deep-sea invasion. Mar. Ecol. Prog. Ser. 314, 109–117.
Vincent, J. F., and Wegst, U. G. (2004). Design and mechanical properties of insect cuticle. Arthropod Struct. Dev. 33 (3), 187–199. doi:10.1016/j.asd.2004.05.006
Walter, W., Sanchez-Cabo, F., and Ricote, M. (2015). GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics 31 (17), 2912–2914. doi:10.1093/bioinformatics/btv300
Wang, H. B., Moriyama, M., Iwanaga, M., and Kawasaki, H. (2010). Ecdysone directly and indirectly regulates a cuticle protein gene, BMWCP10, in the wing disc of Bombyx mori. Insect Biochem. Mol. Biol. 40 (6), 453–459. doi:10.1016/j.ibmb.2010.04.004
Wang, H., Wang, K., Xiao, G., Ma, J., Wang, B., Shen, S., et al. (2015). Molecular mechanisms for high hydrostatic pressure-induced wing mutagenesis in. Drosophila melanogaster. Sci. Rep. 5, 14965. doi:10.1038/srep14965
Wang, J., Itgen, M. W., Wang, H., Gong, Y., Jiang, J., Li, J., et al. (2021). Gigantic genomes provide empirical tests of transposable element dynamics models. Genomics, Proteomics Bioinforma. 19 (1), 123–139. doi:10.1016/j.gpb.2020.11.005
Wang, K., Shen, Y., Yang, Y., Gan, X., Liu, G., Hu, K., et al. (2019). Morphology and genome of a snailfish from the Mariana Trench provide insights into deep-sea adaptation. Nat. Ecol. Evol. 3 (5), 823–833. doi:10.1038/s41559-019-0864-8
Wang, Z., Klipfell, E., Bennett, B. J., Koeth, R., Levison, B. S., Dugar, B., et al. (2011). Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature 472 (7341), 57–63. doi:10.1038/nature09922
Weston, J. N. J., Espinosa-Leal, L., Wainwright, J. A., Stewart, E. C. D., Gonzalez, C. E., Linley, T. D., et al. (2021). Eurythenes atacamensis sp. nov. (Crustacea: Amphipoda) exhibits ontogenetic vertical stratification across abyssal and hadal depths in the atacama trench, eastern south Pacific ocean. Mar. Biodivers. a J. Senckenberg Res. Inst. 51 (3), 51. doi:10.1007/s12526-021-01182-z
Xu, J., Li, Z., Zuo, X., Li, G., Zhang, X., Zhang, B., et al. (2021). Knockdown of NAA25 suppresses breast cancer progression by regulating apoptosis and cell cycle. Front. Oncol. 11, 755267. doi:10.3389/fonc.2021.755267
Yamamoto, E. (1999). Studies on sex-manipulation and production of cloned populations in hirame, Paralichthys olivaceus (Temminck et Schlegel). Aquaculture 173 (1-4), 235–246. doi:10.1016/s0044-8486(98)00448-7
Yancey, P. H. (2020). Cellular responses in marine animals to hydrostatic pressure. J. Exp. Zoology. Part A, Ecol. Integr. Physiology 333 (6), 398–420. doi:10.1002/jez.2354
Yancey, P. H., Gerringer, M. E., Drazen, J. C., Rowden, A. A., and Jamieson, A. (2014). Marine fish may be biochemically constrained from inhabiting the deepest ocean depths. Proc. Natl. Acad. Sci. U. S. A. 111 (12), 4461–4465. doi:10.1073/pnas.1322003111
Yang, G. N., Ahangar, P., Strudwick, X. L., Kopecki, Z., and Cowin, A. J. (2021). Overexpression of flii during murine embryonic development increases symmetrical division of epidermal progenitor cells. Intertional J. Mol. Sci. 22 (15), 8235. doi:10.3390/ijms22158235
Yang, Z. H. (2007). PAML 4:Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24 (8), 1586–1591. doi:10.1093/molbev/msm088
Young, C. M., Tyler, P. A., and Fenaux, L. (1997). Potential for deep sea invasion by mediterranean shallow water echinoids: Pressure and temperature as stage-specific dispersal barriers. Mar. Ecol. Prog. Ser. 154, 197–209. doi:10.3354/meps154197
Zhang, X., Yuan, J., Sun, Y., Li, S., Gao, Y., Yu, Y., et al. (2019). Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nat. Commun. 10 (1), 356. doi:10.1038/s41467-018-08197-4
Zhu, L., Geng, D., Pan, B., Li, W., Jiang, S., and Xu, Q. (2022). Trace elemental analysis of the exoskeleton, leg muscle, and gut of three hadal amphipods. Biol. Trace Elem. Res. 200 (3), 1395–1407. doi:10.1007/s12011-021-02728-9
Zhu, X., Lin, Z., Wu, Z., Li, J., and You, F. (2017). Effect of initiation time of hydrostatic pressure shock on chromosome set doubling of tetraploidization in turbot Scophthalmus maximus. Mar. Biotechnol. (New York, N.Y.) 19 (5), 528–540. doi:10.1007/s10126-017-9771-7
Keywords: hadal trench, amphipod, eggs, high hydrostatic pressure, molecular evolution
Citation: Li W, Wang F, Jiang S, Pan B, Liu Q and Xu Q (2022) Morphological and molecular evolution of hadal amphipod’s eggs provides insights into embryogenesis under high hydrostatic pressure. Front. Cell Dev. Biol. 10:987409. doi: 10.3389/fcell.2022.987409
Received: 06 July 2022; Accepted: 26 August 2022;
Published: 12 September 2022.
Edited by:
Juan Jesus Tena, Spanish National Research Council (CSIC), SpainReviewed by:
José María Santos-Pereira, Andalusian Center for Development Biology, Spanish National Research Council (CSIC), SpainCarlos Santamaria, University of Tampa, United States
Copyright © 2022 Li, Wang, Jiang, Pan, Liu and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qianghua Xu, qhxu@shou.edu.cn