Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps

Baetidae (Insecta: Ephemeroptera) of Aurès Mountains (Algeria): A New Species of the Baetis alpinus Species Group, with Notes on Baetis Laech, 1815 Biogeography within Maghreb

A new species, Baetis (Baetis) dihyaesp. nov., belonging to the Baetis alpinus species group, is described and illustrated based on larval material collected in the Aurès Mountains (northeastern Algeria) in 2020-2021. This new species is closely related to three European species, e.g., Baetis (B.) alpinus (Pictet, 1843); B. (B.) nubecularis Eaton, 1898; and B. (B.) pasquetorum Righetti & Thomas, 2002 by the combination of the following characteristics: (i) more than one short, stout bristle at the tip of segment II of the maxillary palp and (ii) a well-developed paracercus. However, the new species clearly differs from all congeners of the Baetis alpinus species group primarily by the (a) structure of mouthparts-with 14-18 long submarginal setae arranged in a single irregular row on the dorsal surface of the labrum; 2-6 short, stout bristles at the tip of segment II of the maxillary palp; and segment II of the labial palp without a considerably developed apico-internal lobe); (b) setation of abdominal terga, with a few triangular-shaped scales sparsely scattered near the posterior margin only; and (c) a well-developed paracercus, comprised of more than 50 segments. Primary data on the biology and distribution of this new species are provided, and molecular affinities are verified by the analysis of COI (barcode) sequences. Detailed notes on the distribution of mayfly species belonging to the Baetis alpinus species group common in Western Europe and the western part of North Africa are presented. The historical movement of Baetis representatives between Europe, North West Africa, and subsequently Algeria, with the land bridges 'Strait of Gibraltar' and 'Strait of Sicily' as colonization routes, is discussed in detail and identified in the present study as the Western Algeria colonization path and Eastern Algeria colonization path, respectively.

A new shrew mole species of the genus Uropsilus (Eulipotyphla: Talpidae) from northwestern Vietnam

A new medium-sized shrew mole species of the genus Uropsilus from Mount Fansipan, Hoang Lien National Park, Lao Cai Province, northwestern Vietnam is described based on morphological and molecular differences. Uropsilus fansipanensis sp. nov. is distinguished from the other Uropsilus species by the combination of the following features: the dorsum is lightly reddish-brown and venter is dark gray; the dark gray tail is long and slender, with a scattered white base and short bristle hairs; orbital process is oriented upwards posteriorly; lacrimal foramen is well developed and much larger than infraorbital foramen; the lower first premolar is approximately the same size as the lower third premolar. Genetic distances in terms of mitochondrial cytochrome b from other Uropsilus species presented pairwise divergences from 8.63 to 20.70%. To date, the new species is known to exist only in the type locality of Mt. Fansipan, a wet and cold temperate climate area with an upper montane forest at an elevation of approximately 2900 m, forming the southernmost distribution of the genus Uropsilus.

DNA barcoding and phylogeography of the Hoplias malabaricus species complex

Hoplias malabaricus (Bloch, 1794) is a carnivorous fish species widely distributed from northern to southern South America. This taxon is believed to be a good model for the investigation of biogeographic events that shape the ichthyofauna evolution in the Neotropical freshwater systems. However, many studies have revealed that H. malabaricus hides a species complex that hampers its taxonomic identity and limit its practical value for evolutionary and biogeographic studies. In this paper, we used the mitochondrial gene cytochrome c oxidase subunit I (COI) to delimit cryptic species and explore the phylogeography of H. malabaricus sensu stricto. We found genetic evidence for putative new species in the genus Hoplias and showed that H. malabaricus (Bloch, 1794) is a major clade assigned to barcode index number (BIN) BOLD:ABZ3047. This species is structured in six subpopulations differentiated by high Fst values and restricts gene flow. The subpopulations of the São Francisco/East Atlantic/Eastern Northeast Atlantic/Parnaíba/Itapecuru River basins and Tapajós River Basin were the most differentiated and showed demographic fluctuations. The present distributional pattern is most likely explained through a scenario from the Pleistocene.

Phylogeography of the northernmost distributed Anisocentropus caddisflies and their comparative genetic structures based on habitat preferences

Knowledge of the factors that determine the distribution ranges of organisms is necessary to understand their evolutionary and ecological significance and contribution to biodiversity. A very effective mean of studying such factors is to compare the distribution characteristics and genetic structures of closely related species with differing habitat preferences. Freshwater aquatic insects are relatively easy to observe and the basis of their corresponding niche differentiation easier to identify. Freshwater habitats are categorized lotic or lentic water according to flow regime. In Japanese Islands, the genus Anisocentropus of the calamoceratid caddisfly, the target group in this study, was morphologically reconfirmed that three species, that is, Anisocentropus kawamurai, A. pallidus, and A. magnificus. Among these, A. kawamurai prefers lotic environments and A. pallidus is adapted to lentic water habitats. The distribution range of these sister species overlaps within the Japanese Islands. We estimated the phylogeny and the evolutionary history of Anisocentropus caddisflies worldwide. We estimated divergence periods by two methods, a single locus with various specimens and multiple loci with reduced numbers of the specimens. As a result, we elucidated the phylogenetic position of Japanese species within the cosmopolitan genus Anisocentropus, and also revealed their dual origin. In addition, we demonstrated that the contrasting genetic structures between the sister species distributed in widely overlapping areas were due to differentiation in their respective adapted environmental preferences. Although, in general, it is known that species adapted to lentic water have greater dispersal potential and so are associated with wider distribution areas by means of examining their comparative genetic structures, we revealed a new pattern of genetic locality existing in the genetic structures of the species adapted to lentic water. We then present evidence that suggests the ecological preferences of a species are an important factor in understanding the evolutionary history of that species.

DNA Barcoding and geographical scale effect: The problems of undersampling genetic diversity hotspots

DNA barcoding identification needs a good characterization of intraspecific genetic divergence to establish the limits between species. Yet, the number of barcodes per species is many times low and geographically restricted. A poor coverage of the species distribution range may hamper identification, especially when undersampled areas host genetically distinct lineages. If so, the genetic distance between some query sequences and reference barcodes may exceed the maximum intraspecific threshold for unequivocal species assignation. Taking a group of Quercus herbivores (moths) in Europe as model system, we found that the number of DNA barcodes from southern Europe is proportionally very low in the Barcoding of Life Data Systems. This geographical bias complicates the identification of southern query sequences, due to their high intraspecific genetic distance with respect to barcodes from higher latitudes. Pairwise intraspecific genetic divergence increased along with spatial distance, but was higher when at least one of the sampling sites was in southern Europe. Accordingly, GMYC (General Mixed Yule Coalescent) single-threshold model retrieved clusters constituted exclusively by Iberian haplotypes, some of which could correspond to cryptic species. The number of putative species retrieved was more reliable than that of multiple-threshold GMYC but very similar to results from ABGD and jMOTU. Our results support GMYC as a key resource for species delimitation within poorly inventoried biogeographic regions in Europe, where historical factors (e.g., glaciations) have promoted genetic diversity and singularity. Future European DNA barcoding initiatives should be preferentially performed along latitudinal gradients, with special focus on southern peninsulas.

Exploring the potential of metabarcoding to disentangle macroinvertebrate community dynamics in intermittent streams

Taxonomic sufficiency represents the level of taxonomic detail needed to detect ecological patterns to a level that match the requirement of a study. Most bioassessments apply the taxonomic sufficiency concept and assign specimens to the family or genus level given time constraints and the difficulty to correctly identify species. This holds particularly true for stream invertebrates because small and morphologically similar larvae are hard to distinguish. Low taxonomic resolution may hinder detecting true community dynamics, which thus leads to incorrect inferences about community assembly processes. DNA metabarcoding is a new, affordable and cost-effective tool for the identification of multiple species from bulk samples of organisms. As it provides high taxonomic resolution, it can be used to compare results obtained from different identification levels. Measuring the effect of taxonomic resolution on the detection of community dynamics is especially interesting in extreme ecosystems like intermittent streams to test if species at intermittent sites are subsets of those from perennial sources or if independently recruiting taxa exist. Here we aimed to compare the performance of morphological identification and metabarcoding to detect macroinvertebrate community dynamics in the Trebbia River (Italy). Macroinvertebrates were collected from four perennial and two intermittent sites two months after flow resumption and before the next dry phase. The identification level ranged from family to haplotype. Metabarcoding and morphological identifications found similar alpha diversity patterns when looking at family and mixed taxonomic levels. Increasing taxonomic resolution with metabarcoding revealed a strong partitioning of beta diversity in nestedness and turnover components. At flow resumption, beta diversity at intermittent sites was dominated by nestedness when family-level information was employed, while turnover was evidenced as the most important component when using Operational Taxonomic Units (OTUs) or haplotypes. The increased taxonomic resolution with metabarcoding allowed us to detect species adapted to deal with intermittency, like the chironomid Cricotopus bicinctus and the ephemeropteran Cloeon dipterum. Our study thus shows that family and mixed taxonomic level are not sufficient to detect all aspects of macroinvertebrate community dynamics. High taxonomic resolution is especially important for intermittent streams where accurate information about species-specific habitat preference is needed to interpret diversity patterns induced by drying and the nestedness/turnover components of beta diversity are of interest to understand community assembly processes.

Habitat requirements affect genetic variation in three species of mayfly (Ephemeroptera, Baetidae) from South Africa

This study investigates genetic diversity in three species of Ephemeroptera, one eurytopic and therefore widespread (Afroptilumsudafricanum) and two stenotopic and thus endemic (Demoreptusnatalensis and Demoreptuscapensis) species, all of which co-occur in the southern Great Escarpment, South Africa. Mitochondrial DNA was analysed to compare the genetic diversity between the habitat generalist and the two habitat specialists. Afroptilumsudafricanum showed no indication of population genetic structure due to geographic location, while both Demoreptus species revealed clear genetic differentiation between geographic localities and catchments, evident from phylogenetic analyses and high F_ST values from AMOVA. In addition, the phylogenetic analyses indicate some deeper haplotype divergences within A.sudafricanum and Demoreptus that merit taxonomic attention. These results give important insight into evolutionary processes occurring through habitat specialisation and population isolation. Further research and sampling across a wider geographic setting that includes both major mountain blocks of the Escarpment and lowland non-Escarpment sites will allow for refined understanding of biodiversity and associated habitat preferences, and illuminate comparative inferences into gene flow and cryptic speciation.

Microbial assemblages reflect environmental heterogeneity in alpine streams

Alpine streams are dynamic habitats harboring substantial biodiversity across small spatial extents. The diversity of alpine stream biota is largely reflective of environmental heterogeneity stemming from varying hydrological sources. Globally, alpine stream diversity is under threat as meltwater sources recede and stream conditions become increasingly homogeneous. Much attention has been devoted to macroinvertebrate diversity in alpine headwaters, yet to fully understand the breadth of climate change threats, a more thorough accounting of microbial diversity is needed. We characterized microbial diversity (specifically Bacteria and Archaea) of 13 streams in two disjunct Rocky Mountain subranges through 16S rRNA gene sequencing. Our study encompassed the spectrum of alpine stream sources (glaciers, snowfields, subterranean ice, and groundwater) and three microhabitats (ice, biofilms, and streamwater). We observed no difference in regional (γ) diversity between subranges but substantial differences in diversity among (β) stream types and microhabitats. Within-stream (α) diversity was highest in groundwater-fed springs, lowest in glacier-fed streams, and positively correlated with water temperature for both streamwater and biofilm assemblages. We identified an underappreciated alpine stream type-the icy seep-that are fed by subterranean ice, exhibit cold temperatures (summer mean <2°C), moderate bed stability, and relatively high conductivity. Icy seeps will likely be important for combatting biodiversity losses as they contain similar microbial assemblages to streams fed by surface ice yet may be buffered against climate change by insulating debris cover. Our results show that the patterns of microbial diversity support an ominous trend for alpine stream biodiversity; as meltwater sources decline, stream communities will become more diverse locally, but regional diversity will be lost. Icy seeps, however, represent a source of optimism for the future of biodiversity in these imperiled ecosystems.

Ecosystem shifts in Alpine streams under glacier retreat and rock glacier thaw: A review

This review provides a detailed synthesis of the effects of glacier retreat and permafrost thaw on stream ecosystems in the European Alps. As a working framework, we present a conceptual model developed from an integration of current knowledge and understanding of the habitat and ecological shifts in Alpine streams caused by deglaciation. In our work, we depict how climate change and the loss of cryosphere trigger complex cascading effects on Alpine hydrology, as the main water sources shift from snow and glaciers to rock glaciers, groundwater, and precipitation. The associated changes in habitat conditions, such as channel stability, turbidity, temperature, nutrient loadings, and concentrations of legacy pollutants and trace elements are identified. These changes are followed by complex ecological shifts in the stream communities (microbial community, primary producers, invertebrates) and food webs, with a predicted loss of biotic diversity. Corresponding increases in taxa abundances, biomass, functional diversity, and in the complexity of food webs, are predicted to occur in the upper reaches of Alpine catchments in response to ameliorating climatic and habitat conditions. Finally, current knowledge gaps are highlighted as a basis for framing future research agendas. In particular, we call for an improved understanding of permafrost influence on Alpine headwaters, including the ecology of rock-glacier fed streams, as these streams are likely to become increasingly important for water supply in many glacier-free Alpine valleys in the near future.

Hidden genomic evolution in a morphospecies-The landscape of rapidly evolving genes in Tetrahymena

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei—the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])—within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp–exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

Genomic comparison of ten morphologically very similar species of ciliate from the genus Tetrahymena reveals how parallel microevolutionary processes have shaped their genomes and created unique genes through retrotransposition.

Taxonomy of Baetis Leach in Israel (Ephemeroptera, Baetidae)

The taxonomy and systematics of the genus Baetis Leach (Ephemeroptera: Baetidae) in Israel is clarified for the first time as part of an ongoing comprehensive survey of the Israeli mayfly fauna. Six clearly defined species are currently recognized in Israel, four of which are described here as new to science. The validity of all species is supported by a molecular analysis of the mitochondrial COI gene. A key for the identification of the Israeli species based on the nymphal stage is provided, as well as data on their distribution patterns and ecology. The local fauna represents four Palearctic species groups, three of which reach their limit of distribution range in Israel. Baetis species typically inhabit lotic, pristine habitats in northern Israel, with few exceptions for some species that can be found in the Dead Sea area or in ephemeral ponds.

Comparison of environmental DNA and bulk-sample metabarcoding using highly degenerate cytochrome c oxidase I primers

Freshwater biodiversity provides important ecosystem services and is at the core of water quality monitoring worldwide. To assess freshwater biodiversity, genetic methods such as metabarcoding are increasingly used as they are faster and allow better taxonomic resolution than manual identification methods. Either sampled organisms are used directly for "bulk metabarcoding," or water is filtered and the extracted environmental DNA serves as a proxy for biodiversity via "eDNA metabarcoding." Despite the advantages of both methods, questions remain regarding their comparability and applicability for routine biomonitoring and stressor impact assessment. Therefore, we compared metabarcoding results from bulk and eDNA samples taken from 19 streams spanning a wide gradient of farming intensities in New Zealand. We performed PCR with highly degenerate cytochrome c oxidase I primers and sequenced libraries on an Illumina MiSeq. The inferred community composition differed strongly between the two methods. More taxa were captured by eDNA than bulk-sample metabarcoding (5,819 vs. 1,483), but more of the commonly used invertebrate bioindicator taxa (mayflies, stoneflies and caddisflies) were found in bulk (47) than eDNA samples (37). Catchment-wide and local land use impacts on communities were detected better by eDNA metabarcoding, especially for non-metazoan taxa. Our findings imply that bulk-sample metabarcoding resembles classical freshwater biomonitoring approaches better, as more indicator macroinvertebrate taxa are captured. However, eDNA metabarcoding might be better suited to infer the impact of stressors on stream ecosystems at larger scales, as many new and potentially more informative taxa are registered. We therefore suggest exploring both methods in future assessments of stream biodiversity.

The Trouble with MEAM2: Implications of Pseudogenes on Species Delimitation in the Globally Invasive Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species Complex

Molecular species identification using suboptimal PCR primers can over-estimate species diversity due to coamplification of nuclear mitochondrial (NUMT) DNA/pseudogenes. For the agriculturally important whitefly Bemisia tabaci cryptic pest species complex, species identification depends primarily on characterization of the mitochondrial DNA cytochrome oxidase I (mtDNA COI) gene. The lack of robust PCR primers for the mtDNA COI gene can undermine correct species identification which in turn compromises management strategies. This problem is identified in the B. tabaci Africa/Middle East/Asia Minor clade which comprises the globally invasive Mediterranean (MED) and Middle East Asia Minor I (MEAM1) species, Middle East Asia Minor 2 (MEAM2), and the Indian Ocean (IO) species. Initially identified from the Indian Ocean island of Réunion, MEAM2 has since been reported from Japan, Peru, Turkey and Iraq. We identified MEAM2 individuals from a Peruvian population via Sanger sequencing of the mtDNA COI gene. In attempting to characterize the MEAM2 mitogenome, we instead characterized mitogenomes of MEAM1. We also report on the mitogenomes of MED, AUS, and IO thereby increasing genomic resources for members of this complex. Gene synteny (i.e., same gene composition and orientation) was observed with published B. tabaci cryptic species mitogenomes. Pseudogene fragments matching MEAM2 partial mtDNA COI gene exhibited low frequency single nucleotide polymorphisms that matched low copy number DNA fragments (<3%) of MEAM1 genomes, whereas presence of internal stop codons, loss of expected stop codons and poor primer annealing sites, all suggested MEAM2 as a pseudogene artifact and so not a real species.

Climate-induced glacier and snow loss imperils alpine stream insects

Climate warming is causing rapid loss of glaciers and snowpack in mountainous regions worldwide. These changes are predicted to negatively impact the habitats of many range-restricted species, particularly endemic, mountaintop species dependent on the unique thermal and hydrologic conditions found only in glacier-fed and snow melt-driven alpine streams. Although progress has been made, existing understanding of the status, distribution, and ecology of alpine aquatic species, particularly in North America, is lacking, thereby hindering conservation and management programs. Two aquatic insects - the meltwater stonefly (Lednia tumana) and the glacier stonefly (Zapada glacier) - were recently proposed for listing under the U.S. Endangered Species Act due to climate-change-induced habitat loss. Using a large dataset (272 streams, 482 total sites) with high-resolution climate and habitat information, we describe the distribution, status, and key environmental features that limit L. tumana and Z. glacier across the northern Rocky Mountains. Lednia tumana was detected in 113 streams (175 sites) within Glacier National Park (GNP) and surrounding areas. The probability of L. tumana occurrence increased with cold stream temperatures and close proximity to glaciers and permanent snowfields. Similarly, densities of L. tumana declined with increasing distance from stream source. Zapada glacier was only detected in 10 streams (24 sites), six in GNP and four in mountain ranges up to ~600 km southwest. Our results show that both L. tumana and Z. glacier inhabit an extremely narrow distribution, restricted to short sections of cold, alpine streams often below glaciers predicted to disappear over the next two decades. Climate warming-induced glacier and snow loss clearly imperils the persistence of L. tumana and Z. glacier throughout their ranges, highlighting the role of mountaintop aquatic invertebrates as sentinels of climate change in mid-latitude regions.

Baetis (Baetis) cypronyx sp. n., a new species of the Baetis alpinus species-group (Insecta, Ephemeroptera, Baetidae) from Cyprus, with annotated checklist of Baetidae in the Mediterranean islands

A detailed description of the larvae of Baetis (Baetis) cypronyxsp. n., a representative of the Baetisalpinus species-group within the mayfly family Baetidae, is provided, including a differential diagnosis with regard to closely related species of the group, especially Baetismelanonyx (Pictet, 1843) and Baetisbaroukianus Thomas & Dia, 1984. The new species is mainly distinguished by mouthparts (i.e. the shape and setation of labrum, maxillary and labial palps, details of paraglossae and mandibular incisors), setation of legs and abdominal terga, and length of paracercus. All available data on the biology of this putative endemic species of Cyprus are summarized. Annotated distributional data of the 33 species of Baetidae so far recorded from the Mediterranean islands are given, including new records and also including first data from Malta.

Climate change and alpine stream biology: progress, challenges, and opportunities for the future

In alpine regions worldwide, climate change is dramatically altering ecosystems and affecting biodiversity in many ways. For streams, receding alpine glaciers and snowfields, paired with altered precipitation regimes, are driving shifts in hydrology, species distributions, basal resources, and threatening the very existence of some habitats and biota. Alpine streams harbour substantial species and genetic diversity due to significant habitat insularity and environmental heterogeneity. Climate change is expected to affect alpine stream biodiversity across many levels of biological resolution from micro- to macroscopic organisms and genes to communities. Herein, we describe the current state of alpine stream biology from an organism-focused perspective. We begin by reviewing seven standard and emerging approaches that combine to form the current state of the discipline. We follow with a call for increased synthesis across existing approaches to improve understanding of how these imperiled ecosystems are responding to rapid environmental change. We then take a forward-looking viewpoint on how alpine stream biologists can make better use of existing data sets through temporal comparisons, integrate remote sensing and geographic information system (GIS) technologies, and apply genomic tools to refine knowledge of underlying evolutionary processes. We conclude with comments about the future of biodiversity conservation in alpine streams to confront the daunting challenge of mitigating the effects of rapid environmental change in these sentinel ecosystems.

Colonization and diversification of aquatic insects on three Macaronesian archipelagos using 59 nuclear loci derived from a draft genome

The study of processes driving diversification requires a fully sampled and well resolved phylogeny, although a lack of phylogenetic markers remains a limitation for many non-model groups. Multilocus approaches to the study of recent diversification provide a powerful means to study the evolutionary process, but their application remains restricted because multiple unlinked loci with suitable variation for phylogenetic or coalescent analysis are not available for most non-model taxa. Here we identify novel, putative single-copy nuclear DNA (nDNA) phylogenetic markers to study the colonization and diversification of an aquatic insect species complex, Cloeon dipterum L. 1761 (Ephemeroptera: Baetidae), in Macaronesia. Whole-genome sequencing data from one member of the species complex were used to identify 59 nDNA loci (32,213 base pairs), followed by Sanger sequencing of 29 individuals sampled from 13 islands of three Macaronesian archipelagos. Multispecies coalescent analyses established six putative species. Three island species formed a monophyletic clade, with one species occurring on the Azores, Europe and North America. Ancestral state reconstruction indicated at least two colonization events from the mainland (to the Canaries, respectively Azores) and one within the archipelago (between Madeira and the Canaries). Random subsets of the 59 loci showed a positive linear relationship between number of loci and node support. In contrast, node support in the multispecies coalescent tree was negatively correlated with mean number of phylogenetically informative sites per locus, suggesting a complex relationship between tree resolution and marker variability. Our approach highlights the value of combining genomics, coalescent-based phylogeography, species delimitation, and phylogenetic reconstruction to resolve recent diversification events in an archipelago species complex.

discomark: nuclear marker discovery from orthologous sequences using draft genome data

High-throughput sequencing has laid the foundation for fast and cost-effective development of phylogenetic markers. Here we present the program discomark, which streamlines the development of nuclear DNA (nDNA) markers from whole-genome (or whole-transcriptome) sequencing data, combining local alignment, alignment trimming, reference mapping and primer design based on multiple sequence alignments to design primer pairs from input orthologous sequences. To demonstrate the suitability of discomark, we designed markers for two groups of species, one consisting of closely related species and one group of distantly related species. For the closely related members of the species complex of Cloeon dipterum s.l. (Insecta, Ephemeroptera), the program discovered a total of 78 markers. Among these, we selected eight markers for amplification and Sanger sequencing. The exon sequence alignments (2526 base pairs) were used to reconstruct a well-supported phylogeny and to infer clearly structured haplotype networks. For the distantly related species, we designed primers for the insect order Ephemeroptera, using available genomic data from four sequenced species. We developed primer pairs for 23 markers that are designed to amplify across several families. The discomark program will enhance the development of new nDNA markers by providing a streamlined, automated approach to perform genome-scale scans for phylogenetic markers. The program is written in Python, released under a public licence (GNU GPL version 2), and together with a manual and example data set available at: https://github.com/hdetering/discomark.