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

Ngành chăn nuôi gia cầm giữ vai trò kinh tế quan trọng trong cơ cấu nông nghiệp Việt Nam, cung cấp nguồn thực phẩm dồi dào về thịt và trứng với chu kỳ sinh trưởng ngắn, khả năng quay vòng vốn nhanh. Năm 2021, tổng đàn gà cả nước đạt trên 420,2 triệu con (Cục Chăn nuôi, 2021), trong đó các giống gà bản địa và gà đặc sản ngày càng chiếm thị phần lớn. Gà Hắc Phong là giống gà thuốc có nguồn gốc từ Trung Quốc, hiện được nuôi phổ biến tại Quảng Ninh và một số địa phương. Giống gà này sở hữu những đặc điểm ngoại hình và phẩm chất thịt đặc trưng: toàn bộ lông, da, thịt, xương và chân đều có màu đen, hàm lượng mỡ thấp, thịt thơm ngon và có giá trị dinh dưỡng, hỗ trợ bồi bổ sức khỏe (Phạm Công Thiếu và ctv, 2018; Nguyễn Thị Phương Giang và ctv, 2022).

Tuy nhiên, do phương thức chăn nuôi hiện nay chủ yếu mang tính nông hộ và quy mô trang trại nhỏ lẻ, công tác chọn lọc giống chưa được đầu tư bài bản, dẫn đến nguy cơ thoái hóa giống và suy giảm năng suất sinh sản theo thời gian. Sự phát triển của công nghệ sinh học di truyền phân tử mở ra giải pháp chọn lọc giống thông qua chỉ thị phân tử (Marker-Assisted Selection - MAS), giúp rút ngắn thời gian và nâng cao độ chính xác so với phương pháp chọn giống truyền thống dựa trên kiểu hình. Trong số các gen ứng viên liên quan đến tính trạng sinh sản ở gia cầm, gen Neuropeptide Y (NPY) giữ vai trò điều hòa thần kinh - nội tiết quan trọng, kiểm soát sự tiết hormone Gonadotropin từ tuyến yên, tác động trực tiếp đến sự phát triển của buồng trứng và tuổi thành thục sinh dục (Contijoch và ctv, 1993; Dunn và ctv, 2004).

Mục tiêu nghiên cứu: Phân tích đa hình gen Neuropeptide Y bằng kỹ thuật PCR-RFLP và đánh giá ảnh hưởng của đa hình gen này đến một số chỉ tiêu sinh sản ở giống gà Hắc Phong.

Nhiệm vụ nghiên cứu:

  1. Khuếch đại đoạn gen mục tiêu NPY bằng phản ứng chuỗi polymerase (PCR).
  2. Nhận diện và phân tích đa hình gen bằng enzyme cắt giới hạn DraI.
  3. Phân tích ảnh hưởng của đa hình gen NPY đến một số chỉ tiêu sinh sản (khối lượng cơ thể lúc 19 tuần tuổi và năng suất trứng tích lũy đến 38 tuần tuổi) của gà Hắc Phong.

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

  • Đối tượng nghiên cứu: 398 cá thể gà Hắc Phong (gồm 49 gà trống và 349 gà mái) được nuôi giữ tại Trung tâm Nghiên cứu và Phát triển Chăn nuôi Gia cầm VIGOVA (Phân viện Chăn nuôi Nam Bộ - Viện Chăn nuôi). Toàn bộ đàn gà được chăm sóc trong cùng điều kiện chuồng trại, thức ăn và nước uống.
  • Phạm vi không gian: Thu mẫu tại Trung tâm VIGOVA; phân tích phân tử tại Phòng thí nghiệm Công nghệ Phôi Động vật, Viện Nghiên cứu Công nghệ Sinh học và Môi trường, Khoa Khoa học Sinh học, Trường Đại học Nông Lâm TP. Hồ Chí Minh.
  • Phạm vi thời gian: Đề tài được triển khai thực hiện từ tháng 7/2022 đến tháng 11/2022.

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

Cơ sở lý thuyết

Gen Neuropeptide Y (NPY) ở gà định vị trên nhiễm sắc thể số 2, có kích thước khoảng 7 kb, cấu trúc gồm 3 intron và 4 exon mã hóa chuỗi peptide gồm 97 axit amin (NCBI). NPY là một peptide thần kinh phân bố rộng rãi trong hệ thống thần kinh trung ương, đóng vai trò cầu nối điều hòa giữa trạng thái dinh dưỡng và hoạt động sinh sản (Gladysz và Polkowska, 2006). NPY tham gia kiểm soát cảm giác thèm ăn, điều hòa giải phóng Gonadotropin-Releasing Hormone (GnRH) tại vùng dưới đồi, qua đó tác động đến sự rụng trứng và tuổi đẻ quả trứng đầu tiên (Hilal và ctv, 1996; Dunn và ctv, 2004; Xu và ctv, 2011a).

Kỹ thuật PCR-RFLP (Polymerase Chain Reaction - Restriction Fragment Length Polymorphism) kết hợp phản ứng khuếch đại gen đặc hiệu với sự phân cắt của enzyme giới hạn. Khi có sự xuất hiện của đột biến mất đoạn (indel) hoặc đột biến điểm (SNP), vị trí nhận biết đặc hiệu của enzyme bị thay đổi, tạo ra các phân đoạn DNA có kích thước khác nhau sau khi phân cắt, từ đó cho phép nhận diện các kiểu gen đồng hợp và dị hợp trong quần thể nghiên cứu.

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

  1. Thu thập mẫu và tách chiết DNA tổng số:

    • Thu 398 mẫu máu từ tĩnh mạch cánh của gà Hắc Phong vào ống nghiệm chống đông chứa EDTA, bảo quản ở -30°C.
    • Tách chiết DNA tổng số bằng bộ kit TopPURE® blood DNA extraction (ABT - Việt Nam) theo nguyên lý bắt giữ DNA trên màng cột silica qua 7 bước: phá vỡ màng tế bào bằng BL Buffer và Proteinase K ở 72°C trong 10 phút, kết tủa bằng ethanol, rửa sạch tạp chất qua cột silica bằng dung dịch WBI và WB2 (ly tâm 13.000 vòng/phút), giải hấp phụ DNA bằng đệm EB Buffer.
    • Đo nồng độ và độ tinh sạch DNA ở bước sóng 260/280 nm bằng máy NanoBioDrop (BioDrop - Anh); kiểm tra độ toàn vẹn của DNA bằng điện di trên gel agarose 1%.

  2. Khuếch đại đoạn gen mục tiêu bằng phản ứng PCR:

    • Sử dụng cặp mồi xuôi (5'-TCTCAGAGCTCCAACGTATGA-3') và mồi ngược (5'-ATATTTCTGTGCCTGAACAACA-3') tương ứng với vị trí nucleotide 415-666 của gen NPY (Dunn và ctv, 2004; Xu và ctv, 2011a).
    • Thành phần phản ứng PCR (tổng thể tích 12,5 µL): 6,25 µL MyTaq™ Mix 2X (Bioline - Anh), 0,4 µL mồi xuôi (10 pmol/µL), 0,4 µL mồi ngược (10 pmol/µL), 1,0 µL DNA khuôn (10 ng/µL) và 4,45 µL nước khử nuclease.
    • Chu trình nhiệt trên máy luân nhiệt Thermal Cycler Bioer: Tiền biến tính ở 95°C (3 phút); 35 chu kỳ gồm biến tính ở 95°C (30 giây), bắt cặp ở 56°C (30 giây), kéo dài ở 72°C (30 giây); kéo dài lần cuối ở 72°C (7 phút).
    • Kiểm tra sản phẩm PCR trên gel agarose 1,5% nhuộm GelRed 0,6X, điện di ở hiệu điện thế 100 V trong 30 phút. Gửi sản phẩm giải trình tự Sanger để kiểm tra độ đặc hiệu.

  3. Phân tích đa hình gen bằng enzyme cắt giới hạn DraI:

    • Phản ứng phân cắt (tổng thể tích 11 µL) gồm: 2 µL sản phẩm PCR, 1 µL enzyme DraI (BioLab - Anh), 3 µL Buffer và 5 µL nước khử nuclease. Hỗn hợp được ủ ở 37°C trong 16 giờ, sau đó bất hoạt enzyme ở 65°C trong 15 phút.
    • Điện di sản phẩm cắt trên gel agarose 2,5% ở 100 V trong 40 phút.
    • Quy ước kiểu gen tại locus NPY/DraI:
      • Kiểu gen AA (không bị cắt): phân đoạn 252 bp.
      • Kiểu gen AB (dị hợp tử): 3 phân đoạn 252 bp, 167 bp và 81 bp.
      • Kiểu gen BB (đồng hợp tử đột biến): 2 phân đoạn 167 bp và 81 bp.
    • Giải trình tự các mẫu đại diện cho 3 kiểu gen để xác định vị trí đột biến mất 4 nucleotide tạo vị trí nhận biết của enzyme DraI (TTTAAA).

  4. Thu thập dữ liệu kiểu hình và xử lý thống kê:

    • Dữ liệu khối lượng gà ở 19 tuần tuổi (trước khi vào chu kỳ đẻ) và năng suất trứng tích lũy đến tuần tuổi thứ 38 được theo dõi trên 349 cá thể gà mái.
    • Tần số alen và tần số kiểu gen tính theo định luật Hardy-Weinberg.
    • Tần số dị hợp mong đợi ($H_e$) tính theo công thức Nei (1978): $$H_e = 1 - \sum_{i=1}^{n} P_i^2$$
    • Hệ số thông tin đa hình (PIC) tính theo công thức Botstein và ctv (1980): $$PIC = 1 - \sum_{i=1}^{n} P_i^2 - \sum_{i=1}^{n-1}\sum_{j=i+1}^{n} 2P_i^2 P_j^2$$
    • Phân tích phương sai một yếu tố (ANOVA) và kiểm định sai khác giữa các giá trị trung bình bằng phép thử Tukey trên phần mềm thống kê (mức ý nghĩa $P < 0,05$). Số liệu được trình bày dưới dạng Mean ± SEM.

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

Chương 1. Đặt vấn đề

Chương này nêu bật vị trí của ngành chăn nuôi gia cầm tại Việt Nam và tiềm năng kinh tế - dinh dưỡng của giống gà thuốc Hắc Phong. Tác giả chỉ ra thực trạng thiếu hụt dữ liệu chọn giống di truyền phân tử trên đối tượng này, từ đó xác lập tính cấp thiết, mục tiêu và 3 nội dung nghiên cứu trọng tâm xoay quanh việc ứng dụng kỹ thuật PCR-RFLP khảo sát gen NPY.

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

Chương 2 hệ thống hóa cơ sở lý luận về đặc điểm sinh học của giống gà Hắc Phong, cấu trúc và vai trò sinh lý của gen Neuropeptide Y, cũng như nguyên lý kỹ thuật PCR, chỉ thị phân tử RFLP và phương pháp PCR-RFLP. Tác giả tổng hợp các công trình nghiên cứu trong nước (Trần Thị Bình Nguyên và ctv, 2018 trên gà Liên Minh; Nguyễn Trọng Ngữ và ctv, 2015 trên gà Nòi) và nghiên cứu quốc tế (Dunn và ctv, 2004; Xu và ctv, 2011a; Li và ctv, 2009; Fatemi và ctv, 2012) chứng minh mối tương quan giữa các đột biến tại gen NPY với tuổi đẻ quả trứng đầu tiên và năng suất trứng ở các giống gia cầm.

Chương 3. Vật liệu và phương pháp

Chương này mô tả chi tiết địa điểm nghiên cứu (Trường ĐH Nông Lâm TP.HCM và Trung tâm VIGOVA), đối tượng 398 cá thể gà Hắc Phong, toàn bộ quy trình thực nghiệm gồm 7 bước tách chiết DNA bằng kit TopPURE®, thông số cặp mồi, hóa chất, chu trình nhiệt PCR 35 chu kỳ, điều kiện phân cắt của enzyme DraI (37°C trong 16 giờ), các nồng độ gel agarose điện di (1%, 1,5% và 2,5%), cùng hệ thống công thức toán học xác định tần số di truyền ($H_e$, $H_o$, PIC) và mô hình phân tích phương sai ANOVA/Tukey.

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

Chương này trình bày các phát hiện thực nghiệm qua 4 nội dung chính:

  • Đánh giá chất lượng DNA tổng số: Kết quả điện di trên gel agarose 1% cho thấy các giếng chỉ xuất hiện một vệt DNA duy nhất, sáng rõ, không đứt gãy, đạt tiêu chuẩn thực hiện PCR.
  • Khuếch đại đoạn gen mục tiêu: Sản phẩm PCR có kích thước chính xác 252 bp, vệt sáng đồng đều, không tạo dimer mồi, giếng đối chứng âm không xuất hiện vệt. Kết quả giải trình tự mẫu số 790 khẳng định đoạn gen khuếch đại tương đồng hoàn toàn với trình tự nucleotide từ vị trí 415-666 của gen NPY.
  • Nhận diện đa hình gen NPY/DraI: Sau khi phân cắt bằng enzyme DraI, ghi nhận sự phân cắt đa hình với 2 alen (A, B) và 3 kiểu gen (AA, AB, BB). Giải trình tự xác nhận đột biến mất 4 nucleotide (indel) đã hình thành trình tự nhận diện TTTAAA của DraI. Phân tích 398 cá thể cho thấy alen B (0,751) và kiểu gen BB (0,555) chiếm ưu thế. Hệ số PIC đạt 0,304 (đa hình mức độ trung bình), $H_e = 0,374$ và $H_o = 0,392$, phản ánh quần thể có xu hướng giao phối ngẫu nhiên chiếm ưu thế.
  • Ảnh hưởng của đa hình gen đến chỉ tiêu sinh sản: Khối lượng gà mái 19 tuần tuổi giữa 3 kiểu gen không có sự sai khác mang ý nghĩa thống kê ($P > 0,05$). Tuy nhiên, sản lượng trứng tích lũy đến 38 tuần tuổi ở nhóm mang kiểu gen BB đạt 55,9 quả, cao hơn có ý nghĩa thống kê ($P < 0,05$) so với nhóm AB (50,0 quả) và nhóm AA (48,5 quả).

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

Chương 5 tóm tắt các kết quả cốt lõi: khuếch đại thành công đoạn gen 252 bp, xác định 3 kiểu gen tại locus NPY/DraI, chứng minh kiểu gen BB có ưu thế vượt trội về sản lượng trứng ở 38 tuần tuổi. Đồng thời, tác giả đưa ra đề nghị mở rộng nghiên cứu trên nhiều thế hệ tiếp theo để củng cố độ tin cậy của chỉ thị phân tử phục vụ công tác chọn lọc giống.

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

Số liệu và kết quả thực nghiệm chính

Bảng tổng hợp các thông số di truyền quần thể tại locus NPY/DraI trên 398 cá thể gà Hắc Phong:

Đối tượng Số cá thể ($n$) Tần số alen A Tần số alen B Tần số kiểu gen AA Tần số kiểu gen AB Tần số kiểu gen BB Dị hợp mong đợi ($H_e$) Hệ số đa hình (PIC)
Gà trống 49 0,306 0,694 0,082 (4 con) 0,449 (22 con) 0,469 (23 con) 0,425 0,335
Gà mái 349 0,238 0,762 0,049 (17 con) 0,378 (132 con) 0,573 (200 con) 0,363 0,297
Chung 398 0,249 0,751 0,053 (21 con) 0,392 (156 con) 0,555 (221 con) 0,374 0,304

Bảng ảnh hưởng của các kiểu gen NPY/DraI đến khối lượng và năng suất trứng của gà Hắc Phong:

Kiểu gen Số cá thể mái ($n$) Khối lượng ở 19 tuần tuổi (g/con) Năng suất trứng tích lũy đến 38 tuần tuổi (quả)
AA 17 $1.252,7 \pm 10,1^{\text{a}}$ $48,5 \pm 1,5^{\text{a}}$
AB 132 $1.248,3 \pm 5,2^{\text{a}}$ $50,0 \pm 0,6^{\text{a}}$
BB 200 $1.250,5 \pm 4,6^{\text{a}}$ $55,9 \pm 0,5^{\text{b}}$

(Ghi chú: Trong cùng một cột, các giá trị trung bình mang chữ cái khác nhau biểu thị sự sai khác có ý nghĩa thống kê ở mức $P < 0,05$).

Đóng góp học thuật và thực tiễn

  • Đóng góp dữ liệu phân tử: Nghiên cứu là công trình thực nghiệm đầu tiên công bố dữ liệu về tần số alen, tần số kiểu gen và các chỉ số đa hình di truyền của locus NPY/DraI trên quần thể gà Hắc Phong tại Việt Nam.
  • Xác định chỉ thị chọn giống tiềm năng: Phát hiện đột biến mất 4 nucleotide (tạo alen B và kiểu gen BB) có tác động tích cực làm tăng sản lượng trứng (đạt 55,9 quả so với 48,5 - 50,0 quả ở nhóm mang alen A), trong khi không làm thay đổi tiêu cực đến khối lượng cơ thể gà mái ở giai đoạn chuẩn bị vào đẻ (19 tuần tuổi).
  • So sánh với các nghiên cứu trước: Tác giả chỉ ra sự khác biệt về phân bố kiểu gen giữa các giống gà: ở gà Hắc Phong, kiểu gen BB chiếm tỷ lệ cao nhất (55,54%), trái ngược với giống gà Liên Minh ở Việt Nam vốn có kiểu gen không đột biến AA chiếm đa số (76%) (Trần Thị Bình Nguyên và ctv, 2018), hoặc gà Văn Xương ở Trung Quốc có kiểu gen dị hợp AB chiếm ưu thế (50%) (Li và ctv, 2009).

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

  • Hạn chế được ghi nhận từ văn bản:

    • Dữ liệu năng suất trứng mới chỉ được thu thập và theo dõi tích lũy đến 38 tuần tuổi (giai đoạn đầu của chu kỳ sinh sản), chưa bao quát toàn bộ chu kỳ đẻ kéo dài đến 72 tuần tuổi.
    • Nghiên cứu được thực hiện trên một thế hệ cá thể tại một cơ sở chăn nuôi tập trung, chưa khảo sát qua nhiều thế hệ liên tiếp để đánh giá mức độ ổn định di truyền của tính trạng.
    • Đề tài chỉ tập trung phân tích một locus đơn lẻ (NPY/DraI), chưa kết hợp đánh giá đồng thời với các gen ứng viên sinh sản khác (như Prolactin, GnRHR, IGF-I).

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

    • Mở rộng quy mô theo dõi đa hình gen locus NPY/DraI qua nhiều thế hệ tiếp theo nhằm củng cố độ chính xác của cơ sở dữ liệu di truyền.
    • Theo dõi sản lượng trứng trong toàn bộ chu kỳ sinh sản hoàn chỉnh để hỗ trợ xây dựng quy trình chọn giống gà Hắc Phong theo hướng nâng cao năng suất trứng bằng chỉ thị phân tử.

Giá trị tham khảo

  • Đối tượng thụ hưởng:
    • Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Công nghệ Sinh học, Di truyền học động vật, Chăn nuôi và Thú y.
    • Các kỹ thuật viên phòng thí nghiệm sinh học phân tử thực hiện các đề tài ứng dụng chỉ thị PCR-RFLP trong chọn giống vật nuôi bản địa.
  • Nội dung giá trị chuyên sâu:
    • Quy trình chi tiết và điều kiện tối ưu hóa phản ứng PCR (nhiệt độ bắt cặp 56°C, thời gian 35 chu kỳ) và phản ứng cắt enzyme giới hạn DraI (ủ 37°C trong 16 giờ) đối với đoạn gen NPY dài 252 bp.
    • Phương pháp luận phân tích thống kê di truyền quần thể (Hardy-Weinberg, $H_e$, PIC) kết hợp kiểm định mối tương quan giữa kiểu gen và tính trạng số lượng bằng mô hình ANOVA/Tukey.

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

1. Kích thước đoạn gen mục tiêu được khuếch đại bằng PCR là bao nhiêu và tương ứng với vị trí nào trên gen Neuropeptide Y?
Đoạn gen mục tiêu được khuếch đại có kích thước 252 bp (đối với alen không đột biến) hoặc 248 bp (đối với alen xảy ra đột biến mất 4 nucleotide). Trình tự này nằm từ vị trí nucleotide số 415 đến 666 của gen NPY.

2. Enzyme cắt giới hạn DraI phân cắt sản phẩm PCR thành các đoạn có kích thước như thế nào đối với từng kiểu gen?

  • Kiểu gen AA (không đột biến, không có vị trí cắt): giữ nguyên phân đoạn 252 bp.
  • Kiểu gen AB (dị hợp tử, một alen bị cắt): xuất hiện 3 phân đoạn có kích thước 252 bp, 167 bp và 81 bp.
  • Kiểu gen BB (đồng hợp tử đột biến, cả 2 alen bị cắt): tạo thành 2 phân đoạn có kích thước 167 bp và 81 bp.

3. Tần số alen và kiểu gen ghi nhận trên tổng đàn 398 cá thể gà Hắc Phong là bao nhiêu?
Tần số alen A là 0,249 và alen B là 0,751. Tần số 3 kiểu gen lần lượt là: AA chiếm 0,053 (21 cá thể); AB chiếm 0,392 (156 cá thể); BB chiếm 0,555 (221 cá thể). Tần số alen B và kiểu gen BB chiếm ưu thế trong quần thể.

4. Đa hình gen tại locus NPY/DraI ảnh hưởng như thế nào đến khối lượng và năng suất trứng của gà Hắc Phong?
Đa hình gen không làm thay đổi có ý nghĩa thống kê về khối lượng gà mái ở 19 tuần tuổi ($1.252,7$ g ở nhóm AA; $1.248,3$ g ở nhóm AB; $1.250,5$ g ở nhóm BB; $P > 0,05$). Tuy nhiên, đối với năng suất trứng tích lũy đến 38 tuần tuổi, nhóm gà mang kiểu gen BB đạt 55,9 quả, cao hơn có ý nghĩa thống kê ($P < 0,05$) so với nhóm AA (48,5 quả) và AB (50,0 quả).

5. Hệ số thông tin đa hình (PIC) và tần số dị hợp mong đợi (He) của quần thể đạt giá trị bao nhiêu?
Hệ số thông tin đa hình (PIC) của quần thể chung đạt 0,304 (thuộc mức đa hình trung bình). Tần số dị hợp mong đợi ($H_e$) là 0,374, thấp hơn tần số dị hợp quan sát được ($H_o = 0,392$), phản ánh quần thể có xu hướng giao phối ngẫu nhiên chiếm ưu thế.

Kết luận

Khóa luận đã ứng dụng thành công kỹ thuật PCR-RFLP với enzyme cắt giới hạn DraI để khuếch đại đoạn gen NPY kích thước 252 bp và xác định 3 kiểu gen (AA, AB, BB) trên quần thể 398 cá thể gà Hắc Phong. Kết quả phân tích chỉ ra kiểu gen đồng hợp đột biến BB chiếm tỷ lệ cao nhất (55,54%) và có ưu thế vượt trội về sản lượng trứng tích lũy đến 38 tuần tuổi (55,9 quả) so với kiểu gen AA và AB ($P < 0,05$), trong khi không làm ảnh hưởng đến khối lượng cơ thể ở 19 tuần tuổi. Đây là nguồn dữ liệu khoa học thực nghiệm có giá trị, định hướng cho việc sử dụng locus NPY/DraI làm chỉ thị phân tử hỗ trợ chọn lọc nâng cao năng suất trứng ở giống gà Hắc Phong.