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

Article information

Asian-Australas J Anim Sci.. 2014;27(3):316-323
1Akita Prefectural Livestock Experiment Station, Daisen 019-1701, Japan.
2National Livestock Breeding Center, Nishigo, Fukushima 961-8511, Japan.
3Graduate School of Agricultural Science, Tohoku University, Sendai 981-8555, Japan.
*Corresponding Author: Takahashi Hideaki. Tel: +81-29-838-8623, Fax: +81-29-838-8606, E-mail: naoe@affrc.go.jp

National Institute of Livestock and Grassland Science, Tsukuba, Ibaraki 305-0901, Japan

Received 2013 August 13; Revised 2013 November 16; Accepted 2013 October 10.

Abstract

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 F2 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 F2 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 F2 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 F2 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.

INTRODUCTION

The Hinai-dori is a slow-growing breed of chicken native to Akita Prefecture, in northern Honshu Island, Japan. Although Hinai-dori meat has a characteristic taste and the breed has been used for a long time, it has decreased in numbers in recent times owing to the introduction of exotic breeds, and for a while was at risk of extinction. The Hinai-dori breed has been conserved by hobbyists who belong to the Preservation Society (PS) of the Hinai-dori Breed and is now mainly used for exhibition purposes. For effective use of the breed, selection experiments have been performed at the Akita Prefectural Livestock Experiment Station (LES) since 1973 when fertilized eggs were introduced to LES from PS with a view towards increasing growth performance (Rikimaru et al., 2013). At present, the sex-averaged body weight of LES line at 14 weeks of age is approximately 0.8 kg heavier than that of PS line (Rikimaru et al., 2013). F1 chickens produced by crossing the improved LES Hinai-dori sires with Rhode Island Red dams have been commercialized as the Hinai-jidori chicken, which is one of the most popular high-quality chickens on the Japanese market (Rikimaru et al., 2007).

The identification and utilization of potential candidate genes for quantitative trait loci (QTL) with significant effects on economically important traits is becoming increasingly important in animal breeding programs. There are currently two main strategies for detecting QTLs: genome scans based on linkage mapping in a cross population and association tests using candidate genes (Andersson, 2001). In the past decade, QTL mapping for chicken growth traits such as body weight has been widely studied (Abasht et al., 2006; Hu et al., 2007). Moreover, there have been numerous association studies investigating the relationship between growth traits and candidate genes, including pre-B-cell colony enhancing factor I (PBEF1) (Han et al., 2011), high mobility group AT-hook 2 (HMGA2) (Song et al., 2011), insulin-like growth factor I (IGF1) (Zhou et al., 2005), pituitary-specific positive transcription factor 1 (PIT1) (Nie et al., 2008), and thyroid hormone responsive spot 14α (THRSP) (Cao et al., 2007) on chromosome 1; insulin-like growth factor binding protein (IGFBP) 1 and 3 (Qu et al., 2009) on chromosome 2; ornithine decarboxylase (ODC) (Uemoto et al., 2011) on chromosome 3; cholecystokinin type A receptor (CCKAR) (Rikimaru et al., 2012) on chromosome 4; calpain 3 (Zhang et al., 2009) on chromosome 5; IGFBP 2 (Lei et al., 2005; Li et al., 2006; Leng et al., 2009) on chromosome 7; growth hormone secretagogue receptor (GHSR) (Fang et al., 2010) on chromosome 9; and insulin-like growth factor 1 receptor (IGF1R) (Lei et al., 2008) on chromosome 10.

The LES line was developed from PS chickens: hence, these individuals possessed genes that influence growth traits. To identify QTLs associated with growth traits in the Hinai-dori breed, an F2 resource population produced by crossing low- and high-growth lines of the breed was analyzed. A significant QTL for body weight at 10 and 14 weeks of age and average daily gain between 4 and 10 weeks and between 10 and 14 weeks of age were accordingly mapped in a common region between ADL0198 (chr 1: 171.7 Mb) and ABR0287 (chr 1: 173.4 Mb) on chromosome 1 (Rikimaru et al., 2011). Thus, the above-mentioned genes, e.g., IGFBP 1–3, ODC, CCKAR, GHSR, and IGF1R, are not objects of our interest, since they are not located on chromosome 1. Moreover, PBEF1 (chr 1: 15.2 Mb), HMGA2 (chr 1: 36.1 Mb), IGF1 (chr 1: 57.3 Mb), PIT1 (chr 1: 96.2 Mb), THRSP (chr 1: 197.8 Mb) are not located between ADL0198 and ABR0287.

Of all the candidate genes for growth traits on chromosome 1, we herein focused our investigation on the motilin receptor (MLNR, chr 1: 173.4 Mb) because it was located just under the QTL peak adjacent to ABR0287 on chromosome 1, and motilin is a key regulatory hormones for appetite (Garthwaite, 1985; Rosenfeld and Garthwaite, 1987), and thus MLNR is a positional and functional candidate gene. Our main objective in this study was to analyze polymorphism of the MLNR gene and test its association with growth traits in Hinai-dori F2 crossbreds.

MATERIALS AND METHODS

Resource population

PS and LES chickens of the Hinai-dori are shown in Figure 1. F1 chickens were produced by crossing three cocks from a low-growth line (obtained from PS) with nine hens from a high-growth line (obtained from LES). In this cross, one to three hens were randomly selected to mate with each male. The F2 individuals were produced by full-sib mating of 17 F1 cocks and 60 F1 hens. We accordingly obtained 418 F2 individuals comprising 206 cocks and 212 hens. The F2 chickens were hatched on the same day, raised in the same chicken house, and fed the same diet ad libitum for the duration of the experiment. All animals received humane care as outlined in the Guidelines for Proper Conduct of Animal Experiments, Science Council of Japan (Science Council of Japan, 2006).

Figure 1

Hinai-dori. (A) Hinai-dori chickens of the Preservation Society (PS) line. Left, male; right, female. (B) Hinai-dori chickens of the Akita Prefectural Livestock Experiment Station (LES) line. Left, female; right, male.

Traits

Body weight was measured at day 0 (BW 0 d), 4 weeks (BW 4 wk), 10 weeks (BW 10 wk), and 14 weeks (BW 14 wk) of age. Average daily gain between 0 and 4 weeks of age (ADG 0–4 wk), between 4 and 10 weeks of age (ADG 4–10 wk), between 10 and 14 weeks of age (ADG 10–14 wk), and 0 and 14 weeks of age (ADG 0–14 wk) was calculated from BW at each week of age.

Identification of MLNR haplotypes

A draft sequence of the chicken genome (May 2006 assembly), available on the University of California, Santa Cruz (UCSC) Genome Browser (2004) and the Ensembl Genome Browser (2004), was used in this study. The nucleotide sequences of the two exons of MLNR in the parent individuals were determined by polymerase chain reaction (PCR) amplification followed by direct sequencing to determine nucleotide variance of the gene in the resource family. Each genomic DNA of the resource family was purified using a SepaGene kit (EIDIA, Tokyo, Japan). Four PCR primers were designed to separately amplify the two exons of MLNR and the PCRs were performed (Table 1). For the PCRs, we used 15-μL reaction volumes containing the following: 10 pmol of each primer for each marker, 200 μM of each deoxyribonucleotide triphosphate (dNTP), 1.2 mM magnesium sulfate (MgSO4), 0.5 units of KOD-plus polymerase (Toyobo, Tokyo, Japan), 1× reaction buffer provided by the manufacturer, and 10 ng genomic DNA. Reactions were performed in a 96-well plate in an iCycler Thermal Cycler (Bio-Rad Laboratories, Hercules, CA, USA) using the following conditions: initial denaturation at 94°C for 2 min; 30 cycles at x°C (x was the denaturizing temperature shown in Table 1) for 15 s, at y°C (y was the annealing temperature shown in Table 1) for 30 s, at 68°C for 30 s; and a final elongation at 68°C for 4 min 30 s. The PCR products were purified using High Pure 96 UF Cleanup Plates (Roche Diagnostics GmbH, Mannheim, Germany) and subjected to cycle sequencing using a BigDye Terminator Cycle Sequencing FS Ready Reaction kit (Perkin-Elmer Applied Biosystems, Foster City, CA, USA) and the same primers as used for the PCR amplification. Sequencing was performed in both directions and each sequencing product was run in an ABI 3130 Genetic Analyzer (Perkin-Elmer Applied Biosystems). The MLNR gene haplotypes in the F2 intercross population were then identified.

Polymerase chain reaction primers for the amplification of the 2 exons in the gene encoding the motilin receptor

Statistical analyses

A mixed-inheritance animal model was used to evaluate the effects of haplotypes on traits. The ld_fix option of the Qxpak software (Pérez-Enciso and Misztal, 2004) was used for the MLNR haplotypes. For statistical models, the additive effect of each allele and the sex effect were used as fixed effects. In this analysis, the infinitesimal genetic effect was included and treated as a random effect with covariance matrix Aσu2 (where A is the numerator relationship matrix). Likelihood ratio tests were performed by not including the additive effects in the model, and nominal p-values were obtained by assuming a x2 distribution of the likelihood ratio test.

Percentage of haplotype variance explained by the model was calculated as

Variance percentage=100×(R_variance-F_variance)/R_variance

where R_variance is the residual and polygenic variances from the reduced model, omitting haplotype effect but including sex and polygenic effects, and F_variance is the residual variance from the full model, including haplotype, sex, and polygenic effects.

RESULTS

On sequencing of the MLNR gene, we found GGGGC insertion polymorphisms located 54 bp upstream of the 5′-untranslated region and the sequenced region differed by 6 single nucleotide polymorphisms (SNPs) from the corresponding region in the draft sequence of the chicken genome. Resultantly, three MLNR haplotypes (Haplotypes 1–3) were identified in the parent individuals (Table 2). The nucleotide sequences of the three haplotypes have been registered in the DNA data bank of Japan (DDBJ) with the following accession numbers: AB703068, AB703069, and AB703070. The results indicated that these haplotypes might be segregated and individuals with six different genotypes (1/1, /2/2, 3/3, 1/2, 1/3, and 2/3) might appear in the F2 population. Therefore, a mismatch amplification mutation assay (MAMA) PCR protocol was developed that detects the 6 genotypes described by Cha et al. (1992). We designed four PCR primers to distinguish MLNR haplotypes and PCR and genotyping were performed as described in Table 3 and Figure 2.

Genotypes of parent individuals of the F2 populations

Polymerase chain reaction (PCR) primers for the mismatch amplification mutation assay (MAMA)

Figure 2

Typical patterns of 2% agarose gel electrophoresis of MAMA results. (A) MAMA results at the SNP site of AB703068: g.1008 C>T. MLNR-1008-13 and MLNR -1008-2 primer sets were used in lane 1 and 2, and in lane 3 and 4, respectively. Lane M, DNA size marker; lane 1 and 3, patterns of Haplotype 2; lane 2 and 4, patterns of Haplotype 1 and 3. (B) MAMA results at the SNP site of AB703068: g.1029 A>G. MLNR-1029-1 and MLNR-1029-23 primer sets were used in lane 5 and 6, and in lane 7 and 8, respectively. Lane M, DNA size marker; lane 5 and 7; pattern of Haplotype 2 and 3; lane 6 and 8; patterns of Haplotype 1.

We observed an association between three MLNR haplotypes (Haplotype 1, 2, and 3) and growth traits (Table 4). There was no association between these haplotypes and growth traits in the overall F2 population. When data were analyzed separately by sex, a significant (p<0.05) association was found between these haplotypes and BW 10 wk, BW 14 wk, ADG 4–10 wk, ADG 10–14 wk, and ADG 0–14 wk in F2 females, whereas there was no association between these haplotypes and growth traits in F2 males.

Phenotypic values of growth traits and effects of Motilin receptor haplotypes on the traits in the F2 population

DISCUSSION

Motilin, a 22 amino acid polypeptide first discovered in the mucosa of upper small intestine of the hog (Brown et al., 1971), is predominantly recognized for its prokinetic effects on gastrointestinal motility and emptying (Poitras et al., 1994). Motilin has been identified in many mammals (Bond et al., 1988; Banfield et al., 1992; Xu et al., 2001; Strausberg et al., 2002; Xu e al., 2003; Ohshiro et al., 2008) and chickens (DeClercq et al., 1996). Motilin binding sites are present in the gastrointestinal tract and central nervous system, along with a G-protein-coupled receptor, GPR38, which has been identified as MLNR (Feighner et al., 1999). Initially, GPR38/MLNR was identified as an orphan receptor belonging to the GHSR family. This family contains receptors for ghrelin (Kojima et al., 1999), neuromedin U (Howard et al., 2000), and neurotensin (Vincent et al., 1999). Of the receptors, the GHSR is most homologous to the MLNR; the chicken forms share 53% identical amino acids (Yamamoto et al., 2008).

It has been shown that motilin and ghrelin have contractile activities in the gastrointestinal tract and their receptors are presented in the gastrointestinal tract of chickens (Kitazawa et al., 1997; Tanaka et al., 2003), although detail distributions of MLNR and GHSR have not been reported. DeClercq et al. (1996) reported that motilin was most abundant in the mucosa of the colon, followed by the jejunum, duodenum, and ileum; low motilin contents were found in the proventriculus and gizzard. Neglia et al. (2004) reported that the largest numbers of ghrelin-positive cells were located at the base of lobuli of the proventriculus gland, along villi of the intestines and in crypts of the duodenum, while lower numbers of ghrelin-immunostained cells were located in crypts of jejunum and ileum and only few ghrelin-immunostained cells were detected at the base of crypts of the large intestine. Kitazawa et al. (2007) reported that chicken ghrelin caused contraction of chicken gastrointestinal smooth muscles in a region dependent manner (crop>colon>esophagus>proventriculus>jejunum> ileum>duodenum), and chicken motilin caused contraction of chicken gastrointestinal tract in a region dependent manner (duodenum>jejunum>ileum>colon>proventriculus >crop>esophagus). These data suggest that ghrelin mainly regulates contractions of the upper (crop, esophagus, and proventriculus) gastrointestinal tract, motilin principally regulates small intestinal motility in chicken. Since it has been reported that motilin can stimulate the GHSR but ghrelin does not activate MLNR in rabbits (Dass et al., 2003) and chickens (Yamamoto et al., 2008), motilin abundantly presented in colon may activate GHSR and lead gastrointestinal mortality in colon. Taken together, these reports suggest that motilin and MLNR play central roles in the regulation of gastrointestinal motility and gastric emptying with ghrelin and GHSR in chicken.

The association between haplotypes and growth traits cannot be explained by the protein sequences, since the SNPs detected in the coding regions of MLNR do not induce missense mutations, and thus MLNR proteins are functionally normal. Therefore, at present, it is unknown whether MLNR haplotypes have any functional significance. The associations detected in this study may be produced by linkage disequilibrium between these haplotypes and another linked gene directly involved in the regulation of these traits. Recently, Zhang et al. (2011) proposed positional candidate genes for growth traits in the same region on chromosome 1, namely, retinoblastoma 1 (RB1) (chr 1: 173.0 Mb); lysophosphatidic acid receptor 6 (chr 1: 173.0 Mb); fibronectin type III domain containing 3A (chr 1: 173.3 Mb); MLNR (chr 1: 173.4 Mb); and calcium binding protein 39-like (chr 1: 173.5 Mb). Of the genes, Zhang et al. (2011) suggested that two SNPs of RB1 gene might be quantitative trait nucleotides for BW; however they have not proposed any hypothetical mechanisms why the SNPs affect BW.

One likely possibility to explain the association between MLNR haplotypes and growth traits is that presence of GGGGC insertion polymorphism in the 5′-upstream region of MLNR, corresponding to the predicted Sp1 binding site, may affect MLNR gene expression. Sp1 is a transcription factor that plays an important role in many cellular activities such as metabolism, cell growth, differentiation, angiogenesis and apoptosis (Deniaud et al., 2009). Sp1 contains three Cys2His2-type zinc finger motifs and it is generally accepted that its mode of action is mediated mainly by binding to the decanucleotide consensus sequence 5′-(G/T)GGGCGG(G/A)(G/A)(C/T)-3′ in double-stranded DNA (dsDNA) (Song et al., 2001). In the predicted Sp1 binding site of MLNR gene, Haplotype 1, 2, and 3 have (GGGGC)3, (GGGGC)4, and (GGGGC)5, respectively, while the draft sequence of the chicken genome has (GGGGC)2. Since it was reported that Sp1-tandem repeat polymorphisms affect binding capabilities of Sp1 and transcriptional efficiencies (In et al., 1997; Sigurdsson et al., 2008), the repeat polymorphisms in MLNR gene promoter may affect MNLR gene expression. Our data suggest that shorter GGGGC-repeats shows better growth performance rather than longer ones. Therefore, further studies are needed to elucidate the differences in transcriptional and translational efficiencies caused by haplotypes.

It is also unclear the reason that MLNR haplotypes are associated with gender differences in growth traits. Motilin and MLNR genes in mice and rats exist as pseudogenes and are thought to genetically knockdown; however, it is still debated whether or not native MLNR is expressed (He et al., 2010). Recently, Miegueu et al. (2011) demonstrated that motilin stimulates preadipocyte proliferation and differentiation and adipocyte lipid storage in murine 3T3-L1 and primary rat adipocytes. Previous reports indicated that acylated and unacylated ghrelin might influence preadipocyte proliferation (Zhang et al., 2004) and could enhance preadipocyte differentiation (Miegueu et al., 2011), thus, motilin may act via GHSR in the murine and rat adipocytes. Zigman (Zigman et al., 2005) reported that female GHSR-null mice were significantly lighter than wild-type controls between 12 and 19 weeks of age, while male GHSR-null and wild-type mice did not differ in body weight at any weeks of age. After 19 weeks of age, female GHSR-null mice had 35.6% less body fat than wild-type mice. Taken together, these data may imply that MLNR plays a key role in fat deposition and body weight especially in females; however, further studies needed to clarify the reason how MLNR affect fat traits in chicken.

In conclusion, we observed significant associations between MLNR haplotypes and growth traits in an F2 population produced by crossing low- and high-growth lines of the Hinai-dori breed of chicken. Our data suggest that haplotype-assisted selection is a potentially useful approach in the Hinai-dori breed. In the future, the associations between haplotypes and growth traits should be studied in other chicken lines and breeds.

ACKNOWLEDGEMENTS

This work was financially supported by the Integrated Research Project for Plant, Insect, and Animal Using Genome Technology from the Ministry of Agriculture, Forestry and Fisheries of Japan.

References

Abasht B, Dekkers JCM, Lamont SJ. 2006;Review of quantitative trait loci identified in the chicken. Poult Sci 85:2079–2096.
Andersson L. 2001;Genetic dissection of phenotypic diversity in farm animals. Nat Rev Genet 2:130–138.
Banfield DK, MacGillivray RTA, Brown JC, McIntosh CHS. 1992;The isolation and characterization of rabbit motilin precursor cDNA. Biochim Biophys Acta 1131:341–344.
Bond CT, Nilaver G, Godfrey B, Zimmerman EA, Adelman JP. 1988;Characterization of complementary deoxyribonucleic acid for precursor of porcine motilin. Mol Endocrinol 2:175–180.
Brown JC, Mutt V, Dryburgh JR. 1971;The further purification of motilin, a gastric motor activity stimulating polypeptide from the mucosa of the small intestine of hogs. Can J Physiol Pharmacol 49:399–405.
Cao ZP, Wang SZ, Wang QG, Wang YX, Li H. 2007;Association of spot14α gene polymorphisms with body weight in the chicken. Poult Sci 86:1873–1880.
Cha RS, Zarbl H, Keohavong P, Thilly WG. 1992;Mismatch amplification mutation assay (MAMA): application to the c-H-ras gene. PCR Methods Appl 2:14–20.
Dass NB, Hill J, Muir A, Testa T, Wise A, Sanger GJ. 2003;The rabbit motilin receptor: molecular characterization and pharmacology. Br J Pharmacol 140:948–954.
DeClercq P, Depoortere I, Macielag M, Vandermeers A, Vandermeers-Piret MC, Peeters TL. 1996;Isolation, sequence, and bioactivity of chicken motilin. Peptides 17:203–208.
Deniaud E, Baguet J, Chalard R, Blanquier B, Brinza L, Meunier J, Michallet MC, Laugraud A, Ah-Soon C, Wierinckx A, Castellazzi M, Lachuer J, Gautier C, Marvel J, Leverrier Y. 2009;Overexpression of transcription factor Sp1 leads to gene expression perturbations and cell cycle inhibition. PLoS ONE 4:e7035.
Ensembl Genome Browser. 2004. European Bioinformatics Institute; UK: http://ensembl.genomics.org.cn/. Accessed June 2013.
Fang M, Nie Q, Luo C, Zhang D, Zhang X. 2010;Associations of GHSR gene polymorphisms with chicken growth and carcass traits. Mol Biol Rep 37:423–428.
Feighner SD, Tan CP, McKee KK, Palyha OC, Hreniuk DL, Pong SS, Austin CP, Figueroa D, MacNeil D, Cascieri MA, Nargund R, Bakshi R, Abramovitz M, Stocco R, Kargman S, O’Neill G, Van Der Ploeg LH, Evans J, Patchett AA, Smith RG, Howard AD. 1999;Receptor for motilin identified in the human gastrointestinal system. Science 284:2184–2188.
Garthwaite TL. 1985;Peripheral motilin administration stimulates feeding in fasted rats. Peptides 6:41–44.
Han RL, Lia ZJ, Lia MJ, Lib JQ, Lana XY, Sunb GR, Kangb XT, Chena H. 2011;Novel 9-bp indel in visfatin gene and its associations with chicken growth. Br Poult Sci 52:52–57.
He J, Irwin DM, Chen R, Zhang YP. 2010;Stepwise loss of motilin and its specific receptor genes in rodents. J Mol Endocrinol 44:37–44.
Howard AD, Wang R, Pong SS, Mellin TN, Strack A, Guan XM, Zeng Z, Williams DL Jr, Feighner SD, Nunes CN, Murphy B, Stair JN, Yu H, Jiang Q, Clements MK, Tan CP, McKee KK, Hreniuk DL, McDonald TP, Lynch KR, Evans JF, Austin CP, Caskey CT, Van der Ploeg LH, Liu Q. 2000;Identification of receptors for neuromedin U and its role in feeding. Nature 406:70–74.
Hu ZL, Fritz ER, Reecy JM. 2007;AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond. Nucl Acids Res 35:604–609.
In KH, Asano K, Beier D, Grobholz J, Finn PW, Silverman EK, Silverman ES, Collins T, Fischer AR, Keith TP, Serino K, Kim SW, De Sanctis GT, Yandava C, Pillari A, Rubin P, Kemp J, Israel E, Busse W, Ledford D, Murray JJ, Segal A, Tinkleman D, Drazen JM. 1997;Naturally occurring mutations in the human 5-lipoxygenase gene promoter that modify transcription factor binding and reporter gene transcription. J Clin Invest 99:1130–1137.
Kitazawa T, Kaiya H, Taneike T. 2007;Contractile effects of ghrelin-related peptides on the chicken gastrointestinal tract in vitro. Peptides 28:617–624.
Kitazawa T, Taneike T, Ohga A. 1997;Functional characterization of neural and smooth muscle motilin receptors in the chicken proventriculus and ileum. Regul Pept 71:87–95.
Kojima M, Hosoda H, Date Y, Nakazato M, Matsuo H, Kangawa K. 1999;Ghrelin is a growth-hormone-releasing acylated peptide from stomach. Nature 402:656–660.
Lei M, Peng X, Zhou M, Luo C, Nie Q, Zhang X. 2008;Polymorphisms of the IGF1R gene and their genetic effects on chicken early growth and carcass traits. BMC Genet 9:70.
Lei MM, Nie QH, Peng X, Zhang DX, Zhang XQ. 2005;Single nucleotide polymorphisms of the chicken insulin-like factor binding protein 2 gene associated with chicken growth and carcass traits. Poult Sci 84:1191–1198.
Leng L, Wang S, Li Z, Wang Q, Li H. 2009;A polymorphism in the 3′-flanking region of insulin-like growth factor binding protein 2 gene associated with abdominal fat in chickens. Poult Sci 88:938–942.
Li ZH, Li H, Zhang H, Wang SZ, Wang QG, Wang YX. 2006;Identification of a single nucleotide polymorphism of the insulin-like growth factor binding protein 2 gene and its association with growth and body composition traits in the chicken. J Anim Sci 84:2902–2906.
Miegueu P, Cianflone K, Richard D, St-Pierre DH. 2011;Motilin stimulates preadipocyte proliferation and differentiation and adipocyte lipid storage. Am J Physiol Endocrinol Metab 301:E758–766.
Neglia S, Arcamone N, Esposito V, Gargiulo G. 2004;Ghrelin in the gastroenteric tract of birds: immunoreactivity expression. Vet Res Commun 28:213–215.
Nie Q, Fang M, Xie L, Zhou M, Liang Z, Luo Z, Wang G, Bi W, Liang C, Zhang W, Zhang X. 2008;The PIT1 gene polymorphisms were associated with chicken growth traits. BMC Genet 9:20.
Ohshiro H, Nonaka M, Ichikawa K. 2008;Molecular identification and characterization of the dog motilin receptor. Regul Pept 146:80–87.
Pérez-Enciso M, Misztal I. 2004;Qxpak: a versatile mixed model application for genetical genomics and QTL analyses. Bioinformatics 20:2792–2798.
Poitras P, Miller P, Gagnon D, Stpierre S. 1994;Motilin synthetic analogues and motilin receptor antagonists. Biochem Biophys Res Commun 205:449–454.
Qu JT, Tang SQ, Sun DX, Zhang Y. 2009;Polymorphisms of three neuroendocrine-correlated genes associated with growth and reproductive traits in the chicken. Poult Sci 88:722–727.
Rikimaru K, Takahashi H. 2007;A method for discriminating a Japanese brand of chicken, the Hinai-jidori, using microsatellite markers. Poult Sci 86:1881–1886.
Rikimaru K, Sasaki O, Koizumi N, Komatsu M, Suzuki K, Takahashi H. 2011;Mapping of quantitative trait loci affecting growth traits in a Japanese native chicken cross. Asian-Aust J Anim Sci 24:1329–1334.
Rikimaru K, Komatsu M, Suzuki K, Uemoto Y, Takeda H, Takahashi H. 2012;Association between cholecystokinin type A receptor haplotypes and growth traits in Japanese Hinai-dori crossbred chickens. Mol Biol Rep 39:4479–4484.
Rikimaru K, Takeda H, Uemoto Y, Komatsu M, Takahasi K, Suzuki D, Takahashi H. 2013;Effect of a single-nucleotide polymorphism in the cholecystokinin type A receptor gene on growth traits in the Hinai-dori chicken breed. J Poult Sci 50:206–211.
Rosenfeld DJ, Garthwaite TL. 1987;Central administration of motilin stimulates feeding in rats. Physiol Behav 39:753–756.
Science Council of Japan. 2006. Guidelines for proper conduct of animal experiments Tokyo, Japan:
Sigurdsson S, Göring HH, Kristjansdottir G, Milani L, Nordmark G, Sandling JK, Eloranta ML, Feng D, Sangster-Guity N, Gunnarsson I, Svenungsson E, Sturfelt G, Jönsen A, Truedsson L, Barnes BJ, Alm G, Rönnblom L, Syvänen AC. 2008;Comprehensive evaluation of the genetic variants of interferon regulatory factor 5 (IRF5) reveals a novel 5 bp length polymorphism as strong risk factor for systemic lupus erythematosus. Hum Mol Genet 17:872–881.
Song C, Gu X, Feng C, Wang Y, Gao Y, Hu X, Li N. 2011;Evaluation of SNPs in the chicken HMGA2 gene as markers for body weight gain. Anim Genet 42:333–336.
Song J, Ugai H, Kanazawa I, Sun K, Yokoyama KK. 2001;Independent repression of a GC-rich housekeeping gene by Sp1 and MAZ involves the same cis-elements. J Biol Chem 276:19897–19904.
Strausberg RL, Feingold EA, Grouse LH, Derge JG, Klausner RD, Collins FS, Wagner L, Shenmen CM, Schuler GD, Altschul SF, et al. 2002;Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences. Proc. Natl. Acad. Sci USA 99:16899–16903.
Tanaka M, Miyazaki T, Yamamoto I, Nakai N, Ohta Y, Tsushima N, Wakita M, Shimada K. 2003;Molecular characterization of chicken growth hormone secretagogue receptor gene. Gen Comp Endocrinol 134:198–202.
UCSC Genome Browser Home. 2004. University of California; Santa Cruz, USA: http://genome.ucsc.edu. Accessed June 2013.
Uemoto Y, Sato S, Ohtake T, Sato S, Okumura Y, Kobayashi E. 2011;Ornithine decarboxylase gene is a positional candidate gene affecting growth and carcass traits in F2intercross chickens. Poult Sci 90:35–41.
Vincent JP, Mazella J, Kitabgi P. 1999;Neurotensin and neurotensin receptors. Trends Pharmacol Sci 20:302–309.
Xu L, Depoortere I, Tang M, Peeters TL. 2001;Identification and expression of the motilin precursor in the guinea pig. FEBS Lett 490:7–10.
Xu L, Depoortere I, Thielemans L, Huang Z, Tang M, Peeters TL. 2003;Sequence, distribution and quantification of the motilin precursor in the cat. Peptides 24:1387–1395.
Yamamoto I, Kaiya H, Tsutsui C, Sakai T, Tsukada A, Miyazato M, Tanaka M. 2008;Primary structure, tissue distribution, and biological activity of chicken motilin receptor. Gen Comp Endocrinol 156:509–514.
Zhang H, Liu SH, Zhang Q, Zhang YD, Wang SZ, Wang QG, Wang YX, Tang ZQ, Li H. 2011;Fine-mapping of quantitative trait loci for body weight and bone traits and positional cloning of the RB1 gene in chicken. J Anim Breed Genet 128:366–375.
Zhang W, Zhao L, Lin TR, Chai B, Fan Y, Gantz I, Mulholland MW. 2004;Inhibition of adipogenesis by ghrelin. Mol. Biol Cell 15:2484–2491.
Zhang ZR, Liu YP, Yao YG, Jiang XS, Du HR, Zhu Q. 2009;Identification and association of the single nucleotide polymorphisms in calpain3 (CAPN3) gene with carcass traits in chickens. BMC Genet 10:10.
Zhou H, Mitchell AD, McMurtry JP, Ashwell CM, Lamont SJ. 2005;Insulin-like growth factor-I gene polymorphism associations with growth, body composition, skeleton integrity, and metabolic traits in chickens. Poult Sci 84:212–219.
Zigman JM, Nakano Y, Coppari R, Balthasar N, Marcus JN, Lee CE, Jones JE, Deysher AE, Waxman AR, White RD, Williams TD, Lachey JL, Seeley RJ, Lowell BB, Elmquist JK. 2005;Mice lacking ghrelin receptors resist the development of diet-induced obesity. J Clin Invest 115:3564–3572.

Article information Continued

Figure 1

Hinai-dori. (A) Hinai-dori chickens of the Preservation Society (PS) line. Left, male; right, female. (B) Hinai-dori chickens of the Akita Prefectural Livestock Experiment Station (LES) line. Left, female; right, male.

Figure 2

Typical patterns of 2% agarose gel electrophoresis of MAMA results. (A) MAMA results at the SNP site of AB703068: g.1008 C>T. MLNR-1008-13 and MLNR -1008-2 primer sets were used in lane 1 and 2, and in lane 3 and 4, respectively. Lane M, DNA size marker; lane 1 and 3, patterns of Haplotype 2; lane 2 and 4, patterns of Haplotype 1 and 3. (B) MAMA results at the SNP site of AB703068: g.1029 A>G. MLNR-1029-1 and MLNR-1029-23 primer sets were used in lane 5 and 6, and in lane 7 and 8, respectively. Lane M, DNA size marker; lane 5 and 7; pattern of Haplotype 2 and 3; lane 6 and 8; patterns of Haplotype 1.

Table 1

Polymerase chain reaction primers for the amplification of the 2 exons in the gene encoding the motilin receptor

Exon Set Forward primer (5′ → 3′) Reverse primer (5′→3′) Product (bp) Denaturizing temp. (°C) Annealing temp. (°C)
1 1 TCTTTTTCGTCCTTCGTTCCTT CTGAAGTCGGTGCGGTTGT 704 99.5 64.0
1 2 GTCCGACCTGCTCATCCTGCT CTCATCCGTCTCACCTCCCCTT 736 98.0 60.0
2 1 GGGGAGACACCCAAGGTCAGGC CCAAGCATTTCAGCACACAAAGAGT 644 94.0 62.0
2 2 CCAGGACCTACACAGAGACCAG TACAGGGTTGGAGTTAGATGGC 898 94.0 62.0

Table 2

Genotypes of parent individuals of the F2 populations

Individuals No. Genotype
Parent male 1 2/2
2 2/2
3 2/3
Parent female 1 1/2
2 1/2
3 2/2
4 2/2
5 2/2
6 1/2
7 1/2
8 2/2
9 1/1

Table 3

Polymerase chain reaction (PCR) primers for the mismatch amplification mutation assay (MAMA)

Primer sets Primers (5′→3′) PCR product (bp) Target position (AB703068) PCR amplification

Haplotype

1 2 3
MLNR-1008-13 CAGCGTCCTCTACGGCTTCA
GCCGGCGGTCCGTGGtCa
217 g1008. C>T + +
MLNR-1008-2 CAGCGTCCTCTACGGCTTCA
GCCGGCGGTCCGTGGtCg
217 g1008. C>T +
MLNR-1029-1 GCGTCCTCTACGGCTTCATC
AGTCTGTGGAATGGCGCgTc
238 g1029. A>G +
MLNR-1029-23 GCGTCCTCTACGGCTTCATC
AGTCTGTGGAATGGCGCgTt
238 g1029. A>G + +

Bases shown in lower case with an underline represent induced mismatches. Bases shown in lower case at the 3′-end represent target single nucleotide polymorphisms (SNPs). The haplotypes that can or cannot be amplified by PCR for each primer set are shown as ‘+’ or ‘−’, respectively. For the PCRs, we used 15-μL reaction volumes containing the following: 3 pmol of each primer for each marker, 200 μM of each dNTP, 0.5 units of Paq5000DNA Polymerase (Agilent Technologies, La Jolla, CA, USA), 1× reaction buffer (containing 2 mM MgCl2) provided by the manufacturer, and 10 ng genomic DNA, Reactions were performed in a 96-well plate in an iCycler Thermal Cycler (Bio-Rad Laboratories, Hercules, CA, USA) using the following conditions: initial denaturation at 95°C for 2 min; 30 cycles at 94°C for 30 s, at 62°C for 30 s, at 72°C for 1 min (MLNR-1008-13 and −2), or 35 cycles at 94°C for 30 s, at 64°C for 30 s, at 72°C for 1 min (MLNR-1029-1 and −23); and a final elongation at 72°C for 7 min. The PCR products were electrophoresed on a 2.0% agarose gel with 1× Tris-acetate EDTA (TAE) buffer and stained with ethidium bromide. The combination of these results enabled us to identify diplotypes in the F2 population.

Table 4

Phenotypic values of growth traits and effects of Motilin receptor haplotypes on the traits in the F2 population

Trait n Phenotypic values LRT p value Haplotype 1 Haplotype 2 Haplotype 3 Variance (%)


Mean±SD mean±SE mean±SE mean±SE
Total
 BW 0 d (g) 418 34.4±2.4 2.7 0.255 - - - -
 BW 4 wk (g) 417 1 231.1±38.3 0.2 0.888 - - - -
 BW 10 wk (g) 418 960.6±163.5 4.3 0.116 - - - -
 BW 14 wk (g) 418 1,466.9±260.7 5.2 0.075 - - - -
 ADG 0–4 wk (g/d) 417 1 5.2±1.0 0.3 0.856 - - - -
 ADG 4–10 wk (g/d) 417 1 17.4±3.3 5.1 0.077 - - - -
 ADG 10–14 wk (g/d) 418 18.1±4.4 3.6 0.166 - - - -
 ADG 0–14 wk (g/d) 418 14.6±2.7 5.2 0.073 - - - -
Male
 BW 0 d (g) 206 34.6±2.5 5.0 0.041 - - - -
 BW 4 wk (g) 205 1 243.6±38.3 0.3 0.797 - - - -
 BW 10 wk (g) 206 1,053.4±144.8 2.2 0.536 - - - -
 BW 14 wk (g) 206 1,658.8±199.5 1.5 0.670 - - - -
 ADG 0–4 wk (g/d) 205 1 5.5±1.0 0.3 0.835 - - - -
 ADG 4–10 wk (g/d) 205 1 19.3±2.9 2.8 0.485 - - - -
 ADG 10–14wk (g/d) 206 21.6±2.9 0.1 0.970 - - - -
 ADG 0–14 wk (g/d) 206 16.6±2.0 0.8 0.672 - - - -
Female
 BW 0 d (g) 211 34.2±2.3 0.7 0.915 - - - -
 BW 4 wk (g) 211 219.0±34.3 3.2 0.374 - - - -
 BW 10 wk (g) 211 870.8±125.9 8.6 0.021* 462.9±20.2 426.1±15.5 394.4±23.8 7.5
 BW 14 wk (g) 211 1,280.6±157.4 12.0 0.002* 692.6±25.6 629.8±19.5 582.5±30.3 11.0
 ADG 0–4 wk (g/d) 211 4.9±0.9 3.0 0.375
 ADG 4–10wk (g/d) 211 15.5±2.5 8.9 0.041* 8.2±0.4 7.6±0.3 7.0±0.5 5.7
 ADG 10–14 wk (g/d) 211 14.6±2.3 7.5 0.007* 8.2±0.4 7.3±0.3 6.7±0.5 5.5
 ADG 0–14 wk (g/d) 211 12.7±1.6 12.3 0.002* 6.9±0.3 6.3±0.2 5.8±0.3 11.4

BW 0 d, BW 4 wk, BW 10 wk, and BW 14 wk: body weight measured at day 0, 4 weeks, 10 weeks, and 14 weeks of age.

ADG 0–4 wk, ADG 4–10 wk, ADG 10–14 wk, and ADG 0–14 wk: average daily gain between 0 and 4 weeks, between 4 and 10 weeks, between 10 and 14 weeks, and between 0 and 14 weeks.

SD = Standard deviation. SE = Standard error. LRT = Log likelihood ratio test statistics.

1

Since BW 4 wk has a missing datum, ADG 0–4 wk and 4–10 wk have missing data.

*

p<0.05.