26
|
|
27
|
Guang A, Zapata F, Howison M, Lawrence CE, Dunn CW. An Integrated Perspective on Phylogenetic Workflows. Trends Ecol Evol 2016; 31:116-126. [DOI: 10.1016/j.tree.2015.12.007] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2015] [Revised: 12/02/2015] [Accepted: 12/03/2015] [Indexed: 11/29/2022]
|
28
|
Zapata F, Goetz FE, Smith SA, Howison M, Siebert S, Church SH, Sanders SM, Ames CL, McFadden CS, France SC, Daly M, Collins AG, Haddock SHD, Dunn CW, Cartwright P. Phylogenomic Analyses Support Traditional Relationships within Cnidaria. PLoS One 2015; 10:e0139068. [PMID: 26465609 PMCID: PMC4605497 DOI: 10.1371/journal.pone.0139068] [Citation(s) in RCA: 108] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2015] [Accepted: 09/07/2015] [Indexed: 12/04/2022] Open
Abstract
Cnidaria, the sister group to Bilateria, is a highly diverse group of animals in terms of morphology, lifecycles, ecology, and development. How this diversity originated and evolved is not well understood because phylogenetic relationships among major cnidarian lineages are unclear, and recent studies present contrasting phylogenetic hypotheses. Here, we use transcriptome data from 15 newly-sequenced species in combination with 26 publicly available genomes and transcriptomes to assess phylogenetic relationships among major cnidarian lineages. Phylogenetic analyses using different partition schemes and models of molecular evolution, as well as topology tests for alternative phylogenetic relationships, support the monophyly of Medusozoa, Anthozoa, Octocorallia, Hydrozoa, and a clade consisting of Staurozoa, Cubozoa, and Scyphozoa. Support for the monophyly of Hexacorallia is weak due to the equivocal position of Ceriantharia. Taken together, these results further resolve deep cnidarian relationships, largely support traditional phylogenetic views on relationships, and provide a historical framework for studying the evolutionary processes involved in one of the most ancient animal radiations.
Collapse
|
29
|
González VL, Andrade SCS, Bieler R, Collins TM, Dunn CW, Mikkelsen PM, Taylor JD, Giribet G. A phylogenetic backbone for Bivalvia: an RNA-seq approach. Proc Biol Sci 2015; 282:20142332. [PMID: 25589608 PMCID: PMC4308999 DOI: 10.1098/rspb.2014.2332] [Citation(s) in RCA: 85] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022] Open
Abstract
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.
Collapse
|
30
|
Dunn CW, Ryan JF. The evolution of animal genomes. Curr Opin Genet Dev 2015; 35:25-32. [PMID: 26363125 DOI: 10.1016/j.gde.2015.08.006] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2015] [Revised: 08/18/2015] [Accepted: 08/20/2015] [Indexed: 11/18/2022]
Abstract
Genome sequences are now available for hundreds of species sampled across the animal phylogeny, bringing key features of animal genome evolution into sharper focus. The field of animal evolutionary genomics has focused on identifying and classifying the diversity genomic features, reconstructing the history of evolutionary changes in animal genomes, and testing hypotheses about the evolutionary relationships of animals. The grand challenges moving forward are to connect evolutionary changes in genomes with particular evolutionary changes in phenotypes, and to determine which changes are driven by selection. This will require far greater genome sampling both across and within species, extensive phenotype data, a well resolved animal phylogeny, and advances in comparative methods.
Collapse
|
31
|
Siebert S, Pugh PR, Haddock SHD, Dunn CW. Re-evaluation of characters in Apolemiidae (Siphonophora), with description of two new species from Monterey Bay, California. Zootaxa 2015; 3702:201-32. [PMID: 26146720 DOI: 10.11646/zootaxa.3702.3.1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
Siphonophores are polymorphic planktonic marine Cnidarians. The family Apolemiidae is sister to all other species of physonect and calycophoran siphonophores. Although this enigmatic group arguably includes the longest animal species on the planet, their colony-level organization and growth patterns are not well understood. Here we describe two new apolemiid species: Apolemia lanosa sp. nov. and A. rubriversa sp. nov. We provide detailed descriptions of zooid budding and the organization of mature zooids within the siphosome. Our findings reveal that at least two distinct general patterns of siphosomal organization are found in different Apolemia species. In the first pattern, dispersed organization, zooids independently attach directly to the siphosomal stem. In the second pattern, pedunculate organization, only the gastrozooid is attached directly to the stem, and the other zooids of the cormidium branch from its peduncle. This diversity within Apolemia indicates that fundamental aspects of zooid budding and organization are homoplastic within Siphonophora, as both patterns are also found in other siphonophores. The observations presented here greatly clarify the interpretation of diagnostic characters within Apolemiidae, bear on the status of the three previously described species, provide critical detail for understanding the diversity of colony-level organization in siphonophores, and establish a foundation for the description of additional apolemiid species.
Collapse
|
32
|
Haddock SHD, Dunn CW. Fluorescent proteins function as a prey attractant: experimental evidence from the hydromedusa Olindias formosus and other marine organisms. Biol Open 2015; 4:1094-104. [PMID: 26231627 PMCID: PMC4582119 DOI: 10.1242/bio.012138] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
Although proteins in the green fluorescent protein family (GFPs) have been discovered in a wide array of taxa, their ecological functions in these organisms remain unclear. Many hypothesized roles are related to modifying bioluminescence spectra or modulating the light regime for algal symbionts, but these do not explain the presence of GFPs in animals that are non-luminous and non-symbiotic. Other hypothesized functions are unrelated to the visual signals themselves, including stress responses and antioxidant roles, but these cannot explain the localization of fluorescence in particular structures on the animals. Here we tested the hypothesis that fluorescence might serve to attract prey. In laboratory experiments, the predator was the hydromedusa Olindias formosus (previously known as O. formosa), which has fluorescent and pigmented patches on the tips of its tentacles. The prey, juvenile rockfishes in the genus Sebastes, were significantly more attracted (P<1×10−5) to the medusa's tentacles under lighting conditions where fluorescence was excited and tentacle tips were visible above the background. The fish did not respond significantly when treatments did not include fluorescent structures or took place under yellow or white lights, which did not generate fluorescence visible above the ambient light. Furthermore, underwater observations of the behavior of fishes when presented with a brightly illuminated point showed a strong attraction to this visual stimulus. In situ observations also provided evidence for fluorescent lures as supernormal stimuli in several other marine animals, including the siphonophore Rhizophysa eysenhardti. Our results support the idea that fluorescent structures can serve as prey attractants, thus providing a potential function for GFPs and other fluorescent proteins in a diverse range of organisms.
Collapse
|
33
|
Church SH, Ryan JF, Dunn CW. Automation and Evaluation of the SOWH Test with SOWHAT. Syst Biol 2015; 64:1048-58. [PMID: 26231182 PMCID: PMC4604836 DOI: 10.1093/sysbio/syv055] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2014] [Accepted: 07/24/2015] [Indexed: 01/31/2023] Open
Abstract
The Swofford–Olsen–Waddell–Hillis (SOWH) test evaluates statistical support for incongruent phylogenetic topologies. It is commonly applied to determine if the maximum likelihood tree in a phylogenetic analysis is significantly different than an alternative hypothesis. The SOWH test compares the observed difference in log-likelihood between two topologies to a null distribution of differences in log-likelihood generated by parametric resampling. The test is a well-established phylogenetic method for topology testing, but it is sensitive to model misspecification, it is computationally burdensome to perform, and its implementation requires the investigator to make several decisions that each have the potential to affect the outcome of the test. We analyzed the effects of multiple factors using seven data sets to which the SOWH test was previously applied. These factors include a number of sample replicates, likelihood software, the introduction of gaps to simulated data, the use of distinct models of evolution for data simulation and likelihood inference, and a suggested test correction wherein an unresolved “zero-constrained” tree is used to simulate sequence data. To facilitate these analyses and future applications of the SOWH test, we wrote SOWHAT, a program that automates the SOWH test. We find that inadequate bootstrap sampling can change the outcome of the SOWH test. The results also show that using a zero-constrained tree for data simulation can result in a wider null distribution and higher p-values, but does not change the outcome of the SOWH test for most of the data sets tested here. These results will help others implement and evaluate the SOWH test and allow us to provide recommendations for future applications of the SOWH test. SOWHAT is available for download from https://github.com/josephryan/SOWHAT.
Collapse
|
34
|
Zapata F, Wilson NG, Howison M, Andrade SCS, Jörger KM, Schrödl M, Goetz FE, Giribet G, Dunn CW. Phylogenomic analyses of deep gastropod relationships reject Orthogastropoda. Proc Biol Sci 2015; 281:20141739. [PMID: 25232139 DOI: 10.1098/rspb.2014.1739] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Gastropods are a highly diverse clade of molluscs that includes many familiar animals, such as limpets, snails, slugs and sea slugs. It is one of the most abundant groups of animals in the sea and the only molluscan lineage that has successfully colonized land. Yet the relationships among and within its constituent clades have remained in flux for over a century of morphological, anatomical and molecular study. Here, we re-evaluate gastropod phylogenetic relationships by collecting new transcriptome data for 40 species and analysing them in combination with publicly available genomes and transcriptomes. Our datasets include all five main gastropod clades: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia. We use two different methods to assign orthology, subsample each of these matrices into three increasingly dense subsets, and analyse all six of these supermatrices with two different models of molecular evolution. All 12 analyses yield the same unrooted network connecting the five major gastropod lineages. This reduces deep gastropod phylogeny to three alternative rooting hypotheses. These results reject the prevalent hypothesis of gastropod phylogeny, Orthogastropoda. Our dated tree is congruent with a possible end-Permian recovery of some gastropod clades, namely Caenogastropoda and some Heterobranchia subclades.
Collapse
|
35
|
Church SH, Siebert S, Bhattacharyya P, Dunn CW. The histology of Nanomia bijuga (Hydrozoa: Siphonophora). JOURNAL OF EXPERIMENTAL ZOOLOGY PART B-MOLECULAR AND DEVELOPMENTAL EVOLUTION 2015; 324:435-49. [PMID: 26036693 PMCID: PMC5032985 DOI: 10.1002/jez.b.22629] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/08/2014] [Accepted: 04/10/2015] [Indexed: 12/03/2022]
Abstract
The siphonophore Nanomia bijuga is a pelagic hydrozoan (Cnidaria) with complex morphological organization. Each siphonophore is made up of many asexually produced, genetically identical zooids that are functionally specialized and morphologically distinct. These zooids predominantly arise by budding in two growth zones, and are arranged in precise patterns. This study describes the cellular anatomy of several zooid types, the stem, and the gas‐filled float, called the pneumatophore. The distribution of cellular morphologies across zooid types enhances our understanding of zooid function. The unique absorptive cells in the palpon, for example, indicate specialized intracellular digestive processing in this zooid type. Though cnidarians are usually thought of as mono‐epithelial, we characterize at least two cellular populations in this species which are not connected to a basement membrane. This work provides a greater understanding of epithelial diversity within the cnidarians, and will be a foundation for future studies on N. bijuga, including functional assays and gene expression analyses. J. Exp. Zool. (Mol. Dev. Evol.) 324B:435–449, 2015. © 2015 The Authors. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution Published by Wiley Periodicals, Inc.
Collapse
|
36
|
Siebert S, Goetz FE, Church SH, Bhattacharyya P, Zapata F, Haddock SHD, Dunn CW. Stem cells in Nanomia bijuga (Siphonophora), a colonial animal with localized growth zones. EvoDevo 2015; 6:22. [PMID: 26090088 PMCID: PMC4471933 DOI: 10.1186/s13227-015-0018-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2015] [Accepted: 05/11/2015] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Siphonophores (Hydrozoa) have unparalleled colony-level complexity, precision of colony organization, and functional specialization between zooids (i.e., the units that make up colonies). Previous work has shown that, unlike other colonial animals, most growth in siphonophores is restricted to one or two well-defined growth zones that are the sites of both elongation and zooid budding. It remained unknown, however, how this unique colony growth and development is realized at the cellular level. RESULTS To understand the colony-level growth and development of siphonophores at the cellular level, we characterize the distribution of proliferating cells and interstitial stem cells (i-cells) in the siphonophore Nanomia bijuga. Within the colony, we find evidence that i-cells are present at the tip of the horn, the structure within the growth zone that gives rise to new zooids. Co-localized gene expression of vasa-1, pl10, piwi, nanos-1, and nanos-2 suggests that i-cells persist in the youngest zooid buds and that i-cells become progressively restricted to specific regions within the zooids until they are mostly absent from the oldest zooids. The examined genes remain expressed in gametogenic regions. No evidence for i-cells is found in the stem between maturing zooids. Domains of high cell proliferation include regions where the examined genes are expressed, but also include some areas in which the examined genes were not expressed such as the stem within the growth zones. Cell proliferation in regions devoid of vasa-1, pl10, piwi, nanos-1, and nanos-2 expression indicates the presence of mitotically active epithelial cell lineages and, potentially, progenitor cell populations. CONCLUSIONS We provide the first evidence for i-cells in a siphonophore. Our findings suggest maintenance of i-cell populations at the sites of growth zones and that these sites are the main source of i-cells. This restriction of stem cells to particular regions in the colony, in combination with localized budding and spatial patterning during pro-bud subdivision, may play a major role in facilitating the precision of siphonophore growth. Spatially restricted maintenance of i-cells in mature zooids and absence of i-cells along the stem may explain the reduced developmental plasticity in older parts of the colony.
Collapse
|
37
|
Salinas-Saavedra M, Stephenson TQ, Dunn CW, Martindale MQ. Par system components are asymmetrically localized in ectodermal epithelia, but not during early development in the sea anemone Nematostella vectensis. EvoDevo 2015; 6:20. [PMID: 26101582 PMCID: PMC4476184 DOI: 10.1186/s13227-015-0014-6] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2014] [Accepted: 03/12/2015] [Indexed: 11/10/2022] Open
Abstract
Background The evolutionary origins of cell polarity in metazoan embryos are unclear. In most bilaterian animals, embryonic and cell polarity are set up during embryogenesis with the same molecules being utilized to regulate tissue polarity at different life stages. Atypical protein kinase C (aPKC), lethal giant larvae (Lgl), and Partitioning-defective (Par) proteins are conserved components of cellular polarization, and their role in establishing embryonic asymmetry and tissue polarity have been widely studied in model bilaterian groups. However, the deployment and role of these proteins in animals outside Bilateria has not been studied. We address this by characterizing the localization of different components of the Par system during early development of the sea anemone Nematostella vectensis, a member of the clade Cnidaria, the sister group to bilaterian animals. Results Immunostaining using specific N. vectensis antibodies and the overexpression of mRNA-reporter constructs show that components of the N. vectensis Par system (NvPar-1, NvPar-3, NvPar-6, NvaPKC, and NvLgl) distribute throughout the microtubule cytoskeleton of eggs and early embryos without clear polarization along any embryonic axis. However, they become asymmetrically distributed at later stages, when the embryo forms an ectodermal epithelial layer. NvLgl and NvPar-1 localize in the basolateral cortex, and NvaPKC, NvPar-6, and NvPar-3 at the apical zone of the cell in a manner seen in bilaterian animals. Conclusions The cnidarian N. vectensis exhibits clear polarity at all stages of early embryonic development, which appears to be established independent of the Par system reported in many bilaterian embryos. However, in N. vectensis, using multiple immunohistochemical and fluorescently labeled markers in vivo, components of this system are deployed to organize epithelial cell polarity at later stages of development. This suggests that Par system proteins were co-opted to organize early embryonic cell polarity at the base of the Bilateria and that, therefore, different molecular mechanisms operate in early cnidarian embryogenesis. Electronic supplementary material The online version of this article (doi:10.1186/s13227-015-0014-6) contains supplementary material, which is available to authorized users.
Collapse
|
38
|
Helm RR, Tiozzo S, Lilley MKS, Lombard F, Dunn CW. Comparative muscle development of scyphozoan jellyfish with simple and complex life cycles. EvoDevo 2015; 6:11. [PMID: 25932322 PMCID: PMC4415277 DOI: 10.1186/s13227-015-0005-7] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2014] [Accepted: 03/23/2015] [Indexed: 11/10/2022] Open
Abstract
Background Simple life cycles arise from complex life cycles when one or more developmental stages are lost. This raises a fundamental question - how can an intermediate stage, such as a larva, be removed, and development still produce a normal adult? To address this question, we examined the development in several species of pelagiid jellyfish. Most members of Pelagiidae have a complex life cycle with a sessile polyp that gives rise to ephyrae (juvenile medusae); but one species within Pelagiidae, Pelagia noctiluca, spends its whole life in the water column, developing from a larva directly into an ephyra. In many complex life cycles, adult features develop from cell populations that remain quiescent in larvae, and this is known as life cycle compartmentalization and may facilitate the evolution of direct life cycles. A second type of metamorphic processes, known as remodeling, occurs when adult features are formed through modification of already differentiated larval structures. We examined muscle morphology to determine which of these alternatives may be present in Pelagiidae. Results We first examined the structure and development of polyp and ephyra musculature in Chrysaora quinquecirrha, a close relative of P. noctiluca with a complex life cycle. Using phallotoxin staining and confocal microscopy, we verified that polyps have four to six cord muscles that persist in strobilae and discovered that cord muscles is physically separated from ephyra muscle. When cord muscle is removed from ephyra segments, normal ephyra muscle still develops. This suggests that polyp cord muscle is not necessary for ephyra muscle formation. We also found no evidence of polyp-like muscle in P. noctiluca. In both species, we discovered that ephyra muscle arises de novo in a similar manner, regardless of the life cycle. Conclusions The separate origins of polyp and ephyra muscle in C. quinquecirrha and the absence of polyp-like muscle in P. noctiluca suggest that polyp muscle is not remodeled to form ephyra muscle in Pelagiidae. Life cycle stages in Scyphozoa may instead be compartmentalized. Because polyp muscle is not directly remodeled, this may have facilitated the loss of the polyp stage in the evolution of P. noctiluca. Electronic supplementary material The online version of this article (doi:10.1186/s13227-015-0005-7) contains supplementary material, which is available to authorized users.
Collapse
|
39
|
Zapata F, Wilson NG, Howison M, Andrade SCS, Jörger KM, Schrödl M, Goetz FE, Giribet G, Dunn CW. Correction to Phylogenomic analyses of deep gastropod relationships reject Orthogastropoda. Proc Biol Sci 2015. [DOI: 10.1098/rspb.2014.2941] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
|
40
|
Dunn CW, Giribet G, Edgecombe GD, Hejnol A. Animal Phylogeny and Its Evolutionary Implications. ANNUAL REVIEW OF ECOLOGY EVOLUTION AND SYSTEMATICS 2014. [DOI: 10.1146/annurev-ecolsys-120213-091627] [Citation(s) in RCA: 261] [Impact Index Per Article: 26.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
41
|
Ryan JF, Pang K, Schnitzler CE, Nguyen AD, Moreland RT, Simmons DK, Koch BJ, Francis WR, Havlak P, Smith SA, Putnam NH, Haddock SHD, Dunn CW, Wolfsberg TG, Mullikin JC, Martindale MQ, Baxevanis AD. The genome of the ctenophore Mnemiopsis leidyi and its implications for cell type evolution. Science 2013; 342:1242592. [PMID: 24337300 DOI: 10.1126/science.1242592] [Citation(s) in RCA: 451] [Impact Index Per Article: 41.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
An understanding of ctenophore biology is critical for reconstructing events that occurred early in animal evolution. Toward this goal, we have sequenced, assembled, and annotated the genome of the ctenophore Mnemiopsis leidyi. Our phylogenomic analyses of both amino acid positions and gene content suggest that ctenophores rather than sponges are the sister lineage to all other animals. Mnemiopsis lacks many of the genes found in bilaterian mesodermal cell types, suggesting that these cell types evolved independently. The set of neural genes in Mnemiopsis is similar to that of sponges, indicating that sponges may have lost a nervous system. These results present a newly supported view of early animal evolution that accounts for major losses and/or gains of sophisticated cell types, including nerve and muscle cells.
Collapse
|
42
|
Dunn CW, Howison M, Zapata F. Agalma: an automated phylogenomics workflow. BMC Bioinformatics 2013; 14:330. [PMID: 24252138 PMCID: PMC3840672 DOI: 10.1186/1471-2105-14-330] [Citation(s) in RCA: 123] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2013] [Accepted: 11/13/2013] [Indexed: 12/11/2022] Open
Abstract
BACKGROUND In the past decade, transcriptome data have become an important component of many phylogenetic studies. They are a cost-effective source of protein-coding gene sequences, and have helped projects grow from a few genes to hundreds or thousands of genes. Phylogenetic studies now regularly include genes from newly sequenced transcriptomes, as well as publicly available transcriptomes and genomes. Implementing such a phylogenomic study, however, is computationally intensive, requires the coordinated use of many complex software tools, and includes multiple steps for which no published tools exist. Phylogenomic studies have therefore been manual or semiautomated. In addition to taking considerable user time, this makes phylogenomic analyses difficult to reproduce, compare, and extend. In addition, methodological improvements made in the context of one study often cannot be easily applied and evaluated in the context of other studies. RESULTS We present Agalma, an automated tool that constructs matrices for phylogenomic analyses. The user provides raw Illumina transcriptome data, and Agalma produces annotated assemblies, aligned gene sequence matrices, a preliminary phylogeny, and detailed diagnostics that allow the investigator to make extensive assessments of intermediate analysis steps and the final results. Sequences from other sources, such as externally assembled genomes and transcriptomes, can also be incorporated in the analyses. Agalma is built on the BioLite bioinformatics framework, which tracks provenance, profiles processor and memory use, records diagnostics, manages metadata, installs dependencies, logs version numbers and calls to external programs, and enables rich HTML reports for all stages of the analysis. Agalma includes a small test data set and a built-in test analysis of these data. In addition to describing Agalma, we here present a sample analysis of a larger seven-taxon data set. Agalma is available for download at https://bitbucket.org/caseywdunn/agalma. CONCLUSIONS Agalma allows complex phylogenomic analyses to be implemented and described unambiguously as a series of high-level commands. This will enable phylogenomic studies to be readily reproduced, modified, and extended. Agalma also facilitates methods development by providing a complete modular workflow, bundled with test data, that will allow further optimization of each step in the context of a full phylogenomic analysis.
Collapse
|
43
|
Howison M, Zapata F, Dunn CW. Toward a statistically explicit understanding of de novo sequence assembly. Bioinformatics 2013; 29:2959-63. [PMID: 24021385 DOI: 10.1093/bioinformatics/btt525] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION Draft de novo genome assemblies are now available for many organisms. These assemblies are point estimates of the true genome sequences. Each is a specific hypothesis, drawn from among many alternative hypotheses, of the sequence of a genome. Assembly uncertainty, the inability to distinguish between multiple alternative assembly hypotheses, can be due to real variation between copies of the genome in the sample, errors and ambiguities in the sequenced data and assumptions and heuristics of the assemblers. Most assemblers select a single assembly according to ad hoc criteria, and do not yet report and quantify the uncertainty of their outputs. Those assemblers that do report uncertainty take different approaches to describing multiple assembly hypotheses and the support for each. RESULTS Here we review and examine the problem of representing and measuring uncertainty in assemblies. A promising recent development is the implementation of assemblers that are built according to explicit statistical models. Some new assembly methods, for example, estimate and maximize assembly likelihood. These advances, combined with technical advances in the representation of alternative assembly hypotheses, will lead to a more complete and biologically relevant understanding of assembly uncertainty. This will in turn facilitate the interpretation of downstream analyses and tests of specific biological hypotheses.
Collapse
|
44
|
Abstract
Phylogenetic analyses of gene expression have great potential for addressing a wide range of questions. These analyses will, for example, identify genes that have evolutionary shifts in expression that are correlated with evolutionary changes in morphological, physiological, and developmental characters of interest. This will provide entirely new opportunities to identify genes related to particular phenotypes. There are, however, 3 key challenges that must be addressed for such studies to realize their potential. First, data on gene expression must be measured from multiple species, some of which may be field-collected, and parameterized in such a way that they can be compared across species. Second, it will be necessary to develop comparative phylogenetic methods suitable for large multidimensional datasets. In most phylogenetic comparative studies to date, the number n of independent observations (independent contrasts) has been greater than the number p of variables (characters). The behavior of comparative methods for these classic problems is now well understood under a wide variety of conditions. In studies of gene expression, and in studies based on other high-throughput tools, the number n of samples is dwarfed by the number p of variables. The estimated covariance matrices will be singular, complicating their analysis and interpretation, and prone to spurious results. Third, new approaches are needed to investigate the expression of the many genes whose phylogenies are not congruent with species phylogenies due to gene loss, gene duplication, and incomplete lineage sorting. Here we outline general considerations of project design for phylogenetic analyses of gene expression and suggest solutions to these three categories of challenges. These topics are relevant to high-throughput phenotypic data well beyond gene expression.
Collapse
|
45
|
Siebert S, Robinson MD, Tintori SC, Goetz F, Helm RR, Smith SA, Shaner N, Haddock SHD, Dunn CW. Differential gene expression in the siphonophore Nanomia bijuga (Cnidaria) assessed with multiple next-generation sequencing workflows. PLoS One 2011; 6:e22953. [PMID: 21829563 PMCID: PMC3146525 DOI: 10.1371/journal.pone.0022953] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2011] [Accepted: 07/01/2011] [Indexed: 02/02/2023] Open
Abstract
We investigated differential gene expression between functionally specialized feeding polyps and swimming medusae in the siphonophore Nanomia bijuga (Cnidaria) with a hybrid long-read/short-read sequencing strategy. We assembled a set of partial gene reference sequences from long-read data (Roche 454), and generated short-read sequences from replicated tissue samples that were mapped to the references to quantify expression. We collected and compared expression data with three short-read expression workflows that differ in sample preparation, sequencing technology, and mapping tools. These workflows were Illumina mRNA-Seq, which generates sequence reads from random locations along each transcript, and two tag-based approaches, SOLiD SAGE and Helicos DGE, which generate reads from particular tag sites. Differences in expression results across workflows were mostly due to the differential impact of missing data in the partial reference sequences. When all 454-derived gene reference sequences were considered, Illumina mRNA-Seq detected more than twice as many differentially expressed (DE) reference sequences as the tag-based workflows. This discrepancy was largely due to missing tag sites in the partial reference that led to false negatives in the tag-based workflows. When only the subset of reference sequences that unambiguously have tag sites was considered, we found broad congruence across workflows, and they all identified a similar set of DE sequences. Our results are promising in several regards for gene expression studies in non-model organisms. First, we demonstrate that a hybrid long-read/short-read sequencing strategy is an effective way to collect gene expression data when an annotated genome sequence is not available. Second, our replicated sampling indicates that expression profiles are highly consistent across field-collected animals in this case. Third, the impacts of partial reference sequences on the ability to detect DE can be mitigated through workflow choice and deeper reference sequencing.
Collapse
|
46
|
Edgecombe GD, Giribet G, Dunn CW, Hejnol A, Kristensen RM, Neves RC, Rouse GW, Worsaae K, Sørensen MV. Higher-level metazoan relationships: recent progress and remaining questions. ORG DIVERS EVOL 2011. [DOI: 10.1007/s13127-011-0044-4] [Citation(s) in RCA: 206] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
47
|
Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Baguñà J, Bailly X, Jondelius U, Wiens M, Müller WEG, Seaver E, Wheeler WC, Martindale MQ, Giribet G, Dunn CW. Assessing the root of bilaterian animals with scalable phylogenomic methods. Proc Biol Sci 2009; 276:4261-70. [PMID: 19759036 DOI: 10.1098/rspb.2009.0896] [Citation(s) in RCA: 488] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
A clear picture of animal relationships is a prerequisite to understand how the morphological and ecological diversity of animals evolved over time. Among others, the placement of the acoelomorph flatworms, Acoela and Nemertodermatida, has fundamental implications for the origin and evolution of various animal organ systems. Their position, however, has been inconsistent in phylogenetic studies using one or several genes. Furthermore, Acoela has been among the least stable taxa in recent animal phylogenomic analyses, which simultaneously examine many genes from many species, while Nemertodermatida has not been sampled in any phylogenomic study. New sequence data are presented here from organisms targeted for their instability or lack of representation in prior analyses, and are analysed in combination with other publicly available data. We also designed new automated explicit methods for identifying and selecting common genes across different species, and developed highly optimized supercomputing tools to reconstruct relationships from gene sequences. The results of the work corroborate several recently established findings about animal relationships and provide new support for the placement of other groups. These new data and methods strongly uphold previous suggestions that Acoelomorpha is sister clade to all other bilaterian animals, find diminishing evidence for the placement of the enigmatic Xenoturbella within Deuterostomia, and place Cycliophora with Entoprocta and Ectoprocta. The work highlights the implications that these arrangements have for metazoan evolution and permits a clearer picture of ancestral morphologies and life histories in the deep past.
Collapse
|
48
|
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sørensen MV, Haddock SHD, Schmidt-Rhaesa A, Okusu A, Kristensen RM, Wheeler WC, Martindale MQ, Giribet G. Broad phylogenomic sampling improves resolution of the animal tree of life. Nature 2008; 452:745-9. [PMID: 18322464 DOI: 10.1038/nature06614] [Citation(s) in RCA: 1252] [Impact Index Per Article: 78.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2007] [Accepted: 12/20/2007] [Indexed: 11/09/2022]
Abstract
Long-held ideas regarding the evolutionary relationships among animals have recently been upended by sometimes controversial hypotheses based largely on insights from molecular data. These new hypotheses include a clade of moulting animals (Ecdysozoa) and the close relationship of the lophophorates to molluscs and annelids (Lophotrochozoa). Many relationships remain disputed, including those that are required to polarize key features of character evolution, and support for deep nodes is often low. Phylogenomic approaches, which use data from many genes, have shown promise for resolving deep animal relationships, but are hindered by a lack of data from many important groups. Here we report a total of 39.9 Mb of expressed sequence tags from 29 animals belonging to 21 phyla, including 11 phyla previously lacking genomic or expressed-sequence-tag data. Analysed in combination with existing sequences, our data reinforce several previously identified clades that split deeply in the animal tree (including Protostomia, Ecdysozoa and Lophotrochozoa), unambiguously resolve multiple long-standing issues for which there was strong conflicting support in earlier studies with less data (such as velvet worms rather than tardigrades as the sister group of arthropods), and provide molecular support for the monophyly of molluscs, a group long recognized by morphologists. In addition, we find strong support for several new hypotheses. These include a clade that unites annelids (including sipunculans and echiurans) with nemerteans, phoronids and brachiopods, molluscs as sister to that assemblage, and the placement of ctenophores as the earliest diverging extant multicellular animals. A single origin of spiral cleavage (with subsequent losses) is inferred from well-supported nodes. Many relationships between a stable subset of taxa find strong support, and a diminishing number of lineages remain recalcitrant to placement on the tree.
Collapse
|
49
|
Oota H, Dunn CW, Speed WC, Pakstis AJ, Palmatier MA, Kidd JR, Kidd KK. Conservative evolution in duplicated genes of the primate Class I ADH cluster. Gene 2006; 392:64-76. [PMID: 17204375 DOI: 10.1016/j.gene.2006.11.008] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2006] [Revised: 11/11/2006] [Accepted: 11/15/2006] [Indexed: 11/22/2022]
Abstract
Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.
Collapse
|
50
|
Dunn CW, Wagner GP. The evolution of colony-level development in the Siphonophora (Cnidaria:Hydrozoa). Dev Genes Evol 2006; 216:743-54. [PMID: 16983540 DOI: 10.1007/s00427-006-0101-8] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2006] [Accepted: 07/05/2006] [Indexed: 10/24/2022]
Abstract
Evolutionary developmental biology has focused almost exclusively on multicellular organisms, but there are other relevant levels of biological organization that have remained largely neglected. Animal colonies are made up of multiple physiologically integrated and genetically identical units called zooids that are each homologous to solitary, free-living animals. Siphonophores, a group of pelagic hydrozoans (Cnidaria), have the most complex colony-level organization of all animals. Here the colony-level development of five siphonophore species, strategically sampled across the siphonophore phylogeny, is described from specimens collected using deep-sea submersibles and by self-contained underwater breathing apparatus diving. These species include three cystonects, Bathyphysa sibogae, Rhizophysa filiformis, and Rhizophysa eysenhardti, and two "physonects", Agalma elegans and Nanomia bijuga. These data, together with previous findings, are analyzed in a phylogenetic framework to reconstruct key features of the history of colony-level organization and development in the Siphonophora. It is shown that gonodendra and gastrozooids of the examined cystonects arise as independent buds directly on the stem, whereas probud subdivision (the origin of feeding, reproductive, and other zooids from a single bud) is a synapomorphy of the Codonophora. The origin of probud subdivision is associated with the origin of cormidia as integrated units of colony organization, and may have allowed for greater morphological and ecological diversification in the Codonophora relative to the Cystonectae. It is also found that symmetry is labile in siphonophores, with multiple gains and/or losses of directional asymmetry in the group. This descriptive work will enable future mechanistic and molecular studies of colony-level development in the siphonophores.
Collapse
|