Cloning and expression of Iranian Turkmen-thoroughbred horse follicle stimulating hormone in Pichia pastoris

Document Type : Research Paper


1 Department of Genetics, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran Department of Animal Breeding, Genetics and Physiology, Gorgan University of Agricultural Sciences and Natural Resources, Gorgan, Iran

2 Department of Genetics, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran

3 Department of Animal Breeding, Genetics and Physiology, Gorgan University of Agricultural Sciences and Natural Resources, Gorgan, Iran

4 Department of Genetics, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran.

5 Institute of Bioscience, Universiti Putra Malaysia, Serdang, Selangor, Malaysia

6 Department of Genetics, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran,



Background: Follicle stimulating hormone (FSH) plays an essential role in reproductive physiology and follicular development. Objective: A new variant of the equine fsh (efsh) gene was cloned, sequenced, and expressed in Pichia pastoris (P. pastoris) GS115 yeast expression system. Materials and Methods: The full-length cDNAs of the efshα and efshβ chains were amplified by reverse transcription polymerase chain reaction (RT-PCR) using the total RNA isolated from an Iranian Turkmen-thoroughbred horse's anterior pituitary gland. The amplified efsh chains were cloned into the pPIC9 vector and transferred into P. pastoris. The secretion of recombined eFSH using P. pastoris expression system was confirmed by Western blotting and immunoprecipitation (IP) methods. Results: The DNA sequence of the efshβ chain accession number JX861871, predicted two putative differential nucleotide arrays, both of which are located in the 3′ UTR. Western blotting showed a molecular mass of 13 and 18 kDa for eFSHα and eFSHβ subunits, respectively. The expression of desired protein was confirmed by protein G immunoprecipitation kit. Conclusions: eFSH successfully expressed in P. pastoris. These findings lay a foundation to improve ovulation and embryo recovery rates as well as the efficiency of total embryo-transfer process in mares.


Main Subjects

1. Background
    Superovulation in horses is much less successful than in other domestic species. In addition to superovulation rate, the number of embryos are typically negligible (1). Hence, there is a great demand to develop products that induce multiple ovulations in mares due to changes in the governance of the horse reproduction industry (2). Follicle stimulating hormone (FSH), a member of the glycoprotein hormone family, plays an essential role in this process (3). Generally, the anterior pituitary gland is the tissue that is responsible for releasing this hormone (4-8), and granulosa cells in the ovary are the target cells for the numerous actions of FSH (6, 8-11). Similar to other members of the glycoprotein hormone family, including luteinizing hormone (LH), chorionic gonadotropin hormone (CGH), and thyroid stimulating hormone (TSH), FSH consists of two heterodimer non-covalently linked subunits: a common alpha subunit and a hormone-specific beta subunit (6, 12). An increase in the exogenous FSH concentration in the circulatory system during follicular development may leads to increase in the number of ovulatory follicles, which results in an increase in the number of mature follicles, the ovulatory rate, and the rate of embryo transfer (3). Currently, this method is applicable for both human infertility and animal breeding programs(13). Some gonadotropin sources with FSH-like activity, such as pituitary extraction from domestic animals, pregnant mare serum gonadotropin (PMSG), and FSH extraction from the urine of postmenopausal women, have been used for the same purpose. Despite the success of some of these treatments in many domestic species, these sources are less effective in mares (14-18). It has been reported that recombinant eFSH (reFSH) shows a much more specific effect in increasing the number of total preovulatory follicles than other therapeutic methods (3). Many of the recombinant proteins that have been successfully produced using the yeast expression system have been used in basic laboratory research and the therapeutic industry (19). The methylotrophic P. pastoris yeast is a suitable host for the production of many recombinant proteins because of its easy genetic manipulation, expression of intracellular or extracellular foreign protein (20), and capability to perform eukaryotic protein modifications, such as glycosylation, disulfide bond formation, and proteolytic processing (20-22).

2. Objectives
    The purpose of the present research was to express the recombinant horse FSH using the P. pastoris expression system and to detect this protein in the secretion media. This study describes the first production of the horse FSH in the yeast expression system using the Iranian native mare fsh gene.

3. Materials and Methods
    The chemical and biological laboratory reagents and instruments that were used in this study include NucleoSpin RNA kit (Macherey Nagel, Düren, Germany); Revert AidTM H minus first-strand cDNA synthesis kit (Fermentas, Massachusetts, USA); pTZ57R/T vector (Fermentas, Massachusetts, USA); DNA sequencing analysis (MWG, Ebersberg, Germany); Pichia expression kit, (Invitrogen, Massachusetts, USA); P. pastoris cell culture media (Invitrogen, Massachusetts, USA); Gene-Pulser (Bio-Rad, callifornia, USA); sorbitol (Invitrogen, Massachusetts, USA); Amplicon Ultra 3000 MWCO (Millipore Bedford, Massachusetts, USA); Primary antibodies against the eFSHa and eFSHb subunits  (Santacruz, Texas, USA); horse radish peroxidase (HRP)-conjugated secondary antibodies against eFSHa and eFSHb subunits (Sigma Aldrich, Missouri, USA); nitrocellulose membranes (Amersham, london, UK); enhanced chemiluminescence (ECL) reagent (Amersham ECL, London, UK); Aautoradiography GBX developer (Sigma Aldrich, Missouri, USA); Amersham Hyperfilm ECL instrument (Amersham, london, UK); and protein G immunoprecipitation kit (Sigma Aldrich, Missouri, USA).

3.1. Strains and Media
    The E. coli DH5α strain and P. pastoris yeast in a Pichia expression kit were used as the competent cells. Luria Bertani (LB) medium (1% (w/v) tryptone, 0.5% (w/v) yeast extract, and 1% (w/v) NaCl, pH 7.0) was used as the bacterial growth culture at 37°C (25). The P. pastoris cell culture media  was composed of yeast extract peptone dextrose (YEPD) medium (1% (w/v) yeast extract, 2% (w/v) peptone, 2% (w/v) dextrose, and 0.1% (w/v) ampicillin), minimal dextrose (MD) medium (1.34% (w/v) yeast nitrogen base with ammonium sulfate without amino acids (YNB), 4×10-5% (w/v) biotin, 2% (w/v) dextrose, 0.05% (w/v) ampicillin, and 1.5% (w/v) agarose), buffered glycerol-complex (BMGY) medium (1% (w/v) yeast extract, 2% (w/v) Bacto-peptone, 1.34% (w/v) YNB, 1% (v/v) glycerol, 4×10-5% (w/v) biotin, and 100 mM potassium phosphate buffer pH 6.0) and buffered methanol-complex (BMMY) medium, supplemented with 0.5% (v/v) methanol.

3.2. Cloning of the efshα and efshb cDNAs into Expression Vector
    The total mRNA from the Iranian Turkmen-thoroughbred mare’s anterior pituitary gland was freshly extracted using a NucleoSpin RNA kit. The cDNA synthesis was then performed with 1 mg of total mRNA using a Revert AidTM H minus first-strand cDNA synthesis kit and an oligo (dt) primer that, followed by RT-PCR.
    The primers carrying the AvrII site (Table 1) were used to amplify the efsha and efshb cDNA sequences by polymerase chain reaction (PCR).
     The PCR mixture for each subunit consisted of a total volume of 25 mL that included 2 units of Taq DNA polymerase, 2 mL of RT-cDNA, 2.5 mL of 10× PCR buffer, 1 mL of 10 mM dNTP, 1.5 mL of 50 mM MgCl2, and 0.5 mL of 10 pmoles of each forward and reverse primer. The PCR thermocycling steps for efshα and efshb genes were as follow in Table 2.
     The amplified PCR products were separately cloned through the T/A cloning method using the pTZ57R/T vector and frozen DH5a competent cells. The recombinant bacterial colonies were identified by colony PCR, plasmid DNA digestion using AvrII, and DNA sequencing analysis. Ultimately, the gel-purified products were subcloned by ligation into the AvrII site of the linearized pPIC9 plasmid and then transformed into the DH5a competent cells through the standard protocol (23).

3.3. Recombinant P. pastoris Transformants
    As shown in Figure 1, the recombinant pPIC9-eFSHα and pPIC9-eFSHb vectors were linearized with SalI and SacI, respectively (Pichia expression kit). These recombinant vectors were co-transformed (5 mg) into the host P. pastoris GS115 His+ Mut+ cells (containing aox1 gene) by co-electroporation in a Gene-Pulser with a voltage of 1900 V, capacitance of 25 mF, and resistance of 400 ohms (24).
    The yeast transformants were resuspended in 1 M sorbitol, and the cell mixture was then spread and cultured on selective MD medium plates (histidine-restricted media) at 29.5ºC for four days. The single colonies for His+ transformants were inoculated into the YEPD liquid medium and incubated for 48 h at 29.5ºC with 290-rpm agitation. The presence of both efsha and efshb genes in P. pastoris genome was confirmed by PCR using forward and reverse aox1 primers. The pPIC9 vector without any insert and the recombinant pPIC9-eFSHa and pPIC9-eFSHb vectors were used as control groups (Pichia expression kit).

3.4. Secretion of the Recombinant eFSH
    The positive single-recombinant colonies carrying both recombinant pPIC9-eFSHa and pPIC9-eFSHb vectors were selected and cultured in BMGY medium (25 mL) on a rotary incubator for 18 h at 29.5°C with 290-rpm agitation (25) to reach an OD600 value of approximately 3.0 (24). The yeast cells were harvested by centrifugation at 3000 ×g for 5 min and subsequently resuspended and cultivated in BMMY medium (100 mL). The incubation was continued at 29.5ºC for five days. Methanol was added every 24 h to a final concentration of 0.5% to maintain the induction phase. The samples from five consecutive days (days 0 to 5) were collected at 0, 6, 24, 30, 48, 54, 72, 78, 96, 102, and 120 h after induction. The supernatant of the suspension was separated by centrifugation at 13000 ×g for 3 min and stored at -70ºC for use in the biochemical confirmatory tests. The negative controls, which included GS115 (His+ Muts) cells expressing extracellular albumin, GS115 (His+ Mut+) cells expressing intracellular b-galactosidase, and GS115 (His+ Mut+) cells transformed with non-recombinant pPIC9, were cultured and collected in the same manner (Pichia expression kit).

3.5. Western Analysis
    The secretion media from the recombinant P. pastoris GS115 strain (His+ Mut+) and the control groups were concentrated five-fold using an Amplicon Ultra 3000 MWCO at 13000 ×g. The concentrated samples were electrophoresed in 15% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) apparatuses through the Laemmli method and blotted 16 h overnight onto nitrocellulose membranes. The membranes were incubated in the blocking solution   (1 M TBS-Tween-20 plus 3% non-fat milk) for 2 h with shaking at 22ºC. The primary antibodies against the eFSHa (Sc-292422) and eFSHb (Sc-18222) subunits were diluted at ratios of 1:750 and 1:800, respectively. Each membrane was incubated with a specific primary antibody for 2 h at 22ºC. The membranes were then washed three times with 1 M TBS-Tween-20 solution for 15 min. The corresponding horse radish peroxidase (HRP)-conjugated secondary antibodies against eFSHa (A5420) and eFSHb (A0545) subunits were diluted at ratios of 1:125000 and 1:500000, respectively. Each membrane was incubated with the corresponding specific diluted secondary antibody for 1.5 h at 22ºC. The membranes were washed five times for 15 min and treated with an enhanced chemiluminescence (ECL) reagent to reveal the immunoreactive bands through an equilibration technique using. The signals were detected on an Amersham Hyperfilm ECL instrument.
3.6. Immunoprecipitation
    A protein G immunoprecipitation kit was used to detect the expressed recombinant eFSH in the P. pastoris BMMY medium. This enrichment was accomplished by binding the recombinant eFSH with a primary specific b-subunit antibody followed by the precipitation of the immune complexes with protein G. The concentrated sample obtained 54 h after induction was chosen to confirm the previous experiments. First, 600 mL of the sample was mixed with 0.5 to 1 mg of the primary antibody, and the mixture was shaken for 14 h at 4ºC. After the sample was transferred to a filter column, 30 mL of G-protein was added to the column, and the mixture was shaken for 14 h at 4ºC. Subsequently, the obtained purified recombinant eFSHb sample and the control sample were analyzed using the Western blotting.

4. Results

4.1. Nucleotide Sequence Analysis of efsh Subunits
      After the amplification of the efsha and efshb genes from cDNA (Figure 2), the sequence analysis of these amplified gene fragments showed that the efsha gene had the expected size of 387 base pairs, which confirmed the information included in the GenBank. However, the efshb gene was a fragment with 465 bp. Our sequence analysis was revealed that eFSHb reported here was a bit different (accession no. JX861871) to what was in GenBank at its 3¢-UTR (Figure 3) with 3 more nucleotides within the 413-to-418 (GCC ® GCTGCC) and having a nucleotide substitution at 443 base pair region (TAT instead of TAC).

4.2. Genomic Integration of the pPIC9-eFSHa and pPIC9-e FSHb Vectors
    The recombinant pPIC9-eFSHα and pPIC9-eFSHb vectors were co-inserted randomly into the genome of P. pastoris cells (Pichia expression kit). The agarose gel electrophoresis showed that which of the P. pastoris GS115 (His+ Mut+) transformants containing both recombinant pPIC9-eFSHα and pPIC9-eFSHb vectors (Figure 4). Based on the standard protocol (Pichia expression kit), the expected size of the P. pastoris genomic aox1 gene and the integrated efsha and efshb genes were approximately 2.2 kb, 891 bp and 969 bp respectively. These colonies were then selected and cultured in specific media for expression of the recombinant protein.

4.3. Western Analysis
    The presence of reFSH in each collected medium sample was confirmed by Western analysis. The Western blotting with specific polyclonal antibodies for each subunit revealed two bands of approximately 13 and 23 kDa, for eFSHa subunit, and two bands of approximately 18 and 29 kDa, corresponding to the eFSHb subunit (Figure 5). The bands began appearing 24 h after methanol induction, and strong bands were observed 54 h after methanol induction.

4.4. Immunoprecipitation    
    The IP analysis with a rabbit polyclonal antibody against the FSHb subunit (a common antibody that is used for both human and equine species) showed a similar pattern between the bands that appeared in the concentrated sample obtained 54 h after induction and the pure human FSH sample (positive control). As shown in Figure 6, two similar bands were found in each lane at the same position: the first band is approximately 29 kDa, and the second band is 50 kDa. We predicted that the smaller band corresponds to the FSHβ subunit as it could not be detected in the negative control, and the larger band can be corresponded either to the G-protein or our product.

5. Discussion
    The eFSH protein has been used to improve ovulation and embryo recovery rates, which can most likely increase the efficiency of the total embryo-transfer process (26). Jennings et al. (2009) reported that the stimulation of follicles with reFSH may be more potent than the stimulation with purified eFSH, which would lead to an increase in the follicular growth rate in mares (3). Thus, it is obvious that the production of reFSH has a practical application in the equine reproduction industry. Therefore, the current study provides the first report of the utilization of the yeast P. pastoris expression system for the secretion of a new variant of the Iranian Turkmen-thoroughbred horse fsh gene. Saneyoshi et al. (2001) reported a 462 bp sequence for efshb cDNA (27). The amino acid sequence that was deduced from their study is similar to the data that were reported by Fujiki et al. (1978) (3). However, our findings showed a 465 bp sequence for the efshb with some significant variations in 3΄-UTR. The heterologous expression of proteins with disulfide bonds often fails in bacteria (28, 29). As a eukaryotic expression system, P. pastoris was chosen because of its cost-effectiveness, potency for appropriate post translational modifications, secretion of foreign proteins into the culture medium (25), and a respiratory rather than a fermentative mode of growth that lowers the production of toxic waste materials in a high-cell-density environment (30). In fact, the expression of recombinant mammalian GTHs, such as porcine FSH (7), ovine FSH (3), and human chorionic gonadotropin (hCG) (31), in P. pastoris has shown equal biological activity compared with the native GTHs. Thus, P. pastoris was used for the expression of reFSH. We co-integrated the efshα and efshb genes into the P. pastoris genome by separated signal peptide sequences. This secretion pattern is consistent with the native form, in which the two subunits are co-expressed using two independent expression constructs. Yu et al. (2010) reported that the abovementioned form of protein secretion is much more efficient than the use of a common signal peptide leader (25). The integration of a foreign gene into the P. pastoris genome is a preferential method (3) that causes more stable transformants (32) and thus leads to a high expression level of the heterologous gene (3). The recombinant gene constructs were co-electroporated into the P. pastoris genome based on different sites to increase the possibility of the co-insertion of both gene constructs into the genome. The analysis of the Western blotting results revealed that each subunit of the reFSH molecule was successfully synthesized and secreted into the P. pastoris culture medium; however, the results showed two specific bands for each subunit. Based on reports by Kasuto and Levavi-Sivan (2005), it is likely that the larger protein bands (23 and 29 kDa bands in the Figure 5) are highly glycosylated forms of the subunits or multimers (33). The existence of some bands greater than 50 kDa that immunoreacted with the eFSHα and even the eFSHb antisera can be considered nonspecific bands that originated from the yeast (34). Although successful expression of recombinant FSH was reported in CHO by Fachal et al. (2010), the purification of this intracellular protein is often tedious and time consuming. In the current study, production of the recombinant eFSH was conducted in P. pastoris using extracellular secretion system. Application of this system resulted in secretion of the protein into the media which led to simple afterward purification. Post-translational modification is a concern in production of recombinant proteins in P. pastoris. However, this concern could be solved through application of Glyco-Switc strain with the ability of post translational modifications same as mammalian cells (35-37). Meanwhile, further studies are required to compare the functionality of the P. pastoris-derived eFSH in mares or other species.

    This study was substantially funded by the Royan Institute, Tehran, Iran. We express our appreciation to Gorgan University of Agricultural Sciences and Natural Resources and to the Royan proteomics lab for facilitating the research phases. In addition, the authors thank Zahra Mansouri and Najmeh Sadat Masoudi for their kind assistance.

1.    Briant C, Toutain PL, Ottogalli M, Magallon T, Guillaume D. Kinetic studies and production rate of equine (e) FSH in ovariectomized pony mares. Application to the determination of a dosage regimen for eFSH in a superovulation treatment. J Endocrinol. 2004;182(1):43-54. DOI: 10.1677/joe.0.1820043
2.    Fachal MV, Furlan M, Clark R, Card CE, Chedrese PJ. Synthesis and characterization of biologically active recombinant elk and horse FSH. J Anim Reprod Sci. 2010;117(3):331-40. DOI:
3.    Jennings M, Boime I, Daphna-Iken D, Jablonka-Shariff A, Conley A, Colgin M,  Jennings MW, Boime I, Daphna-Iken D, Jablonka-Shariff A, Conley AJ, Colgin M, Bidstrup LA, Meyers-Brown GA, Famula TR, Roser JF. The efficacy of recombinant equine follicle stimulating hormone (reFSH) to promote follicular growth in mares using a follicular suppression model. J Anim Reprod Sci. 2009;116(3):291-307. DOI:
4.    Dias JA, Cohen BD, Lindau-Shepard B, Nechamen CA, Peterson AJ, Schmidt A. Molecular, structural, and cellular biology of follitropin and follitropin receptor. J Vitam Horm. 2002;64:249-322. DOI: 10.1016/S0083-6729(02)64008-7
5.    Allen W, Moor R. The origin of the equine endometrial cups. J Reproduction Fertility. 1972;29(2):313-6. 5.1591907. DOI: 10.1530/jrf.0.0290313
6.    Pierce JG, Parsons TF. Glycoprotein hormones: structure and function. J Annu Rev Biochem. 1981;50(1):465-495. DOI: 10.1146/
7.    Richard F, Robert P, Remy J-J, Martinat N, Bidart J-M, Salesse R, Richard F, Robert P, Remy J-J, Martinat N, Bidart J-M, Salesse R, Combarnous Y. High-Level secretion of biologically active recombinant porcine follicle-stimulating hormone by the methylotrophic yeast Pichia pastoris. Biochem Biophys Res Commun. 1998;245(3):847-852. DOI: 10.1006/ bbrc.1998.8532
8.    Simoni M, Gromoll Jr, Nieschlag E. The follicle-stimulating hormone receptor: biochemistry, molecular biology, physiology, and pathophysiology 1. J Endocrine Rev. 1997;18(6):739-73. DOI:
9.    Licht P, Papkoff H, Farmer S, Muller CH, Tsui HW, Crews D. Evolution of gonadotropin structure and function. Recent Prog Horm Res. 1975;33:169-248. DOI: 10.1016/B978-0-12-571133-3.50012-X
10.    Mapletoft RJ, Steward KB, Adams GP. Recent advances in the superovulation in cattle. J Reprod Nutr Dev. 2002;42(6):601-611. DOI: 10.1051/rnd:2002046
11.    Squires EL. Superovulation in Mares. Veterinary Clinics: J Equine Pract. 2006;22(3):819-830. DOI: 10.1016/j.cveq.2006. 07.005
12.    Gharib SD, Wierman ME, Shupnik MA, Chin WW. Molecular biology of the pituitary gonadotropins. J Endocrine Rev. 1990;11(1):177-199. DOI: 10.1210/edrv-11-1-177
13.    Ben-Chetrit A, Gotlieb L, Wong P, Casper R. Ovarian response to recombinant human follicle-stimulating hormone in luteinizing hormone-depleted women: examination of the two cell, two gonadotropin theory. J Fertil Steril. 1996;65(4):711-717. DOI:
14.    Hofferer S, Lecompte F, Magallon T, Palmer E, Combarnous Y. Induction of ovulation and superovulation in mares using equine LH and FSH separated by hydrophobic interaction chromatography. J F Reprod Fertil. 1993;98(2):597-602. DOI: 10.1530/jrf.0.0980597
15.    McCue P, LeBlanc M, Squires E. eFSH in clinical equine practice. J Theriogenol. 2007;68(3):429-33. DOI:
16.    Aucouturier P, Carp RI, Carnaud C, Wisniewski T. Prion diseases and the immune system. J Clin Immunol. 2000;96(2):79-85. DOI: 10.1006/clim.2000.4875
17.    Tharasanit T, Colenbrander B, Bevers M, Stout T. Effects of recombinant human follicle stimulating hormone on follicle development and ovulation in the mare. J Theriogenol. 2006;65(6):1071-1081. DOI:
18.    Squires E, McCue P. Superovulation in mares. J Anim Reprod Sci. 2007;99(1):1-8. DOI:
19.    Demain AL, Vaishnav P. Production of recombinant proteins by microbes and higher organisms. J Biotechnol Adv. 2009;27(3):297-306. DOI: 10.1016/j.biotechadv.2009.01.008
20.    Cregg JM, Cereghino JL, Shi J, Higgins DR. Recombinant protein expression in Pichia pastoris. J Mol Biotechnol. 2000;16(1):23-52.
21.    Briand L, Eloit C, Nespoulous C, Bézirard V, Huet J-C, Henry C, Briand L, Eloit C, Nespoulous C, Bézirard V, Huet J-C, Henry C, Blon F, Trotier D, Pernollet J-C. Evidence of an odorant-binding protein in the human olfactory mucus: location, structural characterization, and odorant-binding properties. J Biochem. 2002;41(23):7241-7252. DOI: 10.1021/bi 015916c
22.    Briand L, Nespoulous C, Huet JC, Takahashi M, Pernollet JC. Ligand binding and physico-chemical properties of ASP2, a recombinant odorant-binding protein from honeybee (Apis mellifera L.). Eur J Biochem. 2001;268(3):752-760. DOI: 10.1046/j.1432-1327.2001.01927.x.
23.    Sambrook J, Russell DW. Molecular cloning. A laboratory manual. Third. Cold pring Harbor Laboratory Press, New York. 2001.
24.    Du CH, Xiao AF. Secretory expression and purification of the recombinant duck interleukin-2 in Pichia pastoris. J Microbiol Biotechnol. 2011;21(12):1264-1269.
25.    Yu X, Lin S-W, Kobayashi M, Ge W. Expression of recombinant zebrafish follicle-Stimulating Hormone (FSH) in Methylotropic Yeast Pichia pastoris. Fish Physiol Biochem. 2010;36(2):273-281.
26.    Logan NL, McCue PM, Squires EL, editors. How to use equine FSH efficiently to enhance embryo recovery. Proceedings of the annual convention; 2006.
27.    Saneyoshi T, Min K-S, Ma XJ, Nambo Y, Hiyama T, Tanaka S, Saneyoshi T, Min K-S, Ma XJ, Nambo Y, Hiyama T, Tanaka S, Shiota K. Equine follicle-stimulating hormone: molecular cloning of b subunit and biological role of the asparagine-linked oligosaccharide at asparagine56 of a subunit. J Biol Reprod. 2001;65(6):1686-1690. DOI: 10.1095/biolreprod65.6.1686
28.    Almeida MS, Cabral KS, de Medeiros LN, Valente AP, Almeida FC, Kurtenbach E. cDNA cloning and heterologous expression of functional cysteine-rich antifungal protein Psd1 in the yeast Pichia pastoris. J Arch Biochem Biophys. 2001;395(2):199-207. DOI: 10.1006/abbi.2001.2564
29.    Li Z, Moy A, Sohal K, Dam C, Kuo P, Whittaker J, Li Z, Moy A, Sohal K, Dam C, Kuo P, Whittaker J, Whittaker M, Düzgünes N, Konopka K, Franz A. Expression and characterization of recombinant human secretory leukocyte protease inhibitor (SLPI) protein from Pichia pastoris. J Protein Expr Purif. 2009;67(2):175-181. DOI: 10.1016/j.pep.2009.06.001
30.    Cereghino GPL, Cereghino JL, Ilgen C, Cregg JM. Production of recombinant proteins in fermenter cultures of the yeast Pichia pastoris. J Curr Opin Biotechnol. 2002;13(4):329-332. DOI: 10.1016/S0958-1669(02)00330-0
31.    Gupta CS, Dighe R. Hyperexpression of biologically active human chorionic gonadotropin using the methylotropic yeast, Pichia pastoris. J Mol Endocrinol. 1999;22(3):273-283. DOI:  10.1677/jme.0.0220273
32.    Daly R, Hearn MT. Expression of heterologous proteins in Pichia pastoris: a useful experimental tool in protein engineering and production. J Mol Recognit. 2005;18(2):119-138. DOI: 10.1002/jmr.687
33.    Kasuto H, Levavi-Sivan B. Production of biologically active tethered tilapia LHβα by the methylotrophic yeast Pichia pastoris. J Gen Comp Endocrinol. 2005;140(3):222-232. DOI: 10.1016/j.ygcen.2004.10.016
34.    Kamei H, Ohira T, Yoshiura Y, Uchida N, Nagasawa H, Aida K. Expression of a biologically active recombinant follicle stimulating hormone of Japanese eel Anguilla japonica using methylotropic yeast, Pichia pastoris. J Gen Comp Endocrinol. 2003;134(3):244-254. DOI: 10.1016/S0016-6480(03)00259-4
35.    Hamilton SR, Bobrowicz P, Bobrowicz B, Davidson RC, Li H, Mitchell T, Nett JH, Rausch S, Stadheim TA, Wischnewski H. Production of complex human glycoproteins in yeast. Science 2003;301(5637):1244-1246.
36.    Houdebine L-M. Production of pharmaceutical proteins by transgenic animals. J Comp Immunol Microbiol Infect Dis. 2009;32(2):107-121. DOI: 10.1016/j.cimid.2007.11.005
37.    De Schutter K, Lin Y-C, Tiels P, Van Hecke A, Glinka S, Weber-Lehmann J, De Schutter K, Lin Y-C, Tiels P, Van Hecke A, Glinka S, Weber-Lehmann J, Rouzé P, Van de Peer Y, Callewaert N. Genome sequence of the recombinant protein production host Pichia pastoris. J Nat Biotechnol. 2009;27(6):561-566. DOI: 10.1038/nbt.1544