  Neurodegenerative diseases are characterized by numerous markers of
progression and clinical endpoints. For instance, Multiple System Atrophy
(MSA), a rare neurodegenerative synucleinopathy, is characterized by various
combinations of progressive autonomic failure and motor dysfunction, and a very
poor prognosis. Describing the progression of such complex and
multi-dimensional diseases is particularly difficult. One has to simultaneously
account for the assessment of multivariate markers over time, the occurrence of
clinical endpoints, and a highly suspected heterogeneity between patients. Yet,
such description is crucial for understanding the natural history of the
disease, staging patients diagnosed with the disease, unravelling
subphenotypes, and predicting the prognosis. Through the example of MSA
progression, we show how a latent class approach modeling multiple repeated
markers and clinical endpoints can help describe complex disease progression
and identify subphenotypes for exploring new pathological hypotheses. The
proposed joint latent class model includes class-specific multivariate mixed
models to handle multivariate repeated biomarkers possibly summarized into
latent dimensions and class-and-cause-specific proportional hazard models to
handle time-to-event data. Maximum likelihood estimation procedure, validated
through simulations is available in the lcmm R package. In the French MSA
cohort comprising data of 598 patients during up to 13 years, five
subphenotypes of MSA were identified that differ by the sequence and shape of
biomarkers degradation, and the associated risk of death. In posterior
analyses, the five subphenotypes were used to explore the association between
clinical progression and external imaging and fluid biomarkers, while properly
accounting for the uncertainty in the subphenotypes membership.
