Tổng quan nghiên cứu

Bệnh Lở mồm long móng (LMLM - Foot and Mouth Disease) là một trong những bệnh truyền nhiễm nguy hiểm hàng đầu ở động vật móng guốc chẵn, gây thiệt hại kinh tế hàng triệu USD mỗi năm cho ngành chăn nuôi toàn cầu. Tại Việt Nam, dịch bệnh đã tồn tại hơn 100 năm và thường xuyên tái bùng phát theo chu kỳ từ 2 đến 3 năm với diễn biến dịch tễ học vô cùng phức tạp.

Vấn đề nghiên cứu trọng tâm xuất phát từ thực tế giai đoạn 2016 - 2018, khi phạm vi lưu hành của bệnh có xu hướng mở rộng trên nhiều đối tượng gia súc, đặc biệt là sự gia tăng đột biến về tỷ lệ mắc và tử vong trên đàn lợn thịt cùng đàn bò vận chuyển qua biên giới. Mục tiêu cụ thể của luận văn nhằm xác định mối quan hệ và ảnh hưởng của các yếu tố không gian, thời gian đến tần suất phân bố ổ dịch; định danh chính xác các serotype, topotype vi rút lưu hành; đồng thời đánh giá mức độ tương đồng kháng nguyên của vắc xin và đáp ứng miễn dịch bảo hộ trên đàn gia súc thực địa.

Phạm vi nghiên cứu bao quát toàn bộ 63 tỉnh, thành phố trên cả nước trong giai đoạn 3 năm từ 2016 đến 2018, kết hợp giám sát huyết thanh học sau tiêm phòng vắc xin tại các vùng trọng điểm trong năm 2018. Ý nghĩa thực tiễn của công trình được khẳng định qua việc xác lập chỉ số nguy cơ dịch bệnh trung bình đạt 5,95 xã có dịch/100 xã-năm (95% CI 5,51 - 6,39), đồng thời chứng minh tỷ lệ bảo hộ kháng thể sau tiêm phòng đạt từ 74,0% đến 89,3%. Kết quả này cung cấp cơ sở khoa học trực tiếp phục vụ việc điều chỉnh "Chương trình quốc gia phòng, chống bệnh LMLM giai đoạn 2016 - 2020" do Bộ Nông nghiệp và Phát triển nông thôn ban hành.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng thuyết tam giác dịch tễ học truyền nhiễm kinh điển (tương tác mầm bệnh - vật chủ - môi trường) và lý thuyết bệnh học phân tử của vi rút học thú y.

Tác nhân gây bệnh LMLM thuộc giống Aphthovirus, họ Picornaviridae, mang cấu trúc hình khối đa diện 20 mặt đều với kích thước siêu vi từ 25 đến 30 nm. Bộ gen của vi rút là chuỗi đơn dương ARN chứa khoảng 8450 bazơ, mã hóa cho 4 protein cấu trúc gồm VP1, VP2, VP3, VP4 và 10 protein phi cấu trúc. Trong đó, protein bề mặt VP1 giữ vai trò quyết định trong việc gắn kết thụ thể tế bào và kích thích tạo kháng thể trung hòa. Vi rút có sức đề kháng cao ở nhiệt độ lạnh, sống tới 425 ngày ở 0°C nhưng bị bất hoạt sau 5 đến 15 phút ở nhiệt độ 60 - 70°C hoặc môi trường có độ pH dưới 6,5 và trên 11.

Các khái niệm trọng tâm trong nghiên cứu bao gồm:

  • Hệ thống 7 serotype: Phân định các nhóm huyết thanh O, A, C, SAT-1, SAT-2, SAT-3 và Asia1 không có miễn dịch chéo.
  • Tình trạng mang trùng (Carrier state): Hiện tượng vi rút tồn tại dai dẳng tại hạch hầu họng trên 28 ngày sau nhiễm, kéo dài tới 3,5 năm ở loài bò.
  • Hệ số tương quan kháng nguyên (r1 value): Chỉ số đo lường mức độ tương đồng giữa chủng vắc xin và chủng vi rút thực địa, trong đó giá trị r1 ≥ 0,3 đại diện cho khả năng bảo hộ chéo hiệu quả.

Phương pháp nghiên cứu

  • Nguồn dữ liệu và cỡ mẫu: Thu thập toàn bộ dữ liệu dịch tễ học báo cáo từ cấp xã của 63 tỉnh thành giai đoạn 2016 - 2018 do Cục Thú y quản lý. Cỡ mẫu giám sát huyết thanh năm 2018 được thu thập tại các huyện trọng điểm: mỗi huyện chọn ngẫu nhiên 3 xã, mỗi xã lấy 10 mẫu huyết thanh trâu và 10 mẫu huyết thanh bò (tổng số 60 mẫu/huyện), thực hiện lấy mẫu sau mũi tiêm phòng cuối cùng từ 4 tuần trở lên.
  • Phương pháp chọn mẫu: Sử dụng phương pháp chọn mẫu phân tầng kết hợp ngẫu nhiên theo vùng nguy cơ dịch tễ (vùng khống chế và vùng đệm) nhằm đảm bảo tính đại diện cao nhất cho quần thể gia súc cảm nhiễm.
  • Phương pháp phân tích: Áp dụng thiết kế nghiên cứu thuần tập hồi cứu (Retrospective cohort study) để xác định tỷ lệ mới mắc và nguy cơ dịch tễ. Phân tích không gian sử dụng thuật toán ước tính mật độ Kernel (Kernel density estimation) với kích thước ô lưới 200m x 200m trên phần mềm ArcGIS 10.3 và ngôn ngữ R (gói epiR, Spatstat). Diễn biến thời gian được theo dõi qua Tỷ số lây lan ước tính (Estimated Dissemination Ratio - EDR) theo chu kỳ 5 ngày. Phân tích phân tử thực hiện qua kỹ thuật giải trình tự đoạn gen VP1 (kích thước 1300 bp) tại Phòng thí nghiệm tham chiếu OIE Pirbright (Vương quốc Anh) và kiểm tra hiệu giá kháng thể bảo hộ bằng phản ứng LPB-ELISA theo Tiêu chuẩn Quốc gia TCVN 8400-1:2010.

Lý do lựa chọn tổ hợp phương pháp này là nhằm tích hợp đồng bộ giữa dịch tễ học không gian, mô hình hóa toán học và sinh học phân tử, khắc phục triệt để nhược điểm của các nghiên cứu mô tả cắt ngang truyền thống.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

  • Cơ cấu loài mắc bệnh và mức độ cảm nhiễm: Trong giai đoạn 2016 - 2018, bò là loài có tỷ lệ mắc bệnh cao nhất, chiếm 45,9% (95% CI 45,3 - 46,5) tổng số ca bệnh; tiếp theo là lợn chiếm 42,5% (95% CI 42,0 - 43,1); trâu chiếm 11,3% (95% CI 11,0 - 11,7) và thấp nhất là dê, cừu với 0,3% (95% CI 0,2 - 0,3). Đáng chú ý, năm 2017 tỷ lệ bò mắc bệnh lên tới 75,7%, trong khi năm 2018 dịch chuyển biến mạnh sang đàn lợn với tỷ lệ mắc đạt 58,9%.
  • Tỷ lệ tử vong cục bộ: Tổng số gia súc chết do dịch trong 3 năm là 543 con. Lợn là loài chịu tổn thất lớn nhất, chiếm 67,4% (366 con chết, 95% CI 62,8 - 72,0) tổng số gia súc tử vong. Tỷ lệ chết ở bò chiếm 19,9% (108 con), dê cừu chiếm 11,0% (60 con) và trâu chỉ chiếm 1,7% (9 con).
  • Phân bố không gian và động lực thời gian: Nguy cơ trung bình các xã xuất hiện bệnh LMLM là 5,95 xã có dịch/100 xã-năm. Dịch bệnh tập trung mật độ cao tại các tỉnh miền núi phía Bắc, Tây Nguyên và Duyên hải Nam Trung Bộ. Về mặt thời gian, dịch xảy ra quanh năm nhưng bùng phát dữ dội vào 2 tháng cuối năm 2018 với số điểm dịch ghi nhận cao nhất lên tới 15 xã/ngày, đường cong chỉ số EDR vượt ngưỡng 1 báo hiệu sự phát tán diện rộng.
  • Đặc điểm phân tử và đáp ứng miễn dịch: Xác định sự lưu hành của serotype O (thuộc các topotype Mya-98, Cathay, PanAsia, Ind-2001d, Ind-2001e) và serotype A (lineage Sea-97). Đáng lưu ý, các chủng O/SEA/Mya-98 và O/Cathay phân lập trên lợn có độ tương đồng kháng nguyên thấp với chủng vắc xin O1Manisa và O3039 (r1 < 0,3). Ngược lại, vắc xin chứa kháng nguyên A22 Iraq đạt tương đồng cao với vi rút A/ASIA/Sea-97 (r1 > 0,3). Tỷ lệ trâu bò có kháng thể bảo hộ sau tiêm phòng vắc xin serotype O đạt trên 70% (bò sữa đạt đỉnh 89,3%), serotype A đạt 74,0% (95% CI 70,6 - 77,2).

Thảo luận kết quả

Sự gia tăng tỷ lệ mắc ở bò (45,9% so với mức 21,1% của giai đoạn 2006 - 2012) bắt nguồn từ hoạt động vận chuyển gia súc thương phẩm qua đường tiểu ngạch biên giới và sự thiếu hụt miễn dịch bảo hộ tại các ổ dịch cũ. Tỷ lệ tử vong cao đột biến ở lợn (67,4%) phản ánh sự lưu hành của topotype Cathay và Mya-98 thích ứng độc lực cao trên loài lợn, kết hợp với thực tế đàn lợn thịt thương phẩm nuôi nhỏ lẻ hầu như không được tiêm phòng định kỳ.

Về phương diện trực quan hóa, dữ liệu dịch tễ thời gian được mô hình hóa hiệu quả qua biểu đồ đường cong EDR kết hợp biểu đồ cột thể hiện tần suất ca bệnh theo ngày. Không gian nguy cơ dịch bệnh được thể hiện chi tiết trên bản đồ mật độ Kernel, phân tầng các vùng nguy cơ thành các mảng màu nhiệt trực quan tương ứng với mật độ điểm dịch trên từng ô chuẩn hóa 200m x 200m, hỗ trợ đắc lực cho việc khoanh vùng kiểm soát.

Đề xuất và khuyến nghị

  1. Chuẩn hóa và cập nhật chủng kháng nguyên vắc xin: Cục Thú y cùng các cơ quan nghiên cứu chuyên ngành cần rà soát, thay thế hoặc bổ sung các chủng vắc xin serotype O tương thích với topotype Mya-98 và Cathay thực địa (đảm bảo r1 ≥ 0,3); tiếp tục duy trì sử dụng vắc xin serotype A chứa kháng nguyên A22 Iraq. Mục tiêu nâng tỷ lệ đàn gia súc có bảo hộ kháng thể đạt trên 80% tại 100% các huyện thuộc vùng đệm và vùng khống chế trước quý 4 hàng năm.
  2. Mở rộng diện tiêm phòng bắt buộc cho đàn lợn thịt: Ủy ban nhân dân và Chi cục Chăn nuôi và Thú y các tỉnh cần điều chỉnh chính sách, đưa đàn lợn thịt tại các nông hộ nhỏ lẻ và vùng có mật độ nuôi cao vào diện tiêm phòng vắc xin bắt buộc. Mục tiêu giảm tỷ lệ tử vong do LMLM trên đàn lợn xuống dưới 10% trong vòng 12 đến 24 tháng tới.
  3. Ứng dụng hệ thống giám sát và dự báo sớm bằng chỉ số EDR: Trung tâm Chẩn đoán Thú y Trung ương cần thiết lập nền tảng số hóa theo dõi chỉ số EDR theo chu kỳ 5 ngày trên phạm vi toàn quốc. Khi chỉ số EDR vượt ngưỡng 1 trong 2 chu kỳ liên tiếp, cơ quan thú y địa phương phải kích hoạt ngay quy trình ứng phó khẩn cấp và khoanh vùng ổ dịch trước thời điểm lây lan diện rộng từ 10 đến 15 ngày.
  4. Tăng cường kiểm dịch vận chuyển liên tỉnh và biên giới: Lực lượng Thú y phối hợp liên ngành (Quản lý thị trường, Công an giao thông) siết chặt kiểm soát 100% các phương tiện vận chuyển trâu, bò, lợn qua các tuyến quốc lộ huyết mạch tại miền Bắc, miền Trung và Tây Nguyên, đặc biệt trong khoảng thời gian cao điểm từ tháng 10 đến tháng 1 năm sau; kiên quyết xử lý và tiêu hủy động vật nhập lậu không rõ nguồn gốc kiểm dịch.

Đối tượng nên tham khảo luận văn

  • Cơ quan quản lý nhà nước và hoạch định chính sách Thú y: Cán bộ Cục Thú y, Sở Nông nghiệp và Phát triển nông thôn có thể khai thác các chỉ số nguy cơ dịch tễ học không gian (5,95 xã/100 xã-năm) và bản đồ phân vùng Kernel để xây dựng quy hoạch chăn nuôi an toàn dịch bệnh, phân bổ ngân sách dự phòng vắc xin hàng năm.
  • Nhà khoa học, giảng viên và học viên cao học ngành Thú y: Cung cấp tài liệu tham khảo chuẩn mực về quy trình ứng dụng dịch tễ học hồi cứu, phân tích chuỗi gen VP1 kích thước 1300 bp và phương pháp đánh giá tương quan kháng nguyên r1 phục vụ nghiên cứu bệnh truyền nhiễm.
  • Doanh nghiệp sản xuất, nhập khẩu vắc xin thú y: Sử dụng kết quả kiểm tra tương đồng vắc xin (vaccine matching) để định hướng nhập khẩu, đăng ký lưu hành và sản xuất các dòng vắc xin đơn giá, nhị giá O-A đạt hiệu lực bảo hộ chéo cao với các biến chủng đang lưu hành tại Việt Nam.
  • Chủ trang trại và kỹ sư quản lý chăn nuôi quy mô công nghiệp: Tham khảo thông tin về độc lực vi rút trên từng đối tượng vật nuôi để tối ưu hóa quy trình an toàn sinh học, thiết lập lịch tiêm phòng bổ sung cho đàn lợn và trâu bò nhằm duy trì mức kháng thể bảo hộ đàn trên 70%.

Câu hỏi thường gặp

Tại sao bò có tỷ lệ mắc bệnh cao nhất nhưng lợn lại có tỷ lệ tử vong cao nhất? Bò có tỷ lệ mắc cao nhất (45,9%) do niêm mạc hô hấp rất nhạy cảm (chỉ cần liều nhiễm 20 TCID50) và chịu tác động lớn từ việc lưu thông buôn bán liên tỉnh. Ngược lại, lợn chiếm tới 67,4% tổng số gia súc chết (366 con) do vi rút gây thoái hóa hoại tử cơ tim cấp tính (hiện tượng tim rằn ri) và thực tế đàn lợn thịt không được tiêm phòng vắc xin đầy đủ.

Chỉ số tỷ số lây lan ước tính EDR có ý nghĩa như thế nào trong kiểm soát dịch LMLM? Chỉ số EDR đo lường tỷ số giữa số ổ dịch trong chu kỳ 5 ngày ủ bệnh với 5 ngày liền kề trước đó. Khi EDR vượt ngưỡng 1 (như ghi nhận vào tháng 12/2018 với 15 xã có dịch/ngày), dịch đang trong pha phát tán mạnh, là căn cứ khoa học cảnh báo cơ quan chức năng phong tỏa ổ dịch tức thì.

Hiện tượng sai khác kháng nguyên giữa chủng thực địa và vắc xin LMLM thể hiện ra sao? Kết quả phân tích vaccine matching ghi nhận các chủng O/SEA/Mya-98 và O/Cathay lưu hành trên lợn có giá trị r1 < 0,3 khi so sánh với kháng nguyên vắc xin O1Manisa và O3039. Điều này khẳng định vắc xin chứa hai chủng kháng nguyên trên không đủ khả năng bảo hộ chéo cho đàn lợn trước các chủng thực địa này.

Hiệu lực bảo hộ kháng thể sau tiêm phòng vắc xin LMLM tại Việt Nam có đạt chuẩn quốc tế? Có. Đánh giá bằng kỹ thuật LPB-ELISA cho thấy tỷ lệ trâu bò dương tính với kháng thể serotype O đạt trên 70% (bò sữa đạt 89,3%) và serotype A đạt 74,0% (95% CI 70,6 - 77,2). Các chỉ số này đều vượt ngưỡng yêu cầu bảo hộ tối thiểu 70% theo quy chuẩn của Tổ chức Thú y Thế giới (OIE).

Mô hình phân tích không gian Kernel hỗ trợ công tác thú y như thế nào? Thuật toán Kernel chia bản đồ thực địa thành các ô lưới 200m x 200m để tính toán mật độ ổ dịch lặp lại. Phương pháp này giúp nhận diện chính xác các điểm nóng dịch tễ tại miền Bắc, Tây Nguyên và Nam Trung Bộ, loại bỏ các sai số không gian so với phương pháp đếm ca bệnh hành chính thông thường.

Kết luận

  • Xác định nguy cơ dịch bệnh LMLM trung bình cả nước giai đoạn 2016 - 2018 là 5,95 xã có dịch/100 xã-năm, tập trung mật độ cao tại miền Bắc, Tây Nguyên và Duyên hải Nam Trung Bộ.
  • Chứng minh tính chu kỳ thời gian với đỉnh dịch bùng phát nghiêm trọng vào 2 tháng cuối năm 2018, đạt đỉnh 15 xã có dịch/ngày khi chỉ số EDR vượt ngưỡng 1.
  • Làm rõ sự phân hóa vật chủ khi bò dẫn đầu về tỷ lệ mắc bệnh (45,9%) nhưng lợn chịu tổn thất nặng nề nhất về tỷ lệ tử vong (67,4%).
  • Định danh các chủng vi rút lưu hành gồm serotype O (Mya-98, Cathay, PanAsia, Ind-2001d,e) và serotype A (Sea-97); phát hiện sự sai lệch kháng nguyên (r1 < 0,3) giữa chủng vắc xin serotype O hiện hành với chủng thực địa trên lợn.
  • Khẳng định hiệu quả bảo hộ miễn dịch sau tiêm phòng đạt chuẩn OIE với tỷ lệ kháng thể serotype O đạt trên 70% và serotype A đạt 74,0%.

Đóng góp chính của luận văn là cung cấp bức tranh dịch tễ học không gian - thời gian và sinh học phân tử toàn diện, làm cơ sở khoa học cốt lõi cho các chiến lược thú y quốc gia. Bước tiếp theo cần tập trung mở rộng ngân hàng dữ liệu gen vi rút và hoàn thiện quy trình tiêm phòng bao phủ cho đàn lợn thịt giai đoạn 2021 - 2025.

Quý độc giả, nhà quản lý và các chuyên gia chăn nuôi hãy tham khảo toàn văn luận văn thạc sĩ để ứng dụng các mô hình dịch tễ học chuyên sâu vào thực tiễn kiểm soát an toàn dịch bệnh.