Alan Brelsford

Hybridization and speciation

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near‐instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky–Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock‐on effects on speciation both within and outside regions of hybridization.

The biogeography of mitochondrial and nuclear discordance in animals

Combining nuclear (nuDNA) and mitochondrial DNA (mtDNA) markers has improved the power of molecular data to test phylogenetic and phylogeographic hypotheses and has highlighted the limitations of studies using only mtDNA markers. In fact, in the past decade, many conflicting geographic patterns between mitochondrial and nuclear genetic markers have been identified (i.e. mito‐nuclear discordance). Our goals in this synthesis are to: (i) review known cases of mito‐nuclear discordance in animal systems, (ii) to summarize the biogeographic patterns in each instance and (iii) to identify common drivers of discordance in various groups. In total, we identified 126 cases in animal systems with strong evidence of discordance between the biogeographic patterns obtained from mitochondrial DNA and those observed in the nuclear genome. In most cases, these patterns are attributed to adaptive introgression of mtDNA, demographic disparities and sex‐biased asymmetries, with some studies also implicating hybrid zone movement, human introductions and Wolbachia infection in insects. We also discuss situations where divergent mtDNA clades seem to have arisen in the absence of geographic isolation. For those cases where foreign mtDNA haplotypes are found deep within the range of a second taxon, data suggest that those mtDNA haplotypes are more likely to be at a high frequency and are commonly driven by sex‐biased asymmetries and/or adaptive introgression. In addition, we discuss the problems with inferring the processes causing discordance from biogeographic patterns that are common in many studies. In many cases, authors presented more than one explanation for discordant patterns in a given system, which indicates that likely more data are required. Ideally, to resolve this issue, we see important future work shifting focus from documenting the prevalence of mito‐nuclear discordance towards testing hypotheses regarding the drivers of discordance. Indeed, there is great potential for certain cases of mitochondrial introgression to become important natural systems within which to test the effect of different mitochondrial genotypes on whole‐animal phenotypes.

The biogeography of mitochondrial and nuclear discordance in animals.

Combining nuclear (nuDNA) and mitochondrial DNA (mtDNA) markers has improved the power of molecular data to test phylogenetic and phylogeographic hypotheses and has highlighted the limitations of studies using only mtDNA markers. In fact, in the past decade, many conflicting geographic patterns between mitochondrial and nuclear genetic markers have been identified (i.e. mito‐nuclear discordance). Our goals in this synthesis are to: (i) review known cases of mito‐nuclear discordance in animal systems, (ii) to summarize the biogeographic patterns in each instance and (iii) to identify common drivers of discordance in various groups. In total, we identified 126 cases in animal systems with strong evidence of discordance between the biogeographic patterns obtained from mitochondrial DNA and those observed in the nuclear genome. In most cases, these patterns are attributed to adaptive introgression of mtDNA, demographic disparities and sex‐biased asymmetries, with some studies also implicating hybrid zone movement, human introductions and Wolbachia infection in insects. We also discuss situations where divergent mtDNA clades seem to have arisen in the absence of geographic isolation. For those cases where foreign mtDNA haplotypes are found deep within the range of a second taxon, data suggest that those mtDNA haplotypes are more likely to be at a high frequency and are commonly driven by sex‐biased asymmetries and/or adaptive introgression. In addition, we discuss the problems with inferring the processes causing discordance from biogeographic patterns that are common in many studies. In many cases, authors presented more than one explanation for discordant patterns in a given system, which indicates that likely more data are required. Ideally, to resolve this issue, we see important future work shifting focus from documenting the prevalence of mito‐nuclear discordance towards testing hypotheses regarding the drivers of discordance. Indeed, there is great potential for certain cases of mitochondrial introgression to become important natural systems within which to test the effect of different mitochondrial genotypes on whole‐animal phenotypes.

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

INCIPIENT SPECIATION DESPITE LITTLE ASSORTATIVE MATING: THE YELLOW-RUMPED WARBLER HYBRID ZONE

Hybrid zones between recently diverged taxa are natural laboratories for speciation research, allowing us to determine whether there is reproductive isolation between divergent forms and the causes of that isolation. We present a study of a classic avian hybrid zone in North America between two subspecies of the yellow-rumped warbler (Dendroica coronata). Although previous work has shown very little differentiation in mitochondrial DNA across this hybrid zone, we identified two nuclear loci (one sexlinked and one autosomal) that show fixed differences across the hybrid zone, in a close concordance with patterns of plumage variation. Temporal stability and limited width of the hybrid zone, along with substantial linkage disequilibrium between these two diagnostic markers in the center of the zone, indicate that there is moderate reproductive isolation between these populations, with an estimated strength of selection maintaining the zone of 18%. Pairing data indicate that assortative mating is either very weak or absent, suggesting that this reproductive isolation is largely due to postmating barriers. Thus, despite extensive hybridization the two forms are distinct evolutionary groups carrying genes for divergent adaptive peaks, and this situation appears relatively stable.

Hybrid origin of Audubon’s warbler

Several animal species have recently been shown to have hybrid origins, but no avian examples have been documented with molecular evidence. We investigate whether the Audubon’s warbler (Dendroica auduboni), one of four visually distinct species in the yellow‐rumped warbler complex, has originated through hybridization between two other species in this group, the myrtle warbler (D. coronata) and black‐fronted warbler (D. nigrifrons). Analysis of nuclear amplified fragment length polymorphism (AFLP) and sequence markers shows that Audubon’s warblers are genetically intermediate and carry a mixture of alleles otherwise found only in one or the other of their putative parental species. Audubon’s warblers also carry two deeply divergent mitochondrial DNA lineages, each shared with only one putative parental form. Broad clines between Audubon’s and black‐fronted warblers in AFLP markers call into question the validity of these two forms as full species; nevertheless, our results suggest that the Audubon’s warbler probably originated through hybridization between two long‐diverged species. It is likely that more cases of avian species of hybrid origin will be revealed by surveys of variation in nuclear DNA and other traits.

Plumage Genes and Little Else Distinguish the Genomes of Hybridizing Warblers

When related taxa hybridize extensively, their genomes may become increasingly homogenized over time. This mixing via hybridization creates conservation challenges when it reduces genetic or phenotypic diversity and when it endangers previously distinct species via genetic swamping [1]. However, hybridization also facilitates admixture mapping of traits that distinguish each species and the associated genes that maintain distinctiveness despite ongoing gene flow [2]. We address these dual aspects of hybridization in the golden-winged/blue-winged warbler complex, two phenotypically divergent warblers that are indistinguishable using traditional molecular markers and that draw substantial conservation attention [3-5]. Whole-genome comparisons show that differentiation is extremely low: only six small genomic regions exhibit strong differences. Four of these divergence peaks occur in proximity to genes known to be involved in feather development or pigmentation: agouti signaling protein (ASIP), follistatin (FST), ecodysplasin (EDA), wingless-related integration site (Wnt), and beta-carotene oxygenase 2 (BCO2). Throat coloration-the most striking plumage difference between these warblers-is perfectly associated with the promoter region of agouti, and genotypes at this locus obey simple Mendelian recessive inheritance of the black-throated phenotype characteristic of golden-winged warblers. The more general pattern of genomic similarity between these warblers likely results from a protracted period of hybridization, contradicting the broadly accepted hypothesis that admixture results from solely anthropogenic habitat change in the past two centuries [4]. Considered in concert, these results are relevant to both the genetic architecture of avian feather pigmentation and the evolutionary history and conservation challenges associated with these declining songbirds.

HOMOLOGOUS SEX CHROMOSOMES IN THREE DEEPLY DIVERGENT ANURAN SPECIES

Comparative genomic studies are revealing that, in sharp contrast with the strong stability found in birds and mammals, sex determination mechanisms are surprisingly labile in cold‐blooded vertebrates, with frequent transitions between different pairs of sex chromosomes. It was recently suggested that, in context of this high turnover, some chromosome pairs might be more likely than others to be co‐opted as sex chromosomes. Empirical support, however, is still very limited. Here we show that sex‐linked markers from three highly divergent groups of anurans map to Xenopus tropicalis scaffold 1, a large part of which is homologous to the avian sex chromosome. Accordingly, the bird sex determination gene DMRT1, known to play a key role in sex differentiation across many animal lineages, is sex linked in all three groups. Our data provide strong support for the idea that some chromosome pairs are more likely than others to be co‐opted as sex chromosomes because they harbor key genes from the sex determination pathway.

Bioinformatic processing of RAD‐seq data dramatically impacts downstream population genetic inference

1. Restriction site-associated DNA sequencing (RAD-seq) provides high-resolution population genomic data at low cost, and has become an important component in ecological and evolutionary studies. As with all high-throughput technologies, analytic strategies require critical validation to ensure precise and unbiased interpretation. 2. To test the impact of bioinformatic data processing on downstream population genetic inferences, we analysed mammalian RAD-seq data (>100 individuals) with 312 combinations of methodology (de novo vs. mapping to references of increasing divergence) and filtering criteria (missing data, HWE, F-IS, coverage, mapping and genotype quality). In an effort to identify commonalities and biases in all pipelines, we computed summary statistics (nr. loci, nr. SNP, pi, Het(obs), F-IS, F-ST, N-e and m) and compared the results to independent null expectations (isolation-by-distance correlation, expected transition-to-transversion ratio T-s/T-v and Mendelian mismatch rates of known parent-offspring trios). 3. We observed large differences between reference-based and de novo approaches, the former generally calling more SNPs and reducing F-IS and T-s/T-v. Data completion levels showed little impact on most summary statistics, and FST estimates were robust across all pipelines. The site frequency spectrum was highly sensitive to the chosen approach as reflected in large variance of parameter estimates across demographic scenarios (single-population bottlenecks and isolation-with-migration model). Null expectations were best met by reference-based approaches, although contingent on the specific criteria. 4. We recommend that RAD-seq studies employ reference-based approaches to a closely related genome, and due to the high stochasticity associated with the pipeline advocate the use of multiple pipelines to ensure robust population genetic and demographic inferences.

High-density sex-specific linkage maps of a European tree frog (Hyla arborea) identify the sex chromosome without information on offspring sex.

Identifying homology between sex chromosomes of different species is essential to understanding the evolution of sex determination. Here, we show that the identity of a homomorphic sex chromosome pair can be established using a linkage map, without information on offspring sex. By comparing sex-specific maps of the European tree frog Hyla arborea, we find that the sex chromosome (linkage group 1) shows a threefold difference in marker number between the male and female maps. In contrast, the number of markers on each autosome is similar between the two maps. We also find strongly conserved synteny between H. arborea and Xenopus tropicalis across 200 million years of evolution, suggesting that the rate of chromosomal rearrangement in anurans is low. Finally, we show that recombination in males is greatly reduced at the centers of large chromosomes, consistent with previous cytogenetic findings. Our research shows the importance of high-density linkage maps for studies of recombination, chromosomal rearrangement and the genetic architecture of ecologically or economically important traits.