An Excitatory Circuit in the Perioculomotor Midbrain for Non-REM Sleep Control.

The perioculomotor (pIII) region of the midbrain was postulated as a sleep-regulating center in the 1890s but largely neglected in subsequent studies. Using activity-dependent labeling and gene expression profiling, we identified pIII neurons that promote non-rapid eye movement (NREM) sleep. Optrode recording showed that pIII glutamatergic neurons expressing calcitonin gene-related peptide alpha (CALCA) are NREM-sleep active; optogenetic and chemogenetic activation/inactivation showed that they strongly promote NREM sleep. Within the pIII region, CALCA neurons form reciprocal connections with another population of glutamatergic neurons that express the peptide cholecystokinin (CCK). Activation of CCK neurons also promoted NREM sleep. Both CALCA and CCK neurons project rostrally to the preoptic hypothalamus, whereas CALCA neurons also project caudally to the posterior ventromedial medulla. Activation of each projection increased NREM sleep. Together, these findings point to the pIII region as an excitatory sleep center where different subsets of glutamatergic neurons promote NREM sleep through both local reciprocal connections and long-range projections.

Restoration of DNA repair mitigates genome instability and increases productivity of Chinese hamster ovary cells.

Chinese hamster ovary (CHO) cells are the primary host for manufacturing of therapeutic proteins. However, productivity loss is a major problem and is associated with genome instability, as chromosomal aberrations reduce transgene copy number and decrease protein expression. We analyzed whole-genome sequencing data from 11 CHO cell lines and found deleterious single-nucleotide variants in DNA repair genes. Comparison with primary Chinese hamster cells confirmed DNA repair to be compromised in CHO. Correction of key DNA repair genes by single-nucleotide variant reversal or expression of intact complementary DNAs successfully improved DNA repair and mitigated karyotypic instability. Moreover, overexpression of intact copies of LIG4 and XRCC6 in a CHO cell line expressing secreted alkaline phosphatase mitigated transgene copy loss and improved protein titer retention. These results show that correction of DNA repair genes yields improvements in genome stability in CHO, and provide new opportunities for cell line development for sustainable protein expression.

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.

Systematically gap-filling the genome-scale metabolic model of CHO cells.

OBJECTIVE: Chinese hamster ovary (CHO) cells are the leading cell factories for producing recombinant proteins in the biopharmaceutical industry. In this regard, constraint-based metabolic models are useful platforms to perform computational analysis of cell metabolism. These models need to be regularly updated in order to include the latest biochemical data of the cells, and to increase their predictive power. Here, we provide an update to iCHO1766, the metabolic model of CHO cells. RESULTS: We expanded the existing model of Chinese hamster metabolism with the help of four gap-filling approaches, leading to the addition of 773 new reactions and 335 new genes. We incorporated these into an updated genome-scale metabolic network model of CHO cells, named iCHO2101. In this updated model, the number of reactions and pathways capable of carrying flux is substantially increased. CONCLUSIONS: The present CHO model is an important step towards more complete metabolic models of CHO cells.

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.

Proteogenomic Annotation of Chinese Hamsters Reveals Extensive Novel Translation Events and Endogenous Retroviral Elements.

A high-quality genome annotation greatly facilitates successful cell line engineering. Standard draft genome annotation pipelines are based largely on de novo gene prediction, homology, and RNA-Seq data. However, draft annotations can suffer from incorrect predictions of translated sequence, inaccurate splice isoforms, and missing genes. Here, we generated a draft annotation for the newly assembled Chinese hamster genome and used RNA-Seq, proteomics, and Ribo-Seq to experimentally annotate the genome. We identified 3529 new proteins compared to the hamster RefSeq protein annotation and 2256 novel translational events (e.g., alternative splices, mutations, and novel splices). Finally, we used this pipeline to identify the source of translated retroviruses contaminating recombinant products from Chinese hamster ovary (CHO) cell lines, including 119 type-C retroviruses, thus enabling future efforts to eliminate retroviruses to reduce the costs incurred with retroviral particle clearance. In summary, the improved annotation provides a more accurate resource for CHO cell line engineering, by facilitating the interpretation of omics data, defining of cellular pathways, and engineering of complex 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.

A reference genome of the Chinese hamster based on a hybrid assembly strategy.

Accurate and complete genome sequences are essential in biotechnology to facilitate genome-based cell engineering efforts. The current genome assemblies for Cricetulus griseus, the Chinese hamster, are fragmented and replete with gap sequences and misassemblies, consistent with most short-read-based assemblies. Here, we completely resequenced C. griseus using single molecule real time sequencing and merged this with Illumina-based assemblies. This generated a more contiguous and complete genome assembly than either technology alone, reducing the number of scaffolds by >28-fold, with 90% of the sequence in the 122 longest scaffolds. Most genes are now found in single scaffolds, including up- and downstream regulatory elements, enabling improved study of noncoding regions. With >95% of the gap sequence filled, important Chinese hamster ovary cell mutations have been detected in draft assembly gaps. This new assembly will be an invaluable resource for continued basic and pharmaceutical research.

Rapid neurogenesis through transcriptional activation in human stem cells.

Advances in cellular reprogramming and stem cell differentiation now enable ex vivo studies of human neuronal differentiation. However, it remains challenging to elucidate the underlying regulatory programs because differentiation protocols are laborious and often result in low neuron yields. Here, we overexpressed two Neurogenin transcription factors in human-induced pluripotent stem cells and obtained neurons with bipolar morphology in 4 days, at greater than 90% purity. The high purity enabled mRNA and microRNA expression profiling during neurogenesis, thus revealing the genetic programs involved in the rapid transition from stem cell to neuron. The resulting cells exhibited transcriptional, morphological and functional signatures of differentiated neurons, with greatest transcriptional similarity to prenatal human brain samples. Our analysis revealed a network of key transcription factors and microRNAs that promoted loss of pluripotency and rapid neurogenesis via progenitor states. Perturbations of key transcription factors affected homogeneity and phenotypic properties of the resulting neurons, suggesting that a systems-level view of the molecular biology of differentiation may guide subsequent manipulation of human stem cells to rapidly obtain diverse neuronal types.

Occurrence of macroplastic debris in the long-term plastic film-mulched agricultural soil: A case study of Northwest China.

Widespread presence of plastic mulch has led to macroplastic (MaP) pollution. While this issue is widely explored in aquatic ecosystems, MaP pollution on land has been neglected. In 2019, we conducted a large-scale survey of MaPs in Northwest China in 0-30 cm soil with long-term mulching. Samples of MaP debris were collected from 67 sites across Gansu, east Qinghai, and north Shannxi Provinces. All visible MaP pieces for each site were separated and weighed. The mass of each MaP piece was calibrated by size measured in digital images. The MaP mass averaged 47.2 kg ha-1, and the number of MaPs averaged 266.2 pieces ha-1. The mass and number of MaPs varied from site to site. The mean size of MaPs was 19.5 cm2 piece-1 or 28.0 mg piece-1. More importantly, the number of small MaP pieces (<5 and 5-20 cm2 piece-1) accounted for 76.7% of the total number of MaPs detected, and small-sized plastic debris (<10 and 10-25 mg piece-1) were detected in 70.1% of the sampling sites. The percentage of small fragments increased before 15-year of mulching and then declined. However, the amount of medium-large debris (20-50 and >50 cm2 piece-1) showed a trend opposite to that of small fragments. The percentage of MaPs was greater in the small size group than in the medium-large size group. The arid to semi-arid area exhibited higher MaP contamination compared with the semi-arid to the semi-humid area. These observations indicate that plastic debris residing in soil tend to be fragmented, making plastic film recovery more challenging and causing severe soil pollution. Further studies are required to regulate plastic mulch methods and explore the degradation process.

Rare germline variants in individuals diagnosed with schizophrenia within multiplex families.

An extensive catalog of common and rare genetic variants contributes to overall risk for schizophrenia and related disorders. As a complement to population genetics efforts, here we present whole genome sequences of multiple affected probands within individual families to search for possible high penetrance driver variants. From a total of 15 families diagnostically evaluated by a single research psychiatrist, we performed whole genome sequencing of a total of 61 affected individuals, called SNPs, indels, and copy number variants, and compared to reference genomes. In fourteen out of fifteen families, the schizophrenia polygenic risk score for each proband was within the control range defined by the Thousand Genomes cohort. In six families, each affected member carried a very rare or private, predicted-damaging, variant in at least one gene. Among these genes, variants in LRP1 and TENM2 suggest these are candidate disease-related genes when taken into context with existing population genetic studies and biological information. Results add to the number of pedigree sequences reported, suggest pathways for the investigation of biological mechanisms, and are consistent with the overall accumulating evidence that very rare damaging variants contribute to the heritability of schizophrenia.

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.

Erratum: An optimized genome-wide, virus-free CRISPR screen for mammalian cells.

[This corrects the article DOI: 10.1016/j.crmeth.2021.100062.].

An optimized genome-wide, virus-free CRISPR screen for mammalian cells.

Pooled CRISPR screens have been widely applied to mammalian and other organisms to elucidate the interplay between genes and phenotypes of interest. The most popular method for delivering the CRISPR components into mammalian cells is lentivirus based. However, because lentivirus is not always an option, virus-free protocols are starting to emerge. Here, we demonstrate an improved virus-free, genome-wide CRISPR screening platform for Chinese hamster ovary cells with 75,488 gRNAs targeting 15,028 genes. Each gRNA expression cassette in the library is precisely integrated into a genomic landing pad, resulting in a very high percentage of single gRNA insertions and minimal clonal variation. Using this platform, we perform a negative selection screen on cell proliferation that identifies 1,980 genes that affect proliferation and a positive selection screen on the toxic endoplasmic reticulum stress inducer, tunicamycin, that identifies 77 gene knockouts that improve survivability.

Model-based assessment of mammalian cell metabolic functionalities using omics data.

Omics experiments are ubiquitous in biological studies, leading to a deluge of data. However, it is still challenging to connect changes in these data to changes in cell functions because of complex interdependencies between genes, proteins, and metabolites. Here, we present a framework allowing researchers to infer how metabolic functions change on the basis of omics data. To enable this, we curated and standardized lists of metabolic tasks that mammalian cells can accomplish. Genome-scale metabolic networks were used to define gene sets associated with each metabolic task. We further developed a framework to overlay omics data on these sets and predict pathway usage for each metabolic task. We demonstrated how this approach can be used to quantify metabolic functions of diverse biological samples from the single cell to whole tissues and organs by using multiple transcriptomic datasets. To facilitate its adoption, we integrated the approach into GenePattern (www.genepattern.org-CellFie).

Genome-scale reconstructions of the mammalian secretory pathway predict metabolic costs and limitations of protein secretion.

In mammalian cells, >25% of synthesized proteins are exported through the secretory pathway. The pathway complexity, however, obfuscates its impact on the secretion of different proteins. Unraveling its impact on diverse proteins is particularly important for biopharmaceutical production. Here we delineate the core secretory pathway functions and integrate them with genome-scale metabolic reconstructions of human, mouse, and Chinese hamster ovary cells. The resulting reconstructions enable the computation of energetic costs and machinery demands of each secreted protein. By integrating additional omics data, we find that highly secretory cells have adapted to reduce expression and secretion of other expensive host cell proteins. Furthermore, we predict metabolic costs and maximum productivities of biotherapeutic proteins and identify protein features that most significantly impact protein secretion. Finally, the model successfully predicts the increase in secretion of a monoclonal antibody after silencing a highly expressed selection marker. This work represents a knowledgebase of the mammalian secretory pathway that serves as a novel tool for systems biotechnology.

Combating viral contaminants in CHO cells by engineering innate immunity.

Viral contamination in biopharmaceutical manufacturing can lead to shortages in the supply of critical therapeutics. To facilitate the protection of bioprocesses, we explored the basis for the susceptibility of CHO cells to RNA virus infection. Upon infection with certain ssRNA and dsRNA viruses, CHO cells fail to generate a significant interferon (IFN) response. Nonetheless, the downstream machinery for generating IFN responses and its antiviral activity is intact in these cells: treatment of cells with exogenously-added type I IFN or poly I:C prior to infection limited the cytopathic effect from Vesicular stomatitis virus (VSV), Encephalomyocarditis virus (EMCV), and Reovirus-3 virus (Reo-3) in a STAT1-dependent manner. To harness the intrinsic antiviral mechanism, we used RNA-Seq to identify two upstream repressors of STAT1: Gfi1 and Trim24. By knocking out these genes, the engineered CHO cells exhibited activation of cellular immune responses and increased resistance to the RNA viruses tested. Thus, omics-guided engineering of mammalian cell culture can be deployed to increase safety in biotherapeutic protein production among many other biomedical applications.

Dual RNA-seq identifies human mucosal immunity protein Mucin-13 as a hallmark of Plasmodium exoerythrocytic infection.

The exoerythrocytic stage of Plasmodium infection is a critical window for prophylactic intervention. Using genome-wide dual RNA sequencing of flow-sorted infected and uninfected hepatoma cells we show that the human mucosal immunity gene, mucin-13 (MUC13), is strongly upregulated during Plasmodium exoerythrocytic hepatic-stage infection. We confirm MUC13 transcript increases in hepatoma cell lines and primary hepatocytes. In immunofluorescence assays, host MUC13 protein expression distinguishes infected cells from adjacent uninfected cells and shows similar colocalization with parasite biomarkers such as UIS4 and HSP70. We further show that localization patterns are species independent, marking both P. berghei and P. vivax infected cells, and that MUC13 can be used to identify compounds that inhibit parasite replication in hepatocytes. This data provides insights into host-parasite interactions in Plasmodium infection, and demonstrates that a component of host mucosal immunity is reprogrammed during the progression of infection.

Ribosome profiling-guided depletion of an mRNA increases cell growth rate and protein secretion.

Recombinant protein production coopts the host cell machinery to provide high protein yields of industrial enzymes or biotherapeutics. However, since protein translation is energetically expensive and tightly controlled, it is unclear if highly expressed recombinant genes are translated as efficiently as host genes. Furthermore, it is unclear how the high expression impacts global translation. Here, we present the first genome-wide view of protein translation in an IgG-producing CHO cell line, measured with ribosome profiling. Through this we found that our recombinant mRNAs were translated as efficiently as the host cell transcriptome, and sequestered up to 15% of the total ribosome occupancy. During cell culture, changes in recombinant mRNA translation were consistent with changes in transcription, demonstrating that transcript levels influence specific productivity. Using this information, we identified the unnecessary resistance marker NeoR to be a highly transcribed and translated gene. Through siRNA knock-down of NeoR, we improved the production- and growth capacity of the host cell. Thus, ribosomal profiling provides valuable insights into translation in CHO cells and can guide efforts to enhance protein production.

Whole-Genome Sequencing of Invasion-Resistant Cells Identifies Laminin α2 as a Host Factor for Bacterial Invasion.

To understand the role of glycosaminoglycans in bacterial cellular invasion, xylosyltransferase-deficient mutants of Chinese hamster ovary (CHO) cells were created using clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated gene 9 (CRISPR-cas9) gene targeting. When these mutants were compared to the pgsA745 cell line, a CHO xylosyltransferase mutant generated previously using chemical mutagenesis, an unexpected result was obtained. Bacterial invasion of pgsA745 cells by group B Streptococcus (GBS), group A Streptococcus, and Staphylococcus aureus was markedly reduced compared to the invasion of wild-type cells, but newly generated CRISPR-cas9 mutants were only resistant to GBS. Invasion of pgsA745 cells was not restored by transfection with xylosyltransferase, suggesting that an additional mutation conferring panresistance to multiple bacteria was present in pgsA745 cells. Whole-genome sequencing and transcriptome sequencing (RNA-Seq) uncovered a deletion in the gene encoding the laminin subunit α2 (Lama2) that eliminated much of domain L4a. Silencing of the long Lama2 isoform in wild-type cells strongly reduced bacterial invasion, whereas transfection with human LAMA2 cDNA significantly enhanced invasion in pgsA745 cells. The addition of exogenous laminin-α2ß1γ1/laminin-α2ß2γ1 strongly increased bacterial invasion in CHO cells, as well as in human alveolar basal epithelial and human brain microvascular endothelial cells. Thus, the L4a domain in laminin α2 is important for cellular invasion by a number of bacterial pathogens. IMPORTANCE: Pathogenic bacteria penetrate host cellular barriers by attachment to extracellular matrix molecules, such as proteoglycans, laminins, and collagens, leading to invasion of epithelial and endothelial cells. Here, we show that cellular invasion by the human pathogens group B Streptococcus, group A Streptococcus, and Staphylococcus aureus depends on a specific domain of the laminin α2 subunit. This finding may provide new leads for the molecular pathogenesis of these bacteria and the development of novel antimicrobial drugs.