Alexander disease (AxD) is an extremely rare neurodegenerative disorder caused by glial fibrillary acidic protein (GFAP) gene mutations. Compared with the cerebral type, which is characterized by infantile onset, the bulbospinal type and intermediate form are associated with a late onset, spanning from juveniles to the elderly, and more diverse clinical spectrum, suggesting the existence of factors contributing to phenotypic diversity. To build a foundation for future genetic studies of this rare disease, we obtained genomic data by whole exome-sequencing (WES) and DNA microarray derived from thirty-one AxD patients with the bulbospinal type and intermediate form. Using this data, we aimed to identify genetic variations determining the age at onset (AAO) of AxD. As a result, WES- or microarray-based association studies between younger (<45 years; n = 13)- and older (≥45 years; n = 18)-onset patients considering the predicted GFAP-mutation pathogenicity identified no genome-wide significant variant. The candidate gene approach identified several variants likely correlated with AAO (p < 0.05): GAN, SLC1A2, CASP3, HDACs, and PI3K. Although we need to replicate the results using an independent population, this is the first step towards constructing a database, which may serve as an important tool to advance our understanding of AxD.