Neighbor-Joining-Algorithmus

Neighbor-Joining-Algorithmus

Der Neighbor-Joining-Algorithmus ist ein mathematisches Verfahren, um Datensätze zu vergleichen und hierarchisch bifurcal (zweigabelig) anzuordnen. Dieses Verfahren wurde 1987 von Saitou und Nei vorgestellt und 1988 von Studier und Keppler weiterentwickelt und vereinfacht.

Inhaltsverzeichnis

Anwendung

In der Bioinformatik bezeichnet das Neighbor-Joining-Verfahren eine phänetische bottom-up Clustermethode, welche zur Erstellung von phylogenetischen Baumstrukturen verwendet wird. Hiermit soll anhand von variierenden Merkmalen in der Datenmatrix die Wahrscheinlichkeit einer Abstammungs- oder Verwandtschafts-Beziehung in einer Stammbaumartigen Darstellung berechnet werden.

Normalerweise werden damit Bäume aus DNA- oder Proteinsequenzdaten oder klassisch morphologischen Datensätzen erstellt. Der Algorithmus benötigt Wissen über die Distanz zwischen zwei Paaren von Taxa (also beispielsweise Arten oder Sequenzen) in einem Baum.

Algorithmus

Neighbor-joining basiert meist auf dem „Minimum Evolution“-Kriterium für phylogenetische Bäume: Ausgehend von einem zunächst sternförmigen „Baum“, in dem alle Taxa mit einem „Zentrum“ verbunden sind, werden paarweise die DNA- oder Proteinsequenzen mit der geringsten genetischen Distanz ausgewählt und zu einem Ast des Baumes vereinigt. Die genetischen Distanzen der Sequenzen werden neu berechnet und wieder die nächstverwandten zu einem Ast mit zwei Taxa zusammengefügt. Dies erfolgt solange, bis alle Taxa in dem Baum eingefügt wurden und die Sternstruktur des Baumes völlig aufgelöst wurde. Im Unterschied zum UPGMA berücksichtigt Neighbor-Joining, dass die Evolutionsgeschwindigkeit nicht konstant ist: Wenn ein Taxon von allen anderen Taxa weit entfernt ist, so ist dies nicht auf einen entfernten Verwandtschaftsgrad, sondern auf beschleunigte Evolution zurückzuführen.

Beispiel

[1]

Folgend ist eine typische Tabelle von Distanzen zwischen Taxa angegeben, wobei die Werte rein hypothetisch aber realistisch sind:

Mensch Maus Rose Tulpe
Mensch 0 3 14 12
Maus 3 0 13 11
Rose 14 13 0 4
Tulpe 12 11 4 0

Da die Tabelle dreiecks-symmetrisch ist, muss die untere Hälfte nicht unbedingt gespeichert werden. Die Werte in dieser Tabelle werden als di,j benannt.

Schritt 1: Es müssen die Durchschnittlichen Distanzen von jedem Taxon zu jedem anderen berechnet werden. Dies geschieht mit folgender Formel für die Netto-Divergenz ri:

r_i=\frac{1}{N-2}\sum_{k=1}^Nd_{i,k}

Wobei N die Anzahl der Taxa ist.

Mensch r_1 = \frac{0 + 3 + 14 + 12}{4-2} = 14{,}5
Maus r_2 = \frac{3 + 0 + 13 + 11}{4-2} = 13{,}5
Rose r_3 = \frac{14 + 13 + 0 + 4}{4-2} = 15{,}5
Tulpe r_4 = \frac{12 + 11 + 4 + 0}{4-2} = 13{,}5

Interpretation: Die Rose besitzt in unserem Beispiel die größte Netto-Divergenz, hat also im Vergleich mit den anderen Taxa eine größere Evolutionsgeschwindigkeit durchlebt.

Schritt 2: Wir berechnen eine Zwischenmatrix M.

Mi,j = di,j − (ri + rj)

Wie z.B. zwischen Mensch und Maus:

M1,2 = d1,2 − (r1 + r2) = 3 − (14,5 + 13,5) = − 25

M Mensch Maus Rose Tulpe
Mensch -25 -16 -16
Maus -25 -16 -16
Rose -16 -16 -25
Tulpe -16 -16 -25

Schritt 3: In dieser neu berechneten Distanzmatrix M wird nun der kleinste Wert, also die kleinste Distanz zwischen zwei Taxa, gesucht, und die gefundenen zwei Taxa werden zu einem neuen Teilbaum u = (i,j) zusammengefügt. In diesem Beispiel ergeben sich also die zwei Möglichkeiten Mensch und Maus, oder Rose und Tulpe zu einem Teilbaum zusammenzufügen. Wir entscheiden uns zunächst für Mensch und Maus.

Die Kantenlänge des Knotens zu der Verzweigung berechnet sich wie folgt:

v_{i,u} = \frac{d_{i,j} + r_i -r_j}{2}

vj,u = di,jvi,u

Also Mensch zu MeMa ist gleich  \frac{3 + 14,5 - 13,5}{2} = 2

Schritt 4: In der ursprünglichen Distanzmatrix wird der neue Eintrag u = MeMa angefügt:

Mensch Maus Rose Tulpe MeMa
Mensch 0 3 14 12  ?
Maus 3 0 13 11  ?
Rose 14 13 0 4  ?
Tulpe 12 11 4 0  ?
MeMa  ?  ?  ?  ? 0

Um die Distanzen des neuen Eintrages u = (i,j) = (1,2) = (Mensch, Maus) = MeMa zu den restlichen Taxa zu berechnen, wird folgende Formel verwendet:

d_{u,k}=\frac{d_{i,k}+d_{j,k}-d_{i,j}}{2}

Wobei die Einträge i und j zu einem neuen Eintrag u zusammengefügt wird, und die Distanz zum Eintrag k ausgerechnet wird. Die Distanz zwischen Rose und dem neuen Teilbaum ist also:

d_{5,3} = \frac{d_{1,3} + d_{2,3} - d_{1,2}}{2} = \frac{14 + 13 - 3}{2} = 12

Die „alten“, zusammengefügten Einträge, werden aus den Distanzmatrixen gelöscht.

Rose Tulpe MeMa
Rose 0 4 12
Tulpe 4 0 10
MeMa 12 10 0

Danach werden wieder ri und Mi,j berechnet, neu zusammengefügt und wieder von vorne angefangen. Dies wird solange wiederholt, bis nur noch zwei Taxa übrig bleiben, die dann schlussendlich verbunden werden.

Das Ergebnis unseres Beispiels lässt sich als additiven Baum darstellen:

Mensch                  Rose
   \                     /
    \ 2               3 /
     \        9        /
      -----------------
     /                 \
    / 1               1 \
   /                     \
 Maus                   Tulpe

oder die Ausgabe des Phylip Programms

 +----Maus
 !
 !                                         +-------------Rose
 1-----------------------------------------2
 !                                         +----Tulpe
 !
 +--------Mensch

Einordnung

Neighbor-Joining gehört zu den expliziten Methoden. Dies bedeutet, dass bei der Berechnung der genetischen Distanzen unterschiedliche Evolutionsmodelle, d. h. unterschiedliche Wahrscheinlichkeiten für Punktmutationen angenommen werden können. Die Richtigkeit dieser Stammbäume beruht auf der Annahme, dass die Veränderung der betrachteten Merkmale keine unbekannten Zwischenschritte enthält. Es wird also vereinfacht angenommen, dass "die Evolution keine Umwege geht" ("minimum evolution").

Der Neighbor-Joining-Algorithmus berechnet den Stammbaum schrittweise und findet deshalb nicht zwangsläufig die optimale Baum-Topologie mit der geringsten Verzweigungslänge. Dies beruht auf seinem Konstruktionsprinzip, als Greedy Algorithmus. Im Gegensatz zu anderen Algorithmen berechnet dieser nicht alle möglichen Bäume und wählt zum Schluss die optimalen aus, sondern verwirft schon während des Verfahrens einige Rechenwege. Obwohl der Algorithmus suboptimal ist, wurde er ausführlich getestet und findet normalerweise einen Baum, der dem Optimum relativ nahe kommt.

Vorteile

Der größte Vorteil dieses Verfahrens ist seine Geschwindigkeit. Man kann es auf gewaltige Datenmengen anwenden, selbst dort, wo andere Methoden der phylogenetischen Analyse wie maximum parsimony und Maximum-Likelihood) nicht mehr durchführbar sind. Im Gegensatz zum UPGMA-Algorithmus (Unweighted Pair Group Method with Arithmetic mean) zur phylogenetischen Baumrekonstruktion nimmt Neighbor-Joining nicht an, dass die Entwicklung der Abstammungslinien mit derselben Rate (siehe auch Molekulare Uhr) stattfindet und erzeugt daher infolgedessen einen unbalancierten Baum.

Literatur

  • N. Saitou and M. Nei : The neighbor-joining method: a new method for reconstructing phylogenetic trees. In: Molecular Biology and Evolution. 1987, Vol 4(4), p406-425. (full text).
  • J.A. Studier and K.J. Keppler : A note on the neighbor-joining algorithm of Saitou and Nei. In: Molecular Biology and Evolution,1988, Vol 5(6), p729-731. (full text).
  • Knoop, Volker; Müller, Kai: Gene und Stammbäume - Ein Handbuch zur molekularen Phylogenetik. Spektrum Akademischer Verlag, 2006. ISBN 3-8274-1642-6

Quellen

  1. http://www.icp.ucl.ac.be/~opperd/private/neighbor.html

Weblinks


Wikimedia Foundation.

Игры ⚽ Нужно сделать НИР?

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Neighbor-Joining-Verfahren — In diesem Artikel oder Abschnitt fehlen folgende wichtige Informationen: Hier müssen einfach noch viel mehr Infos rein. Das Beispiel muss noch weitergeführt werden. Bilder sollten das Ganze verständlich darstellen (ungewurzelte additive Bäume… …   Deutsch Wikipedia

  • Neighbour-Joining-Algorithmen — In diesem Artikel oder Abschnitt fehlen folgende wichtige Informationen: Hier müssen einfach noch viel mehr Infos rein. Das Beispiel muss noch weitergeführt werden. Bilder sollten das Ganze verständlich darstellen (ungewurzelte additive Bäume… …   Deutsch Wikipedia

  • Neighbour-joining — In diesem Artikel oder Abschnitt fehlen folgende wichtige Informationen: Hier müssen einfach noch viel mehr Infos rein. Das Beispiel muss noch weitergeführt werden. Bilder sollten das Ganze verständlich darstellen (ungewurzelte additive Bäume… …   Deutsch Wikipedia

  • Neighbour joining — In diesem Artikel oder Abschnitt fehlen folgende wichtige Informationen: Hier müssen einfach noch viel mehr Infos rein. Das Beispiel muss noch weitergeführt werden. Bilder sollten das Ganze verständlich darstellen (ungewurzelte additive Bäume… …   Deutsch Wikipedia

  • ClustalW — Clustal Entwickler: Gibson T. (EMBL), Thompson J. (CNRS), Higgins D. (UCD) Aktuelle Version: 2.10 (14. Oktober 2008) Betriebssystem …   Deutsch Wikipedia

  • ClustalX — Clustal Entwickler: Gibson T. (EMBL), Thompson J. (CNRS), Higgins D. (UCD) Aktuelle Version: 2.10 (14. Oktober 2008) Betriebssystem …   Deutsch Wikipedia

  • Clustal W — Clustal Entwickler: Gibson T. (EMBL), Thompson J. (CNRS), Higgins D. (UCD) Aktuelle Version: 2.10 (14. Oktober 2008) Betriebssystem …   Deutsch Wikipedia

  • Clustal — Entwickler Gibson T. (EMBL), Thompson J. (CNRS), Higgins D. (UCD) Aktuelle Version 2.1 (17. November 2010) Betriebssystem Unix, Linux, Mac OS X, Microsoft Win …   Deutsch Wikipedia

  • Phylogenetischer Baum — basierend auf rRNA Genen Ein phylogenetischer Baum ist ein Baum, der die evolutionären Beziehungen zwischen verschiedenen Arten oder anderen Einheiten, von denen man vermutet, dass sie einen gemeinsamen Vorfahren besitzen, darstellt. Damit ist… …   Deutsch Wikipedia

  • Biologische Systematik — Systematik (v. altgriech. συστηματικός systēmatikós „geordnet“) oder Biosystematik ist ein Fachgebiet der Biologie. Die klassische Systematik beschäftigt sich hauptsächlich mit der Bestimmung und Benennung der Lebewesen (Taxonomie). Die moderne… …   Deutsch Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”