The role of isolation on contrasting phylogeographic patterns in two cave crustaceans.

BACKGROUND: The underlying mechanisms and processes that prompt the colonisation of extreme environments, such as caves, constitute major research themes of evolutionary biology and biospeleology. The special adaptations required to survive in subterranean environments (low food availability, hypoxic waters, permanent darkness), and the geographical isolation of caves, nominate cave biodiversity as ideal subjects to answer long-standing questions concerning the interplay amongst adaptation, biogeography, and evolution. The present project aims to examine the phylogeographic patterns exhibited by two sympatric species of surface and cave-dwelling peracarid crustaceans (Asellus aquaticus and Niphargus hrabei), and in doing so elucidate the possible roles of isolation and exaptation in the colonisation and successful adaptation to the cave environment. RESULTS: Specimens of both species were sampled from freshwater hypogean (cave) and epigean (surface) habitats in Hungary, and additional data from neighbouring countries were sourced from Genbank. Sequencing of mitochondrial and nuclear loci revealed, through haplotype network reconstruction (TCS) and phylogenetic inference, the genetic structure, phylogeographic patterns, and divergence-time estimates of A. aquaticus and N. hrabei surface and cave populations. Contrasting phylogeographic patterns were found between species, with A. aquaticus showing strong genetic differentiation between cave and surface populations and N. hrabei lacking any evidence of genetic structure mediated by the cave environment. Furthermore, N. hrabei populations show very low levels of genetic differentiation throughout their range, which suggests the possibility of recent expansion events over the last few thousand years. CONCLUSIONS: Isolation by cave environment, rather than distance, is likely to drive the genetic structuring observed between immediately adjacent cave and surface populations of A. aquaticus, a predominantly surface species with only moderate exaptations to subterranean life. For N. hrabei, in which populations exhibit a fully 'cave-adapted' (troglomorphic) phenotype, the lack of genetic structure suggests that subterranean environments do not pose a dispersal barrier for this surface-cave species.

Phylogenetic and transcriptomic analyses reveal the evolution of bioluminescence and light detection in marine deep-sea shrimps of the family Oplophoridae (Crustacea: Decapoda).

Bioluminescence is essential to the survival of many organisms, particularly in the deep sea where light is limited. Shrimp of the family Oplophoridae exhibit a remarkable mechanism of bioluminescence in the form of a secretion used for predatory defense. Three of the ten genera possess an additional mode of bioluminescence in the form of light-emitting organs called photophores. Phylogenetic analyses can be useful for tracing the evolution of bioluminescence, however, the few studies that have attempted to reconcile the relationships within Oplophoridae have generated trees with low-resolution. We present the most comprehensive phylogeny of Oplophoridae to date, with 90% genera coverage using seven genes (mitochondrial and nuclear) across 30 oplophorid species. We use our resulting topology to trace the evolution of bioluminescence within Oplophoridae. Previous studies have suggested that oplophorid visual systems may be tuned to differentiate the separate modes of bioluminescence. While all oplophorid shrimp possess a visual pigment sensitive to blue-green light, only those bearing photophores have an additional pigment sensitive to near-ultraviolet light. We attempt to characterize opsins, visual pigment proteins essential to light detection, in two photophore-bearing species (Systellaspis debilis and Oplophorus gracilirostris) and make inferences regarding their function and evolutionary significance.

Phylogenetic and transcriptomic characterization of insulin and growth factor receptor tyrosine kinases in crustaceans.

Receptor tyrosine kinases (RTKs) mediate the actions of growth factors in metazoans. In decapod crustaceans, RTKs are implicated in various physiological processes, such molting and growth, limb regeneration, reproduction and sexual differentiation, and innate immunity. RTKs are organized into two main types: insulin receptors (InsRs) and growth factor receptors, which include epidermal growth factor receptor (EGFR), fibroblast growth factor receptor (FGFR), vascular endothelial growth factor receptor (VEGFR), and platelet-derived growth factor receptor (PDGFR). The identities of crustacean RTK genes are incomplete. A phylogenetic analysis of the CrusTome transcriptome database, which included all major crustacean taxa, showed that RTK sequences segregated into receptor clades representing InsR (72 sequences), EGFR (228 sequences), FGFR (129 sequences), and PDGFR/VEGFR (PVR; 235 sequences). These four receptor families were distinguished by the domain organization of the extracellular N-terminal region and motif sequences in the protein kinase catalytic domain in the C-terminus or the ligand-binding domain in the N-terminus. EGFR1 formed a single monophyletic group, while the other RTK sequences were divided into subclades, designated InsR1-3, FGFR1-3, and PVR1-2. In decapods, isoforms within the RTK subclades were common. InsRs were characterized by leucine-rich repeat, furin-like cysteine-rich, and fibronectin type 3 domains in the N-terminus. EGFRs had leucine-rich repeat, furin-like cysteine-rich, and growth factor IV domains. N-terminal regions of FGFR1 had one to three immunoglobulin-like domains, whereas FGFR2 had a cadherin tandem repeat domain. PVRs had between two and five immunoglobulin-like domains. A classification nomenclature of the four RTK classes, based on phylogenetic analysis and multiple sequence alignments, is proposed.

CrusTome: a transcriptome database resource for large-scale analyses across Crustacea.

Transcriptomes from nontraditional model organisms often harbor a wealth of unexplored data. Examining these data sets can lead to clarity and novel insights in traditional systems, as well as to discoveries across a multitude of fields. Despite significant advances in DNA sequencing technologies and in their adoption, access to genomic and transcriptomic resources for nontraditional model organisms remains limited. Crustaceans, for example, being among the most numerous, diverse, and widely distributed taxa on the planet, often serve as excellent systems to address ecological, evolutionary, and organismal questions. While they are ubiquitously present across environments, and of economic and food security importance, they remain severely underrepresented in publicly available sequence databases. Here, we present CrusTome, a multispecies, multitissue, transcriptome database of 201 assembled mRNA transcriptomes (189 crustaceans, 30 of which were previously unpublished, and 12 ecdysozoans for phylogenetic context) as an evolving and publicly available resource. This database is suitable for evolutionary, ecological, and functional studies that employ genomic/transcriptomic techniques and data sets. CrusTome is presented in BLAST and DIAMOND formats, providing robust data sets for sequence similarity searches, orthology assignments, phylogenetic inference, etc. and thus allowing for straightforward incorporation into existing custom pipelines for high-throughput analyses. In addition, to illustrate the use and potential of CrusTome, we conducted phylogenetic analyses elucidating the identity and evolution of the cryptochrome/photolyase family of proteins across crustaceans.

In silico analysis of crustacean hyperglycemic hormone family G protein-coupled receptor candidates.

Ecdysteroid molting hormone synthesis is directed by a pair of molting glands or Y-organs (YOs), and this synthesis is inhibited by molt-inhibiting hormone (MIH). MIH is a member of the crustacean hyperglycemic hormone (CHH) neuropeptide superfamily, which includes CHH and insect ion transport peptide (ITP). It is hypothesized that the MIH receptor is a Class A (Rhodopsin-like) G protein-coupled receptor (GPCR). The YO of the blackback land crab, Gecarcinus lateralis, expresses 49 Class A GPCRs, three of which (Gl-CHHR-A9, -A10, and -A12) were provisionally assigned as CHH-like receptors. CrusTome, a transcriptome database assembled from 189 crustaceans and 12 ecdysozoan outgroups, was used to deorphanize candidate MIH/CHH GPCRs, relying on sequence homology to three functionally characterized ITP receptors (BNGR-A2, BNGR-A24, and BNGR-A34) in the silk moth, Bombyx mori. Phylogenetic analysis and multiple sequence alignments across major taxonomic groups revealed extensive expansion and diversification of crustacean A2, A24, and A34 receptors, designated CHH Family Receptor Candidates (CFRCs). The A2 clade was divided into three subclades; A24 clade was divided into five subclades; and A34 was divided into six subclades. The subclades were distinguished by conserved motifs in extracellular loop (ECL) 2 and ECL3 in the ligand-binding region. Eleven of the 14 subclades occurred in decapod crustaceans. In G. lateralis, seven CFRC sequences, designated Gl-CFRC-A2α1, -A24α, -A24ß1, -A24ß2, -A34α2, -A34ß1, and -A34ß2, were identified; the three A34 sequences corresponded to Gl-GPCR-A12, -A9, and A10, respectively. ECL2 in all the CFRC sequences had a two-stranded ß-sheet structure similar to human Class A GPCRs, whereas the ECL2 of decapod CFRC-A34ß1/ß2 had an additional two-stranded ß-sheet. We hypothesize that this second ß-sheet on ECL2 plays a role in MIH/CHH binding and activation, which will be investigated further with functional assays.

A journey on plate tectonics sheds light on European crayfish phylogeography.

Crayfish can be used as model organisms in phylogeographic and divergence time studies if reliable calibrations are available. This study presents a comprehensive investigation into the phylogeography of the European stone crayfish (Austropotamobius torrentium) and includes samples from previously unstudied sites. Two mitochondrial markers were used to reveal evolutionary relationships among haplogroups throughout the species' distributional range and to estimate the divergence time by employing both substitution rates and geological calibration methods. Our haplotype network reconstruction and phylogenetic analyses revealed the existence of a previously unknown haplogroup distributed in Romania's Apuseni Mountains. This haplogroup is closely related to others that are endemic in the Dinarides, despite their vast geographical separation (~600 km). The separation is best explained by the well-dated tectonic displacement of the Tisza-Dacia microplate, which started in the Miocene (~16 Ma) and possibly carried part of the A. torrentium population to the current location of the Apuseni Mountains. This population may thus have been isolated from the Dinarides for a period of ca. 11 m.y. by marine and lacustrine phases of the Pannonian Basin. The inclusion of this geological event as a calibration point in divergence time analyses challenges currently accepted crayfish evolutionary time frames for the region, constraining the evolution of this area's crayfish to a much earlier date. We discuss why molecular clock calibrations previously employed to date European crayfish species divergences should therefore be reconsidered.

Transcriptomic Insights into the Loss of Vision in Molnár János Cave's Crustaceans.

Animals that inhabit subterranean environments often undergo various distinct phenotypic modifications (referred to as "troglomorphy") as they transition to life in perpetual darkness. However, the molecular basis behind troglomorphy remains poorly understood, particularly in regards to the mechanisms involved in the reduction and/or loss of traits at the transcriptomic level. In this study, we investigate the transcriptional basis behind vision loss in populations of cave-dwelling crustaceans. We employ phylogenetic and transcriptomic methods on surface and cave-adapted populations of an emerging model species for biospeleology, the isopod Asellus aquaticus (Linnaeus, 1758), and the amphipod Niphargus hrabei S. Karaman, 1932. These two species show contrasting directionality in the surface-cave transition, which positions them as ideal study subjects. Asellus aquaticus is common in surface waters and is only occasionally found in caves, where its populations present different degrees of eye reduction and pigmentation. On the other hand, the eyeless N. hrabei has successfully colonized surface environments despite belonging to an almost exclusively cave-dwelling genus. By sequencing and assembling robust de novo transcriptomes we characterized differences in visual genes and pathways among surface and cave populations of the aforementioned species. Our results indicate that despite having reduced eyes, recent cave colonizer A. aquaticus is still capable of expressing functional visual opsins and major components of the phototransduction pathway within the cave. Niphargus hrabei, a species with an ancient cave origin, shows no clear indication of being capable of sight. However, the expression of putative functional visual opsins and other phototransduction genes was maintained, which suggests that this eyeless species might be capable of extraocular photoreception. With the present study, we aim to bring forth the Molnár János Cave system as a promising research avenue to improve our understanding of patterns of reduction and loss of vision in caves and other aphotic environments.

RNA profile diversity across arthropoda: guidelines, methodological artifacts, and expected outcomes.

High-quality RNA is an important precursor for high-throughput RNA sequencing (RNAseq) and subsequent analyses. However, the primary metric used to assess RNA quality, the RNA Integrity Number (RIN), was developed based on model bacterial and vertebrate organisms. Though the phenomenon is not widely recognized, invertebrate 28S ribosomal RNA (rRNA) is highly prone to a form of denaturation known as gap deletion, in which the subunit collapses into two smaller fragments. In many nonmodel invertebrates, this collapse of the 28S subunit appears as a single band similar in size to the 18S rRNA subunit. This phenomenon is hypothesized to be commonplace among arthropods and is often misinterpreted as a "degraded" rRNA profile. The limited characterization of gap deletion in arthropods, a highly diverse group, as well as other nonmodel invertebrates, often biases RNA quality assessments. To test whether the collapse of 28S is a general pattern or a methodological artifact, we sampled more than half of the major lineages within Arthropoda. We found that the 28S collapse is present in ∼90% of the species sampled. Nevertheless, RNA profiles exhibit considerable diversity with a range of banding patterns. High-throughput RNAseq and subsequent assembly of high-quality transcriptomes from select arthropod species exhibiting collapsed 28S subunits further illustrates the limitations of current RIN proxies in accurately characterizing RNA quality in nonmodel organisms. Furthermore, we show that this form of 28S denaturation, which is often mistaken for true "degradation," can occur at relatively low temperatures.