W obliczu bezprecedensowego tempa wymierania gatunków na świecie w różnych taksonach1,2,3,4, dążenie do stworzenia kompleksowego inwentarza różnorodności biologicznej, obejmującego gatunki dotychczas nieodkryte i nieopisane, staje się wyścigiem z czasem2,5. Po pierwsze, ochrona różnorodności biologicznej ma sens jedynie w świetle właściwej wiedzy taksonomicznej6,7. Zatem identyfikacja różnorodności biologicznej na poziomie gatunku jest niezbędna, biorąc pod uwagę, że blisko spokrewnione, lecz odmienne gatunki mogą zajmować unikalne nisze fundamentalne8.
W systematyce chrząszczy taksonomia alfa, czyli taksonomia tradycyjna, stanowi złoty standard w identyfikacji gatunków, podobnie jak w większości, jeśli nie we wszystkich grupach zwierząt9. Podejście to opiera się na cechach morfologicznych w celu klasyfikacji organizmów na wyższych poziomach taksonomicznych oraz zapewnienia identyfikacji na poziomie gatunku10,11. W wielu przypadkach identyfikacja na poziomie rodziny12 lub nawet rodzaju13 może być przeprowadzona wyłącznie na podstawie zewnętrznych cech morfologicznych, takich jak kształt ciała i specjalistyczne adaptacje. Tymczasem identyfikacja gatunkowa chrząszczy jest zazwyczaj przeprowadzana poprzez porównanie szczegółów strukturalnych aedeagusa lub narządów kopulacyjnych samca14,15,16.
Pomimo szerokiej akceptacji tej opartej na morfologii alfa-taksonomii, podejście to do odkrywania gatunków jest utrudnione przez liczne „bariery taksonomiczne”6,17,18, zwłaszcza wśród bezkręgowców19. W systematyce chrząszczy identyfikacja i/lub odkrywanie gatunków za pomocą alfa-taksonomii napotyka trudności związane z ograniczoną liczbą wykwalifikowanych taksonomików20,21, zwięzły i mało informacyjne wcześniejsze opisy22 oraz problemy logistyczne związane z dostępem do literatury i okazów typowych23.
Sytuacja ta jest jeszcze bardziej utrudniona, gdy rozważany takson jest kryptyczny lub mało charakterystyczny. Gatunki kryptyczne to zbiór gatunków z tego samego rodzaju, a nawet z różnych rodzajów, których nie można łatwo rozróżnić na podstawie morfologii porównawczej ze względu na subtelne międzygatunkowe różnice fenetyczne24. Niektóre taksony kryptyczne, określane jako „dark taxa”25,26,27, dodatkowo charakteryzują się wysoką różnorodnością przy jednoczesnym dużym niedoborze badań. Niestety, zjawisko to nie jest rzadkie w przypadku chrząszczy nadrzecznych i wodnych8,28. Brak precyzyjnego wyznaczenia granic gatunkowych wśród kryptycznych „dark taxa” zagraża różnorodności biologicznej, ponieważ jednolite środki ochrony stosuje się wobec organizmów o potencjalnie różnych niszach i wymaganiach ekologicznych8,29.
W obliczu tych wyzwań, „taksonomia integracyjna”23, czyli podejście łączące taksonomię alfa z innym rodzajem dowodów, zyskało popularność w ciągu ostatnich dwóch dekad jako metoda opisywania nowych gatunków30,31,32. Jednym z takich rodzajów dowodów wykorzystywanych w badaniach taksonomii integracyjnej chrząszczy33,34,35,36,37,38,39 oraz innych rzędów owadów40,41,42,43,44,45 jest sekwencjonowanie DNA. W szczególności gen mitochondrialnej podjednostki I cytochromu c oksydazy (COI) jest wykorzystywany do kodowania kreskowego (barcodingu) różnorodnych grup zwierząt ze względu na jego umiarkowany konserwatyzm46,47,48. Podczas gdy jego 658bp-długi koniec 5′ (COI-5′), znany również jako „fragment kodowania kreskowego” lub „fragment Folmera”, został zaproponowany jako sekwencja „diagnostyczna”49, stosowany jest również 723 bp-długi fragment COI-3′ ze względu na dostępność solidnego zestawu starterów dla końca 3′50,51,52.
Podczas gdy kod kreskowy DNA z wykorzystaniem COI został opracowany do identyfikacji na poziomie gatunków, jest on ograniczony przez repozytorium sekwencji a priori 53, takie jak GenBank54 i BOLD55; podejścia do molekularnego wyznaczania granic gatunków określają granice gatunkowe, wykorzystując jedynie dostarczone sekwencje DNA, bez konieczności posiadania wcześniejszego repozytorium56. Istnieją dwie klasy podejść do molekularnego wyznaczania granic gatunków, w zależności od pliku wejściowego. Po pierwsze, stosując dopasowanie sekwencji jako dane wejściowe, podejścia oparte na progach wykorzystują wartości odcięcia do określenia zróżnicowania między gatunkami i wewnątrz gatunków48,57,58. Jako metody oparte na dystansie, podejścia te opierają się na zaobserwowanym maksymalnym dystansie wewnątrzgatunkowym wynoszącym 3% dla owadów, co określa się mianem „luki kodowania” (barcoding gap)51,59. Do podejść opartych na progach należą m.in. SpeciesIdentifier zaimplementowany w TaxonDNA60, dystans genetyczny via model dwuparametrowy Kimury (K2P)61 oraz Assemble Species by Automatic Partitioning (ASAP)62.
Po drugie, wykorzystując drzewo jako plik wejściowy, podejścia oparte na teorii koalescencji identyfikują granice gatunków i sortują niezależne linie poprzez określanie zmian w tempie rozgałęziania linii63,64,66. Często opierając się na filogenetycznej koncepcji gatunku67, podejścia oparte na koalescencji obejmują procesy rozkładu Poissona dla drzew (PTP)68 oraz ich wariant, multi-rate PTP69. W każdym przypadku klastry utworzone przez metody molekularnego delimitowania gatunków określa się zbiorczo jako molekularne operacyjne jednostki taksonomiczne (MOTUs).
W niniejszej pracy przedstawiamy metody potoku analizy służącego do generowania MOTU, obejmujące trzy podejścia oparte na progach oraz trzy podejścia oparte na koalescencji (Rysunek 1). Potok został przetestowany na dwóch zbiorach danych, mianowicie (1) sekwencjach COI-3′ chrząszczy z rodzaju Byrrhinus Motschulsky, 1858 żyjących w strefach nadrzecznych, składającym się z dwóch opublikowanych gatunków35 wyodrębnionych wcześniej za pomocą taksonomii integratywnej oraz dodatkowych sekwencji z nieopisanych gatunków, oraz (2) sekwencjach COI-5′ chrząszczy z rodzaju Anacaena Thomson, 1859, składającym się z sekwencji dwóch opublikowanych gatunków70 wyodrębnionych wcześniej za pomocą alfa-taksonomii oraz dodatkowych sekwencji z nieopisanych gatunków.

Rysunek 1: Schemat blokowy procesu generowania MOTU. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.