Skip to main content

Characterization of metabolic responses, genetic variations, and microsatellite instability in ammonia-stressed CHO cells grown in fed-batch cultures



As bioprocess intensification has increased over the last 30 years, yields from mammalian cell processes have increased from 10’s of milligrams to over 10’s of grams per liter. Most of these gains in productivity can be attributed to increasing cell densities within bioreactors. As such, strategies have been developed to minimize accumulation of metabolic wastes, such as lactate and ammonia. Unfortunately, neither cell growth nor biopharmaceutical production can occur without some waste metabolite accumulation. Inevitably, metabolic waste accumulation leads to decline and termination of the culture. While it is understood that the accumulation of these unwanted compounds imparts a suboptimal culture environment, little is known about the genotoxic properties of these compounds that may lead to global genome instability. In this study, we examined the effects of high and moderate extracellular ammonia on the physiology and genomic integrity of Chinese hamster ovary (CHO) cells.


Through whole genome sequencing, we discovered 2394 variant sites within functional genes comprised of both single nucleotide polymorphisms and insertion/deletion mutations as a result of ammonia stress with high or moderate impact on functional genes. Furthermore, several of these de novo mutations were found in genes whose functions are to maintain genome stability, such as Tp53, Tnfsf11, Brca1, as well as Nfkb1. Furthermore, we characterized microsatellite content of the cultures using the CriGri-PICR Chinese hamster genome assembly and discovered an abundance of microsatellite loci that are not replicated faithfully in the ammonia-stressed cultures. Unfaithful replication of these loci is a signature of microsatellite instability. With rigorous filtering, we found 124 candidate microsatellite loci that may be suitable for further investigation to determine whether these loci may be reliable biomarkers to predict genome instability in CHO cultures.


This study advances our knowledge with regards to the effects of ammonia accumulation on CHO cell culture performance by identifying ammonia-sensitive genes linked to genome stability and lays the foundation for the development of a new diagnostic tool for assessing genome stability.


Biopharmaceutical manufacturing represents nearly 2% of the total US GDP [1] which makes it an important driver of the US economy. Biopharmaceuticals include monoclonal antibodies, recombinant proteins, and assemblies of proteins produced by biological means. Commercial products are used as blood factors, thrombolytic agents, therapeutics, growth factors, interferons, and vaccines [2, 3]. The most common mammalian cell line used is the Chinese hamster ovary (CHO) cell line, due to its ability to produce biopharmaceutical molecules with post-translational modifications required in humans [4]. However, it is well understood that recombinant CHO cell lines are susceptible to genome instability that is often observed after approximately 70 generations [5,6,7,8]. Previous studies have characterized genomic variants across various CHO cell lines that may be a contributing factor to genome instability [9,10,11]. An unstable genome can result in reduced productivity in continuous cultures and fed-batch systems [12, 13]. A common occurrence in both continuous cultures and fed-batch systems is the accumulation of metabolic waste products, such as ammonia and lactate. The role these waste products play in cellular processes, such as glycosylation, metabolism, and productivity have been characterized [14,15,16,17]; however, the effects of these waste products on genome stability have not been directly assessed.

Microsatellite instability (MSI) is described as genetic hypermutability at microsatellite loci where a high frequency of insertion or deletion (indel) mutations accumulate in daughter cells during cell division [18, 19]. MSI results from improperly functioning mismatch repair (MMR) pathways which are key to maintaining genome stability [20]. Rather than correcting DNA mismatch errors that occur spontaneously during DNA replication, cells with impaired MMR systems accumulate these errors over the course of subsequent propagation. The prevalence of these errors allow for MSI loci to be utilized as stable genetic biomarkers that are capable of diagnosing many human cancers [21, 22]. Studies have shown that approximately 15% of human patients with colorectal cancer [20, 23], 20% of patients with stomach cancer [24], and 30% of patients with endometrial cancer [25] could attribute their disease to genome instability that can be diagnosed with MSI biomarkers. The clinical uptake of MSI-based diagnostics, such as the Bethesda Panel, demonstrates the reliability and clinical utility of MSI loci as biomarkers [26].

In this study, we investigated the effects of exogenous ammonia exposure on genome stability during fed-batch cultures of CHO cells. Specifically, the accumulation of DNA mutations in cells exposed to elevated ammonia were compared to cultures grown under standard fed-batch conditions. Ammonia was added to duplicate parallel cultures at 10 mM and 30 mM final concentrations at 12 h of culture time to establish mild and high ammonia stresses respectively. After 72 h of elevated ammonia exposure, samples were taken for whole genome sequencing (WGS). These sequences were then analyzed for MSI, single nucleotide polymorphisms (SNPs), and insertion/deletion (indel) variations. The SNPs and indels were mapped to the Chinese hamster genome and assessed for functional impact in both coding and regulatory genetic regions. Microsatellite regions were analyzed to identify loci with dose-dependent indel mutations that could be used as potential biomarkers.

Materials and methods

Culture conditions

A recombinant CHO-K1 Clone A11 from the Vaccine Research Center at the National Institutes of Health (NIH), which expresses the anti-HIV antibody VRC01 (IgG1) was used. The inoculum train was expanded in 250 mL shake flasks with 70 mL ActiPro media (GE Healthcare) that were maintained at 5% CO2 and 37 °C. The bioreactors were ambr® 250 bioreactors (Sartorius Stedim, Göttingen, Germany) with two pitched blade impellers and an open pipe sparger (vessel part number: 001-5G25). The bioreactors were inoculated at a target cell density of 0.4 × 106 cells/mL in ActiPro batch media and fed daily beginning on Day 3 (3% (v/v) Boost 7a and 0.3% (v/v) Boost 7b (GE Healthcare)). Duplicate cultures were ammonia-stressed at 12-h post inoculation with 0 mM, 10 mM, or 30 mM NH4Cl. The 0 mM and 10 mM cultures used saline to normalize the volume of the 0 mM and 10 mM cultures to the 30 mM cultures. Dissolved oxygen was controlled at 50% of air saturation using PID control that increased the O2 mixture in the gas sparge to 100%, then the stir speed from 300 to 600 rpm. Antifoam (10% solution in media; SH30897.41 – GE Healthcare) was added as needed to control foaming. All gases were supplied through the open pipe sparger; an overlay was not used. The pH was controlled via sparging CO2 and air, and base pump (1 M NaOH). The pH setpoint was 7.0 with a 0.2 deadband. Temperature was controlled at 37 °C. Samples for WGS and MSI analysis were harvested at 84 h culture time (72 h post-stress) and centrifuged at approximately 2000 x g for 15 min at 4 °C. The supernatant was removed, and the pellet was stored at − 80 °C.

DNA extraction, whole genome sequencing, and microsatellite variant discovery

Pellets of approximately 0.5 × 106 cells were pre-washed with 1X phosphate buffered saline (PBS) prior to extraction. Total genomic DNA (gDNA) was purified from 2 replicate samples per condition with the DNAeasy Blood and Tissue Kit (Qiagen), following the manufacturer’s recommended procedures and combined prior to sequencing. Whole genome shotgun sequencing was performed on an Illumina NovaSeq (2 × 150 paired end) through a third- party vendor to approximately 30x genome coverage. Raw sequence data was assessed for quality with the FASTQC software ( Raw sequence data was preprocessed to remove low quality bases and adapter sequences with the Trimmomatic software v.0.38 [27]. Preprocessed reads were aligned to the CriGri-PICR version assembly of the Chinese hamster genome (Cricetulus griseus) (RefSeq assembly accession: GCF_003668045.1) with the Bowtie2 v. short read aligner [28]. Alignments were coordinate sorted and indexed with SamTools v1.3.1 [29]. SNPs and indels were determined with the HaploTypeCaller Walker from the Genome Analysis Toolkit (GATK v.4.0) [30]. Genetic variants in the resulting VCF file (SNP and INDEL) were hard filtered according to the following criteria: read depth (DP 6) and mapping quality (MQ 30). Variant sites were further filtered to remove variant loci that were common between the sample groups but differed from the CHO PICR reference assembly. Finally, variant sites were kept in the final VCF file only if one or both of the treatment samples differed from the control. Functional SNPs were characterized with the SnpEff software, v4.3 [31]. Genome-wide microsatellite loci were determined against the PICR CH assembly with MISA, a microsatellite finder software [32]. Microsatellite loci were intersected with indel coordinates using BedTools Intersect command 2.27.1 [33] to identify indel variants associated with microsatellites.

Text/data mining and functional enrichment analysis

The query “genomic instability [MeSH Terms] “was used to search PubMed to retrieve the abstracts with PMIDs (14,968 PMIDs). The PubTator [34] tool was used to collect genes annotated in these abstracts with Entrez Gene IDs ( released 2/14/2020). Among the 5073 retrieved genes, 3131 genes were human (Homo sapiens), 882 mouse (Mus musculus), 435 yeast (Saccharomyces cerevisiae) representing the three top species with the largest number of genes mapped in those PubMed abstracts. The other top mapped species include seven vertebrates: rat, chicken, zebrafish, frog, Chinese hamster, dog and pig as well as three non-vertebrates: fly, Arabidopsis, and worm. The ortholog pairs between human and the eight other above vertebrates were mapped with NCBI ortholog assignment ( released 07/20/2020), whereas those between human and the four non-vertebrate species were mapped with OMAbrowser ( Altogether, 2897 human genes linked to “genomic instability” were matched with corresponding Chinese hamster orthologs. For the SNP list with high and moderate impact mutations by SnpEff, 273 Chinese hamster genes were mapped with their corresponding human orthologs that had been linked to “genomic instability”. ClusterProfiler [35, 36] was used to obtain the enriched KEGG pathways and GO annotations for the given gene lists.

Identification of candidate MSI loci

Candidate MSI loci were determined with a filtering strategy that leverages several criteria as follows: First, each novel indel-variable genomic locus was assigned a mutation score which is a proportion of the number of variant reads (allelic depth) by the total depth of reads for each site extracted from the .vcf file. Second, the mutation scores of the control cultures were subtracted from the mutation scores of the treated loci in order to generate a mutation score relative to the control. This allowed for the removal of loci that did not exhibit dose-dependent responses to the exogenous ammonia. Concurrently, loci with nonpositive relative mutation scores were also removed. The remaining loci were then intersected with the genome-wide microsatellite coordinates determined by MISA with the Intersect command of BedTools v2.27.1 [33] to identify loci within known microsatellites. The final ranked set of candidate MSI loci contain sites where control samples have fewer to no variant reads in comparison to the treated samples.


Growth and metabolite profiles

Recombinant CHO cells expressing the monoclonal antibody VRC01 were cultured in tightly controlled ambr® 250 bioreactors for 12 h prior to the addition of ammonia to stress the cultures. Up to 1.5 days post-inoculation, there were no observable differences in the viable cell densities (VCD); however, at 2.5 days, the 30 mM ammonia-stressed cultures had substantially lower VCDs compared to the control and 10 mM stressed cultures (Fig. 1a). The 10 mM ammonia-stressed cultures had similar VCDs to the control cultures until Day 7; yet cell viabilities were similar to the control cultures for the entire culture durations. In contrast, the 30 mM ammonia-stressed cultures reached peak VCDs on Day 4 and gradually declined until the cultures were harvested on Day 8.5 due to low viability (< 70%); a cell viability below 70% is a standard harvesting threshold. Samples for genome sequencing were taken 84 h post inoculation (Day 3.5), i.e., 72 h post-stress. At the time of harvest of genome sequencing samples, the viability for all samples was greater than 90% (Fig. 1a). The mildly stressed (10 mM) cultures had no significant change in the ammonia levels between 12 and 84-h, while the high stress (30 mM) cultures had a gradual decline in the ammonia concentration until harvested (Fig. 1b). The glucose and lactate profiles (Fig. 1d, e) confirm that the control and 10 mM ammonia-stressed cultures were paired closely throughout the entire cultures, although the 10 mM stressed culture had slightly lower VCD beginning on Day 6. In contrast, the 30 mM ammonia-stressed cultures began to accumulate glucose and lactate after Day 5, most likely due to the set feeding protocol based on culture volume, and the significantly lower cell growth (Fig. 1a). It is well-known that excessive glucose inevitably leads to lactate accumulation [15], which was observed for the 30 mM stressed cultures. Amino acid profiles were also obtained for these cultures [37]. The amino acid profile that showed the greatest differences between the control and 10 mM cultures was alanine; both the 10 mM and 30 mM ammonia-stressed culture alanine profiles were very similar through Day 6, while the control cultures had profiles that represented a higher consumption rate, as alanine was fed starting on Day 3 (Fig. 1c). The glutamine profiles for the control and 10 mM cultures were similar up to Day 7, until the 10 mM cultures began to accumulate glutamine (Fig. 1g). The glutamine accumulation can be attributed to the feeding of glutamate (Fig. 1h), which when in excess can be aminated to form glutamine [38]. The 30 mM cultures were terminated prior to any significant differences in the glutamine accumulation being observed between the control and 30 mM cultures. Therefore, the glutamate feeding, based on volume, caused glutamate to accumulate due to lower VCD relative to the control cultures for the 10 mM and 30 mM cultures. This in turn impacted the glutamine profile. A global measure of cell health is the overall protein production and cell-specific productivity. The monoclonal antibody titer at the end of the cultures was about 50% lower for the 10 mM cultures compared to the control cultures, whereas the 30 mM cultures had negligible protein productivity in (Fig. 1f). Furthermore, cell-specific productivity (picograms of IgG per cell per day) was found to be substantially higher in the control cultures when compared to the 10 mM sample. It should also be noted that the control and 10 mM cultures had a relatively stable production rate, whereas the 30 mM cultures declined (Fig. 1i). Overall, the samples for the genome sequencing analysis were taken at culture times when there were no substantial VCD, viability, or metabolic differences between the control and 10 mM ammonia-stressed cultures; however, the VCD was significantly lower for the 30 mM ammonia-stressed cultures (Fig. 1).

Fig. 1

Cell growth, ammonia, titer and metabolic profiles for CHO K-1 VRC01 cells cultured in duplicate in the ambr® 250 bioreactor. The ammonia stresses (10 mM and 30 mM) were added at 12 h. Samples for genomic analysis were harvested at 84 h (3.5 days) as shown by the solid grey line. a Viable cell density (VCD) and viability (filled and hollow symbols, respectively), b ammonia, c alanine, d glucose, e lactate, f titer of recombinant monoclonal antibody, g glutamine, h glutamate, and i Cell specific productivity (qp). Due to low levels of the recombinant protein in culture prior to day 5, the qp value is not shown until a significant titer has been reached. In industry, it is common to only measure titers starting at day 7. Control - 0 mM (blue and purple lines); 10 mM (orange and red lines); 30 mM (green and dark green lines)

Whole genome shotgun sequencing and variant discovery in stressed conditions

Whole genome shotgun sequences were collected for the control and treated samples to an approximate depth of 30X coverage to assess the genomic impact of exogenous ammonia exposure. A total of 389,694 variant sites were identified across both stress levels that were composed of 310,597 SNPs and 79,097 indels, (Supplemental Tables S1 and S2 respectively). Of the 389,694 variant sites, a total of 135,913 variant sites reside in protein coding genes (Supplemental Table S3). The variant sites were seemingly randomly distributed in both intergenic and genic positions across the genome. A distribution and density map of variant positions relative to annotated coding genes is depicted in Fig. 2. These variants were further filtered to remove sites annotated as a modifier or low impact variant (e.g. synonymous mutations) predicted by SnpEff. This led to the discovery of 2394 variants within protein coding genes with significant impact variations (high/moderate impact predicted by SnpEff) due to ammonia stress (Supplemental Table S4).

Fig. 2

Genome coverage map of genetic variants and MSI loci in the 10 longest CHO scaffolds. A circos plot of the 10 longest CHO genome scaffolds (parsed into 100 kb windows that represents approximately ~ 20% of the CHO genome) that depicts distribution of genetic variants and MSI loci. The innermost track (light purple) depicts candidate MSI loci (23 out of the 124); Indel density and distribution is depicted in light blue; SNP density and distribution is in light red; and gene density is plotted outside of the CHO ideograms in light orange. Ideogram ticks are scaled in megabases

Functional impact of ammonia-induced variants in genome stability genes

The above described 2394 variants were assigned to 1843 Chinese hamster protein-coding genes with certain functional impact. Through mapping of human orthologs for those Chinese hamster genes, we found 273 genes that are linked to genome instability terms via text mining (Supplemental Table S5). Figure 3a shows the map for KEGG enrichment result of over-representation test of the 273 genes. The five most significant KEGG pathways include breast cancer, cellular senescence, longevity regulating pathway, MAPK signaling pathway, and cell cycle. It is critical to note that the KEGG enrichment analysis (Fig. 3a) combines all variants found in the 10 mM and 30 mM stress cultures, whereas gene lists for variants exclusively detected in 10 mM or 30 mM stress samples generated no enrichment of KEGG pathways. Fig. 3b shows the KEGG comparison between three gene lists: one for all variants in 10 mM sample, another for 30 mM sample, and one for variants from the combined list (i.e., the above mentioned 273 genes). Breast cancer, cellular senescence, longevity regulating pathway are the three KEGG pathways common to the three gene lists (Fig. 3b). These corresponding genes are listed in Table 1. Figure 4 summarizes the significant GO terms enriched among the genes for variants existing at both stress levels. Notable GO biological process terms in Fig. 4 include DNA recombination, cell cycle checkpoint, regulation of response to DNA damage stimulus, telomere organization, and DNA damage checkpoint. Additionally, notable functions of genes include double-strand break repair (Brca1), mismatch repair (Mlh3), and centromere generation (Cenpc) (Supplemental Table S5). More detailed enrichment analysis on KEGG and GO of variant genes can be found in Supplemental Tables S6, S7, S8 and S9.

Fig. 3

KEGG enrichment results from over representation analysis of ammonia-sensitive genes linked to genomic instability. a Network plot of most significant enriched KEGG pathways. Enrichment scores (i.e., adjusted p values) and gene counts (the number of genes enriched in a KEGG term) are depicted by dot color and size. b Comparison of enrichment results of KEGG pathways for genes identified with significant variants in the 10 mM Ammonia-stressed culture (111 genes), in the 30 mM culture (102 genes), and in the union of 10 mM and 30 mM cultures (148 genes). Enrichment scores (i.e., adjusted p values) and gene ratios (the percentage of total genes in the given KEGG term) are depicted by dot color and size. The plots are made with clusterProfiler

Table 1 A summary of select KEGG enrichment genes discovered in ammonia-stressed cultures that can be linked to genome instability in humans via text mining
Fig. 4

Gene Ontology (GO) enrichment results from over representation analysis of ammonia-sensitive genes linked to genomic instability. Terms in the three GO categories are grouped with color bars: biological process (BP, red), cellular component (CC, orange), and molecular function (MF, blue). Fold enrichment is shown in x-axis. Enrichment scores (i.e., adjusted p values) and gene counts (the number of genes enriched in a GO term) are depicted by dot color and size. The plot is made with clusterProfiler

Microsatellite and candidate MSI loci

A whole-genome scan for microsatellites discovered a total of 409,628 loci, with motifs that included di-, tri-, and tetranucleotide repeats (Supplemental Table S10). As expected, the microsatellites composed of dinucleotide repeats were the most prevalent with a total of 287,124. Trinucleotide and tetranucleotide motifs were less abundant with 46,602 and 75,902 occurrences, respectively. An analysis of genome-wide indels in ammonia treated and control samples revealed 1022 microsatellites that were lengthened or shortened due to the ammonia stress (Supplemental Table S11). An example microsatellite locus with desirable length variation resulting from elevated ammonia is shown in Fig. 5. In this example, there is a higher abundance of mapped reads with deletions for the 30 mM ammonia stressed cultures, suggesting a dose-dependent response. Furthermore, we developed a custom mutation score and stringent filtering criteria (see Methods) to identify a candidate set of 124 MSI loci where stable mutations were present in both ammonia-stressed cultures, but were not present in the control cultures. These 124 MSI loci can be used as a foundation for future research as diagnostics for genome instability (Fig. 2). It is important to note that because the mutation score was calculated using the allelic depth; loci with more reads are statistically more significant than those with fewer reads. With this in mind, the 124 candidate MSI loci may not be all inclusive of the optimal loci due to the variation in mapped read depth across the genome. The remaining loci after each filter step is summarized in Table 2. A full list of loci in each step can be found in Supplemental Tables S12, S13, S14 and S15, while the location of all candidate loci are summarized in Table 3.

Fig. 5

An Integrated Genome Viewer (IGV) image of a microsatellite located on scaffold NW_020822544.1 at position 4,160,116. a control, b 10 mM, and c 30 mM cultures. This microsatellite contains five repeats of an AAC motif. They grey lines above each nucleotide are indicative of the total read depth at that location; note that it is significantly lower in regions where deletions were detected

Table 2 Numerical representation of filter progression
Table 3 Location and composition of all candidate microsatellites. More detailed information on candidate loci can be found in Supplemental Table S15


Ammonia is a common metabolic waste product in cell cultures. The accumulation of ammonia most often leads to decreased cell and recombinant protein productivity. Typical fed-batch cultures last for 14 to 20 days, where in the exponential phase, cell division can occur daily. De novo mutations that occur early in culture will be amplified and have the potential to dominate the cell population as the culture approaches harvest. In this study, two ammonia stresses (10 mM) and (30 mM) were used to investigate the genotoxic effects of this byproduct on CHO cell fed-batch cultures. Further, the role of ammonia stress on genome instability was investigated. Despite a relatively short exposure duration of 72 h, MSI loci were identified, which have the potential to be biomarkers for genome instability.

Metabolic response

The VCD, cell viability, and metabolic profiles indicated that the 30 mM ammonia stress significantly impacted the culture health, as the characteristic cell growth and metabolic profiles were significantly different from the control cultures. The effects of the 10 mM ammonia stress were less profound, yet the metabolic profiles and protein productivity were more sensitive to these changes than the cell viability and VCD profiles. At the time of sampling for genetic analysis (3.5 days), VCDs for the control and 10 mM cultures appeared to be matched, whereas the 30 mM ammonia-stressed cultures had lower VCDs. The decreased consumption of alanine observed for the ammonia-stressed cultures was the only metabolic difference observed at sampling for the whole genome sequencing. Alanine metabolic changes are known to occur under ammonia stress [39]. Therefore, the whole genome sequencing would identify changes due to the ammonia stresses, and not due to potential other culture condition differences that might accumulate.

Genome instability

Until now, efforts to characterize ammonia stress effects on CHO cells have mainly focused on transcriptome, proteome, and product characteristic changes [14, 16, 17, 40,41,42,43,44]. In this study, the effects of ammonia stress were further characterized by examining variants within functional genes and microsatellites. Whole genome sequencing allowed for variant SNPs and indels to be identified. Moreover, greater than 2300 high or moderate-impact novel gene variants were identified from the ammonia-stressed cultures that may impact cellular functions of critical pathways. KEGG and GO enrichment analyses confirmed that many of the variant genes affected pathways could lead to suboptimal clone performance. Though thousands of variant genes were identified, this list was narrowed to focus on genes pertaining to pathways involved in genome stability (Figs. 3 and 4, Table 1).

Alterations in critical genes responsible for a wide variety of processes such as transcription regulation, cell cycle regulation, tumor suppression, and signaling pathways may lead to global genome instability. De novo genomic SNPs and indels accumulating is typically the result of replication errors which can result from a variety of mechanisms such as replication stalling [45], replication fork collapse [46, 47], double-strand breaks [48, 49], environmental stressors, transcription regulation errors, or other replication errors [50]. All of these replication mechanisms can be linked to error correction fidelity of DNA repair mechanisms. These DNA repair errors, in turn, can lead to an accelerated variant accumulation rates and loss of genome stability [51]. Mutations, such as synonymous base changes in coding and regulatory regions, normally have little to no effect on gene transcription and translation, however, non-synonymous changes can have functional effects on the subsequent amino acid sequence and folding or function that ultimately can be linked to loss in cell viability.

Through text mining approaches, 273 of the variant genes found in the CHO genome were linked to human orthologs; whose function are related to genome stability maintenance. One gene identified is exceptionally well-known for its role in double-strand break repair and tumor suppression, Brca1 (Table 1); loss of Brca1 function has been associated with increased breast cancer incidence and metastasis, which demonstrates its critical function in maintaining stability [52]. Genome instability can be further exacerbated by the loss of tumor suppressor function. For example, Lin9 (Table 1) is a tumor suppressor that inhibits DNA synthesis and acts synergistically with the well-knownRb1 gene to prevent rapid, uncontrolled cell division [53]. Therefore, loss of function in Lin9 can lead to cancer-like growth of mutant cells that would eventually dominate the culture population.

Some variant genes belonged to three significantly enriched KEGG pathways related to genome instability in humans - cellular senescence, cell cycle, and homologous recombination (Fig. 2). Cellular senescence occurs as a result of multiple stimuli such as DNA damage and oxidative stress [54]. By forcing the cells into a non-replicative state, senescence can severely limit the productivity of cell culture, especially when it occurs before or in the early exponential growth phase. The 30 mM stressed cultures had more genes enriched in the senescence pathway (Supplemental Table S6), which makes variant genes in this pathway a likely contributor to the poor growth observed. The second pathway, cell cycle, was observed to have significant enrichment in union genes of the 10 mM and 30 mM ammonia-stressed cultures. The cell cycle contains multiple checkpoints to ensure daughter cells are healthy and contain undamaged DNA [55, 56]. Significant enrichment in this pathway indicates that damaged or otherwise improperly replicated DNA could be passed on to daughter cells. Finally, the homologous recombination pathway repairs damage caused by double strand breaks by using an identical sequence as a template [57]. This repair method is much more accurate than non-homologous end joining and is less prone to variant generation [58].

It should be noted that while the mismatch repair (MMR) pathway was not found to be significantly enriched, three notable MMR genes accumulated variants: Mlh3 (a MutL homolog), Rpa1, and Abl1 (Supplemental Table S5). An impaired or inefficient MMR system can lead to the accumulation of mutations in functional genes over cell divisions that are critical to the cell’s survival and can lead to loss of genetic stability [59] or disease states, such as cancer [60]. The need for a highly conserved MMR system can be observed by the presence of multiple orthologs of MutS and MutL in eukaryotic genomes [61]. MutS binds to base mismatches or small indels [62, 63] while MutL is responsible for communicating the identification of mismatch events to downstream elements of MMR such as exonucleases [64].

Microsatellite instability

Variants in genes that regulate the MMR pathway may be an origin to the cascade of events that leads to genome instability. When the MMR pathway in a cell is compromised, mistakes can occur and propagate indiscriminately across the genome as cell division occurs [61,62,63]. Unfaithful replication of genomic repeats, such as microsatellite repeats, have been used as effective biomarkers in predicting certain diseases, such as cancer In this study, we found 1022 microsatellites with variable repeat lengths in the WGS reads from ammonia stressed samples. We developed criteria and an approach to identify microsatellite loci that have variant length that can be attributed to ammonia stress. For each site, we considered read depth, mutation type, and frequency in affected and control samples to subset 124 candidate MSI loci that contain indels with a dose dependent response to the ammonia concentrations that were not observed in the control cultures. This set of 124 MSI loci represent a potential biomarker set that could have utility to predict genome instability in CHO cell cultures under stressful culture conditions.


The accumulation of metabolic wastes, such as ammonia, can have a profound effect on CHO cell culture viability, transcriptome, and recombinant protein productivity. Additionally, past work, as well as this study, have observed shifts in growth patterns and metabolic profiles due to the ammonia stress. Further, in this study, it was observed that high levels of exogenous ammonia caused de novo mutations, such as SNPs and indels, within functional genes. More importantly, these mutations persisted throughout the culture population. Variants were identified in the genes that regulate critical cellular processes, such as DNA repair; which is a hallmark of genome instability. In addition to characterizing the microsatellite content of the Chinese hamster genome, potential MSI loci that exhibited unfaithful replication in the presence of exogenous ammonia were identified; these microsatellites could be utilized as a tool to diagnose genome instability in future work.

Availability of data and materials

The genomic sequence data generated and/or analyzed during the current study are available in the NCBI short read archive under BioProject: PRJNA579347 Submission ID: SAMN13108689. All other data generated or analyzed during this study are included in this published article and its supplementary information files.


  1. 1.

    Carlson R. Biodesic 2011 bioeconomy update. In: BioDesic; 2011.

    Google Scholar 

  2. 2.

    Walsh G. Biopharmaceutical benchmarks 2014. Nat Biotechnol. 2014;32(10):992–1000.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    The Pharmaceutical and Biotech Industries in the United States [].

  4. 4.

    Walsh G. Post-translational modifications of protein biopharmaceuticals. Drug Discov Today. 2010;15(17–18):773–80.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Betts Z, Dickson AJ. Assessment of UCOE on recombinant EPO production and expression stability in amplified Chinese hamster ovary cells. Mol Biotechnol. 2015;57(9):846–58.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Dahodwala H, Lee KH. The fickle CHO: a review of the causes, implications, and potential alleviation of the CHO cell line instability problem. Curr Opin Biotechnol. 2019;60:128–37.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Kim M, O'Callaghan PM, Droms KA, James DC. A mechanistic understanding of production instability in CHO cell lines expressing recombinant monoclonal antibodies. Biotechnol Bioeng. 2011;108(10):2434–46.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Veith N, Ziehr H, MacLeod RA, Reamon-Buettner SM. Mechanisms underlying epigenetic and transcriptional heterogeneity in Chinese hamster ovary (CHO) cell lines. BMC Biotechnol. 2016;16:6.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  9. 9.

    Feichtinger J, Hernandez I, Fischer C, Hanscho M, Auer N, Hackl M, Jadhav V, Baumann M, Krempl PM, Schmidl C, et al. Comprehensive genome and epigenome characterization of CHO cells in response to evolutionary pressures and over time. Biotechnol Bioeng. 2016;113(10):2241–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Kaas CS, Kristensen C, Betenbaugh MJ, Andersen MR. Sequencing the CHO DXB11 genome reveals regional variations in genomic stability and haploidy. BMC Genomics. 2015;16:160.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. 11.

    Lewis NE, Liu X, Li Y, Nagarajan H, Yerganian G, O'Brien E, Bordbar A, Roth AM, Rosenbloom J, Bian C, et al. Genomic landscapes of Chinese hamster ovary cell lines as revealed by the Cricetulus griseus draft genome. Nat Biotechnol. 2013;31(8):759–65.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Arnold L, Lee K, Rucker-Pezzini J, Lee JH. Implementation of fully integrated continuous antibody processing: effects on productivity and COGm. Biotechnol J. 2019;14(2):e1800061.

    PubMed  Article  CAS  Google Scholar 

  13. 13.

    Bandyopadhyay AA, O’Brien SA, Zhao L, Fu H-Y, Vishwanathan N, Hu W-S. Recurring genomic structural variation leads to clonal instability and loss of productivity. Biotechnol Bioeng. 2019;116(1):41–53.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Ha TK, Kim YG, Lee GM. Understanding of altered N-glycosylation-related gene expression in recombinant Chinese hamster ovary cells subjected to elevated ammonium concentration by digital mRNA counting. Biotechnol Bioeng. 2015;112(8):1583–93.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Mulukutla BC, Gramer M, Hu WS. On metabolic shift to lactate consumption in fed-batch culture of mammalian cells. Metab Eng. 2012;14(2):138–49.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Kim DY, Chaudhry MA, Kennard ML, Jardon MA, Braasch K, Dionne B, Butler M, Piret JM. Fed-batch CHO cell t-PA production and feed glutamine replacement to reduce ammonia production. Biotechnol Prog. 2013;29(1):165–75.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Pereira AGM, Walther JL, Hollenbach M, Young JD. C-13 Flux Analysis Reveals that Rebalancing Medium Amino Acid Composition can Reduce Ammonia Production while Preserving Central Carbon Metabolism of CHO Cell Cultures. Biotechnol J. 2018;13(10):1700518.

    Article  CAS  Google Scholar 

  18. 18.

    Peltomaki P. Role of DNA mismatch repair defects in the pathogenesis of human cancer. J Clin Oncol. 2003;21(6):1174–9.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Preston BD, Albertson TM, Herr AJ. DNA replication fidelity and cancer. Semin Cancer Biol. 2010;20(5):281–93.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Muzny DM, Bainbridge MN, Chang K, Dinh HH, Drummond JA, Fowler G, Kovar CL, Lewis LR, Morgan MB, Newsham IF, et al. Comprehensive molecular characterization of human colon and rectal cancer. Nature. 2012;487(7407):330–7.

    CAS  Article  Google Scholar 

  21. 21.

    Hause RJ, Pritchard CC, Shendure J, Salipante SJ. Classification and characterization of microsatellite instability across 18 cancer types. Nat Med. 2016;22(11):1342–50.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Wu X, Xu Y, Chai W, Her C. Causal link between microsatellite instability and hMRE11 dysfunction in human cancers. Mol Cancer Res. 2011;9(11):1443–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Vilar E, Gruber SB. Microsatellite instability in colorectal cancer-the stable evidence. Nat Rev Clin Oncol. 2010;7(3):153–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Bass AJ, Thorsson V, Shmulevich I, Reynolds SM, Miller M, Bernard B, Hinoue T, Laird PW, Curtis C, Shen H, et al. Comprehensive molecular characterization of gastric adenocarcinoma. Nature. 2014;513(7517):202–9.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  25. 25.

    Cancer Genome Atlas Research N, Kandoth C, Schultz N, Cherniack AD, Akbani R, Liu Y, Shen H, Robertson AG, Pashtan I, Shen R, et al. Integrated genomic characterization of endometrial carcinoma. Nature. 2013;497(7447):67–73.

    Article  CAS  Google Scholar 

  26. 26.

    Boland CR, Thibodeau SN, Hamilton SR, Sidransky D, Eshleman JR, Burt RW, Meltzer SJ, Rodriguez-Bigas MA, Fodde R, Ranzani GN, et al. A National Cancer Institute workshop on microsatellite instability for cancer detection and familial predisposition: development of international criteria for the determination of microsatellite instability in colorectal cancer. Cancer Res. 1998;58(22):5248–57.

    CAS  PubMed  Google Scholar 

  27. 27.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome project data processing S. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. 30.

    Poplin R, Ruano-Rubio V, DePristo MA, Fennell TJ, Carneiro MO, Van der Auwera GA, Kling DE, Gauthier LD, Levy-Moonshine A, Roazen D, et al. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv. 2018:201178. [Preprints].

  31. 31.

    Cingolani P, Platts A, Wangle L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.

    CAS  Article  Google Scholar 

  32. 32.

    Thiel T, Michalek W, Varshney RK, Graner A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003;106(3):411–22.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Wei CH, Kao HY, Lu Z. PubTator: a web-based text mining tool for assisting biocuration. Nucleic Acids Res. 2013;41(Web Server issue):W518–22.

    PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Yu G. clusterProfiler: An universal enrichment tool for functional and comparative study. bioRxiv. 2018:256784. [Preprints].

  37. 37.

    Elliott K, Anderson J, Gavin C, Blakeman K, Harcum S, Harris G. Spent media analysis with an integrated CE-MS analyzer of Chinese hamster ovary cells grown in an ammonia-stressed parallel microbioreactor platform. BioProcess J. 2020. p. 19.

  38. 38.

    Fan Y, Del Val IJ, Muller C, Sen JW, Rasmussen SK, Kontoravdi C, Weilguny D, Andersen MR. Amino acid and glucose metabolism in fed-batch CHO cell culture affects antibody production and glycosylation. Biotechnol Bioeng. 2015;112(3):521–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Pereira S, Kildegaard HF, Andersen MR. Impact of CHO Metabolism on Cell Growth and Protein Production: An Overview of Toxic and Inhibiting Metabolites and Nutrients. Biotechnol J. 2018;13(3):e1700499.

    PubMed  Article  CAS  Google Scholar 

  40. 40.

    Brodsky AN, Caldwell M, Bae S, Harcum SW. Glycosylation-related genes in NS0 cells are insensitive to moderately elevated ammonium concentrations. J Biotechnol. 2014;187:78–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Chen P, Harcum SW. Effects of elevated ammonium on glycosylation gene expression in CHO cells. Metab Eng. 2006;8(2):123–32.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Freund NW, Croughan MS. A Simple Method to Reduce both Lactic Acid and Ammonium Production in Industrial Animal Cell Culture. Int J Mol Sci. 2018;19(2):385.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  43. 43.

    Yang M, Butler M. Effects of Ammonia and glucosamine on the heterogeneity of erythropoietin Glycoforms. Biotechnol Prog. 2002;18(1):129–38.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Genzel Y, Ritter JB, Konig S, Alt R, Reichl U. Substitution of glutamine by pyruvate to reduce ammonia formation and growth inhibition of mammalian cells. Biotechnol Prog. 2005;21(1):58–69.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Macheret M, Halazonetis TD. DNA replication stress as a hallmark of cancer. Annu Rev Pathol. 2015;10:425–48.

    CAS  Article  Google Scholar 

  46. 46.

    Zeman MK, Cimprich KA. Causes and consequences of replication stress. Nat Cell Biol. 2014;16(1):2–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Petermann E, Orta ML, Issaeva N, Schultz N, Helleday T. Hydroxyurea-stalled replication forks become progressively inactivated and require two different RAD51-mediated pathways for restart and repair. Mol Cell. 2010;37(4):492–502.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Bartkova J, Rezaei N, Liontos M, Karakaidos P, Kletsas D, Issaeva N, Vassiliou LV, Kolettas E, Niforou K, Zoumpourlis VC, et al. Oncogene-induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints. Nature. 2006;444(7119):633–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Di Micco R, Fumagalli M, Cicalese A, Piccinin S, Gasparini P, Luise C, Schurra C, Garre M, Nuciforo PG, Bensimon A, et al. Oncogene-induced senescence is a DNA damage response triggered by DNA hyper-replication. Nature. 2006;444(7119):638–42.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  50. 50.

    Saitoh T, Hiraga S. Studies on Molecular Mechanism of DNA-Replication in Escherichia-Coli .3. Genetic-Analysis of Mutation Causing Resumption of DNA-Replication Sensitive to Rifampicin, Which Exists in DNA Mutant Defective in Initiation of Replication. Jpn J Genet. 1975;49(5):320.

    Google Scholar 

  51. 51.

    Ponder RG, Fonville NC, Rosenberg SM. A switch from high-fidelity to error-prone DNA double-strand break repair underlies stress-induced mutation. Mol Cell. 2005;19(6):791–804.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Motoyama N, Naka K. DNA damage tumor suppressor genes and genomic instability. Curr Opin Genet Dev. 2004;14(1):11–6.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Stein TI, Nudel R, Lieder I, Mazor Y, et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinformatics. 2016;54:1 30 31.

    Article  Google Scholar 

  54. 54.

    van Deursen JM. The role of senescent cells in ageing. Nature. 2014;509(7501):439–46.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  55. 55.

    Nakanishi M, Shimada M, Niida H. Genetic instability in cancer cells by impaired cell cycle checkpoints. Cancer Sci. 2006;97(10):984–9.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Houtgraaf JH, Versmissen J, van der Giessen WJ. A concise review of DNA damage checkpoints and repair in mammalian cells. Cardiovasc Revasc Med. 2006;7(3):165–72.

    PubMed  Article  Google Scholar 

  57. 57.

    Thompson LH, Schild D. Homologous recombinational repair of DNA ensures mammalian chromosome stability. Mutat Res. 2001;477(1–2):131–53.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Rodgers K, McVey M. Error-prone repair of DNA double-Strand breaks. J Cell Physiol. 2016;231(1):15–24.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Harfe BD, Jinks-Robertson S. DNA mismatch repair and genetic instability. Annu Rev Genet. 2000;34:359–99.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Peltomaki P. DNA mismatch repair gene mutations in human cancer. Environ Health Perspect. 1997;105(Suppl 4):775–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Buermeyer AB, Deschenes SM, Baker SM, Liskay RM. Mammalian DNA mismatch repair. Annu Rev Genet. 1999;33:533–64.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Modrich P. Mechanisms and biological effects of mismatch repair. Annu Rev Genet. 1991;25:229–53.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Modrich P, Lahue R. Mismatch repair in replication fidelity, genetic recombination, and cancer biology. Annu Rev Biochem. 1996;65:101–33.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Sancar A, Hearst JE. Molecular matchmakers. Science. 1993;259(5100):1415–20.

    CAS  PubMed  Article  Google Scholar 

Download references


This work was supported in part by the National Science Foundation [IIP-1624641 and CBET-1218345]; and the Advanced Mammalian Biomanufacturing Innovation Center (AMBIC) industrial membership fees.

Author information




CS and SH conceived and designed the study. DC performed analysis of microsatellite data and was the primary contributor in writing the manuscript; QW performed KEGG and GO term enrichment analyses, identified human-Chinese hamster orthologs, and contributed greatly to writing the manuscript; KE led CHO cell culturing; AB assisted with constructing the microsatellite filtering criteria and DNA extraction and library prep; DJ assisted with DNA extraction and library prep, ZL led the DNA extraction and library prep, CW assisted with manuscript review; SH assisted with manuscript review, CS processed raw genomic data and assisted with manuscript review. The author(s) read and approved the final manusrcipt.

Corresponding author

Correspondence to Christopher A. Saski.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Table S1.

Complete list of SNP variants identified in each treatment group. Supplemental Table S2. Complete list of indel variants identified in each treatment group and initial mutation score calculations. Supplemental Table S3. Variant reads identified within functional genes. Supplemental Table S4. Variants identified within protein coding genes. Supplemental Table S5.Human-Chinese hamster gene orthologs that can be linked to genome instability via text mining.

Additional file 2: Supplemental Table S6.

KEGG enrichment of variant genes that can be linked to genome instability. Supplemental Table S7. Statistically significant Biological Process GO terms of variant genes. Supplemental Table S8. Statistically significant Molecular Function GO terms of variant genes. Supplemental Table S9. Statistically significant Cellular Component GO terms of variant genes. Supplemental Table S10. Complete list of microsatelites found in the Chinese hamster genome. Supplemental Table S11. Microsatellites containing insertion and or deletion mutations. Supplemental Table S12. Indel loci with higher mutation frequencies in 30 mM ammonia-stressed cultures. Supplemental Table S13. Indel loci from table S12 where mutation frequency is higher in the 10 mM samples compared to the control. Supplemental Table S14. The intersection of dose-dependent indels (Table S13) and genome-wide microsatellites (S10). Supplemental Table S15. The 124 candidate microsatelite loci that exhibited dose-dependent variation in response to ammonia stress.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chitwood, D.G., Wang, Q., Elliott, K. et al. Characterization of metabolic responses, genetic variations, and microsatellite instability in ammonia-stressed CHO cells grown in fed-batch cultures. BMC Biotechnol 21, 4 (2021).

Download citation


  • CHO
  • Ammonia
  • MSI
  • Genome instability
  • Biomarker