Next Article in Journal
New Records of Powdery Mildews from Taiwan: Erysiphe ipomoeae comb. nov., E. aff. betae on Buckwheat, and E. neolycopersici comb. nov. on Cardiospermum halicacabum
Next Article in Special Issue
Similar Response of a Range Expanding Dragonfly to Low- and High-Elevation Predators
Previous Article in Journal
Response of Diatoms to the Changing Water Quality in the Myristica Swamps of the Western Ghats, India
Previous Article in Special Issue
Evaluating Potential Distribution and Niche Divergence among Populations of the World’s Largest Living Damselfly, Megaloprepus caerulatus (Drury, 1782)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DNA Barcoding and New Records of Odonates (Insecta: Odonata) from Paraíba State, Brazil

by
Ricardo Koroiva
1,2,*,
Vanessa Gabrielle Nóbrega Gomes
2 and
Diogo Silva Vilela
3
1
Departamento de Sistemática e Ecologia, Universidade Federal da Paraíba, João Pessoa 58051-900, Paraíba, Brazil
2
Núcleo de Biodiversidade, Instituto Nacional do Semiárido, Campina Grande 58434-700, Paraíba, Brazil
3
Departamento de Ciências Biológicas, Faculdade de Ciências e Letras de Assis, Universidade Estadual Paulista Júlio de Mesquita Filho, São Paulo 19806-900, Brazil
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(3), 203; https://doi.org/10.3390/d14030203
Submission received: 21 February 2022 / Revised: 7 March 2022 / Accepted: 8 March 2022 / Published: 10 March 2022
(This article belongs to the Special Issue Diversity, Ecology and Evolution of Odonata)

Abstract

:
Odonates (Insecta: Odonata) are important insects in the food chains of freshwater environments around the world, being used as a model species for areas of behavior and analysis of environmental quality. In Brazil, especially in the Northeastern region, both knowledge about the distribution and molecular information of odonate species found in the two main biomes of the region is still limited. Aiming to improve these issues, here, we carried out an Odonata survey in two locations and built a DNA barcode database for species from the state of Paraíba. In total, 15 first records were reported for this Brazilian state and 142 specimens from 27 genera and 45 species had their ‘Folmer’ cytochrome c oxidase subunit I (COI) fragment evaluated. The database we generated includes data for 70% of the Odonata species found in Paraíba state. For 16 species, this is the first DNA barcode available in public sequence repositories. Our results demonstrate that using the COI in the regional scale can help identify and delimit those evaluated. Eight species (17%) showed a low percentage of differentiation (<2%) compared to other species currently deposited in the GenBank or BOLD System; nevertheless, we present morphological traits that reaffirm our identifications. Barcode data provide new insights into Neotropical diversity and deliver basic information for taxonomic analyses.

1. Introduction

Odonates (Insecta: Odonata) are a group of fascinating insects that can be found all over the world (except for Antarctica) [1]. They have been used as model animals in several areas, such as behavioral studies and the analysis of environmental quality (e.g., [2,3,4]). This group plays an important role in the trophic network of freshwater environments, both as an efficient predator of invertebrates and as a prey for several vertebrates [5]. Over the last decade, scientists have increasingly discovered new species, which has improved the taxonomic knowledge of odonates in the Neotropical realm; however, information on their taxonomy (‘Linnean’ shortfall), distribution (‘Wallacean’ shortfall) and genetic diversity remains scarce (see [6,7]).
The traditional identification of odonates based on their phenotypic traits is often assumed to be difficult, considering that the older depictions are based on brief descriptions and rare illustrations (see examples in [8]). Even the most recent definitions present limited information to facilitate the identification of both larvae and adults. In general, in adults, anal appendages, wing venation and genitalia are often used to identify and classify odonates [9,10]. Despite the possibility of identifying some species at the larval level, the number of this type of descriptions is still scarce, especially in the Neotropical region [9].
Since Hebert and colleagues [11,12] suggested using cytochrome oxidase subunit I (COI) sequences as a global bioidentification system for animals in 2003, there has been great progress in using it for species identification and species discovery of odonates. This system, which is called DNA barcoding, has proven to be highly effective for delimiting and identifying different groups (e.g., [13,14]). As a result of an initiative by researchers at the University of Guelph (Ontario, Canada), a worldwide database has been maintained for the deposit and public identification of sequences; this database is called the BOLD system (the DNA Barcode of Life Data System) [15]. Among the advantages of this system compared with other public sequence repositories (e.g., GenBank) is that it is easy to access and download specimen data and sequences. For Neotropical odonates, the use of DNA barcoding is still limited, and there has been some discussion regarding its effectiveness by Koroiva et al. [16]. Moreover, Vilela et al. [17,18] have substantiated its direct application through integrative taxonomy, i.e., the framework to delimit and describe taxa by integrating information from multiple and complementary perspectives [19], using molecular and morphological data.
Currently, there are about 872 species of odonates in Brazil [20]. Several biomes can be found in this country, including the Atlantic Forest and ‘Caatinga’. The former is considered a biodiversity hotspot, and it is one of the most threatened tropical forests in the world [21]. Meanwhile, the latter is the only biome exclusive to Brazil. It is characterized by a semi-arid climate and has been recognized as an affected biome as a result of anthropic actions, such as free-living livestock and fuelwood extraction [22,23]. These two biomes are present in the state of Paraíba, a small territory state localized in Northeastern Brazil. The region’s list of Odonata species was recently published by Koroiva et al. [8]. In this work, the authors identified 49 species living in the region; this provides a good sampling of the diversity of the state, despite the few collections present in the Atlantic Forest. Nonetheless, it is essential to elaborate upon a broad and diverse species collection in order to realize representative genetic databases. In turn, using them to construct a reference database is the next step to enable the use of new molecular tools (e.g., DNA metabarcoding) in this region.
Keeping in mind the difficulties in the morphological identification of dragonflies, molecular tools present a promising way to solve this impediment in Neotropical species [16]. In this study, we present a DNA barcode library for the Odonata identification in the Paraíba state. This work aims (i) to improve the information on the species present in the Atlantic Forest of Paraíba; (ii) to establish DNA barcode libraries for the odonatofauna of Paraíba based on the COI gene and (iii) to evaluate the accuracy of the DNA barcodes in defining species in relation to both the regional and global databases.

2. Materials and Methods

2.1. Ethics Statement

This study was conducted with the appropriate permission (SISBIO license number 74324-5 and JBBM license number 002/2021/JBBM/SUDEMA).

2.2. Data Collection

Between May 2020 and November 2021, 740 specimens were collected from 10 municipalities in the state of Paraíba. Information about climate classification, precipitation, vegetation types and the geology of the sampling area are available in the work of Koroiva et al. [8]. In addition, two sites not considered by Koroiva et al. [8] were sampled during the 10 sampling campaigns between October 2020 and November 2021: João Pessoa Botanical Garden (JBBM), which is located in the municipality of João Pessoa (−7.135867, −34.860025; datum WGS84), and the “Banho do Jair” stream, found in the municipality of Santa Rita (−7.000965, −34.98836; datum WGS84). Notably, both areas are located in the Atlantic Forest fragments. The morphological identification of all specimens was done with the help of experts (see Acknowledgements) in Odonata taxonomy and by using the taxonomic keys of Lencioni [9,24,25] and Garrison et al. [10,26]. Our collection of specimens followed the methodology presented in Vilela et al. [27]. In terms of the classification, we followed Paulson et al. [28]. Voucher specimens were deposited into the Entomological Collection of the Department of Systematics and Ecology at the Federal University of Paraíba (DSEC/UFPB).

2.3. Extracting, Amplifying and Sequencing

All the DNA from the samples was extracted using the Blood & Tissue DNA Mini Kit (Ludwig Biotec, Alvorada, Brazil) from one leg, and it was preserved in ethanol. We amplified the genetic material of 44 specimens of 21 species collected in the two sampled sites mentioned above. We also used another 70 specimens of the 27 species previously collected by Koroiva et al. [8] and deposited in DSEC/UFPB (see Figure 1). For most species, we used specimens from different municipalities in an effort to analyze different populations (see BOLD dataset on http://dx.doi.org/10.5883/DS-ODOPB, accessed on 11 February 2022). In total, 658 bp were amplified from the 5′ region of the Cox1 gene using the M13-tailed primers OdoF1_t1 (5′–TGTAAAACGACGGCCAGTATTCAACHAATCATAARGATATTG G-3’) and OdoR1_t1 (5′–CAGGAAACAGCTATGACTAAACTTCTGGATGYCCRAARAAYCA-3’) (Semotok, unpublished, BOLD Systems http://www.boldsystems.org/index.php/Public_Primer_PrimerSearch, accessed on 27 January 2022). When it was not possible to amplify them, an approximately 421 bp long fragment at the 3′ end of the barcoding region was amplified by using the forward primer BF2 (5’-GCHCCHGAYATRGCHTTYCC-3’) and the reverse primer BR2 (5’-TCDGGRTGNCCRAARAAYCA-3’) described by Elbrecht and Leese [29]. This procedure was performed with consideration that the primers commonly used in COI amplification (OdoF1_t1-OdoR1_t1, HCO2198-LCO1490 [30], HCO2198_t1-LCO1490_t1 [31], LepF1-LepR1 [32] and LepF1_t1-LepR1_t1 [33]) were not successful in the amplification for many species, especially in Zygoptera, with exceptions of Telebasis corallina (Selys, 1876) and Hetaerina rosea Selys, 1853.
The PCR conditions for amplification consisted of 1 × buffer (Colorless GoTaq® Flexi Buffer; Promega Corp., Madison, WI, USA), 0.2 mM dNTP mix, 0.2 μM of each primer, 2 mM MgCl2, 1U Taq polymerase (GoTaq® G2 hot start polymerase, Promega Corp., Madison, WI, USA) and 2 μL of template DNA; these materials were placed in a total reaction volume of 25 μL. The PCR cycling program to OdoF1_t1 and OdoR1_t1 followed Vilela et al. [18]. For the BF2-BR2 primers, the PCR cycling program was run as follows: initial denaturation step with 3 min at 95 °C, 35 cycles of denaturation for 30 s at 95 °C, annealing for 45 s at 50 °C and extension for 1 min at 72 °C, and final extension for 5 min at 72 °C. The PCR products were purified with ethanol/sodium acetate and sequenced in an ABI 3130 Genetic Analyzer (Applied Biosystems, Waltham, MA, USA). The OdoF1_t1 and OdoR1_t1 sequencing were performed using M13 Forward (5’-TGTAAAACGACGGCCAGT-3’) and M13 Reverse primers (5’-CAGGAAACAGCTATGAC-3’), respectively. The sequence data were uploaded to GenBank (accession numbers OL806732 to OL806735 and OL806621 to OL806730) and were made available on BOLD as a dataset (http://dx.doi.org/10.5883/DS-ODOPB, accessed on 11 February 2022).

2.4. Data Analysis

To check the sequence quality of both strands and to assemble and edit them if necessary, we used GENEIOUS v 9.0.5 [34]. Furthermore, we aligned the sequences for each gene loci using Muscle v3.8.425 [35] (module implemented in GENEIOUS v 9.0.5) at the default setting. Five species that had less than three individuals were found to have sequenced in our database (i.e., singletons and doubletons; see on http://dx.doi.org/10.5883/DS-ODOPB, accessed on 11 February 2022). Using the sequences from the Brazilian specimens deposited in the BOLD system (accessed on 27 January 2022), we determined that our database had the incorporation of all sequences of these species (28 specimens); in turn, this allowed the analysis of intraspecific variation, totaling 142 specimens. The genetic distances between and within species were estimated using the Kimura’s two-parameter substitution model (K2P) (but see Srivathsan and Meier [36]); these were calculated using the MEGA X software [37]. To increase the robustness of the homology statement and to elevate the matrix occupancy, long sequences were truncated to cover only the ‘Folmer’ region of the COI gene. This is the most commonly used region for DNA barcoding, as it covers 658 nt of the 5′-end of the gene. For insects, the region can be amplified using the ‘Folmer’ primer pair (HCO2198 and LCO1490; [30]); subsequently, truncation was carried out following the positioning of these primers.
Next, we calculated the mean and maximum genetic divergence values and the lowest genetic distance on our (regional) database to the nearest neighbor in MEGA X [37]. We then plotted the empirical K2P values associated with intra- and interspecific comparisons against each other following the methods detailed by Koroiva and Kvist [38], in order to highlight and visualize any potential “barcode gap” (but see discussion on Wiemers and Fiedler [39]). To evaluate them in the global databases, we used the default settings of Web BLAST (Basic Local Alignment Search Tool; https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 11 February 2022) on GenBank [40] in order to identify the nearest matching sequences and the Ident and E(xpect) values. In addition, we used the Species Level Barcode Records option on BOLD (http://www.boldsystems.org/index.php/IDS_OpenIdEngine; accessed on 11 February 2022, [41]) to obtain a list of the sequences with the highest similarity. With consideration of previous works (e.g., [16]), we examined the species closest to those that showed less than 2% genetic divergence.

2.5. Species Delimitation

We used two methods to delimit the species: distance-based and tree-based methods. For a distance-based method, we performed ABGD (Automatic Barcode Gap Discovery analysis) online (http://wwwabi.snv.jussieu.fr/public/abgd/, accessed on 11 February 2022) [42] because the program is optimized for the COI gene; we used the values of Pmin = 0.001 and Pmax = 0.10, steps = 20, relative gap width (X) = 0.75 and Nb bins = 20 and K2P. The ABGD resulted in a stable genetic group count with a range of prior intraspecific values (P = 0.0113–0.0483) and the results of these grouping are presented. ASAP (Assemble Species by Automatic Partitioning) [43] analyses were performed on the website (https://bioinfo.mnhn.fr/abi/public/asap/, accessed on 11 February 2022) using K2P; we considered only the partition showing the lowest ASAP score.
For tree-based methods, we used the Poisson tree process (PTP) as implemented in the PTP (http://species.h-its.org/, accessed on 11 February 2022) using the default settings [44,45] and the generalized mixed Yule coalescent (GMYC) method. To run the PTP analysis, we first built a tree with RAxML (v 8.2.12) [46] using a GTR GAMMA I model and 1000 bootstrap replicates. The resulting tree was used as the input tree to run on the web server.
A single-threshold GMYC analysis was conducted in the species delimitation service (https://species.h-its.org/gmyc/, accessed on 11 February 2022) [47]. An ultrametric single-locus gene tree was obtained using BEAST v.2.6.6 [48] with 1.5 × 108 MCMC generations under relaxed lognormal clock, the Yule process tree model, and a burn-in of the first 10% generations of the final consensus tree. The posterior distributions (ESS > 200) were examined in Tracer v. 1.6 [49]. The best fitting available model was identified through the jModelTest v. 2.1.7 (AIC & BIC, GTR + I + G) [50,51].
Finally, we created a neighbor-joining (NJ) dendrogram to provide a graphic representation of the divergence pattern between species. An NJ tree was inferred using MEGA X software [37]. In all the tree-based analysis, three species were used as the out group (HQ941355 Baetis adonis, HQ943539 Baetis phoebus and HQ987969 Ephemerella mucronata). It should be noted that the tree presented here is only intended to represent the distance matrix and it should not be interpreted as a phylogenetic hypothesis.

3. Results

3.1. Sampling

Our samples added 15 new species records to the odonatofauna of Paraíba: Dasythemis venosa (Burmeister, 1839), Dythemis nigra Martin, 1897, Epipleoneura metallica Rácenis, 1955, Erythrodiplax cf. fervida (Erichson in Schomburgk, 1848), Idioneura ancilla Selys, 1860, Macrothemis imitans Karsch, 1890, Metaleptobasis bicornis (Selys, 1877), Micrathyria didyma (Selys in Sagra, 1857), Micrathyria mengeri Ris, 1919, Nephepeltia berlai Santos, 1950, Orthemis flavopicta Kirby, 1889, Perithemis lais (Perty, 1834), Tauriphila australis (Hagen, 1867), Telebasis griffinii (Martin, 1896) and Triacanthagyna septima (Selys in Sagra, 1857).
After performing the morphological analysis, all specimens (including those deposited in the DSEC/UFPB) that had previously been identified as Anatya guttata (Erichson in Schomburgk, 1848) were now identified as Anatya januaria Ris, 1911.

3.2. Genetic Variation

A total of 142 mitochondrial COI barcode sequences were obtained from five families, 27 genera and 45 species (mean by species 3.15; max = 15, min = 1, Table 1, Figure 2). All of the analyzed sequences were larger than 353 bp. The average base pairs of the sequences were 524 bp (SD = 116.54) and the median was 592 bp. Thirteen singletons were registered in the database (Table 1). In the regional database (Table 1), the greatest intraspecific variation was the Erythrodiplax fusca (Rambur, 1842) with 1.85% and the smallest inter-specific variation occurred between the Erythrodiplax leticia Machado, 1996 and the Erythrodiplax cf. fervida, with 4.87%.
Figure 3 shows the clear separation of intraspecific and interspecific distances and the so-called “barcoding gap” on the regional DNA barcode library. Considering the global databases, eight species showed close proximity (<2%) to records of other deposited species (Table 1): A. januaria, Erythemis carmelita Williamson, 1923, Erythrodiplax basalis (Kirby, 1897), H. rosea, N. berlai, Perithemis tenera (Say, 1840), Telebasis filiola (Perty, 1834) and Tramea cophysa Hagen, 1867.

3.3. Species Delimitation

Species delimitation analyses provided 45 species defined by morphological delimitation. Using molecular methods to delimit species, the results were 45 for both the ABGD-initial partition and ABGD-recursive partition, 39 for ASAP, 46 for PTP and GMYC (Figure 4).

4. Discussion

This study adds 15 Odonata species to the state of Paraíba. As a result, this state now has 64 species recorded; thus, it is third in number of species in the Northeast region of Brazil, behind only Bahia and Ceará (174 species [52] and 73 species [53], respectively). This library of DNA sequences is the first publication of DNA barcoding for the 16 Odonata species in public sequence repositories: Acanthagrion gracile (Rambur, 1842), A. januaria, Ep. metallica, E. avittata Borror, 1942, E. leticia, E. cf. fervida, E. basalis, E. cf. unimaculata (De Geer, 1773), M. griseofrons Calvert, 1909, Mecistogaster kesselringi Soldati and Machado, 2019, Metaleptobasis bicornis (Selys, 1877), N. berlai, O. flavopicta, Progomphus dorsopallidus Byers, 1934, T. filiola and Zenithoptera lanei Santos, 1941. The database generated for the present study provides data for 70% (45 species) of the odonate species found in Paraíba state. Despite all these numbers, we recognize that this work is a starting point for the study of odonates in Paraíba state, considering that the likely diversity of the region must be far greater than our estimates.
It is important to highlight the methodological limitations that we faced for the amplification of Neotropical odonates. Although many studies indicate the effectiveness of the tested primers on several continents (e.g., [54,55]), they did not have the amplification capacity for all the species amplified in this study, most notably in Zygoptera. Problems with amplifying Neotropical species using primers commonly used in the world is not new, and it is not exclusive to odonates (see [13] for frogs and [14] for fishes). Jennings and collaborators [14] highlight the importance of this type of information considering that the trial and error nature for the choice of primers wastes labor and reagents.
In addition to demonstrating that the two pairs of primers used in this study are capable of amplifying the COI fragment in the two suborders present in the Neotropical region, this study shows that it is also important to consider the usefulness of this marker to discriminate between different species in terms of their DNA barcode and metabarcoding studies. The regional database does not present the overlapping of inter- and intraspecific genetic variation; however, in the global database analysis (using Genbank and the BOLD System), eight species showed close genetic proximity to the species we examined (<2%). Below, we discuss the most likely hypothesis for each taxon as well as the disagreements we noticed in the global dataset evaluation.
Our results showed that the analyzed specimens of A. januaria from Paraíba state are 99.74% similar to the available sequences for A. guttata (Erichson, 1848). Problems with determining the species of Anatya are largely recognized, mainly due to some specific variable characteristics being wrongly interpreted (see [10], p. 224). Our specimens had subtle differences among them, including the size of cerci and body length; however, based on the comparison of these structures drawn by Ris ([56], p. 424) and Garrison et al. [10], the shape of the posterior hamule left almost no doubt that our specimens belong to A. januaria (Figure 5A).
In turn, the E. carmelita we analyzed was 98.67% similar to a sequence identified as the species E. mithroides. We are confident with regard to our morphological identification because, despite their color resemblance, these two species are very easily distinguishable because of their abdominal characteristics. Erythemis carmelita belongs to the group of Erythemis species that possesses greatly swollen basal abdominal segments (which are slightly swollen in E. mithroides) and narrow remainder segments (which are broad in E. mithroides) ([10], p. 240–241).
In addition, we found a high similarity between the sequences of our E. basalis and the deposited sequences of E. paraguayensis (Förster, 1905). Borror [57] states that specimens of E. paraguayensis are “apt to be confused with small individuals of basalis” (p. 153); however, morphologically, these species are difficult to be confused with one another because of multiple characteristics: E. basalis belongs to Borror’s basalis group, which consists of species that have hamules that are usually slender and have an outer branch a little longer than their inner branch (evident in lateral view, Figure 5B), a terminal segment of a slender penis, small and rounded lateral lobes (Figure 5C)—characteristics that we can easily observe in E. basalis specimens [57].
In contrast, E. paraguayensis belongs to the connata group, in which males present hamules that are moderately robust, their outer branch is almost equal in size to inner branch, and the median process of the penis is very prominent and slender, similar to what we can observe in the males of E. paraguayensis. In addition, E. paraguayensis is one of the smallest species of the genus and the extent of the small spots of the hind wings nearly trespasses upon the cubital space. In contrast, E. basalis is a much larger species, its hind wing spots reaching (or in some cases, trespassing upon) the first antenodal vein. Based on these characters, we argue that our morphological identification is correct.
Although our sequences of H. rosea are 98.95% similar to the ones assigned to H. sanguinea Selys, 1853, those two species are highly unlikely to be confused with one another. This is mainly because of the paraprocts (Figure 5D) are long and well developed in H. rosea, while they are vestigial in H. sanguinea (see [58] for a revision). Another character that can easily separate the two species is the median lobe of cercus. It is bilobed in H. rosea and entirely so in H. sanguinea (see [24], pp. 66–67).
Nephepeltia berlai shares several morphological characters with N. aequisetis Calvert, 1909, with the most prominent among them being the tubercle on the venter of thorax, the length and placement of the hind tibiae spurs, and also the vesica spermalis (see [59] for a revision). However, there are some characters on the vesica spermalis (medio-ectal distal process) and the cercus (level of the distal end of ventral toothed carina at about distal fourth of cercus length) that allow for a separation of the two species. Based on those characters, we believe that our specimens are N. berlai.
We based our identifications of the genus Perithemis largely on the study of von Ellenrieder and Muzón [60], which was the first study to separate the species of this genus using characteristics other than coloration and wing venation. The males we identified as P. tenera (Say, 1840) (regarded as P. mooma Kirby, 1889 in [60]) present wings that are uniformly colored (as in P. icteroptera Selys in Sagra, 1857), have a tip of hamuli (Figure 5E) at least 0.40 of the ventral margin (at the level of ventral margin in P. icteroptera), and a penis with a first segment trapezoidal (rounded in P. icteroptera). Our sequences were 99.83% similar to the sequences assigned to P. icteroptera; however, the characteristics described above led us to identify our specimens as P. tenera.
The sequences of our specimens identified as T. filiola were 98.96% identical to the specimens deposited assigned as T. willinki Fraser, 1948. These species are very similar morphologically, as it is stated in the revision of the genus [61]. However, despite the great resemblance of these taxa, Garrison [61] properly diagnosed the two species, showing that the cercus of T. filiola (Figure 5F) is distinctly shorter than the paraproct (it is subequal to the paraproct in T. willinki). In turn, we followed the same diagnosis to assign our specimens as T. filiola.
Lastly, we identified some males of T. cophysa showing agreement with the diagnosis presented in DeMarmels and Rácenis [62]. However, our sequences were 98.34% similar to the sequences assigned to T. binotata (Rambur, 1842) in the BOLD System. These two taxa, although they belong to the same genus, are quite different. Due to their color and morphological differences, they were placed in different groups within Tramea. The cophysa group, to which T. cophysa belongs, is composed of four species that share “only one constant characteristic common to all four species and, peculiar to the “cophysa-group”, are two oblique pale lateral bands on the synthorax” [62]. Such a characteristic is present in our specimens while it is absent in T. binotata, as this is a species with an overall blue-grey coloration, contrasting with the reddish coloration of the species of the cophysa group. Additionally, as a result of our comparison of morphological features, such as the shape of the posterior hamule and the length and shape of the cercus (Figure 5G), we were able to make a safe distinction between our specimens to other taxa within Tramea.
The molecular species delimitation results were identical to the morphological results for most of the species we examined. The presence of divergent results is commonly used to indicate possible cryptic species (e.g., [63]); in turn our identical results in terms of morphological, ABGD, and the occasional divergences suggest that COI has the ability to delimit in the evaluated species.
In summary, our results demonstrate that DNA barcoding can be used to delimit and differentiate odonates on a regional scale. Of the 45 species evaluated in this study, only eight species (17%) showed any disagreement with the global databases and all species (100%) could be identified when we consider only our regional DNA barcode database. Considering the issues with the global databases, our results for the number of species that can be readily identified using their DNA barcoding (83%) are close to those found in other regions of the world. In a genetic database for the Central and North European odonates, the effectiveness was 88% in a set of 103 species [64]. Values between 79% and 89% were also found in datasets in countries such as the Philippines (in a set of 38 species; [54]), Italy (in a set of 88 species; [65]) and Malta (in a set of 10 species, [55]). In Brazil, the only evaluation (in a set of 38 species) performed indicated a success rate between 79% and 94% depending on the analysis criteria [16].
Subsequently, all these results indicate the importance of performing careful morphological analysis (see this kind of problem in [66]). Moreover, as indicated by Koroiva and Kvist [38], it reinforces that the use of the global database does not allow the correct establishment of all molecular identifications to be correct. Many reasons justify the discrepancy between the morphological and molecular results. Among in odonates, three main causes are the presence of cryptic species, rapid and/or recent radiation events and errors in identifying the deposited specimens [38]. Regarding this last issue, the description and presentation of the key structures for identification presented above aimed to facilitate future comparisons in this sense.

5. Conclusions

The establishment of an Odonata DNA barcode library for the Paraíba state is a milestone that will improve the taxonomy and biodiversity conservation for Neotropical species. Despite the difficulties of using traditional primers for amplifying the Neotropical species, our results demonstrate that using the COI in the regional scale can help identify and delimit those evaluated. Our results for the number of species that can be readily identified using their DNA barcoding (83%) are close to the results found in other regions of the world. Keeping in mind the problems of using public genetic databases for identification, here, we present morphological evidence for our identifications in the cases of disagreement. In turn, this facilitates comparisons and allows for new questions to arise about the genetic diversity of tropical species.

Author Contributions

Conceptualization, R.K.; methodology, R.K., V.G.N.G. and D.S.V.; formal analysis, R.K., V.G.N.G. and D.S.V.; data curation, R.K. and D.S.V.; writing—original draft preparation, R.K., V.G.N.G. and D.S.V.; writing—review and editing, R.K., V.G.N.G. and D.S.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by PROPESQ/PRPG/UFPB, grant number PVA13280-2020 and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES), finance Code 001. DSV thanks Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) for a postdoctoral fellowship grant (Proc. 2019/26438-9).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are contained within the article; specimens analyzed in this study are deposited in Entomological Collection of the Department of Systematics and Ecology of the Federal University of Paraíba (DSEC/UFPB) and are available on request to the collection managers. The sequences are available at GenBank (accession numbers OL806732 to OL806735 and OL806621 to OL806730) and BOLD system (http://dx.doi.org/10.5883/DS-ODOPB, accessed on 11 February 2022).

Acknowledgments

We are thankful to Jardim Botânico de João Pessoa (JBBM/SUDEMA) for the assistance in sampling. We are grateful to Frederico Augusto de Atayde Lencioni and Rosser W. Garrison for confirmation of species identifications. We are also grateful to the academic editor and the anonymous reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Trueman, J.W.H. A brief history of the classification and nomenclature of Odonata. Zootaxa 2007, 1668, 381–394. [Google Scholar] [CrossRef]
  2. Cordoba-Aguilar, A. Dragonflies and Damselflies: Model Organisms for Ecological and Evolutionary Research, 1st ed.; Oxford University Press: Oxford, UK, 2009. [Google Scholar]
  3. Miguel, T.B.; Oliveira-Junior, J.M.B.; Ligeiro, R.; Juen, L. Odonata (Insecta) as a tool for the biomonitoring of environmental quality. Ecol. Indic. 2017, 81, 555–566. [Google Scholar] [CrossRef]
  4. Veras, D.S.; Lustosa, G.S.; Moura, L.P.; Ferreira, M.F.R.; Juen, L. Differences in land use modify Odonata assemblages in the Cerrado-Caatinga ecotone. Acta Limnol. Bras. 2020, 32, e15. [Google Scholar] [CrossRef]
  5. May, M. Odonata: Who They Are and What They Have Done for Us Lately: Classification and Ecosystem Services of Dragonflies. Insects 2019, 10, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bried, J.; Ries, L.; Smith, B.; Patten, M.; Abbott, J.; Ball-Damerow, J.; Cannings, R.; Cordero-Rivera, A.; Córdoba-Aguilar, A.; De Marco, P.; et al. Towards Global Volunteer Monitoring of Odonate Abundance. Bioscience 2020, 70, 914–923. [Google Scholar] [CrossRef]
  7. Gómez-Tolosa, M.; Rivera-Velázquez, G.; Rioja-Paradela, T.M.; Mendoza-Cuenca, L.F.; Tejeda-Cruz, C.; López, S. The use of Odonata species for environmental assessment: A meta-analysis for the Neotropical region. Environ. Sci. Pollut. Res. 2021, 28, 1381–1396. [Google Scholar] [CrossRef] [PubMed]
  8. Koroiva, R.; Pereira-Colavite, A.; Batista, F.R.D.C.; Vilela, D.S. Checklist and contribution to the knowledge of the odonatofauna of Paraíba state, Brazil. Biota Neotrop. 2021, 21, e20211196. [Google Scholar] [CrossRef]
  9. Lencioni, F.A.A. Damselflies of Brazil—An Illustrated Identification Guide—Southeast Region; F.A.A. Lencioni: Jacareí, Brazil, 2017. [Google Scholar]
  10. Garrison, R.W.; von Ellenrieder, N.; Louton, J.A. Dragonfly Genera of the New Word: An Illustrated and Annotated Key to the Anisoptera. Available online: http://entomologia.net/L_Odonata/Odonata_Dragonfly_Genera_of_the_New_World__An_Illustrated_and_Annotated_Key_to_the_Anisoptera.pdf (accessed on 5 March 2016).
  11. Hebert, P.D.N.; Cywinska, A.; Ball, S.L.; de Waard, J.R. Biological identifications through DNA barcodes. Proc. Biol. Sci. 2003, 270, 313–321. [Google Scholar] [CrossRef] [Green Version]
  12. Hebert, P.D.N.; Ratnasingham, S.; de Waard, J.R. Barcoding animal life: Cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. B Biol. Sci. 2003, 270, S96–S99. [Google Scholar] [CrossRef] [Green Version]
  13. Lyra, M.L.; Haddad, C.F.B.; de Azeredo-Espin, A.M.L. Meeting the challenge of DNA barcoding Neotropical amphibians: Polymerase chain reaction optimization and new COI primers. Mol. Ecol. Resour. 2017, 17, 966–980. [Google Scholar] [CrossRef]
  14. Jennings, W.B.; Ruschi, P.A.; Ferraro, G.; Quijada, C.C.; Silva-Malanski, A.C.G.; Prosdocimi, F.; Buckup, P.A. Barcoding the Neotropical freshwater fish fauna using a new pair of universal COI primers with a discussion of primer dimers and M13 primer tails. Genome 2019, 62, 77–83. [Google Scholar] [CrossRef] [PubMed]
  15. Ratnasingham, S.; Hebert, P.D.N. BOLD: The Barcode of Life Data System (www.barcodinglife.org). Mol. Ecol. Notes 2007, 7, 355–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Koroiva, R.; Pepinelli, M.; Rodrigues, M.E.; Roque, F.O.; Lorenz-Lemke, A.P.; Kvist, S. DNA barcoding of odonates from the Upper Plata basin: Database creation and genetic diversity estimation. PLoS ONE 2017, 8, e0182283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Vilela, D.S.; Koroiva, R.; Cordero-Rivera, A.; Guillermo-Ferreira, R. A further study on Franciscobasis Machado & Bedê, 2016 (Odonata: Coenagrionidae), a newly described genus from Minas Gerais, Brazil. PLoS ONE 2019, 14, e0223241. [Google Scholar] [CrossRef]
  18. Vilela, D.S.; Anjos-Santos, D.; Koroiva, R.; Cordero-Rivera, A.; Guillermo-Ferreira, R. Revision of the genus Minagrion Santos, 1965 (Odonata: Coenagrionidae). Zootaxa 2020, 4786, 176–198. [Google Scholar] [CrossRef] [PubMed]
  19. Pante, E.; Schoelinck, C.; Puillandre, N. From integrative taxonomy to species description: One step beyond. Syst. Biol. 2015, 64, 152–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. IUCN. The IUCN Red List of Threatened Species: Version 2021-3. Available online: https://www.iucnredlist.org/search?permalink=7d7c85cb-009c-4c4d-83ae-25e4675dcd32 (accessed on 18 February 2022).
  21. Laurance, W.F. Conserving the hottest of the hotspots. Biol. Conserv. 2009, 142, 1137. [Google Scholar] [CrossRef]
  22. Marinho, F.P.; Mazzochini, G.G.; Manhães, A.P.; Weisser, W.W.; Ganade, G. Effects of past and present land use on vegetation cover and regeneration in a tropical dryland forest. J. Arid Environ. 2016, 132, 26–33. [Google Scholar] [CrossRef]
  23. da Silva, J.M.C.; Barbosa, L.C.F. Impact of Human Activities on the Caatinga. In Caatinga; Silva, J.M.C., da Tabarelli, M., Leal, I.R., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 359–368. [Google Scholar]
  24. Lencioni, F.A.A. Damselflies of Brazil: An Illustrated Identification Guide, Volume 1: Non-Coenagrionidae Families; All Print Editora: São Paulo, Brazil, 2005; ISBN 8598310565. [Google Scholar]
  25. Lencioni, F.A.A. Damselflies of Brazil: An Illustrated Identification Guide, Volume 2: Coenagrionidae; All Print Editora: São Paulo, Brazil, 2006; ISBN 8577180344. [Google Scholar]
  26. Garrison, R.W.; von Ellenrieder, N.; Louton, J.A. Damselfly Genera of the New World: An Illustrated and Annotated Key to the Zygoptera; Johns Hopkins University Press: Baltimore, MD, USA, 2010; ISBN 0801896703. [Google Scholar]
  27. Vilela, D.S.; Koroiva, R.; Tosta, T.H.A.; Novaes, M.C.; Guillermo-Ferreira, R. Dragonflies and damselflies from the West of Minas Gerais, Brazil: Checklist and new records. Biota Neotrop. 2020, 20, e20190851. [Google Scholar] [CrossRef] [Green Version]
  28. Paulson, D.; Schorr, M.; Deliry, C. World Odonata List. Available online: https://www2.pugetsound.edu/academics/academic-resources/slater-museum/biodiversity-resources/dragonflies/world-odonata-list2/ (accessed on 11 February 2021).
  29. Elbrecht, V.; Leese, F. Validation and Development of COI Metabarcoding Primers for Freshwater Macroinvertebrate Bioassessment. Front. Environ. Sci. 2017, 5, 11. [Google Scholar] [CrossRef] [Green Version]
  30. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar] [PubMed]
  31. Foottit, R.G.; Maw, H.E.L.; Havill, N.P.; Ahern, R.G.; Montgomery, M.E. DNA barcodes to identify species and explore diversity in the Adelgidae (Insecta: Hemiptera: Aphidoidea). Mol. Ecol. Resour. 2009, 9, 188–195. [Google Scholar] [CrossRef] [PubMed]
  32. Hebert, P.D.N.; Penton, E.H.; Burns, J.M.; Janzen, D.H.; Hallwachs, W. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc. Natl. Acad. Sci. USA 2004, 101, 14812–14817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhang, J.; Hanner, R. Molecular Approach to the Identification of Fish in the South China Sea. PLoS ONE 2012, 7, e30621. [Google Scholar] [CrossRef] [PubMed]
  34. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  35. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  36. Srivathsan, A.; Meier, R. On the inappropriate use of Kimura-2-parameter (K2P) divergences in the DNA-barcoding literature. Cladistics 2012, 28, 190–194. [Google Scholar] [CrossRef] [PubMed]
  37. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  38. Koroiva, R.; Kvist, S. Estimating the barcoding gap in a global dataset of cox1 sequences for Odonata: Close, but no cigar. Mitochondrial DNA Part A 2018, 29, 765–771. [Google Scholar] [CrossRef]
  39. Wiemers, M.; Fiedler, K. Does the DNA barcoding gap exist?—A case study in blue butterflies (Lepidoptera: Lycaenidae). Front. Zool. 2007, 4, 8. [Google Scholar] [CrossRef] [Green Version]
  40. Benson, D.A.; Cavanaugh, M.; Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. Nucleic Acids Res. 2013, 41, D36–D42. [Google Scholar] [CrossRef] [Green Version]
  41. Ratnasingham, S.; Hebert, P.D.N. A DNA-Based Registry for All Animal Species: The Barcode Index Number (BIN) System. PLoS ONE 2013, 8, e66213. [Google Scholar] [CrossRef] [Green Version]
  42. Puillandre, N.; Lambert, A.; Brouillet, S.; Achaz, G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 2012, 21, 1864–1877. [Google Scholar] [CrossRef]
  43. Puillandre, N.; Brouillet, S.; Achaz, G. ASAP: Assemble species by automatic partitioning. Mol. Ecol. Resour. 2021, 21, 609–620. [Google Scholar] [CrossRef]
  44. Zhang, J.; Kapli, P.; Pavlidis, P.; Stamatakis, A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics 2013, 29, 2869–2876. [Google Scholar] [CrossRef] [Green Version]
  45. Kapli, P.; Lutteropp, S.; Zhang, J.; Kobert, K.; Pavlidis, P.; Stamatakis, A.; Flouri, T. Multi-rate Poisson Tree Processes for single-locus species delimitation under Maximum Likelihood and Markov Chain Monte Carlo. Bioinformatics 2017, 33, 1630–1638. [Google Scholar] [CrossRef] [Green Version]
  46. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  47. Fujisawa, T.; Barraclough, T.G. Delimiting species using single-locus data and the generalized mixed yule coalescent approach: A revised method and evaluation on simulated data sets. Syst. Biol. 2013, 62, 707–724. [Google Scholar] [CrossRef] [Green Version]
  48. Bouckaert, R.; Vaughan, T.G.; Barido-Sottani, J.; Duchêne, S.; Fourment, M.; Gavryushkina, A.; Heled, J.; Jones, G.; Kühnert, D.; De Maio, N.; et al. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 2019, 15, e1006650. [Google Scholar] [CrossRef] [Green Version]
  49. Rambaut, A.; Suchard, M.A.; Xie, D.; Drummond, A.J. Tracer, v.1.6. Available online: http://beast.bio.ed.ac.uk/Tracer (accessed on 6 May 2016).
  50. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [Green Version]
  51. Guindon, S.; Gascuel, O. A Simple, Fast, and Accurate Algorithm to Estimate Large Phylogenies by Maximum Likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Ribeiro, C.; Firme, B.; Araujo, S.A.; de Sá, A.; Zander, F.; Teixeira, K.; Santos, L.R.; Rodrigues, M.E. Check-list of Odonata from the state of Bahia, Brazil: Ecological information, distribution, and new state records. Odonatologica 2021, 50, 161–186. [Google Scholar] [CrossRef]
  53. Takiya, D.; Santos, A.P.; Pinto, Â.; Henriques-Oliveira, A.L.; Carvalho, A.; Sampaio, B.; Clarkson, B.; Moreira, F.; Avelino-Capistrano, F.; Gonçalves, I.; et al. Aquatic Insects from the Caatinga: Checklists and diversity assessments of Ubajara (Ceará State) and Sete Cidades (Piauí State) National Parks, Northeastern Brazil. Biodivers. Data J. 2016, 4, e8354. [Google Scholar] [CrossRef] [PubMed]
  54. Casas, P.A.S.; Sing, K.-W.; Lee, P.-S.; Nuñeza, O.M.; Villanueva, R.J.T.; Wilson, J.-J. DNA barcodes for dragonflies and damselflies (Odonata) of Mindanao, Philippines. Mitochondrial DNA Part A 2018, 29, 206–211. [Google Scholar] [CrossRef] [PubMed]
  55. Rewicz, T.; Móra, A.; Tończyk, G.; Szymczak, A.; Grabowski, M.; Calleja, E.J.; Pernecker, B.; Csabai, Z. First records raise questions: DNA barcoding of Odonata in the middle of the Mediterranean. Genome 2021, 64, 196–206. [Google Scholar] [CrossRef]
  56. Ris, F. Libellulinen. In Collections Zoologiques du Baron Edm. De Selys Longchamps—Catalogue Systématique et Descriptif; Hayez: Bruxelles, Belgium, 1911; pp. 385–528. [Google Scholar]
  57. Borror, D.J. A revision of the Libelluline genus Erythrodiplax (Odonata); The Ohio State University: Columbus, OH, USA, 1942. [Google Scholar]
  58. Garrison, R.W. A synopsis of the genus Hetaerina with descriptions of four new species (Odonata: Calopterygidae). Trans. Am. Entomol. Soc. 1990, 116, 175–259. [Google Scholar]
  59. von Ellenrieder, N. A synopsis of the Neotropical genus Nephepeltia (Odonata: Libellulidae), including description of a new species, synonymies, and a key to males. Zootaxa 2014, 3796, 121–146. [Google Scholar] [CrossRef] [Green Version]
  60. von Ellenrieder, N.; Muzón, J. The Argentinean species of the genus Perithemis Hagen (Anisoptera: Libellulidae). Odonatologica 1999, 28, 1–14. [Google Scholar]
  61. Garrison, R.W. A synopsis of the genus Telebasis (Odonata: Coenagrionidae). Int. J. Odonatol. 2009, 12, 1–121. [Google Scholar] [CrossRef]
  62. De Marmels, J.; Rácenis, J. An analysis of the cophysa-group of Tramea Hagen, with descriptions of two new species (Anisoptera: Libellulidae). Odonatologica 1982, 11, 109–128. [Google Scholar]
  63. Koroiva, R.; Rodrigues, L.R.R.; Santana, D.J. DNA barcoding for identification of anuran species in the central region of South America. PeerJ 2020, 8. [Google Scholar] [CrossRef] [PubMed]
  64. Geiger, M.; Koblmüller, S.; Assandri, G.; Chovanec, A.; Ekrem, T.; Fischer, I.; Galimberti, A.; Grabowski, M.; Haring, E.; Hausmann, A.; et al. Coverage and quality of DNA barcode references for Central and Northern European Odonata. PeerJ 2021, 9, e11192. [Google Scholar] [CrossRef] [PubMed]
  65. Galimberti, A.; Assandri, G.; Maggioni, D.; Ramazzotti, F.; Baroni, D.; Bazzi, G.; Chiandetti, I.; Corso, A.; Ferri, V.; Galuppi, M.; et al. Italian odonates in the Pandora’s box: A comprehensive DNA barcoding inventory shows taxonomic warnings at the Holarctic scale. Mol. Ecol. Resour. 2021, 21, 183–200. [Google Scholar] [CrossRef] [PubMed]
  66. Vega-Sánchez, Y.M.; Lorenzo-Carballa, M.O.; Vilela, D.S.; Guillermo-Ferreira, R.; Koroiva, R. Comment on “Molecular identification of seven new Zygopteran genera from South China through partial cytochrome oxidase subunit I (COI) gene”. Meta Gene 2020, 25, 100759. [Google Scholar] [CrossRef]
Figure 1. Geographical location of sampled specimens of odonates. (A) Map of South America (dark grey) highlighting the geopolitical division of Brazil (white); (B) municipality division of Paraíba with the locality of our two sampling sites (black dots) and the other localities for specimens obtained in DSEC/UFPB (grey dots).
Figure 1. Geographical location of sampled specimens of odonates. (A) Map of South America (dark grey) highlighting the geopolitical division of Brazil (white); (B) municipality division of Paraíba with the locality of our two sampling sites (black dots) and the other localities for specimens obtained in DSEC/UFPB (grey dots).
Diversity 14 00203 g001
Figure 2. Examples of odonates (Insecta: Odonata) collected and sequenced from Paraíba State, Brazil. (A) Anatya januaria, (B) Perithemis tenera, (C) Tramea cophysa; (D) Hetaerina rosea, (E) Telebasis filiola, (F) Lestes forficula; (G) Micrathyria hesperis, (H) Erythrodiplax basalis and (I) Erythemis plebeja.
Figure 2. Examples of odonates (Insecta: Odonata) collected and sequenced from Paraíba State, Brazil. (A) Anatya januaria, (B) Perithemis tenera, (C) Tramea cophysa; (D) Hetaerina rosea, (E) Telebasis filiola, (F) Lestes forficula; (G) Micrathyria hesperis, (H) Erythrodiplax basalis and (I) Erythemis plebeja.
Diversity 14 00203 g002
Figure 3. Frequency distribution of intraspecific (black) and interspecific (grey) genetic divergence in the sampled odonates (regional database).
Figure 3. Frequency distribution of intraspecific (black) and interspecific (grey) genetic divergence in the sampled odonates (regional database).
Diversity 14 00203 g003
Figure 4. Molecular species delimitation of odonates from the Paraíba state, Brazil, based on DNA barcodes. Black bars indicate congruent results between the molecular and morphological identifications, and grey bars indicate divergence results from morphological identification. It should be noted that the neighbor-joining tree presented here is only intended to represent the distance matrix, and it should not be interpreted as a phylogenetic hypothesis.
Figure 4. Molecular species delimitation of odonates from the Paraíba state, Brazil, based on DNA barcodes. Black bars indicate congruent results between the molecular and morphological identifications, and grey bars indicate divergence results from morphological identification. It should be noted that the neighbor-joining tree presented here is only intended to represent the distance matrix, and it should not be interpreted as a phylogenetic hypothesis.
Diversity 14 00203 g004
Figure 5. Morphology of male Odonata adults collected in Paraíba state, Brazil (lateral view of): (A) secondary genitalia of Anatya januaria; (B) secondary genitalia and (C) vesica spermalis of Erythrodiplax basalis; (D) caudal appendages of Hetaerina rosea; (E) secondary genitalia of Perithemis tenera; (F) caudal appendages of Telebasis filiola; (G) secondary genitalia of Tramea cophysa.
Figure 5. Morphology of male Odonata adults collected in Paraíba state, Brazil (lateral view of): (A) secondary genitalia of Anatya januaria; (B) secondary genitalia and (C) vesica spermalis of Erythrodiplax basalis; (D) caudal appendages of Hetaerina rosea; (E) secondary genitalia of Perithemis tenera; (F) caudal appendages of Telebasis filiola; (G) secondary genitalia of Tramea cophysa.
Diversity 14 00203 g005
Table 1. Similarity match approach of the generated sequences of collected specimens with the regional database sequences and global database sequences of GenBank and BOLD (Species Level Barcode Records) for species identification. N, number of barcode sequences; Mean(%) = average of intraspecific genetic distance value (expressed as percent); Max(%) = maximum intraspecific genetic distance value (expressed as percent); Nearest neighbor = most closely related species retrieved; DNN = lowest genetic distance to the nearest neighbor (expressed as percent).
Table 1. Similarity match approach of the generated sequences of collected specimens with the regional database sequences and global database sequences of GenBank and BOLD (Species Level Barcode Records) for species identification. N, number of barcode sequences; Mean(%) = average of intraspecific genetic distance value (expressed as percent); Max(%) = maximum intraspecific genetic distance value (expressed as percent); Nearest neighbor = most closely related species retrieved; DNN = lowest genetic distance to the nearest neighbor (expressed as percent).
Similarity Match Regional DatabaseSimilarity Match Global Database
SpeciesNMean (%)Max (%)Nearest NeighborDNN (%)GenBank id. (%)GenBank E-ValueGenBank SpeciesBold SimilarityBold Species
Acanthagrion gracile1--Lestes forficula19.0192.171.00 × 10−149Acanthagrion cuyabaeno matchno match
Anatya januaria30.360.57Erythrodiplax umbrata15.4499.730Anatya sp. 99.74Anatya guttata
Argia reclusa30.160.25Orthemis discolor18.8898.740Argia reclusa98.73Argia reclusa
Brachymesia furcata20.850.85Brachymesia herbida10.5599.670Brachymesia furcata99.67Brachymesia furcata
Brachymesia herbida1--Brachymesia furcata10.5590.050Brachymesia furcata99.33Brachymesia herbida
Diastatops obscura30.861.3Orthemis aequilibris12.299.830Diastatops obscura99.83Diastatops obscura
Dythemis nigra30.160.25Orthemis discolor13.871000Dythemis multipunctata100Dythemis nigra
Epipleoneura metallica1--Dythemis nigra20.5985.40Argia fumipennisno matchno match
Erythemis plebeja600Erythemis carmelita10.6490.950Erythemis peruviana100Erythemis plebeja
Erythemis vesiculosa50.320.50Erythrodiplax umbrata14.6287.500Acisoma attenboroughi99.85Erythemis vesiculosa
Erythemis peruviana20.330.33Erythemis plebeja13.0288.390Agrionoptera insignis100Erythemis peruviana
Erythemis carmelita30.100.16Erythemis plebeja10.6497.900Erythemis sp. 98.67Erythemis mithroides
Erythrodiplax avittata200Erythrodiplax basalis11.0486.930Erythrodiplax paraguayensisno matchno match
Erythrodiplax umbrata31.141.37Orthemis discolor13.7499.840Erythrodiplax umbrata99.83Erythrodiplax umbrata
Erythrodiplax leticia20.160.16Erythrodiplax cf fervida4.8792.530Erythemis sp. no matchno match
Erythrodiplax fusca150.891.85Erythrodiplax avittata15.7594.530Erythrodiplax connata99.85Erythrodiplax fusca
Erythrodiplax cf. fervida200Erythrodiplax leticia4.8791.710Erythemis sp. no matchno match
Erythrodiplax basalis60.110.75Erythrodiplax avittata11.0487.246.00 × 10−118Erythrodiplax paraguayensis99.74Erythrodiplax paraguayensis
Erythrodiplax cf. unimaculata1--Erythrodiplax leticia8.5895.866.00 × 10−163Libellulidae sp.97.11Erythrodiplax kimminsi
Hetaerina rosea30.520.78Ischnura capreolus19.698.960Hetaerina sanguinea98.95Hetaerina sanguinea
Ischnura fluviatilis40.120.26Telebasis filiola15.9699.220Ischnura fluviatilis98.68Ischnura fluviatilis
Ischnura capreolus40.250.49Lestes forficula13.21000Ischnura capreolus100Ischnura capreolus
Lestes forficula30.691.03Ischnura capreolus13.21000Lestes forficula98.78Lestes forficula
Macrothemis griseofrons60.10.17Erythrodiplax umbrata15.3988.070Agrionoptera insignisno matchno match
Mecistogaster kesselringi1--Lestes forficula14.7695.926.00 × 10−178Mecistogaster amaliano matchno match
Metaleptobasis bicornis40.120.25Telebasis corallina13.7593.426.00 × 10−163Metaleptobasis selysino matchno match
Miathyria marcella60.530.92Perithemis lais13.1987.840Agrionoptera insignis100Miathyria marcella
Micrathyria hesperis40.380.51Uracis imbuta14.3787.984.00 × 10−125Micrathyria stawiarskiino matchno match
Micrathyria ocellata30.110.16Nephepeltia belai13.761000Micrathyria ocellata100Micrathyria ocellata
Nephepeltia berlai1--Micrathyria ocellata13.7698.110Nephepeltia aequisetis98.1Nephepeltia aequisetis
Orthemis schmidt40.620.78Orthemis discolor5.0099.480Orthemis schmidti99.48Orthemis schmidti
Orthemis aequilibris1--Orthemis discolor4.9295.250Orthemis discolor99.83Orthemis aequilibris
Orthemis discolor70.510.91Orthemis aequilibris4.9299.030Orthemis sp. 100Orthemis discolor
Orthemis flavopicta1--Orthemis discolor9.6993.794.00 × 10−176Orthemis cultriformisno matchno match
Pantala flavescens30.650.65Tramea cophysa14.999.840Pantala flavescens99.83Pantala flavescens
Perithemis tenera30.871.27Brachymesia furcata13.2891.10Libellulidae sp.99.83Perithemis icteroptera
Perithemis lais1--Brachymesia furcata13.1488.210Libellulinae sp.99.25Perithemis lais
Progomphus dorsopallidus1--Perithemis tenera18.9386.224.00 × 10−171Tanypteryx hagenino matchno match
Tauriphila australis1--Orthemis flavopicta14.641000Tauriphila australis100Tauriphila australis
Telebasis corallina40.470.8Orthemis schmidt13.0888.280Telebasis willinki97.9Telebasis corallina
Telebasis filiola40.260.51Telebasis corallina13.2698.720Telebasis willinki98.96Telebasis willinki
Telebasis griffinii30.520.78Uracis imbuta17.1288.898.00 × 10−127Telebasis digiticollisno matchno match
Tramea cophysa40.831.02Brachymesia furcata12.8899.110Tramea cophysa98.34Tramea binotata
Uracis imbuta1--Nephepeltia berlai14.3594.081.00 × 10−179Uracis imbuta99.04Uracis imbuta
Zenithoptera lanei1--Telebasis corallina19.3398.248.00 × 10−166Zenithoptera sp. no matchno match
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koroiva, R.; Gomes, V.G.N.; Vilela, D.S. DNA Barcoding and New Records of Odonates (Insecta: Odonata) from Paraíba State, Brazil. Diversity 2022, 14, 203. https://doi.org/10.3390/d14030203

AMA Style

Koroiva R, Gomes VGN, Vilela DS. DNA Barcoding and New Records of Odonates (Insecta: Odonata) from Paraíba State, Brazil. Diversity. 2022; 14(3):203. https://doi.org/10.3390/d14030203

Chicago/Turabian Style

Koroiva, Ricardo, Vanessa Gabrielle Nóbrega Gomes, and Diogo Silva Vilela. 2022. "DNA Barcoding and New Records of Odonates (Insecta: Odonata) from Paraíba State, Brazil" Diversity 14, no. 3: 203. https://doi.org/10.3390/d14030203

APA Style

Koroiva, R., Gomes, V. G. N., & Vilela, D. S. (2022). DNA Barcoding and New Records of Odonates (Insecta: Odonata) from Paraíba State, Brazil. Diversity, 14(3), 203. https://doi.org/10.3390/d14030203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop