Genome-wide association study of growth traits in the Jinghai Yellow chicken

Genetic parameter estimation and genetic evaluation of important economic traits in white and yellow broilers

1. This study calculated descriptive statistics for the production traits of two broiler populations: 1) the Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF white broilers), including fat and lean lines; and 2) the Guangxi yellow broilers. Their genetic parameters were estimated, including (co)variance components, heritability (h2) and genetic correlations (rg), using the REML method.2. Heritability estimates (h2) for NEAUHLF white broilers ranged from 0.07 to 0.61. Traits with high heritability (h2 >0.3) included body weight at 3, 5 and 7 weeks of age (BW3, BW5, BW7), carcass weight (CW), metatarsal circumference (MeC), liver weight (LW), gizzard weight (GW), spleen weight (SW) and testis weight (TeW), while in Guangxi yellow broilers, heritability estimates ranged from 0.18 to 0.76, with every trait exhibiting high heritability, except for SW (0.18).3. Positive genetic correlations for NEAUHLF were found (rg >0.3, ranging from 0.31 to 0.84) between BW7 and metatarsal length (MeL), MeC, body oblique length (BoL), chest angle (ChA), LW, GW, heart weight (HW) and SW. Genetic correlations between abdominal fat weight (AFW) and BW1, BW3, BW5, CW, MeL, keel length (KeL), BoL and LW were positive (rg >0.3, ranging from 0.31 to 0.58).4. Among the Guangxi population, BW (125 d of age) showed strong positive genetic correlations with all other traits (rg >0.3, ranging from 0.30 to 0.99), while AFW displayed strong positive genetic correlations with leg muscle weight (LeW), CW, BW and thigh diameter (TD) (rg >0.3, ranging from 0.44 to 0.51).5. It was concluded that the characteristics of the two populations were different, which means there is a need to use different strategies when performing the breeding work to improve productivity and efficiency in both broiler populations.

Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

Genome-Wide Association Studies Provide Insight Into the Genetic Determination for Hyperpigmentation of the Visceral Peritoneum in Broilers

Hyperpigmentation of the visceral peritoneum (HVP) has been becoming one of the most challenging problems in yellow-feathered chicken production, which seriously affected chicken carcass quality traits. Detecting which genes dominantly impact pigmentation in the peritoneum tissues is of great benefit to the genetic improvement of HVP. To investigate the genetic mechanism of HVP in yellow-feathered broilers, genome-wide association studies (GWASs) were conducted in the F₂ generation of a cross broiler population with 395 birds. A total of 115,706 single-nucleotide polymorphisms (SNPs) of 122,415 were retained to identify quantitative trait loci (QTL) associated to HVP in chicken. The GWAS results based on the logistic mixed model (LMM) revealed that a narrow genomic location on chromosomes 1 (49.2–51.3 Mb) was significantly associated (p ≤ 4.32 × 10⁻⁷) with HVP, which contained 23 SNP makers related to 14 functional genes (MFNG, POLDIP3, POLR2F, PICK1, PDXP, SGSM3, RANGAP1, MYH9, RPL3, GALP3, LGALS1, MICALL1, ATF4, and CYP2D6). Four highly associated (p < 10⁻⁵) haplotype blocks of 0.80 kb (two SNPs), 0.06 kb (two SNPs), 0.95 kb (two SNPs), and 0.03 kb (two SNPs) were identified with two, two, four, and four distinct haplotypes, respectively. As a melanoma-associated gene, CYP2D6 were also possibly involved in the development of HVP occurring in chicken with two significant variations (rs314284996 and rs317955795) in the promoter regions. Further tests revealed that the expression of CYP2D6 was obviously higher in the visceral peritoneum tissue of chicken with HVP than that in the normal group (p < 0.05). Our results provide a novel clue to understand the genetic mechanism of HVP generation in chicken, and the mapped QTL or candidate genes might serve for genomic selection to improve carcass quality in the yellow-feathered chicken industry.

Genome-Wide Association Study of Growth Performance and Immune Response to Newcastle Disease Virus of Indigenous Chicken in Rwanda

A chicken genome has several regions with quantitative trait loci (QTLs). However, replication and confirmation of QTL effects are required particularly in African chicken populations. This study identified single nucleotide polymorphisms (SNPs) and putative genes responsible for body weight (BW) and antibody response (AbR) to Newcastle disease (ND) in Rwanda indigenous chicken (IC) using genome-wide association studies (GWAS). Multiple testing was corrected using chromosomal false detection rates of 5 and 10% for significant and suggestive thresholds, respectively. BioMart data mining and variant effect predictor tools were used to annotate SNPs and candidate genes, respectively. A total of four significant SNPs (rs74098018, rs13792572, rs314702374, and rs14123335) significantly (p ≤ 7.6E-5) associated with BW were identified on chromosomes (CHRs) 8, 11, and 19. In the vicinity of these SNPs, four genes such as pre-B-cell leukaemia homeobox 1 (PBX1), GPATCH1, MPHOSPH6, and MRM1 were identified. Four other significant SNPs (rs314787954, rs13623466, rs13910430, and rs737507850) all located on chromosome 1 were strongly (p ≤ 7.6E-5) associated with chicken antibody response to ND. The closest genes to these four SNPs were cell division cycle 16 (CDC16), zinc finger, BED-type containing 1 (ZBED1), myxovirus (influenza virus) resistance 1 (MX1), and growth factor receptor bound protein 2 (GRB2) related adaptor protein 2 (GRAP2). Besides, other SNPs and genes suggestively (p ≤ 1.5E-5) associated with BW and antibody response to ND were reported. This work offers a useful entry point for the discovery of causative genes accountable for essential QTLs regulating BW and antibody response to ND traits. Results provide auspicious genes and SNP-based markers that can be used in the improvement of growth performance and ND resistance in IC populations based on gene-based and/or marker-assisted breeding selection.

Association of growth hormone and insulin-like growth factor I genotype with body weight, dominance of body weight, and mRNA expression in Korat slow-growing chickens

Objective

Growth hormone (GH) and insulin-like growth factor I (IGF-I) play a critical role in animal growth rates. We aimed to investigate the effect of GH and IGF-I genotypes on body weight (BW), dominance, and gene expression in slow-growing chickens at different ages.

Methods

A total of 613 Korat chickens (KRs) were bred and divided into three groups by genotype – A1A1, A1A3, and A3A3 for GH and AA, AC, and CC for IGF-I. Chickens were weighed every two weeks, and liver and breast muscle tissues were collected at 10 weeks of age. Genetic parameters of KRs were estimated using ASReml software. The GH and IGF-I mRNA levels were measured by quantitative polymerase chain reaction. Significant differences between traits were analyzed using the generalized linear model.

Results

A significant effect of GH genotypes on BW was found at most ages, and the A1A1 genotype had the highest value of BW. Compared with the A3A3 genotype, the A1A1 and A1A3 genotypes showed a higher dominance effect at 0 and 2 weeks, and genotype A1A1 had the highest value of dominance at 8 weeks of age. A difference in GH mRNA levels between genotypes was detected in breast muscle at 6 weeks and in the liver tissue at 2 weeks. In the case of IGF-I gene, the AA genotype had the highest BW at the beginning of life. Significant differences in BW dominance were found at 2 weeks. However, IGF-I mRNA levels were not different among genotypes in both breast muscles and liver tissues.

Conclusion

Our results revealed that GH and IGF-I influence growth, but may not be involved in heterosis. GH can be used as a marker gene in selection programs for growth because the homozygous genotype (A1A1) had the highest BW at all ages. The IGF-I is not a useful marker gene for selection programs.

MicroRNA-7 Targets the KLF4 Gene to Regulate the Proliferation and Differentiation of Chicken Primary Myoblasts

The proliferation and differentiation of chicken primary myoblasts (CPMs) play an important role in the development of skeletal muscle. In our previous research, RNA-seq analysis showed that microRNA-7 (miR-7) was relatively highly expressed in the proliferation phase of CPMs, but its expression level decreased significantly after CPM_S-induced differentiation. Meanwhile, the mechanism by which the miR-7 regulates the proliferation and differentiation of CPMs is still unknown. In this study, we found that the expression levels of miR-7 and the Krüppel-like factor 4 (KLF4) gene were negatively correlated during the embryonic phase, and in vitro induced differentiation. A dual-luciferase assay and a rescue experiment show that there is a target relationship between miR-7 and the KLF4 gene. Meanwhile, the results show that overexpression of miR-7 inhibited the proliferation and differentiation of CPMs, while inhibition of miR-7 had the opposite effects. Furthermore, overexpression of the KLF4 gene was found to significantly promote the proliferation and differentiation of CPMs. Conversely, inhibition of the KLF4 gene was able to significantly decrease the proliferation and differentiation of CPMs. Our results demonstrate, for the first time, that miR-7 inhibits the proliferation and differentiation of myoblasts by targeting the KLF4 gene in chicken primary myoblasts.

Transcriptome Analysis of Differentially Expressed Genes Related to the Growth and Development of the Jinghai Yellow Chicken

The growth traits are important traits in chickens. Compared to white feather broiler breeds, Chinese local broiler breeds have a slow growth rate. The main genes affecting the growth traits of local chickens in China are still unclear and need to be further explored. This experiment used fast-growth and slow-growth groups of the Jinghai Yellow chicken as the research objects. Three males and three females with similar body weights were selected from the two groups at four weeks old and eight weeks old, respectively, with a total of 24 individuals selected. After slaughter, their chest muscles were taken for transcriptome sequencing. In the differentially expressed genes screening, all of the genes obtained were screened by fold change ≥ 2 and false discovery rate (FDR) < 0.05. For four-week-old chickens, a total of 172 differentially expressed genes were screened in males, where there were 68 upregulated genes and 104 downregulated genes in the fast-growth group when compared with the slow-growth group. A total of 31 differentially expressed genes were screened in females, where there were 11 upregulated genes and 20 downregulated genes in the fast-growth group when compared with the slow-growth group. For eight-week-old chickens, a total of 37 differentially expressed genes were screened in males. The fast-growth group had 28 upregulated genes and 9 downregulated genes when compared with the slow-growth group. A total of 44 differentially expressed genes were screened in females. The fast-growth group had 13 upregulated genes and 31 downregulated genes when compared with the slow-growth group. Through gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, many genes were found to be related to cell proliferation and differentiation, muscle growth, and cell division such as SNCG, MCL1, ARNTL, PLPPR4, VAMP1, etc. Real-time PCR results were consistent with the RNA-Seq data and validated the findings. The results of this study will help to understand the regulation mechanism of the growth and development of Jinghai Yellow chicken and provide a theoretical basis for improving the growth rate of Chinese local chicken breeds.

Discovery and characterization of functional modules associated with body weight in broilers

Aim of the present study was to investigate whether body weight (BW) in broilers is associated with functional modular genes. To this end, first a GWAS for BW was conducted using 6,598 broilers and the high density SNP array. The next step was to search for positional candidate genes and QTLs within strong LD genomic regions around the significant SNPs. Using all positional candidate genes, a network was then constructed and community structure analysis was performed. Finally, functional enrichment analysis was applied to infer the functional relevance of modular genes. A total number of 645 positional candidate genes were identified in strong LD genomic regions around 11 genome-wide significant markers. 428 of the positional candidate genes were located within growth related QTLs. Community structure analysis detected 5 modules while functional enrichment analysis showed that 52 modular genes participated in developmental processes such as skeletal system development. An additional number of 14 modular genes (GABRG1, NGF, APOBEC2, STAT5B, STAT3, SMAD4, MED1, CACNB1, SLAIN2, LEMD2, ZC3H18, TMEM132D, FRYL and SGCB) were also identified as related to body weight. Taken together, current results suggested a total number of 66 genes as most plausible functional candidates for the trait examined.

Derivation and genome-wide association study of a principal component-based measure of heat tolerance in dairy cattle

Heat stress represents a key factor that negatively affects the productive and reproductive performance of farm animals. In the present work, a new measure of tolerance to heat stress for dairy cattle was developed using principal component analysis. Data were from 590,174 test-day records for milk yield, fat and protein percentages, and somatic cell score of 39,261 Italian Holstein cows. Test-day records adjusted for main systematic factors were grouped into 11 temperature-humidity index (THI) classes. Daughter trait deviations (DTD) were calculated for 1,540 bulls as means of the adjusted test-day records for each THI class. Principal component analysis was performed on the DTD for each bull. The first 2 principal components (PC) explained 42 to 51% of the total variance of the system across the 4 traits. The first PC, a measure of the level at which the curve is located, was interpreted as a measure of the level at which the DTD curve was located. The second PC, which shows the slope of increasing or decreases DTD curves, synthesized the behavior of the DTD pattern. Heritability of the 2 component scores was moderate to high for level across all traits (range = 0.23-0.82) and low to moderate for slope (range = 0.16-0.28). For each trait, phenotypic and genetic correlations between level and slope were equal to zero. A genome-wide association analysis was carried out on a subsample of 423 bulls genotyped with the Illumina 50K bovine bead chip (Illumina, San Diego, CA). Two single nucleotide polymorphisms were significantly associated with slope for milk yield, 4 with level for fat percentage, and 2 with level and slope of protein percentage, respectively. The gene discovery was carried out considering windows of 0.5 Mb surrounding the significant markers and highlighted some interesting candidate genes. Some of them have been already associated with the mechanism of heat tolerance as the heat shock transcription factor (HSF1) and the malonyl-CoA-acyl carrier protein transacylase (MCAT). The 2 PC were able to describe the overall level and the slope of response of milk production traits across increasing levels of THI index. Moreover, they exhibited genetic variability and were genetically uncorrelated. These features suggest their use as measures of thermotolerance in dairy cattle breeding schemes.

Transcriptomic profile of leg muscle during early growth in chicken

The early growth pattern, especially the age of peak growth, of broilers affects the time to market and slaughter weight, which in turn affect the profitability of the poultry industry. However, the underlying mechanisms regulating chicken growth and development have rarely been studied. This study aimed to identify candidate genes involved in chicken growth and investigated the potential regulatory mechanisms of early growth in chicken. RNA sequencing was applied to compare the transcriptomes of chicken muscle tissues at three developmental stages during early growth. In total, 978 differentially expressed genes (DEGs) (fold change ≥ 2; false discovery rate < 0.05) were detected by pairwise comparison. Functional analysis showed that the DEGs are mainly involved in the processes of cell growth, muscle development, and cellular activities (such as junction, migration, assembly, differentiation, and proliferation). Many of the DEGs are well known to be related to chicken growth, such as MYOD1, GH, IGF2BP2, IGFBP3, SMYD1, CEBPB, FGF2, and IGFBP5. KEGG pathway analysis identified that the DEGs were significantly enriched in five pathways (P < 0.1) related to growth and development: extracellular matrix–receptor interaction, focal adhesion, tight junction, insulin signaling pathway, and regulation of the actin cytoskeleton. A total of 42 DEGs assigned to these pathways are potential candidate genes inducing the difference in growth among the three developmental stages, such as MYH10, FGF2, FGF16, FN1, CFL2, MAPK9, IRS1, PHKA1, PHKB, and PHKG1. Thus, our study identified a series of genes and several pathways that may participate in the regulation of early growth in chicken. These results should serve as an important resource revealing the molecular basis of chicken growth and development.