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

Cây đước đôi (Rhizophora apiculata Blume) là loài thực vật chiếm ưu thế tại hệ sinh thái rừng ngập mặn Cần Giờ, giữ vai trò trọng yếu về mặt sinh thái học và kinh tế. Về mặt kinh tế, gỗ đước được sử dụng trong xây dựng, đóng đồ gia dụng, làm tà vẹt, chống lò, làm bột giấy và than đước cho nhiệt lượng cao, ít khói; vỏ cây giàu chất tanin phục vụ công nghiệp thuộc da, dược phẩm, nhuộm và keo dán. Về mặt môi trường, rừng đước đóng vai trò cân bằng sinh thái, điều hòa vi khí hậu, giảm xói lở bờ biển, chắn sóng và tạo bã mùn hữu cơ làm nền tảng cho chuỗi thức ăn cửa sông ven biển, đồng thời cung cấp sinh cảnh cho nhiều loài động vật quý hiếm.

Trong chiến tranh, rừng ngập mặn Cần Giờ bị tàn phá nặng nề bởi bom đạn và chất độc hóa học. Từ năm 1978, TP. Hồ Chí Minh đã chỉ đạo khôi phục hệ sinh thái này, phần lớn nguồn giống đước ban đầu được thu mua từ Năm Căn (Cà Mau). Ngày 21/01/2000, tổ chức UNESCO chính thức công nhận khu vực này là "Khu Dự trữ sinh quyển Rừng ngập mặn Cần Giờ". Tuy nhiên, quần thể đước trồng thuần loài đối mặt với nguy cơ mẫn cảm trước sâu bệnh và biến đổi môi trường. Do đó, việc nghiên cứu đa dạng di truyền bằng công cụ sinh học phân tử là yêu cầu cấp thiết nhằm phục vụ công tác quản lý, bảo tồn và định hướng chọn tạo giống phát triển rừng bền vững.

Khóa luận xác định các mục tiêu và nhiệm vụ cụ thể sau:

  1. Mục tiêu nghiên cứu:
    • Hoàn thiện kỹ thuật RAPD (Random Amplified Polymorphic DNA) trên đối tượng cây đước đôi (Rhizophora apiculata Blume).
    • Đánh giá mức độ đa dạng di truyền của quần thể đước đôi tại Khu dự trữ sinh quyển rừng ngập mặn Cần Giờ.
  2. Nhiệm vụ/Yêu cầu kỹ thuật:
    • Hoàn thiện thao tác và quy trình ly trích DNA tổng số từ mẫu lá đước đạt độ tinh sạch cao.
    • Khảo sát và tối ưu hóa các thông số phản ứng RAPD (nồng độ ion $\text{Mg}^{2+}$, nồng độ primer).
    • Ứng dụng phần mềm NTSYSpc để phân tích hệ số tương đồng di truyền và thiết lập sơ đồ phân nhóm (dendrogram).
  3. Đối tượng và phạm vi nghiên cứu:
    • Đối tượng: Quần thể đước đôi (Rhizophora apiculata Blume) tự nhiên và trồng qua các năm tại Cần Giờ.
    • Phạm vi không gian: Thu thập mẫu tại các tiểu khu thuộc Khu dự trữ sinh quyển rừng ngập mặn Cần Giờ (tổng diện tích khu dự trữ 75.740 ha).
    • Phạm vi thời gian: Đề tài được triển khai thực nghiệm từ ngày 08/05/2007 đến ngày 30/08/2007.

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

Khóa luận dựa trên nền tảng lý thuyết về hệ thống chỉ thị phân tử DNA trong nghiên cứu đa dạng di truyền thực vật, đặc biệt là kỹ thuật RAPD-PCR. RAPD sử dụng các đoạn mồi ngẫu nhiên ngắn (thường là 10 nucleotide) để khuếch đại các đoạn DNA genome mà không cần biết trước trình tự nucleotide. Mức độ đa hình di truyền được ghi nhận dựa trên sự xuất hiện (mã hóa 1) hoặc vắng mặt (mã hóa 0) của các phân đoạn DNA sau điện di. Mối quan hệ di truyền giữa các mẫu được tính toán dựa trên hệ số đồng dạng Nei và Li (1979): $$S_{xy} = \frac{2N_{xy}}{N_x + N_y}$$ Trong đó $N_{xy}$ là số băng chung giữa hai mẫu $x$ và $y$; $N_x$ và $N_y$ lần lượt là tổng số băng của mẫu $x$ và mẫu $y$. Khoảng cách di truyền được xác định qua công thức $D_{xy} = 1 - S_{xy}$. Dữ liệu ma trận nhị phân được xử lý trên phần mềm chuyên dụng NTSYSpc Version 2.1 (Numerical Taxonomy System).

Phương pháp nghiên cứu thực nghiệm được thiết kế qua các nội dung chính:

  • Phương pháp thu thập và xử lý mẫu: Thu thập 40 mẫu lá đước theo tuyến đường chéo (khoảng cách 1500 m thu 1 mẫu) kết hợp định vị tọa độ bằng hệ thống GPS. Thu thập cả lá non (có 2 lá phụ bao bọc) và lá già trên các cây có kiểu hình khác nhau (cây sinh trưởng mạnh, cằn cỗi, cây bị mối mọt, cây có u, quả có màu sắc khác biệt) và các năm tuổi khác nhau.
  • Phương pháp ly trích và kiểm tra DNA: Thử nghiệm so sánh 3 quy trình ly trích: Quy trình 1 (mẫu tươi theo Doyle và Doyle, 1988), Quy trình 2 (cải tiến bởi Lâm Vỹ Nguyên, 2006 có dùng $\text{N}2$ lỏng), và Quy trình 3 (quy trình hoàn thiện bổ sung Polyvidon và điều chỉnh thông số ly tâm). Kiểm tra định tính độ nguyên vẹn của DNA bằng điện di trên gel agarose 1% nhuộm Ethidium Bromide ($0,5\text{ }\mu\text{g/ml}$), quan sát dưới ánh sáng tử ngoại (UV). Đo mật độ quang học (OD) tại bước sóng 260 nm và 280 nm bằng máy quang phổ kế để xác định độ tinh sạch ($\text{OD}{260}/\text{OD}{280}$) và nồng độ DNA theo công thức: $$\text{Nồng độ DNA } (\mu\text{g/ml}) = [(62,9 \times \text{OD}{260}) - (36,0 \times \text{OD}_{280})] \times \text{Độ pha loãng}$$
  • Phương pháp tối ưu hóa RAPD: Thí nghiệm khảo sát nồng độ $\text{Mg}^{2+}$ (3 nghiệm thức: 2,0 mM; 2,5 mM; 3,0 mM), nồng độ mồi OPAC10 (4 nghiệm thức: 0,6; 0,8; 1,0; 1,2 pmol/$\mu\text{l}$), và thử nghiệm so sánh 2 mồi ngẫu nhiên gồm OPAC10 ($5'\text{-AGCAGCGAGG-3'}$, $T_m = 34^\circ\text{C}$) và OPA10 ($5'\text{-AGTGAACGCC-3'}$, $T_m = 30,7^\circ\text{C}$).
  • Nguồn dữ liệu: Dữ liệu thực nghiệm thu thập trực tiếp tại thực địa rừng ngập mặn Cần Giờ và các phân tích phòng thí nghiệm tại Viện Nghiên cứu Công nghệ Sinh học và Công nghệ Môi trường kết hợp Bộ môn Công nghệ Sinh học – Trường Đại học Nông Lâm TP.HCM.

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

Chương 1: Mở đầu

Chương 1 trình bày bối cảnh lịch sử khôi phục rừng Cần Giờ từ sau năm 1978 và việc UNESCO công nhận danh hiệu Khu dự trữ sinh quyển thế giới năm 2000. Tác giả nêu bật sự cần thiết của việc đánh giá đa dạng di truyền của cây đước đôi nhằm khắc phục tính mẫn cảm của rừng trồng thuần loài, hỗ trợ định hướng bảo tồn và tạo giống. Chương này đồng thời xác định rõ ràng mục đích, 3 yêu cầu cụ thể và giới hạn nghiên cứu về không gian, đối tượng, số lượng primer khảo sát và khung thời gian thực hiện.

Chương 2: Tổng quan tài liệu

Chương này tổng hợp các tài liệu khoa học liên quan đến địa bàn và đối tượng nghiên cứu:

  • Khu dự trữ sinh quyển Cần Giờ: Phân tích vị trí địa lý ($10^\circ22'14'' - 10^\circ37'9''\text{B}$, $106^\circ46'12'' - 107^\circ00'59''\text{Đ}$), điều kiện khí hậu, cấu trúc 3 phân vùng (Vùng lõi: 4.721 ha; Vùng đệm: 37.339 ha; Vùng chuyển tiếp: 29.310 ha; tổng diện tích 75.740 ha) và chức năng môi trường (hấp thụ 10.440 tấn $\text{CO}_2$, sinh ra 6.960 tấn $\text{O}_2$, lọc 587 triệu $\text{m}^3$ nước thải/năm).
  • Đặc điểm sinh học loài đước đôi: Phân loại học, đặc điểm hình thái thân, rễ chân nôm, cấu tạo hoa, quả có hiện tượng sinh con (vivipary) và phân bố địa lý.
  • Cơ sở sinh học phân tử: Cơ chế tách chiết DNA thực vật, các yếu tố gây nhiễu (polysaccharide, polyphenol); nguyên lý kỹ thuật PCR và chỉ thị phân tử (RFLP, AFLP, RAPD, SSR, STS); tổng hợp các công trình nghiên cứu sử dụng chỉ thị phân tử trên cây đước trong nước và quốc tế.

Chương 3: Vật liệu và phương pháp tiến hành

Chương 3 mô tả chi tiết quy trình thực nghiệm:

  • Danh mục trang thiết bị, máy luân nhiệt PCR, hệ thống điện di, máy quang phổ và hóa chất sử dụng (CTAB, Polyvidon, Chloroform:Isoamyl alcohol 24:1, Isopropanol, Taq DNA polymerase, mồi OPA10, OPAC10).
  • Trình bày chi tiết các bước tiến hành của 3 quy trình ly trích DNA so sánh.
  • Thiết kế 5 thí nghiệm cụ thể: Thí nghiệm 1 (khảo sát quy trình RAPD cơ bản), Thí nghiệm 2 (khảo sát nồng độ $\text{Mg}^{2+}$ ở 2,0 mM, 2,5 mM, 3,0 mM), Thí nghiệm 3 (khảo sát nồng độ mồi ở 0,6; 0,8; 1,0; 1,2 pmol/$\mu\text{l}$), Thí nghiệm 4 (thử nghiệm mồi OPA10 và OPAC10), Thí nghiệm 5 (thực hiện RAPD trên 10 mẫu đước đại diện theo quy trình tối ưu).
  • Chu kỳ luân nhiệt phản ứng RAPD: Biến tính ban đầu $94^\circ\text{C}$ (5 phút); 45 chu kỳ gồm [$94^\circ\text{C}$ trong 30 giây, $36^\circ\text{C}$ trong 30 giây, $72^\circ\text{C}$ trong 1 phút]; kéo dài kết thúc $72^\circ\text{C}$ (5 phút) và giữ ở $4^\circ\text{C}$.

Chương 4: Kết quả và thảo luận

Chương 4 phân tích các số liệu thực nghiệm thu được:

  • Kết quả thu mẫu và bảo quản: Thu thập 40 mẫu lá đước. Phân tích nguyên nhân lá già bị hóa nâu nhanh chóng khi rã đông từ $-20^\circ\text{C}$ (do tinh thể nước phá vỡ vách tế bào, giải phóng polyphenol), trong khi lá non có 2 lá phụ bảo vệ duy trì chất lượng tốt trên 4 tháng ở $-20^\circ\text{C}$.
  • Kết quả hoàn thiện ly trích DNA: Quy trình 1 cho DNA bị gãy vụn; Quy trình 2 cải thiện độ nguyên vẹn nhưng dịch chiết còn nhớt do dính polysaccharide; Quy trình 3 (bổ sung Polyvidon, nghiền nitơ lỏng) cho kết quả tối ưu. Đo quang phổ 40 mẫu cho thấy 65% số mẫu đạt tỷ lệ $\text{OD}{260}/\text{OD}{280} > 1,8$ và hàm lượng DNA trung bình đạt 243 ng/$\mu\text{l}$ (tất cả các mẫu đều $> 80\text{ ng/}\mu\text{l}$).
  • Kết quả tối ưu hóa RAPD: Xác định nồng độ $\text{Mg}^{2+}$ tối ưu là 3,0 mM và nồng độ mồi OPAC10 tối ưu là 0,6 pmol/$\mu\text{l}$.
  • Phân tích đa dạng di truyền: Phản ứng RAPD với mồi OPAC10 trên 10 mẫu đại diện ghi nhận 6 phân đoạn DNA (kích thước từ 100 đến 600 bp), gồm 5 băng đa hình (83,3%) và 1 băng đồng hình (16,7% ở kích thước 170 bp). Ma trận hệ số đồng dạng cho thấy mức tương đồng giữa các mẫu trồng cùng năm đạt 0,83 – 1,00; giữa các mẫu trồng năm 1978 và 1980 đạt 0,50 – 0,83; giữa nguồn giống tại chỗ (1991, 1996) là 0,50; giữa nguồn giống tại chỗ và nguồn giống gốc Cà Mau dao động từ 0,33 đến 0,67.

Chương 5: Kết luận và đề nghị

Tổng kết toàn bộ các kết quả đạt được về quy trình thu mẫu, bảo quản lá non, quy trình ly trích DNA chuẩn hóa bằng Polyvidon, hệ hóa chất RAPD tối ưu và bức tranh phân hóa di truyền của quần thể đước Cần Giờ. Đưa ra các khuyến nghị nghiên cứu tiếp theo về mở rộng số lượng primer và cỡ mẫu phân tích.

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

Nghiên cứu đã hoàn thiện các giải pháp kỹ thuật phòng thí nghiệm và cung cấp bộ dữ liệu thực nghiệm về cấu trúc di truyền của cây đước đôi tại Cần Giờ:

Thông số / Thành phần Kết quả thực nghiệm xác định trong khóa luận
Bảo quản mẫu lá Lá non có 2 lá phụ bao bọc bảo quản ở $-20^\circ\text{C}$ trên 4 tháng không bị hóa nâu
Quy trình ly trích DNA (Quy trình 3) Bổ sung Polyvidon vào dung dịch đệm EB, nghiền trong $\text{N}2$ lỏng; 65% mẫu đạt $\text{OD}{260}/\text{OD}_{280} > 1,8$; nồng độ DNA trung bình $243\text{ ng/}\mu\text{l}$
Thành phần phản ứng RAPD tối ưu (25 $\mu\text{l}$) PCR buffer 1X, $\text{MgCl}_2$ 3,0 mM, dNTP 0,2 $\mu\text{M}$, Primer 0,6 pmol/$\mu\text{l}$, Taq polymerase 1 U, DNA khuôn 50 ng
Mức độ đa hình DNA (Primer OPAC10) 6 phân đoạn DNA (100 – 600 bp); 5 băng đa hình (83,3%), 1 băng đồng hình (16,7% tại 170 bp)
Nhóm mẫu đước so sánh Hệ số tương đồng di truyền ($S_{xy}$)
Mẫu đước trồng cùng năm 0,83 – 1,00
Mẫu đước trồng năm 1978 so với năm 1980 0,50 – 0,83
Mẫu đước nguồn gốc tại chỗ (1991 so với 1996) 0,50
Mẫu đước nguồn tại chỗ so với nguồn giống Cà Mau 0,33 – 0,67

Các đóng góp chính do tác giả tự xác định:

  • Khắc phục thành công hiện tượng dịch trích bị nhớt do hàm lượng polysaccharide và hợp chất polyphenol cao trong lá đước bằng cách tối ưu hóa Quy trình 3 với Polyvidon.
  • Xây dựng thành công hệ hóa chất và chu kỳ nhiệt RAPD ổn định cho cây đước đôi, giải quyết vấn đề khó khuếch đại DNA ở loài cây rừng ngập mặn này.
  • Cung cấp dẫn liệu khoa học đầu tiên bằng chỉ thị RAPD chỉ ra sự phân hóa di truyền rõ rệt giữa nguồn đước du nhập từ Cà Mau và nguồn đước tái sinh tại chỗ ở Cần Giờ ($S_{xy}$ thấp, từ 0,33 đến 0,67).

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

Khóa luận được hoàn thành trong phạm vi giới hạn nghiên cứu mà tác giả đã chỉ rõ:

  • Hạn chế:
    • Quy mô mẫu phân tích RAPD chuyên sâu mới thực hiện trên 10 mẫu đại diện do giới hạn thời gian (từ tháng 5 đến tháng 8/2007) và kinh phí.
    • Số lượng primer sử dụng để khảo sát đa hình di truyền còn ít (chủ yếu tập trung phân tích sâu trên mồi OPAC10 sau khi thử nghiệm với OPA10).
  • Hướng nghiên cứu tiếp:
    • Mở rộng số lượng mồi RAPD ngẫu nhiên khác để thu nhận thêm nhiều locus đa hình, nâng cao độ chính xác của sơ đồ cây phát sinh di truyền.
    • Tăng cường số lượng mẫu thu thập trên toàn bộ các tiểu khu của cả 3 phân vùng (lõi, đệm, chuyển tiếp).
    • Kết hợp chỉ thị RAPD với các chỉ thị phân tử đồng trội có độ lặp lại cao hơn như SSR (Microsatellite) hoặc AFLP để đánh giá chi tiết tính trạng dị hợp tử và đa dạng di truyền quần thể.

Giá trị tham khảo

Công trình là tài liệu tham khảo hữu ích cho sinh viên, học viên cao học và cán bộ nghiên cứu thuộc các chuyên ngành Công nghệ Sinh học, Lâm nghiệp, Sinh học Thực vật và Quản lý Tài nguyên Môi trường:

  • Phần quy trình tách chiết DNA (Chương 3 và 4): Cung cấp giải pháp kỹ thuật chi tiết để xử lý các mô thực vật chứa hàm lượng cao chất chuyển hóa thứ cấp, polyphenol và polysaccharide.
  • Phần tối ưu hóa phản ứng RAPD: Cung cấp thông số cụ thể về nồng độ ion magie, nồng độ mồi và chu kỳ nhiệt mẫu cho các nghiên cứu đa hình di truyền thực vật rừng ngập mặn.
  • Dữ liệu phân loại học và đa dạng di truyền: Cung cấp cơ sở lý luận và số liệu thực tế phục vụ công tác quy hoạch, quản lý rừng phòng hộ và bảo tồn nguồn gen tại Ban quản lý Rừng phòng hộ Cần Giờ.

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

1. Tại sao lá đước già không thích hợp cho việc bảo quản và ly trích DNA trong nghiên cứu này?
Khi bảo quản ở nhiệt độ $-20^\circ\text{C}$, các tinh thể nước lớn hình thành phá vỡ vách tế bào lá đước già. Khi đưa ra nhiệt độ phòng, tinh thể tan ra khiến các hợp chất polyphenol bên trong phóng thích ra ngoài làm lá hóa nâu và dập nát chỉ trong khoảng 10 giây, gây phân hủy hoặc ức chế quá trình tách chiết DNA. Ngược lại, lá non được bảo vệ bởi 2 lá phụ màu đỏ thẫm giúp giữ trạng thái tươi ổn định trên 4 tháng ở $-20^\circ\text{C}$ và không bị hóa nâu trong 15 phút sau khi rã đông.

2. Hóa chất nào được tác giả bổ sung vào quy trình ly trích DNA để loại bỏ hiện tượng dịch trích bị nhớt?
Tác giả đã bổ sung Polyvidon vào dung dịch đệm ly trích EB (Extraction Buffer) trong Quy trình 3. Chất này có tác dụng liên kết và loại bỏ các hợp chất polysaccharide và polyphenol hòa tan trong mô lá đước, giúp dịch chiết không bị nhớt và tạo điều kiện tủa DNA tinh sạch.

3. Tỷ lệ đa hình di truyền ghi nhận được qua phản ứng RAPD với primer OPAC10 là bao nhiêu?
Phản ứng RAPD với primer OPAC10 ($5'\text{-AGCAGCGAGG-3'}$) trên 10 mẫu đước tạo ra tổng cộng 6 phân đoạn DNA (kích thước 100 – 600 bp), trong đó có 5 băng mang tính đa hình (chiếm tỷ lệ 83,3%) và 1 băng đồng hình (chiếm tỷ lệ 16,7% ở kích thước 170 bp).

4. Mối quan hệ di truyền giữa cây đước có nguồn gốc tại chỗ và cây đước có nguồn gốc từ Cà Mau thể hiện như thế nào?
Kết quả phân tích bằng phần mềm NTSYSpc cho thấy các cây đước nguồn gốc tại chỗ và cây đước có nguồn gốc di thực từ Năm Căn (Cà Mau) có khoảng cách di truyền tương đối xa, với hệ số đồng dạng di truyền ($S_{xy}$) dao động ở mức khá thấp, từ 0,33 đến 0,67.

Kết luận

Khóa luận của tác giả Võ Công Danh đã chuẩn hóa thành công quy trình ly trích DNA tổng số từ lá non cây đước đôi bằng phương pháp sử dụng Polyvidon và nghiền nitơ lỏng, cho hiệu suất và độ tinh sạch cao. Nghiên cứu đã xác lập hệ hóa chất tối ưu cho phản ứng RAPD với primer OPAC10 ($3,0\text{ mM }\text{MgCl}_2$, $0,6\text{ pmol/}\mu\text{l}$ primer) và ghi nhận tỷ lệ đa hình di truyền đạt 83,3%. Phân tích quan hệ di truyền qua phần mềm NTSYSpc đã làm rõ sự biến thiên hệ số tương đồng di truyền giữa các năm trồng và phản ánh sự phân hóa rõ rệt giữa nguồn giống tại chỗ Cần Giờ và nguồn giống di thực từ Cà Mau.