Research Article |
Corresponding author: Omar Torres-Carvajal ( omartorcar@gmail.com ) Academic editor: Sven Bacher
© 2015 Omar Torres-Carvajal.
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:
Torres-Carvajal O (2015) On the origin of South American populations of the common house gecko (Gekkonidae: Hemidactylus frenatus). NeoBiota 27: 69-79. doi: 10.3897/neobiota.27.5437
|
Hemidactylus frenatus is an Asian gecko species that has invaded many tropical regions to become one of the most widespread lizards worldwide. This species has dispersed across the Pacific Ocean to reach Hawaii and subsequently Mexico and other Central American countries. More recently, it has been reported from northwestern South America. Using 12S and cytb mitochondrial DNA sequences I found that South American and Galápagos haplotypes are identical to those from Hawaii and Papua New Guinea, suggesting a common Melanesian origin for both Hawaii and South America. Literature records suggest that H. frenatus arrived in Colombia around the mid-‘90s, dispersed south into Ecuador in less than five years, and arrived in the Galápagos about one decade later.
Galápagos, Gekkonidae, Hemidactylus, invasive species, South America
With more than 120 species occurring in warm regions worldwide, Hemidactylus accounts for nearly 13% of the total number of recognized species in the family Gekkonidae. It is one of the most species-rich and widely distributed reptile clades (
The natural range of many of these widespread species is sometimes not clear. The common house gecko, H. frenatus Duméril & Bibron, 1836, is thought to have originated in tropical Asia and possibly the Indo-Pacific (
Distribution of Hemidactylus frenatus in South America and adjacent islands. Areas above 2000 m are shaded in dark grey. Localities of samples included in this study are shown in red. Locality data was taken from
I obtained nucleotide (nt) sequences of the mitochondrial ribosomal small subunit (12S, 370 nt) and cytochrome b (cytb, 303 nt) genes from 15 specimens collected on both sides of the Andes in Ecuador, as well as the Galápagos islands, and deposited in the herpetological collection of Museo de Zoología, Pontificia Universidad Católica del Ecuador (QCAZ). In addition, I retrieved sequences from GenBank representing samples of Hemidactylus frenatus from Colombia, Hawaii, India, Myanmar and Papua New Guinea. I used H. brookii and H. flaviviridis as outgroup taxa because they are closely related to H. frenatus (
Vouchers, locality data, and GenBank accession numbers of taxa and gene regions included in this study. Geographical coordinates in decimal degrees are provided for new localities sampled in this study.
Taxon | Voucher | Locality | GenBank accession number | |
---|---|---|---|---|
cytb | 12S | |||
Hemidactylus brookii | E1109.10 | India | DQ120276 | DQ120447 |
H. flaviviridis | E912.2 | Yemen | DQ120284 | DQ120455 |
H. frenatus | CES07035 | India | HM595655 | HM595691 |
H. frenatus | E509.5 | India | DQ120282 | DQ120453 |
H. frenatus | E509.2 | Myanmar | DQ120281 | DQ120452 |
H. frenatus | E509.1 | Myanmar | DQ120280 | DQ120451 |
H. frenatus | NV | Papua New Guinea | AY217801 | AY218005 |
H. frenatus | E509.7 | Hawaii | DQ120278 | DQ120449 |
H. frenatus | E509.6 | Hawaii | DQ120277 | DQ120448 |
H. frenatus | E509.3 | Colombia | DQ120279 | DQ120450 |
H. frenatus | QCAZ4524 | Ecuador: Esmeraldas 1.0425; -78.6304 |
KT455016 | KT455031 |
H. frenatus | QCAZ4875 | Ecuador: Manabí -0.9505; -80.7423 |
KT455017 | KT455032 |
H. frenatus | QCAZ5076 | Ecuador: Esmeraldas 0.8740; -79.8450 |
KT455018 | KT455033 |
H. frenatus | QCAZ8124 | Ecuador: Pastaza -1.4529; -77.4425 |
KT455019 | KT455034 |
H. frenatus | QCAZ8130 | Ecuador: Pastaza -1.4529; -77.4425 |
KT455020 | KT455035 |
H. frenatus | QCAZ8472 | Ecuador: Guayas -2.2126; -79.4472 |
KT455021 | KT455036 |
H. frenatus | QCAZ9111 | Ecuador: Guayas -2.1822; -80.0181 |
KT455022 | KT455037 |
H. frenatus | QCAZ10197 | Ecuador: Pichincha -0.5888; -79.3627 |
KT455023 | KT455038 |
H. frenatus | QCAZ10213 | Ecuador: Orellana -0.4720; -76.9807 |
KT455024 | KT455039 |
H. frenatus | QCAZ10215 | Ecuador: Orellana -0.4720; -76.9807 |
KT455025 | KT455040 |
H. frenatus | QCAZ11128 | Ecuador: Galápagos -0.9573; -90.9674 |
KT455026 | KT455041 |
H. frenatus | QCAZ11165 | Ecuador: Galápagos -0.9573; -90.9674 |
KT455027 | KT455042 |
H. frenatus | QCAZ11197 | Ecuador: Galápagos -0.9573; -90.9674 |
KT455028 | KT455043 |
H. frenatus | QCAZ11452 | Ecuador: Manabí 0.0740; -80.0480 |
KT455029 | KT455044 |
H. frenatus | QCAZ11593 | Ecuador: Orellana -0.9167; -75.4000 |
KT455030 | KT455045 |
Genomic DNA was isolated from frozen muscle or liver tissues using a guanidinium isothiocyanate extraction protocol. Polymerase Chain Reaction (PCR) amplification of gene fragments was performed in a final volume of 25 µl reactions using 1X PCR Buffer (– Mg), 3 mM MgCl2, 0.2 mM dNTP mix, 0.2 µM of each primer, 0.1 U/µl of Platinum® Taq DNA Polymerase (Invitrogen, Carlsbad, CA) and 1 µl of extracted DNA. Negative controls were run on all amplifications to check for contamination. Gene fragments were amplified using the primers 12S1L, 12S2H and 12sb for 12S (
Data were aligned in MAFFT under default settings (
A total of 673 aligned sites of 12S (370 nt) and cytb (303 nt) were obtained. Selected models were K80+G, 000010+F, TrN, and 012212+G+F for 1st, 2nd and 3rd codon positions of cytb and 12S, respectively. Of the 25 aligned sequences, 17 contained missing data ranging between 1–65 sites, which together represented 2.6% of the total sites in the matrix.
All 15 new sequences obtained in this study from continental Ecuador and the Galápagos, as well as GenBank sequences of specimens from Colombia, Hawaii and Papua New Guinea were identical (missing data ignored). The clade formed by these sequences (PP=0.95) was recovered with high support (PP=0.91) as sister to a clade with two samples from Myanmar (PP=0.90); samples from India were nested in a clade (PP=0.78) sister to all other samples of Hemidactylus frenatus (Fig.
Intraspecific genetic distances between individuals of H. frenatus from Papua New Guinea/Hawaii/South America and individuals from India and Myanmar included in the analysis varied between 0.069–0.132 and 0.047–0.078 for cytb and 12S, respectively (Table
Majority rule (50%) consensus tree of 36,000 trees obtained from a Bayesian analysis of 25 specimens and 673 nucleotides corresponding to cytb and 12S mitochondrial gene regions. Outgroup taxa Hemidactylus brookii and H. flaviviridis are not shown; all terminals correspond to H. frenatus. Numbers on branches are posterior probability values. Voucher numbers (if available) and country of collection are shown on each terminal. For samples collected in Ecuador, voucher numbers and province names are indicated.
Uncorrected genetic distances among taxa included in this study for cytb (upper diagonal) and 12S (lower diagonal) gene fragments. Taxon name along with voucher number and country of collection are indicated in first column. The sequence of H. frenatus from Colombia represents other sequences from South America, as well as those from Hawaii and Papua New Guinea included in this study.
Taxon sample | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
---|---|---|---|---|---|---|---|
1: H. flaviviridis | |||||||
E912.2 Yemen | 0.251 | 0.201 | 0.205 | 0.191 | 0.218 | 0.187 | |
2: H. brookii | |||||||
E1109.10 India | 0.208 | 0.178 | 0.195 | 0.162 | 0.185 | 0.189 | |
3: H. frenatus | |||||||
E509.1 Myanmar | 0.172 | 0.158 | 0.076 | 0.076 | 0.129 | 0.130 | |
4: H. frenatus | |||||||
E509.2 Myanmar | 0.174 | 0.160 | 0.025 | 0.069 | 0.135 | 0.105 | |
5: H. frenatus | |||||||
E509.3 Colombia | 0.175 | 0.158 | 0.047 | 0.047 | 0.132 | 0.105 | |
6: H. frenatus | |||||||
E509.5 India | 0.168 | 0.151 | 0.056 | 0.064 | 0.078 | 0.105 | |
7: H. frenatus | |||||||
CES07035 India | 0.177 | 0.164 | 0.078 | 0.081 | 0.070 | 0.061 |
The fact that cytb and 12S haplotypes of Hemidactylus frenatus from Hawaii, the Galápagos, and mainland South America are identical to those in Papua New Guinea sheds some light on the origin and dispersal of this species from Melanesia to South America across the Pacific Ocean. The invasive populations in Hawaii and South America most likely originated from a single ‘stock’ in Melanesia; otherwise, we would expect more genetic variation among invasive samples. This is supported by the genetic variation that was observed only among the four samples from India and Myanmar included in this study. Although these four samples come from geographically close localities, their genetic distances vary between 0.076–0.135 and 0.025–0.064 for cytb and 12S, respectively (Table
How Hemidactylus frenatus arrived in South America remains an open question given its ability for massive, human-mediated range expansion (
The common house gecko, as its name suggests, is easy to spot at human settlements feeding around light bulbs at night. Therefore, we can assume that the first time it is found as an invasive species in a certain location corresponds approximately to the colonization time at that location. Thus, based on the year H. frenatus was first reported from each site included in this study, it seems like its general colonization route after leaving Papua New Guinea or somewhere nearby was (first report year follows each site) Hawaii 1940s (
Recent arrival of Hemidactylus frenatus to South America including iconic conservation sites as the Galápagos archipelago should be of concern. Given its great dispersal ability and potential distribution (
Only recently was Hemidactylus frenatus reported for the first time in Galápagos (
I thank S. Báez, S. Ron, I. Tapia and many PUCE biology students for their help in the field; L. Jaramillo and G. Castillo for laboratory work; S. Lobos and S. Quintero for producing the map; W. Tapia for help with logistics in the Galápagos islands. Special thanks to M. Tonione for improving earlier versions of this manuscript with helpful comments. Collection of samples was authorized by permits 008-09 IC-FAU-DNB/MA, 001-10 IC-FAU-DNB/MA, 001-11IC-FAU-DNB/MA and PC-05-10 issued by Ministerio del Ambiente de Ecuador and Galápagos National Park. This research was funded by Pontificia Universidad Católica del Ecuador and Secretaría de Educación Superior, Ciencia, Tecnología e Innovación del Ecuador (SENESCYT).