- Department of Pediatrics, The First Affiliated Hospital, Jinan University, Guangzhou, China
Na+-taurocholate cotransporting polypeptide deficiency (NTCPD) is an autosomal recessive disorder arising from biallelic SLC10A1 mutations. As a newly-described inborn error of bile acid metabolism, the epidemiology of this condition remains largely unclear in Chinese population so far. In this study, a total of 2,828 peripheral blood samples were collected from 12 cities in Guangdong, a province with the largest population in China, and the four prevalent SLC10A1 variants c.800C > T (p.Ser267Phe), c.263T > C (p.Ile88Thr), c.595A > C (p.Ser199Arg) and c.665T > C (p.Leu222Ser) were screened for by using polymerase chain reaction (PCR)- restriction fragment length polymorphism (RFLP). As a result, 663 mutated SLC10A1 alleles were detected, and the mutated allele frequency was calculated to be 11.72% (663/5,656), with a carrier frequency 20.69% (1/5) and a theoretical morbidity rate 1.37% (1/73) of NTCPD in Guangdong province. The variant c.800C > T (p.Ser267Phe) exhibited highest allele frequency among the four prevalent variants (χ2 = 1501.27, p < 0.0001) as well as higher allele frequency in the peripheral region than that within the Pearl River Delta (χ2 = 4.834, p < 0.05). The results suggested that NTCPD might be a disorder rather common in Guangdong province. The findings depicted the molecular epidemiologic features of NTCPD, providing preliminary but significant laboratory evidences for the subsequent NTCPD diagnosis and management in Guangdong population.
Introduction
Na+-Taurocholate Cotransporting polypeptide deficiency (NTCPD) is an autosomal recessive disorder affecting the hepatic uptake of bile acids caused by biallelic variants of the solute carrier family 10 member1 (SLC10A1) gene which encodes NTCP protein. The main clinical presentation of this disease is refractory hypercholanemia, transient cholestatic jaundice in infancy, and indirect hyperbilirubinemia in neonates (Deng et al., 2021; Deng, 2021). The gene SLC10A1 is located at chromosome 14q24.2, contains five exons and has a total length of 23 kb (Hagenbuch and Dawson, 2004; Shiao et al., 2000). The protein product NTCP is expressed in the sinusoidal plasma membrane of the hepatocyte, where it functions to uptake bile salts from plasma in a Na+-dependent manner, playing a crucial role in the enterohepatic circulation of bile acids (Hagenbuch and Dawson, 2004; Anwer and Stieger, 2014).
Although SLC10A1 gene was cloned by (Hagenbuch and Meier, 1994), the first case of NTCPD was just reported by (Vaz et al., 2015), and since then, increasing number of such patients have been reported (Yan et al., 2020). However, some patients were overinvestigated and intervented due to unclear etiology, even undergoing surgical operation such as exploratory laparotomy (Li et al., 2018; Deng, 2021), while pregnant women with NTCPD were sometimes misdiagnosed with intrahepatic cholestasis of pregnancy (ICP) and given cesarean section (Chen et al., 2019). Currently, the SLC10A1 genotypic and phenotypic features of NTCPD still remain far from being well-understood.
Previous studies revealed that the allele frequency of c.800C > T (p.Ser267Phe), a pathogenic variant of SLC10A1 gene (Deng et al., 2016; Ho et al., 2004), varied greatly among different ethnic groups and geographical populations. It was the most prevalent variant in South China and Vietnam, but was not found in European Americans, African Americans or Hispanics (Peng et al., 2015). Of note, in a sample of 50 Chinese Americans, the allele frequency of the variant was 7.5% (Ho et al., 2004), suggesting that NTCPD may not be rare in China. However, as a newly-described inborn error of metabolism, the epidemiology of NTCPD remains largely unclear in Chinese population.
From June 2015 to September 2021, our team diagnosed 318 NTCPD patients by SLC10A1 gene analysis in China, and the four SLC10A1 variants of c.800C > T (p.Ser267Phe), c.263T > C (p.Ile88Thr), c.595A > C (p.Ser199Arg) and c.665T > C (p.Leu222Ser) were at the top of the list, accounting for 96.3, 2.5, 0.6 and 0.4% of all mutated alleles, respectively (Deng, 2021; Deng et al., 2021). In this study, these SLC10A1 variants were screened to investigate the epidemiology of NTCPD in Guangdong, a province with the largest population in China.
Methods
Participants
According to the latest data from the seventh national Census of China’s National Bureau of Statistics, the permanent resident population in Guangdong were 126,012,510 (http://www.stats.gov.cn/tjsj/tjgb/rkpcgb/qgrkpcgb/202106/t20210628_1818822.html). In order to achieve statistical significance, the sample size was calculated to be at least 1,537 on the basis of the estimated carrier rates of 1/50 by using the online sample size calculator EPITOOLS (https://epitools.ausvet.com.au/oneproportion), with a confidence level of 95% and the precision of 0.02. The number of samples from each region was in accordance with the population distribution.
Inclusion criteria: Apparently healthy individuals regardless of their genders and ages. The research subjects in this study were 2,828 used blood samples for health examination collected from 12 different cities in Guangdong Province from April 2011 to March 2013. Among them, there were 476, 276, 189, 169 and 192 samples collected from the five cities Guangzhou, Shenzhen, Foshan, Huizhou and Zhongshan, respectively, within Pearl River Delta region. The remaining samples were from the seven peripheral cities to the Pearl River Delta, with 162, 328, 178, 147, 181, 333, and 197 samples collected from Meizhou, Heyuan, Zhanjiang, Qingyuan, Shaoguan, Yunfu, and Shantou, respectively (Figure 1).
FIGURE 1. Regional division and the distribution of the mutated SLC10A1 alleles in different cities of Guangdong province, including the cities within the Pearl River Delta (blue) and the peripheral region (yellow). The numbers of the mutated SLC10A1 alleles from different cities were presented in parentheses. This figure was generated by means of the software Microsoft PowerPoint 2019. The base map was created by incrementally assembling the outlines of the Chinese administrative regions, which could be downloaded via the URL link http://www.900ppt.com/.
Exclusion criteria: Patients with any positive symptoms or signs including jaundice or hepatosplenomegaly.
Genetic Analysis
As in previous publications, genomic DNA was purified from the peripheral blood samples, and the variants c.800C > T (p.Ser267Phe) (Deng et al., 2016), c.263T > C (p.Ile88Thr) (Qiu et al., 2017) and c.595A > C (p.Ser199Arg) (Li et al., 2019) were detected using established polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) methods. The variants c.665T > C (p. Leu222Ser) was screened using a novel PCR-RFLP procedure as below. The nucleotide sequences of the forward and reverse primers in PCR amplification were 5′-GTGCTTGGCTGAGTTTGTAATAATC-3′ and 5′- GTGTTTGGATACCTTTGGTGTCTG-3′, respectively (Invitrogen; Thermo Fisher Scientific, Inc.). The target fragment was amplified using a PCR kit (Takara Biotechnology Co., Ltd.) and the PCR thermocycling conditions were: 94°C for 5 min, followed by 35 cycles at 94°C for 30 s, 58°C for 40 s and 72°C for 50 s, and 72°C for 10 min. The TaqI restriction enzyme (Thermo Fisher Scientific, inc.) was used to digest the PCR products and the digested DNA products were subsequently separated by electrophoresis in a 4% agarose gel (Figure 2).
FIGURE 2. Novel approach developed for the screening of the SLC10A1 variant c.665T > C (p. Leu222Ser). (A) schematic diagram of the PCR-RFLP approach. The mutated SLC10A1 allele had a TaqI restriction enzyme site and produced the 599 and 485 bp fragments from the 1,084 bp band following enzymatic digestion. (B) Gel electrophoresis showed that the heterozygote of the variant c.665T > C (p. Leu222Ser) had three bands of 1084, 599, and 485 bp, while the wildtype had only one band of 1084 bp.
Calculation of the Allele Frequencies, Carrier Frequencies and NTCPD Morbidity Rate
The mutated SLC10A1 allele frequencies, carrier frequencies and NTCPD morbidity rate in different areas were calculated based on the Hardy-Weinberg equilibrium. The genotypes AA (healthy individual), AB (carrier) and BB (patient) of a biallelic genetic marker were expected to have the relative frequencies of p2, 2pq and q2, with p and q being the A (wild type) and B (mutant) allele frequency, respectively; and thus p + q = 1. The values of 2pq and q2 represented as the carrier frequency of SLC10A1 variants and NTCPD morbidity rate, respectively (Graffelman et al., 2017; Graffelman and Weir, 2018).
Statistical Analysis
By using the statistical software SPSS version 26.0 (IBM Corp., Chicago, IL, United States), the distribution of the four SLC10A1 variants was compared via Chi-square tests among different geographic regions, with a p value < 0.05 indicating statistical significance.
This study was conducted in accordance with the Declaration of Helsinki (as revised in 2013) and was approved by the Medical Ethics Committee of the First Affiliated Hospital, Jinan University, Guangdong, China (No.KY-2019-052). Participants were genotyped retrospectively using blood samples previously collected for the purpose of health examination. The data related to individual identification were anonymized during the entire study process. Therefore, informed consent was waived.
Results
Variants Screening
This study detected 544 heterozygous and 20 homozygous variants of c.800C > T (p.Ser267Phe), 62 heterozygotes of c.263T > C (p. Ile88Thr), four heterozygotes of c.595A > C (p. Ser199Arg) and three heterozygotes of c.665T > C (p. Leu222Ser), along with five compound heterozygote of c.800C > T (p. Ser267Phe) and c.263T > C (p. Ile88Thr). Therefore, a total of 663 mutated SLC10A1 alleles were detected in 5,656 independent alleles (2,828 samples), as shown in Table 1.
TABLE 1. Distribution of the four prevalent SLC10A1 variants in different regions of Guangdong Province.
Carrier Frequencies and Theoretical Morbidity Rates
The mutated SLC10A1 allele frequency in the Guangdong was calculated to be 11.72% (663/5,656), the carrier frequency, 20.69% (1/5), and the theoretical morbidity rate 1.37% (1/73), respectively. According to China’s seventh national census population data and the theoretical morbidity rate, it was estimated that there were at least 1,726,371 NTCPD patients in Guangdong province.
The mutated SLC10A1 allele frequency in the Pearl River Delta was 10.87% (283/2,604), with a carrier frequency 19.38% (1/5), and a theoretical morbidity 1.18% (1/85); in the peripheral region to Pearl River Delta, the mutated SLC10A1 allele frequency was 12.45% (380/3,052), with a carrier frequency 21.80% (1/5) and a theoretical morbidity rate 1.55% (1/65) (Table 1).
Distribution Comparison of the Prevalent SLC10A1 Variants Among Different Geographic Regions
The allele frequency of c.800C > T (p.Ser267Phe) was the highest among the four prevalent variants (Table 2). The allele frequency of the variant c.263T > C (p. Ile88Thr), c.595A > C (p. Ser199Arg) and c.665T > C (p. Leu222Ser) between the Pearl River Delta region and the periphery of the Pearl River Delta had no significant difference; However, the distribution of the variant c.800C > T (p. Ser267Phe) was significantly different between the two regions, and the variant frequency (9.45%, 246/2,604) in the periphery of the Pearl River Delta region was higher than that (11.24%, 343/3,052) in the Pearl River Delta region (χ2 = 4.834, p < 0.05) (Table 3).
TABLE 3. Comparison of the distribution of the four prevalent SLC10A1 variants between the Pearl River Delta and the peripheral region.
Discussion
For the first time, this study reported that the allele frequency of four prevalent SLC10A1 variants in Guangdong was 11.72% (663/5,656), with a carrier frequency 20.69% (1/5) and a theoretical morbidity rate 1.37% (1/73) of NTCPD in this province. According to the latest data from the seventh national population census in 2020, the population in Guangdong reached 126,012,510 (http://www.stats.gov.cn/tjsj/tjgb/rkpcgb/qgrkpcgb/202106/t20210628_1818822. html), and thus it was estimated that at least 1,726,371 NTCPD patients were distributed in Guangdong province. Although no relevant data have been reported in terms of the actual incidence or diagnosis rate of NTCPD, the findings in this study provided preliminary but significant data for the estimation of the NTCPD burden, constituting an important epidemiologic basis for the diagnosis and management of NTCPD in Guangdong, and so high a theoretical morbidity implied that a large number of NTCPD patients might have been overlooked or even misdiagnosed at least in Guangdong population.
Among the four prevalent variants in this study, c.800C > T (p.Ser267Phe) had the highest allele frequency, accounting for 89% (589/663) of the total variants. NTCP was also a functional receptor for human hepatitis B virus (HBV) (Yan et al., 2012), and this variant not only affected the NTCP residue critical for bile salts binding, but also severely impaired viral infection by HBV (Yan et al., 2014). Actually, HBV infection was rather common in Guangdong population. In 1992, the positive rate of HBsAg was as high as 17.85% in this province, while the national average was just 9.75% (Wu et al., 2009). Even in 2006, the positive rate of HBsAg reached 17.55% in Guangdong population aged 15 ∼ 59 years (Cui et al., 2012). Therefore, the high allele frequency of c.800C > T (p. Ser267Phe) in Guangdong might be a result of positive selection in response to HBV infection (Peng et al., 2015). Actually, this variant was associated with a lower incidence of acute-on-chronic liver failure (ACLF) in chronic hepatitis B (CHB) patients (Peng et al., 2015) and reduced the risk for HBV infection and disease progression in human (An et al., 2018). Moreover, the c.800C > T (p.Ser267Phe) variant exhibited relatively higher allele frequency (Table 2), which could be explained by a founder effect of this variant which had originated in a far remote ancestor in the long process of evolution.
This study showed that the allele frequency of c.800C > T (p.Ser267Phe) was higher in the periphery region than that within the Pearl River Delta region. Guangdong is one of the main inflow places of China’s migrant population. In 2010, the migrant population reached 31.28 million, mainly from Hunan, Guangxi, Sichuan, Hubei and other provinces, and nearly 90% of the migrant was distributed in Guangzhou, Shenzhen, Dongguan, Foshan and other cities in the Pearl River Delta region (Chen, 2017). Since the HBsAg carrier rate in these migrant source provinces was lower than that in Guangdong province (Cui et al., 2012; Duan, 2018), it was not surprising that the SLC10A1 variant frequency in the Pearl River Delta region was lower than that in the periphery of the Pearl River Delta region dominated by local residents. In other words, the different geographic distribution of the variant c.800C > T (p. Ser267Phe) might be attributed to the genetic flow which occurred between distinct founding populations.
Interestingly, 20 homozygotes of c.800C > T (p.Ser267Phe) and five compound heterozygotes of this variant with c.263T > C (p. Ile88Thr) were found in this study, who had no clinical signs and symptoms such as jaundice and hepatosplenomegaly. This is not strange since similar findings have been reported that, in affected adults, NTCPD only presented with laboratory abnormalities including hypercholanemia while the clinical presentation was usually negative (Deng et al., 2016; Liu et al., 2017), although this disease caused abnormal bilirubin metabolism in pediatric patients (Yan et al., 2020), and increased the risks of indirect hyperbilirubinemia in affected neonates as well as transient cholestatic jaundice, elevated liver enzymes and 25-hydroxyvitamin D (Vit D) deficiency during early infancy (Deng et al., 2021).
There were still some limitations in this paper. This study only focused on four prevalent SLC10A1 variants in populations from 11 cities in Guangdong Province, and sampling bias might exist with limited sample sizes in relevant cities. Therefore, further epidemiology study focusing on more SLC10A1 variants with larger sample size was in need.
Conclusion
NTCPD, with a theoretical morbidity rate of 1.37%, might be a disorder rather common in Guangdong province, and the prevalent SLC10A1 variant c.800C > T (p. Ser267Phe) exhibited significantly different distribution between Pearl River Delta and the peripheral regions. The findings depicted the molecular epidemiologic features of NTCPD in Guangdong population, providing preliminary but significant laboratory evidences for subsequent NTCPD diagnosis and management in this province.
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 author.
Author Contributions
The author Y-ZS designed the study and directed its implementation, including quality assurance and control. HL helped prepare the introduction, methods and the discussion sections of the text. RC helped conduct the literature review and prepare the methods sections of the text. G-ZL completed the collection of samples, DNA extraction and PCR-RFLP. Author W-XL directed the implementation of DNA extraction and PCR-RFLP. M-RY participated in DNA extraction of some samples.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81974057).
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.
References
An, P., Zeng, Z., and Winkler, C. A. (2018). The Loss-Of-Function S267F Variant in HBV Receptor NTCP Reduces Human Risk for HBV Infection and Disease Progression. J. Infect. Dis. 218, 1404–1410. doi:10.1093/infdis/jiy355
Anwer, M. S., and Stieger, B. (2014). Sodium-dependent Bile Salt Transporters of the Slc10a Transporter Family: More Than Solute Transporters. Pflugers Arch. - Eur. J. Physiol. 466, 77–89. doi:10.1007/s00424-013-1367-0
Chen, R., Deng, M., Rauf, Y.-M., Lin, G.-Z., Qiu, J.-W., Zhu, S.-Y., et al. (2019). Intrahepatic Cholestasis of Pregnancy as a Clinical Manifestation of Sodium-Taurocholate Cotransporting Polypeptide Deficiency. Tohoku J. Exp. Med. 248, 57–61. doi:10.1620/tjem.248.57
Chen, X. Z. (2017). Studies on the Urban Living Willingness of the Floating Population in Guangdong Province. [doctor’s Thesis]. Changchun: Jilin University.
Cui, F. Q., Gong, X. H., Chen, Y. S., Wang, F. Z., Zhen, H., Wu, Z. H., et al. (2012). Vaccination Progress of Hepatitis B Vaccine and Epidemiology Changes of Carrying Rate of Hepatitis B Surface Antigen by Province in China, 1992-2006. Chin J. Vaccin. Immunization 18, 6–13.
Deng, L. J. (2021). Clinical Phenotypic Characterization of Na+-Taurocholate Cotransporting Polypeptide Deficiency (NTCPD) in Humans: A Case-Control Study Based on SLC10A1 Genotyping Analysis. [doctor’s Thesis]. Guangzhou: Jinan University.
Deng, L. J., Ouyang, W. X., Liu, R., Deng, M., Qiu, J. W., Yaqub, M. R., et al. (2021). Clinical Characterization of NTCP Deficiency in Paediatric Patients : A Case‐control Study Based on SLC10A1 Genotyping Analysis. Liver Int. 41, 2720–2728. doi:10.1111/liv.15031
Deng, M., Mao, M., Guo, L., Chen, F.-P., Wen, W.-R., and Song, Y.-Z. (2016). Clinical and Molecular Study of a Pediatric Patient with Sodium Taurocholate Cotransporting Polypeptide Deficiency. Exp. Ther. Med. 12, 3294–3300. doi:10.3892/etm.2016.3752
Duan, R. B. (2018). The Changing Characteristic of Hepatitis B Virus Carrying Rate in China in the Past 25 years: A Meta-Analysis. [master’s Thesis]. [Taiyuan]: Shanxi medical University.
Graffelman, J., Jain, D., and Weir, B. (2017). A Genome-wide Study of hardy-weinberg Equilibrium with Next Generation Sequence Data. Hum. Genet. 136, 727–741. doi:10.1007/s00439-017-1786-7
Graffelman, J., and Weir, B. S. (2018). On the Testing of hardy-weinberg Proportions and equality of Allele Frequencies in Males and Females at Biallelic Genetic Markers. Genet. Epidemiol. 42, 34–48. doi:10.1002/gepi.22079
Hagenbuch, B., and Dawson, P. (2004). The Sodium Bile Salt Cotransport Family Slc10. Pflgers Archiv Eur. J. Physiol. 447, 566–570. doi:10.1007/s00424-003-1130-z
Hagenbuch, B., and Meier, P. J. (1994). Molecular Cloning, Chromosomal Localization, and Functional Characterization of a Human Liver Na+/bile Acid Cotransporter. J. Clin. Invest. 93, 1326–1331. doi:10.1172/JCI117091
Ho, R. H., Leake, B. F., Roberts, R. L., Lee, W., and Kim, R. B. (2004). Ethnicity-dependent Polymorphism in Na+-Taurocholate Cotransporting Polypeptide (SLC10A1) Reveals a Domain Critical for Bile Acid Substrate Recognition. J. Biol. Chem. 279, 7213–7222. doi:10.1074/jbc.M305782200
Li, H., Qiu, J. W., Lin, G. Z., Deng, M., Lin, W. X., Cheng, Y., et al. (2018). Clinical and Genetic Analysis of a Pediatric Patient with Sodium Taurocholate Cotransporting Polypeptide Deficiency. Zhongguo Dang Dai Er Ke Za Zhi 20, 279–284. doi:10.7499/j.issn.1008-8830.2018.04.005
Li, H., Deng, M., Guo, L., Qiu, J. W., Lin, G. Z., Long, X. L., et al. (2019). Clinical and Molecular Characterization of Four Patients with Ntcp Deficiency from Two Unrelated Families Harboring the Novel Slc10a1 Variant c.595a>c (p.Ser199arg). Mol. Med. Rep. 20, 4915–4924. doi:10.3892/mmr.2019.10763
Liu, R., Chen, C., Xia, X., Liao, Q., Wang, Q., Newcombe, P. J., et al. (2017). Homozygous p.Ser267phe in Slc10a1 Is Associated with a New Type of Hypercholanemia and Implications for Personalized Medicine. Sci. Rep. 7, 9214. doi:10.1038/s41598-017-07012-2
Peng, L., Zhao, Q., Li, Q., Li, M., Li, C., Xu, T., et al. (2015). The p.Ser267Phe Variant inSLC10A1is Associated with Resistance to Chronic Hepatitis B. Hepatology 61, 1251–1260. doi:10.1002/hep.27608
Qiu, J.-W., Deng, M., Cheng, Y., Atif, R.-M., Lin, W.-X., Guo, L., et al. (2017). Sodium Taurocholate Cotransporting Polypeptide (Ntcp) Deficiency: Identification of a Novel SLC10A1 Mutation in Two Unrelated Infants Presenting with Neonatal Indirect Hyperbilirubinemia and Remarkable Hypercholanemia. Oncotarget 8, 106598–106607. doi:10.18632/oncotarget.22503
Shiao, T., Iwahashi, M., Fortune, J., Quattrochi, L., Bowman, S., Wick, M., et al. (2000). Structural and Functional Characterization of Liver Cell-specific Activity of the Human Sodium/taurocholate Cotransporter. Genomics 69, 203–213. doi:10.1006/geno.2000.6329
Vaz, F. M., Paulusma, C. C., Huidekoper, H., de Ru, M., Lim, C., Koster, J., et al. (2015). Sodium Taurocholate Cotransporting Polypeptide (SLC10A1) Deficiency: Conjugated Hypercholanemia without a clear Clinical Phenotype. Hepatology 61, 260–267. doi:10.1002/hep.27240
Wu, C. G., Han, K., Shao, X. P., Shu, J., and Lin, Y. J. (2009). Sero-epidemiological Study on Hepatitis B Virus Infection in Guangdong Province. Chin. Prev. Med. 10, 350–354. doi:10.16506/j.1009-6639.2009.05.035
Yan, H., Peng, B., Liu, Y., Xu, G., He, W., Ren, B., et al. (2014). Viral Entry of Hepatitis B and D Viruses and Bile Salts Transportation Share Common Molecular Determinants on Sodium Taurocholate Cotransporting Polypeptide. J. Virol. 88, 3273–3284. doi:10.1128/jvi.03478-13
Yan, H., Zhong, G., Xu, G., He, W., Jing, Z., Gao, Z., et al. (2012). Sodium Taurocholate Cotransporting Polypeptide Is a Functional Receptor for Human Hepatitis B and D Virus. Elife 1, e00049. doi:10.7554/eLife.00049
Keywords: NTCP deficiency, epidemiology, SLC10A1 gene, Guangdong, prevalent variants
Citation: Li H, Chen R, Lin G-Z, Lin W-X, Yaqub M-R and Song Y-Z (2022) Molecular Epidemiology of Na+-Taurocholate Cotransporting Polypeptide Deficiency in Guangdong Province, China: A Pilot Study by Screening for Four Prevalent Variants of the Causative Gene SLC10A1. Front. Genet. 13:874379. doi: 10.3389/fgene.2022.874379
Received: 12 February 2022; Accepted: 24 March 2022;
Published: 27 April 2022.
Edited by:
Yuan Shi, Children’s Hospital of Chongqing Medical University, ChinaReviewed by:
Shaotao Tang, Huazhong University of Science and Technology, ChinaSujun Zheng, Capital Medical University, China
Copyright © 2022 Li, Chen, Lin, Lin, Yaqub and Song. 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: Yuan-Zong Song, songyuanzong@vip.tom.com
†These authors have contributed equally to this work