Epidemiological research has firmly established the association of social context with late-life health and functioning (Cacioppo et al., Reference Cacioppo, Hughes, Waite, Hawkley and Thisted2006; Cohen, Reference Cohen2004). Yet this research unambiguously explicates neither the basis for these associations nor how social context relates to the biological and genetic factors known to contribute to later life functioning (Cohen & Janicki-Deverts, Reference Cohen and Janicki-Deverts2009; Hawkley & Cacioppo, Reference Hawkley and Cacioppo2003). Several longitudinal twin samples, located in different countries and employing somewhat different measures, have individually accumulated substantial data. Because recruiting and maintaining long-running studies of older adults is made especially difficult by participants’ relative frailty, these samples have generally been too small to reliably detect many of the kinds of effects likely to be of interest. Consequently, a collaboration among eight existing longitudinal twin and family studies was initiated, with a central focus on determining how social context is related to physical functioning (health, functional ability), and psychological functioning (well-being, cognition) in mid-life and older ages. Through this, we hope to lay the foundation for future studies of gene-environment interplay in late-life functioning. In the following we describe the rationale and design of the Interplay of Genes and Environment across Multiple Studies (IGEMS) consortium.
With their rich registry systems, the Nordic countries have pioneered many of the major longitudinal twin studies of aging, including the Longitudinal Study of Aging Danish Twins (LSADT) in Denmark (Christensen et al., Reference Christensen, Holm, McGue, Corder and Vaupel1999) the Swedish Adoption Twin Study of Aging (SATSA; Finkel & Pedersen, Reference Finkel and Pedersen2004), Origins of Variance in the Old-Old (OCTO-Twin; McClearn et al., Reference McClearn, Johansson, Berg, Pedersen, Ahern, Petrill and Plomin1997), and Ageing in Women and Men: A Longitudinal Study of Gender Differences in Health Behaviour and Health among Elderly (Gender; Gold et al., Reference Gold, Malmberg, McClearn, Pedersen and Berg2002) projects in Sweden. Moreover, investigators in both countries have extended their investigations to younger adult twin samples (e.g., the Middle-Age Danish Twin [MADT] study; Osler et al., Reference Osler, McGue, Lund and Christensen2008) and Twin and Offspring Study in Sweden (TOSS; Neiderhiser & Lichtenstein, Reference Neiderhiser and Lichtenstein2008), to provide foundations for longitudinal investigations of the influence of earlier developmental periods on aging outcomes. While undertaking similar longitudinal twin research in the United States has been more logistically challenging, members of IGEMS have established analogous studies in the United States, most notably the Minnesota Twin Study of Adult Development and Aging (MTSADA; Finkel & McGue, Reference Finkel and McGue1993) and the Vietnam Era Twin Study of Aging (VETSA; Kremen et al., Reference Kremen, Thompson-Brenner, Leung, Grant, Franz, Eisen and Lyons2006). Collectively, these studies include over 17,500 participants (including nearly 2,600 monogygotic [MZ] and 4,300 dizygotic [DZ] twin pairs and over 1,700 family members) whose ages ranged from 25 years to over 100 years at their intake assessments. Participants in these studies have undergone multiple longitudinal assessments (median follow-up period = 9.5 years) that include a rich variety of measures of aging-relevant outcomes in three broad domains: physical health and functional ability, psychological well-being (emotional stability/depression), and cognitive health. These studies also include multiple indicators of the social environment, spanning early childhood through late adulthood, and encompass participants with several different environmental contexts, such as differences in healthcare and retirement systems. The combined result is a set of longitudinal studies with assessments of outcomes and social contexts spanning mid- to late-adulthood.
The heuristic model that guides our investigation of gene-environment interplay in late-life functioning is presented in Figure 1. The model emphasizes the importance of two gene-environment processes. First, we hypothesize that environmental exposures and social contexts do not occur at random, but rather reflect an individual's genetically influenced behavior and choices (i.e., gene-environment correlation). We know, for example, that there are heritable influences on social factors including social support (Bergeman et al., Reference Bergeman, Neiderhiser, Pedersen and Plomin2001), social engagement (McGue & Christensen, Reference McGue and Christensen2007), and social isolation (Boomsma et al., Reference Boomsma, Cacioppo, Muthen, Asparouhov and Clark2007). We further hypothesize that gene-environment correlations are a major contributor to late-life phenotypic stability. For example, based on Schmalhausen (Reference Schmalhausen1946), we hypothesize that a shift from being socially engaged and physically active in middle age to being socially isolated in late-life will be accompanied by marked discontinuities in behavioral stability. Data in the IGEMS consortium allow us to chart the contributions of changes in social context to behavioral discontinuities; the twin-study design allows us to further explore the nature of gene-environment interplay in these transitions.
Second, we hypothesize that experiential factors can either diminish or amplify the influences of genetic effects on late-life outcomes (genotype-environment interaction, G × E). Two of Shanahan and Hofer's (Reference Shanahan and Hofer2005) four alternative models of G × E are directly relevant. Evolutionary models (Schmalhausen, Reference Schmalhausen1946) predict that the G × E form termed Contextual Triggering by Shanahan and Hofer, but recognized more broadly as the diathesis-stress model, will be more pronounced in late life. Specifically, during most of life, the typical human lives in an environment that is low in stress relative to the environments in which their species evolved (Owens, Reference Owens2002). In old age, biological, psychological, and social stresses are considered to be higher due to more negative change or declines. Thus, at times of relatively increased stress—such as old age—an organism's response to an unfavorable environment will depend on a unique combination of genes, so that genetic effects on physical manifestations of health will be maximized (Schmalhausen, Reference Schmalhausen1946). This pattern of findings is supported by earlier work by IGEMS investigators, showing that the heritability of physical health in mid-life is maximized when current income (Johnson & Krueger, Reference Johnson and Krueger2005a) and perceived life control (Johnson & Krueger, Reference Johnson and Krueger2005b) are low.
Alternatively, environmental factors can suppress genetic effects through a form of G × E Shanahan and Hofer call Social Context as Social Control. We predict that the Social Context as Social Control form of G × E is particularly relevant to extreme forms of physical disability or disease, which diminish the impact of genetic factors by reducing the opportunity for self-selection of activities that maintain stable functioning.
IGEMS Sample
Eight studies are involved in the IGEMS consortium. Brief descriptions of the studies are provided here and sample sizes are summarized in Table 1.
*Swedish and Danish studies have additional follow-ups through national registries.
SATSA = Swedish Adoption Twin Study of Aging; OCTO-Twin = Origins of Variance in the Old-Old; TOSS = Twin and Offspring Study in Sweden; LSDAT = Longitudinal Study of Aging Danish Twins; MADT = Middle-Aged Danish Twins; MTSADA = Minnesota Twin Study of Adult Development and Aging; VETSA = Vietnam Era Twins Study of Aging.
Swedish Studies
Four twin studies based in Sweden are included in IGEMS. The SATSA (Finkel & Pedersen, Reference Finkel and Pedersen2004) is a longitudinal twin study that began in 1984. Multiple measures of physical and psychological functioning and environmental/social context variables were assessed via interviews and questionnaires. Eight questionnaires have been sent to all study participants. Mean age at the first questionnaire assessment was 60.14 (SD = 14.02), 42.1% male, and all are Caucasians. Mean level of attained education was 8 years (SD = 2.38; range 6–16). Nine waves of in-person testing of twins over the age of 50 have been completed; the tenth wave is currently ongoing. The base population comprises all pairs of twins from the population-based Swedish Twin Registry, who indicated that they had been separated before the age of 10 and reared apart; and a sample of twins reared together, matched on the basis of gender, and date and county of birth. Individuals were included in the longitudinal study whether or not both members of the pair responded, and those who did not respond at the first wave were invited at subsequent waves. Zygosity was initially based on questions concerning similarity during childhood, administered when the Swedish Twin Registry was first compiled. These assignments were refined for participants in the first in-person testing-based serological markers, and further refined after genotyping (most recently with the cardio-metabochip).
The OCTO-Twin Study was the first major population-based longitudinal twin study directed to the oldest old (McClearn et al., Reference McClearn, Johansson, Berg, Pedersen, Ahern, Petrill and Plomin1997). The OCTO-Twin Study included five measurement occasions at 2-year intervals between 1991 and 2002 of 351 like-sexed twin pairs drawn from the Swedish Twin Registry who were over the age of 80 at baseline. Mean age at baseline assessment was 83.58 (SD = 3.17), 33.3% male, and all are Caucasians. In-person testing included assessment of physical and psychological functioning, and environmental/social context variables. The assessments were built on those from SATSA. Only pairs where both responded at baseline are included. Mean level of attained education was 7.09 years (SD = 2.24; range 0–25). Like other Swedish twin studies, zygosity was initially based on registry-based responses to questions regarding similarity. For those pairs with undetermined zygosity at the start of OCTO-Twin, zygosity was evaluated through Restriction Fragment Length Polymorphism (RFLP) markers. Like SATSA, some correction in zygosity has been made as further genotyping with the cardio-metabochip has become available.
Gender is a study of unlike-sex twin pairs born between 1906 and 1925 (Gold et al., Reference Gold, Malmberg, McClearn, Pedersen and Berg2002). A survey concerning health and health behavior was mailed in 1994 with responses from 1,210 twins from 605 pairs where both responded. Mean age at baseline questionnaire assessment was 74.43 (SD = 4.28) and all are Caucasians. Mean level of attained education was 7.28 years (SD = 2.4; range 5–25). A baseline in-person evaluation of 498 twins from 249 pairs between 70 and 80 years of age was undertaken between 1995 and 1997, and included an interview and tests of cognitive and physical functioning, similar to SATSA and OCTO-Twin. Two additional in-person waves followed at 4-year intervals. A baseline in-person evaluation of pairs 80 years and older took place between 1998 and 2002 with an interview and tests of cognitive and physical functioning similar to OCTO-Twin. Finally, a second survey was mailed in 2007 to all living twins who participated in the first mailed survey.
The TOSS was designed to examine how family relationships influence the adjustment of adult twins and their family members, how genetic and environmental factors may explain these associations, and to better specify genotype-environment correlation (Neiderhiser & Lichtenstein, Reference Neiderhiser and Lichtenstein2008). TOSS includes 909 pairs of same-sex twins, one adolescent child of each twin, and the spouse/partner of each twin. The average age is 47 (SD = 4.7 years, range 35–59) for the 350 male twin pairs and (SD = 4.6 years, range 32–60) for the 559 female twin pairs. Mean level of attained education was 13.28 years (SD = 2.74; range 9–16). Zygosity was determined by analyses of markers from DNA extracted from mouthwash or saliva samples. For the 54 twins who did not provide DNA, zygosity was determined by a questionnaire assessing twin similarity and by applying an algorithm base on a discriminant analysis of twins with DNA-confirmed zygosity. There were 254 MZ and 285 DZ twin mother pairs and 128 MZ and 183 DZ twin father pairs. The TOSS included a single assessment that indexed psychological functioning and extensive environmental/social context variables; however, longitudinal information is available from STR assessments. Only the twins and their spouses are included in the IGEMS.
Danish Studies
Two Danish studies are included in IGEMS. The LSADT (Christensen et al, Reference Christensen, Holm, McGue, Corder and Vaupel1999) began in 1995 with the assessment of members of like-sex twin pairs born in Denmark prior to 1920 (i.e., at least 75 years old). The surviving members of the initial cohort were followed up every 2 years in 1997, 1999, 2001, 2003, and 2005. Additional cohorts were added in 1997, 1999, and 2001 and the minimum age requirement was progressively reduced to age 70. Mean level of attained education was 8.65 years (SD = 2.87; range 5–18). Mean age at intake was 77.26 years (SD = 5.65; range 70–102) and 41% were males. Multiple measures of physical and psychological functioning and environmental/social context variables were assessed.
The study of MADT includes twins ranging in age from 46 to 68 years at the original assessment from 1998 to 1999 (Osler et al., Reference Osler, McGue, Lund and Christensen2008). Mean level of attained education was 11.10 years (SD = 2.96; range 7–18). Mean age at intake was 56.43 years (SD = 6.33) and 51% were males. A second assessment of MADT participants began in 2009. The MADT assessment closely parallels that used in LSADT.
The participants of LSADT and MDAT surveys are all Caucasians. Determination of zygosity in LSADT as well as MADT is based on the twins’ responses to four questions on physical similarity, a method that has previously been demonstrated to assign the true zygosity in more than 95% of the twin pairs (Christiansen et al., Reference Christiansen, Frederiksen, Schousboe, Skytthe, von Wurmb-Schwark, Christensen and Kyvik2003).
US Studies
The MTSADA (Finkel & McGue, Reference Finkel and McGue1993; McGue et al., Reference McGue, Hirsch and Lykken1993) is a population-based sample drawn from state birth records. Although MTSADA has a broad age range, two-thirds of the sample was age 60 years or older at intake (mean 58.6 years, SD = 10.8). Mean level of attained education was 13.31 years (SD = 2.61; range 0–25), 42% were males, and 98.5% were Caucasians, consistent with the demographics of Minnesota for the birth cohorts sampled. Participation involved two interview sessions separated by 1–3 weeks at which multiple physical functioning, psychological functioning, and environmental/social context variables were assessed. Interviews took place between 1984 and 1994. The MTSADA sample completed a follow-up mail survey approximately 6 years after initial assessment. Zygosity was assessed using five questions dealing with physical similarity, which previous research has shown yields approximately 95% accuracy when compared to genetic-based methods (Lykken et al., Reference Lykken, Bouchard, McGue and Tellegen1990).
The VETSA (Kremen et al., Reference Kremen, Thompson-Brenner, Leung, Grant, Franz, Eisen and Lyons2006) is a community-dwelling sample of male–male twin pairs, all of whom served in some branch of US military service at sometime between 1965 and 1975. The majority were not in combat. They now live throughout the United States and are generally comparable in terms of demographic and health characteristics to American men in their age range (Schoeneborn & Heyman, Reference Schoeneborn and Heyman2009). Wave 1 interviewing took place between 2003 and 2007. The design is somewhat unusual in that the VETSA has a rather narrow age range. At the time of recruitment for the Wave 1 assessment, all participants were in their 50s (although two pairs turned 60 by the time they came to the laboratory). The mean level of formal education completed was 13.84 years (SD = 2.11; range: 8–20). Eighty-six percent were Caucasian. Zygosity was determined on the basis of 25 microsatellite markers for 92% of the sample, and by questionnaire and blood group information for the remaining 8%. The latter method agreed with the DNA-based results in 95% of cases. Both twins in a pair had to agree to participate in order to be enrolled in Wave 1, but there is no such requirement for follow-up assessments. Participants filled out multiple questionnaires at home about physical and psychological functioning, lifestyle, and environmental/social context. They also underwent a full day of in-person assessments (in San Diego or Boston according to participant preference) focused on cognitive and physical function. Participants who could or would not come to one of the laboratories were tested in their hometowns. Assessments are currently planned for every 5 years, and the Wave 2 assessment is nearly complete (as of September 2012) with a retention rate of approximately 85%. A key goal of the VETSA is to make use of its age-homogenous cohort to understand individual differences in aging trajectories in the transition from mid- to later life.
IGEMS Measures
There is a wealth of data available from the eight studies, predominantly from the series of measurement occasions within each longitudinal study. Table 2 shows the key domains from the heuristic model and the measures available across the eight studies to assess those domains. As indicated in the table, these measures were largely obtained by mailed questionnaire or interview, with some physical or functional variables measured during in-person visits. In addition to these data, national registry data available in Sweden and Denmark allow us to obtain supplemental objective measures of health and early environment (such as birth weight). Furthermore, the twin registries also have information provided by the twins when the registries were compiled (in the 1960s and 1970s).
*Limited or available on a subset.
Data Harmonization
Given the differences in study focus, design, and specific measures assessed, we have placed strong emphasis on data harmonization, to ensure the greatest possible consistency of measurement across studies. We developed a common administrative file structure that includes demographic variables, last known vital status, age at each time of assessment, and reasons for non-participation if the twin did not participate at that wave. Because of inter-study differences in how cohorts were defined and schedule of data collection, it was crucial to come to common definitions that would allow combining studies longitudinally.
At the same time, we undertook harmonization of relevant phenotypes and outcomes to enable the combined analysis of the multiple sets of data. We are not simply recoding to a reduced set of common categories, but rather are identifying overlapping item content and response formats, applying psychometric analysis to establish measurement invariance via item response theory (IRT)-factor analytic approaches, and conducting meta-analyses and integrated data analysis of pooled data where warranted (Bauer & Hussong, Reference Bauer and Hussong2009; Curran & Hussong, Reference Curran and Hussong2009; Rosenthal, Reference Rosenthal1991; Rosenthal & Di Matteo, Reference Rosenthal and DiMatteo2001).
As a first step, we created spreadsheets for measures that correspond to the constructs in the model, including the specific questions and response options. Table 2 indicates which variables were available; detailed spreadsheets show how each variable was assessed (see example in Figure 2).
Because integrative data analysis (IDA), where data are pooled across samples to increase power and finely consider between-study heterogeneity, is the preferred method (Curran & Hussong, Reference Curran and Hussong2009), we are creating scale scores that are common across studies. This requires overlapping item content across the studies as well as across time for longitudinal hypotheses. Thus, a second step entailed identifying those measures for which it was possible to identify a common metric; for example, Body Mass Index (BMI), lung function, and blood pressure. For harmonizing education, we have created year equivalents for educational levels among Danish and Swedish studies that did not collect the data in years.
Where a common metric is not already available, IRT techniques are implemented to create harmonized scores across studies (cf. Bauer & Hussong, Reference Bauer and Hussong2009). The premise is that trait scores are estimable based on individual responses to and features of items (e.g., item ‘difficulty’, discrimination and reliability; Embretson & Reise, Reference Embretson and Reise2000). Moreover, in the case of IDA, corresponding indices are available across samples with items that invariantly tap a particular construct(s) (Bauer & Hussong, Reference Bauer and Hussong2009). Multiple common items should be available to test that item functioning is similar across samples (Bauer & Hussong, Reference Bauer and Hussong2009), and to establish ‘crosswalks’ between the different measurement indices (e.g., Velozo et al., Reference Velozo, Byers, Wang and Joseph2007).
For example, all Swedish studies (SATSA, OCTO-Twin, Gender, TOSS) and the US studies (MTSADA and VETSA) measured depressive symptomatology using the Center for Epidemiologic Studies – Depression scale (CES-D) (Radloff, Reference Radloff1977), while the Danish studies used an adapted version of the Cambridge Mental Disorders of the Elderly Examination (CAMDEX; see McGue & Christensen, Reference McGue and Christensen1997). There was no item with the same wording and the same response format on the two scales. However, both scales appeared to encompass the same conceptualization of depression. The CES-D can be analyzed as a single score or factored into four sub-scales (Gatz et al., Reference Gatz, Johansson, Pedersen, Berg and Reynolds1993): depressed mood, somatic complaints, interpersonal difficulties, and lack of well-being. The 21 CAMDEX items can be analyzed as a single score or as two-factor analytically derived sub-scales: an affect sub-scale and a somatic sub-scale. Substantially similar results have been reported across studies with respect to the relative roles of genetic and environmental influences on depressive symptoms (Gatz et al., Reference Gatz, Pedersen, Plomin, Nesselroade and McClearn1992; McGue & Christensen, Reference McGue and Christensen1997) and the possibility that heritability may be greater in women than men (Jansson et al., Reference Jansson, Gatz, Berg, Johansson, Malmberg, McClearn and Pedersen2004; McGue & Christensen, Reference McGue and Christensen1997). Since we could find no published article where the same sample had taken both measures, as a third step in data harmonization, we recruited a separate new sample and administered both scales to that sample, after which it was possible to use IRT methods to establish a ‘crosswalk’ to apply to the scores to create a common scale.
We also used this approach with perceived socio-economic status (SES) and self-rated health items, where the variations on how each was assessed were extensive across studies. Similar to the depression example, we have used IRT modeling to calibrate a common metric for a subset of perceived SES items measured in three Swedish studies (SATSA, OCTO-Twin, and Gender) and the Minnesota study (MTSADA) by using recruitment samples that were administered a subset of items in common among and unique to the Swedish and Minnesotan samples. Specifically, we evaluated a total of six items from the Swedish and Minnesota studies. Examples of perceived SES items included: (1) ‘Do you have difficulty in meeting general expenses like food, rent, bills etc?’ (No/Yes); (SATSA, Gender); (2) ‘I have enough finances to pay monthly bills’ (Yes/No); (MTSADA). Additional recruitment samples will be given full sets of perceived SES items across all available Swedish and US samples, adding items from TOSS and VETSA.
In a third example, we have examined similar self-rated health items included in all eight studies. While item prompts are similar, response metrics varied significantly, ranging from two to as many as seven response options. All combinations of items and response options were collected from an independent sample, and these data were used to develop latent self-rated health factors that could be combined across studies.
At the end of this process, we will have a common set of measures, each from as many studies as possible, to pursue testing the heuristic model presented above. Where integrative analyses cannot be used, we are relying on meta-analytic methods (Rosenthal, Reference Rosenthal1991; Rosenthal & Di Matteo, Reference Rosenthal and DiMatteo2001), using fixed- and random-effects approaches where appropriate.
Conclusion
The expansion of the elderly population has brought with it challenges for healthcare and social delivery systems, which may become overburdened with large numbers of disabled and cognitively impaired individuals (Christensen et al., Reference Christensen, Doblhammer, Rau and Vaupel2009). There is a clear need to understand the basis for healthy aging and especially to identify and characterize the nature of the modifiable risk and protective factors associated with later-life functioning (Christensen et al., Reference Christensen, McGue, Petersen, Jeune and Vaupel2008). IGEMS will support a series of epidemiological, genetic, and biological investigations of one class of potentially modifiable factors, those related to social context.
Acknowledgments
The IGEMS consortium acknowledges the contributions of Chris J. Hahn, Marielena De Araujo-Held, and Catalina Zavala. Most importantly, the authors gratefully acknowledge the continued cooperation and participation of the twins and their families.
IGEMS is supported by the National Institutes of Health grant no. R01 AG037985. SATSA was supported by grants R01 AG04563, R01 AG10175, the MacArthur Foundation Research Network on Successful Aging, the Swedish Council for Working Life and Social Research (FAS) (97:0147:1B, 2009-0795) and Swedish Research Council (825-2007-7460, 825-2009-6141). OCTO-Twin was supported by grant R01 AG08861. Gender was supported by the MacArthur Foundation Research Network on Successful Aging, The Axel and Margaret Ax:son Johnson's Foundation, The Swedish Council for Social Research, and the Swedish Foundation for Health Care Sciences and Allergy Research. TOSS was supported by grant R01 MH54610 from the National Institute of Health. The Danish Twin Registry is supported by grants from The National Program for Research Infrastructure 2007 from the Danish Agency for Science and Innovation, the Velux Foundation and the US National Institute of Health (P01 AG08761). The Minnesota Twin Study of Adult Development and Aging was supported by NIA grant R01 AG 06886. VETSA was supported by National Institute of Health grants R01 AG018384, R01 AG018386, R01 AG022381, and R01 AG022982, and, in part, with resources of the VA San Diego Center of Excellence for Stress and Mental Health. The Cooperative Studies Program of the Office of Research & Development of the United States Department of Veterans Affairs has provided financial support for the development and maintenance of the Vietnam Era Twin (VET) Registry.