Enzymes, vectors, media, and strains
T4 DNA ligase, KOD-Plus-Neo polymerase, EcoRI, NotI, and SacI were purchased from Bao Bioengineering, Co., Ltd. (Dalian, China). PMD18-T vector, pPIC9K vector, and E. coli DH5α were purchased from Shanghai Sangon Biotech Co., Ltd. (Shanghai, China). Minimal dextrose (MD), buffered glycerol-complex (BMGY), and buffered methanol-complex (BMMY) medium were prepared according to the recipe provided in the Invitrogen Pichia pastoris expression kit instruction manual.
Animals and salivary tissue sample collection
All leeches were obtained from an adult H. nipponia colony grown in a medical leech breeding base of the Chongqing Academy of Chinese Materia Medica (Chongqing, China). Fifty healthy leeches were maintained in an aquaculture net cage filled with 15 L dechlorinated tap water at 20–22 °C with a 12h/12h day/night cycle prior to dissection. Every 5 days, half of the water was replaced with fresh water. Salivary tissue masses lying posterior to the three muscular jaws were removed aseptically using a sterilized dissecting tool; subsequently, they were rinsed in 0.5% bleach for 1 min followed by rinsing in deionized water for 1 min48-50. The tissues were then stored in RNAlater (Qiagen, Hilden, Germany) according to the manufacturer’s specifications at –80 °C.
RNA extraction and cDNA preparation
Total RNA was extracted from the aforementioned salivary tissues of H. nipponia using Trizol reagent (Tiangen, Beijing, China) according to the manufacturer instructions. RNA quality was assessed by electrophoresis on a 1.0% agarose gel, and RNA concentration was determined using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). First-strand cDNA was synthesized using the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, USA).
Molecular cloning of the hirudin gene of H. nipponia
The transcript (c16237_g1) that we previously isolated from the salivary gland transcriptome databases of H. nipponia27 was screened for further analysis. Since this transcript showed the highest homology with the Hirudin variant HV1(VV) (ALA22934.1), with an e-value of 3e-13, we designated this as a Hirudin transcript. Amplification primers were designed using Primer Premier 5 software to clone the coding sequence of Hirudin and synthesized by Shanghai Sangon Co., Ltd. (see Table 1). The cycling parameters were 94 °C for 2 min; followed by 35 cycles of 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s; and final extension at 72 °C for 10 min. To obtain full-length cDNA of the Hirudin, rapid amplification of 5′ and 3′ cDNA ends (RACE) was performed using a FirstChoice RLM-RACE Kit (Invitrogen, Waltham, MA, USA) according to the manufacturer’s manual. Pairs of gene-specific primers were designed for 5′ and 3′ RACE polymerase chain reaction (PCR) based on the acquired coding DNA sequence (Table 1).
The PCR products were analyzed on gels, and the bands of interest were excised and purified with MiniBEST Agarose Gel DNA Extraction Kit (Takara Bio Inc., Shiga, Japan). The fragments were subcloned into a pMD-18T vector and sequenced by Shanghai Sangon Co., Ltd.
Sequence analysis and multiple sequence alignment of hirudin
Based on the cloned sequence, the full-length aa sequence of hirudin protein was predicted using ORFfinder (http://www.ncbi.nlm.nih.gov/orffinder/). The molecular weight and theoretical isoelectric point (pI) were predicted using ProtParam (https://www.expasy.org/tools/protparam.html). Protein signal peptides were predicted using the SignalP 6.0 server (https://services.healthtech.dtu.dk/service.php?SignalP). Prediction of the physicochemical properties of the hirudin protein after removal of signal peptides was performed with the ExPASy ProtParam tool (https://web.expasy.org/protparam/).
A protein BLAST (http://www.ncbi.nlm.nih.gov/BLAST) search was performed based on the hirudin aa sequence. The other 17 sequences of hirudin from different leech species were obtained from the National Center for Biotechnology Information (NCBI) website and their relevant information is presented in Table 2. Multiple sequence alignment was accomplished with Clustal W ver. 2.0.10 and GeneDoc ver. 2.7.0.
Construction of pPIC9K-Hirudin recombinant plasmid
The coding sequence of Hirudin was amplified using appropriate primer pairs (Hir-p-F/Hir-p-R; Table 1) and KOD-Plus-Neo polymerase. Thus, appropriate restriction sites for the enzymes EcoRI and NotI as well as protecting bases were added to both ends of the gene. The pPIC9K vector and the PCR product were subjected to 2 h of double enzyme digestion using the same EcoRI and NotI enzymes at 37 °C and ligated with T4 DNA ligase. The recombinant plasmid pPIC9K-Hirudin was then transformed into competent E. coli DH5α cells for amplification. The positive transformants were screened and verified by colony PCR and double enzyme digestion, and the PCR product was analyzed on a 1.5% agarose gel.
Electro-transformation of Pichia pastoris GS115
P. pastoris GS115 competent cells were prepared according to the manual of Invitrogen Corporation. The recombinant expression plasmid was linearized by the SacI restriction enzyme and transformed electronically into P. pastoris GS115. Pichia pastoris GS115 cells harboring pPIC9K-Hirudin were inoculated into 200 μL of MD medium and cultivated for 3 days at 30 °C to screen positive colonies. Single yeast colonies were retrieved from the MD plate using sterile toothpicks. Colony identification was confirmed by PCR using the primers 5′ AOX1 and 3′ AOX1 (Table 1). PCR products were analyzed on a 1.5% agarose gel, and the target bands were excised and purified with MiniBEST Agarose Gel DNA Extraction Kit (Takara). P. pastoris GS115 transformants harboring pPIC9K-Hirudin were designated as GS115/pPIC9K-Hirudin. P. pastoris cells transformed with the empty pPIC9K plasmid were applied as the control and named GS115/pPIC9K.
Induction and expression of recombinant yeast transformants
The yeast clones expressing pPIC9K-Hrudin or pPIC9K were respectively inoculated into a 250-mL bottle containing 50 mL BMGY medium overnight after identified by PCR. The yeast was cultured at 28 ℃ with shaking at 250 rpm until the optical density at 600 nm (A600) reached 2–4. The yeast was then collected and induced in 80 mL BMMY medium until the A600 reached 1.5 and cultured continuously at 28 ℃ for 3 days. During induction, samples were collected every 24 h and methanol was added to a final concentration of 0.05% (v/v). Samples were centrifuged at 12,000 ×g for 2 min at 4 ℃. The supernatant and precipitate were both collected and analyzed by 15% SDS-PAGE as described previously51.
Western blot analysis
Supernatants and precipitates from GS115/pPIC9K and GS115/pPIC9K-Hirudin were respectively loaded on a 150 g/L SDS-PAGE gel for protein separation and transferred to a polyvinylidene difluoride membrane by electrotransfer (54 mA, 80 min). After blocking with 1.5% BSA in blocking buffer at 4 ℃ overnight, the membrane was incubated with 1000× dilutions of Penta-His antibody (Qiagen) at room temperature for 1 h, washed three times with tris-buffered saline, and mixed with 1000× dilutions of rabbit anti-mouse IgG antibody coupled with horseradish peroxidase (Takara) at room temperature for 1 h. The bands were visualized by TrueBlue Peroxidase Substrate (Takara) for 1 min and then with WesternBrightTM ECL (Takara) for 5 min.
Antithrombin activity analysis
The antithrombin activity of hirudin in the supernatants of P. pastoris GS115/pPIC9K-Hirudin after 72 h induction was evaluated according to the Pharmacopoeia of the People’s Republic of China (2020 edition)7 and a previously reported method52. Briefly, antithrombin activity was quantitatively measured by titrating a solution of thrombin and expressed in antithrombin units (ATU): one unit is defined as the neutralization of one NIH unit of thrombin (Sigma, Germany) at 37 ℃. For the activity assay, 200 μL of 0.05% bovine fibrinogen solution prepared with 50 mM Tris-HCl buffer (pH 7.4) was thoroughly mixed with 100 μL of the sample, and then 5 μL of 40 NIH/mL thrombin solution (0.2 NIH unit) was added progressively and mixed gently. The above reaction mixture was incubated at 37 ℃ ± 0.5 ℃ for 1 min. The endpoint of the titration was considered to occur when a fibrin clot formed within 1 min. Otherwise, another 5 μL of thrombin solution was added continuously until a fibrin clot was observed.
Ethical approval
We declare that the experiments described in this paper comply with the current laws in China. This article does not contain any studies with human participants performed by any of the authors.
Data availablity
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.