An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves

The viscoelastic behavior of soft materials, especially cells and tissues, has been extensively investigated due to its importance in many biological and physiological processes that take place during development and even disease.*
Many techniques are used to quantify the mechanical properties of cells, among them micropipette aspiration, optical stretching, deformability cytometry and atomic force microscopy (AFM).*

The AFM, in particular, is still nowadays one of the most popular methods due to its conformity with various material types and geometries and the rather simple analysis process of the material properties.*

For a typical AFM indentation measurement, an AFM cantilever, with a distinct AFM tip shape, moves toward the sample with a predefined velocity and indents it until a prescribed force is reached. The AFM cantilever then moves upwards while detaching from the sample. The deflection and displacement signals of the AFM cantilever are processed further to extract the mechanical properties of the sample. Generally, a Hertzian model is fitted to the approach part of the force-indentation curves to quantify the apparent Young’s modulus.*

When applying the Hertzian model, few assumptions need to be considered, such as the material being homogeneous, isotropic, and linearly elastic. *

Cells and tissues, however, show not only elastic but also viscous behavior that is evident from the hysteresis between the approach and retraction segments of the force-indentation curve. Consequently, assessing this viscoelastic behavior is imperative for understanding the complex nature of biological matter.*

A number of studies utilized AFM to measure the viscoelastic properties of cells in both time and frequency domains.*

Ideally, to investigate the whole range of the viscoelastic behavior one needs to probe the material for a long time and observe its response or apply oscillatory signals and evaluate its phase lag. These approaches require the user to alter the probing method and add several steps to account for the time-dependent drift or the effect of the hydrodynamic drag of the surrounding medium. On top of that, in many of studies, the biological materials were probed with a linear approach followed by immediate retraction. The force-indentation curves from these studies were used to evaluate the apparent elastic modulus of the probed material using the standard Hertzian model. However, additional information concerning energy dissipation can still be extracted from the same curves to evaluate the viscoelasticity of the material.*

In the article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves”, Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland propose a new fitting model to extract the viscoelastic properties of soft materials from AFM force-indentation curves. *

To construct the explicit relation of force and indentation, the authors first use a generalization of Maxwell and Kelvin-Voigt models to describe soft materials, and numerically simulate the indentation of such material with a spherical indenter. *

Shada Abuhattum et al. show that the proposed Kelvin-Voigt-Maxwell (KVM) model adequately captures the force-indentation curves of materials having different mechanical characteristics. *

Based on the simulation results, Shada Abuhattum et al. further propose an explicit force-indentation relation to be fitted to the force-indentation curves. This explicit relation simplifies the association of the mechanical properties with physically meaningful components and processes.
Finally, the authors apply the fitting model to a number of samples, including poroelastic and viscoelastic hydrogels as well as HeLa cells in two different cell cycle phases, interphase and mitotic. *

Shada Abuhattum et al. demonstrate that the distinct nature of the hydrogels, arising from the different crosslinking mechanisms, can be described with the fitting model. For the HeLa cells, the mitotic cells had a higher apparent elasticity and a lower apparent viscosity, implying a stiffer actin cortex and a diluted cytoplasm protein concentration, when compared with interphase cells.*

Their findings demonstrate that the proposed model can reliably extract viscoelastic properties from conventional force-indentation curves. Moreover, the model is able to assess the contribution of the different elastic and viscous elements, and thus allows a direct comparison between the viscoelastic nature of different materials.*

AFM measurements were preformed using a commercially available Atomic Force Microscope. To indent the samples, NanoWorld Pyrex-Nitride tipless AFM cantilevers PNP-TR-TL with a nominal spring constant of 0.08 mN/m were modified by gluing 5 μm diameter polystyrene beads to the underside of the AFM cantilevers using two component glue.*

The AFM cantilevers were calibrated prior to each experiment using the thermal noise method and their accurate spring constant ranged between 0.047-0.059 mN/m. For PAAm and agarose hydrogels, the AFM cantilever was lowered with a constant velocity (5, 10, or 15 μm/s) toward the surface of the sample until a force of 2 nN for agarose and 4 nN for PAAm was reached. These force set points accounted for an indentation in the range of 0.5–1 μm. For HeLa cells, the AFM cantilever was lowered with a constant velocity of 2 μm/s and the cells were indented until a force of 2 nN was reached, which accounted for an indentation depth in the range of 0.5–1.5 μm.*

Graphical abstract for the article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves” by Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland consisting of 4 squares. showing a symbol for numerical simulations in the top left square, an arrow points to the bottom left square showing a graph and a formula as symbols for fitting algorithm a further arrow points to the bottom right square symbolizing the extraction of viscoelastic properties. The pictures in this square show on the left a drawing of the end of a tipless AFM cantilever on which a sphere is glued pressing on a cell, on the right of this picture there is another picture showing the end of a tipless AFM cantilever on which a sphere is glued pressing on a sphere or bead, underneath a graph symbolizing the mechanical properties of hydrogels is shown. Above this square on the top right a graph with a symbol for the mechanical behavior of the indented material is shown.NanoWorld Pyrex-Nitride tipless AFM cantilevers PNP-TR-TL with a nominal spring constant of 0.08 mN/m were modified by gluing 5 μm diameter polystyrene beads to the underside of the AFM cantilevers using two component glue were used for the atomic force microscopy indentation measurements described in the cited article.
Graphical abstract for the article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves” by Shada Abuhattum at al. 2022. NanoWorld Pyrex-Nitride tipless AFM cantilevers PNP-TR-TL with a nominal spring constant of 0.08 mN/m were modified by gluing 5 μm diameter polystyrene beads to the underside of the AFM cantilevers using two component glue were used for the atomic force microscopy indentation measurements described in the cited article.

*Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland
An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves
iScience, Volume 25, ISSUE 4, 104016, April 15, 2022
DOI: https://doi.org/10.1016/j.isci.2022.104016

The article “An explicit model to extract viscoelastic properties of cells from AFM force-indentation curves” by Shada Abuhattum, Dominic Mokbel, Paul Müller, Despina Soteriou, Jochen Guck and Sebastian Aland is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/.

Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis

Cells communicate with their environments via the plasma membrane and various membrane proteins. Clathrin-mediated endocytosis (CME) plays a central role in such communication and proceeds with a series of multiprotein assembly, deformation of the plasma membrane, and production of a membrane vesicle that delivers extracellular signaling molecules into the cytoplasm.*

In the article “Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis”, Aiko Yoshida, Nobuaki Sakai, Yoshitsugu Uekusa, Yuka Imaoka, Yoshitsuna Itagaki, Yuki Suzuki and Shige H. Yoshimura describe how they utilized their home-built correlative imaging system comprising high-speed atomic force microscopy (HS-AFM) and confocal fluorescence microscopy to simultaneously image morphological changes of the plasma membrane and protein localization during CME in a living cell.*

Overlaying AFM and fluorescence images revealed the dynamics of protein assembly and concomitant morphological changes of the plasma membrane with high spatial resolution. In particular, the authors elucidate the role of actin in the closing step of CME.*

The results revealed a tight correlation between the size of the pit and the amount of clathrin assembled. Actin dynamics play multiple roles in the assembly, maturation, and closing phases of the process, and affects membrane morphology, suggesting a close relationship between endocytosis and dynamic events at the cell cortex. Knock down of dynamin also affected the closing motion of the pit and showed functional correlation with actin.*

An AFM tip-scan–type HS-AFM unit combined with an inverted fluorescence/optical microscope equipped with a phase contrast system and a confocal unit was used for this study.*

The modulation method was set to phase modulation mode to detect AFM tip–sample interactions. A customized NanoWorld Ultra-Short AFM cantilever with an electron beam–deposited sharp AFM tip with a spring constant of 0.1 N m−1 (USC-F0.8-k0.1-T12) was used. *

All observations were performed at 28 °C. The AFM tip was aligned with confocal views as described in the Results section of the article. The images from the confocal microscope and AFM were simultaneously acquired at a scanning rate of 10 s/frame. The captured sequential images were overlaid by using AviUTL (http://spring-fragrance.mints.ne.jp/aviutl/) based on the AFM tip position.
The fluorescence intensity was quantified by Image J software (http://rsbweb.nih.gov/ij/). *

Fig 1 from Aiko Yoshida et al 2018 “Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis” :Aligning the confocal image and the AFM image. (A) Schematic illustration of the sample stage. A cross-shaped movable XY-stage (orange) is mounted on the base plate (light green) of the inverted optical microscope (IX83) via a stage guide (gray) equipped at each of the 4 ends of the cross. A 3-point support plate (purple) for mounting the AFM scanner unit is fixed on the base plate with a configuration that does not hinder the sliding of the XY-stage along the x-axis and y-axis. These setups allow the sample stage to move independently of the AFM unit and the optical axis. (B) Side view of the HS-AFM unit mounted on the stage illustrated in panel A. (C) Overlaying a confocal image and an AFM image. COS-7 cells expressing EGFP-CLCa were fixed with 5% paraformaldehyde and subjected to AFM (left) and CLSM (middle) imaging. The x-y position of the probe tip was determined as described in S1 Fig. Two images were overlaid (right) based on the x-y center position. Scale bar: 1 μm. Autofluorescence of the probe was much weaker than clathrin spot and could not be detected during the fast scanning. (D) AFM images of CCP on the cytoplasmic surface of the plasma membrane. COS-7 cells were “unroofed” by mild sonication as described in Materials and methods and then fixed with glutaraldehyde. Scale bar: 0.1 μm. AFM, atomic force microscopy; CCP, clathrin-coated pit; CLSM, confocal laser scanning microscopy; COS-7, CV-1 in origin with SV40 gene line 7; EGFP, enhanced green fluorescent protein; EGFP-CLCa, EGFP-fused clathrin light chain a; HS-AFM, high-speed AFM. https://doi.org/10.1371/journal.pbio.2004786.g001 customized NanoWorld Ultra-Short AFM cantilevers with an electron beam–deposited sharp AFM tip with a spring constant of 0.1 N m−1 (USC-F0.8-k0.1-T12) were used
Fig 1 from Aiko Yoshida et al 2018 “Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis” :
Aligning the confocal image and the AFM image.
(A) Schematic illustration of the sample stage. A cross-shaped movable XY-stage (orange) is mounted on the base plate (light green) of the inverted optical microscope (IX83) via a stage guide (gray) equipped at each of the 4 ends of the cross. A 3-point support plate (purple) for mounting the AFM scanner unit is fixed on the base plate with a configuration that does not hinder the sliding of the XY-stage along the x-axis and y-axis. These setups allow the sample stage to move independently of the AFM unit and the optical axis. (B) Side view of the HS-AFM unit mounted on the stage illustrated in panel A. (C) Overlaying a confocal image and an AFM image. COS-7 cells expressing EGFP-CLCa were fixed with 5% paraformaldehyde and subjected to AFM (left) and CLSM (middle) imaging. The x-y position of the probe tip was determined as described in S1 Fig. Two images were overlaid (right) based on the x-y center position. Scale bar: 1 μm. Autofluorescence of the probe was much weaker than clathrin spot and could not be detected during the fast scanning. (D) AFM images of CCP on the cytoplasmic surface of the plasma membrane. COS-7 cells were “unroofed” by mild sonication as described in Materials and methods and then fixed with glutaraldehyde. Scale bar: 0.1 μm. AFM, atomic force microscopy; CCP, clathrin-coated pit; CLSM, confocal laser scanning microscopy; COS-7, CV-1 in origin with SV40 gene line 7; EGFP, enhanced green fluorescent protein; EGFP-CLCa, EGFP-fused clathrin light chain a; HS-AFM, high-speed AFM.
https://doi.org/10.1371/journal.pbio.2004786.g001

 

Supporting figure 1 from Aiko Yoshida et al 2018 “Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis”:1 Fig. Aligning confocal and AFM images. (A) Scanning electron microscopy (SEM) images of a cantilever equipped with an EBD tip with tilt angle of 12°. Scale bar, 5 μm. Note that the cantilever is held on the AFM head unit with a tilt angle of 102° (from the x-y plane) so that the relative tip–sample angle (θ) is 90°. This setup makes it possible to precisely determine the position of the AFM tip. Scale bar, 2 μm. (B) Determining the position of the AFM probe in a fluorescence image. While the AFM probe was attached on the glass surface without scanning, the autofluorescence signal of the probe was imaged by the confocal scanning unit. The observed fluorescence spot (arrowhead in the middle panel) is defined as an origin of the fluorescence image plane (x = 0, y = 0) and used to define the optical axis (left panel). The position of a fluorescence spot derived from EGFP-CLCa was determined on this axis. On the other hand, the scanning area of the AFM scanner covers the area of (−3, 2.25) (left top), (3, 2.25) (right top), (3, −2.25) (right bottom), and (−3, −2.25) (left bottom) (all right panel). By aligning the axis from both images, the x, y position of the AFM image and that of the confocal fluorescence image could be merged. AFM, atomic force microscopy; EBD, electron beam–deposited; EGFP, enhanced green fluorescent protein; EGFP-CLCa, EGFP-fused clathrin light chain a. https://doi.org/10.1371/journal.pbio.2004786.s001 customized NanoWorld Ultra-Short AFM cantilevers with an electron beam–deposited sharp AFM tip with a spring constant of 0.1 N m−1 (USC-F0.8-k0.1-T12) were used
Supporting figure 1 from Aiko Yoshida et al 2018 “Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis”:
1 Fig. Aligning confocal and AFM images.
(A) Scanning electron microscopy (SEM) images of a cantilever equipped with an EBD tip with tilt angle of 12°. Scale bar, 5 μm. Note that the cantilever is held on the AFM head unit with a tilt angle of 102° (from the x-y plane) so that the relative tip–sample angle (θ) is 90°. This setup makes it possible to precisely determine the position of the AFM tip. Scale bar, 2 μm. (B) Determining the position of the AFM probe in a fluorescence image. While the AFM probe was attached on the glass surface without scanning, the autofluorescence signal of the probe was imaged by the confocal scanning unit. The observed fluorescence spot (arrowhead in the middle panel) is defined as an origin of the fluorescence image plane (x = 0, y = 0) and used to define the optical axis (left panel). The position of a fluorescence spot derived from EGFP-CLCa was determined on this axis. On the other hand, the scanning area of the AFM scanner covers the area of (−3, 2.25) (left top), (3, 2.25) (right top), (3, −2.25) (right bottom), and (−3, −2.25) (left bottom) (all right panel). By aligning the axis from both images, the x, y position of the AFM image and that of the confocal fluorescence image could be merged. AFM, atomic force microscopy; EBD, electron beam–deposited; EGFP, enhanced green fluorescent protein; EGFP-CLCa, EGFP-fused clathrin light chain a.
https://doi.org/10.1371/journal.pbio.2004786.s001

*Aiko Yoshida, Nobuaki Sakai, Yoshitsugu Uekusa, Yuka Imaoka, Yoshitsuna Itagaki, Yuki Suzuki and     Shige H. Yoshimura
Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis
PLoS Biol 16(5) (2018): e2004786
DOI: https://doi.org/10.1371/journal.pbio.2004786

The article “Morphological changes of plasma membrane and protein assembly during clathrin-mediated endocytosis” by Aiko Yoshida, Nobuaki Sakai, Yoshitsugu Uekusa, Yuka Imaoka, Yoshitsuna Itagaki, Yuki Suzuki and Shige H. Yoshimura is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/.

Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer

Flux-closure structures, vortices/antivortices, skyrmions, and merons in oxides, metals and polymers represent non-trivial topologies in which a local polar/magnetic order undergoes quasi-continuous spatial variations in a host crystal lattice. These structures are now extensively studied due to emergent functionalities, but the application of electrical/mechanical fields has so far only served to destroy the polar topologies of interest. *

Topology created by quasi-continuous spatial variations of a local polarization direction represents an exotic state of matter, but field-driven manipulation has been hitherto limited to creation and destruction. *In the article “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer” Mengfan Guo, Erxiang Xu, Houbing Huang, Changqing Guo, Hetian Chen, Shulin Chen, Shan He, Le Zhou, Jing Ma, Zhonghui Shen, Ben Xu, Di Yi, Peng Gao, Ce-Wen Nan, Neil. D. Mathur and Yang Shen report that relatively small electric or mechanical fields can drive the non-volatile rotation of polar spirals in discretized microregions of the relaxor ferroelectric polymer poly(vinylidene fluoride-ran-trifluoroethylene).*

These polar spirals arise from the asymmetric Coulomb interaction between vertically aligned helical polymer chains, and can be rotated in-plane through various angles with robust retention. *

Given also that their manipulation of topological order can be detected via infrared absorption, Mengfan Guo et al.’s work suggests a new direction for the application of complex materials. *

Mengfan Guo et al. produced a 100-nm-thick monolayer of face-on lamellae with vertically aligned polymer chains by melt-recrystallizing spin-coated thin films of P(VDF-TrFE).

The resulting melt-recrystallized thin film of the relaxor ferroelectric polymer was characterized by the authors using a commercial atomic force microscope for in-plane piezo-response force microscopy (IP-PFM).

NanoWorld Platinum Iridium coated Arrow-CONTPt AFM probes (typical resonant frequency: 14 kHz, typical force constant: 0.2 N/m, typical AFM tip radius 25 nm) were used for the in-plane (IP) PFM tests and the PFM lithography tests.

For piezo-response force microscopy (PFM) imaging, Vector mode was used where AFM tips were modulated at around 240 kHz for IP imaging, with the AC voltage set at 2 V. The images obtained by Vector Mode were double checked by using dual AC resonance tracking (DART) mode and the patterns could be reproduced. *

For angle-resolved IP-PFM tests, the rotation of sample was controlled by a protractor. To ensure identical position was imaged after rotating the sample, the authors made cross-scratches as a mark on the sample surface in advance. This method was applied to locate the scanning position in other situations if Mengfan Guo et al. had to move the sample in between the scanning probe microscopic studies.

For electric-field-induced manipulations using PFM lithography, the DC voltage on AFM tip was previously edited in the software. The scan speed was set at 1.95 Hz and no AC voltage was applied during the scanning. The DC voltage was divided by film thickness (100 nm) to obtain the electric field value. And an electric field with downward direction is defined with a positive sign.

For stress-induced manipulations, the deflection value of the PFM cantilever, which is a signal from photodetector, was preset to control the stress/force applied onto the sample. The difference in deflection value between a pressed AFM cantilever and a free AFM cantilever reflects how hard the AFM tip and sample surface are pressed to each other.*

To obtain the force value F, Mengfan Guo et al. first calibrated the AFM tips by the thermal noise method, and obtain the inverse optical lever sensitivity (InvOLS) and the spring constant k of the AFM tips.

The authors also conducted a polarization analysis based on their PFM measurements. *

To obtain the nominal toroidal order evaluated by the local curvature, the obtained IP-PFM amplitude image was firstly divided into 33 × 33 arrays, and each region was then subjected to a recognition of potential domain walls and measurement of an averaged curvature radius. *

To obtain polarization maps, angle-resolved IP-PFM images were first aligned to correct spatial distortion in nanoscale measurement. Positions with specific morphological characteristics were selected as reference points to determine the coordinate. After the correction, improved angle-resolved IP-PFM phase images would be divided into 64 × 64 arrays for deriving polarization maps. *

Fig. 1 from Mengfan Guo et al. (2024) “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer”:Observation of a microregion containing an in-plane polar spiral. a Morphology of a melt-recrystallized thin film of the relaxor ferroelectric polymer. The scale bar is 2 μm. IP-PFM phase (b) and amplitude (c) images of the same area in a exhibiting concentric ring-shaped domains in curly stripe domains. d Distribution of domain wall curvatures in the same area in a–c evidencing nominal toroidal order. It is assumed that the local polarization is parallel to the nearest domain wall so that larger curvature (denoted red) reflects stronger toroidal order. IP-PFM phase images of identical concentric ring-shaped domains with the axis along vertical (e) and horizontal (f) measurement directions. The scale bar is 0.3 μm. The curl (g) and the divergence (h) of local polarization in the same area as e and f, revealing the polar spiral topology. i Schematic stereoscopic view of a CCW polar spiral, arrows represent regions of polarization. The red/blue arrows denote the polar source/sink that spirals in/out. The white arrows represent Néel rotation along the radial direction, as shown in more detail via the inset. NanoWorld Platinum Iridium coated Arrow-CONTPt AFM probes (typical resonant frequency: 14 kHz, typical force constant: 0.2 N/m, typical AFM tip radius 25 nm) were used for the in-plane (IP) PFM tests and the PFM lithography tests.
Fig. 1 from Mengfan Guo et al. (2024) “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer”:
Observation of a microregion containing an in-plane polar spiral.
a Morphology of a melt-recrystallized thin film of the relaxor ferroelectric polymer. The scale bar is 2 μm. IP-PFM phase (b) and amplitude (c) images of the same area in a exhibiting concentric ring-shaped domains in curly stripe domains. d Distribution of domain wall curvatures in the same area in a–c evidencing nominal toroidal order. It is assumed that the local polarization is parallel to the nearest domain wall so that larger curvature (denoted red) reflects stronger toroidal order. IP-PFM phase images of identical concentric ring-shaped domains with the axis along vertical (e) and horizontal (f) measurement directions. The scale bar is 0.3 μm. The curl (g) and the divergence (h) of local polarization in the same area as e and f, revealing the polar spiral topology. i Schematic stereoscopic view of a CCW polar spiral, arrows represent regions of polarization. The red/blue arrows denote the polar source/sink that spirals in/out. The white arrows represent Néel rotation along the radial direction, as shown in more detail via the inset.

*Mengfan Guo, Erxiang Xu, Houbing Huang, Changqing Guo, Hetian Chen, Shulin Chen, Shan He, Le Zhou, Jing Ma, Zhonghui Shen, Ben Xu, Di Yi, Peng Gao, Ce-Wen Nan, Neil. D. Mathur and Yang Shen
Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer
Nature Communications volume 15, Article number: 348 (2024)
DOI: https://doi.org/10.1038/s41467-023-44395-5

The article “Electrically and mechanically driven rotation of polar spirals in a relaxor ferroelectric polymer” by Mengfan Guo, Erxiang Xu, Houbing Huang, Changqing Guo, Hetian Chen, Shulin Chen, Shan He, Le Zhou, Jing Ma, Zhonghui Shen, Ben Xu, Di Yi, Peng Gao, Ce-Wen Nan, Neil. D. Mathur and Yang Shen is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third-party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/.