Abstract
Caspases have functions particularly in apoptosis and inflammation. Increasing evidence indicates novel roles of these proteases in cell differentiation, including those involved in osteogenesis. This investigation provides a complex screening of osteogenic markers affected by pan caspase inhibition in micromass cultures derived from mouse forelimbs. PCR Array analysis showed significant alterations in expression of 49 osteogenic genes after 7 days of inhibition. The largest change was a decrease in CD36 expression, which was confirmed at organ level by caspase inhibition in cultured mouse ulnae followed by CD36 immunohistochemical analysis. So far, available data point to osteogenic potential of pro-apoptotic caspases. Therefore, the expression of pro-apoptotic caspases (-3, -6, -7, -8, -9) within the growth plate of mouse forelimbs at the stage where the individual zones are clearly apparent was studied. Caspase-9 was reported in the growth plate for the first time as well as caspase-6 and -7 in the resting zone, caspase-7 in the proliferation, and caspase-6 and -8 in the ossification zone. For all caspases, there was a gradient increase in activation toward the ossification zone. The distribution of staining varied significantly from that of apoptotic cells, and thus, the results further support non-apoptotic participation of caspases in osteogenesis.
Keywords
Introduction
Caspases are a family of cysteine proteases with well-established roles in apoptosis and inflammation. Recent work has been uncovering their manifold non-apoptotic functions, including those in differentiation of various cell types, cell cycle control, metabolism, and autophagy.1–3
These non-apoptotic roles of caspases appear to be important also in osteogenesis.4–8 The majority of bones develop via endochondral ossification, where condensing mesenchymal cells differentiate into chondrocytes, which gradually proliferate, mature, and hypertrophy in the forming growth plate. 9 Chondrocytes create an extracellular matrix of cartilage that is later invaded by blood vessels, osteoblasts, osteoclasts, and bone marrow cells. Recent evidence shows that some of the chondrocytes undergo apoptosis, some divide but many of them become bone cells. 10 The endochondral ossification is controlled and influenced by plenty of factors, molecules, and hormones,9,11,12 including caspases. Caspases are thought to be involved in chondrocyte terminal differentiation and apoptotic death of osteoblasts.13,14 Osteogenic potential of caspases has so far been described for caspase-3,5,8 caspase-7, 7 and suggested for caspase-8. 4
The objective of this work was to provide a complex screening of the osteogenic potential of caspases related to endochondral ossification. Osteogenic PCR Arrays of 84 candidate genes was performed after pan caspase inhibition in micromass cultures derived from embryonic mouse forelimbs. The greatest change (decrease) in expression in inhibited cultures found for CD36 was corroborated by further immunohistochemical study in caspase-inhibited cultured mouse ulnae. As the so far available data point particularly to osteogenic potential of pro-apoptotic caspases, activation of caspase-8 and -9 as the dominant initiator caspases and caspase-3, -6, -7 as the trio of executioner caspases was investigated across zones of the mouse ulnae growth plate and correlated with distribution of apoptotic cells.
Materials and Methods
Animals, Samples
Mouse (Mus musculus var. alba) forelimbs at embryonic day (E)18 from the mice strain CD1 were collected and decalcified in buffered ethylenediaminetetraacetic acid (EDTA) after immediate fixation in 4% buffered paraformaldehyde. Mice were killed according to the experimental protocol related to the project GACR 16-18430S approved by the University of Veterinary and Pharmaceutical Sciences, Brno, Czech Republic. Decalcified samples were dehydrated in gradient series of ethanol, treated with xylene, and embedded in paraffin. Serial sections (5 µm) were prepared and split over several slides: hematoxylin–eosin (H&E) and Alcian blue staining, cleaved caspases, and TdT-mediated biotin-dUTP nick end labeling (TUNEL) assay.
Micromass Cultures
Micromass cultures (primary mesenchymal cell cultures) were prepared from E12 mouse forelimbs (as, for example, in Mello and Tuan, 15 and Chlastakova et al. 16 ). The forelimbs were disintegrated mechanically in Puck’s saline A buffer (PSA; 0.4 g of KCl/l, 8 g of NaCl/l, 0.35 g of NaHCO3/l, 1 g of glucose/l) with 10% fetal bovine serum (FBS; Invitrogen, Karlsruhe, Germany, 9:1) and enzymatically in 10 U/ml of Dispase II (Sigma-Aldrich; Darmstadt, Germany) in 10% FBS/PSA at 37C, 1–1.5 hr, 1000 rpm, and vortexed every 15 min. After the disintegration, differentiating medium (60% F12/40% Dulbecco’s modified Eagle’s medium [DMEM], 10% FBS, 50 µg/ml ascorbic acid, 10 mM β-glycerol phosphate) was added to neutralize the Dispase activity, and the cell suspension was centrifuged at 1000 × g/5 min. Sediment was resuspended in medium, run through a 40 µm cell strainer (Corning Life Sciences; Falcon) and centrifuged. The final concentration of the cell suspension was 2 × 107 cells/ml. Spots (10 µl) were applied on Nunclon 6-well Delta surface culture dishes (Thermo Fisher Scientific; Waltham, MA) and after 1 hr at 37C in a 5% CO2 incubator, 2 ml of medium was added into each plate. Pan caspase inhibitor FMK001 (R&D Systems; Minneapolis, MN) diluted in dimethyl sulfoxide (DMSO; Sigma-Aldrich, St. Louis, MI) to the final concentration of 100 µM was applied to the micromass medium after stabilization of the micromass cultures per the manufacturer’s recommendation and previous experiments. 17 Differentiating medium with DMSO (control group) and pan caspase inhibitor (experimental group) was changed every second day. After 7 days, cells were collected using 350 µl lysis buffer RLT (Qiagen; Hilden, Germany) with β-ME (Sigma-Aldrich) and prepared for RNA isolation.
Limb Organ Cultures
Limb organ cultures were prepared from dissected mouse E18 ulnae. These were placed on Millipore filters above a metal grid and cultured in differentiating medium with DMSO (control group, inhibitor vehicle) and pan caspases inhibitor (experimental group, 100 µM) at 37C in a 5% CO2 incubator; the medium was changed every day. After 7 days, ulnae were fixed in 4% buffered paraformaldehyde overnight and histologically processed for further analysis. They were dehydrated in gradient series of ethanol, treated with xylene, and embedded in paraffin. Serial sections (5 µm) were prepared and split over several slides: H&E and Alcian blue staining, CD36 antigen immunohistochemistry.
PCR Arrays
RNeasy Mini Kit (Qiagen) was used for RNA isolation, then mRNA was transcribed into cDNA using SuperScript VILO (Invitrogen; Carlsbad, CA), and PCR Arrays were applied to analyze 84 genes connected with osteogenesis (PAMM-026A-24; SA Biosciences, Frederick, MD). Data were statistically evaluated by PCR Array Data Analysis V4 (SA Biosciences). Statistical significance was determined as p<0.05. The threshold of fold change was established as ±1.5. Housekeeping genes used in this analysis were: Actb, B2m, Gapdh, Gusb and Hsp90ab1. Three independent biological repetitions were performed and analyzed.
Immunohistochemistry
Histological sections were deparaffinized in xylene, rehydrated in alcohol series, and pretreated with citrate buffer (5 min at 97C). Endogenous peroxidase activity was eliminated by 3% hydrogen peroxide in phosphate-buffered saline (PBS; 1 min at room temperature), and nonspecific staining was eliminated by incubation in goat serum (45 min, room temperature, ABC kit; Vectastain, Vector Laboratories). Primary antibodies (Table 1) were used overnight at 4C in a humidified chamber. The next day, biotinylated secondary antibody was applied (30 min, room temperature, ABC kit; Vectastain, Vector Laboratories) followed by streptavidin-peroxidase complex (30 min, room temperature, ABC kit; Vectastain, Vector Laboratories). Chromogen substrate 3,3′-diaminobenzidine tetrachloride (DAB, K3466; Dako, Glostrup, Denmark) was used for visualization of positive cells (brown). Hematoxylin was used for counterstaining to clearly distinguish negative cells in blue. Negative controls were performed by omitting the primary antibodies. Interdigital tissue of embryonic mouse limbs was used as a positive control in the case of caspase-3, -6, -7, -9, 18 caspase-8, 19 and hypertrophic chondrocytes in the case of CD36. 20
Information About Primary Antibodies Used in This Study.
TUNEL Assay
After rehydration and elimination of endogenous peroxidase activity by 3% hydrogen peroxide in PBS (10 min, room temperature), samples were pretreated with proteinase K (20 µg/ml, 15 min, room temperature; Merck Millipore, Billerica, MA). ApopTag Peroxidase In Situ Apoptosis Detection Kit was used (Merck Millipore), after equilibration buffer (15 min, room temperature), reaction mix was applied for 45 min at 37C in humidified chamber. After anti-digoxigenin-peroxidase reaction (30 min, room temperature), positive cells were visualized by the color reaction using the DAB substrate (K3466; Dako), and slides were counterstained with hematoxylin.
Results
Pan Caspase Inhibition Significantly Impacts Osteogenic Gene Expression Profile in Micromass Cultures
To examine the role of caspases in osteogenic signaling pathways, osteogenic markers were evaluated by PCR Arrays analysis of micromass cultures treated for 1 week with pan caspase inhibitor. Expression of 49 genes was significantly different in the micromass cultures after pan caspase inhibition, compared with non-inhibited controls (treated with inhibitor vehicle DMSO only), 26 of them fold change above 2 (Fig. 1) and 23 of them with fold change between 2 and 1.5 (Fig. 2). Statistical significance was determined using a t-test and distinguished by color in the figures: p<0.05 (blue color), p<0.01 (green color), p<0.001 (brown color), p<0.0001 (red color). Decreased gene expression was seen in the case of Alkaline phosphatase (AlpI), Annexin A5 (Anxa5), Bone morphogenetic proteins (BMPs), and their receptors (Bmp1, Bmp2, Bmp3, Bmp4, Bmp5, Bmp6, Bmpr1a, Bmpr1b, Bmpr2), CD 36 antigen, Cadherin 11 (Cdh11), Chordin (Chrd), Collagens (Col10a1, Col14a1, Col1a2, Col2a1, Col3a1, Col5a1), Cartilage oligomeric matrix protein (Comp), Cathepsin K (Ctsk), Distal-less homeobox 5 (Dlx5), Epidermal growth factor (Egf), Fibroblast growth factor receptors (Fgfr1, Fgfr2), GLI-Kruppel family member (Gli1), Intercellular adhesion molecule 1 (Icam1), Insulin-like growth factor 1 (Igf1), Integrin alpha 2 (Itga2), Integrin alpha M (Itgam), Matrix metalloproteinase 9 (Mmp9), Nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 (Nfkb1), Noggin (Nog), Phosphate regulating gene with homologies to endopeptidases on the X chromosome (Phex), Runt-related transcription factor 2 (Runx2), Serine (or cysteine) peptidase inhibitor, member 1 (Serpinh1), Smads (Smad2, Smad3, Smad5), Sclerostin (Sost), SRY-box containing gene 9 (Sox9), Sp7 transcription factor 7 (Sp7), Transforming growth factor beta 2 (Tgfb2), and Transforming growth factor beta receptor III (Tgfbr3). On the contrary, increased gene expression was observed in the case of Mmp10, Secreted phosphoprotein 1 (Spp1), and Tumor necrosis factor (ligand) superfamily, member 11 (Tnfsf11).

Pan caspase inhibitor altered expression of 49 osteogenic genes, 26 genes by the fold change above 2. Abbreviations: Spp1, Secreted phosphoprotein 1; Sp7, Sp7 transcription factor 7; Sox9, SRY-box containing gene 9; Sost, Sclerostin; Phex, Phosphate regulating gene with homologies to endopeptidases on the X chromosome; Nog, Noggin; MMP9 and 10, Matrix metalloproteinase 9 and 10; Itgam, Integrin alpha M; Itga2, Integrin alpha 2; Igf1, Insulin-like growth factor 1; Gli1, GLI-Kruppel family member; FgFr2, Fibroblast growth factor receptors; Egf, Epidermal growth factor; DIx5, Distal-less homeobox 5; Ctsk, Cathepsin K; Comp, Cartilage oligomeric matrix protein; Col2a1, Col14a1, Col10a1, Collagens; Cd36, CD 36 antigen; Bmpr1b, 6, 5, 3, Bone morphogenetic proteins, and their receptors; AlpI, Alkaline phosphatase. The p value is distinguished by the color: p<0.05 (blue color), p<0.01 (green color), p<0.001 (brown color), p<0.0001 (red color). *Expression of genes impacted by caspase-3, data published in our previous report by Adamova et al. 8 **The fold change of CD36 reduced 10 times.

Pan caspase inhibitor altered expression of 49 osteogenic genes, 23 genes by the fold change between 2 and 1.5. Abbreviations: Tnfsf11, Tumor necrosis factor (ligand) superfamily, member 11; Tgfbr2, Transforming growth factor beta 2; Tgfbr3, Transforming growth factor beta receptor III; Smad5, 4, 3, 2, Smads; Serpinh1, Serine (or cysteine) peptidase inhibitor, member 1; Runx2, Runt-related transcription factor 2; Nfkb1, Nuclear factor of kappa light polypeptide gene enhancer in B-cells 1; Icam1, Intercellular adhesion molecule 1; Fgfr1, Fibroblast growth factor receptors; Col5a1, Col3a1, Col1a2, Collagens; Chrd, Chordin; Cdh11, Cadherin 11; Bmpr2, 1a, 4, Bone morphogenetic proteins, and their receptors; Bmp2, 1, Bone morphogenetic proteins; Anxa5, Annexin A5. The p value is distinguished by the color: p<0.05 (blue color), p<0.01 (green color), p<0.001 (brown color), p<0.0001 (red color). *Expression of genes impacted by caspase-3, data published in our previous report by Adamova et al. 8
CD36 Is Expressed Throughout the Prenatal Development of Long Bones and Its Expression Is Significantly Decreased After Caspase Inhibition
As CD36 expression was the most prominent change after caspase inhibition, its expression on the protein level was examined in developing mouse forelimbs, with a focus on the ulnae. CD36 was detected in the resting chondrocytes, and was particularly abundant in hypertrophic chondrocytes at E15 (data not shown). E18 ulna was used (Fig. 3A) to evaluate CD36 in the growth plate (Fig. 3B). CD36 antigen was present in the resting (Fig. 3C), proliferating (Fig. 3D), hypertrophic (Fig. 3E), and ossification zones (Fig. 3F).

CD36 is expressed in the growth plate and significantly downregulated after caspase inhibition. Dissected ulna at E18 (A), hematoxylin–eosin and Alcian blue staining of the forelimb at E18 (B). Expression of CD36 in the resting (C), proliferating (D), hypertrophic (E), and ossification (F) zones of the growth plate. E18 ulna + 7 days of cultivation in control medium (G) and expression of CD36 in the ulna cultivated in control medium (H, I, J). E18 ulna + 7 days of cultivation in medium with inhibitor (K) and expression of CD36 in medium with inhibitor (L, M, N). Arrows point to positive cells (brown). Scale bar = 100 µm.
To examine the impact of caspase deficiency on CD36 expression in the intact bone, dissected ulnae were cultured for 7 days under the conditions of caspase inhibition. In non-inhibited control samples (Fig. 3G), CD36 is expressed throughout the growth plate, especially in the resting, proliferation, and hypertrophic zones (Fig. 3H, I, and J). By contrast, in the pan caspase inhibitor treated samples (Fig. 3K), the resting (Fig. 3L), proliferation (Fig. 3M), and hypertrophic (Fig. 3N) zones are mostly CD36 negative.
Pro-Apoptotic Caspases Are Activated in the Growth Plate During Endochondral Ossification, With a Gradient Increase in Activation Toward the Ossification Zone
To follow caspase activation in all major zones of the long bone development, mouse forelimbs at E18 were collected for examination. Resting, proliferation, hypertrophic, and ossification zones were clearly distinguishable at this stage. The resting zone of the cartilage has typical round-shaped chondrocytes containing large nuclei compared with the cytoplasm, in the proliferation zone, where round chondrocytes become flattened and create multicellular clusters. Hypertrophic zone is characterized by 5- to 10-fold increased size of chondrocytes and massive secretion of extracellular matrix and in the ossification zone, osteoblasts are differentiated and the matrix is becoming calcified (Fig. 4A).

Activated pro-apoptotic caspases are expressed throughout the growth plate of developing long bones. Hematoxylin–eosin and Alcian blue staining of the growth plate (A), immunohistochemical detection of individual caspases in resting, proliferation, hypertrophic, and ossification zones of the growth plate: caspase-3 (B, C, D, E), caspase-6 (F, G, H, I), caspase-7 (J, K, L, M), caspase-8 (N, O, P, Q), and caspase-9 (R, S, T, U). Arrows point to positive cells (brown). Scale bar = 50 µm.
Some expression of active caspase-3 was detected in the resting zone of the growth plate (Fig. 4B), and increasing activation was observed in the proliferation (Fig. 4C) and hypertrophic zones (Fig. 4D). Presence of caspase-3 was apparent mostly in the ossification zone (Fig. 4E). Active caspase-6 in the resting (Fig. 4F), proliferating (Fig. 4G), and hypertrophic (Fig. 4H) zones displayed low activation but active caspase-6 was visible particularly in the ossification zone (Fig. 4I). Active caspase-7 was prominent through the growth plate in the resting (Fig. 4J), proliferating (Fig. 4K), hypertrophic (Fig. 4L), as well as ossification zone (Fig. 4M). Similar results were obtained in the case of caspase-8, also activation of this caspase displayed an increasing trend from the resting (Fig. 4N), through proliferating (Fig. 4O), and hypertrophic (Fig. 4P), toward the ossification zone (Fig. 4Q). Some cells positive for caspase-9 were observed in the resting (Fig. 4R), proliferating (Fig. 4S), and hypertrophic zones (Fig. 4T), however, activation of caspase-9 was linked particularly with the ossification zone (Fig. 4U).
Pro-Apoptotic Caspases Are Present in Non-Apoptotic Cells of the Growth Plate as Demonstrated in the Ossification Zone
To evaluate the presence of pro-apoptotic caspases in non-apoptotic cells, the pattern of expression of individual active caspases was correlated with localization of apoptotic (TUNEL positive) cells. The ossification zone, where the expression of all caspases was most prominent, was used for examination. The distribution of staining in the ossification zone for individual caspases (Fig. 5A, C, E, G, and I) was different from that for TUNEL (Fig. 5B, D, F, H, and J), showing presence of pro-apoptotic caspases in non-apoptotic cells.

Pro-apoptotic caspase activation occurs predominantly in non-apoptotic cells. Caspase-3 (A), caspase-6 (C), caspase-7 (E), caspase-8 (G), and caspase-9 (I) expression in the ossification zone, apoptotic cells detected by TUNEL assay on serial sections (B, D, F, H, J). Arrows point to positive cells (brown). Scale bar = 100 µm. Abbreviation: TUNEL, TdT-mediated biotin-dUTP nick end labeling.
Discussion
This work provides novel data related to emerging non-apoptotic roles of caspases in osteogenesis. The general screening of osteogenic impact of caspases related to endochondral ossification identified novel genes whose expression was influenced by active caspases within the osteogenic networks. The evidence about new functions of caspases in endochondral ossification is further supplemented by the activation pattern of caspases-3, -6, -7, -8, and -9 in non-apoptotic cells of the developing mouse growth plate.
In previous studies dealing with caspases and osteogenesis, pan caspase inhibition was performed in limb cultures, and altered expression of Runx2, Alp, Col2a1, and Gdf5 was followed. 21 In caspase-inhibited human chondrocytes, apoptosis was prevented and delayed while chondrocytes remained functional.22,23 In our previous article, changes in gene expression (Bmp1, Bmp5, Bmp6, Col10a1, Col2a1, Comp, Egf, Fgfr2, Gli1, Igf1, Nog, Phex, Sox9, and Spp1) following inhibition of caspase-3 were reported, and corresponding alterations of these genes are obvious also after general caspase inhibition. 8
The most prominent change in gene expression upon caspase inhibition was observed in the case of the CD36 antigen. The apparent downregulation of CD36 demonstrated by PCR Arrays was confirmed by pan caspase inhibition in cultured ulnae, where the expression of CD36 was maintained in controls but was strongly decreased in treated samples. CD36 functions as a receptor for collagens,24,25 trombospondin, 26 and oxidized low-density lipoprotein. 27 CD36 is suggested to be involved in bone resorption by activation of c-src signaling in osteoclasts.28,29 CD36 was also presented as a chondrocyte hypertrophy marker, 20 and it participates as a key regulator in bone mineralization and osteoblast functions.30,31
After pan caspase inhibition, significantly altered gene expression was also detected in the case of proteins involved in the degradation of the extracellular matrix, the matrix proteinases (Mmp9, Mmp10, and Ctsk). Mmp10 is expressed throughout the growth plate and in osteoclasts. 32 Mmp9 is responsible for the invasion of osteoclasts into mineralized cartilage and apoptosis of hypertrophic chondrocytes, and is involved in angiogenesis. 33 Ctsk is abundantly expressed by osteoclasts. 34 Other extracellular matrix molecules (Itga2, Itgam) were significantly impacted as well, such as integrins participating in cell-cell and cell-extracellular matrix interactions. 35
Some BMPs, members of the TGFβ superfamily, were significantly altered by caspase inhibition, and have an important role in cartilage formation (chondrocyte proliferation) during endochondral ossification.36,37 The downregulation of Bmp3 and Bmpr1b receptor following caspase inhibition may indicate a role for caspases in chondrogenesis maintenance, as the upregulation of BMPs in mesenchymal stem cells has osteoinductive activity. 38 Moreover, our observations are in agreement with reduction of Alp activation after inhibition of caspase activity in BMP-4-treated cells. 4 Alkaline phosphatase has a role in mineralization of the extracellular matrix produced by chondrocytes. 39 The BMP antagonist, Sost, is known as a marker of osteocytes 40 involved in the inhibition of osteoblast activity and the reduction of differentiation of osteoprogenitors. 41
Collagen proteins are the main organic components of bone extracellular matrix. 42 The decreased expression of several types of collagens, especially Col14a1 after caspase inhibition indicates changes in the creation of connective tissue in treated micromass cultures. Several transcription factors involved in osteogenesis, particularly Sp7 (Osterix) and Dlx5, were downregulated in the treated cultures. Both molecules are necessary for osteoblast differentiation, and their expression is induced by Bmp2. 43
Osteogenic potential has been associated particularly with pro-apoptotic caspases.4,5,7,8 Our temporospatial investigation within the growth plate supplements previous knowledge (Table 2). In our study, the expression of caspase-6, -7, and -9 in the resting and ossification zones, caspase-7 in the proliferation zone, and caspase-9 in the proliferation and hypertrophic zones were demonstrated for the first time. Moreover, comparison of the activation of individual caspases within the ossification zone with the pattern of apoptosis clearly showed presence of these caspases in non-apoptotic cells. In osteogenesis, non-apoptotic activity of pro-apoptotic caspases was earlier demonstrated in the case of caspase-3 in vitro,4–6 caspase-7,7,44 and caspase-8 in vitro. 4 Non-apoptotic engagement of these caspases, particularly in proliferation and differentiation, is known also from other organs/tissues/cells, such as skin, brain, muscles, keratinocytes, or immune system cells (reviewed in Shalini et al. 3 ).
References Reporting About Activated Caspases in the Growth Plate of Long Bones Supplemented by Our Data.
Individual references are given for published works.
Data investigated and presented here for the first time.
Regarding engagement of examined caspases in apoptosis of hypertrophic chondrocytes, our results are consistent with those of previous research.50,51 A number of possible mechanisms have been posited for apoptotic versus non-apoptotic engagement of caspases, including subcellular localization,6,52,53 co-expressions with inhibitor of apoptosis proteins, 54 levels of activation of caspase-activated DNase, 55 sublethal concentrations of caspases, 56 and substrate specificity of different caspases determined by posttranslational modifications and interactions with other proteins.57–59 Apparently, the complex cellular context determines the consequences of caspase activation. 60
This investigation added missing data to pro-apoptotic caspase activation within the growth plate and, thus, completed a comprehensive overview of pro-apoptotic caspase pattern related to apoptotic and non-apoptotic events in the resting, proliferating, hypertrophic, and osteogenic zones of the growth plate. The increasing gradient of caspase-positive cells from resting to ossification zone pointed to non-apoptotic functions of caspases in osteogenic pathways. PCR Array-based analysis of general osteogenic potential of caspases indicated several genes within osteogenic molecular networks whose expression is likely to interfere with active caspases in a non-apoptotic way and, thus, provided a list of candidates for further research or targeted modulations related to chondrogenesis and osteogenesis.
Footnotes
Competing Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Author Contributions
EJ performed the experiments, analyzed the results, and drafted the manuscript, PB performed the experiments, and EM designed the study and revised the manuscript. All authors have read and approved the final manuscript.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by the Czech Science Foundation (GACR 16-18430S) and P.B. by the IGA UVPS Brno, Czech Republic (117/2017/FVL).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
