ORIGINAL RESEARCH article
Sec. Neuromuscular Disorders and Peripheral Neuropathies
Volume 12 - 2021 | https://doi.org/10.3389/fneur.2021.659922
Platelet Derived Growth Factor-AA Correlates With Muscle Function Tests and Quantitative Muscle Magnetic Resonance in Dystrophinopathies
- 1Neuromuscular Disorders Unit, Neurology Department, Hospital de la Santa Creu i Sant Pau, Departament de Medicina. Universitat Autònoma de Barcelona, Barcelona, Spain
- 2Neurology Department, Neuromuscular Reference Center, University Hospital of Antwerp, Antwerp, Belgium
- 3Biomedical Network Research Centre on Rare Diseases (CIBERER), Barcelona, Spain
- 4John Walton Muscular Dystrophy Research Centre, International Centre for Life, Newcastle University, Newcastle upon Tyne, United Kingdom
- 5Neuromuscular Unit, Neuropediatrics Department, Institut de Recerca Sant Joan de Déu, Hospital Sant Joan de Déu, Barcelona, Spain
- 6Rehabilitation and Physiotherapy Department, Hospital de la Santa Creu i Sant Pau, Universitat Autònoma de Barcelona, Barcelona, Spain
- 7Departamento de Genética, Microbiología y Estadística, Universidad de Barcelona, Barcelona, Spain
- 8Radiology Department, Hospital de la Santa Creu i Sant Pau, Universitat Autònoma de Barcelona, Barcelona, Spain
- 9Philips Healthcare Iberia, Madrid, Spain
Introduction: Duchenne (DMD) and Becker (BMD) muscular dystrophy are X-linked muscular disorders produced by mutations in the DMD gene which encodes the protein dystrophin. Both diseases are characterized by progressive involvement of skeletal, cardiac, and respiratory muscles. As new treatment strategies become available, reliable biomarkers and outcome measures that can monitor disease progression are needed for clinical trials.
Methods: We collected clinical and functional data and blood samples from 19 DMD patients, 13 BMD patients, and 66 healthy controls (8 pediatric and 58 adult controls), and blood samples from 15 patients with dysferlinopathy (DYSF) and studied the serum concentration of 4 growth factors involved in the process of muscle fibrosis. We correlated the serum concentration of these growth factors with several muscle function tests, spirometry results and fat fraction identified by quantitative Dixon muscle MRI.
Results: We found significant differences in the serum concentration of Platelet Derived Growth Factor-AA (PDGF-AA) between DMD patients and pediatric controls, in Connective Tissue Growth Factor (CTGF) between BMD patients and adult controls, and in and Transforming Growth Factor- β1 (TGF-β1) between BMD and DYSF patients. PDGF-AA showed a good correlation with several muscle function tests for both DMD and BMD patients and with thigh fat fraction in BMD patients. Moreover, PDGF-AA levels were increased in muscle biopsies of patients with DMD and BMD as was demonstrated by immunohistochemistry and Real-Time PCR studies.
Conclusion: Our study suggests that PDGF-AA should be further investigated in a larger cohort of DMD and BMD patients because it might be a good biomarker candidate to monitor the progression of these diseases.
Duchenne (DMD) and Becker (BMD) muscular dystrophy are X-linked recessive disorders produced by mutations in the DMD gene which encodes the protein dystrophin (1). Dystrophin is a cytoplasmic protein which links the intracellular cytoskeleton network to transmembrane components of the dystrophin-glycoprotein complex, including dystroglycan, sarcoglycans, and sarcospan (2). Its role is to stabilize the membrane during muscle contraction to prevent contraction-induced damage (3).
DMD has an homogeneous clinical picture characterized by early onset of muscle weakness progressing during childhood and leading to loss of ambulation during adolescence (4). All DMD patients invariably develop cardiac and respiratory involvement, which are the main causes of death. In contrast, BMD is an heterogeneous disease characterized by muscle weakness of variable severity that can start during childhood but also late in adulthood (5). BMD patients can also develop respiratory and cardiac muscle involvement. Despite having substantial phenotypic differences, DMD and BMD share common pathological pathways. Mutations in the DMD gene lead to abnormal dystrophin expression which results in a fragile muscle membrane susceptible to be damaged during contraction (3, 6). Continuous muscle damage leads to persistent inflammation, muscle fiber necrosis and replacement by fat and fibrous tissue, producing permanent weakness and disability (6).
In slow progressive diseases, the identification of functional changes over a short period is not an easy task. In recent years several outcome measures such as the 6 Minutes-Walk Test (6MWT), the Timed Up-and-Go Test (TuGo) or the North Star Ambulatory Assessment (NSAA) have been developed to monitor muscle function in natural history studies and clinical trials (7–10). Additionally, specific tests have been developed to assess upper limb muscle function such as the Brooke Upper Extremity Functional Rating (Brooke scale) or the Performance of the Upper Limb (PUL) (11–13). However, as muscle fiber necrosis and replacement by fat is a slowly progressing process, changes in motor functional tests could take time to be evident. In view of this limitation, muscle magnetic resonance imaging (MRI) and serum biomarkers have been investigated as potential outcome measures that might complement muscle function tests (14–16). Although muscle MRI is becoming more widely accepted as a tool to monitor patients, there is not a serum biomarker accepted so far that correlates with several muscle function tests and that could therefore be considered useful for clinical trials or natural history studies.
The main aim of this study was to analyze the serum concentration of a group of growth factors related to the process of muscle fibrosis in a cohort of 19 patients with DMD and 13 patients with BMD and study whether there was a correlation with the results of muscle function tests and quantitative MRI (qMRI).
Materials and Methods
Study Design and Participants
We report the results of a transversal study of a cohort of 19 DMD and 13 BMD patients seen at Hospital de la Santa Creu i Sant Pau (HSCSP) and Hospital Sant Joan de Déu (HSJD) in Barcelona. All patients included in the study were assessed using muscle function tests and spirometry and filled out a daily life activities questionnaire. qMRI was not a mandatory test and was performed in a subset of patients followed-up in HSCSP only. Blood samples were obtained before the muscle function assessment. The HSCSP and HSJD Ethics Committees approved the study, and all participants signed an informed consent form. All study procedures were performed in accordance with Spanish regulations.
Inclusion criteria for the study were: (1) Diagnosis of DMD or BMD confirmed by the identification of pathogenic variants in the DMD gene and (2) willingness to complete muscle function tests, respiratory assessment, and daily life activities scale.
Muscle function tests were performed by physiotherapists with considerable experience in neuromuscular disorders. We assessed patients using timed tests, functional ability scales, and measurement of muscle strength. Timed tests included the 6MWT, the 10 Meter Run/Walk test (10MWT), the Time to Climb Up (Tup4) and Down (Tdo4) Four Steps and the Time to Rise from the Floor (TRF). Functional ability included the NSAA and the Motor Function Measure 20-item scale (MFM-20) (17). The MFM-20 was also analyzed based on its 3 dimensions: D1 (standing and transfers), D2 (axial and proximal), and D3 (distal). PUL test was used to assess upper extremity function in both ambulant and non-ambulant individuals. Egen Klassification was used to assess mobility in wheelchair-bound patients (18). Muscle strength was studied using the Medical Research Council scale (MRC). For analysis, 5-point MRC power grades were converted to an 11-point scale as previously described (19). We assessed bilateral shoulder flexion, shoulder abduction, elbow flexion, elbow flexion brachioradialis, wrist extension, wrist flexion, hip adduction, hip abduction, hip flexion, hip extension, knee extension, knee flexion, ankle eversion, ankle inversion, ankle plantar flexion, and ankle plantar extension. Grip and pinch strength were studied using a hand dynamometer and the best of 3 attempts was used for analysis. Daily life activities were assessed using the ACTIVity LIMitations scale for patients with upper and/or lower limb impairments (ACTIVLIM) (20). We also assessed respiratory function using the Carefusion Microlab ML 3500 MK8 spirometer (Carefusion, Yorba Linda, CA, USA). We analyzed the Forced Vital Capacity (FVC) while seating and lying and the Forced Expiratory Volume in the first second (FEV1).
Quantification of Growth Factors in Serum Samples
We measured the serum concentration of the following growth factors related with muscle fibrosis: Serum Platelet-Derived Growth Factor BB (PDGF-BB), Transforming Growth Factor β1 (TGF-β1), Platelet-Derived Growth Factor AA (PDGF-AA) and connective tissue growth factor (CTGF). Blood was centrifuged at 1,600 g for 9 min at 4°C. Samples were aliquoted and stored at −80°C until analysis. PDGF-BB and TGF-β1 levels were measured using commercial enzyme-linked immunosorbent assay (ELISA) kits (R&D, Minneapolis, MN, USA), according to the manufacturer's instructions. PDGF-AA human ELISA kit was provided by ThermoFisher (Thermo Fisher Scientific, Nepean, Canada) and CTGF by EIAAB Science Co (Wuhan, China). Minimum detectable cytokine concentrations for these assays were measured to be 1.7 pg/ml for TGF-β1, 15 pg/ml for PDGF-BB, 40 pg/ml for PDGF-AA and 0.18 ng/ml for CTGF. Samples were measured in duplicate and read on a microplate reader Beckman Coulter AD 340 (Beckam-Coulter, Brea, CA, USA) with AD-LD software.
These growth factors were determined in the 19 DMD patients, 13 BMD patients, 15 patients with mutations in the DYSF gene that were being visited in HSCSP (from now on DYSF patients) and in 66 healthy controls. The healthy controls were divided into a pediatric group (<15 years old) and adult group (≥15 years old). Since the average age of DMD cohort is significantly different from the DMB and DYSF cohort, the DMD group was compared to the pediatric control group while the DMB and DYSF group was compared to the adult control group.
7 DMD and 8 BMD patients underwent thigh muscle MRI in a 1.5T Ingenia MR system (Philips Healthcare, Best, the Netherlands) at HSCSP. Axial 3D FFE Dixon sequence was acquired with the following parameters: TR/TE = 5.78/1.84 ms, flip angle = 15°, voxel size = 1 × 1 × 3 mm and FOV 520 × 340 × 300 mm. Water and fat images were automatically obtained from the Dixon acquisition using a single peak model. One investigator (A.A-J) estimated fat content in the images using the PRIDE (Philips Research Image Development Environment) tool developed for this purpose. This tool provides the fat fraction (FF), defined as fat/(fat+water). Regions of interest (ROIs) were manually drawn on five slices on the right thigh of the following muscles: rectus femoris, vastus intermedius, vastus lateralis, vastus medialis, adductor magnus, sartorius, gracilis, biceps femoris long head, semitendinosus, and semimembranosus. We intended to include slices of the proximal, middle and distal part of the thigh. To do so, we identified the head of the femur and set the first slice at 5 centimeters below the head. The remaining four slices were distally distributed in the thighs with a separation of 20 slices between them as previously published (21). In patients in which individual muscles in the posterior or anterior compartment of the thigh could not be analyzed due to a complete fat replacement, a ROI was manually drawn over the whole compartment. Neurovascular bundles were avoided. An average of the fat fraction of the 5 slices was calculated for each muscle, and weighted averages were calculated by normalizing the FF by the area of the muscle (bigger muscles had higher weight) to obtain the global thigh FF.
RNA Extraction and Reverse Transcription
Total RNA was extracted from triceps or biceps biopsy samples from 2 healthy controls and 2 DMD patients using RNeasy® Micro Kit (Qiagen, Hilden, Germany) following manufacturer's instruction. RNA was quantified using a nanodrop ND-1000 spectrophotometer (Nanodrop Technologies Inc., Wil- mington, DE, USA). One μg of total RNA was reverse-transcribed to complementary DNA (cDNA) using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA).
Real-Time Quantitative PCR Analysis
Real-Time PCR (qPCR) was performed using the TaqMan® Universal PCR Master Mix (Applied Biosystems, Foster City, CA, USA) and a 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA). All mRNA- specific FAM-labeled primers/probe were purchased from Applied Biosystems. Relative quantification was performed using the comparative Ct method and all results were compared with the control samples. GAPDH was used as endogenous control.
PDGF-AA immunohistochemistry was performed from paraffinized muscle. The samples were deparaffinized (xylol, absolute ethanol, 95° ethanol, and 70° ethanol), placed in distilled water and pretreated with 10 mM citrate (pH 6) at 100°C. Then samples were washed in PBS and incubated in blocking solution (4% bovine serum albumin). The primary anti-PDGF-AA antibody (Merck Millipore, Darmstadt, Germany) was then added at room temperature for 1 h. After 3 washes, it was incubated with biotinylated anti-rabbit antibody (Vector Laboratory Inc. Burlingame, CA) for 1 h at room temperature and then avidin-biotin peroxidase complex (Dako, Glostrup, Denmark) was added. Finally, the sections were soaked in Mayer's hematoxylin for 10 s, washed under running water and mounted in the aqueous mounting medium Aquatex (Merck Millipore, Darmstadt, Germany).
Non-parametric tests were used for the statistical analysis as we demonstrated that the distribution of the variables was not uniform as studied using Kolmogorov-Smirnov and did not follow a normal distribution as studied using Shapiro-Willis test. The Mann-Whitney U test was used to identify significant differences in variables between 2 groups. Non-parametric Kruskal–Wallis test analysis was used to assess differences among 3 groups followed by Dunn's multiple comparison post-test. We used Spearman's rank correlation (coefficient reported as ρ) to investigate any correlation between the serum concentration of growth factors and the results of the muscle function tests, quality of life scales and the thigh FF obtained using Dixon imaging. Correlation coefficients higher than 0.6 were considered good. Because the objective of our research was to explore the possibilities of these growth factors for further studies, we decided not to use multiple comparisons. The results of all statistical studies were considered significant if P-value was lower than 0.05. Statistical studies were performed using IBM SPSS® Statistics software version 21. Graphs were created using Graphpad software and R software.
A total of 19 patients with DMD (mean age 10 ± 5.6 years), 13 patients with BMD (mean age 31 ± 20.6 years), 15 DYSF patients (mean age 46 ± 18.6 years), 8 pediatric controls (mean age 10 ± 2 years) and 58 adult controls (mean age 43 ± 15.32 years) were enrolled in the study. There was no statistical difference in age between the DMD group and the pediatric controls (Mann-Whitney U-test, p = 0.979) and among the BMD, DYSF group and the adult control group (Kruskall-Wallis test, p = 0.052) We identified deletions in the DMD gene in 71.9% of patients. Duplications accounted for 12.5% of cases and 15.6% were due to point mutations. Among the DMD patients, 13 were ambulant (68.42%) and 6 non-ambulant (31.58%). In the BMD group, 8 patients were ambulant (61.5%) and 5 non-ambulant (38.5%). Five DYSF patients were non-ambulant (33.33%). Table 1 describes the mean demographic and clinical data of all patients included in the cohort.
Growth Factor Serum Levels Comparison Between DMD, BMD, and Age Matched Controls
Our first aim was to study whether there were differences in growth factor serum levels between DMD, BMD and DYSF patients compared to controls. We compared serum levels of PDGF-AA, PDGF-BB, CTGF, and TGF-β1 between DMD and pediatric healthy controls whereas BMD and DYSF patients were compared against adult healthy controls (Figure 1). In the case of DMD patients, we only observed significant differences in the serum levels of PDGF-AA as shown in Figure 1A. The median and IQR of PDGF-AA serum levels were 2245.32 ng/ml (IQR: 1589.94-3832.68) in DMD patients compared to 1242.45 ng/ml (IQR: 915.05-1373.13) in pediatric controls (Mann-Whitney U-test, p = 0.002). In the case of BMD and DYSF patients we observed significant differences in the serum levels of CTGF between BMD and adult controls and of TGF-β1 between BMD and DYSF patients as shown in Figures 1C,D. The median and IQR of CTGF serum levels were 5.66 ng/ml (IQR: 3.82-8.82) in BMD patients compared to 2.37 ng/ml (IQR: 1.35-4.07) in adult controls (Kruskal-Wallis test/Dunn's multiple comparison analysis, p = 0.004). The median and IQR of TGF-β1 serum levels were 59.10 ng/ml (IQR: 51.11-71.51) in BMD patients compared to 40.42 ng/ml (IQR: 28.44-53.75) in DYSF patients (Kruskal-Wallis test/Dunn's multiple comparison analysis, p = 0.01). Although PDGF-AA did not show statistical differences, there was a tendency toward increased levels in BMB patients 2415.60 ng/ml (IQR: 1228.59-3981.10) compared to adult controls 1859.62 (IQR: 1363.08-2477.52) (Kruskal-Wallis test/Dunn's multiple comparison analysis, p = 0.06).
Figure 1. Serum levels of different growth factors in DMD, BMD, DYSF and in pediatric and adult healthy controls. (A) PDGF-AA, (B) PDGF-BB, (C) TGF-β1, and (D) CTGF levels were measured. Variables are represented as median and interquartile range (IQR). Statistical significance of the results by Mann-Whitney test for the DMD compared against pediatric controls and Kruskal-Wallis test/Dunn's multiple comparison analysis for the BMD, DYSF, and adult controls: *p ≤ 0.05 and **p ≤ 0.01.
Correlation of Growth Factor Serum Samples With Muscle Function Test Results
In a second step, we analyzed if there was a correlation between growth factor serum concentration and the results of muscle function tests, spirometry, daily life activities scales and Dixon MRI in DMD and BMD patients. The growth factor that showed a larger number of significant correlations was PDGF-AA as it is shown in Table 2. In the case of DMD, PDGF-AA serum levels correlated with the results of the 10MWT, the 6MWT, the MRC scale. and the dimension 3 of the MFM scale. In the case of BMD, PDGF-AA correlated with the results of the 10MWT, the 6MWT, the Tup4, Tdo4, TRF and the NSAA, FVC sitting, FEV1 and the thigh FF. Figure 2 shows graphs plotting these correlations. Moreover, we analyzed the distribution of PDGF-AA values in relation to age and to the MRC total score and we observed a bell distribution both in DMD and in BMD patients (Figures 2E,F). Younger and stronger patients had low levels of PDGF-AA; patients in an intermediate stage that were losing muscle strength had higher levels of PDGF-AA, while patients in the end-stage of the disease, that were very weak or older, had again low levels of PDGF-AA. Muscle strength measured by MRC significantly correlated with TGF- β1 levels in DMD patients (p = 0.026, Spearman's rho = 0.523) and with CTGF levels in BMD patients (p = 0.012, Spearman's rho = −0.669). We did not find statistically significant correlations between the PDGF-BB serum levels and the function tests.
Table 2. Correlation between PDGF-AA serum levels and different muscle function tests, spirometry values, and muscle MRI.
Figure 2. Correlation between PDGF-AA levels and different muscle tests in patients with DMD and BMD. (A) Correlation between PDGF-AA and the results of the 6MWT; (B) Correlation between PDGF-AA 10 MWT; (C) Correlation between PDGF-AA and results of the NSAA scale; (D) Correlation between PDGF-AA and thigh fat fraction; (E) Distribution of PDGF-AA serum levels depending on the muscle strength quantified using MRC score; (F) distribution of PDGF-AA serum levels depending on the age of the patients. Spearman test was performed, R, Correlation coefficient; p, statistical significance. Red lines: BMD patients; blue lines: DMD patients.
PDGF-AA Expression in Skeletal Muscle Biopsies of DMD and BMD Patients
In a final step we studied PDGF-AA expression patterns in muscle biopsies of patients with DMD and BMD and in controls using immunohistochemistry and Real Time PCR. As it has been previously shown, PDGF-AA is expressed in the capillaries in muscles of healthy controls, where it probably has a role in binding the pericytes to the endothelial cell (Figure 3A) (22, 23). In patients with DMD, we observed that muscle fibers expressed PDGF-AA as it is shown in Figure 3B. qPCR confirmed the results of immunohistochemistry and showed an increased expression of PDGF-AA in muscles from DMD compared to healthy controls (Figure 3C).
Figure 3. PDGF-AA expression is increased in muscles of patients with dystrophinopathy. (A) PDGF-AA expression in skeletal muscles of healthy controls; (B) PDGF-AA expression in a muscle biopsy of a DMD patient and (C) Real Time PCR showing PDGF-AA mRNA expression in skeletal muscle samples of controls and patients with dystrophinopathy. Bars show mean mRNA expression and standard deviation.
The goal of the work behind this manuscript was to identify serum candidate biomarkers to track progression in DMD and BMD.
The current development of new therapies for muscle dystrophies, partly due to the progress in the field of genetics in the latest decades, raises the need for reliable outcome measures. DMD and BMD are slow progressive diseases where weakness is developed over years. However, clinical trials designed for these disorders must be able to prove presence or absence of effect in a shorter period of time. Biomarkers that can closely track the pathophysiological mechanisms of these diseases, and therefore be used as outcome measures, would be desirable. The research in the field has been centered in two different measurements. On the one hand, qMRI has demonstrated to be a reliable tool to analyze changes in muscle structure in patients with muscular diseases. Fat replacement, assessed using Dixon imaging or spectroscopy could be a suitable outcome measure, because it correlates with muscle function and it is sensitive to changes over short periods in several disorders such as DMD, Pompe or LGMD-2I (7, 15, 24). On the other hand, blood biomarkers are especially interesting because they are easy to obtain. Serum concentration of molecules such as cytokines, metalloproteinases or microRNAs have shown promising results in some neuromuscular diseases (16, 25). Creatine kinase (CK) is a widely accepted biomarker for diagnosis of DMD, BMD and other muscular dystrophies. However, its levels fluctuate and are influenced by different conditions such as physical activity, therefore it is not considered a good biomarker to study diseases' progression (16). Other molecules present in muscle and involved in the pathophysiology of muscular dystrophies have been studied. For example, Hathout et al. (26) found 44 proteins which were able to differentiate DMD patients from healthy volunteers. Other studies aimed not only to differentiate between patients and controls, but also to correlate with functional state and treatment (27–29) or muscle function tests (30). Fewer manuscripts aim to provide longitudinal data on the suggested biomarkers (31–33).
We were interested in investigating the serum concentration of PDGF-AA, PDGF-BB, CTGF, and TGF β1 in our patients because of their role in fibrosis. The process of muscle fibrosis in patients with muscular dystrophies is highly complex and not completely understood (34). It is known that chronic muscle degeneration leads to persistent inflammatory infiltration, muscle necrosis and activation of fibroadipogenic progenitor cells (FAPs) and fibroblasts, which release proteins of the extracellular matrix such as collagen-I, leading to the expansion of fibrotic tissue (34–36). The PDGF family in general and PDGF-AA in particular play an important role in this process (22, 37). PDGF-AA is produced by different cell types, such as platelets, muscle fibers or macrophages and acts via the PDGF receptor alpha (PDGFRα) (38). When PDGF-AA binds to its receptor, several signaling pathways involved with cellular proliferation, differentiation and motility are triggered (39). After an acute muscle injury, damaged fibers and inflammatory cells release several cytokines, including PDGF-AA, which activate FAPs inducing the release of extracellular matrix components (40). Extracellular matrix is needed for muscle regeneration as it acts as a scaffold for muscle fiber regeneration (41). When muscle regeneration is completed, inflammatory cells and the new muscle fibers release TNF-α that promotes apoptosis of FAPs and fibroblasts reducing their number and avoiding the development of fibrosis. At that stage, levels of PDGF-AA decrease (42). However, in dystrophinopathies where there is chronic muscle damage, inflammatory cells permanently infiltrate the injured muscle continuously releasing PDGF-AA, which expand the number of activated FAPs leading to an increased non-controlled fibrotic response (40, 43). Consequently, treatment of murine models of DMD with tyrosine kinase inhibitors that block the PDGFRα such as imatinib, crenolanib and nintedanib, decreased fibrosis and improved muscle function of these animals (40, 42, 44). We think that due to this increased pro-fibrotic response in patients with DMD, the levels of PDGF-AA are high in comparison to controls, as shown in our study.
Other primary fibrotic growth factors including the TGF-β1 family and CTGF have also been reported to be involved in the process of muscle fibrosis (36, 45, 46). TGF-β1 and PDGF-BB participate in mesenchymal cell proliferation and fibrotic remodeling (47, 48) and CTGF stimulates chemotaxis and proliferation of fibroblasts and the production of extracellular matrix components (45).
Despite its role in the physiopathology of muscle degeneration, there are only few reports studying the potential utility of growth factors as biomarkers of the disease (49–51).
In our study, the serum concentration of PDGF-AA was increased in patients with DMD in comparison to pediatric controls. There was also an increase in PDGF-AA levels in BMD patients compared to adult controls, although it did not reach statistical significance. Moreover, PDGF-AA serum levels correlated with several muscle function tests both in DMD and BMD, and with the results of Dixon MRI in BMD. Because of the common underlying pathophysiologic mechanism of both DMD and BMD, we speculate that serum levels of PDGF-AA should be further studied in a larger cohort of DMD/BMD patients and considered as a potential biomarker for the disease.
Our results suggest that the expression of PDGF-AA in DMD and BMD patients follows a bell curve distribution. In initial stages of the disease, when histological changes are consistent with acute muscle damage, PDGF-AA expression by muscle fibers is still low, similarly to that observed in controls. As the disease progresses, muscle fibers degenerate leading to increased production of PDGF-AA that potentially activates FAPs involved in the expansion of fibrotic and adipogenic tissue (22). In advanced stages, when most of the muscle has been replaced by fibrotic and fatty tissue, levels drop again as there are no more muscle fibers releasing PDGF-AA. This hypothesis allows us to understand the results obtained in the correlations between PDGF-AA levels and the muscular function of the patients. Patients in early stages of the disease, in whom muscle strength assessed using the MRC scale is only mildly affected, have low levels of PDGF-AA. These patients in early stages are still able to walk, have high scores on the NSAA, need little time to walk 10 m and cover long distances in the 6MWT. Patients in an intermediate stage of the disease have high levels of PDGF-AA and are clearly weak as it is demonstrated using the MRC scales. These patients are able to cover 200–400 m in the 6MWT, but they still have intact respiratory muscle strength. In contrast, patients in advanced stages of the disease have lower levels of serum PDGF-AA, are extremely weak as assessed with the MRC scale and have higher fat replacement in the muscles. Moreover, these patients are not able to walk and have affected respiratory muscles as shown using the spirometry. Although levels of CTGF were also increased in BMD patients in comparison to controls, we only found significant correlations with the MRC scale.
The process of fibrosis is not exclusive to muscular tissue. In fact, fibrotic tissue remodeling can affect virtually every system. Both primary fibrotic diseases such as systemic sclerosis and idiopathic pulmonary fibrosis and secondary fibrotic responses in other organs share many commonalities. An initial injury triggers the reparative process in which lymphocytes and macrophages release profibrotic mediators. These mediators promote the activation of myofibroblasts which liberate extracellular matrix proteins leading to structural changes. Therefore, it is not surprising to find that the growth factors that we have studied here are also implicated in other diseases such as atherosclerosis, liver and pulmonary fibrosis, cirrhosis, myocardial infarction and neoplasms (52–58). However, the primary site of injury, the cells that differentiate into myofibroblasts and the organization of the immune response differs among the different types of disorders (52). Interestingly, in our study TGF-β1 showed a different behavior in BMD than in DYSF patients, since TGF-β1 levels were increased in BMD patients compared to DYSF patients. In our opinion this data suggests that, although the final stage of the disease may be common among muscular dystrophies, the process leading to replacement of muscle fibers by fibrotic and fatty tissue may have some particularities that could be exclusive to each type of muscular dystrophy, such as the type of muscle injury or the molecular pathways that are activated in response to this muscle injury. However, we did not find significant differences in TGF- β1 between patients and controls (only between BMD and DYSF patients) and therefore these results should be taken with caution.
The four growth factors that we have studied have been previously related with fibrosis and with muscle regeneration in muscular dystrophies (59, 60) but have also been studied in other neuromuscular diseases. For instance, our group identified a decrease in PDGF-BB levels in patients with Pompe disease suggesting that there is a deficit in muscle regeneration in this disease, and Sato et al. showed that CTGF levels were high in patients with polymyositis and dermatomyositis (61, 62).
This is a small study with several limitations: firstly, it is a transversal study without longitudinal data. Secondly, only half of the patients underwent Dixon MRI, and the ones who did were already severely affected by the disease and therefore had high FF.
Despite these limitations, our study suggests that serum levels of PDGF-AA can differentiate between DMD patients and controls, and that they correlate with the results of different muscle tests. For BMD patients we have also found correlation between PDGF-AA levels and muscle function tests, including fat muscle content assessed by Dixon MRI. Based on these results, we propose PDGF-AA as a growth factor to be further investigated in in a larger cohort of patients and in longitudinal studies.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
The studies involving human participants were reviewed and approved by Hospital de la Santa Creu i Sant Pau Ethics Committee and Hospital de Sant Joan de Deu Ethics Committee. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
AA-J, EF-S, and JD-M contributed to the study conception and design. Clinical evaluation of patients was performed by AA-J, JD-M, JA-P, DN-d, CO, and CJ-M. Blood samples and data collection was obtained by SS. EG, JD-M, and EF-S designed and performed the experiments. PP-J, AC-R, and XS-C also collaborated with the execution of the experiments. PM, JD-M, CN-P, and JL designed the MRI protocol. The analysis of the MRI images was performed by AA-J. Functional assessments were performed by IB, IP, CG, and EM. AA-J, EF-S, and JD-M analyzed the results and wrote the manuscript. All authors contributed with their comments to the final version.
This work was supported by FIS PI15/01822 to JD-M and FIS (PI18/1525) to JD-M and XS-C, and FIS (PI19/00122) to CJ-M (Fondo Europeo de Desarrollo Regional [FEDER], Instituto de Salud Carlos III, Spain). XS-C was supported by Sara Borrell postdoctoral fellowship (CD18/00195), Fondo Social Europeo (FSE), Instituto de Salud Carlos III (Spain). Funders did not participate in the study design. This work is generated within the European Reference Network for Neuromuscular diseases.
Conflict of Interest
PM is employed by the company Philips Healthcare Iberia.
The remaining 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.
We would like to thank all the patients and families who participated in the study for their support and patience. We are indebted to the Biobanc de l'Hospital Infantil Sant Joan de Déu per a la Investigació integrated in the Spanish Biobank Network of ISCIII for the sample and data procurement. We thank the MRI technician team for their help. We also thank Mr. John Wilkos for his language advice.
1. Birnkrant DJ, Bushby K, Bann CM, Apkon SD, Blackwell A, Brumbaugh D, et al. Diagnosis and management of Duchenne muscular dystrophy, part 1: diagnosis, and neuromuscular, rehabilitation, endocrine, and gastrointestinal and nutritional management. Lancet Neurol. (2018) 17:211–2. doi: 10.1016/S1474-4422(18)30024-3
2. Wallace GQ, McNally EM. Mechanisms of muscle degeneration, regeneration, and repair in the muscular dystrophies. Annu Rev Physiol. (2009) 71:37–57. doi: 10.1146/annurev.physiol.010908.163216
3. Guiraud S, Aartsma-Rus A, Vieira NM, Davies KE, Van Ommen GJB, Kunkel LM. The pathogenesis and therapy of muscular dystrophies. Annu Rev Genomics Hum Genet. (2015) 16:281–308. doi: 10.1146/annurev-genom-090314-025003
4. Flanigan KM. Duchenne and becker muscular dystrophies. Neurol Clin. (2014) 32:671–88. doi: 10.1016/j.ncl.2014.05.002
5. Domingos J, Sarkozy A, Scoto M, Muntoni F. Dystrophinopathies and limb-girdle muscular dystrophies. Neuropediatrics. (2017) 48:262–72. doi: 10.1055/s-0037-1601860
6. Dowling JJD, Gonorazky H, Cohn RD, Campbell C. Treating pediatric neuromuscular disorders: The future is now. Am J Med Genet Part A. (2018) 176:804–41. doi: 10.1002/ajmg.a.38418
7. Barp A, Bello L, Caumo L, Campadello P, Semplicini C, Lazzarotto A, et al. Muscle MRI and functional outcome measures in Becker muscular dystrophy. Sci Rep. (2017) 7:16060. doi: 10.1038/s41598-017-16170-2
8. Arora H, Willcocks RJ, Lott DJ, Harrington AT, Senesac CR, Zilke KL, et al. Longitudinal timed function tests in duchenne muscular dystrophy: Imagingdmd cohort natural history. Muscle Nerve. (2018) 58:631–8. doi: 10.1002/mus.26161
9. Bello L, Campadello P, Barp A, Fanin M, Semplicini C, Sorarù G, et al. Functional changes in Becker muscular dystrophy: Implications for clinical trials in dystrophinopathies. Sci Rep. (2016) 6:32439. doi: 10.1038/srep32439
10. Mazzone ES, Messina S, Vasco G, Main M, Eagle M, D'Amico A, et al. Reliability of the north star ambulatory assessment in a multicentric setting. Neuromuscul Disord. (2009) 19:458–61. doi: 10.1016/j.nmd.2009.06.368
11. Schmidt S, Hafner P, Klein A, Rubino-Nacht D, Gocheva V, Schroeder J, et al. Timed function tests, motor function measure, and quantitative thigh muscle MRI in ambulant children with Duchenne muscular dystrophy. A cross-sectional analysis. Neuromuscul Disord. (2018) 28:16–23. doi: 10.1016/j.nmd.2017.10.003
12. Long C, Amoasii L, Mireault AA, McAnally JR, Li H, Sanchez-Ortiz E, et al. Postnatal genome editing partially restores dystrophin expression in a mouse model of muscular dystrophy. Science. (2016) 351:400–3. doi: 10.1126/science.aad5725
13. Mayhew A, Mazzone ES, Eagle M, Duong T, Ash M, Decostre V, et al. Development of the performance of the upper limb module for duchenne muscular dystrophy. Dev Med Child Neurol. (2013) 55:1038–45. doi: 10.1111/dmcn.12213
14. Godi C, Ambrosi A, Nicastro F, Previtali SC, Santarosa C, Napolitano S, et al. Longitudinal MRI quantification of muscle degeneration in Duchenne muscular dystrophy. Ann Clin Transl Neurol. (2016) 3:607–22. doi: 10.1002/acn3.319
15. Damon BM, Li K, Dortch RD, Welch EB, Park JH, Buck AKW, et al. Quantitative magnetic resonance imaging of skeletal muscle disease. J Vis Exp. (2016) 2016:52352. doi: 10.3791/52352
16. Al-Khalili Szigyarto C, Spitali P. Biomarkers of Duchenne muscular dystrophy: current findings. Degener Neurol Neuromuscul Dis. (2018) 8:1–13. doi: 10.2147/dnnd.s121099
17. Bérard C, Payan C, Hodgkinson I, Fermanian J. A motor function measure scale for neuromuscular diseases. Construction and validation study. Neuromuscul Disord. (2005) 15:463–70. doi: 10.1016/j.nmd.2005.03.004
18. Fagoaga J, Girabent-Farrés M, Bagur-Calafat C, Febrer A, Steffensen BF. Translation and validation of the Egen Klassifikation scale for the Spanish population: functional assessment for non-ambulatory individuals with Duchenne's muscular dystrophy and spinal muscular atrophy. Rev Neurol. (2013) 56:555–61. doi: 10.33588/rn.5611.2013106
19. Harris E, Bladen CL, Mayhew A, James M, Bettinson K, Moore U, et al. The clinical outcome study for dysferlinopathy an international multicenter study. Neurol Genet. (2016) 2:89. doi: 10.1212/NXG.0000000000000089
20. Bleyenheuft Y, Paradis J, Renders A, Thonnard JL, Arnould C. ACTIVLIM-CP a new Rasch-built measure of global activity performance for children with cerebral palsy. Res Dev Disabil. (2017) 60:285–94. doi: 10.1016/j.ridd.2016.10.005
21. Figueroa-Bonaparte S, Llauger J, Segovia S, Belmonte I, Pedrosa I, Montiel E, et al. Quantitative muscle MRI to follow up late onset Pompe patients: a prospective study. Sci Rep. (2018) 8:10898. doi: 10.1038/s41598-018-29170-7
22. Zhao Y, Haginoya K, Sun G, Dai H, Onuma A, Iinuma K. Platelet-derived growth factor and its receptors are related to the progression of human muscular dystrophy: An immunohistochemical study. J Pathol. (2003) 201:149–59. doi: 10.1002/path.1414
23. Andrae J, Gallini R, Betsholtz C. Role of platelet-derived growth factors in physiology and medicine. Genes Dev. (2008) 22:1276–312. doi: 10.1101/gad.1653708
24. Diaz-Manera J, Fernandez-Torron R, LLauger J, James MK, Mayhew A, Smith FE, et al. Muscle MRI in patients with dysferlinopathy: pattern recognition and implications for clinical trials. J Neurol Neurosurg Psychiatry. (2018) 89:1071–81. doi: 10.1136/jnnp-2017-317488
25. Strandberg K, Ayoglu B, Roos A, Reza M, Niks E, Signorelli M, et al. Blood-derived biomarkers correlate with clinical progression in Duchenne muscular dystrophy. J Neuromuscul Dis. (2020) 7:231–46. doi: 10.3233/JND-190454
26. Hathout Y, Brody E, Clemens PR, Cripe L, DeLisle RK, Furlong P, et al. Large-scale serum protein biomarker discovery in Duchenne muscular dystrophy. Proc Natl Acad Sci USA. (2015) 112:7153–8. doi: 10.1073/pnas.1507719112
27. Burch PM, Pogoryelova O, Goldstein R, Bennett D, Guglieri M, Straub V, et al. Muscle-derived proteins as serum biomarkers for monitoring disease progression in three forms of muscular dystrophy. J Neuromuscul Dis. (2015) 2:241–55. doi: 10.3233/JND-140066
28. Ayoglu B, Chaouch A, Lochmüller H, Politano L, Bertini E, Spitali P, et al. Affinity proteomics within rare diseases: a BIO - NMD study for blood biomarkers of muscular dystrophies. EMBO Mol Med. (2014) 6:918–36. doi: 10.15252/emmm.201303724
29. Burch PM, Pogoryelova O, Palandra J, Goldstein R, Bennett D, Fitz L, et al. Reduced serum myostatin concentrations associated with genetic muscle disease progression. J Neurol. (2017) 264:541–53. doi: 10.1007/s00415-016-8379-6
30. Anaya-Segura MA, García-Martínez FA, Montes-Almanza LÁ, Díaz BG, Ávila-Ramírez G, Alvarez-Maya I, et al. Non-invasive biomarkers for duchenne muscular dystrophy and carrier detection. Molecules. (2015) 20:11154–72. doi: 10.3390/molecules200611154
31. Signorelli M, Ayoglu B, Johansson C, Lochmüller H, Straub V, Muntoni F, et al. Longitudinal serum biomarker screening identifies malate dehydrogenase 2 as candidate prognostic biomarker for Duchenne muscular dystrophy. J Cachexia Sarcopenia Muscle. (2020) 11:505–17. doi: 10.1002/jcsm.12517
32. Nadarajah VD, van Putten M, Chaouch A, Garrood P, Straub V, Lochmüller H, et al. Serum matrix metalloproteinase-9 (MMP-9) as a biomarker for monitoring disease progression in Duchenne muscular dystrophy (DMD). Neuromuscul Disord. (2011) 21:569–78. doi: 10.1016/j.nmd.2011.05.011
33. Spitali P, Hettne K, Tsonaka R, Charrout M, van den Bergen J, Koeks Z, et al. Tracking disease progression non-invasively in Duchenne and Becker muscular dystrophies. J Cachexia Sarcopenia Muscle. (2018) 9:715–26. doi: 10.1002/jcsm.12304
34. Bersini S, Gilardi M, Mora M, Krol S, Arrigoni C, Candrian C, et al. Tackling muscle fibrosis: From molecular mechanisms to next generation engineered models to predict drug delivery. Adv Drug Deliv Rev. (2018) 129:64–77. doi: 10.1016/j.addr.2018.02.009
35. Cappellari O, Mantuano P, De Luca A. “The Social Network” and muscular dystrophies: the lesson learnt about the niche environment as a target for therapeutic strategies. Cells. (2020) 9:1659. doi: 10.3390/cells9071659
36. Smith LR, Barton ER. Regulation of fibrosis in muscular dystrophy. Matrix Biol. (2018) 68–9:602–15. doi: 10.1016/j.matbio.2018.01.014
37. Piñol-Jurado P, Gallardo E, de Luna N, Suárez-Calvet X, Sánchez-Riera C, Fernández-Simón E, et al. Platelet-derived growth factor BB influences muscle regeneration in duchenne muscle dystrophy. Am J Pathol. (2017) 187:1814–27. doi: 10.1016/j.ajpath.2017.04.011
38. Fredriksson L, Li H, Eriksson U. The PDGF family: Four gene products form five dimeric isoforms. Cytokine Growth Factor Rev. (2004) 15:197–204. doi: 10.1016/j.cytogfr.2004.03.007
39. Betsholtz C. Insight into the physiological functions of PDGF through genetic studies in mice. Cytokine Growth Factor Rev. (2004) 15:215–28. doi: 10.1016/j.cytogfr.2004.03.005
40. Ieronimakis N, Hays A, Prasad A, Janebodin K, Duffield JS, Reyes M. PDGFRα signalling promotes fibrogenic responses in collagen-producing cells in Duchenne muscular dystrophy. J Pathol. (2016) 240:410–44. doi: 10.1002/path.4801
41. Baghdadi MB, Tajbakhsh S. Regulation and phylogeny of skeletal muscle regeneration. Dev Biol. (2018) 433:200–9. doi: 10.1016/j.ydbio.2017.07.026
42. Lemos DR, Babaeijandaghi F, Low M, Chang CK, Lee ST, Fiore D, et al. Nilotinib reduces muscle fibrosis in chronic muscle injury by promoting TNF-mediated apoptosis of fibro/adipogenic progenitors. Nat Med. (2015) 21:786–94. doi: 10.1038/nm.3869
43. Olson LE, Soriano P. Increased PDGFRα activation disrupts connective tissue development and drives systemic fibrosis. Dev Cell. (2009) 16:303–13. doi: 10.1016/j.devcel.2008.12.003
44. Piñol-Jurado P, Suárez-Calvet X, Fernández-Simón E, Gallardo E, De La Oliva N, Martínez-Muriana A, et al. Nintedanib decreases muscle fibrosis and improves muscle function in a murine model of dystrophinopathy. Cell Death Dis. (2018) 9:776. doi: 10.1038/s41419-018-0792-6
45. Moussad EEDA, Brigstock DR. Connective tissue growth factor: what's in a name? Mol Genet Metab. (2000) 71:276–92. doi: 10.1006/mgme.2000.3059
46. Klinkhammer BM, Floege J, Boor P. PDGF in organ fibrosis. Mol Aspects Med. (2018) 62:44–62. doi: 10.1016/j.mam.2017.11.008
47. Yablonka-Reuveni Z, Balestreri TM, Bowen-Pope DF. Regulation of proliferation and differentiation of myoblasts derived from adult mouse skeletal muscle by specific isoforms of PDGF. J Cell Biol. (1990) 111:1623–9. doi: 10.1083/jcb.111.4.1623
48. Delaney K, Kasprzycka P, Ciemerych MA, Zimowska M. The role of TGF-β1 during skeletal muscle regeneration. Cell Biol Int. (2017) 41:706–15. doi: 10.1002/cbin.10725
49. Tsonaka R, Signorelli M, Sabir E, Seyer A, Hettne K, Aartsma-Rus A, et al. Longitudinal metabolomic analysis of plasma enables modeling disease progression in Duchenne muscular dystrophy mouse models. Hum Mol Genet. (2020) 29:745–55. doi: 10.1093/hmg/ddz309
50. Morales MG, Acuña MJ, Cabrera D, Goldschmeding R, Brandan E. The pro-fibrotic connective tissue growth factor (CTGF/CCN2) correlates with the number of necrotic-regenerative foci in dystrophic muscle. J Cell Commun Signal. (2018) 12:413–21. doi: 10.1007/s12079-017-0409-3
51. Song Y, Yao S, Liu Y, Long L, Yang H, Li Q, et al. Expression levels of TGF-β1 and CTGF are associated with the severity of duchenne muscular dystrophy. Exp Ther Med. (2017) 13:1209–14. doi: 10.3892/etm.2017.4105
52. Distler JHW, Györfi AH, Ramanujam M, Whitfield ML, Königshoff M, Lafyatis R. Shared and distinct mechanisms of fibrosis. Nat Rev Rheumatol. (2019) 15:705–30. doi: 10.1038/s41584-019-0322-7
53. Bujak M, Frangogiannis NG. The role of TGF-β signaling in myocardial infarction and cardiac remodeling. Cardiovasc Res. (2007) 74:184–95. doi: 10.1016/j.cardiores.2006.10.002
54. Kishi M, Aono Y, Sato S, Koyama K, Azuma M, Abe S, et al. Blockade of platelet-derived growth factor receptor-β, not receptor-α ameliorates bleomycin-induced pulmonary fibrosis in mice. PLoS ONE. (2018) 13:e0209786. doi: 10.1371/journal.pone.0209786
55. Roskoski R. The role of small molecule platelet-derived growth factor receptor (PDGFR) inhibitors in the treatment of neoplastic disorders. Pharmacol Res. (2018) 129:65–83. doi: 10.1016/j.phrs.2018.01.021
56. Shamsdin SA, Mehrafshan A, Rakei SM, Mehrabanilare D. Evaluation of VEGF, FGF and PDGF and serum levels of inflammatory cytokines in patients with glioma and meningioma in Southern Iran. Asian Pacific J Cancer Prev. (2019) 20:2883–90. doi: 10.31557/APJCP.2019.20.10.2883
57. Ulloa-Padilla JP, Ghassibi MP, Dubovy SR, Kerr DA. Clinicopathologic correlation of kaposi sarcoma involving the ocular adnexa: immunophenotyping of diagnostic and therapeutic targets. Ophthal Plast Reconstr Surg. (2020) 36:185–90. doi: 10.1097/IOP.0000000000001506
58. Kurys-Denis E, Prystupa A, Luchowska-Kocot D, Krupski W, Bis-Wencel H, Panasiuk L. PDGF-BB homodimer serum level – A good indicator of the severity of alcoholic liver cirrhosis. Ann Agric Environ Med. (2020) 27:80–5. doi: 10.26444/aaem/115997
59. Morales MG, Gutierrez J, Cabello-Verrugio C, Cabrera D, Lipson KE, Goldschmeding R, et al. Reducing CTGF/CCN2 slows down mdx muscle dystrophy and improves cell therapy. Hum Mol Genet. (2013) 22:4938–51. doi: 10.1093/hmg/ddt352
60. Ceco E, Bogdanovich S, Gardner B, Miller T, DeJesus A, Earley JU, et al. Targeting latent TGFβ release in muscular dystrophy. Sci Transl Med. (2014) 6:259ra144. doi: 10.1126/scitranslmed.3010018
61. Fernández-Simón E, Carrasco-Rozas A, Gallardo E, Figueroa-Bonaparte S, Belmonte I, Pedrosa I, et al. PDGF-BB serum levels are decreased in adult onset Pompe patients. Sci Rep. (2019) 9:1–10. doi: 10.1038/s41598-018-38025-0
62. Sato S, Nagaoka T, Hasegawa M, Tamatani T, Nakanishi T, Takigawa M, et al. Serum levels of connective tissue growth factor are elevated in patients with systemic sclerosis: association with extent of skin sclerosis and severity of pulmonary fibrosis. J Rheumatol. (2000) 27:149–54.
Keywords: PDGF, MRI, biomarker, Duchenne, Becker, dystrophinopathy
Citation: Alonso-Jiménez A, Fernández-Simón E, Natera-de Benito D, Ortez C, García C, Montiel E, Belmonte I, Pedrosa I, Segovia S, Piñol-Jurado P, Carrasco-Rozas A, Suárez-Calvet X, Jimenez-Mallebrera C, Nascimento A, Llauger J, Nuñez-Peralta C, Montesinos P, Alonso-Pérez J, Gallardo E, Illa I and Díaz-Manera J (2021) Platelet Derived Growth Factor-AA Correlates With Muscle Function Tests and Quantitative Muscle Magnetic Resonance in Dystrophinopathies. Front. Neurol. 12:659922. doi: 10.3389/fneur.2021.659922
Received: 02 February 2021; Accepted: 17 May 2021;
Published: 11 June 2021.
Edited by:Jens Schmidt, Universitätsklinikum Göttingen, Germany
Reviewed by:Melissa Hooijmans, VU University Medical Center, Netherlands
John Vissing, Rigshospitalet, Denmark
Copyright © 2021 Alonso-Jiménez, Fernández-Simón, Natera-de Benito, Ortez, García, Montiel, Belmonte, Pedrosa, Segovia, Piñol-Jurado, Carrasco-Rozas, Suárez-Calvet, Jimenez-Mallebrera, Nascimento, Llauger, Nuñez-Peralta, Montesinos, Alonso-Pérez, Gallardo, Illa and Díaz-Manera. 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: Jordi Díaz-Manera, email@example.com
†These authors have contributed equally to this work and share first authorship