{'fraction_non_alphanumeric': 0.03160928641560827, 'fraction_numerical': 0.015899661102841037, 'mean_word_length': 5.174032068830661, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 3, 'https://': 1, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Differences in topological properties of functional brain networks between menstrually-related and non-menstrual migraine without aura

Menstrually-related migraine without aura refers to a specific type of migraine that is associated with the female ovarian cycle. Compared with non-menstrual migraine without aura, in menstrually-related migraine without aura, there are additional attacks of migraine outside of the menstrual period. Menstrually-related migraine without aura tends to be less responsive to acute treatment and more prone to relapse than non-menstrual migraine without aura. Currently menstrually-related migraine without aura is treated no differently from any other migraine but, the differences in the central mechanisms underlying menstrually-related migraine without aura and non-menstrual migraine without aura remain poorly understood. Here, using resting-state functional magnetic resonance imaging and graph theory approaches, we aimed to explore the differences in topological properties of functional networks in 51 menstrually-related migraine without aura patients and 47 non-menstrual migraine without aura patients. The major finding of our study was that significant differences in topological properties between the two groups were mainly evident in the nodal centrality of the inferior frontal gyrus and the thalamus. Nodal centrality in inferior frontal gyrus was negatively correlated with Headache Impact Test questionnaire scores in the menstrually-related migraine without aura patients. Partial least squares correlation analysis revealed enhanced correlations of inferior frontal gyrus to pain-related behavior in the non-menstrual migraine without aura group, while within the menstrually-related migraine without aura group these effects were non-significant. These results indicate that the regulatory mechanisms in the central nervous system may differ between the two subtypes of migraine. The results provide novel insights into the pathophysiology of different subtypes of migraine, and could help us to enhance their clinical diagnosis and treatment.


Introduction
Migraine is an idiopathic headache disorder that is characterized as moderate to severe intensity, and is often associated with a combination of symptoms, including nausea, vomiting, tiredness, phonophobia, and photophobia (May 2009). Menstrually-related migraine without aura (MRM), a particular subtype of migraine without aura, includes additional migraines at menstruation, unlike non-menstrual migraine without aura (NMM) (Güven et al. 2017). Currently, the treatment for MRM is similar to that for NMM and mainly uses triptans and nonsteroidal anti-inflammatory drugs (NSAIDs) (Calhoun 2018), but MRM tends to be more resistant to Yutong Zhang and Tao Xu contributed equally to this work. treatment than NMM (Güven et al. 2017). However, the differences in pathophysiological mechanisms between MRM and NMM are not yet fully understood.
Sex hormones have long been linked with MRM as playing a possible mechanistic role in the condition, as they exert profound influences on the central nervous system (CNS) (Calhoun 2018) that may underlie MRM. A decline in estrogen levels in MRM patients may have a facilitating effect on the development of central sensitization, causing changes in the pain threshold of trigeminal neurons (Martin et al. 2007;Welch et al. 2006). These changes may be responsible for the stronger migraine attacks in patients with MRM rather than NMM (Güven et al. 2017). Additionally, the pathophysiology of MRM probably involves the activation of a number of different "pain processing networks" within the CNS, as well as vasodilatation of meningeal arteries secondary to calcitonin gene-related peptide release from the trigeminal nerve (Martin and Behbehani 2006). Ovarian hormones could modulate these pathways to increase the frequency, severity, or duration of MRM (Martin and Behbehani 2006). In light of such differences in MRM and NMM, an investigation of the pathogenic mechanisms in the two different subtypes of migraine may help the development of effective and innovative treatments for clinical application.
Neuroimaging approaches have been used to measure structural and functional brain changes in migraine patients (TJ & DW Schwedt and Dodick 2009), providing feasible, efficient, and noninvasive tools for investigating the pathophysiological mechanisms of migraine (Liu et al. 2015). Several brain regions in patients with migraine have shown structural or functional alterations, including prefrontal cortex (Gao et al. 2016), anterior cingulate cortex , basal ganglia (BG) , and insula (Liu et al. 2015). Further neuroimaging studies have pointed out that brain dysfunctions associated with migraine may not only occur in one or more isolated regions, but are also closely associated with distributed brain networks (Baliki et al. 2008;Liu et al. 2011). Topological alterations can reveal brain network dysfunction or reorganization in patients, which may affect information segregation and integration (Zhang et al. 2017). The investigation of the segregated and integrated processing in the female migraineurs' brain network will significantly advance our understanding the progression of migraine and improve our ability to implement effective treatments.
The incidence and prevalence of migraine are significantly higher in women than men (Stewart et al. 2008), and functional magnetic resonance imaging (fMRI) results show that women have more vulnerable functional networks than men (May 2009). These sex-specific differences suggest that clinical treatment of women should differ to that of men (von Deneen et al. 2019). As a prevalent subtype of migraine in women, the symptoms of MRM are likely to be similar to those of NMM. However, the pathogenesis and clinical characteristics may partially differ between these two subtypes, which might contribute to variations in the underlying topological features of functional brain networks in each.
In the current study, we aimed to explore differences between MRM and NMM in terms of the topological properties of their underlying functional brain networks. To achieve this, resting-state fMRI and graph theory analysis (GTA) were used to investigate differences in the topological properties of brain networks between MRM and NMM patients.

Methods
This study was approved by the Ethics Committee of the Hospital of Chengdu University of Traditional Chinese Medicine (TCM) and was conducted in accordance with the Declaration of Helsinki. All participants gave written, informed consent to participate after the experimental procedures had been fully explained and they were informed that they could stop participating at any time.
Participants 51 MRM patients (age range, 18-50 years; mean, 33.04 ± 6.43 years) and 47 NMM patients (age range, 18-50 years; mean, 35.30 ± 9.43 years) were included in the final data analysis (Table 1). All patients were enrolled from the outpatient clinic of the Departments of Neurology and Gynecology in two clinical centers: (1) the Hospital of Chengdu University of TCM; and (2) Chengdu Integrated TCM &Western Medicine Hospital. Recruitment took place from June 2015 through August 2018. The diagnosis of MRM and NMM was established according to the ICHD III-beta criteria (Headache Classification Committee of the International Head] 2018). Participants who met all the following inclusion criteria were included in the study: (1) female, 18-50 years old, right-handed; (2) fulfilling criteria for migraine without aura; (3) migraine occurred on day 1 ± 2 of menstruation in at least two out of three menstrual cycles, and additionally at other times of the cycle for the MRM group; (4) migraine was experienced outside of the menstrual cycle in the NMM group; (5) history of migraine without aura for 6 months or more; and (6) regular menstrual cycle with a length of 28 ± 7 days. Patients with any of the following conditions were excluded: (1) macroscopic T2-visible brain lesions on magnetic resonance imaging (MRI) scans; (2) suffered from neurological diseases, immunodeficiency, bleeding disorders, or allergies; (3) MRI contraindications such as claustrophobia; (4) alcohol or drug abuse; (5) pregnancy, lactation, or plans to become pregnant within 6 months.
The drugs used for the prophylaxis of migraine were stopped 4 weeks before the experiment; however, participants were allowed to take emergency medicine when their pain was difficult to endure. Detailed information about the participants' drug intake was recorded. All participants had been free from a typical migraine attack for at least 1 week prior to the MRI examination. MRI was performed during the periovulatory phase (days 12-16 of the menstrual cycle). Urine kits for luteinizing hormone were used to verify that the participants were in their periovulatory phase. A positive urine kit test result indicates a surge in the release of luteinizing hormone indicating impending ovulation. After scanning, all participants reported that they did not experience any headaches or migraines and remained awake during the measurement.
This study was approved by the supervision of the Sichuan Regional Ethics Review Committee on TCM (ethical approval number: 2015KL-004) and was registered in the Chinese Clinical Trial Registry (registration number: ChiCTR-IOR-15006648).

Data acquisition
Information collected included age, height, and weight. All patients were required to keep a migraine diary to record visual analogue scale (VAS) scores (scale of 0-10, with 10 being the most intense pain imaginable), migraine attack frequency (number of times), and duration of migraine attack (h). The Headache Impact Test (HIT-6) questionnaire was adopted to assess the impact of migraine on the lives of the patients (Shin et al. 2008), and the Migraine-Specific Quality of Life Questionnaire (MSQ), a 14-item health-related quality of life questionnaire, was used to measure three dimensions of functional status specific to migraine (using preventive, restrictive, and emotional functional subscales) (Bagley et al. 2012).
MRI data were acquired with a GE Discovery MR750 3.0 T system with an eight-channel, phased-array head coil (General Electric, Milwaukee, WI, USA). The functional images were obtained with a single-shot gradient-echo echoplanar imaging (GRE-EPI) sequence with the following parameters: repetition time = 2000 ms; echo time = 25 ms; flip angle = 90°, field of view = 240 × 240 mm, data matrix = 64 × 64, slice thickness = 3 mm. During the whole functional scan, all participants were instructed to keep their eyes closed and stay awake during the entire session.

Data preprocessing
The first 10 time points of functional data were discarded to eliminate non-equilibrium effects of magnetization and allow subjects to get used to the scanning environment. Data preprocessing of the remaining resting-state images was performed using the Data Processing & Analysis for Brain Imaging toolbox version 2.3 (DPABI v. 2.3; http://rfmri.org/dpabi). All data sets were processed using the following steps: 1) slice timing correction; 2) realignment to correct for spatial displacements due to head motion; 3) spatial normalization to the Montreal Neurological Institute echo-planar imaging template image and resampling to 3-mm isotropic voxels; 4) removal of the linear trend of the smoothed images; 5) regression of nuisance covariates, including 24 head movement parameters (Friston et al. 1996), cerebrospinal fluid signals, and white matter signals; 6) temporal band-pass filtering (0.01-0. 08 Hz). Individuals with an estimated maximum displacement in any direction larger than 2 mm or head rotation larger than The data are shown as the mean ± SD. The independent-sample t-test was used to test the difference in each of the demographic characteristics between the two groups. Abbreviations: MRM, menstrually-related migraine without aura; NMM, non-menstrual migraine without aura; VAS, visual analogue scale; HIT-6, Headache Impact Test version 6. MSQ, Migraine-Specific Quality of Life Questionnaire 2°were discarded from the study. Four patients in the NMM group were excluded from further analysis owing to excessive head movement (>2 mm) during fMRI scanning. No subjects exhibited head motions exceeding 2 mm of translation or 2°of rotation in any direction in the MRM group.

Network construction
We constructed the functional network using the Graph Theoretical Network Analysis (GRETNA) toolbox (https:// www.nitrc.org/projects/gretna/) (Wang et al. 2015); this method of network construction and calculation has been used in previous studies of brain networks (Huang et al. 2017).
For each subject, a functional connection matrix was constructed by defining nodes representing brain regions and edges representing interregional resting-state functional connectivity (RSFC) (Liu et al. 2015;Liu et al. 2012). To define the brain nodes of the functional networks, we used the Human Brainnetome Atlas, which divides the whole brain into 246 regions of interest (ROIs) (Fan et al. 2016). To measure the RSFC, the Pearson correlation coefficients between each pair of ROIs were calculated to form an N × N correlation matrix (where N is the number of ROIs, N = 246).

Network metrics
GTA was employed to compute the topological properties of binarized functional connectivity networks at global and regional (nodal) levels. The global metrics examined in this study included the assortativity coefficient, hierarchical organization, global efficiency and small-world properties. The node metrics examined the nodal degree centrality (DC), which emphasizes the impact and significance of a network at the voxel level. It represents the most local and directly quantifiable centrality measure and reflects the properties of the functional brain network "hub" in relation to network information communication (Telesford et al. 2011).
We calculated the area under the curve (AUC) for each metric using the GRETNA toolbox (Wang et al. 2015). The AUC provides a summarized scalar for brain topological properties independent of single threshold selection and can sensitively detect topological alterations in neurological disorders (Huang et al. 2017;Lei et al. 2015). The matrix was thresholded at a value that yielded a binary undirected network of different link densities ranging from 0%-100% in 1% increments.

Statistical analysis
Software (SPSS version 23.0; SPSS, Chicago, IL, USA) was used for demographic analysis. The independent-sample t-test was used to compare all demographic characteristics between the two groups. A significant difference was indicated when the p value was less than 0.05.
To determine whether there were significant differences in the topological properties of the functional brain networks between the MRM and NMM groups, two-sample t tests with a significance level of α = 0.05 were performed on the AUC of each network metric. To correct for multiple comparisons in each of these values, the Bonferroni correction was applied (p = 0.05/246, corresponding to a threshold of p = 0.0002).
To distinguish the potential effect of central mechanisms from that of clinical pathological parameters between the MRM and NMM groups, the demographic and clinical characteristics of the patients (including age, disease duration, attack frequency, VAS score, and duration of migraine attack) were used as covariates. After controlling for the demographic and clinical characteristics, we re-analyzed the AUC of each network metric. The Bonferroni correction was also used to correct for multiple comparisons (p = 0.05/246, corresponding to p = 0.0002).

Partial least squares correlation analysis
To investigate the relationship between nodal degree centrality and pain-related behavioral variables (attack frequency, VAS score, HIT-6 score, and scores on all subscales of the MSQ), partial least squares correlation (PLSC) was used. PLSC is a correlational technique that analyzes associations between two sets of data (Krishnan et al. 2011), which can analyze the association between multiple behavioral measures and multimodal neuroimaging data in a single statistical model. The correlation matrix generated by PLSC is decomposed into latent variables using singular value decomposition to identify unique correlation patterns. PLSC is ideally suited to deriving underlying latent patterns in highly correlated dependent measures and may have higher sensitivity in identifying complex and potentially weak, yet interrelated patterns; thus, it has already been widely adapted and applied in neuroimaging studies (Keresztes et al. 2017).
Pain is a multi-dimensional subjective experience that involves the physical, emotional, and perceptual integration of noxious information (Blackburn-Munro and Blackburn-Munro 2003), and the PLSC has the advantage of being able to identify associations between pain-related behavioral data and brain activity. Overviews of the computational details of the method have been published by McIntosh and Laugh (2004), and Krishnan et al. (2011).
To determine the statistical significance of each latent variable in the PLSC, a permutation test was used (p < 0.001; 5000 permutation iterations). Five-thousand bootstrap samples were also created to assess the reliability of the effect at each node of the functional brain network. Brain and design saliences were recalculated for every bootstrap sample, yielding bootstrap ratios (BSRs; original saliences/bootstrap standard errors). We used a BSR threshold of ±3.3, which is equivalent to a p value of approximately 0.001.

Demographic and clinical characteristics of subjects
The demographic characteristics of all subjects are summarized in Table 1. We found no statistical differences between the MRM and NMM groups in terms of their age, height, weight, disease duration, attack frequency, VAS score, or duration of migraine attack, or in any subscales of the MSQ. (all p > 0.05). Only HIT-6 showed a statistical difference between the two groups (p < 0.05; Table 1).

Differences in network properties between groups
At the level of the global properties, no significant differences were observed in the AUCs of the hierarchical organization, global efficiency, small-world properties, or assortativity coefficient between the MRM and NMM groups (all p > 0.05; Table 2).
Regarding the nodal DC, compared with the NMM group, the MRM group had higher a degree mainly in the insula, inferior frontal gyrus (IFG), and orbital gyrus (OrG), and a lower degree in the thalamus and BG. After controlling for the demographic and clinical characteristics, group differences between MRM and NMM patients were no longer present in the insula, but they persisted within the IFG, thalamus, and BG ( Fig. 1; Table 3). Additionally, we further investigated the association between the above regions and HIT-6 scores using Pearson correlation analysis. We found that the DC in the IFG of MRM patients was negatively correlated with their HIT-6 scores (p = 0.012, r = −0.349; Fig. 2), but an equivalent correlation was not found in the NMM group.
Taken together, after controlling for the demographic and clinical characteristics, these results indicate the group differences in nodal degree between the MRM and NMM patients were no longer present in the insula, while those in the IFG, thalamus, and BG persisted. Nodal centrality in the IFG of MRM patients was negatively correlated with their HIT-6 scores (p < 0.05).

Results of PLSC analysis
In terms of the relationship between nodal degree and painrelated behavior, only one latent variable showed a significant contribution to the covariance according to a 5000-fold permutation test (34.15%, p < 0.001). The behavioral saliences of the significant latent variable (LV) are shown in Fig. 3a. According to our findings, scores from the HIT-6 and from all subscales of the MSQ significantly contributed to the brain-behavior correlation in the NMM group, while equivalent results were not found in the MRM group.
All significant brain regions showing an association with pain-related behavior are shown in Fig. 3b. For the NMM group, positive associations between HIT-6 score and nodal degree were found in the OrG, superior frontal gyrus, fusiform gyrus, parahippocampal gyrus, precuneus, cingulate gyrus. Negative correlations were found in the IFG, middle frontal gyrus, precentral gyrus, inferior temporal gyrus, posterior superior temporal sulcus, superior parietal lobule, inferior parietal lobule, and insula. The above association was opposite for all subscales of the MSQ.

Discussion
In this study, we investigated differences in the topological properties of functional brain networks between MRM and NMM patients, and further explored differences in nodal centrality between the two groups. The major finding of our study is that the topological properties of the functional network show significant differences between the two groups, mainly in terms of the nodal centrality. Combined with regression analysis, our results suggest that the IFG, thalamus, and BG play important roles in the specific pathophysiology of MRM.
The IFG exhibited a higher nodal centrality in MRM patients than NMM patients. Nodal centrality is a powerful measure for assessing the relative importance of a node in a network, and has been applied to evaluating the information integration ability of an individual brain region within cortical networks (Liu et al. 2012). Our results indicate that there may be a difference between the two groups in terms of information transmission of the IFG in the functional brain network at rest. Furthermore, previous studies of pain found that the IFG not only plays a role in the cognitive modulation of pain (Bingel and Tracey 2008), but is also involved in memory retrieval (Denkova et al. 2013) and emotional regulation of pain (Ochsner et al. 2012). In our results, nodal centrality in the IFG was negatively correlated with HIT-6 score only in The two-sample t-test was used to test the difference in global network properties between the two groups. Abbreviations: MRM, menstruallyrelated migraine without aura; NMM, non-menstrual migraine without aura the MRM patients. The HIT-6 was used to assess the impact of headache on the quality of life of the respondents (Martynowicz et al. 2019). The questions within this test cover the content areas found in widely used measures of headache impact, including pain, social-role limitations, cognitive functioning, psychological distress, and vitality (Kosinski et al. 2003). These findings further indicate that the function of the IFG may relate to regulation of pain and emotion processing, and that the IFG would be a crucial component of the central mechanisms of MRM. However, no significant correlation was found for NMM, suggesting this condition may be mediated by other networks.
In the present study, we found the thalamus exhibited lower nodal centrality in MRM patients compared with NMM patients. Circulating sex hormones may modulate pain sensitivity by affecting neural activity (Veldhuijzen et al. 2013); when estrogen levels are declining rapidly or at very low levels, there is the greatest intensity of pain in premenstrual women with migraine (Maurer et al. 2016). The thalamus is a nuclear complex located in the diencephalon, the deep, central brain structure located between the midbrain and cortex, which plays a key role in the transmission of nociceptive information to the cerebral cortex (Yen and Lu 2013). Neuroimaging studies have reported that patients with migraine show white matter lesions (Yu et al. Fig. 1 Significant between-group differences in the nodal degree centrality between the MRM and NMM groups. (a) Significant differences in the nodal degree centrality between the two groups before regression analysis. (b) Significant differences in the nodal degree centrality between the two groups after regression analysis. Regions colored in red have a significantly higher nodal degree for the MRM group, while regions colored in blue indicate a lower degree for the MRM group. Data shown were visualized with the BrainNet viewer (www.nitrc. org/projects/bnv). Abbreviations: BG, basal ganglia; IFG, inferior frontal gyrus; INS, insula; MRM, menstrually-related migraine without aura; NMM, nonmenstrual migraine without aura; OrG, orbital gyrus; Tha, thalamus The data are shown as the mean ± SD. a before regression analysis; b after regression analysis. The two-sample ttest was used to test the difference in degree centrality between the two groups. The Bonferroni correction was used to correct for multiple comparisons. Abbreviations: MRM, menstrually-related migraine without aura; NMM, non-menstrual migraine without aura 2013), gray matter atrophy (Zhao et al. 2014), and abnormal spontaneous functional activity  in thalamus. Some studies have also found that thalamus is activated during the pain stimulation period of the luteal phase, during which the pain threshold is lower for most pain stimulation procedures (Choi et al. 2006). Current evidence suggests that the relationship of the neural mechanisms of MRM to sex hormone levels is likely to be diverse and complex. As a significant betweengroup difference was found in the functional integration of the thalamus, our findings indicate that the MRM and NMM groups may differ in terms of the thalamus-related pain processing circuit. This difference in the thalamus between the two groups may relate to the fact that migraine attacks occur at different phases in the menstrual cycle in MRM compared to NMM. A prior study reported that hormonal factors may have a facilitating effect on the development of central sensitization (Güven et al. 2017), which enhances central sensitization of the trigeminal system after a decline in estrogen levels in MRM patients (Martin et al. 2007). The activation of the trigeminal system has been found to be associated with ascending nociceptive transmission via the trigemino-thalamo-cortical pathway (Sokolov et al. 2014). Therefore, we speculate that estrogen levels may affect functional brain networks in MRM patients through the central sensitization of the thalamus. Previous work suggests that the BG is involved in most aspects of pain processing, including sensory-discriminative, emotional/affective, and cognitive dimensions of pain, and in pain modulation (Borsook et al. 2010). With regard to migraine, brain imaging studies have shown that the BG serves as a sensitive neuroimaging marker, reflecting the disease duration of female migraine patients without aura (Gao et al. 2016) and that BG activity in response to pain is altered by an increasing frequency of migraine (Maniyar and Goadsby 2013). Migraine attacks in MRM not only occur at menstruation, but also additional attacks can occur at other times during the cycle. The altered nodal centrality in BG may further indicate why migraine attacks are more frequent in MRM than in NMM. Moreover, the alterations in nodal centrality may affect pain processing or modulation pathways, and may even affect functional integration or segregation processing within the functional network (Zhang et al. 2017). In our study, the functional network of the MRM group had a lower nodal degree in the BG than the NMM group, suggesting that pain processing or modulation pathways may be disrupted in MRM, which in line with previous research that the BG might be implicated in impaired pain processing and modulatory processes in migraineurs without aura .
With PLSC analysis, it is possible to explore the relationship between multiple predictor variables and multiple response variables, enabling complex relationships within the data to be evaluated (Weaving et al. 2019). In our study, we used PLSC to analyze the relationship between nodal degree centrality and pain-related behavioral variables between the two groups. Our findings indicate the regulatory mechanisms in the CNS may differ between the two subtypes of migraine. The different central mechanisms occurring in NMM may correlate with painrelated behavioral variables, while those in MRM might be affected by various determinants, including sex hormone fluctuations and genetic variance. Thus, further investigations are needed to elucidate the underlying neural mechanisms of MRM.   There are several limitations to this study. First, because hormonal testing was not carried out, we were unable to compare hormone results between patients in the two groups. In view of the important role of hormones in MRM, hormone testing should be introduced in future studies to explore the influence of hormones on the neuroimaging results. Second, currently there are numerous fMRI studies that have investigated differences between migraine patients and healthy controls. To explore the different central mechanisms underlying the two subtypes of migraine, we chose not to include a group of healthy subjects as a control, which may have prevented us from examining the difference between the two subtypes in depth. Third, the MRM group showed a strong trend compared with the NMM group (Table 1), but because of the limited sample size, no statistical difference was found between the two groups. Further longitudinal fMRI studies will be required to establish and confirm the findings of the current study.

Conclusions
Our findings reveal the different topological properties of functional brain networks in MRM and NMM. Group differences in nodal degree between the two subtypes of migraine were found in the IFG, thalamus, and BG. Combined with regression and PLSC analysis, the results indicate that the regulatory mechanisms in the CNS may differ between the two subtypes of migraine. MRM might be affected by the IFG and thalamus, while the pathophysiology of NMM may be associated with pain-related behavioral variables. When considering treatment options in this patient population, it is important to distinguish NMM from MRM, as the latter includes additional migraines at other times of the month, and it can influence the choice of treatment. (20ZDYF1199), and the First-Class Disciplines Development supported by Chengdu University of Traditional Chinese Medicine (CZYJC1901).

Compliance with ethical standards
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Conflict of interest
The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.