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.

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.

Condition-specific genetic interaction maps reveal crosstalk between the cAMP/PKA and the HOG MAPK pathways in the activation of the general stress response.

Cells must quickly respond and efficiently adapt to environmental changes. The yeast Saccharomyces cerevisiae has multiple pathways that respond to specific environmental insults, as well as a generic stress response program. The later is regulated by two transcription factors, Msn2 and Msn4, that integrate information from upstream pathways to produce fast, tunable, and robust response to different environmental changes. To understand this integration, we employed a systematic approach to genetically dissect the contribution of various cellular pathways to Msn2/4 regulation under a range of stress and growth conditions. We established a high-throughput liquid handling and automated flow cytometry system and measured GFP levels in 68 single-knockout and 1,566 double-knockout strains that carry an HSP12-GFP allele as a reporter for Msn2/4 activity. Based on the expression of this Msn2/4 reporter in five different conditions, we identified numerous genetic and epistatic interactions between different components in the network upstream to Msn2/4. Our analysis gains new insights into the functional specialization of the RAS paralogs in the repression of stress response and identifies a three-way crosstalk between the Mediator complex, the HOG MAPK pathway, and the cAMP/PKA pathway.

High-resolution chromatin dynamics during a yeast stress response.

Covalent histone modifications are highly conserved and play multiple roles in eukaryotic transcription regulation. Here, we mapped 26 histone modifications genome-wide in exponentially growing yeast and during a dramatic transcriptional reprogramming-the response to diamide stress. We extend prior studies showing that steady-state histone modification patterns reflect genomic processes, especially transcription, and display limited combinatorial complexity. Interestingly, during the stress response we document a modest increase in the combinatorial complexity of histone modification space, resulting from roughly 3% of all nucleosomes transiently populating rare histone modification states. Most of these rare histone states result from differences in the kinetics of histone modification that transiently uncouple highly correlated marks, with slow histone methylation changes often lagging behind the more rapid acetylation changes. Explicit analysis of modification dynamics uncovers ordered sequences of events in gene activation and repression. Together, our results provide a comprehensive view of chromatin dynamics during a massive transcriptional upheaval.

Systematic dissection of roles for chromatin regulators in a yeast stress response.

Packaging of eukaryotic genomes into chromatin has wide-ranging effects on gene transcription. Curiously, it is commonly observed that deletion of a global chromatin regulator affects expression of only a limited subset of genes bound to or modified by the regulator in question. However, in many single-gene studies it has become clear that chromatin regulators often do not affect steady-state transcription, but instead are required for normal transcriptional reprogramming by environmental cues. We therefore have systematically investigated the effects of 83 histone mutants, and 119 gene deletion mutants, on induction/repression dynamics of 170 transcripts in response to diamide stress in yeast. Importantly, we find that chromatin regulators play far more pronounced roles during gene induction/repression than they do in steady-state expression. Furthermore, by jointly analyzing the substrates (histone mutants) and enzymes (chromatin modifier deletions) we identify specific interactions between histone modifications and their regulators. Combining these functional results with genome-wide mapping of several histone marks in the same time course, we systematically investigated the correspondence between histone modification occurrence and function. We followed up on one pathway, finding that Set1-dependent H3K4 methylation primarily acts as a gene repressor during multiple stresses, specifically at genes involved in ribosome biosynthesis. Set1-dependent repression of ribosomal genes occurs via distinct pathways for ribosomal protein genes and ribosomal biogenesis genes, which can be separated based on genetic requirements for repression and based on chromatin changes during gene repression. Together, our dynamic studies provide a rich resource for investigating chromatin regulation, and identify a significant role for the "activating" mark H3K4me3 in gene repression.

Failure of origin activation in response to fork stalling leads to chromosomal instability at fragile sites.

Perturbed DNA replication in early stages of cancer development induces chromosomal instability preferentially at fragile sites. However, the molecular basis for this instability is unknown. Here, we show that even under normal growth conditions, replication fork progression along the fragile site, FRA16C, is slow and forks frequently stall at AT-rich sequences, leading to activation of additional origins to enable replication completion. Under mild replication stress, the frequency of stalling at AT-rich sequences is further increased. Strikingly, unlike in the entire genome, in the FRA16C region additional origins are not activated, suggesting that all potential origins are already activated under normal conditions. Thus, the basis for FRA16C fragility is replication fork stalling at AT-rich sequences and inability to activate additional origins under replication stress. Our results provide a mechanism explaining the replication stress sensitivity of fragile sites and thus, the basis for genomic instability during early stages of cancer development.

Impaired replication stress response in cells from immunodeficiency patients carrying Cernunnos/XLF mutations.

Non-Homologous End Joining (NHEJ) is one of the two major pathways of DNA Double Strand Breaks (DSBs) repair. Mutations in human NHEJ genes can lead to immunodeficiency due to its role in V(D)J recombination in the immune system. In addition, most patients carrying mutations in NHEJ genes display developmental anomalies which are likely the result of a general defect in repair of endogenously induced DSBs such as those arising during normal DNA replication. Cernunnos/XLF is a recently identified NHEJ gene which is mutated in immunodeficiency with microcephaly patients. Here we aimed to investigate whether Cernunnos/XLF mutations disrupt the ability of patient cells to respond to replication stress conditions. Our results demonstrate that Cernunnos/XLF mutated cells and cells downregulated for Cernunnos/XLF have increased sensitivity to conditions which perturb DNA replication. In addition, under replication stress, these cells exhibit impaired DSB repair and increased accumulation of cells in G2/M. Moreover Cernunnos/XLF mutated and down regulated cells display greater chromosomal instability, particularly at fragile sites, under replication stress conditions. These results provide evidence for the role of Cernunnos/XLF in repair of DSBs and maintenance of genomic stability under replication stress conditions. This is the first study of a NHEJ syndrome showing association with impaired cellular response to replication stress conditions. These findings may be related to the clinical features in these patients which are not due to the V(D)J recombination defect. Additionally, in light of the emerging important role of replication stress in the early stages of cancer development, our findings may provide a mechanism for the role of NHEJ in preventing tumorigenesis.

Homologous recombination and nonhomologous end-joining repair pathways regulate fragile site stability.

Common fragile sites are specific loci that form gaps and constrictions on metaphase chromosomes exposed to replication stress, which slows DNA replication. These sites have a role in chromosomal rearrangements in tumors; however, the molecular mechanism of their expression is unclear. Here we show that replication stress leads to focus formation of Rad51 and phosphorylated DNA-PKcs, key components of the homologous recombination (HR) and nonhomologous end-joining (NHEJ), double-strand break (DSB) repair pathways, respectively. Down-regulation of Rad51, DNA-PKcs, or Ligase IV, an additional component of the NHEJ repair pathway, leads to a significant increase in fragile site expression under replication stress. Replication stress also results in focus formation of the DSB markers, MDC1 and gammaH2AX. These foci colocalized with those of Rad51 and phospho-DNA-PKcs. Furthermore, gammaH2AX and phospho-DNA-PKcs foci were localized at expressed fragile sites on metaphase chromosomes. These findings suggest that DSBs are formed at common fragile sites as a result of replication perturbation. The repair of these breaks by both HR and NHEJ pathways is essential for chromosomal stability at these sites.

Molecular basis for expression of common and rare fragile sites.

Fragile sites are specific loci that form gaps, constrictions, and breaks on chromosomes exposed to partial replication stress and are rearranged in tumors. Fragile sites are classified as rare or common, depending on their induction and frequency within the population. The molecular basis of rare fragile sites is associated with expanded repeats capable of adopting unusual non-B DNA structures that can perturb DNA replication. The molecular basis of common fragile sites was unknown. Fragile sites from R-bands are enriched in flexible sequences relative to nonfragile regions from the same chromosomal bands. Here we cloned FRA7E, a common fragile site mapped to a G-band, and revealed a significant difference between its flexibility and that of nonfragile regions mapped to G-bands, similar to the pattern found in R-bands. Thus, in the entire genome, flexible sequences might play a role in the mechanism of fragility. The flexible sequences are composed of interrupted runs of AT-dinucleotides, which have the potential to form secondary structures and hence can affect replication. These sequences show similarity to the AT-rich minisatellite repeats that underlie the fragility of the rare fragile sites FRA16B and FRA10B. We further demonstrate that the normal alleles of FRA16B and FRA10B span the same genomic regions as the common fragile sites FRA16C and FRA10E. Our results suggest that a shared molecular basis, conferred by sequences with a potential to form secondary structures that can perturb replication, may underlie the fragility of rare fragile sites harboring AT-rich minisatellite repeats and aphidicolin-induced common fragile sites.