Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species.

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Across a phylogeographic break in the Qinling Mountains-Huaihe River Line: Quaternary evolutionary history of a medicinal and edible homologous plant (Allium macrostemon) in China.

Biogeographic barriers to gene flow are central to studies of plant phylogeography. There are many physical and geographic barriers in China, but few studies have used molecular ecological evidence to investigate the natural geographic isolation barrier of the Qinling Mountains-Huaihe River Line (QHL). Allium macrostemon is a precious Chinese perennial herb belonging to the Amaryllidaceae family. It is used as a food and medicine, with a variety of health and healing properties. Five SSR markers, three chloroplast DNA (cpDNA) markers (psbA-trnH, rps16 and trnL-F), one nuclear ribosomal DNA (nrDNA) marker (ITS), and simplified genome GBS sequencing were used to analyse the genetic diversity and structure of A. macrostemon. Combining SSR, cpDNA, nrDNA ITS data and GBS analysis results, we divided A. macrostemon populations into northern and southern groups, with the southern group further divided into southwestern and central-southeastern groups. Niche simulation results reveal that the distribution area of A. macrostemon will reach its maximum in the future. These data indicate that the regional separation of A. macrostemon has been maintained by the combined influence of a geographical barrier and Quaternary climate, and that the back-and-forth fluctuations of QHL and Quaternary climate have played an important role in this process. QHL acts as a north-south dividing line in phylogeography and population genetic structure, promoting physical geographic isolation. This study provides a theoretical basis for the conservation, development, and utilization of A. macrostemon resources. It further provides a reference for understanding the systematic geographical pattern of the large-scale spatial distribution of plants in China and enriches our understanding of Quaternary plant evolution in areas with complex terrain.

The population genomics of Conyza spp. in soybean macroregions suggest the spread of herbicide resistance through intraspecific and interspecific gene flow.

Herbicide-resistant Conyza spp. are a threat to many crops. These widespread weeds are closely related species and often cooccur. To characterize the origins of their resistance and the mechanisms underlying their spread, we assessed the genomic variation in glyphosate-resistant Conyza spp. in Brazil. Twenty populations were sampled from soybean fields across four macroregions (MRSs). A genotyping-by-sequencing study resulted in 2,998 single-nucleotide polymorphisms (SNPs) obtained for C. bonariensis (L.) and the closely related C. sumatrensis (Retz) E. Walker. Higher genomic diversity (π) and heterozygosity (HO/HE) and lower inbreeding coefficient (FIS) values were detected in populations of Conyza spp. from MRS 1 (southern) than in those from other MRSs. Strong genomic structure clustered individuals into three groups (FST = 0.22; p value = 0.000) associated with the MRSs. Thus, resistance to glyphosate originated from independent selection in different MRSs across Brazil. Our dataset supports the occurrence of intraspecific gene flow in Brazil and identified individuals of C. bonariensis that did not group within species. These findings suggest that allelic introgressions within and among species have impacted the evolution and spread of resistance to glyphosate in Conyza spp. We discuss how to mitigate new resistance cases, particularly for the released stacked traits herbicide tolerance in soybeans.

A consequential one-night stand: Episodic historical hybridization leads to mitochondrial takeover in sympatric desert ant-eating spiders.

Disentangling the genomic intricacies underlying speciation and the causes of discordance between sources of evidence can offer remarkable insights into evolutionary dynamics. The ant-eating spider Zodarion nitidum, found across the Middle East and Egypt, displays yellowish and blackish morphs that co-occur sympatrically. These morphs additionally differ in behavioral and physiological features and show complete pre-mating reproductive isolation. In contrast, they possess similar sexual features and lack distinct differences in their mitochondrial DNA. We analyzed both Z. nitidum morphs and outgroups using genome-wide and additional mitochondrial DNA data. The genomic evidence indicated that Yellow and Black are reciprocally independent lineages without signs of recent admixture. Interestingly, the sister group of Yellow is not Black but Z. luctuosum, a morphologically distinct species. Genomic gene flow analyses pinpointed an asymmetric nuclear introgression event, with Yellow contributing nearly 5 % of its genome to Black roughly 320,000 years ago, intriguingly aligning with the independently estimated origin of the mitochondrial DNA of Black. We conclude that the blackish and yellowish morphs of Z. nitidum are long-diverged distinct species, and that the ancient and modest genomic introgression event registered resulted in a complete mitochondrial takeover of Black by Yellow. This investigation underscores the profound long-term effects that even modest hybridization events can have on the genome of organisms. It also exemplifies the utility of phylogenetic networks for estimating historical events and how integrating independent lines of evidence can increase the reliability of such estimations.

Establishing species boundaries in Bornean geckos.

Species delimitation using mitochondrial DNA (mtDNA) remains an important and accessible approach for discovering and delimiting species. However, delimiting species with a single locus (e.g. DNA barcoding) is biased towards overestimating species diversity. The highly diverse gecko genus Cyrtodactylus is one such group where delimitation using mtDNA remains the paradigm. In this study, we use genomic data to test putative species boundaries established using mtDNA within three recognized species of Cyrtodactylus on the island of Borneo. We predict that multi-locus genomic data will estimate fewer species than mtDNA, which could have important ramifications for the species diversity within the genus. We aim to (i) investigate the correspondence between species delimitations using mtDNA and genomic data, (ii) infer species trees for each target species, and (iii) quantify gene flow and identify migration patterns to assess population connectivity. We find that species diversity is overestimated and that species boundaries differ between mtDNA and nuclear data. This underscores the value of using genomic data to reassess mtDNA-based species delimitations for taxa lacking clear species boundaries. We expect the number of recognized species within Cyrtodactylus to continue increasing, but, when possible, genomic data should be included to inform more accurate species boundaries.

Population genomics unravels a lag phase during the global fall armyworm invasion.

The time that elapsed between the initial introduction and the proliferation of an invasive species is referred to as the lag phase. The identification of the lag phase is critical for generating plans for pest management and for the prevention of biosecurity failure. However, lag phases have been identified mostly through retrospective searches of historical records. The agricultural pest fall armyworm (FAW; Spodoptera frugiperda) is native to the New World. FAW invasion was first reported from West Africa in 2016, then it spread quickly through Africa, Asia, and Oceania. Here, using population genomics approaches, we demonstrate that the FAW invasion involved an undocumented lag phase. Invasive FAW populations have negative signs of genomic Tajima's D, and invasive population-specific genetic variations have particularly decreased Tajima's D, supporting a substantial amount of time for the generation of new mutations in introduced FAW populations. Model-based diffusion approximations support the existence of a period with a cessation of gene flow between native and invasive FAW populations. Taken together, these results provide strong support for the presence of a lag phase during the FAW invasion. These results show the usefulness of using population genomics analyses to identify lag phases in biological invasions.

The Possible Earliest Allopolyploidization in Tracheophytes Revealed by Phylotranscriptomics and Morphology of Selaginellaceae.

Selaginellaceae, originated in the Carboniferous and survived the Permian-Triassic mass extinction, is the largest family of lycophyte, which is sister to other tracheophytes. It stands out from tracheophytes by exhibiting extraordinary habitat diversity and lacking polyploidization. The organelle genome-based phylogenies confirmed the monophyly of Selaginella, with six or seven subgenera grouped into two superclades, but the phylogenetic positions of the enigmatic Selaginella sanguinolenta clade remained problematic. Here, we conducted a phylogenomic study on Selaginellaceae utilizing large-scale nuclear gene data from RNA-seq to elucidate the phylogeny and explore the causes of the phylogenetic incongruence of the S. sanguinolenta clade. Our phylogenetic analyses resolved three different positions of the S. sanguinolenta clade, which were supported by the sorted three nuclear gene sets, respectively. The results from the gene flow test, species network inference, and plastome-based phylogeny congruently suggested a probable hybrid origin of the S. sanguinolenta clade involving each common ancestor of the two superclades in Selaginellaceae. The hybrid hypothesis is corroborated by the evidence from rhizophore morphology and spore micromorphology. The chromosome observation and Ks distributions further suggested hybridization accompanied by polyploidization. Divergence time estimation based on independent datasets from nuclear gene sets and plastid genome data congruently inferred that allopolyploidization occurred in the Early Triassic. To our best knowledge, the allopolyploidization in the Mesozoic reported here represents the earliest record of tracheophytes. Our study revealed a unique triad of phylogenetic positions for a hybrid-originated group with comprehensive evidence and proposed a hypothesis for retaining both parental alleles through gene conversion.

Ecological genomics in the Northern krill uncovers loci for local adaptation across ocean basins.

Krill are vital as food for many marine animals but also impacted by global warming. To learn how they and other zooplankton may adapt to a warmer world we studied local adaptation in the widespread Northern krill (Meganyctiphanes norvegica). We assemble and characterize its large genome and compare genome-scale variation among 74 specimens from the colder Atlantic Ocean and warmer Mediterranean Sea. The 19 Gb genome likely evolved through proliferation of retrotransposons, now targeted for inactivation by extensive DNA methylation, and contains many duplicated genes associated with molting and vision. Analysis of 760 million SNPs indicates extensive homogenizing gene-flow among populations. Nevertheless, we detect signatures of adaptive divergence across hundreds of genes, implicated in photoreception, circadian regulation, reproduction and thermal tolerance, indicating polygenic adaptation to light and temperature. The top gene candidate for ecological adaptation was nrf-6, a lipid transporter with a Mediterranean variant that may contribute to early spring reproduction. Such variation could become increasingly important for fitness in Atlantic stocks. Our study underscores the widespread but uneven distribution of adaptive variation, necessitating characterization of genetic variation among natural zooplankton populations to understand their adaptive potential, predict risks and support ocean conservation in the face of climate change.

Continent-wide genomic analysis of the African buffalo (Syncerus caffer).

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus).

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

Patterns of Genetic Diversity and Gene Flow Associated With an Aridity Gradient in Populations of Common Mole-rats, Cryptomys hottentotus hottentotus.

Genetic adaptation is the change of a population toward a phenotype that best fits the present ecological conditions of the environment it inhabits. As environmental conditions change, allele frequencies shift, resulting in different populations of the same species possessing genetic variation and divergent phenotypes. Cooperatively breeding common mole-rats (Cryptomys hottentotus hottentotus) inhabit environments along an aridity gradient in South Africa, which provides an opportunity for local genetic adaptations to occur. Using one mitochondrial gene (cytochrome b) and 3,540 SNP loci across the whole genome, we determined the phylogenetic relationship, population structure and genetic diversity of five populations of C. h. hottentotus located along an aridity gradient. Mitochondrial data identified population-specific clades that were less distinct in the two mesic populations, potentially indicating historical or recent gene flow, or the retention of ancestral haplotypes. Arid and semi-arid populations formed a distinct cluster from the non-arid populations. Genetic diversity and gene flow were higher in arid-dwelling individuals, suggesting greater connectivity and interactions between colonies in arid regions in comparison to mesic ones. Using an Aridity Index, we determined that isolation by environment, rather than isolation by geographical distance, best explains the genetic distance between the populations. Further analyses using target loci may determine if there are differing underlying genetic adaptations among populations of C. h. hottentotus. These analyses could help unravel population differences in response to environmental factors within a subspecies of bathyergid mole-rat and determine the adaptive capacity of this small nonmigratory subterranean rodent species in response to aridification in the face of climate change.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies.

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode.

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years.

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

Impact of patrilocality on contrasting patterns of paternal and maternal heritage in Central-West Africa.

Despite their ancient past and high diversity, African populations are the least represented in human population genetic studies. In this study, uniparental markers (mtDNA and Y chromosome) were used to investigate the impact of sociocultural factors on the genetic diversity and inter-ethnolinguistic gene flow in the three major Nigerian groups: Hausa (n = 89), Yoruba (n = 135) and Igbo (n = 134). The results show a distinct history from the maternal and paternal perspectives. The three Nigerian groups present a similar substrate for mtDNA, but not for the Y chromosome. The two Niger-Congo groups, Yoruba and Igbo, are paternally genetically correlated with populations from the same ethnolinguistic affiliation. Meanwhile, the Hausa is paternally closer to other Afro-Asiatic populations and presented a high diversity of lineages from across Africa. When expanding the analyses to other African populations, it is observed that language did not act as a major barrier to female-mediated gene flow and that the differentiation of paternal lineages is better correlated with linguistic than geographic distances. The results obtained demonstrate the impact of patrilocality, a common and well-established practice in populations from Central-West Africa, in the preservation of the patrilineage gene pool and in the affirmation of identity between groups.

Population dynamics of the Lyme disease bacterium, Borrelia burgdorferi, during rapid range expansion in New York State.

Recent changes in climate and human land-use have resulted in alterations of the geographic range of many species, including human pathogens. Geographic range expansion and population growth of human pathogens increase human disease risk. Relatively little empirical work has investigated the impact of range changes on within-population variability, a contributor to both colonization success and adaptive potential, during the precise time in which populations are colonized. This is likely due to the difficulties of collecting appropriate natural samples during the dynamic phase of migration and colonization. We systematically collected blacklegged ticks (Ixodes scapularis) across New York State (NY), USA, between 2006 and 2019, a time period coinciding with a rapid range expansion of ticks and their associated pathogens including Borrelia burgdorferi, the etiological agent of Lyme disease. These samples provide a unique opportunity to investigate the population dynamics of human pathogens as they expand into novel territory. We observed that founder effects were short-lived, as gene flow from long-established populations brought almost all B. burgdorferi lineages to newly colonized populations within just a few years of colonization. By 7 years post-colonization, B. burgdorferi lineage frequency distributions were indistinguishable from long-established sites, indicating that local B. burgdorferi populations experience similar selective pressures despite geographic separation. The B. burgdorferi lineage dynamics elucidate the processes underlying the range expansion and demonstrate that migration into, and selection within, newly colonized sites operate on different time scales.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States.

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range.

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

The biogeographic and evolutionary processes shaping population divergence in Laupala.

Speciation generates biodiversity and the mechanisms involved are thought to vary across the tree of life and across environments. For example, well-studied adaptive radiations are thought to be fuelled by divergent ecological selection, but additionally are influenced heavily by biogeographic, genomic and demographic factors. Mechanisms of non-adaptive radiations, producing ecologically cryptic taxa, have been less well-studied but should likewise be influenced by these latter factors. Comparing among contexts can help pinpoint universal mechanisms and outcomes, especially if we integrate biogeographic, ecological and evolutionary processes. We investigate population divergence in the swordtail cricket Laupala cerasina, a wide-spread endemic on Hawai'i Island and one of 38 ecologically cryptic Laupala species. The nine sampled populations show striking population genetic structure at small spatio-temporal scales. The rapid differentiation among populations and species of Laupala shows that neither a specific geographical context nor ecological opportunity are pre-requisites for rapid divergence. Spatio-temporal patterns in population divergence, population size change, and gene flow are aligned with the chronosequence of the four volcanoes on which L. cerasina occurs and reveal the composite effects of geological dynamics and Quaternary climate change on population dynamics. Spatio-temporal patterns in genetic variation along the genome reveal the interplay of genetic and genomic architecture in shaping population divergence. In early phases of divergence, we find elevated differentiation in genomic regions harbouring mating song loci. In later stages of divergence, we find a signature of linked selection that interacts with recombination rate variation. Comparing our findings with recent work on complementary systems supports the conclusion that mostly universal factors influence the speciation process.

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish.

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.