Open Access Article
Liangliang
Yang‡
ab,
Qi
Gao‡
c,
Lu
Ge
ab,
Qihui
Zhou
d,
Eliza M.
Warszawik
ab,
Reinier
Bron
ab,
King Wai Chiu
Lai
*c and
Patrick
van Rijn
*ab
aDepartment of Biomedical Engineering-FB40, University of Groningen, University Medical Center Groningen, Groningen, A. Deusinglaan 1, 9713 AV Groningen, The Netherlands. E-mail: p.van.rijn@umcg.nl
bW.J. Kolff Institute for Biomedical Engineering and Materials Science-FB41, University of Groningen, University Medical Center Groningen, Groningen, A. Deusinglaan 1, 9713 AV Groningen, The Netherlands
cDepartment of Biomedical Engineering, City University of Hong Kong, Hong Kong. E-mail: kinglai@cityu.edu.hk
dInstitute for Translational Medicine, Department of Stomatology, The Affiliated Hospital of Qingdao University, Qingdao University, Qingdao 266003, China
First published on 5th April 2020
Topography-driven alterations in cell morphology tremendously influence cell biological processes, particularly stem cell differentiation. Aligned topography is known to alter the cell shape, which we anticipated to also induce altered physical properties of the cell. Here, we show that topography has a significant influence on single cell stiffness of human bone marrow derived-Mesenchymal Stem Cells (hBM-MSCs) and the osteogenic differentiation of these. Aligned topographies were used to control the cell elongation, depicted as the cell aspect ratio (CAR). Intriguingly, an equal CAR elicited from different topographies, resulted in highly altered differentiation behavior and the underlying single cell mechanics was found to be critical. The cell behavior was found to be focal adhesion-mediated and induced stiffness alterations rather than just influencing the cell elongation. The effect was further corroborated by investigations of the transcriptional regulators YAP. Our study provides insight into how mechanical properties of the cell, which are stimulated by topography, modulate the osteogenesis of hBM-MSCs, which is beneficial for improving the understanding of interactions between stem cells and topography for developing applications of tissue engineering and regenerative medicine.
Several types of geometries, namely pits,21 tubes,22 pillars,23 and gratings24 have been previously investigated in terms of osteogenic differentiation. Facilitating osteogenic differentiation by means of non-biochemical cues would have great implications for bone tissue engineering and bone-inductive implant materials.25,26 The topography of biomaterials modulates the differentiation of stem cells via mechanotransduction. Previous studies demonstrate that topography induced acceleration of osteogenic differentiation is caused by focal adhesion, RhoA/Rock signaling pathway,25,27 and the regulation microRNA.28 Furthermore, Liu et al., Wang et al., and Niu et al. revealed that topography-induced N-cadherin via β-catenin signaling,29 integrin-linked kinase (ILK)/β-catenin pathway,30 and crosstalk between focal adhesion kinase/MAPK and ILK/β-catenin pathways31 up-regulated osteogenesis.
Previous studies indicate that cell elongation could be a factor that facilitates differentiation of individual hBM-MSCs,24,32–36 even without the external chemical induction factors.36 The cell elongation can be quantified by determining the cell aspect ratio (CAR, defined as the ratio of the length of the major axis to the length on the minor axis of a single cell) and is an important characterization of cell shape, which has significant influence on the cell fate.34,36 Osteogenic differentiation of MSCs with different CAR (varing between 1–16) was studied by Wagner and co-workers,24 Mrksich and co-workers,32 and Ding and co-workers34 who independently showed that there is a dependency between the CAR and osteogenic differentiation. For instance, Wagner et al.24 found that compared to Flat surface, stem cells grown on microgrooves with higher CAR showed enhanced osteogenic differentiation. Mrksich et al.32 showed that cells cultured in rectangles with increasing aspect ratio, display different osteogenesis profiles. Specifically, patterned cells cultured on rectangles having aspect ratios of 1
:
1, 3
:
2 and 4
:
1 exhibited that the degree of osteogenic differentiation increased with aspect ratio. However, Ding et al.36 studied the effect of CAR on lineage commitment of MSCs in growth medium without the chemical induction factors, and a non-monotonic change of osteogenic differentiation was found with optimal CAR of about 2. This demonstrates that CAR itself is an inherent factor to control stem cell differentiation, let alone with or without chemical induction factors. However, it remains elusive why cell elongation enhances the differentiation and whether or not the morphology is indeed the cause or that other factors are influenced by the topography.
Previous studies demonstrate that topography influences cell behavior by changing focal adhesion assembly, giving rise to variations in cytoskeletal organization and cell mechanical properties.37,38 However, few studies investigated the relationship between cell elongation, cell mechanical properties and osteogenic differentiation of stem cells. We hypothesized that topography-mediated osteogenic differentiation of hBM-MSCs is associated with single cell mechanical properties rather than the cell aspect ratio (CAR). The idea is that cell elongation not only associated with biological properties, for example, focal adhesion and mechanosensing, but also result in altered physical properties (cell stiffness) (Fig. 1), which is rarely investigated and compared to the biological response. To study the effect of surface topography on cell morphology (cell orientation, cell area and CAR), we prepared polydimethylsiloxane (PDMS) substrates with anisotropic wave-like topographies with different wavelengths and amplitudes to controllably affect and direct the CAR. The CAR in connection to single cell mechanical alterations was determined by atomic force microscopy (AFM) nanoindentation and displayed the modulation of osteogenic differentiation exemplifying how influencing mechanical components in cell development is crucial. Topography-mediated mechanosensing was further studied on protein and molecular level using immunostaining of focal adhesion and the transcriptional regulator YAP, respectively. Furthermore, the altering mechanical behavior was studied in combination with osteogenesis and it was found that altering mechanical behavior is of major importance. For the first time, we describe the relationship between CAR, degree of osteogenesis, and mechanical properties of hBM-MSCs during that process. We found that cell stiffness is another important factor for topography-driven CAR guiding osteogenic differentiation of stem cells.
:
1 or 15
:
1 by weight. The mixtures were vigorously stirred with a spatula, degassed under vacuum for 15 min, and deposited onto clean 12 × 12 cm squared polystyrene Petri dish. The poured PDMS was cured overnight in an oven at 70 °C.
| Substrate | Ratio of prepolymer and cross-linker | Operating pressure | Plasma time | Stretch percent (%) |
|---|---|---|---|---|
| W and A are the abbreviation of wavelength and amplitude, respectively, and the unit for W and A is μm. The different surface topographies are further reported as W0.5, W3, W10, and W27. | ||||
| W0.5/A0.05 | 10 : 1 |
14 Torr | 60 s | 30 |
| W3/A0.7 | 10 : 1 |
25 mTorr | 20 s | 30 |
| W10/A3.5 | 10 : 1 |
25 mTorr | 650 s | 20 |
| W27/A4.3 | 15 : 1 |
25 mTorr | 30 min | 10 |
:
1 by weight) was poured, followed by curing at 70 °C overnight. After that, the molds were removed and the freshly prepared PDMS substrates baring the imprint were additionally treated with air plasma at 500 mTorr for 1 min and used for cell seeding. The Flat substrate (control) was prepared and treated the same as those obtained from the imprinting process mentioned above (same ratio for elastomer base and cross-linker (10
:
1), same curing process and condition of air plasma), to make sure all the samples maintain the same surface chemistry and mechanical properties.
:
100, v/v), vinculin (clone hVin-1, Sigma, 1
:
100) or YAP (Santa Cruz Biotechnology, SC-101199, 1
:
500, v/v) for 1 h, followed by staining with Rhodamine Red™-X-labeled goat-anti-mouse antibody (Jackson Immunolab, 1
:
100, v/v) as the secondary antibody. The nucleus and actin were stained using DAPI and TRITC-phalloidin, respectively by incubation for 1 hour. Finally, the cells were imaged with TissueFaxs (TissueGnostics GmbH, Vienna, Austria) at 10× magnification. YAP staining and vinculin staining were observed using a LEICA TCS SP2 confocal laser scanning microscopy (CLSM) equipped with a 40× NA 0.80 water immersion objective. Additionally, image analysis of focal adhesion was done by an online Focal Adhesion Analysis Server,40 and Fiji software was used to measure cell orientation and elongation. Cell elongation was quantified as the aspect ratio between the cell length and width measured via fluorescent F-actin stained cells (at least 60 cells counted for each sample, also for cell area). Directionality analysis for cell orientation was conducted with Fiji using the Orientation J plug-in. Only aligned elements were taken into account for the directionality algorithm in an entire image, and at least 200 cells were taken into consideration for each substrate. For cell elongation and orientation, ≥6 images per sample were used and 3 independent samples were analyzed.
![]() | (1) |
For each cell, approximately 25 force–distance curves were collected. At each condition, at least 25 cells were tested. Data analysis was performed by Nanoscope v1.8.
The surface features after imprinting were characterized by AFM. As shown in Fig. 2A, the non-stretched PDMS results in a Flat surface, while using the preparation conditions as illustrated in Table 1, wave-like topographies are achieved of varied dimensions (wavelength (W; μm) and amplitude (A; μm)) of W0.5A0.05, W3A0.7, W10A3.5 and W27A4.3. The different surface topographies are further reported as Flat, W0.5, W3, W10, and W27. The wavelength and amplitude are dependent and both increase simultaneously as is shown in Fig. 2B.
![]() | ||
| Fig. 2 (A) AFM images taken from middle part of substrate (images from side of the substrate showed in Fig. S1†) and height profiles of the structured PDMS surfaces obtained after imprinting. (B) Wavelength, and amplitude of created wrinkled surface. Data are reported as mean ± standard deviation (SD) (n ≥ 30 wrinkles for each imprint, three independent imprints). Scale bar is 4 μm and applies to all images. | ||
To investigate the influence of the surface topography on cell morphology, hBM-MSCs were seeded on the different topographies and allowed to attach and spread for 24 hours. To visualize the cytoskeleton and nucleus, the cells were stained with phalloidin and DAPI, respectively. As shown in Fig. 3A and B (zoom-in), cell orientation was strongly influenced by the surface topography on which the cells reside. Cells grown on Flat were randomly oriented. Cells cultured on W0.5 promoted slight orientation of cells in the direction of the topography. The cells cultured on the larger topographies resulted in cells with much higher orientation.
Cell alignment, expressed as the percentage of cells that have their main axis within 10° from the direction of the topography,46,47 was quantified with Fiji. The degree of cell orientation (Fig. 3C) increased with the increasing topography dimensions up to W10 (98%), but decreased to 88% upon further increasing the topography dimensions for W27. The results indicate that larger wavelength and amplitude improve the cell orientation, but the trend is a non-monotonic change. Cell area is a well-accepted characterization of cell attachment and spreading, and it is well documented that cell extension and consequent spreading could give rise to altering levels of differentiation markers of stem cells.48Fig. 3D displays that cell area (2 dimensional) gradually decreased with the increase of wrinkle size. hBM-MSCs on Flat and W0.5 surface exhibited the largest average cell area (2300 μm2). In contrast, the cell area for W27 (1580 μm2) was significantly lower compared to Flat and W0.5 substrates. These results suggest that there is an significant influence of surface topography on cell alignment and spreading.
It was found that the CAR is also strongly influenced by surface topography. For cells grown on a Flat surface, the CAR reached 4.4 and further increased with the increasing of wrinkle size. For cells grown on W0.5, W3, and W10, the CAR was 7.0, 11.4, 19.0, respectively, and then decreased to 12.2 for cells cultured on W27 (Fig. 3E). Interestingly, cells grown on both W3 and W27 resulted in a similar CAR despite significant differences in the surface topography. However, cells seem to have a different appearance on both substrates and this is likely attributed to the calculation for cell aspect ratio, as it is defined as the ratio of the length of the major axis to the length on the minor axis of a single cell, so different cell appearance could have the same value of aspect ratio. The general cell elongation here far exceeds the reached elongation factors in other studies and more interestingly, it enables us to compare elongated cells with similar CAR stimulated by different topography allowing to identify if role of cell morphology specifically.
With the culture time increasing, the stiffness of the MSCs cultured on W10 and W27 were still less affected than MSCs on W0.5 and W3. After 14 days, a slight increase (from 3.7 kPa to 4.4 kPa) of Young's modulus and little decrease (from 6.4 kPa to 5.6 kPa) were observed for cells grown on Flat and W3, respectively, while there were no significant changes detected over time for those grown on the other topography substrates. Therefore, the Young's modulus of hBM-MSCs cultured on W0.5 and W3 in GM was significantly higher than for the cells grown on Flat, W10, and W27. These results demonstrate that the same CAR (W3 and W27) does not lead to the same stiffness alteration and that an increasing CAR does not evoke an increasing stiffness. MSCs on W0.5 and W3 have a similar stiffness and on W10 a lower stiffness after 1 day even though the CAR is increased going from W0.5 to W10. These results indicate that the topography induced CAR alterations (altered cell elongation) is most likely not the best predictor for the cell development.
The expression of FA in hBM-MSCs was assessed by immunofluorescence staining for vinculin and visualized by CLSM after 24 hours seeded onto the different substrates. As shown in Fig. 5A, major differences in focal adhesion number and morphology were observed. For substrates of W0.5 and W3, hBM-MSCs had more well defined dash like vinculin spots (typical regarded as mature focal adhesions). In contrast, cells grown on Flat, W10 and W27 showed dot-like (transient) vinculin spots, indicating that W0.5 and W3 could enhance the expression of vinculin. To better understand the focal adhesion formation on different substrates, the FA area per cell was quantitatively analyzed (Fig. 5B). As vinculin is a key FA molecule that controls FA maturation by interactions with other FA proteins such as talin and actin,57 the size of focal adhesion complexes is always characterized by vinculin staining.25 FA area for cells cultured on W0.5 (275 μm2) and W3 (260 μm2) are much larger than Flat, W10 and W27 (about 150 μm2). These results suggest that, compared to Flat and higher wavelength wrinkle substrate (W10 and W27), small wavelength wrinkle substrate (W0.5 and W3) facilitate the formation of FA complexes.
The expression of YAP in hBM-MSCs was assessed by immunofluorescence staining 24 hours after being seeded onto the different topographies. As shown in Fig. 6A, the topography had a substantial effect on the localization of YAP (cells with non-nuclear localization displayed in Fig. S2†). For cells cultured on Flat, W0.5 and W3 substrates, YAP was present in the cell nucleus for 52–60% of the cells (Fig. 6B). For cells cultured on W10 and W27, expression of YAP were predominantly cytoplasmic and the percentage of cells with nuclear positive YAP was only 27% and 30%, respectively. Although, all cells received the same mechanical dosing originating from an oxidized silicone surface, the only difference is the applied topographical stimulus. The percentage of nuclear YAP is conventionally higher for stiff surfaces but it clearly shows the difference in topography regulation.
An immediate difference in fluorescence intensity of immunostained ALP is observed in Fig. 7A and B. Comparing MSCs cultured under growth conditions and osteogenic conditions show that hBM-MSCs cultured in GM for 7 and 14 days do not express ALP, while cells in OM showed increased expression of ALP. The expressed ALP varies among the different topographies that have been applied and illustrates the topography dependency on which the MSCs were cultured whereby W10 and W27 inhibited ALP expression as compared to the Flat control while W0.5 and W3 enhanced the ALP expression.
These effects were quantified by assessing the fluorescence and the results are shown in Fig. 7C and D for the 7 days and 14 days differentiation, respectively. The fluorescence output is corrected for the cell number and it correlated well with the qualitative analysis. The results show that W0.5 (CAR = 7.0) and W3 (CAR = 11.4) significantly facilitate osteogenic differentiation of hBM-MSCs as illustrated by a 1.2 and 1.3 fold increase, respectively as compared to the Flat substrate. Topographies W10 (CAR = 18.96) and W27 (CAR = 12.23) substantially suppressed ALP expression and displayed a 0.6 and 0.3 fold reduction, respectively (Fig. 7C). The difference became even more striking after 14 days of differentiation where a 2.8 and 3.2 fold increase was observed for W0.5 and W3 and a 0.4 and 0.5 fold reduction for W10 and W27 (Fig. 7D). While previously, an increasing CAR was considered to positively affect osteogenic differentiation, our results show that there is a limit to what extent elongation is beneficial. Interestingly, while W3 (CAR = 11.4) and W27 (CAR = 12.2) displayed similar elongation, their differentiation behavior were significantly different, indicating that cell elongation itself is not a determining factor.
hBM-MSCs cultured on W0.5 and W3 in OM exhibited more mineralized calcium nodules than those on the Flat surface, while mineralization on W10 and W27 was rather minimal. To quantify the mineralization degree of hBM-MSC, the stained calcium deposits were de-stained, and the optical density (OD) of the extracted stains was measured at 540 nm. As shown in Fig. 8C, the highest OD540 was obtained for cells grown on W3 followed by those on W0.5, which indicates a higher level of osteogenesis. The OD540 was significantly lower for W10 and W27 while the OD540 obtained from the mineral deposition on Flat was in between W0.5/W3 and W10/W27. This further confirms that cells grown on the W0.5 (CAR = 7.0) and W3 (CAR = 11.4) substrates enhances osteogenic differentiation of hBM-MSC, while culturing on W10 (CAR = 19.0) and W27 (CAR = 12.2) inhibits osteogenesis. The mineralization results are in excellent agreement with the ALP expression determined on 7th and 14th day of differentiation. Overall, these results indicate that the CAR induced by topography regulates the fate of stem cells but not in the fashion as previously considered. A higher CAR does not result in enhanced osteogenic differentiation and inducing the same CAR using different topography does not lead to the same osteogenic differentiation capabilities.
Similar to the cells cultured in growth medium (Fig. 4B) after 1 day, hBM-MSCs grown on W0.5 and W3 displayed enhanced stiffness compared to Flat, W10 and W27. Intriguingly, the Young's modulus of cells grown on W3 in OM displayed the highest at the beginning and the fastest reduction in stiffness after 14 days of differentiation, from the initial value of ∼6.2 kPa on day one to ∼2.1 kPa on day 14. A similar trend was observed for MSCs cultured on W0.5. For cells grown on Flat, the Young's modulus slightly decreased after 14 days of cell culture from about 3.6 kPa to 2.7 kPa. The stiffness of the MSCs cultured on W10 and W27 did not display the same reduction in stiffness and thereby suggests that the mechanical alterations, reduction in stiffness, that is needed for proper differentiation is inhibited.
Previously, different groups studied the effect of CAR on differentiation of mesenchymal stem cells. Mrksich and co-workers32 reported that shapes leading to high contractility, an increase in CAR, increased the osteogenesis. On the other hand, according to Ding and coworkers,34 a CAR of about 2 was reported to be optimal, which is much lower than the CAR of around 7 as reported by Abagnale et al.24 These contradicting results indicate that the differentiation behavior is not directly associated with the cell shape but that other phenomena are affecting it. Based on our results (Fig. 7 and 8), we have shown that the degree of osteogenesis increases with the increasing CAR: 4.4, 7.0, 11.4. The relation between CAR and osteogenesis was non-linear, as the differentiation for cell with a CAR of 19.0 and 12.2 induced by the larger topographies, resulted in lower osteogenesis. Interestingly, despite the similar CAR for W3 and W27, the degree of osteogenesis of hBM-MSCs on those substrates was substantially different. This clearly illustrates that the somewhat contradictive results in previous studies may find its origin not in the cell shape/morphology but in another factor that directs the cellular behavior.
Despite W3 and W27 lead to similar CAR, the cells on W3 showed significantly higher stiffness at day 1 compared to W27 (Fig. 4). It is important to highlight that this is prior to differentiation induction rather than during differentiation, which means that cell stiffness is the cause for affecting stem cell differentiation rather than being an accompanying consequence of the process. More importantly, the differentiation on W3 is more efficient than on W27, which indicates that the higher cell stiffness at day 1 is another important factor compared with CAR. The higher cell mechanical properties on W3 is likely attributed to the more focal adhesions (Fig. 5), which leads to more percentage of YAP located into nuclear (Fig. 6), giving rise to higher degree of osteogenic differentiation. Our results are consistent with Dalby, Furukawa and Fu's study, which demonstrate that osteogenesis is mediated by enhanced focal adhesion,28 actin polymerization25 and YAP activation.64 A schematic representation of the mechanistic phenome at play is illustrated in Fig. 10 and describes the effect of similar CAR, stimulated by different wrinkle dimensions, on the osteogenic differentiation, which is at least partly mediated by the mechanical properties of hBM-MSCs and possibly connected to the focal adhesions providing the contact guidance information and YAP signal pathway. However, further investigation is necessary to fully elucidate focal adhesion, cytoskeleton organization, and components of YAP signaling cascade, such as the associated upstream regulators, in the regulation osteogenesis of hBM-MSCs stimulated by topography. In addition, our results illustrate that the absolute change of cell stiffness between initial and final differentiation is necessary and important for the successful osteogenic differentiation. This finding coincides well with Cho's study,65 who reported a Young's modulus for osteoblasts of about 1.7 kPa and also observed a significant decrease in the elastic modulus for hMSCs exposed to an osteogenic medium for 10 days.66 Therefore, based on the mechanical results, we not only show that the CAR is not the determining factor for osteogenesis, but that modulating the mechanical properties of the cells by means of topography is more important than the actual CAR itself for stem cell differentiation.
Several studies demonstrate that mechanical properties of the cell may be used as novel biological markers of cell phenotypes.67,68 Previous research has shown that the mechanical properties of single cells are strongly connected with their lineage.68 For example, osteoblasts exhibit stiffer properties than chondrocytes, which in turn is stiffer than adipocytes. Importantly, during the differentiation of hBM-MSCs into osteoblasts, the mechanics change accordingly. Previously, it was illustrated that there are striking different characteristics in actin organization between hBM-MSCs and osteoblasts. In hBM-MSCs, the actin is organized as thick bundles or stress fibers while thin dense microfilament meshwork are present in osteoblasts.65,66 The difference in actin organization was shown to be reflected in the cellular elasticity. The Young's modulus for hBM-MSCs is about two times higher than that for osteoblast (1.7 kPa).66 In our study, compared to control (Flat), MSCs cultured on W3 substrate showed substantial higher stiffness at the very beginning (day 1), and then exhibited fast and significant decrease of mechanical properties. In contrast, stem cells on W27 did not displayed the higher stiffness at day 1 nor the decrease during the whole differentiation process. Such detailed characterization of hMSC mechanics is helpful for the quickly and accurately distinguishing the successfully differentiation behavior or not, and could be used along or even replace the analysis of the biological or chemical processes (for example, ALP staining, Alizarin red staining, RT-PCR and western blotting), which is important for the rapidly growing field of stem-cell-based tissue engineering.
Mechano-sensing is often investigated by determining the activity of YAP, which has been correlated with mesenchymal stem cell fate decision via the interaction with key determinants of osteogenic (Runx2) or adipogenic (PPARγ) differentiation.69,70 High stiffness triggers YAP translocation to the nucleus and this process is accompanied by changes in cytoskeletal F-actin and focal adhesions.71,72 Our results showed that YAP localization and activity correlates with the cell differentiation behavior driven by topography size (Fig. 6). A higher percentage of nuclear localization of YAP was determined for cells on W0.5 and W3, which topographies facilitated osteogenic differentiation. However, for planar surface, even though the similar percentage of nuclear localization, the degree of osteogenesis was much lower when compared to W0.5 and W3. This disparity is likely attributable to two aspects. The higher nuclear localization of YAP is possibly due to the stiffness of PDMS material, which could significantly affect the localization of YAP. Previous reports shown that YAP was more localized in the cytoplasm of MSCs if cultured on soft hydrogel (0.7 kPa) whereas a stiff hydrogel (40 kPa) induces nuclear localization.72 In our study, the stiffness of PDMS used for culturing stem cells is within the range of MPa, orders of magnitude larger than kPa, so it is reasonable that cells grown on the planar surface possess high percentage of nuclear localization. hBM-MSCs cultured on the W10 and W27 displayed significantly lower nuclear YAP localization, which indicates that the topography is interfering with the mechano-sensing capabilities as the material properties in all cases is the same. Therefore, material stiffness promotes osteogenic differentiation via mechanotransduction of the material to the cell while topography promotes osteogenic differentiation via mechano-altering effects on the cell by the topography, which can be stimulated as well as inhibited by choosing the correct topography. This effect is not directly related to the cell elongation and suggests that different feature sizes act on different cellular mechanisms.
Although, here mechanical properties and cell morphology are carefully characterized, it has to be noted that the origin of the mechanical deviations is not known. Cytoskeletal arrangements, membrane tension, intracellular fluid pressure, location of larger organelles such as the nucleus below the cell membrane may all attribute to the altering properties. Also the projection of the cells are two dimensional, as is commonly done, but in terms of mechanical properties, volume could play a role as well and future investigations might provide more insights into the subcellular characteristics and the influence of these on cellular mechanical properties.
Footnotes |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d0bm00264j |
| ‡ These authors contributed equally. |
| This journal is © The Royal Society of Chemistry 2020 |