A complete workflow for single cell mtDNAseq in CHO cells, from cell culture to bioinformatic analysis.

Chinese hamster ovary (CHO) cells have a long history in the biopharmaceutical industry and currently produce the vast majority of recombinant therapeutic proteins. A key step in controlling the process and product consistency is the development of a producer cell line derived from a single cell clone. However, it is recognized that genetic and phenotypic heterogeneity between individual cells in a clonal CHO population tends to arise over time. Previous bulk analysis of CHO cell populations revealed considerable variation within the mtDNA sequence (heteroplasmy), which could have implications for the performance of the cell line. By analyzing the heteroplasmy of single cells within the same population, this heterogeneity can be characterized with greater resolution. Such analysis may identify heterogeneity in the mitochondrial genome, which impacts the overall phenotypic performance of a producer cell population, and potentially reveal routes for genetic engineering. A critical first step is the development of robust experimental and computational methods to enable single cell mtDNA sequencing (termed scmtDNAseq). Here, we present a protocol from cell culture to bioinformatic analysis and provide preliminary evidence of significant mtDNA heteroplasmy across a small panel of single CHO cells.

CAR Gene Delivery by T-cell Targeted Lentiviral Vectors is Enhanced by Rapamycin Induced Reduction of Antiviral Mechanisms.

Lentiviral vectors (LV) have become the dominant tool for stable gene transfer into lymphocytes including chimeric antigen receptor (CAR) gene delivery to T cells, a major breakthrough in cancer therapy. Yet, room for improvement remains, especially for the latest LV generations delivering genes selectively into T cell subtypes, a key requirement for in vivo CAR T cell generation. Toward improving gene transfer rates with these vectors, whole transcriptome analyses on human T lymphocytes are conducted after exposure to CAR-encoding conventional vectors (VSV-LV) and vectors targeted to CD8+ (CD8-LV) or CD4+ T cells (CD4-LV). Genes related to quiescence and antiviral restriction are found to be upregulated in CAR-negative cells exposed to all types of LVs. Down-modulation of various antiviral restriction factors, including the interferon-induced transmembrane proteins (IFITMs) is achieved with rapamycin as verified by mass spectrometry (LC-MS). Strikingly, rapamycin enhances transduction by up to 7-fold for CD8-LV and CD4-LV without compromising CAR T cell activities but does not improve VSV-LV. When administered to humanized mice, CD8-LV results in higher rates of green fluorescent protein (GFP) gene delivery. Also in vivo CAR T cell generation is improved in kinetics and tumor control, however to a moderate extent, leaving room for improvement by optimizing the rapamycin administration schedule. The data favor multi-omics approaches for improvements in gene delivery.

Light curves and colours of the ejecta from Dimorphos after the DART impact.

On 26 September 2022, the Double Asteroid Redirection Test (DART) spacecraft struck Dimorphos, a satellite of the asteroid 65803 Didymos1. Because it is a binary system, it is possible to determine how much the orbit of the satellite changed, as part of a test of what is necessary to deflect an asteroid that might threaten Earth with an impact. In nominal cases, pre-impact predictions of the orbital period reduction ranged from roughly 8.8 to 17 min (refs. 2,3). Here we report optical observations of Dimorphos before, during and after the impact, from a network of citizen scientists' telescopes across the world. We find a maximum brightening of 2.29 ± 0.14 mag on impact. Didymos fades back to its pre-impact brightness over the course of 23.7 ± 0.7 days. We estimate lower limits on the mass contained in the ejecta, which was 0.3-0.5% Dimorphos's mass depending on the dust size. We also observe a reddening of the ejecta on impact.

Dissecting insect cell heterogeneity during influenza VLP production using single-cell transcriptomics.

The insect cell-baculovirus expression vector system (IC-BEVS) has been widely used to produce recombinant protein at high titers, including complex virus-like particles (VPLs). However, cell-to-cell variability upon infection is yet one of the least understood phenomena in virology, and little is known about its impact on production of therapeutic proteins. This study aimed at dissecting insect cell population heterogeneity during production of influenza VLPs in IC-BEVS using single-cell RNA-seq (scRNA-seq). High Five cell population was shown to be heterogeneous even before infection, with cell cycle being one of the factors contributing for this variation. In addition, infected insect cells were clustered according to the timing and level of baculovirus genes expression, with each cluster reporting similar influenza VLPs transgenes (i.e., hemagglutinin and M1) transcript counts. Trajectory analysis enabled to track infection progression throughout pseudotime. Specific pathways such as translation machinery, protein folding, sorting and degradation, endocytosis and energy metabolism were identified as being those which vary the most during insect cell infection and production of Influenza VLPs. Overall, this study lays the ground for the application of scRNA-seq in IC-BEVS processes to isolate relevant biological mechanisms during recombinant protein expression towards its further optimization.

Impact of dual-baculovirus infection on the Sf9 insect cell transcriptome during rAAV production using single-cell RNA-seq.

The insect cell-baculovirus expression vector system (IC-BEVS) has shown to be a powerful platform to produce complex biopharmaceutical products, such as recombinant proteins and virus-like particles. More recently, IC-BEVS has also been used as an alternative to produce recombinant adeno-associated virus (rAAV). However, little is known about the variability of insect cell populations and the potential effect of heterogeneity (e.g., stochastic infection process and differences in infection kinetics) on product titer and/or quality. In this study, transcriptomics analysis of Sf9 insect cells during the production of rAAV of serotype 2 (rAAV2) using a low multiplicity of infection, dual-baculovirus system was performed via single-cell RNA-sequencing (scRNA-seq). Before infection, the principal source of variability in Sf9 insect cells was associated with the cell cycle. Over the course of infection, an increase in transcriptional heterogeneity was detected, which was linked to the expression of baculovirus genes as well as to differences in rAAV transgenes (rep, cap and gfp) expression. Noteworthy, at 24 h post-infection, only 29.4% of cells enclosed all three necessary rAAV transgenes to produce packed rAAV2 particles, indicating limitations of the dual-baculovirus system. In addition, the trajectory analysis herein performed highlighted that biological processes such as protein folding, metabolic processes, translation, and stress response have been significantly altered upon infection. Overall, this work reports the first application of scRNA-seq to the IC-BEVS and highlights significant variations in individual cells within the population, providing insight into the rational cell and process engineering toward improved rAAV2 production in IC-BEVS.

Transcriptome analysis of Sf9 insect cells during production of recombinant Adeno-associated virus.

The insect cell-baculovirus expression vector system (IC-BEVS) has emerged as an alternative time- and cost-efficient production platform for recombinant Adeno-associated virus (AAV) for gene therapy. However, a better understanding of the underlying biological mechanisms of IC-BEVS is fundamental to further optimize this expression system toward increased product titer and quality. Here, gene expression of Sf9 insect cells producing recombinant AAV through a dual baculovirus expression system, with low multiplicity of infection (MOI), was profiled by RNA-seq. An 8-fold increase in reads mapping to either baculovirus or AAV transgene sequences was observed between 24 and 48 h post-infection (hpi), confirming a take-over of the host cell transcriptome by the baculovirus. A total of 336 and 4784 genes were identified as differentially expressed at 24 hpi (vs non-infected cells) and at 48 hpi (vs. infected cells at 24 hpi), respectively, including dronc, birc5/iap5, and prp1. Functional annotation found biological processes such as cell cycle, cell growth, protein folding, and cellular amino acid metabolic processes enriched along infection. This work uncovers transcriptional changes in Sf9 in response to baculovirus infection, which provide new insights into cell and/or metabolic engineering targets that can be leveraged for rational bioprocess engineering of IC-BEVS for AAV production.

Gene Expression Analysis of Adapted Insect Cells during Influenza VLP Production Using RNA-Sequencing.

Adaptive laboratory evolution has been used to improve production of influenza hemagglutinin (HA)-displaying virus-like particles (VLPs) in insect cells. However, little is known about the underlying biological mechanisms promoting higher HA-VLP expression in such adapted cell lines. In this article, we present a study of gene expression patterns associated with high-producer insect High Five cells adapted to neutral pH, in comparison to non-adapted cells, during expression of influenza HA-VLPs. RNA-seq shows a decrease in the amount of reads mapping to host cell genomes along infection, and an increase in those mapping to baculovirus and transgenes. A total of 1742 host cell genes were found differentially expressed between adapted and non-adapted cells throughout infection, 474 of those being either up- or down-regulated at both time points evaluated (12 and 24 h post-infection). Interestingly, while host cell genes were found up- and down-regulated in an approximately 1:1 ratio, all differentially expressed baculovirus genes were found to be down-regulated in infected adapted cells. Pathway analysis of differentially expressed genes revealed enrichment of ribosome biosynthesis and carbohydrate, amino acid, and lipid metabolism. In addition, oxidative phosphorylation and protein folding, sorting and degradation pathways were also found to be overrepresented. These findings contribute to our knowledge of biological mechanisms of insect cells during baculovirus-mediated transient expression and will assist the identification of potential engineering targets to increase recombinant protein production in the future.

Erratum: Monitoring CAR T cell generation with a CD8-targeted lentiviral vector by single-cell transcriptomics.

[This corrects the article DOI: 10.1016/j.omtm.2021.09.019.].

Monitoring CAR T cell generation with a CD8-targeted lentiviral vector by single-cell transcriptomics.

Quantifying gene expression in individual cells can substantially improve our understanding about complex genetically engineered cell products such as chimeric antigen receptor (CAR) T cells. Here we designed a single-cell RNA sequencing (scRNA-seq) approach to monitor the delivery of a CD19-CAR gene via lentiviral vectors (LVs), i.e., the conventional vesicular stomatitis virus (VSV)-LV and the CD8-targeted CD8-LV. LV-exposed human donor peripheral blood mononuclear cells (PBMCs) were evaluated for a panel of 400 immune response-related genes including LV-specific probes. The resulting data revealed a trimodal expression for the CAR and CD8A, demanding a careful distribution-based identification of CAR T cells and CD8+ lymphocytes in scRNA-seq analysis. The fraction of T cells expressing high CAR levels was in concordance with flow cytometry results. More than 97% of the cells hit by CD8-LV expressed the CD8A gene. Remarkably, the majority of the potential off-target cells were in fact on-target cells, resulting in a target cell selectivity of more than 99%. Beyond that, differential gene expression analysis revealed the upregulation of restriction factors in CAR-negative cells, thus explaining their protection from CAR gene transfer. In summary, we provide a workflow and subsetting approach for scRNA-seq enabling reliable distinction between transduced and untransduced cells during CAR T cell generation.

Tracing production instability in a clonally derived CHO cell line using single-cell transcriptomics.

A variety of mechanisms including transcriptional silencing, gene copy loss, and increased susceptibility to cellular stress have been associated with a sudden or gradual loss of monoclonal antibody (mAb) production in Chinese hamster ovary (CHO) cell lines. In this study, we utilized single-cell RNA-seq (scRNA-seq) to study a clonally derived CHO cell line that underwent production instability leading to a dramatic reduction of the levels of mAb produced. From the scRNA-seq data, we identified subclusters associated with variations in the mAb transgenes and observed that heavy chain gene expression was significantly lower than that of the light chain across the population. Using trajectory inference, the evolution of the cell line was reconstructed and was found to correlate with a reduction in heavy and light chain gene expression. Genes encoding for proteins involved in the response to oxidative stress and apoptosis were found to increase in expression as cells progressed along the trajectory. Future studies of CHO cell lines using this technology have the potential to dramatically enhance our understanding of the characteristics underpinning efficient manufacturing performance as well as product quality.

From Verbal Account to Written Evidence: Do Written Statements Generated by Officers Accurately Represent What Witnesses Say?

Most countries compile evidence from witnesses and victims manually, whereby the interviewer assimilates what the interviewee says during the course of an interview to produce an evidential statement. This exploratory research examined the quality of evidential statements generated in real world investigations. Transcribed witness/victim interviews (N = 15) were compared to the resultant written statements produced by the interviewing officer and signed as an accurate record by the interviewee. A coding protocol was devised to assess the consistency of information between what was said by the interviewee in the verbal interview and what was reported in the written statement. Statements contained numerous errors including omissions, distortions, and the inclusion of information not mentioned in the verbal interview. This exploratory work highlights an important area for future research focus.

Expanding the Chinese hamster ovary cell long noncoding RNA transcriptome using RNASeq.

Our ability to study Chinese hamster ovary (CHO) cell biology has been revolutionised over the last decade following the development of next generation sequencing technology and publication of reference DNA sequences for CHO cells and the Chinese hamster. RNA sequencing has not only enabled the association of transcript expression with bioreactor conditions and desirable bioprocess phenotypes but played a key role in the characterisation of protein coding and small noncoding RNAs. The annotation of long noncoding RNAs, and therefore our understanding of their role in CHO cell biology, has been limited to date. In this manuscript, we use high-resolution RNASeq data to more than double the number of annotated lncRNA transcripts for the CHO K1 genome. In addition, the utilisation of strand-specific sequencing enabled the identification of more than 1,000 new antisense and divergent lncRNAs. The utility of monitoring lncRNA expression is demonstrated through an analysis of the transcriptomic response to a reduction of cell culture temperature and identification of simultaneous sense/antisense differential expression for the first time in CHO cells. To enable further studies of lncRNAs, the transcripts annotated in this study have been made available for the CHO cell biology community.

Subphysiological temperature induces pervasive alternative splicing in Chinese hamster ovary cells.

RNA sequencing (RNASeq) has been widely used to associate alterations in Chinese hamster ovary (CHO) cell gene expression with bioprocess phenotypes; however, alternative messenger RNA (mRNA) splicing, has thus far, received little attention. In this study, we utilized RNASeq for transcriptomic analysis of a monoclonal antibody (mAb) producing CHO K1 cell line subjected to a temperature shift. More than 2,465 instances of differential splicing were observed 24 hr after the reduction of cell culture temperature. A total of 1,197 of these alternative splicing events were identified in genes where no changes in abundance were detected by standard differential expression analysis. Ten examples of alternative splicing were selected for independent validation using quantitative polymerase chain reaction in the mAb-producing CHO K1 cell line used for RNASeq and a further two CHO K1 cell lines. This analysis provided evidence that exon skipping and mutually exclusive splicing events occur in genes linked to the cellular response to changes in temperature and mitochondrial function. While further work is required to determine the impact of these changes in mRNA sequence on cellular phenotype, this study demonstrates that alternative splicing analysis can be utilized to gain a deeper understanding of post-transcriptional regulation in CHO cells during biopharmaceutical production.

Genetic and Epigenetic Variation across Genes Involved in Energy Metabolism and Mitochondria of Chinese Hamster Ovary Cell Lines.

The increasingdemandfor biopharmaceutical products drives the search for efficient cell factories that are able to sustainably support rapid growth, high productivity, and product quality. As these depend on energy generation, here the genomic variation in nuclear genes associated with mitochondria and energy metabolism and the mitochondrial genome of 14 cell lines is investigated. The variants called enable reliable tracing of lineages. Unique sequence variations are observed in cell lines adapted to grow in protein-free media, enriched in signaling pathways or mitogen-activated protein kinase 3. High-producing cell lines bear unique mutations in nicotinamide adenine dinucleotide (NADH) dehydrogenase (ND2 and ND4) and in peroxisomal acyl-CoA synthetase (ACSL4), involved in lipid metabolism. As phenotypes are determined not only by functional mutations, but also by the exquisite regulation of expression patterns, it is not surprising that ≈50% of the genes investigated here are found to be differentially methylated and thus epigenetically controlled, enabling a clear distinction of high producers, and cells adapted to a minimal, glutamine (Gln)-free medium. Similar pathways are enriched as those identified by genome variation. This strengthens the hypothesis that these phenomena act together to define cell behavior.

Transcriptomic analysis of IgG4 Fc-fusion protein degradation in a panel of clonally-derived CHO cell lines using RNASeq.

In this study, we report an investigation of a panel of clonally-derived Chinese hamster ovary (CHO) cell lines exhibiting variability in the proportion of full-length IgG4 Fc-fusion protein produced. The recombinant protein was found to be degraded during cell culture into four shorter "clipped" species (three of the four cleavage sites occurred at arginine residues) and preliminary analyses suggested that a host cell enzyme was responsible for proteolysis. To identify the specific enzyme responsible, RNA sequencing was used to identify gene expression differences between the cell lines with a "high" and "low" clipping phenotype. From this analysis, six protease-encoding genes were found to be significantly upregulated in those cell lines yielding the lowest proportion of full-length IgG4 Fc-fusion protein. Four of these protease candidates were deprioritized after examination of their cleavage site specificity. The remaining enzymes, Adam19 and Furin, were found to be capable of cleavage at arginine residues, and inhibitors for both proteases were added to cell-free media to determine if the product degradation could be reduced. While the Adam19 inhibitor had no impact, Furin inhibitor I (specific for the proprotein convertase family of enzymes) was found to result in a 33-39% increase in complete IgG4 Fc-fusion protein when compared with untreated samples.

Bicaudal D1 impairs autophagosome maturation in chronic obstructive pulmonary disease.

Bicaudal D1 (BICD1), an adaptor for the dynein-dynactin motor complex, has been identified as a susceptibility gene in chronic obstructive pulmonary disease (COPD). Autophagy, an essential cellular homeostasis process, is defective in COPD, in which oxidative stress-induced misfolded proteins accumulate into toxic aggregates dependent on the accumulation of the autophagic cargo receptor p62. Defective autophagy can be caused by mutations in the dynein and dynactin motor complex suggesting a possible link between BICD1 and defective autophagy in COPD. BICD1 levels were measured in peripheral lung tissue from COPD patients together with markers of autophagy and found to be increased in COPD together with autophagosomes, p62 and p62 oligomers. In vitro exposure of bronchial epithelial cells to cigarette smoke extracts (CSEs) revealed that high concentrations of CSE induced defective autophagosome maturation with accumulation of BICD1, p62 and ubiquitin-associated p62 oligomers. This was confirmed in vivo using CS-exposed mice. Furthermore, we identified that formation of CS-induced p62 oligomers required an interaction with Keap1. Overexpression and ablation of BICD1 confirmed that increased BICD1 negatively regulates autophagosome maturation inducing accumulation of p62 and p62 oligomers and that it can be reversed by cardiac glycosides. We conclude that defective autophagosome maturation in COPD is caused by oxidative stress-mediated BICD1 accumulation.

Understanding biopharmaceutical production at single nucleotide resolution using ribosome footprint profiling.

Biopharmaceuticals such as monoclonal antibodies have revolutionised the treatment of a variety of diseases. The production of recombinant therapeutic proteins, however, remains expensive due to the manufacturing complexity of mammalian expression systems and the regulatory burden associated with administrating these medicines to patients in a safe and efficacious manner. In recent years, academic and industrial groups have begun to develop a greater understanding of the biology of host cell lines, such as Chinese hamster ovary (CHO) cells and utilise that information for process development and cell line engineering. In this review, we focus on ribosome footprint profiling (RiboSeq), an exciting next generation sequencing (NGS) method that provides genome-wide information on translation, and discuss how its application can transform our understanding of therapeutic protein production.

Parallel mRNA, proteomics and miRNA expression analysis in cell line models of the intestine.

AIM: To identify miRNA-regulated proteins differentially expressed between Caco2 and HT-29: two principal cell line models of the intestine. METHODS: Exponentially growing Caco-2 and HT-29 cells were harvested and prepared for mRNA, miRNA and proteomic profiling. mRNA microarray profiling analysis was carried out using the Affymetrix GeneChip Human Gene 1.0 ST array. miRNA microarray profiling analysis was carried out using the Affymetrix Genechip miRNA 3.0 array. Quantitative Label-free LC-MS/MS proteomic analysis was performed using a Dionex Ultimate 3000 RSLCnano system coupled to a hybrid linear ion trap/Orbitrap mass spectrometer. Peptide identities were validated in Proteome Discoverer 2.1 and were subsequently imported into Progenesis QI software for further analysis. Hierarchical cluster analysis for all three parallel datasets (miRNA, proteomics, mRNA) was conducted in the R software environment using the Euclidean distance measure and Ward's clustering algorithm. The prediction of miRNA and oppositely correlated protein/mRNA interactions was performed using TargetScan 6.1. GO biological process, molecular function and cellular component enrichment analysis was carried out for the DE miRNA, protein and mRNA lists via the Pathway Studio 11.3 Web interface using their Mammalian database. RESULTS: Differential expression (DE) profiling comparing the intestinal cell lines HT-29 and Caco-2 identified 1795 Genes, 168 Proteins and 160 miRNAs as DE between the two cell lines. At the gene level, 1084 genes were upregulated and 711 were downregulated in the Caco-2 cell line relative to the HT-29 cell line. At the protein level, 57 proteins were found to be upregulated and 111 downregulated in the Caco-2 cell line relative to the HT-29 cell line. Finally, at the miRNAs level, 104 were upregulated and 56 downregulated in the Caco-2 cell line relative to the HT-29 cell line. Gene ontology (GO) analysis of the DE mRNA identified cell adhesion, migration and ECM organization, cellular lipid and cholesterol metabolic processes, small molecule transport and a range of responses to external stimuli, while similar analysis of the DE protein list identified gene expression/transcription, epigenetic mechanisms, DNA replication, differentiation and translation ontology categories. The DE protein and gene lists were found to share 15 biological processes including for example epithelial cell differentiation [P value ≤ 1.81613E-08 (protein list); P ≤ 0.000434311 (gene list)] and actin filament bundle assembly [P value ≤ 0.001582797 (protein list); P ≤ 0.002733714 (gene list)]. Analysis was conducted on the three data streams acquired in parallel to identify targets undergoing potential miRNA translational repression identified 34 proteins, whose respective mRNAs were detected but no change in expression was observed. Of these 34 proteins, 27 proteins downregulated in the Caco-2 cell line relative to the HT-29 cell line and predicted to be targeted by 19 unique anti-correlated/upregulated microRNAs and 7 proteins upregulated in the Caco-2 cell line relative to the HT-29 cell line and predicted to be targeted by 15 unique anti-correlated/downregulated microRNAs. CONCLUSION: This first study providing "tri-omics" analysis of the principal intestinal cell line models Caco-2 and HT-29 has identified 34 proteins potentially undergoing miRNA translational repression.