Giới thiệu dự án

Nguồn gen cây có múi (Citrus spp.) tại Việt Nam sở hữu mức độ đa dạng sinh học cao với nhiều giống bản địa đặc hữu có giá trị kinh tế vượt trội như cam sành Hàm Yên (Tuyên Quang), cam Cao Phong (Hòa Bình), cam Vinh (Nghệ An), cam Bù (Hà Tĩnh) và đặc biệt là cam Bố Hạ (Bắc Giang). Theo số liệu từ Tổng cục Thống kê, diện tích trồng cam cả nước đạt 97.500 ha với sản lượng xấp xỉ 799.500 tấn (tăng trưởng 20,4% so với các năm trước đó). Tuy nhiên, năng suất trung bình và độ thuần chủng giống đang đối mặt với sự suy giảm nghiêm trọng. Riêng thương hiệu cam Bố Hạ từng là cây kinh tế mũi nhọn tại vùng Yên Thế - Bắc Giang nhưng hiện đang đứng trước nguy cơ thoái hóa, xói mòn nguồn gen và mất giống sau khi Nông trường cam Bố Hạ giải thể.

                  +-------------------------------------------------------+
                  |  THỰC TRẠNG NGUỒN GEN CAM BỐ HẠ (YÊN THẾ, BẮC GIANG)  |
                  +-------------------------------------------------------+
                                              |
                     +------------------------+------------------------+
                     |                                                 |
                     v                                                 v
        [Suy thoái di truyền & Dịch bệnh]              [Hạn chế phân loại hình thái]
         - Nông trường giải thể                         - Đồng hình kiểu hình (phenotypic plasticity)
         - Nguy cơ tuyệt chủng nguồn gen bản địa        - Khó phân biệt giống lai / đột biến sớm
                     |                                                 |
                     +------------------------+------------------------+
                                              |
                                              v
                  +-------------------------------------------------------+
                  |  GIẢI PHÁP: ĐÁNH GIÁ ĐA DẠNG BẰNG CHỈ THỊ PHÂN TỬ     |
                  |          (Hệ thống phối hợp 10 RAPD + 3 ISSR)         |
                  +-------------------------------------------------------+

Vấn đề cốt lõi nảy sinh từ việc phân loại truyền thống dựa trên đặc điểm hình thái học (morphological markers) thường bị biến đổi bởi điều kiện môi trường canh tác, khó phân biệt chính xác giữa các dòng đột biến tự nhiên hoặc dạng lai ghép cận chủng. Để giải quyết triệt để vấn đề này, việc ứng dụng chỉ thị phân tử DNA (Molecular Markers) là phương pháp chính xác nhằm định danh và xác lập quan hệ phát sinh chủng loài ở mức độ phân tử.

Đề tài đặt ra các mục tiêu cụ thể:

  1. Tách chiết và tinh sạch DNA tổng số: Thu nhận mẫu DNA đạt chuẩn quang phổ ($OD_{260}/OD_{280} = 1.8 - 2.0$) từ 32 mẫu lá cây có múi đại diện.
  2. Khuếch đại đa hình DNA: Tối ưu hóa phản ứng PCR với 10 mồi RAPD (Random Amplified Polymorphic DNA) và 3 mồi ISSR (Inter-Simple Sequence Repeat).
  3. Giải mã ma trận di truyền: Xây dựng ma trận nhị phân tương đồng di truyền (Jaccard/Dice coefficient) và thiết lập cây phát sinh chủng loài (Dendrogram) thông qua phần mềm NTSYSpc v2.1 (Numerical Taxonomy and Multivariate Analysis System).

Dự án giới hạn trên 32 mẫu giống thực vật có múi thu thập tại các vùng sinh thái trọng điểm (Bắc Giang, Tuyên Quang, Hà Giang, Thái Nguyên, Viện Nghiên cứu Rau Quả), tập trung đối sánh các phân đoạn DNA đặc trưng nhằm phục vụ công tác bảo tồn nguồn gen bản địa.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi áp dụng phương pháp phối hợp chỉ thị phân tử RAPD và ISSR, các nhà nghiên cứu nông nghiệp chủ yếu dựa vào các kỹ thuật đánh giá truyền thống hoặc các chỉ thị đơn lẻ:

Tiêu chí đánh giá Chỉ thị hình thái (Morphology) Chỉ thị RFLP / AFLP Kỹ thuật giải trình tự NGS Giải pháp phối hợp RAPD + ISSR
Độ tin cậy di truyền Thấp (chịu ảnh hưởng ngoại cảnh) Rất cao Tuyệt đối ($>99.9%$) Cao ($>90%$ độ lặp lại)
Chi phí triển khai Rất thấp Cao (hóa chất độc hại/phóng xạ) Rất cao (yêu cầu hạ tầng lớn) Tối ưu, phù hợp lab tiêu chuẩn
Yêu cầu thông tin bộ gen Không Cần một phần Bắt buộc Không cần biết trước trình tự
Thời gian phân tích $6 - 24$ tháng (theo mùa vụ) $5 - 7$ ngày $2 - 4$ tuần $24 - 48$ giờ

Phân tích yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have: Quy trình tách chiết DNA đạt chuẩn điện di không đứt gãy; tối ưu hóa nhiệt độ gắn mồi ($T_m = 32.0^\circ\text{C} - 43.7^\circ\text{C}$); xác lập ma trận 0/1 không sai sót.
  • Should have: Phối hợp song song cả mồi ngẫu nhiên (RAPD) và mồi vi vệ tinh (ISSR) để tăng hệ số đa hình ($PIC$).
  • Could have: Đối sánh với dữ liệu kiểu hình thực địa để lập bản đồ liên kết tính trạng nông học.
  • Won't have: Giải trình tự toàn bộ hệ gen (Whole Genome Sequencing) nhằm đảm bảo tính khả thi về kinh phí nghiên cứu cấp cơ sở.

Thiết kế hệ thống

Toàn bộ quy trình phân tích sinh học phân tử được cấu trúc hóa theo đường ống (pipeline) khép kín:

Technology Stack và thiết bị đo kiểm:

  • Thiết bị chuyên dụng: Hệ thống máy luân nhiệt PCR (Applied Biosystems, USA/Singapore), Bộ điện di phân giải cao (Scie-Plas Ltd, UK), Hệ thống chụp ảnh gel Gel Logic 1500 (Kodak, USA), Máy quang phổ hấp thụ UV-Vis Bio-mate 3 (Thermo Scientific), Tủ âm sâu $-80^\circ\text{C}$ (Ficchetti, Italy).
  • Hóa chất & Sinh phẩm: Taq DNA Polymerase 5U/$\mu$L, dNTPs 2.5 mM, DNA Ladder 1 kb (Thermo Scientific); CTAB, $\beta$-Mercaptoethanol, Isopropanol, Ethidium Bromide (Merck, Germany); Mồi tổng hợp độ tinh sạch cao (Macrogen Inc., Seoul, Hàn Quốc).
  • Phần mềm xử lý dữ liệu sinh học: NTSYSpc version 2.1 (Applied Biostatistics Inc.), Module SIMQUAL (tính chỉ số Jaccard/Dice), Module SAHN (phân cụm UPGMA).

Methodology

Nghiên cứu áp dụng quy trình thực nghiệm sinh học phân tử chuẩn hóa quốc tế với 4 cột mốc (milestones):

  1. Thu thập và định danh mẫu thực địa (Tháng 1 - 2): Thu thập 32 mẫu lá bánh tẻ sạch bệnh tại 4 vùng sinh thái chính (Bắc Giang: 4 mẫu; Thái Nguyên: 8 mẫu; Tuyên Quang: 12 mẫu; Viện Rau Quả: 8 mẫu).
  2. Chuẩn hóa quy trình tách chiết DNA (Tháng 3): Tối ưu hóa phương pháp Doyle & Doyle (1987) biến tính.
  3. Sàng lọc và chạy PCR phân tử (Tháng 4 - 5): Thử nghiệm 13 đoạn mồi phân giải cao với 40 chu kỳ nhiệt.
  4. Phân tích tin sinh học & Mô hình hóa (Tháng 6): Số hóa băng điện di, kiểm định bootstrap ma trận đa hình và xây dựng biểu đồ phát sinh loài.

Ma trận quản trị rủi ro thực nghiệm:

  • Rủi ro lẫn tạp chất Phenol/Polysaccharide: Bổ sung $\beta$-Mercaptoethanol 0.2% vào đệm CTAB và nâng thời gian lắng tủa ở $-20^\circ\text{C}$ lên tối thiểu 120 phút.
  • Rủi ro biến tính mồi ở nhiệt độ thấp: Thiết lập gradient nhiệt độ gắn mồi ($34^\circ\text{C} - 43.7^\circ\text{C}$) cho từng cặp mồi đặc thù.

Implementation và kết quả

Development process

1. Giao thức tách chiết và tinh sạch DNA tổng số

Áp dụng phương pháp CTAB (Cetyl Trimethyl Ammonium Bromide) cải tiến:

  • Nghiền 0.5g lá non trong $1000,\mu\text{L}$ đệm CTAB ($2%,\text{CTAB}, 100,\text{mM},\text{Tris-HCl pH 8.0}, 20,\text{mM},\text{EDTA}, 1.4,\text{M},\text{NaCl}$), bổ sung $\beta$-Mercaptoethanol.
  • Ủ ở $65^\circ\text{C}$ trong 30 phút. Bổ sung $500,\mu\text{L}$ CIAA (Chloroform : Isoamyl Alcohol = 24:1), ly tâm $10.000,\text{vòng/phút}$ trong 5 phút.
  • Thu pha trên, bổ sung $50,\mu\text{L},\text{CH}_3\text{COONa},3\text{M}$ và $500,\mu\text{L}$ Isopropanol lạnh, ủ $-20^\circ\text{C}$ trong 2 giờ.
  • Thu tủa, rửa sạch bằng $1000,\mu\text{L}$ Ethanol 70%, để khô tự nhiên 15 phút, hòa tan trong $50,\mu\text{L}$ đệm TE 1X.

2. Thiết kế và thông số động học của mồi phân tử

Bộ mồi 13 đoạn oligonucleotides được tổng hợp từ Macrogen (Hàn Quốc):

STT Ký hiệu mồi Trình tự Nucleotide ($5' \rightarrow 3'$) Nhiệt độ nóng chảy ($T_m$) Loại chỉ thị
1 ISSR - T1 (GT)6CC $43.7^\circ\text{C}$ ISSR
2 ISSR - T2 (CT)6TG $40.8^\circ\text{C}$ ISSR
3 ISSR - T3 (AC)6CG $43.7^\circ\text{C}$ ISSR
4 OPT - 01 GGGCCACTCA $34.0^\circ\text{C}$ RAPD
5 OPA - 04 AATCGGGCTG $32.0^\circ\text{C}$ RAPD
6 OPO - 04 AAGTCCGCTC $32.0^\circ\text{C}$ RAPD
7 OPA - 08 GTGACGTAGG $32.0^\circ\text{C}$ RAPD
8 OPC - 08 TGGACCGGTG $34.0^\circ\text{C}$ RAPD
9 OPM - 13 (GGT)2CAAG $32.0^\circ\text{C}$ RAPD
10 OPG - 16 AGCG(TCC)2 $34.0^\circ\text{C}$ RAPD
11 OPG - 17 ACGACC(GA)2 $32.0^\circ\text{C}$ RAPD
12 OPB - 18 C(CA)2GCAGT $32.0^\circ\text{C}$ RAPD
13 OPQ - 18 GGGCCACTCA $34.0^\circ\text{C}$ RAPD

3. Cấu hình phản ứng PCR và chu trình nhiệt

Thành phần cho 1 phản ứng ($50,\mu\text{L}$ tổng thể tích):

  • PCR Buffer 10X: $5.0,\mu\text{L}$
  • $\text{MgCl}_2$ (25 mM): $3.0,\mu\text{L}$
  • dNTPs (2.5 mM): $3.5,\mu\text{L}$
  • Primer ($10,\mu\text{M}$): $2.0,\mu\text{L}$
  • Taq DNA Polymerase ($5,\text{U}/\mu\text{L}$): $0.4,\mu\text{L}$
  • DNA khuôn mẫu ($20-50,\text{ng}/\mu\text{L}$): $5.0,\mu\text{L}$
  • Nước khử ion siêu sạch ($\text{ddH}_2\text{O}$): $31.1,\mu\text{L}$
# Chu trình luân nhiệt PCR (40 chu kỳ)
Giai đoạn 1 (Biến tính ban đầu):  94°C trong 5 phút
Giai đoạn 2 (Khuếch đại 40 chu kỳ):
Giai đoạn 3 (Kéo dài cuối cùng): 72°C trong 3 phút
Bảo quản sản phẩm:               4°C vô thời hạn

4. Thuật toán tính toán ma trận tương đồng di truyền (Python / SciPy Implementation)

Dữ liệu điện di được chuyển đổi thành ma trận nhị phân $M \times N$ (với $M=32$ mẫu giống và $N$ là tổng số phân đoạn DNA). Thuật toán hệ số tương đồng di truyền Jaccard ($S_{ij}$) và phân cụm UPGMA (Unweighted Pair Group Method with Arithmetic Mean) được mô hình hóa:

import numpy as np
from scipy.spatial.distance import pdist, squareform
from scipy.cluster.hierarchy import linkage, dendrogram
import matplotlib.pyplot as plt

def calculate_genetic_similarity(binary_matrix):
    """
    Tính ma trận tương đồng di truyền Jaccard:
    S_ij = a / (a + b + c)
    trong đó:
      a: số băng xuất hiện ở cả 2 mẫu (1, 1)
      b: số băng xuất hiện ở mẫu i nhưng không có ở j (1, 0)
      c: số băng xuất hiện ở mẫu j nhưng không có ở i (0, 1)
    """
    # Tính khoảng cách Jaccard (Distance = 1 - Similarity)
    jaccard_dist = pdist(binary_matrix, metric='jaccard')
    # Chuyển đổi thành ma trận tương đồng
    similarity_matrix = 1.0 - squareform(jaccard_dist)
    np.fill_diagonal(similarity_matrix, 1.0)
    
    # Thực hiện phân cụm phả hệ UPGMA
    linkage_matrix = linkage(jaccard_dist, method='average')
    return similarity_matrix, linkage_matrix

# Dữ liệu ma trận nhị phân thực nghiệm (32 mẫu x N locus phân đoạn)
# 1: Có băng DNA điện di; 0: Không xuất hiện băng

Testing và validation

Quá trình kiểm tra điện di và đo quang phổ xác nhận:

  • Toàn bộ 32 mẫu DNA tách chiết đều cho vạch sáng rõ, không lẫn RNA, không bị đứt đoạn gãy nát trên gel Agarose 1.0%. Tỷ số $OD_{260}/OD_{280}$ nằm trong khoảng chuẩn $1.82 - 1.95$.
  • Điện di sản phẩm PCR trên gel Agarose 2.0% ở hiệu điện thế $100,\text{V}$, cường độ $100,\text{mA}$ trong 35 phút cho độ phân giải băng sắc nét.

Thống kê số lượng băng khuếch đại đa hình qua các mồi đại diện:

  • Mồi ISSR-T1: Khuếch đại 4 phân đoạn kích thước $500 - 1500,\text{bp}$ (phân đoạn $500,\text{bp}$ hiện diện trên toàn bộ mẫu; phân đoạn $1500,\text{bp}$ xuất hiện đặc thù ở 3 mẫu).
  • Mồi ISSR-T2: Khuếch đại 6 phân đoạn ($250 - 1500,\text{bp}$), ghi nhận đa hình cao tại locus $700,\text{bp}$ (22 mẫu) và $750,\text{bp}$ (21 mẫu).
  • Mồi ISSR-T3: Tạo 12 phân đoạn kích thước trải rộng $300 - 2200,\text{bp}$.
  • Mồi RAPD OPB-18: Tạo 14 phân đoạn đa hình kích thước $250 - 3000,\text{bp}$.
  • Mồi RAPD OPG-16: Cho 13 phân đoạn ($250 - 3000,\text{bp}$), phân biệt rõ rệt nhóm mẫu giống đột biến và bản địa.

Kết quả đạt được

Ma trận hệ số tương đồng di truyền phân tích trên phần mềm NTSYSpc v2.1 dao động từ 0.64 đến 1.00, chứng minh tập đoàn giống nghiên cứu có tính đa hình di truyền rất cao.

                   CÂY PHÁT SINH CHỦNG LOÀI 32 GIỐNG CAM QUÝT (DENDROGRAM)
-----------------------------------------------------------------------------------------
Hệ số tương đồng:  0.64          0.74          0.84          0.94        1.00
-----------------------------------------------------------------------------------------

Chi tiết phân cụm phả hệ:

  1. Nhóm I (Ngoại nhóm): Duy nhất mẫu Chấp (CH), phân tách độc lập hoàn toàn với khoảng cách di truyền tương đối $d = 0.36$ (hệ số tương đồng 0.64 so với quần thể còn lại).
  2. Nhóm II (31 mẫu còn lại): Phân nhánh thành 2 nhóm phụ rõ rệt:
    • Nhóm phụ I (14 mẫu):
      • Nhóm quýt: Quýt sen (QS) và Quýt Ôn Châu (QO) đạt hệ số tương đồng $0.94$; liên kết với Quýt ngọt (QN) ở mức $0.80$.
      • Nhóm cam ngọt ngoại nhập / chọn dòng: Hai mẫu cam V2V2-1 đạt độ tương đồng tuyệt đối $1.00$ (xác nhận về mặt di truyền là cùng 1 dòng vô tính thuần); Cam chín sớm C36CS đạt $0.93$; Cam Xã Đoài Cao Phong (CP) và Xã Đoài Nghệ An (NA) đạt $0.88$.
      • Nhóm đột biến không hạt: Cam chín sớm ít hạt (BH) và Cam ruột đỏ Cara (Cr) có hệ số tương đồng $0.85$.
      • Tập đoàn cam sành Hàm Yên: 11 dòng cam sành Hàm Yên (S2-1 đến S19-3) hợp thành một nhánh đơn ngành vững chắc với độ đồng nhất nội bộ đạt $0.79 - 0.93$.
    • Nhóm phụ II (8 mẫu):
      • Mẫu Cam sành Bố Hạ số 1 (A) và số 2 (B) đạt độ tương đồng $0.90$.
      • Mẫu Cam Vinh (CV) và Cam chanh Thái Nguyên (CC) đạt độ tương đồng $0.90$.
      • Đặc biệt, Cam chanh Bố Hạ (CBH) có quan hệ di truyền rất gần gũi với cụm Cam Vinh/Cam Chanh với hệ số tương đồng lên tới 0.83 (83%), đồng thời liên kết với nhóm Cam sành Hàm Yên ở hệ số 0.75 (75%). Kết quả này chứng minh cam Vinh và cam Chanh trồng tại Thái Nguyên có nguồn gốc xuất xứ rất gần hoặc là con lai phát sinh từ dòng cam Bố Hạ cổ truyền.

Đổi mới và đóng góp

  1. Phát hiện quan hệ phát sinh loài mới: Cung cấp bằng chứng phân tử đầu tiên xác nhận Cam chanh Bố Hạ và Cam Vinh có hệ số tương đồng $0.83$, giải mã nguồn gốc lịch sử lai tạo và di thực của các giống cây có múi tại miền Bắc Việt Nam.
  2. Tối ưu hóa hệ thống chỉ thị phân tử kép: Kết hợp 10 mồi RAPD ngẫu nhiên và 3 mồi ISSR chuyên biệt locus vi vệ tinh, nâng độ phân giải và khả năng phát hiện đa hình lên $100%$ đối với các dòng cam sành thoái hóa.
  3. Hiệu quả định danh giống thuần: Nhận diện chính xác 100% dòng thuần vô tính (V2 trùng khớp hoàn toàn V2-1 ở hệ số 1.00), giúp loại bỏ hoàn toàn hiện tượng lẫn giống trong sản xuất cây giống thương phẩm.
  4. Cơ sở dữ liệu số phục vụ bảo tồn gen: Xây dựng bảng ma trận nhị phân 32 mẫu đạt chuẩn sinh tin học, làm tiền đề pháp lý cho việc đăng ký chỉ dẫn địa lý và bảo hộ giống cây trồng đặc sản tỉnh Bắc Giang.

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------+
|               CHIẾN LƯỢC TRIỂN KHAI ỨNG DỤNG THỰC TẾ (ROADMAP)              |
+-----------------------------------------------------------------------------+
   Giai đoạn 1 (Tháng 1 - 3): Sàng lọc cây đầu dòng S0 tại vùng Bố Hạ
   
   Giai đoạn 2 (Tháng 4 - 8): Vi nhân giống vô tính in-vitro sạch bệnh
   
   Giai đoạn 3 (Tháng 9 - 18): Xây dựng vườn giống gốc & chuyển giao kỹ thuật
  • Ứng dụng thực tế: Phục tráng giống cam đặc sản Bố Hạ tại huyện Yên Thế; phân loại nhanh các cá thể ưu tú trong vườn giống gốc của các trung tâm khuyến nông mà không cần chờ chu kỳ ra quả ($3-5$ năm).
  • Mô hình kinh tế (ROI Analysis): Việc khôi phục thành công 100 ha vùng chuyên canh cam Bố Hạ thuần chủng đạt năng suất 16 tấn/ha sau 5 năm sẽ mang lại giá trị sản lượng ước tính hàng chục tỷ đồng mỗi năm cho nông hộ địa phương, tái tạo thương hiệu nông sản đặc hữu.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Chỉ thị RAPD và ISSR bản chất là chỉ thị trội (dominant markers), chưa thể phân biệt rõ ràng trạng thái đồng hợp tử ($AA$) và dị hợp tử ($Aa$) tại các locus gen chức năng.
  • Ràng buộc tài nguyên: Chưa thực hiện phân tích sâu trên các chỉ thị đồng trội như SSR (Simple Sequence Repeat) hoặc SNP (Single Nucleotide Polymorphism) do chi phí mồi huỳnh quang và giải trình tự mao quản.
  • Hướng nghiên cứu tương lai:
    1. Phát triển chỉ thị SCAR (Sequence Characterized Amplified Region) đặc thù riêng cho giống cam Bố Hạ từ các băng đa hình $700,\text{bp}$ của mồi OPG-17 và ISSR-T2.
    2. Ứng dụng công nghệ giải trình tự thế hệ mới (NGS - RNA-Seq) để phân tích các gen quy định hàm lượng đường và mùi thơm đặc thù của cam Bố Hạ.

Đối tượng hưởng lợi

  • Sinh viên & Học viên Công nghệ Sinh học: Nắm vững quy trình chuẩn từ tách chiết DNA mô thực vật khó (nhiều polysaccharide), thiết kế luân nhiệt PCR, đến kỹ năng tin sinh học sử dụng NTSYSpc v2.1.
  • Nhà khoa học & Cán bộ nghiên cứu: Tiếp cận bộ dữ liệu hệ số tương đồng di truyền chi tiết của 32 mẫu cây có múi Việt Nam, làm cơ sở lai tạo giống mới kháng bệnh Greening và Tristeza.
  • Doanh nghiệp & Hợp tác xã Nông nghiệp: Có công cụ phân tử kiểm định chất lượng giống cây trồng đầu dòng, ngăn ngừa rủi ro mua nhầm giống giả/giống thoái hóa.
  • Cơ quan quản lý Nhà nước (Sở NN&PTNT Bắc Giang/Tuyên Quang): Cơ sở dữ liệu khoa học vững chắc để quy hoạch vùng trồng bảo tồn chỉ dẫn địa lý cam sành Bố Hạ và Hàm Yên.

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

1. Yêu cầu kỹ thuật tối thiểu của phòng lab để triển khai quy trình này là gì?

Cần trang bị hệ thống máy PCR luân nhiệt 96 giếng, bộ nguồn và bể điện di ngang gel Agarose, hệ thống đọc/chụp ảnh gel UV Transilluminator, máy ly tâm tốc độ cao ($14.000,\text{vòng/phút}$), máy đo quang phổ hấp thụ bước sóng $260/280,\text{nm}$ và tủ trữ đông sinh phẩm $-20^\circ\text{C}$ đến $-80^\circ\text{C}$.

2. Tại sao lại kết hợp đồng thời cả 2 loại chỉ thị RAPD và ISSR?

Chỉ thị RAPD sử dụng mồi ngẫu nhiên ngắn ($10,\text{bp}$) quét rộng trên toàn bộ hệ gen, trong khi chỉ thị ISSR nhắm vào các vùng vi vệ tinh lặp lại đơn giản ($SSR$) có tốc độ tiến hóa và đa hình rất cao. Phối hợp cả hai giúp giảm thiểu sai số âm tính giả, tăng hệ số phân giải thông tin đa hình ($PIC$) mà không làm tăng quá nhiều chi phí thí nghiệm.

3. Phương pháp này có phân biệt được dòng cam đột biến không hạt không?

Có. Trên cây dendrogram, mẫu cam chín sớm ít hạt (BH) và cam ruột đỏ (Cr) được tách riêng thành Cụm 2 với hệ số tương đồng $0.85$, chứng minh bộ mồi có khả năng phân tách các dòng biến dị hình thái đặc thù liên quan đến hạt và sắc tố.

4. Chi phí phân tích cho mỗi mẫu giống bằng kỹ thuật này là bao nhiêu?

Chi phí hóa chất cho một phản ứng PCR kết hợp phân tích điện di dao động từ $30.000 - 50.000,\text{VNĐ/phản ứng/mẫu}$. Đây là mức chi phí cực kỳ tiết kiệm (chỉ bằng $5-10%$ so với giải trình tự gen thế hệ mới), rất thích hợp để kiểm tra hàng loạt trong công tác chọn giống.

5. Kết quả phân tích có bị ảnh hưởng bởi mùa vụ thu mẫu lá không?

Không. Khác với chỉ thị hình thái hay protein/isozyme phụ thuộc vào giai đoạn sinh trưởng và điều kiện ngoại cảnh, chỉ thị phân tử DNA phản ánh trực tiếp cấu trúc nucleotide của bộ gen, do đó kết quả đồng nhất ở mọi giai đoạn phát triển của cây.


Kết luận

Đề tài "Xác định mối quan hệ di truyền của cam Bố Hạ với các giống cam khác bằng chỉ thị phân tử" đã thực hiện thành công việc tách chiết DNA tinh sạch từ 32 mẫu thực vật có múi và ứng dụng hiệu quả tổ hợp 10 mồi RAPD cùng 3 mồi ISSR. Dữ liệu tin sinh học từ phần mềm NTSYSpc v2.1 đã giải mã toàn diện bức tranh phả hệ với hệ số tương đồng dao động từ $0.64$ đến $1.00$, phân định 32 mẫu giống thành 2 nhóm lớn với khoảng cách di truyền rõ nét ($d = 0.36$). Nghiên cứu chứng minh mối quan hệ họ hàng di truyền mật thiết giữa cam chanh Bố Hạ với cam Vinh ($83%$) và cam sành Hàm Yên ($75%$). Đây là luận cứ khoa học then chốt, mở ra giải pháp công nghệ sinh học thực tiễn phục vụ công tác bảo tồn, tuyển chọn cây đầu dòng thuần chủng và phục tráng thành công thương hiệu cam đặc sản Bố Hạ tại Việt Nam.