Skip to main content

Fast and accurate identification of cryptic and sympatric mayfly species of the Baetis rhodani group



Species of the Baetis rhodani group are among the most widespread mayflies of the Palearctic region. However, frequent occurrence of morphologically cryptic species complicates the identification of sympatric species. Here, we proposed and tested a method for the fast, accurate, and cost-effective assignment of a large number of individuals to their putative species, based on high resolution melting profiles of a standard mitochondrial gene fragment. We tested this method using a system of three recently identified cryptic species inhabiting the Tyrrhenian Islands (western Mediterranean basin).


Highly species-specific high resolution melting profiles were obtained, allowing the unequivocal attribution of each individual to the respective species. This assay provides a convenient and easily customizable alternative to traditional barcoding approaches, provided that the mayfly taxa occurring within the geographic area of interest have been previously identified and their high resolution melting profiles assessed.


Baetis rhodani (Pictet 1843) is a complex of morphologically cryptic mayfly species widespread in the Palearctic region [1]. Despite a considerable research effort in recent years [2,3,4,5], current knowledge of the phylogenetic relationships within this group and the taxonomic status of several of its representatives remain incomplete, with new species emerging as new areas of the Palearctic are investigated.

Methods of species recognition and individual assignment in mayflies have been mostly based on two distinct approaches, each one with its own advantages and drawbacks. The first and more traditional approach is based on the analysis of morphological diagnostic traits, primarily of the larval stages [6, 7]. Although relatively inexpensive and time-saving, this approach relies entirely on an expert knowledge of the taxonomic group under study. Most importantly, with the advent of the second approach, molecular taxonomy, morphological approaches have soon appeared largely to underestimate the amount of species diversity within baetidae mayflies, including B. rhodani [2]. Indeed, under the single morphological species B. rhodani, several deeply divergent species have recently been recognized using molecular markers [2,3,4,5], and even more are likely to be described in the near future. On the other hand, molecular methods, such as those based on DNA barcodes, entail unprecedented processing costs per sample, and their extensive use for the preliminary phase of species delimitation has been largely debated (see [8, 9]).

Recently, three endemic species belonging to the B. rhodani species group have been identified within the Tyrrhenian Islands, based on a multi-locus analysis of genetic variation among 112 individuals from 28 populations [5]. However, to obtain a careful characterization of the respective distribution ranges, patterns of co-occurrence on single islands, relative abundance in sites of co-occurrence, as well as for all subsequent ecological and evolutionary investigations on these species, the characterization of a much larger sample of individuals and populations will be mandatory. Since nuclear and mitochondrial patterns of variation were largely congruent to one another in this group of species [5], the use of a DNA barcoding approach would be fully suitable.

We developed and validated a method for the assignment of a large number of individual mayflies to their putative species, based on the analysis of the high-resolution melting (HRM) curve of a standard molecular marker (NADH dehydrogenase subunit 1 gene). HRM is an effective, still not fully exploited, acquisition in the molecular taxonomy toolbox (see e.g. [10, 11]), which allows to assess the occurrence of sequence differences of diagnostic value between taxa, based on melting temperature profiles and the associated fluorescence peaks. Although based on a DNA barcoding approach, this method is faster, more cost-effective, and equally accurate when compared with traditional sequence-based approaches, provided that the pattern of variation within the focal taxa has already been characterized, as is the case for the B. rhodani group on the Tyrrhenian Islands [5].

Main text


In total, 399 larvae—100 individuals from a previous study [5] and 299 new samples-of the Baetis rhodani species group were collected from 59 localities (Table 1, Fig. 1) spanning the three main Tyrrhenian Islands (Sardinia, Corsica and Elba islands; no ethical approval is needed for sampling or processing these invertebrate species). The larvae were preserved in 95% ETOH prior to DNA extraction. A formal species description of the three species occurring in the study area is still pending. For consistency with the previous work which identified them as distinct biological species, we will refer to them here as Species A, Species B, and Species C. An unequivocal species assignment was made possible by the deep mitochondrial sequence divergence observed between the three species. Indeed, the percent sequence divergence between species (p-distance) was as follows [5]: A vs B = 14%, A vs C = 29%, and B vs C = 26%.

Table 1 Geographic location and sample size (n) of the 59 sites sampled within the Tyrrhenian Islands
Fig. 1

Geographic location of the 59 sites sampled within the Tyrrhenian Islands. Localities are numbered as in Table 1. Pie diagrams show the proportion of individuals belonging to the three species within each sampled locality. (The digital elevation model was downloaded by the WorldClim 1.4 database, freely available at

Total genomic DNA extraction was carried out from the thorax, legs, or entire larval body by proteinase K (Promega) digestion and a standard phenol–chloroform protocol [12]. The amount of genomic DNA of each sample was assessed with the Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA), using the Qubit dsDNA BR Assay kit (Thermo Fisher Scientific Inc.). Each sample was then concentrated to 10 ng/μL for downstream analyses.

Real-time quantitative PCR assays were carried out in 96-well plates through amplification of the mitochondrial NADH dehydrogenase subunit 1 gene (ND1) fragment with the following primers: ND1F (TAAAGTTAGCAGGTTCATACCC) and ND1R (CACCTATATTTGTACTTTGAAGG). Amplifications were run in a 15-μl reaction volume containing MgCl2 (2 mM, Promega), four dNTPs (0.2 mM each, Promega), two primers (0.2 μM each, Sigma-Aldrich), the enzyme Hotstart Taq polymerase (0.5 U, Promega), the reaction buffer (1×, Promega), EvaGreen® dye (1×, Biotium), and ~ 10–30 ng of DNA template.

Real time PCRs were run using a Roche LightCycler® LC480. The PCR cycling conditions were as follows: an initial step of 5 min at 95 °C, 35 cycles of 1 min at 94 °C, 45 s at 57 °C, and 90 s at 72 °C, and an extension step of 10 min at 72 °C. A final melting cycle was performed for 3 min at 95 °C and a melt from 40 to 95 °C collecting fluorescence continuously at a ramping rate of 0.1 °C per second. The HRM peaks obtained from single individuals were examined after each real-time PCR using the LightCycler® 480 Software (release version 1.5.0).

In order to characterize the HRM curve profiles of each species, we included in the analyses all the individuals belonging to the three species analysed in the previous study (n = 100; [5]), where the diagnostic value of sequence variation among the three species was assessed. To further test the accuracy of the proposed method, we randomly selected and sequenced 24 individuals from the pool of newly screened samples (Genebank Accession Numbers: MG581893–MG581916). We then compared their species assignment based on HRM profiles to those achieved through standard sequencing.

In order to evaluate the diagnostic ability of the HRM curve method for the three species studied, we performed a principal component analysis (PCA) using the IBM SPSS v. 23 (IBM SPSS Statistics for Windows, Armonk, NY: IBM Corp), based on peak temperature and fluorescence values of the 100 individuals from the previous paper [5]. The first principal component was then used to perform a receiver operating characteristic (ROC) curve analysis, whose accuracy was evaluated by means of the area under the curve (AUC) values.

Results and discussion

The derived HRM curves obtained for the individuals already sequenced from a previous study clearly showed three distinct profiles (see Fig. 2a). A comparison of these profiles with the standard sequence-based species identification unequivocally indicated the diagnostic value of the HRM curves. This facilitated assigning each unique profile to one of the three species, without misidentifications.

Fig. 2

High resolution melting curve analysis. a A negative first-derivative plot of fluorescence over temperature for three representative individuals of Species A, B, and C (as defined by [5]); redrawn from LightCycler 480 Software. b A 2D scatterplot of melting temperatures at fluorescence peaks for the 100 individuals used to assess the diagnostic value of melting curves for each species

The three HRM curves observed within the negative first-derivative plot were easily distinguishable from one another, based on curve shape and melting peak values (Fig. 2a). The curve identifying Species A showed two peaks, the first at 83 °C with a maximum value of 17.3 and a minimum value of 12.8 (−dF/dT), and the second peak at 86 °C with a maximum value of 20.5 and a minimum value of 8.6 (−dF/dT). The HRM curve for Species B was diagnosable by a single peak at 83–84 °C with a maximum value of 28.4 and a minimum value of 19.7 (−dF/dT). Finally, the HRM profile for Species C showed both the highest and narrowest curve profile, with a melting peak at 82 °C and a maximum value of 41.6 and a minimum value of 29.6 (−dF/dT). As shown by the scatterplot in Fig. 2b (see also Additional file 1), inter-individual variation within each species did not affect the diagnostic significance of the individual peak values.

The analysis of the 399 samples collected throughout the Tyrrhenian Islands (Sardinia, Corsica and Elba islands) did not reveal any new patterns for the HRM curve profiles in addition to the three described above. We assigned all analysed samples to a putative species: 128 to Species A, 170 to Species B, and 101 to Species C. As expected, based on the absence of new HRM curve variants among the newly studied samples, sequence-based assignment of all 24 individuals used as the control perfectly matched species identity as defined by melting profiles. Moreover, the ROC curve analysis based on the first principal component (variance explained: 82%) strongly supported the diagnostic ability of HRM method used, showing extremely high values of the AUC, in all pairwise comparisons between species (range of the AUC values: 0.97–1.00).

Our results clearly demonstrated that the HRM-based method proposed here to assess species identity for mayflies of the B. rhodani group on the Tyrrhenian Islands, proved as accurate and reliable as the standard sequence-based approach, besides allowing to avoid sequencing costs and efforts.

We utilized this method for a fast and accurate evaluation of the actual geographic distribution of the three cryptic species, which is fundamental information for any future study employing these species as target organisms. The results of this evaluation revealed distinct distribution patterns for the three studied species (Fig. 1), allowing us to incorporate them within a long-term study on the evolution of the insular biota (see also [13,14,15,16,17,18,19], to investigate species-specific eco-evolutionary dynamics, and to improve ongoing bio-monitoring programs in the area [20].


  • The method does not provide information about genetic diversity, population structure, or phylogenetic relationships of the study organisms.

  • The method is only applicable when background information on the interspecific genetic variation is available for diagnostic markers.





NADH dehydrogenase subunit 1 gene fragment


  1. 1.

    Brittain JE. Biology of mayflies. Annu Rev Entomol. 1982;27:119–97.

    Article  Google Scholar 

  2. 2.

    Williams HC, Ormerod SJ, Bruford MW. Molecular systematics and phylogeography of the cryptic species complex Baetis rhodani (Ephemeroptera, Baetidae). Mol Phylogenet Evol. 2006;40:370–82.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Lucentini L, Rebora M, Puletti ME, Gigliarelli L, Fontaneto D, Gaino E, et al. Geographical and seasonal evidence of cryptic diversity in the Baetis rhodani complex (Ephemeroptera, Baetidae) revealed by means of DNA taxonomy. Hydrobiologia. 2011;673:215–28.

    CAS  Article  Google Scholar 

  4. 4.

    Rutschmann S, Gattolliat JL, Hughes SJ, Baez M, Sartori M, Monaghan MT. Evolution and island endemism of morphologically cryptic Baetis and Cloeon species (Ephemeroptera, Baetidae) on the Canary Islands and Madeira. Freshw Biol. 2014.

    Google Scholar 

  5. 5.

    Bisconti R, Canestrelli D, Tenchini R, Belfiore C, Buffagni A, Nascetti G. Cryptic diversity and multiple origins of the widespread mayfly species group Baetis rhodani (Ephemeroptera: Baetidae) on northwestern Mediterranean islands. Ecol Evol. 2016;6(21):7901–10.

    Article  Google Scholar 

  6. 6.

    Müller-Liebenau I. Revision der europäischen Arten Gattung Baetis Leach, 1815 (Insecta, Ephemeroptera). Gewässer und Abwässer. 1969;48–49:1–214.

    Google Scholar 

  7. 7.

    Gattolliat JL, Sartori M. What is Baetis rhodani (Pictet 1843)(Insecta, Ephemeroptera, Baetidae)? Designation of a neotype and redescription of the species from its original area. Zootaxa. 2008;1957:69–80.

    Google Scholar 

  8. 8.

    Song H, Buhay JE, Whiting MF, Crandall KA. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci. 2008;105:13486–91.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Hudson RR, Coyne JA. Mathematical consequences of the genealogical species concept. Evolution. 2002;56:1557–65.

    Article  PubMed  Google Scholar 

  10. 10.

    Gopaul KK, Sells J, Lee R, Beckstrom-Sternberg SM, Foster JT, Whatmore AM. Development and assessment of multiplex high resolution melting assay as a tool for rapid single-tube identification of five Brucella species. BMC Res Notes. 2014;7(1):903.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Madesis P, Ganopoulos I, Anagnostis A, Tsaftaris A. The application of Bar-HRM (Barcode DNA-High Resolution Melting) analysis for authenticity testing and quantitative detection of bean crops (Leguminosae) without prior DNA purification. Food Control. 2012;25(2):576–82.

    CAS  Article  Google Scholar 

  12. 12.

    Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual. 2nd ed. New York: Cold Spring Harbor Laboratory Press; 1989.

    Google Scholar 

  13. 13.

    Salvi D, Schembri PJ, Sciberras A, Harris DJ. Evolutionary history of the Maltese wall lizard Podarcis filfolensis: insights on the ‘Expansion-Contraction’ model of Pleistocene biogeography. Mol Ecol. 2014;23:1167–87.

    Article  PubMed  Google Scholar 

  14. 14.

    Salvi D, Bisconti R, Canestrelli D. High phylogeographical complexity within Mediterranean islands: insights from the Corsican fire salamander. J Biogeogr. 2016;43(1):192–203.

    Article  Google Scholar 

  15. 15.

    Bisconti R, Canestrelli D, Nascetti G. Has living on islands been so simple? Insights from the insular endemic frog Discoglossus montalentii. PLoS ONE. 2013;8:e55735.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Bisconti R, Canestrelli D, Salvi D, Nascetti G. A geographic mosaic of evolutionary lineages within the insular endemic newt Euproctus montanus. Mol Ecol. 2013;22:143–56.

    Article  PubMed  Google Scholar 

  17. 17.

    Bauzà-Ribot MM, Jaume D, Fornos JJ, Juan C, Pons J. Islands beneath islands: phylogeography of a groundwater amphipod crustacean in the Balearic archipelago. BMC Evol Biol. 2011;11:211.

    Article  Google Scholar 

  18. 18.

    Bisconti R, Canestrelli D, Nascetti G. Multiple lines of evidence for demographic and range expansion of a temperate species Hyla sarda during the last glaciation. Mol Ecol. 2011;20:5313–27.

    Article  PubMed  Google Scholar 

  19. 19.

    Bisconti R, Canestrelli D, Nascetti G. Genetic diversity and evolutionary history of the Tyrrhenian treefrog Hyla sarda (Anura: Hylidae): adding pieces to the puzzle of Corsica-Sardinia biota. Biol J Linn Soc. 2011;103:159–67.

    Article  Google Scholar 

  20. 20.

    Buffagni A, Tenchini R, Cazzola M, Erba S, Balestrini R, Belfiore C, et al. Detecting the impact of bank and channel modification on invertebrate communities in Mediterranean temporary streams (Sardinia, SW Italy). Sci Total Environ. 2016;565:1138–50.

    CAS  Article  PubMed  Google Scholar 

Download references

Authors’ contributions

DC, CB, and GN conceived and designed the study. CB and RT carried out field sampling and species assignment. RB and RT performed the experiments. DC and RB analysed the data. DC and RB drafted the manuscript. All authors read and approved the manuscript.


We are grateful to Andrea Buffagni and Simone Cardoni for their help with sample collection.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Genbank Accession Numbers: MG581893–MG581916.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


This research was funded by the University of Tuscia.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Roberta Bisconti.

Additional file

Additional file 1.

An excel file format containing data used to draw the scatterplot of Fig. 2b.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bisconti, R., Tenchini, R., Belfiore, C. et al. Fast and accurate identification of cryptic and sympatric mayfly species of the Baetis rhodani group. BMC Res Notes 11, 7 (2018).

Download citation


  • High resolution melting (HRM) curve
  • Species identification
  • Baetis rhodani
  • Tyrrhenian islands
  • Real time PCR