{'fraction_non_alphanumeric': 0.04148964672560952, 'fraction_numerical': 0.025337773184827956, 'mean_word_length': 5.067812354853692, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Insights Into Mechanism of the Naphthalene-Enhanced Biodegradation of Phenanthrene by Pseudomonas sp. SL-6 Based on Omics Analysis

The existence of polycyclic aromatic hydrocarbons (PAHs) in contaminated environment is multifarious. At present, studies of metabolic regulation focus on the degradation process of single PAH. The global metabolic regulatory mechanisms of microorganisms facing coexisting PAHs are poorly understood, which is the major bottleneck for efficient bioremediation of PAHs pollution. Naphthalene (NAP) significantly enhanced the biodegradation of phenanthrene (PHE) by Pseudomonas sp. SL-6. To explore the underlying mechanism, isobaric tags for relative and absolute quantification (iTRAQ) labeled quantitative proteomics was used to characterize the differentially expressed proteins of SL-6 cultured with PHE or NAP + PHE as carbon source. Through joint analysis of proteome and genome, unique proteins were identified and quantified. The up-regulated proteins mainly concentrated in PAH catabolism, Transporters and Electron transfer carriers. In the process, the regulator NahR, activated by salicylate (intermediate of NAP-biodegradation), up-regulates degradation enzymes (NahABCDE and SalABCDEFGH), which enhances the biodegradation of PHE and accumulation of toxic intermediate–1-hydroxy-2-naphthoic acid (1H2Na); 1H2Na stimulates the expression of ABC transporter, which maintains intracellular physiological activity by excreting 1H2Na; the up-regulation of cytochrome C promotes the above process running smoothly. Salicylate works as a trigger that stimulates cell to respond globally. The conjecture was verified at transcriptional and metabolic levels. These new insights contribute to improving the overall understanding of PAHs-biodegradation processes under complex natural conditions, and promoting the application of microbial remediation technology for PAHs pollution.


INTRODUCTION
Polycyclic aromatic hydrocarbons (PAHs) are organic compounds with fused aromatic rings in molecular structure (Alvarez-Ramírez and Ruiz-Morales, 2020). PAHs in the environment are mainly caused by incomplete combustion of coal, oil, wood and macromolecular organic compounds, as well as the exploitation, transportation, use and emission of petroleum and its products . PAHs are one of the main organic pollutants in soil, which have strong mutagenesis, carcinogenesis and teratogenesis effects, thus listed as persistent organic pollutants by the UN Environment Programme (Usman et al., 2016). 16 kinds of PAHs, including naphthalene (NAP) and phenanthrene (PHE), have been identified by written laws or decrees in various countries and regions, the general background value of PAHs in soil is 100-1000 ng/kg, and the lowest published toxic dose (TDLo) of NAP and PHE to mice are 27 and 71 mg/kg, respectively (Fang et al., 2020;Wu et al., 2020;Han et al., 2021). Due to their comparatively simpler structure, NAP and PHE are representative substances in PAHs, and have attracted extensive attentions and researches. It is relatively easier to study the degradation mechanism of NAP and PHE, which is of great significance for bioremediation of PAHs contaminated environment.
In recent years, PAHs degrading bacteria and their degradation mechanism have been intensively studied (Weissenfels et al., 1990;Feitkenhauer et al., 2003), and many genera with degradation ability have been identified (Sàágua et al., 2002;Roy et al., 2011;Lin et al., 2014;Alvarez et al., 2019). However, most of the studies focus on the biodegradation of a single PAH by microorganisms. The actual situation is more complex, many kinds of PAHs coexist in natural environment (Telesiński and Kiepas-Kokot, 2021). When PAHs coexist, the degradation of one PAH could be enhanced or inhibited by others. Sphingomonas yanoikuyae JAR02 (Rentz et al., 2008) grown on salicylate completely removed the solubility limit of 1.2 mg/L benzo(α)pyrene in 20 h. For the efficient pyrenedegrading strain Mycobacterium sp. N12 (Fengchai et al., 2011), it could not grow with benzo(α)pyrene as the sole carbon source and energy source. However, when PHE and pyrene coexisted, the degradation rate of benzo(α)pyrene reached 79% within 9 days. For PHE-degrading strain P-1 (Peng, 2015), the addition of NAP reduced the degradation rate of PHE by 41.7%, anthracene increased the degradation rate by 12.65%, and the removal rate of PHE decreased by 22.82% in the presence of both NAP and anthracene. Only microbial degradation phenomena or characteristics of mixed PAHs have been so far described, and the degradation mechanism has not been reported. The revealing of the underlying mechanism of biodegradation of coexisting PAHs has guiding significance for microbial remediation of PAHs-polluted environment.
Polycyclic aromatic hydrocarbons pollution mostly occurs in extreme saline alkali environment (Zhao et al., 2008). Biodegradation of PAHs is very difficult in saline-alkaline environment due to the inhibition of microbial growth under saline-alkaline stress (Wang et al., 2012). Halotolerant microorganism thus represent promising candidates for the implementation of PAHs bioremediation strategies in hypersaline conditions (González-Abradelo et al., 2019). In this study, a moderately halophilic strain, Pseudomonas sp. SL-6, which has an efficient ability to degrade PAHs, was screened. It is very interesting that the PHE degradability of Pseudomonas sp. SL-6 has been greatly enhanced in the presence of NAP. To understand the mechanism(s) behind the phenomenon, a combined analysis of proteomics and genomics was used.
A preliminary conclusion was reached: Salicylate is the trigger that makes cell respond globally. The addition of NAP was accompanied by the production of salicylate. Salicylateactivated NahR promoted the expression of degradation enzymes, which accelerated the biodegradation of PHE and accumulation of 1-hydroxy-2-naphthoic acid (1H2Na). The upregulation of ABC transporter excretes the substance 1H2Na to maintain intracellular physiological activity. The up-regulation of cytochrome C promotes the above development running smoothly. The conjecture was verified at transcriptional and metabolic levels. These results provide theoretical support for further understanding of microbial degradation characteristics of coexisting PAHs.

Bacterial Preparation
Pseudomonas sp. SL-6 was isolated from oil contaminated soil samples in Shengli Oilfield of Shandong, China. After growing to exponential phase in LB (3% NaCl, w/v), Pseudomonas sp. SL-6 was cultured in minimal salt medium (MSM) with PHE or PHE + NAP as carbon source (10% inoculation, 30 • C, 180 rpm).

Genome Sequencing
Genomic DNA of Pseudomonas sp. SL-6 was extracted using DN12 microbial DNA isolation kit (Aidlab Biotech, Beijing, China) and quality detected by fluorescent dye (Quant-iT PicoGreen dsDNA Assay Kit). The genome sequencing library for Illumina TruSeq Nano DNA LT was constructed in accordance with Illumina TruSeq DNA Sample Preparation Guide. After discarding low quality and short length subreads, a total of 2,093,925,763 bp of clean data were obtained. The clean data were assembled by Canu (version 1.5) (Koren et al., 2017). The sequencing depth was 471.99×, and the total length of the complete genome was 4.43 Mb. The genome of Pseudomonas sp. SL-6 has been deposited at NMDC under the accession number NMDC60018261.

Protein Extraction, Digestion, and Isobaric Tags for Relative and Absolute Quantification-Labeling
In order to identify the differential proteins, according to the degradation characteristics of Pseudomonas sp. SL-6, the 24 h was set as the protein extraction point. The proteins were extracted by total protein extraction kit for gram negative bacteria (BANGFEI BIOSCIENCE, BP0601-50), and its quantitative and SDS-PAGE electrophoresis detection showed that the quality inspection was qualified, the total protein could meet the requirements of downstream experiments. Trypsin (Promega, V5113) was added to 200 µg protein solution at ratio of 1:50, and enzymolysis was carried out at 37 • C for 16 h. 100 µg peptide sample was labeled according to iTRAQ Reagent-8Plex Multiplex Kit (AB Sciex, 4390812) (Unwin et al., 2010;Cox et al., 2014). Then the labeled peptides were pooled together.

Acquisition of Valid Proteomic Data
In order to analyze the mechanism by which NAP enhances the biodegradation of PHE, differential proteomics was performed. The group cultured with PHE as sole carbon was set as control, PHE + NAP as treated group, three parallels for each group. Using p-value ≤ 0.05 and fold change > 1. 2 (i.e., up-regulated > 1.2 and down-regulated < 0.833), 2573 proteins were identified, of which 75 were up-regulated and 101 were down-regulated (Zhou et al., 2016). Through the Accession number of identified differential proteins, the amino acid sequence was obtained from Uniprot 1 , then was blasted with genome sequence of Pseudomonas sp. SL-6 using tblastn 2 of NCBI. The genome was annotated by RAST (Aziz et al., 2008;Weirong et al., 2020) 3 . Through the mutual corroboration of proteome and genome, 63 up-regulated and 80 down-regulated were obtained, which were divided into seven categories due to their function, 33 proteins were discarded (no significant similarity found/function unknown) ( Table 1 and Supplementary Table 2). Based on the existing knowledge, in order to better explain the mechanisms of enhancement, three categories (Table 1) were discussed: PAH catabolism, Transporters, and Electron transfer carriers, in this paper.

Quantitative Real-Time PCR
Primer-BLAST of NCBI 4 was used to design Quantitative Realtime PCR (qPCR) primers for NAP-degradation enzymes genes (nah cluster and sal cluster) and 16S rRNA gene (Supplementary Table 1). The electrophoretic bands of conventional PCR using Pseudomonas sp. SL-6's genomic DNA as template showed that the primers were specific (data not shown). MSM + PHE as the control group, MSM + NAP and MSM + salicylate as treated groups were used to impulse MSM-washed Pseudomonas sp. SL-6 for 8 h that had grown to OD 600 = ∼0.5 in LB. Bacterial RNA Rapid Extraction Kit (GeneBetter, Beijing) was used to extract total RNA. RT Kit With gDNA Eraser (GeneBetter, Beijing) was used to perform reverse transcription. 2 × SYBR Green qPCR Mix Kit (GeneBetter, Beijing) was used for qPCR. qPCR was performed for 15 genes of each group with parallel experiments. 2 − Ct was used for relative quantitative analysis of the data.

Enhanced Effect of Naphthalene on Biodegradation of Phenanthrene
Pseudomonas sp. SL-6, a moderately halophilic bacterium, is able to degrade many PAHs and their derivatives in saline alkali environment. When NAP (100 mg/L) was used as the sole carbon source or PHE (100 mg/L) coexisted, Pseudomonas sp. SL-6 could completely biodegrade NAP in 1 day. When PHE was used as the sole carbon (100 mg/L), more than half of PHE could be degraded in 1 day. However, the addition of NAP greatly enhances the degradation of PHE by Pseudomonas sp. SL-6. The PHE-degradation rate of the group with additional NAP at hour 24 was 24.7% higher than the group without NAP. At hour 48, more than 90% PHE was degraded in the group with additional NAP, while the group with PHE as sole carbon source degraded about 74.4% PHE (Figure 1). Meanwhile, there was no significant difference in biomass between the two group. Many functional microbial populations enter viable but non-culturable (VBNC) state  under unfavorable conditions to deal with environmental stresses, including oligotrophic nutrients and high *Data are from the differential expression proteins analysis. # Data are from the genome annotation analysis by RAST.
Frontiers in Microbiology | www.frontiersin.org concentration of pollutants . Under the conditions of NAP or/and PHE as carbon sources, SL-6 didn't show an obvious growth comparing to growing on LB (Supplementary Figure 1). The reason could be that SL-6 entered VBNC state to cope with the environmental pressure (Fida et al., 2017). As mentioned in Introduction, the phenomenon that one PAH affects the degradation of another has been widely observed. Different characteristics of PHE-biodegradation by Pseudomonas sp. SL-6 between with and without NAP imply great intracellular changes of metabolic regulation.

Metabolic Characteristics of Polycyclic Aromatic Hydrocarbons on Pseudomonas sp. SL-6 Genome
The total length of the genome sequence was 4,436,370 bp with a GC content of 63.8% (Supplementary Figure 2). Through genome annotation (Aziz et al., 2008), 4201 coding sequences were annotated, and Pseudomonas sp. SL-6 harbors 64 PAHs metabolism genes, which constitutes or relates to catechol branch of beta-ketoadipate pathway (eight genes), salicylate and gentisate catabolism (two genes) and so on (Supplementary Figure 3). This indicates that Pseudomonas sp. SL-6 has great degrading potential for PAHs. Pseudomonas sp. SL-6 contains no plasmid and possesses nah and sal cluster (Figure 2). nah and sal cluster encodes enzymes that convert NAP to salicylate and salicylate to pyruvate, respectively. NAH7-like plasmid that contains nah and sal clusters is the most classic NAP-biodegradation plasmid, which was reported to be able to degrades PAHs like anthracene and PHE besides NAP (Ghosal et al., 2016). RHD protein (NahA) was found to easily oxidize NAP and PHE to dihydrodiol form, which had a broad substrate specificity for PAHs with the highest activity to PHE (Wang et al., 2017). The enzymes involved in the conversion of NAP to salicylate can degrade PHE to 1H2Na. 1H2Na is hydroxylated to 1,2-dihydroxynaphthalene, which then enters the NAP-degradation pathway (Yin et al., 2020). In Pseudomonas sp. SL-6, PHE is transferred to 1H2Na by NahABCDE, then 1H2Na is converted to 1,2-dihydroxynaphthalene by salicylate hydroxylase (Balashova et al., 2001;Pathak et al., 2016;Sun et al., 2019), which enters into the salicylate-pathway of NAP biodegradation and finally completely degraded (Figure 3). There are 215 genes annotated as regulator, among which the transcriptional regulator, LysR family (Figure 2, nahR) locates in sal cluster, which may be a regulator for this pathway.

Polycyclic Aromatic Hydrocarbon Catabolism
Polycyclic Aromatic Hydrocarbon catabolism is the most significantly up-regulated category. In PAH catabolism, the proteins are mainly enzymes related to NAP-degrading. There are Salicylate hydroxylase (A0A023WZM3, 2.25 folds) and Naphthalene 1,2-dioxygenase (A0A023WZF5, 2.11 folds) that are two indicative proteins, as well as other degradation enzymes ( Table 1). Locate all the up-regulated proteins from PAH catabolism on genome, they form two clusters with a distance of 10204 bp, nah and sal (Figure 2), which encode integral NAPdegrading salicylate pathway's enzymes. Nah cluster encodes degradation enzymes for ring-open of NAP to salicylate. And sal cluster are genes for degradation enzymes of salicylate to acetaldehyde. Accordingly, the two clusters constitute a complete pathway for PHE catabolism.
The results of genomic annotation and differential expression proteins show that, there is no typical cluster or gene directly related to PHE degradation, such as phn, phd, nid, or nag genes (Singleton et al., 2009). 3,4-dihydroxynaphthalene, 1-hydroxy-2naohthoic acid (1H2Na), salicylate and other metabolites were detected by HPLC when using PHE as sole carbon to culture Pseudomonas sp. SL-6 (data not shown). The results suggested that Pseudomonas sp. SL-6 might utilize the salicylate-pathway of NAP biodegradation to metabolize PHE completely. NahAc and c23o genes (Asmaa et al., 2020) are considered to be the main marker genes for PHE biodegradation by bacteria. The two proteins are significantly up-regulated in proteome data ( Table 1, A0A023WZF5, A0A023WZM3), which are coded by nahA and salA in nah and sal cluster, respectively (Figure 2). Addition of NAP stimulates the up-regulation of degradation enzymes, thus promoting PHE degradation.
The expression of ABC-transporter-related proteins in Pseudomonas putida KT2440 cultured with phenol or pyruvate was much higher than that cultured with glucose, suggesting that toxic stress is the reason for the increased expression of ABC transporters (Leonid, 2006). Some transporter proteins (Vandera et al., 2015), such as Asphe3_22000, was expressed in A. phenanthrenivorans Sphe3 cultured only with PHE or phthalic acid, which shares 88% sequence similarity with a possible 1H2Na transporter in Nocardioides sp. KP7 (Saito et al., 1999). ABC transporter plays an important role for cellular response to PHE (Park et al., 2019). The oxidation products of PAHs are more toxic than themselves, which may react with DNA or other biomolecules (Lundstedt et al., 2007). 1H2Na, for example, is the most easily accumulated oxidation product in biodegradation of PAHs, which is highly toxic to most microorganisms (Knecht et al., 2013), and will completely inhibit the growth of bacteria when the concentration is high (1000 mg/L) (Longmei et al., 2018). In the process of Pseudomonas sp. SL-6 degrading PHE, 1H2Na, which made the medium yellow, accumulated in culture medium (Supplementary Figure 4). The reason could be that 1H2Na was actively excreted through ABC transporter to avoid its accumulation in vivo cells (Figure 3). Compared with PHE alone, the addition of NAP enhanced the biodegradation of PHE, which accelerated the formation of 1H2Na. Up-regulation of ABC transporter was conducive to the efflux of toxic substances. Thus, it was beneficial to maintain the intracellular physiological activity, which in turn was instrumental in biodegrading PHE.

Electron Transfer Carriers
What is noted is that in the electron transfer carriers (Table 1), three up-regulated proteins, biogenesis protein CcsA (A0A172WRB4, 1.55 folds) and Cbb3-type cytochrome c oxidase subunit II (A0A0H3YVF2, 1.57 folds), I (A0A023WTI7, 1.55 folds), are identified as cytochrome C and cytochrome C oxidase, respectively. Similarly, in the process of PHE and benzo[a]pyrene biodegradation by Rhodotorula mucilaginosa EXF-1630, high NADPH-cytochrome c reductase activity was observed (Martínez-Ávila et al., 2021). Cytochrome C is an electron transmitter in biological respiration, which transfers electrons to molecular oxygen, accompanying ATP producing. NAP dioxygenase has been studied and it contains three components, A, B, and C (Ferraro et al., 2017). In the NAP dioxygenase system, component A oxidizes NADH in the presence of cytochrome C, that is, the electron flow direction is NADH-A-C-B, then reductive B combines NAP and oxygen to produce cis-1,2-dihydro-1,2-naphthalenediol. In the process of oxygenating benzene-ring structure, electron transport system is needed to supply electrons. ATP synthesis is coupled in the process of electron transport. ABC transporters are ATP driven pumps (Hara et al., 2010), which obtain energy by binding and hydrolyzing ATP (Khan et al., 2020) and drive the transport of substances across the membrane. When Pseudomonas sp. SL-6 biodegrades PHE, the ATP required for ABC transporter to excrete 1H2Na can be provided by electron transfer. To recap briefly, addition of NAP stimulated the up-regulation of enzymes (NahABCDE and SalABCDEFGH), which promoted the conversion of PHE to 1H2Na. The accumulation of 1H2Na impelled the up-regulation of ABC transporter, which excreted 1H2Na to maintain intracellular physiological activity. Furthermore, the up-regulated cytochrome C guaranteed ATP supply for ABC transporter and the electron demand for enzymes, which ensured the smooth development of the above process (Figure 3). The foregoing analysis shows that, addition of NAP is like pulling one hair, then the whole body is affected, cell

Regulator
No regulatory proteins directly related to PAHs degradation were identified in the up-regulated proteins (data not shown). The transcriptional regulators in PAHs biodegradation pathway has been reported in detail (Tropel and van der Meer, 2004), including LysR-type, IclR family, AraC/XylS family, GntR-type, TetR family, MarR-type, FNR-type, two-component transcription system, XylR/NtrC-type, which exist in different bacteria or plasmids and play roles by following different mechanisms. Therefore, genome of Pseudomonas sp. SL-6 was analyzed. Near sal cluster, there is a transcription regulator of LysR-type (Figure 2), nahR, which locates upstream of salA (salicylate hydroxylase), 156 bp apart, reversely. PAHs biodegradation enzymes gene clusters are usually located in plasmids or genome (Weirong et al., 2020). For Pseudomonas sp. SL-6, it contains no plasmid, the nah cluster and sal cluster are in its genome, and nahR gene is near the sal cluster. NahR is encoded by nahR and activated by salicylate. Activated NahR up-regulates the expression of nah operon and sal operon (Yen and Gunsalus, 1985;Tropel and van der Meer, 2004;Engstrom and Pfleger, 2017) by interacting with promoter sequence of −83 to −45 bp upstream of transcription initiation site (Schell and Wender, 1986). The up-regulation of biodegradation enzymes (nah and sal clusters) in PAH catabolism was observed. The potential reason can be that the increased salicylate that stems from adding NAP activates NahR. Based on the above analysis, one conjecture was deduced that salicylate could be the trigger for the global response of the NAP-enhanced degradation of PHE in Pseudomonas sp. SL-6.

Verification of Salicylate-Enhanced Biodegradation at Transcription Level
In order to validate if the observed up-regulating occurs at the transcriptional level, as well as for further verification of proteomic results, the expression of NAP-degradation enzyme genes was analyzed by qPCR. Compared with the control group (PHE as the sole carbon), the up-regulated trend of genes in nah and sal operon and the up-regulated of proteins in iTRAQ share a high consistency after adding NAP, showing different degrees of up-regulation (Figure 4). For instance, nahAa (11.32 folds) and salA (8.95 folds) of qPCR are most obviously up-regulated, corresponding to the most up-regulated Naphthalene 1,2-dioxygenase (A0A023WZF5, 2.11 folds) and Salicylate hydroxylase (A0A023WZM3, 2.25 folds) in iTRAQ data (Table 1), respectively. The two genes are the initial genes in nah and sal cluster, respectively (Figure 2). Most of other genes are up-regulated 2∼6 folds. In the presence of NAP, the up-regulated fold of RNA and protein of these genes is not strictly corresponding. The possible reason can be that different regulation at translation and post-translation levels occurs in different proteins. Comparing with NAP-added group, the upregulation of salicylate-added group is very similar and more obvious. Especially for nahAb, nahAc, nahD, salA, salC, and salD in salicylate-added group is 1.85 times, 2.38 times, 1.60 times, 1.28 times, 2.23 times, and 1.27 times of NAP-added group (Figure 4), respectively. Addition of salicylate does cause the upregulation of nah and sal cluster, even stronger than addition of NAP. This result confirms the conjecture at transcriptional level.

Validation of Salicylate-Enhanced Biodegradation at Metabolic Level
To validate it at metabolic level, salicylate + PHE and NAP + PHE were set as treated groups, PHE as control group, to determine the biodegradation of PHE by Pseudomonas sp. SL-6 ( Figure 5). Both of additions of NAP or salicylate enhance the biodegradation of PHE. The group with additional salicylate, Pseudomonas sp. SL-6 degraded 81.5% PHE at hour 20, about 5% more than that with additional NAP. Similarly, when adding salicylate, Pseudomonas sp. SL-6 degraded more than 95% of the PHE within 2 days. The degradation character of PHE + salicylate resembles that of NAP + PHE group, even slightly better. Whether NAP or salicylate was added, the enhancement both appeared at hour 8, which indicated the mechanism of enhancement could be the same. From hour 16 on, the enhancement of salicylate became stronger than NAP, indicating that salicylate might enhance PHE-biodegradation more directly. The results verify the conjecture above at metabolic level.

CONCLUSION
In the process of the NAP-enhanced degradation of PHE, a global metabolic regulation occurs in Pseudomonas sp. SL-6 ( Figure 3): NahR, activated by salicylate (intermediate of NAP-biodegradation), promotes the expression of degradation enzymes (NahABCDE and Sal ABCDEFGH), which enhances the biodegradation of PHE and the accumulation of toxic intermediate-1H2Na; 1H2Na stimulates the expression of ABC transporter, which maintains intracellular physiological activity by excreting 1H2Na; the up-regulation of cytochrome C promotes the above process running smoothly. Salicylate works as a trigger that stimulates cell to respond globally. Metabolism of cell is a complex and global physiological process. Mining omics data helps us understand the mechanism behind the phenomenon more comprehensively, especially for the factors that are not directly related. The existence of PAHs in natural environment is usually diverse and complicated. The present results provide new insights into the overall understanding of the mechanism of PAHs degradation by microorganism and promote the application of microbial remediation technology for PAHs pollution.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
HC: conceptualization, methodology, formal analysis, investigation, data curation, and writing-original draft. XZ and JL: conceptualization, validation, resources, and investigation. SW: conceptualization, visualization, and investigation. DH: writing-review and editing and visualization. BZ: writingreviewing and editing. HW: writing-reviewing and editing, supervision, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.