Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci.

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P < 10-9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P < 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Ageing-associated changes in the human DNA methylome: genomic locations and effects on gene expression.

BACKGROUND: Changes in DNA methylation are among the mechanisms contributing to the ageing process. We sought to identify ageing-associated DNA methylation changes at single-CpG-site resolution in blood leukocytes and to ensure that the observed changes were not due to differences in the proportions of leukocytes. The association between DNA methylation changes and gene expression levels was also investigated in the same individuals. RESULTS: We identified 8540 high-confidence ageing-associated CpG sites, 46% of which were hypermethylated in nonagenarians. The hypermethylation-associated genes belonged to a common category: they were predicted to be regulated by a common group of transcription factors and were enriched in a related set of GO terms and canonical pathways. Conversely, for the hypomethylation-associated genes only a limited set of GO terms and canonical pathways were identified. Among the 8540 CpG sites associated with ageing, methylation level of 377 sites was also associated with gene expression levels. These genes were enriched in GO terms and canonical pathways associated with immune system functions, particularly phagocytosis. CONCLUSIONS: We find that certain ageing-associated immune-system impairments may be mediated via changes in DNA methylation. The results also imply that ageing-associated hypo- and hypermethylation are distinct processes: hypermethylation could be caused by programmed changes, whereas hypomethylation could be the result of environmental and stochastic processes.

Molecular mechanisms associated with the strength of the anti-CMV response in nonagenarians.

BACKGROUND: Infection with human cytomegalovirus (CMV) affects the function and composition of the immune system during ageing. In addition to the presence of the pathogen, the strength of the immune response, as measured by the anti-CMV IgG titre, has a significant effect on age-related pathogenesis. High anti-CMV IgG titres have been associated with increased mortality and functional impairment in the elderly. In this study, we were interested in identifying the molecular mechanisms that are associated with the strength of the anti-CMV response by examining the gene expression profiles that are associated with the level of the anti-CMV IgG titre. RESULTS: The level of the anti-CMV IgG titre is associated with the expression level of 663 transcripts in nonagenarians. These transcripts and their corresponding pathways are, for the most part, associated with metabolic functions, cell development and proliferation and other basic cellular functions. However, no prominent associations with the immune system were found, and no associated transcripts were found in young controls. CONCLUSIONS: The lack of defence pathways associated with the strength of the anti-CMV response can indicate that the compromised immune system can no longer defend itself against the CMV infection. Our data imply that the association between high anti-CMV IgG titres and increased mortality and frailty is mediated by basic cellular processes.

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies.

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

Polygenic risk for schizophrenia, social dispositions, and pace of epigenetic aging: Results from the Young Finns Study.

Schizophrenia is often regarded as a disorder of premature aging. We investigated (a) whether polygenic risk for schizophrenia (PRSsch ) relates to pace of epigenetic aging and (b) whether personal dispositions toward active and emotionally close relationships protect against accelerated epigenetic aging in individuals with high PRSsch . The sample came from the population-based Young Finns Study (n = 1348). Epigenetic aging was measured with DNA methylation aging algorithms such as AgeAccelHannum , EEAAHannum , IEAAHannum , IEAAHorvath , AgeAccelHorvath , AgeAccelPheno , AgeAccelGrim , and DunedinPACE. A PRSsch was calculated using summary statistics from the most comprehensive genome-wide association study of schizophrenia to date. Social dispositions were assessed in terms of extraversion, sociability, reward dependence, cooperativeness, and attachment security. We found that PRSsch did not have a statistically significant effect on any studied indicator of epigenetic aging. Instead, PRSsch had a significant interaction with reward dependence (p = 0.001-0.004), cooperation (p = 0.009-0.020), extraversion (p = 0.019-0.041), sociability (p = 0.003-0.016), and attachment security (p = 0.007-0.014) in predicting AgeAccelHannum , EEAAHannum , or IEAAHannum . Specifically, participants with high PRSsch appeared to display accelerated epigenetic aging at higher (vs. lower) levels of extraversion, sociability, attachment security, reward dependence, and cooperativeness. A rather opposite pattern was evident for those with low PRSsch . No such interactions were evident when predicting the other indicators of epigenetic aging. In conclusion, against our hypothesis, frequent social interactions may relate to accelerated epigenetic aging in individuals at risk for psychosis. We speculate that this may be explained by social-cognitive impairments (perceiving social situations as overwhelming or excessively arousing) or ending up in less supportive or deviant social groups.

NMR metabolomic modeling of age and lifespan: A multicohort analysis.

Metabolomic age models have been proposed for the study of biological aging, however, they have not been widely validated. We aimed to assess the performance of newly developed and existing nuclear magnetic resonance spectroscopy (NMR) metabolomic age models for prediction of chronological age (CA), mortality, and age-related disease. Ninety-eight metabolic variables were measured in blood from nine UK and Finnish cohort studies (N ≈31,000 individuals, age range 24-86 years). We used nonlinear and penalized regression to model CA and time to all-cause mortality. We examined associations of four new and two previously published metabolomic age models, with aging risk factors and phenotypes. Within the UK Biobank (N ≈102,000), we tested prediction of CA, incident disease (cardiovascular disease (CVD), type-2 diabetes mellitus, cancer, dementia, and chronic obstructive pulmonary disease), and all-cause mortality. Seven-fold cross-validated Pearson's r between metabolomic age models and CA ranged between 0.47 and 0.65 in the training cohort set (mean absolute error: 8-9 years). Metabolomic age models, adjusted for CA, were associated with C-reactive protein, and inversely associated with glomerular filtration rate. Positively associated risk factors included obesity, diabetes, smoking, and physical inactivity. In UK Biobank, correlations of metabolomic age with CA were modest (r = 0.29-0.33), yet all metabolomic model scores predicted mortality (hazard ratios of 1.01 to 1.06/metabolomic age year) and CVD, after adjustment for CA. While metabolomic age models were only moderately associated with CA in an independent population, they provided additional prediction of morbidity and mortality over CA itself, suggesting their wider applicability.

The relationship of trait-like compassion with epigenetic aging: The population-based prospective Young Finns Study.

Introduction: Helping others within and beyond the family has been related to living a healthy and long life. Compassion is a prosocial personality trait characterized by concern for another person who is suffering and the motivation to help. The current study examines whether epigenetic aging is a potential biological mechanism that explains the link between prosociality and longevity. Methods: We used data from the Young Finns Study that follows six birth-cohorts from age 3-18 to 19-49. Trait-like compassion for others was measured with the Temperament and Character Inventory in the years 1997 and 2001. Epigenetic age acceleration and telomere length were measured with five DNA methylation (DNAm) indicators (DNAmAgeHorvath, IEAA_Hannum, EEAA_Hannum, DNAmPhenoAge, and DNAmTL) based on blood drawn in 2011. We controlled for sex, socioeconomic status in childhood and adulthood, and body-mass index. Results and discussion: An association between higher compassion in 1997 and a less accelerated DNAmPhenoAge, which builds on previous work on phenotypic aging, approached statistical significance in a sex-adjusted model (n = 1,030; b = -0.34; p = 0.050). Compassion in 1997 predicted less accelerated epigenetic aging over and above the control variables (n = 843; b = -0.47; p = 0.016). There was no relationship between compassion in 2001 (n = 1108/910) and any of the other four studied epigenetic aging indicators. High compassion for others might indeed influence whether an individual's biological age is lower than their chronological age. The conducted robustness checks partially support this conclusion, yet cannot rule out that there might be a broader prosocial trait behind the findings. The observed associations are interesting but should be interpreted as weak requiring replication.

NMR metabolomic modelling of age and lifespan: a multi-cohort analysis.

Metabolomic age models have been proposed for the study of biological aging, however they have not been widely validated. We aimed to assess the performance of newly developed and existing nuclear magnetic resonance spectroscopy (NMR) metabolomic age models for prediction of chronological age (CA), mortality, and age-related disease. 98 metabolic variables were measured in blood from nine UK and Finnish cohort studies (N ≈ 31,000 individuals, age range 24-86 years). We used non-linear and penalised regression to model CA and time to all-cause mortality. We examined associations of four new and two previously published metabolomic age models, with ageing risk factors and phenotypes. Within the UK Biobank (N≈ 102,000), we tested prediction of CA, incident disease (cardiovascular disease (CVD), type-2 diabetes mellitus, cancer, dementia, chronic obstructive pulmonary disease) and all-cause mortality. Cross-validated Pearson's r between metabolomic age models and CA ranged between 0.47-0.65 in the training set (mean absolute error: 8-9 years). Metabolomic age models, adjusted for CA, were associated with C-reactive protein, and inversely associated with glomerular filtration rate. Positively associated risk factors included obesity, diabetes, smoking, and physical inactivity. In UK Biobank, correlations of metabolomic age with chronological age were modest (r = 0.29-0.33), yet all metabolomic model scores predicted mortality (hazard ratios of 1.01 to 1.06 / metabolomic age year) and CVD, after adjustment for CA. While metabolomic age models were only moderately associated with CA in an independent population, they provided additional prediction of morbidity and mortality over CA itself, suggesting their wider applicability.

Androgen receptor overexpression alters binding dynamics of the receptor to chromatin and chromatin structure.

BACKGROUND: Castration-resistant prostate cancers (CRPCs) overexpress often androgen receptor (AR). Here, we investigated the effect of AR overexpression on the dynamics of AR loading and RNA polymerase II (RNA Pol II) recruitment to chromatin. Acetylation of histone 3 (AcH3) on lysines 9 and 14 (K9 and K14) was also studied. METHODS: We used an LNCaP-based AR overexpression cell line model that includes a control line and two sublines, LNCaP-ARmo and LNCaP-ARhi, which overexpress AR twofold to threefold and fourfold to fivefold, respectively. Cells were exposed to 1 or 100 nM of dihydrotestosterone (DHT). Chromatin immunoprecipitation (ChIP) on the promoters and enhancers of prostate specific antigen (PSA) and transmembrane protease, serine 2 (TMPRSS2) genes was performed. qRT-PCR was used to measure the levels of PSA and TMPRSS2 transcripts. RESULTS: Upon stimulation with 1 nM DHT, AR and RNA Pol II were recruited onto PSA and TMPRSS2 enhancer regions to a greater extent (P < 0.05) in AR-overexpressing cells compared to control cells. The difference in AR loading between the control and AR-overexpressing cells was abolished by a higher DHT concentration. The ratio of AcH3/H3 was increased in AR-overexpressing cells. The induction of transcription of PSA and TMPRSS2 occurred earlier in the AR-overexpressing cells. CONCLUSIONS: Our findings suggest that the levels of AR potentiate the recruitment of the AR, as well as components of the basic transcription machinery, to chromatin and affect the acetylation of histones in the presence of low levels of androgens. These changes result in enhanced gene transcription of AR target genes.

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia.

BACKGROUND: In humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals. DATA: We obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs. RESULTS: In human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data. CONCLUSIONS: Our results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

Methylation status of VTRNA2-1/nc886 is stable across populations, monozygotic twin pairs and in majority of tissues.

Aims & methods: The aim of this study was to characterize the methylation level of a polymorphically imprinted gene, VTRNA2-1/nc886, in human populations and somatic tissues.48 datasets, consisting of more than 30 tissues and >30,000 individuals, were used. Results: nc886 methylation status is associated with twin status and ethnic background, but the variation between populations is limited. Monozygotic twin pairs present concordant methylation, whereas ∼30% of dizygotic twin pairs present discordant methylation in the nc886 locus. The methylation levels of nc886 are uniform across somatic tissues, except in cerebellum and skeletal muscle. Conclusion: The nc886 imprint may be established in the oocyte, and, after implantation, the methylation status is stable, excluding a few specific tissues.

Aging-associated increase in indoleamine 2,3-dioxygenase (IDO) activity appears to be unrelated to the transcription of the IDO1 or IDO2 genes in peripheral blood mononuclear cells.

BACKGROUND: Old age is associated with increased levels of circulating pro-inflammatory cytokines, a phenomenon termed inflamm-aging. Elevated levels of pro-inflammatory cytokines have been associated with several age-associated diseases and with a shortened lifespan. Indoleamine 2,3-dioxygenase (IDO) has immunomodulatory properties and its activity is elevated in inflammation, autoimmune disorders and malignancies. We have previously shown that IDO activity is increased in nonagenarians compared to young individuals and that high IDO activity is associated with mortality at old age. FINDINGS: In this study our aim was to assess whether this difference in IDO activity in the plasma was due to the differential expression of either the IDO1 or IDO2 gene in peripheral blood mononuclear cells. Our results show that IDO1 and IDO2 are not differently expressed in nonagenarians compared to controls and that the expression of IDO genes is not associated with the level of IDO activity in the plasma. CONCLUSION: The level of IDO activity in the plasma is not regulated through the expression of IDO1 or IDO2 in the peripheral blood mononuclear cells.

Ageing-associated changes in DNA methylation in X and Y chromosomes.

BACKGROUND: Ageing displays clear sexual dimorphism, evident in both morbidity and mortality. Ageing is also associated with changes in DNA methylation, but very little focus has been on the sex chromosomes, potential biological contributors to the observed sexual dimorphism. Here, we sought to identify DNA methylation changes associated with ageing in the Y and X chromosomes, by utilizing datasets available in data repositories, comprising in total of 1240 males and 1191 females, aged 14-92 years. RESULTS: In total, we identified 46 age-associated CpG sites in the male Y, 1327 age-associated CpG sites in the male X, and 325 age-associated CpG sites in the female X. The X chromosomal age-associated CpGs showed significant overlap between females and males, with 122 CpGs identified as age-associated in both sexes. Age-associated X chromosomal CpGs in both sexes were enriched in CpG islands and depleted from gene bodies and showed no strong trend towards hypermethylation nor hypomethylation. In contrast, the Y chromosomal age-associated CpGs were enriched in gene bodies, and showed a clear trend towards hypermethylation with age. CONCLUSIONS: Significant overlap in X chromosomal age-associated CpGs identified in males and females and their shared features suggest that despite the uneven chromosomal dosage, differences in ageing-associated DNA methylation changes in the X chromosome are unlikely to be a major contributor of sex dimorphism in ageing. While age-associated CpGs showed good replication across datasets in the present study, only a limited set of previously reported age-associated CpGs were replicated. One contributor to the limited overlap are differences in the age range of individuals included in each data set. Further study is needed to identify biologically significant age-associated CpGs in the sex chromosomes.

Reproductive history and blood cell DNA methylation later in life: the Young Finns Study.

BACKGROUND: Women with a history of complications of pregnancy, including hypertensive disorders, gestational diabetes or an infant fetal growth restriction or preterm birth, are at higher risk for cardiovascular disease later in life. We aimed to examine differences in maternal DNA methylation following pregnancy complications. METHODS: Data on women participating in the Young Finns study (n = 836) were linked to the national birth registry. DNA methylation in whole blood was assessed using the Infinium Methylation EPIC BeadChip. Epigenome-wide analysis was conducted on differential CpG methylation at 850 K sites. Reproductive history was also modeled as a predictor of four epigenetic age indices. RESULTS: Fourteen significant differentially methylated sites were found associated with both history of pre-eclampsia and overall hypertensive disorders of pregnancy. No associations were found between reproductive history and any epigenetic age acceleration measure. CONCLUSIONS: Differences in epigenetic methylation profiles could represent pre-existing risk factors, or changes that occurred as a result of experiencing these complications.

Adulthood blood levels of hsa-miR-29b-3p associate with preterm birth and adult metabolic and cognitive health.

Preterm birth (PTB) is associated with increased risk of type 2 diabetes and neurocognitive impairment later in life. We analyzed for the first time the associations of PTB with blood miRNA levels in adulthood. We also investigated the relationship of PTB associated miRNAs and adulthood phenotypes previously linked with premature birth. Blood MicroRNA profiling, genome-wide gene expression analysis, computer-based cognitive testing battery (CANTAB) and serum NMR metabolomics were performed for Young Finns Study subjects (aged 34-49 years, full-term n = 682, preterm n = 84). Preterm birth (vs. full-term) was associated with adulthood levels of hsa-miR-29b-3p in a fully adjusted regression model (p = 1.90 × 10-4, FDR = 0.046). The levels of hsa-miR-29b-3p were down-regulated in subjects with PTB with appropriate birthweight for gestational age (p = 0.002, fold change [FC] = - 1.20) and specifically in PTB subjects with small birthweight for gestational age (p = 0.095, FC = - 1.39) in comparison to individuals born full term. Hsa-miR-29b-3p levels correlated with the expressions of its target-mRNAs BCL11A and CS and the gene set analysis results indicated a target-mRNA driven association between hsa-miR-29b-3p levels and Alzheimer's disease, Parkinson's disease, Insulin signaling and Regulation of Actin Cytoskeleton pathway expression. The level of hsa-miR-29b-3p was directly associated with visual processing and sustained attention in CANTAB test and inversely associated with serum levels of VLDL subclass component and triglyceride levels. In conlcusion, adult blood levels of hsa-miR-29b-3p were lower in subjects born preterm. Hsa-miR-29b-3p associated with cognitive function and may be linked with adulthood morbidities in subjects born preterm, possibly through regulation of gene sets related to neurodegenerative diseases and insulin signaling as well as VLDL and triglyceride metabolism.

Modular genome-wide gene expression architecture shared by early traits of osteoporosis and atherosclerosis in the Young Finns Study.

We analysed whole blood genome-wide expression data to identify gene co-expression modules shared by early traits of osteoporosis and atherosclerosis. Gene expression was profiled for the Young Finns Study participants. Bone mineral density and content were measured as early traits of osteoporosis. Carotid and bulbus intima media thickness were measured as early traits of atherosclerosis. Joint association of the modules, identified with weighted co-expression analysis, with early traits of the diseases was tested with multivariate analysis. Among the six modules significantly correlated with early traits of both the diseases, two had significant (adjusted p-values (p.adj) < 0.05) and another two had suggestively significant (p.adj < 0.25) joint association with the two diseases after adjusting for age, sex, body mass index, smoking habit, alcohol consumption, and physical activity. The three most significant member genes from the significant modules were NOSIP, GXYLT2, and TRIM63 (p.adj ≤ 0.18). Genes in the modules were enriched with biological processes that have separately been found to be involved in either bone metabolism or atherosclerosis. The gene modules and their most significant member genes identified in this study support the osteoporosis-atherosclerosis comorbidity hypothesis and can provide new joint biomarkers for both diseases and their dual prevention.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs.

BACKGROUND: Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines. RESULTS: We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood. CONCLUSIONS: These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Ageing-associated changes in the expression of lncRNAs in human tissues reflect a transcriptional modulation in ageing pathways.

Ageing-associated changes in the protein coding transcriptome have been extensively characterised, but less attention has been paid to the non-coding portion of the human genome, especially to long non-coding RNAs (lncRNAs). Only a minority of known lncRNAs have been functionally characterised; however, a handful of these lncRNAs have already been linked to ageing-associated processes. To gain more information on the effects of ageing on lncRNAs, we identified from GTEx data lncRNAs that show ageing-associated expression patterns (age-lncRNAs) in 29 human tissues in 20-79-year-old individuals. The age-lncRNAs identified were highly tissue-specific, but the protein coding genes co-expressed with the age-lncRNAs and the functional categories associated with the age-lncRNAs showed significant overlap across tissues. Functions associated with the age-lncRNAs, including immune system processes and transcription, were similar to what has previously been reported for protein coding genes with ageing-associated expression pattern. As the tissue-specific age-lncRNAs were associated with shared functions across tissues, they may reflect the tissue-specific fine-tuning of the common ageing-associated processes. The present study can be utilised as a resource when selecting and prioritising lncRNAs for further functional analyses.

Epigenome-450K-wide methylation signatures of active cigarette smoking: The Young Finns Study.

Smoking as a major risk factor for morbidity affects numerous regulatory systems of the human body including DNA methylation. Most of the previous studies with genome-wide methylation data are based on conventional association analysis and earliest threshold-based gene set analysis that lacks sensitivity to be able to reveal all the relevant effects of smoking. The aim of the present study was to investigate the impact of active smoking on DNA methylation at three biological levels: 5'-C-phosphate-G-3' (CpG) sites, genes and functionally related genes (gene sets). Gene set analysis was done with mGSZ, a modern threshold-free method previously developed by us that utilizes all the genes in the experiment and their differential methylation scores. Application of such method in DNA methylation study is novel. Epigenome-wide methylation levels were profiled from Young Finns Study (YFS) participants' whole blood from 2011 follow-up using Illumina Infinium HumanMethylation450 BeadChips. We identified three novel smoking related CpG sites and replicated 57 of the previously identified ones. We found that smoking is associated with hypomethylation in shore (genomic regions 0-2 kilobases from CpG island). We identified smoking related methylation changes in 13 gene sets with false discovery rate (FDR) ≤ 0.05, among which is olfactory receptor activity, the flagship novel finding of the present study. Overall, we extended the current knowledge by identifying: (i) three novel smoking related CpG sites, (ii) similar effects as aging on average methylation in shore, and (iii) a novel finding that olfactory receptor activity pathway responds to tobacco smoke and toxin exposure through epigenetic mechanisms.

Aging-associated patterns in the expression of human endogenous retroviruses.

Human endogenous retroviruses (HERV) are relics of ancient retroviral infections in our genome. Most of them have lost their coding capacity, but proviral RNA or protein have been observed in several disease states (e.g. in inflammatory and autoimmune diseases and malignancies). However, their clinical significance as well as their mechanisms of action have still remained elusive. As human aging is associated with several biological characteristics of these diseases, we now analyzed the aging-associated expression of the individual proviruses of two HERV families, HERV-K (91 proviruses) and HERV-W (213 proviruses) using genome-wide RNA-sequencing (RNA-seq). RNA was purified from blood cells derived from healthy young individuals (n = 7) and from nonagenarians (n = 7). The data indicated that in the case of HERV-K (HML-2) 33 proviruses had a detectable expression but in only 3 of those the expression levels were significantly different between the young and old individuals. In the HERV-W family expression was observed in 45 loci and only in one case the young/old difference was significant. However, applying hierarchical clustering on the HERV expression data resulted in the formation of two distinct clusters, one containing the young individuals and another the nonagenarians. This suggests, that even though the aging-associated differences in the expression levels of individual proviruses are minor, there seems to be some underlying aging-related pattern. These data indicate that aging does not have a strong effect on the expression of individual HERV proviruses, but instead several proviruses are affected moderately, leading to age-dependent expression profiles.