Rapid spread of invasive genes into a threatened native species

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species-specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri-specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure), bidirectional estimates of gene flow (migrate-n and IMa2) and phylogeny-based methods (*beast, bucky) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species-specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri-like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.

Fitness dynamics within a poplar hybrid zone: II. Impact of exotic sex on native poplars in an urban jungle

Trees bearing novel or exotic gene components are poised to contribute to the bioeconomy for a variety of purposes such as bioenergy production, phytoremediation, and carbon sequestration within the forestry sector, but sustainable release of trees with novel traits in large-scale plantations requires the quantification of risks posed to native tree populations. Over the last century, exotic hybrid poplars produced through artificial crosses were planted throughout eastern Canada as ornamentals or windbreaks and these exotics provide a proxy by which to examine the fitness of exotic poplar traits within the natural environment to assess risk of exotic gene escape, establishment, and spread into native gene pools. We assessed postzygotic fitness traits of native and exotic poplars within a naturally regenerated stand in eastern Canada (Quebec City, QC). Pure natives (P. balsamifera and P. deltoides spp. deltoides), native hybrids (P. deltoides × P. balsamifera), and exotic hybrids (trees bearing Populus nigra and P. maximowiczii genetic components) were screened for reproductive biomass, yield, seed germination, and fungal disease susceptibility. Exotic hybrids expressed fitness traits intermediate to pure species and were not significantly different from native hybrids. They formed fully viable seed and backcrossed predominantly with P. balsamifera. These data show that exotic hybrids were not unfit and were capable of establishing and competing within the native stand. Future research will seek to examine the impact of exotic gene regions on associated biotic communities to fully quantify the risk exotic poplars pose to native poplar forests.

Hybridization in a warmer world

Climate change is profoundly affecting the evolutionary trajectory of individual species and ecological communities, in part through the creation of novel species assemblages. How climate change will influence competitive interactions has been an active area of research. Far less attention, however, has been given to altered reproductive interactions. Yet, reproductive interactions between formerly isolated species are inevitable as populations shift geographically and temporally as a result of climate change, potentially resulting in introgression, speciation, or even extinction. The susceptibility of hybridization rates to anthropogenic disturbance was first recognized in the 1930s. To date, work on anthropogenically mediated hybridization has focused primarily on either physical habitat disturbance or species invasion. Here, I review recent literature on hybridization to identify how ecological responses to climate change will increase the likelihood of hybridization via the dissolution of species barriers maintained by habitat, time, or behavior. Using this literature, I identify several cases where novel hybrid zones have recently formed, likely as a result of changing climate. Future research should focus on identifying areas and taxonomic groups where reproductive species interactions are most likely to be influenced by climate change. Furthermore, a better understanding of the evolutionary consequences of climate-mediated secondary contact is urgently needed. Paradoxically, hybridization is both a major conservation concern and an important source of novel genetic and phenotypic variation. Hybridization may therefore both contribute to increasing rates of extinction and stimulate the creation of novel phenotypes that will speed adaptation to novel climates. Predicting which result will occur following secondary contact will be an important contribution to conservation for many species.

Genomic replacement of native Cobitis lutheri with introduced C. tetralineata through a hybrid swarm following the artificial connection of river systems

River connections via artificial canals will bring about secondary contacts between previously isolated fish species. Here, we present a genetic consequence of such a secondary contact between Cobitis fish species, C. lutheri in the Dongjin River, and C. tetralineata in the Seomjin River in Korea. The construction of water canals about 80 years ago has unidirectionally introduced C. tetralineata into the native habitat of C. lutheri, and then these species have hybridized in the main stream section of the Dongjin River. According to the divergence population genetic analyses of DNA sequence data, the two species diverged about 3.3 million years ago, which is interestingly coincident with the unprecedented paleoceanographic change that caused isolations of the paleo-river systems in northeast Asia due to sea-level changes around the late Pliocene. Multilocus genotypic data of nine microsatellites and three nuclear loci revealed an extensively admixed structure in the hybrid zone with a high proportion of various post-F1 hybrids. Surprisingly, pure native C. lutheri was absent in the hybrid zone in contrast to the 7% of pure C. tetralineata. Such a biased proportion must have resulted from the dominant influence of continually introducing C. tetralineata on the native C. lutheri which has no supply of natives from other tributaries to the hybrid zone due to numerous low-head dams. In addition, mating experiments indicated that there is no discernible reproductive isolation between them. All the results suggest that the gene pool of native C. lutheri is being rapidly replaced by that of continually introducing C. tetralineata through a hybrid swarm for the last 80 years, which will ultimately lead to the genomic extinction of natives in this hybrid zone.

Intensive removal of signal crayfish (Pacifastacus leniusculus) from rivers increases numbers and taxon richness of macroinvertebrate species

Invasive species are a major cause of species extinction in freshwater ecosystems, and crayfish species are particularly pervasive. The invasive American signal crayfish Pacifastacus leniusculus has impacts over a range of trophic levels, but particularly on benthic aquatic macroinvertebrates. Our study examined the effect on the macroinvertebrate community of removal trapping of signal crayfish from UK rivers. Crayfish were intensively trapped and removed from two tributaries of the River Thames to test the hypothesis that lowering signal crayfish densities would result in increases in macroinvertebrate numbers and taxon richness. We removed 6181 crayfish over four sessions, resulting in crayfish densities that decreased toward the center of the removal sections. Conversely in control sections (where crayfish were trapped and returned), crayfish density increased toward the center of the section. Macroinvertebrate numbers and taxon richness were inversely correlated with crayfish densities. Multivariate analysis of the abundance of each taxon yielded similar results and indicated that crayfish removals had positive impacts on macroinvertebrate numbers and taxon richness but did not alter the composition of the wider macroinvertebrate community. Synthesis and applications: Our results demonstrate that non-eradication-oriented crayfish removal programmes may lead to increases in the total number of macroinvertebrates living in the benthos. This represents the first evidence that removing signal crayfish from riparian systems, at intensities feasible during control attempts or commercial crayfishing, may be beneficial for a range of sympatric aquatic macroinvertebrates.

Changes in the genetic structure of Atlantic salmon populations over four decades reveal substantial impacts of stocking and potential resiliency

While the stocking of captive-bred fish has been occurring for decades and has had substantial immediate genetic and evolutionary impacts on wild populations, its long-term consequences have only been weakly investigated. Here, we conducted a spatiotemporal analysis of 1428 Atlantic salmon sampled from 1965 to 2006 in 25 populations throughout France to investigate the influence of stocking on the neutral genetic structure in wild Atlantic salmon (Salmo salar) populations. On the basis of the analysis of 11 microsatellite loci, we found that the overall genetic structure among populations dramatically decreased over the period studied. Admixture rates among populations were highly variable, ranging from a nearly undetectable contribution from donor stocks to total replacement of the native gene pool, suggesting extremely variable impacts of stocking. Depending on population, admixture rates either increased, remained stable, or decreased in samples collected between 1998 and 2006 compared to samples from 1965 to 1987, suggesting either rising, long-lasting or short-term impacts of stocking. We discuss the potential mechanisms contributing to this variability, including the reduced fitness of stocked fish and persistence of wild locally adapted individuals.

Alternative forms for genomic clines

Understanding factors regulating hybrid fitness and gene exchange is a major research challenge for evolutionary biology. Genomic cline analysis has been used to evaluate alternative patterns of introgression, but only two models have been used widely and the approach has generally lacked a hypothesis testing framework for distinguishing effects of selection and drift. I propose two alternative cline models, implement multivariate outlier detection to identify markers associated with hybrid fitness, and simulate hybrid zone dynamics to evaluate the signatures of different modes of selection. Analysis of simulated data shows that previous approaches are prone to false positives (multinomial regression) or relatively insensitive to outlier loci affected by selection (Barton's concordance). The new, theory-based logit-logistic cline model is generally best at detecting loci affecting hybrid fitness. Although some generalizations can be made about different modes of selection, there is no one-to-one correspondence between pattern and process. These new methods will enhance our ability to extract important information about the genetics of reproductive isolation and hybrid fitness. However, much remains to be done to relate statistical patterns to particular evolutionary processes. The methods described here are implemented in a freely available package “HIest” for the R statistical software (CRAN; http://cran.r-project.org/).

Asymmetric introgression in the horticultural living fossil cycas sect. Asiorientales using a genome-wide scanning approach

The Asian cycads are mostly allopatric, distributed in small population sizes. Hybridization between allopatric species provides clues in determining the mechanism of species divergence. Horticultural introduction provides the chance of interspecific gene flow between allopatric species. Two allopatrically eastern Asian Cycas sect. Asiorientales species, C. revoluta and C. taitungensis, which are widely distributed in Ryukyus and Fujian Province and endemic to Taiwan, respectively, were planted in eastern Taiwan for horticultural reason. Higher degrees of genetic admixture in cultivated samples than wild populations in both cycad species were detected based on multilocus scans by neutral AFLP markers. Furthermore, bidirectional but asymmetric introgression by horticultural introduction of C. revoluta is evidenced by the reanalyses of species associated loci, which are assumed to be diverged after species divergence. Partial loci introgressed from native cycad to the invaders were also detected at the loci of strong species association. Consistent results tested by all neutral loci, and the species-associated loci, specify the recent introgression from the paradox of sharing of ancestral polymorphisms. Phenomenon of introgression of cultivated cycads implies niche conservation among two geographic-isolated cycads, even though the habitats of the extant wild populations of two species are distinct.

Genomic patterns of introgression in rainbow and westslope cutthroat trout illuminated by overlapping paired-end RAD sequencing

Rapid and inexpensive methods for genomewide single nucleotide polymorphism (SNP) discovery and genotyping are urgently needed for population management and conservation. In hybridized populations, genomic techniques that can identify and genotype thousands of species-diagnostic markers would allow precise estimates of population- and individual-level admixture as well as identification of 'super invasive' alleles, which show elevated rates of introgression above the genomewide background (likely due to natural selection). Techniques like restriction-site-associated DNA (RAD) sequencing can discover and genotype large numbers of SNPs, but they have been limited by the length of continuous sequence data they produce with Illumina short-read sequencing. We present a novel approach, overlapping paired-end RAD sequencing, to generate RAD contigs of >300-400 bp. These contigs provide sufficient flanking sequence for design of high-throughput SNP genotyping arrays and strict filtering to identify duplicate paralogous loci. We applied this approach in five populations of native westslope cutthroat trout that previously showed varying (low) levels of admixture from introduced rainbow trout (RBT). We produced 77 141 RAD contigs and used these data to filter and genotype 3180 previously identified species-diagnostic SNP loci. Our population-level and individual-level estimates of admixture were generally consistent with previous microsatellite-based estimates from the same individuals. However, we observed slightly lower admixture estimates from genomewide markers, which might result from natural selection against certain genome regions, different genomic locations for microsatellites vs. RAD-derived SNPs and/or sampling error from the small number of microsatellite loci (n = 7). We also identified candidate adaptive super invasive alleles from RBT that had excessively high admixture proportions in hybridized cutthroat trout populations.

Lethal effects of water quality on threatened California salamanders but not on co-occurring hybrid salamanders

Biological invasions and habitat alteration are often detrimental to native species, but their interactions are difficult to predict. Interbreeding between native and introduced species generates novel genotypes and phenotypes, and human land use alters habitat structure and chemistry. Both invasions and habitat alteration create new biological challenges and opportunities. In the intensively farmed Salinas Valley, California (U.S.A.), threatened California tiger salamanders (Ambystoma californiense) have been replaced by hybrids between California tiger salamander and introduced barred tiger salamanders (Ambystoma tigrinum mavortium). We conducted an enclosure experiment to examine the effects habitat modification and relative frequency of hybrid and native California tiger salamanders have on recruitment of salamanders and their prey, Pacific chorus frogs (Pseudacris regilla). We tested whether recruitment differed among genetic classes of tiger salamanders (hybrid or native) and pond hydroperiod (seasonal or perennial). Roughly 6 weeks into the experiment, 70% (of 378 total) of salamander larvae died in 4 out of 6 ponds. Native salamanders survived (n = 12) in these ponds only if they had metamorphosed prior to the die-offs. During die-offs, all larvae of native salamanders died, whereas 56% of hybrid larvae died. We necropsied native and hybrid salamanders, tested water quality, and queried the California Department of Pesticide Regulation database to investigate possible causes of the die-offs. Salamander die-offs, changes in the abundance of other community members (invertebrates, algae, and cyanobacteria), shifts in salamander sex ratio, and patterns of pesticide application in adjacent fields suggest that pesticide use may have contributed to die-offs. That all survivors were hybrids suggests that environmental stress may promote rapid displacement of native genotypes.

Transgressive Hybrids as Hopeful Monsters

The origin of novelty is a critical subject for evolutionary biologists. Early geneticists speculated about the sudden appearance of new species via special macromutations, epitomized by Goldschmidt’s infamous “hopeful monster”. Although these ideas were easily dismissed by the insights of the Modern Synthesis, a lingering fascination with the possibility of sudden, dramatic change has persisted. Recent work on hybridization and gene exchange suggests an underappreciated mechanism for the sudden appearance of evolutionary novelty that is entirely consistent with the principles of modern population genetics. Genetic recombination in hybrids can produce transgressive phenotypes, “monstrous” phenotypes beyond the range of parental populations. Transgressive phenotypes can be products of epistatic interactions or additive effects of multiple recombined loci. We compare several epistatic and additive models of transgressive segregation in hybrids and find that they are special cases of a general, classic quantitative genetic model. The Dobzhansky-Muller model predicts “hopeless” monsters, sterile and inviable transgressive phenotypes. The Bateson model predicts “hopeful” monsters with fitness greater than either parental population. The complementation model predicts both. Transgressive segregation after hybridization can rapidly produce novel phenotypes by recombining multiple loci simultaneously. Admixed populations will also produce many similar recombinant phenotypes at the same time, increasing the probability that recombinant “hopeful monsters” will establish true-breeding evolutionary lineages. Recombination is not the only (or even most common) process generating evolutionary novelty, but might be the most credible mechanism for sudden appearance of new forms.

Parallel tagged amplicon sequencing reveals major lineages and phylogenetic structure in the North American tiger salamander (Ambystoma tigrinum) species complex

Modern analytical methods for population genetics and phylogenetics are expected to provide more accurate results when data from multiple genome-wide loci are analysed. We present the results of an initial application of parallel tagged sequencing (PTS) on a next-generation platform to sequence thousands of barcoded PCR amplicons generated from 95 nuclear loci and 93 individuals sampled across the range of the tiger salamander (Ambystoma tigrinum) species complex. To manage the bioinformatic processing of this large data set (344 330 reads), we developed a pipeline that sorts PTS data by barcode and locus, identifies high-quality variable nucleotides and yields phased haplotype sequences for each individual at each locus. Our sequencing and bioinformatic strategy resulted in a genome-wide data set with relatively low levels of missing data and a wide range of nucleotide variation. structure analyses of these data in a genotypic format resulted in strongly supported assignments for the majority of individuals into nine geographically defined genetic clusters. Species tree analyses of the most variable loci using a multi-species coalescent model resulted in strong support for most branches in the species tree; however, analyses including more than 50 loci produced parameter sampling trends that indicated a lack of convergence on the posterior distribution. Overall, these results demonstrate the potential for amplicon-based PTS to rapidly generate large-scale data for population genetic and phylogenetic-based research.

Estimating ancestry and heterozygosity of hybrids using molecular markers

Background

Hybridization, genetic mixture of distinct populations, gives rise to myriad recombinant genotypes. Characterizing the genomic composition of hybrids is critical for studies of hybrid zone dynamics, inheritance of traits, and consequences of hybridization for evolution and conservation. Hybrid genomes are often summarized either by an estimate of the proportion of alleles coming from each ancestral population or classification into discrete categories like F1, F2, backcross, or merely “hybrid” vs. “pure”. In most cases, it is not realistic to classify individuals into the restricted set of classes produced in the first two generations of admixture. However, the continuous ancestry index misses an important dimension of the genotype. Joint consideration of ancestry together with interclass heterozygosity (proportion of loci with alleles from both ancestral populations) captures all of the information in the discrete classification without the unrealistic assumption that only two generations of admixture have transpired.

Methods

I describe a maximum likelihood method for joint estimation of ancestry and interclass heterozygosity. I present two worked examples illustrating the value of the approach for describing variation among hybrid populations and evaluating the validity of the assumption underlying discrete classification.

Results

Naively classifying natural hybrids into the standard six line cross categories can be misleading, and false classification can be a serious problem for datasets with few molecular markers. My analysis underscores previous work showing that many (50 or more) ancestry informative markers are needed to avoid erroneous classification.

Conclusion

Although classification of hybrids might often be misleading, valuable inferences can be obtained by focusing directly on distributions of ancestry and heterozygosity. Estimating and visualizing the joint distribution of ancestry and interclass heterozygosity is an effective way to compare the genetic structure of hybrid populations and these estimates can be used in classic quantitative genetic methods for assessing additive, dominant, and epistatic genetic effects on hybrid phenotypes and fitness. The methods are implemented in a freely available package “HIest” for the R statistical software ( http://cran.r-project.org/web/packages/HIest/index.html).

Population genomics as a new tool for wildlife management

Admixture and introgression have varied effects on population viability and fitness. Admixture might be an important source of new alleles, particularly for small, geographically isolated populations. However, admixture might also cause outbreeding depression if populations are adapted to different ecological or climatic conditions. Because of the emerging use of translocation and admixture as a conservation and wildlife management strategy to reduce genetic load (termed genetic rescue), the possible effects of admixture have practical consequences (Bouzat et al. 2009; Hedrick & Fredrickson 2010). Importantly, genetic load and local adaptation are properties of individual loci and epistatic interactions among loci rather than properties of genomes. Likewise, the outcome and consequences of genetic rescue depend on the fitness effects of individual introduced alleles. In this issue of Molecular Ecology, Miller et al. (2012) use model-based, population genomic analyses to document locus-specific effects of a recent genetic rescue in the bighorn sheep population within the National Bison Range wildlife refuge (NBR; Montana, USA). They find a subset of introduced alleles associated with increased fitness in NBR bighorn sheep, some of which experienced accelerated introgression following their introduction. These loci mark regions of the genome that could constitute the genetic basis of the successful NBR bighorn sheep genetic rescue. Although population genomic analyses are frequently used to study local adaptation and selection (e.g. Hohenlohe et al. 2010; Lawniczak et al. 2010), this study constitutes a novel application of this analytical framework for wildlife management. Moreover, the detailed demographic data available for the NBR bighorn sheep population provide a rare and powerful source of information and allow more robust population genomic inference than is often possible.

Understanding admixture patterns in supplemented populations: a case study combining molecular analyses and temporally explicit simulations in Atlantic salmon

Genetic admixture between wild and introduced populations is a rising concern for the management of endangered species. Here, we use a dual approach based on molecular analyses of samples collected before and after hatchery fish introduction in combination with a simulation study to obtain insight into the mechanisms of admixture in wild populations. Using 17 microsatellites, we genotyped pre- and post-stocking samples from four Atlantic salmon populations supplemented with non-native fish to estimate genetic admixture. We also used individual-based temporally explicit simulations based on realistic demographic and stocking data to predict the extent of admixture. We found a low admixture by hatchery stocks within prestocking samples but moderate to high values in post-stocking samples (from 12% to 60%). The simulation scenarios best fitting the real data suggested a 10-25 times lower survival of stocked fish relative to wild individuals. Simulations also suggested relatively high dispersal rates of stocked and wild fish, which may explain some high levels of admixture in weakly stocked populations and the persistence of indigenous genotypes in heavily stocked populations. This study overall demonstrates that combining genetic analyses with simulations can significantly improve the understanding of admixture mechanisms in wild populations.

Genomic consequences of genetic rescue in an insular population of bighorn sheep (Ovis canadensis)

Genetic rescue is a management intervention whereby a small population is supplemented with individuals from other populations in an attempt to reverse the effects of inbreeding and increased genetic load. One such rescue was recently documented in the population of bighorn sheep (Ovis canadensis) within the National Bison Range wildlife refuge (Montana, USA). Here, we examine the locus-specific effects of rescue in this population using a newly developed genome-wide set of 195 microsatellite loci and first-generation linkage map. We found that the rate of introgression varied among loci and that 111 loci, 57% of those examined, deviated from patterns of neutral inheritance. The most common deviation was an excess of homozygous genotypes relative to neutral expectations, indicative of directional selection. As in previous study of this rescue, individuals with more introduced alleles had higher reproductive success and longevity. In addition, we found 30 loci, distributed throughout the genome, which seem to have individual effects on these life history traits. Although the potential for outbreeding depression is a major concern when translocating individuals between populations, we found no evidence of such effects in this population.

Did hybridization save the Norfolk Island boobook owl Ninox novaeseelandiae undulata?

The population of the Norfolk Island boobook owl Ninox novaeseelandiae undulata, a nocturnal bird restricted to the Australian territory of Norfolk Island, was reduced to a single female in 1986. Deliberate introduction of two males of its nearest relative, the New Zealand boobook N. n. novaeseelandiae, as a conservation intervention has allowed the taxon to persist on Norfolk Island, albeit in hybrid form. Although declared Extinct in 2000, a re-examination of this unique situation has concluded there is a strong argument that the taxon should be categorized as Critically Endangered because, on average, approximately half the nuclear genome of the original taxon and all the mitochondrial DNA is conserved in all living owls on the island. This thus represents a special case in which the taxon can be considered to be extant, in hybrid form, even though no pure-bred individuals survive. More generally, we suggest that, in exceptional cases, hybridization may not be a threat to highly threatened species and that guidelines are needed to determine when to consider hybrid populations as extant forms of the original taxon, and when to declare extinction through hybridization.

Species introduction promotes hybridization and introgression in Coregonus: is there sign of selection against hybrids?

Species introductions are considered one of the major drivers of biodiversity loss via ecological interactions and genetic admixture with local fauna. We examined two well-recognized fish species, native whitefish (Coregonus lavaretus) and introduced vendace (Coregonus albula), as well as their morphological hybrids in a single lake to test for selection against hybrids and backcrosses in the wild. A representative random subsample of 693 individuals (27.8%) was taken from the total catch of coregonids. This subsample was examined with the aim to select c. 50 individuals of pure whitefish (n = 52), pure vendace (n = 55) and putative hybrid (n = 19) for genetic analyses. The subsequent microsatellites and mitochondrial (mt) DNA analyses provided compelling evidence of hybridization and introgression. Of the 126 fish examined, four were found to be F(1) , 14 backcrosses to whitefish and seven backcrosses to vendace. The estimates of historical gene flow suggested higher rates from introduced vendace into native whitefish than vice versa, whereas estimates of contemporary gene flow were equal. Mitochondrial introgression was skewed, with 18 backcrosses having vendace mtDNA and only three with whitefish mtDNA. Hybrids and backcrosses had intermediate morphology and niche utilization compared with parental species. No evidence of selection against hybrids or backcrosses was apparent, as both hybrid and backcross growth rates and fecundities were high. Hybrids (F(1) ) were only detected in 2 year-classes, suggesting temporal variability in mating between vendace and whitefish. However, our data show that hybrids reached sexual maturity and reproduced actively, with backcrosses recorded from six consecutive year-classes, whereas no F(2) individuals were found. The results indicate widespread introgression, as 10.8% of coregonids were estimated to be backcrosses.

Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600-km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later-generation hybrid genotypes. Observed clines in ecologically important phenotypic traits-fur coloration and cranial morphology-were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.

Speciation genetics of biological invasions with hybridization

The negative effects of human-induced habitat disturbance and modification on multiple dimensions of biological diversity are well chronicled (Turner 1996; Harding et al. 1998; Lawton et al. 1998; Sakai et al. 2001). Among the more insidious consequences is secondary contact between formerly allopatric taxa (Anderson & Hubricht 1938; Perry et al. 2002; Seehausen 2006). How the secondary contact will play out is unpredictable (Ellstrand et al. 2010), but if the taxa are not fully reproductively isolated, hybridization is likely, and if the resulting progeny are fertile, the eventual outcome is often devastating from a conservation perspective (Rhymer & Simberloff 1996; Wolf et al. 2001; McDonald et al. 2008). In this issue of Molecular Ecology, Steeves et al. (2010) present an analysis of hybridization between two avian species, one of which is critically endangered and the other of which is invasive. Their discovery that the endangered species has not yet been hybridized to extinction is promising and not what one would necessarily expect from theory.

Hybrid fitness across time and habitats

There has been considerable debate about the role of hybrids in the evolutionary process. One question has involved the relative fitness of hybrid versus non-hybrid genotypes. For some, the assumption of lower hybrid fitness continues to be integral to their concept of species and speciation. In contrast, numerous workers have suggested that hybrid genotypes might demonstrate higher relative fitness under various environmental settings. Of particular importance in deciding between these opposing hypotheses are long-term analyses coupling ecological and genetic information. Although currently rare, such analyses have provided a test of the fitness of hybrid genotypes across generations and habitats and their role in adaptation and speciation. Here we discuss examples of these analyses applied to viruses, prokaryotes, plants and Darwin's Finches.

Retention of low-fitness genotypes over six decades of admixture between native and introduced tiger salamanders

BACKGROUND

Introductions of non-native tiger salamanders into the range of California tiger salamanders have provided a rare opportunity to study the early stages of secondary contact and hybridization. We produced first- and second-generation hybrid salamanders in the lab and measured viability among these early-generation hybrid crosses to determine the strength of the initial barrier to gene exchange. We also created contemporary-generation hybrids in the lab and evaluated the extent to which selection has affected fitness over approximately 20 generations of admixture. Additionally, we examined the inheritance of quantitative phenotypic variation to better understand how evolution has progressed since secondary contact.

RESULTS

We found significant variation in the fitness of hybrids, with non-native backcrosses experiencing the highest survival and F2 hybrids the lowest. Contemporary-generation hybrids had similar survival to that of F2 families, contrary to our expectation that 20 generations of selection in the wild would eliminate unfit genotypes and increase survival. Hybrid survival clearly exhibited effects of epistasis, whereas size and growth showed mostly additive genetic variance, and time to metamorphosis showed substantial dominance.

CONCLUSIONS

Based on first- and second- generation cross types, our results suggest that the initial barrier to gene flow between these two species was relatively weak, and subsequent evolution has been generally slow. The persistence of low-viability recombinant hybrid genotypes in some contemporary populations illustrates that while hybridization can provide a potent source of genetic variation upon which natural selection can act, the sorting of fit from unfit gene combinations might be inefficient in highly admixed populations. Spatio-temporal fluctuation in selection or complex genetics has perhaps stalled adaptive evolution in this system despite selection for admixed genotypes within generations.