Genetic analysis of an endemic archipelagic lizard reveals sympatric cryptic lineages and taxonomic discordance

The importance of genetic data in biodiversity conservation is well established, and knowledge of standing genetic variation within and between populations is important for designing conservation strategies. We investigated partitioning of genetic diversity in an endemic lizard (Leiocephalus psammodromus) distributed in the Turks and Caicos archipelago using mtDNA and AFLP data from 259 individuals sampled across 13 islands. Current taxonomy identifies six or more subspecies of L. psammodromus within the archipelago, several of which have undergone recent drastic reductions in range due to extirpation. However, our results indicate the presence of two independent lineages, one on each of the Turks and Caicos banks, and a third sympatric cryptic lineage on both banks. These lineages do not correspond to current taxonomy and alter our understanding of diversity and conservation of this species. Gross morphological data (mass and snout-vent length) indicate some variation in female size among lineages, indicating the possibility of cryptic morphological variation. Instead of initiating separate conservation measures for nominate subspecies, we recommend a more thorough investigation of the morphology and genetics of this group and a more inclusive conservation program. Our surprising results indicate that other endemic squamates in the Bahamas Archipelago might also exhibit sympatric cryptic diversity that does not correspond to current taxonomic understanding and could have significant impacts on our approach to conservation in this region.     

Population genetic diversity and divergence of the halobiotic herb Limonium sinense estimated by AFLP and ISSR, and implications for conservation

Limonium sinense is a halobiotic herb endemic to China that has been traditionally used for hundreds of years for its good restorative function. Genetic variation and population structure of this species were investigated by using amplified fragment length polymorphisms (AFLPs) and inter simple sequence repeats (ISSRs). A high level of genetic diversity was detected [AFLP: H _E = 0.284, percentage of polymorphic loci (PPL) = 92.68 %; ISSR: H _E = 0.257, PPL = 85.71 %] at the species level with POPGENE. Based on analysis of molecular variation (AMOVA), the among-population component accounted for 29.03 % (AFLP) and 28.81 % (ISSR) of the genetic variation, indicating that most of the genetic variation was between individuals within populations. The Shannon diversity index (I) was higher for AFLP (0.432) than for ISSR (0.395). Five main clusters were shown in the unweighted pair-group method with arithmetic mean (UPGMA) dendrogram created using TFPGA, consistent with the result of principal coordinate analysis using NTSYS. In situ conservation is advocated first. Keeping a stable environment for this halobiotic herb is necessary. For ex situ conservation, it is important to establish a germplasm bank. AFLP and ISSR markers were proved to be efficient tools in assessing the genetic variation among populations of L. sinense. The patterns of variation appeared to be consistent for these two marker systems, and they can be used for management of genetic structure, protection of the halobiotic plant, and conservation of germplasm.     

DNA fingerprinting of Leonurus cardiaca L. germplasm in Iran using amplified fragment length polymorphism and inter-retrotransposon amplified polymorphism

In the present study, the extent of inter and intra-population genetic variation was evaluated in Leonurus cardiaca accessions naturally growing in Iran by AFLP and IRAP markers. The fingerprints corresponding to AFLP and IRAP markers revealed high levels of heterozygosity, indicating that L. cardiaca is predominantly an out-crossing species. The average percentage polymorphism was detected as 58% and 90.8% on utilizing AFLP and IRAP data, respectively. Gene diversity values within populations varied 0.14 to 0.20 for AFLP and 0.12 to 0.21 for IRAP. The overall levels of genetic variation present in the L. cardiaca germplasm in Iran were finally determined by combining the AFLP and IRAP datasets to ensure wide genome coverage. The phenogram depicted that the accessions of Dargaz population were genetically distinct from other populations. Based on AFLP and IRAP analysis, it is concluded that L. cardiaca maintains high levels of genetic variation at inter and intra-population level.     

Genetic diversity in traditional Ethiopian highland maize accessions assessed by AFLP markers and morphological traits

In the highland regions of Ethiopia the heterogeneity of the land, the climate, and soil favors the presence of a large number of landraces. We analyzed a representative sample of 62 traditional Ethiopian highland maize accessions, using amplified fragment length polymorphism (AFLP^®) markers and morphological traits with the aim to group the accessions based on the their genetic profiles and morphological traits, to study agroecological variation and to assess the level of correlation between phenotypic and genetic distances. Eight EcoRI/MseI primer combinations and 15 morphological traits were used. The accessions varied significantly for all of the measured morphological traits. Of a total of 650 AFLP markers that were scored, 89.5% were polymorphic. Pair-wise genetic distance estimates based on AFLP data revealed dissimilarity coefficients ranging from 0.32 to 0.69 (mean of 0.57). Cluster analysis of the AFLP data grouped most accessions collected from the Northern highlands into one major cluster. It, however, failed to separate the Western and Southern accessions into different clusters. Regardless of the large variation in environmental conditions between agroecologies, only 9% of the total genetic variation was found between agroecologies, whereas 91% was found within agroecologies in Ethiopia. This finding may be explained by long distance seed exchange, continuous seed introduction and gene flow between agroecologies. The relationship between morphological and AFLP-based distances was significant and positive. Based on this study, three groups of highland accessions, with distinctive genetic profiles and morphological traits were identified. This information will be useful for further collections and conservation of the unique diversity included in the highland maize landraces of Ethiopia.     

Genetic diversity analysis in sorghum germplasm as estimated by AFLP, SSR and morpho-agronomical markers

A comparison of the different methods of the estimation of genetic diversity is important to evaluate their utility as a tool in germplasm conservation and plant breeding. Amplified fragment length polymorphism (AFLP), microsatellites or SSR and morphological traits markers were used to evaluate 45 sorghum germplasm for genetic diversity assessment and discrimination power. The mean polymorphism information content (PIC) values were 0.65 (AFLPs) and 0.46 (SSRs). The average pairwise genetic distance estimates were 0.57 (morphological traits), 0.62 (AFLPs) and 0.60 (SSRs) markers data sets. The Shannon diversity index was higher for morphological traits (0.678) than AFLP (0.487) and SSR (0.539). The correlation coefficients obtained by the Mantel matrix correspondence test, which was used to compare the cophenetic matrices for the different markers, showed that estimated values of genetic relationship given for AFLP and SSR markers, as well as for morphological and SSR markers were significantly related (p <0.001). However, morphological and AFLP data showed non-significant correlation (p >0.05). Both data sets from AFLP and SSR allowed all accessions to be uniquely identified; two accessions could not be distinguished by the morphological data. In summary, AFLP and SSR markers proved to be efficient tools in assessing the genetic variability among sorghum genotypes. The patterns of variation appeared to be consistent for the three marker systems, and they can be used for designing breeding programmes, conservation of germplasm and management of sorghum genetic resources.     

Delineating species along shifting shorelines: <Emphasis Type="Italic">Tropheus</Emphasis> (Teleostei,Cichlidae) from the southern subbasin of Lake Tanganyika

Background

Species delineation is particularly challenging in taxa with substantial intra-specific variation. In systematic studies of fishes, meristics and linear measurements that describe shape are often used to delineate species. Yet, little is known about the taxonomic value of these two types of morphological characteristics. Here, we used Tropheus (Teleostei, Cichlidae) from the southern subbasin of Lake Tanganyika to test which of these types of characters best matched genetic lineages that could represent species in this group of stenotypic rock-dwelling cichlids. We further investigated intra-population variation in morphology. By linking this to a proxy of a population’s age, we could assess the evolutionary stability of different kinds of morphological markers.

Results

Morphological data was collected from 570 specimens originating from 86 localities. An AFLP approach revealed the presence of five lineages in the southern subbasin: T. moorii, T. brichardi, T. sp. ‘maculatus’, T. sp. ‘Mpimbwe’ and T. sp. ‘red’, which we consider to represent distinct species. Although both types of morphological data supported this classification, a comparison of P_ST-values that describe inter-population morphological differentiation, revealed a better correspondence between the taxon delineation based on AFLP data and the patterns revealed by an analysis of meristics than between the AFLP-based taxon delineation and the patterns revealed by an analysis of shape. However, classifying southern populations of Tropheus was inherently difficult as they contained a large amount of clinal variation, both in genetic and in morphological data, and both within and among species. A scenario is put forward to explain the current-day distribution of the species and colour varieties and the observed clinal variation across the subbasin’s shoreline. Additionally, we observed that variation in shape was larger in populations from shallow shores whereas populations from steep shores were more variable in meristics. This difference is explained in terms of the different timescales at which small and large scale lake level fluctuations affected populations of littoral cichlids at steep and shallow shores.

Conclusions

Our results showed meristics to be more evolutionary stable, and of higher taxonomic value for species delimitation in Tropheus, than linear measurements that describe shape. These results should be taken into account when interpreting morphological differences between populations of highly stenotypic species, such as littoral cichlids from the Great East African Lakes.

     

High genetic diversity and differentiation of an extremely narrowly distributed and critically endangered decaploid rose (<Emphasis Type="Italic">Rosa praelucens</Emphasis>): implications for its conservation

Rosa praelucens is a critically endangered decaploid alpine rose with an extremely narrow geographic distribution in Northwestern Yunnan, China. We sampled almost all the extant individuals (527 individuals in 31 natural locations and 56 individuals preserved in three local living collections) to assess the genetic variation and to probe the genetic connectivity among the individuals and populations based on three cpDNA intergenic spacers and six fluorescent amplified fragment length polymorphism (AFLP) markers. The morphological traits from seven populations were also measured. R. praelucens exhibited high levels of morphological variation, genetic diversity, and differentiation. The extant individuals were clustered into eight groups in neighbor-net networks, and subsequent Bayesian analysis assigned them into three larger gene pools, both in accordance with their morphological traits, especially flower color. The living collections embraced two private cpDNA haplotypes and included three out of the species’ total eight AFLP genotypes. Rhizome clonal growth, decaploid, and mixed breeding system may largely contribute to high genetic diversity and differentiation in R. praelucens. We concluded that the endangered status of R. praelucens may mainly be due to habitat fragmentation and loss and inherent reproductive difficulties, rather than low genetic diversity. The populations contributing higher cpDNA genetic diversity, representing more AFLP genotypes, and encompassing private cpDNA haplotypes should be given conservation priority by creating plant-micro reserves. The living collections should also be targeted for further ex situ conservation, population recovery, and reintroduction of R. praelucens plants.     

AFLP and ISSR analysis reveals high genetic variation and inter-population differentiation in fragmented populations of the endangered Litsea szemaois (Lauraceae) from Southwest China

Litsea szemaois (Lauraceae) is an endemic and endangered species from the tropical rain forests of Xishuangbanna, southern Yunnan, SW China, but habitat fragmentation, especially exacerbated by rubber planting, has caused a decline in population size of the species. AFLP and ISSR were used to investigate the genetic diversity and population structure of eight populations from across its known distribution. Three AFLP and ten ISSR primer combinations produced a total of 203 and 77 unambiguous and repeatable bands respectively, of which 164 (80.8%) and 67 (87.0%) were polymorphic for the two markers. These two markers showed that Litsea szemaois exhibits comparatively high genetic diversity at species level (heterozygosity (hs) = 0.2109) relative to some other Lauraceae. Most of the genetic variation was partitioned within populations, but genetic differentiation between populations was significant and relatively high (Φ _st = 0.2420, θ^B = 0.1986) compared with other tropical plants. The genetic characteristics of L. szemaois may be related to its outbreeding system, insect pollination and fragmented distribution. Because L. szemaois is dioecious and slow to mature, ex situ conservation across its genetic diversity is unlikely to succeed, although seedlings grow well under cultivation. Thus, in situ conservation is very important for this endangered species, especially as only 133 adult individuals are known in the wild. In particular, the Nabanhe and Mandian populations should be given a high conservation priority due to their higher genetic diversity, larger population size and better field condition, but wider sampling is required across all populations to determine additional areas with significant genetic conservation value.     

Phylogeography of palo santo trees (Bursera graveolens and Bursera malacophylla; Burseraceae) in the Galápagos archipelago

Two species of palo santo trees are distributed in the Galápagos archipelago, the native Bursera graveolens and the endemic Bursera malacophylla (Burseraceae). However, a zone of individuals morphologically intermediate between the two exists on northern Santa Cruz Island and south‐eastern Santiago Island, suggesting that they may not be reproductively isolated. Here we review the species' distributional and morphological differences and test for evidence of hybridization between the two species using DNA sequence and AFLP data. We find that the species lack distinguishing synapomorphies across the five nuclear and plastid regions examined. Population assignment tests and population genetic analyses of AFLP data indicate that genetically similar palo santo individuals (N = 87), including putative hybrids, partition into two genealogical groups that do not uniformly correspond to island‐ or taxon‐based membership. Furthermore, genotypic admixture levels among morphologically intermediate individuals do not indicate widespread hybridization. Thus, we recommend recognizing the endemic palo santo taxon as B. graveolens subspecies malacophylla (B.L.Rob.) A. Weeks & Tye comb. & stat. nov. in light of its close genetic relationship to B. graveolens subsp. graveolens and its distinctive morphology and distribution. Future research should quantify phenotypic variation in palo santo populations as another means for understanding the basis of morphological differences between the subspecies. © 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 161 , 396–410.     

Genetic structure of the <Emphasis Type="Italic">Anthyllis vulneraria</Emphasis> L. s. l. species complex in Estonia based on AFLPs

Anthyllis vulneraria L. (Fabaceae) s. lato includes many cryptic taxa, ranging from 25 to 60 subspecies according to different authors. The delimitation of intraspecific taxa of A. vulneraria s. lato has always been complicated and inconsistent. Different data sets (multivariate analyses of morphological variation, allozymes, chloroplast SSRs and ITS) have not resolved the existing problem with distinguishing some subspecies. We used the amplified fragment length polymorphism (AFLP) analysis to describe the differentiation in this species complex and to characterize variation on a geographic scale. Some correlation was found between genetic variability and geographic distribution (western-eastern directional variation), but AFLP data analysis did not reveal clear intraspecific structure of the seven analysed taxa. The analysed specimens did not comprise groups correlated with the subspecies.     

Local genetic structure of a montane herb among isolated grassland patches: implications for the preservation of genetic diversity under climate change

Habitat loss, fragmentation of meadow patches, and global climate change (GCC) threaten plant communities of montane grasslands. We analyzed the genetic structure of the montane herb Geranium sylvaticum L. on a local scale in order to understand the effects of habitat fragmentation and potential GCC impacts on genetic diversity and differentiation. Amplified fragment length polymorphism (AFLP) fingerprinting and cpDNA sequencing was performed for 295 individuals of 15 G. sylvaticum populations spanning the entire distribution range of the species in the Taunus mountain range in Germany. We found patterns of substantial genetic differentiation among populations using 150 polymorphic AFLP markers (mean F _ST = 0.105), but no variation in 896 bp of plastid DNA sequences. While populations in the center of their local distribution range were genetically diverse and less differentiated, higher F _ST values and reduced genetic variability was revealed for the populations at the low-altitudinal distribution margins. Projections of GCC effects on the distribution of G. sylvaticum in 2050 showed that GCC will likely lead to the extinction of most edge populations. To maintain regional genetic diversity, conservation efforts should focus on the diverse high-altitude populations, although a potential loss of unique variations in genetically differentiated peripheral populations could lower the overall genetic diversity and potentially the long-term viability in the study region. This study documents the usefulness of fine-scale assessments of genetic population structure in combination with niche modeling to reveal priority regions for the effective long-term conservation of populations and their genetic variation under climate change.     

Norfolk Island Robins are a distinct endangered species: ancient DNA unlocks surprising relationships and phenotypic discordance within the Australo-Pacific Robins

Uncertain taxonomy hinders the effective prioritization of taxa for conservation. This problem is acute for understudied island populations in the southwest Pacific Ocean, which are increasingly threatened by habitat loss, predation and climate change. Here, we offer the first test of taxonomic limits and phylogenetic affinities of the iconic Pacific Robin radiation (Petroica multicolor) in order to prioritize the conservation of its nominotypical subspecies, the endangered Norfolk Island Robin (P. m. multicolor). We integrate phylogenetic analyses of ancient DNA and quantitative measures of plumage and morphometric variation to show that the Norfolk Island Robin should be recognized as a distinct species. Phenotypic and genetic datasets contradict the longstanding treatment of Pacific Robins (including Norfolk Island Robins) and Scarlet Robins (P. boodang) as a single species. Instead, we show that Norfolk Island Robins are deeply divergent from Scarlet Robins and have more genetic similarity to Red-capped Robins (P. goodenovii) than to other Pacific Robins. This finding is unrepresentative of the current taxonomic and conservation status of the Norfolk Island Robin, which we propose should be recognised as an endemic endangered species. Our study clearly shows that in the absence of contemporary tissues, ancient DNA approaches using historical museum specimens can address taxonomic questions that morphological traits are unable to resolve. Further, it highlights the need for similar studies of other threatened Norfolk fauna with uncertain taxonomic status in order to ensure appropriate conservation prioritization.     

Genetic diversity and structure in fragmented populations of the endangered species Ranunculus cabrerensis (Ranunculaceae): implications for conservation

Ranunculus cabrerensis is an endemic and endangered species of the Northwestern Iberian Peninsula. The molecular markers AFLP and ISSR were used to investigate the genetic diversity and population structure of four populations across its known distribution. Fifteen selective primer combinations of AFLP and seventeen ISSR primer combinations produced a total of 2830 and 103 unambiguously repeatable fragments respectively, of which 97.57 and 81.38% were polymorphic for both markers. The genetic diversity of R. cabrerensis at species level was high (H _E = 0.294 by ISSR and H _E = 0.191 by AFLP) and differentiation between sampled locations was also relatively high (G _ST = 0.316 and 0.158 by ISSR and AFLP analysis respectively) compared to other studies of endangered and rare species using the same techniques. The analysis of molecular variance (AMOVA) indicated that the main genetic variation was within sampled locations (73% by AFLP; 52% by ISSR), even though the variation among locations was also significant. Principal Coordinates, NeighborNet and Bayesian analyses revealed a weak but significant relationship between the genetic structures of different populations in R. cabrerensis, with gene flow acting as a homogenizing force that prevents stronger differentiation of populations. Finally, suggestions for conservation strategies to preserve the genetic resources of this species are outlined.     

Ecotypes and genetic structure of Rhinanthus alectorolophus (Orobanchaceae) in southwestern Germany

European annual species of the genus Rhinanthus often exhibit seasonal ecotypic variation, a phenomenon also known from related genera of hemiparasitic Orobanchaceae. Populations with different flowering times exist, correlated with differences in a number of morphological characters. The present study evaluates the correlation of morphological characters and genetic differentiation of populations of Rhinanthus alectorolophus. Thirty-nine populations of three different subspecies from southwestern Germany were sampled. A total of 798 individuals were used for morphological analyses and 187 of these for AFLP analyses. Principal component analysis showed that morphological variation is mostly continuous. In a discriminant analysis based on morphological characters, only 89.7 % of all individuals were correctly assigned to their previously determined subspecies, indicating that subspecies identification is ambiguous for some populations. Using AFLP data and Bayesian assignment analysis, the sampled individuals could be grouped in three genetic clusters which do not correspond to the three subspecies. Instead, the clustering shows a clear geographic pattern and a Mantel test likewise revealed a significant correlation between genetic and geographic distances. Correlations of genetic distances with differences in morphological characters were weak and mostly insignificant. The results indicate that the subspecies of R. alectorolophus do not form discrete entities and that the character combinations distinguishing them are homoplastic.     

Contrasting coarse and fine scale genetic structure among isolated relic populations of Kmeria septentrionalis

Trees of the Magnoliaceae family are of scientific, cultural and socio-economic importance. Kmeria septentrionalis Dandy (Magnoliaceae) is a dioecious tree, found in small, isolated, relic populations in Southern China, and is subject to extensive protection due to its rarity and high economic values. To improve conservation outcomes and in particular, germplasm collection guidelines, information on spatial genetic structure of the species is required. In this study, we investigated the spatial genetic structure and genetic diversity of 161 individuals of K. septentrionalis collected from five natural populations using AFLP molecular markers. Within-population genetic variation was measured, with percentage of polymorphic bands (PPB) ranged from 63% to 87%, while H _S (genetic diversity within population) varied from 0.185 to 0.244 with a mean of 0.215 ± 0.025. Significant genetic differentiations were revealed between pairwise populations, indicating each population existing as an independent evolutionarily significant unit. Mantel test results showed no pattern of isolation-by-distance among populations separated by large distance. Fine scale spatial patterns of genetic variation suggested significant effects of isolation-by-distance within population at distances of 22 m. The results of contrasting genetic structure at coarse and fine scale in K. septentrionalis may indicate restricted pollen flow and seed dispersal at fine scales, and separated evolution in isolated populations over long period of time at coarser scales. Finally, we make several suggestions for improved management practices that may assist in the conservation of this species.     

Conservation genetics of Heritiera littoralis (Sterculiaceae), a threatened mangrove in China,based on AFLP and ISSR markers

AFLP and ISSR markers were used to determine the genetic variations in eight mangrove and non-mangrove populations of Heritiera littoralis (Sterculiaceae), a threatened species in China. Our results showed a moderate to high level of genetic variation in this species (P = 63.69%, H_T = 0.20 for AFLP; P = 76.07%, H_T = 0.22 for ISSR), and a relatively high level of genetic differentiation among populations (G_ST = 0.24 for AFLP and 0.27 for ISSR). Life history traits, long-distance dispersal of floating seeds, and local environments may play important roles in shaping the genetic diversity and genetic structure of this species. Investigation of the plant’s reproductive capacity showed that the natural seed production of H. littoralis was low, as was the germination rate and the transformation rate from juvenile to adult. H. littoralis is seriously threatened and is in urgent need of conservation in China.     

Integrated species delimitation and conservation implications of an endangered weevil Pachyrhynchus sonani (Coleoptera: Curculionidae) in Green and Orchid Islands of Taiwan

Oceanic islands are productive habitats for generating new species and high endemism, which is primarily due to their geographical isolation, smaller population sizes and local adaptation. However, the short divergence times and subtle morphological or ecological divergence of insular organisms may obscure species identity, so the cryptic endemism on islands may be underestimated. The endangered weevil Pachyrhynchus sonani Kôno (Coleoptera: Curculionidae: Entiminae: Pachyrhynchini) is endemic to Green Island and Orchid Island of the Taiwan‐Luzon Archipelago and displays widespread variation in coloration and host range, thus raising questions regarding its species boundaries and degree of cryptic diversity. We tested the species boundaries of P. sonani using an integrated approach that combined morphological (body size and shape, genital shape, coloration and cuticular scale), genetic (four genes and restriction site‐associated DNA sequencing, RAD‐seq) and ecological (host range and distribution) diversity. The results indicated that all the morphological datasets for male P. sonani, except for the colour spectrum, reveal overlapping but statistically significant differences between islands. In contrast, the morphology of the female P. sonani showed minimum divergence between island populations. The populations of P. sonani on the two islands were significantly different in their host ranges, and the genetic clustering and phylogenies of P. sonani established two valid evolutionary species. Integrated species delimitation combining morphological, molecular and ecological characters supported two distinct species of P. sonani from Green Island and Orchid Island. The Green Island population was described as P. jitanasaius sp.n. Chen & Lin, and it is recommended that its threatened conservation status be recognized. Our findings suggest that the inter‐island speciation of endemic organisms inhabiting both islands may be more common than previously thought, and they highlight the possibility that the cryptic diversity of small oceanic islands may still be largely underestimated.     

Erodium maritimum (Geraniaceae), a species with an uneven and fragmented distribution along the Western Mediterranean and European Atlantic coasts,has a weak genetic structure

Erodium maritimum L. is an annual species presenting heterogeneous, sometimes very small, and distant populations, distributed along a discontinuous coastal strip of the European Atlantic and the central and western Mediterranean basin. The aim of this study is to investigate genetic variation and geographic structure changes across its large distribution. Fourteen populations of E. maritimum were studied using AFLP fingerprints, together with their population sizes, reproductive systems and flower visitors. AFLP markers revealed the genetic structure of the species to be weak. Many individuals from one population clustered together with those of other populations, showing a high degree of genetic admixture. Despite having a self‐compatible reproductive system, populations (especially the largest ones) showed high levels of genetic polymorphism, and the majority of genetic variation was contained within populations. The low genetic structure suggests high levels of gene flow, which might be explained through the dispersability of the species’ fruits. Finally, recommendations are provided for management strategies to facilitate the conservation of this endangered species.     

The Role of Geographical and Ecological Factors on Population Divergence of the Neotropical otter <Emphasis Type="Italic">Lontra longicaudis</Emphasis> (Carnivora,Mustelidae)

The geographic distribution of the populations of a species are influenced by the spatial structure of the ecosystems, the environmental factors and the presence of geographic barriers. The Neotropical otter, Lontra longicaudis, is widely distributed throughout the Americas, where a wide range of environmental conditions and geographical features could promote genetic and morphological variation on the three currently recognized subspecies. In this study, we combined phylogeographic, morphometric and environmental niche modelling analyses to examine whether: (1) genetic variation is associated with the presence of barriers to gene flow and/or hydrography; (2) genetic and morphologic variation are associated with environmental variation; and (3) the observed variation in L. longicaudis populations corresponds to the previously defined subspecies. We found strong phylogeographic structure between the northern (L. l. annectens) and the two-southern subspecies (L. l. longicaudis and L. l. enudris), and although shallower, we also detected genetic differentiation between the two South American subspecies. Such genetic differentiation corresponds to the hydrography and to the geographical barriers characteristic of the distributional area of the species. We found a correlation between the shape of the skull and mandible with the environmental variation through the distribution of the species, and we rejected the hypothesis of niche equivalency and similarity between the three identified genetic lineages, suggesting adaptations to different environmental conditions. Our results support that the variation in environmental conditions, in concert with geographical barriers to gene flow and hydrography, have led to population divergence of L. longicaudis along the Neotropics. These results have important taxonomic implications for the species and its conservation.