Rapid fixation of non-native alleles revealed by genome-wide SNP analysis of hybrid tiger salamanders

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.

An Account of the Taxonomy of North American Wolves From Morphological and Genetic Analyses

The available scientific literature was reviewed to assess the taxonomic standing of North American wolves, including subspecies of the gray wolf, Canis lupus. The recent scientific proposal that the eastern wolf, C. l. lycaon, is not a subspecies of gray wolf, but a full species, Canis lycaon, is well-supported by both morphological and genetic data. This species' range extends westward to Minnesota, and it hybridizes with gray wolves where the two species are in contact in eastern Canada and the Upper Peninsula of Michigan, Wisconsin, and Minnesota. Genetic data support a close relationship between eastern wolf and red wolf Canis rufus, but do not support the proposal that they are the same species; it is more likely that they evolved independently from different lineages of a common ancestor with coyotes. The genetic distinctiveness of the Mexican wolf Canis lupus baileyi supports its recognition as a subspecies. The available genetic and morphometric data do not provide clear support for the recognition of the Arctic wolf Canis lupus arctos, but the available genetic data are almost entirely limited to one group of genetic markers (microsatellite DNA) and are not definitive on this question. Recognition of the northern timber wolf Canis lupus occidentalis and the plains wolf Canis lupus nubilus as subspecies is supported by morphological data and extensive studies of microsatellite DNA variation where both subspecies are in contact in Canada. The wolves of coastal areas in southeastern Alaska and British Columbia should be assigned to C. lupus nubilus. There is scientific support for the taxa recognized here, but delineation of exact geographic boundaries presents challenges. Rather than sharp boundaries between taxa, boundaries should generally be thought of as intergrade zones of variable width.

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).

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.

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci

BACKGROUND

Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion.

RESULTS

Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species.

CONCLUSION

This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.

Hybrid fitness, adaptation and evolutionary diversification: lessons learned from Louisiana Irises

Estimates of hybrid fitness have been used as either a platform for testing the potential role of natural hybridization in the evolution of species and species complexes or, alternatively, as a rationale for dismissing hybridization events as being of any evolutionary significance. From the time of Darwin's publication of The Origin, through the neo-Darwinian synthesis, to the present day, the observation of variability in hybrid fitness has remained a challenge for some models of speciation. Yet, Darwin and others have reported the elevated fitness of hybrid genotypes under certain environmental conditions. In modern scientific terminology, this observation reflects the fact that hybrid genotypes can demonstrate genotype × environment interactions. In the current review, we illustrate the development of one plant species complex, namely the Louisiana Irises, into a 'model system' for investigating hybrid fitness and the role of genetic exchange in adaptive evolution and diversification. In particular, we will argue that a multitude of approaches, involving both experimental and natural environments, and incorporating both manipulative analyses and surveys of natural populations, are necessary to adequately test for the evolutionary significance of introgressive hybridization. An appreciation of the variability of hybrid fitness leads to the conclusion that certain genetic signatures reflect adaptive evolution. Furthermore, tests of the frequency of allopatric versus sympatric/parapatric divergence (that is, divergence with ongoing gene flow) support hybrid genotypes as a mechanism of evolutionary diversification in numerous species complexes.

The use of genetics for the management of a recovering population: temporal assessment of migratory peregrine falcons in North America

BACKGROUND

Our ability to monitor populations or species that were once threatened or endangered and in the process of recovery is enhanced by using genetic methods to assess overall population stability and size over time. This can be accomplished most directly by obtaining genetic measures from temporally-spaced samples that reflect the overall stability of the population as given by changes in genetic diversity levels (allelic richness and heterozygosity), degree of population differentiation (F(ST) and D(EST)), and effective population size (N(e)). The primary goal of any recovery effort is to produce a long-term self-sustaining population, and these genetic measures provide a metric by which we can gauge our progress and help make important management decisions.

METHODOLOGY/PRINCIPAL FINDINGS

The peregrine falcon in North America (Falco peregrinus tundrius and anatum) was delisted in 1994 and 1999, respectively, and its abundance will be monitored by the species Recovery Team every three years until 2015. Although the United States Fish and Wildlife Service makes a distinction between tundrius and anatum subspecies, our genetic results based on eleven microsatellite loci suggest limited differentiation that can be attributed to an isolation by distance relationship and warrant no delineation of these two subspecies in its northern latitudinal distribution from Alaska through Canada into Greenland. Using temporal samples collected at Padre Island, Texas during migration (seven temporal time periods between 1985-2007), no significant differences in genetic diversity or significant population differentiation in allele frequencies between time periods were observed and were indistinguishable from those obtained from tundrius/anatum breeding locations throughout their northern distribution. Estimates of harmonic mean N(e) were variable and imprecise, but always greater than 500 when employing multiple temporal genetic methods.

CONCLUSIONS/SIGNIFICANCE

These results, including those from simulations to assess the power of each method to estimate N(e), suggest a stable or growing population, which is consistent with ongoing field-based monitoring surveys. Therefore, historic and continuing efforts to prevent the extinction of the peregrine falcon in North America appear successful with no indication of recent decline, at least from the northern latitude range-wide perspective. The results also further highlight the importance of archiving samples and their use for continual assessment of population recovery and long-term viability.

Birth of a hotspot of intraspecific genetic diversity: notes from the underground

Hotspots of intraspecific diversity have been observed in most species, often within areas of putative Pleistocene refugia. They have thus mostly been viewed as the outcome of prolonged stability of large populations within the refugia. However, recent evidence has suggested that several other microevolutionary processes could also be involved in their formation. Here, we investigate the contribution of these processes to current range-wide patterns of genetic diversity in the Italian endemic mole Talpa romana, using both nuclear (30 allozyme loci) and mitochondrial markers (cytochrome b sequences). Southern populations of this species showed an allozyme variation that is amongst the highest observed in small mammals (most populations had an expected heterozygosity of 0.10 or above), which was particularly unexpected for a subterranean species. Population genetic, phylogeographic and historical demographic analyses indicated that T. romana populations repeatedly underwent allopatric differentiations followed by secondary admixture within the refugial range in southern Italy. A prolonged demographic stability was reliably inferred from the mitochondrial DNA data only for a population group located north and east of the Calabrian peninsula, showing comparatively lower levels of allozyme variability, and lacking evidence of secondary admixture with other groups. Thus, our results point to the admixture between differentiated lineages as the main cause of the higher levels of diversity of refugial populations. When compared with the Pleistocene evolutionary history recently inferred for species from both the same and other geographic regions, these results suggest the need for a reappraisal of the role of gene exchange in the formation of intraspecific hotspots of genetic diversity.

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.

Rapid spread of invasive genes into a threatened native species

When introduced or cultivated plants or animals hybridize with their native relatives, the spread of invasive genes into native populations might have biological, aesthetic, and legal implications. Models suggest that the rate of displacement of native by invasive alleles can be rapid and inevitable if they are favored by natural selection. We document the spread of a few introduced genes 90 km into a threatened native species (the California Tiger Salamander) in 60 years. Meanwhile, a majority of genetic markers (65 of 68) show little evidence of spread beyond the region where introductions occurred. Using computer simulations, we found that such a pattern is unlikely to emerge by chance among selectively neutral markers. Therefore, our results imply that natural selection has favored both the movement and fixation of these exceptional invasive alleles. The legal status of introgressed populations (native populations that are slightly genetically modified) is unresolved by the US Endangered Species Act. Our results illustrate that genetic and ecological factors need to be carefully weighed when considering different criteria for protection, because different rules could result in dramatically different geographic areas and numbers of individuals being protected.

Adaptation by introgression

Both selective and random processes can affect the outcome of natural hybridization. A recent analysis in BMC Evolutionary Biology of natural hybridization between an introduced and a native salamander reveals the mosaic nature of introgression, which is probably caused by a combination of selection and demography.