Abstract
A novel Gram-stain-positive, aerobic bacterial strain, designated AK-R2A1-2 T, was isolated from the surface-sterilized needle leaves of an Abies koreana tree. Strain AK-R2A1-2 T had 97.3% and 96.7% 16S rRNA gene sequence similarities with Subtercola boreus K300T and Subtercola lobariae 9583bT, respectively, but formed a distinct phyletic lineage from these two strains. Growth of strain AK-R2A1-2 T was observed at 4–25 °C at pH 5.0–8.0. Strain AK-R2A1-2 T contained menaquinone 9 (MK-9) and menaquinone 10 (MK-10) as the predominant respiratory quinones. The major cellular fatty acids were anteiso-C15:0 and summed feature 8 (C18:1ω7c or/and C18:1ω6c), and the polar lipids included diphosphatidylglycerol (DPG) and three unknown aminolipids, AKL2, AKL3, and AKL4. The complete genome of strain AK-R2A1-2 T was sequenced to understand the genetic basis of its survival at low temperatures. Multiple copies of cold-associated genes involved in cold-active chaperon, stress response, and DNA repair supported survival of the strain at low temperatures. Strain AK-R2A1-2 T was also able to significantly improve rice seedling growth under low temperatures. Thus, this strain represents a novel species of the genus Subtercola, and the proposed name is Subtercola endophyticus sp. nov. The type strain is AK-R2A1-2 T (= KCTC 49721 T = GDMCC 1.2921 T).
Similar content being viewed by others
Introduction
Abies koreana is a species of native fir that grows in the high mountains (1000–1900 m) of South Korea1. The tree is tolerant to alkaline soils and prefers to grow in cold locations with heavy winter snowfall2. Abies koreana was classified as endangered by the International Union for Conservation of Nature (IUCN) due to many unverified reasons leading to a massive decline in numbers and sudden dried death of these trees. This tree has potential applications in traditional medicine for the treatment of colds, stomachache, indigestion, and rheumatic disease3. Other important secondary metabolites of Abies koreana, such as essential oils, have antioxidant, antimicrobial, and anti-inflammatory properties1,4.
The genus Subtercola, a member of the family Microbacteriaceae within the order Micrococcales, was first characterized by Männistö et al.5 with the type species Subtercola boreus DSM 13056 T. In January 2022, this genus contained four child taxa with published names (https://lpsn.dsmz.de/genus/subtercola). Species of the genus Subtercola have been isolated from extremely high-altitude cold environmental niches, including groundwater5, snowy mountains6, and lakes7. Members of the genus grow optimally at low temperatures5,6,7 and are aerobic, Gram stain-positive, and non-endospore-producing, with pale to yellow, circular, convex, and smooth colonies. The polar lipids of this genus are phosphatidylglycerol, diphosphatidylglycerol, one unknown phospholipid, and two glycolipids5, the major fatty acid is anteiso-C15:0, and MK-9 and MK-10 are the major respiratory quinones5,6,7. The genus Subtercola has a high G + C content that ranges from 64.4 to 69.0% (https://ncbi.nlm.nih.gov/genome/?term=subtercola).
Omics approaches have facilitated the investigation of microbial diversity and plant host-endophyte interactions in multiple ecological communities. Many endophytes are able to synthesize novel compounds to maintain a mutualistic relationship with the host. A recent investigation of the bacterial community of Abies koreana trees sought to identify plant–microbe interactions and define important secondary metabolites. During this analysis, 81 endophytic strains were isolated from surface-sterilized needle leaves and one of the novel strains, designated AK-R2A1-2 T, was selected for further exploration in the present study. The taxonomic status of strain AK-R2A1-2 T was evaluated using phenotypic, phylogenetic, genotypic, and chemotaxonomic data. Genome annotation indicated that this strain has multiple copies of cold-associated genes and the potential to produce secondary metabolites with antioxidant properties. The plant growth promotion effects of strain AK-R2A1-2 T were assessed in rice as a cold-sensitive crop. The results showed that this novel Subtercola species isolated from Abies koreana in Mount Halla has potential value for crop cultivation.
Results and discussion
Phylogenetic analysis
The assembled 16S rRNA gene amplicon of strain AK-R2A1-2 T had 1,457 base pairs (bp). Comparative sequence analysis showed that strain AK-R2A1-2 T belonged to the genus Subtercola and, shared the highest 16S rRNA gene sequence similarity with Subtercola boreus DSM 13056 T (97.3%), followed by 96.7% sequence similarity with Agreia bicolorata VKM Ac-1804 T, Agreia pratensis VKM Ac-2510 T, and Subtercola lobariae KCTC 33586 T (all 96.7%). Sequence similarities with the other two members of the genus Subtercola were 96.2% and 96.1%. Strain AK-R2A1-2 T was regarded as a novel species of the genus Subtercola based on its novel species recognition threshold value of 98.6%8. To identify the phylogenetic location of strain AK-R2A1-2 T, phylogenetic trees from 16S rRNA gene sequences were reconstructed using the MEGA program. The NJ phylogenetic tree showed that strain AK-R2A1-2 T formed a cluster with S. boreus DSM 13056 T, and all Subtercola members comprised a single group. The ML and ME phylogenetic trees supported this result (Fig. 1). Thus, strain AK-R2A1-2 T was confirmed as a member of the genus Subtercola rather than the genus Agreia. Based on the 16S rRNA gene sequence similarity and phylogenetic analysis, strains S. boreus DSM 13056 T, S. vilae DSM 105013 T, Subtercola frigoramans KCTC 49696 T, and Subtercola lobariae KCTC 33586 T were selected for further comparative study under the same conditions.
Neighbor-joining (NJ) phylogenetic tree based on 16S rRNA gene sequences showing the position of strain AK-R2A1-2 T. Bootstrap values (> 50%) were calculated using the NJ, maximum-likelihood (ML), and minimum evolution (ME) algorithms. Filled circles on the nodes indicate that the relationships were also identified using the ML and ME algorithms, whereas open circles indicate nodes were identified by either the ML or the ME algorithm. Scale bar: 0.0100 substitutions per nucleotide position.
Phenotypic characteristics
Phenotypic characteristics based on the optimal growth medium, temperature, pH, and salt tolerance were investigated. Strain AK-R2A1-2 T was shown to grow well on R2A, PDA, MEA, YEP, ISP2, and GYM medium, but did not on MA, LB, NA, and TSA. Cells were short and rod-shaped without flagella (0.2–0.3 µm in width and 0.3–1.6 µm in length, Fig. S1), and were Gram stain-positive, non-motile, catalase-positive, and oxidase-negative. Growth was observed in R2A medium at 4–25 °C (optimal 20 °C), pH 5–8 (optimal pH 5), and 0–1% NaCl (optimal 0%). Differential physiological characteristics that distinguish strain AK-R2A1-2 T from closely related strains are shown in Table 1.
Chemotaxonomic features
The cellular fatty acid profiles (> 1%) of strain AK-R2A1-2 T and closely related strains are shown in Table 2. The major fatty acids (> 10%) of strain AK-R2A1-2 T were anteiso-C15:0 (12.3%) and summed feature 8 (76.7%). This was most similar to S. lobariae KCTC 33586 T, which also contained anteiso-C15:0 (25.8%) and summed feature 8 (52.1%) as its major fatty acids. The major fatty acids of S. boreus DSM 13056 T were anteiso-C15:0 (46.0%), iso-C15:0 (20.7%), 2-OH C14:0 (10.4%), and summed feature 8 (11.8%), while those of S. vilae DSM 105013 T were anteiso-C15:0 (64.1%), iso-C16:0 (12.6%), and 2-OH C14:0 (15.6%), and those of S. frigoramans KCTC 49696 T were summed feature 8 (40.4%) and summed feature 3 (30.2%). Differences in major and minor fatty acids can differentiate strain AK-R2A1-2 T from other closely related members of the genus Subtercola (Table 2). The respiratory quinones of strain AK-R2A1-2 T were menaquinone 9 (MK-9) and menaquinone 10 (MK-10), which is similar to other closely related Subtercola members. The major polar lipids in strain AK-R2A1-2 T were diphosphatidylglycerol (DPG) and three unknown aminolipids, AKL2, AKL3, and AKL4. The polar lipid profile was more simple for strain AK-R2A1-2 T than for strains of other Subtercola species (Fig. S2).
Genome analyses
General and functional features of the AK-R2A1-2 T genome
After assembly using the Canu (version 1.7) de novo assembler, the whole genome of strain AK-R2A1-2 T was shown to comprise a single circular chromosome of 4,318,731 bp, and the N50 value was 12,176 bases, with coverage of 95x. The G + C content calculated based on the respective whole-genome sequence was 65.8%. A comparison between two copies of the 16S rRNA gene fragments with the whole-genome sequence showed that DNA sequence contamination did not occur during the genome assembly of strain AK-R2A1-2 T. The whole genome sequence of the novel strain was deposited into NCBI under accession number CP087997.1. The genome sequence was annotated using RAST9,10, and protein-coding sequences were determined using PGAP11. The NCBI PGAP annotation revealed that the AK-R2A1-2 T genome contained 3,874 protein-coding genes and 56 RNA genes, including two 5S rRNA, two 16S rRNA, two 23S rRNA, three non-coding RNA (ncRNA), and 47 tRNA genes (Table S1). A total of 3,674 protein-coding sequences (CDSs) were annotated using the cluster orthologous group (Fig. 2). The largest group of CDSs was classified as unknown (1,106 of the total CDSs, 30.1%), while the groups with the lowest CDSs numbers were classified as cytoskeleton and RNA processing and modification (1 of the total CDSs, 0.02%). Transcription (8.9%) and carbohydrate transport and metabolism (7.8%) gene had the next highest numbers of CDSs annotated from the whole genome (Fig. 2 and S3).
The cluster of orthologous groups (COG) classification of putative proteins in the whole genome sequence of strain AK-R2A1-2 T.
A whole genome-based phylogenetic tree was reconstructed using UBCG (version 3.0) and showed that strain AK-R2A1-2 T forms a subgroup with S. lobariae 9583bT, and that all members of the genus Subtercola comprise one group (Fig. 3). Concurrently, members of the genus Agreia form a single cluster, consistent with the phylogenetic tree created using 16S rRNA gene sequences. These findings suggested that strain AK-R2A1-2 T was a novel species of the genus Subtercola.
ML phylogenomic tree based on UBCGs (concatenated alignment of 92 core genes). Gene support index (GSI, left) and bootstrap values (right) are indicated at the nodes. Scale bar, 0.050 substitutions per position.
Genomic features associated with the cold-adaption of strain AK-R2A1-2 T
The RAST server, PGAP, and BlastKOALA pipeline were used to perform genome annotation and determine the metabolic pathways of strain AK-R2A1-2 T. The following information is predicted from the genome sequence analyses:
Stress response Members of the genus Subtercola grow well at low temperatures, and while strain AK-R2A1-2 T was isolated at 25 °C, it was shown to grow at 4 °C; the type strains of other Subtercola members can even grow at − 2°C5. The potential mechanism(s) by which bacteria invoke stress responses to adapt to environmental changes have not yet been fully elucidated. Bacterial adaptations to facilitate survival at low temperatures include cell membrane fluidity (increased unsaturated fatty acids) and protein-specific responses (refolding cold-damaged proteins)12. Cold-shock proteins (Csp) are produced by bacteria in response to a rapid decrease in temperature and act as RNA chaperones to help prevent mRNA misfolding13. The PGAP analysis in the current study showed that the AK-R2A1-2 T genome contained two genes encoding a cold-shock protein (UFS59601.1; UFS60376.1) and a cold shock domain-containing protein (UFS59609.1). When whole-genome mining of strain AK-R2A1-2 T was performed using RAST SEED and BlastKOALA, many well-studied, cold-inducible genes14,15,16,17,18,19,20 were annotated by the BlastKOALA pipeline, including aceE (encoding pyruvate dehydrogenase E1 component [EC:1.2.4.1]), aceF (pyruvate dehydrogenase E2 component, dihydrolipoamide acetyltransferase [EC:2.3.1.12]), cspA (cold shock protein), deaD (ATP-dependent RNA helicase DeaD [EC:3.6.4.13]), dnaA (chromosomal replication initiator protein), gyrA (DNA gyrase subunit A [EC:5.6.2.2]), hupB (DNA-binding protein HU-beta), infA (translation initiation factor IF-1), infB (translation initiation factor IF-2), infC (translation initiation factor IF-3), nusA (transcription termination/antitermination protein NusA), otsA (trehalose 6-phosphate synthase [EC:2.4.1.15 2.4.1.347]), otsB (trehalose 6-phosphate phosphatase [EC:3.1.3.12]), pnp (polyribonucleotide nucleotidyltransferase [EC:2.7.7.8]), rbfA (ribosome-binding factor A), recA (recombination protein RecA), and tig (trigger factor) (Table 3). Compatible solutes, a group of compounds that are important for osmotic stress and cold shock, can reduce membrane stress using various mechanisms. Biosynthetic genes of the recognized cold-protective solutes, glycine betaine (opuC, opuBD, opuA), choline (betaA), and proline (proA, proB, proC) were annotated21,22. In addition, genes responsible for the synthesis of the osmoprotectant proline from glutamate (proA, glutamate-5-semialdehyde dehydrogenase [EC:1.2.1.41]; proB, glutamate 5-kinase [EC:2.7.2.11]; proC (pyrroline-5-carboxylate reductase; [EC:1.5.1.2]), two choline dehydrogenases (betA, [EC:1.1.99.1]), and the glycine betaine transport system (opuC, opuBD, opuA) were identified in the AK-R2A1-2 T genome23,24,25,26,27 (Table 3). Using the RAST SEED annotation of strain AK-R2A1-2 T (Fig. S4), 22 genes were associated with the stress response. Of these, three were involved in osmotic stress (osmoregulation), 11 were related to oxidative stress [oxidative stress (6), glutathione: redox cycle (3), glutaredoxins (1), glutathionylspermidine and trypanothione (1)], and eight were related to detoxification. Some freezing-protective solutes, which can contribute to cold stress remission28,29, such as amino sugar and nucleotide sugars (26), amino acids (96), starch, and sucrose (27) were also detected in the genome information from strain AK-R2A1-2 T.
Motility Bacterial motility is critical to its establishment in and colonization of its host, and promotes plant growth30. Two genes related to motility, one for flagellar assembly (rpoD) and one for bacterial chemotaxis (rbsB), were annotated from the whole genome of strain AK-R2A1-2 T. The incomplete motility pathway of strain AK-R2A1-2 T explains why it was non-motile.
Fatty acid metabolism Membrane fatty acids are the major determinants for a sufficiently fluid membrane state required for low-temperature membrane function. The melting point of membrane fatty acids at cold temperatures is reduced by the synthesis of short-chain fatty acids, increases in the relative abundance of mono/polyunsaturated fatty acids, and increases in the anteiso/iso branched-chain fatty acids ratio31,32,33,34. The beta-oxidation-related degradation pathway, including 3-hydroxyacyl-CoA dehydrogenase [EC:1.1.1.35], acyl-CoA oxidase [EC:1.3.3.6], acyl-CoA dehydrogenase [EC:1.3.8.7], acetyl-CoA acyltransferase [EC:2.3.1.16], and enoyl-CoA hydratase [EC:4.2.1.17]35, was detected in the whole genome of strain AK-R2A1-2 T. In addition, fatty acid biosynthesis genes, fabD, fabF, fabG, fabH, fabL, OXSM, OAR1, and unsaturated fatty acid biosynthesis-related genes including tesB were detected in the whole genome.
Genome comparison and biosynthetic potential of the novel strain
The proposed threshold values of dDDH and ANI for delineating novel bacterial species are < 70% and < 95–96%, respectively36. In the current study, the ANI values were 76.8–80% between strain AK-R2A1-2 T and other members of the genus Subtercola, 73.6–74% between AK-R2A1-2 T and members of the genus Agreia, and < 74.7% between AK-R2A1-2 T and other genera in the family Microbacteriaceae. Consistently, strain AK-R2A1-2 T exhibited a dDDH value of 22.5–24.5% with members of the genus Subtercola, 20.1–20.5% with members of the genus Agreia, and < 20.7% with other genera in the family Microbacteriaceae (Table S2). Similar results were obtained for orthoANI values, which were < 80.3% with most closely related strains in the family Microbacteriaceae (Table S2).
Bacteria produce a large number of secondary metabolites that act as a reservoir of bioactive metabolites, with some exhibiting unique functions such as salt resistance37. In the genus Subtercola, there are limited studies on the production of specialized metabolites and the potential source of natural products remains untapped. In recent years, many computational methods have been developed to identify biosynthetic gene clusters (BGCs) in genomic data. AntiSMASH is a platform that is widely used for genome mining of secondary metabolites. The secondary metabolite clusters of the AK-R2A1-2 T genome were searched using antiSMASH (V. 6.0.1). Seven clusters were identified from the whole genome and found to contain non-alpha poly-amino acids like ε-polylysine, microansamycin (polyketide), alkylresorcinol (polyketide), kosinostatin (NPR + polyketide), beta-lactone, carotenoid (terpene), and lipopolysaccharide (saccharide: lipopolysaccharide) (Table S3 and Figs. S5–S11). Some of these compounds have vital antioxidant, antimicrobial, and anti-inflammatory properties. For example, NAPAA is a cationic peptide shown to prevent microbe proliferation and is approved as a food-grade cationic antimicrobial metabolite38, which has natural antioxidant and antimicrobial activities39. Carotenoid can be used to protect against oxidative stress in model systems40,41,42,43, and regulates membrane fluidity at low temperature using cell membrane adaptation44,45. Higher carotenoid levels may also increase resistance to cold shock46. Alkylresorcinol can protect epithelial cells against oxidative damage and has antioxidant47,48,49, antigenotoxic49 and potential antiglioma activity50. Kosinostatin has antitumor, antimicrobial, and antiproliferative functions51,52, while beta-lactone is a natural product with significant clinical applications provided by its antimicrobial, anticancer, and anti-obesity properties53,54. Microansamycin is a member of the macrolactams family, important bioactive compounds with anti-tuberculous and anti-tumor activity55. Similarly, its host, the Korean fir, produces rich polyphenol compounds like quercetin (flavonoid) and carotenoids to provide the strong antioxidant properties of Korean fir56. The essential oil of Korean fir includes the terpene group (monoterpenes, oxygenated monoterpenes, sesquiterpenes, diterpenes, and limonene)57,58 which also have strong antioxidant, antimicrobial, and anti-inflammatory activities1,4. The functional secondary metabolites of strain AK-R2A1-2 T were very similar to its host plant.
Effect of strain AK-R2A1-2 T on rice seedling growth under low temperature
Many endophytic bacteria have a broad host range. For example, Bacillus altitudinis, isolated from the wild plant species, Glyceria chinensis, can promote the growth of Arabidopsis thaliana, N. tabacum, corn, and soybean59. Bacillus megaterium RmBm31, an endophytic bacterium isolated from Retama monosperma, increases the biomass of A. thaliana60. Rice is a cold-sensitive crop, and its exposure to low-temperature stress (< 20 °C) during germination and early seedling growth, can have a negative impact on initial stand establishment61,62. To assess the impact of AK-R2A1-2 T on rice seedling growth under low temperatures, rice seeds were co-cultivated with strain AK-R2A1-2 T for 10 days at 20 °C, and total fresh weight, shoot weight, shoot length, root number, root length, and chlorophyll contents were examined in the rice seedling (Fig. 4). Shoot weight and shoot length were 1.67-fold and 1.36-fold higher, respectively, in the AK-R2A1-2 T-exposed seedlings than in controls (CK), while root length and root number were 1.42-fold and 1.40-fold higher, respectively, in the AK-R2A1-2 T-exposed seedings (Fig. 4). Strain AK-R2A1-2 T significantly improved rice seedling biomass and root morphological parameters under low temperatures. Whole genome analysis showed that these improvements were associated with phytohormone biosynthesis genes for indole-3-acetic acid (trpABCDES) and phosphate solubilization (pstABCS, phoHLU, and phnB). Indole-3-acetic acid can induce the growth of auxin-dependent lateral root formation, root hair development, and primary root growth, which promote plant growth. Meanwhile, phosphate-solubilizing bacteria enhance plant production by solubilizing insoluble phosphorus to make it available and increasing phosphorus nutrition. Secondary metabolites including terpenes were also shown to promote plant growth and some BGCs may produce new secondary metabolites because they display low similarity to known clusters. The function of these secondary metabolites on plant growth will require further study.
Effects of coating strain AK-R2A1-2 T on rice growth under low temperature. (A) Comparison of rice seedling growth in response to a control (CK, R2A broth) or AK-R2A1-2 T bacterial suspension (1 × 108 cfu/mL) on water agar plates at 20 °C for 10 days. (B) Quantification of total fresh weight, shoot weight, shoot length, root number, root length, and chlorophyll contents. Error bars indicate the standard deviation of the mean (n = 15). Asterisks indicate a statistically significant difference between the control and coating bacterial suspension (two-way ANOVA, **p < 0.01, and ****p < 0.0001). Experiments were repeated twice with similar results.
Conclusion
According to phylogenetic, phenotypic, and chemotaxonomic data, strain AK-R2A1-2 T was shown to belong to the genus Subtercola. The differentiating phenotypic properties and chemotaxonomic data presented in Table 1 preliminarily distinguish this strain as a novel species. Genome analyses of strain AK-R2A1-2 T using RAST SEED, BlastKOALA, and antiSMASH identified genes for a cold-active chaperone, osmotic, and oxidative stress responses, and DNA repair mechanisms, as well as the gene clusters responsible for host plant secondary metabolites. These data support the physiological properties of the psychrophilic genus Subtercola, demonstrate bacterial and plant viability at low temperatures, and define the secondary metabolites isolated from the host plant.
Description of Subtercola endophyticus sp. nov
Subtercola endophyticus (en.do.phy’ti.cus. Gr. pref. endo-, within; Gr. neut. n. phyton, plant; L. masc. adj. suff. -icus, adjectival suffix used with the sense of belonging to; N.L. masc. adj. endophyticus, endophytic, within plants, pertaining to the original isolation from plant tissues).
Colonies grown on R2A medium are white, circular (1–4 mm in diameter), smooth, and convex after a 4 day culture at 25 °C. Cells (0.2–0.3 µm in width and 0.3–1.6 µm in length) are short rod-shaped, strictly aerobic, Gram stain-positive, non-motile, lacking flagella, catalase-positive, and oxidase-negative. Growth occurs at 4–25 °C (optimum, 20 °C), at pH 5–8 (optimal 5), and with 0–1% NaCl (optimal 0%). Strain AK-R2A1-2 T grows on R2A, PD, MEA, YEP, ISP2, and GYM (optimum, PD, and R2A), but not on MA, LB, NA, and TSA medium. Strain AK-R2A1-2 T was positive for the following enzyme activities in API ZYM strips: esterase (C4), esterase lipase (C8), lipase (C14), leucine arylamidase, valine arylamidase, cysteine arylamidase, trypsin, acid phosphate, naphthol-AS-BI-phosphohydrolase, β-galactosidase, α-glucosidase, β-glucosidase, and α-mannosidase. In the API 20NE strips, strain AK-R2A1-2 T was positive for aesculin hydrolysis and β-galactosidase, but negative for all other activities. In the API 50CH strips, all tests were negative except aesculin ferric citrate. The major fatty acids detected in strain AK-R2A1-2 T were summed featured 8 and anteiso-C15:0, the respiratory quinones were MK-9 and MK-10 and the major polar lipids were DPG, AKL2, AKL3, and AKL4.
The GenBank accession numbers of the 16S rRNA gene and whole genome sequences of strain AK-R2A1-2 T are OL314543 and CP087997.1, respectively. The strain is available from the Korea Collection for Type Culture (KCTC 49721 T) and the Guangdong Microbial Culture Collection Center (GDMCC 1.2921 T).
Materials and methods
Sample isolation and culture conditions
Plant samples were collected from the top region of Mount Halla, Jeju Island, Republic of Korea (33°21′32.4″ N, 126°31′68″ E) in July 2020. Two grams of needle leaves from Abies koreana tree samples were surface-sterilized with 1% sodium hypochlorite, rinsed five times in sterile distilled water, and crushed with 20 mL of 1 × phosphate-buffered saline (PBS) in a grinder. Ground samples were serially diluted from 10−1 to 10−4 with 1 × PBS, and 100 µL aliquots were spread onto tenfold diluted Reasoner’s 2A agar (R2A, Difco), and incubated at 25 °C for 7 days. Morphologically distinct colonies were selected and subsequently streaked on fresh R2A medium. All strains were preserved in sterile skimmed milk (10%, w/v) at − 80 °C. Among all the isolates, one white, circular, convex, and smooth colony, designated AK-R2A1-2 T, was selected for further study. Subtercola boreus DSM 13056 T(=K300T), Subtercola vilae DSM 105013 T(=DB165T), Subtercola frigoramans KCTC 49696 T(=K265T), and Subtercola lobariae KCTC 33586 T (=9583bT) were obtained from the corresponding culture collections as closely related strains to analyze the taxonomic characteristics under comparable culture conditions. Unless otherwise stated, bacterial strains were grown on R2A or ISP2 for 5 days for the subsequent tests.
Strain identification and phylogenetic analysis
The 16S rRNA gene of strain AK-R2A1-2 T was amplified by PCR using universal primers 27F (5′-AGAGTTTGATCMTGGCTCAG-3′) and 1492R (5′-TACGGYTACCTTGTTACGACTT-3′) with bacterial DNA as the template. The purified PCR product was sequenced by BioFact (Daejeon, Republic of Korea) using primers 27F, 1492R, 518F (5′-CCAGCAGCCGCGGTAATACG-3′), and 800R (5′-TACCAGGGTATCTAATCC-3′). Sequencing reads were assembled using Vector NTI software (1.6.1), and an almost full-length 16S rRNA gene sequence was obtained and subjected to similarity-based BLAST analyses against the GenBank (https://blast.ncbi.nlm.nih.gov/Blast.cgi) and EzBioCloud (http://www.ezbiocloud.net)36 databases. Phylogenetic trees based on 16S rRNA gene sequences were reconstructed with the Molecular Evolutionary Genetics Analysis (MEGA V.7.0) program63, using the neighbor-joining (NJ), minimum evolution (ME), and maximum-likelihood (ML) algorithms with 1,000 bootstrap iterations. Litorihabitans aurantiacus HIsM16-52 T (GenBank accession number LC203065) was included as an outgroup.
Phenotypic characteristics
AK-R2A1-2 T cell morphology was observed by scanning electron microscopy at the Chuncheon Center, Korea Basic Science Institute (KBSI) after growing the strain on R2A agar plates for 4 days. Gram staining was conducted with a Gram staining kit (Difco) according to the manufacturer’s instructions. Cell motility was tested by growing the strain on a semi-solid R2A medium (0.4%). Catalase activity was determined by the production of bubbles after adding 3% (v/v) hydrogen peroxide solution to fresh cells64, while oxidase activity was measured using oxidase test strips (bioMérieux) that changed color to purple. Different media including R2A, potato dextrose agar (PDA, Difco), malt extract agar (MEA, Difco), yeast extract peptone agar (YEP, Difco), ISP medium No.2 (ISP2, Difco), and GYM medium (GYM: glucose 4 g/L, yeast extract 4.0 g/L, malt extract 10 g/L, CaCO3 2 g/L, agar 15 g/L), marine agar 2216 (MA, Difco), Luria–Bertani agar (LB, Difco), nutrient agar (NA, Difco), and trypticase soy agar (TSA, Difco) were screened for optimal growth of the tested strains. Various temperatures (4, 10, 15, 20, 25, 30, 37, 40, 45, and 60 °C) were assayed to identify the growth range of the novel strain. The pH range (3.0 to 12.0, adjusted with 1 N HCl and NaOH) for growth and NaCl tolerance (0 to 15%, w/v; 1% concentration increments) were measured in R2A liquid medium, monitoring the optical density at 600 nm (OD600) using a microplate spectrophotometer (Multiskan skyhigh, Thermo Fisher Scientific). Other biochemical features, such as the utilization of substrates, acid production from carbohydrates, and enzyme activities, were tested using API 20NE (bioMérieux) or API ZYM with NaCl 0.85% medium (bioMérieux) and API 50CH according to the manufacturer’s protocols. All closely related type strains were tested under the same conditions.
Chemotaxonomy features
Biomass for analysis of cellular fatty acids from strain AK-R2A1-2 T and closely related type strains was extracted from cells grown in R2A media. Using the Sherlock Microbial Identification System version 6.0 (MIDI) standard protocol, whole-cell fatty acid methyl esters were extracted and analyzed by gas chromatography (model 6890 N; Agilent) using the Microbial Identification software package65. Quinones were extracted from 100 mg of freeze-dried cells by shaking in chloroform: methanol (2:1, v/v) and purified by thin-layer chromatography as described by Collins et al.66. Analysis of the purified quinones was performed using reverse-phase high-performance liquid chromatography (HPLC) with ultraviolet (UV) absorbance detection at 270 nm. For extraction and analysis of polar lipids, 100 mg of freeze-dried cells were boiled in methanol, mixed and shaken with chloroform, evaporated, and identified by two-dimensional thin-layer chromatography (TLC) on Kieselgel 60 F254 plates (silica gel, 10 × 10 cm; Merck). Lipids spots were visualized by spraying with 0.2% ninhydrin (Sigma-Aldrich), α-naphthol, molybdenum blue (Sigma-Aldrich), 4% phosphomolybdic acid reagent, and Dragendorff’s solution to detect amino group-containing, sugar-containing, phosphorus-containing, total, and quaternary nitrogen-containing lipids, respectively.
Genome sequencing and assembly
Whole genome sequencing was performed by the Macrogen facility (Macrogen, Republic of Korea) using the PacBio Sequel System (Pacific Biosciences, Inc.) and the Illumina sequencing platform. Raw sequencing data were assembled using the Canu (version 1.7) de novo assembler67. Error corrections were performed by Pilon (version 1.21)68. Potential contamination of the assembled genome was checked using the ContEst16S algorithm69 to compare 16S rRNA within the whole genome. To verify that strain AK-R2A1-2 T belonged to the genus Subtercola, a whole genome-based phylogenetic tree was reconstructed using the up-to-date bacterial core gene set and pipeline (UBCG) as described by Na et al.70.
Genome annotation and comparative analysis
The National Center for Biotechnology Information Prokaryotic Genome Annotation Pipeline (PGAP) and Rapid Annotation of microbial genomes using Subsystem Technology (RAST) SEED9,10,11 were used for genome annotation. Metabolic pathways were reconstructed using BlastKOALA, which is based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database71. Genome mining for the presence of secondary metabolite gene clusters was performed using the antiSMASH program (https://antismash.secondarymetabolites.org/) with a “relaxed detection strictness” parameter and a “known cluster blast” feature including non-ribosomal peptide synthetases (NRPSs), type I and type II polyketide synthases (PKSs), lanthipeptides, lasso peptides, sactipeptides, and thiopeptides72. Comparative genome analyses including digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI) values, and orthoANI values between strain AK-R2A1-2 T and the most closely related members of the family Microbacteriaceae were performed using the Genome-to-Genome Distance Calculation (GGDC) webserver (http://ggdc.dsmz.de/)73, ANI calculator (https://www.ezbiocloud.net/), and standalone Orthologous Average Nucleotide Identity (OAT) software74, respectively. A graphical circular map of the AK-R2A1-2 T genome was constructed using CGView (https://cgview.ca/).
Rice growth and low-temperature conditions
Seeds of the Japonica rice cultivar Haepung were surface-sterilized with 1% (v/v) sodium hypochlorite for 10 min followed by repetitive washes with distilled water and soaked in AK-R2A1-2 T suspension (1 × 108 CFU/mL) or R2A broth as a control as described as Jiang et al.37. The control and AK-R2A1-2 T bacteria-coated seeds were incubated in the dark at 14 °C for 4 days in three replicated square plates (245 mm × 245 mm) containing 0.15% water agar. All the plates were incubated at 20 °C with a 16 h light/8 h dark photoperiod for 10 days. To estimate shoot weight, shoot length, root weight, root number, and root length, fifteen seedling plants from each experiment were considered. Chlorophyll contents were calculated as described by Fu et al.75. All experiments were performed in triplicate yielding similar results.
Data availability
The strain is available from the Korean Collection for Type Cultures (KCTC 49721T) and the Guangdong Microbial Culture Collection Center (GDMCC 1.2921T). The GenBank accession number of AK-R2A1-2T for the 16S rRNA gene is OL314543. The whole genome sequence of AK-R2A1-2T accession number is CP087997.1. The associated BioSample and BioProject accession numbers are SAMN23075623 and PRJNA678113, respectively. The taxonomy ID is 2,895,559.
References
Yoon, W.-J. et al. Abies koreana essential oil inhibits drug-resistant skin pathogen growth and LPS-induced inflammatory effects of murine macrophage. Lipids 44(5), 471–476 (2009).
Yumoto, I. et al. Bacillus asahii sp. nov., a novel bacterium isolated from soil with the ability to deodorize the bad smell generated from short-chain fatty acids. Int. J. Syst. Evol. Microbiol. 54(6), 1997–2001 (2004).
Yeşilada, E. et al. Traditional medicine in Turkey. V. folk medicine in the inner Taurus Mountains. J. Ethnopharmacol. 46(3), 133–152 (1995).
Lee, J.-H. & Hong, S.-K. Comparative analysis of chemical compositions and antimicrobial activities of essential oils from Abies holophylla and Abies koreana activities of essential oils from Abies holophylla and Abies koreana. J. Microbiol. Biotechnol. 19(4), 372–377 (2009).
Mannisto, M. K. et al. Subtercola boreus gen. nov., sp. nov. and Subtercola frigoramans sp. nov., two new psychrophilic actinobacteria isolated from boreal groundwater. Int. J. Syst. Evol Microbiol. 50, 1731–1739 (2000).
Si, H.-L. et al. Subtercola lobariae sp. nov., an actinobacterium of the family Microbacteriaceae isolated from the lichen Lobaria retigera. Int. J. Syst. Evol. Microbiol. 67(5), 1516–1521 (2017).
Villalobos, A. S. et al. Subtercola vilae sp. Nov., a novel actinobacterium from an extremely high-altitude cold volcano lake in Chile. Antonie Van Leeuwenhoek 111(6), 955–963 (2017).
Kim, M. et al. Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int. J. Syst. Evol. Microbiol. 64(2), 346–351 (2014).
Overbeek, R. et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42, 206–214 (2014).
Aziz, R. K. et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics 9, 75 (2008).
Tatusova, T. et al. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44(14), 6614–6624 (2016).
Zhang, Y. & Gross, C. A. Cold shock response in bacteria. Annu. Rev. Genet. 55, 377–400 (2021).
Snopkova, K. et al. Genome sequences of two Antarctic strains of Pseudomonas prosekii: insights into adaptation to extreme conditions. Arch. Microbiol. 202(3), 447–454 (2020).
Barria, C., Malecki, M. & Arraiano, C. M. Bacterial adaptation to cold. Microbiology (Reading). 159(Pt 12), 2437–2443 (2013).
Gualerzi, C. O., Giuliodori, A. M. & Pon, C. L. Transcriptional and post-transcriptional control of cold-shock genes. J. Mol. Biol. 331(3), 527–539 (2003).
Moll, I. et al. Effects of ribosomal proteins S1, S2 and the DeaD/CsdA DEAD-box helicase on translation of leaderless and canonical mRNAs in Escherichia coli. Mol. Microbiol. 44(5), 1387–1396 (2002).
Lelivelt, M. J. & Kawula, T. H. Hsc66, an Hsp70 homolog in Escherichia coli, is induced by cold shock but not by heat shock. J. Bacteriol. 177(17), 4900–4907 (1995).
Grainger, D. C. et al. Selective repression by Fis and H-NS at the Escherichia coli dps promoter. Mol Microbiol. 68(6), 1366–1377 (2008).
Giangrossi, M. et al. Selective expression of the beta-subunit of nucleoid-associated protein HU during cold shock in Escherichia coli. Mol Microbiol. 44(1), 205–216 (2002).
Di Pietro, F. et al. Role of the ribosome-associated protein PY in the cold-shock response of Escherichia coli. Microbiologyopen. 2(2), 293–307 (2013).
Hoffmann, T. & Bremer, E. Protection of Bacillus subtilis against cold stress via compatible-solute acquisition. J. Bacteriol. 193(7), 1552–1562 (2011).
Ko, R., Smith, L. T. & Smith, G. M. Glycine betaine confers enhanced osmotolerance and cryotolerance on Listeria monocytogenes. J. Bacteriol. 176(2), 426–431 (1994).
Bhargava, S. Genetically modified cyanobacterium Nostoc muscorum overproducing proline in response to salinity and osmotic stresses. J. Biosci. 31(2), 265–272 (2006).
Hilger, D. et al. Backbone structure of transmembrane domain IX of the Na+/proline transporter PutP of Escherichia coli. Biophys. J. 96(1), 217–225 (2009).
Kohler, C. et al. A comprehensive genomic, transcriptomic and proteomic analysis of a hyperosmotic stress sensitive alpha-proteobacterium. BMC Microbiol. 15, 71 (2015).
Scholz, A. et al. Osmotic stress response in Acinetobacter baylyi: identification of a glycine-betaine biosynthesis pathway and regulation of osmoadaptive choline uptake and glycine-betaine synthesis through a choline-responsive BetI repressor. Environ. Microbiol. Rep. 8(2), 316–322 (2016).
Zaprasis, A. et al. Uptake of amino acids and their metabolic conversion into the compatible solute proline confers osmoprotection to Bacillus subtilis. Appl. Environ. Microbiol. 81(1), 250–259 (2015).
Wargo, M. J. Homeostasis and catabolism of choline and glycine betaine: lessons from Pseudomonas aeruginosa. Appl. Environ. Microbiol. 79(7), 2112–2120 (2013).
Jiang, J. et al. Physiological and TMT-labeled proteomic analyses reveal important roles of sugar and secondary metabolism in Citrus junos under cold stress. J. Proteomics. 237, 104145 (2021).
Allard-Massicotte, R., et al., Bacillus subtilis early colonization of Arabidopsis thaliana roots involves multiple chemotaxis receptors. mBio. 7(6), 01664–16 (2016).
Kralova, S. Role of fatty acids in cold adaptation of Antarctic psychrophilic Flavobacterium spp. Syst. Appl. Microbiol. 40(6), 329–333 (2017).
Chattopadhyay, M. K. & Jagannadham, M. V. A branched chain fatty acid promotes cold adaptation in bacteria. J. Biosci. 28(4), 363–364 (2003).
Weber, M. H. et al. Role of the Bacillus subtilis fatty acid desaturase in membrane adaptation during cold shock. Mol. Microbiol. 39(5), 1321–1329 (2001).
Svensk, E. et al. PAQR-2 regulates fatty acid desaturation during cold adaptation in C elegans. PLoS Genet. 9(9), e1003801 (2013).
Hu, F. et al. Low-temperature effects on docosahexaenoic acid biosynthesis in Schizochytrium sp TIO01 and its proposed underlying mechanism. Biotechnol. Biofuels. 13, 172 (2020).
Yoon, S. H. et al. Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67(5), 1613–1617 (2017).
Jiang, L. M. et al. Genome insights into the novel species Jejubacter calystegiae, a plant growth-promoting bacterium in saline conditions. Diversity 13(1), 24 (2021).
Luz, C. et al. Antimicrobial packaging based on varepsilon-polylysine bioactive film for the control of mycotoxigenic fungi in vitro and in bread. J. Food Process. Preserv. 42(1), e13370 (2018).
Alirezalu, K. et al. Combined effect of natural antioxidants and antimicrobial compounds during refrigerated storage of nitrite-free frankfurter-type sausage. Food Res. Int. 120, 839–850 (2019).
Martorell, P. et al. Evaluation of carotenoids protection against oxidative stress in the animal model caenorhabditis elegans. Methods Mol. Biol. 2083, 387–401 (2020).
Mordi, R. C. et al. Free radical mediated oxidative degradation of carotenes and xanthophylls. Molecules 25(5), 1038 (2020).
Lucas, A., Morales, J. & Velando, A. Differential effects of specific carotenoids on oxidative damage and immune response of gull chicks. J. Exp. Biol. 217(Pt 8), 1253–1262 (2014).
El-Agamey, A. et al. Carotenoid radical chemistry and antioxidant/pro-oxidant properties. Arch. Biochem. Biophys. 430(1), 37–48 (2004).
Jagannadham, M. V. et al. Carotenoids of an Antarctic psychrotolerant bacterium, Sphingobacterium antarcticus, and a mesophilic bacterium, Sphingobacterium multivorum. Arch Microbiol. 173(5–6), 418–424 (2000).
Chattopadhyay, M. K. et al. Carotenoid pigments of an antarctic psychrotrophic bacterium Micrococcus roseus: temperature dependent biosynthesis, structure, and interaction with synthetic membranes. Biochem Biophys. Res. Commun. 239(1), 85–90 (1997).
Seel, W. et al. Carotenoids are used as regulators for membrane fluidity by Staphylococcus xylosus. Sci Rep. 10(1), 330 (2020).
Wang, Z. et al. Wheat alkylresorcinols protect human retinal pigment epithelial cells against H2O2-induced oxidative damage through Akt-dependent Nrf2/HO-1 signaling. Food Funct. 10(5), 2797–2804 (2019).
Elder, A. S., Coupland, J. N. & Elias, R. J. Effect of alkyl chain length on the antioxidant activity of alkylresorcinol homologues in bulk oils and oil-in-water emulsions. Food Chem. 346, 128885 (2021).
Parikka, K. et al. In vitro antioxidant activity and antigenotoxicity of 5-n-alkylresorcinols. J. Agric. Food. Chem. 54(5), 1646–1650 (2006).
Liu, J. et al. Effect of alkylresorcinols on the formation of Nepsilon-(carboxymethyl)lysine and sensory profile of wheat bread. Food Sci. Nutr. 9(1), 489–498 (2021).
Rambabu, V., Suba, S. & Vijayakumar, S. Antimicrobial and antiproliferative prospective of kosinostatin - a secondary metabolite isolated from Streptomyces sp. J. Pharm. Anal. 5(6), 378–382 (2015).
Back, C. R. et al. A new micromonospora strain with antibiotic activity isolated from the microbiome of a mid-Atlantic deep-sea sponge. Mar Drugs. 19(2), 105 (2021).
Robinson, S. L., Christenson, J. K. & Wackett, L. P. Biosynthesis and chemical diversity of beta-lactone natural products. Nat. Prod. Rep. 36(3), 458–475 (2019).
Wang, J., Shi, Y. & Jiang, D. β-lactone derivatives and their anticancer activities: a short review. Curr. Top Med. Chem. 21(18), 1645–1656 (2021).
Wang, J. et al. Pentaketide ansamycin microansamycins A-I from micromonospora sp reveal diverse post-PKS modifications. Org Lett. 20(4), 1058–1061 (2018).
Latos-Brozio, M. et al. Characteristics of the polyphenolic profile and antioxidant activity of cone extracts from conifers determined using electrochemical and spectrophotometric methods. Antioxidants (Basel). 10(11), 1723 (2021).
Wajs-Bonikowska, A. et al. Composition of essential oils from seeds of Abies koreana. Nat. Prod. Commun. 8(2), 227–230 (2013).
Wajs-Bonikowska, A. et al. Chemical composition and biological activity of Abies alba and A. koreana seed and cone essential oils and characterization of their seed hydrolates. Chem. Biodivers. 12(3), 407–418 (2015).
Zhang, D. et al. Endophytic bacillus altitudinis strain uses different novelty molecular pathways to enhance plant growth. Front Microbiol. 12, 692313 (2021).
Dahmani, M. A. et al. Unearthing the plant growth-promoting traits of Bacillus megaterium RmBm31, an endophytic bacterium isolated from root nodules of retama monosperma. Front Plant Sci. 11, 124 (2020).
Gammulla, C. G. et al. Differential metabolic response of cultured rice (Oryza sativa) cells exposed to high- and low-temperature stress. Proteomics 10(16), 3001–3019 (2010).
Cruz, RPd. et al. Avoiding damage and achieving cold tolerance in rice plants. Food Energy Secur. 2(2), 96–119 (2013).
Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33(7), 1870–1874 (2016).
Lee, Y. & Jeon, C. O. Cohnella algarum sp. nov, isolated from a freshwater green alga Paulinella chromatophora. Int. J. Syst. Evol. Microbiol. 67(11), 4767–4772 (2017).
Sasser, M., Bacterial identification by gas chromatographic analysis of fatty acids methyl esters (GC-FAME). MIDI Technical MIDI Technical Note #101. MIDI Inc, Newark, DE, USA. (2006).
Collins, M., Shah, H. & Minnikin, D. A note on the separation of natural mixtures of bacterial menaquinones using reverse phase thin-layer chromatography. J. Appl. Bacteriol. 48(2), 277–282 (1980).
Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27(5), 722–736 (2017).
Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9(11), e112963 (2014).
Lee, I. et al. ContEst16S: an algorithm that identifies contaminated prokaryotic genomes using 16S RNA gene sequences. Int. J. Syst. Evol Microbiol. 67(6), 2053–2057 (2017).
Na, S. I. et al. UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction. J. Microbiol. 56(4), 280–285 (2018).
Kanehisa, M., Sato, Y. & Morishima, K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol. 428(4), 726–731 (2016).
Blin, K. et al. antiSMASH 6.0: improving cluster detection and comparison capabilities. Nucleic Acids Res. 49(1), 29–35 (2021).
Meier-Kolthoff, J. P. et al. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 14, 60 (2013).
Lee, I. et al. OrthoANI: an improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 66(2), 1100–1103 (2016).
Fu, S. F. et al. Plant growth-promoting traits of yeasts isolated from the phyllosphere and rhizosphere of Drosera spatulata Lab. Fungal Biol. 120(3), 433–448 (2016).
Acknowledgements
This work was supported by the Korea Research Institute of Bioscience and Biotechnology (KRIBB) Research Initiative Program (KGM5282223), the National Institute of Ecology (NIE-B-2022-15), and the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2020R111A2072308).
Author information
Authors and Affiliations
Contributions
J.L. and H.C.P. conceptualized the study; J.L. and C.Y.K. acquired funding; L.J. designed the research and carried out the experiment; Y.P., J.S., D.J., M.G.J., J.H.L., J.C.J., and C.Y.K. edited the manuscript; L.J. wrote the manuscript; J.L. and L.J. modified the manuscript. All authors read and approved the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
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/.
About this article
Cite this article
Jiang, L., Peng, Y., Seo, J. et al. Subtercola endophyticus sp. nov., a cold-adapted bacterium isolated from Abies koreana. Sci Rep 12, 12114 (2022). https://doi.org/10.1038/s41598-022-16116-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-16116-3
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
