Global Gene Expression Profiling and Transcription Factor Network Analysis of Cognitive Aging in Monozygotic Twins.

Cognitive aging is one of the major problems worldwide, especially as people get older. This study aimed to perform global gene expression profiling of cognitive function to identify associated genes and pathways and a novel transcriptional regulatory network analysis to identify important regulons. We performed single transcript analysis on 400 monozygotic twins using an assumption-free generalized correlation coefficient (GCC), linear mixed-effect model (LME) and kinship model and identified six probes (one significant at the standard FDR < 0.05 while the other results were suggestive with 0.18 ≤ FDR ≤ 0.28). We combined the GCC and linear model results to cover diverse patterns of relationships, and meaningful and novel genes like APOBEC3G, H6PD, SLC45A1, GRIN3B, and PDE4D were detected. Our exploratory study showed the downregulation of all these genes with increasing cognitive function or vice versa except the SLC45A1 gene, which was upregulated with increasing cognitive function. Linear models found only H6PD and SLC45A1, the other genes were captured by GCC. Significant functional pathways (FDR < 3.95e-10) such as focal adhesion, ribosome, cysteine and methionine metabolism, Huntington's disease, eukaryotic translation elongation, nervous system development, influenza infection, metabolism of RNA, and cell cycle were identified. A total of five regulons (FDR< 1.3e-4) were enriched in a transcriptional regulatory analysis in which CTCF and REST were activated and SP3, SRF, and XBP1 were repressed regulons. The genome-wide transcription analysis using both assumption-free GCC and linear models identified important genes and biological pathways implicated in cognitive performance, cognitive aging, and neurological diseases. Also, the regulatory network analysis revealed significant activated and repressed regulons on cognitive function.

Differential lncRNA expression profiling of cognitive function in middle and old aged monozygotic twins using generalized association analysis.

Cognitive impairment is the most prominent symptom in neurodegenerative disorders affecting quality of life and mortality. However, despite years of research, the molecular mechanism underlying the regulation of cognitive function and its impairment is poorly understood. This study aims to elucidate the role of long non-coding RNAs (lncRNAs) expression and lncRNA-mRNA interaction networks, by analyzing lncRNA expression in whole blood samples of 400 middle and old aged monozygotic twins in association with cognitive function using both linear models and a generalized correlation coefficient (GCC) to capture the diverse patterns of correlation. We detected 13 probes (p < 1e-03) displaying nonlinear and 7 probes (p < 1e-03) showing linear correlations. After combining the results, we identified 20 lncRNA probes with p < 1e-03. The top lncRNA probes were annotated to genes, along with the non-coding MALAT1, that play roles in neurodegenerative diseases. The top lncRNAs were linked to functional clusters including peptidyl-glycine modification, vascular smooth muscle cells, mitotic spindle organization and protein tyrosine phosphatase. In addition, mapping of the top significant lncRNAs to the lncRNA-mRNA interaction network detected significantly enriched biological pathways involving neuroactive ligand-receptor interaction, proteasome and chemokines. We show that GCC served as a complementary approach in detecting lncRNAs missed by the conventional linear models. A combination of GCC and linear models identified lncRNAs of diverse patterns of association enriched for GO biological and molecular functions meaningful in cognitive performance and cognitive decline. The novel lncRNA regulatory network further contributed to detect significant pathways implicated in cognition.

Genetic meta-analysis of twin birth weight shows high genetic correlation with singleton birth weight.

Birth weight (BW) is an important predictor of newborn survival and health and has associations with many adult health outcomes, including cardiometabolic disorders, autoimmune diseases and mental health. On average, twins have a lower BW than singletons as a result of a different pattern of fetal growth and shorter gestational duration. Therefore, investigations into the genetics of BW often exclude data from twins, leading to a reduction in sample size and remaining ambiguities concerning the genetic contribution to BW in twins. In this study, we carried out a genome-wide association meta-analysis of BW in 42 212 twin individuals and found a positive correlation of beta values (Pearson's r = 0.66, 95% confidence interval [CI]: 0.47-0.77) with 150 previously reported genome-wide significant variants for singleton BW. We identified strong positive genetic correlations between BW in twins and numerous anthropometric traits, most notably with BW in singletons (genetic correlation [rg] = 0.92, 95% CI: 0.66-1.18). Genetic correlations of BW in twins with a series of health-related traits closely resembled those previously observed for BW in singletons. Polygenic scores constructed from a genome-wide association study on BW in the UK Biobank demonstrated strong predictive power in a target sample of Dutch twins and singletons. Together, our results indicate that a similar genetic architecture underlies BW in twins and singletons and that future genome-wide studies might benefit from including data from large twin registers.

Novel DNA methylation marker discovery by assumption-free genome-wide association analysis of cognitive function in twins.

Privileged by rapid increase in available epigenomic data, epigenome-wide association studies (EWAS) are to make a profound contribution to understand the molecular mechanism of DNA methylation in cognitive aging. Current statistical methods used in EWAS are dominated by models based on multiple assumptions, for example, linear relationship between molecular profiles and phenotype, normal distribution for the methylation data and phenotype. In this study, we applied an assumption-free method, the generalized correlation coefficient (GCC), and compare it to linear models, namely the linear mixed model and kinship model. We use DNA methylation associated with a cognitive score in 400 and 206 twins as discovery and replication samples respectively. DNA methylation associated with cognitive function using GCC, linear mixed model, and kinship model, identified 65 CpGs (p < 1e-04) from discovery sample displaying both nonlinear and linear correlations. Replication analysis successfully replicated 9 of these top CpGs. When combining results of GCC and linear models to cover diverse patterns of relationships, we identified genes like KLHDC4, PAPSS2, and MRPS18B as well as pathways including focal adhesion, axon guidance, and some neurological signaling. Genomic region-based analysis found 15 methylated regions harboring 11 genes, with three verified in gene expression analysis, also the 11 genes were related to top functional clusters including neurohypophyseal hormone and maternal aggressive behaviors. The GCC approach detects valuable methylation sites missed by traditional linear models. A combination of methylation markers from GCC and linear models enriched biological pathways sensible in neurological function that could implicate cognitive performance and cognitive aging.

Genome-wide association analysis of cognitive function in Danish long-lived individuals.

Cognitive function is a substantially heritable trait related to numerous important life outcomes. Several genome-wide association studies of cognitive function have in recent years led to the identification of thousands of significantly associated loci and genes. Individuals included in these studies have rarely been nonagenarians and centenarians, and since cognitive function is an important component of quality of life for this rapidly expanding demographic group, there is a need to explore genetic factors associated with individual differences in cognitive function at advanced ages. In this study, we pursued this by performing a genome-wide association study of cognitive function in 490 long-lived Danes (age range 90.1-100.8 years). While no genome-wide significant SNPs were identified, suggestively significant SNPs (P < 1 × 10-5) were mapped to several interesting genes, including ZWINT, CELF2, and DNAH5, and the glutamate receptor genes GRID2 and GRM7. Additionally, results from a gene set over-representation analysis indicated potential roles of gene sets related to G protein-coupled receptor (GPCR) signaling, interaction between L1 and ankyrins, mitogen-activated protein kinase (MAPK) signaling, RNA degradation, and cell cycle. Larger studies are needed to shed further light on the possible importance of these suggestive genes and pathways in cognitive function in nonagenarians and centenarians.

Exome-Wide Association Study Identifies FN3KRP and PGP as New Candidate Longevity Genes.

Despite enormous research efforts, the genetic component of longevity has remained largely elusive. The investigation of common variants, mainly located in intronic or regulatory regions, has yielded only little new information on the heritability of the phenotype. Here, we performed a chip-based exome-wide association study investigating 62 488 common and rare coding variants in 1248 German long-lived individuals, including 599 centenarians and 6941 younger controls (age < 60 years). In a single-variant analysis, we observed an exome-wide significant association between rs1046896 in the gene fructosamine-3-kinase-related-protein (FN3KRP) and longevity. Noteworthy, we found the longevity allele C of rs1046896 to be associated with an increased FN3KRP expression in whole blood; a database look-up confirmed this effect for various other human tissues. A gene-based analysis, in which potential cumulative effects of common and rare variants were considered, yielded the gene phosphoglycolate phosphatase (PGP) as another potential longevity gene, though no single variant in PGP reached the discovery p-value (1 × 10E-04). Furthermore, we validated the previously reported longevity locus cyclin-dependent kinase inhibitor 2B antisense RNA 1 (CDKN2B-AS1). Replication of our results in a French longevity cohort was only successful for rs1063192 in CDKN2B-AS1. In conclusion, we identified 2 new potential candidate longevity genes, FN3KRP and PGP which may influence the phenotype through their role in metabolic processes, that is, the reverse glycation of proteins (FN3KRP) and the control of glycerol-3-phosphate levels (PGP).

Polygenic Risk Score Prediction for Endometriosis.

Endometriosis is a major health care challenge because many young women with endometriosis go undetected for an extended period, which may lead to pain sensitization. Clinical tools to better identify candidates for laparoscopy-guided diagnosis are urgently needed. Since endometriosis has a strong genetic component, there is a growing interest in using genetics as part of the clinical risk assessment. The aim of this work was to investigate the discriminative ability of a polygenic risk score (PRS) for endometriosis using three different cohorts: surgically confirmed cases from the Western Danish endometriosis referral Center (249 cases, 348 controls), cases identified from the Danish Twin Registry (DTR) based on ICD-10 codes from the National Patient Registry (140 cases, 316 controls), and replication analysis in the UK Biobank (2,967 cases, 256,222 controls). Patients with adenomyosis from the DTR (25 cases) and from the UK Biobank (1,883 cases) were included for comparison. The PRS was derived from 14 genetic variants identified in a published genome-wide association study with more than 17,000 cases. The PRS was associated with endometriosis in surgically confirmed cases [odds ratio (OR) = 1.59, p = 2.57× 10-7] and in cases from the DTR biobank (OR = 1.50, p = 0.0001). Combining the two Danish cohorts, each standard deviation increase in PRS was associated with endometriosis (OR = 1.57, p = 2.5× 10-11), as well as the major subtypes of endometriosis; ovarian (OR = 1.72, p = 6.7× 10-5), infiltrating (OR = 1.66, p = 2.7× 10-9), and peritoneal (OR = 1.51, p = 2.6 × 10-3). These findings were replicated in the UK Biobank with a much larger sample size (OR = 1.28, p < 2.2× 10-16). The PRS was not associated with adenomyosis, suggesting that adenomyosis is not driven by the same genetic risk variants as endometriosis. Our results suggest that a PRS captures an increased risk of all types of endometriosis rather than an increased risk for endometriosis in specific locations. Although the discriminative accuracy is not yet sufficient as a stand-alone clinical utility, our data demonstrate that genetics risk variants in form of a simple PRS may add significant new discriminatory value. We suggest that an endometriosis PRS in combination with classical clinical risk factors and symptoms could be an important step in developing an urgently needed endometriosis risk stratification tool.

Quickly home again: patients' experiences of early discharge after minor stroke.

BACKGROUND: Patients with minor stroke experience some of the same issues as patients experiencing stroke of increased severity such as fatigue, anxiety and cognitive symptoms. It is current practice that patients with minor stroke receive accelerated treatment and care, yet studies indicate that patients find it difficult to return to their everyday lives after being discharged. We aimed to explore how patients with minor stroke experience the transitional period from the hospital through the first 2-4 weeks after an accelerated care pathway with discharge within 72 hours after stroke onset. METHODS: A qualitative study consisting of semi-structured interviews with 11 patients experiencing first-time stroke 2-4 weeks after discharge. RESULTS: The patients struggled to identify themselves as having had a stroke. They strived to find a new everyday life, but were challenged by existential concerns, mental fatigue and the fear of having a stroke again. Unresolved questions and misunderstandings arose, and the patients expressed a need for health professionals to support them and discuss unclear issues after discharge. Patients searched for others with similar issues in order to find a new sense of self. CONCLUSION: Patients with minor stroke struggle with everyday life after discharge. There is a need for support after discharge from healthcare professionals with specialised knowledge of stroke. Patients also requested an opportunity to meet other patients with minor stroke.

Generalized correlation coefficient for genome-wide association analysis of cognitive ability in twins.

Despite a strong genetic background in cognitive function only a limited number of single nucleotide polymorphisms (SNPs) have been found in genome-wide association studies (GWASs). We hypothesize that this is partially due to mis-specified modeling concerning phenotype distribution as well as the relationship between SNP dosage and the level of the phenotype. To overcome these issues, we introduced an assumption-free method based on generalized correlation coefficient (GCC) in a GWAS of cognitive function in Danish and Chinese twins to compare its performance with traditional linear models. The GCC-based GWAS identified two significant SNPs in Danish samples (rs71419535, p = 1.47e-08; rs905838, p = 1.69e-08) and two significant SNPs in Chinese samples (rs2292999, p = 9.27e-10; rs17019635, p = 2.50e-09). In contrast, linear models failed to detect any genome-wide significant SNPs. The number of top significant genes overlapping between the two samples in the GCC-based GWAS was higher than when applying linear models. The GCC model identified significant genetic variants missed by conventional linear models, with more replicated genes and biological pathways related to cognitive function. Moreover, the GCC-based GWAS was robust in handling correlated samples like twin pairs. GCC is a useful statistical method for GWAS that complements traditional linear models for capturing genetic effects beyond the additive assumption.

Differential long noncoding RNA profiling of BMI in twins.

Aim: Many efforts have been deployed to identify genetic variants associated with BMI. Alternatively, we explore epigenetic contribution to BMI variation by focusing on long noncoding RNAs (lncRNAs) which represents a key layer of epigenetic control. Materials & methods: We analyzed lncRNA expression in whole blood of 229 monozygotic twin pairs in association with BMI using generalized estimating equations. Results & conclusion: Six lncRNA probes were identified as significant (false discovery rate <0.05), with BMI showing causal effects on the expression of the significant lncRNAs. Functional annotation of differential profiles identified Gene ontology biological processes including kidney development, regulations of lipid biosynthetic process, circadian rhythm, notch signaling, etc. Whole blood lncRNAs are significantly expressed in response to BMI variation.

A decade of epigenetic change in aging twins: Genetic and environmental contributions to longitudinal DNA methylation.

BACKGROUND: Epigenetic changes may result from the interplay of environmental exposures and genetic influences and contribute to differences in age-related disease, disability, and mortality risk. However, the etiologies contributing to stability and change in DNA methylation have rarely been examined longitudinally. METHODS: We considered DNA methylation in whole blood leukocyte DNA across a 10-year span in two samples of same-sex aging twins: (a) Swedish Adoption Twin Study of Aging (SATSA; N = 53 pairs, 53% female; 62.9 and 72.5 years, SD = 7.2 years); (b) Longitudinal Study of Aging Danish Twins (LSADT; N = 43 pairs, 72% female, 76.2 and 86.1 years, SD=1.8 years). Joint biometrical analyses were conducted on 358,836 methylation probes in common. Bivariate twin models were fitted, adjusting for age, sex, and country. RESULTS: Overall, results suggest genetic contributions to DNA methylation across 358,836 sites tended to be small and lessen across 10 years (broad heritability M = 23.8% and 18.0%) but contributed to stability across time while person-specific factors explained emergent influences across the decade. Aging-specific sites identified from prior EWAS and methylation age clocks were more heritable than background sites. The 5037 sites that showed the greatest heritable/familial-environmental influences (p < 1E-07) were enriched for immune and inflammation pathways while 2020 low stability sites showed enrichment in stress-related pathways. CONCLUSIONS: Across time, stability in methylation is primarily due to genetic contributions, while novel experiences and exposures contribute to methylation differences. Elevated genetic contributions at age-related methylation sites suggest that adaptions to aging and senescence may be differentially impacted by genetic background.

A Genome-Wide Integrative Association Study of DNA Methylation and Gene Expression Data and Later Life Cognitive Functioning in Monozygotic Twins.

Monozygotic twins are genetically identical but rarely phenotypically identical. Epigenetic and transcriptional variation could influence this phenotypic discordance. Investigation of intra-pair differences in molecular markers and a given phenotype in monozygotic twins controls most of the genetic contribution, enabling studies of the molecular features of the phenotype. This study aimed to identify genes associated with cognition in later life using integrated enrichment analyses of the results of blood-derived intra-pair epigenome-wide and transcriptome-wide association analyses of cognition in 452 middle-aged and old-aged monozygotic twins (56-80 years). Integrated analyses were performed with an unsupervised approach using KeyPathwayMiner, and a supervised approach using the KEGG and Reactome databases. The supervised approach identified several enriched gene sets, including "neuroactive ligand receptor interaction" (p-value = 1.62∗10-2), "Neurotrophin signaling" (p-value = 2.52∗10-3), "Alzheimer's disease" (p-value = 1.20∗10-2), and "long-term depression" (p-value = 1.62∗10-2). The unsupervised approach resulted in a 238 gene network, including the Alzheimer's disease gene APP (Amyloid Beta Precursor Protein) as an exception node, and several novel candidate genes. The strength of the unsupervised method is that it can reveal previously uncharacterized sub-pathways and detect interplay between biological processes, which remain undetected by the current supervised methods. In conclusion, this study identified several previously reported cognition genes and pathways and, additionally, puts forward novel candidates for further verification and validation.

Cohort Differences in the Associations of Selected Candidate Genes With Risk of All-Cause Mortality at Advanced Ages.

Considerable efforts have been made to identify the genetic basis of human longevity, with only limited progress. One important drawback of current genetic studies is the limited knowledge of gene-environment interaction. Using 2 cohorts of long-lived individuals born in 1905 and 1915 in Denmark, we performed survival analysis to estimate risk of mortality for major candidate genes of aging and longevity and their cohort effects. Through statistical modeling that combines individual genetic and survival information with cohort-specific survival data, we estimated the relative risks of mortality from ages 95 to 103 years associated with genetic variants in apolipoprotein E (APOE), forkhead box class O3a, clusterin, and phosphatidylinositol binding clathrin assembly protein. Our analysis estimated a decreased risk of carrying the APOE$\varepsilon $4 allele (change in risk = -0.403, 95% confidence interval (CI): -0.831, 0.021; P = 0.040) in men of the later cohort, although the allele itself was harmful to survival across sexes (relative risk = 1.161, 95% CI: 1.027, 1.345; P = 0.026). We also estimated a cohort effect of increased risk for the minor allele of rs3851179 in phosphatidylinositol binding clathrin assembly protein with borderline significance (change in risk = 0.165, 95% CI: -0.010, 0.331; P = 0.052) in women. Our estimated significant cohort effect on APOE$\varepsilon $4 is indicative of the interplay between the gene and the changing environment that modulates survival at extreme ages.

Inspiratory muscle strength and walking capacity in patients with COPD.

Purpose: It has been suggested that patients with inspiratory muscle weakness could benefit from specific inspiratory muscle training (IMT). We aimed to examine the frequency of patients with inspiratory muscle weakness in a Danish hospital-based outpatient pulmonary rehabilitation program, and to evaluate the association between inspiratory muscle strength and peripheral muscle strength and walking capacity. Methods: Maximal Inspiratory Pressure (MIP) was assessed in 97 patients with COPD (39 men, 58 women, mean age years 70 ± 9, forced expiratory volume in 1 s ((FEV1) = 35 ± 10% pred.). The impact of MIP on knee-extension strength, walking distance, and symptom burden was evaluated using multiple linear regression analyses. Results: The MIP of the patients with COPD was 63 (95% CI 59; 67) cmH2O and it was significantly reduced compared to gender and age-matched reference values 76 (95% CI 73; 79) cmH2O (p < 0.001). Seven patients (7.2%) were under the lower limit of normal. MIP was negatively correlated with increasing age, female gender, decreasing knee-extension strength and lower FEV1% pred. Walking distance was associated with knee-extension strength and it was not associated with MIP. Conclusion: Maximal inspiratory pressure was reduced in patients with COPD but only a few patients had a weak MIP. Whilst MIP was associated with leg muscle strength, it was not associated with walking distance or symptoms.

Global expression profiling of cognitive level and decline in middle-aged monozygotic twins.

Only few studies have investigated the genomewide transcriptome of normative cognitive aging. We therefore aimed at investigating blood gene expression patterns associated with cognitive aging using a population-based sample of 235 middle-aged monozygotic twin pairs with longitudinal data on cognitive function. This unique setup enabled examination of gene expression differences associated with individual and intrapair differences in cognitive level and change while controlling for underlying genetic variation and shared early environment. Overall, increased expression of several gene sets was found to strongly correlate with a lower cognitive level and cognitive decline. The most significantly correlated gene sets were related to protein metabolism, translation, RNA metabolism, infectious disease, and the immune system, which are all processes previously linked to transcription signatures of pathological and normal brain aging, and aging in blood. The results of our study thus suggest that gene expression patterns of cognitive level and decline in our sample mirror those seen in cognitively impaired individuals, which could point toward a more generic response to cognitive aging and aging in general.

A meta-analysis of genome-wide association studies identifies multiple longevity genes.

Human longevity is heritable, but genome-wide association (GWA) studies have had limited success. Here, we perform two meta-analyses of GWA studies of a rigorous longevity phenotype definition including 11,262/3484 cases surviving at or beyond the age corresponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose age at death or at last contact was at or below the age corresponding to the 60th survival percentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE) ε4) is associated with lower odds of surviving to the 90th and 99th percentile age, while rs7412 (ApoE ε2) shows the opposite. Moreover, rs7676745, located near GPR78, associates with lower odds of surviving to the 90th percentile age. Gene-level association analysis reveals a role for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation of the longevity GWA results with that of several disease-related phenotypes points to a shared genetic architecture between health and longevity.

DNA methylome profiling in identical twin pairs discordant for body mass index.

OBJECTIVE: Body mass index (BMI) serves as an important measurement of obesity and adiposity, which are highly correlated with cardiometabolic diseases. Although high heritability has been estimated, the identified genetic variants by genetic association studies only explain a small proportion of BMI variation. As an active effort for further exploring the molecular basis of BMI variation, large-scale epigenome-wide association studies have been conducted but with limited number of loci reported, perhaps due to poorly controlled confounding factors, including genetic factors. Being genetically identical, monozygotic twins discordant for BMI are ideal subjects for analyzing the epigenetic association between DNA methylation and BMI, providing perfect control on their genetic makeups largely responsible for BMI variation. SUBJECTS: We performed an epigenome-wide association study on BMI using 30 identical twin pairs (15 male and 15 female pairs) with age ranging from 39 to 72 years and degree of BMI discordance ranging from 3-7.5 kg/m2. Methylation data from whole blood samples were collected using the reduced representation bisulfite sequencing technique. RESULTS: After adjusting for blood cell composition and clinical variables, we identified 136 CpGs with p-value < 1e-4, 30 CpGs with p < 1e-05 but no CpGs reached genome-wide significance. Genomic region-based analysis found 11 differentially methylated regions harboring coding and non-coding genes some of which were validated by gene expression analysis on independent samples. CONCLUSIONS: Our DNA methylation sequencing analysis on identical twins provides new references for the epigenetic regulation on BMI and obesity.

Epigenome-wide exploratory study of monozygotic twins suggests differentially methylated regions to associate with hand grip strength.

Hand grip strength is a measure of muscular strength and is used to study age-related loss of physical capacity. In order to explore the biological mechanisms that influence hand grip strength variation, an epigenome-wide association study (EWAS) of hand grip strength in 672 middle-aged and elderly monozygotic twins (age 55-90 years) was performed, using both individual and twin pair level analyses, the latter controlling the influence of genetic variation. Moreover, as measurements of hand grip strength performed over 8 years were available in the elderly twins (age 73-90 at intake), a longitudinal EWAS was conducted for this subsample. No genome-wide significant CpG sites or pathways were found, however two of the suggestive top CpG sites were mapped to the COL6A1 and CACNA1B genes, known to be related to muscular dysfunction. By investigating genomic regions using the comb-p algorithm, several differentially methylated regions in regulatory domains were identified as significantly associated to hand grip strength, and pathway analyses of these regions revealed significant pathways related to the immune system, autoimmune disorders, including diabetes type 1 and viral myocarditis, as well as negative regulation of cell differentiation. The genes contributing to the immunological pathways were HLA-B, HLA-C, HLA-DMA, HLA-DPB1, MYH10, ERAP1 and IRF8, while the genes implicated in the negative regulation of cell differentiation were IRF8, CEBPD, ID2 and BRCA1. In conclusion, this exploratory study suggests hand grip strength to associate with differentially methylated regions enriched in immunological and cell differentiation pathways, and hence merits further investigations.