Mouse Genome Research Articles

A Bcl11bN797K variant isolated from an immunodeficient patient inhibits early thymocyte development in mice.

BCL11B is a transcription factor with six C2H2-type zinc-finger domains. Studies in mice have shown that Bcl11b plays essential roles in T cell development. Several germline heterozygous BCL11B variants have been identified in human patients with inborn errors of immunity (IEI) patients. Among these, two de novo mis-sense variants cause asparagine (N) to lysine (K) replacement in distinct zinc-finger domains, BCL11BN441K and BCL11BN807K. To elucidate the pathogenesis of the BCL11BN807K variant, we generated a mouse model of BCL11BN807K by inserting the corresponding mutation, Bcl11bN797K, into the mouse genome. In Bcl11b+/N797K mice, the proportion of immature CD4-CD8+ single-positive thymocytes was increased, and the development of invariant natural killer cells was severely inhibited in a T-cell-intrinsic manner. Under competitive conditions, γδT cell development was outcompeted by control cells. Bcl11bN797K/N797K mice died within one day of birth. Recipient mice reconstituted with Bcl11bN797K/N797K fetal liver cells nearly lacked CD4+CD8+ double-positive thymocytes, which was consistent with the lack of their emergence in culture from Bcl11bN797K/N797K fetal liver progenitors. Interestingly, Bcl11bN797K/N797K progenitors gave rise to aberrant c-Kit+ and CD44+ cells both in vivo and in vitro. The increase in the proportion of immature CD8 single-positive thymocytes in the Bcl11bN797K mutants is caused, in part, by the inefficient activation of the Cd4 gene due to the attenuated function of the two Cd4 enhancers via distinct mechanisms. Therefore, we conclude that immunodeficient patient-derived Bcl11bN797K mutant mice elucidated a novel role for Bcl11b in driving the appropriate transition of CD4-CD8- into CD4+CD8+ thymocytes.

Candidate circulating microRNA biomarkers in dogs with chronic pancreatitis.

Pancreatitis is an important cause of disease and death in dogs. Available circulating biomarkers are not sufficiently sensitive and specific for a definitive diagnosis. Circulating microRNAs would be differentially expressed in dogs with chronic pancreatitis and could have potential as diagnostic biomarkers. Healthy controls (n = 19) and dogs with naturally occurring pancreatitis (n = 17). A retrospective case-control study. Dogs with pancreatitis were included if they satisfied diagnostic criteria for pancreatitis as adjudicated by 3 experts. MicroRNA was extracted from stored serum samples and sequenced. Reads were mapped to mature microRNA sequences in the canine, mouse, and human genomes. Differentially expressed microRNAs were identified and the potential mechanistic relevance explored using Qiagen Ingenuity Pathway Analysis (IPA). Reads mapping to 196 mature microRNA sequences were detected. Eight circulating microRNAs were significantly differentially expressed in dogs with pancreatitis (≥2-fold change and false discovery rate <0.05). Four of these mapped to the canine genome (cfa-miR-221, cfa-miR-222, cfa-miR-23a, and cfa-miR-205). Three mapped to the murine genome (mmu-miR-484, mmu-miR-6240, mmu-miR-101a-3p) and 1 to the human genome (hsa-miR-1290). Expression in dogs with pancreatitis was higher for 7 microRNAs and lower for mmu-miR-101a-3p. Qiagen IPA demonstrated a number of the differently expressed microRNAs are involved in a common pancreatic inflammatory pathway. The significantly differentially expressed microRNAs represent promising candidates for further validation as diagnostic biomarkers for canine pancreatitis.

A common allele increases endometrial Wnt4 expression, with antagonistic implications for pregnancy, reproductive cancers, and endometriosis

The common human SNP rs3820282 is associated with multiple phenotypes including gestational length and likelihood of endometriosis and cancer, presenting a paradigmatic pleiotropic variant. Deleterious pleiotropic mutations cause the co-occurrence of disorders either within individuals, or across population. When adverse and advantageous effects are combined, pleiotropy can maintain high population frequencies of deleterious alleles. To reveal the causal molecular mechanisms of this pleiotropic SNP, we introduced this substitution into the mouse genome by CRISPR/Cas 9. Previous work showed that rs3820282 introduces a high-affinity estrogen receptor alpha-binding site at the Wnt4 locus. Here, we show that this mutation upregulates Wnt4 transcription in endometrial stroma, following the preovulatory estrogen peak. Effects on uterine transcription include downregulation of epithelial proliferation and induction of progesterone-regulated pro-implantation genes. We propose that these changes increase uterine permissiveness to embryo invasion, whereas they decrease resistance to invasion by cancer and endometriotic foci in other estrogen-responsive tissues.

Gene replacement-Alzheimer's disease (GR-AD): Modeling the genetics of human dementias in mice.

Genetic studies conducted over the past four decades have provided us with a detailed catalog of genes that play critical roles in the etiology of Alzheimer's disease (AD) and related dementias (ADRDs). Despite this progress, as a field we have had only limited success in incorporating this rich complexity of human AD/ADRD genetics findings into our animal models of these diseases. Our primary goal for the gene replacement (GR)-AD project is to develop mouse lines that model the genetics of AD/ADRD as closely as possible. To do this, we are generating mouse lines in which the genes of interest are precisely and completely replaced in the mouse genome by their full human orthologs. Each model set consists of a control line with a wild-type human allele and variant lines that precisely match the human genomic sequence in the control line except for a high-impact pathogenic mutation or risk variant.

Donor template delivery by recombinant adeno-associated virus for the production of knock-in mice

BackgroundThe ability of recombinant adeno-associated virus to transduce preimplantation mouse embryos has led to the use of this delivery method for the production of genetically altered knock-in mice via CRISPR-Cas9. The potential exists for this method to simplify the production and extend the types of alleles that can be generated directly in the zygote, obviating the need for manipulations of the mouse genome via the embryonic stem cell route.ResultsWe present the production data from a total of 13 genetically altered knock-in mouse models generated using CRISPR-Cas9 electroporation of zygotes and delivery of donor repair templates via transduction with recombinant adeno-associated virus. We explore the efficiency of gene targeting at a total of 12 independent genetic loci and explore the effects of allele complexity and introduce strategies for efficient identification of founder animals. In addition, we investigate the reliability of germline transmission of the engineered allele from founder mice generated using this methodology. By comparing our production data against genetically altered knock-in mice generated via gene targeting in embryonic stem cells and their microinjection into blastocysts, we assess the animal cost of the two methods.ConclusionsOur results confirm that recombinant adeno-associated virus transduction of zygotes provides a robust and effective delivery route for donor templates for the production of knock-in mice, across a range of insertion sizes (0.9–4.7 kb). We find that the animal cost of this method is considerably less than generating knock-in models via embryonic stem cells and thus constitutes a considerable 3Rs reduction.

CRISPR-TE: a web-based tool to generate single guide RNAs targeting transposable elements

BackgroundThe CRISPR/Cas systems have emerged as powerful tools in genome engineering. Recent studies highlighting the crucial role of transposable elements (TEs) have stimulated research interest in manipulating these elements to understand their functions. However, designing single guide RNAs (sgRNAs) that are specific and efficient for TE manipulation is a significant challenge, given their sequence repetitiveness and high copy numbers. While various sgRNA design tools have been developed for gene editing, an optimized sgRNA designer for TE manipulation has yet to be established.ResultsWe present CRISPR-TE, a web-based application featuring an accessible graphical user interface, available at https://www.crisprte.cn/, and currently tailored to the human and mouse genomes. CRISPR-TE identifies all potential sgRNAs for TEs and provides a comprehensive solution for efficient TE targeting at both the single copy and subfamily levels. Our analysis shows that sgRNAs targeting TEs can more effectively target evolutionarily young TEs with conserved sequences at the subfamily level.ConclusionsCRISPR-TE offers a versatile framework for designing sgRNAs for TE targeting. CRISPR-TE is publicly accessible at https://www.crisprte.cn/ as an online web service and the source code of CRISPR-TE is available at https://github.com/WanluLiuLab/CRISPRTE/.

Genetic influences on motor learning and performance and superperforming mutants revealed by random mutational survey of the mouse genome.

Evolution depends upon genetic variations that influence physiology. As defined in a genetic screen, phenotypic performance may be enhanced or degraded by such mutations. We set out to detect mutations that influence motor function, including motor learning. Thus, we tested the motor effects of 36,444 non-synonymous coding/splicing mutations induced in the germline of C57BL/6J mice with N-ethyl-N-nitrosourea by measuring changes in the performance of repetitive rotarod trials while blinded to genotype. Automated meiotic mapping was used to implicate individual mutations in causation. In total, 32,726 mice bearing all the variant alleles were screened. This was complemented with the simultaneous testing of 1408 normal mice for reference. In total, 16.3% of autosomal genes were thus rendered detectably hypomorphic or nullified by mutations in homozygosity and motor tested in at least three mice. This approach allowed us to identify superperformance mutations in Rif1, Tk1, Fan1 and Mn1. These genes are primarily related, among other less well-characterized functions, to nucleic acid biology. We also associated distinct motor learning patterns with groups of functionally related genes. These functional sets included, preferentially, histone H3 methyltransferase activity for mice that learnt at an accelerated rate relative to the remaining mutant mice. The results allow for an estimation of the fraction of mutations that can modify a behaviour influential for evolution such as locomotion. They may also enable, once the loci are further validated and the mechanisms elucidated, the harnessing of the activity of the newly identified genes to enhance motor ability or to counterbalance disability or disease. KEY POINTS: We studied the effect of chemically induced random mutations on mouse motor performance. An array of mutations influenced the rate of motor learning. DNA regulation genes predominated among these mutant loci. Several mutations in unsuspected genes led to superperformance. Assuming little-biased mutagenicity, the results allow for an estimation of the probability for any spontaneous mutation to influence a behaviour such as motor learning and ultimate performance.

An analysis pipeline for understanding 6-thioguanine effects on a mouse tumour genome

Mouse tumour models are extensively used as a pre-clinical research tool in the field of oncology, playing an important role in anticancer drugs discovery. Accordingly, in cancer genomics research, the demand for next-generation sequencing (NGS) is increasing, and consequently, the need for data analysis pipelines is likewise growing. Most NGS data analysis solutions to date do not support mouse data or require highly specific configuration for their use. Here, we present a genome analysis pipeline for mouse tumour NGS data including the whole-genome sequence (WGS) data analysis flow for somatic variant discovery, and the RNA-seq data flow for differential expression, functional analysis and neoantigen prediction. The pipeline is based on standards and best practices and integrates mouse genome references and annotations. In a recent study, the pipeline was applied to demonstrate the efficacy of low dose 6-thioguanine (6TG) treatment on low-mutation melanoma in a pre-clinical mouse model. Here, we further this study and describe in detail the pipeline and the results obtained in terms of tumour mutational burden (TMB) and number of predicted neoantigens, and correlate these with 6TG effects on tumour volume. Our pipeline was expanded to include a neoantigen analysis, resulting in neopeptide prediction and MHC class I antigen presentation evaluation. We observed that the number of predicted neoepitopes were more accurate indicators of tumour immune control than TMB. In conclusion, this study demonstrates the usability of the proposed pipeline, and suggests it could be an essential robust genome analysis platform for future mouse genomic analysis.

Comparative Gene Signature of Nociceptors Innervating Mouse Molar Teeth, Cranial Meninges, and Cornea.

The trigeminal ganglion (TG) collects afferent sensory information from various tissues. Recent large-scale RNA sequencing of neurons of the TG and dorsal root ganglion has revealed a variety of functionally distinct neuronal subpopulations, but organ-specific information is lacking. To link transcriptomic and tissue-specific information, we labeled small-diameter neurons of 3 specific subpopulations of the TG by local application of lipophilic carbocyanine dyes to their innervation site in the dental pulp, cornea, and meninges (dura mater). We then collected mRNA-sequencing data from fluorescent neurons. Differentially expressed genes (DEGs) were analyzed and subjected to downstream gene set enrichment analysis (GSEA), and ion channel profiling was performed. A total of 10,903 genes were mapped to the mouse genome (>500 reads). DEG analysis revealed 18 and 81 genes with differential expression (log 2 fold change > 2, Padj < .05) in primary afferent neurons innervating the dental pulp (dental primary afferent neurons [DPAN]) compared to those innervating the meninges (meningeal primary afferent neurons [MPAN]) and the cornea (corneal primary afferent neurons [CPAN]). We found 250 and 292 genes differentially expressed in MPAN as compared to DPAN and to CPAN, and 21 and 12 in CPAN as compared to DPAN and MPAN. Scn2b had the highest log 2 fold change when comparing DPAN versus MPAN and Mmp12 was the most prominent DEG when comparing DPAN versus CPAN and, CPAN versus MPAN. GSEA revealed genes of the immune and mitochondrial oxidative phosphorylation system for the DPAN versus MPAN comparison, cilium- and ribosome-related genes for the CPAN versus DPAN comparison, and respirasome, immune cell- and ribosome-related gene sets for the CPAN versus MPAN comparison. DEG analysis for ion channels revealed no significant differences between the neurons set except for the sodium voltage-gated channel beta subunit 2, Scn2b . However, in each tissue a few ion channels turned up with robust number of reads. In DPAN, these were Cacna1b , Trpv2 , Cnga4 , Hcn1 , and Hcn3 , in CPAN Trpa1 , Trpv1 , Cacna1a , and Kcnk13 and in MPAN Trpv2 and Scn11a . Our study uncovers previously unknown differences in gene expression between sensory neuron subpopulations from the dental pulp, cornea, and dura mater and provides the basis for functional studies, including the investigation of ion channel function and their suitability as targets for tissue-specific analgesia.

An ancient polymorphic regulatory region within the BDNF gene associated with obesity modulates anxiety-like behaviour in mice and humans

Obesity and anxiety are morbidities notable for their increased impact on society during the recent COVID-19 pandemic. Understanding the mechanisms governing susceptibility to these conditions will increase our quality of life and resilience to future pandemics. In the current study, we explored the function of a highly conserved regulatory region (BE5.1) within the BDNF gene that harbours a polymorphism strongly associated with obesity (rs10767664; p = 4.69 × 10–26). Analysis in primary cells suggested that the major T-allele of BE5.1 was an enhancer, whereas the obesity-associated A-allele was not. However, CRISPR/CAS9 deletion of BE5.1 from the mouse genome (BE5.1KO) produced no significant effect on the expression of BDNF transcripts in the hypothalamus, no change in weight gain after 28 days and only a marginally significant increase in food intake. Nevertheless, transcripts were significantly increased in the amygdala of female mice and elevated zero maze and marble-burying tests demonstrated a significant increase in anxiety-like behaviour that could be reversed by diazepam. Consistent with these observations, human GWAS cohort analysis demonstrated a significant association between rs10767664 and anxiousness in human populations. Intriguingly, interrogation of the human GTEx eQTL database demonstrated no effect on BDNF mRNA levels associated with rs10767664 but a highly significant effect on BDNF-antisense (BDNF-AS) gene expression and splicing. The subsequent observation that deletion of BE5.1 also significantly reduced BDNF-AS expression in mice suggests a novel mechanism in the regulation of BDNF expression common to mice and humans, which contributes to the modulation of mood and anxiety in both species.

Abstract C034: Genome-wide CRISPR screen identifies novel immunotherapeutic targets in pancreatic cancer

Abstract Background: Cancer-cell-intrinsic properties dictate the immune contexture of tumors, impacting the outcome of immunotherapy. Despite the increasing value of CRISPR/Cas9-based loss-of-function (LOF) genetic screens for systematically assessing gene knockout phenotype in tumor cells, it still present challenges, particularly when it comes to genomic editing in the context of large scale experiments and the identification of immunomodulatory genes in in vivo application. To address this, we have established an AsCpf1 (Cas12a)-based screen platform to enable comprehensive whole genome in vivo LOF screening and uncover novel immunotherapeutic targets for pancreatic cancer, which is notoriously resistant to current immunotherapies. Methods: We designed and constructed the optimized minimal CRISPR/AsCpf1 library that targets the entire mouse genome and is only one-fourth the size of the existing genome-wide Cas9-based library. To identify novel immune-related gene vulnerabilities in pancreatic cancer, we infected cell cultures derived from the classic p48-Cre KrasLSL-G12D/+ p53+/- genetically engineered pancreatic cancer mouse model with our AsCpf1 library for in vivo LOF screen in immune-competent and -incompetent mice. Next generation sequencing analysis was performed with the resulting tumors. The potential immunomodulatory hits were called out through BAGEL analysis, using a false discovery rate (FDR) &amp;lt; 0.05, and by comparing between the two mouse groups. Results: Our AsCpf1-based screen identified highly confident hits with potential immunomodulation properties. Notably, we observed depletions of various components of the MHC complex, including H13, specifically in the immune-competent background, implicating their requirement for the regulation of anti-tumor immunity. Depletion of H13 with CRISPR/Cas9 system significantly suppressed pancreatic tumor growth in T-cell dependent manner. Further immune profiling revealed that the loss of H13 triggered potent immune response which was featured by a substantial increase in effective CD8+ T cells and reduction in immunosuppressive macrophages in both subcutaneous and orthotopic models of pancreatic cancer. These findings are in accordance with our clinically relevant analysis showing that low expression of H13 is correlated with improved patient survival in pancreatic cancer and augmented responsiveness to immune checkpoint therapies in several cancer types. Conclusions: Our AsCpf1-based LOF screening platform enables unbiased exploration of in vivo functional players at whole-genome level, which can be applied to diverse disease context and models. In our lead screen, we have identified H13 as a promising immunomodulatory target for patients with pancreatic cancer. Citation Format: Xiaofei Wang, Jintan Liu, Giulio Draetta, Haoqiang Ying, Wantong Yao. Genome-wide CRISPR screen identifies novel immunotherapeutic targets in pancreatic cancer [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Pancreatic Cancer; 2023 Sep 27-30; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(2 Suppl):Abstract nr C034.

Circular RNA Cdr1as inhibits proliferation and delays injury-induced regeneration of the intestinal epithelium.

Circular RNAs (circRNAs) are highly expressed in the mammalian intestinal epithelium, but their functions remain largely unknown. Here, we identified the circRNA Cdr1as as a repressor of intestinal epithelial regeneration and defense. Cdr1as levels increased in mouse intestinal mucosa after colitis and septic stress, as well as in human intestinal mucosa from patients with inflammatory bowel disease and sepsis. Ablation of the Cdr1as locus from the mouse genome enhanced renewal of the intestinal mucosa, promoted injury-induced epithelial regeneration, and protected the mucosa against colitis. We found approximately 40 microRNAs, including miR-195, differentially expressed between intestinal mucosa of Cdr1as-knockout (Cdr1as-/-) versus littermate mice. Increasing the levels of Cdr1as inhibited intestinal epithelial repair after wounding in cultured cells and repressed growth of intestinal organoids cultured ex vivo, but this inhibition was abolished by miR-195 silencing. The reduction in miR-195 levels in the Cdr1as-/- intestinal epithelium was the result of reduced stability and processing of the precursor miR-195. These findings indicate that Cdr1as reduces proliferation and repair of the intestinal epithelium at least in part via interaction with miR-195 and highlight a role for induced Cdr1as in the pathogenesis of unhealed wounds and disrupted renewal of the intestinal mucosa.

CRISPR-Cas9 base editing of pathogenic CaMKIIδ improves cardiac function in a humanized mouse model.

Cardiovascular diseases are the most common cause of worldwide morbidity and mortality, highlighting the necessity for advanced therapeutic strategies. Ca2+/calmodulin-dependent protein kinase IIδ (CaMKIIδ) is a prominent inducer of various cardiac disorders, which is mediated by 2 oxidation-sensitive methionine residues within the regulatory domain. We have previously shown that ablation of CaMKIIδ oxidation by CRISPR-Cas9 base editing enables the heart to recover function from otherwise severe damage following ischemia/reperfusion (IR) injury. Here, we extended this therapeutic concept toward potential clinical translation. We generated a humanized CAMK2D knockin mouse model in which the genomic sequence encoding the entire regulatory domain was replaced with the human sequence. This enabled comparison and optimization of two different editing strategies for the human genome in mice. To edit CAMK2D in vivo, we packaged the optimized editing components into an engineered myotropic adeno-associated virus (MyoAAV 2A), which enabled efficient delivery at a very low AAV dose into the humanized mice at the time of IR injury. CAMK2D-edited mice recovered cardiac function, showed improved exercise performance, and were protected from myocardial fibrosis, which was otherwise observed in injured control mice after IR. Our findings identify a potentially effective strategy for cardioprotection in response to oxidative damage.

Myotis bat STING attenuates aging-related inflammation in female mice.

Bats, notable as the only flying mammals, serve as natural reservoir hosts for various highly pathogenic viruses in humans (e.g., SARS-CoV and Ebola virus). Furthermore, bats exhibit an unparalleled longevity among mammals relative to their size, particularly the Myotis bats, which can live up to 40 years. However, the mechanisms underlying these distinctive traits remain incompletely understood. In our prior research, we demonstrated that bats exhibit dampened STING-interferon activation, potentially conferring upon them the capacity to mitigate virus- or aging-induced inflammation. To substantiate this hypothesis, we established the first in vivo bat-mouse model for aging studies by integrating Myotis davidii bat STING ( MdSTING) into the mouse genome. We monitored the genotypes of these mice and performed a longitudinal comparative transcriptomic analysis on MdSTING and wild-type mice over a 3-year aging process. Blood transcriptomic analysis indicated a reduction in aging-related inflammation in female MdSTING mice, as evidenced by significantly lower levels of pro-inflammatory cytokines and chemokines, immunopathology, and neutrophil recruitment in aged female MdSTING mice compared to aged wild-type mice in vivo. These results indicated that MdSTING knock-in attenuates the aging-related inflammatory response and may also improve the healthspan in mice in a sex-dependent manner. Although the underlying mechanism awaits further study, this research has critical implications for bat longevity research, potentially contributing to our comprehension of healthy aging in humans.

Enhancer Arrays Regulating Developmental Genes: Sox2 Enhancers as a Paradigm.

Enhancers are the primary regulatory DNA sequences in eukaryotes and are mostly located in the non-coding sequences of genes, namely, intergenic regions and introns. The essential characteristic of an enhancer is the ability to activate proximal genes, e.g., a reporter gene in a reporter assay, regardless of orientation, relative position, and distance from the gene. These characteristics are ascribed to the interaction (spatial proximity) of the enhancer sequence and the gene promoter via DNA looping, discussed in the latter part of this chapter.Developmentally regulated genes are associated with multiple enhancers carrying distinct cell and developmental stage specificities, which form arrays on the genome. We discuss the array of enhancers regulating the Sox2 gene as a paradigm. Sox2 enhancers are the best studied enhancers of a single gene in developmental regulation. In addition, the Sox2 gene is located in a genomic region with a very sparse gene distribution (no other protein-coding genes in ~1.6Mb in the mouse genome), termed a "gene desert," which means that most identified enhancers in the region are associated with Sox2 regulation. Furthermore, the importance of the Sox2 gene in stem cell regulation and neural development justifies focusing on Sox2-associated enhancers.

4mCPred-GSIMP: Predicting DNA N4-methylcytosine sites in the mouse genome with multi-Scale adaptive features extraction and fusion.

The epigenetic modification of DNA N4-methylcytosine (4mC) is vital for controlling DNA replication and expression. It is crucial to pinpoint 4mC's location to comprehend its role in physiological and pathological processes. However, accurate 4mC detection is difficult to achieve due to technical constraints. In this paper, we propose a deep learning-based approach 4mCPred-GSIMP for predicting 4mC sites in the mouse genome. The approach encodes DNA sequences using four feature encoding methods and combines multi-scale convolution and improved selective kernel convolution to adaptively extract and fuse features from different scales, thereby improving feature representation and optimization effect. In addition, we also use convolutional residual connections, global response normalization and pointwise convolution techniques to optimize the model. On the independent test dataset, 4mCPred-GSIMP shows high sensitivity, specificity, accuracy, Matthews correlation coefficient and area under the curve, which are 0.7812, 0.9312, 0.8562, 0.7207 and 0.9233, respectively. Various experiments demonstrate that 4mCPred-GSIMP outperforms existing prediction tools.

Optimal Structured Matrix Approximation for Robustness to Incomplete Biosequence Data.

We propose a general method for optimally approximating an arbitrary matrix M by a structured matrix T (circulant, Toeplitz/Hankel, etc.) and examine its use for estimating the spectra of genomic linkage disequilibrium matrices. This application is prototypical of a variety of genomic and proteomic problems that demand robustness to incomplete biosequence information. We perform a simulation study and corroborative test of our method using real genomic data from the Mouse Genome Database [1]. The results confirm the predicted utility of the method and provide strong evidence of its potential value to a wide range of bioinformatics applications. Our optimal general matrix approximation method is expected to be of independent interest to an even broader range of applications in applied mathematics and engineering.

Identification of rat Vstm1 with conservative anti-inflammatory effect between rat and human homologs

Human VSTM1 (also known as SIRL1) is an inhibitory immune checkpoint receptor involved in leukocyte activation. Identification of the homologous genes in other species, such as mice and rats, will undoubtedly contribute to functional studies and clinical applications. Here, we successfully cloned the Vstm1 gene in rats, as supported by high-throughput sequencing data. However, Vstm1 is degenerated to a pseudogene in the mouse genome. Rat Vstm1 mRNA contains a complete open reading frame (ORF) of 630 nucleotides encoding 209 amino acids. Rat Vstm1 is highly expressed in bone marrow, especially in granulocytes. The expression levels of Vstm1 gradually increase with the development of granulocytes in bone marrow but are downregulated in response to inflammatory stimuli. Rat VSTM1 does not have an immunoreceptor tyrosine-based inhibitory motif (ITIM), however, it shows a conservative function of inflammatory inhibition with human VSTM1, and both are anti-correlated with many inflammatory cytokines, such as IL-1α and TNF-α. In bone marrow-derived macrophages (BMDMs), either rat or human VSTM1 suppressed the secretion of inflammatory cytokines in response to LPS stimulation. Further analysis in lung cancer microenvironment revealed that VSTM1 is mainly expressed in myeloid cells, anti-correlated with inflammatory cytokines and associated with tumor development and metastasis.

Combining Off-flow, a Nextflow-coded program, and whole genome sequencing reveals unintended genetic variation in CRISPR/Cas-edited iPSCs

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas nucleases and human induced pluripotent stem cell (iPSC) technology can reveal deep insight into the genetic and molecular bases of human biology and disease. Undesired editing outcomes, both on-target (at the edited locus) and off-target (at other genomic loci) hinder the application of CRISPR-Cas nucleases. We developed Off-flow, a Nextflow-coded bioinformatic workflow that takes a specific guide sequence and Cas protein input to call four separate off-target prediction programs (CHOPCHOP, Cas-Offinder, CRISPRitz, CRISPR-Offinder) to output a comprehensive list of predicted off-target sites. We applied it to whole genome sequencing (WGS) data to investigate the occurrence of unintended effects in human iPSCs that underwent repair or insertion of disease-related variants by homology-directed repair. Off-flow identified a 3-base-pair-substitution and a mono-allelic genomic deletion at the target loci, KCNQ2, in 2 clones. Unbiased WGS analysis further identified off-target missense variants and a mono-allelic genomic deletion at the targeted locus, GNAQ, in 10 clones. On-target substitution and deletions had escaped standard PCR and Sanger sequencing analysis, while missense variants at other genomic loci were not detected by Off-flow. We used these results to filter out iPSC clones for subsequent functional experiments. Off-flow, which we make publicly available, works for human and mouse genomes currently and can be adapted for other genomes. Off-flow and WGS analysis can improve the integrity of studies using CRISPR/Cas-edited cells and animal models.

ZSCAN4-binding motif-TGCACAC is conserved and enriched in CA/TG microsatellites in both mouse and human genomes.

The Zinc finger and SCAN domain containing 4 (ZSCAN4) protein, expressed transiently in pluripotent stem cells, gametes, and early embryos, extends telomeres, enhances genome stability, and improves karyotypes in mouse embryonic stem (mES) cells. To gain insights into the mechanism of ZSCAN4 function, we identified genome-wide binding sites of endogenous ZSCAN4 protein using ChIP-seq technology in mouse and human ES cells, where the expression of endogenous ZSCAN4 was induced by treating cells with retinoic acids or by overexpressing DUX4. We revealed that both mouse and human ZSCAN4 bind to the TGCACAC motif located in CA/TG microsatellite repeats, which are known to form unstable left-handed duplexes called Z-DNA that can induce double-strand DNA breaks and mutations. These ZSCAN4 binding sites are mostly located in intergenic and intronic regions of the genomes. By generating ZSCAN4 knockout in human ES cells, we showed that ZSCAN4 does not seem to be involved in transcriptional regulation. We also found that ectopic expression of mouse ZSCAN4 enhances the suppression of chromatin at ZSCAN4-binding sites. These results together suggest that some of the ZSCAN4 functions are mediated by binding to the error-prone regions in mouse and human genomes.