The expressivity of Mendelian diseases can be influenced by factors independent from the pathogenic mutation: in Duchenne muscular dystrophy (DMD), for instance, age at loss of ambulation (LoA) varies between individuals whose DMD mutations all abolish dystrophin expression. This suggests the existence of trans-acting variants in modifier genes. Common single nucleotide polymorphisms (SNPs) in candidate genes (SPP1, encoding osteopontin, and LTBP4, encoding latent transforming growth factor β [TGFβ]-binding protein 4) have been established as DMD modifiers. We performed a genome-wide association study of age at LoA in a sub-cohort of European or European American ancestry (n = 109) from the Cooperative International Research Group Duchenne Natural History Study (CINRG-DNHS). We focused on protein-altering variants (Exome Chip) and included glucocorticoid treatment as a covariate. As expected, due to the small population size, no SNPs displayed an exome-wide significant p value (< 1.8 × 10-6). Subsequently, we prioritized 438 SNPs in the vicinities of 384 genes implicated in DMD-related pathways, i.e., the nuclear-factor-κB and TGFβ pathways. The minor allele at rs1883832, in the 5-untranslated region of CD40, was associated with earlier LoA (p = 3.5 × 10-5). This allele diminishes the expression of CD40, a co-stimulatory molecule for T cell polarization. We validated this association in multiple independent DMD cohorts (United Dystrophinopathy Project, Bio-NMD, and Padova, total n = 660), establishing this locus as a DMD modifier. This finding points to cell-mediated immunity as a relevant pathogenetic mechanism and potential therapeutic target in DMD.

OBJECTIVE: Respiratory insufficiency is a major complication of Duchenne muscular dystrophy (DMD). Its progression shows considerable interindividual variability, which has been less thoroughly characterized and understood than in skeletal muscle. We collected pulmonary function testing (PFT) data from a large retrospective cohort followed at Centers collaborating in the Italian DMD Network. Furthermore, we analyzed PFT associations with different DMD mutation types, and with genetic variants in SPP1, LTBP4, CD40, and ACTN3, known to modify skeletal muscle weakness in DMD. Genetic association findings were independently validated in the Cooperative International Neuromuscular Research Group Duchenne Natural History Study (CINRG-DNHS). METHODS AND RESULTS: Generalized estimating equation analysis of 1852 PFTs from 327 Italian DMD patients, over an average follow-up time of 4.5 years, estimated that forced vital capacity (FVC) declined yearly by -4.2%, forced expiratory volume in 1 sec by -5.0%, and peak expiratory flow (PEF) by -2.9%. Glucocorticoid (GC) treatment was associated with higher values of all PFT measures (approximately + 15% across disease stages). Mutations situated 3 of DMD intron 44, thus predicted to alter the expression of short dystrophin isoforms, were associated with lower (approximately -6%) PFT values, a finding independently validated in the CINRG-DNHS. Deletions amenable to skipping of exon 51 and 53 were independently associated with worse PFT outcomes. A meta-analysis of the two cohorts identified detrimental effects of SPP1 rs28357094 and CD40 rs1883832 minor alleles on both FVC and PEF. INTERPRETATION: These findings support GC efficacy in delaying respiratory insufficiency, and will be useful for the design and interpretation of clinical trials focused on respiratory endpoints in DMD.

BACKGROUND: Dilated cardiomyopathy (DCM) is a major complication of, and leading cause of mortality in Duchenne muscular dystrophy (DMD). Its severity, age at onset, and rate of progression display wide variability, whose molecular bases have been scarcely elucidated. Potential DCM-modifying factors include glucocorticoid (GC) and cardiological treatments, DMD mutation type and location, and variants in other genes. METHODS AND RESULTS: We retrospectively collected 3138 echocardiographic measurements of left ventricular ejection fraction (EF), shortening fraction (SF), and end-diastolic volume (EDV) from 819 DMD participants, 541 from an Italian multicentric cohort and 278 from the Cooperative International Neuromuscular Group Duchenne Natural History Study (CINRG-DNHS). Using generalized estimating equation (GEE) models, we estimated the yearly rate of decrease of EF (-0.80%) and SF (-0.41%), while EDV increase was not significantly associated with age. Utilizing a multivariate generalized estimating equation (GEE) model we observed that mutations preserving the expression of the C-terminal Dp71 isoform of dystrophin were correlated with decreased EDV (-11.01 mL/m2, p = 0.03) while for dp116 were correlated with decreased EF (-4.14%, p = <0.001). The rs10880 genotype in the LTBP4 gene, previously shown to prolong ambulation, was also associated with increased EF and decreased EDV (+3.29%, p = 0.002, and -10.62 mL/m2, p = 0.008) with a recessive model. CONCLUSIONS: We quantitatively describe the progression of systolic dysfunction progression in DMD, confirm the effect of distal dystrophin isoform expression on the dystrophin-deficient heart, and identify a strong effect of LTBP4 genotype of DCM in DMD.