Phototactic preference and its genetic basis in the planulae of the colonial Hydrozoan Hydractinia symbiolongicarpus.

Background: Marine organisms with sessile adults commonly possess motile larval stages that make settlement decisions based on integrating environmental sensory cues. Phototaxis, the movement toward or away from light, is a common behavioral characteristic of aquatic and marine metazoan larvae, and of algae, protists, and fungi. In cnidarians, behavioral genomic investigations of motile planulae larvae have been conducted in anthozoans (corals and sea anemones) and scyphozoans (true jellyfish), but such studies are presently lacking in hydrozoans. Here, we examined the behavioral genomics of phototaxis in planulae of the hydrozoan Hydractinia symbiolongicarpus. Results: A behavioral phototaxis study of day 3 planulae indicated preferential phototaxis to green (523 nm) and blue (470 nm) wavelengths of light, but not red (625 nm) wavelengths. A developmental transcriptome study where planula larvae were collected from four developmental time points for RNA-seq revealed that many genes critical to the physiology and development of ciliary photosensory systems are dynamically expressed in planula development and correspond to the expression of phototactic behavior. Microscopical investigations using immunohistochemistry and in situ hybridization demonstrated that several transcripts with predicted function in photoreceptors, including cnidops class opsin, CNG ion channel, and CRX-like transcription factor, localize to ciliated bipolar sensory neurons of the aboral sensory neural plexus, which is associated with the direction of phototaxis and the site of settlement. Conclusions: The phototactic preference displayed by planulae is consistent with the shallow sandy marine habitats they experience in nature. Our genomic investigations add further evidence of similarities between cnidops-mediated photoreceptors of hydrozoans and other cnidarians and ciliary photoreceptors as found in the eyes of humans and other bilaterians, suggesting aspects of their shared evolutionary history.

Authentication of a lophotrochozoan adipokinetic hormone receptor in a Gastropod, Aplysia californica.

Gonadotropin-releasing hormone (GnRH) superfamily comprises multiple families of signaling peptides in both protostomes and deuterostomes. Among this superfamily, vertebrate GnRH stimulates reproduction, but other GnRH superfamily members elicit diverse pleiotropic effects. Within the GnRH superfamily members, adipokinetic hormone (AKH) and its receptor are well described in ecdysozoans but understudied in other lineages. To fill this knowledge gap, we deorphanized a putative receptor for a lophotrochozoan AKH in a gastropod mollusk, Aplysia californica, and named it Aplca-AKHR. Phylogenetic analysis revealed an orthologous relationship of Aplca-AKHR with ecdysozoan AKHRs and other putative lophotrochozoan AKHRs. Aplca-AKHR bound specifically to the previously identified Aplca-AKH with high affinity and activated the inositol phosphate pathway. Aplca-AKHR was expressed widely among central and peripheral tissues, but most prominently in several central ganglia and the heart. The expression of Aplca-AKHR was downregulated by a hyposaline challenge, consistent with a role in volume and fluid regulation previously described for its ligand, Aplca-AKH. In summary, this is the first pairing of a lophotrochozoan AKH with its cognate receptor. Expression data further support diverse central and peripheral roles, including volume and fluid control, of this ligand/receptor pair.

Multisensory integration by polymodal sensory neurons dictates larval settlement in a brainless cnidarian larva.

Multisensory integration (MSI) combines information from more than one sensory modality to elicit behaviours distinct from unisensory behaviours. MSI is best understood in animals with complex brains and specialized centres for parsing different modes of sensory information, but dispersive larvae of sessile marine invertebrates utilize multimodal environmental sensory stimuli to base irreversible settlement decisions on, and most lack complex brains. Here, we examined the sensory determinants of settlement in actinula larvae of the hydrozoan Ectopleura crocea (Cnidaria), which possess a diffuse nerve net. A factorial settlement study revealed that photo-, chemo- and mechanosensory cues each influenced the settlement response in a complex and hierarchical manner that was dependent on specific combinations of cues, an indication of MSI. Additionally, sensory gene expression over development peaked with developmental competence to settle, which in actinulae, requires cnidocyte discharge. Transcriptome analyses also highlighted several deep homological links between cnidarian and bilaterian mechano-, chemo-, and photosensory pathways. Fluorescent in situ hybridization studies of candidate transcripts suggested cellular partitioning of sensory function among the few cell types that comprise the actinula nervous system, where ubiquitous polymodal sensory neurons expressing putative chemo- and photosensitivity interface with mechanoreceptive cnidocytes. We propose a simple multisensory processing circuit, involving polymodal chemo/photosensory neurons and mechanoreceptive cnidocytes, that is sufficient to explain MSI in actinulae settlement. Our study demonstrates that MSI is not exclusive to complex brains, but likely predated and contextualized their evolution.

Epidermal threads reveal the origin of hagfish slime.

When attacked, hagfishes produce a soft, fibrous defensive slime within a fraction of a second by ejecting mucus and threads into seawater. The rapid setup and remarkable expansion of the slime make it a highly effective and unique form of defense. How this biomaterial evolved is unknown, although circumstantial evidence points to the epidermis as the origin of the thread- and mucus-producing cells in the slime glands. Here, we describe large intracellular threads within a putatively homologous cell type from hagfish epidermis. These epidermal threads averaged ~2 mm in length and ~0.5 µm in diameter. The entire hagfish body is covered by a dense layer of epidermal thread cells, with each square millimeter of skin storing a total of ~96 cm threads. Experimentally induced damage to a hagfish's skin caused the release of threads, which together with mucus, formed an adhesive epidermal slime that is more fibrous and less dilute than the defensive slime. Transcriptome analysis further suggests that epidermal threads are ancestral to the slime threads, with duplication and diversification of thread genes occurring in parallel with the evolution of slime glands. Our results support an epidermal origin of hagfish slime, which may have been driven by selection for stronger and more voluminous slime.

Hagfishes are deep-sea animals, and they represent one of the oldest living relatives of animals with backbones. To defend themselves against predators, they produce a remarkable slime that is reinforced with fibers and can clog a predator's gills, thwarting the attack. The slime deploys in less than half a second, exuding from specialized glands on the hagfish's body and expanding up to 10,000 times its ejected volume. The defensive slime is highly dilute, consisting mostly of sea water, with low concentrations of mucus and strong, silk-like threads that are approximately 20 centimeters long. Where and how hagfish slime evolved remains a mystery. Zeng et al. set out to answer where on the hagfish's body the slime glands originated, and how they may have evolved. First, Zeng et al. examined hagfishes and found that cells in the surface layer of their skin (the epidermis) produce threads roughly two millimeters in length that are released when the hagfish's skin is damaged. These threads mix with the mucus that is produced by ruptured skin cells to form a slime that likely adheres to predators' mouths. This slime could be a precursor of the slime produced by the specialized glands. To test this hypothesis, Zeng et al. analyzed which genes are turned on and off both in the hagfishes' skin and in their slime glands. The patterns they found are consistent with the slime glands originating from the epidermis. Based on these results, Zeng et al. propose that ancient hagfishes first evolved the ability to produce slime with anti-predator effects when their skin was damaged in attacks. Over time, hagfishes that could produce and store more slime and eject it actively into a predator's mouth likely had a better chance of surviving. This advantage may have led to the appearance of increasingly specialized glands that could carry out these functions. The findings of Zeng et al. will be of interest to evolutionary biologists, marine biologists, and those studying the ecology of predator-prey interactions. Because of its unique material properties, hagfish slime is also of interest to biophysicists, bioengineers and those engaged in biomimetic research. The origin of hagfish slime glands is an interesting example of how a new trait evolved, and may provide insight into the evolution of other adaptive traits.

Cophylogeny and convergence shape holobiont evolution in sponge-microbe symbioses.

Symbiotic microbial communities of sponges serve critical functions that have shaped the evolution of reef ecosystems since their origins. Symbiont abundance varies tremendously among sponges, with many species classified as either low microbial abundance (LMA) or high microbial abundance (HMA), but the evolutionary dynamics of these symbiotic states remain unknown. This study examines the LMA/HMA dichotomy across an exhaustive sampling of Caribbean sponge biodiversity and predicts that the LMA symbiotic state is the ancestral state among sponges. Conversely, HMA symbioses, consisting of more specialized microorganisms, have evolved multiple times by recruiting similar assemblages, mostly since the rise of scleractinian-dominated reefs. Additionally, HMA symbioses show stronger signals of phylosymbiosis and cophylogeny, consistent with stronger co-evolutionary interaction in these complex holobionts. These results indicate that HMA holobionts are characterized by increased endemism, metabolic dependence and chemical defences. The selective forces driving these patterns may include the concurrent increase in dissolved organic matter in reef ecosystems or the diversification of spongivorous fishes.

Evolution of a remarkable intracellular polymer and extreme cell allometry in hagfishes.

The size of animal cells rarely scales with body size, likely due to biophysical and physiological constraints.1,2 In hagfishes, gland thread cells (GTCs) each produce a silk-like proteinaceous fiber called a slime thread.3,4 The slime threads impart strength to a hagfish's defensive slime and thus are potentially subject to selection on their function outside of the body.5-8 Body size is of fundamental importance in predator-prey interactions, which led us to hypothesize that larger hagfishes produce longer and stronger slime threads than smaller ones.9 Here, by sampling a range of sizes of hagfish from 19 species, we systematically examined the scaling of GTC and slime-thread dimensions with body size within both phylogenetic and ontogenetic contexts. We found that the length of GTCs varied between 40 and 250 µm and scaled positively with body size, exhibiting an allometric exponent greater than those in other animal cells. Correspondingly, larger hagfishes produce longer and thicker slime threads and thus are equipped to defend against larger predators. With diameter and length varying 4-fold (0.7-4 µm and 5-22 cm, respectively) over a body-size range of 10-128 cm, the slime threads characterize the largest intracellular polymers known in biology. Our results suggest selection for stronger defensive slime in larger hagfishes has driven the evolution of extreme size and allometry of GTCs.

Signal, bias, and the role of transcriptome assembly quality in phylogenomic inference.

BACKGROUND: Phylogenomic approaches have great power to reconstruct evolutionary histories, however they rely on multi-step processes in which each stage has the potential to affect the accuracy of the final result. Many studies have empirically tested and established methodology for resolving robust phylogenies, including selecting appropriate evolutionary models, identifying orthologs, or isolating partitions with strong phylogenetic signal. However, few have investigated errors that may be initiated at earlier stages of the analysis. Biases introduced during the generation of the phylogenomic dataset itself could produce downstream effects on analyses of evolutionary history. Transcriptomes are widely used in phylogenomics studies, though there is little understanding of how a poor-quality assembly of these datasets could impact the accuracy of phylogenomic hypotheses. Here we examined how transcriptome assembly quality affects phylogenomic inferences by creating independent datasets from the same input data representing high-quality and low-quality transcriptome assembly outcomes. RESULTS: By studying the performance of phylogenomic datasets derived from alternative high- and low-quality assembly inputs in a controlled experiment, we show that high-quality transcriptomes produce richer phylogenomic datasets with a greater number of unique partitions than low-quality assemblies. High-quality assemblies also give rise to partitions that have lower alignment ambiguity and less compositional bias. In addition, high-quality partitions hold stronger phylogenetic signal than their low-quality transcriptome assembly counterparts in both concatenation- and coalescent-based analyses. CONCLUSIONS: Our findings demonstrate the importance of transcriptome assembly quality in phylogenomic analyses and suggest that a portion of the uncertainty observed in such studies could be alleviated at the assembly stage.

The Genome of the Softshell Clam Mya arenaria and the Evolution of Apoptosis.

Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

The N termini of the inhibitory γ-subunits of phosphodiesterase-6 (PDE6) from rod and cone photoreceptors differentially regulate transducin-mediated PDE6 activation.

Phosphodiesterase-6 (PDE6) plays a central role in both rod and cone phototransduction pathways. In the dark, PDE6 activity is suppressed by its inhibitory γ-subunit (Pγ). Rhodopsin-catalyzed activation of the G protein transducin relieves this inhibition and enhances PDE6 catalysis. We hypothesized that amino acid sequence differences between rod- and cone-specific Pγs underlie transducin's ability to more effectively activate cone-specific PDE6 than rod PDE6. To test this, we analyzed rod and cone Pγ sequences from all major vertebrate and cyclostome lineages and found that rod Pγ loci are far more conserved than cone Pγ sequences and that most of the sequence differences are located in the N-terminal region. Next we reconstituted rod PDE6 catalytic dimer (Pαß) with various rod or cone Pγ variants and analyzed PDE6 activation upon addition of the activated transducin α-subunit (Gtα*-GTPγS). This analysis revealed a rod-specific Pγ motif (amino acids 9-18) that reduces the ability of Gtα*-GTPγS to activate the reconstituted PDE6. In cone Pγ, Asn-13 and Gln-14 significantly enhanced Gtα*-GTPγS activation of cone Pγ truncation variants. Moreover, we observed that the first four amino acids of either rod or cone Pγ contribute to Gtα*-GTPγS-mediated activation of PDE6. We conclude that physiological differences between rod and cone photoreceptor light responsiveness can be partially ascribed to ancient, highly conserved amino acid differences in the N-terminal regions of Pγ isoforms, demonstrating for the first time a functional role for this region of Pγ in the differential activation of rod and cone PDE6 by transducin.

Unraveling inter-species differences in hagfish slime skein deployment.

Hagfishes defend themselves from fish predators by producing defensive slime consisting of mucous and thread components that interact synergistically with seawater to pose a suffocation risk to their attackers. Deployment of the slime occurs in a fraction of a second and involves hydration of mucous vesicles as well as unraveling of the coiled threads to their full length of ∼150 mm. Previous work showed that unraveling of coiled threads (or 'skeins') in Atlantic hagfish requires vigorous mixing with seawater as well as the presence of mucus, whereas skeins from Pacific hagfish tend to unravel spontaneously in seawater. Here, we explored the mechanisms that underlie these different unraveling modes, and focused on the molecules that make up the skein glue, a material that must be disrupted for unraveling to proceed. We found that Atlantic hagfish skeins are also held together with a protein glue, but compared with Pacific hagfish glue, it is less soluble in seawater. Using SDS-PAGE, we identified several soluble proteins and glycoproteins that are liberated from skeins under conditions that drive unraveling in vitro Peptides generated by mass spectrometry of five of these proteins and glycoproteins mapped strongly to 14 sequences assembled from Pacific hagfish slime gland transcriptomes, with all but one of these sequences possessing homologs in the Atlantic hagfish. Two of these sequences encode unusual acidic proteins that we propose are the structural glycoproteins that make up the skein glue. These sequences have no known homologs in other species and are likely to be unique to hagfishes. Although the ecological significance of the two modes of skein unraveling described here are unknown, they may reflect differences in predation pressure, with selection for faster skein unraveling in the Eptatretus lineage leading to the evolution of a glue that is more soluble.

A mosaic of independent innovations involving eyes shut are critical for the evolutionary transition from fused to open rhabdoms.

A fundamental question in evolutionary biology is how developmental processes are modified to produce morphological innovations while abiding by functional constraints. Here we address this question by investigating the cellular mechanism responsible for the transition between fused and open rhabdoms in ommatidia of apposition compound eyes; a critical step required for the development of visual systems based on neural superposition. Utilizing Drosophila and Tribolium as representatives of fused and open rhabdom morphology in holometabolous insects respectively, we identified three changes required for this innovation to occur. First, the expression pattern of the extracellular matrix protein Eyes Shut (EYS) was co-opted and expanded from mechanosensory neurons to photoreceptor cells in taxa with open rhabdoms. Second, EYS homologs obtained a novel extension of the amino terminus leading to the internalization of a cleaved signal sequence. This amino terminus extension does not interfere with cleavage or function in mechanosensory neurons, but it does permit specific targeting of the EYS protein to the apical photoreceptor membrane. Finally, a specific interaction evolved between EYS and a subset of Prominin homologs that is required for the development of open, but not fused, rhabdoms. Together, our findings portray a case study wherein the evolution of a set of molecular novelties has precipitated the origin of an adaptive photoreceptor cell arrangement.

Cloning and characterization of a second lamprey pituitary glycoprotein hormone, thyrostimulin (GpA2/GpB5).

A novel heterodimeric glycoprotein hormone (GpH) comprised of alpha (GpA2) and beta (GpB5) subunits was discovered in 2002 and called thyrostimulin for its ability to activate the TSH receptor in mammals, but its central function in vertebrates has not been firmly established. We report here the cloning and expression of lamprey (l)GpB5, and its ability to heterodimerize with lGpA2 to form a functional l-thyrostimulin. The full-length cDNA of lGpB5 encodes 174 amino acids with ten conserved cysteine residues and one glycosylation site that is conserved with other vertebrate GpB5 sequences. Phylogenetic and synteny analyses support that lGpB5 belongs to the vertebrate GpB5 clade. Heterodimerization of lGpB5 and lGpA2 was shown by nickel pull-down of histidine-tagged recombinant subunits. RNA transcripts of lGpB5 were detected in the pituitary of lampreys during both parasitic and adult life stages. Intraperitoneal injection with lGnRH-III (100 µg/kg) increased pituitary lGpA2, lGpB5, and lGpHß mRNA expression in sexually mature, adult female lampreys. A recombinant l-thyrostimulin produced by expression of a fusion gene in Pichia pastoris activated lamprey GpH receptors I and II as measured by cAMP enzymeimmunoassay. In contrast to jawed vertebrates that have pituitary LH, FSH, and TSH, our data support that lampreys only have two functional pituitary GpHs, lGpH and l-thyrostimulin, which consist of lGpA2 and unique beta subunits. It is hypothesized that lGpH and l-thyrostimulin differentially regulate reproductive and thyroid activities in some unknown way(s) in lampreys.

Comparative Genomics of Pathogenic and Nonpathogenic Beetle-Vectored Fungi in the Genus Geosmithia.

Geosmithia morbida is an emerging fungal pathogen which serves as a model for examining the evolutionary processes behind pathogenicity because it is one of two known pathogens within a genus of mostly saprophytic, beetle-associated, fungi. This pathogen causes thousand cankers disease in black walnut trees and is vectored into the host via the walnut twig beetle. Geosmithia morbida was first detected in western United States and currently threatens the timber industry concentrated in eastern United States. We sequenced the genomes of G. morbida in a previous study and two nonpathogenic Geosmithia species in this work and compared these species to other fungal pathogens and nonpathogens to identify genes under positive selection in G. morbida that may be associated with pathogenicity. Geosmithia morbida possesses one of the smallest genomes among the fungal species observed in this study, and one of the smallest fungal pathogen genomes to date. The enzymatic profile in this pathogen is very similar to its nonpathogenic relatives. Our findings indicate that genome reduction or retention of a smaller genome may be an important adaptative force during the evolution of a specialized lifestyle in fungal species that occupy a specificniche, such as beetle vectored tree pathogens. We also present potential genes under selection in G. morbida that could be important for adaptation to a pathogenic lifestyle.

Recruitment of Rod Photoreceptors from Short-Wavelength-Sensitive Cones during the Evolution of Nocturnal Vision in Mammals.

Vertebrate ancestors had only cone-like photoreceptors. The duplex retina evolved in jawless vertebrates with the advent of highly photosensitive rod-like photoreceptors. Despite cones being the arbiters of high-resolution color vision, rods emerged as the dominant photoreceptor in mammals during a nocturnal phase early in their evolution. We investigated the evolutionary and developmental origins of rods in two divergent vertebrate retinas. In mice, we discovered genetic and epigenetic vestiges of short-wavelength cones in developing rods, and cell-lineage tracing validated the genesis of rods from S cones. Curiously, rods did not derive from S cones in zebrafish. Our study illuminates several questions regarding the evolution of duplex retina and supports the hypothesis that, in mammals, the S-cone lineage was recruited via the Maf-family transcription factor NRL to augment rod photoreceptors. We propose that this developmental mechanism allowed the adaptive exploitation of scotopic niches during the nocturnal bottleneck early in mammalian evolution.

Ancient origins of metazoan gonadotropin-releasing hormone and their receptors revealed by phylogenomic analyses.

The discovery of genes related to gonadotropin-releasing hormones (GnRH) and their receptors from diverse species has driven important advances in comparative endocrinology. However, our view of the evolutionary histories and nomenclature of these gene families has become inconsistent as several different iterations of GnRH and receptor relationships have been proposed. Whole genome sequence data are now available for most of the major lineages of animals, and an exhaustive view of the phylogenies of GnRH and their receptors is now possible. In this paper, we leverage data from publically available whole genome sequences to present a new phylogenomic analysis of GnRH and GnRH receptors and the distant relatives of each across metazoan phylogeny. Our approach utilizes a phylogenomics pipeline that searches data from 36 whole genome sequences and conducts phylogenetic analyses of gene trees. We provide a comprehensive analysis of the major groupings of GnRH peptides, related hormones and their receptors and provide some suggestions for a new nomenclature. Our study provides a framework for understanding the functional diversification of this family of neuromodulatory peptides and their receptors.

Extracting phylogenetic signal and accounting for bias in whole-genome data sets supports the Ctenophora as sister to remaining Metazoa.

BACKGROUND: Understanding the phylogenetic relationships among major lineages of multicellular animals (the Metazoa) is a prerequisite for studying the evolution of complex traits such as nervous systems, muscle tissue, or sensory organs. Transcriptome-based phylogenies have dramatically improved our understanding of metazoan relationships in recent years, although several important questions remain. The branching order near the base of the tree, in particular the placement of the poriferan (sponges, phylum Porifera) and ctenophore (comb jellies, phylum Ctenophora) lineages is one outstanding issue. Recent analyses have suggested that the comb jellies are sister to all remaining metazoan phyla including sponges. This finding is surprising because it suggests that neurons and other complex traits, present in ctenophores and eumetazoans but absent in sponges or placozoans, either evolved twice in Metazoa or were independently, secondarily lost in the lineages leading to sponges and placozoans. RESULTS: To address the question of basal metazoan relationships we assembled a novel dataset comprised of 1080 orthologous loci derived from 36 publicly available genomes representing major lineages of animals. From this large dataset we procured an optimized set of partitions with high phylogenetic signal for resolving metazoan relationships. This optimized data set is amenable to the most appropriate and computationally intensive analyses using site-heterogeneous models of sequence evolution. We also employed several strategies to examine the potential for long-branch attraction to bias our inferences. Our analyses strongly support the Ctenophora as the sister lineage to other Metazoa. We find no support for the traditional view uniting the ctenophores and Cnidaria. Our findings are supported by Bayesian comparisons of topological hypotheses and we find no evidence that they are biased by long-branch attraction. CONCLUSIONS: Our study further clarifies relationships among early branching metazoan lineages. Our phylogeny supports the still-controversial position of ctenophores as sister group to all other metazoans. This study also provides a workflow and computational tools for minimizing systematic bias in genome-based phylogenetic analyses. Future studies of metazoan phylogeny will benefit from ongoing efforts to sequence the genomes of additional invertebrate taxa that will continue to inform our view of the relationships among the major lineages of animals.

A critical appraisal of the use of microRNA data in phylogenetics.

Recent progress in resolving the tree of life continues to expose relationships that resist resolution, which drives the search for novel sources of information to solve these difficult phylogenetic problems. A recent example, the presence and absence of microRNA families, has been vigorously promoted as an ideal source of phylogenetic data and has been applied to several perennial phylogenetic problems. The utility of such data for phylogenetic inference hinges critically both on developing stochastic models that provide a reasonable description of the process that give rise to these data, and also on the careful validation of those models in real inference scenarios. Remarkably, however, the statistical behavior and phylogenetic utility of microRNA data have not yet been rigorously characterized. Here we explore the behavior and performance of microRNA presence/absence data under a variety of evolutionary models and reexamine datasets from several previous studies. We find that highly heterogeneous rates of microRNA gain and loss, pervasive secondary loss, and sampling error collectively render microRNA-based inference of phylogeny difficult. Moreover, our reanalyses fundamentally alter the conclusions for four of the five studies that we reexamined. Our results indicate that the capacity of miRNA data to resolve the tree of life has been overstated, and we urge caution in their application and interpretation.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Pan-metazoan phylogeny of the DMRT gene family: a framework for functional studies.

The family of Doublesex-Mab-3 Related Transcription factors (DMRTs) includes key regulators of sexual differentiation and neurogenesis. To help understand the functional diversification of this gene family, we examined DMRT gene complements from the whole genome sequences and predicted gene models of 32 animal species representing 12 different phyla and from several non-metazoan outgroups. DMRTs are present in all animals except the sponge Amphimedon queenslandica, but are not found in any of the outgroups, indicating that this gene family is specific to animals and has an ancient pre-eumetazoan origin. Our analyses suggest that DMRT genes diversified independently in bilaterian and non-bilaterian animals. Most clades in the DMRT gene tree, including those containing the well-characterized DMRT1 and doublesex genes, have phylogenetically limited distributions.

Gene co-expression modules underlying polymorphic and monomorphic zooids in the colonial hydrozoan, Hydractinia symbiolongicarpus.

Advances in sequencing technology have forced a quantitative revolution in Evolutionary Biology. One important feature of this renaissance is that comprehensive genomic resources can be obtained quickly for almost any taxon, thus speeding the development of new model organisms. Here, we analyze 20 RNA-seq libraries from morphologically, sexually, and genetically distinct polyp types from the gonochoristic colonial hydrozoan, Hydractinia symbiolongicarpus (Cnidaria). Analyses of these data using weighted gene co-expression networks highlight deeply conserved genetic elements of animal spermatogenesis and demonstrate the utility of these methods in identifying modules of genes that correlate with different zooid types across various statistical contrasts. RNA-seq data and analytical scripts described here are deposited in publicly available databases.