- Open Access
Micro-computed tomography of pulmonary fibrosis in mice induced by adenoviral gene transfer of biologically active transforming growth factor-β1
- Thomas Rodt1_1024, 2_1024Email author,
- Christian von Falck1_1024, 2_1024,
- Sabine Dettmer1_1024,
- Roman Halter2_1024,
- Regina Maus3_1024,
- Kjetil Ask4_1024,
- Martin Kolb4_1024,
- Jack Gauldie4_1024,
- Florian Länger5_1024,
- Ludwig Hoy6_1024,
- Tobias Welte7_1024,
- Michael Galanski1_1024,
- Ulrich A Maus3_1024Email author and
- Jürgen Borlak2_1024Email author
© The Author(s) 2010
- Received: 7 July 2010
- Accepted: 22 December 2010
- Published: 22 December 2010
Micro-computed tomography (micro-CT) is a novel tool for monitoring acute and chronic disease states in small laboratory animals. Its value for assessing progressive lung fibrosis in mice has not been reported so far. Here we examined the importance of in vivo micro-CT as non-invasive tool to assess progression of pulmonary fibrosis in mice over time.
Pulmonary fibrosis was induced in mice by intratracheal delivery of an adenoviral gene vector encoding biologically active TGF-ß1 (AdTGF-ß1). Respiratory gated and ungated micro-CT scans were performed at 1, 2, 3, and 4 weeks post pulmonary adenoviral gene or control vector delivery, and were then correlated with respective histopathology-based Ashcroft scoring of pulmonary fibrosis in mice. Visual assessment of image quality and consolidation was performed by 3 observers and a semi-automated quantification algorithm was applied to quantify aerated pulmonary volume as an inverse surrogate marker for pulmonary fibrosis.
We found a significant correlation between classical Ashcroft scoring and micro-CT assessment using both visual assessment and the semi-automated quantification algorithm. Pulmonary fibrosis could be clearly detected in micro-CT, image quality values were higher for respiratory gated exams, although differences were not significant. For assessment of fibrosis no significant difference between respiratory gated and ungated exams was observed.
Together, we show that micro-CT is a powerful tool to assess pulmonary fibrosis in mice, using both visual assessment and semi-automated quantification algorithms. These data may be important in view of pre-clinical pharmacologic interventions for the treatment of lung fibrosis in small laboratory animals.
- Idiopathic Pulmonary Fibrosis
- Pulmonary Fibrosis
- Lung Fibrosis
- Seed Point
- Image Quality Assessment
Pulmonary fibrosis is a severe, chronic lung disease, which is associated with high mortality rates, afflicting more than five million patients worldwide . Its most severe form, idiopathic pulmonary fibrosis (IPF) is typically associated with an average survival of 2-5 years from time of diagnosis [2, 3]. Depending on the underlying cause, progression of pulmonary fibrosis is associated with alveolar epithelial cell injury and hyperplasia, inflammatory leukocyte accumulations within the bronchoalveolar space and lung parenchymal tissue, as well as fibroblast hyperplasia and exuberant deposition of extracellular matrix (ECM) components such as collagen and fibronectin in the lungs, eventually leading to formation of honey comb cysts and scar tissue [2, 4]. Specifically in IPF, immunosuppressive and anti-inflammatory therapies have proved to be ineffective in preventing the fatal consequences of lung tissue remodelling characterizing pulmonary fibrosis.
Over the past few years, it has become evident that animal model systems to study the pathogenesis of lung fibrosis in rodents need to be more carefully evaluated with respect to their applicability to human lung fibrosis . As an example, application of the anti-neoplastic drug bleomycin into the lungs of mice or rats is the most frequently employed model system to study the pathogenesis of this devastating lung disease in small laboratory animals. However, in contrast to human lung fibrosis frequently developing over years and decades, bleomycin treatment of mice is characterized by early alveolar epithelial cell injury and acute lung inflammation followed by a transient lung fibrotic response, thereby not adequately reflecting human disease progression . Another previously described model system to initiate fibrosis in small laboratory animals makes use of somatic cell gene transfer of biologically active transforming growth factor-β1, which, due to its unique fibrogenic characteristics, has been shown to trigger slowly progressing lung fibrosis in rats, which is irreversible and lacks acute inflammatory episodes .
Diagnostic tools to assess the degree of lung fibrosis in humans currently include lung CT scans, lung function tests and histological examinations of lung tissue biopsies. In mice, where complete lungs can be removed, particularly histological and biochemical approaches are employed to assess the degree of fibrosis [6, 7]. In recent years, small animal imaging has emerged as an innovative and powerful, non-invasive tool to perform follow up of both acute and chronic disease processes in laboratory animals over time [8–11]. Various non-invasive techniques such as micro-CT, micro-positron emission tomography (micro-PET), and magnetic resonance imaging (MRI) are innovative tools that may be useful for evaluation of novel anti-fibrogenic therapeutic strategies in mouse models in vivo. However, systematic micro-computed tomography-based analysis of developing lung fibrosis induced by gene transfer of biologically active transforming growth factor-β1 as compared to classical histopathological approaches to monitor fibrotic disease progression in mice using both visual assessment and semi-automated quantification algorithms has not been reported so far.
Therefore, in the current study, pulmonary fibrosis was induced in mice by gene transfer of biologically active transforming growth factor-β1, and fibrotic disease progression was evaluated by micro-computed tomography and then correlated to histological scoring of lung fibrosis as pioneered by Ashcroft .
Two adenoviral vectors were employed, which were constructed as outlined in detail previously . AdTGF-β1223-225 contains the cDNA of the coding region of full-length porcine TGF-β1, while containing a mutation of cysteine to serine at positions 223 and 225, rendering expressed TGF-β1 biologically active . Empty viral vector AdDL70-3 constructed as previously described was used as control .
Induction of pulmonary fibrosis
Lung histology and assessment of pulmonary fibrosis
Complete lungs were inflated in situ with PBS-buffered (4.5%, pH 7.0) formalin (Roth, Darmstadt, Germany), and were then carefully removed and immersed in PBS-buffered formalin. Lung tissue samples were paraffin-embedded, and sections of 5 μm were prepared and stained with haematoxylin/eosin (H/E) and Masson's or Elastica-van-Gieson staining. The degree of lung fibrosis was determined according to the method described by Ashcroft and colleagues , employing a numerical scaling system in lung samples ranging from 0 (normal lung) to 8 (total fibrous obliteration of the field). The mean degree of lung fibrosis according to the definitions by Ashcroft et al.  was calculated from individual scores of ~15 microscopic fields analyzed per mouse lung.
Micro-CT scans of mice
A GE Explore Locus cone-beam micro-CT (GE Healthcare, Chalfont St Giles, GB) was used for data-acquisition in prone position under isoflurane inhalation anaesthesia (tube voltage 80 kV, tube current 450 μA, 0.094 mm effective pixel size) with and without respiratory gating (i.e. synchronization of acquisition of micro-CT projections with a timepoint in the respiratory cycle of the individual mouse). Scanning took approximately 10 and 20 min without and with respiratory gating, respectively. Respiratory monitoring was performed using a pressure transducer pad under the animal's chest. Images were reconstructed and assessed at a constant window width/window level (5000/2000).
Assessment of fibrosis by observer rating scale
Three radiologists blinded to the experimental groups and respiratory gating strategies, independently assessed image qualities of all individual entire micro CT scans, using a rating scale ranging from 1 (very good) to 5 (insufficient) and the amount of lung consolidations using a rating scale ranging from 1 (no consolidation) to 5 (very pronounced consolidation).
Due to the diffuse nature of fibrotic consolidations, direct quantification of the pathology of lung fibrosis by measurement of focal fibrotic lesions was not feasible. Therefore, quantification of aerated lung areas known to be compromized by developing lung fibrosis was used as a surrogate marker for progressive lung fibrosis. A region-growing segmentation algorithm was applied on each individual entire micro CT scan using the MevisLab software package (Mevis Research, Bremen, Germany) . The algorithm uses the seed points as starting points. From these points the segmentation volume "grows", with the structures included in the volume depending on the mathematical parameters. The segmentation algorithm was originally developed for quantification of diffuse lung tumors in transgenic mice, and was slightly modified for its use in the current study. Segmentation was carried out using 20-40 seed points, the segmentation threshold tolerance of the region growing algorithm was set to 2%.
All statistical computations have been performed using the statistical program SPSS 17 (SPSS Inc., Chicago, Ill., USA) The significance level of all tests has been set to 0.05 (5%). Because most of the observed variables were ordinal scaled we used non-parametric tests for statistical analysis. Significance of differences in image quality values were tested for all three observers in respiratory gated and respiratory ungated exams using a Wilcoxon signed ranks test. Furthermore, significance of differences in consolidation assessment values for respiratory gated and respiratory ungated exams was tested using a Wilcoxon signed ranks test for all three observers. Differences between the assessment of the 3 observers were evaluated comparing the 95% confidence interval for the mean for assessment values of the image quality for respiratory gated/ungated exams and the assessment values of the consolidation assessment for respiratory gated/ungated exams.
Correlation of the Ashcroft histology scores (assessing the degree of fibrosis) and the lung consolidation values (obtained by assessment of the images by three observers) was calculated for respiratory ungated/respiratory gated exams and all three observers using Spearman's rank correlation coefficient for ordinary scaled variables.
Correlation of the Ashcroft histology scores and the lung segmentation values (obtained by the semiautomated segmentation routine described above) was calculated for respiratory ungated and respiratory gated exams using a Pearson correlation for continuous variables. To assess differences between the 4 time points for histology scores, visual consolidation assessment and segmentation of aerated lung areas 95% confidence intervals were plotted against the 4 time points and differences were tested for significance.
In vivo micro-CT imaging of mouse lungs after intrapulmonary delivery of control vector or AdTGFβ1
As shown in Figure 1 micro-CT imaging of mouse lungs subjected to treatment with control vector or application of AdTGFβ1 for up to 4 weeks revealed that no major consolidations were detectable in control vector treated mice, whereas AdTGFβ1 treatment of mice provoked strong consolidations easily detectable particularly on weeks 2 and 3 post-treatment (Figure 1).
Assessment of image quality and consolidations in respiratory gated as compared to ungated exams of fibrosing lungs
We next questioned whether detection of consolidations in fibrosing mouse lungs would differ between observers. However, initial comparative image quality assessment between three different observers revealed only slightly increased mean image quality values in respiratory gated (MVSD 2.70, MVCVF 2.33, MVTR 2.27) as compared to ungated exams of fibrosing lungs (MVSD 2.33, MVCVF 2.10, MVTR 1.97). Moreover, all three observers determined similar consolidation assessment values for respiratory gated (MVSD 2.13, MVCVF 2.40, MVTR 2.30) as compared to ungated exams (MVSD 2.13, MVCVF 2.43, MVTR 2.40).
Regarding the assessment variability between the 3 observers, we found that the lower and upper bound of the 95% confidence interval for the mean for the 3 observers showed no relevant differences except for the image quality assessment in gated exams. Here, lower and upper bound of the 95% confidence interval for the mean were distinctly higher for observer SD. The Cohen's kappa statistic also showed significant agreement between nearly all observer scores, but no significant agreement was found for the image quality evaluation scores by observers SD and TR for gated (p = 0.051) and ungated exams (p = 0.699).Overall, these data illustrate that visual assessment of micro-CT based lung fibrosis in mice was independent of the observer.
Correlation of histological grading of lung fibrosis with micro-CT based fibrosis grading
Due to the limited number of study objects within each group (AdTGF-β1 for weeks 1-4: n = 3, n = 6, n = 6, n = 3), we were not able to detect significant differences between different time points for visual consolidation assessment and segmentation of aerated lung areas. However, to better exemplify the time course of fibrosis in mice on days 7, 14, 21, and 28 post treatment, mean values and standard deviations for Ashcroft histology scores and visual assessment of all three observers for ungated exams are given (histology for control group: 0,15 ± 0,03; 0,14 ± 0,01; 0,17 ± 0,07; 0,15 ± 0,05/histology for fibrosis group: 1,79 ± 0,58; 2,13 ± 0,5; 1,67 ± 0,84; 1,64 ± 0,45/visual consolidation assessment scores for control group: 1,11 ± 0,19; 1,11 ± 0,19; 1,22 ± 0,38; 1,11 ± 0,19/visual consolidation assessment scores for fibrosis group: 3,56 ± 0,58; 4,06 ± 1,11; 2,50 ± 0,89; 2,00 ± 0,72).
In this study we evaluated the use of micro-CT for assessment of pulmonary fibrosis induced by adenoviral gene transfer of biologically active TGF-β1 in mice. We found micro-CT imaging to be a highly valuable tool to study pulmonary fibrosis progression in mice using both visual inspection of lung images and a semi-automated segementation algorithm that estimates aerated lung areas as a surrogate marker inversely correlating with developing lung fibrosis in mice.
Micro-CT allows a detailed assessment of lung morphology due to high density differences and well defined borders between air and lung tissue. However, one limitation of this method is that the activity of inflammatory processes can only be assessed by secondary characteristics by micro-CT, such as affection of adjacent strucutre or pleural effusion. Other imaging techniques such as magnetic resonance imaging (MRI), micro-positron emission tomography (PET) and optical imaging could provide functional data, with micro-PET beeing an especially promising candidate to robustly assess inflammatory activity by increase of glucose metabolism as detected by 18F-fluoro-deoxyglucose micro-PET [17–19], which might be particularly relevant in terms of monitoring acute exacerbations of pulmonary fibrosis in small laboratory animals.
Image quality is of equal importance in small animal imaging as in human imaging. However, data on imaging protocols are still not extensive. One major issue addressed in previous studies relates to respiratory gating. Various gating techniques, including prospective and retrospective as well as intrinsic and extrinsic techniques have revealed that mostly, gated imaging resulted in increased image quality [20, 21]. In the current study, we found that mean values for image quality were better in gated exams for all observers, however differences were minor and not significant. Both gated and ungated exams resulted in a highly significant correlation with histology-based Ashcroft scores. Therefore, we conclude that gating should be applied whenever possible, but is not an absolutely critical issue when this setup is not available.
A further important aspect of micro CT imaging of normal or fibrosing mouse lungs is the time-point in the respiratory cycle when images are acquired. Due to the typical respiratory pattern in anesthesia, with a short breath followed by a long expiratory plateau, imaging is mostly performed in expiration. Under those experimental conditions, where intratracheal intubation and ventilation is performed, micro-CT imaging can also be performed in inspiration, as a different and consistent respiratory pattern can be created . Without gating, imaging is effectively performed during expiration, due to the longer period of time as compared to inspiration that only represents a small fraction contributing to the average imaging time-point.
The semi-automated segmentation routine described in this study was found to assess consolidations in a quantitative way, as proven by the significant correlation of segmentation values with Ashcroft-based histological grading. To a minor degree, observer interaction is possible, as the observer has to place its own seed points for the region growing algorithm. As seed points have to be placed carefully within aerated lung areas, inadvertent positioning of seed points within non-aerated lung areas may increase the likelihood that consolidated aereas will be included in the segmentation volume and thus may cause assessment bias, which however, according to the presented data, appears to be very low.
In the current study, segmentation volumes were not correlated with total lung volume or body weight of individual mice. However, when considering that differences in expiratory total lung volumes between individual mice of the same age and strain are neglectable, it may be largely excluded that this aspect may have confounded the reported correlations. Future algorithms might take segmentation assessments of the entire lung into account, thereby allowing the correlation of measured consolidations with total lung volumes, or simply by employing a quotient taking the body weight into account.
Additionally, different techniques of fibrosis assessment could be applied. A promising approach is the measurement of CT densities expressed by Hounsfield Units and assessing fibrosis not only by consolidation but also by evaluation of density distribution . Morphological pattern assessment as applied in human fibrosis assessment has also been described for imaging of lung fibrosis in small laboratory animals [24, 25]. However, such approaches would be dependent on an image quality that is consistently comparable to human imaging, which is technically still difficult to achieve in in vivo imaging of small laboratory animals.
In another study we evaluated the applied radiation dose for the protocols used in this study applying phantom and cadaver thermoluminescence dosimetric measurements (own unpublished observations). The expected mean doses in an average C57BL/6 mouse for the respiratory gated and ungated protocols used in this study were 201 mGy and 194 mGy respectively. The expected mean value for an estimated time of 1 minute fluoroscopy to determine the scan field of view was 22 mGy. Day et al. reported that high radiation doses of 7 - 9 Gy did not result in significant lung fibrosis in mice after 90 days, just few microscopically small collagen-rich foci were detected subpleurally . Therefore, it appears that radiation-induced fibrosis can be largely excluded even when serial examinations are performed.
Human lung fibrosis differs from the animal model reported here, as it is irreversible and shows typical morphological patterns in computed tomography representing definite fibrotic changes of the lung. The application of the methods reported for imaging in human lung fibrosis therefore is not intended. We believe that the currently applied technique will primarily be helpful to better assess e.g. novel therapeutic strategies for the treatment of pulmonary fibrosis in rodent model systems.
In conclusion, the current study shows that respiratory gated/ungated micro-CT allows valid in vivo evaluation of the degree of pulmonary fibrosis in mice using both visual assessment and semi-automated quantification algorithms. Micro-CT thus offers a non-invasive imaging tool for evaluation of novel therapeutic strategies for the treatment of pulmonary fibrosis in mice.
This study has been supported by the German Research Foundation, grant SFB 587 to UAM and TW
- Verma S, Slutsky AS: Idiopathic pulmonary fibrosis--new insights. N Engl J Med 2007, 356:1370–1372.View ArticlePubMedGoogle Scholar
- King TE, Schwarz MI, Brown K, Tooze JA, Colby TV, Waldron JA, Flint A, Thurlbeck W, Cherniack RM: Idiopathic pulmonary fibrosis: relationship between histopathologic features and mortality. Am J Respir Crit Care Med 2001, 164:1025–1032.View ArticlePubMedGoogle Scholar
- Coultas DB, Zumwalt RE, Black WC, Sobonya RE: The epidemiology of interstitial lung diseases. Am J Respir Crit Care Med 1994, 150:967–972.View ArticlePubMedGoogle Scholar
- Gauldie J, Kolb M: Animal models of pulmonary fibrosis: how far from effective reality? Am J Physiol Lung Cell Mol Physiol 2008, 294:L151.View ArticlePubMedGoogle Scholar
- Sime PJ, Xing Z, Graham FL, Csaky KG, Gauldie J: Adenovector-mediated gene transfer of active transforming growth factor-beta1 induces prolonged severe fibrosis in rat lung. J Clin Invest 1997, 100:768–776.View ArticlePubMedPubMed CentralGoogle Scholar
- Srivastava M, Steinwede K, Kiviranta R, Morko J, Hoymann HG, Langer F, Buhling F, Welte T, Maus UA: Overexpression of cathepsin K in mice decreases collagen deposition and lung resistance in response to bleomycin-induced pulmonary fibrosis. Respir Res 2008, 9:54.View ArticlePubMedPubMed CentralGoogle Scholar
- van den Brule S, Misson P, Buhling F, Lison D, Huaux F: Overexpression of cathepsin K during silica-induced lung fibrosis and control by TGF-beta. Respir Res 2005, 6:84.View ArticlePubMedPubMed CentralGoogle Scholar
- Cavanaugh D, Johnson E, Price RE, Kurie J, Travis EL, Cody DD: In vivo respiratory-gated micro-CT imaging in small-animal oncology models. Mol Imaging 2004, 3:55–62.View ArticlePubMedGoogle Scholar
- De Clerck NM, Meurrens K, Weiler H, Van DD, Van HG, Terpstra P, Postnov AA: High-resolution X-ray microtomography for the detection of lung tumors in living mice. Neoplasia 2004, 6:374–379.View ArticlePubMedPubMed CentralGoogle Scholar
- Fushiki H, Kanoh-Azuma T, Katoh M, Kawabata K, Jiang J, Tsuchiya N, Satow A, Tamai Y, Hayakawa Y: Quantification of mouse pulmonary cancer models by microcomputed tomography imaging. Cancer Sci 2009, 100:1544–1549.View ArticlePubMedGoogle Scholar
- Ask K, Labiris R, Farkas L, Moeller A, Froese A, Farncombe T, McClelland GB, Inman M, Gauldie J, Kolb MR: Comparison between conventional and "clinical" assessment of experimental lung fibrosis. J Transl Med 2008, 6:16.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashcroft T, Simpson JM, Timbrell V: Simple method of estimating severity of pulmonary fibrosis on a numerical scale. J Clin Pathol 1988, 41:467–470.View ArticlePubMedPubMed CentralGoogle Scholar
- Bett AJ, Haddara W, Prevec L, Graham FL: An efficient and flexible system for construction of adenovirus vectors with insertions or deletions in early regions 1 and 3. Proc Natl Acad Sci USA 1994, 91:8802–8806.View ArticlePubMedPubMed CentralGoogle Scholar
- Maus UA, Janzen S, Wall G, Srivastava M, Blackwell TS, Christman JW, Seeger W, Welte T, Lohmeyer J: Resident alveolar macrophages are replaced by recruited monocytes in response to endotoxin-induced lung inflammation. Am J Respir Cell Mol Biol 2006, 35:227–235.View ArticlePubMedGoogle Scholar
- Taut K, Winter C, Briles DE, Paton JC, Christman JW, Maus R, Baumann R, Welte T, Maus UA: Macrophage Turnover Kinetics in the Lungs of Mice Infected with Streptococcus pneumoniae. Am J Respir Cell Mol Biol 2008, 38:105–113.View ArticlePubMedGoogle Scholar
- Bitter I, Van UR, Wolf I, Ibanez L, Kuhnigk JM: Comparison of four freely available frameworks for image processing and visualization that use ITK. IEEE Trans Vis Comput Graph 2007, 13:483–493.View ArticlePubMedGoogle Scholar
- Grassi R, Cavaliere C, Cozzolino S, Mansi L, Cirillo S, Tedeschi G, Franchi R, Russo P, Cornacchia S, Rotondo A: Small animal imaging facility: new perspectives for the radiologist. Radiol Med 2009, 114:152–167.View ArticlePubMedGoogle Scholar
- Ambrosini V, Nanni C, Pettinato C, Fini M, D'Errico A, Trepidi S, Spinelli A, Al-Nahhas A, Rubello D, Zompatori M, et al.: Assessment of a chemically induced model of lung squamous cell carcinoma in mice by 18F-FDG small-animal PET. Nucl Med Commun 2007, 28:647–652.View ArticlePubMedGoogle Scholar
- Pinkernelle JG, Stelter L, Hamm B, Teichgraber U: [Small animal MRI: clinical MRI as an interface to basic biomedical research]. Rofo 2008, 180:505–513.View ArticlePubMedGoogle Scholar
- Bartling SH, Dinkel J, Stiller W, Grasruck M, Madisch I, Kauczor HU, Semmler W, Gupta R, Kiessling F: Intrinsic respiratory gating in small-animal CT. Eur Radiol 2008, 18:1375–1384.View ArticlePubMedGoogle Scholar
- Bartling SH, Stiller W, Grasruck M, Schmidt B, Peschke P, Semmler W, Kiessling F: Retrospective motion gating in small animal CT of mice and rats. Invest Radiol 2007, 42:704–714.View ArticlePubMedGoogle Scholar
- Shofer S, Badea C, Qi Y, Potts E, Foster WM, Johnson GA: A micro-CT analysis of murine lung recruitment in bleomycin-induced lung injury. J Appl Physiol 2008, 105:669–677.View ArticlePubMedPubMed CentralGoogle Scholar
- Kauczor HU, Heitmann K, Heussel CP, Marwede D, Uthmann T, Thelen M: Automatic detection and quantification of ground-glass opacities on high-resolution CT using multiple neural networks: comparison with a density mask. AJR Am J Roentgenol 2000, 175:1329–1334.View ArticlePubMedGoogle Scholar
- Lee HJ, Goo JM, Kim NR, Kim MA, Chung DH, Son KR, Kim HC, Lee CH, Park CM, Chun EJ, et al.: Semiquantitative measurement of murine bleomycin-induced lung fibrosis in in vivo and postmortem conditions using microcomputed tomography: correlation with pathologic scores--initial results. Invest Radiol 2008, 43:453–460.View ArticlePubMedGoogle Scholar
- Cavanaugh D, Travis EL, Price RE, Gladish G, White RA, Wang M, Cody DD: Quantification of bleomycin-induced murine lung damage in vivo with micro-computed tomography. Acad Radiol 2006, 13:1505–1512.View ArticlePubMedGoogle Scholar
- Day RM, Barshishat-Kupper M, Mog SR, McCart EA, Prasanna PG, Davis TA, Landauer MR: Genistein protects against biomarkers of delayed lung sequelae in mice surviving high-dose total body irradiation. J Radiat Res (Tokyo) 2008, 49:361–372.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.