È richiesta un abbonamento a JoVE per visualizzare questo contenuto. Accedi o inizia la tua prova gratuita.

Articolo metodologico

Una pipeline di intruglio per la generazione di unità tassonomiche operative molecolari (MOTU) tra coleotteri ripariali e acquatici

914 visualizzazioni

DOI:

10.3791/68323

11 luglio 2025

In questo articolo

Sommario

Il protocollo qui presentato mira a guidare gli utenti nell'elaborazione computerizzata di sequenze geniche della subunità I della citocromo c ossidasi (COI) generate dai coleotteri, in modo tale che l'ipotesi di raggruppamento delle specie chiamata unità tassonomiche operative molecolari (MOTU) possa essere generata dal DNA.

Abstract

Il declino della biodiversità sta avvenendo a ritmi inimmaginabili, ma il monitoraggio della biodiversità è ostacolato da una pletora di fattori, tra cui l'inventario incompleto della diversità di base e la diminuzione della popolazione di tassonomisti qualificati. Colmare questo ostacolo nell'inventario delle specie di "taxa oscuri", o gruppi iperdiversi spesso poco studiati nella tassonomia, è particolarmente cruciale nei tropici che affrontano pressioni ambientali insormontabili. Mentre l'alfa-tassonomia basata sulla morfologia rimane il gold standard nell'identificazione delle specie, i metodi basati sul DNA hanno mostrato un enorme potenziale nell'accelerare la delineazione e la scoperta delle specie. Qui, descriviamo una pipeline di intruglio per la generazione di unità tassonomiche operative molecolari (MOTU) che coinvolge sei approcci di delimitazione delle specie molecolari implementati su due generi di coleotteri altamente poco appariscenti, vale a dire Byrrhinus Motschulsky, 1858 della famiglia Limnichidae e Anacaena Thomson, 1859 della famiglia Hydrophilidae, dalle Filippine. Poiché questo protocollo si concentra sull'elaborazione dei codici a barre del DNA e dato che questi dati genetici possono essere ottenuti sia mediante sequenziamento che scaricando da database online, il punto di partenza per il protocollo qui descritto sono le sequenze di DNA. Dati i codici a barre COI, la pipeline utilizza una combinazione di tre approcci di delimitazione delle specie molecolari basati su soglie, vale a dire TaxonDNA, K2P e ASAP, che utilizzano l'allineamento delle sequenze come dati di input. Inoltre, la pipeline procede con tre approcci basati sulla coalescenza, ovvero PTP-ML, PTP-BI e mPTP, che richiedono un file ad albero come dati di input. Alcuni approcci prevedono l'uso di server web, mentre altri utilizzano software locale. Con la selezione dei metodi algoritmici presentati, i MOTU possono essere identificati, sia per consenso che a maggioranza. Sorprendentemente, questa conduttura ha delineato specie di coleotteri quasi cospicui provenienti dalle Filippine, anche quelle documentate dalla stessa località o da località adiacenti.

Introduzione

Con i tassi senza precedenti di declino delle specie globali in vari taxa 1,2,3,4, l'impulso a creare un inventario completo della biodiversità, che includa specie precedentemente non scoperte e non descritte, diventa una corsa contro il tempo 2,5. Per prima cosa, la conservazione della biodiversità ha senso solo alla luce di un'adeguata conoscenza tassonomica 6,7. Pertanto, l'identificazione della biodiversità a livello di specie è necessaria, dato che specie strettamente correlate ma diverse possono avere nicchie fondamentali uniche8.

Nella sistematica dei coleotteri, l'alfa-tassonomia o tassonomia tradizionale è il gold standard nell'identificazione delle specie simile alla maggior parte, se non a tutti, i gruppi di animali9. Questo approccio si basa sui caratteri morfologici per classificare gli organismi a una tassonomia più alta e per fornire l'identificazione a livello di specie10,11. In molti casi, l'identificazione della famiglia12 o anche a livello di genere13 può essere effettuata utilizzando solo caratteristiche morfologiche esterne, come la forma del corpo e adattamenti specializzati. Nel frattempo, l'identificazione a livello di specie nei coleotteri viene solitamente effettuata attraverso il confronto dei dettagli strutturali dell'aedeagus o dei genitali maschili 14,15,16.

Nonostante l'ampia accettazione di questa alfa-tassonomia basata sulla morfologia, questo approccio alla scoperta delle specie è ostacolato da numerosi "impedimenti tassonomici"6,17,18, specialmente tra gli invertebrati19. Nella sistematica dei coleotteri, l'identificazione e/o la scoperta delle specie utilizzando l'alfa-tassonomia affronta le sfide di avere un numero limitato di tassonomisti qualificati20,21, di descrizioni precedenti brevi e scarsamente informative22 e di problemi logistici legati all'accesso alla letteratura e agli esemplari tipo23.

Questo è anche esacerbato quando il taxon in questione è criptico o altamente poco appariscente. La specie criptica si riferisce all'insieme di specie dello stesso genere, o anche di generi diversi, che non possono essere facilmente delineate dalla morfologia comparativa, date le sottili differenze fenetiche interspecifiche24. Alcuni taxa criptici, considerati come "taxa oscuri"25,26,27, soffrono inoltre di essere iperdiversi ma fortemente poco studiati. Sfortunatamente, questo non è raro nei coleotteri ripariali e acquatici 8,28. L'incapacità di affrontare la delineazione delle specie tra criptici taxa oscuri minaccia la biodiversità poiché viene data una misura di conservazione unica per tutti a organismi di nicchie e requisiti ecologici potenzialmente diversi 8,29.

Alla luce di queste sfide, la "tassonomia integrativa"23, o l'approccio di accoppiare l'alfa-tassonomia con un'altra linea di evidenza, ha guadagnato terreno come mezzo per erigere nuove specie negli ultimi due decenni 30,31,32. Una di queste linee di evidenza utilizzate negli studi tassonomici integrativi sui coleotteri 33,34,35,36,37,38,39 e altri ordini di insetti 40,41,42,43,44,45 è la sequenza del DNA. In particolare, il gene della subunità I (COI) della citocromo c ossidasi mitocondriale viene utilizzato per codificare a barre un insieme diversificato di animali a causa della sua natura moderatamente conservativa 46,47,48. Mentre la sua estremità 5' lunga 658 bp (COI-5'), nota anche come 'frammento di codifica a barre' o 'frammento di Folmer', è stata proposta come sequenza 'diagnostica'49, la COI-3' lunga 723 bp viene utilizzata anche grazie a un robusto set di primer per l'estremità 3' 50,51,52.

Mentre il DNA barcoding che utilizza COI è stato concepito per l'identificazione a livello di specie, è limitato da un repository di sequenze a priori 53, come GenBank54 e BOLD55, gli approcci molecolari di delimitazione delle specie stabiliscono i limiti delle specie utilizzando solo le sequenze di DNA fornite senza alcun repository di prerequisiti56. Esistono due classi di approcci di delimitazione delle specie molecolari basati sul file di input. In primo luogo, utilizzando l'allineamento della sequenza come input, gli approcci basati sulla soglia utilizzano valori di cut-off per determinare la divergenza tra e all'interno delle specie 48,57,58. Come metodi basati sulla distanza, questi approcci si basano su una distanza intraspecifica massima osservata del 3% per gli insetti, che viene definita "barcoding gap"51,59. Alcuni degli approcci basati su soglie sono SpeciesIdentifier implementato in TaxonDNA60, la distanza genetica tramite Kimura 2-parameter (K2P)61 e Assemble Species by Automatic Partitioning (ASAP)62.

In secondo luogo, utilizzando l'albero come file di input, gli approcci basati sulla coalescenza identificano i confini delle specie e ordinano i lignaggi indipendenti determinando il cambiamento nel tasso di ramificazione del lignaggio 63,64,66. Spesso basati sul concetto filogenetico delle specie67, gli approcci basati sulla coalescenza includono i processi dell'albero di Poisson (PTP)68 e la loro variazione, PTP69 multi-rate. In ogni caso, i cluster formati dagli approcci di delimitazione delle specie molecolari sono collettivamente indicati come unità tassonomiche operative molecolari (MOTU).

Qui, presentiamo i metodi per una pipeline di intruglio per la generazione di MOTU che coinvolge tre approcci basati su soglie e tre approcci basati sulla coalescenza (Figura 1). La pipeline è testata su due set di dati, vale a dire (1) sequenze COI-3' del ripariale Byrrhinus Motschulsky, 1858 beetles, che è composto da due specie pubblicate35 precedentemente erette utilizzando la tassonomia integrativa e sequenze aggiuntive da specie non descritte, e (2) sequenze COI-5' di Anacaena Thomson, 1859 beetles, che è composto da sequenze di due specie pubblicate70 precedentemente eretto utilizzando l'alfa-tassonomia e sequenze aggiuntive di specie non descritte.

figure-introduction-1
Figura 1: Diagramma di flusso della pipeline di intruglio per la generazione di MOTU. Clicca qui per visualizzare una versione più grande di questa figura.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Protocollo

I materiali di partenza per questa pipeline di intruglio sono sequenze di DNA, che possono essere scaricate da repository online (ad esempio, BOLD, GenBank) o generate dal sequenziamento Sanger71. La Tabella 1 e la Tabella 2 elencano i numeri di accesso per i set di dati i cui risultati rappresentativi sono presentati in questo articolo.

1. Fase preparatoria per gli approcci basati su soglie: allineamento delle sequenze

  1. Sequenze di carico in MEGA X72.
    NOTA: MEGA X può essere scaricato da https://www.megasoftware.net/
    1. Fare clic su Allinea per aprire un menu a discesa. Fare clic su Modifica/Crea allineamento e selezionare Crea un nuovo allineamento. Fare clic su OK per confermare la selezione. Selezionare DNA come tipo di dati per l'allineamento.
    2. Passate il mouse sulla scheda Modifica e selezionate Inserisci sequenza da file. Naviga nella directory con le sequenze e selezionale per essere caricate in MEGA.
  2. Allinea le sequenze facendo clic su Allineamento, quindi su Allinea con ClustalW73. Continua con le impostazioni predefinite e fai clic su OK.
  3. Modifica manualmente le sequenze tagliando entrambe le estremità e pulendo eventuali inserimenti o eliminazioni.
    1. Per eliminare gli inserimenti, fare clic sulle basi e/o sulle posizioni inserite e premere Elimina sulla tastiera.
    2. Per correggere le cancellazioni, fare clic sulla posizione che si suppone sia stata eliminata come indicato da "-". Elimina il "-" e digita la base desiderata.
    3. Trova la prima posizione in cui tutte le sequenze hanno un carattere. Fare clic sulla casella vuota nell'intestazione della riga della posizione a sinistra della posizione determinata e trascinare fino a selezionare le posizioni in eccesso per tutte le sequenze fino alla fine iniziale. Premere Elimina per tagliare.
    4. Ripetere il passaggio 1.3.3 cercando l'ultima posizione in cui tutte le sequenze hanno un carattere. Fare clic sulla casella vuota sulla riga di intestazione della posizione a destra della posizione determinata, trascinare fino a selezionare le posizioni in eccesso per tutte le sequenze fino alla fine del terminale. Premere Elimina.
      NOTA: La Figura 2A mostra un buon allineamento delle sequenze di DNA. Il set di dati finale è una matrice di sequenze di uguale lunghezza; In tal modo, la lunghezza della sequenza non è un fattore che influenza i risultati dello studio. Le sequenze di riferimento servono come base per la lunghezza finale della matrice di allineamento.
  4. Tradurre la sequenza per verificare la presenza di codoni di stop, dato che il COI è una sequenza codificante proteine.
    1. Seleziona tutte le sequenze, quindi fai clic sulla scheda Sequenze proteiche tradotte .
    2. Quando richiesto, verificare il codice genetico come Mitocondriale invertebrato. Se il codice genetico è diverso, fare clic su No e apparirà un menu che consente di spuntare la casella per il codice genetico mitocondriale degli invertebrati .
    3. Se i codoni di stop, come indicato da asterischi nel corpo di allineamento, sono presenti per un'intera colonna, fare clic su Sequenze DNA ed eliminare la prima posizione per tutte le sequenze.
    4. Ripetere il passaggio precedente se i codoni di stop sono ancora presenti in qualsiasi sequenza.
    5. Se i codoni di stop sono ancora presenti, ripristinare le posizioni precedentemente eliminate. Seleziona tutte le sequenze, quindi passa il mouse sulla scheda Dati e fai clic su Inverti.
    6. Ripetere il passaggio 1.4.1. Se i codoni di stop sono ancora presenti, ripetere i passaggi 1.4.3 e 1.4.4.
      NOTA: La Figura 2B mostra un buon allineamento dei codoni.
  5. Fare clic su Sequenze DNA e salvare l'allineamento della sequenza in .mas/x. Esporta la sequenza in altri tipi di file pertinenti, come il file .meg e il file .fasta.
    NOTA: Alcuni approcci di delimitazione delle specie molecolari funzionano meglio senza l'outgroup. Pertanto, una possibile considerazione consiste nel generare una copia del file di allineamento senza l'outgroup.

figure-protocol-1
Figura 2: Allineamento della sequenza in MEGA X. (A) Allineamento finale della sequenza di DNA pulito manualmente e (B) allineamento della sequenza proteica tradotta del set di dati di Byrrhinus che non mostra inserzioni, delezioni e lacune. Clicca qui per visualizzare una versione più grande di questa figura.

2. Delimitazione per modulo DNA taxon nell'identificatore di specie 1.8

NOTA: L'identificatore di specie può essere scaricato da https://taxondna.sourceforge.net/

  1. Apri TaxonDNA. Fare clic su Importa, quindi fare clic su FASTA. Carica l'allineamento in formato .fasta.
  2. Fare clic su Moduli e fare clic su Cluster. Imposta la soglia al 3%. Quindi, selezionare Genera informazioni sul singolo cluster (Figura 3).
  3. Clicca su Crea cluster ora!. Salva i risultati facendo uno screenshot.

figure-protocol-2
Figura 3: TaxonDNA nella finestra dell'identificatore di specie 1.8 che mostra la soglia del 3% identificata dall'utente per la delimitazione. Clicca qui per visualizzare una versione più grande di questa figura.

3. Delimitazione con Kimura a 2 parametri (K2P) in MEGA X

  1. In MEGA, apri il file .meg facendo clic su File, quindi Apri un file/sessione. Fare clic su Distanza, quindi fare clic su Calcola distanze a coppie. Confermare il file .meg per la delimitazione.
  2. Passa il mouse sopra la casella gialla accanto a Modello/Metodo, quindi fai clic sulla freccia che appare a destra della casella.
  3. Seleziona il modello Kimura a 2 parametri dal menu a discesa. Fare clic su Ok per eseguire il programma. Aprire la finestra di output dei dati (Figura 4).
  4. Fare clic su File, quindi su Esporta/Stampa distanze.
  5. Apri il menu a discesa accanto a Formato di output e seleziona XL: cartella di lavoro di Microsoft Excel. Per Cifre decimali, scegli 4.
  6. Salva il foglio di calcolo risultante in Microsoft Excel. Usa la formattazione condizionale per evidenziare quelli con una distanza inferiore a 0,03 (o 3%).
    NOTA: A differenza degli altri approcci illustrati in questo articolo, K2P non fornisce immediatamente il clustering perché dipende dalla soglia identificata dall'utente. Questa pipeline sta lavorando sull'ipotesi che la soglia sia del 3%, il che significa che qualsiasi coppia di sequenze la cui distanza genetica è inferiore al 3% appartiene allo stesso cluster molecolare.

figure-protocol-3
Figura 4: Finestra risultante in MEGA X che mostra i risultati del parametro Kimura-2. Clicca qui per visualizzare una versione più grande di questa figura.

4. Delimitazione per specie assemblate mediante partizionamento automatico (ASAP)

  1. Accedere al server Web ASAP: https://bioinfo.mnhn.fr/abi/public/asap/asapweb.html
  2. Carica l'allineamento facendo clic sulla casella arancione con l'etichetta Scegli un file... e selezionando il file .fasta. In alternativa, trascina il file .fasta nella casella arancione. Scorri verso il basso e fai clic su Vai.
  3. Per scaricare il clustering, fare clic sull'elenco per la riga con il punteggio ASAP più basso e il rango p-val più alto (Figura 5).
    NOTA: L'output del risultato è un file .txt, che presenta sottoinsiemi dei cluster molecolari.

figure-protocol-4
Figura 5: Finestra risultante in ASAP che mostra le possibili delimitazioni. Il primo risultato (casella blu) è accettato a causa del punteggio ASAP più basso (casella rossa) e del grado P-val più alto (casella gialla). Clicca qui per visualizzare una versione più grande di questa figura.

5. Fase preparatoria per approcci basati sulla coalescenza: stima dell'albero

  1. Apri il file .meg facendo clic su Dati, quindi Apri un file/sessione. Testare il modello di sostituzione più adatto.
    1. Fai clic su Modelli, quindi trova i migliori modelli di DNA/proteine (ML). Confermare il file .meg per l'analisi.
    2. Fare clic su Ok nel menu Preferenze di analisi (Figura 6A). Prendere nota del modello migliore che utilizza il valore BIC e AICc più basso (Figura 6B).
      NOTA: Di solito è il primo nell'elenco dei modelli nella finestra risultante. Fare riferimento ai nomi completi dei nomi dei modelli abbreviati nella parte inferiore della finestra risultante.
  2. Generare l'albero.
    1. Fare clic su Filogenesi e quindi su Costruisci/Testa l'albero della massima verosimiglianza.
    2. Utilizzare il modello di sostituzione nucleotidica più adatto. Fallo cliccando sulla casella gialla accanto a Modello/Metodo. Seleziona il modello dal menu a discesa.
    3. Se il modello di sostituzione nucleotidica più adatto ha '+G' o '+I' o entrambi, fare clic sulla casella accanto a Tassi tra siti e selezionare i parametri (G, I, G + I) dal menu a discesa (Figura 6C).
    4. In Test filogenesi, fai clic sulla casella accanto a Test di filogenesi e seleziona Metodo Bootstrap dal menu a discesa. Fare clic sulla casella accanto a "No. di repliche Bootstrap" e digitare 1000 (Figura 6C). Fare clic su Ok per eseguire l'analisi.
  3. Una volta aperta una finestra contenente l'albero risultante (Figura 6D), salvare la sessione dell'albero come .mts/x. Salvare l'output come file con estensione nwk. Salva anche l'albero come file .png immagine.

figure-protocol-5
Figura 6: Stima dell'albero in MEGA X come fase preparatoria per approcci basati sulla coalescenza. (A) Impostazione per determinare il modello di sostituzione nucleotidica più adatto. (B) Per il set di dati di Byrrhinus , il modello più adatto è GTR + I a causa dei valori BIC e AICc più bassi (riquadro rosso). (C) Nell'esecuzione dell'albero della massima verosimiglianza, GTR (casella gialla) + I (casella blu) sono stati impostati come parametri. (D) L'albero di ML grezzo risultante da utilizzare nelle analisi successive. Clicca qui per visualizzare una versione più grande di questa figura.

6. Delimitazione mediante processi ad albero di Poisson (PTP)

NOTA: I passaggi seguenti produrranno i risultati di PTP-ML e PTP-BI.

  1. Vai al server web PTP: https://species.h-its.org/ptp/
  2. Carica un file newick facendo clic su Scegli file. In Il mio albero è, selezionare Radicato (Figura 7A).
  3. Inserire i parametri per l'analisi. Sotto il n. Generazioni MCMC, inserire 100000. In Diluizione (Thinning), inserite 100. In Burn-in, immettere 0.1. In Seme, inserisci 123.
  4. Inserire l'outgroup digitando il nome delle punte nella casella sotto Nomi taxa outgroup (se presenti). Le istruzioni per la formattazione dell'input si trovano a destra della casella di input.
  5. Inserisci l'indirizzo e-mail nel campo corrispondente e fai clic su Invia. Nella finestra dei risultati della delimitazione delle specie PTP, fare clic qui per aprire la finestra dei risultati.
  6. Nella soluzione di massima verosimiglianza, scaricare i risultati di PTP-ML facendo clic su Scarica risultati di delimitazione (Figura 7B).
  7. Nella soluzione supportata da Bayesian, scaricare i risultati di PTP-BI facendo clic su Scarica risultati della delimitazione.
    NOTA: i risultati indicherebbero gruppi di sequenze nelle partizioni. Questi saranno etichettati come specie numerate seguite da un valore di supporto, in cui 1.000 è l'ideale.

figure-protocol-6
Figura 7: Delimitazione delle specie molecolari mediante processi ad albero di Poisson. (A) Finestre che mostrano i parametri di input. (B) Finestra risultante nel server web PTP, con collegamento contenente i risultati della delimitazione (riquadro rosso). Clicca qui per visualizzare una versione più grande di questa figura.

7. Delimitazione mediante processi ad albero di Poisson multi-rate (mPTP)

  1. Vai al server web mPTP: https://mptp.h-its.org/#/tree
  2. Carica il file newick trascinandolo sul rettangolo grigio o facendo clic sul rettangolo. Una volta caricati i dati sulla stessa pagina, fare clic su Procedi alla selezione dell'outgroup (Figura 8A).
  3. Nella pagina delle specifiche dell'outgroup , selezionare l'outgroup facendo clic sulla casella accanto alle etichette dei taxa degli esemplari dell'outgroup. Fare clic sull'opzione Ritaglia outgroup nella parte inferiore della sequenza. Fare clic su Selezione modello.
  4. Fare clic su MPTP, quindi fare clic su Opzioni di visualizzazione. Nella pagina Opzioni di visualizzazione , accettare le impostazioni predefinite.
    NOTA: In alternativa, la dimensione totale dell'immagine può essere regolata inserendo un valore (in px) sotto la larghezza SVG. Il carattere e la spaziatura del testo possono essere regolati inserendo le dimensioni desiderate (in px) in Dimensione carattere SVG, Spazio tra le punte, Margine sinistro, Margine destro, Margine superiore e Margine inferiore.
  5. Fare clic su Dettagli di contatto. Inserisci l'indirizzo e-mail per il riferimento futuro dei risultati (facoltativo).
  6. Fare clic su Panoramica. Assicurati di rivedere il numero di taxa in input e gli outgroup selezionati.
  7. Fare clic su Invia. Fare clic con il pulsante destro del mouse sui file in File scaricabili (Figura 8B) e selezionare Salva con nome... per salvare i risultati.
    NOTA: Il file di testo di output indicherà il numero totale di specie delimitate all'inizio dei risultati. Presenterà quindi quali sequenze rientrano in ciascuno dei cluster. Questi cluster e le sequenze corrispondenti sono disposti in modo simile all'albero generato nel file .svg. Il file SVG mPTP visualizzerà il clustering delle sequenze. I cluster sono presentati da gruppi di sequenze collegate da rami rossi lungo l'albero. Ogni serie isolata di rami rossi indica un MOTU singolare.

figure-protocol-7
Figura 8: Delimitazione delle specie molecolari mediante processi multi-rate dell'albero di Poisson. (A) Finestre che mostrano le informazioni di lettura mPTP relative all'albero caricato. (B) Finestra risultante nel server web mPTP, con collegamento contenente i risultati della delimitazione (riquadro rosso). Clicca qui per visualizzare una versione più grande di questa figura.

8. Generazione MOTU

  1. Apri l'albero utilizzando un programma di editor di foto o utilizzando PowerPoint.
  2. Creare una barra per rappresentare i risultati di ciascun approccio di delimitazione delle specie molecolari (vedere la Figura 9 e la Figura 10 per i campioni). Assicurati che ogni singola barra sia allineata con le sequenze che include nel suo cluster molecolare.
  3. Determinare il MOTU con uno di questi due approcci74
    1. Se tutti gli approcci generano risultati identici per un dato cluster molecolare, erigere il cluster molecolare come MOTU per consenso.
    2. Se la maggior parte degli approcci genera risultati identici per un dato cluster molecolare, erigere il cluster molecolare come MOTU a maggioranza.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Risultati

Questo articolo ha implementato una pipeline di intruglio per la delimitazione delle specie molecolari, che include TaxonDNA, K2P, ASAP, due varianti di PTP e mPTP, per coleotteri acquatici e ripariali. Il risultato di ciascuno dei sei approcci è indicato come "cluster molecolari", mentre il riassunto è indicato come "unità tassonomiche operative molecolari" o MOTU75.

La pipeline di intruglio è stata testata su due set di dati. Il primo dataset (Tabella 1,

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Discussione

Mentre le sequenze di DNA, in particolare il COI49, possono essere confrontate con le sequenze di riferimento 22,53,78 negli archivi pubblici online 54,55,79,80 per un sistema di bioidentificazione46, tale approccio è limitato dalle sequenze depositate disponibili.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Dichiarazioni

Non abbiamo nulla da rivelare.

Ringraziamenti

Ringraziamo il Bureau of Fisheries and Aquatic Resources (BFAR) e le unità governative locali interessate per averci concesso il consenso informato preventivo per eseguire il campionamento della biodiversità. Ringraziamo il dottor Hendrik Freitag per la guida sul lavoro tassonomico. Ringraziamo anche i revisori e gli editori per i loro commenti utili e costruttivi. Infine, esprimiamo la nostra massima gratitudine all'University Research Council Big Project Grant (URC 2022-09) e all'Open Access Publication Grant dell'Ateneo de Manila University, al SC Johnson Student Research and Development Fund for Environmental Leadership dell'Ateneo Research Institute of Science and Engineering, al LinnéSys: Systematics Research Fund della Linnean Society di Londra e alla Systematics Association per aver finanziato questo studio.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Materiali

Elenco dei materiali utilizzati in questo articolo
NomeAziendaNumero di catalogoCommenti
AL PIÙ PRESTO Esploratore di partihttps://bioinfo.mnhn.fr/abi/public/asap/asapweb.html
ComputerQualsiasi unità di elaborazione
MEGA XMEGA Softwarehttps://www.megasoftware.net/
Server web mPTPhttps://mptp.h-its.org/#/tree
Server web PTPIstituto di studi teorici di Heidelberg (HITS)https://species.h-its.org/ptp/
Identificatore di specie   https://taxondna.sourceforge.net/

Riferimenti

  1. Ceballos, G., Ehrlich, P. R. Mammal population losses and the extinction crisis. Science. 296 (5569), 904-907 (2002).
  2. Dirzo, R., et al. Defaunation in the Anthropocene. Science. 345 (6195), 401-406 (2014).
  3. Cardoso, P., et al. Scientists' warning to humanity on insect extinctions. Biol Conserv. 242, 108426(2020).
  4. Wagner, D. L., Grames, E. M., Forister, M. L., Berenbaum, M. R., Stopak, D. Insect decline in the Anthropocene: Death by a thousand cuts. Proc Natl Acad Sci U S A. 118 (2), 1-10 (2021).
  5. Steffen, W., Crutzen, P. J., McNeill, J. R. Welcome to the Anthropocene. Ambio. 36 (8), 614-621 (2007).
  6. de Carvalho, M. R., et al. Taxonomic impediment or impediment to taxonomy? A commentary on systematics and the cybertaxonomic-automation paradigm. Evol Biol. 34 (3-4), 140-143 (2007).
  7. Satler, J. D., Carstens, B. C., Hedin, M. Multilocus species delimitation in a complex of morphologically conserved trapdoor spiders (Mygalomorphae, Antrodiaetidae, Aliatypus). Syst Biol. 62 (6), 805-823 (2013).
  8. Arribas, P., Andújar, C., Sánchez-Fernández, D., Abellán, P., Millán, A. Integrative taxonomy and conservation of cryptic beetles in the Mediterranean region (Hydrophilidae). Zool Scr. 42 (2), 182-200 (2013).
  9. Baker, K., et al. Carpooling with ecologists, geographers and taxonomists: Perceptions from conducting environmental research in tropical regions. Biodivers Conserv. 28 (4), 975-981 (2019).
  10. Tan, D. S. H., Ang, Y., Lim, G. S., Bin Ismail, M. R., Meier, R. From "cryptic species" to integrative taxonomy: An iterative process involving DNA sequences, morphology, and behaviour leads to the resurrection of Sepsis pyrrhosoma (Sepsidae: Diptera). Zool Scr. 39 (1), 51-61 (2010).
  11. Tan, M. K., Storozhenko, S. Y., Hwang, W. S., Meier, R. Integrative taxonomy reveals two sympatric species of the genus Eucriotettix Hebard, 1930 (Orthoptera: Tetrigidae). Zootaxa. 4268 (3), 377-394 (2017).
  12. Yule, C. M. Freshwater invertebrates. Freshwater Invertebrates of the Malaysian Region. , 23-31 (2004).
  13. Komarek, A., Jäch, M. A., Ji, L. Hydrophilidae: I. Check list and key to Palearctic and Oriental genera of aquatic Hydrophilidae (Coleoptera). Water Beetles of China. 3, 383-395 (2003).
  14. Rudoy, A., Ribera, I. The macroevolution of size and complexity in insect male genitalia. PeerJ. 4, e1882(2016).
  15. Rudoy, A., Beutel, R. G., Ribera, I. Evolution of the male genitalia in the genus Limnebius Leach, 1815 (Coleoptera, Hydraenidae). Zool J Linn Soc. 178, 97-127 (2016).
  16. Shapiro, A. The lock-and-key hypothesis: Evolutionary and biosystematic interpretation of insect genitalia. Annu Rev Entomol. 34 (1), 231-245 (1989).
  17. Scotland, R., Hughes, C., Bailey, D., Wortley, A. The big machine and the much-maligned taxonomist. Syst Biodivers. 1 (2), 139-143 (2003).
  18. Vinarski, M. V. Roots of the taxonomic impediment: Is the "integrativeness" a remedy. Integr Zool. 15 (1), 2-15 (2020).
  19. Cardoso, P., Erwin, T. L., Borges, P. A. V., New, T. R. The seven impediments in invertebrate conservation and how to overcome them. Biol Conserv. 144 (11), 2647-2655 (2011).
  20. Engel, M. S., et al. The taxonomic impediment: A shortage of taxonomists, not the lack of technical approaches. Zool J Linn Soc. XX, 1-7 (2021).
  21. Tancoigne, E., Dubois, A. Cladistics taxonomy: No decline, but inertia. Cladistics. , 1-4 (2013).
  22. Riedel, A., Sagata, K., Suhardjono, Y. R., Tänzler, R., Balke, M. Integrative taxonomy on the fast track: Towards more sustainability in biodiversity research. Front Zool. 10 (15), 1-9 (2013).
  23. Dayrat, B. Towards integrative taxonomy. Biol J Linn Soc. 85 (3), 407-415 (2005).
  24. Bickford, D., et al. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 22 (3), 148-155 (2007).
  25. Page, R. D. M. DNA barcoding and taxonomy: Dark taxa and dark texts. Philos Trans R Soc Lond B Biol Sci. 371 (1702), 20150334(2016).
  26. Morinière, J., et al. A DNA barcode library for 5,200 German flies and midges (Insecta: Diptera) and its implications for metabarcoding-based biomonitoring. Mol Ecol Resour. 19 (4), 900-928 (2019).
  27. Hartop, E., Srivathsan, A., Ronquist, F., Meier, R. Towards large-scale integrative taxonomy (LIT): Resolving the data conundrum for dark taxa. Syst Biol. 71 (6), 1404-1422 (2022).
  28. Fossen, E. I., Ekrem, T., Nilsson, A. N., Bergsten, J. Species delimitation in northern European water scavenger beetles of the genus Hydrobius (Coleoptera, Hydrophilidae). ZooKeys. 564, 71-120 (2016).
  29. Bálint, M., et al. Cryptic biodiversity loss linked to global climate change. Nat Clim Change. 1 (6), 313-318 (2011).
  30. Padial, J. M., De La Riva, I. A response to recent proposals for integrative taxonomy. Biol J Linn Soc. 101 (3), 747-756 (2010).
  31. Andújar, C., Arribas, P., Ruiz, C., Serrano, J., Gómez-Zurita, J. Integration of conflict into integrative taxonomy: Fitting hybridization in species delimitation of Mesocarabus (Coleoptera: Carabidae). Mol Ecol. 23 (17), 4344-4361 (2014).
  32. Schlick-Steiner, B. C., et al. Integrative taxonomy: A multisource approach to exploring biodiversity. Annu Rev Entomol. 55 (1), 421-438 (2010).
  33. Stüben, P. E., Astrin, J. J. Integrative taxonomy, phylogeny, and new species of the weevil genus Onyxacalles Stüben (Coleoptera: Curculionidae: Cryptorhynchinae). Psyche (Lond). 2012, 1-22 (2012).
  34. Hawlitschek, O., et al. Pleistocene climate change promoted rapid diversification of aquatic invertebrates in Southeast Australia. BMC Evol Biol. 12 (142), 1-15 (2012).
  35. Delocado, E. D., Freitag, H. Two new species of Byrrhinus Motschulsky, 1858 (Coleoptera, Limnichidae, Limnichinae) from Negros, Philippines. ZooKeys. 1070, 51-72 (2021).
  36. De Vera, G. C. M., Sia, J. I. L., Freitag, H., Delocado, E. D. The first new species of Leptelmis Sharp, 1888 (Coleoptera, Elmidae) from the Philippines in the last 50 years. Tijdschr Entomol. 166 (2-3), 175-186 (2023).
  37. Seno, C. B. N., Delocado, E. D., Freitag, H. Three new species of Ancyronyx Erichson, 1847 (Coleoptera, Elmidae) from Mindanao, Philippines. Tijdschr Entomol. 72, 61-80 (2022).
  38. López-López, A., Hudson, P., Galián, J. Islands in the desert: Species delimitation and evolutionary history of Pseudotetracha tiger beetles (Coleoptera: Cicindelidae: Megacephalini) from Australian salt lakes. Mol Phylogenet Evol. 101, 279-285 (2016).
  39. Ibarra, J. W. B., Bantiding, F. J. A., Delocado, E. D. First Georissidae (Coleoptera) from the Philippines with description of two new species of Georissus Latreille, 1809. Tijdschr Entomol. , 45-58 (2025).
  40. Boonsoong, B., Sartori, M. Review and integrative taxonomy of the genus Prosopistoma Latreille, 1833 (Ephemeroptera, Prosopistomatidae) in Thailand, with description of a new species. ZooKeys. 825, 123-144 (2019).
  41. Wang, Y., Nansen, C., Zhang, Y. Integrative insect taxonomy based on morphology, mitochondrial DNA, and hyperspectral reflectance profiling. Zool J Linn Soc. 177 (2), 378-394 (2016).
  42. Schutze, M. K., et al. One and the same: Integrative taxonomic evidence that Bactrocera invadens (Diptera: Tephritidae) is the same species as the Oriental fruit fly Bactrocera dorsalis. Syst Entomol. 40 (2), 472-486 (2015).
  43. Schär, S., et al. Integrative biodiversity inventory of ants from a Sicilian archipelago reveals high diversity on young volcanic islands (Hymenoptera: Formicidae). Org Divers Evol. 20, 405-416 (2020).
  44. Ševčík, J., et al. Integrative taxonomy of Central European Planetella (Diptera: Cecidomyiidae) indicates high species diversity, intraspecific variation and low host specificity. Acta Entomol Musei Natl Pragae. 63 (2), 413-450 (2023).
  45. Mengual, X., Ståhls, G., Vujić, A., Marcos-García, M. Á Integrative taxonomy of Iberian Merodon species (Diptera, Syrphidae). Zootaxa. 1377, 1-26 (2006).
  46. Hebert, P. D. N., et al. Biological identifications through DNA barcodes. Proc R Soc B Biol Sci. 270 (1512), 313-321 (2003).
  47. Hebert, P. D. N., et al. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc Natl Acad Sci U S A. 101 (41), 14812-14817 (2004).
  48. Hebert, P. D. N., Gregory, T. R. The promise of DNA barcoding for taxonomy. Syst Biol. 54 (5), 852-859 (2005).
  49. Hebert, P. D. N., Ratnasingham, S., DeWaard, J. R. Barcoding animal life: Cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc Lond B. 270 (Suppl), S96-S99 (2003).
  50. Simon, C., et al. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 87 (6), 651-701 (1994).
  51. Hendrich, L., et al. Mitochondrial Cox1 sequence data reliably uncover patterns of insect diversity but suffer from high lineage-idiosyncratic error rates. PLoS ONE. 5 (12), e14448(2010).
  52. Balke, M., et al. Suggestions for a molecular biodiversity assessment of Southeast Asian freshwater invertebrates: Lessons from the megadiverse beetles (Coleoptera). J Limnol. 72 (Suppl 2), 61-68 (2013).
  53. Virgilio, M., et al. Comparative performances of DNA barcoding across insect orders. BMC Bioinformatics. 11, 206(2010).
  54. Benson, D. A., et al. GenBank. Nucleic Acids Res. 41 (D1), 36-42 (2013).
  55. Ratnasingham, S., Hebert, P. D. N. BOLD: The barcode of life data system. Mol Ecol Notes. 7, 355-364 (2007).
  56. Pons, J., et al. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 55 (4), 595-609 (2006).
  57. Puillandre, N., et al. automatic barcode gap discovery for primary species delimitation. Mol Ecol. 21 (8), 1864-1877 (2012).
  58. Fontaneto, D., Flot, J. F., Tang, C. Q. Guidelines for DNA taxonomy, with a focus on the meiofauna. Mar Biodivers. 45 (3), 433-451 (2015).
  59. Balke, M., et al. New Guinea highland origin of a widespread arthropod supertramp. Proc R Soc B Biol Sci. 276 (1666), 2359-2367 (2009).
  60. Meier, R., et al. DNA barcoding and taxonomy in Diptera: A tale of high intraspecific variability and low identification success. Syst Biol. 55 (5), 715-728 (2006).
  61. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 16 (2), 111-120 (1980).
  62. Puillandre, N., Brouillet, S., Achaz, G. ASAP: Assemble species by automatic partitioning. Mol Ecol Resour. 21 (2), 609-620 (2021).
  63. Wiens, J. J., Penkrot, T. A. Delimiting species using DNA and morphological variation and discordant species limits in spiny lizards (Sceloporus). Syst Biol. 51 (1), 69-91 (2002).
  64. 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. 62 (5), 707-724 (2013).
  65. Yang, Z., Rannala, B. Unguided species delimitation using DNA sequence data from multiple loci. Mol Biol Evol. 31 (12), 3125-3135 (2014).
  66. Leaché, A. D., Fujita, M. K. Bayesian species delimitation in West African forest geckos (Hemidactylus fasciatus). Proc R Soc B Biol Sci. 277 (1697), 3071-3077 (2010).
  67. Luo, A., et al. Comparison of methods for molecular species delimitation across a range of speciation scenarios. Syst Biol. 67 (5), 830-846 (2018).
  68. Zhang, J., et al. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 29 (22), 2869-2876 (2013).
  69. Kapli, P., et al. Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo. Bioinformatics. 33 (11), 1630-1638 (2017).
  70. Sanchez, E. G. S., Delocado, E. D., Freitag, H. Two new species of Anacaena Thomson, 1859 (Coleoptera, Hydrophilidae) from the Philippines. ZooKeys. 1112, 11-25 (2022).
  71. Sanger, F., Coulson, A. R. A rapid method for determining sequences in DNA by primed synthesis with DNA polymerase. J Mol Biol. 94 (3), 441-448 (1975).
  72. Kumar, S., et al. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 35 (6), 1547-1549 (2018).
  73. Thompson, J. D., Higgins, D. G., Gibson, T. J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673-4680 (1994).
  74. Pentinsaari, M., Vos, R., Mutanen, M. Algorithmic single-locus species delimitation: Effects of sampling effort, variation and nonmonophyly in four methods and 1870 species of beetles. Mol Ecol Resour. 17 (3), 393-404 (2017).
  75. Galimberti, A., et al. Integrated operational taxonomic units (IOTUs) in echolocating bats: A bridge between molecular and traditional taxonomy. PLoS ONE. 7 (6), e40122(2012).
  76. Kundrata, R., Jäch, M. A., Bocak, L. Molecular phylogeny of the Byrrhoidea-Buprestoidea complex (Coleoptera, Elateriformia). Zool Scr. 46 (2), 150-164 (2017).
  77. Komarek, A., Freitag, H. Taxonomic revision of Agraphydrus Règimbart, 1903 V. Philippine species and their first DNA barcodes. Koleopt Rundsch. 90, 201-242 (2020).
  78. Tautz, D., et al. A plea for DNA taxonomy. Trends Ecol Evol. 18 (2), 70-74 (2003).
  79. Altschul, S. F., et al. Basic local alignment search tool. J Mol Biol. 215 (3), 403-410 (1990).
  80. Porter, T., Hajibabaei, M. Over 2.5 million COI sequences in GenBank and growing. bioRxiv. , 353904(2018).
  81. Gómez, R. A., Reddell, J., Will, K., Moore, W. Up high and down low: Molecular systematics and insight into the diversification of the ground beetle genus Rhadine LeConte. Mol Phylogenet Evol. 98, 161-175 (2016).
  82. Kodada, J., et al. Ancyronyx clisteri, a new spider riffle beetle species from Borneo, redescription of A. sarawacensis Jäch including a description of the larva and new distribution data for A. procerus Jäch using DNA barcodes (Coleoptera, Elmidae). ZooKeys. 912, 25-64 (2020).
  83. Hajibabaei, M., Janzen, D. H., Burns, J. M., Hallwachs, W., Hebert, P. D. N. DNA barcodes distinguish species of tropical Lepidoptera. Proc Natl Acad Sci U S A. 103 (4), 968-971 (2006).
  84. Yeo, D., Srivathsan, A., Meier, R. Longer is not always better: Optimizing barcode length for large-scale species discovery and identification. Syst Biol. 69 (5), 999-1015 (2020).
  85. Sabordo, M. R., Delocado, E., Freitag, H. Two new species of the genus Ancyronyx Erichson, 1847 from the island of Negros, Philippines (Insecta, Coleoptera, Elmidae). Tijdschr Entomol. 163 (1), 13-30 (2020).
  86. Astrin, J. J., et al. Exploring diversity in cryptorhynchine weevils (Coleoptera) using distance-, character- and tree-based species delineation. Mol Phylogenet Evol. 63 (1), 1-14 (2012).
  87. Pritchard, J. K., Stephens, M., Donnelly, P. Inference of population structure using multilocus genotype data. Genetics. 155, 945-959 (2000).
  88. Huelsenbeck, J. P., Andolfatto, P., Huelsenbeck, E. T. Structurama: Bayesian inference of population structure. Evol Bioinform. 7, 55-59 (2011).
  89. Hausdorf, B., Hennig, C. Species delimitation using dominant and codominant multilocus markers. Syst Biol. 59 (5), 491-503 (2010).
  90. Fujita, M. K., et al. Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol Evol. 27 (9), 480-488 (2012).
  91. Monaghan, M. T., et al. Accelerated species inventory on Madagascar using coalescent-based models of species delineation. Syst Biol. 58 (3), 298-311 (2009).
  92. Yang, Z., Rannala, B. Bayesian species delimitation using multilocus sequence data. Proc Natl Acad Sci U S A. 107 (20), 9264-9269 (2010).
  93. Ence, D. D., Carstens, B. C. SpedeSTEM: A rapid and accurate method for species delimitation. Mol Ecol Resour. 11 (3), 473-480 (2011).
  94. de Queiroz, K. Species concepts and species delimitation. Syst Biol. 56 (6), 879-886 (2007).
  95. Blair, C., Bryson, R. W. Cryptic diversity and discordance in single-locus species delimitation methods within horned lizards (Phrynosomatidae: Phrynosoma). Mol Ecol Resour. 17 (6), 1168-1182 (2017).
  96. Dellicour, S., Flot, J. F. Delimiting species-poor data sets using single molecular markers: A study of barcode gaps, haplowebs and GMYC. Syst Biol. 64 (6), 900-908 (2015).
  97. Dellicour, S., Flot, J. F. The hitchhiker's guide to single-locus species delimitation. Mol Ecol Resour. 18 (6), 1234-1246 (2018).
  98. Liu, L., et al. Estimating phylogenetic trees from genome-scale data. Ann N Y Acad Sci. 1360 (1), 36-53 (2015).
  99. Lim, G. S., Balke, M., Meier, R. Determining species boundaries in a world full of rarity: Singletons, species delimitation methods. Syst Biol. 61 (1), 165-169 (2012).
  100. Zhou, Z., et al. molecular species delimitation based on COI-5P barcode sequences revealed high cryptic/undescribed diversity for Chinese katydids (Orthoptera: Tettigoniidae). BMC Evol Biol. 19 (1), (2019).
  101. Toussaint, E. F. A., et al. Comparative molecular species delimitation in the charismatic Nawab butterflies (Nymphalidae, Charaxinae, Polyura). Mol Phylogenet Evol. 91, 194-209 (2015).
  102. Carstens, B. C., Dewey, T. A. Species delimitation using a combined coalescent and information-theoretic approach: An example from North American Myotis bats. Syst Biol. 59 (4), 400-414 (2010).
  103. Esselstyn, J. A., et al. Single-locus species delimitation: A test of the mixed Yule-coalescent model, with an empirical application to Philippine round-leaf bats. Proc R Soc Lond B. 279 (1743), 3678-3686 (2012).
  104. Stech, M., et al. Molecular species delimitation in the Racomitrium canescens complex (Grimmiaceae) and implications for DNA barcoding of species complexes in mosses. PLoS One. 8 (1), e53134(2013).
  105. Kekkonen, M., Hebert, P. D. N. DNA barcode-based delineation of putative species: Efficient start for taxonomic workflows. Mol Ecol Resour. 14 (4), 706-715 (2014).
  106. Meier, R., et al. "Dark taxonomy": A new protocol for overcoming the taxonomic impediments for dark taxa and broadening the taxon base for biodiversity assessment. Cladistics. 41, 223-238 (2025).
  107. Ebach, M. C., Holdrege, C. DNA barcoding is no substitute for taxonomy. Nature. 434, 697(2005).
  108. Ebach, M. C., Holdrege, C. More taxonomy, not DNA barcoding. BioScience. 55 (10), 822-823 (2005).
  109. Mora, C., et al. How many species are there on Earth and in the ocean. PLoS Biol. 9 (8), e1001127(2011).
  110. Stork, N. E. How many species of insects and other terrestrial arthropods are there on Earth. Annu Rev Entomol. 63, 31-45 (2018).
  111. Scheffers, B. R., et al. What we know and don't know about Earth's missing biodiversity. Trends Ecol Evol. 27 (9), 501-510 (2012).
  112. Carbayo, F., Marques, A. C. The cost of describing the entire animal kingdom. Trends Ecol Evol. 26 (4), 153-154 (2011).
  113. Walsh, J. C., et al. Trends and biases in the listing and recovery planning for threatened species: An Australian case study. Oryx. 47 (1), 134-143 (2013).
  114. Carrizo, S. F., et al. Freshwater megafauna: Flagships for freshwater biodiversity under threat. BioScience. 67 (10), 919-927 (2017).

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Ristampe e permessi

Tag

Pipeline MOTUDNA barcodingdelimitazione delle speciecoleotteri riparialisequenze COIallineamento delle sequenzealbero di massima verosimiglianzamodello di Kimura a 2 parametri