Phylogenetic insights into Central European Chorthippus and Pseudochorthippus (Orthoptera: Acrididae) species using ddRADseq data

The evolution of several orthopteran groups, especially within the grasshopper family Acrididae, remains poorly understood. This is particularly true for the subfamily Gomphocerinae, which comprises cryptic sympatric and syntopic species. Previous mitochondrial studies have highlighted major discrepancies between taxonomic and phylogenetic hypotheses, thereby emphasizing the necessity of genome-wide approaches. In this study, we employ double-digest restriction site-associated DNA sequencing (ddRADseq) to reconstruct the evolution of Central European Chorthippus and Pseudochorthippus species, especially C.smardai, P.tatrae and the C.biguttulus group. Our phylogenomic analyses recovered deep discordance with mitochondrial DNA barcoding, emphasizing its unreliability in Gomphocerinae grasshoppers. Specifically, our data robustly distinguished the C.biguttulus group and confirmed the distinctiveness of C.eisentrauti, also shedding light on its presence in the Berchtesgaden Alps. Moreover, our results support the reclassification of C.smardai to the genus Pseudochorthippus and of P.tatrae to the genus Chorthippus. Our study demonstrates the efficiency of high-throughput genomic methods such as RADseq without prior optimization to elucidate the complex evolution of grasshopper radiations with direct taxonomic implications. While RADseq has predominantly been utilized for population genomics and within-genus phylogenomics, its application extends to resolve relationships between deeply-diverged clades representative of distinct genera.

The Genomics of Isolated Populations of Gampsocleis glabra (Orthoptera: Tettigoniidae) in Central and Western Europe

Habitat destruction and fragmentation are among the major current threats to global biodiversity. Fragmentation may also affect species with good dispersal abilities. We study the heath bushcricket Gampsocleis glabra, a specialist of steppe-like habitats across Europe that are highly fragmented, investigating if these isolated populations can be distinguished using population genomics and if there are any traces of admixture or dispersal among them. We try to answer these questions using genome-wide SNP data generated with ddRAD sequencing. We calculated F-statistics and visualized differentiation using STRUCTURE plots. While limited by the difficulty of sampling this threatened species, our results show that all populations except one that was represented by a singleton were clearly distinct, with pairwise FST values between 0.010 and 0.181. STRUCTURE indicated limited but visible admixture across most populations and probably also an exchange of individuals between populations of Germany and The Netherlands. We conclude that in G. glabra, a certain amount of gene flow has persisted, at least in the past, also among populations that are isolated today. We also detect a possibly more recent dispersal event between a population in The Netherlands and one in Germany, which may be human aided. We suggest that the conservation of larger populations should be maintained, that efforts should be taken to restore abandoned habitat, that the preservation even of small habitat fragments may be beneficial for the conservation of this species, and that these habitats should be regularly monitored for possible (re-)colonization.

Tandem Repeat DNA Provides Many Cytological Markers for Hybrid Zone Analysis in Two Subspecies of the Grasshopper Chorthippus parallelus

Recent advances in next generation sequencing (NGS) have greatly increased our understanding of non-coding tandem repeat (TR) DNA. Here we show how TR DNA can be useful for the study of hybrid zones (HZ), as it serves as a marker to identify introgression in areas where two biological entities come in contact. We used Illumina libraries to analyse two subspecies of the grasshopper Chorthippus parallelus, which currently form a HZ in the Pyrenees. We retrieved a total of 152 TR sequences, and used fluorescent in situ hybridization (FISH) to map 77 families in purebred individuals from both subspecies. Our analysis revealed 50 TR families that could serve as markers for analysis of this HZ, using FISH. Differential TR bands were unevenly distributed between chromosomes and subspecies. Some of these TR families yielded FISH bands in only one of the subspecies, suggesting the amplification of these TR families after the geographic separation of the subspecies in the Pleistocene. Our cytological analysis of two TR markers along a transect of the Pyrenean hybrid zone showed asymmetrical introgression of one subspecies into the other, consistent with previous findings using other markers. These results demonstrate the reliability of TR-band markers for hybrid zone studies.

Genomic evidence of speciation by fusion in a recent radiation of grasshoppers

Postdivergence gene flow can trigger a number of creative evolutionary outcomes, ranging from the transfer of beneficial alleles across species boundaries (i.e., adaptive introgression) to the formation of new species (i.e., hybrid speciation). Although neutral and adaptive introgression has been broadly documented in nature, hybrid speciation is assumed to be rare and the evolutionary and ecological context facilitating this phenomenon still remains controversial. Through combining genomic and phenotypic data, we evaluate the hypothesis that the dual feeding regime (based on both scrub legumes and gramineous herbs) of the taxonomically controversial grasshopper Chorthippus saulcyi algoaldensis resulted from hybridization between the sister taxa C. binotatus (that exclusively feeds on scrub legumes) and C. saulcyi (that only feeds on gramineous herbs). Genetic clustering analyses and inferences from coalescent-based demographic simulations confirm that C. s. algoaldensis represents an independently evolving lineage and support the ancient hybrid origin of this taxon (about 1.4 Ma), which sheds light on its uncertain phylogenetic position and might explain its broader trophic niche. We propose a Pleistocene hybrid speciation model where range shifts resulting from climatic oscillations can promote the formation of hybrid swarms and facilitate their long-term persistence through geographic isolation from parental forms in topographically complex landscapes.

Extensive introgression at late stages of species formation: Insights from grasshopper hybrid zones

The process of species formation is characterised by the accumulation of multiple reproductive barriers. The evolution of hybrid male sterility, or Haldane's rule, typically characterises later stages of species formation, when reproductive isolation is strongest. Yet, understanding how quickly reproductive barriers evolve and their consequences for maintaining genetic boundaries between emerging species remains a challenging task because it requires studying taxa that hybridise in nature. Here, we address these questions using the meadow grasshopper Pseudochorthippus parallelus, where populations that show multiple reproductive barriers, including hybrid male sterility, hybridise in two natural hybrid zones. Using mitochondrial data, we infer that such populations have diverged some 100,000 years ago, at the beginning of the last glacial cycle in Europe. Nuclear data shows that contractions at multiple glacial refugia, and post-glacial expansions have facilitated genetic differentiation between lineages that today interact in hybrid zones. We find extensive introgression throughout the sampled species range, irrespective of current strength of reproductive isolation. Populations exhibiting hybrid male sterility in two hybrid zones show repeatable patterns of genomic differentiation, consistent with shared genomic constraints affecting ancestral divergence or with the role of those regions in reproductive isolation. Together, our results suggest that reproductive barriers that characterise late stages of species formation can evolve relatively quickly, particularly when associated with strong demographic changes. Moreover, we show that such barriers persist in the face of extensive gene flow, allowing future studies to identify associated genomic regions.

Songs and morphology in three species of the Chorthippusbiguttulus group (Orthoptera, Acrididae, Gomphocerinae) in Russia and adjacent countries

Songs and morphology are compared between Chorthippusmiramae (Vorontsovsky, 1928) that was previously named as C.porphyropterus and two other closely related species, C.brunneus (Thunberg, 1815) and C.maritimus Mistshenko, 1951. We compare them because the calling song of C.miramae was previously shown to have song elements similar to those of other two species. One morphological character, the length of stridulatory file, appeared to be the best character to distinguish between all three species. For C.maritimus and C.miramae, we present the morphological descriptions since they are absent in the literature. We also establish the synonymy C.maritimus = C.bornhalmi Harz, 1971, syn. n. = C.biguttuluseximius Mistshenko, 1951, syn. n. In the song analysis, we analyse not only the sound but also the leg-movement pattern, which is very helpful to find a homology between various song elements. We show that the calling song of C.miramae usually contains two elements, one element being similar to the C.brunneus calling song, and another – to the C.maritimus calling song. Despite some similarities, the calling song elements in C.miramae have some peculiarities. The courtship song of C.miramae is similar to the C.brunneus song, whereas the rivalry songs of C.miramae comprise both the maritimus-like elements and the unique ones. C.miramae generally demonstrates a richer song repertoire than the other two species.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae)

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.