The Mammalian Phenotype Ontology as a unifying standard for experimental and high-throughput phenotyping data

Multi-tool copy number detection highlights common body size-associated variants in miniature pig breeds from different geographical regions

BACKGROUND

Copy number variations (CNVs) represent a common and highly specific type of variation in the genome, potentially influencing genetic diversity and mammalian phenotypic development. Structural variants, such as deletions, duplications, and insertions, have frequently been highlighted as key factors influencing traits in high-production pigs. However, comprehensive CNV analyses in miniature pig breeds are limited despite their value in biomedical research.

RESULTS

This study performed whole-genome sequencing in 36 miniature pigs from nine breeds from America, Asia and Oceania, and Europe. By employing a multi-tool approach (CNVpytor, Delly, GATK gCNV, Smoove), the accuracy of CNV identification was improved. In total, 34 homozygous CNVs overlapped with exonic regions in all samples, suggesting a role in expressing specific phenotypes such as uniform growth patterns, fertility, or metabolic function. In addition, 386 copy number variation regions (CNVRs) shared by all breeds were detected, covering 33.6 Mb (1.48% of the autosomal genome). Further, 132 exclusive CNVRs were identified for American breeds, 47 for Asian and Oceanian breeds, and 114 for European breeds. Functional enrichment analysis revealed genes within the common CNVRs involved in body height determination and other growth-related parameters. Exclusive CNVRs were located in the region of genes enriched for lipid metabolism in American minipigs, reproductive traits in Asian and Oceanian breeds, and cardiovascular features and body height in European breeds. In the selected groups, quantitative trait loci associated with body size, meat quality, reproduction, and disease susceptibility were highlighted.

CONCLUSION

This investigation of the CNV landscape of minipigs underlines the impact of selective breeding on structural variants and its role in the development of specific breed phenotypes across geographical areas. The multi-tool approach provides a valuable resource for future studies on the effects of artificial selection on livestock genomes.

Improving prenatal diagnosis through standards and aggregation

Advances in sequencing and imaging technologies enable enhanced assessment in the prenatal space, with a goal to diagnose and predict the natural history of disease, to direct targeted therapies, and to implement clinical management, including transfer of care, election of supportive care, and selection of surgical interventions. The current lack of standardization and aggregation stymies variant interpretation and gene discovery, which hinders the provision of prenatal precision medicine, leaving clinicians and patients without an accurate diagnosis. With large amounts of data generated, it is imperative to establish standards for data collection, processing, and aggregation. Aggregated and homogeneously processed genetic and phenotypic data permits dissection of the genomic architecture of prenatal presentations of disease and provides a dataset on which data analysis algorithms can be tuned to the prenatal space. Here we discuss the importance of generating aggregate data sets and how the prenatal space is driving the development of interoperable standards and phenotype-driven tools.

Mapping novel QTL and fine mapping of previously identified QTL associated with glucose tolerance using the collaborative cross mice

A chronic metabolic illness, type 2 diabetes (T2D) is a polygenic and multifactorial complicated disease. With an estimated 463 million persons aged 20 to 79 having diabetes, the number is expected to rise to 700 million by 2045, creating a significant worldwide health burden. Polygenic variants of diabetes are influenced by environmental variables. T2D is regarded as a silent illness that can advance for years before being diagnosed. Finding genetic markers for T2D and metabolic syndrome in groups with similar environmental exposure is therefore essential to understanding the mechanism of such complex characteristic illnesses. So herein, we demonstrated the exclusive use of the collaborative cross (CC) mouse reference population to identify novel quantitative trait loci (QTL) and, subsequently, suggested genes associated with host glucose tolerance in response to a high-fat diet. In this study, we used 539 mice from 60 different CC lines. The diabetogenic effect in response to high-fat dietary challenge was measured by the three-hour intraperitoneal glucose tolerance test (IPGTT) test after 12 weeks of dietary challenge. Data analysis was performed using a statistical software package IBM SPSS Statistic 23. Afterward, blood glucose concentration at the specific and between different time points during the IPGTT assay and the total area under the curve (AUC0-180) of the glucose clearance was computed and utilized as a marker for the presence and severity of diabetes. The observed AUC0-180 averages for males and females were 51,267.5 and 36,537.5 mg/dL, respectively, representing a 1.4-fold difference in favor of females with lower AUC0-180 indicating adequate glucose clearance. The AUC0-180 mean differences between the sexes within each specific CC line varied widely within the CC population. A total of 46 QTL associated with the different studied phenotypes, designated as T2DSL and its number, for Type 2 Diabetes Specific Locus and its number, were identified during our study, among which 19 QTL were not previously mapped. The genomic interval of the remaining 27 QTL previously reported, were fine mapped in our study. The genomic positions of 40 of the mapped QTL overlapped (clustered) on 11 different peaks or close genomic positions, while the remaining 6 QTL were unique. Further, our study showed a complex pattern of haplotype effects of the founders, with the wild-derived strains (mainly PWK) playing a significant role in the increase of AUC values.

The Human Phenotype Ontology in 2024: phenotypes around the world

The Human Phenotype Ontology (HPO) is a widely used resource that comprehensively organizes and defines the phenotypic features of human disease, enabling computational inference and supporting genomic and phenotypic analyses through semantic similarity and machine learning algorithms. The HPO has widespread applications in clinical diagnostics and translational research, including genomic diagnostics, gene-disease discovery, and cohort analytics. In recent years, groups around the world have developed translations of the HPO from English to other languages, and the HPO browser has been internationalized, allowing users to view HPO term labels and in many cases synonyms and definitions in ten languages in addition to English. Since our last report, a total of 2239 new HPO terms and 49235 new HPO annotations were developed, many in collaboration with external groups in the fields of psychiatry, arthrogryposis, immunology and cardiology. The Medical Action Ontology (MAxO) is a new effort to model treatments and other measures taken for clinical management. Finally, the HPO consortium is contributing to efforts to integrate the HPO and the GA4GH Phenopacket Schema into electronic health records (EHRs) with the goal of more standardized and computable integration of rare disease data in EHRs.

BioKA: a curated and integrated biomarker knowledgebase for animals

Biomarkers play an important role in various area such as personalized medicine, drug development, clinical care, and molecule breeding. However, existing animals' biomarker resources predominantly focus on human diseases, leaving a significant gap in non-human animal disease understanding and breeding research. To address this limitation, we present BioKA (Biomarker Knowledgebase for Animals, https://ngdc.cncb.ac.cn/bioka), a curated and integrated knowledgebase encompassing multiple animal species, diseases/traits, and annotated resources. Currently, BioKA houses 16 296 biomarkers associated with 951 mapped diseases/traits across 31 species from 4747 references, including 11 925 gene/protein biomarkers, 1784 miRNA biomarkers, 1043 mutation biomarkers, 773 metabolic biomarkers, 357 circRNA biomarkers and 127 lncRNA biomarkers. Furthermore, BioKA integrates various annotations such as GOs, protein structures, protein-protein interaction networks, miRNA targets and so on, and constructs an interactive knowledge network of biomarkers including circRNA-miRNA-mRNA associations, lncRNA-miRNA associations and protein-protein associations, which is convenient for efficient data exploration. Moreover, BioKA provides detailed information on 308 breeds/strains of 13 species, and homologous annotations for 8784 biomarkers across 16 species, and offers three online application tools. The comprehensive knowledge provided by BioKA not only advances human disease research but also contributes to a deeper understanding of animal diseases and supports livestock breeding.

The mechanisms underlying alcohol-induced decreased splenic size: A network meta-analysis study

BACKGROUND

Organ weight change is widely accepted as a measure of toxicologic pathology. We and other groups have shown that excessive alcohol exposure leads to decreased spleen weight in rodents. This study explores the mechanisms underlying alcohol-induced splenic injury through a network meta-analysis.

METHODS

QIAGEN Ingenuity Pathway Analysis (IPA) and Mammalian Phenotype (MP) Ontology were used to identify alcohol-related molecules associated with the small spleen phenotype. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and IPA bioinformatics tools were then used to analyze the biologic processes and enriched signaling pathways engaging these molecules. In addition, the "downstream effects analysis" algorithm was used to quantify alcohol's effects.

RESULTS

IPA identified 623 molecules affected by alcohol and a Venn diagram revealed that 26 of these molecules overlapped with those associated with the MP Ontology of small spleen. The 26 molecules are TGFB1, CASP8, MTOR, ESR1, CXCR4, CAMK4, NFKBIA, DRD2, BCL2, FAS, PEBP1, TRAF2, ATM, IGHM, EDNRB, MDM2, GLRA1, PRF1, TLR7, IFNG, ALOX5, FOXO1, IL15, APOE, IKBKG, and RORA. Some of the 26 molecules were also associated with the MP Ontology of abnormal white pulp and red pulp morphology of the spleen, abnormal splenic cell ratio, decreased splenocyte number, abnormal spleen physiology, increased splenocyte apoptosis, and reduced splenocyte proliferation. STRING and IPA "Core Analysis" showed that these molecules were mainly involved in pathways related to cell apoptosis, proliferation, migration, and immune responses. IPA's "Molecular Activity Predictor" tool showed that concurrent effects of activation and inhibition of these molecules led to decreased spleen size by modulating apoptosis, proliferation, and migration of splenocytes.

CONCLUSIONS

Our network meta-analysis revealed that excessive alcohol exposure can damage the spleen through a variety of molecular mechanisms, thereby affecting immune function and human health. We found that alcohol-mediated splenic atrophy is largely mediated by increased apoptosis signaling, migration of cells, and inhibition of splenocyte proliferation.

Mouse phenome database: curated data repository with interactive multi-population and multi-trait analyses

The Mouse Phenome Database continues to serve as a curated repository and analysis suite for measured attributes of members of diverse mouse populations. The repository includes annotation to community standard ontologies and guidelines, a database of allelic states for 657 mouse strains, a collection of protocols, and analysis tools for flexible, interactive, user directed analyses that increasingly integrates data across traits and populations. The database has grown from its initial focus on a standard set of inbred strains to include heterogeneous mouse populations such as the Diversity Outbred and mapping crosses and well as Collaborative Cross, Hybrid Mouse Diversity Panel, and recombinant inbred strains. Most recently the system has expanded to include data from the International Mouse Phenotyping Consortium. Collectively these data are accessible by API and provided with an interactive tool suite that enables users' persistent selection, storage, and operation on collections of measures. The tool suite allows basic analyses, advanced functions with dynamic visualization including multi-population meta-analysis, multivariate outlier detection, trait pattern matching, correlation analyses and other functions. The data resources and analysis suite provide users a flexible environment in which to explore the basis of phenotypic variation in health and disease across the lifespan.

Biochemical pathways represented by Gene Ontology-Causal Activity Models identify distinct phenotypes resulting from mutations in pathways

Gene inactivation can affect the process(es) in which that gene acts and causally downstream ones, yielding diverse mutant phenotypes. Identifying the genetic pathways resulting in a given phenotype helps us understand how individual genes interact in a functional network. Computable representations of biological pathways include detailed process descriptions in the Reactome Knowledgebase and causal activity flows between molecular functions in Gene Ontology-Causal Activity Models (GO-CAMs). A computational process has been developed to convert Reactome pathways to GO-CAMs. Laboratory mice are widely used models of normal and pathological human processes. We have converted human Reactome GO-CAMs to orthologous mouse GO-CAMs, as a resource to transfer pathway knowledge between humans and model organisms. These mouse GO-CAMs allowed us to define sets of genes that function in a causally connected way. To demonstrate that individual variant genes from connected pathways result in similar but distinguishable phenotypes, we used the genes in our pathway models to cross-query mouse phenotype annotations in the Mouse Genome Database (MGD). Using GO-CAM representations of 2 related but distinct pathways, gluconeogenesis and glycolysis, we show that individual causal paths in gene networks give rise to discrete phenotypic outcomes resulting from perturbations of glycolytic and gluconeogenic genes. The accurate and detailed descriptions of gene interactions recovered in this analysis of well-studied processes suggest that this strategy can be applied to less well-understood processes in less well-studied model systems to predict phenotypic outcomes of novel gene variants and to identify potential gene targets in altered processes.

The accelerated organ senescence and proteotoxicity in thyrotoxicosis mice

The mechanism of aging has always been the focus of research, because aging is related to disease susceptibility and seriously affects people's quality of life. The diseases also accelerate the aging process, especially the pathological changes of substantive organs, such as cardiac hypertrophy, severely shortened lifespan. So, lesions in organs are both a consequence and a cause of aging. However, the disease in a given organ is not in isolation but is a systemic problem. Our previous study found that thyrotoxicosis mice model has aging characteristics including immunosenescence, lipotoxicity, malnutrition. But all these characteristics will lead to organ senescence, therefore, this study continued to study the aging changes of important organs such as heart, liver, and kidney in thyrotoxicosis mice using tandem mass tags (TMT) proteomics method. The results showed that the excess thyroxine led to cardiac hypertrophy. In the liver, the ability to synthesize functional proteins, detoxify, and metabolism were declined. The effect on the kidney was the decreased ability of detoxify and metabolism. The main finding of the present study was that the acceleration of organ senescence by excess thyroxine was due to proteotoxicity. The shared cause of proteotoxicity in the three organs included the intensify of oxidative phosphorylation, the redundancy production of ribosomes, and the lack of splicing and ubiquitin proteasome system function. Totally, proteotoxicity was another parallel between thyrotoxicosis and aging in addition to lipotoxicity. Our research provided a convenient and appropriate animal model for exploring aging mechanism and antiaging drugs.

Genic constraint against nonsynonymous variation across the mouse genome

BACKGROUND

Selective constraint, the depletion of variation due to negative selection, provides insights into the functional impact of variants and disease mechanisms. However, its characterization in mice, the most commonly used mammalian model, remains limited. This study aims to quantify mouse gene constraint using a new metric called the nonsynonymous observed expected ratio (NOER) and investigate its relationship with gene function.

RESULTS

NOER was calculated using whole-genome sequencing data from wild mouse populations (Mus musculus sp and Mus spretus). Positive correlations were observed between mouse gene constraint and the number of associated knockout phenotypes, indicating stronger constraint on pleiotropic genes. Furthermore, mouse gene constraint showed a positive correlation with the number of pathogenic variant sites in their human orthologues, supporting the relevance of mouse models in studying human disease variants.

CONCLUSIONS

NOER provides a resource for assessing the fitness consequences of genetic variants in mouse genes and understanding the relationship between gene constraint and function. The study's findings highlight the importance of pleiotropy in selective constraint and support the utility of mouse models in investigating human disease variants. Further research with larger sample sizes can refine constraint estimates in mice and enable more comprehensive comparisons of constraint between mouse and human orthologues.

An allelic series of spontaneous Rorb mutant mice exhibit a gait phenotype, changes in retina morphology and behavior, and gene expression signatures associated with the unfolded protein response

The Retinoid-related orphan receptor beta (RORβ) gene encodes a developmental transcription factor and has 2 predominant isoforms created through alternative first exon usage; one specific to the retina and another present more broadly in the central nervous system, particularly regions involved in sensory processing. RORβ belongs to the nuclear receptor family and plays important roles in cell fate specification in the retina and cortical layer formation. In mice, loss of RORβ causes disorganized retina layers, postnatal degeneration, and production of immature cone photoreceptors. Hyperflexion or "high-stepping" of rear limbs caused by reduced presynaptic inhibition by Rorb-expressing inhibitory interneurons of the spinal cord is evident in RORβ-deficient mice. RORβ variants in patients are associated with susceptibility to various neurodevelopmental conditions, primarily generalized epilepsies, but including intellectual disability, bipolar, and autism spectrum disorders. The mechanisms by which RORβ variants confer susceptibility to these neurodevelopmental disorders are unknown but may involve aberrant neural circuit formation and hyperexcitability during development. Here we report an allelic series in 5 strains of spontaneous Rorb mutant mice with a high-stepping gait phenotype. We show retinal abnormalities in a subset of these mutants and demonstrate significant differences in various behavioral phenotypes related to cognition. Gene expression analyses in all 5 mutants reveal a shared over-representation of the unfolded protein response and pathways related to endoplasmic reticulum stress, suggesting a possible mechanism of susceptibility relevant to patients.

Biochemical Pathways Represented by Gene Ontology Causal Activity Models Identify Distinct Phenotypes Resulting from Mutations in Pathways

Gene inactivation can affect the process(es) in which that gene acts and causally downstream ones, yielding diverse mutant phenotypes. Identifying the genetic pathways resulting in a given phenotype helps us understand how individual genes interact in a functional network. Computable representations of biological pathways include detailed process descriptions in the Reactome Knowledgebase, and causal activity flows between molecular functions in Gene Ontology-Causal Activity Models (GO-CAMs). A computational process has been developed to convert Reactome pathways to GO-CAMs. Laboratory mice are widely used models of normal and pathological human processes. We have converted human Reactome GO-CAMs to orthologous mouse GO-CAMs, as a resource to transfer pathway knowledge between humans and model organisms. These mouse GO-CAMs allowed us to define sets of genes that function in a connected and well-defined way. To test whether individual genes from well-defined pathways result in similar and distinguishable phenotypes, we used the genes in our pathway models to cross-query mouse phenotype annotations in the Mouse Genome Database (MGD). Using GO-CAM representations of two related but distinct pathways, gluconeogenesis and glycolysis, we can identify causal paths in gene networks that give rise to discrete phenotypic outcomes for perturbations of glycolysis and gluconeogenesis. The accurate and detailed descriptions of gene interactions recovered in this analysis of well-studied processes suggest that this strategy can be applied to less well-understood processes in less well-studied model systems to predict phenotypic outcomes of novel gene variants and to identify potential gene targets in altered processes.

Altered hair root gene expression profiles highlight calcium signaling and lipid metabolism pathways to be associated with curly hair initiation and maintenance in Mangalitza pigs

Hair types have been under strong targeted selection in domestic animals for their impact on skin protection, thermoregulation and exterior morphology, and subsequent economic importance. In pigs, a very special hair phenotype was observed in Mangalitza, who expresses a thick coat of curly bristles and downy hair. Two breed-specific missense variants in TRPM2 and CYP4F3 were suggested to be associated with the Mangalitza pig's hair shape due to their role in hair follicle morphogenesis reported for human and mice. However, the mechanism behind this expression of a curly hair type is still unclear and needs to be explored. In our study, hair shafts were measured and investigated for the curvature of the hair in Mangalitza and crossbreeds in comparison to straight-coated pigs. For molecular studies, hair roots underwent RNA sequencing for a differential gene expression analysis using DESeq2. The output matrix of normalized counts was then used to construct weighted gene co-expression networks. The resulting hair root gene expression profiles highlighted 454 genes to be significantly differentially expressed for initiation of curly hair phenotype in newborn Mangalitza piglets versus post-initiation in later development. Furthermore, 2,554 genes showed a significant differential gene expression in curly hair in comparison to straight hair. Neither TRPM2 nor CYP4F3 were identified as differentially expressed. Incidence of the genes in weighted co-expression networks associated with TRPM2 and CYP4F3, and prominent interactions of subsequent proteins with lipids and calcium-related pathways suggested calcium signaling and/or lipid metabolism as essential players in the induction of the curly hair as well as an ionic calcium-dependency to be a prominent factor for the maintenance of this phenotype. Subsequently, our study highlights the complex interrelations and dependencies of mutant genes TRPM2 and CYP4F3 and associated gene expression patterns, allowing the initiation of curly hair type during the development of a piglet as well as the maintenance in adult individuals.

Network modeling helps to tackle the complexity of drug-disease systems

From the (patho)physiological point of view, diseases can be considered as emergent properties of living systems stemming from the complexity of these systems. Complex systems display some typical features, including the presence of emergent behavior and the organization in successive hierarchic levels. Drug treatments increase this complexity scenario, and from some years the use of network models has been introduced to describe drug-disease systems and to make predictions about them with regard to several aspects related to drug discovery. Here, we review some recent examples thereof with the aim to illustrate how network science tools can be very effective in addressing both tasks. We will examine the use of bipartite networks that lead to the important concept of "disease module", as well as the introduction of more articulated models, like multi-scale and multiplex networks, able to describe disease systems at increasing levels of organization. Examples of predictive models will then be discussed, considering both those that exploit approaches purely based on graph theory and those that integrate machine learning methods. A short account of both kinds of methodological applications will be provided. Finally, the point will be made on the present situation of modeling complex drug-disease systems highlighting some open issues. This article is categorized under: Neurological Diseases > Computational Models Infectious Diseases > Computational Models Cardiovascular Diseases > Computational Models.

GWAS Central: an expanding resource for finding and visualising genotype and phenotype data from genome-wide association studies

The GWAS Central resource gathers and curates extensive summary-level genome-wide association study (GWAS) data and puts a range of user-friendly but powerful website tools for the comparison and visualisation of GWAS data at the fingertips of researchers. Through our continued efforts to harmonise and import data received from GWAS authors and consortia, and data sets actively collected from public sources, the database now contains over 72.5 million P-values for over 5000 studies testing over 7.4 million unique genetic markers investigating over 1700 unique phenotypes. Here, we describe an update to integrate this extensive data collection with mouse disease model data to support insights into the functional impact of human genetic variation. GWAS Central has expanded to include mouse gene-phenotype associations observed during mouse gene knockout screens. To allow similar cross-species phenotypes to be compared, terms from mammalian and human phenotype ontologies have been mapped. New interactive interfaces to find, correlate and view human and mouse genotype-phenotype associations are included in the website toolkit. Additionally, the integrated browser for interrogating multiple association data sets has been updated and a GA4GH Beacon API endpoint has been added for discovering variants tested in GWAS. The GWAS Central resource is accessible at https://www.gwascentral.org/.

Mouse Phenome Database: towards a more FAIR-compliant and TRUST-worthy data repository and tool suite for phenotypes and genotypes

The Mouse Phenome Database (MPD; https://phenome.jax.org; RRID:SCR_003212), supported by the US National Institutes of Health, is a Biomedical Data Repository listed in the Trans-NIH Biomedical Informatics Coordinating Committee registry. As an increasingly FAIR-compliant and TRUST-worthy data repository, MPD accepts phenotype and genotype data from mouse experiments and curates, organizes, integrates, archives, and distributes those data using community standards. Data are accompanied by rich metadata, including widely used ontologies and detailed protocols. Data are from all over the world and represent genetic, behavioral, morphological, and physiological disease-related characteristics in mice at baseline or those exposed to drugs or other treatments. MPD houses data from over 6000 strains and populations, representing many reproducible strain types and heterogenous populations such as the Diversity Outbred where each mouse is unique but can be genotyped throughout the genome. A suite of analysis tools is available to aggregate, visualize, and analyze these data within and across studies and populations in an increasingly traceable and reproducible manner. We have refined existing resources and developed new tools to continue to provide users with access to consistent, high-quality data that has translational relevance in a modernized infrastructure that enables interaction with a suite of bioinformatics analytic and data services.

Mouse models of NADK2 deficiency analyzed for metabolic and gene expression changes to elucidate pathophysiology

NADK2 encodes the mitochondrial form of nicotinamide adenine dinucleotide (NAD) kinase, which phosphorylates NAD. Rare recessive mutations in human NADK2 are associated with a syndromic neurological mitochondrial disease that includes metabolic changes, such as hyperlysinemia and 2,4 dienoyl CoA reductase (DECR) deficiency. However, the full pathophysiology resulting from NADK2 deficiency is not known. Here, we describe two chemically induced mouse mutations in Nadk2-S326L and S330P-which cause severe neuromuscular disease and shorten lifespan. The S330P allele was characterized in detail and shown to have marked denervation of neuromuscular junctions by 5 weeks of age and muscle atrophy by 11 weeks of age. Cerebellar Purkinje cells also showed progressive degeneration in this model. Transcriptome profiling on brain and muscle was performed at early and late disease stages. In addition, metabolomic profiling was performed on the brain, muscle, liver and spinal cord at the same ages and on plasma at 5 weeks. Combined transcriptomic and metabolomic analyses identified hyperlysinemia, DECR deficiency and generalized metabolic dysfunction in Nadk2 mutant mice, indicating relevance to the human disease. We compared findings from the Nadk model to equivalent RNA sequencing and metabolomic datasets from a mouse model of infantile neuroaxonal dystrophy, caused by recessive mutations in Pla2g6. This enabled us to identify disrupted biological processes that are common between these mouse models of neurological disease, as well as those processes that are gene-specific. These findings improve our understanding of the pathophysiology of neuromuscular diseases and describe mouse models that will be useful for future preclinical studies.

Flow blockage disrupts cilia-driven fluid transport in the epileptic brain

A carpet of ependymal motile cilia lines the brain ventricular system, forming a network of flow channels and barriers that pattern cerebrospinal fluid (CSF) flow at the surface. This CSF transport system is evolutionary conserved, but its physiological function remains unknown. Here we investigated its potential role in epilepsy with studies focused on CDKL5 deficiency disorder (CDD), a neurodevelopmental disorder with early-onset epilepsy refractory to seizure medications and the most common cause of infant epilepsy. CDKL5 is a highly conserved X-linked gene suggesting its function in regulating cilia length and motion in the green alga Chlamydomonas might have implication in the etiology of CDD. Examination of the structure and function of airway motile cilia revealed both the CDD patients and the Cdkl5 knockout mice exhibit cilia lengthening and abnormal cilia motion. Similar defects were observed for brain ventricular cilia in the Cdkl5 knockout mice. Mapping ependymal cilia generated flow in the ventral third ventricle (v3V), a brain region with important physiological functions showed altered patterning of flow. Tracing of cilia-mediated inflow into v3V with fluorescent dye revealed the appearance of a flow barrier at the inlet of v3V in Cdkl5 knockout mice. Analysis of mice with a mutation in another epilepsy-associated kinase, Yes1, showed the same disturbance of cilia motion and flow patterning. The flow barrier was also observed in the Foxj1^± and FOXJ1Cre^ERT:Cdkl5^y/fl mice, confirming the contribution of ventricular cilia to the flow disturbances. Importantly, mice exhibiting altered cilia-driven flow also showed increased susceptibility to anesthesia-induced seizure-like activity. The cilia-driven flow disturbance arises from altered cilia beating orientation with the disrupted polarity of the cilia anchoring rootlet meshwork. Together these findings indicate motile cilia disturbances have an essential role in CDD-associated seizures and beyond, suggesting cilia regulating kinases may be a therapeutic target for medication-resistant epilepsy.

How much do model organism phenotypes contribute to the computational identification of human disease genes?

Computing phenotypic similarity helps identify new disease genes and diagnose rare diseases. Genotype–phenotype data from orthologous genes in model organisms can compensate for lack of human data and increase genome coverage. In the past decade, cross-species phenotype comparisons have proven valuble, and several ontologies have been developed for this purpose. The relative contribution of different model organisms to computational identification of disease-associated genes is not fully explored. We used phenotype ontologies to semantically relate phenotypes resulting from loss-of-function mutations in model organisms to disease-associated phenotypes in humans. Semantic machine learning methods were used to measure the contribution of different model organisms to the identification of known human gene–disease associations. We found that mouse genotype–phenotype data provided the most important dataset in the identification of human disease genes by semantic similarity and machine learning over phenotype ontologies. Other model organisms' data did not improve identification over that obtained using the mouse alone, and therefore did not contribute significantly to this task. Our work impacts on the development of integrated phenotype ontologies, as well as for the use of model organism phenotypes in human genetic variant interpretation.

This article has an associated First Person interview with the first author of the paper.

Editor's choice: We investigated the use of model organism phenotypes in the computational identification of disease genes, identifying several data biases and concluding that mouse model phenotypes contribute most to computational disease gene identification.

The Vulnerability of the Developing Brain: Analysis of Highly Expressed Genes in Infant C57BL/6 Mouse Hippocampus in Relation to Phenotypic Annotation Derived From Mutational Studies

The hippocampus has been shown to have a major role in learning and memory, but also to participate in the regulation of emotions. However, its specific role(s) in memory is still unclear. Hippocampal damage or dysfunction mainly results in memory issues, especially in the declarative memory but, in animal studies, has also shown to lead to hyperactivity and difficulty in inhibiting responses previously taught. The brain structure is affected in neuropathological disorders, such as Alzheimer's, epilepsy, and schizophrenia, and also by depression and stress. The hippocampus structure is far from mature at birth and undergoes substantial development throughout infant and juvenile life. The aim of this study was to survey genes highly expressed throughout the postnatal period in mouse hippocampus and which have also been linked to an abnormal phenotype through mutational studies to achieve a greater understanding about hippocampal functions during postnatal development. Publicly available gene expression data from C57BL/6 mouse hippocampus was analyzed; from a total of 5 time points (at postnatal day 1, 10, 15, 21, and 30), 547 genes highly expressed in all of these time points were selected for analysis. Highly expressed genes are considered to be of potential biological importance and appear to be multifunctional, and hence any dysfunction in such a gene will most likely have a large impact on the development of abilities during the postnatal and juvenile period. Phenotypic annotation data downloaded from Mouse Genomic Informatics database were analyzed for these genes, and the results showed that many of them are important for proper embryo development and infant survival, proper growth, and increase in body size, as well as for voluntary movement functions, motor coordination, and balance. The results also indicated an association with seizures that have primarily been characterized by uncontrolled motor activity and the development of proper grooming abilities. The complete list of genes and their phenotypic annotation data have been compiled in a file for easy access.

Quantitative phosphoproteomics reveals GSK3A substrate network is involved in the cryodamage of sperm motility

During sperm cryopreservation, the most significant phenotype of cryodamage is the decrease in sperm motility. Several proteomics studies have already been performed to search for key regulators at the protein level. However, sperm functions are known to be highly regulated by phosphorylation signaling. Here, we constructed a quantitative phosphoproteome to investigate the expression change of phosphorylated sites during sperm cryopreservation. A total of 3107 phosphorylated sites are identified and 848 of them are found to be significantly differentially expressed (DE). Bioinformatics analysis showed that the corresponding genes of these regulated sites are highly associated with sperm motility, providing a connection between the molecular basis and the phenotype of cryodamage. We then performed kinase enrichment analysis and successfully identified glycogen synthase kinase-3α (GSK3A) as the key kinase that may play an important role in the regulation of sperm motility. We further constructed a GSK3A centric network that could help us better understand the molecular mechanism of cryodamage in sperm motility. Finally, we also verified that GSK3A was abnormally activated during this process. The presented phosphoproteome and functional associations provide abundant research resources for us to learn the regulation of sperm functions, as well as to optimize the cryoprotectant for sperm cryopreservation.

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 .

Diffusion enables integration of heterogeneous data and user-driven learning in a desktop knowledge-base

Integrating reference datasets (e.g. from high-throughput experiments) with unstructured and manually-assembled information (e.g. notes or comments from individual researchers) has the potential to tailor bioinformatic analyses to specific needs and to lead to new insights. However, developing bespoke analysis pipelines from scratch is time-consuming, and general tools for exploring such heterogeneous data are not available. We argue that by treating all data as text, a knowledge-base can accommodate a range of bioinformatic data types and applications. We show that a database coupled to nearest-neighbor algorithms can address common tasks such as gene-set analysis as well as specific tasks such as ontology translation. We further show that a mathematical transformation motivated by diffusion can be effective for exploration across heterogeneous datasets. Diffusion enables the knowledge-base to begin with a sparse query, impute more features, and find matches that would otherwise remain hidden. This can be used, for example, to map multi-modal queries consisting of gene symbols and phenotypes to descriptions of diseases. Diffusion also enables user-driven learning: when the knowledge-base cannot provide satisfactory search results in the first instance, users can improve the results in real-time by adding domain-specific knowledge. User-driven learning has implications for data management, integration, and curation.

Genetic mapping of novel modifiers for ApcMin induced intestinal polyps' development using the genetic architecture power of the collaborative cross mice

BACKGROUND

Familial adenomatous polyposis is an inherited genetic disease, characterized by colorectal polyps. It is caused by inactivating mutations in the Adenomatous polyposis coli (Apc) gene. Mice carrying a nonsense mutation in the Apc gene at R850, which is designated ApcMin/+ (Multiple intestinal neoplasia), develop intestinal adenomas. Several genetic modifier loci of Min (Mom) were previously mapped, but so far, most of the underlying genes have not been identified. To identify novel modifier loci associated with ApcMin/+, we performed quantitative trait loci (QTL) analysis for polyp development using 49 F1 crosses between different Collaborative Cross (CC) lines and C57BL/6 J-ApcMin/+mice. The CC population is a genetic reference panel of recombinant inbred lines, each line independently descended from eight genetically diverse founder strains. C57BL/6 J-ApcMin/+ males were mated with females from 49 CC lines. F1 offspring were terminated at 23 weeks and polyp counts from three sub-regions (SB1-3) of small intestinal and colon were recorded.

RESULTS

The number of polyps in all these sub-regions and colon varied significantly between the different CC lines. At 95% genome-wide significance, we mapped nine novel QTL for variation in polyp number, with distinct QTL associated with each intestinal sub-region. QTL confidence intervals varied in width between 2.63-17.79 Mb. We extracted all genes in the mapped QTL at 90 and 95% CI levels using the BioInfoMiner online platform to extract, significantly enriched pathways and key linker genes, that act as regulatory and orchestrators of the phenotypic landscape associated with the ApcMin/+ mutation.

CONCLUSIONS

Genomic structure of the CC lines has allowed us to identify novel modifiers and confirmed some of the previously mapped modifiers. Key genes involved mainly in metabolic and immunological processes were identified. Future steps in this analysis will be to identify regulatory elements - and possible epistatic effects - located in the mapped QTL.

CNVxplorer: a web tool to assist clinical interpretation of CNVs in rare disease patients

Copy Number Variants (CNVs) are an important cause of rare diseases. Array-based Comparative Genomic Hybridization tests yield a ∼12% diagnostic rate, with ∼8% of patients presenting CNVs of unknown significance. CNVs interpretation is particularly challenging on genomic regions outside of those overlapping with previously reported structural variants or disease-associated genes. Recent studies showed that a more comprehensive evaluation of CNV features, leveraging both coding and non-coding impacts, can significantly improve diagnostic rates. However, currently available CNV interpretation tools are mostly gene-centric or provide only non-interactive annotations difficult to assess in the clinical practice. Here, we present CNVxplorer, a web server suited for the functional assessment of CNVs in a clinical diagnostic setting. CNVxplorer mines a comprehensive set of clinical, genomic, and epigenomic features associated with CNVs. It provides sequence constraint metrics, impact on regulatory elements and topologically associating domains, as well as expression patterns. Analyses offered cover (a) agreement with patient phenotypes; (b) visualizations of associations among genes, regulatory elements and transcription factors; (c) enrichment on functional and pathway annotations and (d) co-occurrence of terms across PubMed publications related to the query CNVs. A flexible evaluation workflow allows dynamic re-interrogation in clinical sessions. CNVxplorer is publicly available at http://cnvxplorer.com.

Prioritizing Pain-Associated Targets with Machine Learning

While hundreds of genes have been associated with pain, much of the molecular mechanisms of pain remain unknown. As a result, current analgesics are limited to few clinically validated targets. Here, we trained a machine learning (ML) ensemble model to predict new targets for 17 categories of pain. The model utilizes features from transcriptomics, proteomics, and gene ontology to prioritize targets for modulating pain. We focused on identifying novel G-protein-coupled receptors (GPCRs), ion channels, and protein kinases because these proteins represent the most successful drug target families. The performance of the model to predict novel pain targets is 0.839 on average based on AUROC, while the predictions for arthritis had the highest accuracy (AUROC = 0.929). The model predicts hundreds of novel targets for pain; for example, GPR132 and GPR109B are highly ranked GPCRs for rheumatoid arthritis. Overall, gene-pain association predictions cluster into three groups that are enriched for cytokine, calcium, and GABA-related cell signaling pathways. These predictions can serve as a foundation for future experimental exploration to advance the development of safer and more effective analgesics.

Curation of over 10 000 transcriptomic studies to enable data reuse

Vast amounts of transcriptomic data reside in public repositories, but effective reuse remains challenging. Issues include unstructured dataset metadata, inconsistent data processing and quality control, and inconsistent probe-gene mappings across microarray technologies. Thus, extensive curation and data reprocessing are necessary prior to any reuse. The Gemma bioinformatics system was created to help address these issues. Gemma consists of a database of curated transcriptomic datasets, analytical software, a web interface and web services. Here we present an update on Gemma's holdings, data processing and analysis pipelines, our curation guidelines, and software features. As of June 2020, Gemma contains 10 811 manually curated datasets (primarily human, mouse and rat), over 395 000 samples and hundreds of curated transcriptomic platforms (both microarray and RNA sequencing). Dataset topics were represented with 10 215 distinct terms from 12 ontologies, for a total of 54 316 topic annotations (mean topics/dataset = 5.2). While Gemma has broad coverage of conditions and tissues, it captures a large majority of available brain-related datasets, accounting for 34% of its holdings. Users can access the curated data and differential expression analyses through the Gemma website, RESTful service and an R package. Database URL: https://gemma.msl.ubc.ca/home.html.

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.

The Human Phenotype Ontology in 2021

The Human Phenotype Ontology (HPO, https://hpo.jax.org) was launched in 2008 to provide a comprehensive logical standard to describe and computationally analyze phenotypic abnormalities found in human disease. The HPO is now a worldwide standard for phenotype exchange. The HPO has grown steadily since its inception due to considerable contributions from clinical experts and researchers from a diverse range of disciplines. Here, we present recent major extensions of the HPO for neurology, nephrology, immunology, pulmonology, newborn screening, and other areas. For example, the seizure subontology now reflects the International League Against Epilepsy (ILAE) guidelines and these enhancements have already shown clinical validity. We present new efforts to harmonize computational definitions of phenotypic abnormalities across the HPO and multiple phenotype ontologies used for animal models of disease. These efforts will benefit software such as Exomiser by improving the accuracy and scope of cross-species phenotype matching. The computational modeling strategy used by the HPO to define disease entities and phenotypic features and distinguish between them is explained in detail.We also report on recent efforts to translate the HPO into indigenous languages. Finally, we summarize recent advances in the use of HPO in electronic health record systems.

Container-aided integrative QTL and RNA-seq analysis of Collaborative Cross mice supports distinct sex-oriented molecular modes of response in obesity

BACKGROUND

The Collaborative Cross (CC) mouse population is a valuable resource to study the genetic basis of complex traits, such as obesity. Although the development of obesity is influenced by environmental factors, underlying genetic mechanisms play a crucial role in the response to these factors. The interplay between the genetic background and the gene expression pattern can provide further insight into this response, but we lack robust and easily reproducible workflows to integrate genomic and transcriptomic information in the CC mouse population.

RESULTS

We established an automated and reproducible integrative workflow to analyse complex traits in the CC mouse genetic reference panel at the genomic and transcriptomic levels. We implemented the analytical workflow to assess the underlying genetic mechanisms of host susceptibility to diet induced obesity and integrated these results with diet induced changes in the hepatic gene expression of susceptible and resistant mice. Hepatic gene expression differs significantly between obese and non-obese mice, with a significant sex effect, where male and female mice exhibit different responses and coping mechanisms.

CONCLUSION

Integration of the data showed that different genes but similar pathways are involved in the genetic susceptibility and disturbed in diet induced obesity. Genetic mechanisms underlying susceptibility to high-fat diet induced obesity are different in female and male mice. The clear distinction we observed in the systemic response to the high-fat diet challenge and to obesity between male and female mice points to the need for further research into distinct sex-related mechanisms in metabolic disease.

Integration of anatomy ontology data with protein-protein interaction networks improves the candidate gene prediction accuracy for anatomical entities

Background

Identification of genes responsible for anatomical entities is a major requirement in many fields including developmental biology, medicine, and agriculture. Current wet lab techniques used for this purpose, such as gene knockout, are high in resource and time consumption. Protein–protein interaction (PPI) networks are frequently used to predict disease genes for humans and gene candidates for molecular functions, but they are rarely used to predict genes for anatomical entities. Moreover, PPI networks suffer from network quality issues, which can be a limitation for their usage in predicting candidate genes. Therefore, we developed an integrative framework to improve the candidate gene prediction accuracy for anatomical entities by combining existing experimental knowledge about gene-anatomical entity relationships with PPI networks using anatomy ontology annotations. We hypothesized that this integration improves the quality of the PPI networks by reducing the number of false positive and false negative interactions and is better optimized to predict candidate genes for anatomical entities. We used existing Uberon anatomical entity annotations for zebrafish and mouse genes to construct gene networks by calculating semantic similarity between the genes. These anatomy-based gene networks were semantic networks, as they were constructed based on the anatomy ontology annotations that were obtained from the experimental data in the literature. We integrated these anatomy-based gene networks with mouse and zebrafish PPI networks retrieved from the STRING database and compared the performance of their network-based candidate gene predictions.

Results

According to evaluations of candidate gene prediction performance tested under four different semantic similarity calculation methods (Lin, Resnik, Schlicker, and Wang), the integrated networks, which were semantically improved PPI networks, showed better performances by having higher area under the curve values for receiver operating characteristic and precision-recall curves than PPI networks for both zebrafish and mouse.

Conclusion

Integration of existing experimental knowledge about gene-anatomical entity relationships with PPI networks via anatomy ontology improved the candidate gene prediction accuracy and optimized them for predicting candidate genes for anatomical entities.

Mouse Phenome Database: a data repository and analysis suite for curated primary mouse phenotype data

The Mouse Phenome Database (MPD; https://phenome.jax.org) is a widely accessed and highly functional data repository housing primary phenotype data for the laboratory mouse accessible via APIs and providing tools to analyze and visualize those data. Data come from investigators around the world and represent a broad scope of phenotyping endpoints and disease-related traits in naïve mice and those exposed to drugs, environmental agents or other treatments. MPD houses rigorously curated per-animal data with detailed protocols. Public ontologies and controlled vocabularies are used for annotation. In addition to phenotype tools, genetic analysis tools enable users to integrate and interpret genome–phenome relations across the database. Strain types and populations include inbred, recombinant inbred, F1 hybrid, transgenic, targeted mutants, chromosome substitution, Collaborative Cross, Diversity Outbred and other mapping populations. Our new analysis tools allow users to apply selected data in an integrated fashion to address problems in trait associations, reproducibility, polygenic syndrome model selection and multi-trait modeling. As we refine these tools and approaches, we will continue to provide users a means to identify consistent, quality studies that have high translational relevance.

The Monarch Initiative in 2019: an integrative data and analytic platform connecting phenotypes to genotypes across species

In biology and biomedicine, relating phenotypic outcomes with genetic variation and environmental factors remains a challenge: patient phenotypes may not match known diseases, candidate variants may be in genes that haven’t been characterized, research organisms may not recapitulate human or veterinary diseases, environmental factors affecting disease outcomes are unknown or undocumented, and many resources must be queried to find potentially significant phenotypic associations. The Monarch Initiative (https://monarchinitiative.org) integrates information on genes, variants, genotypes, phenotypes and diseases in a variety of species, and allows powerful ontology-based search. We develop many widely adopted ontologies that together enable sophisticated computational analysis, mechanistic discovery and diagnostics of Mendelian diseases. Our algorithms and tools are widely used to identify animal models of human disease through phenotypic similarity, for differential diagnostics and to facilitate translational research. Launched in 2015, Monarch has grown with regards to data (new organisms, more sources, better modeling); new API and standards; ontologies (new Mondo unified disease ontology, improvements to ontologies such as HPO and uPheno); user interface (a redesigned website); and community development. Monarch data, algorithms and tools are being used and extended by resources such as GA4GH and NCATS Translator, among others, to aid mechanistic discovery and diagnostics.

Sickle Cell Disease Ontology Working Group AdekileAdekunleAnieKofi AHamdaCherif BenBrownBiobeleBukiniDaimaCampbellAndrewChaouchMelekChimusaEmileChunda-LiyokaCatherineDennis-AntwiJemimaDerebailVimal KFlor-ParkMiriamGeardAmyGhediraKaisHaendelMelissaHanchardNeil AHotchkissJadeJonasMarioIbrahimMuntaserIngramClairInusaBabaJimohAdijat OzohuJuppSimonKamgaKarenKashimZainab AbimbolaKnight-MaddenJenniferLandouréGuidaLopez-SallPhilomeneMakaniJulieMalasaLeonardMasekoamengTshepisoMazanduGastonMnikaKhuthalaMulderNicolaMunungNchangwi SyntiaMunubeDeogratiasMwitaLiberataNembawareVictoriaNnoduObiageliOfori-AcquahSolomonOhene-FrempongKwakuOsei-AkotoAlexPaintsilVivianPanjiSumirRahimyMohamed CherifRoyalCharmaineSangedaRaphael ZTayoBamideleTiouiriInesTluwayFurahiniTreadwellMarshaTshiloloLeonVasilevskyNicoleWaiswaKasadhakawo MusaWonkamAmbroise. The Sickle Cell Disease Ontology: enabling universal sickle cell-based knowledge representation

Sickle cell disease (SCD) is one of the most common monogenic diseases in humans with multiple phenotypic expressions that can manifest as both acute and chronic complications. Although described more than a century ago, challenges in comprehensive disease management and collaborative research on this disease are compounded by the complex molecular and clinical phenotypes of SCD, environmental and psychosocial factors, limited therapeutic options and ambiguous terminology. This ambiguous terminology has hampered the integration and interoperability of existing SCD knowledge, and SCD research translation. The SCD Ontology (SCDO), which is a community-driven integrative and universal knowledge representation system for SCD, overcomes this issue by providing a controlled vocabulary developed by a group of experts in both SCD and ontology design. SCDO is the first and most comprehensive standardized human- and machine-readable resource that unambiguously represents terminology and concepts about SCD for researchers, patients and clinicians. It is built around the central concept ‘hemoglobinopathy’, allowing inclusion of non-SCD haemoglobinopathies, such as thalassaemias, which may interfere with or influence SCD phenotypic manifestations. This collaboratively developed ontology constitutes a comprehensive knowledge management system and standardized terminology of various SCD-related factors. The SCDO will promote interoperability of different research datasets, facilitate seamless data sharing and collaborations, including meta-analyses within the SCD community, and support the development and curation of data-basing and clinical informatics in SCD.

Expansion of the Human Phenotype Ontology (HPO) knowledge base and resources

The Human Phenotype Ontology (HPO)—a standardized vocabulary of phenotypic abnormalities associated with 7000+ diseases—is used by thousands of researchers, clinicians, informaticians and electronic health record systems around the world. Its detailed descriptions of clinical abnormalities and computable disease definitions have made HPO the de facto standard for deep phenotyping in the field of rare disease. The HPO’s interoperability with other ontologies has enabled it to be used to improve diagnostic accuracy by incorporating model organism data. It also plays a key role in the popular Exomiser tool, which identifies potential disease-causing variants from whole-exome or whole-genome sequencing data. Since the HPO was first introduced in 2008, its users have become both more numerous and more diverse. To meet these emerging needs, the project has added new content, language translations, mappings and computational tooling, as well as integrations with external community data. The HPO continues to collaborate with clinical adopters to improve specific areas of the ontology and extend standardized disease descriptions. The newly redesigned HPO website (www.human-phenotype-ontology.org) simplifies browsing terms and exploring clinical features, diseases, and human genes.

The mouse Gene Expression Database (GXD): 2019 update

The mouse Gene Expression Database (GXD) is an extensive, well-curated community resource freely available at www.informatics.jax.org/expression.shtml. Covering all developmental stages, GXD includes data from RNA in situ hybridization, immunohistochemistry, RT-PCR, northern blot and western blot experiments in wild-type and mutant mice. GXD's gene expression information is integrated with the other data in Mouse Genome Informatics and interconnected with other databases, placing these data in the larger biological and biomedical context. Since the last report, the ability of GXD to provide insights into the molecular mechanisms of development and disease has been greatly enhanced by the addition of new data and by the implementation of new web features. These include: improvements to the Differential Gene Expression Data Search, facilitating searches for genes that have been shown to be exclusively expressed in a specified structure and/or developmental stage; an enhanced anatomy browser that now provides access to expression data and phenotype data for a given anatomical structure; direct access to the wild-type gene expression data for the tissues affected in a specific mutant; and a comparison matrix that juxtaposes tissues where a gene is normally expressed against tissues, where mutations in that gene cause abnormalities.

Bioinformatic framework for analysis of transcription factor changes as the molecular link between replicative cellular senescence signaling pathways and carcinogenesis

Cellular senescence is a natural condition of irreversible cell cycle arrest and apoptotic resistance that occurs in cells exposed to various stress factors, such as replicative stress or overexpression of oncogenes. Unraveling the complex regulation of senescence in cells is essential to strengthen senescence-related therapeutic approaches in cancer, as cellular senescence plays a dual role in tumorigenesis, having both anti- and pro-tumorigenic effects. In our study we created a model of replicative cellular senescence, based on transcriptomic data, including an extra intermediate time-point prior to cells entering senescence, to elucidate the interplay of networks governing cellular senescence with networks involved in tumorigenesis. We reveal specific changes that occur in transcription factor activity at different timepoints before and after cells entering senescence and model the signaling networks that govern these changes.

Data-driven method to enhance craniofacial and oral phenotype vocabularies

BACKGROUND

A significant amount of clinical information captured as free-text narratives could be better used for several applications, such as clinical decision support, ontology development, evidence-based practice, and research. The Human Phenotype Ontology (HPO) is specifically used for semantic comparisons for diagnostic purposes. All these functions require quality coverage of the domain of interest. The authors used natural language processing to capture craniofacial and oral phenotype signatures from electronic health records and then used these signatures for evaluation of existing oral phenotype ontology coverage.

METHODS

The authors applied a text-processing pipeline based on the clinical Text Analysis and Knowledge Extraction System to annotate the clinical notes with Unified Medical Language System codes. The authors extracted the disease or disorder phenotype terms, which were then compared with HPO terms and their synonyms.

RESULTS

The authors retrieved 2,153 deidentified clinical notes from 558 patients. Finally, 2,416 unique diseases or disorders phenotype terms were extracted, which included 210 craniofacial or oral phenotype terms. Twenty-six of these phenotypes were not found in the HPO.

CONCLUSIONS

The authors demonstrated that natural language processing tools could extract relevant phenotype terms from clinical narratives, which could help identify gaps in existing ontologies and enhance craniofacial and dental phenotyping vocabularies.

PRACTICAL IMPLICATIONS

The expansion of terms in the dental, oral, and craniofacial domains in the HPO is particularly important as the dental community moves toward electronic health records.

Network-Based Functional Prediction Augments Genetic Association To Predict Candidate Genes for Histamine Hypersensitivity in Mice

Genetic mapping is a primary tool of genetics in model organisms; however, many quantitative trait loci (QTL) contain tens or hundreds of positional candidate genes. Prioritizing these genes for validation is often ad hoc and biased by previous findings. Here we present a technique for prioritizing positional candidates based on computationally inferred gene function. Our method uses machine learning with functional genomic networks, whose links encode functional associations among genes, to identify network-based signatures of functional association to a trait of interest. We demonstrate the method by functionally ranking positional candidates in a large locus on mouse Chr 6 (45.9 Mb to 127.8 Mb) associated with histamine hypersensitivity (Histh). Histh is characterized by systemic vascular leakage and edema in response to histamine challenge, which can lead to multiple organ failure and death. Although Histh risk is strongly influenced by genetics, little is known about its underlying molecular or genetic causes, due to genetic and physiological complexity of the trait. To dissect this complexity, we ranked genes in the Histh locus by predicting functional association with multiple Histh-related processes. We integrated these predictions with new single nucleotide polymorphism (SNP) association data derived from a survey of 23 inbred mouse strains and congenic mapping data. The top-ranked genes included Cxcl12, Ret, Cacna1c, and Cntn3, all of which had strong functional associations and were proximal to SNPs segregating with Histh. These results demonstrate the power of network-based computational methods to nominate highly plausible quantitative trait genes even in challenging cases involving large QTL and extreme trait complexity.

The JAX Synteny Browser for mouse-human comparative genomics

Visualizing regions of conserved synteny between two genomes is supported by numerous software applications. However, none of the current applications allow researchers to select genome features to display or highlight in blocks of synteny based on the annotated biological properties of the features (e.g., type, function, and/or phenotype association). To address this usability gap, we developed an interactive web-based conserved synteny browser, The Jackson Laboratory (JAX) Synteny Browser. The browser allows researchers to highlight or selectively display genome features in the reference and/or the comparison genome according to the biological attributes of the features. Although the current implementation for the browser is limited to the reference genomes for the laboratory mouse and human, the software platform is intentionally genome agnostic. The JAX Synteny Browser software can be deployed for any two genomes where genome coordinates for syntenic blocks are defined and for which biological attributes of the features in one or both genomes are available in widely used standard bioinformatics file formats. The JAX Synteny Browser is available at: http://syntenybrowser.jax.org/. The code base is available from GitHub: https://github.com/TheJacksonLaboratory/syntenybrowser and is distributed under the Creative Commons Attribution license (CC BY).

Computational resources associating diseases with genotypes, phenotypes and exposures

The causes of a disease and its therapies are not only related to genotypes, but also associated with other factors, including phenotypes, environmental exposures, drugs and chemical molecules. Distinguishing disease-related factors from many neutral factors is critical as well as difficult. Over the past two decades, bioinformaticians have developed many computational resources to integrate the omics data and discover associations among these factors. However, researchers and clinicians are experiencing difficulties in choosing appropriate resources from hundreds of relevant databases and software tools. Here, in order to assist the researchers and clinicians, we systematically review the public computational resources of human diseases related to genotypes, phenotypes, environment factors, drugs and chemical exposures. We briefly describe the development history of these computational resources, followed by the details of the relevant databases and software tools. We finally conclude with a discussion of current challenges and future opportunities as well as prospects on this topic.

Shifts in myeloarchitecture characterise adolescent development of cortical gradients

We studied an accelerated longitudinal cohort of adolescents and young adults (n = 234, two time points) to investigate dynamic reconfigurations in myeloarchitecture. Intracortical profiles were generated using magnetization transfer (MT) data, a myelin-sensitive magnetic resonance imaging contrast. Mixed-effect models of depth specific intracortical profiles demonstrated two separate processes i) overall increases in MT, and ii) flattening of the MT profile related to enhanced signal in mid-to-deeper layers, especially in heteromodal and unimodal association cortices. This development was independent of morphological changes. Enhanced MT in mid-to-deeper layers was found to spatially co-localise specifically with gene expression markers of oligodendrocytes. Interregional covariance analysis revealed that these intracortical changes contributed to a gradual differentiation of higher-order from lower-order systems. Depth-dependent trajectories of intracortical myeloarchitectural development contribute to the maturation of structural hierarchies in the human neocortex, providing a model for adolescent development that bridges microstructural and macroscopic scales of brain organisation.

Potentially repurposable drugs for schizophrenia identified from its interactome

We previously presented the protein-protein interaction network of schizophrenia associated genes, and from it, the drug-protein interactome which showed the drugs that target any of the proteins in the interactome. Here, we studied these drugs further to identify whether any of them may potentially be repurposable for schizophrenia. In schizophrenia, gene expression has been described as a measurable aspect of the disease reflecting the action of risk genes. We studied each of the drugs from the interactome using the BaseSpace Correlation Engine, and shortlisted those that had a negative correlation with differential gene expression of schizophrenia. This analysis resulted in 12 drugs whose differential gene expression (drug versus normal) had an anti-correlation with differential expression for schizophrenia (disorder versus normal). Some of these drugs were already being tested for their clinical activity in schizophrenia and other neuropsychiatric disorders. Several proteins in the protein interactome of the targets of several of these drugs were associated with various neuropsychiatric disorders. The network of genes with opposite drug-induced versus schizophrenia-associated expression profiles were significantly enriched in pathways relevant to schizophrenia etiology and GWAS genes associated with traits or diseases that had a pathophysiological overlap with schizophrenia. Drugs that targeted the same genes as the shortlisted drugs, have also demonstrated clinical activity in schizophrenia and other related disorders. This integrated computational analysis will help translate insights from the schizophrenia drug-protein interactome to clinical research - an important step, especially in the field of psychiatric drug development which faces a high failure rate.

Mouse Phenome Database: an integrative database and analysis suite for curated empirical phenotype data from laboratory mice

The Mouse Phenome Database (MPD; https://phenome.jax.org) is a widely used resource that provides access to primary experimental trait data, genotypic variation, protocols and analysis tools for mouse genetic studies. Data are contributed by investigators worldwide and represent a broad scope of phenotyping endpoints and disease-related traits in naïve mice and those exposed to drugs, environmental agents or other treatments. MPD houses individual animal data with detailed, searchable protocols, and makes these data available to other resources via API. MPD provides rigorous curation of experimental data and supporting documentation using relevant ontologies and controlled vocabularies. Most data in MPD are from inbreds and other reproducible strains such that the data are cumulative over time and across laboratories. The resource has been expanded to include the QTL Archive and other primary phenotype data from mapping crosses as well as advanced high-diversity mouse populations including the Collaborative Cross and Diversity Outbred mice. Furthermore, MPD provides a means of assessing replicability and reproducibility across experimental conditions and protocols, benchmarking assays in users’ own laboratories, identifying sensitized backgrounds for making new mouse models with genome editing technologies, analyzing trait co-inheritance, finding the common genetic basis for multiple traits and assessing sex differences and sex-by-genotype interactions.

Mouse Genome Database (MGD)-2018: knowledgebase for the laboratory mouse

The Mouse Genome Database (MGD; http://www.informatics.jax.org) is the key community mouse database which supports basic, translational and computational research by providing integrated data on the genetics, genomics, and biology of the laboratory mouse. MGD serves as the source for biological reference data sets related to mouse genes, gene functions, phenotypes and disease models with an increasing emphasis on the association of these data to human biology and disease. We report here on recent enhancements to this resource, including improved access to mouse disease model and human phenotype data and enhanced relationships of mouse models to human disease.

eRAM: encyclopedia of rare disease annotations for precision medicine

Rare diseases affect over a hundred million people worldwide, most of these patients are not accurately diagnosed and effectively treated. The limited knowledge of rare diseases forms the biggest obstacle for improving their treatment. Detailed clinical phenotyping is considered as a keystone of deciphering genes and realizing the precision medicine for rare diseases. Here, we preset a standardized system for various types of rare diseases, called encyclopedia of Rare disease Annotations for Precision Medicine (eRAM). eRAM was built by text-mining nearly 10 million scientific publications and electronic medical records, and integrating various data in existing recognized databases (such as Unified Medical Language System (UMLS), Human Phenotype Ontology, Orphanet, OMIM, GWAS). eRAM systematically incorporates currently available data on clinical manifestations and molecular mechanisms of rare diseases and uncovers many novel associations among diseases. eRAM provides enriched annotations for 15 942 rare diseases, yielding 6147 human disease related phenotype terms, 31 661 mammalians phenotype terms, 10,202 symptoms from UMLS, 18 815 genes and 92 580 genotypes. eRAM can not only provide information about rare disease mechanism but also facilitate clinicians to make accurate diagnostic and therapeutic decisions towards rare diseases. eRAM can be freely accessed at http://www.unimd.org/eram/.

Genetic mapping in Diversity Outbred mice identifies a Trpa1 variant influencing late-phase formalin response

Identification of genetic variants that influence susceptibility to pain is key to identifying molecular mechanisms and targets for effective and safe therapeutic alternatives to opioids. To identify genes and variants associated with persistent pain, we measured late-phase response to formalin injection in 275 male and female Diversity Outbred mice genotyped for over 70,000 single nucleotide polymorphisms. One quantitative trait locus reached genome-wide significance on chromosome 1 with a support interval of 3.1 Mb. This locus, Nociq4 (nociceptive sensitivity quantitative trait locus 4; MGI: 5661503), harbors the well-known pain gene Trpa1 (transient receptor potential cation channel, subfamily A, member 1). Trpa1 is a cation channel known to play an important role in acute and chronic pain in both humans and mice. Analysis of Diversity Outbred founder strain allele effects revealed a significant effect of the CAST/EiJ allele at Trpa1, with CAST/EiJ carrier mice showing an early, but not late, response to formalin relative to carriers of the 7 other inbred founder alleles (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, PWK/PhJ, and WSB/EiJ). We characterized possible functional consequences of sequence variants in Trpa1 by assessing channel conductance, TRPA1-TRPV1 interactions, and isoform expression. The phenotypic differences observed in CAST/EiJ relative to C57BL/6J carriers were best explained by Trpa1 isoform expression differences, implicating a splice junction variant as the causal functional variant. This study demonstrates the utility of advanced, high-precision genetic mapping populations in resolving specific molecular mechanisms of variation in pain sensitivity.

Curating gene sets: challenges and opportunities for integrative analysis

Genomic data interpretation often requires analyses that move from a gene-by-gene focus to a focus on sets of genes that are associated with biological phenomena such as molecular processes, phenotypes, diseases, drug interactions or environmental conditions. Unique challenges exist in the curation of gene sets beyond the challenges in curation of individual genes. Here we highlight a literature curation workflow whereby gene sets are curated from peer-reviewed published data into GeneWeaver (GW), a data repository and analysis platform. We describe the system features that allow for a flexible yet precise curation procedure. We illustrate the value of curation by gene sets through analysis of independently curated sets that relate to the integrated stress response, showing that sets curated from independent sources all share significant Jaccard similarity. A suite of reproducible analysis tools is provided in GW as services to carry out interactive functional investigation of user-submitted gene sets within the context of over 150 000 gene sets constructed from publicly available resources and published gene lists. A curation interface supports the ability of users to design and maintain curation workflows of gene sets, including assigning, reviewing and releasing gene sets within a curation project context.

Quantitative phenotype analysis to identify, validate and compare rat disease models

The laboratory rat has been widely used as an animal model in biomedical research. There are many strains exhibiting a wide variety of phenotypes. Capturing these phenotypes in a centralized database provides researchers with an easy method for choosing the appropriate strains for their studies. Existing resources have provided some preliminary work in rat phenotype databases. However, existing resources suffer from problems such as small number of animals, lack of updating, web interface queries limitations and lack of standardized metadata. The Rat Genome Database (RGD) PhenoMiner tool has provided the first step in this effort by standardizing and integrating data from individual studies. Our work, mainly utilizing data curated in RGD, involves the following key steps: (i) we developed a meta-analysis pipeline to automatically integrate data from heterogeneous sources and to produce expected ranges (standardized phenotype ranges) for different strains and phenotypes under different experimental conditions; (ii) we created tools to visualize expected ranges for individual strains and strain groups. We developed a meta-analysis pipeline and an interactive web interface that summarizes and visualizes expected ranges produced from the meta-analysis pipeline. Automation of the pipeline allows for updates as additional data becomes available. The interactive web interface provides curators and researchers with a platform for identifying and validating expected ranges for a variety of quantitative phenotypes. The data analysis result and visualization tools will promote an understanding of rat disease models, guide researchers to choose optimal strains for their research needs and encourage data sharing from different research hubs. Such resources also help to promote research reproducibility. The interactive platforms created in this project will continue to provide a valuable resource for translational research efforts.