Bối cảnh và vấn đề nghiên cứu

Theo số liệu thống kê của Tổ chức Y tế Thế giới (WHO, 2018), ung thư vú là bệnh lý ác tính đứng thứ hai trong nhóm các bệnh ung thư gây tử vong hàng đầu ở phụ nữ trên phạm vi toàn cầu. Hàng năm thế giới ghi nhận khoảng 2,1 triệu ca mắc mới và khoảng 670.000 trường hợp tử vong do căn bệnh này (Ghonche và cộng sự, 2016; WHO, 2018). Riêng tại Việt Nam, số liệu năm 2018 của WHO ghi nhận có 15.229 trường hợp ung thư vú mới được phát hiện. Sự hình thành khối u tuyến vú bắt nguồn từ sự tăng sinh mất kiểm soát của tế bào biểu mô ống dẫn sữa hoặc thùy tuyến vú, chịu sự chi phối của các đột biến gen kết hợp các biến đổi biểu sinh (epigenetics). Trong đó, hiện tượng tăng methyl hóa bất thường (hypermethylation) tại vùng promoter của các gen ức chế khối u (tumor suppressor genes) dẫn đến bất hoạt phiên mã gen đang được ghi nhận là một dấu chứng sinh học phân tử tiềm năng cho tiên lượng và chẩn đoán sớm.

Hai gen giữ vai trò kiểm soát chu kỳ tế bào và sửa chữa tổn thương di truyền được tập trung khảo sát là BRCA1 (nằm trên nhiễm sắc thể 17q21) và p16INK4a (còn gọi là CDKN2A, nằm trên nhiễm sắc thể 9p21). Nhiều nghiên cứu độc lập trên thế giới đã ghi nhận tần số methyl hóa bất thường trên hai gen này dao động trong biên độ rất lớn tùy thuộc vào chủng tộc, quần thể bệnh nhân và loại mẫu xét nghiệm (ví dụ tần số methyl hóa BRCA1 dao động từ 24% đến 68%, trên p16INK4a dao động từ 24,16% đến 40%). Xuất phát từ tính biến thiên đó cùng nhu cầu kế thừa hướng nghiên cứu về dấu chứng sinh học phân tử của Khoa Công nghệ Sinh học - Trường Đại học Mở TP.HCM, đề tài khóa luận "Khảo sát tính chất methyl hóa bất thường trên gen BRCA1 và p16INK4a trên bộ mẫu ung thư vú ở Việt Nam: Phân tích tổng hợp & khảo sát thực nghiệm" được thực hiện.

Mục tiêu nghiên cứu của khóa luận là xác định sự tương quan giữa tính chất methyl hóa bất thường trên vùng promoter gen BRCA1 và p16INK4a ở bệnh ung thư vú trên thế giới và tại Việt Nam, tập trung vào nội dung phân tích tổng hợp và khảo sát thực nghiệm bằng phương pháp MSP với một số mẫu mô vú của người Việt Nam.

Khóa luận triển khai 3 nội dung nghiên cứu chính:

  1. Phân tích tổng hợp: Xác định sự tương quan giữa tính chất methyl hóa vượt mức trên vùng promoter gen BRCA1 và p16INK4a và bệnh ung thư vú, dựa vào phương pháp xác định chỉ số tỷ suất chênh (Odds Ratio - OR).
  2. Khảo sát trên máy tính: Cập nhật thông tin dữ liệu khảo sát bộ mồi phù hợp với quy trình Nested-MSP xác định tính chất methyl hóa trên vùng promoter gen BRCA1 và p16INK4a đối với một số mẫu bệnh phẩm ung thư vú.
  3. Khảo sát thực nghiệm: Thực hiện quy trình Nested-MSP để xác định tính chất methyl hóa trên vùng promoter gen BRCA1 và p16INK4a đối với một số mẫu bệnh phẩm (mẫu mô vú) của người Việt Nam.

Đối tượng nghiên cứu bao gồm: cơ sở dữ liệu các công trình nghiên cứu ca - chứng công bố quốc tế về tính chất methyl hóa gen BRCA1 và p16INK4a; trình tự nucleotide promoter và các bộ mồi PCR đặc hiệu; các mẫu bệnh phẩm mô đúc khối paraffin của bệnh nhân ung thư vú tại Việt Nam. Phạm vi nghiên cứu về mặt thời gian và không gian: dữ liệu phân tích tổng hợp được thu thập, sàng lọc cập nhật đến tháng 4 năm 2019 từ hệ thống cơ sở dữ liệu PubMed/NCBI và Google Scholar; nội dung khảo sát tin sinh học và thực nghiệm được tiến hành tại phòng thí nghiệm chuyên ngành Công nghệ Sinh học Y Dược, Trường Đại học Mở TP.HCM, nghiệm thu vào tháng 5 năm 2019.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của nghiên cứu dựa trên các nguyên lý sinh học phân tử và di truyền học biểu sinh (epigenetics). Hiện tượng methyl hóa DNA là phản ứng gắn nhóm methyl (-CH3) vào vị trí carbon số 5 của vòng cytosine tạo thành 5-methylcytosine (5mC), được xúc tác bởi enzyme DNA methyltransferases (DNMTs) với cơ chất cho nhóm methyl là S-adenosylmethionine (SAM). Vùng tập trung các cặp dinucleotide CpG được gọi là đảo CpG (CpG islands - CGIs), thường có chiều dài từ 300 đến 3.000 cặp base (bp), với tỷ lệ GC trên 50% và tỷ lệ quan sát/kỳ vọng CpG/GpC lớn hơn 0,6. Khi các đảo CpG tại vùng promoter của các gen đè nén u như BRCA1 (tham gia sửa chữa DNA theo cơ chế tái tổ hợp tương đồng) và p16INK4a (ức chế phức hợp kinase phụ thuộc cyclin CDK4/6, duy trì trạng thái khử phosphoryl của protein Rb nhằm dừng chu kỳ tế bào ở pha G1) bị tăng methyl hóa, sự gắn kết của các yếu tố phiên mã bị cản trở, dẫn đến khóa phiên mã gen và thúc đẩy chuyển dạng ác tính.

Phương pháp nghiên cứu được chia làm ba hợp phần:

  1. Phương pháp phân tích tổng hợp (Meta-analysis): Tìm kiếm các bài báo khoa học theo nguyên tắc PICO trên PubMed, PubMed Central, Google Scholar; trích xuất các biến số định lượng gồm cỡ mẫu bệnh, cỡ mẫu lành, số ca có methyl hóa, phương pháp xét nghiệm, chủng tộc và đặc điểm giải phẫu bệnh; tính toán chỉ số tỷ suất chênh (OR) và chỉ số nguy cơ tương đối (RR) kèm khoảng tin cậy 95% (95% CI); đánh giá tính bất đồng nhất giữa các nghiên cứu bằng chỉ số Higgins $I^2$. Áp dụng mô hình ảnh hưởng ngẫu nhiên (Random effects model) khi $I^2$ cao và mô hình ảnh hưởng cố định (Fixed effects model) khi $I^2$ thấp; kiểm tra sai số xuất bản (publication bias) bằng đồ thị hình phễu Funnel plot trên phần mềm MedCalc v18.

  2. Phương pháp khảo sát trên máy tính (In silico): Thu thập trình tự chuẩn của gen BRCA1 (mã số truy cập NC_000017.11) và gen p16INK4a (mã số truy cập NC_000009.12) từ cơ sở dữ liệu NCBI; định vị các đảo CpG bằng công cụ trực tuyến MethPrimer; xác định các vị trí gắn yếu tố điều hòa phiên mã bằng công cụ LASAGNA; tính toán các thông số vật lý của bộ mồi (nhiệt độ nóng chảy $T_m$, tỷ lệ %GC, mức năng lượng tự do $\Delta G$ cho khả năng tạo kẹp tóc hairpin, self-dimer, hetero-dimer) bằng IDT OligoAnalyzer và kiểm tra sự bắt cặp đặc hiệu bằng phần mềm Annhyb.

  3. Phương pháp khảo sát thực nghiệm trong phòng thí nghiệm:

  • Tách chiết DNA tổng số từ mẫu mô đúc paraffin: Khử sáp paraffin bằng dung dịch Xylene (ly tâm 13.000 vòng/phút, 2 phút), hoàn nguyên mô qua chuỗi cồn ethanol giảm dần (100%, 70%, 50%); ly giải tế bào và màng nhân bằng hỗn hợp SDS 10% cùng Proteinase K (20 ng/ml) ủ ở 56°C trong 24-48 giờ; tinh sạch protein bằng hỗn hợp dung môi phenol:chloroform (tỷ lệ 1:1) và chloroform; kết tủa DNA bằng ammonium acetate (NH4OAc) 5M và ethanol tuyệt đối ở 4°C, rửa tủa bằng ethanol 70% và hòa tan trong nước cất hai lần vô trùng.
  • Đo quang phổ kiểm tra độ tinh sạch và nồng độ DNA: Đo độ hấp thụ quang học ở bước sóng 260 nm ($A_{260}$) và 280 nm ($A_{280}$), đánh giá độ tinh sạch dựa trên tỷ số $A_{260}/A_{280}$ đạt ngưỡng 1,8 - 2,0.
  • Xử lý biến đổi DNA bằng sodium bisulfite: Sử dụng bộ sinh phẩm thương mại EZ DNA Methylation-Gold Kit nhằm chuyển hóa cytosine không methyl hóa thành uracil, trong khi cytosine đã methyl hóa (5mC) được bảo toàn nguyên vẹn.
  • Phản ứng Nested-MSP (Nested Methylation-Specific PCR): Tiến hành qua hai vòng nhiệt. Vòng 1 sử dụng cặp mồi ngoài khuếch đại bao quát vùng promoter chứa đảo CpG bất kể trạng thái methyl hóa. Vòng 2 sử dụng sản phẩm vòng 1 làm khuôn, chia thành hai phản ứng song song: phản ứng M sử dụng cặp mồi đặc hiệu cho allele bị methyl hóa, và phản ứng U sử dụng cặp mồi đặc hiệu cho allele không bị methyl hóa.
  • Điện di kiểm tra sản phẩm PCR: Chạy điện di trên bản gel agarose 1,5% trong dung dịch đệm TAE 1X ở hiệu điện thế 70V trong 30 phút, phát hiện các băng DNA dưới tia cực tím (UV) bằng chất nhuộm GelRed trên hệ thống Gel Doc Bio-Rad.

Nguồn dữ liệu của nghiên cứu bao gồm: dữ liệu thứ cấp từ 18 nghiên cứu ca - chứng về BRCA1 và 11 nghiên cứu ca - chứng về p16INK4a công bố trên thế giới; dữ liệu sơ cấp từ các mẫu mô ung thư vú đúc paraffin thu thập tại Việt Nam, kế thừa nguồn mẫu và thông số từ các công trình công bố trước đó của nhóm nghiên cứu ThS. Trương Kim Phượng, PGS.TS Lê Huyền Ái Thúy (2014, 2017) và dữ liệu khóa luận của sinh viên Bùi Thị Mỹ Trang.

Nội dung chính theo từng chương

Chương I: Tổng quan nghiên cứu

Chương này hệ thống hóa các kiến thức bệnh học, dịch tễ học và cơ chế sinh học phân tử liên quan đến ung thư vú. Tác giả trình bày chi tiết phân loại bệnh học gồm ung thư không xâm lấn, ung thư xâm lấn, ung thư biểu mô ống dẫn sữa tại chỗ (DCIS), ung thư tiểu thùy tại chỗ (LCIS), ung thư ống dẫn sữa xâm lấn (IDC), ung thư biểu mô dạng tủy, dạng nhầy, dạng ống, ung thư vú dạng viêm, bệnh Paget núm vú và bướu diệp thể. Các bệnh lý tuyến vú lành tính có nguy cơ chuyển dạng ác tính như xơ nang, u nhú trong ống dẫn sữa và hoại tử mỡ cũng được phân tích.

Tiếp theo, chương luận giải cơ chế di truyền biểu sinh, sự hình thành 5-methylcytosine tại các dinucleotide CpG dưới tác dụng của enzyme DNMTs và nguồn gốc nhóm methyl từ SAM. Tác giả giới thiệu cấu trúc gen BRCA1 tại vị trí 17q21 tham gia cơ chế sửa chữa đứt gãy DNA sợi đôi, và gen p16INK4a tại vị trí 9p21 điều hòa điểm kiểm soát chu kỳ tế bào G1/S thông qua ức chế CDK4/6 ngăn cản bất hoạt protein pRb. Tổng quan ghi nhận các công bố quốc tế về tần số methyl hóa của hai gen này: Wu và cs. (2016) ghi nhận tần số methyl hóa BRCA1 là 24% trên 70 mẫu tại Trung Quốc; Hassan và cs. (2017), Yadav và cs. (2018) ghi nhận tỷ lệ lần lượt là 55% trên 60 mẫu và 68% trên 75 mẫu tại Ấn Độ. Trên gen p16INK4a, Cavusoglu và cs. (2010) ghi nhận 24,16% trên 62 mẫu ở Thổ Nhĩ Kỳ; Lee và cs. (2012) ghi nhận 34,5% trên 200 mẫu ở Hàn Quốc; Wu và cs. (2016) ghi nhận 40% trên 70 mẫu ở Trung Quốc. Cuối chương, tác giả trình bày nguyên lý toán thống kê trong phân tích tổng hợp (chỉ số OR, RR, $I^2$, mô hình Fixed/Random effects) và nguyên lý kỹ thuật phân tích methyl hóa bằng MSP và Nested-MSP.

Chương II: Vật liệu và phương pháp nghiên cứu

Chương II mô tả chi tiết danh mục trang thiết bị, hóa chất, cơ sở dữ liệu tin sinh học (NCBI, GeneCards, MethPrimer, BLAST, LASAGNA, IDT OligoAnalyzer, phần mềm Annhyb, phần mềm thống kê MedCalc v18) và các quy trình thao tác chuẩn:

  • Quy trình phân tích tổng hợp: Xác lập từ khóa tìm kiếm, sàng lọc bài báo theo sơ đồ dòng, trích xuất dữ liệu ca - chứng, chạy mô hình thống kê, vẽ biểu đồ Forest plot và Funnel plot.
  • Quy trình tách chiết DNA từ mô paraffin: Các bước loại nến bằng Xylene, hoàn nguyên mô bằng dãy ethanol nồng độ giảm dần (100%, 70%, 50%), tiêu hóa mô bằng đệm ly giải và proteinase K (20 ng/ml) ủ ở 56°C từ 24 đến 48 giờ; khử tạp protein bằng dung môi hữu cơ phenol:chloroform (1:1) và chloroform; kết tủa DNA bằng NH4OAc 5M và ethanol tuyệt đối, rửa tủa bằng ethanol 70% và hòa tan trong nước cất vô trùng.
  • Quy trình kiểm tra chất lượng DNA: Đo mật độ quang $A_{260}$ và $A_{280}$, tính nồng độ theo công thức $C (\mu g/ml) = 50 \times A_{260} \times d$.
  • Quy trình chuyển đổi hóa học DNA bằng sodium bisulfite: Chi tiết các bước sử dụng bộ kit EZ DNA Methylation-Gold Kit qua chu trình nhiệt (98°C trong 10 phút, 64°C trong 2,5 giờ và giữ ở 4°C), kết hợp làm sạch qua cột ly tâm Zymo-Spin IC.
  • Quy trình phản ứng Nested-MSP: Thiết lập thành phần phản ứng tổng thể tích 50 $\mu l$ gồm Master mix 2X (25 $\mu l$), mồi xuôi 10 $\mu M$ (0,5 $\mu l$), mồi ngược 10 $\mu M$ (0,5 $\mu l$), khuôn DNA (10 - 100 ng/$\mu l$) và nước cất vô trùng. Mô tả chu trình luân nhiệt cho từng phản ứng khuếch đại phổ quát và phản ứng phân biệt trạng thái methyl hóa (M và U).
  • Quy trình điện di gel agarose 1,5%: Sử dụng đệm TAE 1X, chất nhuộm GelRed, điện thế 70V trong 30 phút, đọc kết quả trên hệ thống soi gel UV Gel Doc.

Chương III: Kết quả khai thác dữ liệu và phân tích tổng hợp

Chương này trình bày các phát hiện định lượng từ dữ liệu tổng hợp:

  • Sàng lọc bài báo: Từ 37 bài báo ban đầu về BRCA1, qua các bước loại trừ chọn lọc được 18 công bố khoa học ca - chứng đạt chuẩn. Từ 30 bài báo ban đầu về p16INK4a, chọn lọc được 11 công bố khoa học ca - chứng đạt chuẩn.
  • Phân tích tương quan trên gen BRCA1: Tổng hợp từ 18 nghiên cứu gồm 1.540 mẫu ung thư vú (76 mẫu máu, 1.464 mẫu mô) và 907 mẫu lành (169 mẫu máu, 738 mẫu mô). Giá trị bất đồng nhất $I^2 = 76,36%$ đối với OR và $I^2 = 72,77%$ đối với RR, áp dụng mô hình ảnh hưởng ngẫu nhiên. Kết quả chỉ số tỷ suất chênh gộp $OR = 8,979$ ($95% CI: 4,359 - 18,497; P < 0,0001$), cho thấy tính chất methyl hóa promoter BRCA1 xuất hiện ở nhóm mô bệnh cao gấp 8,979 lần nhóm lành. Chỉ số nguy cơ gộp $RR = 4,973$ ($95% CI: 2,523 - 9,802; P < 0,0001$), phản ánh nguy cơ mắc ung thư vú cao gấp 4,973 lần khi có hiện tượng tăng methyl hóa gen BRCA1.
  • Phân tích tương quan trên gen p16INK4a: Tổng hợp từ 11 nghiên cứu gồm 1.067 mẫu ung thư vú (134 mẫu máu, 933 mẫu mô) và 761 mẫu lành (263 mẫu máu, 498 mẫu mô). Tính bất đồng nhất $I^2 = 66,74%$ đối với OR và $I^2 = 90,13%$ đối với RR, áp dụng mô hình ảnh hưởng ngẫu nhiên. Kết quả chỉ số tỷ suất chênh gộp $OR = 12,925$ ($95% CI: 4,585 - 36,429; P < 0,0001$), cho thấy độ chênh xuất hiện methyl hóa ở nhóm bệnh cao gấp 12,925 lần nhóm lành. Chỉ số nguy cơ gộp $RR = 8,920$ ($95% CI: 2,203 - 36,101; P = 0,002$), phản ánh nguy cơ mắc ung thư vú cao gấp 8,92 lần khi promoter gen p16INK4a bị methyl hóa.
  • Phân tích phân nhóm theo loại mẫu, phương pháp và chủng tộc:
    • Loại mẫu: Trên gen BRCA1, mẫu mô có $OR = 7,100$ ($95% CI: 3,240 - 15,558$), mẫu máu có $OR = 13,549$ ($95% CI: 1,721 - 106,647$). Trên gen p16INK4a, mẫu mô có $OR = 7,908$ ($95% CI: 2,521 - 24,804$), mẫu máu có $OR = 45,972$ ($95% CI: 15,296 - 138,164; P < 0,001$).
    • Phương pháp phân tích: Trên gen BRCA1, kỹ thuật MSP ghi nhận $OR = 7,996$ ($95% CI: 3,556 - 17,979$), các kỹ thuật khác (MS-MLPA, MS-HRM, Pyrosequencing) ghi nhận $OR = 7,013$. Trên gen p16INK4a, kỹ thuật MSP ghi nhận $OR = 12,129$, phương pháp RT-PCR ghi nhận $OR = 16,411$.
    • Chủng tộc: Trên gen BRCA1, người châu Á có $OR = 8,654$ ($95% CI: 3,810 - 19,657; P < 0,001$), người châu Âu có $OR = 3,616$ ($P = 0,331$), nhóm khác (châu Phi, châu Mỹ) có $OR = 15,228$. Trên gen p16INK4a, người châu Á có $OR = 23,573$ ($95% CI: 13,652 - 40,703; P < 0,001$), trong khi người châu Âu chỉ đạt $OR = 1,254$ ($P = 0,838$).
    • Đặc điểm bệnh học: Sự tương quan giữa methyl hóa và độ tuổi, phân độ mô học, di căn hạch lympho đều có giá trị $P > 0,05$. Riêng trên gen BRCA1, bệnh nhân ở giai đoạn muộn (III và IV) có độ chênh methyl hóa cao gấp 1,543 lần so với giai đoạn sớm (I và II) với mức ý nghĩa thống kê ($OR = 1,543; 95% CI: 1,142 - 2,085; P = 0,005$).

Chương IV: Khảo sát trên máy tính và kết quả thực nghiệm

Chương IV trình bày dữ liệu tin sinh học và các kết quả tối ưu hóa thực nghiệm:

  • Dữ liệu vùng promoter và đảo CpG: Vùng promoter gen BRCA1 được định vị tại vị trí 43.044.295 đến 43.125.483 trên mạch âm nhiễm sắc thể số 17 (NC_000017.11) với chiều dài 1.500 bp; công cụ MethPrimer xác định được 2 đảo CpG có kích thước lần lượt là 211 bp và 200 bp; công cụ LASAGNA xác định các vị trí gắn yếu tố phiên mã CACCC-binding factor và FOXJ2. Vùng promoter gen p16INK4a được định vị trên mạch dương nhiễm sắc thể số 9 (NC_000009.12) với chiều dài 2.000 bp; MethPrimer xác định 4 đảo CpG có kích thước lần lượt là 257 bp, 142 bp, 130 bp và 202 bp, cũng chứa các vị trí gắn yếu tố phiên mã CACCC-binding factor và FOXJ2.
  • Đánh giá bộ mồi Nested-MSP: Tác giả kiểm tra các thông số lý hóa của các cặp mồi:
    • Gen BRCA1: Cặp mồi ngoài vòng 1 (BRCA1_F: 27 bp, $T_m = 56,4^\circ C$; BRCA1_R: 24 bp, $T_m = 54,9^\circ C$) cho kích thước sản phẩm 143 bp. Cặp mồi đặc hiệu allele methyl hóa vòng 2 (BRCA1_MF: 20 bp, $T_m = 58,3^\circ C$; BRCA1_MR: 20 bp, $T_m = 63,0^\circ C$) cho kích thước sản phẩm 68 bp. Cặp mồi đặc hiệu allele không methyl hóa vòng 2 (BRCA1_UF: 20 bp, $T_m = 49,7^\circ C$; BRCA1_UR: 20 bp, $T_m = 52,7^\circ C$) cho kích thước sản phẩm 68 bp. Mức độ tương đồng bắt cặp đạt 100%.
    • Gen p16INK4a: Cặp mồi ngoài vòng 1 (P16_F: 21 bp, $T_m = 56,1^\circ C$; P16_R: 24 bp, $T_m = 57,9^\circ C$) cho kích thước sản phẩm 277 bp. Cặp mồi đặc hiệu allele methyl hóa vòng 2 (P16_MF: 24 bp, $T_m = 62,6^\circ C$; P16_MR: 21 bp, $T_m = 63,3^\circ C$) cho kích thước 150 bp. Cặp mồi đặc hiệu allele không methyl hóa vòng 2 (P16_UF: 24 bp, $T_m = 56,8^\circ C$; P16_UR: 22 bp, $T_m = 55,4^\circ C$) cho kích thước 151 bp.
  • Kết quả khảo sát thực nghiệm: Thực hiện thành công quy trình tách chiết DNA từ mẫu mô paraffin ung thư vú và xử lý chuyển đổi bisulfite. Tác giả vận hành thành công phản ứng Nested-MSP với nhiệt độ lai tối ưu của cặp mồi ngoài khuếch đại promoter BRCA1 là 56°C. Kết quả điện di trên gel agarose 1,5% xác nhận sự xuất hiện rõ nét của các phân đoạn DNA đặc hiệu ở giếng M (allele methyl hóa) và giếng U (allele không methyl hóa), khẳng định tính khả thi của quy trình trên nguồn mẫu bệnh nhân ung thư vú tại Việt Nam.

Kết quả và đóng góp

Khóa luận đã tổng hợp và lượng hóa một cách toàn diện mối liên quan giữa tính chất tăng methyl hóa vùng promoter hai gen ức chế khối u BRCA1 và p16INK4a với nguy cơ mắc bệnh ung thư vú.

Bảng số liệu tổng hợp các chỉ số dịch tễ học phân tử từ phân tích tổng hợp:

Gen mục tiêu Số công trình (N) Cỡ mẫu bệnh Cỡ mẫu lành Tính bất đồng nhất $I^2$ (%) Mô hình phân tích Tỷ suất chênh gộp OR (95% CI) Chỉ số nguy cơ gộp RR (95% CI) Mức ý nghĩa thống kê (P)
BRCA1 18 1.540 907 76,36 Ngẫu nhiên (R) 8,979 [4,359 - 18,497] 4,973 [2,523 - 9,802] < 0,0001
p16INK4a 11 1.067 761 66,74 Ngẫu nhiên (R) 12,925 [4,585 - 36,429] 8,920 [2,203 - 36,101] < 0,0001

Bảng tổng hợp chỉ số tỷ suất chênh (OR) theo phân nhóm chủng tộc và loại mẫu:

Phân nhóm khảo sát Chỉ số OR trên gen BRCA1 (95% CI) Mức ý nghĩa P (BRCA1) Chỉ số OR trên gen p16INK4a (95% CI) Mức ý nghĩa P (p16INK4a)
Chủng tộc châu Á 8,654 [3,810 - 19,657] < 0,001 23,573 [13,652 - 40,703] < 0,001
Chủng tộc châu Âu 3,616 [0,271 - 48,175] 0,331 1,254 [0,285 - 4,728] 0,838
Mẫu mô sinh thiết 7,100 [3,240 - 15,558] < 0,001 7,908 [2,521 - 24,804] < 0,001
Mẫu máu ngoại vi 13,549 [1,721 - 106,647] 0,013 45,972 [15,296 - 138,164] < 0,001

Về mặt đóng góp khoa học so với các công bố trước:

  • So với nghiên cứu của Le Huyen Ai Thuy và cộng sự (2017) trên gen BRCA1 ($OR = 4,312; 95% CI: 2,395 - 7,765$ trên 24 nghiên cứu), tác giả đã cập nhật dữ liệu đến năm 2019 và áp dụng tiêu chí sàng lọc chặt chẽ loại trừ các công bố có khoảng tin cậy dưới 1, ghi nhận độ chênh có ý nghĩa cao hơn ($OR = 8,979$), đồng thuận và củng cố vững chắc nhận định về vai trò của methyl hóa BRCA1.
  • Khóa luận làm rõ đặc tính dịch tễ học biểu sinh mang tính đặc thù chủng tộc: Tần số tăng methyl hóa p16INK4a ở bệnh nhân ung thư vú châu Á có độ chênh lệch cực kỳ rõ rệt so với người lành ($OR = 23,573$), cao vượt trội so với quần thể người châu Âu ($OR = 1,254; P = 0,838$). Điều này chứng minh p16INK4a là dấu ấn sinh học đặc biệt phù hợp cho việc tầm soát ung thư vú ở quần thể bệnh nhân châu Á, bao gồm Việt Nam.
  • Khẳng định giá trị của mẫu máu ngoại vi trong xét nghiệm sinh thiết lỏng (liquid biopsy), khi chỉ số OR trên mẫu máu đạt mức rất cao (13,549 trên BRCA1 và 45,972 trên p16INK4a).
  • Chuẩn hóa thành công hệ thống mồi và điều kiện phản ứng Nested-MSP (nhiệt độ lai 56°C cho cặp mồi ngoài BRCA1), chứng minh khả năng ứng dụng thực tế trên nguồn mẫu mô paraffin thu thập tại Việt Nam.

Tác giả đề xuất các kiến nghị:

  • Tiếp tục mở rộng khảo sát thực nghiệm trên tập mẫu bệnh phẩm ung thư vú lớn hơn tại các bệnh viện ở Việt Nam để tăng độ tin cậy thống kê cho các nghiên cứu dịch tễ phân tử trong nước.
  • Triển khai phân tích song song giữa mẫu mô và mẫu máu ngoại vi của cùng một bệnh nhân nhằm đánh giá độ nhạy và độ đặc hiệu của phương pháp chẩn đoán không xâm lấn.
  • Gửi giải trình tự Sanger các sản phẩm khuếch đại Nested-MSP để tái kiểm chứng độ chính xác của các phân đoạn DNA đích sau biến đổi bisulfite.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Tính thiên vị công bố (Publication bias): Đồ thị Funnel plot trên cả hai tập dữ liệu BRCA1 và p16INK4a bộc lộ sự bất đối xứng, phản ánh hiện tượng thiên vị xuất bản đối với các nghiên cứu có kết quả dương tính hoặc cỡ mẫu nhỏ.
  • Dữ liệu lâm sàng trong các công bố quốc tế chưa đồng bộ: Việc thiếu hụt dữ liệu chi tiết về tuổi tác, giai đoạn TNM, phân độ mô học và tình trạng di căn hạch trong các bài báo ca - chứng gốc dẫn đến việc phân tích tương quan giữa tính chất methyl hóa và các đặc điểm giải phẫu bệnh học chưa đạt được ý nghĩa thống kê ($P > 0,05$), ngoại trừ mối tương quan giữa methyl hóa BRCA1 và giai đoạn bệnh tiến triển ($P = 0,005$).
  • Giới hạn phần thực nghiệm: Số lượng mẫu mô vú paraffin khảo sát thực nghiệm tại phòng thí nghiệm còn ít; nghiên cứu chưa có điều kiện gửi giải trình tự nucleotide sản phẩm PCR vòng 2 để tái kiểm chứng trực tiếp từng vị trí cytosine bị methyl hóa.

Hướng nghiên cứu tiếp:

  • Mở rộng quy mô thực nghiệm trên các bộ mẫu mô và mẫu huyết tương/huyết thanh thu thập từ các cơ sở ung bướu tại Việt Nam.
  • Nghiên cứu kết hợp đồng thời nhiều gen ức chế khối u (panel đa gen gồm BRCA1, p16INK4a, RASSF1A, GSTP1...) để xây dựng quy trình chẩn đoán đa dấu ấn sinh học có độ nhạy và độ đặc hiệu cao hơn.
  • Ứng dụng các kỹ thuật định lượng hiện đại như Quantitative Methylation-Specific PCR (Q-MSP) hoặc giải trình tự thế hệ mới xử lý bisulfite (Bisulfite Pyrosequencing/NGS) nhằm định lượng chính xác mức độ methyl hóa.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên khảo hữu ích cho:

  • Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Công nghệ Sinh học Y Dược, Sinh học Phân tử, Di truyền Y học và Ung bướu học.
  • Cán bộ nghiên cứu và kỹ thuật viên xét nghiệm tại các phòng thí nghiệm sinh học phân tử lâm sàng, trung tâm xét nghiệm di truyền và chẩn đoán phân tử.

Các phần nội dung có giá trị tham khảo ứng dụng cao:

  • Quy trình chuẩn bị mẫu và thực nghiệm: Các bước xử lý mẫu mô paraffin (tách sáp bằng xylene, hoàn nguyên mô bằng dãy cồn, ly giải proteinase K kéo dài) và quy trình chuyển đổi bisulfite bằng bộ kit EZ DNA Methylation-Gold Kit là chỉ dẫn kỹ thuật có thể ứng dụng trực tiếp cho các nghiên cứu liên quan đến phân tích DNA biến tính.
  • Hệ thống thông số mồi Nested-MSP: Bảng trình tự mồi, nhiệt độ nóng chảy $T_m$, năng lượng tự do cấu trúc thứ cấp và chu trình luân nhiệt của hai gen BRCA1 và p16INK4a là nguồn dữ liệu chuẩn hóa, giúp tiết kiệm thời gian thiết kế và tối ưu hóa phản ứng PCR cho các đề tài tiếp nối.
  • Phương pháp luận phân tích tổng hợp: Khung nghiên cứu meta-analysis hoàn chỉnh từ chiến lược tìm kiếm theo PICO, trích xuất dữ liệu, kiểm định tính bất đồng nhất ($I^2$), lựa chọn mô hình thống kê (Fixed/Random effects) đến phân tích phân nhóm và đánh giá publication bias trên MedCalc là tài liệu tham khảo phương pháp cho các nghiên cứu dịch tễ học phân tử.

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

1. Chỉ số tỷ suất chênh (OR) và nguy cơ tương đối (RR) phản ánh mối liên quan giữa tính chất methyl hóa promoter gen BRCA1, p16INK4a với ung thư vú là bao nhiêu?
Kết quả phân tích tổng hợp theo mô hình ảnh hưởng ngẫu nhiên ghi nhận: Đối với gen BRCA1, $OR = 8,979$ ($95% CI: 4,359 - 18,497; P < 0,0001$) và $RR = 4,973$ ($95% CI: 2,523 - 9,802; P < 0,0001$). Đối với gen p16INK4a, $OR = 12,925$ ($95% CI: 4,585 - 36,429; P < 0,0001$) và $RR = 8,920$ ($95% CI: 2,203 - 36,101; P = 0,002$). Các số liệu này chứng minh sự hiện diện của hiện tượng tăng methyl hóa làm tăng nguy cơ mắc ung thư vú lên từ 4,973 đến 8,92 lần.

2. Phân tích tổng hợp cho thấy sự khác biệt về tính chất methyl hóa gen p16INK4a giữa quần thể bệnh nhân châu Á và châu Âu như thế nào?
Sự khác biệt thể hiện rất rõ nét qua chỉ số tỷ suất chênh: Ở quần thể người châu Á (tổng hợp từ 9 nghiên cứu), tính chất methyl hóa p16INK4a ở nhóm bệnh cao gấp 23,573 lần so với nhóm lành ($OR = 23,573; 95% CI: 13,652 - 40,703; P < 0,001$). Trong khi đó, ở quần thể người châu Âu (2 nghiên cứu), chỉ số này chỉ đạt $OR = 1,254$ ($95% CI: 0,285 - 4,728$) và không có ý nghĩa thống kê ($P = 0,838$).

3. Kích thước các phân đoạn khuếch đại (amplicon) dự kiến của các cặp mồi Nested-MSP cho hai gen mục tiêu là bao nhiêu?
Theo kết quả khảo sát tin sinh học và tài liệu tham khảo:

  • Gen BRCA1: Cặp mồi ngoài vòng 1 (BRCA1_F/R) tạo sản phẩm kích thước 143 bp; cặp mồi allele methyl hóa vòng 2 (BRCA1_MF/MR) tạo sản phẩm 68 bp; cặp mồi allele không methyl hóa vòng 2 (BRCA1_UF/UR) tạo sản phẩm 68 bp.
  • Gen p16INK4a: Cặp mồi ngoài vòng 1 (P16_F/R) tạo sản phẩm kích thước 277 bp; cặp mồi allele methyl hóa vòng 2 (P16_MF/MR) tạo sản phẩm 150 bp; cặp mồi allele không methyl hóa vòng 2 (P16_UF/UR) tạo sản phẩm 151 bp.

4. Vì sao tác giả lựa chọn kỹ thuật Nested-MSP thay vì kỹ thuật MSP thông thường cho nghiên cứu thực nghiệm?
Tác giả lựa chọn Nested-MSP vì kỹ thuật này gồm hai vòng phản ứng nối tiếp: vòng 1 khuếch đại lượng lớn đoạn trình tự đích phổ quát, sau đó dùng sản phẩm này làm khuôn cho vòng 2 với các mồi đặc hiệu trạng thái methyl hóa. Quy trình cải tiến này giúp tăng đáng kể độ nhạy và độ đặc hiệu khi phân tích mẫu mô đúc paraffin - loại mẫu thường bị phân mảnh DNA và chịu ảnh hưởng tiêu cực từ quá trình cố định mô bằng formalin.

5. Kết quả phân tích có tìm thấy mối tương quan có ý nghĩa thống kê giữa tính chất methyl hóa hai gen mục tiêu với các đặc điểm bệnh học ung thư vú hay không?
Hầu hết các đặc điểm bệnh học gồm độ tuổi (cao so với thấp), phân độ mô học (cao so với thấp) và biểu hiện di căn hạch lympho (có so với không) đều cho chỉ số OR dao động từ 1 đến 2 nhưng giá trị $P > 0,05$ (không có ý nghĩa thống kê). Riêng đối với gen BRCA1, bệnh nhân ở giai đoạn muộn (III và IV) có tỷ lệ methyl hóa cao hơn có ý nghĩa thống kê so với bệnh nhân giai đoạn sớm (I và II) với $OR = 1,543$ ($95% CI: 1,142 - 2,085; P = 0,005$).

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Hoàng Trinh đã hoàn thành toàn diện việc khảo sát tính chất methyl hóa bất thường trên vùng promoter hai gen ức chế khối u BRCA1 và p16INK4a thông qua phân tích tổng hợp dịch tễ học phân tử kết hợp khảo sát tin sinh học và thực nghiệm phòng thí nghiệm. Kết quả meta-analysis định lượng rõ rệt mối tương quan thuận giữa hiện tượng tăng methyl hóa hai gen này với nguy cơ ung thư vú trên quy mô toàn cầu, đặc biệt chỉ ra độ chênh methyl hóa cao vượt trội ở quần thể người châu Á và trên loại mẫu máu ngoại vi. Đồng thời, nghiên cứu đã chuẩn hóa thành công các thông số kỹ thuật của bộ mồi và quy trình thực nghiệm Nested-MSP trên mẫu mô paraffin tại Việt Nam, đặt nền tảng phương pháp luận vững chắc cho các nghiên cứu ứng dụng dấu chứng sinh học phân tử vào công tác chẩn đoán sớm và tiên lượng ung thư vú trong nước.