A DNA methylation atlas of normal human cell types.

DNA methylation is a fundamental epigenetic mark that governs gene expression and chromatin organization, thus providing a window into cellular identity and developmental processes1. Current datasets typically include only a fraction of methylation sites and are often based either on cell lines that underwent massive changes in culture or on tissues containing unspecified mixtures of cells2-5. Here we describe a human methylome atlas, based on deep whole-genome bisulfite sequencing, allowing fragment-level analysis across thousands of unique markers for 39 cell types sorted from 205 healthy tissue samples. Replicates of the same cell type are more than 99.5% identical, demonstrating the robustness of cell identity programmes to environmental perturbation. Unsupervised clustering of the atlas recapitulates key elements of tissue ontogeny and identifies methylation patterns retained since embryonic development. Loci uniquely unmethylated in an individual cell type often reside in transcriptional enhancers and contain DNA binding sites for tissue-specific transcriptional regulators. Uniquely hypermethylated loci are rare and are enriched for CpG islands, Polycomb targets and CTCF binding sites, suggesting a new role in shaping cell-type-specific chromatin looping. The atlas provides an essential resource for study of gene regulation and disease-associated genetic variants, and a wealth of potential tissue-specific biomarkers for use in liquid biopsies.

Universal lung epithelium DNA methylation markers for detection of lung damage in liquid biopsies.

BACKGROUND: Circulating biomarkers for lung damage are lacking. Lung epithelium-specific DNA methylation patterns can potentially report the presence of lung-derived cell-free DNA (cfDNA) in blood, as an indication of lung cell death. METHODS: We sorted human lung alveolar and bronchial epithelial cells from surgical specimens, and obtained their methylomes using whole-genome bisulfite sequencing. We developed a PCR sequencing assay determining the methylation status of 17 loci with lung-specific methylation patterns, and used it to assess lung-derived cfDNA in the plasma of healthy volunteers and patients with lung disease. RESULTS: Loci that are uniquely unmethylated in alveolar or bronchial epithelial cells are enriched for enhancers controlling lung-specific genes. Methylation markers extracted from these methylomes revealed that normal lung cell turnover probably releases cfDNA into the air spaces, rather than to blood. People with advanced lung cancer show a massive elevation of lung cfDNA concentration in blood. Among individuals undergoing bronchoscopy, lung-derived cfDNA is observed in the plasma of those later diagnosed with lung cancer, and to a lesser extent in those diagnosed with other lung diseases. Lung cfDNA is also elevated in patients with acute exacerbation of COPD compared with patients with stable disease, and is associated with future exacerbation and mortality in these patients. CONCLUSIONS: Universal cfDNA methylation markers of normal lung epithelium allow for mutation-independent, sensitive and specific detection of lung-derived cfDNA, reporting on ongoing lung injury. Such markers can find broad utility in the study of normal and pathologic human lung dynamics.

Small RNA sequences derived from pre-microRNAs in the supraspliceosome.

MicroRNAs (miRNAs) are short non-coding RNAs that negatively regulate the expression and translation of genes in healthy and diseased tissues. Herein, we characterize short RNAs from human HeLa cells found in the supraspliceosome, a nuclear dynamic machine in which pre-mRNA processing occurs. We sequenced small RNAs (<200 nt) extracted from the supraspliceosome, and identified sequences that are derived from 200 miRNAs genes. About three quarters of them are mature miRNAs, whereas the rest account for various defined regions of the pre-miRNA, and its hairpin-loop precursor. Out of these aligned sequences, 53 were undetected in cellular extract, and the abundance of additional 48 strongly differed from that in cellular extract. Notably, we describe seven abundant miRNA-derived sequences that overlap non-coding exons of their host gene. The rich collection of sequences identical to pre-miRNAs at the supraspliceosome suggests overlooked nuclear functions. Specifically, the abundant hsa-mir-99b may affect splicing of LINC01129 primary transcript through base-pairing with its exon-intron junction. Using suppression and overexpression experiments, we show that hsa-mir-7704 negatively regulates the level of the lncRNA HAGLR. We claim that in cases of extended base-pairing complementarity, such supraspliceosomal pre-miRNA sequences might have a role in transcription attenuation, maturation and processing.

Spliceosome-Associated microRNAs Signify Breast Cancer Cells and Portray Potential Novel Nuclear Targets.

MicroRNAs (miRNAs) act as negative regulators of gene expression in the cytoplasm. Previous studies have identified the presence of miRNAs in the nucleus. Here we study human breast cancer-derived cell-lines (MCF-7 and MDA-MB-231) and a non-tumorigenic cell-line (MCF-10A) and compare their miRNA sequences at the spliceosome fraction (SF). We report that the levels of miRNAs found in the spliceosome, their identity, and pre-miRNA segmental composition are cell-line specific. One such miRNA is miR-7704 whose genomic position overlaps HAGLR, a cancer-related lncRNA. We detected an inverse expression of miR-7704 and HAGLR in the tested cell lines. Specifically, inhibition of miR-7704 caused an increase in HAGLR expression. Furthermore, elevated levels of miR-7704 slightly altered the cell-cycle in MDA-MB-231. Altogether, we show that SF-miR-7704 acts as a tumor-suppressor gene with HAGLR being its nuclear target. The relative levels of miRNAs found in the spliceosome fractions (e.g., miR-100, miR-30a, and let-7 family) in non-tumorigenic relative to cancer-derived cell-lines was monitored. We found that the expression trend of the abundant miRNAs in SF was different from that reported in the literature and from the observation of large cohorts of breast cancer patients, suggesting that many SF-miRNAs act on targets that are different from the cytoplasmic ones. Altogether, we report on the potential of SF-miRNAs as an unexplored route for cancerous cell state.

DNA Methylation-Based Assessment of Cell Composition in Human Pancreas and Islets.

Assessment of pancreas cell type composition is crucial to the understanding of the genesis of diabetes. Current approaches use immunodetection of protein markers, for example, insulin as a marker of ß-cells. A major limitation of these methods is that protein content varies in physiological and pathological conditions, complicating the extrapolation to actual cell number. Here, we demonstrate the use of cell type-specific DNA methylation markers for determining the fraction of specific cell types in human islet and pancreas specimens. We identified genomic loci that are uniquely demethylated in specific pancreatic cell types and applied targeted PCR to assess the methylation status of these loci in tissue samples, enabling inference of cell type composition. In islet preparations, normalization of insulin secretion to ß-cell DNA revealed similar ß-cell function in pre-type 1 diabetes (T1D), T1D, and type 2 diabetes (T2D), which was significantly lower than in donors without diabetes. In histological pancreas specimens from recent-onset T1D, this assay showed ß-cell fraction within the normal range, suggesting a significant contribution of ß-cell dysfunction. In T2D pancreata, we observed increased α-cell fraction and normal ß-cell fraction. Methylation-based analysis provides an accurate molecular alternative to immune detection of cell types in the human pancreas, with utility in the interpretation of insulin secretion assays and the assessment of pancreas cell composition in health and disease.

Elevated cfDNA after exercise is derived primarily from mature polymorphonuclear neutrophils, with a minor contribution of cardiomyocytes.

Strenuous physical exercise causes a massive elevation in the concentration of circulating cell-free DNA (cfDNA), which correlates with effort intensity and duration. The cellular sources and physiological drivers of this phenomenon are unknown. Using methylation patterns of cfDNA and associated histones, we show that cfDNA in exercise originates mostly in extramedullary polymorphonuclear neutrophils. Strikingly, cardiomyocyte cfDNA concentration increases after a marathon, consistent with elevated troponin levels and indicating low-level, delayed cardiac cell death. Physical impact, low oxygen levels, and elevated core body temperature contribute to neutrophil cfDNA release, while muscle contraction, increased heart rate, ß-adrenergic signaling, or steroid treatment fail to cause elevation of cfDNA. Physical training reduces neutrophil cfDNA release after a standard exercise, revealing an inverse relationship between exercise-induced cfDNA release and training level. We speculate that the release of cfDNA from neutrophils in exercise relates to the activation of neutrophils in the context of exercise-induced muscle damage.

The DNA methylome of human vascular endothelium and its use in liquid biopsies.

BACKGROUND: Vascular endothelial cells (VECs) are an essential component of each tissue, contribute to multiple pathologies, and are targeted by important drugs. Yet, there is a shortage of biomarkers to assess VEC turnover. METHODS: To develop DNA methylation-based liquid biopsies for VECs, we determined the methylome of VECs isolated from freshly dissociated human tissues. FINDINGS: A comparison with a human cell-type methylome atlas yielded thousands of loci that are uniquely unmethylated in VECs. These sites are typically gene enhancers, often residing adjacent to VEC-specific genes. We also identified hundreds of genomic loci that are differentially methylated in organotypic VECs, indicating that VECs feeding specific organs are distinct cell types with a stable epigenetic identity. We established universal and lung-specific VEC markers and evaluated their presence in circulating cell-free DNA (cfDNA). Nearly 2.5% of cfDNA in the plasma of healthy individuals originates from VECs. Sepsis, graft versus host disease, and cardiac catheterization are associated with elevated levels of VEC-derived cfDNA, indicative of vascular damage. Lung-specific VEC cfDNA is selectively elevated in patients with chronic obstructive pulmonary disease (COPD) or lung cancer, revealing tissue-specific vascular turnover. CONCLUSIONS: VEC cfDNA biomarkers inform vascular dynamics in health and disease, potentially contributing to early diagnosis and monitoring of pathologies, and assessment of drug activity. FUNDING: This work was supported by the Beutler Research Program, Helmsley Charitable Trust, JDRF, Grail and the DON Foundation (to Y.D.). Y.D holds the Walter & Greta Stiel Chair in heart studies. B.G., R.S., J.M., D.N., T.K., and Y.D. filed patents on cfDNA analysis.

B cell-derived cfDNA after primary BNT162b2 mRNA vaccination anticipates memory B cells and SARS-CoV-2 neutralizing antibodies.

BACKGROUND: Much remains unknown regarding the response of the immune system to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) vaccination. METHODS: We employed circulating cell-free DNA (cfDNA) to assess the turnover of specific immune cell types following administration of the Pfizer/BioNTech vaccine. FINDINGS: The levels of B cell cfDNA after the primary dose correlated with development of neutralizing antibodies and memory B cells after the booster, revealing a link between early B cell turnover-potentially reflecting affinity maturation-and later development of effective humoral response. We also observed co-elevation of B cell, T cell, and monocyte cfDNA after the booster, underscoring the involvement of innate immune cell turnover in the development of humoral and cellular adaptive immunity. Actual cell counts remained largely stable following vaccination, other than a previously demonstrated temporary reduction in neutrophil and lymphocyte counts. CONCLUSIONS: Immune cfDNA dynamics reveal the crucial role of the primary SARS-CoV-2 vaccine in shaping responses of the immune system following the booster vaccine. FUNDING: This work was supported by a generous gift from Shlomo Kramer. Supported by grants from Human Islet Research Network (HIRN UC4DK116274 and UC4DK104216 to R.S. and Y.D.), Ernest and Bonnie Beutler Research Program of Excellence in Genomic Medicine, The Alex U Soyka Pancreatic Cancer Fund, The Israel Science Foundation, the Waldholtz/Pakula family, the Robert M. and Marilyn Sternberg Family Charitable Foundation, the Helmsley Charitable Trust, Grail, and the DON Foundation (to Y.D.). Y.D. holds the Walter and Greta Stiel Chair and Research Grant in Heart Studies. I.F.-F. received a fellowship from the Glassman Hebrew University Diabetes Center.

Lessons from applied large-scale pooling of 133,816 SARS-CoV-2 RT-PCR tests.

Pooling multiple swab samples before RNA extraction and real-time reverse transcription polymerase chain reaction (RT-PCR) analysis has been proposed as a strategy to reduce costs and increase throughput of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) tests. However, reports on practical large-scale group testing for SARS-CoV-2 have been scant. Key open questions concern reduced sensitivity due to sample dilution, the rate of false positives, the actual efficiency (number of tests saved by pooling), and the impact of infection rate in the population on assay performance. Here, we report an analysis of 133,816 samples collected between April and September 2020 and tested by Dorfman pooling for the presence of SARS-CoV-2. We spared 76% of RNA extraction and RT-PCR tests, despite the frequently changing prevalence (0.5 to 6%). We observed pooling efficiency and sensitivity that exceeded theoretical predictions, which resulted from the nonrandom distribution of positive samples in pools. Overall, our findings support the use of pooling for efficient large-scale SARS-CoV-2 testing.

Remote immune processes revealed by immune-derived circulating cell-free DNA.

Blood cell counts often fail to report on immune processes occurring in remote tissues. Here, we use immune cell type-specific methylation patterns in circulating cell-free DNA (cfDNA) for studying human immune cell dynamics. We characterized cfDNA released from specific immune cell types in healthy individuals (N = 242), cross sectionally and longitudinally. Immune cfDNA levels had no individual steady state as opposed to blood cell counts, suggesting that cfDNA concentration reflects adjustment of cell survival to maintain homeostatic cell numbers. We also observed selective elevation of immune-derived cfDNA upon perturbations of immune homeostasis. Following influenza vaccination (N = 92), B-cell-derived cfDNA levels increased prior to elevated B-cell counts and predicted efficacy of antibody production. Patients with eosinophilic esophagitis (N = 21) and B-cell lymphoma (N = 27) showed selective elevation of eosinophil and B-cell cfDNA, respectively, which were undetectable by cell counts in blood. Immune-derived cfDNA provides a novel biomarker for monitoring immune responses to physiological and pathological processes that are not accessible using conventional methods.