Next Article in Journal
Phylogeny and Biogeography of Branchipolynoe (Polynoidae, Phyllodocida, Aciculata, Annelida), with Descriptions of Five New Species from Methane Seeps and Hydrothermal Vents
Next Article in Special Issue
Three New Lizard Species of the Liolaemus montanus Group from Perú
Previous Article in Journal
Nuclear Orthologs Derived from Whole Genome Sequencing Indicate Cryptic Diversity in the Bemisia tabaci (Insecta: Aleyrodidae) Complex of Whiteflies
Previous Article in Special Issue
Conservation Status of Brachycephalus Toadlets (Anura: Brachycephalidae) from the Brazilian Atlantic Rainforest
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Barcoding Analysis of Paraguayan Squamata

1
Instituto de Investigación Biológica del Paraguay, Del Escudo 1607, 1425 Asunción, Paraguay
2
Guyra Paraguay, Av. Cnel. Carlos Bóveda, Parque Ecológico Capital Verde—Viñas Cué, 1719 Asunción, Paraguay
3
Subtropica, Teniente Rivas 741, 1732 Asunción, Paraguay
4
Senckenberg Forschungsinstitut und Naturmuseum Frankfurt, Senckenberganlage 25, 60325 Frankfurt a.M., Germany
5
Institute for Ecology, Evolution & Diversity, Goethe-University, Biologicum, Building C, Max-von-Laue-Str. 13, 60438 Frankfurt am Main, Germany
*
Authors to whom correspondence should be addressed.
Submission received: 26 June 2019 / Revised: 21 August 2019 / Accepted: 23 August 2019 / Published: 30 August 2019
(This article belongs to the Special Issue Systematics and Conservation of Neotropical Amphibians and Reptiles)

Abstract

:
Paraguay is a key spot in the central region of South America where several ecoregions converge. Its fauna (and specifically its herpetofauna) is getting better studied than years before, but still there is a lack of information regarding molecular genetics, and barcoding analyses have proven to be an excellent tool in this matter. Here, we present results of a barcoding analysis based on 16S rRNA gene sequences, providing valuable data for the scientific community in the region. We based our fieldwork in several areas of Paraguay. We analyzed 249 samples (142 sequenced by us) with a final alignment of 615 bp length. We identified some taxonomic incongruences that can be addressed based on our results. Furthermore, we identify groups, where collecting efforts and research activities should be reinforced. Even though we have some blanks in the geographical coverage of our analysis—and there is still a lot to do towards a better understanding of the taxonomy of the Paraguayan herpetofauna—here, we present the largest genetic dataset for the mitochondrial DNA gene 16S of reptiles (particularly, Squamata) from Paraguay, which can be used to solve taxonomic problems in the region.

1. Introduction

Ideally, profound knowledge of biodiversity is the first step before any conservation action, sustainable management or biological study is carried out in a given area [1,2]. The scientific community is not only aware that the world is facing a major extinction event [3], but also that many lineages are disappearing even before becoming known to science [4,5]. In the last decades, the use of molecular data has helped to improve our knowledge about biological diversity and the use of genetics as a tool for species recognition is now routine. Molecular data are more often used every day and are applied to the species identification [6,7,8] and species delimitation [9,10,11,12] of all kinds of living organisms, with applications even in food control quality [13].
In this context DNA barcoding analysis is the examination and comparison of short and stable fragments of DNA (DNA barcodes), usually mitochondrial, that represent genetic identifiers for a species [14], and have proven to be a reliable technique for taxonomy [15,16,17]. Nevertheless, researchers have to be cautious because despite being a power tool, barcoding analysis potentially can lead to misinterpretations if sequences used for comparison were generated from misidentified specimens [18]. Therefore, it is strongly recommended that molecular genetic tools are complemented with morphological, bioacoustics, and ecological data [19]. In conclusion, molecular genetics open a path for more detailed taxonomic studies [20].
For biological works that are concerned with the central portion of South America, Paraguay is critical since the country is located in a confluence zone of different ecoregions, such as Cerrado, Pantanal, Atlantic Forest, Chaco (Humid and Dry), and Southern Cone Mesopotamian Savanna, each of them having its own distinct origin and evolutionary history [21,22]. Additionally, Paraguay is key for works in northern Argentina, Uruguay, southwestern Brazil, and Bolivia. In spite of this biogeographical importance, Paraguay has been poorly explored, and in the current era of molecular genetics, the investigations that include genetic samples from Paraguay are extremely rare in herpetology. For instance, the natural history museum of Paraguay (Museo Nacional de Historia Natural del Paraguay) started its tissue collection for genetic analyses in this decade, whereas other neighbor countries began cryo tissue collections already some decades ago.
The herpetofauna from Paraguay, and specifically the squamate diversity, is still poorly known, evidenced by the fact that even in the last decade, and without the help of molecular tools, several new records for the country were made (e.g., Ophiodes fragilis, Epictia vellardi, Chironius exoletus, Lygophis paucidens, Philodryas livida, and Micrurus silviae) [23,24,25,26,27,28] and some species new to science were described (Tropidurus lagunablanca Carvalho, 2016; Tropidurus tarara Carvalho, 2016; Tropidurus teyumirim Carvalho, 2016; Ophiodes luciae Cacciali & Scott, 2015; Phalotris normanscotti Cabral & Cacciali, 2015) [29,30,31]. The incorporation of molecular genetics in taxonomy opened new pathways leading to a higher resolution in species delimitation, identifying several cryptic species. Some herpetologists in the region included genetic samples of Squamata from Paraguay (Gamble et al., 2012, Werneck et al., 2012, Morando et al., 2014, Recoder et al., 2014) although only very occasionally [32,33,34,35].
In 2015 we started a project of barcoding the reptile fauna (Squamata specifically) of Paraguay. It is important to mention here that the two genetic markers commonly used for barcoding analyses are the mtDNA genes 16S rRNA and Cytochrome Oxidase Subunit I (COI) [15,16,36,37,38] and both markers seem to work rather equally good for species identification. However, for the South American herpetofauna, 16S was more used than COI, and it is therefore better represented in GenBank for comparison. In addition, studies with 16S have shown not only good results in species recognition but also the systematic relationships among related species [37,39,40]. Thus, we decided to use sequences of the mtDNA gene 16S in our barcoding analysis.
During the project we gathered a lot of information about squamate diversity, and as a product of this work, some papers were published [41,42,43,44,45]. The use of DNA barcoding offers a starting point for recording the number of species that occur in a given region. Our results show how the use of DNA barcode data can augment and increase the accuracy of herpetological inventory surveys. Our barcoding study of the Paraguayan Squamata reveals the depth of taxonomic diversity in this country. Furthermore, our DNA barcode data represent the so far most comprehensive DNA barcode reference library for lizards and snakes of Paraguay. These reference data provide the scientific community with resources of numerous possibilities, ranging from species inventories, species identifications, taxonomic studies to wildlife trafficking.

2. Materials and Methods

2.1. Study Area

Paraguay is located in the center of South America (Figure 1) between parallels 18°18′ and 27°30′ S, and the meridians 54°19′ and 62°38′ W; with a total surface of 406,752 km2. The country is divided by the Paraguay River into two portions: the Occidental Region (commonly known as “Chaco”) with an area of 246,925 km2 (60.7% of the country), and the Eastern Region (or Oriental Region) with a surface of 159,827 km2 (39.9% of the country).
The topography of the Chaco region is a flat savanna with few small isolate hills in the center/north and east. Bad drainage creates vast flooded areas especially south and east, north-west has some dunes formation and is extremely dry. The oriental region is undulated with hills and the highest point in the center is well irrigated by tributaries of the Paraguayan and Parana rivers basins.
The climatic conditions vary in a northwestern–southeastern gradient, being more humid and cooler in the southeast. The mean annual temperature in the whole country is about 23 °C, being 24.5 °C in the western region and 22.5 °C in the eastern region. It is important to note that there are two seasons, the wet season, in which it rains frequently, coincide with the warm period from October to April and the dry season from May to September, where rain is less frequent and coincides with the coldest period.
There is a big difference with respect to the variation in temperature, given that the mean maximum is 25 °C, but the absolute maximum temperature could reach around 50 °C, especially in the northwestern portion of the country in January or February. The coldest month is July, and the absolute minimum temperature can be −6 °C in the south. Nevertheless, in Paraguay the “true” winter usually does not last longer than 16 days each year. Thus, Paraguay has a warm/hot climate during most parts of the year.

2.2. Data Collection

Even though there are some blanks in the areas sampled, the coverage of collecting sites in this study is rather vast (Figure 2). Nevertheless, there are two ecoregions in the Occidental Region (Cerrado Chaqueño and Médanos del Chaco), from where we have no samples. The methods used in the field were the traditional techniques for herpetology: Active searching at different times of the day and night, examining potential shelters (e.g., barks, logs, caves, mud, leaf litter, etc.) [46] (Figure 3). Fast moving lizards (e.g., Ameiva and Teius) were collected using compressed air rifles [47]. Additionally, some habitats, such as ant nests and swamps, were dug looking for hypogeal organisms [48], and floating vegetation was sampled using a trawl net (Figure 3). In total, 147 days of fieldwork were accounted for this project, and about 400 specimens collected.
It is important to highlight that the exotic lizard Hemidactylus mabouia is currently widely distributed in the country and now is part of the Paraguayan herpetofauna, and thus also included in the study. These genetic data may help in future studies about the colonization of the species, which in Paraguay has been recorded in the Concepción, San Pedro, Central, Alto Paraná, and Itapúa departments [49].
Reptiles that were captured alive were euthanized with a pericardial injection of a solution of embutramide, mebezonium iodide, and tetracaine hydrochloride (T-61®, Intervet International GmbH, Unterschleissheim, Germany) or Sodium Thiopental (Tiopental Sódico®, Biosano, Chile). The Secretaría del Ambiente from Paraguay (Currently “Ministerio del Ambiente y Desarrollo Sustentable”) authorized the collecting of specimens though permits SEAM [Secretaría del Ambiente] N° 004/11 and 009/2014. Exportation permits for tissues and specimens were also issued by the same authority through the permits SEAM N° 002/14, 016/2016, and 084/2016.
After euthanasia, tissue samples were taken either from the muscle of the thigh, tongue, finger clips, tail (when regenerated), or liver. Tissues were preserved in vials containing 98% non-denatured ethanol, and stored at −20 °C as soon as possible.
Hemipenes of Squamata were everted after euthanasia, with an injection of 70% ethanol after manually everting the organs. All specimens were fixed with a solution of 36% formalin and 96% ethanol in the proportion of 5:1000 (e.g., 5 mL formalin in 1 L ethanol), injected in the body cavity, thighs, and thickest part of the tail. Following fixation, the specimens were maintained in 70% ethanol.

2.3. Molecular Protocols

We used two different methods of DNA extraction. For sets containing few samples (usually eight or fewer), we used the DNeasy® Blood & Tissue Kit of Qiagen® (Hilden, Germany), whereas for sets of 96 samples we used the fiberglass plate [50]. Both methods are detailed below. The DNA was isolated from tissues whenever possible, or taken from preserved specimens that had been stored for a considerable time in 70% ethanol at room temperature in some cases.
For the DNeasy® method, we used tissue fragments of ~2 mm2. When buffers formed precipitates, they were warmed up at 56 °C before use. All reagents for this protocol are included in the kit. Tissues were digested adding 180 μL (all values are for individual samples) of ATL Buffer and 20 μL of proteinase K. Samples in that mix were incubated in a rocking platform at 56 °C for 4 to 12 h until the tissue was completely lysed.
Following digestion, 200 μL of AL lysis Buffer + 200 μL of ethanol (98%) were added. This mix was centrifuged (8000 rpm) in DNeasy® Mini spin columns, discarding all the flow-through. Then, 500 μL of AW1 washing Buffer was added and centrifuged (8000 rpm) discarding the flow-through. Finally, 500 μL of AW2 washing Buffer was added and centrifuged (14,000 rpm) discarding the flow-through. The final elution was made with 200 μL of AE Buffer, after an incubation of one minute, followed by centrifugation (8000 rpm).
For the fiberglass extraction method, we used tissue fragments of ~1 mm2. Specifications of reagents used in this protocol, are detailed in Table S1. Initially, the samples were washed with 50 μL (values per sample) of a solution of 1× Tris-Ethylenediaminetetraacetic acid (TE) Buffer to remove the remaining ethanol, for ~15 h. Following, the samples were digested with 50 μL of a solution of Vertebrate lysis Buffer and proteinase K (10:1), and incubated in a rocking platform at 56 °C for 12–24 h.
Once the samples were digested, the DNA extraction was made adding 100 μL of Binding Buffer and centrifuging at 2800 rpm. These products were transferred to a Pall® (Cortland, NY, USA) AcroPrep® filter plate, where the plate was vacuumed for 2 min. Then it was added 180 μL of Washing Buffer 1 and vacuumed again for 2 min. Posteriorly, it was added 750 μL of the Washing Buffer 2 and vacuumed for 2. Then TE Buffer was used to elute the DNA, adding 50 μL and incubating it for 2 min at 56 °C.
We amplified fragments of the mtDNA 16S gene with forward and reverse reactions using the following primers respectively: F: L2510 (5′-CGCCTGTTTAACAAAAACAT-3′) and R: H3056 (5′-CGGTCTGAACTCAGATCACGT-3′) [6]. The master mix cocktail used for amplification was of 1 μL of the DNA template, 2.5 μL of Y Buffer, 4 μL of dNTPs, 0.5 μL of TaqPol, 1 μL of MgCl2, 1 μL of the forward and reverse primers, and 14 μL of distilled water, reaching a final volume of 25 μL. Amplification reactions were performed in an Eppendorf Mastercycler® pro (Hamburg, Germany) thermocycler using the following PCR conditions: initial denaturation 2 min (94 °C)—[denaturation 35 s (94 °C)—hybridization 35 s (48.5 °C)—elongation 60 s (72 °C)] × 40—final elongation 10 min (72 °C) [51]. Sequencing was performed using a BigDye® Terminator (ThermoFisher Scientific®, Waltham, MA, USA) with the following cycling conditions: 1 min at 95 °C, 30 × [10 s at 95 °C, 10 s at 50 °C, 2 min at 60 °C], with 10 μL of reaction volume.
Additionally to our own samples, we included non-Paraguayan data from sequences downloaded from GenBank, selecting preferably those sequences associated with museum vouchers, to avoid common problems of misidentifications in that repository [52,53,54]. In most cases, we downloaded only sequences from species represented in our samples, except for Bothrops and Amphisbaena. In these two cases, we downloaded samples from all the species present in Paraguay, because of the difficulty of these taxa for morphological identification. In the case of Micrurus, there is only one sequence of a species present in Paraguay available in GenBank (JQ627286). A particular case was the only available sample of Vanzosaura rubricauda (AF420716) uploaded in the framework of a lizards’ phylogeny [55]. That specimen (MRT 05059) from Vacaria, Estado de Bahia, Brazil, actually is V. multiscutata [35]. Nevertheless, it was included in the analysis to evaluate the clustering with the genetic sample from Paraguay.
Codes of sequences downloaded from GenBank, plus accession numbers of sequences generated in this work, are available in the Table S2.

2.4. Data Analysis

Chromatograms of forward and reverse sequences were assessed, and a consensus sequence for each sample generated in SeqTrace 0.9.0 [56]. For sequences alignment, we employed MAFFT2 [57,58] through the webserver [59], which includes a special search strategy (Q-INS-i) for the secondary structure of the rRNA 16S [60]. No later manual edition was introduced. Results of MAFFT2 were visualized in MSA Viewer [61] and exported as fasta files.
The best scheme for substitution model was explored in PartitionFinder 2.1.1 [62], using linked branch lengths (supported by most of the phylogenetic programs) using a PhyML 3.0 analysis [63]. Given that the analysis is based on Squamata, which can show highly divergent clades, we estimated the relative quality of the statistical models using the Bayesian Information Criterion (BIC) [64] since it penalizes more the number of parameters in the model and then is better for a large degree of heterogeneity [65].
Given that it is not recommended to use both +I (significant proportion of invariable sites) plus +G (rate of variation among sites follows a gamma distribution) together in the same substitution model, we chose the best suggested model using +I or +G in the partition schemes, but never both together [66].
The phylogenetic hypothesis was performed under a Maximum Likelihood (ML) approach, using IQ-Tree [67] through its webserver [68], setting 10,000 non-parametric bootstrap replicates plus 10,000 replicates of Shimodaira-Hasegawa approximate likelihood ratio (SH-aLRT) [69] and 10,000 ultrafast bootstrap (UFBoot) approximation replicates [70]. We used a sequence of Sphenodon punctatus to root the tree [71], which has been proposed as the sister clade to Squamata [72]. Here, it is important to highlight that the phylogenetic hypothesis is used to sorting groups, and we do not seek for a comprehensive evolutionary reconstruction.
For visualization and edition (branch arrangement, colors, font sizes, etc.) of the tree generated through ML analysis, we used FigTree 1. 4.3 [73]. The final alignment plus the ML tree are stored in TreeBASE repository (https://treebase.org/) under the submission number 24616. To do this, we first managed the alignments and trees in nexus format and combined them in a single file (containing one alignment and the ML tree) using Mesquite 3.31 [74].

3. Results

We generated a total of 142 sequences of 64 species of Squamata from Paraguay, including one exotic species: Hemidactylus mabouia. In the Table S2, we present a list of specimens used for genetic analyses based on the field work for this project. For comparison we added 107 sequences from GenBank (Table S2). The final alignment constituted of a dataset of 249 samples of 615 bp length. Sequences are available in GenBank.
The best substitution model for the Barcoding dataset was GTR+G, according to the BIC. The sample of Sphenodon punctatus was retrieved as the sister clade to the Squamata (Figure 4). Deep nodes have low bootstrap values, meaning that the phylogenetic relationships are weakly supported. Nevertheless, the shallowest divergences have higher support values, recovering most of the genera included in the analysis as monophyletic, with the exception of Manciola (Scincidae) and the tribe Xenodontini (Colubridae).
The tribe Xenodontini (Figure S1) of the Subfamily Dipsadinae (Colubridae) contains the samples of Erythrolamprus aesculapii in a monophyletic clade, whereas E. poecilogyrus appears as paraphyletic. Erythrolamprus reginae clusters sister to the above-mentioned taxa. The genus Xenodon seems to be paraphyletic, given that two samples of Xenodon pulcher are sister to Erythrolamprus, whereas Xenodon merremi is sister to Xenodon pulcher + Erythrolamprus. Finally, in this clade a sample of Erythrolamprus sagittifer is sister to a sample of Lygophis dilepis.
Sister to Xenodintini is a clade composed of Phalotris + Philodryas (Figure S2). Both genera are monophyletic in the tree. The genera Psomophis and Dipsas are clustered together, and nested as sister to the above-mentioned snakes (Figure S3). Other genera of Dipsadinae that are rendered as monophyletic are Hydrodynastes, Helicops, and Thamnodynastes (Figure S4). A clade grouping members of the Colubrinae subfamily is composed of three genera (Chironius, Leptophis, and Palusophis) that also show monophyly (Figure S5). The Pseudoboini (Dipsadinae: Colubridae) is shown in its own clade (Figure S6) with four of the five genera used in the analysis being monophyletic, whereas the two species of the genus Phimophis appear in different positions of the gene tree.
The genus Micrurus (Elapidae) is the sister clade of the Colubridae, whereas the genus Bothrops (Viperidae) seems to be the sister to Elapidae + Colubridae (Figure S7). Located in a most basal position among snakes are the two species of Epicrates (Boidae), with Amerotyphlops (Typhlopidae) as the sister clade of the remaining snakes (Figure S8). The genus Amphisbaena (Amphisbaenidae) is monophyletic, where A. alba and A. bolivica are in their own clades, and A. mertensii shows also monophyly (Figure S9). Amphisbaena angustifrons is the sister taxon of the other Amphisbaena, with a sample of Amphisbaena sp. (PCS 314) as the most basal taxon of the clade (Figure S9). In our analysis, Amphisbaena is sister to Teiidae + Gymnophthalmidae. Gymnophthalmidae appears as a monophyletic clade, and the four genera show monophyly as well (Figure S10).
The Family Teiidae is shown as paraphyletic. The only Tupinambinae in our samples was Salvator, which clusters as sister to Gymnophthalmidae (Figure S11). Samples of Teiinae are clustered together showing monophyly, where Ameivula and Kentropyx are sister clades (Figure S11), as are Teius and Ameiva (Figure S12). The Family Scincidae is sister to the all above-mentioned clades (Figure S13). Samples of Manciola show paraphyly (Figure S13). The remaining cluster contains members of the Anguidae, Gekkonidae, Phyllodactylidae, Liolaemidae, Polychrotidae, and Tropiduridae families. The clade composed by geckos shows monophyly in the genera, but not in the families given that Phyllopezus and Homonota are currently placed in Phyllodactylidae, whereas Hemidactylus and Lygodactylus are Gekkonidae (Figure S14). Sister to the Gekkota (Gekkonidae + Phyllodactylidae) is Liolaemus, and Stenocercus is rendered as sister to the Gekkota + Liolaemus.
The last cluster, sister to Gekkota + Liolaemus + Stenocercus, is represented by Ophiodes (Anguidae), Polychrus (Polychrotidae) (Figure S15), and Tropidurus (Tropiduridae) (Figure S16). In this case, the Family Tropiduridae is polyphyletic since the other member of the family (Stenocercus) is sister clade to Gekkota + Liolaemuus.

4. Discussion

Molecular genetics, and in particular barcoding analyses, proved to be a powerful tool to generate preliminary information about the taxonomic status of problematic taxa [20,75]. We present here the most comprehensive analysis of genetic samples of Squamata from Paraguay. The results obtained here will be useful to help identify questionable specimens and in some cases also to clarify some taxonomic issues of the Squamata fauna from the central region of South America. Thus, the data generated here will have a positive impact in a larger geographic context, beyond Paraguay’s borders.
As said before, genetics alone will not yield a well-founded taxonomy. Nevertheless, molecular genetics open a path for defining operational taxonomic units (OTUs), identifying potential undescribed species and pointing to taxonomic problems, and thus have to be seen as a first informative step and a complementary evidence line in the framework of the modern integrative taxonomic approach [20,76].
Some taxonomic results of this project were already published. For instance, the samples of Colobosaura exhibit large genetic distances, and then Colobosaura kraepelini was revalidated [41]. The Tropidurus samples show monophyly in the species of the torquatus group (T. catalanensis and T. etheridgei), but indicate several uncertainties within the spinulosus group (formerly T. guarani, T. lagunablanca, T. spinulosus, T. tarara, and T. teyumirim), that resulted in the synonymization of T. guarani with T. spinulosus, and T. tarara and T. teyumirim with T. lagunablanca [43]. Regarding the Family Phyllodactylidae, there is strong evidence for the recognition of two different Homonota species in the Chaco [42,44] and a highly distinctive Phyllopezus clade, separated from populations from Cerrado and Chaco [45].
The samples of Vanzosaura rubricauda from Cerrado (field number “ALA”) show a high branch distance compared with Vanzosaura rubricauda from Chaco (GK 3801) which is even larger than the distance from V. multiscutata (Figure S10). Integrating molecular and morphological data, a new species of Vanzosaura (V. savanicola) was previously described, and Gymnodactylus multiscutatus was transferred to the genus Vanzosaura [35]. Nevertheless, their genetic tree [35] included only a single sample from Paraguay and none from Argentina. In their map, obviously two divergent populations of Vanzosaura rubriauda are recognized: One west of the Paraguay River in the Dry Chaco, and another east of the Paraguay River in the Cerrado. Keeping a conservative approach, the authors maintained V. rubricauda as a single taxonomic unit, but with our additional samples it might be possible to generate new taxonomic hypotheses.
Furthermore, we recommend further studies on Amphisbaenidae, because one of the major and latest revisions of Amphisbaenidae in the Neotropics concluded that A. mertensi and A. cunhai (not recorded in Paraguay) are the most basal lineages of the genus [77]. Our analysis showed that the most basal sample (Amphisbaena sp. PCS 314) seems to be a different species as those within the remaining clade. Additional analyses, including more samples and a detailed morphological revision, are necessary to assess the specific status of that specimen.
The weakest part of this work was the analysis of snakes. These animals are usually the harder ones to sample. Compared to the actual diversity of Colubridae, our dataset had fewer samples of this family and therefore it was not possible to draw detailed taxonomic conclusions. However, the presence of the genus Xenodon in two different clusters suggests that more taxonomic work with this group of snakes is needed. Several taxonomic modifications occurred within the Colubridae in the last decade, where the genera Lystrophis and Waglerophis were synonymized with Xenodon, based on the analysis of gene sequences 12S and 16S for the genus Lystrophis, and Cytb and bdnf for one sample of Waglerophis merremii [78]. In our analysis, we found the samples of X. pulcher (previously Lystrophis pulcher) separated from X. merremii (previously Waglerophis merremii). It is desirable to perform phylogenies in this group using more nuclear data to get more robust relationships in the deep nodes.
In the clade of the genus Thamnodynastes, the two species used in our analysis (T. chaquensis and T. hypoconia) are nested in the same node. A more specific genetic study of the genus is highly recommendable. Furthermore, the revision of the phylogenetic status of Phimophis is advisable, since here it appears polyphyletic. In a former study, two samples of Phimophis were used: P. guerini (GQ457761) and P. iglesiasi (JQ598891) and due to polyphyly the authors described the genus Rodriguesophis to include the latter species [78]. This genus is characterized by the absence of the loreal scale. Both P. guerini and P. vittatus have a loreal scale, so they cannot be assigned to Rodriguesophis. Thus, a deeper integrative (morphological and molecular) analysis is needed to understand their relationships.
The genus Micrurus is scarcely represented in GenBank, and comparisons are not possible. The only sample from GenBank is M. altirostris, which is differentiated from Paraguayan samples by a rather long branch distance. There is a polytomy with three samples (PCS 310, 334, and 337), and the only identified specimen is Micrurus pyrrhocryptus (PCS 310) from Pantanal (northern part of Paraguay). The other samples are from Concepción at the other side of the river, and with a body color different from the pattern of M. pyrrhocryptus. It is important to highlight here that some of the specimens that we have of Micrurus, were decapitated (killed by rural farmers) and therefore, without genetic samples for comparison (in GenBank) and without cephalic data (which contains important diagnostic characters) [79], its specific taxonomic allocation becomes difficult.
Some of our Paraguayan samples of Bothrops of the neuwiedi complex, from distant parts of Paraguay, are clustered with a sample of B. diporus from GenBank (Samples PCS 302, PCS 318, PCS 331, PCS 504, Figure S7). A thorough analysis of this complex of Bothrops is needed to understand the true diversity in the group.
Regarding the Scincidae, a sample from GenBank of Manciola guaporicola (KX364960) from Mbaracayú Reserve (Paraguay) appears out of the clade of the remaining M. guaporicola from Paraguay and Brazil. This is also a topic that should be further investigated.
Regarding conservation, one of the major problems in Paraguay for several years was habitat loss due to extensive soybean crops in the eastern part of Paraguay [80]. Nevertheless, habitat fragmentation is currently also affecting the landscapes of the Occidental Region of Paraguay [81,82]. Thus, currently, the protected areas are the best strategy for conservation of biodiversity in Paraguay, although many conservation units face legal problems (e.g., lack of official measurements, management plans, forest guards, infrastructure, etc.) and then the long-term maintenance of their biodiversity is not guaranteed [80]. There are some reptile species absent from protected areas in Paraguay; therefore, monitoring and conservation efforts should be intensified for these taxa [83]. One species recently revalidated is Colobosaura kraepelini. This lizard is known only from the holotype from the locality of Puerto Max (San Pedro Department), the neotype from Altos and an additional specimen from San Bernardino, both localities in Cordillera Department [41]. There are no protected areas in the Cordillera Department, but there are some in the northern portion of Central Department (border with Cordillera), located less than 10 km from the known localities of C. kraepelini. The presence of this species in a conservation unit should be confirmed, but is possible that it is protected by “Monumento Natural Cerro Chororí” and “Monumento Natural Cerro Kõi”. It is important to note that the conservation unit closer to the distribution of C. kraepelini is the “Parque Nacional Lago Ypacaraí”, although only the lagoon is protected and not the surroundings. The species Homonota septentrionalis was described from the driest part of Paraguay (northwestern Chaco) and is abundant in the “Parque Nacional Teniente Enciso” [44]. The four species of Tropidurus found in Paraguay [43] are well represented in several protected areas. The last herpetofaunal conservation assessment was published in 2009 [84], and thus a new conservation assessment of Paraguayan reptile fauna including the new taxa, is necessary to provide a sound basis for conservation planning for those species that require special attention.
Finally, given that the sequenced specimens are yet a small portion of the actual diversity of Paraguay (Figure 5), it will be of the utmost importance to continue and expand these studies that will further improve our taxonomic knowledge.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-2818/11/9/152/s1, Table S1: Reagents used for molecular procedures; Table S2: Specimens used in the barcoding analysis; Figures S1–S16: Details of sub-trees sections from the general barcoding analysis.

Author Contributions

Conceptualization, P.C., E.B., and G.K.; data curation, P.C., E.B., and G.K.; methodology, P.C.; formal analysis, P.C.; supervision, G.K.; writing—original draft preparation, P.C.; writing—review and editing, G.K. and E.B.; funding acquisition, P.C. and E.B.

Funding

This research was developed in the framework of the Ph.D. thesis project “Lizards of Paraguay: an integrative approach to solve taxonomic problems in central South America” presented by P.C. in the Faculty of Biosciences of the Johann Wolfgang Goethe-Universität, which had a strong financial support from the Deutscher Akademischer Austauschdienst (DAAD). Additionally, P.C. received a subsidy from the Consejo Nacional de Ciencia y Tecnología (CONACYT) through the Programa Nacional de Incentivo a los Investigadores (PRONII), field equipment provided by Idea Wild, and a grant to perform fieldwork activities in the Mbaracayú Reserve by funded by the GNB Bank through the Fundación Moisés Bertoni (FMB). Additionally, we received financial support from PRESIDENT ENERGY to pay the APC of this work.

Acknowledgments

This was a major project that had the help of many people. For this, we owe gratitude to Norman Scott, Luciano Avila, Mariana Morando, Tony Gamble, Sebastian Lotzkat, Martin Jansen, Arne Schultze, Christian Printzen, Raúl Maneyro, and Abel Batista, who provided advice and technical support at different stages of the project. This project would not have been achieved without the valuable help of many colleagues and friends who offered their help during field work. Hence, thank you so much Jorge Ayala, Frederick Bauer, Aníbal Bogado, Enrique Bragayrac, Diego Bueno, Emilio Buongermini, “Huguitus” Cabral, Paulo Campos Filho, Gloria “Lolex” Céspedes, Jorge A. Céspedez, Julia Coda, Sixto Fernández, Andrea Ferreira, Marcela Ferreira, Celeste “Tita” Gauto, Kevin Guest, Hugo del Castillo, Stefan Harrison (you rock man!), Monica Kozykariski, Victoria Kuntz, Arne Lestheruis, Pamela Marchi, Nicolás Martínez, José Méndez, Martha Motte, Cristian F. Pérez, Pastor Pérez, Rachel Pitts, Sergio Ríos, Lía Romero, Mirtha Ruiz Díaz, Humberto “Ka’umber” Sánchez, Rebecca Sheehan, Nelson Silva, María E. Tedesco, Dulcy Vázquez, Thomas and Sabine Vinke, Akira Yoshikawa, and Víctor Zaracho. Furthermore, we have to acknowledge the hospitality of persons (Rubén Ávila and Nilda Torres de Ávila, Massimo and Angela Coda, Rosario Gabaglio, Ana María Macedo, and Goli Stroessner) and organizations (Fundación Moisés Bertoni, Guyra Paraguay, and Para La Tierra) that provided accommodations during field wok. We also want to express our gratitude to some friends and colleagues that provided tissue samples, literature, or valuable information. They are Luciano J. Avila, Francisco Brusquetti, Arley Camargo, Santiago Carreira, André de Carvalho, Mariana Morando, and Paul Smith. Additionally, we are grateful to the staff (especially Heike Kappes) of the Grunelius-Möllgaard Laboratory (Senckenberg Forschungsinstitut und Naturmuseum Frankfurt - SMF), and to Linda Mogk (SMF) for lab support. P.C. also wants to thank to Tachi, Dante and Rafa for patience and support. Finally, to the CONACYT for providing important tools for scientific research in the country.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Pullin, A.; Sutherland, W.; Gardner, T.; Kapos, V.; Fa, J.E. Conservation priorities: Identifying need, taking action and evaluating success. In Key Topics in Conservation Biology 2; Macdonald, D.W., Willis, K.J., Eds.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2013; pp. 3–22. [Google Scholar] [CrossRef]
  2. Baumgartner, J.; Esselman, R.; Salzer, D.; Young, J. Conservation Action Planning Handbook; The Nature Conservancy: Arlington, VA, USA, 2007. [Google Scholar]
  3. Ripple, W.J.; Wolf, C.; Newsome, T.M.; Galetti, M.; Alamgir, M.; Crist, E.; Mahmoud, M.I.; Laurance, W.F. World scientists’ warning to humanity: A second warning. Bioscience 2017, 67, 1026–1028. [Google Scholar] [CrossRef]
  4. Young, B.E.; Stuart, S.N.; Chanson, J.S.; Cox, N.A.; Boucher, T.M. Disappearing Jewels: The Status of New World Amphibians; NatureServe: Arlington, VA, USA, 2004. [Google Scholar]
  5. Johnson, C.N.; Balmford, A.; Brook, B.W.; Buettel, J.C.; Galetti, M.; Guangchun, L.; Wilmshurst, J.M. Biodiversity losses and conservation responses in the Anthropocene. Science 2017, 356, 270–275. [Google Scholar] [CrossRef] [PubMed]
  6. Palumbi, S.R.; Martin, A.; Romano, S.; Mcmillan, W.O.; Stice, L.; Grabowski, G. The Simple Fool’s Guide to PCR; University of Hawaii Press: Honolulu, HI, USA, 1991. [Google Scholar]
  7. Rudnick, J.A.; Katzner, T.E.; Bragin, E.A.; DeWoody, J.A. Species identification of 548 birds through genetic analysis of naturally shed feathers. Mol. Ecol. Notes 2007, 7, 757–762. [Google Scholar] [CrossRef]
  8. Yang, L.; Tan, Z.; Wang, D.; Xue, L.; Guan, M.X.; Huang, T.; Li, R. Species identification through mitochondrial rRNA genetic analysis. Sci. Rep. 2014, 4, 1–11. [Google Scholar] [CrossRef] [PubMed]
  9. Sites, J.; Marshall, J.C. Delimiting species: A renaissance issue in systematic biology. Trends. Ecol. Evol. 2003, 18, 462–470. [Google Scholar] [CrossRef]
  10. Pinzón, J.H.; LaJaunesse, T.C. Species delimitation of common reef corals in the genus Pocillopora using nucleotide sequence phylogenies, population genetics and symbiosis ecology. Mol. Ecol. 2011, 20, 311–325. [Google Scholar] [CrossRef] [PubMed]
  11. Khodami, S.; Martínez Arbizu, P.; Stöhr, S.; Laakmann, S. Molecular species delimitation of Icelandic brittle stars (Ophiuroidea). Pol. Polar Res. 2014, 35, 243–260. [Google Scholar] [CrossRef]
  12. Leliaert, F.; Verbruggenc, H.; Vanormelingend, P.; Steena, F.; López-Bautistab, J.M.; Zuccarelloe, G.C.; De Clercka, O. DNA-based species delimitation in algae. Eur. J. Phycol. 2014, 49, 179–196. [Google Scholar] [CrossRef] [Green Version]
  13. Barcaccia, G.; Lucchini, M.; Cassandro, M. DNA Barcoding as a molecular tool to track down mislabeling and food piracy. Diversity 2016, 8, 2. [Google Scholar] [CrossRef]
  14. Lane, N. On the origin of bar codes. Nature 2009, 462, 272–274. [Google Scholar] [CrossRef]
  15. Hebert, P.D.N.; Cywinska, A.; Ball, S.; de Waard, J.R. Biological identifications through DNA barcodes. Proc. R. Soc. B 2003, 270, 313–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Murphy, R.W.; Crawford, A.J.; Bauer, A.M.; Che, J.; Donnellan, S.C.; Fritz, U.; Haddad, C.F.B.; Nagy, Z.T.; Poyarkov, N.A.; Vences, M.; et al. Cold Code: The global initiative to DNA barcode amphibians and nonavian reptiles. Mol. Ecol. Resour. 2013, 13, 161–167. [Google Scholar] [CrossRef]
  17. Perl, R.G.B.; Nagy, Z.T.; Sonet, G.; Glaw, F.; Wollenberg, K.C.; Vences, M. DNA barcoding Madagascar’s amphibian fauna. Amphib.-Reptil. 2014, 35, 197–206. [Google Scholar] [CrossRef]
  18. Shen, Y.Y.; Chen, X.; Murphy, R.W. Assessing DNA barcoding as a tool for species identification and data quality control. PLoS ONE 2013, 8, e57125. [Google Scholar] [CrossRef] [PubMed]
  19. Meiri, S.; Mace, G.M. New taxonomy and the origin of species. PLoS Biol. 2007, 5, e194. [Google Scholar] [CrossRef]
  20. Hajibabaei, M.; Singer, G.A.C.; Hebert, P.D.N.; Hickey, D.A. DNA barcoding: How it complements taxonomy, molecular phylogenetics and population genetics. Trends Genet. 2007, 23, 167–172. [Google Scholar] [CrossRef] [PubMed]
  21. Spichiger, R.; Palese, R.; Chautems, A.; Ramella, L. Origin, affinities and diversity hot spots of the Paraguayan dendrofloras. Candollea 1995, 50, 515–537. [Google Scholar]
  22. Morrone, J.J. Biogeografía de América Latina y el Caribe; Manuales & Tesis; SEA/UNESCO/CYTED: Zaragoza, España, 2001. [Google Scholar]
  23. Cacciali, P.; Scott, N. Revisión del Género Ophiodes de Paraguay (Squamata: Anguidae). Bol. Soc. Zool. Urug. 2012, 21, 1–8. [Google Scholar]
  24. Cabral, H.; Netto, F. Epictia vellardi. Geographic distribution. Herp. Rev. 2016, 47, 83. [Google Scholar]
  25. Cacciali, P.; Cabral, H. The genus Chironius (Serpentes, Colubridae) in Paraguay: Composition, distribution, and morphology. Bas. App. Herp. 2015, 29, 51–60. [Google Scholar] [CrossRef]
  26. Cacciali, P.; Smith, P.; Källberg, A.; Pheasey, H.; Atkinson, K. Reptilia, Squamata, Serpentes, Lygophis paucidens Hoge, 1952: First records for Paraguay. Check List 2013, 9, 131–132. [Google Scholar] [CrossRef]
  27. Smith, P.; Cacciali, P.; Scott, N.; del Castillo, H.; Pheasey, H.; Atkinson, K. First record of the globally-threatened Cerrado endemic snake Philodryas livida (Amaral, 1923) (Serpentes, Dipsadidae) from Paraguay, and the importance of the Reserva Natural Laguna Blanca to its conservation. Cuad. Herpetol. 2014, 28, 169–171. [Google Scholar]
  28. Cacciali, P.; Espínola, D.; Centrón Viñales, S.; Gauto Espínola, I.; Cabral, H. Squamata, Serpentes, Micrurus silviae Di-Bernardo, Borges-Martins and Silva, 2007: Presence confirmation in Paraguay. Check List 2011, 7, 809–810. [Google Scholar] [CrossRef]
  29. Carvalho, A.L.G. Three new species of the Tropidurus spinulosus group (Squamata: Tropiduridae) from Eastern Paraguay. Am. Mus. Nov. 2016, 3853, 1–44. [Google Scholar] [CrossRef]
  30. Cacciali, P.; Scott, N. Key to the Ophiodes (Squamata: Sauria: Diploglossidae) of Paraguay with the description of a new species. Zootaxa 2015, 3980, 42–50. [Google Scholar] [CrossRef]
  31. Cabral, H.; Cacciali, P. A new species of Phalotris (Serpentes: Dipsadidae) from the Paraguayan Chaco. Herpetologica 2015, 71, 72–77. [Google Scholar] [CrossRef]
  32. Gamble, T.; Colli, G.R.; Rodrigues, M.T.; Werneck, F.P.; Simons, A.M. Phylogeny and cryptic diversity in geckos (Phyllopezus; Phyllodactylidae; Gekkota) from South America’s open biomes. Mol. Phyl. Evol. 2012, 62, 943–953. [Google Scholar] [CrossRef]
  33. Werneck, F.P.; Gamble, T.; Colli, G.R.; Rodrigues, M.T.; Sites, J. Deep diversification and long-term persistence in the South American “Dry Diagonal”: Integrating continent-wide phylogeography and distribution modeling of geckos. Evolution 2012, 66, 3014–3034. [Google Scholar] [CrossRef]
  34. Morando, M.; Medina, C.D.; Ávila, L.J.; Pérez, C.H.F.; Buxton, A.; Sites, J.W. Molecular phylogeny of the New World gecko genus Homonota (Squamata: Phyllodactylidae). Zool. Scr. 2014, 43, 249–260. [Google Scholar] [CrossRef]
  35. Recoder, R.S.; Werneck, F.; Teixeira, M.; Colli, G.R.; Sites, J.W.; Rodrigues, M.T. Geographic variation and systematic review of the lizard genus Vanzosaura (Squamata, Gymnophthalmidae), with the description of a new species. Zool. J. Lin. Soc. 2014, 171, 206–225. [Google Scholar] [CrossRef]
  36. Sacchi, C.T.; Whitney, A.M.; Mayer, L.W.; Morey, R.; Steigerwalt, A.; Boras, A.; Weyant, R.S.; Popovic, T. Sequencing of 16S rRNA gene: A rapid tool for identification of Bacillus anthracis. Emerg. Infect. Dis. 2002, 8, 1117–1123. [Google Scholar] [CrossRef]
  37. Jansen, M.; Schultze, A. Molecular, morphology and bioacoustic data suggest Bolivian distribution of a large species of the Leptodactylus pentadactylus group (Amphibia: Anura: Leptodactylidae). Zootaxa 2012, 3307, 35–47. [Google Scholar] [CrossRef]
  38. Scherz, M.D.; Vences, M.; Borrell, J.; Ball, L.; Nomenjanahary, D.H.; Parker, D.; Rakotondratsima, M.; Razafimandimby, E.; Starnes, T.; Rabearivony, J.; et al. A new frog species of the subgenus Asperomantis (Anura, Mantellidae, Gephyromantis) from the Bealanana District of northern Madagascar. Zoosyst. Evol. 2017, 93, 451–466. [Google Scholar] [CrossRef]
  39. Batista, A.; Hertz, A.; Köhler, G.; Mebert, K.; Veselý, M. Morphological variation and phylogeography of frogs related to Pristimantis caryophyllaceus (Anura: Terra-rana: Craugastoridae) in Panama. Salamandra 2014, 50, 155–171. [Google Scholar]
  40. Köhler, G.; Townsend, J.H.; Petersen, C.B. A taxonomic revision of the Norops tropidonotus complex (Squamata, Dactyloidae), with the resurrection of N. spilorhipis (Álvarez del Toro and Smith, 1956) and the description of two new species. Mesoam. Herpetol. 2016, 3, 8–41. [Google Scholar]
  41. Cacciali, P.; Martínez, N.; Köhler, G. Revision of the phylogeny and chorology of the tribe Iphisini with the revalidation of Colobosaura kraepelini Werner, 1910 (Reptilia, Squamata, Gymnophthalmidae). ZooKeys 2017, 669, 89–105. [Google Scholar] [CrossRef]
  42. Cacciali, P.; Morando, M.; Medina, C.D.; Köhler, G.; Motte, M.; Avila, L.J. Taxonomic analysis of Paraguayan samples of Homonota fasciata Duméril & Bibron (1836) with the revalidation of Homonota horrida Burmeister (1861) (Reptilia: Squamata: Phyllodactylidae) and the description of a new species. PeerJ 2017, 5, e3523. [Google Scholar] [CrossRef]
  43. Cacciali, P.; Köhler, G. Diversity of Tropidurus (Squamata: Tropiduridae) in Paraguay—An integrative taxonomic approach based on morphological and molecular genetic evidence. Zootaxa 2018, 4375, 511–536. [Google Scholar] [CrossRef]
  44. Cacciali, P.; Morando, M.; Avila, L.J.; Köhler, G. Description of a new species of Homonota (Reptilia, Squamata, Phyllodactylidae) from the central region of northern Paraguay. Zoosyst. Evol. 2018, 94, 147–161. [Google Scholar] [CrossRef]
  45. Cacciali, P.; Lotzkat, S.; Gamble, T.; Köhler, G. Cryptic diversity in the Neotropical gecko genus Phyllopezus Peters, 1878 (Reptilia: Squamata: Phyllodactylidae): A new species from Paraguay. Int. J. Zool 2018. [Google Scholar] [CrossRef]
  46. Cacciali, P. Colecta y Preparación de Anfibios y Reptiles: Manual para Colecta Científica; Editorial Académica Española: Saarbrücken, Germany, 2013. [Google Scholar]
  47. Scrocchi, G.; Kretzschmar, S. Guía de Métodos de Captura y Preparación de Anfibios y Reptiles para Estudios Científicos y Manejo de Colecciones Herpetológicas; Miscelánea N° 102; Fundación Miguel Lillo: San Miguel de Tucumán, Argentina, 1996. [Google Scholar]
  48. Simmons, J.E. Herpetological Collecting and Collections Management; Society for the Study of Amphibians and Reptiles: Salt Lake, UT, USA, 2002. [Google Scholar]
  49. Cacciali, P.; Scott, N.; Aquino, A.L.; Fitzgerald, L.A.; Smith, P. The Reptiles of Paraguay: Literature, distribution, and an annotated taxonomic checklist. Spec. Publ. Mus. Southwest. Biol. 2016, 11, 1–373. [Google Scholar]
  50. Ivanova, N.V.; Dewaard, J.R.; Hebert, P.D. An inexpensive, automation-friendly protocol for recovering high-quality DNA. Mol. Ecol. Notes 2006, 6, 998–1002. [Google Scholar] [CrossRef]
  51. Lotzkat, S.; Hertz, A.; Bienentreu, J.F.; Köhler, G. Distribution and variation of the giant alpha anoles (Squamata: Dactyloidae) of the genus Dactyloa in the highlands of western Panama, with the description of a new species formerly referred to as D. microtus. Zootaxa 2013, 3626, 1–54. [Google Scholar] [CrossRef]
  52. Vilgalys, R. Taxonomic misidentification in public DNA databases. New Phytol. 2003, 160, 4–5. [Google Scholar] [CrossRef]
  53. Zhang, G. Specimens versus sequences. Science 2009, 323, 1672. [Google Scholar] [CrossRef]
  54. Park, K.S.; Ki, C.S.; Kang, C.I.; Kim, Y.J.; Chung, D.R.; Peck, K.R.; Song, J.H.; Lee, N.Y. Evaluation of the GenBank, EzTaxon, and BIBI services for molecular identification of clinical blood culture isolates that were unidentifiable or misidentified by conventional methods. J. Clin. Microbiol. 2012, 50, 1792–1795. [Google Scholar] [CrossRef]
  55. Pellegrino, K.C.M.; Rodrigues, M.T.; Yonenaga-Yassuda, Y.; Sites, J.W. A molecular perspective on the evolution of microteiid lizards (Squamata, Gymnophthalmidae), and a new classification for the family. Biol. J. Lin. Soc. 2001, 74, 315–338. [Google Scholar] [CrossRef] [Green Version]
  56. Stucky, B.J. SeqTrace: A graphical tool for rapidly processing DNA sequencing chromatograms. J. Biomol. Tech. 2012, 23, 90–93. [Google Scholar] [CrossRef]
  57. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef]
  58. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef]
  59. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 2017, bbx108. [Google Scholar] [CrossRef]
  60. Katoh, K.; Toh, H. Recent developments in the MAFFT multiple sequence alignment program. Brief. Bioinform. 2008, 9, 286–298. [Google Scholar] [CrossRef] [Green Version]
  61. Yachdav, G.; Wilzbach, S.; Rauscher, B.; Sheridan, R.; Sillitoe, I.; Procter, J.; Lewis, S.E.; Rost, B.; Goldberg, T. MSAViewer: Interactive JavaScript visualization of multiple sequence alignments. Bioinformatics 2016, 32, 3501–3503. [Google Scholar] [CrossRef]
  62. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. PartitionFinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 2016, 34, 772–773. [Google Scholar] [CrossRef]
  63. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef]
  64. Schwarz, G.E. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  65. Brewer, M.J.; Butler, A.; Cooksley, S. The relative performance of AIC, AICc and BIC in the presence of unobserved heterogeneity. Met. Ecol. Evol. 2016, 7, 679–692. [Google Scholar] [CrossRef]
  66. Yang, Z.; Landry, J.F.; Hebert, P.D.N. A DNA Barcode Library for North American Pyraustinae (Lepidoptera: Pyraloidea: Crambidae). PLoS ONE 2016, 11, e0161449. [Google Scholar] [CrossRef]
  67. Nguyen, L.T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating Maximum Likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  68. Trifinopoulos, J.; Nguyen, L.T.; von Haeseler, A.; Minh, B.Q. W-IQ-TREE: A fast online phylogenetic tool for Maximum Likelihood analysis. Nucleic Acid Res. 2016, 44, W232–W235. [Google Scholar] [CrossRef]
  69. Anisimova, M.; Gil, M.; Dufayard, J.F.; Dessimoz, C.; Gascuel, O. Survey of branch support methods demonstrates accuracy, power, and robustness of fast Likelihood-based approximation schemes. Syst. Biol. 2011, 60, 685–699. [Google Scholar] [CrossRef]
  70. Minh, B.Q.; Thi Nguyen, M.A.; von Haeseler, A. Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 2013, 30, 1188–1195. [Google Scholar] [CrossRef]
  71. Miller, H.C. Cloacal and buccal swabs are a reliable source of DNA for microsatellite genotyping of reptiles. Conserv. Genet. 2006, 7, 1001–1003. [Google Scholar] [CrossRef]
  72. Pyron, R.A.; Burbrink, F.T.; Wiens, J.J. A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. BMC Evol. Biol. 2013, 13, 93. [Google Scholar] [CrossRef]
  73. Bogaardt, C.; Carvalho, L.; Hill, V.; O’Toole, A.; Rambaut, A. FigTree, Version 1.4.3, Molecular Evolution, Phylogenetics and Epidemiology, United Kingdom. 2018. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 14 November 2018).
  74. Madison, W.P.; Madison, D.R. Mesquite: A Modular System for Evolutionary Analysis, Version 3.2, Mesquite Project, United States. 2017. Available online: http://mesquiteproject.org (accessed on 15 June 2018).
  75. Chapple, D.G.; Ritchie, P.A. A retrospective approach to testing DNA barcoding method. PLoS ONE 2013, 8, e77882. [Google Scholar] [CrossRef]
  76. Padial, J.M.; Miralles, A.; De la Riva, I.; Vences, M. The integrative future of taxonomy. Front. Zool. 2010, 7, 16. [Google Scholar] [CrossRef]
  77. Mott, T.; Vieites, D.R. Molecular phylogenetics reveals extreme morphological homoplasy in Brazilian worm lizards challenging current taxonomy. Mol. Phylogenet. Evol. 2009, 51, 190–200. [Google Scholar] [CrossRef]
  78. Grazziotin, F.G.; Zaher, H.; Murphy, R.W.; Scrocchi, G.; Benavides, M.A.; Zhang, Y.P.; Bonatto, S.L. Molecular phylogeny of the New World Dipsadidae (Serpentes: Colubroidea): A reappraisal. Cladistics 2012, 28, 437–459. [Google Scholar] [CrossRef]
  79. Da Silva, N.J.; Sites, J.W. Revision of the Micrurus frontalis complex (Serpentes: Elapidae). Herpetol. Monogr. 1999, 13, 142–194. [Google Scholar] [CrossRef]
  80. Cartes, J.L. Brief history of conservation in the Interior Atlantic Forest. In The Atlantic Forest of South America; Galindo-Leal, C., Gusmão Câmara, I., Eds.; Island Press: London, UK, 2003; pp. 269–287. [Google Scholar]
  81. Huang, C.; Kim, S.; Song, K.; Townshend, J.; Davis, P.; Altstatt, A.; Rodas, O.; Yanosky, A.A.; Clay, R.; Tucker, C.J.; et al. Assessment of Paraguay’s forest cover change using Landsat observations. Glob. Planet. Chang. 2009, 67, 1–12. [Google Scholar] [CrossRef]
  82. Yanosky, A.A. Paraguay’s challenge of conserving natural habitats and biodiversity with global markets demanding for products. In Conservation Biology: Voices from the Tropics; Sodhi, N.S., Gibson, L., Raven, P.H., Eds.; Wiley-Blackwell: Oxford, UK, 2013; pp. 113–119. [Google Scholar]
  83. Cacciali, P.; Cabral, H.; Yanosky, A.A. Conservation implications of protected areas’ coverage for Paraguay’s reptiles. Parks 2015, 21, 101–119. [Google Scholar] [CrossRef]
  84. Motte, M.; Núñez, K.; Cacciali, P.; Brusquetti, F.; Scott, N.; Aquino, A.L. Categorización del estado de conservación de los anfibios y reptiles de Paraguay. Cuad. Herpetol. 2009, 23, 5–18. [Google Scholar]
Figure 1. Administrative divisions of Paraguay, showing the two regions known as “Occidental Region” or Chaco (light brown) and “Oriental Region” (light green), divided by the Paraguay River (highlighted in blue in the map).
Figure 1. Administrative divisions of Paraguay, showing the two regions known as “Occidental Region” or Chaco (light brown) and “Oriental Region” (light green), divided by the Paraguay River (highlighted in blue in the map).
Diversity 11 00152 g001
Figure 2. Ecoregions of Paraguay, showing the collecting sites (brown dots) for this project.
Figure 2. Ecoregions of Paraguay, showing the collecting sites (brown dots) for this project.
Diversity 11 00152 g002
Figure 3. Sampling methods during fieldwork included diverse techniques to search in different environments.
Figure 3. Sampling methods during fieldwork included diverse techniques to search in different environments.
Diversity 11 00152 g003
Figure 4. General view of the barcoding tree performed with a Maximum Likelihood approach. Clades are collapsed (in genera or tribes) for a better visualization. Numbers on nodes represent Bootstrap values. The scale bar represents branch length (substitutions/site).
Figure 4. General view of the barcoding tree performed with a Maximum Likelihood approach. Clades are collapsed (in genera or tribes) for a better visualization. Numbers on nodes represent Bootstrap values. The scale bar represents branch length (substitutions/site).
Diversity 11 00152 g004
Figure 5. Comparison between known diversity of Paraguayan Squamata (blue bars) vs. diversity of sampled taxa for this study (red bars).
Figure 5. Comparison between known diversity of Paraguayan Squamata (blue bars) vs. diversity of sampled taxa for this study (red bars).
Diversity 11 00152 g005

Share and Cite

MDPI and ACS Style

Cacciali, P.; Buongermini, E.; Köhler, G. Barcoding Analysis of Paraguayan Squamata. Diversity 2019, 11, 152. https://0-doi-org.brum.beds.ac.uk/10.3390/d11090152

AMA Style

Cacciali P, Buongermini E, Köhler G. Barcoding Analysis of Paraguayan Squamata. Diversity. 2019; 11(9):152. https://0-doi-org.brum.beds.ac.uk/10.3390/d11090152

Chicago/Turabian Style

Cacciali, Pier, Emilio Buongermini, and Gunther Köhler. 2019. "Barcoding Analysis of Paraguayan Squamata" Diversity 11, no. 9: 152. https://0-doi-org.brum.beds.ac.uk/10.3390/d11090152

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop