Hybridization and postzygotic isolation promote reinforcement of male mating preferences in a diverse group of fishes with traditional sex roles

Riverscape genetics of the orangethroat darter complex

Freshwater darters belonging to the orangethroat darter species complex, or Ceasia, are widely distributed in the Central and Southern United States, with ranges that span both glaciated and unglaciated regions. Up to 15 species have been recognized in the complex, with one, Etheostoma spectabile, having a widespread northern distribution and another, Etheostoma pulchellum, having a sizeable southern distribution. The other species in the complex have much more restricted distributions in unglaciated regions of the Central Highlands. We sampled 384 darters from 52 sites covering much of the range of Ceasia and evaluated patterns of genetic diversity, genetic structure, and pre- and post-glacial patterns of range contraction and expansion. We anticipated finding much stronger signals of genetic differentiation and diversification in unglaciated regions, given the higher species diversity and levels of endemism reported there. Surprisingly, microsatellite genotyping revealed two well-differentiated genetic clusters of E. spectabile in samples from glaciated regions, one confined to the Illinois River basin and another found in the Wabash drainage and Great Lakes tributaries. This suggests that there was expansion from two isolated glacial refugia, with little subsequent post-glacial gene flow. Fish collected from throughout the unglaciated region were less genetically differentiated. Fish assigned to Etheostoma burri and Etheostoma uniporum based on collection sites and morphological characters were not genetically differentiated from E. spectabile samples from the region. Hybridization and introgression occurring in the Central Highlands may confound genetic delineation of species in this region of high endemism and diversity.

Male recognition of conspecific female chemical cues in a diverse clade of freshwater fishes

Where orangethroat darters (Etheostoma: Ceasia) and rainbow darters (Etheostoma caeruleum) co-occur, males prefer conspecific over heterospecific females. The cues males use to identify conspecific females remain unclear. We conducted behavioral trials to ask whether chemical cues function in conspecific recognition. We found that males from three orangethroat darter species preferentially associate with female scent over a control. Our results support the use of olfaction in conspecific identification in the orangethroat clade and contribute to our understanding of signals that may facilitate species recognition and underlie the evolution of behavioral isolation.

Variable rates of hybridization among contact zones between a pair of topminnow species, Fundulus notatus and F. olivaceus

Pairs of species that exhibit broadly overlapping distributions, and multiple geographically isolated contact zones, provide opportunities to investigate the mechanisms of reproductive isolation. Such naturally replicated systems have demonstrated that hybridization rates can vary substantially among populations, raising important questions about the genetic basis of reproductive isolation. The topminnows, Fundulus notatus and F. olivaceus, are reciprocally monophyletic, and co-occur in drainages throughout much of the central and southern United States. Hybridization rates vary substantially among populations in isolated drainage systems. We employed genome-wide sampling to investigate geographic variation in hybridization, and to assess the possible importance of chromosome fusions to reproductive isolation among nine separate contact zones. The species differ by chromosomal rearrangements resulting from Robertsonian (Rb) fusions, so we hypothesized that Rb fusion chromosomes would serve as reproductive barriers, exhibiting steeper genomic clines than the rest of the genome. We observed variation in hybridization dynamics among drainages that ranged from nearly random mating to complete absence of hybridization. Contrary to predictions, our use of genomic cline analyses on mapped species-diagnostic SNP markers did not indicate consistent patterns of variable introgression across linkage groups, or an association between Rb fusions and genomic clines that would be indicative of reproductive isolation. We did observe a relationship between hybridization rates and population phylogeography, with the lowest rates of hybridization tending to be found in populations inferred to have had the longest histories of drainage sympatry. Our results, combined with previous studies of contact zones between the species, support population history as an important factor in explaining variation in hybridization rates.

Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Rapture-ready darters: Choice of reference genome and genotyping method (whole-genome or sequence capture) influence population genomic inference in Etheostoma

Researchers studying nonmodel organisms have an increasing number of methods available for generating genomic data. However, the applicability of different methods across species, as well as the effect of reference genome choice on population genomic inference, remain difficult to predict in many cases. We evaluated the impact of data type (whole-genome vs. reduced representation) and reference genome choice on data quality and on population genomic and phylogenomic inference across several species of darters (subfamily Etheostomatinae), a highly diverse radiation of freshwater fish. We generated a high-quality reference genome and developed a hybrid RADseq/sequence capture (Rapture) protocol for the Arkansas darter (Etheostoma cragini). Rapture data from 1,900 individuals spanning four darter species showed recovery of most loci across darter species at high depth and consistent estimates of heterozygosity regardless of reference genome choice. Loci with baits spanning both sides of the restriction enzyme cut site performed especially well across species. For low-coverage whole-genome data, choice of reference genome affected read depth and inferred heterozygosity. For similar amounts of sequence data, Rapture performed better at identifying fine-scale genetic structure compared to whole-genome sequencing. Rapture loci also recovered an accurate phylogeny for the study species and demonstrated high phylogenetic informativeness across the evolutionary history of the genus Etheostoma. Low cost and high cross-species effectiveness regardless of reference genome suggest that Rapture and similar sequence capture methods may be worthwhile choices for studies of diverse species radiations.

Reinforcement in the banded darter Etheostoma zonale: The effect of sex and sympatry on preferences

Reinforcement occurs when selection against hybrid offspring strengthens behavioral isolation between parental species and may be an important factor in speciation. Theoretical models and experimental evidence indicate that both female and male preferences can be strengthened upon secondary contact via reinforcement. However, the question remains whether this process is more likely to affect the preferences of one sex or the other. Males of polygynous species are often predicted to exhibit weaker preferences than females, potentially limiting the ability for reinforcement to shape male preferences. Yet, in darters (Percidae: Etheostoma), male preference for conspecific mates appears to arise before female preferences during the early stages of allopatric speciation, and research suggests that male, but not female, preferences become reinforced upon secondary contact. In the current study, we aimed to determine whether the geographically widespread darter species Etheostoma zonale exhibits a signature of reinforcement, by comparing the strength of preference for conspecific mates between populations that are sympatric and allopatric with respect to a close congener, E. barrenense. We examined the strength of preference for conspecifics for males and females separately to determine whether the preferences of one or both sexes have been strengthened by reinforcement. Our results show that both sexes of E. zonale from sympatric populations exhibit stronger conspecific preferences than E. zonale from allopatric populations, but that female preferences appear to be more strongly reinforced than male preferences. Results therefore suggest that reinforcement of female preferences may promote behavioral isolation upon secondary contact, even in a genus that is characterized by pervasive male mate choice.

Genomic Resources for Darters (Percidae: Etheostominae) Provide Insight into Postzygotic Barriers Implicated in Speciation

Comparative genomic approaches are increasingly being used to study the evolution of reproductive barriers in nonmodel species. Although numerous studies have examined prezygotic isolation in darters (Percidae), investigations into postzygotic barriers have remained rare due to long generation times and a lack of genomic resources. Orangethroat and rainbow darters naturally hybridize and provide a remarkable example of male-driven speciation via character displacement. Backcross hybrids suffer from high mortality, which appears to promote behavioral isolation in sympatry. To investigate the genomic architecture of postzygotic isolation, we used Illumina and PacBio sequencing to generate a chromosome-level, annotated assembly of the orangethroat darter genome and high-density linkage maps for orangethroat and rainbow darters. We also analyzed genome-wide RADseq data from wild-caught adults of both species and laboratory-generated backcrosses to identify genomic regions associated with hybrid incompatibles. Several putative chromosomal translocations and inversions were observed between orangethroat and rainbow darters, suggesting structural rearrangements may underlie postzygotic isolation. We also found evidence of selection against recombinant haplotypes and transmission ratio distortion in backcross hybrid genomes, providing further insight into the genomic architecture of genetic incompatibilities. Notably, regions with high levels of genetic divergence between species were enriched for genes associated with developmental and meiotic processes, providing strong candidates for postzygotic isolating barriers. These findings mark significant contributions to our understanding of the genetic basis of reproductive isolation between species undergoing character displacement. Furthermore, the genomic resources presented here will be instrumental for studying speciation in darters, the most diverse vertebrate group in North America.

Intergeneric hybrids inform reproductive isolating barriers in the Antarctic icefish radiation

Interspecific hybridization or barriers to hybridization may have contributed to the diversification of Antarctic icefishes (Channichthyidae), but data supporting these hypotheses is scarce. To understand the potential for hybridization and to investigate reproductive isolating mechanisms among icefish species, we performed in vitro fertilization experiments using eggs from a female blackfin icefish Chaenocephalus aceratus and sperm from a male of another genera, the ocellated icefish Chionodraco rastrospinosus. Sequencing of genomic and mitochondrial DNA confirmed the intergeneric hybrid nature of resulting embryos which successfully developed and hatched as active larvae at about four and a half months during the Antarctic winter. This result demonstrates the compatibility of gametes of these two species and the viability of resulting zygotes and larvae. Due to logistic constraints and the slow developmental rate of icefishes, we could not test for long-term hybrid viability, fertility, fitness, or hybrid breakdown. Analysis of our fishing records and available literature, however, suggests that the strongest barriers to hybridization among parapatric icefish species are likely to be behavioral and characterized by assortative mating and species-specific courtship and nesting behaviors. This conclusion suggests that, in long-lived fish species with late sexual maturity and high energetic investment in reproduction like icefishes, pre-mating barriers are energetically more efficient than post-mating barriers to prevent hybridization.