Deciphering the biodesulfurization pathway employing marine mangrove Bacillus aryabhattai strain NM1-A2 according to whole genome sequencing and transcriptome analyses

Abstract

In the biogeochemical cycle, sulfur oxidation plays a vital role and is typically referred to as the elemental sulfur or reductive sulfide oxidation process. This study aimed to characterize a subtropical mangrove-isolated bacterial strain using biochemical, whole-genome, and transcriptome sequencing analyses to enhance our understanding of sulfur metabolism and biodegradation from a molecular genetic perspective. Strain NM1-A2 was characterized as Gram-positive and found to have a close molecular phylogenetic relationship with Bacillus aryabhattai. NM1-A2 efficiently converted dibenzothiophene (DBT) into 2-hydroxybiphenyl (2-HBP) via a 4S pathway with 95% efficiency, using enzymes encoded by the dsz operon (dszA, dszB, and dszC), which determine monooxygenases (DszA & DszC) and desulfinase (DszB). The whole-genome sequence of NM1-A2 had a length of approximately 5,257,678 bp and included 16 sulfur metabolism-related genes, featuring the ABC transport system, small subunit (ssu) and cysteine (cys) gene families, and adenosine 5′-phosphosulfate (APS) and 3′-phosphoadenosine-5′-phosphosulfate (PAPS) biosynthesis-related genes. Transcriptomic analysis of NM1-A2 using three sulfur groups—magnesium sulfate (MS), sulfur powder (SP), and sodium thiosulfate (ST) resulted in a significant number of differentially expressed genes (1200, 2304, and 2001, respectively). This analysis revealed that intracellular cysteine concentration directly regulated the expression of cys and ssu genes. Sulfate did not directly affect cys gene expression but repressed ssu gene expression. The cys gene expression levels decreased during the conversion of sulfate to sulfide and cysteine. The transcriptomic data was validated by analyzing the expression patterns of NM1-A2 using real-time quantitative PCR validation analysis. The expression levels of cysl, mccB, and nrnA were significantly upregulated, while cysH, metB, and sat were downregulated in the SP, ST, and MS groups, respectively. This research contributes to our understanding of marine mangrove microorganisms' bacterial efficiency through characterization, whole-genome, and transcriptome sequencing-based molecular degradation of organic compounds in the mangrove ecosystem, which may enhance nutrient availability.