Mapping systemic lupus erythematosus heterogeneity at the single-cell level.

Patients with systemic lupus erythematosus (SLE) display a complex blood transcriptome whose cellular origin is poorly resolved. Using single-cell RNA sequencing, we profiled ~276,000 peripheral blood mononuclear cells from 33 children with SLE with different degrees of disease activity and 11 matched controls. Increased expression of interferon-stimulated genes (ISGs) distinguished cells from children with SLE from healthy control cells. The high ISG expression signature (ISGhi) derived from a small number of transcriptionally defined subpopulations within major cell types, including monocytes, CD4+ and CD8+ T cells, natural killer cells, conventional and plasmacytoid dendritic cells, B cells and especially plasma cells. Expansion of unique subpopulations enriched in ISGs and/or in monogenic lupus-associated genes classified patients with the highest disease activity. Profiling of ~82,000 single peripheral blood mononuclear cells from adults with SLE confirmed the expansion of similar subpopulations in patients with the highest disease activity. This study lays the groundwork for resolving the origin of the SLE transcriptional signatures and the disease heterogeneity towards precision medicine applications.

The Relationship Between Metabolic Syndrome and Mortality Among Patients With Acute Respiratory Distress Syndrome in Acute Respiratory Distress Syndrome Network and Prevention and Early Treatment of Acute Lung Injury Network Trials.

OBJECTIVES: Metabolic syndrome is known to predict outcomes in COVID-19 acute respiratory distress syndrome (ARDS) but has never been studied in non-COVID-19 ARDS. We therefore aimed to determine the association of metabolic syndrome with mortality among ARDS trial subjects. DESIGN: Retrospective cohort study of ARDS trials' data. SETTING: An ancillary analysis was conducted using data from seven ARDS Network and Prevention and Early Treatment of Acute Lung Injury Network randomized trials within the Biologic Specimen and Data Repository Information Coordinating Center database. PATIENTS: Hospitalized patients with ARDS and metabolic syndrome (defined by obesity, diabetes, and hypertension) were compared with similar patients without metabolic syndrome (those with less than three criteria). INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: The primary outcome was 28-day mortality. Among 4288 ARDS trial participants, 454 (10.6%) with metabolic syndrome were compared with 3834 controls (89.4%). In adjusted analyses, the metabolic syndrome group was associated with lower 28-day and 90-day mortality when compared with control (adjusted odds ratio [aOR], 0.70 [95% CI, 0.55-0.89] and 0.75 [95% CI, 0.60-0.95], respectively). With each additional metabolic criterion from 0 to 3, adjusted 28-day mortality was reduced by 18%, 22%, and 40%, respectively. In subgroup analyses stratifying by ARDS etiology, mortality was lower for metabolic syndrome vs. control in ARDS caused by sepsis or pneumonia (at 28 d, aOR 0.64 [95% CI, 0.48-0.84] and 90 d, aOR 0.69 [95% CI, 0.53-0.89]), but not in ARDS from noninfectious causes (at 28 d, aOR 1.18 [95% CI, 0.70-1.99] and 90 d, aOR 1.26 [95% CI, 0.77-2.06]). Interaction p = 0.04 and p = 0.02 for 28- and 90-day comparisons, respectively. CONCLUSIONS: Metabolic syndrome in ARDS was associated with a lower risk of mortality in non-COVID-19 ARDS. The relationship between metabolic inflammation and ARDS may provide a novel biological pathway to be explored in precision medicine-based trials.

Computational Analysis Concerning the Impact of DNA Accessibility on CRISPR-Cas9 Cleavage Efficiency.

Defining the variables that impact the specificity of CRISPR/Cas9 has been a major research focus. Whereas sequence complementarity between guide RNA and target DNA substantially dictates cleavage efficiency, DNA accessibility of the targeted loci has also been hypothesized to be an important factor. In this study, functional data from two genome-wide assays, genome-wide, unbiased identification of DSBs enabled by sequencing (GUIDE-seq) and circularization for in vitro reporting of cleavage effects by sequencing (CIRCLE-seq), have been computationally analyzed in conjunction with DNA accessibility determined via DNase I-hypersensitive sequencing from the Encyclopedia of DNA Elements (ENCODE) Database and transcriptome from the Sequence Read Archive to determine whether cellular factors influence CRISPR-induced cleavage efficiency. CIRCLE-seq and GUIDE-seq datasets were selected to represent the absence and presence of cellular factors, respectively. Data analysis revealed that correlations between sequence similarity and CRISPR-induced cleavage frequency were altered by the presence of cellular factors that modulated the level of DNA accessibility. The above-mentioned correlation was abolished when cleavage sites were located in less accessible regions. Furthermore, CRISPR-mediated edits were permissive even at regions that were insufficient for most endogenous genes to be expressed. These results provide a strong basis to dissect the contribution of local chromatin modulation markers on CRISPR-induced cleavage efficiency.

Green and Highly-Efficient Microwave Synthesis Route for Sulfur/Carbon Composite for Li-S Battery.

Multiporous carbons (MPCs) are prepared using ZnO as a hard template and biomass pyrolysis oil as the carbon source. It is shown that the surface area, pore volume, and mesopore/micropore ratio of the as-prepared MPCs can be easily controlled by adjusting the ZnO/oil ratio. Sulfur/MPC (S/MPC) composite is prepared by blending sulfur powder with the as-prepared MPCs followed by microwave heating at three different powers (100 W/200 W/300 W) for 60 s. The unique micro/mesostructure characteristics of the resulting porous carbons not only endow the S/MPC composite with sufficient available space for sulfur storage, but also provide favorable and efficient channels for Li-ions/electrons transportation. When applied as the electrode material in a lithium-ion battery (LIB), the S/MPC composite shows a reversible capacity (about 500 mAh g-1) and a high columbic efficiency (>95%) after 70 cycles. Overall, the method proposed in this study provides a simple and green approach for the rapid production of MPCs and S/MPC composite for high-performance LIBs.

Gate voltage impact on charge mobility in end-on stacked conjugated oligomers.

We present a model of the charge transport in thin film organic field-effect transistors with the active channel made of linear conjugated chains stacked on the substrate with end-on-orientation. The transport was simulated in a box consisting of 25 polymer chains, in which the delocalized quantum orbital eigenstates of the on-chain hole distribution were calculated. The inter-chain charge transfer was solved semi-classically. The full self-consistent distribution of charge density and electric field was determined for various applied gate and source-drain voltages. We found that the dependence of charge mobility on gate voltage is not monotonic: it first increases with increasing gate voltage for a limited interval of the latter, otherwise it decreases with the gate voltage. Next, we found formation of the second resonant peak for higher gate voltages. The mobility dependence on the gate voltage confirmed that the current flowing through the active semiconductor layer should be described not only as the hole transfer between adjacent repeat units of the neighbouring chains, but also as the transfer of coherences among on-chain repeat units. The presented model can also give a new insight into the charge transport in organic field-effect transistors with a novel vertical architecture.

Somaclonal variation does not preclude the use of rice transformants for genetic screening.

Rice (Oryza sativa) is one of the world's most important crops. Rice researchers make extensive use of insertional mutants for the study of gene function. Approximately half a million flanking sequence tags from rice insertional mutant libraries are publicly available. However, the relationship between genotype and phenotype is very weak. Transgenic plant assays have been used frequently for complementation, overexpression or antisense analysis, but sequence changes caused by callus growth, Agrobacterium incubation medium, virulence genes, transformation and selection conditions are unknown. We used high-throughput sequencing of DNA from rice lines derived from Tainung 67 to analyze non-transformed and transgenic rice plants for mutations caused by these parameters. For comparison, we also analyzed sequence changes for two additional rice varieties and four T-DNA tagged transformants from the Taiwan Rice Insertional Mutant resource. We identified single-nucleotide polymorphisms, small indels, large deletions, chromosome doubling and chromosome translocations in these lines. Using standard rice regeneration/transformation procedures, the mutation rates of regenerants and transformants were relatively low, with no significant differences among eight tested treatments in the Tainung 67 background and in the cultivars Taikeng 9 and IR64. Thus, we could not conclusively detect sequence changes resulting from Agrobacterium-mediated transformation in addition to those caused by tissue culture-induced somaclonal variation. However, the mutation frequencies within the two publically available tagged mutant populations, including TRIM transformants or Tos17 lines, were about 10-fold higher than the frequency of standard transformants, probably because mass production of embryogenic calli and longer callus growth periods were required to generate these large libraries.

Predicting genome terminus sequences of Bacillus cereus-group bacteriophage using next generation sequencing data.

BACKGROUND: Most tailed bacteriophages (phages) feature linear dsDNA genomes. Characterizing novel phages requires an understanding of complete genome sequences, including the definition of genome physical ends. RESULT: We sequenced 48 Bacillus cereus phage isolates and analyzed Next-generation sequencing (NGS) data to resolve the genome configuration of these novel phages. Most assembled contigs featured reads that mapped to both contig ends and formed circularized contigs. Independent assemblies of 31 nearly identical I48-like Bacillus phage isolates allowed us to observe that the assembly programs tended to produce random cleavage on circularized contigs. However, currently available assemblers were not capable of reporting the underlying phage genome configuration from sequence data. To identify the genome configuration of sequenced phage in silico, a terminus prediction method was developed by means of 'neighboring coverage ratios' and 'read edge frequencies' from read alignment files. Termini were confirmed by primer walking and supported by phylogenetic inference of large DNA terminase protein sequences. CONCLUSIONS: The Terminus package using phage NGS data along with the contig circularity could efficiently identify the proximal positions of phage genome terminus. Complete phage genome sequences allow a proposed characterization of the potential packaging mechanisms and more precise genome annotation.

Dynamic association rules for gene expression data analysis.

BACKGROUND: The purpose of gene expression analysis is to look for the association between regulation of gene expression levels and phenotypic variations. This association based on gene expression profile has been used to determine whether the induction/repression of genes correspond to phenotypic variations including cell regulations, clinical diagnoses and drug development. Statistical analyses on microarray data have been developed to resolve gene selection issue. However, these methods do not inform us of causality between genes and phenotypes. In this paper, we propose the dynamic association rule algorithm (DAR algorithm) which helps ones to efficiently select a subset of significant genes for subsequent analysis. The DAR algorithm is based on association rules from market basket analysis in marketing. We first propose a statistical way, based on constructing a one-sided confidence interval and hypothesis testing, to determine if an association rule is meaningful. Based on the proposed statistical method, we then developed the DAR algorithm for gene expression data analysis. The method was applied to analyze four microarray datasets and one Next Generation Sequencing (NGS) dataset: the Mice Apo A1 dataset, the whole genome expression dataset of mouse embryonic stem cells, expression profiling of the bone marrow of Leukemia patients, Microarray Quality Control (MAQC) data set and the RNA-seq dataset of a mouse genomic imprinting study. A comparison of the proposed method with the t-test on the expression profiling of the bone marrow of Leukemia patients was conducted. RESULTS: We developed a statistical way, based on the concept of confidence interval, to determine the minimum support and minimum confidence for mining association relationships among items. With the minimum support and minimum confidence, one can find significant rules in one single step. The DAR algorithm was then developed for gene expression data analysis. Four gene expression datasets showed that the proposed DAR algorithm not only was able to identify a set of differentially expressed genes that largely agreed with that of other methods, but also provided an efficient and accurate way to find influential genes of a disease. CONCLUSIONS: In the paper, the well-established association rule mining technique from marketing has been successfully modified to determine the minimum support and minimum confidence based on the concept of confidence interval and hypothesis testing. It can be applied to gene expression data to mine significant association rules between gene regulation and phenotype. The proposed DAR algorithm provides an efficient way to find influential genes that underlie the phenotypic variance.

Multilocus approach to clarify species status and the divergence history of the Bemisia tabaci (Hemiptera: Aleyrodidae) species complex.

The sweet potato whitefly, Bemisia tabaci, is a highly differentiated species complex. Despite consisting of several morphologically indistinguishable entities and frequent invasions on all continents with important associated economic losses, the phylogenetic relationships, species status, and evolutionary history of this species complex is still debated. We sequenced and analyzed one mitochondrial and three single-copy nuclear genes from 9 of the 12 genetic groups of B. tabaci and 5 closely related species. Bayesian species delimitation was applied to investigate the speciation events of B. tabaci. The species statuses of the different genetic groups were strongly supported under different prior settings and phylogenetic scenarios. Divergence histories were estimated by a multispecies coalescence approach implemented in (*)BEAST. Based on mitochondrial locus, B. tabaci was originated 6.47 million years ago (MYA). Nevertheless, the time was 1.25MYA based on nuclear loci. According to the method of approximate Bayesian computation, this difference is probably due to different degrees of migration among loci; i.e., although the mitochondrial locus had differentiated, gene flow at nuclear loci was still possible, a scenario similar to parapatric mode of speciation. This is the first study in whiteflies using multilocus data and incorporating Bayesian coalescence approaches, both of which provide a more biologically realistic framework for delimiting species status and delineating the divergence history of B. tabaci. Our study illustrates that gene flow during species divergence should not be overlooked and has a great impact on divergence time estimation.

ASK1 promotes apoptosis of normal and malignant plasma cells.

Although the overproduction of immunoglobulins by short-lived plasma cells accompanying an immune response links with their apoptosis, how long-lived plasma cells adapt to ensure their longevity in this context is obscure. Here, we show that apoptosis signal-regulating kinase 1 (ASK1) contributes to apoptosis of plasma cells because ASK1 activity was induced during differentiation of short-lived plasma cells, and, when produced by ASK1-deficient mice, these cells survived better than those of control mice. Moreover, antigen-specific long-lived plasma cells generated by immunization accumulated in ASK1-deficient mice, suggesting ASK1 also plays a negative role in survival of long-lived plasma cells. In malignant plasma cells, ASK1 transcription was directly suppressed by B lymphocyte-induced maturation protein-1 (Blimp-1). The expression of ASK1 and Blimp-1 showed an inverse correlation between normal human mature B cells and bone marrow plasma cells from patients with multiple myeloma (MM). Suppression of ASK1 is crucial for cell survival because its enforced expression in MM cells caused apoptosis in vitro and lowered MM load in a xenograft animal model; furthermore, alteration of ASK1 activity affected MM cell survival. Our findings indicate a novel mechanism underlying the regulation of survival in normal and malignant plasma cells by ASK1.

STOPTHEBURN: A Randomized Controlled Trial of Death Cafés for Burnout Prevention in Intensive Care Unit Employees.

Rationale: Effective interventions to prevent burnout among intensive care unit (ICU) clinicians are urgently needed. Death cafés, group discussions about death, build a sense of community and create a space for reflection on distressing events. Objective: To assess whether participation in regular death cafés can prevent burnout in ICU clinicians (physicians, nurses, pharmacists, therapists). Methods: A randomized clinical trial was conducted from July 2020 to December 2022 in 10 ICUs in Louisiana. Subjects were randomized to attend four psychotherapist-facilitated virtual death cafés or to a control arm. The primary outcome was burnout defined by the Maslach Burnout Inventory-Human Services Survey at 6 months. Depression and anxiety scores were measured, as were qualitative data on stressors, coping, and death café experience. Results: Among 340 clinicians who were screened and gave consent (171 physicians, 169 nonphysicians), 251 participated (mean age, 31.0 ± 6.8 years; 63% female; 72% White; 37% nurses, 27% residents, 25% interns, 11% other). Burnout prevalence was 19% at baseline. Of 136 participants who completed the 6-month follow-up, no significant differences were found between intervention and control for the primary outcome (18% vs. 25%; unadjusted odds ratio, 0.64; 95% confidence interval, 0.26-1.57; P = 0.33). There were no differences in anxiety or depression. Notably, the study was limited by an inability to achieve target enrollment and a high attrition rate (46%). Conclusions: Virtual death cafés were unable to reduce burnout, although the study was underpowered to detect differences between groups. Clinical trial registered with clinicaltrials.gov (NCT04347811).

Production of rAAV by plasmid transfection induces antiviral and inflammatory responses in suspension HEK293 cells.

Recombinant adeno-associated virus (rAAV) is a clinically proven viral vector for delivery of therapeutic genes to treat rare diseases. Improving rAAV manufacturing productivity and vector quality is necessary to meet clinical and commercial demand. These goals will require an improved understanding of the cellular response to rAAV production, which is poorly defined. We interrogated the kinetic transcriptional response of HEK293 cells to rAAV production following transient plasmid transfection, under manufacturing-relevant conditions, using RNA-seq. Time-series analyses identified a robust cellular response to transfection and rAAV production, with 1,850 transcripts differentially expressed. Gene Ontology analysis determined upregulated pathways, including inflammatory and antiviral responses, with several interferon-stimulated cytokines and chemokines being upregulated at the protein level. Literature-based pathway prediction implicated multiple pathogen pattern sensors and signal transducers in up-regulation of inflammatory and antiviral responses in response to transfection and rAAV replication. Systematic analysis of the cellular transcriptional response to rAAV production indicates that host cells actively sense vector manufacture as an infectious insult. This dataset may therefore illuminate genes and pathways that influence rAAV production, thereby enabling the rational design of next-generation manufacturing platforms to support safe, effective, and affordable AAV-based gene therapies.

Transcriptional activation of Jun and Fos members of the AP-1 complex is a conserved signature of immune aging that contributes to inflammaging.

Diverse mouse strains have different health and life spans, mimicking the diversity among humans. To capture conserved aging signatures, we studied long-lived C57BL/6J and short-lived NZO/HILtJ mouse strains by profiling transcriptomes and epigenomes of immune cells from peripheral blood and the spleen from young and old mice. Transcriptional activation of the AP-1 transcription factor complex, particularly Fos, Junb, and Jun genes, was the most significant and conserved aging signature across tissues and strains. ATAC-seq data analyses showed that the chromatin around these genes was more accessible with age and there were significantly more binding sites for these TFs with age across all studied tissues, targeting pro-inflammatory molecules including Il6. Age-related increases in binding sites of JUN and FOS factors were also conserved in human peripheral blood ATAC-seq data. Single-cell RNA-seq data from the mouse aging cell atlas Tabula Muris Senis showed that the expression of these genes increased with age in B, T, NK cells, and macrophages, with macrophages from old mice expressing these molecules more abundantly than other cells. Functional data showed that upon myeloid cell activation via poly(I:C), the levels of JUN protein and its binding activity increased more significantly in spleen cells from old compared to young mice. In addition, upon activation, old cells produced more IL6 compared to young cells. In sum, we showed that the aging-related transcriptional activation of Jun and Fos family members in AP-1 complex is conserved across immune tissues and long- and short-living mouse strains, possibly contributing to increased inflammation with age.

Assessment of anti-HIV-1 guide RNA efficacy in cells containing the viral target sequence, corresponding gRNA, and CRISPR/Cas9.

The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 gene editing system has been shown to be effective at inhibiting human immunodeficiency virus type 1 (HIV-1). Studies have not consistently used a trackable dual reporter system to determine what cells received the Cas9/gRNA to determine the overall knockdown of HIV. Some studies have used stably transduced cells under drug selection to accomplish this goal. Here a two-color system was used that allows tracking of viral protein expression and which cells received the CRISPR/Cas9 system. These experiments ensured that each gRNA used was a perfect match to the intended target to remove this variable. The data showed that gRNAs targeting the transactivation response element (TAR) region or other highly conserved regions of the HIV-1 genome were effective at stopping viral gene expression, with multiple assays demonstrating greater than 95 percent reduction. Conversely, gRNAs targeting conserved sites of the 5' portion of the U3 region were largely ineffective, demonstrating that the location of edits in the long terminal repeat (LTR) matter with respect to function. In addition, it was observed that a gRNA targeting Tat was effective in a T-cell model of HIV-1 latency. Taken together, these studies demonstrated gRNAs designed to highly conserved functional regions have near 100% efficacy in vitro in cells known to have received the Cas9/gRNA pair.

Computational Design of gRNAs Targeting Genetic Variants Across HIV-1 Subtypes for CRISPR-Mediated Antiviral Therapy.

Clustered regularly interspaced short palindromic repeats (CRISPR)-based HIV-1 genome editing has shown promising outcomes in in vitro and in vivo viral infection models. However, existing HIV-1 sequence variants have been shown to reduce CRISPR-mediated efficiency and induce viral escape. Two metrics, global patient coverage and global subtype coverage, were used to identify guide RNA (gRNA) sequences that account for this viral diversity from the perspectives of cross-patient and cross-subtype gRNA design, respectively. Computational evaluation using these parameters and over 3.6 million possible 20-bp sequences resulted in nine lead gRNAs, two of which were previously published. This analysis revealed the benefit and necessity of considering all sequence variants for gRNA design. Of the other seven identified novel gRNAs, two were of note as they targeted interesting functional regions. One was a gRNA predicted to induce structural disruption in the nucleocapsid binding site (Ψ), which holds the potential to stop HIV-1 replication during the viral genome packaging process. The other was a reverse transcriptase (RT)-targeting gRNA that was predicted to cleave the subdomain responsible for dNTP incorporation. CRISPR-mediated sequence edits were predicted to occur on critical residues where HIV-1 has been shown to develop resistance against antiretroviral therapy (ART), which may provide additional evolutionary pressure at the DNA level. Given these observations, consideration of broad-spectrum gRNAs and cross-subtype diversity for gRNA design is not only required for the development of generalizable CRISPR-based HIV-1 therapy, but also helps identify optimal target sites.

Off-Target Analysis in Gene Editing and Applications for Clinical Translation of CRISPR/Cas9 in HIV-1 Therapy.

As genome-editing nucleases move toward broader clinical applications, the need to define the limits of their specificity and efficiency increases. A variety of approaches for nuclease cleavage detection have been developed, allowing a full-genome survey of the targeting landscape and the detection of a variety of repair outcomes for nuclease-induced double-strand breaks. Each approach has advantages and disadvantages relating to the means of target-site capture, target enrichment mechanism, cellular environment, false discovery, and validation of bona fide off-target cleavage sites in cells. This review examines the strengths, limitations, and origins of the different classes of off-target cleavage detection systems including anchored primer enrichment (GUIDE-seq), in situ detection (BLISS), in vitro selection libraries (CIRCLE-seq), chromatin immunoprecipitation (ChIP) (DISCOVER-Seq), translocation sequencing (LAM PCR HTGTS), and in vitro genomic DNA digestion (Digenome-seq and SITE-Seq). Emphasis is placed on the specific modifications that give rise to the enhanced performance of contemporary techniques over their predecessors and the comparative performance of techniques for different applications. The clinical relevance of these techniques is discussed in the context of assessing the safety of novel CRISPR/Cas9 HIV-1 curative strategies. With the recent success of HIV-1 and SIV-1 viral suppression in humanized mice and non-human primates, respectively, using CRISPR/Cas9, rigorous exploration of potential off-target effects is of critical importance. Such analyses would benefit from the application of the techniques discussed in this review.

Designing Safer CRISPR/Cas9 Therapeutics for HIV: Defining Factors That Regulate and Technologies Used to Detect Off-Target Editing.

Human immunodeficiency virus type-1 (HIV-1) infection has resulted in the death of upward of 39 million people since being discovered in the early 1980s. A cure strategy for HIV-1 has eluded scientists, but gene editing technologies such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) offer a new approach to developing a cure for HIV infection. While the CRISPR/Cas9 system has been used successfully in a number of different types of studies, there remains a concern for off-target effects. This review details the different aspects of the Cas9 system and how they play a role in off-target events. In addition, this review describes the current technologies available for detecting off-target cleavage events and their advantages and disadvantages. While some studies have utilized whole genome sequencing (WGS), this method sacrifices depth of coverage for interrogating the whole genome. A number of different approaches have now been developed to take advantage of next generation sequencing (NGS) without sacrificing depth of coverage. This review highlights four widely used methods for detecting off-target events: (1) genome-wide unbiased identification of double-stranded break events enabled by sequencing (GUIDE-Seq), (2) discovery of in situ Cas off-targets and verification by sequencing (DISCOVER-Seq), (3) circularization for in vitro reporting of cleavage effects by sequencing (CIRCLE-Seq), and (4) breaks labeling in situ and sequencing (BLISS). Each of these technologies has advantages and disadvantages, but all center around capturing double-stranded break (DSB) events catalyzed by the Cas9 endonuclease. Being able to define off-target events is crucial for a gene therapy cure strategy for HIV-1.

Integrated Human Immunodeficiency Virus Type 1 Sequence in J-Lat 10.6.

The full length of HIV/R7/E-/GFP integrated in the J-Lat 10.6 cell line was sequenced in this study. The single copy of the integrated virus, including the breakpoints from the human chromosome to the provirus, was amplified by two separate PCRs. A 10,200-bp genome sequence was acquired, analyzed, and deposited in GenBank.

Safe CRISPR-Cas9 Inhibition of HIV-1 with High Specificity and Broad-Spectrum Activity by Targeting LTR NF-κB Binding Sites.

Viral latency of human immunodeficiency virus type 1 (HIV-1) has become a major hurdle to a cure in the highly effective antiretroviral therapy (ART) era. The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has successfully been demonstrated to excise or inactivate integrated HIV-1 provirus from infected cells by targeting the long terminal repeat (LTR) region. However, the guide RNAs (gRNAs) have classically avoided transcription factor binding sites (TFBSs) that are readily observed and known to be important in human promoters. Although conventionally thought unfavorable due to potential impact on human promoters, our computational pipeline identified gRNA sequences that were predicted to inactivate HIV-1 transcription by targeting the nuclear factor κB (NF-κB) binding sites (gNFKB0, gNFKB1) with a high safety profile (lack of predicted or observed human edits) and broad-spectrum activity (predicted coverage of known viral sequences). Genome-wide, unbiased identification of double strand breaks (DSBs) enabled by sequencing (GUIDE-seq) showed that the gRNAs targeting NF-κB binding sites had no detectable CRISPR-induced off-target edits in HeLa cells. 5' LTR-driven HIV-1 transcription was significantly reduced in three HIV-1 reporter cell lines. These results demonstrate a working model to specifically target well-known TFBSs in the HIV-1 LTR that are readily observed in human promoters to reduce HIV-1 transcription with a high-level safety profile and broad-spectrum activity.