Giới thiệu dự án

Cỏ cú (Cyperus rotundus L.) được Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO) và Hiệp hội Khoa học Cỏ dại Hoa Kỳ (WSSA) xếp vào nhóm "loài cỏ dại nguy hiểm hàng đầu thế giới" (world's worst weed). Theo các thống kê nông nghiệp quốc tế, loài dịch hại này đã xâm lấn trên 92 quốc gia thuộc vùng nhiệt đới và cận nhiệt đới, gây tổn thất nghiêm trọng cho hơn 50 loại cây trồng chủ lực. Tại Nam Á, cỏ cú làm sụt giảm 23% đến 89% năng suất cây trồng trong các vụ mùa chính; làm giảm 42% đến 50% sản lượng lúa nước và chiếm từ 60% đến 80% tổng sinh khối cỏ dại trong các vùng chuyên canh mía đường, gây thiệt hại kinh tế hàng tỷ USD mỗi năm.

Vấn đề cốt lõi (Problem Statement) nằm ở khả năng thích ứng sinh thái cực đoan và hệ thống sinh sản vô tính siêu việt của cỏ cú. Một củ cỏ cú đơn lẻ có thể tạo ra mạng lưới 550 – 1.200 chồi ngầm trong vòng 3 – 4 tháng, đạt mật độ 10 – 30 triệu củ/ha với ngân hàng hạt ngầm vượt mức 100 triệu hạt/ha/năm. Quá trình cơ giới hóa làm đất thường xuyên làm đứt đoạn thân rễ, vô tình kích thích các mắt ngủ trên củ nảy mầm theo cấp số nhân. Đồng thời, sự biến dị di truyền và biến động kiểu hình phát hoa sâu sắc khiến cỏ cú hình thành tính kháng các nhóm hoạt chất trừ cỏ phổ biến (như glyphosate, halosulfuron, lactofen, acetochlor), vô hiệu hóa các biện pháp kiểm soát hóa học truyền thống.

Mạng lưới nhân giống vô tính của Cỏ cú (Cyperus rotundus L.):
[Củ gốc / Initial Tuber] 

Mục tiêu cụ thể của dự án nghiên cứu:

  1. Định lượng sai biệt hình thái: Thu thập và đánh giá 24/25 chỉ tiêu hình thái phát hoa theo hệ thống mô tả chuẩn của Tayyar et al. (2003) trên các mẫu cỏ cú tại địa bàn huyện Củ Chi, TP. Hồ Chí Minh.
  2. Giải mã đa hình phân tử: Ứng dụng kỹ thuật RAPD (Random Amplified Polymorphic DNA) với bộ 13 chỉ thị mồi ngẫu nhiên (theo Molin et al., 2019) để xác lập bản đồ phân tử DNA bộ gen.
  3. Mô hình hóa mối quan hệ di truyền: Xây dựng ma trận khoảng cách Euclid (Euclidean distance), hệ số tương đồng Dice (Dice similarity coefficient), cây phân nhóm UPGMA (Unweighted Pair Group Method with Arithmetic Mean) và phân tích thành phần chính PCA (Principal Component Analysis).

Phương pháp tiếp cận tích hợp hình thái học và sinh học phân tử cho phép đánh giá chính xác tác động của địa hình canh tác (vùng gò đồi 10 – 15 m và vùng chuyển tiếp 5 – 10 m) cùng hệ thống cây trồng (cây ngắn ngày, cây ăn quả, rau ăn lá, rau ăn quả) đến sự phân hóa di truyền.

Phạm vi và giới hạn: Nghiên cứu thực hiện trên 10 mẫu đại diện tại huyện Củ Chi trong khung thời gian từ tháng 02/2022 đến tháng 07/2022, tập trung vào biến động hình thái hoa và đa hình DNA phân tử, không đi sâu vào phân tích lý hóa tính vi lượng của thổ nhưỡng.


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

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

Nghiên cứu cấu trúc quần thể cỏ dại đòi hỏi các công cụ chỉ thị có độ tin cậy cao nhằm phân biệt các dòng sinh thái (ecotypes) và cơ chế kháng thuốc.

Phương pháp chỉ thị Bản chất sinh học Ưu điểm kỹ thuật Hạn chế chính Nhu cầu DNA & Chi phí
Chỉ thị Hình thái (Morphological) Kiểu hình cơ quan sinh dưỡng & phát hoa Trực quan, chi phí hóa chất thấp, quan sát thực địa Bị ảnh hưởng mạnh bởi môi trường, độ biến thiên chủ quan Không cần DNA; Tốn diện tích & thời gian
Chỉ thị Tế bào học (Cytological) Số lượng & cấu trúc nhiễm sắc thể (NST) Phân biệt mức bội thể, chuyển đoạn NST Độ phân giải thấp, quy trình nhuộm NST phức tạp Mô phân sinh sống; Thiết bị quang học cao
Chỉ thị Phân tử RAPD (Molecular RAPD) Đa hình chiều dài đoạn DNA nhân bản ngẫu nhiên Đơn mồi, không cần biết trước trình tự genome, độ bao phủ rộng Chỉ thị tính trạng trội, yêu cầu chuẩn hóa nhiệt độ PCR 2 ng/µL DNA; Chi phí tối ưu, tốc độ nhanh
Chỉ thị SSR / Microsatellite Trình tự lặp lại đơn chuỗi nucleotide (2-6 bp) Chỉ thị đồng trội, độ đa hình cực cao, tái lập tuyệt đối Đòi hỏi thiết kế cặp mồi đặc hiệu loài, chi phí cao 10-50 ng DNA; Chi phí hóa chất đắt đỏ

Nghiên cứu giải quyết bài toán thiếu hụt dữ liệu di truyền học quần thể của Cyperus rotundus tại vùng đệm nông nghiệp TP.HCM, nơi áp lực thuốc diệt cỏ và sự luân canh cây trồng diễn ra với tần suất cao (4 – 12 lần xử lý/năm).

Thiết kế hệ thống

Hệ thống phân tích tích hợp quy trình thực nghiệm phòng Lab và chuỗi xử lý dữ liệu tính toán sinh học (bioinformatics computational pipeline):

Quy trình thực nghiệm & Phân tích tin sinh học:

Công cụ và thông số kỹ thuật thiết bị:

  • Phần mềm xử lý dữ liệu: NTSYSpc phiên bản 2.1 (Numerical Taxonomy and Multivariate Analysis System) kết hợp Microsoft Excel 2019.
  • Hệ thống nhân gen: Máy luân nhiệt 2720 Thermal Cycler (Applied Biosystems, Foster City, CA, USA).
  • Hệ thống điện di: Bể điện di nằm ngang Mupid-One (Nhật Bản) trên nền gel Agarose 2.0% w/v chuẩn độ phân giải cao trong dung dịch đệm 0.5X TAE.
  • Quang phổ kế kiểm chuẩn DNA: Máy quang phổ thể tích vi lượng BioDrop UV-Vis Spectrophotometer (BioDrop, UK).
  • Thiết bị đo trắc lượng hình thái: Thước cặp điện tử Mitutoyo 500-196-30 (Dải đo 0 – 150 mm, độ chính xác ±0.01 mm, Mitutoyo, Nhật Bản) và Kính hiển vi soi nổi kỹ thuật số KTST-978PRO (Độ phóng đại kiểm định 0 – 46 mm, sai số 0.01 mm, KTECH, Đài Loan).

Mô hình toán học ứng dụng trong phân tích di truyền:

  1. Khoảng cách di truyền Euclid (Phân tích hình thái): $$d(j, k) = \sqrt{\sum_{i=1}^{n} (x_{ij} - x_{ik})^2}$$ Trong đó: $x_{ij}$ và $x_{ik}$ lần lượt là giá trị chuẩn hóa của chỉ tiêu hình thái thứ $i$ trên mẫu $j$ và mẫu $k$.

  2. Hệ số tương đồng Dice (Phân tích băng DNA nhị phân): $$S_{Dice} = \frac{2a}{2a + b + c}$$ Trong đó: $a$ là số lượng băng DNA có mặt ở cả hai mẫu; $b$ là số băng chỉ có ở mẫu 1; $c$ là số băng chỉ có ở mẫu 2.

# Thuật toán chuyển đổi băng RAPD và tính ma trận Dice (Mô phỏng quy trình xử lý của NTSYSpc 2.1)
def calculate_dice_similarity(sample_a, sample_b):
    a = sum(1 for band_a, band_b in zip(sample_a, sample_b) if band_a == 1 and band_b == 1)
    b = sum(1 for band_a, band_b in zip(sample_a, sample_b) if band_a == 1 and band_b == 0)
    c = sum(1 for band_a, band_b in zip(sample_a, sample_b) if band_a == 0 and band_b == 1)
    if (2 * a + b + c) == 0:
        return 0.0
    return (2.0 * a) / (2.0 * a + b + c)

Methodology

Phương pháp luận nghiên cứu tuân thủ mô hình thực nghiệm 4 giai đoạn với cấu trúc kiểm soát sai số sinh học nghiêm ngặt:

  • Tiêu chuẩn thu thập mẫu: Thu thập ngẫu nhiên củ ngầm từ các cây đã xuất hiện phát hoa hoàn chỉnh trên 2 tiểu vùng địa hình: Gò đồi (độ cao 10 – 15 m: xã Phú Mỹ Hưng, An Nhơn Tây, Nhuận Đức) và Chuyển tiếp (độ cao 5 – 10 m: Tân Thạnh Đông, Trung An, Tân Phú Trung).
  • Chuẩn hóa giá thể: Toàn bộ 10 mẫu củ (kích thước chuẩn 8 – 12 mm × 6 – 9 mm) được trồng chung điều kiện tiểu khí hậu trong chậu 22 × 17 × 17 cm chứa giá thể đồng nhất (đất sạch : phân bò hoai : phân trùn quế = 2:1:1) trong 10 – 12 tuần nhằm triệt tiêu các biến động hình thái do sai khác dinh dưỡng đất cục bộ.
  • Đánh giá rủi ro và Kiểm soát chất lượng: Loại trừ nguy cơ ức chế phản ứng PCR do hợp chất thứ sinh (polyphenol/polysaccharide đặc thù của họ Cyperaceae) bằng cách bổ sung chất chống oxy hóa 1.0% Polyvinylpyrrolidone (PVP) và 0.2% 2-mercaptoethanol vào dung dịch đệm chiết SDS.

Implementation và kết quả

Development process

Quá trình triển khai kỹ thuật sinh học phân tử và phân tích số liệu diễn ra qua 4 bước thực nghiệm chi tiết:

Quy trình luân nhiệt PCR-RAPD (40 Chu kỳ):

Chi tiết giao thức ly trích DNA tổng số (Giao thức cải tiến SDS-PCI):

  • Bước 1: Nghiền mịn 0.2 g lá tươi trong 500 µL dung dịch chiết SDS (chứa 1% PVP và 0.2% 2-mercaptoethanol).
  • Bước 2: Tinh sạch protein bằng 500 µL PCI (Phenol : Chloroform : Isoamyl alcohol = 25:24:1 v/v), ly tâm 13.000 rpm trong 2 phút ở 4°C.
  • Bước 3: Chiết tách pha hữu cơ với 300 µL dung dịch CI (Chloroform : Isoamyl alcohol = 24:1 v/v), ly tâm 13.000 rpm trong 2 phút.
  • Bước 4 – 7: Kết tủa DNA bằng 120 µL Isopropanol 99% lạnh (-20°C trong 30 phút), ly tâm 12.000 rpm trong 10 phút, rửa cặn DNA 2 lần bằng Ethanol 70%.
  • Bước 8 – 10: Hòa tan cặn DNA trong 50 µL dung dịch TE 1X, đo mật độ quang học ở bước sóng $\lambda = 260\text{ nm}$ và $280\text{ nm}$.

Testing và validation

Nghiên cứu thực hiện giải trình và sàng lọc 13 mồi RAPD ngẫu nhiên 10-mer tổng hợp từ Integrated DNA Technologies (Hoa Kỳ):

Tên mồi Trình tự nucleotide (5' – 3') Hàm lượng GC (%) Khả năng nhân bản Trạng thái đa hình
OPB6 GTTTCGCTCC 60% Thành công Đa hình (Polymorphic)
OPB7 GGTGACGCAG 70% Thành công Đa hình (Polymorphic)
OPB14 TCCGCTCTGG 70% Thành công Đa hình (Polymorphic)
OPB17 AGGGAACGAG 60% Thành công Đa hình (Polymorphic)
OPB18 CCACAGCAGT 60% Thành công Đa hình (Polymorphic)
OPC2 GTGAGGCGTC 70% Thành công Đa hình (Polymorphic)
OPC8 TGGACCGGTG 70% Thành công Đa hình (Polymorphic)
OPG17 CACTCTCCTC 60% Thành công Đa hình (Polymorphic)
OPH2 AGAGGGCACA 60% Thành công Đa hình (Polymorphic)
OPH18 GAATCTCCTC 50% Thành công Đa hình (Polymorphic)
OPI13 CAGCGACAAG 60% Thành công Đa hình (Polymorphic)
OPL7 AGGCGGGAAG 70% Không tạo dải Không nhân bản
OPN6 GAGACGCACA 60% Không tạo dải Không nhân bản

Dải sản phẩm khuếch đại điện di nằm trong khoảng phân tử lượng 250 bp đến 2.500 bp đối chiếu theo thang chuẩn DNA 1 kb Ladder (Meridian Bioscience, Hoa Kỳ).

Thống kê biến động trắc lượng hình thái phát hoa (10 mẫu cỏ cú):

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

  1. Phân tích quan hệ hình thái hoa:

    • Khoảng cách di truyền Euclid biến thiên mạnh từ 1.01 đến 9.61.
    • Cặp mẫu tương đồng hình thái cao nhất là RAL01TRAL02T ($d = 1.01$, cùng thu thập tại vùng chuyên canh rau ăn lá huyện Củ Chi).
    • Cặp mẫu phân kỳ lớn nhất là CNN02TCAQ02D ($d = 9.61$).
    • Ở ngưỡng khoảng cách 5.22, 10 mẫu cỏ cú được phân chia rõ rệt thành 3 nhóm chính:
      • Nhóm I: Gồm 7 mẫu (CNN01D, CNN02D, CAQ01D, CNN01T, RAQ01T, RAL01T, RAL02T).
      • Nhóm II: Gồm 2 mẫu (CAQ02D, RAQ02T) mang đặc trưng phát hoa lớn và mật độ hoa cao (14 – 18 hoa/gié).
      • Nhóm III: Gồm duy nhất mẫu CNN02T mang các đột biến thu nhỏ kích thước hoa (chiều dài hoa thứ 3 chỉ đạt 0.35 mm).
  2. Phân tích đa hình di truyền RAPD:

    • 11/13 mồi biểu hiện tính đa hình sắc nét (tỷ lệ hữu dụng đạt 84.62%).
    • Hệ số tương đồng di truyền Dice dao động trong khoảng 0.15 đến 0.57.
    • Ở ngưỡng tương đồng di truyền 0.37, cây phát sinh chủng loại UPGMA phân lập quần thể cỏ cú thành 3 nhóm phân tử độc lập:
      • Nhóm RAPD 1: CNN01D, CNN02D, CAQ01D, CNN01T, RAQ01T.
      • Nhóm RAPD 2: CAQ02D, CNN02T, RAQ02T.
      • Nhóm RAPD 3: RAL01T, RAL02T.
  3. Tương quan giữa Di truyền và Địa hình Canh tác: Phân tích PCA chỉ ra rằng các yếu tố địa hình (gò đồi vs bưng trũng) và hệ thống luân canh không tạo ra sự cách ly di truyền hoàn toàn. Điều này minh chứng cho tính bảo tồn cao của phương thức nhân giống vô tính qua củ ngầm, đồng thời sự phát tán thụ phấn chéo bằng hạt dù tỷ lệ nảy mầm thấp (1 – 5%) vẫn đóng vai trò duy trì dòng gen (gene flow) đồng nhất qua các vùng canh tác.


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

Nghiên cứu mang lại những đóng góp khoa học và giá trị ứng dụng đột phá so với các công trình trước đây:

So sánh đối chiếu giá trị kỹ thuật với các công trình công bố:

Các đổi mới nổi bật:

  • Chuẩn hóa quy trình tách chiết DNA mô lá cỏ cú nhiều tạp chất: Cải tiến thành công dung dịch đệm SDS bổ sung 1% PVP và 0.2% 2-mercaptoethanol, loại bỏ hoàn toàn hiện tượng hóa nâu dịch chiết do polyphenol, đảm bảo độ tinh sạch $A_{260}/A_{280}$ đạt chuẩn quang phổ 1.8 – 2.0.
  • Tiết kiệm chi phí phân tích phân tử (Giảm >65%): Thiết lập thành công hệ thống phản ứng PCR-RAPD tối ưu thể tích 25 µL với nồng độ mồi cực thấp nhưng vẫn tái tạo 100% các dải băng đa hình từ 250 – 2500 bp mà không cần giải trình tự đắt tiền.
  • Cung cấp cơ sở dữ liệu vi biến dị đầu tiên tại TP.HCM: Xây dựng bộ ma trận kép (Hình thái học & Phân tử học) làm tiền đề phục vụ công tác dự báo tính kháng thuốc bảo vệ thực vật.

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

  1. Quản lý cỏ dại tích hợp (IWM - Integrated Weed Management):

    • Dữ liệu về sự đa dạng hình thái và mức độ biến thiên vòi nhụy/bao phấn (CV 17% – 34%) giải thích nguyên nhân các biện pháp kiểm soát cơ giới đơn lẻ thường thất bại.
    • Ứng dụng đề xuất: Kết hợp cày sâu 35 – 40 cm (làm khô kiệt 70 – 100% sức sống củ) luân phiên với phủ màng Polyethylene và sử dụng các chủng nấm đối kháng tiềm năng như Cercospora caricis hoặc Puccinia romagnoliana.
  2. Chiến lược sử dụng thuốc bảo vệ thực vật:

    • Dựa trên ma trận phân nhóm di truyền, các nhóm cỏ cú thuộc Nhóm RAPD 3 (RAL01T, RAL02T) có tần suất phun thuốc cao (6 – 8 lần/năm) cần được áp dụng luân phiên hoạt chất ức chế ALS (Halosulfuron) với các hoạt chất tiền nảy mầm (Acetochlor, Lactofen) nhằm ngăn chặn triệt để nguy cơ đột biến đơn dòng kháng chéo thuốc trừ cỏ.
  3. Lộ trình nhân rộng và chuẩn hóa:

    • Giai đoạn 1 (0 – 6 tháng): Tích hợp bộ 11 mồi RAPD đặc hiệu vào quy trình giám định nhanh kiểu sinh thái (ecotypes) cỏ cú tại các Chi cục Trồng trọt và Bảo vệ Thực vật.
    • Giai đoạn 2 (6 – 18 tháng): Mở rộng quét bản đồ di truyền trên toàn vùng kinh tế trọng điểm phía Nam (Đồng bằng sông Cửu Long và Đông Nam Bộ).

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

Mặc dù đạt được những kết quả định lượng chuẩn xác, đề tài vẫn ghi nhận một số giới hạn kỹ thuật:

  • Hạn chế chỉ thị trội: Kỹ thuật RAPD chỉ phân biệt được sự hiện diện hoặc vắng mặt của đoạn nhân bản (tính trạng trội), không thể phân biệt trạng thái kiểu gen đồng hợp tử (homozygous) và dị hợp tử (heterozygous).
  • Quy mô quần thể: Khảo sát tập trung trên 10 mẫu tiêu biểu đại diện cho 2 tiểu vùng địa hình huyện Củ Chi, chưa bao phủ hết các vùng thổ nhưỡng đất phèn, đất mặn ven sông Sài Gòn.

Hướng phát triển tiếp theo:

  1. Nâng cấp công nghệ chỉ thị: Ứng dụng công nghệ giải trình tự thế hệ mới GBS (Genotyping-by-Sequencing) hoặc chỉ thị vi vệ tinh SSR đa alen để xây dựng bản đồ liên kết gen chi tiết.
  2. Nghiên cứu biểu hiện gen kháng thuốc: Phân lập và giải trình tự gen mã hóa enzyme đích ALS (Acetolactate Synthase) và ACCase nhằm xác định chính xác các đột biến điểm (SNPs) quy định tính kháng hoạt chất trừ cỏ trên từng nhóm phân tử đã phân loại.

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

  • Sinh viên & Học viên chuyên ngành Nông học / Bảo vệ Thực vật: Tiếp cận tài liệu chuẩn hóa về phương pháp đánh giá 24 tính trạng phát hoa theo tiêu chuẩn quốc tế và cẩm nang quy trình PCR-RAPD chi tiết.
  • Kỹ sư Nông nghiệp & Nhà khoa học Cây trồng: Nắm bắt dữ liệu biến dị di truyền ($S_{Dice} = 0.15 - 0.57$) phục vụ công tác thiết kế thí nghiệm đồng ruộng và quản lý tính kháng cỏ dại.
  • Doanh nghiệp Sản xuất & Phân phối Thuốc BVTV: Tối ưu hóa hướng dẫn sử dụng sản phẩm, xây dựng công thức phối trộn hoạt chất tiền nảy mầm và hậu nảy mầm hiệu quả cho từng vùng sinh thái.
  • Cơ quan Quản lý Nông nghiệp Địa phương: Định hình chiến lược bảo vệ thực vật vùng đô thị và chuyển đổi cơ cấu cây trồng thích ứng với biến đổi môi trường canh tác.

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

1. Yêu cầu kỹ thuật tối thiểu trong phòng thí nghiệm để tái lập quy trình RAPD này là gì?

Phòng thí nghiệm cần trang bị máy nhân gen PCR (Thermal Cycler) kiểm soát nhiệt độ chính xác ±0.2°C, hệ thống điện di ngang gel agarose, buồng đọc UV/Blue LED transilluminator, máy quang phổ vi lượng đo nồng độ DNA (như BioDrop hoặc NanoDrop) và hóa chất Master Mix tiêu chuẩn.

2. Vì sao chỉ thị RAPD vẫn được lựa chọn thay vì giải trình tự toàn bộ hệ gen (WGS)?

RAPD sử dụng các mồi ngẫu nhiên 10-mer, không yêu cầu thông tin trình tự genome tham chiếu trước đó của loài cỏ cú (vốn có bộ gen phức tạp và tính đa bội thể cao), cho phép sàng lọc đa hình nhanh chóng với chi phí hóa chất chỉ bằng 10 – 15% so với giải trình tự thế hệ mới.

3. Sự phân bố địa hình gò đồi và chuyển tiếp tác động như thế nào đến di truyền cỏ cú?

Kết quả phân tích UPGMA và PCA chứng minh địa hình không tạo ra sự phân hóa di truyền triệt để giữa các nhóm. Cỏ cú duy trì tính đa hình cao trong nội bộ quần thể thông qua mạng lưới củ ngầm lan truyền cơ giới và khả năng thụ phấn chéo tự nhiên.

4. Tại sao cần đo đạc đến 24 chỉ tiêu phát hoa trong khi cỏ cú chủ yếu sinh sản vô tính?

Dù tỷ lệ hạt nảy mầm ngoài tự nhiên chỉ từ 1 – 5%, phát hoa là cơ quan phản ánh độ biến thiên kiểu hình trung thực nhất của thực vật có hoa dưới áp lực chọn lọc của môi trường canh tác và chế độ làm đất, giúp phát hiện các nòi sinh thái thích nghi mới.

5. Chi phí triển khai bộ kit giám định 11 mồi RAPD cho một đợt phân tích quần thể là bao nhiêu?

Chi phí hóa chất trung bình ước tính khoảng 45.000 – 60.000 VNĐ cho mỗi phản ứng PCR-RAPD đơn mồi (bao gồm tách chiết, hóa chất PCR, điện di và phân tích ảnh gel), mang lại hiệu quả kinh tế cao cho các cơ sở nghiên cứu và ứng dụng nông nghiệp.


Kết luận

Đồ án khóa luận tốt nghiệp "Đánh giá đa dạng di truyền mẫu cỏ cú (Cyperus rotundus L.) tại Thành phố Hồ Chí Minh" của tác giả Nguyễn Văn Phúc đã hoàn thành toàn diện các mục tiêu khoa học đề ra. Bằng cách tích hợp chuẩn mực trắc lượng 24 tính trạng hình thái phát hoa với kỹ thuật sinh học phân tử RAPD (11 mồi đa hình hiệu quả), công trình đã định lượng chính xác hệ số biến dị hình thái Euclid (1.01 – 9.61) và hệ số tương đồng di truyền Dice (0.15 – 0.57) trên 10 mẫu đại diện tại huyện Củ Chi.

Kết quả nghiên cứu khẳng định tính thích nghi sinh thái rộng và mức độ đa dạng nội tại phong phú của cỏ cú, đồng thời cung cấp luận cứ khoa học vững chắc phục vụ việc xây dựng chiến lược Quản lý cỏ dại tích hợp (IWM) bền vững, giảm thiểu ô nhiễm hóa chất bảo vệ thực vật trong nền nông nghiệp công nghệ cao.