Our analytical sample comprised of 628 individuals, of whom 220 (35%) were men and 408 (64%) women. Women were older than men (women: 83.61 (SD=3.17) yrs. old vs men: 82.95 (SD=2.76) yrs. old, t-test, p-value=0.007) and had fewer years of education than men (women: 6.96 (SD=1.89), men: 7.51 (SD=2.90), t-test, p-value=0.01). More than 73% of the sample (n=463) remained free of dementia over the study follow-up period. See Table 1 for descriptive sample characteristics. The sample included 272 (43%) monozygotic and 356 (57%) dizygotic individuals. Sixty-three individuals (10% of total sample, 24 monozygotic and 39 dizygotic twins) were excluded due to missing data for their co-twin. Monozygotic twins were slightly more educated than dizygotic twins (7.50 (SE=0.16) years of education vs 6.89 (SE=0.10) years of education respectively, t-test, p-value=0.0009) and a larger proportion of monozygotic twins were sedentary compared to dizygotic twins (63% vs 53% respectively, chi square test, p-value=0.02). Sixty per cent (60%) of monozygotic and 68% of dizygotic individuals were women.
Inspection of BIC, SSABIC, and the Lo-Mendel-Rubin likelihood ratio test values resulted in the identification of a 4-class model as the best fitting model. The entropy for this model was 0.74, a value that suggests a good discrimination of individuals into these 4 classes.
Table 1 Descriptive statistics of the OCTO-Twin sample
Variables
|
N
|
Mean (SD)
|
Baseline characteristics
|
N (%)
|
MMSE 1
|
628
|
26.3 (3.9)
|
Monozygotic
|
272 (43.0)
|
MMSE 2
|
524
|
24.9 (6.2)
|
Female
|
408 (65.0)
|
MMSE 3
|
404
|
24.4 (7.1)
|
Smoker
|
248 (39.5)
|
MMSE 4
|
298
|
23.7 (7.5)
|
Stimulates the body
|
266 (42.4)
|
MMSE 5
|
217
|
22.3 (7.8)
|
Stimulates brain
|
416 (66.3)
|
Education (years)
|
|
7.5 (2.3)
|
Diagnosed with dementia
|
165 (26.3)
|
MMSE 1-5 are the MMSE scores at every follow-up. SD standard deviation
MMSE trajectory classes and characteristics
First, we examine the characteristics revealed in the 4 class model (see Table 2). In Figure 1 we show a graphical illustration of the observed and estimated MMSE trajectories in our analytical sample.
Table 2 Results from the 4- class growth mixture model that best fitted the OCTO-Twin sample
|
High performers and stable
N=198
|
High performers with slow decline N=194
|
High performers with fast decline
N=181
|
Impaired very fast decline
N=55
|
|
β (SE)
|
p
|
β (SE)
|
p
|
β (SE)
|
p
|
β (SE)
|
p
|
Fixed effects
|
|
|
|
|
|
|
|
|
MMSE Level
|
28.9 (0.14)
|
0.00
|
27.62 (0.32)
|
0.00
|
26.34 (0.59)
|
0.00
|
19.30 (2.18)
|
0.00
|
Baseline age
|
-0.08 (0.03)
|
0.00
|
-0.16 (0.07)
|
0.02
|
-0.18 (0.12)
|
0.12
|
0.28 (0.23)
|
0.24
|
Education
|
0.04 (0.02)
|
0.04
|
0.09 (0.08)
|
0.26
|
0.30 (0.17)
|
0.09
|
0.73 (0.62)
|
0.23
|
Women
|
0.26 (0.13)
|
0.06
|
0.24 (0.37)
|
0.51
|
-0.15 (0.72)
|
0.83
|
2.71 (2.26)
|
0.23
|
Dementia
|
-0.50 (0.09)
|
0.02
|
0.13 (0.42)
|
0.75
|
-1.18 (0.74)
|
0.11
|
-4.29 (1.28)
|
0.00
|
Linear Slope
|
-0.04 (0.04)
|
0.30
|
-0.31 (0.09)
|
0.00
|
-0.96 (0.14)
|
0.00
|
-2.95 (0.86)
|
0.00
|
Baseline age
|
0.00 (0.00)
|
0.60
|
-0.04 (0.08)
|
0.82
|
-0.02 (0.03)
|
0.47
|
-0.02 (0.09)
|
0.77
|
Education
|
0.00 (0.00)
|
0.42
|
-0.01 (0.02)
|
0.51
|
-0.08 (0.05)
|
0.11
|
-0.84 (0.53)
|
0.10
|
Women
|
0.02 (0.04)
|
0.61
|
0.03 (0.10)
|
0.72
|
0.10 (0.18)
|
0.57
|
-0.06 (1.03)
|
0.95
|
Dementia
|
-0.11 (0.09)
|
0.23
|
-0.76 (0.14)
|
0.00
|
-1.40 (0.18)
|
0.00
|
-0.76 (0.52)
|
0.15
|
Random effects
|
|
|
|
|
|
|
|
|
Level
|
0.06 (0.08)
|
0.42
|
0.06 (0.08)
|
0.42
|
0.06 (0.08)
|
0.42
|
0.06 (0.08)
|
0.42
|
Linear Slope
|
0.02 (0.00)
|
0.00
|
0.02 (0.00)
|
0.00
|
0.02 (0.00)
|
0.00
|
0.02 (0.00)
|
0.00
|
Error
|
0.66 (0.06)
|
0.00
|
5.30 (0.48)
|
0.00
|
22.24 (2.43)
|
0.00
|
27.11 (3.72)
|
0.00
|
High Performers and Stable class
The largest class was comprised of 32% of the sample (n=198). A reference individual in this class is a man aged 83 years at study entry with 7 years of education with an average MMSE score at study entry of 28.9 (SE=0.1), and who remained free of dementia over the study period. The annual rate of MMSE change was only -0.04 (SE=0.04), an estimate that did not reach conventional significance thresholds. In this class, older baseline age was associated with a lower baseline MMSE performance (ß=-0.08 (SE=0.03), while education was associated with higher baseline MMSE score (ß=0.04 (SE=0.02)). Individuals who later received a diagnosis of dementia had poorer MMSE performance at study entry than individuals who remained free of dementia (ß=-0.48, SE=0.10))
High Performers with Slow Decline class
Thirty-one per cent (31%, n=194) of the sample were assigned to a class with an average MMSE score of 27.6 (SE=0.3) at study entry. Their annual decline was -0.3(0.09) MMSE points. Age was associated with poorer baseline MMSE performance, and those who ever received a diagnosis of dementia declined at a faster rate than those who remained free from dementia.
High Performers with Fast Decline class
Twenty-nine per cent (29%, n=181) of the sample was assigned to a class with a slightly lower baseline MMSE score than the other class of high performing individuals, but they declined at a faster rate. Specifically, in this class, the reference individual had an average MMSE score at study entry of 26.3 (SE=0.6) with an annual decline of -0.9 (SE=0.1) MMSE points. Individuals who received a diagnosis of dementia declined 1.4 MMSE points faster per year than those who remained non-demented. None of the risk factors emerged as significantly associated with MMSE baseline performance, nor with rate of change.
Impaired Very Fast Decline class
Finally, 8 % (n=55) of the sample was assigned to a class characterized by low baseline scores and more substantial decline. Their average MMSE score at study entry was 19.3 (SE=2.2), followed by an annual decline of -2.9 (SE=0.8). Individuals who received a diagnosis of dementia had even poorer baseline MMSE performance, compared with those who remained free from dementia. None of the other risk factors reached significance levels.
Risk factors and class assignment probability
Next, we examine the effect of risk factors in relation to the above classification (see Table 3). Compared with the high performing and stable class and as expected, older age at study entry, was associated with higher odds of being in the class of impaired performers with very fast decline. This was also the case in the two classes of high performing individuals. Women and more educated individuals had lower odds of being in the class of High Performers with Fast decline individuals, than in the High Performers and Stable class of individuals. Individuals who did not engage in physical activity were more likely to be in the High Performing with Fast Decline class, compared to the High Performing and Stable class of individuals, a finding that suggests a partial protective effect of education and physical activity.
Table 3 Odds ratio results from multinomial model for class assignment with reference to the High Performers and Stable class
|
High performers with slow decline
|
High performers with fast decline
|
Impaired very fast decline
|
|
OR (SE)
|
p
|
OR (SE)
|
p
|
OR (SE)
|
p
|
Baseline age
|
1.12 (0.06)
|
0.05
|
1.17 (0.07)
|
0.01
|
1.22 (0.09)
|
0.00
|
Education
|
0.92 (0.05)
|
0.11
|
0.75 (0.08)
|
0.00
|
0.75 (0.11)
|
0.03
|
Women
|
0.88 (0.28)
|
0.67
|
0.51 (0.18)
|
0.00
|
1.04 (0.62)
|
0.93
|
Smoking
|
1.18 (0.34)
|
0.60
|
0.96(0.32)
|
0.92
|
1.08 (0.57)
|
0.88
|
Stimulates the body
|
1.12 (0.33)
|
0.71
|
2.28 (0.60)
|
0.03
|
5.20 (2.36)
|
0.07
|
Stimulates the brain
|
1.31 (0.36)
|
0.38
|
1.89 (0.54)
|
0.10
|
6.73 (4.84)
|
0.23
|
OR odds ratio; SE standard error
Distribution of twins across classes
In the next step, we examine the distribution of twins and partners across the four classes. In the high performing and stable class there were 47 twin pairs (i.e., 94 individuals out of 198). In the classes of High Performers with Slow Decline and High Performers with Fast Decline, there were 27 twin pairs (54 individuals out of 194) and 23 pairs (46 individuals out of 181), respectively. In the class of Impaired individuals with Very Fast Decline, there were 4 pairs of twins (8 individuals out of 55). These numbers indicate that 47.5 % of the individuals in the High Performers and Stable Class in fact had a twin partner in the same class. The corresponding proportion of having a twin in the same class was 23.7% among the High Performers with Slow Decline, and 29.8% in the class of High Performers with Fast Decline. Lastly, in the Impaired and Fast Declining class 14.4 % of individuals had a co-twin in the same class.
In the final step we conducted posteriori analyses, to test whether zygosity was associated with the likelihood of being assigned to a certain class (see Figure 2). Our analyses, however, failed to find evidence in support of an association (X2(3 p=0.46), which means that the likelihood for being in a certain class was similar for MZ and same-sex DZ twins.