The Creeper trait, a classical monogenic phenotype of chicken, is controlled by a dominant semi-lethal gene. This trait has been widely cited in the genetics and molecular biology textbooks for illustrating autosomal dominant semi-lethal inheritance over decades. However, the genetic basis of the Creeper trait remains unknown. Here we have utilized ultra-deep sequencing and extensive analysis for targeting causative mutation controlling the Creeper trait. Our results indicated that the deletion of Indian hedgehog (IHH) gene was only found in the whole-genome sequencing data of lethal embryos and Creeper chickens. Large scale segregation analysis demonstrated that the deletion of IHH was fully linked with early embryonic death and the Creeper trait. Expression analysis showed a much lower expression of IHH in Creeper than wild-type chickens. We therefore suggest the deletion of IHH to be the causative mutation for the Creeper trait in chicken. Our findings unravel the genetic basis of the longstanding Creeper phenotype mystery in chicken as the same gene also underlies bone dysplasia in human and mouse, and thus highlight the significance of IHH in animal development and human haploinsufficiency disorders.

Adipose tissues have a central role in organisms, and adipose content is a crucial economic trait of poultry. Pekin duck is an ideal model to study the mechanism of abdominal and subcutaneous adipose deposition for its high ability of adipose synthesis and deposition. Alternative splicing contributes to functional diversity in abdominal and subcutaneous adipose. However, there has been no systematic analysis of the dynamics of differential alternative splicing of abdominal and subcutaneous adipose in Pekin duck. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to explore the transcriptional complexity of abdominal and subcutaneous adipose in Pekin ducks. In total, 143,931 and 111,337 full-length non-chimeric transcriptome sequences of abdominal and subcutaneous adipocytes were obtained from 41.78 GB raw data, respectively. These data led us to identify 19,212 long non-coding RNAs (lncRNAs) and 74,571 alternative splicing events. In addition, combined with the next-generation sequencing technology, we correlated the structure and function annotation with the differential expression profiles of abdominal and subcutaneous adipose transcripts. This study identified lots of novel alternative splicing events and major transcripts of transcription factors related to adipose synthesis. STAT3 was reported as a vital gene for adipogenesis, and we found that its major transcript is STAT3-1, which may play a considerable role in the process of adipose synthesis in Pekin duck. This study greatly increases our understanding of the gene models, genome annotations, genome structures, and the complexity and diversity of abdominal and subcutaneous adipose in Pekin duck. These data provide insights into the regulation of alternative splicing events, which form an essential part of transcript diversity during adipogenesis in poultry. The results of this study provide an invaluable resource for studying alternative splicing and tissue-specific expression.

A considerable number of muscle development-related genes were differentially expressed in the early stage of avian adipocyte differentiation. However, the functions of them in adipocyte differentiation remain largely known. In this study, the myoblast determination protein 1 (MYOD1) was selected as a representative of muscle development. We investigated its expression, function, and regulation in avian adipocyte differentiation.

Fat deposition is very important in pig production, and its mechanism is not clearly understood. MicroRNAs (miRNAs) play critical roles in fat deposition and energy metabolism. In the current study, we investigated the mRNA and miRNA transcriptome in the livers of Landrace pigs with extreme backfat thickness to explore miRNA-mRNA regulatory networks related to lipid deposition and metabolism. A comparative analysis of liver mRNA and miRNA transcriptomes from pigs (four pigs per group) with extreme backfat thickness was performed. We identified differentially expressed genes from RNA-seq data using a Cufflinks pipeline. Seventy-one differentially expressed genes (DEGs), including twenty-eight well annotated on the porcine reference genome genes, were found. The upregulation genes in pigs with higher backfat thickness were mainly involved in fatty acid synthesis, and included fatty acid synthase (FASN), glucokinase (GCK), phosphoglycerate dehydrogenase (PHGDH), and apolipoprotein A4 (APOA4). Cytochrome P450, family 2, subfamily J, polypeptide 34 (CYP2J34) was lower expressed in pigs with high backfat thickness, and is involved in the oxidation of arachidonic acid. Moreover, 13 differentially expressed miRNAs were identified. Seven miRNAs were associated with fatty acid synthesis, lipid metabolism, and adipogenic differentiation. Based on comprehensive analysis of the transcriptome of both mRNAs and miRNAs, an important regulatory network, in which six DEGs could be regulated by differentially expressed miRNAs, was established for fat deposition. The negative correlate in the regulatory network including, miR-545-5p and GRAMD3, miR-338 and FASN, and miR-127, miR-146b, miR-34c, miR-144 and THBS1 indicate that direct suppressive regulation may be involved in lipid deposition and energy metabolism. Based on liver mRNA and miRNA transcriptomes from pigs with extreme backfat thickness, we identified 28 differentially expressed genes and 13 differentially expressed miRNAs, and established an important miRNA-mRNA regulatory network. This study provides new insights into the molecular mechanisms that determine fat deposition in pigs.

A lack of the complete pig proteome has left a gap in our knowledge of the pig genome and has restricted the feasibility of using pigs as a biomedical model. In this study, we developed a tissue-based proteome map using 34 major normal pig tissues. A total of 5841 unknown protein isoforms were identified and systematically characterized, including 2225 novel protein isoforms, 669 protein isoforms from 460 genes symbolized beginning with LOC, and 2947 protein isoforms without clear NCBI annotation in the current pig reference genome. These newly identified protein isoforms were functionally annotated through profiling the pig transcriptome with high-throughput RNA sequencing of the same pig tissues, further improving the genome annotation of the corresponding protein-coding genes. Combining the well-annotated genes that have parallel expression pattern and subcellular witness, we predicted the tissue-related subcellularlocations and potential functions for these unknown proteins. Finally, we mined 3081 orthologous genes for 52.7% of unknown protein isoforms across multiple species, referring to 68 KEGG pathways as well as 23 disease signaling pathways. These findings provide valuable insights and a rich resource for enhancing studies of pig genomics and biology, as well as biomedical model application to human medicine.

The endometrium of sheep consists of plenty of raised aglandular areas called caruncular (C), and intensely glandular intercaruncular areas (IC). In order to better understand the endometrium involved mechanisms of implantation, we used LC-MS/MS technique to profile the proteome of ovine endometrial C areas and IC areas separately during the peri-implantation period, and then compared the proteomic profiles between these two areas. We successfully detected 1740 and 1813 proteins in C areas and IC areas respectively. By comparing the proteome of these two areas, we found 170 differentially expressed proteins (DEPs) (P < 0.05), functional bioinformatics analysis showed these DEPs were mainly involved in growth and remodeling of endometrial tissue, cell adhesion and protein transport, and so on. Our study, for the first time, provided a proteomic reference for elucidating the differences between C and IC areas, as an integrated function unit respectively, during the peri-implantation period. The results could help us to better understand the implantation in the ewes. In addition, we established a relatively detailed protein database of ovine endometrium, which provide a unique reference for further studies.

To select metabolic biomarkers and differentially expressed genes (DEGs) associated with resistant-ascites syndrome (resistant-AS), we used innovative techniques such as metabolomics and transcriptomics to comparatively examine resistant-AS chickens and AS controls. Metabolomic evaluation of chicken serum using ultra-performance liquid chromatography-quadruple time-of-flight high-sensitivity mass spectrometry (UPLC-QTOF/HSMS) showed significantly altered lysoPC(18:1), PE(18:3/16:0), PC(20:1/18:3), DG(24:1/22:6/0:0), PS(18:2/18:0), PI(16:0/16:0), PS(18:0/18:1), PS(14:1/14:0), dihydroxyacetone, ursodeoxycholic acid, tryptophan, L-valine, cycloserine, hypoxanthine, and 4-O-Methylmelleolide concentrations on day 21 and LysoPC(18:0), LysoPE(20:1/0:0), LysoPC(16:0), LysoPE(16:0/0:0), hypoxanthine, dihydroxyacetone, 4-O-Methylmelleolide, LysoPC(18:2), and PC(14:1/22:1) concentrations on day 35, between the susceptible and resistant groups. Compared to the susceptible group, transcriptomic analysis of liver samples using RNA-seq revealed 413 DEGs on day 21 and 214 DEGs on day 35 in the resistant group. Additional evaluations using gene ontology (GO) indicate that significant enrichment occurred in the oxygen transportation, defensive reactions, and protein modifications of the decreased DEGs as well as in the cell morphological formation, neural development, and transforming growth factor (TGF)-beta signalling of the increased DEGs on day 21. Oxygen transportation was also significantly enriched for downregulated DEGs on day 35. The combinatory evaluation of the metabolome and the transcriptome suggests the possible involvement of glycerophospholipid metabolism in the development of resistant-AS in broilers.

The reproductive system of a female bird is responsible for egg production. The genes highly expressed in oviduct are potentially important. From RNA-seq analysis, C2H9orf152 (an orthologous gene of human C9orf152) was identified as highly expressed in chicken uterus. To infer its function, we obtained and characterized its complete cDNA sequence, determined its spatiotemporal expression, and probed its transcription factor(s) through pharmaceutical approach. Data showed that the complete cDNA sequence was 1468bp long with a 789bp of open reading frame. Compared to other tested tissues, this gene was highly expressed in the oviduct and liver tissues, especially uterus. Its expression in uterus was gradually increased during developmental and reproductive periods, which verified its involvement in the growth and maturity of reproductive system. In contrast, its expression was not significant different between active and quiescent uterus, suggesting the role of C2H9orf152 in reproduction is likely due to its long-term effect. Moreover, based on its 5'-flanking sequence, Foxd3 and Hnf4a were predicted as transcription factors of C2H9orf152. Using berberine or retinoic acid (which can regulate the activities of Hnf4a and Foxd3, respectively), we demonstrated suppression of C2H9orf152 by the chemicals in chicken primary hepatocytes. As retinoic acid regulates calcium metabolism, and Hnf4a is a key nuclear factor to liver, these findings suggest that C2H9orf152 is involved in liver function and calcium metabolism of reproductive system. In conclusion, C2H9orf152 may have a long-term effect on chicken reproductive system by regulating calcium metabolism, suggesting this gene has an important implication in the improvement of egg production and eggshell quality.

Japanese quail (Coturnix japonica), a recently domesticated poultry species, is important not only as an agricultural product, but also as a model bird species for genetic research. However, most of the biological questions concerning genomics, phylogenetics, and genetics of some important economic traits have not been answered. It is thus necessary to complete a high-quality genome sequence as well as a series of comparative genomics, evolution, and functional studies.

Selection for cold tolerance in chickens is important for improving production performance and animal welfare. The identification of chicken breeds with higher cold tolerance and production performance will help to target candidates for the selection. The thyroid gland plays important roles in thermal adaptation, and its function is influenced by breed differences and transcriptional plasticity, both of which remain largely unknown in the chicken thyroid transcriptome. In this study, we subjected Bashang Long-tail (BS) and Rhode Island Red (RIR) chickens to either cold or warm environments for 21 weeks and investigated egg production performance, body weight changes, serum thyroid hormone concentrations, and thyroid gland transcriptome profiles. RIR chickens had higher egg production than BS chickens under warm conditions, but BS chickens produced more eggs than RIRs under cold conditions. Furthermore, BS chickens showed stable body weight gain under cold conditions while RIRs did not. These results suggested that BS breed is a preferable candidate for cold-tolerance selection and that the cold adaptability of RIRs should be improved in the future. BS chickens had higher serum thyroid hormone concentrations than RIRs under both environments. RNA-Seq generated 344.3 million paired-end reads from 16 sequencing libraries, and about 90% of the processed reads were concordantly mapped to the chicken reference genome. Differential expression analysis identified 46-1,211 genes in the respective comparisons. With regard to breed differences in the thyroid transcriptome, BS chickens showed higher cell replication and development, and immune response-related activity, while RIR chickens showed higher carbohydrate and protein metabolism activity. The cold environment reduced breed differences in the thyroid transcriptome compared with the warm environment. Transcriptional plasticity analysis revealed different adaptive responses in BS and RIR chickens to cope with the cold, and showed higher responsiveness in BS compared with RIR chickens, suggesting greater adaptability of the thyroid in BS chickens. Moreover, 10,053 differential splicing events were revealed among the groups, with RNA splicing and processing, gene expression, transport, and metabolism being the main affected biological processes, identifying a valuable alternative splicing repertoire for the chicken thyroid. A short isoform of TPO (encoding thyroid peroxidase) containing multiple open reading frames was generated in both breeds by skipping exons 4 and 5 in the cold environment. These findings provide novel clues for future studies of the molecular mechanisms underlying cold adaptation and/or acclimation in chickens.

The fat deposition has a crucial role in animal meat flavor, and fat deposition-related traits are vital for breeding in the commercial duck industry. Avian fat-related traits are typical complex phenotypes, which need a large amount of data to analyze the genetic loci.

Mitochondrial remodeling during the peri-implantation stage is the hallmark event essential for normal embryogenesis. Among the changes, enhanced oxidative phosphorylation is critical for supporting high energy demands of postimplantation embryos, but increases mitochondrial oxidative stress, which in turn threatens mitochondrial DNA (mtDNA) stability. However, how mitochondria protect their own histone-lacking mtDNA, during this stage remains unclear. Concurrently, the mitochondrial genome gain DNA methylation by this stage. Its spatiotemporal coincidence with enhanced mitochondrial stress led us to ask if mtDNA methylation has a role in maintaining mitochondrial genome stability. Herein, we report that mitochondrial genome undergoes de novo mtDNA methylation that can protect mtDNA against enhanced oxidative damage during the peri-implantation window. Mitochondrial genome gains extensive mtDNA methylation during transition from blastocysts to postimplantation embryos, thus establishing relatively hypermethylated mtDNA from hypomethylated state in blastocysts. Mechanistic study revealed that DNA methyltransferase 3A (DNMT3A) and DNMT3B enter mitochondria during this process and bind to mtDNA, via their unique mitochondrial targeting sequences. Importantly, loss- and gain-of-function analyses indicated that DNMT3A and DNMT3B are responsible for catalyzing de novo mtDNA methylation, in a synergistic manner. Finally, we proved, in vivo and in vitro, that increased mtDNA methylation functions to protect mitochondrial genome against mtDNA damage induced by increased mitochondrial oxidative stress. Together, we reveal mtDNA methylation dynamics and its underlying mechanism during the critical developmental window. We also provide the functional link between mitochondrial epigenetic remodeling and metabolic changes, which reveals a role for nuclear-mitochondrial crosstalk in establishing mitoepigenetics and maintaining mitochondrial homeostasis.

The composition of the gut microbiome plays important roles in digestion, nutrient absorption, and health. Here, we analyzed the microbial composition in the duodenum and ileum of yellow broilers. Chickens were grouped based on feed efficiency (high feed efficiency [HFE] and low feed efficiency [LFE] groups; n = 22 each). Microbial samples from the duodenum and ileum were collected, and 16S rRNA sequencing of the V3-V4 region was performed. The dominant bacteria in the duodenum were from the phyla Firmicutes and Cyanobacteria and the genera Lactobacillus, Faecalibacterium, and Ruminococcus. In the ileum, the phyla Firmicutes and Proteobacteria and the genera Lactobacillus, SMB53 and Enterococcus were predominant. Alpha diversity analysis showed that the microbiota diversity was significantly higher in the duodenum than in the ileum. The structure of the ileal microbiota was similar between groups, and the species richness of the microbiota in the HFE group was significantly higher than that in the LFE group. In the HFE and LFE groups, Firmicutes and Cyanobacteria were negatively correlated, and Lactobacillus had medium to high negative correlations with most other genera. Functional prediction analysis showed that the gluconeogenesis I pathway was the most abundant differential metabolic pathway and was significantly altered in the LFE group. Moreover, although the microbial community structures were similar in the duodenum and ileum, the diversity of the microbial community was significantly higher in the duodenum than in the ileum. Pearson correlation analysis revealed that the phylum Chloroflexi and genera Acinetobacter, Pseudomonas, Bacillus and Neisseria were with coefficients <-0.3 or >0.3. In the ileum, Ruminococcus may be associated with HFE whereas Faecalibacterium may be associated with LFE. These findings may provide valuable foundations for future research on composition and diversity of intestinal microbes and provide insights into the roles of intestinal microbes in improving feed efficiency and the industrial economic benefits of yellow broilers.

Chicken is considered an ideal model species to study the synthesis of polyunsaturated fatty acids (PUFAs) due to its appropriate proportions of fatty acids and abundant content of PUFAs, suitable for human consumption. However, the molecular mechanisms regulating poultry PUFA synthesis remain unclear. Here, we systematically explored the transcriptional regulation activity of the gene family related to PUFA synthesis in chicken by carrying out the Dual-Luciferase Reporter Assay. We identified the core promoter regions of members of the chicken PUFA synthesis-related gene family, including ELOVL1, ELOVL2, ELOVL3, ELOVL4, ELOVL5, ELOVL6, ELOVL7, FADS1, FADS2, FADS6, SCD, and SCD5. Additionally, changes in relative fluorescence values of different truncated segments in the upstream regulatory region of these genes indicate the existence of regulatory regions. Furthermore, we predicted the transcription factors that bind to the identified core promoter regions of multiple genes, including Sp1, NF-1, C/EBPalpha, etc. These findings provide a basis for the molecular mechanisms regulating poultry PUFA synthesis and offer new scientific insight into the potential improvement of poultry meat quality in the future.

Fat deposition is highly correlated with the growth, meat quality, reproductive performance and immunity of pigs. Fatty acid synthesis takes place mainly in the adipose tissue of pigs; therefore, in this study, a high-throughput massively parallel sequencing approach was used to generate adipose tissue transcriptomes from two groups of Songliao black pigs that had opposite backfat thickness phenotypes. The total number of paired-end reads produced for each sample was in the range of 39.29-49.36 millions. Approximately 188 genes were differentially expressed in adipose tissue and were enriched for metabolic processes, such as fatty acid biosynthesis, lipid synthesis, metabolism of fatty acids, etinol, caffeine and arachidonic acid and immunity. Additionally, many genetic variations were detected between the two groups through pooled whole-genome resequencing. Integration of transcriptome and whole-genome resequencing data revealed important genomic variations among the differentially expressed genes for fat deposition, for example, the lipogenic genes. Further studies are required to investigate the roles of candidate genes in fat deposition to improve pig breeding programs.

The chorioallantoic placenta is a shared derived feature of "placental" mammals essential for the success of eutherian reproduction. Identifying the genes involved in the emergence of the placenta may provide clues for understanding the biology of this organ. Here we identify among 4960 single copy genes in mammals, 222 that show high expression levels in human placentas at term. Further, we present evidence that 94 of these 222 genes evolved adaptively during human evolutionary history since the time of the last common ancestor of eutherian mammals. Remarkably, the majority of positive selection occurred on the eutherian stem lineage suggesting that ancient adaptations have been retained in the human placenta. Of these positively selected genes, 28 have been shown to play a role in human pregnancy and placental biology, and at least 26 have important pregnancy-related phenotypes in mice. Adaptations in genes highly expressed in human placenta are attractive candidates for functional and clinical studies.

Pekin duck is an important animal model for its ability for fat synthesis and deposition. However, transcriptional dynamic regulation of adipose differentiation driven by complex signal cascades remains largely unexplored in this model. This study aimed to explore adipogenic transcriptional dynamics before (proliferation) and after (differentiation) initial preadipocyte differentiation in ducks.

Early pregnancy failure has a profound impact on both human reproductive health and animal production. 2/3 pregnancy failures occur during the peri-implantation period; however, the underlying mechanism(s) remains unclear. Well-organized modification of the endometrium to a receptive state is critical to establish pregnancy. Aberrant endometrial modification during implantation is thought to be largely responsible for early pregnancy loss.

Fatness traits in animals are important for their growth, meat quality, reproductive performance, and immunity. The liver is the principal organ of the regulation of lipid metabolism, and this study used massive parallelized high-throughput sequencing technologies to determine the porcine liver tissue transcriptome architecture of two full-sibling Songliao black pigs harboring extremely different phenotypes of backfat thickness.