Inferring a spatial code of cell-cell interactions across a whole animal body

Cell-cell interactions shape cellular function and ultimately organismal phenotype. Interacting cells can sense their mutual distance using combinations of ligand-receptor pairs, suggesting the existence of a spatial code, i.e., signals encoding spatial properties of cellular organization. However, this code driving and sustaining the spatial organization of cells remains to be elucidated. Here we present a computational framework to infer the spatial code underlying cell-cell interactions from the transcriptomes of the cell types across the whole body of a multicellular organism. As core of this framework, we introduce our tool cell2cell, which uses the coexpression of ligand-receptor pairs to compute the potential for intercellular interactions, and we test it across the Caenorhabditis elegans' body. Leveraging a 3D atlas of C. elegans' cells, we also implement a genetic algorithm to identify the ligand-receptor pairs most informative of the spatial organization of cells across the whole body. Validating the spatial code extracted with this strategy, the resulting intercellular distances are negatively correlated with the inferred cell-cell interactions. Furthermore, for selected cell-cell and ligand-receptor pairs, we experimentally confirm the communicatory behavior inferred with cell2cell and the genetic algorithm. Thus, our framework helps identify a code that predicts the spatial organization of cells across a whole-animal body.

Spatiotemporal dynamics of human high gamma discriminate naturalistic behavioral states

In analyzing the neural correlates of naturalistic and unstructured behaviors, features of neural activity that are ignored in a trial-based experimental paradigm can be more fully studied and investigated. Here, we analyze neural activity from two patients using electrocorticography (ECoG) and stereo-electroencephalography (sEEG) recordings, and reveal that multiple neural signal characteristics exist that discriminate between unstructured and naturalistic behavioral states such as “engaging in dialogue” and “using electronics”. Using the high gamma amplitude as an estimate of neuronal firing rate, we demonstrate that behavioral states in a naturalistic setting are discriminable based on long-term mean shifts, variance shifts, and differences in the specific neural activity’s covariance structure. Both the rapid and slow changes in high gamma band activity separate unstructured behavioral states. We also use Gaussian process factor analysis (GPFA) to show the existence of salient spatiotemporal features with variable smoothness in time. Further, we demonstrate that both temporally smooth and stochastic spatiotemporal activity can be used to differentiate unstructured behavioral states. This is the first attempt to elucidate how different neural signal features contain information about behavioral states collected outside the conventional experimental paradigm.

Immune response to intravenous immunoglobulin in patients with Kawasaki disease and MIS-C

BACKGROUNDMultisystem inflammatory syndrome in children (MIS-C) is a rare but potentially severe illness that follows exposure to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Kawasaki disease (KD) shares several clinical features with MIS-C, which prompted the use of intravenous immunoglobulin (IVIG), a mainstay therapy for KD. Both diseases share a robust activation of the innate immune system, including the IL-1 signaling pathway, and IL-1 blockade has been used for the treatment of both MIS-C and KD. The mechanism of action of IVIG in these 2 diseases and the cellular source of IL-1β have not been defined.METHODSThe effects of IVIG on peripheral blood leukocyte populations from patients with MIS-C and KD were examined using flow cytometry and mass cytometry (CyTOF) and live-cell imaging.RESULTSCirculating neutrophils were highly activated in patients with KD and MIS-C and were a major source of IL-1β. Following IVIG treatment, activated IL-1β+ neutrophils were reduced in the circulation. In vitro, IVIG was a potent activator of neutrophil cell death via PI3K and NADPH oxidase, but independently of caspase activation.CONCLUSIONSActivated neutrophils expressing IL-1β can be targeted by IVIG, supporting its use in both KD and MIS-C to ameliorate inflammation.FUNDINGPatient Centered Outcomes Research Institute; NIH; American Asthma Foundation; American Heart Association; Novo Nordisk Foundation; NIGMS; American Academy of Allergy, Asthma and Immunology Foundation.

A Chinese hamster transcription start site atlas that enables targeted editing of CHO cells

Chinese hamster ovary (CHO) cells are widely used for producing biopharmaceuticals, and engineering gene expression in CHO is key to improving drug quality and affordability. However, engineering gene expression or activating silent genes requires accurate annotation of the underlying regulatory elements and transcription start sites (TSSs). Unfortunately, most TSSs in the published Chinese hamster genome sequence were computationally predicted and are frequently inaccurate. Here, we use nascent transcription start site sequencing methods to revise TSS annotations for 15 308 Chinese hamster genes and 3034 non-coding RNAs based on experimental data from CHO-K1 cells and 10 hamster tissues. We further capture tens of thousands of putative transcribed enhancer regions with this method. Our revised TSSs improves upon the RefSeq annotation by revealing core sequence features of gene regulation such as the TATA box and the Initiator and, as exemplified by targeting the glycosyltransferase gene Mgat3, facilitate activating silent genes by CRISPRa. Together, we envision our revised annotation and data will provide a rich resource for the CHO community, improve genome engineering efforts and aid comparative and evolutionary studies.

Multiple freeze-thaw cycles lead to a loss of consistency in poly(A)-enriched RNA sequencing

Background

Both RNA-Seq and sample freeze-thaw are ubiquitous. However, knowledge about the impact of freeze-thaw on downstream analyses is limited. The lack of common quality metrics that are sufficiently sensitive to freeze-thaw and RNA degradation, e.g. the RNA Integrity Score, makes such assessments challenging.

Results

Here we quantify the impact of repeated freeze-thaw cycles on the reliability of RNA-Seq by examining poly(A)-enriched and ribosomal RNA depleted RNA-seq from frozen leukocytes drawn from a toddler Autism cohort. To do so, we estimate the relative noise, or percentage of random counts, separating technical replicates. Using this approach we measured noise associated with RIN and freeze-thaw cycles. As expected, RIN does not fully capture sample degradation due to freeze-thaw. We further examined differential expression results and found that three freeze-thaws should extinguish the differential expression reproducibility of similar experiments. Freeze-thaw also resulted in a 3′ shift in the read coverage distribution along the gene body of poly(A)-enriched samples compared to ribosomal RNA depleted samples, suggesting that library preparation may exacerbate freeze-thaw-induced sample degradation.

Conclusion

The use of poly(A)-enrichment for RNA sequencing is pervasive in library preparation of frozen tissue, and thus, it is important during experimental design and data analysis to consider the impact of repeated freeze-thaw cycles on reproducibility.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07381-z.

Awakening dormant glycosyltransferases in CHO cells with CRISPRa

Chinese hamster ovary (CHO) cells are the preferred workhorse for the biopharmaceutical industry, and CRISPR/Cas9 has proven powerful for generating targeted gene perturbations in CHO cells. Here, we expand the CRISPR engineering toolbox with CRISPR activation (CRISPRa) to increase transcription of endogenous genes. We successfully increased transcription of Mgat3 and St6gal1, and verified their activity on a functional level by subsequently detecting that the appropriate glycan structures were produced. This study demonstrates that CRISPRa can make targeted alterations of CHO cells for desired phenotypes.

Reduced apoptosis in Chinese hamster ovary cells via optimized CRISPR interference

Chinese hamster ovary (CHO) cells are widely used for biopharmaceutical protein production. One challenge limiting CHO cell productivity is apoptosis stemming from cellular stress during protein production. Here we applied CRISPR interference (CRISPRi) to downregulate the endogenous expression of apoptotic genes Bak, Bax, and Casp3 in CHO cells. In addition to reduced apoptosis, mitochondrial membrane integrity was improved and the caspase activity was reduced. Moreover, we optimized the CRISPRi system to enhance the gene repression efficiency in CHO cells by testing different repressor fusion types. An improved Cas9 repressor has been identified by applying C-terminal fusion of a bipartite repressor domain, KRAB-MeCP2, to nuclease-deficient Cas9. These results collectively demonstrate that CHO cells can be rescued from cell apoptosis by targeted gene repression using the CRISPRi system.

Replay of large-scale spatio-temporal patterns from waking during subsequent NREM sleep in human cortex

Animal studies support the hypothesis that in slow-wave sleep, replay of waking neocortical activity under hippocampal guidance leads to memory consolidation. However, no intracranial electrophysiological evidence for replay exists in humans. We identified consistent sequences of population firing peaks across widespread cortical regions during complete waking periods. The occurrence of these “Motifs” were compared between sleeps preceding the waking period (“Sleep-Pre”) when the Motifs were identified, and those following (“Sleep-Post”). In all subjects, the majority of waking Motifs (most of which were novel) had more matches in Sleep-Post than in Sleep-Pre. In rodents, hippocampal replay occurs during local sharp-wave ripples, and the associated neocortical replay tends to occur during local sleep spindles and down-to-up transitions. These waves may facilitate consolidation by sequencing cell-firing and encouraging plasticity. Similarly, we found that Motifs were coupled to neocortical spindles, down-to-up transitions, theta bursts, and hippocampal sharp-wave ripples. While Motifs occurring during cognitive task performance were more likely to have more matches in subsequent sleep, our studies provide no direct demonstration that the replay of Motifs contributes to consolidation. Nonetheless, these results confirm a core prediction of the dominant neurobiological theory of human memory consolidation.

The emerging role of systems biology for engineering protein production in CHO cells

To meet the ever-growing demand for effective, safe, and affordable protein therapeutics, decades of intense efforts have aimed to maximize the quantity and quality of recombinant proteins produced in CHO cells. Bioprocessing innovations and cell engineering efforts have improved product titer; however, uncharacterized cellular processes and gene regulatory mechanisms still hinder cell growth, specific productivity, and protein quality. Herein, we summarize recent advances in systems biology and data-driven approaches aiming to unravel how molecular pathways, cellular processes, and extrinsic factors (e.g. media supplementation) influence recombinant protein production. In particular, as the available omics data for CHO cells continue to grow, predictive models and screens will be increasingly used to unravel the biological drivers of protein production, which can be used with emerging genome editing technologies to rationally engineer cells to further control the quantity, quality and affordability of many biologic drugs.

Cellular and neurochemical basis of sleep stages in the thalamocortical network

The link between the combined action of neuromodulators in the brain and global brain states remains a mystery. In this study, using biophysically realistic models of the thalamocortical network, we identified the critical intrinsic and synaptic mechanisms, associated with the putative action of acetylcholine (ACh), GABA and monoamines, which lead to transitions between primary brain vigilance states (waking, non-rapid eye movement sleep [NREM] and REM sleep) within an ultradian cycle. Using ECoG recordings from humans and LFP recordings from cats and mice, we found that during NREM sleep the power of spindle and delta oscillations is negatively correlated in humans and positively correlated in animal recordings. We explained this discrepancy by the differences in the relative level of ACh. Overall, our study revealed the critical intrinsic and synaptic mechanisms through which different neuromodulators acting in combination result in characteristic brain EEG rhythms and transitions between sleep stages.