- 1National Clinical Research Center for Obstetric and Gynecologic Diseases, Department of Obstetrics and Gynecology, Peking Union Medical College Hospital, Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China
- 2Clinical Research Laboratory, Peking Union Medical College Hospital, Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China
- 3Center for Rare Diseases Research, Chinese Academy of Medical Sciences, Beijing, China
Background: Complete androgen insensitivity syndrome (CAIS, OMIM; 300068) is a disorder of sex development with X-linked recessive inheritance. Cases of CAIS usually present as female phenotype, with primary amenorrhea and/or inguinal hernia. Family aggregation is a rare scenario.
Methods: This study is a retrospective analysis of CAIS cases in a three-generation pedigree. The patients’ genomes were determined by sequencing the androgen receptor (AR) gene. The clinical data of the patients, including manifestations, hormone levels, and AR variants, were analyzed.
Results: Sixteen people in this family were involved. A deletion variant (c.1847_1849del; p. Arg616del) was identified in exon 3 of AR, which encodes the DNA binding domain. Until now, four patients and four carriers have been identified in three generations of this family. All the patients live as female, and one has developed gonadal malignancy.
Conclusion: The present study identified a deletion variant in three generations of a family with CAIS, including four carriers and four patients. This study verified the genetic pattern and the corresponding clinical characteristics of CAIS. Furthermore, a case with gonadal malignancy was discovered. The information on diagnosis and treatment in this pedigree is useful for prenatal diagnosis and genetic counseling of similar families.
1 Introduction
Complete androgen insensitivity syndrome (CAIS) is a rare X-linked recessive hereditary disorder, which affects the sexual development of XY embryos through complete androgen resistance (Batista et al., 2018). CAIS involves mutations in the Xq11-q12 region of the androgen receptor (AR) gene (Radpour et al., 2013). The effect of androgens depends mainly on their direct interactions with AR encoded by the AR gene (Verhoeven and Swinnen, 1999), which contains eight exons and encodes 920 amino acid residues (Chamberlain et al., 1996; Sack et al., 2001). AR is a single-strand polypeptide composed of four functional domains: the N-terminal domain (NTD), the DNA-binding domain (DBD), the hinge domain, and the C-terminal ligand-binding domain (LBD) (Tsai and O'Malley, 1994). AR mutations are identified in >95% of patients with CAIS, 70% of which are maternally inherited mutations; the remaining 30% are de novo mutations (Oakes et al., 2008; Zhang et al., 2021).
Androgens are responsible for the proper development of internal and external genitalia during the fetal stage (Alemany, 2022) and promote male secondary sexual characteristics in adolescence and stimulate muscle and bone growth, spermatogenesis, and fertility in adults (Yamada et al., 2006). The most common presentation of CAIS is the occurrence of unilateral or bilateral inguinal hernias in infants or children; however, some individuals are not diagnosed until adolescence due to primary amenorrhea.
In this study, we report a variant in the DBD region of AR across three generations from a single Chinese pedigree with CAIS. A hemizygous three-base deletion was identified in exon 3 of AR (GenBank NM_000044.6: c.1847_1849del), leading to the deletion of Arg616. This study also discusses the genetic variant, hereditary pattern, and genetic consultation in this family.
2 Patients and methods
2.1 Family and patients
The clinical data of this CAIS pedigree were retrospectively collected from medical records and telephone inquiries (Figure 1). The study was approved by the Ethics Committee of Peking Union Medical College Hospital (PUMCH) (IRB Number: JS-2510) and written consent to participate and consent to publication were obtained from each individual or legal guardian if underage. The data used in this study were anonymous or kept confidential. This study did not affect the diagnosis or treatment of the patients.
FIGURE 1. Pedigree analysis of the Chinese family affected by CAIS. The genotypic 46, XX and 46, XY individuals are represented by circles and squares, respectively. Black squares, affected individuals; semi-black circles, carriers; arrows, proband (II-8) in this pedigree.
2.2 Clinical data
Clinical information was summarized, including the chief complaints and signs related to sexual development, serum hormone levels, surgical information, and pathological results. We used the Tanner stage to describe the development of breasts or pubic or axillary hair. Patient height, weight, and body mass index (BMI) were documented at the first visit. The hormone levels, including serum follicle-stimulating hormone (FSH), luteinizing hormone (LH), testosterone (T), estradiol (E2), and progesterone (P), were measured on an automated Elecsys Immunoanalyzer (Beckmann, United States) system.
2.3 AR analysis
Targeted Sanger sequencing of AR was performed to evaluate the variant. Peripheral blood samples from family members including the proband (II-8), her sisters (II-2, II-3, II-4, II-6, and II-7), and her parents (I-1 and I-2) were sent for AR testing. Additionally, blood samples from the husbands and daughters of her sisters (II-1, II-5, III-1, III-2, III-3, III-4, III-5, and III-6) were also tested for AR (Figure 1).
The Lab-Aid 820 automatic DNA extraction kit (Xiamen Zhishan Biotechnology Co., Ltd., China) was used to extract genomic DNA from the peripheral blood samples. The eight exons and intron-exon boundaries of AR (GenBank accession number NM_000044) were amplified by polymerase chain reaction (PCR) with primers designed by Primer3 online software (Table 1). The PCR volume was 25 μL, containing 20–100 ng genomic DNA, 12.5 μL PCR master mix (Beijing Tianyihuiyuan Biological Company), and 10 μM of each primer. The regions of candidate SNVs were amplified according to the manufacturer’s instructions.
For Sanger sequencing, the extracted genomic DNA was sequenced with the BigDye Terminator Cycle Sequencing Ready Reaction Kit, version 3.1 (Applied Biosystems; Thermo Fisher Scientific, Inc., Waltham, MA, United States) according to the manufacturer’s instructions. The sequences were analyzed on an ABI 3130 Genetic Analyzer (http://tools.thermofisher.com/content/sfs/manuals/4477796.pdf). The sequencing results were compared with the reference genomic sequence obtained from the UCSC Genome Browser (Santa Cruz, CA, United States; https://genome.ucsc.edu/). Variants were confirmed based on the Human Gene Mutation Database (HGMD) and the NCBI dbSNP database (http://www.ncbi.nlm.nih.gov/SNP).
3 Results
3.1 Clinical manifestations and hormone profiles
Four patients were identified in this family. Their chromosomal analysis confirmed a 46, XY karyotype. II-8 was the proband in this pedigree. The clinical manifestations and hormone profiles of the four patients (II-8, III-2, III-4, and III-1) are presented in Table 2.
The proband (II-8, 26 years) was the first from this pedigree to visit our hospital after 11 years of inguinal hernia and primary amenorrhea. Gonadal malignancy was identified during further treatment. Her niece (III-2, 16 years) developed the same symptom of primary amenorrhea. The ultrasound at the local hospital showed a congenital absence of a uterus. For further treatment, she was admitted to undergo an operation. After their diagnosis of CAIS, her 8-year-old younger niece (III-4) presented to outpatient clinics with findings of inguinal masses for 1 year. Additionally, a 3 cm mass protruding from the right groin after walking or running was observed in her 2-year-old niece (III-1). She was diagnosed genetically as a patient of CAIS without any secondary sexual developmental abnormalities.
Except for the 2-year-old girl without physical examination, the other three patients all presented a blind-ended vagina and lack of cervix. The absence of a uterus was confirmed by ultrasound in the four patients. Regarding secondary sexual characteristics, neither of the two postpubertal patients (II-8 and III-2) showed well-developed pubic or axillary hair, and Tanner stage IV breast development that was consistent with a typical CAIS, while the prepubertal patients (III-4 and III-1) were consistent with their age.
Preoperative sex hormone testing showed serum levels of E2 and P in II-8 and III-2 in agreement with the lower early follicular phase. The T levels were elevated to the normal male range or significantly higher than the female average. The gonadotropin level was normal in III-2 but was increased significantly to the postmenopausal level in II-8. III-4 and III-1 presented with normal sex hormone levels consistent with their age.
The proband (II-8) developed a malignancy of the gonads and received cytoreductive surgery. Pathology confirmed a yolk sac tumor, seminoma, and testis. The postoperative pathological stage was IIIc according to the National Comprehensive Cancer Network guidelines. After surgery and four courses of chemotherapy with a combination of bleomycin, cisplatin, and etoposide, the cancer antigen 125 (CA125) level plummeted from 455 U/ml to 6.5 U/ml, and the alpha-fetoprotein (AFP) level from 804 ng/ml to 2 ng/ml. Proper hormone replacement therapy was administered. III-2 had received laparoscopic bilateral gonadectomy, in which the histological analysis demonstrated underdeveloped testes of both gonads.
3.2 Genetic analysis of AR
In this family, four hemizygous mutations (II-8, III-1, III-2, and III-4) and four heterozygous mutations (I-2, II-2, II-6, and III-5) were identified (Supplementary Figure S1). Pedigree analysis showed that the variant was from matrilineal inheritance (Figure 1).
In this family, a deletion variant (c.1847_1849del; p. Arg616del) was identified in exon 3 of AR, which encodes the DBD, and has been previously reported (Figure 2). The c.1847_1849del (p. Arg616del) is a known variant included in the HGMD database and had been detected in an androgen-insensitive pedigree (PMID: 8162033). This variant causes a change in protein length (PM4). This variant was not detected in the Genome Aggregation Database (gnomAD) normal population database (PM2_Supporting). Based on available evidence, according to the 2015 American College of Medical Genetics and Genomics (ACMG) guideline, this variant was defined as a variant of unknown clinical significance (PM4+PM2_Supporting). Although this variant was defined as a variant of unknown clinical significance, it caused clinical symptoms in our study, leading to the occurrence of CAIS in this family.
FIGURE 2. Sanger sequencing of AR in the family. (A) Sanger sequencing confirmed that the proband (II-8) and her nieces (III-1, III-2, and III-4) have the c.1847_1849del (p. Arg616del) variant. (B) The proband’s mother (I-2), sisters (II-2 and II-6), and another niece (III-5) are heterozygous carriers, (C) while her father (I-1) and other sisters (II-3, II-4, and II-7) are wild type.
4 Discussion
This study describes three generations of a family with CAIS, including 16 people and four patients. The chief complaints to the medical service were groin mass or primary amenorrhea. The four patients in this family were raised as female due to the lack of male phenotype. Since three members of this family visited the hospital successively and there was a phenomenon of family aggregation, the whole AR gene of the family members was sequenced to explore the pathogenesis. Sanger sequencing showed an AR deletion variant c.1847_1849del in this family, which was first reported in an androgen-insensitive family in 1994 (Beitel et al., 1994). The X-linked recessive inheritance pattern of CAIS was well-validated in this family. As the carrier of the variant, I-2 passed it on to the next generation. This pedigree was valuable and it was rare for the grandparents (I-2) to have raised six daughters in the era of family planning policy in China in the 1970s. Clinical and genetic data from the whole pedigree were collected abundantly and thoroughly. Moreover, AR was sequenced simultaneously in all blood samples of 16 people of this pedigree.
Arg616 is a strictly conserved residue in the DBD in all members of the superfamily, such as steroid, thyroid, and vitamin D receptor proteins. Beitel et al. (1994) reported that Arg616 is an invariant residue in the C-terminal a-helix of the DBD. They predicted that Arg616 contributes to the normal tetrahedral arrangement of the four cysteines at the base of the C-terminal zinc finger, and to its associated α-helices. Therefore, the Arg616 deletion causes CAIS by perturbing the normal androgen receptor-androgen response element (AR-ARE) interaction. Our results showed that the Arg616 deletion mutation in the DBD region, exon 3, is the pathogenic cause of CAIS in this Chinese family. However, functional verification in cell and animal experiments is lacking, and the genotype-phenotype correlation requires further study.
Our genetic analysis revealed that the patients inherited the variant from their maternal line. The 46XY progeny carrying the mutated gene can develop CAIS. As an X-linked recessive genetic disease, the three-level prevention of CAIS is advocated. The primary prevention is carrier screening. It is important to identify heterozygous carriers of gene mutations for genetic counseling. AR testing is recommended for women with a family history of CAIS. In this study, we revealed the existence of the exon 3 p. Arg616del variant in this family. Therefore, it would be more economical to assess only exon 3 instead of sequencing the entire AR gene in future genetic tests in this family. The secondary prevention is prenatal diagnosis. The diagnosis can be confirmed or excluded by amniocentesis, villus sampling for fetal chromosome karyotype analysis, and genetic testing. The phenotype of external genital development can be determined by imaging examinations (ultrasound and magnetic resonance imaging) during pregnancy, to determine the consistency of fetal chromosomal sex and external genital phenotype. Tertiary prevention focuses on postnatal physical examinations. Unclear female vulva or genitalia sex or 46XY chromosome karyotype, combined with CAIS family history, can be considered for the diagnosis of CAIS (罗敏 and 姚凤霞,2021). Further genetic testing can be performed to determine the cause, providing the possibility of genetic counseling for future pregnancies.
II-8 had experienced gonadal malignancy (testicular germ cell tumor), which was first identified in patients with CAIS with this p. Arg616del variant. Chung et al. reported that gonadal malignant transformation was related to a variety of factors, including testicular location, individual genetic susceptibility, and residual androgen activity (Chung et al., 2013). Oncogenetic susceptibility may also be associated with one or more single nucleotide polymorphisms (SNPs) (Kratz et al., 2011). In this study, all patients carried the same variant and had ectopic gonads, but thus far only II-8 had developed gonadal malignancy. We speculated that the tumor was not necessarily associated with the variant; rather, it may be related to the gonad location, residual androgen activity, or the postponement of gonadectomy. The literature indicates that the risk of malignant progression increases with age (Cools et al., 2005; Cools et al., 2006). Therefore, testes excision is of great significance for this family of CAIS cases.
The timing of gonadectomy in patients with CAIS remains controversial. If prophylactic orchiectomy is performed, the occurrence of testicular germ cell tumor (TGCT) in patients with CAIS can be reduced. But the reservation of testicular tissue produces androgens and can be converted to estrogen by aromatase to activate pubertal development (Lanciotti et al., 2019). The reported incidence of TGCT is >22% in adulthood (Deans et al., 2012), while the incidence in childhood and adolescence is very low. In the Morris and Mahesh study of 181 cases of AIS, among patients with gonadal tumors, only one was a teenager and two were patients in their 20s (Hurt et al., 1989). Most malignancies in patients with CAIS after puberty are less invasive (Tack et al., 2018). Therefore, bilateral inguinal hernia in childhood is no longer an indication for early gonadectomy (Döhnert et al., 2017). Most patients who wish to live as female can postpone bilateral gonadectomy to post-puberty after pubertal breast development, followed by estrogen replacement therapy. However, early diagnosis and regular follow-up are still required. Ultrasound remains the first-line assessment method for gonads (Kim et al., 2007). From adolescence, annual ultrasound follow-up is recommended (Looijenga and Cools, 2017). Some classic serum markers (β-HCG, AFP, and LDH) and hormonal assessments (FSH, LH, T, and inhibin B) can also be used as a follow-up in patients with CAIS (Cools et al., 2017). Although current treatment advances advocate postponing gonadectomy until after puberty, in this family, due to the occurrence of gonadal malignancy, we advocate gonadectomy relatively earlier to prevent such malignancy.
5 Conclusion
An X-linked inheritance variant in AR was identified in a three-generation family involving 16 people. The variant c.1847_1849del (p. Arg616del) was highly conserved and was verified to perturb normal androgen–AR interactions. This variant severely disrupted androgen activity, with typical clinical characteristics of CAIS. Furthermore, a case with gonadal malignancy was diagnosed. These findings may be useful for the prenatal diagnosis and genetic counseling of the offspring of this pedigree or relevant families with CAIS.
Data availability statement
The datasets for this article are not publicly available due to concerns regarding participant/patient anonymity. Requests to access the datasets should be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Peking Union Medical College Hospital (PUMCH). Written informed consent to participate in this study was provided by the participant’s legal guardian/next of kin. Written informed consent was obtained from the individual(s), and minor(s) legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author contributions
Conceptualization: QT and LD. Data curation: FY, ML, and SD. Data collection: QT, LD, and DZ. Funding acquisition: QT. Investigation: FY. Methodology: DZ and LD. Writing–original draft: LD and DZ. Writing–review and editing: QT. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (grant number 81671424).
Acknowledgments
The authors are grateful to the families who participated in this study.
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/fgene.2023.1140083/full#supplementary-material
References
Alemany, M. (2022). The roles of androgens in humans: Biology, metabolic regulation and health. Int. J. Mol. Sci. 23 (19), 11952. doi:10.3390/ijms231911952
Batista, R., Costa, E., Rodrigues, A., Gomes, N., Faria, J., Nishi, M., et al. (2018). Androgen insensitivity syndrome: A review. Archives Endocrinol. metabolism 62 (2), 227–235. doi:10.20945/2359-3997000000031
Beitel, L. K., Prior, L., Vasiliou, D. M., Gottlieb, B., Kaufman, M., Lumbroso, R., et al. (1994). Complete androgen insensitivity due to mutations in the probable alpha-helical segments of the DNA-binding domain in the human androgen receptor. Hum. Mol. Genet. 3 (1), 21–27. doi:10.1093/hmg/3.1.21
Chamberlain, N., Whitacre, D., and Miesfeld, R. (1996). Delineation of two distinct type 1 activation functions in the androgen receptor amino-terminal domain. J. Biol. Chem. 271 (43), 26772–26778. doi:10.1074/jbc.271.43.26772
Chung, C., Kanetsky, P., Wang, Z., Hildebrandt, M., Koster, R., Skotheim, R., et al. (2013). Meta-analysis identifies four new loci associated with testicular germ cell tumor. Nat. Genet. 45 (6), 680–685. doi:10.1038/ng.2634
Cools, M., Drop, S., Wolffenbuttel, K., Oosterhuis, J., and Looijenga, L. (2006). Germ cell tumors in the intersex gonad: Old paths, new directions, moving frontiers. Endocr. Rev. 27 (5), 468–484. doi:10.1210/er.2006-0005
Cools, M., van Aerde, K., Kersemaekers, A., Boter, M., Drop, S., Wolffenbuttel, K., et al. (2005). Morphological and immunohistochemical differences between gonadal maturation delay and early germ cell neoplasia in patients with undervirilization syndromes. J. Clin. Endocrinol. metabolism 90 (9), 5295–5303. doi:10.1210/jc.2005-0139
Cools, M., Wolffenbuttel, K., Hersmus, R., Mendonca, B., Kaprová, J., Drop, S., et al. (2017). Malignant testicular germ cell tumors in postpubertal individuals with androgen insensitivity: Prevalence, pathology and relevance of single nucleotide polymorphism-based susceptibility profiling. Hum. Reprod. Oxf. Engl. 32 (12), 2561–2573. doi:10.1093/humrep/dex300
Deans, R., Creighton, S., Liao, L., and Conway, G. (2012). Timing of gonadectomy in adult women with complete androgen insensitivity syndrome (CAIS): Patient preferences and clinical evidence. Clin. Endocrinol. 76 (6), 894–898. doi:10.1111/j.1365-2265.2012.04330.x
Döhnert, U., Wünsch, L., and Hiort, O. (2017). Gonadectomy in complete androgen insensitivity syndrome: Why and when? Sexual development: Genetics, molecular biology, evolution, endocrinology, embryology, and pathology of sex determination and differentiation. Sex. Dev. 11 (4), 171–174. doi:10.1159/000478082
Hurt, W., Bodurtha, J., McCall, J., and Ali, M. (1989). Seminoma in pubertal patient with androgen insensitivity syndrome. Am. J. obstetrics Gynecol. 161 (3), 530–531. doi:10.1016/0002-9378(89)90350-5
Kim, W., Rosen, M. A., Langer, J. E., Banner, M. P., Siegelman, E. S., and Ramchandani, P. (2007). US–MR imaging correlation in pathologic conditions of the scrotum. RadioGraphics 27 (5), 1239–1253. doi:10.1148/rg.275065172
Kratz, C., Greene, M., Bratslavsky, G., and Shi, J. (2011). A stratified genetic risk assessment for testicular cancer. Int. J. Androl. 34, e98–e102. doi:10.1111/j.1365-2605.2011.01156.x
Lanciotti, L., Cofini, M., Leonardi, A., Bertozzi, M., Penta, L., and Esposito, S. (2019). Different clinical presentations and management in complete androgen insensitivity syndrome (CAIS). Int. J. Environ. Res. public health 16 (7), 1268. doi:10.3390/ijerph16071268
Looijenga, L., and Cools, M. (2017). Update on the pathophysiology and risk factors for the development of malignant testicular germ cell tumors in complete androgen insensitivity syndrome. Sexual development: Genetics, molecular biology, evolution, endocrinology, embryology, and pathology of sex determination and differentiation. Sex. Dev. 11 (4), 175–181. doi:10.1159/000477921
Oakes, M., Eyvazzadeh, A., Quint, E., and Smith, Y. (2008). Complete androgen insensitivity syndrome--a review. J. Pediatr. Adolesc. Gynecol. 21 (6), 305–310. doi:10.1016/j.jpag.2007.09.006
Radpour, R., Falah, M., Aslani, A., Zhong, X. Y., and Saleki, A. (2013). Identification of a critical novel mutation in the exon 1 of androgen receptor gene in 2 brothers with complete androgen insensitivity syndrome. J. Androl. 30 (3), 230–232. doi:10.2164/jandrol.108.005520
Sack, J., Kish, K., Wang, C., Attar, R., Kiefer, S., An, Y., et al. (2001). Crystallographic structures of the ligand-binding domains of the androgen receptor and its T877A mutant complexed with the natural agonist dihydrotestosterone. Proc. Natl. Acad. Sci. U. S. A. 98 (9), 4904–4909. doi:10.1073/pnas.081565498
Tack, L., Maris, E., Looijenga, L., Hannema, S., Audi, L., Köhler, B., et al. (2018). Management of gonads in adults with androgen insensitivity: An international survey. Hormone Res. Paediatr. 90 (4), 236–246. doi:10.1159/000493645
Tsai, M., and O'Malley, B. (1994). Molecular mechanisms of action of steroid/thyroid receptor superfamily members. Annu. Rev. Biochem. 63, 451–486. doi:10.1146/annurev.bi.63.070194.002315
Verhoeven, G., and Swinnen, J. (1999). Indirect mechanisms and cascades of androgen action. Mol. Cell. Endocrinol. 151, 205–212. doi:10.1016/s0303-7207(99)00014-3
Yamada, G., Suzuki, K., Haraguchi, R., Miyagawa, S., Satoh, Y., Kamimura, M., et al. (2006). Molecular genetic cascades for external genitalia formation: An emerging organogenesis program. Dev. Dyn. official Publ. Am. Assoc. Anatomists 235 (7), 1738–1752. doi:10.1002/dvdy.20807
Zhang, D., Yao, F., Tian, T., Deng, S., Luo, M., and Tian, Q. (2021). Clinical characteristics and molecular genetics of complete androgen insensitivity syndrome patients: A series study of 30 cases from a Chinese tertiary medical center. Fertil. Steril. 115 (5), 1270–1279. doi:10.1016/j.fertnstert.2020.12.008
Keywords: androgen receptor, complete androgen insensitivity syndrome, gonadal malignancy, variant, X-linked recessive inheritance
Citation: Ding L, Zhang D, Yao F, Luo M, Deng S and Tian Q (2023) A deletion variant Arg616 of androgen receptor in a Chinese family with complete androgen insensitivity syndrome. Front. Genet. 14:1140083. doi: 10.3389/fgene.2023.1140083
Received: 08 January 2023; Accepted: 20 April 2023;
Published: 09 May 2023.
Edited by:
Jared C. Roach, Institute for Systems Biology (ISB), United StatesCopyright © 2023 Ding, Zhang, Yao, Luo, Deng and Tian. 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: Qinjie Tian, qinjietn@163.com
†These authors have contributed equally to this work