Abstract
BACKGROUND:
Low muscle quantity commonly seen in patients undergoing hemodialysis (HD) is the key contributor of declined physical performance and increases the risk of morbidity and mortality. However, how to normalize muscle mass for operational criteria in this population remains unknown.
OBJECTIVE:
To identify the clinically useful whole body lean tissue mass (LTM) and appendicular skeletal muscle mass (ASM) indices pertinent to physical performance in patients undertaking HD.
METHODS:
Whole body LTM in 38 and ASM in 22 patients undergoing HD were measured by body composition monitor (BCM) and dual-energy X-ray absorptiometry (DXA), respectively. Physical performance was assessed by handgrip strength, the incremental shuttle walk test, sit-to-stand tests, gait speed, Short Physical Performance Battery and Duke Activity Status Index. Besides crude LTM and ASM, the other muscle indices were utilized normalizing for height, height squared, weight, body mass index (BMI), fat mass and body fat%.
RESULTS:
Regardless of BCM or DXA being used, the relationships between different muscle mass indices and physical function were not consistent. While the most useful LTM index which was strongly associated with physical function involved adjustment for height, the strongest (and most useful) ASM index was normalization for BMI.
CONCLUSION:
The superiority of adjustment for BMI or height (height2) recommended by international sarcopenia consensus is also suitable for patients undergoing HD. Patients’ BMI or fat mass should be considered in estimating prevalence of sarcopenia and evaluating relationship between muscle mass and physical performance.
Introduction
Sarcopenia is a progressive and age-related skeletal muscle disorder, including low muscle quantity, strength and physical performance, in which muscle quantity and strength are crucial criteria for diagnosis of sarcopenia [1]. Patients undergoing hemodialysis (HD) are characterized by increased muscle degradation and decreased muscle synthesis due to a variety of conditions inherent to end-stage kidney disease, thereby accelerating the loss of muscle mass and leading to the high prevalence of sarcopenia [2, 3, 4]. In light of different methods and cutoff limit, the prevalence of muscle mass loss varied from 4.0 to 73.5% and sarcopenia from 4.0 to 63.0% in patients undergoing HD [4]. Reports showing that muscle loss lowers physical function, quality of life and contributes to the higher mortality risk in this population [5, 6].
In order to obtain a thorough understanding of sarcopenia, it is necessary to develop accurate measures against muscle mass. Anthropometric screening, Magnetic resonance imaging (MRI), computed tomography (CT), Bioelectrical Impedance Analysis (BIA) and Dual Energy X-Ray Absorptiometry (DXA) were commonly used methods providing assessments of body composition for patients undergoing HD. However, doubts about the accuracy of the results remain. Generally, DXA is reviewed as a reliable and precise instrument when determining muscle quantity. However, the cost of DXA is high and it is not widely available, especially in the community and at poverty areas in the low and middle-income countries. Another disadvantage is that the accuracy of DXA measurements can be influenced by the hydration status of patients undergoing HD [1, 7]. Although BIA is also influenced by hydration status of patients, as a multifrequency bioimpedance spectroscopy device, the body composition monitor (BCM) is regarded a convenient and safe method and frequently used for measuring body composition in patients on HD [8] by determining the quantitative amount of overhydration and lean tissue index [9, 10]. However, BCM is only able to quantitatively assess whole body lean tissue mass (LTM) rather than skeletal muscle mass (SMM) or appendicular skeletal muscle mass (ASM) both of which are recommended for diagnosing sarcopenia in clinical practice and in research [1, 11]. Notably, there are emerging tools that offer innovative ways to assess body composition in patients undergoing HD. A study conducted by Yuri Battaglia et al. employed B-mode ultrasound system to measure the thicknesses of quadriceps rectus femoris muscle (QRFM) and subcutaneous fat tissue (SFT), revealing a significant correlation between muscle ultrasound measurements and bioelectrical impedance vector analysis (BIVA) parameters. Furthermore, it was observed that chronic dialysis patients exhibited decreased thigh SFT thicknesses and ultrasound QRFM, highlighting the potential of ultrasound in assessing muscle mass and fat tissue across diverse clinical settings [12].
An intriguing result has been reported that the prevalence of low muscle mass showed a wide range of variation from 8.1 to 32.4%, even when the same muscle quantities were presented [13]. It follows then that the appropriate normalizing methods used for muscle mass is of great significance. Although the European Working Group on Sarcopenia in Older People (EWGSOP) keeps updating the definition and criteria for sarcopenia [1] and the Asian Working Group on Sarcopenia (AWGS) established criteria for Asian population in 2019 [11], many problems remain to be solved. European or Asian consensus, even Chinese guidelines for sarcopenia were specifically developed for the elderly [14]. There is an ongoing debate about the preferred adjustment for quantifying muscle mass and whether the normalizing method recommended to the older populations can be used for patients undergoing HD. Most studies with patients undergoing HD adopted the criteria recommended by the sarcopenia consensus for the older population [15, 16, 17]. In other words, it is still unknown what is the preferred adjustment for muscle quantity for patients undergoing HD.
Recommendations for measuring muscle quantity in sarcopenia in different guidelines
Recommendations for measuring muscle quantity in sarcopenia in different guidelines
IWGS: The International Working Group on Sarcopenia; FNIH: The Foundation for the National Institutes of Health Sarcopenia; EWGSOP: European Working Group on Sarcopenia in Older People; AWGS: The Asian Working Group on Sarcopenia; DXA: Dual-energy X-ray absorptiometry; BIA: Bioelectrical Impedance Analysis; ASM: Appendicular skeletal muscle mass; ALM: Appendicular lean muscle mass.
The absolute level of SMM or ASM adjusting for body size is the commonly employed metric of muscle quantity. While SMM or ASM divided by height, often height squared are supported as the most frequently employed metric of relative muscle mass [18, 19] in older people by the EWGSOP, AWGS and the international Working Group on Sarcopenia (IWGS) [20], Chinese expert consensus on diagnosis and treatment for elderly with sarcopenia (2021) [14] and the Foundation for the National Institutes of Health Sarcopenia Project suggest adjusting muscle mass for BMI in the diagnosis of sarcopenia [21] (Table 1 shows the recommendations for measuring muscle quantity in sarcopenia from different guidelines). However, seldom studies showed that which normalizing method for muscle mass is suitable for the evaluation of muscle quantity when diagnosing sarcopenia in patients undergoing HD. It is well known that BMI is unable to differentiate between fat and muscle mass, the prevalence of sarcopenia may be overestimated when patients with a greater BMI due to a larger amount of fat mass [13, 22]. Obviously, height squared is positively correlated with BMI, and this index may have the same limitation with the one adjusted for BMI [23]. Due to the negative effect of excessive adipose on physical function [24], muscle mass adjusted for fat mass may be a more useful index for diagnosing low muscle quantity. Wilkinson’s study involving 72 chronic kidney disease patients including 11 patients undergoing HD identified that the most useful muscle mass index pertinent to physical performance may be ASM normalized for adiposity [22]. As the key contributor of the declined physical function [25, 26], the muscle mass indices showing the strong association with physical performance may be the useful operational criteria for predicting clinical states. Obviously, it is important to note that the differences in muscle quantity index may cause inconsistent and erroneous associations between muscle mass and physical function [22], leading to misunderstanding of clinical states depending on evaluation of low muscle mass.
The study was aiming to identify the most useful ASM indices pertinent to physical performance for patients undergoing HD. In light of the widespread use of BCM with this population, the purpose of the study was also to obtain the preferred whole body LTM indices related to physical performance in this population.
This is an exploratory secondary-analysis of body composition data from a study conducted to explore the relationship between physical condition and associated factors in patients receiving HD.
Participants
Patients receiving maintenance HD (receiving treatment for
Eligible patients at a HD unit from the affiliated hospital of xx University in China were identified by unit staff. The patients who expressed an interest in participating in the study were provided with the patient information sheets and discussed the study in more details with the researcher. Informed consent was obtained from each participant at the screening prior to any study-related activities being performed. Participants were recruited between September 2016 and March 2018, and September to December in 2021. All aspects of the study were approved by the Ethical committee of the Affiliated Hospital of Nantong University (Ref. 2015-12) and conducted in accordance with the Declaration of Helsinki
Data collection
Participants were asked to complete a series of outcome measure assessments once only. Based on the underlying fatigue or discomfort after the HD session and and the potential impact of water overload on physical performance, four physical performance tests including the Shuttle Walking Test (SWT), Sit to Stand (STS) 60, Handgrip Strength (HGS) and Short Physical Performance Battery (SPPB), were conducted for evaluating physical function with patients undergoing HD before a HD session, which was with one non-HD day interval between two consecutive HD sessions. Participants were given a chance to become familiar with the tests before they were conducted. During the period, patients completed the questionnaire of the Duke Activity Status Index (DASI). Patients’ body composition was assessed by the Fresenius Body Composition Monitor (BCM) before the dialysis session. Patients were requested to lie flat on a bed or a couch and electrodes of the BCM were placed on one hand and foot. The procedure usually took around two minutes. Additionally, patients were scanned on a GE Healthcare Lunar iDXA scanner to obtain body composition at any time based on their own convenience or preference. They had the option to separate all the assessments into several visits and the total assessments took about 3 hours altogether.
Outcome measures
Demographic and clinical parameters
Demographics and clinical parameters were obtained from their medical records with the patient’s consent, including age, gender, vintage (the duration between the initiation of dialysis and the implementation of the study), comorbidities, medication, bicarbonate (mmol/L) and hemoglobin (g/L). The data was collected on the same day as the investigation of the questionnaire in the study. The clinical parameters were recorded with the pre-dialysis results, which was recorded with the average of three available dialysis days that are closest to the survey day. When there was no data available for the required date, the closest date available was accepted. The results available with
Questionnaire
As a surrogate of peak oxygen uptake (VO2peak) for patients undergoing HD, the Duke Activity Status Index (DASI) has been documented as being a powerful marker of physical function [27]. It consists of two responses ‘yes’ and ‘no’ which refer to 12 activities of personal care, ambulation, household tasks, sexual function and recreational activities. A response of ‘no’ indicates a score of zero, and a response of ‘yes’ to each item is marked as a value based on metabolic cost by calculating Metabolic Equivalent of Task (METs), which means the Energy Expenditure (EE) of a specific physical activity. The question is weighted depending on the level of activity and the values of each specific activity are summed up to produce a numerical total score (0.00-58.2 METs).
Physical performance tests
The Shuttle Walking Test (SWT) is a validated test of maximal exercise capacity in patients receiving HD and has been used as a vital assessment tool for physical function [28, 29]. The SWT includes the Incremental Shuttle Walking Test (ISWT) which was conducted in the study. Participants were asked to walk a 10-meter shuttle course at a speed controlled by an external audible bleep signal. For the ISWT, participants walked back and forth continuously at a progressively increasing pace. The test was terminated when participants failed to complete the shuttle course in the allowed time. The distance in ISWT’s performance was recorded.
The Sit to Stand (STS) 60 has been confirmed as a reliable and valid measure of functional outcome in patients receiving HD [30, 31]. The participants started from a seated position on a hard, upright chair, with arms folded across the chest, the feet flat on the floor and the knees bent at 90∘. For the test, the subject simply stood up fully and then sat down again to the starting position, without using their arms. For the STS60, this is the number of STS cycles achieved in 60 seconds.
The Short Physical Performance Battery (SPPB) has been validated [32] and widely used across diverse populations of elderly people by providing objective function measurements [33, 34, 35]. The SPPB includes a gait speed test, muscle strength test (STS5) and a balance test. Participants were asked to walk a 5 m course twice at their usual walking speed with a walking aid if it was normally used. For increased accuracy in measuring patients’ typical walking speed, the duration time was recorded when they completed a 4-meter section of the course. The average time was obtained from the participant’s two performances. The time taken for the participant to complete the STS cycles over 5 times as quickly as possible was then recorded for the measurement of the STS5. In terms of the balance test, participant attempted to stand with 3 different movements for 10 seconds in each one, including standing still with the feet together side by side (side-by-side), the side of the heel of one foot touching the big toe of the other foot (semi-tandem position), and with the heel of one foot in front of and touching the big toe of the other foot (full tandem position). These three separate tests had their own scoring system, ranging from 0 to 4. The scores from the separate tests were used to calculate the overall SPPB score.
The Handgrip strength (HGS) has been reported as being a good simple measure of muscle strength in arms and often used in the diagnosis of sarcopenia [1]. A hand-held dynamometer was used to measure HGS three times and an average reading was obtained by averaging the three readings. Participants were asked to remain seated with the shoulder adducted elbow flexed at 90∘ and forearm in a neutral position [36, 37, 38]. The dynamometer was held with the side opposite to the vascular access arm by the patients. Participants were asked to squeeze the handle as hard as they could. 30 seconds was allowed for the interval between measurements.
Body composition
As a tool for distinguishing over-hydration, lean and adipose tissue mass, the BCM has been extensively used for understanding the body composition of dialysis patients [39, 40], and has played an important role in evaluating nutritional status [41] and predicting protein-energy wasting [42, 43]. However, the BCM only enables to obtain quantitative assessment of whole-body LTM, instead of ASM or SMM. DXA is viewed the gold standard for estimating three principal components including fat mass, fat free mass and bone mineral content by using X-rays passing through bone and soft tissue. However, its measurements can be influenced by overhydration of the patients.
The following are the different LTM or ASM normalization indices.
LTM or ASM/height (m)
Statistical analysis
Continuous data were checked for normality of distribution using Shapiro-Wilk tests. Whilst variables with normally distributed were presented as mean
All statistical analyses were performed using the statistical package of Social Sciences (IBM SPSS Statistics v.26, New York, USA.) and GraphPad Prism (v.7 Graphpad Software Inc, CA, USA). Statistical significance was taken as
Results
43 patients undergoing HD were approached. 2 refused to take part and 1 withdrew due to a serious medical condition which developed during the study. 38 patients completed body composition assessment conducted with the BCM and physical function tests. 22 patients completed the DXA assessment due to the far distance between the DXA assessment centre and their HD unit.
Demographics and characteristics of the patients undergoing HD using BCM and DXA
Demographics and characteristics of the patients undergoing HD using BCM and DXA
HGS: Handgrip strength; ISWT: Incremental shuttle walk test; STS: Sit-to-stand; SPPB:Short physical performance battery; DASI: Duke activity status index; BMI: Body mass index; ASM: Appendicular skeletal muscle; LTM: Lean Tissue Mass.
The median ages of patients who completed the BCM and DXA were 48.6 and 47.0 years, ranging from 30 to 65 and 35 to 65 years old, respectively. The most frequently reported comorbidities of the two groups were hypertension (BCM: 65.8%; DXA: 54.5%) and ischaemic heart disease (BCM: 23.9% and DXA: 18.1%). The main medications patients administered were anti-hypertensives (BCM: 86.8%; DXA: 54.5%) and oral bicarbonate (BCM: 23.7%; DXA: 31.8%). The means of bicarbonate in the BCM and DXA groups were 21.1 mmol/L and 21.8 mmol/L, both of which were lower than the normal range. Table 2 shows the demographics, clinical parameters, physical function and body composition in patients undergoing HD assessed by the BCM and DXA.
|
|
LTM (kg) | LTM/weight (kg) | LTM/height (m) | LTM/height2 (m2) | LTM/BMI (kg/m2) | LTM/fat mass (kg) | LTM/body fat% (kg) |
|---|---|---|---|---|---|---|---|
| Upper extremities function | |||||||
| HGS (Kg) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| Lower extremities function | |||||||
| ISWT (m) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| STS60 (reps) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| SPPB |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| STS5 (secs) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| GS (m/s) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| Self-reported function | |||||||
|
|
|
|
|
|
|
|
|
| DASI |
|
|
|
|
|
|
|
Averaged correlation coefficient value for each LTM index and their association with physical performance tests
All indices of LTM were positively and significantly correlated with HGS of the non-fistula hand. LTM/height, crude LTM and LTM/body fat% were the top three indices strongly related with HGS (
The strongest LTM normalization index with physical performance tests
LTM/height was observed having the overall strongest association with physical function (
|
|
ASM (kg) | ASM/weight (kg) | ASM/height (m) | ASM/height2 (m2) | ASM/BMI (kg/m2) | ASM/fat mass (kg) | ASM/body fat % (kg) |
|---|---|---|---|---|---|---|---|
| Upper extremities function | |||||||
| HGS |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| Lower extremities function | |||||||
| ISWT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| STS60 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| SPPB |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| STS5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| GS |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| Self-reported function | |||||||
|
|
|
|
|
|
|
|
|
| DASI |
|
|
|
|
|
|
|
Averaged correlation coefficient value for each ASM index and their association with physical performance tests
HGS was significantly and strongly associated with all indices of ASM in patients who completed DXA assessment. In particular, ASM/BMI, ASM/height and crude ASM showed the top three strongest correlation with HGS (
The strongest ASM normalization index with physical performance tests
ASM/BMI was identified as the strongest index associated with overall physical performance tests, followed by crude ASM (
Discussion
Muscle mass loss is a crucial hallmark of sarcopenia in patients receiving HD, contributing to declines in physical performance, poor quality of life and high mortality risk [5, 6, 45]. The loss of muscle mass results in poor outcomes mainly through physical impairment. Therefore, it is clinically useful to identify the muscle mass index that is associated with physical function for patients on HD. In light of the different muscle mass measured by DXA and BCM, the aim of the study was to identify the most useful or the strongest normalization methods in regard to physical function performance for whole body LTM (taken from BCM) and ASM (taken from DXA). The main finding is that different relationship was shown between different LTM or ASM normalization indices and physical function. While LTM adjusted for height (either height or height squared) is the strongest LTM normalization index associated with physical function, the strongest ASM adjustment index is ASM adjusted for BMI.
Fundamentally, individual’s muscle mass is influenced by body size. The useful and well-established normalization of muscle mass is crucial for patients undergoing HD to predict physiological status or negative outcomes. It is well-known that the EWGSOP [1], the AWGS [11] and IWGS [20]support skeletal muscle mass or ASM adjusted for height, particularly height-squared as the most frequently used indexing metric of muscle mass. ASM normalized for height squared have also been commonly used in patients undergoing HD [15, 16, 17]. An observational, longitudinal and multi-center study including 170 older adults on HD revealed that the group with sarcopenia showed a hazard ratio 2.65 (95% CI: 0.86–7.05) for mortality when compared to a group no-sarcopenia and no-malnutrition [16]. Sarcopenia in the patients was diagnosed based on low muscle quantity which was defined by ASM adjusted for height squared and adopted the cutoffs recommended in the 2019 EWGSOP [1]. A Korean study using the cut-off points suggested by the AWGSOP showed that ASM/height2 may be the most valuable indices for predicting physical performance or strength compared with the other ones in male patients on HD [46].
Superiority of ASM/height2 is consistent with our finding in the present study. The difference lies in whole-body LTM which was obtained instead of ASM in the study. The result is in agreement with a previous study by Kittiskulnam et al., which revealed that whole-body LTM taken from pre-dialysis bioimpedance spectroscopy adjusted for height2 was a stronger predictor of mortality compared to alternative indexing metrics, such as body weight, body surface area, and body mass index [13]. The whole body LTM is not recommended as a criterion for diagnosing sarcopenia by international sarcopenia working groups. However, numerous studies have shown that BCM is a robust tool for measuring total body fat and lean body mass in patients undergoing HD [47, 48, 49]. It is also extensively applied in patients as it is a reproducible, more feasible and lower cost method comparing with DXA [50, 51, 52]. The most prominent reason behind its wide application in patients undergoing HD is that its measurement is not affected by overhydration. New research further revealed that body composition parameters obtained using BCM have been confirmed as useful predictors of mortality in this population [48, 53, 54].
There is no doubt that height or height squared as indexing metric of muscle mass is recommended by major consensus panels. However, muscle mass relative to height may be unable to clearly identify patients’ sarcopenia if they are obese or overweight with low muscle quantity [13, 23]. Increasing studies have shown that adjusting for height would not recognize patients with higher fat body mass index and lower muscle mass, a condition that is called “obese sarcopenia” [55]. Wilkinson et al. noted that the strongest and most clinically useful muscle mass index in regard to physical function was ASM adjusted for fat mass [22]. The demographics of all 72 patients with renal disease in the study revealed that their mean BMI was up to 27.6 kg/m2. In terms of patients undergoing HD (
The Foundation for the National Institutes of Health (FNIH) Sarcopenia Project recently derived from nine sources of community-dwelling older persons and proposed the definition of sarcopenia in which ASM was adjusted for BMI [21]. Additionally, Chinese expert consensus on diagnosis and treatment for elderly with sarcopenia (2021) also clearly recommends to identify the absolute level of muscle mass by using ASM indexed to BMI [14]. In patients undergoing HD, our study has consistently shown that ASM normalized for BMI had the strongest (and most useful muscle mass index) association with physical performance. To our knowledge, limited previous studies involving patients undergoing HD diagnosed low muscle mass using the ASM/BMI metric, as the overwhelming researchers adopted ASM/height2 as the optimal diagnostic criteria of sarcopenia. In a study of data from 100 patients with chronic kidney disease not yet on dialysis, the inverse correlation was shown between ASM adjusted for BMI and the serum levels of high-sensitivity C-reactive protein (hsCRP) (
Strengths and limitations
As far as we know, few studies focused on the useful muscle mass indices for predicting physical function in patients receiving HD [46]. Moreover, physical performances and muscle strength tests conducted in the current study are all well-established recommended by the various international sarcopenia consensus [1, 11, 14]. Additionally, whilst DXA was used for identifying low muscle mass, BCM which is commonly used for assessing body composition in patients on HD was applied indexing metric of whole-body lean mass as well. A limitation of the study is the restricted sample size. Additionally, participants in this study were not derived from multiple centers and did not represent a diverse sample of Chinese patients undergoing HD. However, this paper only presents an exploratory secondary analysis of body composition data, providing support for the clinical utility of measuring lean body mass or skeletal muscle mass as a predictor of physical functional and performance. Future studies with a larger sample from multi-centers or regions are warranted to confirm the association between muscle mass indices and physical function performance.
Conclusion
Regardless of DAX or BCM being applied for assessing muscle mass in the current study, inconsistent relationship has been appeared between different muscle mass indices and physical function performance. Our study corroborates the evidence to identify the optimal adjustment indices for predicting low muscle strength or physical performance in patients receiving HD. While the most useful LTM index was adjustment for height, ASM index should normalize for BMI. Consequently, superiority of adjustment for height (height2) or BMI recommended by various international sarcopenia working groups have been consistently shown in patients undergoing HD. However, it is required to take into account the influence of individual’s BMI or obesity on prevalence of low muscle mass.
Relevance to clinical practice
Patients undergoing HD suffer from muscle mass loss that leads to high risk of morbidity and mortality. It is clinically useful to identify the preferred muscle mass normalization method by examining its association with physical function. The exploratory study indicates that ASM indexed for BMI may be clinically and physiologically more relevant in patients undergoing HD. Additionally, although BCM only provides whole body LTM which is probably not pertinent to the definition of sarcopenia, it has been widely used to assess body composition in the patients. Obviously, it is equally imperative in clinical practice to figure out the optimal whole body LTM indices. The study suggests that the most clinically useful LTM index involves adjustment for height. The preferred adjustment method recommended in international sarcopenia consensus can be used for patients undergoing HD. However, it is essential for healthcare professionals or researchers to consider patients’ BMI or obesity when understanding their clinical states based on estimation of low muscle mass.
Footnotes
Acknowledgments
We would like to extend our sincere gratitude to our patients who were willing to participate in our study, conduct a variety of physical performance tests and accept different body composition assessments. The study is supported by the Social Livelihood Science and Technology Fund of Nantong (MS22022064) and the National Social Science Fund of China (21BGL302).
Conflict of interest
The authors declare no competing interests.
Author contributions
Each of the authors has made a significant contribution to the work. The individual contributions of each co-author were as follows: Conception; YS. Data acquisition; YS XRZ JXL. Data analysis and interpretation; YS XRZ. Statistical analysis; YS XRZ. Manuscript preparation; YS.
