Omixer: multivariate and reproducible sample randomization to proactively counter batch effects in omics studies

Genetic and environmental contributions to epigenetic aging across adolescence and young adulthood

BACKGROUND

Epigenetic aging estimators commonly track chronological and biological aging, quantifying its accumulation (i.e., epigenetic age acceleration) or speed (i.e., epigenetic aging pace). Their scores reflect a combination of inherent biological programming and the impact of environmental factors, which are suggested to vary at different life stages. The transition from adolescence to adulthood is an important period in this regard, marked by an increasing and, then, stabilizing epigenetic aging variance. Whether this pattern arises from environmental influences or genetic factors is still uncertain. This study delves into understanding the genetic and environmental contributions to variance in epigenetic aging across these developmental stages. Using twin modeling, we analyzed four estimators of epigenetic aging, namely Horvath Acceleration, PedBE Acceleration, GrimAge Acceleration, and DunedinPACE, based on saliva samples collected at two timepoints approximately 2.5 years apart from 976 twins of four birth cohorts (aged about 9.5, 15.5, 21.5, and 27.5 years at first and 12, 18, 24, and 30 years at second measurement occasion).

RESULTS

Half to two-thirds (50-68%) of the differences in epigenetic aging were due to unique environmental factors, indicating the role of life experiences and epigenetic drift, besides measurement error. The remaining variance was explained by genetic (Horvath Acceleration: 24%; GrimAge Acceleration: 32%; DunedinPACE: 47%) and shared environmental factors (Horvath Acceleration: 26%; PedBE Acceleration: 47%). The genetic and shared environmental factors represented the primary sources of stable differences in corresponding epigenetic aging estimators over 2.5 years. Age moderation analyses revealed that the variance due to individually unique environmental sources was smaller in younger than in older cohorts in epigenetic aging estimators trained on chronological age (Horvath Acceleration: 47-49%; PedBE Acceleration: 33-68%). The variance due to genetic contributions, in turn, potentially increased across age groups for epigenetic aging estimators trained in adult samples (Horvath Acceleration: 18-39%; GrimAge Acceleration: 24-43%; DunedinPACE: 42-57%).

CONCLUSIONS

Transition to adulthood is a period of the increasing variance in epigenetic aging. Both environmental and genetic factors contribute to this trend. The degree of environmental and genetic contributions can be partially explained by the design of epigenetic aging estimators.

Tissue-specific methylomic responses to a lifestyle intervention in older adults associate with metabolic and physiological health improvements

Across the lifespan, diet and physical activity profiles substantially influence immunometabolic health. DNA methylation, as a tissue-specific marker sensitive to behavioral change, may mediate these effects through modulation of transcription factor binding and subsequent gene expression. Despite this, few human studies have profiled DNA methylation and gene expression simultaneously in multiple tissues or examined how molecular levels react and interact in response to lifestyle changes. The Growing Old Together (GOTO) study is a 13-week lifestyle intervention in older adults, which imparted health benefits to participants. Here, we characterize the DNA methylation response to this intervention at over 750 thousand CpGs in muscle, adipose, and blood. Differentially methylated sites are enriched for active chromatin states, located close to relevant transcription factor binding sites, and associated with changing expression of insulin sensitivity genes and health parameters. In addition, measures of biological age are consistently reduced, with decreases in grimAge associated with observed health improvements. Taken together, our results identify responsive molecular markers and demonstrate their potential to measure progression and finetune treatment of age-related risks and diseases.

DNA methylation differences between cord and adult white blood cells reflect postnatal immune cell maturation

Epigenetic modifications such as DNA methylation are both cell type and developmental age specific. Here, we show that the immunological maturation of blood cell types influences DNA methylation changes from naive cord blood to fully functional adult blood. Lymphoid cells in adult blood showed more variability than in cord blood suggesting an antigen-dependent maturation of DNA methylation in lymphoid cells throughout the lifespan. Fewer DNA methylation changes between cord and adult blood were observed in myeloid cells, particularly in monocytes, which demonstrated the least number of DNA methylation changes between cord and adult blood. We also noted differences in epigenetic ages by immune cell types within the same individuals, specifically in cord blood where monocytes were epigenetically oldest compared to the other cell types. In addition, we provide a publicly available resource to the community as an R Shiny web application to interactively explore epigenetic patterns between naive cord white blood cells and fully functional adult white blood cells for six immune cell types.

A Comparison of Three Automated Nucleic Acid Extraction Systems for Human Stool Samples

Automated nucleic acid extractors are useful instruments for the high-throughput processing of bio-samples and are expected to improve research throughput in addition to decreased inter-sample variability inherent to manual processing. We evaluated three commercial nucleic acid extractors Bioer GenePure Pro (Bioer Technology, Hangzhou, China), Maxwell RSC 16 (Promega Corporation, Madison, WI, USA), and KingFisher Apex (ThermoFisher Scientific, Waltham, MA, USA) based on their DNA yield, DNA purity, and 16S rRNA gene amplicon results using both human fecal samples and a mock community (ZymoBIOMICS Microbial Community Standard (Zymo Research Corp., Irvine, CA, USA)). Bead-beating provided incremental yield to effectively lyse and extract DNA from stool samples compared to lysis buffer alone. Differential abundance analysis and comparison of prevalent bacterial species revealed a greater representation of Gram-positive bacteria in samples subjected to mechanical lysis, regardless of sample type. All three commercial extractors had differences in terms of yield, inter-sample variability, and subsequent sequencing readouts, which we subsequently share in the paper and believe are significant considerations for all researchers undertaking human fecal microbiota research.

Epigenetic clocks, sex markers and age-class diagnostics in three harvested large mammals

The development of epigenetic clocks, or the DNA methylation-based inference of age, is an emerging tool for ageing in free ranging populations. In this study, we developed epigenetic clocks for three species of large mammals that are the focus of extensive management throughout their range in North America: white-tailed deer, black bear and mountain goat. We quantified differential DNA methylation patterns at over 30,000 cytosine-guanine sites (CpGs) from tissue samples of all three species (black bear n = 49; white-tailed deer n = 47; mountain goat n = 45). We used a penalized regression model (elastic net) to build explanatory (black bear r = .95; white-tailed deer r = .99; mountain goat r = .97) and robust (black bear Median Absolute Error or MAE = 1.33; white-tailed deer MAE = 0.29; mountain goat MAE = 0.61) models of age or clocks. We also characterized individual CpG sites within each species that demonstrated clear differences in methylation levels between age classes and sex, which can be used to develop a suite of accessible diagnostic markers. This tool has the potential to contribute to wildlife monitoring by providing easily obtainable representations of age structure in managed populations.

DNA methylation patterns of FKBP5 regulatory regions in brain and blood of humanized mice and humans

Humanized mouse models can be used to explore human gene regulatory elements (REs), which frequently lie in non-coding and less conserved genomic regions. Epigenetic modifications of gene REs, also in the context of gene x environment interactions, have not yet been explored in humanized mouse models. We applied high-accuracy measurement of DNA methylation (DNAm) via targeted bisulfite sequencing (HAM-TBS) to investigate DNAm in three tissues/brain regions (blood, prefrontal cortex and hippocampus) of mice carrying the human FK506-binding protein 5 (FKBP5) gene, an important candidate gene associated with stress-related psychiatric disorders. We explored DNAm in three functional intronic glucocorticoid-responsive elements (at introns 2, 5, and 7) of FKBP5 at baseline, in cases of differing genotype (rs1360780 single nucleotide polymorphism), and following application of the synthetic glucocorticoid dexamethasone. We compared DNAm patterns in the humanized mouse (N = 58) to those in human peripheral blood (N = 447 and N = 89) and human postmortem brain prefrontal cortex (N = 86). Overall, DNAm patterns in the humanized mouse model seem to recapitulate DNAm patterns observed in human tissue. At baseline, this was to a higher extent in brain tissue. The animal model also recapitulated effects of dexamethasone on DNAm, especially in peripheral blood and to a lesser extent effects of genotype on DNAm. The humanized mouse model could thus assist in reverse translation of human findings in psychiatry that involve genetic and epigenetic regulation in non-coding elements.

Proteomic Evolution from Acute to Post-COVID-19 Conditions

Many COVID-19 survivors have post-COVID-19 conditions, and females are at a higher risk. We sought to determine (1) how protein levels change from acute to post-COVID-19 conditions, (2) whether females have a plasma protein signature different from that of males, and (3) which biological pathways are associated with COVID-19 when compared to restrictive lung disease. We measured protein levels in 74 patients on the day of admission and at 3 and 6 months after diagnosis. We determined protein concentrations by multiple reaction monitoring (MRM) using a panel of 269 heavy-labeled peptides. The predicted forced vital capacity (FVC) and diffusing capacity of the lungs for carbon monoxide (DLCO) were measured by routine pulmonary function testing. Proteins associated with six key lipid-related pathways increased from admission to 3 and 6 months; conversely, proteins related to innate immune responses and vasoconstriction-related proteins decreased. Multiple biological functions were regulated differentially between females and males. Concentrations of eight proteins were associated with FVC, %, and they together had c-statistics of 0.751 (CI:0.732-0.779); similarly, concentrations of five proteins had c-statistics of 0.707 (CI:0.676-0.737) for DLCO, %. Lipid biology may drive evolution from acute to post-COVID-19 conditions, while activation of innate immunity and vascular regulation pathways decreased over that period. (ProteomeXchange identifiers: PXD041762, PXD029437).

Injectiondesign: web service of plate design with optimized stratified block randomization for modern GC/LC-MS-based sample preparation

BACKGROUND

Plate design is a necessary and time-consuming operation for GC/LC-MS-based sample preparation. The implementation of the inter-batch balancing algorithm and the intra-batch randomization algorithm can have a significant impact on the final results. For researchers without programming skills, a stable and efficient online service for plate design is necessary.

RESULTS

Here we describe InjectionDesign, a free online plate design service focused on GC/LC-MS-based multi-omics experiment design. It offers the ability to separate the position design from the sequence design, making the output more compatible with the requirements of a modern mass spectrometer-based laboratory. In addition, it has implemented an optimized block randomization algorithm, which can be better applied to sample stratification with block randomization for an unbalanced distribution. It is easy to use, with built-in support for common instrument models and quick export to a worksheet.

CONCLUSIONS

InjectionDesign is an open-source project based on Java. Researchers can get the source code for the project from Github: https://github.com/CSi-Studio/InjectionDesign . A free web service is also provided: http://www.injection.design .

Epigenetic age acceleration as a biomarker for impaired cognitive abilities in adulthood following early life adversity and psychiatric disorders

Background

Early life adversity and psychiatric disorders are associated with earlier declines in neurocognitive abilities during adulthood. These declines may be preceded by changes in biological aging, specifically epigenetic age acceleration, providing an opportunity to uncover genome-wide biomarkers that identify individuals most likely to benefit from early screening and prevention.

Methods

Five unique epigenetic age acceleration clocks derived from peripheral blood were examined in relation to latent variables of general and speeded cognitive abilities across two independent cohorts: 1) the Female Growth and Development Study (FGDS; n = 86), a 30-year prospective cohort study of substantiated child sexual abuse and non-abused controls, and 2) the Biological Classification of Mental Disorders study (BeCOME; n = 313), an adult community cohort established based on psychiatric disorders.

Results

A faster pace of biological aging (DunedinPoAm) was associated with lower general cognitive abilities in both cohorts and slower speeded abilities in the BeCOME cohort. Acceleration in the Horvath clock was significantly associated with slower speeded abilities in the BeCOME cohort but not the FGDS. Acceleration in the Hannum clock and the GrimAge clock were not significantly associated with either cognitive ability. Accelerated PhenoAge was associated with slower speeded abilities in the FGDS but not the BeCOME cohort.

Conclusions

The present results suggest that epigenetic age acceleration has the potential to serve as a biomarker for neurocognitive decline in adults with a history of early life adversity or psychiatric disorders. Estimates of epigenetic aging may identify adults at risk of cognitive decline that could benefit from early neurocognitive screening.

Early-life exposure to cigarette smoke primes lung function and DNA methylation changes at Cyp1a1 upon exposure later in life

Prenatal and early-life exposure to cigarette smoke (CS) has repeatedly been shown to induce stable, long-term changes in DNA methylation (DNAm) in offspring. It has been hypothesized that these changes might be functionally related to the known outcomes of prenatal and early-life CS exposure, which include impaired lung development, altered lung function, and increased risk of asthma and wheeze. However, to date, few studies have examined DNAm changes induced by prenatal CS in tissues of the lung, and even fewer have attempted to examine the specific influences of prenatal versus early postnatal exposures. Here, we have established a mouse model of CS exposure which isolates the effects of prenatal and early postnatal CS exposures in early life. We have used this model to measure the effects of prenatal and/or postnatal CS exposures on lung function and immune cell infiltration as well as DNAm and expression of Cyp1a1, a candidate gene previously observed to demonstrate DNAm differences on CS exposure in humans. Our study revealed that exposure to CS prenatally and in the early postnatal period causes long-lasting differences in offspring lung function, gene expression, and lung Cyp1a1 DNAm, which wane over time but are reestablished on reexposure to CS in adulthood. This study creates a testable mouse model that can be used to investigate the effects of prenatal and early postnatal CS exposures and will contribute to the design of intervention strategies to mediate these detrimental effects.NEW & NOTEWORTHY Here, we isolated effects of prenatal from early postnatal cigarette smoke and showed that exposure to cigarette smoke early in life causes changes in offspring DNA methylation at Cyp1a1 that last through early adulthood but not into late adulthood. We also showed that smoking in adulthood reestablished these DNA methylation patterns at Cyp1a1, suggesting that a mechanism other than DNA methylation results in long-term memory associated with early-life cigarette smoke exposures at this gene.

Extensive evaluation of DNA methylation of functional elements in the murine Fkbp5 locus using high-accuracy DNA methylation measurement via targeted bisulfite sequencing

FKBP5 is an important stress-regulatory gene implicated in stress-related psychiatric diseases. Single nucleotide polymorphisms of the FKBP5 gene were shown to interact with early life stress to alter the glucocorticoid-related stress response and moderate disease risk. Demethylation of cytosine-phosphate-guanine-dinucleotides (CpGs) in regulatory glucocorticoid-responsive elements was suggested to be the mediating epigenetic mechanism for long-term stress effects, but studies on Fkbp5 DNA methylation (DNAm) in rodents are so far limited. We evaluated the applicability of high-accuracy DNA methylation measurement via targeted bisulfite sequencing (HAM-TBS), a next-generation sequencing-based technology, to allow a more in-depth characterisation of the DNA methylation of the murine Fkbp5 locus in three different tissues (blood, frontal cortex and hippocampus). In this study, we not only increased the number of evaluated sites in previously described regulatory regions (in introns 1 and 5), but also extended the evaluation to novel, possibly relevant regulatory regions of the gene (in intron 8, the transcriptional start site, the proximal enhancer and CTCF-binding sites within the 5'UTR). We here describe the assessment of HAM-TBS assays for a panel of 157 CpGs with possible functional relevance in the murine Fkbp5 gene. DNAm profiles were tissue-specific, with lesser differences between the two brain regions than between the brain and blood. Moreover, we identified DNAm changes in the Fkbp5 locus after early life stress exposure in the frontal cortex and blood. Our findings indicate that HAM-TBS is a valuable tool for broader exploration of the DNAm of the murine Fkbp5 locus and its involvement in the stress response.

Transdiagnostic evaluation of epigenetic age acceleration and burden of psychiatric disorders

Different psychiatric disorders as well as exposure to adverse life events have individually been associated with multiple age-related diseases and mortality. Age acceleration in different epigenetic clocks can serve as biomarker for such risk and could help to disentangle the interplay of psychiatric comorbidity and early adversity on age-related diseases and mortality. We evaluated five epigenetic clocks (Horvath, Hannum, PhenoAge, GrimAge and DunedinPoAm) in a transdiagnostic psychiatric sample using epigenome-wide DNA methylation data from peripheral blood of 429 subjects from two studies at the Max Planck Institute of Psychiatry. Burden of psychiatric disease, represented by a weighted score, was significantly associated with biological age acceleration as measured by GrimAge and DunedinPoAm (R2-adj. 0.22 and 0.33 for GrimAge and DunedinPoAm, respectively), but not the other investigated clocks. The relation of burden of psychiatric disease appeared independent of differences in socioeconomic status and medication. Our findings indicate that increased burden of psychiatric disease may associate with accelerated biological aging. This highlights the importance of medical management of patients with multiple psychiatric comorbidities and the potential usefulness of specific epigenetic clocks for early detection of risk and targeted intervention to reduce mortality in psychiatric patients.

The Developing Human Connectome Project Neonatal Data Release

The Developing Human Connectome Project has created a large open science resource which provides researchers with data for investigating typical and atypical brain development across the perinatal period. It has collected 1228 multimodal magnetic resonance imaging (MRI) brain datasets from 1173 fetal and/or neonatal participants, together with collateral demographic, clinical, family, neurocognitive and genomic data from 1173 participants, together with collateral demographic, clinical, family, neurocognitive and genomic data. All subjects were studied in utero and/or soon after birth on a single MRI scanner using specially developed scanning sequences which included novel motion-tolerant imaging methods. Imaging data are complemented by rich demographic, clinical, neurodevelopmental, and genomic information. The project is now releasing a large set of neonatal data; fetal data will be described and released separately. This release includes scans from 783 infants of whom: 583 were healthy infants born at term; as well as preterm infants; and infants at high risk of atypical neurocognitive development. Many infants were imaged more than once to provide longitudinal data, and the total number of datasets being released is 887. We now describe the dHCP image acquisition and processing protocols, summarize the available imaging and collateral data, and provide information on how the data can be accessed.