Tổng quan về luận án
Ung thư vú (Breast Cancer - BC) là bệnh lý ác tính có tỷ lệ mắc mới và tử vong hàng đầu ở nữ giới trên toàn cầu. Thống kê dịch tễ học của GLOBOCAN ghi nhận riêng năm 2012 có khoảng 1,67 triệu ca mắc mới (chiếm 25% tổng số các ca ung thư ở nữ) và 522.000 trường hợp tử vong. Tại Việt Nam, tỷ lệ chuẩn hóa theo tuổi của ung thư vú tăng nhanh rõ rệt từ 13,8/100.000 dân (năm 2000) lên 28,1/100.000 dân (năm 2010), với đặc thù khởi phát ở độ tuổi trung bình khoảng 50 tuổi – sớm hơn đáng kể so với quần thể các nước phương Tây (60–70 tuổi). Dù các dấu ấn sinh học phân tử như thụ thể estrogen (ER), thụ thể progesterone (PR), thụ thể yếu tố tăng trưởng biểu bì 2 (HER2/neu), CA15-3 và CA27-29 đã được ứng dụng trong phác đồ điều trị, chúng chủ yếu phục vụ giai đoạn xâm lấn và tiên lượng sau can thiệp. Trong khi đó, các kỹ thuật sàng lọc tiêu chuẩn hiện hành như nhũ ảnh (mammography) và thăm khám lâm sàng bộc lộ tỷ lệ dương tính giả tích lũy sau 10 năm theo dõi lần lượt lên tới 49,1% và 22,3%, đặt ra yêu cầu cấp thiết về việc phát hiện các chỉ thị sinh học phân tử mới có độ nhạy và độ đặc hiệu cao hơn.
Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ thực tế: hệ gen ty thể (mitochondrial DNA - mtDNA) vốn đóng vai trò trung tâm trong chuỗi truyền điện tử, quá trình phosphoryl hóa oxy hóa (OXPHOS) và chu trình chết theo chương trình (apoptosis), nhưng các biến đổi di truyền và biểu sinh của mtDNA ở bệnh nhân ung thư vú người Việt Nam chưa từng được khảo sát một cách toàn diện và có hệ thống. Các công bố quốc tế về biến đổi số bản sao mtDNA, mất đoạn lớn $\Delta\text{mtDNA}^{4977}$ và đột biến đơn nucleotide (SNP/somatic mutation) trên các gen mã hóa (ATP6, ND1, ND3) và tARN còn nhiều mâu thuẫn giữa các chủng tộc (như nghiên cứu của Tseng và cs., 2006; Fan và cs., 2009; Dani và cs., 2004; Canter và cs., 2005). Tại Việt Nam, các nghiên cứu trước đây chủ yếu khảo sát đa hình vùng điều khiển D-loop trên quần thể người bình thường hoặc các bệnh lý ty thể thần kinh (LHON, MELAS), ngoại trừ khảo sát sơ bộ trên 2 ca bệnh của Đái Duy Ban và cs. (2003).
Để giải quyết khoảng trống học thuật này, luận án tiến sĩ chuyên ngành Nhân chủng học (Mã số: 62420102) của nghiên cứu sinh Nguyễn Thị Tú Linh, thực hiện dưới sự hướng dẫn khoa học của PGS. Trịnh Hồng Thái và PGS. Tạ Văn Tờ tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội phối hợp với Bệnh viện K Trung ương, đã xác lập các câu hỏi nghiên cứu và giả thuyết khoa học trọng tâm:
- Câu hỏi nghiên cứu 1 (RQ1): Tần suất, mức độ biến đổi số lượng bản sao mtDNA và mất đoạn lớn $\Delta\text{mtDNA}^{4977}$ xuất hiện như thế nào trong mô u ung thư biểu mô ống tuyến vú so với mô lành liền kề, mô u tuyến vú lành tính và máu ngoại vi của người Việt Nam?
- Giả thuyết 1 (H1): Mô ung thư vú có sự suy giảm đáng kể số bản sao mtDNA và mức độ mất đoạn lớn phản ánh tình trạng tái lập trình chuyển hóa tế bào và tổn thương tích lũy qua stress oxy hóa.
- Câu hỏi nghiên cứu 2 (RQ2): Các biến đổi nucleotide trên các gen mã hóa tiểu phần chuỗi hô hấp tế bào (ATP6, ND1, ND3) và các gen tARN ty thể diễn ra với tần suất nào và có mối liên hệ tương quan gì với các chỉ số giải phẫu bệnh học lâm sàng (kích thước khối u T, tình trạng di căn hạch N, độ mô học, giai đoạn TNM)?
- Giả thuyết 2 (H2): Các đột biến điểm và đa hình nucleotide trên các gen ATP6, ND1, ND3, tARN làm biến đổi cấu trúc bậc hai của RNA hoặc thay đổi chuỗi polypeptide, liên quan trực tiếp đến độ tiến triển ác tính của khối u.
- Câu hỏi nghiên cứu 3 (RQ3): Liệu việc tích hợp dữ liệu số bản sao mtDNA và mức độ mất đoạn lớn vào mô hình phân tích hồi quy đa biến có thể tạo ra một công cụ sinh học định lượng có giá trị dự báo nguy cơ ung thư vú hay không?
- Giả thuyết 3 (H3): Mô hình phối hợp đa biến số bản sao và tỷ lệ mất đoạn mtDNA cung cấp độ nhạy, độ đặc hiệu vượt trội và giảm thiểu tỷ lệ dương tính giả/âm tính giả so với các mô hình đơn biến truyền thống.
Khung lý thuyết của công trình được định vị trên Thuyết chuyển hóa ung thư của Otto Warburg (Warburg Effect, 1924/1956), Thuyết gốc tự do và lão hóa ty thể của Denham Harman (1956), kết hợp với Thuyết ngưỡng dị tế bào chất (Heteroplasmy Threshold Effect) của Douglas Wallace (1999, 2005). Nghiên cứu được triển khai trên cỡ mẫu gồm 102 cặp mô u và mô liền kề (cách rìa khối u $\ge 5\text{ cm}$) từ bệnh nhân ung thư biểu mô ống tuyến vú nguyên phát, đối chứng với 20 mẫu mô và máu của bệnh nhân u tuyến vú lành tính cùng 65 mẫu máu từ người hiến máu khỏe mạnh tại Viện Huyết học - Truyền máu Trung ương (thời gian thu thập 12/2012 – 12/2013). Đóng góp mang tính đột phá của luận án là lần đầu tiên tại Việt Nam thiết lập được ngân hàng dữ liệu di truyền ty thể hoàn chỉnh cho ung thư vú và xây dựng thành công mô hình hồi quy logistic đa biến định lượng nguy cơ ác tính dựa trên động học hệ gen ty thể.
Literature Review và Positioning
Nghiên cứu về vai trò của ty thể trong ung thư học hiện đại tập trung vào ba luồng học thuyết và thực nghiệm chính:
-
Động học số bản sao mtDNA (mtDNA Copy Number Dynamics): Hệ gen ty thể người gồm phân tử ADN vòng mạch kép kích thước 16.569 bp, mã hóa 13 polypeptide thuộc chuỗi hô hấp tế bào, 22 tARN và 2 rARN. Mỗi tế bào chứa từ hàng trăm đến hàng nghìn bản sao mtDNA phụ thuộc vào nhu cầu chuyển hóa năng lượng. Trong ung thư, sự biến đổi số bản sao thể hiện tính chất đa chiều: sự gia tăng số bản sao được ghi nhận ở ung thư buồng trứng (Wang và cs., 2006), ung thư đại trực tràng (Cui và cs., 2013), trong khi sự suy giảm rõ rệt lại xuất hiện ở ung thư biểu mô tế bào thận (Meierhofer và cs., 2004) và ung thư biểu mô tế bào gan (Lee và cs., 2004). Trên bệnh nhân ung thư vú, Tseng và cs. (2006) tại Đài Loan ghi nhận 63% bệnh nhân giảm số bản sao mtDNA ở mô u, đồng thuận với nghiên cứu của Fan và cs. (2009) tại Trung Quốc (82% suy giảm). Ngược lại, một số nghiên cứu trên mẫu máu ngoại vi như của Shen và cs. (2009) hay Lemnrau và cs. (2015) lại nhận thấy số bản sao mtDNA cao có tương quan thuận với nguy cơ mắc bệnh ($P_{\text{trend}} = 0,007$), tạo nên sự tranh luận học thuật kéo dài về giá trị tiên lượng của chỉ số này giữa mô đích và dịch lỏng sinh học.
-
Cơ chế mất đoạn lớn $\Delta\text{mtDNA}^{4977}$ và Hiện tượng "Vùng cung mất đoạn chính" (Major Deletion Arc): Mất đoạn 4977 bp xảy ra giữa hai đoạn lặp trực tiếp 13 bp ở vị trí nucleotide 8470–8482 và 13447–13459, loại bỏ hoàn toàn 5 gen tARN và 7 gen mã hóa thuộc phức hệ I, IV và V. Về mặt bệnh học, tồn tại sự bất đồng gay gắt giữa hai luồng quan điểm:
- Quan điểm thứ nhất: Dani và cs. (2004) phát hiện $\Delta\text{mtDNA}^{4977}$ hiện diện ở 24% mô u và 77% mô liền kề ở 87 bệnh nhân ung thư vú, nhưng cho rằng sự có mặt ở mô u là do tạp nhiễm mô bình thường và kết luận khối u ác tính đào thải mất đoạn này. Ngược lại, Ye và cs. (2008) ứng dụng real-time PCR phát hiện mất đoạn này trong 100% mẫu mô (mức độ từ 0,000149% đến 7,0%) nhưng không thấy khác biệt có ý nghĩa thống kê giữa mô u và mô lành.
- Quan điểm thứ hai: Tseng và cs. (2006, 2009) chứng minh tỷ lệ mất đoạn 4977 bp ở mô lành liền kề (47–48,3%) cao hơn vượt trội so với mô u (5%), chứng minh tổn thương oxy hóa xuất hiện từ sớm ở vi môi trường mô xung quanh u ("Hiệu ứng tổn thương vùng" - Field Cancerization), trong khi các tế bào phân chia nhanh trong khối u loại bỏ phân tử mtDNA mang mất đoạn lớn nhằm bảo tồn quá trình sao chép tối thiểu.
-
Đột biến điểm các gen phức hệ I (ND1, ND3), phức hệ V (ATP6) và tARN:
- Đa hình ND3 A10398G (gây thay thế Thr $\rightarrow$ Ala tại codon 114) là tâm điểm tranh luận di truyền học chủng tộc: Canter và cs. (2005) chứng minh alen 10398A làm tăng nguy cơ ung thư vú ở phụ nữ Mỹ gốc Phi nhưng không ảnh hưởng trên người Mỹ da trắng; Darvishi và cs. (2007) và Sultana và cs. (2011) xác nhận nguy cơ liên quan alen 10398A ở phụ nữ Ấn Độ và Bangladesh. Trái lại, Covarrubias và cs. (2008), Bai và cs. (2007) trên người Mỹ da trắng, Czarnecka và cs. (2010a) trên người Ba Lan ($p = 0,0008$), và Nadiah và cs. (2012) trên người Malaysia lại khẳng định alen 10398G mới là yếu tố nguy cơ. Ngược lại, Setiawan và cs. (2008), Mosquera-Miguel và cs. (2008) và Francis và cs. (2013) không tìm thấy mối liên hệ nào.
- Đối với gen ATP6, Grzybowska-Szatkowska và cs. (2014a) ghi nhận 72% bệnh nhân ung thư vú mang biến đổi nucleotide; Ghaffarpour và cs. (2014) tìm thấy biến đổi A8860G ở 100% mẫu u. Ngược lại, Sharp và cs. (1992) hay Chintha và cs. (2014) trên 180 bệnh nhân lại không tìm thấy bất kỳ biến đổi nào của ATP6.
==========================================================================================
BẢN ĐỒ VỊ TRÍ NGHIÊN CỨU TRONG LITERATURE QUỐC TẾ
==========================================================================================
[Tranh luận Copy Number] [Tranh luận Mất đoạn 4977 bp] [Tranh luận Đa hình A10398G]
Tseng (2006) Đài Loan: -63% Dani (2004): 24% U vs 77% Liền kề Canter (2005): Alen A tăng nguy cơ (Mỹ gốc Phi)
Fan (2009) TQ: -82% Tseng (2009): 5% U vs 48% Liền kề Czarnecka (2010): Alen G tăng nguy cơ (Ba Lan)
Shen (2009): Máu ngoại vi tăng Ye (2008): 100% mẫu, ko ý nghĩa Francis (2013): Không liên quan (Nam Ấn Độ)
\ | /
\ | /
\ | /
===================================================
LUẬN ÁN NGUYỄN THỊ TÚ LINH (VIỆT NAM, 2017):
- Chuẩn hóa qPCR nồng độ tuyệt đối plasmid (ND1, ND4, HBB)
- Khảo sát hệ thống 102 cặp mô U / Liền kề + 20 Lành + 65 Máu
- Tích hợp Đa biến (Copy Number + Deletion Level) -> Khung ROC
===================================================
Luận án của Nguyễn Thị Tú Linh đã định vị chính xác vào giao điểm của các tranh luận quốc tế này, cung cấp chứng cứ thực nghiệm vững chắc từ một quần thể nhân chủng học Đông Nam Á chưa từng được phân tích hệ thống.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mở rộng và bổ sung căn bản cho các khung lý thuyết kinh điển:
- Mở rộng Thuyết chuyển hóa năng lượng Warburg (Warburg Metabolic Theory): Otto Warburg (1924, 1956) cho rằng tế bào ung thư chuyển dịch từ hô hấp hiếu khí OXPHOS sang đường phân ái khí (aerobic glycolysis) do khiếm khuyết không hồi phục của hệ thống ty thể. Luận án cung cấp bằng chứng di truyền phân tử định lượng: sự suy giảm số lượng bản sao mtDNA kết hợp với sự ức chế biểu hiện gen tổng hợp ATP (ATP6) và các phức hệ I (ND1, ND3) làm suy giảm hiệu suất tổng hợp ATP qua chuỗi chuyền điện tử, trực tiếp thúc đẩy tế bào u thích nghi với con đường đường phân nhằm tạo tiền chất sinh tổng hợp nucleotide, lipid và protein phục vụ tăng sinh tế bào ác tính.
- Bổ trợ Thuyết ngưỡng dị tế bào chất (Heteroplasmy Threshold Hypothesis): Theo Douglas Wallace (1999) và Chinnery (2000), đột biến mtDNA chỉ gây biểu hiện bệnh lý lâm sàng khi tỷ lệ phân tử đột biến vượt qua ngưỡng tới hạn (thường từ 50% đến 90%). Luận án đã lượng hóa chính xác ngưỡng biểu hiện lâm sàng của tỷ lệ mất đoạn lớn mtDNA trong mô u tuyến vú, chứng minh các mất đoạn lớn dị tế bào chất tích lũy âm thầm và đạt điểm bùng phát gây rối loạn chuỗi hô hấp tế bào trước khi biểu hiện ra kiểu hình giải phẫu bệnh ác tính.
- Củng cố Thuyết tổn thương vùng trong sinh ung thư (Field Cancerization Theory): Khởi xướng bởi Slaughter và cs. (1953), thuyết này chỉ ra rằng các mô lành xung quanh khối u đã trải qua các biến đổi tiền ác tính về mặt di truyền. Bằng chứng thực nghiệm từ luận án cho thấy mức độ mất đoạn $\Delta\text{mtDNA}^{4977}$ ở mô liền kề (cách mép khối u 5 cm) đạt tỷ lệ cao hơn nhiều so với mô u thực thể, khẳng định sự tích lũy tổn thương do các gốc tự do chứa oxy (ROS) trong vi môi trường mô tuyến vú trước khi xuất hiện khối u xâm lấn.
Các giả thuyết và mệnh đề lý thuyết được kiểm chứng cụ thể:
- Mệnh đề 1 ($P_1$): Mức độ sao chép mtDNA tỷ lệ nghịch với độ tiến triển mô học và kích thước khối u (Grade III và T3-4 có số bản sao thấp nhất).
- Mệnh đề 2 ($P_2$): Tỷ số mất đoạn lớn mô u/mô liền kề phản ánh sự chọn lọc tế bào dòng vô tính (clonal selection) trong tiến trình sinh ung thư.
- Mệnh đề 3 ($P_3$): Đột biến thay thế nucleotide trên các vòng kẹp tARN làm mất ổn định năng lượng tự do tối thiểu (Minimum Free Energy - MFE), cản trở quá trình aminoacyl hóa và dịch mã protein ty thể.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp ba trục lý thuyết: Di truyền số lượng bản sao (Copy Number Variation) – Động học cấu trúc bộ gen (Large Structural Deletions) – Đột biến đơn nucleotide chức năng (Functional Point Mutations).
+-----------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH HỆ GEN TY THỂ |
+-----------------------------------------------------------------------------------+
|
+------------------------------------+-----------------------------------+
| | |
v v v
[Trục Số Lượng Bản Sao] [Trục Cấu Trúc Mất Đoạn] [Trục Đột Biến Điểm]
- Đoạn gen chuẩn: ND1 - Đoạn mất: Major Arc - Phức hệ V: ATP6 (G9053A)
- Gen nhân đối chiếu: HBB - Mất đoạn chuẩn: 4977 bp - Phức hệ I: ND1, ND3 (A10398G)
- Phương pháp: Real-time PCR - Đoạn kiểm soát: ND4/ND1 - Cấu trúc: 22 tARN (A5536T)
- Chuẩn hóa: Plasmid curve - Kỹ thuật: PCR định tính/qPCR - Kỹ thuật: PCR-RFLP, Sanger
| | |
+------------------------------------+-----------------------------------+
|
v
+-----------------------------------------------------------------------+
| MÔ HÌNH TOÁN HỌC DỰ BÁO NGUY CƠ TÍCH HỢP |
| - Hồi quy Logistic Đa biến (Multivariate Logistic Regression) |
| - Tối ưu hóa điểm cắt ngưỡng (Cut-off Threshold Optimization) |
| - Kiểm định phân biệt đường cong ROC (Receiver Operating Curve) |
+-----------------------------------------------------------------------+
Điều kiện biên (Boundary Conditions): Khung phân tích này được thiết lập chuẩn hóa trên mô bệnh phẩm ung thư biểu mô ống tuyến vú nguyên phát (Invasive Ductal Carcinoma - IDC), loại trừ các thể ung thư biểu mô tiểu thùy hoặc ung thư di căn thứ phát từ các cơ quan khác đến vú; giới hạn phân tích trên DNA tổng số tách chiết từ mô tươi đông lạnh sâu $-80^\circ\text{C}$ và máu toàn phần chống đông EDTA.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) với thiết kế nghiên cứu bệnh chứng ghép cặp (Paired Case-Control Experimental Design) kết hợp định lượng phân tử đa cấp độ:
- Cấp độ mô đích: 102 cặp mô u ung thư vú và mô liền kề (lấy cách mép khối u $\ge 5\text{ cm}$ tại vùng mô không bị hoại tử), phẫu thuật triệt căn có vét hạch nách và được chẩn đoán mô bệnh học xác định tại Khoa Giải phẫu bệnh – Tế bào, Bệnh viện K Trung ương.
- Cấp độ mô bệnh lành tính: 20 cặp mô u tuyến vú lành tính (Fibroadenoma) và mẫu máu đối chứng tương ứng.
- Cấp độ cá thể đối chứng khỏe mạnh: 65 mẫu máu toàn phần từ người hiến máu tình nguyện tại Viện Huyết học - Truyền máu Trung ương.
Đặc điểm giải phẫu bệnh học của 102 bệnh nhân ung thư vú được phân nhóm chi tiết: độ tuổi $<50$ (38%), $\ge 50$ (62%); kích thước khối u $\le 2\text{ cm}$ (22,5%), $2,1 - 5\text{ cm}$ (60,8%), $>5\text{ cm}$ (16,7%); di căn hạch nách âm tính $N_0$ (39,6%), $N_{1-3}$ (35,6%), $N_{>3}$ (24,8%); mức độ biệt hóa: cao - Grade I (10,1%), vừa - Grade II (66,7%), kém - Grade III (23,2%); giai đoạn TNM: giai đoạn sớm 0–II (79,2%), giai đoạn tiến triển III–IV (20,8%).
Quy trình nghiên cứu rigorous
Quy trình thực nghiệm được chuẩn hóa nghiêm ngặt từ khâu thu thập đến phân tích dữ liệu:
-
Tách chiết và kiểm tra độ tinh sạch DNA:
- DNA tổng số từ mô được tách chiết bằng QIAamp DNA Mini Kit (Qiagen, Hilden, Đức); DNA từ máu tách bằng QIAamp DNA Blood Mini Kit. Quá trình ly giải mô sử dụng Proteinase K và ủ nhiệt ở $56^\circ\text{C}$ đến khi hòa tan hoàn toàn.
- Nồng độ và độ tinh sạch của DNA được đo quang phổ hấp thụ tử ngoại tại bước sóng 260 nm và 280 nm bằng máy NanoDrop 2000c (Thermo Fisher Scientific, Wilmington, USA). Tất cả các mẫu đạt tỷ số $A_{260}/A_{280}$ trong khoảng nghiêm ngặt từ $1,80$ đến $2,00$ mới được đưa vào phản ứng khuếch đại. Kiểm tra tính toàn vẹn của DNA bằng điện di trên gel agarose 0,8–1% nhuộm Ethidium Bromide (EtBr).
-
Thiết kế mồi và hệ thống phản ứng PCR:
- Các cặp mồi đặc hiệu được thiết kế bằng phần mềm Primer-BLAST dựa trên trình tự hệ gen ty thể tham chiếu chuẩn quốc tế NC_012920 (NCBI Recombinant Human Mitochondrial Genome).
- Tổng cộng 10 cặp mồi chuyên biệt được tối ưu hóa nhiệt độ gắn mồi ($T_a$) và chu trình nhiệt trên hệ thống máy GeneAmp PCR System 9700 (Applied Biosystems, Foster City, CA, USA) (xem bảng chi tiết):
| Gen đích |
Tên cặp mồi |
Kích thước đoạn (bp) |
Vị trí nucleotide |
Trình tự mồi xuôi ($5' \rightarrow 3'$) / mồi ngược ($5' \rightarrow 3'$) |
$T_a$ ($^\circ\text{C}$) |
| ND1 |
mt / ND1-1 |
433 / 335 |
3457–3889 / 4015–4349 |
gacgccataaaactcttcac / ggttggtctctgctagtgtg |
54 / 54 |
| $\Delta\text{mtDNA}^{4977}$ |
4977-1 / 4977-2 |
496 / 381 |
8167–13639 / 8196–13553 |
tcaatgctctgaaatctgtgg / gttgacctgttagggtgagaag |
54 / 60 |
| ND3 |
10398 |
246 |
10307–10552 |
cctgccactaatagttatgtc / gatatgaggtgtgagcgata |
54 |
| ATP6 |
ATP6 / 9053 |
1148 / 298 |
8197–9344 / 8802–9099 |
cagtttcatgcccatcgt / gcctagtatgaggagcgtta |
56 / 54 |
| tARN |
TW-1 / TW-2 |
940 / 982 |
4015–4954 / 5152–6133 |
cgaccctacttctaacc / atggataagattgagagagtga |
56 / 56 |
| tARN-Trp |
TW-3 |
173 |
5471–5643 |
gctactcctacctatctcc / taattaaagtggctggtt |
50 |
-
Kỹ thuật PCR-RFLP và Giải trình tự trực tiếp (Sanger Sequencing):
- Sản phẩm PCR được tinh sạch bằng hệ thống ExoSAP-IT (Affymetrix, Santa Clara, CA, USA) hoặc QIAquick Gel Extraction Kit.
- Phân tích RFLP sử dụng các enzyme cắt giới hạn hoạt tính cao FastDigest (Fermentas/Thermo Scientific): Hin6I (phát hiện biến đổi ATP6 G9053A), FspBI (sàng lọc tARN A5536T), DdeI (nhận diện ND3 A10398G) và NlaIII (xác định ND1 A4164G/T4216C). Sản phẩm cắt được điện di phân tách trên gel polyacrylamide 8–10% trong đệm 1X TBE và chụp ảnh phân giải cao trên hệ thống GelDoc (Bio-Rad Laboratories, Hercules, CA, USA).
- Các mẫu mang biến đổi được giải trình tự Sanger hai chiều trên máy tự động ABI PRISM 3100 Genetic Analyzer để xác thực đột biến.
-
Xây dựng đường chuẩn plasmid xác định số bản sao và tỷ lệ mất đoạn bằng Real-time PCR:
- Các đoạn gen ND1 (đại diện cho tổng lượng mtDNA), ND4 (nằm trong vùng mất đoạn chính, bị mất khi xảy ra $\Delta\text{mtDNA}^{4977}$) và HBB (gen Hemoglobin beta trên nhiễm sắc thể nhân số 11, đại diện cho hệ gen nhân) được nhân dòng vào vector pJET1.2/blunt bằng kit CloneJET PCR Cloning Kit và chuyển nạp vào tế bào khả biến Escherichia coli chủng DH5$\alpha$.
- Plasmid tái tổ hợp được tách chiết bằng QIAprep Spin Miniprep Kit, đo mật độ quang học để tính toán số lượng bản sao tuyệt đối theo công thức:
$$\text{Số bản sao}/\mu\text{l} = \frac{C \times 10^{-9} \times 6,022 \times 10^{23}}{M \times 660}$$
(Trong đó: $C$ là nồng độ plasmid đo được tính bằng ng/$\mu$l; $M$ là kích thước phân tử plasmid tái tổ hợp tính bằng bp; $660\text{ Da/bp}$ là trọng lượng phân tử trung bình của một cặp bazơ).
- Lập chuỗi pha loãng bậc 10 (từ $10^7$ đến $10^2$ bản sao/$\mu$l) để thiết lập đường chuẩn tuyến tính có hệ số tương quan $R^2 > 0,99$ và hiệu suất phản ứng ($E$) đạt từ 95% đến 105%. Phản ứng Real-time PCR chạy với hóa chất qPCRBIO SyGreen Mix trên hệ thống MyGo Pro Real-Time PCR (IT-IS Life Science Ltd, Cork, Ireland).
==========================================================================================
QUY TRÌNH THỰC NGHIỆM ĐỒNG BỘ VÀ NGHIÊM NGẶT
==========================================================================================
Bệnh phẩm (102 Cặp Mô U/Liền kề + 20 Lành + 65 Máu)
│
├─► Tách chiết DNA (QIAamp Kit) ──► Kiểm định NanoDrop 2000c (A260/A280 = 1.8-2.0)
│
├─► PCR Nhân bản Gen mục tiêu (ND1, ND4, HBB, ATP6, ND3, tARN)
│ │
│ ├─► Tinh sạch (ExoSAP-IT / Gel Extraction)
│ ├─► PCR-RFLP (FastDigest: Hin6I, FspBI, DdeI, NlaIII) ──► Điện di Polyacrylamide
│ └─► Sanger Sequencing (ABI PRISM 3100) ──► Đối chiếu BLAST (NC_012920)
│
└─► Real-time PCR Tuyệt đối (qPCRBIO SyGreen trên MyGo Pro System)
│
├── Chuẩn nội: Nhân dòng pJET1.2/blunt mang đoạn chèn ND1, ND4, HBB vào E.coli DH5a
├── Thiết lập đường chuẩn pha loãng bậc 10 (R² > 0.99, Efficiency 95-105%)
└── Định lượng:
+ Số bản sao mtDNA = (Số bản sao ND1 / Số bản sao HBB) * 2
+ Mức độ mất đoạn lớn (%) = [1 - (Số bản sao ND4 / Số bản sao ND1)] * 100%
==========================================================================================
Data và phân tích
Số liệu thực nghiệm được xử lý trên phần mềm thống kê sinh học chuyên sâu (SPSS phiên bản 20.0 và MedCalc). Mức ý nghĩa thống kê được xác lập tại $p < 0,05$. Các phép kiểm thống kê bao gồm:
- So sánh biến định lượng giữa hai nhóm độc lập bằng kiểm định Mann-Whitney U (cho phân phối không chuẩn) và kiểm định Student's t-test (cho phân phối chuẩn). So sánh cặp ghép mô u và mô liền kề bằng kiểm định phi tham số Wilcoxon signed-rank test.
- Phân tích bảng tiếp liên giữa các biến định tính bằng kiểm định $\chi^2$ (Chi-square test) hoặc kiểm định chính xác của Fisher (Fisher's exact test) khi tần số kỳ vọng nhỏ hơn 5.
- Phân tích tương quan đơn biến và đa biến bằng mô hình hồi quy Logistic nhị phân (Binary Logistic Regression) để ước lượng tỷ số Odds Ratio (OR) kèm khoảng tin cậy 95% (95% CI).
- Xây dựng đường cong đặc trưng hoạt động máy thu (Receiver Operating Characteristic - ROC), xác định diện tích dưới đường cong (Area Under Curve - AUC), độ nhạy (Sensitivity), độ đặc hiệu (Specificity), giá trị dự báo dương tính (PPV) và âm tính (NPV) nhằm so sánh hiệu năng giữa mô hình đơn biến và mô hình đa biến tích hợp.
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Sự suy giảm sâu sắc số bản sao mtDNA trong mô ung thư và mối tương quan bệnh học:
- Kết quả Real-time PCR cho thấy số bản sao mtDNA trung bình trong mô u ung thư vú giảm có ý nghĩa thống kê ($p < 0,001$) so với mô lành liền kề và mô tuyến vú lành tính.
- Đáng chú ý, sự suy giảm số bản sao tương quan nghịch chặt chẽ với mức độ ác tính lâm sàng: nhóm khối u có kích thước lớn ($T_{3-4}$) và giai đoạn tiến triển (giai đoạn III–IV) có số lượng bản sao mtDNA thấp hơn rõ rệt so với nhóm giai đoạn sớm ($T_{1-2}$, giai đoạn 0–II, $p = 0,028$). Điều này chứng minh quá trình mất dần bản sao mtDNA đồng hành với sự tăng tốc phân chia tế bào bất thường và mất biệt hóa mô học.
-
Động học nghịch lý của mất đoạn $\Delta\text{mtDNA}^{4977}$ – Bằng chứng của quá trình chọn lọc vô tính:
- Tần suất xuất hiện mất đoạn $\Delta\text{mtDNA}^{4977}$ ở mô lành liền kề (44,1%) cao hơn gấp nhiều lần so với mô u ung thư vú thực thể (12,7%, $p < 0,001$). Mức độ mất đoạn định lượng tuyệt đối cũng ghi nhận cao hơn ở mô liền kề.
- Hiện tượng tưởng chừng nghịch lý này được giải thích về mặt cơ chế phân tử: mô liền kề chịu tác động kéo dài của stress oxy hóa vi môi trường dẫn đến tích lũy đột biến mất đoạn lớn nhưng tế bào chưa phân chia nhanh; trong khi tại khối u ác tính, các phân tử mtDNA mất đoạn 4977 bp (bị mất 7 tiểu phần OXPHOS và 5 tARN) bị đào thải chọn lọc hoặc pha loãng bởi áp lực nhân đôi nhanh chóng của các tế bào phân chia vô tính nhằm duy trì mức hô hấp tối thiểu cần thiết cho sự sống của tế bào u.
- Luận án cũng đã phát hiện các mất đoạn lớn mới chưa từng được công bố quốc tế tại vùng cung mất đoạn chính của hệ gen ty thể ở bệnh nhân Việt Nam.
-
Phát hiện đột biến và đa hình trên gen ATP6 và tác động cấu trúc:
- Phân tích giải trình tự phát hiện đột biến thay thế nucleotide G9053A trên gen ATP6 (dẫn đến thay thế axit amin tương ứng trong tiểu phần tổng hợp ATP synthase F0).
- Tần suất mang alen đột biến G9053A ở nhóm bệnh nhân ung thư vú cao hơn có ý nghĩa thống kê so với nhóm người bình thường đối chứng ($p = 0,014$, $\text{OR} = 2,84$, $95%\text{ CI}: 1,21 - 6,65$), định danh đây là một biến thể chỉ thị nguy cơ cao trong quần thể người Việt.
-
Đột biến điểm gen tARN và biến đổi cấu trúc bậc hai:
- Đột biến điểm A5536T trên gen tARN-Trp được phát hiện và sàng lọc thành công bằng kỹ thuật PCR-RFLP (FspBI).
- Phân tích mô hình hóa tin sinh học cấu trúc phân tử (sử dụng thuật toán dự đoán năng lượng tự do tối thiểu MFE) chứng minh đột biến A5536T làm phá vỡ cặp liên kết bazơ bổ sung trên nhánh thân chấp nhận của tARN-Trp, làm tăng mức năng lượng tự do từ $-8,70\text{ kcal/mol}$ lên $-5,20\text{ kcal/mol}$, dẫn đến mất ổn định cấu trúc không gian ba chiều dạng lá chẽ ba (cloverleaf), cản trở trực tiếp phức hợp Aminoacyl-tRNA synthetase gắn kết tryptophan trong quá trình dịch mã ty thể.
-
Đa hình ND3 A10398G và ND1 A4164G:
- Đa hình ND3 A10398G phân bố với tỷ lệ alen 10398G chiếm ưu thế ở bệnh nhân ung thư vú người Việt, tương đồng với các phát hiện trên quần thể người Ba Lan (Czarnecka và cs., 2010a) và Malaysia (Nadiah và cs., 2012), nhưng khác biệt với quần thể người Mỹ da đen (Canter và cs., 2005).
- Biến đổi ND1 A4164G và T4216C ghi nhận sự gia tăng tần suất ở nhóm bệnh nhân có hạch di căn ($N_{>0}$), đóng vai trò chỉ thị tiềm năng cho khả năng xâm lấn hạch vùng.
==========================================================================================
BẢNG TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN
==========================================================================================
Dấu ấn phân tử Mô u (Tumor) Mô liền kề Máu ĐC/Lành Giá trị P / OR (95% CI)
──────────────────────────────────────────────────────────────────────────────────────────
Số bản sao mtDNA Giảm sâu Bình thường Bình thường p < 0,001 (Nghịch biến với T, Stage)
Mất đoạn 4977 bp 12,7% 44,1% 0% - 5% p < 0,001 (Field cancerization)
ATP6 G9053A Cao Trung bình Thấp p = 0,014 / OR = 2,84 (1,21-6,65)
tARN-Trp A5536T Xuất hiện Không Không Tăng MFE (-8,7 -> -5,2 kcal/mol)
ND3 A10398G (Alen G) Ưu thế - Thấp p = 0,032 (Tương quan nguy cơ)
Mô hình Đa biến AUC = 0,868 (Độ nhạy: 82,4%, Độ đặc hiệu: 85,0%) Vượt trội mô hình đơn biến
==========================================================================================
Implications đa chiều
- Về mặt lý luận khoa học: Luận án cung cấp chứng cứ sinh học phân tử khẳng định rằng tái cấu trúc hệ gen ty thể không phải là hiện tượng thứ cấp ngẫu nhiên mà là một cấu phần thúc đẩy tiến trình sinh ung thư vú. Công trình đã liên kết thành công cơ chế suy giảm sao chép mtDNA với hiệu ứng Warburg và quá trình chọn lọc vô tính trong mô đặc.
- Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa định lượng tuyệt đối bản sao mtDNA và mức độ mất đoạn bằng kỹ thuật Real-time PCR dựa trên thang chuẩn plasmid tái tổ hợp (pJET1.2/blunt chứa các đoạn gen chuẩn ND1, ND4, HBB). Quy trình này có thể chuyển giao và áp dụng trực tiếp cho các nghiên cứu trên các bệnh lý ung thư biểu mô khác (ung thư đại trực tràng, ung thư gan, ung thư dạ dày).
- Về mặt ứng dụng lâm sàng và y học dịch chuyển:
- Đề xuất bộ dấu ấn sinh học phân tử kết hợp giữa đột biến điểm (ATP6 G9053A, ND3 A10398G) và chỉ số số lượng bản sao mtDNA phục vụ sàng lọc phát hiện sớm nguy cơ ung thư vú từ sinh thiết mô hoặc định hướng phát triển sinh thiết lỏng (liquid biopsy).
- Tỷ lệ mất đoạn cao ở mô lành liền kề cung cấp cảnh báo giải phẫu bệnh học quan trọng: mép cắt phẫu thuật cần được đánh giá ở cấp độ biến đổi phân tử nhằm tránh tái phát khối u tại chỗ do các tế bào mang tổn thương tiền ác tính còn sót lại trong vi môi trường mô.
Limitations và Future Research
Mặc dù đạt được nhiều phát hiện có tính bước ngoặt, công trình vẫn tồn tại một số hạn chế khách quan cần được hoàn thiện:
- Thiết kế nghiên cứu cắt ngang hồi cứu: Cỡ mẫu 102 bệnh nhân tuy đủ độ mạnh thống kê cho các phép phân tích đơn biến và đa biến nhưng được thu thập từ một trung tâm điều trị đơn lẻ (Bệnh viện K). Cần triển khai các nghiên cứu đoàn hệ tiến cứu đa trung tâm (Multicenter Prospective Cohort Studies) với cỡ mẫu lớn hơn để kiểm chứng độ ổn định của mô hình.
- Thiếu vắng dữ liệu theo dõi dọc (Longitudinal Follow-up): Luận án chưa theo dõi được tỷ lệ sống thêm toàn bộ (Overall Survival - OS) và tỷ lệ sống thêm không bệnh (Disease-Free Survival - DFS) sau 5 năm của bệnh nhân để xác lập giá trị tiên lượng sinh tồn độc lập của số bản sao mtDNA và mất đoạn $\Delta\text{mtDNA}^{4977}$.
- Phân tích mô tổng thể (Bulk Tissue Analysis): DNA được tách chiết từ toàn bộ khối mô sinh thiết, chưa thực hiện kỹ thuật vi phẫu cắt chọn lọc bằng laser (Laser Capture Microdissection - LCM) hoặc giải trình tự đơn tế bào (Single-cell RNA/DNA Sequencing) để tách biệt hoàn toàn tế bào biểu mô ác tính với tế bào chất đệm, tế bào miễn dịch thâm nhiễm và mạch máu.
- Chưa đo lường trực tiếp hoạt tính chuyển hóa sinh hóa: Nghiên cứu tập trung vào biến đổi vật liệu di truyền mà chưa kết hợp đo trực tiếp tốc độ tiêu thụ oxy (Oxygen Consumption Rate - OCR) và tốc độ axit hóa ngoài tế bào (Extracellular Acidification Rate - ECAR) bằng hệ thống phân tích chuyển hóa tế bào sống (như Seahorse XF Analyzer).
Chương trình nghiên cứu tương lai:
- Hướng 1: Ứng dụng công nghệ giải trình tự thế hệ mới toàn bộ hệ gen ty thể (Next-Generation Sequencing - mtDNA NGS) với độ sâu bao phủ cao ($>10.000\text{X}$) nhằm phát hiện toàn diện các đột biến soma dị tế bào chất ở tần số siêu thấp ($<1%$).
- Hướng 2: Nghiên cứu động học số bản sao mtDNA và đột biến mất đoạn trong DNA ty thể tự do lưu hành trong huyết tương (circulating cell-free mtDNA - ccf-mtDNA) nhằm hoàn thiện công cụ sinh thiết lỏng không xâm lấn cho sàng lọc cộng đồng.
- Hướng 3: Thực nghiệm chức năng in vitro và in vivo: sử dụng công nghệ chuyển nhân tế bào (Transmitochondrial Cybrids) đưa ty thể mang đột biến ATP6 G9053A hoặc tARN A5536T vào dòng tế bào ung thư vú cạn kiệt mtDNA ($\rho^0$) để chứng minh mối quan hệ nhân quả trực tiếp đối với tốc độ tăng sinh và kháng hóa chất điều trị.
- Hướng 4: Mở rộng đánh giá tương tác giữa đột biến gen nhân (BRCA1, BRCA2, TP53, PTEN) và biến đổi hệ gen ty thể trong điều hòa trục truyền tín hiệu nhân – ty thể (Mitochondrial-Nuclear Retrograde Signaling).
Tác động và ảnh hưởng
- Tác động học thuật: Công trình là tài liệu tham khảo nền tảng cho chuyên ngành Nhân chủng học, Di truyền học người và Ung thư học phân tử tại Việt Nam. Các bài báo khoa học công bố từ luận án trên các tạp chí chuyên ngành trong nước và quốc tế đã bổ sung mảnh ghép quan trọng về đa dạng di truyền ty thể của chủng tộc người Việt (Kinh) vào cơ sở dữ liệu hệ gen người toàn cầu (MITOMAP, GenBank).
- Chuyển giao công nghệ y học: Bộ quy trình Real-time PCR định lượng chuẩn tuyệt đối plasmid và hệ thống mồi - enzyme giới hạn PCR-RFLP của luận án đã được chuẩn hóa, sẵn sàng chuyển giao cho các phòng xét nghiệm sinh học phân tử tại các bệnh viện ung bướu tuyến trung ương phục vụ công tác nghiên cứu dịch chuyển.
- Định hướng chính sách y tế cộng đồng: Kết quả chỉ ra sự bộc phát ung thư vú ở độ tuổi trẻ tại Việt Nam (trung bình 50 tuổi) cùng với gánh nặng dương tính giả cao của các phương pháp sàng lọc thông thường, cung cấp cơ sở dữ liệu thực nghiệm giá trị cho Bộ Y tế trong việc xây dựng chiến lược sàng lọc phân tầng nguy cơ kết hợp dấu ấn phân tử cho phụ nữ từ 40 tuổi trở lên.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận nghiên cứu hệ gen ty thể chuẩn mực, từ thiết kế mồi, giải trình tự, mô hình hóa cấu trúc RNA đến phân tích toán học đường cong ROC.
- Bác sĩ Giải phẫu bệnh và Ung bướu lâm sàng: Có thêm góc nhìn phân tử về "vùng tổn thương tiền ác tính" tại rìa diện cắt phẫu thuật, hỗ trợ quyết định phạm vi phẫu thuật triệt căn bảo tồn tuyến vú chính xác hơn.
- Các công ty Công nghệ sinh học & Dược phẩm: Khai thác các trình tự gen đột biến (ATP6, ND3, tARN) làm mục tiêu sinh học (biomarker targets) để phát triển các kit chẩn đoán phân tử multiplex PCR hoặc định hướng các hoạt chất tái lập trình chuyển hóa ty thể nhằm tiêu diệt tế bào gốc ung thư (Cancer Stem Cells).
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và nó mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc chứng minh sự tích hợp giữa Thuyết chuyển hóa Warburg và Thuyết ngưỡng dị tế bào chất Wallace trong môi trường vi mô mô ung thư vú người Việt. Luận án chỉ ra rằng sự suy giảm số lượng bản sao mtDNA đi kèm với các đột biến làm bất hoạt phức hệ I (ND1, ND3) và phức hệ V (ATP6) không chỉ làm suy yếu con đường phosphoryl hóa oxy hóa (OXPHOS) mà còn tạo áp lực chọn lọc vô tính đào thải các phân tử mất đoạn lớn $\Delta\text{mtDNA}^{4977}$ trong khối u nhằm duy trì mức sao chép tối thiểu phục vụ quá trình phân chia nhanh chóng của tế bào ác tính.
2. Đổi mới phương pháp luận của luận án vượt trội so với các nghiên cứu trước đây như thế nào?
So với các nghiên cứu của Đái Duy Ban và cs. (2003) vốn chỉ dùng PCR định tính trên 2 mẫu, hoặc các nghiên cứu quốc tế dùng phương pháp $\Delta\Delta Ct$ bán định lượng dễ sai số, luận án đã:
- Tự thiết kế và tạo dòng thành công 3 vector plasmid tái tổ hợp pJET1.2/blunt chứa đoạn gen ND1, ND4 và gen nhân HBB.
- Thiết lập đường chuẩn tuyệt đối pha loãng bậc 10 đạt hiệu suất $95 - 105%$ và $R^2 > 0,99$, cho phép định lượng chính xác số bản sao mtDNA trên mỗi tế bào đơn vị và tính toán tỷ lệ phần trăm mất đoạn lớn trực tiếp từ tỷ số khuếch đại ND4/ND1.
- Kết hợp hoàn hảo giữa kỹ thuật sàng lọc nhanh PCR-RFLP (tối ưu hóa với 4 enzyme FastDigest: Hin6I, FspBI, DdeI, NlaIII), giải trình tự Sanger và mô hình hóa năng lượng tự do MFE cấu trúc bậc 2 của RNA.
3. Phát hiện bất ngờ nhất (counter-intuitive result) thu được từ dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là tỷ lệ và mức độ mất đoạn lớn $\Delta\text{mtDNA}^{4977}$ ở mô u ung thư vú (12,7%) thấp hơn có ý nghĩa thống kê so với mô lành liền kề cách mép u 5 cm (44,1%, $p < 0,001$). Về mặt trực giác ban đầu, mô ác tính thường được dự đoán chứa nhiều tổn thương DNA hơn. Tuy nhiên, kết quả chứng minh mô liền kề là nơi gánh chịu tổn thương oxy hóa mạn tính tích lũy (Field Cancerization), trong khi tế bào ung thư trong khối u có tốc độ phân chia cực cao đã kích hoạt cơ chế chọn lọc thanh lọc (purifying selection), loại bỏ các phân tử mtDNA mang mất đoạn quá lớn (mất tới 4977 bp) để ưu tiên giữ lại các phân tử mtDNA nguyên vẹn hơn phục vụ quá trình nhân đôi DNA ty thể.
4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?
Hoàn toàn đầy đủ. Luận án công bố chi tiết:
- Trình tự 10 cặp mồi xuôi và ngược, kích thước khuếch đại chính xác từng cặp bazơ, vị trí tọa độ nucleotide theo chuẩn tham chiếu NCBI NC_012920.
- Nồng độ hóa chất thành phần cho phản ứng PCR ($6,25,\mu\text{l}$ Maxima Mastermix 2X, $0,25,\mu\text{l}$ mồi $0,2,\mu\text{M}$, khuôn DNA $1 - 2,5,\text{ng}/\mu\text{l}$) và chu trình nhiệt biến tính/gắn mồi/kéo dài chi tiết đến từng giây.
- Điều kiện phản ứng cắt enzyme giới hạn ($10,\mu\text{l}$ PCR, $2,\mu\text{l}$ FastDigest buffer 10X, $1,\mu\text{l}$ enzyme, ủ $37^\circ\text{C}$ trong 5–10 phút, bất hoạt nhiệt $65 - 80^\circ\text{C}$).
- Công thức toán học chuyển đổi nồng độ đo quang sang số lượng bản sao plasmid và phương trình hồi quy tuyến tính Real-time PCR.
5. Chương trình nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào?
Chương trình 10 năm gồm 3 giai đoạn chiến lược:
- Giai đoạn 1 (Năm 1–3): Mở rộng kiểm chứng đoàn hệ tiến cứu đa trung tâm trên 1.000 bệnh nhân ung thư vú tại miền Bắc, miền Trung và miền Nam; giải trình tự toàn bộ hệ gen ty thể (Deep mtDNA NGS).
- Giai đoạn 2 (Năm 4–6): Tối ưu hóa kỹ thuật Digital PCR (dPCR) để định lượng ccf-mtDNA và các đột biến ATP6, ND3, tARN trong huyết tương bệnh nhân, phát triển bộ kit chẩn đoán sinh thiết lỏng sàng lọc sớm không xâm lấn.
- Giai đoạn 3 (Năm 7–10): Thực nghiệm tế bào Cybrid và mô hình chuột PDX (Patient-Derived Xenograft) mang đột biến ty thể đặc thù người Việt để thử nghiệm các phân tử thuốc nhắm trúng đích chuyển hóa ty thể (Mitochondrial-targeted therapeutics).
Kết luận
Công trình nghiên cứu luận án tiến sĩ của NCS. Nguyễn Thị Tú Linh đã đạt được 5 đóng góp học thuật và thực tiễn cốt lõi:
- Thiết lập cơ sở dữ liệu hệ thống đầu tiên: Cung cấp bộ dữ liệu hoàn chỉnh, có tính hệ thống cao nhất tại Việt Nam về biến đổi số bản sao, mất đoạn lớn và đột biến điểm các gen ty thể (ATP6, ND1, ND3, tARN) trên 102 cặp mô ung thư vú và các nhóm đối chứng chuẩn tắc.
- Làm sáng tỏ cơ chế sinh học phân tử: Xác nhận hiện tượng suy giảm sâu sắc số bản sao mtDNA trong mô u tỷ lệ nghịch với giai đoạn TNM và kích thước khối u; đồng thời giải mã hiện tượng nghịch lý của mất đoạn $\Delta\text{mtDNA}^{4977}$ qua cơ chế chọn lọc vô tính tế bào và hiệu ứng tổn thương vùng (Field Cancerization).
- Định danh các chỉ thị phân tử mới: Phát hiện đột biến ATP6 G9053A ($\text{OR} = 2,84$), đa hình ND3 A10398G, đột biến tARN-Trp A5536T (làm tăng năng lượng tự do MFE và phá vỡ cấu trúc không gian tARN) và các mất đoạn lớn mới tại vùng Major Arc ở bệnh nhân Việt Nam.
- Đổi mới công nghệ định lượng phân tử: Chuẩn hóa thành công hệ thống đường chuẩn plasmid tuyệt đối cho Real-time PCR định lượng số bản sao và tỷ lệ mất đoạn mtDNA với độ tin cậy và khả năng tái lặp cao.
- Xây dựng mô hình toán học dự báo nguy cơ: Thiết lập thành công mô hình hồi quy đa biến tích hợp số bản sao và mức độ mất đoạn lớn, đạt diện tích dưới đường cong $\text{AUC} = 0,868$, mở ra hướng tiếp cận định lượng mới trong sàng lọc phân tầng nguy cơ và đánh giá diện cắt phẫu thuật ung thư vú.
Công trình tạo tiền đề vững chắc mở ra các nhánh nghiên cứu chuyên sâu về sinh thiết lỏng ty thể và điều trị trúng đích chuyển hóa ung thư, khẳng định bước tiến quan trọng của khoa học y sinh và nhân chủng học phân tử Việt Nam trên bản đồ học thuật quốc tế.
