Ribosomal DNA copy number is associated with body mass in humans and other mammals.

Body mass results from a complex interplay between genetics and environment. Previous studies of the genetic contribution to body mass have excluded repetitive regions due to the technical limitations of platforms used for population scale studies. Here we apply genome-wide approaches, identifying an association between adult body mass and the copy number (CN) of 47S-ribosomal DNA (rDNA). rDNA codes for the 18 S, 5.8 S and 28 S ribosomal RNA (rRNA) components of the ribosome. In mammals, there are hundreds of copies of these genes. Inter-individual variation in the rDNA CN has not previously been associated with a mammalian phenotype. Here, we show that rDNA CN variation associates with post-pubertal growth rate in rats and body mass index in adult humans. rDNA CN is not associated with rRNA transcription rates in adult tissues, suggesting the mechanistic link occurs earlier in development. This aligns with the observation that the association emerges by early adulthood.

Ribosomal DNA copy number variation associates with hematological profiles and renal function in the UK Biobank.

The phenotypic impact of genetic variation of repetitive features in the human genome is currently understudied. One such feature is the multi-copy 47S ribosomal DNA (rDNA) that codes for rRNA components of the ribosome. Here, we present an analysis of rDNA copy number (CN) variation in the UK Biobank (UKB). From the first release of UKB whole-genome sequencing (WGS) data, a discovery analysis in White British individuals reveals that rDNA CN associates with altered counts of specific blood cell subtypes, such as neutrophils, and with the estimated glomerular filtration rate, a marker of kidney function. Similar trends are observed in other ancestries. A range of analyses argue against reverse causality or common confounder effects, and all core results replicate in the second UKB WGS release. Our work demonstrates that rDNA CN is a genetic influence on trait variance in humans.

Donor whole blood DNA methylation is not a strong predictor of acute graft versus host disease in unrelated donor allogeneic haematopoietic cell transplantation.

Allogeneic hematopoietic cell transplantation (HCT) is used to treat many blood-based disorders and malignancies, however it can also result in serious adverse events, such as the development of acute graft-versus-host disease (aGVHD). This study aimed to develop a donor-specific epigenetic classifier to reduce incidence of aGVHD by improving donor selection. Genome-wide DNA methylation was assessed in a discovery cohort of 288 HCT donors selected based on recipient aGVHD outcome; this cohort consisted of 144 cases with aGVHD grades III-IV and 144 controls with no aGVHD. We applied a machine learning algorithm to identify CpG sites predictive of aGVHD. Receiver operating characteristic (ROC) curve analysis of these sites resulted in a classifier with an encouraging area under the ROC curve (AUC) of 0.91. To test this classifier, we used an independent validation cohort (n = 288) selected using the same criteria as the discovery cohort. Attempts to validate the classifier failed with the AUC falling to 0.51. These results indicate that donor DNA methylation may not be a suitable predictor of aGVHD in an HCT setting involving unrelated donors, despite the initial promising results in the discovery cohort. Our work highlights the importance of independent validation of machine learning classifiers, particularly when developing classifiers intended for clinical use.

Hemin availability induces coordinated DNA methylation and gene expression changes in Porphyromonas gingivalis.

Periodontal disease is a chronic inflammatory disease in which the oral pathogen Porphyromonas gingivalis plays an important role. Porphyromonas gingivalis expresses virulence determinants in response to higher hemin concentrations, but the underlying regulatory processes remain unclear. Bacterial DNA methylation has the potential to fulfil this mechanistic role. We characterized the methylome of P. gingivalis, and compared its variation to transcriptome changes in response to hemin availability. Porphyromonas gingivalis W50 was grown in chemostat continuous culture with excess or limited hemin, prior to whole-methylome and transcriptome profiling using Nanopore and Illumina RNA-Seq. DNA methylation was quantified for Dam/Dcm motifs and all-context N6-methyladenine (6mA) and 5-methylcytosine (5mC). Of all 1,992 genes analyzed, 161 and 268 were respectively over- and under-expressed with excess hemin. Notably, we detected differential DNA methylation signatures for the Dam "GATC" motif and both all-context 6mA and 5mC in response to hemin availability. Joint analyses identified a subset of coordinated changes in gene expression, 6mA, and 5mC methylation that target genes involved in lactate utilization and ABC transporters. The results identify altered methylation and expression responses to hemin availability in P. gingivalis, with insights into mechanisms regulating its virulence in periodontal disease. IMPORTANCE DNA methylation has important roles in bacteria, including in the regulation of transcription. Porphyromonas gingivalis, an oral pathogen in periodontitis, exhibits well-established gene expression changes in response to hemin availability. However, the regulatory processes underlying these effects remain unknown. We profiled the novel P. gingivalis epigenome, and assessed epigenetic and transcriptome variation under limited and excess hemin conditions. As expected, multiple gene expression changes were detected in response to limited and excess hemin that reflect health and disease, respectively. Notably, we also detected differential DNA methylation signatures for the Dam "GATC" motif and both all-context 6mA and 5mC in response to hemin. Joint analyses identified coordinated changes in gene expression, 6mA, and 5mC methylation that target genes involved in lactate utilization and ABC transporters. The results identify novel regulatory processes underlying the mechanism of hemin regulated gene expression in P. gingivalis, with phenotypic impacts on its virulence in periodontal disease.

Genetic variation at mouse and human ribosomal DNA influences associated epigenetic states.

BACKGROUND: Ribosomal DNA (rDNA) displays substantial inter-individual genetic variation in human and mouse. A systematic analysis of how this variation impacts epigenetic states and expression of the rDNA has thus far not been performed. RESULTS: Using a combination of long- and short-read sequencing, we establish that 45S rDNA units in the C57BL/6J mouse strain exist as distinct genetic haplotypes that influence the epigenetic state and transcriptional output of any given unit. DNA methylation dynamics at these haplotypes are dichotomous and life-stage specific: at one haplotype, the DNA methylation state is sensitive to the in utero environment, but refractory to post-weaning influences, whereas other haplotypes entropically gain DNA methylation during aging only. On the other hand, individual rDNA units in human show limited evidence of genetic haplotypes, and hence little discernible correlation between genetic and epigenetic states. However, in both species, adjacent units show similar epigenetic profiles, and the overall epigenetic state at rDNA is strongly positively correlated with the total rDNA copy number. Analysis of different mouse inbred strains reveals that in some strains, such as 129S1/SvImJ, the rDNA copy number is only approximately 150 copies per diploid genome and DNA methylation levels are < 5%. CONCLUSIONS: Our work demonstrates that rDNA-associated genetic variation has a considerable influence on rDNA epigenetic state and consequently rRNA expression outcomes. In the future, it will be important to consider the impact of inter-individual rDNA (epi)genetic variation on mammalian phenotypes and diseases.

An integrative epi-transcriptomic approach identifies the human cartilage chitinase 3-like protein 2 (CHI3L2) as a potential mediator of B12 deficiency in adipocytes.

Vitamin B12 has multiple biochemical functions including in the one-carbon cycle generating a methyl group for DNA methylation, and metabolism of fatty acids and amino acids to generate energy via the citric acid cycle. The aim of our study was to use a combined epigenomic and transcriptomic approach to identify novel genes mediating the effect of B12 on adipogenesis.Human pre-adipocytes (CHUB-S7) were treated with a range of B12 (0-500 nM) concentrations from the day of cell seeding until harvesting in discovery and validation experiments prior to genome-wide methylation analysis using the Illumina HumanMethylation 450Beadchip. For transcriptomic analysis, RNA-seq libraries were run on the Illumina HiSeq 2500. To further investigate the expression of any genes on human adipogenesis, a second human preadipocyte strain was studied (SGBS) by real-time quantitative PCR (qRT-PCR).A combined epigenetic and transcriptomic approach in differentiated human pre-adipocyte cell line, CHUB-S7, identified that the Human cartilage chitinase 3-like protein 2 (CHI3L2) gene was hypo-methylated and had increased expression in low B12 conditions. Furthermore, there was an approximately 1000-fold increase in CHI3L2 expression in the early days of adipocyte differentiation, which paralleled an increase of lipid droplets in differentiated SGBS cells and an increased expression level of markers of mature adipocytes.In summary, we have identified a potential role of the human cartilage chitinase 3-like protein 2 (CHI3L2) in adipocyte function in the presence of low B12 levels.

Comparative epigenetic analysis of tumour initiating cells and syngeneic EPSC-derived neural stem cells in glioblastoma.

Epigenetic mechanisms which play an essential role in normal developmental processes, such as self-renewal and fate specification of neural stem cells (NSC) are also responsible for some of the changes in the glioblastoma (GBM) genome. Here we develop a strategy to compare the epigenetic and transcriptional make-up of primary GBM cells (GIC) with patient-matched expanded potential stem cell (EPSC)-derived NSC (iNSC). Using a comparative analysis of the transcriptome of syngeneic GIC/iNSC pairs, we identify a glycosaminoglycan (GAG)-mediated mechanism of recruitment of regulatory T cells (Tregs) in GBM. Integrated analysis of the transcriptome and DNA methylome of GBM cells identifies druggable target genes and patient-specific prediction of drug response in primary GIC cultures, which is validated in 3D and in vivo models. Taken together, we provide a proof of principle that this experimental pipeline has the potential to identify patient-specific disease mechanisms and druggable targets in GBM.

The DNA methylome of human sperm is distinct from blood with little evidence for tissue-consistent obesity associations.

Epidemiological research suggests that paternal obesity may increase the risk of fathering small for gestational age offspring. Studies in non-human mammals indicate that such associations could be mediated by DNA methylation changes in spermatozoa that influence offspring development in utero. Human obesity is associated with differential DNA methylation in peripheral blood. It is unclear, however, whether this differential DNA methylation is reflected in spermatozoa. We profiled genome-wide DNA methylation using the Illumina MethylationEPIC array in a cross-sectional study of matched human blood and sperm from lean (discovery n = 47; replication n = 21) and obese (n = 22) males to analyse tissue covariation of DNA methylation, and identify obesity-associated methylomic signatures. We found that DNA methylation signatures of human blood and spermatozoa are highly discordant, and methylation levels are correlated at only a minority of CpG sites (~1%). At the majority of these sites, DNA methylation appears to be influenced by genetic variation. Obesity-associated DNA methylation in blood was not generally reflected in spermatozoa, and obesity was not associated with altered covariation patterns or accelerated epigenetic ageing in the two tissues. However, one cross-tissue obesity-specific hypermethylated site (cg19357369; chr4:2429884; P = 8.95 × 10-8; 2% DNA methylation difference) was identified, warranting replication and further investigation. When compared to a wide range of human somatic tissue samples (n = 5,917), spermatozoa displayed differential DNA methylation across pathways enriched in transcriptional regulation. Overall, human sperm displays a unique DNA methylation profile that is highly discordant to, and practically uncorrelated with, that of matched peripheral blood. We observed that obesity was only nominally associated with differential DNA methylation in sperm, and therefore suggest that spermatozoal DNA methylation is an unlikely mediator of intergenerational effects of metabolic traits.

DNA methylation clocks as a predictor for ageing and age estimation in naked mole-rats, Heterocephalus glaber.

The naked mole-rat, Heterocephalus glaber (NMR), the longest-lived rodent, is of significance and interest in the study of biomarkers for ageing. Recent breakthroughs in this field have revealed 'epigenetic clocks' that are based on the temporal accumulation of DNA methylation at specific genomic sites. Here, we validate the hypothesis of an epigenetic clock in NMRs based on changes in methylation of targeted CpG sites. We initially analysed 51 CpGs in NMR livers spanning an age range of 39-1,144 weeks and found 23 to be significantly associated with age (p<0.05). We then built a predictor of age using these sites. To test the accuracy of this model, we analysed an additional set of liver samples, and were successfully able to predict their age with a root mean squared error of 166 weeks. We also profiled skin samples with the same age range, finding a striking correlation between their predicted age versus their actual age (R=0.93), but which was lower when compared to the liver, suggesting that skin ages slower than the liver in NMRs. Our model will enable the prediction of age in wild-caught and captive NMRs of unknown age, and will be invaluable for further mechanistic studies of mammalian ageing.

DNA methylation aging clocks: challenges and recommendations.

Epigenetic clocks comprise a set of CpG sites whose DNA methylation levels measure subject age. These clocks are acknowledged as a highly accurate molecular correlate of chronological age in humans and other vertebrates. Also, extensive research is aimed at their potential to quantify biological aging rates and test longevity or rejuvenating interventions. Here, we discuss key challenges to understand clock mechanisms and biomarker utility. This requires dissecting the drivers and regulators of age-related changes in single-cell, tissue- and disease-specific models, as well as exploring other epigenomic marks, longitudinal and diverse population studies, and non-human models. We also highlight important ethical issues in forensic age determination and predicting the trajectory of biological aging in an individual.

BulkVis: a graphical viewer for Oxford nanopore bulk FAST5 files.

MOTIVATION: The Oxford Nanopore Technologies (ONT) MinION is used for sequencing a wide variety of sample types with diverse methods of sample extraction. Nanopore sequencers output FAST5 files containing signal data subsequently base called to FASTQ format. Optionally, ONT devices can collect data from all sequencing channels simultaneously in a bulk FAST5 file enabling inspection of signal in any channel at any point. We sought to visualize this signal to inspect challenging or difficult to sequence samples. RESULTS: The BulkVis tool can load a bulk FAST5 file and overlays MinKNOW (the software that controls ONT sequencers) classifications on the signal trace and can show mappings to a reference. Users can navigate to a channel and time or, given a FASTQ header from a read, jump to its specific position. BulkVis can export regions as Nanopore base caller compatible reads. Using BulkVis, we find long reads can be incorrectly divided by MinKNOW resulting in single DNA molecules being split into two or more reads. The longest seen to date is 2 272 580 bases in length and reported in eleven consecutive reads. We provide helper scripts that identify and reconstruct split reads given a sequencing summary file and alignment to a reference. We note that incorrect read splitting appears to vary according to input sample type and is more common in 'ultra-long' read preparations. AVAILABILITY AND IMPLEMENTATION: The software is available freely under an MIT license at https://github.com/LooseLab/bulkvis. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Early life diet conditions the molecular response to post-weaning protein restriction in the mouse.

BACKGROUND: Environmental influences fluctuate throughout the life course of an organism. It is therefore important to understand how the timing of exposure impacts molecular responses. Herein, we examine the responses of two key molecular markers of dietary stress, namely variant-specific methylation at ribosomal DNA (rDNA) and small RNA distribution, including tRNA fragments, in a mouse model of protein restriction (PR) with exposure at pre- and/or post-weaning. RESULTS: We first confirm that pre-weaning PR exposure modulates the methylation state of rDNA in a genotype-dependent manner, whereas post-weaning PR exposure has no such effect. Conversely, post-weaning PR induces a shift in small RNA distribution, but there is no effect in the pre-weaning PR model. Intriguingly, mice exposed to PR throughout their lives show neither of these two dietary stress markers, similar to controls. CONCLUSIONS: The results show that the timing of the insult affects the nature of the molecular response but also, critically, that 'matching' diet exposure either side of weaning eliminates the stress response at the level of rDNA methylation and small RNA in sperm.

Integration of human pancreatic islet genomic data refines regulatory mechanisms at Type 2 Diabetes susceptibility loci.

Human genetic studies have emphasised the dominant contribution of pancreatic islet dysfunction to development of Type 2 Diabetes (T2D). However, limited annotation of the islet epigenome has constrained efforts to define the molecular mechanisms mediating the, largely regulatory, signals revealed by Genome-Wide Association Studies (GWAS). We characterised patterns of chromatin accessibility (ATAC-seq, n = 17) and DNA methylation (whole-genome bisulphite sequencing, n = 10) in human islets, generating high-resolution chromatin state maps through integration with established ChIP-seq marks. We found enrichment of GWAS signals for T2D and fasting glucose was concentrated in subsets of islet enhancers characterised by open chromatin and hypomethylation, with the former annotation predominant. At several loci (including CDC123, ADCY5, KLHDC5) the combination of fine-mapping genetic data and chromatin state enrichment maps, supplemented by allelic imbalance in chromatin accessibility pinpointed likely causal variants. The combination of increasingly-precise genetic and islet epigenomic information accelerates definition of causal mechanisms implicated in T2D pathogenesis.

Ageing-associated DNA methylation dynamics are a molecular readout of lifespan variation among mammalian species.

BACKGROUND: Mammalian species exhibit a wide range of lifespans. To date, a robust and dynamic molecular readout of these lifespan differences has not yet been identified. Recent studies have established the existence of ageing-associated differentially methylated positions (aDMPs) in human and mouse. These are CpG sites at which DNA methylation dynamics show significant correlations with age. We hypothesise that aDMPs are pan-mammalian and are a dynamic molecular readout of lifespan variation among different mammalian species. RESULTS: A large-scale integrated analysis of aDMPs in six different mammals reveals a strong negative relationship between rate of change of methylation levels at aDMPs and lifespan. This relationship also holds when comparing two different dog breeds with known differences in lifespans. In an ageing cohort of aneuploid mice carrying a complete copy of human chromosome 21, aDMPs accumulate far more rapidly than is seen in human tissues, revealing that DNA methylation at aDMP sites is largely shaped by the nuclear trans-environment and represents a robust molecular readout of the ageing cellular milieu. CONCLUSIONS: Overall, we define the first dynamic molecular readout of lifespan differences among mammalian species and propose that aDMPs will be an invaluable molecular tool for future evolutionary and mechanistic studies aimed at understanding the biological factors that determine lifespan in mammals.

Roadmap for investigating epigenome deregulation and environmental origins of cancer.

The interaction between the (epi)genetic makeup of an individual and his/her environmental exposure record (exposome) is accepted as a determinant factor for a significant proportion of human malignancies. Recent evidence has highlighted the key role of epigenetic mechanisms in mediating gene-environment interactions and translating exposures into tumorigenesis. There is also growing evidence that epigenetic changes may be risk factor-specific ("fingerprints") that should prove instrumental in the discovery of new biomarkers in cancer. Here, we review the state of the science of epigenetics associated with environmental stimuli and cancer risk, highlighting key developments in the field. Critical knowledge gaps and research needs are discussed and advances in epigenomics that may help in understanding the functional relevance of epigenetic alterations. Key elements required for causality inferences linking epigenetic changes to exposure and cancer are discussed and how these alterations can be incorporated in carcinogen evaluation and in understanding mechanisms underlying epigenome deregulation by the environment.

DNA methylation profiles between airway epithelium and proxy tissues in children.

AIM: Epidemiological studies of deoxyribonucleic acid (DNA) methylation in airway disease have largely been conducted using blood or buccal samples. However, given tissue specificity of DNA methylation, these surrogate tissues may not allow reliable inferences about methylation in the lung. We sought to compare the pattern of DNA methylation in blood, buccal and nasal epithelial cells to that in airway epithelial cells from children. METHODS: Samples of blood, and buccal, nasal and airway epithelium were obtained from six children undergoing elective anaesthesia for adenotonsillectomy. DNA methylation was assessed at 450 000 5'-C-phosphate-G-3' (CpG) sites using the Illumina HumanMethylation450 array. RESULTS: Eighteen samples from all sites were suitable for analysis. Hierarchical clustering demonstrated that the methylation profile in nasal epithelium was most representative of that in airway epithelium; the profile in buccal cells was moderately similar and that in blood was least similar. CONCLUSION: DNA methylation in blood poorly reflects methylation in airway epithelium. Future epidemiological studies of DNA methylation and airway diseases should consider measurement of methylation either in buccal cells or, preferably, in nasal epithelial cells.

eFORGE: A Tool for Identifying Cell Type-Specific Signal in Epigenomic Data.

Epigenome-wide association studies (EWAS) provide an alternative approach for studying human disease through consideration of non-genetic variants such as altered DNA methylation. To advance the complex interpretation of EWAS, we developed eFORGE (http://eforge.cs.ucl.ac.uk/), a new standalone and web-based tool for the analysis and interpretation of EWAS data. eFORGE determines the cell type-specific regulatory component of a set of EWAS-identified differentially methylated positions. This is achieved by detecting enrichment of overlap with DNase I hypersensitive sites across 454 samples (tissues, primary cell types, and cell lines) from the ENCODE, Roadmap Epigenomics, and BLUEPRINT projects. Application of eFORGE to 20 publicly available EWAS datasets identified disease-relevant cell types for several common diseases, a stem cell-like signature in cancer, and demonstrated the ability to detect cell-composition effects for EWAS performed on heterogeneous tissues. Our approach bridges the gap between large-scale epigenomics data and EWAS-derived target selection to yield insight into disease etiology.

Increased DNA methylation variability in type 1 diabetes across three immune effector cell types.

The incidence of type 1 diabetes (T1D) has substantially increased over the past decade, suggesting a role for non-genetic factors such as epigenetic mechanisms in disease development. Here we present an epigenome-wide association study across 406,365 CpGs in 52 monozygotic twin pairs discordant for T1D in three immune effector cell types. We observe a substantial enrichment of differentially variable CpG positions (DVPs) in T1D twins when compared with their healthy co-twins and when compared with healthy, unrelated individuals. These T1D-associated DVPs are found to be temporally stable and enriched at gene regulatory elements. Integration with cell type-specific gene regulatory circuits highlight pathways involved in immune cell metabolism and the cell cycle, including mTOR signalling. Evidence from cord blood of newborns who progress to overt T1D suggests that the DVPs likely emerge after birth. Our findings, based on 772 methylomes, implicate epigenetic changes that could contribute to disease pathogenesis in T1D.