Haplotype analysis of genomic prediction by incorporating genomic pathway information based on high-density SNP marker in Chinese yellow-feathered chicken.

Genomic selection using single nucleotide polymorphism (SNP) markers is now intensively investigated in breeding and has been widely utilized for genetic improvement. Currently, several studies have used haplotype (consisting of multiallelic SNPs) for genomic prediction and revealed its performance advantage. In this study, we comprehensively evaluated the performance of haplotype models for genomic prediction in 15 traits, including 6 growth, 5 carcass, and 4 feeding traits in a Chinese yellow-feathered chicken population. We adopted 3 methods to define haplotypes from high-density SNP panels, and our strategy included combining Kyoto Encyclopedia of Genes and Genomes pathway information and considering linkage disequilibrium (LD) information. Our results showed an increase in prediction accuracy due to haplotypes ranging from -0.04∼27.16% in all traits, where the significant improvements were found in 12 traits. The estimates of haplotype epistasis heritability were strongly correlated with the accuracy increase by haplotype models. In addition, incorporating genomic annotation information could further increase the accuracy of the haplotype model, where the further increase in accuracy is significantly relative to the increase of relative haplotype epistasis heritability. The genomic prediction using LD information for constructing haplotypes has the best prediction performance among the 4 traits. These results uncovered that haplotype methods were beneficial for genomic prediction, and the accuracy could be further increased by incorporating genomic annotation information. Moreover, using LD information would potentially improve the performance of genomic prediction.

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

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

Genome-Wide Association Studies Provide Insight Into the Genetic Determination for Hyperpigmentation of the Visceral Peritoneum in Broilers.

Hyperpigmentation of the visceral peritoneum (HVP) has been becoming one of the most challenging problems in yellow-feathered chicken production, which seriously affected chicken carcass quality traits. Detecting which genes dominantly impact pigmentation in the peritoneum tissues is of great benefit to the genetic improvement of HVP. To investigate the genetic mechanism of HVP in yellow-feathered broilers, genome-wide association studies (GWASs) were conducted in the F2 generation of a cross broiler population with 395 birds. A total of 115,706 single-nucleotide polymorphisms (SNPs) of 122,415 were retained to identify quantitative trait loci (QTL) associated to HVP in chicken. The GWAS results based on the logistic mixed model (LMM) revealed that a narrow genomic location on chromosomes 1 (49.2-51.3 Mb) was significantly associated (p ≤ 4.32 × 10-7) with HVP, which contained 23 SNP makers related to 14 functional genes (MFNG, POLDIP3, POLR2F, PICK1, PDXP, SGSM3, RANGAP1, MYH9, RPL3, GALP3, LGALS1, MICALL1, ATF4, and CYP2D6). Four highly associated (p < 10-5) haplotype blocks of 0.80 kb (two SNPs), 0.06 kb (two SNPs), 0.95 kb (two SNPs), and 0.03 kb (two SNPs) were identified with two, two, four, and four distinct haplotypes, respectively. As a melanoma-associated gene, CYP2D6 were also possibly involved in the development of HVP occurring in chicken with two significant variations (rs314284996 and rs317955795) in the promoter regions. Further tests revealed that the expression of CYP2D6 was obviously higher in the visceral peritoneum tissue of chicken with HVP than that in the normal group (p < 0.05). Our results provide a novel clue to understand the genetic mechanism of HVP generation in chicken, and the mapped QTL or candidate genes might serve for genomic selection to improve carcass quality in the yellow-feathered chicken industry.

Transcriptome Response of Liver and Muscle in Heat-Stressed Laying Hens.

Exposure to high ambient temperature has detrimental effects on poultry welfare and production. Although changes in gene expression due to heat exposure have been well described for broiler chickens, knowledge of the effects of heat on laying hens is still relatively limited. In this study, we profiled the transcriptome for pectoralis major muscle (n = 24) and liver (n = 24), during a 4-week cyclic heating experiment performed on layers in the early phase of egg production. Both heat-control and time-based contrasts were analyzed to determine differentially expressed genes (DEGs). Heat exposure induced different changes in gene expression for the two tissues, and we also observed changes in gene expression over time in the control animals suggesting that metabolic changes occurred during the transition from onset of lay to peak egg production. A total of 73 DEGs in liver were shared between the 3 h heat-control contrast, and the 4-week versus 3 h time contrast in the control group, suggesting a core set of genes that is responsible for maintenance of metabolic homeostasis regardless of the physiologic stressor (heat or commencing egg production). The identified DEGs improve our understanding of the layer's response to stressors and may serve as targets for genetic selection in the future to improve resilience.

Author Correction: 863 genomes reveal the origin and domestication of chicken.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

863 genomes reveal the origin and domestication of chicken.

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world's most numerous domestic animal.

Whole metagenome sequencing of cecum microbiomes in Ethiopian indigenous chickens from two different altitudes reveals antibiotic resistance genes.

We analyzed the whole genomes of cecum microbiomes of Ethiopian indigenous chickens from two distinct geographical zones: Afar (AF) district (Dulecha, 730 m above sea level) and Amhara (AM) district (Menz Gera Midir, 3300 m). Through metagenomic analysis we found that microbial populations were mainly dominated by Bacteroidetes and Firmicutes. We identified 2210 common genes in the two groups. LEfSe showed that the distribution of Coprobacter, Geobacter, Cronobacter, Alloprevotella, and Dysgonomonas were more abundant in AF than AM. Analyses using KEGG, eggNOG, and CAZy databases indicated that the pathways of metabolism, genetic information processing, environmental information processing, and cellular process were significantly enriched. Functional abundance was found to be associated with the nutrient absorption and microbial localization of indigenous chickens. We also investigated antibiotic resistant genes and found antibiotics like LSM, cephalosporin, and tetracycline were significantly more abundant in AF than AM.

Characterization of microRNA and mRNA expression profiles in skin tissue between early-feathering and late-feathering chickens.

BACKGROUND: Early feathering and late feathering in chickens are sex-linked phenotypes, which have commercial application in the poultry industry for sexing chicks at hatch and have important impacts on performance traits. However, the genetic mechanism controlling feather development and feathering patterns is unclear. Here, miRNA and mRNA expression profiles in chicken wing skin tissues were analysed through high-throughput transcriptomic sequencing, aiming to understand the biological process of follicle development and the formation of different feathering phenotypes. RESULTS: Compared to the N1 group with no primary feathers extending out, 2893 genes and 31 miRNAs displayed significantly different expression in the F1 group with primary feathers longer than primary-covert feathers, and 1802 genes and 11 miRNAs in the L2 group displayed primary feathers shorter than primary-covert feathers. Only 201 altered genes and 3 altered miRNAs were identified between the N1 and L2 groups (fold change > 2, q value < 0.01). Both sequencing and qPCR tests revealed that PRLR was significantly decreased in the F1 and L2 groups compared to the N1 group, whereas SPEF2 was significantly decreased in the F1 group compared to the N1 or L2 group. Functional analysis revealed that the altered genes or targets of altered miRNAs were involved in multiple biological processes and pathways related to feather growth and development, such as the Wnt signalling pathway, the TGF-beta signalling pathway, the MAPK signalling pathway, epithelial cell differentiation, and limb development. Integrated analysis of miRNA and mRNA showed that 14 pairs of miRNA-mRNA negatively interacted in the process of feather formation. CONCLUSIONS: Transcriptomic sequencing of wing skin tissues revealed large changes in F1 vs. N1 and L2 vs. N1, but few changes in F1 vs. L2 for both miRNA and mRNA expression. PRLR might only contribute to follicle development, while SPEF2 was highly related to the growth rate of primary feathers or primary-covert feathers and could be responsible for early and late feather formation. Interactions between miR-1574-5p/NR2F, miR-365-5p/JAK3 and miR-365-5p/CDK6 played important roles in hair or feather formation. In all, our results provide novel evidence to understand the molecular regulation of follicle development and feathering phenotype.

Circular RNAs are abundant and dynamically expressed during embryonic muscle development in chickens.

The growth and development of skeletal muscle is regulated by proteins as well as non-coding RNAs. Circular RNAs (circRNAs) are universally expressed in various tissues and cell types, and regulate gene expression in eukaryotes. To identify the circRNAs during chicken embryonic skeletal muscle development, leg muscles of female Xinghua (XH) chicken at three developmental time points 11 embryo age (E11), 16 embryo age (E16) and 1 day post hatch (P1) were performed RNA sequencing. We identified 13,377 circRNAs with 3,036 abundantly expressed and most were derived from coding exons. A total of 462 differentially expressed circRNAs were identified (fold change > 2; q-value < 0.05). Parental genes of differentially expressed circRNAs were related to muscle biological processes. There were 946 exonic circRNAs have been found that harbored one or more miRNA-binding site for 150 known miRNAs. We validated that circRBFOX2s promoted cell proliferation through interacted with miR-206. These data collectively indicate that circRNAs are abundant and dynamically expressed during embryonic muscle development and could play key roles through sequestering miRNAs as well as other functions.

miR-16 controls myoblast proliferation and apoptosis through directly suppressing Bcl2 and FOXO1 activities.

Myogenesis mainly involves several steps including myoblast proliferation, differentiation, apoptosis and fusion. Except for muscle specific regulators, few miRNAs were proved to coordinate this complex process. Here, we reported that miR-16 inhibited myoblast proliferation and promoted myoblast apoptosis by directly targeting Bcl2 and FOXO1. The expression level of miR-16 was significantly decreased in the hypertrophic pectoral muscle compared to the normal pectoral muscle in chicken. In vitro, elevating miR-16 significantly inhibited myoblast proliferation and promoted myoblast apoptosis, resulting in about 11.2% cells arrested in G1 phase and 12.3% apoptotic cells in the early stage. Bioinformatic and biochemical analyses revealed Bcl2 and FOXO1 as direct targets of miR-16. Consist to the effect of miR-16 on myogenesis, specific inhibition of Bcl2 or FOXO1 significantly suppressed myoblast proliferation and induced myoblast apoptosis, indicating that both Bcl2 and FOXO1 contributed to miR-16 regulatory function in myogenesis. Interestingly, FOXO1, as the core target, mediated multiple growth-related pathways induced by miR-16 such as PI3K-AKT-MAPK and PI3K-AKT-mTOR. Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) revealed that 234 annotated genes bound by FOXO1 in the early-differentiated myoblasts, which were significantly enriched in myogenic proliferation, death and hypotrophy. Altogether, we proposed that miR-16 acted as a coordinated mediator to suppress myogenesis in avian through the control of myoblast proliferation and apoptosis. These findings have provided a novel mechanism whereby miR-16 represses Bcl2 and FOXO1 expression to maintain myoblast growth and skeletal muscle mass.

Novel MicroRNA Involved in Host Response to Avian Pathogenic Escherichia coli Identified by Deep Sequencing and Integration Analysis.

Avian pathogenic Escherichia coli (APEC) causes one of the most common bacterial diseases of poultry worldwide. Effective control methods are therefore desirable and will be facilitated by a better understanding of the host response to the pathogen. Currently, microRNAs (miRNAs) involved in host resistance to APEC are unknown. Here, we applied RNA sequencing to explore the changed miRNAs and deregulated genes in the spleen of three groups of broilers: nonchallenged (NC), APEC-challenged with mild pathology (CM), and APEC-challenged with severe pathology (CS). Twenty-seven differentially expressed miRNAs (fold change >1.5; P value <0.01) were identified, including 13 miRNAs between the NC and CM, 17 between the NC and CS, and 14 between the CM and CS groups. Through functional analysis of these miRNA targets, 12 immune-related biological processes were found to be significantly enriched. Based on combined analyses of differentially expressed miRNAs and mRNAs within each of the three groups, 43 miRNA-mRNA pairs displayed significantly negative correlations (r < -0.8). Notably, gga-miR-429 was greatly increased in the CS group compared to levels in both the CM and NC groups. In vitro, gga-miR-429 directly repressed luciferase reporter gene activity via binding to 3' untranslated regions of TMEFF2, NTRK2, and SHISA2. Overexpression of gga-miR-429 in the HD11 macrophage cell line significantly inhibited TMEFF2 and SHISA2 expression, which are involved in the lipopolysaccharide-induced platelet-derived growth factor (PDGF) and Wnt signaling pathways. In summary, we provide the first report characterizing the miRNA changes during APEC infection, which may help to shed light on the roles of these recently identified genetic elements in the mechanisms of host resistance and susceptibility to APEC.

A short insertion mutation disrupts genesis of miR-16 and causes increased body weight in domesticated chicken.

Body weight is one of the most important quantitative traits with high heritability in chicken. We previously mapped a quantitative trait locus (QTL) for body weight by genome-wide association study (GWAS) in an F2 chicken resource population. To identify the causal mutations linked to this QTL, expression profiles were determined on livers of high-weight and low-weight chicken lines by microarray. Combining the expression pattern with SNP effects by GWAS, miR-16 was identified as the most likely potential candidate with a 3.8-fold decrease in high-weight lines. Re-sequencing revealed that a 54-bp insertion mutation in the upstream region of miR-15a-16 displayed high allele frequencies in high-weight commercial broiler line. This mutation resulted in lower miR-16 expression by introducing three novel splicing sites instead of the missing 5' terminal splicing of mature miR-16. Elevating miR-16 significantly inhibited DF-1 chicken embryo cell proliferation, consistent with a role in suppression of cellular growth. The 54-bp insertion was significantly associated with increased body weight, bone size and muscle mass. Also, the insertion mutation tended towards fixation in commercial broilers (Fst > 0.4). Our findings revealed a novel causative mutation for body weight regulation that aids our basic understanding of growth regulation in birds.

Characterization of miR-206 Promoter and Its Association with Birthweight in Chicken.

miRNAs have been widely investigated in terms of cell proliferation and differentiation. However, little is known about their effects on bird growth. Here we characterized the promoter of miR-206 in chicken and found that the preferable promoter was located in 1200 bp upstream of pri-miR-206. In this region, many key transcription factors, including MyoD, c-Myb, CEBPα/ß, AP-4, RAP1, Brn2, GATA-1/2/3, E47, Sn, upstream stimulatory factor (USF) and CdxA, were predicted to bind and interact with miR-206 promoter. Overexpression of MyoD sharply increased miR-206 expression in both fibroblast and myoblast cells, and also the regulation in the myoblast cells was much stronger, indicating that miR-206 was regulated by MyoD combined with other muscle specific transcriptional factors. Aiming to further investigate the relationship between miR-206 mutation and transcriptional expression, total of 23 SNPs were identified in the two distinct bird lines by sequencing. Interestingly, the motif bound by MyoD was individually destroyed by G-to-C mutation located at 419 bp upstream of miR-206 precursor. Co-transfecting MyoD and miR-206 promoter in DF-1 cells, the luciferase activity of promoter containing homozygous GG types was significantly higher than CC ones (p < 0.05). Thus, this mutation caused low expression of miR-206. Consistently, eight variants including G-419C mutation exhibited a great effect on birthweight through maker-trait association analysis in F2 population (p < 0.05). Additionally, the regulation of miR-206 on embryo muscle mass mainly by increasing MyoG and muscle creatine kinase (MCK) expression (p < 0.05) with little change in MyoD, TMEM8C and myosin heavy chain (MHC). In conclusion, our findings provide a novel mutation destroying the promoter activity of miR-206 in birds and shed new light to understand the regulation mechanism of miR-206 on the embryonic muscle growth.

Deep Sequencing Analysis of miRNA Expression in Breast Muscle of Fast-Growing and Slow-Growing Broilers.

Growth performance is an important economic trait in chicken. MicroRNAs (miRNAs) have been shown to play important roles in various biological processes, but their functions in chicken growth are not yet clear. To investigate the function of miRNAs in chicken growth, breast muscle tissues of the two-tail samples (highest and lowest body weight) from Recessive White Rock (WRR) and Xinghua Chickens (XH) were performed on high throughput small RNA deep sequencing. In this study, a total of 921 miRNAs were identified, including 733 known mature miRNAs and 188 novel miRNAs. There were 200, 279, 257 and 297 differentially expressed miRNAs in the comparisons of WRRh vs. WRRl, WRRh vs. XHh, WRRl vs. XHl, and XHh vs. XHl group, respectively. A total of 22 highly differentially expressed miRNAs (fold change > 2 or < 0.5; p-value < 0.05; q-value < 0.01), which also have abundant expression (read counts > 1000) were found in our comparisons. As far as two analyses (WRRh vs. WRRl, and XHh vs. XHl) are concerned, we found 80 common differentially expressed miRNAs, while 110 miRNAs were found in WRRh vs. XHh and WRRl vs. XHl. Furthermore, 26 common miRNAs were identified among all four comparisons. Four differentially expressed miRNAs (miR-223, miR-16, miR-205a and miR-222b-5p) were validated by quantitative real-time RT-PCR (qRT-PCR). Regulatory networks of interactions among miRNAs and their targets were constructed using integrative miRNA target-prediction and network-analysis. Growth hormone receptor (GHR) was confirmed as a target of miR-146b-3p by dual-luciferase assay and qPCR, indicating that miR-34c, miR-223, miR-146b-3p, miR-21 and miR-205a are key growth-related target genes in the network. These miRNAs are proposed as candidate miRNAs for future studies concerning miRNA-target function on regulation of chicken growth.

Identification and characterization of genes that control fat deposition in chickens.

BACKGROUND: Fat deposits in chickens contribute significantly to meat quality attributes such as juiciness, flavor, taste and other organoleptic properties. The quantity of fat deposited increases faster and earlier in the fast-growing chickens than in slow-growing chickens. In this study, Affymetrix Genechip® Chicken Genome Arrays 32773 transcripts were used to compare gene expression profiles in liver and hypothalamus tissues of fast-growing and slow-growing chicken at 8 wk of age. Real-time RT-PCR was used to validate the differential expression of genes selected from the microarray analysis. The mRNA expression of the genes was further examined in fat tissues. The association of single nucleotide polymorphisms of four lipid-related genes with fat traits was examined in a F2 resource population. RESULTS: Four hundred genes in the liver tissues and 220 genes hypothalamus tissues, respectively, were identified to be differentially expressed in fast-growing chickens and slow-growing chickens. Expression levels of genes for lipid metabolism (SULT1B1, ACSBG2, PNPLA3, LPL, AOAH) carbohydrate metabolism (MGAT4B, XYLB, GBE1, PGM1, HKDC1)cholesttrol biosynthesis (FDPS, LSS, HMGCR, NSDHL, DHCR24, IDI1, ME1) HSD17B7 and other reaction or processes (CYP1A4, CYP1A1, AKR1B1, CYP4V2, DDO) were higher in the fast-growing White Recessive Rock chickens than in the slow-growing Xinghua chickens. On the other hand, expression levels of genes associated with multicellular organism development, immune response, DNA integration, melanin biosynthetic process, muscle organ development and oxidation-reduction (FRZB, DMD, FUT8, CYP2C45, DHRSX, and CYP2C18) and with glycol-metabolism (GCNT2, ELOVL 6, and FASN), were higher in the XH chickens than in the fast-growing chickens. RT-PCR validated high expression levels of nine out of 12 genes in fat tissues. The G1257069A and T1247123C of the ACSBG2 gene were significantly associated with abdominal fat weight. The G4928024A of the FASN gene were significantly associated with fat bandwidth, and abdominal fat percentage. The C4930169T of the FASN gene was associated with abdominal fat weight while the A59539099G of the ELOVL 6 was significantly associated with subcutaneous fat. The A8378815G of the DDT was associated with fat band width. CONCLUSION: The differences in fat deposition were reflected with differential gene expressions in fast and slow growing chickens.

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

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

MicroRNA profile analysis on duck feather follicle and skin with high-throughput sequencing technology.

Skin acts an important protection role in animal survival and it evolves with the animal divergence. We identified the conserved miRNA families of skin among duck and other species. Cluster analysis showed that the species with similar skin characteristics were clustered into the same group, indicating miRNAs are important in skin function and skin evolution. The miRNA profiles demonstrated that different miRNA regulation mechanism may exist in contour feather follicles (with the surrounding skin) and down feather follicles (with the surrounding skin). Comparing the highly abundant miRNAs with those of mammalian hair follicles and skins, different miRNAs and miRNA families were found, suggesting the different ways in feather follicles and mammalian hair follicles. Bioinformatics prediction indicated that seven miRNAs probably targeted the genes of Wnt/ß-catenin, Shh/BMP and Notch pathways which were important in feather morphogenesis. Further analysis should be conducted to experimentally validate the relationships between miRNAs and their predicted target genes because the target genes were based exclusively upon the bioinformatics.

Let-7b regulates the expression of the growth hormone receptor gene in deletion-type dwarf chickens.

BACKGROUND: A deletion mutation in the growth hormone receptor (GHR) gene results in the inhibition of skeletal muscle growth and fat deposition in dwarf chickens. We used microarray techniques to determine microRNA (miRNA) and mRNA expression profiles of GHR in the skeletal muscles of 14-day-old embryos as well as 7-week-old deletion-type dwarf and normal-type chickens. Our aim was to elucidate the miRNA regulation of GHR expression with respect to growth inhibition and fat deposition. RESULTS: At the same developmental stages, different expression profiles in skeletal muscles of dwarf and normal chickens occurred for four miRNAs (miR-1623, miR-181b, let-7b, and miR-128). At different developmental stages, there was a significant difference in the expression profiles of a greater number of miRNAs. Eleven miRNAs were up-regulated and 18 down-regulated in the 7-week-old dwarf chickens when compared with profiles in 14-day-old embryos. In 7-week-old normal chickens, seven miRNAs were up-regulated and nine down-regulated compared with those in 14-day-old embryos. In skeletal muscles, 22 genes were up-regulated and 33 down-regulated in 14-day-old embryos compared with 7-week-old dwarf chickens. Sixty-five mRNAs were up-regulated and 108 down-regulated in 14-day-old embryos as compared with 7-week-old normal chickens. Thirty-four differentially expressed miRNAs were grouped into 18 categories based on overlapping seed and target sequences. Only let-7b was found to be complementary to its target in the 3' untranslated region of GHR, and was able to inhibit its expression. Kyoto Encyclopedia of Genes and Genomes pathway analysis and quantitative polymerase chain reactions indicated there were three main signaling pathways regulating skeletal muscle growth and fat deposition of chickens. These were influenced by let-7b-regulated GHR. Suppression of the cytokine signaling 3 (SOCS3) gene was found to be involved in the signaling pathway of adipocytokines. CONCLUSIONS: There is a critical miRNA, let-7b, involved in the regulation of GHR. SOCS3 plays a critical role in regulating skeletal muscle growth and fat deposition via let-7b-mediated GHR expression.

Analysis of muscle and ovary transcriptome of Sus scrofa: assembly, annotation and marker discovery.

Pig (Sus scrofa) is an important organism for both agricultural and medical purpose. This study aims to investigate the S. scrofa transcriptome by the use of Roche 454 pyrosequencing. We obtained a total of 558 743 and 528 260 reads for the back-leg muscle and ovary tissue each. The overall 1 087 003 reads give rise to 421 767 341 bp total residues averaging 388 bp per read. The de novo assemblies yielded 11 057 contigs and 60 270 singletons for the back-leg muscle, 12 204 contigs and 70 192 singletons for the ovary and 18 938 contigs and 102 361 singletons for combined tissues. The overall GC content of S. scrofa transcriptome is 42.3% for assembled contigs. Alternative splicing was found within 4394 contigs, giving rise to 1267 isogroups or genes. A total of 56 589 transcripts are involved in molecular function (40 916), biological process (38 563), cellular component (35 787) by further gene ontology analyses. Comparison analyses showed that 336 and 553 genes had significant higher expression in the back-leg muscle and ovary each. In addition, we obtained a total of 24 214 single-nucleotide polymorphisms and 11 928 simple sequence repeats. These results contribute to the understanding of the genetic makeup of S. scrofa transcriptome and provide useful information for functional genomic research in future.