Abstract
In the present study, we investigated the possible existence of new lineages of peripatids through comparisons between known Neotropical species and specimens obtained from two locations in Pará, a state in eastern Brazilian Amazonia using a molecular approach based on sequences of the mtDNA genes COI, 16Sr RNA, and 18S RNA. The analyses included also sequences of Asian and African taxa for a more systematic understanding of the phylogenetic relationships within the group. The analysis of the COI, 16S rRNA and 18S RNA sequences permitted the identification of three distinct lineages (A, B and C) based on two different phylogenetic approaches (Bayesian methods and ML). The three lineages presented here are completely distinct from all other peripatid taxa so far defined by molecular data. The presence of specimens of three independent onychophoran lineages occurring in sympatry in the Amazon basin was confirmed in all the analyses, providing consistent support for the phylogenies presented in this study. These findings reinforce the importance of the Amazon region in the diversification of Neotropical peripatids, and indicate that onychophoran diversity is much greater than previously thought, given that the number of taxa found at a single site was equivalent to the total number of allopatric species described for the entire region.
Keywords: Amazonia; new lineages; Peripatidae; phylogeny; Onychophora
Velvet worms are soft-bodied terrestrial invertebrates that inhabit humid forests, where they are typically associated with the soil, decaying trunks, and leaf litter, as well as caves (Vasconcellos et al., 2004). These worms are usually confined to stable microhabitats with high levels of humidity, and have limited dispersal capabilities in open environments (New, 1995). Currently, the phylum Onychophora encompasses two widely-distributed but allopatric families: Peripatidae Evans, 1901, found in Central and South America, West Africa, and Southeast Asia; and Peripatopsidae Bouvier, 1905 which is restricted to South Africa, Australasia, and South America (Oliveira et al., 2012a).
The family Peripatidae is distributed in 12 genera, with a total of 85 species currently recognized (Oliveira et al., 2012a,b, 2013, 2014). Four genera of this family are found in the Amazon basin, where only seven species have been identified so far: Epiperipatus brasiliensis (Bouvier, 1889); Epiperipatus edwardsii (Blanchard, 1847); Epiperipatus simoni (Bouvier, 1899); Epiperipatus tucupi (Froehlich, 1968), designated nomen dubium by Oliveira et al. (2012a); Macroperipatus geayi (Bouvier, 1889); Oroperipatus balzani (Camerano, 1897); and Oroperipatus eisenii (Wheeler, 1898), although Sampaio-Costa et al. (2009) have also described a morphospecies of the genus Peripatus. While onychophorans are known to occur in the Amazon basin (Read, 1988; Sampaio-Costa et al., 2009), the diversity of these terrestrial invertebrates in this region – considered “megadiverse” for many other groups of animals (Mittermeier et al., 2003) – is still poorly known, and the last species of velvet worm from the Amazon biome was described more than 50 years ago (Froehlich, 1968).
While there have been recent advances in the description of onychophoran species, more than 70% of the 85 recognized species of Peripatidae are in need of revision. In addition to the new diagnostic characteristics proposed by Oliveira et al. (2012b, 2013, 2014), genetic data obtained over the past few decades have been used increasingly to complement morphological analyses, contributing to the identification of cryptic species, the determination of species limits, phylogeny, and the distribution patterns of this invertebrate group (Lacorte et al., 2011; Oliveira et al., 2011; Murienne et al., 2013).
In the present study, we analyzed the sequences of two mitochondrial markers and one nuclear marker in specimens of peripatids collected from fragments of secondary forest in the eastern Brazilian Amazon basin. Based on the results, we investigated the possible existence of new taxa through comparisons with the Neotropical species already analyzed using this approach, for the definition of phylogenetic relationships. The findings also contribute to our knowledge of the natural history of this poorly-known group of Amazonian organisms.
The tissue samples were obtained from 22 specimens collected between 2006 and 2011 from two localities in Pará, a state in eastern Brazilian Amazonia. Specimen collection was authorized by the Brazilian federal environment institute (Sistema de Autorização e Informação em Biodiversidade - SISBIO / Instituto Chico Mendes de Conservação da Biodiversidade – ICMBio; special license 24714-1). The sample localities shown in Figure 1 and listed in Table S1 were: (7) Outeiro Island - OTI (1°14′56″ S, 48°26′57″ W) (n = 10 specimens) and (8) Bragança - BRG (1°02′04″ S, 46°45'48″ W) (n = 12 specimens), which are separated by a distance of 234 km (Figure 1). The first site (OTI) is a small patch of forest, of approximately 2 km2, located within the urban area of Outeiro Island. The habitat is well-shaded, with many decomposing tree trunks on the ground, and is interspersed with plantations of regional fruit species such as “açaí”, Euterpe oleraceae Mart. and “cupuaçu”, Theobroma grandiflorum (Willd. ex Spreng). The second site (BRG) is a fragment of secondary forest of less than 6 km2, with high levels of anthropogenic disturbance, located within the town's urban zone. Specimens were mostly found in rotten trunks, leaf litter, tree roots, and cavities in the ground. The animals were euthanized using the Monge-Nájera and Bernal (1994) protocol and preserved in 70% ethanol, prior to being deposited in the zoological collection of the Bragança campus of the Federal University of Pará, Brazil.
Map showing the geographic distribution of the species of Peripatidae and sites analyzed in the present study. Numbers correspond to locality records listed in Table S1 (Supplementary material). Colours refer to the clades based on the phylogenetic tree. Topology obtained from the ML analyses of mitochondrial COI. Above: asterisks indicated Bayesian posteriors probabilities > 0.9; below: numbers nodes are bootstrap values > 75%. Abbreviations: OTI = Outeiro Island; BRG = municipality of Bragança.
Total DNA was extracted from fragments of the tegument of the specimens using a phenol-chloroform protocol (Sambrook et al., 1989). Fragments of two mitochondrial loci, Cytochrome Oxidase I (COI) and the large Ribosomal Subunit (16S rRNA), as well as the nuclear ribosomal gene 18S (18S rRNA) locus, were amplified by polymerase chain reaction (PCR). For COI, the primers were L5584 and H6174 (Oliveira et al., 2011), whereas for 16S rRNA, they were L1987 and H2609 (Palumbi et al., 1991), and for the nuclear marker 18S rRNA, the primers were 1F and 5R, from Giribet et al. (1996).
The nucleotide sequences of the resulting gene fragments were determined in an ABI 3500 automatic sequencer (Applied Biosystems). In addition to the sequences obtained in the present study, COI, 16S rRNA and 18S rRNA sequences were obtained from GenBank for inclusion in the analyses (Table S1), representing 13 species of peripatids, 11 from the Neotropical region, one from tropical Africa and one from Asia, with the latter two being used as the outgroups, based on the results of Murienne et al. (2013).
The sequences were aligned in ClustalW (Thompson et al., 1994) using the Bioedit v7.0.5 sequence editor (Hall, 1999). The nucleotide composition and divergence rates (p distances) between and within the observed lineages were calculated in MEGA v6.0 (Tamura et al., 2013). The JModelTest v2.0.2 program (Guindon and Gascuel 2003; Darriba et al., 2012) was used to select the optimal evolutionary model for the phylogenetic analysis of the sequences of both mtDNA and ncDNA regions using the Akaike and Bayesian information criteria (AIC and BIC).
The maximum likelihood (ML) analysis was run in PhyML v3.0 (Guindon and Gascuel, 2003). Support for the groups was evaluated using a bootstrap approach with 1000 replicates. Bootstrap support values of more than 75% were considered to represent a well-defined group. The evolutionary Bayesian inference (BI) models produced by the MrBayes v3.2.0 program (Ronquist and Huelsenbeck, 2003) were selected based on the Bayesian Information Criterion (BIC). The Bayesian Monte Carlo Markov Chain (BMCMC) procedure was based on four chains run simultaneously over 107 generations, with samples being taken every 100 generations. Convergence and effective sample sizes (ESS) were assessed in Tracer v1.6.1 (Rambaut and Drummond, 2007), and the first 100 trees of each run were discarded as burn-in. Bayesian posterior probability values lower than 0.9 were considered to be inconclusive. The trees were visualized in FigTree v.1.4 (Rambaut, 2012).
The best model of nucleotide substitution and the partitioning of the concatenated data set were selected using Partition Finder, v. 1.1.1 (Lanfear et al., 2012), based on the Akaike Information Criterion, or AIC (Akaike, 1973). The phylogenetic trees derived from the ML analysis was estimated using RaxML, v. 8.0 (Stamatakis, 2014), with the optimal partitioning for this analysis and the support for each branch node being calculated using a nonparametric bootstrap analysis, with 1000 pseudo-replicates (Felsenstein, 1985), which also provides an estimate of the confidence for the results. The trees were visualized in FigTree v.1.4 (Rambaut, 2012).
It was not possible to sequence all the genes for each of the specimens collected (Table S1). The sequencing of the COI gene produced fragments of 501 base pairs (bps) with 183 variable sites for 18 of the specimens collected for this study. The 16S rRNA was sequenced in 13 specimens, producing fragments of 392 bps, with 97 variable sites. Fragments of the 18S rRNA were amplified in ten specimens, providing 775 bps and 149 variable sites. There was a predominance of the A+T nucleotides in all the markers, which is typical of both onychophorans and most other invertebrates (Trewick, 2000; Jeon et al., 2012).
The GTR+ I + G model was selected for the COI sequences for both probabilistic approaches (ML and BI). The K81uf+G, TVM+G, and HKY+G models were selected for the ML analysis based on the concatenated COI, 16S rRNA, and 18S rRNA sequences, respectively. The phylogenetic analyses based on the COI, 16S rRNA and 18S rRNA sequences permitted the identification of three distinct lineages from the Amazon region. The first two lineages, denominated A and B, are made up of specimens from both Outeiro Island and Bragança, while the third, lineage C, was formed exclusively by specimens from Outeiro Island (Figure 1; Figure S1). The topology derived did not provide statistically support for any phylogenetic relationship between the species with sequences available and the lineages identified in the present study. Even though, all the species and lineages represented by more than one specimen were validated with high levels of statistical support, confirming the capacity of the markers to distinguish valid taxa in the onychophorans.
We also confirmed that the three clades identified in the present study do not form a single group. Lineage A appears to be close to E. edwardsii, whereas lineages B and C form a sister group to the Caribbean Peripatus dominicae basilensis, although in all cases, with reduced statistical support (Figure 1). The data also indicate that the species of the genera Epiperipatus and Peripatus are non-monophyletic.
The genetic divergence found between the species and the new Neotropical lineages varied from 4.9% to 22.6% when the outgroups were included, and from 4.9% to 20.6% when only the ingroup was considered (Table 1). These values were obtained for the COI marker, which permits the greatest number of comparisons due to the large number of species for which published data are available. The minimum interspecific divergence observed between Amazonian taxa (lineages B and C) was 9.9%, while the maximum was 14.5%, between the A and C lineages. Intraspecific genetic divergence was relatively low, at 0.5% (lineage A), 0.4% (lineage B), and 0.2% (lineage C), similar to that recorded for the species Principapillatus hitoyensis (0.2%), and lower than the divergence found in E. machadoi, E. diadenoproctus, E. paurognostus and E. adenocryptus (1.4%, 1.1%, 1% and 2%) (Oliveira et al., 2011), which indicates low levels of intraspecific variation in the sequences of the specimens analyzed. In the case of the 16S rRNA gene, divergence was between 9% and 12% (including outgroups), with a minimum of 9% between the A and B lineages and 11% between B and C. Intraspecific genetic divergence was 0.4% (lineage A), 0.0% (lineage B), and 0.2% (lineage C). Based on the 18S rRNA gene, the minimum distance was 1.1%, between lineages A and C, while the maximum was 2.2%, between lineages B and C.
The p-distances recorded within and between the lineages and species analyzed in the present study based on sequences of the mitochondrial COI gene.
The comparative analysis of the specimens, together with the other peripatid taxa for which molecular data are available, indicates that they represent three completely distinct lineages. The phylogenies based on the mitochondrial and nuclear sequences, and the marked divergence found between the lineages identified in the analyses (9-13%) indicate the existence of distinct species (Table 1), considering the phylogenetic species concept (sensuMishler and Theriot, 2000). In fact, a genetic distance of only 4.4% was considered diagnostic of the species-level differentiation of the allopatric peripatids E. adenocryptus and E. paurognostus (Oliveira et al., 2011).
As the type localities of three onychophoran species, M. geayi, E. brasiliensis, and E. tucupi (Sampaio-Costa et al., 2009), are located in the eastern Amazon basin, it seems likely that at least one of the lineages identified in the present study (which need to be described formally) may correspond to one of these species. However, as the type-specimen for these species are by now too degraded to provide material for genetic analysis and no other specimens are available from the type localities, it is impossible to provide a direct comparison using molecular tools. In addition, there is no information on the exact geographical location of the type localities, which impedes the collection of new specimens. Clearly, it will be necessary to examine the specimens analyzed in the present study very meticulously for the identification of diagnostic traits in order to confirm their potential species status and avoid synonymy.
The present analysis of mitochondrial and nuclear markers did not provide a well-resolved arrangement (Figure 1; Figure S1). The Amazonian lineages do not form a single clade when analyzed in comparison with the other species for which data are available on the same molecular parameters. This finding contrasts with the situation observed in five allopatric species of Epiperipatus from the state of Minas Gerais, Brazil (Oliveira et al., 2011). While the latter are separated by distances of between 11 km and 155 km, they are nevertheless phylogenetically closely related, forming a single group with a common ancestor.
The sympatric occurrence of these phylogenetically distinct lineages also indicates that this group of animals has undergone distinct spatial-temporal differentiation processes, which have molded species ranges and their diversity in this biome, as observed in other groups of organisms, reflecting the complex zoogeographic and cladogenetic processes that are typical of the Amazon biome. While these animals are restricted to humid habitats, the different lineages are probably adapted to distinct ecological conditions. In the Blue Mountains of Australia, for example, the sympatric species Cephalophovea tomahmontis Ruhberg, Tait, Briscoe and Storch, 1988 and Euperipatoides leuckartii (Sänger, 1871) present quite distinct life history strategies (Leishman and Eldridge, 1990). A high level of interspecific variation (9.22%) has been observed in sympatric peripatopsids, between the Peripatopsis moseleyi (Wood-Mason, 1879) and Peripatopsis balfouri (Sedgwick, 1885) species complexes (Daniels and Ruhberg, 2010; Daniels et al., 2013).
An interesting aspect of the results of our study is the occurrence in Bragança of two of the three lineages found on Outeiro Island, 234 km to the west (Figure 1). A probable scenario is that the lineages were more amply distributed prior to the formation of the island, with the present-day distribution attesting to the ancient connectivity of these environments. Given this situation, the restriction of lineage C to Outeiro Island may reflect a sampling effect rather than the presence of an endemic taxon on the island, which was not identified as a vicariant factor in the establishment of any of the lineages. It would thus be important to expand the number of points sampled on the mainland in order to confirm the more ample distribution of all three lineages.
The small number of available specimens and sample localities are insufficient for a more conclusive analysis of geographic limits or distribution patterns of the three lineages, although they do confirm their occurrence. These findings reinforce the importance of the Amazon region in the diversification of the Neotropical peripatids, and indicate that onychophoran diversity is much greater than previously thought, given that the number of taxa found at a single site was equivalent to the total number of allopatric species described for the entire region (Sampaio-Costa et al., 2009). Similarly, understanding how the velvet worms colonized the region, and which barriers contributed to their diversification, may provide important insights into speciation patterns in the Amazon basin, given that onychophorans have low vagility and are sensitive to environmental impacts. These characteristics may favor the isolation of these organisms, making them an appropriate model for the analysis of biogeographic patterns (Murienne et al., 2013).
Acknowledgments
This study was part of the MS dissertation of WTRC, which was financed by the Amazon Research Foundation - FAPESPA (Scholarship 012/2011) and by the Brazilian National Research Council - CNPq (grants 306233/2009-6, 473341/2010-7 to IS, and 305645/2009-9 to HS). We thank Stephen Ferrari for revision of the manuscript. We also would like thank the anonymous reviewers for their helpful comments.
References
- Akaike H (1973) Information theory and an extension of the maximum likelihood principle. In: Petrov BN and Csaki F (eds) Second International Symposium on Information Theory. Akademiai Kiado, Japan, pp 267-281.
- Daniels SR and Ruhberg H (2010) Molecular and morphological variation in a South African velvet worm Peripatopsis moseleyi (Onychophora, Peripatopsidae): Evidence for cryptic speciation. J Zool 282:171-179.
- Daniels SR, McDonald DN and Picker MD (2013) Evolutionary insight into the Peripatopsis balfouri sensu lato species complex (Onychophora: Peripatopsidae) reveals novel lineages and zoogeographic patterning. Zool Scr 42:656-674.
- Darriba D, Taboada GL, Doallo R and Posada D (2012) jModelTest 2: More models, new heuristics and parallel computing. Nat Methods 9:772.
- Felsenstein J (1985) Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39:783-791.
- Froehlich CG (1968) On some Brazilian onychophores. Beitr Neotrop Fauna 5:160-171.
- Giribet G, Carranza S, Baguñà J, Riutort M and Ribera C (1996) First molecular evidence for the existence of a Tardigrada + Arthropoda clade. Mol Biol Evol 13:76-84.
- Guindon S and Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52:696-704.
- Hall TA (1999) BIOEDIT: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp 41:95-98.
- Jeon MJ, Song JH and Ahn KJ (2012) Molecular phylogeny of the marine littoral genus Cafius (Coleoptera: Staphylinidae: Staphylininae) and implications for classification. Zool Scr 41:150-159.
- Lacorte GA, Oliveira IS and Fonseca CG (2011) Phylogenetic relationships among the Epiperipatus lineages (Onychophora: Peripatidae) from the Minas Gerais State, Brazil. Zootaxa 2755:57-65.
- Lanfear R, Calcott B, Ho SYW and Guindon S (2012) PartitionFinder: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol 29:1695-1701.
- Leishman MR and Eldridge DB (1990) Life history characteristics of two sympatric onychophoran species from the blue mountains, New South Wales. Proc Linn Soc New South Wales 112:173-185.
- Mishler BD and Theriot EC (2000) The phylogenetic species concept (sensu Mishler and Theriot): Monophyly, apomorphy, and phylogenetic species concepts. In: Wheeler QD and Meier R (eds) Species Concepts and Phylogenetic Theory: A Debate. Columbia University Press, New York, pp 44-54.
- Mittermeier RA, Mittermeier CG, Brooks TM, Pilgrim JD, Konstant WR, da Fonseca GAB and Kormos C (2003) Wilderness and biodiversity conservation. Proc Natl Acad Sci U S A 100:10309-10313.
- Monge-Nájera J and Bernal MB (1994) Morphological and physiological characteristics of two species of Epiperipatus from Costa Rica (Onychophora: Peripatidae). Rev Biol Trop 42:181-188.
- Murienne J, Daniels SR, Buckley TR, Mayer G and Giribet G (2013) A living fossil tale of Pangaean biogeography. Proc R Soc Lond B Biol Sci 281:20132648.
- New TR (1995) Onychophora in invertebrate conservation: Priorities, practice and prospects. Zool J Linn Soc 114:77-89.
- Oliveira IS, Lacorte GA, Fonseca CG, Wieloch AH and Mayer G (2011) Cryptic speciation in Brazilian Epiperipatus (Onychophora: Peripatidae) reveals an underestimated diversity among the peripatid velvet worms. PLoS One 6:e19973.
- Oliveira IS, Read VMSJ and Mayer G (2012a) A world checklist of Onychophora (velvet worms), with notes on nomenclature and status of names. Zookeys 211:1-70.
- Oliveira IS, Franke FA, Hering L, Schaffer S, Rowell DM, Weck-Heimann A, Monge-Nájera J, Morera-Brenes B and Mayer G (2012b) Unexplored character diversity in Onychophora (velvet worms): A comparative study of three peripatid species. PLoS One 7:e51220.
- Oliveira IS, Schaffer S, Kvartalnov PV, Galoyan EA, Palko IV, Weck-Heimann A, Geissler P, Ruhberg H and Mayer G (2013) A new species of Eoperipatus (Onychophora) from Vietnam reveals novel morphological characters for the South-East Asian Peripatidae. Zool Anz 254:495-510.
- Oliveira IS, Lacorte GA, Weck-Heimann A, Cordeiro LM, Wieloch AH and Mayer G (2014) A new and critically endangered species and genus of Onychophora (Peripatidae) from the Brazilian savannah - A vulnerable biodiversity hotspot. Syst Biodivers 3:211-233.
- Palumbi S, Martin A, Romano S, McMillian WO, Stice L and Grabowiski G (1991) The Simple Fool's Guide to PCR. University of Hawaii Press, Honolulu, 45 pp.
- Read VMSJ (1988) The application of scanning electron microscopy to the systematics of the Neotropical Peripatidae (Onychophora). Zool J Linn Soc 93:187-223.
- Ronquist F and Huelsenbeck JP (2003) Mrbayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19:1572-1574.
- Sambrook J, Fritsch E and Maniatis T (1989) Molecular Cloning: A Laboratory Manual. 2nd edition. Cold Spring Harbor Laboratory Press, New York, 1626 pp.
- Sampaio-Costa C, Chagas-Junior A and Baptista RLC (2009) Brazilian species of Onychophora with notes on their taxonomy and distribution. Zoologia 26:553-561.
- Stamatakis A (2014) RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312-1313.
- Tamura K, Stecher G, Peterson D, Filipski A and Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30:2725-2729.
- Thompson JD, Higgins DG and Gibson TJ (1994) Clustal W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673-4680.
- Trewick SA (2000) Mitochondrial DNA sequences support allozyme evidence for cryptic radiation of new zealand Peripatoides (Onychophora). Mol Ecol 9:269-281.
- Vasconcellos A, Almeida WO and Eloy ECC (2004) Onychophora de florestas úmidas do complexo da mata atlântica do nordeste brasileiro e sua importância para conservação e estudos sistemáticos. In: Pôrto KC, Cabral JJP and Tabarelli M (eds) Brejos de Altitude: História Natural, Ecologia e Conservação. Ministério do Meio Ambiente, Brasília, pp 139-144.
Internet Resources
-
Rambaut A (2012) FigTree - Tree figure draw tool version 1.4. Available via http://tree.bio.ed.ac.uk/ (September 9, 2014).
» http://tree.bio.ed.ac.uk/ -
Rambaut A and Drummond AJ (2007) Tracer v1.4. Available via http://beast.bio.ed.ac.uk/Tracer (September 9, 2014).
» http://beast.bio.ed.ac.uk/Tracer
Supplementary material
The following online material is available for this article:
Figure S1 - Maximum-likelihood tree for three genes (COI, 16S, 18S).
Publication Dates
-
Publication in this collection
02 Mar 2017 -
Date of issue
Jan-Mar 2017
History
-
Received
21 Feb 2016 -
Accepted
27 Aug 2016