Unveiling the Genetic Blueprint of a Desert Scorpion: A Chromosome-level Genome of Hadrurus arizonensis Provides the First Reference for Parvorder Iurida.

Over 400 million years old, scorpions represent an ancient group of arachnids and one of the first animals to adapt to life on land. Presently, the lack of available genomes within scorpions hinders research on their evolution. This study leverages ultralong nanopore sequencing and Pore-C to generate the first chromosome-level assembly and annotation for the desert hairy scorpion, Hadrurus arizonensis. The assembled genome is 2.23 Gb in size with an N50 of 280 Mb. Pore-C scaffolding reoriented 99.6% of bases into nine chromosomes and BUSCO identified 998 (98.6%) complete arthropod single copy orthologs. Repetitive elements represent 54.69% of the assembled bases, including 872,874 (29.39%) LINE elements. A total of 18,996 protein-coding genes and 75,256 transcripts were predicted, and extracted protein sequences yielded a BUSCO score of 97.2%. This is the first genome assembled and annotated within the family Hadruridae, representing a crucial resource for closing gaps in genomic knowledge of scorpions, resolving arachnid phylogeny, and advancing studies in comparative and functional genomics.

Out from under the wing: reconceptualizing the insect wing gene regulatory network as a versatile, general module for body-wall lobes in arthropods.

Body plan evolution often occurs through the differentiation of serially homologous body parts, particularly in the evolution of arthropod body plans. Recently, homeotic transformations resulting from experimental manipulation of gene expression, along with comparative data on the expression and function of genes in the wing regulatory network, have provided a new perspective on an old question in insect evolution: how did the insect wing evolve? We investigated the metamorphic roles of a suite of 10 wing- and body-wall-related genes in a hemimetabolous insect, Oncopeltus fasciatus. Our results indicate that genes involved in wing development in O. fasciatus play similar roles in the development of adult body-wall flattened cuticular evaginations. We found extensive functional similarity between the development of wings and other bilayered evaginations of the body wall. Overall, our results support the existence of a versatile development module for building bilayered cuticular epithelial structures that pre-dates the evolutionary origin of wings. We explore the consequences of reconceptualizing the canonical wing-patterning network as a bilayered body-wall patterning network, including consequences for long-standing debates about wing homology, the origin of wings and the origin of novel bilayered body-wall structures. We conclude by presenting three testable predictions that result from this reconceptualization.

Something old, something new, something borrowed, something red: the origin of ecologically relevant novelties in Hemiptera.

Comparative transcriptomics, applied in an evolutionary context, has transformed the possibilities for studying phenotypic evolution in non-model taxa. We review recent discoveries about the development of novel, ecologically relevant phenotypes in hemipteran insects. These discoveries highlight the diverse genomic substrates of novelty: 'something old', when novelty results from changes in the regulation of existing genes or gene duplication; 'something new', wherein lineage-restricted genes contribute to the evolution of new phenotypes; and 'something borrowed', showcasing contributions of horizontal gene transfer to the evolution of novelty, including carotenoid synthesis (resulting in 'something red'). These findings show the power and flexibility of comparative transcriptomic approaches for expanding beyond the 'toolkit' model for the evolution of development. We conclude by raising questions about the relationship between new genes and new traits and outlining a research framework for answering them in Hemiptera.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

Slender salamanders (genus Batrachoseps) reveal Southern California to be a center for the diversification, persistence, and introduction of salamander lineages.

BACKGROUND: The southern California biodiversity hotspot has had a complex geological history, with both plate tectonic forces and sea level changes repeatedly reconfiguring the region, and likely driving both lineage splittings and extinctions. Here we investigate patterns of genetic divergence in two species of slender salamanders (Plethodontidae: Batrachoseps) in this region. The complex geological history in combination with several organismal traits led us to predict that these species harbor multiple ancient mitochondrial lineages endemic to southern California. These species belong to a clade characterized by fine-scale mitochondrial structure, which has been shown to track ancient splits. Both focal species, Batrachoseps major and B. nigriventris, are relatively widely distributed in southern California, and estimated to have persisted there across millions of years. Recently several extralimital populations of Batrachoseps were found in the San Joaquin Valley of California, a former desert area that has been extensively modified for agriculture. The origins of these populations are unknown, but based on morphology, they are hypothesized to result from human-mediated introductions of B. major. METHODS: We sequenced the mitochondrial gene cytochrome b from a geographically comprehensive sampling of the mitochondrial lineages of B. major and B. nigriventris that are endemic to southern California. We used phylogenetic analyses to characterize phylogeographic structure and identify mitochondrial contact zones. We also included the San Joaquin Valley samples to test whether they resulted from introductions. We used a bootstrap resampling approach to compare the strength of isolation-by-distance in both Batrachoseps species and four other salamander species with which they co-occur in southern California. RESULTS: The northern lineage of B. major harbors at least eight deeply differentiated, geographically cohesive mitochondrial subclades. We identify geographic contact between many of these mtDNA lineages and some biogeographic features that are concordant with lineage boundaries. Batrachoseps nigriventris also has multiple deeply differentiated clades within the region. Comparative analyses highlight the smaller spatial scales over which mitochondrial divergence accumulates in Batrachoseps relative to most other salamander species in southern California. The extralimital populations of Batrachoseps from the San Joaquin Valley are assigned to B. major and are shown to result from at least two independent introductions from different source populations. We also suggest that B. major on Catalina Island, where it is considered native, may be the result of an introduction. Some of the same traits that facilitate the build-up of deep phylogeographic structure in Batrachoseps likely also contribute to its propensity for introductions, and we anticipate that additional introduced populations will be discovered.

Developmental temperature influences color polymorphism but not hatchling size in a woodland salamander.

Phenotypic plasticity can be an important adaptive response to climate change, particularly for dispersal-limited species. Temperature frequently alters developmental and phenotypic traits including morphology, behavior, and reproductive cycles. We often lack crucial information about if and how thermal conditions during development will interact with genetic responses and facilitate persistence or adaptation under climate change. Polymorphic species offer an ideal test for this, as alternative morphs often confer differential adaptive advantages. However, few studies have examined the effects of incubation temperature on color expression or development in polymorphic taxa. Here we test if developmental temperature mediates morph frequency in the polymorphic salamander Plethodon cinereus. Although previous research suggests geographic variation in morph proportions results from differential climate adaptation, it remains unknown if plasticity also contributes to this variation. We used a split-clutch common garden experiment to determine the effects of developmental temperature on the color and development of P. cinereus. Our results indicate developmental temperature affects coloration in P. cinereus, either via plasticity or differential mortality, with eggs incubated at warmer temperatures yielding a higher proportion of unstriped individuals than those from cooler temperatures. This temperature response may contribute to the spatial variation in morph frequencies in natural populations. Surprisingly, we found neither temperature nor egg size affected hatchling size. Our study provides important insights into the potential for climate-induced responses to preserve diversity in dispersal-limited species, like P. cinereus, and enable time for adaptive evolution.

Co-option of wing-patterning genes underlies the evolution of the treehopper helmet.

Understanding the origin of novelty is a key question in evolutionary developmental biology. In arthropods, the body wall has served as a repeated source of morphological novelty. In treehoppers, an ancestrally flat part of the dorsal body wall (the pronotum) was transformed into a three-dimensional structure (the helmet), which was subsequently moulded by natural selection into diverse shapes. Here, we test three hypotheses for the developmental origin of the helmet by comparing body-region transcriptomes in a treehopper and a leafhopper that retains more ancestral morphology. In leafhoppers, pronotal gene expression is most similar to that of its serial homologue, the mesonotum. By contrast, in treehoppers, helmet gene expression is most similar to that of wings, supporting the wing-patterning network co-option hypothesis for the origin of the helmet. These results suggest that serial homologues may diverge evolutionarily through replacement of, rather than tinkering with, their ancestrally shared patterning network.

New insights into the identities of the elasmobranch fauna of Sri Lanka.

Identities of elasmobranchs from Sri Lanka encountered during collections conducted in an intensive nine-day survey of fish markets and landing sites at 11 localities in the North Western, Northern, and Eastern Provinces in March of 2018 were assessed. In total, 111 specimens representing 34 elasmobranch species were examined. Sequence data for the NADH2 gene were generated for all specimens. Independent Neighbor-Joining analyses, which included data for related taxa, were conducted for 25 subgroups of elasmobranchs to help confirm specific identifications. Five of the 34 species encountered are likely new to science. These consist of one species each of the batoid genera Brevitrygon, Narcine, and Torpedo, and the selachian genera Centrophorus, and Chiloscyllium. The specific identities of 12 species previously known to occur in Sri Lanka are updated to conform to current taxonomy; four of these (Gymnura cf. poecilura 2, Carcharhinus cf. limbatus, Echinorhinus sp. 1, and Iago cf. omanensis 1) represent what appear to be undescribed species reported previously from other localities. Three species (Maculabatis arabica, Acroteriobatus variegatus, and Centroscymnus owstonii) are reported from Sri Lanka for the first time; the latter species also represents the first documented record of this genus and family for the island nation. One of the two specimens on which the recent description of the new species of Planonasus indicus was based was also collected as part of this survey. Although some of the species confirmed to occur in Sri Lanka have also been found in India, others were previously known only from the Persian Gulf, Gulf of Oman, or localities in Southeast Asia. The high amount of novelty discovered as a result of a survey of such short duration emphasizes the importance of more intensive survey efforts in this region now that the civil unrest that precluded such work for nearly three decades has come to an end.

Developmental and Evolutionary Perspectives on the Origin and Diversification of Arthropod Appendages.

Jointed, segmented appendages are a key innovation of arthropods. The subsequent diversification of these appendages, both along the body axis and across taxa, has contributed to the evolutionary success of arthropods. Both developmental and fossil data are informative for understanding how these transitions occurred. Comparative analyses help to pinpoint the developmental novelties that distinguish arthropod appendages from the lobopodous appendages of other panarthropods, and that distinguish different appendage types. The fossil record of stem group arthropods is diverse and preserves intermediate steps in these evolutionary transitions, including some that cannot be directly inferred based on extant taxa. These lead to hypotheses that can be tested with comparative developmental data, as well as to reinterpretations of developmental results. One developmental novelty of arthropods is the reiterated deployment of the joint formation network, which divides the appendages into segments. The fossil record raises questions about how this joint formation network was first deployed, given the contrasting morphologies of appendages in stem group versus extant arthropods. The fossil record supports a character tree for appendage diversification showing progressive individuation of appendages in an anterior-to-posterior sequence. However, to date, developmental evidence provides at best limited support for this character tree. Recent interpretations of the fossil record suggest that the labrum of extant arthropods is a greatly reduced protocerebral appendage pair; this hypothesis is consistent with the extensive shared developmental patterning of the labrum and jointed appendages. Reciprocal illumination from fossils and developmental patterning in a phylogenetic context both makes sense of some results and helps motivates questions for future research.

The Compact Body Plan of Tardigrades Evolved by the Loss of a Large Body Region.

The superphylum Panarthropoda (Arthropoda, Onychophora, and Tardigrada) exhibits a remarkable diversity of segment morphologies, enabling these animals to occupy diverse ecological niches. The molecular identities of these segments are specified by Hox genes and other axis patterning genes during development [1, 2]. Comparisons of molecular segment identities between arthropod and onychophoran species have yielded important insights into the origins and diversification of their body plans [3-9]. However, the relationship of the segments of tardigrades to those of arthropods and onychophorans has remained enigmatic [10, 11], limiting our understanding of early panarthropod body plan diversification. Here, we reveal molecular identities for all of the segments of a tardigrade. Based on our analysis, we conclude that tardigrades have lost a large intermediate region of the body axis-a region corresponding to the entire thorax and most of the abdomen of insects-and that they have lost the Hox genes that originally specified this region. Our data suggest that nearly the entire tardigrade body axis is homologous to just the head region of arthropods. Based on our results, we reconstruct a last common ancestor of Panarthropoda that had a relatively elongate body plan like most arthropods and onychophorans, rather than a compact, tardigrade-like body plan. These results demonstrate that the body plan of an animal phylum can originate by the loss of a large part of the body.

The effects of inference method, population sampling, and gene sampling on species tree inferences: an empirical study in slender salamanders (Plethodontidae: Batrachoseps).

Species tree methods are now widely used to infer the relationships among species from multilocus data sets. Many methods have been developed, which differ in whether gene and species trees are estimated simultaneously or sequentially, and in how gene trees are used to infer the species tree. While these methods perform well on simulated data, less is known about what impacts their performance on empirical data. We used a data set including five nuclear genes and one mitochondrial gene for 22 species of Batrachoseps to compare the effects of method of analysis, within-species sampling and gene sampling on species tree inferences. For this data set, the choice of inference method had the largest effect on the species tree topology. Exclusion of individual loci had large effects in *BEAST and STEM, but not in MP-EST. Different loci carried the greatest leverage in these different methods, showing that the causes of their disproportionate effects differ. Even though substantial information was present in the nuclear loci, the mitochondrial gene dominated the *BEAST species tree. This leverage is inherent to the mtDNA locus and results from its high variation and lower assumed ploidy. This mtDNA leverage may be problematic when mtDNA has undergone introgression, as is likely in this data set. By contrast, the leverage of RAG1 in STEM analyses does not reflect properties inherent to the locus, but rather results from a gene tree that is strongly discordant with all others, and is best explained by introgression between distantly related species. Within-species sampling was also important, especially in *BEAST analyses, as shown by differences in tree topology across 100 subsampled data sets. Despite the sensitivity of the species tree methods to multiple factors, five species groups, the relationships among these, and some relationships within them, are generally consistently resolved for Batrachoseps.

Hox genes require homothorax and extradenticle for body wall identity specification but not for appendage identity specification during metamorphosis of Tribolium castaneum.

The establishment of segment identity is a key developmental process that allows for divergence along the anteroposterior body axis in arthropods. In Drosophila, the identity of a segment is determined by the complement of Hox genes it expresses. In many contexts, Hox transcription factors require the protein products of extradenticle (exd) and homothorax (hth) as cofactors to perform their identity specification functions. In holometabolous insects, segment identity may be specified twice, during embryogenesis and metamorphosis. To glean insight into the relationship between embryonic and metamorphic segmental identity specification, we have compared these processes in the flour beetle Tribolium castaneum, which develops ventral appendages during embryogenesis that later metamorphose into adult appendages with distinct morphologies. At metamorphosis, comparisons of RNAi phenotypes indicate that Hox genes function jointly with Tc-hth and Tc-exd to specify several region-specific aspects of the adult body wall. On the other hand, Hox genes specify appendage identities along the anteroposterior axis independently of Tc-hth/Tc-exd and Tc-hth/Tc-exd specify proximal vs. distal identity within appendages independently of Hox genes during this stage. During embryogenesis, Tc-hth and Tc-exd play a broad role in the segmentation process and are required for specification of body wall identities in the thorax; however, contrasting with results from other species, we did not obtain homeotic transformations of embryonic appendages in response to Tc-hth or Tc-exd RNAi. In general, the homeotic effects of interference with the function of Hox genes and Tc-hth/Tc-exd during metamorphosis did not match predictions based on embryonic roles of these genes. Comparing metamorphic patterning in T. castaneum to embryonic and post-embryonic development in hemimetabolous insects suggests that holometabolous metamorphosis combines patterning processes of both late embryogenesis and metamorphosis of the hemimetabolous life cycle.

Metamorphic labral axis patterning in the beetle Tribolium castaneum requires multiple upstream, but few downstream, genes in the appendage patterning network.

The arthropod labrum is an anterior appendage-like structure that forms the dorsal side of the preoral cavity. Conflicting interpretations of fossil, nervous system, and developmental data have led to a proliferation of scenarios for labral evolution. The best supported hypothesis is that the labrum is a novel structure that shares development with appendages as a result of co-option. Here, we use RNA interference in the red flour beetle Tribolium castaneum to compare metamorphic patterning of the labrum to previously published data on ventral appendage patterning. As expected under the co-option hypothesis, depletion of several genes resulted in similar defects in the labrum and ventral appendages. These include proximal deletions and proximal-to-distal transformations resulting from depletion of the leg gap genes homothorax and extradenticle, large-scale deletions resulting from depletion of the leg gap gene Distal-less, and smaller distal deletions resulting from knockdown of the EGF ligand Keren. However, depletion of dachshund and many of the genes that function downstream of the leg gap genes in the ventral appendages had either subtle or no effects on labral axis patterning. This pattern of partial similarity suggests that upstream genes act through different downstream targets in the labrum. We also discovered that many appendage axis patterning genes have roles in patterning the epipharyngeal sensillum array, suggesting that they have become integrated into a novel regulatory network. These genes include Notch, Delta, and decapentaplegic, and the transcription factors abrupt, bric à brac, homothorax, extradenticle and the paralogs apterous a and apterous b.

A functional genetic analysis in flour beetles (Tenebrionidae) reveals an antennal identity specification mechanism active during metamorphosis in Holometabola.

The antenna was the first arthropod ventral appendage to evolve non-leg identity. Models of antennal evolution have been based on comparisons of antennal and leg identity specification mechanisms in Drosophila melanogaster, a species in which appendages develop from highly derived imaginal discs during the larval period. We test for conservation of the Drosophila antennal identity specification mechanism at metamorphosis in Tribolium castaneum and three other flour beetle species (Tribolium confusum, Tribolium brevicornis and Latheticus oryzae) in the family Tenebrionidae. In Drosophila, loss of function of four transcription factors-homothorax, extradenticle, Distal-less, and spineless-causes large-scale transformations of the antenna to leg identity. Distal-less and spineless function similarly during metamorphosis in T. castaneum. RNA interference (RNAi) targeting homothorax (hth) or extradenticle (exd) caused transformation of the proximal antenna to distal leg identity in flour beetles, but did not affect the identity of the distal antenna. This differs from the functional domain of these genes in early instar Drosophila, where they are required for identity specification throughout the antenna, but matches their functional domain in late instar Drosophila. The similarities between antennal identity specification at metamorphosis in flour beetles and in late larval Drosophila likely reflect the conservation of an ancestral metamorphic developmental mechanism. There were two notable differences in hth/exd loss of function phenotypes between flies and beetles. Flour beetles retained all of their primary segments in both the antenna and legs, whereas flies undergo reduction and fusion of primary segments. This difference in ground state appendage morphology casts doubt on interpretations of developmental ground states as evolutionary atavisms. Additionally, adult Tribolium eyes were transformed to elytron-like structures; we provide a developmental hypothesis for this evolutionarily surprising transformation.

Distal-less and dachshund pattern both plesiomorphic and apomorphic structures in chelicerates: RNA interference in the harvestman Phalangium opilio (Opiliones).

The discovery of genetic mechanisms that can transform a morphological structure from a plesiomorphic (=primitive) state to an apomorphic (=derived) one is a cardinal objective of evolutionary developmental biology. However, this objective is often impeded for many lineages of interest by limitations in taxonomic sampling, genomic resources, or functional genetic methods. In order to investigate the evolution of appendage morphology within Chelicerata, the putative sister group of the remaining arthropods, we developed an RNA interference (RNAi) protocol for the harvestman Phalangium opilio. We silenced the leg gap genes Distal-less (Dll) and dachshund (dac) in the harvestman via zygotic injections of double-stranded RNA (dsRNA), and used in situ hybridization to confirm RNAi efficacy. Consistent with the conserved roles of these genes in patterning the proximo-distal axis of arthropod appendages, we observed that embryos injected with Dll dsRNA lacked distal parts of appendages and appendage-like structures, such as the labrum, the chelicerae, the pedipalps, and the walking legs, whereas embryos injected with dac dsRNA lacked the medial podomeres femur and patella in the pedipalps and walking legs. In addition, we detected a role for these genes in patterning structures that do not occur in well-established chelicerate models (spiders and mites). Dll RNAi additionally results in loss of the preoral chamber, which is formed from pedipalpal and leg coxapophyses, and the ocularium, a dorsal outgrowth bearing the eyes. In one case, we observed that an embryo injected with dac dsRNA lacked the proximal segment of the chelicera, a plesiomorphic podomere that expresses dac in wild-type embryos. This may support the hypothesis that loss of the cheliceral dac domain underlies the transition to the two-segmented chelicera of derived arachnids.

Extent With Modification: Leg Patterning in the Beetle Tribolium castaneum and the Evolution of Serial Homologs.

Serial homologs are similar structures that develop at different positions within a body plan. These structures share some, but not all, aspects of developmental patterning, and their evolution is thought to be constrained by shared, pleiotropic gene functions. Here we describe the functions of 17 developmental genes during metamorphic development of the legs in the red flour beetle, Tribolium castaneum. This study provides informative comparisons between appendage development in Drosophila melanogaster and T. castaneum, between embryonic and adult development in T. castaneum, and between the development of serially homologous appendages. The leg gap genes Distal-less and dachshund are conserved in function. Notch signaling, the zinc-finger transcription factors related to odd-skipped, and bric-à-brac have conserved functions in promoting joint development. homothorax knockdown alters the identity of proximal leg segments but does not reduce growth. Lim1 is required for intermediate leg development but not distal tarsus and pretarsus development as in D. melanogaster. Development of the tarsus requires decapentaplegic, rotund, spineless, abrupt, and bric-à-brac and the EGF ligand encoded by Keren. Metathoracic legs of T. castaneum have four tarsomeres, whereas other legs have five. Patterns of gene activity in the tarsus suggest that patterning in the middle of the tarsal region, not the proximal- or distal-most areas, is responsible for this difference in segment number. Through comparisons with other recent studies of T. castaneum appendage development, we test hypotheses for the modularity or interdependence of development during evolution of serial homologs.