Runs of homozygosity and selection signature analyses reveal putative genomic regions for artificial selection in layer breeding.

BACKGROUND: The breeding of layers emphasizes the continual selection of egg-related traits, such as egg production, egg quality and eggshell, which enhance their productivity and meet the demand of market. As the breeding process continued, the genomic homozygosity of layers gradually increased, resulting in the emergence of runs of homozygosity (ROH). Therefore, ROH analysis can be used in conjunction with other methods to detect selection signatures and identify candidate genes associated with various important traits in layer breeding. RESULTS: In this study, we generated whole-genome sequencing data from 686 hens in a Rhode Island Red population that had undergone fifteen consecutive generations of intensive artificial selection. We performed a genome-wide ROH analysis and utilized multiple methods to detect signatures of selection. A total of 141,720 ROH segments were discovered in whole population, and most of them (97.35%) were less than 3 Mb in length. Twenty-three ROH islands were identified, and they overlapped with some regions bearing selection signatures, which were detected by the De-correlated composite of multiple signals methods (DCMS). Sixty genes were discovered and functional annotation analysis revealed the possible roles of them in growth, development, immunity and signaling in layers. Additionally, two-tailed analyses including DCMS and ROH for 44 phenotypes of layers were conducted to find out the genomic differences between subgroups of top and bottom 10% phenotype of individuals. Combining the results of GWAS, we observed that regions significantly associated with traits also exhibited selection signatures between the high and low subgroups. We identified a region significantly associated with egg weight near the 25 Mb region of GGA 1, which exhibited selection signatures and has higher genomic homozygosity in the low egg weight subpopulation. This suggests that the region may be play a role in the decline in egg weight. CONCLUSIONS: In summary, through the combined analysis of ROH, selection signatures, and GWAS, we identified several genomic regions that associated with the production traits of layers, providing reference for the study of layer genome.

Temperature-induced embryonic diapause in chickens is mediated by PKC-NF-κB-IRF1 signaling.

BACKGROUND: Embryonic diapause (dormancy) is a state of temporary arrest of embryonic development that is triggered by unfavorable conditions and serves as an evolutionary strategy to ensure reproductive survival. Unlike maternally-controlled embryonic diapause in mammals, chicken embryonic diapause is critically dependent on the environmental temperature. However, the molecular control of diapause in avian species remains largely uncharacterized. In this study, we evaluated the dynamic transcriptomic and phosphoproteomic profiles of chicken embryos in pre-diapause, diapause, and reactivated states. RESULTS: Our data demonstrated a characteristic gene expression pattern in effects on cell survival-associated and stress response signaling pathways. Unlike mammalian diapause, mTOR signaling is not responsible for chicken diapause. However, cold stress responsive genes, such as IRF1, were identified as key regulators of diapause. Further in vitro investigation showed that cold stress-induced transcription of IRF1 was dependent on the PKC-NF-κB signaling pathway, providing a mechanism for proliferation arrest during diapause. Consistently, in vivo overexpression of IRF1 in diapause embryos blocked reactivation after restoration of developmental temperatures. CONCLUSIONS: We concluded that embryonic diapause in chicken is characterized by proliferation arrest, which is the same with other spices. However, chicken embryonic diapause is strictly correlated with the cold stress signal and mediated by PKC-NF-κB-IRF1 signaling, which distinguish chicken diapause from the mTOR based diapause in mammals.

De Novo Assembly of 20 Chicken Genomes Reveals the Undetectable Phenomenon for Thousands of Core Genes on Microchromosomes and Subtelomeric Regions.

The gene numbers and evolutionary rates of birds were assumed to be much lower than those of mammals, which is in sharp contrast to the huge species number and morphological diversity of birds. It is, therefore, necessary to construct a complete avian genome and analyze its evolution. We constructed a chicken pan-genome from 20 de novo assembled genomes with high sequencing depth, and identified 1,335 protein-coding genes and 3,011 long noncoding RNAs not found in GRCg6a. The majority of these novel genes were detected across most individuals of the examined transcriptomes but were seldomly measured in each of the DNA sequencing data regardless of Illumina or PacBio technology. Furthermore, different from previous pan-genome models, most of these novel genes were overrepresented on chromosomal subtelomeric regions and microchromosomes, surrounded by extremely high proportions of tandem repeats, which strongly blocks DNA sequencing. These hidden genes were proved to be shared by all chicken genomes, included many housekeeping genes, and enriched in immune pathways. Comparative genomics revealed the novel genes had 3-fold elevated substitution rates than known ones, updating the knowledge about evolutionary rates in birds. Our study provides a framework for constructing a better chicken genome, which will contribute toward the understanding of avian evolution and the improvement of poultry breeding.

Host genetics and gut microbiota jointly regulate blood biochemical indicators in chickens.

Blood biochemical indicators play a crucial role in assessing an individual's overall health status and metabolic function. In this study, we measured five blood biochemical indicators, including total cholesterol (CHOL), low-density lipoprotein cholesterol (LDL-CH), triglycerides (TG), high-density lipoprotein cholesterol (HDL-CH), and blood glucose (BG), as well as 19 growth traits of 206 male chickens. By integrating host whole-genome information and 16S rRNA sequencing of the duodenum, jejunum, ileum, cecum, and feces microbiota, we assessed the contributions of host genetics and gut microbiota to blood biochemical indicators and their interrelationships. Our results demonstrated significant negative phenotypic and genetic correlations (r = - 0.20 ~ - 0.67) between CHOL and LDL-CH with growth traits such as body weight, abdominal fat content, muscle content, and shin circumference. The results of heritability and microbiability indicated that blood biochemical indicators were jointly regulated by host genetics and gut microbiota. Notably, the heritability of HDL-CH was estimated to be 0.24, while the jejunal microbiability for BG and TG reached 0.45 and 0.23. Furthermore, by conducting genome-wide association study (GWAS) with the single-nucleotide polymorphism (SNPs), insertion/deletion (indels), and structural variation (SV), we identified RAP2C, member of the RAS oncogene family (RAP2C), dedicator of cytokinesis 11 (DOCK11), neurotensin (NTS) and BOP1 ribosomal biogenesis factor (BOP1) as regulators of HDL-CH, and glycerophosphodiester phosphodiesterase domain containing 5 (GDPD5), dihydrodiol dehydrogenase (DHDH), and potassium voltage-gated channel interacting protein 1 (KCNIP1) as candidate genes of BG. Moreover, our findings suggest that cecal RF39 and Clostridia_UCG_014 may be linked to the regulation of CHOL, and jejunal Streptococcaceae may be involved in the regulation of TG. Additionally, microbial GWAS results indicated that the presence of gut microbiota was under host genetic regulation. Our findings provide valuable insights into the complex interaction between host genetics and microbiota in shaping the blood biochemical profile of chickens. KEY POINTS: • Multiple candidate genes were identified for the regulation of CHOL, HDL-CH, and BG. • RF39, Clostridia_UCG_014, and Streptococcaceae were implicated in CHOL and TG modulation. • The composition of gut microbiota is influenced by host genetics.

Mechanisms of hepatic steatosis in chickens: integrated analysis of the host genome, molecular phenomics and gut microbiome.

Hepatic steatosis is the initial manifestation of abnormal liver functions and often leads to liver diseases such as nonalcoholic fatty liver disease in humans and fatty liver syndrome in animals. In this study, we conducted a comprehensive analysis of a large chicken population consisting of 705 adult hens by combining host genome resequencing; liver transcriptome, proteome, and metabolome analysis; and microbial 16S ribosomal RNA gene sequencing of each gut segment. The results showed the heritability (h2 = 0.25) and duodenal microbiability (m2 = 0.26) of hepatic steatosis were relatively high, indicating a large effect of host genetics and duodenal microbiota on chicken hepatic steatosis. Individuals with hepatic steatosis had low microbiota diversity and a decreased genetic potential to process triglyceride output from hepatocytes, fatty acid ß-oxidation activity, and resistance to fatty acid peroxidation. Furthermore, we revealed a molecular network linking host genomic variants (GGA6: 5.59-5.69 Mb), hepatic gene/protein expression (PEMT, phosphatidyl-ethanolamine N-methyltransferase), metabolite abundances (folate, S-adenosylmethionine, homocysteine, phosphatidyl-ethanolamine, and phosphatidylcholine), and duodenal microbes (genus Lactobacillus) to hepatic steatosis, which could provide new insights into the regulatory mechanism of fatty liver development.

Hepatic steatosis is associated with dysregulated cholesterol metabolism and altered protein acetylation dynamics in chickens.

BACKGROUND: Hepatic steatosis is a prevalent manifestation of fatty liver, that has detrimental effect on the health and productivity of laying hens, resulting in economic losses to the poultry industry. Here, we aimed to systematically investigate the genetic regulatory mechanisms of hepatic steatosis in laying hens. METHODS: Ninety individuals with the most prominent characteristics were selected from 686 laying hens according to the accumulation of lipid droplets in the liver, and were graded into three groups, including the control, mild hepatic steatosis and severe hepatic steatosis groups. A combination of transcriptome, proteome, acetylome and lipidome analyses, along with bioinformatics analysis were used to screen the key biological processes, modifications and lipids associated with hepatic steatosis. RESULTS: The rationality of the hepatic steatosis grouping was verified through liver biochemical assays and RNA-seq. Hepatic steatosis was characterized by increased lipid deposition and multiple metabolic abnormalities. Integration of proteome and acetylome revealed that differentially expressed proteins (DEPs) interacted with differentially acetylated proteins (DAPs) and were involved in maintaining the metabolic balance in the liver. Acetylation alterations mainly occurred in the progression from mild to severe hepatic steatosis, i.e., the enzymes in the fatty acid oxidation and bile acid synthesis pathways were significantly less acetylated in severe hepatic steatosis group than that in mild group (P < 0.05). Lipidomics detected a variety of sphingolipids (SPs) and glycerophospholipids (GPs) were negatively correlated with hepatic steatosis (r ≤ -0.5, P < 0.05). Furthermore, the severity of hepatic steatosis was associated with a decrease in cholesterol and bile acid synthesis and an increase in exogenous cholesterol transport. CONCLUSIONS: In addition to acquiring a global and thorough picture of hepatic steatosis in laying hens, we were able to reveal the role of acetylation in hepatic steatosis and depict the changes in hepatic cholesterol metabolism. The findings provides a wealth of information to facilitate a deeper understanding of the pathophysiology of fatty liver and contributes to the development of therapeutic strategies.

Genetic and microbiome analysis of feed efficiency in laying hens.

Improving feed efficiency is an important target for poultry breeding. Feed efficiency is affected by host genetics and the gut microbiota, but many of the mechanisms remain elusive in laying hens, especially in the late laying period. In this study, we measured feed intake, body weight, and egg mass of 714 hens from a pedigreed line from 69 to 72 wk of age and calculated the residual feed intake (RFI) and feed conversion ratio (FCR). In addition, fecal samples were also collected for 16S ribosomal RNA gene sequencing (V4 region). Genetic analysis was then conducted in DMU packages by using AI-REML with animal model. Moderate heritability estimates for FCR (h2 = 0.31) and RFI (h2 = 0.52) were observed, suggesting that proper selection programs can directly improve feed efficiency. Genetically, RFI was less correlated with body weight and egg mass than that of FCR. The phenotypic variance explained by gut microbial variance is defined as the microbiability (m2). The microbiability estimates for FCR (m2 = 0.03) and RFI (m2 = 0.16) suggested the gut microbiota was also involved in the regulation of feed efficiency. In addition, our results showed that the effect of host genetics on fecal microbiota was minor in three aspects: 1) microbial diversity indexes had low heritability estimates, and genera with heritability estimates more than 0.1 accounted for only 1.07% of the tested fecal microbiota; 2) the genetic relationship correlations between host genetics and different microbial distance were very weak, ranging from -0.0057 to -0.0003; 3) the microbial distance between different kinships showed no significant difference. Since the RFI has the highest microbiability, we further screened out three genera, including Anaerosporobacter, Candidatus Stoquefichus, and Fournierella, which were negatively correlated with RFI and played positive roles in improving the feed efficiency. These findings contribute to a great understanding of the genetic background and microbial influences on feed efficiency.

Microbial composition of egg component and its association with hatchability of laying hens.

The internal quality of eggs is critical for human consumption and embryonic development. However, microorganisms inside eggs have not been thoroughly investigated for their roles in determining the egg's internal quality. Here, a total of 21 hens were selected from more than 1,000 chickens based on their hatching results and were divided into high- and low-hatchability groups. Then, we collected 72 eggs from these 21 hens to obtain egg whites and yolks, including 54 fresh eggs and 18 eggs after 12 days of incubation. We characterized the microbial composition of egg yolks and whites, the microbial change along incubation, and differences in microbial abundance between the high- and low-hatchability groups. The results indicated that egg whites are not sterile. Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes were the dominant phyla in egg yolk and white. There was a large difference in the microbial composition between egg whites and yolks, and this difference increased after 12 days of incubation. Egg whites have lower microbial diversity than egg yolks owing to the presence of antibacterial substances such as lysozyme in the egg white. After a 12-day incubation, the microbial diversity decreased in egg whites but increased slightly in egg yolks. Meanwhile, the microbes in egg white can migrate to egg yolk during incubation. Additionally, Genus Muribaculaceae was identified as a biomarker in egg yolks incubated for 12 days and was more often detected in healthy groups. On the contrary, more genus Rothia were found in the fresh egg yolk of the low hatchability groups and was considered to have low virulence. These findings shed light on the composition and differences in microbiota between egg yolks and whites and may open new avenues for studying embryonic development in chickens.

Microbiota continuum along the chicken oviduct and its association with host genetics and egg formation.

The microbiota of female reproductive tract have attracted considerable attention in recent years due to their effects on host fitness. However, the microbiota throughout the chicken oviduct and its symbiotic relationships with the host have not been well characterized. Here, we characterized the microbial composition of six segments of the reproductive tract, including the infundibulum, magnum, isthmus, uterus, vagina and cloaca, in pedigreed laying hens with phenotypes of egg quality and quantity. We found that the microbial diversity gradually increased along the reproductive tract from the infundibulum to the cloaca, and the microbial communities were distinct among the cloaca, vagina and four other oviductal segments. The magnum exhibited the lowest diversity, given that the lysozyme and other antimicrobial proteins are secreted at this location. The results of correlation estimated showed that the relationship between host genetic kinship and microbial distance was negligible. Additionally, the genetically related pairwise individuals did not exhibit a more similar microbial community than unrelated pairs. Although the egg might be directly contaminated with potential pathogenic bacteria during egg formation and oviposition, some microorganisms provide long-term benefits to the host. Among these, we observed that increased abundance of vaginal Staphylococcus and Ralstonia was significantly associated with darker eggshells. Meanwhile, vaginal Romboutsia could be used as a predictor for egg number. These findings provide insight into the nature of the chicken reproductive tract microbiota and highlight the effect of oviductal bacteria on the process of egg formation.

The Spatial and Temporal Characterization of Gut Microbiota in Broilers.

The gut microbiota of chickens plays an important role in host physiology. However, the colonization and prevalence of gut microbiota have not been well-characterized. Here, we performed 16S rRNA gene sequencing on the duodenal, cecal and fecal microbiota of broilers at 1, 7, 21, and 35 days of age and characterized the dynamic succession of microbiota across the intestinal tract. Our results showed that Firmicutes was the most abundant phylum detected in each gut site at various ages, while the microbial diversity and composition varied among the duodenum, cecum, and feces at different ages. The microbial diversity and complexity of the cecal microbiota increased with age, gradually achieving stability at 21 days of age. As a specific genus in the cecum, Clostridium_sensu_stricto_1 accounted for 83.50% of the total abundance at 1 day of age, but its relative abundance diminished with age. Regarding the feces, the highest alpha diversity was observed at 1 day of age, significantly separated from the alpha diversity of other ages. In addition, no significant differences were observed in the alpha diversity of duodenal samples among 7, 21, and 35 days of age. The predominant bacterium, Lactobacillus, was relatively low (0.68-6.04%) in the intestinal tract of 1-day-old chicks, whereas its abundance increased substantially at 7 days of age and was higher in the duodenum and feces. Escherichia-Shigella, another predominant bacterium in the chicken intestinal tract, was also found to be highly abundant in fecal samples, and the age-associated dynamic trend coincided with that of Lactobacillus. In addition, several genera, including Blautia, Ruminiclostridium_5, Ruminococcaceae_UCG-014, and [Ruminococcus]_torques_group, which are related to the production of short-chain fatty acids, were identified as biomarker bacteria of the cecum after 21 days of age. These findings shed direct light on the temporal and spatial dynamics of intestinal microbiota and provide new opportunities for the improvement of poultry health and production.