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.

Early sample tagging and pooling enables simultaneous SARS-CoV-2 detection and variant sequencing.

Most severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) diagnostic tests have relied on RNA extraction followed by reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays. Whereas automation improved logistics and different pooling strategies increased testing capacity, highly multiplexed next-generation sequencing (NGS) diagnostics remain a largely untapped resource. NGS tests have the potential to markedly increase throughput while providing crucial SARS-CoV-2 variant information. Current NGS-based detection and genotyping assays for SARS-CoV-2 are costly, mostly due to parallel sample processing through multiple steps. Here, we have established ApharSeq, in which samples are barcoded in the lysis buffer and pooled before reverse transcription. We validated this assay by applying ApharSeq to more than 500 clinical samples from the Clinical Virology Laboratory at Hadassah hospital in a robotic workflow. The assay was linear across five orders of magnitude, and the limit of detection was Ct 33 (~1000 copies/ml, 95% sensitivity) with >99.5% specificity. ApharSeq provided targeted high-confidence genotype information due to unique molecular identifiers incorporated into this method. Because of early pooling, we were able to estimate a 10- to 100-fold reduction in labor, automated liquid handling, and reagent requirements in high-throughput settings compared to current testing methods. The protocol can be tailored to assay other host or pathogen RNA targets simultaneously. These results suggest that ApharSeq can be a promising tool for current and future mass diagnostic challenges.

ChIP-seq of plasma cell-free nucleosomes identifies gene expression programs of the cells of origin.

Cell-free DNA (cfDNA) in human plasma provides access to molecular information about the pathological processes in the organs or tumors from which it originates. These DNA fragments are derived from fragmented chromatin in dying cells and retain some of the cell-of-origin histone modifications. In this study, we applied chromatin immunoprecipitation of cell-free nucleosomes carrying active chromatin modifications followed by sequencing (cfChIP-seq) to 268 human samples. In healthy donors, we identified bone marrow megakaryocytes, but not erythroblasts, as major contributors to the cfDNA pool. In patients with a range of liver diseases, we showed that we can identify pathology-related changes in hepatocyte transcriptional programs. In patients with metastatic colorectal carcinoma, we detected clinically relevant and patient-specific information, including transcriptionally active human epidermal growth factor receptor 2 (HER2) amplifications. Altogether, cfChIP-seq, using low sequencing depth, provides systemic and genome-wide information and can inform diagnosis and facilitate interrogation of physiological and pathological processes using blood samples.

Author Correction: Deciphering eukaryotic gene-regulatory logic with 100 million random promoters.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Deciphering eukaryotic gene-regulatory logic with 100 million random promoters.

How transcription factors (TFs) interpret cis-regulatory DNA sequence to control gene expression remains unclear, largely because past studies using native and engineered sequences had insufficient scale. Here, we measure the expression output of >100 million synthetic yeast promoter sequences that are fully random. These sequences yield diverse, reproducible expression levels that can be explained by their chance inclusion of functional TF binding sites. We use machine learning to build interpretable models of transcriptional regulation that predict ~94% of the expression driven from independent test promoters and ~89% of the expression driven from native yeast promoter fragments. These models allow us to characterize each TF's specificity, activity and interactions with chromatin. TF activity depends on binding-site strand, position, DNA helical face and chromatin context. Notably, expression level is influenced by weak regulatory interactions, which confound designed-sequence studies. Our analyses show that massive-throughput assays of fully random DNA can provide the big data necessary to develop complex, predictive models of gene regulation.

Fine-Resolution Mapping of TF Binding and Chromatin Interactions.

Transcription factor (TF) binding to DNA is crucial for transcriptional regulation. There are multiple methods for mapping such binding. These methods balance between input requirements, spatial resolution, and compatibility with high-throughput automation. Here, we describe SLIM-ChIP (short-fragment-enriched, low-input, indexed MNase ChIP), which combines enzymatic fragmentation of chromatin and on-bead indexing to address these desiderata. SLIM-ChIP reproduces a high-resolution binding map of yeast Reb1 comparable with existing methods, yet with less input material and full compatibility with high-throughput procedures. We demonstrate the robustness and flexibility of SLIM-ChIP by probing additional factors in yeast and mouse. Finally, we show that SLIM-ChIP provides information on the chromatin landscape surrounding the bound transcription factor. We identify a class of Reb1 sites where the proximal -1 nucleosome tightly interacts with Reb1 and maintains unidirectional transcription. SLIM-ChIP is an attractive solution for mapping DNA binding proteins and charting the surrounding chromatin occupancy landscape at a single-cell level.

Mapping the Landscape of a Eukaryotic Degronome.

The ubiquitin-proteasome system (UPS) for protein degradation has been under intensive study, and yet, we have only partial understanding of mechanisms by which proteins are selected to be targeted for proteolysis. One of the obstacles in studying these recognition pathways is the limited repertoire of known degradation signals (degrons). To better understand what determines the susceptibility of intracellular proteins to degradation by the UPS, we developed an unbiased method for large-scale identification of eukaryotic degrons. Using a reporter-based high-throughput competition assay, followed by deep sequencing, we measured a degradation potency index for thousands of native polypeptides in a single experiment. We further used this method to identify protein quality control (PQC)-specific and compartment-specific degrons. Our method provides an unprecedented insight into the yeast degronome, and it can readily be modified to study protein degradation signals and pathways in other organisms and in various settings.

Elucidating Combinatorial Chromatin States at Single-Nucleosome Resolution.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental to our current view of chromatin structure and function. It allows genome-wide mapping of histone marks, which demarcate biologically relevant domains. However, ChIP-seq is an ensemble measurement reporting the average occupancy of individual marks in a cell population. Consequently, our understanding of the combinatorial nature of chromatin states relies almost exclusively on correlation between the genomic distributions of individual marks. Here, we report the development of combinatorial-iChIP to determine the genome-wide co-occurrence of histone marks at single-nucleosome resolution. By comparing to a null model, we show that certain combinations of overlapping marks (H3K36me3 and H3K79me3) co-occur more frequently than would be expected by chance, while others (H3K4me3 and H3K36me3) do not, reflecting differences in the underlying chromatin pathways. We further use combinatorial-iChIP to illuminate aspects of the Set2-RPD3S pathway. This approach promises to improve our understanding of the combinatorial complexity of chromatin.

Microbial-derived lithocholic acid and vitamin K2 drive the metabolic maturation of pluripotent stem cells-derived and fetal hepatocytes.

UNLABELLED: The liver is the main organ responsible for the modification, clearance, and transformational toxicity of most xenobiotics owing to its abundance in cytochrome P450 (CYP450) enzymes. However, the scarcity and variability of primary hepatocytes currently limits their utility. Human pluripotent stem cells (hPSCs) represent an excellent source of differentiated hepatocytes; however, current protocols still produce fetal-like hepatocytes with limited mature function. Interestingly, fetal hepatocytes acquire mature CYP450 expression only postpartum, suggesting that nutritional cues may drive hepatic maturation. We show that vitamin K2 and lithocholic acid, a by-product of intestinal flora, activate pregnane X receptor (PXR) and subsequent CYP3A4 and CYP2C9 expression in hPSC-derived and isolated fetal hepatocytes. Differentiated cells produce albumin and apolipoprotein B100 at levels equivalent to primary human hepatocytes, while demonstrating an 8-fold induction of CYP450 activity in response to aryl hydrocarbon receptor (AhR) agonist omeprazole and a 10-fold induction in response to PXR agonist rifampicin. Flow cytometry showed that over 83% of cells were albumin and hepatocyte nuclear factor 4 alpha (HNF4α) positive, permitting high-content screening in a 96-well plate format. Analysis of 12 compounds showed an R(2) correlation of 0.94 between TC50 values obtained in stem cell-derived hepatocytes and primary cells, compared to 0.62 for HepG2 cells. Finally, stem cell-derived hepatocytes demonstrate all toxicological endpoints examined, including steatosis, apoptosis, and cholestasis, when exposed to nine known hepatotoxins. CONCLUSION: Our work provides fresh insights into liver development, suggesting that microbial-derived cues may drive the maturation of CYP450 enzymes postpartum. Addition of these cues results in the first functional, inducible, hPSC-derived hepatocyte for predictive toxicology.

Genome-wide "re"-modeling of nucleosome positions.

Two recent studies mapped nucleosomes across the yeast and human genomes, teasing apart the relative contributions of DNA sequence and chromatin remodelers to nucleosome organization. These data suggest two emerging models: chromatin remodelers position nucleosomes around transcriptional start sites in yeast, and a few "locked" nucleosomes may serve as barriers from which nucleosome arrays emanate in human genomes.

PRAJA1 is a ubiquitin ligase for the polycomb repressive complex 2 proteins.

Methylation of lysine 27 on histone H3 by the polycomb repressive complex 2 (PRC2) leads to transcriptional repression of genes which are critical to development. PRC2 core complex is composed of the histone methyltransferase EZH2, EED, and SUZ12. Knockdown of any of the PRC2 core subunits results in a concomitant loss of the other subunits which is mediated by the ubiquitin (Ub)-proteasome system (UPS). Inhibition of cellular methyltransferases by 3-deazaneplanocin A (DZNep) also leads to dissociation of the PRC2 complex and rapid degradation of its subunits. Interestingly, the expression of several Ub ligases was induced following DZNep treatment, suggesting that PRC2 might repress the Ub ligase(s) that target its subunits for degradation. Here we confirm that individual PRC2 subunits are ubiquitinated and rapidly degraded by the proteasome. One of the DZNep-induced Ub ligases, PRAJA1, can target PRC2 subunits for proteasomal degradation. PRAJA1 directly ubiquitinates individual PRC2 subunits in a cell free system, which leads to their proteasomal degradation. Expression of PRAJA1 but not of an inactive RING finger mutant of the protein, enhanced the degradation of individual PRC2 subunits in cells. Taken together, our results suggest a role for PRAJA1 in regulating the level of PRC2 by targeting its free subunits for Ub-mediated proteasomal degradation.

The N-terminal domain of MyoD is necessary and sufficient for its nuclear localization-dependent degradation by the ubiquitin system.

A growing number of proteins, including the myogenic transcription factor MyoD, are targeted for proteasomal degradation after N-terminal ubiquitination (NTU) where the first ubiquitin moiety is conjugated to the N-terminal residue rather than to an internal lysine. NTU might be essential in targeting both lysine-containing and naturally occurring lysine-less proteins such as p16(INK4a) and p14(ARF); however, the mechanisms that underlie this process are largely unknown. Specifically, the recognition motif(s) in the target substrates and the ubiquitin ligase(s) that catalyze NTU are still obscure. Here we show that the N-terminal domain of MyoD is critical for its degradation and that its destabilizing effect depends on nuclear localization of the protein. Deletion of the first 15 aa of MyoD blocked completely its lysine-independent degradation. Importantly, transfer of the first 30 N-terminal residues of MyoD to GFP destabilized this otherwise stable protein, and, here too, targeting for degradation depended on localization of the protein to the nucleus. Deletion of the N-terminal domain of lysine-less MyoD did not abolish completely ubiquitination of the protein, suggesting that this domain may be required for targeting the protein also in a postubiquitination step. Interestingly, NTU is evolutionarily conserved: in the yeast Saccharomyces cerevisiae lysine-less (LL) MyoD is degraded in a ubiquitin-, N-terminal domain-, and nuclear localization-dependent manner. Taken together, our data suggest that a short N-terminal segment of MyoD is necessary and sufficient to render MyoD susceptible for ubiquitin- and nuclear-dependent degradation.