Survey of genome size in 28 hydrothermal vent species covering 10 families

Genetic sex determination in three closely related hydrothermal vent gastropods, including one species with intersex individuals

Molluscs have undergone many transitions between separate sexes and hermaphroditism, which is of interest in studying the evolution of sex determination and differentiation. Here, we combined multi-locus genotypes obtained from restriction site-associated DNA (RAD) sequencing with anatomical observations of the gonads of three deep-sea hydrothermal vent gastropods of the genus Alviniconcha living in the southwest Pacific. We found that all three species (Alviniconcha boucheti, Alviniconcha strummeri, and Alviniconcha kojimai) share the same male-heterogametic XY sex-determination system but that the gonads of XX A. kojimai individuals are invaded by a variable proportion of male reproductive tissue. The identification of Y-specific RAD loci (found only in A. boucheti) and the phylogenetic analysis of three sex-linked loci shared by all species suggested that X-Y recombination has evolved differently within each species. This situation of three species showing variation in gonadal development around a common sex-determination system provides new insights into the reproductive mode of poorly known deep-sea species and opens up an opportunity to study the evolution of recombination suppression on sex chromosomes and its association with mixed or transitory sexual systems.

The First Genome Survey of the Snail Provanna glabra Inhabiting Deep-Sea Hydrothermal Vents

The snail P. glabra is an endemic species in deep-sea chemosynthetic ecosystems of the Northwest Pacific Ocean. To obtain more genetic information on this species and provide the basis for subsequent whole-genome map construction, a genome survey was performed on this snail from the hydrothermal vent of Okinawa Trough. The genomic size of P. glabra was estimated to be 1.44 Gb, with a heterozygosity of 1.91% and a repeated sequence content of 69.80%. Based on the sequencing data, a draft genome of 1.32 Gb was assembled. Transposal elements (TEs) accounted for 40.17% of the entire genome, with DNA transposons taking the highest proportion. It was found that most TEs were inserted in the genome recently. In the simple sequence repeats, the dinucleotide motif was the most enriched microsatellite type, accounting for 53% of microsatellites. A complete mitochondrial genome of P. glabra with a total length of 16,268 bp was assembled from the sequencing data. After comparison with the published mitochondrial genome of Provanna sp. from a methane seep, 331 potential single nucleotide polymorphism (SNP) sites were identified in protein-coding genes (PCGs). Except for the cox1 gene, nad2, nad4, nad5, and cob genes are expected to be candidate markers for population genetic and phylogenetic studies of P. glabra and other deep-sea snails. Compared with shallow-water species, three mitochondrial genes of deep-sea gastropods exhibited a higher evolutionary rate, indicating strong selection operating on mitochondria of deep-sea species. This study provides insights into the genome characteristics of P. glabra and supplies genomic resources for further studies on the adaptive evolution of the snail in extreme deep-sea chemosynthetic environments.

Genome sizes and repeatome evolution in zoantharians (Cnidaria: Hexacorallia: Zoantharia)

Across eukaryotes, large variations of genome sizes have been observed even between closely related species. Transposable elements as part of the repeated DNA have been proposed and confirmed as one of the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and transposable element dynamics are not well understood. Together with phenotypic traits, they are commonly referred to as the "C-value enigma". The order Zoantharia are benthic cnidarians found from intertidal zones to the deep sea, and some species are particularly abundant in coral reefs. Despite their high ecological relevance, zoantharians have yet to be largely studied from the genomic point of view. This study aims at investigating the role of the repeatome (total content of repeated elements) in genome size variations across the order Zoantharia. To this end, whole-genomes of 32 zoantharian species representing five families were sequenced. Genome sizes were estimated and the abundances of different repeat classes were assessed. In addition, the repeat overlap between species was assessed by a sequence clustering method. The genome sizes in the dataset varied up to 2.4 fold magnitude. Significant correlations between genome size, repeated DNA content and transposable elements, respectively (Pearson's correlation test R2 = 0.47, p = 0.0016; R2 = 0.22, p = 0.05) were found, suggesting their involvement in the dynamics of genome expansion and reduction. In all species, long interspersed nuclear elements and DNA transposons were the most abundant identified elements. These transposable elements also appeared to have had a recent expansion event. This was in contrast to the comparative clustering analysis which revealed species-specific patterns of satellite elements' amplification. In summary, the genome sizes of zoantharians likely result from the complex dynamics of repeated elements. Finally, the majority of repeated elements (up to 70%) could not be annotated to a known repeat class, highlighting the need to further investigate non-model cnidarian genomes. More research is needed to understand how repeated DNA dynamics relate to zoantharian evolution and their biology.

Distinct genomic routes underlie transitions to specialised symbiotic lifestyles in deep-sea annelid worms

Bacterial symbioses allow annelids to colonise extreme ecological niches, such as hydrothermal vents and whale falls. Yet, the genetic principles sustaining these symbioses remain unclear. Here, we show that different genomic adaptations underpin the symbioses of phylogenetically related annelids with distinct nutritional strategies. Genome compaction and extensive gene losses distinguish the heterotrophic symbiosis of the bone-eating worm Osedax frankpressi from the chemoautotrophic symbiosis of deep-sea Vestimentifera. Osedax's endosymbionts complement many of the host's metabolic deficiencies, including the loss of pathways to recycle nitrogen and synthesise some amino acids. Osedax's endosymbionts possess the glyoxylate cycle, which could allow more efficient catabolism of bone-derived nutrients and the production of carbohydrates from fatty acids. Unlike in most Vestimentifera, innate immunity genes are reduced in O. frankpressi, which, however, has an expansion of matrix metalloproteases to digest collagen. Our study supports that distinct nutritional interactions influence host genome evolution differently in highly specialised symbioses.

Patterns of genome size variation in caridean shrimps: new estimates for non-gambarelloides Synalpheus species

Genome size (GS) or DNA nuclear content is considered a useful index for making inferences about evolutionary models and life history in animals, including taxonomic, biogeographical, and ecological scenarios. However, patterns of GS variation and their causes in crustaceans are still poorly understood. This study aimed to describe the GS of five Neotropical Synalpheus nongambarelloides shrimps (S. apioceros, S. minus, S. brevicarpus, S. fritzmueller, and S. scaphoceris) and compare the C-values of all Caridea Infraorder in terms of geography and phylogenetics. All animals were sampled in the coast of São Paulo State, Brazil and GS was assessed by flow cytometry analysis (FCA). The C-values ranged from 7.89 pg in S. apioceros to 12.24 pg in S. scaphoceris. Caridean shrimps had higher GS than other Decapoda crustaceans. The results reveal a tendency of obtaining larger genomes in species with direct development in Synalpheus shrimps. In addition, a tendency of positive biogeographical (latitudinal) correlation with Caridea Infraorder was also observed. This study provides remarkable and new protocol for FCA (using gating strategy for the analysis), which led to the discovery of new information regarding GS of caridean shrimps, especially for Neotropical Synalpheus, which represents the second-largest group in the Caridea Infraorder.

Novel insights on obligate symbiont lifestyle and adaptation to chemosynthetic environment as revealed by the giant tubeworm genome

The mutualism between the giant tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has been extensively researched over the past 40 years. However, the lack of the host whole genome information has impeded the full comprehension of the genotype/phenotype interface in Riftia. Here we described the high-quality draft genome of Riftia, its complete mitogenome, and tissue-specific transcriptomic data. The Riftia genome presents signs of reductive evolution, with gene family contractions exceeding expansions. Expanded gene families are related to sulphur metabolism, detoxification, anti-oxidative stress, oxygen transport, immune system, and lysosomal digestion, reflecting evolutionary adaptations to the vent environment and endosymbiosis. Despite the derived body plan, the developmental gene repertoire in the gutless tubeworm is extremely conserved with the presence of a near intact and complete Hox cluster. Gene expression analyses establishes that the trophosome is a multi-functional organ marked by intracellular digestion of endosymbionts, storage of excretory products and haematopoietic functions. Overall, the plume and gonad tissues both in contact to the environment harbour highly expressed genes involved with cell cycle, programmed cell death, and immunity indicating a high cell turnover and defence mechanisms against pathogens. We posit that the innate immune system plays a more prominent role into the establishment of the symbiosis during the infection in the larval stage, rather than maintaining the symbiostasis in the trophosome. This genome bridges four decades of physiological research in Riftia, whilst simultaneously provides new insights into the development, whole organism functions and evolution in the giant tubeworm.

Chromosome-level genome assembly of Paralithodes platypus provides insights into evolution and adaptation of king crabs

The blue king crab, Paralithodes platypus, which belongs to the family Lithodidae, is a commercially and ecologically important species. However, a high-quality reference genome for the king crab has not yet been reported. Here, we assembled the first chromosome-level blue king crab genome, which contains 104 chromosomes and an N50 length of 51.15 Mb. Furthermore, we determined that the large genome size can be attributed to the insertion of long interspersed nuclear elements and long tandem repeats. Genome assembly assessment showed that 96.54% of the assembled transcripts could be aligned to the assembled genome. Phylogenetic analysis showed the blue king crab to have a close relationship with the Eubrachyura crabs, from which it diverged 272.5 million years ago. Population history analyses indicated that the effective population of the blue king crab declined sharply and then gradually increased from the Cretaceous and Neogene periods, respectively. Furthermore, gene families related to developmental pathways, steroid and thyroid hormone synthesis, and inflammatory regulation were expanded in the genome, suggesting that these genes contributed substantially to the environmental adaptation and unique body plan evolution of the blue king crab. The high-quality reference genome reported here provides a solid molecular basis for further study of the blue king crab's development and environmental adaptation.

The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest

BACKGROUND

Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture.

RESULTS

Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world.

CONCLUSIONS

The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

Production of WW males lacking the masculine Z chromosome and mining the Macrobrachium rosenbergii genome for sex-chromosomes

The cultivation of monosex populations is common in animal husbandry. However, preselecting the desired gender remains a major biotechnological and ethical challenge. To achieve an efficient biotechnology for all-female aquaculture in the economically important prawn (Macrobrachium rosenbergii), we achieved – for the first time – WW males using androgenic gland cells transplantation which caused full sex-reversal of WW females to functional males. Crossing the WW males with WW females yielded all-female progeny lacking the Z chromosome. We now have the ability to manipulate – by non-genomic means – all possible genotype combinations (ZZ, WZ and WW) to retain either male or female phenotypes and hence to produce monosex populations of either gender. This calls for a study of the genomic basis underlying this striking sexual plasticity, questioning the content of the W and Z chromosomes. Here, we report on the sequencing of a high-quality genome exhibiting distinguishable paternal and maternal sequences. This assembly covers ~ 87.5% of the genome and yielded a remarkable N50 value of ~ 20 × 10⁶ bp. Genomic sex markers were used to initiate the identification and validation of parts of the W and Z chromosomes for the first time in arthropods.

The Genome Sizes of Ostracod Crustaceans Correlate with Body Size and Evolutionary History, but not Environment

Within animals, a positive correlation between genome size and body size has been detected in several taxa but not in others, such that it remains unknown how pervasive this pattern may be. Here, we provide another example of a positive relationship in a group of crustaceans whose genome sizes have not previously been investigated. We analyze genome size estimates for 46 species across the 2 most diverse orders of Class Ostracoda, commonly known as seed shrimps, including 29 new estimates made using Feulgen image analysis densitometry and flow cytometry. Genome sizes in this group range ~80-fold, a level of variability that is otherwise not seen in crustaceans with the exception of some malacostracan orders. We find a strong positive correlation between genome size and body size across all species, including after phylogenetic correction. We additionally detect evidence of XX/XO sex determination in 3 species of marine ostracods where male and female genome sizes were estimated. On average, genome sizes are larger but less variable in Order Myodocopida than in Order Podocopida, and marine ostracods have larger genomes than freshwater species, but this appears to be explained by phylogenetic inertia. The relationship between phylogeny, genome size, body size, and habitat is complex in this system and provides a baseline for future studies examining the interactions of these biological traits.

Genome Sequences of Marine Shrimp Exopalaemon carinicauda Holthuis Provide Insights into Genome Size Evolution of Caridea

Crustacea, particularly Decapoda, contains many economically important species, such as shrimps and crabs. Crustaceans exhibit enormous (nearly 500-fold) variability in genome size. However, limited genome resources are available for investigating these species. Exopalaemoncarinicauda Holthuis, an economical caridean shrimp, is a potential ideal experimental animal for research on crustaceans. In this study, we performed low-coverage sequencing and de novo assembly of the E. carinicauda genome. The assembly covers more than 95% of coding regions. E. carinicauda possesses a large complex genome (5.73 Gb), with size twice higher than those of many decapod shrimps. As such, comparative genomic analyses were implied to investigate factors affecting genome size evolution of decapods. However, clues associated with genome duplication were not identified, and few horizontally transferred sequences were detected. Ultimately, the burst of transposable elements, especially retrotransposons, was determined as the major factor influencing genome expansion. A total of 2 Gb repeats were identified, and RTE-BovB, Jockey, Gypsy, and DIRS were the four major retrotransposons that significantly expanded. Both recent (Jockey and Gypsy) and ancestral (DIRS) originated retrotransposons responsible for the genome evolution. The E. carinicauda genome also exhibited potential for the genomic and experimental research of shrimps.

Comparative transcriptome analysis of Rimicaris sp. reveals novel molecular features associated with survival in deep-sea hydrothermal vent

Shrimp of the family Alvinocarididae are the predominant megafauna of deep-sea hydrothermal vents. However, genome information on this family is currently unavailable. In the present study, by employing Illumina sequencing, we performed the first de novo transcriptome analysis of the gills of the shrimp Rimicaris sp. from the hydrothermal vent in Desmos, Manus Basin. The analysis was conducted in a comparative manner with the shrimp taken directly from the vent (GR samples) and the shrimp that had been maintained for ten days under normal laboratory condition (mGR samples). Among the 128,938 unigenes identified, a large number of differentially expressed genes (DEGs) between the GR and mGR samples were detected, including 2365 and 1607 genes significantly upregulated and downregulated, respectively, in GR. The DEGs covered diverse functional categories. Most of the DEGs associated with immunity were downregulated in GR, while most of the DEGs associated with sulfur metabolism and detoxification were upregulated in GR. These results provide the first comprehensive transcriptomic resource for hydrothermal vent Rimicaris and revealed varied categories of genes likely involved in deep-sea survival.

Patterns of genome size variation in snapping shrimp

Although crustaceans vary extensively in genome size, little is known about how genome size may affect the ecology and evolution of species in this diverse group, in part due to the lack of large genome size datasets. Here we investigate interspecific, intraspecific, and intracolony variation in genome size in 39 species of Synalpheus shrimps, representing one of the largest genome size datasets for a single genus within crustaceans. We find that genome size ranges approximately 4-fold across Synalpheus with little phylogenetic signal, and is not related to body size. In a subset of these species, genome size is related to chromosome size, but not to chromosome number, suggesting that despite large genomes, these species are not polyploid. Interestingly, there appears to be 35% intraspecific genome size variation in Synalpheus idios among geographic regions, and up to 30% variation in Synalpheus duffyi genome size within the same colony.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Phylogeny and New Classification of Hydrothermal Vent and Seep Shrimps of the Family Alvinocarididae (Decapoda)

The paper addresses the phylogeny and classification of the hydrothermal vent shrimp family Alvinocarididae. Two morphological cladistic analyses were carried out, which use all 31 recognized species of Alvinocarididae as terminal taxa. As outgroups, two species were included, both representing major caridean clades: Acanthephyra purpurea (Acanthephyridae) and Alpheus echiurophilus (Alpheidae). For additional support of the clades we utilised available data on mitochondrial Cytochrome c Oxidase I gene (CO1) and 16S ribosomal markers. Both morphological and molecular methods resulted in similar tree topologies and nearly identical clades. We consider these clades as evolutionary units and thus erect two new subfamilies: Rimicaridinae (Alvinocaridinides, Manuscaris, Opaepele, Shinkaicaris, Rimicaris), Alvinocaridinae (Alvinocaris), whilst recognising Mirocaridinae (with genera Mirocaris and Nautilocaris) at subfamily level. One genus, Keldyshicaris could not be assigned to any subfamily and is thus left as incertae sedis. The monophyly of Alvinocardinae was supported by morphological data, but not supported by molecular data (two analyses); the monophyly of all subfamilies was supported both by morphological and molecular data. Chorocaris is herein synonymized with Rimicaris, whilst Opaepele vavilovi is herein transferred to a new genus Keldyshicaris. Morphological trends within Alvinocarididae are discussed and short biogeographical remarks are given. We provide emended diagnoses for all subfamilies and genera along with keys to all recognized species.

LTR-retrotransposons in R. exoculata and other crustaceans: the outstanding success of GalEa-like copia elements

Transposable elements are major constituents of eukaryote genomes and have a great impact on genome structure and stability. They can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution among several genomes is an essential condition to study their dynamics and to better understand their role in species evolution. LTR-retrotransposons have been reported in many diverse eukaryote species, describing a ubiquitous distribution. Given their abundance, diversity and their extended ranges in C-values, environment and life styles, crustaceans are a great taxon to investigate the genomic component of adaptation and its possible relationships with TEs. However, crustaceans have been greatly underrepresented in transposable element studies. Using both degenerate PCR and in silico approaches, we have identified 35 Copia and 46 Gypsy families in 15 and 18 crustacean species, respectively. In particular, we characterized several full-length elements from the shrimp Rimicaris exoculata that is listed as a model organism from hydrothermal vents. Phylogenic analyses show that Copia and Gypsy retrotransposons likely present two opposite dynamics within crustaceans. The Gypsy elements appear relatively frequent and diverse whereas Copia are much more homogeneous, as 29 of them belong to the single GalEa clade, and species- or lineage-dependent. Our results also support the hypothesis of the Copia retrotransposon scarcity in metazoans compared to Gypsy elements. In such a context, the GalEa-like elements present an outstanding wide distribution among eukaryotes, from fishes to red algae, and can be even highly predominant within a large taxon, such as Malacostraca. Their distribution among crustaceans suggests a dynamics that follows a "domino days spreading" branching process in which successive amplifications may interact positively.

Insights into metazoan evolution from Alvinella pompejana cDNAs

Background

Alvinella pompejana is a representative of Annelids, a key phylum for evo-devo studies that is still poorly studied at the sequence level. A. pompejana inhabits deep-sea hydrothermal vents and is currently known as one of the most thermotolerant Eukaryotes in marine environments, withstanding the largest known chemical and thermal ranges (from 5 to 105°C). This tube-dwelling worm forms dense colonies on the surface of hydrothermal chimneys and can withstand long periods of hypo/anoxia and long phases of exposure to hydrogen sulphides. A. pompejana specifically inhabits chimney walls of hydrothermal vents on the East Pacific Rise. To survive, Alvinella has developed numerous adaptations at the physiological and molecular levels, such as an increase in the thermostability of proteins and protein complexes. It represents an outstanding model organism for studying adaptation to harsh physicochemical conditions and for isolating stable macromolecules resistant to high temperatures.

Results

We have constructed four full length enriched cDNA libraries to investigate the biology and evolution of this intriguing animal. Analysis of more than 75,000 high quality reads led to the identification of 15,858 transcripts and 9,221 putative protein sequences. Our annotation reveals a good coverage of most animal pathways and networks with a prevalence of transcripts involved in oxidative stress resistance, detoxification, anti-bacterial defence, and heat shock protection. Alvinella proteins seem to show a slow evolutionary rate and a higher similarity with proteins from Vertebrates compared to proteins from Arthropods or Nematodes. Their composition shows enrichment in positively charged amino acids that might contribute to their thermostability. The gene content of Alvinella reveals that an important pool of genes previously considered to be specific to Deuterostomes were in fact already present in the last common ancestor of the Bilaterian animals, but have been secondarily lost in model invertebrates. This pool is enriched in glycoproteins that play a key role in intercellular communication, hormonal regulation and immunity.

Conclusions

Our study starts to unravel the gene content and sequence evolution of a deep-sea annelid, revealing key features in eukaryote adaptation to extreme environmental conditions and highlighting the proximity of Annelids and Vertebrates.