Abstract
OBJECTIVE: Dual-energy X-ray absorptiometry (DXA)-derived bone mineral density (BMD) often fails to predict fragility fractures. Quantitative textural analysis using magnetic resonance imaging (MRI) may potentially yield useful radiomic features to predict fractures. We aimed to investigate the correlation between BMD and texture attributes (TAs) extracted from MRI scans and the interobserver reproducibility of the analysis.
METHODS: Forty-nine volunteers underwent lumbar spine 1.5-T MRI and DXA. Three-dimensional (3-D) gray-level co-occurrence matrices were measured from routine sagittal T2 fast spin-echo images using the IBEX software. Twenty-two TAs were extracted from 3-D segmented L3 vertebrae. The estimated concordance coefficient was calculated using linear regression analysis. A Pearson correlation coefficient analysis was performed to evaluate the correlation between BMD and the TAs. Interobserver reproducibility was assessed with the concordance coefficient described by Lin.
RESULTS: The results revealed a fair-to-moderate significant correlation between BMD and 13 TAs (r=−0.20 to 0.39; p<0.05). Eight TAs (autocorrelation, energy, homogeneity 1, homogeneity 1.1, maximum probability, sum average, sum variance, and inverse difference normalized) negatively correlated with BMD (r=−0.20 to −0.38; p<0.05), whereas five TAs (dissimilarity, difference entropy, entropy, sum entropy, and information measure corr 1) positively correlated with BMD (r=0.29-0.39; p<0.05). The interobserver agreement was almost perfect for all significant TAs (95% confidence interval, 0.92-1.00; p<0.05).
CONCLUSION: Specific TAs could be reliably extracted from routine MRI and correlated with BMD. Our results encourage future evaluation of the potential usefulness of quantitative texture measurements from MRI scans for predicting fragility fractures.
Magnetic Resonance Imaging; Textural Attribute; Bone Mineral Density
INTRODUCTION
Osteoporosis is the most common metabolic bone disease and is characterized by reduced bone mineral density (BMD), deteriorated bone microarchitecture, and increased bone fragility and susceptibility to fractures (1).
BMD quantified with dual-energy X-ray absorptiometry (DXA) is used as a major determinant of osteoporosis; however, bone mass measurements are insufficient to correctly predict all osteoporotic fragility fractures. In spite of this, DXA continues to be an important tool in the clinical evaluation of osteoporosis. Moreover, DXA is still used as an initial, but insufficient, step in the evaluation of the efficacy of new drugs. This has encouraged the development of other techniques that can reveal bone strength abnormalities. Currently, several bone quality parameters are used with clinical data to improve the prediction of future fragility fractures (2-4).
Together, bone geometry, cortical porosity, collagen properties, bone turnover rates, trabecular microarchitecture, and percentage of microdamage, and bone marrow adiposity (BMAT) are bone strength determinants. In addition, each property may independently contribute to an increased or decreased risk of fracture, even without changes in BMD (5-10). Recent studies have shown that texture-derived features extracted from multidetector computed tomography (MDCT) and DXA scans can indirectly quantify the trabecular bone microarchitecture and thus improve the traditional examination of bone quality and the prediction of fragility fractures (11,12).
The gray-level co-occurrence matrix (GLCM) model is a texture analysis method that uses several statistical parameters to quantify the spatial relationship between diagnostic image pixels, allowing characterization of the morphological properties of underlying tissues (13,14). This model has been validated for characterizing textural attributes (TAs) from MDCT images; however, a few studies have used the same model to extract TAs from magnetic resonance imaging (MRI) scans. Most of these studies extracted texture patterns from the MRI scans of patients with tumors to improve diagnostic accuracy and prognostic prediction (15-17).
MRI is a versatile imaging technique with great potential for enabling the multiparametric analysis of bone quality in a single examination. Advanced MRI sequences such as proton spectroscopy and gradient-echo sequences have been used to quantify BMAT and trabecular bone microarchitecture, respectively (18-25). However, to the best of our knowledge, no study has used routine MRI sequences to extract bone TA, a tool that can be used to indirectly study the microarchitecture of trabecular bone.
This study investigated the interobserver reproducibility of bone texture parameters extracted from routine sagittal T2-weighted fast spin-echo (FSE) lumbar spine (LS) MRI sequences obtained from healthy eutrophic individuals by using GLCM texture features and evaluated their relationship with BMD.
ETHICAL CONSIDERATIONS
Approval by our institutional ethics committee was obtained (47580315.5.0000.5440). All the participants agreed to participate in the study and provided written informed consent. In addition, the procedures were performed in accordance with the Declaration of Helsinki.
MATERIALS AND METHODS
Subjects
The study subjects were comprised 49 volunteers who underwent MRI and DXA between May 2011 and October 2014. Both examinations were performed on the same day. All the individuals were at least 18 years old and had a body mass index (BMI) between 18 and 25 kg/m2. Individuals who used medications that affect bone metabolism, reported histories of smoking and/or alcoholism in the last 10 years, and presented any structural abnormality in the LS that could affect MRI or DXA measurements were excluded from the study.
Dual-energy X-ray absorptiometry
The BMD of the LS (L1-L4) was determined using DXA (Hologic Discovery Wi, QDR series, Hologic, Inc., Waltham, MA, USA). The precision error in the LS was 1.2%.
MRI and TA extraction
All the subjects underwent LS MRI on a 1.5-Tesla system (Philips Achieva, Philips Medical Systems, Best, the Netherlands). Routine sagittal T2-weighted FSE sequences were used to extract TAs (echo time [TE], 120 ms; gap, 4.4 mm; echo-train length, 19; repetition time [TR], 3900 ms; slice thickness, 3.0 mm; scan duration, 2 min 16 s). The L3 vertebrae were manually segmented by two trained raters (JGM, a musculoskeletal radiologist with 3 years of experience in this area) and (LCT, a medical student research fellow). The GLCM characteristics were segmented and extracted using the radiomics analysis platform IBEX v 1.0 (26) as demonstrated in Figure 1.
A 31-year-old female volunteer. Illustrated are the segmentation and positioning of the volume of interest (VOI) in the L3 vertebral body using the IBEX program for volumetric (three-dimensional) extraction of textural attributes (TAs).
The IBEX software characterizes the spatial distribution of gray-level intensities in a region (ROI) or volume of interest (VOI) (27). GLCMs obtained from a VOI or ROI are used to calculate the probability of occurrence of pixel/voxel pairs of gray-level intensities i and j given a distance d from an orientation θ (13,28).
A total of 22 attributes were identified from the L3 segmented VOI. Thus, volumetric three-dimensional (3-D) GLCM (COM3Ds) analyses were performed for the vertebral body. COM3Ds directly computed the occurrence of intensity pairs in the orientation θ for the dimensions x and y, and in the orientation φ for the dimension z. The distances ranged from 1 to 5 voxels, and all 13 orientations of θ and φ (θ=0° and φ=0°, 45°, 90°, and 135°; θ=45° and φ=45°, 90°, and 135°; θ=90° and φ=45°, 90°, and 135°; and θ=135° and φ=45°, 90°, and 135°) were used, resulting in 65 COM3Ds. The 22 GLCM attributes extracted were autocorrelation, cluster prominence, cluster shade, cluster tendency, contrast, correlation, difference entropy, dissimilarity, energy, entropy, two measures of homogeneity (homogeneities 1 and 2), two measures of information correlation (IMC1 and IMC2), inverse difference moment normalized, inverse difference normalized, inverse variance, maximum probability, sum average, sum entropy, sum variance, and variance.
Statistical analyses
The results were analyzed using one-way analysis of variance followed by the Duncan posttest. The confidence interval (CI) was 95%, and the level of significance was set at 0.05. Simple (model 1) and multiple (model 2) linear regression models were tested to verify the association of the parameters. Model 2 was adjusted by age and BMI. In addition, the Pearson correlation coefficient (r) was calculated to measure the strength of the relationship between the variables. The interobserver agreement was evaluated using the concordance test proposed by Lin in 1989 (29). The levels of strength of the agreement associated with kappa statistics were described according to the corresponding kappa ranges as follows: poor, 0.00-0.19; fair, 0.20-0.39; moderate, 0.40-0.59; substantial, 0.60-0.79; and almost perfect, 0.80-1.00 (30). All analyses were performed using the SAS 9.4 software (31).
RESULTS
Patients’ demographics, clinical characteristics, and BMD
The study subjects were comprised 49 healthy individuals with a mean age of 41±14.2 years (range, 20-68 years; median, 39 years) and a mean BMI of 22.6 kg/m2. Of the subjects, 20 were men (mean age, 41±15.5 years; range, 20-68 years) and 29 were women (mean age, 40.9±13.4 years; range, 21-66 years). All nine women aged >50 years were in the postmenopausal stage. The mean BMI was similar between the men and the women (Table 1). The LS BMD results were available for all the subjects. No significant difference in BMD measured with DXA was found between the men and the women. In addition, the differences in the BMDs of the LS (L1-L4) and third lumbar vertebral body (L3) were not significant (p>0.05; Table 1).
Correlation between BMD and bone TAs
The results revealed a fair-to-moderate significant correlation between BMD and 13 TAs derived from MRI scans (r=−0.20 to 0.39; p<0.05). Eight TAs (autocorrelation, energy, homogeneity 1, homogeneity 1.1, maximum probability, sum average, sum variance, and inverse difference normalized) negatively correlated with BMD (r=−0.20 to −0.38; p<0.05), and five TAs (dissimilarity, difference entropy, entropy, sum entropy, and information measure corr 1) positively correlated with BMD (r=0.29-0.39; p<0.05; Table 2).
Estimated coefficients of the association between bone mineral density and bone texture attributes.
Inter-rater reliability
This study evaluated the agreement between two observers (JGM and LCT) for the 13 significant TAs derived from MRI scans. Both researchers performed manual segmentation of the L3 vertebra using the IBEX software and extracted TAs using the software. Both researchers were blinded to the patients’ clinical data and identification. The extracted values were tabulated in Excel, and the agreement coefficient was calculated using the Lin coefficient. The results demonstrated an almost perfect interobserver agreement for all 13 significant TAs extracted (95% CI, 0.92-1.00; p<0.05). The level of agreement of the two observers was greater for the following parameters: autocorrelation, cluster tendency, contrast, difference entropy, dissimilarity, energy, entropy, homogeneity 1, homogeneity 2, inverse different moment normalized, inverse different normalized, inverse variance, maximum probability, sum average, sum entropy, sum variance, and variance.
DISCUSSION
Radiomics is a research area that has shown promising growth and usefulness in diagnostic and prognostic evaluations and therapeutic decision making (32-34). Texture-based quantitative features are used to characterize uniformity, randomness, and repetitive patterns in an image and have been used to detect pathologies and malignancies. Moreover, texture analysis using spectral and structural attributes extracted from MRI scans has been used to differentiate between osteoporotic and malignant compression fractures of vertebral bodies, with promising results reported in the literature (35).
Texture features have been used to study the microarchitectural characteristics of the trabecular bone from high-resolution quantitative computed tomography (HR-pQCT) images. However, no studies have been published using GLCM to extract TAs that may reflect the microarchitecture of the trabecular bone and bone fragility from spinal MRI scans. In our study, we used GLCM to characterize the spatial distribution of the gray-level within volumetric LS images from routine sagittal T2-weighted MRI scans. The choice of sagittal T2-weighted FSE sequences is justified because this sequence is widely available in many routine spinal MRI protocols. The distribution of gray levels is directly impacted by the spatial resolution of MRI sequences. The same attribute is expected to present different values for the same image obtained at different spatial resolutions. In fact, a reported methodology takes advantage of this possible variation and pattern recognition techniques, using attributes extracted from a pyramid of images taken from the same original image at different spatial resolutions. The extent of the effect of this variation on the results depends on the content of the image. We used standardized acquisitions to avoid significant variations in attribute values due to the acquisition parameters.
In this work, the two observers showed good agreement for all TAs extracted from routine MRI, and we found a significant correlation between BMD and 13 TAs. Our results suggest that specific texture parameters could be reliably extracted from routine MRI scans, with clinically acceptable reproducibility for most of the parameters studied.
Our results should encourage future studies to identify isolated attributes or a set of attributes to possibly differentiate subjects at risk of fragility fractures when used with clinical data and other bone quality parameters. Quantitative measurements based on TAs could potentially discriminate subjects with bone mass loss and those at risk of fragility fractures.
This study has some limitations that deserve mention. The cervical and thoracic spinal segments were not studied. This is a relative limitation because currently, only the lumbar vertebrae are used to measure BMD. In addition, our study population was mainly composed of young volunteers with normal bone mass and without fragility fractures. Finally, our study used BMD as a surrogate parameter of fracture susceptibility and did not test the potential usefulness of textural features extracted from MRI scans to directly predict fracture risk. However, our study provides important preliminary data for a more comprehensive study to examine the role of the technique in the evaluation of fracture susceptibility. Future studies are necessary to evaluate the potential clinical application of texture analysis using spine MRI.
CONCLUSION
In summary, several specific TAs could be reliably extracted from routine sagittal T2-weighted MRI scans and showed a strong relationship with BMD in a healthy adult eutrophic population. Our results encourage future evaluations of MRI quantitative textural measurements to potentially discriminate subjects at risk of fragility fractures.
ACKNOWLEDGMENTS
The authors are indebted to Rita C. Varrichio, Luciana Mijoler, and Rodrigo Pessini for their competent technical assistance. Financial Support: JGM: CAPES CODE 1; FJAP: FAPESP 2018/14060-9; MHN-B: FAPESP 2018/04266-9, CNPq 307138/2017-8 and 404795/2018-8; IMA: CNPq 156950/2018-9. FAPESP: São Paulo Research Foundation, CNPq: National Council for Scientific and Technological Development. The funders had no role in the study design.
REFERENCES
-
1 NIH Consensus Development Panel on Osteoporosis Prevention, Diagnosis, and Therapy. Osteoporosis prevention, diagnosis and therapy. JAMA. 2001;285(6):785-95. https://doi.org/10.1001/jama.285.6.785
» https://doi.org/10.1001/jama.285.6.785 -
2 de Paula FJ, de Araújo IM, Carvalho AL, Elias J, Salmon CE, Nogueira-Barbosa MH. The Relationship of Fat Distribution and Insulin Resistance with Lumbar Spine Bone Mass in Women. PLoS One. 2015;10(6):e0129764. https://doi.org/10.1371/journal.pone.0129764
» https://doi.org/10.1371/journal.pone.0129764 - 3 Lang P, Steiger P, Faulkner K, Glüer C, Genant HK. Osteoporosis. Current techniques and recent developments in quantitative bone densitometry. Radiol Clin North Am. 1991;29(1):49-76.
-
4 Bauer JS, Link TM. Advances in osteoporosis imaging. Eur J Radiol. 2009;71(3):440-9. https://doi.org/10.1016/j.ejrad.2008.04.064
» https://doi.org/10.1016/j.ejrad.2008.04.064 -
5 Schuit SC, van der Klift M, Weel AE, de Laet CE, Burger H, Seeman E, et al. Fracture incidence and association with bone mineral density in elderly men and women: the Rotterdam Study. Bone. 2004;34(1):195-202. https://doi.org/10.1016/j.bone.2003.10.001
» https://doi.org/10.1016/j.bone.2003.10.001 -
6 Dall’Ara E, Pahr D, Varga P, Kainberger F, Zysset P. QCT-based finite element models predict human vertebral strength in vitro significantly better than simulated DEXA. Osteoporos Int. 2012;23(2):563-72. https://doi.org/10.1007/s00198-011-1568-3
» https://doi.org/10.1007/s00198-011-1568-3 -
7 Bousson V, Le Bras A, Roqueplan F, Kang Y, Mitton D, Kolta S, et al. Volumetric squantitative computed tomography of the proximal femur: relationships linking geometric and densitometric variables to bone strength. Role for compact bone. Osteoporos Int. 2006;17(6):855-64. https://doi.org/10.1007/s00198-006-0074-5
» https://doi.org/10.1007/s00198-006-0074-5 -
8 Donnelly E. Methods for assessing bone quality: a review. Clin Orthop Relat Res. 2011;469(8):2128-3. https://doi.org/. 10.1007/s11999-010-1702-0
» https://doi.org/10.1007/s11999-010-1702-0 -
9 Maciel JG, de Araújo IM, Carvalho AL, Simão MN, Bastos CM, Troncon LE, et al. Marrow Fat Quality Differences by Sex in Healthy Adults. J Clin Densitom. 2017;20(1):106-13. https://doi.org/10.1016/j.jocd.2016.08.002
» https://doi.org/10.1016/j.jocd.2016.08.002 -
10 Baum T, Grande Garcia E, Burgkart R, Gordijenko O, Liebl H, Jungmann PM, et al. Osteoporosis imaging: effects of bone preservation on MDCT-based trabecular bone microstructure parameters and finite element models. BMC Med Imaging. 2015;15:22. https://doi.org/10.1186/s12880-015-0066-z
» https://doi.org/10.1186/s12880-015-0066-z -
11 Showalter C, Clymer BD, Richmond B, Powell K. Three-dimensional texture analysis of cancellous bone cores evaluated at clinical CT resolutions. Osteoporos Int. 2006;17(2):259-66. https://doi.org/10.1007/s00198-005-1994-1
» https://doi.org/10.1007/s00198-005-1994-1 -
12 Valentinitsch A, Patsch JM, Burghardt AJ, Link TM, Majumdar S, Fischer L, et al. Computational identification and quantification of trabecular microarchitecture classes by 3-D texture analysis-based clustering. Bone. 2013;54(1):133-40. https://doi.org/10.1016/j.bone.2012.12.047
» https://doi.org/10.1016/j.bone.2012.12.047 -
13 Haralick RM, Shanmugam K, Dinstein I. Textural Features for Image Classification. IEEE Trans Syst Man Cybern. 1973;3:610-21. https://doi.org/10.1109/TSMC.1973.4309314
» https://doi.org/10.1109/TSMC.1973.4309314 -
14 Haralick RM. Statistical and Structural Approaches to Texture. Proc IEEE. 1979;67(5):786-804. https://doi.org/10.1109/PROC.1979.11328
» https://doi.org/10.1109/PROC.1979.11328 -
15 Emblem KE, Pinho MC, Zöllner FG, Due-Tonnessen P, Hald JK, Schad LR, et al. A generic support vector machine model for preoperative glioma survival associations. Radiology. 2015;275(1):228-34. https://doi.org/10.1148/radiol.14140770
» https://doi.org/10.1148/radiol.14140770 -
16 Fehr D, Veeraraghavan H, Wibmer A, Gondo T, Matsumoto K, Vargas HA, et al. Automatic classification of prostate cancer Gleason scores from multiparametric magnetic resonance images. Proc Natl Acad Sci U S A. 2015;112(46):E6265-73. https://doi.org/10.1073/pnas.1505935112
» https://doi.org/10.1073/pnas.1505935112 -
17 Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006. https://doi.org/10.1038/ncomms5006
» https://doi.org/10.1038/ncomms5006 -
18 Dooms GC, Fisher MR, Hricak H, Richardson M, Crooks LE, Genant HK. Bone marrow imaging: magnetic resonance studies related to age and sex. Radiology. 1985;155(2):429-32. https://doi.org/10.1148/radiology.155.2.3983394
» https://doi.org/10.1148/radiology.155.2.3983394 -
19 Griffith JF, Yeung DK, Antonio GE, Lee FK, Hong AW, Wong SY, et al. Vertebral bone mineral density, marrow perfusion, and fat content in healthy men and men with osteoporosis: dynamic contrast-enhanced MR imaging and MR spectroscopy. Radiology. 2005;236(3):945-51. https://doi.org/10.1148/radiol.2363041425
» https://doi.org/10.1148/radiol.2363041425 -
20 Patsch JM, Li X, Baum T, Yap SP, Karampinos DC, Schwartz AV, et al. Bone marrow fat composition as a novel imaging biomarker in postmenopausal women with prevalent fragility fractures. J Bone Miner Res. 2013;28(8):1721-8. https://doi.org/10.1002/jbmr.1950
» https://doi.org/10.1002/jbmr.1950 -
21 Yeung DK, Griffith JF, Antonio GE, Lee FK, Woo J, Leung PC. Osteoporosis is associated with increased marrow fat content and decreased marrow fat unsaturation: a proton MR spectroscopy study. J Magn Reson Imaging. 2005;22(2):279-85. https://doi.org/10.1002/jmri.20367
» https://doi.org/10.1002/jmri.20367 -
22 Li X, Kuo D, Schafer AL, Porzig A, Link TM, Black D, et al. Quantification of vertebral bone marrow fat content using 3 Tesla MR spectroscopy: reproducibility, vertebral variation, and applications in osteoporosis. J Magn Reson Imaging. 2011;33(4):974-9. https://doi.org/10.1002/jmri.22489
» https://doi.org/10.1002/jmri.22489 -
23 Chang G, Wang L, Liang G, Babb JS, Wiggins GC, Saha PK, et al. Quantitative assessment of trabecular bone micro-architecture of the wrist via 7 Tesla MRI: preliminary results. MAGMA. 2011;24(4):191-9. https://doi.org/10.1007/s10334-011-0252-0
» https://doi.org/10.1007/s10334-011-0252-0 -
24 Oei L, Koromani F, Rivadeneira F, Zillikens MC, Oei EH. Quantitative imaging methods in osteoporosis. Quant Imaging Med Surg. 2016;6(6):680-98. https://doi.org/10.21037/qims.2016.12.13
» https://doi.org/10.21037/qims.2016.12.13 -
25 Shen Y, Zhang YH, Shen L. Postmenopausal women with osteoporosis and osteoarthritis show different microstructural characteristics of trabecular bone in proximal tibia using high-resolution magnetic resonance imaging at 3 tesla. BMC Musculoskelet Disord. 2013;14:136. https://doi.org/10.1186/1471-2474-14-136
» https://doi.org/10.1186/1471-2474-14-136 -
26 Zhang L, Fried DV, Fave XJ, Hunter LA, Yang J, Court LE. IBEX: an open infrastructure software platform to facilitate collaborative work in radiomics. Med Phys. 2015;42(3):1341-53. https://doi.org/10.1118/1.4908210
» https://doi.org/10.1118/1.4908210 -
27 Echegaray S, Gevaert O, Shah R, Kamaya A, Louie J, Kothary N, et al. Core samples for radiomics features that are insensitive to tumor segmentation: method and pilot study using CT images of hepatocellular carcinoma. J Med Imaging. 2015;2(4):041011. https://doi.org/10.1117/1.JMI.2.4.041011
» https://doi.org/10.1117/1.JMI.2.4.041011 - 28 Haralick R, Shapiro L. Computer and robot vision. Addison-Wesley. 1992.
-
29 Lin LI. A concordance correlation coefficient to evaluate reproducibility. Biometrics. 1989;45(1):255-68. https://doi.org/10.2307/2532051
» https://doi.org/10.2307/2532051 -
30 Landis JR, Koch GG. The measurement of observer agreement for categorical data. Biometrics. 1977;33(1):159-74. https://doi.org/10.2307/2529310
» https://doi.org/10.2307/2529310 - 31 SAS Institute Inc. SAS/STAT user’s guide, version 9.4. Cary, NC: SAS Institute, Inc. 2013.
-
32 Wu W, Parmar C, Grossmann P, Quackenbush J, Lambin P, Bussink J, et al. Exploratory Study to Identify Radiomics Classifiers for Lung Cancer Histology. Front Oncol. 2016;6:71. https://doi.org/10.3389/fonc.2016.00071
» https://doi.org/10.3389/fonc.2016.00071 -
33 Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More than Pictures, They Are Data. Radiology. 2016;278(2):563-77. https://doi.org/10.1148/radiol.2015151169
» https://doi.org/10.1148/radiol.2015151169 -
34 Kumar V, Gu Y, Basu S, Berglund A, Eschrich SA, Schabath MB, et al. Radiomics: the process and the challenges. Magn Reson Imaging. 2012;30(9):1234-48. https://doi.org/10.1016/j.mri.2012.06.010
» https://doi.org/10.1016/j.mri.2012.06.010 -
35 Frighetto-Pereira L, Rangayyan RM, Metzner GA, de Azevedo-Marques PM, Nogueira-Barbosa MH. Shape, texture and statistical features for classification of benign and malignant vertebral compression fractures in magnetic resonance images. Comput Biol Med. 2016;73:147-56. https://doi.org/10.1016/j.compbiomed.2016.04.006
» https://doi.org/10.1016/j.compbiomed.2016.04.006
Publication Dates
-
Publication in this collection
26 Aug 2020 -
Date of issue
2020
History
-
Received
21 Mar 2020 -
Accepted
26 May 2020