Posts mit dem Label Phylogeny werden angezeigt. Alle Posts anzeigen
Posts mit dem Label Phylogeny werden angezeigt. Alle Posts anzeigen

Dienstag, 10. Mai 2016

Optimierung von phylogenetischen Bäumen

Siehe auch...
...Einen Neighbor-Joining Baum erstellen...
...Phylogentische Analysen - Grundlagen von Bäumen...

Nachdem man einen phylogenetischen Baum mit z.B. MEGA berechnet hat, stellt sich die Frage, was man mit dem Baum zeigen will. Phylogenetische Bäume lassen sich in der Regel unterschiedlich darstellen. Zum einen gibt es grundsätzliche Darstellungen (MEGA unterscheidet hier zwischen klassisch, radiär und kreisförmig) und zum anderen ist vielleicht nicht jede Information im Baum relevant für die aktuelle Fragestellung.

Abb. 1: NJ-Phylogeny von Bambus, basierend auf rbcL Sequenzen. Klassisches Phylogram. Rohfassung.
In Abbildung 1 sehen wir die direkte Ausgabe (Rohfassung) eines phylogenetischen Baumes (NJ) der mit 1000 Bootstrap Replikaten getestet wurde. An den Knotenpunkten stehen Prozent Angaben. Diese geben wieder, wie oft die nachfolgenden Einheiten bei den Baum-Replikaten in dieser Zusammensetzung vor kamen. Je höher die Zahl, desto sicherer ist es, dass diese Aufteilung von Bedeutung ist. Im dargestellten Phylogram wird Wert auf die Darstellung der genetischen Distanz (in Form von Veränderungen innerhalb des untersuchten Sequenz-Bereichs) gelegt. Deswegen findet man in der Abbildung auch einen Maßstab. Die Länge von Knotenpunkt zu Knotenpunkt bzw. zum Ende eines Zweiges ist relativ zur genetischen Distanz. Im angeführten Beispiel entspricht die Länge des Maßstabes 0.2 % genetischer Distanz. Bei einer Sequenz Länge von 500 Nucleotiden wäre das 1 Nucleotid, dass im Vergleich anders ist.
Will man mit seinem Baum besonders auf die genetische Distanz hinweisen, würde es sich anbieten, diese für die Zweige anzuzeigen und die Bootstrap Werte auszublenden. Das lässt sich bewerkstelligen indem man über das View-Menü die Optionen auswählt und im Branch-Tab die Anzeige der Statistik/Frequenz abwählt und statt dessen die Zweig-Längen anzeigen lässt.

Abb. 2: Branch Tab, Tree Options Menü von MEGA
Zur Veranschaulichung der 3 Hauptgruppen innerhalb der Bambus Gewächse habe ich diverse Gruppierungen zusammengeklappt.

Abb. 3: Modifizierter Baum aus Abb. 1 mit Schwerpunkt auf die Distanzen zwischen den Bambus Hauptgruppen Arundinarieae, Bambuseae und Olyreae.
Will man mehr auf die Gruppierungen eingehen, empfiehlt es sich den Baum als Cladogram darzustellen und Verzweigungen nur dann zu berücksichtigen, wenn die zugrunde liegenden Daten diese mit einer gewissen Sicherheit unterstützen. In MEGA lässt sich das erreichen indem man einen kondensierten Baum erstellt (Compute - Condensed Tree):

Abb. 4: Modifizierter (kondensierter) Baum aus Abb. 1 mit Schwerpunkt auf Gruppierungen (hier innerhalb der Arundinarieae) die mit mehr als 50% der Bootstrap Replikate unterstützt werden. 
Zuletzt will ich noch auf die Subtree Drawing Options hinweisen, mit denen man auf übergeordnete Gruppierungen hinweisen kann (Abb. 4, Arundinarieae).


Zusammen mit den Tree Options lässt sich so auf besondere Dinge hinweisen, die in der Rohfassung des Baumes nicht direkt ersichtlich sind.

Abb. 5: Modifizierter (kondensierter) Baum aus Abb. 1 mit Schwerpunkt auf Arundinarieae und Bambuseae einschließlich einer Blatt Charakteristik (40 X stereomikroskopische Aufnahme, getrocknertes Blatt, Unterseite).
Mehr Informationen, wie z.B. alternative Software für die Darstellung von phylogenetischen Bäumen findet man hier.

Grüße
Thomas

Samstag, 8. August 2015

Morphologie - Morphometrie

Die Evolution hat viele verschiedene morphologische Phänomene hervorgebracht, die wir sehen, messen und vergleichen können. Im folgenden werden Programme vorgestellt, mit denen man entsprechende Daten erzeugen und auswerten kann.

Übersicht

ProgrammFunktionDaten
ShapeErfassung und AuswertungKonturen
LeafAnalyserErfassung und AuswertungKonturen
LaminaErfassungDimensionen, Konturen
MorphoJAuswertungKonturen
Morphomatica?Konturen
Morpheus et al?Konturen
CLICErfassung und AuswertungKonturen

Shape

Iwata, H. & Ukai, Y., 2002. SHAPE: A Computer Program Package for Quantitative Evaluation of Biological Shapes Based on Elliptic Fourier Descriptors. Journal of Heredity, 93(5), pp.384–385.
Available at: http://jhered.oxfordjournals.org/content/93/5/384.short .
  • Es handelt sich hierbei um ein Programm zur Analyse von geschlossenen Konturen.
  • Die im Paper angegebene URL , welche zum Programm Download führen sollte, funktioniert nicht. Über eine Google Suche nach "SHAPE Iwata" gelangt man aber zur derzeitigen Web Seite  auf der man das Programm runter laden kann.

LeafAnalyser

Weight, C., Parnham, D. & Waites, R., 2008. TECHNICAL ADVANCE: LeafAnalyser: a computational method for rapid and large-scale analyses of leaf shape variation. The Plant Journal, 53(3), pp.578–586.
Available at: http://dx.doi.org/10.1111/j.1365-313X.2007.03330.x .
  • Die im Paper angegebene URL , welche zum Programm Download führen sollte, funktioniert nicht. Es hat den Anschein, als wäre die Seite gehackt worden. Man findet das Programm aber bei SourceForge.

Morphomatica

"Morphomatica ist ein benutzerfreundliches Computerprogramm, dass zur morphometrischen Analyse der Konturen von Ostrakodenschalen entworfen wurde. Es wurde im Zeitraum von 2001 – 2007 von J. Linhart, W. Brauneis und W. Neubauer (Universität Salzburg) unter Verwendung der biologischen Grundlagen von Dan. L. Danielopol (Österreichische Akademie der Wissenschaften) entwickelt. Die Version Morphomatica 1.6.1 eignet sich besonders für vergleichende Studien von merkmalsarmen Ostrakodenschalen."

Collection of Landmark for Identification and Characterization (CLIC)


Freitag, 29. Mai 2015

Warum verwenden wir gerade DEN Marker ?

Nachdem wir uns mit der allgemeinen Frage "Was ist eigentlich ein Marker?" auseinander gesetzt haben, wollen wir nun auf die Frage eingehen, warum man gerade den/die Marker verwendet und nicht andere. In anderen Worten, wir wollen wissen warum ein Marker besser dafür geeignet ist, Arten zu unterscheiden oder die Beziehungen bestimmter taxonomischer Einheiten zu untersuchen als ein anderer.

Marker in der DNA Diagnostik und phylogenetischen Untersuchungen

Für die Identifizierung taxonomischer Einheiten werden entweder bekannte Sequenz-Marker verwendet oder es werden über DNA Fingerprint Techniken neue (anonyme) DNA Bereiche als Marker herangezogen. Wir verwenden dafür bevorzugt bekannte Sequenz-Marker, da diese zum einen gut charakterisiert sind und zum anderen bereits Daten zur Verfügung stehen. Mit Hilfe dieser Informationen kann man sich im Vorfeld ein Bild darüber machen, ob der Bereich zur Unterscheidung der entsprechenden Art bzw. zum Vergleich entsprechender taxonomischer Einheiten geeignet ist oder eher nicht.

Was genau erwarten wir von unserem Marker ?

Wenn wir einen DNA Bereich verschiedener Pflanzen Arten amplifizieren wollen benötigen wir - für entsprechende Primer - Bereiche, die in all diesen Arten konserviert sind. Für die Unterscheidung der Arten benötigen wir einen Bereich, der dafür ausreichend Variation enthält.

Eine Aufstellung von DNA Bereichen (Zellkern, Chloroplast, Mitochondrium), die in verschiedenen Arbeiten für die Unterscheidung bzw. Identifizierung unterschiedlicher taxonomischer Gruppen (Population, Unterart, Art, Gattung, Familie, Ordnung) verwendet wurden. Aus Yip et al. 2007
Neben diesen grundsätzlichen Eigenschaften gibt es weitere, die im Zuge eines Projektes relevant werden können. Ein Beispiel bezieht sich auf den Schritt der Sequenzierung. Nicht jeder Bereiche lässt sich in allen taxonomischen Gruppen erfolgreich sequenzieren. Der Chloroplast Marker psbA-trnH lässt sich mit universellen Primern unter anderem in Dracocephalum problemlos amplifizieren und sequenzieren. In Dianthus dagegen lässt sich dieser Bereich zwar amplifizieren, aber bei der Sequenzierung stößt man auf das Problem von Mononucleotid Wiederholungen (poly-T).


Ein weiteres Beispiel betrifft den Bereich der ribosomalen DNA (rDNA) welcher im Kerngenom liegt und zwar in vielen Kopien an unterschiedlichen Stellen. In manchem Fällen, wie zum Beispiel bei einigen Rheum Arten [Xie et al. 2014], wurde festgestellt, dass diese Kopien sich voneinander unterscheiden, was dazu führt das man verschiedene Fragmente in einer PCR amplifiziert und eine einfache Sequenzierung der PCR nicht möglich ist.

Solche Regionen lassen sich in den entsprechenden taxonomischen Gruppen nur über Umwege und damit verbundenen Mehrkosten sequenzieren. Je nachdem um was es in einem Projekt geht, können solche "Probleme" aber auch eine zusätzliche Quelle für Informationen darstellen und ein Mehraufwand könnte sich lohnen.

Zusammenfassung

Bei der Marker Wahl haben wir zum einen grundsätzliche Faktoren zu beachten, die abhängig davon sind, welche Art von Beziehungen wir untersuchen wollen. Zum anderen müssen wir Faktoren beachten, die auf Erfahrung (vorhandene Messergebnisse - Untersuchungen) beruhen. Dies macht eine ausgiebige Literatur und Datenbank Recherche zu einer wichtigen Grundvoraussetzung für unsere Projekte.

Referenzen

Yip, P. Y., Chau, C. F., Mak, C. Y., & Kwan, H. S. (2007). DNA methods for identification of Chinese medicinal materials. Chinese Medicine, 2, 9. doi:10.1186/1749-8546-2-9

Ma, X., Xie, C., Guan, M., Xu, X., Miki, E., Takeda, O., … Chen, S. (2014). High Levels of Genetic Diversity within One Population of Rheum tanguticum on the Qinghai-Tibet Plateau have Implications for Germplasm Conservation. Pharmaceutical Crops, (5), 1–8.

Samstag, 2. Mai 2015

Einen Neighbor-Joining Baum erstellen...

Die Erstellung eines Neighbor-Joining Baumes lässt sich in drei Schritte gliedern:
  1. Erstellung eines Alignments
  2. Berechnung einer Distanz Matrix
  3. Erstellung eines Baumes
Alle Punkte können mit Hilfe von MEGA ausgeführt werden. Zur Erstellung des Alignments öffnet man eine entsprechende Sammlung von Sequenzen und lässt diese mit Clustal oder Muscle angleichen. Im Anschluss schaut man sich die Seitenbereiche an und schneidet diese entweder entsprechend zu, oder entfernt extrem kurze Sequenzen und schneidet anschließend zu. Hier muss man abwägen zwischen den Informationen, die verloren gehen wenn man die Ränder zuschneidet, und den Informationen die in einer entsprechend kurzen Sequenz stecken.
Die Schritte 2 und 3 werden von MEGA in einem durchgeführt. Wenn man mit dem Alignment zufrieden ist, exportiert man die Daten für eine phylogenetische Analyse.
Danach wechselt man zum Hauptfenster von MEGA und wählt unter Phylogeny "Construct/Test Neighbor-Joining Tree...".
Im nun offenen Fenster (Analysis Preferences) wählt man die gewünschten Einstellungen zur Berechnung der Distanz Matrix (siehe Phylogentische Analysen - Grundlagen von Bäumen) und definiert wie man die Phylogenie testen will. Zweiteres kann man mit gutem Gewissen mit einem Bootstrap Test (500 Replikate) durchführen.
Durch die Bestätigung der Einstellungen wird die Baumkonstruktion und der Test der Phylogenie gestartet.

Das Ergebnis lässt sich als MEGA Tree-Session speichern, um später weiter zu arbeiten, als Bild exportieren um es in einer Präsentation zu verwenden oder man kann an den Details der Darstellung arbeiten. Mehr dazu in einem anderen Post.

Thomas

Phylogentische Analysen - Grundlagen von Bäumen

Ein phylogenetischer Baum (Baum) repräsentiert Verwandtschaftsverhältnisse zwischen biologischen Einheiten basierend auf messbaren Eigenschaften. Anfänglich wurden in der numerischen Taxomonie morphologische Daten verwendet. Später folgten Sequenz-Informationen von Proteinen und DNA sowie Informationen aus DNA Fingerprint Analysen.

Die einfachste und gleichzeitig schnellste Methode, um einen Baum zu erstellen, ist mit Hilfe einer Distanz Matrix. Hierfür werden die einzelnen Sequenzen paarweise verglichen und jeweils ein Wert, der die Distanz beschreibt, festgehalten. In MEGA stehen Neighbor-Joining (NJ) und Unweighted Pair Group Method with Arithmetic Mean (UPGMA) als Distanz basierte Methoden zur Generierung eines phylogentischen Baumes zur Verfügung. Bei der Berechnung der entsprechenden Distanzen gibt es einfache Methoden und komplexere Modelle die man verwenden kann.

Zu den einfachen Methoden gehört, die Distanz als die Anzahl der Unterschiede (No. of differences) oder als proportionale Distanz (p-distance) zu bewerten. Die komplexeren Modelle (z.B. das Kimura 2-parameter Model) beziehen bei der Berechnung der Distanz verschiedene andere Faktoren mit ein. Zum Beispiel werden die DNA Veränderungen Transition und Transversion separat Betrachtet und können mit unterschiedlichen Raten belegt werden, was im Umkehrschluss dazu führt, dass die Distanz abhängig davon ist, wie viele Transversionen und Trasitionen im Vergleich vorkommen.
Neben der Wahl des Models bzw. der Methode, spielt auch der Umgang mit fehlenden Daten und Deletionen (gaps) eine nicht unwesentliche Rolle bei der Erstellung der Distanz Matrix. Zur Wahl steht, alle Positionen mit fehlenden Daten und Deletionen vor dem paarweise Vergleich aus dem Datensatz zu löschen (complete deletion), die entsprechenden Positionen nur beim Vergleich zweier Sequenzen zu ignorieren (pairwise deletion) oder Positionen deren Datengehalt unter einem bestimmten Maß liegt vor dem paarweise Vergleich aus dem Datensatz zu löschen (partial deletion).
Für unsere Zwecke - die DNA Diagnostik - spielen in erster Linie nur die tatsächlichen Unterschiede eine Rolle, weswegen wir auf die Verwendung eines komplexen Modells verzichten und die p-distance als Methode zur Erstellung der Distanz Matrix wählen. Beim Umgang mit fehlenden Date und Deletionen richten wir uns in der Regel nach der Art des Markers. Bei kodierenden Regionen (z.B. rbcL, matK), wo wir in der Regel keine Deletionen finden, verwenden wir complete deletion und bei Spacern (z.B. psbA-trnH) und Introns nutzen wir pairwise deletion um die entsprechende Information (InDels) nicht zu verlieren.

Voraussetzung für ein gutes Ergebnis ist allerdings ein gutes Alignment!

Wenn wir pairwise deletion verwenden, sollten das Alignment möglichst komplett sein (keine Positionen mit fehlenden Daten enthalten).

Simulationsdaten:

Ein kompletter Datensatz...

S1: ACTACTACTACTACT
S2: ACTACTACTACTATT
S3: ACTACTACTACTATT
S4: ACTACTACTACTAAT
S5: ACTACTACTACTAAT

NJ Baum auf S1 gewurzelt, pairwise deletion, kompletter Datensatz
Der selbe Datensatz mit fehlenden Daten...

S1: ACTACTACTACTACT
S2: ACTACTACTACTATT
S3: ACTACTACTACTATT
S4: --------TACTAAT
S5: ACTACTACTACTAAT

NJ Baum auf S1 gewurzelt, pairwise deletion, Datensatz mit fehlenden Informationen
Man sieht hier wie sich die Distanzen im Baum durch unvollständige Daten verändern und die "wahre" Situation verzerrt wird.

Bei complete deletion ist es besonders wichtig, dass keine sehr kurzen oder sehr schlechten Sequenzen (viele N's) enthalten sind, da sich das Löschen an der schlechstesten und kürzesten Sequenz orientiert und die Gefahr relevante Daten zu verlieren hoch ist.

So weit zu den wichtigsten Punkten...

Thomas