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

Ong là nhóm côn trùng đơn ngành (monophyletic) thuộc bộ Cánh màng (Hymenoptera), lớp Côn trùng (Insecta), với hơn 19.000 loài đã được ghi nhận trên toàn thế giới. Trong 7 họ ong hiện diện (bao gồm Andrenidae, Apidae, Colletidae, Halictidae, Megachilidae, Melittidae và Stenotritidae), họ Ong mật (Apidae) là nhóm lớn và đa dạng bậc nhất với hơn 5.700 loài, chia làm 3 phân họ (Apinae, Xylocopinae, Nomadinae). Trong phân họ Apinae, tông ong Dú (Meliponini - stingless bees) là nhóm ong không ngòi đốt phân bố chủ yếu ở các vùng nhiệt đới và cận nhiệt đới như Châu Phi, Châu Úc, Trung và Nam Mỹ, với khoảng 54 chi khác nhau (theo hệ thống phân loại của NCBI tính đến năm 2023).

Ong Dú đóng vai trò sinh thái quan trọng trong việc thụ phấn cây trồng và mang giá trị kinh tế - y học cao từ mật ong và keo ong có hoạt tính sinh học. Tuy nhiên, các loài ong Dú thường có kích thước rất nhỏ (khoảng 2 mm), đặc điểm hình thái tương đồng giữa các chi, dễ dẫn đến nhầm lẫn khi phân loại chỉ dựa trên khóa hình thái truyền thống. Việc kết hợp chỉ thị sinh học phân tử và phân tích phát sinh loài là hướng tiếp cận cần thiết để định danh chính xác các loài ong Dú tại Việt Nam.

Khóa luận tốt nghiệp "Xây dựng mô hình phân tích phả hệ phân tử nhằm hỗ trợ định danh một số loài ong Dú ở Việt Nam" do sinh viên Phạm Văn Nam thực hiện (chuyên ngành Sinh học Phân tử, Khoa Công nghệ Sinh học, Trường Đại học Mở TP. Hồ Chí Minh, dưới sự hướng dẫn của TS. Trương Kim Phượng) nhằm giải quyết các mục tiêu và nhiệm vụ cụ thể:

  1. Mục tiêu tổng quát: Áp dụng phương pháp PCR kết hợp giải trình tự Sanger và các công cụ tin sinh học để xây dựng mô hình phân tích phả hệ phân tử, hỗ trợ định danh một số mẫu ong Dú có giá trị dược tính tại Việt Nam dựa trên hai gen mục tiêu là COI và EF-1α.
  2. Nhiệm vụ nghiên cứu:
    • Khai thác và hệ thống hóa dữ liệu trình tự gen COI và EF-1α của các loài thuộc tông Meliponini từ GenBank (NCBI).
    • Khảo sát in silico để kiểm tra, đánh giá tính đặc hiệu và thông số kỹ thuật của các cặp mồi chuyên biệt cho vùng gen COI và EF-1α.
    • Thực nghiệm tách chiết DNA tổng số, tối ưu hóa điều kiện phản ứng PCR và giải trình tự nucleotide của các sản phẩm khuếch đại.
    • Phân tích độ tương đồng trình tự (BLAST) và tái dựng cây phát sinh loài phân tử bằng phần mềm MEGA X để xác định vị trí phân loại của các mẫu nghiên cứu.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: 7 mẫu ong Dú do Phòng Công nghệ Sinh học Vật liệu và Nano - Trung tâm Công nghệ Sinh học TP. Hồ Chí Minh cung cấp (đã được phân loại sơ bộ qua hình thái).
  • Phạm vi nghiên cứu: Tập trung vào vùng chỉ thị gen ti thể COI (Cytochrome c oxidase subunit I) và gen nhân EF-1α (elongation factor-1α), kết hợp dữ liệu trình tự đối chứng và nhóm ngoại (outgroup) được công bố trên ngân hàng gen quốc tế GenBank.

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

Khung lý thuyết và chỉ thị phân tử

Nghiên cứu ứng dụng hai hệ thống chỉ thị phân tử bổ trợ lẫn nhau:

  • Gen ti thể COI (Cytochrome c oxidase subunit I): Gen bảo tồn trong hệ gen ti thể, đóng vai trò tạo trung tâm oxy hóa - khử của enzyme cytochrome c oxidase trong chuỗi hô hấp tế bào. Vùng gen này có tốc độ tiến hóa thích hợp, được sử dụng phổ biến làm mã vạch DNA (DNA barcoding) và phân tích phả hệ côn trùng.
  • Gen nhân EF-1α (elongation factor-1α): Gen mã hóa yếu tố kéo dài chuỗi polypeptide trong dịch mã, có độ bảo tồn cao ở mức độ nucleotide và amino acid, đặc biệt hữu ích trong việc phân tích các mối quan hệ phát sinh loài ở cấp độ chi và loài trong họ Apidae.

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

Nghiên cứu phối hợp chặt chẽ giữa phân tích tin sinh học (in silico) và thực nghiệm sinh học phân tử (in vitro):

Khai thác dữ liệu GenBank ➔ Thiết kế & Đánh giá mồi in silico (BLAST, IDT, AnnHyb)
       ➔ Tách chiết DNA mẫu ong (Phenol:Chloroform) ➔ Đo OD A260/A280
       ➔ Phản ứng PCR & Điện di gel Agarose ➔ Giải trình tự Sanger 2 chiều
       ➔ Xử lý chuỗi (BioEdit/Chromas) ➔ Tìm mô hình tiến hóa (ModelTest) 
       ➔ Tái dựng cây phả hệ (ML, MP, NJ trong MEGA X với Bootstrap 1.000)
  1. Thu thập và xử lý dữ liệu tin sinh học:

    • Thu thập bộ trình tự tham chiếu và nhóm ngoại từ NCBI/GenBank.
    • Đánh giá thông số lý - hóa của các mồi (nhiệt độ nóng chảy $T_m$, %GC, năng lượng tự bắt cặp, cấu trúc kẹp tóc) bằng phần mềm trực tuyến IDT OligoAnalyzer.
    • Khảo sát độ đặc hiệu mồi bằng công cụ NCBI BLAST và kiểm tra vị trí bắt cặp thực tế trên chuỗi đích qua phần mềm AnnHyb.
  2. Quy trình thực nghiệm phòng thí nghiệm:

    • Tách chiết DNA genome: Áp dụng phương pháp ly giải tế bào với SDS/Proteinase K kết hợp tinh sạch bằng Phenol : Chloroform (tỷ lệ 1:1) và kết tủa bằng ethanol tuyệt đối.
    • Đánh giá chất lượng DNA: Đo mật độ quang phổ kế UV-Vis tại bước sóng 260 nm và 280 nm ($A_{260}/A_{280}$) để xác định độ tinh sạch và nồng độ DNA.
    • Khuếch đại PCR: Thiết lập phản ứng PCR thể tích 50 µL (gồm 25 µL Master mix 2X, 1 µL mồi F, 1 µL mồi R, 1 µL DNA khuôn và 22 µL $H_2O$) với chu trình nhiệt tối ưu hóa.
    • Điện di và giải trình tự: Kiểm tra sản phẩm PCR bằng điện di trên gel agarose; gửi mẫu giải trình tự Sanger hai chiều; xử lý và kiểm tra tín hiệu sắc ký đồ nucleotide bằng BioEdit và Chromas Lite.
  3. Phương pháp phân tích phát sinh loài phân tử:

    • Dóng hàng đa trình tự nucleotide bằng thuật toán tích hợp trong phần mềm MEGA X.
    • Xác định mô hình tiến hóa thay thế nucleotide tối ưu thông qua công cụ ModelTest trong MEGA X.
    • Xây dựng cây phát sinh loài theo 3 phương pháp độc lập: Cực đại khả dĩ (Maximum Likelihood - ML), Tiết kiệm tối đa (Maximum Parsimony - MP) và Gia nhập láng giềng (Neighbor-Joining - NJ).
    • Kiểm định độ tin cậy của các nút nhánh trên cây phả hệ bằng phân tích tái lấy mẫu Bootstrap với 1.000 lần lặp.

Nội dung chính theo từng phần của khóa luận

Phần I: Mở đầu và Tổng quan tài liệu

Tác giả trình bày hệ thống hóa các kiến thức cơ sở về cấu trúc cơ thể côn trùng (chia làm 3 phần: đầu, ngực, bụng; đặc điểm cánh màng và các biến thái hoàn toàn), đặc điểm sinh học và đời sống xã hội của họ Ong mật (Apidae) và tông Ong Dú (Meliponini). Phần này làm rõ đặc điểm làm tổ bằng gôm, nhựa, sáp và cấu trúc không ngòi đốt của ong Dú. Đồng thời, tác giả tổng quan vai trò của hai vùng chỉ thị phân tử COI và EF-1α trong các nghiên cứu phân loại sinh học côn trùng trên thế giới.

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

Trình bày chi tiết nguồn gốc 7 mẫu ong Dú thu thập tại Việt Nam, quy trình 7 bước tách chiết DNA bằng phenol:chloroform, công thức hóa chất đệm ly giải (Tris-HCl 1M, EDTA 0,5M, NaCl 5M, SDS 10%, Proteinase K 20 mg/mL). Tác giả mô tả chi tiết quy trình chạy PCR, nguyên lý phản ứng giải trình tự chuỗi bằng phương pháp kết thúc chuỗi dideoxynucleotide (Sanger) và các bước dựng ma trận khoảng cách di truyền phục vụ phân tích phả hệ trên phần mềm MEGA X.

Phần III: Kết quả và Thảo luận

1. Xây dựng bộ dữ liệu tham chiếu

Từ cơ sở dữ liệu NCBI/GenBank, tác giả đã thu thập hai bộ dữ liệu chuẩn:

  • Dữ liệu gen COI: Gồm 18 trình tự, trong đó có 16 trình tự thuộc các chi ong Dú (Partamona, Melipona, Plebeia, Liotrigona, Trigona, Hypotrigona, Frieseomelitta, Geniotrigona, Tetragonula) và 2 trình tự nhóm ngoại thuộc họ Kiến Formicidae (Leptomyrmex cnemidatus HQ207338.1, Leptomyrmex erythrocephalus HQ207346.1).
  • Dữ liệu gen EF-1α: Gồm 14 trình tự, trong đó có 12 trình tự tông Meliponini (Lepidotrigona, Scaptotrigona, Hypotrigona, Plebeia, Melipona, Geniotrigona, Tetragonula) và 2 trình tự nhóm ngoại thuộc tông Ancylini (Baeopelma foersteri KY983275.1) và chi Apis (Acizzia albizziae KY983277.1).

2. Kết quả khảo sát bộ mồi in silico

Tác giả đánh giá các cặp mồi chuyên biệt qua công cụ tin sinh học:

  • Cặp mồi gen COI (CB2-F và CB2-R):
    • Trình tự: CB2-F (5’-GGATCTCCTGATATAGCATTTCCACG-3’) và CB2-R (3’-GCTCCAGCTAATACTGGTAATG-5’).
    • Kích thước khuếch đại dự kiến: 353 bp; $T_m$ đo được: 56,54 °C.
    • Phân tích BLAST và AnnHyb cho thấy mồi bắt cặp đặc hiệu với độ tương đồng > 90% trên các trình tự Meliponini đối chứng (Partamona bilineata đạt 100%, Mourella caerulea đạt 91,36% - 92,69%, Liotrigona mahafalya đạt 85,38% - 91,36%).
  • Cặp mồi gen EF-1α (EB1-F và EB1-R):
    • Trình tự: EB1-F (5’-CTACTGAACCACCATACTCCGAGACC-3’) và EB1-R (5’-GACCAACTGGTACTGTTCCAATACC-3’).
    • Kích thước khuếch đại dự kiến: 244 bp; $T_m$: 59,49 °C.
    • BLAST cho độ tương đồng 100% với Melipona compressipes (EU163224); phân tích AnnHyb trên Mourella caerulea, Tetragonisca angustula và Oxytrigona mediorufa cho thấy mồi xuôi đạt độ tương đồng 85,38% - 92%, mồi ngược đạt 100%.

3. Kết quả thực nghiệm tách chiết và phản ứng PCR

  • Quy trình tách chiết cho nồng độ và độ tinh sạch DNA tốt. Tỷ số quang phổ $A_{260}/A_{280}$ của các mẫu dao động trong ngưỡng chuẩn từ 1,8 đến 2,0 (ví dụ mẫu 1 đạt tỷ lệ $A_{260}/A_{280} = 1,8$ với nồng độ DNA đạt 12,8 ng/µL).
  • Điện di sản phẩm PCR cho thấy:
    • Với gen COI: Hai mẫu số 1 và số 5 khuếch đại vạch sáng rõ, được ký hiệu lần lượt là C8 và C12.
    • Với gen EF-1α: Bốn mẫu cho sản phẩm khuếch đại rõ nét, được ký hiệu lần lượt là E3, E4, E6 và E7.
  • Chu trình nhiệt PCR chuẩn hóa thành công: Biến tính ban đầu ở 95 °C trong 5 phút; 40 chu kỳ lặp lại gồm [95 °C trong 30 giây, bắt cặp ở 56 °C trong 30 giây, kéo dài ở 72 °C trong 120 giây]; kéo dài kết thúc ở 72 °C trong 7 phút.

4. Kết quả giải trình tự và phân tích độ tương đồng (BLAST)

Tất cả các phân đoạn PCR chọn lọc đều được giải trình tự Sanger 2 chiều thành công:

  • Mẫu gen COI:
    • Mẫu C8: BLAST cho độ tương đồng 83% - 87% (độ bao phủ 95% - 100%) với các loài Tetragonula fuscobalteata, Tetragonula iridipennis, Mourella caerulea.
    • Mẫu C12: BLAST cho độ tương đồng 84% - 87% (độ bao phủ 95% - 99%) với Tetragonula fuscobalteata, Frieseomelitta nigra, Trigona chanchamayoensis.
  • Mẫu gen EF-1α:
    • Mẫu E3, E6, E7: BLAST cho độ tương đồng cao từ 94% đến 99% (độ bao phủ 94% - 99%) với các loài Lepidotrigona terminata, Lepidotrigona ventralis, Heterotrigona itama.
    • Mẫu E4: BLAST cho độ tương đồng 98% - 99% (độ bao phủ ~93%) với Tetragonula sirindhornae, Tetragonula hirashimai, Tetragonula sapiens.

5. Phân tích mô hình tiến hóa và xây dựng cây phát sinh loài

  • Kiểm tra qua công cụ ModelTest xác định mô hình tiến hóa thay thế nucleotide tối ưu cho cả hai tập dữ liệu COI và EF-1α là GTR + G + I (General Time Reversible với phân phối biến thiên tốc độ Gamma và tỷ lệ vị trí bất biến Invariant).
  • Tái dựng cây phả hệ bằng 3 phương pháp (ML, MP, NJ):
    • Địa hình học (topology) của 3 cây phả hệ trên từng gen hoàn toàn tương đồng.
    • Nhóm ngoại phân tách riêng biệt với chỉ số bootstrap cao (đạt 86% - 100% ở cây COI và 100% ở cây EF-1α).
    • Vị trí các mẫu trên cây phả hệ:
      • Mẫu C8 phân nhóm cùng nhánh với loài Tetragonula iridipennis (bootstrap > 60%).
      • Mẫu C12 nằm cùng phân nhóm với các chi Melipona và Partamona (bootstrap đạt 100%).
      • Các mẫu E3, E6, E7 gom cụm cùng chi Lepidotrigona với giá trị bootstrap đạt 85% - 88%.
      • Mẫu E4 gom cụm cùng chi Tetragonula với độ tin cậy bootstrap đạt 100%.

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

Bảng tổng hợp kết quả định danh phân tử các mẫu ong Dú nghiên cứu

Ký hiệu mẫu Gen mục tiêu Cặp mồi sử dụng Kết quả phân tích BLAST (NCBI) Kết quả cây phát sinh loài (MEGA X) Phân loại phân tử kết luận
C8 COI CB2-F / CB2-R Tương đồng 83% - 87% với Tetragonula fuscobalteata, T. iridipennis Phân nhóm cùng Tetragonula iridipennis (Bootstrap > 60%) Chi Tetragonula (gần loài T. iridipennis)
C12 COI CB2-F / CB2-R Tương đồng 84% - 87% với T. fuscobalteata, Frieseomelitta nigra Cùng phân nhóm với Melipona và Partamona (Bootstrap = 100%) Phân nhóm liên quan Melipona / Partamona
E3 EF-1α EB1-F / EB1-R Tương đồng 95% - 99% với Lepidotrigona terminata, L. ventralis Gom cụm cùng chi Lepidotrigona (Bootstrap = 88%) Chi Lepidotrigona
E4 EF-1α EB1-F / EB1-R Tương đồng 98% - 99% với Tetragonula sirindhornae, T. sapiens Gom cụm cùng chi Tetragonula (Bootstrap = 100%) Chi Tetragonula
E6 EF-1α EB1-F / EB1-R Tương đồng 94% - 99% với Lepidotrigona terminata, L. ventralis Gom cụm cùng chi Lepidotrigona (Bootstrap = 87%) Chi Lepidotrigona
E7 EF-1α EB1-F / EB1-R Tương đồng 95% - 99% với Lepidotrigona terminata, L. ventralis Gom cụm cùng chi Lepidotrigona (Bootstrap = 85%) Chi Lepidotrigona

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

  1. Thiết lập quy trình kỹ thuật hoàn chỉnh: Xây dựng quy trình chuẩn từ tách chiết DNA mô côn trùng nhỏ, thiết kế - đánh giá mồi in silico, tối ưu hóa thông số luân nhiệt PCR (56 °C cho cả hai cặp mồi) đến giải trình tự Sanger.
  2. Xác định mô hình tiến hóa chuẩn: Chứng minh tính hiệu quả của mô hình thay thế nucleotide GTR + G + I trong việc phân tích cây phát sinh loài cho các chi ong Dú họ Meliponini.
  3. Cung cấp dữ liệu định danh sinh học phân tử: Hỗ trợ định danh chính xác các mẫu ong Dú thu thập tại Việt Nam thuộc các chi Tetragonula và Lepidotrigona, khắc phục hạn chế của phương pháp hình thái học đối với các loài côn trùng có kích thước nhỏ.

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

Hạn chế

  • Số lượng mẫu thực nghiệm thu nhận được sản phẩm khuếch đại và giải trình tự thành công còn giới hạn (2 mẫu gen COI và 4 mẫu gen EF-1α từ 7 mẫu ban đầu).
  • Bộ dữ liệu tham chiếu trên GenBank đối với các loài ong Dú đặc hữu tại khu vực Đông Nam Á và Việt Nam còn chưa thực sự phong phú, dẫn đến một số mẫu (như C12) mới chỉ dừng lại ở mức độ phân nhóm chi lớn.

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

  • Mở rộng phạm vi thu thập mẫu ong Dú tại nhiều vùng sinh thái khác nhau trên cả nước.
  • Mở rộng phân tích đa gen (kết hợp đồng thời nhiều chỉ thị: COI, COII, 16S rRNA, ND5 và EF-1α) để tăng độ phân giải phân loại học và xây dựng cây phát sinh loài toàn diện hơn.
  • Kết hợp sâu hơn giữa dữ liệu mã vạch phân tử và các phân tích hình thái hiển vi chuyên sâu.

Giá trị tham khảo

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

  • Sinh viên, học viên ngành Công nghệ Sinh học, Sinh học Phân tử, Côn trùng học: Cung cấp tài liệu thực hành chi tiết về quy trình khai thác dữ liệu NCBI, sử dụng các công cụ IDT OligoAnalyzer, AnnHyb, chạy ModelTest và xây dựng cây phả hệ trên MEGA X.
  • Các phòng thí nghiệm và nhà nghiên cứu phân loại sinh học: Cung cấp thông số nhiệt, trình tự mồi chuẩn (CB2-F/R và EB1-F/R) để khuếch đại thành công vùng gen COI (353 bp) và EF-1α (244 bp) trên mẫu ong Dú kích thước nhỏ.

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

1. Cặp mồi nào được tác giả lựa chọn để khuếch đại vùng gen COI và EF-1α trong thực nghiệm?

Tác giả đã chọn lọc và kiểm chứng 2 cặp mồi:

  • Vùng gen COI: Cặp mồi CB2-F (5’-GGATCTCCTGATATAGCATTTCCACG-3’) và CB2-R (3’-GCTCCAGCTAATACTGGTAATG-5’), cho sản phẩm khuếch đại kích thước 353 bp, nhiệt độ nóng chảy $T_m = 56,54\ ^\circ\text{C}$.
  • Vùng gen EF-1α: Cặp mồi EB1-F (5’-CTACTGAACCACCATACTCCGAGACC-3’) và EB1-R (5’-GACCAACTGGTACTGTTCCAATACC-3’), cho sản phẩm khuếch đại kích thước 244 bp, nhiệt độ nóng chảy $T_m = 59,49\ ^\circ\text{C}$.

2. Điều kiện chu trình nhiệt PCR được chuẩn hóa trong khóa luận như thế nào?

Phản ứng PCR 50 µL được tiến hành theo chu trình:

  • Giai đoạn biến tính ban đầu: $95\ ^\circ\text{C}$ trong 5 phút (1 chu kỳ).
  • Giai đoạn chu kỳ (40 chu kỳ): Biến tính ở $95\ ^\circ\text{C}$ trong 30 giây; bắt cặp mồi ở $56\ ^\circ\text{C}$ trong 30 giây; kéo dài ở $72\ ^\circ\text{C}$ trong 120 giây.
  • Giai đoạn kéo dài cuối cùng: $72\ ^\circ\text{C}$ trong 7 phút (1 chu kỳ).

3. Mô hình tiến hóa nucleotide nào được lựa chọn để tái dựng cây phả hệ và trên phần mềm nào?

Thông qua công cụ ModelTest tích hợp trong phần mềm MEGA X, tác giả đã xác định mô hình tiến hóa tối ưu cho cả hai tập dữ liệu gen COI và EF-1α là mô hình GTR + G + I (General Time Reversible với phân phối biến thiên Gamma và tỷ lệ vị trí bất biến).

4. Kết quả phân tích phả hệ đã định danh các mẫu ong nghiên cứu vào những phân nhóm loài nào?

  • Mẫu C8 (COI): Nằm cùng phân nhóm với loài Tetragonula iridipennis (bootstrap > 60%).
  • Mẫu C12 (COI): Nằm cùng phân nhóm với chi Melipona và Partamona (bootstrap 100%).
  • Các mẫu E3, E6, E7 (EF-1α): Nằm cùng phân nhóm với chi Lepidotrigona (bootstrap đạt từ 85% đến 88%).
  • Mẫu E4 (EF-1α): Nằm cùng phân nhóm với chi Tetragonula (bootstrap 100%).

Kết luận

Khóa luận tốt nghiệp của sinh viên Phạm Văn Nam đã xây dựng thành công quy trình phân tích phả hệ phân tử dựa trên hai vùng chỉ thị gen COI và EF-1α nhằm hỗ trợ phân loại các loài ong Dú tại Việt Nam. Bằng việc kết hợp đánh giá mồi in silico, thực nghiệm PCR - giải trình tự Sanger và xây dựng cây phát sinh loài trên MEGA X theo mô hình GTR + G + I, nghiên cứu đã định danh thành công các mẫu ong thuộc các chi Tetragonula và Lepidotrigona. Kết quả của đề tài đóng góp cơ sở dữ liệu phân tử và quy trình kỹ thuật hữu ích phục vụ cho công tác giám định, bảo tồn và khai thác nguồn gen côn trùng có giá trị dược liệu.