A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages

Lipid nanoparticle-mediated base-editing of the Hao1 gene achieves sustainable primary hyperoxaluria type 1 therapy in rats

Primary hyperoxaluria type 1 (PH1) is a severe hereditary disease, leading to the accumulation of oxalate in multiple organs, particularly the kidney. Hydroxyacid oxidase 1 (HAO1), a pivotal gene involved in oxalate production, is an approved target for the treatment of PH1. In this study, we demonstrated the discovery of several novel therapeutic sites of the Hao1 gene and the efficient editing of Hao1 c.290-2 A in vivo with lipid nanoparticles (LNP) delivered adenine base editing (ABE) mRNA. A single infusion of LNP-ABE resulted in a near-complete knockout of Hao1 in the liver, leading to the sustainable normalization of urinary oxalate (for at least 6 months) and complete rescue of the patho-physiology in PH1 rats. Additionally, a significant correlation between Hao1 editing efficiency and urinary oxalate levels was observed and over 60% Hao1 editing efficiency was required to achieve the normalization of urinary oxalate in PH1 rats. These findings suggest that the LNP-mediated base-editing of Hao1 c.290-2 A is an efficient and safe approach to PH1 therapy, highlighting its potential utility in clinical settings.

Cardiac Aging in the Multi-Omics Era: High-Throughput Sequencing Insights

Cardiovascular diseases are a leading cause of mortality worldwide, and the risks of both developing a disease and receiving a poor prognosis increase with age. With increasing life expectancy, understanding the mechanisms underlying heart aging has become critical. Traditional techniques have supported research into finding the physiological changes and hallmarks of cardiovascular aging, including oxidative stress, disabled macroautophagy, loss of proteostasis, and epigenetic alterations, among others. The advent of high-throughput multi-omics techniques offers new perspectives on the molecular mechanisms and cellular processes in the heart, guiding the development of therapeutic targets. This review explores the contributions and characteristics of these high-throughput techniques to unraveling heart aging. We discuss how different high-throughput omics approaches, both alone and in combination, produce robust and exciting new findings and outline future directions and prospects in studying heart aging in this new era.

Expression profiling of cochlear genes uncovers sex-based cellular function in mouse cochleae

Sex is a pivotal biological factor that significantly impacts tissue homeostasis and disease susceptibility. In the auditory system, sex differences have been observed in cochlear physiology and responses to pathological conditions. However, the underlying molecular mechanisms responsible for these differences remain elusive. The current research explores the differences in gene expression profiles in the cochlea between male and female mice, aiming to understand the functional implication of sex-biased gene expression in each sex. Using RNA-sequencing analysis on cochlear tissues obtained from male and female mice, we identified a significant number of genes exhibiting sex-biased expression differences. While some of these differentially expressed genes are located on sex chromosomes, most are found on autosomal chromosomes. Further bioinformatic analysis revealed that these genes are involved in several key cellular functions. In males, these genes are notably linked to oxidative phosphorylation and RNA synthesis and processing, suggesting their involvement in mitochondrial energy production and regulatory control of gene expression. In contrast, sex-biased genes are associated with mechano-transduction and synaptic transmission within female cochleae. Collectively, our study provides valuable insights into the molecular differences between the sexes and emphasizes the need for future research to uncover their functional implications and relevance to auditory health and disease development.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

New mutation within a common haplotype is associated with calf muscle weakness in Holsteins

A recessive haplotype resulting in elevated calf mortality but with apparent incomplete penetrance was previously linked to the end of chromosome 16 (78.7-80.7 Mbp). Genotype analysis of 5.6 million Holsteins indicated that the haplotype was common and traced back to 1952, with a key ancestor born in 1984 (HOUSA1964484, Southwind) identified from chip genotypes as homozygous for the suspect haplotype. Sequence data from Southwind (an affected calf) and the sire of the affected calf were scanned for candidate mutations. A missense mutation with a deleterious projected impact at 79,613,592 bp was homozygous in the affected calf and heterozygous in the calf's sire and Southwind. Sequence data available from the Cooperative Dairy DNA Repository for 299 other Holsteins indicated a 97% concordance with the haplotype and an 89% call rate. The exon amino acid sequence appears to be broadly conserved in the CACNA1S gene, and mutations in humans and mice can cause phenotypes of temporary or permanent paralysis analogous to those in calves with the haplotype causing muscle weakness (HMW). Improved methods for using pedigree to track new mutations within existing haplotypes were developed and applied to the haplotypes for both muscle weakness and Holstein cholesterol deficiency (HCD). For HCD, concordance of the gene test with its haplotype status was greatly improved. For both defects, haplotype status was matched to heifer livability records for 558,000 calves. For HMW, only 46 heifers with livability records were homozygous and traced only to Southwind on both sides. Of those, 52% died before 18 mo at an average age of 1.7 ± 1.6 mo, but that death rate may be underestimated if only healthier calves were genotyped. The death rate was 2.4% for noncarriers. Different reporting methods or dominance effects may be needed to include HMW and other partially lethal effects in selection and mating. Direct tests are needed for new mutations within existing common haplotypes because tracking can be difficult even with accurate pedigrees when the original haplotype has a high frequency.

Effects of age and diet consistency on the expression of myosin heavy-chain isoforms on jaw-closing and jaw-opening muscles in a rat model

BACKGROUND

Skeletal craniofacial morphology can be influenced by changes in masticatory muscle function, which may also change the functional profile of the muscles.

OBJECTIVES

To investigate the effects of age and functional demands on the expression of Myosin Heavy-Chain (MyHC) isoforms in representative jaw-closing and jaw-opening muscles, namely the masseter and digastric muscles respectively.

METHODS

Eighty-four male Wistar rats were divided into four age groups, namely an immature (n = 12; 4-week-old), early adult (n = 24; 16-week-old), adult (n = 24; 26-week-old) and mature adult (n = 24; 38-week-old) group. The three adult groups were divided into two subgroups each based on diet consistency; a control group fed a standard (hard) diet, and an experimental group fed a soft diet. Rats were sacrificed, and masseter and digastric muscles dissected. Real-time quantitative polymerase chain reaction was used to compare the mRNA transcripts of the MyHC isoforms-Myh7 (MyHC-I), Myh2 (MyHC-IIa), Myh4 (MyHC-IIb) and Myh1 (MyHC-IIx)-of deep masseter and digastric muscles.

RESULTS

In the masseter muscle, hypofunction increases Myh1 (26, 38 weeks; p < .0001) but decreases Myh4 (26 weeks; p = .046) and Myh2 (26 weeks; p < .0001) expression in adult rats. In the digastric muscle, hypofunction increases Myh1 expression in the mature adult rats (38 weeks; p < .0001), while Myh2 expression decreases in adult rats (26 weeks; p = .021) as does Myh4 (26 weeks; p = .001). Myh7 expression is increased in the digastric muscle of mature adult rats subjected to hypofunction (38 weeks; p = <.0001), while it is very weakly expressed in the masseter.

CONCLUSION

In jaw-opening and jaw-closing muscles, differences in myosin expression between hard- and soft-diet-fed rats become evident in adulthood, suggesting that long-term alteration of jaw function is associated with changes in the expression of MyHC isoforms and potential fibre remodelling. This may give insight into the role of function on masticatory muscles and the resultant craniofacial morphology.

Four-week inhibition of the renin-angiotensin system in spontaneously hypertensive rats results in persistently lower blood pressure with reduced kidney renin and changes in expression of relevant gene networks

AIMS

Prevention of human hypertension is an important challenge and has been achieved in experimental models. Brief treatment with renin-angiotensin system (RAS) inhibitors permanently reduces the genetic hypertension of the spontaneously hypertensive rat (SHR). The kidney is involved in this fascinating phenomenon, but relevant changes in gene expression are unknown.

METHODS AND RESULTS

In SHR, we studied the effect of treatment between 10 and 14 weeks of age with the angiotensin receptor blocker, losartan, or the angiotensin-converting enzyme inhibitor, perindopril [with controls for non-specific effects of lowering blood pressure (BP)], on differential RNA expression, DNA methylation, and renin immunolabelling in the kidney at 20 weeks of age. RNA sequencing revealed a six-fold increase in renin gene (Ren) expression during losartan treatment (P < 0.0001). Six weeks after losartan, arterial pressure remained lower (P = 0.006), yet kidney Ren showed reduced expression by 23% after losartan (P = 0.03) and by 43% after perindopril (P = 1.4 × 10-6) associated with increased DNA methylation (P = 0.04). Immunolabelling confirmed reduced cortical renin after earlier RAS blockade (P = 0.002). RNA sequencing identified differential expression of mRNAs, miRNAs, and lncRNAs with evidence of networking and co-regulation. These included 13 candidate genes (Grhl1, Ammecr1l, Hs6st1, Nfil3, Fam221a, Lmo4, Adamts1, Cish, Hif3a, Bcl6, Rad54l2, Adap1, Dok4), the miRNA miR-145-3p, and the lncRNA AC115371. Gene ontogeny analyses revealed that these networks were enriched with genes relevant to BP, RAS, and the kidneys.

CONCLUSION

Early RAS inhibition in SHR resets genetic pathways and networks resulting in a legacy of reduced Ren expression and BP persisting for a minimum of 6 weeks.

Sex-specific developmental gene expression atlas unveils dimorphic gene networks in C. elegans

Sex-specific traits and behaviors emerge during development by the acquisition of unique properties in the nervous system of each sex. However, the genetic events responsible for introducing these sex-specific features remain poorly understood. In this study, we create a comprehensive gene expression atlas of pure populations of hermaphrodites and males of the nematode Caenorhabditis elegans across development. We discover numerous differentially expressed genes, including neuronal gene families like transcription factors, neuropeptides, and G protein-coupled receptors. We identify INS-39, an insulin-like peptide, as a prominent male-biased gene expressed specifically in ciliated sensory neurons. We show that INS-39 serves as an early-stage male marker, facilitating the effective isolation of males in high-throughput experiments. Through complex and sex-specific regulation, ins-39 plays pleiotropic sexually dimorphic roles in various behaviors, while also playing a shared, dimorphic role in early life stress. This study offers a comparative sexual and developmental gene expression database for C. elegans. Furthermore, it highlights conserved genes that may underlie the sexually dimorphic manifestation of different human diseases.

Circular RNA expression profiles and functional predication after restraint stress in the amygdala of rats

Prolonged or repeated exposure to stress elevates the risk of various psychological diseases, many of which are characterized by central nervous system dysfunction. Recent studies have demonstrated that circular RNAs (circRNAs) are highly abundant in the mammalian brain. Although their precise expression and function remain unknown, they have been hypothesized to regulate transcriptional and post-transcriptional gene expression. In this investigation, we comprehensively analyzed whether restraint stress for 2 days altered the circRNA expression profile in the amygdala of male rats. The impact of restraint stress on behavior was evaluated using an elevated plus maze and open field test. Serum corticosterone levels were measured using an enzyme-linked immunosorbent assay. A total of 10,670 circRNAs were identified using RNA sequencing. Ten circRNAs were validated by reverse transcription and quantitative polymerase chain reaction analysis. Gene ontology and Kyoto encyclopedia of genes and genomes pathway analyzes supported the notion that genes associated with differentially expressed circRNAs are primarily implicated in neuronal activity and neurotransmitter transport. Moreover, the three differentially expressed circRNAs showed high specificity in the amygdala. Overall, these findings indicate that differentially expressed circRNAs are highly enriched in the amygdala and offer a potential direction for further research on restraint stress.

Comprehensive transcriptome analysis reveals altered mRNA splicing and post-transcriptional changes in the aged mouse brain

A comprehensive understanding of molecular changes during brain aging is essential to mitigate cognitive decline and delay neurodegenerative diseases. The interpretation of mRNA alterations during brain aging is influenced by the health and age of the animal cohorts studied. Here, we carefully consider these factors and provide an in-depth investigation of mRNA splicing and dynamics in the aging mouse brain, combining short- and long-read sequencing technologies with extensive bioinformatic analyses. Our findings encompass a spectrum of age-related changes, including differences in isoform usage, decreased mRNA dynamics and a module showing increased expression of neuronal genes. Notably, our results indicate a reduced abundance of mRNA isoforms leading to nonsense-mediated RNA decay and suggest a regulatory role for RNA-binding proteins, indicating that their regulation may be altered leading to the reshaping of the aged brain transcriptome. Collectively, our study highlights the importance of studying mRNA splicing events during brain aging.

Evaluating the efficacy of protein quantification methods on membrane proteins

Protein quantification is an important tool for a wide range of biological applications. The most common broadscale methods include the Lowry, bicinchoninic acid (BCA), and Coomassie Bradford assays. Despite their wide applicability, the mechanisms of action imply that these methods may not be ideal for large transmembrane proteins due to the proteins' integration in the plasma membrane. Here, we investigate this problem by assessing the efficacy and applicability of these three common protein quantification methods on a candidate transmembrane protein - the Na,K-ATPase (NKA). We compared these methods to an ELISA, which we newly developed and describe here for the quantification of NKA. The use of a relative standard curve allows this ELISA to be easily adapted to other proteins and across the animal kingdom. Our results revealed that the three conventional methods significantly underestimate the concentration of NKA compared to the ELISA. Further, by applying the protein concentrations determined by the different methods to in vitro assays, we found that variation in the resulting data was consistently low when the assay reactions were prepared based on concentrations determined from the ELISA. Thus, when target protein concentrations vary across samples, the conventional quantification methods cannot produce reliable results in downstream applications. In contrast, the ELISA we describe here consistently provides robust results.

Dopamine receptor D2 regulates genes involved in germ cell movement and sperm motility in rat testes†

The function of dopamine receptor D2 (D2R) is well associated with sperm motility; however, the physiological role of D2R present on testicular cells remains elusive. The aim of the present study is to delineate the function of testicular D2R. Serum dopamine levels were found to decrease with age, whereas testicular D2R expression increased. In rat testicular sections, D2R immunolabeling was observed in interstitial cells, spermatogonia, spermatocytes and mature elongated spermatids, whereas tyrosine hydroxylase immunolabeling was selectively detected in Leydig cells. In vitro seminiferous tubule culture following bromocriptine (D2R agonist) treatment resulted in decreased cAMP levels. Microarray identified 1077 differentially expressed genes (511 up-regulated, 566 down-regulated). The majority of differentially expressed genes were present in post-meiotic cells including early and late spermatids, and sperm. Gene ontology elucidated processes related to extra-cellular matrix to be enriched and was supported by differential expression of various collagens and laminins, thereby indicating a role of dopamine in extra-cellular matrix integrity and transport of spermatids across the seminiferous epithelium. Gene ontology and enrichment map also highlighted cell/sperm motility to be significantly enriched. Therefore, genes involved in sperm motility functions were further validated by RT-qPCR. Seven genes (Akap4, Ccnyl1, Iqcf1, Klc3, Prss55, Tbc1d21, Tl18) were significantly up-regulated, whereas four genes (Dnah1, Dnah5, Clxn, Fsip2) were significantly down-regulated by bromocriptine treatment. The bromocriptine-stimulated reduction in seminiferous tubule cyclic AMP and associated changes in spermatid gene expression suggests that dopamine regulates both spermatogenesis and spermiogenesis within the seminiferous epithelium, and spermatozoa motility following spermiation, as essential processes for fertility.

Role of increased IGFBP2 in trophoblast cell proliferation and recurrent spontaneous abortion development: A pilot study

Recurrent spontaneous abortion (RSA) is a serious condition that adversely affects women's health. Differentially expressed proteins (DEPs) in plasma of patients experiencing RSA is helpful to find new therapeutic targets and identified with mass spectrometry. In 57 DEPs, 21 were upregulated and 36 were downregulated in RSA. Gene ontology analyses indicated that identified DEPs were associated with cell proliferation, including significantly downregulated insulin-like growth factor binding protein 2 (IGFBP2). Immunohistochemical result using clinical decidual tissues also showed that IGFBP2 expression was significantly decreased in RSA trophoblasts. Cell proliferation assay indicated that IGFBP2 treatment increased the proliferation and mRNA expressions of PCNA and Ki67 in trophoblast cells. Transcriptome sequencing experiments and Kyoto Encyclopedia of Genes and Genomes analyses revealed that gene expression for components in PI3K-Akt pathway in trophoblasts was significantly upregulated following IGFBP2 treatment. To confirm bioinformatics findings, we did cell-based experiments and found that treatment of inhibitors for insulin-like growth factor (IGF)-1 receptor-PI3K-Akt pathway significantly reduced IGFBP2-induced trophoblast cell proliferation and mRNA expressions of PCNA and Ki67. Our findings suggest that IGFBP2 may increase trophoblast proliferation through the PI3K-Akt signaling pathway to affect pregnancy outcomes and that IGFBP2 may be a new target for future research and treatment of RSA.

The potential impact of CYP and UGT drug-metabolizing enzymes on brain target site drug exposure

Drug metabolism is one of the critical determinants of drug disposition throughout the body. While traditionally associated with the liver, recent research has unveiled the presence and functional significance of drug-metabolizing enzymes (DMEs) within the brain. Specifically, cytochrome P-450 enzymes (CYPs) and UDP-glucuronosyltransferases (UGTs) enzymes have emerged as key players in drug biotransformation within the central nervous system (CNS). This comprehensive review explores the cellular and subcellular distribution of CYPs and UGTs within the CNS, emphasizing regional expression and contrasting profiles between the liver and brain, humans and rats. Moreover, we discuss the impact of species and sex differences on CYPs and UGTs within the CNS. This review also provides an overview of methodologies for identifying and quantifying enzyme activities in the brain. Additionally, we present factors influencing CYPs and UGTs activities in the brain, including genetic polymorphisms, physiological variables, pathophysiological conditions, and environmental factors. Examples of CYP- and UGT-mediated drug metabolism within the brain are presented at the end, illustrating the pivotal role of these enzymes in drug therapy and potential toxicity. In conclusion, this review enhances our understanding of drug metabolism's significance in the brain, with a specific focus on CYPs and UGTs. Insights into the expression, activity, and influential factors of these enzymes within the CNS have crucial implications for drug development, the design of safe drug treatment strategies, and the comprehension of drug actions within the CNS. To that end, CNS pharmacokinetic (PK) models can be improved to further advance drug development and personalized therapy.

Rat microbial biogeography and age-dependent lactic acid bacteria in healthy lungs

The laboratory rat emerges as a useful tool for studying the interaction between the host and its microbiome. To advance principles relevant to the human microbiome, we systematically investigated and defined the multitissue microbial biogeography of healthy Fischer 344 rats across their lifespan. Microbial community profiling data were extracted and integrated with host transcriptomic data from the Sequencing Quality Control consortium. Unsupervised machine learning, correlation, taxonomic diversity and abundance analyses were performed to determine and characterize the rat microbial biogeography and identify four intertissue microbial heterogeneity patterns (P1-P4). We found that the 11 body habitats harbored a greater diversity of microbes than previously suspected. Lactic acid bacteria (LAB) abundance progressively declined in lungs from breastfed newborn to adolescence/adult, and was below detectable levels in elderly rats. Bioinformatics analyses indicate that the abundance of LAB may be modulated by the lung-immune axis. The presence and levels of LAB in lungs were further evaluated by PCR in two validation datasets. The lung, testes, thymus, kidney, adrenal and muscle niches were found to have age-dependent alterations in microbial abundance. The 357 microbial signatures were positively correlated with host genes in cell proliferation (P1), DNA damage repair (P2) and DNA transcription (P3). Our study established a link between the metabolic properties of LAB with lung microbiota maturation and development. Breastfeeding and environmental exposure influence microbiome composition and host health and longevity. The inferred rat microbial biogeography and pattern-specific microbial signatures could be useful for microbiome therapeutic approaches to human health and life quality enhancement.

SpliceProt 2.0: A Sequence Repository of Human, Mouse, and Rat Proteoforms

SpliceProt 2.0 is a public proteogenomics database that aims to list the sequence of known proteins and potential new proteoforms in human, mouse, and rat proteomes. This updated repository provides an even broader range of computationally translated proteins and serves, for example, to aid with proteomic validation of splice variants absent from the reference UniProtKB/SwissProt database. We demonstrate the value of SpliceProt 2.0 to predict orthologous proteins between humans and murines based on transcript reconstruction, sequence annotation and detection at the transcriptome and proteome levels. In this release, the annotation data used in the reconstruction of transcripts based on the methodology of ternary matrices were acquired from new databases such as Ensembl, UniProt, and APPRIS. Another innovation implemented in the pipeline is the exclusion of transcripts predicted to be susceptible to degradation through the NMD pathway. Taken together, our repository and its applications represent a valuable resource for the proteogenomics community.

Synthesis of N-Doped Graphene Photo-Catalyst for Photo-Assisted Charging of Li-Ion Oxygen Battery

In this work, nitrogen (N)-doped graphene film is synthesized, as a photo-catalyst, on one side of the copper foam by chemical vapor deposition and the copper foam is directly used as an electrode after porous Pd@rGO cathode loading to the other side of the foam for the photo-assisted charging of the Li-ion oxygen battery. The amount of urea (CO(NH2)2), which is used as N atom source, is optimized to get maximum photo-anodic currents from the n-type graphene films. The optical band gap and the valance band edge potential of the optimized N-doped graphene film are determined as 2.00 eV and 3.71 VLi+/Li, respectively. X-ray photoelectron spectra provided that the atomic percent of N atoms in the graphene film is 1.34% and the graphitic, pyrrolic and pyridinic N atom percentages are 54.01%, 42.20% and 3.79%, respectively. The photo-assisted charging tests indicated that the N-doped graphene film photo-catalyst reduced the charging potential significantly even at 1000 mA g-1 (0.1 mA cm-2) current density and improved the cyclic discharge-charge performance of the Li-ion oxygen battery considerably.

TAF7L regulates early stages of male germ cell development in the rat

Male germ cell development is dependent on the orchestrated regulation of gene networks. TATA-box binding protein associated factors (TAFs) facilitate interactions of TATA-binding protein with the TATA element, which is known to coordinate gene transcription during organogenesis. TAF7 like (Taf7l) is situated on the X chromosome and has been implicated in testis development. We examined the biology of TAF7L in testis development using the rat. Taf7l was prominently expressed in preleptotene to leptotene spermatocytes. To study the impact of TAF7L on the testis we generated a global loss-of-function rat model using CRISPR/Cas9 genome editing. Exon 3 of the Taf7l gene was targeted. A founder was generated possessing a 110 bp deletion within the Taf7l locus, which resulted in a frameshift and the premature appearance of a stop codon. The mutation was effectively transmitted through the germline. Deficits in TAF7L did not adversely affect pregnancy or postnatal survival. However, the Taf7l disruption resulted in male infertility due to compromised testis development and failed sperm production. Mutant germ cells suffer meiotic arrest at late zygotene/early pachynema stages, with defects in sex body formation. This testis phenotype was more pronounced than previously described for the subfertile Taf7l null mouse. We conclude that TAF7L is essential for male germ cell development in the rat.

Binding and sequestration of poison frog alkaloids by a plasma globulin

Alkaloids are important bioactive molecules throughout the natural world, and in many animals they serve as a source of chemical defense against predation. Dendrobatid poison frogs bioaccumulate alkaloids from their diet to make themselves toxic or unpalatable to predators. Despite the proposed roles of plasma proteins as mediators of alkaloid trafficking and bioavailability, the responsible proteins have not been identified. We use chemical approaches to show that a ~50 kDa plasma protein is the principal alkaloid-binding molecule in blood of poison frogs. Proteomic and biochemical studies establish this plasma protein to be a liver-derived alkaloid-binding globulin (ABG) that is a member of the serine-protease inhibitor (serpin) family. In addition to alkaloid-binding activity, ABG sequesters and regulates the bioavailability of 'free' plasma alkaloids in vitro. Unexpectedly, ABG is not related to saxiphilin, albumin, or other known vitamin carriers, but instead exhibits sequence and structural homology to mammalian hormone carriers and amphibian biliverdin-binding proteins. ABG represents a new small molecule binding functionality in serpin proteins, a novel mechanism of plasma alkaloid transport in poison frogs, and more broadly points toward serpins acting as tunable scaffolds for small molecule binding and transport across different organisms.

Transcriptomic analysis of mRNAs in human whole blood identified age-specific changes in healthy individuals

Older age is one of the most important shared risk factors for multiple chronic diseases, increasing the medical burden to contemporary societies. Current research focuses on identifying aging biomarkers to predict aging trajectories and developing interventions aimed at preventing and delaying the progression of multimorbidity with aging. Here, a transcriptomic changes analysis of whole blood genes with age was conducted. The age-related whole blood gene-expression profiling datasets were downloaded from the Gene Expression Omnibus (GEO) database. We screened the differentially expressed genes (DEGs) between healthy young and old individuals and performed functional enrichment analysis. Cytoscape with Cytohubba and MCODE was used to perform an interaction network of DEGs and identify hub genes. In addition, ROC curves and correlation analysis were used to evaluate the accuracy of hub genes. In total, we identified 29 DEGs between young and old samples that were enriched mainly in immunoglobulin binding and complex, humoral immune response, and immune response-activating signaling pathways. In combination with the PPI network and topological analysis, 4 hub genes (IGLL5, Jchain, POU2AF1, and Bach2) were identified. Pearson analysis showed that the expression changes of these hub genes were highly correlated with age. Among them, 3 hub genes (IGLL5, POU2AF1, and Bach2) were identified with good accuracy (AUC score > 0.7), indicating that these genes were the best indicators of age. Together, our results provided potential biomarkers IGLL5, POU2AF1, and Bach2 to identify individuals at high early risk of age-related disease to be targeted for early interventions and contribute to understanding the molecular mechanisms in the progression of aging.

Changes in glutamic oxaloacetic transaminase 2 during rat physiological and pathological cardiomyocyte hypertrophy

BACKGROUND

Physiological and pathological cardiomyocyte hypertrophy are important pathophysiological processes of adult congenital heart disease-associated ventricular hypertrophy. Glutamic oxaloacetic transaminase (GOT) is a vital marker of myocardial injury. This study aimed to investigate the changes in GOT levels during physiological and pathological cardiomyocyte hypertrophy in rats.

METHODS

RNA-seq analysis and colorimetric methods were used to evaluate the changes in GOT mRNA and activity, respectively. GOT2 protein expression was detected by western blotting and immunofluorescence. Hematoxylin-eosin and wheat germ agglutinin methods were used to observe changes in rat cardiomyocyte morphology.

RESULTS

In juvenile rat hearts, GOT mRNA expression and activity, and GOT2 protein level increased with age-related physiological cardiomyocyte hypertrophy; however, GOT2 protein level was reduced in hypoxia-induced pathological cardiomyocyte hypertrophy.

CONCLUSIONS

GOT2 may regulate physiological and pathological myocardial hypertrophy in rats. We speculated that the low GOT2 level contributed to the rapid occurrence of pathological cardiomyocyte hypertrophy, causing strong plasticity of right ventricular cardiomyocytes in the early postnatal period and heart failure in adulthood.

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.

Long-term α5 GABA A receptor negative allosteric modulator treatment reduces NMDAR-mediated neuronal excitation and maintains basal neuronal inhibition

α5 subunit-containing GABA type-A receptors (α5 GABAARs) are enriched in the hippocampus and play critical roles in neurodevelopment, synaptic plasticity, and cognition. α5 GABAAR preferring negative allosteric modulators (α5 NAMs) show promise mitigating cognitive impairment in preclinical studies of conditions characterized by excess GABAergic inhibition, including Down syndrome and memory deficits post-anesthesia. However, previous studies have primarily focused on acute application or single-dose α5 NAM treatment. Here, we measured the effects of chronic (7-day) in vitro treatment with L-655,708 (L6), a highly selective α5 NAM, on glutamatergic and GABAergic synapses in rat hippocampal neurons. We previously showed that 2-day in vitro treatment with L6 enhanced synaptic levels of the glutamate NMDA receptor (NMDAR) GluN2A subunit without modifying surface α5 GABAAR expression, inhibitory synapse function, or L6 sensitivity. We hypothesized that chronic L6 treatment would further increase synaptic GluN2A subunit levels while maintaining GABAergic inhibition and L6 efficacy, thus increasing neuronal excitation and glutamate-evoked intracellular calcium responses. Immunofluorescence experiments revealed that 7-day L6 treatment slightly increased the synaptic levels of gephyrin and surface α5 GABAARs. Functional studies showed that chronic α5 NAM treatment did not alter inhibition or α5 NAM sensitivity. Surprisingly, chronic L6 exposure decreased surface levels of GluN2A and GluN2B subunits, concurrent with reduced NMDAR-mediated neuronal excitation as seen by faster synaptic decay rates and reduced glutamate-evoked calcium responses. Together, these results show that chronic in vitro treatment with an α5 NAM leads to subtle homeostatic changes in inhibitory and excitatory synapses that suggest an overall dampening of excitability.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

RNA-Seq transcriptomic landscape profiling of spontaneously hypertensive rats treated with a sodium-glucose cotransporter 2 (SGLT2) inhibitor

BACKGROUND

Sodium-glucose co-transporter-2 inhibitor (SGLT2i) are antihyperglycemic medications that reduce cardiovascular disease (CVD) and improve chronic kidney disease prognosis in patients with diabetes mellitus. The specific impact of SGLT2i treatment on hypertensive individuals, however, remains to be established. This underscores the need for systematic efforts to profile the molecular landscape associated with SGLT2i administration.

METHODS

We conducted a detailed RNA-sequencing (RNA-Seq)-based exploration of transcriptomic changes in response to empagliflozin in eight different tissues (i.e., atrium, aorta, ventricle, white adipose, brown adipose, kidney, lung, and brain) from a male rat model of spontaneously hypertension. Corresponding computational analyses (i.e., clustering, differentially-expressed genes [DEG], and functional association) were performed to analyze these data. Blood pressure measurements, tissue staining studies and RT-qPCR were performed to validate our in silico findings.

RESULTS

We discovered that empagliflozin exerted potent transcriptomic effects on various tissues, most notably the kidney, white adipose, and lung in spontaneously hypertension rats (SHR). The functional enrichment of DEGs indicated that empagliflozin may regulate blood pressure, blood glucose and lipid homeostasis in SHR. Consistent with our RNA-Seq findings, immunohistochemistry and qPCR analyses revealed decreased renal expression of mitogen-activated protein kinase 10 (MAPK10) and decreased pulmonary expression of the proinflammatory factors Legumain and cathepsin S (CTSS) at 1 month of empagliflozin administration. Notably, immunofluorescence experiments showed increased expression of the AMP-activated protein kinases Prkaa1 and Prkaa2 in white adipose tissue of SHR following empagliflozin therapy. Furthermore, the transcriptomic signatures of the blood pressure-lowing effect by empagliflozin were experimentally validated in SHR.

CONCLUSIONS

This study provided an important resource of the effects of empagliflozin on various tissues of SHRs. We identified tissue-specific and tissue-enriched transcriptomic signatures, and uncovered the beneficial effects of empagliflozin on hypertension, weight gain and inflammatory response in validated experiments.

CD44 correlates with longevity and enhances basal ATF6 activity and ER stress resistance

The naked mole rat (NMR) is the longest-lived rodent, resistant to multiple age-related diseases including neurodegeneration. However, the mechanisms underlying the NMR's resistance to neurodegenerative diseases remain elusive. Here, we isolated oligodendrocyte progenitor cells (OPCs) from NMRs and compared their transcriptome with that of other mammals. Extracellular matrix (ECM) genes best distinguish OPCs of long- and short-lived species. Notably, expression levels of CD44, an ECM-binding protein that has been suggested to contribute to NMR longevity by mediating the effect of hyaluronan (HA), are not only high in OPCs of long-lived species but also positively correlate with longevity in multiple cell types/tissues. We found that CD44 localizes to the endoplasmic reticulum (ER) and enhances basal ATF6 activity. CD44 modifies proteome and membrane properties of the ER and enhances ER stress resistance in a manner dependent on unfolded protein response regulators without the requirement of HA. HA-independent role of CD44 in proteostasis regulation may contribute to mammalian longevity.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Polygenic genetic variation affecting antibody formation underlies hypertensive renal injury in the stroke-prone spontaneously hypertensive rat

During development of the spontaneously hypertensive rat (SHR), several distinct but closely related lines were generated. Most lines are resistant to hypertensive renal disease. However, the SHR-A3 line (stroke-prone SHR) experiences end-organ injury (EOI) and provides a model of injury susceptibility that can be used to uncover genetic causation. In the present study, we generated a congenic line in which three distinct disease loci in SHR-A3 are concurrently replaced with homologous loci from an injury-resistant SHR line (SHR-B2). Verification that all three loci were homozygously replaced in this triple congenic line [SHR-A3(Trip B2)] while the genetic background of SHR-A3 was fully retained was obtained by whole genome sequencing. Congenic genome substitution was without effect on systolic blood pressure [198.9 ± 3.34 mmHg, mean ± SE, SHR-A3(Trip B2) = 194.7 ± 2.55 mmHg]. Measures of renal injury (albuminuria, histological injury scores, and urinary biomarker levels) were reduced in SHR-A3(Trip B2) animals, even though only 4.5 Mbases of the 2.8 Gbases of the SHR-B2 genome (0.16% of the genome) was transferred into the congenic line. The gene content of the three congenic loci and the functional effects of gene polymorphism within suggest a role of immunoglobulin in EOI pathogenesis. To prove the role of antibodies in EOI in SHR-A3, we generated an SHR-A3 line in which expression from the immunoglobulin heavy chain gene was knocked out (SHR-A3-IGHKO). Animals in the SHR-A3-IGHKO line lack B cells and immunoglobulin, but the hypertensive phenotype is not affected. Renal injury, however, was reduced in this line, confirming a pathogenic role for immunoglobulin in hypertensive EOI in this model of heritable risk.NEW & NOTEWORTHY Here, we used a polygenic animal model of hypertensive renal disease to show that genetic variation affecting antibody formation underlies hypertensive renal disease. We proved the genetic thesis by generating an immunoglobulin knockout in the susceptible animal model.

Research Note: Identification of core promoter region of the polyunsaturated fatty acid synthesis-related gene family in chicken

Chicken is considered an ideal model species to study the synthesis of polyunsaturated fatty acids (PUFAs) due to its appropriate proportions of fatty acids and abundant content of PUFAs, suitable for human consumption. However, the molecular mechanisms regulating poultry PUFA synthesis remain unclear. Here, we systematically explored the transcriptional regulation activity of the gene family related to PUFA synthesis in chicken by carrying out the Dual-Luciferase Reporter Assay. We identified the core promoter regions of members of the chicken PUFA synthesis-related gene family, including ELOVL1, ELOVL2, ELOVL3, ELOVL4, ELOVL5, ELOVL6, ELOVL7, FADS1, FADS2, FADS6, SCD, and SCD5. Additionally, changes in relative fluorescence values of different truncated segments in the upstream regulatory region of these genes indicate the existence of regulatory regions. Furthermore, we predicted the transcription factors that bind to the identified core promoter regions of multiple genes, including Sp1, NF-1, C/EBPalpha, etc. These findings provide a basis for the molecular mechanisms regulating poultry PUFA synthesis and offer new scientific insight into the potential improvement of poultry meat quality in the future.

Amino Acid Catabolism: An Overlooked Area of Metabolism

Amino acids have been extensively studied in nutrition, mainly as key elements for maintaining optimal protein synthesis in the body as well as precursors of various nitrogen-containing compounds. However, it is now known that amino acid catabolism is an important element for the metabolic control of different biological processes, although it is still a developing field to have a deeper understanding of its biological implications. The mechanisms involved in the regulation of amino acid catabolism now include the contribution of the gut microbiota to amino acid oxidation and metabolite generation in the intestine, the molecular mechanisms of transcriptional control, and the participation of specific miRNAs involved in the regulation of amino acid degrading enzymes. In addition, molecules derived from amino acid catabolism play a role in metabolism as they are used in the epigenetic regulation of many genes. Thus, this review aims to examine the mechanisms of amino acid catabolism and to support the idea that this process is associated with the immune response, abnormalities during obesity, in particular insulin resistance, and the regulation of thermogenesis.

Neutrophil extracellular traps involved in the pathogenesis of IgA vasculitis: Confirmed in two IgAV rat models

BACKGROUND

Neutrophil extracellular traps (NETs) have been found to play a role in the development of autoimmune diseases. In the past two years, studies have demonstrated a significantly increase of NETs in skin tissues during the early stages of IgAV, indicating their involvement in disease activity among children with IgAV. However, the presence of NETs in IgAV animal models has not yet been reported. The objective of this study is to investigate whether NETs are involved in the pathogenesis of IgA vasculitis (IgAV) rats.

METHODS

Twenty-four SD rats were randomly divided into three groups: the ovalbumin group, the gliadin group, and the control group. The IgAV rat models were established administering Indian ink with ovalbumin (ovalbumin group) or gliadin (gliadin group) with Freund's complete adjuvant. The cell-free DNA (cf-DNA) was quantified by using dsDNA quantification kit, while the levels of Immunoglobulins, complement C3 and myeloperoxidase-DNA (MPO-DNA) in serum were tested using enzyme linked immunosorbent assay (ELISA). The IgA, complement C3 and NETs in tissues were detected through multiple immunofluorescences.

RESULTS

Both the ovalbumin group and gliadin group showed IgA and C3 deposition in various tissues, including the glomerular mesangial region, skin, and digestive tract, while the control group showed no such deposition. The levels of circulatory cf-DNA and MPO-DNA, which are components of NETs, were significantly elevated in both ovalbumin and gliadin groups compared with the control group. Furthermore, the presence of NETs were found in gastrointestinal and renal tissues of the ovalbumin and gliadin groups, but not in the control group.

CONCLUSIONS

IgAV model rat can be established through the combination of ovalbumin and gliadin with Indian ink and Freund's complete adjuvant. This study provides the first confirmation that NETs are involved in the pathogenesis of IgAV rat.

The role of Zn2+ in shaping intracellular Ca2+ dynamics in the heart

Increasing evidence suggests that Zn2+ acts as a second messenger capable of transducing extracellular stimuli into intracellular signaling events. The importance of Zn2+ as a signaling molecule in cardiovascular functioning is gaining traction. In the heart, Zn2+ plays important roles in excitation-contraction (EC) coupling, excitation-transcription coupling, and cardiac ventricular morphogenesis. Zn2+ homeostasis in cardiac tissue is tightly regulated through the action of a combination of transporters, buffers, and sensors. Zn2+ mishandling is a common feature of various cardiovascular diseases. However, the precise mechanisms controlling the intracellular distribution of Zn2+ and its variations during normal cardiac function and during pathological conditions are not fully understood. In this review, we consider the major pathways by which the concentration of intracellular Zn2+ is regulated in the heart, the role of Zn2+ in EC coupling, and discuss how Zn2+ dyshomeostasis resulting from altered expression levels and efficacy of Zn2+ regulatory proteins are key drivers in the progression of cardiac dysfunction.

Revisiting the expression of BDNF and its receptors in mammalian development

Brain-derived neurotrophic factor (BDNF) promotes the survival and functioning of neurons in the central nervous system and contributes to proper functioning of many non-neural tissues. Although the regulation and role of BDNF have been extensively studied, a rigorous analysis of the expression dynamics of BDNF and its receptors TrkB and p75NTR is lacking. Here, we have analyzed more than 3,600 samples from 18 published RNA sequencing datasets, and used over 17,000 samples from GTEx, and ~ 180 samples from BrainSpan database, to describe the expression of BDNF in the developing mammalian neural and non-neural tissues. We show evolutionarily conserved dynamics and expression patterns of BDNF mRNA and non-conserved alternative 5' exon usage. Finally, we also show increasing BDNF protein levels during murine brain development and BDNF protein expression in several non-neural tissues. In parallel, we describe the spatiotemporal expression pattern of BDNF receptors TrkB and p75NTR in both murines and humans. Collectively, our in-depth analysis of the expression of BDNF and its receptors gives insight into the regulation and signaling of BDNF in the whole organism throughout life.

TransOrGAN: An Artificial Intelligence Mapping of Rat Transcriptomic Profiles between Organs, Ages, and Sexes

Animal studies are required for the evaluation of candidate drugs to ensure patient and volunteer safety. Toxicogenomics is often applied in these studies to gain understanding of the underlying mechanisms of toxicity, which is usually focused on critical organs such as the liver or kidney in young male rats. There is a strong ethical reason to reduce, refine and replace animal use (the 3Rs), where the mapping of data between organs, sexes and ages could reduce the cost and time of drug development. Herein, we proposed a generative adversarial network (GAN)-based framework entitled TransOrGAN that allowed the molecular mapping of gene expression profiles in different rodent organ systems and across sex and age groups. We carried out a proof-of-concept study based on rat RNA-seq data from 288 samples in 9 different organs of both sexes and 4 developmental stages. First, we demonstrated that TransOrGAN could infer transcriptomic profiles between any 2 of the 9 organs studied, yielding an average cosine similarity of 0.984 between synthetic transcriptomic profiles and their corresponding real profiles. Second, we found that TransOrGAN could infer transcriptomic profiles observed in females from males, with an average cosine similarity of 0.984. Third, we found that TransOrGAN could infer transcriptomic profiles in juvenile, adult, and aged animals from adolescent animals with an average cosine similarity of 0.981, 0.983, and 0.989, respectively. Altogether, TransOrGAN is an innovative approach to infer transcriptomic profiles between ages, sexes, and organ systems, offering the opportunity to reduce animal usage and to provide an integrated assessment of toxicity in the whole organism irrespective of sex or age.

Transcriptome Profile of the Rat Adrenal Gland: Parenchymal and Interstitial Cells

The homeostasis of the adrenal gland plays a decisive role in its proper functioning, both in non-stressful conditions and under the influence of various types of stress. This consists of interactions between all types of cells that make up the organ, including parenchymal and interstitial cells. The amount of available information on this subject in the rat adrenal glands under non-stressful conditions is insufficient; the aim of the research was to determine the expression of marker genes for rat adrenal cells depending on their location. The material for the study consisted of adrenal glands taken from intact adult male rats that were separated into appropriate zones. Transcriptome analysis by means of Affymetrix^® Rat Gene 2.1 ST Array was used in the study, followed by real-time PCR validation. Expression analysis of interstitial cell marker genes revealed both the amount of expression of these genes and the zone in which they were expressed. The expression of marker genes for fibroblasts was particularly high in the cells of the ZG zone, while the highest expression of specific macrophage genes was observed in the adrenal medulla. The results of this study, especially with regard to interstitial cells, provide a so far undescribed model of marker gene expression of various cells, both in the cortex and medulla of the sexually mature rat adrenal gland. The interdependence between parenchymal and interstitial cells creates a specific microenvironment that is highly heterogeneous within the gland with respect to some of the interstitial cells. This phenomenon most likely depends on the interaction with the differentiated parenchymal cells of the cortex, as well as the medulla of the gland.

Rat microbial biogeography and age-dependent lactic acid bacteria in healthy lungs

The laboratory rat emerges as a useful tool for studying the interaction between the host and its microbiome. To advance principles relevant to the human microbiome, we systematically investigated and defined a multi-tissue full lifespan microbial biogeography for healthy Fischer 344 rats. Microbial community profiling data was extracted and integrated with host transcriptomic data from the Sequencing Quality Control (SEQC) consortium. Unsupervised machine learning, Spearman's correlation, taxonomic diversity, and abundance analyses were performed to determine and characterize the rat microbial biogeography and the identification of four inter-tissue microbial heterogeneity patterns (P1-P4). The 11 body habitats harbor a greater diversity of microbes than previously suspected. Lactic acid bacteria (LAB) abundances progressively declined in lungs from breastfeed newborn to adolescence/adult and was below detectable levels in elderly rats. LAB's presence and levels in lungs were further evaluated by PCR in the two validation datasets. The lung, testes, thymus, kidney, adrenal, and muscle niches were found to have age-dependent alterations in microbial abundance. P1 is dominated by lung samples. P2 contains the largest sample size and is enriched for environmental species. Liver and muscle samples were mostly classified into P3. Archaea species were exclusively enriched in P4. The 357 pattern-specific microbial signatures were positively correlated with host genes in cell migration and proliferation (P1), DNA damage repair and synaptic transmissions (P2), as well as DNA transcription and cell cycle in P3. Our study established a link between metabolic properties of LAB with lung microbiota maturation and development. Breastfeeding and environmental exposure influence microbiome composition and host health and longevity. The inferred rat microbial biogeography and pattern-specific microbial signatures would be useful for microbiome therapeutic approaches to human health and good quality of life.

The methylimidazolium ionic liquid M8OI is a substrate for OCT1 and p-glycoprotein-1 in rat

The methylimidazolium ionic liquid M8OI was recently found to be present in both the environment and man. In this study, M8OI disposition and toxicity were examined in an established rat progenitor-hepatocyte model. The progenitor B-13 cell was approx. 13 fold more sensitive to the toxic effects of M8OI than the hepatocyte B-13/H cell. However, this difference in sensitivity was not associated with a difference in metabolic capacities. M8OI toxicity was significantly decreased in a dose-dependent manner by co-addition of the OCT1 (SLC22A1) inhibitor clonidine, but not by OCT2 or OCT3 inhibitors in B-13 cells. M8OI toxicity was also dose-dependently increased by the co-addition of p-glycoprotein-1 (ABCB1B, multi drug resistant protein 1 (MDR1)) substrates/inhibitors. Excretion of B-13-loaded fluorophore Hoechst 33342 was also inhibited by the p-glycoproteins substrate cyclosporin A and by M8OI in a dose-dependent manner. Comparing levels of OCT and p-glycoprotein transcripts and proteins in B-13 and B-13/H cells suggest that the lower sensitivity to M8OI in B-13/H cells is predominantly associated with their higher expression of p-glycoprotein-1. These data together therefore suggest that a determinant in M8OI toxicity in rats is the expression and activity of the p-glycoprotein-1 transporter.

Gene Self-Expressive Networks as a Generalization-Aware Tool to Model Gene Regulatory Networks

Self-expressiveness is a mathematical property that aims at characterizing the relationship between instances in a dataset. This property has been applied widely and successfully in computer-vision tasks, time-series analysis, and to infer underlying network structures in domains including protein signaling interactions and social-networks activity. Nevertheless, despite its potential, self-expressiveness has not been explicitly used to infer gene networks. In this article, we present Generalizable Gene Self-Expressive Networks, a new, interpretable, and generalization-aware formalism to model gene networks, and we propose two methods: GXN•EN and GXN•OMP, based respectively on ElasticNet and OMP (Orthogonal Matching Pursuit), to infer and assess Generalizable Gene Self-Expressive Networks. We evaluate these methods on four Microarray datasets from the DREAM5 benchmark, using both internal and external metrics. The results obtained by both methods are comparable to those obtained by state-of-the-art tools, but are fast to train and exhibit high levels of sparsity, which make them easier to interpret. Moreover we applied these methods to three complex datasets containing RNA-seq informations from different mammalian tissues/cell-types. Lastly, we applied our methodology to compare a normal vs. a disease condition (Alzheimer), which allowed us to detect differential expression of genes’ sub-networks between these two biological conditions. Globally, the gene networks obtained exhibit a sparse and modular structure, with inner communities of genes presenting statistically significant over/under-expression on specific cell types, as well as significant enrichment for some anatomical GO terms, suggesting that such communities may also drive important functional roles.

Genomic regulation of Krüppel-like-factor family members by corticosteroid receptors in the rat brain

Hippocampal mineralocorticoid receptors (MRs) and glucocorticoid receptors (GRs) mediate glucocorticoid hormone (GC) action in the hippocampus. These receptors bind to glucocorticoid responsive elements (GREs) within target genes, eliciting transcriptional effects in response to stress and circadian variation. Until recently, little was known about the genome-wide targets of hippocampal MRs and GRs under physiological conditions. Following on from our genome-wide MR and GR ChIP-seq and Ribo-Zero RNA-seq studies on rat hippocampus, we investigated the Krüppel-like factors (KLFs) as targets of MRs and GRs throughout the brain under circadian variation and after acute stress. In particular, Klf2, Klf9 and Klf15 are known to be stress and/or GC responsive and play a role in neurobiological processes including synaptic plasticity and neuronal differentiation. We found increased binding of MR and GR to GREs within Klf2, Klf9 and Klf15 in the hippocampus, amygdala, prefrontal cortex, and neocortex after acute stress and resulting from circadian variation, which was accompanied by upregulation of corresponding hnRNA and mRNA levels. Adrenalectomy abolished transcriptional upregulation of specific Klf genes. These results show that MRs and GRs regulate Klf gene expression throughout the brain following exposure to acute stress or in response to circadian variation, likely alongside other transcription factors.

Developmental programming: adverse sexually dimorphic transcriptional programming of gestational testosterone excess in cardiac left ventricle of fetal sheep

Adverse in-utero insults during fetal life alters offspring's developmental trajectory, including that of the cardiovascular system. Gestational hyperandrogenism is once such adverse in-utero insult. Gestational testosterone (T)-treatment, an environment of gestational hyperandrogenism, manifests as hypertension and pathological left ventricular (LV) remodeling in adult ovine offspring. Furthermore, sexual dimorphism is noted in cardiomyocyte number and morphology in fetal life and at birth. This study investigated transcriptional changes and potential biomarkers of prenatal T excess-induced adverse cardiac programming. Genome-wide coding and non-coding (nc) RNA expression were compared between prenatal T-treated (T propionate 100 mg intramuscular twice weekly from days 30 to 90 of gestation; Term: 147 days) and control ovine LV at day 90 fetus in both sexes. Prenatal T induced differential expression of mRNAs in the LV of female (2 down, 5 up) and male (3 down, 1 up) (FDR < 0.05, absolute log2 fold change > 0.5); pathways analysis demonstrated 205 pathways unique to the female, 382 unique to the male and 23 common pathways. In the male, analysis of ncRNA showed differential regulation of 15 lncRNAs (14 down, 1 up) and 27 snoRNAs (26 down and 1 up). These findings suggest sexual dimorphic modulation of cardiac coding and ncRNA with gestational T excess.

Characterization of the MMP9 Gene and Its Association with Cryptocaryon irritans Resistance Traits in Trachinotus ovatus (Linnaeus, 1758)

The MMPs are endogenous proteolytic enzymes that require zinc and calcium as cofactors. MMP9 is one of the most complex matrix metalloproteinases in the gelatinase family and has many biological functions. In mammals, mmp9 is thought to be closely associated with cancer. However, studies in fish have rarely been reported. In this study, to understand the expression pattern of the ToMMP9 gene and its association with the resistance of Trachinotus ovatus to Cryptocaryon irritans, the sequence of the MMP9 gene was obtained from the genome database. The expression profiles were measured by qRT-PCR, the SNPs were screened by direct sequencing, and genotyping was performed. The ToMMP9 gene contained a 2058 bp ORF encoding a putative amino acid sequence of 685 residues. The homology of the ToMMP9 in teleosts was more than 85%, and the genome structure of ToMMP9 was conserved in chordates. The ToMMP9 gene was expressed in different tissues of healthy individuals and was highly expressed in the fin, the gill, the liver and the skin tissues. The ToMMP9 expression in the skin of the infected site and its adjacent sites increased significantly after C. irritans infection. Two SNPs were identified in the ToMMP9 gene, and the SNP (+400A/G) located in the first intron was found to be significantly associated with the susceptibility/resistance to C. irritans. These findings suggest that ToMMP9 may play an important role in the immune response of T. ovatus against C. irritans.

Expression profiles of lncRNAs and their possible regulatory role in monocrotaline-induced HSOS in rats

Aims: Long non-coding RNAs (lncRNAs) contribute to the regulation of vital physiological processes and play a role in the pathogenesis of many diseases. Monocrotaline (MCT) can cause large-scale outbreaks of toxic liver disease in humans and animals in the form of hepatic sinusoidal obstruction syndrome (HSOS). Although many experiments have been carried out to explain the pathogenesis of Monocrotaline-induced hepatic sinusoidal obstruction syndrome and to develop treatments for it, no studies have examined the role of Long non-coding RNAs in this condition. This study aimed to investigate the Long non-coding RNAs-mRNA regulation network in Monocrotaline-induced hepatic sinusoidal obstruction syndrome in rats. Main methods: We established a model for MCT-induced hepatic sinusoidal obstruction syndrome, and then carried out microarray for liver tissues of SD rats in a model of early hepatic sinusoidal obstruction syndrome (12 h Monocrotaline treatment vs. control group) to investigate the differentially expressed Long non-coding RNAs and mRNAs in early hepatotoxicity. This was followed by RT-PCR analysis of selected Long non-coding RNAs, which were markedly altered. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome analyses were also conducted. Key findings: 176 Long non-coding RNAs (63 downregulated and 113 upregulated) and 4,221 mRNAs (2,385 downregulated and 1836 upregulated) were differentially expressed in the Monocrotaline-treated group compared to the control group. The biological processes identified in GO enrichment analysis as playing a role in hepatotoxicity were positive regulation of guanosine triphosphate phosphohydrolase, liver development, and the oxidation-reduction process. Pathway analysis revealed that the metabolism pathways, gap junction, and ribosome biogenesis in eukaryotes were closely related to Monocrotaline-induced hepatotoxicity. According to these analyses, LOC102552718 might play an essential role in hepatotoxicity mechanisms by regulating the expression of inositol 1,4,5-trisphosphate receptor-1 (Itpr-1). Significance: This study provides a basis for further research on the molecular mechanisms underlying Monocrotaline-induced hepatotoxicity and its treatment, especially in the early stage, when successful treatment is critical before irreversible liver damage occurs.

Preclinical species gene expression database: Development and meta-analysis

The evaluation of toxicity in preclinical species is important for identifying potential safety liabilities of experimental medicines. Toxicology studies provide translational insight into potential adverse clinical findings, but data interpretation may be limited due to our understanding of cross-species biological differences. With the recent technological advances in sequencing and analyzing omics data, gene expression data can be used to predict cross species biological differences and improve experimental design and toxicology data interpretation. However, interpreting the translational significance of toxicogenomics analyses can pose a challenge due to the lack of comprehensive preclinical gene expression datasets. In this work, we performed RNA-sequencing across four preclinical species/strains widely used for safety assessment (CD1 mouse, Sprague Dawley rat, Beagle dog, and Cynomolgus monkey) in ∼50 relevant tissues/organs to establish a comprehensive preclinical gene expression body atlas for both males and females. In addition, we performed a meta-analysis across the large dataset to highlight species and tissue differences that may be relevant for drug safety analyses. Further, we made these databases available to the scientific community. This multi-species, tissue-, and sex-specific transcriptomic database should serve as a valuable resource to enable informed safety decision-making not only during drug development, but also in a variety of disciplines that use these preclinical species.

Age-dependent expression of ion channel genes in rat

Ion channels regulate a large number of cellular functions and their functional role in many diseases makes them potential therapeutic targets. Given their diverse distribution across multiple organs, the roles of ion channels, particularly in age-associated transcriptomic changes in specific organs, are yet to be fully revealed. Using RNA-seq data, we investigated the rat transcriptomic profiles of ion channel genes across 11 organs/tissues and 4 developmental stages in both sexes of Fischer 344 rats and identify tissue-specific and age-dependent changes in ion channel gene expression. Organ-enriched ion channel genes were identified. In particular, the brain showed higher tissue-specificity of ion channel genes, including Gabrd, Gabra6, Gabrg2, Grin2a, and Grin2b. Notably, age-dependent changes in ion channel gene expression were prominently observed in the thymus, including in Aqp1, Clcn4, Hvcn1, Itpr1, Kcng2, Kcnj11, Kcnn3, and Trpm2. Our comprehensive study of ion channel gene expression will serve as a primary resource for biological studies of aging-related diseases caused by abnormal ion channel functions.

Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.

The expression profile of 79 genes from 107 viruses revealed new insights into disease susceptibility in rats, mice, and muskrats

The coronavirus disease 2019 (COVID-19), which is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has rapidly spread over the world, resulting in a global severe pneumonia pandemic. Both the cell receptor angiotensin-converting enzyme 2 (ACE2) and the breakdown of S protein by transmembrane serine protease 2 (TMPRSS2) are required by SARS-CoV-2 to enter the host cells. Similarly, the expression level of viral receptor genes in various organs determines the likelihood of viral infection. Several animal species have been found to be infected by the SARS-CoV-2, such as minks, posing an enormous threat to humans. Because the mice and rats were closely related to human and the fact that rats and mice have a risk of infection by SARS-CoV-2 with specific variants, we investigated the expression patterns of 79 receptor genes from 107 viruses, including SARS-CoV-2, in 14 organs of the rat and mouse, and 5 organs of the muskrat, to find the most likely host organs to become infected with certain viruses. The findings of this study are anticipated to aid in prevention of zoonotic infections spread by rats, mice, muskrats, and other rodents.

Sexually concordant and dimorphic transcriptional responses to maternal trichloroethylene and/or N-acetyl cysteine exposure in Wistar rat placental tissue

Numerous Superfund sites are contaminated with the volatile organic chemical trichloroethylene (TCE). In women, exposure to TCE in pregnancy is associated with reduced birth weight. Our previous study reported that TCE exposure in pregnant rats decreased fetal weight and elevated oxidative stress biomarkers in placentae, suggesting placental injury as a potential mechanism of TCE-induced adverse birth outcomes. In this study, we investigated if co-exposure with the antioxidant N-acetylcysteine (NAC) attenuates TCE exposure effects on RNA expression. Timed-pregnant Wistar rats were exposed orally to 480 mg TCE/kg/day on gestation days 6-16. Exposure of 200 mg NAC/kg/day alone or as a pre/co-exposure with TCE occurred on gestation days 5-16 to stimulate antioxidant genes prior to TCE exposure. Tissue was collected on gestation day 16. In male and female placentae, we evaluated TCE- and/or NAC-induced changes to gene expression and pathway enrichment analyses using false discovery rate (FDR) and fold-change criteria. In female placentae, exposure to TCE caused significant differential expression 129 genes while the TCE+NAC altered 125 genes, compared with controls (FDR< 0.05 + fold-change >1). In contrast, in male placentae TCE exposure differentially expressed 9 genes and TCE+NAC differentially expressed 35 genes, compared with controls (FDR< 0.05 + fold-change >1). NAC alone did not significantly alter gene expression in either sex. Differentially expressed genes observed with TCE exposure were enriched in mitochondrial biogenesis and oxidative phosphorylation pathways in females whereas immune system pathways and endoplasmic reticulum stress pathways were differentially expressed in both sexes (FDR<0.05). TCE treatment was differentially enriched for genes regulated by the transcription factors ATF6 (both sexes) and ATF4 (males only), indicating a cellular condition triggered by misfolded proteins during endoplasmic reticulum stress. This study demonstrates novel genes and pathways involved in TCE-induced placental injury and showed antioxidant co-treatment largely did not attenuate TCE exposure effects.

Gpx4, Selenov, and Txnrd3 Are Three Most Testis-Abundant Selenogenes Resistant to Dietary Selenium Concentrations and Actively Expressed During Reproductive Ages in Rats

Almost all selenogenes are expressed in the testis, and those have the highest and constant expressions will be the primary candidates for functional analysis of selenium (Se) in male reproduction. This study aimed to profile the mRNA expressions of the testis-abundant selenogenes of rat models in responses to growth and dietary Se concentrations. Forty-eight weaning SD male rats were fed Se deficient basal diet (BD) for 5 weeks and then randomly grouped (n = 12/group) for being fed BD or BD plus 0.25, 3, or 5 mg Se/kg for 4 more weeks before sacrifice. Abundances of selenogenomic mRNAs in the liver and testis were determined with relative qPCR and those of the testis-abundant selenogenes in 13 kinds of tissues were assayed with a molecular beacon-based qPCR. Spatiotemporal expressions of rat selenogenome were also analyzed with the RNA-Seq transcriptomic data published by NCBI. mRNA abundances of glutathione peroxidase 4 (Gpx4), nuclear Gpx4 (nGpx4), selenoprotein V (Selenov), and thioredoxin reductase 3 (Txnrd3) in the testis were significantly higher than that in any other tissues (P < 0.05). Moreover, testicular mRNA abundances of Gpx4, Selenov, and Txnrd3 were not affected by levels of dietary Se supplementation (P > 0.05), and much higher at 6-21 weeks old than at 2 and 104 weeks old (P < 0.05). The result showed that Gpx4, Selenov, and Txnrd3 were most highly expressed in the testis of rats especially at reproductive ages and resistant to the impact of dietary Se levels, which suggested their specific importance in male reproduction.

Expression of alternative transcription factor 4 mRNAs and protein isoforms in the developing and adult rodent and human tissues

Transcription factor 4 (TCF4) belongs to the class I basic helix-loop-helix family of transcription factors (also known as E-proteins) and is vital for the development of the nervous system. Aberrations in the TCF4 gene are associated with several neurocognitive disorders such as schizophrenia, intellectual disability, post-traumatic stress disorder, depression, and Pitt-Hopkins Syndrome, a rare but severe autism spectrum disorder. Expression of the human TCF4 gene can produce at least 18 N-terminally distinct protein isoforms, which activate transcription with different activities and thus may vary in their function during development. We used long-read RNA-sequencing and western blot analysis combined with the analysis of publicly available short-read RNA-sequencing data to describe both the mRNA and protein expression of the many distinct TCF4 isoforms in rodent and human neural and nonneural tissues. We show that TCF4 mRNA and protein expression is much higher in the rodent brain compared to nonneural tissues. TCF4 protein expression is highest in the rodent cerebral cortex and hippocampus, where expression peaks around birth, and in the rodent cerebellum, where expression peaks about a week after birth. In human, highest TCF4 expression levels were seen in the developing brain, although some nonneural tissues displayed comparable expression levels to adult brain. In addition, we show for the first time that out of the many possible TCF4 isoforms, the main TCF4 isoforms expressed in the rodent and human brain and other tissues are TCF4-B, -C, -D, -A, and-I. Taken together, our isoform specific analysis of TCF4 expression in different tissues could be used for the generation of gene therapy applications for patients with TCF4-associated diseases.

DPP3: From biomarker to therapeutic target of cardiovascular diseases

Cardiovascular disease is the leading cause of death globally among non-communicable diseases, which imposes a serious socioeconomic burden on patients and the healthcare system. Therefore, finding new strategies for preventing and treating cardiovascular diseases is of great significance in reducing the number of deaths and disabilities worldwide. Dipeptidyl peptidase 3 (DPP3) is the first zinc-dependent peptidase found among DPPs, mainly distributes within the cytoplasm. With the unique HEXXGH catalytic sequence, it is associated with the degradation of oligopeptides with 4 to 10 amino acids residues. Accumulating evidences have demonstrated that DPP3 plays a significant role in almost all cellular activities and pathophysiological mechanisms. Regarding the role of DPP3 in cardiovascular diseases, it is currently mainly used as a biomarker for poor prognosis in patients with cardiovascular diseases, suggesting that the level of DPP3 concentration in plasma is closely linked to the mortality of diseases such as cardiogenic shock and heart failure. Interestingly, it has been reported recently that DPP3 regulates blood pressure by interacting with the renin-angiotensin system. In addition, DPP3 also participates in the processes of pain signaling, inflammation, and oxidative stress. But the exact mechanism by which DPP3 affects cardiovascular function is not clear. Hence, this review summarizes the recent advances in the structure and catalytic activity of DPP3 and its extensive biological functions, especially its role as a therapeutic target in cardiovascular diseases. It will provide a theoretical basis for exploring the potential value of DPP3 as a therapeutic target for cardiovascular diseases.