2025-11-15T01:28:11.271605

Decoding Positive Selection in Mycobacterium tuberculosis with Phylogeny-Guided Graph Attention Models

Wang, Campino, Clark et al.
Positive selection drives the emergence of adaptive mutations in Mycobacterium tuberculosis, shaping drug resistance, transmissibility, and virulence. Phylogenetic trees capture evolutionary relationships among isolates and provide a natural framework for detecting such adaptive signals. We present a phylogeny-guided graph attention network (GAT) approach, introducing a method for converting SNP-annotated phylogenetic trees into graph structures suitable for neural network analysis. Using 500 M. tuberculosis isolates from four major lineages and 249 single-nucleotide variants (84 resistance-associated and 165 neutral) across 61 drug-resistance genes, we constructed graphs where nodes represented isolates and edges reflected phylogenetic distances. Edges between isolates separated by more than seven internal nodes were pruned to emphasise local evolutionary structure. Node features encoded SNP presence or absence, and the GAT architecture included two attention layers, a residual connection, global attention pooling, and a multilayer perceptron classifier. The model achieved an accuracy of 0.88 on a held-out test set and, when applied to 146 WHO-classified "uncertain" variants, identified 41 candidates with convergent emergence across multiple lineages, consistent with adaptive evolution. This work demonstrates the feasibility of transforming phylogenies into GNN-compatible structures and highlights attention-based models as effective tools for detecting positive selection, aiding genomic surveillance and variant prioritisation.
academic

माइकोबैक्टीरियम ट्यूबरकुलोसिस में सकारात्मक चयन को डिकोड करना: फाइलोजेनी-गाइडेड ग्राफ अटेंशन मॉडल

बुनियादी जानकारी

  • पेपर ID: 2510.08703
  • शीर्षक: माइकोबैक्टीरियम ट्यूबरकुलोसिस में सकारात्मक चयन को डिकोड करना फाइलोजेनी-गाइडेड ग्राफ अटेंशन मॉडल के साथ
  • लेखक: लिनफेंग वांग, सुसाना कैम्पिनो, टेन जी. क्लार्क, जोडी ई. फेलान
  • वर्गीकरण: q-bio.PE (जनसंख्या और विकास), cs.LG (मशीन लर्निंग)
  • संस्थान: लंदन स्कूल ऑफ हाइजीन एंड ट्रॉपिकल मेडिसिन
  • पेपर लिंक: https://arxiv.org/abs/2510.08703

सारांश

यह अध्ययन ट्यूबरकुलोसिस बैक्टीरिया में सकारात्मक चयन संकेतों का पता लगाने के लिए फाइलोजेनेटिक ट्री-गाइडेड ग्राफ अटेंशन नेटवर्क (GAT) विधि प्रस्तावित करता है। SNP एनोटेटेड फाइलोजेनेटिक ट्री को तंत्रिका नेटवर्क विश्लेषण के लिए उपयुक्त ग्राफ संरचना में परिवर्तित करके, यह विधि 500 ट्यूबरकुलोसिस बैक्टीरिया आइसोलेट्स और 249 एकल न्यूक्लिओटाइड वेरिएंट्स पर 0.88 की सटीकता प्राप्त करती है, और 41 अनुकूलनीय विकास विशेषताओं वाले उम्मीदवार वेरिएंट्स की सफलतापूर्वक पहचान करती है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

क्षय रोग (TB) अभी भी वैश्विक संक्रामक रोग मृत्यु के प्रमुख कारणों में से एक है, 2024 में 109 लाख लोगों की मृत्यु हुई। दवा प्रतिरोध के विकास ने इस महामारी को बढ़ा दिया है, जिसमें 40 लाख नए TB मामले कम से कम एक प्रथम-पंक्ति दवा राइफैम्पिसिन के प्रति प्रतिरोधी हैं। सकारात्मक चयन ट्यूबरकुलोसिस बैक्टीरिया के विकास का एक मुख्य चालक है, जो अनुकूलनीय उत्परिवर्तन के उद्भव को बढ़ावा देता है, दवा प्रतिरोध, संचरणीयता और विषाणुता को प्रभावित करता है।

अनुसंधान का महत्व

  1. नैदानिक महत्व: सकारात्मक चयन उत्परिवर्तन की पहचान दवा प्रतिरोध तंत्र को समझने और उपचार रणनीति के निर्देशन के लिए महत्वपूर्ण है
  2. विकास जीव विज्ञान मूल्य: ट्यूबरकुलोसिस बैक्टीरिया की कठोर क्लोनल जनसंख्या संरचना और पुनर्संयोजन की कमी इसे अनुकूलनीय विकास के अध्ययन के लिए एक आदर्श मॉडल बनाती है
  3. सार्वजनिक स्वास्थ्य आवश्यकता: जीनोमिक निगरानी को अनुकूलनीय लाभ वाले वेरिएंट्स की तेजी से और सटीक पहचान की आवश्यकता है

मौजूदा विधियों की सीमाएं

  1. पारंपरिक फाइलोजेनेटिक विश्लेषण: मानव व्याख्या पर निर्भर, बड़े पैमाने पर डेटा को संभालना मुश्किल
  2. मानक GNN विधियां: फाइलोजेनेटिक जानकारी और उत्परिवर्तन पैटर्न को प्रभावी ढंग से एकीकृत नहीं कर सकते
  3. मौजूदा वर्गीकरण विधियां: विकास पृष्ठभूमि पर विचार की कमी, महत्वपूर्ण अनुकूलनीय संकेतों को छोड़ सकते हैं

मुख्य योगदान

  1. विधि नवाचार: फाइलोजेनेटिक ट्री को ग्राफ न्यूरल नेटवर्क संगत संरचना में परिवर्तित करने की विधि पहली बार प्रस्तावित करना
  2. आर्किटेक्चर डिजाइन: किनारे की लंबाई की जानकारी को एकीकृत करने वाली ग्राफ अटेंशन नेटवर्क आर्किटेक्चर विकसित करना, जो टोपोलॉजी संरचना और उत्परिवर्तन पैटर्न दोनों को संभाल सकता है
  3. व्यावहारिक अनुप्रयोग: WHO "अनिश्चित" वेरिएंट वर्गीकरण में 41 अभिसारी उपस्थिति पैटर्न वाले उम्मीदवार अनुकूलनीय वेरिएंट्स की पहचान करना
  4. उपकरण विकास: संपूर्ण ओपन-सोर्स कोड और डेटा प्रोसेसिंग पाइपलाइन प्रदान करना

विधि विवरण

कार्य परिभाषा

इनपुट: SNP एनोटेटेड फाइलोजेनेटिक ट्री, जहां नोड्स ट्यूबरकुलोसिस बैक्टीरिया आइसोलेट्स का प्रतिनिधित्व करते हैं, किनारे फाइलोजेनेटिक दूरी को दर्शाते हैं आउटपुट: द्विआधारी वर्गीकरण पूर्वानुमान, यह निर्धारित करना कि क्या विशिष्ट SNP सकारात्मक चयन के अधीन है बाधाएं: फाइलोजेनेटिक संबंधों की पूर्णता बनाए रखते हुए, ग्राफ न्यूरल नेटवर्क इनपुट आवश्यकताओं के अनुकूल

मॉडल आर्किटेक्चर

डेटा संरचना रूपांतरण

  1. ग्राफ निर्माण: फाइलोजेनेटिक ट्री को अप्रत्यक्ष ग्राफ में परिवर्तित करना, नोड्स आइसोलेट्स का प्रतिनिधित्व करते हैं, किनारे वजन आंतरिक नोड गणना दूरी है
  2. किनारा छंटाई: 7 से अधिक आंतरिक नोड्स से अलग किए गए नमूनों के बीच किनारों को हटाना, स्थानीय विकास संरचना को उजागर करना
  3. नोड विशेषताएं: SNP उपस्थिति/अनुपस्थिति स्थिति को एन्कोड करने के लिए बाइनरी संकेतक का उपयोग करना

GAT आर्किटेक्चर डिजाइन

चरण 1: दोहरी-परत ग्राफ अटेंशन नेटवर्क
- पहली परत: 8 अटेंशन हेड्स, प्रत्येक हेड 32 आउटपुट विशेषताएं
- दूसरी परत: एकल अटेंशन हेड, 256-आयामी आउटपुट
- अवशिष्ट कनेक्शन: दोनों परतों के आउटपुट को जोड़ना

चरण 2: वैश्विक पूलिंग और वर्गीकरण
- वैश्विक अटेंशन पूलिंग
- बहु-परत परसेप्ट्रॉन वर्गीकारक (256→32→2)

अटेंशन तंत्र

मुख्य नवाचार किनारा-जागरूक अटेंशन गणना में है:

hi(l+1)=σ(∑j∈N(i)αijWhj(l))h_i^{(l+1)} = \sigma\left(\sum_{j \in N(i)} \alpha_{ij} W h_j^{(l)}\right)

जहां अटेंशन वजन αij\alpha_{ij} नोड विशेषताओं और किनारे की लंबाई की जानकारी दोनों पर विचार करता है: αij=softmax(σ(aT[Whi∥Whj]+b⋅edgeij))\alpha_{ij} = \text{softmax}\left(\sigma\left(\mathbf{a}^T [Wh_i \| Wh_j] + b \cdot edge_{ij}\right)\right)

तकनीकी नवाचार बिंदु

  1. फाइलोजेनेटिक-जागरूक: पहली बार आंतरिक नोड गणना को किनारे वजन के रूप में ग्राफ न्यूरल नेटवर्क में पेश करना
  2. अनुकूली छंटाई: दूरी थ्रेशहोल्ड के माध्यम से स्थानीय पड़ोस संरचना को बनाए रखना, शोर को कम करना
  3. बहु-पैमाने अटेंशन: नोड-स्तर और किनारा-स्तर जानकारी के अटेंशन तंत्र को जोड़ना
  4. अवशिष्ट डिजाइन: गहरे नेटवर्क के प्रशिक्षण स्थिरता को सुनिश्चित करना

प्रायोगिक सेटअप

डेटासेट

  • नमूना आकार: 500 ट्यूबरकुलोसिस बैक्टीरिया नैदानिक नमूने
  • वंशावली कवरेज: चार मुख्य वंशावली (L1-L4), वितरण L1:8, L2:175, L3:109, L4:223
  • वेरिएंट डेटा: 249 SNP वेरिएंट्स, 61 दवा प्रतिरोध जीन में फैले हुए
  • लेबल वितरण: 84 WHO-पुष्टि दवा प्रतिरोध-संबंधित उत्परिवर्तन, 165 तटस्थ वेरिएंट्स

डेटा प्रोसेसिंग प्रवाह

  1. अनुक्रम प्रोसेसिंग: गुणवत्ता नियंत्रण और संरेखण के लिए Trimmomatic और BWA-mem का उपयोग
  2. वेरिएंट कॉलिंग: BCF/VCF उपकरण सेट, >10x कवरेज
  3. फाइलोजेनेटिक पुनर्निर्माण: अधिकतम संभावना ट्री के लिए RAxML
  4. डेटा विभाजन: प्रशिक्षण सेट 149, सत्यापन सेट 50, परीक्षण सेट 50

मूल्यांकन मेट्रिक्स

  • सटीकता (Accuracy): 0.88
  • AUC: 0.89
  • F1 स्कोर: 0.81
  • संवेदनशीलता: 0.76
  • विशिष्टता: 0.94

तुलनात्मक विश्लेषण

हालांकि पेपर पारंपरिक विधियों की प्रत्यक्ष तुलना प्रदान नहीं करता है, लेकिन WHO वर्गीकरण के साथ सामंजस्य सत्यापन के माध्यम से विधि की प्रभावशीलता को सत्यापित किया गया है।

प्रायोगिक परिणाम

मुख्य परिणाम

50 परीक्षण नमूनों के holdout डेटासेट पर:

  • समग्र प्रदर्शन: 0.88 सटीकता, अच्छी सामान्यीकरण क्षमता दिखाता है
  • वर्ग संतुलन: उच्च विशिष्टता (0.94) और मध्यम संवेदनशीलता (0.76), स्क्रीनिंग अनुप्रयोगों के लिए उपयुक्त
  • जैविक तर्कसंगतता: मॉडल लगभग पूरी तरह से समानार्थी उत्परिवर्तन को बाहर करता है, कार्यात्मक अपेक्षाओं के अनुरूप

अटेंशन विश्लेषण

शीर्ष-k अटेंशन गुणवत्ता (TAM) विश्लेषण के माध्यम से खोज:

  • अटेंशन एकाग्रता: शीर्ष 10% किनारे कुल अटेंशन का 44.1% कैप्चर करते हैं
  • जैविक अर्थ: उच्च अटेंशन किनारे मुख्य रूप से उत्परिवर्तन विविधता में समृद्ध केंद्रीय नोड्स को जोड़ते हैं
  • संरचना समझ: मॉडल विकास की दृष्टि से महत्वपूर्ण ग्राफ क्षेत्रों की पहचान और ध्यान केंद्रित कर सकता है

व्यावहारिक अनुप्रयोग सत्यापन

146 WHO "अनिश्चित" वेरिएंट्स में:

  • पूर्वानुमान परिणाम: 27 (18.5%) को सकारात्मक चयन के रूप में पूर्वानुमानित किया गया
  • अभिसारी पैटर्न: 41 उम्मीदवार वेरिएंट्स कई वंशावली में अभिसारी रूप से दिखाई देते हैं
  • कार्यात्मक संबंध: ज्ञात दवा प्रतिरोध उत्परिवर्तन और मुआवजे उत्परिवर्तन की पहचान

महत्वपूर्ण खोजें

  1. embA c.-43G>C: 43 उप-वंशावली में दिखाई देता है, MDR+ आवृत्ति 47.48%
  2. rpoC श्रृंखला उत्परिवर्तन: कई मुआवजे उत्परिवर्तन सफलतापूर्वक पहचाने गए
  3. ubiA वेरिएंट्स: एथमब्यूटोल प्रतिरोध से संबंधित नए उम्मीदवार उत्परिवर्तन

संबंधित कार्य

पारंपरिक फाइलोजेनेटिक विधियां

  • dN/dS अनुपात विश्लेषण: चयन दबाव का पता लगाने की शास्त्रीय विधि
  • फाइलोजेनेटिक अभिसरण विश्लेषण: स्वतंत्र उत्पत्ति घटनाओं की मानव पहचान
  • आणविक घड़ी विश्लेषण: उत्परिवर्तन घटना समय का अनुमान

ग्राफ न्यूरल नेटवर्क अनुप्रयोग

  • जैविक नेटवर्क विश्लेषण: प्रोटीन इंटरैक्शन नेटवर्क में GNN का अनुप्रयोग
  • फाइलोजेनेटिक अनुमान: गहन शिक्षा-आधारित ट्री पुनर्निर्माण विधियां
  • जीनोमिक विश्लेषण: अनुक्रम वर्गीकरण और कार्यात्मक पूर्वानुमान

इस पेपर के लाभ

  1. पहली बार: फाइलोजेनेटिक ट्री को GNN इनपुट में व्यवस्थित रूप से परिवर्तित करने वाला पहला
  2. एकीकृतता: टोपोलॉजी और विशेषता जानकारी दोनों पर विचार करना
  3. व्यावहारिकता: दवा प्रतिरोध निगरानी की वास्तविक आवश्यकताओं पर सीधा अनुप्रयोग

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. तकनीकी व्यवहार्यता: फाइलोजेनेटिक ट्री से ग्राफ न्यूरल नेटवर्क में रूपांतरण की व्यवहार्यता सफलतापूर्वक सिद्ध करना
  2. पूर्वानुमान क्षमता: GAT मॉडल सकारात्मक चयन संकेतों को प्रभावी ढंग से पहचान सकता है
  3. अनुप्रयोग मूल्य: WHO अनिश्चित वेरिएंट वर्गीकरण में कई मूल्यवान उम्मीदवारों की खोज

सीमाएं

  1. नमूना आकार: अपेक्षाकृत छोटा डेटासेट (249 वेरिएंट्स) मॉडल सामान्यीकरण क्षमता को सीमित कर सकता है
  2. लेबल शोर: दवा प्रतिरोध को सकारात्मक चयन प्रॉक्सी के रूप में उपयोग करना वर्गीकरण त्रुटियां पेश कर सकता है
  3. विधि निर्भरता: इनपुट के रूप में उच्च गुणवत्ता की फाइलोजेनेटिक ट्री की आवश्यकता
  4. कम्प्यूटेशनल जटिलता: बड़े पैमाने पर डेटासेट प्रोसेसिंग दक्षता सत्यापन की आवश्यकता है

भविष्य की दिशाएं

  1. विस्तारित अनुप्रयोग: अन्य रोगजनकों के अनुकूलनीय विकास अनुसंधान के लिए लागू करना
  2. विधि सुधार: ग्राफ-अज्ञेयवादी शिक्षण आर्किटेक्चर विकसित करना
  3. बहु-मोडल एकीकरण: फेनोटाइप और जीनोटाइप डेटा को जोड़ना
  4. वास्तविक समय निगरानी: ऑनलाइन दवा प्रतिरोध निगरानी प्रणाली का निर्माण

गहन मूल्यांकन

लाभ

  1. मजबूत नवाचार: पहली बार फाइलोजेनेटिक जानकारी को गहन शिक्षा ढांचे में व्यवस्थित रूप से एकीकृत करना
  2. विधि तर्कसंगतता: किनारा छंटाई रणनीति और अटेंशन तंत्र डिजाइन जैविक अंतर्ज्ञान के अनुरूप है
  3. व्यावहारिक मूल्य: तपेदिक दवा प्रतिरोध निगरानी की वास्तविक आवश्यकताओं को सीधे पूरा करना
  4. ओपन-सोर्स योगदान: संपूर्ण कोड और डेटा प्रदान करना, क्षेत्र विकास को बढ़ावा देना

कमियां

  1. तुलना अपर्याप्त: पारंपरिक फाइलोजेनेटिक विधियों के साथ मात्रात्मक तुलना की कमी
  2. सीमित सत्यापन: पूर्वानुमान परिणामों का प्रायोगिक सत्यापन अभी भी आगे के अनुसंधान की आवश्यकता है
  3. सामान्यीकरण अज्ञात: अन्य रोगजनकों पर प्रयोज्यता अभी तक सत्यापित नहीं है
  4. सैद्धांतिक आधार: GAT इस कार्य के लिए विशेष रूप से उपयुक्त क्यों है इसका सैद्धांतिक विश्लेषण की कमी

प्रभाव

  1. पद्धति विज्ञान योगदान: फाइलोजेनेटिक जीनोमिक्स के लिए नई विश्लेषण उपकरण प्रदान करना
  2. अनुप्रयोग संभावनाएं: संक्रामक रोग निगरानी और विकास जीव विज्ञान में व्यापक अनुप्रयोग संभावनाएं
  3. अंतःविषय मूल्य: विकास जीव विज्ञान, मशीन लर्निंग और सार्वजनिक स्वास्थ्य क्षेत्रों को जोड़ना

लागू परिदृश्य

  1. रोगज़नक निगरानी: नए उभरते दवा प्रतिरोध उत्परिवर्तन की वास्तविक समय पहचान
  2. विकास अनुसंधान: बड़े पैमाने पर अनुकूलनीय विकास संकेत पहचान
  3. दवा विकास: संभावित दवा प्रतिरोध लक्ष्यों का पूर्वानुमान
  4. महामारी विज्ञान: दवा प्रतिरोधी बैक्टीरिया तनाव के प्रसार पैटर्न को ट्रैक करना

संदर्भ

पेपर 26 महत्वपूर्ण संदर्भों का हवाला देता है, जो तपेदिक महामारी विज्ञान, फाइलोजेनेटिक विश्लेषण, ग्राफ न्यूरल नेटवर्क और अन्य कई क्षेत्रों को कवर करते हैं, अनुसंधान के लिए एक मजबूत सैद्धांतिक आधार प्रदान करते हैं।

समग्र मूल्यांकन: यह एक महत्वपूर्ण नवाचार मूल्य वाला अंतःविषय अनुसंधान पेपर है, जो गहन शिक्षा तकनीक को संक्रामक रोग विकास जीनोमिक्स में सफलतापूर्वक लागू करता है, तपेदिक दवा प्रतिरोध निगरानी के लिए नई तकनीकी विधियां प्रदान करता है। कुछ सीमाओं के बावजूद, इसके पद्धति विज्ञान योगदान और व्यावहारिक अनुप्रयोग मूल्य की पुष्टि की जानी चाहिए।