Genome-wide association studies for growth traits in broilers

Non-Targeted Metabolomics of Serum Reveals Biomarkers Associated with Body Weight in Wumeng Black-Bone Chickens

Growth performance is an important economic trait of broilers but the related serum metabolomics remains unclear. In this study, we utilized non-targeted metabolomics using ultra-high-performance liquid phase tandem mass spectrometry (UHPLC-MS/MS) to establish metabolite profiling in the serum of Chinese Wumeng black-bone chickens. The biomarker metabolites in serum associated with growth performance of chickens were identified by comparing the serum metabolome differences between chickens that significantly differed in their weights at 160 days of age when fed identical diets. A total of 766 metabolites were identified including 13 differential metabolite classes such as lipids and lipid-like molecules, organic acids and their derivatives, and organoheterocyclic compounds. The results of difference analysis using a partial least squares discriminant analysis (PLS-DA) model indicated that the low-body-weight group could be differentiated based on inflammatory markers including prostaglandin a2, kynurenic acid and fatty acid esters of hydroxy fatty acids (FAHFA), and inflammation-related metabolic pathways including tryptophan and arachidonic acid metabolism. In contrast, the sera of high-body-weight chickens were enriched for riboflavin and 2-isopropylmalic acid and for metabolic pathways including riboflavin metabolism, acetyl group transfer into mitochondria, and the tricarboxylic acid (TCA) cycle. These results provide new insights into the practical application of improving the growth performance of local chickens.

Single-nucleus RNA sequencing reveals cell types, genes, and regulatory factors influencing melanogenesis in the breast muscle of Xuefeng black-bone chicken

The black-bone chicken, known for its high melanin content, holds significant economic value due to this unique trait. Particularly notable is the prominent melanin deposition observed in its breast muscle. However, the molecular mechanisms governing melanin synthesis and deposition in the breast muscle of black-bone chickens remain largely unknown. This study employed a single-nucleus transcriptome assay to identify genes associated with melanin deposition in the breast muscle of black-bone chickens, which are presumed to influence pigmentation levels. A comprehensive analysis of the nuclear transcriptome was conducted on the breast muscle of Xuefeng black-bone chickens, encompassing 18 distinct cell types, including melanocytes. Our findings revealed that STIMATE, LRRC7, ENSGALG00000049990, and GLDC play pivotal regulatory roles in melanin deposition within the breast muscle. Further exploration into the molecular mechanisms unveiled transcription factors and protein interactions suggesting that RARB, KLF15, and PRDM4 may be crucial regulators of melanin accumulation in the breast muscle. Additionally, HPGDS, GSTO1, and CYP1B1 may modulate melanin production and deposition in the breast muscle by influencing melanocyte metabolism. Our findings also suggest that melanocyte function in the breast muscle may be intertwined with intercellular signaling pathways such as PTPRK-WNT5A, NOTCH1-JAG1, IGF1R-IGF1, IDE-GCG, and ROR2-WNT5A. Leveraging advanced snRNA-seq technology, we generated a comprehensive single-cell nuclear transcriptome atlas of the breast muscle of Xuefeng black-bone chickens. This facilitated the identification of candidate genes, regulatory factors, and cellular signals potentially influencing melanin deposition and melanocyte function. Overall, our study provides crucial insights into the molecular basis of melanin deposition in chicken breast muscle, laying the groundwork for future breeding programs aimed at enhancing black-bone chicken cultivation.

Application of GWAS and mGWAS in Livestock and Poultry Breeding

In recent years, genome-wide association studies (GWAS) and metabolome genome-wide association studies (mGWAS) have emerged as crucial methods for investigating complex traits in animals and plants. These have played pivotal roles in research on livestock and poultry breeding, facilitating a deeper understanding of genetic diversity, the relationship between genes, and genetic bases in livestock and poultry. This article provides a review of the applications of GWAS and mGWAS in animal genetic breeding, aiming to offer reference and inspiration for relevant researchers, promote innovation in animal genetic improvement and breeding methods, and contribute to the sustainable development of animal husbandry.

Genome-wide association analysis identify candidate genes for feed efficiency and growth traits in Wenchang chickens

BACKGROUND

Wenchang chickens are one of the most popular local chicken breeds in the Chinese chicken industry. However, the low feed efficiency is the main shortcoming of this breed. Therefore, there is a need to find a more precise breeding method to improve the feed efficiency of Wenchang chickens. In this study, we explored important candidate genes and variants for feed efficiency and growth traits through genome-wide association study (GWAS) analysis.

RESULTS

Estimates of genomic heritability for growth and feed efficiency traits, including residual feed intake (RFI) of 0.05, average daily food intake (ADFI) of 0.21, average daily weight gain (ADG) of 0.24, body weight (BW) at 87, 95, 104, 113 days of age (BW87, BW95, BW104 and BW113) ranged from 0.30 to 0.44. Important candidate genes related to feed efficiency and growth traits were identified, such as PLCE1, LAP3, MED28, QDPR, LDB2 and SEL1L3 genes.

CONCLUSION

The results identified important candidate genes for feed efficiency and growth traits in Wenchang chickens and provide a theoretical basis for the development of new molecular breeding technology.

Dissecting Selective Signatures and Candidate Genes in Grandparent Lines Subject to High Selection Pressure for Broiler Production and in a Local Russian Chicken Breed of Ushanka

Breeding improvements and quantitative trait genetics are essential to the advancement of broiler production. The impact of artificial selection on genomic architecture and the genetic markers sought remains a key area of research. Here, we used whole-genome resequencing data to analyze the genomic architecture, diversity, and selective sweeps in Cornish White (CRW) and Plymouth Rock White (PRW) transboundary breeds selected for meat production and, comparatively, in an aboriginal Russian breed of Ushanka (USH). Reads were aligned to the reference genome bGalGal1.mat.broiler.GRCg7b and filtered to remove PCR duplicates and low-quality reads using BWA-MEM2 and bcftools software; 12,563,892 SNPs were produced for subsequent analyses. Compared to CRW and PRW, USH had a lower diversity and a higher genetic distinctiveness. Selective sweep regions and corresponding candidate genes were examined based on ZFST, hapFLK, and ROH assessment procedures. Twenty-seven prioritized chicken genes and the functional projection from human homologs suggest their importance for selection signals in the studied breeds. These genes have a functional relationship with such trait categories as body weight, muscles, fat metabolism and deposition, reproduction, etc., mainly aligned with the QTLs in the sweep regions. This information is pivotal for further executing genomic selection to enhance phenotypic traits.

Integrating miRNA and full-length transcriptome profiling to elucidate the mechanism of muscle growth in Muscovy ducks reveals key roles for miR-301a-3p/ANKRD1

BACKGROUND

The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis.

RESULTS

A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p.

CONCLUSIONS

miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms.

Identify Candidate Genes Associated with the Weight and Egg Quality Traits in Wenshui Green Shell-Laying Chickens by the Copy Number Variation-Based Genome-Wide Association Study

Copy number variation (CNV), as an essential source of genetic variation, can have an impact on gene expression, genetic diversity, disease susceptibility, and species evolution in animals. To better understand the weight and egg quality traits of chickens, this paper aimed to detect CNVs in Wenshui green shell-laying chickens and conduct a copy number variation regions (CNVRs)-based genome-wide association study (GWAS) to identify variants and candidate genes associated with their weight and egg quality traits to support related breeding efforts. In our paper, we identified 11,035 CNVRs in Wenshui green shell-laying chickens, which collectively spanned a length of 13.1 Mb, representing approximately 1.4% of its autosomal genome. Out of these CNVRs, there were 10,446 loss types, 491 gain types, and 98 mixed types. Notably, two CNVRs showed significant correlations with egg quality, while four CNVRs exhibited significant associations with body weight. These significant CNVRs are located on chromosome 4. Further analysis identified potential candidate genes that influence weight and egg quality traits, including FAM184B, MED28, LAP3, ATOH8, ST3GAL5, LDB2, and SORCS2. In this paper, the CNV map of the Wenshui green shell-laying chicken genome was constructed for the first time through population genotyping. Additionally, CNVRs can be employed as molecular markers to genetically improve chickens' weight and egg quality traits.

Genome-wide unraveling SNP pairwise epistatic effects associated with sheep body weight

Epistatic effects are an important part of the genetic effect of complex traits in livestock. In this study, we used 218 synthetic ewes from the Xinjiang Academy of Agricultural Reclamation in China to identify interacting paired with genome-wide single nucleotide polymorphisms (SNPs) associated with birth weight, weaning weight, and one-yearling weight. We detected 2 and 66 SNP-SNP interactions of sheep birth weight and weaning weight, respectively. No significant epistatic interaction of one-year-old body weight was detected. The genetic interaction of sheep body weight is dynamic and time-dependent. Most significant interactions of weaning body weight contributed 1% or higher. In the weaning weight trait, 66 significant SNP pairs consisted of 98 single SNPs covering 23 chromosomes, 5 of which were nonsynonymous SNPs (nsSNPs), resulting in single amino acid substitution. We found that genes that interact with transcription factors (TFs) are target genes for the corresponding TFs. Four epitron networks affecting weaning weight, including subnetworks of HIVEP3 and BACH2 transcription factors, constructed using significant SNP pairs, were also analyzed and annotated. These results suggest that transcription factors may play an important role in explaining epistatic effects. It provides a new idea to study the genetic mechanism of weight developing.

Genome-Wide Association Study of Egg Production Traits in Shuanglian Chickens Using Whole Genome Sequencing

Egg production is the most important economic trait in laying hens. To identify molecular markers and candidate genes associated with egg production traits, such as age at first egg (AFE), weight at first egg (WFE), egg weight (EW), egg number (EN), and maximum consecutive egg laying days (MCD), a genome-wide analysis by whole genome sequencing was performed in Shuanglian chickens. Through whole genome sequencing and quality control, a total of 11,006,178 SNPs were obtained for further analysis. Heritability estimates ranged from moderate to high for EW (0.897) and MCD (0.632), and from low to moderate (0.193~0.379) for AFE, WFE, and EN. The GWAS results showed 11 genome-wide significant SNPs and 23 suggestive significant SNPs were identified to be associated with EN, MCD, WFE, and EW. Linkage disequilibrium analysis revealed twenty-seven SNPs associated with EN were located in a 0.57 Mb region on GGA10, and clustered into five blocks. Through functional annotation, three candidate genes NEO1, ADPGK, and CYP11A1, were identified to be associated with EN, while the S1PR4, LDB2, and GRM8 genes was linked to MCD, WFE, and EW, respectively. These findings may help us to better understand the molecular mechanisms underlying egg production traits in chickens and contribute to genetic improvement of these traits.

Metabolomic characterization of Liancheng white and Cherry Valley duck breast meat and their relation to meat quality

Liancheng white duck is a typical local duck breed in Fujian Province famous for its meat traits. To better understand how meat quality varies with breed, the chemical composition of breast meats of Liancheng white ducks (LD) and Cherry Valley ducks (CD) were examined using ultra-high-performance liquid chromatography quadrupole time-of-flight mass spectrometry (UHPLC-QTOF-MS).The correlation between meat quality and the differential metabolites was further analyzed. The results showed that the effects of breed on duck breast meat were significant for pH, color, cooking loss, and shear force. Liancheng white duck breast meat exhibited a higher shear force and pH, and lower cooking loss and lightness (L*24), redness (a*24), and yellowness (b*24) than CD. Metabolomic analysis revealed significant differences between the meat extracts from the 2 duck breeds. A total of 49 and 57 significantly different metabolites were identified in positive and negative ion modes, respectively. These differentially accumulated metabolites (DAMs) could be divided into 28 classes, of which the 4 main categories were carbohydrates, amino acids, fatty acids, and eicosanoids. Liancheng white duck might have better nutritional and medicinal value considering the higher content of (4Z,7Z,10Z,13Z,16Z,19Z)-4,7,10,13,16,19-docosahexaenoic acid (DHA), docosapentaenoic acid (DPA), eicosapentaenoic acid (EPA), and prostaglandinF3α (PGF3α), having anti-inflammatory orantioxidant effects. Carbohydrate concentration negatively correlated with pH24. The 4 metabolites positively correlated with the shear force. These results provide an overall perspective for bridging the gap between variation of duck meat quality and metabolites with respect to breed.

Whole-genome sequencing revealed genetic diversity, structure and patterns of selection in Guizhou indigenous chickens

BACKGROUND

The eight phenotypically distinguishable indigenous chicken breeds in Guizhou province of China are great resources for high-quality development of the poultry industry in China. However, their full value and potential have yet to be understood in depth. To illustrate the genetic diversity, the relationship and population structure, and the genetic variation patterns shaped by selection in Guizhou indigenous chickens, we performed a genome-wide analysis of 240 chickens from 8 phenotypically and geographically representative Guizhou chicken breeds and 60 chickens from 2 commercial chicken breeds (one broiler and one layer), together with 10 red jungle fowls (RJF) genomes available from previous studies.

RESULTS

The results obtained in this present study showed that Guizhou chicken breed populations harbored higher genetic diversity as compared to commercial chicken breeds, however unequal polymorphisms were present within Guizhou indigenous chicken breeds. The results from the population structure analysis markedly reflected the breeding history and the geographical distribution of Guizhou indigenous chickens, whereas, some breeds with complex genetic structure were ungrouped into one cluster. In addition, we confirmed mutual introgression within Guizhou indigenous chicken breeds and from commercial chicken breeds. Furthermore, selective sweep analysis revealed candidate genes which were associated with specific and common phenotypic characteristics evolved rapidly after domestication of Guizhou local chicken breeds and economic traits such as egg production performance, growth performance, and body size.

CONCLUSION

Taken together, the results obtained from the comprehensive analysis of the genetic diversity, genetic relationships and population structures in this study showed that Guizhou indigenous chicken breeds harbor great potential for commercial utilization, however effective conservation measures are currently needed. Additionally, the present study drew a genome-wide selection signature draft for eight Guizhou indigenous chicken breeds and two commercial breeds, as well as established a resource that can be exploited in chicken breeding programs to manipulate the genes associated with desired phenotypes. Therefore, this study will provide an essential genetic basis for further research, conservation, and breeding of Guizhou indigenous chickens.

A Genome-Wide Association Study of the Chest Circumference Trait in Xinjiang Donkeys Based on Whole-Genome Sequencing Technology

Animal genotyping by means of genome-wide association studies is important for connecting phenotypes of interest with their underlying genetics in livestock. However, the use of whole genome sequencing to investigate chest circumference (CC) in donkeys has rarely been reported. We aimed to use the genome-wide association study approach to detect significant single nucleotide polymorphisms (SNPs) and key genes associated with chest circumference traits in Xinjiang donkeys. We assessed 112 Xinjiang donkeys in this study. The chest circumference of each was measured 2 h before milking. We re-sequenced blood samples from the Xinjiang donkeys, and genome-wide association study analyses were performed using a mixed model with the PLINK, GEMMA, and REGENIE programs. We tested 38 donkeys for candidate SNPs for genome-wide association study using three software programs. Additionally, 18 SNP markers reached genome-wide significance (p < 1.61 × 10^-9). On the basis of these, 41 genes were identified. Previously proposed candidate genes for CC traits were supported by this study, including NFATC2 (Nuclear Factor of Activated T Cells 2), PROP1 (PROP Paired-Like Homeobox 1), UBB (Ubiquitin B), and HAND2 (Heart and Neural Crest Derivatives Expressed 2). These promising candidates provide a valuable resource for validating potential meat production genes and will facilitate the development of high-yielding Xinjiang donkey breeds through marker-assisted selection or gene editing.

Analysis of genomic copy number variations through whole-genome scan in Chinese Qaidam cattle

Qaidam cattle (CDM) are indigenous breed inhabiting Northwest China. In the present study, we newly sequenced 20 Qaidam cattle to investigate the copy number variants (CNVs) based on the ARS-UMD1.2 reference genome. We generated the CNV region (CNVR) datasets to explore the genomic CNV diversity and population stratification. The other four cattle breeds (Xizang cattle, XZ; Kazakh cattle, HSK; Mongolian cattle, MG; and Yanbian cattle, YB) from the regions of North China embracing 43 genomic sequences were collected and are distinguished from each of the other diverse populations by deletions and duplications. We also observed that the number of duplications was significantly more than deletions in the genome, which may be less harmful to gene formation and function. At the same time, only 1.15% of CNVRs overlapped with the exon region. Population differential CNVRs and functional annotations between the Qaidam cattle population and other cattle breeds revealed the functional genes related to immunity (MUC6), growth (ADAMTSL3), and adaptability (EBF2). Our analysis has provided numerous genomic characteristics of some Chinese cattle breeds, which are valuable as customized biological molecular markers in cattle breeding and production.

Identification of candidate genomic regions for egg yolk moisture content based on a genome-wide association study

BACKGROUND

Eggs represent important sources of protein and are widely loved by consumers. Egg yolk taste is an important index for egg selection, and the moisture content of the egg yolk affects the taste. To understand the molecular mechanism underlying egg yolk moisture content, this study determined the phenotype and heritability of egg yolk water content and conducted a genome-wide association study (GWAS) using a mixed linear model.

RESULTS

We determined the phenotype and heritability of thermogelled egg yolk water content (TWC) and found that the average TWC was 47.73%. Moreover, significant variations occurred (41.06-57.12%), and the heritability was 0.11, which indicates medium-low heritability. Through the GWAS, 48 single nucleotide polymorphisms (SNPs) related to TWC (20 significantly, 28 suggestively) were obtained, and they were mainly located on chromosomes 10 and 13. We identified 36 candidate genes based on gene function and found that they were mainly involved in regulating fat, protein, and water content and embryonic development. FGF9, PIAS1, FEM1B, NOX5, GLCE, VDAC1, IGFBP7, and THOC5 were involved in lipid formation and regulation; AP3S2, GNPDA1, HSPA4, AP1B1, CABP7, EEF1D, SYTL3, PPP2CA, SKP1, and UBE2B were involved in protein folding and hydrolysis; and CSF2, SOWAHA, GDF9, FSTL4, RAPGEF6, PAQR5, and ZMAT5 were related to embryonic development and egg production. Moreover, MICU2, ITGA11, WDR76, BLM, ANPEP, TECRL, EWSR1, and P4HA2 were related to yolk quality, while ITGA11, WDR76, BLM, and ANPEP were potentially significantly involved in egg yolk water content and thus deserve further attention and research. In addition, this study identified a 19.31-19.92 Mb genome region on GGA10, and a linkage disequilibrium analysis identified strong correlations within this region. Thus, GGA10 may represent a candidate region for TWC traits.

CONCLUSION

The molecular genetic mechanism involved in TWC was revealed through heritability measurements and GWAS, which identified a series of SNPs, candidate genes, and candidate regions related to TWC. These results provide insights on the molecular mechanism of egg yolk moisture content and may aid in the development of new egg traits.

Candidate Genes and Gene Networks Change with Age in Japanese Black Cattle by Blood Transcriptome Analysis

Age is an important physiological factor that affects the metabolism and immune function of beef cattle. While there have been many studies using the blood transcriptome to study the effects of age on gene expression, few have been reported on beef cattle. To this end, we used the blood transcriptomes of Japanese black cattle at different ages as the study subjects and screened 1055, 345, and 1058 differential expressed genes (DEGs) in the calf vs. adult, adult vs. old, and calf vs. old comparison groups, respectively. The weighted co-expression network consisted of 1731 genes. Finally, blue, brown, and yellow age-specific modules were obtained, in which genes were enriched in signaling pathways related to growth and development and immune metabolic dysfunction, respectively. Protein-protein interaction (PPI) analysis showed gene interactions in each specific module, and 20 of the highest connectivity genes were chosen as potential hub genes. Finally, we identified 495, 244, and 1007 genes by exon-wide selection signature (EWSS) analysis of different comparison groups. Combining the results of hub genes, we found that VWF, PARVB, PRKCA, and TGFB1I1 could be used as candidate genes for growth and development stages of beef cattle. CORO2B and SDK1 could be used as candidate marker genes associated with aging. In conclusion, by comparing the blood transcriptome of calves, adult cattle, and old cattle, the candidate genes related to immunity and metabolism affected by age were identified, and the gene co-expression network of different age stages was constructed. It provides a data basis for exploring the growth, development, and aging of beef cattle.

Integrated Transcriptomics Profiling in Chahua and Digao Chickens' Breast for Assessment Molecular Mechanism of Meat Quality Traits

Meat quality traits are an important economic trait and remain a major argument, from the producer to the consumer. However, there are a few candidate genes and pathways of chicken meat quality traits that were reported for chicken molecular breeding. The purpose of the present study is to identify the candidate genes and pathways associated with meat quality underlying variations in meat quality. Hence, transcriptome profiles of breast tissue in commercial Digao (DG, 5 male) and Chahua (CH, 5 male) native chicken breeds were analyzed at the age of 100 days. The results found 3525 differentially expressed genes (DEGs) in CH compared to DG with adjusted p-values of ≤0.05 and log2FC ≥ 0.1 FDR ≤ 0.05. Functional analysis of GO showed that the DEGs are mainly involved in the two types of processes of meat quality, such as positive regulation of the metabolic process, extracellular structure organization, collagen trimer, cellular amino acid metabolic process, cellular amino acid catabolic process, and heme binding. Functional analysis of KEGG showed that the DEGs are mainly involved in the two types of processes of meat quality, such as oxidative phosphorylation, carbon metabolism, valine, leucine, and isoleucine degradation, and fatty acid degradation. Many of the DEGs are well known to be related to meat quality, such as COL28A1, COL1A2, MB, HBAD, HBA1, ACACA, ACADL, ACSL1, ATP8A1, CAV1, FADS2, FASN, DCN, CHCHD10, AGXT2, ALDH3A2, and MORN4. Therefore, the current study detected multiple pathways and genes that could be involved in the control of the meat quality traits of chickens. These findings should be used as an essential resource to improve the accuracy of selection for meat traits in chickens using marker-assisted selection based on differentially expressed genes.