Identification of eczema clusters and their association with filaggrin and atopic comorbidities: analysis of five birth cohorts

Background: Longitudinal modelling of the presence/absence of current eczema through childhood has identified similar phenotypes, but their characteristics often differ between studies.

Objectives: To demonstrate that a more comprehensive description of longitudinal pattern of symptoms may better describe trajectories than binary information on eczema presence.

Methods: We derived six multidimensional variables of eczema spells from birth to 18 years of age (including duration, temporal sequencing and the extent of persistence/recurrence). Spells were defined as consecutive observations of eczema separated by no eczema across 5 epochs in five birth cohorts: infancy (first year); early childhood (age 2-3 years); preschool/early school age (4-5 years); middle childhood (8-10 years); adolescence (14-18 years). We applied Partitioning Around Medoids clustering on these variables to derive clusters of the temporal patterns of eczema. We then investigated the stability of the clusters, within-cluster homogeneity and associated risk factors, including FLG mutations.

Results: Analysis of 7464 participants with complete data identified five clusters: (i) no eczema (51.0%); (ii) early transient eczema (21.6%); (iii) late-onset eczema (LOE; 8.1%); (iv) intermittent eczema (INT; 7.5%); and (v) persistent eczema (PE; 11.8%). There was very-high agreement between the assignment of individual children into clusters when using complete or imputed (n = 15 848) data (adjusted Rand index = 0.99; i.e. the clusters were very stable). Within-individual symptom patterns across clusters confirmed within-cluster homogeneity, with consistent patterns of symptoms among participants within each cluster and no overlap between the clusters. Clusters were characterized by differences in associations with risk factors (e.g. parental eczema was associated with all clusters apart from LOE; sensitization to inhalant allergens was associated with all clusters, with the highest risk in the PE cluster). All clusters apart from LOE were associated with FLG mutations. Of note, the strongest association was for PE [relative risk ratio (RRR) 2.70, 95% confidence interval (CI) 2.24-3.26; P < 0.001] followed by INT (RRR 2.29, 95% CI 1.82-2.88; P < 0.001).

Conclusions: Clustering of multidimensional variables identified stable clusters with different genetic architectures. Using multidimensional variables may capture eczema development and derive stable and internally homogeneous clusters. However, deriving homogeneous symptom clusters does not necessarily mean that these are underpinned by completely unique mechanisms.