- 1Department of Thoracic Surgery, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
- 2Cardiothoracic Surgery Brain injury and brain protection key laboratory of Ganzhou, Jiangxi, China
- 3Department of Respiratory, Shangyou Hospital of Traditional Chinese Medicine, Ganzhou, China
- 4Department of Respiratory, The First Affiliated Hospital of Gannan Medical University, Ganzhou, China
Background: The telomere has long been regarded as a dependable biomarker for cellular senescence. The lung function can reflect the function and status of the lungs. As individuals age beyond adulthood, there is a gradual decline in lung function. However, the existence of a associated between leukocyte telomere length (LTL) and lung function remains uncertain.
Methods: A two-sample Mendelian randomization (MR) analysis was used. The Single-nucleotide polymorphisms (SNPs) of LTL from the genome-wide association (GWAS) study were used as exposure instruments variable, and the lung function indicator including Forced expiratory volume in 1-s (FEV1), FEV1 Best measure, FEV1 predicted and Forced vital capacity (FVC) from the Neale Lab and MRC-IEU were used as outcomes. The associated between the exposures and outcomes was assessed using inverse-variance weighted (IVW), MR-Egger, and weighted median methods. Sensitivity analysis was conducted using Cochran’s Q-test, MR-Egger intercept test, MR-PRESSO, leave-one-out analysis, and Steriger test.
Results: Using the IVW method, a significant association was identified between genetically determined telomere length extension and enhanced lung function in FEV1, with ukb-a-336 (P = 0.127, OR = 1.028,95CI% = 1.003–1.042) and ukb-b-19657 (P = 7.26E-05, OR = 1.051,95CI% = 1.025–1.077),in FEV1 predicted, ukb-a-234 (P = 0.013, OR = 1.029,95CI% = 1.003–1.042), ukb-b-8428 (P = 0.001, OR = 1.032,95CI% = 1.012–1.052), in FEV1 best measure, ukb-a-231 (P = 7.24E-05, OR = 1.050,95CI% = 1.025–1.075), ukb-b-11141 (P = 1.40E-09, OR = 1.067,95CI% = 1.045–1.090).The sensitivity analysis did not reveal heterogeneity or horizontal pleiotropy.Meanwhile, the Steriger test results also indicate that the directionality between exposure and outcome is correct. Therefore, the results indicated robustness.
Conclusion: There is a correlation between longer LTL and better lung function in the European dataset.
Introduction
The telomere, which is a segment located at the terminal end of a chromosome, plays a crucial role in safeguarding genes from harm, thereby ensuring the stability of chromosomes and the integrity of cells (Blackburn, 2001). Telomeres possess significant biological functions, including the stabilization of chromosome functioning, prevention of degradation and fusion at the terminus, protection of chromosome DNA, and regulation of normal cellular processes (Blackburn, 1991). For quite some time, telomeres have been regarded as a dependable biomarker for cellular senescence (López-Otín et al., 2013). In recent years, aside from their correlation with certain age-related diseases (Sanford et al., 2020; Armanios and Blackburn, 2012), telomere shortening has also been linked to an elevated risk of diabetes (Zee et al., 2010) and tumor development (Rode et al., 2016).
The assessment of pulmonary function can serve as an indicator of the respiratory system’s condition in patients. By conducting lung function examinations, healthcare professionals can not only assess the severity of lung diseases and develop targeted treatment strategies but also aid in the early detection of lung and respiratory diseases. A corresponding study has indicated a correlation between prematurely shortened LTL and the presence of idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary disease (COPD) (Dai et al., 2015). The decline of lung function is intricately linked to various lung diseases such as COPD, IPF, and emphysema. Recent studies have indicated a connection between telomere biology and the development of several lung diseases, including IPF and COPD (Alder et al., 2008; Duckworth et al., 2021). However, there is currently a lack of research investigating the potential involvement of LTL in the progression of lung-related diseases through their impact on lung function.
Mendelian randomization (MR) serves as a viable approach to address the common issues of unmeasured confounding and reverse correlation encountered in conventional observational epidemiology (Smith and Ebrahim, 2003). In the present investigation, we employed a dual sample MR design to assess the plausible correlation between LTL and lung function.
Methods
Date source
The summary of LTL GWAS results was obtained in the genome-wide meta-analysis, which included 472174 European populations, including 45.8% of males and 54.2% of females (Codd et al., 2021). In the relevant evaluation indicators of lung function, FEV1, FEV1 Best measure, FEV1 predicted and FVCwere selected as the evaluation indicators. For each indicator was driven by two different research consortiums, Neale Lab and MRC-IEU. The dataset ukb-b-19657 of FEV1 includes 421986 cases, and ukb-a-337 includes 307638 cases In FEV1 Best measure, ukb-b-11141 contains 345665 cases, ukb-a-231 contains 255492 cases, ukb-b-8428 in FEV1 predicted contains 148653 cases, ukb-a-234 contains 110423 cases, and in FVC, ukb-a-336 contains 307638 cases., Ukb-b-7953 contains 421986 cases, all of which are from European men and women.
In our meta-regression analysis, we utilized publicly accessible summary data derived from the open GWAS database pertaining to telomere length and lung function datasets. It is worth noting that moral voting was not deemed relevant or applicable to our study.
Selection of instrumental variables
To be selected as instrumental variables for screening LTL genetic variants, SNPs had to meet three key assumptions: 1) strong association with exposure factors, 2) no association with any confounders of exposure factor-outcome associations (horizontal pleiotropy), and 3) the genetic variant only affects the outcome through exposure factors and not through other pathways.In the process of instrumental variables (IV) selection, the initial criterion employed is whole genome significance (p < 5 × 10−8). Additionally, to eliminate the presence of Linkage disequilibrium, R2 < 0.001 and a window size of 10000 kb are set. Consequently, a total of 154 SNPs are ultimately identified as IVs for telomere length. It is noteworthy that the F-statistic for all these genetic variations surpasses the threshold of 10.
Mendelian randomization analysis
In the course of the analysis, the IVW method served as the primary approach for assessing the correlation between LTL and lung function. Additionally, the MR Egger, Weighted median, Simple mode, and Weighted mode methods were employed for validation purposes. The evaluation of IVs’ heterogeneity was conducted through the utilization of MR Egger regression (MR Egger intercept test) and IVW in Cochran’s Q statistic (Greco et al., 2015). A significance level of p > 0.05 was employed to determine the absence of heterogeneity or horizontal pleiotropy. Simultaneously, we employed the Mendelian randomization Randomized Multivalidity Residual and Outlier (MR-PRESSO) test to detect and eliminate outliers at the level (Verbanck et al., 2018). Throughout this process, we utilized retention analysis to mitigate the potential influence of a single SNP on the overall findings. Furthermore, we employed the Steiger test to compute the variance explained by the instrumental variables estimation on both the exposure and outcome variables, and to determine if the outcome’s variance (c) is indeed lower than that of the exposure. If this condition is met, the direction is deemed accurate (Hemani et al., 2017).
In this study, statistical analyses were performed using R code (version 4.2.2). MR analyses were conducted using the “TwoSampleMR” package (version 0.5.6) (Hemani et al., 2018)and the “MR-PRESSO” package (version 0.5.6).
Results
Based on the results of the MR analysis, a potential correlation was observed between extended LTL and better lung function indicators, FEV1, predicted FEV1, the best measure of FEV1, and FVC. Employing the IVW method, a significant correlation was identified between genetically determined LTL extension and enhanced lung function in FEV1, with ukb-a-336 (P = 0.127, OR = 1.028,95CI% = 1.003–1.042) and ukb-b-19657 (P = 7.26E-05, OR = 1.051,95CI% = 1.025–1.077), in FEV1 predicted, ukb-a-234 (P = 0.013, OR = 1.029,95CI% = 1.003–1.042), ukb-b-8428 (P = 0.001, OR = 1.032,95CI% = 1.012–1.052), in FEV1 best measure, ukb-a-231 (P = 7.24E-05, OR = 1.050,95CI% = 1.025–1.075), ukb-b-11141 (P = 1.40E-09, OR = 1.067,95CI% = 1.045–1.090), in FVC, ukb-a-336 (P = 0.024, OR = 1.022,95CI% = 1.003–1.042), ukb-b-7953 (P = 0.038, OR = 1.018,95CI% = 1.001–1.036) (Table 1).
Table 1. Genetic prediction of telomere length (TL) and its association with lung function indicators (FEV1,FEV1 predicted,FEV1 best measure and FVC.
Furthermore, we employed various methods to assess the results obtained from the IVW method, all of which yielded consistent directions and similar effects (Figure 1). Our analysis of FVC outcomes using Cochran’s Q test and MR Egger test revealed varying degrees of heterogeneity (Table 2). Additionally, we conducted statistical analysis using the intercept term egger of MR egger_ intercept, which yielded a p-value greater than 0.05, indicating the absence of pleiotropy (Table 2). However, we conducted a retest of pleiotropy levels using MR-PRESSO and identified the presence of pleiotropy in FVC-related outcomes (ukb-a-336, P = 0.031, ukb-b-7953, P = 0.0016).
Figure 1. Scatter plot of the association between the telomere length (TL) and lung function indicators (FEV1,FEV1 predicted,FEV1 best measure and FVC).
Table 2. Heterogenity test and Horizontal pleiotropy test of telomere length (TL) and its association with lung function indicators (FEV1,FEV1 predicted,FEV1 best measure and FVC.
In the context of sensitivity analysis, the exclusion of a single SNP does not significantly alter the correlation estimation of genetic predicted LTL on lung function (Figure 2). This implies that the removal of any specific SNP does not fundamentally impact the results, thereby ensuring the stability of our dual sample MR analysis. Additionally, the Steriger test reveals that the variance explained by instrumental variables estimation in the outcome is lower than that in the exposure. This suggests that the instrumental variables estimation influences the outcome primarily through its effect on the exposure, as supported by the data presented in Supplementary Table S1.
Figure 2. Leave-one-out analyses to evaluate whether any single instrumental variable was driving the correlation effect.
Discussion
In this randomized study utilizing a two-sample MR approach, our findings indicate a positive association between longer LTL and improved lung function indicators, encompassing FEV1, predicted FEV1, and the most accurate measure of FEV1. Consequently, acquiring an in-depth comprehension of the interplay between LTL and lung function holds immense importance in the realm of lung disease prevention and treatment.
LTL functions as a crucial indicator reflecting the cellular biological condition and overall health status. The shortening of LTL is associated with an augmented risk of cardiovascular diseases, diabetes, and certain cancers (Pauleck et al., 2023; Haycock et al., 2014; Ahmed and Lingner, 2018). Relevant studies have demonstrated that neutrophils induce telomere dysfunction through reactive oxygen species (ROS) dependence (Lagnado et al., 2021). Simultaneously, LTL is prone to ROS-induced impairment, and oxidation can result in telomere shortening (Ahmed and Lingner, 2018). Nevertheless, the length of LTL can be maintained via the UFMylation of MRE11 (Lee et al., 2021). The regulatory mechanism of LTL encompasses multiple aspects such as telomerase activity, oxidative stress, inflammatory response, and genetic factors, yet the specific one remains ambiguous. LTL is a vital regulator of Cellular senescence and genome stability. The deterioration of lung function indicates the aging process. Research indicates that individuals aged 25 to 40 with decreased lung function (FEV<80% predicted) in early adulthood have a higher risk of developing respiratory and cardiovascular diseases, as well as a greater all-cause mortality rate (Agustí et al., 2017).Numerous studies have verified the correlation between LTL and the lung disease. Some research disclosed that individuals with shorter telomere length presented reduced lung function and an elevated risk of developing COPD (Rode et al., 2013; Savale et al., 2009). At the same time, forced expiratory volume (FEV) decreased significantly more in smokers with shorter LTL than in smokers with longer LTL (Andujar et al., 2018). Alder et al. conducted research on the role of telomere length as a genetic predisposing factor for Emphysema by inducing a chronic cigarette smoke model (Alder et al., 2011). Both LTL and lung function decline with age. At the same time, there was a moderate correlation between shorter telomere length in middle-aged people and compromised airflow parameters, although there was no significant correlation with FVC, which is consistent with our own findings (Nguyen et al., 2019). Studies have shown that telomere shortening has been identified as a characteristic of idiopathic interstitial pneumonia and may lead to idiopathic sexual organ failure, which is clinically manifested by lung and liver problems (Alder et al., 2008; Duckworth et al., 2021). Additionally, it has been observed that individuals with shorter telomere length have an increased risk of death and a poorer prognosis among lung cancer patients (Zhang et al., 2015; Jang et al., 2008; Doherty et al., 2018; Astuti et al., 2017).
It should be emphasized that telomere length can be affected by various genetic and environmental factors, including age, lifestyle choices, psychological factors, and the existence of other disease conditions. Numerous studies have indicated that although there is currently a lack of direct evidence to prove that exercise can prolong LTL, a significant difference in LTL is observed between regular exercisers and non-exercisers, as demonstrated by comparative analyses of exercise and non-exercise groups (Cherkas et al., 2008; Venturelli et al., 2014; Borghini et al., 2015; Silva et al., 2016). Hence, adopting a healthier lifestyle and avoiding smoking may contribute to maintaining longer telomere length and optimal lung function.
In this study, there are indeed certain limitations. Our investigation was solely centered on the relationship between LTL and lung function. While we recognize that other factors such as genetic variations and environmental influences may also play a role in lung function. Additionally, it should be noted that our study only included individuals of European descent. Hence, future data from other populations are required to verify this result.
Conclusion
Our research findings indicate a correlation between increased LTL and enhanced lung function. Integrating LTL measurement into clinical practice holds promise for valuable insights into individual risk assessment, disease prognosis, and potential treatment approaches. Further investigation is warranted to elucidate the underlying mechanisms and assess the potential of LTL-related interventions in ameliorating lung function and health.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
SZ: Data curation, Formal Analysis, Funding acquisition, Software, Writing–original draft, Writing–review and editing. WZ: Methodology, Validation, Writing–review and editing. DR: Validation, Visualization, Writing–review and editing. ZT: Validation, Visualization, Writing–review and editing. XL: Conceptualization, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by Health Commission science and technology plan of Jiangxi Provincial, CHINA (SKJP220229742).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2024.1373064/full#supplementary-material
References
Agustí A., Noell G., Brugada J., Faner R. (2017). Lung function in early adulthood and health in later life: a transgenerational cohort analysis. Lancet Respir. Med. 5 (12), 935–945. doi:10.1016/s2213-2600(17)30434-4
Ahmed W., Lingner J. (2018). PRDX1 and MTH1 cooperate to prevent ROS-mediated inhibition of telomerase. Genes Dev. 32 (9-10), 658–669. doi:10.1101/gad.313460.118
Alder J. K., Chen J. J., Lancaster L., Danoff S., Su S. C., Cogan J. D., et al. (2008). Short telomeres are a risk factor for idiopathic pulmonary fibrosis. Proc. Natl. Acad. Sci. U. S. A. 105 (35), 13051–13056. doi:10.1073/pnas.0804280105
Alder J. K., Guo N., Kembou F., Parry E. M., Anderson C. J., Gorgy A. I., et al. (2011). Telomere length is a determinant of emphysema susceptibility. Am. J. Respir. Crit. Care Med. 184 (8), 904–912. doi:10.1164/rccm.201103-0520OC
Andujar P., Courbon D., Bizard E., Marcos E., Adnot S., Boyer L., et al. (2018). Smoking, telomere length and lung function decline: a longitudinal population-based study. Thorax 73 (3), 283–285. doi:10.1136/thoraxjnl-2017-210294
Armanios M., Blackburn E. H. (2012). The telomere syndromes. Nat. Rev. Genet. 13 (10), 693–704. doi:10.1038/nrg3246
Astuti Y., Wardhana A., Watkins J., Wulaningsih W.PILAR Research Network (2017). Cigarette smoking and telomere length: a systematic review of 84 studies and meta-analysis. Environ. Res. 158, 480–489. doi:10.1016/j.envres.2017.06.038
Blackburn E. H. (1991). Structure and function of telomeres. Nature 350 (6319), 569–573. eng. doi:10.1038/350569a0
Blackburn E. H. (2001). Switching and signaling at the telomere. Cell 106 (6), 661–673. doi:10.1016/s0092-8674(01)00492-5
Borghini A., Giardini G., Tonacci A., Mastorci F., Mercuri A., Mrakic-Sposta S., et al. (2015). Chronic and acute effects of endurance training on telomere length. Mutagenesis 30 (5), 711–716. doi:10.1093/mutage/gev038
Cherkas L. F., Hunkin J. L., Kato B. S., Richards J. B., Gardner J. P., Surdulescu G. L., et al. (2008). The association between physical activity in leisure time and leukocyte telomere length. Arch. Intern Med. 168 (2), 154–158. doi:10.1001/archinternmed.2007.39
Codd V., Wang Q., Allara E., Musicha C., Kaptoge S., Stoma S., et al. (2021). Polygenic basis and biomedical consequences of telomere length variation. Nat. Genet. 53 (10), 1425–1433. doi:10.1038/s41588-021-00944-6
Dai J., Cai H., Li H., Zhuang Y., Min H., Wen Y., et al. (2015). Association between telomere length and survival in patients with idiopathic pulmonary fibrosis. Respirology 20 (6), 947–952. doi:10.1111/resp.12566
Doherty J. A., Grieshober L., Houck J. R., Barnett M. J., Tapsoba J. D., Thornquist M., et al. (2018). Telomere length and lung cancer mortality among Heavy smokers. Cancer Epidemiol. Biomarkers Prev. 27 (7), 829–837. doi:10.1158/1055-9965.Epi-17-1183
Duckworth A., Gibbons M. A., Allen R. J., Almond H., Beaumont R. N., Wood A. R., et al. (2021). Telomere length and risk of idiopathic pulmonary fibrosis and chronic obstructive pulmonary disease: a mendelian randomisation study. Lancet Respir. Med. 9 (3), 285–294. doi:10.1016/s2213-2600(20)30364-7
Greco M. F., Minelli C., Sheehan N. A., Thompson J. R. (2015). Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat. Med. 34 (21), 2926–2940. doi:10.1002/sim.6522
Haycock P. C., Heydon E. E., Kaptoge S., Butterworth A. S., Thompson A., Willeit P. (2014). Leucocyte telomere length and risk of cardiovascular disease: systematic review and meta-analysis. Bmj 349, g4227. doi:10.1136/bmj.g4227
Hemani G., Tilling K., Davey Smith G. (2017). Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 13 (11), e1007081. doi:10.1371/journal.pgen.1007081
Hemani G., Zheng J., Elsworth B., Wade K. H., Haberland V., Baird D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife 7 (7), e34408. doi:10.7554/eLife.34408
Jang J. S., Choi Y. Y., Lee W. K., Choi J. E., Cha S. I., Kim Y. J., et al. (2008). Telomere length and the risk of lung cancer. Cancer Sci. 99 (7), 1385–1389. doi:10.1111/j.1349-7006.2008.00831.x
Lagnado A., Leslie J., Ruchaud-Sparagano M. H., Victorelli S., Hirsova P., Ogrodnik M., et al. (2021). Neutrophils induce paracrine telomere dysfunction and senescence in ROS-dependent manner. Embo J. 40 (9), e106048. doi:10.15252/embj.2020106048
Lee L., Perez Oliva A. B., Martinez-Balsalobre E., Churikov D., Peter J., Rahmouni D., et al. (2021). UFMylation of MRE11 is essential for telomere length maintenance and hematopoietic stem cell survival. Sci. Adv. 7 (39), eabc7371. doi:10.1126/sciadv.abc7371
López-Otín C., Blasco M. A., Partridge L., Serrano M., Kroemer G. (2013). The hallmarks of aging. Cell 153 (6), 1194–1217. doi:10.1016/j.cell.2013.05.039
Nguyen M. T., Saffery R., Burgner D., Lycett K., Vryer R., Grobler A., et al. (2019). Telomere length and lung function in a population-based cohort of children and mid-life adults. Pediatr. Pulmonol. 54 (12), 2044–2052. doi:10.1002/ppul.24489
Pauleck S., Sinnott J. A., Zheng Y. L., Gadalla S. M., Viskochil R., Haaland B., et al. (2023). Association of telomere length with Colorectal cancer risk and prognosis: a systematic review and meta-analysis. Cancers (Basel) 15 (4), 1159. doi:10.3390/cancers15041159
Rode L., Bojesen S. E., Weischer M., Vestbo J., Nordestgaard B. G. (2013). Short telomere length, lung function and chronic obstructive pulmonary disease in 46,396 individuals. Thorax 68 (5), 429–435. doi:10.1136/thoraxjnl-2012-202544
Rode L., Nordestgaard B. G., Bojesen S. E. (2016). Long telomeres and cancer risk among 95 568 individuals from the general population. Int. J. Epidemiol. 45 (5), 1634–1643. doi:10.1093/ije/dyw179
Sanford S. L., Welfer G. A., Freudenthal B. D., Opresko P. L. (2020). Mechanisms of telomerase inhibition by oxidized and therapeutic dNTPs. Nat. Commun. 11 (1), 5288. doi:10.1038/s41467-020-19115-y
Savale L., Chaouat A., Bastuji-Garin S., Marcos E., Boyer L., Maitre B., et al. (2009). Shortened telomeres in circulating leukocytes of patients with chronic obstructive pulmonary disease. Am. J. Respir. Crit. Care Med. 179 (7), 566–571. doi:10.1164/rccm.200809-1398OC
Silva L. C., de Araújo A. L., Fernandes J. R., Matias Mde S., Silva P. R., Duarte A. J., et al. (2016). Moderate and intense exercise lifestyles attenuate the effects of aging on telomere length and the survival and composition of T cell subpopulations. Age (Dordr) 38 (1), 24. doi:10.1007/s11357-016-9879-0
Smith G. D., Ebrahim S. (2003). Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int. J. Epidemiol. 32 (1), 1–22. doi:10.1093/ije/dyg070
Venturelli M., Morgan G. R., Donato A. J., Reese V., Bottura R., Tarperi C., et al. (2014). Cellular aging of skeletal muscle: telomeric and free radical evidence that physical inactivity is responsible and not age. Clin. Sci. (Lond). 127 (6), 415–421. doi:10.1042/cs20140051
Verbanck M., Chen C. Y., Neale B., Do R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50 (5), 693–698. doi:10.1038/s41588-018-0099-7
Zee R. Y., Castonguay A. J., Barton N. S., Germer S., Martin M. (2010). Mean leukocyte telomere length shortening and type 2 diabetes mellitus: a case-control study. Transl. Res. 155 (4), 166–169. doi:10.1016/j.trsl.2009.09.012
Keywords: leukocyte telomere length, lung function, FEV1, mendelian randomization, FVC
Citation: Zhu S, Zheng W, Rao D, Tang Z and Liao X (2024) Leukocyte telomere length and lung function: a mendelian randomization study in European population. Front. Physiol. 15:1373064. doi: 10.3389/fphys.2024.1373064
Received: 19 January 2024; Accepted: 10 October 2024;
Published: 24 October 2024.
Edited by:
Yasumasa Okada, Murayama Medical Center (NHO), JapanReviewed by:
Huibin Tang, Stanford University, United StatesThiago Santos Rosa, Catholic University of Brasilia (UCB), Brazil
Copyright © 2024 Zhu, Zheng, Rao, Tang and Liao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinhui Liao, 408699371@qq.com