Ovary Transcriptome Profiling in Broody and Egg-laying Chahua Chickens.

Broodiness in egg-laying hens (EHs) leads to ovarian atrophy, resulting in reduced egg-laying performance. However, the ovarian regulatory mechanisms in broody hens (BCs) remain elusive. Therefore, ovaries were removed from 300-day-old BCs and EHs for RNA sequencing. Ovarian morphology and histological characteristics of the BC and EH groups were compared and analyzed. The EH group had significantly more hierarchical follicles (HFs) and small yellow follicles (SYFs) than that of the BC group. Although several secondary follicles (SFs) and primary follicles were observed in the ovaries of the EH group, only a few SFs were observed in the ovaries of the BC group. Subsequently, RNA-sequencing analysis was conducted to determine the ovarian expression profiles of the two groups. Transcriptome sequencing identified 259 differentially expressed genes (DEGs) between the BC and EH groups. Of the 259 DEGs, 136 were upregulated and 123 were downregulated. The DEGs were mapped to 22 gene ontology terms and 4 Kyoto Encyclopedia of Genes and Genomes pathways for ovarian tissue. The analysis showed that matrix metalloproteinases 11/13 (MMP11/MMP13) were enriched in the extracellular matrix. The extracellular matrix mediated by MMP13 is affected by follicle-stimulating hormone, prolactin, and estrogen, which are critical signaling pathways that may affect ovarian follicle development to regulate the large yellow follicle reserve process and the ovulation cycle of broody Chahua chickens. These findings indicate that understanding differences in gene expression between the ovarian tissues of BCs and EHs could serve as a valuable reference point for enhancing egg-laying performance in Chahua chickens.

Study on chromatin regulation patterns of expression vectors in the PhiC31 integration site.

The PhiC31 integration system allows for targeted and efficient transgene integration and expression by recognizing pseudo attP sites in mammalian cells and integrating the exogenous genes into the open chromatin regions of active chromatin. In order to investigate the regulatory patterns of efficient gene expression in the open chromatin region of PhiC31 integration, this study utilized Ubiquitous Chromatin Opening Element (UCOE) and activating RNA (saRNA) to modulate the chromatin structure in the promoter region of the PhiC31 integration vector. The study analysed the effects of DNA methylation and nucleosome occupancy changes in the integrated promoter on gene expression levels. The results showed that for the OCT4 promoter with moderate CG density, DNA methylation had a smaller impact on expression compared to changes in nucleosome positioning near the transcription start site, which was crucial for enhancing downstream gene expression. On the other hand, for the SOX2 promoter with high CG density, increased methylation in the CpG island upstream of the transcription start site played a key role in affecting high expression, but the positioning and clustering of nucleosomes also had an important influence. In conclusion, analysing the DNA methylation patterns, nucleosome positioning, and quantity distribution of different promoters can determine whether the PhiC31 integration site possesses the potential to further enhance expression or overcome transgene silencing effects by utilizing chromatin regulatory elements.

Integrated analysis of microRNA and mRNA interactions in ovary of counter-season breeding and egg-ceased geese (Anser cygnoides).

Egg-ceasing is a phenomenon that occurs in most avian species and significantly reduces productivity. Although several factors are reported to regulate the reproduction progress, the underlying molecular mechanism of egg-ceasing remains obscure. Herein, we identified and explored the differentially expressed miRNAs and mRNAs involved in ovarian atrophy via high throughput sequencing. We identified a total of 901 mRNAs and 50 miRNAs that were differentially expressed in egg-laying and atrophic ovaries. Among them, numerous differentially expressed gene (DEG) transcripts and target genes for miRNAs were significantly enriched in Gene Ontology terms such as reproductive processes, cell proliferation, and apoptosis pathways. In addition, an interaction network was constructed by considering target relationships and correlation of the expression levels between ovary development-related genes, miRNAs and pathways. We discovered mRNA and miRNAs transcripts that are candidate regulators of ovary development in egg-ceased geese. Our findings expanded our understanding of the functional of miRNAs in ovarian atrophy and demonstrated that RNA-Seq is a powerful tool for examining the molecular mechanism in regulating egg-ceasing.

miR-122-5p regulates proliferation and apoptosis of chicken granulosa cells of hierarchal follicles by targeting MAPK3.

Chicken follicles plays a crucial role in the reproductive performance, especially in laying period. Recently, miR-122-5p has been found to be differentially expressed in the ovaries of rats with polycystic ovary syndrome and normal rats, indicating the potential role of miR-122-5p in the development of granulosa cells (GCs). In present study, we found that miR-122-5p was highly expressed in chicken atrophic ovaries. Herein, we investigated its function on GC proliferation and apoptosis of chicken in vitro. We found that overexpression of miR-122-5p significantly inhibited proliferation and promoted apoptosis of GCs, whereas the opposite effects were detected in miR-122-5p knockdown GCs. Meanwhile, mitogen-activated protein kinase 3 (MAPK3) was confirmed as a new target gene of miR-122-5p by bioinformatics software prediction and the dual-luciferase reporter assay verification. Furthermore, after knockdown of MAPK3, the function of MAPK3 for GC proliferation and apoptosis was opposite to that of miR-122-5p. Collectively, our results indicated that miR-122-5p impeded chicken GC proliferation and promoted apoptosis through the post-transcriptional downregulation of MAPK3.

Comparative Analysis of Skeletal Muscle DNA Methylation and Transcriptome of the Chicken Embryo at Different Developmental Stages.

DNA methylation is a key epigenetic mechanism involved in embryonic muscle development and plays an important role in early muscle development. In this study, we sought to investigate the effects of genome-wide DNA methylation by combining the expression profiles of the chicken embryonic muscle. Genome-wide DNA methylation maps and transcriptomes of muscle tissues collected from different embryonic development points (E7, E11, E17, and D1) were used for whole-genome bisulfite sequencing (WGBS) and RNA sequencing, respectively. We found that the differentially methylated genes (DMGs) were significantly associated with muscle organ development, regulation of skeletal muscle satellite cell proliferation, and actin filament depolymerization. Furthermore, genes TBX1, MEF2D, SPEG, CFL2, and TWF2 were strongly correlated with the methylation-caused expression switch. Therefore, we chose the CFL2 gene to explore its function in skeletal muscle satellite cells, and the in vitro experiments showed that CFL2 acts as a negative regulator of chicken skeletal muscle satellite cell proliferation and can induce cell apoptosis. These results provide valuable data for future genome and epigenome studies of chicken skeletal muscle and may help reveal the molecular mechanisms of potential economic traits.

miR-21-5p Regulates the Proliferation and Differentiation of Skeletal Muscle Satellite Cells by Targeting KLF3 in Chicken.

The proliferation and differentiation of skeletal muscle satellite cells (SMSCs) play an important role in the development of skeletal muscle. Our previous sequencing data showed that miR-21-5p is one of the most abundant miRNAs in chicken skeletal muscle. Therefore, in this study, the spatiotemporal expression of miR-21-5p and its effects on skeletal muscle development of chickens were explored using in vitro cultured SMSCs as a model. The results in this study showed that miR-21-5p was highly expressed in the skeletal muscle of chickens. The overexpression of miR-21-5p promoted the proliferation of SMSCs as evidenced by increased cell viability, increased cell number in the proliferative phase, and increased mRNA and protein expression of proliferation markers including PCNA, CDK2, and CCND1. Moreover, it was revealed that miR-21-5p promotes the formation of myotubes by modulating the expression of myogenic markers including MyoG, MyoD, and MyHC, whereas knockdown of miR-21-5p showed the opposite result. Gene prediction and dual fluorescence analysis confirmed that KLF3 was one of the direct target genes of miR-21-5p. We confirmed that, contrary to the function of miR-21-5p, KLF3 plays a negative role in the proliferation and differentiation of SMSCs. Si-KLF3 promotes cell number and proliferation activity, as well as the cell differentiation processes. Our results demonstrated that miR-21-5p promotes the proliferation and differentiation of SMSCs by targeting KLF3. Collectively, the results obtained in this study laid a foundation for exploring the mechanism through which miR-21-5p regulates SMSCs.

Whole-genome resequencing reveals loci with allelic transmission ratio distortion in F1 chicken population.

Allelic transmission ratio distortion (TRD) is the significant deviation from the expected ratio under Mendelian inheritance theory, which may be resulted from multiple disrupted biological processes, including germline selection, meiotic drive, gametic competition, imprint error, and embryo lethality. However, it is less known that whether or what extent the allelic TRD is present in farm animals. In this study, whole-genome resequencing technology was applied to reveal TRD loci in chicken by constructing a full-sib F1 hybrid population. Through the whole-genome resequencing data of two parents (30 ×) and 38 offspring (5 ×), we detected a total of 2850 TRD SNPs (p-adj < 0.05) located within 400 genes showing TRD, and all of them were unevenly distributed on macrochromosomes and microchromosomes. Our findings suggested that TRD in the chicken chromosome 16 might play an important role in chicken immunity and disease resistance and the MYH1F with significant TRD and allele-specific expression could play a key role in the fast muscle development. In addition, functional enrichment analyses revealed that many genes (e.g., TGFBR2, TGFBR3, NOTCH1, and NCOA1) with TRD were found in the significantly enriched biological process and InterPro terms in relation to embryonic lethality and germline selection. Our results suggested that TRD is considerably prevalent in the chicken genome and has functional implications.

Genome-wide analysis of spatiotemporal allele-specific expression in F1 hybrids of meat- and egg-type chickens.

In diploid organisms, each gene locus is composed of two parental alleles, which would interact with each other for determining the phenotypic variation. Better understanding of the allele-specific expression (ASE) in farm animals is much important to explore the genetic basis underlying economically important traits, which have been poorly understood yet. In this study, genome-wide analysis was applied to explore the spatiotemporal pattern of ASE in the F1 hybrids of chicken. First, meat- and egg-type chickens were selected for producing a full-sib F1 hybrid population (n = 57). Then, genome resequencing of two parents and 38 offspring were performed and liver and breast muscle samples (n = 38) were subjected to strand-specific RNA sequencing (ssRNA-seq) for ASE detection at 1, 28, and 56 days of age, respectively. The results accurately identified a total of 465 informative genes that could be distinguished with respect to their parental origins. There were 0.4% - 4.1% of informative genes showing ASE, and 57 of them were found across different tissues and time points. Besides, most ASE genes in chickens were tissue-specific, and no matter what the time-point pattern of one ASE gene, the same parental allele of this gene almost showed consistently higher or lower expression across all time points in the same type tissue. In conclusion, this study indicated that most of ASE genes were tissue-specific and time-dependent.

Genotype frequency distributions of 28 SNP markers in two commercial lines and five Chinese native chicken populations.

BACKGROUND: Modern breeding in the poultry industry mainly aims to produce high-performance poultry lines and breeds in two main directions of productivity, meat and eggs. To understand more about the productive potential of lowly selected Chinese native chicken populations, we selected 14 representative SNP markers strongly associated with growth traits or carcass traits and 14 SNP markers strongly associated with egg laying traits through previous reports. By using the MassArray technology, we detected the genotype frequency distributions of these 28 SNP markers in seven populations including four lowly selected as well as one moderately selected Sichuan native chicken populations, one commercial broiler line and one commercial layer line. RESULTS: Based on the genotype frequency distributions of these 28 SNP markers in 5 native chicken populations and 2 commercial lines, the results suggested that these Chinese indigenous chicken populations have a relatively close relationship with the commercial broiler line but a marked distinction from the commercial layer line. Two native chicken breeds, Shimian Caoke Chicken and Daheng Broilers, share similar genetic structure with the broiler line. CONCLUSIONS: Our observations may help us to better select and breed superior domestic chickens and provide new clues for further study of breeding programs in local chicken populations.

Analysis of Expression and Single Nucleotide Polymorphisms of INHA Gene Associated with Reproductive Traits in Chickens.

Inhibin α (INHA) is a candidate gene controlling ovulation in poultry. As the functional center of inhibin, INHA is a molecular marker associated with egg-laying performance. The objective of the current study was to analyze the expression differences of INHA in reproductive system and single nucleotide polymorphisms (SNPs) associations with reproductive traits in chickens. A total of 260 LuHua chickens (barred-feather chicken) were adopted. Twelve SNPs were detected in INHA gene. Among the exonic SNPs, three (g. 22177991A>G, g. 22178249G>C, and g. 22178414G>A) were missense mutations, resulting in the amino acid substitutions Val→Ala, Ala→Gly, and Ala→Gly, respectively. Four SNPs in the 3＇ untranslated region of INHA were predicted to either disturb or create microRNA-target interactions. Five SNPs (g. 22176870T>C, g. 22177100T>C, g. 22177149T>C, g. 22177991A>G, and g. 22178975G>A) were significantly associated with the number of eggs at 300 d of age (EN) (P < 0.05). Birds carrying GA genotype exhibited more EN than those with AA genotype (P < 0.01). In addition, quantitative real-time PCR revealed that INHA is mainly expressed in follicles on d 300 in chickens. Firstly, INHA expression increased and then decreased. The highest INHA mRNA abundance was found in the fifth largest preovulatory follicle (F5) (P < 0.01). In the prehierarchical follicles, INHA mRNA expression increased dramatically in small yellow follicles (SYF) (P < 0.01). Western blotting analysis showed that the INHA protein expression profile in the follicle was similar to its mRNA counterpart with greater expression in F5 and SYF follicles and lowest expression in F1 follicles (P < 0.05). These results suggest that INHA is a potential candidate gene improving reproductive traits in chickens.

Transcriptome analysis reveals differentially expressed genes and pathways for oviduct development and defense in prelaying and laying hens.

PROBLEM: The oviduct plays an indispensable role in the formation of eggs, especially the magnum and uterus. The identification of oviduct development in different stages will help to target candidate genes and pathways in regulation of albumen and eggshell formation, as well as defense mechanism in oviduct and egg. METHODS: To identify the function differences and the molecular defense mechanism of the oviduct and egg, we performed transcriptome sequencing analysis of the magnum and uterus in 120-d-old and 300-d-old Lohmann layers, three birds in each group. RESULTS: With fold changes (log2 ratio) ≥ 2 and false discovery rate (FDR) < 0.01, RNA-Seq revealed 1,040 genes expressed differentially in the magnum and 595 genes in the uterus. By combining GO enrichment and KEGG pathway analysis, it served to show that gene activities of the magnum and uterus in prelaying chickens were more likely to concentrate on growth and development, and after egg-laying, they were mainly inclined to enhance the substances transmembrane transport and secretion activities. We further characterized 1579 new genes, while only 803 of them were functionally annotated. A complex mixture of proteins related to defense was measured in this study. A subset of avian ß-defensins (AvBDs) and ovodefensins (OvoDs), that is, AvBD12, AvBD11, AvBD10, OvoDA1, OvoDB1, OvoDA2, OvoDA3, and OvoDBß, was detected to express in the magnum of laying hens at high levels. CONCLUSION: Collectively, the identification and functional analysis of these differentially expressed genes (DEGs), as well as specific expression of avian defensins, may contribute to understand the development and defense mechanisms of oviduct and eggs.

KLF5 functions in proliferation, differentiation, and apoptosis of chicken satellite cells.

KLF5 is an important regulator of cell proliferation, differentiation, and apoptosis in mammals. Little is known about the function of KLF5 in the regulation of chicken. Hence, qPCR was used to detect the expression of KLF5 in different tissues of chicken. And chicken skeletal muscle satellite cells (SMSCs) were transfected KLF5-specific small interfering RNA (siRNA) to assay SMSCs' proliferation, differentiation, and apoptosis. The results showed that KLF5 expressed higher in skeletal muscle than in the other tissues of chicken. Knockdown of KLF5 significantly inhibited the differentiation and increased apoptosis of chicken SMSCs, but it had no significant effect on proliferation of SMSCs. These results indicate that KLF5 plays an essential role during myogenesis, which will affect muscle repair and muscle regeneration, and may ameliorate muscle aging or sarcopenia.

Polymorphism in MC1R, TYR and ASIP genes in different colored feather chickens.

Coat color genetics successfully adapted and applied to different animal species, which provides a good demonstration of the concept of comparative genetics. In this study, we sequenced 945 bp fragments of melanocortin 1 receptor (MC1R) gene, 421 bp fragments of exon 1 of tyrosinase (TYR) gene and 266 bp fragments of exon 3 of agouti signaling protein (ASIP) gene for 250 individuals with five plumage color patterns. We detected a total of three SNPs (T398A, T637C, and G920C) in MC1R and built six haplotypes (H1-H6) based on the three SNPs. H5 and H6 haplotypes were mainly concentrated in white and grey chicken. And diplotypes H2H3 occurred in white feather and black-speckle feather with the same frequency. Moreover, a total of three SNPs (C47G, T120C, and T172C) in TYR were found and built six haplotypes (P1-P6) based on the three SNPs. Among them, haplotype P2, P3 and P6 were not occurred in black chicken, the diplotypes P1P6 and P4P6 were only distributed in white, gray and black-speckled feather. We only detected one SNP (T168C) in ASIP gene and found that genotype TT was advantage genotype in the different plumage color groups of chickens. Collectively, our study suggested an association between plumage color and genetic variation of MC1R, TYR and ASIP in chicken.

Molecular Characterization, Expression and Functional Analysis of Chicken STING.

Innate immunity is an essential line of defense against pathogen invasion which is gained at birth, and the mechanism involved is mainly to identify pathogen-associated molecular patterns through pattern recognition receptors. STING (stimulator of interferon genes) is a signal junction molecule that hosts the perception of viral nucleic acids and produces type I interferon response, which plays a crucial role in innate immunity. However, relatively few studies have investigated the molecular characterization, tissue distribution, and potential function of STING in chickens. In this study, we cloned the full-length cDNA of chicken STING that is composed of 1341 bp. Sequence analyses revealed that STING contains a 1140-bp open-reading frame that probably encodes a 379-amino acid protein. Multiple sequence alignments showed that the similarity of the chicken STING gene to other birds is higher than that of mammals. Real-time polymerase chain reaction (PCR) assays revealed that STING is highly expressed in the spleen, thymus and bursa of fabricious in chickens. Furthermore, we observed that STING expression was significantly upregulated both in vitro and in vivo following infection with Newcastle disease virus (NDV). STING expression was also significantly upregulated in chicken embryo fibroblasts upon stimulation with poly(I:C) or poly(dA:dT). Taken together, these findings suggest that STING plays an important role in antiviral signaling pathways in chickens.

Genotype frequency contributions of Mx1 gene in eight chicken breeds under different selection pressures.

Chicken Mx1 gene, as a positive antiviral gene, has been reported to provide resistance to several viruses especially avian influenza virus. In present research, the genotype frequency contributions of chicken Mx1 polymorphisms were characterized in five lowly selected as well as one moderately selected Sichuan native chicken populations and two highly selected commercial chicken breeds. Together with two newly-identified mutation sites (r.8A > G and r.1257T > C), a total of 13 single nucleotide polymorphisms (SNPs), including seven nonsynonymous mutation and six synonymous mutation, were found in the coding region of chicken Mx1 gene. Local Chinese chicken populations exhibited higher nucleotide diversity than commercial populations. Moreover, amino acid substitution sites as well as positive selection sites were located only in the domain not determined and GTPase domain, implying that amino acids mutations were likely needed in the modulatory and structural regions to better adapt the environment. Collectively, our results suggest that different selection pressures greatly influenced the genotype frequency contributions of chicken Mx1 gene. Understanding the interaction between genetic diversity and artificial selection may help us to better select and breed superior domestic chickens.

Molecular characterization of a novel ovodefensin gene in chickens.

Host defense peptides (HDPs) represent a large group of diverse small peptides that play important roles in host defense and disease resistance. In vertebrates, one of the main types of HDPs belong to defensins, which are less than 100 amino acid residues and characterized by a highly conserved motif of cysteine residues. Recently, a subfamily of defensins, namely ovodefensins (OvoDs), has been identified in birds and reptiles. However, both their family members and evolutionary relationships remain unclear. In the present study, we cloned and characterized a novel gene namely OvoDBß in chickens. Our results showed that the full length of chicken OvoDBß mRNA contains 344 bp nucleotides and encodes a 61-amino acid protein. We further revealed that the mRNA of OvoDBß is abundant in the oviduct of laying hens but absent in many other tissues. Additionally, sequences comparison and analyses suggested that OvoDBß is orthologous to the gene previously known as zebra finch OvoDB1, albeit it might exhibit specific structures. Furthermore, both OvoDBα and OvoDBß were existent in the genome of each bird, implying that two types of OvoDBs sharing same cysteine motif have already emerged before the species divergence. More importantly, recombinant OvoDBß mature peptide exerted antibacterial activity against Escherischia coli (CICC23657 strain) in vitro. These results collectively indicated that the putative sequence, namely chicken OvoDBß, is a function gene with potential antimicrobial property. Discovery and function characterization of novel HDP genes may help us develop novel antimicrobial agents in the future.

Genotypes of IFIH1 and IFIT5 in seven chicken breeds indicated artificial selection for commercial traits influenced antiviral genes.

Innate immunity is the first line against the invasion of pathogenic microorganisms. Previous reports only demonstrated production traits of commercial importance were often negatively correlated with innate disease resistance. However, whether different purpose of artificial selection influences innate immunity have not been understood. In this study, we cloned exon1, exon6 of IFIH1 and exon2 of IFIT5 by molecular biology techniques in seven different chicken breeds to detect the potential effect of artificial selection for commercial traits on disease resistance for the first time. In total, 8 single nucleotide polymorphisms(SNPs) of IFIH1 gene exon1 and exon6, 19 SNPs of IFIT5 gene exon2 were detected. We found all native chicken breeds had a relatively close relationship to broiler breeds but a remote relationship to layer breed. A great difference between CB and LLH with different selected purpose were observed. The allele frequencies of these two positive antiviral genes were associated with different purpose of artificial selection. Our experiment constituted the foundation for the interaction between commercial traits and immune trait.

Experimental Verification of CAPN1 and CAST Gene Polymorphisms in Different Generations of Da-Heng Broilers.

The micromolar calcium activated neutral protease (CAPN1) and calpastatin (CAST) have been widely regarded as genes related to muscle growth and meat tenderness. The objective of this study was to verify the association of SNPs of CAPN1 and CAST genes with carcass and tenderness traits and search the possible change patterns of SNPs in CAPN1 and CAST genes in six generations of broiler breeding process for growth rate, efficiency, and reproduction, during the third generation and the ninth generation, respectively. We found that, for CAPN1, genetic effects between SNPs (G3535A, C7198A) and meat tenderness were similar in different generations, while SNP3 (G7324A) was a novel polymorphism and had significant association with carcass and tenderness traits (P < 0.05) in this study. Furthermore, there was significant association between SNP4 (G9950A) and carcass indexes instead of tenderness traits (P < 0.05) which was consistent in the two generations. Moreover, although SNP6 (G37868A) of CAST had no relevance to carcass traits or tenderness traits in the third generation, it showed significant association with LW and CW in the ninth generation (P < 0.05).

The Relationship between MC1R Mutation and Plumage Color Variation in Pigeons.

The polymorphisms of MC1R gene play a crucial role in coat color variation in mammals; however, the relationship is still unclear in pigeons. In this study, we sequenced 741 bp fragment of the MC1R for 39 individuals with five plumage color patterns (gray plumage, n = 12; black plumage, n = 9; white plumage, n = 3; spotted plumage, n = 12; red plumage, n = 3). A total of three single nucleotide polymorphisms (SNPs) were detected, including G199A, G225A, and A466G, which subsequently determined four haplotypes (H1-H4). Among them, H1 is the predominant haplotype. Association analysis revealed that H1 and H3 were significantly associated with the black plumage trait (P < 0.05), while the H4 was significantly associated with gray plumage trait (P < 0.05). Furthermore, only diplotype H1H1 was significantly associated with black and gray traits of pigeons. Collectively, our study suggested an association between genetic variation of MC1R and plumage color in pigeon.

MyD88 Polymorphisms and Association with Susceptibility to Salmonella Pullorum.

Myeloid differentiation primary response gene 88 (MYD88), a universal adapter protein, plays an important role in activating the nuclear factor-κB (NF-κB) and regulating the expression of proinflammatory genes like tumor necrosis factor (TNF) and interleukin-1 (IL-1), which were highly involved in Salmonella Pullorum infection. To detect the relationship between polymorphisms of the MyD88 gene and Salmonella Pullorum disease, we screened the coding region (CDS) of the MYD88 gene by DNA pool construction and sequencing based on case-control study. Eight single nucleotide polymorphisms (SNPs) in the sequenced fragment (5 exons), 7 known loci and one novel mutation named G4810372T (SNP8), were found in the fifth exon. In addition, we found 7 nonsynonymous substitutions. The allele frequency of only one SNP, g.4810191C > T (SNP1), was significantly different (P < 0.05) between case and control groups. The genotype frequencies of SNP1 (g.4810191C > T) and SNP3 (g.4810257G > T) were of significant difference between the case and the control groups (P < 0.05). Collectively, SNPs of the MyD88 gene were significantly associated with susceptibility to Salmonella Pullorum infection, which can be used as a disease-resistant marker in chicken. These results provided a theoretical basis for future research on chicken breeding by marker-assisted selection.