Giới thiệu dự án

Bệnh lao (Tuberculosis - TB) do trực khuẩn Mycobacterium tuberculosis gây ra tiếp tục là gánh nặng y tế công cộng nghiêm trọng trên toàn cầu. Theo báo cáo của Tổ chức Y tế Thế giới (WHO), năm 2017 toàn cầu ghi nhận khoảng 10 triệu ca mắc mới và hơn 1,3 triệu ca tử vong. Việt Nam xếp thứ 16/30 quốc gia có gánh nặng bệnh lao cao nhất thế giới (124.000 ca mắc mới, 12.840 ca tử vong năm 2017) và xếp thứ 13/30 về gánh nặng lao đa kháng thuốc (MDR-TB). Trong phác đồ điều trị lao chuẩn hàng 1 của Chương trình Chống lao Quốc gia, Isoniazid (INH) đóng vai trò then chốt nhờ khả năng ức chế sinh tổng hợp acid mycolic ở thành tế bào vi khuẩn.

Tuy nhiên, hiệu quả lâm sàng và độc tính của INH có sự dao động lớn giữa các cá thể. Enzym Arylamine N-Acetyltransferase 2 (NAT2) ở gan chịu trách nhiệm chuyển hóa 50% - 90% lượng INH hấp thu thành acetyl-isoniazid không độc. Các biến thể di truyền đơn nucleotide (Single Nucleotide Polymorphisms - SNP) trên gen NAT2 làm suy giảm hoạt tính hoặc độ bền enzym, dẫn đến tích lũy nồng độ INH tự do và gia tăng chuyển hóa qua con đường phụ thuộc Cytochrome P450 2E1 (CYP2E1). Quá trình này tạo ra các chất trung gian độc hại như hydrazine và acetylhydrazine, gây tổn thương gan do thuốc (Drug-Induced Liver Injury - DILI) ở 20% bệnh nhân và viêm dây thần kinh ngoại biên do cạn kiệt pyridoxine (vitamin B6). Ngược lại, nhóm chuyển hóa nhanh có nguy cơ thất bại điều trị sớm lên tới 15% - 38% khi áp dụng liều chuẩn (Standard Treatment Dose - STD: 5 mg/kg/ngày).

          +--------------------------------------------------+
          |         Isoniazid (INH) Uống (5 mg/kg)           |
          +--------------------------------------------------+
                                   |
                  +----------------+----------------+
                  | (50% - 90%)                     | (Con đường phụ)
                  v                                 v
        [ Enzym NAT2 (Gan) ]               [ Thủy phân trực tiếp ]
                  |                                 |
                  +------------+                    v
                  |            |               Hydrazine
                  v            v                    |
          Acetyl-isoniazid   (NAT2 yếu)             v (CYP2E1 oxi hóa)
                  |            |             +-----------------------+
                  v            +------------>| Gốc tự do độc với gan |
           Thải trừ an toàn                  |       (INH-DILI)      |
             qua nước tiểu                   +-----------------------+

Đề tài tập trung giải quyết các mục tiêu cụ thể:

  1. Xây dựng và chuẩn hóa quy trình phân tích các đa hình di truyền trên gen NAT2 sử dụng kỹ thuật PCR-RFLP và giải trình tự Sanger.
  2. Ứng dụng quy trình đã tối ưu hóa để khảo sát tần số alen, kiểu gen và dự đoán kiểu hình acetyl hóa trên 100 bệnh nhân lao phổi người Việt Nam.
  3. Đánh giá tính khả thi kỹ thuật và đối chiếu mức độ tương đồng giữa phương pháp RFLP (3 SNP) và giải trình tự kết hợp thuật toán NAT2PRED (6 SNP).

Phạm vi nghiên cứu thực hiện trên 100 mẫu máu toàn phần của bệnh nhân lao phổi có AFB dương tính tại Bệnh viện Phổi Trung ương, Bệnh viện 74 Trung ương và Bệnh viện Phổi Hà Nội trong giai đoạn từ tháng 6/2017 đến tháng 9/2018.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Các kỹ thuật phân tích đột biến gen NAT2 phổ biến hiện nay bao gồm Real-time PCR, Sắc ký lỏng cao áp biến tính (DHPLC), Phân tích đa hình chiều dài đoạn cắt giới hạn (PCR-RFLP) và Giải trình tự Sanger (Direct Sequencing).

Phương pháp Ưu điểm Nhược điểm Chi phí / mẫu Độ phức tạp
PCR-RFLP Chi phí rẻ, thiết bị cơ bản, thao tác nhanh, kết quả tin cậy cho các SNP đã biết Cần gel nồng độ cao cho đoạn ngắn, không phát hiện được đột biến mới Thấp (~50.000 - 80.000 VNĐ) Trung bình
Sanger Sequencing Tiêu chuẩn vàng, bao phủ toàn bộ 870 bp exon 2, xác định đồng thời 6+ SNP Chi phí cao hơn, thời gian chờ gửi mẫu dài Cao (~150.000 - 250.000 VNĐ) Cao
Real-time PCR Tốc độ nhanh, khép kín, tránh ngoại nhiễm Chi phí đầu tư máy và probe TaqMan rất đắt Rất cao Trung bình
DHPLC Tự động hóa cao, nhạy với dị sợi kép Thiết bị chuyên biệt đắt tiền, chi phí bảo trì lớn Cao Rất cao

Yêu cầu kỹ thuật được phân loại theo mô hình MoSCoW:

  • Must-have (Bắt buộc): Chiết xuất DNA đạt độ tinh sạch OD260/OD280 từ 1,6 - 2,0; khuếch đại thành công đoạn gen NAT2 kích thước 1093 bp; xác định chính xác 3 SNP cốt lõi (NAT2*5, *6, *7).
  • Should-have (Cần có): Đồng kiểm tra 6 vị trí SNP (NAT2*5, *6, *7, *11, *12, *13) qua giải trình tự Sanger để dự đoán kiểu hình bằng NAT2PRED.
  • Could-have (Có thể mở rộng): Tích hợp phần mềm tự động hóa tính liều INH cá thể hóa dựa trên kiểu gen.
  • Won't-have (Chưa triển khai): Giải trình tự toàn bộ hệ gen (WGS) hoặc giải trình tự thế hệ mới (NGS) do không tối ưu về chi phí thường quy.

Thiết kế hệ thống

Quy trình kỹ thuật được thiết kế khép kín gồm 4 khối chính: Khối xử lý mẫu sinh phẩm, Khối khuếch đại phân tử (PCR), Khối phân tích kiểu gen (RFLP/Sanger) và Khối xử lý tin sinh học.

+-----------------------------------------------------------------------------------+
|                           KIẾN TRÚC QUY TRÌNH PHÂN TÍCH                          |
+-----------------------------------------------------------------------------------+
  [ Mẫu máu EDTA 2 mL ] 
  [ Nhánh 1: PCR-RFLP ]           [ Nhánh 2: Sanger Sequencing ]
  • KpnI   (c.481C>T / NAT2*5)    • Gửi First BASE (Malaysia)
  • TaqI   (c.590G>A / NAT2*6)    • BioEdit v7.0.9 đọc chromatogram
  • BamHI  (c.857G>A / NAT2*7)    • NAT2PRED phân loại 6-SNP
  [ Phân loại kiểu hình Acetyl hóa: Nhanh (25%) | TB (38%) | Chậm (37%) ]
  [ Khuyến nghị lâm sàng: Điều chỉnh liều INH theo Dược di truyền học (PGx) ]
  • Công nghệ và Hóa chất sử dụng:
    • Kit tách chiết: E.Z.N.A. Blood DNA Mini Kit (Omega Bio-tek).
    • Enzym PCR: KAPA2G™ Robust HotStart ReadyMix (2X) (Kapabiosystems).
    • Enzym cắt giới hạn: FastDigest KpnI, FastDigest TaqI, FastDigest BamHI (Thermo Scientific).
    • Phần mềm xử lý: BioEdit Sequence Alignment Editor version 7.0.9 / 7.2.5; Web server NAT2PRED; Chi-square calculator.
  • Tiêu chuẩn an toàn sinh học & Bảo mật: Mẫu bệnh phẩm được mã hóa danh tính (anonymized ID), thu thập và phân tích trong tủ an toàn sinh học cấp 2 (ESCO Class II), đạt chuẩn đạo đức nghiên cứu y sinh học.

Methodology

Nghiên cứu ứng dụng phương pháp thực nghiệm tiến cứu kết hợp phân tích di truyền quần thể:

  • Milestone 1: Thu nhận 100 mẫu máu toàn phần từ bệnh nhân lao, bảo quản ở -20°C.
  • Milestone 2: Tối ưu hóa điều kiện luân nhiệt và nồng độ hóa chất phản ứng PCR.
  • Milestone 3: Thực hiện phản ứng cắt giới hạn RFLP và giải trình tự Sanger 100% mẫu.
  • Milestone 4: Kiểm định cân bằng di truyền Hardy-Weinberg ($\chi^2$-test) và thống kê tỷ lệ kiểu hình.

Implementation và kết quả

Development process

Cặp mồi đặc hiệu được thiết kế để khuếch đại toàn bộ khung đọc mở exon 2 dài 870 bp của gen NAT2 (NCBI Accession ID: NM_000015.2), tạo sản phẩm PCR dài 1093 bp:

  • Mồi xuôi (Forward): 5'-GGA ACA AAT TGG ACT TGG-3'
  • Mồi ngược (Reverse): 5'-TCT AGC ATG AAT CAC TCT GC-3'

Quy trình nhiệt PCR tối ưu được thiết lập trên máy Prime Thermal Cycler:

+-----------------------------------------------------------------------------+
|               CHU TRÌNH NHIỆT PHẢN ỨNG PCR GEN NAT2 (1093 bp)               |
+-----------------------------------------------------------------------------+
1. Biến tính ban đầu : 95°C trong 3 phút
2. Chu kỳ lặp lại (35 chu kỳ):
3. Kéo dài cuối cùng : 72°C trong 10 phút
4. Bảo quản          : 4°C vô hạn
+-----------------------------------------------------------------------------+

Thành phần phản ứng 10 µL gồm: 5,0 µL MasterMix Kapa2G Robust (2X), 0,3 µL mồi xuôi (10 µM), 0,3 µL mồi ngược (10 µM), 1,0 µL DNA khuôn (200 ng/µL) và 3,4 µL nước khử ion siêu sạch.

Quy trình cắt giới hạn FastDigest RFLP: Thể tích 10 µL gồm 3,5 µL DNA sản phẩm PCR (100 ng/µL), 1,0 µL đệm 10X FastDigest Buffer, 1,0 µL enzyme FastDigest tương ứng và 4,5 µL nước khử ion; ủ 5 phút ở 37°C, sau đó bất hoạt enzyme ở 80°C trong 5 phút.

+-----------------------------------------------------------------------------+
|                      SƠ ĐỒ PHÂN ĐOẠN CẮT GIỚI HẠN PCR-RFLP                  |
+-----------------------------------------------------------------------------+
• NAT2*5 (FastDigest KpnI nhận biết c.481C>T):
  - Đồng hợp tử dại (CC)     : Cắt thành 2 băng  --> [ 659 bp ] + [ 434 bp ]
  - Dị hợp tử (CT)           : 3 băng            --> [ 1093 bp ] + [ 659 bp ] + [ 434 bp ]
  - Đồng hợp tử đột biến (TT): Không bị cắt      --> [ 1093 bp ]

• NAT2*6 (FastDigest TaqI nhận biết c.590G>A):
  - Đồng hợp tử dại (GG)     : 4 băng            --> [ 381 bp ] + [ 316 bp ] + [ 226 bp ] + [ 170 bp ]
  - Dị hợp tử (GA)           : 5 băng            --> [ 396 bp ] + [ 381 bp ] + [ 316 bp ] + [ 226 bp ] + [ 170 bp ]
  - Đồng hợp tử đột biến (AA): 3 băng            --> [ 396 bp ] + [ 381 bp ] + [ 316 bp ]

• NAT2*7 (FastDigest BamHI nhận biết c.857G>A):
  - Đồng hợp tử dại (GG)     : Cắt thành 2 băng  --> [ 810 bp ] + [ 283 bp ]
  - Dị hợp tử (GA)           : 3 băng            --> [ 1093 bp ] + [ 810 bp ] + [ 283 bp ]
  - Đồng hợp tử đột biến (AA): Không bị cắt      --> [ 1093 bp ]
+-----------------------------------------------------------------------------+

Testing và validation

  1. Hiệu suất tách chiết: 100/100 mẫu đạt nồng độ từ 20 đến 230 ng/µL; tỷ lệ quang phổ OD260/OD280 đạt từ 1,6 đến 2,2; điện di trên gel agarose 1,2% cho băng DNA tổng số nguyên vẹn tại vị trí 23,1 kb (đối chứng thang $\lambda$/HindIII).
  2. Độ đặc hiệu PCR: Tại nhiệt độ gắn mồi $T_a = 58^\circ\text{C}$ và nồng độ DNA 200 ng/µL, phản ứng đạt độ đặc hiệu 100%, vạch sáng rõ tại kích thước 1093 bp, hoàn toàn không xuất hiện băng phụ hay dimer mồi; chứng âm không xuất hiện tạp nhiễm.
  3. Phân giải đoạn cắt: Điện di trên gel agarose 2,0% (cho KpnI, BamHI) và gel polyacrylamide 10% nhuộm bạc (cho TaqI) cho phép phân tách rõ các đoạn chênh lệch từ 15 bp.

Kết quả đạt được

Phân tích đa hình trên 100 bệnh nhân lao cho thấy sự hiện diện đầy đủ của các biến thể với phân bố tuân theo định luật Hardy-Weinberg ($p > 0,05$):

Đa hình SNP Vị trí nucleotide Thay thế acid amin Đồng hợp tử dại Dị hợp tử Đồng hợp tử đột biến Tần số alen Giá trị $\chi^2$ ($p$-value)
NAT2*5 (rs1801280) 341T>C Ile114Thr TT = 89 TC = 10 CC = 1 T: 0,94 | C: 0,06 1,288 ($p = 0,256$)
NAT2*6 (rs1799930) 590G>A Arg197Gln GG = 50 GA = 37 AA = 13 G: 0,685 | A: 0,32 2,034 ($p = 0,154$)
NAT2*7 (rs1799931) 857G>A Gly286Glu GG = 78 GA = 20 AA = 2 G: 0,88 | A: 0,12 0,281 ($p = 0,596$)
NAT2*11 (rs1799929) 481C>T Leu161Leu CC = 89 CT = 10 TT = 1 C: 0,94 | T: 0,06 1,288 ($p = 0,256$)
NAT2*12 (rs1208) 803A>G Arg268Lys AA = 88 AG = 11 GG = 1 A: 0,93 | G: 0,07 0,903 ($p = 0,342$)
NAT2*13 (rs1041983) 282C>T Tyr94Tyr CC = 29 CT = 44 TT = 27 C: 0,51 | T: 0,49 1,432 ($p = 0,232$)

Phân loại kiểu hình acetyl hóa trên 100 bệnh nhân lao Việt Nam ghi nhận:

  • Acetyl hóa nhanh (Rapid): 25 bệnh nhân (25,0%).
  • Acetyl hóa trung bình (Intermediate): 38 bệnh nhân (38,0%).
  • Acetyl hóa chậm (Slow): 37 bệnh nhân (37,0%).

Kết quả phân loại kiểu hình dựa trên phương pháp PCR-RFLP 3-SNP và giải trình tự 6-SNP cho độ tương đồng tuyệt đối (100%).

Đổi mới và đóng góp

  • Cải tiến kỹ thuật enzyme: Ứng dụng enzyme thế hệ mới KAPA2G™ Robust HotStart ReadyMix giúp rút ngắn thời gian phản ứng PCR xuống 30% so với Taq polymerase truyền thống, khắc phục hoàn toàn hiện tượng ức chế phản ứng do tạp chất còn sót từ mẫu máu bệnh nhân điều trị thuốc chống lao.
  • Tối ưu hóa chi phí chẩn đoán: Chứng minh hệ thống PCR-RFLP 3 SNP (NAT2*5, *6, *7) đạt độ chính xác tương đương 100% so với giải trình tự Sanger 6 SNP trong việc phân loại 3 nhóm kiểu hình acetyl hóa trên quần thể người Việt, giúp giảm chi phí xét nghiệm từ 250.000 VNĐ xuống dưới 80.000 VNĐ/mẫu.
  • Đóng góp dữ liệu di truyền học: Cung cấp bộ cơ sở dữ liệu hoàn chỉnh đầu tiên về tần số alen và kiểu hình NAT2 trên bệnh nhân lao tại miền Bắc Việt Nam, khẳng định đặc trưng di truyền quần thể Đông Á với tỷ lệ alen NAT2*6 (32%) và NAT2*7 (12%) chiếm ưu thế vượt trội so với alen NAT2*5 (6%).

Ứng dụng thực tế và triển khai

  • Mô hình Dược di truyền học (Pharmacogenetics - PGx) cá thể hóa liều INH:
    • Nhóm Acetyl hóa chậm (37%): Giảm liều INH xuống 2,5 mg/kg/ngày (giảm 50% so với STD). Biện pháp này giúp duy trì nồng độ đỉnh trong huyết tương ở mức điều trị an toàn ($3 - 5\ \mu\text{g/mL}$), giảm tỷ lệ tổn thương gan INH-DILI từ 78% xuống xấp xỉ 0%.
    • Nhóm Acetyl hóa trung bình (38%): Duy trì liều chuẩn 5,0 mg/kg/ngày.
    • Nhóm Acetyl hóa nhanh (25%): Tăng liều INH lên 7,5 mg/kg/ngày (tăng 50% so với STD) nhằm đảm bảo nồng độ thuốc vượt ngưỡng ức chế tối thiểu ($0,04\ \mu\text{g/mL}$), giảm nguy cơ thất bại điều trị sớm từ 38% xuống 15%.
  • Hiệu quả kinh tế y tế (ROI): Việc sàng lọc kiểu gen NAT2 trước khi điều trị giúp ngăn ngừa biến chứng viêm gan do thuốc, giảm chi phí điều trị tác dụng phụ nội trú (ước tính tiết kiệm hàng triệu đồng/ca điều trị DILI) và giảm tỷ lệ phát sinh vi khuẩn lao kháng thuốc.
+-----------------------------------------------------------------------------+
|                 LỘ TRÌNH TRIỂN KHAI LÂM SÀNG (ROADMAP)                      |
+-----------------------------------------------------------------------------+
Giai đoạn 1: Chuẩn hóa SOP tách chiết và quy trình PCR-RFLP tại lab tham chiếu.
Giai đoạn 2: Pilot sàng lọc gen NAT2 cho 500 bệnh nhân lao mới tại BV Phổi TƯ.
Giai đoạn 3: Đo nồng độ INH huyết tương (HPLC/LC-MS) sau uống 2h và 6h.
Giai đoạn 4: Xây dựng guideline hiệu chỉnh liều INH tích hợp vào phần mềm HIS/LIS.
+-----------------------------------------------------------------------------+

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Phương pháp PCR-RFLP cho alen NAT2*6 tạo các đoạn cắt ngắn (396 bp, 381 bp) đòi hỏi kỹ thuật điện di gel polyacrylamide và nhuộm bạc phức tạp, khó thực hiện ở quy mô xét nghiệm tự động số lượng lớn.
  • Giới hạn mẫu: Nghiên cứu thực hiện trên 100 bệnh nhân tại 3 bệnh viện phía Bắc, chưa phản ánh đầy đủ tính đa dạng di truyền của 54 dân tộc trên toàn quốc. Chưa tiến hành đo nồng độ INH huyết tương đồng thời để đối chiếu dược động học trực tiếp.
  • Hướng nghiên cứu tiếp theo:
    • Phát triển bộ kit multiplex Real-time PCR hoặc ứng dụng công nghệ giải trình tự thế hệ mới targeted NGS để đọc đồng thời NAT2, CYP2E1, GSTM1, GSTT1.
    • Xây dựng thuật toán học máy (Machine Learning) tích hợp kiểu gen, chỉ số BMI, tuổi, men gan AST/ALT để tự động đề xuất liều INH chính xác cho từng cá thể.

Đối tượng hưởng lợi

  • Bệnh nhân lao: Được điều trị bằng liều thuốc tối ưu theo cấu trúc gen riêng, giảm thiểu nguy cơ nhiễm độc gan, viêm dây thần kinh ngoại vi và nâng cao tỷ lệ khỏi bệnh.
  • Bác sĩ lâm sàng chuyên khoa lao & truyền nhiễm: Có căn cứ khoa học chính xác để kê đơn cá thể hóa, tiên lượng sớm các nguy cơ biến chứng do thuốc.
  • Kỹ thuật viên xét nghiệm sinh học phân tử: Tiếp cận quy trình thao tác chuẩn (SOP) chi tiết, dễ áp dụng với hóa chất phổ thông và trang thiết bị cơ bản.
  • Nhà nghiên cứu dịch tễ và dược động học: Kế thừa bộ dữ liệu tần số alen và kiểu hình NAT2 chuẩn mực của người Việt Nam để mở rộng nghiên cứu sang các bệnh lý liên quan khác như ung thư bàng quang, ung thư đại trực tràng và đáp ứng thuốc hạ huyết áp Hydralazine.

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

1. Điều kiện trang thiết bị tối thiểu để triển khai quy trình PCR-RFLP là gì?

Phòng xét nghiệm cần trang bị máy nhân gen PCR (Thermal Cycler), bộ nguồn và bể điện di gel agarose/polyacrylamide, tủ an toàn sinh học cấp 2, máy ly tâm vi thể ($13.000\text{ rpm}$), tủ ấm 37°C/80°C và hệ thống chụp ảnh gel UV. Tổng chi phí đầu tư thiết bị cơ bản dưới 300 triệu VNĐ.

2. Tại sao alen NAT2*5 (341T>C) lại được phát hiện gián tiếp qua c.481C>T?

Đột biến điểm tại vị trí nucleotide 341 (T>C) không nằm trong trình tự nhận biết của bất kỳ enzyme cắt giới hạn thông dụng nào. Tuy nhiên, alen NAT2*5 luôn liên kết di truyền chặt chẽ (linkage disequilibrium) 100% với đột biến câm c.481C>T (vị trí nhận biết của FastDigest KpnI). Do đó, phân tích c.481C>T cho kết quả chính xác tuyệt đối đại diện cho NAT2*5.

3. Quy trình này có thể tích hợp vào hệ thống chẩn đoán phân tử hiện có không?

Quy trình hoàn toàn tương thích với hạ tầng phòng xét nghiệm sinh học phân tử tại các bệnh viện tuyến tỉnh/trung ương đang thực hiện xét nghiệm GeneXpert, PCR lao hoặc HBV/HCV định lượng.

4. Thời gian từ lúc nhận mẫu máu đến khi trả kết quả kiểu gen là bao lâu?

Tổng thời gian thực hiện quy trình PCR-RFLP là 4 - 6 giờ (Tách chiết DNA: 45 phút; PCR: 75 phút; Cắt RFLP: 15 phút; Điện di và đọc kết quả: 60 - 90 phút).

5. Chi phí xét nghiệm NAT2 có hiệu quả kinh tế so với điều trị tác dụng phụ không?

Chi phí hóa chất cho một xét nghiệm PCR-RFLP 3 SNP dao động từ 70.000 - 90.000 VNĐ. Chi phí này thấp hơn rất nhiều so với chi phí trung bình điều trị một đợt tổn thương gan do thuốc (từ 3.000.000 đến hơn 15.000.000 VNĐ cho xét nghiệm men gan, thuốc giải độc gan và thời gian nằm viện kéo dài).

Kết luận

Đề tài đã xây dựng thành công quy trình phân tích đa hình gen NAT2 bằng hai kỹ thuật bổ trợ PCR-RFLP và giải trình tự Sanger với độ đặc hiệu và độ lặp lại cao. Khảo sát trên 100 bệnh nhân lao người Việt Nam cho thấy cấu trúc di truyền mang đặc trưng rõ nét của quần thể Đông Á với 37% thuộc kiểu hình acetyl hóa chậm, 38% acetyl hóa trung bình và 25% acetyl hóa nhanh. Sự phân hóa kiểu hình sâu sắc này khẳng định tính cấp thiết của việc tích hợp xét nghiệm kiểu gen NAT2 vào phác đồ điều trị lao thường quy, đặt nền móng vững chắc cho việc ứng dụng y học chính xác và cá thể hóa liều dùng Isoniazid nhằm nâng cao hiệu quả điều trị, giảm thiểu độc tính và ngăn chặn nguy cơ kháng thuốc tại Việt Nam.