Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 25 April 2022
Sec. Alzheimer's Disease and Related Dementias
This article is part of the Research Topic Epigenetic and Genetic Marks in Human Stem Cell Models of Aging and Neurodegeneration View all 5 articles

Expression Pattern of p62 in Primary Age-Related Tauopathy: Staging of p62 in PART

\r\nXin Wang,&#x;Xin Wang1,2†Lei Zhang,&#x;Lei Zhang1,2†Hui LuHui Lu1Juanli WuJuanli Wu1Huazheng LiangHuazheng Liang3Bing SunBing Sun1Keqing Zhu,*\r\nKeqing Zhu1,2*
  • 1China Brain Bank and Department of Neurology in Second Affiliated Hospital, Key Laboratory of Medical Neurobiology of Zhejiang Province, and Department of Neurobiology, Zhejiang University School of Medicine, Hangzhou, China
  • 2Department of Pathology, Zhejiang University School of Medicine, Hangzhou, China
  • 3Translational Research Institute of Brain and Brain-Like Intelligence, Shanghai Fourth People’s Hospital, School of Medicine, Tongji University, Shanghai, China

The present study analyzed the distribution pattern of p62 immunoreactivity in brains of primary age-related tauopathy (PART) and Braak NFT matched pre-AD and Alzheimer’s disease (AD) patients using immunohistochemistry in combination with semi-quantitative evaluation. In PART and AD brains, p62 was found positive in seven regions, including the neocortex, thalamus, basal ganglia, hippocampus, brainstem, cerebellar dentate nucleus, and the cervical spinal cord. There was a positive correlation between the Braak NFT stage and the distribution of p62 expression. Six stages of expression of p62 were proposed from the present study. Expression of p62 in the hippocampus of PART and AD was classified stage I, the brainstem stage II, the thalamus stage I _I _I, the basal ganglia stage IV, the neocortex stage V, the cerebellum and the cervical spinal cord stage VI. The hippocampus was the site initially affected by p62, especially the CA1 and the subiculum. They might be the earliest accumulation site of p62.

Introduction

P62, also named sequestosome 1 (SQSTM1), is a selective cargo receptor for autophagy to degrade misfolded proteins (Du et al., 2009). A previous study showed that p62 was associated with neurodegenerative diseases, and one of the underlying mechanisms is p62 induced alteration in autophagy (Liu et al., 2017). p62 induced autophagy failure significantly accelerates aggregation of misfolded proteins. Furthermore, accumulation of misfolded proteins leads to aberrant p62 expression (Liu et al., 2017). Screening with p62 immunohistochemistry has the potential to reveal other types of aggregates that have evaded detection by other stains (Kuusisto et al., 2008). p62-immunoreactivity (IR) is present at a very early stage of inclusion formation. Since it exhibits a high staining intensity of the inclusions with optimal contrast toward the background, it permits reliable evaluation of pathologic lesions (Salminen et al., 2012).

Primary age-related tauopathy (PART) is recently defined by the presence of Alzheimer’s disease (AD)-like neurofibrillary changes without, or with few Aβ plaques. PART can be designated as “Definite” or “Possible” depending on the presence of Aβ plaques. The working classification for definite PART is based on the presence of neurofibrillary tangles (NFT) along with Braak stage ≤ IV and Thal Aβ Phase as 0 (Crary et al., 2014). Based on the neuropathological diagnostic criteria, we found that PART is common in the Chinese population as observed from specimens in our newly established brain bank. It is still unclear whether PART is an early form of AD or a distinct tauopathy different from AD (Duyckaerts et al., 2015; Jellinger et al., 2015). Here, we described the anatomic distribution of p62-IR in PART and AD subjects using immunohistochemistry as a screening stain for diagnostic practice.

Materials

A total of 18 definite PART and 11 AD cases from the Chinese Brain Bank of School of Medicine, Zhejiang University were examined. All subjects fulfilled the following inclusion criteria: age at death ≥ 50 years, and postmortem delay less than 24 h. The diagnosis of PART was confirmed by a neuropathological investigation with NFT Braak stage ≤ IV and Thal Aβ Phase as 0. AD was confirmed by Braak NFT stage ≥ IV and CERAD plaque density C and pre-AD by Braak NFT stage I _I _I and CERAD plaque density C or B. When used in the NIA-AA diagnostic criteria for AD, the CERAD plaque densities scores are: none = 0, sparse = A, moderate = B, frequent = C (Mirra et al., 1991).

Immunohistochemistry

Immunohistochemistry of tissues about 605 pieces in total was performed on paraformaldehyde (PFA)-fixed, paraffin-embedded tissue from all autopsy cases. Small blocks of brain were dissected at autopsy and fixed in 4% PFA in 0.1 M phosphate buffer (pH 7.4) for 2 days. Following cryoprotection in 15% sucrose in 0.01 M phosphate-buffered saline (PBS, pH 7.4), blocks were cut into 3 μm thick sections using a microtome and mounted onto slides. Sections were incubated with 3% H2O2 for 10 min to eliminate endogenous peroxidase activity in the tissue. Prior to immunostaining, sections underwent microwave antigen retrieval for 15 min in the citrate buffer (pH 6.0). After washing with PBS containing 0.3% Triton X-100 (Tx-PBS) for 30 min, sections were blocked with 10% normal serum from the same origin of the secondary antibody, and then incubated with the primary antibody (mouse anti-p62, Biosciences, 610833, 1:100) for 24 h at 4°C. Following incubation with the appropriate secondary antibody, labeling was detected using the avidin–biotinylated HRP complex (ABC) system (Vector Laboratories, Burlingame, CA, United States). The peroxidase reaction was carried out using a developer solution containing 0.4 mg/ml DAB and 0.0006% hydrogen peroxide dissolved in TBS. For the negative control, the primary antibodies were omitted and all other steps carried out as described above. After the staining procedures, sections dehydrated and cleared before being coverslipped with DPX mounting medium. Expression of p62 pathology was semi-quantitatively scored based on a four-point grading scale (-/0, none;+/1, mild; ++/2, moderate;+++/3, severe). Of note, all cases were initially checked using the routine protocol for the Braak staging system (Braak et al., 2006). All images were collected with Virtual Slide System 120 microscope in 10× magnification and analyzed with CellSens Dimension software.

Ethics Approval and Consent to Participate

Use of tissue for this study was approved by the Medical Ethics Committee of School of Medicine, Zhejiang University. Written consent was obtained with all participants.

Statistical Methods

All sections were scanned by Virtual Slide System 120 microscope in 10× magnification and analyzed with CellSens Dimension software. IBM SPSS version statistics software (Armonk, NY, United States) was used for statistical analysis, and the standard error of the means was used when applicable. Correlation between the studied variables was assessed using the Spearman correlation analysis. Mann–Whitney U (two groups) test and Kruskal–Wallis H test (three or more groups) were implemented to test the statistical difference between the studied groups. To calculate the quantitative p62 score in hippocampus, three square (100 μm × 100 μm) with same length were selected in each area according to the organizational structures, and the number of p62 positive cells were counted as p62 score.

Results

1. Expression Pattern of p62 in Primary Age-Related Tauopathy

Expression of p62 in PART was assessed in various brain regions (Figure 1), including the neocortex, thalamus, basal ganglia, hippocampus, brainstem, cerebellar dentate nucleus, and the cervical spinal cord. The neocortex includes seven parts: the central anterior gyrus, central posterior gyrus, temporal pole, frontal lobe, cingulate gyrus, the parietal and the occipital poles; the thalamus includes the hypothalamus and subthalamus; the basal ganglia includes the caudate nucleus, putamen, and the globus pallidus; the hippocampus includes the anterior hippocampus, middle hippocampus, and the posterior hippocampus; the brainstem includes the substantia nigra, colliculus inferior/pons (CIP), locus coeruleus/pons (LCP), and the medulla oblongata.

FIGURE 1
www.frontiersin.org

Figure 1. Immunohistochemistry of the p62 + in PART brain area. (A) Gyrus pre/postcentralis, (B) Hypothalamus, (C) Caudatum/Putamen, (D) Amygdala, (E) Colliculus inferior/pons, (F) Cerebellum, (G) Spinal cord cervical, (H) Granule cell layer of hippocampus dentate gyrus (DG), (I) CA4, (J) CA1, (K) Subiculum (Sub), (L) Entorhinal Cortex(EC). Arrow point positive cell. Scale bars, (A–G) 20 μm, (H–L) 100 μm.

Overall, the positive rate of p62 increased with the increase of Braak NFT stage (Figure 2A) (*p < 0.05; **p < 0.01). Among the four groups, the positive rate of p62 was low in PART with Braak NFT stage I/II, moderate in PART with Braak NFT stage I _I _I/IV, which was close to that of pre-AD, and high in AD with Braak stage V/VI. There were statistical differences between PART with Braak NFT stage I/II and PART Braak NFT stage I _I _I/IV as well as between PART with Braak NFT stage I/II and AD, and between pre-AD and AD. Among the four groups of PART, positive expression of p62 increased with the increase of Braak NFT stage, but there was no statistical difference between these PART groups (Figure 2B).

FIGURE 2
www.frontiersin.org

Figure 2. (A) Quantification of the positive rate of p62 cell in PART and AD. Values represent mean ± SEM (n = 9–16 subjects; **p < 0.01; Kruskal–Wallis H test). (B) Quantification of the positive rate of p62 + cell in PART. Values represent mean ± SEM (n = 5–9 subjects; Kruskal–Wallis H test). (C) Quantification of the positive rate of p62 + cell in PART in brain areas. Values represent mean ± SEM (n = 0–9; *p < 0.05; Kruskal–Wallis H test). (D) Quantification of the positive rate of p62 + cell in PART in brain areas. Values represent mean ± SEM (n = 0–9; **p < 0.01; Kruskal-Wallis H test). (E) Quantification of the positive rate of p62 + cell in AD in brain areas. Values represent mean ± SEM (n = 0–9; *p < 0.05; Kruskal–Wallis H test). (F) Quantification of the positive rate of p62 + cell in AD braak stage III-IV in brain areas. Values represent mean ± SEM (n = 8–18; *p < 0.05; Kruskal–Wallis H test). (G) Quantification of the positive rate of p62 + cell in hippocampus. Values represent mean ± SEM (n = 8–18; **p < 0.01; Mann–Whitney U test).

Three regions with the highest expression of p62 in PART were the hippocampus, brainstem, and the thalamus, and there were significant differences (p < 0.05) between the neocortex and the hippocampus as well as between the brainstem and the neocortex (Figure 2C). Expression of p62 in PART cases with NFT Braak stage I was highest in the hippocampus and the brainstem, and there was a statistical difference (p < 0.01) between these two regions and the neocortex, respectively. P62 was not detectable in the cerebellum in PART cases with NFT Braak stage I (Figure 2D).

TABLE 1
www.frontiersin.org

Table 1. p62-positive structures in primary age-related tauopathy (PART).

2. Expression Pattern of p62 in Alzheimer’s Disease

The distribution pattern of p62-IR in the same brain regions of 11 AD cases with different NFT Braak stages was examined and compared with that of PART cases (Table 2). The positive rate of p62 was higher in AD than in PART. The highest expression was found in the hippocampus, and there was a statistical difference between the hippocampus and the neocortex (Figure 2E) (p < 0.05). When pre-AD (Braak NFT stage I _I _I/IV) was analyzed separately, the highest expression was found in the hippocampus, followed by the brainstem. There were significant differences(p < 0.05) between these two regions and the neocortex, respectively (Figure 2F). The positive rates of p62 in PART and AD were ranked in the order of the hippocampus, brain stem, thalamus, basal ganglia, neocortex, cerebellum, and the spinal cord (Figures 2C,E).

TABLE 2
www.frontiersin.org

Table 2. p62-positive structures in Alzheimer’s disease (AD).

3. Expression Pattern of p62 in the Hippocampus

Five regions of the hippocampus were checked, including the CA1, CA4, dentate gyrus, subiculum, and the entorhinal cortex, and semi-quantitative analysis was performed. The number of p62 positive cases in the AD group was significantly larger than in the pre-AD group, PART with Braak NFT stage I/II and PART with Braak NFT stage I _I _I/IV groups (Figure 2G). In the hippocampus of AD cases, the highest expression of p62 was found in the subiculum, followed by the CA1, the entorhinal cortex, and the CA4. p62 was rarely observed in the dentate gyrus. There were statistical differences between CA4 and the other three regions, respectively (p < 0.01). In pre-AD, the highest expression of p62 was found in the subiculum, followed by the CA1, the entorhinal cortex, and the CA4. There were significant differences between the subiculum and the entorhinal cortex as well as between the subiculum and the CA4 (p < 0.01) (Figure 2G).

Discussion

It remains to be determined whether p62 plays a role in PART. When anatomical regions were taken into consideration, we found that expression of p62 increased with the increase in Braak stages in both PART and AD. We also showed that p62 expression was highest in the hippocampus and the brainstem, especially in the subiculum and the CA1 of the hippocampus, and detected early in these regions. p62 expression was positive in the neocortex at later stages. These data provide the first direct evidence showing p62 expression in different brain regions in PART.

Expression of p62 was detected in multiple brain regions in PART, pre-AD, and AD brains with the highest expression in the hippocampus, followed by the brainstem and the thalamus (Figure 3). We lack the p62 expression data in the normal brain. This is our limitation. Therefore, it is difficult to infer that p62 expression of AD or PART is higher than that of normal controls. Here we focus on p62 expression of AD and PART. There was a positive correlation between expression of p62 and Braak stages in both PART and AD. Within the hippocampus, the subiculum had the highest level of p62, followed by the CA1 and other regions. In our previous research, we demonstrated regional patterns of selective vulnerability as distinguishing features of PART and AD in functionally relevant structures of the hippocampus. In AD cases, tau pathology was high in both CA1 and subiculum, followed by CA2/3, entorhinal cortex (EC), CA4, and dentate gyrus (DG). In PART, the severity of tau pathology in CA1 and subiculum was high, followed by EC, CA2/3, CA4, and DG (Zhang et al., 2020). These suggest that the subiculum of the hippocampus might be the earliest brain region that presents with AD pathology.

FIGURE 3
www.frontiersin.org

Figure 3. The schematic about the location of the tissues required for p62 staining. (A) Lateral view of brain, all lobes. (B) Sagittal cut of brain. (C) Coronal cut at basal ganglia. (D) Coronal cut of brain at hippocampus. Granule cell layer of dentate gyrus (DG), Subiculum (Sub), Entorhinal Cortex (EC).

TABLE 3
www.frontiersin.org

Table 3. The quantitative p62 score in hippocampus.

Six stages of p62 expression are proposed from the present study. At stage I, p62 immunoreactive inclusions are present in neurons of the hippocampus. At stage II, p62 expression expands to involve neurons of the brainstem. At stage I _I _I, the pathological inclusions spread to the thalamus. At stage IV, the basal ganglia show p62 expression. At stage V, the neocortical neurons show p62 positive inclusions. At Stage VI, the cerebellum shows p62 expression. These stages might represent disease progression in a sense. It has a pathological meaning in the clinical practice.

The multifunctional protein p62 is associated with neuropathological inclusions in various neurodegenerative disorders, including frontotemporal lobar degeneration (FTLD), amyotrophic lateral sclerosis and AD (Ma et al., 2019; Deng et al., 2020). Strong evidence has shown that p62 immunoreactivity is associated with neurofibrillary tangles formation (Kuusisto et al., 2002). In general, p62 protein is not a specific marker. Immunohistochemistry against p62 shows potential to facilitate neuropathological diagnosis by displaying a general pattern of stained inclusions. The use of p62 immunodetection was recommended as a sensitive diagnostic screening test for assessment of degenerative dementias, including FTLD-U (Pikkarainen et al., 2008).

Previous studies have shown that p62 participates in the formation of NFT in AD. p62 could be an early marker for neurodegenerative changes in PART, AD, and other neurodegenerative diseases (Piras et al., 2016). In the present study, we found that p62 was positive at the early stage of PART with Braak stage I/II and its expression was positively correlated with the NFT Braak stage. This supports that p62 might be involved in the aggregation of NFT at the early stage of AD (Zhang et al., 2020). In this study, p62 first appeared in the hippocampus. This is different from results demonstrating that LC is the earliest site of NFT formation (Braak et al., 2011; Attems et al., 2012).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by family members of brain donors. This research was approved by the Medical Ethics Committee of Zhejiang University School of Medicine. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

XW, HLu, JW, and BS prepared the tissue for immunohistochemistry and performed statistical analysis. KZ and LZ contributed to the statistical assessment and data processing. KZ designed the study, analyzed the results, and wrote the first version of the manuscript which was circulated among all the contributors for comments and suggestions. HLi contributed to the final version of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the National Science Foundation China (81971184) (to KZ), (91632109) (to Zhong and KZ), by the Zhejiang Provincial Natural Science Foundation (LY16H090013) (to KZ), and by the Zhejiang Medical and Health Science and Technology Plan Project (WKJ2013-2-009) (to KZ).

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.

Acknowledgments

We wish to thank the families of the patients who donated their brains to China Brain Bank in School of Medicine, Zhejiang University, to allow the completion of this study.

References

Attems, J., Thal, D. R., and Jellinger, K. A. (2012). The relationship between subcortical tau pathology and Alzheimer’s disease. Biochem. Soc. Trans. 40, 711–715. doi: 10.1042/BST20120034

PubMed Abstract | CrossRef Full Text | Google Scholar

Braak, H., Alafuzoff, I., Arzberger, T., Kretzschmar, H., and Del Tredici, K. (2006). Staging of Alzheimer disease-associated neurofibrillary pathology using paraffin sections and immunocytochemistry. Acta Neuropathol. 112, 389–404. doi: 10.1007/s00401-006-0127-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Braak, H., Thal, D. R., Ghebremedhin, E., and Del Tredici, K. (2011). Stages of the pathologic process in Alzheimer disease: age categories from 1 to 100 years. J. Neuropathol. Exp. Neurol. 70, 960–969. doi: 10.1097/NEN.0b013e318232a379

PubMed Abstract | CrossRef Full Text | Google Scholar

Crary, J. F., Trojanowski, J. Q., Schneider, J. A., Abisambra, J. F., Abner, E. L., Alafuzoff, I., et al. (2014). Primary age-related tauopathy (PART): a common pathology associated with human aging. Acta Neuropathol. 128, 755–766. doi: 10.1007/s00401-014-1349-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Z., Lim, J., Wang, Q., Purtell, K., Wu, S., Palomo, G. M., et al. (2020). ALS-FTLD-linked mutations of SQSTM1/p62 disrupt selective autophagy and NFE2L2/NRF2 anti-oxidative stress pathway. Autophagy 16, 917–931. doi: 10.1080/15548627.2019.1644076

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Y., Wooten, M. C., and Wooten, M. W. (2009). Oxidative damage to the promoter region of SQSTM1/p62 is common to neurodegenerative disease. Neurobiol. Dis. 35, 302–310. doi: 10.1016/j.nbd.2009.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Duyckaerts, C., Braak, H., Brion, J. P., Buée, L., Del Tredici, K., Goedert, M., et al. (2015). PART is part of Alzheimer disease. Acta Neuropathol. 129, 749–756. doi: 10.1007/s00401-015-1390-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Jellinger, K. A., Alafuzoff, I., Attems, J., Beach, T. G., Cairns, N. J., Crary, J. F., et al. (2015). PART, a distinct tauopathy, different from classical sporadic Alzheimer disease. Acta Neuropathol. 129, 757–762. doi: 10.1007/s00401-015-1407-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuusisto, E., Kauppinen, T., and Alafuzoff, I. (2008). Use of p62/SQSTM1 antibodies for neuropathological diagnosis. Neuropathol. Appl. Neurobiol. 34, 169–180. doi: 10.1111/j.1365-2990.2007.00884.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuusisto, E., Salminen, A., and Alafuzoff, I. (2002). Early accumulation of p62 in neurofibrillary tangles in Alzheimer’s disease: possible role in tangle formation. Neuropathol. Appl. Neurobiol. 28, 228–237. doi: 10.1046/j.1365-2990.2002.00394.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Dai, C., Fan, Y., Guo, B., Ren, K., Sun, T., et al. (2017). From autophagy to mitophagy: the roles of P62 in neurodegenerative diseases. J. Bioenerg. Biomembr. 49, 413–422. doi: 10.1007/s10863-017-9727-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, S., Attarwala, I. Y., and Xie, X. Q. (2019). SQSTM1/p62: a potential target for neurodegenerative disease. ACS Chem. Neurosci. 10, 2094–2114. doi: 10.1021/acschemneuro.8b00516

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirra, S. S., Heyman, A., McKeel, D., Sumi, S. M., Crain, B. J., Brownlee, L. M., et al. (1991). The consortium to establish a registry for Alzheimer’s disease (CERAD). Part II. Standardization of the neuropathologic assessment of Alzheimer’s disease. Neurology 41, 479–486. doi: 10.1212/wnl.41.4.479

PubMed Abstract | CrossRef Full Text | Google Scholar

Pikkarainen, M., Hartikainen, P., and Alafuzoff, I. (2008). Neuropathologic features of frontotemporal lobar degeneration with ubiquitin-positive inclusions visualized with ubiquitin-binding protein p62 immunohistochemistry. J. Neuropathol. Exp. Neurol. 67, 280–298. doi: 10.1097/NEN.0b013e31816a1da2

PubMed Abstract | CrossRef Full Text | Google Scholar

Piras, A., Collin, L., Grüninger, F., Graff, C., and Rönnbäck, A. (2016). Autophagic and lysosomal defects in human tauopathies: analysis of post-mortem brain from patients with familial Alzheimer disease, corticobasal degeneration and progressive supranuclear palsy. Acta Neuropathol. Commun. 4:22. doi: 10.1186/s40478-016-0292-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Salminen, A., Kaarniranta, K., Haapasalo, A., Hiltunen, M., Soininen, H., and Alafuzoff, I. (2012). Emerging role of p62/sequestosome-1 in the pathogenesis of Alzheimer’s disease. Prog. Neurobiol. 96, 87–95.

PubMed Abstract | Google Scholar

Zhang, L., Jiang, Y., Zhu, J., Liang, H., He, X., Qian, J., et al. (2020). Quantitative assessment of hippocampal tau pathology in AD and PART. J. Mol. Neurosci. 70, 1808–1811. doi: 10.1007/s12031-020-01573-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: p62, primary age-related tauopathy, hippocampus, Braak NFT stage, Alzheimier’s disease

Citation: Wang X, Zhang L, Lu H, Wu J, Liang H, Sun B and Zhu K (2022) Expression Pattern of p62 in Primary Age-Related Tauopathy: Staging of p62 in PART. Front. Aging Neurosci. 14:793353. doi: 10.3389/fnagi.2022.793353

Received: 12 October 2021; Accepted: 28 February 2022;
Published: 25 April 2022.

Edited by:

Shruti Arya, University of California, San Francisco, United States

Reviewed by:

Manuela Ceccarelli, National Research Council (CNR), Italy
Nina Schultz, Lund University, Sweden

Copyright © 2022 Wang, Zhang, Lu, Wu, Liang, Sun and Zhu. 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: Keqing Zhu, zhukeqing@zju.edu.cn

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.