Michael S. Brewer

De novo characterization of the gene-rich transcriptomes of two color-polymorphic spiders, Theridion grallator and T. californicum (Araneae: Theridiidae), with special reference to pigment genes

BackgroundA number of spider species within the family Theridiidae exhibit a dramatic abdominal (opisthosomal) color polymorphism. The polymorphism is inherited in a broadly Mendelian fashion and in some species consists of dozens of discrete morphs that are convergent across taxa and populations. Few genomic resources exist for spiders. Here, as a first necessary step towards identifying the genetic basis for this trait we present the near complete transcriptomes of two species: the Hawaiian happy-face spider Theridion grallator and Theridion californicum. We mined the gene complement for pigment-pathway genes and examined differential expression (DE) between morphs that are unpatterned (plain yellow) and patterned (yellow with superimposed patches of red, white or very dark brown).ResultsBy deep sequencing both RNA-seq and normalized cDNA libraries from pooled specimens of each species we were able to assemble a comprehensive gene set for both species that we estimate to be 98-99% complete. It is likely that these species express more than 20,000 protein-coding genes, perhaps 4.5% (ca. 870) of which might be unique to spiders. Mining for pigment-associated Drosophila melanogaster genes indicated the presence of all ommochrome pathway genes and most pteridine pathway genes and DE analyses further indicate a possible role for the pteridine pathway in theridiid color patterning.ConclusionsBased upon our estimates, T. grallator and T. californicum express a large inventory of protein-coding genes. Our comprehensive assembly illustrates the continuing value of sequencing normalized cDNA libraries in addition to RNA-seq in order to generate a reference transcriptome for non-model species. The identification of pteridine-related genes and their possible involvement in color patterning is a novel finding in spiders and one that suggests a biochemical link between guanine deposits and the pigments exhibited by these species.

Millipede Taxonomy after 250 Years: Classification and Taxonomic Practices in a Mega-Diverse yet Understudied Arthropod Group

Background The arthropod class Diplopoda is a mega-diverse group comprising >12,000 described millipede species. The history of taxonomic research within the group is tumultuous and, consequently, has yielded a questionable higher-level classification. Few higher-taxa are defined using synapomorphies, and the practice of single taxon descriptions lacking a revisionary framework has produced many monotypic taxa. Additionally, taxonomic and geographic biases render global species diversity estimations unreliable. We test whether the ordinal taxa of the Diplopoda are consistent with regards to underlying taxonomic diversity, attempt to provide estimates for global species diversity, and examine millipede taxonomic effort at a global geographic scale. Methodology/Principal Findings A taxonomic distinctness metric was employed to assess uniformity of millipede ordinal taxa. We found that ordinal-level taxa are not uniform and are likely overinflated with higher-taxa when compared to related groups. Several methods of estimating global species richness were employed (Bayesian, variation in taxonomic productivity, extrapolation from nearly fully described taxa). Two of the three methods provided estimates ranging from 13,413–16,760 species. Variations in geographic diversity show biases to North America and Europe and a paucity of works on tropical taxa. Conclusions/Significance Before taxa can be used in an extensible way, they must be definable with respect to the diversity they contain and the diagnostic characters used to delineate them. The higher classification for millipedes is shown to be problematic from a number of perspectives. Namely, the ordinal taxa are not uniform in their underlying diversity, and millipedes appear to have a disproportionate number of higher-taxa. Species diversity estimates are unreliable due to inconsistent taxonomic effort at temporal, geographic, and phylogenetic scales. Lack of knowledge concerning many millipede groups compounds these issues. Diplopods are likely not unique in this regard as these issues may persist in many other diverse yet poorly studied groups.

Ordinal-Level Phylogenomics of the Arthropod Class Diplopoda (Millipedes) Based on an Analysis of 221 Nuclear Protein-Coding Loci Generated Using Next-Generation Sequence Analyses

Background The ancient and diverse, yet understudied arthropod class Diplopoda, the millipedes, has a muddled taxonomic history. Despite having a cosmopolitan distribution and a number of unique and interesting characteristics, the group has received relatively little attention; interest in millipede systematics is low compared to taxa of comparable diversity. The existing classification of the group comprises 16 orders. Past attempts to reconstruct millipede phylogenies have suffered from a paucity of characters and included too few taxa to confidently resolve relationships and make formal nomenclatural changes. Herein, we reconstruct an ordinal-level phylogeny for the class Diplopoda using the largest character set ever assembled for the group. Methods Transcriptomic sequences were obtained from exemplar taxa representing much of the diversity of millipede orders using second-generation (i.e., next-generation or high-throughput) sequencing. These data were subject to rigorous orthology selection and phylogenetic dataset optimization and then used to reconstruct phylogenies employing Bayesian inference and maximum likelihood optimality criteria. Ancestral reconstructions of sperm transfer appendage development (gonopods), presence of lateral defense secretion pores (ozopores), and presence of spinnerets were considered. The timings of major millipede lineage divergence points were estimated. Results The resulting phylogeny differed from the existing classifications in a number of fundamental ways. Our phylogeny includes a grouping that has never been described (Juliformia+Merocheta+Stemmiulida), and the ancestral reconstructions suggest caution with respect to using spinnerets as a unifying characteristic for the Nematophora. Our results are shown to have significantly stronger support than previous hypotheses given our data. Our efforts represent the first step toward obtaining a well-supported and robust phylogeny of the Diplopoda that can be used to answer many questions concerning the evolution of this ancient and diverse animal group.

Phylogenetics of the millipede genus Brachycybe Wood, 1864 (Diplopoda: Platydesmida: Andrognathidae): Patterns of deep evolutionary history and recent speciation

The genus Brachycybe Wood is a little known group of millipedes comprising eight nominal species distributed throughout North America, Japan, South Korea, Taiwan, and China. The groups species are relatively morphologically homogenous and have been described primarily on the basis of differences in somatic morphology largely ignoring the often-used characters in millipede taxonomy and systematics - male genitalia (the gonopods). The objectives of this study were to survey male gonopods with the aim of evaluating inter-specific variation, assess existing species boundaries and phylogeny using molecular characters, examine the historical biogeography of the genus, and estimate the timing of lineage divergence using a molecular clock. We surveyed two mitochondrial genes (cytochrome c oxidase I and cytochrome b) and one nuclear protein-coding gene (glutamyl and prolyl-tRNA synthetase) each comprising 952, 746, and 555 aligned base pairs respectively. Phylogenetic inference coupled with an analysis of species delimitation using a generalized mixed Yule coalescent model recovered eight species, two of which were considered to be cryptic. Molecular dating analyses coupled with ancestral range reconstructions indicate that the group is quite ancient (age>50 million years) and its origins are likely traced back geographically to the mountains of Californias Sierra Nevada and Coastal Ranges. The genus Brachycybe appears to have expanded its range at least twice out of present day California into eastern North America and at least once into Asia. This study highlights the need for integrative approaches to describe biodiversity and furthers the evidence for cryptic diversity even in groups where genitalia are generally thought to be rapidly evolving, diagnostic features.

Shifting habitats, morphology, and selective pressures: developmental polyphenism in an adaptive radiation of Hawaiian spiders.

Particularly intriguing examples of adaptive radiation are those in which lineages show parallel or convergent evolution, suggesting utilization of similar genetic or developmental pathways. The current study focuses on an adaptive radiation of Hawaiian “spiny‐leg” spiders in which diversification is associated with repeated convergent evolution leading to similar sets of ecomorphs on each island. However, two species on the oldest islands in the archipelago exhibit variability, occurring as two different ecomorphs. More derived species on the younger islands show much less variability, any one species displaying a single ecomorph. We measured ecomorphological features within individuals over time to determine the nature of the variability. Then, using transcriptomes, we conducted lineage‐based tests for selection under varying models and analyses of gene tree versus species tree incongruencies. Our results provide strong evidence that variability in color in Tetragnatha kauaiensis and T. polychromata is associated with development within individuals (polyphenism). Moreover, a total of 28 loci showed a signature of selection associated with loss of the color‐changing phenotype, and 37 loci showed a signature of selection associated with the colonization of a new environment. The results illustrate how developmental polyphenism might provide an avenue for the repeated evolution of ecomorphs during adaptive radiation.

Phylogeny and Classification of the Trapdoor Spider Genus Myrmekiaphila: An Integrative Approach to Evaluating Taxonomic Hypotheses

Background Revised by Bond and Platnick in 2007, the trapdoor spider genus Myrmekiaphila comprises 11 species. Species delimitation and placement within one of three species groups was based on modifications of the male copulatory device. Because a phylogeny of the group was not available these species groups might not represent monophyletic lineages; species definitions likewise were untested hypotheses. The purpose of this study is to reconstruct the phylogeny of Myrmekiaphila species using molecular data to formally test the delimitation of species and species-groups. We seek to refine a set of established systematic hypotheses by integrating across molecular and morphological data sets. Methods and Findings Phylogenetic analyses comprising Bayesian searches were conducted for a mtDNA matrix composed of contiguous 12S rRNA, tRNA-val, and 16S rRNA genes and a nuclear DNA matrix comprising the glutamyl and prolyl tRNA synthetase gene each consisting of 1348 and 481 bp, respectively. Separate analyses of the mitochondrial and nuclear genome data and a concatenated data set yield M. torreya and M. millerae paraphyletic with respect to M. coreyi and M. howelli and polyphyletic fluviatilis and foliata species groups. Conclusions Despite the perception that molecular data present a solution to a crisis in taxonomy, studies like this demonstrate the efficacy of an approach that considers data from multiple sources. A DNA barcoding approach during the species discovery process would fail to recognize at least two species (M. coreyi and M. howelli) whereas a combined approach more accurately assesses species diversity and illuminates speciation pattern and process. Concomitantly these data also demonstrate that morphological characters likewise fail in their ability to recover monophyletic species groups and result in an unnatural classification. Optimizations of these characters demonstrate a pattern of “Dollo evolution” wherein a complex character evolves only once but is lost multiple times throughout the groups history.

New sequencing technologies, the development of genomics tools, and their applications in evolutionary arachnology

Abstract Molecular genetic tools have been a boon to arachnologists for decades and used to study many unique aspects of arachnid biology including genomics, phylogenetics, population genetics, and biogeography. These tools have evolved over time and now provide myriad methods for exploring evolutionary questions. Early tools, while still useful under the proper circumstances, are giving way to a new generation of DNA sequencing technologies. These new platforms yield impressive amounts of data at a fraction of the cost of traditional techniques. Herein, we discuss the history and future of molecular evolutionary arachnology in terms of available genetic/genomic tools and their potential applications, strengths, weaknesses, and relative costs. Next-generation sequencing (NGS) platforms are varied in their methods and potential uses, making high-throughput sequencing studies focusing on a wide array of questions tractable. To date, relatively few studies have employed NGS technologies using arachnids, but many could benefit from using them. Because no model species exist within the class Arachnida, we have a limited understanding of arachnid genomics. With the ever-advancing nature of sequencing technologies and bioinformatics, arachnologists can relatively easily implement NGS studies to bridge the gaps in our understanding and open avenues for deeper and more powerful experiments. To this end, we discuss examples of applications of NGS technologies focusing on arachnid taxa. Despite the allure of acquiring massive quantities of sequence data, we should recognize the limitations of existing NGS technologies and not forsake pre-NGS methods when these technologies could adequately address our questions.

Arthropod Phylogenetics in Light of Three Novel Millipede (Myriapoda: Diplopoda) Mitochondrial Genomes with Comments on the Appropriateness of Mitochondrial Genome Sequence Data for Inferring Deep Level Relationships

Background Arthropods are the most diverse group of eukaryotic organisms, but their phylogenetic relationships are poorly understood. Herein, we describe three mitochondrial genomes representing orders of millipedes for which complete genomes had not been characterized. Newly sequenced genomes are combined with existing data to characterize the protein coding regions of myriapods and to attempt to reconstruct the evolutionary relationships within the Myriapoda and Arthropoda. Results The newly sequenced genomes are similar to previously characterized millipede sequences in terms of synteny and length. Unique translocations occurred within the newly sequenced taxa, including one half of the Appalachioria falcifera genome, which is inverted with respect to other millipede genomes. Across myriapods, amino acid conservation levels are highly dependent on the gene region. Additionally, individual loci varied in the level of amino acid conservation. Overall, most gene regions showed low levels of conservation at many sites. Attempts to reconstruct the evolutionary relationships suffered from questionable relationships and low support values. Analyses of phylogenetic informativeness show the lack of signal deep in the trees (i.e., genes evolve too quickly). As a result, the myriapod tree resembles previously published results but lacks convincing support, and, within the arthropod tree, well established groups were recovered as polyphyletic. Conclusions The novel genome sequences described herein provide useful genomic information concerning millipede groups that had not been investigated. Taken together with existing sequences, the variety of compositions and evolution of myriapod mitochondrial genomes are shown to be more complex than previously thought. Unfortunately, the use of mitochondrial protein-coding regions in deep arthropod phylogenetics appears problematic, a result consistent with previously published studies. Lack of phylogenetic signal renders the resulting tree topologies as suspect. As such, these data are likely inappropriate for investigating such ancient relationships.

Brain transcriptome of the violet-eared waxbill Uraeginthus granatina and recent evolution in the songbird genome

Songbirds are important models for the study of social behaviour and communication. To complement the recent genome sequencing of the domesticated zebra finch, we sequenced the brain transcriptome of a closely related songbird species, the violet-eared waxbill (Uraeginthus granatina). Both the zebra finch and violet-eared waxbill are members of the family Estrildidae, but differ markedly in their social behaviour. Using Roche 454 RNA sequencing, we generated an assembly and annotation of 11 084 waxbill orthologues of 17 475 zebra finch genes (64%), with an average transcript length of 1555 bp. We also identified 5985 single nucleotide polymorphisms (SNPs) of potential utility for future population genomic studies. Comparing the two species, we found evidence for rapid protein evolution (ω) and low polymorphism of the avian Z sex chromosome, consistent with prior studies of more divergent avian species. An intriguing outlier was putative chromosome 4A, which showed a high density of SNPs and low evolutionary rate relative to other chromosomes. Genome-wide ω was identical in zebra finch and violet-eared waxbill lineages, suggesting a similar demographic history with efficient purifying natural selection. Further comparisons of these and other estrildid finches may provide insights into the evolutionary neurogenomics of social behaviour.

Comparative Transcriptomics of Maturity-Associated Color Change in Hawaiian Spiders

Discontinuous variation within individuals is increasingly recognized as playing a role in diversification and ecological speciation. This study is part of an effort to investigate the molecular genetic underpinnings of adaptive radiation in Hawaiian spiders (genus Tetragnatha). This radiation is found throughout the Hawaiian Islands, showing a common pattern of evolutionary progression from older to younger islands. Moreover, the species are characterized by repeated evolution of similar ecomorphs that can be recognized on the basis of color--Green, Maroon, Large Brown, and Small Brown. However, 2 species (including T. kauaiensis from Kauai) are polyphenic, changing from 1 ecomorph (Green) to another (Maroon) at a specific developmental period. The current study focuses on the age-associated color change in the early stages of the radiation to determine whether this ancestral flexibility in phenotype may have translated into diversification of more derived taxa. We conducted a comparative analysis of transcriptome data (expressed genes) from the Maroon morph of T. kauaiensis and T. perreirai (Oahu), which exhibits a single fixed ecomorph (Maroon). Over 70 million sequence reads were generated using Illumina sequencing of messenger RNA. Using reciprocal best hit BLAST searches, 9027 orthologous genes were identified, of which 32 showed signatures of positive selection between the 2 taxa and may be involved in the loss of the ancestral developmental polyphenism and associated switch to separate monophenic ecomorphs. These results provide critical groundwork that will allow us to advance our understanding of the genomic elements associated with adaptive radiations.: Discontinuous variation within individuals is increasingly recognized as playing a role in diversification and ecological speciation. This study is part of an effort to investigate the molecular genetic underpinnings of adaptive radiation in Hawaiian spiders (genus Tetragnatha). This radiation is found throughout the Hawaiian Islands, showing a common pattern of evolutionary progression from older to younger islands. Moreover, the species are characterized by repeated evolution of similar ecomorphs that can be recognized on the basis of color--Green, Maroon, Large Brown, and Small Brown. However, 2 species (including T. kauaiensis from Kauai) are polyphenic, changing from 1 ecomorph (Green) to another (Maroon) at a specific developmental period. The current study focuses on the age-associated color change in the early stages of the radiation to determine whether this ancestral flexibility in phenotype may have translated into diversification of more derived taxa. We conducted a comparative analysis of transcriptome data (expressed genes) from the Maroon morph of T. kauaiensis and T. perreirai (Oahu), which exhibits a single fixed ecomorph (Maroon). Over 70 million sequence reads were generated using Illumina sequencing of messenger RNA. Using reciprocal best hit BLAST searches, 9027 orthologous genes were identified, of which 32 showed signatures of positive selection between the 2 taxa and may be involved in the loss of the ancestral developmental polyphenism and associated switch to separate monophenic ecomorphs. These results provide critical groundwork that will allow us to advance our understanding of the genomic elements associated with adaptive radiations.