Genome-wide association studies (GWAS) underscore the genetic basis of multiple sclerosis (MS); however, only few of the newly reported genetic variations relevant in MS have been replicated or correlated for clinical/paraclinical phenotypes such as spinal cord atrophy in independent patient cohorts. We genotyped 141 MS patients for 58 variations reported to reach significance in GWAS. Expanded disability status scale (EDSS) and disease duration (DD) are available from regular clinical examinations. MRI included sagittal high-resolution 3D T1-weighted magnetization-prepared rapid acquisition gradient echo of the cervical cord region used for volumetry. Due dependency of mean upper cervical cord area (MUCCA) with EDSS and/or DD, correction operations were performed compensating for EDSS/DD. We assessed each MS risk locus for possible MUCCA association. We identified twelve risk loci that significantly correlated with MUCCA. For nine loci-BATF, CYP27B1, IL12B, NFKB1, IL7, PLEK, EVI5, TAGAP and nrs669607-patients revealed significantly higher degree of atrophy; TYK2, RGS1 and CLEC16A revealed inverse effects. The weighted genetic risk score over the twelve loci showed significant correlation with MUCCA. Our data reveal a risk gene depending paraclinical/clinical phenotype. Since MUCCA clearly correlates with disability, the candidates identified here may serve as prognostic markers for disability progression.