Molecular and morphological patterns of introgression between two large white-headed gull species in a zone of recent secondary contact

Incomplete reproductive isolation promotes gene flow between diverging taxa. However, any gene encoding for traits involved in the reproductive barriers will be less prone to introgression than neutral markers. Comparing introgression rates among loci is thus informative of the number and functions of loci involved in the reproductive barriers. This study aimed at identifying possible mechanisms of restriction to gene flow across a zone of recent secondary contact between Larus argentatus and Larus cachinnans by comparing introgression patterns for nine microsatellite loci, a fragment of mitochondrial DNA and a set of phenotypic traits. The low linkage disequilibrium between neutral nuclear markers indicated introgression without any barrier to gene flow. However, asymmetric introgression of mitochondrial DNA suggested that interspecific crosses may be more successful in one direction. The introgression rate for phenotypic traits was variable and low compared to neutral molecular markers. This was particularly evident in colouration of bare parts: individuals with intermediate colouration were scarcer in sympatry than expected if the genomes recombined freely. We hypothesized that one of these variables, the orbital ring colour, may play a role in mate choice, acting as an incomplete premating barrier through assortative mating. This study emphasizes that multilocus approaches are useful to discriminate among possible mechanisms responsible for the maintenance of hybrid zones.     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Appearance of morphological novelty in a hybrid zone between two species of land snail

On the small oceanic island of Chichijima, two endemic species of land snails, Mandarina mandarina and M. chichijimana, have discrete distributions separated by a hybrid zone. This study investigates the potential of hybridization as a source of morphological novelty in these snails. Mandarina mandarina possesses a shell with a higher whorl expansion rate and a smaller protoconch than M. chichijimana, relative to shell size. The number of whorls and shell size of M. mandarina do not differ from those of M. chichijimana, because the effect of higher expansion rate on number of whorls and size of the former is compensated for by its smaller protoconch. The whorl expansion rate and protoconch diameter of the individuals from the hybrid populations are intermediate or typical of either of the two species, and their average values show clinal changes along the hybrid zone. However, the hybrid populations include exceptionally high shells with many whorls and flat shells with few whorls, which are never found in the pure populations of either species. In addition, gradual increases in variance in shell height and number of whorls were found from the edges to the center of the hybrid zone. A combination of low expansion rate (typical of M. chichijimana) and a small protoconch (typical of M. mandarina) produces a shell with an extremely large number of whorls because of the geometry of shell coiling. However, the combination of high expansion rate and a large protoconch produces a shell with an extremely small number of whorls. Because of the correlation between the number of whorls and shell height, shells with an exceptional number of whorls possess an extraordinarily high or flat spire. Hybrids can inherit a mosaic of characters that, as they play out during growth, lead to novel adult morphologies. These findings emphasize the importance of hybridization as a source of morphological variation and evolutionary novelty in land snails.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves,Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early‐breeding house sparrows and hybrids and rural late‐breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z‐chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land‐use changes in a mosaic landscape.     

Low Levels of Nucleotide Diversity in Crocodylus moreletii and Evidence of Hybridization with C. acutus

Examinations of both population genetic structure and the processes that lead to such structure in crocodilians have been initiated in several species in response to a call by the IUCN Crocodile Specialist Group. A recent study used microsatellite markers to characterize Morelet's crocodile (Crocodylus moreletii) populations in north-central Belize and presented evidence for isolation by distance. To further investigate this hypothesis, we sequenced a portion of the mitochondrial control region for representative animals after including samples from additional locales in Belize, Guatemala and Mexico. While there is limited evidence of subdivision involving other locales, we found that most of the differentiation among populations of C. moreletiican be attributed to animals collected from a single locale in Belize, Banana Bank Lagoon. Furthermore, mitochondrial DNA sequence analysis showed that animals from this and certain other locales display a haplotype characteristic of the American crocodile, C. acutus, rather than C. moreletii. We interpret this as evidence of hybridization between the two species and comment on how these new data have influenced our interpretation of previous findings. We also find very low levels of nucleotide diversity in C. moreletiihaplotypes and provide evidence for a low rate of substitution in the crocodilian mitochondrial control region. Finally, the conservation implications of these findings are discussed.     

A narrow hybrid zone between two crayfish species from a Mexican cave

In Northern Chiapas (Mexico), two newly discovered species of Procambarus crayfish inhabit a subterranean stream. These species can be morphologically distinguished only by comparing extreme phenotypes (dark, thick-eyed, surface dwelling-like individuals vs light, elongated, microphtalmic, cave dwelling-like individuals). Individuals with intermediate phenotypes co-occur with those exhibiting extreme phenotypes. Crayfish were assayed electrophoretically and individual patterns at 23 gene loci were obtained. Unusually high levels of heterozygosity in both species and a clear discrimination between the two gene pools were revealed. The relationships between individuals were investigated by means of multivariate analysis on individual multilocus genotype profiles. Results showed the occurrence of individuals genetically intermediate between the two major clusters, which shared allozymic variants with both species. Due to the occurrence of alternative alleles in the two gene pools, we could quantify patterns of introgression, which revealed asymmetric gene flow between the two species. Moreover, differential levels of introgression in subsamples within the surface-like species were found: most introgressed individuals came from the inner section of the cave, where the two species were greatly mixed. These results are also discussed in reference to the morphometric results from a companion paper. A possible evolutionary pathway, leading to the situation in this cave, and possibly in neighbouring cave systems, is outlined. The hypothesis of a past history of allopatric divergence from a common ancestor and a subsequent secondary contact between these two Procambarus species is supported by geological studies. Crayfish sympatry and competitive exclusion are also discussed.     

Genetic characterization of captive Cuban crocodiles (Crocodylus rhombifer) and evidence of hybridization with the American crocodile (Crocodylus acutus)

There is a surprising lack of genetic data for the Cuban crocodile (Crocodylus rhombifer), especially given its status as a critically endangered species. Samples from captive individuals were used to genetically characterize this species in comparison with other New World crocodilians. Partial mitochondrial sequence data were generated from cyt-b (843 bp) and the tRNA(Pro)- tRNA(Phe)-D-loop region (442 bp). Phylogenetic analyses were performed by generating maximum parsimony, maximum likelihood, and Bayesian-based topologies. In addition, in an effort to identify species-specific alleles, ten polymorphic microsatellite loci were genotyped. Distance and model-based clustering analyses were performed on microsatellite data, in addition to a model-based assignment of hybrid types. Both mitochondrial and nuclear markers identified two distinct C. rhombifer genetic sub-clades (alpha and beta); and microsatellite analyses revealed that most admixed individuals were F(2) hybrids between C. rhombifer-alpha and the American crocodile (C. acutus). All individuals in the C. rhombifer-beta group were morphologically identified as C. acutus and formed a distinct genetic assemblage. J. Exp. Zool. 309A:649-660, 2008. (c) 2008 Wiley-Liss, Inc.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.     

Human disturbance causes the formation of a hybrid swarm between two naturally sympatric fish species

Most evidence for hybrid swarm formation stemming from anthropogenic habitat disturbance comes from the breakdown of reproductive isolation between incipient species, or introgression between allopatric species following secondary contact. Human impacts on hybridization between divergent species that naturally occur in sympatry have received considerably less attention. Theory predicts that reinforcement should act to preserve reproductive isolation under such circumstances, potentially making reproductive barriers resistant to human habitat alteration. Using 15 microsatellites, we examined hybridization between sympatric populations of alewife (Alosa pseudoharengus) and blueback herring (A. aestivalis) to test whether the frequency of hybridization and pattern of introgression have been impacted by the construction of a dam that isolated formerly anadromous populations of both species in a landlocked freshwater reservoir. The frequency of hybridization and pattern of introgression differed markedly between anadromous and landlocked populations. The rangewide frequency of hybridization among anadromous populations was generally 0–8%, whereas all landlocked individuals were hybrids. Although neutral introgression was observed among anadromous hybrids, directional introgression leading to increased prevalence of alewife genotypes was detected among landlocked hybrids. We demonstrate that habitat alteration can lead to hybrid swarm formation between divergent species that naturally occur sympatrically, and provide empirical evidence that reinforcement does not always sustain reproductive isolation under such circumstances.     

Morphological and molecular evidence for hybridization and introgression in a willow (Salix) hybrid zone

Hybrid zones provide biologists with the opportunity to examine genetic and ecological interactions between differentiated populations. Accurate identification of hybrid genealogies is considered a necessary prerequisite to understanding observed patterns of hybridization-related phenomena. We analysed molecular and morphological data from individuals in a hybrid zone between two species of willows (Salix sericea Marshall and S. eriocephala Michaux) and report the use of randomly amplified polymorphic DNA (RAPD), chloroplast DNA (cpDNA), and ribosomal DNA (rDNA) markers, as well as vegetative morphology and foliar chemistry data to identify individuals in terms of hybrid genealogy and to infer the direction and extent of backcrossing and introgression within the hybrid zone. A novel version of a maximum likelihood estimate approach (developed for this study) was used to calculate hybrid index scores from RAPD marker data; this method produced results similar to those obtained using traditional arithmetic methods. Distribution of rDNA, cpDNA, and chemistry data were examined within the graphical context of RAPD-based hybrid index score histograms and principal component analyses (PCA) on RAPD and morphology data. Seven of the 21 plants classified as S. eriocephala in the field were possible introgressants. Another plant presented an unequivocal example of backcrossed S. sericea chemistry and RAPD markers. Inter- and intraspecific chloroplast diversity found within the hybrid zone suggests both historic introgression (perhaps in a glacial refugium), and contemporary hybridization. Patterns of inheritance and expression within the hybrid zone suggest that morphological characters are often not expressed in a simple additive fashion, and problems associated with both morphological and molecular data are considered.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

Genetic variation in a geographically restricted and two widespread species of South American Nothofagus

The amount and distribution of genetic variation is compared using starch gel electrophoresis among populations of the three evergreen species of South American Nothofagus : the geographically restricted Nothofagus nitida (Phil.) Krasser and the widespread Nothofagus betuloides (Mirb.) Oerst. and Nothofagus dombeyi (Mirb.) Oerst. Ten enzyme systems were resolved coding for fifteen putative genetic loci. N. betuloides and N. dombeyi were more genetically variable than N. nitida ; they had higher total number of alleles, mean number of alleles per locus, percent of polymorphic loci, mean expected heterozygosity, and mean total genetic diversity. These results are in agreement with the hypothesis that widespread species are often associated with historically large and more continuously distributed populations which in turn can maintain higher levels of genetic variation. Conversely, the genetic structure of N. nitida including its low genetic identity with the other two species and the presence of four unique alleles support the hypothesis that N. nitida has been isolated from the other two species for a considerable period of time. Genetic data in concert with the fossil record indicate that these taxa may be members of a pre-Pleistocene flora and consequently their genetic structures reflect a more ancient evolutionary history than previously hypothesized.     

Differential introgression of a female competitive trait in a hybrid zone between sex‐role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male–male competition in maintaining or eroding species barriers. Although female–female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex‐role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female–female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex‐role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.     

Discordant patterns of geographic variation between mitochondrial and microsatellite markers in the Mexican black iguana (Ctenosaura pectinata) in a contact zone

Aim To delimit the distribution of matrilines of the black iguana (Ctenosaura pectinata), detect potential contact zones between them, and test the hypothesis of geographic concordance in the structuring of matrilines with respect to genetic structuring across the nuclear genome. Location Pacific coast of Mexico in the states of Nayarit, Jalisco, Colima, Michoacán in the lowlands adjacent to the western end of the Mexican Volcanic Belt. Methods Tissue samples were obtained every 10–20 km along a transect across the range of three previously described mitochondrial DNA (mtDNA) clades of C. pectinata. Iguanas were genotyped with 11 nuclear microsatellites and maternal lineages were inferred based on mtDNA sequences. Geographic structure and geographic concordance between biparentally and maternally inherited markers were analysed with clustering methods, F_ST and N_ST indexes as well as haplotype networks. Results Geographic structure was evident for both markers and the existence of contact zones confirmed. Nevertheless, the distribution of nuclear and mtDNA genetic variation is not geographically concordant. Four matrilines exist in the area, whereas only two nuclear clusters occur. A contact zone between nuclear clusters extends along a distance of c. 60 km, and introgression is evident in individuals from several localities. Main conclusions High levels of genetic diversity were detected on the western coast of Mexico. Historical and contemporaneous processes seem to shape the distribution of genetic variation in C. pectinata. There are no evident geographic, environmental or ecological barriers that coincide with genetic subdivisions, and the observed mtDNA structure is likely to be the result of past climatic changes. The less structured distribution of nuclear genetic variation is consistent with the homogenizing effect of male‐biased dispersal. Our findings have implications for the taxonomy and conservation strategies for this threatened species, and highlight the fact that geographic analyses of both cytoplasmic and nuclear genetic variation are important for the meaningful inference of evolutionary and demographic histories.     

Genetic,morphological, and niche variation in the widely hybridizing Rhus integrifolia-Rhus ovata species complex

Hybridization and introgression are common processes among numerous plant species that present both challenges and opportunities for studies of species delimitation, phylogenetics, taxonomy and adaptation. Rhus integrifolia and R. ovata are two ecologically important shrubs native to the southwestern USA and Mexico, and are known to hybridize frequently, but the morphological, genetic and ecological implications of hybridization in these species are poorly studied on a broad geographic scale. Analyses were conducted using leaf morphology, genetic variation of plastid and nuclear loci, and species distribution models for both species and their putative hybrid introgressants across 19 localities in California and Arizona, USA. These analyses revealed evidence for morphological and genetic distinction among localities comprising putative parental species, but a high degree of morpho-genetic intermediacy among localities with putative hybrids. Comparison of morphological and genetic population structure among localities revealed evidence for putative local adaptation or widespread phenotypic plasticity. Multiple regression models identified a weak but statistically significant negative association between leaf area and precipitation. Finally, species distribution modeling inferred northward range shifts over time, with both species predicted to occupy more coastal regions in the future, possibly increasing the frequency of hybridization among them. These findings underscore the importance of integrative assessment of multiple data sources in the study of hybridizing species and highlight the R. integrifolia-ovata complex as a powerful model for investigating the adaptive implications of hybridization.     

Genetic, morphological, and ecological characterization of a hybrid zone that spans a migratory divide

This study characterizes a hybrid zone that spans a migratory divide between subspecies of the Swainson's thrush (Catharus ustulatus), a long distance migratory songbird, in the Coast Mountains of British Columbia. To assess the potential for a barrier to gene flow between the subspecies, I: (1) analyzed the shape and width of genetic and morphological clines relative to estimates of dispersal distance, (2) assessed the ratio of parental to hybrid genotypes across the hybrid zone, (3) estimated population density across the hybrid zone, and (4) compared the spatial relationship between the hybrid zone and an existing environmental gradient. The results indicate that the hybrid zone is characterized by mostly concordant character clines that are narrow relative to dispersal, the absence of a hybrid swarm, and low population density at the center of the zone. This hybrid zone and additional regions of contact between these subspecies are found on the border between coastal and interior climatic regions throughout the Pacific Northwest. An identified shift in the location, but not the width, of the mtDNA cline relative to the nuclear clines is consistent with asymmetrical hybridization. Neutral diffusion of populations following secondary contact and hybrid superiority within an ecotone are insufficient explanations for the observed patterns. The hypothesis that best fits the data is that the Swainson's thrush hybrid zone is a tension zone maintained by dispersal and ecologically mediated barriers to gene flow.     

Genetic evidence for natural hybridization and apparent introgression between freshwater snail species (Viviparus ater and V. contectus)

We examined six allozyme markers and shell morphology to study gene flow between the naturally sympatric freshwater snail species Viviparus ater and V. contectus in Lake Garda, Italy, using offspring from experimental crosses and snails collected from natural populations. Hybrid offspring obtained from experimental crosses and hybrids collected from natural populations were heterozygous at four or five out of five diagnostic loci. Shell morphology was a poor predictor of the hybrid status of individual snails. Backcrossed offspring from experimentally bred hybrids and either of the parental species morphologically resembled the parental species. There was no evidence of segregation distortion at the six loci. Linkage analysis revealed one pair of linked loci (GPI and PNP). These cross experiments indicated a Mendelian type hybrid system in which gene introgression may occur. F₁ hybrids were found at four out of six sampling sites of Lake Garda, Italy. Local frequencies of F₁ hybrids ranged from 0% to 1.6% (estimated average = 0.74%) of the total population including both species. Alleles typical of V. ater were found at low frequencies in V. contectus at all six sampling sites. Alleles typical of V. contectus were found at low frequencies in V. ater at three out of six sampling sites. This is consistent with the hypothesis of introgression in both directions.     

Effect of nesting environment on incubation temperature and hatching success of Morelet's crocodile (Crocodylus moreletii) in an urban lake of Southeastern Mexico

Incubation temperature is an important aspect in terms of biological performance among crocodiles, and several controlled experiments have demonstrated a significant relationship between incubation temperature, success in hatching and survival of hatchlings. However, a few studies have tested these relationships in the wild. The objective of this study was to determine the relationship of nest characteristics and environment (hatch year, nest basal area and height, clutch size, distance to shore line, and vegetation cover), to incubation temperature and hatching success among Morelet's crocodile (Crocodylus moreletii). The study was carried out during the nesting seasons of Morelet's crocodile, from 2007 to 2009 in the Laguna de Las Ilusiones, an urban lake located in Villahermosa, Tabasco, Mexico. We physically characterized 18 nests and inserted a temperature data logger in each nest chamber. At the end of the nesting season and prior to hatching, we recovered the crocodile eggs and data loggers and calculated hatching success, under laboratory conditions. We related the environmental variables of the nest with the mean and fluctuation (standard deviation) of nest temperature, using linear models. We also related the environmental variables affecting the nest, to mean nest temperature and fluctuation in incubation temperature and to hatching success, using linear models. Although we found differences in incubation temperature between nests, mean incubation temperature did not differ between years, but there were differences in nest thermal fluctuation between years. The mean incubation temperature for 11 nests (61.1%) was lower than the suggested Female–Male pivotal temperature (producing 50% of each sex) for this species, and all hatchlings obtained were males. There were no differences in clutch size between years, but hatching success varied. Our study indicates that hatching success depends on certain environmental variables and nest conditions to which the eggs are subjected, including season, nest size and clutch size. We also discuss the importance of the fluctuation of incubation temperature on hatching success and sex determination.