Tổng quan nghiên cứu

Bệnh Alzheimer (Alzheimer's Disease - AD) hiện là nguyên nhân hàng đầu gây ra hội chứng sa sút trí tuệ ở người cao tuổi, ảnh hưởng đến khoảng 50 triệu người trên toàn cầu vào năm 2020 và được dự báo sẽ tăng lên 152 triệu người vào năm 2050 theo các thống kê y tế quốc tế. Quá trình bệnh lý diễn tiến phức tạp qua nhiều giai đoạn, khởi đầu từ trạng thái nhận thức bình thường (Cognitive Normal - CN), chuyển tiếp qua suy giảm nhận thức nhẹ (Mild Cognitive Impairment - MCI), và tiến triển thành sa sút trí tuệ do AD. Việc chẩn đoán sớm và chính xác đóng vai trò then chốt trong việc can thiệp kịp thời nhằm làm chậm tiến trình thoái hóa thần kinh. Tuy nhiên, thách thức lớn nhất nằm ở tính không đồng nhất giữa các bệnh nhân và mối quan hệ phi tuyến tính phức tạp giữa các nguồn dữ liệu sinh học khác nhau.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Truyền thông và Dữ liệu của tác giả Vũ Duy Thanh (năm 2023) tập trung giải quyết bài toán chẩn đoán tự động AD thông qua việc xây dựng khung học máy tích hợp đa phương thức (Multimodal Data Fusion). Nghiên cứu tích hợp đồng thời 4 nguồn chỉ dấu sinh học: hình ảnh cộng hưởng từ cấu trúc (MRI), cắt lớp phát xạ positron chuyển hóa glucose (FDG-PET), nồng độ protein trong dịch não tủy (CSF) và dữ liệu giải trình tự toàn bộ hệ gene (Whole Genome Sequencing - WGS) với 6.244 biến thể di truyền SNP.

Nghiên cứu được thử nghiệm trên tập dữ liệu chuẩn ADNI (Alzheimer's Disease Neuroimaging Initiative) gồm 331 đối tượng (121 người bình thường CN, 100 người suy giảm nhận thức nhẹ MCI và 110 bệnh nhân AD). Kết quả nghiên cứu đem lại ý nghĩa quan trọng khi giải quyết triệt để hạn chế của các phương pháp kết hợp tuyến tính thông thường, nâng cao độ chính xác phân loại và cải thiện diện tích dưới đường cong ROC (AUC), từ đó cung cấp công cụ hỗ trợ quyết định lâm sàng có độ tin cậy cao.

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 kết hợp giữa lý thuyết hàm hạt nhân (Kernel Methods), phân tích dữ liệu đa tuyến tính (Tensor Decomposition) và mạng học sâu (Deep Learning). Về lý thuyết hàm hạt nhân, ma trận tương đồng được biểu diễn thông qua hàm RBF (Radial Basis Function) hoặc mô hình rừng ngẫu nhiên (Random Forest), cho phép ánh xạ dữ liệu từ không gian đặc trưng ban đầu sang không gian Hilbert nhiều chiều nhằm phân tách phi tuyến tính.

Về lý thuyết phân rã tensor, nghiên cứu áp dụng mô hình phân rã CANDECOMP/PARAFAC (CP) không âm đối với tensor bậc 3. Kỹ thuật này phân tích tensor chứa các ma trận kernel chuẩn hóa từ các phương thức khác nhau thành tổng của các tích ngoài giữa các vector thành phần không âm, giúp biểu diễn cấu trúc tương quan đa chiều mà vẫn duy trì tính diễn giải sinh học.

Khung lý thuyết tích hợp 3 nhóm khái niệm then chốt:

  1. Chỉ dấu sinh học thần kinh: Bao gồm nồng độ peptide Amyloid-beta 42 (Aβ42), Total-tau (T-tau) và Phospho-tau181 (P-tau181) trong dịch não tủy, phản ánh sự tích tụ mảng amyloid và thoái hóa sợi trục thần kinh.
  2. Hình ảnh học chức năng và cấu trúc não bộ: Đo lường mức tiêu thụ glucose của não (CMRglc) qua ảnh FDG-PET và đo lường mức độ teo thể tích mô xám (Grey Matter - GM), mô trắng (White Matter - WM) qua ảnh T1-weighted MRI.
  3. Biến thể di truyền đơn nucleotide (SNP): Mã hóa số lượng alen phụ (minor alleles) nhận các giá trị rời rạc 0, 1 hoặc 2 tương ứng với mức độ biến dị di truyền của từng cá thể.

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

Nguồn dữ liệu được trích xuất từ cơ sở dữ liệu quốc tế ADNI với cỡ mẫu gồm 331 đối tượng được chọn lọc kỹ lưỡng, đảm bảo có đầy đủ dữ liệu ảnh chụp và xét nghiệm sinh học. Quy trình tiền xử lý hình ảnh được chuẩn hóa nghiêm ngặt:

  • Ảnh sMRI gốc (kích thước 176x240x256 voxel, độ phân giải 1,2x1x1 mm³) được phân đoạn mô não, tái tạo bề mặt vỏ não và đưa về không gian chuẩn MNI bằng công cụ CAT12 thuộc nền tảng SPM.
  • Ảnh FDG-PET gốc (kích thước 160x160x96 voxel, độ phân giải 1,5x1,5x1,5 mm³) được đồng đăng ký không gian (co-registration) với ảnh sMRI thông qua bộ công cụ ANTs và chuẩn hóa cường độ hấp thu theo vùng tham chiếu chuẩn, đưa về kích thước đồng nhất 113x137x113 voxel.

Phương pháp phân tích đề xuất gồm hai kỹ thuật học hạt nhân mới:

  • Học hạt nhân sâu (Deep Kernel Learning - DKL): Sử dụng mạng tích chập hoàn toàn với các bộ lọc kích thước 1x1 kết hợp hàm kích hoạt ReLU để học trọng số phi tuyến tính giữa các ma trận kernel.
  • Học hạt nhân tensor (Tensor Kernel Learning - TKL): Xây dựng tensor ma trận tương đồng bậc 3 và tối ưu hóa hàm mục tiêu phân rã CP không âm bằng thuật toán luân phiên bình phương tối thiểu (Alternating Least Squares).

Sau khi thu được kernel hợp nhất, thuật toán học đa tạp UMAP (Uniform Manifold Approximation and Projection) được ứng dụng để giảm số chiều về không gian biểu diễn 10 chiều trước khi đưa vào bộ phân loại SVM (Support Vector Machine) và Random Forest với 500 lần lặp kiểm thử chéo để đánh giá độ ổn định.

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

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

Phân tích hiệu năng phân loại trên các phương thức đơn lẻ qua 500 lần thử nghiệm lặp lại cho thấy sự chênh lệch rõ rệt giữa các chỉ dấu sinh học:

  • Phương thức FDG-PET đạt hiệu quả đơn lẻ cao nhất trong việc phân biệt giữa nhóm bệnh nhân AD và nhóm đối chứng bình thường CN, với độ chính xác trung bình đạt 87,74% và chỉ số diện tích dưới đường cong ROC (AUC) đạt 93,80%.
  • Phương thức MRI cấu trúc (đặc trưng thể tích mô xám GM) đứng thứ hai với độ chính xác đạt 80,54% và AUC đạt 90,38%. Khi kết hợp thêm đặc trưng thể tích mô trắng (GM+WM), khả năng phân biệt giữa giai đoạn MCI và AD được cải thiện đáng kể.
  • Phương thức xét nghiệm dịch não tủy (CSF) đạt độ chính xác 80,20% và AUC đạt 88,40%, thể hiện khả năng nhận diện tốt ở một phân nhóm bệnh nhân nhất định.
  • Ngược lại, dữ liệu 6.244 biến thể di truyền SNP đơn lẻ chỉ đạt độ chính xác khoảng 50,00% trên cả 3 tác vụ phân loại nhị phân (CN vs AD, CN vs MCI, MCI vs AD), gần tương đương với mức phân loại ngẫu nhiên.

Khi áp dụng hai phương pháp hợp nhất DKL và TKL, mô hình đã giải quyết triệt để sự nhiễu loạn của dữ liệu di truyền, kết hợp hài hòa ưu thế của hình ảnh PET và MRI, nâng cao độ chính xác tổng thể và cải thiện rõ rệt chỉ số AUC ở các bài toán chẩn đoán sớm khó khăn như phân loại CN và MCI.

Thảo luận kết quả

Cơ chế mang lại sự vượt trội của phương pháp DKL nằm ở khả năng ánh xạ phi tuyến tính thông qua các lớp tích chập sâu 1x1, giúp phát hiện mối tương tác phức tạp giữa sự suy giảm chuyển hóa glucose ở vùng vỏ não đỉnh-thái dương (trên ảnh FDG-PET) và mức độ teo hồi hải mã (trên ảnh MRI). Trong khi đó, phương pháp TKL dựa trên phân rã CP không âm cho phép bảo toàn cấu trúc đa chiều của dữ liệu, ngăn ngừa tình trạng quá khớp (overfitting) và loại bỏ hiện tượng các phương thức có số lượng biến lớn áp đảo các phương thức có số lượng biến nhỏ.

Dữ liệu ma trận kernel được trực quan hóa thông qua biểu đồ nhiệt (heatmap) và không gian nhúng đa tạp UMAP cho thấy sự tách biệt rõ ràng giữa các nhóm bệnh lý:

  • Trong ma trận kernel có giám sát (xây dựng bằng Random Forest 200 cây), ranh giới giữa cụm bệnh nhân CN và AD xuất hiện sắc nét, phân tách thành hai khối trực giao rõ rệt dọc theo đường chéo chính.
  • Biểu đồ phân tán UMAP 2 chiều cho thấy các đối tượng AD phân hóa thành 2 phân nhóm rõ rệt, phản ánh các mức độ thoái hóa thần kinh và giai đoạn tiến triển bệnh khác nhau.
  • Bảng so sánh giữa các hàm mất mát (loss functions) trong mạng DKL chứng minh rằng các hàm tối ưu hóa căn chỉnh hạt nhân (Kernel Alignment) đạt độ hội tụ nhanh và ổn định hơn so với hàm mất mát MSE truyền thống qua các mốc huấn luyện từ 100 đến 520 chu kỳ lặp.

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

Nhằm chuyển hóa các kết quả nghiên cứu thành các giải pháp ứng dụng thực tiễn trong y tế số và chẩn đoán lâm sàng, nghiên cứu đề xuất 4 khuyến nghị trọng tâm:

  1. Chuẩn hóa quy trình tiền xử lý dữ liệu hình ảnh y tế tự động: Các bệnh viện và trung tâm chẩn đoán hình ảnh cần áp dụng quy trình tiền xử lý đồng bộ sử dụng CAT12 và ANTs để chuyển đổi toàn bộ dữ liệu sMRI và PET về không gian giải phẫu chuẩn MNI kích thước 113x137x113 voxel. Mục tiêu nhằm giảm hơn 40% độ sai lệch do khác biệt giữa các dòng máy quét, hoàn thiện quy chuẩn kỹ thuật trong vòng 6 tháng.

  2. Ứng dụng mô hình Deep Kernel Learning trong sàng lọc giai đoạn sớm: Đội ngũ kỹ sư y sinh và chuyên gia công nghệ thông tin y tế cần tích hợp kiến trúc DKL vào hệ thống chẩn đoán hỗ trợ máy tính (CAD). Mục tiêu nâng cao độ nhạy phát hiện bệnh nhân suy giảm nhận thức nhẹ (MCI) có nguy cơ chuyển thành AD lên trên 85% trong giai đoạn triển khai 12 tháng tại các cơ sở y tế tuyến trung ương.

  3. Xây dựng nền tảng học bán giám sát cá nhân hóa dựa trên TKL: Các nhà phát triển phần mềm y tế nên tận dụng tính năng phân rã tensor không âm của TKL để kết hợp dữ liệu bệnh nhân chưa gán nhãn với dữ liệu chuẩn ADNI. Giải pháp này giúp cập nhật hồ sơ nguy cơ cá nhân hóa với thời gian xử lý dưới 10 giây cho mỗi ca bệnh mà không đòi hỏi huấn luyện lại toàn bộ mô hình.

  4. Thiết lập khung đánh giá đa chỉ dấu kết hợp di truyền và sinh hóa: Hội đồng chuyên môn thần kinh học và tin sinh học cần xây dựng hướng dẫn chẩn đoán phối hợp, tích hợp tỷ lệ sinh hóa CSF (Aβ42/Tau) cùng với 6.244 biến thể SNP để đánh giá nguy cơ tiền lâm sàng, hoàn thành tài liệu hướng dẫn trong khung thời gian 12 đến 18 tháng.

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

Nội dung và phương pháp trong luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chính:

  • Kỹ sư Trí tuệ nhân tạo Y tế và Khoa học Dữ liệu (Medical AI & Data Science Engineers): Nắm bắt kiến trúc học sâu trên ma trận hạt nhân (DKL) và kỹ thuật phân rã tensor bậc 3 (TKL), ứng dụng trực tiếp vào việc xử lý các tập dữ liệu đa phương thức, đa chiều phức tạp trong y học và các ngành kỹ thuật khác.
  • Bác sĩ chuyên khoa Thần kinh và Chẩn đoán hình ảnh: Tiếp cận góc nhìn định lượng mới về mối tương quan giữa mức độ hấp thu FDG-PET vùng đai sau, hồi hải mã và nồng độ protein Aβ42 trong dịch não tủy, hỗ trợ nâng cao độ chính xác trong phân loại giai đoạn bệnh.
  • Nhà nghiên cứu Tin sinh học và Di truyền học (Bioinformatics Researchers): Tham khảo phương pháp mã hóa và xử lý ma trận khoảng cách cho 6.244 biến thể SNP, kết hợp dữ liệu hệ gene với kiểu hình hình ảnh học thần kinh (Imaging Genetics).
  • Đơn vị phát triển Hệ thống Hỗ trợ Ra quyết định Lâm sàng (CDSS Developers): Ứng dụng quy trình tích hợp thuật toán giảm chiều UMAP và bộ phân loại SVM/Random Forest để xây dựng các mô-đun phần mềm chẩn đoán bệnh thoái hóa thần kinh theo thời gian thực.

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

Tại sao mô hình Deep Kernel Learning (DKL) lại hiệu quả hơn phương pháp nối chuỗi đặc trưng truyền thống?

Phương pháp nối chuỗi đặc trưng (Concatenation) làm tăng đột biến số chiều dữ liệu, dễ dẫn đến hiện tượng quá khớp và khiến các phương thức có số lượng biến lớn áp đảo các phương thức nhỏ. Mạng DKL sử dụng tích chập 1x1 và hàm kích hoạt ReLU trên ma trận kernel $N \times N$, giúp nắm bắt trọn vẹn các tương tác phi tuyến tính giữa các phương thức mà vẫn giữ nguyên kích thước ma trận.

Phân rã Tensor CP không âm trong phương pháp TKL mang lại lợi ích gì?

Phân rã CP không âm phân tích tensor kernel bậc 3 thành các thành phần cơ sở mang giá trị dương, phản ánh mức độ đóng góp và tính nhất quán giữa các nguồn dữ liệu MRI, PET, CSF và SNP. Ràng buộc không âm giúp mô hình dễ dàng giải thích về mặt sinh học, đồng thời cho phép triển khai học bán giám sát trên cả tập dữ liệu huấn luyện và kiểm thử.

Tại sao dữ liệu di truyền SNP chỉ đạt độ chính xác khoảng 50% khi phân loại đơn lẻ?

Dữ liệu di truyền gồm 6.244 biến thể SNP có tính nhiễu cao và tác động của từng alen riêng lẻ lên kiểu hình sa sút trí tuệ thường rất nhỏ. Khi đứng đơn lẻ, tín hiệu di truyền không đủ mạnh để bộ phân loại phân tách ranh giới nhị phân, nhưng khi được tích hợp qua kernel DKL, thông tin di truyền đóng vai trò bổ trợ quan trọng cho các chỉ dấu hình ảnh.

Chỉ dấu sinh học nào thể hiện hiệu năng chẩn đoán phân biệt AD và CN cao nhất?

Hình ảnh chuyển hóa FDG-PET đạt hiệu năng đơn lẻ cao nhất với độ chính xác 87,74% và chỉ số AUC đạt 93,80% qua 500 lần thử nghiệm. Điều này phản ánh chính xác cơ chế bệnh sinh khi sự suy giảm tiêu thụ glucose tại các vùng vỏ não đỉnh-thái dương và hồi đai sau là dấu hiệu đặc trưng xuất hiện rất sớm ở bệnh nhân AD.

Việc áp dụng thuật toán UMAP có nguy cơ gây rò rỉ dữ liệu trong quá trình đánh giá không?

Đối với kernel có giám sát được tạo bởi Random Forest, nghiên cứu cố định hoàn toàn chỉ số của tập huấn luyện và tập kiểm tra trong suốt quá trình xây dựng kernel, chiếu UMAP và phân loại cuối cùng. Việc bảo toàn phân chia tập dữ liệu này loại bỏ hoàn toàn nguy cơ rò rỉ dữ liệu (data leakage), đảm bảo tính khách quan của kết quả đánh giá.

Kết luận

  • Luận văn đã xây dựng thành công khung tích hợp dữ liệu đa phương thức y sinh toàn diện kết hợp hình ảnh sMRI, FDG-PET, chỉ dấu dịch não tủy CSF và 6.244 biến thể di truyền SNP.
  • Đề xuất hai phương pháp học hạt nhân tiên tiến gồm Deep Kernel Learning (DKL) giúp khai phá tương tác phi tuyến tính và Tensor Kernel Learning (TKL) tối ưu hóa biểu diễn đa tuyến tính.
  • Đánh giá thực nghiệm quy mô lớn trên 331 đối tượng thuộc tập dữ liệu quốc tế ADNI qua 500 lần lặp kiểm thử, chứng minh sự ổn định vượt trội so với các mô hình đơn phương thức.
  • Xác định rõ vai trò vượt trội của FDG-PET (độ chính xác 87,74%, AUC 93,80%) và thể tích mô xám MRI (độ chính xác 80,54%) trong việc định hình không gian phân loại bệnh lý.
  • Mở ra hướng tiếp cận học bán giám sát và phân tích cá nhân hóa trong việc chẩn đoán sớm và theo dõi tiến trình thoái hóa thần kinh ở người cao tuổi.

Đóng góp cốt lõi của nghiên cứu là cung cấp một giải pháp toán học và học máy vững chắc để dung hòa các nguồn dữ liệu y sinh có cấu trúc và thang đo hoàn toàn khác biệt. Trong lộ trình 6 đến 12 tháng tới, mô hình có thể tiếp tục được mở rộng để đánh giá trên các tập dữ liệu lâm sàng đa trung tâm độc lập. Các nhà nghiên cứu và kỹ sư y sinh quan tâm có thể khai thác mã nguồn và khung phương pháp luận này để mở rộng cho các bài toán phân loại bệnh lý thoái hóa thần kinh phức tạp khác.