Search

IGHMBP2 is a nonessential, superfamily 1 DNA/RNA helicase that is mutated in patients with rare neuromuscular diseases SMARD1 and CMT2S. IGHMBP2 is implicated in translational and transcriptional regulation via biochemical association with ribosomal proteins, pre-rRNA processing factors, and tRNA-related species. To uncover the cellular consequences of perturbing IGHMBP2, we generated full and partial IGHMBP2 deletion K562 cell lines. Using polysome profiling and a nascent protein synthesis assay, we found that IGHMBP2 deletion modestly reduces global translation. We performed Ribo-seq and RNA-seq and identified diverse gene expression changes due to IGHMBP2 deletion, including ATF4 up-regulation. With recent studies showing the integrated stress response (ISR) can contribute to tRNA metabolism-linked neuropathies, we asked whether perturbing IGHMBP2 promotes ISR activation. We generated ATF4 reporter cell lines and found IGHMBP2 knockout cells demonstrate basal, chronic ISR activation. Our work expands upon the impact of IGHMBP2 in translation and elucidates molecular mechanisms that may link mutant IGHMBP2 to severe clinical phenotypes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Prime editing enables the precise modification of genomes through reverse transcription of template sequences appended to the 3' ends of CRISPR-Cas guide RNAs1. To identify cellular determinants of prime editing, we developed scalable prime editing reporters and performed genome-scale CRISPR-interference screens. From these screens, a single factor emerged as the strongest mediator of prime editing: the small RNA-binding exonuclease protection factor La. Further investigation revealed that La promotes prime editing across approaches (PE2, PE3, PE4 and PE5), edit types (substitutions, insertions and deletions), endogenous loci and cell types but has no consistent effect on genome-editing approaches that rely on standard, unextended guide RNAs. Previous work has shown that La binds polyuridine tracts at the 3' ends of RNA polymerase III transcripts2. We found that La functionally interacts with the 3' ends of polyuridylated prime editing guide RNAs (pegRNAs). Guided by these results, we developed a prime editor protein (PE7) fused to the RNA-binding, N-terminal domain of La. This editor improved prime editing with expressed pegRNAs and engineered pegRNAs (epegRNAs), as well as with synthetic pegRNAs optimized for La binding. Together, our results provide key insights into how prime editing components interact with the cellular environment and suggest general strategies for stabilizing exogenous small RNAs therein.

Introduction: The Vel– phenotype is a rare blood group, and it is challenging for identifying this phenotype due to limited available reagents. Moreover, there are relatively few studies on genomic editing of erythroid antigens and generation of knockout (KO) cell lines at present. Methods: To identify the high-efficiency small-guiding RNA (sgRNA) sequence, candidate sgRNAs were transfected into HEK 293T cells and analyzed using Sanger sequencing. Following this, the high-efficiency sgRNA was transfected into K562 cells using lentivirus transduction to generate KO Vel blood group gene cells. The expression of the Vel protein was detected using Western blot on single-cell clones. Additionally, flow cytometry was used to detect the erythroid markers CD235a and CD71. Hemoglobin quantification and Giemsa staining were also performed to evaluate the erythroid differentiation of KO clones induced by hemin. Results: The high-efficiency sgRNA was successfully obtained and used for CRISPR-Cas9 editing in K562 cells. After limiting dilution and screening, two KO clones had either deleted 2 or 4 bases and showed no expression of the Vel protein. In the hemin-induced KO clone, there was a significant difference in erythroid marker and hemoglobin quantification compared to untreated cells. The morphological changes were also observed for the hemin-induced KO clone. Conclusion: In this study, a highly efficient sgRNA was screened out and used to generate Vel erythroid antigen KO single-cell clones in K562 cells. The edited cells could then be induced to undergo erythroid differentiation with the use of hemin. [ABSTRACT FROM AUTHOR]

Objective: To explore the role and mechanism of IL-39 in K562 cells. Methods: The expressions of IL-39R and STAT1/STAT3 in K562 cells were detected via qRT-PCR； the protein levels of phospho-STAT1/STAT3 and STAT1/STAT3 in K562 cells were detected by Western blot； the mRNA regulated by IL-39 in K562 cells was predicted by RNA sequencing. Results: rIL-39 significantly promoted the expression of IL-39R in K562 cells； IL-39 affected K562 cells through STAT1 pathway； IL-39 significantly regulated the expression of mRNA in K562 cells， thus potentially affecting a series of biological processes. Conclusion: IL-39 can directly act on K562 cells， alter gene transcriptional expression and inhibit tumor growth through the STAT1 pathway. IL-39 may similarly modulate gene expression profile in human leukemia cells. [ABSTRACT FROM AUTHOR]

Hydroquinone (HQ), one of the phenolic metabolites of benzene, is widely recognized as an important participant in benzene-induced hematotoxicity. However, there are few relevant proteomics in HQ-induced hematotoxicity and the mechanism hasn't been fully understood yet. In this study, we treated K562 cells with 40 μmol/L HQ for 72 h, examined and validated protein expression changes by Label-free proteomic analysis and Parallel reaction monitoring (PRM), and performed bioinformatics analysis to identify interaction networks. One hundred and eighty-seven upregulated differentially expressed proteins (DEPs) and 279 downregulated DEPs were identified in HQ-exposed K562 cells, which were involved in neutrophil-mediated immunity, blood microparticle, and other GO terms, as well as the lysosome, metabolic, cell cycle, and cellular senescence-related pathways. Focusing on the 23 DEGs and 5 DEPs in erythroid differentiation-related pathways, we constructed the network of protein interactions and determined 6 DEPs (STAT1, STAT3, CASP3, KIT, STAT5B, and VEGFA) as main hub proteins with the most interactions, among which STATs made a central impact and may be potential biomarkers of HQ-induced hematotoxicity. Our work reinforced the use of proteomics and bioinformatic approaches to advance knowledge on molecular mechanisms of HQ-induced hematotoxicity at the protein level and provide a valuable basis for further clarification. [ABSTRACT FROM AUTHOR]

Dinitrosyl iron complexes (DNICs) stabilize nitric oxide in cells and tissues and constitute an important form of its storage and transportation. DNICs may comprise low-molecular-weight ligands, e.g., thiols, imidazole groups in chemical compounds with low molecular weight (LMWDNICs), or high-molecular-weight ligands, e.g., peptides or proteins (HMWDNICs). The aim of this study was to investigate the role of low- and high-molecular-weight ligands in DNIC formation. Lysosomal and proteasomal proteolysis was inhibited by specific inhibitors. Experiments were conducted on human erythroid K562 cells and on K562 cells overexpressing a heavy chain of ferritin. Cell cultures were treated with •NO donor. DNIC formation was monitored by electron paramagnetic resonance. Pretreatment of cells with proteolysis inhibitors diminished the intensity and changed the shape of the DNIC-specific EPR signal in a treatment time-dependent manner. The level of DNIC formation was significantly influenced by the presence of protein degradation products. Interestingly, formation of HMWDNICs depended on the availability of LMWDNICs. The extent of glutathione involvement in the in vivo formation of DNICs is minor yet noticeable, aligning with our prior research findings. [ABSTRACT FROM AUTHOR]

Single nucleotide polymorphisms (SNPs) of HLA‐E are related to the occurrence of many diseases, but their functions remain unclear. In this study, the function of SNPs at HLA‐E rs76971248 and rs1264457 on the myeloid leukemia cells was analyzed by a progressive procedure, included genotyping, mRNA transcription, regulatory element, protein expression, and anti‐tumor effect. The frequencies of rs76971248 G and rs1264457 G were found higher in myeloid leukemia patients than those in healthy blood donors (p < 0.05). For myeloid leukemia, rs76971248 T was protective, while rs1264457 G was susceptible. We also found that rs76971248 affected HLA‐E mRNA transcription and membrane HLA‐E (mHLA‐E) expression in K562 cells through differently binding to transcription factor HOXA5 (p < 0.0001), while rs1264457 affected mHLA‐E expression by changing mRNA transcription and an encoding amino acid (p < 0.01). In contrast, the expression of soluble HLA‐E (sHLA‐E) was not influenced by both rs1264457 and rs76971248. The higher HLA‐E expression was detected among myeloid leukemia patients, and the K562 cells with higher HLA‐E molecules played a significant inhibitory effect on the killing activity of NK‐92MI cells (p < 0.05). In conclusion, the higher HLA‐E expression of myeloid leukemia cells is promoted by rs76971248 G and rs1264457 G, which helps escape from NK‐92MI cells' killing. [ABSTRACT FROM AUTHOR]

Oncological diseases represent a significant global health challenge, with high mortality rates. Early detection is crucial for effective treatment, and aptamers, which demonstrate superior specificity and stability compared to antibodies, offer a promising avenue for diagnostic advancement. This study presents the design, development and evaluation of a quartz crystal microbalance (QCM) sensor functionalized with the T2-KK1B10 aptamer for the sensitive and specific detection of Chronic Myeloid Leukemia (CML) K562 cells. The research focuses on optimizing the biorecognition layer by adjusting the aptamer conditions, demonstrating the sensor's ability to detect these CML cells with high specificity and sensitivity. The aptamer-modified QCM sensor operates on the principle of mass change detection upon binding of target cells. By employing the Langmuir isotherm model, the performance of the sensor was optimized for the capture of CML cells from biological samples with LOD of 263 K562 cells. The sensor was also successfully regenerated multiple times without sensitivity loss. Validation of the sensor's performance was conducted under controlled laboratory settings, followed by extensive testing utilizing human lyophilized plasma and clinical samples from patients. The sensor exhibited high sensitivity and specificity in the detection of CML cells within clinical specimens, thereby illustrating its potential for practical clinical deployment. This research presents a novel approach to the early diagnosis of CML, facilitating timely intervention and enhanced patient outcomes. The developed aptasensor demonstrates potential for broader application in cancer diagnostics and personalized medicine.

Background: Imatinib resistance remains a major obstacle in the treatment of chronic myelogenous leukemia (CML). Crocin (CRC) and astaxanthin (ATX) are phytochemicals with anti-cancer properties. Aims: This study aimed to explore the effects of combination treatment of Imatinib with CRC and ATX on Imatinib-resistant K562 (IR-K562) cells. Methods and results: After the establishment of IR-K562 cells, growth inhibitory activity was determined by the MTT assay. To test the regeneration potential, a colony formation assay was performed. Cell cycle analyses were examined by flow cytometry. Cell injury was evaluated by lactate dehydrogenase (LDH) leakage. Real-time PCR was applied to assess the expression of IL6, TNF-α, STAT3, BAD, CASP3, TP53, and Bcl-2 genes. Caspase-3 activity was determined by a colorimetric assay. Antioxidant activity was measured using a diphenylpicrylhydrazyl (DPPH) assay. After 48 h of treatment, ATX (IC50 = 30µM) and CRC (IC50 = 190µM) significantly inhibited cell proliferation and colony formation ability, induced G1 cell cycle arrest and cell injury, upregulated the expression of apoptosis-associated genes, and downregulated the expression of anti-apoptotic and inflammatory genes. The combination of IM with ATX and/or CRC synergistically reduced cell viability (combination index [CI] < 1). Conclusion: Our data suggest that IM shows better therapeutic efficacy at lower doses when combined with ATX and/or CRC. [ABSTRACT FROM AUTHOR]

To electrochemically measure human myeloid leukemia cells (K562 cells), we constructed a probe consisting of peptide/single‐strand (ss) DNA. Ac‐H6Y4C with an acetylated N‐terminal of peptide was used to enhance the probe to allow electrode responses that could detect target cells. A ss‐DNA was selected as the target cell recognition moiety. The probe exhibits properties that combine the functionalities of both DNA and peptides. The measurement principle is based on changes in the peak currents of the peptide moieties that are caused by interactions between the ss‐DNA and target cells. The peak currents were proportional to the concentration of K 562 cells that ranged from 10 to 2,000 cells/mL with a LOD of 3 cells/mL. [ABSTRACT FROM AUTHOR]

BACKGROUND: Natural killer (NK) cell-based immunotherapy is a promising treatment approach for multiple myeloma (MM), but obtaining a sufficient number of activated NK cells remains challenging. Here, we report an improved method to generate ex vivo expanded NK (eNK) cells from MM patients based on genetic engineering of K562 cells to express OX40 ligand and membrane-bound (mb) IL-18 and IL-21. METHODS: K562-OX40L-mbIL-18/-21 cells were generated by transducing K562-OX40L cells with a lentiviral vector encoding mbIL-18 and mbIL-21, and these were used as feeder cells to expand NK cells from peripheral blood mononuclear cells of healthy donors (HDs) and MM patients in the presence of IL-2/IL-15. Purity, expansion rate, receptor expression, and functions of eNK cells were determined over four weeks of culture. RESULTS: NK cell expansion was enhanced by short exposure of soluble IL-18 and IL-21 with K562-OX40L cells. Co-culture of NK cells with K562-OX40L-mbIL-18/-21 cells resulted in remarkable expansion of NK cells from HDs (9,860-fold) and MM patients (4,929-fold) over the 28-day culture period. Moreover, eNK cells showed increased expression of major activation markers and enhanced cytotoxicity towards target K562, U266, and RPMI8226 cells. CONCLUSIONS: Our data suggest that genetically engineered K562 cells expressing OX40L, mbIL-18, and mbIL-21 improve the expansion of NK cells, increase activation signals, and enhance their cytolytic activity towards MM cells.

Regulation of gene expression requires the combinatorial binding of sequence-specific transcription factors (TFs) at promoters and enhancers. Prior studies showed that alterations in the spacing between TF binding sites can influence promoter and enhancer activity. However, the relative importance of TF spacing alterations resulting from naturally occurring insertions and deletions (InDels) has not been systematically analyzed. To address this question, we first characterized the genome-wide spacing relationships of 73 TFs in human K562 cells as determined by ChIP-seq (chromatin immunoprecipitation sequencing). We found a dominant pattern of a relaxed range of spacing between collaborative factors, including 45 TFs exclusively exhibiting relaxed spacing with their binding partners. Next, we exploited millions of InDels provided by genetically diverse mouse strains and human individuals to investigate the effects of altered spacing on TF binding and local histone acetylation. These analyses suggested that spacing alterations resulting from naturally occurring InDels are generally tolerated in comparison to genetic variants directly affecting TF binding sites. To experimentally validate this prediction, we introduced synthetic spacing alterations between PU.1 and C/EBPβ binding sites at six endogenous genomic loci in a macrophage cell line. Remarkably, collaborative binding of PU.1 and C/EBPβ at these locations tolerated changes in spacing ranging from 5 bp increase to >30 bp decrease. Collectively, these findings have implications for understanding mechanisms underlying enhancer selection and for the interpretation of non-coding genetic variation.

Viral infection triggers activation of the integrated stress response (ISR). In response to viral double-stranded RNA (dsRNA), RNA-activated protein kinase (PKR) phosphorylates the translation initiation factor eIF2, converting it from a translation initiator into a potent translation inhibitor and this restricts the synthesis of viral proteins. Phosphorylated eIF2 (eIF2-P) inhibits translation by binding to eIF2's dedicated, heterodecameric nucleotide exchange factor eIF2B and conformationally inactivating it. We show that the NSs protein of Sandfly Fever Sicilian virus (SFSV) allows the virus to evade the ISR. Mechanistically, NSs tightly binds to eIF2B (KD = 30 nM), blocks eIF2-P binding, and rescues eIF2B GEF activity. Cryo-EM structures demonstrate that SFSV NSs and eIF2-P directly compete, with the primary NSs contacts to eIF2Bα mediated by five 'aromatic fingers'. NSs binding preserves eIF2B activity by maintaining eIF2B's conformation in its active A-State.

RNA molecules are known to fold into specific structures which often play a central role in their functions and regulation. In silico folding of RNA transcripts, especially when assisted with structure profiling (SP) data, is capable of accurately elucidating relevant structural conformations. However, such methods scale poorly to the swaths of SP data generated by transcriptome-wide experiments, which are becoming more commonplace and advancing our understanding of RNA structure and its regulation at global and local levels. This has created a need for tools capable of rapidly deriving structural assessments from SP data in a scalable manner. One such tool we previously introduced that aims to process such data is patteRNA, a statistical learning algorithm capable of rapidly mining big SP datasets for structural elements. Here, we present a reformulation of patteRNA's pattern recognition scheme that sees significantly improved precision without major compromises to computational overhead. Specifically, we developed a data-driven logistic classifier which interprets patteRNA's statistical characterizations of SP data in addition to local sequence properties as measured with a nearest neighbour thermodynamic model. Application of the classifier to human structurome data reveals a marked association between detected stem-loops and RNA binding protein (RBP) footprints. The results of our application demonstrate that upwards of 30% of RBP footprints occur within loops of stable stem-loop elements. Overall, our work arrives at a rapid and accurate method for automatically detecting families of RNA structure motifs and demonstrates the functional relevance of identifying them transcriptome-wide.

Dinitrosyl iron complexes (DNICs) stabilize nitric oxide in cells and tissues and constitute an important form of its storage and transportation. DNICs may comprise low-molecular-weight ligands, e.g., thiols, imidazole groups in chemical compounds with low molecular weight (LMWDNICs), or high-molecular-weight ligands, e.g., peptides or proteins (HMWDNICs). The aim of this study was to investigate the role of low- and high-molecular-weight ligands in DNIC formation. Lysosomal and proteasomal proteolysis was inhibited by specific inhibitors. Experiments were conducted on human erythroid K562 cells and on K562 cells overexpressing a heavy chain of ferritin. Cell cultures were treated with •NO donor. DNIC formation was monitored by electron paramagnetic resonance. Pretreatment of cells with proteolysis inhibitors diminished the intensity and changed the shape of the DNIC-specific EPR signal in a treatment time-dependent manner. The level of DNIC formation was significantly influenced by the presence of protein degradation products. Interestingly, formation of HMWDNICs depended on the availability of LMWDNICs. The extent of glutathione involvement in the in vivo formation of DNICs is minor yet noticeable, aligning with our prior research findings.

Dysregulated pre-mRNA splicing is an emerging Achilles heel of cancers and myelodysplasias. To expand the currently limited portfolio of small-molecule drug leads, we screened for chemical modulators of the U2AF complex, which nucleates spliceosome assembly and is mutated in myelodysplasias. A hit compound specifically enhances RNA binding by a U2AF2 subunit. Remarkably, the compound inhibits splicing of representative substrates and stalls spliceosome assembly at the stage of U2AF function. Computational docking, together with structure-guided mutagenesis, indicates that the compound bridges the tandem U2AF2 RNA recognition motifs via hydrophobic and electrostatic moieties. Cells expressing a cancer-associated U2AF1 mutant are preferentially killed by treatment with the compound. Altogether, our results highlight the potential of trapping early spliceosome assembly as an effective pharmacological means to manipulate pre-mRNA splicing. By extension, we suggest that stabilizing assembly intermediates may offer a useful approach for small-molecule inhibition of macromolecular machines.

Additional sex combs-like 1 (ASXL1) interacts with BRCA1-associated protein 1 (BAP1) deubiquitinase to oppose the polycomb repressive complex 1 (PRC1)-mediated histone H2A ubiquitylation. Germline BAP1 mutations are found in a spectrum of human malignancies, while ASXL1 mutations recurrently occur in myeloid neoplasm and are associated with poor prognosis. Nearly all ASXL1 mutations are heterozygous frameshift or nonsense mutations in the middle or to a less extent the C-terminal region, resulting in the production of C-terminally truncated mutant ASXL1 proteins. How ASXL1 regulates specific target genes and how the C-terminal truncation of ASXL1 promotes leukemogenesis are unclear. Here, we report that ASXL1 interacts with forkhead transcription factors FOXK1 and FOXK2 to regulate a subset of FOXK1/K2 target genes. We show that the C-terminally truncated mutant ASXL1 proteins are expressed at much higher levels than the wild-type protein in ASXL1 heterozygous leukemia cells, and lose the ability to interact with FOXK1/K2. Specific deletion of the mutant allele eliminates the expression of C-terminally truncated ASXL1 and increases the association of wild-type ASXL1 with BAP1, thereby restoring the expression of BAP1-ASXL1-FOXK1/K2 target genes, particularly those involved in glucose metabolism, oxygen sensing, and JAK-STAT3 signaling pathways. In addition to FOXK1/K2, we also identify other DNA-binding transcription regulators including transcription factors (TFs) which interact with wild-type ASXL1, but not C-terminally truncated mutant. Our results suggest that ASXL1 mutations result in neomorphic alleles that contribute to leukemogenesis at least in part through dominantly inhibiting the wild-type ASXL1 from interacting with BAP1 and thereby impairing the function of ASXL1-BAP1-TF in regulating target genes and leukemia cell growth.

Chronic myeloid leukemia (CML) is an uncommon type of white blood cells cancer that originates from bone marrow stem cells. Progesterone (P4) and berberine (BBR) are bioactive compounds that inhibit the growth of tumor cells. The present study aimed to assess the simultaneous effect of P4 and BBR on the level of reactive oxygen species (ROS) in K562 cells. The K562 cells were cultured in a complete cell culture medium and simultaneously exposed to IC50 different concentrations of P4 and BBR at 24 h (P4:102.4 μM; BBR:125 μM), 48 h (P4:78.4 μM; BBR:114 μM), and 72 h (P4:70 μM; BBR:45 μM). Then, the cell viability and cellular ROS level were determined using MTT assay and 2′, 7′-dichlorofluorescin diacetate (DCF) by flow cytometry, respectively. Our results showed that the combination of P4 and BBR inhibited the cells more effectively than P4 and BBR alone at 72 h, as well as P4 and their combination reduced ROS level in the cells compared to BBR and untreated cells at 24 h and 72 h. Taken together, P4 through a ROS-dependent pathway and BBR through a ROS-independent pathway may be inhibited cell growth.

Objective To delineate RNA binding profile of NAT10(N-acetyltransferase 10) in human chronic myelogenous leukemia cell line K562. Methods Enhanced UV crosslinking and immunoprecipitation (eCLIP) were used to capture the transcript collection bound by NAT10 in K562 cells. Data were obtained through high throughput sequencing and then analyzed with bioinformatics. The type and regions of NAT10 binding genes were identified via peak-calling and annotation. The function of the binding RNA was further explored by gene function enrichment analysis. Results NAT10 binding transcripts which were enriched in 3′UTR region, were mainly protein coding genes. Further analysis showed that they were functionally associated with DNA damage and repair. Conclusions NAT10 may regulate gene expression by binding to the 3′UTR region of mRNA associated with DNA damage repair.

Human erythroleukemic K562 cells represent the prototypical cell culture model of chronic myeloid leukemia (CML). The cells are pseudo-triploid and positive for the Philadelphia chromosome. Therefore, K562 cells have been widely used for investigating the BCR/ABL1 oncogene and the tyrosine kinase inhibitor, imatinib-mesylate. Further, K562 cells overexpress transferrin receptors (TfR) and have been used as a model for targeting cytotoxic therapies, via receptor-mediated endocytosis. Here, we have characterized K562 cells focusing on the karyotype of cells in prolonged culture, regulation of expression of TfR in wildtype (WT) and doxorubicin-resistant cells, and responses to histone deacetylase inhibition (HDACi). Karyotype analysis indicates novel chromosomes and gene expression analysis suggests a shift of cultured K562 cells away from patient-derived leukemic cells. We confirm the high expression of TfR on K562 cells using immunofluorescence and cell-surface receptor binding radioassays. Importantly, high TfR expression is observed in patient-derived cells, and we highlight the persistent expression of TfR following doxorubicin acquired resistance. Epigenetic analysis indicates that permissive histone acetylation and methylation at the promoter region regulates the transcription of TfR in K562 cells. Finally, we show relatively high expression of HDAC enzymes in K562 cells and demonstrate the chemotoxic effects of HDACi, using the FDA-approved hydroxamic acid, vorinostat. Together with a description of morphology, infrared spectral analysis, and examination of metabolic properties, we provide a comprehensive characterization of K562 cells. Overall, K562 cell culture systems remain widely used for the investigation of novel therapeutics for CML, which is particularly important in cases of imatinib-mesylate resistance. [ABSTRACT FROM AUTHOR]

Platelet aggregation at the site of atherosclerotic vascular injury is the underlying pathophysiology of myocardial infarction and stroke. To build upon prior GWAS, here we report on 16 loci identified through a whole genome sequencing (WGS) approach in 3,855 NHLBI Trans-Omics for Precision Medicine (TOPMed) participants deeply phenotyped for platelet aggregation. We identify the RGS18 locus, which encodes a myeloerythroid lineage-specific regulator of G-protein signaling that co-localizes with expression quantitative trait loci (eQTL) signatures for RGS18 expression in platelets. Gene-based approaches implicate the SVEP1 gene, a known contributor of coronary artery disease risk. Sentinel variants at RGS18 and PEAR1 are associated with thrombosis risk and increased gastrointestinal bleeding risk, respectively. Our WGS findings add to previously identified GWAS loci, provide insights regarding the mechanism(s) by which genetics may influence cardiovascular disease risk, and underscore the importance of rare variant and regulatory approaches to identifying loci contributing to complex phenotypes.

X-linked agammaglobulinemia (XLA) is a monogenic primary immune deficiency characterized by very low levels of immunoglobulins and greatly increased risks for recurrent and severe infections. Patients with XLA have a loss-of-function mutation in the Bruton's tyrosine kinase (BTK) gene and fail to produce mature B lymphocytes. Gene editing in the hematopoietic stem cells of XLA patients to correct or replace the defective gene should restore B cell development and the humoral immune response. We used the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 platform to precisely target integration of a corrective, codon-optimized BTK complementary DNA (cDNA) cassette into its endogenous locus. This process is driven by homologous recombination and should place the transgenic BTK under transcriptional control of its endogenous regulatory elements. Each integrated copy of this cDNA in BTK-deficient K562 cells produced only 11% as much BTK protein as the wild-type gene. The donor cDNA was modified to include the terminal intron of the BTK gene. Successful integration of the intron-containing BTK donor led to a nearly twofold increase in BTK expression per cell over the base donor. However, this donor variant was too large to package into an adeno-associated viral vector for delivery into primary cells. Donors containing truncated variants of the terminal intron also produced elevated expression, although to a lesser degree than the full intron. Addition of the Woodchuck hepatitis virus posttranscriptional regulatory element led to a large boost in BTK transgene expression. Combining these modifications led to a BTK donor template that generated nearly physiological levels of BTK expression in cell lines. These reagents were then optimized to maximize integration rates into human hematopoietic stem and progenitor cells, which have reached potentially therapeutic levels in vitro. The novel donor modifications support effective gene therapy for XLA and will likely assist in the development of other gene editing-based therapies for genetic disorders.

Overexpressed tumor-associated antigens [for example, epidermal growth factor receptor (EGFR) and human epidermal growth factor receptor 2 (HER2)] are attractive targets for therapeutic T cells, but toxic "off-tumor" cross-reaction with normal tissues that express low levels of target antigen can occur with chimeric antigen receptor (CAR)-T cells. Inspired by natural ultrasensitive response circuits, we engineered a two-step positive-feedback circuit that allows human cytotoxic T cells to discriminate targets on the basis of a sigmoidal antigen-density threshold. In this circuit, a low-affinity synthetic Notch receptor for HER2 controls the expression of a high-affinity CAR for HER2. Increasing HER2 density thus has cooperative effects on T cells-it increases both CAR expression and activation-leading to a sigmoidal response. T cells with this circuit show sharp discrimination between target cells expressing normal amounts of HER2 and cancer cells expressing 100 times as much HER2, both in vitro and in vivo.

T cells experience complex temporal patterns of stimulus via receptor-ligand-binding interactions with surrounding cells. From these temporal patterns, T cells are able to pick out antigenic signals while establishing self-tolerance. Although features such as duration of antigen binding have been examined, our understanding of how T cells interpret signals with different frequencies or temporal stimulation patterns is relatively unexplored. We engineered T cells to respond to light as a stimulus by building an optogenetically controlled chimeric antigen receptor (optoCAR). We discovered that T cells respond to minute-scale oscillations of activation signal by stimulating optoCAR T cells with tunable pulse trains of light. Systematically scanning signal oscillation period from 1 to 150 min revealed that expression of CD69, a T cell activation marker, reached a local minimum at a period of ∼25 min (corresponding to 5 to 15 min pulse widths). A combination of inhibitors and genetic knockouts suggest that this frequency filtering mechanism lies downstream of the Erk signaling branch of the T cell response network and may involve a negative feedback loop that diminishes Erk activity. The timescale of CD69 filtering corresponds with the duration of T cell encounters with self-peptide-presenting APCs observed via intravital imaging in mice, indicating a potential functional role for temporal filtering in vivo. This study illustrates that the T cell signaling machinery is tuned to temporally filter and interpret time-variant input signals in discriminatory ways.

Recent technological advances have expanded the annotated protein coding content of mammalian genomes, as hundreds of previously unidentified, short open reading frame (ORF)-encoded peptides (SEPs) have now been found to be translated. Although several studies have identified important physiological roles for this emerging protein class, a general method to define their interactomes is lacking. Here, we demonstrate that genetic incorporation of the photo-crosslinking noncanonical amino acid AbK into SEP transgenes allows for the facile identification of SEP cellular interaction partners using affinity-based methods. From a survey of seven SEPs, we report the discovery of short ORF-encoded histone binding protein (SEHBP), a conserved microprotein that interacts with chromatin-associated proteins, localizes to discrete genomic loci, and induces a robust transcriptional program when overexpressed in human cells. This work affords a straightforward method to help define the physiological roles of SEPs and demonstrates its utility by identifying SEHBP as a short ORF-encoded transcription factor.

DNA strand asymmetries can have a major effect on several biological functions, including replication, transcription and transcription factor binding. As such, DNA strand asymmetries and mutational strand bias can provide information about biological function. However, a versatile tool to explore this does not exist. Here, we present Asymmetron, a user-friendly computational tool that performs statistical analysis and visualizations for the evaluation of strand asymmetries. Asymmetron takes as input DNA features provided with strand annotation and outputs strand asymmetries for consecutive occurrences of a single DNA feature or between pairs of features. We illustrate the use of Asymmetron by identifying transcriptional and replicative strand asymmetries of germline structural variant breakpoints. We also show that the orientation of the binding sites of 45% of human transcription factors analyzed have a significant DNA strand bias in transcribed regions, that is also corroborated in ChIP-seq analyses, and is likely associated with transcription. In summary, we provide a novel tool to assess DNA strand asymmetries and show how it can be used to derive new insights across a variety of biological disciplines.

Dysfunction of the mitochondrial electron transport chain (mETC) is a major cause of human mitochondrial diseases. To identify determinants of mETC function, we screened a genome-wide human CRISPRi library under oxidative metabolic conditions with selective inhibition of mitochondrial Complex III and identified ovarian carcinoma immunoreactive antigen (OCIA) domain-containing protein 1 (OCIAD1) as a Complex III assembly factor. We find that OCIAD1 is an inner mitochondrial membrane protein that forms a complex with supramolecular prohibitin assemblies. Our data indicate that OCIAD1 is required for maintenance of normal steady-state levels of Complex III and the proteolytic processing of the catalytic subunit cytochrome c1 (CYC1). In OCIAD1 depleted mitochondria, unprocessed CYC1 is hemylated and incorporated into Complex III. We propose that OCIAD1 acts as an adaptor within prohibitin assemblies to stabilize and/or chaperone CYC1 and to facilitate its proteolytic processing by the IMMP2L protease.

Natural killer (NK) cells are key effectors of the innate immune system, but major differences between human and murine NK cells have impeded translation. Outbred dogs offer an important link for studies of NK biology and immunotherapy. We analyzed gene expression of putative NK populations from healthy dogs and dogs with naturally-occurring cancers examining differential gene expression across multiple conditions, including steady-state, in vitro activation with cytokines and co-culture, and in vivo activation with inhaled IL-15 in dogs receiving IL-15 immunotherapy. We also compared dog, mouse and human CD3-NKp46+ NK cells using a novel orthologous transcriptome. Distinct transcriptional profiles between NK populations exist between conditions and in vitro versus in vivo treatments. In cross-species analysis, canine NK cells were globally more similar to human NK cells than mice. These data define canine NK cell gene expression under multiple conditions and across species, filling an important gap in translational NK studies.

Our understanding of the beads-on-a-string arrangement of nucleosomes has been built largely on high-resolution sequence-agnostic imaging methods and sequence-resolved bulk biochemical techniques. To bridge the divide between these approaches, we present the single-molecule adenine methylated oligonucleosome sequencing assay (SAMOSA). SAMOSA is a high-throughput single-molecule sequencing method that combines adenine methyltransferase footprinting and single-molecule real-time DNA sequencing to natively and nondestructively measure nucleosome positions on individual chromatin fibres. SAMOSA data allows unbiased classification of single-molecular 'states' of nucleosome occupancy on individual chromatin fibres. We leverage this to estimate nucleosome regularity and spacing on single chromatin fibres genome-wide, at predicted transcription factor binding motifs, and across human epigenomic domains. Our analyses suggest that chromatin is comprised of both regular and irregular single-molecular oligonucleosome patterns that differ subtly in their relative abundance across epigenomic domains. This irregularity is particularly striking in constitutive heterochromatin, which has typically been viewed as a conformationally static entity. Our proof-of-concept study provides a powerful new methodology for studying nucleosome organization at a previously intractable resolution and offers up new avenues for modeling and visualizing higher order chromatin structure.

BackgroundCircular RNAs (circRNAs) are stable, often highly expressed RNA transcripts with potential to modulate other regulatory RNAs. A few circRNAs have been shown to bind RNA-binding proteins (RBPs); however, little is known about the prevalence and distribution of these interactions in different biological contexts.MethodsWe conduct an extensive screen of circRNA-RBP interactions in the ENCODE cell lines HepG2 and K562. We profile circRNAs in deep-sequenced total RNA samples and analyze circRNA-RBP interactions using a large set of eCLIP data with binding sites of 150 RBPs. We validate interactions for select circRNAs and RBPs by performing RNA immunoprecipitation and functionally characterize our most interesting candidates by conducting knockdown studies followed by RNA-Seq.ResultsWe generate a comprehensive catalog of circRNA-RBP interactions in HepG2 and K562 cells. We show that KHSRP binding sites are enriched in flanking introns of circRNAs and that KHSRP depletion affects circRNA biogenesis. We identify circRNAs that are highly covered by RBP binding sites and experimentally validate individual circRNA-RBP interactions. We show that circCDYL, a highly expressed circRNA with clinical and functional implications in bladder cancer, is almost completely covered with GRWD1 binding sites in HepG2 cells, and that circCDYL depletion counteracts the effect of GRWD1 depletion. Furthermore, we confirm interactions between circCDYL and RBPs in bladder cancer cells and demonstrate that circCDYL depletion affects hallmarks of cancer and perturbs the expression of key cancer genes, e.g., TP53. Finally, we show that elevated levels of circCDYL are associated with overall survival of bladder cancer patients.ConclusionsOur study demonstrates transcriptome-wide and cell-type-specific circRNA-RBP interactions that could play important regulatory roles in tumorigenesis.

Impaired DNA crosslink repair leads to Fanconi anemia (FA), characterized by a unique manifestation of bone marrow failure and pancytopenia among diseases caused by DNA damage response defects. As a germline disorder, why the hematopoietic hierarchy is specifically affected is not fully understood. We find that reprogramming transcription during hematopoietic differentiation results in an overload of genotoxic stress, which causes aborted differentiation and depletion of FA mutant progenitor cells. DNA damage onset most likely arises from formaldehyde, an obligate by-product of oxidative protein demethylation during transcription regulation. Our results demonstrate that rapid and extensive transcription reprogramming associated with hematopoietic differentiation poses a major threat to genome stability and cell viability in the absence of the FA pathway. The connection between differentiation and DNA damage accumulation reveals a novel mechanism of genome scarring and is critical to exploring therapies to counteract the aplastic anemia for the treatment of FA patients.

BackgroundA critical step in uncovering rules of RNA processing is to study the in vivo regulatory networks of RNA binding proteins (RBPs). Crosslinking and immunoprecipitation (CLIP) methods enable mapping RBP targets transcriptome-wide, but methodological differences present challenges to large-scale analysis across datasets. The development of enhanced CLIP (eCLIP) enabled the mapping of targets for 150 RBPs in K562 and HepG2, creating a unique resource of RBP interactomes profiled with a standardized methodology in the same cell types.ResultsOur analysis of 223 eCLIP datasets reveals a range of binding modalities, including highly resolved positioning around splicing signals and mRNA untranslated regions that associate with distinct RBP functions. Quantification of enrichment for repetitive and abundant multicopy elements reveals 70% of RBPs have enrichment for non-mRNA element classes, enables identification of novel ribosomal RNA processing factors and sites, and suggests that association with retrotransposable elements reflects multiple RBP mechanisms of action. Analysis of spliceosomal RBPs indicates that eCLIP resolves AQR association after intronic lariat formation, enabling identification of branch points with single-nucleotide resolution, and provides genome-wide validation for a branch point-based scanning model for 3' splice site recognition. Finally, we show that eCLIP peak co-occurrences across RBPs enable the discovery of novel co-interacting RBPs.ConclusionsThis work reveals novel insights into RNA biology by integrated analysis of eCLIP profiling of 150 RBPs with distinct functions. Further, our quantification of both mRNA and other element association will enable further research to identify novel roles of RBPs in regulating RNA processing.

BACKGROUND:During transcription, numerous transcription factors (TFs) bind to targets in a highly coordinated manner to control the gene expression. Alterations in groups of TF-binding profiles (i.e. "co-binding changes") can affect the co-regulating associations between TFs (i.e. "rewiring the co-regulator network"). This, in turn, can potentially drive downstream expression changes, phenotypic variation, and even disease. However, quantification of co-regulatory network rewiring has not been comprehensively studied. RESULTS:To address this, we propose DiNeR, a computational method to directly construct a differential TF co-regulation network from paired disease-to-normal ChIP-seq data. Specifically, DiNeR uses a graphical model to capture the gained and lost edges in the co-regulation network. Then, it adopts a stability-based, sparsity-tuning criterion -- by sub-sampling the complete binding profiles to remove spurious edges -- to report only significant co-regulation alterations. Finally, DiNeR highlights hubs in the resultant differential network as key TFs associated with disease. We assembled genome-wide binding profiles of 104 TFs in the K562 and GM12878 cell lines, which loosely model the transition between normal and cancerous states in chronic myeloid leukemia (CML). In total, we identified 351 significantly altered TF co-regulation pairs. In particular, we found that the co-binding of the tumor suppressor BRCA1 and RNA polymerase II, a well-known transcriptional pair in healthy cells, was disrupted in tumors. Thus, DiNeR successfully extracted hub regulators and discovered well-known risk genes. CONCLUSIONS:Our method DiNeR makes it possible to quantify changes in co-regulatory networks and identify alterations to TF co-binding patterns, highlighting key disease regulators. Our method DiNeR makes it possible to quantify changes in co-regulatory networks and identify alterations to TF co-binding patterns, highlighting key disease regulators.

Gene expression variability, differences in the number of mRNA per cell across a population of cells, is ubiquitous across diverse organisms with broad impacts on cellular phenotypes. The role of chromatin in regulating average gene expression has been extensively studied. However, what aspects of the chromatin contribute to gene expression variability is still underexplored. Here we addressed this problem by leveraging chromatin diversity and using a systematic investigation of randomly integrated expression reporters to identify what aspects of chromatin microenvironment contribute to gene expression variability. Using DNA barcoding and split-pool decoding, we created a large library of isogenic reporter clones and identified reporter integration sites in a massive and parallel manner. By mapping our measurements of reporter expression at different genomic loci with multiple epigenetic profiles including the enrichment of transcription factors and the distance to different chromatin states, we identified new factors that impact the regulation of gene expression distributions.

Background and Objective: Today, cancer is one of the most important causes of death in the world, and so far, many treatment methods have been used in this field. Immunotherapy is considered one of the newest developments in this science, and it is still being investigated in some forms in different cancers and with a variety of antigens as well. One of the subsets of cancer immunotherapy is its treatment using parasitic antigens. The present study evaluated the effect of using somatic antigens of protoscoleces of Echinococcus granulosus on K562 cancer cells. Methods: In this study, hydatid cysts' protoscolex antigens were extracted, purified, and added to K562 cancer cells at three concentrations (0.1, 1, and 2 mg/ml) and on three times (24, 48, and 72 h). The number of apoptotic cells was compared to the control flask. The antigen concentration of 2 mg/ml was used as a control sample to investigate its cytotoxic effect on the growth of healthy HFF3 cells. Annexin V and PI tests were also performed to differentiate apoptosis from necrosis. Results: In flasks treated with hydatid cyst protoscolex antigen, all three concentrations significantly reduced the growth of cancer cells compared with the control flask, and concentration 2 of crude antigen significantly caused the death of cancer cells. Furthermore, more cancer cells underwent apoptosis by increasing the time of exposure to the antigen. On the other hand, flow cytometry results also showed that the amount of apoptosis has increased compared to the control group. In fact, Protoscolex hydatid cyst somatic antigens induce programmed cell death in K562 cancer cells while not having a cytotoxic effect on normal cells. Conclusion: Therefore, it is suggested to do more research on the anti-cancer and therapeutic properties of the antigens of this parasite. [ABSTRACT FROM AUTHOR]

Objective To delineate RNA binding profile of NAT10 (N-acetyltransferase 10) in human chronic myelogenous leukemia cell line K562. Methods Enhanced UV crosslinking and immunoprecipitation (eCLIP) were used to capture the transcript collection bound by NAT10 in K562 cells. Data were obtained through high throughput sequencing and then analyzed with bioinformatics. The type and regions of NAT10 binding genes were identified via peak-calling and annotation. The function of the binding RNA was further explored by gene function enrichment analysis. Results NAT10 binding transcripts which were enriched in 3′UTR region, were mainly protein coding genes. Further analysis showed that they were functionally associated with DNA damage and repair. Conclusions NAT10 may regulate gene expression by binding to the 3′UTR region of mRNA associated with DNA damage repair. [ABSTRACT FROM AUTHOR]

Nowadays, the search for food products that promote consumers' health has gained interest, and dairy by-products, due to their biological quality, could have a prominent position among products with health benefits. However, little is known about their activity on cancer cells. This study aimed to provide evidence about the effect of ovine colostrum and milk whey on K562 cells, a model of the human chronic myeloid leukemia cell line. The exposure of K562 cells to a single administration of sheep by-products at different concentrations for three days and three treatments for three days was carried out. Using a flow cytometric approach, we found that CD235a expression remained stable in the cells exposed to ovine whey (milk and colostrum) at concentrations ranging from 1 ng/mL to 100 μg/mL, after three days from one or three administrations, respectively. A significant reduction in fluorescent cells was observed in the populations exposed to 1 mg/mL of both milk and colostrum at the same time points. In these conditions, the size and granularity of the leukemic cells also changed, with a substantial reduction in the number of actively dividing cells in the S phase of the cell cycle. This phenomenon was highlighted by the Annexin V/PI cytofluorimetric test, which is able to provide quantitative results regarding the population of cells in early or late apoptosis or necrotic cells after exposure to a single dose or three doses of colostrum or sheep whey for three days, respectively. This report showed that both colostrum and milk whey were able to modify the phenotypic profile and cell cycle of the K562 cell line, inducing apoptosis at the highest concentration. [ABSTRACT FROM AUTHOR]

Quiescence is a hallmark of CD4+ T cells latently infected with human immunodeficiency virus 1 (HIV-1). While reversing this quiescence is an effective approach to reactivate latent HIV from T cells in culture, it can cause deleterious cytokine dysregulation in patients. As a key regulator of T-cell quiescence, FOXO1 promotes latency and suppresses productive HIV infection. We report that, in resting T cells, FOXO1 inhibition impaired autophagy and induced endoplasmic reticulum (ER) stress, thereby activating two associated transcription factors: activating transcription factor 4 (ATF4) and nuclear factor of activated T cells (NFAT). Both factors associate with HIV chromatin and are necessary for HIV reactivation. Indeed, inhibition of protein kinase R-like ER kinase, an ER stress sensor that can mediate the induction of ATF4, and calcineurin, a calcium-dependent regulator of NFAT, synergistically suppressed HIV reactivation induced by FOXO1 inhibition. Thus, our studies uncover a link of FOXO1, ER stress and HIV infection that could be therapeutically exploited to selectively reverse T-cell quiescence and reduce the size of the latent viral reservoir.

Disrupted energy metabolism drives cell dysfunction and disease, but approaches to increase or preserve ATP are lacking. To generate a comprehensive metabolic map of genes and pathways that regulate cellular ATP-the ATPome-we conducted a genome-wide CRISPR interference/activation screen integrated with an ATP biosensor. We show that ATP level is modulated by distinct mechanisms that promote energy production or inhibit consumption. In our system HK2 is the greatest ATP consumer, indicating energy failure may not be a general deficiency in producing ATP, but rather failure to recoup the ATP cost of glycolysis and diversion of glucose metabolites to the pentose phosphate pathway. We identify systems-level reciprocal inhibition between the HIF1 pathway and mitochondria; glycolysis-promoting enzymes inhibit respiration even when there is no glycolytic ATP production, and vice versa. Consequently, suppressing alternative metabolism modes paradoxically increases energy levels under substrate restriction. This work reveals mechanisms of metabolic control, and identifies therapeutic targets to correct energy failure.

Sprouty-related, EVH1 domain-containing (SPRED) proteins negatively regulate RAS/mitogen-activated protein kinase (MAPK) signaling following growth factor stimulation. This inhibition of RAS is thought to occur primarily through SPRED1 binding and recruitment of neurofibromin, a RasGAP, to the plasma membrane. Here, we report the structure of neurofibromin (GTPase-activating protein [GAP]-related domain) complexed with SPRED1 (EVH1 domain) and KRAS. The structure provides insight into how the membrane targeting of neurofibromin by SPRED1 allows simultaneous interaction with activated KRAS. SPRED1 and NF1 loss-of-function mutations occur across multiple cancer types and developmental diseases. Analysis of the neurofibromin-SPRED1 interface provides a rationale for mutations observed in Legius syndrome and suggests why SPRED1 can bind to neurofibromin but no other RasGAPs. We show that oncogenic EGFR(L858R) signaling leads to the phosphorylation of SPRED1 on serine 105, disrupting the SPRED1-neurofibromin complex. The structural, biochemical, and biological results provide new mechanistic insights about how SPRED1 interacts with neurofibromin and regulates active KRAS levels in normal and pathologic conditions.

Therapeutic targeting of epigenetic modulators offers a novel approach to the treatment of multiple diseases. The cellular consequences of chemical compounds that target epigenetic regulators (epi-drugs) are complex. Epi-drugs affect global cellular phenotypes and cause local changes to gene expression due to alteration of a gene chromatin environment. Despite increasing use in the clinic, the mechanisms responsible for cellular changes are unclear. Specifically, to what degree the effects are a result of cell-wide changes or disease related locus specific effects is unknown. Here we developed a platform to systematically and simultaneously investigate the sensitivity of epi-drugs at hundreds of genomic locations by combining DNA barcoding, unique split-pool encoding, and single cell expression measurements. Internal controls are used to isolate locus specific effects separately from any global consequences these drugs have. Using this platform we discovered wide-spread loci specific sensitivities to epi-drugs for three distinct epi-drugs that target histone deacetylase, DNA methylation and bromodomain proteins. By leveraging ENCODE data on chromatin modification, we identified features of chromatin environments that are most likely to be affected by epi-drugs. The measurements of loci specific epi-drugs sensitivities will pave the way to the development of targeted therapy for personalized medicine.

Chimeric antigen receptor (CAR)-T cell therapy has shown remarkable clinical efficacy against B-cell malignancies, yet marked vulnerability to antigen escape and tumor relapse exists. Here we report the rational design and optimization of bispecific CAR-T cells with robust activity against heterogeneous multiple myeloma (MM) that is resistant to conventional CAR-T cell therapy targeting B-cell maturation antigen (BCMA). We demonstrate that BCMA/CS1 bispecific CAR-T cells exhibit superior CAR expression and function compared to T cells that co-express individual BCMA and CS1 CARs. Combination therapy with anti-PD-1 antibody further accelerates the rate of initial tumor clearance in vivo, while CAR-T cell treatment alone achieves durable tumor-free survival even upon tumor re-challenge. Taken together, the BCMA/CS1 bispecific CAR presents a promising treatment approach to prevent antigen escape in CAR-T cell therapy against MM, and the vertically integrated optimization process can be used to develop robust cell-based therapy against novel disease targets.

DISCOVER-seq (discovery of in situ Cas off-targets and verification by sequencing) is a broadly applicable approach for unbiased CRISPR-Cas off-target identification in cells and tissues. It leverages the recruitment of DNA repair factors to double-strand breaks (DSBs) after genome editing with CRISPR nucleases. Here, we describe a detailed experimental protocol and analysis pipeline with which to perform DISCOVER-seq. The principle of this method is to track the precise recruitment of MRE11 to DSBs by chromatin immunoprecipitation followed by next-generation sequencing. A customized open-source bioinformatics pipeline, BLENDER (blunt end finder), then identifies off-target sequences genome wide. DISCOVER-seq is capable of finding and measuring off-targets in primary cells and in situ. The two main advantages of DISCOVER-seq are (i) low false-positive rates because DNA repair enzyme binding is required for genome edits to occur and (ii) its applicability to a wide variety of systems, including patient-derived cells and animal models. The whole protocol, including the analysis, can be completed within 2 weeks.

Repair of double strand DNA breaks (DSBs) can result in gene disruption or gene modification via homology directed repair (HDR) from donor DNA. Altering cellular responses to DSBs may rebalance editing outcomes towards HDR and away from other repair outcomes. Here, we utilize a pooled CRISPR screen to define host cell involvement in HDR between a Cas9 DSB and a plasmid double stranded donor DNA (dsDonor). We find that the Fanconi Anemia (FA) pathway is required for dsDonor HDR and that other genes act to repress HDR. Small molecule inhibition of one of these repressors, CDC7, by XL413 and other inhibitors increases the efficiency of HDR by up to 3.5 fold in many contexts, including primary T cells. XL413 stimulates HDR during a reversible slowing of S-phase that is unexplored for Cas9-induced HDR. We anticipate that XL413 and other such rationally developed inhibitors will be useful tools for gene modification.

A lack of tools to precisely control gene expression has limited our ability to evaluate relationships between expression levels and phenotypes. Here, we describe an approach to titrate expression of human genes using CRISPR interference and series of single-guide RNAs (sgRNAs) with systematically modulated activities. We used large-scale measurements across multiple cell models to characterize activities of sgRNAs containing mismatches to their target sites and derived rules governing mismatched sgRNA activity using deep learning. These rules enabled us to synthesize a compact sgRNA library to titrate expression of ~2,400 genes essential for robust cell growth and to construct an in silico sgRNA library spanning the human genome. Staging cells along a continuum of gene expression levels combined with single-cell RNA-seq readout revealed sharp transitions in cellular behaviors at gene-specific expression thresholds. Our work provides a general tool to control gene expression, with applications ranging from tuning biochemical pathways to identifying suppressors for diseases of dysregulated gene expression.

Reactive oxygen species (ROS) play critical roles in metabolism and disease, yet a comprehensive analysis of the cellular response to oxidative stress is lacking. To systematically identify regulators of oxidative stress, we conducted genome-wide Cas9/CRISPR and shRNA screens. This revealed a detailed picture of diverse pathways that control oxidative stress response, ranging from the TCA cycle and DNA repair machineries to iron transport, trafficking, and metabolism. Paradoxically, disrupting the pentose phosphate pathway (PPP) at the level of phosphogluconate dehydrogenase (PGD) protects cells against ROS. This dramatically alters metabolites in the PPP, consistent with rewiring of upper glycolysis to promote antioxidant production. In addition, disruption of peroxisomal import unexpectedly increases resistance to oxidative stress by altering the localization of catalase. Together, these studies provide insights into the roles of peroxisomal matrix import and the PPP in redox biology and represent a rich resource for understanding the cellular response to oxidative stress.

The analysis of samples from unsequenced and/or understudied species as well as samples where the proteome is derived from multiple organisms poses two key questions. The first is whether the proteomic data obtained from an unusual sample type even contains peptide tandem mass spectra. The second question is whether an appropriate protein sequence database is available for proteomic searches. We describe the use of automated de novo sequencing for evaluating both the quality of a collection of tandem mass spectra and the suitability of a given protein sequence database for searching that data. Applications of this method include the proteome analysis of closely related species, metaproteomics, and proteomics of extinct organisms.