Carrier Burden and Birth Prevalence of Recessive Inborn Errors of Metabolism in South Asians: A gnomAD and GenomeIndia Study

ملخص تقني: عبء الحامل للمرض وانتشار الولادات لأخطاء التمثيل الغذائي الوراثية المتنحية لدى سكان جنوب آسيا

بيان المشكلة
تعد أخطاء التمثيل الغذراي الموروثة (IEMs)، رغم كونها نادرة بشكل فردي، مسببات هامة جماعياً للأمراض في مرحلة الطفولة المبكرة، والإعاقة الذهنية، والوفاة، وتنتقل غالباً بطريقة وراثية متنحية (AR). ويعد تقدير انتشار هذه الاضطرابات في جنوب آسيا أمراً صعباً بسبب عدم تجانس تصميمات الدراسات، ومحدودية تغطية الفحص، واستبعاد المتغيرات المتوقع أن تكون ممرضة. علاوة على ذلك، فإن التقديرات الحالية غالباً ما تفشل في مراعاة معدلات التزاوج الداخلي وزواج الأقارب المرتفعة في المنطقة، مما يزيد من حالة "التماثل الجيني الذاتي" (autozygosity) واحتمالية حالات المتماثل الجيني المتنحي. وقد اقترحت النماذج السابقة معدلات انتشار تتراوح بين 1 في 1,000 إلى 1 في 3,100 من المواليد الأحياء، ولكن قد تكون هذه الأرقام أقل من الواقع. تعالج هذه الدراسة الحاجة إلى تقدير شامل وقائم على الأسس الجينية لعبء الحامل للمرض وانتشار الولادات لأخطاء التمثيل الغذراي الوراثية المتنحية (AR-IEMs) عبر سكان جنوب آسيا (SAS)، مع التحقق من النتائج مقابل بيانات مجموعات مستقلة وتصحيح انحيازات التقصي.

المنهجية
استخدمت الدراسة مساراً حوسبياً متعدد الخطوات يدمج بيانات gnomAD v4 ومشروع GenomeIndia:

1. تنقية لوحة الجينات: تم اشتقاق لوحة تضم 696 جيناً من جينات أخطاء التمثيل الغذراي الوراثية المتنحية من Genomics England PanelApp (R98)، ورسم خرائطها وفق نظام GRCh38، ومواءمتها مع تسميات gnomAD.
2. مصادر البيانات:
   • المجموعة الأولية: بيانات gnomAD v4 لسكان جنوب آسيا (SAS) والتي تضم 86,254 فرداً.
   • مجموعة التحقق: الإحصاءات الملخصة لـ GenomeIndia لـ 9,768 فرداً غير مرتبطين عبر 83 مجموعة من المجموعات ذات التزاوج الداخلي.
3. تصنيف المتغيرات: تم تصنيف المتغيرات الممرضة إلى ثلاث مجموعات:
   • متغيرات ClinVar الممرضة/المحتمل أن تكون ممرضة (P/LP) دون وجود تفسيرات متضاربة.
   • متغيرات فقدان الوظيفة (LoF) عالية الثقة (HC) المستمدة من LOFTEE.
   • متغيرات مغلوطة (missense) يُتوقع أنها ضارة (AlphaMissense "محتمل أن تكون ممرضة" مع REVEL ≥ 0.75).
   • تم تطبيق فلاتر جودة صارمة، بما في ذلك فلتر للآثار الناتجة عن التماثل الجيني فقط، وسقف للتردد الأليلي (AF) عند 0.5% لاستبعاد المتغيرات الشائعة ضعيفة التأثير (مثل MTHFR c.677C>T).
4. النمذجة الجينية:
   • التزاوج العشوائي: تم حساب الانتشار باستخدام توازن هاردي-واينبرغ (HWE)، مع مراعاة كل من المتماثلين جينياً (homozygotes) والمتباينين جينياً (compound heterozygotes) ($GP = q^2$).
   • نمذجة زواج الأقوب (Consanguinity): إدراكاً للتغاير في زواج الأقارب في جنوب آسيا، قامت الدراسة بنمذجة الانتشار عبر طيف من معاملات وفيض رايت (Wright's inbreeding coefficients - $F$): 0 (عشوائي)، 0.01، 0.02 (متوسط السكان)، 0.03، و0.0625 (اتحاد أبناء العمومة المباشرين).
   • المعايرة: تمت معايرة التقديرات مقابل سجلات فحص حديثي الولادة (NBS) وتعديلها لتصحيح فجوات اكتشاف المتغيرات الهيكلية (CNVs/indels) التي لا يمكن اكتشافها بواسطة تسلسل القراءة القصيرة.
5. المقارنة عبر الأصول العرقية: تم تكرار المسار عبر ثماني مجموعات عرقية عالمية من gnomAD لتقييم انحياز التقصي ومقارنة ترددات الحامل للمرض.

النتائج الرئيسية

• عبء الحامل للمرض: في ظل التزاوج العشوائي، يحمل 84.2% من سكان جنوب آسيا أليلاً واحداً على الأقل من أصل 696 جيناً. ويمثل هذا أعلى تردد للحامل للمرض بين ثمانية أصول عرقية عالمية، يليه مباشرة الأفارقة (84.0%) والشرق آسيويين (83.5%)، وهو أعلى بكثير من المجموعات المؤسسة مثل الفنلنديين (64.8%) واليهود الأشكناز (63.6%).
• انتشار الولادات (التزاوج العشوائي): يبلغ الانتشار الجيني للولادات المقدر 255 لكل 100,000 مولود حي (1 من 392).
  • تفصيل حسب نوع المتغير: 40.6 لكل 100,000 (ClinVar P/LP)، و41.3 لكل 100,000 (HC LoF)، و62.7 لكل 100,000 (المتغيرات المغلوطة المتوقع ضررها).
  • التباين الجيني المركب (Compound Heterozygosity): ينشأ حوالي 43% من إجمالي العبء (110.5 لكل 100,000) من التباين الجيني المركب، مما يسلط الضوء على أهمية نمذجة المصفوفة ثنائية الأليل بدلاً من مجرد عد المتماثلين جينياً.
• تأثير زواج الأقارب: يزداد الانتشار بشكل كبير مع زيادة التزاوج الداخلي. تحت نموذج متوسط السكان ($F=0.02$)، يرتفع الانتشار إلى 2,072 لكل 100,000 (1 من 48). وبالنسبة لذرية اتحاد أبناء العمومة المباشرين ($F=0.0625$)، يصل المعدل إلى 5,824 لكل 100,000 (1 من 17).
• الجينات الأكثر مساهمة: لا يوجد جين واحد يستحوذ على معظم العبء. تشمل المساهمات الرئيسية: MAN2B2، DPYD، PYGL، ATP7B، PAH، MMUT، PYGM، DHCR24، SPG7، SACS، GLDC، وMLYCD.
• التحقق (GenomeIndia): في مجموعة GenomeIndia، بلغ الانتشار للمتغيرات المشتركة 139.6 لكل 100,000، وارتفع إلى 151.5 لكل 100,000 عند تضمين أليلات ClinVar P/LP الخاصة بـ GenomeIndia. ورغم أنه أقل من تقدير gnomAD لجنوب آسيا (ربما بسبب صغر حجم العينة وانخفاض اكتشاف الأليلات شديدة الندرة)، إلا أن أحمال الجينات لكل جين كانت مرتبطة بقوة (معامل سبيرمان $\rho = 0.83$). وقد قدم التقدير المجمع 189.1 لكل 100,000.
• انحياز التقصي: أظهر تحليل يعتمد على ClinVar فقط أن التقدير يقل بشكل كبير عن عبء جنوب آسيا (41 لكل 100,000) بينما يبالغ في تقدير عبء غير الفنلنديين الأوروبيين، مما يظهر انحياز تقصي مركزي أوروبي. وقد صحح إدراج المتغيرات المتوقع ضررها وفلاتر HC LoF هذا الاختلال.
• المعايرة: وفر تطبيق عامل قياس سريري-إلى-جيني مشتق من فحص حديثي الولادة (0.65×) مع مراعاة فجوات المتغيرات الهيكلية (1.59×) حدوداً استكشافية، حيث بلغ معدل الحدوث المعاير 165 لكل 100,000، والحد الأعلى 406 لكل 100,000.

الأهمية والادعاءات
يزعم البحث أن العبء الجماعي لأخطاء التمثيل الغذراي الوراثية المتنحية في جنوب آسيا كبير، حيث يحمل ما يقرب من 85% من السكان أليلاً واحداً على الأقل ممرضاً. وتؤكد الدراسة أنه بينما يكون الانتشار الجيني الأساسي تحت التزاوج العشوائي مشابهاً للمجموعات العالمية الأخرى، فإن الانتشار السريري يُقاد بشكل كبير بواسطة زواج الأقارب والتزاوج الداخلي، مما يمكن أن يزيد من خطر الولادات المصابة بمقدار عشرة أضعاف.

تشمل المساهمات الرئيسية ما يلي:

1. تصحيح الانحياز: إثبات أن الاعتماد على ClinVar وحده يؤدي إلى انحياز مركزي أوروبي، مما يقلل من تقدير مخاطر جنوب آسيا، وأن دمج المتنبئات الحاسوبية (AlphaMissense, REVEL) وفلاتر HC LoF يوفر نموذجاً جينياً أكثر إنصافاً ودقة.
2. النمذجة التفصيلية: تجاوز تقديم تقدير واحد للانتشار عبر توفير نطاق من السيناريوهات بناءً على معاملات التزاوج الداخلي ($F$)، مع الاعتراف بتغاير المجموعات السكانية في جنوب آسيا.
3. الفائدة السريرية: تحديد الجينات ذات العبء العالي (مثل DPYD، ATP7B، PAH) التي يجب منحها الأولوية لبرامج فحص الحامل للمرض الموسعة وبرامج فحص حديثي الولادة المصممة خصيصاً لسكان جنوب آسيا، مع ملاحظة أن اللوحات المركزية الأوروبية قد تكون ذات حساسية منخفضة.
4. قابلية التكرار: توفير إطار عمل قابل للتكرار بالكامل، بما في ذلك الكود والبيانات وجداول مستوى المتغيرات، لتسهيل المزيد من الأبحاث وتطوير السياسات.

يخلص المؤلف إلى أن هذه النتائج توفر قاعدة أدلة حاسمة للسياسة الجينية، واستراتيجيات فحص الحامل للمرض، وتوسيع فحص حديثي الولادة عبر جنوب آسيا، مؤكداً على الحاجة إلى نهج خاص بالسكان لإدارة العبء المرتفع للاضطرابات الأيضية الوراثية المتنحية.