Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 30 May 2024
Sec. Microbiotechnology

Study on biodegradation mechanism of Fusarium solani NK-NH1 on the hull wood of the Nanhai No. 1 shipwreck

Yu Wang,&#x;Yu Wang1,2Yeqing Han&#x;Yeqing Han3Naisheng Li&#x;Naisheng Li4Cen WangCen Wang3Kaixuan MaKaixuan Ma3Xinduo HuangXinduo Huang3Jing DuJing Du4Hong Guo,
Hong Guo1,2*Jiao Pan,
Jiao Pan1,2*
  • 1Key Laboratory of Archaeomaterials and Conservation, Ministry of Education, University of Science and Technology Beijing, Beijing, China
  • 2Institute for Cultural Heritage and History of Science and Technology, University of Science and Technology Beijing, Beijing, China
  • 3Department of Microbiology, College of Life Sciences, Nankai University, Tianjin, China
  • 4National Centre for Archaeology, Beijing, China

The Nanhai No. 1 shipwreck is an ancient wooden ship in the Southern Song Dynasty. Currently, serious challenges of microbial diseases exist on the hull wood. This study aimed to obtain microbial samples from the ship hull in December 2021 and analyze the microbial diseases through scanning electron microscopy and high-throughput sequencing to preserve the Nanhai No. 1 shipwreck. The biodegradation mechanism of diseased microorganisms was explored through whole genome sequencing and the detection of enzyme activity and gene expression levels of diseased microorganisms under different conditions. The results showed that there was obvious fungal colonization on the surface of the hull wood and Fusarium solani NK-NH1 was the dominant disease fungus on the surface. NK-NH1 has strong cellulose and lignin degradation ability. Its whole genome size is 52,389,955 bp, and it contains 17,402 genes. It has a variety of key enzyme genes involved in cellulose and lignin degradation. The NK-NH1 dominant degrading enzyme lignin peroxidase has the highest enzyme activity at pH = 4, NaCl concentration of 30%, and FeSO4 concentration of 50 mg/L, while laccase has the highest enzyme activity at pH = 4, NaCl concentration of 10%, and FeSO4 concentration of 100 mg/L. The above research results prove that NK-NH1 is a key fungus to the biodegradation of ship hull wood when it is exposed to air, low pH, high salt, and rich in sulfur iron compounds. This study provides a theoretical basis for the preservation of the Nanhai No. 1 shipwreck.

1 Introduction

Cultural relics’ biological corrosion is caused by biological activities, mainly including microorganisms such as bacteria, fungi, and algae, as well as insects such as moths and beetles. The rapid growth of organisms and secretion of secondary metabolites are the most important reasons for the destruction of cultural relics, which may cause irreversible damage. Among them, microbial corrosion, mainly caused by bacteria and fungi, is a recent research focus in cultural relic preservation. Cultural relics of various materials may be corroded by microorganisms (Mazzoli et al., 2018), including organic cultural relics such as wooden relics (Liu et al., 2017), textiles (Gutarowska et al., 2017), silk (Wang et al., 2023), and leather (Liu et al., 2018a). It also includes inorganic cultural relics, such as stone relics (Li et al., 2018), earth ruins (Sun et al., 2020), ancient buildings (Piñar et al., 2009; Gaylarde et al., 2011), and murals (Ma et al., 2020). In particular, organic cultural relics, whose material composition is natural and complex, can provide rich nutrients for microorganism’s growth (Kip and Veen, 2015; Liu et al., 2017). Factors such as airflow and human activities can lead to the deposition and colonization of a large number of microorganisms on the surface of cultural relics, causing changes in the physical and chemical properties of cultural relics, ultimately leading to their damage (Singh, 2012). Bacteria can secrete a variety of secondary metabolites to corrode cultural relics. The hyphae and spores of fungi and actinomycetes will spread and grow on the surface of cultural relics, affecting their artistic and aesthetic value. Simultaneously, they will also invade the interiors of the material, causing damage to the structure of the cultural relic material (Diaz-Herraiz et al., 2013; Mazzoli et al., 2018).

Marine effluent wooden cultural relics are a typical organic matter cultural relic, mainly carried by wooden sunken ships, which are important cultural relics with large volumes, severe corrosion, and difficult preservation. Owing to long-term immersion in seawater, they all have the following characteristics: (1) The moisture content of wood is extremely high; (2) The wood has a high degree of decay, serious corrosion and degradation, and the content of cellulose and hemicellulose in wood is far lower than that of fresh wood; (3) Wood contains high salt content, including soluble salts dominated by sodium chloride and insoluble salts dominated by sulfur iron compounds. High concentrations of salt will promote the degradation and fracture of cellulose in crystallization and dissolution; (4) Ships are often equipped with iron tools such as nails and pans. Sulfate-reducing bacteria and iron tools on the sea floor will form sulfur iron compounds and penetrate wood under long-term action. Sulfur iron compounds will oxidize to generate sulfuric acid and sulfate under the action of oxygen, leading to wood acidification (Fors, 2008; Zhang et al., 2014). Therefore, it is extremely difficult to preserve wooden artifacts in marine water. Famous wooden artifacts from marine water include the Swedish “Vasa” sunken ship (Hocker, 2006), the British “Mary Rose” sunken ship (Preston et al., 2014), the Australian “Batavia” sunken ship, the Chinese “Xiaobaijiao I,” and the Chinese “Huaguangjiao I,” among others. The salvage methods of these sunken ships and the conservation status are different. However, after the sunken ships were salvaged and brought out of the seawater, microbial corrosion needed to be highly valued.

The Nanhai No. 1 shipwreck was an ancient wooden ship that was wrecked and sunk during the Southern Song Dynasty (1127–1279 AD) when it was transporting porcelain out of the Maritime Silk Road. It is a shining pearl on the Maritime Silk Road, key evidence of the ancient Maritime Silk Road, and has critical significance and value for studying marine wood relics (Ma et al., 2022). Among all the cultural relics of the Nanhai No. 1 shipwreck, the wooden hull relic is the most precious and difficult to protect. This is owing to the huge volume and uneven degree of decay of wooden hulls, the long conservation period, and continuous excavation. The average moisture content of the Nanhai No. 1 shipwreck hull wood is 300–700%, which belongs to moderate and severe corrosion. The pH value of the hull wood is stable at 5–7, indicating a neutral to slightly acidic nature. The hull wood contains 6.07–46.13% cellulose, 0.95–8.08% hemicellulose, and 2.63–6.55% ash. Compared with 75% holocellulose content (including cellulose and hemicellulose) and less than 1% ash content in fresh pine (Passialis, 1997), the holocellulose content in wood is far lower than that in fresh wood, indicating that the wood has a high degree of decay and high salt content in wood.

Our research group has been involved in the scientific and technological preservation project of the Nanhai No. 1 shipwreck since April 2015, mainly researching the microbial diseases of the Nanhai No. 1 shipwreck. The preliminary research results show that the dominant disease fungus on the ship is Fusarium solani NK-NH1, which has a strong ability to degrade wood (Liu et al., 2018b; Han et al., 2021). Based on previous research, this study comprehensively analyzes the microbial diseases of the Nanhai No. 1 shipwreck by combining traditional microbial separation and detection techniques with cutting-edge sequencing techniques. Simultaneously, the entire genome of the dominant disease fungus F. solani NK-NH1 was sequenced, with a focus on studying its biodegradation mechanism. This study has the characteristic of interdisciplinary integration, integrating biology, chemistry, cultural relic protection and other disciplines, reflecting the research characteristics of common orientation and cross-integration. This study provides certain guiding suggestions for the future scientific and technological conservation work of the Nanhai No. 1 shipwreck. Simultaneously, it also provides a reference for the study of microbial diseases in cultural relics and a theoretical basis and data support for the long-term conservation of wooden cultural relics.

2 Materials and methods

2.1 Investigation and sample collection of microbial diseases on the Nanhai No. 1 shipwreck

This study investigated the microbial diseases of the Nanhai No. 1 shipwreck in December 2021. White microbial disease plaques were detected on some areas of the hull surface, mainly in the blind area of the antibacterial agent spraying system. In December 2021, four sampling sites (NH.SH10–NH.SH13) on the hull surface were collected. White plaques were on the surface of these four sampling sites (Figure 1). The sampling sites were all on the deck of Nanhai No. 1 Shipwreck, which was placed near the ship. The spraying system could only spray on some areas of the hull, while the sampling sites could not be sprayed. The average annual temperature of the preservation environment was 25.6°C and average annual humidity was 84.1%.

Figure 1
www.frontiersin.org

Figure 1. Sampling of microbial diseases.

2.2 Scanning electron microscope (SEM) observation

Carbon conductive adhesive was used to adhere microbial samples to the surface of the Nanhai No. 1 shipwreck, which were then dried in a drying dish. After drying, the sample was affixed to the SEM sample stage. Gold was sprayed with a 24-mA current for 300 s, the sample was observed using SEM, and images were recorded. The measurement conditions were as follows: EHT: 15.0 kV, WD: 9.6–10.2 mm, Mag: 0.3KX-5KX.

2.3 High throughput sequencing analysis

The DNeasy PowerSoil Pro Kit (QIAGEN, Germany, Cat. No. 47014) was used to extract total DNA from solid samples. The total extracted DNA from samples collected was sent to NovoMagic Technology Co., Ltd. The ITS rRNA genes of distinct regions (ITS1-5F) were amplified using specific primers (ITS5-1737F, ITS2-2043R) with barcodes. Sequencing libraries were prepared with the TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, United States, Cat. No. 15032317) following the manufacturer’s instructions, and index codes were added. Library quality was assessed using a Qubit@ 2.0 Fluorometer (Thermo Scientific) and an Agilent Bioanalyzer 2,100 system. The library was then sequenced on an Illumina NovaSeq platform, generating 250 bp paired-end reads. Paired-end reads were matched to samples based on their unique barcodes and then trimmed to remove the barcodes and primer sequences. The paired-end reads were merged using FLASH (V1.2.7, http://ccb.jhu.edu/software/FLASH). Quality filtering of raw tags was done under specific conditions to obtain high-quality clean tags following the QIIME (V1.9.1, http://qiime.org/scripts/split_libraries_fastq.html) quality control process. The tags were compared to the Silva database using the UCHIME algorithm1 to identify and remove chimeric sequences. The effective tags were then obtained. Sequence analysis was conducted using the Uparse software (Uparse v7.0.1001, http://drive5.com/uparse/). Sequences with ≥97% similarity were assigned to the same operational OTUs. A representative sequence from each OTU was selected for further annotation, using the Silva Database2 and the Mothur algorithm to annotate the taxonomic information. The abundance information of OTUs was normalized to a standard sequence number corresponding to the sample with the lowest number of sequences. The raw sequencing data are accessible at the NCBI Sequence Read Archive (SRA) with the study accession number PRJNA1015580.

2.4 Fungal whole genome sequencing analysis

NK-NH1 was cultured in a PD liquid medium at 28°C for 3 days. It was centrifuged at 4°C using a 50 mL centrifuge tube at low speed, rinsed with distilled water three times, and the fungal body was collected through centrifugation. It was frozen in liquid nitrogen and sent to BGI Gene (China). The data volume of fungal whole genome sequencing is large and the parsing is complex. Currently, a combination of short-read long and single molecule sequencing is used to achieve de novo assembly of fungal genome fine maps. The raw sequencing data are accessible at the NCBI WGS with the study accession number PRJNA1016167.

2.5 Detection of the activity and specific activity of cellulose and lignin-degrading enzymes in microorganisms

The improved Lowry protein quantification kit was used to measure the total protein content in the supernatant, and then measure the activities of the corresponding enzymes according to the following method.

2.5.1 Detection of cellulase and specific activities

2.5.1.1 Preparation of glucose standard curve

Multiple gradient glucose solutions were prepared. Approximately 1 mL of 1% CM-Na buffer solution was taken, soaked in 50°C water for 30 min, 2 mL of DNS was added, followed by 0.5 mL of glucose gradient solution. It was boiled for 10 min, cooled, and the absorbance at 550 nm was measured. A standard curve was drawn with glucose concentration as the horizontal axis and absorbance value as the vertical axis.

2.5.1.2 Determination of glucose concentration in the supernatant

The experimental group consisted of taking 0.5 mL of the supernatant, adding 1 mL of 1% CMC-Na buffer, bathing at 50°C for 30 min, adding 2 mL of DNS, and boiling for 10 min. The blank control was prepared by taking 1 mL of 1% CMC-Na buffer solution, water bathing at 50°C for 30 min, adding 2 mL of DNS, adding 0.5 mL of supernatant, and boiling for 10 min. The absorbance at 550 nm was measured after cooling in a cold-water bath. Finally, the glucose concentration in the supernatant was obtained based on the glucose standard curve.

2.5.1.3 Determination of cellulase and specific activities

The unit of cellulase activity (U) is defined as the amount of enzyme required to produce 1 mmol of glucose per minute. The enzyme activity unit U/mL per milliliter of enzyme solution = (X × V1)/(M × T × V2); the specific activity of cellulase U/ng = (X × V1)/(M × T × V2 × P). In the above equation, X is the glucose concentration (mg/mL) calculated based on the standard curve; V1 is the volume of the reaction solution (mL); M is the molar mass of glucose; T is the reaction time (min); V2 is the volume of the supernatant (mL); and P is the total protein concentration (ng/mL) calculated based on the standard curve.

2.5.2 LiP activity and specific activity detection

Approximately 1 mL of 125 mmol/L tartaric acid sodium buffer solution was taken, 0.5 mL of 0.16 mmol/L aniline blue solution, 0.5 mL of supernatant, and 0.5 mL of 2 mmol/L H2O2 solutions were added to start the reaction. The absorbance change value at 651 nm was determined in the first 3 min. The LiP activity unit (U) is defined as the amount of enzyme required to decrease the OD value by 0.1 per min. The enzyme activity unit per milliliter of enzyme solution (U/mL) = N/(0.1 × T × V2). Specific activity of LiP (U/ng) = N/(0.1 × T × V2 × P). In the above equation, N is the absorbance change value at 651 nm in the first 3 min; T is the reaction time (min); V2 is the volume of the supernatant (mL); and P is the total protein concentration (ng/mL) calculated based on the standard curve.

2.5.3 MnP activity and specific activity detection

Approximately 3.4 mL of 50 mmol/L sodium lactate buffer and 0.1 mL of 1.6 mmol/L MnSO4 solutions were taken, 0.4 mL of supernatant was added, and preheated at 37°C for 10 min. Following this, 0.1 mL of 1.6 mmol/L H2O2 solution was added and the absorbance change at 240 nm was measured in the first 3 min. The MnP activity unit (U) is defined as the amount of enzyme required to increase the OD value by 0.1 per min. The enzyme activity unit per milliliter of enzyme solution (IU/mL) = N/(0.1 × T × V2). Specific activity of MnP (U/ng) = N/(0.1 × T × V2× P). In the above equation, N is the absorbance change value at 240 nm in the first 3 min; T is the reaction time (min); V2 is the volume of the supernatant (mL); and P is the total protein concentration (ng/mL) calculated based on the standard curve.

2.5.4 Lac activity and specific activity detection

Approximately 3 mL of 200 mmol/L acetic acid buffer was taken, 0.5 mL of 7 mmol/L ABTS solution and 0.5 mL of supernatant were added, and the absorbance change at 420 nm was measured in the first 3 min. The Lac activity unit (U) is defined as the amount of enzyme required to increase the OD value by 0.01 per min. The enzyme activity units contained in each milliliter of enzyme solution (IU/mL) = N/(0.01 × T × V2). Specific activity of Lac (U/ng) = N/(0.01 × T × V2 × P). In the above equation, N is the absorbance change value at 420 nm in the first 3 min; T is the reaction time (min); V2 is the volume of the supernatant (mL); and P is the total protein concentration (ng/mL) calculated based on the standard curve.

2.6 Detection of gene expression levels of dominant degrading enzymes under different culture conditions using qPCR

2.6.1 Setting of different cultivation conditions

Based on the hull wood powder liquid culture medium, four single variables were set successively:

2.6.1.1 The pH value

A pH gradient was set to detect the activity and specific activity of dominant enzymes and expression level of enzyme genes, to obtain the optimal pH value.

2.6.1.2 NaCl concentration

A NaCl concentration gradient was set to detect the activity and specific activity of dominant enzymes, as well as the expression of enzyme genes, under the optimal pH conditions to obtain the optimal NaCl concentration.

2.6.1.3 FeSO4 concentration

The FeSO4 concentration gradient was set to detect the activity, specific activity, and expression level of dominant degrading enzymes under the optimal pH value and NaCl concentration conditions, respectively, to obtain the optimal FeSO4 concentration.

2.6.1.4 PEG400 concentration

The PEG400 concentration gradient was set to detect the activity, specific activity, and expression level of dominant degrading enzymes under the optimal pH value, NaCl concentration, and FeSO4 concentration conditions, respectively, to obtain the optimal PEG400 concentration.

2.6.2 qPCR primer design

The data obtained from whole gene sequencing was analyzed, and combined with relevant literature reports, and Primer 5 software was used to design qPCR primers.

2.6.3 Trizol method for extracting fungal RNA

The culture precipitate was collected by centrifugation at 4°C, ground with liquid nitrogen, and then an appropriate amount of dry powder was taken into a 1.5 mL EP tube. Approximately 1 mL of Trizol was added to the EP tube and shaken for 2 min. Approximately 200 μL chloroform was added, left to stand for 5 min, centrifuged, and the supernatant was aspirated into a new EP tube. Approximately 500 μL isopropanol was added, left to stand at −20°C for 10 min, centrifuged, and the supernatant was discarded. Approximately 500 μL 70% ethanol was added, centrifuged, the supernatant was discarded, and dried for 5 min. For the preparation of the digestive DNA reaction system (RNA, 10× DNase I buffer, recombinant DNase I, RNase inhibitor, RNase-free H2O), a total of 50 μL reacts in a 37°C incubator for 20–30 min. Approximately 50 μL RNase-free H2O and 100 μL phenol/chloroform/Isoamyl alcohol (25:24:1) were added and completely mixed. After centrifugation, the upper layer was taken to a new EP tube and same amount of chloroform/isoamyl alcohol (24:1) was added. After centrifugation, the upper layer was transferred to a new EP tube. NaAC and anhydrous ethanol were added and left at −80°C for 20 min. The supernatant was discarded after centrifugation. Approximately 500 μL 75% ethanol was added, centrifuged, and the supernatant was discarded. After drying, it was dissolved with an appropriate amount of RNase-free H2O for RNA reverse transcription.

2.6.4 RNA reverse transcription to obtain cDNA

The template RNA/primer system was prepared in a PCR tube [template RNA 2 μg, Oligo (dT) 2 μL]. RNase-free H2O replenishment was performed to 12 μL. After 10 min at 70°C, it was quickly placed on ice for 3 min. A reaction system was prepared in the PCR tube (system mentioned above 12 μL, buffer 4 μL, dNTP mixture 1 μL, RNase inhibitor 0.5 μL. RTase 1 μL, and RNase-free H2O 1.5 μL) 20 μL in total. It was held at 42°C for 1 h, cooled at 70°C for 15 min on ice, and the obtained cDNA was used for qPCR detection of enzyme gene expression.

2.6.5 qPCR detection of enzyme gene expression levels

The qPCR reaction system comprised cDNA 1 μL, primers each 0.5 μL, SYBR 10 μL, and RNase-free H2O 8 μL. The prepared reaction system was placed on ice and stored away from light. The following program was set on the qPCR instrument for gene amplification. Finally, the Ct values during the amplification process were obtained and expression levels of different groups were compared. Three parallel samples were prepared for each sample.

2.7 Statistical analysis

This research was carried out in triplicate with completely randomized experimental designs. Data obtained in this investigation were depicted as averages ± standard deviation. The statistical significance between the results of different sample groups were analyzed by T-test. Statistical significance was set at p < 0.05, *p < 0.05, **p < 0.01, ***p < 0.001.

3 Results

3.1 SEM observation and high-throughput sequencing analysis of microorganisms on the hull of the Nanhai No. 1 shipwreck in December 2021

In December 2021, microbial disease samples were collected from four sites (NH.SH10–NH.SH13) on the hull surface. The SEM results indicated that the typical structure of filamentous fungi still adheres to the surface of the ship, with a high abundance (Figure 2A). The high-throughput sequencing results of fungi show that at the phylum level, the dominant phylum was Ascomycota, with an average relative abundance of 92.48% (Figure 2B). At the level of fungi, Fusarium accounted for the largest proportion, with an average relative abundance of 54.23%. Additionally, Harposporium, Coniochaeta, Talaromyces, Acremonium, and Scedosporium were present (Figure 2C). Compared with previous results (Liu et al., 2018b; Han et al., 2021), the dominant disease fungus on the hull surface was Fusarium since 2015.

Figure 2
www.frontiersin.org

Figure 2. Composition analysis of fungal community on the hull of the Nanhai No. 1 shipwreck (December 2021). (A) SEM observation results; (B) relative abundance of fungi at the phylum level; (C) relative abundance of fungi at the genus level.

3.2 Whole genome sequencing analysis of F. solani NK-NH1

3.2.1 Genomic component analysis

The whole genome size of F. solani NK-NH1 is 52,389,955 bp and contains 17,402 genes. The total length of the CDS sequence is 25,625,398 bp, accounting for 48.91% of the total genome length. Non-coding RNAs (including tRNA, rRNA, sRNA, snRNA, and miRNA) collectively account for 0.1534% of the genome. Repetitive sequences (including DNA transposable elements, tandem repeats, and transposable elements) account for 7.2366% of the genome. According to the genome sequence of NK-NH1, Circos software is used to analyze the genome, gene density, non-coding RNA, repeat sequence, GC-content, GC skew, and other information, and draw the whole genome circle of NK-NH1 (Figure 3).

Figure 3
www.frontiersin.org

Figure 3. Whole genome circle graph of F. solani NK-NH1. The circle graph from outside to inside is as follows: (1) genome (sorted by length); (2) gene density (number of genes in 50 kp non-overlapping window); (3) non-coding RNA density (number of non-coding RNA in 100 kp non-overlapping window); (4) repeat sequence (repeat sequence in 50 kp non-overlapping window); (5) GC (GC rate in 20 kp non-overlapping window); (6) GC skew (GC offset in 20 kp non-overlapping window).

3.2.2 Gene function annotation analysis

Database comparison and annotation were conducted on the entire genome of F. solani NK-NH1. A total of 16 databases were annotated, with a total of 16,950 (97.4%) genes annotated. The statistical results are shown in Supplementary Table S1. An analysis was conducted on the annotation of CAZy, eggNOG, and KEGG databases, and the results showed that carbohydrate metabolism accounted for a large proportion of 12.32% in the secondary classification of KEGG databases (Figure 4). The proportion of carbohydrate transport and metabolism annotations obtained in the functional annotations of the eggNOG database is the highest, reaching 48.47% (Figure 5). Among all the gene functional annotations, the genes related to the cellulose and lignin degradation ability of NK-NH1 were explored. After extensive data screening, they were mainly found to be annotated in CWDE, KEGG, NR, eggNOG, GO, Swiss-Prot, KOG, and CAZy databases. The key genes annotated are shown in Table 1.

Figure 4
www.frontiersin.org

Figure 4. Functional annotation distribution of secondary classification of KEGG database in the whole genome analysis of F. solani NK-NH1.

Figure 5
www.frontiersin.org

Figure 5. Functional annotation distribution of eggNOG database in the whole genome analysis of F. solani NK-NH1.

Table 1
www.frontiersin.org

Table 1. Key genes of cellulose and lignin degradation in the whole genome sequencing of F. solani NK-NH1.

3.3 Effects of different culture conditions on enzyme activity of F. solani NK-NH1

Under the cultivation conditions of using the experimental wood powder of the Nanhai No. 1 shipwreck as the sole carbon source, the dominant degrading enzymes of NK-NH1 are LiP and Lac, with the highest enzyme activity occurring for 8.5 days (Supplementary Table S2). Currently, the Nanhai No. 1 shipwreck is in an acidic, high salt, and sulfur-rich iron compounds preservation environment, and archeological workers are using the protective material PEG400 to dewater and strengthen it. By adjusting pH value, and NaCl, FeSO4, and PEG400 concentrations, the effects of various factors on the LiP and Lac activity of NK-NH1 were studied, and the reason for NK-NH1 becoming a dominant fungus was explored. The research results found that the highest enzyme activity of LiP of NK-NH1 was generated under the following conditions: pH = 4, NaCl concentration of 30%, and FeSO4 concentration of 50 mg/L, without the addition of PEG400 (Figure 6). The maximum enzyme activity of Lac is generated under the following conditions: pH = 4, NaCl concentration of 10%, and FeSO4 concentration of 100 mg/L, without the addition of PEG400 (Figure 7). The detection results of enzyme gene expression are consistent with those of enzyme activity detection.

Figure 6
www.frontiersin.org

Figure 6. The variation trend of LiP enzyme activity, specific activity, and enzyme gene expression of F. solani NK-NH1 under different culture conditions. (A) Effects of pH on LiP enzyme and specific activities. (B) Effect of pH on LiP enzyme gene expression. (C) Effect of NaCl concentration on LiP enzyme and specific activities. (D) Effect of NaCl concentration on LiP enzyme gene expression. (E) Effect of FeSO4 concentration on LiP enzyme and specific activities. (F) Effect of FeSO4 concentration on LiP enzyme gene expression. (G) Effect of PEG400 concentration on LiP enzyme and specific activities. (H) Effect of PEG400 concentration on LiP enzyme gene expression. Statistical significance was set at p < 0.05, *p < 0.05,**p < 0.01, ***p < 0.001.

Figure 7
www.frontiersin.org

Figure 7. The variation trend of Lac enzyme activity, specific activity, and enzyme gene expression of F. solani NK-NH1 under different culture conditions. (A) Effect of pH on Lac enzyme and specific activities. (B) Effect of pH on Lac enzyme gene expression. (C) Effect of NaCl concentration on Lac enzyme and specific activities. (D) Effect of NaCl concentration on Lac enzyme gene expression. (E) Effect of FeSO4 concentration on Lac enzyme and specific activities. (F) Effect of FeSO4 concentration on Lac enzyme gene expression. (G) Effect of PEG400 concentration on Lac enzyme and specific activities. (H) Effect of PEG400 concentration on Lac enzyme gene expression. Statistical significance was set at p < 0.05, *p < 0.05,**p < 0.01, ***p < 0.001.

4 Discussion

The Nanhai No. 1 shipwreck, as one of the rare large and precious wooden cultural relics in China, has extremely important historical significance and conservation value. The microbial diseases on the hull are the biggest challenge in its conservation. Although the effective antibacterial agent K100 has been regularly and quantitatively sprayed on-site through the installed spraying system, owing to the complex structure of the ship, there are blind spots in the spraying system, and the loss of antibacterial agents exists, resulting in microbial contamination on some ship surfaces.

Fungi are the most effective and abundant wood decomposers known to date. Their developed hyphae give them the advantage of spreading growth over most prokaryotes. Many fungi isolated from wooden artifacts can produce cellulase and lignin-degrading enzymes. For example, Penicillium and Cladosporium (Zhang et al., 2019) separated from the canoe of the Tang Dynasty in the National Marine Museum of China, and Aspergillus and Penicillium (Osman et al., 2015) separated from the Islamic Art Museum and other cultural relics. Fusarium can be widely distributed in various environments and can degrade wood, significantly damaging wood (Rodriguez et al., 1996; Panagiotou et al., 2003; M’Barek et al., 2019). There is limited research on the mechanism of wood degradation by Fusarium; however, F. solani is currently known to degrade cellulose and lignin, producing cellulase, Lac, LiP, and MnP (Saparrat et al., 2000; Wu et al., 2010; Obruca et al., 2012). After being infected with F. solani, the cellulose and hemicellulose content in wood will decrease (Rahmawati et al., 2020). Benaddou et al. (2024) experiment revealed that F. oxysporum and F. solani were the most selective lignin-degrading fungi among the tested fungi. The cell wall degrading enzyme (CWDE) in Fusarium is the main pathogenic factor for its invasion and infection of plants, mainly including pectinase, cellulase, and protease (Erazo et al., 2021). Among numerous CWDE, the deletion of a single gene does not significantly affect the virulence of Fusarium; however, if two or more coding genes are knocked out simultaneously and an upstream gene regulating CWDE is knocked out, the pathogenicity will be weakened (Ruiz et al., 2016; Paccanaro et al., 2017). F. solani can decompose aromatic acids (Norris, 1980), and the CWDE produced during spore germination can enhance its infectivity (Köller et al., 1982). Currently, F. solani NK-NH1 has colonized the Nanhai No. 1 shipwreck surface and has become a dominant fungus. This study sequenced its entire genome, and through database annotation, it was found that genes related to carbohydrate transport and metabolism account for a large proportion in various major databases, and multiple key genes involved in cellulose and lignin degradation have been compared in various databases.

The hull wood of the Nanhai No. 1 shipwreck has the characteristics of weak acidity, high salinity, and high sulfur iron compounds, and the surface has the protective material PEG. pH is a critical element that affects enzymatic reactions because it regulates the ionization state of acidic or basic amino acids on the enzyme site and/or affinity to substrate (Xu, 1997). It was reported that NaCl at relatively low concentrations can bind active sites to enhance the affinity of enzyme to substrate and facilitate oxygen reduction (Si et al., 2021a). Additionally, some enzymes require metallic ions to retain the conformation of their active sites (Si et al., 2021b). Therefore, we tested the impact of these factors on NK-NH1 and found that under the cultivation condition of experimental shipwreck wood powder as the only carbon source, the highest enzyme activity of the dominant degrading enzyme LiP was produced under the following conditions: pH = 4, NaCl concentration 30%, FeSO4 concentration 50 mg/L, and no PEG400 added. The highest enzyme activity conditions for the dominant degrading enzyme Lac are pH = 4, NaCl concentration of 10%, FeSO4 concentration of 100 mg/L, and no PEG400 added. This result indicates that the current state of the hull with low pH, high salt, and rich sulfur iron compounds is conducive to the degradation of wood by F. solani NK-NH1. This also explains from another perspective the dominant colonization of NK-NH1 on the hull of the Nanhai No. 1 shipwreck. This is consistent with some previous literature reports (Barclay et al., 1998; Firdous and Shahzad, 2001; Mandeel, 2006; Merlin et al., 2013; Zheng et al., 2017; Dutta et al., 2018).

In addition to cellulose and lignin-degrading enzymes, fungi can also accelerate the degradation of cellulose and lignin through the Fenton reaction (Baldrian and Valášková, 2008). Fe2+, H2O2, and acidic conditions are necessary conditions for the occurrence of Fenton reaction. There are multiple metabolic pathways in organisms that can produce H2O2. The hull is rich in sulfur and iron compounds. Therefore, it is highly likely that F. solani NK-NH1 will produce a Fenton reaction to degrade the hull wood. Through whole genome sequencing, key genes that may be involved in sulfur and iron metabolism in NK-NH1 are found, which are mainly annotated in NR, GO-IPR, KOG, KEGG, NOG, GO, TCDB, and Swiss-Prot databases (Supplementary Table S3 illustrates the key genes annotated). It mainly includes cysteine desulfurizing enzyme, sulfite oxidase, FeS cluster-related protein, sulfite transport and efflux-related protein, and thiosulfate sulfur transferase.

In summary, the current dominant disease fungus, F. solani NK-NH1 promoted the degradation of hull wood. Moreover, the existing low pH, high salt, and sulfur-rich iron compounds preservation state of the Nanhai No. 1 shipwreck will, to a certain extent, promote its degradation of hull wood. Therefore, regular monitoring, physical control, and biochemical antibacterial methods must be combined to address the microbial diseases of the Nanhai No. 1 shipwreck.

5 Conclusion

This study analyzed the microbial diseases of the Nanhai No. 1 shipwreck and biodegradation mechanism of the dominant fungus. The high-throughput sequencing results showed that Fusarium is the dominant fungus, and whole genome sequencing results of the isolated F. solani NK-NH1 showed that it can degrade wood. The highest enzyme activity of the NK-NH1 dominant degrading enzyme LiP is generated under conditions of pH = 4, NaCl concentration of 30%, FeSO4 concentration of 50 mg/L, and no addition of PEG400. The highest enzyme activity of Lac is generated under conditions of pH = 4, NaCl concentration of 10%, FeSO4 concentration of 100 mg/L, and no PEG400 added. This indicates that the current preservation state of the hull with low pH, high salt, and sulfur-rich iron compounds is conducive to the biodegradation of hull wood by F. solani NK-NH1.

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

YW: Formal analysis, Investigation, Software, Writing – original draft, Writing – review & editing. YH: Investigation, Methodology, Writing – original draft. NL: Investigation, Writing – original draft. CW: Writing – original draft. KM: Writing – original draft. XH: Writing – original draft. JD: Supervision, Writing – original draft. HG: Supervision, Validation, Visualization, Writing – review & editing. JP: Conceptualization, Funding acquisition, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Key R&D Program of China (2023YFF0906400) and the Fundamental Research Funds for the Central Universities from the University of Science and Technology Beijing.

Acknowledgments

We acknowledge the assistance of Xingzhong Liu, Dongsheng Wei, Haikun Ma of Nankai University, and Zhiguo Zhang and Yue Chen of the National Centre for Archaeology.

Conflict of interest

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1382653/full#supplementary-material

Footnotes

References

Baldrian, B., and Valášková, V. (2008). Degradation of cellulose by basidiomycetous fungi. FEMS Microbiol. Rev. 32, 501–521. doi: 10.1111/j.1574-6976.2008.00106.x

Crossref Full Text | Google Scholar

Barclay, M., Tett, V. A., and Knowles, C. J. (1998). Metabolism and enzymology of cyanide/metallo cyanide biodegradation by Fusarium solani under neutral and acidic conditions. Enzym. Microb. Technol. 23, 321–330. doi: 10.1016/S0141-0229(98)00055-6

Crossref Full Text | Google Scholar

Benaddou, M., Hajjaj, H., Allali, A., Moubchir, T., M’Barek, H. N., Nafidi, H. A., et al. (2024). Optimizing fungal treatment of lignocellulosic agro-industrial by-products to enhance their nutritional value. Food Sci. Nutr. doi: 10.1002/fsn3.4131

Crossref Full Text | Google Scholar

Diaz-Herraiz, M., Jurado, V., Cuezva, S., Laiz, L., Pallecchi, P., Tiano, P., et al. (2013). The actinobacterial colonization of Etruscan paintings. Sci. Rep. 3, 1440–1445. doi: 10.1038/srep01440

PubMed Abstract | Crossref Full Text | Google Scholar

Dutta, S. D., Tarafder, M., Islam, R., and Datta, B. (2018). Characterization of cellulolytic enzymes of Fusarium soil isolates. Biocatal. Agric. Biotechnol. 14, 279–285. doi: 10.1016/j.bcab.2018.03.011

Crossref Full Text | Google Scholar

Erazo, J. G., Palacios, S. A., Pastor, N., Giordano, F. D., Rovera, M., Reynoso, M. M., et al. (2021). Biocontrol mechanisms of Trichoderma harzianum ITEM 3636 against peanut brown root rot caused by Fusarium solani RC 386. Biol. Control 164:104774. doi: 10.1016/j.biocontrol.2021.104774

Crossref Full Text | Google Scholar

Firdous, H., and Shahzad, S. (2001). Effect of some salts on in vitro growth of Fusarium solani. Pak. J. Bot. 33, 117–124,

Google Scholar

Fors, Y. (2008). Sulfur-related conservation concerns in marine archaeological wood: the origin, speciation and distribution of accumulated sulfur with some remedies for the vasa Stockholm: University of Stockholm.

Google Scholar

Gaylarde, C. C., Morton, L. H. G., Loh, K., and Shirakawa, M. A. (2011). Biodeterioration of external architectural paint films – a review. Int. Biodeterior. Biodegradation 65, 1189–1198. doi: 10.1016/j.ibiod.2011.09.005

Crossref Full Text | Google Scholar

Gutarowska, B., Pietrzak, K., Machnowski, W., and Miczarek, J. M. (2017). Historical textiles – a review of microbial deterioration analysis and disinfection methods. Text. Res. J. 87, 2388–2406. doi: 10.1177/0040517516669076

Crossref Full Text | Google Scholar

Han, Y., Huang, X., Wang, Y., Du, J., Ma, K., Chen, Y., et al. (2021). Fungal community and biodeterioration analysis of hull wood and its storage environment of the Nanhai no. 1 shipwreck. Front. Microbiol. 11:609475. doi: 10.3389/fmicb.2020.609475

PubMed Abstract | Crossref Full Text | Google Scholar

Hocker, E. (2006). From the Micro- to the macro-: managing the conservation of the warship, Vasa. Macromolecular Symposia 238, 16–21. doi: 10.1002/masy.200650603

Crossref Full Text | Google Scholar

Kip, N., and Veen, J. A. (2015). The dual role of microbes in corrosion. ISME J. 9, 542–551. doi: 10.1038/ismej.2014.169

PubMed Abstract | Crossref Full Text | Google Scholar

Köller, W., Allan, C. R., and Kolattukudy, P. E. (1982). Role of cutinase and cell wall degrading enzymes in infection of Pisum sativum by Fusarium solani f. sp. pisi. Physiological Plant Pathol. 20, 47–60. doi: 10.1016/0048-4059(82)90023-6

Crossref Full Text | Google Scholar

Li, Q., Zhang, B., Yang, X., and Ge, Q. (2018). Deterioration-associated microbiome of stone monuments: structure, variation, and assembly. Appl. Environ. Microbiol. 84, e02680–e02617. doi: 10.1128/AEM.02680-17

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Z., Fu, T., Hu, C., Shen, D., Macchioni, N., Sozzi, L., et al. (2018b). Microbial community analysis and biodeterioration of waterlogged archaeological wood from the Nanhai no. 1 shipwreck during storage. Sci. Rep. 8:7170. doi: 10.1038/s41598-018-25484-8

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Z., Wang, Y., Pan, X., Ge, Q., Ma, Q., Li, Q., et al. (2017). Identification of fungal communities associated with the biodeterioration of waterlogged archeological wood in a Han dynasty tomb in China. Front. Microbiol. 8:1633. doi: 10.3389/fmicb.2017.01633

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Z., Zhang, Y., Zhang, F., Hu, C., Liu, G., and Pan, J. (2018a). Microbial community analyses of the deteriorated storeroom objects in the Tianjin museum using culture-independent and culture-dependent approaches. Front. Microbiol. 9:802. doi: 10.3389/fmicb.2018.00802

Crossref Full Text | Google Scholar

M’Barek, H. N., Taidi, B., Smaoui, T., Aziz, M. B., Mansouri, A., and Hajjaj, H. (2019). Isolation. Screening and identification of ligno-cellulolytic fungi from northern Central Morocco. Biotechnol. Agron. Soc. Environ. 23, 207–217,

Google Scholar

Ma, K., Wang, Y., Yang, X., Wang, C., Han, Y., Huang, X., et al. (2022). Analysis of the composition of culturable airborne microorganisms in the archaeological excavation protection site of the Nanhai no. 1 ancient shipwreck. Front. Microbiol. 13:958914. doi: 10.3389/fmicb.2022.958914

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, W., Wu, F., Tian, T., He, D., Zhang, Q., Gu, J., et al. (2020). Fungal diversity and its contribution to the biodeterioration of mural paintings in two 1700-year-old tombs of China. Int. Biodeterior. Biodegradation 152:104972. doi: 10.1016/j.ibiod.2020.104972

Crossref Full Text | Google Scholar

Mandeel, Q. A. (2006). Biodiversity of the genus Fusarium in saline soil habitats. J. Basic Microbiol. 46, 480–494. doi: 10.1002/jobm.200510128

Crossref Full Text | Google Scholar

Mazzoli, R., Giuffrida, M. G., and Pessione, E. (2018). Back to the past: “find the guilty bug-microorganisms involved in the biodeterioration of archeological and historical artifacts”. Appl. Microbiol. Biotechnol. 102, 6393–6407. doi: 10.1007/s00253-018-9113-3

PubMed Abstract | Crossref Full Text | Google Scholar

Merlin, J. N., Christhudas, I. V. S. N., Kumar, P. P., and Agastian, P. (2013). Optimization of growth and bioactive metabolite production: Fusarium solani. Asian J. Pharmaceut. Clin. Res. 6, 98–103,

Google Scholar

Norris, D. M. (1980). Degradation of 14C-labeled lignins and 14C-labeled aromatic acids by Fusarium solani. Appl. Environ. Microbiol. 40, 376–380. doi: 10.1128/aem.40.2.376-380.1980

PubMed Abstract | Crossref Full Text | Google Scholar

Obruca, S., Marova, I., Matoušková, P., and Hároniková, A. (2012). Production of lignocellulose-degrading enzymes employing Fusarium solani F-552. Folia Microbiol. 57, 221–227. doi: 10.1007/s12223-012-0098-5

PubMed Abstract | Crossref Full Text | Google Scholar

Osman, M. E. S., El-Shaphy, A. A. E. N., Meligy, D. A., and Ayid, M. M. (2015). Survey for fungal decaying archaeological wood and their enzymatic activity. Int. J. Conserv. Sci. 5, 295–308,

Google Scholar

Paccanaro, M. C., Sella, L., Castiglioni, C., Giacomello, F., Martínez-Rocha, A. L., D'Ovidio, R., et al. (2017). Synergistic effect of different plant cell wall degrading enzymes is important for virulence of Fusarium graminearum. Mol. Plant-Microbe Interact. 30, 886–895. doi: 10.1094/MPMI-07-17-0179-R

PubMed Abstract | Crossref Full Text | Google Scholar

Panagiotou, G., Kekos, D., Macris, B. J., and Christakopoulos, P. (2003). Production of cellulolytic and xylanolytic enzymes by Fusarium oxysporum grown on corn Stover in solid state fermentation. Ind. Crop. Prod. 18, 37–45. doi: 10.1016/S0926-6690(03)00018-9

Crossref Full Text | Google Scholar

Passialis, C. N. (1997). Physico-chemical characteristics of waterlogged archaeological wood. Environ. Sci. 51, 111–113. doi: 10.1515/hfsg.1997.51.2.111

Crossref Full Text | Google Scholar

Piñar, G., Ripka, K., Weber, J., and Sterflinger, K. (2009). The micro-biota of a sub-surface monument the medieval chapel of St. Virgil (Vienna, Austria). Int. Biodeterior. Biodegradation 63, 851–859. doi: 10.1016/j.ibiod.2009.02.004

Crossref Full Text | Google Scholar

Preston, J., Smith, A. D., Schofield, E. J., Chadwick, A. V., Jones, M. A., and Watts, J. E. M. (2014). The effects of Mary rose conservation treatment on iron oxidation processes and microbial communities contributing to acid production in marine archaeological timbers. PLoS One 9, 1–8. doi: 10.1371/journal.pone.0084169

PubMed Abstract | Crossref Full Text | Google Scholar

Rahmawati, N., Sumardi, I., and Dungani, R. (2020). Isolation and identification of fungi inhabiting rubber-impregnated wood, and their role of quality changing the impregnated wood. BioRes. 15, 2839–2849. doi: 10.15376/biores.15.2.2839-2849

Crossref Full Text | Google Scholar

Rodriguez, A., Perestelo, F., Carnicero, A., Regalado, V., Perez, R., Fuente, G., et al. (1996). Degradation of natural lignins and lignocellulosic substrates by soil-inhabiting fungi imperfecti. FEMS Microbiol. Ecol. 21, 213–219. doi: 10.1111/j.1574-6941.1996.tb00348.x

Crossref Full Text | Google Scholar

Ruiz, G. B., Pietro, A. D., and Roncero, M. I. G. (2016). Combined action of the major secreted exo- and endopolygalacturonases is required for full virulence of Fusarium oxysporum. Mol. Plant Pathol. 17, 339–353. doi: 10.1111/mpp.12283

PubMed Abstract | Crossref Full Text | Google Scholar

Saparrat, M. C. N., Martínez, M. J., Tournier, H. A., Cabello, M. N., and Arambarri, A. M. (2000). Production of ligninolytic enzymes by Fusarium solani strains isolated from different substrata. World J. Microbiol. Biotechnol. 16, 799–803. doi: 10.1023/A:1008923804859

Crossref Full Text | Google Scholar

Si, J., Ma, H. F., Cao, Y. J., Cui, B. K., and Dai, Y. C. (2021b). Introducing a thermo-alkali-stable, metallic ion-tolerant laccase purified from white rot fungus Trametes hirsuta. Front. Microbiol. 12:670163. doi: 10.3389/fmicb.2021.670163

PubMed Abstract | Crossref Full Text | Google Scholar

Si, J., Wu, Y., Ma, H. F., Cao, Y. J., Sun, Y. F., and Cui, B. K. (2021a). Selection of a pH-and temperature-stable laccase from Ganoderma australe and its application for bioremediation of textile dyes. J. Environ. Manag. 299:113619. doi: 10.1016/j.jenvman.2021.113619

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, A. P. (2012). A review of microbial decay types found in wooden objects of cultural heritage recovered from buried and waterlogged environments. J. Cult. Herit. 13, S16–S20. doi: 10.1016/j.culher.2012.04.002

Crossref Full Text | Google Scholar

Sun, M., Zhang, F., Huang, X., Han, Y., Jiang, N., Cui, B., et al. (2020). Analysis of microbial community in the archaeological ruins of Liangzhu city and study on protective materials. Front. Microbiol. 11:684. doi: 10.3389/fmicb.2020.00684

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, B., Zhu, C., Hu, Y., Zhang, B., and Wang, J. (2023). Dynamics of microbial community composition during degradation of silks in burial environment. Sci. Total Environ. 883:163694. doi: 10.1016/j.scitotenv.2023.163694

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, Y., Luo, Z., Chow, R. K., and Vrijmoed, L. L. P. (2010). Purification and characterization of an extracellular laccase from the anthracene-degrading fungus Fusarium solani MAS2. Bioresour. Technol. 101, 9772–9777. doi: 10.1016/j.biortech.2010.07.091

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, F. (1997). Effects of redox potential and hydroxide inhibition on the pH activity profile of fungal laccases. J. Biol. Chem. 272, 924–928. doi: 10.1074/jbc.272.2.924

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, F., Li, L., Sun, M., Hu, C., Zhang, Z., Liu, Z., et al. (2019). Fungal community analyses of a pirogue from the tang dynasty in the national maritime museum of China. Appl. Sci. 9:4129. doi: 10.3390/app9194129

Crossref Full Text | Google Scholar

Zhang, Z., Li, N., Tian, X., Liu, J., and Shen, D. (2014). Research on the removal of the iron sulfides in the Qing dynasty marine shipwreck, Ningbo Xiaobaijiao no. 1. Sci. Conservat. Archaeol. 26, 30–38,

Google Scholar

Zheng, F., An, Q., Meng, G., Wu, X. J., Dai, Y. C., Si, J., et al. (2017). A novel laccase from white rot fungus Trametes orientalis: purification, characterization, and application. Int. J. Biol. Macromol. 102, 758–770. doi: 10.1016/j.ijbiomac.2017.04.089

Crossref Full Text | Google Scholar

Keywords: the Nanhai No. 1 shipwreck, Fusarium solani, biodegradation mechanism, microbial community analysis, cultural relics preservation

Citation: Wang Y, Han Y, Li N, Wang C, Ma K, Huang X, Du J, Guo H and Pan J (2024) Study on biodegradation mechanism of Fusarium solani NK-NH1 on the hull wood of the Nanhai No. 1 shipwreck. Front. Microbiol. 15:1382653. doi: 10.3389/fmicb.2024.1382653

Received: 06 February 2024; Accepted: 29 April 2024;
Published: 30 May 2024.

Edited by:

Hao Zhou, Dalian University of Technology, China

Reviewed by:

Federica Spina, University of Turin, Italy
Jing Si, Beijing Forestry University, China

Copyright © 2024 Wang, Han, Li, Wang, Ma, Huang, Du, Guo and Pan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jiao Pan, cGFuamlhb0B1c3RiLmVkdS5jbg==; Hong Guo, Z3VvaG9uZzU0MzIxQDE2My5jb20=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.