Questionnaire and Interview Studies

Earlier questionnaires were mailed in 1975 and 1981 to all twins born before 1958 and living in Finland, identified from the population register of Finland; the 1990 third-wave questionnaire was targeted to twins born in 1930–1957 who had responded to either earlier questionnaire. Numbers and details are shown in the first three rows of Table 1.

Table 1. Questionnaires and interviews of the adult twin cohort in four waves of data collection. The fourth wave was spread over two decades and used either interviews or questionnaires for data collection in three separate substudies defined by birth year

Note: I = Interview, Q = Questionnaires

^a At wave 1, these numbers include nonbiological twins who satisfied the original selection criteria (pairs of persons with the same surname at birth, same birth date, same local community of birth and same sex) in 1974.

The fourth wave of cohortwide data collection has been conducted in three separate phases and with separate data collection methods. We invited the participation of twins based on their birth years (those born before 1938, 1938–1944 and 1945–1957) for studies conducted between 1999 and 2017. In addition, selected subsets of twins have been approached for more intensive clinical studies.

Cognitive Screening Studies

Starting in 1999, twins from the FTC over 65 years old, and eventually including all those born before 1938, were invited to participate in a telephone interview assessing cognitive functioning. Two validated instruments were used: a telephone assessment of dementia (TELE; Gatz et al., Reference Gatz, Reynolds, John, Johansson, Mortimer and Pedersen2002) and Telephone Interview for Cognitive Status (TICS; Brandt et al., Reference Brandt, Spencer and Folstein1988). Also, a 1-min category (animals) fluency task was administered. The first part of the data collection took place in 1999–2007 and included twins born before 1938. More specifically, monozygotic (MZ) twin-pairs (both members of a pair alive) were interviewed in 1999–2001 and dizygotic (DZ) twins and those with uncertain zygosity, irrespective of the co-twin’s vital status, were interviewed in 2003–2007. All participants were also asked to provide a blood sample for extracting DNA. There were a total of 2606 twins who participated in the telephone interview at a mean age of 74 years. The study collection is termed MEMTWIN I in Table 1.

The second part of the data collection was performed in 2013–2017 and included cohorts born in 1938–1944 (MEMTWIN II in Table 1). All twins — regardless of co-twin’s vital status — were asked to participate. Cognitive measures were similar to earlier data collection, with the exception of using the modified version of the TICS (TICS-m, including a delayed free recall of 10-word list). In addition to cognitive instruments, the Center for Epidemiological Studies Depression scale (CES-D) was also included (Radloff, Reference Radloff1977). All individuals were also asked to participate in a saliva sample collection for extraction of DNA. A total of 1813 twins participated in the telephone interview at a mean age of 74 years. Twins born in 1940–1944 were invited to take part in objective measures of physical activity. See below for details and first results.

These two data collections, collectively called MEMTWIN, have cognitive data available from a total of 4429 twins with a participation rate of 67% (Lindgren et al., Reference Lindgren, Kaprio, Rinne and Vuoksimaa2019). Of these, 3614 have genotype data. Informed consent was obtained from participants, and the study protocol was approved by the Ethics Committee of the Hospital District of Southwest Finland. We are currently planning to administer postal questionnaires and telephone interviews for all FTC twins who reach 90 years of age. Selected pairs potentially discordant for cognitive impairment were invited to the national PET Centre, University of Turku for neuropsychological (NP) assessments, magnetic resonance imaging (MRI) and positron emission tomography (PET) imaging.

Lifestyle, Health and Behavioral Changes Over 35 Years

In the 2013 review (Kaprio, Reference Kaprio2013), the mailed survey conducted among twins born in 1945–1957 was described in detail. In brief, the questionnaire had 85 questions (and 240 items). The topics covered items about contact with their co-twin, health and risk factors (blood pressure [BP] and cholesterol measurements, diabetes, breathlessness, back pain, bruxism, chronic diseases and use of common medications). Women were asked about hormone replacement therapy (HRT) use. Six questions on sleep were included. Eight items on work and pension intentions were asked next. Twelve questions asked a detailed smoking history, including two nicotine dependence scales, while seven questions were directed to alcohol quantity, frequency and problem drinking. Six items asked about leisure-time physical activity (LTPA) and sedentary behavior. Life events and traumatic experiences in childhood and adulthood were questioned, as well as items on social interactions. We repeated the four-item life satisfaction scale, and the brief assessment of extraversion and neuroticism (19 items together) asked in previous surveys. We included the CES-D depression scale. Finally, we asked about their current height and weight, as well as weight from 1 year and 5 years earlier, as well as their weight at age 50 and largest weight ever. Finally, the participants were asked to self-measure their waist circumference with a tape measure included with the questionnaire. We found that this had acceptably high validity when measured by health professionals in the Epigenetics of Hypertension study (Tuomela et al., Reference Tuomela, Kaprio, Sipilä, Silventoinen, Wang, Ollikainen and Piirtolain press). This has enabled longitudinal analyses spanning up to four measurements and 35 years of follow-up (see below for more detail on analyses related to alcohol, weight and physical inactivity). The details of the three subsets of wave 4 are also shown in Table 1, while Figure 1 gives the numbers of twins who have answered various combinations of surveys, which helps when planning longitudinal analyses.

Fig. 1. Overlap in responses to waves 1–4 studies for the FTC. Panel A: The entire cohort; Panel B: twins born before 1938 and assessed for cognition in wave 4 (older part of MEMTWIN); Panel C: twins born 1938–1944 and assessed for cognition in wave 4 (young part of MEMTWIN) and Panel D: twins born in 1945–1957.

Epigenetics of Hypertension Study

Based on the 2011 survey, a few years later, we undertook an in-person study on the epigenetics of hypertension (EH-Epi study, see below) as part of a consortium project led by Dr Xiaoling Wang (Georgia Prevention Institute, Augusta University, GA, USA). A total of 447 twins were examined and measured, including four measurements of BP at two separate time points, weight, height as well as waist and hip circumferences, as well as fasting glucose and cell counts (60 items). A lifestyle interview had items on current medications, use of alcohol, smoking, occupation, physical activity at leisure and handedness (160 items). The participants also completed a 56-item brief food frequency questionnaire.

Blood samples were taken for genetic and epigenetic studies (DNA methylation using 450-k arrays, next-generation sequencing for miRNA and histone modification analyses), as well as for metabolomics on several platforms and transcriptomics. See below for more details and initial results.

Inheritance of Drinking Behavior Cohort

While the older FTC is large and comprehensive, it is not the first in Finland. There is a well-known study of male twins born in the 1920s who took part in a study of genetics of alcohol use and abuse, summarized in the book Inheritance of Drinking Behavior, by Juha Partanen, Kettil Bruun and Touko Markkanen (Partanen et al., Reference Partanen, Bruun and Markkanen1966). They studied in detail 902 male twin-pairs aged 28–37 years in 1958, and we have now been able to access this interview data. We are currently re-entering the paper interviews as the previous electronic records no longer exist. The interview was extensive with 182 items on biographical data, relationships between co-twins, external circumstances, drinking behavior, illnesses, food and taste preferences and anthropological measures. In addition, intelligence and personality tests were conducted and blood samples drawn. As many of the pairs have also taken part in our 1975 and 1981 surveys, this dataset will add data from one more data point to the longitudinal study.

Register Linkages

Through regular register linkages, we have updated vital status and address information as needed, while cancer incidence and cause of death data are updated from the Cancer Registry and Statistics Finland, respectively. For more details on ongoing cancer research, see the review on the Nordic Twin Study of Cancer (Mucci et al., Reference Mucci, Hjelmborg, Harris, Czene, Havelick and Scheike2016). In addition, some analyses of predictors of cancer incidence have been done on the Finnish data alone. More extensive register linkages have been waiting for new legislation on accessing health and social data. The number of deaths and number of cancer cases (all kinds of malignancies) are shown in Table 2 for the like-sex twins who have replied to the specified surveys. This enables planning of analyses of exposure data collected at the surveys in relation to cause of death or cancer outcomes.

Table 2. Number of deaths, cancer cases in the twin cohort by survey year and the numbers genotyped or with methylation data. The number of pairs genotyped on a genomewide array

The Finnish Parliament has also approved the Act on the Secondary Use of Health and Social Data in 2019. It enables wider and speedier use of the data not only from national social and health care registers but also from the operational client and patient systems in primary care, specialist health care and social services. Decisions on accessing the data will be taken by the new Data Permit Authority as a centralized national actor. For sensitive data, it will be required to be handled in a safe and secure environment, such that only the results of the analytics can be used externally. Potentially, this will permit linkage of the twin cohort to multiple new registers for novel analyses. The new regulation implements the GDPR (the EU’s General Data Protection Regulation) in Finland. It is one of the first implementations of the GDPR for the secondary use of data in Europe.

Biobanking

As part of the Finnish Biobanking Legislation (law 688/2012), all biological samples associated with the FinnTwin12 study, the FinnTwin16 study and the older FTC have been deposited in 2018 with the nationwide Biobank of the National Institute for Health and Welfare, which has a license to operate from the National Supervisory Authority for Welfare and Health. Biobanked data are publicly available for use by qualified researchers worldwide following a standardized application procedure. Access to coded data is managed and controlled by the Biobank in accordance with the limitations of the consent forms signed by participants and Finnish law. A unique feature of the Biobank is the ability to link study records to Finland’s national health and social welfare registries, which includes extensive information regarding medical diagnoses and receipt of social services (e.g., disability claims). For details on accessing the biobanked data on the twins, see https://thl.fi/en/web/thl-biobank/for-researchers/application-process.

There are also hospital-based biobanks in Finland, where the pathology archives from the past 30 years or more are accessible. We are in the process of linking the older FTC to pathology archives to identify tissue samples from twin-pairs in which one or both twins have selected cancers, such as prostate, breast, colorectal and lung cancer. This has now become easier as there is a single portal, called FINGENIOUS, run by the Finnish Biobank cooperative (https://finbb.fi/), by which one can reach all six Finnish hospital biobanks by completing just one feasibility and access request.

Many DNA samples have been collected during targeted studies of selected pairs, conducted from the early 1990s onward, and all available samples have now been genotyped, imputed and available for genomewide association studies (GWASs) and other analyses. Table 2 also gives numbers of persons with genomewide genotypes or methylation data. Please note that the methylation data have generally been collected in later years, and not necessarily at the same time as the latest surveys but rather at the time of the respective substudy during an on-site visit. We also have some information on opposite-sex twins, some of which is register-based (twins born in 1950–1957) and some through questionnaire contact in 1996–1997 for twins born in the 1940s. Across all three cohorts, we have DNA samples on more than 15,000 twins and family members, who have taken part in family studies. Notably, the opposite-sex twins and first-degree relatives in the Nicotine Addiction Genetics (NAG) family study are not tabled here.

Almost 2000 methylation arrays have been done overall with 526 in the old cohort from 442 individuals and 83 with repeated measures. We are currently undertaking more methylation array analyses, which will add some 400 more in the old FTC. The more detailed breakdown of Ns will be given in separate review papers for Finntwin12 and Finntwin16 in the forthcoming special issue of twin registries in this journal. In addition, there are transcriptomics, histone modifications, miRNA analyses and metabolomics in selected subsets of twins.

Hypertension Epigenetics

Both genetic and environmental factors contribute to the pathogenesis of essential hypertension (EH). The GWASs have shown that the associated single-nucleotide polymorphisms (SNPs) explain less than 6% of the BP variation in the population (Evangelou et al., Reference Evangelou, Warren, Mosen-Ansorena, Mifsud, Pazoki, Gao and Caulfield2018; Hoffmann et al., Reference Hoffmann, Ehret, Nandakumar, Ranatunga, Schaefer, Kwok and Risch2017; Wain et al., Reference Wain, Vaez, Jansen, Joehanes, van der Most, Erzurumluoglu and Ehret2017; Warren et al., Reference Warren, Evangelou, Cabrera, Gao, Ren, Mifsud and Wain2017). Emerging evidence suggests that epigenetic variants play an important role in the pathogenesis of EH (Kato et al., Reference Kato, Loh, Takeuchi, Verweij, Wang, Zhang and Chambers2015; Wang et al., Reference Wang, Falkner, Zhu, Shi, Su, Xu and Snieder2013). Altogether, 199 twin-pairs and four singletons (age 55–69 years; 60.7% monozygous) from the EH-Epi substudy were included in epigenome and transcriptomewide analysis of hypertension. We identified 13 novel CpG methylation sites (article in preparation) and one new gene, MOK with its expression associated with BP, and replicated several previously identified signals in both gene expression and DNA methylation level. Results from bivariate twin modeling further provided new insights into the genetic and environmental sources of BP-related DNA methylation and gene expression signatures (Huang et al., Reference Huang, Ollikainen, Sipilä, Mustelin, Wang, Su and Wang2018).

Smoking and Mental Health Disorders

Depression and smoking are strongly associated, but the causal nature of their relationship is debated. Depression is a major risk factor for suicide and suicidality, but the role of smoking in suicide is unclear. We have examined longitudinal associations between smoking, depression and suicidality using the data from the older FTC.

We studied smoking behavior in 1975–1981 as a predictor of depression in 1990 (The Beck Depression Inventory, BDI score > 9) among 4164 men and 4934 women. Among the men, persistent smoking and smoking in 1975 but quitting by 1981 were associated with a higher risk of depression, while among the women only the quitters had an elevated risk. When family and genetic background were controlled, smoking remained a predictor of depression. Genetic modeling among men suggested a modest correlation (r _g = .25) between the genetic components of smoking and depression (Korhonen et al., Reference Korhonen, Broms, Varjonen, Romanov, Koskenvuo, Kinnunen and Kaprio2007).

We further examined which depression dimensions are linked to smoking patterns (4980 men and 5997 women). Longitudinal cigarette smoking patterns in 1975–1981 included multiple categories describing consistency and change. The three BDI dimensions were negative attitudes toward self (NATS), performance impairment and weight loss assessed in 1990. Persistent smoking and inconsistent former smoking predicted all the depression dimensions, although such associations with the NATS dimension were independent among men only. Long-term abstinence (consistent former smoking) did not predict risk for any depression dimension. Controlling for familial confounding, the association of persistent smoking with later NATS was replicated within discordant twin-pairs, suggesting that smoking is most robustly associated with the ‘core dimension’ of depression, that is, NATS (Korhonen et al., Reference Korhonen, Koivumaa-Honkanen, Varjonen, Broms, Koskenvuo and Kaprio2011).

While these previous studies were based on self-reported depressive symptoms, the last study investigated whether smoking predicts register-based outcome, that is, antidepressant medication (n = 10,652). A questionnaire was administered to Finnish adult twins in 1990. Antidepressant prescription data during 1995–2004 were obtained from the register of the Finnish Social Insurance Institution and linked to the survey data. Cox Proportional Hazard Models among 1075 cases and 9577 controls assessed the risk for depression in the cohort, whereas within-pair comparisons of smoking twins with their nonsmoking co-twins controlled for shared familial influences. Based on the analysis among those without baseline depression, heavy daily smokers had a significantly elevated likelihood for antidepressant prescription when adjusted for all confounders. A similar analysis using pairs discordant for antidepressant medication confirmed that daily smoking twins had a higher likelihood for antidepressant prescriptions compared with their nonsmoking co-twins (Korhonen et al., Reference Korhonen, Ranjit, Tuulio-Henriksson and Kaprio2017).

Within the older FTC, we investigated the relationship between smoking and suicide among 16,282 twin-pairs born before 1958. Smoking status and dose, marital status, employment and socioeconomic status and indicators of psychiatric and somatic illness were assessed at both 1975 and 1981 surveys. Emergent psychiatric and medical illness and vital status, including suicide determined by forensic autopsy, were evaluated over 35-year follow-up through national registries. The association between smoking and suicide was determined in competing risks hazard models. In twin-pairs discordant for smoking and suicide, the prospective association between smoking and suicide was determined using a matched case-control design. Smokers had a higher cumulative suicide incidence than former or never smokers. Heavy smokers had significantly higher suicide risk than light smokers. Compared with never smokers, smokers, but not former smokers, had increased suicide risk, adjusting for depressive symptoms, alcohol and sedative-hypnotic use, and excluding those who developed a serious somatic or psychiatric illness. In twin-pairs discordant for smoking and suicide, suicide was more likely in smokers. Adults who smoked tobacco were more likely to die by suicide, with a large, dose-dependent effect. This effect remained after consideration of many known predictors of suicide and shared familial effects, consistent with the hypothesis that exposure to tobacco smoke increases the risk of suicide (Evins et al., Reference Evins, Korhonen, Kinnunen and Kaprio2017).

Sleep

Sleep is also a major determinant of wellbeing and mental health, and over the past decades we have reported on the genetic and environmental contributions to sleep length and quality, as well as many other sleep-related traits, such as snoring, parasomnias and sleep bruxism. Recently, we have examined longitudinal associations of sleep quality across the four surveys among twins born in 1945–1957, showing that the large majority of respondents sleep well or very well across a 35-year period, while very few are consistently poor sleepers, indicating that more or less short-term circumstances, such as life events or illnesses, probably lie behind episodes of poor sleep, and poor sleep is persistent only among very few individuals (Hublin et al., Reference Hublin, Lehtovirta, Partinen, Koskenvuo and Kaprio2018). A longitudinal analysis of sleep length is ongoing.

Genetics and Epigenetics of Smoking, Nicotine Metabolism and Nicotine Dependence

Smoking behavior is influenced by both genetic and environmental factors, and epigenetic factors may link these two. FTC twins who have participated in the NAG, Schizophrenia and EH-Epi studies have been used to investigate both genetics and epigenetics and their interaction in smoking-related traits. Data from the FinnTwin12 and FinnTwin16 cohorts have also contributed. Genetics has been assessed by SNPs, while epigenetics has been assessed by measuring DNA methylation using genomewide array technologies.

Using the FinnTwin12 and FinnTwin16 together with other Finnish data (N = 1518), we performed the first GWAS on the rate of nicotine metabolism (Loukola et al., Reference Loukola, Buchwald, Gupta, Palviainen, Hällfors, Tikkanen and Kaprio2015), which influences smoking behavior such as response to smoking cessation treatments (Lerman et al., Reference Lerman, Tyndale, Patterson, Wileyto, Shields, Pinto and Benowitz2006). We found an association locus on chromosome 19, harboring the CYP2A6 gene. Based on our results, we constructed a genetic risk score explaining ~30% of the variance observed in the nicotine metabolite ratio (NMR) to study the association of the NMR and smoking in the older FTC data and other Finnish data. The genetic risk score was positively associated with cigarettes smoked per day (CPD) and smoking cessation, suggesting that individuals with faster nicotine metabolism smoke more, but are more likely to quit. The latter association was no longer statistically significant after accounting for potential confounders, such as quitting due to adverse health consequences or having major depression (Loukola et al., Reference Loukola, Buchwald, Gupta, Palviainen, Hällfors, Tikkanen and Kaprio2015).

Subsequently, we performed GWAS meta-analyses of the NMR and five other smoking-related traits in 5185 Finnish, Canadian and Australian current smokers, including data from the FTC (Buchwald et al., Reference Buchwald, Chenoweth, Palviainen, Zhu, Benner, Gordon and Tyndale2019). This was the largest GWAS to date of the NMR, and the first of the nicotine intake biomarker, cotinine + 3hydroxycotinine, and the smoking intensity variable, cotinine/CPD. We replicated previous findings, such as the chromosome 19 association locus for the NMR, as well as revealed new association loci and shed light on the genetic overlap between the studied traits.

We have also contributed to multiple other meta-analyses of GWAS data on various smoking and alcohol phenotypes, notably the GSCAN consortium GWAS and exome analyses (Brazel et al., Reference Brazel, Jiang, Hughey, Turcot, Zhan, Gong and Vrieze2018; Erzurumluoglu et al., Reference Erzurumluoglu, Liu, Jackson, Barnes, Datta, Melbourne and Howson2019; Liu et al., Reference Liu, Jiang, Wedow, Li, Brazel, Chen and Vrieze2019), for nicotine dependence (Hancock et al., Reference Hancock, Guo, Reginsson, Gaddis, Lutz, Sherva and Johnson2018) and alcohol dependence (Walters et al, Reference Walters, Polimanti, Johnson, McClintick, Adams, Adkins and Agrawal2018). These analyses have been based on alcohol and smoking behavior data from the surveys, as well as from the structured psychiatric interviews conducted as part of the NAG study. Additionally, we used the FTC data in our targeted hypothesis-driven study of the neuregulin signaling pathway (Gupta et al., Reference Gupta, Qaiser, He, Hiekkalinna, Zheutlin, Therman and Loukola2017). Our neuregulin signaling pathway-targeted analyses resulted in identification of SNPs with regulatory potential and provide evidence for the involvement of the pathway in smoking behavior and comorbid diseases. We have also done a GWAS on detailed smoking phenotypes covering the time dependency of progression of a smoker from initiation, to becoming a regular smoker, to smoking cessation based on detailed phenotyping in the NAG study (He et al., Reference He, Pitkäniemi, Heikkilä, Chou, Madden, Korhonen and Loukola2016).

The epigenome has been screened by epigenomewide association studies to identify genomic regions that are affected by or may influence smoking behavior (Gupta et al., Reference Gupta, van Dongen, Fu, Abdellaoui, Tyndale, Velagapudi and Ollikainen2019). Genetic and epigenetic data have also been integrated to uncover the regulatory potential of those variants that associate with smoking behavior phenotypes. We have observed multiple cis and trans methylation quantitative trait loci at genomewide significant SNPs or CpG sites and showed that DNA methylation is a molecular mechanism mediating the effects of genotype on smoking behavior phenotypes at some genetic loci (Gupta et al., Reference Gupta, van Dongen, Fu, Abdellaoui, Tyndale, Velagapudi and Ollikainen2019; Loukola et al., Reference Loukola, Buchwald, Gupta, Palviainen, Hällfors, Tikkanen and Kaprio2015).

As is well known, self-reported smoking status is prone to misclassification due to under-reporting, while biomarkers such as cotinine can measure only recent exposure. It is also known that smoking strongly influences DNA methylation, with current, former and never smokers exhibiting different methylation profiles. To meet the needs for accurate smoking status estimation, we have developed a robust DNA methylation-based predictor, EpiSmokEr, providing an objective measure of smoking status (https://github.com/sailalithabollepalli/EpiSmokEr; Bollepalli et al., Reference Bollepalli, Korhonen, Kaprio, Ollikainen and Anders2019). DNA methylation profiles of 405 EH-Epi twins were used in this project as a test dataset.

Alcohol Use and Its Consequences Over Time

Using measures of alcohol use behaviors at each study wave, we investigated birth cohort effects on the quantity and heritability of alcohol use. As expected, twins in earlier birth cohorts were more often abstainers and consumed less alcohol. However, we also found the contribution of genetic and environmental factors to variation in alcohol consumption to vary by birth cohort: in earlier cohorts, environmental factors shared by co-twins were more important, whereas the heritability of alcohol use behaviors was higher in more recent cohorts (Virtanen et al., Reference Virtanen, Kaprio, Viken, Rose and Latvala2018). The FTC offers unusual opportunities to study long-term associations of alcohol use and test for familial/genetic confounds in comparisons to exposure discordant MZ twin-pairs. While earlier analyses of alcohol use discordant pairs had not shown differences in mortality, longer follow-up and more refined phenotypes show that even in MZ twin-pairs discordant for alcohol use, the heavier user has a higher risk of death (Sipilä et al., Reference Sipilä, Rose and Kaprio2015).

We have completed an interview study of twin-pairs discordant for alcohol use in the 1975 baseline. From all pairs born in 1951–1957 in which both co-twins completed the 2011–2012 questionnaire follow-up, 171 pairs, confirmed or assumed to be MZ, were selected for their drinking discordance in 1975 when aged 18–25 years. Questionnaire measures of drinking over the 37-year follow-up confirm substantial stability for both individual differences and intrapair differences for consumption and alcohol problems across the nearly four-decade follow-up. The signed differences in consumption patterns on which pairs were selected in 1975 remained significant (albeit attenuated in magnitude) in the 2011–2012 questionnaire. Semistructured telephone interviews were completed for 128 of these 171 drinking-discordant twin-pairs, and intrapair differences in symptom counts of alcohol use disorder significantly correlated with measures of psychological wellbeing, general health, self-esteem and satisfaction with life. A smaller subset of interviewed pairs from the FinnTwin16 study shows similar findings from earlier ages and over a shorter duration of follow-up; a report of the combined results is in preparation awaiting zygosity confirmation for some twin-pairs.

Leisure-Time Physical Activity, Weight Gain, Health and Mortality — A Prospective Follow-up in Twins

We have continued the analyses on whether LTPA, adjusted for genetic factors and childhood environment, protects against mortality and different adverse health outcomes. All questionnaires (1975, 1981, 1990 and 2011) included questions on physical activity, allowing for a total metabolic equivalent (MET) index calculation (volume) by multiplying general intensity and the average duration and frequency of physical activity during leisure time. Different Finnish registers were used to study mortality, type 2 diabetes and use of antidepressants (measured as number of purchases). Follow-up telephone interview with detailed questions on the continuation of LTPA, self-measured weight and waist circumference and occurrence of chronic disease was carried out in 2005 (mean age 58.5 years) for physical activity discordant twin-pairs (n = 146 twin-pairs).

During follow-up (1976–2004), 1082 type 2 diabetes cases were observed among 20,487 individuals. Among all individuals, participants in MET quintiles (Q) III–V had significantly decreased risk for type 2 diabetes compared with sedentary individuals (QI). The pairwise analyses showed that the body mass index (BMI)-adjusted hazard ratio (HR) for the active (QII–V) compared with sedentary (QI) co-twins at follow-up was decreased by about one-half (0.54) (Waller, Kaprio et al., Reference Waller, Kaprio, Lehtovirta, Silventoinen, Koskenvuo and Kujala2010). Physical activity volume over 15 years (1975–1990) was used as the predictor of subsequent use of antidepressants for 1995–2004 among 642 physical activity discordant twin-pairs (Waller et al., Reference Waller, Kaprio, Korhonen, Tuulio-Henriksson and Kujala2016). Altogether, 229 persons had used at least one prescribed antidepressant during the study period. Active co-twins had a 20% lower risk for using any quantity of antidepressants than their inactive co-twins; trends were similar for DZ (20% decrease) and MZ pairs (22% decrease). The lowest odds ratio (OR) (0.51, 0.26–0.98) was seen among MZ pairs after adjusting for BMI, smoking and drinking.

Among the 146 LTPA discordant pairs, 24 co-twins (16 inactive and 8 active) had died between 1981 and 2004. The active co-twins had a strongly reduced risk of all-cause mortality as social class-adjusted HR was significantly reduced (HR = 0.39; Waller, Kujala, Rantanen et al., Reference Waller, Kujala, Rantanen, Kauppinen, Silventoinen, Koskenvuo and Kaprio2010). This was not found among the small number of MZ pairs. Mortality follow-up was extended for a period of 1981–2013 when twin-pairs discordant for vigorous LTPA were studied (discordance in both 1975 and 1981). This dataset included 1074 activity discordant twin-pairs (of which 231 were MZ pairs). Among 278 activity and death discordant pairs, 180 inactive and 98 active co-twins had died during follow-up. The unadjusted pairwise HR for all pairs was significantly reduced for all pairs (0.64) and for DZ pairs (0.58) but not for MZ pairs (0.85; Karvinen et al., Reference Karvinen, Waller, Silvennoinen, Koch, Britton, Kaprio and Kujala2015). When this analysis was repeated in the baseline-healthy subgroup, the HRs were further attenuated for DZ pairs and for MZ pairs (Karvinen et al., Reference Karvinen, Waller, Silvennoinen, Koch, Britton, Kaprio and Kujala2015).

On the other hand, an increased mortality by 41% was observed among twins with persistent dyspnea (breathlessness) at baseline (dyspnea in 1975 and 1981) compared with asymptomatic individuals. The within-pair analyses of death discordant pairs showed a similar and significant increase of 47% for all twin-pairs, with a point estimate of 164% for healthy MZ pairs (Waller et al., Reference Waller, Kaprio and Kujala2014).

Among the 95 interviewed physical activity discordant twin-pairs, the risk of type 2 diabetes or glucose intolerance (OR = 0.09, p = .022) and incident elevated BP (OR = 0.46, p = .039) was lower among the active co-twins. The active co-twins were more satisfied with their life at follow-up (p = .047). In contrast, the active co-twins showed a tendency toward more sports-related injuries (OR = 1.9, p = .051; Waller, Kujala, Kaprio et al., Reference Waller, Kujala, Kaprio, Koskenvuo and Rantanen2010). Within the subgroup of 42 pairs discordant for LTPA over 30 years, mean weight gain from 1975 through 2005 was 5.4 kg less (p = .003) and waist circumference 8.4 cm smaller (p < .001) at follow-up among the active compared with inactive co-twins (Waller et al., Reference Waller, Kaprio and Kujala2008).

A smaller subgroup of twins participated in an accelerometer-measured physical activity study called MOBILETWIN, which is a subset of the MEMTWIN younger twins, as it included all available same-sex twin-pairs born in 1940–1944. Twins completed a telephone interview and answered the questionnaire, and altogether 787 individuals, 404 women (283 complete pairs; 119 MZ pairs, mean age 72.9 years) wore a hip triaxial accelerometer for at least 4 days in 2014–2016.

High baseline mean MET (mean from 1975, 1981 and 1990) index predicted higher amounts of objectively measured moderate-to-vigorous physical activity (MVPA), with an increase in regression model R ² of 6.9% after age- and sex-adjustment (p < .001) at follow-up among individuals (Waller, Vähä-Ypyä, Lindgren et al., Reference Waller, Vähä-Ypyä, Lindgren, Kaprio, Sievänen and Kujala2018). Self-reported fitness at follow-up explained moderately the variation in measured MVPA (age- and sex-adjusted R ² 20%). MZ pairs were more similar (r = .57) in objectively measured MVPA at follow-up compared to DZ pairs (r = .28) (Waller, Vähä-Ypyä, Törmäkangas et al., Reference Waller, Vähä-Ypyä, Törmäkangas, Hautasaari, Lindgren, Iso-Markku and Kujala2018). Using mediation analysis, it was shown that shared genetic factors account for a large share of the associations between midlife and later physical activity variables (r _G = .59; 95% CI [0.44, 0.79]; Waller, Vähä-Ypyä, Lindgren et al., Reference Waller, Vähä-Ypyä, Lindgren, Kaprio, Sievänen and Kujala2018). In later years, low objectively measured physical activity was associated with the occurrence of chronic diseases (Kujala et al., Reference Kujala, Hautasaari, Vähä-Ypyä, Waller, Lindgren, Iso-Markku and Sievänen2019).

Three sets of analyses have examined causes and consequences of long-term LTPA (>1.5 MET-h/day) and inactivity (≤1.5 MET-h/day) by comparing 35-year LTPA profiles among twins born in 1945–1957 who replied to all four surveys in 1975, 1981, 1990 and 2011. The associations of 15-year smoking profile, education attainment by 1981 and 35-year BMI change with physical inactivity were studied (Piirtola et al., Reference Piirtola, Kaprio, Kujala, Heikkilä, Koskenvuo, Svedberg and Ropponen2016; Piirtola, Kaprio, Silventoinen et al., Reference Piirtola, Kaprio, Silventoinen, Svedberg, Korhonen and Ropponen2017; Piirtola, Kaprio, Waller et al., Reference Piirtola, Kaprio, Waller, Heikkilä, Koskenvuo, Svedberg and Ropponen2017). We found that long-term persistent daily smoking increased the likelihood of either remaining or becoming physically inactive, whereas quitting smoking decreased this likelihood (Piirtola, Kaprio, Waller et al., Reference Piirtola, Kaprio, Waller, Heikkilä, Koskenvuo, Svedberg and Ropponen2017). In analyzing long-term smoking discordant twin-pairs, any kind of smoking history was associated with physical inactivity. Further confirmation came from the observation that smoking at midlife also predicted less objective measured MVPA at follow-up in pairwise analysis among MZ twin-pairs (Waller, Vähä-Ypyä, Törmäkangas et al., Reference Waller, Vähä-Ypyä, Törmäkangas, Hautasaari, Lindgren, Iso-Markku and Kujala2018). Regarding another major risk factor, we found that higher education protected from physical inactivity because each additional year of education was associated with lower likelihood for inactivity both in the cross-sectional and longitudinal analyses (Piirtola et al., Reference Piirtola, Kaprio, Kujala, Heikkilä, Koskenvuo, Svedberg and Ropponen2016). These associations held true also for within-pair analyses, indicating an independent relationship between education attainment and LTPA. Persistent physical activity was protective against weight gain, also for within-pair analyses, though somewhat attenuated coefficients indicating confounding of familiar factors (Piirtola, Kaprio, Waller et al., Reference Piirtola, Kaprio, Waller, Heikkilä, Koskenvuo, Svedberg and Ropponen2017).

In summary, physical activity helps in maintaining overall health by decreasing the rate of weight gain, lowering waist circumference, reducing the risk for clinical type 2 diabetes and use of antidepressants. However, genetic factors may play a bigger role in explaining the associations between physical activity and mortality, as some of the findings were clearer among the DZ than MZ twin-pairs discordant for LTPA.

Causes and Consequences of Physical Activity: The TWINACTIVE Study

Sixteen twin-pairs (of whom seven were MZ pairs) with a persistent discordance in their LTPA levels were comprehensively identified from the FTC (Leskinen, Waller et al., Reference Leskinen, Waller, Mutikainen, Aaltonen, Ronkainen, Alen and Kujala2009). The discordance between the co-twins of the twin-pairs was initially defined from baseline questionnaires (1975 and 1981). A follow-up interview (from 1980 to 2005) was set to assess discordance in LTPA thereafter. The activity level assessment was based on calculation of the volume of activity, the MET index, from the self-reported physical activity.

Among the 16 pairs, the same twin had higher activity MET index than his or her co-twin throughout the follow-up period up to year 2007 when we conducted the multidimensional health-related measurements. The co-twins did not differ in other lifestyle factors, such as smoking, alcohol or work-related physical activity habits, either at baseline assessment (1975) or at the end of the follow-up (2007). The mean intrapair difference in LTPA MET indexes between the co-twins for the follow-up period was extensive, 8.8 MET h/day. As a result, the physically active co-twins had higher peak exercise capacity (9.3 METs vs. 7.5 METs, p < .001) and 20% higher knee extensor force compared to their inactive co-twins (Leskinen et al., Reference Leskinen, Waller, Mutikainen, Aaltonen, Ronkainen, Alen and Kujala2009, Reference Leskinen, Sipilä, Kaprio, Kainulainen, Alen and Kujala2013).

At the end of the follow-up, the inactive twins had a higher overall body fat content, 1.5 times more visceral and extramyocellular thigh fat and a significantly higher liver fat score compared to their active co-twins (Leskinen, Sipilä et al., Reference Leskinen, Sipilä, Alen, Cheng, Pietiläinen, Usenius and Kujala2009). These differences were also seen among the MZ co-twins only (Leskinen & Kujala, Reference Leskinen and Kujala2015). In skeletal muscle tissue, gene expression of the central pathways of energy metabolism and supportive metabolic pathways, especially those genes related to the processes of oxidative phosphorylation and branched-chain amino acid (BCAA) degradation, were upregulated among the physically active compared to inactive co-twins (Leskinen et al., Reference Leskinen, Rinnankoski-Tuikka, Rintala, Seppänen-Laakso, Pöllänen, Alen and Kujala2010). In subcutaneous fat tissue samples, the upregulated pathways in active compared to inactive co-twins included BCAA degradation and polyunsaturated fatty acid synthesis.

TWINACTIVE study shows that physically active lifestyle maintains aerobic capacity and muscle strength. Regular physical activity versus inactivity lowers the rate of ectopic fat accumulation even when the genetic background is considered. These findings, among our metabolomics finding (Kujala et al., Reference Kujala, Mäkinen, Heinonen, Soininen, Kangas, Leskinen and Ala-Korpela2013), indicate lower cardio-metabolic risk among persistently active versus inactive twins.

Predictors and Early Detection of Old Age Cognitive Impairment and Dementia — The MEMTWIN and ADTWIN Studies

Genetically informed studies of dementia risk factors with assessment of future cognition or dementia are rare; for that reason, the twin studies in our older twins were initiated as described above. On the one hand, questionnaire and register-based information were used (the MEMTWIN study); on the other hand, detailed assessments among selected pairs using neuropsychology and neuroimaging were done (ADTWIN study).

Future Prospects — Cognitive Aging in the Oldest-Old

Given that the older FTC included all Finnish same-sex twin-pairs born before 1958, many of these individuals have reached the age of 90 years. People who are 90 years or older are commonly referred to as the oldest-old, and they are the fastest growing population segment in many countries, including Finland. A big challenge is that we know very little about cognitive and brain aging in the oldest-old. Still, the prevalence and incidence of dementia are highest in this age group (Brumback-Peltz et al., Reference Brumback-Peltz, Balasubramanian, Corrada and Kawas2011). On the other hand, some of the oldest-old maintain their cognitive functioning and do have good memory even when compared to younger age groups.

The effects of risk factors of dementia may differ as a function of age. In addition to low education attainment, middle-age cardiovascular risk factors are the most robust risk factors of dementia. However, it is not clear whether the risk factor-dementia associations are similar across old age. For example, hypertension with onset after 80 years of age was related to decreased risk of dementia in the oldest-old (Corrada et al., Reference Corrada, Hayden, Paganini-Hill, Bullain, DeMoss, Aguirre and Kawas2017). Likewise, metabolic syndrome, a cluster of risk factors (hypertension, obesity, dyslipidemia and impaired glucose metabolism) of cardiovascular disease, has been reported to be associated with attenuated cognitive decline in the oldest-old (van den Berg et al., Reference van den Berg, Biessels, de Craen, Gussekloo and Westendorp2007). Many studies in the oldest-old have a short follow-up time, which may lead to reverse causation when investigating the risk factors after the disease process has started. Results based on several cohorts, including the older FTC, suggest that higher BMI is associated with higher dementia risk when the follow-up is over 20 years, whereas no association is observed with a follow-up of 10–20 years and there is an inverse association between BMI and dementia with less than 10 years of follow-up (Kivimäki et al., Reference Kivimäki, Luukkonen, Batty, Nyberg, Alfredsson, Fransson and Jokela2018). Genetic risk of dementia may also differ as a function of age: having an APOE ϵ4-allele confers risk for AD in people ≥65 years, but this effect seems to be attenuated among the oldest-old (Corrada et al., Reference Corrada, Paganini-Hill, Berlau and Kawas2013; Valerio et al., Reference Valerio, Raventos, Schmeidler, Beeri, Villalobos, Bolanos-Palmieri and Silverman2014). Taken together, the middle-age factors that predict dementia or successful cognitive aging in the oldest-old are not well known. Survival bias in studies of ≥90 year olds may result in having participants with high levels of resilience/cognitive reserve, as indexed by protective factors such as higher general cognitive ability and higher education or protective genetic background.

To learn more about the risk factors for dementia in the oldest-old, we are going to administer the telephone interview protocol to collect cognitive data from the FTC participants who are 90 years or older. Every year, more than 100 individuals are expected to reach the age of 90 years. These data together with the questionnaire data from 1975 will allow investigating the longitudinal association with a follow-up time of 45 years or more. Moreover, we can utilize the longitudinal FTC questionnaire data from the 1981 and 1990 questionnaires. Establishing a new database in 90-year-old twins will also make it possible to conduct more in-depth laboratory assessment, for example, by using the ADTWIN protocol of studying twin-pairs discordant for cognitive functioning. Also, registry-based data in the FTC can be used to study whether the risk factor-dementia associations are similar across the old age.

Finnish Twin Study on Aging

The purpose of the Finnish Twin Study on Aging (FITSA) is to study the genetic and environmental effects on the functional phenotypes that are relevant for the disablement process among older women. As FITSA includes three study waves, it is a longitudinal study embedded within the larger cohort.

The FITSA participants were recruited from the older FTC, including 1260 female twin-pairs born in 1924–1937 (Kaprio & Koskenvuo, Reference Kaprio and Koskenvuo2002; Kaprio et al., Reference Kaprio, Sarna, Koskenvuo and Rantasalo1978; Sarna et al., Reference Sarna, Kaprio, Sistonen and Koskenvuo1978). Of those women who were 63–76 years old in 2000, and who had participated in the FTC surveys in 1975 and 1981, invitations to take part in the FITSA were sent to 414 female twin-pairs living in Finland and 103 MZ and 114 DZ twin-pairs participated. The day-long protocol included the following assessments: biosampling, clinical measurements, anthropometrics, body composition, physical performance and sensory functions. Information about health habits and outcomes such as functional abilities and disabilities were collected in a questionnaire. Surveillance data on falls were collected for 1 year after the laboratory examination using a calendar method. Mortality has been followed up from the national registers.

First follow-up measurements of FITSA were conducted after 3 years, in 2003–2004. Altogether, 419 women participated in the follow-up study, of whom 313 participated in the laboratory measurements and filled in a questionnaire, while 106 women responded solely to the postal questionnaire. The second follow-up round for FITSA took place in 2011–2012. Altogether, 344 (91%) of the 377 surviving participants answered a structured, postal questionnaire, which was supplemented with telephone interviews when needed.

There are multiple dissertations and publications about FITSA pheno- and genotypes; for example, according to quantitative genetic analyses, muscle strength, power and walking speed have a genetic effect in common (Tiainen et al., Reference Tiainen, Pajala, Sipilä, Kaprio, Koskenvuo, Alén and Rantanen2007). Studies about falls have revealed that familial factors underlie the risk of falling (Pajala et al., Reference Pajala, Era, Koskenvuo, Kaprio, Viljanen and Rantanen2006) and that falls, as well as fear of falling, are associated with vision and hearing acuity (Kulmala et al., Reference Kulmala, Viljanen, Sipilä, Pajala, Pärssinen, Kauppinen and Rantanen2009; Viljanen et al., Reference Viljanen, Kaprio, Pyykkö, Sorri, Pajala, Kauppinen and Rantanen2009, Reference Viljanen, Kulmala, Rantakokko, Koskenvuo, Kaprio and Rantanen2013). Hearing and vision are both largely accounted for by genetic differences between individuals (Pärssinen et al., Reference Pärssinen, Kauppinen, Kaprio, Koskenvuo and Rantanen2013, Reference Pärssinen, Kauppinen, Kaprio, Koskenvuo and Rantanen2015, Reference Pärssinen, Kauppinen, Kaprio and Rantanen2016, Reference Pärssinen, Kauppinen, Halekoh, Kaprio and Rantanen2019; Viljanen, Era et al., Reference Viljanen, Era, Kaprio, Pyykkö, Koskenvuo and Rantanen2007; Viljanen, Kaprio et al., Reference Viljanen, Kaprio, Pyykkö, Sorri, Kauppinen, Koskenvuo and Rantanen2007). The FITSA study has contributed new knowledge especially in the area of genetic and environmental effects on phenotypes relevant for sensory and physical functioning underlying the disablement process in old age.

Because FITSA participants have also taken part in the FTC surveys in years 1975 and 1981, it has been possible to do longitudinal analyses in a genetically controlled sample with over 35-year follow-up from persons’ middle age to old age. It has been discovered, for example, that extraverted personality at middle age protected from, while neuroticism increased, the risk of later life depressive symptoms (Pakkala et al., Reference Pakkala, Read, Kaprio, Koskenvuo, Kauppinen and Rantanen2010). Middle-age living habits had no effect on depressive symptoms at old age. According to genetic analyses, neuroticism at middle age and depressive symptoms at old age shared an additive genetic component, while such a common genetic component was not observed with extraversion (Pakkala et al., Reference Pakkala, Read, Kaprio, Koskenvuo, Kauppinen and Rantanen2010). We have also reported that during the 29-year follow-up, BMI increased by 17%. Midlife BMI was a significant predictor of mobility 29 years later. Genetic influences on BMI level and its rate of change accounted for 37% and 25%, respectively, of the genetic influences on mobility later in life. The corresponding environmental influences on BMI level and its rate of change were 35% and 22% (Ortega-Alonso et al., Reference Ortega-Alonso, Sipilä, Kujala, Kaprio and Rantanen2009). FITSA data have also been utilized in several genetic studies. For example, vision, hearing and forced vital capacity have been analyzed as a part of larger multisite consortia identifying several new associated loci/genes (Loth et al., Reference Loth, Artigas, Gharib, Wain, Franceschini, Koch and London2014; Vuckovic et al., Reference Vuckovic, Dawson, Scheffer, Rantanen, Morgan, Di Stazio and Girotto2015).

Estrogen and Sarcopenia: A Co-Twin Study With HRT Discordant Postmenopausal Women (The SAWES Study)

Hormonal aging has notable effects on the functions of different organs and body systems as well as on functional capacity and health among middle-aged and older people. The most notable and well-characterized age-induced change in hormonal systems is the decline in sex hormone production in women during the menopausal transition. Low sex hormone concentrations are among the possible key mechanisms for sarcopenia. Sarcopenia is characterized by progressive loss of skeletal muscle mass and strength with an increased risk of adverse outcomes such as physical disability, poor health and poor quality of life and death. The underlying biological mechanisms by which age-induced sex hormone deprivation affects muscle mass and function are largely unknown.

The Sarcopenia — Skeletal Muscle Adaptation to Postmenopausal Hypogonadism and Effects of Hormone Replacement Therapy and Physical Activity in Older Women: a Genetic and Molecular Biological Study on Estrogen-related Pathways (SAWEs) study was launched in 2007 with the aim of identifying the key genetic factors needed for maintaining proper muscle mass, composition and function after menopause and to understand their connection with estrogen status and physical activity.

An invitation to participate the SAWEs study was sent to all female MZ twin-pairs born in 1943–1952 (537 pairs). Of these, 114 pairs responded and 21 pairs identified themselves as being HRT discordant, that is, one sister was a current HRT user, while the other had never used HRT. Finally, after confirming zygosity, 15 MZ twin-pairs from those who met the inclusion criteria formed the study group. At time of the study, enrolled participants were 54–62 years old and their HRT discordance ranged from 2 to 16 years, with a mean (SD) duration of 6.9 (4.1) years. Of the users, five women used estradiol preparations, six used estradiol with progestogen preparation and four women used tibolone.

This co-twin control study found postmenopausal HRT use to be beneficial for skeletal muscle. In particular, the greater relative area of the thigh muscles and the lower relative fat area within the muscle compartment among the sisters on HRT than their sisters with no history of HRT were found (Ronkainen et al., Reference Ronkainen, Kovanen, Alén, Pöllänen, Palonen, Ankarberg-Lindgren and Sipilä2009). No significant differences between the sisters were found in maximal voluntary isometric force (Finni et al., Reference Finni, Noorkoiv, Pöllänen, Ronkainen, Alén, Kaprio and Sipilä2011; Ronkainen et al., Reference Ronkainen, Kovanen, Alén, Pöllänen, Palonen, Ankarberg-Lindgren and Sipilä2009). However, when the effect of the central nervous system (voluntary activation) was blocked using electrical stimulation to the peripheral nerve, the force output of the plantar flexors was 32% greater among the HRT users compared to nonusers (Finni et al., Reference Finni, Noorkoiv, Pöllänen, Ronkainen, Alén, Kaprio and Sipilä2011). Accordingly, the specific tension of single type I and IIa myofibres was significantly greater among the HRT users than nonusers (Qaisar et al., Reference Qaisar, Renaud, Hedstrom, Pöllänen, Ronkainen, Kaprio and Larsson2013). In fast and more complex voluntary neuromuscular functions, such as the vertical countermovement jump and maximal walking speed, the sisters on HRT performed better than the nonusing sisters (Ronkainen et al., Reference Ronkainen, Kovanen, Alén, Pöllänen, Palonen, Ankarberg-Lindgren and Sipilä2009).

Estrogens have been suggested to play a role in tendon abnormalities, injuries and tendon diameter. Achilles tendon thickness and cross-sectional area did not significantly differ between HRT using and nonusing twin sisters. However, in five physically active pairs, the sister receiving HRT had a smaller tendon cross-sectional area than their sisters (Finni et al., Reference Finni, Kovanen, Ronkainen, Pöllänen, Bashford, Kaprio and Sipilä2009). In order to investigate the influence of long-term HRT use on bone strength and its determinants, the peripheral computed tomography data on tibial shaft, distal tibia and distal radius acquired from the HRT discordant twin-pairs at the FITSA (9 pairs) and SAWEs (15 pairs) studies were combined (Mikkola et al., Reference Mikkola, Heinonen, Kovanen, Cheng, Kujala, Suominen and Sipilä2011). The study found no difference between HRT users and nonusers in the total bone cross-sectional area in any of the measured bone sites but instead showed HRT users to have higher bending strength of the tibial shaft to be accounted from the larger cortical area and higher cortical bone mineral density. Furthermore, the estimated compressive strength of the distal tibia and distal radius was also higher among the HRT users compared to their nonusing co-twins.

Several approaches have been taken in the SAWEs to investigate the underlying molecular mechanisms. The investigation of the global muscle gene expression pointed toward estrogenic regulation of cell structure, cell-matrix interactions, energy metabolism and utilization of nutrients pathways (Ronkainen et al., Reference Ronkainen, Pöllänen, Alén, Pitkänen, Puolakka, Kujala and Kovanen2010). Further studies showed that the whole-body level differences favoring lower adiposity, higher lean mass and better muscle performance among HRT users were also apparent in the adiposity markers measured from blood circulation (Ahtiainen et al., Reference Ahtiainen, Alén, Pöllänen, Pulkkinen, Ronkainen, Kujala and Kovanen2012) and in the pathways identified from muscle proteomic study (Laakkonen et al., Reference Laakkonen, Soliymani, Karvinen, Rintala, Kaprio, Kujala and Lalowski2017). The focused inspection of the intramuscular sex steroid hormones (including estradiol, testosterone, dihydrotestosterone and dehydroepiandrosterone) revealed each hormone to be an independent predictor of female muscle strength and power regardless of age, systemic estradiol and body adiposity (Pöllänen et al., Reference Pöllänen, Kangas, Horttanainen, Niskala, Kaprio, Butler-Browne and Kovanen2015).

Several microRNA-related studies have also been conducted using samples from SAWEs. MicroRNAs are small regulatory RNA molecules that function as repressors of the protein synthesis and thereby regulate the activity of multitude of biological processes. MicroRNAs may directly affect the cells they have been produced but they can also be transported to the blood circulation and function as endocrine regulatory elements. Within the SAWEs study, we were first to identify estradiol-responsive microRNAs from skeletal muscle (Olivieri et al., Reference Olivieri, Ahtiainen, Lazzarini, Pöllänen, Capri, Lorenzi and Procopio2014) and later expanded these studies to include adipose tissue (Kangas, Morsiani et al., Reference Kangas, Morsiani, Pizza, Lanzarini, Aukee, Kaprio and Capri2017) and circulating microRNAs (Kangas et al., Reference Kangas, Pöllänen, Rippo, Lanzarini, Prattichizzo, Niskala and Kovanen2014; Kangas, Törmäkangas et al., Reference Kangas, Törmäkangas, Fey, Pursiheimo, Miinalainen, Alen and Laakkonen2017). The muscle-expressed estradiol-responding microRNAs, miR-182, miR-223 and miR-142–3p, had common targets in the insulin/IGF-1 pathway, which is known to participate in the regulation of muscle mass (Olivieri et al., Reference Olivieri, Ahtiainen, Lazzarini, Pöllänen, Capri, Lorenzi and Procopio2014). The circulating miRs showing differential expression due to use of HRT were associated with apoptosis, inflammation and body metabolic pathways (Kangas et al., Reference Kangas, Pöllänen, Rippo, Lanzarini, Prattichizzo, Niskala and Kovanen2014; Kangas, Morsiani et al., Reference Kangas, Morsiani, Pizza, Lanzarini, Aukee, Kaprio and Capri2017). These results are in line with our previous results indicating that the use of HRT affects muscle signaling at transcriptome (Ronkainen et al., Reference Ronkainen, Pöllänen, Alén, Pitkänen, Puolakka, Kujala and Kovanen2010) and proteome (Laakkonen et al., Reference Laakkonen, Soliymani, Karvinen, Rintala, Kaprio, Kujala and Lalowski2017) levels, as well as white blood cell signaling at methylome and transcriptome levels (Bahl et al., Reference Bahl, Pöllänen, Ismail, Sipilä, Mikkola, Berglund and Ollikainen2015).

Taken together, results from the SAWEs study indicate that menopausal changes in circulating steroid hormones, especially in female sex steroid estradiol, have multilevel effects on skeletal muscle gene expression, which can be seen at the phenotype level as increased adiposity, decreased lean mass, poorer bone properties and reduced muscle performance. The use of HRT partially aids maintaining better body composition and muscle function.

Biological Aging Clocks and Physical Activity and Functioning

Biological clocks, of which the most promising have been telomere length and DNA methylation age, also known as the ‘Epigenetic clock’, present an exciting opportunity to track and help understand the human biological aging process (Horvath, Reference Horvath2013; Jylhävä et al., Reference Jylhävä, Pedersen and Hägg2017). We have investigated genetic and environmental predictors of telomere length and DNA methylation age, and how these clocks are associated with not only determinants of older age physical functioning and LTPA at the total body level but also in muscle tissue, which has a vital role for health and functioning in old age. The aim of these studies has been to increase knowledge about the biological mechanisms underlying the development of age-related functional limitations and diseases.

Telomere length measurements by quantitative polymerase chain reaction (PCR) have been performed at baseline to subsample FITSA (n = 386). For SAWEs twins, telomere length measurements have also been conducted from muscle samples with the Southern blot method (Sillanpää, Niskala et al., Reference Sillanpää, Niskala, Laakkonen, Ponsot, Alén, Kaprio and Sipilä2017). We collaborated in a large European meta-analysis, which showed that telomere length is related to lung function and lung diseases (Albrecht et al., Reference Albrecht, Sillanpää, Karrasch, Alves, Codd, Hovatta and Schulz2014). Following that, we showed that the associations between telomere length and lung phenotypes are explained by shared genetic factors underlying both phenotypes (Sillanpää, Sipilä et al., Reference Sillanpää, Sipilä, Törmäkangas, Kaprio and Rantanen2017). Our analysis also revealed that telomere length predicts a decline in physical functioning over an 11-year follow-up (Sillanpää et al., Reference Sillanpää, Törmäkangas, Rantanen, Kaprio and Sipilä2016) and that telomere length is causally associated with cognitive functions (Hägg et al., Reference Hägg, Zhan, Karlsson, Gerritsen, Ploner, van der Lee and Pedersen2017). These results indicate that telomeres are genetically regulated biomarkers that can be used to predict the development of disability on older age. However, as predictors of life or health span, their value is rather low and overall, associations between telomeres and physical functioning are weak. Telomeres measured from muscle tissue may be more closely related to, for example, age-related loss in muscle mass and strength than telomeres measured from white blood cells (Sillanpää, Niskala et al., Reference Sillanpää, Niskala, Laakkonen, Ponsot, Alén, Kaprio and Sipilä2017).

Recently, we have expanded our biological aging research into DNA methylation clocks. DNA methylation age estimates can be produced from data generated by DNA methylation arrays. Currently, about 2000 DNA samples from individuals in the twin cohorts have been analyzed by Illumina HumanMethylation 450k or EPIC BeadChip. In comparison to telomeres, DNA methylation age may be more closely associated with physical functioning phenotypes in older people (Sillanpää et al., Reference Sillanpää, Kaprio, Laakkonen, Vaara, Rantanen, Kovanen and Ollikainen2018). Quantitative genetic modeling in older and younger FTC revealed that the relative contribution of nonshared environmental factors was larger among older compared with younger twin-pairs (Sillanpää et al., Reference Sillanpää, Ollikainen, Kaprio, Wang, Leskinen, Kujala and Törmäkangas2019). We tested the hypothesis that long-term LTPA is one of the environmental factors that affect variation in DNAm age in older age. A co-twin control analysis with older same-sex twin-pairs who had persistent discordance in physical activity for 32 years according to reported/interviewed physical-activity data (TWINACTIVE) showed no differences in DNA methylation age among active and inactive co-twins (Sillanpää et al., Reference Sillanpää, Ollikainen, Kaprio, Wang, Leskinen, Kujala and Törmäkangas2019). The nonexistent or at most very minor associations between LTPA and DNA methylation age were confirmed also in cross-sectional analysis in a sample of FinnTwin16 twin-pairs.