An intercross population study reveals genes associated with body size and plumage color in ducks

Identification of SNPs in the second intron of IGF2BP1 and their Association with growth traits in Nanjiang Yellow goat

Insulin-like Growth Factor 2 mRNA-binding Protein 1 (IGF2BP1) is a candidate gene of significant interest for modulating economically important traits in livestock and poultry. The second intron of IGF2BP1 has been implicated in growth-related traits, though its precise mechanistic role remains elusive. Initial resequencing analyses in our laboratory indicated strong selective pressures on the IGF2BP1 genomic region, prompting the selection and identification of several single nucleotide polymorphisms (SNPs). Seven SNPs were mapped to the conserved region of the second intron, necessitating further investigation into their functional relevance and association with growth traits. In this study, 348 Nanjiang Yellow goats were analyzed, and the association analysis via the GLM program in SAS 9.4 identified five SNPs significantly correlated with growth traits. Notably, rs652062749(A > G) emerged as a critical locus influencing later-stage growth traits. Furthermore, strong linkage disequilibrium was observed among three SNPs, with the rs638185407 (T > A) variant markedly enhancing luciferase activity in H293T cells. Combination genotypes TTAACT, TTCCCC, and ATCACT were identified as superior for growth traits, offering theoretical insights for genetic co-breeding. This study underscores the potential utility of IGF2BP1 as a functional genetic marker in Nanjiang Yellow goat breeding programs.

Genomic analyses reveal a lack of widespread strong selection in indigenous chickens

The study of domestication has been revolutionized with the advent of molecular genetics. Chickens, with their clear domestication history, emerge as an excellent model for study into the paths of evolution in domestication and improvement. Here we used genomic data from wild, indigenous, and commercial chickens to better understand how genetic drift and selection translate into their differentiations. Our investigation into the patterns of allelic change and divergence reveals a polygenic architecture governing genetic differentiation during domestication and improvement. We uncover distinctive population-specific differentiations in terms of genes and functions among wild, indigenous, and commercial chickens. Using Runs Of Homozygosity (ROH) based mixed model approach developed in this study, we identified only directional selection signatures occurring in wild and commercial chickens. Notably, our findings suggest that indigenous chickens serve as reservoirs of genetic diversity, necessary for rapid adaptation to new environments or subsequent modern breeding. This work provides unprecedented insights into the chicken domestication and improvement, and it illuminates our understanding of the domestication of other animal species.

Genetic and metabolic mechanisms underlying webbed feet pigmentation in geese: Insights from histological, transcriptomic, and metabolomic analyses

This study systematically investigated the genetic and metabolic mechanisms underlying pigmentation in goose webbed feet by integrating histological, transcriptomic, and metabolomic analyses. Histological examinations revealed significant differences in melanin deposition among webbed feet of varying colors. Dark black webbed feet exhibited the highest melanin content, light black webbed feet showed moderate levels, and colorless webbed feet lacked detectable melanin. Transcriptomic analysis identified substantial variations in the expression levels of key genes involved in melanin biosynthesis, including TYRP1, PMEL, DCT, TYR, OCA2, MC1R, RAB38, WNT16, CAMK2A, and MLANA, between pigmented and colorless webbed feet. Notably, the OCA2 gene exhibited significantly higher expression in dark black webbed feet compared to light black webbed feet, underscoring its pivotal role in regulating pigmentation intensity. Enrichment analysis emphasized the importance of pathways related to tyrosine metabolism, melanin production, and amino acid biosynthesis in determining pigmentation differences. Metabolomic profiling supported these findings, revealing that L-tyrosine and 5,6-dihydroxyindole-2-carboxylic acid are critical metabolites in the melanin biosynthesis pathway. Specifically, elevated levels of L-tyrosine were detected in colorless webbed feet, likely due to inhibited melanin synthesis, whereas 5,6-dihydroxyindole-2-carboxylic acid levels were highest in dark black webbed feet, reflecting active melanin production. Correlation analysis between transcriptomic and metabolomic data further validated the central role of tyrosine metabolism and melanin biosynthesis pathways in pigmentation. In conclusion, this study employed multi-omics approaches to elucidate the critical role of the OCA2-centered genetic-metabolic regulatory network in melanin deposition of goose webbed feet, providing important insights into the molecular mechanisms of avian pigmentation and valuable references for poultry breeding.

Nine high-quality Anas genomes provide insights into Anas evolution and domestication

Evolutionary studies of wild and domestic organisms have yielded fascinating discoveries, while the species diversity and the domestication of ducks remain unclear. Here, we assembled eight chromosome-level Anas genomes, combined with the Pekin duck genome, to investigate Anas evolution and domestication. We found that, compared to autosomes, the Z chromosome was less affected by introgression and exhibited relatively stable local phylogenies. From the Z chromosome perspective, we proposed that the speciation of Anas platyrhynchos and Anas zonorhyncha was accompanied by continuous female-biased gene flow and remodeled duck domestication history. Moreover, we constructed an Anas pan-genome and identified several differentiated structural variations (SVs) between domestic and wild ducks. These SVs likely regulate their neighboring genes (i.e., GHR and FER), which represented the promising "domestication genes." Furthermore, a long terminal repeat (LTR) retrotransposon burst was found to have accelerated duck domestication, specifically contributing to functional shifts of the notable MITF and IGF2BP1 genes. These findings presented a live example for understanding animal evolutionary processes.

Metabolome genome-wide association analyses identify a splice mutation in AADAT affects lysine degradation in duck skeletal muscle

Metabolites in skeletal muscles play an important role in their growth, development, immunity and other physiological activities. However, the genetic basis of metabolites in skeletal muscle remains poorly understood. Here, we identified 247 candidate divergent regions containing 905 protein-coding genes closely related to metabolic pathways, including lysine degradation and fatty acid biosynthesis. We then profiled 3,060 metabolites in 246 skeletal muscle samples from F2 segregating population generated by mallard×Pekin duck crosses using metabolomic approaches. We identified 2,044 significant metabolome-based GWAS signals and 21 candidate genes potentially modulating metabolite contents in skeletal muscle. Among them, the levels of 2-aminoadipic acid in skeletal muscle were significantly correlated with body weight and intramuscular fat content, determined by a 939-bp CR1 LINE insertion in AADAT. We further found that the CR1 LINE insertion most possibly led to a splice mutation in AADAT, resulting in the downregulation of the lysine degradation pathway in skeletal muscle. Moreover, intramuscular fat content and fatty acids biosynthesis pathway was significantly increased in individuals with CR1 LINE insertion. This study enhances our understanding of the genetic basis of skeletal muscle metabolic traits and promotes the efficient utilization of metabolite traits in the genetic improvement of animals.

RNA sequencing analysis reveals key genes and pathways associated with feather pigmentation in mule ducks

Feather color is an important morphological trait of poultry. At present, the reports on the inheritance of plumage color of mule ducks at the molecular level are few, and the regulatory mechanism in white plumage rates of different mule ducks remains unclear. This study aimed to broaden the understanding of the white plumage rates in mule ducks to improve their production value. We used RNA sequencing to analyze and compare the mRNA expression profiles in hair follicle tissues from 10-week-old mule ducks with black and white plumages, thereby revealing the temporal gene expression patterns and pathways associated with plumage color regulation. In total, 1672 annotated differentially expressed genes (DEGs) were identified in black and white plumages from different databases between mule ducks with the 2 plumage colors. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment in hair follicle tissues indicated that the aforementioned DEGs were mainly involved in the melanin signaling pathway. Concurrently, we use weighted gene co-expression network analysis to detect core modules and hub genes associated with melanin biosynthesis in feathers. The green module exhibited the strongest correlation with the phenotypic traits, encompassing a total of 1049 genes. Subsequent Kyoto Encyclopedia of Genes and Genomes enrichment analysis identified 11 genes as pivotal in the melanin biosynthetic pathway, including EDRNB2, TYR, KIT, EDNRB, and MC1R. The differential expression of eight selected DEGs was verified using quantitative reverse transcription-PCR, and the results were consistent with RNA-seq data. This study provides a basis for understanding the differences in plumage color development in mule ducks.

Mining of candidate genes related to body size in Chinese native pig breeds based on public data

To study the key genes that influence the body size of local pig breeds in China. Genome-wide SNP chip data from a total of 129 pigs from eight breeds, consisting of four large body size breeds (MZ, HT, ST, RC) and four small body size breeds (XI, BX, WZ, DN) were analyzed. Principal Component Analysis (PCA) was employed to assess the genetic clustering of the eight breeds. Fst and XP-CLR were used to detect selective signals between the large ans small body size breeds groups. The PCA results indicated a clear clustering of small breeds and a dispersion distribution among large breeds. Fst and XP-CLR identified 142 overlapping regions within a 500 kb up & down stream of significant loci. These regions encompassed 520 annotated genes, which were enriched in 34 biological pathways. Gene network analysis highlighted nine key genes, of which five (NPR3, TNFSF11, TBC1D7, FGF2, IGF1R) are known to be associated with bone growth and body size traits in animals. Additionally, four novel candidate genes (IKBKB, SFRP1, LRP6, SPRY1) were identified that might be related to pig body size. Our findings provide a theoretical basis for further revealing the genetic mechanism of pig body size traits.

An mRNA expression atlas for the duck with public RNA-seq datasets

BACKGROUND

Ducks are globally important poultry species and a major source of farm animal products, including meat, eggs, and feathers. A thorough understanding of the functional genomic and transcriptomic sequences is crucial for improving production efficiency.

RESULT

This study constructed the largest duck mRNA expression atlas among all waterfowl species to date. The atlas encompasses 1,257 tissue samples across 30 tissue types, representing all major organ systems. Using advanced clustering analysis, we established co-expression network clusters to describe the transcriptional features in the duck mRNA expression atlas and, when feasible, assign these features to unique tissue types or pathways. Additionally, we identified 27 low-variance, highly expressed housekeeping genes suitable for gene expression experiments. Furthermore, in-depth analysis revealed potential sex-biased gene expression patterns within tissues and specific gene expression profiles in meat-type and egg-type ducks, providing valuable resources to understand the genetic basis of sex differences and particular phenotypes. This research elucidates the biological processes affecting duck productivity.

CONCLUSION

This study presents the most extensive gene expression atlas for any waterfowl species to date. These findings are of significant value for advancing duck biological research and industrial applications.

Standardized ileal digestible lysine requirements based on growth performance of White Pekin ducks for 21 days after hatch

The purpose of this study was to assess the optimal standardized ileal digestible (SID) lysine (Lys) requirement for male White Pekin ducklings with a specific focus on growth performance for the 3 weeks following hatching. A total of 384 one-day-old male White Pekin ducklings were allocated to six different dietary treatments, each containing varying levels of digestible Lys content ranging from 0.72% to 1.12%. All amino acids in the diets remained consistent except for Lys. The ducklings were randomly distributed into 24-floor pens, with each treatment group comprising eight pens, and each pen housing eight ducklings. The diets were offered ad-libitum throughout the study. Weekly measurements of body weight and feed intake were recorded to calculate the feed conversion ratio. The SID Lys requirement was determined by analyzing the data using both linear-plateau and quadratic-plateau models and calculating the mean value. The results demonstrated a significant linear (p < 0.001) and quadratic (p < 0.001) improvement in body weight gain and feed efficiency with increasing SID Lys content in the diet. According to the linear-plateau regression analysis, the estimated SID Lys requirements for final body weight, weight gain, and feed efficiency were 1.00%, 1.00%, and 0.98%, respectively. Conversely, the quadratic-plateau regression analysis yielded estimated SID Lys requirements of 1.11%, 1.11%, and 1.10%, respectively, for the same parameters. In summary, this study established that the recommended SID Lys levels for White Pekin ducklings for the 3 wk period after hatching were found to be 1.05%, 1.05%, and 1.04% for achieving the finest final body weight, daily gain, and feed efficiency, respectively.

The Impact of Mutant EDNRB on the Two-End Black Coat Color Phenotype in Chinese Local Pigs

Endothelin Receptor Type B (EDNRB) is expressed in a variety of cells during embryonic stage, including melanocyte precursors cells. Our previous studies found that 11 bp deletion of EDNRB caused the two-end black (TEB) coat color in Chinese pigs. In this study, we aimed to explore the mutant EDNRB on the formation of TEB coat color in Chinese pigs. We constructed recombinant plasmid for wild and mutant EDNRB and EDN1, respectively, and transfected the recombinant plasmid into mouse B16 melanoma cells in groups. Real-time fluorescent quantitative PCR (RT-qPCR) was performed to detect expression of genes that participate in melanin pathway, including PLCγ, Raf, MITF. Comparing to the wild-type EDNRB cells, expression of the three genes in the cell line expressing mutant EDNRB cells was significantly reduced. We measured the melanin content produced by transfected recombinant granulocytes of wild and mutant EDNRB and found that the amount of melanin in mutant EDNRB cells was significantly lower than that of the wild. Wound-healing assay confirmed that the migration and mobility rate of mutant EDNRB cells were significantly lower than the wild. Co-immunoprecipitation further confirmed that mutant EDNRB could not interact with the EDN1 protein. In conclusion, this study revealed that the 11 bp deletion of EDNRB reduced the melanin production, which may be caused by inhibiting the expression of PLCγ, Raf, and MITF. The mutant EDNRB reduced melanocyte migration and could not interact with the EDN1 protein. We explored the effect of mutant EDNRB in Chinese pigs with TEB coat color, and the results provided a reference for exploring molecular mechanism of mutant EDNRB on the formation of TEB coat color pigs.

Genomic variation responding to artificial selection on different lines of Pekin duck

Understanding the genomic variation in Pekin duck under artificial selection is important for improving the utilization of duck genetic resources. Here, the genomic changes in Pekin duck were analyzed by using the genome resequencing data from 96 individual samples, including 2 conservation populations and 4 breeding populations with different breeding backgrounds. The population structure, runs of homozygosity (ROH), effective population number (Ne), and other genetic parameters were analyzed. The breeding populations showed lower genetic diversity compared to the conservation populations. Maple Leaf duck and Cherry Valley duck retained low genetic diversity compared to other breeding populations, with Cherry Valley duck showing the lowest diversity and the highest inbreeding coefficient. This suggested that Cherry Valley and Maple Leaf ducks have undergone intensive selection compared to other breeding populations. By the analysis of runs of homozygosity (ROHs), some genes (e.g., IGF1R) associated with growth traits were identified. By the analysis of the selection signal, strong selection characteristics in certain genomic regions during the breeding of Peking duck across different selected lines were observed. In addition, copy number variations (CNVs) in Pekin duck populations were analyzed. Six regions of interest were identified, containing RPA1, DOT1L, SLC25A42, RALYL, TRPA1, and IGFBP2. Furthermore, the allele frequency distribution of these genes showed significant differences between breeding populations and conservation populations, indicating that these candidate genes could have undergone strong selection pressure during long-term selection for improved production. These findings contribute to a deeper understanding of the distinct evolutionary processes in Pekin ducks under artificial selection and provide valuable insights for future breeding strategies.

Identification of quantitative trait loci and candidate genes associated with growth curve parameters in chinese wenshang barred chickens

The growth curve is a vital instrument for assessing and forecasting weight and developmental shifts in livestock and poultry, which reflects the changes of bodyweight traits with time and plays a key role in guiding breeding and production approaches. This study performed a genome-wide association study (GWAS) for growth curve parameters generated by nonlinear models which fit original weight-age records, to discover the SNPs and candidate genes correlated with growth traits. Data from 362 Chinese Wenshang Barred Chickens weighed at the age of 0, 2, 4, 6, 8, 10, 12, 14, 16, and 18 weeks were used to fit the Gompertz, Logistic, and von Bertalanffy growth models. The Gompertz model showed the highest coefficient of determination (R2 = 0.974). The mature body weight (A), time scale (b), and maturity rate (k) were treated as phenotypes for single-trait GWAS. The GWAS identified 44, 103, and 5 significant single nucleotide polymorphisms (SNPs) associated with A, b, and K, respectively. Among them, several candidate genes, including LDB2, TOB2, RCBTB1, KPNA3, SLIT2, LCORL, LAP3, and TPRA1, were previously reported to be associated with growth and development. Two lead SNPs (4:76022389, 4:76070237) on the LDB2 gene were significantly associated with the growth curve. Further research of these candidate genes could help explore the full genetic architecture underlying growth and development traits in poultry.

A high-quality assembly revealing the PMEL gene for the unique plumage phenotype in Liancheng ducks

BACKGROUND

Plumage coloration is a distinctive trait in ducks, and the Liancheng duck, characterized by its white plumage and black beak and webbed feet, serves as an excellent subject for such studies. However, academic comprehension of the genetic mechanisms underlying duck plumage coloration remains limited. To this end, the Liancheng duck genome (GCA_039998735.1) was hereby de novo assembled using HiFi reads, and F2 segregating populations were generated from Liancheng and Pekin ducks. The aim was to identify the genetic mechanism of white plumage in Liancheng ducks.

RESULTS

In this study, 1.29 Gb Liancheng duck genome was de novo assembled, involving a contig N50 of 12.17 Mb and a scaffold N50 of 83.98 Mb. Beyond the epistatic effect of the MITF gene, genome-wide association study analysis pinpointed a 0.8-Mb genomic region encompassing the PMEL gene. This gene encoded a protein specific to pigment cells and was essential for the formation of fibrillar sheets within melanosomes, the organelles responsible for pigmentation. Additionally, linkage disequilibrium analysis revealed 2 candidate single-nucleotide polymorphisms (Chr33: 5,303,994A>G; 5,303,997A>G) that might alter PMEL transcription, potentially influencing plumage coloration in Liancheng ducks.

CONCLUSIONS

Our study has assembled a high-quality genome for the Liancheng duck and has presented compelling evidence that the white plumage characteristic of this breed is attributable to the PMEL gene. Overall, these findings offer significant insights and direction for future studies and breeding programs aimed at understanding and manipulating avian plumage coloration.

Genomic evidence for hybridization and introgression between blue peafowl and endangered green peafowl and molecular foundation of leucistic plumage of blue peafowl

INTRODUCTION

The blue peafowl (Pavo cristatus) and the green peafowl (Pavo muticus) have garnered significant public affection due to their stunning appearance, although the green peafowl is currently endangered. The causative mutation that causes the leucistic plumage of the blue peafowl (also called white peafowl) remains unknown.

RESULTS

In this study, we generated a chromosome-level reference genome of the blue peafowl with a contig N50 of 30.6 Mb, including the autosomes, Z and W sex chromosomes, and a complete mitochondria DNA sequence. Data from 77 peafowl whole genomes, 76 peafowl mitochondrial genomes, and 33 peafowl W chromosomes genomes provided the first substantial genetic evidence for recent hybridization between green peafowls and blue peafowls. We found 3 hybrid green peafowls in zoo samples rather than in the wild samples, with a blue peafowl genomic content of 16-34%. Maternal genetic analysis showed that 2 of the hybrid female green peafowls contained complete blue peafowl mitochondrial genomes and W chromosomes. Some animal protection agencies release captive green peafowls in order to maintain the wild population of green peafowls. Therefore, to better protect the endangered green peafowl, we suggest that purebred identification must be carried out before releasing green peafowls from zoos into the wild in order to prevent the hybrid green peafowl from contaminating the wild green peafowl. In addition, we also found that there were historical introgression events of green peafowl to blue peafowl in 4 zoo blue peafowl individuals. The introgressed genomic regions contain IGFBP1 and IGFBP3 genes that could affect blue peafowl body size. Finally, we identified that the nonsense mutation (g.4:12583552G>A) in the EDNRB2 gene is the genetic causative mutation for leucistic plumage of blue peafowl, preventing melanocytes from being transported into plumage, thereby inhibiting melanin deposition.

CONCLUSION

Our research provides both theoretical and empirical support for the conservation of the endangered green peafowl. The high-quality genome and genomic data also provide a valuable resource for blue peafowl genomics-assisted breeding.

Genomic insights into the specialisation and selection of the Jinding duck

The domestication of ducks represents a pivotal evolutionary shift from the unguided propagation of wild species to deliberate human-mediated selection, culminating in distinct behavioural, morphological, and physiological traits that differentiate domesticated ducks from their wild counterparts. This transition has yielded breeds with traits fine-tuned to specific economic roles, such as egg production, meat yield, or dual-purpose functionality. Duck domestication plays a significant role in poultry production globally, meeting the growing demand for eggs and meat in various regions. Here, we focus on the Jinding Duck (JDD), a breed renowned for its prolific egg-laying traits. Employing whole-genome resequencing data from 325 individuals across five Chinese indigenous duck breeds, we aimed to dissect the unique population structure and assess the genetic diversity within the JDD cohort. The findings reveal the distinct genetic heritage of JDD, diverged from other domesticated breeds, and show a relative paucity of genetic diversity. A salient discovery was a 200 kb genomic interval containing three genes (NCF2, SMG7, and ARPC5) with almost exclusive haplotypes, which were inherited from Anas platyrhynchos or Anas zonorhyncha, impacting the morphological attributes of JDD. The study highlights a c.28G>A non-synonymous mutation in the first exon of the LAMC1 gene, which is potentially influencing feather morphology in JDD. Our findings suggest that unique blue eggshell colouration in JDD is likely attributable to variations within the promoter element of the ABCG2 gene, distinguishing it from other breeds. Moreover, the MAP7 and FHL1 genes emerge as significant factors in the laying performance of JDD. These genetic insights are not only crucial for improving the JDD breed but also provide valuable information that could be applied to duck breeding programmes worldwide, helping enhance productivity and meet international demands for high-efficiency poultry breeds.

Insight Into Body Size Evolution in Aves: Based on Some Body Size-Related Genes

Birds exhibit remarkable variations in body size, making them an ideal group for the study of adaptive evolution. However, the genetic mechanisms underlying body size evolution in avian species remain inadequately understood. This study investigates the evolutionary patterns of avian body size by analyzing 15 body-size-related genes, including GHSR, IGF2BP1, and IGFBP7 from the growth hormone/insulin-like growth factor axis, EIF2AK3, GALNS, NCAPG, PLOD1, and PLAG1 associated with tall stature, and ACAN, OBSL1, and GRB10 associated with short stature, four genes previously reported in avian species: ATP11A, PLXDC2, TNS3, and TUBGCP3. The results indicate significant adaptive evolution of body size-related genes across different avian lineages. Notably, in the IGF2BP1 gene, a significant positive correlation was observed between the evolutionary rate and body size, suggesting that larger bird species exhibit higher evolutionary rates of the IGF2BP1 gene. Furthermore, the IGFBP7 and PLXDC2 genes demonstrated accelerated evolution in large- and medium-sized birds, respectively, indicating distinct evolutionary patterns for these genes among birds of different sizes. The branch-site model analysis identified numerous positively selected sites, primarily concentrated near functional domains, thereby reinforcing the critical role of these genes in body size evolution. Interestingly, extensive convergent evolution was detected in lineages with larger body sizes. This study elucidates the genetic basis of avian body size evolution for the first time, identifying adaptive evolutionary patterns of body size-related genes across birds of varying sizes and documenting patterns of convergent evolution. These findings provide essential genetic data and novel insights into the adaptive evolution of body size in birds.

A novel codominant plumage color pattern of white breast patches in WugangTong geese was controlled by EDNRB2

Two basic plumage color patterns are observed in adult geese: solid grey (G) or colorless white (W). However, a Chinese indigenous breed, the Wugangtong goose (WGT), continues to be subject to selective breeding efforts as it displays segregation of plumage colors, including G, W, and a novel color pattern designated Wb (G with white breast circles). The underlying genetic mechanisms responsible for the Wb phenotype are yet to be determined. The current study employed the population differentiation index (FST) to analyze 90 geese exhibiting diverse plumage colors, identifying the fifth intron of EDNRB2 as a particularly noteworthy region with the highest FST values. Sanger sequencing of the region surrounding the EDNRB2 gene identified a 14-bp insertion within exon 3 as the causal mutation. The heterozygosity of this 14-bp insertion and wild-type alleles was completely associated with the Wb phenotype, thereby substantiating the codominant nature of the G and W phenotypes. An inter-species corroborated this finding cross between the graylag (no 14-bp insertion) and the swan goose (homozygous for the 14-bp insertion) breeds, as hybrids from this cross exhibited the Wb phenotype. Transcriptomes from white breast patches and gray dorsal skins of 4 Wb geese were compared. A significant downregulation of genes involved in melanin synthesis and melanocyte development was observed, including EDRNB2 and MLANA. The downregulation of MLANA indicated that the mutated EDNRB2 resulted in melanocyte loss in specific body regions, as MLANA is a marker gene for melanocytes. The findings were corroborated by melanin staining using the Mansson-Fontana method, which revealed no melanin particles deposited in the white breast patches. In summary, the gray plumage color was codominant to the white color in WGT geese, and plumage color variations were controlled by EDNRB2. The findings of our study offer valuable and practical guidance for the purification of plumage colors among WGT, whether through traditional phenotype selection or molecular breeding methods.

Identification of SNPs and INDELS associated with duck egg quality traits through a genome-wide association analysis

Egg quality traits are economically important in the poultry industry. To explore the genetic architecture and identify potential candidate genes, a genome-wide association study (GWAS) was performed for 13 egg quality traits using data from whole-genome sequencing of 299 Longyan Shan-ma female ducks, including 12 quantitative traits and one qualitative trait, eggshell color (ESC; white, light green, green). From estimation of pedigree genetic parameters, heritability (h2) ranged from 0.022 to 0.996 for the 12 quantitative traits, with the highest h2 (0.996) for eggshell color a* value (ESCA) and the lowest h2 (0.022) for egg yolk percentage relative to EW. A total of 8,874 single nucleotide polymorphism (SNP)-based significant associations (1.0 × 10-6) and 247 insertion-deletion (indel)-based significant associations (1.00 × 10-5) were identified, including 5,980 SNPs and 159 indel markers. From 5,924 SNPs and 143 indels associated with ESC traits, 181 potential candidate genes were identified, and most significant SNPs and indels (P < 1.0 × 10-20) were located at 1.86 Mb (44.29-46.15 Mb) on chromosome 4. The top SNP (chr4:45325309:C>A; P = 7.97 × 10-43) and the top indel (chr4:45299595:delTTCCACTCCAC; P = 4.20 × 10-36) for the ESC a* value were within two known ESC candidate genes; ATP-binding cassette subfamily G member 2 (ABCG2) and protein kinase cGMP-dependent 2 (PRKG2). Of 56 SNPs and 16 indels associated with other egg quality traits, 46 potential candidate genes were identified including synapse differentiation-inducing 1-like (SYNDIG1L) for EW, and core histone macro-H2A.1 (LOC101795967) and neurogenin 1 (NEUROG1) for egg shape index; and four genes including collagen type VI alpha 3 chain (COL6A3), lysine demethylase 7A (KDM7A), LOC101802169, and sperm-associated antigen 16 (SPAG16) for egg yolk weight and the percentage of yolk to total egg weight. Of the 46 genes, the molecular functions of 22 are related to protein binding, indicating important roles in the formation of egg quality traits. Our findings provide new insight into the genetic basis of egg quality traits in ducks.

Genome-wide association studies reveal the genetic basis of growth and carcass traits in Sichuan Shelduck

China has abundant local duck resource populations, and evaluating the characteristics of these breeds will help improve development and utilization. In this study, we conducted the first investigations of growth and slaughter performance on Sichuan Shelduck (n = 240), an endangered duck local breed. The average body weight is 1497.91 g at 90 d of age. According to the growth curve through data recorded every 2 wk, we observed a low relative growth rate (RGR) for the early growth stage. The RGR shows a decreasing trend with age increasing in the stage from 0 to 56 d of age. The SNP-based heritability estimation showed the growth rate has a relatively high heritability, indicating high genetic stability for this trait. In the correlation analysis, the percentage of leg muscle is positively correlated with the absolute growth rate (AGR) at 28 to 42 d of age, whereas it is negatively correlated with the earlier stages, exhibiting a time-specific correlation result. Additionally, genome-wide association studies (GWAS) identified PCSK6, TOX2, and TOMM7 as potential candidate genes influencing AGR (42-56) and AGR (56-90), while the candidate genes of slaughter traits were PTP4A2, FAM110B, TOX, UBXN2B, and FCHSD2. These results provide an important reference for further understanding the genetic basis of growth and meat production performance of Sichuan Shelduck.

Joint Genomic and Transcriptomic Analysis Reveals Candidate Genes Associated with Plumage Color Traits in Matahu Ducks

Plumage color is a key trait for identifying waterfowl breeds with significant economic importance. A white-feathered group has recently emerged within the native Matahu duck population, presenting an opportunity for breeding new lines. However, the genetic basis for this plumage variation is still unknown, necessitating further research. This study aims to identify the genetic mechanisms underlying the emergence of white-feathered individuals in the Matahu duck population through combined genome and transcriptome analysis, providing insights for selective breeding and the development of new white-feathered lines. In this study, a total of 1344 selected genes and 1406 significantly differentially expressed genes were identified through selection signal analysis and transcriptomic analysis, respectively. The functional enrichment of these genes revealed several key signaling pathways, including those related to cGMP-PKG, cAMP, PI3K-Akt, and MAPK. Furthermore, important candidate genes involved in melanin biosynthesis, such as MITF, MC1R, TYR, TYRP1, and ABCB6, were identified. Notably, 107 genes were detected by both methods, and, among these, DGKI, GPRC5B, HMX1, STS, ADGRA1, PRKAR2B, and HOXB9 are suggested to play a role in melanin formation and potentially influence plumage traits. Through the integrative approach combining genomic selection signals and transcriptomic analyses, we identified several candidate genes directly associated with plumage color, including MITF, TYR, TYRP1, and MC1R, along with multiple signaling pathways linked to melanin formation. We hypothesize that the expression of DGKI, GPRC5B, HMX1, STS, ADGRA1, PRKAR2B, and HOXB9, detected by both methods, may be closely related to the regulation of plumage color traits. These findings provide a foundational basis for further research aimed at elucidating the genetic mechanisms governing plumage color variation in ducks.

Integrative 3D genomics with multi-omics analysis and functional validation of genetic regulatory mechanisms of abdominal fat deposition in chickens

Chickens are the most abundant agricultural animals globally, with controlling abdominal fat deposition being a key objective in poultry breeding. While GWAS can identify genetic variants associated with abdominal fat deposition, the precise roles and mechanisms of these variants remain largely unclear. Here, we use male chickens from two lines divergently selected for abdominal fat deposition as experimental models. Through the integration of genomic, epigenomic, 3D genomic, and transcriptomic data, we build a comprehensive chromatin 3D regulatory network map to identify the genetic regulatory mechanisms that influence abdominal fat deposition in chickens. Notably, we find that the rs734209466 variant functions as an allele-specific enhancer, remotely enhancing the transcription of IGFBP2 and IGFBP5 by the binding transcription factor IRF4. This interaction influences the differentiation and proliferation of preadipocytes, which ultimately affects phenotype. This work presents a detailed genetic regulatory map for chicken abdominal fat deposition, offering molecular targets for selective breeding.

Genome-wide association study reveals 2 copy number variations associated with the variation of plumage color in the white duck hybrid population

Plumage color is an intuitive external poultry characteristic with rich manifestations and complex genetic mechanisms. In our previous study, we observed that there were more dark variations in plumage color in the F2 population derived from the hybridization of 2 white duck varieties. Therefore, based on the statistics of plumage color of 308 F2 populations, we further used the resequencing data of these individuals to detect copy number variations (CNVs) in the whole genome and conducted genome-wide association studies (GWAS) to determine the genetic basis related to plumage color traits. The CNV detection revealed 9,337 CNVs, with an average length of 15,950 bp and a total length of 142.02 MB, accounting for approximately 12.91% of the reference genome. The CNV distribution on the chromosomes was relatively uniform, and the number of CNVs on each chromosome positively correlated with the length of the chromosome. In the pure black plumage group, 2,101 CNVs were only identified, and 1,714 were specifically identified in the pure white plumage group. Ten CNVs were randomly selected for validation using quantitative real-time PCR, and 9 CNVs had the same CNV types as predicted, with an accuracy of 90%. Based on GWAS, we identified 2 CNVs potentially associated with plumage color variations, with the associated CNV regions covering 9 genes. Enrichment analysis of these 9 candidate genes showed significant enrichment of 3 pathways (ribosome biogenesis in eukaryotes, RNA transport, and protein export) and 17 gene ontology terms. Among these, VWA5A can downregulate MITF by binding to the regulatory factors SOX10. The occurrence of CNV may indirectly contribute to duck plumage color variation by affecting the regulatory factors of the switch gene MITF in the melanogenesis pathway. These findings have improved the understanding of the genetic basis of duck plumage color variation and have been beneficial for developing and using plumage color traits in subsequent poultry breeding.

Genome-Wide Analysis of Genetic Diversity and Selection Signatures in Zaobei Beef Cattle

This investigation provides a comprehensive analysis of genomic diversity and selection signatures in Zaobei beef cattle, an indigenous breed known for its adaptation to hot and humid climates and superior meat quality. Whole-genome resequencing was conducted on 23 Zaobei cattle, compared with 46 Simmental cattle to highlight genetic distinctions. Population structure analysis confirmed the genetic uniqueness of Zaobei cattle. Using methods such as DASDC v1.01, XPEHH, and θπ ratio, we identified 230, 232, and 221 genes through DASDC, including hard sweeps, soft sweeps, and linkage sweeps, respectively. Coincidentally, 109 genes were identified when using XPEHH and θπ ratio methods. Together, these analyses revealed eight positive selection genes (ARHGAP15, ZNF618, USH2A, PDZRN4, SPATA6, ROR2, KCNIP3, and VWA3B), which are linked to critical traits such as heat stress adaptation, fertility, and meat quality. Moreover, functional enrichment analyses showed pathways related to autophagy, immune response, energy metabolism, and muscle development. The comprehensive genomic insights gained from this study provide valuable knowledge for breeding programs aimed at enhancing the beneficial traits in Zaobei cattle.

Genome-wide variation study and inter-tissue communication analysis unveil regulatory mechanisms of egg-laying performance in chickens

Egg-laying performance is of great economic importance in poultry, but the underlying genetic mechanisms are still elusive. In this work, we conduct a multi-omics and multi-tissue integrative study in hens with distinct egg production, to detect the hub candidate genes and construct hub molecular networks contributing to egg-laying phenotypic differences. We identifiy three hub candidate genes as egg-laying facilitators: TFPI2, which promotes the GnRH secretion in hypothalamic neuron cells; CAMK2D, which promotes the FSHβ and LHβ secretion in pituitary cells; and OSTN, which promotes granulosa cell proliferation and the synthesis of sex steroid hormones. We reveal key endocrine factors involving egg production by inter-tissue crosstalk analysis, and demonstrate that both a hepatokine, APOA4, and an adipokine, ANGPTL2, could increase egg production by inter-tissue communication with hypothalamic-pituitary-ovarian axis. Together, These results reveal the molecular mechanisms of multi-tissue coordinative regulation of chicken egg-laying performance and provide key insights to avian reproductive regulation.

Genome-Wide Association Analysis Identifies Important Haplotypes and Candidate Gene XKR4 for Body Size Traits in Pekin Ducks

Body size is an important growth indicator in ducks and is a primary selection criterion for physical improvement. An excessively rapid growth rate in meat ducks can result in excessive body size, which may hinder subsequent processing and slaughter operations. However, only a few molecular markers related to body size have been studied in meat ducks. In this study, we performed a genome-wide association study (GWAS) to identify candidate genes and QTLs affecting body length (BL), keel bone length (KBL), neck length (NL), and breast width (BrW) in Pekin ducks (Anas platyrhynchos domestica). Our results indicate the significant SNP for NL is located within a pseudogene, whereas the significant SNP for BrW is located in an intergenic region. More importantly, our analysis identified a haplotype that was significantly associated with both BL and KBL. This haplotype, containing 48 single-nucleotide polymorphisms (SNPs), is localized within the XKR4 gene. The identification of this haplotype suggests that XKR4 may be a key candidate gene influencing BL and KBL in Pekin ducks. These findings have important implications for the breeding and genetic improvement of Pekin ducks, and provide valuable insights into the genetic architecture of body size traits in this species.

Duck gut metagenome reveals the microbiome signatures linked to intestinal regional, temporal development, and rearing condition

The duck gastrointestinal tract (GIT) harbors an abundance of microorganisms that play an important role in duck health and production. Here, we constructed the first relatively comprehensive duck gut microbial gene catalog (24 million genes) and 4437 metagenome-assembled genomes using 375 GIT metagenomic samples from four different duck breeds across five intestinal segments under two distinct rearing conditions. We further characterized the intestinal region-specific microbial taxonomy and their assigned functions, as well as the temporal development and maturation of the duck gut microbiome. Our metagenomic analysis revealed the similarity within the microbiota of the foregut and hindgut compartments, but distinctive taxonomic and functional differences between distinct intestinal segments. In addition, we found a significant shift in the microbiota composition of newly hatched ducks (3 days), followed by increased diversity and enhanced stability across growth stages (14, 42, and 70 days), indicating that the intestinal microbiota develops into a relatively mature and stable community as the host duck matures. Comparing the impact of different rearing conditions (with and without water) on duck cecal microbiota communities and functions, we found that the bacterial capacity for lipopolysaccharide biosynthesis was significantly increased in ducks that had free access to water, leading to the accumulation of pathogenic bacteria and antibiotic-resistance genes. Taken together, our findings expand the understanding of the microbiome signatures linked to intestinal regional, temporal development, and rearing conditions in ducks, which highlight the significant impact of microbiota on poultry health and production.

Comparative genome-wide association study on body weight in Chinese native ducks using four models

The Jinling White duck represents a newly developed breed characterized by a rapid growth rate and a superior meat quality, offering significant economic value and research potential; however, the genetic basis underlying their body weight traits remains less understood. Here, we performed whole-genome resequencing for 201 diverse Jinling White male ducks and conducted population genomic analyses, suggesting a rich genetic diversity within the Jinling White duck population. Equipped with our genomic resources, we applied genome-wide association analysis for body weight on birth (BWB), body weight on 1 wk (BW1), body weight on 3 wk (BW3), body weight on 5 wk (BW5) and body weight on 7 wk (BW7) using 4 statistical models. Comparative studies indicated that factored spectrally transformed linear mixed models (FaST-LMM) demonstrated the most superior efficiency, yielding more results with the minimal false positives. We discovered that PUS7, FBXO11, FOXN2, MSH6, and SLC4A4 were associated with BWB. RAG2, and TMEFF2 were candidate genes for BW1, and STARD13, Klotho, ZAR1L are likely candidates for BW3 and BW5. PLXNC1, ATP1A1, CD58, FRYL, OCIAD1, and OCIAD2 were linked to BW7. These findings provide a genetic reference for the selection and breeding of Jinling White ducks, while also deepened our understanding of Growth and development phenotypic in ducks.

A comparative analysis reveals the genomic diversity among 8 Muscovy duck populations

The Muscovy duck (Cairina moschata) is a waterfowl indigenous to the neotropical regions of Central and South America. It has low demand for concentrated feed and strong adaptability to different rearing conditions. After introduced to China through Eurasian commercial trade, Muscovy ducks have a domestication history of around 300 years in the Fujian Province of China. In the 1990s, the commodity Muscovy duck breed "Crimo," cultivated in Europe, entered the Chinese market for consumption and breeding purposes. Due to the different selective breeding processes, Muscovy ducks have various populational traits and lack transparency of their genetic background. To remove this burden in the Muscovy duck breeding process, we analyzed genomic data from 8 populations totaling 83 individuals. We identify 11.24 million single nucleotide polymorphisms (SNPs) and categorized these individuals into the Fujian-bred and the Crimo populations according to phylogenetic analyses. We then delved deeper into their evolutionary relationships through assessing population structure, calculating fixation index (FST) values, and measuring genetic distances. Our exploration of runs of homozygosity (ROHs) and homozygous-by-descent (HBD) uncovered genomic regions enriched for genes implicated in fatty acid metabolism, development, and immunity pathways. Selective sweep analyses further indicated strong selective pressures exerted on genes including TECR, STAT2, and TRAF5. These findings provide insights into genetic variations of Muscovy ducks, thus offering valuable information regarding genetic diversity, population conservation, and genome associated with the breeding of Muscovy ducks.

Multi-omic characterization of allele-specific regulatory variation in hybrid pigs

Hybrid mapping is a powerful approach to efficiently identify and characterize genes regulated through mechanisms in cis. In this study, using reciprocal crosses of the phenotypically divergent Duroc and Lulai pig breeds, we perform a comprehensive multi-omic characterization of regulatory variation across the brain, liver, muscle, and placenta through four developmental stages. We produce one of the largest multi-omic datasets in pigs to date, including 16 whole genome sequenced individuals, as well as 48 whole genome bisulfite sequencing, 168 ATAC-Seq and 168 RNA-Seq samples. We develop a read count-based method to reliably assess allele-specific methylation, chromatin accessibility, and RNA expression. We show that tissue specificity was much stronger than developmental stage specificity in all of DNA methylation, chromatin accessibility, and gene expression. We identify 573 genes showing allele specific expression, including those influenced by parent-of-origin as well as allele genotype effects. We integrate methylation, chromatin accessibility, and gene expression data to show that allele specific expression can be explained in great part by allele specific methylation and/or chromatin accessibility. This study provides a comprehensive characterization of regulatory variation across multiple tissues and developmental stages in pigs.

Unique feather color characteristics and transcriptome analysis of hair follicles in Liancheng White ducks

Avian feather color is a fascinating trait, and the genetic mechanism of duck plumage formation is still in the preliminary stage. In this study, feather color of Liancheng White ducks was analyzed by determination of melanin content and RNA-seq analysis. In this research, 9 ducks from Mallards (n = 3), Liancheng White (n = 3) and Pekin ducks (n = 3) were used by high performance liquid chromatography (HPLC) and Masson-Fontana staining to reveal the difference of feather melanin content. RNA-seq from 11 hair follicle tissues (1- and 8-wk-old) of Liancheng White ducks (n = 5) and Pekin ducks (n = 7) was used to analyze the candidate genes for the feather melanin synthesis, and Immunofluorescence experiment was used to show the protein expression in 6 black- and white-feathered ducks. Pectorale, skin, liver, fat, brain, heart, kidney, lung, spleen of an 8-wk-old black-feathered Mallard were collected for candidate gene expression. The results showed that the contents of feathers, beak, web melanin in Liancheng White ducks were higher than in Pekin ducks (p < 0.05). Melanin within hair follicles was located in the barb ridge and hair matrix of black feather duck, also we found that TYRP1, TYR, SOX10 genes were differentially expressed between Liancheng White and Pekin ducks (p < 0.05), and these genes were mainly expressed showed in duck skin tissues. This study revealed the unique feather color phenotype of Liancheng White duck shedding light on the transcriptome that underlies it.

Genome-wide association study reveals the candidate genes of humerus quality in laying duck

Bone plays a crucial role in poultry's health and production. However, during the selection and cage farming, there has been a decline in bone quality. As the development of breeding theory, researchers find that it's possible to enhance bone quality through selective breeding.This study measure 8 humerus quality in 260 samples of the 350-day-old female duck. By descripting the basic characteristic traits, mechanical property traits we found that all the bone quality traits had a large variable coefficient, especially mechanical properties trait (20-70%), indicating that there was a large difference in bone health status among laying ducks. The phenotypic correlations showed a high correlation between weight and density, diameter and perimeter, breaking and toughness (r = 0.52-0.68). And then, we performed the Genome-wide association study (GWAS) to reveal the candidate genes of humerus quality in ducks. Seven candidate protein-coding genes were identified with perimeter trait, and 52 protein-coding genes were associated with toughness trait. We also analysed the candidate region and performed KEGG and GO analyse for 75 candidate genes. Furthermore, the expression analyse of the above candidate genes in different stage of humerus and different tissues were performed. Finally, AP2A2, SMAD3, SMNDC1, NFIA, EPHB2, PMEPA1, UNC5C, ESR1, VAV3, NFATC2 deserve further focus. The obtained results can contribute to new insight into bone quality and provide new genetic biomarkers for application in duck breeding programs.

Whole-genome SNP allele frequency differences between Tibetan and Large white pigs reveal genes associated with skeletal muscle growth

BACKGROUND

The skeletal muscle growth rate and body size of Tibetan pigs (TIB) are lower than Large white pigs (LW). However, the underlying genetic basis attributing to these differences remains uncertain. To address this knowledge gap, the present study employed whole-genome sequencing of TIB (slow growth) and LW (fast growth) individuals, and integrated with existing NCBI sequencing datasets of TIB and LW individuals, enabling the identification of a comprehensive set of genetic variations for each breed. The specific and predominant SNPs in the TIB and LW populations were detected by using a cutoff value of 0.50 for SNP allele frequency and absolute allele frequency differences (△AF) between the TIB and LW populations.

RESULTS

A total of 21,767,938 SNPs were retrieved from 44 TIB and 29 LW genomes. The analysis detected 2,893,106 (13.29%) and 813,310 (3.74%) specific and predominant SNPs in the TIB and LW populations, and annotated to 24,560 genes. Further GO analysis revealed 291 genes involved in biological processes related to striated and/or skeletal muscle differentiation, proliferation, hypertrophy, regulation of striated muscle cell differentiation and proliferation, and myoblast differentiation and fusion. These 291 genes included crucial regulators of muscle cell determination, proliferation, differentiation, and hypertrophy, such as members of the Myogenic regulatory factors (MRF) (MYOD, MYF5, MYOG, MYF6) and Myocyte enhancer factor 2 (MEF2) (MEF2A, MEF2C, MEF2D) families, as well as muscle growth inhibitors (MSTN, ACVR1, and SMAD1); KEGG pathway analysis revealed 106 and 20 genes were found in muscle growth related positive and negative regulatory signaling pathways. Notably, genes critical for protein synthesis, such as MTOR, IGF1, IGF1R, IRS1, INSR, and RPS6KA6, were implicated in these pathways.

CONCLUSION

This study employed an effective methodology to rigorously identify the potential genes associated with skeletal muscle development. A substantial number of SNPs and genes that potentially play roles in the divergence observed in skeletal muscle growth between the TIB and LW breeds were identified. These findings offer valuable insights into the genetic underpinnings of skeletal muscle development and present opportunities for enhancing meat production through pig breeding.

Whole genome resequencing reveals genomic regions related to red plumage in ducks

Plumage color is a characteristic trait of ducks that originates as a result of natural and artificial selection. As a conspicuous phenotypic feature, it is a breed characteristic. Previous studies have identified some genes associated with the formation of black and white plumage in ducks. However, studies on the genetic basis underlying the red plumage phenotype in ducks are limited. Here, genome-wide association analysis (GWAS) and selection signal detection (Fst, θπ ratio, and cross-population composite likelihood ratio [XP-CLR]) were conducted to identify candidate regions and genes underlying duck plumage color phenotype. Selection signal detection revealed 29 overlapping genes (including ENPP1 and ULK1) significantly associated with red plumage color in Ji'an Red ducks. ENSAPLG00000012679, ESRRG, and SPATA5 were identified as candidate genes associated with red plumage using GWAS. Selection signal detection revealed that 19 overlapping genes (including GMDS, PDIA6, and ODC1) significantly correlated with light brown plumage in Brown Tsaiya ducks. GWAS to narrow down the significant regions further revealed nine candidate genes (AKT1, ATP6V1C2, GMDS, LRP4, MAML3, PDIA6, PLD5, TMEM63B, and TSPAN8). Notably, in Brown Tsaiya ducks, GMDS, ODC1, and PDIA6 exhibit significantly differentiated allele frequencies among other feather-colored ducks, while in Ji'an Red ducks, ENSAPLG00000012679 has different allele frequency distributions compared with that in other feather-colored ducks. This study offers new insights into the variation and selection of the red plumage phenotype using GWAS and selective signals.

The relationship of MITF gene expression and promoter methylation with plumage colour in quail

1. This study focused on the relationship between MITF mRNA expression and plumage colour in quail and the effect of promoter methylation on the expression of MITF mRNA.2. The CDS region of MITF mRNA was cloned by RT-PCR, followed by DNA sequencing. The RT-qPCR method was used to analyse the expression levels of MITF mRNA in dorsal skin tissue in Korean quail and Beijing white quail. The promoter region of the MITF gene was cloned, and the CpG island was predicted by the CpGplot program. The methylation levels of the CpG island were analysed using BS-PCR technology.3. Quail MITF mRNA contains a 1,476 bp complete ORF, which encodes a 492 amino acid residue protein. The MITF protein has no signal peptide or transmembrane region. The expression of MITF mRNA in dorsal tissue of Korean quail was significantly higher than that in Beijing white quail (p < 0.01). Abundant cis-elements and a 346 bp CpG island were found in the promoter region of the MITF gene. The average methylation level of the CpG island was 22 (22%) in Korean quail, and 46 (30%) in Beijing white quail (p < 0.05).4. The hypermethylation of the MITF gene promoter region in Beijing white quail resulted in a decrease in expression level, which was related to white feather colour.

How structural variants shape avian phenotypes: Lessons from model systems

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Divergent Regulatory Roles of Transcriptional Variants of the Chicken LDB3 Gene in Muscle Shaping

LIM domain binding 3 (LDB3) serves as a striated muscle-specific Z-band alternatively spliced protein that plays an important role in mammalian skeletal muscle development, but its regulatory role and molecular mechanism in avian muscle development are still unclear. In this study, we reanalyzed RNA sequencing data sets of 1415 samples from 21 chicken tissues published in the NCBI GEO database. First, three variants (LDB3-X, LDB3-XN1, and LDB3-XN2) generated by alternative splicing of the LDB3 gene were identified in chicken skeletal muscle, among which LDB3-XN1 and LDB3-XN2 are novel variants. LDB3-X and LDB3-XN1 are derived from exon skipping in chicken skeletal muscle at the E18-D7 stage and share three LIM domains, but LDB3-XN2 lacks a LIM domain. Our results preliminarily suggest that the formation of three variants of LDB3 is regulated by RBM20. The three splice isomers have divergent functions in skeletal muscle according to in vitro and in vivo assays. Finally, we identified the mechanism by which different variants play different roles through interactions with IGF2BP1 and MYHC, which promote the proliferation and differentiation of chicken myoblasts, in turn regulating chicken myogenesis. In conclusion, this study revealed the divergent roles of three LDB3 variants in chicken myogenesis and muscle remodeling and demonstrated their regulatory mechanism through protein-protein interactions.

Whole-genome selection signature differences between Chaohu and Ji'an red ducks

Assessing the genetic structure of local varieties and understanding their genetic data are crucial for effective management and preservation. However, the genetic differences among local breeds require further explanation. To enhance our understanding of their population structure and genetic diversity, we conducted a genome-wide comparative study of Chaohu and Ji'an Red ducks using genome sequence and restriction site-associated DNA sequencing technology. Our analysis revealed a distinct genetic distinction between the two breeds, leading to divided groups. The phylogenetic tree for Chaohu duck displayed two branches, potentially indicating minimal impact from artificial selection. Additionally, our ROH (runs of homozygosity) analysis revealed that Chaohu ducks had a lower average inbreeding coefficient than Ji'an Red ducks. We identified several genomic regions with high genetic similarity in these indigenous duck breeds. By conducting a selective sweep analysis, we identified 574 candidate genes associated with muscle growth (BMP2, ITGA8, MYLK, and PTCH1), fat deposits (ELOVL1 and HACD2), and pigmentation (ASIP and LOC101797494). These results offer valuable insights for the further enhancement and conservation of Chinese indigenous duck breeds.

Transcriptome analysis of sexual dimorphism in dorsal down coloration in goslings

BACKGROUND

In day-old Hungarian white goose goslings, there is a noticeable difference in dorsal down coloration between males and females, with females having darker dorsal plumage and males having lighter plumage. The ability to autosex day-old goslings based on their dorsal down coloration is important for managing them efficiently and planning their nutrition in the poultry industry. The aim of this study was to determine the biological and genetic factors underlying this difference in dorsal down colorationthrough histological analysis, biochemical assays, transcriptomic profiling, and q‒PCR analysis.

RESULTS

Tissue analysis and biochemical assays revealed that compared with males, 17-day-old embryos and day-old goslings of female geese exhibited a greater density of melanin-containing feather follicles and a greater melanin concentration in these follicles during development. Both female and male goslings had lower melanin concentrations in their dorsal skin compared to 17-day-old embryos. Transcriptome analysis identified a set of differentially expressed genes (DEGs) (MC1R, TYR, TYRP1, DCT and MITF) associated with melanogenesis pathways that were downregulated or silenced specifically in the dorsal skin of day-old goslings compared to 17-day-old embryos, affecting melanin synthesis in feather follicles. Additionally, two key genes (MC1R and MITF) associated with feather coloration showed differences between males and females, with females having higher expression levels correlated with increased melanin synthesis and darker plumage.

CONCLUSION

The expression of multiple melanogenesis genes determines melanin synthesis in goose feather follicles. The dorsal down coloration of day-old Hungarian white goose goslings shows sexual dimorphism, likely due to differences in the expression of the MC1R and MITF genes between males and females. These results could help us better understand why male and female goslings exhibit different plumage patterns.

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet

BACKGROUND

Goose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended.

RESULTS

In this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake.

CONCLUSIONS

We examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits

BACKGROUND

Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits.

RESULTS

In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits.

CONCLUSIONS

Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production 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.

Elucidation of the genetic determination of body weight and size in Chinese local chicken breeds by large-scale genomic analyses

BACKGROUND

Body weight and size are important economic traits in chickens. While many growth-related quantitative trait loci (QTLs) and candidate genes have been identified, further research is needed to confirm and characterize these findings. In this study, we investigate genetic and genomic markers associated with chicken body weight and size. This study provides new insights into potential markers for genomic selection and breeding strategies to improve meat production in chickens.

METHODS

We performed whole-genome resequencing of and Wenshang Barred (WB) chickens (n = 596) and three additional breeds with varying body sizes (Recessive White (RW), WB, and Luxi Mini (LM) chickens; (n = 50)). We then used selective sweeps of mutations coupled with genome-wide association study (GWAS) to identify genomic markers associated with body weight and size.

RESULTS

We identified over 9.4 million high-quality single nucleotide polymorphisms (SNPs) among three chicken breeds/lines. Among these breeds, 287 protein-coding genes exhibited positive selection in the RW and WB populations, while 241 protein-coding genes showed positive selection in the LM and WB populations. Genomic heritability estimates were calculated for 26 body weight and size traits, including body weight, chest breadth, chest depth, thoracic horn, body oblique length, keel length, pelvic width, shank length, and shank circumference in the WB breed. The estimates ranged from 0.04 to 0.67. Our analysis also identified a total of 2,522 genome-wide significant SNPs, with 2,474 SNPs clustered around two genomic regions. The first region, located on chromosome 4 (7.41-7.64 Mb), was linked to body weight after ten weeks and body size traits. LCORL, LDB2, and PPARGC1A were identified as candidate genes in this region. The other region, located on chromosome 1 (170.46-171.53 Mb), was associated with body weight from four to eighteen weeks and body size traits. This region contained CAB39L and WDFY2 as candidate genes. Notably, LCORL, LDB2, and PPARGC1A showed highly selective signatures among the three breeds of chicken with varying body sizes.

CONCLUSION

Overall this study provides a comprehensive map of genomic variants associated with body weight and size in chickens. We propose two genomic regions, one on chromosome 1 and the other on chromosome 4, that could helpful for developing genome selection breeding strategies to enhance meat yield in chickens.

Genome-wide association study identifies multiple loci influencing duck serum biochemical indicators in the laying period

1. Laying performance is an important economic trait in poultry. The blood is essential in transporting nutrients to the yolk and albumen and is necessary for egg formation.2. This study calculated the phenotypic relationships of duck egg quality, egg production efficiency and 22 serum parameters in the egg-laying stage. Using a variety of methodologies, a genome-wide association study (GWAS) was carried out to uncover the genetic foundations of the 22 serum biochemical markers of laying ducks.3. Spearman correlation coefficients between the egg production (226-329 per day) and the serum parameters were all weak, being less than 0.3. This analysis was done on 22 serum parameters, with total protein (TP), total triglycerides (TG), calcium (Ca) and phosphorous (P) having the highest correlation coefficients (r = 0.56-0.88). The coefficients for blood markers, such as total cholesterol (CHOL), total bilirubin (TBIL), low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) varied from 0.70-0.94.4. Based on single-marker single-trait genome-wide analyses by a mixed linear model program of EMMAX, nine candidate genes were associated with enzyme traits (AST/ALT aspartate transaminase/glutamic-pyruvic transaminase, creatine kinase) and 19 candidate genes were associated with metabolism and protein-related serum parameters (glucose, total bile acid, uric acid (UA), albumin (ALB).5. The mvLMM (multivariate linear mixed model) of GEMMA software was used to carry out multiple trait integrated GWAS. Two candidate genes were found in the TP-TG-CA-P analysis and seven candidate genes in the CHOL_LDL-C_HDL-C_TBIL study. There was a high genetic correlation between the two groups.

Transcriptome-based comparison reveals key genes regulating allometry growth of forelimb and hindlimb bone in duck embryos

Allometric growth of the forelimb and hindlimb is a widespread phenomenon observed in vertebrates. As a typical precocial bird, ducks exhibit more advanced development of their hindlimbs compared to their forelimbs, enabling them to walk shortly after hatching. This phenomenon is closely associated with the development of long bones in the embryonic stage. However, the molecular mechanism governing the allometric growth of duck forelimb and hindlimb bones is remains elusive. In this study, we employed phenotypic, histological, and gene expression analyses to investigate developmental differences between the humerus (forelimb bone) and tibia/femur (hindlimb bones) in duck embryos. Our results revealed a gradual increase in weight and length disparity between the tibia and humerus from E12 to E28 (embryo age). At E12, endochondral ossification was observed solely in the tibia but not in the humerus. The number of differentially expressed genes (DEGs) gradually increased at H12 vs. T12, H20 vs. T20, and H28 vs. T28 stages consistent with phenotypic variations. A total of 38 DEGs were found across all 3 stages. Protein-protein interaction network analysis demonstrated strong interactions among members of HOXD gene family (HOXD3/8/9/10/11/12), HOXB gene family (HOXB8/9), TBX gene family (TBX4/5/20), HOXA11, SHOX2, and MEIS2. Gene expression profiling indicated higher expression levels for all HOXD genes in the humerus compared to tibia while opposite trends were observed for HOXA/HOXB genes with low or no expression detected in the humerus. These findings suggest distinct roles played by different clusters within HOX gene family during skeletal development regulation of duck embryo's forelimbs versus hind limbs. Notably, TBX4 exhibited high expression levels specifically in tibia whereas TBX5 showed similar patterns exclusively within humerus as seen previously across other species' studies. In summary, this study identified key regulatory genes involved in allometric growth of duck forelimb and hindlimb bones during embryonic development. Skeletal development is a complex physiological process, and further research is needed to elucidate the regulatory role of candidate genes in endochondral ossification.

Duck pan-genome reveals two transposon insertions caused bodyweight enlarging and white plumage phenotype formation during evolution

Structural variations (SVs) are a major source of domestication and improvement traits. We present the first duck pan-genome constructed using five genome assemblies capturing ∼40.98 Mb new sequences. This pan-genome together with high-depth sequencing data (∼46.5×) identified 101,041 SVs, of which substantial proportions were derived from transposable element (TE) activity. Many TE-derived SVs anchoring in a gene body or regulatory region are linked to duck's domestication and improvement. By combining quantitative genetics with molecular experiments, we, for the first time, unraveled a 6945 bp Gypsy insertion as a functional mutation of the major gene IGF2BP1 associated with duck bodyweight. This Gypsy insertion, to our knowledge, explains the largest effect on bodyweight among avian species (27.61% of phenotypic variation). In addition, we also examined another 6634 bp Gypsy insertion in MITF intron, which triggers a novel transcript of MITF, thereby contributing to the development of white plumage. Our findings highlight the importance of using a pan-genome as a reference in genomics studies and illuminate the impact of transposons in trait formation and livestock breeding.

Genetic Diversity, Population Structure and Selection Signature in Begait Goats Revealed by Whole-Genome Sequencing

Goats belong to a group of animals called small ruminants and are critical sources of livelihood for rural people. Genomic sequencing can provide information ranging from basic knowledge about goat diversity and evolutionary processes that shape genomes to functional information about genes/genomic regions. In this study, we exploited a whole-genome sequencing data set to analyze the genetic diversity, population structure and selection signatures of 44 individuals belonging to 5 Ethiopian goat populations: 12 Aberegalle (AB), 5 Afar (AF), 11 Begait (BG), 12 Central highlands (CH) and 5 Meafure (MR) goats. Our results revealed the highest genetic diversity in the BG goat population compared to the other goat populations. The pairwise genetic differentiation (FST) among the populations varied and ranged from 0.011 to 0.182, with the closest pairwise value (0.003) observed between the AB and CH goats and a distant correlation (FST = 0.182) between the BG and AB goats, indicating low to moderate genetic differentiation. Phylogenetic tree, ADMIXTURE and principal component analyses revealed a classification of the five Ethiopian goat breeds in accordance with their geographic distribution. We also found three top genomic regions that were detected under selection on chromosomes 2, 5 and 13. Moreover, this study identified different candidate genes related to milk characteristics (GLYCAM1 and SRC), carcass (ZNF385B, BMP-7, PDE1B, PPP1R1A, FTO and MYOT) and adaptive and immune response genes (MAPK13, MAPK14, SCN7A, IL12A, EST1 DEFB116 and DEFB119). In conclusion, this information could be helpful for understanding the genetic diversity and population structure and selection scanning of these important indigenous goats for future genetic improvement and/or as an intervention mechanism.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Genome-wide detection of selection signatures in Jianli pigs reveals novel cis-regulatory haplotype in EDNRB associated with two-end black coat color

BACKGROUND

Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown.

RESULTS

In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function.

CONCLUSION

Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs.

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens

INTRODUCTION

Investigating the genetic markers and genomic signatures related to chicken meat production by combing multi-omics methods could provide new insights into modern chicken breeding technology systems.

OBJECT

Chicken is one of the most efficient and environmentally friendly livestock, especially the fast-growing white-feathered chicken (broiler), which is well known for high meat yield, but the underlying genetic basis is poorly understood.

METHOD

We generated whole-genome resequencing of three purebred broilers (n = 748) and six local breeds/lines (n = 114), and sequencing data of twelve chicken breeds (n = 199) were obtained from the NCBI database. Additionally, transcriptome sequencing of six tissues from two chicken breeds (n = 129) at two developmental stages was performed. A genome-wide association study combined with cis-eQTL mapping and the Mendelian randomization was applied.

RESULT

We identified > 17 million high-quality SNPs, of which 21.74% were newly identified, based on 21 chicken breeds/lines. A total of 163 protein-coding genes underwent positive selection in purebred broilers, and 83 genes were differentially expressed between purebred broilers and local chickens. Notably, muscle development was proven to be the major difference between purebred broilers and local chickens, or ancestors, based on genomic and transcriptomic evidence from multiple tissues and stages. The MYH1 gene family showed the top selection signatures and muscle-specific expression in purebred broilers. Furthermore, we found that the causal gene SOX6 influenced breast muscle yield and also related to myopathy occurrences. A refined haplotype was provided, which had a significant effect on SOX6 expression and phenotypic changes.

CONCLUSION

Our study provides a comprehensive atlas comprising the typical genomic variants and transcriptional characteristics for muscle development and suggests a new regulatory target (SOX6-MYH1s axis) for breast muscle yield and myopathy, which could aid in the development of genome-scale selective breeding aimed at high meat yield in broiler chickens.