Meine Merkliste
my.bionity.com  
Login  

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 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 Neighbour-Joining, daß die Evolutionsgeschwindigkeit nicht konstant ist: Wenn ein Taxon von allen anderen Taxa weit entfernt ist, so ist dies nicht auf einen entfernten Verwandschaftsgrad, 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 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 zusammengefügt zu einem neuen Teilbaum u = (i,j). 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, daß "die Evolution keine Umwege geht" ("minimum evolution").

Der Neighbour-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
 
Dieser Artikel basiert auf dem Artikel Neighbor-Joining-Algorithmus aus der freien Enzyklopädie Wikipedia und steht unter der GNU-Lizenz für freie Dokumentation. In der Wikipedia ist eine Liste der Autoren verfügbar.
Ihr Bowser ist nicht aktuell. Microsoft Internet Explorer 6.0 unterstützt einige Funktionen auf ie.DE nicht.