Genealogical discordance and patterns of introgression and selection across a cricket hybrid zone

Phylomitogenomics of two Neotropical species of long-legged crickets Endecous Saussure, 1878 (Orthoptera: Phalangopsidae)

Mitochondrial genomes have provided significant insights into the evolution of several insects. A typical mitogenome contains 37 genes, and variations in gene order can indicate evolutionary relationships between species. In this study, we have assembled the first complete mitogenomes of Endecous chape and E. onthophagus and analyzed the phylogenetic implications for the Gryllidea infraorder. We performed DNA extractions and genome sequencing for both Endecous species. Subsequently, we searched for raw data in the Sequence Read Archive (SRA) in NCBI. Using the SRA data, we assembled the partial mitogenome of Dianemobius nigrofasciatus and annotated the protein-coding genes (PCGs) for nine species. Phylogenomic relationships were reconstructed using Maximum Likelihood (ML) and Bayesian Inference (BI), utilizing the PCGs from 49 Gryllidea species. The mitogenome lengths of E. chape and E. onthophagus are 16,266 bp and 16,023 bp, respectively, while D. nigrofasciatus has a length of 15,359 bp. Our results indicate that species within the infraorder exhibit four types of gene order arrangements that align with their phylogenetic relationships. Both phylogenomic trees displayed strong support, and the ML corroborated with the literature. Gryllidea species have significantly contributed to various fields, and studying their mitogenomes can provide valuable insights into this infraorder evolution.

The evolutionary outcomes of climate-change-induced hybridization in insect populations

Rapid range shifts are one of the most frequent responses to climate change in insect populations. Climate-induced range shifts can lead to the breakdown of isolation barriers, and thus, to an increase in hybridization and introgression. Long-term evolutionary consequences such as the formation of hybrid zones, introgression, speciation, and extinction have been predicted as a result of climate-induced hybridization. Our review shows that there has been an increase in the number of published cases of climate-induced hybridization in insects, and that the formation of hybrid zones and introgression seems to be, at the moment, the most frequent outcomes. Although introgression is considered positive, since it increases species' genetic diversity, in the long term, it could lead to negative outcomes such as species fusion or genetic swamping.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Gene flow, divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.

Mapping reduced introgression loci to the X chromosome of the hybridizing field crickets, Gryllus firmus and G. pennsylvanicus

The genomic architecture of barriers to gene exchange during the speciation process is poorly understood. The genomic islands model suggests that loci associated with barriers to gene exchange prevent introgression of nearby genomic regions via linkage disequilibrium. But few analyses of the actual genomic location of non-introgressing loci in closely related species exist. In a previous study Maroja et al. showed that in the hybridizing field crickets, Gryllus firmus and G. pennsylvanicus, 50 non-introgressing loci are localized on two autosomal regions and the X chromosome, but they were not able to map the loci along the X chromosome because they used a male informative cross. Here, we localize the introgressing and non-introgressing loci on the X chromosome, and reveal that all X-linked non-introgressing loci are restricted to a 50-cM region with 10 of these loci mapped to a single location. We discuss the implications of this finding to speciation.

Speciation genes are more likely to have discordant gene trees

Speciation genes are responsible for reproductive isolation between species. By directly participating in the process of speciation, the genealogies of isolating loci have been thought to more faithfully represent species trees. The unique properties of speciation genes may provide valuable evolutionary insights and help determine the true history of species divergence. Here, we formally analyze whether genealogies from loci participating in Dobzhansky-Muller (DM) incompatibilities are more likely to be concordant with the species tree under incomplete lineage sorting (ILS). Individual loci differ stochastically from the true history of divergence with a predictable frequency due to ILS, and these expectations-combined with the DM model of intrinsic reproductive isolation from epistatic interactions-can be used to examine the probability of concordance at isolating loci. Contrary to existing verbal models, we find that reproductively isolating loci that follow the DM model are often more likely to have discordant gene trees. These results are dependent on the pattern of isolation observed between three species, the time between speciation events, and the time since the last speciation event. Results supporting a higher probability of discordance are found for both derived-derived and derived-ancestral DM pairs, and regardless of whether incompatibilities are allowed or prohibited from segregating in the same population. Our overall results suggest that DM loci are unlikely to be especially useful for reconstructing species relationships, even in the presence of gene flow between incipient species, and may in fact be positively misleading.

Rapid movement and instability of an invasive hybrid swarm

Unstable hybrid swarms that arise following the introduction of non‐native species can overwhelm native congeners, yet the stability of invasive hybrid swarms has not been well documented over time. Here, we examine genetic variation and clinal stability across a recently formed hybrid swarm involving native blacktail shiner (Cyprinella venusta) and non‐native red shiner (C. lutrensis) in the Upper Coosa River basin, which is widely considered to be a global hot spot of aquatic biodiversity. Examination of phenotypic, multilocus genotypic, and mitochondrial haplotype variability between 2005 and 2011 revealed that the proportion of hybrids has increased over time, with more than a third of all sampled individuals exhibiting admixture in the final year of sampling. Comparisons of clines over time indicated that the hybrid swarm has been rapidly progressing upstream, but at a declining and slower pace than rates estimated from historical collection records. Clinal comparisons also showed that the hybrid swarm has been expanding and contracting over time. Additionally, we documented the presence of red shiner and hybrids farther downstream than prior studies have detected, which suggests that congeners in the Coosa River basin, including all remaining populations of the threatened blue shiner (Cyprinella caerulea), are at greater risk than previously thought.

Heterogeneous genome divergence, differential introgression, and the origin and structure of hybrid zones

Hybrid zones have been promoted as windows on the evolutionary process and as laboratories for studying divergence and speciation. Patterns of divergence between hybridizing species can now be characterized on a genomewide scale, and recent genome scans have focused on the presence of 'islands' of divergence. Patterns of heterogeneous genomic divergence may reflect differential introgression following secondary contact and provide insights into which genome regions contribute to local adaptation, hybrid unfitness and positive assortative mating. However, heterogeneous genome divergence can also arise in the absence of any gene flow, as a result of variation in selection and recombination across the genome. We suggest that to understand hybrid zone origins and dynamics, it is essential to distinguish between genome regions that are divergent between pure parental populations and regions that show restricted introgression where these populations interact in hybrid zones. The latter, more so than the former, reveal the likely genetic architecture of reproductive isolation. Mosaic hybrid zones, because of their complex structure and multiple contacts, are particularly good subjects for distinguishing primary intergradation from secondary contact. Comparisons among independent hybrid zones or transects that involve the 'same' species pair can also help to distinguish between divergence with gene flow and secondary contact. However, data from replicate hybrid zones or replicate transects do not reveal consistent patterns; in a few cases, patterns of introgression are similar across independent transects, but for many taxa, there is distinct lack of concordance, presumably due to variation in environmental context and/or variation in the genetics of the interacting populations.

Differential introgression suggests candidate beneficial and barrier loci between two parapatric subspecies of Pearson's horseshoe bat Rhinolophus pearsoni

Observations that rates of introgression between taxa can vary across loci are increasingly common. Here, we test for differential locus-wise introgression in 2 parapatric subspecies of Pearson′s horseshoe bat (Rhinolophus pearsoni chinensis and R. p. pearsoni). To efficiently identify putative speciation genes and/or beneficial genes in our current system, we used a candidate gene approach by including loci from X chromosome that are suggested to be more likely involved in reproductive isolation in other organisms and loci underlying hearing that have been suggested to spread across the hybrid zone in another congeneric species. Phylogenetic and coalescent analyses were performed at 2 X-linked, 4 hearing genes, as well as 2 other autosomal loci individually. Likelihood ratio tests could not reject the model of zero gene flow at 2 X-linked and 2 autosomal genes. In contrast, gene flow was supported at 3 of 4 hearing genes. While this introgression could be adaptive, we cannot rule out stochastic processes. Our results highlight the utility of the candidate gene approach in searching for speciation genes and/or beneficial genes across the species boundary in natural populations.

Sorted gene genealogies and species-specific nonsynonymous substitutions point to putative postmating prezygotic isolation genes in Allonemobius crickets

In the Allonemobius socius complex of crickets, reproductive isolation is primarily accomplished via postmating prezygotic barriers. We tested seven protein-coding genes expressed in the male ejaculate for patterns of evolution consistent with a putative role as postmating prezygotic isolation genes. Our recently diverged species generally lacked sequence variation. As a result, ω-based tests were only mildly successful. Some of our genes showed evidence of elevated ω values on the internal branches of gene trees. In a couple of genes, these internal branches coincided with both species branching events of the species tree, between A. fasciatus and the other two species, and between A. socius and A. sp. nov. Tex. In comparison, more successful approaches were those that took advantage of the varying degrees of lineage sorting and allele sharing among our young species. These approaches were particularly powerful within the contact zone. Among the genes we tested we found genes with genealogies that indicated relatively advanced degrees of lineage sorting across both allopatric and contact zone alleles. Within a contact zone between two members of the species complex, only a subset of genes maintained allelic segregation despite evidence of ongoing gene flow in other genes. The overlap in these analyses was arginine kinase (AK) and apolipoprotein A-1 binding protein (APBP). These genes represent two of the first examples of sperm maturation, capacitation, and motility proteins with fixed non-synonymous substitutions between species-specific alleles that may lead to postmating prezygotic isolation. Both genes express ejaculate proteins transferred to females during copulation and were previously identified through comparative proteomics. We discuss the potential function of these genes in the context of the specific postmating prezygotic isolation phenotype among our species, namely conspecific sperm precedence and the superior ability of conspecific males to induce oviposition in females.

Population Genetic Structure of the Malaria Vector Anopheles sinensis (Diptera: Culicidae) Sensu Stricto and Evidence for Possible Introgression in the Republic of Korea

Anopheles sinensis Wiedemann sensu stricto (s.s.) is a dominant mosquito and considered a secondary malaria vector in the Republic of Korea (ROK). Despite the potential significance for malaria control, population genetics studies have been conducted using only mitochondrial DNA (mtDNA), and studies of the genetics of hybridization have never been attempted. In this study, 346 specimens from 23 localities were subject to experiments. Among them, 305 An. sinensis s.s. specimens from 20 localities were used for mtDNA analysis, and 346 specimens comprising 341 An. sinensis s.s. from 22 localities and five Anopheles kleini Rueda from one locality were examined in the microsatellite study. Neighbor-joining analysis of pairwise FST and RST based on microsatellite results showed that the populations are divided into two groups, as did the mtDNA results. However, the Bayesian analysis and factorial correspondence analysis plots showed three distinct clusters. Among the mtDNA and microsatellite results, only microsatellites represented small but positive and significant isolation-by-distance patterns. Both molecular markers show the Taebaek and Sobaek Mountain ranges as barriers between the northern and southern parts of the ROK. The newly recognized third group suggests possible introgressive hybridization of An. sinensis s.s. with closely related species. The slightly different composition of populations in each group based on different markers is probably because of different population dynamics in each group. These results imply that there is restricted gene flow of epidemiologically important malaria-related genes between the northern and southern parts of the ROK.

Genes with Restricted Introgression in a Field Cricket (Gryllus firmus/Gryllus pennsylvanicus) Hybrid Zone Are Concentrated on the X Chromosome and a Single Autosome

Characterizing the extent of genomic differentiation between recently diverged lineages provides an important context for understanding the early stages of speciation. When such lineages form discrete hybrid zones, patterns of differential introgression allow direct estimates of which genome regions are likely involved in speciation and local adaptation. Here we use a backcross experimental design to construct a genetic linkage map for the field crickets Gryllus firmus and Gryllus pennsylvanicus, which interact in a well-characterized hybrid zone in eastern North America. We demonstrate that loci with major allele frequency differences between allopatric populations are not randomly distributed across the genome. Instead, most are either X-linked or map to a few small autosomal regions. Furthermore, the subset of those highly differentiated markers that exhibit restricted introgression across the cricket hybrid zone are also concentrated on the X chromosome (39 of 50 loci) and in a single 7-cM region of one autosome. Although the accumulation on the sex chromosome of genes responsible for postzygotic barriers is a well-known phenomenon, less attention has been given to the genomic distribution of genes responsible for prezygotic barriers. We discuss the implications of our results for speciation, both in the context of the role of sex chromosomes and also with respect to the likely causes of heterogeneous genomic divergence. Although we do not yet have direct evidence for the accumulation of ecological, behavioral, or fertilization prezygotic barrier genes on the X chromosome, faster-X evolution could make these barriers more likely to be X-linked.

Hybridization and back-crossing in giant petrels (Macronectes giganteus and M. halli) at Bird Island, South Georgia, and a summary of hybridization in seabirds

Hybridization in natural populations provides an opportunity to study the evolutionary processes that shape divergence and genetic isolation of species. The emergence of pre-mating barriers is often the precursor to complete reproductive isolation. However, in recently diverged species, pre-mating barriers may be incomplete, leading to hybridization between seemingly distinct taxa. Here we report results of a long-term study at Bird Island, South Georgia, of the extent of hybridization, mate fidelity, timing of breeding and breeding success in mixed and conspecific pairs of the sibling species, Macronectes halli (northern giant petrel) and M. giganteus (southern giant petrel). The proportion of mixed-species pairs varied annually from 0.4-2.4% (mean of 1.5%), and showed no linear trend with time. Mean laying date in mixed-species pairs tended to be later than in northern giant petrel, and always earlier than in southern giant petrel pairs, and their breeding success (15.6%) was lower than that of conspecific pairs. By comparison, mixed-species pairs at both Marion and Macquarie islands always failed before hatching. Histories of birds in mixed-species pairs at Bird Island were variable; some bred previously or subsequently with a conspecific partner, others subsequently with a different allospecific partner, and some mixed-species pairs remained together for multiple seasons. We also report the first verified back-crossing of a hybrid giant petrel with a female northern giant petrel. We discuss the potential causes and evolutionary consequences of hybridization and back-crossing in giant petrels and summarize the incidence of back-crossing in other seabird species.

Models of selection, isolation, and gene flow in speciation

Many marine ecologists aspire to use genetic data to understand how selection and demographic history shape the evolution of diverging populations as they become reproductively isolated species. I propose combining two types of genetic analysis focused on this key early stage of the speciation process to identify the selective agents directly responsible for population divergence. Isolation-with-migration (IM) models can be used to characterize reproductive isolation between populations (low gene flow), while codon models can be used to characterize selection for population differences at the molecular level (especially positive selection for high rates of amino acid substitution). Accessible transcriptome sequencing methods can generate the large quantities of data needed for both types of analysis. I highlight recent examples (including our work on fertilization genes in sea stars) in which this confluence of interest, models, and data has led to taxonomically broad advances in understanding marine speciation at the molecular level. I also highlight new models that incorporate both demography and selection: simulations based on these theoretical advances suggest that polymorphisms shared among individuals (a key source of information in IM models) may lead to false-positive evidence of selection (in codon models), especially during the early stages of population divergence and speciation that are most in need of study. The false-positive problem may be resolved through a combination of model improvements plus experiments that document the phenotypic and fitness effects of specific polymorphisms for which codon models and IM models indicate selection and reproductive isolation (such as genes that mediate sperm-egg compatibility at fertilization).

Differential introgression among loci across a hybrid zone of the intermediate horseshoe bat (Rhinolophus affinis)

Background

Hybrid zones formed by the secondary contact of divergent lineages represent natural laboratories for studying the genetic basis of speciation. Here we tested for patterns of differential introgression among three X-linked and 11 autosomal regions to identify candidate loci related to either reproductive isolation or adaptive introgression across a hybrid zone between two Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis: R. a. himalayanus and R. a. macrurus.

Results

Our results support the previous suggestion that macrurus formed when a third subspecies (R. a. hainanus) recolonized the mainland from Hainan Island, and that himalayanus is the ancestral taxon. However, this overall evolutionary history was not reflected in all loci examined, with considerable locus-wise heterogeneity seen in gene tree topologies, levels of polymorphism, genetic differentiation and rates of introgression. Coalescent simulations suggested levels of lineage mixing seen at some nuclear loci might result from incomplete lineage sorting. Isolation with migration models supported evidence of gene flow across the hybrid zone at one intronic marker of the hearing gene Prestin.

Conclusions

We suggest that phylogenetic discordance with respect to the species tree seen here is likely to arise via a combination of incomplete lineage sorting and a low incidence of introgression although we cannot rule out other explanations such as selection and recombination. Two X-linked loci and one autosomal locus were identified as candidate regions related to reproductive isolation across the hybrid zone. Our work highlights the importance of including multiple genomic regions in characterizing patterns of divergence and gene flow across a hybrid zone.

Selection and demographic history shape the molecular evolution of the gamete compatibility protein bindin in Pisaster sea stars

Reproductive compatibility proteins have been shown to evolve rapidly under positive selection leading to reproductive isolation, despite the potential homogenizing effects of gene flow. This process has been implicated in both primary divergence among conspecific populations and reinforcement during secondary contact; however, these two selective regimes can be difficult to discriminate from each other. Here, we describe the gene that encodes the gamete compatibility protein bindin for three sea star species in the genus Pisaster. First, we compare the full-length bindin-coding sequence among all three species and analyze the evolutionary relationships between the repetitive domains of the variable second bindin exon. The comparison suggests that concerted evolution of repetitive domains has an effect on bindin divergence among species and bindin variation within species. Second, we characterize population variation in the second bindin exon of two species: We show that positive selection acts on bindin variation in Pisaster ochraceus but not in Pisaster brevispinus, which is consistent with higher polyspermy risk in P. ochraceus. Third, we show that there is no significant genetic differentiation among populations and no apparent effect of sympatry with congeners that would suggest selection based on reinforcement. Fourth, we combine bindin and cytochrome c oxidase 1 data in isolation-with-migration models to estimate gene flow parameter values and explore the historical demographic context of our positive selection results. Our findings suggest that positive selection on bindin divergence among P. ochraceus alleles can be accounted for in part by relatively recent northward population expansions that may be coupled with the potential homogenizing effects of concerted evolution.

Barriers to gene exchange in hybridizing field crickets: the role of male courtship effort and cuticular hydrocarbons

Background

Pre-zygotic barriers often involve some form of sexual selection, usually interpreted as female choice, as females are typically the choosier sex. However, males typically show some mate preferences, which are increasingly reported. Here we document previously uncharacterized male courtship behavior (effort and song) and cuticular hydrocarbon (CHC) profiles in the hybridizing crickets Gryllus firmus and G. pennsylvanicus. These two species exhibit multiple barriers to gene exchange that act throughout their life history, including a behavioral barrier that results in increased time to mate in heterospecific pairs.

Results

We demonstrated that male mate choice (as courtship effort allocation) plays a more important role in the prezygotic behavioral barrier than previously recognized. In gryllids females ultimately decide whether or not to mate, yet we found males were selective by regulating courtship effort intensity toward the preferred (conspecific) females. Females were also selective by mating with more intensely courting males, which happened to be conspecifics. We report no differences in courtship song between the two species and suggest that the mechanism that allows males to act differentially towards conspecific and heterospecific females is the cuticular hydrocarbon (CHC) composition. CHC profiles differed between males and females of both species, and there were clear differences in CHC composition between female G. firmus and G. pennsylvanicus but not between the males of each species.

Conclusion

Although many barriers to gene exchange are known in this system, the mechanism behind the mate recognition leading to reduced heterospecific mating remains unknown. The CHC profiles might be the phenotypic cue that allow males to identify conspecifics and thus to adjust their courtship intensity accordingly, leading to differential mating between species.

Climate-Driven Reshuffling of Species and Genes: Potential Conservation Roles for Species Translocations and Recombinant Hybrid Genotypes

Comprising 50%-75% of the world's fauna, insects are a prominent part of biodiversity in communities and ecosystems globally. Biodiversity across all levels of biological classifications is fundamentally based on genetic diversity. However, the integration of genomics and phylogenetics into conservation management may not be as rapid as climate change. The genetics of hybrid introgression as a source of novel variation for ecological divergence and evolutionary speciation (and resilience) may generate adaptive potential and diversity fast enough to respond to locally-altered environmental conditions. Major plant and herbivore hybrid zones with associated communities deserve conservation consideration. This review addresses functional genetics across multi-trophic-level interactions including "invasive species" in various ecosystems as they may become disrupted in different ways by rapid climate change. "Invasive genes" (into new species and populations) need to be recognized for their positive creative potential and addressed in conservation programs. "Genetic rescue" via hybrid translocations may provide needed adaptive flexibility for rapid adaptation to environmental change. While concerns persist for some conservationists, this review emphasizes the positive aspects of hybrids and hybridization. Specific implications of natural genetic introgression are addressed with a few examples from butterflies, including transgressive phenotypes and climate-driven homoploid recombinant hybrid speciation. Some specific examples illustrate these points using the swallowtail butterflies (Papilionidae) with their long-term historical data base (phylogeographical diversity changes) and recent (3-decade) climate-driven temporal and genetic divergence in recombinant homoploid hybrids and relatively recent hybrid speciation of Papilio appalachiensis in North America. Climate-induced "reshuffling" (recombinations) of species composition, genotypes, and genomes may become increasingly ecologically and evolutionarily predictable, but future conservation management programs are more likely to remain constrained by human behavior than by lack of academic knowledge.

Detecting adaptive trait loci in nonmodel systems: divergence or admixture mapping?

Mapping adaptive trait loci (ATL) underlying ecological divergence is an essential step towards understanding the processes that generate phenotypic diversity. Technological advances have made it possible to sequence exomes in nonmodel systems, providing an efficient means of analysing functional genetic variants. Divergence scans of genetic markers for outlier loci, or 'divergence mapping', have been used to map locally adapted genes, but this approach is likely to be underpowered when background divergence is elevated. Genotype-phenotype association tests in admixed populations, or 'admixture mapping', may provide a useful approach for mapping locally adapted loci when neutral divergence is high. To determine the power and limits of divergence mapping, we simulated exomes containing a single ATL across two parental populations of varying neutral divergence, estimated divergence and quantified the power to identify the ATL. We found that divergence mapping had very high power when background FST is <0.2, but decreased dramatically above this level. To evaluate the utility of admixture mapping, we simulated exomes from admixed populations, then simulated phenotypes, conducted genotype-phenotype association tests and found that even two generations of random mating after admixture could provide high mapping power in scenarios where pure divergence mapping was ineffective (FST = 0.35). Moreover, admixture mapping had high power across all levels of divergence after 20 generations since admixture. Together with high-throughput exome sequencing, admixture mapping could be used to map ATL in systems such as Heliconius butterflies or Gryllus crickets when experimental design and analytical approach are chosen accordingly.

Population genetics and a study of speciation using next-generation sequencing: an educational primer for use with "Patterns of transcriptome divergence in the male accessory gland of two closely related species of field crickets"

Understanding evidence for the genetic basis of reproductive isolation is imperative for supporting students' understanding of mechanisms of speciation in courses such as Genetics and Evolutionary Biology. An article by Andrés et al. in the February 2013 issue of GENETICS illustrates how advances in DNA sequencing are accelerating studies of population genetics in species with limited genetic and genomic resources. Andrés et al. use the latest sequencing technologies to systematically identify and characterize sites in the DNA that vary within, and have diverged between, species to explore speciation in crickets. This primer, coupled with that article, will help instructors introduce and reinforce important concepts in genetics and evolution while simultaneously introducing modern methodology in the undergraduate classroom. Related article in

GENETICS

Andrés, J. A., E. L. Larson, S. M. Bogdanowicz, and R. G. Harrison, 2013 Patterns of transcriptome divergence in the male accessory gland of two closely related species of field crickets. Genetics 193: 501-513.

Differential introgression in a mosaic hybrid zone reveals candidate barrier genes

Hybrid zones act as genomic sieves. Although globally advantageous alleles will spread throughout the zone and neutral alleles can be freely exchanged between species, introgression will be restricted for genes that contribute to reproductive barriers or local adaptation. Seminal fluid proteins (SFPs) are known to contribute to reproductive barriers in insects and have been proposed as candidate barrier genes in the hybridizing field crickets Gryllus pennsylvanicus and Gryllus firmus. Here, we have used 125 single nucleotide polymorphisms to characterize patterns of differential introgression and to identify genes that may contribute to prezygotic barriers between these species. Using a transcriptome scan of the male cricket accessory gland (the site of SFP synthesis), we identified genes with major allele frequency differences between the species. We then compared patterns of introgression for genes encoding SFPs with patterns for genes expressed in the same tissue that do not encode SFPs. We find no evidence that SFPs have reduced gene exchange across the cricket hybrid zone. However, a number of genes exhibit dramatically reduced introgression, and many of these genes encode proteins with functional roles consistent with known barriers.

Structure of a mosaic hybrid zone between the field crickets Gryllus firmus and G. pennsylvanicus

Hybrid zones provide insight into the nature of species boundaries and the evolution of barriers to gene exchange. Characterizing multiple regions within hybrid zones is essential for understanding both their history and current dynamics. Here, we describe a previously uncharacterized region of a well-studied hybrid zone between two species of field crickets, Gryllus pennsylvanicus and G. firmus. We use a combination of mitochondrial DNA sequencing, morphological data, and modeling of environmental variables to identify the ecological factors structuring the hybrid zone and define patterns of hybridization and introgression. We find an association between species distribution and natural habitat; Gryllus pennsylvanicus occupies natural habitat along forest edges and natural clearings, whereas G. firmus occupies more disturbed areas in agricultural and suburban environments. Hybridization and introgression occur across patch boundaries; there is evidence of substantial admixture both in morphological characters and mtDNA, over a broad geographic area. Nonetheless, the distribution of morphological types is bimodal. Given that F1 hybrids are viable and fertile in the lab, this suggests that strong pre-zygotic barriers are operating in this portion of the hybrid zone.

Patterns of transcriptome divergence in the male accessory gland of two closely related species of field crickets

One of the central questions in evolutionary genetics is how much of the genome is involved in the early stages of divergence between populations, causing them to be reproductively isolated. In this article, we investigate genomic differentiation in a pair of closely related field crickets (Gryllus firmus and G. pennsylvanicus). These two species are the result of allopatric divergence and now interact along an extensive hybrid zone in eastern North America. Genes encoding seminal fluid proteins (SFPs) are often divergent between species, and it has been hypothesized that these proteins may play a key role in the origin and maintenance of reproductive isolation between diverging lineages. Hence, we chose to scan the accessory gland transcriptome to enable direct comparisons of differentiation for genes known to encode SFPs with differentiation in a much larger set of genes expressed in the same tissue. We have characterized differences in allele frequency between two populations for >6000 SNPs and >26,000 contigs. About 10% of all SNPs showed nearly fixed differences between the two species. Genes encoding SFPs did not have significantly elevated numbers of fixed SNPs per contig, nor did they seem to show larger differences than expected in their average allele frequencies. The distribution of allele frequency differences across the transcriptome is distinctly bimodal, but the relatively high proportion of fixed SNPs does not necessarily imply "ancient" divergence between these two lineages. Further studies of linkage disequilibrium and introgression across the hybrid zone are needed to direct our attention to those genome regions that are important for reproductive isolation.

The language of speciation

The literature on speciation has expanded dramatically in recent years, catalyzed by the emergence of new conceptual frameworks, new theoretical approaches, and new methods for characterizing pattern and inferring process. As a consequence, the language used to describe the speciation process has become more complex. Increasing complexity may be an accurate reflection of current thinking with respect to how phenotypic differences limit gene flow, how selection results in the evolution of reproductive isolation, and genetic changes that contribute to speciation. However, increased language complexity has come at a cost; old definitions have been reconfigured and new terms have been introduced. In some instances, the introduction of new terminology has failed to recognize historical usage, leading to unnecessary ambiguity and redundancy. Although the writings of Mayr and Dobzhansky remain a reference point in the language of speciation, the last decades of the 20th century saw substantial changes in our thinking about the speciation process. During that period, the language of speciation remained relatively stable. In contrast, the first decade of the 21st century has witnessed a remarkable expansion of the language of speciation. Here, the origin and evolution of ideas about speciation are viewed through the lens of changing language use.

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito-nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent-based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex-biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.

Discordant introgression in a rapidly expanding hybrid swarm

The erosion of species boundaries can involve rapid evolutionary change. Consequently, many aspects of the process remain poorly understood, including the formation, expansion, and evolution of hybrid swarms. Biological invasions involving hybridization present exceptional opportunities to study the erosion of species boundaries because timelines of interactions and outcomes are frequently well known. Here, we examined clinal variation across codominant and maternally inherited genetic markers as well as phenotypic traits to characterize the expansion and evolution of a hybrid swarm between native Cyprinella venusta and invasive Cyprinella lutrensis minnows. Discordant introgression of phenotype, microsatellite multilocus genotype, and mtDNA haplotype indicates that the observable expansion of the C. venusta × C. lutrensis hybrid swarm is a false invasion front. Both parental and hybrid individuals closely resembling C. lutrensis are numerically dominant in the expansion wake, indicating that the non-native parental phenotype may be selectively favored. These findings show that cryptic introgression can extend beyond the phenotypic boundaries of hybrid swarms and that hybrid swarms likely expand more rapidly than can be documented from phenotypic variation alone. Similarly, dominance of a single parental phenotype following an introduction event may lead to instances of species erosion being mistaken for species displacement without hybridization.

The Tegula tango: a coevolutionary dance of interacting, positively selected sperm and egg proteins

Reproductive proteins commonly show signs of rapid divergence driven by positive selection. The mechanisms driving these changes have remained ambiguous in part because interacting male and female proteins have rarely been examined. We isolate an egg protein the vitelline envelope receptor for lysin (VERL) from Tegula, a genus of free-spawning marine snails. Like VERL from abalone, Tegula VERL is a major component of the VE surrounding the egg, includes a conserved zona pellucida (ZP) domain at its C-terminus, and possesses a unique, negatively charged domain of about 150 amino acids implicated in interactions with the positively charged lysin. Unlike for abalone VERL, where this unique VERL domain occurs in a tandem array of 22 repeats, Tegula VERL has just one such domain. Interspecific comparisons show that both lysin and the VERL domain diverge via positive selection, whereas the ZP domain evolves neutrally. Rates of nonsynonymous substitution are correlated between lysin and the VERL domain, consistent with sexual antagonism, although lineage-specific effects, perhaps owing to different ecologies, may alter the relative evolutionary rates of sperm- and egg-borne proteins.

Post-mating prezygotic barriers to gene exchange between hybridizing field crickets

Studies of sexual selection in speciation have traditionally focused on mate preference, with less attention given to traits that act between copulation and fertilization. However, recent work suggests that post-mating prezygotic barriers may play an important role in speciation. Here, we evaluate the role of such barriers in the field crickets, Gryllus firmus and Gryllus pennsylvanicus. Gryllus pennsylvanicus females mated with G. firmus males produce viable, fertile offspring, but when housed with both species produce offspring sired primarily by conspecifics. We evaluate patterns of sperm utilization in doubly mated G. pennsylvanicus females and find no evidence for conspecific sperm precedence. The reciprocal cross (G. firmus female × G. pennsylvanicus male) produces no progeny. Absence of progeny reflects a barrier to fertilization rather than reduced sperm transfer, storage or motility. We propose a classification scheme for mechanisms underlying post-mating prezygotic barriers similar to that used for premating barriers.

Higher differentiation among subspecies of the house mouse (Mus musculus) in genomic regions with low recombination

In the early stages of reproductive isolation, genomic regions of reduced recombination are expected to show greater levels of differentiation, either because gene flow between species is reduced in these regions or because the effects of selection at linked sites within species are enhanced in these regions. Here, we study the patterns of DNA sequence variation at 27 autosomal loci among populations of Mus musculus musculus, M. m. domesticus, and M. m. castaneus, three subspecies of house mice with collinear genomes. We found that some loci exhibit considerable shared variation among subspecies, while others exhibit fixed differences. We used an isolation-with-gene-flow model to estimate divergence times and effective population sizes (N(e) ) and to disentangle ancestral variation from gene flow. Estimates of divergence time indicate that all three subspecies diverged from one another within a very short period of time approximately 350,000 years ago. Overall, N(e) for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest N(e) and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest N(e) and the smallest proportion of monophyletic gene genealogies. M. m. domesticus and M. m. musculus were more differentiated from each other than either were from M. m. castaneus, consistent with greater reproductive isolation between M. m. domesticus and M. m. musculus. F(ST) was significantly greater at loci experiencing low recombination rates compared to loci experiencing high recombination rates in comparisons between M. m. castaneus and M. m. musculus or M. m. domesticus. These results provide evidence that genomic regions with less recombination show greater differentiation, even in the absence of chromosomal rearrangements.

Phylogeny estimation of the radiation of western North American chipmunks (Tamias) in the face of introgression using reproductive protein genes

The causes and consequences of rapid radiations are major unresolved issues in evolutionary biology. This is in part because phylogeny estimation is confounded by processes such as stochastic lineage sorting and hybridization. Because these processes are expected to be heterogeneous across the genome, comparison among marker classes may provide a means of disentangling these elements. Here we use introns from nuclear-encoded reproductive protein genes expected to be resistant to introgression to estimate the phylogeny of the western chipmunks (Tamias: subgenus: Neotamias), a rapid radiation that has experienced introgressive hybridization of mitochondrial DNA (mtDNA). We analyze the nuclear loci using coalescent-based species-tree estimation methods and concatenation to estimate a species tree and we use parametric bootstraps and coalescent simulations to differentiate between phylogenetic error, coalescent stochasticity and introgressive hybridization. Results indicate that the mtDNA gene tree reflects several introgression events that have occurred between taxa of varying levels of divergence and at different time points in the tree. T. panamintinus and T. speciosus appear to be fixed for ancient mitochondrial introgressions from T. minimus. A southern Rocky Mountains clade appears well sorted (i.e., species are largely monophyletic) at multiple nuclear loci, while five of six taxa are nonmonophyletic based on cytochrome b. Our simulations reject phylogenetic error and coalescent stochasticity as causes. The results represent an advance in our understanding of the processes at work during the radiation of Tamias and suggest that sampling reproductive-protein genes may be a viable strategy for phylogeny estimation of rapid radiations in which reproductive isolation is incomplete. However, a genome-scale survey that can statistically compare heterogeneity of genealogical process at many more loci will be necessary to test this conclusion.

Speciation in the European rabbit (Oryctolagus cuniculus): islands of differentiation on the X chromosome and autosomes

Studies of gene flow between closely related taxa can provide insight into the genetic basis of speciation. To evaluate the importance of the X chromosome in reproductive isolation between subspecies of the European rabbit and to study the genomic scale over which islands of differentiation extend, we resequenced a total of 34 loci distributed along the X chromosome and chromosome 14. Previous studies based on few markers suggested that loci in centromeric regions were highly differentiated between rabbit subspecies, whereas loci in telomeric regions were less differentiated. Here, we confirmed this finding but also discovered remarkable variation in levels of differentiation among loci, with F(ST) values from nearly 0 to 1. Analyses using isolation-with-migration models suggest that this range appears to be largely explained by differential levels of gene flow among loci. The X chromosome was significantly more differentiated than the autosomes. On chromosome 14, differentiation decayed very rapidly at increasing distances from the centromere, but on the X chromosome distinct islands of differentiation encompassing several megabases were observed both at the centromeric region and along the chromosome arms. These findings support the idea that the X chromosome plays an important role in reproductive isolation between rabbit subspecies. These results also demonstrate the mosaic nature of the genome at species boundaries.

Evidence for positive selection in the gene fruitless in Anastrepha fruit flies

Background

Many genes involved in the sex determining cascade have indicated signals of positive selection and rapid evolution across different species. Even though fruitless is an important gene involved mostly in several aspects of male courtship behavior, the few studies so far have explained its high rates of evolution by relaxed selective constraints. This would indicate that a large portion of this gene has evolved neutrally, contrary to what has been observed for other genes in the sex cascade.

Results

Here we test whether the fruitless gene has evolved neutrally or under positive selection in species of Anastrepha (Tephritidae: Diptera) using two different approaches, a long-term evolutionary analysis and a populational genetic data analysis. The first analysis was performed by using sequences of three species of Anastrepha and sequences from several species of Drosophila using the ratio of nonsynonymous to synonymous rates of evolution in PAML, which revealed that the fru region here studied has evolved by positive selection. Using Bayes Empirical Bayes we estimated that 16 sites located in the connecting region of the fruitless gene were evolving under positive selection. We also investigated for signs of this positive selection using populational data from 50 specimens from three species of Anastrepha from different localities in Brazil. The use of standard tests of selection and a new test that compares patterns of differential survival between synonymous and nonsynonymous in evolutionary time also provide evidence of positive selection across species and of a selective sweep for one of the species investigated.

Conclusions

Our data indicate that the high diversification of fru connecting region in Anastrepha flies is due at least in part to positive selection, not merely as a consequence of relaxed selective constraint. These conclusions are based not only on the comparison of distantly related taxa that show long-term divergence time, but also on recently diverged lineages and suggest that episodes of adaptive evolution in fru may be related to sexual selection and/or conflict related to its involvement in male courtship behavior.

The divergence of two independent lineages of an endemic Chinese gecko, Gekko swinhonis, launched by the Qinling orogenic belt

The genetic structure and demographic history of an endemic Chinese gecko, Gekko swinhonis, were investigated by analysing the mitochondrial cytochrome b gene and 10 microsatellite loci for samples collected from 27 localities. Mitochondrial DNA data provided a detailed distribution of two highly divergent evolutionary lineages, between which the average pairwise distance achieved was 0.14. The geographic division of the two lineages coincided with a plate boundary consisting of the Qinling and Taihang Mts, suggesting a historical vicariant pattern. The orogeny of the Qinling Mts, a dispersal and major climatic barrier of the region, may have launched the independent lineage divergence. Both lineages have experienced recent expansion, and the current sympatric localities comprised the region of contact between the lineages. Individual-based phylogenetic analyses of nucDNA and Bayesian-clustering approaches revealed a deep genetic structure analogous to mtDNA. Incongruence between nucDNA and mtDNA at the individual level at localities outside of the contact region can be explained by the different inheritance patterns and male-biased dispersal in this species. High genetic divergence, long-term isolation and ecological adaptation, as well as the morphological differences, suggest the presence of a cryptic species.

Population genomics of secondary contact

One common form of reticulate evolution arises as a consequence of secondary contact between previously allopatric populations. Using extensive coalescent simulations, we describe the conditions for, and extent of, the introgression of genetic material into the genome of a colonizing population from an endemic population. The simulated coalescent histories are sampled from models that describe the evolution of entire chromosomes, thereby allowing the expected length of introgressed haplotypes to be estimated. The results indicate that our ability to identify reticulate evolution from genetic data is highly variable and depends critically upon the duration of the period of allopatry, the timing of the secondary contact event, as well as the sizes of the populations at the time of contact. One particularly interesting result arises when secondary contact occurs close to the time of a severe founder event, in this case, genetic introgression can be substantially more difficult to detect. However, if secondary contact occurs after such a founding event, when the range of the colonizing population increases, introgression is more readily detectable across the genome. This result may have important implications for our ability to detect introgression between ancestrally bottlenecked modern human populations and archaic hominin species, such as Neanderthals.