Chromosome-level genome assembly of the coastal horseshoe crab (Tachypleus gigas)

Saving the overlooked mangrove horseshoe crabs-A perspective from enhancing mangrove ecosystem conservation

Despite being widely distributed in Asia, Carcinoscorpius rotundicauda is often overlooked and, its population status remains unclear. Moreover, it is threatened by illegal harvesting and degradation of mangrove ecosystems. Protecting its habitat is essential for population and biodiversity conservation, as mangroves provide nursery grounds and food supply for C. rotundicauda. This review discusses the biological characteristics of C. rotundicauda, including ecology, nutrition, life history, toxicology, and immunology. It also presents information about its distribution and population status. The review emphasizes the challenges faced by C. rotundicauda and proposes a conservation framework that involves the participation of local residents to facilitate conservation efforts. Collaboration between local residents and communities is proposed to protect and monitor the mangrove ecosystem. Additionally, this framework can support field research, protect C. rotundicauda juveniles and other species, and ensure the livelihood of local residents through participation in carbon trading markets and eco-industries such as eco-farming and eco-tourism.

Assessing the relative performance of fast molecular dating methods for phylogenomic data

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

Weighted gene co-expression network analysis of embryos and first instar larvae of the horseshoe crab Tachypleus tridentatus uncovers development gene networks

Horseshoe crabs are marine chelicerates that have existed on Earth for about 450 million years, and they are often used as an experimental model for studying marine invertebrate embryology. In this study, we performed transcriptome gene expression profiling of four continuous embryonic stages (Stages 18-21) and first instar larvae of Tachypleus tridentatus. A mean of 50,742,995 high-quality clean reads was obtained from each library. We then conducted weighted gene co-expression network analysis (WGCNA) for 13,698 genes with fragments per kilobase of exon per million mapped fragments values >5. We identified 17 modules, six of which likely play critical roles in development. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differentially expressed genes was performed on the biologically significant modules. We found that several pathways, such as hedgehog signaling pathway, VEGF signaling pathway, dorso-ventral axis formation, may be involved in the embryonic development process of T. tridentatus. We also identified hub genes that were highly connected in the six critical modules. This is the first study to apply WGCNA to horseshoe crabs to identify hub genes that may play critical roles in development, and our results provide new insight into the mechanisms underlying early development in horseshoe crabs.

Eggs to long-legs: embryonic staging of the harvestman Phalangium opilio (Opiliones), an emerging model arachnid

BACKGROUND

The comparative embryology of Chelicerata has greatly advanced in recent years with the integration of classical studies and genetics, prominently spearheaded by developmental genetic works in spiders. Nonetheless, the understanding of the evolution of development and polarization of embryological characters in Chelicerata is presently limited, as few non-spider species have been well studied. A promising focal species for chelicerate evo-devo is the daddy-long-legs (harvestman) Phalangium opilio, a member of the order Opiliones. Phalangium opilio, breeds prolifically and is easily accessible in many parts of the world, as well as tractable in a laboratory setting. Resources for this species include developmental transcriptomes, a draft genome, and protocols for RNA interference, but a modern staging system is critically missing for this emerging model system.

RESULTS

We present a staging system of P. opilio embryogenesis that spans the most important morphogenetic events with respect to segment formation, appendage elongation and head development. Using time-lapse imaging, confocal microscopy, colorimetric in situ hybridization, and immunohistochemistry, we tracked the development of synchronous clutches from egg laying to adulthood. We describe key events in segmentation, myogenesis, neurogenesis, and germ cell formation.

CONCLUSION

Considering the phylogenetic position of Opiliones and the unduplicated condition of its genome (in contrast to groups like spiders and scorpions), this species is poised to serve as a linchpin for comparative studies in arthropod development and genome evolution. The staging system presented herein provides a valuable reference for P. opilio that we anticipate being useful to the arthropod evo-devo community, with the goal of revitalizing research in the comparative development of non-spider arachnids.

Assessing species coverage and assembly quality of rapidly accumulating sequenced genomes

BACKGROUND

Ambitious initiatives to coordinate genome sequencing of Earth's biodiversity mean that the accumulation of genomic data is growing rapidly. In addition to cataloguing biodiversity, these data provide the basis for understanding biological function and evolution. Accurate and complete genome assemblies offer a comprehensive and reliable foundation upon which to advance our understanding of organismal biology at genetic, species, and ecosystem levels. However, ever-changing sequencing technologies and analysis methods mean that available data are often heterogeneous in quality. To guide forthcoming genome generation efforts and promote efficient prioritization of resources, it is thus essential to define and monitor taxonomic coverage and quality of the data.

FINDINGS

Here we present an automated analysis workflow that surveys genome assemblies from the United States NCBI, assesses their completeness using the relevant BUSCO datasets, and collates the results into an interactively browsable resource. We apply our workflow to produce a community resource of available assemblies from the phylum Arthropoda, the Arthropoda Assembly Assessment Catalogue. Using this resource, we survey current taxonomic coverage and assembly quality at the NCBI, examine how key assembly metrics relate to gene content completeness, and compare results from using different BUSCO lineage datasets.

CONCLUSIONS

These results demonstrate how the workflow can be used to build a community resource that enables large-scale assessments to survey species coverage and data quality of available genome assemblies, and to guide prioritizations for ongoing and future sampling, sequencing, and genome generation initiatives.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

What Is an “Arachnid”? Consensus, Consilience, and Confirmation Bias in the Phylogenetics of Chelicerata

The basal phylogeny of Chelicerata is one of the opaquest parts of the animal Tree of Life, defying resolution despite application of thousands of loci and millions of sites. At the forefront of the debate over chelicerate relationships is the monophyly of Arachnida, which has been refuted by most analyses of molecular sequence data. A number of phylogenomic datasets have suggested that Xiphosura (horseshoe crabs) are derived arachnids, refuting the traditional understanding of arachnid monophyly. This result is regarded as controversial, not least by paleontologists and morphologists, due to the widespread perception that arachnid monophyly is unambiguously supported by morphological data. Moreover, some molecular datasets have been able to recover arachnid monophyly, galvanizing the belief that any result that challenges arachnid monophyly is artefactual. Here, we explore the problems of distinguishing phylogenetic signal from noise through a series of in silico experiments, focusing on datasets that have recently supported arachnid monophyly. We assess the claim that filtering by saturation rate is a valid criterion for recovering Arachnida. We demonstrate that neither saturation rate, nor the ability to assemble a molecular phylogenetic dataset supporting a given outcome with maximal nodal support, is a guarantor of phylogenetic accuracy. Separately, we review empirical morphological phylogenetic datasets to examine characters supporting Arachnida and the downstream implication of a single colonization of terrestrial habitats. We show that morphological support of arachnid monophyly is contingent upon a small number of ambiguous or incorrectly coded characters, most of these tautologically linked to adaptation to terrestrial habitats.

The genome of a daddy-long-legs (Opiliones) illuminates the evolution of arachnid appendages

Chelicerate arthropods exhibit dynamic genome evolution, with ancient whole-genome duplication (WGD) events affecting several orders. Yet, genomes remain unavailable for a number of poorly studied orders, such as Opiliones (daddy-long-legs), which has hindered comparative study. We assembled the first harvestman draft genome for the species Phalangium opilio, which bears elongate, prehensile appendages, made possible by numerous distal articles called tarsomeres. Here, we show that the genome of P. opilio exhibits a single Hox cluster and no evidence of WGD. To investigate the developmental genetic basis for the quintessential trait of this group-the elongate legs-we interrogated the function of the Hox genes Deformed (Dfd) and Sex combs reduced (Scr), and a homologue of Epidermal growth factor receptor (Egfr). Knockdown of Dfd incurred homeotic transformation of two pairs of legs into pedipalps, with dramatic shortening of leg segments in the longest leg pair, whereas homeosis in L3 is only achieved upon double Dfd + Scr knockdown. Knockdown of Egfr incurred shortened appendages and the loss of tarsomeres. The similarity of Egfr loss-of-function phenotypic spectra in insects and this arachnid suggest that repeated cooption of EGFR signalling underlies the independent gains of supernumerary tarsomeres across the arthropod tree of life.

The different fates of two Asian horseshoe crab species with different dispersal abilities

Impending anthropogenic climate change will severely impact coastal organisms at unprecedented speed. Knowledge on organisms' evolutionary responses to past sea-level fluctuations and estimation of their evolutionary potential is therefore indispensable in efforts to mitigate the effects of future climate change. We sampled tens of thousands of genomic markers of ~300 individuals in two of the four extant horseshoe crab species across the complex archipelagic Singapore Straits. Carcinoscorpius rotundicauda Latreille, a less mobile mangrove species, has finer population structure and lower genetic diversity compared with the dispersive deep-sea Tachypleus gigas Müller. Even though the source populations of both species during the last glacial maximum exhibited comparable effective population sizes, the less dispersive C. rotundicauda seems to lose genetic diversity much more quickly because of population fragmentation. Contra previous studies' results, we predict that the more commonly sighted C. rotundicauda faces a more uncertain conservation plight, with a continuing loss in evolutionary potential and higher vulnerability to future climate change. Our study provides important genomic baseline data for the redirection of conservation measures in the face of climate change and can be used as a blueprint for assessment and mitigation of the adverse effects of impending sea-level rise in other systems.

The endangered White Sands pupfish (Cyprinodon tularosa) genome reveals low diversity and heterogenous patterns of differentiation

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08 Gb draft genome with a contig N50 of 1.4 Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193 genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts.

Hemolymph proteins: An overview across marine arthropods and molluscs

In this compilation we collect information about the main protein components in hemolymph and stress the continued interest in their study. The reasons for such an attention span several areas of biological, veterinarian and medical applications: from the notions for better dealing with the species - belonging to phylum Arthropoda, subphylum Crustacea, and to phylum Mollusca - of economic interest, to the development of 'marine drugs' from the peptides that, in invertebrates, act as antimicrobial, antifungal, antiprotozoal, and/or antiviral agents. Overall, the topic most often on focus is that of innate immunity operated by classes of pattern-recognition proteins. SIGNIFICANCE: The immune response in invertebrates relies on innate rather than on adaptive/acquired effectors. At a difference from the soluble and membrane-bound immunoglobulins and receptors in vertebrates, the antimicrobial, antifungal, antiprotozoal and/or antiviral agents in invertebrates interact with non-self material by targeting some common (rather than some highly specific) structural motifs. Developing this paradigm into (semi) synthetic pharmaceuticals, possibly optimized through the modeling opportunities offered by computational biochemistry, is one of the lessons today's science may learn from the study of marine invertebrates, and specifically of the proteins and peptides in their hemolymph.