Developmental biology and potential use of Alboglossiphonia lata (Annelida: Hirudinea) as an “Evo-Devo” model organism
- Brenda Irene Medina Jiménez†1,
- Hee-Jin Kwak†1,
- Jong-Seok Park1,
- Jung-Woong Kim2Email author and
- Sung-Jin Cho1Email author
© The Author(s). 2017
Received: 5 September 2017
Accepted: 27 October 2017
Published: 28 December 2017
The need for the adaptation of species of annelids as “Evo-Devo” model organisms of the superphylum Lophotrochozoa to refine the understanding of the phylogenetic relationships between bilaterian organisms, has promoted an increase in the studies dealing with embryonic development among related species such as leeches from the Glossiphoniidae family. The present study aims to describe the embryogenesis of Alboglossiphonia lata (Oka, 1910), a freshwater glossiphoniid leech, chiefly distributed in East Asia, and validate standard molecular biology techniques to support the use of this species as an additional model for “Evo-Devo” studies.
A. lata undergoes direct development, and follows the highly conserved clitellate annelid mode of spiral cleavage development; the duration from the egg laying to the juvenile stage is ~7.5 days, and it is iteroparous, indicating that it feeds and deposits eggs again after the first round of brooding, as described in several other glossiphoniid leech species studied to date. The embryos hatch only after complete organ development and proboscis retraction, which has not yet been observed in other glossiphoniid genera. The phylogenetic position of A. lata within the Glossiphoniidae family has been confirmed using cytochrome c oxidase subunit 1 (CO1) sequencing. Lineage tracer injections confirmed the fates of the presumptive meso- and ectodermal precursors, and immunostaining showed the formation of the ventral nerve system during later stages of development. Further, the spatiotemporal expression of an EF-hand calcium-binding protein Calsensin ortholog was characterized, which showed a specific pattern in both the ventral and peripheral nervous systems during the later stages.
Our description of the embryonic development of A. lata under laboratory conditions provides new data for further comparative studies with other leech and lophotrochozoa model organisms. Moreover, it offers a basis for the establishment of this species as a model for future “Evo-Devo” studies.
The study of new non-model organisms such as annelids has gained more attention in recent years. Of the three bilaterian clades, namely, Deuterostomia, Ecdysozoa, and Lophotrochozoa [1, 2], the latter remains the least represented clade because of the preference for classic, genetic model organisms. This has led to gaps in understanding the evolutionary history of the bilaterians . Developmental studies reveal a mixture of conserved and derived features, which are interpreted in light of the underlying phylogenetic relationships that are established independently by molecular phylogenetic analysis. Understanding of the actual mechanisms that shape development and evolution requires detailed knowledge of the cellular processes occurring during embryogenesis, a more highly resolved phylogenetic tree for annelids and their allies, and the inclusion of more species in comparative studies . Earthworms and leeches have been studied to address this concern. This allows the establishment of new model organisms for the members of Lophotrochozoa such as Helobdella austinensis (Kutschera et al. 2013) , which provides a reference for studying satellite species .
Alboglossiphonia lata (Oka, 1910)  belongs to the Glossiphoniidae family, which is among the more species-rich leech families in terms of described numbers of species . Worldwide, the presence of A. lata is primarily recorded in East Asia, including China, Japan, Taiwan, and South Korea, as well as in Hawaii [9, 10]. Glossiphoniidae leeches are characterized for having dorso-ventrally flattened and dorsally convex bodies, and for bearing a proboscis. These leeches usually feed on the blood of turtle or amphibians in clean, non-organic polluted streams, irrigation ditches, and open sewers. However, some glossiphoniids, like those belonging to the genera Helobdella (Blanchard, 1896) and Glossiphonia (Johnson, 1816), feed on the hemolymph of aquatic oligochaetes and snails . Regarding parental care, all known Glossiphoniidae have evolved the habit of brooding the eggs and juveniles .
Studies on the embryonic development of the East Asian freshwater leech, A. lata have not been conducted yet. The present study aims to describe the embryonic development of the glossiphoniid A. lata under laboratory conditions. In addition, lineage tracer injection, immunostaining and gene expression experiments were performed to support the use of this species as an “Evo-Devo” model organisms in the future.
Comparison between brooding periods of several Glossiphoniidae species
Alboglossiphonia hyalina 
Alboglossiphonia polypompholyx 
Batracobdella algira 
Glossiphonia complanata 
Oligobdella biannulata 
Theromyzon tessulatum 
Helobdella striata 
Helobdella austinensis 
Helobdella robusta 
Helobdella stagnalis 
Cleavage (stages 1 to 6)
Germinal band formation (stages 7 to 8)
Organogenesis and hatching (stages 9 to 11)
Post-embryonic stage and parental care
The beginning of juvenile stage is marked by the exhaustion of the yolk from within the crop . Starting from mid Stage 11, A. lata grows distinctively flat (Fig. 5a-c). Body coloration of the growing adults is fawn to pale, slightly translucent with a body length of 10–22 mm . Parental care in A. lata ends after the juveniles that have exhausted their yolk start leaving the adult. After brooding is finished, the parent adults feed on snails and reproduce at least two times more before dying.
Calsensin expression patterns
DNA markers from mitochondrial genomes, like cytochrome C oxidase subunit 1 (CO1), are widely used for estimating phylogenetic relationships among closely allied taxa . In the present study, the consensus tree generated by Neighbor-Joining method has confirmed the phylogenetic relationships within the Glossiphoniidae family established by Siddall and collaborators , in which A. lata species cluster alongside G. complanata in a monophyletic group that shares the common trail of presenting three pairs of eyespots. Our results are supported by a pairwise-distance matrix (Additional file 2) .
The developmental process observed in A. lata embryos is overall similar to that observed in the model organisms H. austinensis and other studied glossiphoniidae leech species. Aggregated groups of mature individuals were similar to those observed in A. hyalina , suggesting that hypodermic insemination  occurs in A. lata, although this was not directly observed. Nagao (1958) observed that the cocoon of A.lata is secreted from the clitellar glands surrounding the female gonopore , which is also the case for at least one Helobdella species . Like all glossiphoniids, A. lata develops an embryonic attachment organ [22, 31]. This organ, formed at the later part of stage 8, appears to be more prominent than that of other glossiphoniids such as Helobdella. Although microinjection experiments in A.lata embryos were ultimately successful, their vitelline membrane proved to be more difficult to puncture than that of H. austinensis embryos. To overcome this issue, better injection skills and the making of more resistant and sharper needles were required. At the same time, it was observed that injected A. lata embryos showed more resistance to bacterial infection in comparison to injected Helobdella embryos. The vitelline membrane of A. lata being thicker than that of H. austinensis could have helped decreasing the risk of infection of A. lata embryos after injection. Hatching of A. lata appears to be delayed relative to other glossiphoniid species outside of Alboglossiphonia . The eye spots, midgut diverticula and at least the posterior sucker are well differentiated prior to hatching (Fig. 43Be, Bf). Allowing the embryo to fully develop inside the vitelline membrane could potentially increase the chances of survival in case the embryos were forcefully or accidentally detached from the ventral surface of the parent adult before hatching, in comparison to those species that hatch earlier. We speculate that this delay could have been an evolutionary advantage for the Alboglossiphonia genus. Regarding reproduction strategy, A. lata leeches have shown to be iteroparous, with at least two reproductive cycles before dying. This strategy is known in other glossiphoniid genera . Interestingly, several studied Alboglossiphonia species [17, 18, 33, 34] and at least one Helobdella  and Batracobdella  species are, in contrast, semelparous. It is important to note that our observations occurred under experimental conditions, and therefore, it is recommended to elaborate a methodology for future confirmatory studies in A. lata specimens in situ.
The present study, for detailed comparative purposes with current leech model organisms species (Additional file 3, Table 1), followed the stage division known for Helobdella , in which eleven stages are established. However, in similar embryogenesis studies conducted for other leech species, this stage division varied [17, 34].
Calsensin is a EF-hand calcium-binding protein that was first found in the leech Haemopis marmorata, and is thought to mediate calcium-dependant signal transduction events in growth cones and axones of the developing nervous system . The spatial expression of an EF-hand calcium-binding protein Calsensin ortholog in A. lata (Ala-calsensin) has been characterized, appearing to be expressed in the segmental ganglia and peripheral neurons in the body wall during organogenesis (10 and 11). Considering that Calsensin expression has been detected in central and peripheral nerves of other hirudinid species [36, 37], our results give further support for a potential physiological role of Calsensin in the formation and maintenance of nerve pathways in leech species.
Successful injection of lineage tracing, visualization of neurogenesis during later stages by immunostaining using anti-acetyl tubulin antibody, and spatial expression pattern-based characterization of Calsensin by chemical in situ hybridization support the use of A. lata as a model organisms for “Evo-Devo” studies.
Description of the embryonic development of A. lata in vitro provides new data for further comparative studies involving other leech species. In addition, successful use of molecular biological techniques, such as microinjection of embryos for lineage tracing, in situ hybridization for spatial gene expression, and immunostaining for neurogenesis offers a basis for the development of this leech species as an “Evo-Devo” model organisms in the future.
Adult specimens, collected in nearby ponds and purchased online from Yeosu Aquarium, were bred in the Laboratory of Cellular and Developmental Biology (LCDB) of the Department of Biology of Chungbuk National University, Republic of Korea. Following the Protocol for Handling of Helobdella (Leech) embryos , Alboglossiphonia lata adult specimens (body length at rest: 10 – 18 mm) were deposited in Petri dishes with lid containing Working Solution. The specimens were cleaned once a day by changing the culture medium and scrubbed manually to get rid of any residual waste, and kept in an incubator at 22 °C. Their diet consisted of red snails purchased online, which were bred in fish bowls with Working Solution at room temperature.
Gravid adults were carefully manipulated with blunt forceps in order to remove the cocoons adhered to their ventral body wall. With the help of a sterilized pipette, each cocoon was transported to a separate smaller petri dish containing HTR medium for further examination and culture.
Using sterilized insect pins. Embryos were detached from their respective cocoons and observed through a Leica ZOOM 2000 stereomicroscope to identify their current developmental stage. Expecting a similarity with the timing of developmental stages in Helobdella species, A. lata embryos were checked every half an hour since deposition until it reached stage 4, then they were checked at least twice a day for the following stages. Embryos were imaged at each developmental stage using a Nikon SMZ18 stereomicroscope.
In order to elaborate a timeline for A. lata developmental stages that can be compared to the existing ones for Helobdella robusta and Helobdella austinensis, each one of the stages was indicated in terms of the number of hours after zygote deposition (AZD).
CO1 gene cloning and sequencing
Total RNA from Alboglossiphonia lata embryos was isolated using TRIzol (Ambion). Then mRNA from RNA using Oligo (dT) primer (Promega) was selected, and reverse transcription into cDNA (SuperScript II First-Strand Synthesis System for RT-PCR, Invitrogen) was conducted. CO1 protein coding homologous sequences were searched on sequenced RNA database available in our laboratory of Cellular and Developmental Biology (LCDB). Specific primers (CO1_Contig2 Forward: 5′ GCAGTGAAATATGCTCGGGT 3′, CO1_Contig2 Reverse: 5′ GAGTTAGCACAACCAGGCTCA 3′) were designed in order to amplify CO1 from A. lata cDNA. We followed the TaKaRa protocol for PCR according to standard procedure.
In order to acquire CO1 sequences, we cited the Leech gene sequencing articles available [29, 39] and used the accession numbers from said articles, with the exception of the Hirudo medicinalis CO1 sequence, because it had many gaps and did not qualify. Instead, we used the nucleotide sequence under the accession number AY786458 from NCBI that does not present gaps and is shorter than the previously referenced sequence. We searched additional accession numbers not present in these articles in NCBI. The sequences were aligned and trimmed using biological sequence editor BioEdit (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). Aligned sequences were analyzed in MEGA7.
Calsensin gene identification, gene cloning, probe synthesis, whole-mount in situ hybridization and nuclear staining
Total RNA was isolated from Alboglossiphonia species (lata) embryos of different developmental stages using TRIzol (Ambion). We selected mRNA from RNA using Oligo (dT) primer (Promega), and then conducted reverse transcription into cDNA (SuperScript II First Synthesis System for RT-PCR, Invitrogen). To demonstrate the feasibility of molecular approaches to A. lata, we isolated an EF-hand motif Calsensin gene, conducted alignment using alignment tool ClustalW, and phylogenetic analysis using tool MEGA7. After confirmation, the investigated leech gene was isolated by means of PCR, using gene-specific primers (Ala-calsensin Forward: 5′ GCCAAACGTTACCGAACCTCG 3′; Ala-calsensin Reverse: 5′ GAGAAGGTCCGCGTTGGCG 3′) based on sequenced RNA database available in our laboratory of Cellular and Developmental Biology (LCDB). The amplified fragments were cloned into pGEM T vector (Promega). Dioxigenin-labelled RNA probe were synthesized from the cloned fragments. Then in situ Hybridization (ISH) was performed as previously described . Pre-hybridization was performed at 64.7 °C for one day in hybridization buffer (50% Formamide, 5× SSC, 1× Denhardt’s Solution, 0.1% CHAPS, 100 mg/ml Heparin, 0.1% Tween20, 100 mg/ml tRNA). The prehybridized buffer was replaced with fresh hybridization buffer containing 2 ng/ml of the corresponding probe and embryos were hybridized at 64.7 °C for 2 days. Washed embryos were incubated at room temperature for 2 h in 1% Blocking Regent dissolved in PBT (1× PBS plus 0.1% Tween20) then incubated at 4 °C for 16 h with 1/1000 Anti-DIG/AP in 1% Blocking Reagent. After incubation, the color reaction was carried out using nitro blue tetrazolium chloride/ 5-bromo-4-chloro-3-indoyl-phosphate (Roche) by standard procedures. Stained embryos were dehydrated in ethanol, mounted in plastic embedding solution (PolyBed, Roche), and examined by bright field microscopy on a Nikon SMZ18 stereomicroscope.
Main cells were injected using dextran,tetramethylrhodamine (Molecular probes, D1817) for colorizing red. To visualize the different lineage, we injected green color fluorescence dye, the dextran Alexa fluor 488 (Molecular probes, D22910) in left OPQ cell at stage 6a when cleaved NOPQ to N and OPQ. After injection, the embryo were incubated at 22 °C in antibiotic (Gibco, 15,240,062) treated Helobdella triserialis media. After reaching the desired stage, embryos were fixed and treated with DAPCO glycerol. We labeled the embryo’s bright field color in pseudo blue color to pinpoint the exact location of the embryo. Embryos were imaged by fluorescent microscopy on a Nikon SMZ18 Stereomicroscope.
After rehydrating the embryos (stage 9 to 11), they were pre-incubated in 5% mercaptoethanol and 1% Triton in 0.1 M Tris-HCl (pH 7.5) at 37 °C on shaking incubator (rpm60) for an hour. Following three washes with PBT, the embryos were incubated in Block solution (1:9 10X Roche Western Blocking Reagent in PBT) for two hours. Then, embryos were incubated with a monoclonal anti-acetylated-α-Tubulin antibody (Sigma, T-7451) in Blocking Solution (1:500) at 4 °C for 72 h. After three consecutive washes with PBT, embryos were incubated with a secondary antibody (Abcam, ab150113) in Blocking Solution (1:1000) at 4 °C for 48 h. Consequently, embryos were washed overnight with PBT and then dyed with DAPI in PBT (1:1000) at room temperature in the dark for overnight. After washing with PBT three more times, embryos were finally embedded in 30%, 50% 20 min and 87% glycerol and 2.5 mg/ml of DABCO in 1xPBS. Embryos were imaged by fluorescence microscopy on a Nikon SMZ18 Stereomicroscope and LEICA DM6 B.
We thank members of the Cho Laboratory for valuable comment.
This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Development of Environmental Exosome-reactor and the relationship between Bio/Chemical hazards and exoxomal metabolites, PJ012653)” Rural Development Administration, Republic of Korea.
Availability of data and materials
BIMJ and SJC conceived the project and designed the experiments. BIMJ documented the embryonic development of A. lata embryos by microscopy, performed cloning and sequencing of CO1, elaborated the embryogenesis time table for this species, and was a major contributor in writing the manuscript. HJK performed lineage tracer injection in meso and ectodermal precursor cells, immunostaining in late stage embryos, in situ hybridization and phylogenetic tree construction for Ala-Calsensin gene, and analyzed the results obtained. JSP performed phylogenetic tree analysis. BIMJ and SJC wrote the manuscript, JWK improve the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Aguinaldo AM, Turbeville JM, Linford LS, Rivera MC, Garey JR, Raff RA, Lake JA. Evidence for a clade of nematodes, arthropods and other moulting animals. Nature. 1997;387:489–93.View ArticlePubMedGoogle Scholar
- Ruiz-Trillo I, Riutort M, Littlewood DT, Herniou EA, Baguna J. Acoel flatworms: earliest extant bilaterian metazoans, not members of Platyhelminthes. Science. 1999;283:1919–23.View ArticlePubMedGoogle Scholar
- Weisblat DA, Huang FZ. An overview of glossiphoniid leech development. Can J Zool. 2001;79:218–32.View ArticleGoogle Scholar
- Kuo DH. Annelids in evolutionary developmental biology. In: Shain DH, editor. Annelids in Modern Biology. New Jersey: John Wiley & Sons Inc; 2009. p. 65–87.Google Scholar
- Kutschera U, Langguth H, Kuo DH, Weisblat DA, Shankland M. Description of a new leech species from North America, Helobdella austinensis n. Sp. (Hirudinea: Glossiphoniidae), with observations on its feeding behaviour. Zoosyst Evol. 2013;89:239–46.View ArticleGoogle Scholar
- Kutschera U, Weisblat DA. Leeches of the genus Helobdella as model organisms for Evo-Devo studies. Theory Biosci. 2016;134(3–4):93–104.Google Scholar
- Oka A. Synopsis der japanischen Hirudineen, mit Diagnosen der neuen Spezies. Annot Zool Jpn. 1910;7:165–83.Google Scholar
- Sawyer RT. Leech biology and behavior, vol. 1–3. Oxford: Clarendon Press; 1986.Google Scholar
- Lai YT, Chen JH. Leech fauna of Taiwan. Taipei: National Taiwan University Press; 2010.Google Scholar
- Lee DJ, Yoon CS, Lee JC, Sung SH, Park DR, Cheong SW. Fundamental investigation for long-term ecological monitoring on community of benthic macroinvertebrates in wetland Woopo. J Environ Sci Int. 2009;18(12):1399–410.View ArticleGoogle Scholar
- Sawyer RT. The phylogenetic development of brooding behaviour in the Hirudinea. Hydrobiologia. 1971;37:197–204.View ArticleGoogle Scholar
- Kutschera U, Wirtz P. The evolution of parental care in freshwater leeches. Theory Biosci. 2001;120:115–37.View ArticleGoogle Scholar
- Encyclopedia of Life. http://www.eol.org. Accessed 5 Sept 2017.
- Bissen ST, Weisblat DA. The durations and compositions of cell cycles in embryo of the leech Helobdella triseralis. Development. 1989;106:105–18.PubMedGoogle Scholar
- Fernandez J, Stent GS. Embryonic development of the glossiphoniid leech Theromyzon rude: structure and development of the germinal bands. Dev Biol. 1980;78:407–34.View ArticlePubMedGoogle Scholar
- Huang FZ, Kang D, Ramirez-Weber FA, Bissen ST, Weisblat DA. Micromere lineage in the glossiphoniid leech Helobdella. Development. 2002;129:719–32.PubMedGoogle Scholar
- Ben Ahmed R, Gammoudi M, Tekaya S. Life cycle data and development of the freshwater leech Alboglossiphonia hyalina (Annelida, Clitellata, Hirudinea) in Tunisia. Acta Zool Bulgar. 2015;67(2):265–70.Google Scholar
- El-Shimy NA, Davies RW. The life-cycle, ecology and host specificity of the freshwater leech Alboglossiphonia polypompholyx (Glossiphoniidae) in Egypt. Hydrobiologia. 1991;222:173–8.View ArticleGoogle Scholar
- Romdhane Y, Ben Ahmed R, Tekaya S. Insemination and embryonic development in the leech: Batracobdella algira (Hirudinea, Annelida). Invertebr Reprod Dev. 2014;59(1):1–9.Google Scholar
- William E, Moser R Van, Devender W, Klemm DJ. Life history and distribution of the Leech Oligobdella biannulate (Moore, 1900) (Euhirudinea: Glossiphoniidae). Comp Parasitol. 2005;72:17–21.Google Scholar
- Tarin AM, Charlene WS, Tracie LP, Jonelle C, Daniel HS. Cocoon deposition and hatching in the Aquatic Leech; Theromyzon tessulatum (Annelida, Hirudinea, Glossiphoniidae). Am Midl Nat. 2005;154:78–87.View ArticleGoogle Scholar
- Kutschera U, Wirtz P: Reproduction behaviour and parental care of Helobdella striata (Hirudinea: Glossiphoniidae): a leech that feeds its young. Ethology. 1986b;72:132–42.Google Scholar
- Romdhane Y, Ben Ahmed R, Tekaya S. Sexual behavior, insemination and development of the freshwater leech Helobdella stagnalis (Annelida, Hirudinea, Glossiphoniidae). Invertebr Reprod Dev. 2017; doi:10.1080/07924259.2017.1340354.
- Mann KH. Leeches (Hirudinea): Their structure, physiology, ecology, and embryology. New York: Pergamon Press; 1962.Google Scholar
- Whitman CO. The embryology of Clepsine. Q J Microsc Sci. 1878;18:215–315.Google Scholar
- Clifford HF. Aquatic invertebrates of Alberta. Edmonton: University of Alberta Press; 1991.Google Scholar
- Weisblat DA, Kuo DH. Developmental biology of the leech Helobdella. Int J Dev Biol. 2014;58:429–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Marshall E. Taxonomy. Will DNA bar codes breathe life into classification? Science. 2005;307:1037.Google Scholar
- Siddall ME, Budinoff RB, Borda E. Phylogenetic evaluation of systematics and biogeography of the leech family Glossiphoniidae. Invertebr Syst. 2005;19:105–12.View ArticleGoogle Scholar
- Hall BG. Phylogenetic tress made easy: A how-to manual for molecular biologists. Maryland: Sinauer Association Inc.; 2011.Google Scholar
- Nagao Z. Some observations on the breeding habits in a freshwater leech, Glossiphonia lata Oka. Jap J Zool. 1958;12:219–28.Google Scholar
- Coleman J, Shain DH. Clitellate cocoons and their secretion. In: Shain DH, editor. Annelids in Modern Biology. New Jersey: Wiley; 2009. p. 328–346.View ArticleGoogle Scholar
- Proctor RM, Young JO. The life history, diet and migration of a lake-dwelling population of the leech Alboglossiphonia heteroclita. Hydrobiologia. 1987;150:133–9.View ArticleGoogle Scholar
- Gouda HA. A new Alboglossiphonia species (Hirudinea: Glossiphoniidae) from Egypt: Description and life history data. Zootaxa. 2010;2361:46–56.View ArticleGoogle Scholar
- Briggs KK, Silvers AJ, Johansen KM, Johansen J. Calsensin: A novel calcium-binding protein expressed in a subset of peripheral leech neurons fasciculating in a single axon tract. J Cell Biol. 1995;129(5):1355–62.View ArticlePubMedGoogle Scholar
- Macagno ER, Stewart RR, Zipser. The expression of antigens by embryonic neurons and glia in segmental ganglia of the leech Haemopis marmorata. J Neurosci. 1983;3:1746–59.PubMedGoogle Scholar
- Veldman M, Huang Y, Jellies J, Johansen KM, Johansen J. Differential expression of the EF-hand calcium-binding protein calsensin in the central nervous system of hirudinid leeches. Cell Tissue Res. 1996;286(3):357–64.View ArticlePubMedGoogle Scholar
- Weisblat DA, Kuo DH. Helobdella (Leech): a model for developmental studies. In: Gann A, Crotty D, editors. Emerging Model Organisms: a Laboratory Manual, Volume 1. New York: Cold Spring Harbor Laboratory Press; 2009. p. 245–267.Google Scholar
- Apakupakul K, Siddall ME, Burreson EM. Higher level relationships of leeches (Annelida: Clitellata: Euhirudinea) based on morphology and gene sequences. Mol Phylogenet Evol. 1999;12(3):350–9.View ArticlePubMedGoogle Scholar
- Cho SJ, Valles Y, Giani VC Jr, Seaver EC, Weisblat DA. Evolutionary dynamics of the wnt gene family: a lophotrochozoan perspective. Mol Biol Evol. 2010;27:1645–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.PubMedGoogle Scholar
- Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.View ArticlePubMedGoogle Scholar
- Zuckerkandl E, Pauling L. Evolutionary divergence and convergence in proteins. In: Edited by Bryson V, Vogel HJ, editors. Evolving genes and proteins. New York: Academic Press; 1965. p. 97-166. Google Scholar
- Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci. U.S.A. 2004;101:11030–5.Google Scholar
- Kumar S, Stecher G, Tamura K. MEGA7. Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.View ArticlePubMedGoogle Scholar