- 1Shanghai Clinical Research Center of Bone Diseases, Department of Osteoporosis and Bone Diseases, Shanghai Jiao Tong University Affiliated Sixth People’s Hospital, Shanghai, China
Spondylo-epi-metaphyseal dysplasia (SEMD) is a heterogeneous group of disorders with different modes of inheritance and is characterized by disproportionate or proportionate short stature. To date, more than 30 disease-causing genes have been identified, and different types of SEMD exhibit greatly overlapping clinical features, which usually complicate the diagnosis. This study was performed to expand the clinical and molecular spectrum of SEMD among Chinese subjects and to explore their potential phenotype–genotype relations. We enrolled seven families including 11 affected patients with SEMD, and their clinical, radiographic, and genetic data were carefully analyzed. All the seven probands showed different degrees of short stature, and each of them exhibited additional specific skeletal manifestations; four probands had extraosseous manifestations. X-rays of the seven probands showed common features of SEMD, including vertebral deformities, irregular shape of the epiphysis, and disorganization of the metaphysis. Seven variants were identified in TRPV4 (c.694C> T, p.Arg232Cys), COL2A1 (c.654 + 1G > C; c.3266_3268del, p.Gly1089del), CCN6 (c.396 T> G, p.Cys132Trp; c.721 T>C, p.Cys241Arg), SBDS (c.258 + 2T> C), and ACAN (c.1508C> A, p.Thr503Lys) genes, and two of them were novel. Two families with TRPV4 variants showed considerable intrafamily and interfamily heterogeneities. In addition, we reported one case of SEMD with a severe phenotype caused by ACAN gene mutation. Our study expands the phenotype and genetic spectrum of SEMD and provides evidence for the phenotype–genotype relations, aiding future molecular and clinical diagnosis as well as procreative management of SEMD.
Introduction
Spondylo-epi-metaphyseal dysplasia (SEMD) is a heterogeneous genetic disorder, involving vertebral, epiphyseal, and metaphyseal dysplasia, and is diagnosed based on clinical phenotype, radiographic examination, and molecular sequencing. The primary clinical feature is a different degree of short stature (may present with short limbs or short trunk), combined with specific orthopedic symptoms (such as developmental coxa vara and scoliosis). Epiphyseal dysplasia usually leads to early-onset osteoarthritis, mostly in weight-bearing joints (Briggs et al., 1995). Odontoid hypoplasia causes atlantoaxial instability with severe spinal compression problems (Miyoshi et al., 2004). Radiological features included platyspondyly, vertebral body irregularity, destruction of articular cartilage, and dysplasia of epiphysis and metaphysis (Spranger, 1989).
This category of diseases is usually given the nomenclature according to the site name of manifest radiographic abnormalities (Alanay and Lachman, 2011; Mortier et al., 2019) (Figure 1), including spondyloepiphyseal dysplasia (SED), spondylometaphyseal dysplasia (SMD), spondylo-epi-metaphyseal dysplasia (SEMD), metaphyseal dysplasia (MD), and epiphyseal dysplasia (ED). With the in-depth study, SEMD is further classified according to specific clinical manifestations. For example, SED with different clinical manifestations can be divided into spondyloepiphyseal dysplasia congenita (SEDC, OMIM# 183900), spondyloepiphyseal dysplasia-Kimberley type (SED-KT, OMIM# 608361), and spondyloepiphyseal dysplasia-Maroteaux type (SED-MT, OMIM# 184095).
FIGURE 1. Nomenclature according to the involved area name of dysplasia. A + D: epiphyseal dysplasia, A + E: metaphyseal dysplasia, B + D: spondyloepiphyseal dysplasia, B + E: spondylometaphyseal dysplasia, and B + D + E: spondylo-epi-metaphyseal dysplasia.
So far, more than 30 pathogenic genes have been identified to cause SEMD. These disease-causing genes are involved in encoding various types and functions of proteins (Cormier-Daire, 2008). It is difficult to define them as a common signal pathway among all the SEMD. In such clinical and genetic heterogeneous diseases, different gene–disease associations and overlapping clinical characteristics of different genes complicate the differential diagnosis. According to the nosology and classification of genetic skeletal disorders 2019 revision, groups 10–13 are classified as spondylo-epi-metaphyseal abnormalities, but these diseases still exist in other groups (Mortier et al., 2019). The more common SEMD were COL2A1-related dysplasia and pseudoachondroplasia (OMIM# 177170).
At present, most of the reports of SEMD are single case reports, and there are few studies focusing on exploring the phenotype of SEMD caused by different disease-causing genes. We reported seven families with SEMD caused by TRPV4, COL2A1, CCN6, SBDS, and ACAN genes in order to explore the relationship between phenotype and genotype of them. Several studies have reported the phenotype–genotype relations of skeletal disorders caused by COL2A1, COMP, and CCN6 genes (Barat-Houari et al., 2016a; Madhuri et al., 2016; Liang et al., 2021), but we still need more families to prove it, especially Chinese families. Therefore, we summarized the clinical manifestations, radiological data, and molecular features of patients with SEMD, hoping to improve our understanding of Chinese families with SEMD.
Methods
Patients with clinical assessment
Seven non-consanguineous families (11 affected individuals) with features of spinal, epiphyseal, and metaphyseal dysplasia participated in this study from 2014–2020. Clinical data were collected, including basic information, family history, age of onset, physical examinations [included height, weight, arm, span, upper/lower segment (U/L) ratio, and standard deviation score (SDS) (Li, 2009)], and skeletal and other system manifestations. Radiographs of the spine, pelvises, and affected joints were taken. Laboratory examinations including liver and kidney functions, blood routine, and bone turnover markers were collected. The study was reviewed and approved by the Ethics Committee of the Shanghai Jiao Tong University Affiliated Sixth People’s Hospital and conducted in accordance with the Declaration of Helsinki. Informed written consents were obtained from all study participants or their legal guardians.
Targeted exome sequencing
Peripheral blood samples were collected in EDTA tubes from seven probands and 20 available family members. Informed consent was obtained from all 27 participants before the genetic analysis. The QuickGene DNA whole blood kits (Kurabo Industries Ltd., Osaka, Japan) and the Nucleic Acid Isolation System (QuickGene-610L; Autogen, Inc., Holliston, MA, United States) were used to extract genomic DNA. We used the TES to detect candidate genes of seven probands’ DNA samples. TES was performed using the Agilent SureSelect (SureSelect Reagent Kit; Agilent Technologies, CA, United States) target enrichment kit, which consisted of 322 pathogenic genes known to cause skeletal disease (Lv et al., 2021). Sequencing was performed on the Illumina Hiseq platform using paired-end 150-bp reads.
Evaluation of variants
We prioritized the variants that were known to cause disease in the Human Gene Mutation Database (HGMD) or other databases. The variants had lower allele frequency (< 0.001) in healthy control population databases (dbSNP, 1,000 Genomes, ESP6500, and gnomAD). The conservation and pathogenicity of variants were predicted by in silico tools (UniProt, Mutation Taster, Polyphen2, SIFT, and Human Splicing Finder). The pathogenicity of variants was also determined through appraisal of patients’ clinical phenotypes and scientific literature studies. The interpretation criteria of variants consisted of the American College of Medical Genetics and Genomics (ACMG) practice guidelines (Richards et al., 2015). Sanger sequencing was used to validate the candidate variants in all family members. Novel variants were recognized by using the Human Gene Mutation Database (HGMD) and ClinVar (Landrum et al., 2018; Stenson et al., 2020).
Results
This study included 11 patients derived from seven unrelated families who suffered from spinal, epiphyseal, and metaphyseal abnormalities. They are all of Han nationality from eastern Asia. Also, the ratio of females to males was 2/5, with 2 females and 5 males. The mean age of onset occurred at 4.9 ± 4.7 years. Detailed clinical characteristics of seven probands are shown in Table 1. Two families had an obvious positive family history, and the specific pedigrees of these seven families are shown in Figure 2.
FIGURE 2. Pedigrees of seven families with SEMD. Almost all accessible family members were sequenced and examined. Mutations are shown below each subject, and the corresponding sequence diagrams are displayed next to the pedigrees. Circles and squares indicate females and males, respectively. Arrows identify the proband in the families. Slashes indicate deceased individuals.
TRPV4-related skeletal dysplasia (OMIM# 605427)
In family 1, proband 1 (Ⅱ-1) was an 8.6-year-old girl with normal pregnancy, delivery, and family history. Her height was 125 cm (−1.1 SDS), and her arm span was 118 cm. She came to medical attention at the age of 8 when her mother noticed her scoliosis. Physical examinations revealed mild scoliosis, increased lumbar lordosis, and normal gait. She had myopia and normal hearing and intelligence. X-rays demonstrated mild scoliosis, ovoid vertebral bodies, rough edges of the vertebral bodies, and disorganization of the metaphysis (Figure 3).
FIGURE 3. Photographs and radiographs of families 1 and 2 with TRPV4-related skeletal dysplasia. (A) Radiographs of the spine and pelvis of proband 1 (Ⅱ-1) in family 1. Mild scoliosis, ovoid vertebral bodies, and rough edges of the vertebral bodies were significant in the spine. Disorganization of metaphysis was in the pelvis. (B,C) Radiographs of the spine and pelvis (B), and photographs of feet (C) of proband 2 (Ⅲ-1) in family 2. (B) Rough edges of the vertebral bodies and lumbar lordosis were visible in the spine. Disorganization of the metaphysis and shortening of the femoral neck were in the pelvis. (C) Fourth metatarsal dysplasia was in the left foot. (D) Radiographs of hands, spine, pelvis, and feet of proband 2’s father (Ⅱ-1). Mild contracture of the distal joints was shown in the hands. Scoliosis and depression of the anterior edge of vertebral bodies were displayed in the spine. Flattening of the acetabular roof, shortening of the femoral neck, and stenosis of joint spaces were in the pelvis. Hallux valgus deformity was seen in feet. (E) Radiographs of the spine and pelvis of proband 2’s brother (Ⅲ-2). Irregular shape of the vertebral bodies in the spine. Radiographs of the pelvis were normal.
Family 2 included three affected individuals, but their clinical phenotypes were highly variable. The 13-year-old boy proband 2 (Ⅲ-1) was 155.4 cm (−0.5 SDS) tall. Since the age of 10, there was a significant kyphosis in his lumbar without low back pain and joint movement limitation. Physical examinations revealed severe lumbar lordosis and fourth metatarsal dysplasia. His father (Ⅱ-1, 35 years old) had mild contracture of the distal joints of both hands and hallux valgus deformity of both feet since the age of 10 and difficulty in walking now. His brother (Ⅲ-1, 7 years old) did not show any abnormalities. Their heights were within the normal range; his father was 172 cm tall, and his brother was 123.5 cm tall. X-rays of the spine and pelvis in proband 2 showed rough edges of the vertebral bodies, lumbar lordosis, disorganization of the metaphysis, and shortening of the femoral neck. X-rays of the spine, pelvis, and feet of his father manifested scoliosis, depression of the anterior edge of vertebral bodies, flattening of the acetabular roof, shortening of the femoral neck, stenosis of joint spaces, and hallux valgus deformity (Figure 3). However, his brother’s X-ray showed a mild irregularity in the shape of the vertebral bodies.
COL2A1-related skeletal dysplasia Spondyloepiphyseal dysplasia congenita (OMIM# 183900)
In family 3, proband 3 (Ⅱ-1) was a 9-year-old boy, with a height of 110 cm (−5.6 SDS). At 33 weeks of gestation, his mother’s color doppler ultrasound showed that the fetal femur was shortened by 1 cm. He had visible thoracic abnormalities (pectus excavatum) at birth, learned to walk at the age of 2, and underwent thoracic correction surgery at the age of 3. Since he was 6 years old, he had hip stiffness with pain and limited mobility. At present, he has been unable to straighten his waist. Physical examinations showed that both maxillary lateral incisors were absent. X-rays showed scoliosis, platyspondyly, flattening and irregularity of the acetabular roof, femoral head necrosis, and disappearance of the femoral neck (Figure 4).
FIGURE 4. Photographs and radiographs of the families 3 and 4 with COL2A1-related skeletal dysplasia. (A,B) Radiographs of the spine and pelvis (A) and photographs of overall body and face (B) of proband 3 (Ⅱ-1) in family 3. (A) Scoliosis and platyspondyly were shown in the spine. Flattening and irregularity of the acetabular roof, femoral head necrosis, and disappearance of the femoral neck were visible in the pelvis. (B) Barrel chested and limited hip movement were performed in the overall body. No eruption of bilateral incisors was displayed in the face. (C,D) Radiographs of the spine, pelvis, and upper limb (C) and photographs of hands, elbows, knees, and feet (D) of proband 4 (Ⅱ-1) in family 4. (C) Severe multiple compressions of the vertebral bodies was revealed in the spine. Flattening of the acetabulum, shortening of the femoral neck, enlargement of the trochanter and femoral head, and narrow joint space of bilateral hips were shown in the pelvis. Widened epiphysis and metaphysis and joint cavity stenosis were visible in the elbows. (D) Bilateral elbow, knee, ankle, multiple metacarpophalangeal and interphalangeal joints were enlarged.
Kniest dysplasia (OMIM# 156550)
In family 4, the 29-year-old male proband 4 (Ⅱ-1) was 146 cm (−4.4 SDS) tall. He was born with enlarged knee joints, and thereafter, bilateral interphalangeal, elbow, knee, and ankle joints were enlarged gradually. At the age of 5, knee joint pain began and gradually aggravated, resulting in difficulty in walking. Proband 4 had poor visual acuity (retinal detachment in the left eye at the age of 24, cataract in the right eye at the age of 25, and both underwent surgery) and normal hearing. X-rays displayed severe multiple compressions of the vertebral bodies, flattening of the acetabulum, shortening of the femoral neck, narrow joint space of bilateral hips, and enlargement of the trochanter and femoral head (Figure 4).
Progressive pseudorheumatoid dysplasia (OMIM# 208230)
In family 5, the female proband 5 (Ⅱ-1) was 15 years old, with a height of 154.9 cm (−0.8 SDS). Intermittent large joint pain occurred at the age of 10, and no obvious abnormality was found in physical examinations. At the age of 13, the extension of both elbows was limited, and the elbow, knee, ankle, metacarpophalangeal, and interphalangeal joints gradually enlarged. Physical examinations showed that extension of multiple joints was limited, and no extraosseous manifestations were found. Symptoms aggravated with age. X-rays revealed mild platyspondyly, expansion of the epiphysis, and joint cavity stenosis. The pelvis was normal (Figure 5).
FIGURE 5. Photographs and radiographs of the families 5 and 6. (A,B) Radiographs of the spine, pelvis, hands, knees (A) and photographs (B) of proband 5 (Ⅱ-1) in family 5. (A) Mild platyspondyly was revealed in the spine. The pelvis was normal. Expansion of metaphysis and stenosis of joint spaces were visible in hands. Widened epiphysis and joint cavity stenosis were displayed in the knees. (B) Bilateral elbow, knee, ankle, and multiple metacarpophalangeal and interphalangeal joints were enlarged. (C) Radiographs of the spine, pelvis, and knees of proband (Ⅱ-1) in family 6. Irregular shape of the vertebral bodies was significant in the spine. Flattening of the acetabular roof, shortening of the femoral neck, mildly compressed femoral head, and thickening of metaphyseal texture were displayed in the pelvis. Widened epiphysis and metaphysis, irregular shape of the epiphysis, and disorganization of the metaphysis were obvious in the knees.
Shwachman–Bodian–Diamond syndrome (OMIM# 260400)
In family 6, the male proband 6 (Ⅱ-1) was 12 years old and had a height of 116 cm (−4.8 SDS). He was born prematurely, and growth retardation occurred at birth. He developed knee pain at the age of 8, which was aggravated after exercise. Physical examinations revealed mild intellectual disability, tooth dysplasia, and normal gait. X-rays demonstrated irregular shape of the vertebral bodies and epiphysis, expansion of metaphysis, and disorganization of metaphysis (Figure 5).
Aggrecan-related spondyloepiphyseal dysplasia (OMIM# 608361)
In family 7, the 48-year-old male proband 7 (Ⅱ-2) had a height of 148 cm (−4.0SDS) and an arm span of 135 cm, with a positive family history. Brachydactyly, clubfoot, and growth retardation were found at birth. When learning to walk at the age of 1, he had mild knee valgus and a swinging gait. Pain in the spine, hip, and knee joints appeared at the age of 16, and symptoms worsened after 20 years. Physical examinations revealed stubby fingers and toes, mild knee valgus, and limited hip mobility. X-rays displayed mild flattening and irregular shape of the vertebral bodies, 4–5 lumbar spondylolisthesis, flattening of the acetabulum, hip subluxation, femoral head necrosis, shortening of the femoral neck, osteoarthritis, widened epiphysis, and joint cavity stenosis (Figure 6). Both his mother (Ⅰ-2, 73 years old, 150 cm, −2.0 SDS) and daughter (Ⅲ-1, 22 years old, 148 cm, −2.3 SDS) had similar symptoms: short stature with short limbs. Her daughter’s X-ray reports showed similar findings, but her pelvis was normal.
FIGURE 6. Photographs and radiographs of the family 7. (A,B) Radiographs of the spine, pelvis, hands, and feet (A) and photographs (B) of proband 7 (Ⅱ-2) in family 7. (A) Mild flattening and irregular shape of the vertebral bodies and 4–5 lumbar spondylolisthesis were shown in the spine. Flattening of the acetabulum, hip subluxation, femoral head necrosis, and shortening of the femoral neck were revealed in the pelvis. Widened epiphysis and joint cavity stenosis were displayed in the hands and feet. (B) Short stature with short limb was visible in the full body. Thick and short fingers and toes were shown in hands and feet. (C,D) Radiographs of the spine, pelvis, hands, and feet (C) and photographs (D) of patient (Ⅲ-1) in family 7. (C) Mild flattening and irregular shape of the vertebral bodies were displayed in the spine. The pelvis was normal. Widened epiphysis and joint cavity stenosis were displayed in the hands and feet. (D) Thick and short fingers and toes were shown in hands and feet.
Genetic analysis
We have identified TRPV4, COL2A1, CCN6, SBDS, and ACAN gene mutations in 8 p.R232C (c.694C> T) were detected in two probands and three family members. In family 3, a heterozygous deletion mutation in COL2A1, p.G1089del (c.3266_3268del), was detected in proband 3. In family 4, a heterozygous splicing mutation in COL2A1, c.654+1G > C, was detected in proband 4. In family 5, two missense mutations in CCN6, p.C132W (c.396 T>G) and p.C241R (c.721 T>C), were detected in proband 5. Furthermore, proband 5’s mother (c.721 T>C) and father (c.396 T>G) carried a mutation, respectively. In family 6, a homozygous splicing mutation in SBDS, c.258+2 T>C, was detected in proband 6. Also, we found a heterozygous splicing mutation (c.258+2 T>C) in proband 6’s mother. In family 7, a heterozygous missense mutation in ACAN, p.T503K (c.1508C> A), was detected in proband 7 and two affected family members. A total of seven variants were found, and the pathogenicity prediction and classification are shown in Table 2. Combined with literature studies, two variants (c.654 + 1G > C and p.T503K) were novel, and five variants (p.R232C, p.G1089del, p.C132W, p.C241R, and c.258+2 T>C) were reported previously.
Discussion
This study focused on exploring phenotype and genotype in seven families with SEMD caused by TRPV4, SBDS, COL2A1, CCN6, and ACAN gene mutations, including two TRPV4-related skeletal dysplasia, one spondyloepiphyseal dysplasia, one Kniest dysplasia, one progressive pseudorheumatoid dysplasia, one Shwachman–Bodian–Diamond syndrome, and one aggrecan-related SED. The main clinical feature was short stature, and four probands had severe short stature (<-3 SDS). In addition, four probands had extraosseous manifestations. We have reported seven different mutation sites, including two novel sites and five reported sites. These findings expanded the phenotypic and genetic spectrum of SEMD in the Chinese population.
TRPV4 gene encodes a nonselective calcium-permeable ion channel, widely distributed in bone, nerve, lung, heart, and other tissues (Everaerts et al., 2010). So far, more than 70 variants in the TRPV4 gene had been reported to be related to autosomal-dominant skeletal dysplasia and motor and sensory neuropathies, such as Charcot-Marie-Tooth type 2C (CMT2C, OMIM#606482), scapuloperoneal spinal muscular atrophy (SPSMA, OMIM# 181405), spondylometaphyseal dysplasia-Kozlowski type (SMDK, OMIM#184252), and spondylo-epi-metaphyseal dysplasia Maroteaux pseudo-Morquio type 2 (SEDM-PM2, OMIM# 184095). The R232C variant had been reported repeatedly. Two studies reported that patients with R232C had a CMT2C-SPSMA overlap syndrome (Klein et al., 2011; Koutsis et al., 2015). They all had vocal cord paralysis, which was not found in this study. The probands in families 1 and 2 with R232C manifested slight short stature and scoliosis. Moreover, X-rays showed dysplasia of the vertebral bodies. No evidence of distal purely motor neuropathy was found in two probands, but these symptoms progressed very slowly. Therefore, we should follow up closely. Although only a small number of patients have been identified, they had comorbidity of neuromuscular and skeletal dysplasia (Cho et al., 2012; Faye et al., 2019). Interestingly, the clinical manifestations of family members who found the same variants showed intrafamily and interfamily variabilities. We detected that proband 1′s father harbored the R232C variant unexpectedly. He does not seem to be affected. Proband 2′s father and brother harbored the same variant. The phenotype of his father was the most serious, while his brother had no obvious abnormality, and only a mild irregularity in the shape of the vertebral bodies was found on the X-ray. These performances reinforced the notion that clinical features are heterogeneous in TRPV4-related neuropathies and skeletal dysplasia (Echaniz-Laguna et al., 2014). Meanwhile, the R232C variant exhibited reduced penetrance, which is the same as other reports (Chen et al., 2010; Koutsis et al., 2015). The R236C and R315W variants had also been reported to have reduced penetrance (Auer-Grumbach et al., 2010; Jędrzejowska et al., 2019). However, this phenomenon has not been reported in families with TRPV4-associated skeletal dysplasia (Nishimura et al., 2012).
The COL2A1 gene encodes the alpha-1 chain of type II procollagen, which is the main component of the nucleus pulposus of intervertebral discs, the vitreous humor of the eyes, and hyaline cartilage extracellular matrix (Li et al., 1995). The phenotypic spectrum of COL2A1-related diseases is very broad and includes mainly spondyloepiphyseal dysplasia congenita (SEDC, OMIM#183900), Kniest dysplasia, spondyloperipheral dysplasia (SPPD, OMIM#271700), osteoarthritis with mild chondrodysplasia (OSCDP, OMIM# 604864), and spondyloepiphyseal dysplasia-Strudwick type (SED-ST, OMIM#184250). It encompasses a diverse group of clinical phenotypes characterized by short stature and ocular manifestations, whereas later it is manifested as isolated arthritis. SEDC is characterized by a short trunk, flattened vertebral bodies, and abnormal epiphyses. In our study, proband 3 had a short trunk (U/P = 0.91), severe hip joint mobility limitation, and femoral head necrosis. He had a deletion of glycine in Gly-X-Y repeats (p.Gly1089del), which led to severe impairment of protein assembly and stability (Barat-Houari et al., 2016b). Patients with C-propeptide glycine substitutions have been reported to be shorter than those with N-propeptide substitutions (Terhal et al., 2012). Our center has reported that glycine to serine substitution caused milder phenotypes than glycine to nonserine substitutions (Xu et al., 2020). In Kniest dysplasia, the classical phenotypes are short-trunk dwarfism due to severely affected skeletal growth, scoliosis, platyspondyly, and joint enlargement, while extraskeletal features mainly include myopia, prominent eyes, conductive hearing loss, and mid-face hypoplasia (Barat-Houari et al., 2016a). Proband 4 with Kniest dysplasia had severe short stature (−4.4 SDS, U/L = 0.92), severe skeletal deformity, retinal detachment and cataract, and normal face and hearing. It is worth mentioning that proband 4 was misdiagnosed as PPD before, but he was detected with a COL2A1 mutation. This reminds us that we should pay attention to the differential diagnosis of PPD and COL2A1-related SEDC. The onset age of SEDC is usually at birth, while PPD mostly occurs at the age of 3–8 years (Garcia Segarra et al., 2012). At least 33 mutations were reported to cause Kniest dysplasia, and we reported a novel splice-site mutation (c.654 + 1G > C) which led to exon skipping.
Progressive pseudorheumatoid dysplasia (PPD) is a rare autosomal recessive disease caused by the functional loss or abnormality of cellular communication network factor 6 (CCN6). The clinical features are progressive joint stiffness and enlargement without inflammation (Dalal et al., 2012). Patients with PPD are usually normal at birth, with an onset ranging from 1 to 16 years of age. (Delague et al., 2005; Dalal et al., 2012). In our study, proband 5 developed symptoms at an average age of 10 years and were misdiagnosed as rheumatoid arthritis, but there was no inflammation. PPD initially presented with swelling and stiffness of the interphalangeal joint, gradually involving all large joints. According to literature statistics, the height of adults is usually less than the 3rd percentile (Garcia Segarra et al., 2012). In the early stages of the disease, the stature and proportion of patients may still be normal. Although patients had enlarged joints and limited mobility, X-rays showed no evidence of osteoarthritis. As cartilage destruction increases, secondary osteoarthritis may develop in adulthood. In our study, we found two reportedly missense mutations (C132W and C241R). The C241R variant has been reported many times in Chinese families and seems to be a hot site in Chinese (Ye et al., 2012; Yu et al., 2015; Hu et al., 2017). One research has reported that most PPD patients carry nonsense mutations in both WISP3 alleles. They observed that the phenotype of patients homozygous for a missense mutation was not milder than patients with one or both alleles carrying a nonsense mutation (Garcia Segarra et al., 2012).
Shwachman–Bodian–Diamond syndrome (SBDS, OMIM#260400) is an autosomal recessive disease characterized by hematologic abnormalities, exocrine pancreatic dysfunction, immune deficiency, and skeletal abnormalities (Ginzberg et al., 1999). More than 90% of SBDS were caused by SBDS gene mutations (Boocock et al., 2003), and the most common variants were c.184A > T and c.258 + 2 T > C (Kuijpers et al., 2005; Nelson and Myers, 2018). The diagnosis of SBDS relied on evidence of exocrine pancreatic dysfunction, hematologic abnormalities, and recurrent infections (Huang and Shimamura, 2011; Nelson and Myers, 2018). However, some patients did not have this classic combination of manifestations (Myers et al., 2014). One study evaluated 102 genetically diagnosed patients with SBDS and found 40% had hematologic complications (Donadieu et al., 2012). According to research in China, the onset age of hemocytopenia ranges from 0 to 12 years old (An et al., 2020). In our study, proband 6 had normal hematological function. Therefore, it is necessary to conduct long-term follow-up on the hematologic complications of proband 6. Skeletal dysplasia is also a common manifestation, Mäkitie reported that the skeletal abnormalities in patients with SBDS varied not only with age but also with individuals (Mäkitie et al., 2004). Proband 6 had obvious skeletal manifestations, severe short stature (< −3.0 SDS), and epiphyseal and metaphyseal dysplasia. One study found that patients with SBDS ranged widely in their neurocognitive impairment compared with controls (Kerr et al., 2010) and had been shown to be related to structural brain alterations (Toiviainen-Salo et al., 2008). It shows that the phenotypic spectrum of SBDS is wider than previously thought. One study summarized the literature for large SBDS cohorts and found that no obvious genotype–phenotype correlation was observed (Mäkitie et al., 2004; Donadieu et al., 2012; Myers et al., 2020).
The cartilage aggrecan proteoglycan is encoded by the ACAN gene, the main structural component of the cartilage growth plate, which is essential for skeletal growth and articular cartilage function (Roughley and Mort, 2014). The mutations in the ACAN gene can be classified into four categories: spondyloepiphyseal dysplasia, Kimberley type (SEDK, OMIM#608361); spondylo-epi-metaphyseal dysplasia, aggrecan type (SEMD, OMIM#612813); osteochondritis dissecans (OCD, OMIM#165800); and short stature and advanced bone age. Previously, ACAN mutations had been reported as a cause of short stature with a frequency of 1.4–37.5% in the short stature population (Hauer et al., 2017; Stavber et al., 2020). In 2021, Li reported a large cohort of Chinese short-stature children caused by ACAN mutations (Lin et al., 2021). Patients with SEDK mostly have brachydactyly, platyspondyly, irregular femoral epiphyses, and precocious arthropathy, similar to the patients in this study. However, other features such as facial dysmorphisms and discopathy were not found. In contrast to previous data, brachydactyly was observed in all probands (Sentchordi-Montané et al., 2018). Nevertheless, there are few reports of SEDK caused by ACAN mutations (Gleghorn et al., 2005; Sentchordi-Montané et al., 2018). The new variant we reported expanded the clinical phenotype and genotype of this disease. At present, many clinical trials have proved that combined treatment with growth hormone (GH) or gonadotropin-releasing hormone (GnRH) is used to improve the final height of patients (Lin et al., 2021; Wei et al., 2021). Gkourogianni reported that the average height of patients treated with GH was lower than normal at 2.5 SD, while that of untreated patients was lower than 3 SD (Gkourogianni et al., 2017). The efficacy of GH and GnRH analogs in the treatment of patients with ACAN mutations needs more studies to be proved.
This study had some limitations. First, we were unable to perform electromyogram (EMG) testing in patients from family 1 and family 2 because the patients refused. We failed to determine whether the patients had manifestations of distal purely motor neuropathy. Second, as this was a single-center study, the number of probands with SEMD was small. Therefore, we did not have sufficient evidence to elucidate the phenotype–genotype relationship in SEMD. In the future, we hope to accumulate more patient data to illustrate the correlation between the two.
In the present study, we summarized seven families with SEMD caused by TRPV4, COL2A1, CCN6, SBDS, and ACAN gene mutations. Our results further demonstrated that phenotypic and genotypic spectra of SEMD are very wide. We should pay attention to the complications of SEMD and reduce or prevent their effects. Genetic testing is important for diagnosis, genetic counseling, and prenatal diagnosis. The new findings reported herein would contribute to the investigation of the phenotype–genotype relations among patients with SEMD.
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 Shanghai Jiao Tong University Affiliated Sixth People’s Hospital. Written informed consent to participate in this study was provided by the participants’ 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: ZZ and SL; data curation: SL, JZ, and LL; formal analysis: ZZ and SL; funding acquisition: ZZ, SL, HY, and CW; methodology: ZZ; project administration: ZZ and SL; resources: HZ and HY; supervision: ZZ, SL, and CW; writing—original draft: SL; writing—review and editing: ZZ and SL.
Funding
The study was supported by the National Key Research and Development Program of China (2018YFA0800801), the National Natural Science Foundation of China (NSFC) (81974123 and 81900807), and the Clinical Science and Technology Innovation Project of Shanghai Shenkang Hospital Development Center (SHDC12018120).
Acknowledgments
We are grateful to Genesky Biotechnologies Inc., (Shanghai) for its technical support and to the patients and their families for their participation 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.
References
Alanay, Y., and Lachman, R. S. (2011). A review of the principles of radiological assessment of skeletal dysplasias. Jcrpe 3 (4), 163–178. doi:10.4274/jcrpe.463
An, W. B., Liu, C., Wan, Y., Chang, L. X., Chen, X. Y., and Zhu, X. F. (2020). Clinical features and gene mutations of children with Shwachman-Diamond syndrome and malignant myeloid transformation. Zhongguo Dang Dai Er Ke Za Zhi 22 (5), 460–465. doi:10.7499/j.issn.1008-8830.2001133
Auer-Grumbach, M., Olschewski, A., Papić, L., Kremer, H., McEntagart, M. E., Uhrig, S., et al. (2010). Alterations in the ankyrin domain of TRPV4 cause congenital distal SMA, scapuloperoneal SMA and HMSN2C. Nat. Genet. 42 (2), 160–164. doi:10.1038/ng.508
Barat-Houari, M., Dumont, B., Fabre, A., Them, F. T., Alembik, Y., Alessandri, J. L., et al. (2016a). The expanding spectrum of COL2A1 gene variants IN 136 patients with a skeletal dysplasia phenotype. Eur. J. Hum. Genet. 24 (7), 992–1000. doi:10.1038/ejhg.2015.250
Barat-Houari, M., Sarrabay, G., Gatinois, V., Fabre, A., Dumont, B., Genevieve, D., et al. (2016b). Mutation update forCOL2A1Gene variants associated with type II collagenopathies. Hum. Mutat. 37 (1), 7–15. doi:10.1002/humu.22915
Boocock, G. R., Morrison, J. A., Popovic, M., Richards, N., Ellis, L., Durie, P. R., et al. (2003). Mutations in SBDS are associated with Shwachman-Diamond syndrome. Nat. Genet. 33 (1), 97–101. doi:10.1038/ng1062
Briggs, M. D., Hoffman, S. M., King, L. M., Olsen, A. S., Mohrenweiser, H., Leroy, J. G., et al. (1995). Pseudoachondroplasia and multiple epiphyseal dysplasia due to mutations in the cartilage oligomeric matrix protein gene. Nat. Genet. 10 (3), 330–336. doi:10.1038/ng0795-330
Chen, D. H., Sul, Y., Weiss, M., Hillel, A., Lipe, H., Wolff, J., et al. (2010). CMT2C with vocal cord paresis associated with short stature and mutations in the TRPV4 gene. Neurology 75 (22), 1968–1975. doi:10.1212/WNL.0b013e3181ffe4bb
Cho, T. J., Matsumoto, K., Fano, V., Dai, J., Kim, O. H., Chae, J. H., et al. (2012). TRPV4-pathy manifesting both skeletal dysplasia and peripheral neuropathy: A report of three patients. Am. J. Med. Genet. 158a (4), 795–802. doi:10.1002/ajmg.a.35268
Cormier-Daire, V. (2008). Spondylo-epi-metaphyseal dysplasia. Best Pract. Res. Clin. Rheumatology 22 (1), 33–44. doi:10.1016/j.berh.2007.12.009
Dalal, A., Bhavani G, G. S., Togarrati, P. P., Bierhals, T., Nandineni, M. R., Danda, S., et al. (2012). Analysis of theWISP3gene in Indian families with progressive pseudorheumatoid dysplasia. Am. J. Med. Genet. 158a (11), 2820–2828. doi:10.1002/ajmg.a.35620
Delague, V., Chouery, E., Corbani, S., Ghanem, I., Aamar, S., Fischer, J., et al. (2005). Molecular study ofWISP3 in nine families originating from the Middle-East and presenting with progressive pseudorheumatoid dysplasia: Identification of two novel mutations, and description of a founder effect. Am. J. Med. Genet. 138a (2), 118–126. doi:10.1002/ajmg.a.30906
Donadieu, J., Fenneteau, O., Beaupain, B., Beaufils, S., Bellanger, F., Mahlaoui, N., et al. (2012). Classification of and risk factors for hematologic complications in a French national cohort of 102 patients with Shwachman-Diamond syndrome. Haematologica 97 (9), 1312–1319. doi:10.3324/haematol.2011.057489
Echaniz-Laguna, A., Dubourg, O., Carlier, P., Carlier, R. Y., Sabouraud, P., Pereon, Y., et al. (2014). Phenotypic spectrum and incidence of TRPV4 mutations in patients with inherited axonal neuropathy. Neurology 82 (21), 1919–1926. doi:10.1212/wnl.0000000000000450
Everaerts, W., Nilius, B., and Owsianik, G. (2010). The vanilloid transient receptor potential channel TRPV4: From structure to disease. Prog. Biophysics Mol. Biol. 103 (1), 2–17. doi:10.1016/j.pbiomolbio.2009.10.002
Faye, E., Modaff, P., Pauli, R., and Legare, J. (2019). Combined phenotypes of spondylometaphyseal dysplasia-kozlowski type and charcot-marie-tooth disease type 2C secondary to a TRPV4 pathogenic variant. Mol. Syndromol. 10 (3), 154–160. doi:10.1159/000495778
Garcia Segarra, N., Mittaz, L., Campos-Xavier, A. B., Bartels, C. F., Tuysuz, B., Alanay, Y., et al. (2012). The diagnostic challenge of progressive pseudorheumatoid dysplasia (PPRD): A review of clinical features, radiographic features, and WISP3 mutations in 63 affected individuals. Am. J. Med. Genet. 160c (3), 217–229. doi:10.1002/ajmg.c.31333
Ginzberg, H., Shin, J., Ellis, L., Morrison, J., Ip, W., Dror, Y., et al. (1999). Shwachman syndrome: Phenotypic manifestations of sibling sets and isolated cases in a large patient cohort are similar. J. Pediatr. 135 (1), 81–88. doi:10.1016/s0022-3476(99)70332-x
Gkourogianni, A., Andrew, M., Tyzinski, L., Crocker, M., Douglas, J., Dunbar, N., et al. (2017). Clinical characterization of patients with autosomal dominant short stature due to aggrecan mutations. J. Clin. Endocrinol. Metab. 102 (2), 460–469. doi:10.1210/jc.2016-3313
Gleghorn, L., Ramesar, R., Beighton, P., and Wallis, G. (2005). A mutation in the variable repeat region of the aggrecan gene (AGC1) causes a form of spondyloepiphyseal dysplasia associated with severe, premature osteoarthritis. Am. J. Hum. Genet. 77 (3), 484–490. doi:10.1086/444401
Hauer, N. N., Sticht, H., Boppudi, S., Büttner, C., Kraus, C., Trautmann, U., et al. (2017). Genetic screening confirms heterozygous mutations in ACAN as a major cause of idiopathic short stature. Sci. Rep. 7 (1), 12225. doi:10.1038/s41598-017-12465-6
Hu, Q., Liu, J., Wang, Y., Wang, J., Shi, H., Sun, Y., et al. (2017). Delayed-onset of progressive pseudorheumatoid dysplasia in a Chinese adult with a novel compound WISP3 mutation: A case report. BMC Med. Genet. 18 (1), 149. doi:10.1186/s12881-017-0507-3
Huang, J. N., and Shimamura, A. (2011). Clinical spectrum and molecular pathophysiology of Shwachman-Diamond syndrome. Curr. Opin. Hematol. 18 (1), 30–35. doi:10.1097/MOH.0b013e32834114a5
Jędrzejowska, M., Dębek, E., Kowalczyk, B., Halat, P., Kostera-Pruszczyk, A., Ciara, E., et al. (2019). The remarkable phenotypic variability of the p.Arg269HiS variant in the TRPV4 gene. Muscle Nerve 59 (1), 129–133. doi:10.1002/mus.26346
Kerr, E. N., Ellis, L., Dupuis, A., Rommens, J. M., and Durie, P. R. (2010). The behavioral phenotype of school-age children with shwachman diamond syndrome indicates neurocognitive dysfunction with loss of Shwachman-Bodian-Diamond syndrome gene function. J. Pediatr. 156 (3), 433–438. doi:10.1016/j.jpeds.2009.09.026
Klein, C. J., Shi, Y., Fecto, F., Donaghy, M., Nicholson, G., McEntagart, M. E., et al. (2011). TRPV4 mutations and cytotoxic hypercalcemia in axonal Charcot-Marie-Tooth neuropathies. Neurology 76 (10), 887–894. doi:10.1212/WNL.0b013e31820f2de3
Koutsis, G., Lynch, D., Manone, A., Karadima, G., Reilly, M. M., Houlden, H., et al. (2015). Charcot-Marie-Tooth disease type 2C and scapuloperoneal muscular atrophy overlap syndrome in a patient with the R232C TRPV4 mutation. J. Neurol. 262 (8), 1972–1975. doi:10.1007/s00415-015-7800-x
Kuijpers, T. W., Alders, M., Tool, A. T., Mellink, C., Roos, D., and Hennekam, R. C. (2005). Hematologic abnormalities in shwachman diamond syndrome: Lack of genotype-phenotype relationship. Blood 106 (1), 356–361. doi:10.1182/blood-2004-11-4371
Landrum, M. J., Lee, J. M., Benson, M., Brown, G. R., Chao, C., Chitipiralla, S., et al. (2018). ClinVar: Improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 46 (D1), D1062–d1067. doi:10.1093/nar/gkx1153
Li, H. (2009). Growth standardized values and curves based on weight, length/height and head circumference for Chinese children under 7 years of age. Zhonghua Er Ke Za Zhi 47 (3), 173–178.
Li, S. W., Prockop, D. J., Helminen, H., Fässler, R., Lapveteläinen, T., Kiraly, K., et al. (1995). Transgenic mice with targeted inactivation of the Col2 alpha 1 gene for collagen II develop a skeleton with membranous and periosteal bone but no endochondral bone. Genes Dev. 9 (22), 2821–2830. doi:10.1101/gad.9.22.2821
Liang, H., Hou, Y., Pang, Q., Jiang, Y., Wang, O., Li, M., et al. (2021). Clinical, biochemical, radiological, genetic and therapeutic analysis of patients with COMP gene variants. Calcif. Tissue Int. 110, 313–323. doi:10.1007/s00223-021-00920-6
Lin, L., Li, M., Luo, J., Li, P., Zhou, S., Yang, Y., et al. (2021). A high proportion of novel ACAN mutations and their prevalence in a large cohort of Chinese short stature children. J. Clin. Endocrinol. Metab. 106 (7), e2711–e2719. doi:10.1210/clinem/dgab088
Lv, S., Zhao, J., Xi, L., Lin, X., Wang, C., Yue, H., et al. (2021). Genetics evaluation of targeted exome sequencing in 223 Chinese probands with genetic skeletal dysplasias. Front. Cell Dev. Biol. 9, 715042. doi:10.3389/fcell.2021.715042
Madhuri, V., Santhanam, M., Rajagopal, K., Sugumar, L. K., and Balaji, V. (2016). WISP3 mutational analysis in Indian patients diagnosed with progressive pseudorheumatoid dysplasia and report of a novel mutation at p.Y198*. Bone & Jt. Res. 5 (7), 301–306. doi:10.1302/2046-3758.57.2000520
Mäkitie, O., Ellis, L., Durie, P. R., Morrison, J. A., Sochett, E. B., Rommens, J. M., et al. (2004). Skeletal phenotype in patients with Shwachman-Diamond syndrome and mutations in SBDS. Clin. Genet. 65 (2), 101–112. doi:10.1111/j.0009-9163.2004.00198.x
Miyoshi, K., Nakamura, K., Haga, N., and Mikami, Y. (2004). Surgical treatment for atlantoaxial subluxation with myelopathy in spondyloepiphyseal dysplasia congenita. Spine 29 (21), E488–E491. doi:10.1097/01.brs.0000143621.37688.f3
Mortier, G. R., Cohn, D. H., Cormier‐Daire, V., Hall, C., Krakow, D., Mundlos, S., et al. (2019). Nosology and classification of genetic skeletal disorders: 2019 revision. Am. J. Med. Genet. 179 (12), 2393–2419. doi:10.1002/ajmg.a.61366
Myers, K. C., Bolyard, A. A., Otto, B., Wong, T. E., Jones, A. T., Harris, R. E., et al. (2014). Variable clinical presentation of shwachman-diamond syndrome: Update from the north American shwachman-diamond syndrome registry. J. Pediatr. 164 (4), 866–870. doi:10.1016/j.jpeds.2013.11.039
Myers, K. C., Furutani, E., Weller, E., Siegele, B., Galvin, A., Arsenault, V., et al. (2020). Clinical features and outcomes of patients with shwachman-diamond syndrome and myelodysplastic syndrome or acute myeloid leukaemia: A multicentre, retrospective, cohort study. Lancet Haematol. 7 (3), e238–e246. doi:10.1016/s2352-3026(19)30206-6
Nelson, A. S., and Myers, K. C. (2018). Diagnosis, treatment, and molecular pathology of shwachman-diamond syndrome. Hematology/Oncology Clin. N. Am. 32 (4), 687–700. doi:10.1016/j.hoc.2018.04.006
Nishimura, G., Lausch, E., Savarirayan, R., Shiba, M., Spranger, J., Zabel, B., et al. (2012). TRPV4-associated skeletal dysplasias. Am. J. Med. Genet. 160c (3), 190–204. doi:10.1002/ajmg.c.31335
Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of medical genetics and genomics and the association for molecular pathology. Genet. Med. 17 (5), 405–424. doi:10.1038/gim.2015.30
Roughley, P. J., and Mort, J. S. (2014). The role of aggrecan in normal and osteoarthritic cartilage. J. Exp. Ortop. 1 (1), 8. doi:10.1186/s40634-014-0008-7
Sentchordi-Montané, L., Aza-Carmona, M., Benito-Sanz, S., Barreda- Bonis, A. C., Sánchez-Garre, C., Prieto-Matos, P., et al. (2018). Heterozygous aggrecan variants are associated with short stature and brachydactyly: Description of 16 probands and a review of the literature. Clin. Endocrinol. 88 (6), 820–829. doi:10.1111/cen.13581
Spranger, J. (1989). Radiologic nosology of bone dysplasias. Am. J. Med. Genet. 34 (1), 96–104. doi:10.1002/ajmg.1320340117
Stavber, L., Hovnik, T., Kotnik, P., Lovrečić, L., Kovač, J., Tesovnik, T., et al. (2020). High frequency of pathogenic ACAN variants including an intragenic deletion in selected individuals with short stature. Eur. J. Endocrinol. 182 (3), 243–253. doi:10.1530/eje-19-0771
Stenson, P. D., Mort, M., Ball, E. V., Chapman, M., Evans, K., Azevedo, L., et al. (2020). The human gene mutation database (HGMD): Optimizing its use in a clinical diagnostic or research setting. Hum. Genet. 139 (10), 1197–1207. doi:10.1007/s00439-020-02199-3
Terhal, P. A., van Dommelen, P., Le Merrer, M., Zankl, A., Simon, M. E., Smithson, S. F., et al. (2012). Mutation-based growth charts for SEDC and other COL2A1 related dysplasias. Am. J. Med. Genet. 160c (3), 205–216. doi:10.1002/ajmg.c.31332
Toiviainen-Salo, S., Mäkitie, O., Mannerkoski, M., Hämäläinen, J., Valanne, L., and Autti, T. (2008). Shwachman-Diamond syndrome is associated with structural brain alterations on MRI. Am. J. Med. Genet. 146a (12), 1558–1564. doi:10.1002/ajmg.a.32354
Wei, M., Ying, Y., Li, Z., Weng, Y., and Luo, X. (2021). Identification of novel ACAN mutations in two Chinese families and genotype-phenotype correlation in patients with 74 pathogenic ACAN variations. Mol. Genet. Genomic Med. 9 (11), e1823. doi:10.1002/mgg3.1823
Xu, Y., Li, L., Wang, C., Yue, H., Zhang, H., Gu, J., et al. (2020). Clinical and molecular characterization and discovery of novel genetic mutations of Chinese patients with col2a1-related dysplasia. Int. J. Biol. Sci. 16 (5), 859–868. doi:10.7150/ijbs.38811
Ye, J., Zhang, H. W., Qiu, W. J., Han, L. S., Zhang, Y. F., Gong, Z. W., et al. (2012). Patients with progressive pseudorheumatoid dysplasia: From clinical diagnosis to molecular studies. Mol. Med. Rep. 5 (1), 190–195. doi:10.3892/mmr.2011.619
Keywords: spondylo-epi-metaphyseal dysplasia, phenotype-genotype relation, TRPV4, COL2A1, CCN6, SBDS, ACAN
Citation: Lv S, Zhao J, Liu L, Wang C, Yue H, Zhang H, Li S and Zhang Z (2022) Exploring and expanding the phenotype and genotype diversity in seven Chinese families with spondylo-epi-metaphyseal dysplasia. Front. Genet. 13:960504. doi: 10.3389/fgene.2022.960504
Received: 03 June 2022; Accepted: 01 August 2022;
Published: 31 August 2022.
Edited by:
Rong Qiang, Northwest Women’s and Children’s Hospital, ChinaReviewed by:
Asmat Ullah, University of Copenhagen, DenmarkChao Xu, Shandong Provincial Hospital, China
Sajid Malik, Quaid-i-Azam University, Pakistan
Copyright © 2022 Lv, Zhao, Liu, Wang, Yue, Zhang, Li and Zhang. 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: Shanshan Li, lissane0508@126.com; Zhenlin Zhang, zhangzl@sjtu.edu.cn