High levels of population genetic differentiation in the American crocodile (Crocodylus acutus)

Abstract

The American crocodile (Crocodylus acutus) is a widely distributed species across coastal and brackish areas of the Neotropical region of the Americas and the Greater Antilles. Available information on patterns of genetic differentiation in C. acutus shows a complex structuring influenced by interspecific interactions (mainly hybridization) and anthropogenic actions (mostly historical hunting, recent poaching, habitat loss and fragmentation, and unintentional translocation of individuals). In this study, we used data on mitochondrial DNA control region and 11 nuclear polymorphic microsatellite loci to assess the degree of population structure of C. acutus in South America, North America, Central America and the Greater Antilles. We used traditional genetic differentiation indices, Bayesian clustering and multivariate methods to create a more comprehensive picture of the genetic relationships within the species across its range. Analyses of mtDNA and microsatellite loci show evidence of a strong population genetic structure in the American crocodile, with unique populations in each sampling locality. Our results support previous findings showing large degrees of genetic differentiation between the continental and the Greater Antillean C. acutus. We report three new haplotypes unique to Venezuela, which are considerably less distant from the Central and North American haplotypes than to the Greater Antillean ones. Our findings reveal genetic population differentiation between Cuban and Jamaican C. acutus and offer the first evidence of strong genetic differentiation among the populations of Greater Antillean C. acutus.

About the Authors:

Natalia A. Rossi

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft

* E-mail: [email protected]

Affiliations Wildlife Conservation Society, Bronx, New York, United States of America, Sackler Institute of Comparative Genomics, American Museum of Natural History, Manhattan, New York, United States of America, Department of Ecology, Evolution and Environmental Biology, Columbia University, Columbia, New York, United States of America

ORCID logo http://orcid.org/0000-0002-9026-0740

Angelica Menchaca-Rodriguez

Roles Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – review & editing

Affiliation: School of Biological Sciences, the University of Bristol, Bristol, England, United Kingdom

Rafael Antelo

Roles Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – review & editing

Affiliations Fundación Palmarito Casanare, Bogotá, Colombia, Estación Biológica El Frío, Apure, Venezuela

Byron Wilson

Roles Data curation, Investigation, Methodology, Writing – review & editing

Affiliation: Department of Life Sciences, University of the West Indies, Mona, Jamaica

Kurt McLaren

Roles Data curation, Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Life Sciences, University of the West Indies, Mona, Jamaica

Frank Mazzotti

Roles Data curation, Formal analysis, Investigation, Writing – review & editing

Affiliation: Fort Lauderdale Research and Education Center, University of Florida, Fort Lauderdale, Florida, United States of America

Rafael Crespo

† Deceased.

Roles Data curation, Formal analysis, Writing – review & editing

Affiliation: Fort Lauderdale Research and Education Center, University of Florida, Fort Lauderdale, Florida, United States of America

Joe Wasilewski

Roles Data curation, Writing – review & editing

Affiliation: Fort Lauderdale Research and Education Center, University of Florida, Fort Lauderdale, Florida, United States of America

Fernando Alda

Roles Data curation, Software, Writing – review & editing

Affiliation: Department of Biology, Geology and Environmental Science, University of Tennessee at Chattanooga, Chattanooga, Tennessee, United States of America

Ignacio Doadrio

Roles Data curation, Formal analysis, Writing – review & editing

Affiliation: Departamento de Biodiversidad y Biología Evolutiva, Museo Nacional de Ciencias Naturales, Madrid, Spain

Tito R. Barros

Roles Data curation, Validation, Writing – review & editing

Affiliation: Museo de Biología de la Universidad del Zulia, Maracaibo, Venezuela

Evon Hekkala

Roles Data curation, Methodology, Writing – review & editing

Affiliation: Department of Biological Sciences, Fordham University, Bronx, New York, United States of America

Manuel Alonso-Tabet

Roles Conceptualization, Data curation, Investigation, Methodology, Writing – review & editing

Affiliation: Refugio de Fauna Monte Cabaniguán-Ojo de Agua, Las Tunas, Cuba

Yairen Alonso-Giménez

Roles Conceptualization, Data curation, Writing – review & editing

Affiliation: Refugio de Fauna Monte Cabaniguán-Ojo de Agua, Las Tunas, Cuba

Manuel Lopez

Roles Data curation, Methodology, Writing – review & editing

Affiliation: Refugio de Fauna Monte Cabaniguán-Ojo de Agua, Las Tunas, Cuba

Georgina Espinosa-Lopez

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

Affiliation: Facultad de Biología, Universidad de la Habana, La Habana, Cuba

Joe Burgess

Roles Data curation, Writing – review & editing

Affiliation: US Forest Service, Idleyld Park, Oregon, United States of America

John B. Thorbjarnarson

† Deceased.

Roles Conceptualization, Investigation, Project administration, Supervision, Validation

Affiliation: Wildlife Conservation Society, Bronx, New York, United States of America

Joshua R. Ginsberg

Roles Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Resources, Writing – review & editing

Affiliation: Cary Institute of Ecosystem Studies, Millbrook, New York, United States of America

Kent A. Vliet

Roles Conceptualization, Data curation, Investigation, Methodology, Project administration, Resources, Writing – review & editing

Affiliation: Department of Biology, University of Florida, Gainesville, Florida, United States of America

George Amato

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – review & editing

Affiliation: Sackler Institute of Comparative Genomics, American Museum of Natural History, Manhattan, New York, United States of America

Introduction

The American crocodile (Crocodylus acutus) is widely distributed along coastal swamps, estuarine rivers, and lakes of the Neotropical region of the Americas and the Insular Caribbean. It is primarily a coastal crocodilian inhabiting brackish waters and the saltwater sections of rivers, coastal lagoons, and mangrove swamps [1,2]. The species also occurs in freshwater areas located inland [3–6]. This wide-ranging species overlaps with three other new world crocodilians: the Cuban crocodile (Crocodylus rhombifer) in Zapata, Cuba [7], Morelet’s crocodile (Crocodylus moreletii) in the Yucatan Peninsula, Mexico [8,9] and Belize [10], and the Orinoco crocodile (Crocodylus intermedius) in Venezuela [11]. It hybridizes in the wild with C. moreletii [12–14] and C. rhombifer [7,13]. The ability to tolerate and thrive in highly saline environments [15–17] allowed C. acutus to occupy a unique ecological niche and extend geographically more than any other new world crocodilian.

Like many crocodilians, populations of American crocodiles have declined due to intense hunting pressures since the early 1800s. The American crocodile was listed in the Convention on International Trade in Endangered Species of Flora and Fauna (CITES) at the Convention's inception in 1975 [18]. Since then, a combination of sustainable farming [18], habitat protection [19], and conservation and monitoring programs [20] have aided the recovery of the species. However, this recovery has not been even throughout its range, and new conservation challenges have arisen. New threats have been driven by changes in water regimes and water quality of coastal wetlands, illegal hunting [21], retaliatory killing, destruction of nests, inter- and intraspecific competition [14,22], human-mediated movements of animals across the range [23], and habitat loss and fragmentation due to coastal development [10,24]. Habitat reduction and degradation from coastal development is the leading threat to the survival of the American crocodile, as it decreases the availability of nesting sites, food resources, nursery habitat and hiding sites [8,25–27], and even contributing to increased hybridization through the breakdown of reproductive barriers and increase potential of panmixia among distinctive evolutionary lineages, leading to a decline of genetic diversity and the extinction of local populations [28]. The species is categorized as vulnerable by the International Union for Conservation of Nature (IUCN) Red List of Endangered Species [26] and listed in CITES Appendix I [19], which prohibits trade, except for the populations of Cuba [1], the Integrated Management District of Mangroves of the Bay of Cispata in Colombia, and the population of Mexico (with no export quota allowed of wild specimens for commercial purposes), all listed in Appendix II [29].

An important challenge in the conservation of any widespread species is to identify unique clusters that need to be considered as independent management units [30]. These clusters are usually defined considering habitat types within the species range, local conservation status, threats, enforcement, and legislation [19,31], as well as patterns of genetic variation within and among subpopulations. The majority of the work on C. acutus genetic variation has focused on zones where it can hybridize in the wild with C. rhombifer and C. moreletti [7,9,12,13,32,33]. Though information about genetic differentiation has been used as a tool for delineating local management units on some species of Crocodylus [34], there is still a need to incorporate population-level genetic information into conservation planning for many species of crocodilians, including the American crocodile.

Current information on patterns of genetic differentiation of C. acutus shows a complex population structuring influenced by interspecific interactions (mainly hybridization) and anthropogenic actions (mainly poaching, habitat destruction and fragmentation, and unintentional translocation of individuals). For example, an analysis of variation in the mtDNA control region revealed a minimum of 11 haplotypes in the American crocodile [23,32]. Nine of these haplotypes are also present in Central and North American populations [23], and only two are exclusive to the Greater Antilles [7,32,35]. C. acutus populations from South America, however, have not been incorporated in these analyses, and thus, the relationship of Antillean crocodiles with South American lineages is unclear. Moreover, hybridization with other crocodilians and human-mediated migration across the species range may affect the distribution and frequency of the reported haplotypes [23]. A study based on mitochondrial (mtDNA) and nuclear DNA markers revealed interspecific hybridization, admixture, and significant patterns of population substructure of crocodiles in Southern Florida as a result of human-mediated migrations from Latin America and the Greater Antilles [23]. More recently, a study on populations of C. acutus in the Greater Antilles and the Americas recognized distinguishable phylogenetic relationships and high genetic divergence between these two groups [36]. C. acutus populations along the Pacific coast of Costa Rica exhibited small to moderate levels of inbreeding and significant levels of population differentiation potentially attributed to the lack of connectivity between some localities and the occurrence of population bottlenecks in the past [37]. In addition, Pacheco et al [33] extensive sampling throughout Mexico revealed several unique lineages as a result of historical and present hybridization between C. acutus and C moreletti. In parallel to C. acutus’ species designation debates [33,36], it is imperative to strengthen the comparative research of extant populations to inform local and regional management and incorporate lacking information into the broader species assessment.

Defining conservation management units based on molecular data has proven more effective when combining a wide range of genetic markers [38–40]. The hypervariable control region of the mitochondrial DNA has been one of the most informative and widely used markers to uncover population subdivision [38,41]. In addition, nuclear polymorphic microsatellite markers have proven useful to detect subtle patterns of population structure [42–44], and to accurately assign population origin [34]. Information on population differentiation of American crocodiles derived from mtDNA markers has been useful to detect hierarchies of structuring but has not resolved subtle patterns of population subdivision [23].

To better understand patterns of genetic differentiation in the American crocodile, we used data on mtDNA and nuclear DNA to assess the degree of population structuring between and among localities in South America, North America, Central America, and the Greater Antilles. We incorporated previously unsampled areas (two sites in Jamaica, and five sites in Venezuela) and populations previously studied to create a more comprehensive picture of the genetic relationships within the species across its range. Our work aims at understanding the population genetic structuring of C. acutus, and to provide relevant information for the conservation and recovery of the species.

Materials and methods

Sample collection

As part of crocodilian conservation and monitoring programs, we collected tissue samples of C. acutus from five countries (United States, Belize, Venezuela, Southeastern Cuba, and Jamaica). The dataset was supplemented with mtDNA sequences obtained from Genbank from Mexico, Costa Rica, and Southwestern Cuba (Fig 1A; Table 1).

[Figure omitted. See PDF.]

Fig 1. Geographic distribution, haplotype network, discriminant analysis of principal components, and population structure of the American crocodile.

A) C. acutus distribution range based on most recent data of presence, and sampling localities color-coded by country (Base map credit: NASA Earth Observatory); B) Neighbor-joining haplotype network of the mtDNA control region, circle size is proportional to haplotype frequency, and hatch marks indicate mutation steps; C) scatterplot of DAPC with PC1 and PC2; D) estimated Bayesian genetic structures based on 11 microsatellite loci, with K = 3 and K = 4. Samples are colored according to geographic location.

https://doi.org/10.1371/journal.pone.0235288.g001

[Figure omitted. See PDF.]

Table 1. Geographic distribution and description of the American crocodile samples used in this study.

https://doi.org/10.1371/journal.pone.0235288.t001

Tissue samples were removed from the dorsal section of the base of the tail and stored at -20°C in 95% ethanol until DNA extraction. The sex of the individuals was unknown. Protocols for animal handling and collecting biological samples were approved by the Field Veterinary Program of the Wildlife Conservation Society and conducted following guidelines of the IUCN Crocodile Specialist Group [45]. Samples were obtained in accordance with collection permits from the National Environment and Planning Agency of Jamaica (Exemption Certificate N0 70), the Center for Inspection and Environmental Control of Cuba (License 18/15), the National Office of Biological Diversity of Venezuela (Permit number 5–0358), the Forest Department of Belize, and the Institutional Animal Care and Use Committee (IACUC) of the United States (Permits A037-2006 and D684-2007). Samples from all sites except Venezuela were transported to the Sackler Institute for Comparative Genomics at the American Museum of Natural History, New York, USA (CITES export permits C0001733 and JM2320, and USFWS permit 1018–0093) for laboratory analyses including DNA extraction, amplification and sequencing of the mtDNA control region, and amplification and sequencing of 11 microsatellite loci. Samples from Venezuela were transported to the Museo Nacional de Ciencias Naturales of Madrid, Spain (CITES permit 1395/VE9120190), where laboratory analysis, including mtDNA extraction, amplification, and sequencing were conducted. Nuclear DNA was not recovered from samples originating from Venezuela, and thus, microsatellite loci could not be included in our analyses.

DNA extraction, sequencing and genotyping

We extracted genomic DNA from fresh scale tissue using the DNeasy Blood and Tissue kit (Qiagen, Valencia, CA) following manufacturer recommendations. The final DNA concentration of each sample was measured using a NanoDrop 2000c. A fragment of 525 base pairs (bp) of the mitochondrial control region was amplified with primers drL15459 and CR2HA [46]. Polymerase chain reaction (PCR) was performed at a final volume of 25μl with 1.25 units of AmpliTaq Gold (Applied Biosystems, Foster City, California, USA), 1x PCR Buffer, 0.2mM each dNTP, 0.2 μM each primer, 1mM MgCl2 and 50ng/μl DNA. An initial 2-minute denaturation cycle at 94°C was followed by 33 cycles of denaturation at 94°C for 30 seconds, annealing at 56°C for 1min, and elongation at 72°C for 7 minutes. Amplicons were purified and concentrated via ethanol precipitation. Double-stranded DNA was sequenced with the dideoxy chain termination method using an ABI 3730 automated sequencer.

To genotype samples from Belize, Cuba, Jamaica, and the Everglades, we used a panel of 11 microsatellite markers previously developed for Crocodylus spp., including C391, Cj16, Cj18, Cj35, Cj127, Cj119, Cj131, CUJ131, Cu5123, Cj20, Cj109 [47]. PCR reactions were prepared at a final volume of 12.5μl with 1.0 unit of AmpliTaq Gold (Applied Biosystems, Foster City, California, USA), 1x PCR Buffer, 0.4mM each dNTP, 0.1 μM fluorescently labeled forward primer, 0.1 μM reverse primer, 3.5mM MgCl2 and 50ng/μl DNA. Microsatellites were amplified in single PCRs for 33 cycles and three different annealing temperatures (51 ºC for Cj35 and Cj127, 62 ºC for Cj16, Cj20, and Cj109, and 58°C for all the rest). An initial 2-minute denaturation cycle at 94°C was followed by 33 cycles of denaturation at 94°C for 30 seconds, different annealing temperature for each microsatellite at 51–62°C for 1-minute, elongation at 72°C for 45 seconds and a final elongation step at 72°C for 5 minutes. We visualized PCR products using an ABI 3730 automated sequencer. Genotypes were identified using GeneMapper v5.0 software (Applied Biosystems, Foster City, California, USA). Allelic sizes were scored against the size standard GS500 LIZ.

Haplotypic and genotypic variation

We aligned mtDNA sequences in GENEIOUS 8.1.7 (Biomatters, Ltd., San Francisco, CA, USA) using the Muscle Aligner 3.8.425 with default settings. From the 525 bp mtDNA control region fragment, we selected a 458 bp consensus region containing about 95% of the variation within our samples. We compared the mtDNA sequences obtained with those downloaded from Genbank corresponding to Mexico (AY568308-17)[48] Costa Rica (GU064561-65) [23], and Southwestern Cuba (EU034586 KM577701) [32]. We then matched sequences to previously described haplotypes [23] using DnaSP v5.10.1 software [49]. We used Rodriguez et al [23] haplotype definitions and nomenclature for our analysis. Haplotypic diversity (Hd) [50], nucleotide diversity (π) [51,52], and the mean number of pairwise differences among sequences (K) [51–53] were calculated in Arlequin 3.5 [54,55] and DnaSP. We used TCS 1.21 software to construct statistical parsimony haplotype networks [56] in order to depict an intraspecific genealogy for our mtDNA sequences. We then built a median-joining haplotype network [57] to depict patterns of genetic variation among haplotypes using the PopART 1.7 software [58].

We used Arlequin to calculate the number of alleles, observed (HO) and expected heterozygosities (HE) in our microsatellite marker set. We used GENEPOP 4.3 [59] to test Hardy-Weinberg (HW) expectations per locus as well as genotypic linkage disequilibrium (LD) between loci. Sequential Bonferroni corrections [60] were used to adjust departures for HW and LD (p≤0.05). We took into account that genotyping errors could occur as a result of primer-site mutations or contamination and could produce null alleles and/or allele drop out. Accordingly, genotyping inconsistencies were assessed using Micro-checker 2.2.3 software [61].

Analysis of population structure

Mitochondrial sequence data.

The mtDNA dataset included sequences from 265 individuals corresponding to thirteen locations in seven countries (Fig 1A; Table 1): 1) Everglades National Park, the United States, 2) Quintana Roo, Mexico, 3) Turneffe Atoll, Belize, 4) Costa Rica, 5) Chico River, Cuare, Burro Negro, Negro River, and Santa Rosa River, Venezuela, 6) Zapata Swamp National Park, and the Wildlife Refuge Monte Cabaniguan, Cuba, 7) Black River Lower Morass, and Portland Bight Protected Area, Jamaica.

The partitioning of genetic variation among putative populations was assessed through a nested analysis of molecular variance (AMOVA) implemented in Arlequin. This allows for testing hypotheses of among-group and within-group differences at several hierarchical levels. We computed genetic distances in Arlequin, including pairwise FST (haplotype frequencies only) and ΦST statistics (pairwise differences between haplotypes). The significance of values (Fst and ΦST) was tested using 10,000 nonparametric random permutations. Chi-square analysis was conducted in DnaSP to test for significant differences in haplotype frequency distributions among sampling localities.

Microsatellite data.

We analyzed microsatellite data from 225 individuals from four countries: 1) Turneffe Atoll, Belize, 2) Wildlife Refuge Monte Cabaniguan, Cuba, 3) Black River Lower Morass, and Portland Bight Protected Area, Jamaica, and 4) Everglades National Park, United States (Fig 1; Table 1, and S1 Table). We assessed spatial structure through AMOVA analysis, and estimated pairwise FST statistics in Arlequin, using Weir and Cockerham’s estimator [62], which assumes an infinite allele model of mutation. We evaluated the degree of genetic partitioning among putative populations using spatially explicit Bayesian clustering methods. Such statistical methods describe and quantify the geographic patterns of intraspecific genetic variation.

The package Geneland 3.1.4 [63,64] implemented in R 3.2.3 [65] uses an algorithm based on a spatial model which assigns individuals into a number of genetic clusters (K) making use of genotypes as well as spatial coordinates of sampled individuals. As is common to all explicit clustering methods, Geneland weights information of an individual's location in the search for the most likely K instead of assuming that all clustering solutions are equally possible. Whereas, STRUCTURE 2.3.4 [66–68] uses the number of populations (K) as a fixed parameter to estimate the log-likelihood of the data for the pre-defined K values, and assigns memberships for all individuals in the total sample set. Both programs infer unknown parameters through Markov Chain Monte-Carlo (MCMC) computations, assume Hardy-Weinberg equilibrium with linkage equilibrium between loci, and do not require an a priori definition of putative populations. STRUCTURE, too, allows for the incorporation of sampling location priors. Location priors allow STRUCTURE to assign individuals into genetically similar clusters considering a priori the geographic origin of each individual [69].

For Geneland, we determined K across 10 independent runs with 1,000,000 MCMC iterations and allowing K to vary from 1 to 8. We used a correlated allele model and set the maximum rate of the Poisson process at 225 (the number of individuals) and the maximum number of nuclei in the Poisson-Voronoi tessellation at 675 (three times the number of individuals) as suggested by Guillot et al. [63]. The uncertainty of spatial coordinates of the collection was set at 20 meters. In our study, the uncertainty of spatial coordinates accounted for the recording error. We inferred the most likely number of clusters as the modal K with the highest posterior probability.

We ran STRUCTURE without sampling location priors, through 20 independent runs for K = 1–8. We set a burn-in period of 100,000 and 1,000,000 MCMC iterations to identify the genetic clusters. We assumed an admixture model with correlated allele frequencies, which allows individuals to have mixed ancestry. We determined the optimal value of K according to the ΔK method [70] using STRUCTURE HARVESTER 0.6.94 online application [71]. We performed consensus analyses for the average scores for the inferred K value in CLUMPP 1.1.2 software [71].

Finally, we performed a Discriminant Analysis of Principal Components (DAPC) using the Adegenet 2.0.1 package [72] to further explore the structure of the population without making assumptions of panmixia and to complement Bayesian inference analyses. This method is useful to identify population clusters without prior knowledge and provides a robust alternative to traditional Bayesian approaches [73].

Results

Genetic diversity

Mitochondrial sequence data.

We identified a total of ten mtDNA haplotypes across populations of C. acutus in the Americas and Greater Antilles (Table 2). Four haplotypes had not been previously described (GeneBank accession numbers MT416744, MT416745, MT416746, and MT416747), while the other six matched previously identified ones [23,74]. Haplotype β, described by Cedeño-Vazquez et al. [74], was observed in the Greater Antillean population. The most common haplotype in the Everglades was Ca1; noteworthy, we also found two individuals with haplotype β and one with a newly reported one (Ca10). Crocodiles from Belize exhibited the two more common haplotypes for the species in Central and North America (Ca1 and Ca2). We found three haplotypes for South American crocodiles (Ca11, Ca12, and Ca13).

[Figure omitted. See PDF.]

Table 2. Genetic diversity indices for mtDNA data.

https://doi.org/10.1371/journal.pone.0235288.t002

Crocodiles from Jamaica and one locality in Cuba (WRMC) showed the lowest values of mtDNA genetic diversity (Table 2). Maternal lineages sampled from these populations seemed to be identical in terms of molecular distances, with the same mtDNA haplotype (β) and values of zero for pairwise differences among sequences and nucleotide diversity. Likewise, samples from Costa Rica had values of zero in all genetic diversity estimators and only one haplotype. However, as GeneBank available sequences encompassed only five samples, this may not represent the true number of haplotypes present in Costa Rica. Samples from Mexico exhibited the highest values of haplotypic diversity, followed by samples from Venezuela, Cuba (ZAP_CU), and the Everglades in the United States. Molecular distances, however, were relatively similar among Mexican samples and were higher for Zapata, the Everglades, and Venezuela.

We detected strong patterns of haplotype distribution across geographic locations. The median-joining haplotype network (Fig 1B) reveals differences in haplotype identity between the Americas and the Greater Antilles, in accordance with recent findings [36]. Cuba and Jamaica comprise a cluster separated from all other populations. Haplotypes present in C. acutus from the Everglades, Central America, and Venezuela are closely related than those in the Antilles.

Microsatellite data.

The analysis of microsatellite data did not find evidence of LD, and the null hypothesis of HW equilibrium could not be rejected (p > 0.05). We found no significant differences between the expected heterozygosity under HW and the observed heterozygosity in the data for any of the putative populations (Table 3). Loci Cj109 and Cj127 had the highest allele counts, while CUJ131 and CUD68 showed the lowest.

[Figure omitted. See PDF.]

Table 3. Genetic diversity indices of the microsatellite data.

https://doi.org/10.1371/journal.pone.0235288.t003

The results from Micro-checker showed no evidence of scoring error due to stuttering or large allele dropout and crocodile populations from the four countries sampled (United States, Cuba, Jamaica, and Belize). There was evidence of shared null alleles (Cj20, Cj35, Cj109, Cj131) in all localities. Microsatellite null alleles are commonly encountered in population genetics studies [75] but appear to have little effect on the outcome in Bayesian assignment analyses [76].

Population structure

Mitochondrial sequence data.

The among-groups component of the AMOVA analysis was significant when both the haplotype frequencies and the molecular distances were considered (FST: 0.96, p<0.0001; Φst = 0.75, p<0.001). Similarly, the χ2 showed significant differentiation among populations (χ2 = 659.03, p<0.001, df = 72). We found significant structure between the same pairs of populations when performing the exact test of population differentiation (Table 4). Pairwise comparisons showed significant structuring for all pairs except for Mexico—Belize, Mexico—Costa Rica, and Cuba—Jamaica. Greater Antillean C. acutus populations exhibited strong differentiation when compared with all other populations with fixation indices ranging from 0.45 to 1 (FST).

[Figure omitted. See PDF.]

Table 4. Pairwise estimates of FST and ΦST between sampling localities based on mtDNA.

https://doi.org/10.1371/journal.pone.0235288.t004

Microsatellite data.

The among-groups component of the AMOVA analysis was significant (FST: 0.20, p<0.0001). In contrast to the data from the mtDNA, pairwise comparisons between Greater Antillean populations from Cuba and Jamaica did reveal population structure; congruent to the mtDNA analysis, all remaining comparisons between population pairs were also significant (Table 5).

[Figure omitted. See PDF.]

Table 5. Pairwise estimates of FST between sampling localities based on mtDNA.

https://doi.org/10.1371/journal.pone.0235288.t005

Spatially explicit Bayesian clustering analysis implemented in Geneland yielded a modal number of four populations. All ten independent runs converged in four populations. The comparison of the posterior probability of assignment of individuals to populations allocated each individual to the population from which it originated. Non-spatial Bayesian clustering implemented in STRUCTURE with no prior distribution revealed subdivision of the samples into three populations (K = 3) with Delta K = 950.33 (Fig 1D); clustering individuals from the United States and Belize and separating Cuba from Jamaica. The second-best configuration assigned individuals into four populations: Belize, Cuba, Jamaica, and the United States with Delta K = 147.37 (Fig 1D).

The DACP (Fig 1C) suggested four distinct groups: (1) Belize, (2) Cuba (3) Jamaica, and (4) the United States. This grouping was based predominantly on the first principal component (PC1), which explained 19.23% of the variance in allele frequencies among samples and was augmented by PC2 (12.50%) and, to a lesser extent by PC3 (5.76%).

Discussion

We present new and relevant information to understanding C. acutus evolutionary history and reveal patterns of genetic partitioning important for its conservation. As a widespread species living across a diversity of fresh and brackish water habitats, able to move long distances across saltwater, the American crocodile has been considered "homogenous" across its range before population studies emerged in the early 2000s. With most research focusing on hybrid zones [7,9,19,32,74], our study contributes to current debates on the taxonomic complexity of C. acutus [33,35,36]. Overall, our findings improve our understanding of populations of C. acutus across seven countries and provide sound evidence of genetic structuring with direct consequences into the management of local, unique populations.

Analyses for both mtDNA and nuclear markers show evidence of population genetic structure in the American crocodile with unique populations in North America, Central America, South America, and the Greater Antilles. In accord with other studies [7,23,35,36], our results show the greatest degree of genetic differentiation between the continental and Greater Antillean C. acutus. The suite of genetic differentiation estimators used in our mtDNA analysis indicated that this differentiation is strong, reaching pairwise comparisons estimator values (FST and ΦST) as high as one between Greater Antillean and continental populations. The same estimators found a lack of differentiation between populations in Mexico and Costa Rica, suggesting deep matrilineal phylogenetic divergences between Central America populations and their conspecifics in the Insular Caribbean.

We report four newly discovered haplotypes: three in Venezuela (Ca11, Ca12, and Ca13) and one in the Everglades (Ca10). Our results show that Venezuelan haplotypes are considerably less distant from Central and North American haplotypes than Greater Antillean haplotypes (Fig 1B). Notably, the addition of South American sampling locations into haplotype designation was useful to improve the understanding of C. acutus intraspecific population genealogies.

A study by Bloor et al. [77], using mtDNA cytochrome b and cytochrome oxidase I gene sequences in captive C. acutus from Colombia, revealed two distinct lineages: one closely related to Central American haplotypes and a second one unique to Colombian C. acutus. Our findings strengthen such evidence and highlight the importance of future research to better understand haplotype relationships within South American C. acutus.

In addition, all but one sample from the Greater Antilles share the same haplotype (β). Haplotype β is closer to Cuban crocodile haplotype α than to any other haplotypes from the Americas [7,32,36]; mitochondrial capture may have occurred during an ancient hybridization event between Greater Antillean C. acutus and C. rhombifer [23]. An alternative scenario postulates multiple colonization events to Cuba by C. rhombifer and later by C. acutus, with sustained periods of hybridization and dispersion across the island [32]. In fact, hybridization between Cuban and American crocodiles in Cuba seems to have taken place both historically and in recent times [7]. Similarly, ancestral and present hybridization has occurred between the American and Morelet’s crocodiles in the Yucatan Peninsula [9,33]. Regardless of potential evolutionary explanations for haplotype origins, the geographic distribution of C. acutus haplotypes suggests that strong genetic structuring shaped haplotype identities, followed by differences in haplotype frequencies.

The new haplotype (Ca10) found in the Everglades seems to be closer to South and Central American haplotypes. In addition, we found two individuals in the Everglades with haplotype β. Rodriguez et al. [23] reported the presence of several haplotypes from Latin America and the Caribbean within Southern Florida and attributed genetic structuring as a result of the admixture of local haplotypes with those of foreign and captive American crocodiles. In addition, the authors reported only one haplotype (Ca1) in core nesting areas in the Everglades. In this context, the two additional haplotypes present in our samples from the Everglades (β and Ca10) are most likely introduced or admixed individuals, suggesting that individuals with foreign haplotypes are being released and/or moved from areas nearby, where non-local haplotypes have been previously reported.

Our analysis of microsatellite data further uncovered patterns of genetic subdivision between populations of the Greater Antilles, Belize, and the United States with at least three populations. It is more likely that four populations occur within these sampling locations, as found by Geneland and the DAPC, as the Evanno method may underestimate K when a hierarchical structure is present. Spatial explicit models applied to other studies have detected biologically meaningful clusters where STRUCTURE failed to detect any population subdivision [78,79] and might be more accurate for populations exhibiting some degree of isolation by distance [80]. Overall, we found sound evidence supporting that Cuban and Jamaican C. acutus represent two distinct populations. The presence of admixture in Jamaica and Cuba could be explained by natural or human-mediated migration of individuals between localities.

Conclusion

Although sampling was limited and extremely challenging as a result of the current political and administrative environment where C. acutus is distributed, our study provides a thorough comparison among populations across the species range and robust evidence of genetic differentiation among populations of the Greater Antilles. With this evidence and that of previous studies [36], we propose the incorporation of independent conservation management units for Cuban and Jamaican C. acutus. Our uneven allocation of sampling effort towards Portland Bright Protected Area versus the area around Black River Lower Morass prevented a more even sample across localities in Jamaica. We strongly suggest that future studies expand the sampling area to further clarify potential patterns of subtle structuring between American crocodiles at these two locations in Jamaica.

Information at the regional-scale is crucial for regional planning and conservation of the species. Importantly, in the context of increasing coastal development throughout the range of the American crocodile, with habitat loss and fragmentation jeopardizing local populations [24,81,82], having readily-available genetic information at the population level will be critical to inform country and site-specific management plans, aid ex-situ conservation efforts, and support the implementation of reintroduction programs.

Future research should also incorporate comparative analysis for other sampling localities in South America and Central America. A broader assessment of C. acutus will further aid conservation efforts and population management decisions to achieve an effective range-wide conservation strategy for the species. Additionally, future studies should aim at using whole-genome sequencing to resolve population structuring further and to understand the phylogenetic evolution of the American crocodile and its relationships with other crocodilian species in the Americas and the Insular Caribbean. Finally, as research looks deeper into C. acutus lineages across its range and suggests the occurrence of potential cryptic species as a result of ancient hybridization and other evolutionary processes [33,36], our study aims at providing population-level data to inform current management and conservation as we continue our debate into C. acutus species designations.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Table of microsatellite genotypes for C. acutus samples collected in this study.

https://doi.org/10.1371/journal.pone.0235288.s001

(XLS)

Acknowledgments

We are grateful to all the colleagues across sampling localities who decided to join this regional effort. We are most grateful to all field staff involved in sample collection, which constituted the first and most critical step for our research project. This project was carried out in collaboration with the Wildlife Conservation Society, the American Museum of Natural History, the Fort Lauderdale Research and Education Center at the University of Florida, the Department of Ecology, Evolution and Environmental Biology at Columbia University, Facultad de Biología de la Universidad de la Habana, the Department of Life Sciences at the University of the West Indies, Museo de Historia Natural La Salle in Venezuela, and the Museo Nacional de Ciencias Naturales (CSIC) in Madrid, Spain. We are especially grateful to our co-authors, the late John Thorbjarnarson and Rafael Crespo, who remain in our hearts and guide our passion for conserving crocodiles in Cuba and around the world.

Citation: Rossi NA, Menchaca-Rodriguez A, Antelo R, Wilson B, McLaren K, Mazzotti F, et al. (2020) High levels of population genetic differentiation in the American crocodile (Crocodylus acutus). PLoS ONE 15(7): e0235288. https://doi.org/10.1371/journal.pone.0235288

References

1. Thorbjarnarson J. American Crocodile Crocodylus acutus. Crocodiles Status Survey and Conservation Action Plan. 3rd ed. Manolis S.C. and Stevenson C., editor. Crocodile Specialist Group; 2010. 46–53 p.

2. Thorbjarnarson JB. Ecology of the American Crocodile, Crocodylus acutus. In: Crocodiles, their Ecology, Management and Conservation. Switzerland: IUCN Publications New Series Gland; 1989. p. 228–58.

3. Thorbjarnarson JB. The status and ecology of the American crocodile in Haiti. Bull Florida State Museum Biol Sci. 1988;33:1–86.

4. Schubert A G. Santana. Conservation of the American crocodile (Crocodylus acutus) in the Dominican Republic. Contrib to Herpetol. 1996;12:425–33.

5. Schubert A. Lagos Enriquillo y Azuéi: Donde la Naturaleza siempre tiene una Sorpresa. Editorial Académica Española; 2012. 88 p.

6. Ortega-León AM, Santos-Morales AH, Zamora-Abrego JG, Pérez-Mendoza HA. Analysis of the population dynamics of the endangered American crocodile, Crocodylus acutus in Paramillo National Natural Park. Mar Freshw Res [Internet]. 2020; Available from: https://doi.org/10.1071/MF19026

7. Milián-García Y, Venegas-Anaya M, Frias-Soler R, Crawford AJ, Ramos-Targarona R, Rodríguez-Soberón R, et al. Evolutionary history of Cuban crocodiles Crocodylus rhombifer and Crocodylus acutus inferred from multilocus markers. J Exp Zool Part A Ecol Genet Physiol. 2011;315 A:358–75.

8. MacHkour-M’Rabet S, Hénaut Y, Charruau P, Gevrey M, Winterton P, Legal L. Between introgression events and fragmentation, islands are the last refuge for the American crocodile in Caribbean Mexico. Mar Biol. 2009;156:1321–33.

9. Pacheco-Sierra G, Gompert Z, Domínguez-Laso J, Vázquez-Domínguez E. Genetic and morphological evidence of a geographically widespread hybrid zone between two crocodile species, Crocodylus acutus and Crocodylus moreletii. Mol Ecol. 2016;25(14):3484–98. pmid:27164458

10. Platt SG, Thorbjarnarson JB. Nesting Ecology of the American Crocodile in the coastal zone of Belize. Copeia. 2006;2000:869–73.

11. Barros T, Jiménez-Oraá M, Heredia H, Seijas A. Artificial incubation of wild-collected eggs of American and Orinoco crocodiles (Crocodylus acutus and C. intermedius), Guárico and Zulia, Venezuela. Conserv Evid. 2010;7:111–5.

12. Cedeño-Vázquez JR, Rodriguez D, Calmé S, Ross JP, Densmore LD, Thorbjarnarson AJB. Hybridization between Crocodylus acutus and Crocodylus moreletii in the Yucatan Peninsula: I. Evidence from mitochondrial DNA and morphology. J Exp Zool Part A Ecol Genet Physiol. 2008;309A:661–73.

13. Rodriguez D, Cedeño-Vázquez JR, Forstner MRJ, Densmore LD. Hybridization between Crocodylus acutus and Crocodylus moreletii in the Yucatan Peninsula: II. Evidence from microsatellites. J Exp Zool Part A Ecol Genet Physiol. 2008;309A:674–86.

14. Cissell JR, Steinberg MK. Human Landscape Modification in Placencia, Stann Creek District, Belize: Possible Implications for Crocodile Hybridization. J Lat Am Geogr. 2020;19(2):218–42.

15. Mazzotti FJ, Dunson WA. Adaptations of Crocodylus acutus and Alligator for life in saline water. Comp Biochem Physiol—Part A Physiol. 1984;79:641–6.

16. Wheatley P V., Peckham H, Newsome SD, Koch PL. Estimating marine resource use by the American crocodile Crocodylus acutus in southern Florida, USA. Mar Ecol Prog Ser. 2012;447:211–29.

17. Black J. Assessment of crocodile abundance and seasonal effects of salinity on distribution using both boat based and aerial drone surveys. Master Thesis. Purdue University Graduate School; 2019.

18. Thorbjarnarson J. Crocodile tears and skins: International trade, economic constraints, and limits to the sustainable use of crocodilians. Conserv Biol. 1999;13:465–70.

19. Thorbjarnarson J, Mazzotti F, Sanderson E, Buitrago F, Lazcano M, Minkowski K, et al. Regional habitat conservation priorities for the American crocodile. Biol Conserv. 2006;128:25–36.

20. Kushlan JA. Conservation and management of the American crocodile. Environ Manage. 1988;12:777–90.

21. Platt SG, Thorbjarnarson JB. Status and conservation of the American crocodile, Crocodylus acutus, in Belize. Biol Conserv. 2000;96:13–20.

22. García-Grajales J, Buenrostro-Silva A. Assessment of human–crocodile conflict in Mexico: patterns, trends and hotspots areas. Mar Freshw Res. 2019;70(5):708–20.

23. Rodriguez D, Forstner MRJJ, Moler PE, Wasilewski JA, Cherkiss MS, Densmore LD. Effect of human-mediated migration and hybridization on the recovery of the American crocodile in Florida (USA). Conserv Genet. 2011;12:449–59.

24. Cox K, Brumund B. Crocodile Smiles: How listed species cope in an age of urbanization. Nat Resour Environ. 2018;32(3):44–7.

25. Tellez M, Boucher M. The lessons of history and the future of American Crocodile conservation in Belize. Herpetol Rev. 2018;49(3):492–8.

26. Ponce-Campos, P., Thorbjarnarson, J. & Velasco A (IUCN SCSG. Crocodylus acutus. IUCN Red List Threat Species. 2012;e.T5659A30(http://dx.doi.org/10.2305/IUCN.UK.2012.RLTS.T5659A3043244.en):Downloaded on 26 July 2019.

27. LeBuff C. Historical review of American crocodiles (Crocodylus acutus) along the Florid Gulf Coast. History. 2016;54(3):50–7.

28. Todesco M, Pascual MA, Owens GL, Ostevik KL, Moyers BT, Hübner S, et al. Hybridization and extinction. Evol Appl. 2016;9(7):892–908. pmid:27468307

29. UNEP (2020). The Species+ Website. Nairobi, Kenya. Compiled by UNEP-WCMC, Cambridge, UK. Available at: www.speciesplus.net. [Accessed (04/29/2020)].

30. Palsbøll PJ, Bérubé M, Allendorf FW. Identification of management units using population genetic data. Trends Ecol Evol. 2007;22(1):11–6. pmid:16982114

31. Sanderson EW, Redford KH, Chetkiewicz CLB, Medellin RA, Rabinowitz AR, Robinson JG, et al. Planning to save a species: The jaguar as a model. Conserv Biol. 2002;16:58–72.

32. Milián-García Y, Ramos-Targarona R, Pérez-Fleitas E, Sosa-Rodríguez G, Guerra-Manchena L, Alonso-Tabet M, et al. Genetic evidence of hybridization between the critically endangered Cuban crocodile and the American crocodile: Implications for population history and in situ/ex situ conservation. Heredity (Edinb). 2015;114:272–280.

33. Pacheco-Sierra G, Vázquez-Domínguez E, Pérez-Alquicira J, Suárez-Atilano M, Domínguez-Laso J. Ancestral hybridization yields evolutionary distinct hybrids lineages and species boundaries in crocodiles, posing unique conservation conundrums. Front Ecol Evol. 2018;6(138).

34. Hekkala ER, Amato G, DeSalle R, Blum MJ. Molecular assessment of population differentiation and individual assignment potential of Nile crocodile (Crocodylus niloticus) populations. Conserv Genet. 2010;11:1435–43.

35. Milián-García Y, Castellanos-Labarcena J, Russello MA, Amato G. Mitogenomic investigation reveals a cryptic lineage of Crocodylus in Cuba. Bull Mar Sci. 2018;94(2):329–43.

36. Milián-García Y, Russello MA, Castellanos-Labarcena J, Cichon M, Kumar V, Espinosa G, et al. Genetic evidence supports a distinct lineage of American crocodile (Crocodylus acutus) in the Greater Antilles. PeerJ. 2018;6:e5836. pmid:30473930

37. Mauger LA, Velez E, Cherkiss MS, Brien ML, Mazzotti FJ, Spotila JR. Conservation genetics of American crocodile, Crocodylus acutus, populations in Pacific Costa Rica. Nat Conserv. 2017;17:1.

38. Moritz C. Applications of mitochondrial DNA analysis in conservation: A critical review. Mol Ecol. 1994;3:401–11.

39. Schlötterer C, Pemberton J. The use of microsatellites for genetic analysis of natural populations—a critical review. In: Molecular Approaches to Ecology and Evolution. 1998. p. 71–86.

40. Mable BK. Conservation of adaptive potential and functional diversity: integrating old and new approaches. Conserv Genet. 2019;20(1):89–100.

41. Clarke CR, Karl SA, Horn RL, Bernard AM, Lea JS, Hazin FH, et al. Global mitochondrial DNA phylogeography and population structure of the silky shark, Carcharhinus falciformis. Mar Biol. 2015;162:945–55.

42. Goldstein DB, Roemer GW, Smith DA, Reich DE, Bergman A, Wayne RK. The use of microsatellite variation to infer population structure and demographic history in a natural model system. Genetics. 1999;

43. Balloux F, Lugon-Moulin N. The estimation of population differentiation with microsatellite markers. Mol Ecol. 2002;11:155–65. pmid:11856418

44. Allendorf FW. Genetics and the conservation of natural populations: allozymes to genomes. Mol Ecol. 2017;26(2):420–30. pmid:27933683

45. IUCN CSG Crocodilian Capacity Building Manual. Available at: http://www.iucncsg.org/pages/Crocodilian-Capacity-Building-Manual-Home.html.

46. Weaver JP, Rodriguez D, Venegas-Anaya M, Cedeño-VÁZQUEZ JR, Forstner MRJJ, Densmore LD, et al. Genetic characterization of captive Cuban crocodiles (Crocodylus rhombifer) and evidence of hybridization with the American crocodile (Crocodylus acutus). J Exp Zool Part A Ecol Genet Physiol. 2008;309(10):649–60.

47. Fitzsimmons NN, Tanksley S, Forstner MRJ, Louis EE, Daglish R, Gratten J, et al. Microsatellite markers for Crocodylus: new genetic tools for population genetics, mating system studies and forensics. In: Crocodilian Biology and Evolution. Surrey Beauty and Sons; 2001. p. 51–7.

48. Ray DA, Dever JA, Platt SG, Rainwater TR, Finger AG, McMurry ST, et al. Low levels of nucleotide diversity in Crocodylus moreletii and evidence of hybridization with C. acutus. Conserv Genet. 2004;5:449–62.

49. Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2. pmid:19346325

50. Nei M. Molecular Evolutionary Genetics. New York: Columbia University Press; 1987. 512 p.

51. Nei M. Molecular population genetics and evolution. Amsterdam, Oxford: North-Holland Publishing Company; 1975. 288 p.

52. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20. pmid:7463489

53. Tajima F. Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983;105:437–60. pmid:6628982

54. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7. pmid:21565059

55. Excoffier L, Laval G, Schneider S. Arlequin: An integrated software package for population genetics data analysis. Evol Bioinforma. 2005;1:47–50.

56. Posada D, Crandall KA. Intraspecific gene genealogies: Trees grafting into networks. Trends Ecol Evol. 2001;16:37–45. pmid:11146143

57. Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48. pmid:10331250

58. Leigh JW, Bryant D. POPART: Full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6(9):1110–6.

59. Rousset F. GENEPOP’007: A complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008;8(1):103–6. pmid:21585727

60. Rice WR. Analyzing tables of statistical tests. Evolution (N Y). 1989;43:223–5.

61. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.

62. Weir BS, Cockerham CC. Estimating F-Statistics for the Analysis of Population Structure. Evolution (N Y). 1984;38(6):1358–70.

63. Guillot G, Mortier F, Estoup A. GENELAND: A computer package for landscape genetics. Mol Ecol Notes. 2005;5(3):712–5.

64. Guillot G, Santos F, Estoup A. Analysing georeferenced population genetics data with Geneland: A new algorithm to deal with null alleles and a friendly graphical user interface. Bioinformatics. 2008;24:1406–7. pmid:18413327

65. Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2015.

66. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59. pmid:10835412

67. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164:1567–87. pmid:12930761

68. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: Dominant markers and null alleles. Mol Ecol Notes. 2007;7:574–8. pmid:18784791

69. Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9(5):1322–32. pmid:21564903

70. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Molecular Ecology. 2005 Jul;14:2611–20. pmid:15969739

71. Earl DA, vonHoldt BM. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.

72. Jombart T. Adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5. pmid:18397895

73. Reich D, Price AL, Patterson N. Principal component analysis of genetic data. Nat Genet. 2008;40(5):1–2.

74. Cedeño-Vázquez JR, Rodriguez D, Calmé S, Ross JP, Densmore LD, Thorbjarnarson AJB, et al. Hybridization between Crocodylus acutus and Crocodylus moreletii in the Yucatan Peninsula: I. Evidence from mitochondrial DNA and morphology. Journal of experimental zoology Part A, Ecological genetics and physiology. 2008 Dec 1;309(10):661–73. pmid:18626922

75. Chapuis MP, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24:621–631. pmid:17150975

76. Carlsson J. Effects of microsatellite null alleles on assignment testing. J Hered. 2008;99:616–623. pmid:18535000

77. Bloor P, Ibáñez C, Viloria-Lagares TA. Mitochondrial DNA analysis reveals hidden genetic diversity in captive populations of the threatened American crocodile (Crocodylus acutus) in Colombia. Ecol Evol. 2015;5(1):130–40. pmid:25628870

78. Lavandero B, Miranda M, Ramrez CC, Fuentes-Contreras E. Landscape composition modulates population genetic structure of Eriosoma lanigerum (Hausmann) on Malus domestica Borkh in central Chile. Bull Entomol Res. 2009;99:97–105. pmid:18947448

79. Barr KR, Lindsay DL, Athrey G, Lance RF, Hayden TJ, Tweddale SA, et al. Population structure in an endangered songbird: Maintenance of genetic differentiation despite high vagility and significant population recovery. Mol Ecol. 2008;17:3628–39. pmid:18643883

80. Renner SC, Suarez-Rubio M, Wiesner KR, Drögemüller C, Gockel S, Kalko EKV, et al. Using multiple landscape genetic approaches to test the validity of genetic clusters in a species characterized by an isolation-by-distance pattern. Biol J Linn Soc. 2016;118:292–303.

81. Somaweera R, Brien ML, Platt SG, Manolis C, Webber BL. Direct and indirect interactions with vegetation shape crocodylian ecology at multiple scales. Freshw Biol. 2019;64(2):257–68.

82. Fitzgerald LA, Walkup D, Chyn K, Buchholtz E, Angeli N, Parker M. The future for reptiles: advances and challenges in the Anthropocene. Encycl Anthr DellaSala, D, Goldstein M (Eds) Elsevier Sci Ltd, Oxford, UK. 2018;163–74.