Riassunto Tecnico: Entropia dell'Eredità negli Alberi di Lineaggio Cellulare

1. Definizione del Problema

Quantificare l'impatto della trasmissione ereditaria all'interno degli alberi di lineaggio cellulare rimane una sfida fondamentale in biologia. Sebbene il profiling a singola cellula permetta la caratterizzazione degli stati cellulari, combinare questo approccio con il classico tracciamento del lineaggio per monitorare le ramificazioni ereditarie è difficile. Una sfida specifica sorge con l'inattività (cellule che smettono di dividersi, inclusi senescenza, differenziamento o quiescenza). A differenza della proliferazione attiva, dove le mutazioni ereditarie amplificano il fenotipo nella progenie, l'inattività è auto-limitante: una volta che una cellula diventa inattiva, il suo ramo termina, cancellando le informazioni a valle. Ciò rende difficile determinare se l'emergere dell'inattività sia guidato da eventi stocastici casuali o da fattori ereditari (epigenetici) che regolano l'uscita dal ciclo cellulare. Inoltre, distinguere tra mutazioni genetiche e eredità epigenetica è cruciale, particolarmente nelle colonie di cellule staminali dove l'eterogeneità inter-colonia influenza il potenziale terapeutico.

2. Metodologia

Gli autori introducono un metodo indipendente dal modello basato sull'entropia di Shannon applicata alla topologia degli alberi di lineaggio, concentrandosi specificamente sulla distribuzione delle cellule inattive.

A. Acquisizione dei Dati

Lo studio utilizza il tracciamento del lineaggio a singola cellula ad alta definizione di 32 colonie clonali di cellule stromali del midollo osseo umano (BMSC). Gli esperimenti sono stati condotti tramite microscopia a contrasto di fase in time-lapse fino al raggiungimento della settima generazione ($k_{max}=7$) da parte di tutte le cellule replicanti. Le cellule sono state classificate come:

• Attive: Mitosi osservata.
• Inattive (G0/Morte): Nessuna divisione per 84 ore o morte cellulare.
• Foglie (Leaves): Cellule a $k=7$ tracciate solo fino alla nascita.

B. Sbilanciamento dell'Inattività

Per quantificare l'impatto topologico dell'inattività, gli autori definiscono lo sbilanciamento dell'inattività ($I_m$) per ogni mitosi $m$ nell'albero:
$$I_m = |M_{m}^{sinistra} - M_{m}^{destra}|$$
dove $M$ rappresenta il numero di foglie mancanti a $k_{max}$ nei rami sinistro e destro. Un grande sbilanciamento in un nodo specifico suggerisce una divergenza nella probabilità di inattività tra i due rami figli, indicando potenzialmente un evento mutazionale a monte.

C. Entropia dell'Eredità

La metrica centrale è l'entropia dell'eredità ($S$), calcolata dai pesi normalizzati ($w_m$) degli sbilanciamenti attraverso tutte le mitosi del lineaggio:
$$w_m = \frac{I_m}{\sum_n I_n}, \quad S = -\sum_m w_m \log w_m$$

• Alta Entropia: Gli sbilanciamenti sono dispersi omogeneamente attraverso l'albero, suggerendo fluttuazioni casuali e non ereditarie.
• Bassa Entropia: Gli sbilanciamenti sono concentrati su pochi nodi, indicando un evento ereditario specifico (mutazione) che ha alterato la probabilità di inattività per un determinato sotto-ramo.

D. Insiemi Nulli e Validazione

Per determinare la significatività statistica, l'entropia dell'albero biologico viene confrontata con tre ensemble nulli non ereditari:

1. Galton-Watson (GW): Probabilità costante di inattività ($q$) per tutte le cellule.
2. GW Non Omogeneo (NHGW): Probabilità di inattività dipendente dalla generazione ($q_k$).
3. Rimescolamento Casuale (RR): Il metodo preferito. Questo consiste nello scambiare le posizioni delle cellule all'interno di ogni generazione preservando il numero totale di cellule inattive per generazione. Ciò recide tutti i legami ereditari madre-figlia pur mantenendo i vincoli strutturali dell'albero.

Il P-value è definito come la frazione di alberi dell'ensemble nullo con entropia $S \leq S_{biologica}$. Un basso P-value indica una forte evidenza di eredità.

E. Fitting del Modello (SKH)

Gli autori propongono un modello Stocastico Kirkwood-Holliday (SKH) per descrivere i dati. In questo modello:

• Le cellule "wild" mutano in uno stato "mutato" con probabilità $p$.
• Le cellule mutate (e la loro progenie) hanno una probabilità $q$ di diventare inattive.
• Crucialmente, il modello incorpora un ritardo di mutazione ($l$), ovvero il numero di generazioni tra l'evento di mutazione e l'espressione fenotipica dell'inattività.

3. Contributi Chiave

1. Entropia dell'Eredità: Una nuova metrica, agnostica rispetto al modello, per rilevare la struttura ereditaria negli alberi di lineaggio basata esclusivamente sulla topologia dei rami troncati (inattivi).
2. Validazione degli Modelli Nulli: Uno studio di benchmarking rigoroso che dimostra come, mentre i modelli GW e NHGW falliscono nel distinguere alberi ereditari da alberi non ereditari non omogenei (producendo falsi positivi), il metodo del Rimescolamento Casuale (RR) è robusto, sensibile e indipendente dal modello.
3. Quantificazione del Ritardo di Mutazione: La capacità di localizzare la specifica mitosi in cui la probabilità di inattività è cambiata e di misurare il ritardo di mutazione (il ritardo tra il cambiamento epigenetico e l'espressione fenotipica).
4. Modello Stochastic Kirkwood-Holliday: Un modello stocastico generalizzato che cattura con successo le caratteristiche ereditarie dei lineaggi BMSC, incluso il ritardo tra mutazione e inattività.

4. Risultati

• Struttura Ereditaria: Nel dataset sperimentale di 32 colonie BMSC, 21 colonie (75%) hanno mostrato un carattere ereditario significativo ($P < 0.05$) utilizzando il test RR. Le restanti colonie avevano troppe cellule inattive per produrre un risultato statisticamente significativo, il che non implica necessariamente una mancanza di eredità.
• Analisi dell'Entropia: Gli alberi biologici hanno costantemente mostrato un'entropia significativamente inferiore rispetto ai loro controparti rimescolati casualmente, provando che la distribuzione delle cellule inattive non è casuale ma segue un pattern ereditario.
• Ritardo di Mutazione: Il ritardo di mutazione medio è stato trovato essere di circa 2 generazioni. Ciò conferma che il cambiamento epigenetico che porta all'inattività avviene a monte dell'effettiva espressione fenotipica.
• Fit del Modello: Il modello SKH, adattato ai dati, ha fornito parametri $p \approx 0.1$ (probabilità di mutazione) e $q \approx 0.46$ (probabilità di inattività per le cellule mutate). Gli alberi sintetici generati da SKH con questi parametri somigliavano strettamente alla struttura topologica dei lineaggi BMSC sperimentali.
• Spaziale vs. Topologico: L'analisi ha confermato che il clustering delle cellule inattive non è un artefatto della prossimità fisica (correlazione spaziale) sul piatto di coltura, ma è una caratteristica topologica del lineaggio.

5. Significato e Rivendicazioni

L'articolo sostiene che l'eterogeneità osservata nelle colonie BMSC umane è guidata da fattori ereditari epigenetici che regolano l'uscita dal ciclo cellulare, piuttosto che da meccanismi puramente stocastici o genetici.

• Robustezza: Il test basato sull'entropia fornisce un modo semplice, robusto e agnostico rispetto al modello per valutare l'origine ereditaria della crescita differenziale in qualsiasi albero di lineaggio.
• Rilevanza Clinica: I risultati suggeriscono che l'eterogeneità inter-colonia osservata nei trapianti di BMSC (un ostacolo maggiore per le terapie di rigenerazione ossea) ha un'origine ereditaria epigenetica. Comprendere e controllare questi meccanismi ereditari potrebbe migliorare la consistenza delle terapie con cellule staminali.
• Avanzamento Metodologico: Stabilendo un "ritardo di mutazione", il metodo consente ai ricercatori di tornare indietro all'origine dei cambiamenti fenotipici anche quando il fenotipo (inattività) è auto-cancellante.
• Direzioni Future: Gli autori notano che, sebbene l'attuale metodo si concentri sull'inattività, la correlazione tra la topologia del lineaggio e i tempi di divisione suggerisce che anche la durata del ciclo cellulare possa essere regolata da meccanismi ereditari epigenetici, sebbene ciò richieda ulteriori studi.

Gli autori rimangono modesti, riconoscendo che il modello SKH è una semplificazione e che modelli più complessi potrebbero catturare meglio la varianza dei ritardi di mutazione. Essi sottolineano che il metodo rileva i cambiamenti di probabilità; se un lineaggio non presenta variazioni nell'espressione dell'inattività, il test non rileverà l'eredità, anche se esiste un meccanismo ereditario.