Polymorphisms in Wnt signaling pathway genes are significantly associated with chicken carcass traits

Association of variants and expression levels of MYOD1 gene with carcass and muscle characteristic traits in domestic pigeons

This study was to investigate the correlations of myogenic differentiation 1 (MYOD1) gene polymorphisms with carcass traits and its expression with breast muscle development in pigeons. Four SNPs were found in the pigeon MYOD1 gene. Correlation analysis showed that individuals with AA genotype at both SNPs g.2967A > G (p < .01) and g.3044G > A (p < .05) have significantly higher live weight (LW), carcass weight (CW), semi-eviscerated weight (SEW), eviscerated weight (EW) and breast muscle weight (BMW). Moreover, the two SNPs also had the same significant effects on MYOD1 mRNA expression levels in breast muscle of pigeons, ie, the AA genotype showed higher MYOD1 mRNA expression levels. The diameter and cross-section area of muscle fibers continuously increased from 0w to 4w (p < .05), accompanied with the increasing expression of MYOD1 gene, while the density decreased (p < .05) dramatically from 0w to 1w and continuously fell over in the next few weeks (p > .05). What's more, the expression level of MYOD1 gene was positively correlated with a diameter (r = 0.937, p < .05) and cross-sectional area (r = 0.956, p < .01) of myofiber, and negatively correlated with density (r = -0.769, p < .01). The results showed that individuals with AA genotype at both SNPs g.2967A > G and g.3044G > A have showed higher carcass traits (LW, CW, SEW, EW, and BMW) and higher MYOD1 mRNA expression level in breast muscle than AB and BB genotypes. Moreover, the expression level of MYOD1 gene was closely correlated with muscle characteristic traits, indicating variants of MYOD1 gene was closely related to muscle development and could be a potential candidate gene in marker-assisted selection of pigeons.

Integration of Transcriptomics and Non-Targeted Metabolomics Reveals the Underlying Mechanism of Skeletal Muscle Development in Duck during Embryonic Stage

Skeletal muscle is an important economic trait in duck breeding; however, little is known about the molecular mechanisms of its embryonic development. Here, the transcriptomes and metabolomes of breast muscle of Pekin duck from 15 (E15_BM), 21 (E21_BM), and 27 (E27_BM) days of incubation were compared and analyzed. The metabolome results showed that the differentially accumulated metabolites (DAMs), including the up-regulated metabolites, l-glutamic acid, n-acetyl-1-aspartylglutamic acid, l-2-aminoadipic acid, 3-hydroxybutyric acid, bilirubin, and the significantly down-regulated metabolites, palmitic acid, 4-guanidinobutanoate, myristic acid, 3-dehydroxycarnitine, and s-adenosylmethioninamine, were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites, biosynthesis of cofactors, protein digestion and absorption, and histidine metabolism, suggesting that these pathways may play important roles in the muscle development of duck during the embryonic stage. Moreover, a total of 2142 (1552 up-regulated and 590 down-regulated), 4873 (3810 up-regulated and 1063 down-regulated), and 2401 (1606 up-regulated and 795 down-regulated) DEGs were identified from E15_BM vs. E21_BM, E15_BM vs. E27_BM and E21_BM vs. E27_BM in the transcriptome, respectively. The significantly enriched GO terms from biological processes were positive regulation of cell proliferation, regulation of cell cycle, actin filament organization, and regulation of actin cytoskeleton organization, which were associated with muscle or cell growth and development. Seven significant pathways, highly enriched by FYN, PTK2, PXN, CRK, CRKL, PAK, RHOA, ROCK, INSR, PDPK1, and ARHGEF, were focal adhesion, regulation of actin cytoskeleton, wnt signaling pathway, insulin signaling pathway, extracellular matrix (ECM)-receptor interaction, cell cycle, and adherens junction, which participated in regulating the development of skeletal muscle in Pekin duck during the embryonic stage. KEGG pathway analysis of the integrated transcriptome and metabolome indicated that the pathways, including arginine and proline metabolism, protein digestion and absorption, and histidine metabolism, were involved in regulating skeletal muscle development in embryonic Pekin duck. These findings suggested that the candidate genes and metabolites involved in crucial biological pathways may regulate muscle development in the Pekin duck at the embryonic stage, and increased our understanding of the molecular mechanisms underlying the avian muscle development.

Polymorphisms of CRELD1 and DNAJC30 and their relationship with chicken carcass traits

Carcass traits play important roles in the broiler industry and single nucleotide polymorphism (SNP) can be efficient molecular markers for marker-assisted breeding of chicken carcass traits. Based on our previous RNA-seq data (accession number GSE58755), cysteine rich with epidermal growth factor like domains 1 (CRELD1) and DnaJ heat shock protein family member C30 (DNAJC30) are differentially expressed in breast muscle between white recessive rock chicken (WRR) and Xinghua chicken (XH). In this study, we further characterize the potential function and SNP mutation of CRELD1 and DNAJC30 in chicken for the first time. According to protein interaction network and enrichment analysis, CRELD1 and DNAJC30 may play some roles in chicken muscle development and fat deposition. In WRR and XH, the results of the relative tissue expression pattern demonstrated that CRELD1 and DNAJC30 are not only differentially expressed in breast muscle but also leg muscle and abdominal fat. Therefore, we identified 5 SNP sites of CRELD1 and 7 SNP sites of DNAJC30 and genotyped them in an F₂ chicken population. There are 4 sites of CRELD1 and 3 sites of DNAJC30 are associated with chicken carcass traits like breast muscle weight, body weight, dressed weight, leg weight percentage, eviscerated weight with giblet percentage, intermuscular adipose width, shank length, and girth. These results suggest that the SNP sites of CRELD1 and DNAJC30 can be potential molecular markers to improve the chicken carcass traits and lay the foundation for marker-assisted selection.

Assessing genomic diversity and signatures of selection in Jiaxian Red cattle using whole-genome sequencing data

Background

Native cattle breeds are an important source of genetic variation because they might carry alleles that enable them to adapt to local environment and tough feeding conditions. Jiaxian Red, a Chinese native cattle breed, is reported to have originated from crossbreeding between taurine and indicine cattle; their history as a draft and meat animal dates back at least 30 years. Using whole-genome sequencing (WGS) data of 30 animals from the core breeding farm, we investigated the genetic diversity, population structure and genomic regions under selection of Jiaxian Red cattle. Furthermore, we used 131 published genomes of world-wide cattle to characterize the genomic variation of Jiaxian Red cattle.

Results

The population structure analysis revealed that Jiaxian Red cattle harboured the ancestry with East Asian taurine (0.493), Chinese indicine (0.379), European taurine (0.095) and Indian indicine (0.033). Three methods (nucleotide diversity, linkage disequilibrium decay and runs of homozygosity) implied the relatively high genomic diversity in Jiaxian Red cattle. We used θπ, CLR, F_ST and XP-EHH methods to look for the candidate signatures of positive selection in Jiaxian Red cattle. A total number of 171 (θπ and CLR) and 17 (F_ST and XP-EHH) shared genes were identified using different detection strategies. Functional annotation analysis revealed that these genes are potentially responsible for growth and feed efficiency (CCSER1), meat quality traits (ROCK2, PPP1R12A, CYB5R4, EYA3, PHACTR1), fertility (RFX4, SRD5A2) and immune system response (SLAMF1, CD84 and SLAMF6).

Conclusion

We provide a comprehensive overview of sequence variations in Jiaxian Red cattle genomes. Selection signatures were detected in genomic regions that are possibly related to economically important traits in Jiaxian Red cattle. We observed a high level of genomic diversity and low inbreeding in Jiaxian Red cattle. These results provide a basis for further resource protection and breeding improvement of this breed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07340-0.

Genome-wide DNA methylation profiles reveal novel candidate genes associated with meat quality at different age stages in hens

Poultry meat quality is associated with breed, age, tissue and other factors. Many previous studies have focused on distinct breeds; however, little is known regarding the epigenetic regulatory mechanisms in different age stages, such as DNA methylation. Here, we compared the global DNA methylation profiles between juvenile (20 weeks old) and later laying-period (55 weeks old) hens and identified candidate genes related to the development and meat quality of breast muscle using whole-genome bisulfite sequencing. The results showed that the later laying-period hens, which had a higher intramuscular fat (IMF) deposition capacity and water holding capacity (WHC) and less tenderness, exhibited higher global DNA methylation levels than the juvenile hens. A total of 2,714 differentially methylated regions were identified in the present study, which corresponded to 378 differentially methylated genes, mainly affecting muscle development, lipid metabolism, and the ageing process. Hypermethylation of the promoters of the genes ABCA1, COL6A1 and GSTT1L and the resulting transcriptional down-regulation in the later laying-period hens may be the reason for the significant difference in the meat quality between the juvenile and later laying-period hens. These findings contribute to a better understanding of epigenetic regulation in the skeletal muscle development and meat quality of chicken.

Whole genome scan reveals the genetic signature of African Ankole cattle breed and potential for higher quality beef

BACKGROUND

Africa is home to numerous cattle breeds whose diversity has been shaped by subtle combinations of human and natural selection. African Sanga cattle are an intermediate type of cattle resulting from interbreeding between Bos taurus and Bos indicus subspecies. Recently, research has asserted the potential of Sanga breeds for commercial beef production with better meat quality as compared to Bos indicus breeds. Here, we identified meat quality related gene regions that are positively selected in Ankole (Sanga) cattle breeds as compared to indicus (Boran, Ogaden, and Kenana) breeds using cross-population (XP-EHH and XP-CLR) statistical methods.

RESULTS

We identified 238 (XP-EHH) and 213 (XP-CLR) positively selected genes, of which 97 were detected from both statistics. Among the genes obtained, we primarily reported those involved in different biological process and pathways associated with meat quality traits. Genes (CAPZB, COL9A2, PDGFRA, MAP3K5, ZNF410, and PKM2) involved in muscle structure and metabolism affect meat tenderness. Genes (PLA2G2A, PARK2, ZNF410, MAP2K3, PLCD3, PLCD1, and ROCK1) related to intramuscular fat (IMF) are involved in adipose metabolism and adipogenesis. MB and SLC48A1 affect meat color. In addition, we identified genes (TIMP2, PKM2, PRKG1, MAP3K5, and ATP8A1) related to feeding efficiency. Among the enriched Gene Ontology Biological Process (GO BP) terms, actin cytoskeleton organization, actin filament-based process, and protein ubiquitination are associated with meat tenderness whereas cellular component organization, negative regulation of actin filament depolymerization and negative regulation of protein complex disassembly are involved in adipocyte regulation. The MAPK pathway is responsible for cell proliferation and plays an important role in hyperplastic growth, which has a positive effect on meat tenderness.

CONCLUSION

Results revealed several candidate genes positively selected in Ankole cattle in relation to meat quality characteristics. The genes identified are involved in muscle structure and metabolism, and adipose metabolism and adipogenesis. These genes help in the understanding of the biological mechanisms controlling beef quality characteristics in African Ankole cattle. These results provide a basis for further research on the genomic characteristics of Ankole and other Sanga cattle breeds for quality beef.

Gene set based association analyses for the WSSV resistance of Pacific white shrimp Litopenaeus vannamei

White Spot Syndrome Virus (WSSV) is regarded as a virus with the strongest pathogenicity to shrimp. For the threshold trait such as disease resistance, marker assisted selection (MAS) was considered to be a more effective approach. In the present study, association analyses of single nucleotide polymorphisms (SNPs) located in a set of immune related genes were conducted to identify markers associated with WSSV resistance. SNPs were detected by bioinformatics analysis on RNA sequencing data generated by Illimina sequencing platform and Roche 454 sequencing technology. A total of 681 SNPs located in the exons of immune related genes were selected as candidate SNPs. Among these SNPs, 77 loci were genotyped in WSSV susceptible group and resistant group. Association analysis was performed based on logistic regression method under an additive and dominance model in GenABEL package. As a result, five SNPs showed associations with WSSV resistance at a significant level of 0.05. Besides, SNP-SNP interaction analysis was conducted. The combination of SNP loci in TRAF6, Cu/Zn SOD and nLvALF2 exhibited a significant effect on the WSSV resistance of shrimp. Gene expression analysis revealed that these SNPs might influence the expression of these immune-related genes. This study provides a useful method for performing MAS in shrimp.

Identification of quantitative trait loci for body temperature, body weight, breast yield, and digestibility in an advanced intercross line of chickens under heat stress

Background

Losses in poultry production due to heat stress have considerable negative economic consequences. Previous studies in poultry have elucidated a genetic influence on response to heat. Using a unique chicken genetic resource, we identified genomic regions associated with body temperature (BT), body weight (BW), breast yield, and digestibility measured during heat stress. Identifying genes associated with a favorable response during high ambient temperature can facilitate genetic selection of heat-resilient chickens.

Methods

Generations F18 and F19 of a broiler (heat-susceptible) × Fayoumi (heat-resistant) advanced intercross line (AIL) were used to fine-map quantitative trait loci (QTL). Six hundred and thirty-one birds were exposed to daily heat cycles from 22 to 28 days of age, and phenotypes were measured before heat treatment, on the 1st day and after 1 week of heat treatment. BT was measured at these three phases and BW at pre-heat treatment and after 1 week of heat treatment. Breast muscle yield was calculated as the percentage of BW at day 28. Ileal feed digestibility was assayed from digesta collected from the ileum at day 28. Four hundred and sixty-eight AIL were genotyped using the 600 K Affymetrix chicken SNP (single nucleotide polymorphism) array. Trait heritabilities were estimated using an animal model. A genome-wide association study (GWAS) for these traits and changes in BT and BW was conducted using Bayesian analyses. Candidate genes were identified within 200-kb regions around SNPs with significant association signals.

Results

Heritabilities were low to moderate (0.03 to 0.35). We identified QTL for BT on Gallus gallus chromosome (GGA)14, 15, 26, and 27; BW on GGA1 to 8, 10, 14, and 21; dry matter digestibility on GGA19, 20 and 21; and QTL of very large effect for breast muscle yield on GGA1, 15, and 22 with a single 1-Mb window on GGA1 explaining more than 15 % of the genetic variation.

Conclusions

This is the first study to estimate heritabilities and perform GWAS using this AIL for traits measured during heat stress. Significant QTL as well as low to moderate heritabilities were found for each trait, and these QTL may facilitate selection for improved animal performance in hot climatic conditions.

SNP discovery from transcriptome of the swimbladder of Takifugu rubripes

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies in many species. High-throughput sequencing of RNA was developed primarily to analyze global gene expression, while it is an efficient way to discover SNPs from the expressed genes. In this study, we conducted transcriptome sequencing of the swimbladder of Takifugu rubripes using Illumina HiSeq2000 platform to identify gene-associated SNPs in the swimbladder. A total of 30,312,181 unique-mapped-reads were obtained from 44,736,850 raw reads. A total of 62,270 putative SNPs were discovered, which were located in 11,306 expressed genes and 2,246 scaffolds. The average minor allele frequency (MAF) of the SNPs was 0.26. GO and KEGG pathway analysis were conducted to analyze the genes containing SNPs. Validation of selected SNPs revealed that 54% of SNPs (26/48) were true SNPs. The results suggest that RNA-Seq is an efficient and cost-effective approach to discover gene-associated SNPs. In this study, a large number of SNPs were identified and these data will be useful resources for population genetic study, evolution analysis, resource assessment, genetic linkage analysis and genome-wide association studies.

The sodium channel gene family is specifically expressed in hen uterus and associated with eggshell quality traits

BACKGROUND

Eggshell quality is important for the poultry industry. During eggshell formation a mass of inorganic minerals is deposited. The Sodium Channel (SCNN1) gene family plays an essential role in cation transportation. The objective of this study was to investigate the pattern of expression of members of the SCNN1 gene family, their variation and their effects on eggshell quality.

RESULT

The highest expression of SCNN1a, SCNN1b, and SCNN1g genes were in the active uterus during eggshell mineralization, while SCNN1d showed its highest expression level in the quiescent uterus (no egg present). Nineteen candidate SNPs from the four genes were genotyped in a population of 338 White Leghorn layers. Association analysis between SNPs (haplotypes/diplotypes) and eggshell traits was performed. Among seven significant SNPs, five SNPs were associated with eggshell strength, eggshell thickness, eggshell percentage or/and egg weight, while the other two SNPs within SCNN1d were only associated with eggshell percentage. These SNPs had a 0.25-6.99% contribution to phenotypic variance, depending on the trait. In haplotype analysis, SCNN1b and SCNN1d were associated with egg weight. The SCNN1b and SCNN1g were significantly associated with eggshell weight while only SCNN1g explained 2.04% of phenotypic variance. All the alleles of the members of SCNN1 gene family were associated with eggshell percentage and eggshell thickness, and others members had an association with eggshell strength except for SCNN1a. The contribution of different haplotypes of the SCNN1 gene family to eggshell phenotypic variance ranged from 0.09% to 5.74%.

CONCLUSIONS

Our study indicated that the SCNN1 gene family showed tissue expression specificity and was significantly associated with eggshell traits in chicken. This study provides evidence that genetic variation in members of the sodium channel can influence eggshell quality.

Comparison of the genome-wide DNA methylation profiles between fast-growing and slow-growing broilers

INTRODUCTION

Growth traits are important in poultry production, however, little is known for its regulatory mechanism at epigenetic level. Therefore, in this study, we aim to compare DNA methylation profiles between fast- and slow-growing broilers in order to identify candidate genes for chicken growth. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) was used to investigate the genome-wide DNA methylation pattern in high and low tails of Recessive White Rock (WRR(h); WRR(l)) and that of Xinhua Chickens (XH(h); XH(l)) at 7 weeks of age. The results showed that the average methylation density was the lowest in CGIs followed by promoters. Within the gene body, the methylation density of introns was higher than that of UTRs and exons. Moreover, different methylation levels were observed in different repeat types with the highest in LINE/CR1. Methylated CGIs were prominently distributed in the intergenic regions and were enriched in the size ranging 200-300 bp. In total 13,294 methylated genes were found in four samples, including 4,085 differentially methylated genes of WRR(h) Vs. WRR(l), 5,599 of XH(h) Vs. XH(l), 4,204 of WRR(h) Vs. XH(h), as well as 7,301 of WRR(l) Vs. XH(l). Moreover, 132 differentially methylated genes related to growth and metabolism were observed in both inner contrasts (WRR(h) Vs. WRR(l) and XH(h) Vs. XH(l)), whereas 129 differentially methylated genes related to growth and metabolism were found in both across-breed contrasts (WRR(h) Vs. XH(h) and WRR(l) Vs. XH(l)). Further analysis showed that overall 75 genes exhibited altered DNA methylation in all four contrasts, which included some well-known growth factors of IGF1R, FGF12, FGF14, FGF18, FGFR2, and FGFR3. In addition, we validate the MeDIP-seq results by bisulfite sequencing in some regions.

CONCLUSIONS

This study revealed the global DNA methylation pattern of chicken muscle, and identified candidate genes that potentially regulate muscle development at 7 weeks of age at methylation level.