THE EXTENT OF INTROGRESSION OUTSIDE THE CONTACT ZONE BETWEEN NOTROPIS CORNUTUS AND NOTROPIS CHRYSOCEPHALUS (TELEOSTEI: CYPRINIDAE)

The cyprinid fishes, Notropis cornutus and N. chrysocephalus, hybridize in a long, narrow zone in the midwestern United States. To quantify the extent of introgression of genetic markers outside of this zone, samples were collected along transects starting near the region of contact (as defined by morphological characters), followed by samples progressively more distant. Diagnostic allozymic and mitochondrial DNA (mtDNA) restriction site markers were used to estimate the extent of introgression outside of the zone, while polymorphic allozyme and mtDNA markers were used to evaluate the potential for gene flow among populations within transects. Analysis of populations from the northern transect provided evidence for differentiation of populations for some of the markers; however, on average, enough gene flow has occurred to overcome substantial differentiation. Introgressed mtDNA and allozyme haplotypes were rare and found only in the population closest to the contact zone. The rarity of introgressed alleles in the more northern populations is consistent with the recent origin of these populations after the Wisconsin glaciation (less than 12,000 years bp) and/or selection maintaining the northern boundary of the contact zone. Analysis of populations from the southern transect revealed evidence for population subdivision but no evidence for introgression at the diagnostic allozyme loci; however, nearly all individuals from this transect possessed introgressed mtDNA haplotypes, with samples furthest from the contact zone exhibiting the highest frequencies of introgression. Patterns of variation for one of the polymorphic allozyme markers (Est-A) and introgressed mtDNAs were highly correlated, suggesting that allozymic heterogeneity at this locus is also the result of introgression. The most likely explanation for these data is that these introgressed haplotypes are indicators of a more southern position of the contact zone during the Pleistocene, with the contact zone shifting northward with the recession of the glacial front. Such movement implicates selection in the maintenance of distributional limits of these species, and hence, the width and position of the contact zone.     

PATTERNS OF HYBRIDIZATION IN THE PIRIQUETA CAROLINIANA COMPLEX IN CENTRAL FLORIDA: EVIDENCE FOR AN EXPANDING HYBRID ZONE

We have identified a broad zone of hybridization between two morphologically and ecologically distinct herbaceous perennial taxa (morphotypes) within the Piriqueta caroliniana complex, which extends more than 300 km across the central Florida peninsula. Phylogeographic analyses indicate that the caroliniana morphotype has been present in north and central Florida since the early Pleistocene and that the viridis morphotype has immigrated into southern Florida much more recently. We examine the distribution of diagnostic morphological characters and nuclear genetic markers to assess the extent and patterns of introgression in this system. The results are consistent with the hypothesis that this hybrid zone has expanded north in recent history as viridis alleles have introgressed into regions that were previously occupied by populations of caroliniana. Genetic markers diagnostic for caroliniana have consistently high frequencies across the hybrid zone, whereas markers for viridis are extremely variable among populations with frequency reversals in adjacent populations. The latter pattern is probably the result of the combined stochastic effects of dispersal and drift on viridis alleles as they introgressed northward. Additional evidence for the recent expansion of this hybrid zone comes from patterns of variation for morphological and genetic markers. As expected for an expanding hybrid zone, within-population morphological variation was greatest toward the advancing front of introgression and levels of genetic variation for neutral diagnostic markers were greatest in the region of initial contact and lower in areas of recent expansion. The observed patterns of variation suggest that at least some hybrid genotypes have high fitnesses, which has led to the expansion of the hybrid zone via the displacement of parental genotypes in central Florida.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Large-scale asymmetric introgression of cytoplasmic DNA reveals Holocene range displacement in a North American boreal pine complex

Jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta var. latifolia) are two North American boreal hard pines that hybridize in their zone of contact in western Canada. The main objective of this study was to characterize their patterns of introgression resulting from past and recent gene flow, using cytoplasmic markers having maternal or paternal inheritance. Mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) diversity was assessed in allopatric populations of each species and in stands from the current zone of contact containing morphological hybrids. Cluster analyses were used to identify genetic discontinuities among groups of populations. A canonical analysis was also conducted to detect putative associations among cytoplasmic DNA variation, tree morphology, and site ecological features. MtDNA introgression was extensive and asymmetric: it was detected in P. banksiana populations from the hybrid zone and from allopatric areas, but not in P. contorta populations. Very weak cpDNA introgression was observed, and only in P. banksiana populations. The mtDNA introgression pattern indicated that central Canada was first colonized by migrants from a P. contorta glacial population located west of the Rocky Mountains, before being replaced by P. banksiana migrating westward during the Holocene. In contrast, extensive pollen gene flow would have erased the cpDNA traces of this ancient presence of P. contorta. Additional evidence for this process was provided by the results of canonical analysis, which indicated that the current cpDNA background of trees reflected recent pollen gene flow from the surrounding dominant species rather than historical events that took place during the postglacial colonization.     

Detecting signatures of inter‐regional and inter‐specific hybridization among the Chinese rhesus macaque specific pathogen‐free (SPF) population using single nucleotide polymorphic (SNP) markers

Background While rates of gene flow between rhesus and longtail macaque populations near their hybrid zone in Indochina have been quantified elsewhere, this study demonstrates that the inter‐specific introgression is not limited to the Indochinese hybrid zone but is more geographically widespread. Methods Twelve rhesus and longtail macaque populations were analyzed using single nucleotide polymorphic (SNP) loci. Results There is evidence for inter‐specific admixture between Chinese rhesus and mainland longtails, with implications for genetic diversity both in the Chinese super‐SPF population at the California National Primate Research Center and in other primate facilities. Eastern Chinese rhesus appeared more highly derived than western Chinese rhesus, and allele sharing between longtails and Chinese rhesus was not random with regard to geographic distance, but no significant nuclear genetic differences between eastern and western Chinese rhesus were detected among the 245 genic SNPs assayed. Conclusion The implications of this inter‐specific admixture for the use of Chinese rhesus and mainland longtail in biomedical research should be considered.     

Evaluation of introgressive hybridization among Cervidae in Japan's Kinki District via two novel genetic markers developed from public NGS data

Hybridization and backcrossing of native populations with introduced species can lead to introgression and genetic alteration. In this study, we evaluated introgression in 43 deer from a potential hybrid zone around Okinoshima Island, Kinki District, Japan. This region witnessed the migration of a hybrid population (cross between the Formosan sika deer [Cervus nippon taiouanus] and other deer species) that could potentially breed with the native Japanese sika deer (C. n. centralis). We used an existing genetic marker for the mitochondrial cytochrome b gene and two novel markers for nuclear DNA, developed using publicly available next‐generation sequencing data. We identified one mainland deer with a mitochondrial haplotype identical to that of the Formosan sika deer as well as nuclear heterozygous sequences identical to those of Formosan and Japanese sika deer. This suggests that the mainland deer is a hybrid offspring of the Okinoshima population and native deer. However, only Japanese sika deer sequences were found in the other 42 samples, indicating limited introgression. Nevertheless, hybridization pre‐ and postintroduction in the Okinoshima population could cause multispecies introgression among Japanese sika deer, negatively affecting genetic integrity. We developed a simple test based on polymerase chain reaction–restriction fragment length polymorphism to detect introgression in natural populations. Our method can accelerate genetic monitoring of Japanese sika deer in Kinki District. In conclusion, to prevent further introgression and maintain genetic integrity of Japanese sika deer, we recommend establishing fences around Okinoshima Island to limit migration, besides a continued genetic monitoring of the native deer.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

Phylogeography and the maternal origin of the tetraploid form of the Japanese spined loach, <Emphasis Type="Italic">Cobitis biwae</Emphasis>, revealed by mitochondrial DNA analysis

The Japanese spined loach Cobitis biwae includes the tetraploid form, which has 96 chromosomes, as well as the diploid form, which has 48 chromosomes. In the present study, we analyzed the mitochondrial DNA (mtDNA) to examine the genetic relationships among 82 populations of diploid–tetraploid complexes of C. biwae. Restriction fragment length polymorphism (RFLP) analysis of the ND1 region on 202 individuals revealed that C. biwae contains three genetically divergent major groups corresponding to geographical proximities (the Eastern, Western, and Kochi groups). Phylogenetic analyses (neighbor-joining, NJ and maximum-parsimony, MP) of a part of the cytochrome b gene sequence (748bp) in 31 individuals supported the three major groups recognized by RFLP, and indicated considerable genetic differentiation between the Western group and the other two groups (average, 15.2%). The genetic relationship and distribution pattern of the three major groups hypothesized two major dispersions of C. biwae during the middle Miocene: first, the ancestor of the Eastern and Kochi groups had spread from the West through wide regions of present-day Honshu and Shikoku Islands, and following that period, the ancestor of the Western group with a different mtDNA composition, probably resulting from the mtDNA introgression from C. takatsuensis, moved into western Japan and pushed the former north eastward. All tetraploid form populations were included in the Western group and treated as a monophyletic cluster with low genetic divergence. It is notable that two diploid populations geographically adjacent to the tetraploid range were genetically closely related to the tetraploid forms. This result suggests that these diploid populations were directly related to the maternal origin of the tetraploid form.     

Phylogeography of the mountain chickadee (Poecile gambeli): diversification, introgression, and expansion in response to Quaternary climate change

Since the late 1990s, molecular techniques have fuelled debate about the role of Pleistocene glacial cycles in structuring contemporary avian diversity in North America. The debate is still heated; however, there is widespread agreement that the Pleistocene glacial cycles forced the repeated contraction, fragmentation, and expansion of the North American biota. These demographic processes should leave genetic 'footprints' in modern descendants, suggesting that detailed population genetic studies of contemporary species provide the key to elucidating the impact of the late Quaternary (late Pleistocene-Holocene). We present an analysis of mitochondrial DNA (mtDNA) variation in the mountain chickadee (Poecile gambeli) in an attempt to examine the genetic evidence of the impact of the late Quaternary glacial cycles. Phylogenetic analyses reveal two strongly supported clades of P. gambeli: an Eastern Clade (Rocky Mountains and Great Basin) and a Western Clade (Sierra Nevada and Cascades). Post-glacial introgression is apparent between these two clades in the Mono Lake region of Central California. Within the Eastern Clade there is evidence of isolation-by-distance in the Rocky Mountain populations, and of limited gene flow into and around the Great Basin. Coalescent analysis of genetic variation in the Western Clade indicates that northern (Sierra Nevada/Cascades) and southern (Transverse/Peninsular Ranges) populations have been isolated and evolving independently for nearly 60,000 years.     

Evidences of introgression from cultivated rice to Oryza rufipogon (Poaceae) populations based on SSR fingerprinting: implications for wild rice differentiation and conservation

Crop-to-wild introgression may play an important role in evolution of wild species. Asian cultivated rice (Oryza sativa L.) is of a particular concern because of its cross-compatibility with the wild ancestor, O. rufipogon Griff. The distribution of cultivated rice and O. rufipogon populations is extensively sympatric, particularly in Asia where many wild populations are surrounded by rice fields. Consequently, gene flow from cultivated rice may have a potential to alter genetic composition of wild rice populations in close proximity. In this study, we estimated introgression of cultivated rice with O. rufipogon based on analyses of 139 rice varieties (86 indica and 53 japonica ecotypes) and 336 wild individuals from 11 O. rufipogon populations in China. DNA fingerprinting based on 17 selected rice simple sequence repeat (SSR) primer pairs was adopted to measure allelic frequencies in rice varieties and O. rufipogon samples, and to estimate genetic associations between wild and cultivated rice through cluster analysis. We detected consanguinity of cultivated rice in O. rufipogon populations according to the admixture model of the STRUCTURE program. The analyses showedz that four wild rice populations, DX-P1, DX-P2, GZ-P2, and HL-P, contained some rare alleles that were commonly found in the rice varieties examined. In addition, the four wild rice populations that scattered among the rice varieties in the cluster analysis showed a closer affinity to the cultivars than the other wild populations. This finding supports the contention of substantial gene flow from crop to wild species when these species occur close to each other. The introgressive populations had slightly higher genetic diversity than those that were isolated from rice. Crop-to-wild introgression may have accumulative impacts on genetic variations in wild populations, leading to significant differentiation in wild species. Therefore, effective measure should be taken to avoid considerable introgression from cultivated rice, which may influence the effective in-situ conservation of wild rice species.     

Genetic structure in the Mediterranean seagrass Posidonia oceanica: disentangling past vicariance events from contemporary patterns of gene flow

The Mediterranean Sea is a two‐basin system, with the boundary zone restricted to the Strait of Sicily and the narrow Strait of Messina. Two main population groups are recognized in the Mediterranean endemic seagrass Posidonia oceanica, corresponding to the Western and the Eastern basins. To address the nature of the East–West cleavage in P. oceanica, the main aims of this study were: (i) to define the genetic structure within the potential contact zone (i.e. the Strait of Sicily) and clarify the extent of gene flow between the two population groups, and (ii) to investigate the role of present water circulation patterns vs. past evolutionary events on the observed genetic pattern. To achieve these goals, we utilized SSR markers and we simulated, with respect to current regime, the possible present‐day dispersal pattern of Posidonia floating fruits using 28‐day numerical Lagrangian trajectories. The results obtained confirm the presence of the two main population groups, without any indices of reproductive isolation, with the break zone located at the level of the Southern tip of Calabria. The populations in the Strait of Sicily showed higher affinity with Western than with Eastern populations. This pattern of genetic structure probably reflects historical avenues of recolonization from relict glacial areas and past vicariance events, but seems to persist as a result of the low connectivity among populations via marine currents, as suggested by our dispersal simulation analysis.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Molecular genetic variation of boll weevil populations in North America estimated with microsatellites: Implications for patterns of dispersal

The boll weevil (Anthonomus grandis Boheman) is an insect pest of cotton that underwent a well-documented range expansion across the southeastern U.S. from Mexico beginning about 110 years ago. Eleven microsatellite loci were surveyed to infer the magnitude and pattern of genetic differentiation among boll weevil populations from 18 locations across eight U.S. states and northeast Mexico. Estimates of genetic diversity (allelic diversity and heterozygosity) were greater in Southern than Northern populations, and were greater in the west than the east among Northern populations. Boll weevil populations were genetically structured as a whole across the geographic range sampled, with a global F _ST of 0.241. South-central populations exhibit classic isolation by distance, but evidence suggests that populations within the Eastern and Western regions have not yet reached genetic equilibrium. Gene flow appears to be relatively high among populations within the Eastern region. Population assignment data and estimates of gene flow indicate that migration between locations separated by < 300 km is frequent. The database of microsatellite genotypes generated in this study now makes it possible, through population assignment techniques, to identify the most likely geographic source of a boll weevil reintroduced to an eradication zone, which will help action agencies decide the most appropriate mitigation response. Kyung Seok Kim - Formerly: USDA-ARS, Areawide Pest Management Research Unit, 2771 F & B Rd., College Station,TX 77845, USA.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

A genetically distinct hybrid zone occurs for two globally invasive mosquito fish species with striking phenotypic resemblance

Hybrid zones allow for the investigation of incipient speciation and related evolutionary processes of selection, gene flow, and migration. Interspecific dynamics, like competition, can impact the size, shape, and directional movement of species in hybrid zones. Hybrid zones contribute to a paradox for the biological species concept because interbreeding between species occurs while parental forms remain distinct. A long‐standing zone of intergradation or introgression exists for eastern and western mosquito fish (Gambusia holbrooki and G. affinis) around Mobile Bay, AL. The region has been studied episodically, over decades, making it perfect for addressing temporal dynamics and for providing a deeper understanding of the genetics of these periodically reclassified fishes (as species or subspecies). We used six microsatellite markers to assess the current population structure and gene flow patterns across 19 populations of mosquito fish and then compared our results with historical data. Genetic evidence demonstrates that the current hybrid zone is located in a similar geographic region as the historical one, even after three decades. Hybrid fish, however, demonstrate relatively low heterozygosity and are genetically distinct from western and eastern mosquito fish populations. Fin ray counts, sometimes used to distinguish the two species from one another, demonstrate more eastern (G. holbrooki) phenotype fish within the molecular genetic hybrid zone today. Mosquito fish are globally invasive, often found on the leading edge of flooded waters that they colonize, so the impact of hurricanes in the wake of climate change was also evaluated. An increase in the frequency and intensity of hurricanes in the hybrid region has occurred, and this point warrants further attention since hurricanes are known to move these aggressive, invasive species into novel territory. This work contributes to our classical understanding of hybrid zone temporal dynamics, refines our understanding of mosquito fish genetics in their native range, evaluates important genotype–phenotype relationships, and identifies a potential new impact of climate change.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

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

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

PATTERNS OF VARIATION AND LINKAGE DISEQUILIBRIUM IN A FIELD CRICKET HYBRID ZONE

The distribution of multilocus genotypes found within a natural hybrid zone is determined by the sample of genotypes present when the hybrid zone first formed, by subsequent patterns of genetic exchange between the hybridizing taxa, and by drift and selection within each of the hybrid zone populations. We have used anonymous nuclear DNA restriction fragment polymorphisms (RFLPs) to characterize the array of multilocus genotypes present within a well-studied hybrid zone between two eastern North American field crickets, Gryllus pennsylvanicus and Gryllus firmus. These crickets hybridize along a zone of contact that extends from New England to Virginia. Previous studies have shown that both premating and postmating barriers exist between the two cricket species, but the absence of diagnostic morphological and allozyme markers has made it difficult to assess the consequences of these barriers for genetic exchange. Analyses based on four diagnostic anonymous nuclear markers indicate that hybrid zone populations in Connecticut contain few F₁ hybrids, and that nonrandom associations persist among nuclear gene markers, between nuclear and cytoplasmic markers, and between molecular markers and morphology. Field cricket populations within the hybrid zone are not “hybrid swarms” but consist primarily of crickets that are very much like one or the other of the parental species. Despite ample opportunity for genetic exchange and evidence for introgression at some loci, the two species remain quite distinct. Such a pattern appears to be characteristic of many natural hybrid zones.