TEST0

PMC:4460873 / 21032-21151 JSON TXT

Unraveling adaptation of Pontibacter korlensis to radiation and infertility in desert through complete genome and comparative transcriptomic analysis Abstract The desert is a harsh habitat for flora and microbial life due to its aridness and strong radiation. In this study, we constructed the first complete and deeply annotated genome of the genus Pontibacter (Pontibacter korlensis X14-1T = CCTCC AB 206081T, X14-1). Reconstruction of the sugar metabolism process indicated that strain X14-1 can utilize diverse sugars, including cellulose, starch and sucrose; this result is consistent with previous experiments. Strain X14-1 is also able to resist desiccation and radiation in the desert through well-armed systems related to DNA repair, radical oxygen species (ROS) detoxification and the OstAB and TreYZ pathways for trehalose synthesis. A comparative transcriptomic analysis under gamma radiation revealed that strain X14-1 presents high-efficacy operating responses to radiation, including the robust expression of catalase and the manganese transport protein. Evaluation of 73 novel genes that are differentially expressed showed that some of these genes may contribute to the strain’s adaptation to radiation and desiccation through ferric transport and preservation. Approximately 10% of the Earth’s terrestrial surface is covered by desert with arid environments, which are characterized as environments with nutrient limitation, desiccation, cycles of extreme temperatures and intense radiation1. Nevertheless, diverse bacterial species have been identified and isolated from this extreme biotope2345 and have been found to be tolerant to solar radiation through various mechanisms, such as DNA repair, ROS detoxification and protein protection6. Since the radiation-resistant strain Deinococcus radiodurans R1 was isolated 50 years ago, studies of bacterial resistance and tolerance to solar radiation have been mainly performed on the genus Deinococcus7. D. radiodurans is 200-fold and 20-fold more resistant to ionizing radiation and UV irradiation, respectively, than Escherichia coli8, and the complete genome of D. radiodurans was first published in 19999. To elucidate the extreme resistance phenotype of D. radiodurans R1, various research strategies have been combined10, and three hypotheses regarding DNA repair have been proposed11. The lack of novelty in DNA repair-related genes/proteins and the greater efficiency of specific bacteria to use conventional repair pathways are partially supported by the findings from previous studies121314151617. An in-depth analysis of the D. radiodurans R1 genome and its gene expression profile revealed that many undefined genes, including ddrA, ddrB, ddrC, ddrD and pprA, are involved in DNA repair18192021, suggesting that repair functions are encoded by these hypothetical genes. The last hypothesis is that ring-like nucleoids (RNs) contribute to DNA repair22. There are also three assumptions regarding the maintenance of a low ROS concentration in bacteria10, most of which are detoxifying and scavenging ROS, including small catalase, superoxide dismutase, and antioxidant molecules, and exhibit an increased Mn(II)/Fe ratio intermediated by manganese complexes111323. Daly and Krisko found that molecules smaller than 3 kDa in the extract of Deinococcus radiodurans R1 can impose antioxidant protection on E. coli proteins2324. The promotion of metabolic activities with decreased ROS production (e.g., glyoxylate bypass of the TCA cycle2125) is an alternative to the response to oxidative damage and single antioxidant pathways through high ROS production, which could be inactivated due to redundant ROS-tolerance mechanisms. To lower the ROS, it is also helpful to reduce proteins with Fe-S clusters and the number of respiratory chain enzymes25. In addition to maintaining a low ROS concentration, many other metabolic activities, such as proteolysis and glucose metabolism, contribute to the robustness of D. radiodurans R115232627. In addition to D. radiodurans R1, additional genome sequences of the genus Deinococcus have been published262829303132, and comparative analyses have been performed to elucidate the diverse molecular mechanisms and physiological determinants underlying the extreme resistance phenotype3334. We isolated the strain Pontibacter korlensis X14-1T (X14-1) from the surface layer of a desert in Xinjiang, China, and identified it as a new species of the genus Pontibacter2. This study provides the first complete genome of the genus Pontibacter and attempted to delineate genomic components related to radiation and desiccation resistance in comparison with other species from the genus Pontibacter. A comparative analysis of the gene expression profile under radiation was also conducted to unravel the complicated mechanisms of strain X14-1 involved in its adaptation to the arid environment of the desert. This work will provide referable information for the comprehensive understanding of the evolution and adaptation of the genus Pontibacter as well as various radiation and desiccation resistances. Results Genomic characteristics and phylogeny of strain X14-1 The complete genome sequence of strain X14-1 was produced based on high-quality reads and corrected by read mapping and PCR verification. Strain X14-1 has a larger genome size (5.46 MB) and a lower GC content (47.3%) than three other Pontibacter strains distributed in different species (summarized in Table 1). We found that most of the transposase-related genes are located near genomic islands and next to DNA repair- and ROS detoxification-related genes (Fig. 1), which implied that mobile genetic elements (MGEs) play an important role in the adaptation to radiation and desiccation in the desert. Differences in genome size and MGEs between strain X14-1 and other Pontibacter strains may be attributed to the genomic evolution or gapped assembly of P. actiniarum DSM 19842, Pontibacter sp. BAB1700 and P. roseus DSM 17521. To confirm the phylotype of strain X14-1, we downloaded 40 genomes of the family Cytophagaceae (higher taxonomic classification of the genus Pontibacter) available in the NCBI database. Phylogeny analysis indicated the same results as those previously reported based on the 16S rDNA sequence2, and P. actiniarum DSM 19842 was found to be the most homologous to strain X14-1 (Fig. 2), a finding that is also supported by the following functional analysis. Sugar metabolism in strain X14-1 In comparison with three other genomes from the genus Pontibacter, we found that only strain X14-1 harbors genes encoding D-fructokinase, which is essential for sucrose and fructose utilization, and this is consistent with previous experimental results2. Beta-galactosidase, which is essential for strain X14-1 to use lactose as an alternative carbon source by catalyzing lactose to galactose and glucose, is specific to strain X14-1 compared with other Pontibacter strains. Although comparative analysis revealed a common dispersion of cellobiose glucohydrolase in Pontibacter, enzymes responsible for degrading cellulose to cellobiose are only distributed in strain X14-1. Starch could be degraded to amylose and alpha-D-galactose-1-phosphate, which is an intermediate in the production of UDP-glucose that could link pentose and glucuronate interconversion. This is important for the utilization of D-galactose as a carbon resource. Mannose can enter glycolysis through beta-D-fructose-6-phosphate with the help of hexokinase and mannose-6-phosphate isomerase, which could be encoded by genes in strain X14-1. The ability to utilize versatile sugars as described above (Fig. 3) could partly explain how strain X14-1 survives in an infertile desert. Determinants in the genome for the adaptation of strain X14-1 to the radiation and aridness of the desert Genes related to DNA repair and the stress response, including ROS detoxification and the osmotic response, were analyzed to understand whether and how strain X14-1 protect itself against to desiccation and radiation. In strain X14-1, recombination repair-related genes are the most abundant, followed by base excision repair (BER) and nucleotide excision repair (NER) (Fig. 4, detailed in Supplementary File 1). There are no specific corresponding genetic determinants for the radiation resistance of strain X14-1, implying that the robustness of the resistance of strain X14-1 to radiation and desiccation is not attributed to new genes but to the high-efficacy operation of systems relevant to recovery from radiation. Trehalose is a natural product that can form a protective film outside the cell under low and high temperatures, osmotic pressure and aridness, preventing proteins from becoming inactivated. The analyses conducted in this study elucidated that intermediate products from the metabolism of several sugars, such as starch and lactose, in strain X14-1 could be transferred to trehalose through the OstAB and TreYZ pathways, as summarized in a previous publication35 (Fig. 5). Trehalase-encoding genes, which are homologous to the trehalase-like protein-encoding genes of Gramella forsetii KT0803, were also identified, suggesting the trehalose consumption of strain X 14-1. Responses of strain X14-1 to gamma radiation revealed by RNA-Seq After radiation, the average survival rate of strain X14-1 ranged from 43.4% to 55.1%, as determined from triplicate analyses. To better understand the specific responses of strain X14-1 to gamma radiation, we determined the expression level of all possible determinants responsible for radiation adaptation and ROS detoxification. BER, NER and homologous recombination-related genes, which are mainly required for recovery from gamma radiation, are apparently upregulated (Supplementary File 2). Catalase-encoding genes exhibit intensive upregulation after radiation, but a similar finding was not obtained for superoxide dismutase, which is also important for ROS detoxification (Supplementary File 2). The expression of the manganese transport protein is 10-fold higher after radiation, and a similar trend was identified for subunit A-F of the multicomponent Na+:H+ antiporter as well as other ion-coupled transporters (Supplementary File 2). NADH dehydrogenase I, with the exception of except subunit C, exhibited robust expression in the radiated samples compared with the controls. Ferric uptake regulator and LysR family transcriptional regulator, both of which are transcription factors and have been identified as determinants for oxidative detoxification in the genome, also exhibited increased expression after radiation (Supplementary File 2). Initial evaluation of unannotated genes with differential expression Apparent differential expression was found for 394 of 893 unclassified genes as a result of gamma radiation. To explore whether these unknown genes contribute to the adaptation of strain X14-1 to radiation, an alignment of these genes to the conserved domain database (CDD) was conducted. Of the 394 new genes, 73 could be assigned to specific domains in CDD, and the corresponding hypothesized functions are summarized in Supplementary File 3, indicating that these 73 genes are commonly found in the FecR, DUF and FlgD_ig (associated with ferric transport and preservation) superfamily. Discussion The desert is characterized by a low diversity of inhabitants due to fluctuations in temperature, aridness, strong radiation and poor nutrition136. Elucidating the microbial life in the desert could be important because it could help us understand the dry limit of life and improve our research on extremophiles37. Contributors to radiation resistance and adaptation to poor nutrition in microbes could be applied to improve plant culturing in the desert. This will decrease sand storms mainly induced by the desert and increase arable lands to aid the avoidance of hunger worldwide. Determinants to radiation resistance in microbes residing in the desert could also be industrially applied in cosmetic and soil remediation under radiation. Strain X14-1, which was isolated from the sand surface in the desert, belongs to the poorly understood genus Pontibacter, and there are no publications on an in-depth analysis of the Pontibacter genome. In this study, we performed a high-coverage sequencing of strain X14-1 and produced the first complete map, which may be used as a valuable reference to promote research on the genus Pontibacter. A detailed annotation of the strain X14-1 genome in combination with a comparative transcriptomic analysis was conducted to understand how strain X14-1 utilizes diverse alternative nutrition sources and recovers from desiccation and strong radiation. Previous studies have revealed that strain X14-1 can utilize versatile carbon sources2, and this finding is supported by the reconstruction of the metabolic pathways of different sugars. D-fructokinase, beta-galactosidase, endoglucanase and beta-glucosidase are key for the utilization of sucrose, fructose, lactose, cellulose and glycoside, respectively. Enzymes that are important for metabolizing starch, galactose and mannose were also found in strain X14-1. In addition to their robust competence for utilizing diverse available sugars, another interesting finding is that only strain X14-1 and P. actiniarum DSM19842 possess phosphoenolpyruvate (PEP) carboxylase, which is key for CO2 fixation, mainly in plants, through the C4 cycle and may be used to resemble the autotrophic lifestyle in environments with extremely poor nutrition. The fact that most of the hypothetical proteins in strain X14-1 are homologous to those in P. actiniarum DSM19842 also demonstrates the homology of these two Pontibacter strains revealed by the phylogenetic construction (Fig. 2). In addition to the various metabolic pathways required for the utilization of limited resources, strain X14-1 also harbors an intensive arsenal of DNA repair- and stress response-related determinants that allow it to survive in the desert. An in-depth analysis of the genomic components of strain X14-1 revealed possible broad-spectrum contributors to resistance to radiation (Supplementary File 1). This is consistent with the genetic factors found in the genus Deinococcus, and we also found new genes that encode resistance agents, indicating supplementary mechanisms to those found in Deinococcus. Based on changes in the annotated genes after gamma radiation, we found that catalases but not dismutase exhibited a marked increase, indicating the robustness of catalase in oxidative detoxification in strain X14-1. An increased Mn(II)/Fe ratio has been demonstrated to play an important role in detoxification from oxidative damage11, and this is also supported by our findings: upregulation of the manganese transport protein, ferric uptake regulator and LysR family transcriptional regulator. The expression of gene gshB encoding glutathione synthase was high during radiation, but glutathione peroxidase and reductase, which are responsible for the circulation of GSS (reduced glutathione) and GSSG (oxidized glutathione) exhibited no change, implying that glutathione is not the intermediate in strain X14-1during radiation resistance. An interesting finding is the robustly elevated expression of transposase during gamma radiation (Supplementary File 2), which implies the potently active role of widely distributed transposases in strain X14-1 for desert adaptation. Some of the conserved domains of upregulated unknown genes could be assigned to the FecR and ferritin-like superfamily. This finding may imply that some unknown genes could contribute to oxidative detoxification because these proteins are related with ferric transport or preservation, but further analysis or experiments are required to confirm this hypothesis. Materials and Methods Strains and culture conditions All of the strains were obtained from the China Center for Type Culture Collection (CCTCC). P. korlensis X14-1 was grown at 30 °C on marine broth 2216 (Difco). To determine the tolerance of the culture to gamma radiation, the strains were grown in the appropriate liquid medium to the exponential phase, and 2.5 ml of the cell culture was then subjected to 60Co radiation with a continuous dose rate of 1000 Gray/h at room temperature. The control culture cells were not treated with radiation. Samples were collected 3 h after culture with (radiated samples) or without (controls) radiation treatment. Biological triplicates were harvested from each treatment. These samples were immediately frozen in RNAlater (Qiagen, German) and stored at –80 °C for RNA analysis. DNA extraction, sequencing and assembly Total DNA was extracted through the traditional CTAB method, and its quality was tested by Qubit. Sequencing was performed for a 470-bp insert size library (1.65-Gb clean reads with a length of 90 bp) and mate-paired libraries with insert sizes of 2222 bp and 6200 bp (477-Mb clean reads, each with a length of 49 bp), respectively. Clean reads were used to assemble the complete genome sequence as previously described38. Gene prediction, annotation and comparative analysis Genes were predicted based on the assembled and confirmed sequence using GLIMMER394041, which was developed for microorganisms including bacteria, archaea, and viruses. This was also used for the prediction of genes from other three Pontibacter strains downloaded from the NCBI ftp site. The annotations were conducted by assigning the predicted genes to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/)42, Cluster of Orthologous Groups of proteins (COG)4344, and the protein database in NCBI (http://www.ncbi.nlm.nih.gov/protein/) through a BLASTP search (e-value ≤1e-5, query coverage ≥40%). The strategy described by Tian et al.45 was used to identify the orthologous proteins in four Pontibacter genomes. In the phylogenetic analysis, homologous genes from 41 strains were used to construct the gene family and phylogenetic tree: (i) the gene set was first generated using BLAST (V2.2.23) with the genes from 41 strains; (ii) the gene family was constructed using OrthoMCL (V1.4) with the gene set data; (iii) to construct the phylogenetic tree, single-copy genes in the gene family were screened, and the genes were aligned according to the strains using Muscle (V3.8.31); and (IV) the phylogenetic tree was drawn using TreeBeST (V1.9.2) with 1000 bootstrap replicates. Identification of MGEs The prophages in the genome sequences were identified using the Phage Search Tool (PHAST)46. Genomic islands were predicted with multiple methods (IslandPath-DIOMB, SIGI-HMM, and IslandPicker using Islandviewer474849. We searched repeat sequences with the RepeatMasker program50, and transposases were selected from the annotation results. RNA-Seq analysis Three replicates of each group were used for total RNA extraction and mRNA enrichment using the Ribo-ZeroTM kit. The qualified mRNA was sequenced on the HiSeq2000 platform, and the expression level of each gene was calculated using the reads per kilo bases per million reads (RPKM)51. The genes that were differentially expressed between the two groups were identified as described by Audicet et al.52. The analysis results with a p-value <0.05 were corrected by the FDR (false discovery rate), which was set to ≤0.001. Additional Information How to cite this article: Dai, J. et al. Unraveling adaptation of Pontibacter korlensis to radiation and infertility in desert through complete genome and comparative transcriptomic analysis. Sci. Rep. 5, 10929; doi: 10.1038/srep10929 (2015). Supplementary Material Supplementary File 1 Supplementary File 2 Supplementary File 3 We would like to thank Xiaohui Li, Hongjuan Wang and Yongfeng Liu for coordination during this project. This research was supported by the National Natural Science Foundation of China (No: 31300003), the Hubei Provincial Natural Science Foundation (2015CFB679) and Research Fund for the Doctoral Program of Hubei University of Technology (No: BSQD13003). Author Contributions X.C. managed the project. C.F. and Z.W. offered the strain and designed the experiments. Y.Z., L.Z., Q.Y., X.L., Z.W., M.Z. and Z.J. performed experiments and prepared the DNA and RNA samples. Z.Y. and Q.G. performed sequencing and assembly. C.Q. and W.D. performed other bioinformatics analysis in this work. W.D., C.Q. and J.D. wrote the paper. All authors reviewed this manuscript. Accession codes: Assembly result has been deposited in GenBank nucleotide core database under the accession number CP009621. Clean reads for six transcriptome samples could be found through accession number SRP049974. Figure 1 Genomic components of P. korlensis X14-1. Distribution of genes, genomic islands, transposases, DNA repair and oxidative response related determinants in the genome of strain X14-1. It indicates that strain X14-1 harbors abundant MGEs, suggestive of high genome plasticity. Figure 2 Phylogenetic tree of family Cytophagaceae based on 40 available reference genomes and strain X14-1(highlighted in red). The phylogenetic tree was produced on single copy gene within gene family and drawn by TreeBeST version1.9.2 under bootstrap 1000, implicating taxonomic position of strain X14-1. Figure 3 Pathways of sugar utilization for strain X14-1. PGI: Pentose and Glucuronate Interconversion; PPP: Pentose Phosphate Pathway. Arrows with different colors represent utilizing distinct carbon sources: Blue-lactose and starch, orange-glycoside, gray-sucrose, brown-cellulose, black-glucose, red-galactose, green-mannose, purple-fructose. Figure 4 Genes related to DNA repair and stress response in genus Pontibacter. BER: base excision repair; NER: nucleotide excision repair; DR: direct reversal; RER: recombination repair; ODR: other DNA repair; ORT: other radiation tolerance. Figure 5 Trehalose synthesis of strain X14-1. UDP-Glucose from lactose and galactose (Fig. 3) could be catalyzed to trehalose through trehalose 6-phosphate synthase and trehalose 6-phosphate phosphatase. Glycogen operon protein could degrade starch synthesized from alpha-D-Glucose-6P available from many metabolic pathways to trehalose. Table 1 Comparison of strain X14-1 with other species in genus Pontibacter. P. korlensis X14-1 P. actiniarum DSM19842 Pontibacter sp. BAB1700 P. roseus DSM17521 Habitat Desert Aquatic Multiple NA Genome size(Mb) 5.46 4.95 4.54 4.58 GC content(%) 47.3 53.1 50.0 52.6 Total genes 5,037 4,689 4,849 4,260 Coding regions(%) 86.58 86.42 84.97 87.01 Unannotated genes 893 901 507 328 Insertion sequence 4 6 1 0 Prophage 1 0 0 0 Transposase/Integrase 114 28 10 5

Document structure show

article-title	Unraveling adaptation of Pontibacter korlensis to radiation and infertility in desert through complete genome and comparative transcriptomic analysis
abstract	The desert is a harsh habitat for flora and microbial life due to its aridness and strong radiation. In this study, we constructed the first complete and deeply annotated genome of the genus Pontibacter (Pontibacter korlensis X14-1T = CCTCC AB 206081T, X14-1). Reconstruction of the sugar metabolism process indicated that strain X14-1 can utilize diverse sugars, including cellulose, starch and sucrose; this result is consistent with previous experiments. Strain X14-1 is also able to resist desiccation and radiation in the desert through well-armed systems related to DNA repair, radical oxygen species (ROS) detoxification and the OstAB and TreYZ pathways for trehalose synthesis. A comparative transcriptomic analysis under gamma radiation revealed that strain X14-1 presents high-efficacy operating responses to radiation, including the robust expression of catalase and the manganese transport protein. Evaluation of 73 novel genes that are differentially expressed showed that some of these genes may contribute to the strain’s adaptation to radiation and desiccation through ferric transport and preservation.
p	The desert is a harsh habitat for flora and microbial life due to its aridness and strong radiation. In this study, we constructed the first complete and deeply annotated genome of the genus Pontibacter (Pontibacter korlensis X14-1T = CCTCC AB 206081T, X14-1). Reconstruction of the sugar metabolism process indicated that strain X14-1 can utilize diverse sugars, including cellulose, starch and sucrose; this result is consistent with previous experiments. Strain X14-1 is also able to resist desiccation and radiation in the desert through well-armed systems related to DNA repair, radical oxygen species (ROS) detoxification and the OstAB and TreYZ pathways for trehalose synthesis. A comparative transcriptomic analysis under gamma radiation revealed that strain X14-1 presents high-efficacy operating responses to radiation, including the robust expression of catalase and the manganese transport protein. Evaluation of 73 novel genes that are differentially expressed showed that some of these genes may contribute to the strain’s adaptation to radiation and desiccation through ferric transport and preservation.
body	Approximately 10% of the Earth’s terrestrial surface is covered by desert with arid environments, which are characterized as environments with nutrient limitation, desiccation, cycles of extreme temperatures and intense radiation1. Nevertheless, diverse bacterial species have been identified and isolated from this extreme biotope2345 and have been found to be tolerant to solar radiation through various mechanisms, such as DNA repair, ROS detoxification and protein protection6. Since the radiation-resistant strain Deinococcus radiodurans R1 was isolated 50 years ago, studies of bacterial resistance and tolerance to solar radiation have been mainly performed on the genus Deinococcus7. D. radiodurans is 200-fold and 20-fold more resistant to ionizing radiation and UV irradiation, respectively, than Escherichia coli8, and the complete genome of D. radiodurans was first published in 19999. To elucidate the extreme resistance phenotype of D. radiodurans R1, various research strategies have been combined10, and three hypotheses regarding DNA repair have been proposed11. The lack of novelty in DNA repair-related genes/proteins and the greater efficiency of specific bacteria to use conventional repair pathways are partially supported by the findings from previous studies121314151617. An in-depth analysis of the D. radiodurans R1 genome and its gene expression profile revealed that many undefined genes, including ddrA, ddrB, ddrC, ddrD and pprA, are involved in DNA repair18192021, suggesting that repair functions are encoded by these hypothetical genes. The last hypothesis is that ring-like nucleoids (RNs) contribute to DNA repair22. There are also three assumptions regarding the maintenance of a low ROS concentration in bacteria10, most of which are detoxifying and scavenging ROS, including small catalase, superoxide dismutase, and antioxidant molecules, and exhibit an increased Mn(II)/Fe ratio intermediated by manganese complexes111323. Daly and Krisko found that molecules smaller than 3 kDa in the extract of Deinococcus radiodurans R1 can impose antioxidant protection on E. coli proteins2324. The promotion of metabolic activities with decreased ROS production (e.g., glyoxylate bypass of the TCA cycle2125) is an alternative to the response to oxidative damage and single antioxidant pathways through high ROS production, which could be inactivated due to redundant ROS-tolerance mechanisms. To lower the ROS, it is also helpful to reduce proteins with Fe-S clusters and the number of respiratory chain enzymes25. In addition to maintaining a low ROS concentration, many other metabolic activities, such as proteolysis and glucose metabolism, contribute to the robustness of D. radiodurans R115232627. In addition to D. radiodurans R1, additional genome sequences of the genus Deinococcus have been published262829303132, and comparative analyses have been performed to elucidate the diverse molecular mechanisms and physiological determinants underlying the extreme resistance phenotype3334. We isolated the strain Pontibacter korlensis X14-1T (X14-1) from the surface layer of a desert in Xinjiang, China, and identified it as a new species of the genus Pontibacter2. This study provides the first complete genome of the genus Pontibacter and attempted to delineate genomic components related to radiation and desiccation resistance in comparison with other species from the genus Pontibacter. A comparative analysis of the gene expression profile under radiation was also conducted to unravel the complicated mechanisms of strain X14-1 involved in its adaptation to the arid environment of the desert. This work will provide referable information for the comprehensive understanding of the evolution and adaptation of the genus Pontibacter as well as various radiation and desiccation resistances. Results Genomic characteristics and phylogeny of strain X14-1 The complete genome sequence of strain X14-1 was produced based on high-quality reads and corrected by read mapping and PCR verification. Strain X14-1 has a larger genome size (5.46 MB) and a lower GC content (47.3%) than three other Pontibacter strains distributed in different species (summarized in Table 1). We found that most of the transposase-related genes are located near genomic islands and next to DNA repair- and ROS detoxification-related genes (Fig. 1), which implied that mobile genetic elements (MGEs) play an important role in the adaptation to radiation and desiccation in the desert. Differences in genome size and MGEs between strain X14-1 and other Pontibacter strains may be attributed to the genomic evolution or gapped assembly of P. actiniarum DSM 19842, Pontibacter sp. BAB1700 and P. roseus DSM 17521. To confirm the phylotype of strain X14-1, we downloaded 40 genomes of the family Cytophagaceae (higher taxonomic classification of the genus Pontibacter) available in the NCBI database. Phylogeny analysis indicated the same results as those previously reported based on the 16S rDNA sequence2, and P. actiniarum DSM 19842 was found to be the most homologous to strain X14-1 (Fig. 2), a finding that is also supported by the following functional analysis. Sugar metabolism in strain X14-1 In comparison with three other genomes from the genus Pontibacter, we found that only strain X14-1 harbors genes encoding D-fructokinase, which is essential for sucrose and fructose utilization, and this is consistent with previous experimental results2. Beta-galactosidase, which is essential for strain X14-1 to use lactose as an alternative carbon source by catalyzing lactose to galactose and glucose, is specific to strain X14-1 compared with other Pontibacter strains. Although comparative analysis revealed a common dispersion of cellobiose glucohydrolase in Pontibacter, enzymes responsible for degrading cellulose to cellobiose are only distributed in strain X14-1. Starch could be degraded to amylose and alpha-D-galactose-1-phosphate, which is an intermediate in the production of UDP-glucose that could link pentose and glucuronate interconversion. This is important for the utilization of D-galactose as a carbon resource. Mannose can enter glycolysis through beta-D-fructose-6-phosphate with the help of hexokinase and mannose-6-phosphate isomerase, which could be encoded by genes in strain X14-1. The ability to utilize versatile sugars as described above (Fig. 3) could partly explain how strain X14-1 survives in an infertile desert. Determinants in the genome for the adaptation of strain X14-1 to the radiation and aridness of the desert Genes related to DNA repair and the stress response, including ROS detoxification and the osmotic response, were analyzed to understand whether and how strain X14-1 protect itself against to desiccation and radiation. In strain X14-1, recombination repair-related genes are the most abundant, followed by base excision repair (BER) and nucleotide excision repair (NER) (Fig. 4, detailed in Supplementary File 1). There are no specific corresponding genetic determinants for the radiation resistance of strain X14-1, implying that the robustness of the resistance of strain X14-1 to radiation and desiccation is not attributed to new genes but to the high-efficacy operation of systems relevant to recovery from radiation. Trehalose is a natural product that can form a protective film outside the cell under low and high temperatures, osmotic pressure and aridness, preventing proteins from becoming inactivated. The analyses conducted in this study elucidated that intermediate products from the metabolism of several sugars, such as starch and lactose, in strain X14-1 could be transferred to trehalose through the OstAB and TreYZ pathways, as summarized in a previous publication35 (Fig. 5). Trehalase-encoding genes, which are homologous to the trehalase-like protein-encoding genes of Gramella forsetii KT0803, were also identified, suggesting the trehalose consumption of strain X 14-1. Responses of strain X14-1 to gamma radiation revealed by RNA-Seq After radiation, the average survival rate of strain X14-1 ranged from 43.4% to 55.1%, as determined from triplicate analyses. To better understand the specific responses of strain X14-1 to gamma radiation, we determined the expression level of all possible determinants responsible for radiation adaptation and ROS detoxification. BER, NER and homologous recombination-related genes, which are mainly required for recovery from gamma radiation, are apparently upregulated (Supplementary File 2). Catalase-encoding genes exhibit intensive upregulation after radiation, but a similar finding was not obtained for superoxide dismutase, which is also important for ROS detoxification (Supplementary File 2). The expression of the manganese transport protein is 10-fold higher after radiation, and a similar trend was identified for subunit A-F of the multicomponent Na+:H+ antiporter as well as other ion-coupled transporters (Supplementary File 2). NADH dehydrogenase I, with the exception of except subunit C, exhibited robust expression in the radiated samples compared with the controls. Ferric uptake regulator and LysR family transcriptional regulator, both of which are transcription factors and have been identified as determinants for oxidative detoxification in the genome, also exhibited increased expression after radiation (Supplementary File 2). Initial evaluation of unannotated genes with differential expression Apparent differential expression was found for 394 of 893 unclassified genes as a result of gamma radiation. To explore whether these unknown genes contribute to the adaptation of strain X14-1 to radiation, an alignment of these genes to the conserved domain database (CDD) was conducted. Of the 394 new genes, 73 could be assigned to specific domains in CDD, and the corresponding hypothesized functions are summarized in Supplementary File 3, indicating that these 73 genes are commonly found in the FecR, DUF and FlgD_ig (associated with ferric transport and preservation) superfamily. Discussion The desert is characterized by a low diversity of inhabitants due to fluctuations in temperature, aridness, strong radiation and poor nutrition136. Elucidating the microbial life in the desert could be important because it could help us understand the dry limit of life and improve our research on extremophiles37. Contributors to radiation resistance and adaptation to poor nutrition in microbes could be applied to improve plant culturing in the desert. This will decrease sand storms mainly induced by the desert and increase arable lands to aid the avoidance of hunger worldwide. Determinants to radiation resistance in microbes residing in the desert could also be industrially applied in cosmetic and soil remediation under radiation. Strain X14-1, which was isolated from the sand surface in the desert, belongs to the poorly understood genus Pontibacter, and there are no publications on an in-depth analysis of the Pontibacter genome. In this study, we performed a high-coverage sequencing of strain X14-1 and produced the first complete map, which may be used as a valuable reference to promote research on the genus Pontibacter. A detailed annotation of the strain X14-1 genome in combination with a comparative transcriptomic analysis was conducted to understand how strain X14-1 utilizes diverse alternative nutrition sources and recovers from desiccation and strong radiation. Previous studies have revealed that strain X14-1 can utilize versatile carbon sources2, and this finding is supported by the reconstruction of the metabolic pathways of different sugars. D-fructokinase, beta-galactosidase, endoglucanase and beta-glucosidase are key for the utilization of sucrose, fructose, lactose, cellulose and glycoside, respectively. Enzymes that are important for metabolizing starch, galactose and mannose were also found in strain X14-1. In addition to their robust competence for utilizing diverse available sugars, another interesting finding is that only strain X14-1 and P. actiniarum DSM19842 possess phosphoenolpyruvate (PEP) carboxylase, which is key for CO2 fixation, mainly in plants, through the C4 cycle and may be used to resemble the autotrophic lifestyle in environments with extremely poor nutrition. The fact that most of the hypothetical proteins in strain X14-1 are homologous to those in P. actiniarum DSM19842 also demonstrates the homology of these two Pontibacter strains revealed by the phylogenetic construction (Fig. 2). In addition to the various metabolic pathways required for the utilization of limited resources, strain X14-1 also harbors an intensive arsenal of DNA repair- and stress response-related determinants that allow it to survive in the desert. An in-depth analysis of the genomic components of strain X14-1 revealed possible broad-spectrum contributors to resistance to radiation (Supplementary File 1). This is consistent with the genetic factors found in the genus Deinococcus, and we also found new genes that encode resistance agents, indicating supplementary mechanisms to those found in Deinococcus. Based on changes in the annotated genes after gamma radiation, we found that catalases but not dismutase exhibited a marked increase, indicating the robustness of catalase in oxidative detoxification in strain X14-1. An increased Mn(II)/Fe ratio has been demonstrated to play an important role in detoxification from oxidative damage11, and this is also supported by our findings: upregulation of the manganese transport protein, ferric uptake regulator and LysR family transcriptional regulator. The expression of gene gshB encoding glutathione synthase was high during radiation, but glutathione peroxidase and reductase, which are responsible for the circulation of GSS (reduced glutathione) and GSSG (oxidized glutathione) exhibited no change, implying that glutathione is not the intermediate in strain X14-1during radiation resistance. An interesting finding is the robustly elevated expression of transposase during gamma radiation (Supplementary File 2), which implies the potently active role of widely distributed transposases in strain X14-1 for desert adaptation. Some of the conserved domains of upregulated unknown genes could be assigned to the FecR and ferritin-like superfamily. This finding may imply that some unknown genes could contribute to oxidative detoxification because these proteins are related with ferric transport or preservation, but further analysis or experiments are required to confirm this hypothesis. Materials and Methods Strains and culture conditions All of the strains were obtained from the China Center for Type Culture Collection (CCTCC). P. korlensis X14-1 was grown at 30 °C on marine broth 2216 (Difco). To determine the tolerance of the culture to gamma radiation, the strains were grown in the appropriate liquid medium to the exponential phase, and 2.5 ml of the cell culture was then subjected to 60Co radiation with a continuous dose rate of 1000 Gray/h at room temperature. The control culture cells were not treated with radiation. Samples were collected 3 h after culture with (radiated samples) or without (controls) radiation treatment. Biological triplicates were harvested from each treatment. These samples were immediately frozen in RNAlater (Qiagen, German) and stored at –80 °C for RNA analysis. DNA extraction, sequencing and assembly Total DNA was extracted through the traditional CTAB method, and its quality was tested by Qubit. Sequencing was performed for a 470-bp insert size library (1.65-Gb clean reads with a length of 90 bp) and mate-paired libraries with insert sizes of 2222 bp and 6200 bp (477-Mb clean reads, each with a length of 49 bp), respectively. Clean reads were used to assemble the complete genome sequence as previously described38. Gene prediction, annotation and comparative analysis Genes were predicted based on the assembled and confirmed sequence using GLIMMER394041, which was developed for microorganisms including bacteria, archaea, and viruses. This was also used for the prediction of genes from other three Pontibacter strains downloaded from the NCBI ftp site. The annotations were conducted by assigning the predicted genes to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/)42, Cluster of Orthologous Groups of proteins (COG)4344, and the protein database in NCBI (http://www.ncbi.nlm.nih.gov/protein/) through a BLASTP search (e-value ≤1e-5, query coverage ≥40%). The strategy described by Tian et al.45 was used to identify the orthologous proteins in four Pontibacter genomes. In the phylogenetic analysis, homologous genes from 41 strains were used to construct the gene family and phylogenetic tree: (i) the gene set was first generated using BLAST (V2.2.23) with the genes from 41 strains; (ii) the gene family was constructed using OrthoMCL (V1.4) with the gene set data; (iii) to construct the phylogenetic tree, single-copy genes in the gene family were screened, and the genes were aligned according to the strains using Muscle (V3.8.31); and (IV) the phylogenetic tree was drawn using TreeBeST (V1.9.2) with 1000 bootstrap replicates. Identification of MGEs The prophages in the genome sequences were identified using the Phage Search Tool (PHAST)46. Genomic islands were predicted with multiple methods (IslandPath-DIOMB, SIGI-HMM, and IslandPicker using Islandviewer474849. We searched repeat sequences with the RepeatMasker program50, and transposases were selected from the annotation results. RNA-Seq analysis Three replicates of each group were used for total RNA extraction and mRNA enrichment using the Ribo-ZeroTM kit. The qualified mRNA was sequenced on the HiSeq2000 platform, and the expression level of each gene was calculated using the reads per kilo bases per million reads (RPKM)51. The genes that were differentially expressed between the two groups were identified as described by Audicet et al.52. The analysis results with a p-value <0.05 were corrected by the FDR (false discovery rate), which was set to ≤0.001. Additional Information How to cite this article: Dai, J. et al. Unraveling adaptation of Pontibacter korlensis to radiation and infertility in desert through complete genome and comparative transcriptomic analysis. Sci. Rep. 5, 10929; doi: 10.1038/srep10929 (2015). Supplementary Material Supplementary File 1 Supplementary File 2 Supplementary File 3
p	Approximately 10% of the Earth’s terrestrial surface is covered by desert with arid environments, which are characterized as environments with nutrient limitation, desiccation, cycles of extreme temperatures and intense radiation1. Nevertheless, diverse bacterial species have been identified and isolated from this extreme biotope2345 and have been found to be tolerant to solar radiation through various mechanisms, such as DNA repair, ROS detoxification and protein protection6.
p	Since the radiation-resistant strain Deinococcus radiodurans R1 was isolated 50 years ago, studies of bacterial resistance and tolerance to solar radiation have been mainly performed on the genus Deinococcus7. D. radiodurans is 200-fold and 20-fold more resistant to ionizing radiation and UV irradiation, respectively, than Escherichia coli8, and the complete genome of D. radiodurans was first published in 19999. To elucidate the extreme resistance phenotype of D. radiodurans R1, various research strategies have been combined10, and three hypotheses regarding DNA repair have been proposed11. The lack of novelty in DNA repair-related genes/proteins and the greater efficiency of specific bacteria to use conventional repair pathways are partially supported by the findings from previous studies121314151617. An in-depth analysis of the D. radiodurans R1 genome and its gene expression profile revealed that many undefined genes, including ddrA, ddrB, ddrC, ddrD and pprA, are involved in DNA repair18192021, suggesting that repair functions are encoded by these hypothetical genes. The last hypothesis is that ring-like nucleoids (RNs) contribute to DNA repair22.
p	There are also three assumptions regarding the maintenance of a low ROS concentration in bacteria10, most of which are detoxifying and scavenging ROS, including small catalase, superoxide dismutase, and antioxidant molecules, and exhibit an increased Mn(II)/Fe ratio intermediated by manganese complexes111323. Daly and Krisko found that molecules smaller than 3 kDa in the extract of Deinococcus radiodurans R1 can impose antioxidant protection on E. coli proteins2324. The promotion of metabolic activities with decreased ROS production (e.g., glyoxylate bypass of the TCA cycle2125) is an alternative to the response to oxidative damage and single antioxidant pathways through high ROS production, which could be inactivated due to redundant ROS-tolerance mechanisms. To lower the ROS, it is also helpful to reduce proteins with Fe-S clusters and the number of respiratory chain enzymes25. In addition to maintaining a low ROS concentration, many other metabolic activities, such as proteolysis and glucose metabolism, contribute to the robustness of D. radiodurans R115232627. In addition to D. radiodurans R1, additional genome sequences of the genus Deinococcus have been published262829303132, and comparative analyses have been performed to elucidate the diverse molecular mechanisms and physiological determinants underlying the extreme resistance phenotype3334.
p	We isolated the strain Pontibacter korlensis X14-1T (X14-1) from the surface layer of a desert in Xinjiang, China, and identified it as a new species of the genus Pontibacter2. This study provides the first complete genome of the genus Pontibacter and attempted to delineate genomic components related to radiation and desiccation resistance in comparison with other species from the genus Pontibacter. A comparative analysis of the gene expression profile under radiation was also conducted to unravel the complicated mechanisms of strain X14-1 involved in its adaptation to the arid environment of the desert. This work will provide referable information for the comprehensive understanding of the evolution and adaptation of the genus Pontibacter as well as various radiation and desiccation resistances.
sec	Results Genomic characteristics and phylogeny of strain X14-1 The complete genome sequence of strain X14-1 was produced based on high-quality reads and corrected by read mapping and PCR verification. Strain X14-1 has a larger genome size (5.46 MB) and a lower GC content (47.3%) than three other Pontibacter strains distributed in different species (summarized in Table 1). We found that most of the transposase-related genes are located near genomic islands and next to DNA repair- and ROS detoxification-related genes (Fig. 1), which implied that mobile genetic elements (MGEs) play an important role in the adaptation to radiation and desiccation in the desert. Differences in genome size and MGEs between strain X14-1 and other Pontibacter strains may be attributed to the genomic evolution or gapped assembly of P. actiniarum DSM 19842, Pontibacter sp. BAB1700 and P. roseus DSM 17521. To confirm the phylotype of strain X14-1, we downloaded 40 genomes of the family Cytophagaceae (higher taxonomic classification of the genus Pontibacter) available in the NCBI database. Phylogeny analysis indicated the same results as those previously reported based on the 16S rDNA sequence2, and P. actiniarum DSM 19842 was found to be the most homologous to strain X14-1 (Fig. 2), a finding that is also supported by the following functional analysis. Sugar metabolism in strain X14-1 In comparison with three other genomes from the genus Pontibacter, we found that only strain X14-1 harbors genes encoding D-fructokinase, which is essential for sucrose and fructose utilization, and this is consistent with previous experimental results2. Beta-galactosidase, which is essential for strain X14-1 to use lactose as an alternative carbon source by catalyzing lactose to galactose and glucose, is specific to strain X14-1 compared with other Pontibacter strains. Although comparative analysis revealed a common dispersion of cellobiose glucohydrolase in Pontibacter, enzymes responsible for degrading cellulose to cellobiose are only distributed in strain X14-1. Starch could be degraded to amylose and alpha-D-galactose-1-phosphate, which is an intermediate in the production of UDP-glucose that could link pentose and glucuronate interconversion. This is important for the utilization of D-galactose as a carbon resource. Mannose can enter glycolysis through beta-D-fructose-6-phosphate with the help of hexokinase and mannose-6-phosphate isomerase, which could be encoded by genes in strain X14-1. The ability to utilize versatile sugars as described above (Fig. 3) could partly explain how strain X14-1 survives in an infertile desert. Determinants in the genome for the adaptation of strain X14-1 to the radiation and aridness of the desert Genes related to DNA repair and the stress response, including ROS detoxification and the osmotic response, were analyzed to understand whether and how strain X14-1 protect itself against to desiccation and radiation. In strain X14-1, recombination repair-related genes are the most abundant, followed by base excision repair (BER) and nucleotide excision repair (NER) (Fig. 4, detailed in Supplementary File 1). There are no specific corresponding genetic determinants for the radiation resistance of strain X14-1, implying that the robustness of the resistance of strain X14-1 to radiation and desiccation is not attributed to new genes but to the high-efficacy operation of systems relevant to recovery from radiation. Trehalose is a natural product that can form a protective film outside the cell under low and high temperatures, osmotic pressure and aridness, preventing proteins from becoming inactivated. The analyses conducted in this study elucidated that intermediate products from the metabolism of several sugars, such as starch and lactose, in strain X14-1 could be transferred to trehalose through the OstAB and TreYZ pathways, as summarized in a previous publication35 (Fig. 5). Trehalase-encoding genes, which are homologous to the trehalase-like protein-encoding genes of Gramella forsetii KT0803, were also identified, suggesting the trehalose consumption of strain X 14-1. Responses of strain X14-1 to gamma radiation revealed by RNA-Seq After radiation, the average survival rate of strain X14-1 ranged from 43.4% to 55.1%, as determined from triplicate analyses. To better understand the specific responses of strain X14-1 to gamma radiation, we determined the expression level of all possible determinants responsible for radiation adaptation and ROS detoxification. BER, NER and homologous recombination-related genes, which are mainly required for recovery from gamma radiation, are apparently upregulated (Supplementary File 2). Catalase-encoding genes exhibit intensive upregulation after radiation, but a similar finding was not obtained for superoxide dismutase, which is also important for ROS detoxification (Supplementary File 2). The expression of the manganese transport protein is 10-fold higher after radiation, and a similar trend was identified for subunit A-F of the multicomponent Na+:H+ antiporter as well as other ion-coupled transporters (Supplementary File 2). NADH dehydrogenase I, with the exception of except subunit C, exhibited robust expression in the radiated samples compared with the controls. Ferric uptake regulator and LysR family transcriptional regulator, both of which are transcription factors and have been identified as determinants for oxidative detoxification in the genome, also exhibited increased expression after radiation (Supplementary File 2). Initial evaluation of unannotated genes with differential expression Apparent differential expression was found for 394 of 893 unclassified genes as a result of gamma radiation. To explore whether these unknown genes contribute to the adaptation of strain X14-1 to radiation, an alignment of these genes to the conserved domain database (CDD) was conducted. Of the 394 new genes, 73 could be assigned to specific domains in CDD, and the corresponding hypothesized functions are summarized in Supplementary File 3, indicating that these 73 genes are commonly found in the FecR, DUF and FlgD_ig (associated with ferric transport and preservation) superfamily.
title	Results
sec	Genomic characteristics and phylogeny of strain X14-1 The complete genome sequence of strain X14-1 was produced based on high-quality reads and corrected by read mapping and PCR verification. Strain X14-1 has a larger genome size (5.46 MB) and a lower GC content (47.3%) than three other Pontibacter strains distributed in different species (summarized in Table 1). We found that most of the transposase-related genes are located near genomic islands and next to DNA repair- and ROS detoxification-related genes (Fig. 1), which implied that mobile genetic elements (MGEs) play an important role in the adaptation to radiation and desiccation in the desert. Differences in genome size and MGEs between strain X14-1 and other Pontibacter strains may be attributed to the genomic evolution or gapped assembly of P. actiniarum DSM 19842, Pontibacter sp. BAB1700 and P. roseus DSM 17521. To confirm the phylotype of strain X14-1, we downloaded 40 genomes of the family Cytophagaceae (higher taxonomic classification of the genus Pontibacter) available in the NCBI database. Phylogeny analysis indicated the same results as those previously reported based on the 16S rDNA sequence2, and P. actiniarum DSM 19842 was found to be the most homologous to strain X14-1 (Fig. 2), a finding that is also supported by the following functional analysis.
title	Genomic characteristics and phylogeny of strain X14-1
p	The complete genome sequence of strain X14-1 was produced based on high-quality reads and corrected by read mapping and PCR verification. Strain X14-1 has a larger genome size (5.46 MB) and a lower GC content (47.3%) than three other Pontibacter strains distributed in different species (summarized in Table 1). We found that most of the transposase-related genes are located near genomic islands and next to DNA repair- and ROS detoxification-related genes (Fig. 1), which implied that mobile genetic elements (MGEs) play an important role in the adaptation to radiation and desiccation in the desert. Differences in genome size and MGEs between strain X14-1 and other Pontibacter strains may be attributed to the genomic evolution or gapped assembly of P. actiniarum DSM 19842, Pontibacter sp. BAB1700 and P. roseus DSM 17521.
p	To confirm the phylotype of strain X14-1, we downloaded 40 genomes of the family Cytophagaceae (higher taxonomic classification of the genus Pontibacter) available in the NCBI database. Phylogeny analysis indicated the same results as those previously reported based on the 16S rDNA sequence2, and P. actiniarum DSM 19842 was found to be the most homologous to strain X14-1 (Fig. 2), a finding that is also supported by the following functional analysis.
sec	Sugar metabolism in strain X14-1 In comparison with three other genomes from the genus Pontibacter, we found that only strain X14-1 harbors genes encoding D-fructokinase, which is essential for sucrose and fructose utilization, and this is consistent with previous experimental results2. Beta-galactosidase, which is essential for strain X14-1 to use lactose as an alternative carbon source by catalyzing lactose to galactose and glucose, is specific to strain X14-1 compared with other Pontibacter strains. Although comparative analysis revealed a common dispersion of cellobiose glucohydrolase in Pontibacter, enzymes responsible for degrading cellulose to cellobiose are only distributed in strain X14-1. Starch could be degraded to amylose and alpha-D-galactose-1-phosphate, which is an intermediate in the production of UDP-glucose that could link pentose and glucuronate interconversion. This is important for the utilization of D-galactose as a carbon resource. Mannose can enter glycolysis through beta-D-fructose-6-phosphate with the help of hexokinase and mannose-6-phosphate isomerase, which could be encoded by genes in strain X14-1. The ability to utilize versatile sugars as described above (Fig. 3) could partly explain how strain X14-1 survives in an infertile desert.
title	Sugar metabolism in strain X14-1
p	In comparison with three other genomes from the genus Pontibacter, we found that only strain X14-1 harbors genes encoding D-fructokinase, which is essential for sucrose and fructose utilization, and this is consistent with previous experimental results2. Beta-galactosidase, which is essential for strain X14-1 to use lactose as an alternative carbon source by catalyzing lactose to galactose and glucose, is specific to strain X14-1 compared with other Pontibacter strains. Although comparative analysis revealed a common dispersion of cellobiose glucohydrolase in Pontibacter, enzymes responsible for degrading cellulose to cellobiose are only distributed in strain X14-1. Starch could be degraded to amylose and alpha-D-galactose-1-phosphate, which is an intermediate in the production of UDP-glucose that could link pentose and glucuronate interconversion. This is important for the utilization of D-galactose as a carbon resource. Mannose can enter glycolysis through beta-D-fructose-6-phosphate with the help of hexokinase and mannose-6-phosphate isomerase, which could be encoded by genes in strain X14-1. The ability to utilize versatile sugars as described above (Fig. 3) could partly explain how strain X14-1 survives in an infertile desert.
sec	Determinants in the genome for the adaptation of strain X14-1 to the radiation and aridness of the desert Genes related to DNA repair and the stress response, including ROS detoxification and the osmotic response, were analyzed to understand whether and how strain X14-1 protect itself against to desiccation and radiation. In strain X14-1, recombination repair-related genes are the most abundant, followed by base excision repair (BER) and nucleotide excision repair (NER) (Fig. 4, detailed in Supplementary File 1). There are no specific corresponding genetic determinants for the radiation resistance of strain X14-1, implying that the robustness of the resistance of strain X14-1 to radiation and desiccation is not attributed to new genes but to the high-efficacy operation of systems relevant to recovery from radiation. Trehalose is a natural product that can form a protective film outside the cell under low and high temperatures, osmotic pressure and aridness, preventing proteins from becoming inactivated. The analyses conducted in this study elucidated that intermediate products from the metabolism of several sugars, such as starch and lactose, in strain X14-1 could be transferred to trehalose through the OstAB and TreYZ pathways, as summarized in a previous publication35 (Fig. 5). Trehalase-encoding genes, which are homologous to the trehalase-like protein-encoding genes of Gramella forsetii KT0803, were also identified, suggesting the trehalose consumption of strain X 14-1.
title	Determinants in the genome for the adaptation of strain X14-1 to the radiation and aridness of the desert
p	Genes related to DNA repair and the stress response, including ROS detoxification and the osmotic response, were analyzed to understand whether and how strain X14-1 protect itself against to desiccation and radiation. In strain X14-1, recombination repair-related genes are the most abundant, followed by base excision repair (BER) and nucleotide excision repair (NER) (Fig. 4, detailed in Supplementary File 1). There are no specific corresponding genetic determinants for the radiation resistance of strain X14-1, implying that the robustness of the resistance of strain X14-1 to radiation and desiccation is not attributed to new genes but to the high-efficacy operation of systems relevant to recovery from radiation.
p	Trehalose is a natural product that can form a protective film outside the cell under low and high temperatures, osmotic pressure and aridness, preventing proteins from becoming inactivated. The analyses conducted in this study elucidated that intermediate products from the metabolism of several sugars, such as starch and lactose, in strain X14-1 could be transferred to trehalose through the OstAB and TreYZ pathways, as summarized in a previous publication35 (Fig. 5). Trehalase-encoding genes, which are homologous to the trehalase-like protein-encoding genes of Gramella forsetii KT0803, were also identified, suggesting the trehalose consumption of strain X 14-1.
sec	Responses of strain X14-1 to gamma radiation revealed by RNA-Seq After radiation, the average survival rate of strain X14-1 ranged from 43.4% to 55.1%, as determined from triplicate analyses. To better understand the specific responses of strain X14-1 to gamma radiation, we determined the expression level of all possible determinants responsible for radiation adaptation and ROS detoxification. BER, NER and homologous recombination-related genes, which are mainly required for recovery from gamma radiation, are apparently upregulated (Supplementary File 2). Catalase-encoding genes exhibit intensive upregulation after radiation, but a similar finding was not obtained for superoxide dismutase, which is also important for ROS detoxification (Supplementary File 2). The expression of the manganese transport protein is 10-fold higher after radiation, and a similar trend was identified for subunit A-F of the multicomponent Na+:H+ antiporter as well as other ion-coupled transporters (Supplementary File 2). NADH dehydrogenase I, with the exception of except subunit C, exhibited robust expression in the radiated samples compared with the controls. Ferric uptake regulator and LysR family transcriptional regulator, both of which are transcription factors and have been identified as determinants for oxidative detoxification in the genome, also exhibited increased expression after radiation (Supplementary File 2).
title	Responses of strain X14-1 to gamma radiation revealed by RNA-Seq
p	After radiation, the average survival rate of strain X14-1 ranged from 43.4% to 55.1%, as determined from triplicate analyses. To better understand the specific responses of strain X14-1 to gamma radiation, we determined the expression level of all possible determinants responsible for radiation adaptation and ROS detoxification. BER, NER and homologous recombination-related genes, which are mainly required for recovery from gamma radiation, are apparently upregulated (Supplementary File 2).
p	Catalase-encoding genes exhibit intensive upregulation after radiation, but a similar finding was not obtained for superoxide dismutase, which is also important for ROS detoxification (Supplementary File 2). The expression of the manganese transport protein is 10-fold higher after radiation, and a similar trend was identified for subunit A-F of the multicomponent Na+:H+ antiporter as well as other ion-coupled transporters (Supplementary File 2). NADH dehydrogenase I, with the exception of except subunit C, exhibited robust expression in the radiated samples compared with the controls. Ferric uptake regulator and LysR family transcriptional regulator, both of which are transcription factors and have been identified as determinants for oxidative detoxification in the genome, also exhibited increased expression after radiation (Supplementary File 2).
sec	Initial evaluation of unannotated genes with differential expression Apparent differential expression was found for 394 of 893 unclassified genes as a result of gamma radiation. To explore whether these unknown genes contribute to the adaptation of strain X14-1 to radiation, an alignment of these genes to the conserved domain database (CDD) was conducted. Of the 394 new genes, 73 could be assigned to specific domains in CDD, and the corresponding hypothesized functions are summarized in Supplementary File 3, indicating that these 73 genes are commonly found in the FecR, DUF and FlgD_ig (associated with ferric transport and preservation) superfamily.
title	Initial evaluation of unannotated genes with differential expression
p	Apparent differential expression was found for 394 of 893 unclassified genes as a result of gamma radiation. To explore whether these unknown genes contribute to the adaptation of strain X14-1 to radiation, an alignment of these genes to the conserved domain database (CDD) was conducted. Of the 394 new genes, 73 could be assigned to specific domains in CDD, and the corresponding hypothesized functions are summarized in Supplementary File 3, indicating that these 73 genes are commonly found in the FecR, DUF and FlgD_ig (associated with ferric transport and preservation) superfamily.
sec	Discussion The desert is characterized by a low diversity of inhabitants due to fluctuations in temperature, aridness, strong radiation and poor nutrition136. Elucidating the microbial life in the desert could be important because it could help us understand the dry limit of life and improve our research on extremophiles37. Contributors to radiation resistance and adaptation to poor nutrition in microbes could be applied to improve plant culturing in the desert. This will decrease sand storms mainly induced by the desert and increase arable lands to aid the avoidance of hunger worldwide. Determinants to radiation resistance in microbes residing in the desert could also be industrially applied in cosmetic and soil remediation under radiation. Strain X14-1, which was isolated from the sand surface in the desert, belongs to the poorly understood genus Pontibacter, and there are no publications on an in-depth analysis of the Pontibacter genome. In this study, we performed a high-coverage sequencing of strain X14-1 and produced the first complete map, which may be used as a valuable reference to promote research on the genus Pontibacter. A detailed annotation of the strain X14-1 genome in combination with a comparative transcriptomic analysis was conducted to understand how strain X14-1 utilizes diverse alternative nutrition sources and recovers from desiccation and strong radiation. Previous studies have revealed that strain X14-1 can utilize versatile carbon sources2, and this finding is supported by the reconstruction of the metabolic pathways of different sugars. D-fructokinase, beta-galactosidase, endoglucanase and beta-glucosidase are key for the utilization of sucrose, fructose, lactose, cellulose and glycoside, respectively. Enzymes that are important for metabolizing starch, galactose and mannose were also found in strain X14-1. In addition to their robust competence for utilizing diverse available sugars, another interesting finding is that only strain X14-1 and P. actiniarum DSM19842 possess phosphoenolpyruvate (PEP) carboxylase, which is key for CO2 fixation, mainly in plants, through the C4 cycle and may be used to resemble the autotrophic lifestyle in environments with extremely poor nutrition. The fact that most of the hypothetical proteins in strain X14-1 are homologous to those in P. actiniarum DSM19842 also demonstrates the homology of these two Pontibacter strains revealed by the phylogenetic construction (Fig. 2). In addition to the various metabolic pathways required for the utilization of limited resources, strain X14-1 also harbors an intensive arsenal of DNA repair- and stress response-related determinants that allow it to survive in the desert. An in-depth analysis of the genomic components of strain X14-1 revealed possible broad-spectrum contributors to resistance to radiation (Supplementary File 1). This is consistent with the genetic factors found in the genus Deinococcus, and we also found new genes that encode resistance agents, indicating supplementary mechanisms to those found in Deinococcus. Based on changes in the annotated genes after gamma radiation, we found that catalases but not dismutase exhibited a marked increase, indicating the robustness of catalase in oxidative detoxification in strain X14-1. An increased Mn(II)/Fe ratio has been demonstrated to play an important role in detoxification from oxidative damage11, and this is also supported by our findings: upregulation of the manganese transport protein, ferric uptake regulator and LysR family transcriptional regulator. The expression of gene gshB encoding glutathione synthase was high during radiation, but glutathione peroxidase and reductase, which are responsible for the circulation of GSS (reduced glutathione) and GSSG (oxidized glutathione) exhibited no change, implying that glutathione is not the intermediate in strain X14-1during radiation resistance. An interesting finding is the robustly elevated expression of transposase during gamma radiation (Supplementary File 2), which implies the potently active role of widely distributed transposases in strain X14-1 for desert adaptation. Some of the conserved domains of upregulated unknown genes could be assigned to the FecR and ferritin-like superfamily. This finding may imply that some unknown genes could contribute to oxidative detoxification because these proteins are related with ferric transport or preservation, but further analysis or experiments are required to confirm this hypothesis.
title	Discussion
p	The desert is characterized by a low diversity of inhabitants due to fluctuations in temperature, aridness, strong radiation and poor nutrition136. Elucidating the microbial life in the desert could be important because it could help us understand the dry limit of life and improve our research on extremophiles37. Contributors to radiation resistance and adaptation to poor nutrition in microbes could be applied to improve plant culturing in the desert. This will decrease sand storms mainly induced by the desert and increase arable lands to aid the avoidance of hunger worldwide. Determinants to radiation resistance in microbes residing in the desert could also be industrially applied in cosmetic and soil remediation under radiation.
p	Strain X14-1, which was isolated from the sand surface in the desert, belongs to the poorly understood genus Pontibacter, and there are no publications on an in-depth analysis of the Pontibacter genome. In this study, we performed a high-coverage sequencing of strain X14-1 and produced the first complete map, which may be used as a valuable reference to promote research on the genus Pontibacter. A detailed annotation of the strain X14-1 genome in combination with a comparative transcriptomic analysis was conducted to understand how strain X14-1 utilizes diverse alternative nutrition sources and recovers from desiccation and strong radiation.
p	Previous studies have revealed that strain X14-1 can utilize versatile carbon sources2, and this finding is supported by the reconstruction of the metabolic pathways of different sugars. D-fructokinase, beta-galactosidase, endoglucanase and beta-glucosidase are key for the utilization of sucrose, fructose, lactose, cellulose and glycoside, respectively. Enzymes that are important for metabolizing starch, galactose and mannose were also found in strain X14-1. In addition to their robust competence for utilizing diverse available sugars, another interesting finding is that only strain X14-1 and P. actiniarum DSM19842 possess phosphoenolpyruvate (PEP) carboxylase, which is key for CO2 fixation, mainly in plants, through the C4 cycle and may be used to resemble the autotrophic lifestyle in environments with extremely poor nutrition. The fact that most of the hypothetical proteins in strain X14-1 are homologous to those in P. actiniarum DSM19842 also demonstrates the homology of these two Pontibacter strains revealed by the phylogenetic construction (Fig. 2).
p	In addition to the various metabolic pathways required for the utilization of limited resources, strain X14-1 also harbors an intensive arsenal of DNA repair- and stress response-related determinants that allow it to survive in the desert. An in-depth analysis of the genomic components of strain X14-1 revealed possible broad-spectrum contributors to resistance to radiation (Supplementary File 1). This is consistent with the genetic factors found in the genus Deinococcus, and we also found new genes that encode resistance agents, indicating supplementary mechanisms to those found in Deinococcus. Based on changes in the annotated genes after gamma radiation, we found that catalases but not dismutase exhibited a marked increase, indicating the robustness of catalase in oxidative detoxification in strain X14-1. An increased Mn(II)/Fe ratio has been demonstrated to play an important role in detoxification from oxidative damage11, and this is also supported by our findings: upregulation of the manganese transport protein, ferric uptake regulator and LysR family transcriptional regulator. The expression of gene gshB encoding glutathione synthase was high during radiation, but glutathione peroxidase and reductase, which are responsible for the circulation of GSS (reduced glutathione) and GSSG (oxidized glutathione) exhibited no change, implying that glutathione is not the intermediate in strain X14-1during radiation resistance. An interesting finding is the robustly elevated expression of transposase during gamma radiation (Supplementary File 2), which implies the potently active role of widely distributed transposases in strain X14-1 for desert adaptation. Some of the conserved domains of upregulated unknown genes could be assigned to the FecR and ferritin-like superfamily. This finding may imply that some unknown genes could contribute to oxidative detoxification because these proteins are related with ferric transport or preservation, but further analysis or experiments are required to confirm this hypothesis.
sec	Materials and Methods Strains and culture conditions All of the strains were obtained from the China Center for Type Culture Collection (CCTCC). P. korlensis X14-1 was grown at 30 °C on marine broth 2216 (Difco). To determine the tolerance of the culture to gamma radiation, the strains were grown in the appropriate liquid medium to the exponential phase, and 2.5 ml of the cell culture was then subjected to 60Co radiation with a continuous dose rate of 1000 Gray/h at room temperature. The control culture cells were not treated with radiation. Samples were collected 3 h after culture with (radiated samples) or without (controls) radiation treatment. Biological triplicates were harvested from each treatment. These samples were immediately frozen in RNAlater (Qiagen, German) and stored at –80 °C for RNA analysis. DNA extraction, sequencing and assembly Total DNA was extracted through the traditional CTAB method, and its quality was tested by Qubit. Sequencing was performed for a 470-bp insert size library (1.65-Gb clean reads with a length of 90 bp) and mate-paired libraries with insert sizes of 2222 bp and 6200 bp (477-Mb clean reads, each with a length of 49 bp), respectively. Clean reads were used to assemble the complete genome sequence as previously described38. Gene prediction, annotation and comparative analysis Genes were predicted based on the assembled and confirmed sequence using GLIMMER394041, which was developed for microorganisms including bacteria, archaea, and viruses. This was also used for the prediction of genes from other three Pontibacter strains downloaded from the NCBI ftp site. The annotations were conducted by assigning the predicted genes to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/)42, Cluster of Orthologous Groups of proteins (COG)4344, and the protein database in NCBI (http://www.ncbi.nlm.nih.gov/protein/) through a BLASTP search (e-value ≤1e-5, query coverage ≥40%). The strategy described by Tian et al.45 was used to identify the orthologous proteins in four Pontibacter genomes. In the phylogenetic analysis, homologous genes from 41 strains were used to construct the gene family and phylogenetic tree: (i) the gene set was first generated using BLAST (V2.2.23) with the genes from 41 strains; (ii) the gene family was constructed using OrthoMCL (V1.4) with the gene set data; (iii) to construct the phylogenetic tree, single-copy genes in the gene family were screened, and the genes were aligned according to the strains using Muscle (V3.8.31); and (IV) the phylogenetic tree was drawn using TreeBeST (V1.9.2) with 1000 bootstrap replicates. Identification of MGEs The prophages in the genome sequences were identified using the Phage Search Tool (PHAST)46. Genomic islands were predicted with multiple methods (IslandPath-DIOMB, SIGI-HMM, and IslandPicker using Islandviewer474849. We searched repeat sequences with the RepeatMasker program50, and transposases were selected from the annotation results. RNA-Seq analysis Three replicates of each group were used for total RNA extraction and mRNA enrichment using the Ribo-ZeroTM kit. The qualified mRNA was sequenced on the HiSeq2000 platform, and the expression level of each gene was calculated using the reads per kilo bases per million reads (RPKM)51. The genes that were differentially expressed between the two groups were identified as described by Audicet et al.52. The analysis results with a p-value <0.05 were corrected by the FDR (false discovery rate), which was set to ≤0.001.
title	Materials and Methods
sec	Strains and culture conditions All of the strains were obtained from the China Center for Type Culture Collection (CCTCC). P. korlensis X14-1 was grown at 30 °C on marine broth 2216 (Difco). To determine the tolerance of the culture to gamma radiation, the strains were grown in the appropriate liquid medium to the exponential phase, and 2.5 ml of the cell culture was then subjected to 60Co radiation with a continuous dose rate of 1000 Gray/h at room temperature. The control culture cells were not treated with radiation. Samples were collected 3 h after culture with (radiated samples) or without (controls) radiation treatment. Biological triplicates were harvested from each treatment. These samples were immediately frozen in RNAlater (Qiagen, German) and stored at –80 °C for RNA analysis.
title	Strains and culture conditions
p	All of the strains were obtained from the China Center for Type Culture Collection (CCTCC). P. korlensis X14-1 was grown at 30 °C on marine broth 2216 (Difco). To determine the tolerance of the culture to gamma radiation, the strains were grown in the appropriate liquid medium to the exponential phase, and 2.5 ml of the cell culture was then subjected to 60Co radiation with a continuous dose rate of 1000 Gray/h at room temperature. The control culture cells were not treated with radiation. Samples were collected 3 h after culture with (radiated samples) or without (controls) radiation treatment. Biological triplicates were harvested from each treatment. These samples were immediately frozen in RNAlater (Qiagen, German) and stored at –80 °C for RNA analysis.
sec	DNA extraction, sequencing and assembly Total DNA was extracted through the traditional CTAB method, and its quality was tested by Qubit. Sequencing was performed for a 470-bp insert size library (1.65-Gb clean reads with a length of 90 bp) and mate-paired libraries with insert sizes of 2222 bp and 6200 bp (477-Mb clean reads, each with a length of 49 bp), respectively. Clean reads were used to assemble the complete genome sequence as previously described38.
title	DNA extraction, sequencing and assembly
p	Total DNA was extracted through the traditional CTAB method, and its quality was tested by Qubit. Sequencing was performed for a 470-bp insert size library (1.65-Gb clean reads with a length of 90 bp) and mate-paired libraries with insert sizes of 2222 bp and 6200 bp (477-Mb clean reads, each with a length of 49 bp), respectively. Clean reads were used to assemble the complete genome sequence as previously described38.
sec	Gene prediction, annotation and comparative analysis Genes were predicted based on the assembled and confirmed sequence using GLIMMER394041, which was developed for microorganisms including bacteria, archaea, and viruses. This was also used for the prediction of genes from other three Pontibacter strains downloaded from the NCBI ftp site. The annotations were conducted by assigning the predicted genes to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/)42, Cluster of Orthologous Groups of proteins (COG)4344, and the protein database in NCBI (http://www.ncbi.nlm.nih.gov/protein/) through a BLASTP search (e-value ≤1e-5, query coverage ≥40%). The strategy described by Tian et al.45 was used to identify the orthologous proteins in four Pontibacter genomes. In the phylogenetic analysis, homologous genes from 41 strains were used to construct the gene family and phylogenetic tree: (i) the gene set was first generated using BLAST (V2.2.23) with the genes from 41 strains; (ii) the gene family was constructed using OrthoMCL (V1.4) with the gene set data; (iii) to construct the phylogenetic tree, single-copy genes in the gene family were screened, and the genes were aligned according to the strains using Muscle (V3.8.31); and (IV) the phylogenetic tree was drawn using TreeBeST (V1.9.2) with 1000 bootstrap replicates.
title	Gene prediction, annotation and comparative analysis
p	Genes were predicted based on the assembled and confirmed sequence using GLIMMER394041, which was developed for microorganisms including bacteria, archaea, and viruses. This was also used for the prediction of genes from other three Pontibacter strains downloaded from the NCBI ftp site. The annotations were conducted by assigning the predicted genes to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/)42, Cluster of Orthologous Groups of proteins (COG)4344, and the protein database in NCBI (http://www.ncbi.nlm.nih.gov/protein/) through a BLASTP search (e-value ≤1e-5, query coverage ≥40%). The strategy described by Tian et al.45 was used to identify the orthologous proteins in four Pontibacter genomes.
p	In the phylogenetic analysis, homologous genes from 41 strains were used to construct the gene family and phylogenetic tree: (i) the gene set was first generated using BLAST (V2.2.23) with the genes from 41 strains; (ii) the gene family was constructed using OrthoMCL (V1.4) with the gene set data; (iii) to construct the phylogenetic tree, single-copy genes in the gene family were screened, and the genes were aligned according to the strains using Muscle (V3.8.31); and (IV) the phylogenetic tree was drawn using TreeBeST (V1.9.2) with 1000 bootstrap replicates.
sec	Identification of MGEs The prophages in the genome sequences were identified using the Phage Search Tool (PHAST)46. Genomic islands were predicted with multiple methods (IslandPath-DIOMB, SIGI-HMM, and IslandPicker using Islandviewer474849. We searched repeat sequences with the RepeatMasker program50, and transposases were selected from the annotation results.
title	Identification of MGEs
p	The prophages in the genome sequences were identified using the Phage Search Tool (PHAST)46. Genomic islands were predicted with multiple methods (IslandPath-DIOMB, SIGI-HMM, and IslandPicker using Islandviewer474849. We searched repeat sequences with the RepeatMasker program50, and transposases were selected from the annotation results.
sec	RNA-Seq analysis Three replicates of each group were used for total RNA extraction and mRNA enrichment using the Ribo-ZeroTM kit. The qualified mRNA was sequenced on the HiSeq2000 platform, and the expression level of each gene was calculated using the reads per kilo bases per million reads (RPKM)51. The genes that were differentially expressed between the two groups were identified as described by Audicet et al.52. The analysis results with a p-value <0.05 were corrected by the FDR (false discovery rate), which was set to ≤0.001.
title	RNA-Seq analysis
p	Three replicates of each group were used for total RNA extraction and mRNA enrichment using the Ribo-ZeroTM kit. The qualified mRNA was sequenced on the HiSeq2000 platform, and the expression level of each gene was calculated using the reads per kilo bases per million reads (RPKM)51. The genes that were differentially expressed between the two groups were identified as described by Audicet et al.52. The analysis results with a p-value <0.05 were corrected by the FDR (false discovery rate), which was set to ≤0.001.
sec	Additional Information How to cite this article: Dai, J. et al. Unraveling adaptation of Pontibacter korlensis to radiation and infertility in desert through complete genome and comparative transcriptomic analysis. Sci. Rep. 5, 10929; doi: 10.1038/srep10929 (2015).
title	Additional Information
p	How to cite this article: Dai, J. et al. Unraveling adaptation of Pontibacter korlensis to radiation and infertility in desert through complete genome and comparative transcriptomic analysis. Sci. Rep. 5, 10929; doi: 10.1038/srep10929 (2015).
sec	Supplementary Material Supplementary File 1 Supplementary File 2 Supplementary File 3
title	Supplementary Material
caption	Supplementary File 1
title	Supplementary File 1
caption	Supplementary File 2
title	Supplementary File 2
caption	Supplementary File 3
title	Supplementary File 3
back	We would like to thank Xiaohui Li, Hongjuan Wang and Yongfeng Liu for coordination during this project. This research was supported by the National Natural Science Foundation of China (No: 31300003), the Hubei Provincial Natural Science Foundation (2015CFB679) and Research Fund for the Doctoral Program of Hubei University of Technology (No: BSQD13003). Author Contributions X.C. managed the project. C.F. and Z.W. offered the strain and designed the experiments. Y.Z., L.Z., Q.Y., X.L., Z.W., M.Z. and Z.J. performed experiments and prepared the DNA and RNA samples. Z.Y. and Q.G. performed sequencing and assembly. C.Q. and W.D. performed other bioinformatics analysis in this work. W.D., C.Q. and J.D. wrote the paper. All authors reviewed this manuscript. Accession codes: Assembly result has been deposited in GenBank nucleotide core database under the accession number CP009621. Clean reads for six transcriptome samples could be found through accession number SRP049974.
ack	We would like to thank Xiaohui Li, Hongjuan Wang and Yongfeng Liu for coordination during this project. This research was supported by the National Natural Science Foundation of China (No: 31300003), the Hubei Provincial Natural Science Foundation (2015CFB679) and Research Fund for the Doctoral Program of Hubei University of Technology (No: BSQD13003).
p	We would like to thank Xiaohui Li, Hongjuan Wang and Yongfeng Liu for coordination during this project. This research was supported by the National Natural Science Foundation of China (No: 31300003), the Hubei Provincial Natural Science Foundation (2015CFB679) and Research Fund for the Doctoral Program of Hubei University of Technology (No: BSQD13003).
footnote	Author Contributions X.C. managed the project. C.F. and Z.W. offered the strain and designed the experiments. Y.Z., L.Z., Q.Y., X.L., Z.W., M.Z. and Z.J. performed experiments and prepared the DNA and RNA samples. Z.Y. and Q.G. performed sequencing and assembly. C.Q. and W.D. performed other bioinformatics analysis in this work. W.D., C.Q. and J.D. wrote the paper. All authors reviewed this manuscript. Accession codes: Assembly result has been deposited in GenBank nucleotide core database under the accession number CP009621. Clean reads for six transcriptome samples could be found through accession number SRP049974.
p	Author Contributions X.C. managed the project. C.F. and Z.W. offered the strain and designed the experiments. Y.Z., L.Z., Q.Y., X.L., Z.W., M.Z. and Z.J. performed experiments and prepared the DNA and RNA samples. Z.Y. and Q.G. performed sequencing and assembly. C.Q. and W.D. performed other bioinformatics analysis in this work. W.D., C.Q. and J.D. wrote the paper. All authors reviewed this manuscript.
p	Accession codes: Assembly result has been deposited in GenBank nucleotide core database under the accession number CP009621. Clean reads for six transcriptome samples could be found through accession number SRP049974.
figure	Figure 1 Genomic components of P. korlensis X14-1. Distribution of genes, genomic islands, transposases, DNA repair and oxidative response related determinants in the genome of strain X14-1. It indicates that strain X14-1 harbors abundant MGEs, suggestive of high genome plasticity.
label	Figure 1
caption	Genomic components of P. korlensis X14-1. Distribution of genes, genomic islands, transposases, DNA repair and oxidative response related determinants in the genome of strain X14-1. It indicates that strain X14-1 harbors abundant MGEs, suggestive of high genome plasticity.
title	Genomic components of P. korlensis X14-1.
p	Distribution of genes, genomic islands, transposases, DNA repair and oxidative response related determinants in the genome of strain X14-1. It indicates that strain X14-1 harbors abundant MGEs, suggestive of high genome plasticity.
figure	Figure 2 Phylogenetic tree of family Cytophagaceae based on 40 available reference genomes and strain X14-1(highlighted in red). The phylogenetic tree was produced on single copy gene within gene family and drawn by TreeBeST version1.9.2 under bootstrap 1000, implicating taxonomic position of strain X14-1.
label	Figure 2
caption	Phylogenetic tree of family Cytophagaceae based on 40 available reference genomes and strain X14-1(highlighted in red). The phylogenetic tree was produced on single copy gene within gene family and drawn by TreeBeST version1.9.2 under bootstrap 1000, implicating taxonomic position of strain X14-1.
title	Phylogenetic tree of family Cytophagaceae based on 40 available reference genomes and strain X14-1(highlighted in red).
p	The phylogenetic tree was produced on single copy gene within gene family and drawn by TreeBeST version1.9.2 under bootstrap 1000, implicating taxonomic position of strain X14-1.
figure	Figure 3 Pathways of sugar utilization for strain X14-1. PGI: Pentose and Glucuronate Interconversion; PPP: Pentose Phosphate Pathway. Arrows with different colors represent utilizing distinct carbon sources: Blue-lactose and starch, orange-glycoside, gray-sucrose, brown-cellulose, black-glucose, red-galactose, green-mannose, purple-fructose.
label	Figure 3
caption	Pathways of sugar utilization for strain X14-1. PGI: Pentose and Glucuronate Interconversion; PPP: Pentose Phosphate Pathway. Arrows with different colors represent utilizing distinct carbon sources: Blue-lactose and starch, orange-glycoside, gray-sucrose, brown-cellulose, black-glucose, red-galactose, green-mannose, purple-fructose.
title	Pathways of sugar utilization for strain X14-1.
p	PGI: Pentose and Glucuronate Interconversion; PPP: Pentose Phosphate Pathway. Arrows with different colors represent utilizing distinct carbon sources: Blue-lactose and starch, orange-glycoside, gray-sucrose, brown-cellulose, black-glucose, red-galactose, green-mannose, purple-fructose.
figure	Figure 4 Genes related to DNA repair and stress response in genus Pontibacter. BER: base excision repair; NER: nucleotide excision repair; DR: direct reversal; RER: recombination repair; ODR: other DNA repair; ORT: other radiation tolerance.
label	Figure 4
caption	Genes related to DNA repair and stress response in genus Pontibacter. BER: base excision repair; NER: nucleotide excision repair; DR: direct reversal; RER: recombination repair; ODR: other DNA repair; ORT: other radiation tolerance.
title	Genes related to DNA repair and stress response in genus Pontibacter.
p	BER: base excision repair; NER: nucleotide excision repair; DR: direct reversal; RER: recombination repair; ODR: other DNA repair; ORT: other radiation tolerance.
figure	Figure 5 Trehalose synthesis of strain X14-1. UDP-Glucose from lactose and galactose (Fig. 3) could be catalyzed to trehalose through trehalose 6-phosphate synthase and trehalose 6-phosphate phosphatase. Glycogen operon protein could degrade starch synthesized from alpha-D-Glucose-6P available from many metabolic pathways to trehalose.
label	Figure 5
caption	Trehalose synthesis of strain X14-1. UDP-Glucose from lactose and galactose (Fig. 3) could be catalyzed to trehalose through trehalose 6-phosphate synthase and trehalose 6-phosphate phosphatase. Glycogen operon protein could degrade starch synthesized from alpha-D-Glucose-6P available from many metabolic pathways to trehalose.
title	Trehalose synthesis of strain X14-1.
p	UDP-Glucose from lactose and galactose (Fig. 3) could be catalyzed to trehalose through trehalose 6-phosphate synthase and trehalose 6-phosphate phosphatase. Glycogen operon protein could degrade starch synthesized from alpha-D-Glucose-6P available from many metabolic pathways to trehalose.
table-wrap	Table 1 Comparison of strain X14-1 with other species in genus Pontibacter. P. korlensis X14-1 P. actiniarum DSM19842 Pontibacter sp. BAB1700 P. roseus DSM17521 Habitat Desert Aquatic Multiple NA Genome size(Mb) 5.46 4.95 4.54 4.58 GC content(%) 47.3 53.1 50.0 52.6 Total genes 5,037 4,689 4,849 4,260 Coding regions(%) 86.58 86.42 84.97 87.01 Unannotated genes 893 901 507 328 Insertion sequence 4 6 1 0 Prophage 1 0 0 0 Transposase/Integrase 114 28 10 5
label	Table 1
caption	Comparison of strain X14-1 with other species in genus Pontibacter.
title	Comparison of strain X14-1 with other species in genus Pontibacter.
table	P. korlensis X14-1 P. actiniarum DSM19842 Pontibacter sp. BAB1700 P. roseus DSM17521 Habitat Desert Aquatic Multiple NA Genome size(Mb) 5.46 4.95 4.54 4.58 GC content(%) 47.3 53.1 50.0 52.6 Total genes 5,037 4,689 4,849 4,260 Coding regions(%) 86.58 86.42 84.97 87.01 Unannotated genes 893 901 507 328 Insertion sequence 4 6 1 0 Prophage 1 0 0 0 Transposase/Integrase 114 28 10 5
tr	P. korlensis X14-1 P. actiniarum DSM19842 Pontibacter sp. BAB1700 P. roseus DSM17521
th
th	P. korlensis X14-1
th	P. actiniarum DSM19842
th	Pontibacter sp. BAB1700
th	P. roseus DSM17521
tr	Habitat Desert Aquatic Multiple NA
td	Habitat
td	Desert
td	Aquatic
td	Multiple
td	NA
tr	Genome size(Mb) 5.46 4.95 4.54 4.58
td	Genome size(Mb)
td	5.46
td	4.95
td	4.54
td	4.58
tr	GC content(%) 47.3 53.1 50.0 52.6
td	GC content(%)
td	47.3
td	53.1
td	50.0
td	52.6
tr	Total genes 5,037 4,689 4,849 4,260
td	Total genes
td	5,037
td	4,689
td	4,849
td	4,260
tr	Coding regions(%) 86.58 86.42 84.97 87.01
td	Coding regions(%)
td	86.58
td	86.42
td	84.97
td	87.01
tr	Unannotated genes 893 901 507 328
td	Unannotated genes
td	893
td	901
td	507
td	328
tr	Insertion sequence 4 6 1 0
td	Insertion sequence
td	4
td	6
td	1
td	0
tr	Prophage 1 0 0 0
td	Prophage
td	1
td	0
td	0
td	0
tr	Transposase/Integrase 114 28 10 5
td	Transposase/Integrase
td	114
td	28
td	10
td	5

Annnotations TAB TSV DIC JSON TextAE Lectin_function IAV-Glycan

Denotations: 0
Blocks: 0
Relations: 0

PMC:4460873 / 21032-21151 JSONTXT

Document structure show

Annnotations TAB TSV DIC JSON TextAE Lectin_function IAV-Glycan

PMC:4460873 / 21032-21151 JSON TXT