Design of 240,000 orthogonal 25mer DNA barcode probes

Rapid identification of bacterial isolates using microfluidic adaptive channels and multiplexed fluorescence microscopy

We demonstrate the rapid capture, enrichment, and identification of bacterial pathogens using Adaptive Channel Bacterial Capture (ACBC) devices. Using controlled tuning of device backpressure in polydimethylsiloxane (PDMS) devices, we enable the controlled formation of capture regions capable of trapping bacteria from low cell density samples with near 100% capture efficiency. The technical demands to prepare such devices are much lower compared to conventional methods for bacterial trapping and can be achieved with simple benchtop fabrication methods. We demonstrate the capture and identification of seven species of bacteria with bacterial concentrations lower than 1000 cells per mL, including common Gram-negative and Gram-positive pathogens such as Escherichia coli and Staphylococcus aureus. We further demonstrate that species identification of the trapped bacteria can be undertaken in the order of one-hour using multiplexed 16S rRNA-FISH with identification accuracies of 70-98% with unsupervised classification methods across 7 species of bacteria. Finally, by using the bacterial capture capabilities of the ACBC chip with an ultra-rapid antimicrobial susceptibility testing method employing fluorescence imaging and convolutional neural network (CNN) classification, we demonstrate that we can use the ACBC chip as an imaging flow cytometer that can predict the antibiotic susceptibility of E. coli cells after identification.

Protocol for simultaneous detection of mRNAs and (phospho-)proteins with ARTseq-FISH in mouse embryonic stem cells

Understanding the molecular signatures of individual cells within complex biological systems is crucial for deciphering cellular heterogeneity and uncovering regulatory mechanisms. Here, we present a protocol for simultaneous multiplexed detection of selected mRNAs and (phospho-)proteins in mouse embryonic stem cells using spatial single-cell profiling. We describe steps for employing single-stranded DNA (ssDNA)-labeled antibo'dies, padlock probes, and rolling circle amplification to achieve simultaneous visualization of mRNAs and (phospho-)proteins at subcellular resolution. This protocol has potential application in identifying cells in heterogeneous biological microenvironments. For complete details on the use and execution of this protocol, please refer to Hu et al.1.

Advances and Challenges in Random Access Techniques for In Vitro DNA Data Storage

With digital transformation and the general application of new technologies, data storage is facing new challenges with the demand for high-density loading of massive information. In response, DNA storage technology has emerged as a promising research direction. Efficient and reliable data retrieval is critical for DNA storage, and the development of random access technology plays a key role in its practicality and reliability. However, achieving fast and accurate random access functions has proven difficult for existing DNA storage efforts, which limits its practical applications in industry. In this review, we summarize the recent advances in DNA storage technology that enable random access functionality, as well as the challenges that need to be overcome and the current solutions. This review aims to help researchers in the field of DNA storage better understand the importance of the random access step and its impact on the overall development of DNA storage. Furthermore, the remaining challenges and future research trends in random access technology of DNA storage are discussed, with the goal of providing a solid foundation for achieving random access in DNA storage under large-scale data conditions.

Generation of densely labeled oligonucleotides for the detection of small genomic elements

The genome contains numerous regulatory elements that may undergo complex interactions and contribute to the establishment, maintenance, and change of cellular identity. Three-dimensional genome organization can be explored with fluorescence in situ hybridization (FISH) at the single-cell level, but the detection of small genomic loci remains challenging. Here, we provide a rapid and simple protocol for the generation of bright FISH probes suited for the detection of small genomic elements. We systematically optimized probe design and synthesis, screened polymerases for their ability to incorporate dye-labeled nucleotides, and streamlined purification conditions to yield nanoscopy-compatible oligonucleotides with dyes in variable arrays (NOVA probes). With these probes, we detect genomic loci ranging from genome-wide repetitive regions down to non-repetitive loci below the kilobase scale. In conclusion, we introduce a simple workflow to generate densely labeled oligonucleotide pools that facilitate detection and nanoscopic measurements of small genomic elements in single cells.

Scalable design of orthogonal DNA barcode libraries

Orthogonal DNA barcode library design is an essential task in bioengineering. Here we present seqwalk, an efficient method for designing barcode libraries that satisfy a sequence symmetry minimization (SSM) heuristic for orthogonality, with theoretical guarantees of maximal or near-maximal library size under certain design constraints. Seqwalk encodes SSM constraints in a de Bruijn graph representation of sequence space, enabling the application of recent advances in discrete mathematics1 to the problem of orthogonal sequence design. We demonstrate the scalability of seqwalk by designing a library of >106 SSM-satisfying barcode sequences in less than 20 s on a standard laptop.

An Effective DNA-Based File Storage System for Practical Archiving and Retrieval of Medical MRI Data

DNA-based data storage is a new technology in computational and synthetic biology, that offers a solution for long-term, high-density data archiving. Given the critical importance of medical data in advancing human health, there is a growing interest in developing an effective medical data storage system based on DNA. Data integrity, accuracy, reliability, and efficient retrieval are all significant concerns. Therefore, this study proposes an Effective DNA Storage (EDS) approach for archiving medical MRI data. The EDS approach incorporates three key components (i) a novel fraction strategy to address the critical issue of rotating encoding, which often leads to data loss due to single base error propagation; (ii) a novel rule-based quaternary transcoding method that satisfies bio-constraints and ensure reliable mapping; and (iii) an indexing technique designed to simplify random search and access. The effectiveness of this approach is validated through computer simulations and biological experiments, confirming its practicality. The EDS approach outperforms existing methods, providing superior control over bio-constraints and reducing computational time. The results and code provided in this study open new avenues for practical DNA storage of medical MRI data, offering promising prospects for the future of medical data archiving and retrieval.

Scalable search of massively pooled nucleic acid samples enabled by a molecular database query language

The surge in nucleic acid analytics requires scalable storage and retrieval systems akin to electronic databases used to organize digital data. Such a system could transform disease diagnosis, ecological preservation, and molecular surveillance of biothreats. Current storage systems use individual containers for nucleic acid samples, requiring single-sample retrieval that falls short compared with digital databases that allow complex and combinatorial data retrieval on aggregated data. Here, we leverage protective microcapsules with combinatorial DNA labeling that enables arbitrary retrieval on pooled biosamples analogous to Structured Query Languages. Ninety-six encapsulated pooled mock SARS-CoV-2 genomic samples barcoded with patient metadata are used to demonstrate queries with simultaneous matches to sample collection date ranges, locations, and patient health statuses, illustrating how such flexible queries can be used to yield immunological or epidemiological insights. The approach applies to any biosample database labeled with orthogonal barcodes, enabling complex post-hoc analysis, for example, to study global biothreat epidemiology.

Simultaneous analysis of pMHC binding and reactivity unveils virus-specific CD8 T cell immunity to a concise epitope set

CD8 T cells provide immunity to virus infection through recognition of epitopes presented by peptide major histocompatibility complexes (pMHCs). To establish a concise panel of widely recognized T cell epitopes from common viruses, we combined analysis of TCR down-regulation upon stimulation with epitope-specific enumeration based on barcode-labeled pMHC multimers. We assess CD8 T cell binding and reactivity for 929 previously reported epitopes in the context of 1 of 25 HLA alleles representing 29 viruses. The prevalence and magnitude of CD8 T cell responses were evaluated in 48 donors and reported along with 137 frequently recognized virus epitopes, many of which were underrepresented in the public domain. Eighty-four percent of epitope-specific CD8 T cell populations demonstrated reactivity to peptide stimulation, which was associated with effector and long-term memory phenotypes. Conversely, nonreactive T cell populations were associated primarily with naive phenotypes. Our analysis provides a reference map of epitopes for characterizing CD8 T cell responses toward common human virus infections.

Evaluation of restriction and Cas endonuclease kinetics using matrix-insensitive magnetic biosensors

The enzymatic actions of endonucleases in vivo can be altered due to bound substrates and differences in local environments, including enzyme concentration, pH, salinity, ionic strength, and temperature. Thus, accurate estimation of enzymatic reactions in vivo using matrix-dependent methods in solution can be challenging. Here, we report a matrix-insensitive magnetic biosensing platform that enables the measurement of endonuclease activity under different conditions with varying pH, salinity, ionic strength, and temperature. Using biosensor arrays and orthogonal pairs of oligonucleotides, we quantitatively characterized the enzymatic activity of EcoRI under different buffer conditions and in the presence of inhibitors. To mimic a more physiological environment, we monitored the sequence-dependent star activity of EcoRI under unconventional conditions. Furthermore, enzymatic activity was measured in cell culture media, saliva, and serum. Last, we estimated the effective cleavage rates of Cas12a on anchored single-strand DNAs using this platform, which more closely resembles in vivo settings. This platform will facilitate precise characterization of restriction and Cas endonucleases under various conditions.

Spatial proteomics in neurons at single-protein resolution

To understand biological processes, it is necessary to reveal the molecular heterogeneity of cells by gaining access to the location and interaction of all biomolecules. Significant advances were achieved by super-resolution microscopy, but such methods are still far from reaching the multiplexing capacity of proteomics. Here, we introduce secondary label-based unlimited multiplexed DNA-PAINT (SUM-PAINT), a high-throughput imaging method that is capable of achieving virtually unlimited multiplexing at better than 15 nm resolution. Using SUM-PAINT, we generated 30-plex single-molecule resolved datasets in neurons and adapted omics-inspired analysis for data exploration. This allowed us to reveal the complexity of synaptic heterogeneity, leading to the discovery of a distinct synapse type. We not only provide a resource for researchers, but also an integrated acquisition and analysis workflow for comprehensive spatial proteomics at single-protein resolution.

Tunable Dynamics in a Multistrain Transcriptional Pulse Generator

One challenge in synthetic biology is the tuning of regulatory components within gene circuits to elicit a specific behavior. This challenge becomes more difficult in synthetic microbial consortia since each strain's circuit must function at the intracellular level and their combination must operate at the population level. Here we demonstrate that circuit dynamics can be tuned in synthetic consortia through the manipulation of strain fractions within the community. To do this, we construct a microbial consortium comprised of three strains of engineered Escherichia coli that, when cocultured, use homoserine lactone-mediated intercellular signaling to create a multistrain incoherent type-1 feedforward loop (I1-FFL). Like naturally occurring I1-FFL motifs in gene networks, this engineered microbial consortium acts as a pulse generator of gene expression. We demonstrate that the amplitude of the pulse can be easily tuned by adjusting the relative population fractions of the strains. We also develop a mathematical model for the temporal dynamics of the microbial consortium. This model allows us to identify population fractions that produced desired pulse characteristics, predictions that were confirmed for all but extreme fractions. Our work demonstrates that intercellular gene circuits can be effectively tuned simply by adjusting the starting fractions of each strain in the consortium.

Cell-type-specific 3D-genome organization and transcription regulation in the brain

3D organization of the genome plays a critical role in regulating gene expression. However, it remains unclear how chromatin organization differs among different cell types in the brain. Here we used genome-scale DNA and RNA imaging to investigate 3D-genome organization in transcriptionally distinct cell types in the primary motor cortex of the mouse brain. We uncovered a wide spectrum of differences in the nuclear architecture and 3D-genome organization among different cell types, ranging from the physical size of the cell nucleus to the active-inactive chromatin compartmentalization and radial positioning of chromatin loci within the nucleus. These cell-type-dependent variations in nuclear architecture and chromatin organization exhibited strong correlation with both total transcriptional activity of the cell and transcriptional regulation of cell-type-specific marker genes. Moreover, we found that the methylated-DNA-binding protein MeCP2 regulates transcription in a divergent manner, depending on the nuclear radial positions of chromatin loci, through modulating active-inactive chromatin compartmentalization.

Hybridization Protection Reaction for Sensitive and Robust Gene Expression Profiling of Clinical Formalin-Fixed Paraffin-Embedded Samples

BACKGROUND

RNA profiling of formalin-fixed paraffin-embedded (FFPE) tumor tissues for the molecular diagnostics of disease prognosis or treatment response is often irreproducible and limited to a handful of biomarkers. This has led to an unmet need for robust multiplexed assays that can profile several RNA biomarkers of interest using a limited amount of specimen. Here, we describe hybridization protection reaction (HPR), which is a novel RNA profiling approach with high reproducibility.

METHODS

HPR assays were designed for multiple genes, including 10 radiosensitivity-associated genes, and compared with TaqMan assays. Performance was tested with synthetic RNA fragments, and the ability to analyze RNA was investigated in FPPE samples from 20 normal lung tissues, 40 lung cancer, and 30 esophageal cancer biopsies.

RESULTS

Experiments performed on 3 synthetic RNA fragments demonstrated a linear dynamic range of over 1000-fold with a replicate correlation coefficient of 0.99 and high analytical sensitivity between 3.2 to 10 000 pM. Comparison of HPR with standard quantitative reverse transcription polymerase chain reaction on FFPE specimens shows nonsignificant differences with > 99% confidence interval between 2 assays in transcript profiling of 91.7% of test transcripts. In addition, HPR was effectively applied to quantify transcript levels of 10 radiosensitivity-associated genes.

CONCLUSIONS

Overall, HPR is an alternative approach for RNA profiling with high sensitivity, reproducibility, robustness, and capability for molecular diagnostics in FFPE tumor biopsy specimens of lung and esophageal cancer.

Perturb-tracing enables high-content screening of multiscale 3D genome regulators

Three-dimensional (3D) genome organization becomes altered during development, aging, and disease1-23, but the factors regulating chromatin topology are incompletely understood and currently no technology can efficiently screen for new regulators of multiscale chromatin organization. Here, we developed an image-based high-content screening platform (Perturb-tracing) that combines pooled CRISPR screen, a new cellular barcode readout method (BARC-FISH), and chromatin tracing. We performed a loss-of-function screen in human cells, and visualized alterations to their genome organization from 13,000 imaging target-perturbation combinations, alongside perturbation-paired barcode readout in the same single cells. Using 1.4 million 3D positions along chromosome traces, we discovered tens of new regulators of chromatin folding at different length scales, ranging from chromatin domains and compartments to chromosome territory. A subset of the regulators exhibited 3D genome effects associated with loop-extrusion and A-B compartmentalization mechanisms, while others were largely unrelated to these known 3D genome mechanisms. We found that the ATP-dependent helicase CHD7, the loss of which causes the congenital neural crest syndrome CHARGE24 and a chromatin remodeler previously shown to promote local chromatin openness25-27, counter-intuitively compacts chromatin over long range in different genomic contexts and cell backgrounds including neural crest cells, and globally represses gene expression. The DNA compaction effect of CHD7 is independent of its chromatin remodeling activity and does not require other protein partners. Finally, we identified new regulators of nuclear architectures and found a functional link between chromatin compaction and nuclear shape. Altogether, our method enables scalable, high-content identification of chromatin and nuclear topology regulators that will stimulate new insights into the 3D genome functions, such as global gene and nuclear regulation, in health and disease.

Spatial atlas of the mouse central nervous system at molecular resolution

Spatially charting molecular cell types at single-cell resolution across the 3D volume is critical for illustrating the molecular basis of brain anatomy and functions. Single-cell RNA sequencing has profiled molecular cell types in the mouse brain1,2, but cannot capture their spatial organization. Here we used an in situ sequencing method, STARmap PLUS3,4, to profile 1,022 genes in 3D at a voxel size of 194 × 194 × 345 nm3, mapping 1.09 million high-quality cells across the adult mouse brain and spinal cord. We developed computational pipelines to segment, cluster and annotate 230 molecular cell types by single-cell gene expression and 106 molecular tissue regions by spatial niche gene expression. Joint analysis of molecular cell types and molecular tissue regions enabled a systematic molecular spatial cell-type nomenclature and identification of tissue architectures that were undefined in established brain anatomy. To create a transcriptome-wide spatial atlas, we integrated STARmap PLUS measurements with a published single-cell RNA-sequencing atlas1, imputing single-cell expression profiles of 11,844 genes. Finally, we delineated viral tropisms of a brain-wide transgene delivery tool, AAV-PHP.eB5,6. Together, this annotated dataset provides a single-cell resource that integrates the molecular spatial atlas, brain anatomy and the accessibility to genetic manipulation of the mammalian central nervous system.

Cytogenetic screening of a Canadian swine breeding nucleus using a newly developed karyotyping method named oligo-banding

BACKGROUND

The frequency of chromosomal rearrangements in Canadian breeding boars has been estimated at 0.91 to 1.64%. These abnormalities are widely recognized as a potential cause of subfertility in livestock production. Since artificial insemination is practiced in almost all intensive pig production systems, the use of elite boars carrying cytogenetic defects that have an impact on fertility can lead to major economic losses. To avoid keeping subfertile boars in artificial insemination centres and spreading chromosomal defects within populations, cytogenetic screening of boars is crucial. Different techniques are used for this purpose, but several issues are frequently encountered, i.e. environmental factors can influence the quality of results, the lack of genomic information outputted by these techniques, and the need for prior cytogenetic skills. The aim of this study was to develop a new pig karyotyping method based on fluorescent banding patterns.

RESULTS

The use of 207,847 specific oligonucleotides generated 96 fluorescent bands that are distributed across the 18 autosomes and the sex chromosomes. Tested alongside conventional G-banding, this oligo-banding method allowed us to identify four chromosomal translocations and a rare unbalanced chromosomal rearrangement that was not detected by conventional banding. In addition, this method allowed us to investigate chromosomal imbalance in spermatozoa.

CONCLUSIONS

The use of oligo-banding was found to be appropriate for detecting chromosomal aberrations in a Canadian pig nucleus and its convenient design and use make it an interesting tool for livestock karyotyping and cytogenetic studies.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights

BACKGROUND

Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining.

RESULTS

We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids.

CONCLUSIONS

Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

High-throughput functional profiling of the human fungal pathogen Candida albicans genome

Candida albicans is a major fungal pathogen of humans. Although its genome has been sequenced more than two decades ago, there are still over 4300 uncharacterized C. albicans genes. We previously generated an ORFeome as well as a collection of destination vectors to facilitate overexpression of C. albicans ORFs. Here, we report the construction of ∼2500 overexpression mutants and their evaluation by in vitro spotting on rich medium and in a liquid pool experiment in rich medium, allowing the identification of genes whose overexpression has a fitness cost. The candidates were further validated at the individual strain level. This new resource allows large-scale screens in different growth conditions to be performed routinely. Altogether, based on the concept of identifying functionally related genes by cluster analysis, the availability of this overexpression mutant collection will facilitate the characterization of gene functions in C. albicans.

High-multiplex tissue imaging in routine pathology-are we there yet?

High-multiplex tissue imaging (HMTI) approaches comprise several novel immunohistological methods that enable in-depth, spatial single-cell analysis. Over recent years, studies in tumor biology, infectious diseases, and autoimmune conditions have demonstrated the information gain accessible when mapping complex tissues with HMTI. Tumor biology has been a focus of innovative multiparametric approaches, as the tumor microenvironment (TME) contains great informative value for accurate diagnosis and targeted therapeutic approaches: unraveling the cellular composition and structural organization of the TME using sophisticated computational tools for spatial analysis has produced histopathologic biomarkers for outcomes in breast cancer, predictors of positive immunotherapy response in melanoma, and histological subgroups of colorectal carcinoma. Integration of HMTI technologies into existing clinical workflows such as molecular tumor boards will contribute to improve patient outcomes through personalized treatments tailored to the specific heterogeneous pathological fingerprint of cancer, autoimmunity, or infection. Here, we review the advantages and limitations of existing HMTI technologies and outline how spatial single-cell data can improve our understanding of pathological disease mechanisms and determinants of treatment success. We provide an overview of the analytic processing and interpretation and discuss how HMTI can improve future routine clinical diagnostic and therapeutic processes.

Activating and repressing gene expression between chromosomes during stochastic fate specification

DNA elements act across long genomic distances to regulate gene expression. During transvection in Drosophila, DNA elements on one allele of a gene act between chromosomes to regulate expression of the other allele. Little is known about the biological roles and developmental regulation of transvection. Here, we study the stochastic expression of spineless (ss) in photoreceptors in the fly eye to understand transvection. We determine a biological role for transvection in regulating expression of naturally occurring ss alleles. We identify DNA elements required for activating and repressing transvection. Different enhancers participate in transvection at different times during development to promote gene expression and specify cell fates. Bringing a silencer element on a heterologous chromosome into proximity with the ss locus "reconstitutes" the gene, leading to repression. Our studies show that transvection regulates gene expression via distinct DNA elements at specific timepoints in development, with implications for genome organization and architecture.

HT-smFISH: a cost-effective and flexible workflow for high-throughput single-molecule RNA imaging

The ability to visualize RNA in its native subcellular environment by using single-molecule fluorescence in situ hybridization (smFISH) has reshaped our understanding of gene expression and cellular functions. A major hindrance of smFISH is the difficulty to perform systematic experiments in medium- or high-throughput formats, principally because of the high cost of generating the individual fluorescent probe sets. Here, we present high-throughput smFISH (HT-smFISH), a simple and cost-efficient method for imaging hundreds to thousands of single endogenous RNA molecules in 96-well plates. HT-smFISH uses RNA probes transcribed in vitro from a large pool of unlabeled oligonucleotides. This allows the generation of individual probes for many RNA species, replacing commercial DNA probe sets. HT-smFISH thus reduces costs per targeted RNA compared with many smFISH methods and is easily scalable and flexible in design. We provide a protocol that combines oligo pool design, probe set generation, optimized hybridization conditions and guidelines for image acquisition and analysis. The pipeline requires knowledge of standard molecular biology tools, cell culture and fluorescence microscopy. It is achievable in ~20 d. In brief, HT-smFISH is tailored for medium- to high-throughput screens that image RNAs at single-molecule sensitivity.

Temporally dynamic antagonism between transcription and chromatin compaction controls stochastic photoreceptor specification in flies

Stochastic mechanisms diversify cell fates during development. How cells randomly choose between two or more fates remains poorly understood. In the Drosophila eye, the random mosaic of two R7 photoreceptor subtypes is determined by expression of the transcription factor Spineless (Ss). We investigated how cis-regulatory elements and trans factors regulate nascent transcriptional activity and chromatin compaction at the ss gene locus during R7 development. The ss locus is in a compact state in undifferentiated cells. An early enhancer drives transcription in all R7 precursors, and the locus opens. In differentiating cells, transcription ceases and the ss locus stochastically remains open or compacts. In SsON R7s, ss is open and competent for activation by a late enhancer, whereas in SsOFF R7s, ss is compact, and repression prevents expression. Our results suggest that a temporally dynamic antagonism, in which transcription drives large-scale decompaction and then compaction represses transcription, controls stochastic fate specification.

A transcriptomic axis predicts state modulation of cortical interneurons

Transcriptomics has revealed that cortical inhibitory neurons exhibit a great diversity of fine molecular subtypes1-6, but it is not known whether these subtypes have correspondingly diverse patterns of activity in the living brain. Here we show that inhibitory subtypes in primary visual cortex (V1) have diverse correlates with brain state, which are organized by a single factor: position along the main axis of transcriptomic variation. We combined in vivo two-photon calcium imaging of mouse V1 with a transcriptomic method to identify mRNA for 72 selected genes in ex vivo slices. We classified inhibitory neurons imaged in layers 1-3 into a three-level hierarchy of 5 subclasses, 11 types and 35 subtypes using previously defined transcriptomic clusters3. Responses to visual stimuli differed significantly only between subclasses, with cells in the Sncg subclass uniformly suppressed, and cells in the other subclasses predominantly excited. Modulation by brain state differed at all hierarchical levels but could be largely predicted from the first transcriptomic principal component, which also predicted correlations with simultaneously recorded cells. Inhibitory subtypes that fired more in resting, oscillatory brain states had a smaller fraction of their axonal projections in layer 1, narrower spikes, lower input resistance and weaker adaptation as determined in vitro7, and expressed more inhibitory cholinergic receptors. Subtypes that fired more during arousal had the opposite properties. Thus, a simple principle may largely explain how diverse inhibitory V1 subtypes shape state-dependent cortical processing.

Sensitive visualization of SARS-CoV-2 RNA with CoronaFISH

The current COVID-19 pandemic is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The positive-sense single-stranded RNA virus contains a single linear RNA segment that serves as a template for transcription and replication, leading to the synthesis of positive and negative-stranded viral RNA (vRNA) in infected cells. Tools to visualize vRNA directly in infected cells are critical to analyze the viral replication cycle, screen for therapeutic molecules, or study infections in human tissue. Here, we report the design, validation, and initial application of FISH probes to visualize positive or negative RNA of SARS-CoV-2 (CoronaFISH). We demonstrate sensitive visualization of vRNA in African green monkey and several human cell lines, in patient samples and human tissue. We further demonstrate the adaptation of CoronaFISH probes to electron microscopy. We provide all required oligonucleotide sequences, source code to design the probes, and a detailed protocol. We hope that CoronaFISH will complement existing techniques for research on SARS-CoV-2 biology and COVID-19 pathophysiology, drug screening, and diagnostics.

Machine Learning-Guided Systematic Search of DNA Sequences for Sorting Carbon Nanotubes

The prerequisite of utilizing DNA in sequence-dependent applications is to search specific sequences. Developing a strategy for efficient DNA sequence screening represents a grand challenge due to the countless possibilities of sequence combination. Herein, relying on sequence-dependent recognition between DNA and single-wall carbon nanotubes (SWCNTs), we demonstrate a method for systematic search of DNA sequences for sorting single-chirality SWCNTs. Different from previously documented empirical search, which has a low efficiency and accuracy, our approach combines machine learning and experimental investigation. The number of resolving sequences and the success rate of finding them are improved from ∼10² to ∼10³ and from ∼10% to >90%, respectively. Moreover, the resolving sequence patterns determined from 5-mer and 6-mer short sequences can be extended to sequence search in longer DNA subspaces.

Spatial transcriptomics using combinatorial fluorescence spectral and lifetime encoding, imaging and analysis

Multiplexed mRNA profiling in the spatial context provides new information enabling basic research and clinical applications. Unfortunately, existing spatial transcriptomics methods are limited due to either low multiplexing or complexity. Here, we introduce a spatialomics technology, termed Multi Omic Single-scan Assay with Integrated Combinatorial Analysis (MOSAICA), that integrates in situ labeling of mRNA and protein markers in cells or tissues with combinatorial fluorescence spectral and lifetime encoded probes, spectral and time-resolved fluorescence imaging, and machine learning-based decoding. We demonstrate MOSAICA's multiplexing scalability in detecting 10-plex targets in fixed colorectal cancer cells using combinatorial labeling of five fluorophores with facile error-detection and removal of autofluorescence. MOSAICA's analysis is strongly correlated with sequencing data (Pearson's r = 0.96) and was further benchmarked using RNAscopeTM and LGC StellarisTM. We further apply MOSAICA for multiplexed analysis of clinical melanoma Formalin-Fixed Paraffin-Embedded (FFPE) tissues. We finally demonstrate simultaneous co-detection of protein and mRNA in cancer cells.

Neoantigen-reactive CD8+ T cells affect clinical outcome of adoptive transfer with tumor-infiltrating lymphocytes in melanoma

BACKGROUND

Neoantigen-driven recognition and T cell-mediated killing contribute to tumor clearance following adoptive cell therapy (ACT) with Tumor-Infiltrating Lymphocytes (TILs). Yet, how diversity, frequency, and persistence of expanded neoepitope-specific CD8+ T cells derived from TIL infusion products affect patient outcome is not fully determined.

METHODS

Using barcoded pMHC multimers, we provide a comprehensive mapping of CD8+ T cells recognizing neoepitopes in TIL infusion products and blood samples from 26 metastatic mela-noma patients who received ACT.

RESULTS

We identified 106 neoepitopes within TIL infusion products corresponding to 1.8% of all predicted neoepitopes. We observed neoepitope-specific recognition to be virtually devoid in TIL infusion products given to patients with progressive disease outcome. Moreover, we found that the frequency of neoepitope-specific CD8+ T cells in TIL infusion products correlated with in-creased survival, and that detection of engrafted CD8+ T cells in post-treatment (i.e. originating from the TIL infusion product) were unique to responders of TIL-ACT. Finally, we found that a transcriptional signature for lymphocyte activity within the tumor microenvironment was associated with a higher frequency of neoepitope-specific CD8+ T cells in the infusion product.

CONCLUSIONS

These data support previous case studies of neoepitope-specific CD8+ T cells in melanoma, and indicate that successful TIL-ACT is associated with an expansion of neoepitope-specific CD8+ T cells.

FUNDING

NEYE Foundation; European Research Council; Lundbeck Foundation Fellowship; Carlsberg Foundation.

Scalable Nucleic Acid Storage and Retrieval Using Barcoded Microcapsules

Rapid advances in nucleic acid sequencing and synthesis technologies have spurred a major need to collect, store, and sequence the DNA and RNA from viral, bacterial, and mammalian sources and organisms. However, current approaches to storing nucleic acids rely on a low-temperature environment and require robotics for access, posing challenges for scalable and low-cost nucleic acid storage. Here, we present an alternative method for storing nucleic acids, termed Preservation and Access of Nucleic aciDs using barcOded micRocApsules (PANDORA). Nucleic acids spanning kilobases to gigabases and from different sources, including animals, bacteria, and viruses, are encapsulated into silica microcapsules to protect them from environmental denaturants at room temperature. Molecular barcodes attached to each microcapsule enable sample pooling and subsequent identification and retrieval using fluorescence-activated sorting. We demonstrate quantitative storage and rapid access to targeted nucleic acids from a pool emulating standard retrieval operations implemented in conventional storage systems, including recovery of 100,000-200,000 samples and Boolean logic selection using four unique barcodes. Quantitative polymerase chain reaction and short-read sequencing of the retrieved samples validated the sorting experiments and the integrity of the released nucleic acids. Our proposed approach offers a scalable long-term, room-temperature storage and retrieval of nucleic acids with high sample fidelity.

CRISPR-based peptide library display and programmable microarray self-assembly for rapid quantitative protein binding assays

CRISPR-inspired systems have been extensively developed for applications in genome editing and nucleic acid detection. Here, we introduce a CRISPR-based peptide display technology to facilitate customized, high-throughput in vitro protein interaction studies. We show that bespoke peptide libraries fused to catalytically inactive Cas9 (dCas9) and barcoded with unique single guide RNA (sgRNA) molecules self-assemble from a single mixed pool to programmable positions on a DNA microarray surface for rapid, multiplexed binding assays. We develop dCas9-displayed saturation mutagenesis libraries to characterize antibody-epitope binding for a commercial anti-FLAG monoclonal antibody and human serum antibodies. We also show that our platform can be used for viral epitope mapping and exhibits promise as a multiplexed diagnostics tool. Our CRISPR-based peptide display platform and the principles of complex library self-assembly using dCas9 could be adapted for rapid interrogation of varied customized protein libraries or biological materials assembly using DNA scaffolding.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Random access DNA memory using Boolean search in an archival file storage system

DNA is an ultrahigh-density storage medium that could meet exponentially growing worldwide demand for archival data storage if DNA synthesis costs declined sufficiently and if random access of files within exabyte-to-yottabyte-scale DNA data pools were feasible. Here, we demonstrate a path to overcome the second barrier by encapsulating data-encoding DNA file sequences within impervious silica capsules that are surface labelled with single-stranded DNA barcodes. Barcodes are chosen to represent file metadata, enabling selection of sets of files with Boolean logic directly, without use of amplification. We demonstrate random access of image files from a prototypical 2-kilobyte image database using fluorescence sorting with selection sensitivity of one in 10⁶ files, which thereby enables one in 10^6N selection capability using N optical channels. Our strategy thereby offers a scalable concept for random access of archival files in large-scale molecular datasets.

Direct oligonucleotide sequencing with nanopores

Third-generation DNA sequencing has enabled sequencing of long, unamplified DNA fragments with minimal steps. Direct sequencing of ssDNA or RNA gives valuable insights like base-level modifications, phosphoramidite synthesis yield estimates and strand quality analysis, without the need to add the complimentary strand. Direct sequencing of single-stranded nucleic acid species is challenging as they are non-compatible to the double-stranded sequencing adapters used by manufacturers. The MinION platform from Oxford Nanopore Technologies performs sequencing by passing single-strands of DNA through a layer of biological nanopore sensors; although sequencing is performed on single-strands, the recommended template by the manufacturer is double-stranded. We have identified that the MinION platform can perform sequencing of short, single-strand oligonucleotides directly without amplification or second-strand synthesis by performing a single annealing step before library preparation. Short 5' phosphorylated oligos when annealed to an adapter sequence can be directly sequenced in the 5' to 3' direction via nanopores. Adapter sequences were designed to bind to the 5' end of the oligos and to leave a 3' adenosine overhang after binding to their target. The 3' adenosine overhang of the adapter and the terminal phosphate makes the 5' end of the oligo analogous to an end-prepared dsDNA, rendering it compatible with ligation-based library preparation for sequencing. An oligo-pool containing 42,000, 120 nt orthogonal sequences was phosphorylated and sequenced using this method and ~90% of these sequences were recovered with high accuracy using BLAST. In the nanopore raw data, we have identified that empty signals can be wrongly identified as a valid read by the MinION platform and sometimes multiple signals containing several strands can be fused into a single raw sequence file due to segmentation faults in the software. This direct oligonucleotide sequencing method enables novel applications in DNA data storage systems where short oligonucleotides are the primary information carriers.

Direct oligonucleotide sequencing with nanopores

Third-generation DNA sequencing has enabled sequencing of long, unamplified DNA fragments with minimal steps. Direct sequencing of ssDNA or RNA gives valuable insights like base-level modifications, phosphoramidite synthesis yield estimates and strand quality analysis, without the need to add the complimentary strand. Direct sequencing of single-stranded nucleic acid species is challenging as they are non-compatible to the double-stranded sequencing adapters used by manufacturers. The MinION platform from Oxford Nanopore Technologies performs sequencing by passing single-strands of DNA through a layer of biological nanopore sensors; although sequencing is performed on single-strands, the recommended template by the manufacturer is double-stranded. We have identified that the MinION platform can perform sequencing of short, single-strand oligonucleotides directly without amplification or second-strand synthesis by performing a single annealing step before library preparation. Short 5' phosphorylated oligos when annealed to an adapter sequence can be directly sequenced in the 5' to 3' direction via nanopores. Adapter sequences were designed to bind to the 5' end of the oligos and to leave a 3' adenosine overhang after binding to their target. The 3' adenosine overhang of the adapter and the terminal phosphate makes the 5' end of the oligo analogous to an end-prepared dsDNA, rendering it compatible with ligation-based library preparation for sequencing. An oligo-pool containing 42,000, 120 nt orthogonal sequences was phosphorylated and sequenced using this method and ~90% of these sequences were recovered with high accuracy using BLAST. In the nanopore raw data, we have identified that empty signals can be wrongly identified as a valid read by the MinION platform and sometimes multiple signals containing several strands can be fused into a single raw sequence file due to segmentation faults in the software. This direct oligonucleotide sequencing method enables novel applications in DNA data storage systems where short oligonucleotides are the primary information carriers.

Resolving cell state in iPSC-derived human neural samples with multiplexed fluorescence imaging

Human induced pluripotent stem cell-derived (iPSC) neural cultures offer clinically relevant models of human diseases, including Amyotrophic Lateral Sclerosis, Alzheimer’s, and Autism Spectrum Disorder. In situ characterization of the spatial-temporal evolution of cell state in 3D culture and subsequent 2D dissociated culture models based on protein expression levels and localizations is essential to understanding neural cell differentiation, disease state phenotypes, and sample-to-sample variability. Here, we apply PRobe-based Imaging for Sequential Multiplexing (PRISM) to facilitate multiplexed imaging with facile, rapid exchange of imaging probes to analyze iPSC-derived cortical and motor neuron cultures that are relevant to psychiatric and neurodegenerative disease models, using over ten protein targets. Our approach permits analysis of cell differentiation, cell composition, and functional marker expression in complex stem-cell derived neural cultures. Furthermore, our approach is amenable to automation, offering in principle the ability to scale-up to dozens of protein targets and samples.

Tomov et al. utilize DNA-PRISM to allow for multiplexed imaging of cultured cells using antibodies modified with oligonucleotide probes. The differentiation of iPSCs to cortical and motor neurons is characterized in model cultures, relevant for use in disease research and drug screening.

Transcriptionally inactive hepatitis B virus episome DNA preferentially resides in the vicinity of chromosome 19 in 3D host genome upon infection

The hepatitis B virus (HBV) infects 257 million people worldwide. HBV infection requires establishment and persistence of covalently closed circular (ccc) DNA, a viral episome, in nucleus. Here, we study cccDNA spatial localization in the 3D host genome by using chromosome conformation capture-based sequencing analysis and fluorescence in situ hybridization (FISH). We show that transcriptionally inactive cccDNA is not randomly distributed in host nucleus. Rather, it is preferentially accumulated at specialized areas, including regions close to chromosome 19 (chr.19). Activation of the cccDNA is apparently associated with its re-localization, from a pre-established heterochromatin hub formed by 5 regions of chr.19 to transcriptionally active regions formed by chr.19 and nearby chromosomes including chr.16, 17, 20, and 22. This active versus inactive positioning at discrete regions of the host genome is primarily controlled by the viral HBx protein and by host factors including the structural maintenance of chromosomes protein 5/6 (SMC5/6) complex.

Dimerization and oligomerization of DNA-assembled building blocks for controlled multi-motion in high-order architectures

In living organisms, proteins are organized prevalently through a self-association mechanism to form dimers and oligomers, which often confer new functions at the intermolecular interfaces. Despite the progress on DNA-assembled artificial systems, endeavors have been largely paid to achieve monomeric nanostructures that mimic motor proteins for a single type of motion. Here, we demonstrate a DNA-assembled building block with rotary and walking modules, which can introduce new motion through dimerization and oligomerization. The building block is a chiral system, comprising two interacting gold nanorods to perform rotation and walking, respectively. Through dimerization, two building blocks can form a dimer to yield coordinated sliding. Further oligomerization leads to higher-order structures, containing alternating rotation and sliding dimer interfaces to impose structural twisting. Our hierarchical assembly scheme offers a design blueprint to construct DNA-assembled advanced architectures with high degrees of freedom to tailor the optical responses and regulate multi-motion on the nanoscale.

Creation of high-order architectures using DNA devices is of interest for increasing the complexity of synthetic systems. Here, the authors, inspired by biological oligomers, create DNA dimers and oligomers that combining rotation and walking to make high-order systems with more complex conformational changes.

Genome oligopaint via local denaturation fluorescence in situ hybridization

Cas9 in complex with a programmable guide RNA targets specific double-stranded DNA for cleavage. By harnessing Cas9 as a programmable loader of superhelicase to genomic DNA, we report a physiological-temperature DNA fluorescence in situ hybridization (FISH) method termed genome oligopaint via local denaturation (GOLD) FISH. Instead of global denaturation as in conventional DNA FISH, loading a superhelicase at a Cas9-generated nick allows for local DNA denaturation, reducing nonspecific binding of probes and avoiding harsh treatments such as heat denaturation. GOLD FISH relies on Cas9 cleaving target DNA sequences and avoids the high nuclear background associated with other genome labeling methods that rely on Cas9 binding. The excellent signal brightness and specificity enable us to image nonrepetitive genomic DNA loci and analyze the conformational differences between active and inactive X chromosomes. Finally, GOLD FISH could be used for rapid identification of HER2 gene amplification in patient tissue.

A large-scale survey of pairwise epistasis reveals a mechanism for evolutionary expansion and specialization of PDZ domains

Deep mutational scanning (DMS) facilitates data-driven models of protein structure and function. Here, we adapted Saturated Programmable Insertion Engineering (SPINE) as a programmable DMS technique. We validate SPINE with a reference single mutant dataset in the PSD95 PDZ3 domain and then characterize most pairwise double mutants to study epistasis. We observe wide-spread proximal negative epistasis, which we attribute to mutations affecting thermodynamic stability, and strong long-range positive epistasis, which is enriched in an evolutionarily conserved and function-defining network of "sector" and clade-specifying residues. Conditional neutrality of mutations in clade-specifying residues compensates for deleterious mutations in sector positions. This suggests that epistatic interactions between these position pairs facilitated the evolutionary expansion and specialization of PDZ domains. We propose that SPINE provides easy experimental access to reveal epistasis signatures in proteins that will improve our understanding of the structural basis for protein function and adaptation.

A comprehensive phenotypic CRISPR-Cas9 screen of the ubiquitin pathway uncovers roles of ubiquitin ligases in mitosis

The human ubiquitin proteasome system, composed of over 700 ubiquitin ligases (E3s) and deubiquitinases (DUBs), has been difficult to characterize systematically and phenotypically. We performed chemical-genetic CRISPR-Cas9 screens to identify E3s/DUBs whose loss renders cells sensitive or resistant to 41 compounds targeting a broad range of biological processes, including cell cycle progression, genome stability, metabolism, and vesicular transport. Genes and compounds clustered functionally, with inhibitors of related pathways interacting similarly with E3s/DUBs. Some genes, such as FBXW7, showed interactions with many of the compounds. Others, such as RNF25 and FBXO42, showed interactions primarily with a single compound (methyl methanesulfonate for RNF25) or a set of related compounds (the mitotic cluster for FBXO42). Mutation of several E3s with sensitivity to mitotic inhibitors led to increased aberrant mitoses, suggesting a role for these genes in cell cycle regulation. Our comprehensive CRISPR-Cas9 screen uncovered 466 gene-compound interactions covering 25% of the interrogated E3s/DUBs.

Expansion sequencing: Spatially precise in situ transcriptomics in intact biological systems

Methods for highly multiplexed RNA imaging are limited in spatial resolution and thus in their ability to localize transcripts to nanoscale and subcellular compartments. We adapt expansion microscopy, which physically expands biological specimens, for long-read untargeted and targeted in situ RNA sequencing. We applied untargeted expansion sequencing (ExSeq) to the mouse brain, which yielded the readout of thousands of genes, including splice variants. Targeted ExSeq yielded nanoscale-resolution maps of RNAs throughout dendrites and spines in the neurons of the mouse hippocampus, revealing patterns across multiple cell types, layer-specific cell types across the mouse visual cortex, and the organization and position-dependent states of tumor and immune cells in a human metastatic breast cancer biopsy. Thus, ExSeq enables highly multiplexed mapping of RNAs from nanoscale to system scale.

HyPR-seq: Single-cell quantification of chosen RNAs via hybridization and sequencing of DNA probes

Single-cell quantification of RNAs is important for understanding cellular heterogeneity and gene regulation, yet current approaches suffer from low sensitivity for individual transcripts, limiting their utility for many applications. Here we present Hybridization of Probes to RNA for sequencing (HyPR-seq), a method to sensitively quantify the expression of hundreds of chosen genes in single cells. HyPR-seq involves hybridizing DNA probes to RNA, distributing cells into nanoliter droplets, amplifying the probes with PCR, and sequencing the amplicons to quantify the expression of chosen genes. HyPR-seq achieves high sensitivity for individual transcripts, detects nonpolyadenylated and low-abundance transcripts, and can profile more than 100,000 single cells. We demonstrate how HyPR-seq can profile the effects of CRISPR perturbations in pooled screens, detect time-resolved changes in gene expression via measurements of gene introns, and detect rare transcripts and quantify cell-type frequencies in tissue using low-abundance marker genes. By directing sequencing power to genes of interest and sensitively quantifying individual transcripts, HyPR-seq reduces costs by up to 100-fold compared to whole-transcriptome single-cell RNA-sequencing, making HyPR-seq a powerful method for targeted RNA profiling in single cells.

RNA structure-wide discovery of functional interactions with multiplexed RNA motif library

Biochemical assays and computational analyses have discovered RNA structures throughout various transcripts. However, the roles of these structures are mostly unknown. Here we develop folded RNA element profiling with structure library (FOREST), a multiplexed affinity assay system to identify functional interactions from transcriptome-wide RNA structure datasets. We generate an RNA structure library by extracting validated or predicted RNA motifs from gene-annotated RNA regions. The RNA structure library with an affinity enrichment assay allows for the comprehensive identification of target-binding RNA sequences and structures in a high-throughput manner. As a proof-of-concept, FOREST discovers multiple RNA-protein interaction networks with quantitative scores, including translational regulatory elements that function in living cells. Moreover, FOREST reveals different binding landscapes of RNA G-quadruplex (rG4) structures-binding proteins and discovers rG4 structures in the terminal loops of precursor microRNAs. Overall, FOREST serves as a versatile platform to investigate RNA structure-function relationships on a large scale.

Structured RNA motifs can be obtained by structure probing, duplex capture, and motif prediction. Here the authors develop a multiplexed affinity assay system to identify functional protein interactors from an RNA structure library with validated or predicted RNA motifs.

Antigen discovery tools for adaptive immune receptor repertoire research

The adaptive immune system has evolved to recognize with incredible precision a large diversity of molecules. Innovations in high-throughput sequencing and bioinformatics have accelerated large-scale immune repertoire analyses and given us important insights into the behavior of the adaptive immune system. However, establishing a connection between receptor sequence and its antigen-specificity remains a challenge despite its central role in determining T and B cell fate. We discuss recent large-scale antigen discovery technologies which can be combined with adaptive immune receptor repertoire (AIRR) studies. We highlight important discoveries made using repertoire analyses in the field of host-microbe interactions.

Ribo-Pop: simple, cost-effective, and widely applicable ribosomal RNA depletion

The measurement of RNA abundance derived from massively parallel sequencing experiments is an essential technique. Methods that reduce ribosomal RNA levels are usually required prior to sequencing library construction because ribosomal RNA typically comprises the vast majority of a total RNA sample. For some experiments, ribosomal RNA depletion is favored over poly(A) selection because it offers a more inclusive representation of the transcriptome. However, methods to deplete ribosomal RNA are generally proprietary, complex, inefficient, applicable to only specific species, or compatible with only a narrow range of RNA input levels. Here, we describe Ribo-Pop (ribosomal RNA depletion for popular use), a simple workflow and antisense oligo design strategy that we demonstrate works over a wide input range and can be easily adapted to any organism with a sequenced genome. We provide a computational pipeline for probe selection, a streamlined 20-min protocol, and ready-to-use oligo sequences for several organisms. We anticipate that our simple and generalizable "open source" design strategy would enable virtually any laboratory to pursue full transcriptome sequencing in their organism of interest with minimal time and resource investment.

Highly Multiplexed Single-Cell In Situ RNA and DNA Analysis by Consecutive Hybridization

The ability to comprehensively profile nucleic acids in individual cells in their natural spatial contexts is essential to advance our understanding of biology and medicine. Here, we report a novel method for spatial transcriptomics and genomics analysis. In this method, every nucleic acid molecule is detected as a fluorescent spot at its natural cellular location throughout the cycles of consecutive fluorescence in situ hybridization (C-FISH). In each C-FISH cycle, fluorescent oligonucleotide probes hybridize to the probes applied in the previous cycle, and also introduce the binding sites for the next cycle probes. With reiterative cycles of hybridization, imaging and photobleaching, the identities of the varied nucleic acids are determined by their unique color sequences. To demonstrate the feasibility of this method, we show that transcripts or genomic loci in single cells can be unambiguously quantified with 2 fluorophores and 16 C-FISH cycles or with 3 fluorophores and 9 C-FISH cycles. Without any error correction, the error rates obtained using the raw data are close to zero. These results indicate that C-FISH potentially enables tens of thousands (2¹⁶ = 65,536 or 3⁹ = 19,683) of different transcripts or genomic loci to be precisely profiled in individual cells in situ.

T cell antigen discovery

T cells respond to threats in an antigen-specific manner using T cell receptors (TCRs) that recognize short peptide antigens presented on major histocompatibility complex (MHC) proteins. The TCR-peptide-MHC interaction mediated between a T cell and its target cell dictates its function and thereby influences its role in disease. A lack of approaches for antigen discovery has limited the fundamental understanding of the antigenic landscape of the overall T cell response. Recent advances in high-throughput sequencing, mass cytometry, microfluidics and computational biology have led to a surge in approaches to address the challenge of T cell antigen discovery. Here, we summarize the scope of this challenge, discuss in depth the recent exciting work and highlight the outstanding questions and remaining technical hurdles in this field.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

Targeted insertional mutagenesis libraries for deep domain insertion profiling

Domain recombination is a key principle in protein evolution and protein engineering, but inserting a donor domain into every position of a target protein is not easily experimentally accessible. Most contemporary domain insertion profiling approaches rely on DNA transposons, which are constrained by sequence bias. Here, we establish Saturated Programmable Insertion Engineering (SPINE), an unbiased, comprehensive, and targeted domain insertion library generation technique using oligo library synthesis and multi-step Golden Gate cloning. Through benchmarking to MuA transposon-mediated library generation on four ion channel genes, we demonstrate that SPINE-generated libraries are enriched for in-frame insertions, have drastically reduced sequence bias as well as near-complete and highly-redundant coverage. Unlike transposon-mediated domain insertion that was severely biased and sparse for some genes, SPINE generated high-quality libraries for all genes tested. Using the Inward Rectifier K⁺ channel Kir2.1, we validate the practical utility of SPINE by constructing and comparing domain insertion permissibility maps. SPINE is the first technology to enable saturated domain insertion profiling. SPINE could help explore the relationship between domain insertions and protein function, and how this relationship is shaped by evolutionary forces and can be engineered for biomedical applications.

The COMET toolkit for composing customizable genetic programs in mammalian cells

Engineering mammalian cells to carry out sophisticated and customizable genetic programs requires a toolkit of multiple orthogonal and well-characterized transcription factors (TFs). To address this need, we develop the COmposable Mammalian Elements of Transcription (COMET)-an ensemble of TFs and promoters that enable the design and tuning of gene expression to an extent not, to the best of our knowledge, previously possible. COMET currently comprises 44 activating and 12 inhibitory zinc-finger TFs and 83 cognate promoters, combined in a framework that readily accommodates new parts. This system can tune gene expression over three orders of magnitude, provides chemically inducible control of TF activity, and enables single-layer Boolean logic. We also develop a mathematical model that provides mechanistic insights into COMET performance characteristics. Altogether, COMET enables the design and construction of customizable genetic programs in mammalian cells.