Evolution, Expression Patterns, and Distribution of Novel Ribbon Worm Predatory and Defensive Toxins.

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.

Loss of S-nitrosoglutathione reductase disturbs phytohormone homeostasis and regulates shoot side branching and fruit growth in tomato.

S-Nitrosoglutathione plays a central role in nitric oxide (NO) homeostasis, and S-nitrosoglutathione reductase (GSNOR) regulates the cellular levels of S-nitrosoglutathione across kingdoms. Here, we investigated the role of endogenous NO in shaping shoot architecture and controlling fruit set and growth in tomato (Solanum lycopersicum). SlGSNOR silencing promoted shoot side branching and led to reduced fruit size, negatively impacting fruit yield. Greatly intensified in slgsnor knockout plants, these phenotypical changes were virtually unaffected by SlGSNOR overexpression. Silencing or knocking out of SlGSNOR intensified protein tyrosine nitration and S-nitrosation and led to aberrant auxin production and signaling in leaf primordia and fruit-setting ovaries, besides restricting the shoot basipetal polar auxin transport stream. SlGSNOR deficiency triggered extensive transcriptional reprogramming at early fruit development, reducing pericarp cell proliferation due to restrictions on auxin, gibberellin, and cytokinin production and signaling. Abnormal chloroplast development and carbon metabolism were also detected in early-developing NO-overaccumulating fruits, possibly limiting energy supply and building blocks for fruit growth. These findings provide new insights into the mechanisms by which endogenous NO fine-tunes the delicate hormonal network controlling shoot architecture, fruit set, and post-anthesis fruit development, emphasizing the relevance of NO-auxin interaction for plant development and productivity.

Multifaceted roles of nitric oxide in tomato fruit ripening: NO-induced metabolic rewiring and consequences for fruit quality traits.

Nitric oxide (NO) has been implicated as part of the ripening regulatory network in fleshy fruits. However, very little is known about the simultaneous action of NO on the network of regulatory events and metabolic reactions behind ripening-related changes in fruit color, taste, aroma and nutritional value. Here, we performed an in-depth characterization of the concomitant changes in tomato (Solanum lycopersicum) fruit transcriptome and metabolome associated with the delayed-ripening phenotype caused by NO supplementation at the pre-climacteric stage. Approximately one-third of the fruit transcriptome was altered in response to NO, including a multilevel down-regulation of ripening regulatory genes, which in turn restricted the production and tissue sensitivity to ethylene. NO also repressed hydrogen peroxide-scavenging enzymes, intensifying nitro-oxidative stress and S-nitrosation and nitration events throughout ripening. Carotenoid, tocopherol, flavonoid and ascorbate biosynthesis were differentially affected by NO, resulting in overaccumulation of ascorbate (25%) and flavonoids (60%), and impaired lycopene production. In contrast, the biosynthesis of compounds related to tomato taste (sugars, organic acids, amino acids) and aroma (volatiles) was slightly affected by NO. Our findings indicate that NO triggers extensive transcriptional and metabolic rewiring at the early ripening stage, modifying tomato antioxidant composition with minimal impact on fruit taste and aroma.

C4 and crassulacean acid metabolism within a single leaf: deciphering key components behind a rare photosynthetic adaptation.

Although biochemically related, C4 and crassulacean acid metabolism (CAM) systems are expected to be incompatible. However, Portulaca species, including P. oleracea, operate C4 and CAM within a single leaf, and the mechanisms behind this unique photosynthetic arrangement remain largely unknown. Here, we employed RNA-seq to identify candidate genes involved exclusively or shared by C4 or CAM, and provided an in-depth characterization of their transcript abundance patterns during the drought-induced photosynthetic transitions in P. oleracea. Data revealed fewer candidate CAM-specific genes than those recruited to function in C4 . The putative CAM-specific genes were predominantly involved in night-time primary carboxylation reactions and malate movement across the tonoplast. Analysis of gene transcript-abundance regulation and photosynthetic physiology indicated that C4 and CAM coexist within a single P. oleracea leaf under mild drought conditions. Developmental and environmental cues were shown to regulate CAM expression in stems, whereas the shift from C4 to C4 -CAM hybrid photosynthesis in leaves was strictly under environmental control. Moreover, efficient starch turnover was identified as part of the metabolic adjustments required for CAM operation in both organs. These findings provide insights into C4 /CAM connectivity and compatibility, contributing to a deeper understanding of alternative ways to engineer CAM into C4 crop species.

Role of neuroendocrine modulation and biochemistry in the sepsis in Piaractus mesopotamicus.

Sepsis is a systemic process with multifactorial pathophysiology that affects most animal species. It is responsible for high rates of morbidity and mortality. This work aimed to study the biochemical and neuroendocrine changes of the sepsis process in Piaractus mesopotamicus after Aeromonas hydrophila inoculation analyzing changes in blood leukocyte and differences in neuroendocrine-biochemical modulation using RNA-seq. Fish showed hypercortisolemia, inhibition of glucose absorption, followed by hypocortisolemia and then hyperglycemia. Thyroid hormones (T3 and T4) showed immediate decrease in serum and T4 increased 6 h post-inoculation (HPI). Sepsis-induced hormonal alterations triggered changes in the metabolic pathways increasing protein and lipid catabolism, use of transient anaerobic glycolysis and liver injury. A reference transcriptome was constructed based on blood leukocytes from P. mesopotamicus. The assembly resulted in total 266,272 contigs with a N50 of 2786 bp. There was a reorganization of plasma membrane of leukocytes at the beginning of the septic process with increased expression of neuroendocrine receptors and with continuous flow of neurotransmitters, hormones and solutes with compensatory regulation at 6 HPI. Three and nine HPI seemed to be critical, the expression of a number of transcription factors was increased, including the modulatory DEGs related to glucocorticoid and thyroid hormones induced and suppressed (FDR < 0.05). Neuroendocrine modulation can regulate leukocytes and biochemical parameters of peripheral blood, being important sources for the study of the pathophysiology of sepsis. These finding highlights the importance of further studies focusing on biochemical-neuroendocrine changes in blood leukocytes and systemic sepsis.

Revisiting metazoan phylogeny with genomic sampling of all phyla.

Proper biological interpretation of a phylogeny can sometimes hinge on the placement of key taxa-or fail when such key taxa are not sampled. In this light, we here present the first attempt to investigate (though not conclusively resolve) animal relationships using genome-scale data from all phyla. Results from the site-heterogeneous CAT + GTR model recapitulate many established major clades, and strongly confirm some recent discoveries, such as a monophyletic Lophophorata, and a sister group relationship between Gnathifera and Chaetognatha, raising continued questions on the nature of the spiralian ancestor. We also explore matrix construction with an eye towards testing specific relationships; this approach uniquely recovers support for Panarthropoda, and shows that Lophotrochozoa (a subclade of Spiralia) can be constructed in strongly conflicting ways using different taxon- and/or orthologue sets. Dayhoff-6 recoding sacrifices information, but can also reveal surprising outcomes, e.g. full support for a clade of Lophophorata and Entoprocta + Cycliophora, a clade of Placozoa + Cnidaria, and raising support for Ctenophora as sister group to the remaining Metazoa, in a manner dependent on the gene and/or taxon sampling of the matrix in question. Future work should test the hypothesis that the few remaining uncertainties in animal phylogeny might reflect violations of the various stationarity assumptions used in contemporary inference methods.

Biological and molecular characterization of a putative new potexvirus infecting Senna occidentalis.

In this work, we report the complete genome sequence of, production of polyclonal antibodies against, and development of biological assays for a putative new potexvirus, named senna mosaic virus (SenMV), found infecting Senna occidentalis in the state of São Paulo, Brazil. The complete genome sequence of SenMV comprises 6775 nucleotides excluding the poly(A) tail. The genome organization is similar to those of other potexviruses, with five open reading frames coding for RNA-dependent RNA polymerase (RdRp), the triple gene block (TGB 1, 2, and 3) proteins, and coat protein (CP). The virus was transmitted to S. occidentalis by mechanical inoculation and trimming scissors, but not by seeds.

Resolving the evolutionary relationships of molluscs with phylogenomic tools.

Molluscs (snails, octopuses, clams and their relatives) have a great disparity of body plans and, among the animals, only arthropods surpass them in species number. This diversity has made Mollusca one of the best-studied groups of animals, yet their evolutionary relationships remain poorly resolved. Open questions have important implications for the origin of Mollusca and for morphological evolution within the group. These questions include whether the shell-less, vermiform aplacophoran molluscs diverged before the origin of the shelled molluscs (Conchifera) or lost their shells secondarily. Monoplacophorans were not included in molecular studies until recently, when it was proposed that they constitute a clade named Serialia together with Polyplacophora (chitons), reflecting the serial repetition of body organs in both groups. Attempts to understand the early evolution of molluscs become even more complex when considering the large diversity of Cambrian fossils. These can have multiple dorsal shell plates and sclerites or can be shell-less but with a typical molluscan radula and serially repeated gills. To better resolve the relationships among molluscs, we generated transcriptome data for 15 species that, in combination with existing data, represent for the first time all major molluscan groups. We analysed multiple data sets containing up to 216,402 sites and 1,185 gene regions using multiple models and methods. Our results support the clade Aculifera, containing the three molluscan groups with spicules but without true shells, and they support the monophyly of Conchifera. Monoplacophora is not the sister group to other Conchifera but to Cephalopoda. Strong support is found for a clade that comprises Scaphopoda (tusk shells), Gastropoda and Bivalvia, with most analyses placing Scaphopoda and Gastropoda as sister groups. This well-resolved tree will constitute a framework for further studies of mollusc evolution, development and anatomy.

Differences in the skeletal muscle transcriptome profile associated with extreme values of fatty acids content.

BACKGROUND: Lipids are a class of molecules that play an important role in cellular structure and metabolism in all cell types. In the last few decades, it has been reported that long-chain fatty acids (FAs) are involved in several biological functions from transcriptional regulation to physiological processes. Several fatty acids have been both positively and negatively implicated in different biological processes in skeletal muscle and other tissues. To gain insight into biological processes associated with fatty acid content in skeletal muscle, the aim of the present study was to identify differentially expressed genes (DEGs) and functional pathways related to gene expression regulation associated with FA content in cattle. RESULTS: Skeletal muscle transcriptome analysis of 164 Nellore steers revealed no differentially expressed genes (DEGs, FDR 10%) for samples with extreme values for linoleic acid (LA) or stearic acid (SA), and only a few DEGs for eicosapentaenoic acid (EPA, 5 DEGs), docosahexaenoic acid (DHA, 4 DEGs) and palmitic acid (PA, 123 DEGs), while large numbers of DEGs were associated with oleic acid (OA, 1134 DEGs) and conjugated linoleic acid cis9 trans11 (CLA-c9t11, 872 DEGs). Functional annotation and functional enrichment from OA DEGs identified important genes, canonical pathways and upstream regulators such as SCD, PLIN5, UCP3, CPT1, CPT1B, oxidative phosphorylation mitochondrial dysfunction, PPARGC1A, and FOXO1. Two important genes associated with lipid metabolism, gene expression and cancer were identified as DEGs between animals with high and low CLA-c9t11, specifically, epidermal growth factor receptor (EGFR) and RNPS. CONCLUSION: Only two out of seven classes of molecules of FA studied were associated with large changes in the expression profile of skeletal muscle. OA and CLA-c9t11 content had significant effects on the expression level of genes related to important biological processes associated with oxidative phosphorylation, and cell growth, survival, and migration. These results contribute to our understanding of how some FAs modulate metabolism and may have protective health function.

Articulating "Archiannelids": Phylogenomics and Annelid Relationships, with Emphasis on Meiofaunal Taxa.

Annelid disparity has resulted in morphological-based classifications that disagree with phylogenies based on Sanger sequencing and phylogenomic analyses. However, the data used for the latter studies came from various sources and technologies, involved poorly occupied matrices and lacked key lineages. Here, we generated a new Illumina-based data set to address annelid relationships from a fresh perspective, independent from previously generated data and with nearly fully occupied matrices. Our sampling reflects the span of annelid diversity, including two symbiotic annelid groups (Myzostomida and Spinther) and five meiofaunal groups once referred to as part of Archiannelida (three from Protodrilida, plus Dinophilus and Polygordius). As well as the placement of these unusual annelids, we sought to address the overall phylogeny of Annelida, and provide a new perspective for naming of major clades. Our results largely corroborate the phylogenomic results of Weigert et al. (2014; Illuminating the base of the annelid tree using transcriptomics. Mol Biol Evol. 31:1391-1401), with "Magelona + Owenia" and Chaetopteridae forming a grade with respect to all other annelids. Echiura and Sipuncula are supported as being annelid groups, with Sipuncula closest to amphinomids as sister group to Sedentaria and Errantia. We recovered the three Protodrilida terminals as sister clade to Phyllodocida and Eunicida (=clade Aciculata). We therefore place Protodrilida as part of Errantia. Polygordius was found to be sister group to the scaleworm terminal and the possibility that it is a simplified scaleworm clade, as has been shown for the former family Pisionidae, is discussed. Our results were equivocal with respect to Dinophilus, Myzostomida, and Spinther possibly owing to confounding long-branch effects.

Phylogenomic analyses of a Mediterranean earthworm family (Annelida: Hormogastridae).

Earthworm taxonomy and evolutionary biology remain a challenge because of their scarce distinct morphological characters of taxonomic value, the morphological convergence by adaptation to the uniformity of the soil where they inhabit, and their high plasticity when challenged with stressful or new environmental conditions. Here we present a phylogenomic study of the family Hormogastridae, representing also the first piece of work of this type within earthworms. We included seven transcriptomes of the group representing the main lineages as previously-described, analysed in a final matrix that includes twelve earthworms and eleven outgroups. While there is a high degree of gene conflict in the generated trees that obscure some of the internal relationships, the origin of the family is well resolved: the hormogastrid Hemigastrodrilus appears as the most ancestral group, followed by the ailoscolecid Ailoscolex, therefore rejecting the validity of the family Ailoscolecidae. Our results place the origin of hormogastrids in Southern France, as previously hypothesised.

A transcriptomic approach to ribbon worm systematics (nemertea): resolving the pilidiophora problem.

Resolving the deep relationships of ancient animal lineages has proven difficult using standard Sanger-sequencing approaches with a handful of markers. We thus reassess the relatively well-studied phylogeny of the phylum Nemertea (ribbon worms)-for which the targeted gene approaches had resolved many clades but had left key phylogenetic gaps-by using a phylogenomic approach using Illumina-based de novo assembled transcriptomes and automatic orthology prediction methods. The analysis of a concatenated data set of 2,779 genes (411,138 amino acids) with about 78% gene occupancy and a reduced version with 95% gene occupancy, under evolutionary models accounting or not for site-specific amino acid replacement patterns results in a well-supported phylogeny that recovers all major accepted nemertean clades with the monophyly of Heteronemertea, Hoplonemertea, Monostilifera, being well supported. Significantly, all the ambiguous patterns inferred from Sanger-based approaches were resolved, namely the monophyly of Palaeonemertea and Pilidiophora. By testing for possible conflict in the analyzed supermatrix, we observed that concatenation was the best solution, and the results of the analyses should settle prior debates on nemertean phylogeny. The study highlights the importance, feasibility, and completeness of Illumina-based phylogenomic data matrices.

A phylogenetic backbone for Bivalvia: an RNA-seq approach.

Bivalves are an ancient and ubiquitous group of aquatic invertebrates with an estimated 10 000-20 000 living species. They are economically significant as a human food source, and ecologically important given their biomass and effects on communities. Their phylogenetic relationships have been studied for decades, and their unparalleled fossil record extends from the Cambrian to the Recent. Nevertheless, a robustly supported phylogeny of the deepest nodes, needed to fully exploit the bivalves as a model for testing macroevolutionary theories, is lacking. Here, we present the first phylogenomic approach for this important group of molluscs, including novel transcriptomic data for 31 bivalves obtained through an RNA-seq approach, and analyse these data with published genomes and transcriptomes of other bivalves plus outgroups. Our results provide a well-resolved, robust phylogenetic backbone for Bivalvia with all major lineages delineated, addressing long-standing questions about the monophyly of Protobranchia and Heterodonta, and resolving the position of particular groups such as Palaeoheterodonta, Archiheterodonta and Anomalodesmata. This now fully resolved backbone demonstrates that genomic approaches using hundreds of genes are feasible for resolving phylogenetic questions in bivalves and other animals.

Re-evaluating the phylogeny of Sipuncula through transcriptomics.

Sipunculans (also known as peanut worms) are an ancient group of exclusively marine worms with a global distribution and a fossil record that dates back to the Early Cambrian. The systematics of sipunculans, now considered a distinct subclade of Annelida, has been studied for decades using morphological and molecular characters, and has reached the limits of Sanger-based approaches. Here, we reevaluate their family-level phylogeny by comparative transcriptomic analysis of eight species representing all known families within Sipuncula. Two data matrices with alternative gene occupancy levels (large matrix with 675 genes and 62% missing data; reduced matrix with 141 genes and 23% missing data) were analysed using concatenation and gene-tree methods, yielding congruent results and resolving each internal node with maximum support. We thus corroborate prior phylogenetic work based on molecular data, resolve outstanding issues with respect to the familial relationships of Aspidosiphonidae, Antillesomatidae and Phascolosomatidae, and highlight the next area of focus for sipunculan systematics.

Phylogenomic analyses of deep gastropod relationships reject Orthogastropoda.

Gastropods are a highly diverse clade of molluscs that includes many familiar animals, such as limpets, snails, slugs and sea slugs. It is one of the most abundant groups of animals in the sea and the only molluscan lineage that has successfully colonized land. Yet the relationships among and within its constituent clades have remained in flux for over a century of morphological, anatomical and molecular study. Here, we re-evaluate gastropod phylogenetic relationships by collecting new transcriptome data for 40 species and analysing them in combination with publicly available genomes and transcriptomes. Our datasets include all five main gastropod clades: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia. We use two different methods to assign orthology, subsample each of these matrices into three increasingly dense subsets, and analyse all six of these supermatrices with two different models of molecular evolution. All 12 analyses yield the same unrooted network connecting the five major gastropod lineages. This reduces deep gastropod phylogeny to three alternative rooting hypotheses. These results reject the prevalent hypothesis of gastropod phylogeny, Orthogastropoda. Our dated tree is congruent with a possible end-Permian recovery of some gastropod clades, namely Caenogastropoda and some Heterobranchia subclades.

Continent-island boundary and environment-shaped evolution in the marine amphipod Ampithoe marcuzzii complex (Crustacea: Eumalacostraca: Ampithoidae).

Marine amphipods are crustaceans that lack a larval phase and consequently have low dispersion rates. Despite that, these crustaceans present a remarkable ability to be transported by rafting on natural floating substrata, especially macroalgae, where they find shelter, food and a mating ground. The species Ampithoe marcuzzii is widely distributed throughout the western Atlantic Ocean. Here, it was used as a model to study seascape genomics and phylogeography in invertebrates with low dispersion capacities. We anticipated that the lineages would present isolation-by-distance patterns. However, surface currents and other abiotic variables could facilitate connectivity among distant sites. Based on mitochondrial and nuclear genes, SNPs, and environmental associations, we observed the presence of a species complex within A. marcuzzii, separating mainland and insular populations. Each species showed an independent evolutionary history, with a strong latitudinal population structure and evidence of isolation-by-distance and isolation-by-environment, characterizing the 'continent' species. Historical expansion and environmental variables were observed associated with the southeastern population, and ecological niche modeling corroborated the region as a paleorefuge. Conversely, populations from 'islands' presented complicated evolutionary histories, with closer localities genetically isolated and distant localities connected. These findings indicate that insular populations with low dispersion capacity might be more susceptible to spatial connectivity by floating substrata and to changes in surface currents. In contrast, mainland populations might be more vulnerable to local climate changes due to lack of gene flow.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication.

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Comparing iDNA from mosquitoes and flies to survey mammals in a semi-controlled Neotropical area.

Ingested-derived DNA (iDNA) from insects represents a powerful tool for assessing vertebrate diversity because insects are easy to sample, have a diverse diet and are widely distributed. Because of these advantages, the use of iDNA for detecting mammals has gained increasing attention. Here we aimed to compare the effectiveness of mosquitoes and flies to detect mammals with a small sampling effort in a semi-controlled area, a zoo that houses native and non-native species. We compared mosquitoes and flies regarding the number of mammal species detected, the amount of mammal sequence reads recovered, and the flight distance range for detecting mammals. We also verified if the combination of two mini-barcodes (12SrRNA and 16SrRNA) would perform better than either mini-barcode alone to inform local mammal biodiversity from iDNA. To capture mosquitoes and flies, we distributed insect traps in eight sampling points during 5 days. We identified 43 Operational Taxonomic Units from 10 orders, from the iDNA of 17 mosquitoes and 46 flies. There was no difference in the number of species recovered per individual insect between mosquitoes and flies, but the number of flies captured was higher, resulting in more mammal species recovered by flies. Eight species were recorded exclusively by mosquitoes and 20 by flies, suggesting that using both samplers would allow a more comprehensive screening of the biodiversity. The maximum distance recorded was 337 m for flies and 289 m for mosquitoes, but the average range distance did not differ between insect groups. Our assay proved to be efficient for mammal detection, considering the high number of species detected with a reduced sampling effort.

Coevolutionary Implications of Obligate Commensalism in Sea Turtles: the Case of the Genus Hyachelia Barnard, 1967 (Crustacea, Amphipoda).

Obligate commensalism in the marine environment and its evolutionary role are still poorly understood. Although sea turtles may serve as ideal substrates for epibionts, within amphipods, only the genus Hyachelia evolved in obligate commensalism with turtles. Here, we report a new host record for Hyachelia lowryi on the hawksbill turtle and describe a larger distribution of the genus in the Atlantic Ocean on green and loggerhead turtles. Hyachelia spp. were sampled from nesting sites of Caretta caretta and feeding grounds of Eretmochelys imbricata and Chelonia mydas along the Brazilian coast. Insights regarding the coevolution of this remarkable genus with its hosts based on molecular analyses are inferred based on mitochondrial (COI) and nuclear (18SrRNA) genes using new and previously available sequences from the infraorder Talitrida. Divergence times for Hyachelia are around the Cretaceous (~127.66 Mya), corresponding to an ancient origin and in agreement with modern green turtle (Chelonioidea) radiation. Later, diversification of Hyachelia species is dated at about 26 Mya, suggesting a coevolutionary association between amphipods and Carettini/Chelonini sea turtles.