Assessing the potential of RAD-sequencing to resolve phylogenetic relationships within species radiations: The fly genus Chiastocheta (Diptera: Anthomyiidae) as a case study

Species delimitation, discovery and conservation in a tiger beetle species complex despite discordant genetic data

In an age of species declines, delineating and discovering biodiversity is critical for both taxonomic accuracy and conservation. In recent years, there has been a movement away from using exclusively morphological characters to delineate and describe taxa and an increase in the use of molecular markers to describe diversity or through integrative taxonomy, which employs traditional morphological characters, as well as genetic or other data. Tiger beetles are charismatic, of conservation concern, and much work has been done on the morphological delineation of species and subspecies, but few of these taxa have been tested with genetic analyses. In this study, we tested morphologically based taxonomic hypotheses of polymorphic tiger beetles in the Eunota circumpicta (LaFerté-Sénectère, 1841) species complex using multilocus genomic and mtDNA analyses. We find multiple cryptic species within the previous taxonomic concept of Eunota circumpicta, some of which were historically recognized as subspecies. We found that the mtDNA and genomic datasets did not identify the same taxonomic units and that the mtDNA was most at odds with all other genetic and morphological patterns. Overall, we describe new cryptic diversity, which raises important conservation concerns, and provide a working example for testing species and subspecies validity despite discordant data.

Integrated phylogenomic approaches in insect systematics

The increased accessibility of genomic and imaging methods, and the improved access to ecological, spatial, and other natural history-related data is allowing for insect systematics to grow and find answers to central evolutionary and taxonomic questions. Today, integrated studies in insect phylogenomics and systematics are combining natural history, behavior, developmental biology, morphology, fossils, geographic range data, and ecological interactions. This integration is contributing to the clarification of evolutionary relationships, and the recognition of the role played by these factors on the evolution of insects. Future work should continue to build on these advances, seeking to further increase open-access databasing and support for natural history research, as well as expand its analytical palettes.

Phylogenetic assessment of the Metamasius hemipterus species complex (Coleoptera, Curculionidae, Dryophthorinae)

Metamasius is a large genus of dryophthorine weevils, with nearly 85 species. Among the economically important pests in the genus, M. hemipterus is currently separated in three subspecies, based largely on color patterns of the elytra, pronotum, and sternum. The tenuous limits of M. hemipterus subspecies were created over fifty years ago and never tested under a phylogenetic framework. Here, for the first time, we address the M. hemipterus species boundaries applying a molecular approach. We constructed a reduced genome representation of a few species using restriction site-associated DNA sequencing (RADseq). Phylogenetic analysis using either a complete supermatrix or only SNPs revealed a clear separation of Metamasius species. We suggest that M. h. carbonarius syn. nov. and M. h. sericeus be treated as the same species, M. sericeus (Oliver) stat. n., and elevate M. h. hemipterus as a separate species M. hemipterus (Linnaeus). We updated Vaurie's identification key to reflect the new species status. This systematic reassessment reflects a more natural classification for these remarkable and economically significant weevils.

DNA Barcoding versus Morphological Variability of Pterostichus brevicornis brevicornis (Kirby, 1837) (Coleoptera, Carabidae) in the Arctic and Subarctic

Simple Summary

Taxonomic studies on a polymorphic species inhabiting a region with relatively uniform environmental conditions (e.g., the Arctic) should involve an integrative approach. Ground beetles such as the subgenus Cryobius of the genus Pterostichus are a successful group in expansion in the tundra biome. The current taxonomy of Cryobius species is unclear and could be considered an obstacle to ecological studies; knowledge of their distribution patterns in the Arctic is rather limited. In this study, the first report on the phylogeography and phylogeny of the most abundant tundra subspecies, P. (Cryobius) brevicornis brevicornis (Kirby, 1837), within its continuous range throughout northern Eurasia and North America is presented. The results indicated that the male genitalia morphology of P. b. brevicornis from Eurasian populations shared a higher geographic variability compared with the pronotum shape and the mitochondrial DNA sequences.

Abstract

The geographic patterns of genetic and morphological variability in ground beetles were examined throughout Northern Eurasia and North America using the most abundant circumpolar tundra subspecies, Pterostichus (Cryobius) brevicornis brevicornis (Kirby, 1837), as a model. Phylogenetic structure was assessed on the basis of a Bayesian approach using two DNA markers (partial sequences of the COI and 28S rRNA genes), while phylogeographic patterns and population genetic diversity were estimated using the COI gene only. Morphological patterns were analysed using elliptical Fourier coefficients that were calculated based on the pronotum and male genitalia shape outlines. The subspecies shares 23 COI haplotypes throughout its entire circumpolar range, while eight haplotypes of 28S rRNA were detected in Northern Eurasia. Phylogenetic analysis did not reveal subdivided species lineages with strict geographical imprint. The network, F_ST and uncorrected pairwise divergence analyses showed that the genetic distances between populations increase by longitude from Northeastern Asia to Europe. The genetic variability among the five studied geographical population groups of P. b. brevicornis was relatively high. The MANOVA showed significant regional divergence between local populations in Northern Eurasia based on both morphological markers, but only male genitalia variability was geographically structured. Neither the pronotum shape nor the male genitalia shape aligned with the phylogeographic patterns discovered on the basis of COI sequences. The genetic (COI) marker had more variation within, rather than among, population groups in addition to morphology of pronotum but not male genitalia.

DiscoSnp-RAD: de novo detection of small variants for RAD-Seq population genomics

Restriction site Associated DNA Sequencing (RAD-Seq) is a technique characterized by the sequencing of specific loci along the genome that is widely employed in the field of evolutionary biology since it allows to exploit variants (mainly Single Nucleotide Polymorphism-SNPs) information from entire populations at a reduced cost. Common RAD dedicated tools, such as STACKS or IPyRAD, are based on all-vs-all read alignments, which require consequent time and computing resources. We present an original method, DiscoSnp-RAD, that avoids this pitfall since variants are detected by exploiting specific parts of the assembly graph built from the reads, hence preventing all-vs-all read alignments. We tested the implementation on simulated datasets of increasing size, up to 1,000 samples, and on real RAD-Seq data from 259 specimens of Chiastocheta flies, morphologically assigned to seven species. All individuals were successfully assigned to their species using both STRUCTURE and Maximum Likelihood phylogenetic reconstruction. Moreover, identified variants succeeded to reveal a within-species genetic structure linked to the geographic distribution. Furthermore, our results show that DiscoSnp-RAD is significantly faster than state-of-the-art tools. The overall results show that DiscoSnp-RAD is suitable to identify variants from RAD-Seq data, it does not require time-consuming parameterization steps and it stands out from other tools due to its completely different principle, making it substantially faster, in particular on large datasets.

Simulation with RADinitio improves RADseq experimental design and sheds light on sources of missing data

Restriction-site associated DNA sequencing (RADseq) has become a powerful and versatile tool in modern population genomics, enabling large-scale evolutionary and genomic analyses in otherwise inaccessible biological systems. With its widespread use, different variants on the protocol have been developed to suit specific experimental needs. Researchers face the challenge of choosing the optimal molecular and sequencing protocols for their reduced representation experimental design, an often-complicated process. Strategic errors can lead to biased data generation that has reduced power to answer biological questions. Here, we present RADinitio, simulation software for the selection and optimization of RADseq experiments via the generation of sequencing data that behave similarly to empirical sources. RADinitio provides an evolutionary simulation of populations, implementation of various RADseq protocols with customizable parameters, and thorough assessment of missing data. We test the efficacy of the software using different RAD protocols across several organisms, highlighting the importance of protocol selection on the magnitude and quality of data acquired. Additionally, we test the effects of RAD library preparation and sequencing on allelic dropout, observing that library preparation and sequencing often contributes more to missing alleles than population-level variation.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Genomics overrules mitochondrial DNA, siding with morphology on a controversial case of species delimitation

Species delimitation is a major quest in biology and is essential for adequate management of the organismal diversity. A challenging example comprises the fish species of red snappers in the Western Atlantic. Red snappers have been traditionally recognized as two separate species based on morphology: Lutjanus campechanus (northern red snapper) and L. purpureus (southern red snapper). Recent genetic studies using mitochondrial markers, however, failed to delineate these nominal species, leading to the current lumping of the northern and southern populations into a single species ( L. campechanus). This decision carries broad implications for conservation and management as red snappers have been commercially over-exploited across the Western Atlantic and are currently listed as vulnerable. To address this conflict, we examine genome-wide data collected throughout the range of the two species. Population genomics, phylogenetic and coalescent analyses favour the existence of two independent evolutionary lineages, a result that confirms the morphology-based delimitation scenario in agreement with conventional taxonomy. Despite finding evidence of introgression in geographically neighbouring populations in northern South America, our genomic analyses strongly support isolation and differentiation of these species, suggesting that the northern and southern red snappers should be treated as distinct taxonomic entities.

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build-up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double-digest restriction site-associated DNA sequencing (ddRADseq) to obtain and analyse genomic-scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.