Rice seeds (Variety: Pusa Basmati 1121) were collected from rice field located in Kharar, Punjab (30.6755° N, 76.6723° E) for 3 years and processed in three independent batches for seed sterilization and isolation of endophytic bacteria. Endophytic bacteria were isolated from surface sterilized rice seeds by following the process: The hulls were removed from rice seeds using sterilized forceps, and the seeds (5 g) were put in sterile falcon tubes and washed with sterilized water for 1 min and then with 1% sodium hypochlorite solution for 5 min. The seeds were again washed with 75% ethanol for 1 min. After another wash with sterilized water five times, the surface sterilized rice seeds were crushed in sterile mortar and pestle and suspended in sterile saline solution (0.85% NaCl). The seeds suspension was incubated for 2 h at 28°C under shaking condition. After that, 100 μl of each of Direct, 10⁻¹,10⁻², 10⁻³, and 10⁻⁴ dilution in sterile saline was plated in duplicates onto Nutrient agar (NA); King's medium B (KMB); Glucose yeast chalk agar (GYCA); Tryptic soy agar (TSA); Peptone sucrose agar (PSA) supplemented with 0.01% cycloheximide. The confirmation of surface sterilization was conducted by spreading the last water wash as well as placing the washed rice seeds onto different media plates. Further, the isolates were identified by 16S rRNA gene PCR using 16S rRNA universal primers (27F and 1492R) and 16S rRNA gene sequence was used for characterization using EzTaxon server (http://www.ezbiocloud.net/eztaxon; Kim et al., 2012) prior to whole genome sequencing. For preservation, 15% glycerol stocks of the pure culture of each isolate was prepared and maintained at −80°C.

Endophytic bacteria belong to Firmicutes phylum were revived from −80°C stocks onto NA plates and genomic DNA was isolated using Zymo ZR- Fungal/bacterial DNA isolation kit as per the instruction manual. The quality of the genomic DNA was assessed using agarose gel electrophoresis. DNA concentrations were estimated using a Nanodrop spectrophotometer ND-100 (Thermo Fisher Scientific, USA) and Qubit 2.0 fluorometer (Invitrogen, USA), with a double-stranded DNA High Sensitivity (dsDNA HS) Assay Kit. Both the ratio A260/280 and gel electrophoresis were used to ascertain the quality and purity of DNA samples. The input of 1 μg of genomic DNA from each sample was taken and standard protocol for the Nextera XT DNA (Illumina, San Diego, CA) sample preparation kit was used for library construction. The purified fragmented DNA was used as a template for a limited cycle PCR using Nextera primers and index adaptors. Cluster generation and sequencing of libraries were performed on the Illumina MiSeq platform (Illumina, San Diego, CA) with a 2 × 250 paired-end run.

Demultiplexing, FASTQ generation and adapter trimming in raw sequence reads were automatically performed by Illumina- MiSeq software. The paired-end raw reads containing FASTQ files were assembled into contigs. Total number of library reads, number of contigs, genome size, G+C content, and total coverage were analyzed using CLC Genomics Workbench software version 7.0.3. Complete 16S rRNA sequences were extracted from the assembled genomes using RNAmmer 1.2 server (Lagesen et al., 2007) and characterized by EzTaxon server. tRNA was calculated by tRNAScan-SE. Sequences were annotated using NCBI Prokaryotic Genome Annotation Pipeline (http://www.ncbi.nlm.nih.gov/genome/annotation_prok/) and Rapid Annotation System Technology (RAST) pipeline (Aziz et al., 2008).

Taxonomic relationship among the endophytic bacterial isolates was deduced using 16S rRNA gene which is considered as the important tool for taxonomic and phylogenetic analysis (Mizrahi-Man et al., 2013). In the era of genome based taxonomy, the 16S rRNA gene is still used for preliminary bacterial typing because it is present in at least one copy in every bacterial genome and its conserved regions enable simple sample identification (Land et al., 2015). Therefore, firstly, the 16S rRNA sequences were used for characterization using EzTaxon server. Further, for genome similarity assessment, BLAST-based average nucleotide identity (ANIb) and Genome to Genome Distance calculator or digital DNA-DNA hybridization (dDDH) values were used. Pairwise ANI was calculated using JSpecies (Richter and Rosselló-Móra, 2009) and dDDH (Auch et al., 2010) was calculated using web tool GGDC 2.0.

The assembled genome sequences (in FASTA format) of all 21 Rice seed bacterial endophytes is deposited in NCBI GenBank. Assembly statistics and accession numbers of the sequenced genomes are mentioned in Table 1. Figshare link to download genomes in FASTA format for single point access: https://figshare.com/s/6deae778bbdb90c6563c.

RAST annotated genomes was used to search important genes responsible for plant adaptation and plant protection.

A total of 21 bacterial endophytes isolated from rice seeds (Variety: Pusa Basmati 1121). The isolates were identified based on 16SrRNA gene for initial characterization prior to whole genome sequencing and their identification using EzTaxon server is summarized in Table S1.

Sequence data generated for each strain (Table 1) was de novo assembled with coverage ranging from 30x to 76x. Analysis based on complete 16SrRNA gene sequence extracted from the whole genome sequences assigned the bacterial isolates into 2 distinct genera and 6 species (Table S1). To gain better taxonomic resolution, all the 21 bacterial isolates have been validated by species delineation cut-off of >95% ANI values and >70% dDDH values (Richter and Rosselló-Móra, 2009; Auch et al., 2010). Genome sequence of reference strains were taken from NCBI database. The ANI value heatmap and dDDH values of 21 endophytes genomes with Type strain/reference strain genome are depicted in Figure 1 and Table S2, respectively. Description of the strains assigned to 2 different genera is provided as follows:

Staphylococcus is cocci shaped gram-positive member of phylum Firmicutes and 14 different strains were selected for sequencing from this genus belonging to 3 different year rice seeds sample lots. 16S rRNA sequences have assigned them to four different species with more than 99% identity to Staphylococcus warneri (RE1.51 and RE1.52), Staphylococcus hominis (RE2.7, RE2.8, and RE2.10), Staphylococcus equorum (RE2.24, RE2.35, and RE2.40), and Staphylococcus cohnii (SE3.10, SE4.1, SE4.2, SE4.3, SE4.4, and SE4.5). Moreover, ANI values also supported the species distinction between the four groups as they are above the cut-off (>95%) for species delineation (Figure 1). Further dDDH values of all the strains were more than 70% which further confirms their species identity (Table S2). Both ANI and dDDH resolved Staphylococcus species up to subspecies/strain level which were not possible with 16S rRNA gene sequences based analysis. This suggests that atleast 4 species of Staphylococcus inhabits rice seeds.

Bacillus is a genus of gram-positive, rod-shaped bacteria and a member of the phylum Firmicutes. Seven isolates from different years of rice seed lots belong to two different species, Bacillus subtilis (RE1.1, RE1.3, and RE1.4), and Bacillus licheniformis (RE1.58, RE1.59, RE1.60, and RE1.61) with more than 99% identity of 16SrRNA gene sequence identity. ANI and dDDH values of RE1.58, RE1.59, RE1.60, and RE1.61 of >99 and >97% with type strain genomes unequivocally establish their species identity as B. licheniformis (Figure 1; Table S2). However, in case of RE1.1, RE1.3, and RE1.4, ANI value with type strain of B. subtilis was 95.2% that is close to borderline of species demarcation, and dDDH values of 48.9% also showed that it is closely related to B. subtilis and may belongs to separate species.

From RAST analysis, it was observed that more than 100 genes in the genomes of B. subtilis and B. licheniformis strains whereas more than 70 genes were present in genomes of S. warneri, S. hominis, S. equorum, and S. cohnii were identified that are responsible for stress tolerance. Majority of these genes were responsible for protection against oxidative stress and rest were encoded to cope with heat and osmotic stresses. Moreover, all the endophytes genomes of Staphylococcus and Bacillus species contain genes responsible for auxin biosynthesis and siderophores production. Full details of these genes are available in Table S3. It was reported that members of genus Bacillus are most commonly found as endophytic bacteria in plants. Moreover, they play important role as a biocontrol agent against phytopathogens, stimulates plant growth, and also produce plant growth hormones such as auxin and gibberellin, as well as able to ameliorate drought stress (Forchetti et al., 2007). There are also several reports in literature where members of Staphylococcus were documented as endophytes such as in rice seed (Chaudhry and Patil, 2016), maize kernels (Liu et al., 2012), grapevine, and hybrid spruce (Collins et al., 2004).

In the present study, our goal was to generate genomic data resource of endophytic bacteria associated with healthy rice seeds (var. Pusa Basmati 1121) from India. As representatives, dominant members belonging to Firmicutes phylum were selected. Our genomic dataset which comprised of Staphylococcus and Bacillus genera and their phylogenomic analysis reported here will serve to catalyze future studies by providing a new lens to study their endophytic lifestyle and will help in deciphering novel biological insights of endophytic bacteria-host plant relationships.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.