
==== Front
Data Brief
Data Brief
Data in Brief
2352-3409
Elsevier

S2352-3409(24)00874-6
10.1016/j.dib.2024.110911
110911
Data Article
Whole genome sequence data of ZEN-degrading strain Levilactobacillus brevis PYN10_6_2 isolated from Tenebrio molitor larval feces
Zhao Mengru ab
Guo Baoyuan ab
Jiang Shuzhen c
Wang Yang wy@ags.ac.cn
b⁎
a School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
b Institute of Grain and Oil Quality and Safety, Academy of National Food and Strategic Reserves Administration, Beijing 100834, China
c College of Animal Science and Technology, Shandong Agricultural University, Tai'an 271018, China
⁎ Corresponding author. wy@ags.ac.cn
04 9 2024
12 2024
04 9 2024
57 11091121 6 2024
27 8 2024
30 8 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by-nc/4.0/ This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/).
Levilactobacillus brevis PYN10_6_2, a lactic acid bacterial strain previously isolated from Tenebrio molitor larval feces, possesses the ability to convert zearalenone (ZEN) to α-/β-Zearalenol (α-/β-ZEL). However, the genes involved in the ZEN reduction reaction and the biosafety of this strain remain unknown. In this study, we sequenced, assembled, and annotated the whole genome of L. brevis PYN10_6_2. Genomic sequencing was conducted using short-read sequencing on the Illumina HiSeq X Ten platform and long-read sequencing on the PacBio RS II Single Molecule Real-Time (SMRT) platform. The assembled genome consisted of one circular chromosome, four circular plasmids, with a total size of 2,745,725 bp and a G + C content of 45.52 %. Annotation identified 2,660 coding sequences, 5 rRNAs, 66 tRNAs, and a single CRISPR locus. Average nucleotide identity (ANI) between L. brevis PYN10_6_2 and L. brevis DSM 20054T yielded a value of 98.94 %. Further in-depth analysis revealed 182 antibiotic resistance genes, 237 putative virulence genes, 2 prophages, and 10 genomic islands. Additionally, functional annotation through COG and KEGG databases revealed the presence of three genes encoding 3α- and 3β-hydroxysteroid dehydrogenase (3α-/3β-HSD) within the bacterial chromosome. This comprehensive genomic characterization provides valuable insights into the genetic basis of L. brevis PYN10_6_2′s ZEN-reducing ability and its biosafety profile.

Keywords

Levilactobacillus brevis
Genome sequencing
Zearalenone
Hydroxysteroid dehydrogenase
Biosafety
==== Body
pmcSpecifications TableSubject	Microbiology: Applied Microbiology.	
Specific subject area	Molecular microbiology and Genomics	
Type of data	Table, Figures	
Data collection	Genomic DNA was extracted using the Wizard® Genomic DNA Purification Kit (Promega). Short-read sequencing libraries were prepared with the NEXTflex™ Rapid DNA-Seq Kit and sequenced on an Illumina HiSeq X Ten platform (2 × 150 bp paired-end reads). Long-read sequencing was performed on a PacBio RS II Single Molecule Real Time (SMRT) platform. Using Unicycler v0.4.8 for long-read assembly, followed by sequence polishing with Pilon v1.22. Genome annotation was performed with the Prokaryotic Genome Annotation Pipeline (PGAP) and functionally characterized using Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Antibiotic resistance genes were identified using the Comprehensive Antibiotic Resistance Database (CARD) v3.2.6. Virulence factors were predicted based on homology searches against the Virulence Factors Database (VFDB). Genomic islands were identified using IslandViewer v1.2 and IslandPath-DIMOB v1.0.0, while prophages were predicted using Phigaro v2.3.0.	
Data source location	Institution: Academy of National Food and Strategic Reserves Administration
City/Town/Region: Beijing
Country: China	
Data accessibility	Repository name: NCBI (National Center for Biotechnology Information) GenBank
Nucleotide database
Data identification number: BioProject accession number: PRJNA1111238, NCBI SRA Accession Number: SRR29015671 and SRR29015672, GenBank Accession Number: CP156920-CP156924.
Direct URL to data: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1111238
https://www.ncbi.nlm.nih.gov/sra/SRR29015671
https://www.ncbi.nlm.nih.gov/sra/SRR29015672
https://www.ncbi.nlm.nih.gov/nuccore/CP156920:CP156924[accn]	
Related research article	None	

1 Value of the Data

• The genomic data reveals the potential genes responsible for converting ZEN to α-/β-ZEL in L. brevis PYN10_6_2.

• The genomic data serves as a valuable resource for comparative genomic studies of Levilactobacillus spp. and related taxa, and can be used to assess the evolutionary relevance of short-promoted Levilactobacillus species.

• The genomic data can contribute to the evaluation of the biosafety of L. brevis PYN10_6_2 by identifying genetic elements associated with pathogenicity, resistance, and other risk factors.

2 Background

Zearalenone (ZEN), a prevalent mycotoxin produced by Fusarium species, is frequently found contaminating cereals [1]. Microbial degradation as a promising and environmentally friendly approach for ZEN mitigation [2]. One of the known mechanisms of ZEN degradation is the reduction of the ketocarbonyl groups, which reduces ZEN to α-/β-ZEL [3].

In a previous study, we isolated a strain, L. brevis 10_6_2, from T. molitor larval feces, demonstrating its ability to reduce ZEN to α-/β-ZEL [4]. Studies have shown that one of the metabolic pathways of zearalenone in rat liver is the reduction of the ketocarbonyl groups, which is thought to be catalyzed by 3α-/3β-HSD [5]. However, the genetic basis of this biotransformation process and the potential biosafety of the strain remain unexplored. The present study aims to address these knowledge gaps through a comprehensive genomic analysis of L. brevis PYN10_6_2. Specifically, we aim to: (1) identify potential genes responsible for ZEN reduction, elucidating the molecular mechanisms underlying this biotransformation process; (2) confirm the taxonomic classification of the strain through whole-genome sequencing; (3) evaluate its biosafety profile by screening for pathogenicity islands, virulence factors, and antibiotic resistance genes, a preliminary assessment for its potential application in food and feed safety.

3 Data Description

This study presents the whole genome sequencing information of L. brevis PYN10_6_2, which is able to convert ZEN to α-/β-ZEL. Genome sequencing was performed using the PacBio RS II and Illumina HiSeq X Ten sequencing platforms. L. brevis PYN10_6_2 generated 3,276,526 paired-end reads on the Illumina platform with a genome sequencing depth of 368.47×. The L. brevis PYN10_6_2 genome was sequenced to 104.05× depth using the PacBio platform, which generated 32,302 reads. The assembly coverage was calculated to be 99.95× based on the overall bases mapped to the genome. Genome assembly revealed that L. brevis PYN10_6_2 possesses a 2,745,725 bp, consisting of one circular chromosome and four circular plasmids, with a G + C content of 45.52 % (Fig. 1). Raw read length data can be found in the Sequence Read Archive (SRA), where the PacBio RS II and Illumina HiSeq X Ten reads are available under the accession numbers SRR29015671 and SRR29015672, respectively. The assembled genome can be found in Whole Gun Shotgun (WGS) under the accession numbers CP156920-CP156924. The four plasmids found in L. brevis PYN10_6_2 with the accession numbers CP156921, CP156922, CP156923, and CP156924. The genome was annotated using NCBI PGAP yielded a total of 2,660 CDSs including 2,623 CDSs containing coding proteins, 15 rRNAs, 66 tRNAs and 1 CRISPR (Table 1). The ANI between typical Levilactobacillus strains was obtained using OAT based on the complete genome sequence, and the ANI of L. brevis PYN10_6_2 with L. brevis DSM 20054T was 98.94 % (Fig. 2).Fig. 1 Construction of a circular genome map of L. brevis PYN10_6_2 strain using CGView server (https://proksee.ca/ assessed on 26th May 2024). The first and fourth circles of the circle map from outside to inside are CDS, tRNA and rRNA on the positive and negative strands; the third and second circles are GC skew and GC content on the positive and negative strands, respectively; red color represents CDS; the green peaks represent GC-skew+; purple color represents GC-skew-; and the black peaks are G + C content. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Fig 1

Table 1 General genomic features of L. brevis PYN10_6_2.

Table 1Feature	Value	
Total sequence length (bp)	2,745,725 bp	
G+C content (%)	45.52 %	
Plasmids	4	
CDSs (total)	2,660	
CDSs (with protein)	2,623	
Genes assigned to COG	2,053	
Genes assigned to KEGG	1,856	
rRNA	15	
tRNA	66	
Number of CRISPR locus	1	
Number of Prophages	2	

Fig. 2 Heatmap generated by OAT software showing ANI values for L. brevis PYN10_6_2 and related Levilactobacillus reference strains.

Fig 2

When combined with the annotation results, Functional annotation of CDS by gene evolutionary genealogy, COG annotation results showed that a total of 2,053 gene annotations were obtained, accounting for 76.52 % of the total number of genes. In contrast, KEGG annotation results showed 1,856 genes assigned, representing 69.18 % of the total number of genes (Table 1). The pathway of Steroid hormone biosynthesis in the Lipid metabolism hierarchy classification of KEGG and the Lipid transport and metabolism type of COG are known to have three coding sequences capable of producing 3α-/3β-HSD (ABE785_00920, ABE785_01750, ABE785_01835) present in the bacterial chromosome (Fig. 3). This genome annotation data is valuable for elucidating the molecular mechanism of the biotransformation process of L. brevis PYN10_6_2 to ZEN to α-/β-ZEL.Fig. 3 COG and KEGG annotations for the coding sequence of L. brevis PYN10_6_2. A, COG annotation for the coding sequence; B, KEGG annotation for the sequence.

Fig 3

Information on resistance genes contained in each genome was obtained by CARD comparison prediction, and a total of 182 resistance genes appeared on chromosomes containing resistance genes to tetracyclines, fluoroquinolones, macrolides, and other antibiotics. An annotation map of the virulence factor genes was obtained by comparison with the VFDB, and a total of 237 virulence factors were statistically obtained, of which 8.43 % were present in the plasmids, 17 virulence factors were present in plasmid CP156921, and 3 virulence factors were present in plasmids CP156923 and CP156924, respectively. Ten genomic islands were obtained by genomic island prediction using IslandViewer and IslandPath-DIMOB, four of which were present on plasmids. Two prophages present on the chromosome were obtained by Phigaro's prediction. This genome annotation data is valuable for the preliminary assessment of the biosafety of L. brevis PYN10_6_2 in food and feed applications.

4 Experimental Design, Materials and Methods

4.1 Sample collection

L. brevis PYN10_6_2, previously isolated from T. molitor larval feces and shown to convert ZEN into α-/β-ZEL [4], was subjected to whole-genome sequencing. The microbial strain is inoculated into fresh MRS broth and incubated at 37 °C with a shaking speed of 220 revolutions per minute (rpm) for a duration of 9 h to reach the exponential growth phase and collection of bacterial cells. Bacterial cells were harvested and submitted to Shanghai Majorbio Bio-Pharm Technology Co., Ltd. for whole-genome sequencing.

4.2 Genomic DNA preparation

Genomic DNA was extracted using Wizard® Genomic DNA Purification Kit (Promega) according to manufacturer's protocol. The concentration and purity of genomic DNA were assessed by a TBS-380 fluorometer (Turner BioSystems Inc., Sunnyvale, CA). High quality DNA (OD260/280=1.8–2.0, >20 µg) was used for subsequent analyses.

4.3 Whole genome sequencing and assembly

The genome of L. brevis PYN10_6_2 was sequenced using a hybrid approach combining short-read Illumina and long-read PacBio SMRT sequencing technologies. For Illumina sequencing, genomic DNA was fragmented to 400–500 bp using a Covaris M220 Focused ultrasonicator and libraries were prepared with the NEXTflex™ Rapid DNA-Seq Kit. Sequencing was performed on an Illumina HiSeq X Ten platform (2 × 150 bp paired-end reads). For PacBio sequencing, genomic DNA was sheared, purified, end-repaired, and ligated to SMRTbell adapters following manufacturer's recommendations. The resulting sequencing library was purified three times using 0.45× volumes of Agencourt AMPure XP beads (Beckman Coulter Genomics, MA) according to the manufacturer's instructions. A ∼10 kb insert library was prepared and sequenced on a SMRT cell following standard methods. Raw reads from both platforms underwent quality assessment and trimming. Genome assembly was conducted using Unicycler v0.4.8[6], and the resulting assembly was polished with Pilon v1.22[7]. Complete circular chromosomes and plasmids were identified based on overlapping contig ends.

4.4 Whole genome prediction and annotation

Genome annotation was performed using the PGAP v6.7. Functional annotation of CDS was performed through the EggNOG and KEGG databases using the Diamond v0.8.35 Sequence Comparison Tool with an E-value threshold of 10−5 [[8], [9], [10]]. The assembled genome was visualized as a circular map using the CGView web service (https://proksee.ca/) [11]. ANI of strain with closely related species was determined using the Homologous Average Nucleotide Identity Tool v0.93.1 (OAT), resulting in the construction of a phylogenetic tree and a heatmap [12].

Genomic island prediction was performed using IslandViewer 4 and IslandPath-DIMOB v1.0.0 [13]. Prophage prediction was performed using Phigaro v2.3.0 [14]. The virulence factors were obtained from the VFDB using the Diamond with an E-value threshold of 10−5 [15]. Antibiotic resistance genes were obtained by Diamond with an E-value threshold of 10−5 via CARD v3.2.6 [16].

Limitations

None.

Ethics Statement

The authors have read and follow the ethical requirements for publication in Data in Brief and confirming that the current work does not involve human subjects, animal experiments, or any data collected from social media platforms.

CRediT Author Statement

Mengru Zhao: Conceptualization, Data curation, Formal analysis, Methodology, Validation, Writing – original draft; Baoyuan Guo: Project administration, Funding acquisition; Shuzhen Jiang: Supervision; Yang Wang: Conceptualization, Investigation, Supervision, Writing – review and editing, Funding acquisition.

Data Availability

Levilactobacillus brevis strain PYN10-6-2 plasmid pLB1106_4, complete sequence (Original data) (NCBI GenBank Nucleotide database).

Levilactobacillus brevis strain PYN10-6-2 plasmid pUCCLB521_C, complete sequence (Original data) (NCBI GenBank Nucleotide database).

Levilactobacillus brevis strain PYN10-6-2 plasmid pLB1174-1, complete sequence (Original data) (NCBI GenBank Nucleotide database).

Levilactobacillus brevis strain:PYN10_6_2 Raw sequence reads (Original data) (NCBI GenBank Nucleotide database).

Levilactobacillus brevis strain PYN10-6-2 chromosome, complete genome (Original data) (NCBI GenBank Nucleotide database).

Levilactobacillus brevis strain PYN10-6-2 plasmid p1, complete sequence (Original data) (NCBI GenBank Nucleotide database).

SRX24542287: Bridge PCR of Levilactobacillus brevis PYN10-6-2:Tenebrio Molitor feces (Original data) (NCBI GenBank Nucleotide database).

SRX24542286: Bridge PCR of Levilactobacillus brevis PYN10_6_2:Tenebrio Molitor feces (Original data) (NCBI GenBank Nucleotide database).

Acknowledgments

This research was supported by grants from the Key Research and Development Program of Ningxia, China (2022BBF03031) and the 10.13039/501100001809 National Natural Science Foundation of China (31972605 ).

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
==== Refs
References

1 Kong L. Zhao A.-H. Wang Q.-W. Feng Y.-Q. Yan Z.-H. Li M.-H. Zhang F.-L. Wang H. Shen K.-Y. Liu Y. Sun Y.-J. Shen W. Li L. Maternal Zearalenone exposure impacted ovarian follicle formation and development of suckled offspring Sci. Total Environ. 788 2021 147792 10.1016/j.scitotenv.2021.147792
2 Murtaza B. Jin B. Wang L. Li X. Saleemi M.K. Majeed S. Khatoon A. Li G. Xu Y. Mitigation of zearalenone in vitro using probiotic strains LWT–Food Sci. Technol. 186 2023 115265 10.1016/j.lwt.2023.115265
3 Cai P. Liu S. Tu Y. Shan T. Toxicity, biodegradation, and nutritional intervention mechanism of zearalenone Sci. Total Environ. 911 2024 168648 10.1016/j.scitotenv.2023.168648
4 Yuetong P. Degradation and Metabolic Behavior of Zearalenone in Tenebrio Molitor [D] 2022 Beijing: Beijing Union University
5 Olsen M. Pettersson H. Kiessling K.H. Reduction of zearalenone to zearalenol in female rat liver by 3α-hydroxysteroid dehydrogenase Acta Pharmacol. Toxicol. 48 1981 157 161 10.1111/j.1600-0773.1981.tb01602.x
6 Wick R.R. Judd L.M. Gorrie C.L. Holt K.E. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads PLoS Comput. Biol. 13 2017 e1005595 10.1371/journal.pcbi.1005595
7 Walker B.J. Abeel T. Shea T. Priest M. Abouelliel A. Sakthikumar S. Cuomo C.A. Zeng Q. Wortman J. Young S.K. Earl A.M. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement PLoS ONE 9 2014 e112963 10.1371/journal.pone.0112963
8 Jensen L.J. Julien P. Kuhn M. von Mering C. Muller J. Doerks T. Bork P. eggNOG: automated construction and annotation of orthologous groups of genes Nucleic Acids Res. 36 2008 D250 D254 10.1093/nar/gkm796 17942413
9 Kanehisa M. Goto S. KEGG: kyoto encyclopedia of genes and genomes Nucleic Acids Res. 28 2000 27 30 10.1093/nar/28.1.27 10592173
10 Buchfink B. Xie C. Huson D.H. Fast and sensitive protein alignment using DIAMOND Nat. Methods 12 2015 59 60 10.1038/nmeth.3176 25402007
11 Grant J.R. Enns E. Marinier E. Mandal A. Herman E.K. Chen C.Y. Graham M. Van Domselaar G. Stothard P. Proksee: in-depth characterization and visualization of bacterial genomes Nucleic Acids Res. 51 2023 10.1093/nar/gkad326 W484-w492
12 Lee I. Ouk Kim Y. Park S.-C. Chun J. OrthoANI: an improved algorithm and software for calculating average nucleotide identity Int. J. Syst. Evol. Microbiol. 66 2016 1100 1103 10.1099/ijsem.0.000760 26585518
13 Bertelli C. Laird M.R. Williams K.P. Lau B.Y. Hoad G. Winsor G.L. Brinkman F.S.L. IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets Nucleic Acids Res. 45 2017 10.1093/nar/gkx343 W30-w35
14 Fouts D.E. Phage_Finder: automated identification and classification of prophage regions in complete bacterial genome sequences Nucleic Acids Res. 34 2006 5839 5851 10.1093/nar/gkl732 17062630
15 Chen L. Zheng D. Liu B. Yang J. Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on Nucleic Acids Res. 44 2016 D694 D697 10.1093/nar/gkv1239 26578559
16 Jia B. Raphenya A.R. Alcock B. Waglechner N. Guo P. Tsang K.K. Lago B.A. Dave B.M. Pereira S. Sharma A.N. Doshi S. Courtot M. Lo R. Williams L.E. Frye J.G. Elsayegh T. Sardar D. Westman E.L. Pawlowski A.C. Johnson T.A. Brinkman F.S. Wright G.D. McArthur A.G. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database Nucleic Acids Res. 45 2017 10.1093/nar/gkw1004 D566-d573
