The mouse Gene Expression Database (GXD): 2021 update

Fine mapping and candidate gene analysis of Dravet syndrome modifier loci on mouse chromosomes 7 and 8

Dravet syndrome is a developmental and epileptic encephalopathy (DEE) characterized by intractable seizures, comorbidities related to developmental, cognitive, and motor delays, and a high mortality burden due to sudden unexpected death in epilepsy (SUDEP). Most Dravet syndrome cases are attributed to SCN1A haploinsufficiency, with genetic modifiers and environmental factors influencing disease severity. Mouse models with heterozygous deletion of Scn1a recapitulate key features of Dravet syndrome, including seizures and premature mortality; however, severity varies depending on genetic background. Here, we refined two Dravet survival modifier (Dsm) loci, Dsm2 on chromosome 7 and Dsm3 on chromosome 8, using interval-specific congenic (ISC) mapping. Dsm2 was complex and encompassed at least two separate loci, while Dsm3 was refined to a single locus. Candidate modifier genes within these refined loci were prioritized based on brain expression, strain-dependent differences, and biological relevance to seizures or epilepsy. High priority candidate genes for Dsm2 include Nav2, Ptpn5, Ldha, Dbx1, Prmt3 and Slc6a5, while Dsm3 has a single high priority candidate, Psd3. This study underscores the complex genetic architecture underlying Dravet syndrome and provides insights into potential modifier genes that could influence disease severity and serve as novel therapeutic targets.

Decoding Molecular Bases of Rodent Social Hetero-Grooming Behavior Using in Silico Analyses and Bioinformatics Tools

Highly prevalent in laboratory rodents, 'social' hetero-grooming behavior is translationally relevant to modeling a wide range of neuropsychiatric disorders. Here, we comprehensively evaluated all known to date mouse genes linked to aberrant hetero-grooming phenotype, and applied bioinformatics tools to construct a network of their established protein-protein interactions (PPI). We next identified several distinct molecular clusters within this complex network, including neuronal differentiation, cytoskeletal, WNT-signaling and synapsins-associated pathways. Using additional bioinformatics analyses, we further identified 'central' (hub) proteins within these molecular clusters, likely key for mouse hetero-grooming behavior. Overall, a more comprehensive characterization of intricate molecular pathways linked to aberrant rodent grooming may markedly advance our understanding of underlying cellular mechanisms and related neurological disorders, eventually helping discover novel targets for their pharmacological or gene therapy interventions.

Therapeutic targets for muscle weakness in older adults: proteome-wide Mendelian randomization and colocalization analyses

BACKGROUND

Recent research highlights the importance of muscular strength as a key factor in physical fitness, a strong indicator of overall mortality risk, and a vital target for preventing chronic diseases. This study used a proteome-wide Mendelian randomization analysis plus colocalization analysis for low hand grip strength to explore potential therapeutic targets for muscle weakness.

METHODS

We conducted two two-sample Mendelian randomization analyses from four cohorts to identify and validate the causal relationship between plasma proteins and low grip strength. We also employed bidirectional Mendelian randomization analysis with Steiger filtering, Bayesian co-localization, and phenotype scanning to detect reverse causality, thereby consolidating our Mendelian randomization findings. Downstream analyses were also undertaken of identified proteins, including knockout models, enrichment analyses, and protein-protein interaction networks. Finally, we assessed the druggability of the identified proteins.

RESULTS

At Bonferroni significance (P < 6.82 × 10-5), Mendelian randomization analysis revealed that three proteins were causally associated with low grip strength. Increased MGP (OR = 0.85) and HP (OR = 0.96) decreased the risk of low grip strength, whereas elevated ART4 (OR = 1.06) increased the risk of low grip strength. None of the three proteins had reverse causality with low grip strength. Bayesian co-localization suggested that MGP shared the same variant with low grip strength (coloc.abf-PPH4 = 0.826). Further downstream analyses showed that MGP, which is highly expressed in musculoskeletal system, is a potential novel target for muscle weakness.

CONCLUSIONS

The proteome-wide Mendelian randomization investigation identified three proteins associated with the risk of muscle weakness. MGP, HP, and ART4 deserve further investigation as potential therapeutic targets for muscle weakness.

SPAG7 deletion causes intrauterine growth restriction, resulting in adulthood obesity and metabolic dysfunction

From a forward mutagenetic screen to discover mutations associated with obesity, we identified mutations in the Spag7 gene linked to metabolic dysfunction in mice. Here, we show that SPAG7 KO mice are born smaller and develop obesity and glucose intolerance in adulthood. This obesity does not stem from hyperphagia, but a decrease in energy expenditure. The KO animals also display reduced exercise tolerance and muscle function due to impaired mitochondrial function. Furthermore, SPAG7-deficiency in developing embryos leads to intrauterine growth restriction, brought on by placental insufficiency, likely due to abnormal development of the placental junctional zone. This insufficiency leads to loss of SPAG7-deficient fetuses in utero and reduced birth weights of those that survive. We hypothesize that a 'thrifty phenotype' is ingrained in SPAG7 KO animals during development that leads to adult obesity. Collectively, these results indicate that SPAG7 is essential for embryonic development and energy homeostasis later in life.

The changing metabolic landscape of bile acids - keys to metabolism and immune regulation

Bile acids regulate nutrient absorption and mitochondrial function, they establish and maintain gut microbial community composition and mediate inflammation, and they serve as signalling molecules that regulate appetite and energy homeostasis. The observation that there are hundreds of bile acids, especially many amidated bile acids, necessitates a revision of many of the classical descriptions of bile acids and bile acid enzyme functions. For example, bile salt hydrolases also have transferase activity. There are now hundreds of known modifications to bile acids and thousands of bile acid-associated genes, especially when including the microbiome, distributed throughout the human body (for example, there are >2,400 bile salt hydrolases alone). The fact that so much of our genetic and small-molecule repertoire, in both amount and diversity, is dedicated to bile acid function highlights the centrality of bile acids as key regulators of metabolism and immune homeostasis, which is, in large part, communicated via the gut microbiome.

Insights into digit evolution from a fate map study of the forearm using Chameleon, a new transgenic chicken line

The cellular and genetic networks that contribute to the development of the zeugopod (radius and ulna of the forearm, tibia and fibula of the leg) are not well understood, although these bones are susceptible to loss in congenital human syndromes and to the action of teratogens such as thalidomide. Using a new fate-mapping approach with the Chameleon transgenic chicken line, we show that there is a small contribution of SHH-expressing cells to the posterior ulna, posterior carpals and digit 3. We establish that although the majority of the ulna develops in response to paracrine SHH signalling in both the chicken and mouse, there are differences in the contribution of SHH-expressing cells between mouse and chicken as well as between the chicken ulna and fibula. This is evidence that, although zeugopod bones are clearly homologous according to the fossil record, the gene regulatory networks that contribute to their development and evolution are not fixed.

Systems Genetics Analyses Reveals Key Genes Related to Behavioral Traits in the Striatum of CFW Mice

The striatum plays a central role in directing many complex behaviors ranging from motor control to action choice and reward learning. In our study, we used 55 male CFW mice with rapid decay linkage disequilibrium to systematically mine the striatum-related behavioral functional genes by analyzing their striatal transcriptomes and 79 measured behavioral phenotypic data. By constructing a gene coexpression network, we clustered the genes into 13 modules, with most of them being positively correlated with motor traits. Based on functional annotations as well as Fisher's exact and hypergeometric distribution tests, brown and magenta modules were identified as core modules. They were significantly enriched for striatal-related functional genes. Subsequent Mendelian randomization analysis verified the causal relationship between the core modules and dyskinesia. Through the intramodular gene connectivity analysis, Adcy5 and Kcnma1 were identified as brown and magenta module hub genes, respectively. Knock outs of both Adcy5 and Kcnma1 lead to motor dysfunction in mice, and KCNMA1 acts as a risk gene for schizophrenia and smoking addiction in humans. We also evaluated the cellular composition of each module and identified oligodendrocytes in the striatum to have a positive role in motor regulation.

Mouse Genome Informatics: an integrated knowledgebase system for the laboratory mouse

Mouse Genome Informatics (MGI) is a federation of expertly curated information resources designed to support experimental and computational investigations into genetic and genomic aspects of human biology and disease using the laboratory mouse as a model system. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are core MGI databases that share data and system architecture. MGI serves as the central community resource of integrated information about mouse genome features, variation, expression, gene function, phenotype, and human disease models acquired from peer-reviewed publications, author submissions, and major bioinformatics resources. To facilitate integration and standardization of data, biocuration scientists annotate using terms from controlled metadata vocabularies and biological ontologies (e.g. Mammalian Phenotype Ontology, Mouse Developmental Anatomy, Disease Ontology, Gene Ontology, etc.), and by applying international community standards for gene, allele, and mouse strain nomenclature. MGI serves basic scientists, translational researchers, and data scientists by providing access to FAIR-compliant data in both human-readable and compute-ready formats. The MGI resource is accessible at https://informatics.jax.org. Here, we present an overview of the core data types represented in MGI and highlight recent enhancements to the resource with a focus on new data and functionality for MGD and GXD.

Cross-Domain Text Mining of Pathophysiological Processes Associated with Diabetic Kidney Disease

Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease worldwide. This study's goal was to identify the signaling drivers and pathways that modulate glomerular endothelial dysfunction in DKD via artificial intelligence-enabled literature-based discovery. Cross-domain text mining of 33+ million PubMed articles was performed with SemNet 2.0 to identify and rank multi-scalar and multi-factorial pathophysiological concepts related to DKD. A set of identified relevant genes and proteins that regulate different pathological events associated with DKD were analyzed and ranked using normalized mean HeteSim scores. High-ranking genes and proteins intersected three domains-DKD, the immune response, and glomerular endothelial cells. The top 10% of ranked concepts were mapped to the following biological functions: angiogenesis, apoptotic processes, cell adhesion, chemotaxis, growth factor signaling, vascular permeability, the nitric oxide response, oxidative stress, the cytokine response, macrophage signaling, NFκB factor activity, the TLR pathway, glucose metabolism, the inflammatory response, the ERK/MAPK signaling response, the JAK/STAT pathway, the T-cell-mediated response, the WNT/β-catenin pathway, the renin-angiotensin system, and NADPH oxidase activity. High-ranking genes and proteins were used to generate a protein-protein interaction network. The study results prioritized interactions or molecules involved in dysregulated signaling in DKD, which can be further assessed through biochemical network models or experiments.

Fine Mapping and Candidate Gene Analysis of Dravet Syndrome Modifier Loci on Mouse Chromosomes 7 and 8

Dravet syndrome is a developmental and epileptic encephalopathy (DEE) characterized by intractable seizures, comorbidities related to developmental, cognitive, and motor delays, and a high mortality burden due to sudden unexpected death in epilepsy (SUDEP). Most Dravet syndrome cases are attributed to SCN1A haploinsufficiency, with genetic modifiers and environmental factors influencing disease severity. Mouse models with heterozygous deletion of Scn1a recapitulate key features of Dravet syndrome, including seizures and premature mortality; however, severity varies depending on genetic background. Here, we refined two Dravet survival modifier (Dsm) loci, Dsm2 on chromosome 7 and Dsm3 on chromosome 8, using interval-specific congenic (ISC) mapping. Dsm2 was complex and encompassed at least two separate loci, while Dsm3 was refined to a single locus. Candidate modifier genes within these refined loci were prioritized based on brain expression, strain-dependent differences, and biological relevance to seizures or epilepsy. High priority candidate genes for Dsm2 include Nav2, Ptpn5, Ldha, Dbx1, Prmt3 and Slc6a5, while Dsm3 has a single high priority candidate, Psd3. This study underscores the complex genetic architecture underlying Dravet syndrome and provides insights into potential modifier genes that could influence disease severity and serve as novel therapeutic targets.

The underappreciated diversity of bile acid modifications

The repertoire of modifications to bile acids and related steroidal lipids by host and microbial metabolism remains incompletely characterized. To address this knowledge gap, we created a reusable resource of tandem mass spectrometry (MS/MS) spectra by filtering 1.2 billion publicly available MS/MS spectra for bile-acid-selective ion patterns. Thousands of modifications are distributed throughout animal and human bodies as well as microbial cultures. We employed this MS/MS library to identify polyamine bile amidates, prevalent in carnivores. They are present in humans, and their levels alter with a diet change from a Mediterranean to a typical American diet. This work highlights the existence of many more bile acid modifications than previously recognized and the value of leveraging public large-scale untargeted metabolomics data to discover metabolites. The availability of a modification-centric bile acid MS/MS library will inform future studies investigating bile acid roles in health and disease.

GFRAL Is Widely Distributed in the Brain and Peripheral Tissues of Mice

In 2017, four independent publications described the glial cell-derived neurotrophic factor (GDNF) receptor alpha-like (GFRAL) as receptor for the growth differentiation factor 15 (GDF15, also MIC-1, NAG-1) with an expression exclusively in the mice brainstem area postrema (AP) and nucleus tractus solitarii (NTS) where it mediates effects of GDF15 on reduction of food intake and body weight. GDF15 is a cell stress cytokine with a widespread expression and pleiotropic effects, which both seem to be in contrast to the reported highly specialized localization of its receptor. This discrepancy prompts us to re-evaluate the expression pattern of GFRAL in the brain and peripheral tissues of mice. In this detailed immunohistochemical study, we provide evidence for a more widespread distribution of this receptor. Apart from the AP/NTS region, GFRAL-immunoreactivity was found in the prefrontal cortex, hippocampus, nucleus arcuatus and peripheral tissues including liver, small intestine, fat, kidney and muscle tissues. This widespread receptor expression, not taken into consideration so far, may explain the multiple effects of GDF-15 that are not yet assigned to GFRAL. Furthermore, our results could be relevant for the development of novel pharmacological therapies for physical and mental disorders related to body image and food intake, such as eating disorders, cachexia and obesity.

Scrutinizing the human TEX genes in the context of human male infertility

BACKGROUND

Infertility affects around 15% of all couples worldwide and is increasingly linked to variants in genes specifically expressed in the testis. Well-established causes of male infertility include pathogenic variants in the genes TEX11, TEX14, and TEX15, while few studies have recently reported variants in TEX13B, TEX13C, FAM9A (TEX39A), and FAM9B (TEX39B).

OBJECTIVES

We aimed at screening for novel potential candidate genes among the human TEX ("testis expressed") genes as well as verifying previously described disease associations in this set of genes.

MATERIALS AND METHODS

To this end, we screened the exome sequencing data of 1305 men, including 1056 crypto- and azoospermic individuals, and determined cell-specific expression by analyzing testis-specific single-cell RNA sequencing data for genes with identified variants. To investigate the overarching role in male fertility, we generated testis-specific knockdown (KD) models of all 10 orthologous TEX genes in Drosophila melanogaster.

RESULTS

We detected rare potential disease-causing variants in TEX10, TEX13A, TEX13B, TEX13C, TEX13D, ZFAND3 (TEX27), TEX33, FAM9A (TEX39A), and FAM9B (TEX39B), in 28 infertile men, of which 15 men carried variants in TEX10, TEX27, and TEX33. The KD of TEX2, TEX9, TEX10, TEX13, ZFAND3 (TEX27), TEX28, TEX30, NFX1 (TEX42), TEX261, and UTP4 (TEX292) in Drosophila resulted in normal fertility.

DISCUSSION

Based on our findings, the autosomal dominant predicted genes TEX10 and ZFAND3 (TEX27) and the autosomal recessive predicted gene TEX33, which all three are conceivably required for germ cell maturation, were identified as novel potential candidate genes for human non-obstructive azoospermia. We additionally identified hemizygous loss-of-function (LoF) variants in TEX13B, TEX13C, and FAM9A (TEX39A) as unlikely monogenic culprits of male infertility as LoF variants were also found in control men.

CONCLUSION

Our findings concerning the X-linked genes TEX13B, TEX13C, and FAM9A (TEX39A) contradict previous reports and will decrease false-positive reports in genetic diagnostics of azoospermic men.

Cis-regulatory control of mammalian Trps1 gene expression

TRPS1 serves as the causative gene for tricho-rhino phalangeal syndrome, known for its craniofacial and skeletal abnormalities. The Trps1 gene encodes a protein that represses Wnt signaling through strong interactions with Wnt signaling inhibitors. The identification of genomic cis-acting regulatory sequences governing Trps1 expression is crucial for understanding its role in embryogenesis. Nevertheless, to date, no investigations have been conducted concerning these aspects of Trps1. To identify deeply conserved noncoding elements (CNEs) within the Trps1 locus, we employed a comparative genomics approach, utilizing slowly evolving fish such as coelacanth and spotted gar. These analyses resulted in the identification of eight CNEs in the intronic region of the Trps1 gene. Functional characterization of these CNEs in zebrafish revealed their regulatory potential in various tissues, including pectoral fins, heart, and pharyngeal arches. RNA in-situ hybridization experiments revealed concordance between the reporter expression pattern induced by the identified set of CNEs and the spatial expression pattern of the trps1 gene in zebrafish. Comparative in vivo data from zebrafish and mice for CNE7/hs919 revealed conserved functions of these enhancers. Each of these eight CNEs was further investigated in cell line-based reporter assays, revealing their repressive potential. Taken together, in vivo and in vitro assays suggest a context-dependent dual functionality for the identified set of Trps1-associated CNE enhancers. This functionally characterized set of CNE-enhancers will contribute to a more comprehensive understanding of the developmental roles of Trps1 and can aid in the identification of noncoding DNA variants associated with human diseases.

Identification and characterization of Hand2 upstream genomic enhancers active in developing stomach and limbs

BACKGROUND

The bHLH transcription factor HAND2 plays important roles in the development of the embryonic heart, face, limbs, and sympathetic and enteric nervous systems. To define how and when HAND2 regulates these developmental systems, requires understanding the transcriptional regulation of Hand2.

RESULTS

Remarkably, Hand2 is flanked by an extensive upstream gene desert containing a potentially diverse enhancer landscape. Here, we screened the regulatory interval 200 kb proximal to Hand2 for putative enhancers using evolutionary conservation and histone marks in Hand2-expressing tissues. H3K27ac signatures across embryonic tissues pointed to only two putative enhancer regions showing deep sequence conservation. Assessment of the transcriptional enhancer potential of these elements using transgenic reporter lines uncovered distinct in vivo enhancer activities in embryonic stomach and limb mesenchyme, respectively. Activity of the identified stomach enhancer was restricted to the developing antrum and showed expression within the smooth muscle and enteric neurons. Surprisingly, the activity pattern of the limb enhancer did not overlap Hand2 mRNA but consistently yielded a defined subectodermal anterior expression pattern within multiple transgenic lines.

CONCLUSIONS

Together, these results start to uncover the diverse regulatory potential inherent to the Hand2 upstream regulatory interval.

The unusual gene architecture of polyubiquitin is created by dual-specific splice sites

BACKGROUND

The removal of introns occurs through the splicing of a 5' splice site (5'ss) with a 3' splice site (3'ss). These two elements are recognized by distinct components of the spliceosome. However, introns in higher eukaryotes contain many matches to the 5' and 3' splice-site motifs that are presumed not to be used.

RESULTS

Here, we find that many of these sites can be used. We also find occurrences of the AGGT motif that can function as either a 5'ss or a 3'ss-previously referred to as dual-specific splice sites (DSSs)-within introns. Analysis of the Sequence Read Archive reveals a 3.1-fold enrichment of DSSs relative to expectation, implying synergy between the ability to function as a 5'ss and 3'ss. Despite this suggested mechanistic advantage, DSSs are 2.7- and 4.7-fold underrepresented in annotated 5' and 3' splice sites. A curious exception is the polyubiquitin gene UBC, which contains a tandem array of DSSs that precisely delimit the boundary of each ubiquitin monomer. The resulting isoforms splice stochastically to include a variable number of ubiquitin monomers. We found no evidence of tissue-specific or feedback regulation but note the 8.4-fold enrichment of DSS-spliced introns in tandem repeat genes suggests a driving role in the evolution of genes like UBC.

CONCLUSIONS

We find an excess of unannotated splice sites and the utilization of DSSs in tandem repeats supports the role of splicing in gene evolution. These findings enhance our understanding of the diverse and complex nature of the splicing process.

Identification of ancestral gnathostome Gli3 enhancers with activity in mammals

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

Transcriptome analyses of murine right and left maxilla-mandibular complex

OBJECTIVE

The objective of the study was to investigate differential gene expression between murine right and left maxilla-mandibular (MxMn) complexes.

SETTING AND SAMPLE POPULATION

Wild-type (WT) C57BL/6 embryonic (E) day 14.5 (n = 3) and 18.5 (n = 3) murine embryos.

METHODS

The E14.5 and 18.5 embryos were harvested and hemi-sectioned the MxMn complexes into right and left halves in the mid-sagittal plane. We isolated total RNA using Trizol reagent and further purified using the RNA-easy kit (QIAGEN). We confirmed equal expression of house-keeping genes in right and left halves using RT-PCR and then performed paired-end whole mRNA sequencing in LC Sciences (Houston, TX) followed by differential transcript analyses (>1 or <-1 log fold change; p < .05; q < .05; and FPKM >0.5 in 2/3 samples). The Mouse Genome Informatics and Online Mendelian Inheritance in Man databases as well as gnomAD constraint scores were used to prioritize differentially expressed transcripts.

RESULTS

There were 19 upregulated and 19 downregulated transcripts at E14.5 and 8 upregulated and 17 downregulated transcripts at E18.5 time-points. These differentially expressed transcripts were statistically significant and shown to be associated with craniofacial phenotypes in mouse models. These transcripts also have significant gnomAD constraint scores and are enriched in biological processes critical for embryogenesis.

CONCLUSIONS

We identified significant differential expression of transcripts between E14.5 and 18.5 murine right and left MxMn complexes. These findings when extrapolated to humans, they may provide a biological basis for facial asymmetry. Further experiments are required to validate these findings in murine models with craniofacial asymmetry.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Thalamic Neuron Resilience during Osmotic Demyelination Syndrome (ODS) Is Revealed by Primary Cilium Outgrowth and ADP-ribosylation factor-like protein 13B Labeling in Axon Initial Segment

A murine osmotic demyelinating syndrome (ODS) model was developed through chronic hyponatremia, induced by desmopressin subcutaneous implants, followed by precipitous sodium restoration. The thalamic ventral posterolateral (VPL) and ventral posteromedial (VPM) relay nuclei were the most demyelinated regions where neuroglial damage could be evidenced without immune response. This report showed that following chronic hyponatremia, 12 h and 48 h time lapses after rebalancing osmolarity, amid the ODS-degraded outskirts, some resilient neuronal cell bodies built up primary cilium and axon hillock regions that extended into axon initial segments (AIS) where ADP-ribosylation factor-like protein 13B (ARL13B)-immunolabeled rod-like shape content was revealed. These AIS-labeled shaft lengths appeared proportional with the distance of neuronal cell bodies away from the ODS damaged epicenter and time lapses after correction of hyponatremia. Fine structure examination verified these neuron abundant transcriptions and translation regions marked by the ARL13B labeling associated with cell neurotubules and their complex cytoskeletal macromolecular architecture. This necessitated energetic transport to organize and restore those AIS away from the damaged ODS core demyelinated zone in the murine model. These labeled structures could substantiate how thalamic neuron resilience occurred as possible steps of a healing course out of ODS.

Late fetal hematopoietic failure results from ZBTB11 deficiency despite abundant HSC specification

Hematopoiesis produces diverse blood cell lineages to meet the basal needs and sudden demands of injury or infection. A rapid response to such challenges requires the expansion of specific lineages and a prompt return to balanced steady-state levels, necessitating tightly coordinated regulation. Previously we identified a requirement for the zinc finger and broad complex, tramtrak, bric-a-brac domain-containing 11 (ZBTB11) transcription factor in definitive hematopoiesis using a forward genetic screen for zebrafish myeloid mutants. To understand its relevance to mammalian systems, we extended these studies to mice. When Zbtb11 was deleted in the hematopoietic compartment, embryos died at embryonic day (E) 18.5 with hematopoietic failure. Zbtb11 hematopoietic knockout (Zbtb11hKO) hematopoietic stem cells (HSCs) were overabundantly specified from E14.5 to E17.5 compared with those in controls. Overspecification was accompanied by loss of stemness, inability to differentiate into committed progenitors and mature lineages in the fetal liver, failure to seed fetal bone marrow, and total hematopoietic failure. The Zbtb11hKO HSCs did not proliferate in vitro and were constrained in cell cycle progression, demonstrating the cell-intrinsic role of Zbtb11 in proliferation and cell cycle regulation in mammalian HSCs. Single-cell RNA sequencing analysis identified that Zbtb11-deficient HSCs were underrepresented in an erythroid-primed subpopulation and showed downregulation of oxidative phosphorylation pathways and dysregulation of genes associated with the hematopoietic niche. We identified a cell-intrinsic requirement for Zbtb11-mediated gene regulatory networks in sustaining a pool of maturation-capable HSCs and progenitor cells.

The RELT Family of Proteins: An Increasing Awareness of Their Importance for Cancer, the Immune System, and Development

This review highlights Receptor Expressed in Lymphoid Tissues (RELT), a Tumor Necrosis Factor Superfamily member, and its two paralogs, RELL1 and RELL2. Collectively, these three proteins are referred to as RELTfms and have gained much interest in recent years due to their association with cancer and other human diseases. A thorough knowledge of their physiological functions, including the ligand for RELT, is lacking, yet emerging evidence implicates RELTfms in a variety of processes including cytokine signaling and pathways that either promote cell death or survival. T cells from mice lacking RELT exhibit increased responses against tumors and increased inflammatory cytokine production, and multiple lines of evidence indicate that RELT may promote an immunosuppressive environment for tumors. The relationship of individual RELTfms in different cancers is not universal however, as evidence indicates that individual RELTfms may be risk factors in certain cancers yet appear to be protective in other cancers. RELTfms are important for a variety of additional processes related to human health including microbial pathogenesis, inflammation, behavior, reproduction, and development. All three proteins have been strongly conserved in all vertebrates, and this review aims to provide a clearer understanding of the current knowledge regarding these interesting proteins.

Acyl coenzyme A binding protein (ACBP): An aging- and disease-relevant "autophagy checkpoint"

Acyl coenzyme A binding protein (ACBP), also known as diazepam-binding inhibitor (DBI), is a phylogenetically ancient protein present in some eubacteria and the entire eukaryotic radiation. In several eukaryotic phyla, ACBP/DBI transcends its intracellular function in fatty acid metabolism because it can be released into the extracellular space. This ACBP/DBI secretion usually occurs in response to nutrient scarcity through an autophagy-dependent pathway. ACBP/DBI and its peptide fragments then act on a range of distinct receptors that diverge among phyla, namely metabotropic G protein-coupled receptor in yeast (and likely in the mammalian central nervous system), a histidine receptor kinase in slime molds, and ionotropic gamma-aminobutyric acid (GABA)A receptors in mammals. Genetic or antibody-mediated inhibition of ACBP/DBI orthologs interferes with nutrient stress-induced adaptations such as sporulation or increased food intake in multiple species, as it enhances lifespan or healthspan in yeast, plant leaves, nematodes, and multiple mouse models. These lifespan and healthspan-extending effects of ACBP/DBI suppression are coupled to the induction of autophagy. Altogether, it appears that neutralization of extracellular ACBP/DBI results in "autophagy checkpoint inhibition" to unleash the anti-aging potential of autophagy. Of note, in humans, ACBP/DBI levels increase in various tissues, as well as in the plasma, in the context of aging, obesity, uncontrolled infection or cardiovascular, inflammatory, neurodegenerative, and malignant diseases.

A non-coding insertional mutation of Grhl2 causes gene over-expression and multiple structural anomalies including cleft palate, spina bifida and encephalocele

Orofacial clefts, including cleft lip and palate (CL/P) and neural tube defects (NTDs) are among the most common congenital anomalies, but knowledge of the genetic basis of these conditions remains incomplete. The extent to which genetic risk factors are shared between CL/P, NTDs and related anomalies is also unclear. While identification of causative genes has largely focused on coding and loss of function mutations, it is hypothesized that regulatory mutations account for a portion of the unidentified heritability. We found that excess expression of Grainyhead-like 2 (Grhl2) causes not only spinal NTDs in Axial defects (Axd) mice but also multiple additional defects affecting the cranial region. These include orofacial clefts comprising midline cleft lip and palate and abnormalities of the craniofacial bones and frontal and/or basal encephalocele, in which brain tissue herniates through the cranium or into the nasal cavity. To investigate the causative mutation in the Grhl2Axd strain, whole genome sequencing identified an approximately 4 kb LTR retrotransposon insertion that disrupts the non-coding regulatory region, lying approximately 300 base pairs upstream of the 5' UTR. This insertion also lies within a predicted long non-coding RNA, oriented on the reverse strand, which like Grhl2 is over-expressed in Axd (Grhl2Axd) homozygous mutant embryos. Initial analysis of the GRHL2 upstream region in individuals with NTDs or cleft palate revealed rare or novel variants in a small number of cases. We hypothesize that mutations affecting the regulation of GRHL2 may contribute to craniofacial anomalies and NTDs in humans.

Microtubule-Associated Serine/Threonine (MAST) Kinases in Development and Disease

Microtubule-Associated Serine/Threonine (MAST) kinases represent an evolutionary conserved branch of the AGC protein kinase superfamily in the kinome. Since the discovery of the founding member, MAST2, in 1993, three additional family members have been identified in mammals and found to be broadly expressed across various tissues, including the brain, heart, lung, liver, intestine and kidney. The study of MAST kinases is highly relevant for unraveling the molecular basis of a wide range of different human diseases, including breast and liver cancer, myeloma, inflammatory bowel disease, cystic fibrosis and various neuronal disorders. Despite several reports on potential substrates and binding partners of MAST kinases, the molecular mechanisms that would explain their involvement in human diseases remain rather obscure. This review will summarize data on the structure, biochemistry and cell and molecular biology of MAST kinases in the context of biomedical research as well as organismal model systems in order to provide a current profile of this field.

Plp1 in the enteric nervous system is preferentially expressed during early postnatal development in mouse as DM20, whose expression appears reliant on an intronic enhancer

Recently, the myelin proteolipid protein gene (Plp1) was shown to be expressed in the glia of the enteric nervous system (ENS) in mouse. However, beyond this, not much is known about its expression in the intestine. To address this matter, we investigated Plp1 expression at the mRNA and protein levels in the intestine of mice at different ages (postnatal days 2, 9, 21, and 88). In this study, we show that Plp1 expression preferentially occurs during early postnatal development, primarily as the DM20 isoform. Western blot analysis indicated that DM20 migrated according to its formula weight when isolated from the intestine. However, mobilities of both PLP and DM20 were faster than expected when procured from the brain. The 6.2hPLP(+)Z/FL transgene, which uses the first half of the human PLP1 gene to drive expression of a lacZ reporter gene, recapitulated the developmental pattern observed with the native gene in the intestine, indicating that it can be used as a proxy for Plp1 gene expression. As such, the relative levels of β-galactosidase (β-gal) activity emanating from the 6.2hPLP(+)Z/FL transgene suggest that Plp1 expression is highest in the duodenum, and decreases successively along the segments, toward the colon. Moreover, removal of the wmN1 enhancer region from the transgene (located within Plp1 intron 1) resulted in a dramatic reduction in both transgene mRNA levels and β-gal activity in the intestine, throughout development, suggesting that this region contains a regulatory element crucial for Plp1 expression. This is consistent with earlier studies in both the central and peripheral nervous systems, indicating that it may be a common (if not universal) means by which Plp1 gene expression is governed.

Sohlh1 loss of function male and female infertility model impacts overall health beyond gonadal dysfunction in mice†

Reproductive longevity is associated with health outcomes. Early menopause, loss of ovarian function, and male infertility are linked to shorter lifespan and increased adverse health outcomes. Here we examined the extragonadal effects of whole animal loss of spermatogenesis and oogenesis specific basic helix-loop-helix 1 (Sohlh1) gene in mice, a well-described mouse model of female and male infertility. Sohlh1 encodes a transcription factor that is primarily expressed in the male and female germline and regulates germline differentiation. The Sohlh1 knockout mouse model, just like human individuals with SOHLH1 loss of function, presents with hypergonadotropic hypogonadism and loss of ovarian function in females and impaired spermatogenesis in males, with a seemingly gonad restricted phenotype in both sexes. However, extragonadal phenotyping revealed that Sohlh1 deficiency leads to abnormal immune profiles in the blood and ovarian tissues of female animals, sex-specific alterations of metabolites, and behavior and cognition changes. Altogether, these results show that Sohlh1 deficiency impacts overall health in both male and female mice.

The IgG4 hinge with CD28 transmembrane domain improves VHH-based CAR T cells targeting a membrane-distal epitope of GPC1 in pancreatic cancer

Heterogeneous antigen expression is a key barrier influencing the activity of chimeric antigen receptor (CAR) T cells in solid tumors. Here, we develop CAR T cells targeting glypican-1 (GPC1), an oncofetal antigen expressed in pancreatic cancer. We report the generation of dromedary camel V_HH nanobody (D4)-based CAR T cells targeting GPC1 and the optimization of the hinge (H) and transmembrane domain (TM) to improve activity. We find that a structurally rigid IgG4H and CD28TM domain brings the two D4 fragments in proximity, driving CAR dimerization and leading to enhanced T-cell signaling and tumor regression in pancreatic cancer models with low antigen density in female mice. Furthermore, single-cell-based proteomic and transcriptomic analysis of D4-IgG4H-CD28TM CAR T cells reveals specific genes (e.g., HMGB1) associated with high T-cell polyfunctionality. This study demonstrates the potential of V_HH-based CAR T for pancreatic cancer therapy and provides an engineering strategy for developing potent CAR T cells targeting membrane-distal epitopes.

Determination of sex in cDNA samples from mice

The need to take sex into account in biomedical research is now recognized and mandated by funding institutions. In laboratory rodents, such as mice, sexing is usually performed anatomically or by genotyping using multiplex or simplex PCR techniques on genomic DNA. Here we present a simple RT-PCR-based method targeting Kdm5c and Kdm5d to determine genetic sex in mouse cDNA samples, allowing for retrospective sex determination.

Computational drug repurposing by exploiting large-scale gene expression data: Strategy, methods and applications

De novo drug development is an extremely complex, time-consuming and costly task. Urgent needs for therapies of various diseases have greatly accelerated searches for more effective drug development methods. Luckily, drug repurposing provides a new and effective perspective on disease treatment. Rapidly increased large-scale transcriptome data paints a detailed prospect of gene expression during disease onset and thus has received wide attention in the field of computational drug repurposing. However, how to efficiently mine transcriptome data and identify new indications for old drugs remains a critical challenge. This review discussed the irreplaceable role of transcriptome data in computational drug repurposing and summarized some representative databases, tools and strategies. More importantly, it proposed a practical guideline through establishing the correspondence between three gene expression data types and five strategies, which would facilitate researchers to adopt appropriate strategies to deeply mine large-scale transcriptome data and discover more effective therapies.

Angiopoietin-like 2 is essential to aortic valve development in mice

Aortic valve (AoV) abnormalities during embryogenesis are a major risk for the development of aortic valve stenosis (AVS) and cardiac events later in life. Here, we identify an unexpected role for Angiopoietin-like 2 (ANGPTL2), a pro-inflammatory protein secreted by senescent cells, in valvulogenesis. At late embryonic stage, mice knocked-down for Angptl2 (Angptl2-KD) exhibit a premature thickening of AoV leaflets associated with a dysregulation of the fine balance between cell apoptosis, senescence and proliferation during AoV remodeling and a decrease in the crucial Notch signalling. These structural and molecular abnormalities lead toward spontaneous AVS with elevated trans-aortic gradient in adult mice of both sexes. Consistently, ANGPTL2 expression is detected in human fetal semilunar valves and associated with pathways involved in cell cycle and senescence. Altogether, these findings suggest that Angptl2 is essential for valvulogenesis, and identify Angptl2-KD mice as an animal model to study spontaneous AVS, a disease with unmet medical need.

Delayed puberty, gonadotropin abnormalities and subfertility in male Padi2/Padi4 double knockout mice

BACKGROUND

Peptidylarginine deiminase enzymes (PADs) convert arginine residues to citrulline in a process called citrullination or deimination. Recently, two PADs, PAD2 and PAD4, have been linked to hormone signaling in vitro and the goal of this study was to test for links between PAD2/PAD4 and hormone signaling in vivo.

METHODS

Preliminary analysis of Padi2 and Padi4 single knockout (SKO) mice did not find any overt reproductive defects and we predicted that this was likely due to genetic compensation. To test this hypothesis, we created a Padi2/Padi4 double knockout (DKO) mouse model and tested these mice along with wild-type FVB/NJ (WT) and both strains of SKO mice for a range of reproductive defects.

RESULTS

Controlled breeding trials found that male DKO mice appeared to take longer to have their first litter than WT controls. This tendency was maintained when these mice were mated to either DKO or WT females. Additionally, unsexed 2-day old DKO pups and male DKO weanlings both weighed significantly less than their WT counterparts, took significantly longer than WT males to reach puberty, and had consistently lower serum testosterone levels. Furthermore, 90-day old adult DKO males had smaller testes than WT males with increased rates of germ cell apoptosis.

CONCLUSIONS

The Padi2/Padi4 DKO mouse model provides a new tool for investigating PAD function and outcomes from our studies provide the first in vivo evidence linking PADs with hormone signaling.

HIV-1 Tat Upregulates the Receptor for Advanced Glycation End Products and Superoxide Dismutase-2 in the Heart of Transgenic Mice

Cardiovascular disorder (CVD) is a common comorbidity in people living with HIV (PLWH). Although the underlying mechanisms are unknown, virotoxic HIV proteins, such as the trans-activator of transcription (Tat), likely contribute to CVD pathogenesis. Tat expression in mouse myocardium has been found to induce cardiac dysfunction and increase markers of endothelial toxicity. However, the role that Tat may play in the development of CVD pathogenesis is unclear. The capacity for Tat to impact cardiac function was assessed using AC16 human cardiomyocyte cells and adult male and female transgenic mice that conditionally expressed Tat [Tat(+)], or did not [Tat(-)]. In AC16 cardiomyocytes, Tat increased intracellular calcium. In Tat(+) mice, Tat expression was detected in both atrial and ventricular heart tissue. Tat(+) mice demonstrated an increased expression of the receptor for advanced glycation end products and superoxide dismutase-2 (SOD-2) in ventricular tissues compared to Tat(-) controls. No changes in SOD-1 or α-smooth muscle actin were observed. Despite Tat-mediated changes at the cellular level, no changes in echocardiographic measures were detected. Tat(+) mice had a greater proportion of ventricular mast cells and collagen; however, doxycycline exposure offset the latter effect. These data suggest that Tat exposure promotes cellular changes that can precede progression to CVD.

An anti-ACVR1 antibody exacerbates heterotopic ossification by fibro-adipogenic progenitors in fibrodysplasia ossificans progressiva mice

Fibrodysplasia ossificans progressiva (FOP) is a rare genetic disease characterized by progressive and catastrophic heterotopic ossification (HO) of skeletal muscle and associated soft tissues. FOP is caused by dominantly acting mutations in the gene encoding the bone morphogenetic protein (BMP) type I receptor, ACVR1 (ALK2), the most prevalent of which results in an arginine to histidine substitution at position 206[ACVR1(R206H)]. The fundamental pathological consequence of FOP-causing ACVR1 receptor mutations is to enable activin A to initiate canonical BMP signaling in fibro-adipogenic progenitors (FAPs), which drives HO. We developed a monoclonal blocking antibody (JAB0505) to the extracellular domain of ACVR1 and tested its effect on HO in two independent FOP mouse models. Although JAB0505 inhibited BMP-dependent gene expression in wild-type and ACVR1(R206H)-overexpressing cell lines, JAB0505 treatment profoundly exacerbated injury-induced HO. JAB0505-treated mice exhibited multiple, distinct foci of heterotopic lesions, suggesting an atypically broad anatomical domain of FAP recruitment to endochondral ossification. This was accompanied by dysregulated FAP population growth and an abnormally sustained immunological reaction following muscle injury. JAB0505 drove injury-induced HO in the absence of activin A, indicating that JAB0505 has receptor agonist activity. These data raise serious safety and efficacy concerns for the use of bivalent anti-ACVR1 antibodies to treat patients with FOP.

OmicsView: omics data analysis through interactive visual analytics

With advances in NGS technologies, transcriptional profiling of human tissue across many diseases is becoming more routine, leading to the generation of petabytes of data deposited in public repositories. There is a need for bench scientists with little computational expertise to be able to access and mine this data to understand disease pathology, identify robust biomarkers of disease and the effect of interventions (in vivo or in vitro). To this end we release an open source analytics and visualization platform for expression data called OmicsView, http://omicsview.org.

This platform comes preloaded with 1000 s of samples across many disease areas and normal tissue, including the GTEx database, all processed with a harmonized pipeline. We demonstrate the power and ease-of-use of the platform by means of a Crohn’s disease data mining exercise where we can quickly uncover disease pathology and identify strong biomarkers of disease and response to treatment.

MOET: a web-based gene set enrichment tool at the Rat Genome Database for multiontology and multispecies analyses

Biological interpretation of a large amount of gene or protein data is complex. Ontology analysis tools are imperative in finding functional similarities through overrepresentation or enrichment of terms associated with the input gene or protein lists. However, most tools are limited by their ability to do ontology-specific and species-limited analyses. Furthermore, some enrichment tools are not updated frequently with recent information from databases, thus giving users inaccurate, outdated or uninformative data. Here, we present MOET or the Multi-Ontology Enrichment Tool (v.1 released in April 2019 and v.2 released in May 2021), an ontology analysis tool leveraging data that the Rat Genome Database (RGD) integrated from in-house expert curation and external databases including the National Center for Biotechnology Information (NCBI), Mouse Genome Informatics (MGI), The Kyoto Encyclopedia of Genes and Genomes (KEGG), The Gene Ontology Resource, UniProt-GOA, and others. Given a gene or protein list, MOET analysis identifies significantly overrepresented ontology terms using a hypergeometric test and provides nominal and Bonferroni corrected P-values and odds ratios for the overrepresented terms. The results are shown as a downloadable list of terms with and without Bonferroni correction, and a graph of the P-values and number of annotated genes for each term in the list. MOET can be accessed freely from https://rgd.mcw.edu/rgdweb/enrichment/start.html.

A homozygous missense variant in laminin subunit beta 1 as candidate causal mutation of hemifacial microsomia in Romagnola cattle

Hemifacial microsomia (HFM) was diagnosed in a 9‐day‐old Romagnola calf. The condition was characterized by microtia of the left ear, anotia of the right ear, asymmetry of the face, and deafness. Magnetic resonance imaging revealed agenesis of the right pinna and both tympanic bullae, asymmetry of the temporal bones and temporomandibular joints, and right pontine meningocele. Brainstem auditory evoked responses confirmed the impaired auditory capacity. At gross post mortem examination, there was agenesis and hypoplasia of the right and the left external ear, respectively. No histological abnormalities were detected in the inner ears. A trio whole‐genome sequencing approach was carried out and identified a private homozygous missense variant in LAMB1 affecting a conserved residue (p.Arg668Cys). Genotyping of 221 Romagnola bulls revealed a carrier prevalence <2%. This represents a report of a LAMB1‐related autosomal recessive inherited disorder in domestic animals and adds LAMB1 to the candidate genes for HFM.

The Human Disease Ontology 2022 update

The Human Disease Ontology (DO) (www.disease-ontology.org) database, has significantly expanded the disease content and enhanced our userbase and website since the DO’s 2018 Nucleic Acids Research DATABASE issue paper. Conservatively, based on available resource statistics, terms from the DO have been annotated to over 1.5 million biomedical data elements and citations, a 10× increase in the past 5 years. The DO, funded as a NHGRI Genomic Resource, plays a key role in disease knowledge organization, representation, and standardization, serving as a reference framework for multiscale biomedical data integration and analysis across thousands of clinical, biomedical and computational research projects and genomic resources around the world. This update reports on the addition of 1,793 new disease terms, a 14% increase of textual definitions and the integration of 22 137 new SubClassOf axioms defining disease to disease connections representing the DO’s complex disease classification. The DO’s updated website provides multifaceted etiology searching, enhanced documentation and educational resources.

Mouse Genome Informatics (MGI): latest news from MGD and GXD

The Mouse Genome Informatics (MGI) database system combines multiple expertly curated community data resources into a shared knowledge management ecosystem united by common metadata annotation standards. MGI's mission is to facilitate the use of the mouse as an experimental model for understanding the genetic and genomic basis of human health and disease. MGI is the authoritative source for mouse gene, allele, and strain nomenclature and is the primary source of mouse phenotype annotations, functional annotations, developmental gene expression information, and annotations of mouse models with human diseases. MGI maintains mouse anatomy and phenotype ontologies and contributes to the development of the Gene Ontology and Disease Ontology and uses these ontologies as standard terminologies for annotation. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are MGI's two major knowledgebases. Here, we highlight some of the recent changes and enhancements to MGD and GXD that have been implemented in response to changing needs of the biomedical research community and to improve the efficiency of expert curation. MGI can be accessed freely at http://www.informatics.jax.org .

The Functions of Cytochrome P450 ω-hydroxylases and the Associated Eicosanoids in Inflammation-Related Diseases

The cytochrome P450 (CYP) ω-hydroxylases are a subfamily of CYP enzymes. While CYPs are the main metabolic enzymes that mediate the oxidation reactions of many endogenous and exogenous compounds in the human body, CYP ω-hydroxylases mediate the metabolism of multiple fatty acids and their metabolites via the addition of a hydroxyl group to the ω- or (ω-1)-C atom of the substrates. The substrates of CYP ω-hydroxylases include but not limited to arachidonic acid, docosahexaenoic acid, eicosapentaenoic acid, epoxyeicosatrienoic acids, leukotrienes, and prostaglandins. The CYP ω-hydroxylases-mediated metabolites, such as 20-hyroxyleicosatrienoic acid (20-HETE), 19-HETE, 20-hydroxyl leukotriene B4 (20-OH-LTB₄), and many ω-hydroxylated prostaglandins, have pleiotropic effects in inflammation and many inflammation-associated diseases. Here we reviewed the classification, tissue distribution of CYP ω-hydroxylases and the role of their hydroxylated metabolites in inflammation-associated diseases. We described up-regulation of CYP ω-hydroxylases may be a pathogenic mechanism of many inflammation-associated diseases and thus CYP ω-hydroxylases may be a therapeutic target for these diseases. CYP ω-hydroxylases-mediated eicosanods play important roles in inflammation as pro-inflammatory or anti-inflammatory mediators, participating in the process stimulated by cytokines and/or the process stimulating the production of multiple cytokines. However, most previous studies focused on 20-HETE,and further studies are needed for the function and mechanisms of other CYP ω-hydroxylases-mediated eicosanoids. We believe that our studies of CYP ω-hydroxylases and their associated eicosanoids will advance the translational and clinal use of CYP ω-hydroxylases inhibitors and activators in many diseases.

TISMO: syngeneic mouse tumor database to model tumor immunity and immunotherapy response

Syngeneic mouse models are tumors derived from murine cancer cells engrafted on genetically identical mouse strains. They are widely used tools for studying tumor immunity and immunotherapy response in the context of a fully functional murine immune system. Large volumes of syngeneic mouse tumor expression profiles under different immunotherapy treatments have been generated, although a lack of systematic collection and analysis makes data reuse challenging. We present Tumor Immune Syngeneic MOuse (TISMO), a database with an extensive collection of syngeneic mouse model profiles with interactive visualization features. TISMO contains 605 in vitro RNA-seq samples from 49 syngeneic cancer cell lines across 23 cancer types, of which 195 underwent cytokine treatment. TISMO also includes 1518 in vivo RNA-seq samples from 68 syngeneic mouse tumor models across 19 cancer types, of which 832 were from immune checkpoint blockade (ICB) studies. We manually annotated the sample metadata, such as cell line, mouse strain, transplantation site, treatment, and response status, and uniformly processed and quality-controlled the RNA-seq data. Besides data download, TISMO provides interactive web interfaces to investigate whether specific gene expression, pathway enrichment, or immune infiltration level is associated with differential immunotherapy response. TISMO is available at http://tismo.cistrome.org.

Antigen dominance hierarchies shape TCF1+ progenitor CD8 T cell phenotypes in tumors

CD8 T cell responses against different tumor neoantigens occur simultaneously, yet little is known about the interplay between responses and its impact on T cell function and tumor control. In mouse lung adenocarcinoma, we found that immunodominance is established in tumors, wherein CD8 T cell expansion is predominantly driven by the antigen that most stably binds MHC. T cells responding to subdominant antigens were enriched for a TCF1+ progenitor phenotype correlated with response to immune checkpoint blockade (ICB) therapy. However, the subdominant T cell response did not preferentially benefit from ICB due to a dysfunctional subset of TCF1+ cells marked by CCR6 and Tc17 differentiation. Analysis of human samples and sequencing datasets revealed that CCR6+ TCF1+ cells exist across human cancers and are not correlated with ICB response. Vaccination eliminated CCR6+ TCF1+ cells and dramatically improved the subdominant response, highlighting a strategy to optimally engage concurrent neoantigen responses against tumors.

Annotated expression and activity data for murine recombinase alleles and transgenes: the CrePortal resource

Recombinase alleles and transgenes can be used to facilitate spatio-temporal specificity of gene disruption or transgene expression. However, the versatility of this in vivo recombination system relies on having detailed and accurate characterization of recombinase expression and activity to enable selection of the appropriate allele or transgene. The CrePortal (http://www.informatics.jax.org/home/recombinase) leverages the informatics infrastructure of Mouse Genome Informatics to integrate data from the scientific literature, direct data submissions from the scientific community at-large, and from major projects developing new recombinase lines and characterizing recombinase expression and specificity patterns. Searching the CrePortal by recombinase activity or specific recombinase gene driver provides users with a recombinase alleles and transgenes activity tissue summary and matrix comparison of gene expression and recombinase activity with links to generation details, a recombinase activity grid, and associated phenotype annotations. Future improvements will add cell type-based activity annotations. The CrePortal provides a comprehensive presentation of recombinase allele and transgene data to assist researchers in selection of the recombinase allele or transgene based on where and when recombination is desired.

Multiple genome viewer (MGV): a new tool for visualization and comparison of multiple annotated genomes

The assembled and annotated genomes for 16 inbred mouse strains (Lilue et al., Nat Genet 50:1574–1583, 2018) and two wild-derived strains (CAROLI/EiJ and PAHARI/EiJ) (Thybert et al., Genome Res 28:448–459, 2018) are valuable resources for mouse genetics and comparative genomics. We developed the multiple genome viewer (MGV; http://www.informatics.jax.org/mgv) to support visualization, exploration, and comparison of genome annotations within and across these genomes. MGV displays chromosomal regions of user-selected genomes as horizontal tracks. Equivalent features across the genome tracks are highlighted using vertical ‘swim lane’ connectors. Navigation across the genomes is synchronized as a researcher uses the scroll and zoom functions. Researchers can generate custom sets of genes and other genome features to be displayed in MGV by entering genome coordinates, function, phenotype, disease, and/or pathway terms. MGV was developed to be genome agnostic and can be used to display homologous features across genomes of different organisms.

VDAC Genes Expression and Regulation in Mammals

VDACs are pore-forming proteins, coating the mitochondrial outer membrane, and playing the role of main regulators for metabolites exchange between cytosol and mitochondria. In mammals, three isoforms have evolutionary originated, VDAC1, VDAC2, and VDAC3. Despite similarity in sequence and structure, evidence suggests different biological roles in normal and pathological conditions for each isoform. We compared Homo sapiens and Mus musculus VDAC genes and their regulatory elements. RNA-seq transcriptome analysis shows that VDAC isoforms are expressed in human and mouse tissues at different levels with a predominance of VDAC1 and VDAC2 over VDAC3, with the exception of reproductive system. Numerous transcript variants for each isoform suggest specific context-dependent regulatory mechanisms. Analysis of VDAC core promoters has highlighted that, both in a human and a mouse, VDAC genes show features of TATA-less ones. The level of CG methylation of the human VDAC genes revealed that VDAC1 promoter is less methylated than other two isoforms. We found that expression of VDAC genes is mainly regulated by transcription factors involved in controlling cell growth, proliferation and differentiation, apoptosis, and bioenergetic metabolism. A non-canonical initiation site termed "the TCT/TOP motif," the target for translation regulation by the mTOR pathway, was identified in human VDAC2 and VDAC3 and in every murine VDACs promoter. In addition, specific TFBSs have been identified in each VDAC promoter, supporting the hypothesis that there is a partial functional divergence. These data corroborate our experimental results and reinforce the idea that gene regulation could be the key to understanding the evolutionary specialization of VDAC isoforms.

A universal framework for detecting cis-regulatory diversity in DNA regulatory regions

High-throughput sequencing-based assays measure different biochemical activities pertaining to gene regulation, genome-wide. These activities include transcription factor (TF)–DNA binding, enhancer activity, open chromatin, and more. A major goal is to understand underlying sequence components, or motifs, that can explain the measured activity. It is usually not one motif but a combination of motifs bound by cooperatively acting proteins that confers activity to such regions. Furthermore, regions can be diverse, governed by different combinations of TFs/motifs. Current approaches do not take into account this issue of combinatorial diversity. We present a new statistical framework, cisDIVERSITY, which models regions as diverse modules characterized by combinations of motifs while simultaneously learning the motifs themselves. Because cisDIVERSITY does not rely on knowledge of motifs, modules, cell type, or organism, it is general enough to be applied to regions reported by most high-throughput assays. For example, in enhancer predictions resulting from different assays—GRO-cap, STARR-seq, and those measuring chromatin structure—cisDIVERSITY discovers distinct modules and combinations of TF binding sites, some specific to the assay. From protein–DNA binding data, cisDIVERSITY identifies potential cofactors of the profiled TF, whereas from ATAC-seq data, it identifies tissue-specific regulatory modules. Finally, analysis of single-cell ATAC-seq data suggests that regions open in one cell-state encode information about future states, with certain modules staying open and others closing down in the next time point.

Thousands of induced germline mutations affecting immune cells identified by automated meiotic mapping coupled with machine learning

Forward genetic studies use meiotic mapping to adduce evidence that a particular mutation, normally induced by a germline mutagen, is causative of a particular phenotype. Particularly in small pedigrees, cosegregation of multiple mutations, occasional unawareness of mutations, and paucity of homozygotes may lead to erroneous declarations of cause and effect. We sought to improve the identification of mutations causing immune phenotypes in mice by creating Candidate Explorer (CE), a machine-learning software program that integrates 67 features of genetic mapping data into a single numeric score, mathematically convertible to the probability of verification of any putative mutation-phenotype association. At this time, CE has evaluated putative mutation-phenotype associations arising from screening damaging mutations in ∼55% of mouse genes for effects on flow cytometry measurements of immune cells in the blood. CE has therefore identified more than half of genes within which mutations can be causative of flow cytometric phenovariation in Mus musculus The majority of these genes were not previously known to support immune function or homeostasis. Mouse geneticists will find CE data informative in identifying causative mutations within quantitative trait loci, while clinical geneticists may use CE to help connect causative variants with rare heritable diseases of immunity, even in the absence of linkage information. CE displays integrated mutation, phenotype, and linkage data, and is freely available for query online.

A Multi-Strategy Sequencing Workflow in Inherited Retinal Dystrophies: Routine Diagnosis, Addressing Unsolved Cases and Candidate Genes Identification

The management of unsolved inherited retinal dystrophies (IRD) cases is challenging since no standard pipelines have been established. This study aimed to define a diagnostic algorithm useful for the diagnostic routine and to address unsolved cases. Here, we applied a Next-Generation Sequencing-based workflow, including a first step of panel sequencing (PS) followed by clinical-exome sequencing (CES) and whole-exome sequencing (WES), in 46 IRD patients belonging to 42 families. Twenty-six likely causal variants in retinal genes were found by PS and CES. CES and WES allowed proposing two novel candidate loci (WDFY3 and a X-linked region including CITED1), both abundantly expressed in human retina according to RT-PCR and immunohistochemistry. After comparison studies, PS showed the best quality and cost values, CES and WES involved similar analytical efforts and WES presented the highest diagnostic yield. These results reinforce the relevance of panels as a first step in the diagnostic routine and suggest WES as the next strategy for unsolved cases, reserving CES for the simultaneous study of multiple conditions. Standardizing this algorithm would enhance the efficiency and equity of clinical genetics practice. Furthermore, the identified candidate genes could contribute to increase the diagnostic yield and expand the mutational spectrum in these disorders.