BIOLOGIJA. 2019. Vol. 65. No. 4. P. 262–272 © Lietuvos mokslų akademija, 2019
Due to a number of factors in the anthropogenic load on ecosystems, environmental changes, and competition between native and invasive species, the loss of species diversity has currently risen to an alarming scale. One of the prospective approaches to remedy the current state and halt the loss is the international project “Barcode of Life”. The report describes the results of molecular genetic research on the insect representatives of such orders as Coleoptera (beetles), Lepidoptera (butterflies), and Trichoptera (caddisflies) using the cytochrome c oxidase subunit 1 of the mitochondrial region (COI). We proposed a technique for the sampling of biological material to isolate DNA from the hind legs of individuals. The technique prevents removal of the individuals of the most valuable species collected on the territory of the Republic of Belarus from their habitats. The DNA sample collection of order representatives was investigated. For the subfamily Cetoniinae (flower chafer), significant differences in haplotypes among the representatives of the Belarusian and European parts of their areal were found. For Trichoptera, it was shown that the COI gene has high variability to differentiate species. It was also revealed that some species of caddisflies, which Belarusian researchers believe to be synonymous ones, have sequences with big differences according to the BoldSystem database. This fact should be explored in the future. So, our analysis allows considering the COI gene as a satisfactory marker for the species identification in the Belarusian insects’ taxa studied.
Keywords: DNA barcoding, butterflies, beetles, caddisflies, Belarusian populations
The loss of species diversity has currently risen to an alarming scale, due to a number of factors in anthropogenic load on ecosystems, environmental changes, and competition between native and invasive species (Adamowicz, 2015; Hebert et al., 2003). One of the prospective approaches to remedy the current state and halt the loss is the international project “Barcode of Life” initiated by Dr. Paul D. N. Hebert, University of Guelph, Canada. This project is the largest initiative in biodiversity genomics ever undertaken. Hundreds of biodiversity scientists, specialists in genomics, technologists, and ethicists representing 25 nations are working together to develop a digital identification system for all multi-cellular life based on the DNA barcoding technique (Adamowicz, 2015; Vernooy et al., 2010).
DNA barcoding is a system for fast and accurate species identification that makes an ecological system more accessible using a short DNA sequence instead of the whole genome and is used for eukaryotes (Kaur et al., 2015). A short DNA sequence generated from the standard genome region is known as a marker. This marker is different for various species:
COI gene (cytochrome c oxidase 1) is emerging as a standard barcode region for higher animals (Hebert et al., 2003). It is 658 nucleotide base pairs long in most groups, a very short sequence relative to 3 billion base pairs in the human genome. This site is located at the 5’ end and surrounded by conserved DNA sequences. Due to this, it is possible to develop conservative primers to carry out PCR (Mutanen et al., 2016). It is important that the difference in genetic distances within the species in multicellular organisms is less than 3%, while between the species this figure varies in the range. However, species boundaries are often hard to discern. The potential possibility for interbreeding creates challenges for the ascertainment of species boundaries (Mutanen et al., 2016). Despite this, the species rank has retained its status in biological science.
First studies on DNA barcoding began in the Republic of Belarus at the Belarusian State University. The staff of the Zoology Department actively cooperates with the Consortium for the Barcode of Life. At the National Academy of Sciences of Belarus, such studies, apart from the Institute of Genetics and Cytology, were undertaken at the Scientific and Practical Center for Bioresources (Balashenko et al., 2017).
The report describes the results of molecular genetic research on the insect representatives of such orders as Coleoptera (beetles), Lepidoptera (butterflies), and Trichoptera (caddisflies) using the cytochrome c oxidase subunit I of the mitochondrial region (COI). Among those taxa, 31 species of the Coleoptera and 34 species of the Lepidoptera were included into the 4th edition of the Red Book of the Republic of Belarus (2015), which significantly exceeds the number of species listed in the 3rd edition (2006).
The field material was collected on the territory of the state nature protection institution “Berezinsky Biosphere Reserve”, the state landscape reserve “Srednya Pripyat”, the state landscape reserve “Ozyory”, the state landscape reserve “Grodno Forest”, the biological reserve “Svisloch”, the national landscape reserve “Staritsa”, the national park “Belovezhskaya Pushcha”, and in other sites (Fig. 1). Imago representatives of the subfamily Cetoniinae were mostly collected manually, as well as using special traps. Manual collection was performed on flowering plants and tree trunks with flowing juice. Traps of two types were used: a trap with a lure and a pallet-like trap with a suspended load.
We propose a technique for the sampling of biological material by isolating DNA from the hind legs of individuals (Balashenko et al., 2017), which allowed us not to remove the individuals of the most valuable species collected on the territory of the Republic of Belarus from their habitat. The collected hind limb was fixed in 96% alcohol. The sample/ethanol ratio was maintained mostly at 1 : 5.
Two techniques were tested for isolating and purifying DNA samples: phenol-chloroform extraction and DNA isolation by the NucleoSpin DNA Insect kit from MACHEREY-NAGEL GmbH & Co. Most samples were obtained by phenol-chloroform extraction (Graham, 1978).
Specific oligonucleotide primers LCO1490 (5’GGT CAA CAA ATC ATA AAG ATA TTG G3’)/HCO2198 (5’TAA ACT TCA GGG TGA CCA AAA AAT CA3’) were used to amplify the target regions of the COI gene (Folmer et al., 1994; Simon, 1994). Sequencing was carried out on the genetic analyzer THERMO FISHER SCIENTIFIC ABI PRISM. For the analysis of DNA sequences, software packages SeqScanner (SeqScanner, 2018), MEGA 7, and MEGA X (Kumar et al., 2016) were used: SeqScanner mainly for inspection of the quality of sequencing and the generation of graphic reports, while MEGA 7 and MEGA X were used to align sequences, analyse nucleotide replacements, and construct a phylogenetic tree.
Some DNA samples were deposited with the State DNA Bank of a Human, Plants, Animals and Microorganisms of the Institute of Genetics and Cytology, the National Academy of Sciences of Belarus (Dromashko et al., 2018).
Lepidoptera. Nine butterfly species listed in the Red Book of the Republic of Belarus were genotyped:
This butterfly is one of the rarest butterflies of Europe, where it is mostly considered a postglacial relict species (Fig. 2). The state of its population is critical in some countries. L. helle is listed as endangered (EN) in Europe (Van Swaay et al., 2010) and is treated as a Species of Community Interest included in both Annex II and IV of the Habitats Directive of the European Commission (92/43/EEC). It should be noted that L. helle is listed as endangered (EN) in the IUCN Red List of Threatened Species (2019). Owing to strong population declines throughout its European range, targeted conservation measures are needed to avoid further local and regional extinctions. In addition, violet copper (L. helle) is considered to be an endangered species in Belarus (the Red Book of the Republic of Belarus, 2015). The violet copper is a Palearctic butterfly, widely distributed from central Europe, central and northern Fennoscandia throughout Russia to the Amur region and Mongolia (Bozano et al., 2004). Nine populations are known to have existed in Belarus. Some of those populations are thought to have become extinct due to anthropogenic transformation (Setrakova et al., 2019). For that reason, it is critical to save remaining stable populations in Belarus.
It should be noted that L. helle is particularly sensitive to environmental changes due to its (1) specialized habitat requirements, (2) low dispersal ability, and (3) comparably low fecundity (Fischer, 1997; Fischer et al., 2014). L. helle can be characterized as a hygrophilous habitat specialist, which inhabits wetlands such as wet, swampy meadows with the occurrence of larval host plant – bistort (Bistorta officinalis). Such places are very vulnerable, and taking into consideration that L. helle cannot live in biotopes with a high level of anthropogenic transformation (Habel, Meyer, 2014), we suppose that in Europe the number of favorable habitats is decreasing. At the same time, Belarus has plenty of favorable habitats for L. helle: over 14.1% of its territory is covered with bogs and the host plant (Bistorta officinalis) is widely distributed in Belarus (Setrakova et al., 2019). The populations of this species are typically small and often severely restricted in their geographical ranges (Habel, Meyer, 2014). Taking into account the low dispersal ability of L. helle, we suppose that with growing anthropogenic stress this species does not have enough time to occupy new habitats and this is one of the reasons for the critical population reduction.
Field data and genetic material of L. helle were collected in Minsk Region, Belarus. Stolbtsy District of Minsk Region is re presented by 36 samples.
Analysis of its COI gene sequence revealed low variability of the studied gene area of this species. There were only three haplotypes of the COI gene revealed in the Violet Copper in Belarus. The occurrence of one haplotype was about 95%. The other two haplotypes had two nucleotide substitutions: T/C309 and C/T651. European populations demonstrated similar trends according to the BoldSystems Data Portal (Barcode of Life Data System v 4, 2019). Sixtyseven published records with specimens from nine countries (Germany: 14, France: 14, Romania: 8, Poland: 7, Austria: 7, Russia: 6, Sweden: 5, Switzerland: 4, Finland: 2) were found; of these, 50 were similar to the most common haplotype in Belarus making 74.6%. Thus, the study revealed a limited number of mother lines in L. helle species.
Coleoptera. Thirteen species of the subfamily Cetoniinae live on the territory of Belarus (Alexandrovich et al., 1996; Tsinkevich et al., 2012). Representatives of the subfamily Cetoniinae are biologically associated with wood (especially decaying) at the stages of larvae; adults mainly live on glades and edges of forest biotopes, clearings, and other open spaces (Ryndevich, Parfianovich, 2016). The largest number of species (6) live in deciduous forests. The most favorable biotope for the habitat of bronzes is oak forests with old-growth trees, which correspond to the main environmental preferences of various types of bronzes (the presence of hollow trees for the development of larvae, a sufficient number of fodder flowering plants for adults, etc.). Original reserves of biodiversity in the anthropogenic landscape are forest belts along fields and roads (Ryndevich, Parfianovich, 2016). Due to their attachment to old-growth forests, family Cetoniinae members can be indicators of an ecosystem with low anthropogenic stress. Related to this, the interest in this subfamily in Belarus is growing. In the 4th edition of the Red Book of the Republic of Belarus (2015), six species of the subfamily Cetoniinae are included, which exceeds the number of species listed in the 3rd edition in 2004 (it was only 1 species).
During the period of field research, 11 species of beetles of the subfamily Cetoniinae from Belarus were collected. A general DNA fraction was isolated from the collected material. Data on the number of DNA samples included in the collection are presented in Table 1.
Sequencing of the COI gene region (cytochrome c oxidase I) was carried out for two species Oxythyrea funesta (Poda, 1761) (6 sequences) and Protaetia aeruginosa (Linnaeus, 1767) (3 sequences) (Fig. 3). Oxythyrea funesta is a common species, not a protected one. Protaetia aeruginosa is listed in the Red Book of the Republic of Belarus (category IV of the national conservation value). The species was included in the European list of protected saproxyl species (NT) and the Red Book of Russia (category II: species declining in numbers). Despite that, 66 and 15 samples of O. funesta and P. aeruginosa were taken for DNA extraction; qualitative sequences of the COI gene were obtained only for six and three samples, respectively.
Species | Number of samples |
---|---|
Oxythyrea funesta (Poda, 1761) | 66 |
Tropinota hirta (Poda, 1761) | 4 |
Cetonia aurata (Linnaeus, 1761) | 47 |
Protaetia lugubris (Herbst, 1786) | 3 |
P. metallica (Herbst, 1782) | 5 |
P. fieberi (Kraatz, 1880) | 2 |
P. aeruginosa (Linnaeus, 1767) | 15 |
Valgus hemipterus (Linnaeus, 1758) | 3 |
Osmoderma barnabita (Motsch., 1845) | 2 |
Gnorimus nobilis (Linnaeus, 1758) | 1 |
Trichius fasciatus (Linnaeus, 1758) | 4 |
DNA sequences were analysed using SeqScanner programs, MEGA 7. Differences in the genetic sequences of the habitat representatives from the Belarusian part were found for Oxythyrea funesta and Protaetia aeruginosa. When comparing the obtained genetic sequences with the sequences taken from the international BOLD database (22 samples collected in Russia, Finland, Austria, Germany, France, and Italy), two nucleotide replacements were found for O. funesta. They were found in the individuals of the Belarusian part of the habitat only: A/G (three sequences with genotype G/G, one sequence with genotype A/G) replacement of G/A (three sequences with genotype G/G, one sequence with genotype A/A). In the representatives of P. aeruginosa, eight variants of nucleotide substitutions were found. However, it is difficult to talk about the differences in sequences in the homozygous state that may be due to the differences in the state the populations of those species. There are only three sequences of the COI gene of this species listed under synonymic species names (Cetonischema aeruginosa, Potosia aeruginosa) in the individuals collected in different parts of the habitat (according to the BOLD database). It should also be noted that nucleotides were found to be heterozygous in the sequence of the COI O. funesta region, while in P. aeruginosa, all nucleotides of the resulting sequences appeared.
A sequence analysis of the COI gene region was also performed using 175 variable-length sequences. This number included seven sequences collected on the territory of Belarus and 168 sequences of species belonging to the Cetoniinae subfamily obtained from the BOLD Database. For the analysis, only the COI gene sequences of the specimens collected in Europe were chosen. On the basis of the analysis, a philogram was constructed. Its section is shown in Fig. 4. It can be seen that the analyzed DNA sequences of the P. aeruginosa individuals collected in Belarus and Russia are united in a separate cluster and this clearly shows the presence of significant differences in these sequences as compared to the DNA sequences of the specimens collected in France.
It should be noted that the ciphers of the samples are given in accordance with the BOLD Database, where some species are listed under synonymous names, but they are combined as one species in the construction of a philogram.
Trichoptera. The interest in the aquatic ecosystem is worldwide because it is severely altered or destroyed at a rate greater than that in any other times in the human history and way faster than they are being restored (Cairns et al., 1992). The Water Framework Directive 2000/60/EC (WFD) was adopted to help resolve this issue. One of the concepts of WFD is the use of applied ecology of water ecosystems: classification of water quality is based on direct biological measurements. It means that real biological communities should be sampled and the taxa identified. The significance of the order Trichoptera is that counting of species of this particular group is actively used in biomonitoring (EPT index). The structure of the Ephemeroptera, Plecoptera and Trichoptera (EPT) community is significantly sensitive to water parameters (Hamid et al., 2017). While morphological identification is very difficult, especially of their early stages of development, it is very important to identify the larvae of the EPT group correctly. Monitoring based on morpho-taxonomy has many drawbacks. The key step of species identification is often not taken, as it requires a high level of taxonomic expertise, which is lacking in most organisations conducting bioassesment, or species cannot be identified as they are morphologically cryptic or represent little known groups (Carew et al., 2013). Thereby, a routine method using DNA fragments for identifying representatives of the order Trichoptera is promising for future activities. Together with other methods, it could be a part of the bioassessment toolbox. To use new molecular approach for bioassessment, building a DNA barcode reference library of aquatic invertebrates is essential.
We started gap analysis for Trichoptera order registered in Belarus. Fifty-seven species were registered in Belarus (Moroz, Lipinskaya, 2014). The analysis was carried out using only the public data portal “The Barcode of Life Data Systems”. We found 4411 sequences related to the species listed as native to Belarus. For seven of the species, there are no sequences at all and 45 species have less than ten sequences. Using DNA barcoding method implies a full referenced database in the first place. Improvement of the database will allow DNA sequencing to become feasible in routine monitoring. After dealing with the priority species, we will continue with other EPT species from Belarus, which will make the use of the method of DNA barcoding more applicable for our territory.
To analyze the applicability of the DNA sequence of the COI gene for the specific identification of the caddis of Belarus, the data obtained were supplemented with the data from the BoldSystem Database, and, using the MEGA Х program, cladograms of the evolutionary history of caddis species were constructed. Only DNA sequences of the species recorded on the territory of Belarus in accordance with the Belarusian Catalog of Mayflies (Ephemeroptera), Vesnjans (Plecoptera), and Caddisflies (Trichoptera) were chosen for the analysis. Table 2 shows the genetic distances between the species of the family Hydropsychidae found on the territory of Belarus.
Cheumatopsyche lepida | |||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Cyrnus crenaticornis | 0.33 | ||||||||||||||||||
Cyrnus flavidus | 0.33 | 0.15 | |||||||||||||||||
Cyrnus insolutus | 0.33 | 0.12 | 0.16 | ||||||||||||||||
Cyrnus trimaculatus | 0.33 | 0.09 | 0.16 | 0.12 | |||||||||||||||
Holocentropus dubius | 0.32 | 0.18 | 0.21 | 0.18 | 0.2 | ||||||||||||||
Holocentropus picicornis | 0.33 | 0.14 | 0.17 | 0.19 | 0.16 | 0.14 | |||||||||||||
Holocentropus stagnalis | 0.35 | 0.17 | 0.17 | 0.2 | 0.19 | 0.19 | 0.18 | ||||||||||||
Hydropsyche angustipennis | 0.21 | 0.28 | 0.31 | 0.31 | 0.27 | 0.32 | 0.28 | 0.32 | |||||||||||
Hydropsyche bulgaromanorus | 0.21 | 0.28 | 0.29 | 0.32 | 0.28 | 0.31 | 0.26 | 0.33 | 0.13 | ||||||||||
Hydropsyche contubernalis | 0.23 | 0.27 | 0.33 | 0.32 | 0.31 | 0.31 | 0.27 | 0.33 | 0.16 | 0.12 | |||||||||
Hydropsyche guttata | 0.2 | 0.27 | 0.33 | 0.32 | 0.31 | 0.32 | 0.28 | 0.31 | 0.14 | 0.08 | 0.14 | ||||||||
Hydropsyche pellucidula | 0.25 | 0.33 | 0.33 | 0.33 | 0.35 | 0.33 | 0.31 | 0.34 | 0.18 | 0.18 | 0.18 | 0.19 | |||||||
Hydropsyche saxonica | 0.24 | 0.32 | 0.32 | 0.33 | 0.33 | 0.31 | 0.31 | 0.35 | 0.18 | 0.18 | 0.18 | 0.19 | 0.02 | ||||||
Hydropsyche siltalai | 0.25 | 0.26 | 0.35 | 0.33 | 0.32 | 0.32 | 0.33 | 0.35 | 0.16 | 0.17 | 0.17 | 0.15 | 0.02 | 0.19 | |||||
Neureclipsis bimaculata | 0.33 | 0.19 | 0.17 | 0.18 | 0.17 | 0.23 | 0.2 | 0.2 | 0.29 | 0.29 | 0.3 | 0.31 | 0.35 | 0.37 | 0.32 | ||||
Plectrocnemia conspersa | 0.32 | 0.14 | 0.17 | 0.16 | 0.17 | 0.18 | 0.17 | 0.17 | 0.27 | 0.27 | 0.29 | 0.28 | 0.32 | 0.34 | 0.3 | 0.19 | |||
Polycentropus flavomaculatus | 0.33 | 0.14 | 0.14 | 0.18 | 0.14 | 0.18 | 0.17 | 0.19 | 0.27 | 0.27 | 0.3 | 0.3 | 0.33 | 0.35 | 0.3 | 0.21 | 0.14 | ||
Polycentropus irroratus | 0.33 | 0.17 | 0.19 | 0.2 | 0.2 | 0.24 | 0.19 | 0.23 | 0.3 | 0.28 | 0.3 | 0.29 | 0.32 | 0.34 | 0.32 | 0.23 | 0.18 | 0.12 | |
Polycentropus kingi | 0.32 | 0.14 | 0.14 | 0.18 | 0.15 | 0.2 | 0.18 | 0.2 | 0.27 | 0.26 | 0.3 | 0.29 | 0.32 | 0.34 | 0.31 | 0.23 | 0.15 | 0.05 | 0.1 |
As an example, Fig. 5 shows a circular cladogram of the evolutionary history of the Hydropsychidae family of the order Trichoptera. The cladogram shows that the COI gene makes it possible to determine the points of the species divergence within the family. It was also found that the COI gene has a sufficiently high variability within the Trichoptera order, which makes it possible to clearly separate the types of this order. A total of 329 variable nucleotide positions from 659 nucleotides of the COI gene were detected, which is 49.9%.
In the course of the analysis, it was found that some species considered synonymous in Belarus are listed in the BoldSystem Database (BoldSystems, 2018) with the corresponding sequences, which differ quite significantly among themselves. These include Hydropsyche angustipennis (Curtis, 1834)/Hydropsyche fulvipes (Stephens, 1836); Neureclipsis bimaculata (Linnaeus, 1761)/Philopotamus variegatus (Scopoli, 1763). At the same time, synonymous names of species are often found in the BoldSystem; however, the analysis of genetic sequences allows combining them into one species. As an example, refer to Fig. 4, where a cladogram section is presented based on the analysis of the partial nucleotide sequence of the COI gene of the Cetoniinae subfamily. It shows that synonymous species of Protaetia marmorata/Protaetia lugubris are combined into one species.
1. We proposed a technique for biological material sampling to isolate DNA from the hind legs of individuals, which prevented the removal of the individuals of the most valuable species collected on the territory of the Republic of Belarus from their habitats.
2. A collection of DNA samples of the studied order representatives (i.e., Lepidoptera, Coleoptera and Trichoptera) was compiled. Some samples were deposited with the Republican DNA Bank of a Human, Plants, Animals and Microorganisms of the Institute of Genetics and Cytology, the National Academy of Sciences of Belarus.
3. For some rare species of the order Lepidoptera, low haplotype diversity was shown. As an example, only three haplotypes in the Violet Copper (Lycaena helle) were revealed and one of them has about 95% of occurrence. In the subfamily Cetoniinae (flower chafer) of the order Coleoptera, significant differences were found in haplotypes among representatives of the Belarusian and European parts of their areal.
4. For Trichoptera, it was shown that the COI gene has high variability to differentiate species. It was also revealed that some caddisfly species, which Belarusian researchers believe to be synonymous ones, have sequences with serious differences according to the BoldSystem Database. This fact should be explored in the future.
5. Our analysis allows us to consider the COI gene as a pretty good marker for the species identification in the Belarusian insects’ taxa studied. So, the white spot Belarus occupies in the BoldSystem has become smaller.
Received 28 June 2019
Accepted 25 September 2019
* Corresponding author. Email: s.dromashko@igc.by
Santrauka
Rūšių įvairovės nykimas, susijęs su antropogenine apkrova ekosistemoms, aplinkos pokyčiais, vietinių bei invazinių rūšių konkurencija ir kitais veiksniais, šiuo metu tapo nerimą keliančiu ženklu. Vienas iš perspektyvių būdų ištaisyti dabartinę būklę ir sustabdyti rūšių įvairovės nykimą yra tarptautinis projektas „Gyvenimo barkodingas“. Jo ataskaitoje aprašomi vabzdžių, tokių kaip Coleoptera (vabalai), Lepidoptera (drugiai) ir Trichoptera (apsiuvos), molekulinių genetinių tyrimų, naudojant mitochondrinio regiono (COI) citochromo c oksidazės 1 subvienetą, rezultatai. Mes pasiūlėme biologinės medžiagos mėginių ėmimo metodą, padedantį išskirti DNR iš užpakalinių galūnių. Buvo ištirta skirtingų vabzdžių klasės būrių atstovų DNR mėginių kolekcija. Cetoniinae, skarabėjų pošeimyje, rasta reikšmingų haplotipų skirtumų tarp Baltarusijos ir Europos arealo dalių individų. Trichoptera atveju nustatytas didelis COI geno kintamumas atskirose rūšyse. Atskleista, kad kai kurios apsiuvų rūšių sekos, Baltarusijos tyrinėtojų laikytos identiškomis, pagal „BoldSystem“ duomenų bazę, turi rimtų skirtumų. Šį faktą ateityje reikėtų ištirti. Taigi, mūsų analizė leidžia apsvarstyti COI geną kaip tinkamą žymenį rūšims identifikuoti Baltarusijos vabzdžių taksonuose.
Raktažodžiai: DNR barkodingas, drugiai, vabalai, apsiuvos, Baltarusijos vabzdžių populiacija