Deep learning models incorporating endogenous factors beyond DNA sequences improve the prediction accuracy of base editing outcomes.

Adenine base editors (ABEs) and cytosine base editors (CBEs) enable the single nucleotide editing of targeted DNA sites avoiding generation of double strand breaks, however, the genomic features that influence the outcomes of base editing in vivo still remain to be characterized. High-throughput datasets from lentiviral integrated libraries were used to investigate the sequence features affecting base editing outcomes, but the effects of endogenous factors beyond the DNA sequences are still largely unknown. Here the base editing outcomes of ABE and CBE were evaluated in mammalian cells for 5012 endogenous genomic sites and 11,868 genome-integrated target sequences, with 4654 genomic sites sharing the same target sequences. The comparative analyses revealed that the editing outcomes of ABE and CBE at endogenous sites were substantially different from those obtained using genome-integrated sequences. We found that the base editing efficiency at endogenous target sites of both ABE and CBE was influenced by endogenous factors, including epigenetic modifications and transcriptional activity. A deep-learning algorithm referred as BE_Endo, was developed based on the endogenous factors and sequence information from our genomic datasets, and it yielded unprecedented accuracy in predicting the base editing outcomes. These findings along with the developed computational algorithms may facilitate future application of BEs for scientific research and clinical gene therapy.

Efficient CRISPR/Cas9-mediated gene editing in mammalian cells by the novel selectable traffic light reporters.

CRISPR/Cas9 is a powerful tool for gene editing in various cell types and organisms. However, it is still challenging to screen genetically modified cells from an excess of unmodified cells. Our previous studies demonstrated that surrogate reporters can be used for efficient screening of genetically modified cells. Here, we developed two novel traffic light screening reporters, puromycin-mCherry-EGFP (PMG) based on single-strand annealing (SSA) and homology-directed repair (HDR), respectively, to measure the nuclease cleavage activity within transfected cells and to select genetically modified cells. We found that the two reporters could be self-repaired coupling the genome editing events driven by different CRISPR/Cas nucleases, resulting in a functional puromycin-resistance and EGFP selection cassette that can be afforded to screen genetically modified cells by puromycin selection or FACS enrichment. We further compared the novel reporters with different traditional reporters at several endogenous loci in different cell lines, for the enrichment efficiencies of genetically modified cells. The results indicated that the SSA-PMG reporter exhibited improvements in enriching gene knockout cells, while the HDR-PMG system was very useful in enriching knock-in cells. These results provide robust and efficient surrogate reporters for the enrichment of CRISPR/Cas9-mediated editing in mammalian cells, thereby advancing basic and applied research.

Cytosine base editors induce off-target mutations and adverse phenotypic effects in transgenic mice.

Base editors have been reported to induce off-target mutations in cultured cells, mouse embryos and rice, but their long-term effects in vivo remain unknown. Here, we develop a Systematic evaluation Approach For gene Editing tools by Transgenic mIce (SAFETI), and evaluate the off-target effects of BE3, high fidelity version of CBE (YE1-BE3-FNLS) and ABE (ABE7.10F148A) in ~400 transgenic mice over 15 months. Whole-genome sequence analysis reveals BE3 expression generated de novo mutations in the offspring of transgenic mice. RNA-seq analysis reveals both BE3 and YE1-BE3-FNLS induce transcriptome-wide SNVs, and the numbers of RNA SNVs are positively correlated with CBE expression levels across various tissues. By contrast, ABE7.10F148A shows no detectable off-target DNA or RNA SNVs. Notably, we observe abnormal phenotypes including obesity and developmental delay in mice with permanent genomic BE3 overexpression during long-time monitoring, elucidating a potentially overlooked aspect of side effects of BE3 in vivo.

Optimization of C-to-G base editors with sequence context preference predictable by machine learning methods.

Efficient and precise base editors (BEs) for C-to-G transversion are highly desirable. However, the sequence context affecting editing outcome largely remains unclear. Here we report engineered C-to-G BEs of high efficiency and fidelity, with the sequence context predictable via machine-learning methods. By changing the species origin and relative position of uracil-DNA glycosylase and deaminase, together with codon optimization, we obtain optimized C-to-G BEs (OPTI-CGBEs) for efficient C-to-G transversion. The motif preference of OPTI-CGBEs for editing 100 endogenous sites is determined in HEK293T cells. Using a sgRNA library comprising 41,388 sequences, we develop a deep-learning model that accurately predicts the OPTI-CGBE editing outcome for targeted sites with specific sequence context. These OPTI-CGBEs are further shown to be capable of efficient base editing in mouse embryos for generating Tyr-edited offspring. Thus, these engineered CGBEs are useful for efficient and precise base editing, with outcome predictable based on sequence context of targeted sites.

A Universal Surrogate Reporter for Efficient Enrichment of CRISPR/Cas9-Mediated Homology-Directed Repair in Mammalian Cells.

CRISPR/Cas9-mediated homology-directed repair (HDR) can be leveraged to precisely engineer mammalian genomes. However, the inherently low efficiency of HDR often hampers to identify the desired modified cells. Here, we developed a novel universal surrogate reporter system that efficiently enriches for genetically modified cells arising from CRISPR/Cas9-induced HDR events (namely, the "HDR-USR" system). This episomally based reporter can be self-cleaved and self-repaired via HDR to create a functional puromycin selection cassette without compromising genome integrity. Co-transfection of the HDR-USR system into host cells and transient puromycin selection efficiently achieves enrichment of HDR-modified cells. We tested the system for precision point mutation at 16 loci in different human cell lines and one locus in two rodent cell lines. This system exhibited dramatic improvements in HDR efficiency at a single locus (up to 20.7-fold) and two loci at once (42% editing efficiency compared to zero in the control), as well as greatly improved knockin efficiency (8.9-fold) and biallelic deletion (35.9-fold) at test loci. Further increases were achieved by co-expression of yeast Rad52 and linear single-/double-stranded DNA donors. Taken together, our HDR-USR system provides a simple, robust and efficient surrogate reporter for the enrichment of CRISPR/Cas9-induced HDR-based precision genome editing across various targeting loci in different cell lines.

sgRNA-shRNA Structure Mediated SNP Site Editing on Porcine IGF2 Gene by CRISPR/StCas9.

The SNP within intron 3 of the porcine IGF2 gene (G3072A) plays an important role for muscle growth and fat deposition in pigs. In this study, the StCas9 derived from Streptococcus thermophilus together with the Drosha-mediated sgRNA-shRNA structure were combined to boost the G to A base editing on the IGF2 SNP site, which we called "SNP editing." The codon-humanized StCas9 as we previously reported was firstly compared with the prevalently used SpCas9 derived from Streptococcus pyogenes using our idiomatic surrogate report assay, and the StCas9 demonstrated a comparable targeting activity. On the other hand, by combining shRNA with sgRNA, simultaneous gene silencing and genome targeting can be achieved. Thus, the novel IGF2.sgRNA-LIG4.shRNA-IGF2.sgRNA structure was constructed to enhance the sgRNA/Cas9-mediated HDR-based IGF2 SNP editing by silencing the LIG4 gene, which is a key molecule of the HDR's competitive NHEJ pathway. The sgRNA-shRNA/StCas9 all-in-one expression vector and the IGF2.sgRNA/StCas9 as control were separately used to transfect porcine PK15 cells together with an ssODNs donor for the IGF2 SNP editing. The editing events were detected by the RFLP assay, Sanger sequencing as well as Deep-sequencing, and the Deep-sequencing results finally demonstrated a significant higher HDR-based editing efficiency (16.38%) for our sgRNA-shRNA/StCas9 strategy. In short, we achieved effective IGF2 SNP editing by using the combined sgRNA-shRNA/StCas9 strategy, which will facilitate the further production of base-edited animals and perhaps extend for the gene therapy for the base correction of some genetic diseases.

The Effects of Matrix Metalloproteinase-9 on Dairy Goat Mastitis and Cell Survival of Goat Mammary Epithelial Cells.

Matrix metalloproteinase-9 (MMP-9) is a zinc-dependent enzyme, and plays a crucial role in extracellular matrix degeneration, inflammation and tissue remodeling. However, the relationship between MMP-9 and somatic cell count (SCC) in goat milk and the role of MMP-9 in the regulation of mastitis are still unknown. In this study, we found MMP-9 was predominantly expressed in the spleen, intestine and mammary gland. The SCC in goat milk was positively correlated with MMP-9 expression, and staphylococcus aureus could markedly increase MMP-9 expression in goat mammary epithelial cells (GMEC) in dosage and time dependent manner. We also demonstrated that SB-3CT, an inhibitor of MMP-9, promoted apoptosis and inhibited proliferation in GMEC. Thus, MMP-9 may emerge as an easily measurable and sensitive parameter that reflects the number of somatic cells present in milk and a regulatory factor of apoptosis in GMEC.

Parathyroid hormone-related protein overexpression protects goat mammary gland epithelial cells from calcium-sensing receptor activation-induced apoptosis.

Normal mammary gland epithelial cells and breast cancer cells express the calcium-sensing receptor (CaSR), which is the master regulator of systemic calcium metabolism. During lactation, activation of the CaSR in mammary epithelial cells downregulates parathyroid hormone-related protein (PTHrP) levels in milk and in the circulation, and increases calcium transport into milk. However, very little information is available on the role of CaSR in goat mammary gland epithelial cells (GMECs) apoptosis. In this investigation, the full-length cDNA of CaSR from Xinong Saanen dairy goats was cloned, which contains an open-reading frame of 3,258 bp encoding 1,085 amino acids with a predicted molecular weight of 121.0 kDa and an isoelectric point of 5.65. The amino acid sequence is highly homologous with sheep, and the goat CaSR gene is mapped to chromosome 1. Quantitative real-time PCR suggested that CaSR was predominantly expressed in the heart, kidney and mammary gland. Then, we found the stimulation of CaSR with its activator gadolinium chloride (GdCl3) contributed to increase CaSR mRNA levels in GMECs and simultaneously promoted cell apoptosis, and these effects were abrogated partially by NPS2390 which is an inhibitor of CaSR. We also demonstrated that Ca(2+) increased CaSR mRNA levels and induced GMECs apoptosis and restrained cell proliferation. In contrast, PTHrP overexpression protected GMECs from calcium-induced apoptosis, and promoted cell proliferation. In conclusion, these results suggest that PTHrP overexpression protects GMECs from CaSR activation-induced apoptosis.

Signal transducer and activator of transcription 5a inhibited by pimozide may regulate survival of goat mammary gland epithelial cells by regulating parathyroid hormone-related protein.

The signal transducer and activator of transcription 5a (Stat5a) modulates genes involved in proliferation and survival and plays pivotal roles in regulating the function of the mammary gland during pregnancy, lactation, and involution. However, there is little information about the effects of Stat5a on apoptosis of goat mammary gland epithelial cells (GMECs). In addition, parathyroid hormone-related protein (PTHrP) is a key regulator in cellular calcium transport, mammary gland development and breast tumor biology. This study aimed to explore the interaction of Stat5a and PTHrP in GMEC apoptosis. Quantitative real time PCR (qRT-PCR) suggested that Stat5a was predominantly expressed in the mammary gland, lung, liver and spleen of goats. Treating the GMECs with pimozide, an inhibitor of Stat5a that decreases Stat5a tyrosine phosphorylation, increased PTHrP levels in GMECs in a dose-dependent manner and simultaneously promoted apoptosis of the GMECs. We also demonstrated that PTHrP inhibition induced GMEC apoptosis and restrained cell proliferation. In contrast, PTHrP overexpression protected GMECs from pimozide- and calcium-induced apoptosis, and promoted cell proliferation. Furthermore, pimozide and CaCl2 downregulated the antiapoptotic protein Bcl-2 mRNA expression, respectively, and these effects were protected by PTHrP overexpression. Interestingly, we also found that Stat5a suppressed the expression of matrix metalloproteinase 9 (MMP-9) which can induce goat mammary epithelial cell migration, but PTHrP increased MMP-9 mRNA level. Thus, Stat5a may regulate GMEC survival by regulating the expression of PTHrP.