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

Vô sinh là một vấn đề sức khỏe toàn cầu gây ra những ảnh hưởng tâm lý và xã hội sâu sắc. Tại Việt Nam, nghiên cứu toàn quốc của Bệnh viện Phụ sản Trung ương và Đại học Y Hà Nội ghi nhận tỷ lệ vô sinh ở các cặp vợ chồng trong độ tuổi sinh đẻ là 7,7%. Trên phạm vi thế giới, vô sinh ảnh hưởng đến khoảng 7% nam giới, và các nguyên nhân xuất phát từ "yếu tố nam" chiếm khoảng 40% đến 50% tổng số ca bệnh. Trong đó, rối loạn sinh tinh là dạng phổ biến nhất dẫn đến suy giảm chức năng tinh hoàn nguyên phát vô căn.

Sinh tinh là quá trình biệt hóa tế bào phức tạp từ tinh nguyên bào thành tinh trùng trưởng thành, chịu sự điều hòa của khoảng 2.000 gen, với khoảng 20% trong số đó chỉ biểu hiện đặc hiệu ở dòng tế bào mầm nam. Các đột biến hoặc biến đổi di truyền ở những gen này có thể gây gián đoạn quá trình tạo tinh trùng. Trong cơ chế biểu sinh, các RNA tương tác với Piwi (piRNA) đóng vai trò bảo vệ hệ gen tế bào mầm khỏi các yếu tố chuyển vị (transposon elements). Gen TDRD1 (Tudor domain-containing 1), nằm trên cánh dài của nhiễm sắc thể số 10 (vị trí 10q25), mã hóa protein chứa domain Tudor có chức năng phối hợp với protein Piwi nhằm ức chế retrotransposon (như LINE1) trong quá trình giảm phân và sinh tinh. Các mô hình chuột biến đổi gen bất hoạt TDRD1 cho thấy tình trạng vô sinh nam hoàn toàn do khiếm khuyết sinh tinh ở giai đoạn hợp tử - tiếp hợp.

Khóa luận tốt nghiệp của tác giả Nguyễn Minh Nguyệt (ngành Công nghệ Sinh học, Khoa Công nghệ Sinh học, Học viện Nông nghiệp Việt Nam, thực hiện tại Phòng Thí nghiệm Di truyền hệ gen người – Viện Nghiên cứu hệ gen) được triển khai nhằm thực hiện hai mục tiêu chính:

  1. Xác định sự phân bố kiểu gen và tần số alen của đa hình đơn nucleotide (SNP) TDRD1 rs541192490 ở quần thể người Việt Nam.
  2. Đánh giá mối liên quan giữa biến thể TDRD1 rs541192490 với nguy cơ vô sinh nam ở quần thể người Việt Nam.

Đối tượng và phạm vi nghiên cứu bao gồm 310 nam giới Việt Nam:

  • Nhóm bệnh (cases): Gồm 160 bệnh nhân nam được chẩn đoán vô sinh vô căn (thiểu tinh, vô tinh không do tắc nghẽn, tinh trùng bất động hoặc dị dạng), có bộ nhiễm sắc thể (karyotype) và vùng AZF trên nhiễm sắc thể Y bình thường, không có tiền sử viêm nhiễm hay bệnh lý khác ảnh hưởng đến khả năng sinh sản, được thu thập ngẫu nhiên tại Bệnh viện Đại học Y Hà Nội năm 2019.
  • Nhóm đối chứng (controls): Gồm 150 nam giới khỏe mạnh đến từ một số tỉnh thành tại Việt Nam, đã có ít nhất một con bằng phương pháp thụ thai tự nhiên và có các chỉ số tinh dịch đồ bình thường theo tiêu chuẩn của Tổ chức Y tế Thế giới (WHO).

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

Khung lý thuyết và tiêu chuẩn đánh giá

Nghiên cứu dựa trên cơ sở sinh học phân tử của quá trình sinh tinh và cơ chế ức chế yếu tố chuyển vị qua con đường piRNA:

  • Tiêu chuẩn tinh dịch đồ WHO: Đánh giá các thông số cơ bản gồm thể tích tinh dịch, mật độ tinh trùng, tổng số tinh trùng, độ di động tiến tới, độ sống và hình thái bình thường qua các ấn bản hướng dẫn của WHO (từ ấn bản 1 năm 1980 đến ấn bản 5 năm 2010).
  • Cơ chế hoạt động của TDRD1: Protein TDRD1 chứa 4 domain Tudor kinh điển và 1 domain ngón tay kẽm MYND, định vị tại thể nuage (xi măng gian ty thể và thể nhiễm sắc) trong tế bào sinh tinh, liên kết đặc hiệu với protein MILI (PIWIL2) để methyl hóa DNA de novo tại promoter của LINE1, ngăn ngừa tái hoạt hóa transposon trong giảm phân.
Tiêu chí tinh dịch đồ theo WHO Ấn bản 1 (1980) Ấn bản 2 (1987) Ấn bản 3 (1992) Ấn bản 4 (1999) Ấn bản 5 (2010)
Thể tích tinh dịch Không xác định $\ge 2{,}0\text{ mL}$ $\ge 2{,}0\text{ mL}$ $\ge 2{,}0\text{ mL}$ $\ge 1{,}5\text{ mL}$
Mật độ tinh trùng $20 - 200 \times 10^6/\text{mL}$ $\ge 20 \times 10^6/\text{mL}$ $\ge 20 \times 10^6/\text{mL}$ $\ge 20 \times 10^6/\text{mL}$ $\ge 15 \times 10^6/\text{mL}$
Tổng số tinh trùng $\ge 40 \times 10^6/\text{mL}$ $\ge 40 \times 10^6/\text{mL}$ $\ge 40 \times 10^6/\text{mL}$ $\ge 40 \times 10^6/\text{mL}$ $\ge 39 \times 10^6/\text{lần xuất tinh}$
Di động tiến tới $\ge 60%$ $\ge 50%$ $\ge 50%$ $\ge 50%$ $\ge 32%$
Tỷ lệ sống Không xác định $\ge 50%$ $\ge 75%$ $\ge 75%$ $\ge 58%$
Hình thái bình thường $\ge 80{,}5%$ $\ge 50%$ $\ge 30%$ $\ge 15%$ $\ge 4%$

Phương pháp nghiên cứu thực nghiệm

Quy trình thực nghiệm được tiến hành theo các bước chuẩn mực trong sinh học phân tử:

  • Tách chiết DNA tổng số: Tách chiết từ 200 µL máu toàn phần chống đông bằng EDTA của 310 đối tượng nghiên cứu bằng kit thương mại GeneJET Whole Blood Genomic DNA Purification Kit (Thermo Scientific). Kiểm tra nồng độ và độ tinh sạch của DNA bằng máy quang phổ NanoDrop Lite và điện di trên gel agarose 1%.
  • Thiết kế mồi PCR: Sử dụng công cụ trực tuyến Primer-BLAST và OligoAnalyzer để thiết kế cặp mồi đặc hiệu khuếch đại đoạn gen chứa SNP TDRD1 rs541192490 (vị trí biến thể NM_198795.1009G>A):
    • Mồi xuôi (F): 5'-TTGAGTGCCTCAGAACAGGAG-3'
    • Mồi ngược (R): 5'-GTTTCTCACCCTGGGCAGTAG-3'
  • Kỹ thuật PCR: Phản ứng khuếch đại thể tích 5 µL sử dụng enzyme DreamTaq DNA Polymerase (Thermo Fisher) và 10 ng DNA khuôn, tạo ra sản phẩm PCR có kích thước 442 bp.
  • Kỹ thuật PCR-RFLP (Restriction Fragment Length Polymorphism): Cắt sản phẩm PCR bằng enzyme giới hạn Psp1406I (AclI, Thermo Scientific) nhận biết trình tự cắt đặc hiệu. Sản phẩm cắt ủ ở 37°C trong 5 giờ và phân tách trên gel agarose 1,5%:
    • Kiểu gen đồng hợp hoang dại (GG): Bị cắt thành 2 đoạn (145 bp và 297 bp).
    • Kiểu gen dị hợp tử (GA): Cắt không hoàn toàn tạo 3 đoạn (145 bp, 297 bp và 435 bp).
    • Kiểu gen đồng hợp đột biến (AA): Mất vị trí nhận biết của enzyme, giữ nguyên 1 đoạn 435 bp.
  • Giải trình tự Sanger kiểm chứng: Lựa chọn ngẫu nhiên khoảng 5% tổng số mẫu (15 mẫu gồm 8 bệnh nhân và 7 đối chứng), tinh sạch sản phẩm PCR bằng GeneJET PCR Purification Kit, phản ứng giải trình tự chu kỳ với BigDye Terminator v3.1 trên hệ thống ABI PRISM 3500-Avant Genetic Analyzer và phân tích kết quả bằng phần mềm SnapGene Viewer.
  • Xử lý số liệu thống kê: Sử dụng phần mềm thống kê R (phiên bản 4.2) và Microsoft Excel. Kiểm định tính cân bằng Hardy-Weinberg (HWE) bằng phép thử Chi bình phương ($\chi^2$). Đánh giá mối liên quan giữa biến thể di truyền và nguy cơ vô sinh nam qua tỷ số chênh (Odds Ratio - OR), khoảng tin cậy 95% (95% CI) theo 3 mô hình di truyền: cộng gộp (additive), trội (dominant), và lặn (recessive). Ngưỡng ý nghĩa thống kê xác định khi $p < 0{,}05$.

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

Part 1: Introduction (Mở đầu)

Chương mở đầu trình bày tổng quan về thực trạng vô sinh nam trên thế giới và tại Việt Nam, nêu bật tầm quan trọng của các yếu tố di truyền trong rối loạn sinh tinh. Tác giả nhấn mạnh vai trò của hệ thống piRNA và giới thiệu vị trí, chức năng của gen ứng viên TDRD1. Từ đó, tác giả đưa ra lý do lựa chọn nghiên cứu biến thể TDRD1 rs541192490 và xác định 2 mục tiêu nghiên cứu cụ thể.

Part 2: Literature Review (Tổng quan tài liệu)

Chương này hệ thống hóa các kiến thức chuyên sâu về:

  • Định nghĩa, phân loại vô sinh (vô sinh nguyên phát và thứ phát), các dạng tổn thương sinh tinh (oligozoospermia, asthenozoospermia, teratozoospermia, non-obstructive azoospermia - NOA).
  • Quy trình thăm khám lâm sàng và các tiêu chuẩn phân tích tinh dịch đồ theo WHO qua các giai đoạn 1980, 1987, 1992, 1999 và 2010.
  • Các nhóm nguyên nhân gây vô sinh nam: trục nội tiết hạ đồi - tuyến yên - tinh hoàn, bệnh lý thực thể (giãn tĩnh mạch thừng tinh, tinh hoàn ẩn, nhiễm trùng đường sinh dục), lối sống và môi trường (hút thuốc, rượu bia, hóa chất, phóng xạ), và các yếu tố di truyền (bất thường nhiễm sắc thể như hội chứng Klinefelter 47,XXY; vi mất đoạn vùng AZF trên NST Y; đột biến gen CFTR, SHBG, ESR1, ESR2, FSHR).
  • Cấu trúc, sinh tổng hợp piRNA và chức năng sinh học của họ protein Tudor (TDRDs), đặc biệt là TDRD1 trong kiểm soát transposon bảo vệ tế bào mầm.
  • Nguyên lý kỹ thuật PCR-RFLP và các công nghệ định kiểu gen SNP khác.
  • Lược sử các nghiên cứu liên quan đến đa hình gen sinh tinh trên thế giới (các nghiên cứu GWAS, nghiên cứu gen thuộc con đường piRNA như HIWI, TDRD1 rs77559927) và các nghiên cứu di truyền vô sinh nam đã thực hiện tại Việt Nam (gen NAT2, GSTs, AHRR).

Part 3: Materials and Methods (Đối tượng và phương pháp nghiên cứu)

Trình bày chi tiết tiêu chuẩn chọn mẫu bệnh và mẫu chứng (310 đối tượng); danh mục hóa chất, sinh phẩm, cặp mồi PCR, enzyme cắt giới hạn Psp1406I; thông số thiết bị sử dụng; các bước chi tiết trong quy trình tách chiết DNA từ máu toàn phần, kiểm tra chất lượng DNA, thiết lập chu nhiệt phản ứng PCR, phản ứng cắt enzyme RFLP, tinh sạch sản phẩm PCR, giải trình tự Sanger và phương pháp phân tích thống kê y sinh học.

Part 4: Results (Kết quả nghiên cứu)

Chương kết quả công bố các dữ liệu thực nghiệm thu được:

  • Kết quả kiểm tra DNA tổng số: Băng điện di DNA sắc nét, không bị thoái hóa; chỉ số đo quang học nồng độ đạt chuẩn.
  • Kết quả khuếch đại PCR: Sản phẩm PCR duy nhất có kích thước chính xác 442 bp.
  • Kết quả phân tích kiểu gen bằng PCR-RFLP: Xác định rõ ràng 3 kiểu gen GG, GA, AA trên gel agarose 1,5%.
  • Tần số kiểu gen và alen: Quần thể nghiên cứu tuân theo cân bằng Hardy-Weinberg ($p > 0{,}05$). Tần số alen đột biến A (minor allele frequency) ở nhóm bệnh là 0,041 và nhóm đối chứng là 0,037.
  • Kết quả giải trình tự Sanger: 15 mẫu giải trình tự cho kết quả hoàn toàn trùng khớp với kết quả định kiểu gen bằng PCR-RFLP.
  • Phân tích mối liên quan: Các mô hình di truyền cộng gộp, trội, lặn đều cho giá trị $p > 0{,}05$ ($p = 0{,}952$ ở mô hình trội; $p = 0{,}561$ ở mô hình lặn).

Part 5: Discussion (Bàn luận)

Tác giả thảo luận về chức năng của protein TDRD1 gồm 4 domain Tudor và 1 domain MYND định vị tại thể nuage tế bào mầm, phân tích các bằng chứng thực nghiệm trên chuột đột biến TDRD1 dẫn đến bất hoạt sinh tinh hoàn toàn. So sánh kết quả thu được với các công bố quốc tế (như nghiên cứu của Zhu và cộng sự tại Trung Quốc về biến thể rs77559927 trên vùng 5'UTR của TDRD1 có tác dụng bảo vệ giảm nguy cơ tổn thương sinh tinh). Tác giả lý giải nguyên nhân không tìm thấy mối liên quan có ý nghĩa thống kê của SNP rs541192490 ở người Việt Nam có thể do tính đặc thù quần thể/chủng tộc hoặc biến thể này có thể liên kết với các đa hình khác chưa được khảo sát.

Part 6: Conclusion (Kết luận)

Tóm tắt ngắn gọn các phát hiện chính của khóa luận về sự phân bố kiểu gen, tần số alen của SNP rs541192490 và kết luận về mối liên quan với vô sinh nam ở quần thể người Việt Nam; đồng thời khuyến nghị hướng nghiên cứu tiếp theo đối với các biến thể khác trên gen TDRD1.

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

Số liệu thực nghiệm và thống kê di truyền

Kết quả phân bố kiểu gen và tần số alen của SNP TDRD1 rs541192490 trên 310 mẫu nghiên cứu được tổng hợp trong bảng sau:

Nhóm nghiên cứu Số lượng ($n$) Kiểu gen GG Kiểu gen GA Kiểu gen AA Tần số alen G Tần số alen A Giá trị $p$ (HWE)
Nhóm bệnh (Cases) 160 148 11 1 0,959 0,041 $> 0{,}05$
Nhóm chứng (Controls) 150 139 11 0 0,963 0,037 $> 0{,}05$

Kết quả phân tích mối liên quan giữa đa hình TDRD1 rs541192490 và nguy cơ vô sinh nam theo các mô hình di truyền:

Mô hình di truyền Kiểu gen Nhóm bệnh ($n=160$) Nhóm chứng ($n=150$) Tỷ số chênh (OR) 95% CI Giá trị $p$ ($\chi^2$)
Cộng gộp (Additive) GG / GA / AA 148 / 11 / 1 139 / 11 / 0 - - $> 0{,}05$
Trội (Dominant) GG vs. GA+AA 148 vs. 12 139 vs. 11 1,02 0,44 – 2,39 0,952
Lặn (Recessive) GG+GA vs. AA 159 vs. 1 150 vs. 0 - - 0,561

Ghi chú: $n$ là số lượng người tham gia; OR: Odds Ratio; 95% CI: Khoảng tin cậy 95%; giá trị $p$ tính bằng kiểm định Chi bình phương ($\chi^2$).

Đóng góp khoa học của đề tài

  • Tính tiên phong: Đây là công trình đầu tiên tại Việt Nam khảo sát biến thể sai nghĩa (missense variant) rs541192490 (NM_198795.1009G>A) trên gen TDRD1 ở bệnh nhân vô sinh nam.
  • Đóng góp dữ liệu quần thể: Cung cấp dữ liệu gốc về tần số phân bố kiểu gen và alen của biến thể TDRD1 rs541192490 ở người Việt Nam, khẳng định quần thể tuân theo định luật cân bằng Hardy-Weinberg.
  • Xây dựng quy trình kỹ thuật: Hoàn thiện và chuẩn hóa quy trình phân tích định kiểu gen SNP rs541192490 bằng phương pháp PCR-RFLP sử dụng enzyme cắt Psp1406I, có độ chính xác tương đương 100% với phương pháp giải trình tự gen Sanger chuẩn thức, giúp giảm chi phí và thời gian xét nghiệm.

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

  • Hạn chế:
    • Quy mô mẫu nghiên cứu dừng lại ở 310 người (160 bệnh nhân và 150 đối chứng), đây là cỡ mẫu ở mức trung bình đối với các nghiên cứu di truyền quần thể bệnh lý phức tạp.
    • Nghiên cứu chỉ tập trung phân tích đơn lẻ một biến thể rs541192490 trên gen TDRD1, chưa đánh giá toàn diện các biến thể khác hoặc phân tích haplotype liên kết trong vùng gen.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng quy mô mẫu nghiên cứu trên các nhóm bệnh nhân phân loại chi tiết theo từng thể tổn thương tinh dịch đồ (như riêng nhóm vô tinh không do tắc nghẽn NOA hoặc thiểu tinh nặng).
    • Khảo sát thêm các đa hình đơn nucleotide khác trên gen TDRD1 (ví dụ vùng promoter, 5'UTR, 3'UTR) và các gen khác thuộc con đường piRNA sinh tinh (PIWIL1, PIWIL2, PIWIL4, TDRD9) để làm rõ cơ chế tương tác đa gen trong vô sinh nam tại Việt Nam.

Giá trị tham khảo

Khóa luận có giá trị tham khảo hữu ích cho:

  • Sinh viên và học viên ngành Công nghệ Sinh học, Sinh học Y học, Di truyền học: Nguồn tài liệu thực nghiệm rõ ràng về cách thiết kế mồi chuyên biệt, tối ưu hóa điều kiện PCR, ứng dụng enzyme giới hạn trong kỹ thuật PCR-RFLP và phương pháp giải trình tự Sanger.
  • Các nghiên cứu viên y sinh học: Cung cấp mô hình phương pháp luận nghiên cứu bệnh - chứng (case-control study) trong dịch tễ học di truyền, bao gồm cách thu thập mẫu, sàng lọc tiêu chuẩn lâm sàng, kiểm tra cân bằng Hardy-Weinberg và tính toán các chỉ số OR, 95% CI trên phần mềm thống kê R.
  • Các bác sĩ và chuyên gia hỗ trợ sinh sản: Cung cấp thêm dẫn liệu khoa học về đặc điểm di truyền phân tử của quần thể nam giới vô sinh vô căn tại Việt Nam.

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

1. Khóa luận khảo sát đa hình trên gen nào và sử dụng phương pháp định kiểu gen chính nào?
Khóa luận nghiên cứu biến thể đơn nucleotide rs541192490 (NM_198795.1009G>A) trên gen TDRD1 (Tudor domain-containing 1) định vị tại nhiễm sắc thể 10q25. Phương pháp định kiểu gen chính được sử dụng là PCR-RFLP (khuếch đại PCR đoạn gen 442 bp và cắt bằng enzyme giới hạn Psp1406I), kết hợp kiểm chứng độ chính xác bằng giải trình tự Sanger.

2. Quy mô và tiêu chuẩn lựa chọn mẫu trong nghiên cứu này như thế nào?
Nghiên cứu thực hiện trên 310 nam giới Việt Nam, gồm 160 bệnh nhân vô sinh vô căn (thiểu tinh, vô tinh không do tắc nghẽn, tinh trùng dị dạng/bất động, karyotype và vùng AZF bình thường) thu thập tại Bệnh viện Đại học Y Hà Nội năm 2019 và 150 đối chứng khỏe mạnh đã có con tự nhiên, có tinh dịch đồ đạt chuẩn theo hướng dẫn của WHO.

3. Enzyme cắt giới hạn nào được dùng trong kỹ thuật PCR-RFLP và kích thước các phân đoạn cắt tương ứng với từng kiểu gen là bao nhiêu?
Nghiên cứu sử dụng enzyme cắt giới hạn Psp1406I (AclI). Kích thước các phân đoạn cắt tương ứng:

  • Kiểu gen GG (đồng hợp hoang dại): 2 đoạn (145 bp và 297 bp).
  • Kiểu gen GA (dị hợp tử): 3 đoạn (145 bp, 297 bp và 435 bp).
  • Kiểu gen AA (đồng hợp đột biến): 1 đoạn (435 bp, không bị cắt).

4. Kết quả kiểm định giải trình tự Sanger trên các mẫu ngẫu nhiên cho thấy điều gì?
Tổng cộng 15 mẫu (8 bệnh nhân và 7 đối chứng, chiếm khoảng 5% tổng số mẫu) được giải trình tự Sanger trên máy ABI PRISM 3500-Avant. Kết quả giải trình tự hoàn toàn trùng khớp 100% với kết quả đọc kiểu gen qua kỹ thuật PCR-RFLP, xác nhận độ tin cậy của quy trình thực nghiệm.

5. Kết luận chính về mối liên quan giữa biến thể TDRD1 rs541192490 và vô sinh nam tại Việt Nam là gì?
Phân tích thống kê theo cả 3 mô hình di truyền (cộng gộp, trội, lặn) đều cho giá trị $p > 0{,}05$ (mô hình trội có $\text{OR} = 1{,}02$, $95%\text{ CI} = 0{,}44 - 2{,}39$, $p = 0{,}952$; mô hình lặn có $p = 0{,}561$). Do đó, nghiên cứu kết luận chưa phát hiện mối liên quan có ý nghĩa thống kê giữa biến thể TDRD1 rs541192490 với nguy cơ vô sinh nam ở quần thể người Việt Nam.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Minh Nguyet đã thiết lập thành công quy trình kỹ thuật PCR-RFLP kết hợp giải trình tự Sanger để định kiểu gen biến thể TDRD1 rs541192490 trên 310 mẫu nam giới Việt Nam. Tần số kiểu gen và alen của biến thể này trong quần thể nghiên cứu tuân theo định luật cân bằng Hardy-Weinberg, với tần số alen đột biến A ở nhóm bệnh là 0,041 và nhóm chứng là 0,037. Kết quả phân tích thống kê cho thấy biến thể TDRD1 rs541192490 không có mối liên quan có ý nghĩa thống kê với nguy cơ vô sinh nam ở quần thể nghiên cứu, đóng góp thêm dữ liệu cơ sở cho các nghiên cứu di truyền học sinh sản tiếp theo tại Việt Nam.