II. Phương pháp nghiên cứu
2.5. Phương pháp khảo sát độ đặc hiệu của mồi –mẫu dò
Độ đặc hiệu cho biết khả năng bắt cặp và khuếch đại của mồi và mẫu dò vào trình tự đích mong muốn mà không bắt cặp với trình tự DNA của loài khác có trong mẫu, cụ thể trình tự đích ở đây là ORF 28 của VZV.
Cách thực hiện:
Dùng hệ mồi và mẫu dò đã thiết kế chạy phản ứng Real-time PCR với các mẫu DNA của CMV, HSV, HBV, EBV, HPV, MTB (đây là những chủng có khả năng xuất hiện trong mẫu bệnh phẩm cùng với VZV), chứng dương (DNA của VZV) và một chứng âm (nước cất vô trùng). Đồng thời đặt phản ứng với hỗn hợp phản ứng Real-time PCR chuyên biệt phát hiện DNA của từng loài sử dụng mồi –mẫu dò đặc hiệu cho từng DNA của các loài CMV, HSV, HBV, EBV, HPV, MTB nhằm mục đích xác định xem hệ mồi và mẫu dò của VZV có bắt cặp với trình tự của các loài trên hay không.
CHƯƠNG III: KẾT QUẢ VÀ THẢO LUẬN 3.1. Kết quả thiết kế mồi và mẫu dò trên lý thuyết
Sau khi thực hiện các bước được nêu ở phần phương pháp. Kết quả cho thấy có vùng bảo tồn cao ở gene ORF28. Chúng tôi tiến hành download 18 trình tự trình tự ORF28 của VZV từ trang web: http://www.ncbi.nlm.nih.gov/ (phụ lục 1)
3.1.1. Kết quả thiết kế hệ mồi và mẫu dò
Hình 3.1: Kết quả sắp gióng các trình tự của ORF28 trên BioEdit
Sau khi sắp gióng các trình tự, chúng tôi tiến hành thiết kế mới hoàn toàn mồi và mẫu dòcho ORF28 dựa trên các đoạn bảo tồn của chúng. Kết quả như sau:
Bảng 3.1: Trình tự mồi và mẫu dò trên gene ORF28 của VZV
Tên Trình tự 5’-3’ Chiều dài
đoạn khuếch đại
F-VZV GCAGATTATCCGACATGCAGTCAA 164 bp
P-VZV ATTGGCGACTTGCCGGTCGAACTCG
Kết quả chạy clustal hệ mồi và mẫu dò vừa thiết kế:
• Với mồi xuôi:
Hình 3.2: Kết quả của việc tìm độ tương đồng trên gene ORF28- VZV và mồi xuôi bằng chương trình clustal.
Nhận xét: Tỷ lệ tương đồng giữa mồi xuôi và các đoạn ORF28 -VZV là 100%.
• Với mồi ngược:
Hình 3.3: Kết quả của việc tìm độ tương đồng trên gene ORF28 và mồi ngược bằng chương trình clustal.
Nhận xét: Tỷ lệ tương đồng giữa mồi ngược và các đoạn ORF28-VZV cũng là 100%.
• Với mẫu dò:
Hình 3.4 : Kết quả của việc tìm độ tương đồng trên gene ORF28 và mẫu dò bằng chương trình clustal.
Nhận xét: Tỉ lệ tương đồng của mẫu dò với 18 đoạn gene ORF28 là 100%.
Kết quả chạy Annhyb cho hệ mồi – mẫu dò:
Đối với GENE DQ479957.1: trình tự mồi xuôi, mẫu dò, và mồi ngược được thể hiện như sau:
Nhận xét: qua kết quả thu được có thể khằng định trên lý thuyết hệ mồi và mẫu dò của gene VZV có độ đặc hiệu cao với trình tự đích.
3.1.2. Kết quả khảo sát các đặc tính của hệ mồi – mẫu dò
Sau khi khảo sát các đặc tính của hệ hệ mồi và mẫu dò trên phần mềm trực tuyến Oligo Analyzer 3.0 (IDT, Hoa Kỳ), chúng tôi thu được các kết quả ở bảng
3.2 và phụ lục 3, 4, 5.
Bảng 3.2: Các thông số của hệ mồi và mẫu dò của gene ORF28
Dài (nu) %GC Tm (0C) Hairpin ΔG (kcal/mole) Self-dimer
ΔG(kcal/mole) ΔG(kcal/mole) Hetero-dimer F- VZV 24 45,8 54,7 -1,69 -7,05 Với P:-6,68 Với R:-3,61 P- VZV 25 60 65,2 -3,56 Tm: 53,9 -9,75 Với R:-6,75 R- VZV 19 57,9 56,1 -2,27 Tm: 46 -6,97
Dựa trên kết quả khảo sát đặc tính mồi nhận thấy:
Về mặt kích thước, cả hệ mồi và mẫu dò đạt yêu cầu về chiều dài tốt nhất (18 – 30bp).
Về thành phần %GC, các mồi trên có %GC nằm trong tỉ lệ tốt nhất (40% - 60%).
Nhiệt độ lai của hệ mồi và mẫu dò thoả mãn điều kiện thiết kế [Taqman probe có nhiệt độ nóng chảy cao hơn mồi 5-100
C,( Vì ở giai đọan 600
C là nhiệt độ bắt cặp tối ưu của Taqman probe, probe sẽ bắt cặp vào sợi đích trước, rồi mồi mới bắt cặp)].
Cấu trúc thứ cấp (thể hiện ở phần phụ lục): bao gồm hairpin loop (cấu trúc kẹp tóc), self dimer (oligonucleotide tự bắt cặp với chính nó), hetero dimer (sự bắt cặp giữa hai oligonucleotide khác nhau)]. Trong phản ứng PCR, nếu cấu
trúc bậc hai này càng bền vững bao nhiêu thì hiệu quả nhân bản càng giảm đi bấy nhiêu.
Độ bền vững của các cấu trúc thứ cấp được đo bằng chỉ số năng lượng tự do (∆G, kcal/mole), chỉ số này phải lớn hơn -9kcal (theo Livak,1999) vì ∆G càng lớn bao nhiêu thì cấu trúc càng kém bền bấy nhiêu, do đó sẽ không ảnh hưởng đến hoạt động của PCR.
Kết quả khảo sát trên phần mềm trực tuyến IDT oligo analyzer cho thấy, tất cả cấu trúc thứ cấp của hệ hệ mồi và mẫu dò đều có ∆G > -9kcal, chỉ trừ mẫu dò là có ∆G < -9kcal, nhưng khoảng chênh lệch này không đáng kể và nằm trong khoảng cho phép có thể chấp nhận được.
3.1.3. Kết quả khảo sát độ đặc hiệu của hệ mồi – mẫu dò
Song song với việc kiểm tra các đặc tính của mồi – mẫu dò, chúng tôi cũng tiến hành khảo sát độ đặc hiệu của mồi – mẫu dò trên lý thuyết. Để khảo sát, chúng tôi tiến hành BLAST từng trình tự trên NCBI. Chương trình này cho phép tìm kiếm các trình tự tương đồng với oligonucleotide đã có từ toàn bộ ngân hàng gene nổi tiếng trên thế giới như GENBANK, EMBL, DDBJ. Phân tích kết quả nhận được có thể kết luận cặp mồi – mẫu dò có đặc hiệu với VZV hay không. (Kết quả ở bảng phụ lục 2)
Kết quả BLAST hệ mồi – mẫu dò cho thấy mồi xuôi, mẫu dò, mồi ngược có trình tự tương đồng rất cao với duy nhất VZV. Điều này chứng minh rằng hệ mồi – mẫu dò chúng tôi sử dụng chỉ bắt cặp đặc hiệu với bộ gene của VZV.
3.2. Kết quả khảo sát hệ mồi –mẫu dò trên thực nghiệm 3.2.1. Khảo sát khả năng khuếch đại của hệ mồi –mẫu dò 3.2.1. Khảo sát khả năng khuếch đại của hệ mồi –mẫu dò
Để biết được hệ mồi –mẫu dò chúng tôi thiết kế có khả năng khuếch đại phát tín hiệu huỳnh quang với mẫu VZV trong phản ứng real-time PCR hay không. Chúng tôi tiến hành khảo sát như sau: sử dụng một mẫu chứng dương VZV và một mẫu chứng âm (với những thành phần phản ứng không đổi nhưng thay nguồn DNA bằng nước cất vô trùng).
Cho lần lượt 5µl từ mỗi mẫu trên vào mỗi mix riêng biệt (mix có chứa thành phần của hệ mồi –mẫu dò chúng tôi sử dụng). Ly tâm nhẹ và chạy phản ứng Real-time PCR.
Thí nghiệm được lặp lại ba lần với kết quả thu được:
Mẫu thử Chu kỳ ngưỡng (Ct)
Lần 1 Lần 2 Lần 3 TB
Mẫu chứng dương 25,39 25,77 24,41 25,19
Mẫu chứng âm N/A N/A N/A N/A
Kết quả minh họa bằng biểu đồ như sau:
Hình 3.6: Kết quả khảo sát khả năng khuếch đại của hệ mồi –mẫu dò (lần 2)
Nhận xét:
Biểu đồ cho thấy hệ mồi và mẫu dò đã hoạt động, đường biểu diễn của chứng dương bắt đầu ở chu kỳ trung bình là 25,19. Đường biểu diễn chứng âm không vượt qua tín hiệu nền, đường cong PCR lên đều, đẹp và không có hiện tượng nhiễu tín hiệu huỳnh quang. Chứng tỏ thí nghiệm không có hiện tượng ngoại nhiễm xảy ra.
Kết luận: Hệ mồi –mẫu dò được thiết kế có khuếch đại bản mẫu. Vì vậy, chúng ta tiến hành các khảo sát tiếp theo trên thực nghiệm với hệ mồi – mẫu dò này.
Đường nền Chứng dương
Tiếp theo, chúng tôi kiểm tra tính đặc hiệu của sản phẩm PCR nhận được với hệ mồi – mẫu dò bằng kỹ thuật giải trình tự. Mồi sử dụng cho giải trình tự là mồi đặc hiệu cho VZV đã thiết kế. Chúng tôi lấy sản phẩm PCR từ chứng dương để giải trình tự (mẫu mã số 5), gởi đến công ty Macrogenn (Hàn quốc) nhằm để kiểm tra tính đặc hiệu của sản phẩm PCR khuếch đại. Kết quả được trình bày qua bảng 3.3 và phụ lục 6.
Bảng 3.3: Kết quả giải trình tự gene ORF 28
Tên gene Trình tự nucleotid (5’ – 3’) Kích thước giải được (bp) Kích thước thiết kế (bp) ORF 28 TATCCGACATGCAGTCAATTTCAACGTCGC TTAACGTTAATTGGCGACTTGCCGGTCGAA CTCGAACACGTTCCCCATCAACTCCAGGTT TTAGTTGATACCAACCAAAA CT AAC AAA GCCGGG ATT ATCCAT TAGA AAACGA GTGGT AGCG TCT ACC
158 164
Qua kết quả bảng 3.3 cho thấy kích thước giải trình tự thu được của sản phẩm thấp hơn so với kích thước thiết kế là 6bp. Như vậy, sự chênh lệch này có thể do trong kỹ thuật giải trình tự một số nucleotide đầu tiên trong quá trình giải thường bị nhiễu tín hiệu nên không xác định được. Do vậy, sản phẩm giải trình tự sẽ cho kích thước nhỏ hơn so với sản phẩm trong thiết kế.
Kết quả giải trình tự được kiểm tra bằng phần mềm ClustalX 1.81 cho thấy: trình tự nucleotide giải được của gene ORF28 nằm trong vùng thiết kế khuếch đại sản phẩm và thuộc gene ORF28 (vị trí từ nuleotie 2592 tới 2755) của VZV(phụ lục 6).
Ngoài ra, nhằm đánh giá mức độ tương đồng của sản phẩm khuếch đại với các trình tự gene được công bố, đề tài sử dụng phần mềm BLAST để tìm
kiếm các trình tự gene được công bố trên NCBI tương đồng với sản phẩm giải trình tự. Kết quả cho thấy sản phẩm khuếch đại tương đồng với các trình tự gene ORF28 của VZV với mức độ tương đồng là 100% (phụ lục 6).
Kết luận: Như vậy, có thể khẳng định hệ mồi – mẫu dò hoạt động tốt, đặc hiệu và chuyên biệt, có thể ứng dụng được vào quy trình.
3.2.2. Khảo sát nhiệt độ lai của hệ mồi – mẫu dò
Nhiệt độ lai là một chỉ tiêu quan trọng ảnh hưởng lớn đến kết quả của phản ứng Real-time PCR. Nhiệt độ lai quá cao hoặc quá thấp đều ảnh hưởng không tốt đến kết quả phản ứng Real-rime PCR. Nếu nhiệt độ lai quá thấp so với nhiệt độ biến tính của mồi và mẫu dòsử dụng thì dễ tạo ra các sản phẩm kí sinh không mong muốn. Nếu nhiệt độ lai quá cao sẽ cản trở sự bắt cặp của mồi và mẫu dòvào DNA bản mẫu. Do đó, hiệu quả nhân bản sẽ không cao.
Từ nhiệt độ biến tính của các mồi và mẫu dò tiến hành khảo sát nhiệt độ lai cho phản ứng. Thông thường nhiệt độ lai so với nhiệt độ biến tính của mồi sai khác ±5. Do đó, ta dùng mồi và mẫu dò đã khảo sát khả năng hoạt động ở trên để xác định nhiệt độ lai, chạy trên 8 mix chứng dương, một chứng âm với gradient nhiệt độ được xác định bằng cách lấy nhiệt độ trung bình giữa mồi có nhiệt độ thấp nhất và mồi có nhiệt độ cao nhất ±5 (Ta = (Tm1 + Tm2)/2 ±50
C). Từ đó tính được khoảng gradient nhiệt độ từ 550 đến 650
C trong khi các thành phần (DNA mẫu, mồi, mẫu dò,…) và các điều kiện khác (điều kiện về thiết bị, thời gian phản ứng…) đều được giữ nguyên và đồng nhất.
Bảng 3.4: Kết quả khảo sát nhiệt độ lai Giếng phản ứng Nhiệt độ lai (0C) Chu kỳ ngưỡng (Ct) (Ct) TB lần 1 lần 2 lần 3 A08 650C 26,56 24,26 26,19 25,67 B08 64,50C 25,51 25,05 26,65 25,74 C08 63,30C 25,09 25,16 25,19 25,15 D08 61,40C 26,09 24,51 25,50 25,36 E08 58,90C 25,51 24,42 25,58 25,17 F08 57,10C 25,09 23,85 25,37 24,77 G08 55,80C 25,49 23,56 25,18 24,74 H08 55,00C 23,04 23,01 23,13 23,06
(-) (550C) N/A N/A N/A N/A
Hình 3.7: Kết quả khảo sát nhiệt độ lai (lần 1)
550C 55,80C 61,40C 650C Chứng âm 57,10C 63,30C 57,10C 63,30C
Nhận xét:
Từ kết quả bảng 3.4 cho thấy ở tất cả các nhiệt độ điều cho tín hiệu dương tính và nhiệt độ 550
C có giá trị Ct thấp nhất nghĩa là ở nhiệt độ 550
C cho hiệu quả nhân bản sớm nhất so với các nhiệt độ khác.
Kết luận: Mồi hoạt động tốt ở nhiệt độ lai là 550
C.
3.2.3. Khảo sát nồng độ mẫu dò
Nồng độ mẫu dò cũng là một yếu tố quan trọng trong phản ứng Real-time PCR. Do đó, cần có một nồng độ mẫu dò nhất định để sử dụng vừa đủ trong quá trình thực hiện phản ứng, tránh trường hợp thiếu hoặc thừa mẫu dò, như vậy sẽ không đem lại kết quả như mong muốn, đồng thời tốn nhiều chi phí. Vì thế, chúng tôi tiến hành bước khảo sát này nhằm tìm ra nồng độ mẫu dò thích hợp. Theo các nghiên cứu, nồng độ mẫu dò trong phản ứng Real-time PCR dao động trong khoảng 100nM đến 500nM. Chúng tôi chọn phổ nồng độ từ 100-500nM với đơn vị tăng là 100nM nhằm tìm ra nồng độ cho giá trị Ct nhỏ nhất, thí nghiệm được lặp lại ba lần và kết quả được trình bày ở bảng 3.5.
Bảng 3.5: Kết quả khảo sát nồng độ mẫu dò
Nồng độ mẫu dò (nM) Chu kỳ ngưỡng (Ct) Ct TB Lần 1 Lần 2 Lần 3 100 25,79 25,88 25,84 25,84 200 25,06 25,45 25,18 25,23 300 24,72 24,87 24,84 24,81 400 26,42 26,45 26,42 26,43 500 26,58 26,53 26,57 26,56
Hình 3.8: Kết quả khảo sát nồng độ mẫu dò (lần 2)
Nhận xét:
Kết quả bảng 3.5 cho thấy nồng độ mẫu dò 300nM cho giá trị Ct nhỏ nhất (24,81). Nồng độ này được sử dụng cho nghiên cứu tiếp theo.
Kết luận: Mẫu dò hoạt động tốt ở nồng độ 300nM.
3.2.4. Khảo sát nồng độ Mg2+
MgCl2 là đồng yếu tố (co-factor) cần thiết cho hoạt động tổng hợp của hầu hết các DNA polymerase, trong đó có Taq polymerase. Nồng độ MgCl2 ảnh hưởng đến hiệu quả và tính đặc hiệu của phản ứng. Nếu nồng độ MgCl2 quá thấp làm giảm hoặc mất hẳn hiệu quả nhân bản. Ngược lại, nếu nồng độ MgCl2 quá cao sẽ ức chế phản ứng do tăng tính bền vững cấu trúc mạch đôi của DNA, ổn định cấu trúc bắt cặp không đặc hiệu của mồi, làm giảm tính đặc hiệu và tạo sản phẩm kí sinh.
Do đó, chúng tôi khảo sát nồng độ MgCl2 tìm nồng độ MgCl2 của phản ứng có chu kì ngưỡng Ct nhỏ nhất, nghĩa là phản ứng có độ nhạy cao nhất. Nồng độ cần khảo sát từ 1,5 đến 5mM, với độ tăng là 0,5 mM và chứng âm.
Khi tiến hành khảo sát nồng độ MgCl2, chúng tôi tiến hành các phản ứng ở cùng nhiệt độ lai 550C, thành phần các chất tham gia phản ứng là như nhau, chỉ khác biệt về nồng độ MgCl2. Chúng tôi lặp lại thí nghiệm ba lần và kết quả thu được bảng 3.6. 300nM 500nM 200nM 400nM 100nM
Bảng 3.6: Kết quả khảo sát nồng độ MgCl2 Nồng độ MgCl2 (mM) (Ct) lần 1 (Ct) lần 2 (Ct) lần 3 (Ct) TB 1,5 (+) / 25,14 23,57 / 2,0 26,89 26,64 23,96 25,83 2,5 27,47 26,06 27,54 27,02 3,0 26,75 26,31 25,52 26,19 3,5 24,66 24,20 24,89 24,58 4,0 27,51 27,37 27,71 27,53 4,5 26,18 27,27 25,62 26,36 5,0 28,08 28,16 25,52 27,25 1,5 (-) Âm tính Âm tính Âm tính
Hình 3.9: Kết quả khảo sát nồng độ MgCl2 (lần 2)
Nhận xét:
Qua kết quả bảng 3.6 cho thấy nồng độ MgCl2 = 3,5mM ở chu kỳ ngưỡng thấp nhất (Ct = 24,58) đã có tín hiệu huỳnh quang.
Kết luận: Vậy chúng tôi chọn nồng độ MgCl2= 3,5mM là tối ưu cho phản ứng. MgCl2=3,0 MgCl2=5,0 MgCl2=1,5(-) MgCl2=2,0 MgCl2=2,5 MgCl2=4,5 MgCl2=3,5 MgCl2=1,5(+) MgCl2=4,0
3.2.5. Khảo sát sát khả năng nhân bản chọn lọc của hệ mồi và mẫu dò
Trên lý thuyết, bằng phần mềm Annhyb 4,942 và trình tự genome của các tác nhân HSV1, HSV2, CMV, EBV, HBV, HPV, vi khẩn lao từ Genbank để khảo sát khả năng bắt cặp của hệ mồi – mẫu dò. Kết quả khảo sát bảng 3.7.
Bảng 3.7: Kết quả khảo sát độ đặc hiệu của hệ mồi – mẫu dò thiết kế
Tên chủng HSV1 HSV2 CMV EBV HBV MTB HPV Phần trăm bắt cặp MỒI XUÔI 46,6% 25,8% 52,5% 35% 42,9% 55,4% 42,5% MỒI NGƯỢC 53,6% 57,3% 59% 33,6% 45% 74,7% 48,4% MẪU DÒ 56,8% 33,2% 50% 49% 43,7% 69,2% 33,6%
Qua kết quả bảng 3.7, chúng tôi thấy phần trăm bắt cặp giữa hệ mồi – mẫu dò của VZV trên genome của các loài HSV-1, HSV-2, CMV, EBV, HBV, HCV, HPV, MTB đều không cao, đối với mồi xuôi dao động trong khoảng 25,8% đến 55,4%, với mồi ngược trong khoảng 33,6%-74,7%, còn với mẫu dò nằm trong khoảng 33,2%-69,2% khẳng định hệ mồi-mẫu dò được thiết kế là đặc hiệu cho VZV
Trên thực nghiệm, chúng tôi tiến hành song song hai phản ứng Real-time PCR với các chứng dương CMV, HPV, HSV, MTB, HBV,EBV và một mẫu chứng dương với VZV. Và một mẫu chứng âm thực hiện trên phản ứng mồi – mẫu dò của VZV (đồng thời các mẫu chứng dương sẽ thực hiện phản ứng trên hệ