Identification of candidate genes and mutations in QTL regions for chicken growth using bioinformatic analysis of NGS and SNP-chip data

A primer for Fibroblast Growth Factor 16 (FGF16)

During the discovery of the Fibroblast Growth Factor superfamily, scientists were determined to uncover all the genes that encoded FGF proteins. In 1998, FGF16 was discovered with classical cloning techniques in human and rat heart samples. FGF16 loss- and gain-of-function experiments in several organisms demonstrated a conserved function in vertebrates, and as a component of the FGF9 subfamily of ligands (FGF-E/-9/-20), is functionally conserved and sufficient to rescue loss-of-function phenotypes in invertebrates, like C. elegans. FGF16 has a broad expression pattern, predominantly expressed in brown adipose tissue, heart, with low but detectable levels in the brain, olfactory bulb, inner ear, muscle, thymus, pancreas, spleen, stomach, small intestine, and gonads (testis and ovary). FGF16 is also expressed moderately in the late developing limb bud. Despite its expression levels, this ligand plays notable roles in autopod metacarpal development; loss of one allele causes congenital metacarpal 4-5 fusion and hand deformities in humans. The broad expression pattern of FGF16 in several tissues underscores its multifaceted roles in stem cell maintenance, proliferation, cell fate specification, and metabolism.

Crim1 inhibits angiotensin II-induced hypertrophy and preserves Kv4.2 expression in cardiomyocytes

Objectives

Angiotensin II (Ang II) plays a key role in the regulation of myocardial hypertrophy via downstream cysteine-rich transmembrane bone morphogenetic protein regulator 1 (Crim1). However, it is still unclear whether Crim1 is involved in ionic channel remodeling. The study aimed to explore the effects of Crim1 on transient outward potassium current (Ito) and Kv4.2 (the main subunit of I_to channel) expression in hypertrophic ventricular cardiomyocytes.

Materials and Methods

The ventricular cardiomyocytes were isolated from the neonatal rats. Hypertrophy was induced by Ang II. Crim1 expression was modulated by using adenovirus transfection. The expression of myosin heavy chain beta (β-MHC), Crim1, and Kv4.2 was determined by RT-qPCR and western blot. The cellular surface area was assessed using Image J software. I_to was recorded by the whole-cell patch clamp technique.

Results

Ang II-induced hypertrophy in cardiomyocytes was identified by their larger cellular surface area and higher mRNA expression of β-MHC. Ang II significantly decreased the expression of Crim1 and Kv4.2 and reduced I_to current density. However, Crim1 overexpression abolished the Ang II-induced hypertrophy and preserved the expression of Kv4.2 and I_to current density.

Conclusion

Crim1 overexpression inhibits Ang II-induced hypertrophy and preserves I_to current density via up-regulating Kv4.2 in ventricular cardiomyocytes from neonatal rats. Crim1 could have a role in the development of ventricular arrhythmia in hypertrophic hearts.

Super-assembly of integrated gold magnetic assay with loop-mediated isothermal amplification for point-of-care testing

With the increasing global threat of various diseases and infections, it is essential to develop a fast, low-cost, and easy-to-use point-of-care testing (POCT) system for inspections at all levels of medical institutions and self-examination at home. In this work, gold magnetic nanoparticles (GMNPs) are used as the key material, and a rapid visual detection method is designed through integrating loop-mediated isothermal amplification (LAMP) and lateral flow assay (LFA) biosensor for detecting a variety of analytes which includes whole blood, buccal swabs, and DNA. It is worth to note that the proposed method does not need DNA extraction. Furthermore, uracil DNA glycosylase (UDG) is employed to eliminate carrier contamination for preventing false positive results. The whole detection process can be finished within 25 min. The accuracy of detection is measured by assessing the polymorphisms of the methylenetetrahydrofolate reductase (MTHFR) C677T. The detection limit of the newly developed extraction-free detection system for MTHFR C677T is 0.16 ng/μL. A preliminary clinical study of the proposed method is carried out by analyzing 600 clinical samples (including 200 whole blood samples, 100 buccal swabs, and 300 genomic DNA samples). The results indicate that the proposed method is 100% consistent with the sequencing results which provides a new choice for POCT and shows a broad application prospect in all levels of medical clinics and at home.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material (details for MTHFR C677T primer sequences, the cell count results of samples at different dilution ratios, genotyping results and frequency samples, a Hardy-Weinberg equilibrium test, the sensitivity of the system, detection results of multiple samples, and optimization of the system) is available in the online version of this article at 10.1007/s12274-022-4692-9.

Transcriptomic profile of leg muscle during early growth and development in Haiyang yellow chicken

Skeletal muscle growth and development from embryo to adult consists of a series of carefully regulated changes in gene expression. This study aimed to identify candidate genes involved in chicken growth and development and to investigate the potential regulatory mechanisms of early growth in Haiyang yellow chicken. RNA sequencing was used to compare the transcriptomes of chicken muscle tissues at four developmental stages. In total, 6150 differentially expressed genes (DEGs) (|fold change| ≥ 2; false discovery rate (FDR) ≤ 0.05) were detected by pairwise comparison in female chickens. Functional analysis showed that the DEGs were mainly involved in the processes of muscle growth and development and cell differentiation. Many of the DEGs, such as MSTN, MYOD1, MYF6, MYF5, and IGF1, were related to chicken growth and development. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the DEGs were significantly enriched in four pathways related to growth and development: extracellular matrix (ECM)–receptor interaction, focal adhesion, tight junction, and insulin signalling pathways. A total of 42 DEGs assigned to these pathways are potential candidate genes for inducing the differences in growth among the four development stages, such as MYH1A, EGF, MYLK2, MYLK4, and LAMB3. This study identified a range of genes and several pathways that may be involved in regulating early growth.

A Missense Mutation in the MYBPH Gene Is Associated With Abdominal Fat Traits in Meat-Type Chickens

Chicken is an important source of protein for human nutrition and a model system for growth and developmental biology. Although the genetic architecture of quantitative traits in meat-type chickens has been the subject of ongoing investigation, the identification of mutations associated with carcass traits of economic interest remains challenging. Therefore, our aim was to identify predicted deleterious mutation, which potentially affects protein function, and test if they were associated with carcass traits in chickens. For that, we performed a genome-wide association analysis (GWAS) for breast, thigh and drumstick traits in meat-type chickens and detected 19 unique quantitative trait loci (QTL). We then used: (1) the identified windows; (2) QTL for abdominal fat detected in a previous study with the same population and (3) previously obtained whole genome sequence data, to identify 18 predicted deleterious single nucleotide polymorphisms (SNPs) in those QTL for further association with breast, thigh, drumstick and abdominal fat traits. Using the additive model, a predicted deleterious SNP c.482C > T (SIFT score of 0.4) was associated (p-value < 0.05) with abdominal fat weight and percentage. This SNP is in the second exon of the MYBPH gene, and its allele frequency deviates from Hardy–Weinberg equilibrium. In conclusion, our study provides evidence that the c.482C > T SNP in the MYBPH gene is a putative causal mutation for fat deposition in meat-type chickens.

Identification of new genes and quantitative trait locis associated with growth curve parameters in F2 chicken population using genome-wide association study

The markers which are correlated with the growth curve parameters help in understanding the characteristics of individual growth during the rearing of livestock. This study aimed to identify a set of biomarkers through a GWAS for growth curve parameters in crossbred chickens using the Illumnia 60K chicken SNP Beadchip. Growth data were collected from a total of 301 birds from cross of a broiler line and native chickens. Using the Gompertz-Laird model, two growth curve parameters, the instantaneous growth rate per day (L) and the coefficient of relative growth or maturing index (k), were estimated. The L and k were used to estimate five derived parameters, namely asymptotic (mature) body weight, body weight at inflection point, age at the inflection point, average growth rate and maximum growth rate. These parameters were considered as phenotypic values in the GWAS based on generalized linear models. The results of the GWAS indicated 21 significant markers, which were located near or within 46 genes. A number of these genes, such as GH, RET, GRB14, FTSJ3 and CCK, are important for growth and meat quality in chickens, and some of them are growth related in other species such as sheep and cattle (GPI, XIRP2, GALNTL6, BMS1, THSD4, TRHDE, SHISA9, ACSL6 and DYNC1LI2). The other genes are associated with developmental biological pathways. These genes are particuarly related to body weight, average daily gain and growth QTL. The results of this study can shed light on the genetic mechanism of biological functions of growth factors in broiler chickens, which is useful for developing management practices and accelerating genetic progress in breeding programs.

Regulating appetite in broilers for improving body and muscle development - A review

Appetite is the desire for feed and water and the voluntary intake of feed and is an important regulator of livestock productivity and animal health. Economic traits such as growth rate and muscle development (meat deposition) in broilers are directly correlated to appetite. Factors that may influence appetite include environmental factors, such as stress and temperature variation, and animal‐specific factors, such as learning period, eating capacity and preferences. Feed preferences have been reported to be determined in early life, and this period is important in broilers due to their fast growth and relatively short growth trajectories. This may be of importance when contemplating the use of more circular and sustainable feeds and the optimization of appetite for these feeds. The objective of this review was to review the biological mechanisms underlying appetite using data from human, animal and bird models and to consider the option for modulating appetite particularly as it relates to broiler chickens.

Dissecting Epistatic QTL for Blood Pressure in Rats: Congenic Strains versus Heterogeneous Stocks, a Reality Check

Advances in molecular genetics have provided well-defined physical genetic maps and large numbers of genetic markers for both model organisms and humans. It is now possible to gain a fundamental understanding of the genetic architecture underlying quantitative traits, of which blood pressure (BP) is an important example. This review emphasizes analytical techniques and results obtained using the Dahl salt-sensitive (S) rat as a model of hypertension by presenting results in detail for three specific chromosomal regions harboring genetic elements of increasing complexity controlling BP. These results highlight the critical importance of genetic interactions (epistasis) on BP at all levels of structure, intragenic, intergenic, intrachromosomal, interchromosomal, and across whole genomes. In two of the three examples presented, specific DNA structural variations leading to biochemical, physiological, and pathological mechanisms are well defined. This proves the usefulness of the techniques involving interval mapping followed by substitution mapping using congenic strains. These classic techniques are compared to newer approaches using sophisticated statistical analysis on various segregating or outbred model-organism populations, which in some cases are uniquely useful in demonstrating the existence of higher-order interactions. It is speculated that hypertension as an outlier quantitative phenotype is dependent on higher-order genetic interactions. The obstacle to the identification of genetic elements and the biochemical/physiological mechanisms involved in higher-order interactions is not theoretical or technical but the lack of future resources to finish the job of identifying the individual genetic elements underlying the quantitative trait loci for BP and ascertaining their molecular functions. © 2019 American Physiological Society. Compr Physiol 9:1305-1337, 2019.

Bidirectional Selection for Body Weight on Standing Genetic Variation in a Chicken Model

Experimental populations of model organisms provide valuable opportunities to unravel the genomic impact of selection in a controlled system. The Virginia body weight chicken lines represent a unique resource to investigate signatures of selection in a system where long-term, single-trait, bidirectional selection has been carried out for more than 60 generations. At 55 generations of divergent selection, earlier analyses of pooled genome resequencing data from these lines revealed that 14.2% of the genome showed extreme differentiation between the selected lines, contained within 395 genomic regions. Here, we report more detailed analyses of these data exploring the regions displaying within- and between-line genomic signatures of the bidirectional selection applied in these lines. Despite the strict selection regime for opposite extremes in body weight, this did not result in opposite genomic signatures between the lines. The lines often displayed a duality of the sweep signatures, where an extended region of homozygosity in one line, in contrast to mosaic pattern of heterozygosity in the other line. These haplotype mosaics consisted of short, distinct haploblocks of variable between-line divergence, likely the results of a complex demographic history involving bottlenecks, introgressions and moderate inbreeding. We demonstrate this using the example of complex haplotype mosaicism in the growth1 QTL. These mosaics represent the standing genetic variation available at the onset of selection in the founder population. Selection on standing genetic variation can thus result in different signatures depending on the intensity and direction of selection.

A PCR-lateral flow assay system based on gold magnetic nanoparticles for CYP2C19 genotyping and its clinical applications

Presently, clopidogrel is the standard therapeutic drugs for antiplatelet therapy. Variants in the CYP2C19 gene influence the clinical response of clopidogrel. Thus, the US Food and Drug Administration suggested CYP2C19 genotyping needs to identify before taking medicine. Due to high cost, time consuming, and sophisticated instruments, current single nucleotide polymorphism detection methods are limited in clinical application. In the present study, we established a genotyping method for CYP2C19, which combines amplification refractory mutation system (ARMS)-PCR with a lateral flow assay used gold magnetic nanoparticles (GMNPs) named as PCR-gold magnetic lateral flow assay system (PCR-GoldMag LFA). The PCR products with specific genotype can be explained within 5 minutes, either through visual or by a magnetic reader automatically according to the captured GMNPs probes on the test lines of the LFA strips. The limit of detection of this method is 5 ng of genomic DNA. The PCR-GoldMag LFA system was applied in a clinical trial with 1356 samples of Han Chinese. The concordance rate between the LFA system and sequencing is 99.93%. The allele frequency of CYP2C19*2 and CYP2C19*3 are 30.38 and 7.08% in Han Chinese, respectively. This method provides a new way in the clinical application of CYP2C19 genotyping to guide the clopidogrel medication.

A direct isothermal amplification system adapted for rapid SNP genotyping of multifarious sample types

Genotyping of single nucleotide polymorphisms (SNPs) in point-of-care (POC) settings could be further improved through simplifying the treatment of samples. In this study, we devised an accurate, rapid and easy-to-use SNP detection system based on direct loop-mediated isothermal amplification (LAMP) without DNA extraction, known as Direct-LAMP. Samples from various sources (including whole blood, dried blood spot, buccal swab and saliva), treated with NaOH, can be used directly in amplification. The turnaround time was about 30 min from sample collection to provision of results. The accuracy was evaluated by assessing the polymorphisms of methylenetetrahydrofolate reductase (MTHFR) C677T and aldehyde dehydrogenase-2 (ALDH2) Glu504Lys, which are better known for their critical role in folate and ethanol metabolism, respectively. Completely consistent genotyping results reveal that Direct-LAMP is generally concordant with sequencing. This system can serve as a very promising platform in the fields of disease predisposition, drug metabolism and personalized medicine.

Multiple SNPs Detection Based on Lateral Flow Assay for Phenylketonuria Diagnostic

Single nucleotide polymorphisms (SNPs) are closely related to genetic diseases, but current SNP detection methods, such as DNA microarrays that include tedious procedures and expensive, sophisticated instruments, are unable to perform rapid SNPs detection in clinical practice, especially for those multiple SNPs related to genetic diseases. In this study, we report a sensitive, low cost, and easy-to-use point-of-care testing (POCT) system formed by combining amplification refractory mutation system (ARMS) polymerase chain reaction with gold magnetic nanoparticles (GMNPs) and lateral flow assay (LFA) noted as the ARMS-LFA system, which allow us to use a uniform condition for multiple SNPs detection simultaneously. The genotyping results can be explained by a magnetic reader automatically or through visual interpretation according to the captured GMNPs probes on the test and control lines of the LFA device. The high sensitivity (the detection limit of 0.04 pg/μL with plasmid) and specificity of this testing system were found through genotyping seven pathogenic SNPs in phenylalanine hydroxylase gene ( PAH, the etiological factor of phenylketonuria). This system can also be applied in DNA quantification with a linear range from 0.02 to 2 pg/μL of plasmid. Furthermore, this ARMS-LFA system was applied to clinical trials for screening the seven pathogenic SNPs in PAH of 23 families including 69 individuals. The concordance rate of the genotyping results detected by the ARMS-LFA system was up to 97.8% compared with the DNA sequencing results. This method is a very promising POCT in the detection of multiple SNPs caused by genetic diseases.

Transcriptional analysis of abdominal fat in chickens divergently selected on bodyweight at two ages reveals novel mechanisms controlling adiposity: validating visceral adipose tissue as a dynamic endocrine and metabolic organ

Background

Decades of intensive genetic selection in the domestic chicken (Gallus gallus domesticus) have enabled the remarkable rapid growth of today’s broiler (meat-type) chickens. However, this enhanced growth rate was accompanied by several unfavorable traits (i.e., increased visceral fatness, leg weakness, and disorders of metabolism and reproduction). The present descriptive analysis of the abdominal fat transcriptome aimed to identify functional genes and biological pathways that likely contribute to an extreme difference in visceral fatness of divergently selected broiler chickens.

Methods

We used the Del-Mar 14 K Chicken Integrated Systems microarray to take time-course snapshots of global gene transcription in abdominal fat of juvenile [1-11 weeks of age (wk)] chickens divergently selected on bodyweight at two ages (8 and 36 wk). Further, a RNA sequencing analysis was completed on the same abdominal fat samples taken from high-growth (HG) and low-growth (LG) cockerels at 7 wk, the age with the greatest divergence in body weight (3.2-fold) and visceral fatness (19.6-fold).

Results

Time-course microarray analysis revealed 312 differentially expressed genes (FDR ≤ 0.05) as the main effect of genotype (HG versus LG), 718 genes in the interaction of age and genotype, and 2918 genes as the main effect of age. The RNA sequencing analysis identified 2410 differentially expressed genes in abdominal fat of HG versus LG chickens at 7 wk. The HG chickens are fatter and over-express numerous genes that support higher rates of visceral adipogenesis and lipogenesis. In abdominal fat of LG chickens, we found higher expression of many genes involved in hemostasis, energy catabolism and endocrine signaling, which likely contribute to their leaner phenotype and slower growth. Many transcription factors and their direct target genes identified in HG and LG chickens could be involved in their divergence in adiposity and growth rate.

Conclusions

The present analyses of the visceral fat transcriptome in chickens divergently selected for a large difference in growth rate and abdominal fatness clearly demonstrate that abdominal fat is a very dynamic metabolic and endocrine organ in the chicken. The HG chickens overexpress many transcription factors and their direct target genes, which should enhance in situ lipogenesis and ultimately adiposity. Our observation of enhanced expression of hemostasis and endocrine-signaling genes in diminished abdominal fat of LG cockerels provides insight into genetic mechanisms involved in divergence of abdominal fatness and somatic growth in avian and perhaps mammalian species, including humans.

Electronic supplementary material

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

A Chromosome-Scale Assembly of the Bactrocera cucurbitae Genome Provides Insight to the Genetic Basis of white pupae

Genetic sexing strains (GSS) used in sterile insect technique (SIT) programs are textbook examples of how classical Mendelian genetics can be directly implemented in the management of agricultural insect pests. Although the foundation of traditionally developed GSS are single locus, autosomal recessive traits, their genetic basis are largely unknown. With the advent of modern genomic techniques, the genetic basis of sexing traits in GSS can now be further investigated. This study is the first of its kind to integrate traditional genetic techniques with emerging genomics to characterize a GSS using the tephritid fruit fly pest Bactrocera cucurbitae as a model. These techniques include whole-genome sequencing, the development of a mapping population and linkage map, and quantitative trait analysis. The experiment designed to map the genetic sexing trait in B. cucurbitae, white pupae (wp), also enabled the generation of a chromosome-scale genome assembly by integrating the linkage map with the assembly. Quantitative trait loci analysis revealed SNP loci near position 42 MB on chromosome 3 to be tightly linked to wp Gene annotation and synteny analysis show a near perfect relationship between chromosomes in B. cucurbitae and Muller elements A-E in Drosophila melanogaster This chromosome-scale genome assembly is complete, has high contiguity, was generated using a minimal input DNA, and will be used to further characterize the genetic mechanisms underlying wp Knowledge of the genetic basis of genetic sexing traits can be used to improve SIT in this species and expand it to other economically important Diptera.

Genome-wide association study for performance traits in chickens using genotype by sequencing approach

Performance traits are economically important and are targets for selection in breeding programs, especially in the poultry industry. To identify regions on the chicken genome associated with performance traits, different genomic approaches have been applied in the last years. The aim of this study was the application of CornellGBS approach (134,528 SNPs generated from a PstI restriction enzyme) on Genome-Wide Association Studies (GWAS) in an outbred F₂ chicken population. We have validated 91.7% of these 134,528 SNPs after imputation of missed genotypes. Out of those, 20 SNPs were associated with feed conversion, one was associated with body weight at 35 days of age (P < 7.86E-07) and 93 were suggestively associated with a variety of performance traits (P < 1.57E-05). The majority of these SNPs (86.2%) overlapped with previously mapped QTL for the same performance traits and some of the SNPs also showed novel potential QTL regions. The results obtained in this study suggests future searches for candidate genes and QTL refinements as well as potential use of the SNPs described here in breeding programs.

Imputation-Based Fine-Mapping Suggests That Most QTL in an Outbred Chicken Advanced Intercross Body Weight Line Are Due to Multiple, Linked Loci

The Virginia chicken lines have been divergently selected for juvenile body weight for more than 50 generations. Today, the high- and low-weight lines show a >12-fold difference for the selected trait, 56-d body weight. These lines provide unique opportunities to study the genetic architecture of long-term, single-trait selection. Previously, several quantitative trait loci (QTL) contributing to weight differences between the lines were mapped in an F₂-cross between them, and these were later replicated and fine-mapped in a nine-generation advanced intercross of them. Here, we explore the possibility to further increase the fine-mapping resolution of these QTL via a pedigree-based imputation strategy that aims to better capture the genetic diversity in the divergently selected, but outbred, founder lines. The founders of the intercross were high-density genotyped, and then pedigree-based imputation was used to assign genotypes throughout the pedigree. Imputation increased the marker density 20-fold in the selected QTL, providing 6911 markers for the subsequent analysis. Both single-marker association and multi-marker backward-elimination analyses were used to explore regions associated with 56-d body weight. The approach revealed several statistically and population structure independent associations and increased the mapping resolution. Further, most QTL were also found to contain multiple independent associations to markers that were not fixed in the founder populations, implying a complex underlying architecture due to the combined effects of multiple, linked loci perhaps located on independent haplotypes that still segregate in the selected lines.

Exploiting Natural Variation to Discover Candidate Genes Involved in Photosynthesis-Related Traits

Naturally occurring genetic variation in plants can be very useful to dissect the complex regulation of primary metabolism as well as of physiological traits such as photosynthesis and photorespiration. The physiological and genetic mechanisms underlying natural variation in closely related species or accessions may provide important information that can be used to improve crop yield. In this chapter we describe in detail the use of a population of introgression lines (ILs), with the Solanum pennellii IL population as a study case, as a tool for the identification of genomic regions involved in the control of photosynthetic efficiency.

A novel lateral flow assay based on GoldMag nanoparticles and its clinical applications for genotyping of MTHFR C677T polymorphisms

Current techniques for single nucleotide polymorphism (SNP) detection require tedious experimental procedures and expensive and sophisticated instruments. In this study, a visual genotyping method has been successfully established via combining ARMS-PCR with gold magnetic nanoparticle (GoldMag)-based lateral flow assay (LFA) and applied to the genotyping of methylenetetrahydrofolate reductase (MTHFR) C677T. C677T substitution of the gene MTHFR leads to an increased risk of diseases. The genotyping result is easily achievable by visual observation within 5 minutes after loading of the PCR products onto the LFA device. The system is able to accurately assess a broad detection range of initial starting genomic DNA amounts from 5 ng to 1200 ng per test sample. The limit of detection reaches 5 ng. Furthermore, our PCR-LFA system was applied to clinical trials for screening 1721 individuals for the C677T genotypes. The concordance rate of the genotyping results detected by PCR-LFA was up to 99.6% when compared with the sequencing results. Collectively, our PCR-LFA has been proven to be rapid, accurate, sensitive, and inexpensive. This new method is highly applicable for C677T SNP screening in laboratories and clinical practices. More promisingly, it could also be extended to the detection of SNPs of other genes.

Selection for divergent body size alters rates of embryonic skeletal muscle formation and muscle gene expression patterns

The impact of divergent selection for body size on embryogenesis is poorly understood. The objective of this experiment was to document skeletal muscle development during embryogenesis in two lines of chickens that display divergent growth as adults. Results reveal that after 54 generations of opposing selection from a common founder population, the embryos from the low weight select (LWS) line develop more rapidly during early embryogenesis than those from the high weight select (HWS) line. Muscle formation during the late embryonic period is more rapid and extensive in the HWS embryo than in the LWS contemporary. Isolated muscle progenitors from embryonic day 10 HWS embryos proliferated more rapidly, forming fibers sooner with a larger size than the LWS cells. The limited myogenic capacity of the LWS progenitor cells is not attributed to altered patterns of expression of Pax7, Pax3 or the myogenic regulatory factor genes. Members of the fibroblast growth factor family are potent mitogens and inhibitors of myoblast differentiation. Transcript abundance of FGF2 and FGF4 was measured in cultures of HWS and LWS progenitors as a function of time. The pattern of expression of FGF4 was similar between HWS and LWS with a large increase between days 1 and 3 followed by a reduction at day 5 of culture. Expression of FGF2 in LWS muscle cells did not change while a significant reduction in FGF2 expression was observed by day 5 in the HWS. Our results indicate that divergent selection for postnatal growth has altered embryonic development.

A Genome-Wide mRNA Screen and Functional Analysis Reveal FOXO3 as a Candidate Gene for Chicken Growth

Chicken growth performance provides direct economic benefits to the poultry industry. However, the underlying genetic mechanisms are unclear. The objective of this study was to identify candidate genes associated with chicken growth and investigate their potential mechanisms. We used RNA-Seq to study the breast muscle transcriptome in high and low tails of Recessive White Rock (WRR_h, WRR_l) and Xinghua chickens (XH_h, XH_l). A total of 60, 23, 153 and 359 differentially expressed genes were detected in WRR_h vs. WRR_l, XH_h vs. XH_l, WRR_h vs. XH_h and WRR_l vs. XH_l, respectively. GO, KEGG pathway and gene network analyses showed that CEBPB, FBXO32, FOXO3 and MYOD1 played key roles in growth. The functions of FBXO32 and FOXO3 were validated. FBXO32 was predominantly expressed in leg muscle, heart and breast muscle. After decreased FBXO32 expression, growth-related genes such as PDK4, IGF2R and IGF2BP3 were significantly down-regulated (P < 0.05). FBXO32 was significantly (P < 0.05) associated with carcass and meat quality traits, but not growth traits. FOXO3 was predominantly expressed in breast and leg muscle. In both of these tissues, the FOXO3 mRNA level in XH was significantly higher than that in WRR chickens with normal body weight (P < 0.05). In DF-1 cells, siRNA knockdown of FOXO3 significantly (P < 0.01) inhibited the MYOD expression and significantly up-regulated (P < 0.01 or P < 0.05) the expression of growth-related genes including CEBPB, FBXO32, GH, GHR, IGF1R, IGF2R, IGF2BP1, IGF2BP3, INSR, PDK1 and PDK4. Moreover, 18 SNPs were identified in FOXO3. G66716193A was significantly (P < 0.05) associated with growth traits. The sites C66716002T, C66716195T and A66716179G were significantly (P < 0.05) associated with growth or carcass traits. These results demonstrated that FOXO3 is a candidate gene influencing chicken growth. Our observations provide new clues to understand the molecular basis of chicken 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.

SNP and INDEL detection in a QTL region on chicken chromosome 2 associated with muscle deposition

Genetic improvement is important for the poultry industry, contributing to increased efficiency of meat production and quality. Because breast muscle is the most valuable part of the chicken carcass, knowledge of polymorphisms influencing this trait can help breeding programs. Therefore, the complete genome of 18 chickens from two different experimental lines (broiler and layer) from EMBRAPA was sequenced, and SNPs and INDELs were detected in a QTL region for breast muscle deposition on chicken chromosome 2 between microsatellite markers MCW0185 and MCW0264 (105,849-112,649 kb). Initially, 94,674 unique SNPs and 10,448 unique INDELs were identified in the target region. After quality filtration, 77% of the SNPs (85,765) and 60% of the INDELs (7828) were retained. The studied region contains 66 genes, and functional annotation of the filtered variants identified 517 SNPs and three INDELs in exonic regions. Of these, 357 SNPs were classified as synonymous, 153 as non-synonymous, three as stopgain, four INDELs as frameshift and three INDELs as non-frameshift. These exonic mutations were identified in 37 of the 66 genes from the target region, three of which are related to muscle development (DTNA, RB1CC1 and MOS). Fifteen non-tolerated SNPs were detected in several genes (MEP1B, PRKDC, NSMAF, TRAPPC8, SDR16C5, CHD7, ST18 and RB1CC1). These loss-of-function and exonic variants present in genes related to muscle development can be considered candidate variants for further studies in chickens. Further association studies should be performed with these candidate mutations as should validation in commercial populations to allow a better explanation of QTL effects.

Variant discovery in a QTL region on chromosome 3 associated with fatness in chickens

Abdominal fat content is an economically important trait in commercially bred chickens. Although many quantitative trait loci (QTL) related to fat deposition have been detected, the resolution for these regions is low and functional variants are still unknown. The current study was conducted aiming at increasing resolution for a region previously shown to have a QTL associated with fat deposition, to detect novel variants from this region and to annotate those variants to delineate potentially functional ones as candidates for future studies. To achieve this, 18 chickens from a parental generation used in a reciprocal cross between broiler and layer lines were sequenced using the Illumina next-generation platform with an initial coverage of 18X/chicken. The discovery of genetic variants was performed in a QTL region located on chromosome 3 between microsatellite markers LEI0161 and ADL0371 (33,595,706-42,632,651 bp). A total of 136,054 unique SNPs and 15,496 unique INDELs were detected in this region, and after quality filtering, 123,985 SNPs and 11,298 INDELs were retained. Of these variants, 386 SNPs and 15 INDELs were located in coding regions of genes related to important metabolic pathways. Loss-of-function variants were identified in several genes, and six of those, namely LOC771163, EGLN1, GNPAT, FAM120B, THBS2 and GGPS1, were related to fat deposition. Therefore, these loss-of-function variants are candidate mutations for conducting further studies on this important trait in chickens.