Tsallis-Entropie verfeinert den Fixationsindex für Populationsgenetik


Forscher haben mit $F_q$ eine informationstheoretische Verallgemeinerung des klassischen Fixationsindex $F_{\text{ST}}$ vorgestellt. Das neue Maß quantifiziert, welcher Anteil der Tsallis-$q$-Entropie innerhalb von Subpopulationen im Vergleich zur gepoolten Gesamtpopulation verloren geht. Bemerkenswert: Für $q=2$ reduziert sich $F_q$ exakt auf den varianzbasierten $F_{\text{ST}}$; für $q=1$ ergibt sich eine Shannon-Entropie-Variante, deren absolute Form der Transinformation zwischen Allelen und Populationszuordnungen entspricht.

Durch Variation von $q$ wirkt $F_q$ als spektraler Differenzierer: Kleine $q$-Werte gewichten seltene genetische Varianten stärker, während $q>1$ zunehmend häufige Varianten betont. Damit liefert $F_q$ eine feinere Auflösung der populationsgenetischen Differenzierung als $F_{\text{ST}}$, insbesondere wenn Allelfrequenzspektren verzerrt sind.

An realen Daten — 865 ozeanische Genome mit 1.823.000 untersuchten Stellen — sowie an kontrollierten genealogischen Simulationen (basierend auf 1.432 Gründern der HGDP- und 1000-Genomes-Panels mit 322.216 Stellen) zeigten die Autoren, dass $F_q$ in One-vs-Rest- und Leave-One-Out-Modi klar zuordnet, welche Subpopulationen eine regionale Populationsstruktur prägen. Zudem lassen sich Isolations- und Migrationsereignisse sowie Gründereffekte zeitlich empfindlich eingrenzen. $F_q$ eignet sich damit als ergänzendes Analyseinstrument höherer Auflösung für Simulationsaudits und die Charakterisierung von Populationsstrukturen — und demonstriert zugleich, wie physikalisch geprägte Entropiekonzepte neue Anwendungen in der Genomik finden.

Hintergrund: Die Tsallis-Entropie wurde 1988 vom Physiker Constantino Tsallis als Verallgemeinerung der Boltzmann-Gibbs-Shannon-Entropie für nicht-extensive Systeme eingeführt und wird vor allem in der statistischen Physik und Komplexitätsforschung eingesetzt. Ihre Übertragung auf die Populationsgenetik verbindet informationstheoretische Maße mit klassischen Evolutionsmodellen.

Lernimpuls

Kernidee: Populationsgenetische Differenzierung misst man klassisch mit dem Fixationsindex $F_{\text{ST}}$ — dem Anteil genetischer Vielfalt, der verloren geht, wenn man von der gepoolten Gesamtpopulation auf die Subpopulationen herabblickt. Die neue $F_q$-Statistik tauscht dabei nur die Basisgröße aus: Statt der Shannon-Entropie $H = -\sum_i p_i \ln p_i$ verwendet sie die Tsallis-Entropie $S_q = \frac{1}{q-1}\left(1 - \sum_i p_i^q\right)$, in der der Entropie-Index $q$ die Allelfrequenzen nichtlinear umgewichtet. $F_q$ ist dann wieder eine relative Entropiedifferenz — Tsallis-Entropie des Pools minus der gewichteten mittleren Entropie der Subpopulationen, normiert auf die Gesamtentropie. Für $q = 2$ gilt die bemerkenswerte Identität: Weil $S_2 = 1 - \sum_i p_i^2$ über die Homozygotie mit der Allelfrequenzvarianz verknüpft ist (bei einem biallelen Marker ist $1 - \sum_i p_i^2 = 2\,p(1-p)$), fällt $F_2$ exakt mit dem varianzbasierten $F_{\text{ST}}$ zusammen. Im Grenzfall $q \to 1$ geht die Tsallis- in die Shannon-Entropie über, und der absolute (nicht normierte) Wert der Differenz entspricht der Transinformation zwischen Allelidentität und Populationszugehörigkeit — also der Information, die die Kenntnis der Population über das Allel liefert. Der praktische Gewinn liegt in den Zwischenwerten: Kleine $q$-Werte ($q < 1$) heben seltene Varianten überproportional hervor, große ($q > 1$) betonen die frequenten. Variation von $q$ wirkt damit wie ein spektraler Differenzierer, der das Allelfrequenzspektrum unter verschiedenen „Vergrößerungen" abtastet und verrät, ob eine beobachtete Populationsstruktur von häufigen Kernvarianten oder von seltenen Randvarianten getragen wird — genau dann wertvoll, wenn Gründereffekte oder Migration das Spektrum verzerrt haben.

Häufige Fehlvorstellung: „$F_q$ ist nur ein neuer Name für $F_{\text{ST}}$." Falsch — nur für $q = 2$ deckt sich $F_q$ exakt mit dem varianzbasierten Index; für jedes andere $q$ werden die Allelfrequenzen anders gewichtet, und Populationen können sogar in eine andere Rangfolge kommen. Zweite Fehlvorstellung: „Tsallis-Entropie ist einfach Shannon-Entropie mit einem anderen Logarithmus." Falsch — ein Logarithmenwechsel ändert die Entropie nur um einen konstanten Faktor und bliebe derselbe Funktional; die Tsallis-Form $S_q = \frac{1}{q-1}(1 - \sum_i p_i^q)$ ist dagegen eine echte nichtlineare Deformation, die Shannon-Entropie nur im Grenzfall $q \to 1$ wiederherstellt. Dritte Fehlvorstellung: „$q$ ist eine messbare Eigenschaft der Daten, etwa eine Temperatur." Falsch — $q$ ist ein vom Forscher gewählter Analyseparameter (ein Deformations- bzw. Fokusparameter), kein geschätzter Populationsparameter; man variiert ihn bewusst, um Sensitivitäten zu vergleichen. Vierte Fehlvorstellung: „Ein höheres $F_q$ bedeutet automatisch stärkere Differenzierung." Nicht ungeprüft — die Aussage gilt nur bei festgehaltenem $q$ und identischer Datengrundlage; über verschiedene $q$-Werte hinweg kann derselbe Datensatz bei kleinem $q$ (strukturgebende seltene Varianten) hohe und bei großem $q$ moderate Werte liefern, und gerade dieses Muster ist die eigentliche diagnostische Information.

Prüfe dich:

  1. Zwei Subpopulationen gleicher Größe haben an einem biallelen Marker die Allelfrequenzen $p_1 = 0{,}8$ und $p_2 = 0{,}2$. Berechnen Sie die Tsallis-Entropie $S_2$ für jede Subpopulation und für den gepoolten Datensatz ($\bar{p} = 0{,}5$) und zeigen Sie, dass $F_2 = 0{,}36$ dem varianzbasierten $F_{\text{ST}}$ mit $H_T = 2\bar{p}(1-\bar{p})$ und $H_S$ als mittlerer Heterozygotie entspricht.
  2. Zeigen Sie mit der Regel von L’Hôpital (oder per Ableitung nach $q$ an der Stelle $q = 1$), dass $\lim_{q \to 1} \frac{1 - \sum_i p_i^q}{q-1} = -\sum_i p_i \ln p_i$ gilt — die Tsallis-Entropie also kontinuierlich in die Shannon-Entropie übergeht.
  3. Erklären Sie, warum $q < 1$ seltene Varianten relativ aufwertet: Vergleichen Sie dazu den Beitrag eines raren Allels mit $p = 0{,}001$ unter $p^{0{,}5}$ und $p^{2}$ und diskutieren Sie, was dies für die Detektion von Gründer-Effekten bedeutet.
  4. Begründen Sie, warum der absolute Wert der Entropiedifferenz für $q \to 1$ der Transinformation $I(\text{Allel};\text{Population}) = H(\text{Allel}) - H(\text{Allel}\mid\text{Population})$ entspricht. Warum ist diese Größe ein natürliches, nicht normiertes Differenzierungsmaß?
  5. Ein Datensatz stammt aus einem Gründungsereignis mit wenigen Gründern. Skizzieren Sie, wie dies das Allelfrequenzspektrum verzerrt (viele seltene Varianten gehen verloren, einige werden per Zufall hochfrequent), und erklären Sie, warum ein einziger fester $q$-Wert solche Strukturen übersehen kann, während das $q$-Spektrum von $F_q$ sie sichtbar macht.
  6. In der Studie werden One-vs-Rest- und Leave-One-Out-Modi verwendet. Erläutern Sie, wie man mit $F_q$ in diesen Modi identifizieren kann, welche Subpopulationen die regionale Populationsstruktur prägen — und warum der Entzug einer strukturbildenden Population das $F_q$-Profil besonders stark verändert.

Bei Schwierigkeiten: Stellen Sie sich zwei Dosen mit bunten Glaskugeln vor — die Farben sind die Allele, die Dosen sind die Subpopulationen. Die Entropie misst die „Farbvielfalt": Mischt man beide Dosen in einem großen Behälter, ist die Vielfalt dort meist größer als in jeder Einzeldose, und genau dieser Zuwachs (relativ gemessen) ist der Fixationsindex. Die Tsallis-Variante ändert nur die Brille, mit der Sie Vielfalt bewerten: Mit $q < 1$ tragen Sie eine Weitwinkel-„Sensationsbrille", die auch der einzelnen Kugel einer seltenen Farbe fast so viel Gewicht gibt wie einer Massenfarbe; mit $q > 1$ tragen Sie ein Teleobjektiv, das nur die dominanten Farben scharfstellt und Raritäten praktisch ausblendet. Für $q = 2$ wird die Brille zur Varianzbrille — Vielfalt wird zu Quadratsummen, und das ist genau die alte $F_{\text{ST}}$-Rechnung. Warum die Rechnung für $q \to 1$ in Shannon übergeht: Wenn Sie den Exponenten $q$ sanft auf 1 zubewegen, geht jeder Term $p_i^q$ in $p_i$ über, und die Korrektur $\frac{1}{q-1}$ verwandelt die kleine Krümmung dieses Übergangs exakt in den Logarithmus von $-\sum_i p_i \ln p_i$ — die Tsallis-Entropie ist gewissermaßen eine gespannte Feder, die in Ruheposition die Shannon-Entropie ist. Und die Transinformation? Denken Sie an ein Ratespiel: Wie viel weniger unsicher sind Sie über die Kugelfarbe, wenn ich Ihnen verrate, aus welcher Dose sie stammt? Genau diese Unsicherheitsreduktion ist die Differenzierung zwischen den Populationen.

Entitäten: tsallis-entropie | konzept fixationsindex | konzept f_q-statistik | konzept shannon-entropie | konzept transinformation | konzept gründereffekt | konzept allelfrequenzspektrum | konzept one-vs-rest-analyse | methode hgdp-panel | person 1000-genomes-projekt | person


📄 Quelle

Nachrichten-Artikel

📄 Original-Publikation 📄 Original-Publikation 📄 Original-Publikation


🧪 Verwandte Rechner

Mit diesen interaktiven Werkzeugen können Sie das Thema vertiefen:

🔬 ph-rechner → 🔬 saeuren-basen-gleichgewicht →