A polymorphism in the 3'-flanking region of insulin-like growth factor binding protein 2 gene associated with abdominal fat in chickens

Comprehensive analysis of changes in expression of lncRNA, microRNA and mRNA in liver tissues of chickens with high or low abdominal fat deposition

1. The liver of chickens is a dominant lipid biosynthetic tissue and plays a vital role in fat deposition, particularly in the abdomen. To determine the molecular mechanisms involved in its lipid metabolism, the livers of chickens with high (H) or low (L) abdominal fat content were sampled and sequencing on long non-coding RNA (lncRNA), messenger RNA (mRNA) and small RNA (microRNA) was performed.2. In total, 351 expressed protein-coding genes for long non-coding RNA (DEL; 201 upregulated and 150 downregulated), 400 differentially expressed genes (DEG; 223 upregulated and 177 downregulated) and 10 differentially expressed miRNA (DEM; four upregulated and six downregulated) were identified between the two groups. Multiple potential signalling pathways related to lipogenesis and lipid metabolism were identified via pathway enrichment analysis. In addition, 173 lncRNA - miRNA - mRNA interaction regulatory networks were identified, including 30 lncRNA, 27 mRNA and seven miRNA.3. These networks may help regulate lipid metabolism and fat deposition. Five promising candidate genes and two lncRNA may play important roles in the regulation of adipogenesis and lipid metabolism in chickens.

Genome-wide association study of growth curve parameters reveals novel genomic regions and candidate genes associated with metatarsal bone traits in chickens

The growth and development of chicken bones have an enormous impact on the health and production performance of chickens. However, the development pattern and genetic regulation of the chicken skeleton are poorly understood. This study aimed to evaluate metatarsal bone growth and development patterns in chickens via non-linear models, and to identify the genetic determinants of metatarsal bone traits using a genome-wide association study (GWAS) based on growth curve parameters. Data on metatarsal length (MeL) and metatarsal circumference (MeC) were obtained from 471 F2 chickens (generated by crossing broiler sires, derived from a line selected for high abdominal fat, with Baier layer dams) at 4, 6, 8, 10, and 12 weeks of age. Four non-linear models (Gompertz, Logistic, von Bertalanffy, and Brody) were used to fit the MeL and MeC growth curves. Subsequently, the estimated growth curve parameters of the mature MeL or MeC (A), time-scale parameter (b), and maturity rate (K) from the non-linear models were utilized as substitutes for the original bone data in GWAS. The Logistic and Brody models displayed the best goodness-of-fit for MeL and MeC, respectively. Single-trait and multi-trait GWASs based on the growth curve parameters of the Logistic and Brody models revealed 4 618 significant single nucleotide polymorphisms (SNPs), annotated to 332 genes, associated with metatarsal bone traits. The majority of these significant SNPs were located on Gallus gallus chromosome (GGA) 1 (167.433-176.318 Mb), GGA2 (96.791-103.543 Mb), GGA4 (65.003-83.104 Mb) and GGA6 (64.685-95.285 Mb). Notably, we identified 12 novel GWAS loci associated with chicken metatarsal bone traits, encompassing 35 candidate genes. In summary, the combination of single-trait and multi-trait GWASs based on growth curve parameters uncovered numerous genomic regions and candidate genes associated with chicken bone traits. The findings benefit an in-depth understanding of the genetic architecture underlying metatarsal growth and development in chickens.

A Functional Variant Alters the Binding of Bone morphogenetic protein 2 to the Transcription Factor NF-κB to Regulate Bone morphogenetic protein 2 Gene Expression and Chicken Abdominal Fat Deposition

In this study, we employed a dual-luciferase reporter assay and electrophoretic mobility shift analysis (EMSA) in vitro to explore whether a 12-base pair (bp) insertion/deletion (InDel) variant (namely g.14798187_14798188insTCCCTGCCCCCT) within intron 2 of the chicken BMP2 gene, which was significantly associated with chicken abdominal fat weight and abdominal fat percentage, is a functional marker and its potential regulatory mechanism. The reporter analysis demonstrated that the luciferase activity of the deletion allele was extremely significantly higher than that of the insertion allele (p < 0.01). A bioinformatics analysis revealed that compared to the deletion allele, the insertion allele created a transcription factor binding site of nuclear factor-kappa B (NF-κB), which exhibited an inhibitory effect on fat deposition. A dual-luciferase reporter assay demonstrated that the inhibitory effect of NF-κB on the deletion allele was stronger than that on the insertion allele. EMSA indicated that the binding affinity of NF-κB for the insertion allele was stronger than that for the deletion allele. In conclusion, the 12-bp InDel chicken BMP2 gene variant is a functional variant affecting fat deposition in chickens, which may partially regulate BMP2 gene expression by affecting the binding of transcription factor NF-κB to the BMP2 gene.

Co-expression analysis of long non-coding RNAs and mRNAs involved in intramuscular fat deposition in Muchuan black-bone chicken

The intramuscular fat (IMF) content in meat products is positively correlated with meat quality, making it an important consumer trait. Long non-coding RNAs (lncRNAs) play central roles in regulating various biological processes, but little is currently known about the mechanisms by which they regulate IMF deposition in chickens. This study sampled the breast muscles of chickens with high (H) and low (L) IMF content and constructed six small RNA libraries. High-throughput sequencing technology was used to profile the breast muscle transcriptome (lncRNA and mRNA) and to identify the differentially expressed lncRNAs (DELs) and mRNAs (DEGs) between the H and L groups. In total, 263 DELs (118 up-regulated and 145 down-regulated lncRNAs) and 443 DEGs (203 up-regulated and 240 down-regulated genes) were identified between the two groups. To analyse the DELs-DEGs interaction network, co-expression analysis was conducted to identify lncRNA-mRNA pairs. In total, 19,270 lncRNA/mRNA pairs were identified, including 16,398 significant correlation pairs that presented as positive and 2872 pairs that presented as negative. The lncRNA-mRNA network comprised 263 lncRNA nodes and 440 mRNA nodes. Pathway analysis, using the Kyoto Encyclopedia of Genes and Genomes, indicated that pathways associated with fat deposition and lipid metabolism such as the MAPK, PPAR, GnRH, ErbB and calcium signalling pathways, fatty acid elongation and fatty acid metabolism. Overall, the study identified potential candidate lncRNAs, genes and regulatory networks associated with chicken IMF deposition. These findings provide new insights to help clarify the regulatory mechanisms of IMF deposition in chickens which can be used to improve the IMF content in poultry.

A combination of genome-wide association study and selection signature analysis dissects the genetic architecture underlying bone traits in chickens

The bones of chicken play an important role in supporting and protecting the body. The growth and development of bones have a substantial influence on the health and production performance in chickens. However, genetic architecture underlying chicken bone traits is not well understood. The objectives of this study are to dissect the genetic basis of bone traits in chickens and to identify valuable genes and genetic markers for chicken breeding. We performed a combination of genome-wide association study (GWAS) and selection signature analysis (fixation index values and nucleotide diversity ratios) in an F₂ crossbred experimental population with different genetic backgrounds (broiler × layer) to identify candidate genes and significant variants related to femur, shank, keel length, chest width, metatarsal claw weight, metatarsal length, and metatarsal circumference. A total of 545 individuals were genotyped based on the whole genome re-sequencing method (26 F₀ individuals were re-sequenced at 10 × coverage; 519 F₂ individuals were re-sequenced at 3 × coverage). A total of 2 028 112 single-nucleotide polymorphisms (SNPs) remained to carry out analysis after quality control and imputation. The integration of GWAS and selection signature analysis indicated that all significant SNPs responsible for bone traits were mainly localized on chicken chromosomes 1, 4, and 27. Finally, we identified 21 positional candidate genes that might regulate chicken bone growth and development, including LRCH1, RB1, FNDC3A, MLNR, CAB39L, FOXO1, LHFP, TRPC4, POSTN, SMAD9, RBPJ, PPARGC1A, SLIT2, NCAPG, NKX3-2, CPZ, SPOP, NGFR, SOST, ZNF652, and HOXB3. Additionally, an array of uncharacterized genes was identified. The findings provide an in-depth understanding of the genetic architecture of chicken bone traits and offer a molecular basis for applying genomics in practical chicken breeding.

Chicken C/EBPζ gene: Expression profiles, association analysis, and identification of functional variants for abdominal fat

CCAAT enhancer binding protein ζ (C/EBPζ) plays an important role in adipose proliferation and differentiation in humans. However, very little is known about the effect of C/EBPζ on the growth and development of adipose tissues in domesticated animals. The present study attempted to investigate the mRNA expression profiles of chicken C/EBPζ in a variety of tissues; analyze the association of its variants with abdominal fat; and identify the functional variants for abdominal fat. The tissue expression profiles revealed that C/EBPζ was highly expressed in 19 tissues obtained from broilers. The expression level of C/EBPζ in fat broilers was significantly lower than that in lean broilers in the duodenum, ileum, cecum, kidney, pectoral muscle, and liver (P < 0.05). Among 170 polymorphic loci of C/EBPζ, 9 single nucleotide polymorphisms (SNPs) demonstrated a significant association with chicken abdominal fat traits (P < 0.05) as well as significant discrepancies in their allelic frequencies between fat and lean birds. Particularly, only C/EBPζ g.7085A>C exhibited significant correlation with abdominal fat traits (P < 0.00015) using the Bonferroni method. The results revealed that, in preadipocyte immortalized cells (ICPI), the luciferase activity of the A allele of g.7085A>C locus was remarkably stronger than that of the C allele (P < 0.05). In silico analysis showed that g.7085A>C locus was located in the binding region of the transcription factor SOX5, which possesses the ability to transform C/EBPζ transcription efficiency through binding with SOX5. In summary, the data obtained from this study suggested that C/EBPζ is a potential candidate gene responsible for abdominal fat deposition in chicken and that g.7085A>C is a functional SNP that can be promisingly leveraged for marker assisted selection (MAS) in future chicken breeding programs.

A functional variant in the promoter region of IGF1 gene is associated with chicken abdominal fat deposition

Insulin-like growth factor 1 (IGF1) plays an important role in the regulation of cell growth, proliferation, differentiation, and apoptosis. Previously several studies revealed that genotypes of chicken IGF1 c.-366A > C were significantly associated with abdominal fat weight and body weight in chickens. But the underlying mechanism is still unknown. To investigate the mechanism underlying the association, herein, we performed IGF1 gene mRNA expression profiling, a dual-luciferase reporter assay and electrophoretic mobility shift assay (EMSA). Quantitative real-time PCR results showed that IGF1 gene was widely expressed in 14 tissues. The mRNA expression levels of IGF1 gene in both abdominal fat and jejunum were significantly higher in fat broilers than in lean broilers. However, the opposite results were observed in the pancreas. The reporter gene assay showed that the promoter luciferase activity of allele A was significantly higher than that of allele C (P < 0.05). In addition, the luciferase activity of allele A promoted by the transcription factor AP1 and OCT1 was higher than that of allele C (P < 0.05). Electrophoretic mobility shift assay result showed that allele A binding to the transcription factor AP1 and OCT1 was stronger than that of allele C. All in all, our data indicated that the IGF1 gene c.-366A > C is a functional SNP responsible for chicken adipose deposition.

Profiling of RNA N 6 -Methyladenosine Methylation Reveals the Critical Role of m6A in Chicken Adipose Deposition

One of the main objectives of broiler breeding is to prevent excessive abdominal adipose deposition. The role of RNA modification in adipose deposition is not clear. This study was aimed to map m⁶A modification landscape in chicken adipose tissue. MeRIP-seq was performed to compare the differences in m⁶A methylation pattern between fat and lean broilers. We found that start codons, stop codons, coding regions, and 3′-untranslated regions were generally enriched for m⁶A peaks. The high m⁶A methylated genes (fat birds vs. lean birds) were primarily associated with fatty acid biosynthesis and fatty acid metabolism, while the low m⁶A methylated genes were mainly involved in processes associated with development. Furthermore, we found that the mRNA levels of many genes may be regulated by m⁶A modification. This is the first comprehensive characterization of m⁶A patterns in the chicken adipose transcriptome, and provides a basis for studying the role of m⁶A modification in fat deposition.

Polymorphisms of KLF3 gene coding region and identification of their functionality for abdominal fat in chickens

KLF3 is a member of the Kruppel‐like factor (KLF) family of transcription factors, and plays an important role in several biological processes, including adipogenesis, erythropoiesis and B‐cell development. The purposes of this study are to search for polymorphisms of KLF3 coding region and to provide functional evidence for abdominal fat in chickens. A total of 168 SNPs in KLF3 coding region were detected in a unique chicken population, the Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF). Of which three single nucleotide polymorphisms (g.3452T > C, g.8663A > G and g.10751G > A) were significantly correlated with abdominal fat weight (AFW) and abdominal fat percentage (AFP) of 329 birds from the 19th generation of NEAUHLF (FDR < 0.05). The reporter gene assay was performed to verify functionality of these three SNPs in both ICP‐1 and DF1 cells. Results showed that the luciferase activity of G allele was significantly higher than that of A allele in g.10751G > A (p < 0.05). However, there were no significant differences between different alleles of others two SNPs in luciferase activity. Overall, KLF3 is an important candidate gene that affects chicken abdominal fat content, and the g.10751G > A is a functional variant that potential would be applied to marker‐assisted selection (MAS) for selective breeding programme.

Associations between IGF1, IGFBP2 and TGFß3 Genes Polymorphisms and Growth Performance of Broiler Chicken Lines

Marker-assisted selection based on fast and accurate molecular analysis of individual genes is considered an acceptable tool in the speed-up of the genetic improvement of production performance in chickens. The objective of this study was to detect the single nucleotide polymorphisms (SNPs) in the IGF1, IGFBP2 and TGFß3 genes, and to investigate their associations with growth performance (body weight (BW) and average daily gain (ADG) at 14, 21, 28, 35 and 42 days of age) and carcass traits in broilers. Performance (carcass) data (weight before slaughter; weights of the trunk, giblets, abdominal fat, breast muscle and thigh muscle; slaughter value and slaughter percentage), as well as blood samples for DNA extraction and SNP analysis, were obtained from 97 chickens belonging to two different lines (Hubbard F15 and Cobb E) equally divided between the two sexes. The genotypes were detected using polymerase chain reaction- restriction fragment length polymorphism (PCR-RFLP) methods with specific primers and restrictase for each gene. The statistical analysis discovered significant associations (p < 0.05) between the TGFβ3 SNP and the following parameters: BW at 21, 28 and 35 days, trunk weight and slaughter value. Association analysis of BWs (at 21, 28 and 35 days) and SNPs was always significant for codominant, dominant and overdominant genetic models, showing a possible path for genomic selection in these chicken lines. Slaughter value was significant for codominant, recessive and overdominant patterns, whereas other carcass traits were not influenced by SNPs. Based on the results of this study, we suggested that the TGFβ3 gene could be used as a candidate gene marker for chicken growth traits in the Hubbard F15 and Cobb E population selection programs, whereas for carcass traits further investigation is needed.

Effect of IGFBP2 Overexpression on the Expression of Fatty Acid Synthesis Genes in Primary Cultured Chicken Hepatocytes

The effects of insulin-like growth factor binding protein 2 (IGFBP2) on the expression of fatty acid synthesis regulators and triglyceride production were investigated in primary cultured chicken hepatocytes. The full-length chicken IGFBP2 coding region was synthesized by overlap extension PCR and cloned into the pcDNA3.1 vector. An in situ digestion method was used to prepare the chicken hepatocytes. Primary chicken hepatocytes were maintained in monolayer culture. Real-time PCR was used to detect changes in the expression of IGFBP2, PPARG, IGF1, IGF1R, APOAI, and LFABP, after the overexpression of IGFBP2 in chicken hepatocytes. Triglyceride production and glucose content were also evaluated using triglyceride and glucose analysis methods. The expression level of IGFBP2 increased after transfection of the IGFBP2-containing vector. The expression levels of PPARG, IGF1, and IGF1R also increased in cultured chicken hepatocytes after the overexpression of IGFBP2, whereas the expression of LFABP and APOAI decreased. Triglyceride production in primary cultured chicken hepatocytes increased after the overexpression of IGFBP2. These results suggest that IGFBP2 is involved in lipogenesis, increasing both the expression of fatty acid synthesis regulators, and triglyceride production in primary cultured chicken hepatocytes.

Identification of selection signatures involved in performance traits in a paternal broiler line

BACKGROUND

Natural and artificial selection leads to changes in certain regions of the genome resulting in selection signatures that can reveal genes associated with the selected traits. Selection signatures may be identified using different methodologies, of which some are based on detecting contiguous sequences of homozygous identical-by-descent haplotypes, called runs of homozygosity (ROH), or estimating fixation index (FST) of genomic windows that indicates genetic differentiation. This study aimed to identify selection signatures in a paternal broiler TT line at generations 7th and 16th of selection and to investigate the genes annotated in these regions as well as the biological pathways involved. For such purpose, ROH and FST-based analysis were performed using whole genome sequence of twenty-eight chickens from two different generations.

RESULTS

ROH analysis identified homozygous regions of short and moderate size. Analysis of ROH patterns revealed regions commonly shared among animals and changes in ROH abundance and size between the two generations. Results also suggest that whole genome sequencing (WGS) outperforms SNPchip data avoiding overestimation of ROH size and underestimation of ROH number; however, sequencing costs can limited the number of animals analyzed. FST-based analysis revealed genetic differentiation in several genomic windows. Annotation of the consensus regions of ROH and FST windows revealed new and previously identified genes associated with traits of economic interest, such as APOB, IGF1, IGFBP2, POMC, PPARG, and ZNF423. Over-representation analysis of the genes resulted in biological terms of skeletal muscle, matrilin proteins, adipose tissue, hyperglycemia, diabetes, Salmonella infections and tyrosine.

CONCLUSIONS

Identification of ROH and FST-based analyses revealed selection signatures in TT line and genes that have important role in traits of economic interest. Changes in the genome of the chickens were observed between the 7th and 16th generations showing that ancient and recent selection in TT line may have acted over genomic regions affecting diseases and performance traits.

Genomic Insights Into the Multiple Factors Controlling Abdominal Fat Deposition in a Chicken Model

Genetic selection for an increased growth rate in meat-type chickens has been accompanied by excessive fat accumulation particularly in abdominal cavity. These progressed to indirect and often unhealthy effects on meat quality properties and increased feed cost. Advances in genomics technology over recent years have led to the surprising discoveries that the genome is more complex than previously thought. Studies have identified multiple-genetic factors associated with abdominal fat deposition. Meanwhile, the obesity epidemic has focused attention on adipose tissue and the development of adipocytes. The aim of this review is to summarize the current understanding of genetic/epigenetic factors associated with abdominal fat deposition, or as it relates to the proliferation and differentiation of preadipocytes in chicken. The results discussed here have been identified by different genomic approaches, such as QTL-based studies, the candidate gene approach, epistatic interaction, copy number variation, single-nucleotide polymorphism screening, selection signature analysis, genome-wide association studies, RNA sequencing, and bisulfite sequencing. The studies mentioned in this review have described multiple-genetic factors involved in an abdominal fat deposition. Therefore, it is inevitable to further study the multiple-genetic factors in-depth to develop novel molecular markers or potential targets, which will provide promising applications for reducing abdominal fat deposition in meat-type chicken.

Effects of Fasting and Refeeding on the mRNA levels of Insulin-like Growth Factor-binding Proteins in Chick Liver and Brain

The physiological functions of insulin-like growth factor-binding proteins (IGFBPs) in mammals have been evaluated in several studies. However, the physiological roles of IGFBPs in chickens have not yet been elucidated. In this study, we examined the effects of short-term (6 h) fasting and refeeding on the mRNA levels of IGFBPs in chick liver and brain. Eighteen 8-day-old chicks were weighed and allocated to three groups on the basis of body weight, and subjected to ad libitum feeding, 6 h of fasting, or 6 h of fasting followed by 6 h of refeeding. After the chicks were euthanized by decapitation, the liver and brain were excised, and the brain was dissected into six segments (telencephalon, optic lobes, cerebellum, rostral part of the brainstem, middle part of the brainstem, and caudal part of the brainstem). IGFBP mRNA levels were determined by qRT-PCR. Fasting significantly increased the mRNA levels of IGFBP-1 and -2 in the chick liver, and these changes were reversed by 6 h of refeeding. The mRNA levels of IGFBP-3 in the middle part of the brainstem and IGFBP-5 in the optic lobes were decreased by 6 h of fasting and were not reversed after 6 h of refeeding. These findings suggest that IGFBP-1 and -2 in the liver, IGFBP-3 in the middle part of the brainstem, and IGFBP-5 in the optic lobes may play physiological roles in response to short-term changes in the nutritional status of chicks.

Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines

Background

Meat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection.

Results

A total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development.

Conclusions

In this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4444-0) contains supplementary material, which is available to authorized users.

IGFBP2 enhances adipogenic differentiation potentials of mesenchymal stem cells from Wharton's jelly of the umbilical cord via JNK and Akt signaling pathways

Mesenchymal stem cell (MSC)-mediated tissue engineering represents a promising strategy to address adipose tissue defects. MSCs derived from Wharton’s jelly of the umbilical cord (WJCMSCs) may serve as an ideal source for adipose tissue engineering due to their abundance, safety profile, and accessibility. How to activate the directed differentiation potentials of WJCMSCs is the core point for their clinical applications. A thorough investigation of mechanisms involved in WJCMSC adipogenic differentiation is necessary to support their application in adipose tissue engineering and address shortcomings. Previous study showed, compared with periodontal ligament stem cells (PDLSCs), WJCMSCs had a weakened adipogenic differentiation potentials and lower expression of insulin-like growth factor binding protein 2 (IGFBP2). IGFBP2 may be involved in the adipogenesis of MSCs. Generally, IGFBP2 is involved in regulating biological activity of insulin-like growth factors, however, its functions in human MSCs are unclear. Here, we found IGFBP2 expression was upregulated upon adipogenic induction, and that IGFBP2 enhanced adipogenic differentiation of WJCMSCs and BMSCs. Moreover, IGFBP2 increased phosphorylation of c-Jun N-terminal kinase (p-JNK) and p-Akt, and activated JNK or Akt signaling significantly promoted adipogenic differentiation of MSCs. Furthermore, inhibitor-mediated blockage of either JNK or Akt signaling dramatically reduced IGFBP2-mediated adipogenic differentiation. And the JNK inhibitor, SP600125 markedly blocked IGFBP2-mediated Akt activation. Moreover, IGFBP2 was negatively regulated by BCOR, which inhibited adipogenic differentiation of WJCMSCs. Overall, our results reveal a new function of IGFBP2, providing a novel insight into the mechanism of adipogenic differentiation and identifying a potential target mediator for improving adipose tissue engineering based on WJCMSCs.

Differential fecal microbiota are retained in broiler chicken lines divergently selected for fatness traits

Our study combined 16S rRNA-pyrosequencing and whole genome sequencing to analyze the fecal metagenomes of the divergently selected lean (LL) and fat (FL) line chickens. Significant structural differences existed in both the phylogenic and functional metagenomes between the two chicken lines. At phylum level, the FL group had significantly less Bacteroidetes. At genus level, fourteen genera of different relative abundance were identified, with some known short-chain fatty acid producers (including Subdoligranulum, Butyricicoccus, Eubacterium, Bacteroides, Blautia) and a potentially pathogenic genus (Enterococcus). Redundancy analysis identified 190 key responsive operational taxonomic units (OTUs) that accounted for the structural differences between the phylogenic metagenome of the two groups. Four Cluster of Orthologous Group (COG) categories (Amino acid transport and metabolism, E; Nucleotide transport and metabolism, F; Coenzyme transport and metabolism, H; and Lipid transport and metabolism, I) were overrepresented in LL samples. Fifteen differential metabolic pathways (Biosynthesis of amino acids, Pyruvate metabolism, Nitrotoluene degradation, Lipopolysaccharide biosynthesis, Peptidoglycan biosynthesis, Pantothenate and CoA biosynthesis, Glycosaminoglycan degradation, Thiamine metabolism, Phosphotransferase system, Two-component system, Bacterial secretion system, Flagellar assembly, Bacterial chemotaxis, Ribosome, Sulfur relay system) were identified. Our data highlighted interesting variations between the gut metagenomes of these two chicken lines.

Association between Motilin Receptor Gene Haplotypes and Growth Traits in Japanese Hinai-dori Crossbred Chickens

We previously identified quantitative trait loci (QTL) for body weight and average daily gain in a common region between ADL0198 (chr 1: 171.7 Mb) and ABR0287 (chr 1: 173.4 Mb) on chicken chromosome 1 in an F₂ resource population produced by crossing low- and high-growth lines of the Hinai-dori breed. Motilin receptor (MLNR) is a candidate gene affecting growth traits in the region. In this study, we genotyped polymorphisms of the MLNR gene and investigated its association with growth traits in a Hinai-dori F₂ intercross population. All the exons of the MLNR gene in the parental population were subjected to PCR amplification, nucleotide sequenced and haplotypes identified. To distinguish resultant diplotype individuals in the F₂ population, a mismatch amplification mutation assay was performed. Three haplotypes (Haplotypes 1–3) were accordingly identified. Six genotypes produced by the combination of three haplotypes (Haplotype 1, 2, and 3) were examined in order to identify associations between MLNR haplotypes and growth traits. The data showed that Haplotype 1 was superior to Haplotype 2 and 3 in body weight at 10 and 14 weeks of age, average daily gain between 4 and 10 weeks, 10 and 14 weeks, and 0 and 14 weeks of age in female in F₂ females. It was concluded that MLNR is a useful marker of growth traits and could be used to develop strategies for improving growth traits in the Hinai-dori breed.

Population genomic analyses based on 1 million SNPs in commercial egg layers

Identifying signatures of selection can provide valuable insight about the genes or genomic regions that are or have been under selective pressure, which can lead to a better understanding of genotype-phenotype relationships. A common strategy for selection signature detection is to compare samples from several populations and search for genomic regions with outstanding genetic differentiation. Wright's fixation index, F_ST, is a useful index for evaluation of genetic differentiation between populations. The aim of this study was to detect selective signatures between different chicken groups based on SNP-wise F_ST calculation. A total of 96 individuals of three commercial layer breeds and 14 non-commercial fancy breeds were genotyped with three different 600K SNP-chips. After filtering a total of 1 million SNPs were available for F_ST calculation. Averages of F_ST values were calculated for overlapping windows. Comparisons of these were then conducted between commercial egg layers and non-commercial fancy breeds, as well as between white egg layers and brown egg layers. Comparing non-commercial and commercial breeds resulted in the detection of 630 selective signatures, while 656 selective signatures were detected in the comparison between the commercial egg-layer breeds. Annotation of selection signature regions revealed various genes corresponding to productions traits, for which layer breeds were selected. Among them were NCOA1, SREBF2 and RALGAPA1 associated with reproductive traits, broodiness and egg production. Furthermore, several of the detected genes were associated with growth and carcass traits, including POMC, PRKAB2, SPP1, IGF2, CAPN1, TGFb2 and IGFBP2. Our approach demonstrates that including different populations with a specific breeding history can provide a unique opportunity for a better understanding of farm animal selection.

Expression of Potential Regulatory Genes in Abdominal Adipose Tissue of Broiler Chickens during Early Development

The identities of genes that underlie population variation in adipose tissue development in farm animals are poorly understood. Previous studies in our laboratory have suggested that increased fat tissue involves the expression modulation of an array of genes in broiler chickens. Of special interest are eight genes, FGFR3, EPHB2, IGFBP2, GREM1, TNC, COL3A1, ACBD7, and SCD. To understand their expression regulation and response to dietary manipulation, we investigated their mRNA levels after dietary manipulation during early development. Chickens were fed either a recommended standard or a high caloric diet from hatch to eight weeks of age (WOA). The high caloric diet markedly affected bodyweight of the broiler birds. mRNA levels of the eight genes in the abdominal adipose tissue were assayed at 2, 4, 6, and 8 WOA using RT-qPCR. Results indicate that (1) FGFR3 mRNA level was affected significantly by diet, age, and diet:age interaction; (2) COL3A mRNA level was repressed by high caloric diet; (3) mRNA levels of EPHB2, ACBD7, and SCD were affected by age; (4) mRNA level of TNC was modulated by age:diet interaction; (5) changes in GREM1 and IGFBP2 mRNA levels were not statistically different.

Lipoprotein metabolism differs between Marek's disease susceptible and resistant chickens

Marek's disease (MD) is a lymphoproliferative disease of chickens caused by MD virus and has an important impact on the poultry industry worldwide. There have been reports showing different physiological characteristics between MD susceptible and resistant chickens. However, little is known about whether there are differences in lipid metabolism between MD susceptible and resistant lines of chickens. In this study, we examined the BW and the weight of tissues (abdominal fat, breast muscle with bone, leg muscle with bone, liver, and heart), the lipoprotein-cholesterol concentrations and distributions, and the plasma and tissue levels of adiponectin and its receptors in the highly resistant and susceptible lines during chicken growth. Our data showed that the increase in total cholesterol during growth was mainly due to the elevation of cholesterol in the low-density/very low-density lipoprotein fraction in MD susceptible chickens, whereas the increase of total cholesterol was mainly attributable to the increase in high-density lipoprotein-cholesterol in MD resistant chickens. Meanwhile, the MD resistant line appeared to have increased plasma adiponectin levels compared with MD susceptible chickens during growth. Taken together, our data suggested that lipoprotein-cholesterol and adiponectin metabolism are different between MD susceptible and resistant chickens.

Association between cholecystokinin type A receptor haplotypes and growth traits in Japanese Hinai-dori crossbred chickens

We previously identified quantitative trait loci for body weight and average daily gain in a common region between MCW0240 (chr 4: 69.9 Mb) and ABR0622 (chr 4: 86.3 Mb) on chicken chromosome 4 in an F(2) resource population produced by crossing low- and high-growth lines of the Hinai-dori breed. Cholecystokinin type A receptor (CCKAR) is a candidate gene affecting growth traits in the region. In this study, we genotyped polymorphisms of the CCKAR gene and investigated its association with growth traits in a Hinai-dori F(2) intercross population. All the exons of the CCKAR gene in the parental population were subjected to PCR amplification, nucleotide sequenced and haplotypes identified. To distinguish resultant diplotype individuals in the F(2) population, a mismatch amplification mutation assay was performed. Five haplotypes (Haplotypes 1-5) were accordingly identified. Six genotypes produced by the combination of three haplotypes (Haplotype 1, 3, and 4) were examined in order to identify associations between CCKAR haplotypes and growth traits. The data indicate that Haplotype 1 was superior to Haplotype 3 and 4 in body weight at 10 and 14 weeks of age, average daily gain between 4 and 10 weeks, 10 and 14 weeks, and 0 and 14 weeks of age. It was concluded that CCKAR is a useful marker of growth traits and could be used to develop strategies for improving growth traits in the Hinai-dori breed.

Genetic epistasis analysis of 10 peroxisome proliferator-activated receptor γ-correlated genes in broiler lines divergently selected for abdominal fat content

Chicken peroxisome proliferator-activated receptor γ (PPARγ), which is highly expressed in adipose tissues, is a key factor in fat accumulation in the abdominal fat pad. In this study, association and pairwise epistasis analyses were performed for all the polymorphisms detected in PPARγ and for 9 genes from PPARγ-correlated lipid metabolic pathways for abdominal fat weight (AFW) in 10th-generation populations of Northeast Agricultural University broiler lines divergently selected for abdominal fat content. Epistatic networks were then reconstructed with the identified epistatic effects. Single-marker association analyses showed that 5 of the 20 screened polymorphisms were significantly associated with AFW (P < 0.05), and CCAAT/enhancer-binding protein α (C/EBPα) c.552G>A was 1 of the 5 significant loci. Pairwise interaction analyses showed that 15 pairs of polymorphisms reached a significance level of P < 2.64 × 10(-4) (adjusted by Bonferroni correction) in the lean line, 41 pairs reached significance in the fat line, and 7 pairs reached significance in both lines. Interestingly, many other loci interacted with C/EBPα c.552G>A in both lines. In epistatic network analyses, C/EBPα c.552G>A seemed to behave as a hub for the epistatic network in both lines. All these results revealed that the genetic architecture of C/EBPα c.552G>A for AFW seemed to be an apparent individual main-effect QTL but that it could be dissected into a genetic epistatic network. Our results suggest that C/EBPα c.552G>A might be the most important locus contributing to phenotypic variation in AFW among all the polymorphisms detected in this study.