Correlated Protein Modules Revealing Functional Coordination of Interacting Proteins Are Detected by Single-Cell Proteomics.

Single-cell proteomics has attracted a lot of attention in recent years because it offers more functional relevance than single-cell transcriptomics. However, most work to date has focused on cell typing, which has been widely accomplished by single-cell transcriptomics. Here we report the use of single-cell proteomics to measure the correlation between the translational levels of a pair of proteins in a single mammalian cell. In measuring pairwise correlations among ∼1000 proteins in a population of homogeneous K562 cells under a steady-state condition, we observed multiple correlated protein modules (CPMs), each containing a group of highly positively correlated proteins that are functionally interacting and collectively involved in certain biological functions, such as protein synthesis and oxidative phosphorylation. Some CPMs are shared across different cell types while others are cell-type specific. Widely studied in omics analyses, pairwise correlations are often measured by introducing perturbations into bulk samples. However, some correlations of gene or protein expression under the steady-state condition would be masked by perturbation. The single-cell correlations probed in our experiment reflect intrinsic steady-state fluctuations in the absence of perturbation. We note that observed correlations between proteins are experimentally more distinct and functionally more relevant than those between corresponding mRNAs measured in single-cell transcriptomics. By virtue of single-cell proteomics, functional coordination of proteins is manifested through CPMs.

Highly sensitive single-cell chromatin accessibility assay and transcriptome coassay with METATAC.

Recent advances in single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) and its coassays have transformed the field of single-cell epigenomics and transcriptomics. However, the low detection efficiency of current methods has limited our understanding of the true complexity of chromatin accessibility and its relationship with gene expression in single cells. Here, we report a high-sensitivity scATAC-seq method, termed multiplexed end-tagging amplification of transposase accessible chromatin (METATAC), which detects a large number of accessible sites per cell and is compatible with automation. Our high detectability and statistical framework allowed precise linking of enhancers to promoters without merging single cells. We systematically investigated allele-specific accessibility in the mouse cerebral cortex, revealing allele-specific accessibility of promotors of certain imprinted genes but biallelic accessibility of their enhancers. Finally, we combined METATAC with our high-sensitivity single-cell RNA sequencing (scRNA-seq) method, multiple annealing and looping based amplification cycles for digital transcriptomics (MALBAC-DT), to develop a joint ATAC-RNA assay, termed METATAC and MALBAC-DT coassay by sequencing (M2C-seq). M2C-seq achieved significant improvements for both ATAC and RNA compared with previous methods, with consistent performance across cell lines and early mouse embryos.

Protein and lipid mass concentration measurement in tissues by stimulated Raman scattering microscopy.

Cell mass and chemical composition are important aggregate cellular properties that are especially relevant to physiological processes, such as growth control and tissue homeostasis. Despite their importance, it has been difficult to measure these features quantitatively at the individual cell level in intact tissue. Here, we introduce normalized Raman imaging (NoRI), a stimulated Raman scattering (SRS) microscopy method that provides the local concentrations of protein, lipid, and water from live or fixed tissue samples with high spatial resolution. Using NoRI, we demonstrate that protein, lipid, and water concentrations at the single cell are maintained in a tight range in cells under the same physiological conditions and are altered in different physiological states, such as cell cycle stages, attachment to substrates of different stiffness, or by entering senescence. In animal tissues, protein and lipid concentration varies with cell types, yet an unexpected cell-to-cell heterogeneity was found in cerebellar Purkinje cells. The protein and lipid concentration profile provides means to quantitatively compare disease-related pathology, as demonstrated using models of Alzheimer's disease. This demonstration shows that NoRI is a broadly applicable technique for probing the biological regulation of protein mass, lipid mass, and water mass for studies of cellular and tissue growth, homeostasis, and disease.

Coordination of two enhancers drives expression of olfactory trace amine-associated receptors.

Olfactory sensory neurons (OSNs) are functionally defined by their expression of a unique odorant receptor (OR). Mechanisms underlying singular OR expression are well studied, and involve a massive cross-chromosomal enhancer interaction network. Trace amine-associated receptors (TAARs) form a distinct family of olfactory receptors, and here we find that mechanisms regulating Taar gene choice display many unique features. The epigenetic signature of Taar genes in TAAR OSNs is different from that in OR OSNs. We further identify that two TAAR enhancers conserved across placental mammals are absolutely required for expression of the entire Taar gene repertoire. Deletion of either enhancer dramatically decreases the expression probabilities of different Taar genes, while deletion of both enhancers completely eliminates the TAAR OSN populations. In addition, both of the enhancers are sufficient to drive transgene expression in the partially overlapped TAAR OSNs. We also show that the TAAR enhancers operate in cis to regulate Taar gene expression. Our findings reveal a coordinated control of Taar gene choice in OSNs by two remote enhancers, and provide an excellent model to study molecular mechanisms underlying formation of an olfactory subsystem.

Accurate SNV detection in single cells by transposon-based whole-genome amplification of complementary strands.

Single-nucleotide variants (SNVs), pertinent to aging and disease, occur sporadically in the human genome, hence necessitating single-cell measurements. However, detection of single-cell SNVs suffers from false positives (FPs) due to intracellular single-stranded DNA damage and the process of whole-genome amplification (WGA). Here, we report a single-cell WGA method termed multiplexed end-tagging amplification of complementary strands (META-CS), which eliminates nearly all FPs by virtue of DNA complementarity, and achieved the highest accuracy thus far. We validated META-CS by sequencing kindred cells and human sperm, and applied it to other human tissues. Investigation of mature single human neurons revealed increasing SNVs with age and potentially unrepaired strand-specific oxidative guanine damage. We determined SNV frequencies along the genome in differentiated single human blood cells, and identified cell type-dependent mutational patterns for major types of lymphocytes.

Changes in genome architecture and transcriptional dynamics progress independently of sensory experience during post-natal brain development.

Both transcription and three-dimensional (3D) architecture of the mammalian genome play critical roles in neurodevelopment and its disorders. However, 3D genome structures of single brain cells have not been solved; little is known about the dynamics of single-cell transcriptome and 3D genome after birth. Here, we generated a transcriptome (3,517 cells) and 3D genome (3,646 cells) atlas of the developing mouse cortex and hippocampus by using our high-resolution multiple annealing and looping-based amplification cycles for digital transcriptomics (MALBAC-DT) and diploid chromatin conformation capture (Dip-C) methods and developing multi-omic analysis pipelines. In adults, 3D genome "structure types" delineate all major cell types, with high correlation between chromatin A/B compartments and gene expression. During development, both transcriptome and 3D genome are extensively transformed in the first post-natal month. In neurons, 3D genome is rewired across scales, correlated with gene expression modules, and independent of sensory experience. Finally, we examine allele-specific structure of imprinted genes, revealing local and chromosome (chr)-wide differences. These findings uncover an unknown dimension of neurodevelopment.

MINERVA: A Facile Strategy for SARS-CoV-2 Whole-Genome Deep Sequencing of Clinical Samples.

Analyzing the genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from clinical samples is crucial for understanding viral spread and evolution as well as for vaccine development. Existing RNA sequencing methods are demanding on user technique and time and, thus, not ideal for time-sensitive clinical samples; these methods are also not optimized for high performance on viral genomes. We developed a facile, practical, and robust approach for metagenomic and deep viral sequencing from clinical samples. We demonstrate the utility of our approach on pharyngeal, sputum, and stool samples collected from coronavirus disease 2019 (COVID-19) patients, successfully obtaining whole metatranscriptomes and complete high-depth, high-coverage SARS-CoV-2 genomes with high yield and robustness. With a shortened hands-on time from sample to virus-enriched sequencing-ready library, this rapid, versatile, and clinic-friendly approach will facilitate molecular epidemiology studies during current and future outbreaks.

Structurally Resolved SARS-CoV-2 Antibody Shows High Efficacy in Severely Infected Hamsters and Provides a Potent Cocktail Pairing Strategy.

Understanding how potent neutralizing antibodies (NAbs) inhibit SARS-CoV-2 is critical for effective therapeutic development. We previously described BD-368-2, a SARS-CoV-2 NAb with high potency; however, its neutralization mechanism is largely unknown. Here, we report the 3.5-Å cryo-EM structure of BD-368-2/trimeric-spike complex, revealing that BD-368-2 fully blocks ACE2 recognition by occupying all three receptor-binding domains (RBDs) simultaneously, regardless of their "up" or "down" conformations. Also, BD-368-2 treats infected adult hamsters at low dosages and at various administering windows, in contrast to placebo hamsters that manifested severe interstitial pneumonia. Moreover, BD-368-2's epitope completely avoids the common binding site of VH3-53/VH3-66 recurrent NAbs, evidenced by tripartite co-crystal structures with RBDs. Pairing BD-368-2 with a potent recurrent NAb neutralizes SARS-CoV-2 pseudovirus at pM level and rescues mutation-induced neutralization escapes. Together, our results rationalized a new RBD epitope that leads to high neutralization potency and demonstrated BD-368-2's therapeutic potential in treating COVID-19.

Potent Neutralizing Antibodies against SARS-CoV-2 Identified by High-Throughput Single-Cell Sequencing of Convalescent Patients' B Cells.

The COVID-19 pandemic urgently needs therapeutic and prophylactic interventions. Here, we report the rapid identification of SARS-CoV-2-neutralizing antibodies by high-throughput single-cell RNA and VDJ sequencing of antigen-enriched B cells from 60 convalescent patients. From 8,558 antigen-binding IgG1+ clonotypes, 14 potent neutralizing antibodies were identified, with the most potent one, BD-368-2, exhibiting an IC50 of 1.2 and 15 ng/mL against pseudotyped and authentic SARS-CoV-2, respectively. BD-368-2 also displayed strong therapeutic and prophylactic efficacy in SARS-CoV-2-infected hACE2-transgenic mice. Additionally, the 3.8 Å cryo-EM structure of a neutralizing antibody in complex with the spike-ectodomain trimer revealed the antibody's epitope overlaps with the ACE2 binding site. Moreover, we demonstrated that SARS-CoV-2-neutralizing antibodies could be directly selected based on similarities of their predicted CDR3H structures to those of SARS-CoV-neutralizing antibodies. Altogether, we showed that human neutralizing antibodies could be efficiently discovered by high-throughput single B cell sequencing in response to pandemic infectious diseases.

RNA sequencing by direct tagmentation of RNA/DNA hybrids.

Transcriptome profiling by RNA sequencing (RNA-seq) has been widely used to characterize cellular status, but it relies on second-strand complementary DNA (cDNA) synthesis to generate initial material for library preparation. Here we use bacterial transposase Tn5, which has been increasingly used in various high-throughput DNA analyses, to construct RNA-seq libraries without second-strand synthesis. We show that Tn5 transposome can randomly bind RNA/DNA heteroduplexes and add sequencing adapters onto RNA directly after reverse transcription. This method, Sequencing HEteRo RNA-DNA-hYbrid (SHERRY), is versatile and scalable. SHERRY accepts a wide range of starting materials, from bulk RNA to single cells. SHERRY offers a greatly simplified protocol and produces results with higher reproducibility and GC uniformity compared with prevailing RNA-seq methods.

Three-dimensional genome structures of single sensory neurons in mouse visual and olfactory systems.

Sensory neurons in the mouse eye and nose have unusual chromatin organization. Here we report their three-dimensional (3D) genome structure at 20-kilobase (kb) resolution, achieved by applying our recently developed diploid chromatin conformation capture (Dip-C) method to 409 single cells from the retina and the main olfactory epithelium of adult and newborn mice. The 3D genome of rod photoreceptors exhibited inverted radial distribution of euchromatin and heterochromatin compared with that of other cell types, whose nuclear periphery is mainly heterochromatin. Such genome-wide inversion is not observed in olfactory sensory neurons (OSNs). However, OSNs exhibited an interior bias for olfactory receptor (OR) genes and enhancers, in clear contrast to non-neuronal cells. Each OSN harbored multiple aggregates of OR genes and enhancers from different chromosomes. We also observed structural heterogeneity of the protocadherin gene cluster. This type of genome organization may provide the structural basis of the 'one-neuron, one-receptor' rule of olfaction.

Multiregion Sequencing Reveals the Genetic Heterogeneity and Evolutionary History of Osteosarcoma and Matched Pulmonary Metastases.

Osteosarcoma is the most common primary bone malignancy, and the lung is the most frequent site of metastasis. The limited understanding of the tumoral heterogeneity and evolutionary process of genomic alterations in pulmonary metastatic osteosarcoma impedes development of novel therapeutic strategies. Here we systematically illustrate the genomic disparities between primary tumors and corresponding pulmonary metastatic tumors by multiregional whole-exome and whole-genome sequencing in 86 tumor regions from 10 patients with osteosarcoma. Metastatic tumors exhibited a significantly higher mutational burden and genomic instability compared with primary tumors, possibly due to accumulation of mutations caused by a greater number of alterations in DNA damage response genes in metastatic tumors. Integrated analysis of the architecture and relationships of subclones revealed a dynamic mutational process and diverse dissemination patterns of osteosarcoma during pulmonary metastasis (6/10 with linear and 4/10 with parallel evolutionary patterns). All patients demonstrated more significant intertumoral rather than intratumoral heterogeneity between primary tumors and metastatic tumors. Mutated genes were enriched in the PI3K-Akt pathway at both the early and late stages of tumor evolution and in the MAPK pathway at the metastatic stage. Conversely, metastatic tumors showed improved immunogenicity, including higher neoantigen load, elevated PD-L1 expression, and tumor-infiltrating lymphocytes than the corresponding primary tumors. Our study is the first to report the dynamic evolutionary process and temporospatial tumor heterogeneity of pulmonary metastatic osteosarcoma, providing new insights for diagnosis and potential therapeutic strategies for pulmonary metastasis. SIGNIFICANCE: High-throughput sequencing of primary and metastatic osteosarcoma provides new insights into the diagnosis of and potential clinical therapeutic strategies for pulmonary metastasis.

Label-free imaging of amyloid plaques in Alzheimer's disease with stimulated Raman scattering microscopy.

One of the key pathological features of Alzheimer's disease (AD) is the existence of extracellular deposition of amyloid plaques formed with misfolded amyloid-ß (Aß). The conformational change of proteins leads to enriched contents of ß sheets, resulting in remarkable changes of vibrational spectra, especially the spectral shifts of the amide I mode. Here, we applied stimulated Raman scattering (SRS) microscopy to image amyloid plaques in the brain tissue of an AD mouse model. We have demonstrated the capability of SRS microscopy as a rapid, label-free imaging modality to differentiate misfolded from normal proteins based on the blue shift (~10 cm-1) of amide I SRS spectra. Furthermore, SRS imaging of Aß plaques was verified by antibody staining of frozen thin sections and fluorescence imaging of fresh tissues. Our method may provide a new approach for studies of AD pathology, as well as other neurodegenerative diseases associated with protein misfolding.

Integrated mapping of pharmacokinetics and pharmacodynamics in a patient-derived xenograft model of glioblastoma.

Therapeutic options for the treatment of glioblastoma remain inadequate despite concerted research efforts in drug development. Therapeutic failure can result from poor permeability of the blood-brain barrier, heterogeneous drug distribution, and development of resistance. Elucidation of relationships among such parameters could enable the development of predictive models of drug response in patients and inform drug development. Complementary analyses were applied to a glioblastoma patient-derived xenograft model in order to quantitatively map distribution and resulting cellular response to the EGFR inhibitor erlotinib. Mass spectrometry images of erlotinib were registered to histology and magnetic resonance images in order to correlate drug distribution with tumor characteristics. Phosphoproteomics and immunohistochemistry were used to assess protein signaling in response to drug, and integrated with transcriptional response using mRNA sequencing. This comprehensive dataset provides simultaneous insight into pharmacokinetics and pharmacodynamics and indicates that erlotinib delivery to intracranial tumors is insufficient to inhibit EGFR tyrosine kinase signaling.

Three-dimensional genome structures of single diploid human cells.

Three-dimensional genome structures play a key role in gene regulation and cell functions. Characterization of genome structures necessitates single-cell measurements. This has been achieved for haploid cells but has remained a challenge for diploid cells. We developed a single-cell chromatin conformation capture method, termed Dip-C, that combines a transposon-based whole-genome amplification method to detect many chromatin contacts, called META (multiplex end-tagging amplification), and an algorithm to impute the two chromosome haplotypes linked by each contact. We reconstructed the genome structures of single diploid human cells from a lymphoblastoid cell line and from primary blood cells with high spatial resolution, locating specific single-nucleotide and copy number variations in the nucleus. The two alleles of imprinted loci and the two X chromosomes were structurally different. Cells of different types displayed statistically distinct genome structures. Such structural cell typing is crucial for understanding cell functions.

Genotyping single-sperm cells by universal MARSALA enables the acquisition of linkage information for combined pre-implantation genetic diagnosis and genome screening.

PURPOSE: This paper aims to investigate the feasibility of performing pre-implantation genetic diagnosis (PGD) and pre-implantation genetic screening (PGS) simultaneously by a universal strategy without the requirement of genotyping relevant affected family members or lengthy preliminary work on linkage analysis. METHODS: By utilizing a universal Mutated Allele Revealed by Sequencing with Aneuploidy and Linkage Analyses (MARSALA) strategy based on low depth whole genome sequencing (~3x), not involving specific primers' design nor the enrichment of SNP markers for haplotype construction. Single-sperm cells and trephectoderm cells from in vitro fertilized embryos from a couple carrying HBB mutations were genotyped. Haplotypes of paternal alleles were constructed and investigated in embryos, and the chromosome copy number profiles were simultaneously analyzed. RESULTS: The universal MARSALA strategy allows the selection of a euploid embryo free of disease mutations for in uterus transfer and successful pregnancy. A follow-up amniocentesis was performed at 17 weeks of gestation to confirm the PGD/PGS results. CONCLUSION: We present the first successful PGD procedure based on genotyping multiple single-sperm cells to obtain SNP linkage information. Our improved PGD/PGS procedure does not require genotyping the proband or relevant family members and therefore can be applicable to a wider population of patients when conducting PGD for monogenic disorders.

Highly accurate fluorogenic DNA sequencing with information theory-based error correction.

Eliminating errors in next-generation DNA sequencing has proved challenging. Here we present error-correction code (ECC) sequencing, a method to greatly improve sequencing accuracy by combining fluorogenic sequencing-by-synthesis (SBS) with an information theory-based error-correction algorithm. ECC embeds redundancy in sequencing reads by creating three orthogonal degenerate sequences, generated by alternate dual-base reactions. This is similar to encoding and decoding strategies that have proved effective in detecting and correcting errors in information communication and storage. We show that, when combined with a fluorogenic SBS chemistry with raw accuracy of 98.1%, ECC sequencing provides single-end, error-free sequences up to 200 bp. ECC approaches should enable accurate identification of extremely rare genomic variations in various applications in biology and medicine.

Mapping allele with resolved carrier status of Robertsonian and reciprocal translocation in human preimplantation embryos.

Reciprocal translocations (RecT) and Robertsonian translocations (RobT) are among the most common chromosomal abnormalities that cause infertility and birth defects. Preimplantation genetic testing for aneuploidy using comprehensive chromosome screening for in vitro fertilization enables embryo selection with balanced chromosomal ploidy; however, it is normally unable to determine whether an embryo is a translocation carrier. Here we report a method named "Mapping Allele with Resolved Carrier Status" (MaReCs), which enables chromosomal ploidy screening and resolution of the translocation carrier status of the same embryo. We performed MaReCs on 108 embryos, of which 96 were from 13 RecT carriers and 12 were from three RobT carriers. Thirteen of the sixteen patients had at least one diploid embryo. We have confirmed the accuracy of our carrier status determination in amniotic fluid karyotyping of seven cases as well as in the live birth we have thus far. Therefore, MaReCs accurately enables the selection of translocation-free embryos from patients carrying chromosomal translocations. We expect MaReCs will help reduce the propagation of RecT/RobT in the human population.

Rapid intraoperative histology of unprocessed surgical specimens via fibre-laser-based stimulated Raman scattering microscopy.

Conventional methods for intraoperative histopathologic diagnosis are labour- and time-intensive, and may delay decision-making during brain-tumour surgery. Stimulated Raman scattering (SRS) microscopy, a label-free optical process, has been shown to rapidly detect brain-tumour infiltration in fresh, unprocessed human tissues. Here, we demonstrate the first application of SRS microscopy in the operating room by using a portable fibre-laser-based microscope and unprocessed specimens from 101 neurosurgical patients. We also introduce an image-processing method - stimulated Raman histology (SRH) - which leverages SRS images to create virtual haematoxylin-and-eosin-stained slides, revealing essential diagnostic features. In a simulation of intraoperative pathologic consultation in 30 patients, we found a remarkable concordance of SRH and conventional histology for predicting diagnosis (Cohen's kappa, κ > 0.89), with accuracy exceeding 92%. We also built and validated a multilayer perceptron based on quantified SRH image attributes that predicts brain-tumour subtype with 90% accuracy. Our findings provide insight into how SRH can now be used to improve the surgical care of brain tumour patients.