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

AI-driven Discovery of Morphomolecular Signatures in Toxicology

Early identification of drug toxicity is essential yet challenging in drug development. At the preclinical stage, toxicity is assessed with histopathological examination of tissue sections from animal models to detect morphological lesions. To complement this analysis, toxicogenomics is increasingly employed to understand the mechanism of action of the compound and ultimately identify lesion-specific safety biomarkers for which in vitro assays can be designed. However, existing works that aim to identify morphological correlates of expression changes rely on qualitative or semi-quantitative morphological characterization and remain limited in scale or morphological diversity. Artificial intelligence (AI) offers a promising approach for quantitatively modeling this relationship at an unprecedented scale. Here, we introduce GEESE, an AI model designed to impute morphomolecular signatures in toxicology data. Our model was trained to predict 1,536 gene targets on a cohort of 8,231 hematoxylin and eosin-stained liver sections from Rattus norvegicus across 127 preclinical toxicity studies. The model, evaluated on 2,002 tissue sections from 29 held-out studies, can yield pseudo-spatially resolved gene expression maps, which we correlate with six key drug-induced liver injuries (DILI). From the resulting 25 million lesion-expression pairs, we established quantitative relations between up and downregulated genes and lesions. Validation of these signatures against toxicogenomic databases, pathway enrichment analyses, and human hepatocyte cell lines asserted their relevance. Overall, our study introduces new methods for characterizing toxicity at an unprecedented scale and granularity, paving the way for AI-driven discovery of toxicity biomarkers.

Effect of Short-Term Restraint Stress on the Hypothalamic Transcriptome Profiles of Rats with Inherited Stress-Induced Arterial Hypertension (ISIAH) and Normotensive Wistar Albino Glaxo (WAG) Rats

Emotional stress is one of the health risk factors in the modern human lifestyle. Stress exposure can provoke the manifestation of various pathological conditions, one of which is a sharp increase in the blood pressure level. In the present study, we analyzed changes in the transcriptome profiles of the hypothalamus of hypertensive ISIAH and normotensive WAG rats exposed to a single short-term restraint stress (the rat was placed in a tight wire-mesh cage for 2 h). This type of stress can be considered emotional stress. The functional annotation of differentially expressed genes allowed us to identify the most significantly altered biological processes in the hypothalamus of hypertensive and normotensive rats. The study made it possible to identify a group of genes that describe a general response to stress, independent of the rat genotype, as well as a hypothalamic response to stress specific to each strain. The alternatively changing expression of the Npas4 (neuronal PAS domain protein 4) gene, which is downregulated in the hypothalamus of the control WAG rats and induced in the hypothalamus of hypertensive ISIAH rats, is suggested to be the key event for understanding inter-strain differences in the hypothalamic response to stress. The stress-dependent ISIAH strain-specific induction of Fos and Jun gene transcription may play a crucial role in neuronal activation in this rat strain. The data obtained can be potentially useful in the selection of molecular targets for the development of pharmacological approaches to the correction of stress-induced pathologies related to neuronal excitability, taking into account the hypertensive status of the patients.

Identifying Molecular Pathophysiology and Potential Therapeutic Options in Iatrogenic Tracheal Stenosis

INTRODUCTION

While most patients with iatrogenic tracheal stenosis (ITS) respond to endoscopic ablative procedures, approximately 15% experience a recalcitrant, recurring disease course that is resistant to conventional management. We aimed to explore genetic profiles of patients with recalcitrant ITS to understand underlying pathophysiology and identify novel therapeutic options.

METHODS

We collected 11 samples of granulation tissue from patients with ITS and performed RNA sequencing. We identified the top 10 most highly up- and down-regulated genes and cellular processes that these genes corresponded to. For the most highly dysregulated genes, we identified potential therapeutic options that favorably regulate their expression.

RESULTS

The dysregulations in gene expression corresponded to hyperkeratinization (upregulation of genes involved in keratin production and keratinocyte differentiation) and cellular proliferation (downregulation of cell cycle regulating and pro-apoptotic genes). Genes involved in retinoic acid (RA) metabolism and signaling were dysregulated in a pattern suggesting local cellular RA deficiency. Consequently, RA also emerged as the most promising potential therapeutic option for ITS, as it favorably regulated seven of the ten most highly dysregulated genes.

CONCLUSION

This is the first study to characterize the role of hyperkeratinization and dysregulations in RA metabolism and signaling in the disease pathophysiology. Given the ability of RA to favorably regulate key genes involved in ITS, future studies must explore its efficacy as a potential therapeutic option for patients with recalcitrant ITS.

DNA methylation and type 2 diabetes: a systematic review

OBJECTIVE

DNA methylation influences gene expression and function in the pathophysiology of type 2 diabetes mellitus (T2DM). Mapping of T2DM-associated DNA methylation could aid early detection and/or therapeutic treatment options for diabetics.

DESIGN

A systematic literature search for associations between T2DM and DNA methylation was performed. Prospero registration ID: CRD42020140436.

METHODS

PubMed and ScienceDirect databases were searched (till October 19, 2023). Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines and New Castle Ottawa scale were used for reporting the selection and quality of the studies, respectively.

RESULT

Thirty-two articles were selected. Four of 130 differentially methylated genes in blood, adipose, liver or pancreatic islets (TXNIP, ABCG1, PPARGC1A, PTPRN2) were reported in > 1 study. TXNIP was hypomethylated in diabetic blood across ethnicities. Gene enrichment analysis of the differentially methylated genes highlighted relevant disease pathways (T2DM, type 1 diabetes and adipocytokine signaling). Three prospective studies reported association of methylation in IGFBP2, MSI2, FTO, TXNIP, SREBF1, PHOSPHO1, SOCS3 and ABCG1 in blood at baseline with incident T2DM/hyperglycemia. Sex-specific differential methylation was reported only for HOOK2 in visceral adipose tissue (female diabetics: hypermethylated, male diabetics: hypomethylated). Gene expression was inversely associated with methylation status in 8 studies, in genes including ABCG1 (blood), S100A4 (adipose tissue), PER2 (pancreatic islets), PDGFA (liver) and PPARGC1A (skeletal muscle).

CONCLUSION

This review summarizes available evidence for using DNA methylation patterns to unravel T2DM pathophysiology. Further validation studies in diverse populations will set the stage for utilizing this knowledge for identifying early diagnostic markers and novel druggable pathways.

Updates to the Alliance of Genome Resources central infrastructure

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, Caenorhabditis elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and application programming interfaces (APIs). Here, we focus on developments over the last 2 years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse software. We describe our progress toward a central persistent database to support curation, the data modeling that underpins harmonization, and progress toward a state-of-the-art literature curation system with integrated artificial intelligence and machine learning (AI/ML).

The genetic architecture of complete blood counts in lactating Holstein dairy cows

Complete blood counts (CBCs) measure the abundance of individual immune cells, red blood cells, and related measures such as platelets in circulating blood. These measures can indicate the health status of an animal; thus, baseline circulating levels in a healthy animal may be related to the productive life, resilience, and production efficiency of cattle. The objective of this study is to determine the heritability of CBC traits and identify genomic regions that are associated with CBC measurements in lactating Holstein dairy cattle. The heritability of CBCs was estimated using a Bayes C0 model. The study population consisted of 388 cows with genotypes at roughly 75,000 markers and 16 different CBC phenotypes taken at one to three time points (n = 33, 131, and 224 for 1, 2, and 3 time points, respectively). Heritabilities ranged from 0.00 ± 0.00 (red cell distribution width) to 0.68 ± 0.06 (lymphocytes). A total of 96 different 1-Mb windows were identified that explained more than 1% of the genetic variance for at least one CBC trait, with 10 windows explaining more than 1% of the genetic variance for two or more traits. Multiple genes in the identified regions have functions related to immune response, cell differentiation, anemia, and disease. Positional candidate genes include RAD52 motif-containing protein 1 (RDM1), which is correlated with the degree of immune infiltration of immune cells, and C-X-C motif chemokine ligand 12 (CXCL12), which is critically involved in neutrophil bone marrow storage and release regulation and enhances neutrophil migration. Since animal health directly impacts feed intake, understanding the genetics of CBCs may be useful in identifying more disease-resilient and feed-efficient dairy cattle. Identification of genes responsible for variation in CBCs will also help identify the variability in how dairy cattle defend against illness and injury.

Species-specific responses during Seoul orthohantavirus infection in human and rat lung microvascular endothelial cells

Seoul orthohantavirus (SEOV) is a rat-borne zoonotic virus that is transmitted via inhalation of aerosolized infectious excreta, and can cause hemorrhagic fever with renal syndrome (HFRS) in humans worldwide. In rats, SEOV predominantly exists as a persistent infection in the absence of overt clinical signs. Lack of disease in rats is attributed to downregulation of pro-inflammatory and upregulation of regulatory host responses. As lung microvascular endothelial cells (LMECs) represent a primary target of infection in both human and rats, infections in these cells provide a unique opportunity to study the central role of LMECs in the dichotomy between pathogenicity in both species. In this study, host responses to SEOV infection in primary human and rat LMECs were directly compared on a transcriptional level. As infection of rat LMECs was more efficient than human LMECs, the majority of anti-viral defense responses were observed earlier in rat LMECs. Most prominently, SEOV-induced processes in both species included responses to cytokine stimulus, negative regulation of innate immune responses, responses to type I and II interferons, regulation of pattern recognition receptor signaling and MHC-I signaling. However, over time, in the rat LMECs, responses shifted from an anti-viral state towards a more immunotolerant state displayed by a PD-L1, B2M-, JAK2-focused interaction network aiding in negative regulation of cytotoxic CD8-positive T cell activation. This suggests a novel mechanism by which species-specific orthohantavirus-induced endothelium and T cell crosstalk may play a crucial role in the development of acute disease in humans and persistence in rodents.

Hepatotropic activity of a betulonic acid based compound

Using the model of cyclophosphamide (CP)-induced immunosuppression in C57BL/6 mice, the hepatotropic effects of a conjugate of betulonic acid with 9-(4-methylpiperazin-1-ylmethyl)-2-(1,2,3-triazolyl) oreozelone (BABC) have been studied. In the liver of treated animals the expression of genes for cytochromes (CYP 1A1, CYP 1A2, CYP 3A44, CYP 2B10, CYP 2C29, CYP 17A1), PPARA, and cytokines (TNF-α, IL-1β, IL-12α, IL-10) and the relative levels of NF-κB p65, GST-π, and NAT-1 proteins were determined. On day six after administration of the compound and CP to animals a significant (3.2-fold) increase in the expression of the CYP 2B10 as compared to the control group was observed. Treatment of mice with the compound and CP also caused a 2.4-fold increase in the mRNA level of the pro-inflammatory TNF-α gene as compared to the group of animals receiving CP. Administration of the studied compound to intact animals was accompanied by a 2.5-fold increase in the IL-1β expression and a 1.8-fold decrease in the IL-10 expression as compared to the control group. An increase in the expression of pro-inflammatory cytokine genes in the liver of animals treated with the compound was accompanied by an increase in the content of NF-κB p65 (by 1.6 times), as well as an increase in the relative amount of NAT-1 protein (by 2.7 times) as compared to control animals.

Identification of sleep and circadian alternative polyadenylation sites associated with APA-linked human brain disorders

Sleep and circadian rhythm disruptions are comorbid features of many pathologies and can negatively influence numerous health conditions, including degenerative diseases, metabolic illnesses, cancer, and various neurological disorders. Genetic association studies linking sleep and circadian disturbances with disease susceptibility have mainly focused on changes in gene expression due to mutations, such as single-nucleotide polymorphisms. Thus, associations between sleep and/or circadian rhythm and alternative polyadenylation (APA), particularly in the context of other health challenges, are largely undescribed. APA is a process that generates various transcript isoforms from the same gene, resulting in effects on mRNA translation, stability, localization, and subsequent function. Here, we have identified unique APAs in rat brain that exhibit time-of-day-dependent oscillations in expression as well as APAs that are altered by sleep deprivation and the subsequent recovery period. Genes affected by APA usage include Mapt/Tau, Ntrk2, Homer1A, Sin3band Sorl. Sorl1 has two APAs which cycle with a 24 h period, one additional APA cycles with a 12 h period and one more that is reduced during recovery sleep. Finally, we compared sleep- or circadian-associated APAs with recently described APA-linked brain disorder susceptibility genes and found 46 genes in common.

A3D Model Organism Database (A3D-MODB): a database for proteome aggregation predictions in model organisms

Protein aggregation has been associated with aging and different pathologies and represents a bottleneck in the industrial production of biotherapeutics. Numerous past studies performed in Escherichia coli and other model organisms have allowed to dissect the biophysical principles underlying this process. This knowledge fuelled the development of computational tools, such as Aggrescan 3D (A3D) to forecast and re-design protein aggregation. Here, we present the A3D Model Organism Database (A3D-MODB) http://biocomp.chem.uw.edu.pl/A3D2/MODB, a comprehensive resource for the study of structural protein aggregation in the proteomes of 12 key model species spanning distant biological clades. In addition to A3D predictions, this resource incorporates information useful for contextualizing protein aggregation, including membrane protein topology and structural model confidence, as an indirect reporter of protein disorder. The database is openly accessible without any need for registration. We foresee A3D-MOBD evolving into a central hub for conducting comprehensive, multi-species analyses of protein aggregation, fostering the development of protein-based solutions for medical, biotechnological, agricultural and industrial applications.

The alliance of genome resources: transforming comparative genomics

Comparing genomic and biological characteristics across multiple species is essential to using model systems to investigate the molecular and cellular mechanisms underlying human biology and disease and to translate mechanistic insights from studies in model organisms for clinical applications. Building a scalable knowledge commons platform that supports cross-species comparison of rich, expertly curated knowledge regarding gene function, phenotype, and disease associations available for model organisms and humans is the primary mission of the Alliance of Genome Resources (the Alliance). The Alliance is a consortium of seven model organism knowledgebases (mouse, rat, yeast, nematode, zebrafish, frog, fruit fly) and the Gene Ontology resource. The Alliance uses a common set of gene ortholog assertions as the basis for comparing biological annotations across the organisms represented in the Alliance. The major types of knowledge associated with genes that are represented in the Alliance database currently include gene function, phenotypic alleles and variants, human disease associations, pathways, gene expression, and both protein-protein and genetic interactions. The Alliance has enhanced the ability of researchers to easily compare biological annotations for common data types across model organisms and human through the implementation of shared programmatic access mechanisms, data-specific web pages with a unified "look and feel", and interactive user interfaces specifically designed to support comparative biology. The modular infrastructure developed by the Alliance allows the resource to serve as an extensible "knowledge commons" capable of expanding to accommodate additional model organisms.

Y and Mitochondrial Chromosomes in the Heterogeneous Stock Rat Population

Genome-wide association studies typically evaluate the autosomes and sometimes the X Chromosome, but seldom consider the Y or mitochondrial Chromosomes. We genotyped the Y and mitochondrial chromosomes in heterogeneous stock rats (Rattus norvegicus), which were created in 1984 by intercrossing eight inbred strains and have subsequently been maintained as an outbred population for 100 generations. As the Y and mitochondrial Chromosomes do not recombine, we determined which founder had contributed these chromosomes for each rat, and then performed association analysis for all complex traits (n=12,055; intersection of 12,116 phenotyped and 15,042 haplotyped rats). We found the eight founders had 8 distinct Y and 4 distinct mitochondrial Chromosomes, however only two of each were observed in our modern heterogeneous stock rat population (Generations 81-97). Despite the unusually large sample size, the p-value distribution did not deviate from expectations; there were no significant associations for behavioral, physiological, metabolome, or microbiome traits after correcting for multiple comparisons. However, both Y and mitochondrial Chromosomes were strongly associated with expression of a few genes located on those chromosomes, which provided a positive control. Our results suggest that within modern heterogeneous stock rats there are no Y and mitochondrial Chromosomes differences that strongly influence behavioral or physiological traits. These results do not address other ancestral Y and mitochondrial Chromosomes that do not appear in modern heterogeneous stock rats, nor do they address effects that may exist in other rat populations, or in other species.

Updates to the Alliance of Genome Resources Central Infrastructure Alliance of Genome Resources Consortium

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively-studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, C. elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and APIs. Here we focus on developments over the last two years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse. We describe our progress towards a central persistent database to support curation, the data modeling that underpins harmonization, and progress towards a state-of-the art literature curation system with integrated Artificial Intelligence and Machine Learning (AI/ML).

A Cell Atlas of Thoracic Aortic Perivascular Adipose Tissue: a focus on mechanotransducers

Perivascular adipose tissue (PVAT) is increasingly recognized for its function in mechanotransduction. To examine the cell-specificity of recognized mechanotransducers we used single nuclei RNA sequencing (snRNAseq) of the thoracic aorta PVAT (taPVAT) from male Dahl SS rats compared to subscapular brown adipose tissue (BAT). Approximately 30,000 nuclei from taPVAT and BAT each were characterized by snRNAseq, identifying 8 major cell types expected and one unexpected (nuclei with oligodendrocyte marker genes). Cell-specific differential gene expression analysis between taPVAT and BAT identified up to 511 genes (adipocytes) with many (≥20%) being unique to individual cell types. Piezo1 was the most highly, widely expressed mechanotransducer. Presence of PIEZO1 in the PVAT was confirmed by RNAscope® and IHC; antagonism of PIEZO1 impaired the PVAT's ability to hold tension. Collectively, the cell compositions of taPVAT and BAT are highly similar, and PIEZO1 is likely a mechanotransducer in taPVAT.

Advances in cellular resolution microscopy for brain imaging in rats

Rats are used in neuroscience research because of their physiological similarities with humans and accessibility as model organisms, trainability, and behavioral repertoire. In particular, rats perform a wide range of sophisticated social, cognitive, motor, and learning behaviors within the contexts of both naturalistic and laboratory environments. Further progress in neuroscience can be facilitated by using advanced imaging methods to measure the complex neural and physiological processes during behavior in rats. However, compared with the mouse, the rat nervous system offers a set of challenges, such as larger brain size, decreased neuron density, and difficulty with head restraint. Here, we review recent advances in in vivo imaging techniques in rats with a special focus on open-source solutions for calcium imaging. Finally, we provide suggestions for both users and developers of in vivo imaging systems for rats.