Iodine is an essential trace element in the human body that plays a vital role in the neurodevelopment during fetal development and early childhood(Reference Delange1,Reference Gulaboglu, Akgul and Akgul2) . Iodine deficiency during pregnancy and infancy may impair the growth and neurodevelopment of the infant and increase infant mortality(Reference Zimmermann3). The cerebral cortex of the fetus depends on the mother’s thyroxine, which produces triiodothyronine to bind nuclear receptors and exert biological effects(Reference de Escobar, Obregón and del Rey4). A study conducted in Sicily, Italy, found that even mild iodine deficiency during pregnancy can have long-term adverse effects on fetal nerves and cognition, and even if there is subsequent adequate iodine intake during childhood, the effect cannot be improved(Reference Hynes, Otahal and Hay5). Malnutrition during pregnancy is directly associated with goitre and thyroid dysfunction(Reference Alexander, Pearce and Brent6). Similarly, iodine nutrition in breast-feeding women has received attention, because breastmilk is the primary nutrient for infants(Reference Velasco, Bath and Rayman7,Reference Nazeri, Kabir and Dalili8) . Therefore, it is necessary to understand iodine nutrition in pregnant and lactating women and take effective and targeted measures to prevent iodine deficiency.
In 2007, the WHO, UNICEF and International Council for the Control of Iodine Deficiency Diseases issued guidelines for reliably measuring iodine deficiency disorders(9). Four indicators, including the median urinary iodine concentration (MUIC), thyroid swelling rate, thyroid-stimulating hormone (TSH) level and thyroglobulin (Tg) level, have been recommended as indicators of iodine deficiency disorders. However, the MUIC represents only the iodine nutritional indicator of the population and does not reflect an individual’s iodine nutritional status. During pregnancy, the thyroid volume in iodine-sufficient countries increases by 10 %, whereas that in iodine-deficient areas increases by 20 % to 40 %(Reference Alexander, Pearce and Brent6). Studies have found that thyroid volume in pregnant women in iodine-sufficient areas is mainly related to genetics and not to other factors(Reference Soares, Vanacor and Manica10). TSH is not a sensitive biomarker for iodine deficiency or excess iodine and its value in determining mild iodine deficiency or excess iodine is uncertain(Reference Farebrother, Zimmermann and Andersson11). Tg is an indicator of sensitivity to low or high iodine intake(Reference Farebrother, Zimmermann and Andersson11,Reference Zimmermann, Aeberli and Andersson12) . However, the laboratory testing costs for TSH and Tg are high; special laboratories and equipment are required and the subsequent costs are relatively high. Therefore, it is particularly important to select appropriate indicators for evaluating the nutritional status of iodine.
Pregnancy and lactation are special periods because both the growth and development of fetuses and babies require nutrients from their mothers. As an important element in regulating metabolism and synthesising thyroid hormones, iodine plays an important role in the development of fetuses and babies before the age of 2 years. Our research aimed to identify economical, non-invasive and easily accessible biomarkers of individual iodine nutrition to provide clinical guidance for pregnant and lactating women.
Materials and methods
Population and study design
We conducted a cross-sectional study of pregnant and lactating women in northern Xinjiang. This study aimed to identify the best indicators for evaluating iodine nutrition at an individual level. We collected data on the iodine intake, serum total iodine concentration (StI), SIC, UIC, urine creatinine ratio (UIC/Cr), breastmilk iodine concentration (BMIC), thyroid function indicators and thyroid volume. Data were also recorded on socio-economic characteristics, pre-pregnancy BMI, gestational week, number of pregnancies and number of abortions.
Our study population included pregnant and lactating women from the northern part of Xinjiang Uyghur Autonomous Region. The southern part of the region had implemented oral lipiodol pills (twice a year, 100 mg/pill), which would have affected the calculation of iodine intake. Therefore, this region was excluded from analysis. This study included pregnant and lactating women from four counties in the Yili region of northern Xinjiang. Forty pregnant and twenty lactating women participated in the study after providing voluntary informed consent. The researcher informed each participant of their test results.
In this study, the rate of excess iodine intake (iodine intake > 500 µg/l) among pregnant women in northern Xinjiang in 2023 in the ‘Iodine Deficiency Disease Surveillance’ of the whole Xinjiang was 41·52 % as the prevalence of excess iodine intake, and the sample size of the cross-sectional survey was calculated using the formula: When α = 0·05 and the permissible error d = 0·1P, then n 400 × (1 − P)/P. The sample size, n, for the population of pregnant women in the northern border region was calculated to be approximately 564.
Between April and August 2023, 607 pregnant and 171 lactating women were included in the study. The exclusion criteria were as follows: (1) a history of smoking or alcohol abuse; (2) patients who had used iodine-containing lotions in the previous 3 d, taken Huasu amiodarone hydrochloride tablets, or had recently undergone angiography and (3) occupational exposure to iodine (such as medical or personnel use of iodine disinfectants and iodine contrast agents). Our questionnaire consisted of two parts with independent questionnaires: (1) socio-economic characteristics and (2) a 3-day 24-hour dietary recall. The flow chart is shown in Fig. 1. This study was conducted in accordance with the guidelines of the Declaration of Helsinki. All procedures involving human subjects were approved by the Ethics Committee of the Center for Disease Control and Prevention of the Xinjiang Uygur Autonomous Region (number: 2021-07). Written informed consent was obtained from all participants.
Iodine intake survey and calculation
Before starting the survey, our researchers visited the survey area to collect data on foods eaten with high frequency by pregnant and lactating women and the eating habits of the local residents to determine the food’s iodine content. The high-frequency foods included eleven major food categories: staple foods, beans and their products, livestock and poultry meat, eggs, milk and their products, seafood, vegetables, fruits and dried fruits, snacks, condiments and salt. The dietary guidelines of Wang(Reference Zhishuo, Man and Jieshu13) were then used with other dietary maps for 3 d. A retrospective dietary survey was conducted, and the data were collated and analysed. Dietary iodine intake comes mainly from food, drinking water and iodised salt and was calculated using the following formula:
The 20 % was to make allowance for the cooking loss rate of iodised salt defined by WHO. The drinking water iodine concentration data were obtained from the 2017 Xinjiang Water Iodine Concentration Survey Project, and the edible salt iodine concentration data were obtained from the Xinjiang Iodine Deficiency Area Disease Monitoring Project in 2023.
Food, blood, saliva, breastmilk and urine sample collections and laboratory measurement methods
Food sample collection included unprocessed foods such as vegetables, fruits, nuts, meat, dairy products and snacks. Unprocessed flour and edible oil were collected in sealed bags and transported to the laboratory to determine the iodine concentration. Fasting venous blood samples (5 ml) were collected from each participant. The blood sample was allowed to stand upright for 15 min and then centrifuged at 3000 rpm. After centrifugation, the serum was divided into two 1·5 ml centrifuge tubes (Axygen brand). One part was used to determine the total serum iodine concentration, and the other was used to determine thyroid function. Before saliva collection, the participants were asked to rinse their mouths with distilled water. Five minutes after rinsing, the subjects were asked to spit their saliva into a 5 ml sterile plastic tube and seal it. Before milk collection, lactating women were asked to rinse their nipples with distilled water and then squeeze the milk into a 5 ml sterile plastic tube that was then sealed for storage. For the collection of urine, pregnant and breast-feeding women were asked to collect a single dose of urine in the morning before the ultrasound scan and were asked not to eat or drink before collecting the urine. Each participant was given a sterile urine cup and a 5 ml sterile sealed plastic tube. Participants were asked to collect the middle part of the random urine sample into the urine cup and then transfer it to a 5 ml sterile sealed plastic tube. Serum, saliva, breastmilk and urine were transported to the laboratory under cold-chain conditions for measurement. The cold-chain temperature was monitored during this process and was maintained at 2–5°C.
The iodine concentration in food was measured according to the national standard method: National Standard of the People’s Republic of China – National Standard for Food Safety Determination of Iodine in Food Determination of Iodine in Foodstuffs (National Standard No. GB5009·267–2020) Inductively Coupled Plasma Mass Spectrometry (ICP-MS). This method was also used to determine serum, breastmilk and SIC. The UIC was determined in accordance with the Health Industry Standard of the People’s Republic of China – Determination of Iodine in Urine Part 2 (WS/T107·2–2016): Electro-inductive Coupled Isotope Mass Spectrometry (ICP-MS). The urine creatinine (Cr) concentration was determined in accordance with the People’s Republic of China health industry standard ‘urine creatinine spectrophotometry photometric method’. A five-item checklist was used to determine thyroid function. A Roche E411 electrochemiluminescence immunoassay analyser was used to measure the levels of free triiodothyronine (FT3), free thyroxine (FT4), TSH, thyroglobulin antibody and thyroid peroxidase antibody (TPOAb).
Standards used and variable definitions
Pregnancy was classified according to trimester, as follows: first trimester, 0–13 weeks; second trimester, 14–27 weeks and third trimester, 28–40 weeks. The recommended daily nutrient intake for iodine intake by pregnant and lactating women and the daily intake that should not be exceeded were considered. The WHO standard was adopted. Iodine intake of less than 250 µg/d was defined as insufficient iodine intake, whereas iodine intake of more than 500 µg/d was defined as excess iodine intake(Reference Andersson, de Benoist and Delange14). The pre-pregnancy BMI classification standard adopted by the Chinese adult standard in 2007 classified a BMI between 18·5 and 23·9 kg/m2 as normal, less than 18·5 kg/m2 as underweight, between 24 and 27·9 kg/m2 as overweight, and 28 kg/m2 or more as obese(Reference Wang, Mi and Shan15). The WHO, UNICEF and International Council for the Control of Iodine Deficiency Diseases ICCIDD have evaluated the iodine intake categories in pregnant women. An MUIC < 150 µg/l was classified as insufficient iodine intake; An MUIC between 150 and 250 µg/l was classified as adequate iodine intake and a MUIC between 250 and 500 µg/l was classified as more than adequate. An MUIC > 500 µg/l meant that no added health benefit could be expected. The standard for evaluating iodine intake by lactating women was an MUIC < 100 µg/l was classified as insufficient iodine intake. An MUIC > 100 µg/l was classified as adequate iodine intake(9). The evaluation criteria for thyroid hormones and autoantibodies were based on the guidelines for diagnosing and treating thyroid diseases during pregnancy and the postpartum period(16). The reference range for normal thyroid function in pregnant women was established using the Roche electrochemiluminescence immunoassay. The standard established by the commonly used instrument kits was as follows: The reference range of FT3 was 3·1–6·8 pmol/l. The reference ranges of TSH by trimester were 0·09–4·52 mU/l in the first trimester, 0·45–4·32 mU/l in the second trimester and 0·30–4·98 mU/l in the third trimester. The reference range of FT4 was 13·15–20·78 pmol/l in the first trimester, 9·77–18·89 pmol/l in the second trimester and 9·04–15·22 pmol/l in the third trimester. The reference range of thyroglobulin antibody was 0 − 115 IU/ml. The reference range of TPOAb was 0–34 IU/ml. The criteria for diagnosing thyroid disease are listed in Table 1.
TSH, thyroid-stimulating hormone.
Statistical analysis
The Kolmogorov–Smirnov method was used to verify the normality of the data distribution. Age, gestational week, thyroid volume, FT3, FT4, pre-pregnancy BMI, MUIC, UIC/Cr, StI, SIC, TSH, TgAb and TPOAb were all non-normally distributed data and were therefore reported as medians and interquartile ranges. Medcalc 20·010 software was used to draw receiver operating characteristic (ROC) curves to test SIC, BMIC, UIC, StI, UIC/Cr and other indicators for diagnosing excessive iodine intake, and abnormal thyroid function. The Spearman correlation coefficient was used to analyse the relationships between SIC and height, weight, BMI, gestational age, parity, UIC and UIC/Cr. The relationship between SIC and age, BMI, UIC, UIC/Cr, thyroid volume and thyroid function was explored using quantile regression, with α = 0·05 as the threshold for statistical significance. The statistical analysis and plotting of the quantile regression graph were performed using R version 4.03 (The R Foundation for Statistical Computing, Vienna, Austria).
Results
The socio-economic characteristics of the study groups
Of the pregnant participants in this study, 17 % were in their first trimester, and 42·2 % were in their second and 40·9 % were in their third trimester, respectively. The median age was 28 years, and the 25–35-year age group was the largest age group. The mean weight of the women before pregnancy was 57–58 kg, and the mean weight of the lactating women was 63·5 kg. Most participants had a BMI in the normal range (between 18·5 and 23·9 kg/m2). In the lactating women, 39·8 % were overweight and 12·3 % were obese, which was higher than the pre-pregnancy BMI of the pregnant women. The proportion of first birth or first pregnancy was lower than that of previously pregnant women. Among women in the second trimester, 16·2 % were experiencing their first pregnancy or birth, whereas 26·0 % had a previous pregnancy. In the lactating women, the proportion of first delivery or first pregnancy was 46·8 %. The annual income was relatively concentrated in the range of 20 000–30 000 CNY. The level of education was mainly at the junior high school level. The number of pregnant women who had left their jobs was 52·9 %, whereas the proportion of lactating women who had left their jobs was 43·3 %. The proportion of women who had been exposed to second-hand smoke was 18·7 % in lactating women and 9·4 % among women in the second trimester. All pregnant and lactating women consumed iodised salt, and none consumed nutritional iodine supplements (Table 2).
IQR, interquartile range; SIC, salivary iodine concentrations; UIC, urinary iodine concentrations; BMIC, breastmilk iodine concentration.
The SIC in pregnant women was lowest in the first trimester, with a median concentration of 128·10 μg/l. The SIC in pregnant women in the third trimester was higher, with a median concentration of 144·19 μg/l. The median concentration of SIC in lactating women was 131·22 μg/l (Fig. 2). The median UIC of pregnant women in the first trimester was higher than 207·21 μg/l, and the median UIC of lactating women (144·80 μg/l) was lower than that of pregnant women. The median serum total iodine concentration in lactating women was 64·21 μg/l lower than that of the pregnant women. The median UIC and UIC to creatinine (Cr) ratio showed the same trend as that of the median UIC. The UIC/Cr ratio was higher in women in the first trimester and lowest in lactating women. The mean BMIC of the lactating women was 380·26 μg/l. The median dietary iodine intake was over 410 μg/d (Table 2).
Comparison of different indexes in diagnosing excessive iodine intake and abnormal thyroid function
Iodine overdose was investigated using different indicators including UIC, SIC, StI, UIC/Cr and BMIC. If the upper limit of iodine intake was set as 500 μg/d, the cutoff values of UIC, SIC, StI and UIC/Cr for evaluating excessive iodine intake would be 339·67 μg/l, 131·63 μg/l, 73·62 μg/l and 282·04 μg/l, respectively. The sensitivity and specificity of these measures for diagnosis are shown in Table 3. Among them, the UIC of pregnant women had a higher sensitivity in diagnosing iodine excess but lower specificity. The sensitivity and specificity of SIC for diagnosing iodine overdose were similar, with the highest specificity and maximum area under the ROC curve (AUC) of 0·62 (P < 0·01). In lactating women, SIC had the highest specificity for diagnosing iodine overdose (85·00 %), and the highest AUC was 0·63 (P = 0·03). BMIC had the highest sensitivity (82·41 %), but the specificity was low (35·00 %). Similarly, when using different indicators to diagnose abnormal thyroid function in pregnant and lactating women, the highest specificity of the different indicators was SIC, with an AUC of 0·61 and 0·78 in pregnant and lactating women, respectively (P = 0·02, < 0·01). Values greater than 148·12 μg/l and 101·73 μg/l were the best cut-off values for diagnosing abnormal thyroid function in pregnant and lactating women, respectively. UIC had the highest sensitivity for diagnosing abnormal thyroid function in pregnant women, whereas BMIC had the highest sensitivity for diagnosing abnormal thyroid function in lactating women. However, the specificities of UIC and BMIC were low (Table 3 and Fig. 3).
UIC, urinary iodine concentrations; SIC, salivary iodine concentrations; BMIC, breastmilk iodine concentration; AUC, area under the curve.
Quantile regression model of salivary iodine concentrations and related factors
Spearman rank correlation analysis (Table 4) showed that the SIC of the pregnant women was positively correlated with the number of previous deliveries (ρ = 0·0855, P = 0·04) and was positively correlated with UIC and UIC/Cr (ρ = 0·1869, 0·1367, P < 0·01), which were positively correlated with dietary iodine intake and thyroid-related indexes FT3, FT4 and TPOAb. Among the relevant factors for lactating women, only StI was correlated with SIC (ρ = 0·1980, P = 0·009) (Table 5).
SIC, salivary iodine concentrations; UIC, urinary iodine concentrations; TSH, thyroid-stimulating hormone.
UIC, urinary iodine concentrations; TSH, thyroid-stimulating hormone.
By incorporating meaningful variables from the Spearman correlation analysis into the quantile regression model, the results for pregnant women showed that the effect of UIC on SIC increased gradually with an increase in UIC below the 35th percentile; however, the effect on SIC changed little after the 35th percentile. The effect of UIC/Cr on SIC below the 35th percentile gradually decreased as the value increased; however, above the 35th percentile, the effect on SIC gradually increased and exceeded the CI of the mean regression. The effect of 3-day dietary iodine intake on salivary iodine was not significantly different from the mean regression. When the FT3 level was low, it was negatively correlated with SIC, but positively correlated after the 30th percentile. With an increase in FT3, the influence on SIC and its regression coefficient gradually increased. When the concentration of FT4 was low, it was positively correlated with SIC and mildly negatively correlated with SIC with an increase in FT4. Below the 50th percentile, the effect of TPOAb on SIC gradually decreased and then increased from a positive correlation to a negative correlation; however, above the 50th percentile and blow the 65th percentile, the effect gradually increased. Above the 65th percentile, the effect on SIC gradually decreased with increasing TPOAb levels. The effect of StI on SIC below the 65th percentile fluctuated slightly; however, above the 65th percentile, the effect on SIC gradually increased with an increase in StI and changed from negative to positive. Quantile regression was not performed for lactating women because of the limited number of relevant variables (Fig. 4).
Discussion
This study shows that salivary iodine can be used as a biomarker to evaluate excess iodine intake and predict thyroid dysfunction in pregnant and lactating women. Moreover, it is more convenient to obtain materials that are easy to operate and detect. Few studies have been conducted on biomarkers for iodine nutrition in pregnant and lactating women. The Netherlands-based Generation R (iodine sufficient) and the Spain-based INMA (mildly-to-moderately iodine-deficient) cohorts used thyroglobulin as a biomarker of iodine status in iodine-sufficient and mildly iodine-deficient pregnant women(Reference Dineva, Rayman and Levie17). The results of a cross-sectional survey of lactating women in Taiwan and China suggested that measuring maternal UIC alone may not be sufficient, as BMIC, UIC/Cr, and 24-h UIE were all important biomarkers. Ingestion of dairy products and multivitamins was independently associated with BMIC(Reference Huang, Li and Hwu18). In a multicentre Chinese population-based mother and child cohort in Zhejiang, SIC was shown to be a potential biomarker for assessing individual iodine nutrition and thyroid dysfunction in pregnant women(Reference Li, Tu and Gu19). A cross-sectional survey conducted in Tianjin, China, showed that serum iodine levels can be used as an indicator of the nutritional iodine status and thyroid dysfunction(Reference Xu, Guo and Ren20). One study, which started in June 2012 in Liaoning Province in an iodine-adequate area in China, found that UI/Cr reflected the 24-hour iodine excretion and circulating iodine levels during pregnancy and postpartum(Reference Li, Peng and Zhang21). A cohort study in Zhejiang, China, found thyroglobulin to be a sensitive biomarker of iodine status in mildly and moderately iodine-deficient pregnant women(Reference Jiang, Gu and Mo22).
Currently, UIC in a single random urine sample cannot solve the problem of unbalanced iodine intake caused by dietary habits among different individuals. Researchers need to identify an evaluation index suitable for assessing individual iodine nutrition in the population of pregnant and lactating women(Reference Eriksson, Barregard and Sallsten23,Reference Wu, Yang and Su24) . The physiological structure and hormonal profile of this population undergo significant alterations during pregnancy(Reference Alexander, Pearce and Brent6).
Existing research needs to identify a biomarker suitable for individual iodine nutrition in pregnant and lactating women. A meta-analysis found that, although Tg could be used as an evaluation index for iodine deficiency, there are still limitations in adults and pregnant women. In addition, it was not sensitive to changes in iodine intake(Reference Ma and Skeaff25). Ji et al. (Reference Ji, Wu and Wu26) found that StI truly reflects iodine nutritional status in pregnant women. However, venous blood extraction is an invasive process that making it less accepted by pregnant women. SIC is an economical, easy-to-detect, and non-invasive biomarker. The salivary and thyroid glands have similar iodine-enriching functions. Dekker et al. (Reference Dekker, Touw and van der Horst-Schrivers27) used SIC to estimate the iodine status of adults in clinical practice and achieved significant results. Guo et al. (Reference Guo, Pan and Zhang28) found a good correlation between SIC and UIC and concluded that SIC was useful for diagnosing iodine deficiency and excess. SIC can also be used as an indicator of the nutritional iodine status of individuals and populations. Furthermore, it does not have seasonal changes or diurnal variation(Reference Guo, Dong and Jin29). However, if SIC is to be used as a biomarker of iodine excess or iodine deficiency, it needs to be validated in different populations, and appropriate cut-off values for each population should be calculated. Guo et al. (Reference Guo, Pan and Zhang28) also found significant reference values for the diagnosis of iodine deficiency with SIC of less than 105 μg/l and for the diagnosis of iodine excess and thyroid nodules with SIC of more than 273 μg/l. We also found that SIC can be used to diagnose abnormal thyroid function in pregnant and lactating women. SIC values greater than 148·12 μg/l and 101·73 μg/l were the best cut-off values in pregnant and lactating women, respectively. Dietary deficiency or iodine (I) excess plays an important role in the oral mucosa and salivary gland physiology. Absorption occurs through the same sodium iodide symporter (NIS) receptor as that in the thyroid, and this receptor is expressed in the salivary glands, breast, stomach and placenta. At the beginning of gastrointestinal circulation, saliva recovers iodine that is not absorbed by the thyroid or peripheral tissues from the blood. Iodine is transported and accumulated by NIS through the membrane of the ductal salivary basolateral (BL) surface of the salivary cells of the ductal epithelium(Reference De la Vieja and Santisteban30). Gastric salivary clearance and secretion of iodine is an important part of the gastrointestinal iodine cycle, accounting for approximately 23 % of the iodine reservoir in the human body. The salivary and thyroid glands share iodine concentration capacity through sodium iodide homology (NIS) and peroxidase activity(Reference Venturi and Venturi31,Reference Geiszt, Witta and Baffi32) . A correlation between SIC, UIC and UIC ratio was found in pregnant women. Guo et al. (Reference Guo, Pan and Zhang28) found a strong correlation between salivary and urinary iodine(Reference Guo, Pan and Zhang28). Gulaboglul et al. (Reference Gulaboglu, Akgul and Akgul2) also found a close relationship between salivary and urine iodine. SIC determination should be considered as a non-invasive measurement method if it can replace the use of serum for detecting thyroid dysfunction, especially for screening neonatal thyroid dysfunction. The non-invasive saliva collection method is preferable to the collection of neonatal heel blood to screen for hypothyroidism.
This study has some limitations. First, only individuals from iodine-deficient areas were included. Studies should also be conducted among pregnant and lactating women from iodine-adequate and iodine-excess areas. Second, due to the cross-sectional nature of the survey, we did not continuously collect biological samples from pregnant women during the third trimester. As serial samples were not collected during the first, second and third trimester and during lactation, the sample collection process was not dynamic or continuous. Moreover, we did not collect information on the iodine nutritional status of the pregnant women before pregnancy. As with urine, individual variability in saliva was not considered and corrected for. For example, creatinine levels are measured in to correct for urinary iodine levels, and saliva should be measured for salivary amylase to correct for salivary iodine levels(Reference Porcheri and Mitsiadis33). Lactating women cooperated poorly and could not complete the collection of blood, urine, saliva and milk. Most of the participants cared for their children alone at home. This resulted in a low number of breastfeeding women participating in the survey.
Acknowledgements
We thank all the women who participated in this study and the physicians involved in the data collection and testing of the samples. We also thank the Department of Endemic Diseases, the Center for Disease Control and Prevention of Yili Region and the College of Medical Engineering and Technology at Xinjiang Medical University for their help with the analysis data.
This research was funded by the National Natural Science Foundation of China (grant number: 82260648).
Methodology, C. W.; software, R. T.; validation, J. N., H. M. and Y. Z.; formal analysis, C. W.; investigation, C. W. and Y. Z.; resources, C. W.; data curation, Q. L.; writing – original draft preparation, C. W.; writing – review and editing, C. W.; visualisation, D. M.; supervision, D. M.; project administration, C. W.; funding acquisition, C. W. All authors have read and agreed to the published version of the manuscript.
The authors declare no conflicts of interest.
The data presented in this study are available on request from the corresponding author. The data are not publicly available because they contain personal information.