- 1College of Coastal Agricultural Sciences, Guangdong Ocean University, Zhanjiang, China
- 2Maoming Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Maoming, China
In order to explore the impact of antibiotics (enrofloxacin) on microbial community in aquatic environment, an indoor aquatic ecological model was built, and different concentrations of enrofloxacin (0.05, 0.5, 5, and 50 mg/L) were added in the aquatic ecological model. In addition, the water and sediment samples were collected on the 0, 7, 30, and 60 days, and the changes in microbial community were studied through 16S rDNA high-throughput sequencing. The results showed that when the concentration of enrofloxacin was 50 mg/L, the relative abundance of Actinomycetes was increased. In the water, the bacterial richness and diversity communities first decreased and then gradually recovered with the passage of time; On the 7th day, the diversity and richness index of species in the treatment groups with enrofloxacin at 5 and 50 mg/L decreased to the lowest; On the 30th day, the diversity and richness index of species began to rise; On the 60th day, the diversity index and richness index of water species began to increase, while the diversity index and richness index of sediment species decreased. In conclusion, the addition of enrofloxacin negatively affected the microbial community structure in an indoor aquatic ecological model, 50 mg/L enrofloxacin could increase the relative abundance of Actinomycetes, and decrease the diversity and richness index of water and sediment.
Introduction
Antibiotics are widely used in the fields of medicine, aquaculture, and animal husbandry as pathogenic microorganism inhibitors and/or growth promoters during recent decades (1). However, about 30–90% of antibiotics in humans and animals are discharged into the body through urine and feces, which are often not fully utilized and treated, resulting in a large number of residues of antibiotics and their metabolites in the soil, water, and vegetation of the environment (2–4). It has been shown that the concentration of antibiotic residues directly affects the structure and community of microorganisms in the ecosystem (5). Enrofloxacin is a kind of fluoroquinolone antibiotics, which has strong antibacterial activity against most gram-negative bacteria, some gram-positive bacteria and mycoplasma (6, 7). After the non-covalent combination of two enrofloxacin molecules with DNA-topoisomerase complex (II or IV) near the tyrosine residue of the active site of bacteria, its induced enzyme underwent conformational changes, which led to the formation of enrofloxacin-rotatase/topoisomerase IV-DNA complex, and directly inhibited bacterial DNA replication (8, 9). Low concentration of enrofloxacin would inhibit bacterial DNA rotatase and topoisomerase IV, but high concentration of enrofloxacin can directly destroy bacterial chromosomes, leading to cell death (10, 11).
Because of the excellent and broad-spectrum antibacterial effect, enrofloxacin is frequently used as a therapeutic antibiotic and growth promoter in aquaculture, resulting in its metabolites remaining in the soil and aquaculture environment. According to the report of Li et al. (12), they detected nine kinds of quinolones in the seawater of aquaculture areas, and the concentration range of enrofloxacin residue was 6.91–17.49 ng/L. The main place where norfloxacin accumulates in aquatic environments is surface water, such as rivers, where the level is between 12 ng/L and 4.24 μg/L (13), and the concentrations of streams and ponds are about 17–216 ng/L (14) and 0.50 μg/L (15) respectively. Besides, it is reported that enrofloxacin with concentration ranging from 2.0 to 4.0 ng/L was detected in the tap water of Guangzhou and Macao (16). The largest residues of enrofloxacin are domestic wastewater, municipal wastewater and hospital wastewater, and the concentration of enrofloxacin may be as high as 100 μg/L. (17) Therefore, the toxicity of enrofloxacin to the ecological environment and its impact on the microbial community and structure in the environment are of great concern. The abuse of enrofloxacin in aquaculture has a negative impact on the aquatic ecological environment and may affect human public health through the ecological cycle (18). Microorganisms are an important part of the aquatic environment, which play an important role in the response to antibiotics, the spread of drug resistance, and environmental remediation (19). It has been found that the antibiotics could directly affect the drug resistance of environmental microflora, thus causing environmental toxicity (20, 21). However, the relevant reports about the impact of enrofloxacin on the microbial community in the aquatic environment are limited. Therefore, this study evaluated the effects of different concentrations of enrofloxacin on bacterial community in an indoor aquatic ecological model.
Materials and methods
Reagents
The lake surface water (1–10 cm) and lake surface sediments (1–10 cm) in Guangdong Ocean University (Zhanjiang, China) were collected, and the debris was removed to construct an aquatic microcosm. Enrofloxacin (99.95%) was purchased from Zhejiang Guobang Pharmaceutical Co., Ltd. (190720-2, Hangzhou, China). Water DNA Kit (D5525), Soil DNA Kit (D5625), and MicroElute Gel Extraction Kit (D6294) were obtained from Omega Bio-Tek Company (Guangzhou, China).
Experimental design and sampling
The indoor aquatic ecological model was constructed according to our previous study, and the details of the aquatic ecological model have been reported previously (20). Briefly, the water and sediment used for constructing aquatic models are all taken from the same location of the artificial lake on the campus of Guangdong Ocean University (Zhanjiang, China). The 60 (length) × 50 (width) × 30 (height) cm transparent plastic boxes were used for constructing aquatic model. First, lay the sediment at the bottom of the plastic box, so that the height of the sediment is 10 cm, and then add 45 L of water. The illustration of the indoor aquatic ecological model is presented in Supplementary Figure S1. The experiment was divided into five groups, with three replicates in each group, a total of 15 aquatic model boxes were used. Enrofloxacin was added to the water at 0, 0.05, 0.5, 5, and 50 mg/L (water and sediment were W0–W4 and S0–S4, respectively), and was cultured at room temperature of 25 ± 3°C for 60 days. Samples of water (Group W) and sediment (Group S) were collected at 0, 7, 30, and 60 days (D0, D7, D30, and D60). The 5% sodium hydroxide solution was used as the solvent for enrofloxacin. The specific method for adding different concentrations of enrofloxacin is as follows: take 5 g of enrofloxacin, add an appropriate amount of 5% sodium hydroxide solution until it is completely dissolved, then add double distilled water to a constant volume until the concentration of enrofloxacin is 5 g/L. Then dilute it with double distilled water to several concentrations of 5, 50, and 500 mg/L, and then 30 L of water was taken from each aquatic ecological model. After discarding 450 mL of water, 450 mL of corresponding concentration of enrofloxacin solution was added. For the control group, 450 mL of double distilled water without enrofloxacin was added. So that the enrofloxacin concentrations in the water bodies of each model were 0, 0.05, 0.5, 5, and 50 mg/L.
16S rDNA sequencing and statistical analysis
The supernatant of water and sediment samples was removed by centrifugation, and the total genomic DNA of water and sediment microorganisms was extracted according to the instructions of the kit. The concentration and purity of extracted DNA were determined by Nanodrop UV–Vis spectrophotometer. The extracted DNA was sent to Suzhou Golden Wisdom Biotechnology Co., Ltd. for Illumina MiSeq sequencing. QIIME was used to optimize the original sequencing data: splicing the overlapping regions at the end of the sequence; Remove the sequence with length less than 200 bp; Remove the chimera sequence and get valid data. OTU cluster analysis was performed according to 97% similarity. The Silva 138 16S rDNA database was used for species classification annotation. The statistical analysis of the 16S rDNA sequencing data was reported in our earlier study (20).
Results
Microbial diversity analysis
The operational taxonomic units (OTUs) of bacteria were presented in Figure S2, the rarefaction curve of the OTUs was showed in Figure S3, and the unweighted unifrac distance matrix heat map was showed in Figure S4, and the principal coordinate (PCoA) analysis and principal component (PCA) analysis were presented in Figure S5 and S6. In addition, The community richness index mainly includes Chao1 index and ACE index, and the flora diversity index mainly includes Shannon index and Simpson index. As shown in Supplementary Table S1, among the sediment samples, the diversity index values of S4D60 sample are the lowest; As shown in Supplementary Table S2, in the water, the diversity index values of W3D7 samples are the lowest; It shows that enrofloxacin can inhibit the richness and diversity of bacterial species in S4D60 and W3D7 samples. With the increase of drug concentration in the water body, both the richness and diversity of the colony showed a downward trend. At a higher drug concentration, with the prolongation of the action time, the richness and diversity of the colony showed a downward trend and then gradually rose. At 30 days, the diversity index began to rise, and the higher the drug concentration, the smaller the index recovery trend.
Analysis of bacterial abundance and distribution
The analysis of non-metric multidimensional scale method (NMDS) and unweighted pair group method with arithmetical mean (UPGMA) was presented in Figure S7 and S8. The relative abundance of microbials was analyzed at the phylum level, and the information of the first 30 species is selected to draw a histogram as shown in Figure 1. The dominant bacteria in the water are Proteobacteria, Actinobacteriota, and Bacteroidota, and the dominant bacteria in the sediment are Proteobacteria, Bacteroidota, Desulfobacterota, Acidobacter, and Myxococcota. The water and sediment are the largest species of Proteobacteria, and the relative abundance of Proteobacteria in the water is 25–73%, and the relative abundance of Proteobacteria in the sediment is 24–44%. With the change of time, the relative abundance of Proteobacteria increased first and then decreased; the relative abundance of Actinomycetes and Bacteroides decreased first and then increased with time. In water, the relative abundance of Bacteroides in W3D7 samples increased; the relative abundance of Actinomycetes in W4D60 samples increased, while the relative abundance of Proteobacteria decreased. In the sediment, the relative abundance of Myxococcus in S4D60 samples decreased.
Figure 1. Effects of different concentrations of enrofloxacin on microbial abundance at the phylum level in an aquatic ecological model. W, water; S, sediments.
At the genus level, a total of 413 genera were obtained from sediments, with the number of each group being 275–349, and 248 genera were obtained from water samples, with the number of each group being 117–181. The first 30 genera of horizontal species were selected to draw the Heatmap map, as shown in Figure 2. The main dominant genera in the water without enrofloxacin were Proteobacteria C39, Chlorophyllaceae and Comamandaceae; HgCI-class, CL500-29-marine-group, Myxococcota, PeM15 of Actinobacteriota; Pseudoricicella of Bacteroidota, etc. In water, the relative abundance of Actinomycete PeM15 in W4D60 samples increased. The main dominant bacteria in the sediment are Dechromonas, Anaeromyxoactor, Bacteroidetes-vadinHA17, Geobactaceae, Geothrix, Steroidobactaceae, etc.
Figure 2. Effects of different concentrations of enrofloxacin on microbial abundance at the genus level in an aquatic ecological model. W, water; S, sediments.
Discussion
There are a variety of microorganisms in the aquatic ecosystem, which play a key role in the aquatic environment, and carry out the circulation between material exchanges and maintain the stability of the water environment (22, 23). In recent years, the abuse of antibiotics has left a large number of residues in the aquatic environment, and the microorganisms in the aquatic environment will obtain the exogenous resistance genes through gene mutation, generation inheritance, transduction, conjugation, and transformation (24, 25). Some studies have also described that under the low concentration of antibiotics, microorganisms change the microbial community as a signal molecule so as to adapt to the residues of different concentrations of antibiotics in the aquatic environment (26–28).
In this experiment, the aquatic environment with enrofloxacin concentration less than 0.5 mg/L has no effect on the diversity and abundance of bacterial, while the concentration of enrofloxacin at 5 mg/L has a greater impact on bacterial community. Wei et al. (29) found that the increased stress concentration of enrofloxacin added to soil would significantly affect the soil bacterial community and change the ability of microorganisms to use carbon sources. In addition, fluoroquinolone antibiotics show resistance to hydrolysis and thermal coupling, and have good chemical stability, so there is a high degree of persistent pollution in the complex aquatic environment (30). In aquaculture, the antibiotics in the water will be degraded quickly due to the effect of light (31). However, when enrofloxacin and its metabolite ciprofloxacin enter the water body, 65% will be absorbed by the sediment (32). It has been suggested that the degradation of enrofloxacin is mainly photolysis, so enrofloxacin entering the sediment can be maintained for a long time, and can be slowly released into the water body through the desorption of the sediment (33). This partly explained that enrofloxacin had a greater impact on the diversity and abundance of microbial communities in water and a smaller impact on microbial communities in sediment in this experiment. In this study, we found that when the concentration of enrofloxacin was 50 mg/L, the relative abundance of Actinomycetes increased. According to previous studies, the presence of antibiotics in the aquatic environment will change the microbial community structure and diversity of the biofilm, and with the increase of antibiotic pollution concentration, the content of Actinomycetes in the aquatic environment will also increase (34). This indicates that high concentration of enrofloxacin can increase the relative abundance of drug-resistant bacteria. Therefore, in the aquaculture ecosystem, the residues of antibiotics will induce the microbial community to produce resistance genes, and the number of drug-resistant bacteria will increase significantly, polluting the drinking water source and animal-derived food, and posing a potential threat to human health (35).
In summary, with the increase of enrofloxacin concentration in the water, the richness and diversity of microbial community showed a downward trend. At a higher enrofloxacin concentration, with the prolongation of action time, the richness and diversity of bacterial community showed a downward trend and then gradually rose. The diversity and richness of microbial community in sediment are higher than that in water. The higher concentration of enrofloxacin added, the greater impact on the community and structure of microorganisms in the aquatic environment, and will increase the number of drug-resistant bacteria.
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 at: https://www.ncbi.nlm.nih.gov/, bioproject/926010.
Author contributions
YM and W-CL: funding acquisition and writing—review and editing. YD: investigation. T-YZ and J-JP: methodology. X-PX and YD: project administration. S-SL: software. YM: supervision. YD and J-JP: writing—original draft. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by Guangdong Natural Science Foundation Project (2023A1515012181); Guangdong Provincial Department of Education 2021 Special Project for Key Fields of Ordinary Colleges and Universities (2021ZDZX4003); Zhanjiang Science and Technology Bureau 2021 Provincial Science and Technology Special Funds (“Big Special + Task List”) Competitive Allocation Project (2021A05231); and Start-up Research Project of Maoming Laboratory, Guangdong Laboratory for Lingnan Modern Agriculture (2021TDQD002).
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/fvets.2023.1151988/full#supplementary-material
References
1. Zuccato, E, Castiglioni, S, Bagnati, R, Melis, M, and Fanelli, R. Source, occurrence and fate of antibiotics in the Italian aquatic environment. J Hazard Mater. (2010) 179:1042–8. doi: 10.1016/j.jhazmat.2010.03.110
2. Zhang, Q-Q, Ying, G-G, Pan, C-G, Liu, Y-S, and Zhao, J-L. Comprehensive evaluation of antibiotics emission and fate in the river basins of China: source analysis, multimedia modeling, and linkage to bacterial resistance. Environ Sci Technol. (2015) 49:6772–82. doi: 10.1021/acs.est.5b00729
3. Li, H, Tian, Y, Liu, W, Long, Y, Ye, J, Li, B, et al. Impact of electrokinetic remediation of heavy metal contamination on antibiotic resistance in soil. Chem Eng J. (2020) 400:125866. doi: 10.1016/j.cej.2020.125866
4. Kim, K-R, Owens, G, Kwon, S-I, So, K-H, Lee, D-B, and Ok, YS. Occurrence and environmental fate of veterinary antibiotics in the terrestrial environment. Water Air Soil Pollut. (2011) 214:163–74. doi: 10.1007/s11270-010-0412-2
5. Chow, LK, Ghaly, TM, and Gillings, MR. A survey of sub-inhibitory concentrations of antibiotics in the environment. J Environ Sci. (2021) 99:21–7. doi: 10.1016/j.jes.2020.05.030
6. Liu, X, Yang, Q, Fan, Y, Du, Y, Lei, L, Wang, D, et al. Pharmacokinetics and pharmacodynamics of enrofloxacin treatment of Escherichia coli in a murine thigh infection modeling. BMC Vet Res. (2021) 17:1–13. doi: 10.1186/s12917-021-02908-8
8. Hooper, DC, and Jacoby, GA. Topoisomerase inhibitors: fluoroquinolone mechanisms of action and resistance. Cold Spring Harb Perspect Med. (2016) 6:a025320. doi: 10.1101/cshperspect.a025320
9. Bush, NG, Diez-Santos, I, Abbott, LR, and Maxwell, A. Quinolones: mechanism, lethality and their contributions to antibiotic resistance. Molecules. (2020) 25:5662. doi: 10.3390/molecules25235662
10. Hooper, DC, and Jacoby, GA. Mechanisms of drug resistance: quinolone resistance. Ann N Y Acad Sci. (2015) 1354:12–31. doi: 10.1111/nyas.12830
11. Wang, X, Zhao, X, Malik, M, and Drlica, K. Contribution of reactive oxygen species to pathways of quinolone-mediated bacterial cell death. J Antimicrob Chemother. (2010) 65:520–4. doi: 10.1093/jac/dkp486
12. Li, Z, Dong, X, and Wu, M. Quinolone residues in seawater of aquaculture area, Sanggou bay, Yellow Sea. Mar Environ Sci. (2018) 2:182–6. doi: 10.13634/j.cnki.mes.2018.02.004
13. Nguyen Dang Giang, C, Sebesvari, Z, Renaud, F, Rosendahl, I, Hoang Minh, Q, and Amelung, W. Occurrence and dissipation of the antibiotics sulfamethoxazole, sulfadiazine, trimethoprim, and enrofloxacin in the Mekong Delta, Vietnam. PLoS One. (2015) 10:e0131855. doi: 10.1371/journal.pone.0131855
14. Ikem, A, Lin, C-H, Broz, B, Kerley, M, and Le Thi, H. Occurrence of enrofloxacin in overflows from animal lot and residential sewage lagoons and a receiving-stream. Heliyon. (2017) 3:e00409. doi: 10.1016/j.heliyon.2017.e00409
15. Wei, R, Ge, F, Chen, M, and Wang, R. Occurrence of ciprofloxacin, enrofloxacin, and florfenicol in animal wastewater and water resources. J Environ Qual. (2012) 41:1481–6. doi: 10.2134/jeq2012.0014
16. Wang, Q-J, Mo, C-H, Li, Y-W, Gao, P, Tai, Y-P, Zhang, Y, et al. Determination of four fluoroquinolone antibiotics in tap water in Guangzhou and Macao. Environ Pollut. (2010) 158:2350–8. doi: 10.1016/j.envpol.2010.03.019
17. Babić, S, Periša, M, and Škorić, I. Photolytic degradation of norfloxacin, enrofloxacin and ciprofloxacin in various aqueous media. Chemosphere. (2013) 91:1635–42. doi: 10.1016/j.chemosphere.2012.12.072
18. Jin, X, Liu, S, Zhang, Z, Liu, T, Li, N, Liang, Y, et al. Enrofloxacin-induced transfer of multiple-antibiotic resistance genes and emergence of novel resistant bacteria in red swamp crayfish guts and pond sediments. J Hazard Mater. (2023) 443:130261. doi: 10.1016/j.jhazmat.2022.130261
19. Zhao, F, Yang, L, Li, G, Fang, L, Yu, X, Tang, Y-T, et al. Veterinary antibiotics can reduce crop yields by modifying soil bacterial community and earthworm population in agro-ecosystems. Sci Total Environ. (2022) 808:152056. doi: 10.1016/j.scitotenv.2021.152056
20. Zhang, T, Ding, Y, Peng, J, Dai, Y, Luo, S, Liu, W, et al. Effects of broad-Spectrum antibiotic (Florfenicol) on resistance genes and bacterial community structure of water and sediments in an aquatic microcosm model. Antibiotics. (2022) 11:1299. doi: 10.3390/antibiotics11101299
21. Liao, C-Y, Balasubramanian, B, Peng, J-J, Tao, S-R, Liu, W-C, and Ma, Y. Antimicrobial resistance of Escherichia coli from aquaculture farms and their environment in Zhanjiang, China. Front Vet Sci. (2021) 8:1534. doi: 10.3389/fvets.2021.806653
22. Yitayew, B, Woldeamanuel, Y, Asrat, D, Rahman, A, Mihret, A, Aseffa, A, et al. Antimicrobial resistance genes in microbiota associated with sediments and water from the Akaki river in Ethiopia. Environ Sci Pollut Res. (2022) 29:70040–55. doi: 10.1007/s11356-022-20684-2
23. Shu, W-S, and Huang, L-N. Microbial diversity in extreme environments. Nat Rev Microbiol. (2022) 20:219–35. doi: 10.1038/s41579-021-00648-y
24. Hossain, A, Habibullah-Al-Mamun, M, Nagano, I, Masunaga, S, Kitazawa, D, and Matsuda, H. Antibiotics, antibiotic-resistant bacteria, and resistance genes in aquaculture: risks, current concern, and future thinking. Environ Sci Pollut Res. (2022) 29:11054–75. doi: 10.1007/s11356-021-17825-4
25. Ogunlaja, A, Ogunlaja, OO, Olukanni, OD, Taylor, K, Olorunnisola, CG, Dougnon, VT, et al. Antibiotic resistomes and their chemical residues in aquatic environments in Africa. Environ Pollut. (2022) 312:119783. doi: 10.1016/j.envpol.2022.119783
26. Li, S, Zhang, C, Li, F, Hua, T, Zhou, Q, and Ho, S-H. Technologies towards antibiotic resistance genes (ARGs) removal from aquatic environment: a critical review. J Hazard Mater. (2021) 411:125148. doi: 10.1016/j.jhazmat.2021.125148
27. Dong, H, Chen, Y, Wang, J, Zhang, Y, Zhang, P, Li, X, et al. Interactions of microplastics and antibiotic resistance genes and their effects on the aquaculture environments. J Hazard Mater. (2021) 403:123961. doi: 10.1016/j.jhazmat.2020.123961
28. Binh, VN, Dang, N, Anh, NTK, and Thai, PK. Antibiotics in the aquatic environment of Vietnam: sources, concentrations, risk and control strategy. Chemosphere. (2018) 197:438–50. doi: 10.1016/j.chemosphere.2018.01.061
29. Wei, Z, Wang, J, Zhu, L, Wang, J, and Zhu, G. Toxicity of enrofloxacin, copper and their interactions on soil microbial populations and ammonia-oxidizing archaea and bacteria. Sci Rep. (2018) 8:5828. doi: 10.1038/s41598-018-24016-8
30. Turiel, E, Martín-Esteban, A, Bordin, G, and Rodríguez, AR. Stability of fluoroquinolone antibiotics in river water samples and in octadecyl silica solid-phase extraction cartridges. Anal Bioanal Chem. (2004) 380:123–8. doi: 10.1007/s00216-004-2730-y
31. Zhang, H, Xie, H, Chen, J, and Zhang, S. Prediction of hydrolysis pathways and kinetics for antibiotics under environmental pH conditions: a quantum chemical study on cephradine. Environ Sci Technol. (2015) 49:1552–8. doi: 10.1021/es505383b
32. Kümmerer, K . Drugs in the environment: emission of drugs, diagnostic aids and disinfectants into wastewater by hospitals in relation to other sources–a review. Chemosphere. (2001) 45:957–69. doi: 10.1016/S0045-6535(01)00144-8
33. Slana, M, Žigon, D, and Sollner-Dolenc, M. Enrofloxacin degradation in broiler chicken manure under field conditions and its residuals effects to the environment. Environ Sci Pollut Res. (2017) 24:13722–31. doi: 10.1007/s11356-017-8722-1
34. Proia, L, Lupini, G, Osorio, V, Pérez, S, Barceló, D, Schwartz, T, et al. Response of biofilm bacterial communities to antibiotic pollutants in a Mediterranean river. Chemosphere. (2013) 92:1126–35. doi: 10.1016/j.chemosphere.2013.01.063
Keywords: aquatic environment, aquatic ecological model, bacterial community structure, enrofloxacin, 16S rDNA
Citation: Dai Y, Peng J-J, Zhang T-Y, Xie X-P, Luo S-S, Liu W-C and Ma Y (2023) Effects of antibiotics (enrofloxacin) on microbial community of water and sediment in an aquatic ecological model. Front. Vet. Sci. 10:1151988. doi: 10.3389/fvets.2023.1151988
Edited by:
Enrico Fiore, University of Padua, ItalyReviewed by:
Yongtao Liu, Chinese Academy of Fishery Sciences (CAFS), ChinaHector Sumano, National Autonomous University of Mexico, Mexico
Copyright © 2023 Dai, Peng, Zhang, Xie, Luo, Liu and Ma. 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: Wen-Chao Liu, bGl1d2NAZ2RvdS5lZHUuY24=; Yi Ma, bWF5aTc2MUBnZG91LmVkdS5jbg==
†These authors have contributed equally to this work and share first authorship