Recombinant Dermatophagoides farinae Type 1 Allergen Protein and Its Preparation Method and Application

Abstract

Provided are an optimized proDer f1 gene, a proDer f1 protein encoded thereby, a vector comprising said gene, and a Pichia pastoris strain. Also provided are an expression method and a purification method of the proDer f1 protein.

Claims

1. A DNA sequence encoding proDer f1 protein, having a base sequence as shown in SEQ ID NO: 1.

2. A DNA sequence of claim 1, wherein the DNA sequence is comprised in the vector pAO815, pPIC9, pPIC9K, pPIC3.5, pPIC3.5K, pPICZ A, B, C or pGAPZ A, B, C.

3. A DNA sequence of claim 2, wherein the vector is comprised in the Pichia pastoris strain SMD1168, GS115, KM71, X33 or KM71H.

4. A DNA sequence of claim 3, wherein there is 242 bp interval between the DNA sequence encoding proDer f1 protein and the ATG of AOX1 on Pichia pastoris; and the DNA sequence encoding the proDer f1 protein is preceded by an alpha-factor signal peptide and Kozak sequence GCCACCATGG.

5. Recombinant protein, having an amino acid sequence as shown in SEQ ID NO: 3 or 4.

6. Recombinant protein, of claim 5, wherein the recombinant protein is encoded by the base sequence as shown in SEQ ID NO: 1.

7. Recombinant protein of claim 6, wherein the DNA sequence is comprised in the vector pAO815, pPIC9, pPIC9K, pPIC3.5, pPIC3.5K, pPICZ A, B, C or pGAPZ A, B, C.

8. Recombinant protein of claim 6, wherein the vector is comprised in the Pichia pastoris strain SMD1168, GS115, KM71, X33 or KM71H.

9. Recombinant protein of claim 6, wherein there is 242 bp interval between the DNA sequence encoding proDer f1 protein and the ATG of AOX1 on Pichia pastoris; and the DNA sequence encoding the proDer f1 protein is preceded by an alpha-factor signal peptide and Kozak sequence GCCACCATGG.

10. The use of the recombinant protein of claim 5 in the preparation of a medicament for treating a dust mite allergic disease.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0029] FIG. 1 is a comparison plot of sequences of the recombinant proDer f1 gene before and after optimization.

[0030] The sequence before optimization corresponds to the nucleotide sequence of the natural proDer f1 gene; the sequence after optimization corresponds to the nucleotide sequence of the recombinant proDer f1 gene of the present invention, that is, the codon-optimized sequence.

[0031] FIGS. 2A and 2B show the CAI indices of the proDer f1 gene in the Pichia pastoris expression system before and after optimization.

[0032] FIG. 2A shows that the CAI index of the nucleotide sequence of the natural proDer f1 gene in the Pichia pastoris expression system was calculated by a program to be 0.76. FIG. 2B shows that the CAI index of the optimized proDer f1 codon of the present invention in the Pichia pastoris expression system is calculated by a program to be 0.85.

[0033] FIGS. 3A and 3B are optimal codon frequency distribution region plots of the proDer f1 gene in the Pichia pastoris expression host before and after codon optimization.

[0034] FIG. 3A is an optimal codon frequency distribution region plot of the nucleotide sequence of natural proDer f1 gene in the Pichia pastoris system, and it can be seen from the figure that the occurrence percentage of low-utilization codon in the nucleotide sequence of natural proDer f1 gene is 10%. FIG. 3B shows an optimal codon frequency distribution region plot of the optimized proDer f1 codon of the present invention in the Pichia pastoris system, and the occurrence rate of low-utilization codon in the sequence of optimized proDer f1 codon of the present invention is 0.

[0035] FIGS. 4A and 4B are average GC base content distribution region plots of the proDer f1 gene in the Pichia pastoris expression system before and after codon optimization.

[0036] FIG. 4A shows that the average GC base content of the nucleotide sequence of the natural proDer f1 gene in the Pichia pastoris expression system is 41.85%. FIG. 4B shows that the average GC base content of optimized proDer f1 codon of the present invention in the Pichia pastoris expression system is 42.10%.

[0037] FIG. 5 is an agarose gel electrophoretogram of a PCR product of the codon-optimized proDer f1 gene.

[0038] Lane 1 represents 200 bp DNA ladder; lane 2 represents a PCR product of the recombinant proDer f1 gene containing XhoI and NotI restriction sites at both ends.

[0039] FIG. 6 is a diagram showing a construction process of the expression plasmid pPICZ-proDer f1 for codon-optimized proDer f1.

[0040] FIGS. 7A-7B are diagrams showing the identification of expression of the codon-optimized proDer f1 gene in the host engineering bacteria.

[0041] FIG. 7A is a SDS-PAGE gel electrophoretogram of a supernatant of a solution of the host engineering strain containing the codon-optimized proDer f1 gene, after methanol-induced expression for one week. Lane 1 represents pre-stained protein loading markers in the range of 10-250 KD; and other lanes represent supernatants of cultured solutions of proDer f1 gene-positive monoclonal host engineering strains screened by Zeocin.

[0042] FIG. 7B is a western blot plot of a supernatant of a solution of the host engineering strain containing the codon-optimized proDer f1 gene, after methanol-induced expression for one week. Lane 1 represents pre-stained protein loading markers in the range of 10-250 KD; and lanes 2-10 represent supernatants of proDerf1 monoclonal-induced expression

[0043] FIGS. 8A-8B show a chromatogram of the supernatant of proDer f1 fermentation broth by cation chromatography of the first step, and a gel electrophoretogram.

[0044] FIG. 8A is a chromatogram of the supernatant of proDer f1 fermentation broth by cation chromatography of the first step. FIG. 8B is the identification result of cation chromatography purification of the supernatant of proDer f1 fermentation broth, wherein lane 1 represents 11-100 KD non-pre-stained protein markers, lane 2 represents the supernatant of the proDer f1 fermentation broth before purification, lane 3 represents the flow-through liquid, and lanes 4-8 represent the elution of each tube.

[0045] FIG. 9A is a chromatogram of Der f1 protein by anion chromatography of the second step, and FIG. 9B is a gel electrophoretogram.

[0046] FIG. 9A is a chromatogram of Der f1 protein by anion chromatography in step 2. FIG. 9B is the identification result of anion chromatography in the second step for proDerf1 protein purification, wherein lane 1 represents 11-100 KD non-pre-stained protein markers, lane 2 represents the supernatant before purification of Der f1 protein, lane 3 represents the flow-through liquid, and lane 4 represents elution peaks.

[0047] FIG. 10A shows a chromatogram of Der f1 protein by hydrophobic chromatography of the third step and FIG. 10B is a gel electrophoretogram.

[0048] FIG. 10A is a chromatogram of Der f1 protein by hydrophobic chromatography of the third step. FIG. 10B is a purification result of Der f1 protein by hydrophobic chromatography, wherein lane 1 represents 11-100 KD non-pre-stained protein markers, and lane 2-10 represent respective elution tubes.

[0049] FIG. 11 is a comparison of reactivity to positive serum between recombinant Der f1 and natural Der f1, wherein nDerf1 represents the natural Der f1 protein, rDerf1 represents the recombinant Der f1 protein, and NC represents a PBS solution at pH 7.4.

[0050] FIG. 12 is an agarose gel electrophoretogram of a PCR-amplified GAP gene, wherein lane 1 represents 250 bp DNA ladder and lane 2 represents the GAP gene.

[0051] FIG. 13 is an agarose gel electrophoretogram of positive clone of GAP gene T-vector identified by PCR, wherein lane 1 represents 250 bp DNA ladder, lanes 2-11 represent positive clones obtained by blue-white screening, and lane 12 represents a negative clone obtained by blue-white screening.

[0052] FIG. 14 is an agarose gel electrophoretogram of proDer f1 gene identified by PCR, wherein lane 1 represents 500 bp DNA ladder and lane 2 represents the proDer f1 gene.

[0053] FIG. 15 is an agarose gel electrophoretogram for identifying T-vector positive clone of proDer f1 gene through PCR, wherein lane 1 and 2 represent negative clones obtained by blue-white screening, lanes 3 represents 500 bp DNA ladder, lanes 4-13 represent positive clones obtained by blue-white screening, and lane 14 represents a positive control (proDer f1 gene), in which line 4, 6, 7 are positive clones containing proDer f1 gene, other lines are false positive clones.

[0054] FIGS. 16A and 16B show amplification curves of a standard plasmid.

[0055] FIG. 16A shows amplification curves of the standard plasmid T-GAP, and FIG. 16B shows amplification curves of the standard plasmid T-proDer f1.

[0056] FIGS. 17A and 17B show melting curves of a standard plasmid.

[0057] FIG. 17A shows melting curves of the standard plasmid T-GAP, and FIG. 17B shows melting curves of the standard plasmid T-proDer f1.

[0058] FIGS. 18A and 18B show a standard curve of a standard plasmid.

[0059] FIG. 18A shows a standard curve of the standard plasmid T-GAP, and FIG. 18B shows a standard curve of the standard plasmid T-proDer f1.

[0060] FIGS. 19A-19B show amplification curves of samples to be tested.

[0061] FIG. 19A shows amplification curves obtained when the samples to be tested are amplified with GAP-1 and GAP-2 as primers, and FIG. 19B shows amplification curves obtained when the samples to be tested are amplified with 5AOX and 3AOX as primers.

[0062] FIGS. 20A and 20B show melting curves of samples to be tested.

[0063] FIG. 20A shows melting curves obtained when the samples to be tested are amplified with GAP-1 and GAP-2 as primers, and FIG. 20B shows melting curves obtained when the samples to be tested are amplified with 5AOX and 3AOX as primers.

DETAILED DESCRIPTION OF THE INVENTION

[0064] The invention is further illustrated below in conjunction with specific examples. It should be understood that the examples referred to are merely illustrative of the invention and are not intended to limit the scope of the present invention.

Example 1 Codon Optimization of Recombinant ProDer f1

[0065] Based on the DNA sequence of proDer f1 disclosed in GenBank (GenBank accession no. AB034946.1), as shown in SEQ ID No: 2, the inventors performed codon optimization of the gene to obtain the proDer f1 gene of the present invention of which the nucleotide sequence is as shown in SEQ ID No: 1 and the amino acid sequence is as shown in SEQ ID No: 3. Comparison of each parameter before and after codon optimization of the proDer f1 is as follows:

[0066] 1. Codon Adaptation Index (CAI)

[0067] As can be seen from FIG. 2A, the codon adaptation index (CAI) of the original proDer f1 gene in the Pichia pastoris expression system before codon optimization is 0.76. As can be seen from FIG. 2B, the proDer f1 gene has a CAI index of 0.85 in the Pichia pastoris expression system after codon optimization. Usually, when CAI=1, it is considered that the gene is in the most ideal expression state in the expression system. The lower the CAI index, the worse the expression level of the gene in the host. Thus, it can be seen the gene sequence obtained by codon optimization can increase the expression level of the proDer f1 gene in the Pichia pastoris expression system.

[0068] 2. Optimal Codon Usage Frequency (POP)

[0069] As can be seen from FIG. 3A, based on the Pichia pastoris expression vector, the occurrence percentage of the low-utilization codon (codon with a utilization rate less than 40%) of the proDer f1 gene sequence is 10% before codon optimization. This unoptimized gene uses tandem rare codons that may reduce translation efficiency and even disintegrate a translation assembly. As can be seen from FIG. 3B, the proDer f1 gene has a low utilization codon frequency of 0 in the Pichia pastoris system after codon optimization.

[0070] 3. GC Base Content (GC Curve)

[0071] The ideal distribution region of GC content is 30%-70%, and any peak outside this region will affect transcription and translation efficiency to varying degrees. As can be seen from the comparison of the average GC base content distribution region plots of the proDer f1 gene in FIG. 4A and FIG. 4B, FIG. 4A shows the average GC base content of the proDer f1 gene being 41.85%, and FIG. 4B shows that the peaks of GC content appearing outside the 30%-70% region are removed after optimization, and finally the average GC base content of optimized proDer f1 is 42.80%.

Example 2: Construction of an Expression Plasmid Containing the ProDer f1 Gene

[0072] A sequence of XhoI restriction site was introduced at the 5 end, and a sequence of NotI restriction site was introduced at the 3 end of the codon-optimized proDer f1, and then full gene synthesis was performed. The synthesized gene fragment was constructed into the pUC57 plasmid supplied by GenScript (Nanjing) Co., Ltd., thereby obtaining a plasmid for long-term preservation, denoted as pUC57-proDer f1 plasmid.

[0073] PCR amplification was performed using the pUC57-proDer f1 plasmid as a template, and primers of following sequences:

TABLE-US-00001 upstreamprimer(SEQIDNo:5): M13F: TGTAAAACGACGGCCAGT downstreamprimer(SEQIDNo:6): M13R: CAGGAAACAGCTATGAC

[0074] The total volume of the reaction was 50 L, in which 2.5 L of each primer at a concentration of 10 mol/L was added, 1 L of dNTP at a concentration of 10 mmol/L was added, and 0.5 L DNA polymerase being Q5 (#M0491L, purchased from New England BioLabs) at 2 U/pt was added. The reaction conditions were 98 C. for 5 seconds, 55 C. for 45 seconds, and 72 C. for 30 seconds. After 25 cycles, the product was analyzed by 1.0% agarose gel electrophoresis. The results showed that the product size was consistent with the expected size (915 bp) (results as shown in FIG. 5). The product was digested with XhoI (#R0146S, purchased from New England BioLabs) and NotI (#R0189S, purchased from New England BioLabs), respectively, and electrophoresed on 1% agarose gel to obtain a gene product, which was purified using a DNA gel recovery kit (DP214, purchased from Tiangen Biotech (Beijing) Co., Ltd.). The purified product was ligated to pPICZA plasmid (V173-20, purchased from Invitrogen) with T4 ligase (#M0202S, purchased from New England BioLabs), and transformed into DH5a competent cells (CB101, purchased from Tiangen Biotech (Beijing) Co., Ltd.) and cultured in an LB solid medium containing bleomycin (purchased from Invitrogen) at 37 C. overnight. On the second day, the positive clones were picked and sequenced, and the sequence was found identical to the expected sequence by alignment, thereby obtaining the expression plasmid of codon-optimized proDer f1, denoted as pPICZ-proDer f1 (the plasmid construction was as shown in FIG. 6).

Example 3: Construction of a Pichia pastoris Host Engineering Strain Containing a Recombinant ProDer f1 Gene

[0075] Formulation of YPDS solid medium: the medium was formulated according to the instructions of Easy SelectPichia Expression Kit, Invitrogen, comprising 10 g/L yeast extract, 20 g/L peptone, 20 g/L glucose, 15 g/L agarose, and 182 g/L sorbitol.

[0076] 1. Construction of a Host Engineering Strain Containing Codon-Optimized proDer f1

[0077] Electrocompetent cells were prepared according to the method of instructions of Easy SelectPichia Expression Kit, Invitrogen. The plasmid pPICZ-proDer f1 obtained in Example 2 was linearized with Sac I restriction endonuclease (#R0156S, purchased from New England Biolabs), and precipitated with ethanol. The linearized vector was electrotransformed into competent cells of Pichia pastoris X33. The cells were plated on YPDS solid media and cultured at 30 C. until the transformants grew.

Example 4: Inducible Expression and Identification of Engineering Strains Containing Codon-Optimized ProDer f1 Gene

[0078] Formulation of BMGY medium: the medium was formulated according to the instructions of Easy SelectPichia Expression Kit, Invitrogen, comprising 10 g/L yeast extract, 20 g/L peptone, 3 g/L K.sub.2HPO.sub.4, 11.8 g/L KH.sub.2PO.sub.4, 13.4 g/L YNB, 410.sup.4 g/L biotin, and 10 g/L glycerin.

[0079] Formulation of BMMY medium: the medium was formulated according to the instructions of Easy SelectPichia Expression Kit, Invitrogen, comprising 10 g/L yeast extract, 20 g/L peptone, 3 g/L K.sub.2HPO.sub.4, 11.8 g/L KH.sub.2PO.sub.4, 13.4 g/L YNB, 410.sup.4 g/L biotin, and 5 mL/L methanol.

[0080] 1. Methanol-Induced Expression of an Engineering Strain of Codon-Optimized proDer f1

[0081] The host monoclonal engineering strain obtained in Example 3 was picked into a 5 mL BMGY medium and cultured in a 50 mL sterile centrifuge tube at 30 C. and 220 rpm until OD.sub.600 reaches 1.0-2.0. 1 mL of the culture was stored, and the remaining strain solution was resuspended and transferred to BMMY for induced expression at a small scale, and methanol was supplemented every 24 hours to a final concentration of 1%. One week later, the supernatant of the strain solution was collected by centrifugation, and analyzed by SDS-PAGE gel electrophoresis and Western blotting. Brightness of expressed product bands was observed. FIGS. 7A and 7B are plots of identification of induced expression of gene engineering strains containing proDer f1. As seen from FIGS. 7A and 7B, the proDer f1 protein was significantly expressed in the engineering strain.

Example 5: Purification of Recombinant ProDer f1 Protein

[0082] The Der f1 constructed in this invention is obtained mainly by ion exchange and hydrophobic chromatography purification methods. HiTrap SP FF, HiTrap Q FF, and HiTrap Phenyl HP were selected as the chromatographic packings. The specific steps are as follows:

1. Pretreatment of the Fermentation Broth by Impurity Removal

[0083] The fermentation broth of host engineering strain containing proDer f1 obtained according to Example 4 was centrifuged at a low temperature at 12000 rpm for 15 minutes to collect a supernatant, and the supernatant was dialyzed in a 5 KD dialysis bag against a 25 mM sodium acetate buffer at pH 4.5 for 48 h, and filtered through a 0.45 m filter membrane to obtain a supernatant of the treated fermentation broth.

2. Cation Exchange Chromatography

[0084] The treated fermentation broth of the previous step was loaded on a SPFF cation exchange chromatographic column, wherein the equilibration buffer was 50 mM NaAc at pH 4.5, the elution buffer was 50 mMNaAc and 1.0 M NaCl at pH 4.5, isocratic elution was performed at 12%, 25% and 100%, and the sample peaks were mainly concentrated at the 25% elution peak. FIG. 8A is an ion exchange purification chromatogram of Der f1, and FIG. 8B is an SDS-PAGE analysis plot of Der f1 after ion exchange chromatography.

3. Anion Exchange Chromatography

[0085] The Der f1 protein peak purified in the previous step was collected, and the sample was ultrafiltrated with a 20 mM NaH.sub.2PO.sub.4 solution at pH 6.0, and loaded on a HiTrap Q FF chromatography packing. The equilibration buffer was 20 mM NaH.sub.2PO.sub.4 at pH 6.0, and the elution buffer was 20 mM NaH.sub.2PO.sub.4 and 1.0 M NaCl at pH 6.0. The flow-through peak of Der f1 was collected. The flow-through peak of Der f1 protein was as shown in FIGS. 9A and 9B.

4. Hydrophobic Chromatography

[0086] The flow-through peak of Der f1 from the anion chromatography was collected, and ammonium sulfate was added to a final concentration of 1.5 M. The fermentation broth supernatant treated as above was loaded on a Phenyl HP chromatographic column. The equilibration buffer was 20 mM NaH.sub.2PO.sub.4 and 1.5 M (NH.sub.4).sub.2SO.sub.4 at pH 6.0; the elution buffer was 20 mM NaH.sub.2PO.sub.4 at pH 6.0, isocratic elution was performed at 25%, 50%, 70%, and 100%, and the Der f1 protein is mainly concentrated at the 75% elution peak. FIG. 10A is hydrophobic chromatography purification chromatogram of Der f1, and FIG. 10B is an SDS-PAGE analysis plot of Der f1 after hydrophobic chromatography. The yield of target protein per liter of fermentation broth is as high as 200 mg or more.

Example 6: Analysis of Der f1 Protein Activity

[0087] The purified Der f1 protein was dialyzed against a PBS buffer at pH 7.4, and the protein concentration was determined by a Pierce BCA protein concentration assay kit (Cat No: 23225, purchased from Pierce), and fold-diluted to 250 ng, 125 ng, 62.5 ng, 31.25 ng, and 15.625 ng. The obtained solution was detected for the reactivity with sera of patients allergic to Dermatophagoides farinae by comparing with natural Der f1. FIG. 11 shows that the recombinant Der f1 has substantially identical reactivity with the sera as compared with the natural Der f1, showing that the recombinant Der f1 has a similar biological activity as the natural Der f1.

Example 7: Determination of Gene Copy Number of Recombinant ProDer f1 Engineering Strain

[0088] 1. Inoculation in X33 strain: the strains were cultured in YPD media for 24 h, the X33 genome was extracted by a genomic extraction kit (purchased from Tiangen Biotech (Beijing) Co., Ltd.), and GAP gene was amplified using the X33 genome as a template, and GAP-1 and GAP-2 as primers of which the sequences are as follows:

TABLE-US-00002 upstreamprimer (SEQIDNo:7) GAP-1: GGTATTAACGGTTTCGGACGTATTG downstreamprimer (SEQIDNo:8) GAP-2: GATGTTGACAGGGTCTCTCTCTTGG

[0089] The total volume of the reaction was 50 L, in which 2.5 L of each primer at a concentration of 10 mol/L was added, 1 L of dNTP at a concentration of 10 mmol/L was added, and 0.5 L DNA polymerase being Taq DNA Polymerase (M0267S, purchased from New England BioLabs) at 2 U/L was added. The reaction conditions were 94 C. for 10 minutes, 94 C. for 30 seconds, 55 C. for 30 seconds, 68 C. for 60 seconds, and 68 C. for 5 minutes. After 30 cycles, the product was analyzed by 1.0% agarose gel electrophoresis. The results showed that the product size was consistent with the expected size (400 bp) (results as shown in FIG. 12). The obtained gene product was purified by DNA gel recovery kit (DP214, purchased from Tiangen Biotech (Beijing) Co., Ltd.) and 2 buffer ligated into pGM-T vector kit (VT202-01, purchased from Tiangen Biotech (Beijing) Co., Ltd.). The vector was transformed into the Top10 competent cells (CB104, purchased from Tiangen Biotech (Beijing) Co., Ltd.), and cultured at 37 C. overnight on blue-white screening media. On the next day, white clones were picked and identified by PCR for which the primers used were GAP-1 and GAP-2. The PCR reaction conditions were consistent with the above-mentioned conditions. The obtained product was analyzed by 1.0% agarose gel electrophoresis, and the results showed that the product size is consistent with the expected size (400 bp) (results as shown in FIG. 13). The positive clones were sent to GenScript (Nanjing) Co., Ltd. for sequencing, and the sequence was found completely identical to the expected sequence by alignment, thereby obtaining the T vector clone of GAP gene, denoted as T-GAP. The T-GAP clone having a correct sequence was inoculated in an LB liquid medium at 37 C. overnight, and the plasmid was extracted (using a plasmid mini-extract kit DP103, purchased from Tiangen Biotech (Beijing) Co., Ltd.) to obtain a standard plasmid for real-time quantitative PCR.

[0090] 2. The proDer f1 gene was amplified using the pPICZ-proDer f1 plasmid of Example 2 as a template, and 5 AOX and 3 AOX as primers with the following sequences:

TABLE-US-00003 upstreamprimer(SEQIDNo:9): 5 AOX: GACTGGTTCCAATTGACAAGC downstreamprimer(SEQIDNo:10): 3 AOX: GGCAAATGGCATTCTGACAT

[0091] The total volume of the reaction was 50 in which 2.5 L of each primer at a concentration of 10 mol/L was added, 1 L of dNTP at a concentration of 10 mmol/L was added, and 0.5 L DNA polymerase being Taq DNA Polymerase (# M0267S, purchased from New England BioLabs) at 2 U/L was added. The reaction conditions were 94 C. for 10 minutes, 94 C. for 30 seconds, 49 C. for 30 seconds, and 68 C. for 60 seconds, and 68 C. for 5 minutes. After 30 cycles, the product was analyzed by 1.0% agarose gel electrophoresis. The results showed that the product size was consistent with the expected size (1500 bp) (results as shown in FIG. 14). The obtained gene product was purified by DNA gel recovery kit (DP214, purchased from Tiangen Biotech (Beijing) Co., Ltd.) and ligated into pGM-T vector kit (VT202-01, purchased from Tiangen Biotech (Beijing) Co., Ltd.). The vector was transformed into the Top10 competent cells (CB104, purchased from Tiangen Biotech (Beijing) Co., Ltd.), and cultured at 37 C. overnight on blue-white screening media. On the next day, white clones were picked and identified by PCR for which the primers used were 5AOX and 3AOX. The PCR reaction conditions were consistent with the above-mentioned conditions. The obtained product was analyzed by 1.0% agarose gel electrophoresis, and the results showed that the product size is consistent with the expected size (1500 bp) (results as shown in FIG. 15). The positive clones were sent to GenScript (Nanjing) Co., Ltd. for sequencing, and the sequence was found completely identical to the expected sequence by alignment, thereby obtaining the T vector clone of proDer f1, denoted as T-proDer f1. The T-proDer f1 clone having a correct sequence was inoculated in an LB liquid medium at 37 C. overnight, and the plasmid was extracted using a plasmid mini-extract kit (DP103, purchased from Tiangen Biotech (Beijing) Co., Ltd.) to obtain a standard plasmid for real-time quantitative PCR.

3. Calculation of Gene Copy Number:

[0092] The concentration (ng/L) of the standard plasmid was determined by a nucleic acid microanalyzer (Nanodrop2000, ThermoFisher). Copy numbers of GAP and proDer f1 were calculated according to the following formula:

Copies/u=(6.0210.sup.23)(ng/l10.sup.9)/(DNA length660)

4. Processing Samples to be Tested

[0093] The pPICZ-proDer f1-X33 engineering strain was inoculated in YPD liquid media at 30 C. overnight; and the genome was extracted the next day, and its concentration (ng/L) and purity were determined by a nucleic acid quantitative microanalyzer.

5. Establishment of a Standard Curve

[0094] The standard plasmids of T-GAP and T-proDer f1 with known copy numbers were gradiently diluted to 10.sup.8, 10.sup.7, 10.sup.6, 10.sup.5, 10.sup.4, and 10.sup.3 copies/l, respectively. The fluorescent quantitative PCR were performed using GAP-1 and GAP-2, 5 AOX and 3 AOX as primers, respectively. FIG. 16A shows amplification curves of the standard plasmid T-GAP, FIG. 16B shows amplification curves of the standard plasmid T-proDer f1, FIG. 17A shows melting curves of the standard plasmid T-GAP, and FIG. 17B shows melting curves of the standard plasmid T-proDer f1. Each gradient was assayed 3 times to verify the repeatability of the standard curve. Standard curves were established with the Ct values as the ordinate and the starting template copy numbers as the abscissa. FIG. 18A shows a standard curve of the standard plasmid T-GAP, and FIG. 18B shows a standard curve of the standard plasmid T-proDer f1.

6. Determination of Copy Number of ProDer f1 Gene in Recombinant Strains

[0095] The genome sample of extracted pPICZ-proDer f1-X33 was serially 10-fold-diluted to obtain four gradients of stock solution, 10.sup.1, 10.sup.2, and 10.sup.3. Fluorescent quantitative PCR was performed using GAP-1 and GAP-2, 5 AOX and 3 AOX as primers, and each gradient was assayed three times. FIG. 19A shows amplification curves of the samples to be tested with GAP-1 and GAP-2 as primers, FIG. 19B shows amplification curves of the samples to be tested with 5 AOX and 3 AOX as primers, FIG. 20A shows melting curves of the samples to be tested with GAP-1 and GAP-2 as primers, and FIG. 20B shows melting curves of the samples to be tested with 5 AOX and 3 AOX as primers. The GAP gene exists in Pichia pastoris in a single copy. Therefore, the copy number of the GAP gene can be used to characterize the initial copy number of the genome in the template. The ratio of the copy number of the proDer f1 gene to the copy number of the GAP gene is the copy number of proDer f1 gene in the Pichia pastoris genome. Table 1 shows the detection results of copy number of the proDer f1 gene in the Pichia pastoris gene engineering strain, the detected copy number is between 4.85 and 6.02, and finally the copy number of the proDer f1 gene in the recombinant strain was averaged to eliminate the system error and determined to be 5.

TABLE-US-00004 TABLE 1 Results of copy number of proDer f1 in the genome detected by real-time fluorescent quantitative PCR Average Ct value gene copy number (10.sup.N)Copy number of proDer f1 gene in Pichia pastoris genome Copy number of the proDer f1 gene/copy DNA GAP proDer f1 GAP proDer f1 number of the concentration gene gene gene gene GAP gene Stock 18.802 23.000 6.84 5.96 6.02 solution 10.sup.1 19.939 24.382 6.29 5.57 5.46 10.sup.2 23.650 24.704 5.17 5.29 5.07 10.sup.3 27.966 24.876 3.62 5.19 4.85

Example 8: Analysis of the Acting Elements in the ProDer f1 Genome

[0096] There is no stable additional plasmid in Pichia pastoris, the expression vector is homologously recombined with the host chromosome, and the exogenous gene expression framework is fully integrated into the chromosome to realize the expression of the exogenous gene; the typical Pichia pastoris expression vector contains a regulatory sequence of alcohol oxidase gene, and contains the main structures comprising AOX promoter, multiple cloning site, transcription termination and polyA formation gene sequence (TT), screening markers and the like. The promoter is a cis-element for gene expression regulation and an important element for the genetically engineered expression vector. The important role of the promoter at the transcriptional level determines the gene expression level.

[0097] The proDer f1 genome was extracted according to the method of Example 7, and the proDer f1 gene was amplified from the genome using 5 AOX and 3 AOX as primers. The obtained samples were sent to GenScript (Nanjing) Co., Ltd. to detect the acting element before and after the proDer f1 gene which was inserted into the genome. The results of genome sequencing indicated that the proDer f1 gene expression framework was integrated into the chromosome of Pichia pastoris by a single cross-insertion, which enabled the proDer f1 gene to express the gene using the AOX promoter on the yeast chromosome, and thus the expression level was higher.

[0098] Generally, the closer the first ATG of the exogenous coding sequence to the ATG of AOX1, the better the expression effect. In the gene construction, the inventors chose an enzyme cleavage site closest to the ATG of AOX1, and found that the proDer f1 gene was away from ATG of AOX1 only by 242 bp. In addition, the alpha-factor signal peptide and Kozak sequence GCCACCATGG were added in front of proDer f1 gene, and the signal peptide and the sequence can greatly improve transcription and translation efficiency and increase expression efficiency of proDer f1 gene in eukaryotes.