- Open Access
Molecular cloning and expression patterns of the cholesterol side chain cleavage enzyme (CYP11A1) gene during the reproductive cycle in goose (Anas cygnoides)
Journal of Animal Science and Biotechnology volume 6, Article number: 54 (2015)
CYP11A1, a gene belonging to the family 11 of cytochrome P450, encodes a crucial steroidogenic enzyme that catalyzes the initial step in the production of all classes of steroids. Many studies show that CYP11A1 plays a role in ovary function. However, the role of CYP11A1 in goose reproductive cycle remains largely unknown.
In this study, full-length CYP11A1 cDNA of Zhedong goose was obtained using reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE). The cDNA consisted of a 96-base pair (bp) 5′untranslated region (UTR), a 179-bp 3′UTR and a 1509-bp open reading frame. The open reading frame encodes a putative 503 amino acid protein that shares high homology with CYP11A1 of other birds. The amino acid sequence possesses conserved domains of the P450 superfamily, which include the steroid-binding domain and the heme-binding region. Real-time quantitative polymerase chain reaction (qPCR) analysis revealed CYP11A1 mRNA was expressed ubiquitously in every Zhedong goose tissue analyzed, including the heart, liver, glandular stomach, lung, spleen, kidney, intestinum tenue, intestinum crassum, cerebrum, cerebellum, muscle, oviduct, pituitary, hypothalamus and ovary.. The relatively low levels of CYP11A1 mRNA were detected in pituitary, ovary and oviduct tissues at ovulation when compared with levels at oviposition. Interestingly, higher expression was observed in ovary and oviduct tissues during brooding. Lastly, higher mRNA expression of Yangzhou geese was detected during the ovulation period than that of Zhedong geese.
Our findings reveal the sequence characterization and expression patterns of the CYP11A1 gene during the goose reproductive cycle, which may provides correlative evidence that CYP11A1 expression is important in reproduction activity.
Cytochrome P450 side chain cleavage (P450scc, encoded by the CYP11A1 gene) plays a major role in the regulation of steroidogenesis by mediating the conversion of cholesterol to pregnenolone . Because of its physiological importance, CYP11A1 has been the focus of many studies. Recently, more and more evidence shows that CYP11A1 plays a role in ovary function. In fowl, many investigations have found that CYP11A1 expression changes dynamically in developing follicles [2–4] and that its expression in ovarian granulosa cells is influenced largely by hormonal regulation [5–9]. Another study found that the hormonally and developmentally regulated expression of CYP11A1 was principally driven by multiple trans-acting factors, like Sp-1 [10, 11], SF-1/LRH-1, GATA4, CREB-1,AP-1  AP-2, LBP-1b/LBP-9  and Ff1b .
The goose (Anas cygnoides) is a commercially important food source that is widely cultivated in China. It is an ideal animal model for characterization of fowl reproduction because of obvious reproductive stages and strong broodiness . In a previous study, we identified CYP11A1 as an important candidate gene through transcriptome profiling of ovaries from laying and brooding geese . To further characterize the role of CYP11A1 in the goose reproductive cycle, we cloned Zhedong goose (Anas cygnoides) CYP11A1 and characterized its spatiotemporal expression patterns by Real-time qPCR. We also investigated the differences in CYP11A1 expression in high egg production breed(Yangzhou geese, Anas cygnoides) and low egg production breed(Zhedong goose). These data may facilitate a better understanding of expression patterns of the CYP11A1 gene during the reproductive cycle in goose.
All animal experiments were reviewed and approved by the Institutional Animal Care and Use Committee of Yangzhou University. Procedures were performed in accordance with the Regulations for the Administration of Affairs Concerning Experimental Animals (Yangzhou University, China, 2012) and the Standards for the Administration of Experimental Practices (Jiangsu, China, 2008).
Animals and tissue sample collection
One hundred female Zhedong geese and Yangzhou geese were selected from the breeding farm of Jiangsu Lihua Animal Husbandry Co. Ltd (Changzhou, China)), which laid about thirty and seventy eggs, respectively. They were both housed semi-enclosed building and raised the same conditions according to the farm’s standard practice. During the experiment, geese were fed ad libitum with rice grain supplemented with green grass or water plants whenever possible. The feed was provided during the daytime when the geese were released to an open area. The geese were exposed to natural light and ambient temperature throughout this study. The five 120-day-old young (prelaying) Zhedong female geese were sacrificed to investigate CYP11A1 expression patterns in different tissues. The Zhedong geese were anesthetized with sodium pentobarbital and the heart, liver, glandular stomach, lung, spleen, kidney, intestinum tenue, intestinum crissum, cerebrum, cerebellum, muscle, oviduct(infundibulum), pituitary, hypothalamus and ovary(stroma) were removed and immediately frozen in liquid nitrogen before storage at−80 °C for RNA isolation. Tissue samples from the oviduct(infundibulum), pituitary gland, hypothalamus, and ovary(stroma) were obtained from twelve sacrificed 380-day-old adult geese, including three laying Zhedong geese with an egg in the oviduct (ovulation, the release of an ovum from a ruptured follicle), three laying Zhedong geese without an egg in the oviduct (oviposition, the laying of the egg), three brooding Zhedong geese, and three laying Yangzhou geese with an egg in the oviduct (ovulation). The tissue samples were obtained promptly, as described above, to characterize developmental expression patterns.
Zhedong goose CYP11A1 cDNA cloning and sequencing
Total RNA was extracted from collected tissue samples using TRIzol reagent according to the manufacturer’s instruction (TaKaRa, China). The RNA was resuspended in RNase-free water, and the concentration and purity were determined using a NanoDrop Spectrophotometer (NanoDrop, USA). After purification, total RNA (2 μg) was reverse transcribed using M-MLV reverse transcriptase (Promega, USA) according to the manufacturer’s protocol. BLAST analysis of one unigene revealed that it was highly similar to the chicken CYP11A1 gene . Primers were designed according the unigene (Table 1) and RT-PCR was performed using ovarian cDNA from geese. The PCR product was purified, cloned into the pMD19-T vector (TaKaRa, China), and subjected to sequence analysis. The 5′ and 3′ ends of CYP11A1 were amplified via rapid amplification of cDNA ends (RACE) using the SMART RACE cDNA amplification protocol (Clontech, USA) and the 3′-Full RACE Kit (TaKaRa, China), respectively. RACE primers (Table 1) were designed using the partial CYP11A1 nucleotide sequence obtained from RT-PCR. Touchdown and nested PCRs were performed according to the manufacturer’s instructions. Amplicons were then cloned into a plasmid vector for nucleotide sequencing by Sangon Biotech (Shanghai, China).
The Zhedong goose CYP11A1 cDNA and the deduced amino acid sequences were analyzed using DNAssist (version 2.2) and the Expasy search program (http://au.expasy.org/tools/), respectively. Homology analyses were carried out using Clustal W (http://www.ebi.ac.uk/Tools/msa/). The phylogenetic analyses and statistical neighbor-joining bootstrap tests of the phylogenies were performed using the MEGA package (version 6.0) and the bootstrap method was used to indicate confidence values for tree nodes based on 1000 replicates. The conserved domain analysis was predicted by SMART (http://smart.embl-heidelberg.de/).
CYP11A1 mRNA expression patterns in Zhedong goose
To study CYP11A1 mRNA expression, we performed Real-time qPCR on total RNA isolated from the tissues of prelaying Zhedong geese and the kidney tissue served as a calibrator tissue. Assays were conducted in 20 μL reactions using the SYBR Premix Ex TaqTM (TaKaRa, China) and performed on an ABI two-step RT-PCR system (Applied Biosystems 7500, USA) with diluted first-strand cDNA. qPCR programs for CYP11A1 and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were: 1 cycle of 95 °C for 5 min,40 cycles of 95 °C for 10 s, 60 °C for 34 s data collection, followed by 1 cycle for melting curve analysis. All cDNA synthesis reactions were carried out using 100 ng of total RNA per reaction and assayed in three to four technical replicates for each set of biological samples. The same methods were used to determine the Zhedong CYP11A1 mRNA expression profile during the reproductive cycle. The GAPDH gene served as an internal reference gene and the hypothalamus tissue from prelaying Zhedong geese served as a calibrator tissue for differential expression analyses. To analyze the difference of expression patterns between different egg production goose breeds(Zhedong geese and Yangzhou geese), the mean ΔCt value of the hypothalamus tissue of Zhedong geese within each group was used as the calibrator.
Statistical analyses were performed using t-tests independent of group and SPSS version 13.0 software.
Zhedong goose CYP11A1 cDNA cloning and sequence analyses
To obtain the full-length cDNA of goose CYP11A1, RT-PCR and RACE were carried out. The CYP11A1 cDNA was only found to be 1784 nucleotides in length, which included a 96-nucleotide 5′ UTR, a 179-nucleotide 3′ UTR, and a 1509-nucleotide open reading frame putatively encoding a single 503 amino acid protein. The other transcript variant did not been found among the tissues tested in Zhedong goose.
Phylogenetic analysis of the putative CYP11A1
To evaluate the relationship between goose CYP11A1 and that of other birds, sequence alignment and phylogenetic analysis were carried out. Alignment of the amino acid residues of goose CYP11A1 with those of other birds forms is shown in Fig. 1. Goose CYP11A1 shared fairly high identity with birds orthologs (over 69 % identity), including budgerigar (86.3 % identity), chicken (85.6 % identity), penguin (85.2 % identity), turkey (84.8 % identity), rock pigeon (84.8 % identity), peregrine (84.7 % identity), crested ibis (84.4 % identity), peregrine (83.7 % identity), Tibetan ground-tit (77.0 % identity), Zebra Finch (76.8 % identity), and Medium ground finch (69.7 % identity). Furthermore, the p450 domain (aa 38–495 in geese) was conserved in all birds sequences analyzed.
Phylogenetic analysis of CYP11A1 resulted in clear segregation into three groups, one branch containing mammalian CYP11A1, one branch of birds CYP11A1, another branch of fish CYP11A1(Fig. 2). Goose CYP11A1 clustered with other Phasianidae forms and appears to be significantly diverse from Paridae forms.
Expression pattern of CYP11A1 in different tissues and reproductive cycle stages in Zhedong goose
The qPCR demonstrated that CYP11A1 was ubiquitously expressed in fifteen tissues tested, but the expression levels were distinctly different (Fig. 3). High levels of CYP11A1 transcript were detected in ovary, oviduct, pituitary, hypothalamus, lung, and spleen tissues in Zhedong goose, and lower negligible expression levels were found in heart and muscle tissues (Fig. 3).
To furtherly determine temporal expression patterns we characterized CYP11A1 expression in the ovary, oviduct, pituitary, and hypothalamus tissues during different stages of the goose reproductive cycle (prelaying, ovulation, oviposition, and broody periods). We observed different CYP11A1 expression profiles in the different tissues at different times of the reproductive cycle (Figs. 4). From an overall perspective, it was relatively low during the ovulation period (Fig. 4). In the pituitary tissue, CYP11A1 mRNA levels are relatively high during the prelaying period, and steadily decline through the egg laying period. In the ovary, CYP11A1 mRNA levels dropped to an approximate tenfold decrease during the ovulation period followed by a significant increase during the oviposition period. In the oviduct, CYP11A1 mRNA levels were relatively high during the broody period and relatively low during the prelaying period. In the hypothalamus, CYP11A1 mRNA levels remained low throughout the entire reproductive cycle.
Comparison on CYP11A1 tissue-differential expression between Zhedong geese and Yangzhou geese in ovulation period
Additionally, to explore the differences in CYP11A1 expression in different goose breeds, its expression levels were compared in Zhedong geese (low egg production) and Yangzhou geese (high egg production) in ovulation periods . The expression of CYP11A was significantly higher in the examined reproductive tissues of Yangzhou geese in the ovulation period in pituitary and ovary tissues than that of Zhedong geese (Fig. 5).
In this study, cDNAs encoding CYP11A1 were cloned from Chinese goose ovaries. The deduced amino acid sequence of CYP11A1 shows high identity with that of other birds and contains some structural domains typical of the cytochrome P450 superfamily, specifically a steroid-binding domain and a heme-binding domain . Both domains appear to be well conserved among the different species, indicating that they may be critical for proper enzymatic activity of all P450scc.
As expected on the basis of the predicted goose CYP11A1 amino acid sequence, phylogenetic analysis clearly clustered goose CYP11A1 with other Phasianidae forms and showed that birds CYP11A1 appear to be significantly diverged from mammalian forms and fish forms.
CYP11A1 mRNA was expressed ubiquitously in every tissue analyzed except heart and muscle, albeit to different degrees. As expected, CYP11A1 mRNA was easily detected in reproductive tissues and endocrine tissues. These results are in agreement with recent molecular studies aimed at determining the temporal expression of CYP11A1 in mammals [19, 20]. The expression of the CYP11A1 was previously well documented in tissues such as the adrenals, gonads, placenta [21, 22], pancreas , skin , gut , kidney  and certain areas of the brain , where its expression has been linked to steroidogenic activity. The relative abundance of specific transcripts at these various sites suggests that CYP11A1 could potentially play an important role in regulating local steroid hormone synthesis [27, 28]. We also found CYP11A1 to be expressed in various nonendocrine tissues, including the lung, spleen, stomach, liver, brain, gut and kidney. During the last few years, CYP11A1 has been shown to be expressed in various tissues that are not involved in steroid synthesis. Sigel A et al. have elaborated CYP11A1 own ubiquitous roles .
The qPCR analysis revealed that CYP11A1 transcript abundance drastically changed in ovary, oviduct, and pituitary tissues during the different stages of the reproductive cycle. As expected, CYP11A1 mRNA levels were relatively low at ovulation compared with oviposition. CYP11A1 catalyzes the conversion of cholesterol to pregnenolone in the first step of steroid biosynthesis . It has been reported in geese and hens that the levels of progesterone in the plasma peak 2 to 3 h before ovulation and then decrease [30, 31] which might be due to low CYP11A1 mRNA level. The higher CYP11A1 expression was observed in ovary and oviduct tissues during broodiness, which is in concordance with previous studies that showed that gonadal steroid levels are maintained at a high level in broody birds . Interestingly, The CYP11A1 expression was high in pituitary than the oviduct during prelaying time. In contrast, the expression of gene was high in ovary/oviduct during oviposition, ovulation and brooding, but low in hypothalamus/pituitary (Fig. 4). There may be a feedback mechanism to control this gene/enzyme expression. CYP11A1 is an important enzyme to produce the progestin precursor required for the synthesis of estradiol . We speculate the estradiol might effect on steroid hormone synthesis by hypothalamus-pituitary-gonad (HPG) axis. On the contrary, high concentrations of steroids might regulate CYP11A1 mRNA expression.
We also observed higher CYP11A1 expression in Yangzhou geese than in Zhedong geese during the ovulation period. Yangzhou geese are known for high egg production, while Zhedong geese show low egg production. Yangzhou geese might need to release more pregnenolone to ovulate, so the observed higher CYP11A1 expression levels may be required to ensure sufficient hormone secretion.
In summary, we present the molecular cloning and characterization of the goose CYP11A1, and analyze its expression during the goose reproductive cycle. The relatively low levels of CYP11A1 mRNA were detected at ovulation when compared with levels at oviposition. The higher mRNA expression was investigated during the ovulation period of Yangzhou geese than that of Zhedong geese. Our findings provides correlative evidence that CYP11A1 expression is important in reproduction activity.
side-chain cleavage cytochrome P450
cholesterol side chain cleavage enzyme
reverse transcription polymerase chain reaction
rapid amplification of cDNA ends
quantitative polymerase chain reaction
O’Hara L, York JP, Zhang P, Smith LB. Targeting of GFP-Cre to the mouse Cyp11a1 locus both drives cre recombinase expression in steroidogenic cells and permits generation of Cyp11a1 knock out mice. PLoS One. 2014;9(1), e84541.
Kowalski KI, Tilly JL, Johnson AL. Cytochrome P450 side-chain cleavage (P450scc) in the hen ovary. I. Regulation of P450scc messenger RNA levels and steroidogenesis in theca cells of developing follicles. Biol Reprod. 1991;45(6):955–66.
Tilly JL, Kowalski KI, Johnson AL. Cytochrome P450 side-chain cleavage (P450scc) in the hen ovary. II. P450scc messenger RNA, immunoreactive protein, and enzyme activity in developing granulosa cells. Biol Reprod. 1991;45(6):967–74.
Johnson A, Woods DC. Dynamics of avian ovarian follicle development: cellular mechanisms of granulosa cell differentiation. Gen Comp Endocrinol. 2009;163(1):12–7.
Sechman A, Pawlowska K, Hrabia A. Effect of 3,3′,5-triiodothyronine and 3,5-diiodothyronine on progesterone production, cAMP synthesis, and mRNA expression of STAR, CYP11A1, and HSD3B genes in granulosa layer of chicken preovulatory follicles. Domest Anim Endocrinol. 2011;41(3):137–49.
Rangel PL, Rodriguez A, Rojas S, Sharp PJ, Gutierrez CG. Testosterone stimulates progesterone production and STAR, P450 cholesterol side-chain cleavage and LH receptor mRNAs expression in hen (Gallus domesticus) granulosa cells. Reproduction. 2009;138(6):961–9.
Chiu Chih-Hsien LC-W, Jong D-S, Kuo W-W, Ting LI, Wu L-S. Expression of Steroidogenic Enzymes and Related Protein in Follicles of Tsaiya Duck. J Agr Assoc China. 2007;8(5):456–69.
Chiu CH, Fei CY, Srinivasan R, Wu LS. Inhibitory effects of epidermal growth factor on progesterone production of ovarian granulosa cells in Tsaiya duck (Anas platyrhynchos var. domestica). Br Poult Sci. 2010;51(6):821–7.
Li Z, Johnson A. Regulation of P450 cholesterol side-chain cleavage messenger ribonucleic acid expression and progesterone production in hen granulosa cells. Biol Reprod. 1993;49(3):463–9.
Hammer F, Compagnone NA, Vigne JL, Bair SR, Mellon SH. Transcriptional regulation of P450scc gene expression in the embryonic rodent nervous system. Endocrinology. 2004;145(2):901–12.
Shih M-CM, Chiu Y-N, Hu M-C, Guo C, Chung B-c. Regulation of steroid production: analysis of Cyp11a1 promoter. Mol Cell Endocrinol. 2011;336(1):80–4.
LaVoie HA, King SR. Transcriptional regulation of steroidogenic genes: STARD1, CYP11A1 and HSD3B. Exp Biol Med. 2009;234(8):880–907.
Guo IC, Shih MC, Lan HC, Hsu NC, Hu MC, Chung BC. Transcriptional regulation of human CYP11A1 in gonads and adrenals. J Biomed Sci. 2007;14(4):509–15.
Quek SI, Chan WK. Transcriptional activation of zebrafish cyp11a1 promoter is dependent on the nuclear receptor Ff1b. J Mol Endocrinol. 2009;43(3):121–30.
Xu Q, Zhang Y, Chen Y, Tong YY, Rong GH, Huang ZY, et al. Identification and Differential Expression of microRNAs in Ovaries of Laying and Broody Geese (Anser cygnoides) by Solexa Sequencing. PLoS One. 2014;9(2):e87920.
Xu Q, Zhao WM, Chen Y, Tong YY, Rong GH, Huang ZY, et al. Transcriptome Profiling of the Goose (Anser cygnoides) Ovaries Identify Laying and Broodiness Phenotypes. PLoS One. 2013;8(2):e55496.
Nomura O, Nakabayashi O, Nishimori K, Mizuno S. The cDNA cloning and transient expression of a chicken gene encoding cytochrome P-450scc. Gene. 1997;185(2):217–22.
Ni Y, Zhou Y, Lu L, Grossmann R, Zhao R. Developmental changes of FSH-R, LH-R, ER-beta and GnRH-I expression in the ovary of prepubertal ducks (Anas platyrhynchos). Anim Reprod Sci. 2007;100(3–4):318–28.
Vilchis F, Chavez B, Larrea F, Timossi C, Montiel F. The cDNA cloning and tissue expression of the cytochrome P450scc from Syrian hamster (Mesocricetus auratus). Gen Comp Endocrinol. 2002;126(3):279–86.
Watanabe KP, Kawai YK, Ikenaka Y, Kawata M, Ikushiro S, Sakaki T, et al. Avian cytochrome P450 (CYP) 1–3 family genes: isoforms, evolutionary relationships, and mRNA expression in chicken liver. PLoS One. 2013;8(9), e75689.
Schwartz JR, Roy SK. Developmental expression of cytochrome P450 side-chain cleavage and cytochrome P450 17 alpha-hydroxylase messenger ribonucleic acid and protein in the neonatal hamster ovary. Biol Reprod. 2000;63(6):1586–93.
Slominski A, Semak I, Zjawiony J, Wortsman J, Gandy MN, Li J, et al. Enzymatic metabolism of ergosterol by cytochrome p450scc to biologically active 17alpha,24-dihydroxyergosterol. Chem Biol. 2005;12(8):931–9.
Morales A, Cuellar A, Ramirez J, Vilchis F, Diaz-Sanchez V. Synthesis of steroids in pancreas: evidence of cytochrome P-450scc activity. Pancreas. 1999;19(1):39–44.
Slominski A, Zjawiony J, Wortsman J, Semak I, Stewart J, Pisarchik A, et al. A novel pathway for sequential transformation of 7-dehydrocholesterol and expression of the P450scc system in mammalian skin. Eur J Biochem. 2004;271(21):4178–88.
Cima I, Corazza N, Dick B, Fuhrer A, Herren S, Jakob S, et al. Intestinal epithelial cells synthesize glucocorticoids and regulate T cell activation. J Exp Med. 2004;200(12):1635–46.
Dalla Valle L, Toffolo V, Vianello S, Belvedere P, Colombo L. Expression of cytochrome P450scc mRNA and protein in the rat kidney from birth to adulthood. J Steroid Biochem Mol Biol. 2004;88(1):79–89.
Xiao Y, Ni Y, Huang Y, Wu J, Grossmann R, Zhao R. Effects of kisspeptin-10 on progesterone secretion in cultured chicken ovarian granulosa cells from preovulatory (F1–F3) follicles. Peptides. 2011;32(10):2091–7.
Ahumada-Solórzano SM, Carranza ME, Pedernera E, Rodríguez-Méndez AJ, Luna M, Arámburo C. Local expression and distribution of growth hormone and growth hormone receptor in the chicken ovary: Effects of GH on steroidogenesis in cultured follicular granulosa cells. Gen Comp Endocrinol. 2012;175(2):297–310.
Sigel A, Sigel H, Sigel RK. The ubiquitous roles of cytochrome P450 proteins: metal ions in life sciences. New York: John Wiley & Sons; 2007.
Celebi F, Guven B. Plasma concentrations of 13,14-dihydro-15-keto PGF2alpha and progesterone during the oviposition cycle of the domestic goose (Anser anser domesticus). Poult Sci. 2001;80(2):225–7.
Sechman A, Staruszczak M, Rząsa J. Comparison of Sex Steroid Concentration in Blood Plasma and Ovarian Follicles of White Leghorn and Greenleg Partridge Laying Hens. Ann Ani Sci. 2011;11(4):507–20.
Bedrak E, Harvey S, Chadwick A. Concentrations of pituitary, gonadal and adrenal hormones in serum of laying and broody white rock hens (Gallus domesticus). J Endocrinol. 1981;89(2):197–204.
We are grateful to Jiangsu Lihua Animal Husbandry CO. LTD for offering the animal sample. This work received financial support from the earmarked fund for Modern Agro-industry Technology Research System (CARS-43-3), the Agricultural Science and Technology Independent Innovation Fund of Jiangsu Province (CX(13)2034) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (2011–137).
The authors declare that they have no competing interests.
Conceived and designed the experiments: GC and QX. Performed the experiments:. QX, YS and YC. Analyzed the data: QX, YS and GC. Contributed reagents/materials/analysis tools: RL, YZ, YL, ZH and WZ. Wrote the paper: QX and GC. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Xu, Q., Song, Y., Chen, Y. et al. Molecular cloning and expression patterns of the cholesterol side chain cleavage enzyme (CYP11A1) gene during the reproductive cycle in goose (Anas cygnoides). J Animal Sci Biotechnol 6, 54 (2015). https://doi.org/10.1186/s40104-015-0053-9
- Gene expression