High homozygosity of inversions in sunflower species largely averts accumulation of deleterious mutations

Recombination is critical both for accelerating adaptation and for the purging of deleterious mutations. Chromosomal inversions can act as recombination modifiers that suppress local recombination and, thus, are predicted to accumulate such mutations. In this study, we investigated patterns of recombination, transposable element abundance and coding sequence evolution across the genomes of 1,445 individuals from three sunflower species, as well as within nine inversions segregating within species. We also analyzed the effects of inversion genotypes on 87 phenotypic traits to test for overdominance. We found significant negative correlations of long terminal repeat retrotransposon abundance and deleterious mutations with recombination rates across the genome in all three species. However, we failed to detect an increase in these features in the inversions, except for a modest increase in the proportion of stop codon mutations in several very large or rare inversions. Moreover, there was little evidence of phenotypic overdominance in inversion heterozygotes, consistent with observations of minimal deleterious load. On the other hand, significantly greater load was observed for inversions in populations polymorphic for a given inversion compared to populations monomorphic for one of the arrangements, suggesting that the local state of inversion polymorphism affects deleterious load. These seemingly contradictory results can be explained by the geographic structuring and consequent excess homozygosity of inversions in wild sunflowers. Inversions contributing to local adaptation often exhibit geographic structure; such inversions represent ideal recombination modifiers, acting to facilitate adaptive divergence with gene flow, while largely averting the accumulation of deleterious mutations due to recombination suppression.

Genomic signatures of isolation, hybridization, and selection during speciation of island finches

Sister species occurring sympatrically on islands are rare and offer unique opportunities to understand how speciation can proceed in the face of gene flow. The São Tomé grosbeak is a massive-billed, 'giant' finch endemic to the island of São Tomé in the Gulf of Guinea, where it has diverged from its co-occurring sister species the Príncipe seedeater, an average-sized finch that also inhabits two neighbouring islands. Here, we show that the grosbeak carries a large number of unique alleles different from all three Príncipe seedeater populations, but also shares many alleles with the sympatric São Tomé population of the seedeater, a genomic signature signifying divergence in isolation as well as subsequent introgressive hybridization. Furthermore, genomic segments that remain unique to the grosbeak are situated close to genes, including genes that determine bill morphology, suggesting the preservation of adaptive variation through natural selection during divergence with gene flow. This study reveals a complex speciation process whereby genetic drift, introgression, and selection during periods of isolation and secondary contact all have shaped the diverging genomes of these sympatric island endemic finches.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca−nimia, crecca−carolinensis, and carolinensis−flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ~0.1 - 0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca−nimia), parapatric (crecca−carolinensis), and (mostly) allopatric (carolinensis−flavirostris). Ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems like this.

Faster-X adaptive differentiation under divergence-with-gene-flow in diploids and haplodiploids

Empirical data from diverse taxa indicate that the hemizygous portions of the genome (X/Z chromosomes) evolve more rapidly than their diploid counterparts. Faster-X theory predicts increased rates of adaptive substitutions between isolated species, yet little is known about species experiencing gene flow. Here we investigate how hemizygosity impacts genome-wide patterns of differentiation during adaptive divergence with gene flow, combining simulations under isolation-with-migration models, a meta-analysis of autosomes and sex-chromosomes from diverse taxa, and analysis of haplodiploid species. First, using deterministic and stochastic simulations, we show that elevated differentiation at hemizygous loci occurs when there is gene flow, irrespective of dominance. This faster-X adaptive differentiation stems from more efficient selection resulting in reduced probability of losing the beneficial allele, greater migration-selection threshold, greater allele frequency differences at equilibrium, and a faster time to equilibrium. Second, by simulating neutral variation linked to selected loci, we show that faster-X differentiation affects linked variation due to reduced opportunities for recombination between locally adaptive and maladaptive immigrant haplotypes. Third, after correcting for expected differences in effective population size, we find that most taxon pairs (24 out of 28) exhibit faster-X differentiation in the meta-analysis. Finally, using a novel approach combining demographic modeling and simulations, we found evidence for faster-X differentiation in haplodiploid pine-feeding hymenopteran species adapted to different host plants. Together, our results indicate that divergent selection with gene flow can lead to higher differentiation at selected and linked variation in hemizygous loci (i.e., faster-X adaptive differentiation), both in X/Z-chromosomes and haplodiploid species.

Repeated colonization of caves leads to phenotypic convergence in catfishes (Siluriformes: Trichomycterus) at a small geographical scale

Across various animal groups, adaptation to the challenging conditions of cave environments has resulted in convergent evolution. We document a Neotropical cavefish system with ample potential to study questions related to convergent adaptation to cave environments at the population level. In the karstic region of the Andes of Santander, Colombia, cave-dwelling catfish in the genus Trichomycterus exhibit variable levels of reduction of eyes and body pigmentation relative to surface congeners. We tested whether cave-dwelling, eye-reduced, depigmented Trichomycterus from separate caves in Santander were the result of a single event of cave colonization and subsequent dispersal, or of multiple colonizations to caves by surface ancestors followed by phenotypic convergence. Using mitochondrial DNA sequences to reconstruct phylogenetic relationships, we found that caves in this region have been colonized independently by two separate clades. Additional events of cave colonization – and possibly recolonization of surface streams – may have occurred in one of the clades, where surface and cave-dwelling populations exhibit shallow differentiation, suggesting recent divergence or divergence with gene flow. We also identify potentially undescribed species and likely problems with the circumscription of named taxa. The system appears promising for studies on a wide range of ecological and evolutionary questions.

Selection on a pleiotropic color gene block underpins early differentiation between two warbler species

When one species gradually splits into two, divergent selection on specific traits can cause peaks of differentiation in the genomic regions encoding those traits. Whether speciation is initiated by strong selection on a few genomic regions with large effects or by more diffused selection on many regions with small effects remains controversial. Differentiated phenotypes between differentiating lineages are commonly involved in reproductive isolation, thus their genetic underpinnings are key to the genomics architecture of speciation. When two species hybridize, recombination over multiple generations can help reveal the genetic regions responsible for the differentiated phenotypes against a genomic background that has been homogenized via backcrossing and introgression. We used admixture mapping to investigate genomic differentiation and the genetic basis of differentiated plumage features (relative melanin and carotenoid pigment) between hybridizing sister species in the early stage of speciation: Townsend’s (Setophaga townsendi) and Hermit warblers (S. occidentalis). We found a few narrow and dispersed divergent regions between allopatric parental populations, consistent with the ‘divergence with gene flow’ model of speciation. One of the divergent peaks involves three genes known to affect pigmentation: ASIP, EIF2S2, and RALY (the ASIP-RALY gene block). After controlling for population substructure, we found that a single nucleotide polymorphism (SNP) inside the intron of RALY displays a strong pleiotropic association with cheek, crown, and breast coloration. In addition, we detect selection on the ASIP-RALY gene block, as the geographic cline of the RALY marker of this gene block has remained narrower than the plumage cline, which remained narrower than expected under neutral diffusion over two decades. Despite extensive gene flow between these species across much of the genome, the selection on ASIP-RALY gene block maintains stable genotypic and plumage difference between species allowing further differentiation to accumulate via linkage to its flanking genetic region or linkage-disequilibrium genome-wide.

Shallow genetic divergence and distinct phenotypic differences between two Andean hummingbirds: Speciation with gene flow?

Ecological speciation can proceed despite genetic interchange when selection counteracts the homogenizing effects of migration. We tested predictions of this divergence-with-gene-flow model in Coeligena helianthea and C. bonapartei, 2 parapatric Andean hummingbirds with marked plumage divergence. We sequenced putatively neutral markers (mitochondrial DNA [mtDNA] and nuclear ultraconserved elements [UCEs]) to examine genetic structure and gene flow, and a candidate gene (MC1R) to assess its role underlying divergence in coloration. We also tested the prediction of Gloger’s rule that darker forms occur in more humid environments, and examined morphological variation to assess adaptive mechanisms potentially promoting divergence. Genetic differentiation between species was low in both ND2 and UCEs. Coalescent estimates of migration were consistent with divergence with gene flow, but we cannot reject incomplete lineage sorting reflecting recent speciation as an explanation for patterns of genetic variation. MC1R variation was unrelated to phenotypic differences. Species did not differ in macroclimatic niches but were distinct in morphology. Although we reject adaptation to variation in macroclimatic conditions as a cause of divergence, speciation may have occurred in the face of gene flow driven by other ecological pressures or by sexual selection. Marked phenotypic divergence with no neutral genetic differentiation is remarkable for Neotropical birds, and makes C. helianthea and C. bonapartei an appropriate system in which to search for the genetic basis of species differences employing genomics.

The multiple population genetic and demographic routes to islands of genomic divergence

Abstract1. The way that organisms diverge into reproductively isolated species is a major question in biology. The recent accumulation of genomic data provides promising opportunities to understand the genomic landscape of divergence, which describes the distribution of differences across genomes. Genomic areas of unusually high differentiation have been called genomic islands of divergence. Their formation has been attributed to a variety of mechanisms, but a prominent hypothesis is that they result from divergent selection over a small portion of the genome, with surrounding areas homogenised by gene flow. Such islands have often been interpreted as being associated with divergence with gene flow. However other mechanisms related to genetic architecture and population history can also contribute to the formation of genomic islands of divergence.2. We currently lack a quantitative framework to examine the dynamics of genomic landscapes under the complex and nuanced conditions that are found in natural systems. Here, we develop an individual-based simulation to explore the dynamics of diverging genomes under various scenarios of gene flow, selection and genotype-phenotype maps.3. Our modelling results are consistent with empirical observations demonstrating the formation of genomic islands under genetic isolation. Importantly, we have quantified the range of conditions that produce genomic islands. We demonstrate that the initial level of genetic diversity, drift, time since divergence, linkage disequilibrium, strength of selection and gene flow are all important factors that can influence the formation of genomic islands. Because the accumulation of genomic differentiation over time tends to erode the signal of genomic islands, genomic islands are more likely to be observed in recently divergent taxa, although not all recently diverged taxa will necessarily exhibit islands of genomic divergence. Gene flow primarily slows the swamping of islands of divergence with time.4. By using this framework, further studies may explore the relative influence of particular suites of events that contribute to the emergence of genomic islands under sympatric, parapatric and allopatric conditions. This approach represents a novel tool to explore quantitative expectations of the speciation process, and should prove useful in elucidating past and projecting future genomic evolution of any taxa.