Phylogeographic reconstruction of the marbled crayfish origin

Global invasion and biosecurity risk from the online trade in ornamental crayfish

Trade in undomesticated ornamental animals has rapidly expanded beyond brick-and-mortar retail stores to now include growing numbers of internet marketplaces. The growing volume, diversity, and origins of invasive non-native species in trade challenge already weak national biosecurity policies. Despite widespread focus on vertebrates, many knowledge gaps exist regarding the online global trade of ornamental invertebrates. We conducted the first global assessment of the online trade in and associated invasion risk of freshwater crayfishes, which are increasingly popular aquarium animals. We systematically examined e-commerce marketplaces in multiple languages, scrapping information that included species identity, price, quantity, and shipping designation from each crayfish listing. Next, we combined geographic pathways of shipping associated with online trade (introduction risk) and environmental suitability modeling (establishment risk) to estimate global risk of non-native crayfish invasion risk. We identified hundreds of online marketplaces and thousands of sale listings in 33 countries (5 continents) involving 60 species and representing a selling value of ∼US$1.5 million. Invasion risk of non-native crayfish in trade was widespread, with geographic hotspots coinciding with both elevated opportunities for introduction (greater shipping offerings) and establishment. Precise characterization of the online species trade is fundamental to support new and reformed biosecurity policies, build industry partnerships, and design educational campaigns to prevent species invasions through trade. We found that the taxonomy, geography, and economics of the global online ornamental crayfish trade are vast and require greater attention.

Genomics resources for the Rapa Nui (Eastern Island) spiny lobster Panulirus pascuensis (Crustacea: Decapoda: Achelata)

Background

The Easter Island spiny lobster Panulirus pascuensis (Reed, 1954) or ‘Ura’ in the Rapa Nui language, is a little known species native to the south eastern Pacific Ocean, distributed along the coasts of Easter Island, Pitcairn Island, and the Salas y Gómez Ridge. In Easter Island, P. pascuensis is the target of a small and profitable and probably overexploited fishery. In this study, we profited from a series of bioinformatic analyses to mine biological insight from low-pass short-read next generation sequencing datasets; we have estimated genome size and ploidy in P. pascuensis using a k-mer strategy, discovered, annotated, and quantified mobile elements in the nuclear genome, assembled the 45S rRNA nuclear DNA cassette and mitochondrial chromosome, and explored the phylogenetic position of P. pascuensis within the genus Panulirus using the signal retrieved from translated mitochondrial protein coding genes.

Results

K-mer analyses predicted P. pascuensis to be diploid with a haploid genome size ranging between 2.75 Gbp (with k-mer = 51) and 3.39 Gbp (with k-mer = 18). In P. pascuensis, repetitive elements comprise at least a half and a maximum of three fourths of the nuclear genome. Almost a third (64.94%) of the repetitive elements present in the studied nuclear genome were not assigned to any known family of transposable elements. Taking into consideration only annotated repetitive elements, the most abundant were classified as Long Interspersed Nuclear Elements (22.81%). Less common repetitive elements included Long Terminal Repeats (2.88%), Satellite DNA (2.66%), and DNA transposons (2.45%), among a few others. The 45S rRNA DNA cassette of P. pascuensis was partially assembled into two contigs. One contig, 2,226 bp long, encoded a partially assembled 5′ ETS the entire ssrDNA (1,861 bp), and a partial ITS1. A second contig, 6,714 bp long, encoded a partially assembled ITS1, the entire 5.8S rDNA (158 bp), the entire ITS2, the entire lsrDNA (4,938 bp), and a partial 3′ ETS (549 bp). The mitochondrial genome of P. pascuensis was 15,613 bp long and contained 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and two ribosomal RNA genes (12S ribosomal RNA [rrnS] and 16S ribosomal RNA [rrnL]). A phylomitogenomic analysis based on PCGs retrieved Panulirus pascuensis as sister to a fully supported clade comprising P. cygnus and P. longipes.

Conclusion

We expect that the information generated in this study will guide the assembly of a chromosome-level nuclear genome for P. pascuensis in the near future. The newly assembled 45S rRNA nuclear DNA cassette and mitochondrial chromosome can support bioprospecting and biomonitoring of P. pascuensis using environmental DNA. The same elements can help to survey the public market place and detect mislabelling of this and other spiny lobsters. Overall, the genomic resources generated in this study will aid in supporting fisheries management and conservation strategies in this iconic spiny lobster that is likely experiencing overexploitation.

Considerations for protein and amino acids in standardized reference diet for parthenogenetic marbled crayfish Procambarus virginalis model organism

The concept of a standardized reference diet (SRD) is used in laboratory model organisms to ensure nutritional control between studies and laboratories. Although models using the genetically identical, all female parthenogenetic marbled crayfish (Procambarus virginalis) are growing in popularity, research into nutrition in this species still has many knowledge gaps. To fast track the development of a SRD in terms of protein and amino acids (SRDprotein) for this species, we first analyzed the composition of its body amino acids to determine the ideal protein concept (IPC) of indispensable amino acids in wild-caught P. virginalis (which had an unusually high preponderance of leucine and arginine). Then, we strategically evaluated three common clusters of types of fish feed: (1) ornamental fish feed (SER) fortified with a naturally occurring alga (Spirulina). This type of feed was protein-high in arginine and leucine (SER + SPI) that fulfils the species' IPC for iso-protein (~ 40%), iso-phosphorus (~ 0.8%) and near iso-energetic (~ 475 kcal 100 g-1); (2) freeze-dried live feed consisting of chironomid larvae (CHI) fortified with Spirulina (CHI + SPI) that fulfils the IPC for iso-protein (~ 46%), iso-phosphorus (~ 0.7%) and near iso-energetic (~ 405 kcal 100 g-1); and (3) a commercially standardized 'starter diet' for carnivorous fish larvae (FISH) and post-larval shrimps (SHRIMP) with iso-protein (~ 56%) and iso-phosphorus (~ 1.6%). A total of six diets, embracing a diverse range of proteinaceous feeds, were used in a 100-day ad libitum feeding and growth trial. The FISH group outperformed all the other groups (p < 0.05) and our exploratory multivariate analysis revealed an ideal demand of > 44% protein (tailored to deliver high arginine 3% and leucine 4%, followed by the usual lysine > 3.5% and methionine 1.2%) but also the lowest carbohydrate level (21%). For SRDprotein, our findings show that the FISH diet is ideal and suggest the possibilities of using a CHI + SPI diet for further optimization (more economic use of protein and phosphorus).

Considerations for fatty acids in standardized reference diet for parthenogenetic marbled crayfish Procambarus virginalis model organism

Fatty acid accumulation was studied in the parthenogenetic all-female marbled crayfish Procambarus virginalis using six arbitrarily designed experimental feeds and related to individuals with glair glands (sexual maturity) after 100 days of ad libitum feeding at 21 °C, including gravid females from the wild as a reference. Fatty acids 16:0 and 18:1n-9 comprised 40% of the total amount of fatty acids and tended to up-concentrate in bodies. Shorter chain 14:0 depleted from feed to body. Across diets, there was a concomitant decrease in precursor fatty acid and increase in product fatty acid, such as reinforcements in monounsaturated fatty acid (18:1n-9), eicosanoid precursors 20:4n-6 (arachidonic acid, ARA) and 20:5n-3 (eicosapentaenoic acid, EPA) in-vivo, but not 22:6n-3 (docosahexaenoic acid, DHA) except when deficient in CHI or CHI + SPI diets. Saturation kinetics modeling (R2 0.7-0.9, p < 0.05) showed that when the ARA share is ~ 1%, the EPA share is ~ 8%, and the DHA share is ~ 2% in the food lipids, the accumulation of fatty acids in body lipids levels off. The lowest DHA in the CHI (0% glair glands) or CHI + SPI (0-3.9% glair glands) diets, and the lowest ARA in SER (0% glair glands) or SER + SPI (0-3% glair glands) diets, were synchronous with negligible sexual maturity despite a wide range of observed specific growth rates (2.77-3.60% per day), body size (0.44-0.84 g), ≤ 5% crude lipid and 40-46% crude protein feed. The FISH and SHRIMP diets (56% protein, 11-14% lipid) with the highest ARA, EPA, and DHA together seem to be the most conducive diets for sexual maturity (up to 20% of individuals with glair glands). We propose a fatty acid profile mimicking the FISH or SHRIMP diets as a starting point for designing the lipid content required in the marbled crayfish standardized reference diet.

Temperature and hypoxia trigger developmental phenotypic plasticity of cardiorespiratory physiology and growth in the parthenogenetic marbled crayfish, Procambarus virginalis Lyko, 2017

Attempting to differentiate phenotypic variation caused by environmentally-induced alterations in gene expression from that caused by actual allelic differences can be experimentally difficult. Environmental variables must be carefully controlled and then interindividual genetic differences ruled out as sources of phenotypic variation. We investigated phenotypic variability of cardiorespiratory physiology as well as biometric traits in the parthenogenetically-reproducing marbled crayfish Procambarus virginalis Lyko, 2017, all offspring being genetically identical clones. Populations of P. virginalis were reared from eggs tank-bred at four different temperatures (16, 19, 22 and 25 °C) or two different oxygen levels (9.5 and 20 kPa). Then, at Stage 3 and 4 juvenile stages, physiological (heart rate, oxygen consumption) and morphological (carapace length, body mass) variables were measured. Heart rate and oxygen consumption measured at 23 °C showed only small effects of rearing temperature in Stage 3 juveniles, with larger effects evident in older, Stage 4 juveniles. Additionally, coefficients of variation were calculated to compare our data to previously published data on P. virginalis as well as sexually-reproducing crayfish. Comparison revealed that carapace length, body mass and heart rate (but not oxygen consumption) indeed showed lower, yet notable coefficients of variation in clonal crayfish. Yet, despite being genetically identical, significant variation in their morphology and physiology in response to different rearing conditions nonetheless occurred in marbled crayfish. This suggests that epigenetically induced phenotypic variation might play a significant role in asexual but also sexually reproducing species.

Time-resolved, integrated analysis of clonally evolving genomes

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

Temperature affects predation of schistosome-competent snails by a novel invader, the marbled crayfish Procambarus virginalis

The human burden of environmentally transmitted infectious diseases can depend strongly on ecological factors, including the presence or absence of natural enemies. The marbled crayfish (Procambarus virginalis) is a novel invasive species that can tolerate a wide range of ecological conditions and colonize diverse habitats. Marbled crayfish first appeared in Madagascar in 2005 and quickly spread across the country, overlapping with the distribution of freshwater snails that serve as the intermediate host of schistosomiasis-a parasitic disease of poverty with human prevalence ranging up to 94% in Madagascar. It has been hypothesized that the marbled crayfish may serve as a predator of schistosome-competent snails in areas where native predators cannot and yet no systematic study to date has been conducted to estimate its predation rate on snails. Here, we experimentally assessed marbled crayfish consumption of uninfected and infected schistosome-competent snails (Biomphalaria glabrata and Bulinus truncatus) across a range of temperatures, reflective of the habitat range of the marbled crayfish in Madagascar. We found that the relationship between crayfish consumption and temperature is unimodal with a peak at ~27.5°C. Per-capita consumption increased with body size and was not affected either by snail species or their infectious status. We detected a possible satiation effect, i.e., a small but significant reduction in per-capita consumption rate over the 72-hour duration of the predation experiment. Our results suggest that ecological parameters, such as temperature and crayfish weight, influence rates of consumption and, in turn, the potential impact of the marbled crayfish invasion on snail host populations.

Phenotypic plasticity in the monoclonal marbled crayfish is associated with very low genetic diversity but pronounced epigenetic diversity

Clonal organisms are particularly useful to investigate the contribution of epigenetics to phenotypic plasticity, because confounding effects of genetic variation are negligible. In the last decade, the apomictic parthenogenetic marbled crayfish, Procambarus virginalis, has been developed as a model to investigate the relationships between phenotypic plasticity and genetic and epigenetic diversity in detail. This crayfish originated about 30 years ago by autotriploidy from a single slough crayfish Procambarus fallax. As the result of human releases and active spreading, marbled crayfish has established numerous populations in very diverse habitats in 22 countries from the tropics to cold temperate regions. Studies in the laboratory and field revealed considerable plasticity in coloration, spination, morphometric parameters, growth, food preference, population structure, trophic position, and niche width. Illumina and PacBio whole-genome sequencing of marbled crayfish from representatives of 19 populations in Europe and Madagascar demonstrated extremely low genetic diversity within and among populations, indicating that the observed phenotypic diversity and ability to live in strikingly different environments are not due to adaptation by selection on genetic variation. In contrast, considerable differences were found between populations in the DNA methylation patterns of hundreds of genes, suggesting that the environmentally induced phenotypic plasticity is mediated by epigenetic mechanisms and corresponding changes in gene expression. Specific DNA methylation fingerprints persisted in local populations over successive years indicating the existence of epigenetic ecotypes, but there is presently no information as to whether these epigenetic signatures are transgenerationally inherited or established anew in each generation and whether the recorded phenotypic plasticity is adaptive or nonadaptive.

Environmental Adaptation of Genetically Uniform Organisms with the Help of Epigenetic Mechanisms-An Insightful Perspective on Ecoepigenetics

Organisms adapt to different environments by selection of the most suitable phenotypes from the standing genetic variation or by phenotypic plasticity, the ability of single genotypes to produce different phenotypes in different environments. Because of near genetic identity, asexually reproducing populations are particularly suitable for the investigation of the potential and molecular underpinning of the latter alternative in depth. Recent analyses on the whole-genome scale of differently adapted clonal animals and plants demonstrated that epigenetic mechanisms such as DNA methylation, histone modifications and non-coding RNAs are among the molecular pathways supporting phenotypic plasticity and that epigenetic variation is used to stably adapt to different environments. Case studies revealed habitat-specific epigenetic fingerprints that were maintained over subsequent years pointing at the existence of epigenetic ecotypes. Environmentally induced epimutations and corresponding gene expression changes provide an ideal means for fast and directional adaptation to changing or new conditions, because they can synchronously alter phenotypes in many population members. Because microorganisms inclusive of human pathogens also exploit epigenetically mediated phenotypic variation for environmental adaptation, this phenomenon is considered a universal biological principle. The production of different phenotypes from the same DNA sequence in response to environmental cues by epigenetic mechanisms also provides a mechanistic explanation for the "general-purpose genotype hypothesis" and the "genetic paradox of invasions".

Host-pathogen coevolution drives innate immune response to Aphanomyces astaci infection in freshwater crayfish: transcriptomic evidence

BACKGROUND

For over a century, scientists have studied host-pathogen interactions between the crayfish plague disease agent Aphanomyces astaci and freshwater crayfish. It has been hypothesised that North American crayfish hosts are disease-resistant due to the long-lasting coevolution with the pathogen. Similarly, the increasing number of latent infections reported in the historically sensitive European crayfish hosts seems to indicate that similar coevolutionary processes are occurring between European crayfish and A. astaci. Our current understanding of these host-pathogen interactions is largely focused on the innate immunity processes in the crayfish haemolymph and cuticle, but the molecular basis of the observed disease-resistance and susceptibility remain unclear. To understand how coevolution is shaping the host's molecular response to the pathogen, susceptible native European noble crayfish and invasive disease-resistant marbled crayfish were challenged with two A. astaci strains of different origin: a haplogroup A strain (introduced to Europe at least 50 years ago, low virulence) and a haplogroup B strain (signal crayfish in lake Tahoe, USA, high virulence). Here, we compare the gene expression profiles of the hepatopancreas, an integrated organ of crayfish immunity and metabolism.

RESULTS

We characterised several novel innate immune-related gene groups in both crayfish species. Across all challenge groups, we detected 412 differentially expressed genes (DEGs) in the noble crayfish, and 257 DEGs in the marbled crayfish. In the noble crayfish, a clear immune response was detected to the haplogroup B strain, but not to the haplogroup A strain. In contrast, in the marbled crayfish we detected an immune response to the haplogroup A strain, but not to the haplogroup B strain.

CONCLUSIONS

We highlight the hepatopancreas as an important hub for the synthesis of immune molecules in the response to A. astaci. A clear distinction between the innate immune response in the marbled crayfish and the noble crayfish is the capability of the marbled crayfish to mobilise a higher variety of innate immune response effectors. With this study we outline that the type and strength of the host immune response to the pathogen is strongly influenced by the coevolutionary history of the crayfish with specific A. astaci strains.

Studying phenotypic variation and DNA methylation across development, ecology and evolution in the clonal marbled crayfish: a paradigm for investigating epigenotype-phenotype relationships in macro-invertebrates

Animals can produce different phenotypes from the same genome during development, environmental adaptation and evolution, which is mediated by epigenetic mechanisms including DNA methylation. The obligatory parthenogenetic marbled crayfish, Procambarus virginalis, whose genome and methylome are fully established, proved very suitable to study this issue in detail. Comparison between developmental stages and DNA methylation revealed low expression of Dnmt methylation and Tet demethylation enzymes from the spawned oocyte to the 256 cell embryo and considerably increased expression thereafter. The global 5-methylcytosine level was 2.78% at mid-embryonic development and decreased slightly to 2.41% in 2-year-old adults. Genetically identical clutch-mates raised in the same uniform laboratory setting showed broad variation in morphological, behavioural and life history traits and differences in DNA methylation. The invasion of diverse habitats in tropical to cold-temperate biomes in the last 20 years by the marbled crayfish was associated with the expression of significantly different phenotypic traits and DNA methylation patterns, despite extremely low genetic variation on the whole genome scale, suggesting the establishment of epigenetic ecotypes. The evolution of marbled crayfish from its parent species Procambarus fallax by autotriploidy a few decades ago was accompanied by a significant increase in body size, fertility and life span, a 20% reduction of global DNA methylation and alteration of methylation in hundreds of genes, suggesting that epigenetic mechanisms were involved in speciation and fitness enhancement. The combined analysis of phenotypic traits and DNA methylation across multiple biological contexts in the laboratory and field in marbled crayfish may serve as a blueprint for uncovering the role of epigenetic mechanisms in shaping of phenotypes in macro-invertebrates.