Surprisingly rich repertoire of Wnt genes in the demosponge Halisarca dujardini
- Ilya Borisenko†1,
- Marcin Adamski†2, 3,
- Alexander Ereskovsky1, 4 and
- Maja Adamska2, 3Email authorView ORCID ID profile
© The Author(s). 2016
Received: 6 April 2016
Accepted: 5 June 2016
Published: 10 June 2016
Wnt proteins are secreted signalling molecules found in all animal phyla. In bilaterian animals, including humans, Wnt proteins play key roles in development, maintenance of homeostasis and regeneration. While Wnt gene repertoires and roles are strongly conserved between cnidarians and bilaterians, Wnt genes from basal metazoans (sponges, ctenophores, placozoans) are difficult or impossible to assign to the bilaterian + cnidarian orthologous groups. Moreover, dramatic differences in Wnt numbers among basal metazoan exist, with only three present in the genome of Amphimedon queenslandica, a demosponge, and 21 in the genome of Sycon ciliatum, a calcisponge. To gain insight into the ancestral Wnt repertoire and function, we have chosen to investigate Wnt genes in Halisarca dujardini, a demosponge with relatively well described development and regeneration, and a very distant phylogenetic relationship to Amphimedon.
Here we describe generation of a eukaryotic contamination-free transcriptome of Halisarca dujardini, and analysis of Wnt genes repertoire and expression in this species. We have identified ten Wnt genes, with only one orthologous to Amphimedon Wnt, and six appearing to be a result of a lineage specific expansion. Expression analysis carried out by in situ hybridization of adults and larvae revealed that two Halisarca Wnts are expressed in nested domains in the posterior half of the larvae, and six along the adult body axis, with two specific to the osculum. Strikingly, expression of one of the Wnt genes was elevated in the region undergoing regeneration.
Our results demonstrated that the three Poriferan lineages (Demospongiae, Calcarea and Homoloscleromorpha) are characterized by highly diverse Wnt gene repertoires which do not display higher similarity to each other than they do to the non-sponge (i.e. ctenophore, cnidarian and bilaterian) repertoires. This is in striking contrast to the uniform Wnt repertoires in Cnidarians and Bilaterians, suggesting that the Wnt family composition became “fixed” only in the last common ancestor of Cnidarians and Bilaterians. In contrast, expression of Wnt genes in the apical region of sponge adults and the posterior region of sponge larvae suggests conservation of the Wnt role in axial patterning across the animal kingdom.
Wnt genes encode secreted glycoproteins acting as signalling molecules to direct cell proliferation, migration, differentiation and survival during animal development, maintenance of homeostasis and regeneration [1–6]. While some Wnt pathway components have been identified outside of the animal kingdom, Wnt genes themselves are a conserved metazoan innovation [7, 8].
Representatives of the Wnt family have been identified in all animals studied so far, including so-called “basally branching” or non-bilaterian clades: cnidarians [9, 10], placozoans , ctenophores [12, 13] and sponges [14–18]. Wnt repertoires are surprisingly conserved between cnidarians and bilaterians, with 12 of 13 bilaterian orthologs present in the sea anemone, Nematostella vectensis . This conservation appears to extend to function, as demonstrated by involvement of Wnts in segregation of germ layers during gastrulation, in embryonic and adult axial patterning and in restoration of lost body parts in both cnidarians and bilaterians [19–22]. Conservation of the blastoporal Wnt expression in cnidarians and chordates is particularly striking . In cnidarian polyps such as Hydra, this expression persists in the oral region located in the apical part of the adult body . In line with this, over-activation of Wnt signalling in Hydra results in formation of additional structures with head identity . In chordates, where the blastopore gives rise to the anus, Wnt expression and activity confers posterior identity to developing structures [24, 25]. In line with this, over-activation of the Wnt signalling in vertebrates results in loss of anterior structures, while loss of Wnt function results in posterior truncation [26, 27].
Wnt genes identified in placozoans, ctenophores and sponges are difficult or impossible to assign to the bilaterian + cnidarian orthologous groups [12, 16–18]. Yet, Wnt expression in ctenophores and sponges is consistent with conserved involvement in axial patterning [13, 15, 18]. In particular, Wnt genes have been found to be expressed in the larval posterior pole of two major sponge model species: the demosponge Amphimedon queenslandica  and the calcisponge Sycon ciliatum . In addition, Wnt expression is associated with osculum (the major exhalant opening of adult sponges, located at the apical pole) of Sycon ciliatum . Such expression is consistent with homology of the larval and adult body axes between sponges and cnidarians, supporting homologous relationship between the cnidarian mouth and the sponge osculum [18, 28, 29].
While Wnt expression in adult demosponges has not been reported, pharmacological over-activation of the Wnt pathway in a freshwater species, Ephydatia mulleri, resulted in multiplication of the body axis, as evidenced by formation of multiple oscula . This outcome is strikingly similar to Wnt over-activation experiments in cnidarians, resulting in formation of ectopic head structures . Moreover, experiments involving transplantation of oscula demonstrated their organizer properties, in line with organizer properties of the cnidarian head, and animal blastopores in general [10, 31].
Halisarca dujardini (Chondrillida) is a marine demosponge which is very distantly related to Amphimedon queenslandica (Haplosclerida) [32, 33]. Halisarca embryonic development, metamorphosis and regeneration are well described at morphological level [34–36], but sequence resources have been lacking. Here we report generation of a transcriptome dataset and identification of a surprisingly rich Wnt repertoire (ten genes, in contrast to only three present in the genome of Amphimedon). Two of these genes are expressed in nested domains in the posterior half of the larvae, and six along the adult body axis, with two specific to the osculum. Moreover, Wnt expression is elevated in the region undergoing regeneration, suggesting conservation of the Wnt role in axial patterning and restoration of lost body parts across the animal kingdom.
No permits were required to collect sponge specimens in Norwegian waters. Total RNA and gDNA were isolated from wild-collected adult sponges and several hundred larvae freshly released in laboratory conditions. To avoid eukaryotic contaminations, the larvae were washed in sterile-filtered sea water and visually inspected under dissecting microscope. Nucleic acids were isolated using Allprep Mini kit (Qiagen) following manufacturer’s instructions, and the RNA yield and quality were determined using the NanoDrop spectrophotometer (Thermo Scientific) and the Agilent 2100 BioAnalyzer RNA 6000 Nano chip (Agilent Technologies).
Two RNA-Seq libraries were prepared using Illumina TruSeq RNA Library Prep Kit: one from the wild-collected adult specimen and another one from eukaryotic contaminations-free larvae. An additional gDNA library was prepared from the same larvae using Illumina TruSeq DNA Library Prep Kit. The libraries were paired-end sequenced on Illumina HiSeq 2000 with read length of 100.
The transcriptome was assembled de-novo from the two RNA-Seq libraries. The assembly was performed with Trinity 2.1.1  including reads’ pre-processing with Trimmomatic . We have modified Trinity’s final step called Butterfly to use read pairing information: The fasta sequence files prepared for Butterfly runs were supplemented to include both ends for all the fragments (missing-pair reads were added) and option ‘run_as_paired’ was added to all Butterfly commands. Assembled transcriptome was screened to exclude eukaryotic contaminations by aligning reads from the clean juvenile gDNA library. The alignments were done using bowtie with default parameters. Transcriptome contigs not aligned to any of the clean read were removed from the assembly. Assembly was screened for sequencing vectors using blastn against UniVec database. Transcripts of Wnt ligands were identified by sequence homology using tblastn and Wnt proteins from other organisms and are available in TSA under ids HADA01000001 – HADA01000010.
Wnt protein sequences were aligned with Mafft v7.123 using option L-INS-i. Alignment was then manually trimmed to remove poorly aligned and divergent regions. Phylogenetic tree was built using Mr Bayes 3.1  which we modified to incorporate the LG model (as LG was selected as best fit substitution model by ProTest 3) . Mr Bayes was run with two sets of 4 Markov chains each, till standard deviation of split frequencies dropped below 0.01.
In situ hybridization
In situ hybridization has been carried as described for Sycon ciliatum , except that proteinase treatment was 10 min at 37 °C.
Results and discussion
Ten Wnt genes are present in Halisarca dujardini
We have generated transcriptome dataset for Halisarca dujardini representing genes expressed in adult specimens and free-swimming larvae (see Methods for details). Using a variety of sponge, cnidarian and bilaterian sequences we have BLAST-searched this dataset for Wnt genes and recovered ten complete coding protein sequences (Additional file 1). This stands in contrast with only three Wnt genes present in another demosponge, Amphimedon queenslandica, and also differs from the number of 21 genes identified in Sycon ciliatum (Calcarea, Calcaronea) [16, 18]. For comparison, at least eight Wnt genes are present in Homoscleromorph sponges, e.g. Oscarella sp. .
Thus, all so-far studied sponges, representing three Poriferan lineages (Demospongiae, Calcarea and Homoloscleromorpha) are characterized by highly diverse Wnt gene repertoires which do not display higher similarity to each other than they do to the non-sponge (i.e. ctenophore, cnidarian and bilaterian) repertoires. This is in striking contrast to the uniform Wnt repertoires in Cnidarians and Bilaterians, suggesting that the Wnt family composition became “fixed” only in the last common ancestor of Cnidarians and Bilaterians.
Halisarca Wnts are expressed along the adult and larval axes and during regeneration
In ctenophores, cnidarians and calcareous sponges Wnt genes are expressed along the major (oral-aboral or apical-basal) body axis in sets of nestled domains, suggesting existence of a “Wnt code” possibly conveying positional information [4, 9, 10, 13, 18]. In Oscarella lobularis, a homoscleromorph sponge, two Wnt genes are expressed in a complementary fashion with domains in the ostia (multiple openings in the inhalant canals on the surface of the body) and in exopinacoderm surrounding the ostia of adult specimens . It is important to note here that at least six other Wnt genes are present in Oscarella sp. , expression of which has not been reported so far.
We have attempted cloning and expression analysis of all ten Halisarca Wnts. Of these, six genes were expressed in the adult specimens in four unique patterns:
Detection of WntJ expression revealed positive large cells within the mesohyl of some specimens, which were identified as young oocytes upon sectioning (Fig. 2h, i). While during our collections we have not found any specimens with embryos, we have been able to carry out in situ hybridization on larvae released from adults briefly maintained in laboratory conditions. Two Halisarca Wnt genes revealed robust expression in the larvae: WntK (Fig. 2j) throughout most of the posterior hemisphere, except of the polar cells themselves, and WntJ in a band of equatorial cells (Fig. 2k).
Thus, the identified Wnt expression domains encompass the entire apical-basal axis of the adult Halisarca body, with majority of the genes expressed uniquely or predominantly in the osculum. At the same time, Wnt expression is associated with the posterior region of the larvae. These nested patterns, and the prevalence of apical and posterior expression are consistent with the postulated conservation of Wnt role in axial patterning throughout the metazoans.
In addition to the conserved role in axial patterning, Wnt genes are also known to be involved in wound healing and regeneration in many animal lineages [3, 6, 10, 19–22]. We have recently described cellular processes leading to regeneration of the ectosome in H. dujardini , and we wanted to know whether Wnt genes might be involved in these processes. While majority of the Wnt genes did not display detectable expression changes in the regeneration zone, HduWntK expression was prominent in the exopinacocytes surrounding the wound at 12 h after wounding (Fig. 2l, m). These cells are actively involved in the regeneration, as they temporarily dedifferentiate, phagocyte the debris and contribute to restoration of the ectosome by migration and re-differentiation . Thus, as in other animal lineages, the demosponge Wnt pathway is implicated in the regeneration processes.
Transcriptome sequencing of Halisarca dujardini allowed us the first insight into gene repertoire of a demosponge from a previously unexplored order. We have identified ten Wnt genes, nine of which are without orthologs in any previously reported species. While the diversity of Wnt subfamilies is striking, expression of the identified genes suggests conservation of roles in axial patterning and regeneration. We wonder what mechanisms are responsible – or permissive – for the apparent lack of constraints on Wnt protein sequences in sponges (as well as ctenophores and possibly placozoans) as opposed to cnidarians and bilaterians.
Amq, Amphimedon queenslandica; BLAST, Basic Local Alignment Search Tool; ENA, European Nucleotide Archive; gDNA, genomic deoxyribonucleic acid; Hdu, Halisarca dujardini; LG model, Le-Gascuel substitution model; Ml, Mnemiopsis leidyi; Nv, Nematostella vectensis; Oca, Oscarella carmela; Olo, Oscarella lobularis; RNA, ribonucleic acid; Sci, Sycon ciliatum; SRA, Sequence Read Archive; TSA, Transcriptome Shotgun Assembly.
We thank Sven Leininger for providing adult Halisarca RNA. We thank the “Molecular and cell technologies” research resource center of Saint-Petersburg State University for technical support. The authors are grateful to Dr. A. A. Sukhotin, the director of the White Sea Biological Station of the Zoological Institute of the Russian Academy of Sciences, and The Sven Lovén Centre for Marine Sciences for providing facilities for field research.
We acknowledge funding from the Sars International Centre for Marine Molecular Biology to M. Adamska and grants from the Russian Foundation for Basic Research (n° 13-04-01084 and 16-04-00084), research grant № 22.214.171.1244 awarded by Saint-Petersburg State University and European Community grant ASSEMBLE, agreement no. 227799. I. Borisenko was supported by funds from the Master student research training program at the Sars Centre in 2010 and 2011, and a fellowship of French Government for master students in 2013.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the ENA-SRA and ENA-TSA repositories. Reads were deposited in under study ERP013253, biosamples: SAMEA3672117 for RNA-Seq from larvae, SAMEA3672118 for the RNA-Seq from the adult specimen and SAMEA3672119 for gDNA from larvae. Final transcriptome was deposited in ENA-TSA database under accession numbers HADA01000001 – HADA01138992 http://www.ebi.ac.uk/ena/data/view/HADA01000001-HADA01138992.
I.B. carried out database searches and laboratory work and participated in data analysis. Mar. A. participated in sample collection and carried out sequence assembly and phylogenetic analyses. A.E. participated in study design and data analysis. Maj. A. participated in study design, sample collection and data analysis and drafted the manuscript. All authors contributed to the final version of the manuscript.
The authors declare that they have no competing interests.
Consent to publish
Ethical approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Willert K, Nusse R. Wnt proteins. Cold Spring Harb Perspect Biol. 2012;4:a007864.View ArticlePubMedPubMed CentralGoogle Scholar
- Clevers H, Nusse R. Wnt/β-catenin signaling and disease. Cell. 2012;149:1192–205.View ArticlePubMedGoogle Scholar
- De Robertis EM. Wnt signaling in axial patterning and regeneration: lessons from planaria. Sci Signal. 2010;3(127):e21.View ArticleGoogle Scholar
- Guder C, Philipp I, Lengfeld T, Watanabe H, Hobmayer B, Holstein TW. The Wnt code: cnidarians signal the way. Oncogene. 2006;25(57):7450–60.View ArticlePubMedGoogle Scholar
- Petersen CP, Reddien PW. Wnt Signaling and the Polarity of the Primary Body Axis. Cell. 2009;139:1056–68.View ArticlePubMedGoogle Scholar
- Clevers H, Loh KM, Nusse R. Stem cell signaling. An integral program for tissue renewal and regeneration: Wnt signaling and stem cell control. Science. 2014;346:1248012.View ArticlePubMedGoogle Scholar
- Holstein TW. The evolution of the Wnt pathway. Cold Spring Harb Perspect Biol. 2012;4(7):a007922.View ArticlePubMedPubMed CentralGoogle Scholar
- Adamska M. Developmental signalling and emergence of animal multicellularity. In: Ruiz-Trillo I, Nedelcu AM, editorss. Evolutionary Transitions to Multicellular Life, Advances in Marine Genomics 2, Springer Science + Business Media. ISBN: 978-94-017-9641-5.Google Scholar
- Kusserow A, Pang K, Sturm C, Hrouda M, Lentfer J, Schmidt HA, Technau U, von Haeseler A, Hobmayer B, Martindale MQ, Holstein TW. Unexpected complexity of the Wnt gene family in a sea anemone. Nature. 2005;433:156–60.View ArticlePubMedGoogle Scholar
- Lengfeld T, Watanabe H, Simakov O, Lindgens D, Gee L, Law L, Schmidt HA, Ozbek S, Bode H, Holstein TW. Multiple Wnts are involved in Hydra organizer formation and regeneration. Dev Biol. 2009;330:186–99.View ArticlePubMedGoogle Scholar
- Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, Kuo A, Mitros T, Salamov A, Carpenter ML, Signorovitch AY, Moreno MA, Kamm K, Grimwood J, Schmutz J, Shapiro H, Grigoriev IV, Buss LW, Schierwater B, Dellaporta SL, Rokhsar DS. The Trichoplax genome and the nature of placozoans. Nature. 2008;454(7207):955–60.View ArticlePubMedGoogle Scholar
- Pang K, Ryan JF, Comparative Sequencing Program NISC, Mullikin JC, Baxevanis AD, Martindale MQ. Genomic insights into Wnt signaling in an early diverging metazoan, the ctenophore Mnemiopsis leidyi. Evodevo. 2010;1:10.View ArticlePubMedPubMed CentralGoogle Scholar
- Jager M, Dayraud C, Mialot A, Quéinnec E, le Guyader H, Manuel M. Evidence for involvement of Wnt signalling in body polarities, cell proliferation, and the neuro-sensory system in an adult ctenophore. PLoS One. 2013;8:e84363.View ArticlePubMedPubMed CentralGoogle Scholar
- Nichols SA, Dirks W, Pearse JS, King N. Early evolution of animal cell signaling and adhesion genes. Proc Natl Acad Sci U S A. 2006;103(33):12451–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Adamska M, Degnan SM, Green KM, Adamski M, Craigie A, Larroux C, Degnan BM. Wnt and TGF-beta expression in the sponge Amphimedon queenslandica and the origin of metazoan embryonic patterning. P L o S One. 2007;2:e1031.View ArticleGoogle Scholar
- Adamska M, Larroux C, Adamski M, Green K, Lovas E, Koop D, Richards GS, Zwafink C, Degnan BM. Structure and expression of conserved Wnt pathway components in the demosponge Amphimedon queenslandica. Evol Dev. 2010;12(5):494–518.View ArticlePubMedGoogle Scholar
- Lapébie P, Gazave E, Ereskovsky A, Derelle R, Bézac C, Renard E, Houliston E, Borchiellini C. WNT/beta-catenin signalling and epithelial patterning in the homoscleromorph sponge Oscarella. PLoS One. 2009;4:e5823.View ArticlePubMedPubMed CentralGoogle Scholar
- Leininger S, Adamski M, Bergum B, Guder C, Liu J, Laplante M, Bråte J, Hoffmann F, Fortunato S, Jordal S, Rapp HT and Adamska M. Developmental gene expression provides clues to relationships between sponge and eumetazoan body plans. Nature Comm. 2014;5:3905.View ArticleGoogle Scholar
- Whyte JL, Smith AA, Helms JA. Wnt signaling and injury repair. Cold Spring Harb Perspect Biol. 2012;4(8):a008078.View ArticlePubMedPubMed CentralGoogle Scholar
- Chera S, Ghila L, Dobretz K, Wenger Y, Bauer C, Buzgariu W, Martinou J-C, Galliot B. Apoptotic cells provide an unexpected source of Wnt3 signaling to drive Hydra head regeneration. Dev Cell. 2009;17:279–89.View ArticlePubMedGoogle Scholar
- Trevino M, Stefanik DJ, Rodriguez R, Harmon S, Burton PM. Induction of canonical Wnt signaling by alsterpaullone is sufficient for oral tissue fate during regeneration and embryogenesis in Nematostella vectensis. Dev Dyn. 2011;240:2673–9.View ArticlePubMedPubMed CentralGoogle Scholar
- DuBuc TQ, Nikki Traylor-Knowles N, Martindale MQ. Initiating a regenerative response; cellular and molecular features of wound healing in the cnidarian Nematostella vectensis. BMC Biol. 2014;12:24.View ArticlePubMedPubMed CentralGoogle Scholar
- Broun M, Gee L, Reinhardt B, Bode HR. Formation of the head organizer in hydra involves the canonical Wnt pathway. Development. 2005;132(12):2907–16.View ArticlePubMedGoogle Scholar
- Holland LZ, Holland NN, Schubert M. Developmental expression of AmphiWnt1, an amphioxus gene in the Wnt1/wingless subfamily. Dev Genes Evol. 2000;210:522–4.View ArticlePubMedGoogle Scholar
- Onai T, Lin HC, Schubert M, Koop D, Osborne PW, Alvarez S, Alvarez R, Holland ND, Holland LZ. Retinoic acid and Wnt/beta-catenin have complementary roles in anterior/posterior patterning embryos of the basal chordate amphioxus. Dev Biol. 2009;332:223–33.View ArticlePubMedGoogle Scholar
- Mukhopadhyay M, Shtrom S, Rodriguez-Esteban C, Chen L, Tsukui T, Gomer L, Dorward DW, Glinka A, Grinberg A, Huang SP, Niehrs C, Izpisúa Belmonte JC, Westphal H. Dickkopf1 is required for embryonic head induction and limb morphogenesis in the mouse. Dev Cell. 2001;1:423–34.View ArticlePubMedGoogle Scholar
- Pinson KI, Brennan J, Monkley S, Avery BJ, Skarnes WC. An LDL-receptor-related protein mediates Wnt signalling in mice. Nature. 2000;407(6803):535–8.View ArticlePubMedGoogle Scholar
- Haeckel E. On the organization of sponges and their relationship to the corals. Ann Mag Nat Hist. 1870;5(1–13):107–20.View ArticleGoogle Scholar
- Adamska M. Sponges as the Rosetta Stone of Colonial-to-Multicellular Transition. In: Niklas KJ, Newman SA, editors. Multicellularity. Origins and Evolution. Cambridge, Massachusetts; London, England: The MIT Press; 2016. ISBN 978-0-262-03415-9.Google Scholar
- Windsor PJ, Leys SP. Wnt signaling and induction in the sponge aquiferous system: evidence for an ancient origin of the organizer. Evol Dev. 2010;12:481–590.View ArticleGoogle Scholar
- Gilbert SF. Developmental Biology. 6th edition. Sunderland: Sinauer Associates; 2000. Axis Formation in Amphibians: The Phenomenon of the Organizer. Available from: http://www.ncbi.nlm.nih.gov/books/NBK10101/.
- Ereskovsky AV, Lavrov DV, Boury-Esnault N, Vacelet J. Molecular and morphological description of a new species of Halisarca (Demospongiae: Halisarcida) from Mediterranean Sea with redescription of the type species Halisarca dujardini. Zootaxa. 2011;2768:5–31.Google Scholar
- Morrow C, Cárdenas P. Proposal for a revised classification of the Demospongiae (Porifera). Front Zool. 2015;12:7. doi:10.1186/s12983-015-0099-8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ereskovsky AV, Gonobobleva EL. New data on embryonic development of Halisarca dujardini Johnston, 1842 (Demospongiae: Halisarcida). Zoosystema. 2000;22:355–68.Google Scholar
- Gonobobleva EL, Ereskovsky AV. Metamorphosis of the larva of Halisarca dujardini (Demospongiae, Halisarcida). Bull.inst. Roy. Sci. nat. Belgique. Biol. 2004;74:101–15.Google Scholar
- Borisenko IE, Adamska M, Tokina DB, Ereskovsky AV. Transdifferentiation is a driving force of regeneration in Halisarca dujardini (Demospongiae, Porifera). PeerJ. 2015;3:e1211.View ArticlePubMedPubMed CentralGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Huelsenbeck JP, Ronquist F. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.View ArticlePubMedGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27:1164–5.View ArticlePubMedGoogle Scholar
- Fortunato S, Adamski M, Bergum B, Guder C, Jordal S, Leininger S, Zwafink C, Rapp HT, Adamska M. Genome-wide analysis of the Sox family in the calcareous sponge Sycon ciliatum: multiple genes with unique expression patterns. EvoDevo. 2012;3:14.View ArticlePubMedPubMed CentralGoogle Scholar