True Parthenogenesis and Female-Biased Sex Ratios in Cicadomorpha and Fulgoromorpha (Hemiptera, Auchenorrhyncha)

Insects are renowned for their remarkable diversity of reproductive modes. Among these, the largest non-holometabolous order, Hemiptera, stands out with one of the most diversified arrays of parthenogenesis modes observed among insects. Although there are extensive reviews on reproduction without fertilization in some hemipteran higher taxa, no such analysis has been conducted for the large suborders Fulgoromorpha (planthoppers) and Cicadomorpha (leafhoppers). In both groups, there are species that reproduce by true parthenogenesis, specifically thelytoky, and in Fulgoromorpha, there are species that reproduce by pseudogamy or, more specifically, sperm-dependent parthenogenesis. In this review paper, we give and discuss the only currently known examples of true parthenogenesis in Fulgoromorpha and Cicadomorpha, mainly from the planthopper family Delphacidae and the leafhopper family Cicadellidae. We analyze patterns of distribution, ecology, mating behavior, acoustic communication, and cytogenetic and genetic diversity of parthenoforms and discuss hypotheses about the origin of parthenogenesis in each case. We also highlight examples in which natural populations show a shift in sex ratio toward females and discuss possible causes of this phenomenon, primarily the influence of endosymbiotic bacteria capable of altering the reproductive strategies of the hosts. Our review is mainly based on studies in which the authors have participated.

First annotated draft genomes of nonmarine ostracods (Ostracoda, Crustacea) with different reproductive modes

Ostracods are one of the oldest crustacean groups with an excellent fossil record and high importance for phylogenetic analyses but genome resources for this class are still lacking. We have successfully assembled and annotated the first reference genomes for three species of nonmarine ostracods; two with obligate sexual reproduction (Cyprideis torosa and Notodromas monacha) and the putative ancient asexual Darwinula stevensoni. This kind of genomic research has so far been impeded by the small size of most ostracods and the absence of genetic resources such as linkage maps or BAC libraries that were available for other crustaceans. For genome assembly, we used an Illumina-based sequencing technology, resulting in assemblies of similar sizes for the three species (335-382 Mb) and with scaffold numbers and their N50 (19-56 kb) in the same orders of magnitude. Gene annotations were guided by transcriptome data from each species. The three assemblies are relatively complete with BUSCO scores of 92-96. The number of predicted genes (13,771-17,776) is in the same range as Branchiopoda genomes but lower than in most malacostracan genomes. These three reference genomes from nonmarine ostracods provide the urgently needed basis to further develop ostracods as models for evolutionary and ecological research.

Skipping sex: A nonrecombinant genomic assemblage of complementary reproductive modules

The unusual occurrence and developmental diversity of asexual eukaryotes remain a puzzle. De novo formation of a functioning asexual genome requires a unique assembly of sets of genes or gene states to disrupt cellular mechanisms of meiosis and gametogenesis, and to affect discrete components of sexuality and produce clonal or hemiclonal offspring. We highlight two usually overlooked but essential conditions to understand the molecular nature of clonal organisms, that is, a nonrecombinant genomic assemblage retaining modifiers of the sexual program, and a complementation between altered reproductive components. These subtle conditions are the basis for physiologically viable and genetically balanced transitions between generations. Genomic and developmental evidence from asexual animals and plants indicates the lack of complementation of molecular changes in the sexual reproductive program is likely the main cause of asexuals' rarity, and can provide an explanatory frame for the developmental diversity and lability of developmental patterns in some asexuals as well as for the discordant time to extinction estimations.

Karyotype Variability and Inter-Population Genomic Differences in Freshwater Ostracods (Crustacea) Showing Geographical Parthenogenesis

Transitions from sexual to asexual reproduction are often associated with polyploidy and increased chromosomal plasticity in asexuals. We investigated chromosomes in the freshwater ostracod species Eucypris virens (Jurine, 1820), where sexual, asexual and mixed populations can be found. Our initial karyotyping of multiple populations from Europe and North Africa, both sexual and asexual, revealed a striking variability in chromosome numbers. This would suggest that chromosomal changes are likely to be accelerated in asexuals because the constraints of meiosis are removed. Hence, we employed comparative genomic hybridization (CGH) within and among sexual and asexual populations to get insights into E. virens genome arrangements. CGH disclosed substantial genomic imbalances among the populations analyzed, and three patterns of genome arrangement between these populations: 1. Only putative ribosomal DNA (rDNA)-bearing regions were conserved in the two populations compared indicating a high sequence divergence between these populations. This pattern is comparable with our findings at the interspecies level of comparison; 2. Chromosomal regions were shared by both populations to a varying extent with a distinct copy number variation in pericentromeric and presumable rDNA-bearing regions. This indicates a different rate of evolution in repetitive sequences; 3. A mosaic pattern of distribution of genomic material that can be explained as non-reciprocal genetic introgression and evidence of a hybrid origin of these individuals. We show an overall increased chromosomal dynamics in E. virens that is complementary with available phylogenetic and population genetic data reporting highly differentiated diploid sexual and asexual lineages with a wide variety of genetic backgrounds.

Aquatic long-distance dispersal and vicariance shape the evolution of an ostracod species complex (Crustacea) in four major Brazilian floodplains

Cladogenesis is often driven by the interplay of dispersal and vicariance. The importance of long-distance dispersal in biogeography and speciation is increasingly recognised, but still ill-understood. Here, we study faunal interconnectivity between four large Brazilian floodplains, namely the Amazon, Araguaia, Pantanal (on Paraguay River) and Upper Paraná River floodplains, investigating a species complex of the non-marine ostracod genus Strandesia. We use DNA sequence data from the mitochondrial COI and the nuclear Elongation Factor 1 alpha genes to construct molecular phylogenies and minimum spanning networks, to identify genetic species, analyse biogeographic histories and provide preliminary age estimates of this species complex. The Strandesia species complex includes five morphological and eleven genetic species, which doubles the known diversity in this lineage. The evolutionary history of this species complex appears to comprise sequences of dispersal and vicariance events. Faunal and genetic patterns of connectivity between floodplains in some genetic species are mirrored in modern hydrological connections. This could explain why we find evidence for (aquatic) long-distance dispersal between floodplains, thousands of kilometres apart. Our phylogenetic reconstructions seem to mostly indicate recent dispersal and vicariance events, but the evolution of the present Strandesia species complex could span up to 25 Myr, which by far exceeds the age of the floodplains and the rivers in their current forms.

Karyotypic variation and pollen stainability in resynthesized allopolyploids Tragopogon miscellus and T. mirus

PREMISE OF THE STUDY

Polyploidy has extensively shaped the evolution of plants, but the early stages of polyploidy are still poorly understood. The neoallopolyploid species Tragopogon mirus and T. miscellus are both characterized by widespread karyotypic variation, including frequent aneuploidy and intergenomic translocations. Our study illuminates the origins and early impacts of this variation by addressing two questions: How quickly does karyotypic variation accumulate in Tragopogon allopolyploids following whole-genome duplication (WGD), and how does the fertility of resynthesized Tragopogon allopolyploids evolve shortly after WGD?

METHODS

We used genomic in situ hybridization and lactophenol-cotton blue staining to estimate the karyotypic variation and pollen stainability, respectively, of resynthesized T. mirus and T. miscellus during the first five generations after WGD.

KEY RESULTS

Widespread karyotypic variation developed quickly in synthetics and resembled that of naturally occurring T. mirus and T. miscellus by generation S₄ . Pollen stainability in resynthesized allopolyploids was consistently lower than that of natural T. mirus and T. miscellus, as well as their respective diploid progenitor species. Logistic regression showed that mean pollen stainability increased slightly over four generations in resynthesized T. mirus but remained at equivalent levels in T. miscellus.

CONCLUSIONS

Our results clarify some of the changes that occur in T. mirus and T. miscellus immediately following their origin, most notably the rapid onset of karyotypic variation within these species immediately following WGD.

The thermodynamics of protein aggregation reactions may underpin the enhanced metabolic efficiency associated with heterosis, some balancing selection, and the evolution of ploidy levels

Identifying the physical basis of heterosis (or "hybrid vigor") has remained elusive despite over a hundred years of research on the subject. The three main theories of heterosis are dominance theory, overdominance theory, and epistasis theory. Kacser and Burns (1981) identified the molecular basis of dominance, which has greatly enhanced our understanding of its importance to heterosis. This paper aims to explain how overdominance, and some features of epistasis, can similarly emerge from the molecular dynamics of proteins. Possessing multiple alleles at a gene locus results in the synthesis of different allozymes at reduced concentrations. This in turn reduces the rate at which each allozyme forms soluble oligomers, which are toxic and must be degraded, because allozymes co-aggregate at low efficiencies. The model developed in this paper can explain how heterozygosity impacts the metabolic efficiency of an organism. It can also explain why the viabilities of some inbred lines seem to decline rapidly at high inbreeding coefficients (F > 0.5), which may provide a physical basis for truncation selection for heterozygosity. Finally, the model has implications for the ploidy level of organisms. It can explain why polyploids are frequently found in environments where severe physical stresses promote the formation of soluble oligomers. The model can also explain why complex organisms, which need to synthesize aggregation-prone proteins that contain intrinsically unstructured regions (IURs) and multiple domains because they facilitate complex protein interaction networks (PINs), tend to be diploid while haploidy tends to be restricted to relatively simple organisms.

Seasonal niche partitioning and coexistence of amphimictic and parthenogenetic lineages of Heterocypris barbara (Crustacea: Ostracoda)

Sympatry of amphimictic and parthenogenetic lineages in species with mixed reproductive systems is rarely observed in nature. On Lampedusa Island (Pelagie Islands, Italy), amphimictic and parthenogenetic lineages of Heterocypris barbara (Gauthier and Brehm, 1928) co-occur in a temporary pond. Their sympatric persistence calls for an ecological differentiation. We investigated the role of seasonal variation of temperature and photoperiod conditions by two different approaches: microcosms set up by inundation of dry sediments from the temporary pond and life-table experiments. Microcosms recreate conditions similar to the field and in their sediments random samples of resting eggs of both amphimictic and parthenogenetic females are stored. Life-table experiments supplied individual-based estimates of survivorship, adult life span, fecundity, and sex ratio in the progeny. We carried out the experiments at 24 °C and a photoperiod of 12 h light (L) : 12 h dark (D) (simulating fall conditions) and at 16 °C and a photoperiod of 10 h L : 14 h D (simulating winter conditions). Males and amphimictic females were the most numerous forms at 24 °C and 12 h L : 12 h D; parthenogenetic females were dominant at 16 °C and 10 h L : 14 h D. Life-table experiments showed that amphimictic forms do not complete development at 16 °C and 10 h L : 14 h D. Our results suggest that sympatry of amphimictic and parthenogenetic females in the field depends on seasonal niche partitioning and the storage effect of resting eggs that allows survival through adverse-season conditions.

Ostracod (Ostracoda, Crustacea) genomics - Promises and challenges

Ostracods are well-suited model organisms for evolutionary research. Classic genetic techniques have mostly been used for phylogenetic studies on Ostracoda and were somewhat affected by the lack of large numbers of suitable markers. Genomic methods with their huge potential have so far rarely been applied to this group of crustaceans. We provide relevant examples of genomic studies on other organisms to propose future avenues of genomic ostracod research. At the same time, we suggest solutions to the potential problems in ostracods that the application of genomic techniques might present.

Chromosome numbers in three species groups of freshwater flatworms increase with increasing latitude

Polyploidy in combination with parthenogenesis offers advantages for plasticity and the evolution of a broad ecological tolerance of species. Therefore, a positive correlation between the level of ploidy and increasing latitude as a surrogate for environmental harshness has been suggested. Such a positive correlation is well documented for plants, but examples for animals are still rare. Species of flatworms (Platyhelminthes) are widely distributed, show a remarkably wide range of chromosome numbers, and offer therefore good model systems to study the geographical distribution of chromosome numbers. We analyzed published data on counts of chromosome numbers and geographical information of three flatworm "species" (Phagocata vitta, Polycelis felina and Crenobia alpina) sampled across Europe (220 populations). We used the mean chromosome number across individuals of a population as a proxy for the level of ploidy within populations, and we tested for relationships of this variable with latitude, mode of reproduction (sexual, asexual or both) and environmental variables (annual mean temperature, mean diurnal temperature range, mean precipitation and net primary production). The mean chromosome numbers of all three species increased with latitude and decreased with mean annual temperature. For two species, chromosome number also decreased with mean precipitation and net primary production. Furthermore, high chromosome numbers within species were accompanied with a loss of sexual reproduction. The variation of chromosome numbers within individuals of two of the three species increased with latitude. Our results support the hypothesis that polyploid lineages are able to cope with harsh climatic conditions at high latitudes. Furthermore, we propose that asexual reproduction in populations with high levels of polyploidization stabilizes hybridization events. Chromosomal irregularities within individuals tend to become more frequent at the extreme environments of high latitudes, presumably because of mitotic errors and downsizing of the genome.

Extensive variation in chromosome number and genome size in sexual and parthenogenetic species of the jumping-bristletail genus Machilis (Archaeognatha)

Parthenogenesis in animals is often associated with polyploidy and restriction to extreme habitats or recently deglaciated areas. It has been hypothesized that benefits conferred by asexual reproduction and polyploidy are essential for colonizing these habitats. However, while evolutionary routes to parthenogenesis are manifold, study systems including polyploids are scarce in arthropods. The jumping-bristletail genus Machilis (Insecta: Archaeognatha) includes both sexual and parthenogenetic species, and recently, the occurrence of polyploidy has been postulated. Here, we applied flow cytometry, karyotyping, and mitochondrial DNA sequencing to three sexual and five putatively parthenogenetic Eastern-Alpine Machilis species to investigate whether (1) parthenogenesis originated once or multiply and (2) whether parthenogenesis is strictly associated with polyploidy. The mitochondrial phylogeny revealed that parthenogenesis evolved at least five times independently among Eastern-Alpine representatives of this genus. One parthenogenetic species was exclusively triploid, while a second consisted of both diploid and triploid populations. The three other parthenogenetic species and all sexual species were diploid. Our results thus indicate that polyploidy can co-occur with parthenogenesis, but that it was not mandatory for the emergence of parthenogenesis in Machilis. Overall, we found a weak negative correlation of monoploid genome size (Cx) and chromosome base number (x), and this connection is stronger among parthenogenetic species alone. Likewise, monoploid genome size decreased with elevation, and we therefore hypothesize that genome downsizing could have been crucial for the persistence of alpine Machilis species. Finally, we discuss the evolutionary consequences of intraspecific chromosomal rearrangements and the presence of B chromosomes. In doing so, we highlight the potential of Alpine Machilis species for research on chromosomal and genome-size alterations during speciation.

Genetic uniformity and long-distance clonal dispersal in the invasive androgenetic Corbicula clams

The clam genus Corbicula is an interesting model system to study the evolution of reproductive modes as it includes both sexual and asexual (androgenetic) lineages. While the sexual populations are restricted to the native Asian areas, the androgenetic lineages are widely distributed being also found in America and Europe where they form a major aquatic invasive pest. We investigated the genetic diversity of native and invasive Corbicula populations through a worldwide sampling. The use of mitochondrial and nuclear (microsatellite) markers revealed an extremely low diversity in the invasive populations with only four, undiversified, genetic lineages distributed across Europe and America. On the contrary, in the native populations, both sexual and androgenetic lineages exhibited much higher genetic diversity. Remarkably, the most abundant and widely distributed invasive forms, the so-called form A and form R found in America and Europe respectively, are fixed for the same single COI (cytochrome c oxydase subunit I) haplotype and same multilocus genotype. This suggests that form R, observed in Europe since the 1980s, derived directly from form A found in America since the 1920s. In addition, this form shares alleles with some Japanese populations, indicating a Japanese origin for this invasive lineage. Finally, our study suggests that few androgenetic Corbicula individuals successfully invaded the non-native range and then dispersed clonally. This is one striking case of genetic paradox raising the issue of invasive and evolutionary success of genetically undiversified populations.

Population genetics of bisexual and unisexual populations of the scaly-winged bark louse Echmepteryx hageni (Insecta: Psocoptera)

The scaly-winged bark louse, Echmepteryx hageni, exhibits a unique pattern of co-existence of apparently differnt reproductive modes. Unisexuality is widespread in eastern North America, while sexual populations are restricted to isolated rock out-croppings in southern Illinois and eastern Kentucky. Three of the four nuclear loci examined show greater genetic diversity in the unisexual form compared to the sexual form of E. hageni, in accordance with the pattern previously shown in mitochondrial genetic data. Neutrality tests of the nuclear loci indicate a consistent signal of demographic expansion in asexual populations, but not in sexual populations. There was evidence of inbreeding in the isolated sexual populations at three of the nuclear loci, and one locus had signs of gene specific balancing selection. However, there is no significant genetic differentiation between bisexual and unisexual populations, possibly due to the greater effective population size of nuclear loci relative to mitochondrial loci. The mitochondrial differentiation of E. hageni populations in the northwestern part of their range (Minnesota and Wisconsin) was also not reflected in the nuclear data. We present three hypotheses that may explain the disparity in observed nuclear and mitochondrial genetic diversity between the reproductive forms of E. hageni.

Alteration of sexual reproduction and genetic diversity in the kelp species Laminaria digitata at the southern limit of its range

Adaptation to marginal habitats at species range-limits has often been associated with parthenogenetic reproduction in terrestrial animals and plants. Laboratory observations have shown that brown algae exhibit a high propensity for parthenogenesis by various mechanisms. The kelp Laminaria digitata is an important component of the ecosystem in Northern European rocky intertidal habitats. We studied four L. digitata populations for the effects of marginality on genetic diversity and sexual reproduction. Two populations were marginal: One (Locquirec, in Northern Brittany) was well within the geographic range, but was genetically isolated from other populations by large stretches of sandy beaches. Another population was at the range limits of the species (Quiberon, in Southern Brittany) and was exposed to much higher seasonal temperature changes. Microsatellite analyses confirmed that these populations showed decreased genetic and allelic diversity, consistent with marginality and genetic isolation. Sporophytes from both marginal populations showed greatly diminished spore-production compared to central populations, but only the southern-limit population (Quiberon) showed a high propensity for producing unreduced (2N) spores. Unreduced 2N spores formed phenotypically normal gametophytes with nuclear area consistent with ≥2N DNA contents, and microsatellite studies suggested these were produced at least in part by automixis. However, despite this being the dominant path of spore production in Quiberon sporophyte individuals, the genetic evidence indicated the population was maintained mostly by sexual reproduction. Thus, although spore production and development showed the expected tendency of geographical parthenogenesis in marginal populations, this appeared to be a consequence of maladaptation, rather than an adaptation to, life in a marginal habitat.

Morphotypes of Cytheridella Ilosvayi (Ostracoda) Detected by Soft and Hard Part Analyses

Quantitative approaches dealing with soft- and hard-part morphological variability of ostracods are very rare for most species. Here we present intraspecific limb and shell variability of the Neotropical freshwater ostracod species assumed as Cytheridella ilosvayi. For the first time such an analysis considered adult females and males, as well as juveniles. Variability of the appendages included the antennule, second antenna as well as the first and second thoracopod and was analysed in terms of ratios of podomeres, setae and limbs. Limb variability of adult and juvenile individuals (down to A-3) is generally low. The highest variation is shown in podomere proportions of the antennas, while thoracopods and setae provide minor and/or non-significant influence on the variability. Based on discrimination analyses shell parameters (i.e., shell length, position of the transversal sulcus) emerge to be more important for differentiation of groups than limb ratios. Adult females exhibit a large size range in which two clearly separated morphotypes exist. Although differences in limbs are small their existence can be proven also on the basis of proportions between podomeres of the antennas. Males and juveniles, however, do not display a separation into corresponding subgroups or morphotypes. The presence of two morphologically similar females and only one type of males indicates the coexistence of female morphotypes which may represent either two (cryptic) species or a mixed reproduction population in which parthenogenetic and sexual reproduction coexists.

The distribution of sexual reproduction of the geographic parthenogen Eucypris virens (Crustacea: Ostracoda) matches environmental gradients in a temporary lake

Niche segregation may prevent competitive exclusion and promote local coexistence. This typically results in the occupation of different habitats. In the freshwater ostracod Eucypris virens (Jurine, 1820), the distribution of sexual and parthenogenetic populations in the temporary Lake Caracuel, central Spain, was not homogeneous. Parthenogens were found everywhere including the littoral, whereas sexuals were restricted to the centre. We investigated the hypothesis that spatial distribution responded to ecological differences. We studied the ecological significance of this segregation by linking environmental data to male presence, sexual fraction, ploidy, and genetic structure of our model organism in 12 ponds in the lake basin. Hydro-chemical and biological data indicated that the observed segregation is not only spatial but also ecological, with sexual E. virens occurring in ponds with shorter and probably more unpredictable hydroperiod. The correlations between environmental gradients and E. virens population structure suggest that sexual and parthenogenetic lineages are ecologically segregated. Sexual populations appeared restricted to sufficiently unpredictable environments, while parthenogens dominated environments with longer hydroperiods. Local coexistence seems mediated by spatial heterogeneity in habitat stability.

Can resource costs of polyploidy provide an advantage to sex?

The predominance of sexual reproduction despite its costs indicates that sex provides substantial benefits, which are usually thought to derive from the direct genetic consequences of recombination and syngamy. While genetic benefits of sex are certainly important, sexual and asexual individuals, lineages, or populations may also differ in physiological and life history traits that could influence outcomes of competition between sexuals and asexuals across environmental gradients. Here, we address possible phenotypic costs of a very common correlate of asexuality, polyploidy. We suggest that polyploidy could confer resource costs related to the dietary phosphorus demands of nucleic acid production; such costs could facilitate the persistence of sex in situations where asexual taxa are of higher ploidy level and phosphorus availability limits important traits like growth and reproduction. We outline predictions regarding the distribution of diploid sexual and polyploid asexual taxa across biogeochemical gradients and provide suggestions for study systems and empirical approaches for testing elements of our hypothesis.

Geographic parthenogenesis and plant-enemy interactions in the common dandelion

BACKGROUND

Many species with sexual and asexual variants show a pattern of geographic parthenogenesis where asexuals have broader and higher-latitude distribution than sexuals. Because sexual reproduction is often considered a costly evolutionary strategy that is advantageous in the face of selection by coevolving pests and pathogens, one possible explanation for geographic parthenogenesis is that populations at higher latitudes are exposed to fewer pests and pathogens. We tested this hypothesis in the common dandelion (Taraxacum officinale), a species with well-established geographic parthenogenesis, by screening prevalence and effects of several specialized pests and pathogens in natural dandelion populations.

RESULTS

We did a population survey of 18 dandelion populations along a geographic transect that ranged from the area where sexual and asexual dandelions co-occur northward into the area where only asexuals occur. In addition we used four southern and four northern populations in a 8x8 cross-inoculation greenhouse experiment in which plants were exposed experimentally to each other's natural field soil microbial communities. The cross-inoculation experiment indicated a higher pathogenicity of soil microbial communities from the southern, mostly sexual, populations compared to soil microbial communities from the northern asexual populations. Northern dandelion populations also showed reduced infestation by a specialized seed-eating weevil. A similar trend of reduced rust fungus infection in northern populations was observed but this trend was not statistically significant.

CONCLUSIONS

The prevalence of pests and pathogens decreased along the south-to-north axis of geographic parthenogenesis. This highlights the potential of biotic interactions in shaping patterns of geographic parthenogenesis.

Toward a DNA taxonomy of Alpine Rhithrogena (Ephemeroptera: Heptageniidae) using a mixed Yule-coalescent analysis of mitochondrial and nuclear DNA

Aquatic larvae of many Rhithrogena mayflies (Ephemeroptera) inhabit sensitive Alpine environments. A number of species are on the IUCN Red List and many recognized species have restricted distributions and are of conservation interest. Despite their ecological and conservation importance, ambiguous morphological differences among closely related species suggest that the current taxonomy may not accurately reflect the evolutionary diversity of the group. Here we examined the species status of nearly 50% of European Rhithrogena diversity using a widespread sampling scheme of Alpine species that included 22 type localities, general mixed Yule-coalescent (GMYC) model analysis of one standard mtDNA marker and one newly developed nDNA marker, and morphological identification where possible. Using sequences from 533 individuals from 144 sampling localities, we observed significant clustering of the mitochondrial (cox1) marker into 31 GMYC species. Twenty-one of these could be identified based on the presence of topotypes (expertly identified specimens from the species' type locality) or unambiguous morphology. These results strongly suggest the presence of both cryptic diversity and taxonomic oversplitting in Rhithrogena. Significant clustering was not detected with protein-coding nuclear PEPCK, although nine GMYC species were congruent with well supported terminal clusters of nDNA. Lack of greater congruence in the two data sets may be the result of incomplete sorting of ancestral polymorphism. Bayesian phylogenetic analyses of both gene regions recovered four of the six recognized Rhithrogena species groups in our samples as monophyletic. Future development of more nuclear markers would facilitate multi-locus analysis of unresolved, closely related species pairs. The DNA taxonomy developed here lays the groundwork for a future revision of the important but cryptic Rhithrogena genus in Europe.