Identifying genetic variants that influence the abundance of cell states in single-cell data.

Disease risk alleles influence the composition of cells present in the body, but modeling genetic effects on the cell states revealed by single-cell profiling is difficult because variant-associated states may reflect diverse combinations of the profiled cell features that are challenging to predefine. We introduce Genotype-Neighborhood Associations (GeNA), a statistical tool to identify cell-state abundance quantitative trait loci (csaQTLs) in high-dimensional single-cell datasets. Instead of testing associations to predefined cell states, GeNA flexibly identifies the cell states whose abundance is most associated with genetic variants. In a genome-wide survey of single-cell RNA sequencing peripheral blood profiling from 969 individuals, GeNA identifies five independent loci associated with shifts in the relative abundance of immune cell states. For example, rs3003-T (P = 1.96 × 10-11) associates with increased abundance of natural killer cells expressing tumor necrosis factor response programs. This csaQTL colocalizes with increased risk for psoriasis, an autoimmune disease that responds to anti-tumor necrosis factor treatments. Flexibly characterizing csaQTLs for granular cell states may help illuminate how genetic background alters cellular composition to confer disease risk.

The influence of HLA genetic variation on plasma protein expression.

Genetic variation in the human leukocyte antigen (HLA) loci is associated with risk of immune-mediated diseases, but the molecular effects of HLA polymorphism are unclear. Here we examined the effects of HLA genetic variation on the expression of 2940 plasma proteins across 45,330 Europeans in the UK Biobank, with replication analyses across multiple ancestry groups. We detected 504 proteins affected by HLA variants (HLA-pQTL), including widespread trans effects by autoimmune disease risk alleles. More than 80% of the HLA-pQTL fine-mapped to amino acid positions in the peptide binding groove. HLA-I and II affected proteins expressed in similar cell types but in different pathways of both adaptive and innate immunity. Finally, we investigated potential HLA-pQTL effects on disease by integrating HLA-pQTL with fine-mapped HLA-disease signals in the UK Biobank. Our data reveal the diverse effects of HLA genetic variation and aid the interpretation of associations between HLA alleles and immune-mediated diseases.

Use of Diagnostic Codes for Primary Open-Angle Glaucoma Polygenic Risk Score Construction in Electronic Health Record-Linked Biobanks.

PURPOSE: Polygenic risk scores (PRSs) likely predict risk and prognosis of glaucoma. We compared the PRS performance for primary open-angle glaucoma (POAG), defined using International Classification of Diseases (ICD) codes vs manual medical record review. DESIGN: Retrospective cohort study. METHODS: We identified POAG cases in the Mount Sinai BioMe and Mass General Brigham (MGB) biobanks using ICD codes. We confirmed POAG based on optical coherence tomograms and visual fields. In a separate 5% sample, the absence of POAG was confirmed with intraocular pressure and cup-disc ratio criteria. We used genotype data and either self-reported glaucoma diagnoses or ICD-10 codes for glaucoma diagnoses from the UK Biobank and the lassosum method to compute a genome-wide POAG PRS. We compared the area under the curve (AUC) for POAG prediction based on ICD codes vs medical records. RESULTS: We reviewed 804 of 996 BioMe and 367 of 1006 MGB ICD-identified cases. In BioMe and MGB, respectively, positive predictive value was 53% and 55%; negative predictive value was 96% and 97%; sensitivity was 97% and 97%; and specificity was 44% and 53%. Adjusted PRS AUCs for POAG using ICD codes vs manual record review in BioMe were not statistically different (P ≥.21) by ancestry: 0.77 vs 0.75 for African, 0.80 vs 0.80 for Hispanic, and 0.81 vs 0.81 for European. Results were similar in MGB (P ≥.18): 0.72 vs 0.80 for African, 0.83 vs 0.86 for Hispanic, and 0.74 vs 0.73 for European. CONCLUSIONS: A POAG PRS performed similarly using either manual review or ICD codes in 2 electronic health record-linked biobanks; manual assessment of glaucoma status might not be necessary for some PRS studies. However, caution should be exercised when using ICD codes for glaucoma diagnosis given their low specificity (44%-53%) for manually confirmed cases of glaucoma.

Acute Steroid-Refractory Gastrointestinal Graft-Versus-Host Disease Is Not Associated With Significant Differences in Gut Taxonomic Composition Compared to Steroid-Sensitive Gastrointestinal Graft-Versus-Host Disease Immediately Before Onset of Disease.

Taxonomic composition of the gut microbiota at the time of neutrophil engraftment is associated with the development of acute gastrointestinal graft-versus-host disease (GI GVHD) in patients undergoing allogeneic hematopoietic stem cell transplantation. However, less is known about the relationship between the gut microbiota and development of steroid-refractory GI GVHD immediately before the onset of disease. Markers of steroid-refractory GI GVHD are needed to identify patients who may benefit from the early initiation of non-corticosteroid-based GVHD treatment. Our aim was to identify differences in taxonomic composition in stool samples from patients without GVHD, with steroid-responsive GVHD and with steroid-refractory GI GVHD to identify predictive microbiome biomarkers of steroid-refractory GI GVHD. We conducted a retrospective case-control, single institution study, performing shotgun metagenomic sequencing on stool samples from patients with (n = 36) and without GVHD (n = 34) matched for time since transplantation. We compared the taxonomic composition of the gut microbiome in those with steroid-sensitive GI GVHD (n = 17) and steroid-refractory GI GVHD (n = 19) to each other and to those without GVHD. We also performed associations between steroid-refractory GI GVHD, gut taxonomic composition, and fecal calprotectin, a marker of GI GVHD to develop composite fecal markers of steroid-refractory GVHD before the onset of GI disease. We found that fecal samples within 30 days of GVHD onset from patients with and without GVHD or with and without steroid-refractory GI GVHD did not differ significantly in Shannon diversity (alpha-diversity) or in overall taxonomic composition (beta-diversity). Although those patients without GVHD had higher relative abundance of Clostridium spp., those with and without steroid-refractory GI GVHD did not significantly differ in taxonomic composition between one another. In our study, fecal calprotectin before disease onset was significantly higher in patients with GVHD compared to those without GVHD and higher in patients with steroid-refractory GI GVHD compared to steroid-sensitive GI GVHD. No taxa were significantly associated with higher levels of calprotectin.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Deconstruction of rheumatoid arthritis synovium defines inflammatory subtypes.

Rheumatoid arthritis is a prototypical autoimmune disease that causes joint inflammation and destruction1. There is currently no cure for rheumatoid arthritis, and the effectiveness of treatments varies across patients, suggesting an undefined pathogenic diversity1,2. Here, to deconstruct the cell states and pathways that characterize this pathogenic heterogeneity, we profiled the full spectrum of cells in inflamed synovium from patients with rheumatoid arthritis. We used multi-modal single-cell RNA-sequencing and surface protein data coupled with histology of synovial tissue from 79 donors to build single-cell atlas of rheumatoid arthritis synovial tissue that includes more than 314,000 cells. We stratified tissues into six groups, referred to as cell-type abundance phenotypes (CTAPs), each characterized by selectively enriched cell states. These CTAPs demonstrate the diversity of synovial inflammation in rheumatoid arthritis, ranging from samples enriched for T and B cells to those largely lacking lymphocytes. Disease-relevant cell states, cytokines, risk genes, histology and serology metrics are associated with particular CTAPs. CTAPs are dynamic and can predict treatment response, highlighting the clinical utility of classifying rheumatoid arthritis synovial phenotypes. This comprehensive atlas and molecular, tissue-based stratification of rheumatoid arthritis synovial tissue reveal new insights into rheumatoid arthritis pathology and heterogeneity that could inform novel targeted treatments.

Participant Experience Using Novel Perimetry Tests to Monitor Glaucoma Progression.

PRCIS: Participant surveys taken after using tablet-based and smart visual function analyzer (SVFA) perimetry tests suggest that patients may prefer novel perimetry tests over traditional visual field machines. PURPOSE: Compare patient experience using the IMOvifa SVFA and the tablet-based Melbourne Rapid Fields visual field (VF) tests to the Humphrey Field 24-2 Swedish Interactive Threshold Algorithm Standard. PATIENTS AND METHODS: Prospective observational cohort study on adult participants with diagnoses of glaucoma suspect, ocular hypertension, or glaucoma. Participants attended 2 study visits ~3 months apart. During the first visit, participants were trained to use the 2 novel perimeters, took 1 test on both new devices and the Humphrey Field Analyzer, then were surveyed. Participants received tablets and performed weekly tablet VF tests at home between study visits. At the final study visit, participants re-took the VF tests and completed the same surveys. RESULTS: Eighty-one participants were surveyed twice. At the baseline survey, participants preferred the SVFA (71.7%) and tablet tests (69.2%) over the Humphrey VF. Most were willing to perform weekly monitoring at home on the SVFA (69.1%) and tablet (75.4%). Participants generally had a "very good" overall experience when testing on the SVFA (71.6%) and tablet (90.1%). At the final visit, fewer participants were willing to test on the tablet daily (23.5% to 9.9%; P = 0.02 for change) and more were willing to test monthly (18.5% to 33.3%; P = 0.03 for change). CONCLUSION: Users reported a preference for novel VF devices. Overall participant experience using these devices was positive, supporting the feasibility of home monitoring of VFs from an experience perspective.

Tutorial: a statistical genetics guide to identifying HLA alleles driving complex disease.

The human leukocyte antigen (HLA) locus is associated with more complex diseases than any other locus in the human genome. In many diseases, HLA explains more heritability than all other known loci combined. In silico HLA imputation methods enable rapid and accurate estimation of HLA alleles in the millions of individuals that are already genotyped on microarrays. HLA imputation has been used to define causal variation in autoimmune diseases, such as type I diabetes, and in human immunodeficiency virus infection control. However, there are few guidelines on performing HLA imputation, association testing, and fine mapping. Here, we present a comprehensive tutorial to impute HLA alleles from genotype data. We provide detailed guidance on performing standard quality control measures for input genotyping data and describe options to impute HLA alleles and amino acids either locally or using the web-based Michigan Imputation Server, which hosts a multi-ancestry HLA imputation reference panel. We also offer best practice recommendations to conduct association tests to define the alleles, amino acids, and haplotypes that affect human traits. Along with the pipeline, we provide a step-by-step online guide with scripts and available software ( https://github.com/immunogenomics/HLA_analyses_tutorial ). This tutorial will be broadly applicable to large-scale genotyping data and will contribute to defining the role of HLA in human diseases across global populations.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci.

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Comparison of Perimetric Outcomes from a Tablet Perimeter, Smart Visual Function Analyzer, and Humphrey Field Analyzer.

PURPOSE: The tablet-based Melbourne Rapid Fields (MRF) visual field (VF) test and the IMOvifa Smart Visual Function Analyzer (SVFA) are portable perimeters that may allow for at-home monitoring and more frequent testing. We compared tablet and SVFA results with outputs from the Humphrey Field Analyzer (HFA) 24-2 Swedish Interactive Threshold Algorithm Standard program. DESIGN: Observational cross-sectional study. SUBJECTS: Adult participants with a diagnosis of glaucoma, suspected glaucoma, or ocular hypertension seen in the Massachusetts Eye and Ear glaucoma clinic were enrolled. All participants were reliable and experienced HFA testers. METHODS: Participants were tested with the SVFA and HFA. The study staff also trained participants on the MRF tablet with instructions to take weekly tests at home for 3 months. Visual field results from the 3 devices were compared. MAIN OUTCOME MEASURES: Mean deviation (MD), pattern standard deviation (PSD), reliability parameters, and point sensitivity. RESULTS: Overall, 79 participants (133 eyes) with a mean age of 61 ± 13 years (range, 26-79 years) were included; 59% of the participants were female, and the mean HFA MD was -2.7 ± 3.9 dB. The global indices of MD and PSD did not significantly vary between HFA and the 2 novel devices, except that the tablet VF reported a 0.6 dB higher PSD compared with HFA. However, tablet and SVFA sensitivities significantly differed from those of the HFA at 36 and 39 locations, respectively, out of 52 locations. Relative to HFA, the tablet overestimated light sensitivity in the nasal field while underestimating the temporal field. The SVFA generally underestimated light sensitivity, but its results were more similar to HFA results compared with the tablet. CONCLUSIONS: Although average MD values from the 2 novel devices suggest that they provide similar results to the HFA, point-by-point comparisons highlight notable deviations. Differences in specific point sensitivity values were significant, especially between the tablet and the other 2 devices. These differences may in part be explained by differences in the devices' normative databases as well as how MD is calculated. However, the tablet had substantial differences based on location, indicating that the tablet design itself may be responsible for differences in local sensitivities. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Relationship between Claims-Based Frailty Index and Eye Care Utilization among Medicare Beneficiaries with Glaucoma.

PURPOSE: To determine differences in eye care utilization by frailty levels among Medicare beneficiaries with glaucoma. DESIGN: Retrospective cohort study. PARTICIPANTS: Medicare fee-for-service beneficiaries over 65 years of age with glaucoma, identified using International Classification of Diseases codes before July 1, 2014. METHODS: By using a validated claims-based frailty index (range, 0-1), beneficiaries were classified as nonfrail/prefrail (0-0.19), mildly frail (0.20-0.29), and moderate-to-severely frail (≥ 0.30). Negative binomial regression analyses were used to estimate incident rate ratios (IRRs) of eye care utilization by frailty levels between July 1, 2014, and December 31, 2016. MAIN OUTCOME MEASURES: Current Procedural Terminology codes for eye examinations and eye care-related office visits; eye care-related inpatient and emergency department (ED) encounters; eye care-related nursing facility and home-visit encounters; visual field (VF) and retinal nerve fiber layer (RNFL) OCT tests; and selective laser trabeculoplasties (SLTs) and glaucoma surgeries. RESULTS: Among 76 260 Medicare beneficiaries with glaucoma, the mean age was 78.9 years (standard deviation, 7.8), female beneficiaries constituted 60.5%, and 78.7% of beneficiaries self-identified as non-Hispanic White. According to a claims-based frailty index, 79.5% of beneficiaries were nonfrail/prefrail, 17.1% were mildly frail, and 3.4% were moderate-to-severely frail. Moderate-to-severely frail beneficiaries were less likely than nonfrail/prefrail beneficiaries to have outpatient encounters (IRR, 0.85, 95% confidence interval [CI], 0.83-0.88); VF tests (IRR, 0.64, 95% CI, 0.60-0.67); RNFL OCT tests (IRR, 0.77, 95% CI, 0.73-0.81); SLT (IRR, 0.74, 95% CI, 0.60-0.92); and glaucoma surgery (IRR, 0.74, 95% CI 0.55-0.99), after adjusting for age, gender, glaucoma severity, race, and socioeconomic status. Compared with nonfrail/prefrail beneficiaries, moderate-to-severely frail beneficiaries had higher rates of inpatient/ED encounters (IRR, 5.03, 95% CI, 2.36-10.71) and nursing facility/home-visit encounters (IRR, 34.89, 95% CI, 14.82-82.13). CONCLUSIONS: Compared with nonfrail/prefrail Medicare beneficiaries with glaucoma, beneficiaries with moderate-to-severe frailty had lower rates of eye care utilization in the outpatient setting and higher rates of utilization in acute care settings. This suggests that frail patients may receive less disease monitoring and fewer interventions for their glaucoma management. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

Single-cell multiome of the human retina and deep learning nominate causal variants in complex eye diseases.

Genome-wide association studies (GWASs) of eye disorders have identified hundreds of genetic variants associated with ocular disease. However, the vast majority of these variants are noncoding, making it challenging to interpret their function. Here we present a joint single-cell atlas of gene expression and chromatin accessibility of the adult human retina with more than 50,000 cells, which we used to analyze single-nucleotide polymorphisms (SNPs) implicated by GWASs of age-related macular degeneration, glaucoma, diabetic retinopathy, myopia, and type 2 macular telangiectasia. We integrate this atlas with a HiChIP enhancer connectome, expression quantitative trait loci (eQTL) data, and base-resolution deep learning models to predict noncoding SNPs with causal roles in eye disease, assess SNP impact on transcription factor binding, and define their known and novel target genes. Our efforts nominate pathogenic SNP-target gene interactions for multiple vision disorders and provide a potentially powerful resource for interpreting noncoding variation in the eye.

Cross-tissue, single-cell stromal atlas identifies shared pathological fibroblast phenotypes in four chronic inflammatory diseases.

BACKGROUND: Pro-inflammatory fibroblasts are critical for pathogenesis in rheumatoid arthritis, inflammatory bowel disease, interstitial lung disease, and Sjögren's syndrome and represent a novel therapeutic target for chronic inflammatory disease. However, the heterogeneity of fibroblast phenotypes, exacerbated by the lack of a common cross-tissue taxonomy, has limited our understanding of which pathways are shared by multiple diseases. METHODS: We profiled fibroblasts derived from inflamed and non-inflamed synovium, intestine, lungs, and salivary glands from affected individuals with single-cell RNA sequencing. We integrated all fibroblasts into a multi-tissue atlas to characterize shared and tissue-specific phenotypes. FINDINGS: Two shared clusters, CXCL10+CCL19+ immune-interacting and SPARC+COL3A1+ vascular-interacting fibroblasts, were expanded in all inflamed tissues and mapped to dermal analogs in a public atopic dermatitis atlas. We confirmed these human pro-inflammatory fibroblasts in animal models of lung, joint, and intestinal inflammation. CONCLUSIONS: This work represents a thorough investigation into fibroblasts across organ systems, individual donors, and disease states that reveals shared pathogenic activation states across four chronic inflammatory diseases. FUNDING: Grant from F. Hoffmann-La Roche (Roche) AG.

Repertoire analyses reveal T cell antigen receptor sequence features that influence T cell fate.

T cells acquire a regulatory phenotype when their T cell antigen receptors (TCRs) experience an intermediate- to high-affinity interaction with a self-peptide presented via the major histocompatibility complex (MHC). Using TCRß sequences from flow-sorted human cells, we identified TCR features that promote regulatory T cell (Treg) fate. From these results, we developed a scoring system to quantify TCR-intrinsic regulatory potential (TiRP). When applied to the tumor microenvironment, TiRP scoring helped to explain why only some T cell clones maintained the conventional T cell (Tconv) phenotype through expansion. To elucidate drivers of these predictive TCR features, we then examined the two elements of the Treg TCR ligand separately: the self-peptide and the human MHC class II molecule. These analyses revealed that hydrophobicity in the third complementarity-determining region (CDR3ß) of the TCR promotes reactivity to self-peptides, while TCR variable gene (TRBV gene) usage shapes the TCR's general propensity for human MHC class II-restricted activation.

Co-varying neighborhood analysis identifies cell populations associated with phenotypes of interest from single-cell transcriptomics.

As single-cell datasets grow in sample size, there is a critical need to characterize cell states that vary across samples and associate with sample attributes, such as clinical phenotypes. Current statistical approaches typically map cells to clusters and then assess differences in cluster abundance. Here we present co-varying neighborhood analysis (CNA), an unbiased method to identify associated cell populations with greater flexibility than cluster-based approaches. CNA characterizes dominant axes of variation across samples by identifying groups of small regions in transcriptional space-termed neighborhoods-that co-vary in abundance across samples, suggesting shared function or regulation. CNA performs statistical testing for associations between any sample-level attribute and the abundances of these co-varying neighborhood groups. Simulations show that CNA enables more sensitive and accurate identification of disease-associated cell states than a cluster-based approach. When applied to published datasets, CNA captures a Notch activation signature in rheumatoid arthritis, identifies monocyte populations expanded in sepsis and identifies a novel T cell population associated with progression to active tuberculosis.

Efficient and precise single-cell reference atlas mapping with Symphony.

Recent advances in single-cell technologies and integration algorithms make it possible to construct comprehensive reference atlases encompassing many donors, studies, disease states, and sequencing platforms. Much like mapping sequencing reads to a reference genome, it is essential to be able to map query cells onto complex, multimillion-cell reference atlases to rapidly identify relevant cell states and phenotypes. We present Symphony ( https://github.com/immunogenomics/symphony ), an algorithm for building large-scale, integrated reference atlases in a convenient, portable format that enables efficient query mapping within seconds. Symphony localizes query cells within a stable low-dimensional reference embedding, facilitating reproducible downstream transfer of reference-defined annotations to the query. We demonstrate the power of Symphony in multiple real-world datasets, including (1) mapping a multi-donor, multi-species query to predict pancreatic cell types, (2) localizing query cells along a developmental trajectory of fetal liver hematopoiesis, and (3) inferring surface protein expression with a multimodal CITE-seq atlas of memory T cells.

Maximizing statistical power to detect differentially abundant cell states with scPOST.

To estimate a study design's power to detect differential abundance, we require a framework that simulates many multi-sample single-cell datasets. However, current simulation methods are challenging for large-scale power analyses because they are computationally resource intensive and do not support easy simulation of multi-sample datasets. Current methods also lack modeling of important inter-sample variation, such as the variation in the frequency of cell states between samples that is observed in single-cell data. Thus, we developed single-cell POwer Simulation Tool (scPOST) to address these limitations and help investigators quickly simulate multi-sample single-cell datasets. Users may explore a range of effect sizes and study design choices (such as increasing the number of samples or cells per sample) to determine their effect on power, and thus choose the optimal study design for their planned experiments.

Gliotoxin, identified from a screen of fungal metabolites, disrupts 7SK snRNP, releases P-TEFb, and reverses HIV-1 latency.

A leading pharmacological strategy toward HIV cure requires "shock" or activation of HIV gene expression in latently infected cells with latency reversal agents (LRAs) followed by their subsequent clearance. In a screen for novel LRAs, we used fungal secondary metabolites as a source of bioactive molecules. Using orthogonal mass spectrometry (MS) coupled to latency reversal bioassays, we identified gliotoxin (GTX) as a novel LRA. GTX significantly induced HIV-1 gene expression in latent ex vivo infected primary cells and in CD4+ T cells from all aviremic HIV-1+ participants. RNA sequencing identified 7SK RNA, the scaffold of the positive transcription elongation factor b (P-TEFb) inhibitory 7SK small nuclear ribonucleoprotein (snRNP) complex, to be significantly reduced upon GTX treatment of CD4+ T cells. GTX directly disrupted 7SK snRNP by targeting La-related protein 7 (LARP7), releasing active P-TEFb, which phosphorylated RNA polymerase II (Pol II) C-terminal domain (CTD), inducing HIV transcription.