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.

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.

Lysophospholipid (LPA) receptors in GtoPdb v.2021.3

Lysophosphatidic acid (LPA) receptors (nomenclature as agreed by the NC-IUPHAR Subcommittee on Lysophospholipid Receptors [55, 19, 82, 129]) are activated by the endogenous phospholipid LPA. The first receptor, LPA1, was identified as ventricular zone gene-1 (vzg-1) [40], This discovery represented the beginning of the de-orphanisation of members of the endothelial differentiation gene (edg) family, as other LPA and sphingosine 1-phosphate (S1P) receptors were found. Five additional LPA receptors (LPA2,3,4,5,6) have since been identified [82] and their gene nomenclature codified for human LPAR1, LPAR2, etc. (HUGO Gene Nomenclature Committee, HGNC) and Lpar1, Lpar2, etc. for mice (Mouse Genome Informatics Database, MGI) to reflect species and receptor function of their corresponding proteins. The crystal structure of LPA1 is solved and indicates that LPA accesses the extracellular binding pocket, consistent with its proposed delivery via autotaxin [13]. These studies have also implicated cross-talk with endocannabinoids via phosphorylated intermediates that can also activate these receptors. The binding affinities to LPA1 of unlabeled, natural LPA and anandamide phosphate (AEAp) were measured using backscattering interferometry (pKd = 9) [83, 104]. Utilization of this method indicated affinities that were 77-fold lower than when measured using radioactivity-based protocols [128]. Targeted deletion of LPA receptors has clarified signalling pathways and identified physiological and pathophysiological roles. Multiple groups have independently published validation of all six LPA receptors described in these tables, and further validation was achieved using a distinct read-out via a novel TGFα "shedding* assay [48]. LPA has been proposed to be a ligand for GPR35 [94], supported by a study revealing that LPA modulates macrophage function through GPR35 [54]. However chemokine (C-X-C motif) ligand 17 (CXCL17) is reported to be a ligand for GPR35/CXCR8 [76]. Moreover, LPA has also been described as an agonist for the transient receptor potential (Trp) ion channels TRPV1 [87] and TRPA1 [58]. All of these proposed non-GPCR receptor identities require confirmation and are not currently recognized as bona fide LPA receptors.

Murine allele and transgene symbols: ensuring unique, concise, and informative nomenclature

In addition to naturally occurring sequence variation and spontaneous mutations, a wide array of technologies exist for modifying the mouse genome. Standardized nomenclature, including allele, transgene, and other mutation nomenclature, as well as persistent unique identifiers (PUID) are critical for effective scientific communication, comparison of results, and integration of data into knowledgebases such as Mouse Genome Informatics (MGI), Alliance for Genome Resources, and International Mouse Strain Resource (IMSR). As well as being the authoritative source for mouse gene, allele, and strain nomenclature, MGI integrates published and unpublished genomic, phenotypic, and expression data while linking to other online resources for a complete view of the mouse as a valuable model organism. The International Committee on Standardized Genetic Nomenclature for Mice has developed allele nomenclature rules and guidelines that take into account the number of genes impacted, the method of allele generation, and the nature of the sequence alteration. To capture details that cannot be included in allele symbols, MGI has further developed allele to gene relationships using sequence ontology (SO) definitions for mutations that provide links between alleles and the genes affected. MGI is also using (HGVS) variant nomenclature for variants associated with alleles that will enhance searching for mutations and will improve cross-species comparison. With the ability to assign unique and informative symbols as well as to link alleles with more than one gene, allele and transgene nomenclature rules and guidelines provide an unambiguous way to represent alterations in the mouse genome and facilitate data integration among multiple resources such the Alliance of Genome Resources and International Mouse Strain Resource.

Lysophospholipid (LPA) receptors in GtoPdb v.2021.2

Lysophosphatidic acid (LPA) receptors (nomenclature as agreed by the NC-IUPHAR Subcommittee on Lysophospholipid Receptors [55, 19, 82, 129]) are activated by the endogenous phospholipid LPA. The first receptor, LPA1, was identified as ventricular zone gene-1 (vzg-1) [40], This discovery represented the beginning of the de-orphanisation of members of the endothelial differentiation gene (edg) family, as other LPA and sphingosine 1-phosphate (S1P) receptors were found. Five additional LPA receptors (LPA2,3,4,5,6) have since been identified [82] and their gene nomenclature codified for human LPAR1, LPAR2, etc. (HUGO Gene Nomenclature Committee, HGNC) and Lpar1, Lpar2, etc. for mice (Mouse Genome Informatics Database, MGI) to reflect species and receptor function of their corresponding proteins. The crystal structure of LPA1 is solved and indicates that LPA accesses the extracellular binding pocket, consistent with its proposed delivery via autotaxin [13]. These studies have also implicated cross-talk with endocannabinoids via phosphorylated intermediates that can also activate these receptors. The binding affinities to LPA1 of unlabeled, natural LPA and anandamide phosphate (AEAp) were measured using backscattering interferometry (pKd = 9) [83, 104]. Utilization of this method indicated affinities that were 77-fold lower than when measured using radioactivity-based protocols [128]. Targeted deletion of LPA receptors has clarified signalling pathways and identified physiological and pathophysiological roles. Multiple groups have independently published validation of all six LPA receptors described in these tables, and further validation was achieved using a distinct read-out via a novel TGFα "shedding* assay [48]. LPA LPA has been proposed to be a ligand for GPCR35 [94], supported by a recent study revealing that LPA modulates macrophage function through GPR35 [54]. However chemokine (C-X-C motif) ligand 17 (CXCL17) is reported to be a ligand for GPR35/CXCR8 [76]. Moreover, LPA has also been described as an agonist for the transient receptor potential (Trp) ion channels TRPV1 [87] and TRPA1 [58]. All of these proposed non-GPCR receptor identities require confirmation and are not currently recognized as bona fide LPA receptors.

Lysophospholipid (S1P) receptors in GtoPdb v.2021.2

Sphingosine 1-phosphate (S1P) receptors (nomenclature as agreed by the NC-IUPHAR Subcommittee on Lysophospholipid receptors [89]) are activated by the endogenous lipid sphingosine 1-phosphate (S1P). Originally cloned as orphan members of the endothelial differentiation gene (edg) family [16, 112], the receptors are currently designated as S1P1R through S1P5R [69, 16, 112]. Their gene nomenclature has been codified as human S1PR1, S1PR2, etc. (HUGO Gene Nomenclature Committee, HGNC) and S1pr1, S1pr2, etc. for mice (Mouse Genome Informatics Database, MGI) to reflect species and receptor function. All S1P receptors have been knocked-out in mice constitutively and in some cases, conditionally. S1PRs, particularly S1P1, are expressed throughout all mammalian organ systems. Ligand delivery occurs via two known carriers (or "chaperones"): albumin and HDL-bound apolipoprotein M (ApoM), the latter of which elicits biased agonist signaling by S1P1 in multiple cell types [18, 49]. The five S1PRs, two chaperones, and active cellular metabolism have complicated analyses of receptor ligand binding in native systems. Signaling pathways and physiological roles have been characterized through radioligand binding in heterologous expression systems, targeted deletion of the different S1PRs, and most recently, mouse models that report in vivo S1P1R activation [94, 96]. A crystal structure of an S1P1-T4 fusion protein confirmed aspects of ligand binding, specificity, and receptor activation, determined previously through biochemical and genetic studies [65, 17]. fingolimod (FTY720), the first FDA-approved drug to target any of the lysophospholipid receptors, binds as a phosphorylated metabolite to four of the five S1PRs, and was the first oral therapy for multiple sclerosis (MS) [33]. siponimod and ozanimod that target S1P1 and S1P5 are also FDA approved for the treatment of various MS forms [16, 112]. The mechanisms of action of fingolimod and other S1PR-modulating drugs now in development include binding S1PRs in multiple organ systems, e.g., immune and nervous systems, although the precise nature of their receptor interactions requires clarification [129, 35, 59, 60].

PICKLE 3.0: enriching the human meta-database with the mouse protein interactome extended via mouse–human orthology

Summary The PICKLE 3.0 upgrade refers to the enrichment of this human protein–protein interaction (PPI) meta-database with the mouse protein interactome. Experimental PPI data between mouse genetic entities are rather limited; however, they are substantially complemented by PPIs between mouse and human genetic entities. The relational scheme of PICKLE 3.0 has been amended to exploit the Mouse Genome Informatics mouse–human ortholog gene pair collection, enabling (i) the extension through orthology of the mouse interactome with potentially valid PPIs between mouse entities based on the experimental PPIs between mouse and human entities and (ii) the comparison between mouse and human PPI networks. Interestingly, 43.5% of the experimental mouse PPIs lacks a corresponding by orthology PPI in human, an inconsistency in need of further investigation. Overall, as primary mouse PPI datasets show a considerably limited overlap, PICKLE 3.0 provides a unique comprehensive representation of the mouse protein interactome. Availability and implementation PICKLE can be queried and downloaded at http://www.pickle.gr. Supplementary information Supplementary data are available at Bioinformatics online.

HGNChelper: identification and correction of invalid gene symbols for human and mouse

Gene symbols are recognizable identifiers for gene names but are unstable and error-prone due to aliasing, manual entry, and unintentional conversion by spreadsheets to date format. Official gene symbol resources such as HUGO Gene Nomenclature Committee (HGNC) for human genes and the Mouse Genome Informatics project (MGI) for mouse genes provide authoritative sources of valid, aliased, and outdated symbols, but lack a programmatic interface and correction of symbols converted by spreadsheets. We present HGNChelper, an R package that identifies known aliases and outdated gene symbols based on the HGNC human and MGI mouse gene symbol databases, in addition to common mislabeling introduced by spreadsheets, and provides corrections where possible. HGNChelper identified invalid gene symbols in the most recent Molecular Signatures Database (mSigDB 7.0) and in platform annotation files of the Gene Expression Omnibus, with prevalence ranging from ~3% in recent platforms to 30-40% in the earliest platforms from 2002-03. HGNChelper is installable from CRAN.

The mouse Gene Expression Database (GXD): 2021 update

The Gene Expression Database (GXD; www.informatics.jax.org/expression.shtml) is an extensive and well-curated community resource of mouse developmental gene expression information. For many years, GXD has collected and integrated data from RNA in situ hybridization, immunohistochemistry, RT-PCR, northern blot, and western blot experiments through curation of the scientific literature and by collaborations with large-scale expression projects. Since our last report in 2019, we have continued to acquire these classical types of expression data; developed a searchable index of RNA-Seq and microarray experiments that allows users to quickly and reliably find specific mouse expression studies in ArrayExpress (https://www.ebi.ac.uk/arrayexpress/) and GEO (https://www.ncbi.nlm.nih.gov/geo/); and expanded GXD to include RNA-Seq data. Uniformly processed RNA-Seq data are imported from the EBI Expression Atlas and then integrated with the other types of expression data in GXD, and with the genetic, functional, phenotypic and disease-related information in Mouse Genome Informatics (MGI). This integration has made the RNA-Seq data accessible via GXD’s enhanced searching and filtering capabilities. Further, we have embedded the Morpheus heat map utility into the GXD user interface to provide additional tools for display and analysis of RNA-Seq data, including heat map visualization, sorting, filtering, hierarchical clustering, nearest neighbors analysis and visual enrichment.

HGNChelper: identification and correction of invalid gene symbols for human and mouse

Gene symbols are recognizable identifiers for gene names but are unstable and error-prone due to aliasing, manual entry, and unintentional conversion by spreadsheets to date format. Official gene symbol resources such as HUGO Gene Nomenclature Committee (HGNC) for human genes and the Mouse Genome Informatics project (MGI) for mouse genes provide authoritative sources of valid, aliased, and outdated symbols, but lack a programmatic interface and correction of symbols converted by spreadsheets. We present HGNChelper, an R package that identifies known aliases and outdated gene symbols based on the HGNC human and MGI mouse gene symbol databases, in addition to common mislabeling introduced by spreadsheets, and provides corrections where possible. HGNChelper identified invalid gene symbols in the most recent Molecular Signatures Database (mSigDB 7.0) and in platform annotation files of the Gene Expression Omnibus, with prevalence ranging from ∼3% in recent platforms to 30-40% in the earliest platforms from 2002-03. HGNChelper is installable from CRAN, with open development and issue tracking on GitHub and an associated pkgdown site https://waldronlab.io/HGNChelper/.