Multi-omics blood atlas reveals unique features of immune and platelet responses to SARS-CoV-2 Omicron breakthrough infection.

Although host responses to the ancestral SARS-CoV-2 strain are well described, those to the new Omicron variants are less resolved. We profiled the clinical phenomes, transcriptomes, proteomes, metabolomes, and immune repertoires of >1,000 blood cell or plasma specimens from SARS-CoV-2 Omicron patients. Using in-depth integrated multi-omics, we dissected the host response dynamics during multiple disease phases to reveal the molecular and cellular landscapes in the blood. Specifically, we detected enhanced interferon-mediated antiviral signatures of platelets in Omicron-infected patients, and platelets preferentially formed widespread aggregates with leukocytes to modulate immune cell functions. In addition, patients who were re-tested positive for viral RNA showed marked reductions in B cell receptor clones, antibody generation, and neutralizing capacity against Omicron. Finally, we developed a machine learning model that accurately predicted the probability of re-positivity in Omicron patients. Our study may inspire a paradigm shift in studying systemic diseases and emerging public health concerns.

SARS-CoV-2 RNA stabilizes host mRNAs to elicit immunopathogenesis.

SARS-CoV-2 RNA interacts with host factors to suppress interferon responses and simultaneously induces cytokine release to drive the development of severe coronavirus disease 2019 (COVID-19). However, how SARS-CoV-2 hijacks host RNAs to elicit such imbalanced immune responses remains elusive. Here, we analyzed SARS-CoV-2 RNA in situ structures and interactions in infected cells and patient lung samples using RIC-seq. We discovered that SARS-CoV-2 RNA forms 2,095 potential duplexes with the 3' UTRs of 205 host mRNAs to increase their stability by recruiting RNA-binding protein YBX3 in A549 cells. Disrupting the SARS-CoV-2-to-host RNA duplex or knocking down YBX3 decreased host mRNA stability and reduced viral replication. Among SARS-CoV-2-stabilized host targets, NFKBIZ was crucial for promoting cytokine production and reducing interferon responses, probably contributing to cytokine storm induction. Our study uncovers the crucial roles of RNA-RNA interactions in the immunopathogenesis of RNA viruses such as SARS-CoV-2 and provides valuable host targets for drug development.

Structural insight into H4K20 methylation on H2A.Z-nucleosome by SUV420H1.

DNA replication ensures the accurate transmission of genetic information during the cell cycle. Histone variant H2A.Z is crucial for early replication origins licensing and activation in which SUV420H1 preferentially recognizes H2A.Z-nucleosome and deposits H4 lysine 20 dimethylation (H4K20me2) on replication origins. Here, we report the cryo-EM structures of SUV420H1 bound to H2A.Z-nucleosome or H2A-nucleosome and demonstrate that SUV420H1 directly interacts with H4 N-terminal tail, the DNA, and the acidic patch in the nucleosome. The H4 (1-24) forms a lasso-shaped structure that stabilizes the SUV420H1-nucleosome complex and precisely projects the H4K20 residue into the SUV420H1 catalytic center. In vitro and in vivo analyses reveal a crucial role of the SUV420H1 KR loop (residues 214-223), which lies close to the H2A.Z-specific residues D97/S98, in H2A.Z-nucleosome preferential recognition. Together, our findings elucidate how SUV420H1 recognizes nucleosomes to ensure site-specific H4K20me2 modification and provide insights into how SUV420H1 preferentially recognizes H2A.Z nucleosome.

The Transcriptome and DNA Methylome Landscapes of Human Primordial Germ Cells.

Germ cells are vital for transmitting genetic information from one generation to the next and for maintaining the continuation of species. Here, we analyze the transcriptome of human primordial germ cells (PGCs) from the migrating stage to the gonadal stage at single-cell and single-base resolutions. Human PGCs show unique transcription patterns involving the simultaneous expression of both pluripotency genes and germline-specific genes, with a subset of them displaying developmental-stage-specific features. Furthermore, we analyze the DNA methylome of human PGCs and find global demethylation of their genomes. Approximately 10 to 11 weeks after gestation, the PGCs are nearly devoid of any DNA methylation, with only 7.8% and 6.0% of the median methylation levels in male and female PGCs, respectively. Our work paves the way toward deciphering the complex epigenetic reprogramming of the germline with the aim of restoring totipotency in fertilized oocytes.

Complementary Alu sequences mediate enhancer-promoter selectivity.

Enhancers determine spatiotemporal gene expression programs by engaging with long-range promoters1-4. However, it remains unknown how enhancers find their cognate promoters. We recently developed a RNA in situ conformation sequencing technology to identify enhancer-promoter connectivity using pairwise interacting enhancer RNAs and promoter-derived noncoding RNAs5,6. Here we apply this technology to generate high-confidence enhancer-promoter RNA interaction maps in six additional cell lines. Using these maps, we discover that 37.9% of the enhancer-promoter RNA interaction sites are overlapped with Alu sequences. These pairwise interacting Alu and non-Alu RNA sequences tend to be complementary and potentially form duplexes. Knockout of Alu elements compromises enhancer-promoter looping, whereas Alu insertion or CRISPR-dCasRx-mediated Alu tethering to unregulated promoter RNAs can create new loops to homologous enhancers. Mapping 535,404 noncoding risk variants back to the enhancer-promoter RNA interaction maps enabled us to construct variant-to-function maps for interpreting their molecular functions, including 15,318 deletions or insertions in 11,677 Alu elements that affect 6,497 protein-coding genes. We further demonstrate that polymorphic Alu insertion at the PTK2 enhancer can promote tumorigenesis. Our study uncovers a principle for determining enhancer-promoter pairing specificity and provides a framework to link noncoding risk variants to their molecular functions.

Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm.

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.

Interpretation of 10 years of Alzheimer's disease genetic findings in the perspective of statistical heterogeneity.

Common genetic variants and susceptibility loci associated with Alzheimer's disease (AD) have been discovered through large-scale genome-wide association studies (GWAS), GWAS by proxy (GWAX) and meta-analysis of GWAS and GWAX (GWAS+GWAX). However, due to the very low repeatability of AD susceptibility loci and the low heritability of AD, these AD genetic findings have been questioned. We summarize AD genetic findings from the past 10 years and provide a new interpretation of these findings in the context of statistical heterogeneity. We discovered that only 17% of AD risk loci demonstrated reproducibility with a genome-wide significance of P < 5.00E-08 across all AD GWAS and GWAS+GWAX datasets. We highlighted that the AD GWAS+GWAX with the largest sample size failed to identify the most significant signals, the maximum number of genome-wide significant genetic variants or maximum heritability. Additionally, we identified widespread statistical heterogeneity in AD GWAS+GWAX datasets, but not in AD GWAS datasets. We consider that statistical heterogeneity may have attenuated the statistical power in AD GWAS+GWAX and may contribute to explaining the low repeatability (17%) of genome-wide significant AD susceptibility loci and the decreased AD heritability (40-2%) as the sample size increased. Importantly, evidence supports the idea that a decrease in statistical heterogeneity facilitates the identification of genome-wide significant genetic loci and contributes to an increase in AD heritability. Collectively, current AD GWAX and GWAS+GWAX findings should be meticulously assessed and warrant additional investigation, and AD GWAS+GWAX should employ multiple meta-analysis methods, such as random-effects inverse variance-weighted meta-analysis, which is designed specifically for statistical heterogeneity.

The ATF4-RPS19BP1 axis modulates ribosome biogenesis to promote erythropoiesis.

ABSTRACT: Hematopoietic differentiation is controlled by intrinsic regulators and the extrinsic hematopoietic niche. Activating transcription factor 4 (ATF4) plays a crucial role in the function of fetal and adult hematopoietic stem cell maintenance. However, the precise function of ATF4 in the bone marrow (BM) niche and the mechanism by which ATF4 regulates adult hematopoiesis remain largely unknown. Here, we used 4 cell-type-specific mouse Cre lines to achieve conditional knockout of Atf4 in Cdh5+ endothelial cells, Prx1+ BM stromal cells, Osx+ osteoprogenitor cells, and Mx1+ hematopoietic cells and uncovered the role of Atf4 in niche cells and hematopoiesis. Intriguingly, depletion of Atf4 in niche cells did not affect hematopoiesis; however, Atf4-deficient hematopoietic cells exhibited erythroid differentiation defects, leading to hypoplastic anemia. Mechanistically, ATF4 mediated direct regulation of Rps19bp1 transcription, which is, in turn, involved in 40 S ribosomal subunit assembly to coordinate ribosome biogenesis and promote erythropoiesis. Finally, we demonstrate that under conditions of 5-fluorouracil-induced stress, Atf4 depletion impedes the recovery of hematopoietic lineages, which requires efficient ribosome biogenesis. Taken together, our findings highlight the indispensable role of the ATF4-RPS19BP1 axis in the regulation of erythropoiesis.

Suppression of PKR promotes network excitability and enhanced cognition by interferon-γ-mediated disinhibition.

The double-stranded RNA-activated protein kinase (PKR) was originally identified as a sensor of virus infection, but its function in the brain remains unknown. Here, we report that the lack of PKR enhances learning and memory in several behavioral tasks while increasing network excitability. In addition, loss of PKR increases the late phase of long-lasting synaptic potentiation (L-LTP) in hippocampal slices. These effects are caused by an interferon-γ (IFN-γ)-mediated selective reduction in GABAergic synaptic action. Together, our results reveal that PKR finely tunes the network activity that must be maintained while storing a given episode during learning. Because PKR activity is altered in several neurological disorders, this kinase presents a promising new target for the treatment of cognitive dysfunction. As a first step in this direction, we show that a selective PKR inhibitor replicates the Pkr(-/-) phenotype in WT mice, enhancing long-term memory storage and L-LTP.

Atf7ip and Setdb1 interaction orchestrates the hematopoietic stem and progenitor cell state with diverse lineage differentiation.

Hematopoietic stem and progenitor cells (HSPCs) are a heterogeneous group of cells with expansion, differentiation, and repopulation capacities. How HSPCs orchestrate the stemness state with diverse lineage differentiation at steady condition or acute stress remains largely unknown. Here, we show that zebrafish mutants that are deficient in an epigenetic regulator Atf7ip or Setdb1 methyltransferase undergo excessive myeloid differentiation with impaired HSPC expansion, manifesting a decline in T cells and erythroid lineage. We find that Atf7ip regulates hematopoiesis through Setdb1-mediated H3K9me3 modification and chromatin remodeling. During hematopoiesis, the interaction of Atf7ip and Setdb1 triggers H3K9me3 depositions in hematopoietic regulatory genes including cebpß and cdkn1a, preventing HSPCs from loss of expansion and premature differentiation into myeloid lineage. Concomitantly, loss of Atf7ip or Setdb1 derepresses retrotransposons that instigate the viral sensor Mda5/Rig-I like receptor (RLR) signaling, leading to stress-driven myelopoiesis and inflammation. We find that ATF7IP or SETDB1 depletion represses human leukemic cell growth and induces myeloid differentiation with retrotransposon-triggered inflammation. These findings establish that Atf7ip/Setdb1-mediated H3K9me3 deposition constitutes a genome-wide checkpoint that impedes the myeloid potential and maintains HSPC stemness for diverse blood cell production, providing unique insights into potential intervention in hematological malignancy.

The cryo-EM structure of the SF3b spliceosome complex bound to a splicing modulator reveals a pre-mRNA substrate competitive mechanism of action.

Somatic mutations in spliceosome proteins lead to dysregulated RNA splicing and are observed in a variety of cancers. These genetic aberrations may offer a potential intervention point for targeted therapeutics. SF3B1, part of the U2 small nuclear RNP (snRNP), is targeted by splicing modulators, including E7107, the first to enter clinical trials, and, more recently, H3B-8800. Modulating splicing represents a first-in-class opportunity in drug discovery, and elucidating the structural basis for the mode of action opens up new possibilities for structure-based drug design. Here, we present the cryogenic electron microscopy (cryo-EM) structure of the SF3b subcomplex (SF3B1, SF3B3, PHF5A, and SF3B5) bound to E7107 at 3.95 Å. This structure shows that E7107 binds in the branch point adenosine-binding pocket, forming close contacts with key residues that confer resistance upon mutation: SF3B1R1074H and PHF5AY36C The structure suggests a model in which splicing modulators interfere with branch point adenosine recognition and supports a substrate competitive mechanism of action (MOA). Using several related chemical probes, we validate the pose of the compound and support their substrate competitive MOA by comparing their activity against both strong and weak pre-mRNA substrates. Finally, we present functional data and structure-activity relationship (SAR) on the PHF5AR38C mutation that sensitizes cells to some chemical probes but not others. Developing small molecule splicing modulators represents a promising therapeutic approach for a variety of diseases, and this work provides a significant step in enabling structure-based drug design for these elaborate natural products. Importantly, this work also demonstrates that the utilization of cryo-EM in drug discovery is coming of age.

Structural basis of nucleosome dynamics modulation by histone variants H2A.B and H2A.Z.2.2.

Nucleosomes are dynamic entities with wide-ranging compositional variations. Human histone variants H2A.B and H2A.Z.2.2 play critical roles in multiple biological processes by forming unstable nucleosomes and open chromatin structures, but how H2A.B and H2A.Z.2.2 confer these dynamic features to nucleosomes remains unclear. Here, we report cryo-EM structures of nucleosome core particles containing human H2A.B (H2A.B-NCP) at atomic resolution, identifying large-scale structural rearrangements in the histone octamer in H2A.B-NCP. H2A.B-NCP compacts approximately 103 bp of DNA wrapping around the core histones in approximately 1.2 left-handed superhelical turns, in sharp contrast to canonical nucleosome encompassing approximately 1.7 turns of DNA. Micrococcal nuclease digestion assay reveals that nineteen H2A.B-specific residues, including a ROF ("regulating-octamer-folding") sequence of six consecutive residues, are responsible for loosening of H2A.B-NCPs. Unlike H2A.B-NCP, the H2A.Z.2.2-containing nucleosome (Z.2.2-NCP) adopts a less-extended structure and compacts around 125 bp of DNA. Further investigation uncovers a crucial role for the H2A.Z.2.2-specific ROF in both H2A.Z.2.2-NCP opening and SWR1-dependent histone replacement. Taken together, these first high-resolution structure of unstable nucleosomes induced by histone H2A variants elucidate specific functions of H2A.B and H2A.Z.2.2 in enhancing chromatin dynamics.

DNA methylation safeguards the generation of hematopoietic stem and progenitor cells by repression of Notch signaling.

The earliest hematopoietic stem and progenitor cells (HSPCs) are generated from the ventral wall of the dorsal aorta, through endothelial-to-hematopoietic transition during vertebrate embryogenesis. Notch signaling is crucial for HSPC generation across vertebrates; however, the precise control of Notch during this process remains unclear. In the present study, we used multi-omics approaches together with functional assays to assess global DNA methylome dynamics during the endothelial cells to HSPCs transition in zebrafish, and determined that DNA methyltransferase 1 (Dnmt1) is essential for HSPC generation via repression of Notch signaling. Depletion of dnmt1 resulted in decreased DNA methylation levels and impaired HSPC production. Mechanistically, we found that loss of dnmt1 induced hypomethylation of Notch genes and consequently elevated Notch activity in hemogenic endothelial cells, thereby repressing the generation of HSPCs. This finding deepens our understanding of HSPC specification in vivo, which will provide helpful insights for designing new strategies for HSPC generation in vitro.

A strategy combining denoising and cryo-EM single particle analysis.

In cryogenic electron microscopy (cryo-EM) single particle analysis (SPA), high-resolution three-dimensional structures of biological macromolecules are determined by iteratively aligning and averaging a large number of two-dimensional projections of molecules. Since the correlation measures are sensitive to the signal-to-noise ratio, various parameter estimation steps in SPA will be disturbed by the high-intensity noise in cryo-EM. However, denoising algorithms tend to damage high frequencies and suppress mid- and high-frequency contrast of micrographs, which exactly the precise parameter estimation relies on, therefore, limiting their application in SPA. In this study, we suggest combining a cryo-EM image processing pipeline with denoising and maximizing the signal's contribution in various parameter estimation steps. To solve the inherent flaws of denoising algorithms, we design an algorithm named MScale to correct the amplitude distortion caused by denoising and propose a new orientation determination strategy to compensate for the high-frequency loss. In the experiments on several real datasets, the denoised particles are successfully applied in the class assignment estimation and orientation determination tasks, ultimately enhancing the quality of biomacromolecule reconstruction. The case study on classification indicates that our strategy not only improves the resolution of difficult classes (up to 5 Å) but also resolves an additional class. In the case study on orientation determination, our strategy improves the resolution of the final reconstructed density map by 0.34 Å compared with conventional strategy. The code is available at https://github.com/zhanghui186/Mscale.

Lipoxin A4 modulates the PKA/CREB and NF-κB signaling pathway to mitigate chondrocyte catabolism and apoptosis in temporomandibular joint osteoarthritis.

Temporomandibular joint osteoarthritis (TMJ-OA) is characterized by the degradation of the extracellular matrix (ECM) in cartilage and the apoptosis of chondrocytes, which is caused by inflammation and disruptions of chondrocyte metabolism and inflammation. Lipoxin A4 (LXA4), a specialized pro-resolving mediator, has been shown to inhibit inflammation and regulate the balance between ECM synthesis and degradation. However, the therapeutic effects of LXA4 on TMJ-OA and its underlying mechanisms remain unclear. Interleukin-1 beta (IL-1ß)-induced chondrocyte and surgically induced TMJ-OA rat models were established in this study. The viability of chondrocytes treated with LXA4 was evaluated with the cell counting kit-8 (CCK-8) assay, while protein levels were assessed by western blot analysis, and the apoptosis rate was evaluated with terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick-end labelling (TUNEL) staining. Histological analysis was conducted to evaluate the impact of LXA4 on cartilage degradation in TMJ-OA rat models. In vitro, the qRT-PCR and western blot analysis demonstrated that LXA4 facilitated the upregulation of collagen proteins (Collagen II) and decreased expression of matrix metalloproteinases (MMP-3, and MMP-13) associated with ECM modulation. LXA4 enhanced the TMJ-OA chondrocyte viability and decreased apoptotic rate. In vivo, histology and immunohistochemistry (IHC) analysis revealed that intraperitoneal injection of LXA4 contributed to the amelioration of chondrocyte injuries and deceleration of TMJ-OA. Transcriptomic sequencing revealed that cAMP signaling pathway was up-regulated and NF-κB signaling pathway was down-regulated in LXA4 treated group. LXA4 inhibited the phosphorylation of P65 and inhibitor of nuclear factor kappa B (IκBα) proteins while enhancing the phosphorylation PKA and CREB. This study demonstrates the potential of LXA4 as a therapeutic agent for suppressing chondrocyte catabolism and apoptosis by increasing PKA/CREB activity and decreasing NF-κB signaling.

Boolean map and object reconcile as the unit of visual working memory.

The unit of visual working memory is a fundamental issue under debate in the fields of cognitive psychology and neuroscience, with some traditional research suggesting that it is an object, while other recent studies demonstrating that a Boolean map offers a better account. The controversy surrounding the unit of visual working memory often centers on the representation of objects consist of same dimensional features (e.g. bicolor objects). For 2 colors in a bicolor object, some behavioral studies have suggested that they need to be represented by separate units, while some other studies using electrophysiological measures have found that they can be represented within a single unit. This apparent conflict hints that Boolean map and object may reconcile as the unit of visual working memory. Adopting the contralateral delay activity as an electrophysiological marker of visual working memory, experiments 1 and 2 consistently found that the contralateral delay activity amplitude for memorizing bicolor circles at P7/P8 conformed the Boolean map-based storage throughout the whole maintenance, while the contralateral delay activity amplitude at P3/P4 just conformed the object-based storage during the early period. It suggests though Boolean map got stronger supporting evidence than object, they 2 may coexist as the unit of visual working memory.

Mechanical cloak via data-driven aperiodic metamaterial design.

SignificanceAn invisibility cloak to conceal objects from an outside observer has long been a subject of interest in metamaterial design. While cloaks have been manufactured for optical, thermal, and electric fields, limited progress has been made for mechanical cloaks. Most existing designs rely on mapping-based methods, which have so far been limited to special base cells and a narrow selection of voids with simple shapes. In this study, we develop a fundamentally different approach by exploiting data-driven designs to offer timely, customized solutions to mechanical cloaking that were previously difficult to obtain. Through simulations and experimental validations, we show that excellent cloaking performance can be achieved for various boundary conditions, shapes of voids, base cells, and even multiple voids.

In situ structures of polymerase complex of mammalian reovirus illuminate RdRp activation and transcription regulation.

Mammalian reovirus (reovirus) is a multilayered, turreted member of Reoviridae characterized by transcription of dsRNA genome within the innermost capsid shell. Here, we present high-resolution in situ structures of reovirus transcriptase complex in an intact double-layered virion, and in the uncoated single-layered core particles in the unloaded, reloaded, pre-elongation, and elongation states, respectively, obtained by cryo-electron microscopy and sub-particle reconstructions. At the template entry of RNA-dependent RNA polymerase (RdRp), the RNA-loading region gets flexible after uncoating resulting in the unloading of terminal genomic RNA and inactivity of transcription. However, upon adding transcriptional substrates, the RNA-loading region is recovered leading the RNAs loaded again. The priming loop in RdRp was found to play a critical role in regulating transcription, which hinders the elongation of transcript in virion and triggers the rearrangement of RdRp C-terminal domain (CTD) during elongation, resulting in splitting of template-transcript hybrid and opening of transcript exit. With the integration of these structures, a transcriptional model of reovirus with five states is proposed. Our structures illuminate the RdRp activation and regulation of the multilayered turreted reovirus.

tRNA-derived fragments: Unveiling new roles and molecular mechanisms in cancer progression.

tRNA-derived fragments (tRFs) are novel small noncoding RNAs (sncRNAs) that range from approximately 14 to 50 nt. They are generated by the cleavage of mature tRNAs or precursor tRNAs (pre-tRNAs) at specific sites. Based on their origin and length, tRFs can be classified into three categories: (1) tRF-1 s; (2) tRF-3 s, tRF-5 s, and internal tRFs (i-tRFs); and (3) tRNA halves. They play important roles in stress response, signal transduction, and gene expression processes. Recent studies have identified differential expression of tRFs in various tumors. Aberrantly expressed tRFs have critical clinical value and show promise as new biomarkers for tumor diagnosis and prognosis and as therapeutic targets. tRFs regulate the malignant progression of tumors via various mechanisms, primarily including modulation of noncoding RNA biogenesis, global chromatin organization, gene expression regulation, modulation of protein translation, regulation of epigenetic modification, and alternative splicing regulation. In conclusion, tRF-mediated regulatory pathways could present new avenues for tumor treatment, and tRFs could serve as promising therapeutic targets for cancer therapy.

Alignment of single-cell RNA-seq samples without overcorrection using kernel density matching.

Single-cell RNA sequencing (scRNA-seq) technology is poised to replace bulk cell RNA sequencing for many biological and medical applications as it allows users to measure gene expression levels in a cell type-specific manner. However, data produced by scRNA-seq often exhibit batch effects that can be specific to a cell type, to a sample, or to an experiment, which prevent integration or comparisons across multiple experiments. Here, we present Dmatch, a method that leverages an external expression atlas of human primary cells and kernel density matching to align multiple scRNA-seq experiments for downstream biological analysis. Dmatch facilitates alignment of scRNA-seq data sets with cell types that may overlap only partially and thus allows integration of multiple distinct scRNA-seq experiments to extract biological insights. In simulation, Dmatch compares favorably to other alignment methods, both in terms of reducing sample-specific clustering and in terms of avoiding overcorrection. When applied to scRNA-seq data collected from clinical samples in a healthy individual and five autoimmune disease patients, Dmatch enabled cell type-specific differential gene expression comparisons across biopsy sites and disease conditions and uncovered a shared population of pro-inflammatory monocytes across biopsy sites in RA patients. We further show that Dmatch increases the number of eQTLs mapped from population scRNA-seq data. Dmatch is fast, scalable, and improves the utility of scRNA-seq for several important applications. Dmatch is freely available online.