Inactivation of uncoupling protein 1 is linked to shifts in metabolic rate, body size, and species richness of eight mammalian lineages.

Arachnida is an ancient, diverse and ecologically important animal group that contains a number of species of interest for medical, agricultural and engineering applications. Despite their importance, many aspects of the arachnid tree of life remain unresolved, hindering comparative approaches to arachnid biology. Biologists have made considerable efforts to resolve the arachnid phylogeny; yet, limited and challenging morphological characters, as well as a dearth of genetic resources, have hindered progress. Here, we present a genomic toolkit for arachnids featuring hundreds of conserved DNA regions (ultraconserved elements or UCEs) that allow targeted sequencing of any species in the arachnid tree of life. We used recently developed capture probes designed from conserved regions of available arachnid genomes to enrich a sample of loci from 32 diverse arachnids. Sequence capture returned an average of 487 UCE loci for all species, with a range from 170 to 722. Phylogenetic analysis of these UCEs produced a highly resolved arachnid tree with relationships largely consistent with recent transcriptome-based phylogenies. We also tested the phylogenetic informativeness of UCE probes within the spider, scorpion and harvestman orders, demonstrating the utility of these markers at shallower taxonomic scales and suggesting that these loci will be useful for species-level differences. This probe set will open the door to phylogenomic and population genomic studies across the arachnid tree of life, enabling systematics, species delimitation, species discovery and conservation of these diverse arthropods.

Silk spinning is essential to spider ecology and has had a key role in the expansive diversification of spiders. Silk is composed primarily of proteins called spidroins, which are encoded by a multi-gene family. Spidroins have been studied extensively in the derived clade, Orbiculariae (orb-weavers), from the suborder Araneomorphae (‘true spiders’). Orbicularians produce a suite of different silks, and underlying this repertoire is a history of duplication and spidroin gene divergence. A second class of silk proteins, Egg Case Proteins (ECPs), is known only from the orbicularian species, Lactrodectus hesperus (Western black widow). In L. hesperus, ECPs bond with tubuliform spidroins to form egg case silk fibers. Because most of the phylogenetic diversity of spiders has not been sampled for their silk genes, there is limited understanding of spidroin gene family history and the prevalence of ECPs. Silk genes have not been reported from the suborder Mesothelae (segmented spiders), which diverged from all other spiders >380 million years ago, and sampling from Mygalomorphae (tarantulas, trapdoor spiders) and basal araneomorph lineages is sparse. In comparison to orbicularians, mesotheles and mygalomorphs have a simpler silk biology and thus are hypothesized to have less diversity of silk genes. Here, we present cDNAs synthesized from the silk glands of six mygalomorph species, a mesothele, and a non-orbicularian araneomorph, and uncover a surprisingly rich silk gene diversity. In particular, we find ECP homologs in the mesothele, suggesting that ECPs were present in the common ancestor of extant spiders, and originally were not specialized to complex with tubuliform spidroins. Furthermore, gene-tree/species-tree reconciliation analysis reveals that numerous spidroin gene duplications occurred after the split between Mesothelae and Opisthothelae (Mygalomorphae plus Araneomorphae). We use the spidroin gene tree to reconstruct the evolution of amino acid compositions of spidroins that perform different ecological functions.

Antrodiaetus riversi is a dispersal-limited, habitat-specialized mygalomorph spider species endemic to mesic woodlands of northern and central California. Here, we build upon prior phylogeographic research using a much larger geographic sample and include additional nuclear genes, providing more detailed biogeographic insights throughout the range of this complex. Of particular interest is the uncovering of unexpected and replicated trans-valley biogeographic patterns, where in two separate genetic clades western haplotypes in the California south Coast Ranges are phylogenetically closely related to eastern haplotypes from central and northern Sierran foothills. In both instances, these trans-valley phylogenetic patterns are strongly supported by multiple genes. These western and eastern populations are currently separated by the Central Valley, a well-recognized modern-day and historical biogeographic barrier in California. For one clade, the directionality is clearly northeast to southwest, and all available evidence is consistent with a jump dispersal event estimated at 1.2-1.3 Ma. During this time period, paleogeographic data indicate that northern Sierran rivers emptied to the ocean in the south Coast Ranges, rather than at the San Francisco Bay. For the other trans-valley clade genetic evidence is less conclusive regarding the mechanism and directionality of biogeographic exchange, although the estimated timeframe is similar (approximately 1.8 Ma). Despite the large number of biogeographic studies previously conducted in central California, to the best of our knowledge no prior studies have discussed or revealed a northern Sierran to south Coast Range biogeographic connection. This uniqueness may reflect the low-dispersal biology of mygalomorph spiders, where 'post-event' gene exchange rarely erases historical biogeographic signal.

Antrodiaetus riversi (Araneae, Antrodiaetidae) is a dispersal-limited, habitat specialized mygalomorph spider species endemic to mesic woodlands of northern and central California. This species occupies a disjunct distribution, with populations in the Sierra Nevada and Coast Ranges, separated by the inhospitable Central Valley. Previous studies of morphological and allozyme variation have suggested that these populations may constitute cryptic species. We investigated the phylogeography of A. riversi using both nuclear and mitochondrial DNA sequences, collected for a comprehensive population sample. These data reveal the presence of at least five species in the A. riversi complex - these species are deeply diverged, and genealogically exclusive in both nuclear and mitochondrial genomes. Each of these species is characterized by extreme population subdivision and deep phylogeographical structuring, consistent with minimal gene flow across the dissected Californian landscape. Three species are restricted to the Coast Ranges, one to high altitudes of the central Sierran Nevada, and one species is found in both ranges. These species have allopatric distributions, although species parapatry is hypothesized to occur in several areas. Species diversification appears to have pulsed in the Late Miocene/Early Pliocene, a timing consistent with biogeographical reconstructions for many Californian taxa, and a time of turbulent geological activity in the region.

Next-generation sequencing technologies are rapidly transforming molecular systematic studies of non-model animal taxa. The arachnid order Opiliones (commonly known as “harvestmen”) includes more than 6,400 described species placed into four well-supported lineages (suborders). Fossil plus molecular clock evidence indicates that these lineages were diverging in the late Silurian to mid-Carboniferous, with some fossil harvestmen representing the earliest known land animals. Perhaps because of this ancient divergence, phylogenetic resolution of subordinal interrelationships within Opiliones has been difficult. We present the first phylogenomics analysis for harvestmen, derived from comparative RNA-Seq data for eight species representing all suborders. Over 30 gigabases of original Illumina short-read data were used in de novo assemblies, resulting in 50–80,000 transcripts per taxon. Transcripts were compared to published scorpion and tick genomics data, and a stringent filtering process was used to identify over 350 putatively single-copy, orthologous protein-coding genes shared among taxa. Phylogenetic analyses using various partitioning strategies, data coding schemes, and analytical methods overwhelmingly support the “classical” hypothesis of Opiliones relationships, including the higher-level clades Palpatores and Phalangida. Relaxed molecular clock analyses using multiple alternative fossil calibration strategies corroborate ancient divergences within Opiliones that are possibly deeper than the recorded fossil record indicates. The assembled data matrices, comprising genes that are conserved, highly expressed, and varying in length and phylogenetic informativeness, represent an important resource for future molecular systematic studies of Opiliones and other arachnid groups.

31Arachnida is an ancient, diverse, and ecologically important animal group that contains a number 32 of species of interest for medical, agricultural, and engineering applications. Despite this applied 33 importance, many aspects of the arachnid tree of life remain unresolved, hindering comparative 34 approaches to arachnid biology. Biologists have made considerable efforts to resolve the 35 arachnid phylogeny; yet, limited and challenging morphological characters, as well as a dearth of 36 genetic resources, have confounded these attempts. Here, we present a genomic toolkit for 37 arachnids featuring hundreds of conserved DNA regions (ultraconserved elements or UCEs) that 38 allow targeted sequencing of any species in the arachnid tree of life. We used recently developed 39 capture probes designed from conserved genomic regions of available arachnid genomes to 40 enrich a sample of loci from 32 diverse arachnids. Sequence capture returned an average of 487 41 UCE loci for all species, with a range from 170 to 722. Phylogenetic analysis of these UCEs 42 produced a highly resolved arachnid tree with relationships largely consistent with recent 43 transcriptome-based phylogenies. We also tested the phylogenetic informativeness of UCE 44 probes within the spider, scorpion, and harvestman orders, demonstrating the utility of these 45 markers at shallower taxonomic scales, even down to the level of species differences. This probe 46 set will open the door to phylogenomic and population genomic studies across the arachnid tree 47 of life, enabling systematics, species delimitation, species discovery, and conservation of these 48 diverse arthropods. 49 50

Previous studies have reported inactivated copies of six enamel-related genes (AMBN, AMEL, AMTN, ENAM, KLK4, MMP20) and one dentin-related gene (DSPP) in one or more toothless vertebrates and/or vertebrates with enamelless teeth, thereby providing evidence that these genes are enamel or tooth-specific with respect to their critical functions that are maintained by natural selection. Here, we employ available genome sequences for edentulous and enamelless mammals to evaluate the enamel specificity of four genes (WDR72, SLC24A4, FAM83H, C4orf26) that have been implicated in amelogenesis imperfecta, a condition in which proper enamel formation is abrogated during tooth development. Coding sequences for WDR72, SCL24A4, and FAM83H are intact in four edentulous taxa (Chinese pangolin, three baleen whales) and three taxa (aardvark, nine-banded armadillo, Hoffmann's two-toed sloth) with enamelless teeth, suggesting that these genes have critical functions beyond their involvement in tooth development. By contrast, genomic data for C4orf26 reveal inactivating mutations in pangolin and bowhead whale as well as evidence for deletion of this gene in two minke whale species. Hybridization capture of exonic regions and PCR screens provide evidence for inactivation of C4orf26 in eight additional baleen whale species. However, C4orf26 is intact in all three species with enamelless teeth that were surveyed, as well as in 95 additional mammalian species with enamel-capped teeth. Estimates of selection intensity suggest that dN/dS ratios on branches leading to taxa with enamelless teeth are similar to the dN/dS ratio on branches leading to taxa with enamel-capped teeth. Based on these results, we conclude that C4orf26 is tooth-specific, but not enamel-specific, with respect to its essential functions that are maintained by natural selection. A caveat is that an alternative splice site variant, which translates exon 3 in a different reading frame, is putatively functional in Catarrhini and may have evolved an additional role in this primate clade.