Cost-Effective Method of Gene Synthesis by Sequencing from Microchip-Derived Oligos for Droplet Cloning

Gene synthesis has provided important contributions in various fields including genomics and medicine. Current genes are 7 - 30 cents depending on the assembly and sequencing methods performed. Demand for gene synthesis has been increasing for the past few decades, yet available methods remain expensive. A solution to this problem involves microchip-derived oligonucleotides (oligos), an oligo pool with a substantial number of oligo fragments. Microchips have been proposed as a tool for gene synthesis, but this approach has been criticized for its high error rate during sequencing. This study tests a possible cost-effective method for gene synthesis utilizing fragment assembly and golden gate assembly, which can be employed for quicker manufacturing and efficient execution of genes in the near future. The droplet method was tested in two trials to determine the viability of the method through the accuracy of the oligos sequenced. A preliminary research experiment was performed to determine the efficacy of oligo lengths ranging from two to four overlapping oligos through Gibson assembly. Of the three oligo lengths tested, only two fragment oligos were correctly sequenced. Two fragment oligos were used for the second experiment, which determined the efficacy of the droplet method in reducing gene synthesis cost and speed. The first trial utilized a high-fidelity polymerase and resulted in 3% correctly sequenced oligos, so the second trial utilized a non-high-fidelity polymerase, resulting in 8% correctly sequenced oligos. After calculating, the cost of gene synthesis lowers down to 0.8 cents/base. The final calculated cost of 0.8 cents/base is significantly cheaper than other manufacturing costs of 7 - 30 cents/base. Reducing the cost of gene synthesis provides new insight into the cost-effectiveness of present technologies and protocols and has the potential to benefit the fields of bioengineering and gene therapy.

Gene Cloning and Bioinformatics Analysis of phoR Gene from Vibrio alginolyticus HY9901

PhoR is a histidine kinase in a two-component regulatory system that regulates phosphorus metabolic pathways and undertakes the key mission of information transmission in pathogenic bacteria.The full-length phoR gene was successfully cloned from the Vibrio alginolyticus HY9901 strain.A comprehensive analysis of the cloned gene was conducted using bioinformatics.Sequence analysis revealed that the total length of the phoR gene(GenBank accession No.:KJ958404.1)is 1299 bp,with the coding region containing a total of 432 amino acid residues.The phylogenetic tree of PhoR revealed that it belongs to the same subclade as V.diabolicus.The SMART program was employed for the purpose of functional domain prediction,which revealed that PhoR possesses three major functional domains:PAS(amino acids 98-166),HisKA(amino acids 205-272),and HATPase_c(amino acids 317-429).

Molecular Cloning of sodB Gene from Vibrio alginolyticus HY9901 and Its Bioinformatics Analysis

Vibrio alginolyticus is a zoonotic bacterium.A pair of specific primers was designed using the sodB gene sequence of Vibrio alginolyticus HY9901 in order to amplify the full length of the gene by PCR.The results indicated that the total length of the sodB gene was 585 bp and that it could encode 194 amino acids.The predicted amino acid sequence derivation indicated that the molecular weight of the protein was approximately 21.56 kDa,with an isoelectric point of 4.95.Upon prediction of the N-terminal signal peptide structure of the protein,no significant signal peptide cleavage site was observed,indicating that the protein lacked both a signal peptide and a transmembrane region.The amino acid sequence contained an N-glycosylation site,a casein kinase II phosphorylation site,a microsomal C-terminal target signal site,and a manganese and iron superoxide dismutase signal site.The probability of intracytoplasmic localization of the SodB protein was 56.5%,which was analyzed according to the subcellular localization of the protein.The amino acid sequence of the sodB gene of V.alginolyticus exhibited 98%-100%homology to other Vibrio species,clustering into the same subfamily with V.parahaem,indicating a relatively close relationship between them.In the prediction of protein structure,the proportions ofα-helix,random coil,β-sheet,and extended strand were 48.45%,30.41%,5.67%,and 15.46%,respectively.The similarity to template 1dt0.1.A reached 71.58%.A PTM site analysis revealed the presence of phosphorylation,glycosylation,ubiquitination,sumoylation,acetylation,and methylation modification sites,as well as the absence of lactylation modification sites.

Molecular Cloning and Bioinformatics Analysis of msrA Gene from Vibrio alginolyticus Strain HY9901

[Objectives]This study was conducted to understand the structure and function of MsrA protein.[Methods]With Vibrio alginolyticus HY9901 as the object of study,primers were designed to amplify the full-length gene of msrA,and its bioinformatics analysis was carried out.[Results]The full length of msrA gene was 639 bp,encoding 212 amino acids,and its theoretical molecular weight was about 23729.60 Da.The protein had a stable structure,and it was hydrophobic overall.The structure of signal peptides at the N terminal of the amino acid sequence was predicted,and it was found that there was no signal peptide cleavage site and no transmembrane region.The amino acid sequence of MsrA contained multiple signal binding sites.Protein subcellular localization showed that MsrA protein was most likely located in the cytoplasm.Homology analysis showed that MsrA of V.alginolyticus had high homology with other Vibrio species,and the highest homology with V.alginolyticus.In the prediction of functional domains,MsrA had the function of methionine sulfoxide reduction.In secondary structure prediction,MsrA contained random coils at a proportion of 46.70%,which was the highest.The similarity between the tertiary structure model of MsrA and template Q87SW6.1.A was 89.15%.PTM analysis showed that MsrA protein had many PTM modification sites such as phosphorylation and glycosylation sites.[Conclusions]This study provides some reference value for further study on the role of MsrA in bacterial antioxidant stress.

Cloning and Functional Validation of Mung Bean VrPR Gene

For the purpose of functional validation,the mung bean(Vigna radiata)VrPR gene was cloned and overexpressed in Arabidopsis thaliana.Thefindings revealed that the ORF of VrPR contained 1200 bp,in which 399 amino acids were encoded.Bioinformatics analysis showed that the VrPR protein belonged to the NADB Rossmann superfamily,which was one of the non-transmembrane hydrophilic proteins.VrPR was assumed to have 44 amino acid phosphorylation sites and be contained in chloroplasts.The VrPR secondary structure comprised of random coil,αhelix,βangle,and extended chain,all of which were quite compatible with the anticipated tertiary structure.Moreover,analysis of the phylogenetic tree indicated that the soybean PR(Glyma.12G222200)and VrPR were closely related.Furthermore,chlorophyll content in leaves is markedly increased in Arabidopsis when VrPR is overexpressed.Ourfindings will serve as a reference for more functional studies on the PR genes in mung bean.

Construction and Secretory Expression of β-Galactosidase Gene from Lactobacillus Bulgaricus in Lactococcus Lactis

Objective This study is to examine the secretion effects of β-galactosidase in Lactococcus lactis.Methods The usp45 and β-galactosidase genes were cloned and inserted into plasmid pMG36e to obtain the recombinant plasmid pMG36e-usp-lacZ.This recombinant plasmid was transformed into both Escherichia coli DH5α and L.lactis MG1363.The enzyme activity,gene sequencing,SDS-PAGE and hereditary stability were assessed and studied.Results The lacZ gene inserted into plasmids pMG36e-usp-lacZ was 99.37% similar to the GenBank sequence,and SDS-PAGE revealed an evident idio-strap at 116 KDa between L.lactis MG1363/pMG36eusp-lacZ in both supernatant and cell samples.β-Galactosidase activity measured 0.225 U/mL in L.lactis pMG36e-usp-lacZ transformants,and its secretion rate was 10%.The plasmid pMG36e-usp-lacZ appeared more stable in MG1363.Conclusion The authors concluded that these new recombinant bacteria well expressed and secreted β-galactosidase,indicating that the β-galactosidase expression system was successfully constructed,and this might provide a new solution for management of lactose intolerance specifically and promote the use of gene-modified organisms as part of the food-grade plasmid in general.

Cloning and Bioinformatics Analysis of vscB Gene of T3SS Chaperone of Vibrio alginolyticus

[Objectives]To clone and analyze the vscB gene of Vibrio alginolyticus HY9901 by bioinformatics.[Methods]A pair of specific primers were designed according to the vscB gene sequence of Vibrio alginolyticus HY9901.The full length of the primers was cloned by PCR and analyzed by bioinformatics.[Results]The vscB gene was 429 bp long,encoding 142 amino acids,with a theoretical molecular weight of 16.4 kDa and a pI value of 5.48.Amino acid sequence analysis of VscB showed that VscB was not a secretory protein,without signal peptide and transmembrane region,and there were protein kinase C phosphorylation site and casein kinase II phosphorylation site in the sequence.Homologous comparison of amino acid sequences showed that VscB of V.alginolyticus had the highest protein similarity with Vibrio Parahaemolyticus,reaching 91%.Phylogenetic tree analysis showed that the corresponding proteins of V.alginolyticus VscB,Vibrio Parahaemolyticus and Vibrio diabolicus were clustered in the same subfamily.Functional domain analysis showed that it had CesT family domain.Tertiary structure prediction showed that there were 3α-helices and 5β-turns in VscB protein.[Conclusions]This study provided a theoretical basis for further study on the function of chaperone of V.alginolyticus.

Construction of a Food-Grade Expression Vector Based on pMG36e by Using an α-Galactosidase Gene as a Selectable Marker

Construction of a food-grade expression vector for application to lactic acid bacteria（LAB） is of importance for dairy fermentation system. An α-galactosidase（aga） gene encoding an enzyme degrading melibiose was amplified by PCR from the plasmid p RAF800 of Lactococcus lactis NZ9000. The aga gene was introduced into pMG36 e to substitute the p rimary antibiotic selectable marker of pMG36 e, resulting in construction of a new food-grade expression vector pMG36-aga. To testify the expression efficiency of exogenous gene in pMG36-aga, a 1.5 kb long α-amylase（amy） gene from Ba cillus li cheniformis was cloned by PCR and introduced into the plasmid pMG36-aga. The resultant plasimd pMG36-aga-amy was transformed into L. lactis ML23 by electroporation. The positive clones were selected with the medium containing melibiose as the sole carbon source. Th e selection efficiency of aga was 8.71×103 CFU with a standard deviation of 9.1×102 CFU ？g-1 DNA of pMG36-aga. Furthermore, the SDS-PAGE analysis showed that the pMG36-aga-amy expressed a 56.4 kDa protein which was the same as the putati ve molecular weight of α-amylase. The starch plate assay also indicated that L. lactis ML23 displayed high activity of α-amylase by expressing of amy gene of pMG36-aga-amy.

Defect in an immune regulator gene BrSRFR1 leads to premature leaf senescence in Chinese cabbage

Leaf senescence is the final stage of leaf development, where the nutrients and energy of senescent leaves are redistributed to developing tissues or organs for plant growth, reproduction, and defense. Outer leaves are photosynthetic organs that usually senesce at the late heading stage in Chinese cabbage, and premature leaf senescence often reduces leafy head yield and quality. In this study, 11 premature leaf senescence mutants were screened from an ethyl methanesulfonate-mutagenized population of the double haploid line ‘FT' in Chinese cabbage. At the early heading stage, the mutants exhibited edge yellowing within its outer leaves, and at the mature stage, its leafy head weight decreased significantly. Genetic analysis revealed that the mutated trait of all 11 mutants corresponds to single gene recessive inheritance. Semi-diallel cross tests showed that 5 of the 11 were allelic mutants. MutMap and Kompetitive Allele Specific PCR genotyping revealed that BraA01g001400.3C was the candidate gene, which is orthologous of Arabidopsis SUPPRESSOR OF rps4-RLD 1, encoding an immune regulator, so we named it as BrSRFR1. All the BrSRFR1 in the five allelic mutants exhibited single nucleotide polymorphisms at different positions on their exons and led to premature translation termination, which confirmed that defect in BrSRFR1 led to premature leaf senescence. These results verify the role of Br SRFR1 on leaf senescence and provide a new insight into the mechanisms of leaf senescence in Chinese cabbage, which reveals a novel function of SRFR1 in plant development.

EXPRESSION OF HUMAN α-GALACTOSIDASE AND α1,2-FUCOSYL-TRANSFERASE GENES MODIFIES THE CELL SURFACE GALα1,3GAL ANTIGEN AND CONFERS RESISTANCETO HUMAN SERUM-MEDIATED CYTOLYSIS

Objective To explore the strategies which reduce the amount of xenoantigen Galα1, 3 Gal. Methods Human α-galactosidase gene and α1,2-fucosyltransferase gene were transferred into cul-tured porcine vascular endothelial cells PEDSV.15 and human α-galactosidase transgenic mice were produced. The Galα1,3Gal on the cell surface and susceptibility of cells to human antibody-mediated lysis were analyzed. Results Human α-galactosidase gene alone reduced 78% of Galα1,3Gal on PEDSV.15 cell surface while human α-galactosidase combined with α1,2-fucosyltransferase genes removed Galα1,3Gal completely. Decrease of Galα1,3Gal could reduce susceptibility of cells to human antibody-mediated lysis, especially during co-expression of α-galactosidase gene and α1,2-fucosyltransferase gene. RT-PCR indicated positive human α-galactosidase gene expression in all organs of positive human α-galacto-sidase transgenic F1 mice including heart, liver, kidney, lung, and spleen, the amount of Galα1,3Gal antigens on which was reduced largely. 58% of spleen cells from F1 mice were destroyed by comp-lement-mediated lysis compared with 24% of those from normal mice. Conclusions Human α-galactosidase gene and α1,2-fucosyltransferase gene effectively reduce the expression of Galα1,3Gal antigens on endothelial cell surface and confers resistance to human serum-mediated cytolysis. The expression of human α-galactosidase in mice can also eliminate the Galα1,3Gal antigens in most tissues and decrease the susceptibility of spleen cells to human serum-mediated cytolysis.

Cloning and Prokaryotic Expression of FnBP Ligand Binding Gene of Staphylococcus aureus

[Objective] The study aimed to clone the FnBP ligand binding gene of Staphylococcus aureus and run prokaryotic expression by constructing a prokaryotic expression vector. [Method] The gene encoding FnBP ligand binding gene was amplified from S.aureus chromosomal DNA by PCR technique. After T-A cloning, plasmid pMD18- FnBP was constructed. pMD18- FnBP and pET28a(+)were digested by BamH Ⅰ and EcoR Ⅰ double enzymes, then the purified FnBP ligand binding gene was subcloned into the expression vector pET28a(+), and the prokaryotic expression vector pET28a-FnBP was thus constructed. The constructed plasmid pET28a-FnBP was transformed into Escherichia coli BL21(DE3) competent cells. The bacterium was induced by IPTG and the expressed products were analyzed by SDS-PAGE and Western blot. [Result] The gene fragment with the length of 370 bp was amplified by PCR approach. One approximately 30 kD exogenous protein was observed in SDS-PAGE analysis. Western blot analysis indicates the protein has antigenicity of S.aureus. [Conclusion] The FnBP ligand binding gene of S.aureus was successfully cloned and expressed in prokaryotic cells.

Cloning and Bioinformatics Analysis of P23 Gene from Theileria sergenti

[Objective] The aim of this study is to provide basis for developing genetic engineering vaccine and diagnostic kit for Theileria sergenti infection. [Objective] P23 gene of Theileria sergenti was amplified from its genomic DNA by PCR amplification, and cloned into the pGEM-Easy vector; then the sequencing result was analyzed with bioinformatics methods. [Result] Whole length of the P23 gene from Theileria sergenti is 684 bp containing a 672 bp open reading frame. The deduced amino acid sequence (223 amino acid residues) contains a signal peptide of 19 amino acid residues and two fragments of transmembrane domains, with relative molecular weight of the 25.886 kD and with the pI of 9.22. The homology between the yielded sequence and Chitose of Theileria sergenti P23 gene(TS-Chitose type, D84446), Ikeda of Theileria sergenti P23 gene(TS-Ikeda type, D84447) reached 99% and 90%, respectively. The sequence has been accessed in GenBank(EU573168). [Conclusion] The protein encoded by the P23 gene has better stability and immunogenicity, thus can be used as the antigen candidate for preparing genetic engineering vaccine for Theileria sergenti.

Cloning of Potato POTHR-1 Gene and Its Expression in Response to Infection by Phytophthora infestans and Other Abiotic Stimuli

A complete cDNA of potato Phytophthora infestans-induced hypersensitive response-related protein gene (POTHR-1) was cloned using rapid amplification of cDNA ends (RACE) strategy according to a fragment sequence which we had cloned using suppression subtractive hybridization (SSH) technique. The potato POTHR-1 gene encodes a protein of 225 amino acids, which shares 81% identity with tobacco hin1 gene-enoded protein (harpin-induced protein). Southern blot revealed that there are two to three copies of POTHR-1 in potato genome. The POTHR-1 gene expression in potato leaves showed that its transcripts accumulated remarkably in leaves after 36 h inoculation with P. infestans. Mechanical wounding and jasmonic acid (JA) could induce the POTHR-1 gene expression and osmotic stress just induce a slight accumulation of POTHR-1 gene mRNA, while salicylic acid (SA) had no detectable function on the induction accumulation of POTHR-1 gene transcripts. The potato POTHR-1 gene may preferentially associate with hypersensitive response (HR) or biotic cell death during interaction between host and pathogen.

Molecular Cloning and Expression of RSSG58 Gene in Rice Sperm Cells

Myosins, a large family of structurally diverse mechanoenzymes, which, upon interaction with actin filaments, convert energy from ATP hydrolysis into mechanical force, play an important role in male reproductive processes. In this study we report the rice ( Oryza sativa L.) RSSG58 gene, which was cloned from the cDNA library of rice sperm cells by using sperm cell mainly expression subtractive clone as probe. This gene encodes a putative 66.7 W polypeptide, which shows similarity to the myosin heavy chain of Arabidopsis thaliana, and consists of 579 amino acids with an isoelectric point (pI) of 4.885. RSSG58, which is a member of a divergent gene family, generates transcripts of 2 278 bp and 2 437 bp that differ only in their polyadenylation sites. Southern hybridization showed that RSSG58 has only one copy in rice genome and RSSG58 transcripts are most abundant in sperm cells, with two distinct signals. The RT-PCR analysis indicated that the transcriptions of the RSSG58 gene were various in the different development stages and tissues. The greatest accumulation of RSSG58 mRNA was detected in sperm cells, while weaker expression was detected in leaves, microspore mother cells, unicellular microspore pollen stage, two-cell stage pollens, mature pollens and pollinated ovaries. These results suggest that RSSG58 is especially abundantly expressed in rice sperm cells.

Cloning of Rabbit Bone Morphogenetic Protein 15 and Its Expression During in vitro Maturation of Rabbit Oocytes

Partial cDNA sequence of rabbit BMP15 was cloned by RT-PCR from rabbit ovaries, showing a similarity of 83%-90% with the BMP15 nucleotide sequences in humans, mice, ovine, sheep, cows and pigs. The expression of BMP15 in rabbit cumulus-oocyte complexs during oocytes in vitro maturation （IVM） was measured by fluorescent quantitative RT-PCR method. BMP 15 was expressed at low levels in immature oocytes and increased to the highest level at 16h of IVM, which coincides with the time of cumulus cell expansion, then declined slowly under IVM cultivation. The expression pattern of BMP 15 suggested that it might be important in cumulus expansion in rabbits.

Cloning and Analysis of WRKY Gene of Rice Induced by Rhizoctonia solani Kuhn

[Objective] The aim was to clone the up-regulated expression gene of rice induced by Rhizoctonia solani.[Method] The EST fragment K16 obtained by suppression subtraction hybridization（SSH）was cloned and confirmed by reverse transcription-polymerase chain reaction（RT-PCR）.Then RT-PCR products were cloned into the PMD18-T vector and sequenced.The functions of the sequence were predicted with bioinformatics method.[Result] A 1 079 bp gene was obtained.The gene encoded a protein with 236 amino acids.The protein contains many motif sites,two WRKY domains and a C2H2 zinc finger motif.The gene showed high identities with WRKY8,WRKY24 and WRKY30 gene of rice.[Conclusion] The up-regulated expression gene induced by R.solani was representative WRKY family gene.The gene could play an important role on rice sheath blight resistance.

Cloning, Characterization and Chromosome Localization of Two Powdery Mildew Resistance-Related Gene Sequences from Wheat

Reverse_transcription Polymerase Chain Reaction (RT_PCR) was performed using cDNAs as templates from wheat_ Haynaldia villosa 6VS/6AL translocation line and 'Yangmai 5' induced with fungus Erysiphe graminis , and degenerate primers designed based on the conserved amino acid sequences of known plant disease_resistance genes. The cDNA sequences encoding cyclophilin_like and H +_ATPase_like genes were first isolated and characterized in wheat. The putative amino acid sequences of the two clones showed that they were highly homologous to those of cyclophilin proteins and H +_ATPases isolated from other plants. Thus they were designated as Ta_Cyp and Ta_MAH . The obvious expression differences could be observed between wheat_ H. villosa 6VS/6AL translocation line and susceptible wheat cultivar 'Yangmai 5', implying that the two genes may be related with the resistance of wheat_ H. villosa 6VS/6AL translocation line to disease. Southern blot indicated that the wheat genome contained 2-3 copies of Ta_Cyp gene and one copy of the Ta_MAH gene. Chinese Spring nulli_tetrasomic line analysis located the Ta_Cyp homologous genes on wheat chromosome 6A, 6B and 6D. Southern blot using Ta_Cyp clone as a probe showed that the polymorphic bands existed among the H. villosa , amphiploid of Triticum durum _ H. villosa , wheat_ H. villosa 6VS/6AL translocation line and 'Yangmai 5', suggesting that Ta_Cyp homologies exist in wheat genome as well as on the short arm of chromosome 6V in H. villosa .

Molecular Cloning and Sequence Analysis of Class Ⅱ Chitinase Gene in Leymus chinensis

[ Objective] The aim of this study was to clone Class Ⅱ chitinase gene in Leymus chinensis grown in saline land in Heilongjiang Province and analyze its sequence, which provided a foundation for further study on the biological function and application of chitinasa gene. [ Method] cDNA library of Leymus chinensis leaves were constructed, and its DNA sequence was determined or analyzed, while the homology of chitinasa gene and amino acid sequence was compared with that in GenBank. [ Result] One full length cDNA fragment with length of 996 bp was cloned from cDNA library of Leymus chinensis leaves. The length of ORF was 768 bp encoding 225 amino acids （GenBank accession number： EU344908）. The encoding products lacked CBD and C-terminal extension region from the view of structure, but had structural characters of Class Ⅱ chitinase gene, which indicated that amino acid sequence had high homology compared with Class Ⅱ chitinase gene of rye and wheat. The constructed recombinant vector pQE-LcChi2 could express a protein of 27 kD through induction, which was consistent with the deduced encoding product of pQE-LcChi2 gene. [ Conclusion] LcChl2 gene is an expression gene, which can express in E. coll.

Cloning and Bioinformatics Analysis of Lfcin Gene of Datong Yak

[Objective] This study was to clone Lfcin gene from Datong yak, so as to provide reference for applying this gene in feed industry and breeding industry. [Method] Using PCR technology, the lactoferricin(Lfcin)-encoding gene was obtained from genome of Datong yak; then it was cloned into pGEM-T easy vector, and then sequenced; the sequencing results were subsequently aligned with the sequences of dairy cow accessed in GenBank. Moreover, amino acid sequences of Lfcin gene from various species including yak, dairy cow, human and mouse were used for sequence alignment and phylogenesis analysis. [Result] The second exon of lactoferrin(LF) from Datong yak, which is 778 bp in length, was obtained, within which the coding region of Lfcin gene is 75 bp (25 amino acid residues); sequence analysis showed that there is discrepancy of eleven bases between Datong yak and dairy cow; Lfcin proteins from various species shared high homeology, of which that from Datong yak and dairy cow were completely identical; phylogenesis analysis showed that cladogram based on Lfcin was consistent with species evolutionary law. [Conclusion] This study laid a foundation for the prokaryotic or eukaryotic expression of Lfcin gene and further understanding the activity of Lfcin protein.

Study on the Cloning and Isolation of sus scrofa GPX2 Gene by RACE Method

[Objective] Using molecular biotechnology to clone the sus scrofa GPX2 gene. [Method] Using total RNA of sus scrofa duodenum as template, degenerated primer pairs were designed according to the homology alignment analysis of GPX2 gene of human, rat, mouse, dog and cattle. A sus scrofa GPX2 gene sequence of 330 bp was obtained by RT-PCR application method. Primes were designed respectively according to the known sequence, sus scrofa GPX2 gene was isolated and cloned by 3-RACE and 5-RACE method and analyzed the gene sequence. [Result] A mRNA sequence of 924 bp was successfully cloned and isolated in this research. This sequence contained complete 3＇end and had higher sequence homology with human,mouse,cattle and dog GPX2 gene, and there was codon called TGA which encoding Sec on the position of No. 114-116 gene. [Conclusion] Sequence alignment analysis showed that the cloned gene was sus scrofa GPX2 gene （ NCBI GenBank database, the sequence number was D098982）.