Research Article |
Corresponding author: Riccardo Castiglia ( castiglia@uniroma1.it ) Academic editor: Larissa Kupriyanova
© 2020 Riccardo Castiglia, Oscar Alberto Flores-Villela, Alexandra M. R. Bezerra, Ekaterina Gornung, Flavia Annesi, Luis Antonio Muñoz-Alonso, Emanuela Solano.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Castiglia R, Flores-Villela OA, Bezerra AMR, Gornung E, Annesi F, Muñoz-Alonso LA, Solano E (2020) Detection of cryptic diversity in lizards (Squamata) from two Biosphere Reserves in Mesoamerica. Comparative Cytogenetics 14(4): 613-638. https://doi.org/10.3897/CompCytogen.v14.i4.57765
|
A combined approach based on karyology and DNA taxonomy allowed us to characterize the taxonomic peculiarities in 10 Mesoamerican lizard species, belonging to six genera and five families, inhabiting two Biosphere Reserve in Chiapas, Mexico: La Sepultura Biosphere Reserve, and Montes Azules Biosphere. The karyotypes of four species, Phyllodactylus sp. 3 (P. tuberculosus species group) (2n = 38), Holcosus festivus (Lichtenstein et von Martens, 1856) (2n = 50), Anolis lemurinus Cope, 1861 (2n = 40), and A. uniformis Cope, 1885 (2n = 29–30) are described for the first time, the last one showing a particular X1X1X2X2/X1X2Y condition. In Aspidoscelis deppii (Wiegmann, 1834) (2n = 50) and Anolis capito Peters, 1863 (2n = 42), we found a different karyotype from the ones previously reported for these species. Moreover, in A. capito, the cytogenetic observation is concurrent with a considerable genetic divergence (9%) at the studied mtDNA marker (MT-ND2), which is indicative of a putative new cryptic species. The skink Scincella cherriei (Cope, 1893), showed high values of genetic divergence (5.2% at 16S gene) between the specimens from Montes Azules and those from Costa Rica and Nicaragua, comparable to the values typical of sister species in skinks. A lower level of genetic divergence, compatible with an intraspecific phylogeographic structure, has been identified in Lepidophyma flavimaculatum Duméril, 1851. These new data identify taxa that urgently require more in-depth taxonomic studies especially in these areas where habitat alteration is proceeding at an alarming rate.
Cytotaxomy, DNA, herpetofauna, taxonomy
The Mesoamerican biota, with its number of endemics in different groups of taxa is one of the most diverse and interesting on the planet (for revision see
Saurians are one of the most representative group in terms of karyotypic diversification among reptiles (
Our previous studies aimed to genetically characterize the lizard community of a tropical dry forest in the Chamela-Cuixmala Biosphere Reserve (Jalisco state, Mexico) by means of DNA and chromosome analysis (
This study aims to extend the genetic characterization of lizard species in two additional Biosphere Reserves in Mesoamerica: La Sepultura and Montes Azules Biosphere Reserves, both in Chiapas state, Mexico. From a biotic perspective, Chiapas is an area of transition between the herpetofauna of Mexico and that of Central America, along with the one of the Yucatan Peninsula (
In the present study, karyotypes of the sampled species have been characterized. Then, in conjunction with karyotype data, mtDNA genes for different species, sequenced here and available from GeneBank, were used as molecular markers to identify new putative cryptic species and/or new evolutionarily significant units (ESU) (
Lizard specimens here analyzed were sampled in two localities: La Sepultura Biosphere Reserve, during September 2009, and Montes Azules Biosphere Reserve during 2012, Chiapas state, Mexico (Fig.
Map showing the collection localities of specimens used in this study, in La Sepultura and Montes Azules Biosphere Reserves, Chiapas state, Mexico.
The specimens were captured by hand in active searching in random walks along the surveyed localities. Details on voucher numbers, genes sequenced, chromosome complements and sampled localities, for each species are shown in Table
For karyotyping, specimens were injected with a 1:1000 solution of Velbe (Lilly) for one hour. The femurs, vertebral column, and testes were removed, crushed and left in hypotonic solution (0.075 M KCl) for 40 minutes at room temperature. Cells were collected by centrifugation and fixed with a methanol-acetic acid solution (3:1). Metaphase plates were prepared by standard air-drying method and slides were stained with Giemsa (pH = 7). Metaphases images were captured with a Photometrics Sensys 1600 digital camera (Roper Scientific Photometrics, Tucson, AZ). For each species, we identified the diploid number (2n), the number of macro- and microchromosomes, and the morphology of macrochromosomes. In some species, it was also possible to assess the morphology of the largest microchromosomes.
For molecular analyses, tissues were extracted from liver and body muscle, and preserved in 100% ethanol. A fragment of the mtDNA genome was sequenced for each species, and the sequenced genes were either cytochrome b (MT-CYB), NADH-ubiquinone oxidoreductase core subunit 2 (MT-ND2) or mitochondrially encoded 16S rRNA (16S) (Table
The QIAmp tissue extraction kit (Qiagen) was used for DNA extraction. We used the universal primers L14841 and H15149 (
Molecular identification of the specimen was performed with the BLAST algorithm (https://blast.ncbi.nlm.nih.gov/Blast.cgi) using, the newly obtained sequences and searching for highly similar sequences (Mega BLAST) on the entire nucleotide collection database. When sequence identity was below 98% the sequences were aligned with the sequences from the same species and/or same genus downloaded from GenBank. Phylogenetic relationships were evaluated with Bayesian inference (BI) and the BI tree was built with the software MrBayes v3.2.1 (
For some species a TCS Parsimony Network (
We obtained karyological and molecular data for 10 species (Fig.
Some lizard species analyzed in present study (Photos by Riccardo Castiglia) A Anolis capito B Anolis lemurinus C Holcosus festivus D Lepidophyma flavimaculatum E Anolis uniformis.
Details of gene sequenced, chromosome complement and sampled localities, for each species studied in this work. In voucher numbers, OFV indicated those specimens held in the Museum of Zoology “Alfonso L. Herrera”, Mexico City; every other specimen is held in the Museum of Comparative Anatomy of Vertebrates “Battista Grassi” of the Rome University “La Sapienza”, Rome, Italy.
Taxon | Voucher | Gene sequenced | GenBank accession numbers | Karyotype | Locality |
---|---|---|---|---|---|
Squamata | |||||
Scincidae | |||||
Scincella assata | RCMX 85 | 16S | – | 2n = 28 (7M + 14 m) | La Sepultura Biosphere Reserve |
RCMX 86 | – | ||||
RCMX 92 | MW265933 | ||||
Scincella cherriei | RCMX 219 (OFV 1197) | 16S | MW265931 | 2n = 30 (7M + 16 m) | Montes Azules Biosphere Reserve |
RCMX 235 | MW265932 | ||||
Phyllodactylidae | |||||
Phyllodactylus sp.3 | RCMX 67 RCMX 69 RCMX 93 |
MT-CYB |
MW275909 MW275910 MW275911 |
2n = 38 | La Sepultura Biosphere Reserve |
Xantusiidae | |||||
Lepidophyma flavimaculatum | RCMX 207 (OFV 1177) | MT-CYB | – | 2n = 38 (18M + 20m) | Montes Azules Biosphere Reserve |
RCMX 208 (OFV 1178) | – | ||||
RCMX 212 (OFV 1179) | MW275912 | ||||
RCMX 213 (OFV 1180) | MW275913 | ||||
RCMX 232 (OFV 1255) | MW275914 | ||||
Teiidae | |||||
Aspidoscelis deppii | RCMX 76 | MT-CYB | MW275915 | 2n = 52 (28M + 24m) | La Sepultura Biosphere Reserve |
Holcosus festivus | RCMX 223 (OFV 1213) | MT-ND2 | MW275916 | 2n = 50 (26M + 24m) | Montes Azules Biosphere Reserve |
RCMX 224 (OFV 1214) | – | ||||
RCMX 233 | MW275917 | ||||
Holcosus undulatus | RCMX 77 | MT-ND2 | MW275918 | 2n = 50 (26M + 24m) | La Sepultura Biosphere Reserve |
Dactyloidae | |||||
Anolis capito | RCMX 217 (OFV 1203) | MT-ND2 | MW275927 | 2n = 42 (24M + 18m) | Montes Azules Biosphere Reserve |
RCMX 218 (OFV 1204) | MW275928 | ||||
Anolis lemurinus | RCMX 214 (OFV1186) | MT-ND2 | MW275930 | 2n = 40 (24M + 16m) | Montes Azules Biosphere Reserve |
RCMX 225 (OFV 1215) | MW275929 | ||||
Anolis uniformis | RCMX 201 (OFV 1160) | MT-ND2 | MW275919 | 2n = 29/30 (14M + 15/16m) | Montes Azules Biosphere Reserve |
RCMX 203 | MW275925 | ||||
RCMX 205 (OFV 1164) | MW275926 | ||||
RCMX 206 (OFV 1176) | MW275920 | ||||
RCMX 209 (OFV 1183) | MW275921 | ||||
RCMX 210 (OFV 1173) | MW275922 | ||||
RCMX 215 (OFV 1182) | MW275923 | ||||
RCMX 226 (OFV 1211) | MW275924 |
The Mexican herpetofauna includes seven Scincella species that formerly belonged to the genus Sphenomorphus Fitzinger, 1843. They were reassigned to Scincella based on molecular phylogenetic analyses (
This species is distributed from Colima state, Mexico, southwards to Chiapas state, on the Pacific coast, and towards the southwest to Guatemala and Honduras.
RCMX85 (male*), RCMX86 (female*) and RCMX92 (female*) from La Sepultura, Chiapas, Mexico.
See below under S. cherriei (Cope, 1893) account.
The karyotype, described in
This species inhabits Mexico, from central Veracruz to extreme southeastern Puebla, northern Oaxaca state, southwards to Central America on the Atlantic coast, including the Yucatan Peninsula in México, reaching the eastern Panama.
RCMX219 (male) and RCMX235 (male*) from Estación Chajul, Selva Lacandona, Montes Azules, Chiapas state, Mexico.
The BI phylogenetic tree has been performed on 448-bp alignment of the 16S gene for four individuals of Scincella cherriei [RCMX219 and RCMX235 from the Montes Azules, one from Costa Rica (JF498076) and one from Nicaragua (AB057392)] and three individuals of Scincella assata [RCMX92 from La Sepultura, and two from El Salvador (JF498074 and JF498075)]. Scincella lateralis (Say, 1822) (AB057402 and JF498077) and S. reevesii (Gray, 1838) (JF498078) were used as outgroups. The tree (Fig.
The genetic divergence between the two specimens of S. cherriei from the Montes Azules and S. cherriei from other localities is high (5.2%), comparable to the divergence between S. assata and S. cherriei (6.6%-6.2%). The nominal subspecies S. c. cherriei (Cope, 1893), was described from Palmar, Costa Rica, which is far from from the Montes Azules. The lineage of S. cherriei from the Montes Azules may represent a different taxon worthy of additional detailed morphological and genetic studies.
The karyotype, described in
The genus Phyllodactylus is now constrained to the New World (
provisional distribution of this lineage, probably representing an undescribed species, is restricted to Pacific coast of eastern Oaxaca and western Chiapas states, Mexico (
RCMX67 (female*), RCMX69 (male*) and RCMX93 (female*) from La Sepultura, Chiapas state, Mexico.
The first description of the karyotype of one species of the P. tuberculosus complex is reported here (Fig.
As previously reported, 2n = 38 is the most common karyotype found in species of the genus Phyllodactylus from the Pacific coast of Mexico (
The genus Lepidophyma comprises 20 recognized species and is particularly speciose in Mexico, where 15 species are endemic and, in some cases, restricted to a particular mountain landscape (
Found on the Gulf of Mexico coast from Veracruz and Oaxaca, crossing the base of the Yucatan peninsula, through Central America to Panama.
RCMX207 (female*), RCMX208 (male*), RCMX212 (female*), RCMX213 (male*), and RCMX232 (female*) from Montes Azules, Chiapas state, Mexico.
Our samples have been identified on a morphological basis as Lepidophyma flavimaculatum, a species already reported for Chiapas. We aligned our 309 bp MT-CYB sequences to the 14 haplotypes of the same species published in
Diploid chromosome complements vary from 2n = 24 to 2n = 40 in Xantusiidae (
Species of the genus Aspidoscelis were previously included in Cnemidophorus Wagler, 1830, but based upon divergent morphological, molecular, and enzymatic characters the two genera were separated (
The species groups differ also in their karyotypes. 2n = 52 is observed in the deppii group, 2n = 46 in the sexlineata group, and 2n = 46 with XY sex chromosomal system in the tigris group.
The species has a wide distribution from Morelos and Michoacan (Mexico) south to Guatemala, El Salvador, Honduras, Nicaragua and Costa Rica.
RCMX76 (female*) from La Sepultura, Chiapas, Mexico.
The MT-CYB sequence (294-bp) is 4% divergent from GenBank sequences of Aspidoscelis deppii (KF555517-21) from Mexico (Playa Miramar, Tabasco). Despite the wide distribution, there are no studies on the intraspecific genetic variability of this species. It is a pity because this slight divergence in the MT-CYB could match with a different karyotype (see below).
In the genus Aspidoscelis chromosomal number ranges from 2n = 44 to 2n = 56, with some species showing triploid numbers, such as Aspidoscelis tesselatus (Say, 1823), with 69 chromosomes (
Ten species formerly assigned to the genus Ameiva F. Meyer, 1795 have been reassigned to the genus Holcosus and reorganized in three species groups (
This species is found in the lowlands of Tabasco and Mexico down to Colombia; it does not enter in the Yucatan Peninsula.
RCMX223 (female*), RCMX224 (female*), and RCMX233 (female) from Estación Chajul, Selva Lacandona, Montes Azules, Chiapas, Mexico.
The 600-bp PCR-amplified fragments of the MT-ND2 gene were identical in the two specimens (RCMX223 and RCMX233). The BLASTn search showed that this sequence belongs to Holcosus festivus, with 99.8% – 100% identity to H. festivus (KR058107, Montes Azules) and 96% identity to the other two H. festivus samples (KR058105 and KR058106, Costa Rica).
Here we report the first karyotype description for H. festivus (Fig.
The species is distributed along both coasts of Mexico from southern Nayarit to northern Costa Rica Pacific coast) and from southern Tamaulipas to central Nicaragua (Atlantic coast) including the peninsula of Yucatan.
RCMX77 (female*) from La Sepultura, Chiapas, Mexico.
The MT-ND2 sequence (556-bp) obtained from the individual from Chiapas has a 99% match to two GenBank sequences of H. undulatus parvus Barbour et Noble, 1915 (KR058051 and KR058063). According to Meza-Lazaro and Nieto-Montes de Oca (2015), this subspecies, distributed in the Pacific coast region of Southern Mexico and Northern Guatemala, should be elevated to species rank.
The specimen analyzed shows a 2n = 50 chromosome number (Fig.
Anolis (sensu lato) is the most speciose genus among the reptiles, with about 380 recognized species that have been all enclosed in a complete molecular phylogenetic tree by
The ancestral karyotype of “beta” Anolis (Norops) consists of 2n = 28 or 2n = 30 chromosomes subdivided in 14 macro- and 14 or 16 microchromosomes without evident sex chromosome heteromorphism (
Anolis capito has been found from Tabasco and northern Chiapas south to Central America on the Atlantic coast, to Costa Rica and Panama, where it is found on both coasts.
RCMX217 (female*), RCMX218 (female*) from Montes Azules, Chiapas, Mexico. The specimens were collected close to the northern part of species range and morphologically assigned to Anolis capito. Based on morphological studies from populations of almost all the species range, there is no evidence of cryptic species in A. capito (
We obtained a 685-bp MT-ND2 sequence showing 9% genetic divergence respect to an A. capito sequence collected in Costa Rica (GenBank AY909744). Such a high genetic divergence spurred us to perform a complete phylogenetic analysis with the MT-ND2 gene of Anolis species available in GenBank (not shown). The sequences from our samples cluster with the GenBank A. capito sequence, and together were sister to A. tropidonotus Peters, 1863. This tree topology has been already reported by
The specimens presently studied show, along with Anolis nebuloides Bocourt, 1973, the highest diploid number within the genus Anolis. The macrochromosomes include one pair of metacentric, six pairs of submetacentric, and five pairs of subtelocentric/acrocentric chromosomes. The chromosome shape of two pairs of microchromosomes appears to be biarmed. No heteromorphic sex chromosomes are discernible (unfortunately, no males have been analyzed).
The lack of description of chromosome morphology in Gorman’s study (
Occurs on the Atlantic coast from central Veracruz to central Panama, and on the Pacific coast from Costa Rica to central Panama.
RCMX214 (male*), RCMX225 (male*) Estación Chajul, Selva Lacandona, Montes Azules, Chiapas, Mexico.
BLAST analysis of the 630-bp MT-ND2 gene sequences from both individuals show 99.5% – 100% of identity with a sequence of A. lemurinus from Oaxaca (GenBank KT724761).
No previous chromosomal data are available for A. lemurinus and its karyotype is here described for the first time. Both male specimens from Montes Azules have a 2n = 40 (24M + 16m) karyotype (Fig.
This karyotype has the same composition in micro- and macrochromosomes as all Anolis species with 2n = 40 so far described. Molecular phylogenetics (
Occurs from southern Tamaulipas to north-central Honduras on the Atlantic coast.
RCMX201 (male), RCMX203 (male), RCMX205 (male*), RCMX206 (female*), RCMX209 (female), RCMX210 (male*), RCMX215 (male*) and RCMX226 (female*) from Estación Chajul, Selva Lacandona, Montes Azules, Chiapas, Mexico.
We report here the first description of the karyotype of this species (Fig.
Among the species of the genus Anolis with a known karyotype, this species is phylogenetically close to A. aquaticus Taylor, 1956 and A. biporcatus. Furthermore, A. biporcatus has also a similar composition of the sex chromosomes system, even if the morphology of sex chromosomes is different. In fact, the so-called 2n = 30 karyotype is one of the most common karyotypes in Anolis. However, three variants of this karyotype, based on the number and shape of macro- and microchromosomes, have been described. Among them, two types of 2n = 29–30 are present, type-A and type-B (
The type-A, typical of A. biporcatus, presents a multiple sex chromosomes system where X1 is an acrocentric chromosome, X2 is a microchromosome, and Y is metacentric similar in size to X1.
In our case, the Y is a small acrocentric chromosome, which might have been derived from a pericentric inversion in the submetacentric Y chromosome of the 2n = 29–30 type-A karyotype.Thus, although it is believed that the onset of multiple sex chromosomes in Anolis occurs independently (
Combined karyological and DNA taxonomic approaches have allowed us to highlight some interesting taxonomic peculiarities in 10 Mesoamerican lizard species belonging to six genera and five families. The karyotypes of four species, Phyllodactylus sp. 3 (P. tuberculosus species group), Holcosus festivus, Anolis lemurinus, and A. uniformis are here described for the first time. In Aspidoscelis deppii and Anolis capito, we found different karyotypes from those previously reported for these species. Moreover, in A. capito, the cytogenetic observation is consistent with the considerable genetic divergence at the studied mtDNA marker (MT-ND2), which is indicative of a putative new cryptic species. The anole species here studied exhibited different sex chromosomes configurations including a X1X1X2X2/X1X2Y condition in A. uniformis that should be in future studied by molecular cytogenetic techniques.
Another species that may include cryptic taxa is the skink Scincella cherriei, for which we found high values of genetic divergence among the specimens from Montes Azules and those from Costa Rica and Nicaragua, comparable to the divergence typical of sister species in skinks. A lower level of genetic divergence, compatible with an intraspecific phylogeographic structure, has been identified for L. flavimaculatum. In fact, the studied specimens belong to a mtDNA lineage that is sister with respect to the remaining haplotypes from other populations. However, it should be noted that the novel data represent only the first step in the identification of cryptic species and more efforts are necessary to investigate our assumptions. Both taxonomic revision and the notions related to the chromosome evolution in this hyper-diversified group of reptiles will be worthy of note.
We are grateful to Julia Carabias, Javier de la Maza, Rodrigo León, Diego Noriega, and to all the staff of the Estación Chajul Biological Research Station for their precious help during the permanence of RC and OFV in Montes Azules Biosphere Reserve. Secretaria de Medio Ambiente y Recursos Naturales of Mexico – SEMARNAT provided collection permit (N° FAUT-0015). Thanks are extended to Aaron Bauer and to the anonymous reviewer for their valuable comments on the first version of the manuscript, and to Giulia Trauzzi for help in the laboratory. RC received support grant from the University of Rome “La Sapienza” (Progetti di ricerca di Università); OFV received support grant from Dirección General de Asuntos del Personal Academico – UNAM for sabbatical period; AMRB received research fellowships from Conselho Nacional de Desenvolvimento Científico e Tecnológico – CNPq (PDJ 150599/2008-0, PCI-DA 300670/2019-2, and 301429/2020-0). Part of this manuscript was conceived and written while OFV was in Rome, founded by University Sapienza as Visiting Professor (bando 2018/2019).