A healthy plant of Piper crocatum was collected from Bogor, Indonesia. The plant was identified by Dr. Joeni Setijo Rahajoe (Herbarium Bogoriense, BRIN) by morphological comparison with an authenticated reference specimen of the same species held at the Herbarium Bogoriense (BO), Directorate of Scientific Collection Management, National Research and Innovation Agency (BRIN), Cibinong, Indonesia (specimen accession no. BO-1475717). Plant tissues (leaves, stems, and roots) were thoroughly washed under running tap water, followed by surface sterilisation using 70% ethanol for two minutes, 5% sodium hypochlorite for five minutes, 70% ethanol for 30 seconds, and rinsed four times with sterile distilled water. The sterilised tissues were used for further processing. Tissues were macerated and segmented prior to enrichment in tryptic soy broth (TSB) for four hours in a shaker incubator with 180 RPM at 37 °C. The suspension was then spread onto tryptic soy agar (TSA) supplemented with amphotericin B, and incubated in a 37 °C incubator for 20 hours for bacterial isolation. Distinct bacterial colonies were subcultured onto fresh TSA plates several times until pure isolates were obtained. The isolate used in this study was designated PCL31, preserved at −80 °C in 25% ( v/v) glycerol, and deposited in the i3L University microbial culture collection under the strain designation PCL31. The strain is available from the corresponding author (Mega Safithri; safithri@apps.ipb.ac.id) upon reasonable request for non-commercial research use, subject to institutional approval and a completed material transfer agreement. To confirm endophytic origin, 100 μL of the final rinse water was plated on TSA in parallel; no colonies were recovered after 48 h at 37 °C.

The pure isolates were cultured in tryptic soy broth (TSB) for 24 hours in a shaker incubator at 180 RPM, 37 °C. Then, the pellet of the pure isolate obtained from the cultured broth was centrifuged at 4000 × g for 5 minutes. The pellet was further processed through lysis and purification to isolate genomic DNA using the spin-column method from the gSYNC™ DNA extraction kit (Geneaid Biotech Ltd., Taiwan). Extracted DNA samples were evaluated using gel electrophoresis (1% tris borate-EDTA agarose w/v) for integrity, and NanoDrop Lite Plus microvolume spectrophotometer (Thermo Fisher Scientific, USA) as well as Quantus fluorometer (Promega, USA) equipped with the QuantiFluor ONE dsDNA system, for concentration and purity.

The DNA library was prepared for PromethION ONT using the SQK-NBD114.24 ligation sequencing kit according to the manufacturer’s recommendation (Oxford Nanopore Technology, UK). Kits from NEBNext, FFPE Repair Mix (NEB M6630, UK) and Ultra II End Repair Module (NEB E7546, UK) were employed in the library preparation. In brief, total gDNA was repaired using an end-prep enzyme mix (NEB M6630 & NEB E7546), generating DNA with 5′ phosphorylated, 3′ dA-tailed ends. Barcodes were ligated with an ONT-compatible adapter. The library was quantified with a Quantus Fluorometer (Promega, USA) before loading onto the FLO-PRO114M flow cell. Sequencing control was conducted in MinKNOW v23.04.6 following the default sequencing run from ONT (72 h). The output raw reads were set to POD5 for a separate basecalling process.

The raw POD5 data generated was basecalled using Dorado v0.8.2 super accuracy model and the resulting BAM output was converted to FASTQ using Samtools v1.21 ( Danecek et al., 2021). Next, the raw FASTQ file was assessed using NanoPlot v1.42.0 to analyze long-read ONT whole genome sequences ( De Coster & Rademakers, 2023). A pre-filtering step with Filtlong v0.2.1 ( https://github.com/rrwick/Filtlong) with min-length 1 kb and min-Q 10 did not improve the assembly N50, BUSCO completeness, or estimated coverage in pilot tests, and was therefore omitted from the final pipeline. Reads were then assembled de novo using Flye v2.9.5 ( Kolmogorov et al., 2019) to generate a whole-genome sequencing draft of the isolate in FASTA format. Parameters included --nano-hq --genome-size 5 m --iterations 3 --threads 4. The draft assembly was then polished using Medaka consensus v2.0.1 to enhance sequence accuracy and correct residual errors. The quality of the polished assembly, including total length, GC content, and N50, was evaluated using the quality assessment tool for genome assembly (QUAST) v5.3.0 ( Gurevich et al., 2013). Genome completeness of the polished assembly was assessed using the benchmarking universal single-copy orthologs (BUSCO) v5.8.0 ( Manni et al., 2021) against the enterobacterales_odb10 lineage dataset (n = 440 single-copy orthologs), where assemblies with >95% completeness and < 5% contamination were considered high-quality.

Assembly graph topology was inspected from the Flye assembly_info.txt output. Plasmid replicons were screened with PlasmidFinder v2.1 against the Enterobacterales database (≥95% identity, ≥60% coverage; Carattoli et al., 2014), and molecule classification (chromosome vs. plasmid) and mobility typing were performed with MOB-suite v3.1.9 using mob_recon defaults ( Robertson & Nash, 2018). Per-contig statistics (length, GC content, coverage, circularity, and replicon/molecule assignment) are provided in Underlying Data ( Sidhartha et al., 2026).

Whole genome-based taxonomic analysis was carried out on the type (strain) genome server (TYGS) ( Meier-Kolthoff & Göker, 2019) by submitting the polished genome assembly and utilizing an unrestricted query. The TYGS webserver performed pairwise comparisons of genome sequences against the TYGS database to yield species identification, DNA-DNA hybridization (dDDH) score, and phylogenetic inference. A genome BLAST distance phylogeny (GBDP) tree was also generated.

The result was further confirmed using FastANI v1.34 ( Jain et al., 2018) against the reference Kluyvera ascorbata complete-genome assembly NZ_CP096201.1 (strain SK) retrieved from NCBI RefSeq; the species type strain is ATCC 33433 (assembly JMPL00000000), used for the TYGS dDDH comparison above. The average nucleotide identity (ANI) was calculated by the mean percentage of identical nucleotides across all shared regions between the obtained and reference sequences. Results were interpreted based on the established ANI threshold (>90%).

Genome annotation was performed using Bakta v1.9.4 ( Schwengers et al., 2021) utilizing Bakta database v5.0_2023-03-20 and AMRFinderPlus database v3.12–2024-05-02.2.

BGCs were identified using antiSMASH 8.0 ( Blin et al., 2025) in relaxed detection mode with KnownClusterBlast, ClusterBlast, SubClusterBlast, MIBiG comparison, and ClusterCompare modules enabled; references were drawn from MIBiG v4.0. Similarity scores represent the aggregated region-level score from antiSMASH (range 0–1), where higher values reflect greater gene content and synteny conservation with the reference cluster; scores ≥0.80 typically indicate clusters encoding homologues of the reference compound, whereas scores ≤0.50 were interpreted as putative, low-confidence matches.

Antimicrobial resistance gene (ARG) profiling was done using the Comprehensive Antibiotic Resistance Database (CARD) v4.0.1 Resistance Gene Identifier (RGI) v6.0.5( Alcock et al., 2023), selecting for Perfect and Strict hits, and specifying the input as high quality/coverage assembly.

The assembled genome was screened for CRISPR arrays and Cas operons using the CRISPRCasFinder ( Couvin et al., 2018) utilizing the default settings.

Genome map was visualized using Proksee ( Grant et al., 2023) to generate a circular genome map depicting sequence features, contig boundaries, and genomic architecture. BGCs and ARGs coordinates from their respective tools were manually converted into a location map as.bed file and included as separate forward- and reverse-strand tracks.

TYGS species identification using the whole genome assembly suggested Kluyvera ascorbata, as indicated by the dDDH value of 88.6% with the reference strain Kluyvera ascorbata ATCC 33433 ( Figure 1). The identity was confirmed with the ANI score of 98.6% between the sample and Kluyvera ascorbata reference sequence (NCBI RefSeq NZ_CP096201.1). The next-closest TYGS match was Kluyvera cryocrescens (dDDH = 54.3%), well below the 70% species-delineation threshold, confirming species-level assignment to Kluyvera ascorbata.

Genome mining with antiSMASH predicted four putative biosynthetic gene clusters (BGCs) in the assembly ( Table 3). The 53.8 kb NRP-metallophore/NRPS cluster on contig_1 corresponded to the canonical enterobactin biosynthesis locus of Escherichia coli K-12 (MIBiG BGC0002476; similarity = 0.83), consistent with the presence of this deeply conserved siderophore pathway across Enterobacteriaceae. The remaining three clusters — a RiPP-like, a terpene-precursor, and an azole-containing RiPP region — returned only low-similarity matches (scores ≤0.50) to MIBiG references from phylogenetically distant taxa, and therefore represent clusters with no close match in MIBiG. Low similarity may reflect either true biosynthetic novelty or under-representation of Kluyvera-associated clusters in current reference databases, requiring further metabolomic validation.

CARD RGI screening of the Kluyvera ascorbata assembly returned 26 AMR determinants (1 Perfect, 25 Strict). The sole Perfect hit was bla CTX-M-2, located on the chromosomal contig_4 and with a gene length of 873 bp, which is consistent with Kluyvera ascorbata being the recognised chromosomal progenitor of the CTX-M-2 β-lactamase lineage. Strict hits were dominated by efflux-pump components, including ten resistance-nodulation-cell division (RND)-family genes ( acrA, acrB, acrD, adeF, baeR, CRP, marA, mdtB, mdtC, rsmA) and a MarR-overexpression variant of the AcrAB–TolC system, three major facilitator superfamily (MFS) pumps ( emrR, KpnH, leuO), two small multidrug resistance (SMR) pumps ( KpnE, KpnF), and one ATP-binding cassette (ABC) pump ( msbA). Target-alteration determinants comprised vanL, ArnT, PmrF, two paralogues of EF-Tu carrying R234F, UhpT (E350Q), and PBP3 (D350N/S357N), and an additional inactivation enzyme, FosA3, was also recovered. Collectively, the predicted resistome spans β-lactams, fluoroquinolones, aminoglycosides, macrolides, tetracyclines, glycopeptides, fosfomycin, polymyxins, aminocoumarins, phenicols, rifamycins, and disinfectants. The Strict-hit complement comprises predominantly intrinsic Enterobacteriaceae efflux systems and chromosomal target-modification genes whose contribution to clinical resistance phenotype cannot be inferred from sequence alone. Phenotypic resistance testing (MIC) was outside the scope of this Genome Note and is identified here as a limitation.

No CRISPR/Cas component was detected using the CRISPR/Cas Finder tool. The schematic of a circular map of the genome, including the annotated BGCs and ARGs, is presented in Figure 2.

In summary, this Genome Note describes a high-quality draft genome of an endophytic Kluyvera ascorbata strain PCL31 isolated from the medicinal plant Piper crocatum, together with its annotation, predicted biosynthetic gene clusters, and predicted antimicrobial resistance gene complement. The genome is made publicly available as a resource to support future functional, metabolomic, and antimicrobial resistance studies; detailed analysis and interpretation of these features are beyond the scope of this Genome Note.