The Year of the Rat: The Rat Genome Database at 20: a multi-species knowledgebase and analysis platform

The impact of exercise on gene regulation in association with complex trait genetics

Endurance exercise training is known to reduce risk for a range of complex diseases. However, the molecular basis of this effect has been challenging to study and largely restricted to analyses of either few or easily biopsied tissues. Extensive transcriptome data collected across 15 tissues during exercise training in rats as part of the Molecular Transducers of Physical Activity Consortium has provided a unique opportunity to clarify how exercise can affect tissue-specific gene expression and further suggest how exercise adaptation may impact complex disease-associated genes. To build this map, we integrate this multi-tissue atlas of gene expression changes with gene-disease targets, genetic regulation of expression, and trait relationship data in humans. Consensus from multiple approaches prioritizes specific tissues and genes where endurance exercise impacts disease-relevant gene expression. Specifically, we identify a total of 5523 trait-tissue-gene triplets to serve as a valuable starting point for future investigations [Exercise; Transcription; Human Phenotypic Variation].

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

Interaction models matter: an efficient, flexible computational framework for model-specific investigation of epistasis

PURPOSE

Epistasis, the interaction between two or more genes, is integral to the study of genetics and is present throughout nature. Yet, it is seldom fully explored as most approaches primarily focus on single-locus effects, partly because analyzing all pairwise and higher-order interactions requires significant computational resources. Furthermore, existing methods for epistasis detection only consider a Cartesian (multiplicative) model for interaction terms. This is likely limiting as epistatic interactions can evolve to produce varied relationships between genetic loci, some complex and not linearly separable.

METHODS

We present new algorithms for the interaction coefficients for standard regression models for epistasis that permit many varied models for the interaction terms for loci and efficient memory usage. The algorithms are given for two-way and three-way epistasis and may be generalized to higher order epistasis. Statistical tests for the interaction coefficients are also provided. We also present an efficient matrix based algorithm for permutation testing for two-way epistasis. We offer a proof and experimental evidence that methods that look for epistasis only at loci that have main effects may not be justified. Given the computational efficiency of the algorithm, we applied the method to a rat data set and mouse data set, with at least 10,000 loci and 1,000 samples each, using the standard Cartesian model and the XOR model to explore body mass index.

RESULTS

This study reveals that although many of the loci found to exhibit significant statistical epistasis overlap between models in rats, the pairs are mostly distinct. Further, the XOR model found greater evidence for statistical epistasis in many more pairs of loci in both data sets with almost all significant epistasis in mice identified using XOR. In the rat data set, loci involved in epistasis under the XOR model are enriched for biologically relevant pathways.

CONCLUSION

Our results in both species show that many biologically relevant epistatic relationships would have been undetected if only one interaction model was applied, providing evidence that varied interaction models should be implemented to explore epistatic interactions that occur in living systems.

The Born in Guangzhou Cohort Study enables generational genetic discoveries

Genomic research that targets large-scale, prospective birth cohorts constitutes an essential strategy for understanding the influence of genetics and environment on human health1. Nonetheless, such studies remain scarce, particularly in Asia. Here we present the phase I genome study of the Born in Guangzhou Cohort Study2 (BIGCS), which encompasses the sequencing and analysis of 4,053 Chinese individuals, primarily composed of trios or mother-infant duos residing in South China. Our analysis reveals novel genetic variants, a high-quality reference panel, and fine-scale local genetic structure within BIGCS. Notably, we identify previously unreported East Asian-specific genetic associations with maternal total bile acid, gestational weight gain and infant cord blood traits. Additionally, we observe prevalent age-specific genetic effects on lipid levels in mothers and infants. In an exploratory intergenerational Mendelian randomization analysis, we estimate the maternal putatively causal and fetal genetic effects of seven adult phenotypes on seven fetal growth-related measurements. These findings illuminate the genetic links between maternal and early-life traits in an East Asian population and lay the groundwork for future research into the intricate interplay of genetics, intrauterine exposures and early-life experiences in shaping long-term health.

Genetic analysis of the blood transcriptome of young healthy pigs to improve disease resilience

BACKGROUND

Disease resilience is the ability of an animal to maintain productive performance under disease conditions and is an important selection target. In pig breeding programs, disease resilience must be evaluated on selection candidates without exposing them to disease. To identify potential genetic indicators for disease resilience that can be measured on selection candidates, we focused on the blood transcriptome of 1594 young healthy pigs with subsequent records on disease resilience. Transcriptome data were obtained by 3'mRNA sequencing and genotype data were from a 650 K genotyping array.

RESULTS

Heritabilities of the expression of 16,545 genes were estimated, of which 5665 genes showed significant estimates of heritability (p < 0.05), ranging from 0.05 to 0.90, with or without accounting for white blood cell composition. Genes with heritable expression levels were spread across chromosomes, but were enriched in the swine leukocyte antigen region (average estimate > 0.2). The correlation of heritability estimates with the corresponding estimates obtained for genes expressed in human blood was weak but a sizable number of genes with heritable expression levels overlapped. Genes with heritable expression levels were significantly enriched for biological processes such as cell activation, immune system process, stress response, and leukocyte activation, and were involved in various disease annotations such as RNA virus infection, including SARS-Cov2, as well as liver disease, and inflammation. To estimate genetic correlations with disease resilience, 3205 genotyped pigs, including the 1594 pigs with transcriptome data, were evaluated for disease resilience following their exposure to a natural polymicrobial disease challenge. Significant genetic correlations (p < 0.05) were observed with all resilience phenotypes, although few exceeded expected false discovery rates. Enrichment analysis of genes ranked by estimates of genetic correlations with resilience phenotypes revealed significance for biological processes such as regulation of cytokines, including interleukins and interferons, and chaperone mediated protein folding.

CONCLUSIONS

These results suggest that expression levels in the blood of young healthy pigs for genes in biological pathways related to immunity and endoplasmic reticulum stress have potential to be used as genetic indicator traits to select for disease resilience.

How is Big Data reshaping preclinical aging research?

The exponential scientific and technological progress during the past 30 years has favored the comprehensive characterization of aging processes with their multivariate nature, leading to the advent of Big Data in preclinical aging research. Spanning from molecular omics to organism-level deep phenotyping, Big Data demands large computational resources for storage and analysis, as well as new analytical tools and conceptual frameworks to gain novel insights leading to discovery. Systems biology has emerged as a paradigm that utilizes Big Data to gain insightful information enabling a better understanding of living organisms, visualized as multilayered networks of interacting molecules, cells, tissues and organs at different spatiotemporal scales. In this framework, where aging, health and disease represent emergent states from an evolving dynamic complex system, context given by, for example, strain, sex and feeding times, becomes paramount for defining the biological trajectory of an organism. Using bioinformatics and artificial intelligence, the systems biology approach is leading to remarkable advances in our understanding of the underlying mechanism of aging biology and assisting in creative experimental study designs in animal models. Future in-depth knowledge acquisition will depend on the ability to fully integrate information from different spatiotemporal scales in organisms, which will probably require the adoption of theories and methods from the field of complex systems. Here we review state-of-the-art approaches in preclinical research, with a focus on rodent models, that are leading to conceptual and/or technical advances in leveraging Big Data to understand basic aging biology and its full translational potential.

Omic horizon expression: a database of gene expression based on RNA sequencing data

BACKGROUND

Gene expression profiles have important significance for gene expression characteristics and further functional studies. More attention has been given to the expression databases in humans and mice, but less attention has been given to rats, while rat models also play an irreplaceable role in biomedical experiments.

RESULTS

To depict the rat gene expression profiles in mRNA expression levels, we analyzed over 2,700 RNA sequencing (RNA-Seq) samples from 48 tissues, 40 primary cell types and 25 cell lines; and then mapped them to the latest version of the rat genome reference, mRatBN7.2. Based on these datasets and reanalysis, we constructed a new database, the Omic Horizon Expression Database ( http://immudb.bjmu.edu.cn/expression.html ), which allows expressional profile query of over 25,000 rat genes based on non-redundant gene symbols. The database supports requests using gene symbols (or alias), Ensemble and Entrez gene IDs. Gene expression profiles can be queried in three categories: tissues, primary cells and cell lines. Application examples including expression profiling and comparison, as well as identification of novel rat genes, were illustrated to show the utility of the database.

CONCLUSIONS

As an omic resource, the Omic Horizon Expression Database provides horizons of gene expression profiles across various tissues and cells, which greatly facilitates the identification of rat genes as well as functional clues.

Targeting microRNA-145-mediated progressive phenotypes of early bladder cancer in a molecularly defined in vivo model

A progressive subclass of early-stage non-muscle-invasive bladder cancer (NMIBC) frequently recurs and progress into invasive carcinoma, thus decreasing the overall survival rate of NMIBC. However, therapeutic development for progressive NMIBC has been challenging due to the lack of molecularly validated in vivo models and agents targeting its genetic vulnerability. We herein molecularly characterized an interventional model of progressive NMIBC and revealed the principal functions and therapeutic potential of microRNA-145 (miR-145) in early bladder tumorigenesis. N-butyl-N-(4-hydroxybutyl)nitrosamine-induced premalignant lesions (BiPLs) in rats exhibited downregulated expression of miR-145 as well as highly similar mutation/expression profiles to those of the human progressive NMIBC subclass with the worst prognosis. The expression patterns of miR-145 inversely correlated with those of BC-related oncogenes in BiPLs. We also demonstrated that miR-145 dominantly regulated interferon pathways and c-Myc expression, which play a crucial role in the pathogenesis of progressive NMIBC. Furthermore, we demonstrated that miR-145 replacement with a novel miR-145-based intravesical agent (miR-145S1) significantly inhibited the progression of BiPLs in vivo. These results provide insights into the essential role of miR-145 as the earliest-acting oncogenic driver of bladder tumorigenesis as well as a validated interventional model and novel miR-145-based nucleic acid therapeutic agent for progressive NMIBC.

Rare disease research resources at the Rat Genome Database

Rare diseases individually affect relatively few people, but as a group they impact considerable numbers of people. The Rat Genome Database (https://rgd.mcw.edu) is a knowledgebase that offers resources for rare disease research. This includes disease definitions, genes, quantitative trail loci (QTLs), genetic variants, annotations to published literature, links to external resources, and more. One important resource is identifying relevant cell lines and rat strains that serve as models for disease research. Diseases, genes, and strains have report pages with consolidated data, and links to analysis tools. Utilizing these globally accessible resources for rare disease research, potentiating discovery of mechanisms and new treatments, can point researchers toward solutions to alleviate the suffering of those afflicted with these diseases.

The Rat Genome Database: Genetic, Genomic, and Phenotypic Data Across Multiple Species

The laboratory rat, Rattus norvegicus, is an important model of human health and disease, and experimental findings in the rat have relevance to human physiology and disease. The Rat Genome Database (RGD, https://rgd.mcw.edu) is a model organism database that provides access to a wide variety of curated rat data including disease associations, phenotypes, pathways, molecular functions, biological processes, cellular components, and chemical interactions for genes, quantitative trait loci, and strains. We present an overview of the database followed by specific examples that can be used to gain experience in employing RGD to explore the wealth of functional data available for the rat and other species. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Navigating the Rat Genome Database (RGD) home page Basic Protocol 2: Using the RGD search functions Basic Protocol 3: Searching for quantitative trait loci Basic Protocol 4: Using the RGD genome browser (JBrowse) to find phenotypic annotations Basic Protocol 5: Using OntoMate to find gene-disease data Basic Protocol 6: Using MOET to find gene-ontology enrichment Basic Protocol 7: Using OLGA to generate gene lists for analysis Basic Protocol 8: Using the GA tool to analyze ontology annotations for genes Basic Protocol 9: Using the RGD InterViewer tool to find protein interaction data Basic Protocol 10: Using the RGD Variant Visualizer tool to find genetic variant data Basic Protocol 11: Using the RGD Disease Portals to find disease, phenotype, and other information Basic Protocol 12: Using the RGD Phenotypes & Models Portal to find qualitative and quantitative phenotype data and other rat strain-related information Basic Protocol 13: Using the RGD Pathway Portal to find disease and phenotype data via molecular pathways.

The VGNC: expanding standardized vertebrate gene nomenclature

The Vertebrate Gene Nomenclature Committee (VGNC) was established in 2016 as a sister project to the HUGO Gene Nomenclature Committee, to approve gene nomenclature in vertebrate species without an existing dedicated nomenclature committee. The VGNC aims to harmonize gene nomenclature across selected vertebrate species in line with human gene nomenclature, with orthologs assigned the same nomenclature where possible. This article presents an overview of the VGNC project and discussion of key findings resulting from this work to date. VGNC-approved nomenclature is accessible at https://vertebrate.genenames.org and is additionally displayed by the NCBI, Ensembl, and UniProt databases.

The Gene Ontology knowledgebase in 2023

The Gene Ontology (GO) knowledgebase (http://geneontology.org) is a comprehensive resource concerning the functions of genes and gene products (proteins and noncoding RNAs). GO annotations cover genes from organisms across the tree of life as well as viruses, though most gene function knowledge currently derives from experiments carried out in a relatively small number of model organisms. Here, we provide an updated overview of the GO knowledgebase, as well as the efforts of the broad, international consortium of scientists that develops, maintains, and updates the GO knowledgebase. The GO knowledgebase consists of three components: (1) the GO-a computational knowledge structure describing the functional characteristics of genes; (2) GO annotations-evidence-supported statements asserting that a specific gene product has a particular functional characteristic; and (3) GO Causal Activity Models (GO-CAMs)-mechanistic models of molecular "pathways" (GO biological processes) created by linking multiple GO annotations using defined relations. Each of these components is continually expanded, revised, and updated in response to newly published discoveries and receives extensive QA checks, reviews, and user feedback. For each of these components, we provide a description of the current contents, recent developments to keep the knowledgebase up to date with new discoveries, and guidance on how users can best make use of the data that we provide. We conclude with future directions for the project.

2022 updates to the Rat Genome Database: a Findable, Accessible, Interoperable, and Reusable (FAIR) resource

The Rat Genome Database (RGD, https://rgd.mcw.edu) has evolved from simply a resource for rat genetic markers, maps, and genes, by adding multiple genomic data types and extensive disease and phenotype annotations and developing tools to effectively mine, analyze, and visualize the available data, to empower investigators in their hypothesis-driven research. Leveraging its robust and flexible infrastructure, RGD has added data for human and eight other model organisms (mouse, 13-lined ground squirrel, chinchilla, naked mole-rat, dog, pig, African green monkey/vervet, and bonobo) besides rat to enhance its translational aspect. This article presents an overview of the database with the most recent additions to RGD's genome, variant, and quantitative phenotype data. We also briefly introduce Virtual Comparative Map (VCMap), an updated tool that explores synteny between species as an improvement to RGD's suite of tools, followed by a discussion regarding the refinements to the existing PhenoMiner tool that assists researchers in finding and comparing quantitative data across rat strains. Collectively, RGD focuses on providing a continuously improving, consistent, and high-quality data resource for researchers while advancing data reproducibility and fulfilling Findable, Accessible, Interoperable, and Reusable (FAIR) data principles.

The progression of radiation injury in a Wistar rat model of partial body irradiation with ∼5% bone marrow shielding

PURPOSE

To describe the dose response relationship and natural history of radiation injury in the Wistar rat and it's suitability for use in medical countermeasures (MCM) testing.

MATERIALS & METHODS

In two separate studies, male and female rats were exposed to partial body irradiation (PBI) with 5% bone marrow sparing. Animals were X-ray irradiated from 7 to 12 Gy at 7-10 weeks of age. Acute radiation syndrome (ARS) survival at 30 days and delayed effects of acute radiation exposure (DEARE) survival at 182 days were assessed. Radiation effects were determined by clinical observations, body weights, hematology, clinical chemistry, magnetic resonance imaging of lung, whole-body plethysmography, and histopathology.

RESULTS

Rats developed canonical ARS responses of hematopoietic atrophy and gastrointestinal injury resulting in mortality at doses ≥8Gy in males and ≥8.5 Gy in females. DEARE mortality occurred at doses ≥8Gy for both sexes. Findings indicate lung, kidney, and/or liver injury, and persistent hematological dysregulation, revealing multi-organ injury as a DEARE.

CONCLUSION

The Wistar rat PBI model is suitable for testing MCMs against hematopoietic and gastrointestinal ARS. DEARE multi-organ injury occurred in both sexes irradiated with 8-9Gy, also suggesting suitability for polypharmacy studies addressing the combination of ARS and DEARE injury.

Identifying multiscale translational safety biomarkers using a network-based systems approach

Animal testing is the current standard for drug and chemicals safety assessment, but hazards translation to human is uncertain. Human in vitro models can address the species translation but might not replicate in vivo complexity. Herein, we propose a network-based method addressing these translational multiscale problems that derives in vivo liver injury biomarkers applicable to in vitro human early safety screening. We applied weighted correlation network analysis (WGCNA) to a large rat liver transcriptomic dataset to obtain co-regulated gene clusters (modules). We identified modules statistically associated with liver pathologies, including a module enriched for ATF4-regulated genes as associated with the occurrence of hepatocellular single-cell necrosis, and as preserved in human liver in vitro models. Within the module, we identified TRIB3 and MTHFD2 as a novel candidate stress biomarkers, and developed and used BAC-eGFPHepG2 reporters in a compound screening, identifying compounds showing ATF4-dependent stress response and potential early safety signals.

Conservation and divergence of canonical and non-canonical imprinting in murids

BACKGROUND

Genomic imprinting affects gene expression in a parent-of-origin manner and has a profound impact on complex traits including growth and behavior. While the rat is widely used to model human pathophysiology, few imprinted genes have been identified in this murid. To systematically identify imprinted genes and genomic imprints in the rat, we use low input methods for genome-wide analyses of gene expression and DNA methylation to profile embryonic and extraembryonic tissues at allele-specific resolution.

RESULTS

We identify 14 and 26 imprinted genes in these tissues, respectively, with 10 of these genes imprinted in both tissues. Comparative analyses with mouse reveal that orthologous imprinted gene expression and associated canonical DNA methylation imprints are conserved in the embryo proper of the Muridae family. However, only 3 paternally expressed imprinted genes are conserved in the extraembryonic tissue of murids, all of which are associated with non-canonical H3K27me3 imprints. The discovery of 8 novel non-canonical imprinted genes unique to the rat is consistent with more rapid evolution of extraembryonic imprinting. Meta-analysis of novel imprinted genes reveals multiple mechanisms by which species-specific imprinted expression may be established, including H3K27me3 deposition in the oocyte, the appearance of ZFP57 binding motifs, and the insertion of endogenous retroviral promoters.

CONCLUSIONS

In summary, we provide an expanded list of imprinted loci in the rat, reveal the extent of conservation of imprinted gene expression, and identify potential mechanisms responsible for the evolution of species-specific imprinting.

Transcriptomic Mapping of Neurotoxicity Pathways in the Rat Brain in Response to Intraventricular Polymyxin B

Intraventricular or intrathecal administration of polymyxins are increasingly used to treat multidrug-resistant (MDR) Gram-negative bacteria caused infections in the central nervous system (CNS). However, our limited knowledge of the mechanisms underpinning polymyxin-induced neurotoxicity significantly hinders the development of safe and efficacious polymyxin dosing regimens. To this end, we conducted transcriptomic analyses of the rat brain and spinal cord 1 h following intracerebroventricular administration of polymyxin B into rat lateral ventricle at a clinically relevant dose (0.5 mg/kg). Following the treatment, 66 differentially expressed genes (DEGs) were identified in the brain transcriptome while none for the spinal cord (FDR ≤ 0.05, fold-change ≥ 1.5). DEGs were enriched in signaling pathways associated with hormones and neurotransmitters, including dopamine and (nor)epinephrine. Notably, the expression levels of Slc6a3 and Gabra6 were decreased by 20-fold and 4.3-fold, respectively, likely resulting in major perturbations of dopamine and γ-aminobutyric acid signaling in the brain. Mass spectrometry imaging of brain sections revealed a distinct pattern of polymyxin B distribution with the majority accumulating in the injection-side lateral ventricle and subsequently into third and fourth ventricles. Polymyxin B was not detectable in the left lateral ventricle or brain tissue. Electrophysiological measurements on primary cultured rat neurons revealed a large inward current and significant membrane leakage following polymyxin B treatment. Our work demonstrates, for the first time, the key CNS signaling pathways associated with polymyxin neurotoxicity. This mechanistic insight combined with pharmacokinetic/pharmacodynamic dosing strategies will help guide the design of safe and effective intraventricular/intrathecal polymyxin treatment regimens for CNS infections caused by MDR Gram-negative pathogens.

Cardiac copper content and its relationship with heart physiology: Insights based on quantitative genetic and functional analyses using BXD family mice

Background

Copper (Cu) is essential for the functioning of various enzymes involved in important cellular and physiological processes. Although critical for normal cardiac function, excessive accumulation, or deficiency of Cu in the myocardium is detrimental to the heart. Fluctuations in cardiac Cu content have been shown to cause cardiac pathologies and imbalance in systemic Cu metabolism. However, the genetic basis underlying cardiac Cu levels and their effects on heart traits remain to be understood. Representing the largest murine genetic reference population, BXD strains have been widely used to explore genotype-phenotype associations and identify quantitative trait loci (QTL) and candidate genes.

Methods

Cardiac Cu concentration and heart function in BXD strains were measured, followed by QTL mapping. The candidate genes modulating Cu homeostasis in mice hearts were identified using a multi-criteria scoring/filtering approach.

Results

Significant correlations were identified between cardiac Cu concentration and left ventricular (LV) internal diameter and volumes at end-diastole and end-systole, demonstrating that the BXDs with higher cardiac Cu levels have larger LV chamber. Conversely, cardiac Cu levels negatively correlated with LV posterior wall thickness, suggesting that lower Cu concentration in the heart is associated with LV hypertrophy. Genetic mapping identified six QTLs containing a total of 217 genes, which were further narrowed down to 21 genes that showed a significant association with cardiac Cu content in mice. Among those, Prex1 and Irx3 are the strongest candidates involved in cardiac Cu modulation.

Conclusion

Cardiac Cu level is significantly correlated with heart chamber size and hypertrophy phenotypes in BXD mice, while being regulated by multiple genes in several QTLs. Prex1 and Irx3 may be involved in modulating Cu metabolism and its downstream effects and warrant further experimental and functional validations.

Magel2 truncation alters select behavioral and physiological outcomes in a rat model of Schaaf-Yang syndrome

Previous studies in mice have utilized Magel2 gene deletion models to examine the consequences of its absence. We report the generation, molecular validation and phenotypic characterization of a novel rat model with a truncating Magel2 mutation modeling variants associated with Schaaf-Yang syndrome-causing mutations. Within the hypothalamus, a brain region in which human MAGEL2 is paternally expressed, we demonstrated, at the level of transcript and peptide detection, that rat Magel2 exhibits a paternal, parent-of-origin effect. In evaluations of behavioral features across several domains, juvenile Magel2 mutant rats displayed alterations in anxiety-like behavior and sociability measures. Moreover, the analysis of peripheral organ systems detected alterations in body composition, cardiac structure and function, and breathing irregularities in Magel2 mutant rats. Several of these findings are concordant with reported mouse phenotypes, indicating the conservation of MAGEL2 function across rodent species. Our comprehensive analysis revealing impairments across multiple domains demonstrates the tractability of this model system for the study of truncating MAGEL2 mutations.

Multiple solutions at the genomic level in response to selective breeding for high locomotor activity

Replicate lines under uniform selection often evolve in different ways. Previously, analyses using whole-genome sequence data for individual mice (Mus musculus) from 4 replicate High Runner lines and 4 nonselected control lines demonstrated genomic regions that have responded consistently to selection for voluntary wheel-running behavior. Here, we ask whether the High Runner lines have evolved differently from each other, even though they reached selection limits at similar levels. We focus on 1 High Runner line (HR3) that became fixed for a mutation at a gene of major effect (Myh4Minimsc) that, in the homozygous condition, causes a 50% reduction in hindlimb muscle mass and many pleiotropic effects. We excluded HR3 from SNP analyses and identified 19 regions not consistently identified in analyses with all 4 lines. Repeating analyses while dropping each of the other High Runner lines identified 12, 8, and 6 such regions. (Of these 45 regions, 37 were unique.) These results suggest that each High Runner line indeed responded to selection somewhat uniquely, but also that HR3 is the most distinct. We then applied 2 additional analytical approaches when dropping HR3 only (based on haplotypes and nonstatistical tests involving fixation patterns). All 3 approaches identified 7 new regions (as compared with analyses using all 4 High Runner lines) that include genes associated with activity levels, dopamine signaling, hippocampus morphology, heart size, and body size, all of which differ between High Runner and control lines. Our results illustrate how multiple solutions and "private" alleles can obscure general signatures of selection involving "public" alleles.

Positional cloning of rat mutant genes reveals new functions of these genes

The laboratory rat (Rattus norvegicus) is a key model organism for biomedical research. Rats can be subjected to strict genetic and environmental controls. The rat's large body size is suitable for both surgical operations and repeated measurements of physiological parameters. These advantages have led to the development of numerous rat models for genetic diseases. Forward genetics is a proven approach for identifying the causative genes of these disease models but requires genome resources including genetic markers and genome sequences. Over the last few decades, rat genome resources have been developed and deposited in bioresource centers, which have enabled us to perform positional cloning in rats. To date, more than 100 disease-related genes have been identified by positional cloning. Since some disease models are more accessible in rats than mice, the identification of causative genes in these models has sometimes led to the discovery of novel functions of genes. As before, various mutant rats are also expected to be discovered and developed as disease models in the future. Thus, the forward genetics continues to be an important approach to find genes involved in disease phenotypes in rats. In this review, I provide an overview the development of rat genome resources and describe examples of positional cloning in rats in which novel gene functions have been identified.

Hunting for Genes Underlying Emotionality in the Laboratory Rat: Maps, Tools and Traps

Scientists have systematically investigated the hereditary bases of behaviors since the 19th century, moved by either evolutionary questions or clinically-motivated purposes. The pioneer studies on the genetic selection of laboratory animals had already indicated, one hundred years ago, the immense complexity of analyzing behaviors that were influenced by a large number of small-effect genes and an incalculable amount of environmental factors. Merging Mendelian, quantitative and molecular approaches in the 1990s made it possible to map specific rodent behaviors to known chromosome regions. From that point on, Quantitative Trait Locus (QTL) analyses coupled with behavioral and molecular techniques, which involved in vivo isolation of relevant blocks of genes, opened new avenues for gene mapping and characterization. This review examines the QTL strategy applied to the behavioral study of emotionality, with a focus on the laboratory rat. We discuss the challenges, advances and limitations of the search for Quantitative Trait Genes (QTG) playing a role in regulating emotionality. For the past 25 years, we have marched the long journey from emotionality-related behaviors to genes. In this context, our experiences are used to illustrate why and how one should move forward in the molecular understanding of complex psychiatric illnesses. The promise of exploring genetic links between immunological and emotional responses are also discussed. New strategies based on humans, rodents and other animals (such as zebrafish) are also acknowledged, as they are likely to allow substantial progress to be made in the near future.

Novel peanut-specific human IgE monoclonal antibodies enable screens for inhibitors of the effector phase in food allergy

Background

10% of US residents have food allergies, including 2% with peanut allergy. Mast cell mediators released during the allergy effector phase drive allergic reactions. Therefore, targeting sensitized mast cells may prevent food allergy symptoms.

Objective

We used novel, human, allergen-specific, IgE monoclonal antibodies (mAbs) created using human hybridoma techniques to design an in vitro system to evaluate potential therapeutics targeting sensitized effector cells.

Methods

Two human IgE mAbs specific for peanut, generated through human hybridoma techniques, were used to sensitize rat basophilic leukemia (RBL) SX-38 cells expressing the human IgE receptor (FcϵRI). Beta-hexosaminidase release (a marker of degranulation), cytokine production, and phosphorylation of signal transduction proteins downstream of FcϵRI were measured after stimulation with peanut. Degranulation was also measured after engaging inhibitory receptors CD300a and Siglec-8.

Results

Peanut-specific human IgE mAbs bound FcϵRI, triggering degranulation after stimulation with peanut in RBL SX-38 cells. Sensitized RBL SX-38 cells stimulated with peanut increased levels of phosphorylated SYK and ERK, signal transduction proteins downstream of FcϵRI. Engaging inhibitory cell surface receptors CD300a or Siglec-8 blunted peanut-specific activation.

Conclusion

Allergen-specific human IgE mAbs, expressed from human hybridomas and specific for a clinically relevant food allergen, passively sensitize allergy effector cells central to the in vitro models of the effector phase of food allergy. Peanut reproducibly activates and induces degranulation of RBL SX-38 cells sensitized with peanut-specific human IgE mAbs. This system provides a unique screening tool to assess the efficacy of therapeutics that target allergy effector cells and inhibit food allergen-induced effector cell activation.

Trace Element Interactions, Inflammatory Signaling, and Male Sex Implicated in Reduced Growth Following Excess Oral Iron Supplementation in Pre-Weanling Rats

Iron supplements are frequently provided to infants in high-income countries despite low incidence of iron deficiency. There is growing concern regarding adverse health and development outcomes of excess iron provision in early life. Excess iron may directly damage developing organs through the formation of reactive oxygen species, alter systemic inflammatory signaling, and/or dysregulate trace mineral metabolism. To better characterize the in vivo effects of excess iron on development, we utilized a pre-weanling rat pup model. Lewis rat litters were culled to eight pups (four males and four females) and randomly assigned to daily supplementation groups receiving either vehicle control (CON; 10% w/v sucrose solution) or ferrous sulfate (FS) iron at one of the following doses: 10, 30, or 90 mg iron/kg body weight—FS-10, FS-30, and FS-90, respectively—from postnatal day (PD) 2 through 9. FS-90 litters, but not FS-30 or FS-10, failed to thrive compared to CON litters and had smaller brains on PD 10. Among the groups, FS-90 liver iron levels were highest, as were white blood cell counts. Compared to CON, circulating MCP-1 and liver zinc were increased in FS-90 pups, whereas liver copper was decreased. Growth defects due to excess FS provision in pre-weanling rats may be related to liver injury, inflammation, and altered trace mineral metabolism.

Targeting OXPHOS de novo purine synthesis as the nexus of FLT3 inhibitor-mediated synergistic antileukemic actions

Using a genome-wide CRISPR screen, we identified CDK9, DHODH, and PRMT5 as synthetic lethal partners with gilteritinib treatment in fms-like tyrosine kinase 3 (FLT3)-internal tandem duplication (ITD) acute myeloid leukemia (AML) and genetically and pharmacologically validated their roles in gilteritinib sensitivity. The presence of FLT3-ITD is associated with an increase in anaerobic glycolysis, rendering leukemia cells highly sensitive to inhibition of glycolysis. Supportive of this, our data show the enrichment of single guide RNAs targeting 28 glycolysis-related genes upon gilteritinib treatment, suggesting that switching from glycolysis to oxidative phosphorylation (OXPHOS) may represent a metabolic adaption of AML in gilteritinib resistance. CDK9i/FLT3i, DHODHi/FLT3i, and PRMT5i/FLT3i pairs mechanistically converge on OXPHOS and purine biosynthesis blockade, implying that targeting the metabolic functions of these three genes and/or proteins may represent attractive strategies to sensitize AML to gilteritinib treatment. Our findings provide the basis for maximizing therapeutic impact of FLT3-ITD inhibitors and a rationale for a clinical trial of these novel combinations.

Identification of Hypothalamic Long Noncoding RNAs Associated with Hypertension and the Behavior/Neurological Phenotype of Hypertensive ISIAH Rats

Long noncoding RNAs (lncRNAs) play an important role in the control of many physiological and pathophysiological processes, including the development of hypertension and other cardiovascular diseases. Nonetheless, the understanding of the regulatory function of many lncRNAs is still incomplete. This work is a continuation of our earlier study on the sequencing of hypothalamic transcriptomes of hypertensive ISIAH rats and control normotensive WAG rats. It aims to identify lncRNAs that may be involved in the formation of the hypertensive state and the associated behavioral features of ISIAH rats. Interstrain differences in the expression of seven lncRNAs were validated by quantitative PCR. Differential hypothalamic expression of lncRNAs LOC100910237 and RGD1562890 between hypertensive and normotensive rats was shown for the first time. Expression of four lncRNAs (Snhg4, LOC100910237, RGD1562890, and Tnxa-ps1) correlated with transcription levels of many hypothalamic genes differentially expressed between ISIAH and WAG rats (DEGs), including genes associated with the behavior/neurological phenotype and hypertension. After functional annotation of these DEGs, it was concluded that lncRNAs Snhg4, LOC100910237, RGD1562890, and Tnxa-ps1 may be involved in the hypothalamic processes related to immune-system functioning and in the response to various exogenous and endogenous factors, including hormonal stimuli. Based on the functional enrichment analysis of the networks, an association of lncRNAs LOC100910237 and Tnxa-ps1 with retinol metabolism and an association of lncRNAs RGD1562890 and Tnxa-ps1 with type 1 diabetes mellitus are proposed for the first time. Based on a discussion, it is hypothesized that previously functionally uncharacterized lncRNA LOC100910237 is implicated in the regulation of hypothalamic processes associated with dopaminergic synaptic signaling, which may contribute to the formation of the behavioral/neurological phenotype and hypertensive state of ISIAH rats.

PhenoGeneRanker: Gene and Phenotype Prioritization Using Multiplex Heterogeneous Networks

Uncovering genotype-phenotype relationships is a fundamental challenge in genomics. Gene prioritization is an important step for this endeavor to make a short manageable list from a list of thousands of genes coming from high-throughput studies. Network propagation methods are promising and state of the art methods for gene prioritization based on the premise that functionally related genes tend to be close to each other in the biological networks. Recently, we introduced PhenoGeneRanker, a network-propagation algorithm for multiplex heterogeneous networks. PhenoGeneRanker allows multi-layer gene and phenotype networks. It also calculates empirical p values of gene and phenotype ranks using random stratified sampling of seeds of genes and phenotypes based on their connectivity degree in the network. In this study, we introduce the PhenoGeneRanker Bioconductor package and its application to multi-omics rat genome datasets to rank hypertension disease-related genes and strains. We showed that PhenoGeneRanker performed better to rank hypertension disease-related genes using multiplex gene networks than aggregated gene networks. We also showed that PhenoGeneRanker performed better to rank hypertension disease-related strains using multiplex phenotype network than single or aggregated phenotype networks. We performed a rigorous hyperparameter analysis and, finally showed that Gene Ontology (GO) enrichment of statistically significant top-ranked genes resulted in hypertension disease-related GO terms.

Nutrigenetic Interaction of Spontaneously Hypertensive Rat Chromosome 20 Segment and High-Sucrose Diet Sensitizes to Metabolic Syndrome

Several corresponding regions of human and mammalian genomes have been shown to affect sensitivity to the manifestation of metabolic syndrome via nutrigenetic interactions. In this study, we assessed the effect of sucrose administration in a newly established congenic strain BN.SHR20, in which a limited segment of rat chromosome 20 from a metabolic syndrome model, spontaneously hypertensive rat (SHR), was introgressed into Brown Norway (BN) genomic background. We mapped the extent of the differential segment and compared the genomic sequences of BN vs. SHR within the segment in silico. The differential segment of SHR origin in BN.SHR20 spans about 9 Mb of the telomeric portion of the short arm of chromosome 20. We identified non-synonymous mutations e.g., in ApoM, Notch4, Slc39a7, Smim29 genes and other variations in or near genes associated with metabolic syndrome in human genome-wide association studies. Male rats of BN and BN.SHR20 strains were fed a standard diet for 18 weeks (control groups) or 16 weeks of standard diet followed by 14 days of high-sucrose diet (HSD). We assessed the morphometric and metabolic profiles of all groups. Adiposity significantly increased only in BN.SHR20 after HSD. Fasting glycemia and the glucose levels during the oral glucose tolerance test were higher in BN.SHR20 than in BN groups, while insulin levels were comparable. The fasting levels of triacylglycerols were the highest in sucrose-fed BN.SHR20, both compared to the sucrose-fed BN and the control BN.SHR20. The non-esterified fatty acids and total cholesterol concentrations were higher in BN.SHR20 compared to their respective BN groups, and the HSD elicited an increase in non-esterified fatty acids only in BN.SHR20. In a new genetically defined model, we have isolated a limited genomic region involved in nutrigenetic sensitization to sucrose-induced metabolic disturbances.

GlioMarker: An integrated database for knowledge exploration of diagnostic biomarkers in gliomas

Gliomas are the most frequent malignant and aggressive tumors in the central nervous system. Early and effective diagnosis of glioma using diagnostic biomarkers can prolong patients' lives and aid in the development of new personalized treatments. Therefore, a thorough and comprehensive understanding of the diagnostic biomarkers in gliomas is of great significance. To this end, we developed the integrated and web-based database GlioMarker (http://gliomarker.prophetdb.org/), the first comprehensive database for knowledge exploration of glioma diagnostic biomarkers. In GlioMarker, accurate information on 406 glioma diagnostic biomarkers from 1559 publications was manually extracted, including biomarker descriptions, clinical information, associated literature, experimental records, associated diseases, statistical indicators, etc. Importantly, we integrated many external resources to provide clinicians and researchers with the capability to further explore knowledge on these diagnostic biomarkers based on three aspects. (1) Obtain more ontology annotations of the biomarker. (2) Identify the relationship between any two or more components of diseases, drugs, genes, and variants to explore the knowledge related to precision medicine. (3) Explore the clinical application value of a specific diagnostic biomarker through online analysis of genomic and expression data from glioma cohort studies. GlioMarker provides a powerful, practical, and user-friendly web-based tool that may serve as a specialized platform for clinicians and researchers by providing rapid and comprehensive knowledge of glioma diagnostic biomarkers to subsequently facilitates high-quality research and applications.

Genome-wide linkage analysis combined with genome sequencing in large families with intracranial aneurysms

Rupture of an intracranial aneurysm (IA) leads to aneurysmal subarachnoid haemorrhage (ASAH), a severe type of stroke. Some rare variants that cause IA in families have been identified, but still, the majority of genetic causes, as well as the biological mechanisms of IA development and rupture, remain unknown. We aimed to identify rare, damaging variants for IA in three large Dutch families with multiple affected members with IA (N = 9, 11, and 6). By combining linkage analysis and genome sequencing (GS), we identified six rare and damaging variants for which all cases within one of the families were heterozygous. These variants were p.Tyr87Cys in SYCP1, p.Phe1077Leu in FMNL2, p.Thr754Lys in TBC1D2, p.Arg321His in ZNF782, p.Arg979Trp in CCDC180, and p.Val125Met in NCBP1. None of the variants showed association with IA status in a large cohort of 937 patients from the general IA patient population and 1046 controls. Gene expression in IA and cerebral artery tissue further prioritized FMNL2 and TBC1D2 as potential important players in IA pathophysiology. Further studies are needed to characterize the functional consequences of the identified variants and their role in the biological mechanisms of IA.

Body Composition and Metabolic Changes in a Lyon Hypertensive Congenic Rat and Identification of Ercc6l2 as a Positional Candidate Gene

Central obesity is genetically complex, and its exponential increase in the last decades have made it a critical public health issue. The Lyon Hypertensive (LH) rat is a well-characterized hypertensive model that also exhibits spontaneous and profound differences in body weight and adiposity, relative to its metabolically healthy control, the Lyon Normotensive (LN) rat. The mechanisms underlying the body weight differences between these strains are not well-understood, thus a congenic model (LH17LNa) was developed where a portion of the proximal arm of LN chromosome 17 is introgressed on the LH genomic background to assess the contribution of LN alleles on obesity features. Male and female LH17LNa rats were studied, but male congenics did not significantly differ from LH in this study. Female LH17LNa rats exhibited decreases in total body growth, as well as major alterations to their body composition and adiposity. The LH17LNa female rats also showed decreases in metabolic rate, and a reduction in food intake. The increased adiposity in the female LH17LNa rats was specific to abdominal white adipose tissue, and this phenomenon was further explained by significant hypertrophy in those adipocytes, with no evidence of adipocyte hyperplasia. Sequencing of the parental strains identified a novel frameshift mutation in the candidate gene Ercc6l2, which is involved in transcription-coupled DNA repair, and is implicated in premature aging. The discovery of the significance of Ercc6l2 in the context of female-specific adipocyte biology could represent a novel role of DNA repair failure syndromes in obesity pathogenesis.

Protective effect of pterostilbene in a streptozotocin-induced mouse model of Alzheimer's disease by targeting monoamine oxidase B

Alzheimer's disease (AD) is a neurodegenerative disease in elderly population. Pterostilbene (PTS) is a resveratrol analogue with neuroprotective activity. However, the biological mechanisms of PTS in AD progression are largely uncertain. An animal model of AD was established using streptozotocin (STZ)-treated C57BL/6J mice. Monoamine oxidase B (MAOB) expression was analyzed by bioinformatics analysis and detected by western blotting assay. The memory impairment was investigated by Morris water maze test. The levels of Tau hyperphosphorylation and death-related proteins were detected by western blotting analysis. The levels of amyloid β (Aβ)_1-42 accumulation, oxidative stress-related markers (ROS, MDA, SOD and GSH), and inflammation-relative markers (TNF-α, IL-1β, IL-6 and p-NF-κB) were measured by ELISA. MAOB expression was increased in hippocampus of AD mice, and it was decreased by PTS. PTS attenuated STZ-induced body weight loss and memory impairment by regulating MAOB. PTS mitigated Aβ_1-42 accumulation and Tau hyperphosphorylation by regulating MAOB in STZ-treated mice. PTS attenuated neuronal death by decreasing cleaved caspase-3 and Bax levels and increasing Bcl2 expression in hippocampus by regulating MAOB in STZ-treated mice. PTS weakened STZ-induced oxidative stress in hippocampus by decreasing ROS and MDA levels and increasing SOD and GSH levels by regulating MAOB. PTS protected against STZ-induced neuroinflammation in hippocampus by inhibiting TNF-α, IL-1β, IL-6 and p-NF-κB levels through regulating MAOB. In conclusion, PTS alleviates STZ-induced memory impairment, Aβ_1-42 accumulation, Tau hyperphosphorylation, neuronal death, oxidative stress and inflammation by decreasing MAOB in AD mice, proving anti-AD potential of PTS.

Preventive aerobic training preserves sympathovagal function and improves DNA repair capacity of peripheral blood mononuclear cells in rats with cardiomyopathy

To evaluate the effect of preventive aerobic exercise training on sympathovagal function, cardiac function, and DNA repair capacity in a preclinical model of doxorubicin (DOX)-induced cardiomyopathy. Forty male Wistar-Kyoto rats were allocated into four groups (n = 10/group): D (DOX-treated) and C (controls) remained sedentary, and DT (DOX-trained) and CT (control-trained) performed aerobic training 4 days/week, during 4 weeks before exposure to DOX (4 mg/kg/week during 4 weeks) or saline solution. We evaluated cardiac function (echocardiography), hemodynamic and sympathovagal modulation (artery-femoral cannulation), cardiac troponin T levels, and DNA repair capacity (comet assay). Exercise training preserved ejection fraction (D: − 14.44% vs. DT: − 1.05%, p < 0.001), fractional shortening (D: − 8.96% vs. DT: − 0.27%, p = 0.025) and troponin T levels (D: 6.4 ± 3.6 vs. DT: 2.8 ± 1.7 ng/mL, p = 0.010). DOX increased heart rate variability (C: 27.7 ± 7.9 vs. D: 7.5 ± 2.2 ms², p < 0.001) and induced sympathovagal dysfunction (LF/HF, C: 0.37 ± 0.15 vs. D: 0.15 ± 0.15, p = 0.036) through exacerbation of sympathetic function (LF, C: 0.22 ± 0.01 vs. D: 0.48 ± 0.24 Hz, p = 0.019). Peripheral mononuclear blood cells of DT animals presented lower residual DNA damage (D: 43.4 ± 8.4% vs. DT: 26 ± 3.4%, p = 0.003 after 1 h). Cardioprotective effects of preventive aerobic exercise training are mediated by preservation of sympathovagal function and improvement of DNA repair capacity of peripheral blood mononuclear cells.

TNF-α Induced Myotube Atrophy in C2C12 Cell Line Uncovers Putative Inflammatory-Related lncRNAs Mediating Muscle Wasting

Background: Muscle atrophy is a complex catabolic condition developing under different inflammatory-related systemic diseases resulting in wasting of muscle tissue. While the knowledge of the molecular background of muscle atrophy has developed in recent years, how the atrophic conditions affect the long non-coding RNA (lncRNAs) machinery and the exact participation of the latter in the mediation of muscle loss are still unknown. The purpose of the study was to assess how inflammatory condition developing under the tumor necrosis factor alpha (TNF-α) treatment affects the lncRNAs’ expression in a mouse skeletal muscle cell line. Materials and method: A C2C12 mouse myoblast cell line was treated with TNF-α to develop atrophy, and inflammatory-related lncRNAs mediating muscle loss were identified. Bioinformatics was used to validate and analyze the discovered lncRNAs. The differences in their expression under different TNF-α concentrations and treatment times were investigated. Results: Five lncRNAs were identified in a discovery set as atrophy related and then validated. Three lncRNAs, Gm4117, Ccdc41os1, and 5830418P13Rik, were selected as being significant for inflammatory-related myotube atrophy. Dynamics changes in the expression of lncRNAs depended on both TNF-α concentration and treatment time. Bioinformatics analysis revealed the mRNA and miRNA target for selected lncRNAs and their putative involvement in the molecular processes related to muscle atrophy. Conclusions: The inflammatory condition developing in the myotube under the TNF-α treatment affects the alteration of lncRNAs’ expression pattern. Experimental and bioinformatics testing suggested the prospective role of lncRNAs in the mediation of muscle loss under an inflammatory state.

Analysis of Genomic Alternative Splicing Patterns in Rat under Heat Stress Based on RNA-Seq Data

Heat stress is one of the most severe challenges faced in livestock production in summer. Alternative splicing as an important post-transcriptional regulation is rarely studied in heat-stressed animals. Here, we performed and analyzed RNA-sequencing assays on the liver of Sprague-Dawley rats in control (22 °C, n = 5) and heat stress (4 °C for 120 min, H120; n = 5) groups, resulting in the identification of 636 differentially expressed genes. Identification analysis of the alternative splicing events revealed that heat stress-induced alternative splicing events increased by 20.18%. Compared with other types of alternative splicing events, the alternative start increased the most (43.40%) after heat stress. Twenty-eight genes were differentially alternatively spliced (DAS) between the control and H120 groups, among which Acly, Hnrnpd and mir3064 were also differentially expressed. For DAS genes, Srebf1, Shc1, Srsf5 and Ensa were associated with insulin, while Cast, Srebf1, Tmem33, Tor1aip2, Slc39a7 and Sqstm1 were enriched in the composition of the endoplasmic reticulum. In summary, our study conducts a comprehensive profile of alternative splicing in heat-stressed rats, indicating that alternative splicing is one of the molecular mechanisms of heat stress response in mammals and providing reference data for research on heat tolerance in mammalian livestock.

Exercise-Induced Extracellular Vesicles Delay the Progression of Prostate Cancer

Increasing evidence suggests that regular physical exercise not only reduces the risk of cancer but also improves functional capacity, treatment efficacy and disease outcome in cancer patients. At least partially, these effects are mediated by the secretome of the tissues responding to exercise. The secreted molecules can be released in a carrier-free form or enclosed into extracellular vesicles (EVs). Several recent studies have shown that EVs are actively released into circulation during physical exercise. Here, we for the first time investigated the effects of exercise-induced EVs on the progression of cancer in an F344 rat model of metastatic prostate cancer. Although we did not observe a consistent increase in the circulating EV levels, RNA sequencing analysis demonstrated substantial changes in the RNA content of EVs collected before and immediately after forced wheel running exercise as well as differences between EVs from runners at resting state and sedentary rats. The major RNA biotype in EVs was mRNA, followed by miRNA and rRNA. Molecular functions of differentially expressed RNAs reflected various physiological processes including protein folding, metabolism and regulation of immune responses triggered by the exercise in the parental cells. Intravenous administration of exercise-induced EVs into F344 rats with orthotopically injected syngeneic prostate cancer cells PLS10, demonstrated reduction of the primary tumor volume by 35% and possibly—attenuation of lung metastases. Hence, our data provide the first evidence that exercise-induced EVs may modulate tumor physiology and delay the progression of cancer.

The Rat Genome Database (RGD) facilitates genomic and phenotypic data integration across multiple species for biomedical research

Model organism research is essential for discovering the mechanisms of human diseases by defining biologically meaningful gene to disease relationships. The Rat Genome Database (RGD, ( https://rgd.mcw.edu )) is a cross-species knowledgebase and the premier online resource for rat genetic and physiologic data. This rich resource is enhanced by the inclusion and integration of comparative data for human and mouse, as well as other human disease models including chinchilla, dog, bonobo, pig, 13-lined ground squirrel, green monkey, and naked mole-rat. Functional information has been added to records via the assignment of annotations based on sequence similarity to human, rat, and mouse genes. RGD has also imported well-supported cross-species data from external resources. To enable use of these data, RGD has developed a robust infrastructure of standardized ontologies, data formats, and disease- and species-centric portals, complemented with a suite of innovative tools for discovery and analysis. Using examples of single-gene and polygenic human diseases, we illustrate how data from multiple species can help to identify or confirm a gene as involved in a disease and to identify model organisms that can be studied to understand the pathophysiology of a gene or pathway. The ultimate aim of this report is to demonstrate the utility of RGD not only as the core resource for the rat research community but also as a source of bioinformatic tools to support a wider audience, empowering the search for appropriate models for human afflictions.

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

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

The human hepatocyte TXG-MAPr: gene co-expression network modules to support mechanism-based risk assessment

Mechanism-based risk assessment is urged to advance and fully permeate into current safety assessment practices, possibly at early phases of drug safety testing. Toxicogenomics is a promising source of mechanisms-revealing data, but interpretative analysis tools specific for the testing systems (e.g. hepatocytes) are lacking. In this study, we present the TXG-MAPr webtool (available at https://txg-mapr.eu/WGCNA_PHH/TGGATEs_PHH/ ), an R-Shiny-based implementation of weighted gene co-expression network analysis (WGCNA) obtained from the Primary Human Hepatocytes (PHH) TG-GATEs dataset. The 398 gene co-expression networks (modules) were annotated with functional information (pathway enrichment, transcription factor) to reveal their mechanistic interpretation. Several well-known stress response pathways were captured in the modules, were perturbed by specific stressors and showed preservation in rat systems (rat primary hepatocytes and rat in vivo liver), with the exception of DNA damage and oxidative stress responses. A subset of 87 well-annotated and preserved modules was used to evaluate mechanisms of toxicity of endoplasmic reticulum (ER) stress and oxidative stress inducers, including cyclosporine A, tunicamycin and acetaminophen. In addition, module responses can be calculated from external datasets obtained with different hepatocyte cells and platforms, including targeted RNA-seq data, therefore, imputing biological responses from a limited gene set. As another application, donors' sensitivity towards tunicamycin was investigated with the TXG-MAPr, identifying higher basal level of intrinsic immune response in donors with pre-existing liver pathology. In conclusion, we demonstrated that gene co-expression analysis coupled to an interactive visualization environment, the TXG-MAPr, is a promising approach to achieve mechanistic relevant, cross-species and cross-platform evaluation of toxicogenomic data.

Integrative Map of HIF1A Regulatory Elements and Variations

Hypoxia-inducible factor (HIF) family of transcription factors (HIF1A, EPAS1, and HIF3A) are regulators of the cellular response to hypoxia. They have been shown to be involved in development of various diseases such as cancer, diabetes, and erythrocytosis. A complete map of connections between HIF family of genes with various omics types has not yet been developed. The main aim of the present analysis was to construct the integrative map of genomic elements associated with HIF1A gene and prioritize potentially deleterious variants. Various genomic databases and bioinformatics tools were used, including Ensembl, MirTarBase, STRING, Cytoscape, MethPrimer, CADD, SIFT, and UALCAN. Integrative HIF1A gene map was visualized and includes transcriptional and post-transcriptional regulators, downstream targets, and genetic variants. One CpG island overlaps transcription start site of the HIF1A gene. Out of over 450 missense variants, four have predicted deleterious effect on protein function by at least five bioinformatics tools. Currently there are 85 miRNAs reported to target HIF1A. HIF1A downstream targets include protein-coding genes, long noncoding RNAs, and microRNAs (hypoxamiRs). The study presents the first integration of heterogeneous molecular interactions associated with HIF1A gene enabling a holistic view of the gene and lays the groundwork for supplementing the data in the future.

Domestication Explains Two-Thirds of Differential-Gene-Expression Variance between Domestic and Wild Animals; The Remaining One-Third Reflects Intraspecific and Interspecific Variation

Belyaev's concept of destabilizing selection during domestication was a major achievement in the XX century. Its practical value has been realized in commercial colors of the domesticated fox that never occur in the wild and has been confirmed in a wide variety of pet breeds. Many human disease models involving animals allow to test drugs before human testing. Perhaps this is why investigators doing transcriptomic profiling of domestic versus wild animals have searched for breed-specific patterns. Here we sequenced hypothalamic transcriptomes of tame and aggressive rats, identified their differentially expressed genes (DEGs), and, for the first time, applied principal component analysis to compare them with all the known DEGs of domestic versus wild animals that we could find. Two principal components, PC1 and PC2, respectively explained 67% and 33% of differential-gene-expression variance (hereinafter: log2 value) between domestic and wild animals. PC1 corresponded to multiple orthologous DEGs supported by homologs; these DEGs kept the log2 value sign from species to species and from tissue to tissue (i.e., a common domestication pattern). PC2 represented stand-alone homologous DEG pairs reversing the log2 value sign from one species to another and from tissue to tissue (i.e., representing intraspecific and interspecific variation).

Rodent Estrous Cycle Monitoring Utilizing Vaginal Lavage: No Such Thing As a Normal Cycle

The current methodology establishes a reproducible, standardized, and cost-effective approach to monitoring the estrous cycle of female Sprague Dawley (SD) adolescent rats. This study demonstrates the complexity of hormonal cycles and the broad spectrum of understanding required to construct a reliable and valid monitoring technique. Through an in-depth examination of principal experimental design and procedural elements, this description of the cycle and its fundamental principles provides a framework for further understanding and deconstructs misconceptions for future replication. Along with an outline of the sample collection process employing vaginal lavage, the procedure describes the mechanism of data categorization into the four-stage model of proestrus, estrus, metestrus, and diestrus. These stages are characterized by a new proposed approach, utilizing the 4 categorizing determinants of vaginal fluid condition, cell type(s) present, cell arrangement, and cell quantity at the time of collection. Variations of each stage, favorable and unfavorable samples, the distinction between cyclicity and acyclicity, and graphic depictions of the collected categorizing components are presented alongside effective interpretive and organizational practices of the data. Overall, these tools allow for the publication of quantifiable data ranges for the first time, leading to the standardization of categorization factors upon replication.

Development of mammary cancer in γ-irradiated F1 hybrids of susceptible Sprague-Dawley and resistant Copenhagen rats, with copy-number losses that pinpoint potential tumor suppressors

Copenhagen rats are highly resistant to mammary carcinogenesis, even after treatment with chemical carcinogens and hormones; most studies indicate that this is a dominant genetic trait. To test whether this trait is also dominant after radiation exposure, we characterized the susceptibility of irradiated Copenhagen rats to mammary carcinogenesis, as well as its inheritance, and identified tumor-suppressor genes that, when inactivated or mutated, may contribute to carcinogenesis. To this end, mammary cancer-susceptible Sprague-Dawley rats, resistant Copenhagen rats, and their F1 hybrids were irradiated with 4 Gy of γ-rays, and tumor development was monitored. Copy-number variations and allelic imbalances of genomic DNA were studied using microarrays and PCR analysis of polymorphic markers. Gene expression was assessed by quantitative PCR in normal tissues and induced mammary cancers of F1 rats. Irradiated Copenhagen rats exhibited a very low incidence of mammary cancer. Unexpectedly, this resistance trait did not show dominant inheritance in F1 rats; rather, they exhibited intermediate susceptibility levels (i.e., between those of their parent strains). The susceptibility of irradiated F1 rats to the development of benign mammary tumors (i.e., fibroadenoma and adenoma) was also intermediate. Copy-number losses were frequently observed in chromosome regions 1q52-54 (24%), 2q12-15 (33%), and 3q31-42 (24%), as were focal (38%) and whole (29%) losses of chromosome 5. Some of these chromosomal regions exhibited allelic imbalances. Many cancer-related genes within these regions were downregulated in mammary tumors as compared with normal mammary tissue. Some of the chromosomal losses identified have not been reported previously in chemically induced models, implying a novel mechanism inherent to the irradiated model. Based on these findings, Sprague-Dawley × Copenhagen F1 rats offer a useful model for exploring genes responsible for radiation-induced mammary cancer, which apparently are mainly located in specific regions of chromosomes 1, 2, 3 and 5.

Coding variants in mouse and rat model organisms: mousepost and ratpost

Mice and rats are the most commonly used vertebrate model organisms in biomedical research. The availability of a reference genome in both animals combined with the deep sequencing of several doze of popular inbred lines also provides rich sequence variation data in these species. In some cases, such sequence variants can be linked directly to a distinctive phenotype. In previous work, we created the mouse and rat online searchable databases ("Mousepost" and "Ratpost") where small variant information for protein coding transcripts in mouse and rat inbred strains can be easily retrieved at the amino acid level. These tools are directly useful in forward genetics strategies or as a repository of existing sequence variations. Here, we perform a comparison between the "Mousepost" and "Ratpost" databases and we couple these two tools to a database of human sequence variants ClinVar. We investigated the level of redundancy and complementarity of known variants in protein coding transcripts and found that the large majority of variants is species-specific. However, a small set of positions is conserved in an inbred line between both species. We conclude that both databases are highly complementary, but this may change with further sequencing efforts in both species.

Prefrontal Cortex Metabolome Is Modified by Opioids, Anesthesia, and Sleep

Obtundation of wakefulness caused by opioids and loss of wakefulness caused by anesthetics and sleep significantly alter concentrations of molecules comprising the prefrontal cortex (PFC) metabolome. Quantifying state-selective changes in the PFC metabolome is essential for advancing functional metabolomics. Diverse functions of the PFC suggest the PFC metabolome as a potential therapeutic entry point for countermeasures to state-selective autonomic dysfunction.

The Cocaine and Oxycodone Biobanks, Two Repositories from Genetically Diverse and Behaviorally Characterized Rats for the Study of Addiction

The rat oxycodone and cocaine biobanks contain samples that vary by genotypes (by using genetically diverse genotyped HS rats), phenotypes (by measuring addiction-like behaviors in an advanced SA model), timepoints (samples are collected longitudinally before, during, and after SA, and terminally at three different timepoints in the addiction cycle: intoxication, withdrawal, and abstinence or without exposure to drugs through age-matched naive rats), samples collected (organs, cells, biofluids, feces), preservation (paraformaldehyde-fixed, snap-frozen, or cryopreserved) and application (proteomics, transcriptomics, microbiomics, metabolomics, epigenetics, anatomy, circuitry analysis, biomarker discovery, etc.Substance use disorders (SUDs) are pervasive in our society and have substantial personal and socioeconomical costs. A critical hurdle in identifying biomarkers and novel targets for medication development is the lack of resources for obtaining biological samples with a detailed behavioral characterization of SUD. Moreover, it is nearly impossible to find longitudinal samples. As part of two ongoing large-scale behavioral genetic studies in heterogeneous stock (HS) rats, we have created two preclinical biobanks using well-validated long access (LgA) models of intravenous cocaine and oxycodone self-administration (SA) and comprehensive characterization of addiction-related behaviors. The genetic diversity in HS rats mimics diversity in the human population and includes individuals that are vulnerable or resilient to compulsive-like responding for cocaine or oxycodone. Longitudinal samples are collected throughout the experiment, before exposure to the drug, during intoxication, acute withdrawal, and protracted abstinence, and include naive, age-matched controls. Samples include, but are not limited to, blood plasma, feces and urine, whole brains, brain slices and punches, kidney, liver, spleen, ovary, testis, and adrenal glands. Three preservation methods (fixed in formaldehyde, snap-frozen, or cryopreserved) are used to facilitate diverse downstream applications such as proteomics, metabolomics, transcriptomics, epigenomics, microbiomics, neuroanatomy, biomarker discovery, and other cellular and molecular approaches. To date, >20,000 samples have been collected from over 1000 unique animals and made available free of charge to non-profit institutions through https://www.cocainebiobank.org/ and https://www.oxycodonebiobank.org/.

Dynamic miRNA changes during the process of epileptogenesis in an infantile and adult-onset model

Temporal lobe epilepsy (TLE) is the most common epilepsy type. TLE onset in infancy aggravates features like severity, drug responsiveness, or development of comorbidities. These aggravations may arise from altered micro RNA (miRNA) expression specific to the early onset of the disease. Although the miRNA involvement in TLE is widely studied, the relationship between the onset-age and miRNA expression has not been addressed. Here, we investigated the miRNA profile of infantile and adult-onset TLE in rats combining sequencing and PCR. Since miRNA expression changes with the disease progression, we scrutinized miRNA dynamics across three stages: acute, latent, and chronic. We report that infantile-onset TLE leads to changes in the expression of fewer miRNAs across these stages. Interestingly, the miRNA profile in the acute stage of infantile-onset TLE overlaps in dysregulation of miR-132-5p, -205, and -211-3p with the chronic stage of the disease starting in adulthood. The analysis of putative targets linked the majority of dysregulated miRNAs with pathways involved in epilepsy. Our profiling uncovered miRNA expression characteristic for infantile and adulthood-onset epileptogenesis, suggesting the distinct biology underlying TLE in the onset age-dependent matter. Our results indicate the necessity of addressing the onset age as an important parameter in future epilepsy research.

Advances in Genome Editing and Application to the Generation of Genetically Modified Rat Models

The rat has been extensively used as a small animal model. Many genetically engineered rat models have emerged in the last two decades, and the advent of gene-specific nucleases has accelerated their generation in recent years. This review covers the techniques and advances used to generate genetically engineered rat lines and their application to the development of rat models more broadly, such as conditional knockouts and reporter gene strains. In addition, genome-editing techniques that remain to be explored in the rat are discussed. The review also focuses more particularly on two areas in which extensive work has been done: human genetic diseases and immune system analysis. Models are thoroughly described in these two areas and highlight the competitive advantages of rat models over available corresponding mouse versions. The objective of this review is to provide a comprehensive description of the advantages and potential of rat models for addressing specific scientific questions and to characterize the best genome-engineering tools for developing new projects.

Network-Based Analysis of Cognitive Impairment and Memory Deficits from Transcriptome Data

Aging is an inevitable process that negatively affects all living organisms and their vital functions. The brain is one of the most important organs in living beings and is primarily impacted by aging. The molecular mechanisms of learning, memory and cognition are altered over time, and the impairment in these mechanisms can lead to neurodegenerative diseases. Transcriptomics can be used to study these impairments to acquire more detailed information on the affected molecular mechanisms. Here we analyzed learning- and memory-related transcriptome data by mapping it on the organism-specific protein-protein interactome network. Subnetwork discovery algorithms were applied to discover highly dysregulated subnetworks, which were complemented with co-expression-based interactions. The functional analysis shows that the identified subnetworks are enriched with genes having roles in synaptic plasticity, gliogenesis, neurogenesis and cognition, which are reported to be related to memory and learning. With a detailed analysis, we show that the results from different subnetwork discovery algorithms or from different transcriptomic datasets can be successfully reconciled, leading to a memory-learning network that sheds light on the molecular mechanisms behind aging and memory-related impairments.