Baoping Hu1,2,3,4Yuhe Wang1,2,3,5Zhongtao Wang6Xue He1,2,3
Li Wang1,2,3Dongya Yuan1,2,3Yongjun He1,2,3
Tianbo Jin1,2,3*Shumei He1,2,3*- 1Key Laboratory of Molecular Mechanism and Intervention Research for Plateau Diseases of Tibet Autonomous Region, School of Medicine, Xizang Minzu University, Xianyang, China
- 2Engineering Research Center of Tibetan Medicine Detection Technology, Ministry of Education, Xianyang, China
- 3Key Laboratory of High Altitude Hypoxia Environment and Life Health, School of Medicine, Xizang Minzu University, Xianyang, China
- 4Department of Anesthesiology, The Affiliated Hospital of Xizang Minzu University, Xianyang, China
- 5Department of Clinical Laboratory, The Affiliated Hospital of Xizang Minzu University, Xianyang, China
- 6Department of Infectious Diseases, General Hospital of Tibet Military Command, Lhasa, China
Tuberculosis (TB) is an important health issue in the world. Although the relation of SLC11A1 polymorphisms with TB risk has been extensively studied, it has not been reported in the northwest Chinese Han population. Therefore, this study aimed to investigate the relationships between five polymorphisms in or near the SLC11A1 gene and susceptibility to TB. The Agena MassARRAY platform was conducted for genotyping from 510 TB patients and 508 healthy controls. Odds ratios (ORs) and 95% confidence intervals (CIs) were analyzed through logistic regression adjustment age and gender to assess the relationships between polymorphisms and TB risk. Our results identified that rs7608307 was related to increased TB risk in males (CT vs. CC: OR = 1.69, 95%CI: 1.12–2.56, p = 0.013; CT-TT vs. CC: OR = 1.61, 95%CI: 1.08–2.41, p = 0.020) and age ≤41 group (CT vs. CC: OR = 1.66, 95%CI: 1.04–2.65, p = 0.035), respectively. The SNP rs13062 was associated with the TB risk both in males (p = 0.012) and age >41 group (p = 0.021). In addition, we observed that the CC genotype of rs4674301 was correlated with increased TB risk in females (p = 0.043). Our results demonstrated the relationships between polymorphisms (rs7608307, rs4674301, and rs13062) in or near the SLC11A1 gene and age- and sex-specific TB risk in the northwest Chinese Han population.
Introduction
Tuberculosis (TB) is a common infectious disease with high public health impact globally. TB is predominantly caused by Mycobacterium tuberculosis (MTB); approximately 10% of MTB infected individuals develop clinical disease during their lifetime (Zumla et al., 2013). An estimated 1.7 billion (23%) of the world’s population is infected with MTB, leading to more than 10 million new TB cases each year (Daley, 2019). Several different factors including air pollution, smoking, drinking, poverty, undernutrition, and diabetes may affect TB development (Silva et al., 2018; Chakaya et al., 2021; Davis and Checkley, 2021). Recently, studies have demonstrated that host genetic factors also play a crucial role in the occurrence and development in TB (Ghanavi et al., 2020). Previous genome-wide association studies (GWAS) have identified multiple TB susceptibility gene loci, such as FOXP1 (rs6786408), AGMO (rs916943), ESRRB (rs12437118), and TGM6 (rs6114027) (Grant et al., 2016; Zheng et al., 2018). However, TB susceptibility is affected by age, sex-specific, and population-specific factors (Moller and Kinnear, 2020).
The human solute carrier family 11 member 1 (SLC11A1) was formerly known as natural resistance-associated macrophage protein (NRAMP1). SLC11A1, located on the human chromosome 2q35, is an important regulator of macrophage responses to MTB (McHenry et al., 2020). Genetic variants of SLC11A1 have been reported to be significantly associated with both autoimmune and infectious disease susceptibility (Archer et al., 2015; Shahzad et al., 2022). Although numerous studies have shown that the single-nucleotide polymorphisms (SNPs) of SLC11A1 were associated with the risk of TB in different populations (Wu et al., 2013; Harishankar et al., 2018; Asante-Poku et al., 2021), whether polymorphisms in and near the SLC11A1 gene are associated with susceptibility to TB in the northwest Chinese Han population has not yet been reported.
Therefore, we performed this study including 510 TB patients and 508 controls to investigate the relationships between five polymorphisms (rs11695562, rs7608307, rs4674301, rs2695343, and rs13062) in or near the SLC11A1 gene and TB risk in the northwest Chinese Han population. These results provided insights into the molecular basis for the SLC11A1 in the occurrence and development of TB.
Materials and Methods
Subjects
A total of 1,018 subjects, including 510 TB patients and 508 unrelated controls, from the Xi’an Chest Hospital, Shaanxi province, were enrolled between August 2019 and September 2021. Inclusion criteria included the following: the TB patients were diagnosed, following the TB diagnosis criteria developed by the Chinese Ministry of Health; the controls were confirmed to have no history of TB and no abnormal chest X-ray; all participants were unrelated Chinese Han population. Exclusion criteria included the following: both TB patients and controls with familial hereditary diseases, chronic inflammatory, HIV infection, diabetes mellitus, other respiratory diseases, autoimmune diseases, tumors, transplants, and long-term use of hormones were excluded.
Ethical Approval
The present study was approved by the Ethics Committee of the Xi’an Chest Hospital (No. 2012–15) and was conducted according to the principles of the Declaration of Helsinki. All the participants signed written informed consents for blood sample collection and molecular analysis.
Sample Size Calculation
We used G*Power software to estimate the sample size of the case group and the control group by independent sample t-test, and the parameters were set as follows: Tail = two, Effect size = 0.20, α = 0.05, Power = 0.89, and Allocation ratio = 1.
Single-Nucleotide Polymorphism Selection and Genotyping
The GoldMag DNA Purification Kit (GoldMag, China) was used to isolate the genomic DNA from samples, according to the manufacturer’s instructions. DNA purity and concentration were detected using the NanoDrop 2000 spectrophotometer (Thermo Fisher, United States). We selected the SLC11A1 gene as a candidate gene for the study, according to previously published studies (Meilang et al., 2012; Harishankar et al., 2018; Shahzad et al., 2022) and randomly selected five unreported SNPs (rs11695562, rs7608307, rs4674301, rs2695343, and rs13062) in or near SLC11A1 with the global minor allele frequency (MAF) greater than 0.05 from the 1000 Genomes Project data (http://www.internationalgenome.org/). The amplification and extension primers of each SNP were designed by Agena Bioscience Assay Design Suite V2.0 software. Moreover, the SNP genotype and data management were performed by the MassARRAY iPLEX platform and Agena Bioscience TYPER 4.0 software, respectively.
Statistical Analysis
The demographic characteristics of subjects involving age and sex were performed by SPSS 20.0 software. The t-test was used to detect the difference between the two groups for continuous variables, and the chi-squared test was used to detect the difference of the categorical variables between the two groups. Hardy–Weinberg equilibrium (HWE) p-values of the control group and differences in frequencies of allele and genotype distribution between the two groups were obtained from the χ2 test. Logistic regression analysis was used to calculate the odds ratios (ORs) and 95% confidence interval (CI) values for estimating the association between SNPs and TB risk adjusting age and gender. The PLINK software (version 1.07) was used to evaluate the relationship between SCL11A1 polymorphisms and TB sensibility under four genetic models. Furthermore, the stratified analyses by the mean age and gender were also performed. In this study, p < 0.05 in all statistical tests was thought to be statistically significant.
Results
Demographics
As summarized in Table 1, the basic demographic characteristics of subjects were described. The mean age of TB patients was 41.90 ± 14.83 years old including 318 (62.35%) males, 192 (37.65%) females, 233 (45.7%) age >41, and 277 (54.3%) age ≤41. For healthy controls, the mean age was 41.14 ± 18.42 years old containing 316 (62.20%) males, 192 (37.80%) females, 213 (41.9%) age >41, and 295 (58.1%) age ≤41. No difference was observed in mean age (p = 0.469) and gender (p = 0.961) distribution. Statistical analysis suggested that the age and gender between cases and controls matched.
TABLE 1. Basic characteristics of samples.
Overall Analysis
The detailed characteristics of five SNPs (rs11695562, rs7608307, rs4674301, rs2695343, and rs13062) in or near the SLC11A1 gene are displayed in Table 2. All SNPs were in accordance with HWE (p > 0.05) among healthy controls. The allele and genotype frequency distribution of all SNPs were analyzed by the χ2 test, and the results are shown in Table 2 and 3. However, the five SNPs did not present any difference in the allele and genotype frequencies among TB patients and healthy controls (p > 0.05).
TABLE 2. Allele frequency distribution and association with the TB risk.
TABLE 3. Association between SLC11A1 polymorphisms and the risk of TB.
Stratified Analysis
We performed the stratification analysis by the gender, as shown in Table 4. Analysis results suggested that rs7608307 was significantly related to increase the TB risk in males under the co-dominant model (CT vs. CC: OR = 1.69, 95%CI: 1.12–2.56, and p = 0.013), dominant model (CT-TT vs. CC: OR = 1.61, 95%CI: 1.08–2.41, and p = 0.020), and log-additive model (OR = 1.47, 95%CI: 1.01–2.13, and p = 0.043). The SNP rs13062 was also associated with increased TB risk in males under the co-dominant model (CA vs. CC: OR = 1.52, 95%CI: 1.10–2.12, and p = 0.012). The SNPs were observed to be associated with an increased risk of TB in females under the co-dominant model (CC vs. TT: OR = 3.82, 95%CI: 1.04–4.03, and p = 0.043) and the recessive model (CC vs. TT-TC: OR = 3.85, 95%CI: 1.06–4.02, and p = 0.041).
TABLE 4. Association between SLC11A1 polymorphisms and the TB risk stratified by gender.
The results of age-stratified analysis showed that the CT genotype of rs7608307 was significantly associated with the increased TB risk compared with the CC genotype in the age ≤41 group (OR = 1.66, 95%CI: 1.04–2.65, and p = 0.035). However, the AA genotype of rs13062 was associated with the reduced risk of TB in the old group under the co-dominant genetic model (OR = 0.44, 95%CI: 0.20–0.98, and p = 0.043) and the recessive model (OR = 0.40, 95%CI: 0.18–0.87, and p = 0.021) (Table 5). However, there was no significant difference between the other SNPs and TB susceptibility (p > 0.05); all data were not shown.
TABLE 5. Association between SLC11A1 polymorphisms and the TB risk stratified by age.
Discussion
In this case-control study, we explored the relationships between five SNPs in or near the SLC11A1 gene and susceptibility to TB in the northwestern Chinese Han population. The overall analysis results did not find differences in the allele and genotype frequency distributions of these five SNPs (rs11695562, rs7608307, rs4674301, rs2695343, and rs13062) in or near the SLC11A1 gene between the TB patients and healthy control groups. However, the stratified analysis by age and gender showed that rs7608307, rs13062, and rs4674301 polymorphisms were associated with TB susceptibility (p < 0.05).
SLC11A1 was identified as a proton cation antiporter, which localizes to lysosomes or late endosome. It has been reported that SLC11A1 regulates macrophage activation and plays an important role in host innate immune response against infections (Correa et al., 2017). Numerous studies have demonstrated that SLC11A1 polymorphisms were associated with the risk of TB. A study reported that rs3731865 (SLC11A1) was associated with TB risk in African-Americans, and two SNPs (rs3731863 and rs17221959) in SLC11A1 were associated with TB risk in Caucasians (Velez et al., 2009). A significant association was found between the polymorphisms (Asn543Asp and rs17235409) in the 3′ untranslated region (UTR) of the NRAMP1 gene and TB in Venezuelan (Fernandez-Mestre et al., 2015) and India population (Medapati et al., 2017). In addition, the results of three meta-analysis studies showed that SLC11A1 polymorphisms ((D543N, 3′UTR TGTG ins/del, INT4, [GT]n), rs3731865, and rs17235416)) were significantly associated with TB risk (Li et al., 2011; Meilang et al., 2012; Archer et al., 2015). Recently, rs17235409 (SLC11A1) has been observed with an even stronger protective effect against Mtb infection in Ghana (Asante-Poku et al., 2021). SLC11A1 genetic variation and low expressions have been reported to cause immune response impairment in TB patients (Shahzad et al., 2022).
Several studies about the associations between SLC11A1 polymorphisms and TB risk have also been reported in the Chinese populations. There was a statistical association between 3′UTR polymorphism in NRAMP1 (TGTG deletion) and the risk of TB in the Chinese Kazak population (Wu et al., 2013). The D543NG/A and 3′UTR TGTG+/del were related to TB susceptibility in the Chinese Han population (Wu et al., 2015). The rs17235409 (SLC11A1) was associated with the spinal TB risk in the southern Han Chinese population, and the NRAMP1 protein expression was increased in spinal TB patients (Li et al., 2022). In our study, we only observed that the CT and CT-TT genotypes of rs7608307 and the CA genotype of rs13062 were associated with the increased risk of TB in the younger group (age ≤41) and males, respectively, and the CC genotype of rs4674301 was correlated with the increased TB risk in females. Similar to our results, D543N in SLC11A1 was correlated with the TB risk in age ≤65 years group and the females of the Hong Kong Chinese population (Leung et al., 2007). These findings may suggest the genetic susceptibility to TB differences by age and sex and emphasize the importance of considering heterogeneity in genetic and TB association studies. TB epidemiology is a sex-specific disease characterized by significant differences in prevalence between men and women worldwide (Hertz and Schneider, 2019). Sex differences in TB may be due to genetics, sex hormones, lifestyle, and other factors. Previous research found that the odds of Mtb infection increased significantly with age in young girls, and the increase in odds was borderline significant (Fernandes et al., 2018). Aging has significant effects on both the innate and adaptive immune systems, which may contribute to the increased risk of TB (Byng-Maddick and Noursadeghi, 2016). Our results are needed to verify in further studies with large samples.
There are some limitations to this study that cannot be ignored. First, this study found for the first time that rs7608307, rs4674301, and rs13062 were associated with age- and sex-specific susceptibility to TB in the northwestern Chinese Han population. The results of this study need to be replicated in more Chinese populations. Second, information of some risk factors is missing (e.g., smoking, alcohol consumption, living conditions, and pollution). We will consider these factors in future studies. Finally, this study is only a very preliminary study, and further experiments are needed to explore gene expressions and protein functions.
Conclusion
Our results revealed the significant association between polymorphisms (rs7608307, rs4674301, and rs13062) in or near the SLC11A1 gene with age- and sex-specific susceptibility to TB in the northwestern Chinese Han population. A larger sample size will be needed to confirm our results in further studies. The study provides an important direction to understand the occurrence and development mechanism of TB in Chinese population.
Data Availability Statement
The original contributions presented in the study are included in the article. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the Xi’an Chest Hospital (No. 2012–15). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
TJ conceived and designed the study; BH revised the article; SH and YW wrote the manuscript and analyzed the data; XH, ZW, and LW performed research; DY and YH collected and checked information. All authors have read and agreed with the final manuscript.
Funding
This work was supported by the Natural Science Foundation of Tibet Autonomous Region (No. XZ2019ZRG-27(Z) and the Natural Science Foundation of Xizang (Tibet) Autonomous Region (No. XZ2019ZRG-37(Z)).
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.
Acknowledgments
We are grateful to all participants for providing blood samples. We also appreciate the reviewers and editors for their efforts and patience.
References
Archer, N. S., Nassif, N. T., and O'Brien, B. A. (2015). Genetic Variants of SLC11A1 Are Associated with Both Autoimmune and Infectious Diseases: Systematic Review and Meta-Analysis. Genes Immun. 16 (4), 275–283. doi:10.1038/gene.2015.8
Asante-Poku, A., Morgan, P., Osei-Wusu, S., Aboagye, S. Y., Asare, P., and Otchere, I. D. (2021). Genetic Analysis of TB Susceptibility Variants in Ghana Reveals Candidate Protective Loci in SORBS2 and SCL11A1 Genes. Front. Genet. 12, 729737. doi:10.3389/fgene.2021.729737
Byng-Maddick, R., and Noursadeghi, M. (2016). Does Tuberculosis Threaten Our Ageing Populations? BMC Infect. Dis. 16, 119. doi:10.1186/s12879-016-1451-0
Chakaya, J., Khan, M., Ntoumi, F., Aklillu, E., Fatima, R., Mwaba, P., et al. (2021). Global Tuberculosis Report 2020 - Reflections on the Global TB Burden, Treatment and Prevention Efforts. Int. J. Infect. Dis. 113 (Suppl. 1), S7–S12. doi:10.1016/j.ijid.2021.02.107
Correa, M. A., Canhamero, T., Borrego, A., Katz, I. S. S., Jensen, J. R., Guerra, J. L., et al. (2017). Slc11a1 (Nramp-1) Gene Modulates Immune-Inflammation Genes in Macrophages during Pristane-Induced Arthritis in Mice. Inflamm. Res. 66 (11), 969–980. doi:10.1007/s00011-017-1077-8
Daley, C. L. (2019). The Global Fight against Tuberculosis. Thorac. Surg. Clin. 29 (1), 19–25. doi:10.1016/j.thorsurg.2018.09.010
Davis, J. L., and Checkley, W. (2021). Characterization of Air Pollution Exposures as Risk Factors for Tuberculosis Infection. Am. J. Respir. Crit. Care Med. 204 (10), 1130–1131. doi:10.1164/rccm.202107-1795ED
Fernandes, P., Ma, Y., Gaeddert, M., Tsacogianis, T., Marques-Rodrigues, P., Fregona, G., et al. (2018). Sex and Age Differences in Mycobacterium tuberculosis Infection in Brazil. Epidemiol. Infect. 146 (12), 1503–1510. doi:10.1017/S0950268818001450
Fernandez-Mestre, M., Villasmil, A., Takiff, H., and Alcala, Z. F. (2015). NRAMP1 and VDR Gene Polymorphisms in Susceptibility to Tuberculosis in Venezuelan Population. Dis. Markers 2015, 860628. doi:10.1155/2015/860628
Ghanavi, J., Farnia, P., Farnia, P., and Velayati, A. A. (2020). Human Genetic Background in Susceptibility to Tuberculosis. Int. J. Mycobacteriol 9 (3), 239–247. doi:10.4103/ijmy.ijmy_118_20
Grant, A. V., Sabri, A., Abid, A., Abderrahmani Rhorfi, I., Benkirane, M., Souhi, H., et al. (2016). A Genome-wide Association Study of Pulmonary Tuberculosis in Morocco. Hum. Genet. 135 (3), 299–307. doi:10.1007/s00439-016-1633-2
Harishankar, M., Selvaraj, P., and Bethunaickan, R. (2018). Influence of Genetic Polymorphism towards Pulmonary Tuberculosis Susceptibility. Front. Med. (Lausanne) 5, 213. doi:10.3389/fmed.2018.00213
Hertz, D., and Schneider, B. (2019). Sex Differences in Tuberculosis. Semin. Immunopathol. 41 (2), 225–237. doi:10.1007/s00281-018-0725-6
Leung, K. H., Yip, S. P., Wong, W. S., Yiu, L. S., Chan, K. K., Lai, W. M., et al. (2007). Sex- and Age-dependent Association of SLC11A1 Polymorphisms with Tuberculosis in Chinese: a Case Control Study. BMC Infect. Dis. 7, 19. doi:10.1186/1471-2334-7-19
Li, T., Wang, L., Guo, C., Zhang, H., Xu, P., Liu, S., et al. (2022). Polymorphisms of SLC11A1(NRAMP1) Rs17235409 Associated with and Susceptibility to Spinal Tuberculosis in a Southern Han Chinese Population. Infect. Genet. Evol. 98, 105202. doi:10.1016/j.meegid.2021.105202
Li, X., Yang, Y., Zhou, F., Zhang, Y., Lu, H., Jin, Q., et al. (2011). SLC11A1 (NRAMP1) Polymorphisms and Tuberculosis Susceptibility: Updated Systematic Review and Meta-Analysis. PLoS One 6 (1), e15831. doi:10.1371/journal.pone.0015831
McHenry, M. L., Bartlett, J., Igo, R. P., Wampande, E. M., Benchek, P., Mayanja-Kizza, H., et al. (2020). Interaction between Host Genes and Mycobacterium tuberculosis Lineage Can Affect Tuberculosis Severity: Evidence for Coevolution? PLoS Genet. 16 (4), e1008728. doi:10.1371/journal.pgen.1008728
Medapati, R. V., Suvvari, S., Godi, S., and Gangisetti, P. (2017). NRAMP1 and VDR Gene Polymorphisms in Susceptibility to Pulmonary Tuberculosis Among Andhra Pradesh Population in India: a Case-Control Study. BMC Pulm. Med. 17 (1), 89. doi:10.1186/s12890-017-0431-5
Meilang, Q., Zhang, Y., Zhang, J., Zhao, Y., Tian, C., Huang, J., et al. (2012). Polymorphisms in the SLC11A1 Gene and Tuberculosis Risk: a Meta-Analysis Update. Int. J. Tuberc. Lung Dis. 16 (4), 437–446. doi:10.5588/ijtld.10.0743
Moller, M., and Kinnear, C. J. (2020). Human Global and Population-specific Genetic Susceptibility to Mycobacterium tuberculosis Infection and Disease. Curr. Opin. Pulm. Med. 26 (3), 302–310. doi:10.1097/MCP.0000000000000672
Shahzad, F., Bashir, N., Ali, A., Nadeem, A., Ammar, A., Kashif, M., et al. (2022). SLC11A1 Genetic Variation and Low Expression May Cause Immune Response Impairment in TB Patients. Genes Immun. 23 (2), 85–92. doi:10.1038/s41435-022-00165-9
Silva, D. R., Munoz-Torrico, M., Duarte, R., Galvao, T., Bonini, E. H., Arbex, F. F., et al. (2018). Risk Factors for Tuberculosis: Diabetes, Smoking, Alcohol Use, and the Use of Other Drugs. J. Bras. Pneumol. 44 (2), 145–152. doi:10.1590/s1806-37562017000000443
Velez, D., Hulme, W., Myers, J., Stryjewski, M., Abbate, E., Estevan, R., et al. (2009). Association of SLC11A1 with Tuberculosis and Interactions with. NOS2A TLR2 African-Americans Caucasians 13 (9), 1068–1076.
Wu, F., Zhang, W., Zhang, L., Wu, J., Li, C., Meng, X., et al. (20132013). NRAMP1, VDR, HLA-DRB1, and HLA-DQB1 Gene Polymorphisms in Susceptibility to Tuberculosis Among the Chinese Kazakh Population: a Case-Control Study. Biomed. Res. Int., 484535. doi:10.1155/2013/484535
Wu, L., Deng, H., Zheng, Y., Mansjo, M., Zheng, X., Hu, Y., et al. (2015). An Association Study of NRAMP1, VDR, MBL and Their Interaction with the Susceptibility to Tuberculosis in a Chinese Population. Int. J. Infect. Dis. 38, 129–135. doi:10.1016/j.ijid.2015.08.003
Zheng, R., Li, Z., He, F., Liu, H., Chen, J., Chen, J., et al. (2018). Genome-wide Association Study Identifies Two Risk Loci for Tuberculosis in Han Chinese. Nat. Commun. 9 (1), 4072. doi:10.1038/s41467-018-06539-w
Keywords: tuberculosis, SLC11A1, polymorphism, susceptibility, case-control
Citation: Hu B, Wang Y, Wang Z, He X, Wang L, Yuan D, He Y, Jin T and He S (2022) Association of SLC11A1 Polymorphisms With Tuberculosis Susceptibility in the Chinese Han Population. Front. Genet. 13:899124. doi: 10.3389/fgene.2022.899124
Received: 18 March 2022; Accepted: 02 May 2022;
Published: 22 July 2022.
Edited by:
Ramcés Falfán-Valencia, Instituto Nacional de Enfermedades Respiratorias-México (INER), MexicoReviewed by:
Khuthala Mnika, University of Cape Town, South AfricaKarol J. Nava-Quiroz, National Autonomous University of Mexico, Mexico
Tjip S. van der Werf, University of Groningen, Netherlands
Marco Antonio Ponce-Gallegos, Civil Hospital of Guadalajara, Mexico
Copyright © 2022 Hu, Wang, Wang, He, Wang, Yuan, He, Jin and He. 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: Tianbo Jin, jintianbo@gmail.com, Shumei He, 573447802@qq.com