ORIGINAL RESEARCH article

Front. Microbiol., 09 June 2022

Sec. Extreme Microbiology

Volume 13 - 2022 | https://doi.org/10.3389/fmicb.2022.923038

Genomic Analysis Provides New Insights Into Biotechnological and Industrial Potential of Parageobacillus thermantarcticus M1

  • 1. Department of Bioengineering, Istanbul Medeniyet University, Istanbul, Turkey

  • 2. Institute of Biomolecular Chemistry (ICB), National Research Council, Naples, Italy

  • 3. Department of Bioengineering, Industrial Biotechnology and Systems Biology (IBSB), Marmara University, Istanbul, Turkey

Abstract

Parageobacillus thermantarcticus strain M1 is a Gram-positive, motile, facultative anaerobic, spore forming, and thermophilic bacterium, isolated from geothermal soil of the crater of Mount Melbourne (74°22′ S, 164°40′ E) during the Italian Antarctic Expedition occurred in Austral summer 1986–1987. Strain M1 demonstrated great biotechnological and industrial potential owing to its ability to produce exopolysaccharides (EPSs), ethanol and thermostable extracellular enzymes, such as an xylanase and a β-xylosidase, and intracellular ones, such as xylose/glucose isomerase and protease. Furthermore, recent studies revealed its high potential in green chemistry due to its use in residual biomass transformation/valorization and as an appropriate model for microbial astrobiology studies. In the present study, using a systems-based approach, genomic analysis of P. thermantarcticus M1 was carried out to enlighten its functional characteristics. The elucidation of whole-genome organization of this thermophilic cell factory increased our understanding of biological mechanisms and pathways, by providing valuable information on the essential genes related to the biosynthesis of nucleotide sugar precursors, monosaccharide unit assembly, as well as the production of EPSs and ethanol. In addition, gene prediction and genome annotation studies identified genes encoding xylanolytic enzymes that are required for the conversion of lignocellulosic materials to high-value added molecules. Our findings pointed out the significant potential of strain M1 in various biotechnological and industrial applications considering its capacity to produce EPSs, ethanol and thermostable enzymes via the utilization of lignocellulosic waste materials.

Introduction

Parageobacillus thermantarcticus strain M1 (DSM 9572T) is a thermophilic rod-shaped bacterium isolated from geothermal soil of the crater of Mount Melbourne (74°22′ S, 164°40′ E) during the Italian Antarctic Expedition in Austral summer of 1986–1987 ().

The taxonomical assignment of strain M1 underwent several changes over the years. Originally, at the time of the first taxonomic description, indicated it by the name of Bacillus thermoantarcticus, which was subsequently corrected to “thermantarcticus” in the Validation List N. 84 published in IJSEM (). Based on the 16S rRNA gene sequence analysis, strain M1 was placed as a peripheral member of the thermophilic Bacillus rRNA group 5. Afterward, illustrated the phenotypic and phylogenetic distance of strain M1 from other Bacillus species and transferred B. thermantarcticus to the genus Geobacillus as Geobacillus thermantarcticus comb. nov.

When the whole genome sequences of strain M1 and other Geobacillus species became available, analyzed the phylogenetic relatedness, applying the Average Amino acid Identity, the Average Nucleotide Identity and the digital DNA-DNA hybridization. Phylogenomic comparisons indicated that the genus Geobacillus clustered into two monophyletic clades that differed in terms of nucleotide base composition. Considering this premise, it was found that the Geobacillus species (comprising G. thermantarcticus), currently residing within clade II, were reconsidered as belonging to a new genus named Parageobacillus.

Several interesting enzymatic pathways were found in P. thermantarcticus, making this Antarctic species of great potential in biotechnological applications, ranging from bioprocessing tools to industrial biocatalysis. Thermostable enzymatic activities were studied, such as an intracellular xylose/glucose isomerase and a protease, as well as an extracellular xylanase (1,4-β-D-xylanxylanohydrolase; E.C. 3.2.1.8) and a β-xylosidase (1,4-β-D-xylanxylohydrolase; E.C. 3.2.1.37) (). Presence of xylanolytic enzymes was exploited to develop and optimize processes for the conversion of lignocellulosic biomass into high value added molecules like monosaccharides and valuable prebiotics (xylooligosaccharides) (; ). In particular, enzymatic conversions involved the hemicellulolytic fractions extracted by eco-sustainable processes from the stems and leaves of Cynara cardunculus and from rhizome biomass of Arundo donax, where only the seeds of the former and the stems of the latter are employed for bioethanol production when the enzymatic performances are compared with those described for the other Bacillus species, P. thermantarcticus enzymes stand out with remarkable optimum working conditions (80°C and pH 5.6 for xylanase and 70°C and pH 6.0 for β-xylosidase). Moreover, during the late exponential phase of growth, P. thermantarcticus was able to synthesize cytoplasmatic D-xylose (glucose) isomerase, which transforms both xylose and glucose into xylulose and fructose, respectively, and was remarkably stable at high temperatures (optimal temperature 90°C, pH 7.0) (). Different carbon sources were tested to improve the enzyme production and the best yield was obtained when xylose was employed as unique substrate in the growth media, both for xylose isomerase and glucose–isomerase activities.

reported the presence of an extracellular proteolytic enzyme in P. thermantarcticus and its maximum yield was detected at the end of the exponential growth phase. Skim milk in the growth medium was found to enhance the protease production up to 19-fold.

Besides the biotechnological potentials associated with these thermostable enzymes, successful studies have also demonstrated its employment as a suitable microbial model for Astrobiology, i.e., the multidisciplinary approach to the study of origin and evolution of life on Earth and in the Universe. Several experiments were carried out either in space or in laboratories, by simulating the extreme environmental parameters that are typical of space and exoplanets, such as irradiation with UV and γ-rays, desiccation, low temperature values, etc. (; , ; ). P. thermantarcticus viable cells and its spores were found to survive in a space simulated environment, which in turn lead to a new scenario called as the Panspermia theory that states that life on Earth could have originated from bacterial species from other places transported by solar pressure.

The Antarctic cell factory strain M1 was also studied for its ability to produce two exopolysaccharides (EPSs) using mannose as a sole carbon and energy source, with a yield of 400 mg/L, and spectroscopic analyses have revealed their chemical nature as xanthan and mannan-type polymers (). As with other EPS producing thermophiles, these biopolymers represent an ecological strategy to resist the extreme conditions: due to their protective nature, EPSs, surrounding the cells, represent a defense mechanism against extreme values of temperature and salinity; they are essential in retention of nutrients of water (; ). Without any doubt, EPSs are rapidly emerging as new and industrially important biomaterials, essentially thanks to their unique and complex chemical structures and properties, with commercial applications in various fields, from agriculture to medicine, microbial enhanced oil recovery (MEOR) and wastewater treatment, passing through the food and pharmaceutical fields, as well as packaging, textile and cosmetics use. Another important advantage is their high degree of their biocompatibility, biodegradability and both environmental and human compatibility, which make them viable alternatives to petroleum-based polymers (). Despite the excellent qualities of microbial EPSs, only a few have achieved commercial use due to their high production costs. The interrelations between metabolic pathways and EPS biosynthesis need to be elucidated in order to enhance the microbial productivity, and it is the mandatory step for the maximum exploitation and use of these biopolymers (). Therefore, a system-based approach becomes essential to pursue the objectives of industrial applications and for understanding the interplay between metabolism and bacterial EPS biosynthesis. Since microbial EPS biosynthesis is a result of a complex system of many metabolic pathways, the interplay of these networks needs to be enlightened to control, improve chemical characteristics and optimize production in order to achieve the formerly reported yields. In the view of this, the whole genome of M1 strain was sequenced1 by the U.S. Department of Energy Joint Genome Institute, via Illumina HiSeq 2500-1TB, and the NCBI Bioproject Accession number PRJNA323262 was assigned. For the data. The total number of sequenced bases was 3,448,881, of which 85.19% were coding bases and the DNA C + G% was found to be 43.63%. The total number of genes was 3,714, of which 96.85% were protein-coding and the remaining 3.15% were RNA genes. The number of protein-coding genes with function prediction was 2,783 (79.93%), while the number of those without function prediction were 814.

The genome sequencing and annotation can be a potent tool for understanding the genes utilized for biopolymer production. Moreover, elucidation of EPS biosynthetic pathways could lead to metabolic and genomic engineering strategies for the overproduction and tailor-made synthesis of biopolymers with desired activities and properties.

Hence in the present study, genome analysis of P. thermantarcticus M1 was performed with the aim of understanding the real potential of this organism. Elucidation of the whole genome organization of this thermophilic cell factory increased our understanding of biological mechanisms and pathways, by providing valuable information on the essential genes related to the biosynthesis of nucleotide sugar precursors, monosaccharide unit assembly, as well as the production of EPSs and ethanol. Presence of genes encoding xylanolytic enzymes that are required for the hydrolysis of hemicellulose indicated the potential use of this organism in the production of high-value added industrial products, such as ethanol, from lignocellulosic biomass. Our findings pointed out the significant potential of strain M1 in various biotechnological and industrial applications, considering its capacity to produce EPSs, ethanol and thermostable enzymes via the utilization of lignocellulosic waste materials.

Materials and Methods

Bacterial Strain and Culture Conditions

P. thermantarcticus (DSM 9572T) strain M1 (), a facultative anaerobic thermophilic Gram-positive bacterium, isolated from geothermal soil of Mount Melbourne in Antarctica, was cultured in a flask at 60°C for 12 h. Cells were grown in YN standard complex medium containing 6.0 g/L yeast extract and 3.0 g/L NaCl at pH 6.0 (; ). For the DNA analysis, the cells of strains M1 were collected at the exponential phase, via centrifugation at 5,000 × g for 15 min. Genomic DNA was isolated using DNAzol (Molecular Research Center, Inc., Cincinnati, OH, United States) according to the manufacturer’s instructions. The DNA was quantified spectrophotometrically using a Qubit™ fluorimeter (Thermo Fisher Scientific, MA, United States). The amount of 5 μg/μl of gDNA was considered for the construction of Illumina 300 bp library. The genome was sequenced at Department of Energy Joint Genome Institute under the project entitled: “The Genomic Encyclopedia of Bacteria and Archaea (GEBA).”

Whole Genome Sequence and Genome Assembly

The whole genome of P. thermantarcticus strain M1 was sequenced via Illumina HiSeq 2000 1TB technology. The library method was Illumina Regular Fragment, 300 bp, Plates with insert sizes of 238 ± 34 bp. The genome was assembled from a total of 6,360,534 sequence reads, with an average read length of 151 bp. The reads were assembled into 106 contigs with an average contig length of about 32,536 bp.

Genome Annotation

Genome annotation and gene prediction of P. thermantarcticus M1 were carried out using auto-annotation servers, such as Rapid Annotations using Subsystems Technology (RAST),2 and Pathosystems Resource Integration Center (PATRIC).3 The results obtained from the RAST Server was controlled by matching to that from PATRIC. The predictions of protein-encoding, rRNA and tRNA genes were made by RAST Server. The predictions of genes involved in important subsystems were manually verified via BLASTp4 against protein databases (i) the Universal Protein Resource (UniProt)5 and (ii) National Center for Biotechnology Information (NCBI).6 Subsystem annotation obtained from the RAST Server was used for function and category assignments to genes in the genome. Information on genes that encode enzymes was attained from Kyoto Encyclopedia of Genes and Genomes (KEGG)7 and ExPASy8 databases.

Phylogenomic Analysis

Phylogenomic comparison was performed between the whole genome sequence (WGS) of P. thermantarcticus M1 (GenBank accession number FOJS00000000.1), and WGS of the other species of the Parageobacillus genus [P. thermoglucosidasius DSM 2542 (GenBank accession number NZ_CP012712.1), P. toebii DSM 14590 (GenBank accession number RCWX00000000.1), P. caldoxylosilyticus CIC9 (GenBank accession number AMRO00000000.1), P. genome sp.1 NUB3621 2542 (GenBank accession number AOTZ00000000.1)]. NCBI database was used to download the WGS of the species and MicroScope Microbial Genome Annotation and Analysis Platform9 was used for the comparison of the genomes. MicroScope Microbial Genome Annotation and Analysis Platform assigned putative functions to the proteins encoded and these functions were used in order to compare the selected genomes. Comparison was performed with permissive parameters (50% amino-acid identity, 80% amino-acid alignment coverage).

Comparative Genomics

WGS of 18 strains belonging to 5 different Parageobacillus species [P. thermantarcticus M1 (GenBank accession number FOJS00000000.1), P. thermoglucosidasius ZCTH02-B4 (GenBank accession number LZRS00000000.1), P. thermoglucosidasius DSM 2542 (GenBank accession number NZ_CP012712.1), P. thermoglucosidasius C56-YS93 (GenBank accession numbers of chromosome CP002835.1 and plasmids CP002836.1/CP002837.1), P. thermoglucosidasius NCIMB 11955 (GenBank accession numbers of chromosome CP016622.1 and plasmids CP016623.1/CP016624.1), P. thermoglucosidasius TM242 (GenBank accession numbers of chromosome CP016916.1 and plasmids CP016917.1/CP016918.1), P. thermoglucosidasius TNO-09.020 (GenBank accession number CM001483.1), P. thermoglucosidasius NBRC 107763 (GenBank accession number BAWP00000000.1), P. thermoglucosidasius W-2 (GenBank accession number LXMA00000000.1), P. thermoglucosidasius GT23 (GenBank accession number LUCT00000000.1), P. toebii NBRC 107807 (GenBank accession number BDAQ00000000.1), P. toebii DSM 14590 (GenBank accession number RCWX00000000.1), P. toebii B4110 (GenBank accession number LQYW00000000.1), P. toebii PW12 (GenBank accession number QREZ00000000.1), P. caldoxylosilyticus NBRC 107762 (GenBank accession number BAWO00000000.1), P. caldoxylosilyticus CIC9 (GenBank accession number AMRO00000000.1), P. caldoxylosilyticus B4119 (GenBank accession number LQYS00000000.1), P. genomosp.1 NUB3621 2542 (GenBank accession number AOTZ00000000.1)] were obtained from the NCBI database and annotated automatically by the RAST Server and PATRIC. Encoding proteins that played roles in carbohydrate utilization, and industrially or biotechnologically important pathways, such as xylan hydrolysis, ethanol and, EPS biosynthesis, were found out for 4 genomes of different species and comparative analysis was performed according to the presence/absence of proteins.

Exopolymer Production and Recovery

The precultures of P. thermantarcticus strain M1 were incubated in YN standard complex medium for 8 h; then 100 ml were used to inoculate the culture. The cells were grown in Biostat C-plus Fermenter 10–3 (Sartorius Stedim; Melsungen, Germany) with a working volume of 4 L at 60°C with an agitation speed of 150 rpm. The fermentations were carried out under aerobic conditions with a constant air flow of 1.5 L min–1. pH was adjusted to 6.5 with 2M H2SO4 and 1M NaOH. Fermentation medium containing 0.6% (w/v) carbon source, 0.1% (w/v) yeast extract, and 0.3% (w/v) NaCl, was sterilized in situ at 120°C for 20 min, except that carbon source solution which was added to a sterile medium after filtration (0.22 μm size-pore). Cells were grown on various carbon sources, including fructose, galactose, maltose, xylose, mannose, trehalose, glycerol, lactose, glucose, sucrose. Growth was monitored by measuring the absorbance at 540 nm until the stationary phase was reached.

Samples taken from the cultures were centrifuged at 10,000 rpm for 30 min at 4°C. Equal volume of cold absolute ethanol was added to the polymeric substances in the extracellular fraction and kept at −20°C overnight. In order to precipitate exopolymer, the mixture was centrifuged for 30 min at 10,000 rpm at 4°C. Warm distilled water was used to dissolve the pellet and it was dialyzed against water in dialysis tubes (Spectra/Por MWCO; molecular weight cut-off 12–14 kDa) for 3 days, then lyophilized (Hetodrywinner) and weighed (). Phenol/sulfuric acid method () was used to assay the total carbohydrate content of exopolymers. All experiments were carried out in duplicate.

Ethanol Production and Quantification

The preculture of P. thermantarcticus strain M1 was incubated in anaerobic condition for determining the relation between microbial oxygen and the ethanol production capability. Therefore, after growth in strictly anaerobic medium, the broth was centrifuged at 10,000 rpm for 30 min at 4°C. The amount of ethanol produced by P. thermantarcticus strain M1 was determined from the cell-free supernatant using the Megazyme K-ETOH kit according to the manufacturer’s instructions.

Accession Numbers of Nucleotide Sequence

The genome project of P. thermoantarcticus M1 has been deposited at DDBJ/EMBL/GenBank under the accession FOJS00000000.1. Accession number of the first version was FOJS00000000.1. The gene sequence of 16S rRNA was available in DDBJ/EMBL/GenBank with the accession number KF192950.

Results

General Features of the Parageobacillus thermantarcticus M1 Genome

The gene prediction and genome annotation of P. thermantarcticus M1 revealed 3,957 coding sequences and 86 RNA genes (Table 1). A total of 2,580 protein coding genes were assigned with putative functions, whereas 1,377 hypothetical proteins couldn’t be matched to any known protein. Coding sequences assigned to any subsystem made up twenty-ninth percent of the total coding sequences (Supplementary Materials 1, 2). Furthermore, the gene annotation analysis resulted in the presence of the complete glycolysis, gluconeogenesis, Entner-Doudoroff, pentose phosphate pathways, TCA cycle, and genes encoding enzymes needed for ethanol fermentation. Moreover, pathways for the biosynthesis of folate, thiamine, and riboflavin were found to be present in P. thermantarcticus genome.

TABLE 1

GenomeP. thermantarcticus M1
DomainBacteria
TaxonomyBacteria; Terrabacteria group; Firmicutes; Bacilli; Bacillales; Bacillaceae; Parageobacillus; Parageobacillusthermantarcticus
Size3,448,881 bp
G + C content43.7%
Number of subsystems318
Number of coding sequences3,957
Proteins with functional assignments2,580
N50 value61,729
L50 value15
CRISPR repeats67
CRISPR spacer64
CRISPR array3
Number of RNAs
CheckM Completeness
CheckM Contamination
86
100
0.7

General features of P. thermantarcticus M1 draft genome.

In order to enlighten the functional characteristics of P. thermantarcticus M1, the association of genes with the general cluster of orthologous group (COG) functional categories was analyzed, and 65% of the total protein-coding genes were found to be associated with cellular processing and signaling, information storage and processing, and metabolism related functional categories (Table 2).

TABLE 2

COG codeNumber of genesPercentageDescription
Cellular process and signaling
D621.6004Cell cycle control, cell division, chromosome partitioning
M1233.1750Cell wall/membrane/envelope biogenesis
N792.0390Cell motility
O1173.0201Post-translational modification, protein turnover, and chaperones
T1543.9752Signal transduction mechanisms
U581.4972Intracellular trafficking, secretion, and vesicular transport
V481.2390Defense mechanisms
W20.0516Extracellular structures
Y00Nuclear structure
Z00Cytoskeleton
Information storage and processing
A00RNA processing and modification
B10.0258Chromatin structure and dynamics
J1674.3108Translation, ribosomal structure and biogenesis
K2315.9628Transcription
L2376.1177Replication, recombination and repair
Metabolism
C1794.6205Energy production and conversion
E2947.5891Amino acid transport and metabolism
F701.8069Nucleotide transport and metabolism
G2155.5498Carbohydrate transport and metabolism
H1353.4848Coenzyme transport and metabolism
I972.5039Lipid transport and metabolism
P1874.8271Inorganic ion transport and metabolism
Q571.4713Secondary metabolites biosynthesis, transport, and catabolism
Poorly characterized
R40010.3252General function prediction only
S2857.3567Function unknown
-1,16917.4501Not in COGs

Number of genes associated with the general cluster of orthologous group (COG) functional categories.

Phylogenomic Analysis

The relationships between P. thermantarcticus M1 and four other species of Parageobacillus genus were analyzed by using putative function assignments of encoded proteins, which is called a phylogenomics approach (Figure 1; ). The pan-genome describes the entire set of genes known to exist in all members of a given genus, consisting of a core genome (a gene set commonly shared between all individuals of the species) and the dispensable or accessory genome (a gene set shared within only one or some strains) (). Family/gene numbers responsible for pan-genome, core-genome and variable genome were 6,838/18,697, 1,851/11,030, and 4,987/7,667, respectively. P. thermantarcticus M1 has 900 unique proteins assigned with putative functions that form 24.207% of the total CDS. P. toebii DSM 14590 has lowest variable CDS (1216 CDS that form 35.913% of total CDS) and strain specific CDS (592 CDS that form 17.484% of total CDS). Highest pan CDS (3951 CDS) and variable CDS (1717 CDS) as well as strain specific CDS (928 CDS) were assigned to P. thermoglucosidasius DSM2542. Numbers of genes for five different species are presented in Table 3. P. thermantarcticus M1 shared the highest number of gene family with P. thermoglucosidasius DSM 2542 (2192 CDS) and the lowest number of gene family with P. toebii DSM 14590 (2045 CDS).

FIGURE 1

TABLE 3

OrganismPan CDSCore CDSVar CDSStrain specific CDSCore CDS (%)Var CDS (%)Strain specific CDS (%)
P. thermantarcticus M13,7182,1551,56390057.96142.03924.207
P. toebii DSM 145903,3862,1701,21659264.08735.91317.484
P. genomosp NUB36213,7312,2121,51967059.28740.71317.958
P. caldoxylosilyticus CIC93,9112,2591,65269757.76042.24017.822
P. thermoglucosidasius DSM25423,9512,2341,71792856.59343.45723.488

Pan/Core genome analysis of Parageobacillus genus.

Stress Tolerance/Response

Cells are generally exposed to many environmental changes during industrial production processes and therefore, they develop stress response mechanisms to sense and combat the deleterious effects of various stresses. A total of 34 genes encoding proteins associated with stress response were found in P. thermantarcticus genome.

Under stress conditions, the expressions of stress responsive genes are re-organized. The transcription factors that interact with RNA polymerase are activated to coordinate the gene expression (). In bacteria, sigma factors are required for transcription initiation. When a sigma factor associates with RNA polymerase, it directs promoter recognition. Since one subgroup of sigma factors comprises proteins that are activated in the presence of environmental stress (), the general stress regulons are controlled by the activity state of sigma factors. Therefore, the presence of a serine phosphatase RsbU (Pta.peg.737, Pta.peg.1300, Pta.peg.1512, Pta.peg.2721), which is a regulator of sigma subunit, in P. thermantarcticus genome was notable.

P. thermantarcticus showed two transporter genes encoding proteins that play roles in osmoregulation: glycerol uptake facilitator protein (Pta.peg.2109) and glycine betaine transporter OpuD (Pta.peg.3108) that is known from Bacillus subtilis. The OpuD system is a member of a small family of transport proteins involved in the accumulation of trimethylammonium compounds. It is a single component betaine glycine uptake system that actively participates in stress response of B. subtilis under high osmolarity environment (). The accumulation of glycine betaine at high concentrations was also reported in various bacterial species in response to osmotic stress (; ; ; ; ).

Oxidative stress response is one of the key responses found in P. thermantarcticus. Genes encoding proteins involved in the regulation of oxidative stress response, such as nitrite-sensitive transcriptional repressor NsrR (Pta.peg.2043), peroxide stress regulator PerR, FUR family (Pta.peg.107), ferric uptake regulation protein FUR (Pta.peg.1861), alkyl hydroperoxide reductase subunit C-like protein (Pta.peg.957) were identified in the genome. Moreover, genes encoding a cytosolic enzyme glutathione peroxidase GPX (EC 1.11.1.9, Pta.peg.3266, and Pta.peg.3324) that catalyzes the reduction of hydrogen peroxide to water and oxygen and the reduction of peroxide radicals to alcohols and oxygen were detected. The protection of glutathione against oxidative stress was reported for several bacteria (; ; ). In addition, three superoxide dismutase (SOD) proteins, including superoxide dismutase [Fe] (EC 1.15.1.1, Pta.peg.1893), superoxide dismutase [Cu-Zn] precursor (EC 1.15.1.1, Pta.peg.428), and superoxide dismutase [Mn] (EC 1.15.1.1, Pta.peg.3553) were present in the genome. SOD proteins protect the cells from reactive oxygen species by catalyzing the dismutation of superoxide radicals (O2) to either molecular oxygen (O2) or hydrogen peroxide (H2O2). Therefore, the presence of SOD proteins supported the existence of a possible defense mechanism against oxidative stress in P. thermantarcticus.

Other stress-related genes included a periplasmic stress related gene (intramembrane protease RasP/YluC, implicated in cell division based on FtsL cleavage, Pta.peg.3173), three bacterial hemoglobin genes (hemoglobin-like protein HbO, Pta.peg.2877, and di-guanylate cyclase/phosphodiesterase (GGDEF and EAL domains) with PAS/PAC sensors, Pta.peg.258 and Pta.peg.2815), two genes involved in carbon starvation (carbon starvation protein A, CstA, Pta.peg.2476, and carbon storage regulator, CsrA, Pta.peg.2671), and a hfl operon gene (RNA binding protein Hfq, Pta.peg.3104).

Finally, thermophilic microorganisms are able to not only tolerate high temperatures, but also exploit them to their advantage. Understanding high temperature tolerance of these microorganisms is essential for the analysis of thermophilic enzymes that can work at high temperatures. P. thermantarcticus M1 was found to have spermidine synthase (EC 2.5.1.16, Pta.peg.2044, and Pta.peg.2994) gene in its genome. The presence of the genes related to polyamine metabolism is common in the genome of thermophilic bacteria and considered to be associated with the thermophilicity of thermophiles (, ). Caseinolytic protease (Clp) family proteins play important role in heat stress response. These proteins are high molecular weight chaperones that are part of the Heat Shock Protein 100 (HSP100) family (). ClpB (Pta.peg.2905), ClpX (Pta.peg.3746), ClpC (Pta.peg.1716), ClpE (Pta.peg.1801) and ATP-dependent Clp protease proteolytic subunit (EC 3.4.21.92, Pta.peg.2605) were assigned in the genome of M1 strain. Moreover, the genome included genes encoding heat shock protein GrpE (Pta.peg.3602), chaperone proteins DnaK (Pta.peg.3601), DnaJ (Pta.peg.3600) and GroEL (Pta.peg.713), co-chaperone GroES (Pta.peg.712), and heat-inducible transcription repressor HrcA (Pta.peg.3603).

Transport

Transport is a vitally important system for all organisms on the Earth. P. thermantarcticus M1 showed 272 open reading frames (ORFs) (6.9% of total genome) related with various transport systems. These transport systems of the bacterium take role in the accumulation of necessary nutrients, extrusion of undesirable by-products and maintenance of cytoplasmic content of ions, salts, etc. A total of 113 ORFs were assigned to transport associated subsystems, whereas 159 ORFs were not categorized under a subsystem despite their roles in transport. The products of these ORFs consist of transporters for sugars and carbohydrates, ions (F,S2–, Mg2+, Zn2+, Mn2+, Fe2+/3+, Co2+, Ca2+, and Mo2+), amino acids (arginine, glutamine), and several further molecules (β-xyloside vitamin B12, spermidine/putrescine, glycerol-3-phosphate, polyamine, siderophore and subtilin) (Supplementary Materials 1, 2). Furthermore, the genome included genes responsible for quorum sensing, for instance those encoding ATP-binding cassette (ABC) transporters (Pta.peg. 211, 2887, 2888, 2889, 2890, 3462, 3463, 3464, 3465, 3466) and two-component sensor histidine kinase proteins (Pta.peg.516, 799, 850, 1450, 1634, 2229, 3454, 3504) for detection of the autoinducers.

Uptake and Utilization of Carbon Sources

Genes coding the ABC transporters and transporters for carbon sources (ribose, xylose, arabinose, galactose, glycerol, maltose/maltodextrin, and lactose), and utilization systems for sucrose, fructose, raffinose, xylose, glycerol, lactose, galactose, ribose, and chitin were identified in P. thermantarcticus genome. Moreover, the phosphoenolpyruvate-dependent sugar phosphotransferase system (sugar PTS) that is a major carbohydrate active transport system were detected for trehalose (Pta.peg.2053), maltose (Pta.peg.2362), glucose (Pta.peg.2366), and fructose (Pta.peg.3377). In addition, genes encoding enzymes associated with glucose, mannose and trehalose metabolisms were detected (Supplementary Materials 1, 2). In order to validate sugar uptake and utilization, P. thermantarcticus were cultured in the presence of 10 different sugars as carbon sources. Growth profiles of cultures indicated that glucose, fructose, sucrose, galactose, lactose, maltose, xylose, glycerol, mannose, and trehalose were found to be utilized for growth and EPS production (Figure 2).

FIGURE 2

Exopolysaccharide Biosynthesis

P. thermantarcticus was previously reported to synthesize two different EPSs: (i) EPS1 is a sulfate heteropolysaccharide composed of mannose and glucose (1.0:0.7 molar ratio), (ii) EPS2 is a sulfate homopolysaccharide composed of mannose as the major component (). In the aforementioned study, the EPS fraction was obtained with all sugars tested, while a higher yield was obtained with mannose as the source of carbon and energy. Additionally, in this study, the EPS yields achieved from cultures grown in 10 different carbohydrates were compared (Figure 2 and Table 4). Although undetectable levels of EPSs were obtained in the presence of maltose and glycerol, the highest yields were reached by cultures grown on galactose (18.24 mg/l) and trehalose (14.12 mg/l).

TABLE 4

C-sourceTransport systemAssociated gene(s)
FructosePTS systemPta.peg.3377
GalactoseABC transporterPta.peg.687
MaltoseABC transporter/PTS systemPta.peg.416, 608, 790, 791, 793, 890, 891, 892, 3427, 3429, 3441, 3451/Pta.peg.2362
XyloseABC transporterPta.peg.687, 3450, 3452, 1449
MannosePTS systemPta.peg.3394, 3395, 3396
TrehalosePTS systemPta.peg.2053
GlycerolGlycerol uptake facilitator proteinPta.peg.2109
LactoseABC transporterPta.peg.687
GlucosePTS systemPta.peg.2366
SucroseABC transporterPta.peg.1261, 2473, 3022

Carbohydrate uptake and utilization by P. thermantarcticus M1.

Ethanol Biosynthesis

P. thermantarcticus produced ethanol only in strictly anaerobic conditions, by using glucose as a carbon source, the maximum amount reached was 1.85 mg/l. The comparison of annotated P. thermantarcticus metabolic pathways to the glycolysis and oxidative decarboxylation of pyruvate deposited in KEGG resulted in the existence of the genes encoding enzymes required for ethanol synthesis from glucose. Firstly, glucose was converted to pyruvate by glycolysis pathway and then pyruvate is further used to form Acetyl-CoA through oxidative decarboxylation (Supplementary Material 2). Then, acetate was formed from Acetyl-CoA by acetyl-CoA synthetase (EC 6.2.1.1, Pta.peg.3869), and aldehyde dehydrogenase (EC 1.2.1.3, Pta.peg.671, Pta.peg.1191, Pta.peg.1463, Pta.peg.3297, Pta.peg.3404) catalyzed the acetaldehyde formation from acetate. Finally, the product, ethanol, was synthesized by alcohol dehydrogenase (EC 1.1.1.1, Pta.peg.820, Pta.peg.1274, Pta.peg.1800, Pta.peg.3497). The presence of genes encoding these enzymes showed the metabolic capacity of P. thermantarcticus to produce ethanol by utilizing glucose.

Comparative Parageobacillus Genome Analysis

At time of writing, Parageobacillus includes 7 species and to date, a total of 4 WGS for the strains of the thermophilic species P. caldoxylosilyticus (; ), P. thermoglucosidasius (; ; ; ; ; ), P. toebii (; ; ), as well as P. thermantarcticus M1, were reported. Comparative analysis was performed for these 4 genomes on the basis of carbohydrate uptake/utilization and EPS biosynthesis mechanisms (Supplementary Material 3). All subsystems and proteins that are responsible for sugar uptake/utilization or EPS biosynthesis, but couldn’t be assigned to any subsystem, were also evaluated. Published literature and genome information of these 7 species indicated biological, biotechnological and industrial significance of Parageobacillus genus, as well as similarities and differences within the genus itself.

Parageobacillus species were able to utilize D-xylose. The xylA gene encoding D-xylose isomerase (EC 5.3.1.5), and xylB gene encoding D-xylulose kinase (EC 2.7.1.17) have been found in the genome of P. caldoxylosilyticus, P. thermoglucosidasius, P. toebii and P. thermantarcticus. The xyl gene cluster also comprised genes encoding the D-xylose transport function. Coding genes for D-xylose transportation, such as ATP-binding protein (xylG), substrate-binding protein (xylF) and transport permease protein (xylH) were present in the genomes of all species. On the other hand, genes responsible for the utilization of L-arabinose were present only in the genome of P. caldoxylosilyticus and P. thermoglucosidasius. L-arabinose transporter genes araG (ATP-binding protein) araF (substrate-binding protein) araH (permease protein), as well as araA (arabinose isomerase) were present in these two genomes. However, Parageobacillus species were unable to ferment D-allose, and the utilization pathway of D-allose was completely absent in the genome of all species.

Chitin, a β-1,4-linked polymer of N-acetylglucosamine (GlcNAc), is the second most abundant polysaccharide on Earth, and mainly originates from the cuticle of arthropods and the cell wall of filamentous fungi (). The ability to metabolize chitin provides a competitive advantage to soil-dwelling microorganisms. Indeed, soils are generally carbon-rich and nitrogen-poor environments (), and GlcNAc-containing polymers are a rich source for both essential nutrients. Subsystem of “Chitin and N-acetylglucosamine utilization” is present in the genome of four species.

D-Galacturonate and D-Glucuronate utilization and dTDP-rhamnose synthesis pathways were present in the genomes of P. caldoxylosilyticus and P. thermatarcticus. A total of 7 genes assigned to subsystem of “D-Galacturonate and D-Glucuronate utilization.” The gene cluster required for the degradation of D-glucuronate, including uxaC encoding uronate isomerase (EC 5.3.1.12), uxuB encoding D-mannonate oxidoreductase (EC 1.1.1.57) and uxuA encoding mannonate dehydratase (EC 4.2.1.8), was present in both P. caldoxylosilyticus and P. thermatarcticus. However, uxaB encoding altronate oxidoreductase (EC 1.1.1.58) was absent in all genomes and uxaA encoding altronate dehydratase (EC 4.2.1.7) was present only in the genome of P. caldoxylosilyticus. Pathway of dTDP-L-rhamnose synthesis from D-glucose-1P was complete in both P. caldoxylosilyticus and P. thermatarcticus. RmlA encoding glucose-1-phosphate thymidylyltransferase (EC 2.7.7.24), rmlB encoding dTDP-glucose 4,6-dehydratase (EC 4.2.1.46), rmlC encoding dTDP-4-dehydrorhamnose 3,5-epimerase (EC 5.1.3.13), and rmlD encoding dTDP-4-dehydrorhamnose reductase (EC 1.1.1.133) gene cluster was annotated under “dTDP-rhamnose synthesis” subsystem.

There are three main routes for biosynthesis of microbial EPSs: (i) Wzx/Wzy-dependent pathway; (ii) ATP-binding cassette (ABC) transporter-dependent pathway, and (iii) synthase-dependent pathway (). Fructan or glucan based homopolysaccharides are synthesized by sucrase enzymes (). These glucan-sucrases and fructan-sucrases act on sucrose, hydrolyze glycosidic bonds, and then polymer chain is formed by transglycosylation reactions. The absence of genes encoding these sucrose enzymes in these four Parageobacillus species showed the lack of glucan or fructan biosynthesis, indicating that the EPSs are synthesized intracellularly like most of the EPSs from extremophilic bacteria (). The other two pathways for the synthesis of heteropolysaccharides and some homopolysaccharides are more complicated and requires the uptake of monosaccharides like glucose, xylose, fructose, mannose, or glycerol. Then, these monosaccharides are converted into nucleoside diphosphate sugars (NDP-sugars), and isoprenoid lipid–phosphate is formed in the cytoplasmic membrane assembly of NDP-sugars.

At this stage, repeating sugar monomer units are transferred sequentially from sugar nucleotides by glycosyltransferases, and then an acyl group is added for the modification, and lastly polymerization is completed. As a final point, the synthesized polysaccharide is secreted from the membrane of the cell into the extracellular environment (). Genome analysis showed the existence of a number of genes that might have roles in four different stages of EPS biosynthesis mechanism (Table 5). Comparative analysis of the genomes revealed that most of the genes responsible for biosynthesis of NDP-sugars existed in all Parageobacillus species. Genes encoding enzymes/proteins in the biosynthesis of isoprenoid lipid–phosphate and several glycosyltransferases were found in all compared genomes. All species were found to have deacetylase enzymes; however, chitooligosaccharide deacetylase (EC 3.5.1.-) was annotated for all species except P. thermantarcticus. On the other hand, polysaccharide pyruvyl transferase was annotated for P. thermantarcticus and P. caldoxylosilyticus (Supplementary Material 3). Two main mechanisms exist for the polysaccharide transport system. In the first scenario, a polysaccharide can be exported by ABC transporters across the cytoplasmic membrane. The existence of the related genes in the genome of P. thermantarcticus M1 supported the presence of this mechanism for the EPS export (Table 5). Second scenario is the polymerization by a wzx/wzy-dependent pathway after isoprenoid-linked intermediates are flipped through the cytoplasmic membrane. Comparative genome analysis indicated that Parageobacillus species used first mechanism for the transport of polysaccharides (Supplementary Material 3).

TABLE 5

Steps of EPS biosynthesis/MetabolitesEssential gene(s)Encoded protein(s)
Biosynthesis of NDP-sugars
UDP-glucosePta.peg.1692/
Pta.peg.2221
UTP-glucose-1-phosphate uridylytransferase (EC 2.7.7.9)/Galactose-1-phosphate uridylyltransferase (EC 2.7.7.10)
UDP-galactosePta.peg.392, 393, 1131, 2222UDP-glucose 4-epimerase (EC 5.1.3.2)
GDP-mannosePta.peg.3277Mannose-1-phosphate guanylyltransferase (EC 2.7.7.13)
UDP- N-acetylglucosaminePta.peg.1569N-acetylglucosamine-1-phosphate uridyltransferase (EC 2.7.7.23)/Glucosamine-1-phosphate N-acetyltransferase (EC 2.3.1.157)
UDP- N-acetylgalactosaminePta.peg.392, 393, 1131, 2222UDP-N-acetylglucosamine 4-epimerase (EC 5.1.3.7)/UDP-glucose 4-epimerase (EC 5.1.3.2)
UDP- N-acetylmannosaminePta.peg.876UDP-N-acetylglucosamine 2-epimerase (EC 5.1.3.14)
Assembly on lipid-acceptor
Di-trans,octa-cis-undecaprenyl diphosphatePta.peg.3176Undecaprenyl pyrophosphate synthetase (EC 2.5.1.31)
Di-trans,octa-cis-undecaprenyl phosphatePta.peg.408, 858Undecaprenyl-diphosphatase (EC 3.6.1.27)
Glycosyltransferases
Pta.peg.1414, 1614, 1708, 2637, 2692, 3017, 3276Glycosyltransferase (2.4.1.-)
Pta.peg.1707Glycosyltransferase family 1 (2.4.1.-)
Pta.peg.1694Glycosyltransferase family 2 (2.4.1.-)
Pta.peg.1360Glycosyltransferase family 8 (2.4.1.-)
Pta.peg.864Poly (glycerol-phosphate) α-glucosyltransferase (EC 2.4.1.52)
Pta.peg.3275Mannosyltransferase (2.4.1.-)
Pta.peg.1693Undecaprenyl-phosphate galactosephosphotransferase (EC 2.7.8.6)
Pta.peg.628UDP-N-acetylglucosamine–N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase (EC 2.4.1.227)
Pta.peg.2689Undecaprenyl-phosphate α-N-acetylglucosaminyl 1-phosphate transferase (EC 2.7.8.33)
Pta.peg.2256UDP-N-acetylglucosamine:L-malate glycosyltransferase
Pta.peg.2995, 3868Multimodular transpeptidase-transglycosylase (EC 2.4.1.129) (EC 3.4.-.-)
Pta.peg.865N-acetylmannosaminyltransferase (EC 2.4.1.187)
Modification
Pta.peg.3945Glycogen phosphorylase (EC 2.4.1.1)
Pta.peg.867Polysaccharide pyruvyl transferase
Pta.peg.743,1786Polysaccharide deacetylase
Transport of EPS
Pta.peg.2269, 2612, 3543, 3658, 3748TPR repeat protein
Pta.peg.3037, 3865GAF domain/HD domain protein/GGDEF domain
Pta.peg.78, 116, 401, 788, 891, 1167, 1387, 2833ABC transporter, ATP binding protein
Pta.peg.416, 687, 793, 890, 1257ABC transporter, sugar binding protein
Pta.peg.609, 791, 892, 1258, 1259, 3451ABC transporter, permease protein
Pta.peg.91, 92, 254, 784, 851, 986, 1022, 1337, 1521, 1632Lipoprotein

Predicted genes and proteins involved in EPS biosynthesis.

Essential genes associated with EPS biosynthesis were detected by genome annotation and together with experimental evidences, a hypothetical mechanism for EPS biosynthesis was generated (Figure 3). In this mechanism, the pathways for sugar uptake as well as the biosynthesis of NDP-sugars were shown. However, the annotated glycosyltransferases (Table 5) should be characterized more specifically such that their role in the whole process can be identified more precisely.

FIGURE 3

Features of Industrial and Biotechnological Interest

Thermophilic microorganisms are of the utmost importance considering the production of thermostable enzymes (). WGS of P. thermantarcticus M1 revealed the presence of many industrially important enzymes, including lipases, esterases, proteases, xylanases, β-xylosidases, and xylose isomerases (Figure 4 and Table 6).

FIGURE 4

TABLE 6

Enzyme nameEC numberAssociated gene(s)
LipaseEC 3.1.1.-Pta.peg.2190, 2553, 2708, 3304, 3412, 3912
EsteraseEC 3.1.1.-/ EC 3.1.4.-Pta.peg.258, 379, 606, 621, 916, 1406, 1612, 1849, 2472, 2520, 2594, 2596, 2792, 2815, 2847, 2848, 2937, 3058, 3188, 3350, 3389, 3757
β-XylanaseEC 3.2.1.8Pta.peg.3417, 3446
β-XylosidaseEC 3.2.1.37Pta.peg.3425, 3447
Xylose isomeraseEC 5.3.1.5Pta.peg.3410
PullulanaseEC 3.2.1.41Pta.peg.3892
NeopullulanaseEC 3.2.1.135Pta.peg.1260
Serine proteaseEC 3.4.21.-Pta.peg.161, 2372, 3049
MetalloproteaseEC 3.4.24.-Pta.peg.2997
NitrilaseEC 3.5.5.4Pta.peg.3281
AmidaseEC 3.5.1.-Pta.peg.219, 595, 684, 874, 1193, 1782, 2131, 3134, 3357, 3515, 3517, 3593
DipeptidaseEC 3.4.13.19Pta.peg.3127, 3841, 3896
Aldo-Keto reductaseEC 1.1.1.-Pta.peg.164, 1080
N-acetylglucosamine-6-phosphate deacetylaseEC 3.5.1.25Pta.peg.1904
Glucosamine-6-phosphate deaminaseEC 3.5.99.6Pta.peg.1905

Industrially important enzymes predicted by genome analysis of P. thermantarcticus M1.

β-xylanases and β-xylosidases are required for complete hydrolysis and assimilation of β-1,4-xylan, which is the most important hemicellulose component of lignocellulosic biomass. Xylanases, which have been mostly reported in bacteria and fungi, play a crucial role in pulp and paper industry by hydrolyzing xylan, and consequently releasing lignin from paper pulp, and reducing the use of chlorine as the bleaching agent (). In industrial processes, xylanases obtained from thermophilic organisms are preferred and Geobacillus genus come to prominence due its ability to utilize various carbon sources. reported the production of a highly thermostable xylanase from Geobacillus sp. strain DUSELR13, which had the potential in bioethanol production processes. Genome annotation analysis showed the presence of xylanolytic enzymes, β-xylanase and β-xylosidase, in P. thermantarcticus M1 genome. Moreover, the production of both enzymes by P. thermantarcticus M1 was previously reported and the characterization of these enzymes showed their potential use in various biotechnological applications (). In addition, the presence of the gene encoding alcohol dehydrogenase (EC 1.1.1.1, Pta.peg.820, Pta.peg.1274, Pta.peg.1800, Pta.peg.3497), that converts acetaldehyde to ethanol, endowed the strain with application in bioethanol production from lignocellulosic biomass.

Hemicelluloses can be used as a sole carbon source for growth of various microorganisms. Initially, hemicelluloses are degraded to D-xylose by exoenzymes and then D-xylose is transported into the cell. Xylose isomerase is required for the isomerization of D-xylose to D-xylulose. reported the production of a thermostable xylose isomerase by P. thermantarcticus. Genome annotation analysis showed the presence of the gene encoding xylulose kinase (EC 2.7.1.17, Pta.peg.3409) that catalyzes the phosphorylation reaction of D-xylulose to D-xylulose 5-phosphate, and the gene encoding ribulose-phosphate 3-epimerase (EC 5.1.3.1, Pta.peg.3249) that is required for the conversion of D-xylulose 5-phosphate to D-ribulose 5-phosphate, which enters the pentose phosphate pathway in P. thermantarcticus M1 genome.

P. thermantarcticus M1 genome carried a gene encoding pullulanase, a pullulan degrading enzyme. Pullulanase is known to play a crucial role in starch processing, leading to its widespread use in various applications in food, pharmaceutical, material, and bioenergy industries (). Moreover, the analysis indicated the presence of a gene encoding neopullulanase, a type I pullulan hydrolase (). These pullulan degrading enzymes can hydrolyze the glycosidic bonds of pullulan and starch, and forms reducing sugars. Therefore, the presence of pullulan degrading enzymes in P. thermantarcticus M1 might lead to its potential application in bioethanol production from starchy biomass.

P. thermantarcticus M1 genome also carried multiple genes encoding proteases, specifically serine and metalloproteases, that constitute one of the most important group of industrial enzymes. The production of protease enzyme from P. thermantarcticus M1 was previously reported and its characterization revealed its potential use in biotechnological applications (). However, future studies are required for experimental validation of the other enzymes presented in Table 6.

Discussion

The whole genome sequencing of P. thermantarcticus M1 was performed to elucidate the genetic and metabolic interrelations of this thermophilic bacteria and to enhance our understanding on its potential in industrial applications. P. thermantarcticus M1 was able to produce highly thermostable enzymes and degrade hemicellulose. In particular, it possessed the genes encoding both exo- and endo-xylanolytic enzymes, that are required for the decomposition of lignocellulosic biomass. Together with the previous experimental data, the whole genome analysis made this microorganism a promising cell factory for the development of bioprocesses that convert lignocellulosic waste into biofuels and value-added bioproducts according to biorefinery approach ().

P. thermantarcticus M1 genome annotation showed the presence of genes encoding protease and glucose (xylose) isomerase, underlying the great potential of this Antarctic bacteria in industrially important biotechnological applications. Furthermore, the identification of genes encoding many further valuable catalysts, such as lipase, esterase, pullulanase, aldo-keto reductase, and enzymes that plays roles in chitin utilization could contribute to the development of important industrial processes. In recent years, microbial esterases attracted much attention for their potential use in the degradation of plastics and microplastics due to their presence in all kind of environments (). Today, the existence of microplastics reach up to the food chain and affects human health (). Microplastic contamination is commonly caused by anthropogenic activities, such that plastics are used daily by humans and come from a wide range of sources from industry sectors to domestic activities. Increasing the yields of enzymatic degradation plays great importance in microplastics degradation. Moreover, the identification of novel enzymes through proteomic approach can be an outstanding technique to solve this environmental problem.

The aldo-ketoreductase superfamily comprises enzymes involved in detoxification, and act on substrates such as environmental pollutants, a large number of pharmaceuticals, drugs, and xenobiotics (). Therefore, the presence of the gene encoding an aldo-keto reductase is of fundamental importance in the discovery of new sources for biodegradation of pollutants and drugs from microbial origin. Moreover, the genome analysis indicated the presence of “Chitin and N-acetylglucosamine Utilization” subsystem. The ability to metabolize chitin, which is the second most abundant polysaccharide on Earth, provides a competitive advantage to soil-dwelling microorganisms ().

The identification of the gene encoding alcohol dehydrogenase that is required for ethanol production, could make this thermophile an ideal candidate for bioethanol production. P. thermantarcticus M1 genome was also found to encode for xylanolytic enzymes, and pullulanase involved in lignocellulose and starch degradation, respectively. These features unveiled the potential of this strain as a cell factory for bioethanol production from lignocellulosic biomass and/or starchy waste products.

The genes involved in quorum sensing mechanism, were also detected by genome annotation. At the time of writing, quorum sensing mechanism of P. thermantarcticus M1 has not yet studied and therefore, our findings open new insights for the exploitation of cell-cell communication, and communications between cells of different species, that is important for the coordination of community behavior and for the inhibition of competing species.

Apart from the microbial capability, many cellular mechanisms that allow P. thermantarcticus M1 to survive under stressful conditions, were elucidated by genome annotation. Cells encounter with stressful conditions during industrial production processes, and develop stress response mechanisms to combat the deleterious effects of various stresses. Genome annotation analysis revealed the genes specifically associated with osmotic and oxidative stresses P. thermantarcticus M1 genome.

Another interesting aspect addressed was the ability of P. thermantarcticus M1 to produce EPS from a variety of sugars. Our results indicated the developed versatility of the bacterium that enables it to utilize different types of carbohydrates as carbon and energy source for both growth and EPS production. Combined with its strong lignocellulolytic activity, the wide spectrum of carbon source utilization made P. thermantarcticus M1 to recruit the use of mixed sugars as substrate for EPS production through thermophilic bioprocessing.

Essential genes associated with EPS biosynthesis were detected by genome annotation, and together with the experimental evidences showed that the Antarctic microorganism transported the required EPS monomers across the cytoplasmic membrane by ABC transporters and produce EPS outside the cells. In microbial EPS production, the knowledge of mechanism is needed for process optimization, and improvement of the product quality. Since the genetic information of P. thermantarcticus, such as genome position, coding region, gene product function, and Enzyme Commission (EC) numbers, was identified by genome annotation, this study would play a significant role in the reconstruction of a metabolic model of this extremophilic microorganism. The elucidation of the pathways and molecular mechanisms leading to EPS production would help the development of metabolic engineering strategies for P. thermantarcticus in order to optimize medium, to enhance biopolymer production yields, to lower the production costs, and to improve the quality and efficiency of the process for industrial-scale applications.

This study clearly established the biotechnological and industrial potential of P. thermantarcticus M1. Although further experimental evidences are required to elucidate the underlying molecular mechanisms and to clarify the metabolism of this thermophilic bacterium, the genome annotation information would be crucial for further in silico and in vivo studies. Moreover, this understanding would help to design engineering strategies that lead to the optimized production processes resulted in improved properties of value-added industrial products, especially EPSs and biofuels.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, (FOJS00000000.1, NZ_CP012712.1, RCWX00000000.1, AMRO00000000.1, and AOTZ00000000.1).

Author contributions

IF, AP, ET, CK, and SY conceived the study. CK and SY carried out the genome analysis. IF and AP designed the experiments. IF, IR, and LLe carried out the experiments. LLa and BN analyzed the experimental data. CK, SY, ET, IF, and AP wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was partially supported by the grants from the Programma Nazionale di Ricerche in Antartide (PNRA; grant no. PNRA18_00232).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.923038/full#supplementary-material

References

Summary

Keywords

Parageobacillus thermantarcticus, thermophiles, genome, next-generation sequencing, exopolysaccharides

Citation

Yasar Yildiz S, Finore I, Leone L, Romano I, Lama L, Kasavi C, Nicolaus B, Toksoy Oner E and Poli A (2022) Genomic Analysis Provides New Insights Into Biotechnological and Industrial Potential of Parageobacillus thermantarcticus M1. Front. Microbiol. 13:923038. doi: 10.3389/fmicb.2022.923038

Received

18 April 2022

Accepted

18 May 2022

Published

09 June 2022

Volume

13 - 2022

Edited by

Rajesh K. Sani, South Dakota School of Mines and Technology, United States

Reviewed by

Jia Wang, Oak Ridge National Laboratory, United States; Tanvi Govil, South Dakota School of Mines and Technology, United States

Updates

Copyright

*Correspondence: Ilaria Finore, , orcid.org/0000-0003-2970-387XCeyda Kasavi, , orcid.org/0000-0003-3606-5140

†These authors have contributed equally to this work

This article was submitted to Extreme Microbiology, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics