Modifying the amino acids in conformational motion pathway of the α-amylase of Geobacillus stearothermophilus improved its activity and stability.

Amino acids along the conformational motion pathway of the enzyme molecule correlated to its flexibility and rigidity. To enhance the enzyme activity and thermal stability, the motion pathway of Geobacillus stearothermophilus α-amylase has been identified and molecularly modified by using the neural relational inference model and deep learning tool. The significant differences in substrate specificity, enzymatic kinetics, optimal temperature, and thermal stability were observed among the mutants with modified amino acids along the pathway. Mutants especially the P44E demonstrated enhanced hydrolytic activity and catalytic efficiency (kcat/KM) than the wild-type enzyme to 95.0% and 93.8% respectively, with the optimum temperature increased to 90°C. This mutation from proline to glutamic acid has increased the number and the radius of the bottleneck of the channels, which might facilitate transporting large starch substrates into the enzyme. The mutation could also optimize the hydrogen bonding network of the catalytic center, and diminish the spatial hindering to the substrate entry and exit from the catalytic center.

The Motion Paradigm of Pre-Dock Zearalenone Hydrolase Predictions with Molecular Dynamics and the Docking Phase with Umbrella Sampling.

Zearalenone (ZEN) is one of the most prevalent estrogenic mycotoxins, is produced mainly by the Fusarium family of fungi, and poses a risk to the health of animals. Zearalenone hydrolase (ZHD) is an important enzyme capable of degrading ZEN into a non-toxic compound. Although previous research has investigated the catalytic mechanism of ZHD, information on its dynamic interaction with ZEN remains unknown. This study aimed to develop a pipeline for identifying the allosteric pathway of ZHD. Using an identity analysis, we identified hub genes whose sequences can generalize a set of sequences in a protein family. We then utilized a neural relational inference (NRI) model to identify the allosteric pathway of the protein throughout the entire molecular dynamics simulation. The production run lasted 1 microsecond, and we analyzed residues 139-222 for the allosteric pathway using the NRI model. We found that the cap domain of the protein opened up during catalysis, resembling a hemostatic tape. We used umbrella sampling to simulate the dynamic docking phase of the ligand-protein complex and found that the protein took on a square sandwich shape. Our energy analysis, using both molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MMPBSA) and Potential Mean Force (PMF) analysis, showed discrepancies, with scores of -8.45 kcal/mol and -1.95 kcal/mol, respectively. MMPBSA, however, obtained a similar score to that of a previous report.

High-level secretive expression of a novel achieved Talaromyces cellulolyticus endo-polygalacturonase in Pichia pastoris by improving gene dosage for hydrolysis of natural pectin.

Pectin is a type of complex hydrophilic polysaccharide widely distributed in plant resources. Thermal stable pectinase has its advantage in bioapplication in the fields of food processing, brewing, and papermaking, etc. In this study, we enzymatically characterized a putative endo-polygalacturonase TcPG from a Talaromyces cellulolyticus, realized its high-level expression in Pichia pastoris by in vitro constructing of a series of multi-copy expression cassettes and real time quantitative PCR screening. The secretive expression level of TcPG was nonlinear correlated to the gene dosage. Recombinants with five-copy TcPG gene in the host genome showed the highest expression. After cultivation in a bioreactor for about 96 h, the enzyme activity reached 7124.8 U/mL culture. TcPG has its optimal temperature of 70 °C. Under the optimized parameters, the pectin could be efficiently hydrolyzed into oligosaccharides.

Gene dosage and coexpression with endoplasmic reticulum secretion-associated factors improved the secretory expression of α-galactosidase.

The α-galactosidases, which can catalyze the removal of α-1,6-linked terminal galactose residues from galactooligosaccharide materials, have good potential for industrial applications. The high-level and efficient secretion of the α-galactosidases into the extracellular space has greatly simplified the downstream bioengineering process, facilitating their bioapplications. In this study, the effects of gene dosage and endoplasmic reticulum secretion-associated factors (ERSAs) on the secretory expression of an α-galactosidase gene derived from a Aspergillus oryzae strain were investigated by constructing multicopy expression cassettes and coexpressing the α-galactosidase gene with ERSAs. With the increase in the gene copy-number in the host genome, the expression of GalA was improved. However, the secretory expression level was not linearly related to the copy number. When the number was higher than four copies, the expression level of GalA gene declined. The ERSAs factors HAC1, PDI, and Ero1 improved the secretory expression of α-galactosidase, while Hsp40 inhibited its secretion. After methanol-induced expression in a bench-top bioreactor, Pichia recombinants carrying four copies of GalA genes reached 3520 U/mL in the supernatant of the culture. We further optimized the parameters for α-galactosidase to hydrolyze two types of galactooligosaccharides: raffinose and stachyose. This study has fulfilled the scale-up production of α-galactosidase, thus facilitating its industrial applications.

Vertical stratification of bacteria and the chemical compounds in crude oil-contaminated soil layers of the semi-deserted Dzungharian Basin.

The largely semi-deserted and deserted Dzungharian Basin sites in the northwest of China geologically represent an extension of the Paleozoic Kazakhstan Block and were once part of an independent continent. For reasons of overdevelopment and unreasonable operations during the process of exploitation and transportation, oil pollutants that were discharged into the soil environment caused serious pollution in this weak ecosystem. To explore the bacterial community composition in detail and their possible origination and potential during the natural attenuation of petroleum contaminants in this type of ecologic niche, GC-MS and high-throughput sequencing techniques were used to resolve the organic compounds and bacterial communities in vertical soil layers. The degradation of petroleum contaminants in semi-deserted and deserted soils mainly occurred in the layer at a depth of 45-55 cm. During this process, aromatic and heterocyclic compounds were significantly enriched in soils. The bacterial communities in this basin exhibited a distinct vertical stratification from the surface layer down to the bottom soil layer. Considering the interaction between the community composition and the geochemical properties, we speculate that the degradation of petroleum contaminants in this semi-deserted and deserted soil might represent a microorganism-mediated process and mainly occur in the deeper soil layer.

Aromatic amino acids in the cellulose binding domain of Penicillium crustosum endoglucanase EGL1 differentially contribute to the cellulose affinity of the enzyme.

The cellulose binding domain (CBD) of cellulase binding to cellulosic materials is the initiation of a synergistic action on the enzymatic hydrolysis of the most abundant renewable biomass resources in nature. The binding of the CBD domain to cellulosic substrates generally relies on the interaction between the aromatic amino acids structurally located on the flat face of the CBD domain and the glucose rings of cellulose. In this study, we found the CBD domain of a newly cloned Penicillium crustosum endoglucanase EGL1, which was phylogenetically related to Aspergillus, Fusarium and Rhizopus, and divergent from the well-characterized Trichoderma reeseis cellulase CBD domain, contain two conserved aromatic amino acid-rich regions, Y451-Y452 and Y477-Y478-Y479, among which three amino acids Y451, Y477, and Y478 structurally sited on a flat face of this domain. Cellulose binding assays with green fluorescence protein as the marker, adsorption isotherm assays and an isothermal titration calorimetry assays revealed that although these three amino acids participated in this process, the Y451-Y452 appears to contribute more to the cellulose binding than Y477-Y478-Y479. Further glycine scanning mutagenesis and structural modelling revealed that the binding between CBD domain and cellulosic materials might be multi-amino-acids that participated in this process. The flexible poly-glucose molecule could contact Y451, Y477, and Y478 which form the contacting flat face of CBD domain as the typical model, some other amino acids in or outside the flat face might also participate in the interaction. Thus, it is possible that the conserved Y451-Y452 of CBD might have a higher chance of contacting the cellulosic substrates, contributing more to the affinity of CBD than the other amino acids.

The catalytic domain of Penicillium crustosum endoglucanase EGL1 has cellulose-binding capacity and cellulolytic activity.

The cellulase-mediated degradation of cellulosic materials, which is initiated by endoglucanases by the random cleavage of the glycosidic bonds between glucose units to break long cellulose molecules into shorter ones, represents a major carbon flow in the global carbon cycle. The structure of a typical endoglucanase contains a classical (α/ß)8 barrel fold catalytic domain, a linker region and a cellulose-binding domain. In this study, we found that both the full-length enzyme and the catalytic domain of endoglucanase EGL1 cloned from Penicillium crustosum strain 601 have CMCase and FPase activity. A cellulose-binding assay using green fluorescent protein as a marker further showed that the catalytic domain could also bind the cellulose substrate. The three-dimensional structure of the catalytic domain of EGL1 revealed that this cellulose substrate-binding capacity of the catalytic domain may come from the hydrophobic core formed by aromatic amino acids distributed in or outside the (α/ß)8 barrel fold. A glycine scanning mutagenesis assay further found that the aromatic amino acids at the bottom of the barrel fold and those adjacent to the catalytic site significantly affect the cellulolytic activity and the cellulose binding affinity of the catalytic domain. Thus, it could be speculated that the aromatic amino acids in the bottom of the barrel fold might be the main contributors in the binding capacity of the catalytic domain with the cellulose substrate, and those distributed around the active sites on the top of the enzyme might participate in moving the cellulose substrate to the active site in the barrel fold or releasing the hydrolysis products.

Constitutive expression of Campylobacter jejuni truncated hemoglobin CtrHb improves the growth of Escherichia coli cell under aerobic and anaerobic conditions.

Bacteria hemoglobin could bind to the oxygen, transfer it from the intracellular microenvironment to the respiration process and sustain the energy for the metabolism and reproduction of cells. Heterologous expression of bacteria hemoglobin gene could improve the capacity of the host on oxygen-capturing and allow it to grow even under microaerophilic condition. To develop a system based on hemoglobin to help bacteria cells overcome the oxygen shortage in fermentation, in this study, Campylobacter jejuni truncated hemoglobin (CtrHb) gene was synthesized and expressed under the control of constitutive expression promoters P2 and P(SPO1-II) in Escherichia coli. As showed by the growth curves of the two recombinants P2-CtrHb and P(SPO1-II)-CtrHb, constitutive expression of CtrHb improved cell growth under aerobic shaking-flasks, anaerobic capped-bottles and bioreactor conditions. According to the NMR analysis, this improvement might come from the expression of hemoglobin which could boost the metabolism of cells by supplying more oxygen to the respiratory chain processes. Through semi-quantitative RT-PCR and CO differential spectrum assays, we further discussed the connection between the growth patterns of the recombinants, the expression level of CtrHb and oxygen binding capacity of CtrHb in cells. Based on the growth patterns of these recombinants in bioreactor, a possible choice on different type of recombinants under specific fermentation conditions was also suggested in this study.

Recombinant Escherichia coli strains with inducible Campylobacter jejuni single domain hemoglobin CHb expression exhibited improved cell growth in bioreactor culture.

Maintaining an appropriate concentration of dissolved oxygen in aqueous solution is critical for efficient operation of a bioreactor, requiring sophisticated engineering design and a system of regulation to maximize oxygen transfer from the injected air bubbles to the cells. Bacterial hemoglobins are oxygen-binding proteins that transfer oxygen from the environment to metabolic processes and allow bacteria to grow even under microaerophilic conditions. To improve the oxygen utilization efficiency of cells and overcome the oxygen shortage in bioreactors, the gene coding for the Campylobacter jejuni single domain hemoglobin (CHb) gene was artificially synthesized and functionally expressed under the control of inducible expression promoters PT7 and Pvgh in Escherichia coli. The effects of the recombinants PT7-CHb and Pvgh-CHb on cell growth were evaluated in aerobic shake flasks, anaerobic capped bottles and a 5-L bioreactor, and a pronounced improvement in cell biomass was observed for CHb-expressing cells. To determine the growth curves, CHb gene expression, and CHb oxygen-binding capacity of specific recombinants with different promoters, we determined the time course of CHb gene expression in the two recombinants by semi-quantitative RT-PCR and CO differential spectrum assays. Based on the growth patterns of the two recombinants in the bioreactor, we proposed different recombinant types with optimal performance under specific culture conditions.

Codon optimization, promoter and expression system selection that achieved high-level production of Yarrowia lipolytica lipase in Pichia pastoris.

Lipase (EC 3.1.1.3) stands amongst the most important and promising biocatalysts for industrial applications. In this study, in order to realize a high-level expression of the Yarrowia lipolytica lipase gene in Pichia pastoris, we optimized the codon of LIP2 by de novo gene design and synthesis, which significantly improved the lipase expression when compared to the native lip2 gene. We also comparatively analyzed the effects of the promoter types (PAOX1 and PFLD1) and the Pichia expression systems, including the newly developed PichiaPink system, on lipase production and obtained the optimal recombinants. Bench-top scale fermentation studies indicated that the recombinant carrying the codon-optimized lipase gene syn-lip under the control of promoter PAOX1 has a significantly higher lipase production capacity in the fermenter than other types of recombinants. After undergoing methanol inducible expression for 96h, the wet cell weight of Pichia, the lipase activity and the protein content in the fermentation broth reached their highest values of 262g/L, 38,500U/mL and 2.82g/L, respectively. This study has not only greatly facilitated the bioapplication of lipase in industrial fields but the strategies utilized, such as de novo gene design and synthesis, the comparative analysis among promoters and different generations of Pichia expression systems will also be useful as references for future work in this field.

Community composition and cellulase activity of cellulolytic bacteria from forest soils planted with broad-leaved deciduous and evergreen trees.

Cellulolytic bacteria in forest soil provide carbon sources to improve the soil fertility and sustain the nutrient balance of the forest ecological system through the decomposition of cellulosic remains. These bacteria can also be utilized for the biological conversion of biomass into renewable biofuels. In this study, the community compositions and activities of cellulolytic bacteria in the soils of forests planted with broad-leaved deciduous (Chang Qing Garden, CQG) and broad-leaved evergreen (Forest Park, FP) trees in Wuhan, China were resolved through restriction fragment length polymorphism (RFLP) and sequencing analysis of the 16S rRNA gene. All of the isolates exhibited 35 RFLP fingerprint patterns and were clustered into six groups at a similarity level of 50 %. The phylogeny analysis based on the 16S rRNA gene sequence revealed that these RFLP groups could be clustered into three phylogenetic groups and further divided into six subgroups at a higher resolution. Group I consists of isolates from Bacillus cereus, Bacillus subtilis complex (I-A) and from Paenibacillus amylolyticus-related complex (I-B) and exhibited the highest cellulase activity among all of the cellulolytic bacteria isolates. Cluster II consists of isolates belonging to Microbacterium testaceum (II-A), Chryseobacterium indoltheticum (II-B), and Flavobacterium pectinovorum and the related complex (II-C). Cluster III consists of isolates belonging to Pseudomonas putida-related species. The community shift with respect to the plant species and the soil properties was evidenced by the phylogenetic composition of the communities. Groups I-A and I-B, which account for 36.0 % of the cellulolytic communities in the CQG site, are the dominant groups (88.4 %) in the FP site. Alternatively, the ratio of the bacteria belonging to group III (P. putida-related isolates) shifted from 28.0 % in CQG to 4.0 % in FP. The soil nutrient analysis revealed that the CQG site planted with deciduous broad-leaved trees has a richer organic nutrient (total organic carbon and total nitrogen) than the FP site planted with evergreen broad-leaved trees. Against this background, the population density and the diversity of cellulolytic bacteria in the CQG site are clearly higher than those in the FP site, and the latter was dominated with high-cellulase-activity Bacillus- and Paenibacillus-related bacteria. The canonical correspondence analysis further indicated that the distribution of these groups is correlated with the FP site, whereas groups II and III are correlated with the organic nutrient-rich CQG site.

Community composition of nirS-type denitrifier in a shallow eutrophic lake.

Denitrification is a major biological process to reduce nitrate to molecular nitrogen (N2). In shallow eutrophic lakes, this process can remove the largest portion of fixed nitrogen and plays an important role in self-purification of this ecosystem. To understand the structure of denitrifying communities in a shallow eutrophic lake, denitrifier communities in four sub-lakes of East Lake in Wuhan, China, were explored by restriction fragment length polymorphisms (RFLP) analysis and sequencing of nirS gene clone libraries. nirS is a functional marker gene for denitrification encoding cytochrome cd 1-containing nitrite reductase, which catalyzes the reduction of nitrite to nitric oxide. Both RFLP fingerprints clustering analysis and phylogeny analysis based on the amino acid sequences of NirS revealed that NirS-type communities in East Lake sediment could be roughly divided into three clusters. Cluster I accounted for 74-82 % of clones from the moderately eutrophic sub-lakes Tuan, Tang Ling, and Guo Zheng. Cluster II accounted for 76 % of the communities in hypertrophic sub-lake Miao Lake and cluster III as a minor group (7 % of the total), mainly presented in Miao Lake. Phylogenetic analysis revealed that cluster I was related to the reference clones from a broad range of ecological environments, and clusters II and III were more phylogenetically related to the reference clones from entrophic environments. Canonical correspondence analysis indicated that total nitrogen, total phosphate, total organic carbon, and NH4-N and NO2-N were important environmental factors affecting the dispersion of NirS-type denitrifier in the sediments. Cluster I showed a weak relationship with the nutrient content, while cluster II and III were positively related with the nutrient content. Principal coordinates analysis indicated that NirS-type communities from Tuan Lake, Tang Ling Lake, and Guo Zheng Lake sediments were divergent from those found in river, estuary sediment, and forest soil but similar to communities in constructed wetland sediment despite large geographic distances. The communities from the hypertrophic sub-lake Miao Lake deviated from other sub-lakes and the reference communities and clustered independently. Our results support the argument that environmental factors regulate the composition and distribution of the functional bacterial groups.

Artifactual pyrosequencing reads in multiple-displacement-amplified sediment metagenomes from the Red Sea.

The Multiple Displacement Amplification (MDA) protocol is reported to introduce different artifacts into DNA samples with impurities. In this study, we report an artifactual effect of MDA with sediment DNA samples from a deep-sea brine basin in the Red Sea. In the metagenomes, we showed the presence of abundant artifactual 454 pyrosequencing reads over sizes of 50 to 220 bp. Gene fragments translocated from neighboring gene regions were identified in these reads. Occasionally, the translocation occurred between the gene fragments from different species. Reads containing these gene fragments could form a strong stem-loop structure. More than 60% of the artifactual reads could fit the structural models. MDA amplification is probably responsible for the massive generation of the artifactual reads with the secondary structure in the metagenomes. Possible sources of the translocations and structures are discussed.

de novo design and synthesis of Candida antarctica lipase B gene and α-factor leads to high-level expression in Pichia pastoris.

Candida antarctica lipase B (CALB) is one of the most widely used and studied enzymes in the world. In order to achieve the high-level expression of CALB in Pichia, we optimized the codons of CALB gene and α-factor by using a de novo design and synthesis strategy. Through comparative analysis of a series of recombinants with different expression components, we found that the methanol-inducible expression recombinant carrying the codon-optimized α-factor and mature CALB gene (pPIC9KαM-CalBM) has the highest lipase production capacity. After fermentation parameters optimization, the lipase activity and protein content of the recombinant pPIC9KαM-CalBM reached 6,100 U/mL and 3.0 g/L, respectively, in a 5-L fermentor. We believe this strategy could be of special interest due to its capacity to improve the expression level of target gene, and the Pichia transformants carrying the codon-optimized gene had great potential for the industrial-scale production of CALB lipase.

A simple and accurate two-step long DNA sequences synthesis strategy to improve heterologous gene expression in pichia.

In vitro gene chemical synthesis is a powerful tool to improve the expression of gene in heterologous system. In this study, a two-step gene synthesis strategy that combines an assembly PCR and an overlap extension PCR (AOE) was developed. In this strategy, the chemically synthesized oligonucleotides were assembled into several 200-500 bp fragments with 20-25 bp overlap at each end by assembly PCR, and then an overlap extension PCR was conducted to assemble all these fragments into a full length DNA sequence. Using this method, we de novo designed and optimized the codon of Rhizopus oryzae lipase gene ROL (810 bp) and Aspergillus niger phytase gene phyA (1404 bp). Compared with the original ROL gene and phyA gene, the codon-optimized genes expressed at a significantly higher level in yeasts after methanol induction. We believe this AOE method to be of special interest as it is simple, accurate and has no limitation with respect to the size of the gene to be synthesized. Combined with de novo design, this method allows the rapid synthesis of a gene optimized for expression in the system of choice and production of sufficient biological material for molecular characterization and biotechnological application.

Bacterial niche-specific genome expansion is coupled with highly frequent gene disruptions in deep-sea sediments.

The complexity and dynamics of microbial metagenomes may be evaluated by genome size, gene duplication and the disruption rate between lineages. In this study, we pyrosequenced the metagenomes of microbes obtained from the brine and sediment of a deep-sea brine pool in the Red Sea to explore the possible genomic adaptations of the microbes in response to environmental changes. The microbes from the brine and sediments (both surface and deep layers) of the Atlantis II Deep brine pool had similar communities whereas the effective genome size varied from 7.4 Mb in the brine to more than 9 Mb in the sediment. This genome expansion in the sediment samples was due to gene duplication as evidenced by enrichment of the homologs. The duplicated genes were highly disrupted, on average by 47.6% and 70% for the surface and deep layers of the Atlantis II Deep sediment samples, respectively. The disruptive effects appeared to be mainly due to point mutations and frameshifts. In contrast, the homologs from the Atlantis II Deep brine sample were highly conserved and they maintained relatively small copy numbers. Likely, the adaptation of the microbes in the sediments was coupled with pseudogenizations and possibly functional diversifications of the paralogs in the expanded genomes. The maintenance of the pseudogenes in the large genomes is discussed.

Homologous overexpression of a lipase from Burkholderia cepacia using the lambda Red recombinase system.

Red recombinase system of the lambda phage is widely used for recombination of short linear DNA fragments and genome. Using this system, we obtained T7 RNA polymerase (RNAP) substitution mutants in Burkholderia cepacia. To test the expression abilities of the T7 mutants, four different lipase expression vectors were transformed and the lipase activity of these recombinants was evaluated. Our results suggest that 500 nt homology between the unit and the genome is sufficient to generate mutations and this strategy enables the rapid establishment of mutant strains with efficiencies of 85%. After expression and purification, the highest purified lipase activity obtained was 3,990 U/l, nearly triple that of the wild-type organism.

Aspergillus niger lipase: gene cloning, over-expression in Escherichia coli and in vitro refolding.

From the N-terminal amino acid sequence of the lipase from Aspergillus niger F044, a potential homologous gene A84689 to the lipanl (the gene encoding the lipase from Aspergillus niger F044) was identified. A pair of primers were designed according to the nucleotide sequence of A84689, and the lipanl was cloned by PCR. Nucleotide sequencing revealed that the lipanl has an ORF of 1,044 bp, containing three introns. The deduced amino acid sequence corresponds to 297 amino acid residues. The cloned cDNA fragment encoding the mature lipase from Aspergillus niger F044 was over-expressed in Escherichia coli BL21(De3) and the recombinant protein was refolded in vitro by dilution followed by DEAE Sepharose Fast Flow chromatography.

[Purification and characterization of a lipase from Aspergillus niger F044].

A lipase from Aspergillus niger F044 was purified to homogeneity using ammonium sulfate precipitation, dialysis, DEAE-Sepharose Fast Flow anion exchange chromatography and Sephadex G-75 gel filtration chromatography. This purification protocol resulted in a 73.71-fold purification of lipase with 33.99% final yield, and the relative molecular weight of the enzyme was determined to be approximately 35-40kD using SDS-PAGE. The optimum pH and temperature for lipolytic activity of the lipase was 7.0 and 45 degrees C , respectively. It was extremely stable at 60 degrees C and retained 98.70% of its original activity for 30min. The stability declined rapidly as soon as the temperature rose over 65 degrees C . The lipase was highly stable in the pH range from 2.0 to 9.0 for 4h. Ca2+ and Mg2+ ions stimulated lipolytic activity, whereas Mn2+ , Fe2+ and Zn2+ ions caused inhibition. The values of Km and Vmax calculated from the Lineweaver-Burk plot using pNPP as hydrolysis substrate were 7.37mmol/L and 25.91 micromol/(min x mg), respectively. The N-terminal sequence of the lipase was Ser/Glu/His-Val-Ser-Thr-Ser-Thr-Leu-Asp-Glu-Leu-Gln-Leu-Phe-Ala-Gln, which is highly homogeneity with that of lipase, as reported by Torossian.

Genotypic characteristics of the rrn operon and genome of indigenous soybean Bradyrhizobia in cropping zones of China.

Four genetic assays, 16S rRNA restriction fragment length polymorphism (RFLP), 16S rRNA sequencing, 16S-23S rRNA intergenetic spacer (IGS) RFLP, and amplified fragment length polymorphism (AFLP), were conducted to determine the genotypic characteristics of 44 indigenous strains of Bradyrhizobium from soybean (Glycine max L.) cropping zones of China. The results generated from different assays showed that soybean bradyrhizobial isolates comprised four genomic groups. Group I was composed of strains mainly isolated from the North and Northeast plains of China. All four assays confirmed this group as phylogenetically divergent from all the reference strains. Strains of the group may represent a new species. Strains in Group II isolated from a variety of geographic regions were ascribed to B. liaoningense. Strains in Group III, mainly isolated from Central and East China, were closely related to the reference strains of B. japonicum. Strains in Group IV belonged to B. elkanii.