Marker-free characterization of full-length transcriptomes of single live circulating tumor cells.

The identification and characterization of circulating tumor cells (CTCs) are important for gaining insights into the biology of metastatic cancers, monitoring disease progression, and medical management of the disease. The limiting factor in the enrichment of purified CTC populations is their sparse availability, heterogeneity, and altered phenotypes relative to the primary tumor. Intensive research both at the technical and molecular fronts led to the development of assays that ease CTC detection and identification from peripheral blood. Most CTC detection methods based on single-cell RNA sequencing (scRNA-seq) use a mix of size selection, marker-based white blood cell (WBC) depletion, and antibodies targeting tumor-associated antigens. However, the majority of these methods either miss out on atypical CTCs or suffer from WBC contamination. We present unCTC, an R package for unbiased identification and characterization of CTCs from single-cell transcriptomic data. unCTC features many standard and novel computational and statistical modules for various analyses. These include a novel method of scRNA-seq clustering, named deep dictionary learning using k-means clustering cost (DDLK), expression-based copy number variation (CNV) inference, and combinatorial, marker-based verification of the malignant phenotypes. DDLK enables robust segregation of CTCs and WBCs in the pathway space, as opposed to the gene expression space. We validated the utility of unCTC on scRNA-seq profiles of breast CTCs from six patients, captured and profiled using an integrated ClearCell FX and Polaris workflow that works by the principles of size-based separation of CTCs and marker-based WBC depletion.

Single-cell genome sequencing of human neurons identifies somatic point mutation and indel enrichment in regulatory elements.

Accurate somatic mutation detection from single-cell DNA sequencing is challenging due to amplification-related artifacts. To reduce this artifact burden, an improved amplification technique, primary template-directed amplification (PTA), was recently introduced. We analyzed whole-genome sequencing data from 52 PTA-amplified single neurons using SCAN2, a new genotyper we developed to leverage mutation signatures and allele balance in identifying somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) in PTA data. Our analysis confirms an increase in nonclonal somatic mutation in single neurons with age, but revises the estimated rate of this accumulation to 16 SNVs per year. We also identify artifacts in other amplification methods. Most importantly, we show that somatic indels increase by at least three per year per neuron and are enriched in functional regions of the genome such as enhancers and promoters. Our data suggest that indels in gene-regulatory elements have a considerable effect on genome integrity in human neurons.

Online monitoring and control of upstream cell culture process using 1D and 2D-LC with SegFlow interface.

The biopharmaceutical industry is transitioning from currently deployed batch-mode bioprocessing to a highly efficient and agile next-generation bioprocessing with the adaptation of continuous bioprocessing, which reduces capital investment and operational costs. Continuous bioprocessing, aligned with FDA's quality-by-design platform, is designed to develop robust processes to deliver safe and effective drugs. With the deployment of knowledge-based operations, product quality can be built into the process to achieve desired critical quality attributes (CQAs) with reduced variability. To facilitate next-generation continuous bioprocessing, it is essential to embrace a fundamental shift-in-paradigm from "quality-by-testing" to "quality-by-design," which requires the deployment of process analytical technologies (PAT). With the adaptation of PAT, a systematic approach of process and product understanding and timely process control are feasible. Deployment of PAT tools for real-time monitoring of CQAs and feedback control is critical for continuous bioprocessing. Given the current deficiency in PAT tools to support continuous bioprocessing, we have integrated Infinity 2D-LC with a post-flow-splitter in conjunction with the SegFlow autosampler to the bioreactors. With this integrated system, we have established a platform for online measurements of titer and CQAs of monoclonal antibodies as well as amino acid analysis of bioreactor cell culture.

Design and Structure-Activity Relationships of Isothiocyanates as Potent and Selective N-Acylethanolamine-Hydrolyzing Acid Amidase Inhibitors.

N-Acylethanolamines are signaling lipid molecules implicated in pathophysiological conditions associated with inflammation and pain. N-Acylethanolamine acid amidase (NAAA) favorably hydrolyzes lipid palmitoylethanolamide, which plays a key role in the regulation of inflammatory and pain processes. The synthesis and structure-activity relationship studies encompassing the isothiocyanate pharmacophore have produced potent low nanomolar inhibitors for hNAAA, while exhibiting high selectivity (>100-fold) against other serine hydrolases and cysteine peptidases. We have followed a target-based structure-activity relationship approach, supported by computational methods and known cocrystals of hNAAA. We have identified systemically active inhibitors with good plasma stability (t1/2 > 2 h) and microsomal stability (t1/2 ∼ 15-30 min) as pharmacological tools to investigate the role of NAAA in inflammation, pain, and drug addiction.

Process analytical technology for on-line monitoring of quality attributes during single-use ultrafiltration/diafiltration.

Process analytical technology (PAT) is a fast-growing field within bioprocessing that enables innovation in biological drug manufacturing. This study demonstrates novel PAT methods for monitoring multiple quality attributes simultaneously during the ultrafiltration and diafiltration (UF/DF) process operation, the final step of monoclonal antibody (mAb) purification. Size exclusion chromatography (SEC) methods were developed to measure excipients arginine, histidine, and high molecular weight (HMW) species using a liquid chromatography (LC) system with autosampler for both on-line and at-line PAT modes. The methods were applied in UF/DF studies for the comparison of single-use tangential flow filtration (TFF) cassettes to standard reusable cassettes to achieve very high concentration mAb drug substance (DS) in the order of 100-200 g/L. These case studies demonstrated that single-use TFF cassettes are a functionally equivalent, low-cost alternative to standard reusable cassettes, and that the on-line PAT measurement of purity and excipient concentration was comparable to orthogonal offline methods. These PAT applications using an on-line LC system equipped with onboard sample dilution can become a platform system for monitoring of multiple attributes over a wide dynamic range, a potentially valuable tool for biological drug development and manufacturing.

Erratum: Iyer, A., et al. Integrative Analysis and Machine Learning Based Characterization of Single Circulating Tumor Cells. J. Clin. Med. 2020, 9, 1206.

In the published manuscript [...].

Online/at-line measurement, analysis and control of product titer and critical product quality attributes (CQAs) during process development.

Process analytical technology (PAT) has been defined by the Food and Drug Administration as a system for designing, analyzing, and controlling manufacturing through timely measurements to ensure final product quality. Based on quality-by-design (QbD) principles, real-time or near-real-time data monitoring is essential for timely control of critical quality attributes (CQAs) to keep the process in a state of control. To facilitate next-generation continuous bioprocessing, deployment of PAT tools for real-time monitoring is integral for process understanding and control. Real-time monitoring and control of CQAs are essential to keep the process within the design space and align with the guiding principles of QbD. The contents of this manuscript are pertinent to the online/at-line monitoring of upstream titer and downstream product quality with timely process control. We demonstrated that an ultra-performance liquid chromatography (UPLC) system interfaced with a UPLC-process sample manager (UPLC-PSM) can be utilized to measure titer and CQAs directly from bioreactors and downstream unit operations, respectively. We established online titer measurements from fed-batch and perfusion-based alternating tangential flow bioreactors as well as product quality assessments of downstream operations for real-time peak collection. This integrated, fully automated system for online data monitoring with feedback control is designed to achieve desired product quality.

Integrative Analysis and Machine Learning based Characterization of Single Circulating Tumor Cells.

We collated publicly available single-cell expression profiles of circulating tumor cells (CTCs) and showed that CTCs across cancers lie on a near-perfect continuum of epithelial to mesenchymal (EMT) transition. Integrative analysis of CTC transcriptomes also highlighted the inverse gene expression pattern between PD-L1 and MHC, which is implicated in cancer immunotherapy. We used the CTCs expression profiles in tandem with publicly available peripheral blood mononuclear cell (PBMC) transcriptomes to train a classifier that accurately recognizes CTCs of diverse phenotype. Further, we used this classifier to validate circulating breast tumor cells captured using a newly developed microfluidic system for label-free enrichment of CTCs.

Design and synthesis of cyanamides as potent and selective N-acylethanolamine acid amidase inhibitors.

N-acylethanolamine acid amidase (NAAA) inhibition represents an exciting novel approach to treat inflammation and pain. NAAA is a cysteine amidase which preferentially hydrolyzes the endogenous biolipids palmitoylethanolamide (PEA) and oleoylethanolamide (OEA). PEA is an endogenous agonist of the nuclear peroxisome proliferator-activated receptor-α (PPAR-α), which is a key regulator of inflammation and pain. Thus, blocking the degradation of PEA with NAAA inhibitors results in augmentation of the PEA/PPAR-α signaling pathway and regulation of inflammatory and pain processes. We have prepared a new series of NAAA inhibitors exploring the azetidine-nitrile (cyanamide) pharmacophore that led to the discovery of highly potent and selective compounds. Key analogs demonstrated single-digit nanomolar potency for hNAAA and showed >100-fold selectivity against serine hydrolases FAAH, MGL and ABHD6, and cysteine protease cathepsin K. Additionally, we have identified potent and selective dual NAAA-FAAH inhibitors to investigate a potential synergism between two distinct anti-inflammatory molecular pathways, the PEA/PPAR-α anti-inflammatory signaling pathway,1-4 and the cannabinoid receptors CB1 and CB2 pathways which are known for their antiinflammatory and antinociceptive properties.5-8 Our ligand design strategy followed a traditional structure-activity relationship (SAR) approach and was supported by molecular modeling studies of reported X-ray structures of hNAAA. Several inhibitors were evaluated in stability assays and demonstrated very good plasma stability (t1/2 > 2 h; human and rodents). The disclosed cyanamides represent promising new pharmacological tools to investigate the potential role of NAAA inhibitors and dual NAAA-FAAH inhibitors as therapeutic agents for the treatment of inflammation and pain.

Multimodal Single-Cell Analysis Reveals Physiological Maturation in the Developing Human Neocortex.

In the developing human neocortex, progenitor cells generate diverse cell types prenatally. Progenitor cells and newborn neurons respond to signaling cues, including neurotransmitters. While single-cell RNA sequencing has revealed cellular diversity, physiological heterogeneity has yet to be mapped onto these developing and diverse cell types. By combining measurements of intracellular Ca2+ elevations in response to neurotransmitter receptor agonists and RNA sequencing of the same single cells, we show that Ca2+ responses are cell-type-specific and change dynamically with lineage progression. Physiological response properties predict molecular cell identity and additionally reveal diversity not captured by single-cell transcriptomics. We find that the serotonin receptor HTR2A selectively activates radial glia cells in the developing human, but not mouse, neocortex, and inhibiting HTR2A receptors in human radial glia disrupts the radial glial scaffold. We show highly specific neurotransmitter signaling during neurogenesis in the developing human neocortex and highlight evolutionarily divergent mechanisms of physiological signaling.

Establishing Cerebral Organoids as Models of Human-Specific Brain Evolution.

Direct comparisons of human and non-human primate brains can reveal molecular pathways underlying remarkable specializations of the human brain. However, chimpanzee tissue is inaccessible during neocortical neurogenesis when differences in brain size first appear. To identify human-specific features of cortical development, we leveraged recent innovations that permit generating pluripotent stem cell-derived cerebral organoids from chimpanzee. Despite metabolic differences, organoid models preserve gene regulatory networks related to primary cell types and developmental processes. We further identified 261 differentially expressed genes in human compared to both chimpanzee organoids and macaque cortex, enriched for recent gene duplications, and including multiple regulators of PI3K-AKT-mTOR signaling. We observed increased activation of this pathway in human radial glia, dependent on two receptors upregulated specifically in human: INSR and ITGB8. Our findings establish a platform for systematic analysis of molecular changes contributing to human brain development and evolution.

C1 CAGE detects transcription start sites and enhancer activity at single-cell resolution.

Single-cell transcriptomic profiling is a powerful tool to explore cellular heterogeneity. However, most of these methods focus on the 3'-end of polyadenylated transcripts and provide only a partial view of the transcriptome. We introduce C1 CAGE, a method for the detection of transcript 5'-ends with an original sample multiplexing strategy in the C1TM microfluidic system. We first quantifiy the performance of C1 CAGE and find it as accurate and sensitive as other methods in the C1 system. We then use it to profile promoter and enhancer activities in the cellular response to TGF-ß of lung cancer cells and discover subpopulations of cells differing in their response. We also describe enhancer RNA dynamics revealing transcriptional bursts in subsets of cells with transcripts arising from either strand in a mutually exclusive manner, validated using single molecule fluorescence in situ hybridization.

Regulation of cell-type-specific transcriptomes by microRNA networks during human brain development.

MicroRNAs (miRNAs) regulate many cellular events during brain development by interacting with hundreds of mRNA transcripts. However, miRNAs operate nonuniformly upon the transcriptional profile with an as yet unknown logic. Shortcomings in defining miRNA-mRNA networks include limited knowledge of in vivo miRNA targets and their abundance in single cells. By combining multiple complementary approaches, high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation with an antibody to AGO2 (AGO2-HITS-CLIP), single-cell profiling and computational analyses using bipartite and coexpression networks, we show that miRNA-mRNA interactions operate as functional modules that often correspond to cell-type identities and undergo dynamic transitions during brain development. These networks are highly dynamic during development and over the course of evolution. One such interaction is between radial-glia-enriched ORC4 and miR-2115, a great-ape-specific miRNA, which appears to control radial glia proliferation rates during human brain development.

Vitamin B12 association with mAbs: Mechanism and potential mitigation strategies.

Process control for manufacturing biologics is critical for ensuring product quality, safety, and lot to lot consistency of therapeutic proteins. In this study, we investigated the root cause of the pink coloration observed for various in-process pools and drug substances in the antibody manufacturing process. Vitamin B12 is covalently bound to mAbs via a cobalt-sulfur coordinate bond via the cysteine residues. The vitamin B12 was identified to attach to an IgG4 molecule at cysteine residues on light chain (Cys-214), and heavy chain (Cys-134, Cys-321, Cys-367, and Cys-425). Prior to attachment to mAbs, the vitamin B12 needs to be in its active form of hydroxocobalamin. During culture media preparation, storage and cell culture processing, cyanocobalamin, the chemical form of vitamin B12 added to media, is converted to hydroxocobalamin by white fluorescence light (about 50% degradation in 11-14 days at room temperature and with room light intensity about 500-1,000 lux) and by short-wavelength visible light (400-550 nm). However, cyanocobalamin is stable under red light (wavelength >600 nm) exposure and does not convert to hydroxocobalamin. Our findings suggests that the intensity of pink color depends on concentrations of both free sulfhydryl groups on reduced mAb and hydroxocobalamin, the active form of vitamin B12 . Both reactants are necessary and neither one of them is sufficient to generate pink color, therefore process control strategy can consider limiting either one or both factors. A process control strategy to install red light (wavelength >600 nm) in culture media preparation, storage and culture processing areas is proposed to provide safe light for biologics and to prevent light-induced color variations in final products.

Secretion, isotopic labeling and deglycosylation of N-acylethanolamine acid amidase for biophysical studies.

N-acylethanolamine acid amidase (NAAA) is an N-terminal nucleophile (Ntn) enzyme with a catalytic cysteine residue that has highest activity at acidic pH. The most prominent substrate hydrolyzed is palmitoylethanolamine (PEA), which regulates inflammation. Inhibitors of NAAA have been shown to increase endogenous levels of PEA, and are of interest as potential treatments for inflammatory disorders and other maladies. Currently, there are no X-ray or NMR structures of NAAA available to inform medicinal chemistry. Additionally, there are a limited number of enzyme structures available that are within the Ntn-hydrolase family, have a catalytic cysteine residue, and have a high sequence homology. For these reasons, we developed expression and purification methods for the production of enzyme samples amenable to structural characterization. Mammalian cells are necessary for post-translational processing, including signal sequence cleavage and glycosylation, that are required for a correctly folded zymogen before conversion to active, and mature enzyme. We have identified an expression construct, mammalian cell line, specific media and additives to express and secrete hNAAA zymogen and we further optimized propagation conditions and show this secretion method is suitable for isotopic labeling of the protein. We refined purification methods to achieve a high degree of protein purity potentially suited to crystallography. Glycosylated proteins can present challenges to biophysical methods. Therefore we deglycosylate the enzyme and show that the activity of the mature enzyme is not affected by deglycosylation.

Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex.

Systematic analyses of spatiotemporal gene expression trajectories during organogenesis have been challenging because diverse cell types at different stages of maturation and differentiation coexist in the emerging tissues. We identified discrete cell types as well as temporally and spatially restricted trajectories of radial glia maturation and neurogenesis in developing human telencephalon. These lineage-specific trajectories reveal the expression of neurogenic transcription factors in early radial glia and enriched activation of mammalian target of rapamycin signaling in outer radial glia. Across cortical areas, modest transcriptional differences among radial glia cascade into robust typological distinctions among maturing neurons. Together, our results support a mixed model of topographical, typological, and temporal hierarchies governing cell-type diversity in the developing human telencephalon, including distinct excitatory lineages emerging in rostral and caudal cerebral cortex.

Corrigendum: Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

[This corrects the article on p. 70 in vol. 4, PMID: 27709111.].

The nature and nurture of cell heterogeneity: accounting for macrophage gene-environment interactions with single-cell RNA-Seq.

BACKGROUND: Single-cell RNA-Seq can be a valuable and unbiased tool to dissect cellular heterogeneity, despite the transcriptome's limitations in describing higher functional phenotypes and protein events. Perhaps the most important shortfall with transcriptomic 'snapshots' of cell populations is that they risk being descriptive, only cataloging heterogeneity at one point in time, and without microenvironmental context. Studying the genetic ('nature') and environmental ('nurture') modifiers of heterogeneity, and how cell population dynamics unfold over time in response to these modifiers is key when studying highly plastic cells such as macrophages. RESULTS: We introduce the programmable Polaris™ microfluidic lab-on-chip for single-cell sequencing, which performs live-cell imaging while controlling for the culture microenvironment of each cell. Using gene-edited macrophages we demonstrate how previously unappreciated knockout effects of SAMHD1, such as an altered oxidative stress response, have a large paracrine signaling component. Furthermore, we demonstrate single-cell pathway enrichments for cell cycle arrest and APOBEC3G degradation, both associated with the oxidative stress response and altered proteostasis. Interestingly, SAMHD1 and APOBEC3G are both HIV-1 inhibitors ('restriction factors'), with no known co-regulation. CONCLUSION: As single-cell methods continue to mature, so will the ability to move beyond simple 'snapshots' of cell populations towards studying the determinants of population dynamics. By combining single-cell culture, live-cell imaging, and single-cell sequencing, we have demonstrated the ability to study cell phenotypes and microenvironmental influences. It's these microenvironmental components - ignored by standard single-cell workflows - that likely determine how macrophages, for example, react to inflammation and form treatment resistant HIV reservoirs.

Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

The study of single cells has evolved over the past several years to include expression and genomic analysis of an increasing number of single cells. Several studies have demonstrated wide spread variation and heterogeneity within cell populations of similar phenotype. While the characterization of these populations will likely set the foundation for our understanding of genomic- and expression-based diversity, it will not be able to link the functional differences of a single cell to its underlying genomic structure and activity. Currently, it is difficult to perturb single cells in a controlled environment, monitor and measure the response due to perturbation, and link these response measurements to downstream genomic and transcriptomic analysis. In order to address this challenge, we developed a platform to integrate and miniaturize many of the experimental steps required to study single-cell function. The heart of this platform is an elastomer-based integrated fluidic circuit that uses fluidic logic to select and sequester specific single cells based on a phenotypic trait for downstream experimentation. Experiments with sequestered cells that have been performed include on-chip culture, exposure to various stimulants, and post-exposure image-based response analysis, followed by preparation of the mRNA transcriptome for massively parallel sequencing analysis. The flexible system embodies experimental design and execution that enable routine functional studies of single cells.

Molecular identity of human outer radial glia during cortical development.

Radial glia, the neural stem cells of the neocortex, are located in two niches: the ventricular zone and outer subventricular zone. Although outer subventricular zone radial glia may generate the majority of human cortical neurons, their molecular features remain elusive. By analyzing gene expression across single cells, we find that outer radial glia preferentially express genes related to extracellular matrix formation, migration, and stemness, including TNC, PTPRZ1, FAM107A, HOPX, and LIFR. Using dynamic imaging, immunostaining, and clonal analysis, we relate these molecular features to distinctive behaviors of outer radial glia, demonstrate the necessity of STAT3 signaling for their cell cycle progression, and establish their extensive proliferative potential. These results suggest that outer radial glia directly support the subventricular niche through local production of growth factors, potentiation of growth factor signals by extracellular matrix proteins, and activation of self-renewal pathways, thereby enabling the developmental and evolutionary expansion of the human neocortex.