Update of bioRxiv. 2023 Apr 08;:

Update of bioRxiv. 2023 Jun 25;:

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.

Chronic spinal cord injury (SCI) results in an increased predisposition to various metabolic problems that can be exacerbated by consuming a diet rich in calories and saturated fat. In addition, gastrointestinal symptoms have been reported after SCI, including intestinal dysbiosis of the gut microbiome. The effects of both diet and SCI on the gut microbiome of adult male Long Evans rats euthanized 16 wk after injury were investigated. The rats were either thoracic spinal contused or received sham procedures. After 12 wk of either a low-fat or high-fat diet, cecal contents were analyzed, revealing significant microbial changes to every taxonomic level below the kingdom level. Shannon α diversity analyses demonstrated a significant difference in diversity between the groups based on the surgical condition of the rats. SCI produced a unique signature of changes in commensal bacteria that were significantly different than Sham. Specific changes in commensal bacteria as a result of diet manipulation had high fidelity with reports in the literature, such as Clostridia, Thiohalorhabdales, and Pseudomonadales. In addition, novel changes in commensal bacteria were identified that are unique dietary influences on SCI. Linear regression analysis on body fat and lean mass showed that a consequence of chronic SCI produces uncoupled associations between some commensal bacteria and body composition. In conclusion, despite tightly controlling the protein content and varying the carbohydrate and fat contents, Sham and SCI rats respond uniquely to diet. These data provide potential direction for therapeutic modulation of the microbiome to improve health and wellness following SCI.

Military Deployment to Southwest Asia and Afghanistan and exposure to toxic airborne particulates have been associated with an increased risk of developing respiratory disease, collectively termed deployment-related respiratory diseases (DRRDs). Our knowledge about how particulates mediate respiratory disease is limited, precluding the appropriate recognition or management. Central to this limitation is the lack of understanding of how exposures translate into dysregulated cell identity with dysregulated transcriptional programs. The small airway epithelium is involved in both the pathobiology of DRRD and fine particulate matter deposition. To characterize small airway epithelial cell epigenetic and transcriptional responses to Afghan desert particulate matter (APM) and investigate the functional interactions of transcription factors that mediate these responses, we applied two genomics assays, the assay for transposase accessible chromatin with sequencing (ATAC-seq) and Precision Run-on sequencing (PRO-seq). We identified activity changes in a series of transcriptional pathways as candidate regulators of susceptibility to subsequent insults, including signal-dependent pathways, such as loss of cytochrome P450 or P53/P63, and lineage-determining transcription factors, such as GRHL2 loss or TEAD3 activation. We further demonstrated that TEAD3 activation was unique to APM exposure despite similar inflammatory responses when compared with wood smoke particle exposure and that P53/P63 program loss was uniquely positioned at the intersection of signal-dependent and lineage-determining transcriptional programs. Our results establish the utility of an integrated genomics approach in characterizing responses to exposures and identifying genomic targets for the advanced investigation of the pathogenesis of DRRD.

Comment in Physiol Genomics. 2023 Feb 1;55(2):75-78.

miRNA155 (miR155) has emerged as an important regulator of breast cancer (BrCa) development. Studies have consistently noted an increase in miR155 levels in serum and/or tissues in patients with BrCa. However, what is less clear is whether this increase in miR155 is a reflection of oncogenic or tumor suppressive properties. To study the effects of miR155 in a transgenic model of BrCA, we developed an MMTV-PyMT mouse deficient in miR155 (miR155-/- PyMT). miR155-/- mice (n = 11) exhibited reduced tumor number and volume palpations at ∼14-18 wk of age compared with miR155 sufficient littermates (n = 12). At 19 wk, mammary glands were excised from tumors for RT-PCR, and tumors were counted, measured, and weighed. miR155-/- PyMT mice exhibited reduced tumor volume, number, and weight, which was confirmed by histopathological analysis. There was an increase in apoptosis with miR155 deficiency and a decrease in proliferation. As expected, miR155 deficiency resulted in upregulated gene expression of suppressor of cytokine signaling 1 (Socs1)-its direct target. There was a reduction in gene expression of macrophage markers (CD68, Adgre1, Itgax, Mrc1) with miR-155-/- and this was confirmed with immunofluorescence staining for F4/80. miR155-/- increased expression of M1 macrophage marker Nos2 and reduced expression of M2 macrophage markers IL-10, IL-4, Arg1, and MMP9. Overall, miR155 deficiency reduced BrCA and improved the tumor microenvironment through the reduction of genes associated with protumorigenic processes. However, given the inconsistencies in the literature, additional studies are needed before any attempts are made to harness miR155 as a potential oncogenic or tumor suppressive miRNA.NEW & NOTEWORTHY To examine the effects of miR155 in a transgenic model of breast cancer, we developed an MMTV-PyMT mouse-deficient in miR155. We demonstrate that global loss of miR155 resulted in blunted tumor growth through modulating the tumor microenvironment. Specifically, miR155-deficient mice had smaller and less invasive tumors, an increase in apoptosis and a decrease in proliferation, a reduction in tumor-associated macrophages, and the expression of genes associated with protumoral processes.

Blood pressure (BP) is influenced by genetic variation and sodium intake with sex-specific differences; however, studies to identify renal molecular mechanisms underlying the influence of sodium intake on BP in nonhuman primates (NHP) have focused on males. To address the gap in our understanding of molecular mechanisms regulating BP in female primates, we studied sodium-naïve female baboons (n = 7) fed a high-sodium (HS) diet for 6 wk. We hypothesized that in female baboons variation in renal transcriptional networks correlates with variation in BP response to a high-sodium diet. BP was continuously measured for 64-h periods throughout the study by implantable telemetry devices. Sodium intake, blood samples for clinical chemistries, and ultrasound-guided kidney biopsies were collected before and after the HS diet for RNA-Seq and bioinformatic analyses. We found that on the LS diet but not the HS diet, sodium intake and serum 17 β-estradiol concentration correlated with BP. Furthermore, kidney transcriptomes differed by diet-unbiased weighted gene coexpression network analysis revealed modules of genes correlated with BP on the HS diet but not the LS diet. Our results showed variation in BP on the HS diet correlated with variation in novel kidney gene networks regulated by ESR1 and MYC; i.e., these regulators have not been associated with BP regulation in male humans or rodents. Validation of the mechanisms underlying regulation of BP-associated gene networks in female NHP will inform better therapies toward greater precision medicine for women.

The vast majority of studies focusing on the effects of endurance exercise on hematological parameters and leukocyte gene expression were performed in adult men, so our aim was to investigate these changes in young females. Four young (age 15.3 ± 1.3 yr) elite female athletes completed an exercise session, in which they accomplished the cycling and running disciplines of a junior triathlon race. Blood samples were taken immediately before the exercise, right after the exercise, and then 1, 2, and 7 days later. Analysis of cell counts and routine biochemical parameters were complemented by RNA sequencing (RNA-seq) to whole blood samples. The applied exercise load did not trigger remarkable changes in either cardiovascular or biochemical parameters; however, it caused a significant increase in the percentage of neutrophils and a significant reduction in the ratio of lymphocytes immediately after exercise. Furthermore, endurance exercise induced a characteristic gene expression pattern change in the blood transcriptome. Gene set enrichment analysis (GSEA) using the Reactome database revealed that the expression of genes involved in immune processes and neutrophil granulocyte activation was upregulated, whereas the expression of genes important in translation and rRNA metabolism was downregulated. Comparison of a set of immune cell gene signatures (ImSig) and our transcriptomic data identified 15 overlapping genes related to T-cell functions and involved in podosome formation and adhesion to the vessel wall. Our results suggest that RNA-seq to whole blood together with ImSig analysis are useful tools for the investigation of systemic responses to endurance exercise.

Human-generated negative impacts on aquatic environments are rising. Despite wild fish playing a key role in aquatic ecologies and comprising a major global food source, physiological consequences of these impacts on them are poorly understood. Here we address the issue through the lens of interrelationship between wild fish and their gut microbiota, hypothesizing that fish microbiota are reporters of the aquatic environs. Two geographically separate teleost wild-fish species were studied (Lake Erie, Ohio, and Caribbean Sea, US Virgin Islands). At each geolocation, fresh fecal samples were collected from fish in areas of presence or absence of known aquatic compromise. Gut microbiota was assessed via microbial 16S-rRNA gene sequencing and represents the first complete report for both fish species. Despite marked differences in geography, climate, water type, fish species, habitat, diet, and gut microbial compositions, the pattern of shifts in microbiota shared by both fish species was nearly identical due to aquatic compromise. Next, these data were subjected to machine learning (ML) to examine reliability of using the fish-gut microbiota as an ecomarker for anthropogenic aquatic impacts. Independent of geolocation, ML predicted aquatic compromise with remarkable accuracy (>90%). Overall, this study represents the first multispecies stress-related comparison of its kind and demonstrates the potential of artificial intelligence via ML as a tool for biomonitoring and detecting compromised aquatic conditions.

Atherosclerosis in the carotid artery is a major cause of ischemic stroke and has a strong genetic component. The aim of this study was to identify genetic factors contributing to carotid atherosclerosis. One hundred fifty-four female F2 mice were generated from an intercross between LP/J and BALB/cJ Apoe-null (Apoe-/-) mice and fed 12 wk of Western diet. Atherosclerotic lesions, body weight, and coat color were measured and genotyping was performed using miniMUGA genotyping arrays. A significant quantitative trait locus (QTL) on chromosome (Chr) 7, named Cath20, and five suggestive QTL on Chr 6, 12, 13, 15, and X were identified for carotid lesions. Three significant QTL, Bwfq2, Bw1n, Bwtq6, on Chr 2, 7, and 15 were identified for body weight. Two significant QTL, Chop2 and Albc2, on Chr 4 and 7 were identified for coat color, with Tyr, encoding tyrosinase, being the causal gene of Albc2. Cath20 overlapped with or was close to QTL Bw1n and Albc2 on Chr7. Carotid lesion sizes were significantly correlated with body weight and graded coat color in F2 mice. Cath20 on Chr7 disappeared after adjustment for coat color but remained after adjustment for body weight. Tyr was abundantly expressed in atherosclerotic lesions. These results demonstrate genetic connections of carotid atherosclerosis with body weight and coat color in hyperlipidemic mice and suggest a potential role for Tyr in carotid atherosclerosis.

Increased levels of oxidative stress have been found with heart failure. Whether failing hearts express antioxidant and detoxification enzymes have not been addressed systematically. Nrf2 gene encodes a transcription factor that regulates the expression of antioxidant and detoxification genes. Using RNA-Seq data set from explanted hearts of 37 patients with dilated cardiomyopathy (DCM), 13 patients with ischemic cardiomyopathy (ICM), and 14 nonfailure (NF) donors as a control, we addressed whether failing hearts change the expression of Nrf2, its negative regulator Keap1, and antioxidant or detoxification genes. Significant increases in the ratio of Nrf2 to Keap1 were found to associate with DCM or ICM. Antioxidant genes showed decreased expression in both types of heart failure, including NQO1, SOD1, GPX3, GPX4, GSR, PRDX1, and TXNRD1. Detoxification enzymes, GCLM and EPHX1, also showed decreased expression, whereas the CYP1B1 transcript was elevated in both DCM and ICM. The genes encoding metal-binding protein ferritin were decreased, whereas five out of 12 metallothionein genes showed elevated expression. Our finding on Nrf2 gene expression has been validated by meta-analysis of seven independent data sets of microarray or RNA-Seq for differential gene expression in DCM and ICM from NF controls. In conclusion, minor elevation of Nrf2 gene expression is not coupled to increases in antioxidant and detoxification genes, supporting an impairment of Nrf2 signaling in patients with heart failure. Decreases in multiple antioxidant and detoxification genes are consistent with the observed increases of oxidative stress in failing hearts.

Transcriptomic analysis in metabolically active tissues allows a systems genetics approach to identify causal genes and networks involved in metabolic disease. Outbred heterogeneous stock (HS) rats are used for genetic mapping of complex traits, but to-date, a systems genetics analysis of metabolic tissues has not been done. We investigated whether adiposity-associated genes and gene coexpression networks in outbred heterogeneous stock (HS) rats overlap those found in humans. We analyzed RNAseq data from adipose tissue of 415 male HS rats, correlated these transcripts with body weight (BW) and compared transcriptome signatures to two human cohorts: the "African American Genetics of Metabolism and Expression" and "Metabolic Syndrome in Men." We used weighted gene coexpression network analysis to identify adiposity-associated gene networks and mediation analysis to identify genes under genetic control whose expression drives adiposity. We identified 554 orthologous "consensus genes" whose expression correlates with BW in the rat and with body mass index (BMI) in both human cohorts. Consensus genes fell within eight coexpressed networks and were enriched for genes involved in immune system function, cell growth, extracellular matrix organization, and lipid metabolic processes. We identified 19 consensus genes for which genetic variation may influence BW via their expression, including those involved in lipolysis (e.g., Hcar1), inflammation (e.g., Rgs1), adipogenesis (e.g., Tmem120b), or no previously known role in obesity (e.g., St14 and Ms4a6a). Strong concordance between HS rat and human BW/BMI associated transcripts demonstrates translational utility of the rat model, while identification of novel genes expands our knowledge of the genetics underlying obesity.

The brain renin-angiotensin system (RAS) is implicated in control of blood pressure (BP), fluid intake, and energy expenditure (EE). Angiotensin II (ANG II) within the arcuate nucleus of the hypothalamus contributes to control of resting metabolic rate (RMR) and thereby EE through its actions on Agouti-related peptide (AgRP) neurons, which also contribute to EE control by leptin. First, we determined that although leptin stimulates EE in control littermates, mice with transgenic activation of the brain RAS (sRA) exhibit increased EE and leptin has no additive effect to exaggerate EE in these mice. These findings led us to hypothesize that leptin and ANG II in the brain stimulate EE through a shared mechanism. Because AgRP signaling to the melanocortin MC4R receptor contributes to the metabolic effects of leptin, we performed a series of studies examining RMR, fluid intake, and BP responses to ANG II in mice rendered deficient for expression of MC4R via a transcriptional block (Mc4r-TB). These mice were resistant to stimulation of RMR in response to activation of the endogenous brain RAS via chronic deoxycorticosterone acetate (DOCA)-salt treatment, whereas fluid and electrolyte effects remained intact. These mice were also resistant to stimulation of RMR via acute intracerebroventricular (ICV) injection of ANG II, whereas BP responses to ICV ANG II remained intact. Collectively, these data demonstrate that the effects of ANG II within the brain to control RMR and EE are dependent on MC4R signaling, whereas fluid homeostasis and BP responses are independent of MC4R signaling.

Rat genomic tools have been slower to emerge than for those of humans and mice and have remained less thorough and comprehensive. The arrival of a new and improved rat reference genome, mRatBN7.2, in late 2020 is a welcome event. This assembly, like predecessor rat reference assemblies, is derived from an inbred Brown Norway rat. In this "user" survey we hope to provide other users of this assembly some insight into its characteristics and some assessment of its improvements as well as a few caveats that arise from the unique aspects of this assembly. mRatBN7.2 was generated by the Wellcome Sanger Institute as part of the large Vertebrate Genomes Project. This rat assembly has now joined human, mouse, chicken, and zebrafish in the National Center for Biotechnology Information (NCBI)'s Genome Reference Consortium, which provides ongoing curation of the assembly. Here we examine the technical procedures by which the assembly was created and assess how this assembly constitutes an improvement over its predecessor. We also indicate the technical limitations affecting the assembly, providing illustrations of how these limitations arise and the impact that results for this reference assembly.

Hypertension (HTN) is a complex disease influenced by heritable genetic elements and environmental interactions. Dietary salt is among the most influential modifiable factors contributing to increased blood pressure (BP). It is well established that men and women develop BP impairment in different patterns and a recent emphasis has been placed on identifying mechanisms leading to the differences observed between the sexes in HTN development. The current work reported here builds on an extensive genetic mapping experiment that sought to identify genetic determinants of salt-sensitive (SS) HTN using the Dahl SS rat. BTG antiproliferation factor 2 (Btg2) was previously identified by our group as a candidate gene contributing to SS HTN in female rats. In the current study, Btg2 was mutated using transcription activator-like effector nuclease (TALEN)-targeted gene disruption on the SSBN congenic rat background. The Btg2 mutated rats exhibited impaired BP and proteinuria responses to a high-salt diet compared with wild-type rats. Differences in body weight, mutant pup viability, skeletal morphology, and adult nephron density suggest a potential role for Btg2 in developmental signaling pathways. Subsequent cell cycle gene expression assessment provides several additional signaling pathways that Btg2 may function through during salt handling in the kidney. The expression analysis also identified several potential upstream targets that can be explored to further isolate therapeutic approaches for SS HTN.

Limited reports exist regarding adeno-associated virus (AAV) biodistribution in swine. This study assessed biodistribution following antegrade intracoronary and intravenous delivery of two self-complementary serotype 9 AAV (AAV9sc) biologics designed to target signaling in the cardiomyocyte considered important for the development of heart failure. Under the control of a cardiomyocyte-specific promoter, AAV9sc.shmAKAP and AAV9sc.RBD express a small hairpin RNA for the perinuclear scaffold protein muscle A-kinase anchoring protein β (mAKAPβ) and an anchoring disruptor peptide for p90 ribosomal S6 kinase type 3 (RSK3), respectively. Quantitative PCR was used to assess viral genome (vg) delivery and transcript expression in Ossabaw and Yorkshire swine tissues. Myocardial viral delivery was 2-5 × 105 vg/µg genomic DNA (gDNA) for both infusion techniques at a dose ∼1013 vg/kg body wt, demonstrating delivery of ∼1-3 viral particles per cardiac diploid genome. Myocardial RNA levels for each expressed transgene were generally proportional to dose and genomic delivery, and comparable with levels for moderately expressed endogenous genes. Despite significant AAV9sc delivery to other tissues, including the liver, neither biologic induced toxic effects as assessed using functional, structural, and circulating cardiac and systemic markers. These results indicate successful targeted delivery of cardiomyocyte-selective viral vectors in swine without negative side effects, an important step in establishing efficacy in a preclinical experimental setting.

Current knowledge of the link between microbiota and hypertension is limited to the gut. Besides the gut, oral cavity and skin are other locations where sodium chloride (NaCl) is in direct contact with microbiota. Although oral nitrate-reducing bacteria generate nitric oxide, which leads to vasodilation and lowering of blood pressure (BP), the skin excretes sodium via sweat glands and is an important site for sodium and BP homeostasis. However, knowledge on the contributions of oral and skin microbiota to BP regulation, is limited. Therefore, the current study was conducted to compare the tripartite relationship between site, sex, and genetic effects on the composition of oral, skin, and gut microbiota impacting hypertension. Microbiota were profiled from the oral cavity, skin, and feces of both male and female hypertensive Dahl salt-sensitive (S) and congenic rats with genomic substitutions on rat chromosomes (RNO) 1, 5, 9, and 10, demonstrating disparate BP effects. Sex-specific differences in β-diversity were observed only in skin microbiota. The most abundant taxa of the oral and skin microbiota were Actinobacteria and Cyanobacteria, respectively. Oral Actinobacteria were inversely associated with BP. Although the abundance of oral Actinobacteria was upregulated by the BP locus on RNO10 in both sexes, depletion of skin Cyanobacteria decreased the protection from hypertension in the RNO5 female, but not male, congenic strain. In conclusion, to our knowledge this is the first study to identify specific microbiota in sites other than gut as contributors to BP regulation. Notably, both oral Actinobacteria and skin Cyanobacteria were beneficial for lowering BP.

Broad cellular functions and diseases including muscular dystrophy, arrhythmogenic right ventricular cardiomyopathy (ARVC5) and cancer are associated with transmembrane protein43 (TMEM43/LUMA). The study aimed to investigate biological roles of TMEM43 through genetic regulation, gene pathways and gene networks, candidate interacting genes, and up- or downstream regulators. Cardiac transcriptomes from 40 strains of recombinant inbred BXD mice and two parental strains representing murine genetic reference population (GRP) were applied for genetic correlation, functional enrichment, and coexpression network analysis using systems genetics approach. The results were validated in a newly created knock-in Tmem43-S358L mutation mouse model (Tmem43S358L) that displayed signs of cardiac dysfunction, resembling ARVC5 phenotype seen in humans. We found high Tmem43 levels among BXDs with broad variability in expression. Expression of Tmem43 highly negatively correlated with heart mass and heart rate among BXDs, whereas levels of Tmem43 highly positively correlated with plasma high-density lipoproteins (HDL). Through finding differentially expressed genes (DEGs) between Tmem43S358L mutant and wild-type (Tmem43WT) lines, 18 pathways (out of 42 found in BXDs GRP) that are involved in ARVC, hypertrophic cardiomyopathy, dilated cardiomyopathy, nonalcoholic fatty liver disease, Alzheimer's disease, Parkinson's disease, and Huntington's disease were verified. We further constructed Tmem43-mediated gene network, in which Ctnna1, Adcy6, Gnas, Ndufs6, and Uqcrc2 were significantly altered in Tmem43S358L mice versus Tmem43WT controls. Our study defined the importance of Tmem43 for cardiac- and metabolism-related pathways, suggesting that cardiovascular disease-relevant risk factors may also increase risk of metabolic and neurodegenerative diseases via TMEM43-mediated pathways.

People living with HIV (PLWH) have increased prevalence of comorbid conditions including insulin resistance and at-risk alcohol use. Circulating microRNAs (miRs) may serve as minimally invasive indicators of pathophysiological states. We aimed to identify whether alcohol modulates circulating miR associations with measures of glucose/insulin dynamics in PLWH. PLWH (n = 96; 69.8% males) enrolled in the Alcohol & Metabolic Comorbidities in PLWH: Evidence-Driven Interventions (ALIVE-Ex) study were stratified into negative phosphatidylethanol (PEth < 8 ng/mL, n = 42) and positive PEth (PEth ≥ 8 ng/mL, n = 54) groups. An oral glucose tolerance test (OGTT) was administered, and total RNA was isolated from fasting plasma to determine absolute miR expression. Circulating miRs were selected based on their role in skeletal muscle (miR-133a and miR-206), pancreatic β-cell (miR-375), liver (miR-20a), and adipose tissue (miR-let-7b, miR-146a, and miR-221) function. Correlation and multiple regression analyses between miR expression and adiponectin, 2 h glucose, insulin, and C-peptide values were performed adjusting for body mass index (BMI) category, age, sex, and viral load. miR-133a was negatively associated with adiponectin (P = 0.002) in the negative PEth group, and miR-20a was positively associated with 2 h glucose (P = 0.013) in the positive PEth group. Regression analyses combining miRs demonstrated that miR-133a (P < 0.001) and miR-221 (P = 0.010) together predicted adiponectin in the negative PEth group. miR-20a (P < 0.001) and miR-375 (P = 0.002) together predicted 2 h glucose in the positive PEth group. Our results indicate that associations between miRs and measures of glucose/insulin dynamics differed between PEth groups, suggesting that the pathophysiological mechanisms contributing to altered glucose homeostasis in PLWH are potentially modulated by alcohol use.