BARRIERS TO SYMPATRY BETWEEN AVIAN SIBLING SPECIES (PARIDAE: BAEOLOPHUS) IN LOCAL SECONDARY CONTACT

Comparative Population Genomics of Cryptic Speciation and Adaptive Divergence in Bicknell's and Gray-Cheeked Thrushes (Aves: Catharus bicknelli and Catharus minimus)

Cryptic speciation may occur when reproductive isolation is recent or the accumulation of morphological differences between sister lineages is slowed by stabilizing selection preventing phenotypic differentiation. In North America, Bicknell's Thrush (Catharus bicknelli) and its sister species, the Gray-cheeked Thrush (Catharus minimus), are parapatrically breeding migratory songbirds, distinguishable in nature only by subtle differences in song and coloration, and were recognized as distinct species only in the 1990s. Previous molecular studies have estimated that the species diverged approximately 120,000-420,000 YBP and found very low levels of introgression despite their similarity and sympatry in the spring (prebreeding) migration. To further clarify the history, genetic divergence, genomic structure, and adaptive processes in C. bicknelli and C. minimus, we sequenced and assembled high-coverage reference genomes of both species and resequenced genomes from population samples of C. bicknelli, C. minimus, and two individuals of the Swainson's Thrush (Catharus ustulatus). The genome of C. bicknelli exhibits markedly higher abundances of transposable elements compared with other Catharus and chicken. Demographic and admixture analyses confirm moderate genome-wide differentiation (Fst ≈ 0.10) and limited gene flow between C. bicknelli and C. minimus, but suggest a more recent divergence than estimates based on mtDNA. We find evidence of rapid evolution of the Z-chromosome and elevated divergence consistent with natural selection on genomic regions near genes involved with neuronal processes in C. bicknelli. These genomes are a useful resource for future investigations of speciation, migration, and adaptation in Catharus thrushes.

Environmental (in)dependence of a hybrid zone: Insights from molecular markers and ecological niche modeling in a hybrid zone of Origanum (Lamiaceae) on the island of Crete

The role of environment and the relative significance of endogenous versus exogenous selection in shaping hybrid zones have been crucial issues in the studies of hybridization. Recent advances in ecological niche modeling (ENM) offer new methodological tools, especially in combination with the genotyping of individuals in the hybrid zone. Here, we study the hybrid zone between the widely known spices Origanum onites and Origanum vulgare ssp. hirtum in Crete. We analyze the genetic structure of both parental taxa and their hybrid Origanum × intercendens using AFLP markers on 15 sympatric and 12 allopatric populations and employ ecological niche modeling and niche similarity tests to study their niche patterns. We complement these analyses with seed viability measurements. Our study revealed that the hybridizing taxa O. onites and O. vulgare ssp. hirtum and the resulting genotypic classes showed geographical and environmental niche similarities based on the predictions of ENMs and the subsequent similarity tests. The occurrence of the hybrid zone is not directly dependent on environmental factors which favor the fitness of the hybrid compared to the parental taxa, but rather on aspects such as historical factors and management practices, which may contribute to the localization and maintenance of the contact zone between parental species. Our results suggest that if a minimum required niche differentiation between genotypic classes is not achieved, environmental dependence might not have a prominent role on the outcome of the hybridization.

Potential Implications of Climate Change on Aegilops Species Distribution: Sympatry of These Crop Wild Relatives with the Major European Crop Triticum aestivum and Conservation Issues

Gene flow from crop to wild relatives is a common phenomenon which can lead to reduced adaptation of the wild relatives to natural ecosystems and/or increased adaptation to agrosystems (weediness). With global warming, wild relative distributions will likely change, thus modifying the width and/or location of co-occurrence zones where crop-wild hybridization events could occur (sympatry). This study investigates current and 2050 projected changes in sympatry levels between cultivated wheat and six of the most common Aegilops species in Europe. Projections were generated using MaxEnt on presence-only data, bioclimatic variables, and considering two migration hypotheses and two 2050 climate scenarios (RCP4.5 and RCP8.5). Overall, a general decline in suitable climatic conditions for Aegilops species outside the European zone and a parallel increase in Europe were predicted. If no migration could occur, the decline was predicted to be more acute outside than within the European zone. The potential sympatry level in Europe by 2050 was predicted to increase at a higher rate than species richness, and most expansions were predicted to occur in three countries, which are currently among the top four wheat producers in Europe: Russia, France and Ukraine. The results are also discussed with regard to conservation issues of these crop wild relatives.

Differential effects of climate and species interactions on range limits at a hybrid zone: potential direct and indirect impacts of climate change

The relative contributions of climate versus interspecific interactions in shaping species distributions have important implications for closely related species at contact zones. When hybridization occurs within a contact zone, these factors regulate hybrid zone location and movement. While a hybrid zone's position may depend on both climate and interactions between the hybridizing species, little is known about how these factors interact to affect hybrid zone dynamics. Here, we utilize SDM (species distribution modeling) both to characterize the factors affecting the current location of a moving North American avian hybrid zone and to predict potential direct and indirect effects of climate change on future distributions. We focus on two passerine species that hybridize where their ranges meet, the Black‐capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee. Our contemporary climate models predict the occurrence of climatically suitable habitat extending beyond the hybrid zone for P. atricapillus only, suggesting that interspecific interactions primarily regulate this range boundary in P. atricapillus, while climatic factors regulate P. carolinensis. Year 2050 climate models predict a drastic northward shift in suitable habitat for P. carolinensis. Because of the greater importance of interspecific interactions for regulating the southern range limit of P. atricapillus, these climate‐mediated shifts in the distribution of P. carolinensis may indirectly lead to a range retraction in P. atricapillus. Together, our results highlight the ways climate change can both directly and indirectly affect species distributions and hybrid zone location. In addition, our study lends support to the longstanding hypothesis that abiotic factors regulate species' poleward range limits, while biotic factors shape equatorial range limits.

Patterns of differential introgression in a salamander hybrid zone: inferences from genetic data and ecological niche modelling

Hybrid zones have yielded considerable insight into many evolutionary processes, including speciation and the maintenance of species boundaries. Presented here are analyses from a hybrid zone that occurs among three salamanders -Plethodon jordani, Plethodon metcalfi and Plethodon teyahalee- from the southern Appalachian Mountains. Using a novel statistical approach for analysis of non-clinal, multispecies hybrid zones, we examined spatial patterns of variation at four markers: single-nucleotide polymorphisms (SNPs) located in the mtDNA ND2 gene and the nuclear DNA ILF3 gene, and the morphological markers of red cheek pigmentation and white flecks. Concordance of the ILF3 marker and both morphological markers across four transects is observed. In three of the four transects, however, the pattern of mtDNA is discordant from all other markers, with a higher representation of P. metcalfi mtDNA in the northern and lower elevation localities than is expected given the ILF3 marker and morphology. To explore whether climate plays a role in the position of the hybrid zone, we created ecological niche models for P. jordani and P. metcalfi. Modelling results suggest that hybrid zone position is not determined by steep gradients in climatic suitability for either species. Instead, the hybrid zone lies in a climatically homogenous region that is broadly suitable for both P. jordani and P. metcalfi. We discuss various selective (natural selection associated with climate) and behavioural processes (sex-biased dispersal, asymmetric reproductive isolation) that might explain the discordance in the extent to which mtDNA and nuclear DNA and colour-pattern traits have moved across this hybrid zone.

Quantifying ecological, morphological, and genetic variation to delimit species in the coast horned lizard species complex (Phrynosoma)

Lineage separation and divergence form a temporally extended process whereby populations may diverge genetically, morphologically, or ecologically, and these contingent properties of species provide the operational criteria necessary for species delimitation. We inferred the historical process of lineage formation in the coast horned lizard (Phrynosoma coronatum) species complex by evaluating a diversity of operational species criteria, including divergence in mtDNA (98 specimens; 2,781 bp) and nuclear loci (RAG-1, 1,054 bp; BDNF 529 bp), ecological niches (11 bioclimatic variables; 285 unique localities), and cranial horn shapes (493 specimens; 16 landmarks). A phylogenetic analysis of mtDNA recovers 5 phylogeographic groups arranged latitudinally along the Baja California Peninsula and in California. The 2 southern phylogeographic groups exhibit concordance between genetic, morphological, and ecological divergence; however, differentiation is weak or absent at more recent levels defined by phylogeographic breaks in California. Interpreting these operational species criteria together suggests that there are 3 ecologically divergent and morphologically diagnosable species within the P. coronatum complex. Our 3-species taxonomic hypothesis invokes a deep coalescence event when fitting the mtDNA genealogy into the species tree, which is not unexpected for populations that have diverged recently. Although the hypothesis that the 3 phylogeographic groups distributed across California each represent distinctive species is not supported by all of the operational species criteria evaluated in this study, the conservation status of the imperiled populations represented by these genealogical units remains critical.

Identification and dynamics of a cryptic suture zone in tropical rainforest

Suture zones, shared regions of secondary contact between long-isolated lineages, are natural laboratories for studying divergence and speciation. For tropical rainforest, the existence of suture zones and their significance for speciation has been controversial. Using comparative phylogeographic evidence, we locate a morphologically cryptic suture zone in the Australian Wet Tropics rainforest. Fourteen out of 18 contacts involve morphologically cryptic phylogeographic lineages, with mtDNA sequence divergences ranging from 2 to 15 per cent. Contact zones are significantly clustered in a suture zone located between two major Quaternary refugia. Within this area, there is a trend for secondary contacts to occur in regions with low environmental suitability relative to both adjacent refugia and, by inference, the parental lineages. The extent and form of reproductive isolation among interacting lineages varies across species, ranging from random admixture to speciation, in one case via reinforcement. Comparative phylogeographic studies, combined with environmental analysis at a fine-scale and across varying climates, can generate new insights into suture zone formation and to diversification processes in species-rich tropical rainforests. As arenas for evolutionary experimentation, suture zones merit special attention for conservation.

Phylogeographic history and gene flow among giant Galápagos tortoises on southern Isabela Island

Volcanic islands represent excellent models with which to study the effect of vicariance on colonization and dispersal, particularly when the evolution of genetic diversity mirrors the sequence of geological events that led to island formation. Phylogeographic inference, however, can be particularly challenging for recent dispersal events within islands, where the antagonistic effects of land bridge formation and vicariance can affect movements of organisms with limited dispersal ability. We investigated levels of genetic divergence and recovered signatures of dispersal events for 631 Galápagos giant tortoises across the volcanoes of Sierra Negra and Cerro Azul on the island of Isabela. These volcanoes are among the most recent formations in the Galápagos (<0.7 million years), and previous studies based on genetic and morphological data could not recover a consistent pattern of lineage sorting. We integrated nested clade analysis of mitochondrial DNA control region sequences, to infer historical patterns of colonization, and a novel Bayesian multilocus genotyping method for recovering evidence of recent migration across volcanoes using eleven microsatellite loci. These genetic studies illuminate taxonomic distinctions as well as provide guidance to possible repatriation programs aimed at countering the rapid population declines of these spectacular animals.