ACR: metagenome-assembled prokaryotic and eukaryotic genome refinement tool.

Microbial genome recovery from metagenomes can further explain microbial ecosystem structures, functions and dynamics. Thus, this study developed the Additional Clustering Refiner (ACR) to enhance high-purity prokaryotic and eukaryotic metagenome-assembled genome (MAGs) recovery. ACR refines low-quality MAGs by subjecting them to iterative k-means clustering predicated on contig abundance and increasing bin purity through validated universal marker genes. Synthetic and real-world metagenomic datasets, including short- and long-read sequences, evaluated ACR's effectiveness. The results demonstrated improved MAG purity and a significant increase in high- and medium-quality MAG recovery rates. In addition, ACR seamlessly integrates with various binning algorithms, augmenting their strengths without modifying core features. Furthermore, its multiple sequencing technology compatibilities expand its applicability. By efficiently recovering high-quality prokaryotic and eukaryotic genomes, ACR is a promising tool for deepening our understanding of microbial communities through genome-centric metagenomics.

A case study on the distribution of the environmental resistome in Korean shrimp farms.

Hundreds of tons of antibiotics are widely used in aquaculture to prevent microbial infections and promote fish growth. However, the overuse of antibiotics and chemical products can lead to the selection and spreading of antibiotic-resistant bacteria (ARB) and antimicrobial resistance genes (ARGs), which are of great concern considering the threat to public health worldwide. Here, in-depth metagenome sequencing was performed to explore the environmental resistome and ARB distribution across farming stages in shrimp farms and examine anthropogenic effects in nearby coastal waters. A genome-centric analysis using a metagenome binning approach allowed us to accurately investigate the distribution of pathogens and ARG hosts in shrimp farms. The diversity of resistomes was higher in shrimp farms than in coastal waters, and the distribution of resistomes was dependent on the farming stage. In particular, the tetracycline resistance gene was found mainly at the early post-larval stage regardless of the farm. The metagenome-assembled genomes of Vibrio spp. were dominant at this stage and harbored tet34, which is known to confer resistance to oxytetracycline. In addition, opportunistic pathogens such as Francisella, Mycoplasma, Photobacterium, and Vibrio were found in abundance in shrimp farms, which had multiple virulence factors. This study highlights the increased resistance diversity and environmental selection of pathogens in shrimp farms. The use of environmental pollutants on farms may cause an increase in resistome diversity/abundance and the transmission of pathogens to the surrounding environment, which may pose future risks to public health and aquatic organisms.

Genomic Information on Linezolid-Resistant Sequence-Type 398 Livestock-Associated Methicillin-Resistant Staphylococcus aureus Isolated from a Pig.

The frequent occurrence of sequence-type 398 (ST398) livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) in pigs has become a major public health concern owing to the increased zoonotic potential of the pathogen. Recently, a novel oxazolidinone resistance gene, chloramphenicol-florfenicol resistant (cfr), conferring multiresistance phenotypes to phenicols, lincosamides, oxazolidinones, pleuromutilins, and streptogramin A (PhLOPSA), has been found among ST398 LA-MRSA strains isolated from pigs. In this study, we report the first in silico genome analysis of a linezolid-resistant ST398 LA-MRSA strain, designated PJFA-521M, recovered from a pig in Korea. Genomic analyses revealed that the presence of the cfr gene was responsible for the observed linezolid resistance in the PJFA-521M strain. Moreover, newer antimicrobial resistance genes, such as the dfrG, aadE, spw, lsa(E), lnu(B), and fexA genes, were found in the PJFA-521M strain. In addition to the genetic elements for antimicrobial resistance, the carriage of various virulence genes for adherence, invasion, and immunomodulation was identified in the genome, especially within several mobile genetic elements (MGEs). The presence of multiple antimicrobial resistance genes and virulence genes on MGEs in the genome of a linezolid-resistant ST398 LA-MRSA should raise awareness regarding the use of other antimicrobial agents in pig farms and may also provide selective pressure for the prevalence of the cfr gene and the associated multidrug-resistant phenotype.

Metagenomic analysis of relationships between the denitrification process and carbon metabolism in a bioaugmented full-scale tannery wastewater treatment plant.

The goal of this study was to investigate the relationship between the denitrification process and carbon metabolism in a full-scale tannery wastewater treatment plant bioaugmented with the microbial consortium BM-S-1. The metagenomic analysis of the microbial community showed that Brachymonas denitrificans, a known denitrifier, was present at a high level in the treatment stages of buffering (B), primary aeration (PA), and sludge digestion (SD). The occurrences of the amino acid-degrading enzymes alpha ketoglutarate dehydrogenase (α-KGDH) and tryptophan synthase were highly correlated with the presence of denitrification genes, such as napA, narG, nosZ and norB. The occurrence of glutamate dehydrogenase (GDH) was also highly paralleled with the occurrence of denitrification genes such as napA, narG, and norZ. The denitrification genes (nosZ, narG, napA, norB and nrfA) and amino acid degradation enzymes (tryptophan synthase, α-KGDH and pyridoxal phosphate dependent enzymes) were observed at high levels in B. This indicates that degradation of amino acids and denitrification of nitrate may potentially occur in B. The high concentrations of the fatty acid degradation enzyme groups (enoyl-CoA hydratase, 3-hydroxyacyl-CoA dehydrogenase and ß-ketothiolase) were observed together with the denitrification genes, such as napA, narG and nosZ. Phospholipase/carboxylesterase, enoyl-CoA hydratase/isomerase, acyl-CoA dehydrogenase, phenylacetate degradation enzyme and 3-hydroxyacyl-CoA dehydrogenase 2 were also dominant in B. All these results clearly indicate that the denitrification pathways are potentially linked to the degradation pathways of amino acids and fatty acids whose degradation products go through the TCA cycle, generating the NADH that is used as electron donors for denitrification.

Gut microbiome and metabolic pathways linked to sleep quality.

Sleep quality is a vital determinant of human health as sleep disorders are associated with cognitive deficits, and chronic sleep deprivation is associated with a broad range of health complications. Previous studies on the association between the gut microbiome and sleep quality have been constrained by small sample sizes or have focused on specific sleep disorders, thus yielding inconsistent results. Herein, we investigated the relationship between microbial composition and sleep quality in a cohort of 159 Koreans. Sleep quality was measured using the Pittsburgh Sleep Quality Index (PSQI), determined through a self-administered questionnaire. Gut microbiome analyses were performed using 16S rRNA amplicons. We found no direct correlation between microbial alpha diversity metrics and sleep; however, we identified differences in beta diversity among sleep quality groups (with a PSQI score > 5 indicating poor sleep quality and PSQI ≤5 indicating good sleep quality). We also found differential microbial signatures (Bacteroides, Prevotella 9, and Faecalibacterium) among the groups. Furthermore, functional metabolic pathway profiles revealed significant linear correlations of the L-arginine and L-tryptophan biosynthetic pathways as well as 4-aminobutanoate degradation with sleep status. In particular, Faecalibacterium prausnitzii, which harbors these metabolic pathways, showed differences between sleep quality groups and a linear association with sleep quality scores and was thus identified as the species most strongly associated with sleep status. This study provides a significant advance in our understanding of the relationship between gut microbiota and sleep regulation. The current findings provide a basis for further research into potential therapeutic strategies for sleep disorders targeting the gut microbiome.

Recovery of 240 metagenome-assembled genomes from coastal mariculture environments in South Korea.

The mariculture industry has seen a rapid expansion in recent years due to the increasing global demand for seafood. However, the industry faces challenges from climate change and increased pathogen pressure. Additionally, the chemicals used to enhance mariculture productivity are changing ocean ecosystems. This study analyzed 36 surface-water metagenomes from South Korean mussel, oyster, scallop, and shrimp farms to expand our understanding of aquaculture microbial genetic resources and the potential impacts of these anthropogenic inputs. We recovered 240 non-redundant species-level metagenome-assembled genomes (MAGs), comprising 224 bacteria, 13 archaea, and three eukaryotes. Most MAGs were assigned to Proteobacteria, Bacteroidota, and Actinobacteriota, with 40.7% remaining unclassified at the species level. Among the three eukaryotic MAGs, one was identified as a novel lineage of green algae, highlighting the uncharacterized genetic diversity in mariculture environments. Additionally, 22 prokaryotic MAGs harbored 26 antibiotic and metal resistance genes, with MAGs carrying beta-lactamases being particularly prevalent in most farms. The obtained microbiome data from mariculture environments can be utilized in future studies to foster healthy, sustainable mariculture practices.

Persistence of antibiotic resistance from animal agricultural effluents to surface water revealed by genome-centric metagenomics.

Concerns about antibiotic resistance genes (ARGs) released from wastewaters of livestock or fish farming into the natural environment are increasing, but studies on unculturable bacteria related to the dissemination of antibiotic resistance are limited. Here, we reconstructed 1100 metagenome-assembled genomes (MAGs) to assess the impact of microbial antibiotic resistome and mobilome in wastewaters discharged to Korean rivers. Our results indicate that ARGs harbored in the MAGs were disseminated from wastewater effluents into downstream rivers. Moreover, it was found that ARGs are more commonly co-localized with mobile genetic elements (MGEs) in agricultural wastewater than in river water. Among the effluent-derived phyla, uncultured members of the superphylum Patescibacteria possessed a high number of MGEs, along with co-localized ARGs. Our findings suggest that members of the Patesibacteria are a potential vector for propagating ARGs into the environmental community. Therefore, we propose that the dissemination of ARGs by uncultured bacteria should be further investigated in multiple environments.

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics.

BACKGROUND: DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life. RESULTS: We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark. CONCLUSIONS: Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome. Video Abstract.

Transcriptomic Analysis of Degradative Pathways for Azo Dye Acid Blue 113 in Sphingomonas melonis B-2 from the Dye Wastewater Treatment Process.

BACKGROUND: Acid Blue 113 (AB113) is a typical azo dye, and the resulting wastewater is toxic and difficult to remove. METHODS: The experimental culture was set up for the biodegradation of the azo dye AB113, and the cell growth and dye decolorization were monitored. Transcriptome sequencing was performed in the presence and absence of AB113 treatment. The key pathways and enzymes involved in AB113 degradation were found through pathway analysis and enrichment software (GO, EggNog and KEGG). RESULTS: S. melonis B-2 achieved more than 80% decolorization within 24 h (50 and 100 mg/L dye). There was a positive relationship between cell growth and the azo dye degradation rate. The expression level of enzymes involved in benzoate and naphthalene degradation pathways (NADH quinone oxidoreductase, N-acetyltransferase and aromatic ring-hydroxylating dioxygenase) increased significantly after the treatment of AB113. CONCLUSIONS: Benzoate and naphthalene degradation pathways were the key pathways for AB113 degradation. NADH quinone oxidoreductase, N-acetyltransferase, aromatic ring-hydroxylating dioxygenase and CYP450 were the key enzymes for AB113 degradation. This study provides evidence for the process of AB113 biodegradation at the molecular and biochemical level that will be useful in monitoring the dye wastewater treatment process at the full-scale treatment.

Distinct Bacterial and Fungal Communities Colonizing Waste Plastic Films Buried for More Than 20 Years in Four Landfill Sites in Korea.

Plastic pollution has been recognized as a serious environmental problem, and microbial degradation of plastics is a potential, environmentally friendly solution to this. Here, we analyzed and compared microbial communities on waste plastic films (WPFs) buried for long periods at four landfill sites with those in nearby soils to identify microbes with the potential to degrade plastics. Fourier-transform infrared spectroscopy spectra of these WPFs showed that most were polyethylene and had signs of oxidation, such as carbon-carbon double bonds, carbon-oxygen single bonds, or hydrogen-oxygen single bonds, but the presence of carbonyl groups was rare. The species richness and diversity of the bacterial and fungal communities on the films were generally lower than those in nearby soils. Principal coordinate analysis of the bacterial and fungal communities showed that their overall structures were determined by their geographical locations; however, the microbial communities on the films were generally different from those in the soils. For the pulled data from the four landfill sites, the relative abundances of Bradyrhizobiaceae, Pseudarthrobacter, Myxococcales, Sphingomonas, and Spartobacteria were higher on films than in soils at the bacterial genus level. At the species level, operational taxonomic units classified as Bradyrhizobiaceae and Pseudarthrobacter in bacteria and Mortierella in fungi were enriched on the films. PICRUSt analysis showed that the predicted functions related to amino acid and carbohydrate metabolism and xenobiotic degradation were more abundant on films than in soils. These results suggest that specific microbial groups were enriched on the WPFs and may be involved in plastic degradation.

Prokaryotic DNA methylation and its functional roles.

DNA methylation is known as a universal mechanism of epigenetic regulation in all kingdoms of life. Particularly, given that prokaryotes lack key elements such as histones and nucleosomes that can structurally modify DNA, DNA methylation is considered a major epigenetic regulator in these organisms. However, because DNA methylation studies have focused primarily on eukaryotes, the mechanism of prokaryotic DNA methylation has been less studied than in eukaryotes. DNA methylation in prokaryotes plays an important role in regulating not only the host defense system, but also the cell cycle, gene expression, and virulence that can respond directly to the environment. Recent advances in sequencing techniques capable of detecting methylation signals have allowed for the characterization of prokaryotic genome-wide epigenetic regulation. In this review, we describe representative examples of cellular events regulated by DNA methylation in prokaryotes, from early studies to current applications.

Two DNA Methyltransferases for Site-Specific 6mA and 5mC DNA Modification in Xanthomonas euvesicatoria.

Xanthomonas euvesicatoria (Xe) is a gram-negative phytopathogenic bacterium that causes bacterial spot disease in tomato/pepper leading to economic losses in plantations. DNA methyltransferases (MTases) are critical for the survival of prokaryotes; however, their functions in phytopathogenic bacteria remain unclear. In this study, we characterized the functions of two putative DNA MTases, XvDMT1 and XvDMT2, in Xe by generating XvDMT1- and XvDMT2-overexpressing strains, Xe(XvDMT1) and Xe(XvDMT2), respectively. Virulence of Xe(XvDMT2), but not Xe(XvDMT1), on tomato was dramatically reduced. To postulate the biological processes involving XvDMTs, we performed a label-free shotgun comparative proteomic analysis, and results suggest that XvDMT1 and XvDMT2 have distinct roles in Xe. We further characterized the functions of XvDMTs using diverse phenotypic assays. Notably, both Xe(XvDMT1) and Xe(XvDMT2) showed growth retardation in the presence of sucrose and fructose as the sole carbon source, with Xe(XvDMT2) being the most severely affected. In addition, biofilm formation and production of exopolysaccharides were declined in Xe(XvDMT2), but not Xe(XvDMT1). Xe(XvDMT2) was more tolerant to EtOH than Xe(XvDMT1), which had enhanced tolerance to sorbitol but decreased tolerance to polymyxin B. Using single-molecule real-time sequencing and methylation-sensitive restriction enzymes, we successfully predicted putative motifs methylated by XvDMT1 and XvDMT2, which are previously uncharacterized 6mA and 5mC DNA MTases, respectively. This study provided new insights into the biological functions of DNA MTases in prokaryotic organisms.

Metagenomic Insight of a Full Scale Eco-Friendly Treatment System of Textile Dye Wastewater Using Bioaugmentation of the Composite Culture CES-1.

Effects of bioaugmentation of the composite microbial culture CES-1 on a full scale textile dye wastewater treatment process were investigated in terms of water quality, sludge reduction, dynamics of microbial community structures and their functional genes responsible for degradation of azo dye, and other chemicals. The removal efficiencies for Chemical Oxygen Demand (COD), Total Nitrogen (T-N), Total Phosphorus (T-P), Suspended Solids (SS), and color intensity (96.4%, 78.4, 83.1, 84.4, and 92.0, respectively) 300-531 days after the augmentation were generally improved after bioaugmentation. The denitrification linked to T-N removal appeared to contribute to the concomitant COD removal that triggered a reduction of sludge (up to 22%) in the same period of augmentation. Azo dye and aromatic compound degradation and other downstream pathways were highly metabolically interrelated. Augmentation of CES-1 increased microbial diversity in the later stages of augmentation when a strong microbial community selection of Acinetobacterparvus, Acinetobacterjohnsonii, Marinobacter manganoxydans, Verminephrobacter sp., and Arcobacter sp. occurred. Herein, there might be a possibility that the CES-1 augmentation could facilitate the indigenous microbial community successions so that the selected communities made the augmentation successful. The metagenomic analysis turned out to be a reasonable and powerful tool to provide with new insights and useful biomarkers for the complex environmental conditions, such as the full scale dye wastewater treatment system undergoing bioaugmentation.

Calf Diarrhea Caused by Prolonged Expansion of Autochthonous Gut Enterobacteriaceae and Their Lytic Bacteriophages.

Neonatal calf diarrhea is a common disease leading to a major economic loss for cattle producers worldwide. Several infectious and noninfectious factors are implicated in calf diarrhea, but disease control remains problematic because of the multifactorial etiology of the disease. Here, we conducted diagnostic multiplex PCR assay and meta-omics analysis (16S rRNA gene-based metataxonomics and untargeted transcriptional profiling) of rectal content of normal and diarrheic beef calves (n = 111). In the diarrheic calf gut, we detected both microbial compositional dysbiosis (i.e., increased abundances of the family Enterobacteriaceae members and their lytic bacteriophages) and functional dysbiosis (i.e., elevated levels of aerobic respiration and virulence potential). The calf diarrheic transcriptome mirrored the gene expression of the bovine host and was enriched in cellular pathways of sulfur metabolism, innate immunity, and gut motility. We then isolated 12 nontoxigenic Enterobacteriaceae strains from the gut of diarrheic calves. Feeding a strain mixture to preweaning mice resulted in a significantly higher level of fecal moisture content, with decreased body weight gain and shortened colon length. The presented findings suggest that gut inflammation followed by a prolonged expansion of nontoxigenic autochthonous Enterobacteriaceae contributes to the onset of diarrhea in preweaning animals.IMPORTANCE Calf diarrhea is the leading cause of death of neonatal calves worldwide. Several infectious and noninfectious factors are implicated in calf diarrhea, but disease control remains problematic because of the multifactorial etiology of the disease. The major finding of the current study centers around the observation of microbial compositional and functional dysbiosis in rectal samples from diarrheic calves. These results highlight the notion that gut inflammation followed by a prolonged expansion of autochthonous Enterobacteriaceae contributes to the onset of calf diarrhea. Moreover, this condition possibly potentiates the risk of invasion of notorious enteric pathogens, including Salmonella spp., and the emergence of inflammation-resistant (or antibiotic-resistant) microbiota via active horizontal gene transfer mediated by lytic bacteriophages.

Mobile resistome of human gut and pathogen drives anthropogenic bloom of antibiotic resistance.

BACKGROUND: The impact of human activities on the environmental resistome has been documented in many studies, but there remains the controversial question of whether the increased antibiotic resistance observed in anthropogenically impacted environments is just a result of contamination by resistant fecal microbes or is mediated by indigenous environmental organisms. Here, to determine exactly how anthropogenic influences shape the environmental resistome, we resolved the microbiome, resistome, and mobilome of the planktonic microbial communities along a single river, the Han, which spans a gradient of human activities. RESULTS: The bloom of antibiotic resistance genes (ARGs) was evident in the downstream regions and distinct successional dynamics of the river resistome occurred across the spatial continuum. We identified a number of widespread ARG sequences shared between the river, human gut, and pathogenic bacteria. These human-related ARGs were largely associated with mobile genetic elements rather than particular gut taxa and mainly responsible for anthropogenically driven bloom of the downstream river resistome. Furthermore, both sequence- and phenotype-based analyses revealed environmental relatives of clinically important proteobacteria as major carriers of these ARGs. CONCLUSIONS: Our results demonstrate a more nuanced view of the impact of anthropogenic activities on the river resistome: fecal contamination is present and allows the transmission of ARGs to the environmental resistome, but these mobile genes rather than resistant fecal bacteria proliferate in environmental relatives of their original hosts. Video abstract.

Antarctic tundra soil metagenome as useful natural resources of cold-active lignocelluolytic enzymes.

Lignocellulose composed of complex carbohydrates and aromatic heteropolymers is one of the principal materials for the production of renewable biofuels. Lignocellulose-degrading genes from cold-adapted bacteria have a potential to increase the productivity of biological treatment of lignocellulose biomass by providing a broad range of treatment temperatures. Antarctic soil metagenomes allow to access novel genes encoding for the cold-active lignocellulose-degrading enzymes, for biotechnological and industrial applications. Here, we investigated the metagenome targeting cold-adapted microbes in Antarctic organic matter-rich soil (KS 2-1) to mine lignolytic and celluloytic enzymes by performing single molecule, real-time metagenomic (SMRT) sequencing. In the assembled Antarctic metagenomic contigs with relative long reads, we found that 162 (1.42%) of total 11,436 genes were annotated as carbohydrate-active enzymes (CAZy). Actinobacteria, the dominant phylum in this soil's metagenome, possessed most of candidates of lignocellulose catabolic genes like glycoside hydrolase families (GH13, GH26, and GH5) and auxiliary activity families (AA7 and AA3). The predicted lignocellulose degradation pathways in Antarctic soil metagenome showed synergistic role of various CAZyme harboring bacterial genera including Streptomyces, Streptosporangium, and Amycolatopsis. From phylogenetic relationships with cellular and environmental enzymes, several genes having potential for participating in overall lignocellulose degradation were also found. The results indicated the presence of lignocellulose-degrading bacteria in Antarctic tundra soil and the potential benefits of the lignocelluolytic enzymes as candidates for cold-active enzymes which will be used for the future biofuel-production industry.

Characterization of a nifH-Harboring Bacterial Community in the Soil-Limited Gotjawal Forest.

Using a high-throughput metagenomic approach, we evaluated nifH-harboring bacterial communities and their assembly in the Gotjawal forest, which was naturally formed on basalt rocks with thin layer of soil. Significant differences in soil properties and community structure were observed in comparison with similar communities in various habitats, including other lava-formed forests (on Jeju Island and in Hawaii) and in regions with high humidity (Florida) or low temperatures (Alaska). nifH-harboring bacterial communities were found to assemble along gradients of environmental factors, particularly cation-exchange capacity. Unlike in other regions, in the Gotjawal forest, Paenibacillus and Clostridium, which belong to the phylum Firmicutes, were present in significantly higher proportion than in other regions. Network analysis suggested that much fewer co-occurrence relationships occurred in the Gotjawal forest than in other lava-formed forests. Our results indicate that the unique nifH-harboring bacterial community and its assembly in the Gotjawal forest are due to its distinctive soil properties, which has implications for microbial interactions and functional potentials.

Complete genome sequence of antibiotic and anticancer agent violacein producing Massilia sp. strain NR 4-1.

Massilia sp. NR 4-1 was a violacein producing strain newly isolated from topsoil under nutmeg tree, Torreya nucifera in Korean national monument Bijarim Forest. Violacein is a novel class of drug exhibiting anticancer and antibiotic activities originated from l-tryptophan. Here, we present the complete genome of Massilia sp. strain NR 4-1 of 6,361,416bp and total 5285 coding sequences (CDSs) including a complete violacein biosynthesis pathway, vioABCDE. The genome sequence of Massilia sp. NR 4-1 will provide stable and efficient biotechnological applications of violacein production.

Methylome Analysis of Two Xanthomonas spp. Using Single-Molecule Real-Time Sequencing.

Single-molecule real-time (SMRT) sequencing allows identification of methylated DNA bases and methylation patterns/motifs at the genome level. Using SMRT sequencing, diverse bacterial methylomes including those of Helicobacter pylori, Lactobacillus spp., and Escherichia coli have been determined, and previously unreported DNA methylation motifs have been identified. However, the methylomes of Xanthomonas species, which belong to the most important plant pathogenic bacterial genus, have not been documented. Here, we report the methylomes of Xanthomonas axonopodis pv. glycines (Xag) strain 8ra and X. campestris pv. vesicatoria (Xcv) strain 85-10. We identified N6-methyladenine (6mA) and N4-methylcytosine (4mC) modification in both genomes. In addition, we assigned putative DNA methylation motifs including previously unreported methylation motifs via REBASE and MotifMaker, and compared methylation patterns in both species. Although Xag and Xcv belong to the same genus, their methylation patterns were dramatically different. The number of 4mC DNA bases in Xag (66,682) was significantly higher (29 fold) than in Xcv (2,321). In contrast, the number of 6mA DNA bases (4,147) in Xag was comparable to the number in Xcv (5,491). Strikingly, there were no common or shared motifs in the 10 most frequently methylated motifs of both strains, indicating they possess unique species- or strain-specific methylation motifs. Among the 20 most frequent motifs from both strains, for 9 motifs at least 1% of the methylated bases were located in putative promoter regions. Methylome analysis by SMRT sequencing technology is the first step toward understanding the biology and functions of DNA methylation in this genus.