يُعدّ الرشّ الموضعي داخل المنازل (IRS) الركيزة الأساسية لجهود مكافحة نواقل داء الليشمانيات الحشوي (VL) في الهند. ولا تزال المعلومات شحيحة حول تأثير هذه المكافحة على مختلف أنواع الأسر. في هذه الدراسة، نقوم بتقييم ما إذا كان للرشّ الموضعي باستخدام المبيدات الحشرية نفس التأثيرات المتبقية والتدخلية على جميع أنواع الأسر في إحدى القرى. كما قمنا بتطوير خرائط مخاطر مكانية مُدمجة ونماذج تحليل كثافة البعوض، استنادًا إلى خصائص الأسر، وحساسية البعوض للمبيدات، وحالة الرشّ الموضعي، وذلك لدراسة التوزيع المكاني والزماني للنواقل على المستوى الجزئي.
أُجريت الدراسة في قريتين تابعتين لبلدة ماهنار في مقاطعة فايشالي بولاية بيهار. وتم تقييم مكافحة ناقلات داء الليشمانيا الحشوي (P. argentipes) باستخدام الرش الموضعي للمبيدات الحشرية [ثنائي كلورو ثنائي فينيل ثلاثي كلورو الإيثان (DDT 50%) والبيرثرويدات الاصطناعية (SP 5%)]. وقُيِّمَت فعالية المبيدات الحشرية المتبقية على أنواع مختلفة من الجدران باستخدام طريقة الفحص الحيوي المخروطي وفقًا لتوصيات منظمة الصحة العالمية. كما تم فحص حساسية سمكة الفضة المحلية للمبيدات الحشرية باستخدام فحص حيوي مخبري. ورُصِدَت كثافة البعوض قبل وبعد الرش الموضعي في المنازل وملاجئ الحيوانات باستخدام مصائد ضوئية نصبتها مراكز مكافحة الأمراض والوقاية منها من الساعة 6:00 مساءً إلى 6:00 صباحًا. وتم تطوير النموذج الأمثل لتحليل كثافة البعوض باستخدام تحليل الانحدار اللوجستي المتعدد. تم استخدام تقنية التحليل المكاني القائمة على نظم المعلومات الجغرافية لرسم خريطة لتوزيع حساسية ناقلات الآفات للمبيدات حسب نوع الأسرة، وتم استخدام حالة الرش المتبقي داخل المنزل لشرح التوزيع المكاني والزماني للروبيان الفضي.
تُعدّ بعوضة الفضة شديدة الحساسية لمبيد SP (100%)، لكنها تُظهر مقاومة عالية لمبيد DDT، حيث بلغت نسبة الوفيات 49.1%. وقد أُفيد بأنّ رشّ SP الموضعي يحظى بقبولٍ عامّ أفضل من رشّ DDT الموضعي في جميع أنواع الأسر. وتفاوتت فعالية المبيدات المتبقية باختلاف أسطح الجدران؛ ولم يستوفِ أيٌّ من المبيدات الحشرية مدة التأثير الموصى بها من قِبل منظمة الصحة العالمية للرشّ الموضعي. وفي جميع نقاط القياس بعد الرشّ الموضعي، كان انخفاض أعداد بقّ الرائحة الكريهة الناتج عن رشّ SP الموضعي أكبر بين مجموعات الأسر (أي الرشاشين والمراقبين) مقارنةً برشّ DDT الموضعي. وتُظهر خريطة المخاطر المكانية المُجمّعة أنّ رشّ SP الموضعي يُحقق تأثيرًا أفضل في مكافحة البعوض من رشّ DDT الموضعي في جميع مناطق المخاطر لأنواع الأسر المختلفة. وحدّد تحليل الانحدار اللوجستي متعدد المستويات خمسة عوامل خطر مرتبطة ارتباطًا وثيقًا بكثافة روبيان الفضة.
ستوفر النتائج فهمًا أفضل لممارسات الرش المتبقي داخل المنازل في مكافحة داء الليشمانيات الحشوي في بيهار، مما قد يساعد في توجيه الجهود المستقبلية لتحسين الوضع.
داء الليشمانيات الحشوي، المعروف أيضًا باسم الكالازار، هو مرض استوائي مهمل متوطن ينتقل عن طريق الحشرات، ويسببه طفيليات أولية من جنس الليشمانيا. في شبه القارة الهندية، حيث يُعد الإنسان المضيف الوحيد للطفيلي، ينتقل الطفيل (الليشمانيا الدونوفانية) إلى الإنسان عن طريق لدغات إناث البعوض المصابة (فليبوتوموس أرجنتيبس) [1، 2]. في الهند، ينتشر داء الليشمانيات الحشوي بشكل رئيسي في أربع ولايات وسطى وشرقية: بيهار، وجارخاند، والبنغال الغربية، وأوتار براديش. كما تم الإبلاغ عن بعض حالات التفشي في ماديا براديش (وسط الهند)، وغوجارات (غرب الهند)، وتاميل نادو، وكيرالا (جنوب الهند)، بالإضافة إلى مناطق سفوح جبال الهيمالايا في شمال الهند، بما في ذلك هيماشال براديش وجامو وكشمير [3]. من بين الولايات الموبوءة، تُعدّ ولاية بيهار الأكثر توطناً للمرض، حيث تُصاب 33 مقاطعة بداء الليشمانيا الحشوي، ما يُمثل أكثر من 70% من إجمالي الحالات في الهند سنوياً [4]. ويُقدّر عدد المُعرّضين للخطر في المنطقة بنحو 99 مليون شخص، بمتوسط إصابة سنوي يبلغ 6752 حالة (2013-2017).
في ولاية بيهار ومناطق أخرى من الهند، تعتمد جهود مكافحة داء الليشمانيا الحشوي على ثلاث استراتيجيات رئيسية: الكشف المبكر عن الحالات، والعلاج الفعال، ومكافحة النواقل باستخدام الرش الداخلي بالمبيدات الحشرية في المنازل وملاجئ الحيوانات [4، 5]. وكنتيجة ثانوية لحملات مكافحة الملاريا، نجح الرش الداخلي بالمبيدات الحشرية في السيطرة على داء الليشمانيا الحشوي في ستينيات القرن الماضي باستخدام ثنائي كلورو ثنائي فينيل ثلاثي كلورو الإيثان (DDT 50% مسحوق قابل للبلل، 1 غرام من المادة الفعالة/م²)، كما نجحت المكافحة البرنامجية في السيطرة على داء الليشمانيا الحشوي في عامي 1977 و1992 [5، 6]. ومع ذلك، أكدت دراسات حديثة أن الروبيان الفضي البطن قد طور مقاومة واسعة النطاق لمادة DDT [4، 7، 8]. في عام 2015، استبدل البرنامج الوطني لمكافحة الأمراض المنقولة بالنواقل (NVBDCP، نيودلهي) الرش الداخلي بالمبيدات الحشرية من DDT إلى البيريثرويدات الاصطناعية (SP؛ ألفا-سيبرمثرين 5% مسحوق قابل للبلل، 25 ملغ من المادة الفعالة/م²) [7، 9]. حددت منظمة الصحة العالمية هدفًا للقضاء على داء الليشمانيا الحشوي بحلول عام 2020 (أي أقل من حالة واحدة لكل 10,000 شخص سنويًا على مستوى الشارع/الحي) [10]. وقد أظهرت العديد من الدراسات أن الرش المتبقي داخل المنازل أكثر فعالية من طرق مكافحة النواقل الأخرى في تقليل كثافة ذباب الرمل [11، 12، 13]. ويتوقع نموذج حديث أيضًا أنه في المناطق ذات معدلات الإصابة المرتفعة (أي معدل الإصابة قبل المكافحة 5 حالات لكل 10,000 شخص)، يمكن للرش المتبقي الفعال الذي يغطي 80% من الأسر أن يحقق أهداف القضاء على المرض قبل ذلك بسنة إلى ثلاث سنوات [14]. يؤثر داء الليشمانيا الحشوي على أفقر المجتمعات الريفية في المناطق الموبوءة، وتعتمد مكافحة نواقل المرض فيها كليًا على الرش المتبقي داخل المنازل، ولكن لم تتم دراسة الأثر المتبقي لهذا الإجراء على مختلف أنواع الأسر ميدانيًا في مناطق التدخل [15، 16]. إضافة إلى ذلك، وبعد الجهود المكثفة لمكافحة داء الليشمانيا الحشوي، استمر الوباء في بعض القرى لعدة سنوات وتحول إلى بؤر ساخنة [17]. لذا، من الضروري تقييم الأثر المتبقي للرش المتبقي داخل المنازل على رصد كثافة البعوض في مختلف أنواع المساكن. إضافةً إلى ذلك، يُسهم رسم خرائط المخاطر الجغرافية المكانية على نطاق دقيق في فهم أعداد البعوض والسيطرة عليها بشكل أفضل حتى بعد التدخل. تُعدّ نظم المعلومات الجغرافية (GIS) مزيجًا من تقنيات رسم الخرائط الرقمية التي تُمكّن من تخزين وتراكب ومعالجة وتحليل واسترجاع وعرض مجموعات مختلفة من البيانات الجغرافية والبيئية والاجتماعية الديموغرافية لأغراض متنوعة [18، 19، 20]. يُستخدم نظام تحديد المواقع العالمي (GPS) لدراسة الموقع المكاني لمكونات سطح الأرض [21، 22]. وقد طُبّقت أدوات وتقنيات النمذجة المكانية القائمة على نظم المعلومات الجغرافية ونظام تحديد المواقع العالمي على العديد من الجوانب الوبائية، مثل التقييم المكاني والزماني للأمراض والتنبؤ بتفشيها، وتنفيذ استراتيجيات المكافحة وتقييمها، وتفاعلات مسببات الأمراض مع العوامل البيئية، ورسم خرائط المخاطر المكانية [20، 23، 24، 25، 26]. يمكن للمعلومات التي يتم جمعها واستخلاصها من خرائط المخاطر الجغرافية المكانية أن تسهل اتخاذ تدابير تحكم فعالة وفي الوقت المناسب.
قيّمت هذه الدراسة الفعالية المتبقية وتأثير تدخلات DDT وSP-IRS على مستوى الأسر المعيشية في إطار البرنامج الوطني لمكافحة نواقل داء الليشمانيات الحشوي في ولاية بيهار بالهند. وشملت الأهداف الإضافية تطوير نموذج مُدمج لخريطة المخاطر المكانية وتحليل كثافة البعوض، استنادًا إلى خصائص المساكن، وحساسية نواقل الأمراض للمبيدات الحشرية، وحالة الرش الموضعي للأسر المعيشية، وذلك لدراسة التسلسل الهرمي للتوزيع المكاني والزماني للبعوض على نطاق صغير.
أُجريت الدراسة في منطقة ماهنار التابعة لمقاطعة فايشالي على الضفة الشمالية لنهر الغانج (الشكل 1). تُعدّ ماهنار منطقةً شديدة التوطن، حيث يبلغ متوسط حالات داء الليشمانيا الحشوي فيها 56.7 حالة سنويًا (170 حالة في الفترة 2012-2014)، ويبلغ معدل الإصابة السنوي 2.5-3.7 حالة لكل 10,000 نسمة. تم اختيار قريتين: تشاكيسو كموقع مرجعي (الشكل 1د1؛ لم تُسجّل أي حالات داء الليشمانيا الحشوي خلال السنوات الخمس الماضية)، ولافابور ماهنار كموقع متوطن (الشكل 1د2؛ شديد التوطن، حيث سُجّلت 5 حالات أو أكثر لكل 1000 نسمة سنويًا خلال السنوات الخمس الماضية). تم اختيار القريتين بناءً على ثلاثة معايير رئيسية: الموقع وسهولة الوصول (أي تقع على نهر يسهل الوصول إليه على مدار السنة)، والخصائص الديموغرافية، وعدد الأسر (أي 200 أسرة على الأقل؛ تضم تشاكيسو 202 أسرة و204 أسر بمتوسط حجم أسرة). بلغ عدد سكان قريتي الدراسة 4.9 و5.1 نسمة على التوالي، بينما بلغ عدد سكان قرية لافابور ماهانار 4.9 و5.1 نسمة على التوالي، وتم تحديد نوع الأسرة (أسرة ذات سكن مشترك) وطبيعة توزيعها (أي توزيع عشوائي لأسرة ذات سكن مشترك مختلطة). تقع كلتا القريتين ضمن مسافة 500 متر من مدينة ماخنار ومستشفى المقاطعة. أظهرت الدراسة أن سكان القريتين كانا يشاركان بنشاط كبير في الأنشطة البحثية. تتكون المنازل في قرية التدريب [من غرفة نوم واحدة إلى غرفتين مع شرفة ملحقة، ومطبخ واحد، وحمام واحد، وحظيرة واحدة (ملحقة أو منفصلة)] من جدران من الطوب/الطين وأرضيات من الطوب اللبن، وجدران من الطوب مع طلاء من الجص الجيري الإسمنتي، وأرضيات من الإسمنت، وجدران من الطوب غير مطلية وغير مطلية، وأرضيات من الطين، وسقف من القش. تتمتع منطقة فايشالي بأكملها بمناخ شبه استوائي رطب مع موسم ممطر (من يوليو إلى أغسطس) وموسم جاف (من نوفمبر إلى ديسمبر). يبلغ متوسط هطول الأمطار السنوي 720.4 ملم (يتراوح بين 736.5 و1076.7 ملم)، والرطوبة النسبية 65±5% (تتراوح بين 16 و79%)، ومتوسط درجة الحرارة الشهرية 17.2-32.4 درجة مئوية. ويُعد شهرا مايو ويونيو الأكثر دفئًا (درجات الحرارة 39-44 درجة مئوية)، بينما يُعد شهر يناير الأبرد (7-22 درجة مئوية).
تُظهر خريطة منطقة الدراسة موقع ولاية بيهار على خريطة الهند (أ)، وموقع مقاطعة فايشالي على خريطة بيهار (ب). تم اختيار قريتين للدراسة في منطقة ماخنار (ج): تشاكيسو كموقع مرجعي، ولافابور ماخنار كموقع تجريبي.
في إطار البرنامج الوطني لمكافحة داء الليشمانيا الحشوي، نفّذ مجلس الصحة المجتمعية في بيهار جولتين من الرش المتبقي داخل المنازل خلال عامي 2015 و2016 (الجولة الأولى: فبراير-مارس، والجولة الثانية: يونيو-يوليو)[4]. ولضمان التنفيذ الفعال لجميع أنشطة الرش المتبقي داخل المنازل، أعدّ معهد راجندرا التذكاري الطبي (بيهار)، باتنا، وهو معهد تابع للمجلس الهندي للأبحاث الطبية (نيودلهي)، خطة عمل تفصيلية. وتم اختيار القرى المشاركة في الرش المتبقي داخل المنازل بناءً على معيارين رئيسيين: وجود حالات سابقة من داء الليشمانيا الحشوي وداء الليشمانيا الحشوي الرجعي في القرية (أي القرى التي سُجّلت فيها حالة واحدة أو أكثر خلال أي فترة زمنية في السنوات الثلاث الماضية، بما في ذلك سنة التنفيذ). القرى غير الموبوءة المحيطة بـ"البؤر الساخنة" (أي القرى التي سجلت حالات إصابة بشكل مستمر لمدة سنتين أو أكثر، أو حالتين أو أكثر لكل 1000 نسمة)، والقرى الموبوءة حديثًا (التي لم تسجل أي حالات إصابة خلال السنوات الثلاث الماضية)، والقرى التي تم الإبلاغ عنها في السنة الأخيرة من سنة التنفيذ المذكورة في المرجع [17]. كما تم إدراج القرى المجاورة التي طبقت الجولة الأولى من نظام الضرائب الوطني، والقرى الجديدة أيضًا، في الجولة الثانية من خطة عمل الضرائب الوطنية. في عام 2015، نُفذت جولتان من الرش المتبقي داخل المنازل باستخدام مادة DDT (50% مسحوق قابل للبلل، 1 غرام من المادة الفعالة/م²) في قرى الدراسة التدخلية. ومنذ عام 2016، يُنفذ الرش المتبقي داخل المنازل باستخدام البيريثرويدات الاصطناعية (SP؛ ألفا-سيبرمثرين 5% مسحوق قابل للبخار، 25 ملغ من المادة الفعالة/م²). تم الرش باستخدام مضخة هادسون إكسبرت (13.4 لترًا) مزودة بفلتر ضغط وصمام تدفق متغير (1.5 بار) وفوهة رش مسطحة 8002 للأسطح المسامية [27]. راقب المجلس الهندي للأبحاث الطبية - معهد البحوث الطبية الإقليمي للعلوم الطبية، باتنا (بيهار) برنامج الرش المتبقي داخل المنازل على مستوى الأسر والقرى، وقدم معلومات أولية عنه لسكان القرى عبر مكبرات الصوت خلال اليومين الأولين. تم تجهيز كل فريق رش بجهاز مراقبة (مقدم من معهد البحوث الطبية الإقليمي للعلوم الطبية) لمتابعة أداء الفريق. وتم نشر أمناء المظالم، إلى جانب فرق الرش، في جميع المنازل لإعلام رب الأسرة وطمأنته بشأن الآثار الإيجابية للرش. خلال جولتين من مسوحات الرش، بلغت التغطية الإجمالية للأسر في القرى المدروسة 80% على الأقل [4]. تم تسجيل حالة الرش (أي عدم الرش، والرش الجزئي، والرش الكامل؛ كما هو موضح في الملف الإضافي 1: الجدول S1) لجميع الأسر في القرية المستهدفة خلال جولتي الرش.
أُجريت الدراسة في الفترة من يونيو 2015 إلى يوليو 2016. استخدمت عملية الرش الموضعي للمبيدات مراكز مكافحة الأمراض لرصد كثافة ذباب الرمل، والسيطرة عليها، والوقاية منها في كل جولة من جولات الرش، وذلك قبل التدخل (أي قبل أسبوعين من التدخل؛ مسح أساسي) وبعده (أي بعد أسبوعين، وأربعة أسابيع، واثني عشر أسبوعًا من التدخل؛ مسوحات متابعة). وتم تركيب مصائد ضوئية في غرف النوم وملاجئ الحيوانات في كل منزل. في القرية التي أُجريت فيها الدراسة، تم فحص كثافة ذباب الرمل في 48 منزلًا قبل بدء الرش (12 منزلًا يوميًا لمدة أربعة أيام متتالية حتى اليوم السابق ليوم الرش). تم اختيار 12 منزلًا لكل مجموعة من المجموعات الأربع الرئيسية للمنازل (أي المنازل ذات الجدران الطينية العادية، والمنازل ذات الجدران الإسمنتية المغطاة بالجير، والمنازل المبنية من الطوب غير المطلي وغير المكسو بالجص، والمنازل ذات الأسقف المصنوعة من القش). بعد ذلك، تم اختيار 12 أسرة فقط (من أصل 48 أسرة قبل تطبيق الرش الموضعي للمبيدات) لمواصلة جمع بيانات كثافة البعوض بعد اجتماع الرش. ووفقًا لتوصيات منظمة الصحة العالمية، تم اختيار 6 أسر من مجموعة التدخل (الأسر التي تلقت معالجة الرش الموضعي) ومجموعة المراقبة (الأسر في قرى التدخل، أي تلك التي رفضت السماح بالرش) [28]. أما في مجموعة الضبط (الأسر في القرى المجاورة التي لم تخضع للرش لعدم وجود داء الليشمانيا الحشوي)، فقد تم اختيار 6 أسر فقط لمراقبة كثافة البعوض قبل وبعد جلستين من الرش. بالنسبة لجميع مجموعات مراقبة كثافة البعوض الثلاث (أي التدخل، والمراقبة، والضبط)، تم اختيار الأسر من ثلاث مجموعات لمستويات الخطورة (منخفضة، ومتوسطة، وعالية؛ أسرتان من كل مستوى خطورة)، وتم تصنيف خصائص مخاطر ارتفاع درجة الحرارة (الوحدات والهياكل موضحة في الجدول 1 والجدول 2 على التوالي) [29، 30]. تم اختيار أسرتين لكل مستوى خطورة لتجنب التحيز في تقديرات كثافة البعوض والمقارنات بين المجموعات. في المجموعة التجريبية، تم رصد كثافة البعوض بعد الرش المتبقي في نوعين من الأسر التي خضعت للرش المتبقي: الأسر المعالجة بالكامل (ن = 3؛ أسرة واحدة لكل مستوى من مستويات مجموعة المخاطر) والأسر المعالجة جزئيًا (ن = 3؛ أسرة واحدة لكل مستوى من مستويات مجموعة المخاطر).
نُقلت جميع البعوضات التي جُمعت في أنابيب الاختبار من الحقل إلى المختبر، حيث تم قتلها باستخدام قطن مبلل بالكلوروفورم. حُدد جنس ذباب الرمل الفضي وفُصل عن الحشرات الأخرى والبعوض بناءً على خصائصه المورفولوجية باستخدام رموز تعريف قياسية [31]. ثم عُلبت جميع ذكور وإناث ذباب الرمل الفضي بشكل منفصل في كحول بنسبة 80%. حُسبت كثافة البعوض لكل مصيدة/ليلة باستخدام الصيغة التالية: إجمالي عدد البعوض الذي جُمع / عدد المصائد الضوئية الموضوعة في الليلة. قُدرت النسبة المئوية للتغير في وفرة البعوض (SFC) نتيجة الرش المتبقي داخل المنازل باستخدام DDT وSP باستخدام الصيغة التالية [32]:
حيث A هو متوسط SFC الأساسي للأسر المتدخلة، وB هو متوسط SFC الخاص بمصلحة الضرائب الأمريكية للأسر المتدخلة، وC هو متوسط SFC الأساسي للأسر الضابطة/المراقبة، وD هو متوسط SFC للأسر الضابطة/المراقبة التابعة لمصلحة الضرائب الأمريكية.
تشير نتائج تأثير التدخل، المسجلة كقيم سالبة وموجبة، إلى انخفاض وزيادة في مستوى السكر في الدم بعد الرش الموضعي، على التوالي. إذا بقي مستوى السكر في الدم بعد الرش الموضعي مماثلاً لمستوى السكر في الدم الأساسي، يُحسب تأثير التدخل على أنه صفر.
وفقًا لبرنامج تقييم المبيدات التابع لمنظمة الصحة العالمية (WHOPES)، تم تقييم حساسية روبيان الفضة المحلي لمبيدي DDT وSP باستخدام اختبارات حيوية معيارية في المختبر [33]. تم تعريض إناث روبيان الفضة السليمة وغير المتغذية (18-25 روبيان لكل مجموعة) لمبيدات تم الحصول عليها من جامعة العلوم الماليزية (USM، ماليزيا؛ بتنسيق من منظمة الصحة العالمية) باستخدام مجموعة اختبار حساسية المبيدات التابعة لمنظمة الصحة العالمية [4، 9، 33، 34]. تم اختبار كل مجموعة من الاختبارات الحيوية للمبيدات ثماني مرات (أربع تكرارات، كل منها يتم تشغيله بالتزامن مع المجموعة الضابطة). أُجريت اختبارات المجموعة الضابطة باستخدام ورق مشبع مسبقًا بمادة ريسيلا (لمبيد DDT) وزيت السيليكون (لمبيد SP) مقدم من جامعة العلوم الماليزية. بعد 60 دقيقة من التعرض، وُضعت البعوضات في أنابيب منظمة الصحة العالمية، ووُضعت عليها قطنات ماصة مشبعة بمحلول سكري بنسبة 10%. تمت ملاحظة عدد البعوضات النافقة بعد ساعة واحدة، ومعدل النفوق النهائي بعد 24 ساعة. يُصنَّف وضع المقاومة وفقًا لإرشادات منظمة الصحة العالمية: تشير نسبة الوفيات من 98% إلى 100% إلى قابلية الإصابة، وتشير نسبة الوفيات من 90% إلى 98% إلى احتمال وجود مقاومة تتطلب تأكيدًا، وتشير نسبة الوفيات الأقل من 90% إلى وجود مقاومة [33، 34]. ولأن نسبة الوفيات في المجموعة الضابطة تراوحت بين 0% و5%، لم يتم إجراء أي تعديل على معدل الوفيات.
تم تقييم الفعالية البيولوجية والآثار المتبقية للمبيدات الحشرية على النمل الأبيض المحلي في ظروف الحقل. في ثلاث منازل تجريبية (منزل واحد لكل نوع من المنازل ذات الجدران المطلية بالجص الطيني العادي أو PMP، ومنزل واحد لكل نوع من المنازل المطلية بالجص الإسمنتي والجير أو CPLC، ومنزل واحد من الطوب غير المطلي وغير المدهون أو BUU) بعد أسبوعين، وأربعة أسابيع، واثني عشر أسبوعًا من الرش. أُجري اختبار بيولوجي معياري لمنظمة الصحة العالمية على مخاريط تحتوي على مصائد ضوئية. [27، 32]. تم استبعاد التدفئة المنزلية نظرًا لعدم استواء الجدران. في كل تحليل، استُخدم 12 مخروطًا في جميع المنازل التجريبية (أربعة مخاريط لكل منزل، مخروط واحد لكل نوع من أنواع أسطح الجدران). تم تثبيت المخاريط على كل جدار من جدران الغرفة على ارتفاعات مختلفة: مخروط واحد على مستوى الرأس (من 1.7 إلى 1.8 متر)، ومخروطان على مستوى الخصر (من 0.9 إلى 1 متر)، ومخروط واحد أسفل الركبة (من 0.3 إلى 0.5 متر). وُضِعَت عشر بعوضات أنثى غير متغذية (عشر بعوضات لكل مخروط؛ جُمِعَت من قطعة أرض ضابطة باستخدام جهاز شفط) في كل حجرة مخروطية بلاستيكية تابعة لمنظمة الصحة العالمية (مخروط واحد لكل نوع من أنواع المنازل) كعناصر ضابطة. بعد 30 دقيقة من التعرض، أُزيلت البعوضات بعناية من الحجرة المخروطية باستخدام جهاز شفط ذي رأس كوعي، ونُقِلَت إلى أنابيب تابعة لمنظمة الصحة العالمية تحتوي على محلول سكري بنسبة 10% للتغذية. سُجِّلَ معدل النفوق النهائي بعد 24 ساعة عند درجة حرارة 27 ± 2 درجة مئوية ورطوبة نسبية 80 ± 10%. عُدِّلَت معدلات النفوق التي تتراوح درجاتها بين 5% و20% باستخدام معادلة أبوت [27] كما يلي:
حيث P هي نسبة الوفيات المعدلة، وP1 هي نسبة الوفيات المرصودة، وC هي نسبة الوفيات في المجموعة الضابطة. تم استبعاد التجارب التي تجاوزت فيها نسبة الوفيات في المجموعة الضابطة 20% وإعادة إجرائها [27، 33].
أُجري مسح شامل للأسر في القرية المستهدفة. سُجِّل موقع كل أسرة باستخدام نظام تحديد المواقع العالمي (GPS)، بالإضافة إلى تصميمها ونوع مواد بنائها، وموقعها، وحالة التدخل. طوّرت منصة نظم المعلومات الجغرافية (GIS) قاعدة بيانات جغرافية رقمية تتضمن طبقات حدودية على مستوى القرية والمنطقة والولاية. رُسمت علامات جغرافية لجميع مواقع الأسر باستخدام طبقات نقاط نظم المعلومات الجغرافية على مستوى القرية، ورُبطت معلوماتها الوصفية وحُدِّثت. في كل موقع أسرة، قُيِّمَت المخاطر بناءً على درجة الحرارة العالية، وحساسية نواقل المبيدات الحشرية، وحالة الرش المتبقي داخل المنازل (الجدول 1) [11، 26، 29، 30]. ثم حُوِّلت جميع نقاط مواقع الأسر إلى خرائط موضوعية باستخدام تقنية الاستيفاء المكاني (الترجيح العكسي للمسافة، ورسم الخرائط الالتفافية المكعبة) [35]. تم إنشاء نوعين من خرائط المخاطر المكانية الموضوعية: خرائط موضوعية قائمة على تقنية HT، وخرائط موضوعية لحساسية ناقلات المبيدات وحالة الرش المتبقي داخل المنازل (ISV وIRSS). ثم دُمجت خريطتا المخاطر الموضوعية باستخدام تحليل التراكب الموزون [36]. خلال هذه العملية، أُعيد تصنيف طبقات البيانات النقطية إلى فئات تفضيل عامة لمستويات المخاطر المختلفة (أي، مخاطر عالية، ومتوسطة، ومنخفضة/معدومة). بعد ذلك، ضُربت كل طبقة بيانات نقطية مُعاد تصنيفها بالوزن المُخصص لها بناءً على الأهمية النسبية للمعايير التي تدعم وفرة البعوض (استنادًا إلى الانتشار في قرى الدراسة، ومواقع تكاثر البعوض، وسلوك الراحة والتغذية) [26، 29، 30، 37]. وُزنت خريطتا المخاطر الموضوعية بنسبة 50:50 نظرًا لمساهمتهما المتساوية في وفرة البعوض (الملف الإضافي 1: الجدول S2). بجمع خرائط التراكب الموضوعية الموزونة، يتم إنشاء خريطة مخاطر مركبة نهائية وعرضها على منصة نظم المعلومات الجغرافية (GIS). يتم عرض خريطة المخاطر النهائية ووصفها من حيث قيم مؤشر مخاطر ذبابة الرمل (SFRI) المحسوبة باستخدام الصيغة التالية:
في الصيغة، يُمثل P قيمة مؤشر المخاطر، وL قيمة المخاطر الإجمالية لموقع كل أسرة، وH أعلى قيمة مخاطر لأسرة في منطقة الدراسة. قمنا بإعداد طبقات نظم المعلومات الجغرافية وتحليلها باستخدام برنامج ESRI ArcGIS الإصدار 9.3 (ريدلاندز، كاليفورنيا، الولايات المتحدة الأمريكية) لإنشاء خرائط المخاطر.
أجرينا تحليلات انحدار متعددة لدراسة التأثيرات المشتركة لتقنيات HT وISV وIRSS (كما هو موضح في الجدول 1) على كثافة البعوض داخل المنازل (ن = 24). واعتُبرت خصائص المساكن وعوامل الخطر المستندة إلى تدخل الرش المتبقي داخل المنازل، والمسجلة في الدراسة، متغيرات تفسيرية، بينما استُخدمت كثافة البعوض كمتغير استجابة. وأُجريت تحليلات انحدار بواسون أحادية المتغير لكل متغير تفسيري مرتبط بكثافة ذباب الرمل. وخلال التحليل الأحادي المتغير، أُزيلت المتغيرات غير الدالة إحصائيًا والتي كانت قيمة P لها أكبر من 15% من تحليل الانحدار المتعدد. ولدراسة التفاعلات، أُدرجت حدود التفاعل لجميع التوليفات الممكنة للمتغيرات الدالة إحصائيًا (الموجودة في التحليل الأحادي المتغير) في تحليل الانحدار المتعدد، ثم أُزيلت الحدود غير الدالة إحصائيًا من النموذج تدريجيًا للوصول إلى النموذج النهائي.
أُجري تقييم المخاطر على مستوى الأسر بطريقتين: تقييم المخاطر على مستوى الأسر، وتقييم مكاني مُدمج لمناطق الخطر على الخريطة. قُدّرت تقديرات المخاطر على مستوى الأسر باستخدام تحليل الارتباط بين تقديرات المخاطر وكثافة ذباب الرمل (التي جُمعت من 6 أسر مراقبة و6 أسر خاضعة للتدخل؛ قبل أسابيع من تطبيق الرش المتبقي داخل المنازل وبعدها). قُدّرت مناطق الخطر المكانية باستخدام متوسط عدد البعوض الذي جُمع من مختلف الأسر، وقورنت بين مجموعات الخطر (أي مناطق الخطر المنخفض والمتوسط والعالي). في كل جولة من جولات الرش المتبقي داخل المنازل، اختيرت 12 أسرة عشوائيًا (4 أسر في كل مستوى من مستويات مناطق الخطر الثلاثة؛ تُجرى عمليات جمع ليلية كل أسبوعين و4 أسابيع و12 أسبوعًا بعد الرش المتبقي داخل المنازل) لجمع البعوض لاختبار خريطة المخاطر الشاملة. استُخدمت بيانات الأسر نفسها (أي درجة الحرارة، ومؤشر شدة الإصابة، ومؤشر شدة الإصابة داخل المنازل، ومتوسط كثافة البعوض) لاختبار نموذج الانحدار النهائي. أُجري تحليل ارتباط بسيط بين الملاحظات الميدانية وكثافة البعوض المتوقعة من النموذج في الأسر.
تم حساب الإحصاءات الوصفية، مثل المتوسط والحد الأدنى والحد الأقصى وفترات الثقة 95% والنسب المئوية، لتلخيص البيانات المتعلقة بعلم الحشرات والرش المتبقي داخل المنازل. وتم حساب متوسط عدد/كثافة ومعدل نفوق حشرات البق الفضي (بقايا المبيدات الحشرية) باستخدام اختبارات بارامترية (اختبار t للعينات المزدوجة للبيانات ذات التوزيع الطبيعي) واختبارات لا بارامترية (اختبار ويلكوكسون للإشارة المرتبة) لمقارنة فعالية المبيدات بين أنواع الأسطح في المنازل (مثل: BUU مقابل CPLC، وBUU مقابل PMP، وCPLC مقابل PMP للبيانات غير ذات التوزيع الطبيعي). أُجريت جميع التحليلات باستخدام برنامج SPSS الإصدار 20 (شركة SPSS، شيكاغو، إلينوي، الولايات المتحدة الأمريكية).
تم حساب نسبة تغطية الأسر في القرى المستهدفة خلال جولات الرش الموضعي بمبيد DDT وSP. تلقى ما مجموعه 205 أسر الرش الموضعي في كل جولة، منها 179 أسرة (87.3%) في جولة DDT و194 أسرة (94.6%) في جولة SP لمكافحة ناقل داء الليشمانيا الحشوي. وكانت نسبة الأسر التي عولجت بالكامل بالمبيدات أعلى خلال جولة SP-IRS (86.3%) مقارنةً بجولة DDT-IRS (52.7%). بلغ عدد الأسر التي رفضت الرش الموضعي خلال جولة DDT 26 أسرة (12.7%)، بينما بلغ عددها خلال جولة SP 11 أسرة (5.4%). خلال جولتي DDT وSP، بلغ عدد الأسر المسجلة التي عولجت جزئيًا 71 أسرة (34.6% من إجمالي الأسر المعالجة) و17 أسرة (8.3% من إجمالي الأسر المعالجة) على التوالي.
وفقًا لإرشادات منظمة الصحة العالمية بشأن مقاومة المبيدات، كان جمبري الفضة في موقع التدخل شديد الحساسية لمادة ألفا-سيبرمثرين (0.05%)، حيث بلغ متوسط معدل النفوق المُسجل خلال التجربة (24 ساعة) 100%. وبلغ معدل الإغماء المُلاحظ 85.9% (95% CI: 81.1–90.6%). أما بالنسبة لمادة DDT، فقد بلغ معدل الإغماء بعد 24 ساعة 22.8% (95% CI: 11.5–34.1%)، وبلغ متوسط معدل النفوق في الاختبار الإلكتروني 49.1% (95% CI: 41.9–56.3%). وأظهرت النتائج أن جمبري الفضة قد طور مقاومة كاملة لمادة DDT في موقع التدخل.
يلخص الجدول 3 نتائج التحليل البيولوجي للأقماع لأنواع مختلفة من الأسطح (فترات زمنية مختلفة بعد الرش المتبقي) المعالجة بمادة DDT و SP. أظهرت بياناتنا أنه بعد 24 ساعة، انخفضت معدلات النفوق بشكل مطرد مع مرور الوقت لكلا المبيدين الحشريين (BUU مقابل CPLC: t(2) = -6.42، P = 0.02؛ BUU مقابل PMP: t(2) = 0.25، P = 0.83؛ CPLC مقابل PMP: t(2) = 1.03، P = 0.41 (لـ DDT-IRS وBUU) CPLC: t(2) = -5.86، P = 0.03 وPMP: t(2) = 1.42، P = 0.29؛ IRS وCPLC وPMP: t(2) = 3.01، P = 0.10 وSP: t(2) = 9.70، P = 0.01). بالنسبة لـ SP-IRS: بعد أسبوعين من الرش لجميع أنواع الجدران (أي بلغت نسبة الوفيات الإجمالية 95.6%، وبعد 4 أسابيع من الرش لجدران CPLC فقط (أي 82.5%). في مجموعة DDT، كانت نسبة الوفيات أقل من 70% باستمرار لجميع أنواع الجدران في جميع الأوقات بعد الاختبار الحيوي للرش المتبقي داخل المباني. بلغ متوسط معدلات الوفيات التجريبية لـ DDT وSP بعد 12 أسبوعًا من الرش 25.1% و63.2% على التوالي. بالنسبة لأنواع الأسطح الثلاثة، كانت أعلى متوسطات معدلات الوفيات مع DDT هي 61.1% (لـ PMP بعد أسبوعين من الرش المتبقي داخل المباني)، و36.9% (لـ CPLC بعد 4 أسابيع من الرش المتبقي داخل المباني)، و28.9% (لـ CPLC بعد 4 أسابيع من الرش المتبقي داخل المباني). أما أدنى المعدلات فهي 55% (لـ BUU، بعد أسبوعين من الرش المتبقي داخل المباني)، و32.5% (لـ PMP، بعد 4 أسابيع من الرش المتبقي داخل المباني)، و20% (لـ PMP، بعد 4 أسابيع من الرش المتبقي داخل المباني)؛ الرش المتبقي داخل المباني في الولايات المتحدة. بالنسبة لمبيد SP، بلغت أعلى معدلات الوفيات المتوسطة لجميع أنواع الأسطح 97.2% (لمبيد CPLC، بعد أسبوعين من الرش المتبقي)، و82.5% (لمبيد CPLC، بعد 4 أسابيع من الرش المتبقي)، و67.5% (لمبيد CPLC، بعد 4 أسابيع من الرش المتبقي). أما أدنى المعدلات فكانت 94.4% (لمبيد BUU، بعد أسبوعين من الرش المتبقي)، و75% (لمبيد PMP، بعد 4 أسابيع من الرش المتبقي)، و58.3% (لمبيد PMP، بعد 12 أسبوعًا من الرش المتبقي). بالنسبة لكلا المبيدين، تباينت معدلات الوفيات على الأسطح المعالجة بمبيد PMP بشكل أسرع مع مرور الوقت مقارنةً بالأسطح المعالجة بمبيدي CPLC وBUU.
يلخص الجدول 4 آثار التدخل (أي التغيرات في وفرة البعوض بعد الرش الموضعي) لجولات الرش الموضعي القائمة على مادة DDT ومادة SP (الملف الإضافي 1: الشكل S1). بالنسبة للرش الموضعي باستخدام DDT، بلغت نسب الانخفاض في أعداد خنافس الأرجل الفضية بعد فترة الرش 34.1% (بعد أسبوعين)، و25.9% (بعد 4 أسابيع)، و14.1% (بعد 12 أسبوعًا). أما بالنسبة للرش الموضعي باستخدام SP، فكانت نسب الانخفاض 90.5% (بعد أسبوعين)، و66.7% (بعد 4 أسابيع)، و55.6% (بعد 12 أسبوعًا). وسُجلت أكبر انخفاضات في وفرة الروبيان الفضي في المنازل المرجعية خلال فترات الإبلاغ عن الرش الموضعي باستخدام DDT وSP، حيث بلغت 2.8% (بعد أسبوعين) و49.1% (بعد أسبوعين) على التوالي. خلال فترة الرش الموضعي باستخدام مبيد SP، كان انخفاض أعداد طيور التدرج أبيض البطن (قبل وبعد الرش) متقاربًا في المنازل التي خضعت للرش (t² = -9.09، P < 0.001) وفي المنازل التي خضعت للمراقبة (t² = -1.29، P = 0.33). وكانت أعدادها أعلى مقارنةً بالرش الموضعي باستخدام مبيد DDT في جميع الفترات الزمنية الثلاث بعد الرش. بالنسبة لكلا المبيدين، ازداد عدد حشرات البق الفضي في المنازل التي خضعت للمراقبة بعد 12 أسبوعًا من الرش (أي بنسبة 3.6% و9.9% لمبيدي SP وDDT على التوالي). خلال فترتي الرش الموضعي باستخدام مبيدي SP وDDT، تم جمع 112 و161 جمبري فضي من المزارع التي خضعت للمراقبة على التوالي.
لم تُلاحظ فروقٌ ذات دلالة إحصائية في كثافة الروبيان الفضي بين مجموعات الأسر (أي الرش مقابل المراقبة: t(2) = -3.47، P = 0.07؛ الرش مقابل المجموعة الضابطة: t(2) = -2.03، P = 0.18؛ المراقبة مقابل المجموعة الضابطة: خلال أسابيع الرش المتبقي بعد استخدام مادة DDT، t(2) = -0.59، P = 0.62). في المقابل، لوحظت فروقٌ ذات دلالة إحصائية في كثافة الروبيان الفضي بين مجموعة الرش والمجموعة الضابطة (t(2) = -11.28، P = 0.01) وبين مجموعة الرش والمجموعة الضابطة (t(2) = -4.42، P = 0.05). أُجري الرش المتبقي بعد أسابيع قليلة من الرش الوقائي. بالنسبة للرش الوقائي بعد الرش الوقائي، لم تُلاحظ فروقٌ ذات دلالة إحصائية بين الأسر المراقبة والأسر الضابطة (t(2) = -0.48، P = 0.68). يوضح الشكل 2 متوسط كثافة طيور التدرج ذات البطن الفضي في المزارع التي عولجت كليًا وجزئيًا بعجلات الرش الموضعي. لم تُلاحظ فروق ذات دلالة إحصائية في كثافة طيور التدرج المُدارة كليًا بين الأسر التي عولجت كليًا وتلك التي عولجت جزئيًا (بمعدل 7.3 و2.7 طائر لكل مصيدة/ليلة). كما رُشّت بعض الأسر بمبيدي DDT وSP معًا (بمعدل 7.5 و4.4 طائر في الليلة لكل من DDT وSP على التوالي) (t² ≤ 1.0، P > 0.2). مع ذلك، اختلفت كثافة الروبيان الفضي في المزارع التي رُشّت كليًا وجزئيًا اختلافًا ذا دلالة إحصائية بين جولات الرش الموضعي باستخدام SP وDDT (t² ≥ 4.54، P ≤ 0.05).
متوسط الكثافة التقديرية لحشرات البق ذات الأجنحة الفضية في المنازل المعالجة كليًا وجزئيًا في قرية ماهانار، لافابور، خلال الأسبوعين السابقين لجولات الرش المتكامل للمبيدات الحشرية، و2 و4 و12 أسبوعًا بعد جولات الرش المتكامل للمبيدات الحشرية، وDDT، وSP.
تم إعداد خريطة شاملة للمخاطر المكانية (قرية لافابور ماهانار؛ المساحة الإجمالية: 26,723 كم²) لتحديد مناطق المخاطر المكانية المنخفضة والمتوسطة والعالية، وذلك لرصد ظهور وعودة ظهور الروبيان الفضي قبل وبعد عدة أسابيع من تطبيق الرش المتبقي داخل المنازل (الشكلان 3 و4). وقد سُجلت أعلى درجة مخاطر للأسر عند إعداد خريطة المخاطر المكانية بـ "12" (أي "8" لخرائط المخاطر القائمة على مؤشر المخاطر المكانية، و"4" لخرائط المخاطر القائمة على مؤشري VSI وIRSS). أما أدنى درجة مخاطر محسوبة فهي "صفر" أو "لا يوجد خطر"، باستثناء خرائط DDT-VSI وIRSS التي يبلغ حدها الأدنى 1. وأظهرت خريطة المخاطر القائمة على مؤشر المخاطر المكانية أن مساحة كبيرة (أي 19,994.3 كم²؛ 74.8%) من قرية لافابور ماهانار تُعد منطقة عالية المخاطر، حيث يُرجح أن يواجه السكان البعوض ويعاود الظهور. تتفاوت نسبة التغطية بين المناطق عالية الخطورة (DDT 20.2%؛ SP 4.9%)، والمتوسطة الخطورة (DDT 22.3%؛ SP 4.6%)، والمنخفضة/المعدومة الخطورة (DDT 57.5%؛ SP 90.5%) (t(2) = 12.7، P < 0.05) بين مخططات المخاطر لـ DDT وSP-IS وIRSS (الشكلان 3 و4). أظهرت خريطة المخاطر المركبة النهائية أن SP-IRS يتمتع بقدرات وقائية أفضل من DDT-IRS في جميع مستويات مناطق خطر الإصابة بالتسمم الدموي. انخفضت مساحة المنطقة عالية الخطورة للتسمم الدموي إلى أقل من 7% (1837.3 كم²) بعد تطبيق SP-IRS، وأصبحت معظم المساحة (أي 53.6%) منطقة منخفضة الخطورة. خلال فترة الرش الموضعي باستخدام مبيد DDT، بلغت نسبة المناطق عالية الخطورة والمناطق منخفضة الخطورة التي تم تقييمها بواسطة خريطة المخاطر المدمجة 35.5% (9498.1 كم²) و16.2% (4342.4 كم²) على التوالي. تم رسم كثافة ذباب الرمل المقاسة في المنازل المعالجة والمنازل المرجعية قبل وبعد عدة أسابيع من تطبيق الرش الموضعي، وعرضها على خريطة مخاطر مدمجة لكل جولة من جولات الرش (أي DDT وSP) (الشكلان 3 و4). كان هناك توافق جيد بين درجات مخاطر المنازل ومتوسط كثافة الروبيان الفضي المسجلة قبل وبعد الرش الموضعي (الشكل 5). كانت قيم معامل التحديد (R²) (P < 0.05) لتحليل الاتساق المحسوبة من جولتي الرش المتبقي داخل المنازل (IRS) كما يلي: 0.78 قبل أسبوعين من استخدام مبيد DDT، و0.81 بعد أسبوعين من استخدام مبيد DDT، و0.78 بعد أربعة أسابيع من استخدام مبيد DDT، و0.83 بعد 12 أسبوعًا من استخدام مبيد DDT. وبلغ إجمالي استخدام مبيد DDT بعد الرش الوقائي 0.85، و0.82 قبل أسبوعين من الرش الوقائي، و0.38 بعد أسبوعين من الرش الوقائي، و0.56 بعد أربعة أسابيع من الرش الوقائي، و0.81 بعد 12 أسبوعًا من الرش الوقائي، و0.79 بعد أسبوعين من الرش الوقائي إجمالًا (الملف الإضافي 1: الجدول S3). أظهرت النتائج أن تأثير تدخل الرش الوقائي داخل المنازل (SP-IRS) على جميع أنواع الحشرات قد تعزز خلال الأسابيع الأربعة التي تلت الرش المتبقي داخل المنازل. بينما ظل الرش المتبقي داخل المنازل باستخدام مبيد DDT غير فعال لجميع أنواع الحشرات في جميع الأوقات بعد تطبيق الرش المتبقي داخل المنازل. تُلخّص نتائج التقييم الميداني لمنطقة خريطة المخاطر المتكاملة في الجدول 5. بالنسبة لجولات الرش المتبقي داخل المنازل، كان متوسط وفرة روبيان البطن الفضي ونسبة وفرته الإجمالية في المناطق عالية الخطورة (أي > 55%) أعلى منه في المناطق منخفضة ومتوسطة الخطورة في جميع نقاط القياس بعد الرش. تم تحديد مواقع العائلات الحشرية (أي تلك المختارة لجمع البعوض) وتوضيحها في الملف الإضافي 1: الشكل S2.
ثلاثة أنواع من خرائط المخاطر المكانية القائمة على نظم المعلومات الجغرافية (أي HT و IS و IRSS ومزيج من HT و IS و IRSS) لتحديد مناطق خطر حشرة البق النتن قبل وبعد DDT-IRS في قرية ماهنار، لافابور، مقاطعة فايشالي (بيهار).
ثلاثة أنواع من خرائط المخاطر المكانية القائمة على نظم المعلومات الجغرافية (أي HT و IS و IRSS ومزيج من HT و IS و IRSS) لتحديد مناطق خطر الروبيان المرقط الفضي (مقارنة بـ Kharbang)
تم حساب تأثير DDT-(أ، ج، هـ، ز، ط) وSP-IRS (ب، د، و، ح، ي) على مستويات مختلفة من مجموعات المخاطر المنزلية من خلال تقدير معامل التحديد (R²) بين المخاطر المنزلية. تم تقدير مؤشرات المنازل ومتوسط كثافة P. argentipes قبل أسبوعين من تطبيق الرش الموضعي للمبيدات، وبعد أسبوعين، وأربعة أسابيع، واثني عشر أسبوعًا من تطبيقه في قرية لافابور ماهنار، مقاطعة فايشالي، بيهار.
يلخص الجدول 6 نتائج التحليل الأحادي لجميع عوامل الخطر المؤثرة على كثافة رقائق البعوض. وقد تبين أن جميع عوامل الخطر (عددها 6) مرتبطة ارتباطًا وثيقًا بكثافة البعوض في المنازل. ولوحظ أن مستوى الدلالة الإحصائية لجميع المتغيرات ذات الصلة كان أقل من 0.15. وبناءً على ذلك، تم الاحتفاظ بجميع المتغيرات التفسيرية لإجراء تحليل الانحدار المتعدد. وقد تم إنشاء أفضل نموذج نهائي بناءً على خمسة عوامل خطر: TF، وTW، وDS، وISV، وIRSS. ويسرد الجدول 7 تفاصيل المعلمات المختارة في النموذج النهائي، بالإضافة إلى نسب الأرجحية المعدلة، وفترات الثقة 95%، وقيم P. ويُعد النموذج النهائي ذا دلالة إحصائية عالية، حيث بلغت قيمة R² 0.89 (F(5)=27.9، P<0.001).
تم استبعاد متغير TR من النموذج النهائي نظرًا لانخفاض دلالته الإحصائية (P = 0.46) مقارنةً بالمتغيرات التفسيرية الأخرى. استُخدم النموذج المُطوَّر للتنبؤ بكثافة ذباب الرمل استنادًا إلى بيانات من 12 أسرة مختلفة. أظهرت نتائج التحقق وجود ارتباط قوي بين كثافة البعوض المرصودة ميدانيًا وكثافة البعوض المتوقعة بواسطة النموذج (r = 0.91، P < 0.001).
الهدف هو القضاء على داء الليشمانيا الحشوي في الولايات الهندية الموبوءة بحلول عام 2020 [10]. ومنذ عام 2012، أحرزت الهند تقدمًا ملحوظًا في خفض معدلات الإصابة والوفيات الناجمة عن هذا المرض [10]. وشكّل التحوّل من استخدام مادة DDT إلى SP في عام 2015 نقلة نوعية في تاريخ الرش الموضعي للمبيدات في ولاية بيهار الهندية [38]. ولفهم المخاطر المكانية لداء الليشمانيا الحشوي ووفرة نواقله، أُجريت العديد من الدراسات على المستوى الكلي. ومع ذلك، ورغم تزايد الاهتمام بالتوزيع المكاني لانتشار داء الليشمانيا الحشوي في جميع أنحاء البلاد، إلا أن الأبحاث التي أُجريت على المستوى الجزئي قليلة. علاوة على ذلك، تتسم البيانات على المستوى الجزئي بقلة الاتساق وصعوبة التحليل والفهم. وحسب علمنا، تُعد هذه الدراسة أول تقرير يُقيّم الفعالية المتبقية وتأثير التدخل للرش الموضعي للمبيدات باستخدام مبيدي DDT وSP بين حشرات HT ضمن البرنامج الوطني لمكافحة نواقل داء الليشمانيا الحشوي في ولاية بيهار (الهند). وتُعد هذه أيضًا المحاولة الأولى لتطوير خريطة المخاطر المكانية ونموذج تحليل كثافة البعوض للكشف عن التوزيع المكاني والزماني للبعوض على المستوى الميكروسكوبي في ظل ظروف تدخل الرش المتبقي داخل المنازل.
أظهرت نتائجنا أن اعتماد نظام الرش الموضعي بالمبيدات الحشرية (SP-IRS) كان مرتفعًا في جميع الأسر، وأن معظمها قد استكملت جميع مراحل الرش. وأظهرت نتائج الاختبارات البيولوجية أن ذباب الرمل الفضي في القرية محل الدراسة كان شديد الحساسية لمادة بيتا-سيبرمثرين، ولكنه كان أقل حساسية لمادة دي دي تي. ويقل متوسط معدل وفيات الروبيان الفضي الناتج عن دي دي تي عن 50%، مما يشير إلى مستوى عالٍ من المقاومة لهذه المادة. ويتوافق هذا مع نتائج دراسات سابقة أُجريت في أوقات مختلفة في قرى متفرقة من ولايات هندية موبوءة بداء الليشمانيا الحشوي، بما في ذلك ولاية بيهار [8، 9، 39، 40]. إلى جانب حساسية المبيدات، تُعد فعالية المبيدات المتبقية وتأثيرات التدخل معلومات مهمة أيضًا. وتُعد مدة التأثيرات المتبقية مهمة لدورة البرمجة، إذ تُحدد الفترات الزمنية بين جولات الرش الموضعي بالمبيدات، بحيث يبقى السكان محميين حتى الرشة التالية. وكشفت نتائج الاختبارات البيولوجية باستخدام المخاريط عن اختلافات كبيرة في معدل الوفيات بين أنواع أسطح الجدران في نقاط زمنية مختلفة بعد الرش الموضعي بالمبيدات. كان معدل النفوق على الأسطح المعالجة بمادة DDT دائمًا أقل من المستوى المُرضي لمنظمة الصحة العالمية (أي ≥80%)، بينما ظل معدل النفوق على الجدران المعالجة بمادة SP مُرضيًا حتى الأسبوع الرابع بعد الرش الموضعي. من هذه النتائج، يتضح أنه على الرغم من حساسية روبيان الأرجل الفضية الموجود في منطقة الدراسة لمادة SP، فإن فعاليتها المتبقية تختلف باختلاف درجة الحرارة. ومثل مادة DDT، لا تفي مادة SP أيضًا بمدة الفعالية المحددة في إرشادات منظمة الصحة العالمية [41، 42]. قد يعود هذا القصور في الفعالية إلى سوء تطبيق الرش الموضعي (أي تحريك المضخة بالسرعة المناسبة، والمسافة المناسبة من الجدار، ومعدل التدفق، وحجم قطرات الماء، وترسبها على الجدار)، بالإضافة إلى الاستخدام غير الرشيد للمبيدات (أي تحضير المحلول) [11، 28، 43]. ومع ذلك، نظرًا لأن هذه الدراسة أُجريت تحت رقابة صارمة، فقد يكون سبب آخر لعدم الالتزام بتاريخ انتهاء الصلاحية الموصى به من قِبل منظمة الصحة العالمية هو جودة مادة SP (أي نسبة المادة الفعالة) التي تُشكل عنصر مراقبة الجودة.
من بين أنواع الأسطح الثلاثة المستخدمة لتقييم ثبات المبيدات، لوحظت اختلافات كبيرة في معدل النفوق بين سطحي BUU وCPLC لمبيدين. ومن النتائج الجديدة الأخرى أن سطح CPLC أظهر أداءً أفضل في الحفاظ على المبيدات المتبقية في جميع الفترات الزمنية تقريبًا بعد الرش، يليه سطحا BUU وPMP. مع ذلك، بعد أسبوعين من الرش الموضعي، سجل سطح PMP أعلى معدل نفوق لمبيد DDT، وثاني أعلى معدل نفوق لمبيد SP. تشير هذه النتيجة إلى أن المبيد المترسب على سطح PMP لا يدوم طويلًا. قد يُعزى هذا الاختلاف في فعالية بقايا المبيدات بين أنواع الجدران إلى أسباب متعددة، مثل تركيب المواد الكيميائية في الجدار (زيادة الرقم الهيدروجيني تؤدي إلى تحلل بعض المبيدات بسرعة)، ومعدل الامتصاص (أعلى على الجدران الترابية)، وتوافر التحلل البكتيري، ومعدل تحلل مواد الجدار، بالإضافة إلى درجة الحرارة والرطوبة [44، 45، 46، 47، 48، 49]. تدعم نتائجنا العديد من الدراسات الأخرى حول الفعالية المتبقية للأسطح المعالجة بالمبيدات الحشرية ضد مختلف نواقل الأمراض [45، 46، 50، 51].
أظهرت تقديرات انخفاض أعداد البعوض في المنازل المعالجة أن الرش الموضعي باستخدام مبيد SP كان أكثر فعالية من الرش الموضعي باستخدام مبيد DDT في مكافحة البعوض في جميع الفترات الزمنية اللاحقة للرش (P < 0.001). وبالنسبة لجولات الرش الموضعي باستخدام مبيد SP وDDT، تراوحت معدلات الانخفاض في أعداد البعوض في المنازل المعالجة بين 55.6% و90.5% خلال الفترة من أسبوعين إلى 12 أسبوعًا، وبين 14.1% و34.1% على التوالي. كما أظهرت هذه النتائج وجود تأثيرات ملحوظة على وفرة بعوضة P. argentipes في المنازل المرجعية خلال 4 أسابيع من تطبيق الرش الموضعي؛ حيث ازدادت أعدادها في كلتا جولتي الرش بعد 12 أسبوعًا من تطبيقه. ومع ذلك، لم يكن هناك فرق معنوي في أعداد البعوض في المنازل المرجعية بين جولتي الرش الموضعي (P = 0.33). أظهرت نتائج التحليلات الإحصائية لكثافة روبيان الفضة بين مجموعات الأسر في كل جولة عدم وجود فروق ذات دلالة إحصائية في تركيز مادة DDT بين المجموعات الأربع (أي، الرش مقابل المراقبة؛ الرش مقابل عدم الرش؛ المراقبة مقابل عدم الرش؛ الرش الكامل مقابل الرش الجزئي). كما تم رصد فروق ذات دلالة إحصائية في كثافة روبيان الفضة بين جولتي الرش باستخدام DDT والرش باستخدام SP-IRS (أي، المراقبة مقابل عدم الرش، والرش الكامل مقابل الرش الجزئي). ومع ذلك، لوحظت فروق ذات دلالة إحصائية في كثافة روبيان الفضة بين جولتي الرش باستخدام DDT والرش باستخدام SP-IRS في المزارع التي رُشّت جزئيًا وكليًا. تشير هذه الملاحظة، بالإضافة إلى حساب تأثيرات التدخل عدة مرات بعد الرش، إلى أن SP فعال في مكافحة البعوض في المنازل المعالجة جزئيًا أو كليًا، ولكنه غير فعال في المنازل غير المعالجة. ومع ذلك، على الرغم من عدم وجود فروق ذات دلالة إحصائية في عدد البعوض في منازل المراقبة بين جولتي الرش باستخدام DDT والرش باستخدام SP-IRS، إلا أن متوسط عدد البعوض الذي جُمع خلال جولة DDT-IRS كان أقل مقارنةً بجولة SP-IRS. الكمية تتجاوز الكمية. تشير هذه النتيجة إلى أن المبيد الحشري الحساس للنواقل، والذي حقق أعلى تغطية في الرش المتبقي داخل المنازل، قد يكون له تأثير على مكافحة البعوض في المنازل التي لم تُرش. ووفقًا للنتائج، كان لمبيد SP تأثير وقائي أفضل ضد لدغات البعوض من مبيد DDT في الأيام الأولى بعد الرش المتبقي. إضافةً إلى ذلك، ينتمي ألفا-سيبرمثرين إلى مجموعة SP، وله تأثير تهيجي عند التلامس وسمية مباشرة للبعوض، وهو مناسب للرش المتبقي [51، 52]. قد يكون هذا أحد الأسباب الرئيسية لضعف تأثير ألفا-سيبرمثرين في المناطق النائية. وجدت دراسة أخرى [52] أنه على الرغم من أن ألفا-سيبرمثرين أظهر استجابات موجودة ومعدلات قتل عالية في الاختبارات المعملية وفي الأكواخ، إلا أن المركب لم يُنتج استجابة طاردة للبعوض في ظل ظروف معملية مضبوطة.
في هذه الدراسة، تم تطوير ثلاثة أنواع من خرائط المخاطر المكانية؛ حيث تم تقييم تقديرات المخاطر المكانية على مستوى الأسر وعلى مستوى المنطقة من خلال الملاحظات الميدانية لكثافة روبيان الأرجل الفضية. أظهر تحليل مناطق الخطر بناءً على مؤشر المخاطر المكانية أن غالبية مناطق قرية لافابور-ماهانارا (>78%) تقع ضمن أعلى مستوى من خطر ظهور ذبابة الرمل وعودتها. ولعل هذا هو السبب الرئيسي لشعبية ذبابة الرمل في قرية راوالبور ماهانارا. وقد وُجد أن مؤشر المخاطر المكانية الإجمالي ومؤشر المخاطر المكانية للرش المتكامل، بالإضافة إلى خريطة المخاطر المدمجة النهائية، ينتجان نسبة أقل من المناطق المصنفة ضمن المناطق عالية الخطورة خلال جولة الرش المتكامل للروبيان (ولكن ليس خلال جولة الرش المتكامل للروبيان). بعد جولة الرش المتكامل للروبيان، تم تحويل مساحات كبيرة من مناطق المخاطر العالية والمتوسطة بناءً على مؤشر المخاطر المكانية إلى مناطق منخفضة الخطورة (أي 60.5%؛ تقديرات خريطة المخاطر المدمجة)، وهو ما يقل بنحو أربعة أضعاف (16.2%) عن جولة الرش المتكامل للروبيان. - يوضح الرسم البياني أعلاه الوضع في مخطط مخاطر محفظة الرش المتكامل للروبيان. تشير هذه النتيجة إلى أن الرش المتبقي داخل المنازل هو الخيار الصحيح لمكافحة البعوض، ولكن درجة الحماية تعتمد على جودة المبيد الحشري، والحساسية (للناقل المستهدف)، والمقبولية (في وقت الرش المتبقي داخل المنازل) وتطبيقه؛
أظهرت نتائج تقييم المخاطر المنزلية توافقًا جيدًا (P < 0.05) بين تقديرات المخاطر وكثافة روبيان الفضة الذي جُمع من منازل مختلفة. يشير هذا إلى أن معايير المخاطر المنزلية المحددة ودرجات المخاطر الفئوية الخاصة بها مناسبة تمامًا لتقدير وفرة روبيان الفضة محليًا. بلغت قيمة R² لتحليل التوافق بعد الرش الموضعي لمبيد DDT ≥ 0.78، وهي مساوية أو أكبر من القيمة قبل الرش (أي 0.78). أظهرت النتائج أن الرش الموضعي لمبيد DDT كان فعالًا في جميع مناطق المخاطر العالية (أي العالية والمتوسطة والمنخفضة). بالنسبة لجولة الرش الموضعي لمبيد SP، وجدنا أن قيمة R² تذبذبت في الأسبوعين الثاني والرابع بعد تنفيذ الرش الموضعي، وكانت القيم قبل أسبوعين من تنفيذ الرش وبعد 12 أسبوعًا من تنفيذه متقاربة جدًا؛ تعكس هذه النتيجة التأثير الكبير للتعرض للرش الموضعي لمبيد SP على البعوض، والذي أظهر اتجاهًا تنازليًا مع مرور الوقت بعد الرش الموضعي. وقد تم تسليط الضوء على تأثير الرش الموضعي لمبيد SP ومناقشته في فصول سابقة.
أظهرت نتائج التدقيق الميداني لمناطق الخطر في الخريطة المجمعة أنه خلال جولة الرش المتبقي داخل المنازل، جُمعت أعلى أعداد من الروبيان الفضي في المناطق عالية الخطورة (أي أكثر من 55%)، تليها المناطق متوسطة الخطورة ثم منخفضة الخطورة. وباختصار، أثبت التقييم المكاني للمخاطر القائم على نظم المعلومات الجغرافية فعاليته كأداة لاتخاذ القرارات، حيث يجمع طبقات مختلفة من البيانات المكانية، منفردة أو مجتمعة، لتحديد مناطق خطر ذبابة الرمل. وتوفر خريطة المخاطر المطورة فهمًا شاملًا لظروف ما قبل وما بعد التدخل (أي نوع الأسرة، وحالة الرش المتبقي داخل المنازل، وآثار التدخل) في منطقة الدراسة، والتي تتطلب اتخاذ إجراءات فورية أو تحسينات، لا سيما على المستوى الجزئي. وهذا وضع شائع جدًا. في الواقع، استخدمت العديد من الدراسات أدوات نظم المعلومات الجغرافية لرسم خرائط مخاطر مواقع تكاثر النواقل والتوزيع المكاني للأمراض على المستوى الكلي [24، 26، 37].
تم تقييم خصائص المساكن وعوامل الخطر للتدخلات القائمة على الرش المتبقي داخل المنازل إحصائيًا لاستخدامها في تحليلات كثافة الروبيان الفضي. على الرغم من أن جميع العوامل الستة (أي TF، TW، TR، DS، ISV، وIRSS) ارتبطت ارتباطًا وثيقًا بوفرة الروبيان الفضي محليًا في التحليلات أحادية المتغير، إلا أنه تم اختيار عامل واحد فقط من بين خمسة عوامل في نموذج الانحدار المتعدد النهائي. تُظهر النتائج أن خصائص إدارة الأسر وعوامل التدخل للرش المتبقي داخل المنازل (TF، TW، DS، ISV، IRSS، إلخ) في منطقة الدراسة مناسبة لرصد ظهور الروبيان الفضي وتعافيه وتكاثره. في تحليل الانحدار المتعدد، لم يُثبت أن TR ذو دلالة إحصائية، وبالتالي لم يتم اختياره في النموذج النهائي. كان النموذج النهائي ذا دلالة إحصائية عالية، حيث فسرت المعلمات المختارة 89% من كثافة الروبيان الفضي. أظهرت نتائج دقة النموذج وجود ارتباط قوي بين كثافات الروبيان الفضي المتوقعة والملاحظة. تدعم نتائجنا أيضًا الدراسات السابقة التي ناقشت عوامل الخطر الاجتماعية والاقتصادية والإسكانية المرتبطة بانتشار داء الليشمانيا الحشوي والتوزيع المكاني للناقل في المناطق الريفية في بيهار [15، 29].
لم تُقيّم هذه الدراسة ترسب المبيدات على الجدران المرشوشة، ولا جودة المبيدات المستخدمة في الرش الموضعي. إذ يُمكن أن تؤثر الاختلافات في جودة المبيدات وكميتها على معدل نفوق البعوض وفعالية تدخلات الرش الموضعي. وبالتالي، قد يختلف معدل النفوق المُقدّر بين أنواع الأسطح وتأثيرات التدخل بين مجموعات الأسر عن النتائج الفعلية. وبناءً على هذه النقاط، يُمكن التخطيط لدراسة جديدة. يشمل تقييم المساحة الإجمالية المعرضة للخطر (باستخدام خرائط المخاطر الجغرافية) في قرى الدراسة المناطق المفتوحة بينها، مما يؤثر على تصنيف مناطق الخطر (أي تحديد المناطق) ويمتد ليشمل مناطق خطر مختلفة. مع ذلك، أُجريت هذه الدراسة على مستوى جزئي، لذا فإن للأراضي غير المبنية تأثيرًا طفيفًا فقط على تصنيف مناطق الخطر. إضافةً إلى ذلك، يُمكن أن يُتيح تحديد وتقييم مناطق الخطر المختلفة ضمن المساحة الإجمالية للقرية فرصةً لاختيار مناطق لبناء مساكن جديدة في المستقبل (وخاصةً اختيار المناطق منخفضة الخطر). إجمالًا، تُقدّم نتائج هذه الدراسة معلومات متنوعة لم تُدرس من قبل على المستوى الجزئي. والأهم من ذلك، أن التمثيل المكاني لخريطة مخاطر القرية يساعد في تحديد وتجميع الأسر في مناطق المخاطر المختلفة، وبالمقارنة مع المسوحات الأرضية التقليدية، فإن هذه الطريقة بسيطة ومريحة وفعالة من حيث التكلفة وأقل كثافة في العمل، مما يوفر المعلومات لصناع القرار.
تشير نتائجنا إلى أن حشرات السمك الفضي المحلية في القرية محل الدراسة قد طورت مقاومة (أي مقاومة عالية) لمادة DDT، ولوحظ ظهور البعوض مباشرة بعد الرش الموضعي. ويبدو أن ألفا-سيبرمثرين هو الخيار الأمثل لمكافحة نواقل داء الليشمانيا الحشوي بالرش الموضعي نظرًا لفعاليته العالية في القضاء على حشرات السمك الفضي، وقبوله المجتمعي الأفضل مقارنةً بالرش الموضعي باستخدام DDT. ومع ذلك، وجدنا أن معدل نفوق البعوض على الجدران المعالجة بمادة SP يختلف باختلاف نوع السطح، ولوحظت فعالية متبقية ضعيفة، ولم يتم الالتزام بالوقت الموصى به من قبل منظمة الصحة العالمية بعد الرش الموضعي. توفر هذه الدراسة نقطة انطلاق جيدة للنقاش، وتتطلب نتائجها مزيدًا من البحث لتحديد الأسباب الجذرية الحقيقية. وقد أظهرت دقة التنبؤ لنموذج تحليل كثافة ذباب الرمل أنه يمكن استخدام مزيج من خصائص المساكن، وحساسية النواقل للمبيدات الحشرية، وحالة الرش الموضعي لتقدير كثافة ذباب الرمل في القرى الموبوءة بداء الليشمانيا الحشوي في ولاية بيهار. تُظهر دراستنا أيضًا أن رسم خرائط المخاطر المكانية باستخدام نظم المعلومات الجغرافية (على المستوى الكلي) يُمكن أن يكون أداةً فعّالة لتحديد المناطق المعرضة للخطر، ورصد ظهور الكتل الرملية وعودتها قبل وبعد اجتماعات الرش الموضعي. إضافةً إلى ذلك، تُوفر خرائط المخاطر المكانية فهمًا شاملًا لمدى وطبيعة المناطق المعرضة للخطر على مختلف المستويات، وهو ما لا يُمكن دراسته من خلال المسوحات الميدانية التقليدية وأساليب جمع البيانات المعتادة. يُمكن لمعلومات المخاطر المكانية الدقيقة التي يتم جمعها عبر خرائط نظم المعلومات الجغرافية أن تُساعد العلماء وباحثي الصحة العامة على تطوير وتنفيذ استراتيجيات مكافحة جديدة (مثل التدخل الفردي أو المكافحة المتكاملة للنواقل) للوصول إلى فئات مختلفة من الأسر، وذلك بحسب طبيعة مستويات المخاطر. علاوةً على ذلك، تُساعد خريطة المخاطر على تحسين تخصيص موارد المكافحة واستخدامها في الوقت والمكان المناسبين لتعزيز فعالية البرنامج.
منظمة الصحة العالمية. الأمراض المدارية المهملة، نجاحات خفية، وفرص جديدة. 2009. http://apps.who.int/iris/bitstream/10665/69367/1/WHO_CDS_NTD_2006.2_eng.pdf. تاريخ الوصول: 15 مارس/آذار 2014
منظمة الصحة العالمية. مكافحة داء الليشمانيات: تقرير اجتماع لجنة خبراء منظمة الصحة العالمية المعنية بمكافحة داء الليشمانيات. 2010. http://apps.who.int/iris/bitstream/10665/44412/1/WHO_TRS_949_eng.pdf. تاريخ الوصول: 19 مارس 2014
سينغ إس. الاتجاهات المتغيرة في علم الأوبئة والعرض السريري وتشخيص عدوى الليشمانيا وفيروس نقص المناعة البشرية المشتركة في الهند. المجلة الدولية للأمراض المعدية. 2014؛29:103-112.
البرنامج الوطني لمكافحة الأمراض المنقولة بالنواقل (NVBDCP). تسريع برنامج القضاء على داء الليشمانيا الحشوي. 2017. https://www.who.int/leishmaniasis/resources/Accelerated-Plan-Kala-azar1-Feb2017_light.pdf. تاريخ الوصول: 17 أبريل 2018
مونياراج م. مع قلة الأمل في القضاء على داء الليشمانيات الحشوي (الكالازار) بحلول عام 2010، والذي تتكرر نوباته دوريًا في الهند، هل ينبغي إلقاء اللوم على تدابير مكافحة النواقل أو العدوى المصاحبة بفيروس نقص المناعة البشرية أو العلاج؟ توباراسيتول. 2014؛4:10-9.
ثاكور كيه بي: استراتيجية جديدة للقضاء على داء الليشمانيا الحشوي في المناطق الريفية في بيهار. المجلة الهندية للبحوث الطبية. 2007؛ 126: 447-51.
تاريخ النشر: 20 مايو 2024



