Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.     

Hybridization occurs between Drosophila simulans and D. sechellia in the Seychelles archipelago

Drosophila simulans and D. sechellia are sister species that serve as a model to study the evolution of reproductive isolation. While D. simulans is a human commensal that has spread all over the world, D. sechellia is restricted to the Seychelles archipelago and is found to breed exclusively on the toxic fruit of Morinda citrifolia. We surveyed the relative frequency of males from these two species in a variety of substrates found on five islands of the Seychelles archipelago. We sampled different fruits and found that putative D. simulans can be found in a variety of substrates, including, surprisingly, M. citrifolia. Putative D. sechellia was found preferentially on M. citrifolia fruits, but a small proportion was found in other substrates. Our survey also shows the existence of putative hybrid males in areas where D. simulans is present in Seychelles. The results from this field survey support the hypothesis of current interbreeding between these species in the central islands of Seychelles and open the possibility for fine measurements of admixture between these two Drosophila species to be made.     

Genetic divergence of the Mesoamerican azure‐crowned hummingbird (Amazilia cyanocephala,Trochilidae) across the Motagua‐Polochic‐Jocotán fault system

The cloud forests of Mesoamerica are highly endangered habitats and the existence of narrowly distributed cryptic endemics will increase the number of taxa at potential risk of extinction. Here, we investigate genetic divergence between populations of the azure‐crowned hummingbird (Amazilia cyanocephala), a species complex of endemic hummingbirds to the montane forests of Mesoamerica, by analysing DNA sequences of four mitochondrial markers, morphological data and ecological niche modelling. Our results revealed the presence of two mtDNA lineages corresponding to subspecies A. c. cyanocephala distributed from Tamaulipas to Chiapas in Mexico and Amazilia c. guatemalensis distributed from southern Chiapas to Guatemala. The lineage split can be explained as a consequence of relative isolation of the populations in the different mountain ranges separated by the Motagua‐Polochic‐Jocotán fault system and corresponds to differences in morphology and to the lack of overlap in environmental space between subspecies. The divergence time estimates do not support the proposed model of a highly constrained temporal window at the end of the Pliocene as divergence at this barrier between cyanocephala and guatemalensis and splits of other bird taxa occurred during the Pleistocene.     

1.688 g/cm3 satellite‐related repeats: a missing link to dosage compensation and speciation

Despite the important progress that has been made on dosage compensation (DC), a critical link in our understanding of the X chromosome recognition mechanisms is still missing. Recent studies in Drosophila indicate that the missing link could be a family of DNA repeats populating the euchromatin of the X chromosome. In this opinion article, I discuss how these findings add a new fresh twist on the DC problem. In the following sections, I first summarize our understanding of DC in Drosophila and integrate these recent discoveries into our knowledge of the X chromosome recognition problem. Next, I introduce a model according to which, 1.688 g/cm³ satellite‐related (SR) repeats would be the primary recognition elements for the dosage compensation complex. Contrary to the current belief, I suggest that the DC system in Drosophila is not conserved and static, but it is continuously co‐evolving with the target SR repeats. The potential role of the SR repeats in hybrid incompatibilities and speciation is also discussed.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Interploidy hybridization in sympatric zones: the formation of Epidendrum fulgens × E. puniceoluteum hybrids (Epidendroideae,Orchidaceae)

Interspecific hybridization is a primary cause of extensive morphological and chromosomal variation and plays an important role in plant species diversification. However, the role of interploidal hybridization in the formation of hybrid swarms is less clear. Epidendrum encompasses wide variation in chromosome number and lacks strong premating barriers, making the genus a good model for clarifying the role of chromosomes in postzygotic barriers in interploidal hybrids. In this sense, hybrids from the interploidal sympatric zone between E. fulgens (2n = 2x = 24) and E. puniceoluteum (2n = 4x = 56) were analyzed using cytogenetic techniques to elucidate the formation and establishment of interploidal hybrids. Hybrids were not a uniform group: two chromosome numbers were observed, with the variation being a consequence of severe hybrid meiotic abnormalities and backcrossing with E. puniceoluteum. The hybrids were triploids (2n = 3x = 38 and 40) and despite the occurrence of enormous meiotic problems associated with triploidy, the hybrids were able to backcross, producing successful hybrid individuals with broad ecological distributions. In spite of the nonpolyploidization of the hybrid, its formation is a long‐term evolutionary process rather than a product of a recent disturbance, and considering other sympatric zones in Epidendrum, these events could be recurrent.     

Secondary contact zones of closely‐related Erebia butterflies overlap with narrow phenotypic and parasitic clines

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.     

Reticulate evolution and sea‐level fluctuations together drove species diversification of slipper orchids (Paphiopedilum) in South‐East Asia

South‐East Asia covers four of the world's biodiversity hotspots, showing high species diversity and endemism. Owing to the successive expansion and contraction of distribution and the fragmentation by geographical barriers, the tropical flora greatly diversified in this region during the Tertiary, but the evolutionary tempo and mode of species diversity remain poorly investigated. Paphiopedilum, the largest genus of slipper orchids comprising nearly 100 species, is mainly distributed in South‐East Asia, providing an ideal system for exploring how plant species diversity was shaped in this region. Here, we investigated the evolutionary history of this genus with eight cpDNA regions and four low‐copy nuclear genes. Discordance between gene trees and network analysis indicates that reticulate evolution occurred in the genus. Ancestral area reconstruction suggests that vicariance and long‐distance dispersal together led to its current distribution. Diversification rate variation was detected and strongly correlated with the species diversity in subg. Paphiopedilum (~80 species). The shift of speciation rate in subg. Paphiopedilum was coincident with sea‐level fluctuations in the late Cenozoic, which could have provided ecological opportunities for speciation and created bridges or barriers for gene flow. Moreover, some other factors (e.g. sympatric distribution, incomplete reproductive barriers and clonal propagation) might also be advantageous for the formation and reproduction of hybrid species. In conclusion, our study suggests that the interplay of reticulate evolution and sea‐level fluctuations has promoted the diversification of the genus Paphiopedilum and sheds light into the evolution of Orchidaceae and the historical processes of plant species diversification in South‐East Asia.     

Developmental effects of visual environment on species‐assortative mating preferences in Lake Victoria cichlid fish

Local adaptation can be a potent force in speciation, with environmental heterogeneity leading to niche specialization and population divergence. However, local adaption often requires nonrandom mating to generate reproductive isolation. Population divergence in sensory properties can be particularly consequential in speciation, affecting both ecological adaptation and sexual communication. Pundamilia pundamila and Pundamilia nyererei are two closely related African cichlid species that differ in male coloration, blue vs. red. They co‐occur at rocky islands in southern Lake Victoria, but inhabit different depth ranges with different light environments. The species differ in colour vision properties, and females exert species‐specific preferences for blue vs. red males. Here, we investigated the mechanistic link between colour vision and preference, which could provide a rapid route to reproductive isolation. We tested the behavioural components of this link by experimentally manipulating colour perception – we raised both species and their hybrids under light conditions mimicking shallow and deep habitats – and tested female preference for blue and red males under both conditions. We found that rearing light significantly affected female preference: shallow‐reared females responded more strongly to P. pundamilia males and deep‐reared females favoured P. nyererei males – implying that visual development causally affects mate choice. These results are consistent with sensory drive predictions, suggesting that the visual environment is key to behavioural isolation of these species. However, the observed plasticity could also make the species barrier vulnerable to environmental change: species‐assortative preferences were weaker in females that were reared in the other species’ light condition.     

Ecosystem size matters: the dimensionality of intralacustrine diversification in Icelandic stickleback is predicted by lake size

Cases of evolutionary diversification can be characterized along a continuum from weak to strong genetic and phenotypic differentiation. Several factors may facilitate or constrain the differentiation process. Comparative analyses of replicates of the same taxon at different stages of differentiation can be useful to identify these factors. We estimated the number of distinct phenotypic groups in three‐spine stickleback populations from nine lakes in Iceland and in one marine population. Using the inferred number of phenotypic groups in each lake, genetic divergence from the marine population, and physical lake and landscape variables, we tested whether ecosystem size, approximated by lake size and depth, or isolation from the ancestral marine gene pool predicts the occurrence and the extent of phenotypic and genetic diversification within lakes. We find intralacustrine phenotypic diversification to be the rule rather than the exception, occurring in all but the youngest lake population and being manifest in ecologically important phenotypic traits. Neutral genetic data further indicate nonrandom mating in four of nine studied lakes, and restricted gene flow between sympatric phenotypic groups in two. Although neither the phenotypic variation nor the number of intralacustrine phenotypic groups was associated with any of our environmental variables, the number of phenotypic traits that were differentiated was significantly positively related to lake size, and evidence for restricted gene flow between sympatric phenotypic groups was only found in the largest lakes where trait specific phenotypic differentiation was highest.     

Taxonomy of the long‐tailed mouse Oligoryzomys destructor (Sigmodontinae: Oryzomyini) with the designation of neotypes for Hesperomys destructor Tschudi, 1844 and Hesperomys melanostoma Tschudi, 1844

The genus Oligoryzomys, distributed from southern South America to southern North America, is the most diverse of the tribe Oryzomyini of sigmodontine rodents. Even when 22 species are currently recognized, species boundaries are unclear for several forms. The species Oligoryzomys destructor is one of the least studied species of the genus and is the one with the largest distribution along the Andes (from southern Colombia to northern Bolivia). The species was described without the selection of a holotype and indication of its type locality. In addition, several taxa are regarded as synonyms of O. destructor. These facts are relevant because previous analysis of DNA sequences has shown that O. destructor represents a species complex. Herein, in addition to test the phylogenetic position of O. destructor within the genus Oligoryzomys, we assess patterns of morphological and molecular variation of O. destructor and its associated nominal forms aimed to assess the boundaries of the species. As part of the study, we selected neotypes for Hesperomys destructor and H. melanostoma. At the light of our results, we recognized O. destructor as a species with two subspecies, O. d. destructor and O. d. spodiurus. Also, we discuss the role of Andean rivers, and their different permeability, as allopatric barriers molding the structure of O. destructor.     

A test for within‐lake niche differentiation in the nine‐spined sticklebacks (Pungitius pungitius)

Specialization for the use of different resources can lead to ecological speciation. Accordingly, there are numerous examples of ecologically specialized pairs of fish “species” in postglacial lakes. Using a polymorphic panel of single nucleotide variants, we tested for genetic footprints of within‐lake population stratification in nine‐spined sticklebacks (Pungitius pungitius) collected from three habitats (viz. littoral, benthic, and pelagic) within a northern Swedish lake. Analyses of admixture, population structure, and relatedness all supported the conclusion that the fish from this lake form a single interbreeding unit.     

The accumulation of reproductive isolation in early stages of divergence supports a role for sexual selection

Models of speciation by sexual selection propose that male–female coevolution leads to the rapid evolution of behavioural reproductive isolation. Here, we compare the strength of behavioural isolation to ecological isolation, gametic incompatibility and hybrid inviability in a group of dichromatic stream fishes. In addition, we examine whether any of these individual barriers, or a combined measure of total isolation, is predicted by body shape differences, male colour differences, environmental differences or genetic distance. Behavioural isolation reaches the highest values of any barrier and is significantly greater than ecological isolation. No individual reproductive barrier is associated with any of the predictor variables. However, marginally significant relationships between male colour and body shape differences with ecological and behavioural isolation are discussed. Differences in male colour and body shape predict total reproductive isolation between species; hierarchical partitioning of these two variables' effects suggests a stronger role for male colour differences. Together, these results suggest an important role for divergent sexual selection in darter speciation but raise new questions about the mechanisms of sexual selection at play and the role of male nuptial ornaments.     

Paternity analyses in wild‐caught and laboratory‐reared Caribbean cricket females reveal the influence of mating environment on post‐copulatory sexual selection

Polyandry is ubiquitous in insects and provides the conditions necessary for male‐ and female‐driven forms of post‐copulatory sexual selection to arise. Populations of Amphiacusta sanctaecrucis exhibit significant divergence in portions of the male genitalia that are inserted directly into the female reproductive tract, suggesting that males may exercise some post‐copulatory control over fertilization success. We examine the potential for male–male and male–female post‐copulatory interactions to influence paternity in wild‐caught females of A. sanctaecrucis and contrast our findings with those obtained from females reared in a high‐density laboratory environment. We find that female A. sanctaecrucis exercise control by mating multiple times (females mount males), but that male–male post‐copulatory interactions may influence paternity success. Moreover, post‐copulatory interactions that affect reproductive success of males are not independent of mating environment: clutches of wild‐caught females exhibit higher sire diversity and lower paternity skew than clutches of laboratory‐reared females. There was no strong evidence for last male precedence in either case. Most attempts at disentangling the contributions of male–male and male–female interactions towards post‐copulatory sexual selection have been undertaken in a laboratory setting and may not capture the full context in which they take place – such as the relationship between premating and post‐mating interactions. Our results reinforce the importance of designing studies that can capture the multifaceted nature of sexual selection for elucidating the role of post‐copulatory sexual selection in driving the evolution of male and female reproductive traits, especially when different components (e.g. precopulatory and post‐copulatory interactions) do not exert independent effects on reproductive outcomes.     

‘Becoming a species by becoming a pest’ or how two maize pests of the genus Ostrinia possibly evolved through parallel ecological speciation events

New agricultural pest species attacking introduced crops may evolve from pre‐existing local herbivores by ecological speciation, thereby becoming a species by becoming a pest. We compare the evolutionary pathways by which two maize pests (the Asian and the European corn borers, ACB and ECB) in the genus Ostrinia (Lepidoptera, Crambidae) probably diverged from an ancestral species close to the current Adzuki bean borer (ABB). We typed larval Ostrinia populations collected on maize and dicotyledons across China and eastern Siberia, at microsatellite and mitochondrial loci. We found only two clusters: one on maize (as expected) and a single one on dicotyledons despite differences in male mid‐tibia morphology, suggesting that all individuals from dicotyledons belonged to the ABB. We found evidence for migrants and hybrids on both host plant types. Hybrids suggest that field reproductive isolation is incomplete between ACB and ABB. Interestingly, a few individuals with an ‘ABB‐like’ microsatellite profile collected on dicotyledons had ‘ACB’ mtDNA rather than ‘ABB‐like’ mtDNA, whereas the reverse was never found on maize. This suggests asymmetrical gene flow directed from the ACB towards the ABB. Hybrids and backcrosses in all directions were obtained in no‐choice tests. In laboratory conditions, they survived as well as parental strain individuals. In Xinjiang, we found ACB and ECB in sympatry, but no hybrids. Altogether, our results suggest that reproductive isolation between ACB and ABB is incomplete and mostly prezygotic. This points to ecological speciation as a possible evolutionary scenario, as previously found for ECB and ABB in Europe.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

A test of genomic modularity among life‐history adaptations promoting speciation with gene flow

Speciation with gene flow may require adaptive divergence of multiple traits to generate strong ecologically based reproductive isolation. Extensive negative pleiotropy or physical linkage of genes in the wrong phase affecting these diverging traits may therefore hinder speciation, while genetic independence or “modularity” among phenotypic traits may reduce constraints and facilitate divergence. Here, we test whether the genetics underlying two components of diapause life history, initial diapause intensity and diapause termination timing, constrain differentiation between sympatric hawthorn and apple‐infesting host races of the fly Rhagoletis pomonella through analysis of 10,256 SNPs measured via genotyping‐by‐sequencing (GBS). Loci genetically associated with diapause termination timing were mainly observed for SNPs mapping to chromosomes 1–3 in the genome, most notably for SNPs displaying higher levels of linkage disequilibrium (LD), likely due to inversions. In contrast, selection on initial diapause intensity affected loci on all five major chromosomes of the genome, specifically those showing low levels of LD. This lack of overlap in genetically associated loci suggests that the two diapause phenotypes are largely modular. On chromosome 2, however, intermediate level LD loci and a subgroup of high LD loci displayed significant negative relationships between initial diapause intensity and diapause termination time. These gene regions on chromosome 2 therefore affected both traits, while most regions were largely independent. Moreover, loci associated with both measured traits also tended to exhibit highly divergent allele frequencies between the host races. Thus, the presence of nonoverlapping genetic modules likely facilitates simultaneous, adaptive divergence for the measured life‐history components.     

A life‐cycle approach to species barriers

What maintains reproductive barriers between closely related species is, of course, of fundamental interest to a closer understanding of the mechanisms that generate new biodiversity. One important dichotomy is to separate barriers evolved from divergent selection over environmental gradients (extrinsic barriers) from barriers caused by incompatibilities between different genetic arrangements that may have evolved in isolation (intrinsic barriers). This dichotomy also reflects an important applied consequence. As the extrinsic barriers are associated with specific environmental contexts, they may be partly or completely erased if the environment changes. In contrast, intrinsic barriers are inert to the environmental context and resistant to environmental changes. From a conservation biology perspective, it may thus be important to be able to separate extrinsic and intrinsic species barriers, but this may in many organisms be a complex matter. In this issue of Molecular Ecology, Montecinos et al. ( 2017 ) found a tractable approach that works for species with life cycles that include two reproductive but ecologically similar generations, one haploid and the other diploid. What they demonstrate is that using a life‐cycle approach offers a unique possibility to separate between prezygotic and postzygotic barriers. Indeed, in the case of an isomorphic life cycle, there is even a possibility to suggest whether postzygotic barriers are more likely to be intrinsic or extrinsic. In this way, their approach may be useful both to increase our understanding of the basic mechanisms of speciation and to single out when species barriers will better resist environmental changes.