Temporal variation in production performance, biochemical and oxidative stress markers, and gut microbiota in Pekin ducks during the late growth stage.

In the late growth stage of commercial Pekin ducks, a significant increase in feed intake and a decline in body weight gain have been observed, leading to impaired feed conversion efficiency. To address this issue, we investigated alterations in production performance, blood biochemical indices, ileum tissue architecture, and microbial community structure in Pekin ducks. The primary objective was to provide robust data supporting the improvement of meat duck production efficiency during the late growth stage (28-42-days-old). Forty 28-day-old Pekin ducks were randomly assigned to 8 replicates, with five ducks per replicate. The rearing period lasted 14 days, with feed and water provided ad libitum. Our findings indicated a significant increase in Pekin duck body and heart weights with advancing age (P < 0.05). Moreover, serum antioxidant enzyme and high-density lipoprotein concentrations significantly increased, whereas triglyceride levels decreased (P < 0.05). Notably, the height of the ileal villi was significantly reduced (P < 0.05). The microbial community structure of the ileum exhibited significant changes as ducks aged, accompanied by a substantial increase in microbial flora diversity, particularly with the formation of more tightly connected microbial network modules. Time-dependent enrichment was observed in microbial gene functions related to energy metabolism pathways. At the genus level, Sphingomonas and Subdoligranulum have emerged as crucial players in microbial differential functional pathways and network formation. These bacteria likely serve as the key driving factors in the dynamic microbial changes that occur in Pekin ducks over time. Overall, our findings suggest a potential decline in the absorption function of the small intestine and fat deposition performance of Pekin ducks during later growth stages, which may be attributed to the maturation and proliferation of the gut microbial community.

Genetic Foundation of Male Spur Length and Its Correlation with Female Egg Production in Chickens.

Spurs, which mainly appear in roosters, are protrusions near the tarsometatarsus on both sides of the calves of chickens, and are connected to the tarsometatarsus by a bony core. As a male-biased morphological characteristic, the diameter and length of spurs vary significantly between different individuals, mainly related to genetics and age. As a specific behavior of hens, egg-laying also varies greatly between individuals in terms of traits such as age at first egg (AFE), egg weight (EW), and so on. At present, there are few studies on chicken spurs. In this study, we investigated the inheritance pattern of the spur trait in roosters with different phenotypes and the correlations between spur length, body weight at 18 weeks of age (BW18), shank length at 18 weeks of age (SL18), and the egg-laying trait in hens (both hens and roosters were from the same population and were grouped according to their family). These traits related to egg production included AFE, body weight at first egg (BWA), and first egg weight (FEW). We estimated genetic parameters based on pedigree and phenotype data, and used variance analysis to calculate broad-sense heritability for correcting the parameter estimation results. The results showed that the heritability of male left and right spurs ranged from 0.6 to 0.7. There were significant positive correlations between left and right spur length, BW18, SL18, and BWA, as well as between left and right spur length and AFE. We selected 35 males with the longest spurs and 35 males with the shortest spurs in the population, and pooled them into two sets to obtain the pooled genome sequencing data. After genome-wide association and genome divergency analysis by FST, allele frequency differences (AFDs), and XPEHH methods, we identified 7 overlapping genes (CENPE, FAT1, FAM149A, MANBA, NFKB1, SORBS2, UBE2D3) and 14 peak genes (SAMD12, TSPAN5, ENSGALG00000050071, ENSGALG00000053133, ENSGALG00000050348, CNTN5, TRPC6, ENSGALG00000047655,TMSB4X, LIX1, CKB, NEBL, PRTFDC1, MLLT10) related to left and right spur length through genome-wide selection signature analysis and a genome-wide association approach. Our results identified candidate genes associated with chicken spurs, which helps to understand the genetic mechanism of this trait and carry out subsequent research around it.

Increased Gremlin1 Expression in Pancreatic Ductal Adenocarcinoma Promotes a Fibrogenic Stromal Microenvironment.

OBJECTIVE: The pancreatic ductal adenocarcinoma (PDAC) microenvironment is primarily composed of cancer-associated fibroblasts (CAFs) and immune cells. Gremlin1 (Grem1) is a profibrogenic factor that promotes tumorigenesis in several cancers. However, the role of Grem1 in the PDAC microenvironment is not adequately defined. METHODS: We correlated Grem1 levels with activated stroma and immune cells in human PDAC using The Cancer Genome Atlas (TCGA) RNA-sequencing data and characterized the expression of Grem1 transcripts and isoforms in pancreatic cell lines and PDAC tissues. We assessed the role of Grem1 in the microenvironment by in vitro studies. RESULTS: Grem1 expression is associated with an activated stroma and increased M1 and M2 macrophages. Only full length Grem1 variant 1 and isoform 1 were detectable in human pancreatic cells, and remarkably high levels of Grem1 were observed in pancreatic fibroblasts (P < 0.05). Immunohistochemistry detected Grem1 protein in PDAC tumor cells and stromal cells, which correlated with infiltrating macrophages in PDAC tumors. Grem1 knockdown in CAFs suppressed transforming growth factor (TGF)-ß-induced extracellular matrix proteins (P < 0.05). Grem1 recombinant protein treatment in vitro increased M1 and M2 macrophages (P < 0.05). CONCLUSIONS: Grem1 acts as a profibrogenic factor in the PDAC microenvironment via modulation of fibroblasts and macrophages. Grem1 may have the potential to be developed as a therapeutic target for PDAC.

Genetic parameter estimation and molecular foundation of chicken beak shape.

The bird beak is mainly functioned as feeding and attacking, and its shape has extremely important significance for survival and reproduction. In chickens, since beak shape could lead to some disadvantages including pecking and waste of feed, it is important to understand the inheritance of chicken beak shape. In the present study, we firstly established 4 indicators to describe the chicken beak shapes, including upper beak length (UL), lower beak length (LL), distance between upper and lower beak tips (DB) and upper beak curvature (BC). And then, we measured the 4 beak shape indicators as well as some production traits including body weight (BW), shank length (SL), egg weight (EW), eggshell strength (ES) of a layer breed, Rhode Island Red (RIR), in order to estimate genetic parameters of chicken beak shape. The heritabilities of UL and LL were 0.41 and 0.37, and the heritabilities of DB and BC were 0.22 and 0.21, indicating that beak shape was a highly or mediumly heritable. There were significant positive genetic and phenotypic correlations among UL, LL, and DB. And UL was positively correlated with body weight (BW18) and shank length (SL18) at 18 weeks of age in genetics, and DB was positively correlated with BC in terms of genetics and phenotype. We also found that layers of chicken cages played a role on beak shape, which could be attributed to the difference of lightness in different cage layers. By a genome-wide association study (GWAS) for the chicken UL, we identified 9 significant candidate genes associated with UL in RIR. For the variants with low minor allele frequencies (MAF <0.01) and outside of high linkage disequilibrium (LD) regions, we also conducted rare variant association studies (RVA) and GWAS to find the association between genotype and phenotype. We also analyzed transcriptomic data from multiple tissues of chicken embryos and revealed that all of the 9 genes were highly expressed in beak of chicken embryos, indicating their potential function for beak development. Our results provided the genetic foundation of chicken beak shape, which could help chicken breeding on beak related traits.

Genetic parameter estimation and molecular foundation of chicken egg-laying trait.

The age of first egg (AFE) in chicken can affect early and even life-time egg production performance to some extent, and therefore is an important economic trait that affects production efficiency. To better understand the genetic patterns of AFE and other production traits including body weight at first egg (BWA), first egg weight (FEW), and total egg number from AFE to 58 wk of age (total-EN), we recorded the production performance of 2 widely used layer breeds, white leghorn (WL) and Rhode Island Red (RIR) and estimated genetic parameters based on pedigree and production data. The results showed that the heritability of AFE in both breeds ranged from 0.4 to 0.6, and AFE showed strong positive genetic and phenotypic correlations to BWA as well as FEW, while showing strong negative genetic and phenotypic correlations with total-EN. Furtherly, by genome-wide association analysis study (GWAS), we identified 12 and 26 significant SNPs to be related to AFE in the 2-layer breeds, respectively. A total of 18 genes were identified that could affect AFE based on the significant SNP annotations obtained, but there were no gene overlapped in the 2 breeds indicating the genetic foundation of AFE could differ from breed to breed. Our results provided a deeper understanding of genetic patterns and molecular basement of AFE in different breeds and could help in the selection of egg production traits.

Molecular genetic foundation of a sex-linked tailless trait in Hongshan chicken by whole genome data analysis.

As a Chinese local chicken breed, Hongshan chickens have 2 kinds of tail feather phenotypes, normal and taillessness. Our previous studies showed that taillessness was a sex-linked dominant trait. Abnormal development of the tail vertebrae could be explained this phenomenon in some chicken breeds. However, the number of caudal vertebrae in rumpless Hongshan chickens was normal, so rumplessness in Hongshan chicken was not related to the development of the caudal vertebrae. Afterwards, we found that rumplessness in Hongshan was due to abnormal development of tail feather rather than abnormal development of caudal vertebrae. In order to understand the genetic foundation of the rumplessness of Hongshan chickens, we compared and reanalyzed 2 sets of data in normal and rumpless Hongshan chickens from our previous studies. By joint analysis of genome-wide selection signature analysis and genome-wide association approach, we found that 1 overlapping gene (EDIL3) and 16 peak genes (ENSGALG00000051843, ENSGALG00000053498, ENSGALG00000054800, KIF27, PTPRD, ENSGALG00000047579, ENSGALG00000041052, ARHGEF28, CAMK4, SERINC5, ENSGALG00000050776, ERCC8, MCC, ADAMTS19, ENSGALG00000053322, CHRNA8) located on the Z chromosome was associated with the rumpless trait. The results of this study furtherly revealed the molecular mechanism of the rumpless trait in Hongshan chickens, and identified the candidate genes associated with this trait. Our results will help to improve the shape of chicken tail feathers and to rise individual economic value in some specific market in China.

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.

Significant genomic introgression from grey junglefowl (Gallus sonneratii) to domestic chickens (Gallus gallus domesticus).

BACKGROUND: Chicken is one of the most numerous and widely distributed species around the world, and many studies support the multiple ancestral origins of domestic chickens. The research regarding the yellow skin phenotype in domestic chickens (regulated by BCO2) likely originating from the grey junglefowl serves as crucial evidence for demonstrating the multiple origins of chickens. However, beyond the BCO2 gene region, much remains unknown about the introgression from the grey junglefowl into domestic chickens. Therefore, in this study, based on whole-genome data of 149 samples including 4 species of wild junglefowls and 13 local domestic chicken breeds, we explored the introgression events from the grey junglefowl to domestic chickens. RESULTS: We successfully detected introgression regions besides BCO2, including two associated with growth trait (IGFBP2 and TKT), one associated with angiogenesis (TIMP3) and two members of the heat shock protein family (HSPB2 and CRYAB). Our findings suggest that the introgression from the grey junglefowl may impact the growth performance of chickens. Furthermore, we revealed introgression events from grey junglefowl at the BCO2 region in multiple domestic chicken breeds, indicating a phenomenon where the yellow skin phenotype likely underwent strong selection and was retained. Additionally, our haplotype analysis shed light on BCO2 introgression event from different sources of grey junglefowl into domestic chickens, possibly suggesting multiple genetic flows between the grey junglefowl and domestic chickens. CONCLUSIONS: In summary, our findings provide evidences of the grey junglefowl contributing to the genetic diversity of domestic chickens, laying the foundation for a deeper understanding of the genetic composition within domestic chickens, and offering new perspectives on the impact of introgression on domestic chickens.

Dynamic Imaging of Blood Coagulation Within the Hematoma of Patients With Acute Hemorrhagic Stroke.

BACKGROUND: The dynamics of blood clot (combination of Hb [hemoglobin], fibrin, and a higher concentration of aggregated red blood cells) formation within the hematoma of an intracerebral hemorrhage is not well understood. A quantitative neuroimaging method of localized coagulated blood volume/distribution within the hematoma might improve clinical decision-making. METHODS: The deoxyhemoglobin of aggregated red blood cells within extravasated blood exhibits a higher magnetic susceptibility due to unpaired heme iron electrons. We propose that coagulated blood, with higher aggregated red blood cell content, will exhibit (1) a higher positive susceptibility than noncoagulated blood and (2) increase in fibrin polymerization-restricted localized diffusion, which can be measured noninvasively using quantitative susceptibility mapping and diffusion tensor imaging. In this serial magnetic resonance imaging study, we enrolled 24 patients with acute intracerebral hemorrhage between October 2021 to May 2022 at a stroke center. Patients were 30 to 70 years of age and had a hematoma volume >15 cm3 and National Institutes of Health Stroke Scale score >1. The patients underwent imaging 3×: within 12 to 24 (T1), 36 to 48 (T2), and 60 to 72 (T3) hours of last seen well on a 3T magnetic resonance imaging system. Three-dimensional anatomic, multigradient echo and 2-dimensional diffusion tensor images were obtained. Hematoma and edema volumes were calculated, and the distribution of coagulation was measured by dynamic changes in the susceptibilities and fractional anisotropy within the hematoma. RESULTS: Using a coagulated blood phantom, we demonstrated a linear relationship between the percentage coagulation and susceptibility (R2=0.91) with a positive red blood cell stain of the clot. The quantitative susceptibility maps showed a significant increase in hematoma susceptibility (T1, 0.29±0.04 parts per millions; T2, 0.36±0.04 parts per millions; T3, 0.45±0.04 parts per millions; P<0.0001). A concomitant increase in fractional anisotropy was also observed with time (T1, 0.40±0.02; T2, 0.45±0.02; T3, 0.47±0.02; P<0.05). CONCLUSIONS: This quantitative neuroimaging study of coagulation within the hematoma has the potential to improve patient management, such as safe resumption of anticoagulants, the need for reversal agents, the administration of alteplase to resolve the clot, and the need for surgery.

Clearance of Neutrophils From ICH-Affected Brain by Macrophages Is Beneficial and Is Assisted by Lactoferrin and CD91.

BACKGROUND: Within hours after intracerebral hemorrhage (ICH) onset, masses of polymorphonuclear neutrophils (PMNs) infiltrate the ICH-affected brain. After degranulation involving controlled release of many toxic antimicrobial molecules, the PMNs undergo rapid apoptosis and then are removed by phagocytic microglia/macrophages (MΦ) through a process called efferocytosis. Effective removal of PMNs may limit secondary brain damage and inflammation; however, the molecular mechanisms governing these cleanup activities are not well understood. We propose that scavenger receptor CD91 on myeloid phagocytes especially in presence of CD91 ligand, LTF (lactoferrin, protein abundant in PMNs), plays an important role in clearance of dead apoptotic PMNs (ANs). METHODS: Mice/rats were subjected to an autologous blood injection model of ICH. Primary cultured microglia were used to assess phagocytosis of ANs. Immunohistochemistry was employed to assess CD91 expression and PMN infiltration. CD91 knockout mice selectively in myeloid phagocytes (Mac-CD91-KO) were used to establish the CD91/LTF function in phagocytosis and in reducing ICH-induced injury, as assessed using behavioral tests, hematoma resolution, and oxidative stress. RESULTS: Masses of PMNs are found in ICH-affected brain, and they contain LTF. MΦ at the outer border of hematoma are densely packed, expressing CD91 and phagocytosing ANs. Microglia deficient in CD91 demonstrate defective phagocytosis of ANs, and mice deficient in CD91 (Mac-CD91-KO) subjected to ICH injury have increased neurological dysfunction that is associated with impaired hematoma resolution (hemoglobin and iron clearance) and elevated oxidative stress. LTF that normally ameliorates ICH injury in CD91-proficient control mice shows reduced therapeutic effects in Mac-CD91-KO mice. CONCLUSIONS: Our study suggests that CD91 plays a beneficial role in improving ANs phagocytosis and ultimately post-ICH outcome and that the beneficial effect of LTF in ICH is in part dependent on presence of CD91 on MΦ.

Genomic Analysis Reveals Candidate Genes Underlying Sex-Linked Eyelid Coloboma, Feather Color Traits, and Climatic Adaptation in Huoyan Geese.

Driven by natural and artificial selection, the domestic Huoyan geese from Northern China have gradually generated specific phenotypes and climatic adaptations. To understand the genetic basis of the two specific phenotypes that are sex linked, including upper eyelid coloboma and gosling feather color, as well as the climatic adaptations of the Huoyan goose, which can contribute to the artificial selection and breeding of geese. We selected Huoyan geese and nine Southern Chinese goose breeds and identified their divergence on the genomic level. Using selective sweep analysis, we found that PTPRM on chromosome Z influences the upper eyelid coloboma phenotype of the Huoyan goose, and TYRP1 is a plausible candidate gene for the Huoyan gosling feather color. We obtained a number of genes related to cold adaptation in Huoyan geese, mainly involved in physiological functions such as metabolism, angiogenesis contraction and circulatory system, apoptosis, immunity, stress, and neural system. The most interesting candidates for cold adaptation are PIP5K1B and NMNAT3 that are associated with energy metabolism and stress. We also obtained some genes related to heat adaptation, including AGTPBP1, associated with neurology; GDA, associated with skin pigmentation; and NAA35, associated with apoptosis. These findings deepen our understanding of the genetics of specific phenotypes and climate adaptation in local geese and provide insights for the selection of goose breeds.

Evolution, Inheritance, and Strata Formation of the W Chromosome in Duck (Anas platyrhynchos).

The nonrecombining female-limited W chromosome is predicted to experience unique evolutionary processes. Difficulties in assembling W chromosome sequences have hindered the identification of duck W-linked sequences and their evolutionary footprint. To address this, we conducted three initial contig-level genome assemblies and developed a rigorous pipeline by which to successfully expand the W-linked data set, including 11 known genes and 24 newly identified genes. Our results indicate that the W chromosome expression may not be subject to female-specific selection; a significant convergent pattern of upregulation associated with increased female-specific selection was not detected. The genetic stability of the W chromosome is also reflected in the strong evolutionary correlation between it and the mitochondria; the complete consistency of the cladogram topology constructed from their gene sequences proves the shared maternal coevolution. By detecting the evolutionary trajectories of W-linked sequences, we have found that recombination suppression started in four distinct strata, of which three were conserved across Neognathae. Taken together, our results have revealed a unique evolutionary pattern and an independent stratum evolutionary pattern for sex chromosomes.

Phylogenetic analysis reveals multiple origins of Chinese gamecocks.

Cockfighting is popular worldwide, dating back to 2,800 BC. Primarily, 5 modern Chinese gamecock breeds exist, located in the northeast (Luxi and Henan), west (Turpan), south (Xishuangbanna), and southeast (Zhangzhou) of China. However, whether Chinese gamecocks were derived from a single origin or multiple origins remains controversial. Therefore, this study used next-generation resequencing data to elucidate the origin of Chinese gamecocks by constructing genome-wide and SRY-box transcription factor 5 (SOX5) gene phylogenetic trees. Data from 161 chickens from 27 breeds, including 9 gamecock breeds, were included. Before constructing the SOX5 gene tree, we validated that the pea-comb phenotype mutation in all gamecock breeds was attributed to copy number variation in intron 1 of the SOX5 gene, as previously reported. The specific region was chr1: 65,838,000 to 65,846,000. The phylogenetic tree results suggested that Zhangzhou and Xishuangbanna gamecocks have a monophyletic origin, while Luxi, Henan, and Turpan gamecocks have a common ancestor. Our study provides genome-wide evidence that Chinese gamecocks have multiple origins and advances the understanding of the genetic mechanisms of the pea-comb characteristic.

Comparative population genomics analysis uncovers genomic footprints and genes influencing body weight trait in Chinese indigenous chicken.

Body weight of chicken is a typical quantitative trait, which shows phenotypic variations due to selective breeding. Despite some QTL loci have been obtained, the body weight of native chicken breeds in different geographic regions varies greatly, its genetic basis remains unresolved questions. To address this issue, we analyzed 117 Chinese indigenous chickens from 10 breeds (Huiyang Bearded, Xinhua, Hotan Black, Baicheng You, Liyang, Yunyang Da, Jining Bairi, Lindian, Beijing You, Tibetan). We applied fixation index (FST) analysis to find selected genomic regions and genes associated with body weight traits. Our study suggests that NELL1, XYLT1, and NCAPG/LCORL genes are strongly selected in the body weight trait of Chinese indigenous chicken breeds. In addition, the IL1RAPL1 gene was strongly selected in large body weight chickens, while the PCDH17 and CADM2 genes were strongly selected in small body weight chickens. This result suggests that the patterns of genetic variation of native chicken and commercial chicken, and/or distinct local chicken breeds may follow different evolutionary mechanisms.

Extensive intra- and inter-genetic admixture of Chinese gamecock and other indigenous chicken breeds revealed by genomic data.

Genomic admixture is a widespread phenomenon among domestic animal breeds, including chickens. However, reports on admixture within Chinese gamecocks or other indigenous chickens are limited. This study focuses on the population genetic structure and admixture of 5 Chinese gamecock breeds and the admixture with 9 other indigenous Chinese chicken breeds. Our results showed that Turpan and Henan gamecocks were grouped into one cluster, whereas Luxi, Zhangzhou, and Xishuangbanna gamecocks were grouped into the other cluster. Gene flow occurred between Xishuangbanna and Turpan and Turpan and Luxi gamecocks. Simultaneously, gene flow was observed between gamecocks and indigenous chickens, such as Xishuangbanna and Wenchang. Ancestral component analysis indicated that modern domestic chickens in southern China played an important role in the history of the domestication of modern Chinese gamecock. Our study will be helpful in better understanding the domestication and evolution of Chinese gamecock.

Systematic Selection Signature Analysis of Chinese Gamecocks Based on Genomic and Transcriptomic Data.

Selection pressures driven by natural causes or human interference are key factors causing genome variants and signatures of selection in specific regions of the genome. Gamecocks were bred for cockfighting, presenting pea-combs, larger body sizes, stronger limbs, and higher levels of aggression than other chickens. In this study, we aimed to explore the genomic differences between Chinese gamecocks and commercial, indigenous, foreign, and cultivated breeds by detecting the regions or sites under natural or artificial selection using genome-wide association studies (GWAS), genome-wide selective sweeps based on the genetic differentiation index (FST), and transcriptome analyses. Ten genes were identified using GWAS and FST: gga-mir-6608-1, SOX5, DGKB, ISPD, IGF2BP1, AGMO, MEOX2, GIP, DLG5, and KCNMA1. The ten candidate genes were mainly associated with muscle and skeletal development, glucose metabolism, and the pea-comb phenotype. Enrichment analysis results showed that the differentially expressed genes between the Luxi (LX) gamecock and Rhode Island Red (RIR) chicken were mainly related to muscle development and neuroactive-related pathways. This study will help to understand the genetic basis and evolution of Chinese gamecocks and support the further use of gamecocks as an excellent breeding material from a genetic perspective.

Candidate pigmentation genes related to feather color variation in an indigenous chicken breed revealed by whole genome data.

Chicken plumage color is an inheritable phenotype that was naturally and artificially selected for during domestication. The Baicheng You chicken is an indigenous Chinese chicken breed presenting three main feather colors, lavender, black, and yellow plumages. To explore the genetic mechanisms underlying the pigmentation in Baicheng You chickens, we re-sequenced the whole genome of Baicheng You chicken with the three plumage colors. By analyzing the divergent regions of the genome among the chickens with different feather colors, we identified some candidate genomic regions associated with the feather colors in Baicheng You chickens. We found that EGR1, MLPH, RAB17, SOX5, and GRM5 genes were the potential genes for black, lavender, and yellow feathers. MLPH, GRM5, and SOX5 genes have been found to be related to plumage colors in birds. Our results showed that EGR1 is a most plausible candidate gene for black plumage, RAB17, MLPH, and SOX5 for lavender plumage, and GRM5 for yellow plumage in Baicheng You chicken.

Comprehensive analysis of structural variants in chickens using PacBio sequencing.

Structural variants (SVs) are one of the main sources of genetic variants and have a greater impact on phenotype evolution, disease susceptibility, and environmental adaptations than single nucleotide polymorphisms (SNPs). However, SVs remain challenging to accurately type, with several detection methods showing different limitations. Here, we explored SVs from 10 different chickens using PacBio technology and detected 49,501 high-confidence SVs. The results showed that the PacBio long-read detected more SVs than Illumina short-read technology genomes owing to some SV sites on chromosomes, which are related to chicken growth and development. During chicken domestication, some SVs beneficial to the breed or without any effect on the genomic function of the breed were retained, whereas deleterious SVs were generally eliminated. This study could facilitate the analysis of the genetic characteristics of different chickens and provide a better understanding of their phenotypic characteristics at the SV level, based on the long-read sequencing method. This study enriches our knowledge of SVs in chickens and improves our understanding of chicken genomic diversity.

A ∼ 4.1 kb deletion in IRX1 gene upstream is completely associated with rumplessness in Piao chicken.

Piao chicken, a Chinese indigenous rumpless chicken breed, lacks pygostyle, caudal vertebra, uropygial gland and tail feathers. The rumplessness in Piao chicken presents an autosomal dominant inheritance pattern. However, the molecular genetic mechanisms underlying the rumplessness in Piao chicken remains unclear. In this study, whole-genome resequencing was performed for 146 individuals from 10 chicken breeds, including 9 tailed chicken breeds and Piao rumpless breed. Tailbone CT scan for Piao chickens and WL chickens, revealed that some Piao chicken tails were normal in number, and for a few Piao chickens tail length and tail bone numbers were between the rumpless and the normal tailed chickens. The results showed that the rumpless phenotype has not been completely fixed in Piao chicken breed. Using selection signature analysis and structural variation detection, we found a 4174 bp deletion located in the upstream region of IRX1 gene on chromosome 2 related to rumpless phenotype. Structural variation genotyping showed that the deletion was present in all 32 rumpless Piao chickens (del/del, wild/del) and absent from all 112 tailed chickens included in the dataset for the other 9 breeds and 2 tailed Piao chickens (wild/wild). In summary, all rumpless Piao chickens tested here carry this deletion mutation, to show a complete linkage association with rumplessness trait. We suggested that the 4174 bp deletion could be causative for rumpless phenotype in Piao chicken since this is the only mutation to show the complete linkage disequilibrium with rumplessness on whole genome level across all of 146 chickens from the 10 breeds. This study could facilitate a better understanding of the genetic characteristics of Piao chicken.

Genomic and transcriptomic analyses reveal genetic adaptation to cold conditions in the chickens.

Under the pressure of natural and artificial selection, domestic animals, including chickens, have evolved unique mechanisms of genetic adaptations such as high-altitude adaptation, hot and arid climate adaptation, and desert adaptation. Here, we investigated the genetic basis of cold tolerance in chicken by integrating whole-genome and transcriptome sequencing technologies. Genome-wide comparative analyses of 118 chickens living in different latitudes showed 46 genes and several pathways that may be involved in cold adaptation. The results of the functional enrichment analysis of differentially expressed genes proved the important role of metabolic pathways and immune-related pathways in cold tolerance in chickens. The subsequent integration of whole genome and transcriptome sequencing technology further identified six genes - dnah5 (dynein axonemal heavy chain 5), ptgs2 (prostaglandin-endoperoxide synthase 2), inhba (inhibin beta A subunit), irx2 (iroquois homeobox 2), ensgalg00000054917, and ensgalg00000046652 - requiring more detailed studies. In addition, we also discovered different allele frequency distributions of five SNPs (single nucleotide polymorphisms) within ptgs2 and nine SNPs within dnah5 in chickens in different latitudes, suggesting strong selective pressure of these two genes in chickens. We provide a novel insight into the genetic adaptation in chickens to cold environments, and provide a reference for evaluating and developing adaptive chicken breeds in cold environments.