Transcriptomic analysis to elucidate the molecular mechanisms that underlie feed efficiency in meat-type chickens

Characterisation of runs of homozygosity and inbreeding coefficients in the red-brown Korean native chickens

OBJECTIVE

The analysis of runs of homozygosity (ROH) has been applied to assess the level of inbreeding and identify selection signatures in various livestock species. The objectives of this study were to characterize the ROH pattern, estimate the rate of inbreeding, and identify signatures of selection in the red-brown Korean native chickens.

METHODS

The Illumina 60K single nucleotide polymorphism chip data of 651 chickens was used in the analysis. Runs of homozygosity were analysed using the PLINK v1.9 software. Inbreeding coefficients were estimated using the GCTA software and their correlations were examined. Genomic regions with high levels of ROH were explored to identify selection signatures.

RESULTS

A total of 32,176 ROH segments were detected in this study. The majority of the ROH segments were shorter than 4 Mb. The average ROH inbreeding coefficients (FROH) varied with the length of ROH segments. The means of inbreeding coefficients calculated from different methods were also variable. The correlations between different inbreeding coefficients were positive and highly variable (r = 0.18-1). Five ROH islands harbouring important quantitative trait loci were identified.

CONCLUSION

This study assessed the level of inbreeding and patterns of homozygosity in Red-brown native Korean chickens. The results of this study suggest that the level of recent inbreeding is low which indicates substantial progress in the conservation of red-brown Korean native chickens. Additionally, Candidate genomic regions associated with important production traits were detected in homozygous regions.

Genetic assessment and candidate genes identification for breed-specific characteristics of Qingyuan partridge chicken based on runs of homozygosity

BACKGROUND

Several core breeding and supporting lines of the Qingyuan partridge chicken, a representative local chicken breed in China, have been developed over 20 years. Consequently, its economic traits related to growth and reproduction have been significantly improved by breeding selection and commercial utilization, but some characteristic traits, such as partridge feathers, high meat quality and sufficient flavor, have always been retained. However, effective methods for genetic assessment and functional gene exploration of similar trait groups are lacking. The presence of identical haplotype fragments transmitted from parent to offspring results in runs of homozygosity (ROH), which offer an efficient solution. In this study, genomes of 134 Qingyuan partridge chickens representing two breeding populations and one preserved population were re-sequenced to evaluate the genetic diversity and explore functional genes by analyzing the diversity, distribution, and frequency of ROH.

RESULTS

The results showed a low level of genomic linkage and degree of inbreeding within both the bred and preserved populations, suggesting abundant genetic diversity and an adequate genetic potential of the Qingyuan partridge chicken. Throughout the long-term selection process, 21 genes, including GLI3, ANO5, BLVRA, EFNB2, SLC5A12, and SVIP, associated with breed-specific characteristics were accumulated within three ROH islands, whereas another 21 genes associated with growth traits including IRX1, IRX2, EGFR, TPK1, NOVA1, BDNF and so on were accumulated within five ROH islands.

CONCLUSIONS

These findings provide new insights into the genetic assessment and identification of genes with breed-specific and selective characteristics, offering a solid genetic basis for breeding and protection of Qingyuan partridge chickens.

OMICs approaches and technologies for understanding low-high feed efficiency traits in chicken: implication to breeding

In poultry production, there has been a trend of continuous increase in cost of feed ingredients which represents the major proportion of the production costs. Feed costs can be reduced by improving feed efficiency traits which increase the possibility of using various indigestible feed sources and decrease the environmental impact of the enhanced poultry production. Therefore, feed efficiency has been used as one of the most important economic traits of selection in the breeding program of chickens. Recently, many OMICs experimental studies have been designed to characterize biological differences between the high and low feed efficiency chicken phenotypes. Biological complexity cannot be fully captured by main individual OMICs such as genomics, transcriptomics, proteomics and metabolomics. Therefore, researchers have combined multiple assays from the same set of samples to create multi-OMICs datasets. OMICs findings are crucial in improving existing approaches to poultry breeding. The current review aimed to highlight the components of feed efficiency and general OMICs approaches and technologies. Besides, individual and multi-OMICs based understanding of chicken feed efficiency traits and the application of the acquired knowledge in the chicken breeding program were addressed.

Conservation priority and run of homozygosity pattern assessment of global chicken genetic resources

The conservation of genetic resources is becoming increasingly important for the sustainable development of the poultry industry. In the present study, we systematically analyzed the population structure, conservation priority, runs of homozygosity (ROH) of chicken breeds globally, and proposed rational conservation strategies. We used a 600K Affymetrix Axiom HD genotyping SNP array dataset of 2,429 chickens from 134 populations. The chickens were divided into 5 groups based on their country of origin and sampling location: Asian chickens (AS-LOC), African chickens (AF), European local chickens (EU-LOC), Asian breeds sampled in Germany (AS-DE), and European breeds sampled in Germany (EU-DE). The results indicated that the population structure was consistent with the actual geographical distribution of the populations. AS-LOC had the highest positive contribution to the total gene (HT, 1.00%,) and allelic diversity (AT, 0.0014%), the lowest inbreeding degree and the fastest linkage disequilibrium (LD) decay rate; the lowest contribution are derived by European ex situ chicken breeds (EU-DE:HT = -0.072%, AT = -0.0014%), which showed the highest inbreeding and slowest LD decay. Breeds farmed in ex situ (AS-DE, EU-DE) conditions exhibited reduced genetic diversity and increased inbreeding due to small population size. Given limited funds, it is a better choice for government to conserve the breeds with the highest contribution to genetic diversity in each group. Therefore, we evaluated the contribution of each breed to genetic and allelic diversity in 5 groups. Among each group, KUR(AF), BANG(AS-LOC), ALxx(EU-LOC), BHwsch(AS-DE), and ARw(EU-DE) had the highest contribution to gene diversity in the order of the above grouping. Similarly, according to the allelic diversity standard (in the same order), ZIMxx, PIxx, ALxx, SHsch, and ARsch had the highest contribution. After analyzing ROH, we found a total of 144,708 fragments and 27 islands. The gene and genome regions identified by the ROH islands and QTLs indicate that chicken breeds have potential for adaptation to different production systems. Based on these findings, it is recommended to prioritize the conservation of breeds with the highest genetic diversity in each group, while paying more attention to the conservation of Asian and African breeds. Furthermore, providing a valuable reference for the conservation and utilization of chicken.

Thigh muscle metabolic response is linked to feed efficiency and meat characteristics in slow-growing chicken

The Korat chicken (KR) is a slow-growing Thai chicken breed with relatively poor feed efficiency (FE) but very tasty meat with high protein and low fat contents, and a unique texture. To enhance the competitiveness of KR, its FE should be improved. However, selecting for FE has an unknown effect on meat characteristics. Thus, understanding the genetic basis underlying FE traits and meat characteristics is needed. In this study, 75 male KR birds were raised up to 10 wk of age. For each bird, the feed conversion ratio (FCR), residual feed intake (RFI), and physicochemical properties, flavor precursors, and biological compounds in the thigh meat were evaluated. At 10 wk of age, thigh muscle samples from 6 birds (3 with high FCR and 3 with low FCR values) were selected, and their proteomes were investigated using a label-free proteomic method. Weighted gene coexpression network analysis (WGCNA) was used to screen the key protein modules and pathways. The WGCNA results revealed that FE and meat characteristics significantly correlated with the same protein module. However, the correlation was unfavorable; improving FE may result in a decrease in meat quality through the alteration in biological processes including glycolysis/gluconeogenesis, metabolic pathway, carbon metabolism, biosynthesis of amino acids, pyruvate metabolism, and protein processing in the endoplasmic reticulum. The hub proteins of the significant module (TNNT1, TNNT3, TNNI2, TNNC2, MYLPF, MYH10, GADPH, PGK1, LDHA, and GPI) were also identified to be associated with energy metabolism, and muscle growth and development. Given that the same proteins and pathways are present in FE and meat characteristics but in opposite directions, selection practices for KR should simultaneously consider both trait groups to maintain the high meat quality of KR while improving FE.

Analysis of Conservation Priorities and Runs of Homozygosity Patterns for Chinese Indigenous Chicken Breeds

To achieve sustainable development of the poultry industry, the effective conservation of genetic resources has become increasingly important. In the present study, we systematically elucidated the population structure, conservation priority, and runs of homozygosity (ROH) patterns of Chinese native chicken breeds. We used a high-density genotyping dataset of 157 native chickens from eight breeds. The population structure showed different degrees of population stratification among the breeds. Chahua chicken was the most differentiated breed from the other breeds (Nei = 0.0813), and the Wannan three-yellow chicken (WanTy) showed the lowest degree of differentiation (Nei = 0.0438). On the basis of contribution priority, Xiaoshan chicken had the highest contribution to the total gene diversity (1.41%) and the maximum gene diversity of the synthetic population (31.1%). WanTy chicken showed the highest contribution to the total allelic diversity (1.31%) and the maximum allelic diversity of the syntenic population (17.0%). A total of 5242 ROH fragments and 5 ROH island regions were detected. The longest ROH fragment was 41.51 Mb. A comparison of the overlapping genomic regions between the ROH islands and QTLs in the quantitative trait loci (QTL) database showed that the annotated candidate genes were involved in crucial economic traits such as immunity, carcass weight, drumstick and leg muscle development, egg quality and egg production, abdominal fat precipitation, body weight, and feed intake. In conclusion, our findings revealed that Chahua, Xiaoshan, and WanTy should be the priority conservation breeds, which will help optimize the conservation and breeding programs for Chinese indigenous chicken breeds.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Basal diet composition contributes to differential performance, intestinal health, and immunological responses to a microalgae-based feed ingredient in broiler chickens

Novel feed ingredients may improve poultry health, but functionality of these ingredients may vary across basal diet formulations. This study evaluated a proprietary algae ingredient's effects on broiler performance, intestinal health, systemic immunity, and metabolic/immune kinotypes between corn- or wheat-based diets. Ross 308 broilers were housed in 80 floor pens (14 birds/pen) and assigned to 1 of 4 corn or wheat-based diets ± 0.175% algae ingredient for 42 d. At the end of each 14 d starter, grower, and finisher period, 10 birds/treatment were euthanized for tissue collection to assess intestinal histomorphology, systemic immune cell populations by flow cytometry and kinotypes by peptide arrays. On d 28 and 29, forty-three birds/treatment underwent a 12 h feed restriction challenge followed by a fluorescein isothiocyanate-dextran intestinal permeability assay. For the entire 42 d study, wheat-based diets improved feed conversion rate (FCR) by 5 points compared to corn-based diets (P < 0.0001). Performance benefits related to algae inclusion were diet dependent, with algae inclusion improving 42 d FCR by 6 points only in corn-based diets (P = 0.006). Birds fed wheat-based diets had reduced splenic monocyte/macrophage, CD1.1⁺, and T cell populations in the first 14 d (P < 0.0001) and reduced serum fluorescence on d 28/29 (P = 0.0002). Algae inclusion in the corn-based diet increased villus height in the duodenum on d 28 and jejunum on d 42, while reducing splenic CD3⁺CD8α⁺ cytotoxic T cells 13.4 to 27.5% compared to the corn-based control at the same timepoints (P < 0.0001). Kinome results showed a significant innate immune toll-like receptor (TLR) response via MyD88 at d 14 in the small intestine of birds fed corn-based diets with algae that shifted to a more growth factor and adaptive immune-oriented response by d 42. Concurrent with immune changes, signaling changes indicative of lipid metabolism in the small intestine, ceca, and liver were seen in birds fed the corn-based diet with algae. The observed differential responses to basal diet composition and algae inclusion emphasize the need to comparatively evaluate feed ingredients in various diet formulations.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

Microbiota and Transcriptomic Effects of an Essential Oil Blend and Its Delivery Route Compared to an Antibiotic Growth Promoter in Broiler Chickens

This study evaluated the effect of the delivery of a commercial essential oil blend containing the phytonutrients star anise, cinnamon, rosemary, and thyme oil (via different routes) on broiler chickens’ ileal and ceca microbiota and liver transcriptome compared to an antibiotic growth promoter. Eggs were incubated and allocated into three groups: non-injected, in ovo saline, and in ovo essential oil. On day 18 of incubation, 0.2 mL of essential oil in saline (dilution ratio of 2:1) or saline alone was injected into the amnion. At hatch, chicks were assigned to post-hatch treatment combinations: (A) a negative control (corn-wheat-soybean diet), (B) in-feed antibiotics, (C) in-water essential oil (250 mL/1000 L of drinking water), (D) in ovo saline, (E) in ovo essential oil, and (F) in ovo essential oil plus in-water essential oil in eight replicate cages (six birds/cage) and raised for 28 days. On days 21 and 28, one and two birds per cage were slaughtered, respectively, to collect gut content and liver tissues for further analysis. Alpha and beta diversity differed significantly between ileal and ceca samples but not between treatment groups. In-feed antibiotic treatment significantly increased the proportion of specific bacteria in the family Lachnospiraceae while reducing the proportion of bacteria in the genus Christensenellaceae in the ceca, compared to other treatments. Sex-controlled differential expression of genes related to cell signaling and tight junctions were recorded. This study provides data that could guide the use of these feed additives and a foundation for further research.

Comparative proteomics revealed duodenal metabolic function associated with feed efficiency in slow-growing chicken

The Korat chicken (KR), developed in Thailand, is a slow-growing breed developed as an alternative breed for Thai chicken producers. The growing interest in slow-growing chicken meat, due to its unique taste, distinct texture, health benefits, and higher broiler welfare have led to higher market demand for KR. However, its low feed efficiency (FE) has a significant negative impact on farm profitability. Understanding the molecular mechanism regulating FE allows for designing a suitable selection program and contributing to breeding more efficient chicken for poultry production. Thus, the objective of our study was to investigate the proteome differences and possible pathways associated with FE in male KR using a label-free quantitative proteomic approach. Seventy-five KR males were individually evaluated for FE, and duodenum samples from 6 animals (3 high-FE and 3 low-FE chickens) were collected at 10 wk of age for differential abundant proteins (DAPs), protein networks, functional enrichment, and pathway analyses. In this study, we found 40 DAPs significantly associated with FE pathways, including glycolysis/gluconeogenesis, peroxisome, oxidative phosphorylation, tight junction, and cysteine and methionine metabolism. Thus, variations in observed DAPs or genes related to DAPs could be interesting biomarker candidates for selection for higher feed utilization efficiency in chicken.

Genome-wide association and transcriptome studies identify candidate genes and pathways for feed conversion ratio in pigs

Background

The feed conversion ratio (FCR) is an important productive trait that greatly affects profits in the pig industry. Elucidating the genetic mechanisms underpinning FCR may promote more efficient improvement of FCR through artificial selection. In this study, we integrated a genome-wide association study (GWAS) with transcriptome analyses of different tissues in Yorkshire pigs (YY) with the aim of identifying key genes and signalling pathways associated with FCR.

Results

A total of 61 significant single nucleotide polymorphisms (SNPs) were detected by GWAS in YY. All of these SNPs were located on porcine chromosome (SSC) 5, and the covered region was considered a quantitative trait locus (QTL) region for FCR. Some genes distributed around these significant SNPs were considered as candidates for regulating FCR, including TPH2, FAR2, IRAK3, YARS2, GRIP1, FRS2, CNOT2 and TRHDE. According to transcriptome analyses in the hypothalamus, TPH2 exhibits the potential to regulate intestinal motility through serotonergic synapse and oxytocin signalling pathways. In addition, GRIP1 may be involved in glutamatergic and GABAergic signalling pathways, which regulate FCR by affecting appetite in pigs. Moreover, GRIP1, FRS2, CNOT2, and TRHDE may regulate metabolism in various tissues through a thyroid hormone signalling pathway.

Conclusions

Based on the results from GWAS and transcriptome analyses, the TPH2, GRIP1, FRS2, TRHDE, and CNOT2 genes were considered candidate genes for regulating FCR in Yorkshire pigs. These findings improve the understanding of the genetic mechanisms of FCR and may help optimize the design of breeding schemes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07570-w.

The significant influence of residual feed intake on flavor precursors and biomolecules in slow-growing Korat chicken meat

Objective

This study investigated the association between feed efficiency, physicochemical properties, flavor precursors and biomolecules in the thigh meat of Korat (KR) chickens.

Methods

The feed intake and body weight of individual male KR chickens were recorded from 1 to 10 weeks old to calculate the individual residual feed intake (RFI) of 75 birds. At 10 weeks of age, chickens with the 10 highest (HRFI) and lowest RFI (LRFI) were slaughtered to provide thigh meat samples. The physicochemical properties (ultimate pH, water holding capacity [WHC], drip loss) and flavor precursors (guanosine monophosphate, inosine monophosphate (IMP), adenosine monophosphate and inosine) were analyzed conventionally, and Fourier transform infrared spectroscopy was used to identify the composition of biomolecules (lipids, ester lipids, amide I, amide II, amide III, and carbohydrates) and the secondary structure of the proteins. A group t-test was used to determine significant differences between mean values and principal component analysis to classify thigh meat samples into LRFI and HRFI KR chickens.

Results

The physicochemical properties of thigh meat samples from LRFI and HRFI KR chickens were not significantly different but the IMP content, ratios of lipid, lipid ester, protein (amide I, amide II) were significantly different (p<0.05). The correlation loading results showed that the LRFI group was correlated with high ratios of lipids, lipid esters, collagen content (amide III) and beta sheet protein (rg loading >0.5) while the HRFI group was positively correlated with protein (amide I, amide II), alpha helix protein, IMP content, carbohydrate, ultimate pH and WHC (rg loading >0.5).

Conclusion

The thigh meat from chickens with different RFI differed in physiochemical properties affecting meat texture, and in the contents of flavor precursors and biomolecules affecting the nutritional value of meat. This information can help animal breeders to make genetic improvements by taking more account of traits related to RFI.

Genome-wide analysis reveals the patterns of genetic diversity and population structure of 8 Italian local chicken breeds

The aim of this study was to conduct a genome-wide comparative analysis of 8 local Italian chicken breeds (Ermellinata di Rovigo, Millefiori di Lonigo [PML], Polverara Bianca, Polverara Nera, Padovana, Pepoi [PPP], Robusta Lionata, and Robusta Maculata), all under a conservation plan, to understand their genetic diversity and population structure. A total of 152 animals were analyzed using the Affymetrix Axiom 600 K Chicken Genotyping Array. The levels of genetic diversity were highest and lowest in PML and PPP, respectively. The results of genomic inbreeding based on runs of homozygosity (ROH; FROH) showed marked differences among breeds and ranged from 0.161 (PML) to 0.478 (PPP). Furthermore, in all breeds, short ROH (<4 Mb in length) were more frequent than long segments. Patterns of genetic differentiation, model-based clustering, and neighbor networks showed that most breeds formed nonoverlapping clusters and were clearly separate populations. The 2 Polverara breeds shared a similar genetic background and showed the lowest genetic differentiation in comparison with purebred lines; the local populations showed separated groups. PPP and PML were closer to the group of the purebred broiler lines (BRSA, BRSB, BRDA, and BRDB). Six genomic regions are presented as hotspots of autozygosity among the Italian chicken breeds, with candidate genes involved in multiple morphological phenotypes as breast muscle, muscle dry matter content, and body weight. This study is the first exhaustive genome-wide analysis of the diversity of these Italian local chickens from Veneto region. We conclude that breeds have conserved authentic genetic patterns. The results are of significant importance because they will help design and implement conservation strategies. In fact, the conservation of these breeds may also have positive impacts on the local economy, niche traditional markets, and offering a source of high-quality products to consumers. In this context, genomic information may play a crucial role in the management of local breeds.

Genome-Wide SNP Analysis Reveals the Population Structure and the Conservation Status of 23 Italian Chicken Breeds

The genomic variability of local Italian chicken breeds, which were monitored under a conservation plan, was studied using single nucleotide polymorphisms (SNPs) to understand their genetic diversity and population structure. A total of 582 samples from 23 local breeds and four commercial stocks were genotyped using the Affymetrix 600 K Chicken SNP Array. In general, the levels of genetic diversity, investigated through different approaches, were lowest in the local chicken breeds compared to those in the commercial stocks. The level of genomic inbreeding, based on runs of homozygosity (FROH), was markedly different among the breeds and ranged from 0.121 (Valdarnese) to 0.607 (Siciliana). In all breeds, short runs of homozygosity (ROH) (<4 Mb in length) were more frequent than long segments. The patterns of genetic differentiation, model-based clustering, and neighbor networks showed that most breeds formed non-overlapping clusters and were clearly separate populations, which indicated the presence of gene flow, especially among breeds that originated from the same geographical area. Four genomic regions were identified as hotspots of autozygosity (islands) among the breeds, where the candidate genes are involved in morphological traits, such as body weight and feed conversion ratio. We conclude that the investigated breeds have conserved authentic genetic patterns, and these results can improve conservation strategies; moreover, the conservation of local breeds may play an important role in the local economy as a source of high-quality products for consumers.

Genome-wide Analyses Identifies Known and New Markers Responsible of Chicken Plumage Color

Simple Summary

In order to assess sources of variation related to Polverara breed plumage color (black vs. white), we carried out genome-wide analyses to identify the genomic regions involved in this trait. The present work has revealed new candidate genes involved in the phenotypic variability in local chicken populations. These results also contribute insights into the genetic basis for plumage color in poultry, and confirm the great complexity of the mechanisms that control this trait.

Abstract

Through the development of the high-throughput genotyping arrays, molecular markers and genes related to phenotypic traits have been identified in livestock species. In poultry, plumage color is an important qualitative trait that can be used as phenotypic marker for breed identification. In order to assess sources of genetic variation related to the Polverara chicken breed plumage colour (black vs. white), we carried out a genome-wide association study (GWAS) and a genome-wide fixation index (F_ST) scan to uncover the genomic regions involved. A total of 37 animals (17 white and 20 black) were genotyped with the Affymetrix 600 K Chicken single nucleotide polymorphism (SNP) Array. The combination of results from GWAS and F_ST revealed a total of 40 significant markers distributed on GGA 01, 03, 08, 12 and 21, and located within or near known genes. In addition to the well-known TYR, other candidate genes have been identified in this study, such as GRM5, RAB38 and NOTCH2. All these genes could explain the difference between the two Polverara breeds. Therefore, this study provides the basis for further investigation of the genetic mechanisms involved in plumage color in chicken.

Assessing selection signatures within and between selected lines of dual-purpose black and white and German Holstein cattle

The aim of this study was to detect selection signatures considering cows from the German Holstein (GH) and the local dual-purpose black and white (DSN) population, as well as from generated sub-populations. The 4654 GH and 261 DSN cows were genotyped with the BovineSNP50 Genotyping BeadChip. The geographical herd location was used as an environmental descriptor to create the East-DSN and West-DSN sub-populations. In addition, two further sub-populations of GH cows were generated, using the extreme values for solutions of residual effects of cows for the claw disorder dermatitis digitalis. These groups represented the most susceptible and most resistant cows. We used cross-population extended haplotype homozygosity methodology (XP-EHH) to identify the most recent selection signatures. Furthermore, we calculated Wright's fixation index (F_ST ). Chromosomal segments for the top 0.1 percentile of negative or positive XP-EHH scores were studied in detail. For gene annotations, we used the Ensembl database and we considered a window of 250 kbp downstream and upstream of each core SNP corresponding to peaks of XP-EHH. In addition, functional interactions among potential candidate genes were inferred via gene network analyses. The most outstanding XP-EHH score was on chromosome 12 (at 77.34 Mb) for DSN and on chromosome 20 (at 36.29-38.42 Mb) for GH. Selection signature locations harbored QTL for several economically important milk and meat quality traits, reflecting the different breeding goals for GH and DSN. The average F_ST value between GH and DSN was quite low (0.068), indicating shared founders. For group stratifications according to cow health, several identified potential candidate genes influence disease resistance, especially to dermatitis digitalis.

The association between microbial community and ileal gene expression on intestinal wall thickness alterations in chickens

The dynamic development of the animal intestine with a concurrent succession of microbiota and changes in microbial community and metabolite spectrum can exert far-reaching effects on host physiology. However, the precise mechanism of mutual response between microbiota and the gut is yet to be fully elucidated. Broilers with varying developmental degrees of intestinal wall thickness were selected, and they were divided into the thick group (H type) and the thin group (B type), using multiomics data integration analysis to reveal the fundamental regulatory mechanisms of gut–microbiota interplay. Our data showed, in broilers with similar body weight, the intestinal morphological parameters were improved in H type and the diversity of microbial communities is distinguishable from each other. The beneficial bacteria such as Bifidobacterium breve was increased whereas avian endogenous retrovirus EAV-HP was decreased in the H type compared with the B type. Furthermore, microbial metabolic potentials were more active, especially the biosynthesis of folate was improved in the H type. Similarly, the consolidation of absorption, immunity, metabolism, and development was noticed in the thick group. Correlation analysis indicated that the expression levels of material transport and immunomodulatory-related genes were positively correlated with the relative abundance of several probiotics such as B. breve, Lactobacillus saerimneri, and Butyricicoccus pullicaecorum. Our findings suggest that the chickens with well-developed ileal thickness own exclusive microbial composition and metabolic potential, which is closely related to small intestinal morphogenesis and homeostasis.

Effect of divergence in residual feed intake on expression of lipid metabolism-related genes in the liver of meat-type ducks1

Lipid metabolism is considered one of the important factors affecting residual feed intake (RFI). However, the relationship between RFI and expression of lipid metabolism-related genes is unknown in meat-type ducks. To address this issue, a total of 1,000 male meat-type ducks with similar body weight were randomly selected to measure body weight gain and feed intake from 21 to 42 d of age to estimate RFI. The 8 greatest- (high RFI [HRFI]) and lowest- (low RFI [LRFI]) ranking birds were then selected for the present study. Relative expressions of key genes, namely sirtuin 1 (Sirt1), forkhead box O1 (Foxo1), peroxisome proliferator-activated receptor gamma (PPARγ), sterol regulatory element-binding transcription factor 1c (SREBP-1c), fas cell surface death receptor (FAS), acetyl-CoA carboxylase alpha (ACC), carnitine palmitoyltransferase 1A (CPT1A), and acyl-CoA oxidase 1 (ACOX1), were then determined in the HRFI and LRFI ducks by quantitative PCR. The results showed that RFI, feed conversion ratio (FCR), and average daily feed intake (ADFI) were significantly lower (P < 0.05) in LRFI ducks than in HRFI ducks. In addition, expression of Sirt1, Foxo1, CPT1A, and ACOX1 were significantly higher in LRFI ducks than in HRFI ducks (P < 0.05), whereas PPARγ and FAS expression levels were significantly lower in LRFI ducks than in HRFI ducks (P < 0.01). Correlation analysis showed that Sirt1, CPT1A, and ACOX1 expressions were significantly negatively correlated with FCR (r = -0.81 to -0.93; P < 0.01), whereas PPARγ and FAS expressions were significantly positively correlated with FCR (r = 0.74 to 0.87; P < 0.01). PPARγ expression was significantly positively correlated with RFI (r = 0.83; P < 0.01), whereas CPT1A and ACOX1 expressions were significantly negatively correlated with RFI (r = -0.84 to -0.89; P < 0.01). Sirt1 mRNA expression was positively correlated with Foxo1, CPT1A, and ACOX1 mRNA expression (r = 0.78 to 0.92; P < 0.01). Association of Foxo1 with CPT1A and ACOX1 was positive (r = 0.88 to 0.96; P < 0.01). These results suggest that genes related to fatty acid oxidation are upregulated in the liver of ducks with high feed efficiency, while genes associated with lipid synthesis are downregulated. Furthermore, the inclusion of lipid metabolism-related genes in future breeding programs might be beneficial for selecting ducks with greater feed efficiency phenotype.

Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle

BACKGROUND

Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions.

RESULTS

The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15.

CONCLUSIONS

Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.

Distinct cell proliferation, myogenic differentiation, and gene expression in skeletal muscle myoblasts of layer and broiler chickens

Myoblasts play a central role during skeletal muscle formation and growth. Precise understanding of myoblast properties is thus indispensable for meat production. Herein, we report the cellular characteristics and gene expression profiles of primary-cultured myoblasts of layer and broiler chickens. Broiler myoblasts actively proliferated and promptly differentiated into myotubes compared to layer myoblasts, which corresponds well with the muscle phenotype of broilers. Transcriptomes of layer and broiler myoblasts during differentiation were quantified by RNA sequencing. Ontology analyses of the differentially expressed genes (DEGs) provided a series of extracellular proteins as putative markers for characterization of chicken myogenic cells. Another ontology analyses demonstrated that broiler myogenic cells are rich in cell cycle factors and muscle components. Independent of these semantic studies, principal component analysis (PCA) statistically defined two gene sets: one governing myogenic differentiation and the other segregating layers and broilers. Thirteen candidate genes were identified with a combined study of the DEGs and PCA that potentially contribute to proliferation or differentiation of chicken myoblasts. We experimentally proved that one of the candidates, enkephalin, an opioid peptide, suppresses myoblast growth. Our results present a new perspective that the opioids present in feeds may influence muscle development of domestic animals.

Host transcriptome and microbiome interaction modulates physiology of full-sibs broilers with divergent feed conversion ratio

Efficient livestock production relies on effective conversion of feed into body weight gain (BWG). High levels of feed conversion are especially important in production of broiler chickens, birds reared for meat, where economic margins are tight. Traits associated with improved broiler growth and feed efficiency have been subjected to intense genetic selection, but measures such as feed conversion ratio (FCR) remain variable, even between full siblings (sibs). Non-genetic factors such as the composition and function of microbial populations within different enteric compartments have been recognized to influence FCR, although the extent of interplay between hosts and their microbiomes is unclear. To examine host–microbiome interactions we investigated variation in the composition and functions of host intestinal-hepatic transcriptomes and the intestinal microbiota of full-sib broilers with divergent FCR. Progeny from 300 broiler families were assessed for divergent FCR set against shared genetic backgrounds and exposure to the same environmental factors. The seven most divergent full-sib pairs were chosen for analysis, exhibiting marked variation in transcription of genes as well as gut microbial diversity. Examination of enteric microbiota in low FCR sibs revealed variation in microbial community structure and function with no difference in feed intake compared to high FCR sibs. Gene transcription in low and high FCR sibs was significantly associated with the abundance of specific microbial taxa. Highly intertwined interactions between host transcriptomes and enteric microbiota are likely to modulate complex traits like FCR and may be amenable to selective modification with relevance to improving intestinal homeostasis and health.

Expression of genes related to lipid transport in meat-type ducks divergent for low or high residual feed intake

Objective

This study examined the effects of divergence in residual feed intake (RFI) on expression profiles of key genes related to lipid transport in the liver and duodenal epithelium and their associations with feed efficiency traits in meat-type ducks.

Methods

A total of 1,000 male ducks with similar body weight (1,042.1±87.2 g) were used in this study, and their individual RFI was calculated from 21 to 42 d of age. Finally, the 10 highest RFI (HRFI) and 10 lowest RFI (LRFI) ducks were chosen for examining the expression of key genes related to lipid transport in the liver and duodenal epithelium using quantitative polymerase chain reaction.

Results

In the liver, expression levels of albumin (ALB), CD36 molecule (CD36), fatty acid hydroxylase domain containing 2 (FAXDC2), and choline kinase alpha (CHKA) were significantly higher in LRFI ducks than in HRFI ducks (p<0.01); negative correlations (p<0.05) between expression levels of ALB, CD36, FAXDC2, and CHKA and RFI were detected in the liver. Additionally, ALB expression was strongly positively correlated (p<0.05) with CD36, FAXDC2, CHKA, and apolipoprotein H (APOH) expression in the liver. In duodenal epithelium, we found that mRNA levels of ALB, CD36, FAXDC2, and APOH were significantly higher in LRFI ducks than in HRFI ducks (p<0.01); RFI was strongly negatively correlated (p<0.05) with ALB, FAXDC2, and APOH expression, while ALB expression was strongly positively correlated with APOH expression (p<0.01) in duodenal epithelium. Furthermore, expression levels of both ALB and FAXDC2 genes were significantly associated with feed conversion ratio and RFI in both liver and duodenal epithelium (p<0.05).

Conclusion

Our findings therefore suggest that ALB and FAXDC2 genes might be used as potential gene markers designed to improve feed efficiency in future meat-type duck breeding programs.

Association analysis between feed efficiency and expression of key genes of the avTOR signaling pathway in meat-type ducks

Genes involved in the target of rapamycin (TOR) signaling pathway are implicated in nutrient translation, cell proliferation and differentiation, and anabolism, which can affect both growth and feed intake. However, the role of TOR signaling in the regulation of feed intake and feed efficiency in poultry is not clear. In the present study, a total of 1000 ducks, of similar initial weight, were chosen and transferred to individual cages to determine their residual feed intake (RFI) from the age of 21 to 42 days. Subsequently, 60 ducks, which were divided into high (HRFI) and low (LRFI) groups according to their RFI, were chosen to analyze the TOR signaling activities in the liver. The differential expression level of genes involved in the TOR signaling pathway was assayed by the real-time polymerase chain reaction. In the liver, the expression of AKT, avTOR, avLST8, and S6K1 was significantly higher in LRFI ducks than in HRFI ducks; avTOR and AKT were negatively associated with the feed conversion ratio and RFI. Furthermore, PI3K was moderately positively associated with AKT; AKT was strongly positively associated with PI3K, avTOR, avLST8, and S6K1; and avTOR was strongly positively associated with S6K1. In conclusion, the activation of avTOR signaling in the liver of LRFI ducks might be ascribed to higher energy state or more active nutrient transport (amino acids), or both, than those in the liver of HRFI ducks. The results of the present study indicate that AKT and avTOR of TOR signaling might be used as candidate genes to assess molecular regulation of feed efficiency.

Application of omics technologies for a deeper insight into quali-quantitative production traits in broiler chickens: A review

The poultry industry is continuously facing substantial and different challenges such as the increasing cost of feed ingredients, the European Union's ban of antibiotic as growth promoters, the antimicrobial resistance and the high incidence of muscle myopathies and breast meat abnormalities. In the last decade, there has been an extraordinary development of many genomic techniques able to describe global variation of genes, proteins and metabolites expression level. Proper application of these cutting-edge omics technologies (mainly transcriptomics, proteomics and metabolomics) paves the possibility to understand much useful information about the biological processes and pathways behind different complex traits of chickens. The current review aimed to highlight some important knowledge achieved through the application of omics technologies and proteo-genomics data in the field of feed efficiency, nutrition, meat quality and disease resistance in broiler chickens.

Transcriptional shifts account for divergent resource allocation in feed efficient broiler chickens

Considerable variation in feed efficiency (FE) has been observed in indigenous and selected meat-type chicken populations. Although this variation could be partially linked to extrinsic factors like diet, housing environment and microbiota, it further illustrates the existence of strong molecular mechanisms enabling the differential allocation of resources for various physiological processes. To further deepen the molecular basis of individual allocation capacity in male and female broilers, an RNA-seq experiment was conducted which based on a phenotyped chicken population divergent in FE. Transcriptional differences linked to FE were pronounced in intestinal and muscular tissue sites of male animals. Specifically, signalling pathways of farnesoid X receptor (FXR) and retinoid X receptor (RXR) might contribute to mediate individual FE. The transcriptional profiles suggested ACSBG2 (muscular lipid utilisation), ASBT (intestinal bile salt transport), CLEC2B (natural killer cell activation), HMGCS2 (jejunal, duodenal and muscular ketogenesis), and SCARB1 (jejunal lipid uptake) as potential mediators driving FE. Results indicate that improvements in FE exploit shifts in resource allocation which might occur at the expense of general immune responsiveness in high efficient male chickens. Consequently, to further improve FE traits and to explore causative molecular patterns, effects originating from sex-dimorphism in chickens need to be taken into consideration.

Differences in intestinal size, structure, and function contributing to feed efficiency in broiler chickens reared at geographically distant locations

The contribution of the intestinal tract to differences in residual feed intake (RFI) has been inconclusively studied in chickens so far. It is also not clear if RFI-related differences in intestinal function are similar in chickens raised in different environments. The objective was to investigate differences in nutrient retention, visceral organ size, intestinal morphology, jejunal permeability and expression of genes related to barrier function, and innate immune response in chickens of diverging RFI raised at 2 locations (L1: Austria; L2: UK). The experimental protocol was similar, and the same dietary formulation was fed at the 2 locations. Individual BW and feed intake (FI) of chickens (Cobb 500FF) were recorded from d 7 of life. At 5 wk of life, chickens (L1, n = 157; L2 = 192) were ranked according to their RFI, and low, medium, and high RFI chickens were selected (n = 9/RFI group, sex, and location). RFI values were similar between locations within the same RFI group and increased by 446 and 464 g from low to high RFI in females and males, respectively. Location, but not RFI rank, affected growth, nutrient retention, size of the intestine, and jejunal disaccharidase activity. Chickens from L2 had lower total body weight gain and mucosal enzyme activity but higher nutrient retention and longer intestines than chickens at L1. Parameters determined only at L1 showed increased crypt depth in the duodenum and jejunum and enhanced paracellular permeability in low vs. high RFI females. Jejunal expression of IL1B was lower in low vs. high RFI females at L2, whereas that of TLR4 at L1 and MCT1 at both locations was higher in low vs. high RFI males. Correlation analysis between intestinal parameters and feed efficiency metrics indicated that feed conversion ratio was more correlated to intestinal size and function than was RFI. In conclusion, the rearing environment greatly affected intestinal size and function, thereby contributing to the variation in chicken RFI observed across locations.

Feed conversion ratio, residual feed intake and cholecystokinin type A receptor gene polymorphisms are associated with feed intake and average daily gain in a Chinese local chicken population

Background

The feed conversion ratio (FCR) and residual feed intake (RFI) are common indexes in measuring feed efficiency for livestock. RFI is a feed intake adjusted for requirements for maintenance and production so these two traits are related. Similarly, FCR is related to feed intake and weight gain because it is their ratio. Cholecystokinin type A receptor (CCKAR) plays an important role in animal digestive process. We examined the interplay of these three parameters in a local Chinese chicken population.

Results

The feed intake (FI) and body weights (BW) of 1,841 individuals were monitored on a daily basis from 56 to 105 d of age. There was a strong correlation between RFI and average daily feed intake (ADFI) and a negative correlation between the FCR and daily gain (r_g = − 0.710). Furthermore, we identified 51 single nucleotide polymorphisms (SNPs) in the CCKAR and 4 of these resulted in amino acid mutations. The C334A mutation was specifically associated with FI and the expected feed intake (EFI) (P <  0.01) and significantly associated with the average daily gain (ADG) (P <  0.05). G1290A was significantly associated with FI and EFI (P < 0.05).

Conclusion

FCR is apply to weight selecting, and RFI is more appropriate if the breeding focus is feed intake. And C334A and G1290A of the CCKAR gene can be deemed as candidate markers for feed intake and weight gain.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0261-1) contains supplementary material, which is available to authorized users.

Identifying artificial selection signals in the chicken genome

Identifying the signals of artificial selection can contribute to further shaping economically important traits. Here, a chicken 600k SNP-array was employed to detect the signals of artificial selection using 331 individuals from 9 breeds, including Jingfen (JF), Jinghong (JH), Araucanas (AR), White Leghorn (WL), Pekin-Bantam (PB), Shamo (SH), Gallus-Gallus-Spadiceus (GA), Rheinlander (RH) and Vorwerkhuhn (VO). Per the population genetic structure, 9 breeds were combined into 5 breed-pools, and a 'two-step' strategy was used to reveal the signals of artificial selection. GA, which has little artificial selection, was defined as the reference population, and a total of 204, 155, 305 and 323 potential artificial selection signals were identified in AR_VO, PB, RH_WL and JH_JF, respectively. We also found signals derived from standing and de-novo genetic variations have contributed to adaptive evolution during artificial selection. Further enrichment analysis suggests that the genomic regions of artificial selection signals harbour genes, including THSR, PTHLH and PMCH, responsible for economic traits, such as fertility, growth and immunization. Overall, this study found a series of genes that contribute to the improvement of chicken breeds and revealed the genetic mechanisms of adaptive evolution, which can be used as fundamental information in future chicken functional genomics study.

Integrative approach using liver and duodenum RNA-Seq data identifies candidate genes and pathways associated with feed efficiency in pigs

This study aims identifying candidate genes and pathways associated with feed efficiency (FE) in pigs. Liver and duodenum transcriptomes of 37 gilts showing high and low residual feed intake (RFI) were analysed by RNA-Seq. Gene expression data was explored through differential expression (DE) and weighted gene co-expression network analyses. DE analysis revealed 55 and 112 differentially regulated genes in liver and duodenum tissues, respectively. Clustering genes according to their connectivity resulted in 23 (liver) and 25 (duodenum) modules of genes with a co-expression pattern. Four modules, one in liver (with 444 co-expressed genes) and three in duodenum (gathering 37, 126 and 41 co-expressed genes), were significantly associated with FE indicators. Intra-module analyses revealed tissue-specific candidate genes; 12 of these genes were also identified as DE between individuals with high and low RFI. Pathways enriched by the list of genes showing DE and/or belonging to FE co-expressed modules included response to oxidative stress, inflammation, immune response, lipid metabolism and thermoregulation. Low overlapping between genes identified in duodenum and liver tissues was observed but heat shock proteins were associated to FE in both tissues. Our results suggest tissue-specific rather than common transcriptome regulatory processes associated with FE in pigs.

Assessing serum metabolite profiles as predictors for feed efficiency in broiler chickens reared at geographically distant locations

1. The objective of this study was to investigate differences in growth performance, serum intermediary metabolites, acute-phase proteins and white blood cells in low, medium and high-residual feed intake (RFI) chickens. It was also assessed if the environment affects the feed efficiency (FE) and FE-related performance and serum profiles of chickens. 2. Individual body weight (BW) and feed intake (FI) were recorded from d 7 of life. At 5 weeks of age, female and male broiler chickens (Cobb 500) were selected according to their RFI (L1: Austria; L2: UK; n = 9/RFI group, sex and locatity -45on) and blood samples were collected. 3. Chickens at L1 had similar FI but a 15% higher BW gain compared to chickens at L2. The RFI values of female chickens were -231, 8 and 215 g and those of male chickens -197, 0 and 267 g for low, medium and high RFI, respectively. 4. Location affected serum glucose, urea, cholesterol, non-esterified fatty acids (NEFA) and ovotransferrin in females, and serum glucose and triglycerides in male chickens. Serum uric acid and NEFA linearly increased from low to high RFI in females, whereas in males, cholesterol showed the same linear response from low to high RFI. Serum alpha-1-acid glycoprotein and blood heterophil-to-lymphocyte ratio linearly increased by 35% and 68%, respectively, from low to high RFI but only in male chickens at L1. 5. Regression analysis showed significant positive relationships between RFI and serum uric acid (R² = 0.49) and cholesterol (R² = 0.13). 6. It was concluded that RFI-related variation in serum metabolites of chickens was largely similar for the two environments and that serum metabolite patterns could be used to predict RFI in chickens.

RNA sequencing for global gene expression associated with muscle growth in a single male modern broiler line compared to a foundational Barred Plymouth Rock chicken line

Background

Modern broiler chickens exhibit very rapid growth and high feed efficiency compared to unselected chicken breeds. The improved production efficiency in modern broiler chickens was achieved by the intensive genetic selection for meat production. This study was designed to investigate the genetic alterations accumulated in modern broiler breeder lines during selective breeding conducted over several decades.

Methods

To identify genes important in determining muscle growth and feed efficiency in broilers, RNA sequencing (RNAseq) was conducted with breast muscle in modern pedigree male (PeM) broilers (n = 6 per group), and with an unselected foundation broiler line (Barred Plymouth Rock; BPR). The RNAseq analysis was carried out using Ilumina Hiseq (2 x 100 bp paired end read) and raw reads were assembled with the galgal4 reference chicken genome. With normalized RPM values, genes showing >10 average read counts were chosen and genes showing <0.05 p-value and >1.3 fold change were considered as differentially expressed (DE) between PeM and BPR. DE genes were subjected to Ingenuity Pathway Analysis (IPA) for bioinformatic functional interpretation.

Results

The results indicate that 2,464 DE genes were identified in the comparison between PeM and BPR. Interestingly, the expression of genes encoding mitochondrial proteins in chicken are significantly biased towards the BPR group, suggesting a lowered mitochondrial content in PeM chicken muscles compared to BPR chicken. This result is inconsistent with more slow muscle fibers bearing a lower mitochondrial content in the PeM. The molecular, cellular and physiological functions of DE genes in the comparison between PeM and BPR include organismal injury, carbohydrate metabolism, cell growth/proliferation, and skeletal muscle system development, indicating that cellular mechanisms in modern broiler lines are tightly associated with rapid growth and differential muscle fiber contents compared to the unselected BPR line. Particularly, PDGF (platelet derived growth factor) signaling and NFE2L2 (nuclear factor, erythroid 2-like 2; also known as NRF2) mediated oxidative stress response pathways appear to be activated in modern broiler compared to the foundational BPR line. Upstream and network analyses revealed that the MSTN (myostatin) –FST (follistatin) interactions and inhibition of AR (androgen receptor) were predicted to be effective regulatory factors for DE genes in modern broiler line. PRKAG3 (protein kinase, AMP-activated, gamma 3 non-catalytic subunit) and LIPE (lipase E) are predicted as core regulatory factors for myogenic development, nutrient and lipid metabolism.

Conclusion

The highly upregulated genes in PeM may represent phenotypes of subclinical myopathy commonly observed in the commercial broiler breast tissue, that can lead to muscle hardening, named as woody breast. By investigating global gene expression in a highly selected pedigree broiler line and a foundational breed (Barred Plymouth Rock), the results provide insight into cellular mechanisms that regulate muscle growth, fiber composition and feed efficiency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3471-y) contains supplementary material, which is available to authorized users.

Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus carpio)

Muscle fat content is an important phenotypic trait in fish, as it affects the nutritional, technical and sensory qualities of flesh. To identify loci and candidate genes associated with muscle fat content and abdominal fat traits, we performed a genome-wide association study (GWAS) using the common carp 250 K SNP assay in a common carp F₂ resource population. A total of 18 loci surpassing the genome-wide suggestive significance level were detected for 4 traits: fat content in dorsal muscle (MFdo), fat content in abdominal muscle (MFab), abdominal fat weight (AbFW), and AbFW as a percentage of eviscerated weight (AbFP). Among them, one SNP (carp089419) affecting both AbFW and AbFP reached the genome-wide significance level. Ten of those loci were harbored in or near known genes. Furthermore, relative expressions of 5 genes related to MFdo were compared using dorsal muscle samples with high and low phenotypic values. The results showed that 4 genes were differentially expressed between the high and low phenotypic groups. These genes are, therefore, prospective candidate genes for muscle fat content: ankyrin repeat domain 10a (ankrd10a), tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (tanc2), and four jointed box 1 (fjx1) and choline kinase alpha (chka). These results offer valuable insights into the complex genetic basis of fat metabolism and deposition.

Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria

As feed represents 60 to 70% of the cost of raising an animal to market weight, feed efficiency (the amount of dry weight intake to amount of wet weight gain) remains an important genetic trait in animal agriculture. To gain greater understanding of cellular mechanisms of feed efficiency (FE), shotgun proteomics was conducted using in-gel trypsin digestion and tandem mass spectrometry on breast muscle samples obtained from pedigree male (PedM) broilers exhibiting high feed efficiency (FE) or low FE phenotypes (n = 4 per group). The high FE group had greater body weight gain (P = 0.004) but consumed the same amount of feed (P = 0.30) from 6 to 7 wk resulting in higher FE (P < 0.001). Over 1800 proteins were identified, of which 152 were different (P < 0.05) by at least 1.3 fold and ≤ 15 fold between the high and low FE phenotypes. Data were analyzed for a modified differential expression (DE) metric (Phenotypic Impact Factors or PIF) and interpretation of protein expression data facilitated using the Ingenuity Pathway Analysis (IPA) program. In the entire data set, 228 mitochondrial proteins were identified whose collective expression indicates a higher mitochondrial expression in the high FE phenotype (binomial probability P < 0.00001). Within the top up and down 5% PIF molecules in the dataset, there were 15 mitoproteome proteins up-regulated and only 5 down-regulated in the high FE phenotype. Pathway enrichment analysis also identified mitochondrial dysfunction and oxidative phosphorylation as the number 1 and 5 differentially expressed canonical pathways (up-regulated in high FE) in the proteomic dataset. Upstream analysis (based on DE of downstream molecules) predicted that insulin receptor, insulin like growth receptor 1, nuclear factor, erythroid 2-like 2, AMP activated protein kinase (α subunit), progesterone and triiodothyronine would be activated in the high FE phenotype whereas rapamycin independent companion of target of rapamycin, mitogen activated protein kinase 4, and serum response factor would be inhibited in the high FE phenotype. The results provide additional insight into the fundamental molecular landscape of feed efficiency in breast muscle of broilers as well as further support for a role of mitochondria in the phenotypic expression of FE.

Funding provided by USDA-NIFA (#2013–01953), Arkansas Biosciences Institute (Little Rock, AR), McMaster Fellowship (AUS to WB) and the Agricultural Experiment Station (Univ. of Arkansas, Fayetteville).

QTL variations for growth-related traits in eight distinct families of common carp (Cyprinus carpio)

BACKGROUND

Comparing QTL analyses of multiple pair-mating families can provide a better understanding of important allelic variations and distributions. However, most QTL mapping studies in common carp have been based on analyses of individual families. In order to improve our understanding of heredity and variation of QTLs in different families and identify important QTLs, we performed QTL analysis of growth-related traits in multiple segregating families.

RESULTS

We completed a genome scan for QTLs that affect body weight (BW), total length (TL), and body thickness (BT) of 522 individuals from eight full-sib families using 250 microsatellites evenly distributed across 50 chromosomes. Sib-pair and half-sib model mapping identified 165 QTLs on 30 linkage groups. Among them, 10 (genome-wide P <0.01 or P < 0.05) and 28 (chromosome-wide P < 0.01) QTLs exhibited significant evidence of linkage, while the remaining 127 exhibited a suggestive effect on the above three traits at a chromosome-wide (P < 0.05) level. Multiple QTLs obtained from different families affect BW, TL, and BT and locate at close or identical positions. It suggests that same genetic factors may control variability in these traits. Furthermore, the results of the comparative QTL analysis of multiple families showed that one QTL was common in four of the eight families, nine QTLs were detected in three of the eight families, and 26 QTLs were found common to two of the eight families. These common QTLs are valuable candidates in marker-assisted selection.

CONCLUSION

A large number of QTLs were detected in the common carp genome and associated with growth-related traits. Some of the QTLs of different growth-related traits were identified at similar chromosomal regions, suggesting a role for pleiotropy and/or tight linkage and demonstrating a common genetic basis of growth trait variations. The results have set up an example for comparing QTLs in common carp and provided insights into variations in the identified QTLs affecting body growth. Discovery of these common QTLs between families and growth-related traits represents an important step towards understanding of quantitative genetic variation in common carp.

Post-weaning blood transcriptomic differences between Yorkshire pigs divergently selected for residual feed intake

Background

Improving feed efficiency (FE) of pigs by genetic selection is of economic and environmental significance. An increasingly accepted measure of feed efficiency is residual feed intake (RFI). Currently, the molecular mechanisms underlying RFI are largely unknown. Additionally, to incorporate RFI into animal breeding programs, feed intake must be recorded on individual pigs, which is costly and time-consuming. Thus, convenient and predictive biomarkers for RFI that can be measured at an early age are greatly desired. In this study, we aimed to explore whether differences exist in the global gene expression profiles of peripheral blood of 35 to 42 day-old pigs with extremely low (more efficient) and high RFI (less efficient) values from two lines that were divergently selected for RFI during the grow-finish phase, to use such information to explore the potential molecular basis of RFI differences, and to initiate development of predictive biomarkers for RFI.

Results

We identified 1972 differentially expressed genes (DEGs) (q ≤ 0.15) between the low (n = 15) and high (n = 16) RFI groups of animals by using RNA sequencing technology. We validated 24 of 37 selected DEGs by reverse transcription-quantitative PCR (RT-qPCR) in a joint analysis of 24 (12 per line) of the 31 samples already used for RNA-seq plus 24 (12 per line) novel samples from the same contemporary group of pigs. Using an analysis of the 24 novel samples alone, only nine of the 37 selected DEGs were validated. Genes involved in small molecule biosynthetic process, antigen processing and presentation of peptide antigen via major histocompatibility complex (MHC) class I, and steroid biosynthetic process were overrepresented among DEGs that had higher expression in the low versus high RFI animals. Genes known to function in the proteasome complex or mitochondrion were also significantly enriched among genes with higher expression in the low versus high RFI animals. Alternatively, genes involved in signal transduction, bone mineralization and regulation of phosphorylation were overrepresented among DEGs with lower expression in the low versus high RFI animals. The DEGs significantly overlapped with genes associated with disease, including hyperphagia, eating disorders and mitochondrial diseases (q < 1E-05). A weighted gene co-expression network analysis (WGCNA) identified four co-expression modules that were differentially expressed between the low and high RFI groups. Genes involved in lipid metabolism, regulation of bone mineralization, cellular immunity and response to stimulus were overrepresented within the two modules that were most significantly differentially expressed between the low and high RFI groups. We also found five of the DEGs and one of the co-expression modules were significantly associated with the RFI phenotype of individual animals (q < 0.05).

Conclusions

The post-weaning blood transcriptome was clearly different between the low and high RFI groups. The identified DEGs suggested potential differences in mitochondrial and proteasomal activities, small molecule biosynthetic process, and signal transduction between the two RFI groups and provided potential new insights into the molecular basis of RFI in pigs, although the observed relationship between the post-weaning blood gene expression and RFI phenotype measured during the grow-finish phase was not strong. DEGs and representative genes in co-expression modules that were associated with RFI phenotype provide a preliminary list for developing predictive biomarkers for RFI in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2395-x) contains supplementary material, which is available to authorized users.

Gene map of large yellow croaker (Larimichthys crocea) provides insights into teleost genome evolution and conserved regions associated with growth

The genetic map of a species is essential for its whole genome assembly and can be applied to the mapping of important traits. In this study, we performed RNA-seq for a family of large yellow croakers (Larimichthys crocea) and constructed a high-density genetic map. In this map, 24 linkage groups comprised 3,448 polymorphic SNP markers. Approximately 72.4% (2,495) of the markers were located in protein-coding regions. Comparison of the croaker genome with those of five model fish species revealed that the croaker genome structure was closer to that of the medaka than to the remaining four genomes. Because the medaka genome preserves the teleost ancestral karyotype, this result indicated that the croaker genome might also maintain the teleost ancestral genome structure. The analysis also revealed different genome rearrangements across teleosts. QTL mapping and association analysis consistently identified growth-related QTL regions and associated genes. Orthologs of the associated genes in other species were demonstrated to regulate development, indicating that these genes might regulate development and growth in croaker. This gene map will enable us to construct the croaker genome for comparative studies and to provide an important resource for selective breeding of croaker.