Assortative mating and persistent reproductive isolation in hybrids

Ongoing introgression of a secondary sexual plumage trait in a stable avian hybrid zone

Hybrid zones are dynamic systems where natural selection, sexual selection, and other evolutionary forces can act on reshuffled combinations of distinct genomes. The movement of hybrid zones, individual traits, or both are of particular interest for understanding the interplay between selective processes. In a hybrid zone involving two lek-breeding birds, secondary sexual plumage traits of Manacus vitellinus, including bright yellow collar and olive belly color, have introgressed ~50 km asymmetrically across the genomic center of the zone into populations more genetically similar to Manacus candei. Males with yellow collars are preferred by females and are more aggressive than parental M. candei, suggesting that sexual selection was responsible for the introgression of male traits. We assessed the spatial and temporal dynamics of this hybrid zone using historical (1989-1994) and contemporary (2017-2020) transect samples to survey both morphological and genetic variation. Genome-wide single nucleotide polymorphism data and several male phenotypic traits show that the genomic center of the zone has remained spatially stable, whereas the olive belly color of male M. vitellinus has continued to introgress over this time period. Our data suggest that sexual selection can continue to shape phenotypes dynamically, independent of a stable genomic transition between species.

Differentiating mechanism from outcome for ancestry-assortative mating in admixed human populations

Population genetic theory, and the empirical methods built upon it, often assume that individuals pair randomly for reproduction. However, natural populations frequently violate this assumption, which may potentially confound genome-wide association studies, selection scans, and demographic inference. Within several recently admixed human populations, empirical genetic studies have reported a correlation in global ancestry proportion between spouses, referred to as ancestry-assortative mating. Here, we use forward genomic simulations to link correlations in ancestry between mates to the underlying mechanistic mate-choice process. We consider the impacts of two types of mate-choice model, using either ancestry-based preferences or social groups as the basis for mate pairing. We find that multiple mate-choice models can produce the same correlations in ancestry proportion between spouses; however, we also highlight alternative analytic approaches and circumstances in which these models may be distinguished. With this work, we seek to highlight potential pitfalls when interpreting correlations in empirical data as evidence for a particular model of human mating practices, as well as to offer suggestions toward development of new best practices for analysis of human ancestry-assortative mating.

Selection leads to remarkable variability in the outcomes of hybridisation across replicate hybrid zones

Hybrid zones have been viewed as an opportunity to see speciation in action. When hybrid zones are replicated, it is assumed that if the same genetic incompatibilities are maintaining reproductive isolation across all instances of secondary contact, those incompatibilities should be identifiable by consistent patterns in the genome. In contrast, changes in allele frequencies due to genetic drift should be idiosyncratic for each hybrid zone. To test this assumption, we simulated 20 replicates of each of 12 hybrid zone scenarios with varied genetic incompatibilities, rates of migration, selection and different starting population size ratios of parental species. We found remarkable variability in the outcomes of hybridisation in replicate hybrid zones, particularly with Bateson-Dobzhansky-Muller incompatibilities and strong selection. We found substantial differences among replicates in the overall genomic composition of individuals, including admixture proportions, inter-specific ancestry complement and number of ancestry junctions. Additionally, we found substantial variation in genomic clines among replicates at focal loci, regardless of locus-specific selection. We conclude that processes other than selection are responsible for some consistent outcomes of hybridisation, whereas selection on incompatibilities can lead to genomically widespread and highly variable outcomes. We highlight the challenge of mapping between pattern and process in hybrid zones and call attention to how selection against incompatibilities will commonly lead to variable outcomes. We hope that this study informs future research on replicate hybrid zones and encourages further development of statistical techniques, theoretical models and exploration of additional axes of variation to understand reproductive isolation.

Pervasive gene flow despite strong and varied reproductive barriers in swordtails

One of the mechanisms that can lead to the formation of new species occurs through the evolution of reproductive barriers. However, recent research has demonstrated that hybridization has been pervasive across the tree of life even in the presence of strong barriers. Swordtail fishes (genus Xiphophorus) are an emerging model system for studying the interface between these barriers and hybridization. We document overlapping mechanisms that act as barriers between closely related species, X. birchmanni and X. cortezi, by combining genomic sequencing from natural hybrid populations, artificial crosses, behavioral assays, sperm performance, and developmental studies. We show that strong assortative mating plays a key role in maintaining subpopulations with distinct ancestry in natural hybrid populations. Lab experiments demonstrate that artificial F1 crosses experience dysfunction: crosses with X. birchmanni females were largely inviable and crosses with X. cortezi females had a heavily skewed sex ratio. Using F2 hybrids we identify several genomic regions that strongly impact hybrid viability. Strikingly, two of these regions underlie genetic incompatibilities in hybrids between X. birchmanni and its sister species X. malinche. Our results demonstrate that ancient hybridization has played a role in the origin of this shared genetic incompatibility. Moreover, ancestry mismatch at these incompatible regions has remarkably similar consequences for phenotypes and hybrid survival in X. cortezi × X. birchmanni hybrids as in X. malinche × X. birchmanni hybrids. Our findings identify varied reproductive barriers that shape genetic exchange between naturally hybridizing species and highlight the complex evolutionary outcomes of hybridization.

Genetic divergences and hybridization within the Sebastes inermis complex

The Sebastes inermis complex includes three sympatric species (Sebastes cheni, viz Sebastes inermis, and Sebastes ventricosus) with clear ecomorphological differences, albeit incomplete reproductive isolation. The presence of putative morphological hybrids (PMH) with plausibly higher fitness than the parent species indicates the need to confirm whether hybridization occurs within the complex. In this sense, we assessed the dynamics of genetic divergence and hybridization within the species complex using a panel of 10 microsatellite loci, and sequences of the mitochondrial control region (D-loop) and the intron-free rhodopsin (RH1) gene. The analyses revealed the presence of three distinct genetic clusters, large genetic distances using D-loop sequences, and distinctive mutations within the RH1 gene. These results are consistent with the descriptions of the three species. Two microsatellite loci had signatures of divergent selection, indicating that they are linked to genomic regions that are crucial for speciation. Furthermore, nonsynonymous mutations within the RH1 gene detected in S. cheni and "Kumano" (a PMH) suggest dissimilar adaptations related to visual perception in dim-light environments. The presence of individuals with admixed ancestry between two species confirmed hybridization. The presence of nonsynonymous mutations within the RH1 gene and the admixed ancestry of the "Kumano" morphotype highlight the potential role of hybridization in generating novelties within the species complex. We discuss possible outcomes of hybridization within the species complex, considering hybrid fitness and assortative mating. Overall, our findings indicate that the genetic divergence of each species is maintained in the presence of hybridization, as expected in a scenario of speciation-with-gene-flow.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

A complex genomic architecture underlies reproductive isolation in a North American oriole hybrid zone

Natural hybrid zones provide powerful opportunities for identifying the mechanisms that facilitate and inhibit speciation. Documenting the extent of genomic admixture allows us to discern the architecture of reproductive isolation through the identification of isolating barriers. This approach is particularly powerful for characterizing the accumulation of isolating barriers in systems exhibiting varying levels of genomic divergence. Here, we use a hybrid zone between two species-the Baltimore (Icterus galbula) and Bullock's (I. bullockii) orioles-to investigate this architecture of reproductive isolation. We combine whole genome re-sequencing with data from an additional 313 individuals amplityped at ancestry-informative markers to characterize fine-scale patterns of admixture, and to quantify links between genes and the plumage traits. On a genome-wide scale, we document several putative barriers to reproduction, including elevated peaks of divergence above a generally high genomic baseline, a large putative inversion on the Z chromosome, and complex interactions between melanogenesis-pathway candidate genes. Concordant and coincident clines for these different genomic regions further suggest the coupling of pre- and post-mating barriers. Our findings of complex and coupled interactions between pre- and post-mating barriers suggest a relatively rapid accumulation of barriers between these species, and they demonstrate the complexities of the speciation process.

Assortative Mating in an Ecological Context: Effects of Mate Choice Errors and Relative Species Abundance on the Frequency and Asymmetry of Hybridization

AbstractThe frequency and asymmetry of mixed-species mating set the initial stage for the ecological and evolutionary implications of hybridization. How such patterns of mixed-species mating, in turn, are influenced by the combination of mate choice errors and relative species abundance remains largely unknown. We develop a mathematical model that generates predictions for how relative species abundances and mate choice errors affect hybridization patterns. When mate choice errors are small (<5%), the highest frequency of hybridization occurs when one of the hybridizing species is at low abundance, but when mate choice errors are high (>5%), the highest hybridization frequency occurs when species occur in equal proportions. Furthermore, females of the less abundant species are overrepresented in mixed-species matings. We compare our theoretical predictions with empirical data on naturally hybridizing Ficedula flycatchers and find that hybridization is highest when the two species occur in equal abundance, implying rather high mate choice errors. We discuss ecological and evolutionary implications of our findings and encourage future work on hybrid zone dynamics that take demographic aspects, such as relative species abundance, into account.

The strength of reproductive isolating barriers in seed plants: Insights from studies quantifying premating and postmating reproductive barriers over the past 15 years

Speciation is driven by the evolution of reproductive isolating barriers that reduce, and ultimately prevent, substantial gene flow between lineages. Despite its central role in evolutionary biology, the process can be difficult to study because it proceeds differently among groups and may occur over long timescales. Due to this complexity, we typically rely on generalizations of empirical data to describe and understand the process. Previous reviews of reproductive isolation (RI) in flowering plants have suggested that prezygotic or extrinsic barriers generally have a stronger effect on reducing gene flow compared to postzygotic or intrinsic barriers. Past conclusions have rested on relatively few empirical estimates of RI; however, RI data have become increasingly abundant over the past 15 years. We analyzed data from recent studies quantifying multiple pre- and postmating barriers in plants and compared the strengths of isolating barriers across 89 taxa pairs using standardized RI metrics. Individual prezygotic barriers were on average stronger than individual postzygotic barriers, and the total strength of prezygotic RI was approximately twice that of postzygotic RI. These findings corroborate that ecological divergence and extrinsic factors, as opposed to solely the accumulation of genetic incompatibilities, are important to speciation and the maintenance of species boundaries in plants. Despite an emphasis in the literature on asymmetric postmating and postzygotic RI, we found that prezygotic barriers acted equally asymmetrically. Overall, substantial variability in the strengths of 12 isolating barriers highlights the great diversity of mechanisms that contribute to plant diversification.

Genetic structure and trait variation within a maple hybrid zone underscore North China as an overlooked diversity hotspot

Tertiary relict flora in East Asia can be divided into northern and southern regions. North China is a diversity hotspot because it can be the secondary contact zone of ancient lineages from the two regions. To test the extent of ancient lineages hybridization and distinguish between the putative species pair Acer pictum subsp. mono and Acer truncatum, we conducted genetic and ecological studies within a maple hybrid zone in North China. Our results suggest that the two lineages of Acer coexist in the hybrid zone and that adult and offspring populations show typical bimodal genetic patterns. Hybrid individuals are established at intermediate altitudes between the two parental lineages. Flowering phenology is divergent between lineages, whereas the complex sexual system of Acer may ensure pollination among lineages. Leaf and fruit morphologies are different between the northern and southern origin lineages, corresponding to A. pictum subsp. mono and A. truncatum, respectively. Reduced gene flow between lineages suggests that they should be considered as two species. However, large morphological variations within each species and the existence of hybrids offer low reliability of species identification based solely on morphological traits. Our study underscores North China as an overlooked diversity hotspot that requires further study in the future.

Assortative mating enhances postzygotic barriers to gene flow via ancestry bundling

Hybridization and subsequent genetic introgression are now known to be common features of the histories of many species, including our own. Following hybridization, selection often purges introgressed DNA genome-wide. While assortative mating can limit hybridization in the first place, it is also known to play an important role in postzygotic selection against hybrids and, thus, the purging of introgressed DNA. However, this role is usually thought of as a direct one: a tendency for mates to be conspecific reduces the sexual fitness of hybrids, reducing the transmission of introgressed ancestry. Here, we explore a second, indirect role of assortative mating as a postzygotic barrier to gene flow. Under assortative mating, parents covary in their ancestry, causing ancestry to be "bundled" in their offspring and later generations. This bundling effect increases ancestry variance in the population, enhancing the efficiency with which postzygotic selection purges introgressed DNA. Using whole-genome simulations, we show that the bundling effect can comprise a substantial portion of mate choice's overall effect as a postzygotic barrier to gene flow. We then derive a simple method for estimating the impact of the bundling effect from standard metrics of assortative mating. Applying this method to data from a diverse set of hybrid zones, we find that the bundling effect increases the purging of introgressed DNA by between 1.2-fold (in a baboon system with weak assortative mating) and 14-fold (in a swordtail system with strong assortative mating). Thus, assortative mating's bundling effect contributes substantially to the genetic isolation of species.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities

Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

Predictability and parallelism in the contemporary evolution of hybrid genomes

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

Assortative mating and epistatic mating-trait architecture induce complex movement of the crow hybrid zone

Hybrid zones provide a window into the evolutionary processes governing species divergence. Yet, the contribution of mate choice to the temporal and spatial stability of hybrid zones remains poorly explored. Here, we investigate the effects of assortative mating on hybrid-zone dynamics by means of a mathematical model parameterized with phenotype and genotype data from the hybrid zone between all-black carrion and gray-coated hooded crows. In the best-fit model, narrow clines of the two mating-trait loci were maintained by a moderate degree of assortative mating inducing pre- and postzygotic isolation via positive frequency-dependent selection. Epistasis between the two loci induced hybrid-zone movement in favor of alleles conveying dark plumage followed by a shift in the opposite direction favoring gray-coated phenotypes ∼ 1 200 generations after secondary contact. Unlinked neutral loci diffused near-unimpeded across the zone. These results were generally robust to the choice of matching rule (self-referencing or parental imprinting) and effects of genetic drift. Overall, this study illustrates under which conditions assortative mating can maintain steep clines in mating-trait loci without generalizing to genome-wide reproductive isolation. It further emphasizes the importance of the genetic mating-trait architecture for spatio-temporal hybrid-zone dynamics.

Two new hybrid populations expand the swordtail hybridization model system

Natural hybridization events provide unique windows into the barriers that keep species apart as well as the consequences of their breakdown. Here, we characterize hybrid populations formed between the northern swordtail fish Xiphophorus cortezi and Xiphophorus birchmanni from collection sites on two rivers. We use simulations and new genetic reference panels to develop sensitive and accurate local ancestry calling in this novel system. Strikingly, we find that hybrid populations on both rivers consist of two genetically distinct subpopulations: a cluster of pure X. birchmanni individuals and one of phenotypically intermediate hybrids that derive ∼85-90% of their genome from X. cortezi. Simulations suggest that initial hybridization occurred ∼150 generations ago at both sites, with little evidence for contemporary gene flow between subpopulations. This population structure is consistent with strong assortative mating between individuals of similar ancestry. The patterns of population structure uncovered here mirror those seen in hybridization between X. birchmanni and its sister species, Xiphophorus malinche, indicating an important role for assortative mating in the evolution of hybrid populations. Future comparisons will provide a window into the shared mechanisms driving the outcomes of hybridization not only among independent hybridization events between the same species but also across distinct species pairs.

A narrow window for geographic cline analysis using genomic data: Effects of age, drift, and migration on error rates

The use of genomic and phenotypic data to scan for outliers is a mainstay for studies of hybridization and speciation. Geographic cline analysis of natural hybrid zones is widely used to identify putative signatures of selection by detecting deviations from baseline patterns of introgression. As with other outlier-based approaches, demographic histories can make neutral regions appear to be under selection and vice versa. In this study, we use a forward-time individual-based simulation approach to evaluate the robustness of geographic cline analysis under different evolutionary scenarios. We modelled multiple stepping-stone hybrid zones with distinct age, deme sizes, and migration rates, and evolving under different types of selection. We found that drift distorts cline shapes and increases false positive rates for signatures of selection. This effect increases with hybrid zone age, particularly if migration between demes is low. Drift can also distort the signature of deleterious effects of hybridization, with genetic incompatibilities and particularly underdominance prone to spurious typing as adaptive introgression. Our results suggest that geographic clines are most useful for outlier analysis in young hybrid zones with large populations of hybrid individuals. Current approaches may overestimate adaptive introgression and underestimate selection against maladaptive genotypes.

Exploring Klebsiella pneumoniae in Healthy Poultry Reveals High Genetic Diversity, Good Biofilm-Forming Abilities and Higher Prevalence in Turkeys Than Broilers

Klebsiella pneumoniae is a well-studied human pathogen for which antimicrobial resistant and hypervirulent clones have emerged globally. K. pneumoniae is also present in a variety of environmental niches, but currently there is a lack of knowledge on the occurrence and characteristics of K. pneumoniae from non-human sources. Certain environmental niches, e.g., animals, may be associated with high K. pneumoniae abundance, and these can constitute a reservoir for further transmission of strains and genetic elements. The aim of this study was to explore and characterize K. pneumoniae from healthy broilers and turkeys. A total of 511 cecal samples (broiler n = 356, turkey n = 155), included in the Norwegian monitoring program for antimicrobial resistance (AMR) in the veterinary sector (NORM-VET) in 2018, were screened for K. pneumoniae by culturing on SCAI agar. K. pneumoniae was detected in 207 (40.5%) samples. Among the broiler samples, 25.8% were positive for K. pneumoniae, in contrast to turkey with 74.2% positive samples (p < 0.01). Antibiotic susceptibility testing was performed, in addition to investigating biofilm production. Whole genome sequencing was performed on 203 K. pneumoniae isolates, and analysis was performed utilizing comparative genomics tools. The genomes grouped into 66 sequence types (STs), with ST35, ST4710 and ST37 being the most prevalent at 13.8%, 7.4%, and 5.4%, respectively. The overall AMR occurrence was low, with only 11.3% of the isolates showing both pheno- and genotypic resistance. Genes encoding aerobactin, salmochelin or yersiniabactin were detected in 47 (23.2%) genomes. Fifteen hypervirulent genomes belonging to ST4710 and isolated from turkey were identified. These all encoded the siderophore virulence loci iuc5 and iro5 on an IncF plasmid. Isolates from both poultry species displayed good biofilm-forming abilities with an average of OD₅₉₅ 0.69 and 0.64. To conclude, the occurrence of K. pneumoniae in turkey was significantly higher than in broiler, indicating that turkey might be an important zoonotic reservoir for K. pneumoniae compared to broilers. Furthermore, our results show a highly diverse K. pneumoniae population in poultry, low levels of antimicrobial resistance, good biofilm-forming abilities and a novel hypervirulent ST4710 clone circulating in the turkey population.

Mixed Mating in a Multi-Origin Population Suggests High Potential for Genetic Rescue in North Island Brown Kiwi, Apteryx mantelli

Reinforcement translocations are increasingly utilised in conservation with the goal of achieving genetic rescue. However, concerns regarding undesirable results, such as genetic homogenisation or replacement, are widespread. One factor influencing translocation outcomes is the rate at which the resident and the introduced individuals interbreed. Consequently, post-release mate choice is a key behaviour to consider in conservation planning. Here we studied mating, and its consequences for genomic admixture, in the North Island brown kiwi Apteryx mantelli population on Ponui Island which was founded by two translocation events over 50 years ago. The two source populations used are now recognised as belonging to two separate management units between which birds differ in size and are genetically differentiated. We examined the correlation between male and female morphometrics for 17 known pairs and quantified the relatedness of 20 pairs from this admixed population. In addition, we compared the genetic similarity and makeup of 106 Ponui Island birds, including 23 known pairs, to birds representing the source populations for the original translocations. We found no evidence for size-assortative mating. On the contrary, genomic SNP data suggested that kiwi of one feather did not flock together, meaning that mate choice resulted in pairing between individuals that were less related than expected by random chance. Furthermore, the birds in the current Ponui Island population were found to fall along a gradient of genomic composition consistent with non-clustered representation of the two parental genomes. These findings indicate potential for successful genetic rescue in future Apteryx reinforcement translocations, a potential that is currently under utilised due to restrictive translocation policies. In light of our findings, we suggest that reconsideration of these policies could render great benefits for the future diversity of this iconic genus in New Zealand.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

Ongoing production of low-fitness hybrids limits range overlap between divergent cryptic species

Contact zones between recently diverged taxa provide opportunities to examine the causes of reproductive isolation and the processes that determine whether two species can coexist over a broad region. The Pacific wren (Troglodytes pacificus) and winter wren (Troglodytes hiemalis) are two morphologically similar songbirds that started diverging about 4 million years ago, older than most sister species pairs of temperate songbirds. The ranges of these species come into narrow contact in western Canada, where the two species remain distinct. To assess evidence for differentiation, hybridization and introgression in this system, we examined variation in over 250,000 single nucleotide polymorphism markers distributed across the genome. The two species formed highly divergent genetic clusters, consistent with long-term differentiation. In a set of 75 individuals, two first-generation hybrids (i.e., F₁ 's) were detected, indicating only moderate levels of assortative mating between these taxa. We found no recent backcrosses or other evidence of recent breeding success of F₁ 's, indicating very low or zero fitness of F₁  hybrids. Examination of genomic variation shows evidence for only a single backcrossing event many generations ago. The moderate rate of hybridization combined with very low F₁  hybrid fitness is expected to result in a population sink in the contact zone, largely explaining the narrow overlap of the two species. If such dynamics are common in nature, they could explain the narrow range overlap often observed between pairs of closely related species.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

The Genetic Architecture of Variation in the Sexually Selected Sword Ornament and Its Evolution in Hybrid Populations

Biologists since Darwin have been fascinated by the evolution of sexually selected ornaments, particularly those that reduce viability. Uncovering the genetic architecture of these traits is key to understanding how they evolve and are maintained. Here, we investigate the genetic architecture and evolutionary loss of a sexually selected ornament, the "sword" fin extension that characterizes many species of swordtail fish (Xiphophorus). Using sworded and swordless sister species of Xiphophorus, we generated a mapping population and show that the sword ornament is polygenic-with ancestry across the genome explaining substantial variation in the trait. After accounting for the impacts of genome-wide ancestry, we identify one major-effect quantitative trait locus (QTL) that explains ~5% of the overall variation in the trait. Using a series of approaches, we narrow this large QTL interval to several likely candidate genes, including genes involved in fin regeneration and growth. Furthermore, we find evidence of selection on ancestry at one of these candidates in four natural hybrid populations, consistent with selection against the sword in these populations.

Evolutionary novelty in communication between the sexes

The diversity of signalling traits within and across taxa is vast and striking, prompting us to consider how novelty evolves in the context of animal communication. Sexual selection contributes to diversification, and here we endeavour to understand the initial conditions that facilitate the maintenance or elimination of new sexual signals and receiver features. New sender and receiver variants can occur through mutation, plasticity, hybridization and cultural innovation, and the initial conditions of the sender, the receiver and the environment then dictate whether a novel cue becomes a signal. New features may arise in the sender, the receiver or both simultaneously. We contend that it may be easier than assumed to evolve new sexual signals because sexual signals may be arbitrary, sexual conflict is common and receivers are capable of perceiving much more of the world than just existing sexual signals. Additionally, changes in the signalling environment can approximate both signal and receiver changes through a change in transmission characteristics of a given environment or the use of new environments. The Anthropocene has led to wide-scale disruption of the environment and may thus generate opportunity to directly observe the evolution of new signals to address questions that are beyond the reach of phylogenetic approaches.

Contact-Chemosensory Evolution Underlying Reproductive Isolation in Drosophila Species

The main theme of the review is how changes in pheromone biochemistry and the sensory circuits underlying pheromone detection contribute to mate choice and reproductive isolation. The review focuses primarily on gustatory and non-volatile signals in Drosophila. Premating isolation is prevalent among closely related species. In Drosophila, preference for conspecifics against other species in mate choice underlies premating isolation, and such preference relies on contact chemosensory communications between a female and male along with other biological factors. For example, although D. simulans and D. melanogaster are sibling species that yield hybrids, their premating isolation is maintained primarily by the contrasting effects of 7,11-heptacosadiene (7,11-HD), a predominant female pheromone in D. melanogaster, on males of the two species: it attracts D. melanogaster males and repels D. simulans males. The contrasting preference for 7,11-HD in males of these two species is mainly ascribed to opposite effects of 7,11-HD on neural activities in the courtship decision-making neurons in the male brain: 7,11-HD provokes both excitatory and inhibitory inputs in these neurons and differences in the balance between the two counteracting inputs result in the contrasting preference for 7,11-HD, i.e., attraction in D. melanogaster and repulsion in D. simulans. Introduction of two double bonds is a key step in 7,11-HD biosynthesis and is mediated by the desaturase desatF, which is active in D. melanogaster females but transcriptionally inactivated in D. simulans females. Thus, 7,11-HD biosynthesis diversified in females and 7,11-HD perception diversified in males, yet it remains elusive how concordance of the changes in the two sexes was attained in evolution.

The importance of intrinsic postzygotic barriers throughout the speciation process

Intrinsic postzygotic barriers can play an important and multifaceted role in speciation, but their contribution is often thought to be reserved to the final stages of the speciation process. Here, we review how intrinsic postzygotic barriers can contribute to speciation, and how this role may change through time. We outline three major contributions of intrinsic postzygotic barriers to speciation. (i) reduction of gene flow: intrinsic postzygotic barriers can effectively reduce gene exchange between sympatric species pairs. We discuss the factors that influence how effective incompatibilities are in limiting gene flow. (ii) early onset of species boundaries via rapid evolution: intrinsic postzygotic barriers can evolve between recently diverged populations or incipient species, thereby influencing speciation relatively early in the process. We discuss why the early origination of incompatibilities is expected under some biological models, and detail how other (and often less obvious) incompatibilities may also serve as important barriers early on in speciation. (iii) reinforcement: intrinsic postzygotic barriers can promote the evolution of subsequent reproductive isolation through processes such as reinforcement, even between relatively recently diverged species pairs. We incorporate classic and recent empirical and theoretical work to explore these three facets of intrinsic postzygotic barriers, and provide our thoughts on recent challenges and areas in the field in which progress can be made. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Leveraging eDNA to expand the study of hybrid zones

Hybrid zones are important windows into ecological and evolutionary processes. Our understanding of the significance and prevalence of hybridization in nature has expanded with the generation and analysis of genome‐spanning data sets. That said, most hybridization research still has restricted temporal and spatial resolution, which limits our ability to draw broad conclusions about evolutionary and conservation related outcomes. Here, we argue that rapidly advancing environmental DNA (eDNA) methodology could be adopted for studies of hybrid zones to increase temporal sampling (contemporary and historical), refine and geographically expand sampling density, and collect data for taxa that are difficult to directly sample. Genomic data in the environment offer the potential for near real‐time biological tracking of hybrid zones, and eDNA provides broad, but as yet untapped, potential to address eco‐evolutionary questions.

Natural hybridization reveals incompatible alleles that cause melanoma in swordtail fish

The establishment of reproductive barriers between populations can fuel the evolution of new species. A genetic framework for this process posits that "incompatible" interactions between genes can evolve that result in reduced survival or reproduction in hybrids. However, progress has been slow in identifying individual genes that underlie hybrid incompatibilities. We used a combination of approaches to map the genes that drive the development of an incompatibility that causes melanoma in swordtail fish hybrids. One of the genes involved in this incompatibility also causes melanoma in hybrids between distantly related species. Moreover, this melanoma reduces survival in the wild, likely because of progressive degradation of the fin. This work identifies genes underlying a vertebrate hybrid incompatibility and provides a glimpse into the action of these genes in natural hybrid populations.

Assortative mating by population of origin in a mechanistic model of admixture

Populations whose mating pairs have levels of similarity in phenotypes or genotypes that differ systematically from the level expected under random mating are described as experiencing assortative mating. Excess similarity in mating pairs is termed positive assortative mating, and excess dissimilarity is negative assortative mating. In humans, empirical studies suggest that mating pairs from various admixed populations - whose ancestry derives from two or more source populations - possess correlated ancestry components that indicate the occurrence of positive assortative mating on the basis of ancestry. Generalizing a two-sex mechanistic admixture model, we devise a model of one form of ancestry-assortative mating that occurs through preferential mating based on source population. Under the model, we study the moments of the admixture fraction distribution for different assumptions about mating preferences, including both positive and negative assortative mating by population. We demonstrate that whereas the mean admixture under assortative mating is equivalent to that of a corresponding randomly mating population, the variance of admixture depends on the level and direction of assortative mating. We consider two special cases of assortative mating by population: first, a single admixture event, and second, constant contributions to the admixed population over time. In contrast to standard settings in which positive assortment increases variation within a population, certain assortative mating scenarios allow the variance of admixture to decrease relative to a corresponding randomly mating population: with the three populations we consider, the variance-increasing effect of positive assortative mating within a population might be overwhelmed by a variance-decreasing effect emerging from mating preferences involving other pairs of populations. The effect of assortative mating is smaller on the X chromosome than on the autosomes because inheritance of the X in males depends only on the mother's ancestry, not on the mating pair. Because the variance of admixture is informative about the timing of admixture and possibly about sex-biased admixture contributions, the effects of assortative mating are important to consider in inferring features of population history from distributions of admixture values. Our model provides a framework to quantitatively study assortative mating under flexible scenarios of admixture over time.

Genomic regions underlying metabolic and neuronal signaling pathways are temporally consistent in a moving avian hybrid zone

The study of hybrid zones can provide insight into the genetic basis of species differences that are relevant for the maintenance of reproductive isolation. Hybrid zones can also provide insight into climate change, species distributions, and evolution. The hybrid zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) is shifting northward in response to increasing winter temperatures but is not increasing in width. This pattern indicates strong selection against chickadees with admixed genomes. Using high-resolution genomic data, we identified regions of the genomes that are outliers in both time points and do not introgress between the species; these regions may be involved in the maintenance of reproductive isolation. Genes involved in metabolic regulation processes were overrepresented in this dataset. Several gene ontology categories were also temporally consistent-including glutamate signaling, synaptic transmission, and catabolic processes-but the nucleotide variants leading to this pattern were not. Our results support recent findings that hybrids between black-capped and Carolina chickadees have higher basal metabolic rates than either parental species and suffer spatial memory and problem-solving deficits. Metabolic breakdown, as well as spatial memory and problem-solving, in hybrid chickadees may act as strong postzygotic isolation mechanisms in this moving hybrid zone.

Gene transfer across species boundaries in bryophytes: evidence from major life cycle stages in Homalothecium lutescens and H. sericeum

BACKGROUND AND AIMS

The mosses Homalothecium lutescens and H. sericeum are genetically, morphologically and ecologically differentiated; mixed populations sometimes occur. In sympatric populations, intermediate character states among gametophytes and sporophytes have been observed, suggesting hybridization and introgression in such populations.

METHODS

We determined genotypes using bi-allelic co-dominant single nucleotide polymorphism (SNP) markers, specific to either H. lutescens or H. sericeum, to estimate the degree of genetic mixing in 449 moss samples collected from seven sympatric and five allopatric populations on the island of Öland, south Sweden. The samples represented three generations: haploid maternal gametophytes; diploid sporophytes; and haploid sporelings.

KEY RESULTS

Admixture analyses of SNP genotypes identified a majority as pure H. lutescens or H. sericeum, but 76 samples were identified as mildly admixed (17 %) and 17 samples (3.8 %) as strongly admixed. Admixed samples were represented in all three generations in several populations. Hybridization and introgression were bidirectional.

CONCLUSIONS

Our results demonstrate that admixed genomes are transferred between the generations, so that the populations behave as true hybrid zones. Earlier studies of sympatric bryophyte populations with admixed individuals have not been able to show that admixed alleles are transferred beyond the first generation. The presence of true hybrid zones has strong evolutionary implications because genetic material transferred across species boundaries can be directly exposed to selection in the long-lived haploid generation of the bryophyte life cycle, and contribute to local adaptation, long-term survival and speciation.

Divergent neurogenomic responses shape social learning of both personality and mate preference

Behavior plays a fundamental role in shaping the origin and fate of species. Mating decisions can act to promote or restrict gene flow, as can personality traits that influence dispersal and niche use. Mate choice and personality are often both learned and therefore influenced by an individual's social environment throughout development. Likewise, the molecular pathways that shape these behaviors may also be co-expressed. In this study on swordtail fish (Xiphophorus birchmanni), we show that female mating preferences for species-typical pheromone cues are entirely dependent on social experience with adult males. Experience with adults also shapes development along the shy-bold personality axis, with shy behaviors arising from exposure to risk-averse heterospecifics as a potential stress-coping strategy. In maturing females, conspecific exposure results in a strong upregulation of olfaction and vision genes compared with heterospecific exposure, as well as immune response genes previously linked to anxiety, learning and memory. Conversely, heterospecific exposure involves an increased expression of genes important for neurogenesis, synaptic plasticity and social decision-making. We identify subsets of genes within the social decision-making network and with known stress-coping roles that may be directly coupled to the olfactory processes females rely on for social communication. Based on these results, we conclude that the social environment affects the neurogenomic trajectory through which socially sensitive behaviors are learned, resulting in adult phenotypes adapted for specific social groupings.

Establishment, Hybridization, Dispersal, Impact, and Decline of Diorhabda spp. (Coleoptera: Chrysomelidae) Released for Biological Control of Tamarisk in Texas and New Mexico

Three Diorhabda spp. tamarisk beetles (Coleoptera: Chrysomelidae) were established in Texas from 2003 to 2010 for biological control of tamarisk (Tamarix spp.): Mediterranean tamarisk beetles, D. elongata (Brullé) from Greece, also established in New Mexico; subtropical tamarisk beetles, D. sublineata (Lucas) from Tunisia; and larger tamarisk beetles, D. carinata (Faldermann) from Uzbekistan. More than one million tamarisk beetles were released at 99 sites. Species establishment success ranged from 52 to 83%. All three species now co-occur in New Mexico with the northern tamarisk beetles, D. carinulata (Desbrochers). A phenotypic hybrid scoring system was developed to assess Diorhabda phenotype distributions and character mixing in hybrid zones. Widespread field populations of bispecific hybrid phenotypes for D. carinata/D. elongata and D. sublineata/D. elongata rapidly appeared following contact of parental species. Initial distributions and dispersal of Diorhabda spp. and hybrids are mapped for Texas, New Mexico, Oklahoma, and Kansas, where they produced large-scale tamarisk defoliation and localized dieback for 3-4 yr. However, populations subsequently severely declined, now producing only isolated defoliation and allowing tamarisk to recover. Diorhabda sublineata and D. elongata temporarily produced nontarget spillover defoliation of ornamental athel, Tamarix aphylla (L.) Karst, along the Rio Grande. Hybrid phenotypes were generally bimodally distributed, indicating some degree of reproductive isolation. Additional diagnostic phenotypic characters in males allowed more precise hybrid scoring. Character mixing in some hybrid populations approached or reached that of a hybrid swarm. The significance of hybridization for tamarisk biocontrol is discussed.

Eukaryote hybrid genomes

Interspecific hybridization is the process where closely related species mate and produce offspring with admixed genomes. The genomic revolution has shown that hybridization is common, and that it may represent an important source of novel variation. Although most interspecific hybrids are sterile or less fit than their parents, some may survive and reproduce, enabling the transfer of adaptive variants across the species boundary, and even result in the formation of novel evolutionary lineages. There are two main variants of hybrid species genomes: allopolyploid, which have one full chromosome set from each parent species, and homoploid, which are a mosaic of the parent species genomes with no increase in chromosome number. The establishment of hybrid species requires the development of reproductive isolation against parental species. Allopolyploid species often have strong intrinsic reproductive barriers due to differences in chromosome number, and homoploid hybrids can become reproductively isolated from the parent species through assortment of genetic incompatibilities. However, both types of hybrids can become further reproductively isolated, gaining extrinsic isolation barriers, by exploiting novel ecological niches, relative to their parents. Hybrids represent the merging of divergent genomes and thus face problems arising from incompatible combinations of genes. Thus hybrid genomes are highly dynamic and undergo rapid evolutionary change, including genome stabilization in which selection against incompatible combinations results in fixation of compatible ancestry block combinations within the hybrid species. The potential for rapid adaptation or speciation makes hybrid genomes a particularly exciting subject of in evolutionary biology. Here we summarize how introgressed alleles or hybrid species can establish and how the resulting hybrid genomes evolve.

micRocounter: Microsatellite Characterization in Genome Assemblies

Microsatellites are repetitive DNA sequences usually found in non-coding regions of the genome. Their quantification and analysis have applications in fields from population genetics to evolutionary biology. As genome assemblies become commonplace, the need for software that can facilitate analyses has never been greater. In particular, R packages that can analyze genomic data are particularly important since this is one of the most popular software environments for biologists. We created an R package, micRocounter, to quantify microsatellites. We have optimized our package for speed, accessibility, and portability, making the automated analysis of large genomic data sets feasible. Computationally intensive algorithms were built in C++ to increase speed. Tests using benchmark datasets show a 200-fold improvement in speed over existing software. A moderately sized genome of 500 Mb can be processed in under 50 sec. Results are output as an object in R increasing accessibility and flexibility for practitioners.

Caspase-8-dependent control of NK- and T cell responses during cytomegalovirus infection

Caspase-8 (CASP8) impacts antiviral immunity in expected as well as unexpected ways. Mice with combined deficiency in CASP8 and RIPK3 cannot support extrinsic apoptosis or RIPK3-dependent programmed necrosis, enabling studies of CASP8 function without complications of unleashed necroptosis. These extrinsic cell death pathways are naturally targeted by murine cytomegalovirus (MCMV)-encoded cell death suppressors, showing they are key to cell-autonomous host defense. Remarkably, Casp8^-/-Ripk3^-/-, Ripk1^-/-Casp8^-/-Ripk3^-/- and Casp8^-/-Ripk3^K51A/K51A mice mount robust antiviral T cell responses to control MCMV infection. Studies in Casp8^-/-Ripk3^-/- mice show that CASP8 restrains expansion of MCMV-specific natural killer (NK) and CD8 T cells without compromising contraction or immune memory. Infected Casp8^-/-Ripk3^-/- or Casp8^-/-Ripk3^K51A/K51A mice have higher levels of virus-specific NK cells and CD8 T cells compared to matched RIPK3-deficient littermates or WT mice. CASP8, likely acting downstream of Fas death receptor, dampens proliferation of CD8 T cells during expansion. Importantly, contraction proceeds unimpaired in the absence of extrinsic death pathways owing to intact Bim-dependent (intrinsic) apoptosis. CD8 T cell memory develops in Casp8^-/-Ripk3^-/- mice, but memory inflation characteristic of MCMV infection is not sustained in the absence of CASP8 function. Despite this, Casp8^-/-Ripk3^-/- mice are immune to secondary challenge. Interferon (IFN)γ is recognized as a key cytokine for adaptive immune control of MCMV. Ifngr^-/-Casp8^-/-Ripk3^-/- mice exhibit increased lifelong persistence in salivary glands as well as lungs compared to Ifngr^-/- and Casp8^-/-Ripk3^-/- mice. Thus, mice deficient in CASP8 and RIPK3 are more dependent on IFNγ mechanisms for sustained T cell immune control of MCMV. Overall, appropriate NK- and T cell immunity to MCMV is dependent on host CASP8 function independent of RIPK3-regulated pathways.

Mate Choice versus Mate Preference: Inferences about Color-Assortative Mating Differ between Field and Lab Assays of Poison Frog Behavior

Codivergence of mating traits and mate preferences can lead to behavioral isolation among lineages in early stages of speciation. However, mate preferences limit gene flow only when expressed as mate choice, and numerous factors might be more important than preferences in nature. In the extremely color polytypic strawberry poison frog (Oophaga pumilio), female mate preferences have codiverged with color in most allopatric populations tested. Whether these lab-assayed preferences predict mating (gene flow) in the wild remains unclear. We observed courting pairs in a natural contact zone between red and blue lineages until oviposition or courtship termination. We found color-assortative mating in a disturbed habitat with high population density but not in a secondary forest with lower density. Our results suggest color-assortative O. pumilio mate choice in the wild but also mating patterns that do not match those predicted by lab-assayed preferences.

Evolution favours aging in populations with assortative mating and in sexually dimorphic populations

Since aging seems omnipresent, many authors regard it as an inevitable consequence of the laws of physics. However, recent research has conclusively shown that some organisms do not age, or at least do not age on a scale comparable with other aging organisms. This begets the question why aging evolved in some organisms yet not in others. Here we present a simulation model of competition between aging and non-aging individuals in a sexually reproducing population. We find that the aging individuals may outcompete the non-aging ones if they have a sufficiently but not excessively higher initial fecundity or if individuals mate assortatively with respect to their own phenotype. Furthermore, the aging phenotype outcompetes the non-aging one or resists dominance of the latter for a longer period in populations composed of genuine males and females compared to populations of simultaneous hermaphrodites. Finally, whereas sterilizing parasites promote non-aging, the effect of mortality-enhancing parasites is to enable longer persistence of the aging phenotype relative to when parasites are absent. Since the aging individuals replace the non-aging ones in diverse scenarios commonly found in nature, our study provides important insights into why aging has evolved in most, but not all organisms.

Role of sexual imprinting in assortative mating and premating isolation in Darwin's finches

Global biodiversity is being degraded at an unprecedented rate, so it is important to preserve the potential for future speciation. Providing for the future requires understanding speciation as a contemporary ecological process. Phylogenetically young adaptive radiations are a good choice for detailed study because diversification is ongoing. A key question is how incipient species become reproductively isolated from each other. Barriers to gene exchange have been investigated experimentally in the laboratory and in the field, but little information exists from the quantitative study of mating patterns in nature. Although the degree to which genetic variation underlying mate-preference learning is unknown, we provide evidence that two species of Darwin's finches imprint on morphological cues of their parents and mate assortatively. Statistical evidence of presumed imprinting is stronger for sons than for daughters and is stronger for imprinting on fathers than on mothers. In combination, morphology and species-specific song learned from the father constitute a barrier to interbreeding. The barrier becomes stronger the more the species diverge morphologically and ecologically. It occasionally breaks down, and the species hybridize. Hybridization is most likely to happen when species are similar to each other in adaptive morphological traits, e.g., body size and beak size and shape. Hybridization can lead to the formation of a new species reproductively isolated from the parental species as a result of sexual imprinting. Conservation of sufficiently diverse natural habitat is needed to sustain a large sample of extant biota and preserve the potential for future speciation.

Digest: The evolution of sexual imprinting as an assortative mating mechanism

In this issue, Yeh et al. (2018) investigated whether sexual imprinting could act as an assortative mating mechanism, reducing hybridization and increasing premating isolation. While they indeed find that imprinting leads to assortative mating and reduced hybridization, the strength at which imprinting evolves is usually intermediate, because it is counterbalanced by the costs of imprinting and the benefits of adaptive hybridization. Thus, while sexual imprinting can act as an assortative mating mechanism, it is often not the sole element of female mate choice.

Natural selection interacts with recombination to shape the evolution of hybrid genomes

To investigate the consequences of hybridization between species, we studied three replicate hybrid populations that formed naturally between two swordtail fish species, estimating their fine-scale genetic map and inferring ancestry along the genomes of 690 individuals. In all three populations, ancestry from the "minor" parental species is more common in regions of high recombination and where there is linkage to fewer putative targets of selection. The same patterns are apparent in a reanalysis of human and archaic admixture. These results support models in which ancestry from the minor parental species is more likely to persist when rapidly uncoupled from alleles that are deleterious in hybrids. Our analyses further indicate that selection on swordtail hybrids stems predominantly from deleterious combinations of epistatically interacting alleles.

Female Choice Undermines the Emergence of Strong Sexual Isolation between Locally Adapted Populations of Atlantic Mollies (Poecilia mexicana)

Divergent selection between ecologically dissimilar habitats promotes local adaptation, which can lead to reproductive isolation (RI). Populations in the Poecilia mexicana species complex have independently adapted to toxic hydrogen sulfide and show varying degrees of RI. Here, we examined the variation in the mate choice component of prezygotic RI. Mate choice tests across drainages (with stimulus males from another drainage) suggest that specific features of the males coupled with a general female preference for yellow color patterns explain the observed variation. Analyses of male body coloration identified the intensity of yellow fin coloration as a strong candidate to explain this pattern, and common-garden rearing suggested heritable population differences. Male sexual ornamentation apparently evolved differently across sulfide-adapted populations, for example because of differences in natural counterselection via predation. The ubiquitous preference for yellow color ornaments in poeciliid females likely undermines the emergence of strong RI, as female discrimination in favor of own males becomes weaker when yellow fin coloration in the respective sulfide ecotype increases. Our study illustrates the complexity of the (partly non-parallel) pathways to divergence among replicated ecological gradients. We suggest that future work should identify the genomic loci involved in the pattern reported here, making use of the increasing genomic and transcriptomic datasets available for our study system.