Original Research ARTICLE
The Genetic Architecture of the Chickens Dropping Moisture by Genetic Parameter Estimation and Genome-Wide Association
- 1State Key Laboratory of Animal Nutrition, Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
- 2Hebei Dawu Poultry Breeding Co., Ltd., Hebei, China
- 3Beijing Animal Husbandry and Veterinary Station, Beijing, China
- 4Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
Dropping moisture (DM) refers to the water content of feces. High DM in chickens could be disadvantageous to pathogen control and fecal treatment in chicken farms. DM can be affected by environment, nutrition, disease, and genetics. In the present study, significant individual differences were presented in the DM of Rhode Island Red (RIR) chicken population, indicating that genetics could contribute to DM in the chickens. Subsequently, we estimated the genetic parameters of DM and conducted a genome-wide association study (GWAS) to find the potential genomic regions related to DM. The results showed that the heritability of DM ranged from 0.25 to 0.32. Furthermore, 11 significant loci on chromosome 7 were found to be associated with DM levels by the GWAS. The SNP rs15833816 within the COL6A3 gene was the most significant SNP related to DM. Hens carrying the G allele including GA and GG produced higher DM (P < 0.01) levels than those carrying the other genotype AA. Our results showed that DM is a medium-inheritable trait and that COL6A3 could be a potential candidate gene that regulates DM level in chickens.
The cloaca in laying hens is a passage at the end of the digestive, urinary, and reproductive systems and includes the urethra, ejaculation canal, and anal canal. Because of the unique excretory organ of chickens, feces and urine are mixed together and excreted through the cloaca, and the moisture content of hen feces is typically higher than that of some other animals.
Excess dropping moisture (DM) can severely affect the economic benefits of the laying hen industry. First, high DM can reduce the digestibility of chicken. Next, a higher DM will cause problems for fecal cleanup and fecal composting, thus eventually increasing labor costs (Ostrander and Hart, 1964). In addition, water-like feces attract more insects and bacteria in the summer, and high-DM feces are prone to releasing NH3, H2S, and other harmful gases that threaten the biological safety of laying hens (Jongebreur and Monteny, 2001; Achiano and Giliomee, 2005).
The excessive moisture content of the layer’s feces could be caused by many factors including genetics, disease, and environment. Anatomically, a chicken’s digestive system is relatively short and less time is required to digest feed in the gastrointestinal tract. Owing to the special digestive system of chickens, the DM is common in chickens. In addition, bacteria such as Escherichia coli (Li et al., 2012) and Pasteurella (Thanasarasakulpong et al., 2015) and viruses such as rotavirus (Kuroki et al., 1993; Wickelgren, 2000; Lim et al., 2005; Alam et al., 2011) and Newcastle disease virus (Kapczynski et al., 2013) can cause severe diarrhea in chickens.
In this study, we analyzed DM genetically. To understand the genetics of DM in chickens, we estimated the genetic parameters of DM in a layer population, i.e., Rhode Island Red (RIR) chickens, and performed a genome-wide association study (GWAS) to identify the potential genomic regions that regulate DM in the chickens.
Materials and Methods
Animals and Data Collection
Our experimental RIR hens were obtained from a layer breeding company in China. All hens were raised in three-tier H-shaped single cages. In total, 2,500 hens from 112 half-sib families were used to record the DM for each chicken. Observers started to record the DM when the RIR hens were 45 weeks old. Because Salmonella Pullorum (SP) and avian leukocyte virus (ALV) are common in China, we checked for the SP antibody and ALV p27 antigen titers in our birds. We confirmed that our chickens were all free of the two pathogens.
Because it is difficult to detect the real water content of the feces for each of the chickens, we developed an easy and efficient method to measure the relative DM levels. By this method, the DM levels were evaluated by the appearance of the feces. The lower the DM levels, the drier were the feces. Therefore, we recorded the DM levels by numbers from 1 to 4, each corresponding to the water content of the feces, i.e., normal, slight, medium, and severe DM, respectively (Figure 1). The grading standards of our method in assessing the DM levels are shown in Table 1. The DM levels of the chickens were observed every 3 days for five times by two observers. To keep the consistency of evaluation, two observers were trained with 100 chickens until the observation results are highly consistent, and the data were reviewed during data collection.
In order to study the relationship of DM to cage height and chicken performance, cage height, egg number (EN), and egg weight (EW) at age 48 weeks were recorded for further analysis.
Genetics Parameters Estimation
Two strategies were used in our study to estimate the genetic component variance. Strategy 1: The variance component of DM level was estimated based on the daily record; observer was included in module as a fixed variable. The variance component was estimated by Eq. 1; Strategy 2: To avoid the artificial bias, DM level numbers were added for the five times to represent the severity of DM (SDM) for each chicken, consider that SDM is not normal distribution as SDM is an integer between 5 and 20, SDM was divided into norm (SDM ≤ 6), slight (7 ≤ SDM ≤ 9), medium (10 ≤ SDM ≤ 14), and severe (SDM ≥ 15) according to the SDM distribution (Figure 2), and then the variance component was estimated Eq. 2.
Figure 2. The distribution of the total SDM. The DM state was divided into four categories according to the red dotted line.
where Y is the vector of the daily DM level; b1 represents the fixed effects including the observer and cage height; b2 is the fixed effects including the cage height; e is the vector of random residuals; [X] and [Z] are incidence matrices for the fixed and random effects, respectively. We assumed that [a]∼N([O,Aσ2]), [pe]∼N(0,[Iσ2]), and [e]∼N(0,[Iσ2]). The component estimation was performed using the DMU (v6) package (Madsen et al., 2006).
Genome-Wide Association Study
We performed GWAS in order to identify the DM related genes. Individuals with SDMs of five and 20 were grouped into case and control groups, respectively. Finally, we selected 48 hens, i.e., 25 cases and 23 controls, for genome-wide genotyping. Genomic DNA was extracted from blood samples using the TIANamp Blood DNA Kit [Tiangen Biotech (Beijing) Co., Ltd.], DNA quality was tested by NanoDrop 2000 (Thermo Scientific, United States), and all the DNA samples were tested as qualified (Supplementary Table S1). Then the DNA samples were genotyped on the 600K Affymetrix Axiom chicken panel (Kranis et al., 2013), which is designed based on reference genome Gallus_gallus-5.0 and contains 580,961 SNPs with an average probe spacing of 1.79 kb. Calling and quality control was performed in the Axiom Analysis Suite; only autosomal markers with clear physical position information were used in the analysis. The subsequently quality control was performed by PLINK v1.90. SNPs were filtered based on the call rate (geno > 0.90), minor allele frequency (maf > 0.01), and Hardy–Weinberg equilibrium (hwe > 10–5). The samples were filtered based on the call rate (mind > 0.9). In filtration, 356,008 variants and all the 48 chicken samples passed quality control. Then, a chi-square-based genome association analysis was performed for all samples. Bonferroni correction was applied to account for multiple testing. Consider that many SNPs fall within regions of strong linkage disequilibrium (LD) (“blocks”), the suggestive and significant threshold was set at 1/LD blocks number and 0.05/LD blocks (Duggal et al., 2008). The LD blocks number was predicted by PLINK v1.90.
DM Is a Medium-Heritable Trait
A total of 2,500 hens were phenotyped as candidate experimental individuals in our study (Table 2). According to the DMU instruction, the Gibbs sampling method was used to estimate the genetic parameter of DM. For strategy 1, the heritability of DM ranged from 0.25 to 0.31 (Table 3). For strategy 2, the heritability of DM was 0.32 (Table 3). Our results show that DM is a heritable trait with medium heritability.
DM Correlates With Cage Height
In our study, the chi-square test was used to test the effect of cage height on chicken DM, and the results showed that DM was related to cage height (P < 0.01). The first-level cages not only exhibited a lower proportion of slight DM but also higher proportions of middle and severe DM (Figure 3).
Figure 3. The distribution of DM state in different cages, where L1 is first floor, L2 is second floor, and L3 is third floor.
On the other hand, we counted the EN and randomly selected the eggs to measure the average EW at the age of 48 weeks. The EN ranged from 160 to 180, and the EW ranged from 52 to 66 g. One-way ANOVA analysis showed that EN nor EW was correlated with DM (P > 0.1).
Candidate Regions Identified to Be Related to DM by GWAS
We selected 25 chickens with a total SDM of five in the case group and 23 with total SDM of 20 in the control group to perform GWA. The top-3 principal components explained more than 60% of the variance between the individuals. No significant population stratification was found in our experiment (Figure 4). A total of 18,650 LD blocks were found in the autosomes; therefore, the thresholds and significant threshold of the association analysis were set to 5.35 × 10–5 (1/18,650) and 2.68 × 10–6 (0.05/18,650), respectively. The case-control GWAS was performed for all individuals, and the results are presented in Figure 5A. Eleven regions (SNPs) were beyond the suggested line; however, only one SNP on chromosome 7 remained significant after Bonferroni correction. There is a peak around the SNP, 10 SNPs beyond the suggested line (Figure 5B). This peak overlapped with three genes, COL6A3, MLPH, and RAB17 gene (Figure 5C). The detail SNPs information is presented in Table 4. The genotype frequencies of rs15833816 presented the most significant difference between the case and control groups, which is an intron variation of the COL6A3 (collagen type VI alpha 3 chain) gene. Then gene ontology (GO) was employed for other autosome gene which overlap with significant SNP, but no GO item was enriched, and the gene set and SNPs set was present in Supplementary Table S2.
Figure 5. The result of genome wide association. (A) Manhattan plot of genome-wide association study. (B) Manhattan plot from 2M to 6M in chromosome 7. (C) The overlapping genes of significant peak on chromosome 7.
Heritability of DM
The estimation of variance component is largely dependent on the quality and quantity of data. DM level can be quantified by measuring the water content of feces of each bird. However, it was not feasible to measure the feces of each bird in our experiment. Therefore, we regarded DM as a threshold trait by assigning numbers from 1 to 4 according to the feces shape of each chicken. A similar strategy has been used in other animal studies, such as feather peaking (Brinker et al., 2014) and body condition score of cows (Edmonson et al., 1989). In the present study, the heritability in different models ranged from 0.25 to 0.32, indicating that the additive effect explained at least 25% of the total phenotypic variance. Genetic selection could be an effective method for changing the DM level in the chickens.
A previous study reported that cage height affected the behavior of hens leading to head scratching, and body shaking (Nicol, 1987; Albentosa et al., 2007). In this study, we found that cage height was strongly related to DM level. Chickens in lower cages presented higher DM levels, suggesting that sanitary conditions caused by cage height affected DM levels. There are many environmental factors decided by cage height, such as air condition, number of microorganisms, and temperature. Improving housing conditions is also an effective way to eliminate DM.
Genome-Wide Association Study
With the development of sequencing technology and microarray technology, the cost of genome-wide genotyping is becoming lower. GWAS, an effective method to identify important molecular markers, has been widely used in animal breeding.
In our study, we identified a marker that is significantly related to DM level and three candidate regulatory genes (CLO6A3, MLPH, and RAB17). Considering that the most significant SNPs occurred in the COL6A3 gene, this gene is most likely associated with DM in chickens. The COL6A3 gene starts at 4,808,221 bp on 7 chromosome and ends at 4,861,524 bp in Gallus_gallus-5.0, and its full length is 53,303 bp. As reported in a previous study, mutations in the COL6A3 gene can lead to poor muscle development in humans and cause Down’s syndrome (Dey et al., 2013) and muscle dystonia (Demir et al., 2002). The full name of MLPH gene is melanophilin, which has been found to be related to the feather color in chickens (Vaez et al., 2008). RAB17 gene, full name is member RAS oncogene family, which is a cancer- and immunity-related gene.
Dropping moisture is a medium heritable trait and a genetic marker exhibited significant association with it. Furthermore, we found that the main candidate gene COL6A3 could affect the DM level in chickens. This provides theoretical basis for subsequent functional verification and insight into genetic selection for the DM level to help improve the economic efficiency of layer farms.
Data Availability Statement
The datasets that support the findings of this study are available on FigShare https://figshare.com/articles/GENETICS_OF_DROPPING_MOISTURE/12027936.
The animal study was reviewed and approved by the Committee for Animal Care and Use of China Agricultural University. Written informed consent was obtained from the owners for the participation of their animals in this study.
TZ and T-YZ performed the experiments and data analysis and wrote the manuscript. JW and ZG contributed to sample collection. XZ provided test poultry farms. YJ, YC, LW, XL, WY, and LQ designed the experiments and coordinated the study. LQ and ZN reviewed and revised the manuscript. All authors approved the submitted version.
This work was supposed by Beijing Innovation Team of the Modern Agro-industry Technology Research System (BAIC04-2019).
Conflict of Interest
XZ was employed by the company Hebei Dawu Poultry Breeding Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We sincerely appreciate Xiaoyu Zhao for access to her chickens, and staff at Hebei Dawu Poultry Breeding Co., Ltd., for their assistance, particularly Hao Zhou and Tianzeng Li.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00806/full#supplementary-material
Alam, M. N., Alam, M. M., Nahar, A., and Kobayashi, N. (2011). Molecular epidemiological studies on rotavirus infection causing severe diarrhea in human, animals and poultry. ID - 20133202532, Bangladesh. J Vet Med. 9, 167–175. doi: 10.3329/bjvm.v9i2.13473
Albentosa, M. J., Cooper, J. J., Luddem, T., Redgate, S. E., Elson, H. A., and Walker, A. W. (2007). Evaluation of the effects of cage height and stocking density on the behaviour of laying hens in furnished cages. Br. Poultry Sci. 48, 1–11. doi: 10.1080/00071660601156479
Brinker, T., Bijma, P., Visscher, J., Rodenburg, T. B., and Ellen, E. D. (2014). Plumage condition in laying hens: genetic parameters for direct and indirect effects in two purebred layer lines. Genet. Sel. Evol. 46:33. doi: 10.1186/1297-9686-46-33
Demir, E., Sabatelli, P., Allamand, V., Ferreiro, A., Moghadaszadeh, B., Makrelouf, M., et al. (2002). Mutations in COL6A3 cause severe and mild phenotypes of ullrich congenital muscular dystrophy. Am. J. Hum. Genet. 70, 1446–1458. doi: 10.1086/340608
Dey, A., Bhowmik, K., Chatterjee, A., Sinha, S., and Mukhopadhyay, K. (2013). Down syndrome related muscle hypotonia: association with COL6A3 functional SNP rs2270669. Front. Genet. 4:57. doi: 10.3389/fgene.2013.00057
Duggal, P., Gillanders, E. M., Holmes, T. N., and Bailey-Wilson, J. E. (2008). Establishing an adjusted p-value threshold to control the family-wide type 1 error in genome wide association studies. BMC Genomics 9:516. doi: 10.1186/1471-2164-9-516
Edmonson, A. J., Lean, I. J., Weaver, L. D., Farver, T., and Webster, G. (1989). A body condition scoring chart for holstein dairy cows. J. Dairy Sci. 72, 68–78. doi: 10.3168/jds.S0022-0302(89)79081-0
Jongebreur, A. A., and Monteny, G. J. (2001). Prevention and control of losses of gaseous nitrogen compounds in livestock operations: a review. Sci. World J. 1(Suppl. 2), 844–851. doi: 10.1100/tsw.2001.339
Kranis, A., Gheyas, A. A., Boschiero, C., Turner, F., Yu, L., Smith, S., et al. (2013). Development of a high density 600K SNP genotyping array for chicken. BMC Genomics 14:59. doi: 10.1186/1471-2164-14-59
Kuroki, M., Ikemori, Y., Yokoyama, H., Peralta, R. C., Icatlo, F. C., and Kodama, Y. (1993). Passive protection against bovine rotavirus-induced diarrhea in murine model by specific immunoglobulins from chicken egg yolk. Vet. Microbiol. 37, 135–146. doi: 10.1016/0378-1135(93)90188-D
Li, D. F., Liu, W. B., Liu, J. F., Yi, G. Q., Lian, L., Qu, L. J., et al. (2012). Whole-genome scan for signatures of recent selection reveals loci associated with important traits in White Leghorn chickens. Poultry Sci. 91, 1804–1812. doi: 10.3382/ps.2012-02275
Lim, I. S., Lee, H. S., Kim, W., Choi, E. S., Jung, D. H., Jung, H. K., et al. (2005). Therapeutic effect of anti-rotavirus chicken egg yolk immunoglobulin(IgY) on diarrhea by infection of rotavirus. Korean J. Pediatr. 48, 1354–1361.
Madsen, P., Sørensen, P., Su, G., Damgaard, L. H., Thomsen, H., and Labouriau, R. (2006). DMU - A Package for Analyzing Multivariate Mixed Models. ID - 20063170093. Minas Gerais: Instituto Prociência, 11–27.
Thanasarasakulpong, A., Poolperm, P., Tankaew, P., Sawada, T., and Sthitmatee, N. (2015). Protectivity conferred by immunization with intranasal recombinant outer membrane protein H from Pasteurella multocida serovar A:1 in chickens. J. Vet. Med. Sci. 77, 321–326. doi: 10.1292/jvms.14-0532
Vaez, M., Follett, S. A., Bed’hom, B., Gourichon, D., Tixier-Boichard, M., and Burke, T. (2008). A single point-mutation within the melanophilin gene causes the lavender plumage colour dilution phenotype in the chicken. BMC Genet. 9:7. doi: 10.1186/1471-2156-9-7
Keywords: chicken, dropping moisture, genetic parameters, GWAS, COL6A3 gene
Citation: Zhu T, Zhang T-Y, Wen J, Zhao X, Chen Y, Jia Y, Wang L, Lv X, Yang W, Guan Z, Ning Z and Qu L (2020) The Genetic Architecture of the Chickens Dropping Moisture by Genetic Parameter Estimation and Genome-Wide Association. Front. Genet. 11:806. doi: 10.3389/fgene.2020.00806
Received: 13 March 2020; Accepted: 06 July 2020;
Published: 29 July 2020.
Edited by:Xianyao Li, Shandong Agricultural University, China
Reviewed by:Guirong Sun, Henan Agricultural University, China
Xiaozhu Wang, Auburn University, United States
Copyright © 2020 Zhu, Zhang, Wen, Zhao, Chen, Jia, Wang, Lv, Yang, Guan, Ning and Qu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work