Skip to main content

Complete genome sequence of Arthrobacter sp. PAMC25564 and its comparative genome analysis for elucidating the role of CAZymes in cold adaptation

Abstract

Background

The Arthrobacter group is a known set of bacteria from cold regions, the species of which are highly likely to play diverse roles at low temperatures. However, their survival mechanisms in cold regions such as Antarctica are not yet fully understood. In this study, we compared the genomes of 16 strains within the Arthrobacter group, including strain PAMC25564, to identify genomic features that help it to survive in the cold environment.

Results

Using 16 S rRNA sequence analysis, we found and identified a species of Arthrobacter isolated from cryoconite. We designated it as strain PAMC25564 and elucidated its complete genome sequence. The genome of PAMC25564 is composed of a circular chromosome of 4,170,970 bp with a GC content of 66.74 % and is predicted to include 3,829 genes of which 3,613 are protein coding, 147 are pseudogenes, 15 are rRNA coding, and 51 are tRNA coding. In addition, we provide insight into the redundancy of the genes using comparative genomics and suggest that PAMC25564 has glycogen and trehalose metabolism pathways (biosynthesis and degradation) associated with carbohydrate active enzyme (CAZymes). We also explain how the PAMC26654 produces energy in an extreme environment, wherein it utilizes polysaccharide or carbohydrate degradation as a source of energy. The genetic pattern analysis of CAZymes in cold-adapted bacteria can help to determine how they adapt and survive in such environments.

Conclusions

We have characterized the complete Arthrobacter sp. PAMC25564 genome and used comparative analysis to provide insight into the redundancy of its CAZymes for potential cold adaptation. This provides a foundation to understanding how the Arthrobacter strain produces energy in an extreme environment, which is by way of CAZymes, consistent with reports on the use of these specialized enzymes in cold environments. Knowledge of glycogen metabolism and cold adaptation mechanisms in Arthrobacter species may promote in-depth research and subsequent application in low-temperature biotechnology.

Peer Review reports

Background

The Arthrobacter genus is a member of the family Micrococcaceae, which belongs to the phylum Actinobacteria [1, 2]. Arthrobacter species are often isolated from soil, where they contribute to biochemical cycles and decontamination [3]. Additionally, these species have been isolated worldwide from a variety of environments, including sediments [4], human clinical specimens [5], water [6], glacier cryoconite [7], sewage [8], and glacier ice [9]. Cold environments represent about 75 % of the earth, and their study provides information about new microorganisms and their evolution in cold environments [10]. Psychrophilic microorganisms have colonized all permanently cold environments, from the deep sea to mountains and polar regions [11]. Cold-adapted microorganisms utilize a wide range of metabolic strategies to grow in diverse environments. In general, the ability to adapt to low temperatures requires that microorganisms sense a decrease in temperature, which induces upregulation of cold-associated genes [12]. Arthrobacter is a gram-positive obligate aerobe that requires oxygen to grow in a variety of environments. Obligate aerobes grow through cellular respiration and use oxygen to metabolize substances like sugars, carbohydrates, or fat to obtain energy [13, 14]. However, there is a still lack of research on how obligate aerobes acquire adequate energy in cold environments.

Carbohydrate active enzymes (CAZymes) have functions associated with biosynthesis, binding, and catabolism of carbohydrates. This classification system is based on amino acid sequence similarity, protein folds, and enzymatic mechanism. Thus, one can understand overall enzyme function and carbohydrate metabolism through CAZymes [15]. These enzymes are classified based on their catalytic activity: glycoside hydrolase (GH), carbohydrate esterase (CE), polysaccharide lyase (PL), glycosyltransferase (GT), and auxiliary activity (AA). In addition, CAZymes may have non-catalytic subunits like a carbohydrate-binding module (CBM). CAZymes are already well known in biotechnology, and their industrial applications are of interest to many researchers because they produce precursors for bio-based products such as food, paper, textiles, animal feed, and various chemicals, including biofuels [16, 17].

Most bacteria can use glycogen as an energy storage compound, and the enzymes involved in its metabolism are well known. A recent study showed the physiological impact of glycogen metabolism on the survival of bacteria living in extreme environments [18]. Some microorganisms can adapt quickly to continuously changing environmental conditions by accumulating energy storage compounds to cope with transient starvation periods. These strategies use glycogen-like structures such as polysaccharides composed of α-D-glycosyl units connected by α-1,4 linkages and branched by α-1,6 glycosidic linkages. Such biopolymers differ in their chain length and branching occurrence. To be used as carbon and energy sources, their glucose units are released by specific enzymes [19].

Microorganisms have synergistic enzymes capable of decomposing plant cell walls to release glucose. This phenomenon can be used for energy supply to maintain microbial growth [20]. Starch is an excellent source of carbon and energy for microbes that produce proteins responsible for extracellular hydrolysis of starch, in-cell absorption of fructose, and further decomposition into glucose [21]. In addition, strains that metabolize glycogen show important physiological functions, including use of energy storage compounds for glycogen metabolism. These pathways act as carbon pools that regulate carbon fluxes [22], and partly, this ability is attributed to CAZymes. Using comparative genome analysis of bacteria isolated from cold environments and the genetic patterns of CAZymes within them, this study provides an understanding of how survival adaptation can be achieved in extremely low-temperature environments.

Results and discussion

Profile of the complete genome of Arthrobacter sp. PAMC25564

As shown in Table 1, the complete genome of Arthrobacter sp. PAMC25564 is composed of a circular chromosome of 4,170,970 bp with a 66.74 % GC content. 3,829 genes were predicted on the chromosome of which 3,613 protein-encoding genes were functionally assigned, whereas the remaining genes were predicted as hypothetical proteins. We annotated 147 pseudogenes, 15 rRNA genes, and 51 tRNA genes distributed through the genome. From the predicted genes, 3,449 (90.08 %) were classified into 20 functional Clusters of Orthologous Groups (COG) categories, whereas the remaining 380 (9.92 %) remained un-classified. The most numerous COG categories were S genes with unknown function (705 genes), transcription (category K, 298 genes), amino acid transport and metabolism (category E, 280 genes), carbohydrate transport and metabolism (category G, 276 genes), and energy production and conversion (category C, 259 genes) (Fig. 1). Many of these genes are related to amino acid transport, transcription, carbohydrate transport, and energy production/conversion, which suggests that this strain utilizes CAZymes for energy storage and carbohydrate metabolism. Most bacteria rely on cell respiration to catabolize carbohydrates to obtain the energy used during photosynthesis for converting carbon dioxide into carbohydrates. The energy is stored temporarily in the form of high-energy molecules such as ATP and used in several cell processes [23, 24]. However, Arthrobacter is already known as a genus of bacteria that is commonly found in cold environments. All species in this genus are gram-positive obligate aerobes and as such require oxygen to grow. These organisms use oxygen to metabolize substances like sugars, polysaccharides, or fats, to obtain energy as cellular respiration [14]. Therefore, we predicted that the PAMC25564 strain could also utilize carbohydrate degradation to obtain energy through these results.

Table 1 Genome features of Arthrobacter sp. PAMC25564
Fig. 1
figure1

 A. Circular map of Arthrobacter sp. PAMC25564 genome, B. COG functional categories for forward coding sequences. Metabolism: C, energy production and conversion; G, carbohydrate transport and metabolism; E, amino acid transport and metabolism; F, nucleotide transport and metabolism; H, coenzyme transport and metabolism; I, lipid transport and metabolism; P, inorganic ion transport and metabolism; and Q, secondary metabolites biosynthesis, transport, and catabolism. Cell processing and signaling: D, cell cycle control, cell division, and chromosome partitioning; V, defense mechanisms; T, signal transduction mechanisms; M, cell wall/membrane/envelope biogenesis; N, cell motility; Z, mobilome, prophages, and transposons; and O, posttranslational modification, protein turnover, and chaperones. Information storage and processing: J, translation, ribosomal structure, and biogenesis; A, RNA processing and modification; K, transcription; L, replication, recombination and repair; and B, chromatin structure and dynamics

16 S rRNA phylogenetic analysis and average nucleotide identity (ANI) values

The identification of A. sp. PAMC25564 was verified using 16 S rRNA sequence analysis (Fig. 2). This strain is phylogenetically placed among Arthrobacter and Pseudarthrobacter species. The results from phylogenetic analysis, Basic Local Alignment Search Tool, and EzBio Cloud revealed closely related strains such as P. sulfonivirans ALL (T) (99.09 %), P. siccitolerans 4J27 (T) (98.48 %), A. ginsengisoli DCY81 (T) (98.23 %), and P. phenanthrenivorans Sphe3 (T) (98.13 %). These results confirmed that isolate PAMC25564 belongs to the family Micrococcaceae, phylum Actinobacteria. Recently, several Arthrobacter species have been reclassified into new genera, based on 16 S rRNA sequence similarities and chemotaxonomic traits such as peptidoglycan types, quinone systems, and/or polar lipid profiles [25]. It has been proposed to reclassify five genera within the genus Arthrobacter: Paenarthrobacter gen. nov., Pseudarthrobacter gen. nov., Glutamicibacter gen. nov., Paeniglutamicibacter gen. nov., and Pseudoglutamicibacter gen. nov. Among them, the genus Arthrobacter would be reclassified into the Pseudarthrobacter group as: A. chlorophenolicus, A. defluvii, A. equi, A. niigatensis, A. oxydans, A. phenanthrenivorans, A. polychromogenes, A. scleromae, A. siccitolerans, and A. sulfonivorans [26].

Fig. 2
figure2

Phylogenetic tree of Arthrobacter sp. PAMC25564. The phylogenetic tree was generated using the Maximum Likelihood method and Tamura-Nei model in MEGA X, based on 16 S rRNA sequences. This study used the 16 S rRNA sequence of 40 strains. The tree shows the relationship between three Arthrobacter strains and six Pseudarthrobacter strains and their phylogenetic positions as compared with that of our strain. The capital T in bold indicates our strain

In general, bacterial comparative genome analysis uses the ANI methods. As shown in Fig. 3, each ANI value ranged from 70.67 to 98.46 %. So we see that comparative genome results are much lower than the common ANI values of 92–94. The ANI analysis shows the average nucleotide identity of all bacterial orthologous genes that are shared between any two genomes and offers a robust resolution between bacterial strains of the same or closely related species (i.e., species showing 80–100 % ANI) [27]. However, ANI values do not represent genome evolution, because orthologous genes can vary widely between the genomes being compared. Nevertheless, ANI closely reflects the traditional microbiological concept of DNA-DNA hybridization relatedness for defining species, so many researchers use this method, since it takes into account the fluid nature of the bacterial gene pool and hence implicitly considers shared functions [28]. The results mean the PAMC25564 strain could either belong to the species from which Arthrobacter diverged, or could be a Pseudarthrobacter closely related new species. However, this study classified the strain and allocated a species through 16 S rRNA sequencing and ANI. While our classification is not conclusive, the PAMC25564 strain will probably be reclassified into the genus Pseudarthrobacter in a follow-up study.

Fig. 3
figure3

Orthologous Average Nucleotide Identity of Arthrobacter sp. PAMC 25564 and 25 other genomes calculated using OrthoANI. ANI results are colored yellow to red according to their value mean from 70 to 100 %. They were calculated using an OrthoANI in OAT (Orthologous Average Nucleotide Identity Tool). Strains belonging to the same species are marked with strong color. Arthrobacter sp.: PAMC25564 (NZ_CP039290.1), 24S4-2 (NZ_CP040018.1), YN (NZ_CP022436.1), QXT-31 (NZ_CP019304.1), Rue61a (NC_018531.1/CP003203.1), FB24 (NC_008541.1/CP000454.1), PAMC25486 (NZ_CP007595.1), ZXY-2 (NZ_CP017421.1), U41 (NZ_CP015732.1), DCT-5 (NZ_CP029642.1), PGP41 (NZ_CP026514.1), ERGS1:01 (NZ_CP012479.1), YC-RL1 (NZ_CP013297.1), Hiyo4 (AP014718.1), KBS0702 (NZ_CP042172.1), UKPF54-2 (NZ_CP040174.1), MN05-02 (AP018697.1), Hiyo8 (AP014719.1), and ATCC21022 (NZ_CP014196.1); Arthrobacter crystallopoietes: DSM 20117 (NZ_CP018863.1); Arthrobacter alpinus: R3.8 (NZ_CP012677.1), ERGS4:06 (NZ_CP013200.1), and A3 (NZ_CP013745.1); Pseudarthrobacter phenanthrenivorans: Sphe3 (NC_015145.1/CP002379.1); Pseudarthrobacter chlorophenolicus: A6 (NC_011886.1/CP001341.1); and Pseudarthrobacter sulfonivorans: Ar51 (NZ_CP013747.1). The black circle indicates our strain

CAZyme-encoding genes in Arthrobacter sp. PAMC25564

Among the 3,613 identified protein-encoding genes in PAMC25564, 108 were significantly annotated and classified into CAZyme groups (GH, GT, CE, AA, CBM, and PL) using dbCAN2. The results provided an insight into the carbohydrate utilization mechanisms of PAMC25564. Signal peptides gene retention predicted that 11 genes contained in CAZyme of strain PAMC26554 through Signal P tool. We found that proteins were distributed as follows: 33 GHs, 45 GTs, 23 CEs, 5 AAs, and 2 CBMs. However, no protein was assigned to the PL group. The GH gene annotations revealed that the PAMC25564 genome contains genes involved in glycogen and trehalose metabolism pathways such as β-glucosidase (GH1), glycogen debranching proteins (CBM48 and GH13_11), (1→4)-α-D-glucan 1-α-D-glucosylmutase (GH13_26), α-glucosyltransferase (GH13), α-trehalose phosphorylase (GH65), and 4-α-glucanotransferase (GH77) (Table 2). Previous studies showed the complex interplay of glycogen metabolism in colony development of Streptomycetes (in Actinomycetes species was only reported), showing that spore germination is followed by an increase in glycogen metabolism [29]. The underlying genetic and physiological mechanisms of spore germination remain unknown, but some mechanisms associated with the accumulation of nutrients such as biomass and storage materials in the substrate mycelium during morphological phases of development have been reported [30]. However, not much research has been done yet about whether gram-positive obligate aerobes have glycogen metabolism mechanisms. Recently, Shigella sp. PAMC28760, a pathogen isolated from Antarctica, was also reported to be able to adapt and survive in cold environments through glycogen metabolism [31]. Also, Bacillus sp. TK-2 has been reported to possess cold evolution adaptability through CAZyme genes related to degradation of polysaccharides including cellulose and hemicellulose [32]. These complete genome analyses uncover genomic information and evolutionary insights regarding diverse strains and species from cold environments. However, characteristics of glycogen metabolism in prokaryotes remain less well-studied than those in eukaryotes, and the metabolism of microorganisms isolated from cold environments are not well understood [33]. This study predicts the role of CAZymes in cold adaptation, specifically as being those PAMC25564 genes involved in glycogen and trehalose metabolism.

Table 2 List of CAZyme GH enzymes from Arthrobacter sp. PAMC25564

Comparison of Arthrobacter sp. PAMC25564 genome characteristics with those from closely related species

We compared CAZyme genes from Arthrobacter species to speculate about their bacterial lifestyles and identified relevant CAZymes for potential applications in biotechnology. Considering the accessibility of available genome data, the complete genomes of 26 strains were chosen for the comparative analysis of CAZymes: 19 genomes of Arthrobacter spp., one genome of A. crystallopoietes, three genomes of A. alpinus, and three genomes of Pseudarthrobacter spp. (Table 3). Our results showed that the number of total CAZymes in each genome ranged from a minimum of 56 (A. sp. YC-RL1) to a maximum of 166 (P. chlorophenolicus A6) (Fig. 4). We predicted that common CAZyme genes such as CE14, CE9, GH23, GH65, GT2, GT20, GT28, GT39, GT4, and GT51 would appear in each of the 26 genomes. However, when we compared strains isolated from cold environments, we found CAZyme genes were more common than what is found in the 26 comparison genomes. They include CE1, CE4, CE9, CE10, CE14, AA3, AA7, CBM48, GH1, GH3, GH13, GH15, GH23, GH25, GH38, GH65, GH76, GT2, GT4, GT20, GT28, GT39, and GT51. In particular, CAZyme members GH13, GH65, GH77, GT5, and GT20 (glycogen and trehalose-related genes) are involved in energy storage. This study focuses on those genes related to adaptations in metabolism that allow the species to withstand cold environments. These genes are involved in glycogen degradation and trehalose pathways and were found in strains PAMC25564, 24S4-2, FB24, Hiyo8, KBS0702, MN05-02, PGP41, QXT-31, U41, UKPF54-2, A6, Ar51, and sphe3. These Arthrobacter species isolated from extreme environments have a family of CAZymes and the related genes for proteins with a strong ability to store and release energy and this permits them to survive in such cold areas. We found that strain PAMC25564 had the largest number of CAZyme genes related to glycogen metabolism and the trehalose pathway. In general, CAZymes are a large group of proteins that are mainly responsible for the degradation and biosynthesis/modification of polysaccharides but not all the members of this group are secreted proteins. This study confirms small differences in the gene pattern of CAZymes between species (Additional file 1: Figure S1).

Table 3 Genome information of 26 strains including Arthrobacter sp. PAMC25564
Fig. 4
figure4

Comparative CAZyme-encoding genes found in the genome of Arthrobacter species. On the left side is the CAZyme family and on the right side is the total CAZyme number: GT, glycosyltransferase; GH, glycoside hydrolase; CE, carbohydrate esterase; CBM, carbohydrate-binding modules; and AA, auxiliary activity. The CAZyme family is colored, as indicated below the figure. Arthrobacter sp.: PAMC25564, 24S4-2, YN, QXT-31, Rue61a, FB24, PAMC25486, ZXY-2, U41, DCT-5, PGP41, ERGS1:01, YC-RL1, Hiyo4, KBS0702, UKPF54-2, MN05-02, Hiyo8, and ATCC21022 Arthrobacter crystallopoietes: DSM 20117 Arthrobacter alpinus: R3.8, ERGS4:06, and A3 Pseudarthrobacter phenanthrenivorans: Sphe3 Pseudarthrobacter chlorophenolicus: A6 and Pseudarthrobacter sulfonivorans: Ar51

Bacterial glycogen metabolism in a cold environment

Glycogen is an energy source for plants, animals, and bacteria and is one of the most common carbohydrates. Glycogen consists of D-glucose residues joined by α (1→4) links; and it is a structural part of cellulose and dextran [43]. Glycogen is a polymer with approximately 95 % of α-1, 4 linkages, and 5 % of α-1, 6 branching linkages. In bacteria, glycogen metabolism includes five essential enzymes: ADP-glucose pyrophosphorylase (GlgC), glycogen synthase (GlgA), glycogen branching enzyme (GlgB), glycogen phosphorylase (GlgP), and glycogen debranching enzyme (GlgX) [44]. To adapt and survive in a cold environment, organisms need well-developed functional energy storage systems, one of which is glycogen synthesis. Bacteria have a passive energy saving strategy to adapt to cold environmental conditions such as nutrient deprivation, by using slow glycogen degradation. Glycogen is hypothesized to function as long durability energy reserves, which have been reported as a Durable Energy Storage Mechanism (DESM) to account for the long-term survival of some bacteria in cold environments [45]. Metabolism of maltodextrin has been linked with osmoregulation and sensitivity of bacterial endogenous induction to hyperosmolarity, which is related to glycogen metabolism. Glycogen-generated maltotetraose is dynamically metabolized by maltodextrin phosphorylase (MalP) and maltodextrin glucosidase (MalZ), while 4-α-glucanotransferase (MalQ) is responsible for maltose recycling to maltodextrins [46]. Maltotetraose is produced using GlgB, MalZ, MalQ, and glucokinase (Glk), which act on maltodextrin and glucose. On the other hand, glucose-1-phosphate can be formed by MalP for glycogen synthesis or glycolysis [47]. Thus, glycogen degradation can play an essential role in bacterial adaptation to the environment. Additionally, maltose may form capsular α-glucan, which plays a role in environmental adaptation through the (TreS)-Pep2-GlgE-GlgB pathway [48, 49]. Previous studies indicate that trehalose is involved in bacterial adaptation to temperature fluctuation, hyperosmolarity, and desiccation resistance. Recently, the accumulation of trehalose and glycogen under cold conditions in Propionibacterium freudenreichii has been reported [50, 51]. Therefore, the role of glycogen in bacterial energy metabolism is closely linked to several metabolic pathways associated with bacterial persistence under environmental stresses such as starvation, drying, temperature fluctuations, and hyperosmolarity. Maltodextrin and trehalose pathways are examples of the relationship between glycogen and other metabolic pathways, as shown in Fig. 5. However, further exploration is needed to elucidate the relationship of glycogen with other compounds, and the mechanisms involved in bacterial persistence strategies [46]. The comparative analysis of predicted pathways for glycogen metabolism in Arthrobacter isolates (Additional file 2: Table S1), showed that in PAMC25564 the trehalose biosynthesis follows three metabolic pathways (OtsAB, TreYZ, and TreS) as in Mycobacterium [52]. The trehalose biosynthesis pathway is well known in numerous bacteria, for example, as a defense strategy involving the accumulation of trehalose. Three metabolic pathways to regulate osmotic stress have been reported in Corynebacterium glutamicum [53]. These three metabolic pathways are used for producing trehalose in C. glutamicum, where the gene galU/otsAB allows the increase of trehalose levels up to six times [54, 55]. This pathway was found in A. sp. PAMC25564, and it was predicted that such an isolate could produce energy in cold environments.

Fig. 5
figure5

Predicted pathways for glycogen and trehalose metabolism in Arthrobacter sp. PAMC25564 as cold adaptation response. Kyoto Encyclopedia of Genes and Genomes KEGG-predicted enzyme pathways (yellow squares), and dbCAN2-predicted CAZyme family (blue squares)

Glycogen metabolism and the trehalose pathway in Arthrobacter species

We investigated the glycogen metabolic pathways in each Arthrobacter species strain (Fig. 6). To determine the three pathways of glycogen metabolism and trehalose biosynthesis in Arthrobacter species coming from diverse environments, the level of dissimilarity was analyzed based on the composition of GH, GT, and other major enzymes from the 16 genomes. The analysis showed that only QXT-31, U41, and PGP41 shared the same genes and pathway as our strain, while other strains have a slightly different pattern. This study assumed that the PAMC25564 strain uses different pathways to obtain energy or degrade polysaccharides. Based on the above-mentioned pathway-related genes, we confirmed that strains YN, Rue61, PAMC25486, ZXY-2, ERGS1:01, YC-RL1, ATCC21022, R.3.8, and A3 lack the malQ gene. These strains of Arthrobacter species showed a low number of genes for the three main pathways of trehalose biosynthesis and these are responsible for maltose recycling to maltodextrins. Therefore, the energy supply may be compromised in such isolates. Although most strains showed galU/otsAB genes, strains Hyo8 and ERGS1:01 lack the otsB gene (Fig. 6; Additional file 2: Table S1). This result suggests that these strains would produce a significantly lower amount of trehalose than the isolates having the otsB gene. Additionally, we investigated the phosphotransferase system-related genes in strains R.3.8 and A3. These enzymes constitute another method used by bacteria for sugar uptake when the source of energy is phosphoenolpyruvate. As a result, these two strains probably produce polysaccharides by themselves or from an external source using phosphoenolpyruvate rather than consuming energy. Most of the compared strains were isolated in low-temperature (-18 to 15 ℃) and low-contamination environments. These species have been reported as highly anticipated strains due to their fast adaptability to the environment, and these results suggest this adaptability is related to glycogen metabolism, trehalose, and maltodextrin pathways, which may have an impact on industrial applications. Microorganisms with related genes can make trehalose production economical and are able to draw on their own energy. This result also predicted that trehalose metabolism in microorganisms depends on the requirement of bacterial metabolism in the given environmental conditions, and this is one of the characteristics of bacteria that grow in extreme environments such as the cold.

Fig. 6
figure6

Comparative analysis of predicted glycogen and trehalose metabolic pathways in Arthrobacter species. Glycogen and trehalose metabolism-associated genes are colored, as indicated below the figure. The figure shows the differences in direction, presence, and location of genes among the strains. GH65: α,α-trehalose phosphorylase; GH13_30: α-1,4-glucan-maltose-1-phosphate maltosyltransferase; TreY/GH13_26: 1→4;-α-D-glucan 1-α-D-glucosylmutase; MalQ: 4-α-glucanotransferase; 2.4.1.21: glycogen synthase, and ADP-glucose transglucosylase; OtsA/ GT20: trehalose 6-phosphate synthase; 3.2.1.28: α-trehalase; GalU: UTP-glucose-1-phosphate uridylyltransferase; GlgX/CBM + GH13_11: glycogen debranching protein; TreZ/GH13_10: malto-oligosyltrehalose trehalohydrolase, α-glucosidase, and α-trehalase; GlgA/GT20: glycogen synthase, and ADP-glucose transglucosylase; GlgC: glucose-1-phosphate adenylyltransferase; OtsB/GT20: trehalose 6-phosphate phosphatase; 2.7.1.201: PTS system, and sugar-specific IIA component; GlgB/CBM + GH13_9: 1,4-α-glucan branching enzyme; TreS/GH13_16: maltose α-D-glucosyltransferase/α-amylase; GlgE/GH13_3: α-1,4-glucan-maltose-1-phosphate maltosyltransferase, PTS system, and sugar-specific IIA component; GlgP: glycogen phosphorylase; and 3.2.1.10: oligo-1,6-glucosidase. The black circle indicates our strain. / mean one gene deletion and // mean two more than genes deletion

Conclusions

In this study, we elucidated the complete genome sequence of Arthrobacter sp. PAMC25564 and conducted a comparative genome analysis with other species for studying CAZyme patterns. We isolated bacteria from cryoconite under laboratory conditions and confirmed that the isolate is an Arthrobacter species, based on the analysis of 16 S rRNA sequences. Although the isolation of this species from extreme or contaminated environments has been reported previously, there are no reports on the use of CAZymes in cold environments. We predicted that Arthrobacter sp. PAMC25564 could produce energy autonomously as fast as it can adapt to the environment. The PAMC25564 strain genome is 4.17 Mb in size with a GC content of 66.74 %. The analysis of its complete genome suggested that the isolate has glycogen, trehalose, and maltodextrin pathways associated to CAZyme genes. We confirmed that PAMC25564 has 108 active CAZyme genes from the following groups, 5 AA, 2 CBM, 23 CE, 33 GH, and 45 GT. In addition, a comparative genome analysis of Arthrobacter species revealed that they adapt quickly to the environment. In conclusion, we expect the genome sequence analysis to provide valuable information regarding novel functional enzymes, especially CAZymes, which are active at low temperatures and can be used for biotechnological applications and fundamental research purposes. This study provides a foundation to understand how the PAMC25564 strain produces energy in an extreme environment.

Methods

Isolation and genomic DNA extraction of Arthrobacter sp. PAMC25564

The strain Arthrobacter sp. PAMC25564 was isolated from the Cryoconite of Wurmkogel, Ötztaler Alps, Austria (47°04’ N, 12°41’ E, 2820 m to a maximum elevation) using 0.1 X R2A agar (MB cell Ltd., Seoul, Korea). The strain was isolated at an environmental temperature of 20 ℃. The bacterial sample for DNA analysis was isolated at 15 ℃ by using a pure R2A agar. DNA from A. sp. PAMC25564 was extracted using a QIAamp DNA Mini Kit (Qiagen Inc., Valencia, CA, USA). The quantity and purity of genomic DNA were determined using a spectrophotometer (Biochrome, Libra S35PC, UK). The extracted DNA was checked by agarose gel electrophoresis to evaluate its quality. DNA was stored at -20 ℃ until use.

Genome sequencing and assembly of the whole genome of Arthrobacter sp. PAMC25564

Genome sequencing was performed using PacBio sequel single-molecule real-time (SMRT) sequencing technology (Pacific Biosciences, Menlo Park, CA, USA). SMRTbell library inserts (20 kb) were sequenced using SMRT cells. Raw sequence data were generated from 7,698 reads and 4,170,970 bp that were assembled de novo by using the hierarchical genome-assembly process (HGAP v.4) protocol [56] and HGAP4 assembly using SMRT analysis software (ver. 2.3; Pacific Biosciences, https://github.com/PacificBiosciences/SMRT-Analysis). The complete genome sequence was deposited in the GenBank database under the GenBank accession number NZ_CP039290.1.

Genome annotation of Arthrobacter sp. PAMC25564

The genome of strain PAMC25564 was annotated using the rapid annotation subsystem technology (RAST) server [57]. The predicted gene sequences were translated and searched in the National Center for Biotechnology Information (NCBI) non-redundant database, the COG from the eggnog v.4.5.1 database [58], and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database by cutoff value 0.01 [59]. A circular map of the PAMC25564 genome was prepared using the CGView comparison tool [60]. CAZyme gene analyses were carried out by running dbCAN tool [61] scans using hidden Markov model (HMM) profile downloaded from dbCAN2 HMMdb (version 7.0). At the same time, we obtained results from Signal P (version 4.0) about the presence of CAZyme genes [62]. The e-value cutoff was 1e-15 and the coverage cutoff was > 0.35. In addition, we used DIAMOND [63] (e-value < 1e102) and Hotpep [64] (frequency > 2.6, hits > 6) to improve the prediction accuracy.

Phylogenetic analysis

Strain PAMC25564 was compared with other Arthrobacter species using 16 S rRNA phylogenetic analysis. Alignments were performed using Basic Local Alignment Search Tool from the NCBI database and analyzed using EzBio Cloud (www.ezbiocloud.com). 16 S rRNA sequences were aligned using MUSCLE [65, 66] and MEGA X [67] to reconstruct a neighbor-joining tree and maximum likelihood tree with 1,000 bootstrap replications.

Comparative genomics of Arthrobacter sp. PAMC25564

We used all complete genome sequences of Arthrobacter species available in GenBank (https://www.ncbi.nlm.nih.gov). Firstly, we determined the relationship of PAMC25564 with other strains from the same species using complete genome sequences and checked their similarity by comparing values of ANI, calculated using an OrthoANI in OAT (the Orthologous Average Nucleotide Identity Tool) [68]. The genome information of several Arthrobacter species is available in GenBank, and we compared the CAZymes from registered species referenced in CAZy (http://www.cazy.org). First, based on the 16 S rRNA sequence, strains with similarity were selected, and then 25 strains with a complete genome were chosen. All those sequences were downloaded from the database and all CAZymes were reannotated using the dbCAN2 server. In fact, some kinds of strains (PAMC 25,486, ERGS1:01, and ERGS4:06) reported isolated from cold environments such as the Spitsbergen of Arctic, glacial of Himalaya.

Availability of data and materials

The datasets analyzed during the current study are available in the NCBI repository, accession numbers: NZ_CP039290.1for Arthrobacter sp. strain PAMC25564, complete genome; NZ_CP040018.1 for Arthrobacter sp. strain 24S4-2, complete genome; NZ_CP022436.1 for Arthrobacter sp. strain YN, complete genome; NZ_CP019304.1 for Arthrobacter sp. strain QXT-31, complete genome; NZ_CP003203.1 for Arthrobacter sp. strain Rue61a; NZ_CP000454.1 for Arthrobacter sp. strain FB24, complete genome; NZ_CP007595.1 for Arthrobacter sp. strain PAMC25486, complete genome; NZ_CP017421.1 for Arthrobacter sp. strain ZXY-2, complete genome; NZ_CP015732.1 for Arthrobacter sp. strain U41, complete genome; NZ_CP029642.1 for Arthrobacter sp. strain DCT-5, complete genome; NZ_CP026514.1 for Arthrobacter sp. strain PGP41, complete genome; NZ_CP012479.1 for Arthrobacter sp. strain ERGS1:01, complete genome; NZ_CP013297.1 for Arthrobacter sp. strain YC-RL1, complete genome; AP014718.1 for Arthrobacter sp. strain Hiyo4, complete genome; NZ_CP042172.1 for Arthrobacter sp. strain KBS0702, complete genome; NZ_CP040174.1 for Arthrobacter sp. strain UKPF54-2, complete genome; AP018697.1 for Arthrobacter sp. strain MN05-02, complete genome; AP014719.1 for Arthrobacter sp. strain Hiyo8, complete genome; NZ_CP014196.1 for Arthrobacter sp. strain ATCC21022, complete genome. NZ_CP018863.1 for Arthrobacter crystallopoietes DSM 20117, complete genome; NZ_CP012677.1 for Arthrobacter alpinus strain R3.8, complete genome; NZ_CP013200.1 for Arthrobacter alpinus strain ERGS4:06, complete genome; NZ_CP013745.1 for Arthrobacter alpinus strain A3, complete genome; CP002379.1 for Pseudarthrobacter phenanthrenivorans strain Sphe3, complete genome; CP001341.1 for Pseudarthrobacter chlorophenolicus strain A6, complete genome; NZ_CP013747.1 for Pseudarthrobacter sulfonivorans strain Ar51, complete genome.

Abbreviations

GH:

Glycoside Hydrolase

CE:

Carbohydrate Esterase

PL:

Polysaccharide Lyase

GT:

Glycosyltransferase

AA:

Auxiliary Activity

CBM:

Carbohydrate-Binding Module

ANI:

Average Nucleotide Identity

NCBI:

National Center for Biotechnology Information

References

  1. 1.

    H.J. Conn, I. Dimmick, Soil bacteria similar in morphology to Mycobacterium and Corynebacterium, J. Bacteriol. 1947; 54: 291–303. https://doi.org/10.1128/JB.54.3.291-303.1947.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    C. Koch, P. Schumann, E. Stackebrandt, Reclassification of Micrococcus agilis Ali-Cohen 1889; to the genus Arthrobacter as Arthrobacter agilis comb. nov. and emendation of the genus Arthrobacter, Int. J. Syst. Bacteriol. 1995; 45: 837–839. https://doi.org/10.1099/00207713-45-4-837.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    M. Unell, K. Nordin, C. Jernberg, J. Stenstrom, J.K. Jansson, Degradation of mixtures of phenolic compounds by Arthrobacter chlorophenolicus A6, Biodegradation 2008; 19: 495–505. https://doi.org/10.1007/s10532-007-9154-2.

  4. 4.

    S.G. Dastager, Q. Liu, S.K. Tang, S. Krishnamurthi, J.C. Lee, W.J. Li, Arthrobacter enclensis sp. nov., isolated from sediment sample, Arch. Microbiol. 2014; 196: 775–782. https://doi.org/10.1007/s00203-014-1016-9.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Y. Huang, N. Zhao, L. He, L. Wang, Z. Liu, M. You, F. Guan, Arthrobacter scleromae sp. nov., isolated from human clinical specimens, J. Clin. Microbiol. 2005; 43: 1451–1455. https://doi.org/10.1128/JCM.43.3.1451-1455.2005.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    H.W. Chang, J.W. Bae, Y.D. Nam, H.Y. Kwon, J.R. Park, K.S. Shin, K.H. Kim, Z.X. Quan, S.K. Rhee, K.G. An, Y.H. Park, Arthrobacter subterraneus sp. nov., isolated from deep subsurface water of the South Coast of Korea, J. Microbiol. Biotechnol. 2007; 17: 1875–1879.

    CAS  PubMed  Google Scholar 

  7. 7.

    R. Margesin, P. Schumann, D.C. Zhang, M. Redzic, Y.G. Zhou, H.C. Liu, F. Schinner, Arthrobacter cryoconiti sp. nov., a psychrophilic bacterium isolated from alpine glacier cryoconite, Int. J. Syst. Evol. Microbiol. 2012; 62: 397–402. https://doi.org/10.1099/ijs.0.031138-0.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    K.K. Kim, K.C. Lee, H.M. Oh, M.J. Kim, M.K. Eom, J.S. Lee, Arthrobacter defluvii sp. nov., 4-chlorophenol-degrading bacteria isolated from sewage, Int. J. Syst. Evol. Microbiol. 2008; 58: 1916–1921. https://doi.org/10.1099/ijs.0.65550-0.

  9. 9.

    Q. Liu, Y.H. Xin, X.L. Chen, H.C. Liu, Y.G. Zhou, W.X. Chen, Arthrobacter ruber sp. nov., isolated from glacier ice, Int. J. Syst. Evol. Microbiol. 2018; 68: 1616–1621. https://doi.org/10.1099/ijsem.0.002719.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    R. Cavicchioli, Cold-adapted archaea, Nat. Rev. Microbiol. 2006; 4: 331–343.https://doi.org/10.1038/nrmicro1390.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    R.Y. Morita, Psychrophilic bacteria, Bacteriol. Rev. 1975; 39: 144–167.

    CAS  Article  Google Scholar 

  12. 12.

    S. D’Amico, T. Collins, J.C. Marx, G. Feller, C. Gerday, Psychrophilic microorganisms: challenges for life, EMBO Rep. 2006; 7: 385–389. https://doi.org/10.1038/sj.embor.7400662.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    L.M. Prescott, J.P. Harley, D.A. Klein, Microbiology (3rd ed.). Wm. C. Brown Publishers. 1996; pp. 130–131. ISBN 0-697-29390-4.

  14. 14.

    M. Berney, C. Greening, R. Conrad, W.R Jacobs, G.M Cook, An obligately aerobic soil bacterium activates fermentative hydrogen production to survive reductive stress during hypoxia, Proc. Natl. Acad. Sci. USA 2014; 111: 11479–11484. https://doi.org/10.1073/pnas.1407034111.

  15. 15.

    A. Levasseur, E. Drula, V. Lombard, P.M. Coutinho, B. Henrissat, Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes, Biotechnol. Biofuels 2013; 6: 41. https://doi.org/10.1186/1754-6834-6-41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    J. Rytioja, K. Hildén, J. Yuzon, A. Hatakka, R.P. de Vries, M.R. Mäkelä, Plant-polysaccharide-degrading enzymes from Basidiomycetes, Microbiol. Mol. Biol. Rev. 2014; 78: 614–649. https://doi.org/10.1128/MMBR.00035-14.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    A.K. Sista Kameshwar, W. Qin, Comparative study of genome-wide plant biomass-degrading CAZymes in white rot, brown rot and soft rot fungi, Mycology 2017; 9: 93–105. https://doi.org/10.1080/21501203.2017.1419296.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    A. Klotz, J. Georg, L. Bucinská, S. Watanabe, V. Reimann, W. Januszewski, R. Sobotka, D. Jendrossek, W.R. Hess, K. Forchhammer, Awakening of a dormant cyanobacterium from nitrogen chlorosis reveals a genetically determined program, Curr. Biol. 2016; 26: 2862–2872. https://doi.org/10.1016/j.cub.2016.08.054.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    A. Klotz, K. Forchhammer, Glycogen, a major player for bacterial survival and awakening from dormancy, Future Microbiol. 2017; 12: 101–104. https://doi.org/10.2217/fmb-2016-0218.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    W.R. de Souza, Microbial Degradation of Lignocellulosic Biomass, In: Chandel DA (ed), InTech. 2013; pp. 207–247.

  21. 21.

    V. Valk, W. Eeuwema, F.D. Sarian, R.M. van der Kaaij, L. Dijkhuizen, Degradation of granular starch by the bacterium Microbacterium aurum strain B8. A involves a modular α-amylase enzyme system with FNIII and CBM25 domains, Appl. Environ. Microbiol. 2015; 81: 6610–6620. https://doi.org/10.1128/AEM.01029-15.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Y.J. Goh, T.R. Klaenhammer, Insights into glycogen metabolism in Lactobacillus acidophilus: impact on carbohydrate metabolism, stress tolerance and gut retention, Microb. Cell Fact. 2014; 13: 94. https://doi.org/10.1186/s12934-014-0094-3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    D.L. Nelson, M.M. Cox, A.L. Lehninger, Lehninger principles of biochemistry 6th (ed), W. H. Freeman & Co., New York 2013.

  24. 24.

    L.M. Sanders, Carbohydrate: Digestion, absorption and metabolism. Encyclopedia of food and health, 2016; pp. 643–650.

  25. 25.

    H.J. Busse, P. Schumann, Reclassification of Arthrobacter enclensis as Pseudarthrobacter enclensis comb. nov., and emended descriptions of the genus Pseudarthrobacter, and the species Pseudarthrobacter phenanthrenivorans and Pseudarthrobacter scleromae, Int. J. Syst. Evol. Microbiol. 2019; 69: 3508–3511. https://doi.org/10.1099/ijsem.0.003652.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    H.J. Busse, Review of the taxonomy of the genus Arthrobacter, emendation of the genus Arthrobacter sensu lato, proposal to reclassify selected species of the genus Arthrobacter in the novel genera Glutamicibacter gen. nov., Paeniglutamicibacter gen. nov., Pseudoglutamicibacter gen. nov., Paenarthrobacter gen. nov. and Pseudarthrobacter gen. nov., and emended description of Arthrobacter roseus, Int. J. Syst. Evol. Microbiol. 2016; 66: 9–37. https://doi.org/10.1099/ijsem.0.000702.

  27. 27.

    J. Goris, K.T. Konstantinidis, J.A. Klappenbach, T. Coenye, P. Vandamme, J.M. Tiedje, DNA-DNA hybridization values and their relationship to whole-genome sequence similarities, Int. J. Syst. Evol. Microbiol. 2007; 57: 81–91. https://doi.org/10.1099/ijs.0.64483-0.

  28. 28.

    C. Jain, L.M. Rodriguez-R, A.M. Phillippy, K.T. Konstantinidis, S. Aluru, High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries, Nat. Commun. 2018; 9: 5114. https://doi.org/10.1038/s41467-018-07641-9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    M. Yeo, K. Chater, The interplay of glycogen metabolism and differentiation provides an insight into the developmental biology of Streptomyces coelicolor, Microbiology 2005; 151: 855–861. https://doi.org/10.1099/mic.0.27428-0.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    D. Schneider, C.J. Bruton, K.F. Chater, Duplicated gene clusters suggest an interplay of glycogen and trehalose metabolism during sequential stages of aerial mycelium development in Streptomyces coelicolor A32, Mol. Gen. Genet. 2000; 263: 543–553. https://doi.org/10.1007/s004380051200.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    S.R. Han, D.W. Kim, B. Kim, Y.M. Chi, S. Kang, H. Park, S.H. Jung, J.H. Lee, T.J. Oh, Complete genome sequencing of Shigella sp. PAMC28760: Identification of CAZyme genes and analysis of their potential role in glycogen metabolism for cold survival adaptation, Microb. Pathog. 2019; 137: 103759. https://doi.org/10.1016/j.micpath.2019.103759.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    L. Shen, X. Zang, K. Sun, H. Chen, X. Che, Y. Sun, G. Wang, S. Zhang, G. Chen, Complete genome sequencing of Bacillus sp. TK-2, analysis of its cold evolution adaptability, Sci. Rep. 2021; 11: 4836. https://doi.org/10.1038/s41598-021-84286-7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    D.H.D. Nguyen, S.H. Park, P.L. Tran, J.W. Kim, Q.T. Le, W. Boos, J.T. Park, Characterization of the transglycosylation reaction of 4-α-glucanotransferase MalQ; and its role in glycogen breakdown in Escherichia coli, J. Microbiol. Biotechnol. 2019; 29: 357–366. https://doi.org/10.4014/jmb.1811.11051.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    H. Niewerth, J. Schuldes, K. Parschat, P. Kiefer, J.A. Vorholt, R. Daniel, S. Fetzner, Complete genome sequence and metabolic potential of the quinaldine-degrading bacterium Arthrobacter sp. Rue61a, BMC Genomics 2012; 13: 534. https://doi.org/10.1186/1471-2164-13-534.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    C.H. Nakatsu, R. Barabote, S. Thompson, D. Bruce, C. Detter, T. Brettin, C. Han, F. Beasley, W. Chen, A. Konopka, G. Xie, Complete genome sequence of Arthrobacter sp. strain FB24, Stand. Genomic Sci. 2013; 9: 106–116. https://doi.org/10.4056/sigs.4438185.

    CAS  Article  Google Scholar 

  36. 36.

    R. Kumar, D. Singh, M.K. Swarnkar, A.K. Singh, S. Kumar, Complete genome sequence of Arthrobacter sp. ERGS1: 01, a putative novel bacterium with prospective cold active industrial enzymes, isolated from East Rathong glacier in India, J. Biotechnol. 2015; 214: 139–140. https://doi.org/10.1016/j.jbiotec.2015.09.025.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    L. Ren, Y. Jia, N. Ruth, B. Zhao, Y. Yan, Complete genome sequence of an aromatic compound degrader Arthrobacter sp. YC-RL1, J. Biotechnol. 2016; 219: 34–35. https://doi.org/10.1016/j.jbiotec.2015.12.008.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    S. Hiraoka, A. Machiyama, M. Ijichi, K. Inoue, K. Oshima, M. Hattori, S. Yoshizawa, K. Kogure, W. Iwasaki, Genomic and metagenomic analysis of microbes in a soil environment affected by the 2011 Great East Japan Earthquake tsunami, BMC Genomics 2016; 17: 53. https://doi.org/10.1186/s12864-016-2380-4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    D.A. Russell, G.F. Hatfull, Complete genome sequence of Arthrobacter sp. ATCC 21022, a host for bacteriophage discovery, Genome Announc. 2016; 4: e00168-16. https://doi.org/10.1128/genomeA.00168-16.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    R. Kumar, D. Singh, M.K. Swarnkar, A.K. Singh, S. Kumar, Complete genome sequence of Arthrobacter alpinus ERGS4: 06, a yellow pigmented bacterium tolerant to cold and radiations isolated from Sikkim Himalaya, J. Biotechnol. 2016; 220: 86–87. https://doi.org/10.1016/j.jbiotec.2016.01.016.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    A. Kallimanis, K.M. LaButti, A. Lapidus, A. Clum, A. Lykidis, K. Mavromatis, I. Pagani, K. Liolios, N. Ivanova, L. Goodwin, S. Pitluck, A. Chen, K. Palaniappan, V. Markowitz, J. Bristow, A.D. Velentzas, A. Perisynakis, C.C. Ouzounis, N.C. Kyrpides, A.I. Koukkou, C. Drainas, Complete genome sequence of Arthrobacter phenanthrenivorans type strain Sphe3, Stand. Genomic Sci. 2011; 4: 123–130. https://doi.org/10.4056/sigs.1393494.

    CAS  Article  Google Scholar 

  42. 42.

    J.A. Miranda-Ríos, J.A. Ramirez-Trujillo, B. Nova-Franco, L.F. Lozano-Aguirre Beltrán, G. Iturriaga, R. Suárez-Rodriguez, Draft genome sequence of Arthrobacter chlorophenolicus strain Mor30.16, isolated from the bean rhizosphere, Genome Announc. 2015; 3: e00360-15. https://doi.org/10.1128/genomeA.00360-15.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    S. Ball, C. Colleoni, U. Cenci, J.N. Raj, C. Tirtiaux, The evolution of glycogen and starch metabolism in eukaryotes gives molecular clues to understand the establishment of plastid endosymbiosis, J. Exp. Bot. 2011; 62: 1775–1801. https://doi.org/10.1093/jxb/erq411.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    L. Wang, M.J. Wise, Glycogen with short average chain length enhances bacterial durability, Naturwissenschaften 2011; 98: 719–729. https://doi.org/10.1007/s00114-011-0832-x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    L. Wang, Q. Liu, X. Tan, T. Yang, D. Tang, W. Wang, M.J. Wise, Bacterial glycogen as a durable energy reserve contributing to persistence: an updated bibliography and mathematical model, bioRxiv 2019. https://doi.org/10.1101/536110.

  46. 46.

    R. Dippel, W. Boos, The maltodextrin system of Escherichia coli: metabolism and transport, J. Bacteriol. 2005; 187: 8322–8331. https://doi.org/10.1128/JB.187.24.8322-8331.2005.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    K.H. Park, Roles of enzymes in glycogen metabolism and degradation in Escherichia coli, J. Appl. Glycosci. 2015; 62: 37–45. https://doi.org/10.5458/jag.jag.JAG-2015_005.

  48. 48.

    S. Bornemann, α-Glucan biosynthesis and the GlgE pathway in Mycobacterium tuberculosis, Biochem. Soc. Trans. 2016; 44: 68–73. https://doi.org/10.1042/BST20150181.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    J.T. Park, J.H. Shim, P.L. Tran, I.H. Hong, H.U. Yong, E.F. Oktavina, H.D. Nguyen, J.W. Kim, T.S. Lee, S.H. Park, W. Boos, K.H. Park, Role of maltose enzymes in glycogen synthesis by Escherichia coli, J. Bacteriol. 2011; 193: 2517–2526. https://doi.org/10.1128/JB.01238-10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    M. Dalmasso, J. Aubert, S. Even, H. Falentin, M.-B. Maillard, S. Parayre, V. Loux, J. Tanskanen, A. Thierry, Accumulation of intracellular glycogen and trehalose by Propionibacterium freudenreichii under conditions mimicking cheese ripening in the cold, Appl. Environ. Microbiol. 2012; 78: 6357–6364. https://doi.org/10.1128/AEM.00561-12.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    M. Reina-Bueno, M. Argandoña, J.J. Nieto, A. Hidalgo-Garcia, F. Iglesias-Guerra, M.J. Delgado, C. Vargas, Role of trehalose in heat and desiccation tolerance in the soil bacterium Rhizobium etli, BMC Microbiol. 2012; 12: 207. https://doi.org/10.1186/1471-2180-12-207.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    K.A.L. De Smet, A. Weston, I.N. Brown, D.B. Young, B.D. Robertson, Three pathways for trehalose biosynthesis in Mycobacteria, Microbiology 2000; 146: 199–208. https://doi.org/10.1099/00221287-146-1-199.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    A. Wolf, R. Kramer, S. Morbach, Three pathways for trehalose metabolism in Corynebacterium glutamicum ATCC13032 and their significance in response to osmotic stress, Mol. Microbiol. 2003; 49: 1119–1134. https://doi.org/10.1046/j.1365-2958.2003.03625.x.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    J. Carpinelli, R. Kramer, E. Agosin, Metabolic engineering of Corynebacterium glutamicum for trehalose overproduction: role of the TreYZ trehalose biosynthetic pathway, Appl. Environ. Microbiol. 2006; 72: 1949–1955. https://doi.org/10.1128/AEM.72.3.1949-1955.2006.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    L. Padilla, S. Morbach, R. Kramer, E. Agosin, Impact of heterologous expression of E. coli UDP-glucose pyrophosphorylase on trehalose and glycogen synthesis in Corynebacterium gluticum, Appl. Environ. Microbiol. 2004; 70: 3845–3854. https://doi.org/10.1128/AEM.70.7.3845-3854.2004.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    C.S. Chin, D.H. Alexander, P. Marks, A.A. Klammer, J. Drake, C. Heiner, A. Clum, A. Copeland, J. Huddleston, E.E. Eichler, S.W. Turner, J. Korlach, Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data, Nat. Methods 2013; 10: 563–569. https://doi.org/10.1038/nmeth.2474.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    R.K. Aziz, D. Bartels, A.A. Best, M. DeJongh, T. Disz, R.A. Edwards, K. Formsma, S. Gerdes, E.M. Glass, M. Kubal, F. Meyer, G.J. Olsen, R. Olson, A.L. Osterman, R.A. Overbeek, L.K. McNeil, D. Paarmann, T. Paczian, B. Parrello, G.D. Pusch, C. Reich, R. Stevens, O. Vassieva, V. Vonstein, A. Wilke, O. Zagnitko, The RAST Server: rapid annotations using subsystems technology, BMC Genomics 2008; 9: 75. https://doi.org/10.1186/1471-2164-9-75.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    J. Huerta-Cepas, D. Szklarczyk, K. Forslund, H. Cook, D. Heller, M.C. Walter, T. Rattei, D.R. Mende, S. Sunagawa, M. Kuhn, L.J. Jensen, C. von Mering, P. Bork, eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences, Nucleic Acids Res. 2016; 44: D286–D293. https://doi.org/10.1093/nar/gkv1248.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    M. Kanehisa, S. Goto, KEGG: Kyoto encyclopaedia of genes and genomes, Nucleic Acids Res. 2000; 28: 27–30. https://doi.org/10.1093/nar/28.1.27.

  60. 60.

    J.R. Grant, A.S. Arantes, P. Stothard, Comparing thousands of circular genomes using the CGView Comparison Tool, BMC Genomics 2012; 13: 202. https://doi.org/10.1186/1471-2164-13-202.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    H. Zhang, T. Yohe, L. Huang, S. Entwistle, P. Wu, Z. Yang, P.K. Busk, Y. Xu, Y. Yin, dbCAN2: a meta server for automated carbohydrate-active enzyme annotation, Nucleic Acids Res. 2018; 46: W95–W101. https://doi.org/10.1093/nar/gky418.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    T.N. Petersen, S. Brunak, G. von Heijne, H. Nielsen, SignalP 4.0: discriminating signal peptides from transmembrane regions, Nat. Methods 2011; 8: 785–786. https://doi.org/10.1038/nmeth.1701.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    B. Buchfink, C. Xie, D.H. Huson, Fast and sensitive protein alignment using DIAMOND, Nat. Methods 2015; 12: 59–60. https://doi.org/10.1038/nmeth.3176.

    CAS  Article  Google Scholar 

  64. 64.

    P.K. Busk, B. Pilgaard, M.J. Lezyk, A.S. Meyer, L. Lange, Homology to peptide pattern for annotation of carbohydrate-active enzymes and prediction of function, BMC Bioinformatics 2017; 18: 214. https://doi.org/10.1186/s12859-017-1625-9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    R.C. Edgar, MUSCLE: multiple sequence alignment with high accuracy and high throughput, Nucleic Acids Res. 2004; 32: 1792–1797. https://doi.org/10.1093/nar/gkh340.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    R.C. Edgar, MUSCLE: a multiple sequence alignment method with reduced time and space complexity, BMC Bioinformatics 2004; 5: 113. https://doi.org/10.1186/1471-2105-5-113.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    S. Kumar, G. Stecher, M. Li, C. Knyaz, K. Tamura, MEGA X: Molecular evolutionary genetics analysis across computing platforms, Mol. Biol. Evol. 2018; 35: 1547–1549. https://doi.org/10.1093/molbev/msy096.

  68. 68.

    I. Lee, Y.O. Kim, S.C. Park, J. Chun, OrthoANI: an improved algorithm and software for calculating average nucleotide identity, Int. J. Syst. Evol. Microbiol. 2016; 66: 1100–1103. https://doi.org/10.1099/ijsem.0.000760.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was a part of the project titled “Development of potential antibiotic compounds using polar organism resources (15250103, KOPRI Grant PM20030)”, funded by the Ministry of Oceans and Fisheries, Korea. In addition, this work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (NRF-2017R1A2B4012865).

Author information

Affiliations

Authors

Contributions

H. Park and T.-J. Oh designed and supervised the project. S.-R. Han and B. Kim performed the experiments; S.-R. Han, B. Kim, J.H. Jang, H. Park, and T.-J. Oh wrote the manuscript. All authors discussed the results, commented on the manuscript, and approved the manuscript.

Corresponding authors

Correspondence to Hyun Park or Tae-Jin Oh.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have declared that no conflict of interest exists.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Figure S1

. Comparative CAZyme-encoding genes found in the genome of Arthrobacter species. GT, glycosyl transferase; GH, glycoside hydrolase; CE, carbohydrate esterase; CBM, carbohydrate binding; and AA, auxiliary activity. CAZyme-encoding genes are colored, as indicated below the figure. Arthrobacter sp.: PAMC25564, 24S4-2, YN, QXT-31, Rue61a, FB24, PAMC25486, ZXY-2, U41, DCT-5, PGP41, ERGS1:01, YC-RL1, Hiyo4, KBS0702, UKPF54-2, MN05-02, Hiyo8, and ATCC21022; Arthrobacter crystallopoietes: DSM 20117; Arthrobacter alpinus: R3.8, ERGS4:06, and A3; Pseudarthrobacter phenanthrenivorans: Sphe3; Pseudarthrobacter chlorophenolicus: A6; and Pseudarthrobacter sulfonivorans: Ar51.

Additional file 2: Supplementary Table 1.

Comparative analysis of predicted pathways for glycogen and trehalose metabolism in Arthrobacter species. The symbol + indicates that the isolate produces the enzyme but symbol - indicates that the not produces the enzyme.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Han, SR., Kim, B., Jang, J.H. et al. Complete genome sequence of Arthrobacter sp. PAMC25564 and its comparative genome analysis for elucidating the role of CAZymes in cold adaptation. BMC Genomics 22, 403 (2021). https://doi.org/10.1186/s12864-021-07734-8

Download citation

Keywords

  • Arthrobacter species
  • CAZyme
  • Cold-adapted bacteria
  • Genetic patterns
  • Glycogen metabolism
  • Trehalose pathway