Uncovering microbiomes of the rice phyllosphere using long-read metagenomic sequencing.

The plant microbiome is crucial for plant growth, yet many important questions remain, such as the identification of specific bacterial species in plants, their genetic content, and location of these genes on chromosomes or plasmids. To gain insights into the genetic makeup of the rice-phyllosphere, we perform a metagenomic analysis using long-read sequences. Here, 1.8 Gb reads are assembled into 26,067 contigs including 142 circular sequences. Within these contigs, 669 complete 16S rRNA genes are clustered into 166 bacterial species, 121 of which show low identity (<97%) to defined sequences, suggesting novel species. The circular contigs contain novel chromosomes and a megaplasmid, and most of the smaller circular contigs are defined as novel plasmids or bacteriophages. One circular contig represents the complete chromosome of a difficult-to-culture bacterium Candidatus Saccharibacteria. Our findings demonstrate the efficacy of long-read-based metagenomics for profiling microbial communities and discovering novel sequences in plant-microbiome studies.

Bacteria can maintain rRNA operons solely on plasmids for hundreds of millions of years.

It is generally assumed that all bacteria must have at least one rRNA operon (rrn operon) on the chromosome, but some strains of the genera Aureimonas and Oecophyllibacter carry their sole rrn operon on a plasmid. However, other related strains and species have chromosomal rrn loci, suggesting that the exclusive presence of rrn operons on a plasmid is rare and unlikely to be stably maintained over long evolutionary periods. Here, we report the results of a systematic search for additional bacteria without chromosomal rrn operons. We find that at least four bacterial clades in the phyla Bacteroidota, Spirochaetota, and Pseudomonadota (Proteobacteria) lost chromosomal rrn operons independently. Remarkably, Persicobacteraceae have apparently maintained this peculiar genome organization for hundreds of millions of years. In our study, all the rrn-carrying plasmids in bacteria lacking chromosomal rrn loci possess replication initiator genes of the Rep_3 family. Furthermore, the lack of chromosomal rrn operons is associated with differences in copy numbers of rrn operons, plasmids, and chromosomal tRNA genes. Thus, our findings indicate that the absence of rrn loci in bacterial chromosomes can be stably maintained over long evolutionary periods.

Characterization of a carbapenem- and colistin-resistant Enterobacter cloacae carrying Tn6901 in bla NDM-1 genomic context.

We report a clinical strain of Enterobacter cloacae, PIMB10EC27, isolated in Vietnam in 2010 that was resistant to 21 of 26 tested antibiotics, including carbapenems (MICs >64 µg/mL) and colistin (MIC >128 µg/mL). The complete genome of strain PIMB10EC27 was sequenced by PacBio RSII and the Illumina Miseq system. Whole-genome analysis revealed that PIMB10EC27 contains a chromosome of the ST513 group (PIMBEC27, length 5,272,177 bp) and two plasmids, pEC27-1 of the IncX3 group (length 62,470 bp) and pEC27-2 of the IncHI1 group (length 84,602 bp). It also revealed that strain PIMB10EC27 carries 15 genes that confer resistance to at least 10 antibiotic groups. Particularly, the insertion of ISKpn19 and Tn6901 into the genomic context of bla NDM-1 was first identified and described. In another context, amino acid mutations G273D in PmrB and F515S in PmrC were first identified on the chromosome of PIMB10EC27, which may confer resistance to colistin in this strain.

Metaepigenomic analysis reveals the unexplored diversity of DNA methylation in an environmental prokaryotic community.

DNA methylation plays important roles in prokaryotes, and their genomic landscapes-prokaryotic epigenomes-have recently begun to be disclosed. However, our knowledge of prokaryotic methylation systems is focused on those of culturable microbes, which are rare in nature. Here, we used single-molecule real-time and circular consensus sequencing techniques to reveal the 'metaepigenomes' of a microbial community in the largest lake in Japan, Lake Biwa. We reconstructed 19 draft genomes from diverse bacterial and archaeal groups, most of which are yet to be cultured. The analysis of DNA chemical modifications in those genomes revealed 22 methylated motifs, nine of which were novel. We identified methyltransferase genes likely responsible for methylation of the novel motifs, and confirmed the catalytic specificities of four of them via transformation experiments using synthetic genes. Our study highlights metaepigenomics as a powerful approach for identification of the vast unexplored variety of prokaryotic DNA methylation systems in nature.

Complete Genome Sequence of Methylobacterium sp. Strain AMS5, an Isolate from a Soybean Stem.

Nonrhizobial Methylobacterium spp. inhabit the phyllosphere of a wide variety of plants. We report here the complete genome sequence of Methylobacterium sp. AMS5, which was isolated from a soybean stem. The information is useful for understanding the molecular mechanisms of the interaction between nonrhizobial Methylobacterium spp. and plants.

Sulfur Fertilization Changes the Community Structure of Rice Root-, and Soil- Associated Bacteria.

Under paddy field conditions, biological sulfur oxidation occurs in the oxidized surface soil layer and rhizosphere, in which oxygen leaks from the aerenchyma system of rice plants. In the present study, we examined community shifts in sulfur-oxidizing bacteria associated with the oxidized surface soil layer and rice roots under different sulfur fertilization conditions based on the 16S ribosomal RNA (rRNA) gene in order to explore the existence of oligotrophic sulfur-oxidizing bacteria in the paddy rice ecosystem. Rice plants were grown in pots with no fertilization (control) or CaCO3 or CaSO4 fertilization. A principal-coordinates analysis (PCoA) showed that CaSO4 fertilization markedly affected bacterial communities associated with rice roots and soil, whereas no significant differences were observed in plant growth among the fertilizer treatments examined. In rice roots, the relative abundance of Acidobacteria, Alphaproteobacteria, Gammaproteobacteria, and TM7 was significantly higher in CaSO4-fertilized pots than in control pots. Alphaproteobacteria, Bradyrhizobiaceae, and Methylocystaceae members were significantly more abundant in CaSO4-fertilized roots than in control roots. On the other hand, the abundance of Actinobacteria and Proteobacteria was lower in CaSO4-fertilized soil than in control soil. These results indicate that the bacteria associated with rice roots and soil responded to the sulfur amendment, suggesting that more diverse bacteria are involved in sulfur oxidation in the rice paddy ecosystem than previously considered.

Bacterial clade with the ribosomal RNA operon on a small plasmid rather than the chromosome.

rRNA is essential for life because of its functional importance in protein synthesis. The rRNA (rrn) operon encoding 16S, 23S, and 5S rRNAs is located on the "main" chromosome in all bacteria documented to date and is frequently used as a marker of chromosomes. Here, our genome analysis of a plant-associated alphaproteobacterium, Aureimonas sp. AU20, indicates that this strain has its sole rrn operon on a small (9.4 kb), high-copy-number replicon. We designated this unusual replicon carrying the rrn operon on the background of an rrn-lacking chromosome (RLC) as the rrn-plasmid. Four of 12 strains close to AU20 also had this RLC/rrn-plasmid organization. Phylogenetic analysis showed that those strains having the RLC/rrn-plasmid organization represented one clade within the genus Aureimonas. Our finding introduces a previously unaddressed viewpoint into studies of genetics, genomics, and evolution in microbiology and biology in general.

Symbiosis island shuffling with abundant insertion sequences in the genomes of extra-slow-growing strains of soybean bradyrhizobia.

Extra-slow-growing bradyrhizobia from root nodules of field-grown soybeans harbor abundant insertion sequences (ISs) and are termed highly reiterated sequence-possessing (HRS) strains. We analyzed the genome organization of HRS strains with the focus on IS distribution and symbiosis island structure. Using pulsed-field gel electrophoresis, we consistently detected several plasmids (0.07 to 0.4 Mb) in the HRS strains (NK5, NK6, USDA135, 2281, USDA123, and T2), whereas no plasmids were detected in the non-HRS strain USDA110. The chromosomes of the six HRS strains (9.7 to 10.7 Mb) were larger than that of USDA110 (9.1 Mb). Using MiSeq sequences of 6 HRS and 17 non-HRS strains mapped to the USDA110 genome, we found that the copy numbers of ISRj1, ISRj2, ISFK1, IS1632, ISB27, ISBj8, and IS1631 were markedly higher in HRS strains. Whole-genome sequencing showed that the HRS strain NK6 had four small plasmids (136 to 212 kb) and a large chromosome (9,780 kb). Strong colinearity was found between 7.4-Mb core regions of the NK6 and USDA110 chromosomes. USDA110 symbiosis islands corresponded mainly to five small regions (S1 to S5) within two variable regions, V1 (0.8 Mb) and V2 (1.6 Mb), of the NK6 chromosome. The USDA110 nif gene cluster (nifDKENXSBZHQW-fixBCX) was split into two regions, S2 and S3, where ISRj1-mediated rearrangement occurred between nifS and nifB. ISs were also scattered in NK6 core regions, and ISRj1 insertion often disrupted some genes important for survival and environmental responses. These results suggest that HRS strains of soybean bradyrhizobia were subjected to IS-mediated symbiosis island shuffling and core genome degradation.

Metaproteomic identification of diazotrophic methanotrophs and their localization in root tissues of field-grown rice plants.

In a previous study by our group, CH4 oxidation and N2 fixation were simultaneously activated in the roots of wild-type rice plants in a paddy field with no N input; both processes are likely controlled by a rice gene for microbial symbiosis. The present study examined which microorganisms in rice roots were responsible for CH4 oxidation and N2 fixation under the field conditions. Metaproteomic analysis of root-associated bacteria from field-grown rice (Oryza sativa Nipponbare) revealed that nitrogenase complex-containing nitrogenase reductase (NifH) and the alpha subunit (NifD) and beta subunit (NifK) of dinitrogenase were mainly derived from type II methanotrophic bacteria of the family Methylocystaceae, including Methylosinus spp. Minor nitrogenase proteins such as Methylocella, Bradyrhizobium, Rhodopseudomonas, and Anaeromyxobacter were also detected. Methane monooxygenase proteins (PmoCBA and MmoXYZCBG) were detected in the same bacterial group of the Methylocystaceae. Because these results indicated that Methylocystaceae members mediate both CH4 oxidation and N2 fixation, we examined their localization in rice tissues by using catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH). The methanotrophs were localized around the epidermal cells and vascular cylinder in the root tissues of the field-grown rice plants. Our metaproteomics and CARD-FISH results suggest that CH4 oxidation and N2 fixation are performed mainly by type II methanotrophs of the Methylocystaceae, including Methylosinus spp., inhabiting the vascular bundles and epidermal cells of rice roots.

Impact of Azospirillum sp. B510 inoculation on rice-associated bacterial communities in a paddy field.

Rice seedlings were inoculated with Azospirillum sp. B510 and transplanted into a paddy field. Growth in terms of tiller numbers and shoot length was significantly increased by inoculation. Principal-coordinates analysis of rice bacterial communities using the 16S rRNA gene showed no overall change from B510 inoculation. However, the abundance of Veillonellaceae and Aurantimonas significantly increased in the base and shoots, respectively, of B510-inoculated plants. The abundance of Azospirillum did not differ between B510-inoculated and uninoculated plants (0.02-0.50%). These results indicate that the application of Azospirillum sp. B510 not only enhanced rice growth, but also affected minor rice-associated bacteria.

Complete genome sequence of Bradyrhizobium sp. S23321: insights into symbiosis evolution in soil oligotrophs.

Bradyrhizobium sp. S23321 is an oligotrophic bacterium isolated from paddy field soil. Although S23321 is phylogenetically close to Bradyrhizobium japonicum USDA110, a legume symbiont, it is unable to induce root nodules in siratro, a legume often used for testing Nod factor-dependent nodulation. The genome of S23321 is a single circular chromosome, 7,231,841 bp in length, with an average GC content of 64.3%. The genome contains 6,898 potential protein-encoding genes, one set of rRNA genes, and 45 tRNA genes. Comparison of the genome structure between S23321 and USDA110 showed strong colinearity; however, the symbiosis islands present in USDA110 were absent in S23321, whose genome lacked a chaperonin gene cluster (groELS3) for symbiosis regulation found in USDA110. A comparison of sequences around the tRNA-Val gene strongly suggested that S23321 contains an ancestral-type genome that precedes the acquisition of a symbiosis island by horizontal gene transfer. Although S23321 contains a nif (nitrogen fixation) gene cluster, the organization, homology, and phylogeny of the genes in this cluster were more similar to those of photosynthetic bradyrhizobia ORS278 and BTAi1 than to those on the symbiosis island of USDA110. In addition, we found genes encoding a complete photosynthetic system, many ABC transporters for amino acids and oligopeptides, two types (polar and lateral) of flagella, multiple respiratory chains, and a system for lignin monomer catabolism in the S23321 genome. These features suggest that S23321 is able to adapt to a wide range of environments, probably including low-nutrient conditions, with multiple survival strategies in soil and rhizosphere.

Isolation and genetic characterization of Aurantimonas and Methylobacterium strains from stems of hypernodulated soybeans.

The aims of this study were to isolate Aurantimonas and Methylobacterium strains that responded to soybean nodulation phenotypes and nitrogen fertilization rates in a previous culture-independent analysis (Ikeda et al. ISME J. 4:315-326, 2010). Two strategies were adopted for isolation from enriched bacterial cells prepared from stems of field-grown, hypernodulated soybeans: PCR-assisted isolation for Aurantimonas and selective cultivation for Methylobacterium. Thirteen of 768 isolates cultivated on Nutrient Agar medium were identified as Aurantimonas by colony PCR specific for Aurantimonas and 16S rRNA gene sequencing. Meanwhile, among 187 isolates on methanol-containing agar media, 126 were identified by 16S rRNA gene sequences as Methylobacterium. A clustering analysis (>99% identity) of the 16S rRNA gene sequences for the combined datasets of the present and previous studies revealed 4 and 8 operational taxonomic units (OTUs) for Aurantimonas and Methylobacterium, respectively, and showed the successful isolation of target bacteria for these two groups. ERIC- and BOX-PCR showed the genomic uniformity of the target isolates. In addition, phylogenetic analyses of Aurantimonas revealed a phyllosphere-specific cluster in the genus. The isolates obtained in the present study will be useful for revealing unknown legume-microbe interactions in relation to the autoregulation of nodulation.

Autoregulation of nodulation interferes with impacts of nitrogen fertilization levels on the leaf-associated bacterial community in soybeans.

The diversities leaf-associated bacteria on nonnodulated (Nod(-)), wild-type nodulated (Nod(+)), and hypernodulated (Nod(++)) soybeans were evaluated by clone library analyses of the 16S rRNA gene. To analyze the impact of nitrogen fertilization on the bacterial leaf community, soybeans were treated with standard nitrogen (SN) (15 kg N ha(-1)) or heavy nitrogen (HN) (615 kg N ha(-1)) fertilization. Under SN fertilization, the relative abundance of Alphaproteobacteria was significantly higher in Nod(-) and Nod(++) soybeans (82% to 96%) than in Nod(+) soybeans (54%). The community structure of leaf-associated bacteria in Nod(+) soybeans was almost unaffected by the levels of nitrogen fertilization. However, differences were visible in Nod(-) and Nod(++) soybeans. HN fertilization drastically decreased the relative abundance of Alphaproteobacteria in Nod(-) and Nod(++) soybeans (46% to 76%) and, conversely, increased those of Gammaproteobacteria and Firmicutes in these mutant soybeans. In the Alphaproteobacteria, cluster analyses identified two operational taxonomic units (OTUs) (Aurantimonas sp. and Methylobacterium sp.) that were especially sensitive to nodulation phenotypes under SN fertilization and to nitrogen fertilization levels. Arbuscular mycorrhizal infection was not observed on the root tissues examined, presumably due to the rotation of paddy and upland fields. These results suggest that a subpopulation of leaf-associated bacteria in wild-type Nod(+) soybeans is controlled in similar ways through the systemic regulation of autoregulation of nodulation, which interferes with the impacts of N levels on the bacterial community of soybean leaves.

Community- and genome-based views of plant-associated bacteria: plant-bacterial interactions in soybean and rice.

Diverse microorganisms are living as endophytes in plant tissues and as epiphytes on plant surfaces in nature. Questions about driving forces shaping the microbial community associated with plants remain unanswered. Because legumes developed systems to attain endosymbioses with rhizobia as well as mycorrhizae during their evolution, the above questions can be addressed using legume mutants relevant to genes for symbiosis. Analytical methods for the microbial community have recently been advanced by enrichment procedures of plant-associated microbes and culture-independent analyses targeting the small subunit of rRNA in microbial ecology. In this review, we first deal with interdisciplinary works on the global diversity of bacteria associated with field-grown soybeans with different nodulation genotypes and nitrogen application. A subpopulation of Proteobacteria in aerial parts of soybean shoots was likely to be regulated through both the autoregulation system for plant-rhizobium symbiosis and the nitrogen signaling pathway, suggesting that legumes accommodate a taxonomically characteristic microbial community through unknown plant-microbe communications. In addition to the community views, we then show multiphasic analysis of a beneficial rice endophyte for comparative bacterial genomics and plant responses. The significance and perspectives of community- and genome-based approaches are discussed to achieve a better understanding of plant-microbe interactions.