Cloning of a novel lipase gene, <Emphasis Type="Italic">lipJ08,</Emphasis> from <Emphasis Type="Italic">Candida rugosa</Emphasis> and expression in <Emphasis Type="Italic">Pichia pastoris</Emphasis> by codon optimization

A novel lipase gene, lipJ08, was cloned from Candida rugosa ATCC14830, along with the already reported five lipase genes (lip1–lip5). Nucleotide sequencing indicated that the lipJ08 gene contains a 1650 bp open reading frame (ORF) without introns. The deduced amino acid sequence corresponds to 534 amino acid residues, including a putative signal sequence of 15 amino acid residues. Seventeen of the non-universal serine codons (CTG) of lipJ08 were converted into universal serine codons (TCT) by PCR-based mutagenesis. The native and codon-optimized lipJ08 genes were expressed in Pichia pastoris. The hydrolytic activity of the recombinant LIPJ08 was 4.7 U/ml, whereas the activity of the recombinant wild-type lipase could not be detected.     

Cloning,expression and characterization of a new lipase from <Emphasis Type="Italic">Yarrowia lipolytica</Emphasis>

Bioinformatic analysis of the Yarrowia lipolytica CLIB122 genome has revealed 18 putative lipase genes all of which were expressed in Escherichia coli and screened for hydrolyzing activities against p-nitrophenyl-palmitate. One positive transformant containing an ORF of 1,098 bp encoding a protein of 365 amino acids was obtained. To characterize its enzymatic properties, the lipase gene was functionally expressed in Pichia pastoris. The resulting lipase exhibited the highest activity towards p-NP-decanoate at pH 7 and 35°C. In addition, the new lipase had a lower optimal temperature and pH compared to other Y. lipolytica lipases. It was noticeably enhanced by Ca²⁺, but was inhibited by PMSF, Hg²⁺ and Ni²⁺. The new lipase displayed the 1,3-specificity for triolein.     

Secretion of pro- and mature Rhizopus arrhizus lipases by Pichia pastoris and properties of the proteins

The lipases of Rhizopus spp. share a high 1,3-regiospecificity toward triacylglycerols, which makes them important enzymes in lipid modification. In the present study, the extracellularly active production of recombinant Rhizopus arrhizus lipase was carried out with genes encoding the mature region (mRAL) and the mRAL having the prosequence (ProRAL) in Pichia pastoris. Two transformed P. pastoris clones containing the multicopy of mRAL and ProRAL genes were separately selected for the production of recombinant enzymes. In a fed-batch cultivation, where methanol feeding was controlled by an on-line methanol analyzer, the supernatant contained 91 mg/L recombinant pro-form lipase (rProRAL) and 80 mg/L recombinant mature lipase (rRAL) after 92 h of cultivation. rProRAL and rRAL were purified by ultrafiltration, SP-Sepharose Rast Flow chromatography, and Butyl-Sepharose Fast Flow chromatography. Molecular weights of rProRAL and rRAL are 32 kDa and 29 kDa, respectively. The amino-terminal analysis showed that the 32-kDa protein was mRAL attached with 28 amino acids of the carboxy-terminal part of the prosequence (rPro28RAL). The specific lipase activities of mRAL attached with 28 amino acids of the carboxy-terminal part of the prosequence (rPro28RAL) and rRAL were 1543 U/mg and 2437 U/mg. The rPro28RAL was more stable than rRAL at pH 4.0–7.0, whereas rRAL was more stable at pH 7.0–10.0. The rPro28RAL had the highest lipase activity toward tributyrin (C₄), whereas rRAL had the highest lipase activity toward tricaprylin (C₈).     

Cloning,expression, and characterization of serine protease from thermophilic fungus <Emphasis Type="Italic">Thermoascus aurantiacus</Emphasis> var. <Emphasis Type="Italic">levisporus</Emphasis>

The serine protease gene from a thermophilic fungus Thermoascus aurantiacus var. levisporus, was cloned, sequenced, and expressed in Pichia pastoris and the recombinant protein was characterized. The full-length cDNA of 2,592 bp contains an ORF of 1,482 bp encoding 494 amino acids. Sequence analysis of the deduced amino acid sequence revealed high homology with subtilisin serine proteases. The putative enzyme contained catalytic domain with active sites formed by three residues of Aspl83, His215, and Ser384. The molecular mass of the recombinant enzyme was estimated to be 59.1 kDa after overexpression in P. pastoris. The activity of recombinant protein was 115.58 U/mg. The protease exhibited its maximal activity at 50°C and pH 8.0 and kept thermostable at 60°C, and retained 60% activity after 60 min at 70° C. The protease activity was found to be inhibited by PMSF, but not by DTT or EDTA. The enzyme has broad substrate specificity such as gelatin, casein and pure milk, and exhibiting highest activity towards casein.     

Cloning of exoinulinase gene from Penicillium janthinellum strain B01 and its high-level expression in Pichia pastoris

Aims: The aim of this study is to improve exoinulinase production by expression of a cloned exoinulinase gene inuA1 (GenBank accession no. JF961344 ) from Penicillium janthinellum strain B01 in Pichia pastoris. Methods and Results: A full‐length cDNA of exoinulinase gene (inuA1) was cloned from P. janthinellum strain B01 using RACE PCR. An open reading frame (ORF) of 2115 bp is interrupted by a single intron of 67 bp. The fragment encodes a signal peptide with 20 amino acids and a mature protein with 684 amino acids. The inuA1 was subcloned to the pPICZαC expression vector and succesfully over‐expressed in Pichia pastoris X‐33. The highest activity of exoinlinase reached 272·8 U ml^?1 in the fermentation liquid. It was c. 11‐fold of that produced by wild‐strain B01. A large amount of fructose was identified after the hydrolysis of inulin with the crude recombinant exoinulinase. The recombinant exoinulinase was purified and characterized. The molecular weight of the purified recombinant exoinulianse was 100 kDa. The mass spectrometry result indicated that the purified protein was indeed recombinant exoinulinase. The optimal pH and temperature of the purified recombinant exoinulianse were 4·5 and 50°C, respectively. Conclusions: An exoinulinase gene of P. janthinellum strain B01 was cloned, sequenced and over‐expressed successfully in P. pastoris. Significance and Impact of the Study: Only a few genes have been cloned from P. janthinellum because its molecular biology is poorly understood. In this study, we cloned and over‐expressed inuA1 gene of P. janthinellum in P. pastoris. This recombinant exoinulinase can be used to hydrolyse inulin to produce fructose and facilitate the biofuel production from inulin resources.     

Cloning,Expression and Characterization of a Lipase Encoding Gene from Human Oral Metagenome

The human oral metagenomic DNA cloned into plasmid pUC19 was used to construct a DNA library in Escherichia coli. Functional screening of 40,000 metagenomic clones led to identification of a clone LIP2 that exhibited halo on tributyrin agar plate. Sequence analysis of LIP2 insert DNA revealed a 939 bp ORF (omlip1) which showed homology to lipase 1 of Acinetobacter junii SH205. The omlip1 ORF was cloned and expressed in E. coli BL21 (DE3) using pET expression system. The recombinant enzyme was purified to homogeneity and the biochemical properties were studied. The purified OMLip1 hydrolyzed p-nitrophenyl esters and triacylglycerol esters of medium and long chain fatty acids, indicating the enzyme is a true lipase. The purified protein exhibited a pH and temperature optima of 7 and 37 °C respectively. The lipase was found to be stable at pH range of 6–7 and at temperatures lower than 40 °C. Importantly, the enzyme activity was unaltered, by the presence or absence of many divalent cations. The metal ion insensitivity of OMLip1offers its potential use in industrial processes.     

Cloning,expression, and characterization of Aureobasidium melanogenum lipase in Pichia pastoris

cDNA of Aureobasidium melanogenum lipase comprises 1254 bp encoding 417 amino acids, whereas genomic DNA of lipase comprises 1311 bp with one intron (57 bp). The lipase gene contains a putative signal peptide encoding 26 amino acids. The A. melanogenum lipase gene was successfully expressed in Pichia pastoris. Recombinant lipase in an inducible expression system showed the highest lipase activity of 3.8 U/mL after six days of 2% v/v methanol induction. The molecular mass of purified recombinant lipase was estimated as 39 kDa using SDS-PAGE. Optimal lipase activity was observed at 35–37 °C and pH 7.0 using p-nitrophenyl laurate as the substrate. Lipase activity was enhanced by Mg²⁺, Mn²⁺, Li⁺, Ca²⁺, Ni²⁺, CHAPS, DTT, and EDTA and inhibited by Hg²⁺, Ag⁺, SDS, Tween 20, and Triton X-100. The addition of 10% v/v acetone, DMSO, p-xylene, and octanol increased lipase activity, whereas that of propanol and butanol strongly inhibited it.     

Cloning,characterization and expression of the extracellular lipase gene from <Emphasis Type="Italic">Aureobasidium </Emphasis><Emphasis Type="Italic">pullulans</Emphasis> HN2-3 isolated from sea saltern

The extracellular lipase structural gene was isolated from cDNA of Aureobasidium pullulans HN2-3 by using SMART^TM RACE cDNA amplification kit. The gene had an open reading frame of 1245 bp long encoding a lipase. The coding region of the gene was interrupted by only one intron (55 bp). It encodes 414 amino acid residues of a protein with a putative signal peptide of 26 amino acids. The protein sequence deduced from the extracellular lipase structural gene contained the lipase consensus sequence (G-X-S-X-G) and three conserved putative N-glycosylation sites. According to the phylogenetic tree of the lipases, the lipase from A. pullulans was closely related to that from Aspergillus fumigatus (XP_750543) and Neosartorya fischeri (XP_001257768) and the identities were 50% and 52%, respectively. The mature peptide encoding cDNA was subcloned into pET-24a (+) expression vector. The recombinant plasmid was expressed in Escherichia coli BL21(DE3). The expressed fusion protein was analyzed by SDS-PAGE and western blotting and a specific band with molecular mass of about 47 kDa was found. Enzyme activity assay verified the recombinant protein as a lipase. A maximum activity of 0.96 U/mg was obtained from cellular extract of E. coli BL21(DE3) harboring pET-24a(+)LIP1. Optimal pH and temperature of the crude recombinant lipase were 8.0 and 35 °C, respectively and the crude recombinant lipase had the highest hydrolytic activity towards peanut oil.     

Isolation of a Novel Lipase Gene from <Emphasis Type="Italic">Serratia liquefaciens</Emphasis> S33 DB-1, Functional Expression in <Emphasis Type="Italic">Pichia pastoris</Emphasis> and its Properties

A new lipase gene designated as SlLipA was isolated from Serratia liquefaciens S33 DB-1 by the genomic-walking method. The cloned gene contained an open reading frame (ORF) of 1,845 bp encoding 615 amino acids with a conserved GXSXG motif. Genome sequence analysis showed that an aldo/keto reductase gene closed to the SlLipA gene. The lipase gene was cloned into the expression vector pPICZαA and successfully integrated into the heterologous host, methylotrophic yeast Pichia pastoris GS115. Five transformants could be expressed as secreted recombinant proteins with the high activity on Triglyceride–Agarose plate and as candidates to produce the recombinant enzyme. A C-terminal His tag was used for its purification. The lipase activity of different transformants against substrate para-nitrophenyl laurate (p-NPL) varied from 14 to 16 U ml⁻¹. For the substrates para-nitrophenyl caprate (p-NPC), p-NPL, para-nitrophenyl myristate (p-NPM), para-nitrophenyl palmitate (p-NPP), and para-nitrophenyl stearate (p-NPS), the specific activity was shown to be preferred to long acyl chain length of p-NPS.     

Isolation and characterization of a new alkali-thermostable lipase cloned from a metagenomic library

The construction of a cosmid library from the biomass produced in an enriched Sequencing Fed-Batch Reactor allowed the isolation of a new lipase by functional screening. The open reading frame of 928 bp encoded a polypeptide of 308 amino acids with a molecular mass of 32.6 kDa. The amino acid sequence analysis revealed the presence of the conserved pentapeptide GXSXG essential for lipase activity. Alignment with known sequences of proteins showed no more than 52% identity with different lipases, confirming the discovery of a novel gene sequence. The lipase was cloned and expressed in Streptomyces lividans and further purified by a combination of hydrophobic interaction and size-exclusion chromatography. Spectrophotometric assays with different p-nitrophenyl esters demonstrated a preference for long-length acyl chains, especially p-nitrophenylmyristate (C14). Moreover, the enzyme presented an optimal activity at 60°C and at alkaline pH of 10.5.     

Characterization of an extracellular lipase in Burkholderia sp. HY-10 isolated from a longicorn beetle

Burkholderia sp. HY-10 isolated from the digestive tracts of the longicorn beetle, Prionus insularis, produced an extracellular lipase with a molecular weight of 33.5 kDa estimated by SDS-PAGE. The lipase was purified from the culture supernatant to near electrophoretic homogenity by a one-step adsorption-desorption procedure using a polypropylene matrix followed by a concentration step. The purified lipase exhibited highest activities at pH 8.5 and 60 degrees . A broad range of lipase substrates, from C4 to C18 rho-nitrophenyl esters, were hydrolyzed efficiently by the lipase. The most efficient substrate was rho-nitrophenyl caproate (C6). A 2485 bp DNA fragment was isolated by PCR amplification and chromosomal walking which encoded two polypeptides of 364 and 346 amino acids, identified as a lipase and a lipase foldase, respectively. The N-terminal amino acid sequence of the purified lipase and nucleotide sequence analysis predicted that the precursor lipase was proteolytically modified through the secretion step and produced a catalytically active 33.5 kDa protein. The deduced amino acid sequence for the lipase shared extensive similarity with those of the lipase family I.2 of lipases from other bacteria. The deduced amino acid sequence contained two Cystein residues forming a disulfide bond in the molecule and three, well-conserved amino acid residues, Ser131, His330, and Asp308, which composed the catalytic triad of the enzyme.     

Characterization and gene cloning of a novel β-mannanase from alkaliphilic <Emphasis Type="Italic">Bacillus</Emphasis> sp. N16-5

An alkaline -mannanase was purified to homogeneity from a culture broth of alkaliphilic Bacillus sp. N16-5. The enzyme had optimum activity at pH 9.5 and 70°C. It was composed of a single polypeptide chain with a molecular weight of 55 kDa deduced from SDS-PAGE, and its isoelectric point was around pH 4.3. The enzyme efficiently hydrolyzed galactomannan and glucomannan, producing a series of oligosaccharides and monosaccharides. The -mannanase gene (manA) contained an open reading frame (ORF) of 1,479 bp, encoding a 32-amino acids signal peptide, and a mature protein of 461 amino acids, with a calculated molecular mass of 50,743 Da. Strain N16-5 ManA, deduced from the manA ORF, exhibited relatively high amino acid similarity to the members of the glycosyl hydrolase family 5. The eight conserved active-site amino acids in family 5 glycosyl hydrolase were found in the deduced amino acid sequence of strain N16-5 ManA.     

Secretory expression and characterization of a highly Ca-activated thermostable L2 lipase

Thermostable lipases are important biocatalysts, showing many interesting properties with industrial applications. Previously, a thermophilic Bacillus sp. strain L2 that produces a thermostable lipase was isolated. In this study, the gene encoding for mature thermostable L2 lipase was cloned into a Pichia pastoris expression vector. Under the control of the methanol-inducible alcohol oxidase (AOX) promoter, the recombinant L2 lipase was secreted into the culture medium driven by the Saccharomyces cerevisiae α-factor signal sequence. After optimization the maximum recombinant lipase activity achieved in shake flasks was 125 U/ml. The recombinant 44.5 kDa L2 lipase was purified 1.8-fold using affinity chromatography with 63.2% yield and a specific activity of 458.1 U/mg. Its activity was maximal at 70 °C and pH 8.0. Lipase activity increased 5-fold in the presence of Ca²⁺. L2 lipase showed a preference for medium to long chain triacylglycerols (C₁₀–C₁₆), corn oil, olive oil, soybean oil, and palm oil. Stabilization at high temperature and alkaline pH as well as its broad substrate specificity offer great potential for application in various industries that require high temperature operations.     

Cloning and characterization of lactate dehydrogenase C4 from Ochotona curzoniae

Lactate dehydrogenase C4 (LDH-C4) is considered to be a good target protein for the development of contraceptive drugs. To develop contraceptive rodenticide against pika (Ochotona curzoniae) LDH-C4, the pika LDH-C gene was cloned and expressed in Escherichia coli. The recombinant protein was purified and characterized. The cDNA of pika LDH-C gene was cloned by the RACE method. The cDNA was 1498 bp in length and contained an ORF of 996 bp, which encoded a polypeptide of 332 amino acids. The ORF of pika LDH-C was introduced in E. coli and expressed with no fusion tags added. The recombinant LDH-C4 protein was purified by heating, affinity chromatography and ion-exchange chromatography. The recombinant pika LDH-C4 was a tetramer with a molecular weight of approximately 140 kDa, and it had temperature-dependent catalytic activity, as it was thermally stable up to 60°C. The optimal pH values in the forward and backward reactions were around 7.48 and 10.28, respectively. The apparent Michaelis constants for pyruvate and lactate were 51.2 ± 3.8 and 8568.8 ± 409 μM respectively. The inhibition constant for oxalic acid was 11.8 ± 3.5 mM. This study laid a solid foundation for contraceptive rodenticide development against pika LDH-C4.     

Molecular and enzymatic characterization of a subfamily I.4 lipase from an edible oil-degrader <Emphasis Type="Italic">Bacillus</Emphasis> sp. HH-01

An edible-oil degrading bacterial strain HH-01 was isolated from oil plant gummy matter and was classified as a member of the genus Bacillus on the basis of the nucleotide sequence of the 16S rRNA gene. A putative lipase gene and its flanking regions were cloned from the strain based on its similarity to lipase genes from other Bacillus spp. The deduced product was composed of 214 amino acids and the putative mature protein, consisting of 182 amino acids, exhibited 82% amino acid sequence identity with the subfamily I.4 lipase LipA of Bacillus subtilis 168. The recombinant product was successfully overproduced as a soluble form in Escherichia coli and showed lipase activity. The gene was, therefore, designated as lipA of HH-01. HH-01 LipA was stable at pH 4–11 and up to 30°C, and its optimum pH and temperature were 8–9 and 30°C, respectively. The enzyme showed preferential hydrolysis of the 1(3)-position ester bond in trilinolein. The activity was, interestingly, enhanced by supplementing with 1 mM CoCl₂, in contrast to other Bacillus lipases. The lipA gene seemed to be constitutively transcribed during the exponential growth phase, regardless of the presence of edible oil.     

Complete Sequence of Plasmid pMP1 from the Marine Environmental Vibrio vulnificus and Location of Its Replication Origin

A novel cryptic plasmid, pMP1, from an environmental Vibrio vulnificus MP-4 isolated from Mai Po Nature Reserve in Hong Kong, has been characterized. The 7.6-kb plasmid had guanine–cytosine content of 40.03% and encoded four open reading frames (ORFs) with >100 amino acids. The predicted protein of ORF1 contained 478 amino acids showing 29% identity and 50% similarity over 309 amino acids to the integrase of Vibrio cholerae phage VP2. ORF2 encoded a putative protein of 596 amino acids, which were 23% identity and 42% similarity over 455 amino acids to the tail tape measure protein TP901 of Chromohalobacter salexigens phage. ORF3 and ORF4 encoded putative proteins of 103 and 287 amino acids, respectively, but showed no homologies to any known proteins. Further experiments indicated that a 3.2-kb fragment from EcoRI digestion could self-replicate. Analysis indicated that a sequence upstream of ORF4 had the features characteristic of theta-type replicons: AT-rich region, six potential direct repeats (iterons) spaced approximately two DNA helical turn apart (about 23 bp), two copies of 9 bp dnaA boxes, three Dam methylation sites, and five inverted repeats. Complementation experiments confirmed that the protein encoded by ORF4 was required for plasmid replication. We propose that ORF4 encode a new type of Rep protein and pMP1 is a new type of theta plasmid.     

Gene cloning,overexpression and characterization of a novel organic solvent tolerant and thermostable lipase from Galactomyces geotrichum Y05

Although the lipase of Geotrichum candidum has been extensively reported, little attention has been focused on molecular genetic and biochemical characterizations of Galactomyces geotrichum lipases. A lipase gene from G. geotrichum Y05 was cloned from both genomic DNA and cDNA sources. Nucleotide sequencing revealed that the ggl gene has an ORF of 1692 bp without any introns, encoding a protein of 563 amino acid residues, including a potential signal sequence of 19 amino acid residues. The amino acid sequence of this lipase showed 86% identity to lipase of Trichosporon fermentans WU-C12. The mature lipase gene was subcloned into pPIC9K vector, and overexpressed in methylotrophic Pichia pastoris GS115. Active lipase was accumulated to the level of 100.0 U/ml (0.4 mg/ml) in the shake-flask culture, 10.4-fold higher than the activity of the original strain (9.6 U/ml). This yield dramatically exceeds that previously reported with 23–50 U/ml, 0.06 mg/ml and 0.2 mg/ml. The purified lipase exhibited several properties of significant industrial importance, such as pH and temperature stability, wide organic solvent tolerance and broad hydrolysis on vegetable oils. Such a combination of properties makes it a promising candidate for its application in non-aqueous biocatalysis, such as biodiesel production, selective hydrolysis or esterification for enrichment of PUFAs and oil-contaminated biodegradation, which have been drawn considerable attention currently.     

Isolation and biochemical characterization of two lipases from a metagenomic library of China Holstein cow rumen

Two novel lipase genes RlipE1 and RlipE2 which encoded 361- and 265-amino acid peptides, respectively, were recovered from a metagenomic library of the rumen microbiota of Chinese Holstein cows. A BLAST search revealed a high similarity (90%) between RlipE2 and a carboxylesterase from Thermosinus carboxydivorans Nor1, while there was a low similarity (below 50%) between RlipE1 and other lipases. Phylogenetic analysis indicated that RlipE2 clustered with the lipolytic enzymes from family V while RlipE1 clustered with six other putative bacterial lipases which might constitute a new subfamily. The recombinant lipases were thermally unstable and retained 60% activity over a pH range of 6.5-8.5. Substrate specificity assay indicated that both enzymes had higher hydrolytic activity toward laurate (C₁₂), palmitate (C₁₆) and stearate (C₁₈). The novel phylogenetic affiliation and high specificity of both enzymes for long-chain fatty acid make them interesting targets for manipulation of rumen lipid metabolism.     

Characterization of a novel lipolytic enzyme from Aspergillus oryzae

In this study, we report the characterization of a protein from Aspergillus oryzae, exhibiting sequence identity with paraben esterase from the genus Aspergillus. The coding region of 1,586 bp, including a 77-bp intron, encoded a protein of 502 amino acids. The gene without the signal peptide of 19 amino acids was cloned into a vector, pPICZαC, and expressed successfully in Pichia pastoris as an active extracellular protein. The purified recombinant protein had pH and temperature optima of 7.0–8.0 and 30 °C, respectively, and was stable at the pH range of 7.0–10.0 and up to 40 °C. The optimal substrate for hydrolysis by the purified recombinant protein, among a panel of α-naphthyl esters (C2–C16), was α-naphthyl butyrate (C4), with activity of 0.16 units/mg protein. The considerable hydrolytic activity of the purified recombinant enzyme toward tributyrin was determined. However, no paraben esterase activity was detected toward the ethyl, propyl, and butyl esters of 4-hydroxybenzoic acid. In addition, no activity was detected toward the methyl esters of ferulic, p-coumaric, caffeic, and sinapic acids that would indicate feruloyl esterase activity.     

Isolation and characterization of two novel plasmids pCYM01 and pCYM02 of Cylindrospermum stagnale

Cyanobacteria play a vital role in supplying nitrogen into the soil and aquatic ecosystem. It has an extra chromosomal DNA, whose role is not yet defined well. Isolation and characterization of extra chromosomal DNA in cyanobacteria might help to understand its survival mechanism. Cylindrospermum stagnale isolated (and deposited in NRMCF 3001) from soil showed presence of four plasmids namely pCYLM01, pCYLM02, pCYLM03, and pCYLM04. The following plasmids pCYLM01 and pCYLM02 were subjected to restriction digestion using HindIII restriction enzyme and cloned into pBlueScriptSK(-) vector. The sequence of pCYLM01 contained 4 potential open reading frames (ORFs) that have amino acids in the range of 59–299. Among them, ORF1 shows high sequence homology to the bacterial replication initiator family protein as evident from BLASTP analysis. The analysis of 4359 bp plasmid pCYLM02 sequence revealed 7 ORFs which are longer than 50 amino acids in length. The ORF2 of pCYLM02 has 243 amino acids and is represented in the plasmid sequence from 3045 to 3776 bp. The ORF3 of pCYLM02 corresponds to the plasmid sequence from 2323 to 2976 and codes for a putative protein of 217 amino acids long. A number of small ORFs below 50 bp were also found in the sequence analysis.