Tổng quan nghiên cứu

Khiếm thính là khuyết tật giác quan phổ biến nhất ở người, ảnh hưởng sâu sắc đến sự phát triển nhận thức, ngôn ngữ và khả năng hòa nhập xã hội. Theo thống kê của Tổ chức Y tế Thế giới công bố năm 2018, toàn cầu có khoảng 466 triệu người bị suy giảm thính lực, chiếm hơn 5% dân số thế giới. Tỷ lệ mắc bệnh ở trẻ sơ sinh dao động từ 1 đến 4 trên 1.000 ca sinh, trong đó có khoảng 50% đến 70% trường hợp khiếm thính bẩm sinh bắt nguồn từ các đột biến gen di truyền. Trong nhóm khiếm thính di truyền, khoảng 70% thuộc dạng không hội chứng và 30% thuộc dạng có hội chứng đi kèm nhiều tổn thương phức tạp.

Phương pháp chẩn đoán truyền thống chủ yếu dựa trên đo âm ốc tai hoặc giải trình tự Sanger từng đoạn gen đơn lẻ, bộc lộ nhiều hạn chế về chi phí và thời gian đối với căn bệnh có tính đa hình cao với hơn 150 gen liên quan. Xuất phát từ nhu cầu thực tiễn đó, nghiên cứu được triển khai nhằm mục tiêu xây dựng quy trình kỹ thuật giải trình tự DNA thế hệ mới phát hiện chính xác các đột biến trên 7 gen đích phổ biến gồm GJB2, GJB3, GJB6, KCNQ4, MYO7A, MYO15A và SLC26A4. Đề tài được thực hiện tại Công ty TNHH Công nghệ Sinh học Khoa Thương từ ngày 13/08/2018 đến ngày 02/12/2018.

Công trình mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp nền tảng xét nghiệm di truyền có độ chính xác trên 99,9% và độ phủ phân tích đạt trên 100X. Dữ liệu từ quy trình này đóng vai trò then chốt trong việc xác định căn nguyên bệnh lý, hỗ trợ chỉ định cấy ốc tai điện tử sớm cho trẻ nhỏ và phục vụ hiệu quả cho công tác tư vấn di truyền tiền hôn nhân tại Việt Nam.

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 lý thuyết di truyền phân tử về bệnh khiếm thính và các hướng dẫn chẩn đoán của Hiệp hội Di truyền Y khoa Hoa Kỳ. Về mặt phân tử, hệ thống thính giác của người hoạt động dựa trên mạng lưới tương tác phức tạp của nhiều nhóm protein cấu trúc, kênh ion và protein liên kết khe:

  • Họ protein Connexin: Gồm các gen GJB2, GJB3 và GJB6 mã hóa cho các protein connexin 26, connexin 31 và connexin 30. Chúng tham gia cấu tạo các kênh liên lạc giữa các tế bào ở ốc tai, trong đó các sai hỏng trên gen GJB2 chiếm từ 50% đến 80% tổng số ca khiếm thính lặn trên nhiễm sắc thể thường.
  • Kênh ion màng tế bào: Gen KCNQ4 mã hóa cho kênh kali điện thế màng tế bào, giữ vai trò quan trọng trong việc tái phục hồi ion kali ở tai trong và là nguyên nhân phổ biến của dạng khiếm thính trội.
  • Họ protein Motor Myosin: Gồm gen MYO7A và gen MYO15A mã hóa cho các phân tử myosin vận chuyển ở lông ốc tai, liên quan trực tiếp đến hội chứng Usher và khiếm thính thần kinh giác quan sâu.
  • Protein vận chuyển Anion: Gen SLC26A4 mã hóa cho protein pendrin, liên quan đến hội chứng Pendred và dị tật mở rộng cống tiền đình.

Song song đó, nghiên cứu áp dụng nguyên lý giải trình tự DNA thế hệ mới trên hệ thống Illumina dựa trên cơ chế tổng hợp huỳnh quang kết hợp bắt cặp hai đầu, giúp giải mã hàng triệu đoạn trình tự đồng thời với độ chính xác và thông lượng vượt trội so với kỹ thuật Sanger cổ điển.

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

  • Nguồn dữ liệu và mẫu nghiên cứu: Nghiên cứu sử dụng mẫu chuẩn quốc tế từ Viện Coriell bao gồm mẫu chuẩn dương NA23835 mang đột biến gây bệnh c.101T>C trên gen GJB2, mẫu chuẩn âm NA12878 và 4 mẫu máu toàn phần từ các cá thể người Việt Nam khỏe mạnh không mắc bệnh khiếm thính (mã hiệu HL001, HL002, HL003, HL004). Phương pháp chọn mẫu kiểm chuẩn có chủ đích nhằm thiết lập và thẩm định độ nhạy, độ đặc hiệu của quy trình phân tích.
  • Kỹ thuật Long-range PCR: Thiết kế 28 cặp mồi chuyên biệt bao phủ toàn bộ các exon mã hóa và vùng intron kế cận tối thiểu 15 cặp base của 7 gen đích. Phản ứng khuếch đại phân đoạn lớn được thực hiện bằng enzyme PrimeSTAR GXL DNA Polymerase trên máy luân nhiệt với nhiệt độ lai được tối ưu hóa từ 58°C đến 68°C.
  • Chuẩn bị thư viện và giải trình tự: Sản phẩm PCR được phân mảnh bằng enzyme dsDNA Fragmentase, gắn adapter, làm giàu bằng phản ứng PCR kèm chỉ thị barcode và tinh sạch bằng hạt từ AMPure XP để chọn kích thước thư viện từ 300 đến 600 bp. Quá trình giải trình tự vận hành trên hệ máy Illumina MiniSeq với bộ kit MiniSeq Mid-output 300 chu kỳ.
  • Xử lý và phân tích dữ liệu: Dữ liệu FASTQ thô được kiểm tra chất lượng qua FastQC và MultiQC, gióng hàng lên hệ gen tham chiếu người GRCh38 bằng phần mềm Geneious. Các biến thể được xác định và đối chiếu trên cơ sở dữ liệu quốc tế như ClinVar, 1000 Genomes, dbSNP, phần mềm Alamut Visual và Varsome. Các biến thể gây bệnh hoặc vùng có độ phủ dưới 100X được kiểm chứng lại bằng phương pháp giải trình tự Sanger. Timeline nghiên cứu được hoàn thành tập trung trong 4 tháng cuối năm 2018.

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

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

Nghiên cứu đã xây dựng thành công quy trình hoàn chỉnh từ khâu khuếch đại gen mục tiêu đến giải trình tự và phát hiện biến thể với các kết quả cụ thể:

  • Hiệu quả thiết kế mồi và khuếch đại Long-range PCR: Thiết lập thành công 28 cặp mồi đặc hiệu khuếch đại toàn bộ các vùng exon của 7 gen mục tiêu với kích thước sản phẩm dao động từ 1.117 bp đến 11.031 bp. Các phản ứng PCR cho băng điện di sáng rõ, đặc hiệu, không xuất hiện hiện tượng dimer hay khuếch đại không đặc hiệu.
  • Chất lượng thư viện và thông số chạy máy MiniSeq: Thư viện giải trình tự đạt kích thước chuẩn từ 300 đến 600 bp với nồng độ nạp máy tối ưu 1,4 pM. Quá trình giải trình tự ghi nhận mật độ cluster đạt 219 k/mm², tỷ lệ cụm vượt qua bộ lọc đạt 89,0%, và tỷ lệ nucleotide đạt điểm chất lượng Q30 trở lên đạt 89,9%, vượt xa ngưỡng tiêu chuẩn tối thiểu là 80%.
  • Độ phủ và tính toàn vẹn của dữ liệu tin sinh học: Tỷ lệ sắp xếp gióng cột của các đoạn đọc lên hệ gen tham chiếu GRCh38 đạt trên 85%. Độ sâu giải trình tự trung bình tại các vùng gen mục tiêu đạt trên 100X, tỷ lệ độ phủ ngang lớn hơn 30X đạt trên 80%, đảm bảo độ tin cậy để gọi biến thể di truyền với độ chính xác đạt 99,9%.
  • Độ nhạy và độ đặc hiệu phát hiện biến thể: Quy trình đã phát hiện chính xác đột biến đồng hợp tử gây bệnh c.101T>C trên gen GJB2 ở mẫu chuẩn dương NA23835. Đồng thời, không phát hiện thấy bất kỳ đột biến gây bệnh nào ở mẫu chuẩn âm NA12878 cũng như 4 mẫu máu người bình thường, chứng minh tỷ lệ dương tính giả đạt mức 0%.

Thảo luận kết quả

Các kết quả thực nghiệm khẳng định ưu thế vượt trội của chiến lược giải trình tự nhóm gen mục tiêu kết hợp Long-range PCR so với giải trình tự toàn bộ hệ gen hoặc giải trình tự toàn bộ vùng mã hóa. Bằng cách tập trung vào 7 gen then chốt, chi phí xét nghiệm giảm hơn 100 lần so với phương pháp Sanger truyền thống khi xét trên cùng khối lượng dữ liệu, đồng thời rút ngắn thời gian phân tích từ vài tháng xuống chỉ còn khoảng 1 đến 2 ngày làm việc.

Dữ liệu kích thước thư viện được kiểm soát nghiêm ngặt và thể hiện rõ nét qua biểu đồ điện di mao quản BioAnalyzer 2100 với đỉnh phân bố tập trung chuẩn xác tại khoảng 350 đến 450 bp. Sự phân bố đều của 4 loại nucleotide quanh mức 25% cùng các bảng thống kê độ phủ dọc trên 100X cho thấy quá trình khuếch đại không bị thiên lệch bởi hàm lượng GC cao tại các vùng exon phức tạp.

Tại Việt Nam, các nghiên cứu trước đây chỉ ra rằng đột biến gen GJB2 và SLC26A4 chiếm tỷ lệ phổ biến nhất ở trẻ khiếm thính bẩm sinh. Việc phát hiện chính xác biến thể c.101T>C trên gen GJB2 trong nghiên cứu này mang giá trị lâm sàng đặc biệt quan trọng. Bệnh nhân mang đột biến trên gen GJB2 thường có tổn thương khu trú tại ốc tai trong khi dây thần kinh thính giác vẫn bảo tồn nguyên vẹn. Do đó, phát hiện sớm đột biến này sẽ cung cấp chỉ định vàng cho can thiệp cấy ốc tai điện tử ở giai đoạn vàng từ 1 đến 3 tuổi, mang lại cơ hội phục hồi khả năng nghe và phát triển ngôn ngữ tối ưu cho trẻ.

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

Nhằm hoàn thiện và đưa quy trình xét nghiệm di truyền khiếm thính vào ứng dụng thực tiễn rộng rãi, nghiên cứu đưa ra 4 khuyến nghị hành động trọng tâm:

  1. Mở rộng bảng gen mục tiêu: Bổ sung thêm từ 10 đến 15 gen gây khiếm thính phổ biến khác như STRC, CDH23, TMC1 và OTOF vào bảng xét nghiệm trong vòng 12 tháng tới. Mục tiêu nâng tỷ lệ phát hiện căn nguyên di truyền ở các bệnh nhân khiếm thính không hội chứng lên trên mức 75%. Đơn vị thực hiện: Nhóm nghiên cứu di truyền phân tử phối hợp với các phòng thí nghiệm y sinh.
  2. Chuẩn hóa và thẩm định lâm sàng mở rộng: Triển khai thử nghiệm quy trình trên cỡ mẫu tối thiểu 200 bệnh nhân khiếm thính tại các bệnh viện chuyên khoa Tai Mũi Họng và Nhi đồng trong giai đoạn 2019-2021. Mục tiêu hoàn thiện bộ hồ sơ kỹ thuật đạt tiêu chuẩn xét nghiệm chẩn đoán y khoa ISO 15189. Đơn vị thực hiện: Các bệnh viện chuyên khoa và viện nghiên cứu y sinh học.
  3. Thiết lập ngân hàng dữ liệu biến thể di truyền người Việt: Xây dựng cơ sở dữ liệu hệ gen về các biến thể thính lực bản địa với quy mô thu thập ít nhất 1.000 hồ sơ trong thời gian 2 năm. Cơ sở dữ liệu này giúp các nhà tin sinh học phân loại chính xác các biến thể chưa rõ ý nghĩa chức năng (VUS), giảm thiểu nguy cơ chẩn đoán nhầm lẫn. Đơn vị thực hiện: Các trung tâm tin sinh học và trường đại học.
  4. Tích hợp xét nghiệm gen vào chương trình tầm soát sơ sinh: Đề xuất cơ quan quản lý y tế kết hợp xét nghiệm bảng gen NGS vào quy trình sàng lọc thính lực sau khi kiểm tra âm ốc tai (OAE/ABR) thất bại trong vòng 3 tháng đầu đời. Mục tiêu giúp các bác sĩ đưa ra phác đồ can thiệp cấy ốc tai điện tử sớm cho trẻ từ 1 đến 3 tuổi, đạt tỷ lệ phục hồi giao tiếp trên 90%. Đơn vị thực hiện: Bộ Y tế và các cơ sở y tế dự phòng.

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

Công trình nghiên cứu cung cấp khối lượng kiến thức thực nghiệm và dữ liệu phân tích giá trị, đặc biệt hữu ích cho 4 nhóm đối tượng:

  • Bác sĩ chuyên khoa Tai Mũi Họng và Chuyên gia Di truyền Y học: Sử dụng tài liệu làm cơ sở tham chiếu để chỉ định xét nghiệm phân tử, đánh giá căn nguyên gây điếc bẩm sinh và xây dựng phác đồ tư vấn di truyền tiền hôn nhân, tiền làm tổ cho các gia đình có tiền sử mang gen bệnh.
  • Nhà nghiên cứu và Chuyên viên Công nghệ Sinh học: Tham khảo chi tiết thiết kế 28 cặp mồi Long-range PCR, quy trình tối ưu hóa phản ứng nhiệt và phương pháp chuẩn bị thư viện NGS trên hệ thống Illumina phục vụ phát triển các bộ kit chẩn đoán phân tử.
  • Kỹ sư Tin sinh học (Bioinformatics): Khai thác quy trình xử lý dữ liệu giải trình tự thô dạng FASTQ, các tiêu chí đánh giá kiểm soát chất lượng dữ liệu bằng Picard, QualiMap, và quy trình phân loại biến thể theo chuẩn dữ liệu quốc tế ClinVar, Varsome.
  • Giảng viên, Học viên Cao học và Sinh viên ngành Y Sinh: Sử dụng làm tài liệu học tập mẫu mực về phương pháp nghiên cứu ứng dụng công nghệ giải trình tự thế hệ mới trong y học chính xác và dịch tễ học phân tử.

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

1. Tại sao nghiên cứu lại kết hợp kỹ thuật Long-range PCR với giải trình tự NGS thay vì giải trình tự toàn bộ hệ gen (WGS)?
Giải trình tự toàn bộ hệ gen (WGS) đòi hỏi phân tích khoảng 3,2 tỷ cặp base, làm phát sinh chi phí rất cao và khối lượng dữ liệu khổng lồ. Việc kết hợp Long-range PCR cho 7 gen mục tiêu giúp tập trung làm giàu đúng các vùng exon gây bệnh chính, giảm hơn 100 lần chi phí trên mỗi base, rút ngắn thời gian phân tích xuống dưới 2 ngày và duy trì độ sâu giải trình tự cao trên 100X.

2. Đột biến trên gen GJB2 có ý nghĩa lâm sàng then chốt như thế nào đối với người bệnh?
Gen GJB2 mã hóa cho protein connexin 26 chiếm từ 50% đến 80% trường hợp khiếm thính lặn trên nhiễm sắc thể thường. Các đột biến trên gen này gây tổn thương tế bào ốc tai nhưng vẫn giữ nguyên vẹn chức năng dây thần kinh thính giác. Bệnh nhân mang đột biến gen GJB2 có tiên lượng hồi phục thính lực và phát triển ngôn ngữ xuất sắc khi được cấy ốc tai điện tử sớm trong giai đoạn từ 1 đến 3 tuổi.

3. Tiêu chuẩn chất lượng dữ liệu NGS trong luận văn được đánh giá qua những chỉ số nào?
Dữ liệu giải trình tự phải đáp ứng các tiêu chuẩn kỹ thuật nghiêm ngặt gồm: điểm chất lượng Q30 đạt trên 80% (kết quả thực nghiệm đạt 89,9%), mật độ cluster nằm trong khoảng 160 đến 220 k/mm² (đạt 219 k/mm²), tỷ lệ cụm vượt qua bộ lọc đạt trên 80% (đạt 89,0%), tỷ lệ gióng cột đạt trên 85% và độ sâu phân tích vùng mục tiêu đạt trên 100X.

4. Quy trình này có thể thay thế hoàn toàn kỹ thuật giải trình tự Sanger truyền thống hay không?
Quy trình NGS giúp sàng lọc đồng thời nhiều gen với độ nhạy cao nhưng không thay thế hoàn toàn Sanger. Phương pháp Sanger vẫn được sử dụng như tiêu chuẩn vàng để thẩm định và tái khẳng định lại các biến thể có khả năng gây bệnh hoặc các vùng trình tự có độ phủ đọc dưới 100X nhằm đảm bảo độ chính xác tuyệt đối trước khi trả kết quả lâm sàng.

5. Mẫu bệnh phẩm đầu vào cần đáp ứng những yêu cầu kỹ thuật gì để thực hiện quy trình?
Mẫu xét nghiệm là DNA hệ gen tách chiết từ máu toàn phần bằng cột silica, yêu cầu nồng độ đạt trên 16 ng/µL, tỷ số đo độ sạch OD 260/280 nm nằm trong khoảng 1,6 đến 2,2. Sau khi khuếch đại Long-range PCR, thư viện DNA được chuẩn hóa về nồng độ 10 nM và biến tính pha loãng về nồng độ nạp máy 1,4 pM để bảo đảm phản ứng tạo cụm tối ưu.

Kết luận

  • Xây dựng và tối ưu hóa thành công 28 cặp mồi Long-range PCR khuếch đại toàn bộ các vùng exon và vùng intron kế cận của 7 gen đích liên quan đến khiếm thính di truyền (GJB2, GJB3, GJB6, KCNQ4, MYO7A, MYO15A, SLC26A4).
  • Vận hành hoàn chỉnh quy trình giải trình tự thế hệ mới trên hệ thống Illumina MiniSeq với các chỉ số kỹ thuật xuất sắc gồm Q30 đạt 89,9%, mật độ cụm 219 k/mm², độ phủ dọc trên 100X và độ chính xác phân tích đạt 99,9%.
  • Chứng minh độ tin cậy tuyệt đối của quy trình thông qua việc phát hiện chính xác đột biến gây bệnh c.101T>C trên gen GJB2 ở mẫu chuẩn dương NA23835 và không ghi nhận bất kỳ trường hợp dương tính giả nào trên các mẫu âm tính.
  • Thiết lập bộ thông số kiểm soát chất lượng và quy trình xử lý dữ liệu tin sinh học chuẩn mực từ file FASTQ thô đến bước phân loại biến thể dựa trên các ngân hàng dữ liệu di truyền quốc tế.
  • Công trình là tiền đề khoa học vững chắc cho việc phát triển các xét nghiệm di truyền sàng lọc khiếm thính sơ sinh chi phí thấp tại Việt Nam trong lộ trình 12 đến 24 tháng tới.

Các cơ sở y tế, viện nghiên cứu và phòng xét nghiệm di truyền có thể ứng dụng ngay quy trình kỹ thuật này vào công tác nghiên cứu và triển khai dịch vụ chẩn đoán phân tử, góp phần nâng cao chất lượng dân số và phát hiện sớm khuyết tật thính giác trong cộng đồng.