(Luận án tiến sĩ) chọn tạo hai dòng vịt biển trên cơ sở giống vịt biển 15 đại xuyên

HỌC VIỆN NƠNG NGHIỆP VIỆT NAM CHU HỒNG NGA CHỌN TẠO HAI DÒNG VỊT BIỂN TRÊN CƠ SỞ GIỐNG VỊT BIỂN 15 - ĐẠI XUYÊN LUẬN ÁN TIẾN SĨ NHÀ XUẤT BẢN HỌC VIỆN NÔNG NGHIỆP - 2021 HỌC VIỆN NƠNG NGHIỆP VIỆT NAM CHU HỒNG NGA CHỌN TẠO HAI DÒNG VỊT BIỂN TRÊN CƠ SỞ GIỐNG VỊT BIỂN 15 - ĐẠI XUYÊN Ngành: Chăn nuôi Mã số: 9.62.01.05 Người hướng dẫn khoa học: TS Nguyễn Thanh Sơn GS.TS Đặng Vũ Bình HÀ NỘI - 2021 LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng tôi, kết nghiên cứu trình bày luận án trung thực, khách quan chưa để bảo vệ học vị Tôi xin cam đoan giúp đỡ cho việc thực luận án cảm ơn, thơng tin trích dẫn luận án rõ nguồn gốc Hà Nội, ngày tháng năm 2021 Tác giả luận án Chu Hoàng Nga i LỜI CẢM ƠN Để hồn thành luận án này, tơi xin bày tỏ lịng kính trọng biết ơn sâu sắc đến thầy hướng dẫn khoa học: TS Nguyễn Thanh Sơn, GS.TS Đặng Vũ Bình tận tình hướng dẫn, dành nhiều công sức, thời gian để hướng dẫn tơi suốt q trình học tập, thực đề tài viết luận án Tôi xin bày tỏ lòng biết ơn chân thành Ban Giám đốc, Ban Quản lý đào tạo, Bộ môn Di truyền Giống vật nuôi, Khoa Chăn nuôi - Học viện Nông nghiệp Việt Nam tận tình giúp đỡ tơi q trình học tập, thực đề tài hoàn thành luận án Tôi xin chân thành cảm ơn tập thể lãnh đạo, cán bộ, công nhân viên Trung tâm Nghiên cứu Vịt Đại Xuyên thành viên tham gia Đề tài khoa học cấp Bộ: “Nghiên cứu chọn tạo dịng vịt Biển phục vụ chăn ni vùng xâm ngập mặn” tạo điều kiện để tơi hồn thành luận án Xin chân thành cảm ơn Ban Giám đốc, Khoa Quân nhu, Bộ môn Sản xuất, cán bộ, nhà giáo Học viện Hậu cần Ban Lãnh đạo, Phịng Quản lý học viên cán Đồn 871 - Tổng cục Chính trị - Bộ Quốc Phịng, tạo điều kiện, thời gian để tham gia học tập hoàn thành luận án Xin chân thành cảm ơn nhà khoa học, chuyên gia lĩnh vực chăn nuôi gia cầm, tạo điều kiện, ủng hộ giúp đỡ suốt q trình học tập, nghiên cứu để hồn thành luận án Xin chân thành cảm ơn gia đình, người thân, bạn bè cổ vũ, giúp đỡ mặt, động viên tạo điều kiện thuận lợi để tơi hồn thành luận án./ Hà Nội, ngày tháng năm 2021 Nghiên cứu sinh Chu Hoàng Nga ii MỤC LỤC Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục chữ viết tắt vi Danh mục bảng vii Danh mục hình ix Trích yếu luận án x Thesis abstract xii Phần Mở đầu 1.1 Tính cấp thiết đề tài 1.2 Mục tiêu đề tài 1.2.1 Mục tiêu chung 1.2.2 Mục tiêu cụ thể 1.3 Phạm vi nghiên cứu 1.4 Những đóng góp đề tài 1.5 Ý nghĩa khoa học thực tiễn đề tài 1.5.1 Ý nghĩa khoa học 1.5.2 Ý nghĩa thực tiễn Phần Tổng quan tài liệu 2.1 Cơ sở khoa học vấn đề nghiên cứu 2.1.1 Tính trạng số lượng 2.1.2 Các tham số di truyền 2.1.3 Giá trị giống 2.1.4 Hiệu chọn lọc 10 2.1.5 Các tính trạng vịt yếu tố ảnh hưởng 10 2.2 Tình hình nghiên cứu ngồi nước 22 2.2.1 Tình hình nghiên cứu ngồi nước 22 2.2.2 Tình hình nghiên cứu nước 27 2.3 Đánh giá nghiên cứu vịt biển 35 Phần Nội dung phƣơng pháp nghiên cứu 36 3.1 Nội dung nghiên cứu 36 iii 3.2 Vật liệu phương pháp nghiên cứu 36 3.2.1 Vật liệu nghiên cứu 36 3.2.2 Phương pháp nghiên cứu 37 3.2.3 Xử lý thống kê 47 Phần Kết thảo luận 48 4.1 Chọn tạo dòng trống HY1 48 4.1.1 Một số đặc điểm ngoại hình 48 4.1.2 Ảnh hưởng yếu tố cố định tham số di truyền khối lượng thể 49 4.1.3 Khối lượng thể HY1 qua hệ tuần tuổi 52 4.1.4 Khảo sát sinh trưởng hàm toán học 55 4.1.5 Kích thước chiều đo thể HY1 61 4.1.6 Khả cho thịt HY1 qua hệ 63 4.1.7 Tiêu tốn thức ăn/kg tăng khối lượng HY1 qua hệ 68 4.1.8 Tỷ lệ đẻ suất trứng 20 tuần đẻ vịt HY1 qua hệ .71 4.1.9 Chất lượng trứng vịt HY1 qua hệ 72 4.1.10 Một số tiêu ấp nở vịt HY1 qua hệ 75 4.2 Kết chọn tạo dòng mái HY2 77 4.2.1 Một số đặc điểm ngoại hình dịng HY2 77 4.2.2 Ảnh hưởng số yếu tố đến cố định tham số di truyền khối lượng thể vịt lúc tuần tuổi suất trứng 20 tuần đẻ 78 4.2.3 Tỷ lệ đẻ suất trứng 20 tuần đẻ vịt HY2 qua hệ chọn lọc 81 4.2.4 Tiêu tốn thức ăn 20 tuần đẻ HY2 qua hệ 85 4.2.5 Tỷ lệ đẻ suất trứng vịt HY2 52 tuần đẻ 87 4.2.6 Chất lượng trứng vịt HY2 qua hệ 95 4.2.7 Một số tiêu ấp nở trứng vịt HY2 99 4.2.8 Khối lượng thể vịt HY2 qua tuần tuổi 101 Phần Kết luận kiến nghị 104 5.1 Kết luận 104 5.2 Kiến nghị 105 Danh mục cơng trình cơng bố liên quan đến luận án 106 iv Tài liệu tham khảo 107 Phụ lục 119 v DANH MỤC CHỮ VIẾT TẮT Từ viết tắt Nghĩa tiếng Việt BLUP Dự đốn tuyến tính khơng thiên vị tốt ĐBSCL Đồng sông Cửu Long Cs Cộng DT Dài thân HY1 Dòng trống vịt Biển 15 - Đại Xuyên HY2 Dòng mái vịt Biển 15 - Đại Xuyên EBV Giá trị giống ước tính NT Ngày tuổi NST Năng suất trứng Hệ số xác định SD Độ lệch chuẩn TB Trung bình THXP Thế hệ xuất phát TL Tỷ lệ TTTĂ Tiêu tốn thức ăn VN Vòng ngực R vi DANH MỤC BẢNG TT Tên bảng Trang 3.1 Giá trị dinh dưỡng mức cho ăn vịt HY1 HY2 theo giai đoạn ni 38 4.1 Một số đặc điểm ngoại hình vịt HY1 48 4.2 Các yếu tố ảnh hưởng tới khối lượng thể HY1 49 4.3 Phương sai hiệp phương sai di truyền kiểu hình khối lượng thể HY1 hệ 50 4.4 Hệ số di truyền, tương quan di truyền tương quan kiểu hình khối lượng thể HY1 hệ 50 4.5 Phương sai hiệp phương sai di truyền kiểu hình khối lượng thể HY1 hệ 50 4.6 Hệ số di truyền, tương quan di truyền tương quan kiếu hình khối lượng thể HY1 hệ 51 4.7 Khối lượng vịt mái HY1 qua hệ 52 4.8 Khối lượng vịt trống HY1 qua hệ 54 4.9 Hàm Gompertz vịt mái trống HY1 hệ 55 4.10 Hàm Richards vịt mái trống HY1 hệ 55 4.11 Khối lượng thể tiệm cận, tuổi khối lượng thể HY1 điểm uốn hệ theo hàm Gompertz 58 4.12 Khối lượng thể tiệm cận, tuổi khối lượng thể HY1 điểm uốn hệ theo hàm Richards 59 4.13 Khối lượng chiều đo thể lúc tuần tuổi vịt mái HY1 qua hệ 61 4.14 Khối lượng chiều đo thể lúc tuần tuổi vịt trống HY1 qua hệ 62 4.15 Năng suất thịt xẻ lúc tuần tuổi vịt mái HY1 qua hệ 63 4.16 Năng suất thịt xẻ lúc tuần tuổi vịt trống HY1 qua hệ 64 4.17 Năng suất thịt xẻ lúc tuần tuổi vịt mái HY1 qua hệ 64 4.18 Năng suất thịt xẻ lúc tuần tuổi vịt trống HY1 qua hệ 65 4.19 Năng suất thịt xẻ lúc tuần tuổi vịt mái HY1 qua hệ 65 4.20 Năng suất thịt xẻ lúc tuần tuổi vịt trống HY1 qua hệ 66 4.21 Tiêu tốn thức ăn/kg tăng khối lượng HY1 qua hệ 68 vii 4.22 Tỷ lệ đẻ, suất trứng HY1 qua hệ 71 4.23 Chất lượng trứng HY1 qua hệ 73 4.24 Kết ấp nở HY1 qua hệ 75 4.25 Một số đặc điểm ngoại hình vịt HY2 77 4.26 Các yếu tố ảnh hưởng tới khối lượng thể tuần tuổi suất trứng/20 tuần đẻ HY2 78 4.27 Phương sai, hiệp phương sai di truyền kiểu hình khối lượng thể lúc tuần tuổi suất trứng/20 tuần đẻ HY2 hệ 79 4.28 Các tham số di truyền khối lượng thể lúc tuần tuổi suất trứng/20 tuần đẻ HY2 hệ 79 4.29 Phương sai, hiệp phương sai di truyền kiểu hình khối lượng thể lúc tuần tuổi suất trứng/20 tuần đẻ HY2 hệ 79 4.30 Các tham số di truyền khối lượng thể lúc tuần tuổi suất trứng/20 tuần đẻ HY2 hệ 79 4.31 Tỷ lệ đẻ suất trứng/20 tuần đẻ HY2 qua hệ 81 4.32 So sánh tỷ lệ đẻ, suất trứng/20 tuần đẻ HY2 qua hệ 82 4.33 Tiêu tốn thức ăn/10 trứng 20 tuần đẻ HY2 qua hệ 85 4.34 Tỷ lệ đẻ, suất trứng từ 21 đến 52 tuần đẻ HY2 qua hệ 88 4.35 So sánh tỷ lệ đẻ suất trứng 21 - 52 tuần đẻ HY2 qua hệ 90 4.36 So sánh tỷ lệ đẻ suất trứng/52 tuần đẻ HY2 qua hệ 93 4.37 Chất lượng trứng HY2 qua hệ 95 4.38 Kết ấp nở HY2 qua hệ 99 4.39 Khối lượng thể vịt mái HY2 qua hệ 101 4.40 Khối lượng thể vịt trống HY2 qua hệ 101 viii Function to be estimated: a*exp(-b*exp(-k*t)) Initial parameter estimates: a = 2800.0 b = 4.18 k = 0.35 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 18 Estimation Results Parameter a b k Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 98.0386 percent R-Squared (adjusted for d.f.) = 98.0371 percent Standard Error of Est = 141.972 Mean absolute error = 85.9154 Durbin-Watson statistic = 1.70109 Lag residual autocorrelation = 0.149352 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3724.96*exp(-4.33335*exp(-0.35643*t)) 128 Hàm Gompertz – Trống HY1-2 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*exp(-b*exp(-k*t)) Initial parameter estimates: a = 2800.0 b = 4.18 k = 0.35 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 18 Estimation Results Parameter a b k Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 98.0859 percent R-Squared (adjusted for d.f.) = 98.0829 percent Standard Error of Est = 132.136 Mean absolute error = 80.5345 Durbin-Watson statistic = 1.91204 Lag residual autocorrelation = 0.0439104 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3895.52*exp(-4.16733*exp(-0.327043*t)) 129 Hàm Richard, Mái HY1-0 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*(1-b*exp(-k*t))^(-1/n) Initial parameter estimates: a = 2500.0 b = 0.06 k = 0.36 n = -0.015 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 35 Estimation Results Parameter a b k n Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 97.5749 percent R-Squared (adjusted for d.f.) = 97.5714 percent Standard Error of Est = 117.73 Mean absolute error = 87.8815 Durbin-Watson statistic = 1.59412 Lag residual autocorrelation = 0.202017 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3714.03*(1-0.0473722*exp(-0.299021*t))^(-1/-0.0116363) 130 Hàm Richard, Trống HY1-0 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*(1-b*exp(-k*t))^(-1/n) Initial parameter estimates: a = 2800.0 b = 0.04 k = 0.34 n = -0.01 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 36 Estimation Results Parameter a b k n Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 98.3957 percent R-Squared (adjusted for d.f.) = 98.3888 percent Standard Error of Est = 104.832 Mean absolute error = 78.1499 Durbin-Watson statistic = 1.38229 Lag residual autocorrelation = 0.306917 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3895.52*(1-0.0336886*exp(-0.295065*t))^(-1/-0.00853496) 131 Hàm Richard, Mái HY1-1 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*(1-b*exp(-k*t))^(-1/n) Initial parameter estimates: a = 2500.0 b = 0.06 k = 0.36 n = -0.015 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 22 Estimation Results Parameter a b k n Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 97.5967 percent R-Squared (adjusted for d.f.) = 97.5941 percent Standard Error of Est = 150.77 Mean absolute error = 89.9118 Durbin-Watson statistic = 1.6654 Lag residual autocorrelation = 0.166816 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3756.57*(1-0.0566409*exp(-0.324053*t))^(-1/-0.0136542) 132 Hàm Richard, Trống HY1-1 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*(1-b*exp(-k*t))^(-1/n) Initial parameter estimates: a = 2800.0 b = 0.04 k = 0.34 n = -0.01 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 23 Estimation Results Parameter a b k n Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 97.7375 percent R-Squared (adjusted for d.f.) = 97.7334 percent Standard Error of Est = 145.448 Mean absolute error = 92.3114 Durbin-Watson statistic = 1.71327 Lag residual autocorrelation = 0.142065 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3898.13*(1-0.0298767*exp(-0.312091*t))^(-1/-0.00748724) 133 Hàm Richard, Mái HY1-2 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*(1-b*exp(-k*t))^(-1/n) Initial parameter estimates: a = 2500.0 b = 0.06 k = 0.36 n = -0.015 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 23 Estimation Results Parameter a b k n Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 98.0383 percent R-Squared (adjusted for d.f.) = 98.0362 percent Standard Error of Est = 142.006 Mean absolute error = 86.1081 Durbin-Watson statistic = 1.7009 Lag residual autocorrelation = 0.149449 Residual Analysis n MSE MAE MAPE ME MPE 134 The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3759.85*(1-0.0631561*exp(-0.34921*t))^(-1/-0.0149845) Hàm Richard, Trống HY1-2 Nonlinear Regression - KL Dependent variable: KL Independent variables: t Function to be estimated: a*(1-b*exp(-k*t))^(-1/n) Initial parameter estimates: a = 2800.0 b = 0.04 k = 0.34 n = -0.01 Estimation method: Marquardt Estimation stopped due to convergence of residual sum of squares Number of iterations: Number of function calls: 23 Estimation Results Parameter a b k n Analysis of Variance Source Model Residual Total Total (Corr.) R-Squared = 98.086 percent R-Squared (adjusted for d.f.) = 98.0816 percent Standard Error of Est = 132.183 Mean absolute error = 80.4556 Durbin-Watson statistic = 1.91219 Lag residual autocorrelation = 0.0438358 135 Residual Analysis n MSE MAE MAPE ME MPE The StatAdvisor The output shows the results of fitting a nonlinear regression model to describe the relationship between KL and independent variables The equation of the fitted model is KL = 3922.22*(1-0.039639*exp(-0.32259*t))^(-1/-0.00967747) Kết tính EBV dịng HY1 .PEST UIUC 4.2.3 Monday April 10:27:26 2019 PREDICT BV HY1 TH2 P7 page Dimension of equation system : 1731 Rank 1731 of equation system : General I n f o r m a t i o n: _ Usable memory : Memory used # of equations : : 32.000 MB # of nonzero elements : # of data records Parameter 21.98 % 1731 6063 : file 350000 925 : hy1-2-7 Results file : EBVHY1-2-7 Input data file(raw) : /data/dat(hy1-2-7).dat Input rel file(raw) : /data/ped(hy1-2-7).ped Input data file(bin) : MEMORY Input rel file(bin) : MEMORY _ R u n T i m e I n f o r m a t i o n: _ CPU time spent in hour:min:sec data preparation : setting up and solving Total wall clock time / memory : 00:00:00 00:00:00 : 00:00:08 memory(kb) 793 7196 7202 _ Data File I n f o r m a t i o n: COVs/TRAITs n mean std.dev CV max scaled by _ number of number of total number of animals _ .PEST UIUC 4.2.3 Monday April 10:27:26 2019 PREDICT BV HY1 TH2 P7 136 PEST UIUC 4.2.3 Monday April 10:27:26 2019 SOLVERS # Monday April 10:27:26 2019 SEX ANIMAL - EBV vịt HY1 hệ + Chọn mái theo EBV(P7): 20 vịt có giá trị cao 137 1217 1479 1477 1264 1118 1117 209 133 1079 1136 1307 887 1265 1530 579 886 ………………………………………… + Chọn trống theo EBV(P7): 20 vịt có giá trị cao SỐ CÁNH 384 332 708 605 674 145 141 321 544 59 197 917 389 425 682 397 58 140 146 136 - EBV HY1 hệ + Chọn mái theo giá trị EBV(P7): 20 vịt có giá trị cao SỐ CÁNH 709 383 213 388 603 139 711 713 138 143 1217 1479 1218 356 1477 387 710 712 1264 138 + Chọn trống theo giá trị EBV(P7): 20 vịt có giá trị cao SỐ CÁNH 384 332 708 605 674 145 141 321 544 59 197 917 389 425 682 397 58 140 146 136 Kết tính EBV suất trứng HY2 ‫ﯯ‬ .PEST UIUC 4.2.3 Friday October 08:25:32 2019 PREDICT BV NST HY2 TH2 TM Dimension Rank G e n e r a l I n f o r m a t i o n: _ Usable memory Memory used # of equations # of nonzero elements # of data records Parameter Results Input data Input rel Input data Input rel _ R u n T i m e I n f o r m a t i o n: CPU time spent in data preparation setting up and solving Total wall clock time _ D a t a F i l e I n f o r m a t i o n: COVs/TRAITs 139 NST R E L A T I O N S H I P I n f o r m a t i o n: _ number of number of total number of animals _ .PEST UIUC 4.2.3 Friday October 08:25:32 2019 PREDICT BV NST HY2 TH2 TM .PEST UIUC 4.2.3 Friday October 08:25:32 2019 PREDICT BV NST HY2 TH2 TM ANIMAL 1601 1602 1603 1604 1605 - EBV NST 20 vịt HY2 đạt cao Số cánh 698 697 691 849 692 700 682 699 681 696 693 690 683 687 863 679 850 817 762 685 908 140 MỘT SỐ HÌNH ẢNH Khay ấp trứng kẹp thẻ nhôm vịt lúc ngày tuổi Ghép trống vào gia đình Dịng HY2 Khay trứng chuồng đẻ cá thể 141 Bấm thẻ nhựa thu trứng Mổ Khảo sát Đo kích thước trứng cân khối lượng lịng đỏ 142 ... phương pháp chọn lọc quần thể tách nhóm vịt HY1 để tạo dòng trống HY2 để tạo dòng mái Để nâng cao suất giống vịt Biển 15 – Đại Xuyên, cần tiến hành chọn tạo dòng vịt Biển: dòng HY1 - dòng trống,... protein nuôi vịt thịt thương phẩm hiệu kinh tế chăn nuôi vịt Biển Về vấn đề chọn lọc nhân giống vịt Biển: Trung tâm Nghiên cứu Vịt Đại Xuyên thực việc chọn lọc hệ vịt Biển 15 – Đại Xuyên theo phương... qua hệ 92 4 .15 Năng suất trứng 52 tuần đẻ HY2 qua hệ 92 ix TRÍCH YẾU LUẬN ÁN Tên tác giả: Chu Hoàng Nga Tên Luận án: Chọn tạo hai dòng vịt Biển sở giống vịt Biển 15 - Đại Xuyên Ngành: Chăn