Transgressive gene expression and expression plasticity under thermal stress in a stable hybrid zone.

Interspecific hybridization can lead to myriad outcomes, including transgressive phenotypes in which the hybrids are more fit than either parent species. Such hybrids may display important traits in the context of climate change, able to respond to novel environmental conditions not previously experienced by the parent populations. While this has been evaluated in an agricultural context, the role of transgressive hybrids under changing conditions in the wild remains largely unexplored; this is especially true regarding transgressive gene expression. Using the blue mussel species complex (genus Mytilus) as a model system, we investigated the effects of hybridization on temperature induced gene expression plasticity by comparing expression profiles in parental species and their hybrids following a 2-week thermal challenge. Hybrid expression plasticity was most often like one parent or the other (50%). However, a large fraction of genes (26%) showed transgressive expression plasticity (i.e. the change in gene expression was either greater or lesser than that of both parent species), while only 2% were intermediately plastic in hybrids. Despite their close phylogenetic relationship, there was limited overlap in the differentially expressed genes responding to temperature, indicating interspecific differences in the responses to high temperature in which responses from hybrids are distinct from both parent species. We also identified differentially expressed long non-coding RNAs (lncRNAs), which we suggest may contribute to species-specific differences in thermal tolerance. Our findings provide important insight into the impact of hybridization on gene expression under warming. We propose transgressive hybrids may play an important role in population persistence under future warming conditions.

A transcriptome-based phylogeny for Polynoidae (Annelida: Aphroditiformia).

Polynoidae is the most diverse radiation of Aphroditiformia and one of the most successful groups of all Annelida in terms of diversity and habitats colonized. With such an unmatched diversity, phylogenetic investigations have struggled to understand their evolutionary relationships. Previous phylogenetic analyses have slowly increased taxon sampling and employed methodologies, but despite their diversity and biological importance, large genomic sampling is limited. To investigate the internal relationships within Polynoidae, we conducted the first phylogenomic analyses of the group based on 12 transcriptomes collected from species inhabiting a broad array of habitats, including shallow and deep waters, as well as hydrothermal vents, anchialine caves and the midwater. Our phylogenomic analyses of Polynoidae recovered congruent tree topologies representing the clades Polynoinae, Macellicephalinae and Lepidonotopodinae. Members of Polynoinae and Macellicephalinae clustered in well-supported and independent clades. In contrast, Lepidonotopodinae taxa were always recovered nested within Macellicephalinae. Though our sampling only covers a small proportion of the species known for Polynoidae, our results provide a robust phylogenomic framework to build from, emphasizing previously hypothesized relationships between Macellicephalinae and Lepidonotopodinae taxa, while providing new insights on the origin of enigmatic cave and pelagic lineages.

Sexual Signals Persist over Deep Time: Ancient Co-option of Bioluminescence for Courtship Displays in Cypridinid Ostracods.

Although the diversity, beauty, and intricacy of sexually selected courtship displays command the attention of evolutionists, the longevity of these traits in deep time is poorly understood. Population-based theory suggests sexual selection could either lower or raise extinction risk, resulting in high or low persistence of lineages with sexually selected traits. Furthermore, empirical studies that directly estimate the longevity of sexually selected traits are uncommon. Sexually selected signals-including bioluminescent courtship-originated multiple times during evolution, allowing the empirical study of their longevity after careful phylogenetic and divergence time analyses. Here, we estimate the first transcriptome-based molecular phylogeny and divergence times of Cypridinidae. We report extreme longevity of bioluminescent courtship, a trait important in mate choice and probably under sexual selection. Our relaxed-clock estimates of divergence times coupled with stochastic character mapping show luminous courtship evolved only once in Cypridinidae-in a Sub-Tribe, we name Luxorina-at least 151 millions of years ago from cypridinid ancestors that used bioluminescence only in antipredator displays, defining a Tribe we name Luminini. This time-calibrated molecular phylogeny of cypridinids will serve as a foundation for integrative and comparative studies on the biochemistry, molecular evolution, courtship, diversification, and ecology of cypridinid bioluminescence. The persistence of luminous courtship for hundreds of millions of years suggests that sexual selection did not cause a rapid loss of associated traits, and that rates of speciation within the group exceeded extinction risk, which may contribute to the persistence of a diverse clade of signaling species. [Ancestral state reconstruction; Biodiversity; co-option; divergence time estimates; macroevolution; Ostracoda; phylogenomics; sexual selection.].

A Phylogenomic Backbone for Gastropod Molluscs.

Gastropods have survived several mass extinctions during their evolutionary history resulting in extraordinary diversity in morphology, ecology, and developmental modes, which complicate the reconstruction of a robust phylogeny. Currently, gastropods are divided into six subclasses: Caenogastropoda, Heterobranchia, Neomphaliones, Neritimorpha, Patellogastropoda, and Vetigastropoda. Phylogenetic relationships among these taxa historically lack consensus, despite numerous efforts using morphological and molecular information. We generated sequence data for transcriptomes derived from 12 taxa belonging to clades with little or no prior representation in previous studies in order to infer the deeper cladogenetic events within Gastropoda and, for the first time, infer the position of the deep-sea Neomphaliones using a phylogenomic approach. We explored the impact of missing data, homoplasy, and compositional heterogeneity on the inferred phylogenetic hypotheses. We recovered a highly supported backbone for gastropod relationships that is congruent with morphological and mitogenomic evidence, in which Patellogastropoda, true limpets, are the sister lineage to all other gastropods (Orthogastropoda) which are divided into two main clades 1) Vetigastropoda $s.l.$ (including Pleurotomariida $+$ Neomphaliones) and 2) Neritimorpha $+$ (Caenogastropoda $+$ Heterobranchia). As such, our results support the recognition of five subclasses (or infraclasses) in Gastropoda: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda, and Heterobranchia. [Compositional heterogeneity; fast-evolving; long-branch attraction; missing data; Mollusca; phylogenetics; systematic error.].

Mangrove diversity is more than fringe deep.

Mangroves form coastal tropical forests in the intertidal zone and are an important component of shoreline protection. In comparison to other tropical forests, mangrove stands are thought to have relatively low genetic diversity with population genetic structure gradually increasing with distance along a coastline. We conducted genetic analyses of mangrove forests across a range of spatial scales; within a 400 m2 parcel comprising 181 Rhizophora mangle (red mangrove) trees, and across four sites ranging from 6-115 km apart in Honduras. In total, we successfully genotyped 269 R. mangle trees, using a panel of 677 SNPs developed with 2b-RAD methodology. Within the 400 m2 parcel, we found two distinct clusters with high levels of genetic differentiation (FST = 0.355), corresponding to trees primarily located on the seaward fringe and trees growing deeper into the forest. In contrast, there was limited genetic differentiation (FST = 0.027-0.105) across the sites at a larger scale, which had been predominantly sampled along the seaward fringe. Within the 400 m2 parcel, the cluster closest to the seaward fringe exhibited low genetic differentiation (FST = 0.014-0.043) with the other Honduran sites, but the cluster further into the forest was highly differentiated from them (FST = 0.326-0.414). These findings contradict the perception that genetic structure within mangroves forests occurs mainly along a coastline and highlights that there is greater genetic structure at fine spatial scales.

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.

Draft Genome Assemblies and Annotations of Agrypnia vestita Walker, and Hesperophylax magnus Banks Reveal Substantial Repetitive Element Expansion in Tube Case-Making Caddisflies (Insecta: Trichoptera).

Trichoptera (caddisflies) play an essential role in freshwater ecosystems; for instance, larvae process organic material from the water and are food for a variety of predators. Knowledge on the genomic diversity of caddisflies can facilitate comparative and phylogenetic studies thereby allowing scientists to better understand the evolutionary history of caddisflies. Although Trichoptera are the most diverse aquatic insect order, they remain poorly represented in terms of genomic resources. To date, all long-read based genomes have been sequenced from individuals in the retreat-making suborder, Annulipalpia, leaving ∼275 Ma of evolution without high-quality genomic resources. Here, we report the first long-read based de novo genome assemblies of two tube case-making Trichoptera from the suborder Integripalpia, Agrypnia vestita Walker and Hesperophylax magnus Banks. We find that these tube case-making caddisflies have genome sizes that are at least 3-fold larger than those of currently sequenced annulipalpian genomes and that this pattern is at least partly driven by major expansion of repetitive elements. In H. magnus, long interspersed nuclear elements alone exceed the entire genome size of some annulipalpian counterparts suggesting that caddisflies have high potential as a model for understanding genome size evolution in diverse insect lineages.

Evidence of two deeply divergent co-existing mitochondrial genomes in the Tuatara reveals an extremely complex genomic organization.

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

Next-generation museum genomics: Phylogenetic relationships among palpimanoid spiders using sequence capture techniques (Araneae: Palpimanoidea).

Historical museum specimens are invaluable for morphological and taxonomic research, but typically the DNA is degraded making traditional sequencing techniques difficult to impossible for many specimens. Recent advances in Next-Generation Sequencing, specifically target capture, makes use of short fragment sizes typical of degraded DNA, opening up the possibilities for gathering genomic data from museum specimens. This study uses museum specimens and recent target capture sequencing techniques to sequence both Ultra-Conserved Elements (UCE) and exonic regions for lineages that span the modern spiders, Araneomorphae, with a focus on Palpimanoidea. While many previous studies have used target capture techniques on dried museum specimens (for example, skins, pinned insects), this study includes specimens that were collected over the last two decades and stored in 70% ethanol at room temperature. Our findings support the utility of target capture methods for examining deep relationships within Araneomorphae: sequences from both UCE and exonic loci were important for resolving relationships; a monophyletic Palpimanoidea was recovered in many analyses and there was strong support for family and generic-level palpimanoid relationships. Ancestral character state reconstructions reveal that the highly modified carapace observed in mecysmaucheniids and archaeids has evolved independently.

Open access genomic resources for terrestrial arthropods.

Genome sequencing initiatives like the Arthropod i5k project and other biodiversity genomics research rely on access to high quality DNA and/or tissue. Global collection initiatives such as the Smithsonian Global Genome Initiative (GGI) and its partner network, the Global Genome Biodiversity Network (GGBN) aim to provide access to these resources at high-quality standards. Here, we review progress toward providing genomic resources (tissues, DNA, genomes) for terrestrial arthropods, a megadiverse animal group, and compare progress in genome sequencing to all other animals.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution.

BACKGROUND: The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum. RESULTS: We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication. CONCLUSIONS: Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

The complete mitochondrial genome of Cerion uva uva (Gastropoda: Panpulmonata: Stylommatophora: Cerionidae).

We report the complete mitochondrial genome sequence of Cerion uva uva (Linnaeus 1758), the type species of the type genus of the family Cerionidae. The mitogenome is 15,043 bp in length, has a base composition of A (28.3%), T (34.4%), C (17.3%) and G (20.0%), and contains 13 protein-coding genes, 2 ribosomal RNA genes, as well as 22 transfer RNA genes. Gene order is the same as in Cerion incanum (Leidy 1851), but differs from those of all other Panpulmonata. This is the second mitochondrial genome sequenced within the family Cerionidae and will contribute to the assessment of the phylogeography of this family throughout the islands of the tropical western Atlantic.

Cementing mussels to oysters in the pteriomorphian tree: a phylogenomic approach.

Mussels (Mytilida) are a group of bivalves with ancient origins and some of the most important commercial shellfish worldwide. Mytilida consists of approximately 400 species found in various littoral and deep-sea environments, and are part of the higher clade Pteriomorphia, but their exact position within the group has been unstable. The multiple adaptive radiations that occurred within Pteriomorphia have rendered phylogenetic classifications difficult and uncertainty remains regarding the relationships among most families. To address this phylogenetic uncertainty, novel transcriptomic data were generated to include all five orders of Pteriomorphia. Our results, derived from complex analyses of large datasets from 41 transcriptomes and evaluating possible pitfalls affecting phylogenetic reconstruction (matrix occupancy, heterogeneity, evolutionary rates, evolutionary models), consistently recover a well-supported phylogeny of Pteriomorphia, with the only exception of the most complete but smallest data matrix (Matrix 3: 51 genes, 90% gene occupancy). Maximum-likelihood and Bayesian mixture model analyses retrieve strong support for: (i) the monophyly of Pteriomorphia, (ii) Mytilida as a sister group to Ostreida, and (iii) Arcida as sister group to all other pteriomorphians. The basal position of Arcida is congruent with its shell microstructure (solely composed of aragonitic crystals), whereas Mytilida and Ostreida display a combination of a calcitic outer layer with an aragonitic inner layer composed of nacre tablets, the latter being secondarily lost in Ostreoidea.

Phylogenomics and Divergence Dating of Fungus-Farming Ants (Hymenoptera: Formicidae) of the Genera Sericomyrmex and Apterostigma.

Fungus-farming ("attine") ants are model systems for studies of symbiosis, coevolution, and advanced eusociality. A New World clade of nearly 300 species in 15 genera, all attine ants cultivate fungal symbionts for food. In order to better understand the evolution of ant agriculture, we sequenced, assembled, and analyzed transcriptomes of four different attine ant species in two genera: three species in the higher-attine genus Sericomyrmex and a single lower-attine ant species, Apterostigma megacephala, representing the first genomic data for either genus. These data were combined with published genomes of nine other ant species and the honey bee Apis mellifera for phylogenomic and divergence-dating analyses. The resulting phylogeny confirms relationships inferred in previous studies of fungus-farming ants. Divergence-dating analyses recovered slightly older dates than most prior analyses, estimating that attine ants originated 53.6-66.7 million of years ago, and recovered a very long branch subtending a very recent, rapid radiation of the genus Sericomyrmex. This result is further confirmed by a separate analysis of the three Sericomyrmex species, which reveals that 92.71% of orthologs have 99% - 100% pairwise-identical nucleotide sequences. We searched the transcriptomes for genes of interest, most importantly argininosuccinate synthase and argininosuccinate lyase, which are functional in other ants but which are known to have been lost in seven previously studied attine ant species. Loss of the ability to produce the amino acid arginine has been hypothesized to contribute to the obligate dependence of attine ants upon their cultivated fungi, but the point in fungus-farming ant evolution at which these losses occurred has remained unknown. We did not find these genes in any of the sequenced transcriptomes. Although expected for Sericomyrmex species, the absence of arginine anabolic genes in the lower-attine ant Apterostigma megacephala strongly suggests that the loss coincided with the origin of attine ants.

A phylogenetic backbone for Bivalvia: an RNA-seq approach.

Bivalves are an ancient and ubiquitous group of aquatic invertebrates with an estimated 10 000-20 000 living species. They are economically significant as a human food source, and ecologically important given their biomass and effects on communities. Their phylogenetic relationships have been studied for decades, and their unparalleled fossil record extends from the Cambrian to the Recent. Nevertheless, a robustly supported phylogeny of the deepest nodes, needed to fully exploit the bivalves as a model for testing macroevolutionary theories, is lacking. Here, we present the first phylogenomic approach for this important group of molluscs, including novel transcriptomic data for 31 bivalves obtained through an RNA-seq approach, and analyse these data with published genomes and transcriptomes of other bivalves plus outgroups. Our results provide a well-resolved, robust phylogenetic backbone for Bivalvia with all major lineages delineated, addressing long-standing questions about the monophyly of Protobranchia and Heterodonta, and resolving the position of particular groups such as Palaeoheterodonta, Archiheterodonta and Anomalodesmata. This now fully resolved backbone demonstrates that genomic approaches using hundreds of genes are feasible for resolving phylogenetic questions in bivalves and other animals.

Re-evaluating the phylogeny of Sipuncula through transcriptomics.

Sipunculans (also known as peanut worms) are an ancient group of exclusively marine worms with a global distribution and a fossil record that dates back to the Early Cambrian. The systematics of sipunculans, now considered a distinct subclade of Annelida, has been studied for decades using morphological and molecular characters, and has reached the limits of Sanger-based approaches. Here, we reevaluate their family-level phylogeny by comparative transcriptomic analysis of eight species representing all known families within Sipuncula. Two data matrices with alternative gene occupancy levels (large matrix with 675 genes and 62% missing data; reduced matrix with 141 genes and 23% missing data) were analysed using concatenation and gene-tree methods, yielding congruent results and resolving each internal node with maximum support. We thus corroborate prior phylogenetic work based on molecular data, resolve outstanding issues with respect to the familial relationships of Aspidosiphonidae, Antillesomatidae and Phascolosomatidae, and highlight the next area of focus for sipunculan systematics.

Phylogenomic interrogation of arachnida reveals systemic conflicts in phylogenetic signal.

Chelicerata represents one of the oldest groups of arthropods, with a fossil record extending to the Cambrian, and is sister group to the remaining extant arthropods, the mandibulates. Attempts to resolve the internal phylogeny of chelicerates have achieved little consensus, due to marked discord in both morphological and molecular hypotheses of chelicerate phylogeny. The monophyly of Arachnida, the terrestrial chelicerates, is generally accepted, but has garnered little support from molecular data, which have been limited either in breadth of taxonomic sampling or in depth of sequencing. To address the internal phylogeny of this group, we employed a phylogenomic approach, generating transcriptomic data for 17 species in combination with existing data, including two complete genomes. We analyzed multiple data sets containing up to 1,235,912 sites across 3,644 loci, using alternative approaches to optimization of matrix composition. Here, we show that phylogenetic signal for the monophyly of Arachnida is restricted to the 500 slowest-evolving genes in the data set. Accelerated evolutionary rates in Acariformes, Pseudoscorpiones, and Parasitiformes potentially engender long-branch attraction artifacts, yielding nonmonophyly of Arachnida with increasing support upon incrementing the number of concatenated genes. Mutually exclusive hypotheses are supported by locus groups of variable evolutionary rate, revealing significant conflicts in phylogenetic signal. Analyses of gene-tree discordance indicate marked incongruence in relationships among chelicerate orders, whereas derived relationships are demonstrably robust. Consistently recovered and supported relationships include the monophyly of Chelicerata, Euchelicerata, Tetrapulmonata, and all orders represented by multiple terminals. Relationships supported by subsets of slow-evolving genes include Ricinulei + Solifugae; a clade comprised of Ricinulei, Opiliones, and Solifugae; and a clade comprised of Tetrapulmonata, Scorpiones, and Pseudoscorpiones. We demonstrate that outgroup selection without regard for branch length distribution exacerbates long-branch attraction artifacts and does not mitigate gene-tree discordance, regardless of high gene representation for outgroups that are model organisms. Arachnopulmonata (new name) is proposed for the clade comprising Scorpiones + Tetrapulmonata (previously named Pulmonata).

Cocaine-induced vasculitis: is this a new trend?

Cocaine-induced vasculitis is a rare complication found in drug abusers. It occurs due to cocaine adulterated with levamisole. Levamisole was once used as a chemotherapy and immunomodulator for different conditions. One of the side effects of this medication is necrotizing vasculitis which has been reported in the US and Puerto Rico. Here we present another case of cocaine induced vasculitis in Puerto Rico. We describe a 43-year-old female with past medical history of bronchial asthma, migraine, and crack smoking who presented to the emergency room due to blood in her urine for 5 days. She also reported fever, chills, and fatigue. At the physical exam she had a right knee ulcer with swelling erythema, warmth, and pain. Also, she had retiform purpuric plaque lesions in her ears, bilaterally. Eroded plaques with elevated borders at left foot and finger dorsum were also present. Laboratory workup was positive for cocaine. The patient showed leucopenia and microcytic anemia with a normal absolute neutrophil count in her cell blood count. Blood cultures, urine cultures, and ulcer cultures were negative. Urinalysis was positive for proteinuria and hematuria. Also, the patient had positive perinuclear anti-neutrophil cytoplasmic antibody, cytoplasmic anti-neutrophil cytoplasmic antibody, and antinuclear antibody tests and elastase specificity. She showed negative anticardiolipin and lupus anticoagulant antibodies. Her complement levels were decreased. The punch biopsy of her ear showed superficial thrombosis of superficial vascular plexus with perivascular lymphocytic infiltrates and deeper sections showed epidermal necrosis and necrotizing vasculitis. She was started on a high dose of steroids, but could not complete the treatment because she escaped from the hospital before finishing her treatment.

Into the deep: a phylogenetic approach to the bivalve subclass Protobranchia.

A molecular phylogeny of Protobranchia, the subclass of bivalve mollusks sister to the remaining Bivalvia, has long proven elusive, because many constituent lineages are deep-sea endemics, which creates methodological challenges for collecting and preserving genetic material. We obtained 74 representatives of all 12 extant protobranch families and investigated the internal phylogeny of this group using sequence data from five molecular loci (16S rRNA, 18S rRNA, 28S rRNA, cytochrome c oxidase subunit I, and histone H3). Model-based and dynamic homology parsimony approaches to phylogenetic reconstruction unanimously supported four major clades of Protobranchia, irrespective of treatment of hypervariable regions in the nuclear ribosomal genes 18S rRNA and 28S rRNA. These four clades correspond to the superfamilies Nuculoidea (excluding Sareptidae), Nuculanoidea (including Sareptidae), Solemyoidea, and Manzanelloidea. Salient aspects of the phylogeny include (1) support for the placement of the family Sareptidae with Nuculanoidea; (2) the non-monophyly of the order Solemyida (Solemyidae+Nucinellidae); (3) and the non-monophyly of most nuculoid and nuculanoid genera and families. In light of this first family-level phylogeny of Protobranchia, we present a revised classification of the group. Estimation of divergence times in concert with analyses of diversification rates demonstrate the signature of the end-Permian mass extinction in the phylogeny of extant protobranchs.

Comparative description of ten transcriptomes of newly sequenced invertebrates and efficiency estimation of genomic sampling in non-model taxa.

INTRODUCTION: Traditionally, genomic or transcriptomic data have been restricted to a few model or emerging model organisms, and to a handful of species of medical and/or environmental importance. Next-generation sequencing techniques have the capability of yielding massive amounts of gene sequence data for virtually any species at a modest cost. Here we provide a comparative analysis of de novo assembled transcriptomic data for ten non-model species of previously understudied animal taxa. RESULTS: cDNA libraries of ten species belonging to five animal phyla (2 Annelida [including Sipuncula], 2 Arthropoda, 2 Mollusca, 2 Nemertea, and 2 Porifera) were sequenced in different batches with an Illumina Genome Analyzer II (read length 100 or 150 bp), rendering between ca. 25 and 52 million reads per species. Read thinning, trimming, and de novo assembly were performed under different parameters to optimize output. Between 67,423 and 207,559 contigs were obtained across the ten species, post-optimization. Of those, 9,069 to 25,681 contigs retrieved blast hits against the NCBI non-redundant database, and approximately 50% of these were assigned with Gene Ontology terms, covering all major categories, and with similar percentages in all species. Local blasts against our datasets, using selected genes from major signaling pathways and housekeeping genes, revealed high efficiency in gene recovery compared to available genomes of closely related species. Intriguingly, our transcriptomic datasets detected multiple paralogues in all phyla and in nearly all gene pathways, including housekeeping genes that are traditionally used in phylogenetic applications for their purported single-copy nature. CONCLUSIONS: We generated the first study of comparative transcriptomics across multiple animal phyla (comparing two species per phylum in most cases), established the first Illumina-based transcriptomic datasets for sponge, nemertean, and sipunculan species, and generated a tractable catalogue of annotated genes (or gene fragments) and protein families for ten newly sequenced non-model organisms, some of commercial importance (i.e., Octopus vulgaris). These comprehensive sets of genes can be readily used for phylogenetic analysis, gene expression profiling, developmental analysis, and can also be a powerful resource for gene discovery. The characterization of the transcriptomes of such a diverse array of animal species permitted the comparison of sequencing depth, functional annotation, and efficiency of genomic sampling using the same pipelines, which proved to be similar for all considered species. In addition, the datasets revealed their potential as a resource for paralogue detection, a recurrent concern in various aspects of biological inquiry, including phylogenetics, molecular evolution, development, and cellular biochemistry.