Intra-Genomic Variation in the Ribosomal Repeats of Nematodes

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.     

Use of 16S rRNA and rpoB Genes as Molecular Markers for Microbial Ecology Studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such “heterogeneity hot spots” occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Multilocus patterns of genetic variation across Cryptosporidium species suggest balancing selection at the gp60 locus

Cryptosporidium is an apicomplexan protozoan that lives in most vertebrates, including humans. Its gp60 gene is functionally involved in its attachment to host cells, and its high level of genetic variation has made it the reference marker for sample typing in epidemiological studies. To understand the origin of such high diversity and to determine the extent to which this classification applies to the rest of the genome, we analysed the patterns of variation at gp60 and nine other nuclear loci in isolates of three Cryptosporidium species. Most loci showed low genetic polymorphism (π_S <1%) and similar levels of between‐species divergence. Contrastingly, gp60 exhibited very different characteristics: (i) it was nearly ten times more variable than the other loci; (ii) it displayed a significant excess of polymorphisms relative to between‐species differences in a maximum‐likelihood Hudson–Kreitman–Aguadé test; (iii) gp60 subtypes turned out to be much older than the species they were found in; and (iv) showed a significant excess of polymorphic variants shared across species from random expectations. These observations suggest that this locus evolves under balancing selection and specifically under negative frequency‐dependent selection (FDS). Interestingly, genetic variation at the other loci clusters very well within the groups of isolates defined by gp60 subtypes, which may provide new tools to understand the genome‐wide patterns of genetic variation of the parasite in the wild. These results suggest that gp60 plays an active and essential role in the life cycle of the parasite and that genetic variation at this locus might be essential for the parasite's long‐term success.     

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome‐scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well‐defined categories of locus types (e.g., neutral vs. highly selected). Our multi‐locus analysis reveals cross‐genomic variation in selective constraints on gene flow and locus‐specific flexibility in the permeability of the interspecies membrane.     

Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such "heterogeneity hot spots" occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Phylogenetics of New World Astragalus: Screening of novel nuclear loci for the reconstruction of phylogenies at low taxonomic levels

This study explores methods to use information gathered from genomics technology to understand evolutionary relationships in the hyperdiverse legume group Neo-Astragalus. These species inhabit deserts and mountains of North and South America, and even though the monophyly of the group is well established, relationships within it are still poorly understood. Plastid genes, commonly used to infer phylogenies in plants, are usually not useful for closely related taxa because of low levels of genetic variation. TheMedicago truncatula genome project provided a suite of candidate nuclear loci with high levels of variation that might prove suitable for low-level phylogenetics. This paper reports the development of methods for screening a large number of these nuclear loci, and detailed analysis of four of them. Four different patterns of phylogenetic diversification occur in the loci sampled from these genomes ofAstragalus species. One locus (CNGC4) was single copy and could be directly used in phylogenetic analyses. Two loci (ARG10 and FENR) showed patterns strongly suggestive of duplication events in some taxa, and one locus (tRALS) has apparently undergone a cryptic duplication, making it very difficult to diagnose. Potential methods for using the information provided by these loci are discussed.     

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom‐feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre‐ and post‐mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi‐locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.     

Buchnera has changed flatmate but the repeated replacement of co‐obligate symbionts is not associated with the ecological expansions of their aphid hosts

Symbiotic associations with bacteria have facilitated important evolutionary transitions in insects and resulted in long‐term obligate interactions. Recent evidence suggests that these associations are not always evolutionarily stable and that symbiont replacement, and/or supplementation of an obligate symbiosis by an additional bacterium, has occurred during the history of many insect groups. Yet, the factors favouring one symbiont over another in this evolutionary dynamic are not well understood; progress has been hindered by our incomplete understanding of the distribution of symbionts across phylogenetic and ecological contexts. While many aphids are engaged into an obligate symbiosis with a single Gammaproteobacterium, Buchnera aphidicola, in species of the Lachninae subfamily, this relationship has evolved into a ‘ménage à trois’, in which Buchnera is complemented by a cosymbiont, usually Serratia symbiotica. Using deep sequencing of 16S rRNA bacterial genes from 128 species of Cinara (the most diverse Lachninae genus), we reveal a highly dynamic dual symbiotic system in this aphid lineage. Most species host both Serratia and Buchnera but, in several clades, endosymbionts related to Sodalis, Erwinia or an unnamed member of the Enterobacteriaceae have replaced Serratia. Endosymbiont genome sequences from four aphid species confirm that these coresident symbionts fulfil essential metabolic functions not ensured by Buchnera. We further demonstrate through comparative phylogenetic analyses that cosymbiont replacement is not associated with the adaptation of aphids to new ecological conditions. We propose that symbiont succession was driven by factors intrinsic to the phenomenon of endosymbiosis, such as rapid genome deterioration or competitive interactions between bacteria with similar metabolic capabilities.     

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole‐genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Conservatism and novelty in the genetic architecture of adaptation in Heliconius butterflies

Understanding the genetic architecture of adaptive traits has been at the centre of modern evolutionary biology since Fisher; however, evaluating how the genetic architecture of ecologically important traits influences their diversification has been hampered by the scarcity of empirical data. Now, high-throughput genomics facilitates the detailed exploration of variation in the genome-to-phenotype map among closely related taxa. Here, we investigate the evolution of wing pattern diversity in Heliconius, a clade of neotropical butterflies that have undergone an adaptive radiation for wing-pattern mimicry and are influenced by distinct selection regimes. Using crosses between natural wing-pattern variants, we used genome-wide restriction site-associated DNA (RAD) genotyping, traditional linkage mapping and multivariate image analysis to study the evolution of the architecture of adaptive variation in two closely related species: Heliconius hecale and H. ismenius. We implemented a new morphometric procedure for the analysis of whole-wing pattern variation, which allows visualising spatial heatmaps of genotype-to-phenotype association for each quantitative trait locus separately. We used the H. melpomene reference genome to fine-map variation for each major wing-patterning region uncovered, evaluated the role of candidate genes and compared genetic architectures across the genus. Our results show that, although the loci responding to mimicry selection are highly conserved between species, their effect size and phenotypic action vary throughout the clade. Multilocus architecture is ancestral and maintained across species under directional selection, whereas the single-locus (supergene) inheritance controlling polymorphism in H. numata appears to have evolved only once. Nevertheless, the conservatism in the wing-patterning toolkit found throughout the genus does not appear to constrain phenotypic evolution towards local adaptive optima.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Phylogenetic Utility of Histone H3 Intron Sequences in the Perennial Relatives of Soybean (Glycine:Leguminosae)

Histone H3 loci form a large multigene family in most plant species. InGlycine,some of these loci possess introns, whose sequences can provide characters for assessing phylogenetic relationships among species of the genus. Phylogenetic analyses of two closely related H3-B loci revealed a complex evolutionary pattern, producing trees from which species relationships could not be inferred readily. The single H3-D locus, in contrast, provided data suitable for the construction of gene trees whose topologies were sufficiently similar to other hypotheses of relationships within the subgenusGlycineto give confidence that evolution at this locus is tracking species phylogenies. H3-D topologies identified several of the same groupings found in previous phylogenetic studies using the chloroplast genome. However, histone H3-D and chloroplast genome data sets were in other respects incongruent, as revealed by both topological differences and numerical measures of congruence. The principal difference involvedGlycine falcata,whose chloroplast genome belongs to one of the three strongly supported clades in the subgenus, but whose histone H3-D allele was sister to those of the remaining members of the subgenus. The H3-D topology is more in keeping with the morphologically, ecologically, and genetically divergent nature of this species. The H3-D locus appears to be a useful source of phylogenetic characters for interspecific studies inGlycine,providing resolution among taxa whose relationships were unresolved in previous studies.     

Characterization of the Two Intra-Individual Sequence Variants in the 18S rRNA Gene in the Plant Parasitic Nematode,Rotylenchulus reniformis

The 18S rRNA gene is fundamental to cellular and organismal protein synthesis and because of its stable persistence through generations it is also used in phylogenetic analysis among taxa. Sequence variation in this gene within a single species is rare, but it has been observed in few metazoan organisms. More frequently it has mostly been reported in the non-transcribed spacer region. Here, we have identified two sequence variants within the near full coding region of 18S rRNA gene from a single reniform nematode (RN) Rotylenchulus reniformis labeled as reniform nematode variant 1 (RN_VAR1) and variant 2 (RN_VAR2). All sequences from three of the four isolates had both RN variants in their sequences; however, isolate 13B had only RN variant 2 sequence. Specific variable base sites (96 or 5.5%) were found within the 18S rRNA gene that can clearly distinguish the two 18S rDNA variants of RN, in 11 (25.0%) and 33 (75.0%) of the 44 RN clones, for RN_VAR1 and RN_VAR2, respectively. Neighbor-joining trees show that the RN_VAR1 is very similar to the previously existing R. reniformis sequence in GenBank, while the RN_VAR2 sequence is more divergent. This is the first report of the identification of two major variants of the 18S rRNA gene in the same single RN, and documents the specific base variation between the two variants, and hypothesizes on simultaneous co-existence of these two variants for this gene.     

Concerted and birth-and-death evolution of 26S ribosomal DNA in Camellia L.

Background and AimsThe ribosomal DNA (rDNA) gene family, encoding ribosomal RNA (rRNA), has long been regarded as an archetypal example illustrating the model of concerted evolution. However, controversy is arising, as rDNA in many eukaryotic species has been proved to be polymorphic. Here, a metagenomic strategy was applied to detect the intragenomic polymorphism as well as the evolutionary patterns of 26S rDNA across the genus Camellia.MethodsDegenerate primer pairs were designed to amplify the 26S rDNA fragments from different Camellia species. The amplicons were then paired-end sequenced on the Illumina MiSeq platform.Key ResultsAn extremely high level of rDNA polymorphism existed universally in Camellia. However, functional rDNA was still the major component of the family, and was relatively conserved among different Camellia species. Sequence variations mainly came from rRNA pseudogenes and favoured regions that are rich in GC. Specifically, some rRNA pseudogenes have existed in the genome for a long time, and have even experienced several expansion events, which has greatly enriched the abundance of rDNA polymorphism.Conclusions Camellia represents a group in which rDNA is subjected to a mixture of concerted and birth-and-death evolution. Some rRNA pseudogenes may still have potential functions. Conversely, when released from selection constraint, they can evolve in the direction of decreasing GC content and structural stability through a methylation-induced process, and finally be eliminated from the genome.     

Ready on arrival: Standing variation at a chromosomal inversion contributes to rapid adaptation in an invasive marine crab

In this issue of Molecular Ecology, Tepolt et al. (2021) illustrate how the genetic architecture of adaptation and life history influence invasive success. A marvel of many invasive species is that they are incredibly successful despite evolutionary expectations that they will have low adaptive potential and suffer inbreeding depression due to initially small founding population sizes. Determining the combinations of ecoevolutionary factors that permit this apparent “genetic paradox of invasions” is an ongoing endeavour of invasive species research. Tepolt et al. (2021) study the European green crab in its invasive range on the North American west coast. Following a single introduction into California, this crab quickly spread across a wide latitude gradient, despite low diversity in the original founding population. Adaptation of this crab to clinal variation in temperature appeared largely driven by an inferred chromosomal inversion. This inversion exists as a balanced polymorphism in the European home range of green crabs and is associated with thermal tolerance. Tepolt et al. (2021) therefore demonstrate that adaptive evolution post introduction need not be impeded by bottlenecks if variation at key parts of the genome is available and can be maintained in introduced populations. Moreover, Tepolt et al. (2021) show how chromosomal inversions acting as large-effect loci might facilitate rapid responses to selection in introduced populations.     

Differential divergence in autosomes and sex chromosomes is associated with intra‐island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome‐wide markers to investigate the diversification of the Reunion grey white‐eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping‐by‐sequencing and pooled RAD‐seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z‐linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.     

Highly variable reproductive isolation among pairs of Catostomus species

Hybridization between diverged taxa tests the strength of reproductive isolation and can therefore reveal mechanisms of reproductive isolation. However, it remains unclear how consistent reproductive isolation is across species' ranges and to what extent reproductive isolation might remain polymorphic as species diverge. To address these questions, we compared outcomes of hybridization across species pairs of Catostomus fishes in three rivers in the Upper Colorado River basin, where an introduced species, C. commersoni, hybridizes with at least two native species, C. discobolus and C. latipinnis. We observed substantial heterogeneity in outcomes of hybridization, both between species pairs and across geographically separate rivers within each species pair. We also observed hybridization of additional related species with our focal species, suggesting that reproductive isolation in this group involves interactions of multiple evolutionary and ecological factors. These findings suggest that a better understanding of the determinants of variation in reproductive isolation is needed and that studies of reproductive isolation in hybrids should consider how the dynamics and mechanisms of reproductive isolation vary over ecological space and over evolutionary time. Our results also have implications for the conservation and management of native catostomids in the Colorado River basin. Heterogeneity in outcomes of hybridization suggests that the threat posed by hybridization and genetic introgression to the persistence of native species probably varies with extent of reproductive isolation, both across rivers and across species pairs.     

Lack of gene flow: Narrow and dispersed differentiation islands in a triplet of Leptidea butterfly species

Genome scans in recently separated species can inform on molecular mechanisms and evolutionary processes driving divergence. Large‐scale polymorphism data from multiple species pairs are also key to investigate the repeatability of divergence—whether radiations tend to show parallel responses to similar selection pressures and/or underlying molecular forces. Here, we used whole‐genome resequencing data from six wood white (Leptidea sp.) butterfly populations, representing three closely related species with karyomorph variation, to infer the species' demographic history and characterize patterns of genomic diversity and differentiation. The analyses supported previously established species relationships, and there was no evidence for postdivergence gene flow. We identified significant intraspecific genetic structure, in particular between karyomorph extremes in the wood white (L. sinapis)—a species with a remarkable chromosome number cline across the distribution range. The genomic landscapes of differentiation were erratic, and outlier regions were narrow and dispersed. Highly differentiated (F_ST) regions generally had low genetic diversity (θ_π), but increased absolute divergence (D_XY) and excess of rare frequency variants (low Tajima's D). A minority of differentiation peaks were shared across species and population comparisons. However, highly differentiated regions contained genes with overrepresented functions related to metabolism, response to stimulus and cellular processes, indicating recurrent directional selection on a specific set of traits in all comparisons. In contrast to the majority of genome scans in recently diverged lineages, our data suggest that divergence landscapes in Leptidea have been shaped by directional selection and genetic drift rather than stable recombination landscapes and/or introgression.