Genome size variation in deep-sea amphipods

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Phylogenomic analyses reveal a single deep-water colonisation in Patellogastropoda

Patellogastropoda, the true limpets, is a major group of gastropods widely distributed in marine habitats from the intertidal to deep sea. Though important for understanding their evolutionary radiation, the phylogenetic relationships among the patellogastropod families have always been challenging to reconstruct, with contradictory results likely due to insufficient sampling. Here, we obtained mitogenomic and phylogenomic data (transcriptomic or genomic) from six species representing the three predominantly deep-water patellogastropod families: Lepetidae, Neolepetopsidae, and Pectinodontidae. By using various phylogenetic methods, we show that mitogenome phylogeny recovers monophyly of eight families in most of the trees, though the relationships among families remain contentious. Meanwhile, a more robust family-level topology consistent with morphology was achieved by phylogenomics. This also reveals that these mainly deep-water families are monophyletic, suggesting a single colonisation of the deep water around the Jurassic. We also found a lack of significant correlation between genome size and habitat depth, despite some deep-water species exhibiting larger genome sizes. Our phylogenomic tree provides a stable phylogenetic backbone for Patellogastropoda that includes seven of the nine recognized families and paves the way for future evolutionary analyses in this major group of molluscs.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles

BACKGROUND

Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles.

RESULTS

The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species.

CONCLUSIONS

Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

Genomic Analysis of a Scale Worm Provides Insights into Its Adaptation to Deep-Sea Hydrothermal Vents

Deep-sea polynoid scale worms endemic to hydrothermal vents have evolved an adaptive strategy to the chronically hypoxic environment, but its underlying molecular mechanisms remain elusive. Here, we assembled a chromosome-scale genome of the vent-endemic scale worm Branchipolynoe longqiensis (the first annotated genome in the subclass Errantia) and annotated two shallow-water polynoid genomes, aiming to elucidate the adaptive mechanisms. We present a genome-wide molecular phylogeny of Annelida which calls for extensive taxonomy revision by including more genomes from key lineages. The B. longqiensis genome with a genome size of 1.86 Gb and 18 pseudochromosomes is larger than the genomes of two shallow-water polynoids, possibly due to the expansion of various transposable elements (TEs) and transposons. We revealed two interchromosomal rearrangements in B. longqiensis when compared with the two shallow-water polynoid genomes. The intron elongation and interchromosomal rearrangement can influence a number of biological processes, such as vesicle transport, microtubules, and transcription factors. Furthermore, the expansion of cytoskeleton-related gene families may favor the cell structure maintenance of B. longqiensis in the deep ocean. The expansion of synaptic vesicle exocytosis genes has possibly contributed to the unique complex structure of the nerve system in B. longqiensis. Finally, we uncovered an expansion of single-domain hemoglobin and a unique formation of tetra-domain hemoglobin via tandem duplications, which may be related to the adaptation to a hypoxic environment.

Large effective population size masks population genetic structure in Hirondellea amphipods within the deepest marine ecosystem, the Mariana Trench

The examination of genetic structure in the deep-ocean hadal zone has focused on divergence between tectonic trenches to understand how environment and geography may drive species divergence and promote endemism. There has been little attempt to examine localized genetic structure within trenches, partly because of logistical challenges associated with sampling at an appropriate scale, and the large effective population sizes of species that can be sampled adequately may mask underlying genetic structure. Here we examine genetic structure in the superabundant amphipod Hirondellea gigas in the Mariana Trench at depths of 8126-10,545 m. RAD sequencing was used to identify 3182 loci containing 43,408 single nucleotide polymorphisms (SNPs) across individuals after stringent pruning of loci to prevent paralogous multicopy genomic regions being erroneously merged. Principal components analysis of SNP genotypes resolved no genetic structure between sampling locations, consistent with a signature of panmixia. However, discriminant analysis of principal components identified divergence between all sites driven by 301 outlier SNPs in 169 loci and significantly associated with latitude and depth. Functional annotation of loci identified differences between singleton loci used in analysis and paralogous loci pruned from the data set and also between outlier and nonoutlier loci, all consistent with hypotheses explaining the role of transposable elements driving genome dynamics. This study challenges the traditional perspective that highly abundant amphipods within a trench form a single panmictic population. We discuss the findings in relation to eco-evolutionary and ontogenetic processes operating in the deep sea, and highlight key challenges associated with population genetic analysis in nonmodel systems with inherent large effective population sizes and genomes.

Barriers to gene flow in the deepest ocean ecosystems: Evidence from global population genomics of a cosmopolitan amphipod

The deepest marine ecosystem, the hadal zone, hosts endemic biodiversity resulting from geographic isolation and environmental selection pressures. However, the pan-ocean distribution of some fauna challenges the concept that the hadal zone is a series of isolated island-like habitats. Whether this remains true at the population genomic level is untested. We investigated phylogeographic patterns of the amphipod, Bathycallisoma schellenbergi, from 12 hadal features across the Pacific, Atlantic, Indian, and Southern oceans and analyzed genome-wide single-nucleotide polymorphism markers and two mitochondrial regions. Despite a cosmopolitan distribution, populations were highly restricted to individual features with only limited gene flow between topographically connected features. This lack of connectivity suggests that populations are on separate evolutionary trajectories, with evidence of potential cryptic speciation at the Atacama Trench. Together, this global study demonstrates that the shallower ocean floor separating hadal features poses strong barriers to dispersal, driving genetic isolation and creating pockets of diversity to conserve.

Experimental Crossing Confirms Reproductive Isolation between Cryptic Species within Eulimnogammarus verrucosus (Crustacea: Amphipoda) from Lake Baikal

Ancient lakes are known speciation hotspots. One of the most speciose groups in the ancient Lake Baikal are gammaroid amphipods (Crustacea: Amphipoda: Gammaroidea). There are over 350 morphological species and subspecies of amphipods in Baikal, but the extent of cryptic variation is still unclear. One of the most common species in the littoral zone of the lake, Eulimnogammarus verrucosus (Gerstfeldt, 1858), was recently found to comprise at least three (pseudo)cryptic species based on molecular data. Here, we further explored these species by analyzing their mitogenome-based phylogeny, genome sizes with flow cytometry, and their reproductive compatibility. We found divergent times of millions of years and different genome sizes in the three species (6.1, 6.9 and 8 pg), further confirming their genetic separation. Experimental crossing of the western and southern species, which are morphologically indistinguishable and have adjacent ranges, showed their separation with a post-zygotic reproductive barrier, as hybrid embryos stopped developing roughly at the onset of gastrulation. Thus, the previously applied barcoding approach effectively indicated the separate biological species within E. verrucosus. These results provide new data for investigating genome evolution and highlight the need for precise tracking of the sample origin in any studies in this morphospecies.

Ecological, physiological and life-history traits correlate with genome sizes in decapod crustaceans

Crustaceans are characterized by some of the most variable genome sizes among animals. Significant relationships between genome size and specific eco-physiological and morphological features have been described in many crustacean taxa, such as Amphipoda, Ostracoda, Cladocera, and Copepoda. A consistent pattern of genome size variation is yet to be found, however, in one of the most specious order of crustaceans, the Decapoda. We investigated how genome size is related to selected eco-physiological and life-history traits in species belonging to this order. We analyzed the respiratory, excretory, and developmental adaptations of 179 decapod species inhabiting shallow marine, deep marine, freshwater, intertidal, and terrestrial habitats. Our results show that the genome size is significantly larger in decapods that have a direct development. Moreover, in Anomura, Astacidea, and Brachyura we found larger genome sizes in species (i) living in freshwater habitats, (ii) using gills as a strictly water-breathing organs, and (iii) presenting a direct development. While species that (i) live in shallow waters, and intertidal and terrestrial habitats, (ii) have some degrees of air-breathing adaptations, namely gills that can exchange oxygen with air if wet or branchiostegal lungs, and (iii) possess an indirect or extended development are characterized by significantly smaller genomes. Our analyses show that developmental complexity affects genome size in decapods, and that multiple eco-physiological and life-history traits correlate with genome size in Anomura, Astacidea, and Brachyura.

Adaptive Bird-like Genome Miniaturization During the Evolution of Scallop Swimming Lifestyle

Genome miniaturization drives key evolutionary innovations of adaptive traits in vertebrates, such as the flight evolution of birds. However, whether similar evolutionary processes exist in invertebrates remains poorly understood. Derived from the second-largest animal phylum, scallops are a special group of bivalve molluscs and acquire the evolutionary novelty of the swimming lifestyle, providing excellent models for investigating the coordinated genome and lifestyle evolution. Here, we show for the first time that genome sizes of scallops exhibit a generally negative correlation with locomotion activity. To elucidate the co-evolution of genome size and swimming lifestyle, we focus on the Asian moon scallop (Amusium pleuronectes) that possesses the smallest known scallop genome while being among scallops with the highest swimming activity. Whole-genome sequencing of A. pleuronectes reveals highly conserved chromosomal macrosynteny and microsynteny, suggestive of a highly contracted but not degenerated genome. Genome reduction of A. pleuronectes is facilitated by significant inactivation of transposable elements, leading to reduced gene length, elevated expression of genes involved in energy-producing pathways, and decreased copy numbers and expression levels of biomineralization-related genes. Similar evolutionary changes of relevant pathways are also observed for bird genome reduction with flight evolution. The striking mimicry of genome miniaturization underlying the evolution of bird flight and scallop swimming unveils the potentially common, pivotal role of genome size fluctuation in the evolution of novel lifestyles in the animal kingdom.

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.

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.

A high-quality carabid genome assembly provides insights into beetle genome evolution and cold adaptation

The hyperdiverse order Coleoptera comprises a staggering ~25% of known species on Earth. Despite recent breakthroughs in next generation sequencing, there remains a limited representation of beetle diversity in assembled genomes. Most notably, the ground beetle family Carabidae, comprising more than 40,000 described species, has not been studied in a comparative genomics framework using whole genome data. Here we generate a high-quality genome assembly for Nebria riversi, to examine sources of novelty in the genome evolution of beetles, as well as genetic changes associated with specialization to high-elevation alpine habitats. In particular, this genome resource provides a foundation for expanding comparative molecular research into mechanisms of insect cold adaptation. Comparison to other beetles shows a strong signature of genome compaction, with N. riversi possessing a relatively small genome (~147 Mb) compared to other beetles, with associated reductions in repeat element content and intron length. Small genome size is not, however, associated with fewer protein-coding genes, and an analysis of gene family diversity shows significant expansions of genes associated with cellular membranes and membrane transport, as well as protein phosphorylation and muscle filament structure. Finally, our genomic analyses show that these high-elevation beetles have endosymbiotic Spiroplasma, with several metabolic pathways (e.g., propanoate biosynthesis) that might complement N. riversi, although its role as a beneficial symbiont or as a reproductive parasite remains equivocal.

Eurythenes atacamensis sp. nov. (Crustacea: Amphipoda) exhibits ontogenetic vertical stratification across abyssal and hadal depths in the Atacama Trench, eastern South Pacific Ocean

Eurythenes S.I. Smith in Scudder, 1882 (Crustacea: Amphipoda) are prevalent scavengers of the benthopelagic community from bathyal to hadal depths. While a well-studied genus, molecular systematic studies have uncovered cryptic speciation and multiple undescribed lineages. Here, we apply an integrative taxonomic approach and describe the tenth species, Eurythenes atacamensis sp. nov., based on specimens from the 2018 Atacamex and RV Sonne SO261 Expeditions to the southern sector of the Peru-Chile Trench, the Atacama Trench (24-⁠21°S). Eurythenes atacamensis sp. nov. is a large species, max. observed length 83.2 mm, possesses diagnostic features, including a short gnathopod 1 palm and a chelate gnathopod 2 palm, and a distinct genetic lineage based on a 16S rRNA and COI phylogeny. This species is a dominant bait-attending fauna with an extensive bathymetric range, spanning from 4974 to 8081 m. The RV Sonne SO261 specimens were recovered along a 10-station transect from abyssal to hadal depths and further examined for demographic and bathymetric-related patterns. Ontogenetic vertical stratification was evident across the trench axis, with only juveniles present at abyssal depths (4974-6025 m). Total length-depth analysis revealed that the size of females was unrelated to depth, whereas juveniles followed a sigmoidal relationship with a step-up in size at depths >7200 m. Thus, these bathymetric trends suggest that juveniles and females employ differing ecological strategies in subduction trench environments. This study highlights that even dominant and ecologically important species are still being discovered within the abyssal and hadal environments. Continued systematic expeditions will lead to an improved understanding of the eco-evolutionary drivers of speciation in the world's largest ecosystem.

Genome Size Covaries More Positively with Propagule Size than Adult Size: New Insights into an Old Problem

The body size and (or) complexity of organisms is not uniformly related to the amount of genetic material (DNA) contained in each of their cell nuclei ('genome size'). This surprising mismatch between the physical structure of organisms and their underlying genetic information appears to relate to variable accumulation of repetitive DNA sequences, but why this variation has evolved is little understood. Here, I show that genome size correlates more positively with egg size than adult size in crustaceans. I explain this and comparable patterns observed in other kinds of animals and plants as resulting from genome size relating strongly to cell size in most organisms, which should also apply to single-celled eggs and other reproductive propagules with relatively few cells that are pivotal first steps in their lives. However, since body size results from growth in cell size or number or both, it relates to genome size in diverse ways. Relationships between genome size and body size should be especially weak in large organisms whose size relates more to cell multiplication than to cell enlargement, as is generally observed. The ubiquitous single-cell 'bottleneck' of life cycles may affect both genome size and composition, and via both informational (genotypic) and non-informational (nucleotypic) effects, many other properties of multicellular organisms (e.g., rates of growth and metabolism) that have both theoretical and practical significance.

Molecular Cytogenetic Analysis in Freshwater Prawns of the Genus Macrobrachium (Crustacea: Decapoda: Palaemonidae)

Freshwater prawns of the genus Macrobrachium are one of the important components of circumtropical marine, estuarine, and freshwater environments. They have been extensively exploited for human consumption for many years. More than 250 species reflect the evolutionary success of this highly diversified group, with a complex and challenging taxonomy due to morphological variations and vast geographical distribution. Although genetic approaches have been used to clarify phylogenetic and taxonomic aspects of Macrobrachium species, cytogenetic information is still very scarce and mostly focused on chromosome number and morphology. Here, we present chromosome data for three species from the Neotropical region, M. carcinus,M. acanthurus, and M. amazonicum, and one species from the Oriental region, M. rosenbergii. Using conventional cytogenetic approaches and chromosome mapping of repetitive DNAs by fluorescence in situ hybridization (FISH), we identified numerical diversification of the diploid set, within and between both zoogeographic regions. These included M. acanthurus and M. amazonicum sharing diploid chromosomes of 98, while M. carcinus has 94, and M. rosenbergii has 118 chromosomes. Argentophilic sites are also variable in number, but they occur in a much higher number than 18S rDNA, representing two to 10 sites within the study species. Microsatellites repeat motifs are also abundant in the chromosomes, with a co-localization and uniform distribution along the chromosome arms, but completely absent in the AT-rich centromeric regions. As a whole, our study suggests that the 2n divergence was followed by a considerable rDNA diversification. The abundance of the exceptional amount of microsatellite sequences in the chromosomes also suggests that they are essential components of the Macrobrachium genome and, therefore, maintained as a shared feature by the species, the reason for which is yet unknown.

Comparative Proteomics on Deep-Sea Amphipods after in Situ Copper Exposure

The interest in deep-sea mining increased along with the environmental concerns of these activities to the deep-sea fauna. The discovery of optimal biomarkers of deep-sea mining activities in deep-sea species is a crucial step toward the supply of important ecological information for environmental impact assessment. In this study, an in situ copper exposure experiment was performed on deep-sea scavenging amphipods. Abyssorchomene distinctus individuals were selected among all the exposed amphipods for molecular characterization. Copper concentration within the gut was assessed, followed by a tandem mass tag-based coupled with two-dimensional liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) applied to identify and quantify the protein expression changes after 48 h of exposure. 2937 proteins were identified and annotated, and 1918 proteins among all identified proteins were assigned by at least two nonambiguous peptides. The screening process was performed based on the differences in protein abundance and the specific correlation between the proteins and copper in previous studies. These differentially produced proteins include Na⁺/K⁺ ATPase, cuticle, chitinase, and proteins with unknown function. Their abundances showed correlation with copper and had high sensitivity to indicate the copper level, being here proposed as biomarker candidates for deep-sea mining activities in the future. This is a key step in the development of environmental impact assessment of deep-sea mining activities integrating ecotoxicological data.