THE GENETIC STRUCTURE OF A HYBRID ZONE BETWEEN TWO CHROMOSOME RACES OF THE SCELOPORUS GRAMMICUS COMPLEX (SAURIA, PHRYNOSOMATIDAE) IN CENTRAL MEXICO

Allometric analysis of sexual dimorphism and morphological variation in two chromosome races of the Sceloporus grammicus complex (Squamata: Phrynosomatidae) from Mexico

Sexual dimorphism is a widespread feature in the Animal Kingdom. In lizards of the Sceloporus grammicus complex, studies of sexual dimorphism that analyze the allometric trajectories of body traits remain unexplored. Here we investigate sexual dimorphism in key phenotypic traits, including body size (snout-vent length, SVL) as well as head length (HL), head width (HW), and forearm length (FL). We use an allometric approach to detect differences in scale relationships among body parts in the S. grammicus complex in Mexico. We focus on two chromosomal races within this complex, F5 (2n = 34) and FM2 (2n = 46). In the complex, we found that males are larger than females in all morphological variables, and this pattern was confirmed in both races. We determined negative allometric trajectories (SVLvs.HL and HW), isometry (SVLvs.FL) and intersexual differences in the slopes of the SVLvs.HL and HW; the males showed steeper slopes. Thus, the growth of the head is more pronounced in males than females. Additionally, we found between-race differences in these trajectories (SVLvs.FL) and in all morphological variables (F5 lizards are larger than those of the FM2 race), which correlate with their chromosomal divergence. We discuss biological implications of our findings in relation to sexual selection and natural selection.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Female preference for sympatric vs. allopatric male throat color morphs in the mesquite lizard (Sceloporus grammicus) species complex

Color polymorphic sexual signals are often associated with alternative reproductive behaviors within populations, and the number, frequency, or type of morphs present often vary among populations. When these differences lead to assortative mating by population, the study of such polymorphic taxa may shed light on speciation mechanisms. We studied two populations of a lizard with polymorphic throat color, an important sexual signal. Males in one population exhibit orange, yellow, or blue throats; whereas males in the other exhibit orange, yellow, or white throats. We assessed female behavior when choosing between allopatric and sympatric males. We asked whether females discriminated more when the allopatric male was of an unfamiliar morph than when the allopatric male was similar in coloration to the sympatric male. We found that female rejection of allopatric males relative to sympatric males was more pronounced when males in a pair were more different in throat color. Our findings may help illuminate how behavioral responses to color morph differences between populations with polymorphic sexual signals contribute to reproductive isolation.

Asymmetric reproductive isolation between terminal forms of the salamander ring species Ensatina eschscholtzii revealed by fine-scale genetic analysis of a hybrid zone

Background

Ring species, exemplified by salamanders of the Ensatina eschscholtzii complex, represent a special window into the speciation process because they allow the history of species formation to be traced back in time through the geographically differentiated forms connecting the two terminal forms of the ring. Of particular interest is the nature and extent of reproductive isolation between the geographically terminal forms, in this case E. e. eschscholtzii and E. e. klauberi. Previous studies have documented infrequent hybridization at the end of the ring. Here, we report the first fine-scale genetic analysis of a hybrid zone between the terminal forms in southern California using individual-based Bayesian analyses of multilocus genetic data to estimate levels and direction of hybridization and maximum-likelihood analysis of linkage disequilibrium and cline shape to make inferences about migration and selection in the hybrid zone.

Results

The center of the hybrid zone has a high proportion of hybrids, about half of which were classified as F1s. Clines are narrow with respect to dispersal, and there are significant deviations from Hardy-Weinberg equilibrium as well as nonrandom associations (linkage disequilibria) between alleles characteristic of each parental type. There is cytonuclear discordance, both in terms of introgression and the geographic position of mitochondrial versus nuclear clines. Genetic disequilibrium is concentrated on the eschscholtzii side of the zone. Nearly all hybrids possess klauberi mtDNA, indicating that most hybrids are formed from female klauberi mating with male eschscholtzii or male hybrids (but not vice versa).

Conclusions

Our results are consistent with a tension zone trapped at an ecotone, with gene combinations characteristic of klauberi showing up on the eschscholtzii side of the zone due to asymmetric hybridization. We suggest that the observed asymmetry is best explained by increased discriminatory power of eschscholtzii females, or asymmetric postzygotic isolation. The relatively high frequency of hybrids, particularly F1s, contrasts with other contacts between the terminal forms, and with other contacts between other divergent Ensatina lineages, highlighting the diverse outcomes of secondary contact within a single species complex.

Two behavioural traits promote fine-scale species segregation and moderate hybridisation in a recovering sympatric fur seal population

BACKGROUND

In systems where two or more species experience secondary contact, behavioural factors that regulate interspecific gene flow may be important for maintaining species boundaries and reducing the incidence of hybridisation. At subantarctic Macquarie Island, two species of fur seal breed in close proximity to one another, hybridise at very high levels (up to 21% of hybrid pups are born annually), yet retain discrete gene pools. Using spatial and genetic information collected for pups and adults over twelve years, we assessed two behavioural traits - inter-annual site fidelity and differences in habitat use between the species - as possible contributors to the maintenance of this species segregation. Further, we explored the breakdown of these traits in pure-species individuals and hybrids.

RESULTS

We found virtually complete spatial segregation of the parental species, with only one exception; a single territory that contained adults of both species and also the highest concentration of hybrid pups. The spatial distribution of each species was closely linked to habitat type (pebbled vs boulder beaches), with members of each species breeding almost exclusively on one type or the other but hybrids breeding on both or at the junction between habitats. Inter-annual site fidelity was high for both sexes of pure-species adults, with 66% of females and all males returning to the same territory or a neighbouring one in different years. An important consequence for pure females of breeding on the 'wrong' habitat type, and thus in a heterospecific aggregation, was the production of hybrid pups. Low habitat fidelity of hybrid females facilitated bi-directional backcrossing, resulting in more diverse hybrid offspring.

CONCLUSION

In a disturbed system where two sympatric fur seal species breed in close proximity, discrete gene pools are retained by extremely fine-scale and strong spatial segregation of the species. Two behavioural traits were found to be important in maintaining this stable population structure, and habitat type was a strong indicator of where species locate and a potentially powerful predictor of future directions of hybridisation. A direct consequence of the breakdown of this trait was the production of hybrid offspring, which may have severe implications if hybrids have reduced fitness.

Cryptic diversity in vertebrates: molecular data double estimates of species diversity in a radiation of Australian lizards (Diplodactylus, Gekkota)

A major problem for biodiversity conservation and management is that a significant portion of species diversity remains undocumented (the 'taxonomic impediment'). This problem is widely acknowledged to be dire among invertebrates and in developing countries; here, we demonstrate that it can be acute even in conspicuous animals (reptiles) and in developed nations (Australia). A survey of mtDNA, allozyme and chromosomal variation in the Australian gecko, genus Diplodactylus, increases overall species diversity estimates from 13 to 29. Four nominal species each actually represent multi-species complexes; three of these species complexes are not even monophyletic. The high proportion of cryptic species discovered emphasizes the importance of continuing detailed assessments of species diversity, even in apparently well-known taxa from industrialized countries.

The roles of allopatric divergence and natural selection in quantitative trait variation across a secondary contact zone in the lizard Anolis roquet

Populations of the Caribbean lizard, Anolis roquet, are thought to have experienced long periods of allopatry before recent secondary contact. To elucidate the effects of past allopatry on population divergence in A. roquet, we surveyed parallel transects across a secondary contact zone in northeastern Martinique. We used diagnostic molecular mitochondrial DNA markers to test fine-scale association of mitochondrial DNA lineage and geological region, multivariate statistical techniques to explore quantitative trait pattern, and cline fitting techniques to model trait variation across the zone of secondary contact. We found that lineages were strongly associated with geological regions along both transects, but quantitative trait patterns were remarkably different. Patterns of morphological and mitochondrial DNA variation were consistent with a strong barrier to gene flow on the coast, whereas there were no indications of barriers to gene flow in the transitional forest. Hence, the coastal populations behaved as would be predicted by an allopatric model of divergence in this complex, while those in the transitional forest did not, despite the close proximity of the transects and their shared geological history. Patterns of geographical variation in this species complex, together with environmental data, suggest that on balance, selection regimes on either side of the secondary contact zone in the transitional forest may be more convergent, while those either side of the secondary contact zone on the coast are more divergent. Hence, the evolutionary consequences of allopatry may be strongly influenced by local natural selection regimes.

Extensive Natural Hybridization Between Two Geckos, Gekko tawaensis and Gekko japonicus (Reptilia: Squamata), Throughout Their Broad Sympatric Area

The status of natural hybridization between the two Japanese geckos, Gekko tawaensis and Gekko japonicus, was surveyed by use of 15 diagnostic allozyme loci. Of 438 specimens examined, 9 were identified as F(l), 1 was a first backcross with G. japonicus, and 15 were identified as more successive generations. Hybridizations were detected at 7 of the 15 localities interspersed throughout a broad sympatric area of the two species, forming a mosaic hybrid zone. A comparison of species-hybrid composition between 2-year samples from a single locality and a 5-year interval showed no evidence for progressive introgression or establishment of a hybrid swarm, despite constant emergences of F(l) hybrids. Nonprevalence of the hybrid genotypes was also indicated by the analysis using hybrid index scores for all other localities examined. These results suggest that strong selection acts against hybrid genotypes. Morphological features of hybrid individuals were also provided.

Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies

We describe extensions to the method of Pritchard et al. for inferring population structure from multilocus genotype data. Most importantly, we develop methods that allow for linkage between loci. The new model accounts for the correlations between linked loci that arise in admixed populations ("admixture linkage disequilibium"). This modification has several advantages, allowing (1) detection of admixture events farther back into the past, (2) inference of the population of origin of chromosomal regions, and (3) more accurate estimates of statistical uncertainty when linked loci are used. It is also of potential use for admixture mapping. In addition, we describe a new prior model for the allele frequencies within each population, which allows identification of subtle population subdivisions that were not detectable using the existing method. We present results applying the new methods to study admixture in African-Americans, recombination in Helicobacter pylori, and drift in populations of Drosophila melanogaster. The methods are implemented in a program, structure, version 2.0, which is available at http://pritch.bsd.uchicago.edu.

Mating signal variation and bimodality in a mosaic hybrid zone between Chorthippus grasshopper species

The grasshoppers Chorthippus brunneus and Chorthippus jacobsi are easily distinguished by male calling song and the number of stridulatory pegs on the hind femur, and form a mosaic hybrid zone in northern Spain. In this paper, we fit a two-dimensional cline to variation in male calling song characters, which are of particular interest as they are likely to be involved in mate choice by females. As with variation in peg number, local habitat makes only a small contribution in explaining deviations in mean song score from clinal expectations. However, the fitted width of the cline for song characters is significantly narrower than for peg number, suggesting that mating signals may be associated with reduced hybrid fitness in the field and that recombination rates are sufficient to allow clines for different characters to diverge in width. Despite this, estimates for the overall elevation in linkage disequilibrium at the zone center, based on covariance between peg and song characters, reveal a substantial overrepresentation of parental genotypes at the cline center relative to the expectations of a tension zone of similar width. Examination of covariance at individual sites reveals that this inflated estimate of linkage disequilibrium is caused by several sites where the distribution of phenotypes is effectively bimodal. This substantial variation in linkage disequilibrium at the cline center could result from local variation in the strength of assortative mating or selection against hybrids, or may reflect the long-distance colonization of empty habitat from outside the hybrid zone, which would continually create new contacts between parental genotypes at the cline center. Hybrid zones like this, in which strong linkage disequilibrium occurs in some situations but not in others, are of particular relevance to speciation research and allow investigation of the spread of combinations of alleles through different genetic and ecological backgrounds.

Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

A comparison of nuclear and mitochondrial cline shapes in a hybrid zone in the Sceloporus grammicus complex (Squamata; Phrynosomatidae)

The F5 and FM2 chromosome races of the Sceloporus grammicus complex form a hybrid zone in the Mexican state of Hidalgo. Previous studies of this zone have assessed genetic structure by averaging estimates of shape and width across three diagnostic chromosome markers. This approach is likely to mask subtle differences in cline shape among loci (e.g. selected vs. neutral), and obscure any displacement of cline centres (if present). Here we use maximum likelihood methods to construct the best fitting individual clines for three chromosomal markers, and also add two new markers; the mitochondrial DNA (mtDNA) locus, and the nuclear ribosomal DNA (rDNA) repeat. For each locus, hybrid zone models were fitted by cline shape and width, and the position and number of segments describing the centre of the zone. Pairwise comparisons between all clines revealed concordance between chromosomes 2 and 6, but significant discordance in cline structure among all other paired combinations. The concordance of chromosomes 2 and 6 suggests that these clines are maintained by genome-wide forces. The discordance of the chromosome 1 cline suggests an influence of asymmetric introgression, while the mtDNA cline is probably influenced by selection and drift. The rDNA locus reveals a pattern best explained by either extreme asymmetric introgression or gene conversion. The structure of zone indicates that genome-wide processes and locus specific selective forces as well as drift, are operating to different degrees on different loci. The locus-by-locus approach used here permits a finer discrimination among possible mechanisms responsible for the maintenance of the individual clines.

A comparison of multilocus clines maintained by environmental adaptation or by selection against hybrids

There has recently been considerable debate over the relative importance of selection against hybrids ("endogenous" selection) vs. adaptation to different environments ("exogenous") in maintaining stable hybrid zones and hence in speciation. Single-locus models of endogenous and exogenous viability selection generate clines of similar shape, but the comparison has not been extended to multilocus systems, which are both quantitatively and qualitatively very different from the single-locus case. Here we develop an analytical multilocus model of differential adaptation across an environmental transition and compare it to previous heterozygote disadvantage models. We show that the shape of clines generated by exogenous selection is indistinguishable from that generated by endogenous selection. A stochastic simulation model is used to test the robustness of the analytical description to the effects of drift and strong selection, and confirms the prediction that pairwise linkage disequilibria are predominantly generated by migration. However, although analytical predictions for the width of clines maintained by heterozygote disadvantage fit well with the simulation results, those for environmental adaptation are consistently too narrow; reasons for the discrepancy are discussed. There is a smooth transition between a system in which a set of loci effectively act independently of each other and one in which they act as a single nonrecombining unit.

Tracking the evolution of insecticide resistance in the mosquito Culex pipiens

The evolution of pesticide resistance provides some of the most striking examples of darwinian evolution occurring over a human life span. Identification of resistance alleles opens an outstanding framework in which to study the evolution of adaptive mutations from the beginning of pesticide application, the evolution of interactions between alleles (dominance) or between loci (epistasis). Here we show that resistance alleles can also be used as markers to dissect population processes at a microevolutionary scale. We have focused on the antagonistic roles of selection and migration involved in the dynamics of local adaptation with reference to allelic frequencies at two resistance loci in the mosquito Culex pipiens. We find that their frequencies follow an annual cycle of large amplitude (25%), and we precisely unravel the seasonal variation of migration and selection underlying this cycle. Our results provide a firm basis on which to devise an insecticide treatment strategy that will better control the evolution of resistance genes and the growth of mosquito populations.

Chromosomal heterozygosity and fertility in house mice (Mus musculus domesticus) from Northern Italy

Following the discovery of over 40 Robertsonian (Rb) races of Mus musculus domesticus in Europe and North Africa, the house mouse has been studied extensively as an ideal model to determine the chromosomal changes that may cause or accompany speciation. Current models of chromosomal speciation are based on the assumption that heterozygous individuals have a particularly low fertility, although recent studies indicate otherwise. Despite their importance, fertility estimates for the house mouse are incomplete because traditional measurements, such as anaphase I nondisjunction and germ cell death, are rarely estimated in conjunction with litter size. In an attempt to bridge this gap, we have taken advantage of the house mouse hybrid zone in Upper Valtellina (Lombardy, Italy) in which five Rb races interbreed. We present data on the fertility of naturally occurring ("wild-caught") hybrids and of offspring from laboratory crosses of wild-caught mice ("laboratory-reared"), using various measurements. Wild-caught mice heterozygous for one fusion were more infertile than predicted from past studies, possibly due to genic hybridity; laboratory-reared heterozygotes carrying seven or eight trivalents at meiosis I and heterozygotes carrying one pentavalent also had low fertilities. These low fertilities are especially significant given the probable occurrence of a reinforcement event in Upper Valtellina.

Evaluating gene flow using selected markers: a case study

The extent to which an organism is locally adapted in an environmental pocket depends on the selection intensities inside and outside the pocket, on migration, and on the size of the pocket. When two or more loci are involved in this local adaptation, measuring their frequency gradients and their linkage disequilbria allows one to disentangle the forces-migration and selection-acting on the system. We apply this method to the case of a local adaptation to organophosphate insecticides in the mosquito Culex pipiens pipiens in southern France. The study of two different resistance loci allowed us to estimate with support limits gene flow as well as selection pressure on insecticide resistance and the fitness costs associated with each locus. These estimates permit us to pinpoint the conditions for the maintenance of this pocket of adaptation as well as the effect of the interaction between the two resistance loci.