Original Research ARTICLE
Genotypes and Phenotypes of DMD Small Mutations in Chinese Patients With Dystrophinopathies
- 1Department of Neurology, National Key Clinical Department and Key Discipline of Neurology, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 2Department of Dermatology, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 3Department of Dermatology, The Second Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Dystrophinopathies are a group of neuromuscular disorders resulting from mutations in DMD, including Duchenne muscular dystrophy (DMD), intermediate muscular dystrophy (IMD), and Becker muscular dystrophy (BMD). Herein, we present the characteristics of small mutations in Chinese patients with dystrophinopathies, and explore genotype–phenotype correlations. In our cohort, 115 patients with small mutations (18.49% of all patients) were included and DMD mutations were detected by either Sanger (53.91%) or next generation sequencing (46.09%). In total, 106 small mutations were detected, 28 of which (26.42%) had not been reported previously. The most common mutations were nonsense mutations (52.17%), followed by splicing (24.35%), frameshift (17.39%), and missense mutations (5.22%), in addition to a single untranslated region mutation (0.87%). We discovered distinct mutation characteristics in our patients, such as different positional distributions, indicating different exon skipping therapy strategies for small mutations in Chinese patients. Almost all patients (96.51%) with truncating or missense mutations, were covered by triple/double/single-exon skipping therapy; the most frequent single-exon skipping strategy was skipping exon 32, applicable for 8.51% of patients. Furthermore, splicing classification grades were correlated with phenotypes in nonsense mutations (P < 0.001), and serum creatinine levels differed significantly between DMD/IMD and BMD for patients ≤ 16 years old (P = 0.002). These observations can further aid prognostic judgment and guide treatment. In conclusion, the mutation characteristics and genotype–phenotype correlations in Chinese patients with dystrophinopathies and small mutations could provide insights into the molecular mechanisms of pathogenesis, diagnosis, and treatment designs.
Dystrophinopathies are X chromosome-linked recessive neuromuscular disorders with a pooled prevalence of 6.31 per 100,000 males (Mah et al., 2014), and are characterized by progressive muscle weakness, elevated levels of serum creatine kinase (CK), and death as a result of cardiorespiratory failure (Flanigan, 2014). These disorders can be classified as Duchenne muscular dystrophy (DMD; #310200), Becker muscular dystrophy (BMD; #300376), and intermediate muscular dystrophy (IMD). DMD has a heavy disease course with loss of ambulation before 13 years of age, whereas BMD exhibits mild progression and better prognosis (Bushby et al., 2010); the severity of IMD lies between that of DMD and BMD (Flanigan, 2014).
Dystrophinopathies result from mutations in DMD (OMIM: ∗300377). DMD encodes dystrophin, a structural protein with four domains: actin-binding domain (ABD), rod domain (RD), cysteine-rich domain (CRD), and C-terminal domain (CTD). Dystrophin is located in the myofiber membrane where it assists in membrane stability (Le Rumeur et al., 2010). The most common mutation type in dystrophinopathies involves large rearrangements, which account for 77.7% of mutations (Tuffery-Giraud et al., 2009). The remainder are DMD small mutations, mainly consisting of point mutations, small deletions, and small insertions (Tuffery-Giraud et al., 2009). Uncovering the mutation spectrum observed in specific diseases can lead to the development of genetic medicine, such as exon skipping therapy and nonsense read through therapy (Jarmin et al., 2014). The ‘reading-frame rule’ can explain approximately 91% of the genotype–phenotype correlations observed in patients with large rearrangements in DMD (Aartsma-Rus et al., 2006). However, the genotype–phenotype correlations of small mutations are more complicated (Aartsma-Rus et al., 2006); for example, nonsense mutations introduce premature stop codons (PTCs), which may result in truncated proteins and correlate with severe phenotypes (Ginjaar et al., 2000; Richards et al., 2015). However, an increasing number of recent studies have reported BMD cases with nonsense mutations (Ginjaar et al., 2000; Aartsma-Rus et al., 2006; Flanigan et al., 2011). Therefore, genotype–phenotype correlation with respect to DMD small mutations warrants further research.
Several studies have reported the mutation characteristics of small mutations but not focused on genotype–phenotype correlations in Chinese patients with dystrophinopathies (Chen et al., 2013; Yang et al., 2013; Guo et al., 2015). Therefore, we studied the characteristics of DMD small mutations and assessed genotype–phenotype correlations. Additionally, we examined the correlation between serum creatinine (SCRN) and phenotypes of patients with small mutations, as we had previously verified the correlation between SCRN and phenotypes of patients whose main mutations were large rearrangements (Zhang et al., 2015; Wang et al., 2017).
Materials and Methods
A total of 622 male patients with dystrophinopathies were diagnosed through clinical manifestations, biochemical changes, and genetic examinations between September 2012 and February 2018 at the First Affiliated Hospital, Sun Yat-sen University, China. Of these, 115 unrelated male patients with small mutations were included in this study. All of included patients were East Asian, with 104 of Han ethnicity, 1 Hui, and 10 of unknown ethnicities. Some subjects had also participated in our previous studies (Zhang et al., 2015; Wang et al., 2017).
Four phenotypic subgroups were defined as follows: DMD, IMD, BMD, and pending. These classifications were mainly based on clinical severity, such as age at which the patient lost ambulation (DMD < 13 years; BMD ≥ 16 years; and 13 years ≤ IMD < 16 years; Tuffery-Giraud et al., 2009). Some patients could not be classified in this manner, and thus those with an onset of weakness by the age of 5 were recognized as DMD, while those that retained nearly normal motor function, or had very mild motor dysfunction after the age of 5 were classified as BMD (Brooke et al., 1981; Marden et al., 2005).
Patients whose motor dysfunction severity fell between DMD and BMD were classified as IMD. This study included seven patients with IMD; two of whom lost ambulation at the age of 14 and another two (12 and 13 years old) had severe difficulty walking, manifested obvious waddling gaits, and lost the ability of standing up from seated positions. The remaining three patients were 11–12 years old and presented with obvious waddling gait as well as difficulties climbing stairs and, getting up off the ground or from a chair. They were classified as IMD in this study because their motor dysfunctions were milder than those with DMD and more severe than BMD patients at the same age. We could not completely exclude the possibility that the three IMD patients between 11–12 years old would lose ambulation before the age of 13, which would change their classifications to DMD. However, this does not affect our analysis because DMD and IMD were grouped as DMD/IMD during this study.
Patients whose phenotypes could not be determined or with missing disease history information (only one patient had missing clinical data) were classified as ‘pending.’
SCRN and Genetic Examination
Blood collection was performed as previously described (Wang et al., 2017). Sera were immediately tested for CK and SCRN using enzymatic reactions in a Beckman Coulter AU5800 clinical chemistry analyser (Beckman Coulter, Brea, CA, United States). Generally, results of CK and SCRN levels obtained from the first test were included for analysis of patients with several follow-up serum tests. For DMD mutation analysis, blood (2 mL) was collected and then multiplex ligation-dependent probe amplification (MLPA) was used to detect large sequence rearrangements. Blood samples with negative MLPA results were then checked for small-scale mutations by next-generation sequencing (NGS) using an Illumina HiSeq 2000 platform (Illumina, San Diego, CA, United States) or Sanger sequencing as previously described (Chen et al., 2013; Wang et al., 2017). Sanger sequencing was used to detect whether the probands’ mothers carried the DMD mutation. The DMD transcript number is NM_004006.2.
Twenty-five patients underwent muscle biopsy and then the muscle samples were used for immunohistochemistry. The quadriceps femoris and gastrocnemius muscles were extracted for analysis as previously (Yu et al., 2017). The antibody against dystrophin was purchased from Novocastra (Newcastle, United Kingdom). We analyzed expression levels of the dystrophin C-terminus as using antibodies against this C-terminus is a highly reliable diagnostic approach; furthermore, the dystrophin C-terminus is translated in patients with nonsense mutations if beneficial spontaneous exon skipping occurred (Doriguzzi et al., 1997; Flanigan et al., 2011).
The Leiden Open Variation Database (LOVD), Universal Mutation Database (UMD)-DMD Database, ClinVar database, and Google were used to detect whether mutations had been previously reported. The pathogenicity of a mutation was evaluated according to the guidelines of the American College of Medical Genetics and Genomics (Richards et al., 2015); the pathogenic or likely pathogenic mutations were included in this study. Domains where mutations were located were defined as ABD (exons 1–8), RD (exons 9–62), CRD (exons 63–69), and CTD (exons 70–79; Engel, 2004). The Leiden DMD reading frame checker1 was used to determine whether skipping the exon in which the mutation was located would lead to an open reading frame (ORF) shift [ORF shift (out-of-frame) was set as 0, otherwise it was 1 (in-frame)] (Aartsma-Rus et al., 2006). Mechanisms explaining spontaneous exon skipping of nonsense mutations varied; however, a splicing change resulting from the disruption or creation of a splicing regulatory element (SRE) appeared to be the primary mechanism that mainly occurred via exonic splicing silencer (ESS) and exonic splicing enhancer (ESE) elements (Flanigan et al., 2011). Human Splicing Finder (HSF2) was used to determine changes in ESS and ESE elements in mutated DMD (Desmet et al., 2009). The HSF score was 0 when there were no SRE changes and 1 when a new ESS or broken ESE was present. Splicing classification grades were obtained by multiplying ORF changes by the HSF score as beneficial exon skipping is based on in-frame changes after exon skipping. Grade 0 represents an absence of exon skipping or unbeneficial exon skipping that leads to changes in the reading frame, both of which produce PTCs; grade 1 represents beneficial exon skipping.
Data were analyzed using SPSS 24.0 (IBM, Armonk, NY, United States) and GraphPad Prism 6 (GraphPad Software, San Diego, CA, United States) software. Microsoft Excel 2016 (Microsoft, Redmond, Washington, DC, United States) was also used for plotting information. Contingency tables and Chi square tests were used for categorical variables. When any expected value of contingency tables was less than 5, the results were corrected for continuity, when the value was less than 1 or the total sample number was less than 40, the Fisher exact test was used. For statistical significance, odds ratios (ORs) were reported with 95% confidence intervals (CIs) in 2 × 2 contingency tables. Correlation between phenotypes and factors (age of examination, mutation type and SCRN, or CK levels) were analyzed using logistic regression. P < 0.05 was considered statistically significant.
Phenotypes and Mutation Characteristics
In our cohort of patients with dystrophinopathies (622), 115 (18.49%) presented with small mutations. Among them, DMD was observed in 59 subjects (51.30%), BMD in 31 (26.96%), and IMD in 7 (6.09%), while 18 were classified as pending (15.65%). We categorized DMD and IMD patients into one DMD/IMD group due to prognostic similarities and the low number of IMD.
Out of 115 patients, 106 mutations were detected, 28 of which (26.42%) had not been reported previously (search date: 2018/5/2). Small mutations in 62 (53.91%) and 53 (46.09%) subjects were detected by Sanger sequencing and NGS, respectively. There were seven common mutations; furthermore, the phenotypes were not noticeably different between patients with the same mutations (Table 1 and see Supplementary Table S1 for detailed patient information). As our study patients were unrelated, we analyzed genotype and phenotype characteristics based on the number of patients. The most common types of mutations were nonsense mutations (60/115, 52.17%), followed by splicing (28/115, 24.35%), frameshift (20/115, 17.39%), and missense mutations (6/115, 5.22%). A single untranslated region (UTR) mutation (1/115, 0.87%) was also detected. Of all patients with truncating mutations, those with nonsense mutations (60/80, 75.00%) were more common than those with frameshift mutations (20/80, 25.00%). Differences in mutation types were observed between DMD/IMD and BMD (P = 0.032; Figure 1). Additionally, carrier status of the probands’ mothers was determined for 78 out of 98 females (79.59%) whose DNA was submitted for Sanger sequencing; carrier status was not significantly different between DMD/IMD and BMD (P = 0.433) or patients with different mutations (P = 0.543; Figure 2).
Figure 1. Distribution of mutation types. The ratios of mutation types in (A) all patients, (B) DMD/IMD, and (C) BMD.
Figure 2. Carrier status of probands’ mothers. Carrier status ratios for mothers of (A) all probands, (B) DMD/IMD, and (C) BMD as well as probands with the following different mutation types, (D) nonsense, (E) frameshift, (F) splicing, and (G) missense mutations.
We also determined the distribution of mutations in DMD (Figure 3A) and found that nearly 50% of patients had mutations located in the 5′ region upstream of exon 30 (4071 bp; 36.82% of the gene) in the DMD coding sequence. Additionally, the distribution of mutation types in DMD domains differed significantly (P = 0.002; Figures 3B,C); no splicing and frameshift mutations were observed in the CRD, and no missense mutations were observed in the RD or CTD.
Figure 3. Positional distributions of small mutations. (A) Distribution of exons in which mutations are located. (B) The ratios of mutation types in DMD domains. (C) Distribution of domains in which mutations are located. Distributions of splicing mutation were classified to consensus exons.
Nonsense mutations (60/115) were the most common across all mutation types (52.17%) and when truncating mutations were isolated (60/80, 75.00%). Patients with transition mutations (41/60, 68.33%) were more common than those with transversion mutations (19/60, 31.67%). C-to-T (37/41, 90.24%) and G-to-A (4/41, 9.76%) were the only transition types observed. Similar to the results of a previous study (Tuffery-Giraud et al., 2009), the frequency of patients with mutations located in CpG dinucleotides was 41.67% (25/60; consisting of 23 C-to-T, 1 C-to-A, and 1 G-to-A). These results indicate that CpG islands are a mutational hot spot, likely due to oxidative deamination of 5-methyl cytosine (Tuffery-Giraud et al., 2009). However, the distribution of CpG mutations observed in this study differed from what was previously reported (Figure 4). The most common CpG mutation we observed was c. 10141C > T (4/25, 16.00%) in exon 70, followed by c. 5287 C > T (2/25, 8.00%) in exon 37, and c. 6292 C > T (2/25, 8.00%) in exon 44. Furthermore, we observed mutations at CpG dinucleotides located in exon 9 (c.883C > T), exon 32 (c.4375C > T), and exon 61 (c.9100C > T), which were not reported in a previous study (Tuffery-Giraud et al., 2009).
Figure 4. Positional distribution of nonsense mutations located in CpG dinucleotides and other regions.
Nonsense mutation sites were primarily clustered at the first codon position (49/60, 81.67%), followed by the third (7/60, 11.67%), and second position (4/60, 6.66%). No G-to-T substitutions at the first codon position were located at the CpG dinucleotide; however, C-to-T at the first codon position were more often observed at CpG dinucleotides than at any other site [23 mutations (62.16%) vs. 14 mutations (37.84%), respectively]. Analysis of PTC types observed in nonsense mutations revealed that TGA was the most common (29/60, 48.33%), followed by TAG (15/60, 25.00%), and TAA (16/60, 26.67%). Moreover, there was a statistically significant difference between these PTC types among the different substitutions observed at the first codon position (P < 0.001; Table 2).
Frameshift mutations accounted for 25.00% of all truncating mutations. Within this group, there were 14 small deletions (70.00%) and 6 small insertions (30.00%; Figure 5A). Twelve mutations (60.00%) resulted from either small deletions or insertions of a single base; however, no small insertions of several bases were observed (Figure 5B). Apart from a small insertion (c.1376dupA) and deletion (c.1533_1536del) in exon 12, no exon contained both a small insertion and deletion. The mRNA of transcripts containing frameshift mutations were degraded by nonsense-mediated decay due to the presence of a hidden downstream PTC (Miller and Pearce, 2014). Analysis of the first hidden PTC downstream of the mutated sites identified 11 TGA (55.00%), 5 TAA (25.00%), and 4 TAG (20.00%) codons. These results did not differ from the PTC distribution results in the nonsense mutation group (P > 0.05; Figure 6). However, distribution of the first hidden PTC in small insertions and deletions differed but was not statistically significant (P = 0.361; Figure 5C). The positional distributions were significantly different between small insertions and deletions (P = 0.003); all small deletions were located in the RD, while two small insertions were observed in each of the ABD, RD, and CTD.
Figure 5. Information regarding small insertions and deletions in frameshift mutations. (A) Proportions of small insertions and deletions in frameshift mutations. (B) The number of deleted or inserted bases in small insertions and deletions. (C) Proportions of first hidden premature termination codons in small insertions and deletions.
Figure 6. Ratios of premature termination codons. The constituent ratios of premature termination codons in (A) truncating, (B) nonsense, and (C) frameshift mutations.
In this study, we defined splicing mutations as those within introns that led to DMD splicing changes. Out of 28 patients with splicing mutations, 21 mutations (75.00%) were located at a highly conserved splice site (±1–2 bp from the exon-intron boundary), whereas, the others (5/28, 17.86%) were located 3–5 bp from this boundary. In addition, we identified two splicing mutations, c.9164-3_9164-1delCAG and c.10329-11_10330del, that spanned highly conserved splice regions. Of all splicing mutations, substitutions were the most common mutation type observed (24/28, 85.72%), followed by small insertions (2/28, 7.14%), and deletions (2/28, 7.14%). Furthermore, there were 14 splicing mutations each (50.00%) in the splice acceptor and donor sites. The distribution of splicing mutations is represented in Figure 7.
There were six missense mutations (5.22%) identified in this study, which is consistent with the percentages reported in previous studies (1–5%; Sedlácková et al., 2009; Tuffery-Giraud et al., 2009; Juan-Mateu et al., 2015). Four of these mutations were located in the ABD and the other two in the CRD, consistent with the clustering phenomenon of missense mutation in ABD and CRD demonstrated in previous reports (Tuffery-Giraud et al., 2009; Juan-Mateu et al., 2015). Although the RD accounts for 75.92% of DMD, missense mutations were rare in this domain. Compared with the rate of missense mutations, the proportions of other mutations in the RD correlate with its length (nonsense mutations, 76.67%; frameshift mutations, 80.00%; and splicing mutations, 71.43%; Figure 3C). This phenomenon may be associated with the spectrin-like repeat structure in RD.
Dystrophin Expression in the Muscles
Twenty-five patients underwent muscle biopsy; eight presented with patchy dystrophin staining in muscles, whereas the others showed absence of dystrophin expression. Patchy staining represents dystrophin expression levels > 10%. Among the 8 patients with patchy staining, 7 manifested BMD and 1 was pending. For patients with absence of dystrophin expression, 15 manifested DMD/IMD, 1 BMD, and 1 was pending (see Supplementary Table S1 for detailed information).
Levels of SCRN Differ in Patients With DMD/IMD and BMD
To detect whether SCRN and CK levels differed between phenotypes, we set phenotypes as dependent variables while age of examination, mutation type and SCRN, or CK were independent variables and analyzed via logistic regression. Because all patients above the age of 16 presented with BMD, only ≤16 years old (n = 86) were included in this analysis. Furthermore, 10 subjects (11.63%) with missing SCRN data were excluded from this analysis. The results indicated that SCRN (P = 0.002, OR = 1.934, CI: 1.266–2.956) levels differed between DMD/IMD and BMD, whereas CK levels did not (P > 0.05). Moreover, SCRN levels did not differ between mutation types (P > 0.05; Figure 8).
Figure 8. Correlations between phenotypes and serum creatinine (SCRN) or creatine kinase (CK) levels. (A) SCRN levels in patients with Duchenne muscular dystrophy (DMD)/intermediate muscular dystrophy (IMD) or Becker muscular dystrophy (BMD). (B) CK levels in patients with DMD/IMD or BMD. The median of SCRN levels at each age is represented by short lines.
Dystrophinopathy phenotypes are dependent on the degree to which the dystrophin gene, DMD, is affected by the mutation. Complete loss of DMD function leads to DMD, while relatively mild disruption of DMD function leads to BMD (Monaco et al., 1988). Our results indicate that fewer DMD/IMD patients had non-truncating mutations (15/66, 22.73%) compared with those with BMD (15/31, 48.39%; P = 0.011, OR = 0.314, CI: 0.126–0.779). These results support the theory that mutations with a greater effect on protein function result in more severe phenotypes because truncating mutations usually lead to degeneration of mRNA and complete absence of the protein (Miller and Pearce, 2014).
We further analyzed the correlation between phenotypes and the protein domains in which the mutations were located (Table 3). We found no correlation between phenotypes and domains in which mutations are located for splicing and missense mutations (P > 0.05). However, all the patients with frameshift mutations in RD presented with DMD/IMD, which was significant (P = 0.033).
All nonsense mutations identified in BMD were located in the RD (Table 3). Given the high proportion of exons whose deletions lead to in-frame mutations (30/54, 55.56%) and the probable mild phenotypes resulting from in-frame deletions at RD (Yokota et al., 2012), spontaneous exon skipping may account for BMD with nonsense mutations. To assess the applicability of splicing classification grades for predicting spontaneous exon skipping in nonsense mutations, we analyzed the correlation between splicing classification grades and expression levels of the dystrophin C-terminus. All nonsense mutations with splicing classification grade 0 exhibited absent dystrophin levels, whereas half of nonsense mutations with splicing classification grade 1 presented with patchy dystrophin C-terminus staining (Figure 9). Thus, our results suggest that nonsense mutations with splicing classification grade 0 did not lead to beneficial spontaneous exon skipping.
Figure 9. Correlation between splicing classification grades and dystrophin levels in patients with nonsense mutations. Expression levels of dystrophin C-terminus were evaluated using muscle immunohistochemistry with antibody against this terminus. Two staining conditions were observed in the muscle sections: (A) absence staining and, (B) patchy staining. (C) Bar chart representing percentage of patients versus splicing classification grade.
Furthermore, we found that nonsense mutations with splicing classification grade 0 most frequently represented DMD/IMD (26/28, 92.86%, P < 0.001, OR = 12.000, CI: 2.332–61.758). All nonsense mutations with splicing classification grade 0 at the ABD, CRD, and CTD, were associated with DMD/IMD (Table 4). Analysis of nonsense mutations at the RD revealed that those with splicing classification grade 0 also most frequently represented DMD/IMD (15/17, 88.24%, P = 0.011, OR = 7.500, CI: 1.400–40.178). Thus, splicing classification grades were helpful in distinguishing DMD/IMD from BMD with nonsense mutations.
Table 4. Correlation between splicing classification grade and phenotypes in nonsense mutation cases.
Exon Skipping Strategy for Patients With Dystrophinopathies and Small Mutations
Exon skipping is a promising treatment for dystrophinopathies and can target small mutations by skipping the exons in which the mutations are located (Aartsma-Rus et al., 2009). If skipping a single exon containing the mutation resulted in a frameshift, adjacent exons can also be skipped to maintain the ORF (Yokota et al., 2012). Previous studies confirmed that exon skipping therapy can restore dystrophin in dystrophic mice and in muscle cells of patients with DMD small mutations (Aartsma-Rus et al., 2003; Fletcher et al., 2007; Spitali et al., 2009). In this study, we included truncating (nonsense and frameshift mutations) and missense mutations for our analysis of exon skipping therapy. The rationale for this is that splicing mutations can lead to spontaneous exon skipping or intron retention, which makes exon skipping strategies unreliable depending on the specific DNA changes (Thi Tran et al., 2005; Hagenbuch et al., 2010).
We found that single-exon skipping, single or double-exon (double/single-exon) skipping, and single, double, or triple-exon (triple/double/single-exon) skipping can cover 54.65, 87.21, and 96.51% of patients with truncating and missense mutations, respectively (Figure 10A). The aim of exon skipping treatment is to change severe phenotypes to milder ones. Theoretically, as BMD already represents the mildest phenotype, patients with BMD should not receive this treatment (Aartsma-Rus et al., 2009). Thus, after excluding BMD and pending subjects from our analysis, we found that single-exon, double/single-exon, and triple/double/single-exon skipping can cover 50.00, 92.31, and 98.08% of DMD/IMD patients with truncating and missense mutations (n = 52), respectively (Figure 10B). Excluding mutations representing BMD resulted in different exon skipping strategies with more coverage than that observed when BMD was included (Figures 10C,D).
Figure 10. Exon skipping strategies for truncating and missense mutations. The covering proportions of potentially applicable patients among (A) all patients, and (B) those with Duchenne muscular dystrophy (DMD)/intermediate muscular dystrophy (IMD). Exon skipping strategies for (C) all patients and (D) those with DMD/IMD.
The DMD gene is the largest known human gene, containing 79 exons and is 2.4 Mb in length (Kole and Krieg, 2015). MLPA for detecting DMD large rearrangements is an already proven practical and reliable method (Aartsma-Rus et al., 2016). On the other hand, sanger sequencing for detecting small mutations is costly and time-consuming due to the DMD gene length, making muscle biopsy important for diagnosis of dystrophinopathies (Bushby et al., 2010). In recent years, however, the development of NGS made the detection of DMD small mutations more cost-effective and accurate than before. Furthermore, since obtaining information about mutations is critical for prenatal diagnosis, preimplantation genetic diagnosis and administration of gene therapies (Jarmin et al., 2014), genetic testing is now performed prior to muscle biopsy for diagnosing dystrophinopathies, as stated in the latest guideline (Birnkrant et al., 2018). Considering the expensive cost of directly combining MPLA and sequencing, we initially perform MLPA for large rearrangements, which covers about 80% of patients. If the results are negative, we perform sequencing, thereby cutting down on cost. After the confirmation of probands’ mutations, the probands’ mothers should be verified for the carrier state. As sequencing technology develops and becomes cheaper, large rearrangements and small mutations can be detected simultaneously in the future.
Although muscle biopsies for detecting dystrophin expression are not critical for the diagnosis of dystrophinopathies nowadays, quantification of dystrophin expression levels is helpful in classifying phenotypes. Duchenne patients have little to no dystrophin in muscle tissue, while dystrophin levels in IMD and BMD are higher (Hoffman et al., 1988; Fanin et al., 1992, 1995; Angelini et al., 1994). Furthermore, low levels of dystrophin expression in female carriers suggests a poor prognosis (Pegoraro et al., 1995). Thus, dystrophin expression levels are correlated with clinical phenotypes. Our results were compatible with this theory, with the exception of a BMD patient that lacked dystrophin expression. Previous studies that included patients and DMD dog models also reported some exceptions (Castro-Gago, 2015; Vieira et al., 2015). It is difficult to explain this phenomenon, but it is possible that a supplementary mechanism is responsible for the dystrophin defect, which should be investigated further and may be a novel target for treating patients with dystrophinopathies.
Creatinine is a metabolic product of the creatine pathway, which is important for myofiber energy metabolism (Wyss and Kaddurah-Daouk, 2000). We found that SCRN levels in BMD patients younger than 16 years were higher than those in DMD/IMD, which was independent of mutation types. Furthermore, the difference of SCRN levels between DMD/IMD and BMD showed an increasing trend in an age-dependent manner. Therefore, SCRN levels correlate with clinical phenotypes, rather than genotypes, confirming the results of our previous small sample study (Wang et al., 2017). CK is a typical biomarker in neuromuscular disorders and whose levels are extremely elevated in dystrophinopathies (Pernice et al., 1986). Although CK is a good biomarker for screening and diagnosis, analysis of its levels could not distinguish DMD/IMD from BMD, which may result from the large variability in CK levels.
Regarding dystrophinopathy mutations, the characteristics revealed in this study are consistent with those of previous studies, including proportion of mutation and PTC types, frequency of nonsense mutations in CpG dinucleotides, and the ratios of small deletions and insertions (Sedlácková et al., 2009; Tuffery-Giraud et al., 2009; Magri et al., 2011; Juan-Mateu et al., 2015). The most common mutation type was nonsense mutations, which suggests that nonsense read-through therapies hold promise in China. However, we identified some special characteristics in our cohort. For example, base changes in the first codon position were correlated with PTC types in nonsense mutations, and positional distributions differed between small deletions and insertions.
Furthermore, the positional distribution of small mutations observed in this study differed from those reported previously (Flanigan et al., 2015; Juan-Mateu et al., 2015). In general, small mutations tended to cluster in the 5′ region of the DMD ORF. Additionally, the most common positions for small mutations were exon/intron 32, 40, and 44, which were different from those reported in other studies (Sedlácková et al., 2009; Flanigan et al., 2015; Guo et al., 2015; Juan-Mateu et al., 2015; Li et al., 2015). This difference can be observed between different countries, such as China, Spain, America, and the Czechia (Sedlácková et al., 2009; Flanigan et al., 2015; Juan-Mateu et al., 2015). As for the positional distribution of nonsense mutation, the top three in our cohort were exon 70, 32, and 44, similarly, different characteristics of distribution can be observed between different countries (Sedlácková et al., 2009; Tuffery-Giraud et al., 2009; Flanigan et al., 2015; Guo et al., 2015; Juan-Mateu et al., 2015; Li et al., 2015). However, the proportion of patients with nonsense mutations in exon 70 was the highest in almost all studies, except for a Spanish cohort (Sedlácková et al., 2009; Tuffery-Giraud et al., 2009; Flanigan et al., 2015; Guo et al., 2015; Juan-Mateu et al., 2015; Li et al., 2015). Thus, as the distribution characteristics of small mutations are distinguished in different regions, different strategies should be tailored for targeted genetic treatment.
We examined the utility of exon skipping therapy as an example of such targeted treatment. Our results indicated that exon skipping therapy can be used in truncating and missense mutation cases because single-exon skipping can cover half of all patients with these mutations while triple/double/single-exon skipping can cover 96.51%; this is consistent with previous reports (Aartsma-Rus et al., 2009; Yokota et al., 2012). On the other hand, the most frequent single-exon skipping strategy was exon 32 skipping (covering 8.51% of total patients and 11.54% of DMD/IMD patients with truncating and missense mutations), which was different from a previous study that only focused on nonsense mutations (Yokota et al., 2012). In theory, multiexon skipping can cover many more patients than can single-exon skipping; however, its application is challenging owing to toxicity and non-specificity (Aartsma-Rus et al., 2014). Thus, if strategies were developed to eliminate the associated toxicity and increase specificity, a multi-exon skipping strategy would be the preferred choice in the treatment of small mutations. All mutations that could not be treated by triple/double/single-exon skipping were located in exon 2, which might require exon 2–19 skipping (Yokota et al., 2012).
Theoretically, patients with BMD should not receive exon skipping therapy because the treatment aim is to reduce the phenotype severity of DMD to that of BMD (Aartsma-Rus et al., 2009). Therefore, it is important to evaluate the specific phenotype of patients before initiating treatment. The classical method of phenotypic classification involves evaluating motor dysfunction and measuring dystrophin levels in the muscles (Arahata et al., 1989). However, these classical methods have some disadvantages, such as inaccuracies for young patients and invasiveness, which limit their application (Wang et al., 2017). Our results suggest two supplementary methods that can be used to more accurately determine dystrophinopathy phenotypes. First, we showed that patients with truncating mutations are more likely to present as DMD/IMD. Moreover, patients with nonsense mutations with a splicing classification grade 0 have an increased likelihood of presenting as DMD/IMD (92.86%); this is because nonsense mutations with beneficial spontaneous exon skipping can convert severe phenotypes to milder phenotypes (Flanigan et al., 2011). We believe these findings may promote the clinical application of bioinformatic tools. Second, SCRN levels can be helpful in determining phenotypic classifications (Wang et al., 2017), and remains practicable for small mutations.
In some circumstances, such as for very young patients, phenotypes cannot be determined precisely, which makes the exon skipping strategy of initially analyzing total patients with dystrophinopathies. However, once the patient can be classified into a specific phenotype with the help of clinical manifestations, biomarkers, bioinformatic analysis, and other methods, BMD patients should be excluded from exon skipping therapy. Thus, we analyzed exon skipping strategies in all patients as well as in DMD/IMD patients only. The results indicated that the covering rates of exon skipping strategies for DMD/IMD were similar to those calculated when taking into account all patients. However, the exon skipping strategies with the greatest covering frequencies differed when patients with BMD and pending were excluded from the analysis. These strategies are important and thus close attention is required when medications are designed. After comprehensive analysis, our results indicate that exon skipping strategies as treatment options for small mutations in Chinese patients should focus on exon 32 for single-exon skipping, exons 69–70 for double-exon skipping, and exons 69–71 for triple-exon skipping. All three strategies were the most common across all patients and patients with DMD/IMD.
Finally, there are some limitations to this study. First, deep regions in introns could not be detected due to the limitations of current sequencing technologies. Therefore, the frequency of splicing mutations may be underestimated. Second, although we performed a preliminary investigation of spontaneous exon skipping for nonsense mutations using immunohistochemistry staining of muscle dystrophin, the number of patients that underwent biopsies was limited due to its invasiveness; thus, this topic must be studied further with a larger sample size. Furthermore, the quantification of muscle dystrophin levels detected by western blotting is helpful for classifying phenotypes, which was not applied in this study due to the few number of available muscle samples. Third, the number of patients under the age of 3 subjected to the CK and SCRN test was limited because patients at these ages are pauci-symptomatic and doctors are rarely consulted during this time.
In this study, we identified 106 small mutations in 115 unrelated patients with dystrophinopathies and reported the mutation characteristics and genotype–phenotype correlations. These characteristics suggested a mutation pattern in the DMD gene. Furthermore, compared with other non-splicing mutations, nonsense mutations tended to result in splicing changes, which provided evidence for spontaneous exon skipping in BMD patients with nonsense mutation. One BMD patient exhibited a lack of dystrophin expression, suggesting a novel mechanism for phenotype modification that should be investigated further.
The novel mutations uncovered here will be helpful in the diagnosis of dystrophinopathies and the use of splicing classification grades for predicting phenotypes, indicating a practical clinical application for bioinformatic tools. Moreover, the classification of phenotypes with the help of splicing classification grades and SCRN levels can aid the judgment of prognosis and in the design of treatment strategies. Data regarding mutation distribution and the investigation of exon skipping strategies may offer information for drug design, as 96–98% of patients with truncating or missense mutations were covered by triple/double/single-exon skipping, supporting the potential of exon skipping therapy for dystrophinopathy treatment.
Thus, our study findings have expanded our knowledge and understanding of mutation characteristics and genotype–phenotype correlations in dystrophinopathies, which can provide further insights regarding its pathogenesis, diagnosis, and treatment design.
The study was conducted according to the ethical guidelines of the Declaration of Helsinki (1975), and a waiver of informed parental or personal consent was approved by the ICE for Clinical Research and Animal Trials of the First Affiliated Hospital of Sun Yat-sen University.
LW and MX designed the study, analyzed the data, and drafted the manuscript. LW, HL, YZ, RH, and JL recorded clinical and genetic data. YZ and CZ supervised the study and assisted in drafting and revising the manuscript.
This study was sponsored by grants from the Natural Science Foundation of China (Grant Nos. 81771359, 81471280, 81601087, and 81271401), the Guangdong Provincial Science and Technology Plan (Grant No. 2014A020212130), the Guangzhou City Science and Technology Plan (Grant No. 1561000153/201508020012), the China Scholarship Council (Grant No. 201706380141), the National Key Clinical Department and Key Discipline of Neurology, the Guangdong Provincial Key Laboratory for Diagnosis and Treatment of Major Neurological Diseases (Grant No. 2014B030301035), the Southern China International Cooperation Base for Early Intervention and Functional Rehabilitation of Neurological Diseases (Grant No. 2015B050501003), and the Guangdong Provincial Engineering Center for Major Neurological Disease Treatment, Science, and Technology Planning Project of Guangzhou (Grant No. 201604020010).
Conflict of Interest Statement
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.
The authors would like to thank Prof. Victor Wei Zhang (Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, United States) and Shaoyuan Li (AmCare Genomics Lab, Guangzhou, China) for their assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00114/full#supplementary-material
Aartsma-Rus, A., Ferlini, A., Goemans, N., Pasmooij, A. M., Wells, D. J., Bushby, K., et al. (2014). Translational and regulatory challenges for exon skipping therapies. Hum. Gene. Ther. 25, 885–892. doi: 10.1089/hum.2014.086
Aartsma-Rus, A., Fokkema, I., Verschuuren, J., Ginjaar, I., van Deutekom, J., van Ommen, G. J., et al. (2009). Theoretic applicability of antisense-mediated exon skipping for duchenne muscular dystrophy mutations. Hum. Mutat. 30, 293–299. doi: 10.1002/humu.20918
Aartsma-Rus, A., Janson, A. A., Kaman, W. E., Bremmer-Bout, M., den Dunnen, J. T., Baas, F., et al. (2003). Therapeutic antisense-induced exon skipping in cultured muscle cells from six different DMD patients. Hum. Mol. Genet. 12, 907–914. doi: 10.1093/hmg/ddg100
Aartsma-Rus, A., Van Deutekom, J. C., Fokkema, I. F., Van Ommen, G. J., and Den Dunnen, J. T. (2006). Entries in the Leiden Duchenne muscular dystrophy mutation database: an overview of mutation types and paradoxical cases that confirm the reading-frame rule. Muscle Nerve 34, 135–144. doi: 10.1002/mus.20586
Angelini, C., Fanin, M., Pegoraro, E., Freda, M. P., Cadaldini, M., and Martinello, F. (1994). Clinical-molecular correlation in 104 mild X-linked muscular dystrophy patients: characterization of sub-clinical phenotypes. Neuromuscul Disord. 4, 349–358. doi: 10.1016/0960-8966(94)90071-X
Arahata, K., Hoffman, E. P., Kunkel, L. M., Ishiura, S., Tsukahara, T., Ishihara, T., et al. (1989). Dystrophin diagnosis: comparison of dystrophin abnormalities by immunofluorescence and immunoblot analyses. Proc. Natl. Acad. Sci. U.S.A. 86, 7154–7158. doi: 10.1073/pnas.86.18.7154
Birnkrant, D. J., Bushby, K., Bann, C. M., Apkon, S. D., Blackwell, A., Brumbaugh, D., et al. (2018). Diagnosis and management of Duchenne muscular dystrophy, part 1: diagnosis, and neuromuscular, rehabilitation, endocrine, and gastrointestinal and nutritional management. Lancet Neurol. 17, 251–267. doi: 10.1016/S1474-4422(18)30024-3
Brooke, M. H., Griggs, R. C., Mendell, J. R., Fenichel, G. M., Shumate, J. B., Pellegrino, R. J., et al. (1981). Clinical trial in Duchenne dystrophy. I. The design of the protocol. Muscle Nerve 4, 186–197. doi: 10.1002/mus.880040304
Bushby, K., Finkel, R., Birnkrant, D. J., Case, L. E., Clemens, P. R., Cripe, L., et al. (2010). Diagnosis and management of Duchenne muscular dystrophy, part 1: diagnosis, and pharmacological and psychosocial management. Lancet Neurol. 9, 77–93. doi: 10.1016/S1474-4422(09)70271-6
Chen, W. J., Lin, Q. F., Zhang, Q. J., He, J., Liu, X. Y., Lin, M. T., et al. (2013). Molecular analysis of the dystrophin gene in 407 Chinese patients with Duchenne/Becker muscular dystrophy by the combination of multiplex ligation-dependent probe amplification and Sanger sequencing. Clin. Chim. Acta 423, 35–38. doi: 10.1016/j.cca.2013.04.006
Desmet, F. O., Hamroun, D., Lalande, M., Collod-Béroud, G., Claustres, M., and Béroud, C. (2009). Human splicing finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 37:e67. doi: 10.1093/nar/gkp215
Doriguzzi, C., Palmucci, L., Mongini, T., Chiadò-Piat, L., Maniscalco, M., and Restagno, G. (1997). Systematic use of dystrophin testing in muscle biopsies results in 201 cases. Eur. J. Clin. Invest. 27, 352–358. doi: 10.1046/j.1365-2362.1997.1090650.x
Fanin, M., Danieli, G. A., Cadaldini, M., Miorin, M., Vitiello, L., and Angelini, C. (1995). Dystrophin-positive fibers in Duchenne dystrophy: origin and correlation to clinical course. Muscle Nerve 18, 1115–1120. doi: 10.1002/mus.880181007
Fanin, M., Danieli, G. A., Vitiello, L., Senter, L., and Angelini, C. (1992). Prevalence of dystrophin-positive fibers in 85 Duchenne muscular dystrophy patients. Neuromuscul. Disord. 2, 41–45. doi: 10.1016/0960-8966(92)90025-2
Flanigan, K. M., Dunn, D. M., von Niederhausern, A., Soltanzadeh, P., Gappmaier, E., Howard, M. T., et al. (2015). Mutational spectrum of DMD mutations in dystrophinopathy patients: application of modern diagnostic techniques to a large cohort. Hum. Mutat. 30, 1657–1666. doi: 10.1002/humu.21114
Flanigan, K. M., Dunn, D. M., von Niederhausern, A., Soltanzadeh, P., Howard, M. T., Sampson, J. B., et al. (2011). Nonsense mutation-associated Becker muscular dystrophy: interplay between exon definition and splicing regulatory elements within the DMD gene. Hum. Mutat. 32, 299–308. doi: 10.1002/humu.21426
Fletcher, S., Honeyman, K., Fall, A. M., Harding, P. L., Johnsen, R. D., Steinhaus, J. P., et al. (2007). Morpholino oligomer-mediated exon skipping averts the onset of dystrophic pathology in the mdx mouse. Mol. Ther. 15, 1587–1592. doi: 10.1038/sj.mt.6300245
Ginjaar, I. B., Kneppers, A. L., v d Meulen, J. D., Anderson, L. V., Bremmer-Bout, M., van Deutekom, J. C., et al. (2000). Dystrophin nonsense mutation induces different levels of exon 29 skipping and leads to variable phenotypes within one BMD family. Eur. J. Hum. Genet. 8, 793–796. doi: 10.1038/sj.ejhg.5200535
Guo, R., Zhu, G., Zhu, H., Ma, R., Peng, Y., Liang, D., et al. (2015). DMD mutation spectrum analysis in 613 Chinese patients with dystrophinopathy. J. Hum. Genet. 60, 435–442. doi: 10.1038/jhg.2015.43
Hagenbuch, S. C., Gottliebson, W. M., Wansapura, J., Mazur, W., Fleck, R., Benson, D. W., et al. (2010). Detection of progressive cardiac dysfunction by serial evaluation of circumferential strain in patients with Duchenne muscular dystrophy. Am. J. Cardiol. 105, 1451–1455. doi: 10.1016/j.amjcard.2009.12.070
Hoffman, E. P., Fischbeck, K. H., Brown, R. H., Johnson, M., Medori, R., Loike, J. D., et al. (1988). Characterization of dystrophin in muscle-biopsy specimens from patients with Duchenne’s or Becker’s muscular dystrophy. N. Engl. J. Med. 318, 1363–1368. doi: 10.1056/NEJM198805263182104
Jarmin, S., Kymalainen, H., Popplewell, L., and Dickson, G. (2014). New developments in the use of gene therapy to treat Duchenne muscular dystrophy. Expert. Opin. Biol. Ther. 14, 209–230. doi: 10.1517/14712598.2014.866087
Juan-Mateu, J., Gonzalez-Quereda, L., Rodriguez, M. J., Baena, M., Verdura, E., Nascimento, A., et al. (2015). DMD mutations in 576 dystrophinopathy families: a step forward in genotype-phenotype correlations. PLoS One 10:e0135189. doi: 10.1371/journal.pone.0135189
Li, X., Zhao, L., Zhou, S., Hu, C., Shi, Y., Shi, W., et al. (2015). A comprehensive database of Duchenne and Becker muscular dystrophy patients (0-18 years old) in East China. Orphanet. J. Rare. Dis. 10:5. doi: 10.1186/s13023-014-0220-7
Magri, F., Del Bo, R., D’Angelo, M. G., Govoni, A., Ghezzi, S., Gandossini, S., et al. (2011). Clinical and molecular characterization of a cohort of patients with novel nucleotide alterations of the Dystrophin gene detected by direct sequencing. BMC Med. Genet. 12:37. doi: 10.1186/1471-2350-12-37
Mah, J. K., Korngut, L., Dykeman, J., Day, L., Pringsheim, T., and Jette, N. (2014). A systematic review and meta-analysis on the epidemiology of Duchenne and Becker muscular dystrophy. Neuromuscul. Disord. 24, 482–491. doi: 10.1016/j.nmd.2014.03.008
Marden, F. A., Connolly, A. M., Siegel, M. J., and Rubin, D. A. (2005). Compositional analysis of muscle in boys with Duchenne muscular dystrophy using MR imaging. Skeletal. Radiol. 34, 140–148. doi: 10.1007/s00256-004-0825-3
Monaco, A. P., Bertelson, C. J., Liechti-Gallati, S., Moser, H., and Kunkel, L. M. (1988). An explanation for the phenotypic differences between patients bearing partial deletions of the DMD locus. Genomics 2, 90–95. doi: 10.1016/0888-7543(88)90113-9
Pegoraro, E., Schimke, R. N., Garcia, C., Stern, H., Cadaldini, M., Angelini, C., et al. (1995). Genetic and biochemical normalization in female carriers of Duchenne muscular dystrophy: evidence for failure of dystrophin production in dystrophin-competent myonuclei. Neurology 45, 677–690. doi: 10.1212/WNL.45.4.677
Pernice, W., Guggolz, M. A., Guggolz, M., Beckmann, R., and Wais, U. (1986). A mathematical analysis of creatine kinase activity in the course of Duchenne muscular dystrophy. Muscle Nerve 9, 333–340. doi: 10.1002/mus.880090409
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, 405–424. doi: 10.1038/gim.2015.30
Sedlácková, J., Vondrácek, P., Hermanová, M., Zámecník, J., Hrubá, Z., Haberlová, J., et al. (2009). Point mutations in czech DMD/BMD patients and their phenotypic outcome. Neuromuscul. Disord. 19, 749–753. doi: 10.1016/j.nmd.2009.08.011
Spitali, P., Rimessi, P., Fabris, M., Perrone, D., Falzarano, S., Bovolenta, M., et al. (2009). Exon skipping-mediated dystrophin reading frame restoration for small mutations. Hum. Mutat. 30, 1527–1534. doi: 10.1002/humu.21092
Thi Tran, H. T., Takeshima, Y., Surono, A., Yagi, M., Wada, H., and Matsuo, M. (2005). A G-to-A transition at the fifth position of intron-32 of the dystrophin gene inactivates a splice-donor site both in vivo and in vitro. Mol. Genet. Metab. 85, 213–219. doi: 10.1016/j.ymgme.2005.03.006
Tuffery-Giraud, S., Beroud, C., Leturcq, F., Yaou, R. B., Hamroun, D., Michel-Calemard, L., et al. (2009). Genotype-phenotype analysis in 2,405 patients with a dystrophinopathy using the UMD-DMD database: a model of nationwide knowledgebase. Hum. Mutat. 30, 934–945. doi: 10.1002/humu.20976
Vieira, N. M., Guo, L. T., Estrela, E., Kunkel, L. M., Zatz, M., and Shelton, G. D. (2015). Muscular dystrophy in a family of labrador retrievers with no muscle dystrophin and a mild phenotype. Neuromuscul. Disord. 25, 363–370. doi: 10.1016/j.nmd.2015.02.012
Wang, L., Chen, M., He, R., Sun, Y., Yang, J., Xiao, L., et al. (2017). Serum creatinine distinguishes Duchenne muscular dystrophy from Becker muscular dystrophy in patients aged < = 3 years: a retrospective study. Front. Neurol. 8:196. doi: 10.3389/fneur.2017.00196
Yang, J., Li, S. Y., Li, Y. Q., Cao, J. Q., Feng, S. W., Wang, Y. Y., et al. (2013). MLPA-based genotype-phenotype analysis in 1053 Chinese patients with DMD/BMD. BMC Med. Genet. 14:29. doi: 10.1186/1471-2350-14-29
Yokota, T., Duddy, W., Echigoya, Y., and Kolski, H. (2012). Exon skipping for nonsense mutations in Duchenne muscular dystrophy: too many mutations, too few patients? Expert. Opin. Biol. Ther. 12, 1141–1152. doi: 10.1517/14712598.2012.693469
Yu, M., Zheng, Y., Jin, S., Gang, Q., Wang, Q., Yu, P., et al. (2017). Mutational spectrum of Chinese LGMD patients by targeted next-generation sequencing. PLoS One 12:e0175343. doi: 10.1371/journal.pone.0175343
Keywords: dystrophinopathies, small mutation, mutation spectrum, genotype–phenotype correlation, exon skipping therapy
Citation: Wang L, Xu M, Li H, He R, Lin J, Zhang C and Zhu Y (2019) Genotypes and Phenotypes of DMD Small Mutations in Chinese Patients With Dystrophinopathies. Front. Genet. 10:114. doi: 10.3389/fgene.2019.00114
Received: 24 June 2018; Accepted: 30 January 2019;
Published: 18 February 2019.
Edited by:Jordi Pérez-Tur, Instituto de Biomedicina de Valencia (IBV), Spain
Reviewed by:Muhammad Jawad Hassan, National University of Medical Sciences (NUMS), Pakistan
Corrado Angelini, University of Padova, Italy
Copyright © 2019 Wang, Xu, Li, He, Lin, Zhang and Zhu. 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.
†These authors have contributed equally to this work