Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Systematic replication of smoking disease associations using survey responses and EHR data in the All of Us Research Program

OBJECTIVE

The All of Us Research Program (All of Us) aims to recruit over a million participants to further precision medicine. Essential to the verification of biobanks is a replication of known associations to establish validity. Here, we evaluated how well All of Us data replicated known cigarette smoking associations.

MATERIALS AND METHODS

We defined smoking exposure as follows: (1) an EHR Smoking exposure that used International Classification of Disease codes; (2) participant provided information (PPI) Ever Smoking; and, (3) PPI Current Smoking, both from the lifestyle survey. We performed a phenome-wide association study (PheWAS) for each smoking exposure measurement type. For each, we compared the effect sizes derived from the PheWAS to published meta-analyses that studied cigarette smoking from PubMed. We defined two levels of replication of meta-analyses: (1) nominally replicated: which required agreement of direction of effect size, and (2) fully replicated: which required overlap of confidence intervals.

RESULTS

PheWASes with EHR Smoking, PPI Ever Smoking, and PPI Current Smoking revealed 736, 492, and 639 phenome-wide significant associations, respectively. We identified 165 meta-analyses representing 99 distinct phenotypes that could be matched to EHR phenotypes. At P < .05, 74 were nominally replicated and 55 were fully replicated. At P < 2.68 × 10-5 (Bonferroni threshold), 58 were nominally replicated and 40 were fully replicated.

DISCUSSION

Most phenotypes found in published meta-analyses associated with smoking were nominally replicated in All of Us. Both survey and EHR definitions for smoking produced similar results.

CONCLUSION

This study demonstrated the feasibility of studying common exposures using All of Us data.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

150 risk variants for diverticular disease of intestine prioritize cell types and enable polygenic prediction of disease susceptibility

We conducted a genome-wide association study (GWAS) analysis of diverticular disease (DivD) of intestine within 724,372 individuals and identified 150 independent genome-wide significant DNA variants. Integration of the GWAS results with human gut single-cell RNA sequencing data implicated gut myocyte, mesothelial and stromal cells, and enteric neurons and glia in DivD development. Ninety-five genes were prioritized based on multiple lines of evidence, including SLC9A3, a drug target gene of tenapanor used for the treatment of the constipation subtype of irritable bowel syndrome. A DivD polygenic score (PGS) enables effective risk prediction (area under the curve [AUC], 0.688; 95% confidence interval [CI], 0.645-0.732) and the top 20% PGS was associated with ∼3.6-fold increased DivD risk relative to the remaining population. Our statistical and bioinformatic analyses suggest that the mechanism of DivD is through colon structure, gut motility, gastrointestinal mucus, and ionic homeostasis. Our analyses reinforce the link between gastrointestinal disorders and the enteric nervous system through genetics.

Sex Differences in Health Conditions Associated with Sexual Assault in a Large Hospital Population

Introduction

Sexual assault is an urgent public health concern with both immediate and long-lasting health consequences, affecting 44% of women and 25% of men during their lifetimes. Large studies are needed to understand the unique healthcare needs of this patient population.

Methods

We mined clinical notes to identify patients with a history of sexual assault in the electronic health record (EHR) at Vanderbilt University Medical Center (VUMC), a large university hospital in the Southeastern USA, from 1989 to 2021 (N = 3,376,424). Using a phenome-wide case-control study, we identified diagnoses co-occurring with disclosures of sexual assault. We performed interaction tests to examine whether sex modified any of these associations. Association analyses were restricted to a subset of patients receiving regular care at VUMC (N = 833,185).

Results

The phenotyping approach identified 14,496 individuals (0.43%) across the VUMC-EHR with documentation of sexual assault and achieved a positive predictive value of 93.0% (95% confidence interval = 85.6-97.0%), determined by manual patient chart review. Out of 1,703 clinical diagnoses tested across all subgroup analyses, 465 were associated with sexual assault. Sex-by-trauma interaction analysis revealed 55 sex-differential associations and demonstrated increased odds of psychiatric diagnoses in male survivors.

Discussion

This case-control study identified associations between disclosures of sexual assault and hundreds of health conditions, many of which demonstrated sex-differential effects. The findings of this study suggest that patients who have experienced sexual assault are at risk for developing wide-ranging medical and psychiatric comorbidities and that male survivors may be particularly vulnerable to developing mental illness.

Examining Sex-Differentiated Genetic Effects Across Neuropsychiatric and Behavioral Traits

BACKGROUND

The origin of sex differences in prevalence and presentation of neuropsychiatric and behavioral traits is largely unknown. Given established genetic contributions and correlations, we tested for a sex-differentiated genetic architecture within and between traits.

METHODS

Using European ancestry genome-wide association summary statistics for 20 neuropsychiatric and behavioral traits, we tested for sex differences in single nucleotide polymorphism (SNP)-based heritability and genetic correlation (rg < 1). For each trait, we computed per-SNP z scores from sex-stratified regression coefficients and identified genes with sex-differentiated effects using a gene-based approach. We calculated correlation coefficients between z scores to test for shared sex-differentiated effects. Finally, we tested for sex differences in across-trait genetic correlations.

RESULTS

We observed no consistent sex differences in SNP-based heritability. Between-sex, within-trait genetic correlations were high, although <1 for educational attainment and risk-taking behavior. We identified 4 genes with significant sex-differentiated effects across 3 traits. Several trait pairs shared sex-differentiated effects. The top genes with sex-differentiated effects were enriched for multiple gene sets, including neuron- and synapse-related sets. Most between-trait genetic correlation estimates were not significantly different between sexes, with exceptions (educational attainment and risk-taking behavior).

CONCLUSIONS

Sex differences in the common autosomal genetic architecture of neuropsychiatric and behavioral phenotypes are small and polygenic and unlikely to fully account for observed sex-differentiated attributes. Larger sample sizes are needed to identify sex-differentiated effects for most traits. For well-powered studies, we identified genes with sex-differentiated effects that were enriched for neuron-related and other biological functions. This work motivates further investigation of genetic and environmental influences on sex differences.

Clinical laboratory test-wide association scan of polygenic scores identifies biomarkers of complex disease

BACKGROUND

Clinical laboratory (lab) tests are used in clinical practice to diagnose, treat, and monitor disease conditions. Test results are stored in electronic health records (EHRs), and a growing number of EHRs are linked to patient DNA, offering unprecedented opportunities to query relationships between genetic risk for complex disease and quantitative physiological measurements collected on large populations.

METHODS

A total of 3075 quantitative lab tests were extracted from Vanderbilt University Medical Center's (VUMC) EHR system and cleaned for population-level analysis according to our QualityLab protocol. Lab values extracted from BioVU were compared with previous population studies using heritability and genetic correlation analyses. We then tested the hypothesis that polygenic risk scores for biomarkers and complex disease are associated with biomarkers of disease extracted from the EHR. In a proof of concept analyses, we focused on lipids and coronary artery disease (CAD). We cleaned lab traits extracted from the EHR performed lab-wide association scans (LabWAS) of the lipids and CAD polygenic risk scores across 315 heritable lab tests then replicated the pipeline and analyses in the Massachusetts General Brigham Biobank.

RESULTS

Heritability estimates of lipid values (after cleaning with QualityLab) were comparable to previous reports and polygenic scores for lipids were strongly associated with their referent lipid in a LabWAS. LabWAS of the polygenic score for CAD recapitulated canonical heart disease biomarker profiles including decreased HDL, increased pre-medication LDL, triglycerides, blood glucose, and glycated hemoglobin (HgbA1C) in European and African descent populations. Notably, many of these associations remained even after adjusting for the presence of cardiovascular disease and were replicated in the MGBB.

CONCLUSIONS

Polygenic risk scores can be used to identify biomarkers of complex disease in large-scale EHR-based genomic analyses, providing new avenues for discovery of novel biomarkers and deeper understanding of disease trajectories in pre-symptomatic individuals. We present two methods and associated software, QualityLab and LabWAS, to clean and analyze EHR labs at scale and perform a Lab-Wide Association Scan.

Epidemiology of Functional Seizures Among Adults Treated at a University Hospital

Importance

Functional seizures (formerly psychogenic nonepileptic seizures), paroxysmal episodes that are often similar to epileptic seizures in their clinical presentation and display no aberrant brain electrical patterns, are understudied. Patients experience a long diagnostic delay, few treatment modalities, a high rate of comorbidities, and significant stigma due to the lack of knowledge about functional seizures.

Objective

To characterize the clinical epidemiology of a population of patients with functional seizures observed at Vanderbilt University Medical Center (VUMC).

Design, Setting, and Participants

This case-control study included patients with functional seizures identified in the VUMC electronic health record (VUMC-EHR) system from October 1989 to October 2018. Patients with epilepsy were excluded from the study and all remaining patients in the VUMC medical center system were used as controls. In total, the study included 1431 patients diagnosed with functional seizures, 2251 with epilepsy and functional seizures, 4715 with epilepsy without functional seizures, and 502 200 control patients who received treatment at VUMC for a minimum of a 3 years. Data were analyzed from November 2018 to March 2020.

Exposure

Diagnosis of functional seizures, as identified from the VUMC-EHR system by an automated phenotyping algorithm that incorporated International Classification of Diseases, Ninth Revision (ICD-9) codes, International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10) codes, Current Procedural Terminology codes, and natural language processing.

Main Outcomes and Measures

Associations of functional seizures with comorbidities and risk factors, measured in odds ratios (ORs).

Results

Of 2 346 808 total patients in the VUMC-EHR aged 18 years or older, 3341 patients with functional seizures were identified (period prevalence, 0.14%), 1062 (74.2%) of whom were women and for which the median (interquartile range) age was 49.3 (39.4-59.9) years. This assessment replicated previously reported associations with psychiatric disorders including posttraumatic stress disorder (PTSD) (OR, 1.22; 95% CI, 1.21-1.24; P < 3.02 × 10-5), anxiety (OR, 1.14; 95% CI, 1.13-1.15; P < 3.02 × 10-5), and depression (OR, 1.14; 95% CI, 1.13-1.15; P < 3.02 × 10-5), and identified novel associations with cerebrovascular disease (OR, 1.08; 95% CI, 1.06-1.09; P < 3.02 × 10-5). An association was found between functional seizures and the known risk factor sexual assault trauma (OR, 10.26; 95% CI, 10.09-10.44; P < 3.02 × 10-5), and sexual assault trauma was found to mediate nearly a quarter of the association between female sex and functional seizures in the VUMC-EHR.

Conclusions and Relevance

This case-control study found evidence to support previously reported associations, discovered new associations between functional seizures and PTSD, anxiety, and depression. An association between cerebrovascular disease and functional seizures was also found. Results suggested that sexual trauma may be a mediating factor in the association between female sex and functional seizures.

Potential positive and negative consequences of ZnT8 inhibition

SLC30A8 encodes the zinc transporter ZnT8. SLC30A8 haploinsufficiency protects against type 2 diabetes (T2D), suggesting that ZnT8 inhibitors may prevent T2D. We show here that, while adult chow fed Slc30a8 haploinsufficient and knockout (KO) mice have normal glucose tolerance, they are protected against diet-induced obesity (DIO), resulting in improved glucose tolerance. We hypothesize that this protection against DIO may represent one mechanism whereby SLC30A8 haploinsufficiency protects against T2D in humans and that, while SLC30A8 is predominantly expressed in pancreatic islet beta cells, this may involve a role for ZnT8 in extra-pancreatic tissues. Consistent with this latter concept we show in humans, using electronic health record-derived phenotype analyses, that the 'C' allele of the non-synonymous rs13266634 SNP, which confers a gain of ZnT8 function, is associated not only with increased T2D risk and blood glucose, but also with increased risk for hemolytic anemia and decreased mean corpuscular hemoglobin (MCH). In Slc30a8 KO mice, MCH was unchanged but reticulocytes, platelets and lymphocytes were elevated. Both young and adult Slc30a8 KO mice exhibit a delayed rise in insulin after glucose injection, but only the former exhibit increased basal insulin clearance and impaired glucose tolerance. Young Slc30a8 KO mice also exhibit elevated pancreatic G6pc2 gene expression, potentially mediated by decreased islet zinc levels. These data indicate that the absence of ZnT8 results in a transient impairment in some aspects of metabolism during development. These observations in humans and mice suggest the potential for negative effects associated with T2D prevention using ZnT8 inhibitors.

Pancreatic islet beta cell-specific deletion of G6pc2 reduces fasting blood glucose

The G6PC1, G6PC2 and G6PC3 genes encode distinct glucose-6-phosphatase catalytic subunit (G6PC) isoforms. In mice, germline deletion of G6pc2 lowers fasting blood glucose (FBG) without affecting fasting plasma insulin (FPI) while, in isolated islets, glucose-6-phosphatase activity and glucose cycling are abolished and glucose-stimulated insulin secretion (GSIS) is enhanced at submaximal but not high glucose. These observations are all consistent with a model in which G6PC2 regulates the sensitivity of GSIS to glucose by opposing the action of glucokinase. G6PC2 is highly expressed in human and mouse islet beta cells however, various studies have shown trace G6PC2 expression in multiple tissues raising the possibility that G6PC2 also affects FBG through non-islet cell actions. Using real-time PCR we show here that expression of G6pc1 and/or G6pc3 are much greater than G6pc2 in peripheral tissues, whereas G6pc2 expression is much higher than G6pc3 in both pancreas and islets with G6pc1 expression not detected. In adult mice, beta cell-specific deletion of G6pc2 was sufficient to reduce FBG without changing FPI. In addition, electronic health record-derived phenotype analyses showed no association between G6PC2 expression and phenotypes clearly unrelated to islet function in humans. Finally, we show that germline G6pc2 deletion enhances glycolysis in mouse islets and that glucose cycling can also be detected in human islets. These observations are all consistent with a mechanism by which G6PC2 action in islets is sufficient to regulate the sensitivity of GSIS to glucose and hence influence FBG without affecting FPI.

An Antisense Oligonucleotide Leads to Suppressed Transcription of Hdac2 and Long-Term Memory Enhancement

Knockout of the memory suppressor gene histone deacetylase 2 (Hdac2) in mice elicits cognitive enhancement, and drugs that block HDAC2 have potential as therapeutics for disorders affecting memory. Currently available HDAC2 catalytic activity inhibitors are not fully isoform specific and have short half-lives. Antisense oligonucleotides (ASOs) are drugs that elicit extremely long-lasting, specific inhibition through base pairing with RNA targets. We utilized an ASO to reduce Hdac2 messenger RNA (mRNA) in mice and determined its longevity, specificity, and mechanism of repression. A single injection of the Hdac2-targeted ASO in the central nervous system produced persistent reduction in HDAC2 protein and Hdac2 mRNA levels for 16 weeks. It enhanced object location memory for 8 weeks. RNA sequencing (RNA-seq) analysis of brain tissues revealed that the repression was specific to Hdac2 relative to related Hdac isoforms, and Hdac2 reduction caused alterations in the expression of genes involved in extracellular signal-regulated kinase (ERK) and memory-associated immune signaling pathways. Hdac2-targeted ASOs also suppress a nonpolyadenylated Hdac2 regulatory RNA and elicit direct transcriptional suppression of the Hdac2 gene through stalling RNA polymerase II. These findings identify transcriptional suppression of the target gene as a novel mechanism of action of ASOs.

Polygenic burden in focal and generalized epilepsies

Rare genetic variants can cause epilepsy, and genetic testing has been widely adopted for severe, paediatric-onset epilepsies. The phenotypic consequences of common genetic risk burden for epilepsies and their potential future clinical applications have not yet been determined. Using polygenic risk scores (PRS) from a European-ancestry genome-wide association study in generalized and focal epilepsy, we quantified common genetic burden in patients with generalized epilepsy (GE-PRS) or focal epilepsy (FE-PRS) from two independent non-Finnish European cohorts (Epi25 Consortium, n = 5705; Cleveland Clinic Epilepsy Center, n = 620; both compared to 20 435 controls). One Finnish-ancestry population isolate (Finnish-ancestry Epi25, n = 449; compared to 1559 controls), two European-ancestry biobanks (UK Biobank, n = 383 656; Vanderbilt biorepository, n = 49 494), and one Japanese-ancestry biobank (BioBank Japan, n = 168 680) were used for additional replications. Across 8386 patients with epilepsy and 622 212 population controls, we found and replicated significantly higher GE-PRS in patients with generalized epilepsy of European-ancestry compared to patients with focal epilepsy (Epi25: P = 1.64×10-15; Cleveland: P = 2.85×10-4; Finnish-ancestry Epi25: P = 1.80×10-4) or population controls (Epi25: P = 2.35×10-70; Cleveland: P = 1.43×10-7; Finnish-ancestry Epi25: P = 3.11×10-4; UK Biobank and Vanderbilt biorepository meta-analysis: P = 7.99×10-4). FE-PRS were significantly higher in patients with focal epilepsy compared to controls in the non-Finnish, non-biobank cohorts (Epi25: P = 5.74×10-19; Cleveland: P = 1.69×10-6). European ancestry-derived PRS did not predict generalized epilepsy or focal epilepsy in Japanese-ancestry individuals. Finally, we observed a significant 4.6-fold and a 4.5-fold enrichment of patients with generalized epilepsy compared to controls in the top 0.5% highest GE-PRS of the two non-Finnish European cohorts (Epi25: P = 2.60×10-15; Cleveland: P = 1.39×10-2). We conclude that common variant risk associated with epilepsy is significantly enriched in multiple cohorts of patients with epilepsy compared to controls-in particular for generalized epilepsy. As sample sizes and PRS accuracy continue to increase with further common variant discovery, PRS could complement established clinical biomarkers and augment genetic testing for patient classification, comorbidity research, and potentially targeted treatment.

A salt-sensing kinase in T lymphocytes, SGK1, drives hypertension and hypertensive end-organ damage

We previously showed that angiotensin II (Ang II) increases T cell production of IL-17A, and that mice deficient in IL-17A have blunted hypertension and attenuated renal and vascular dysfunction. It was recently shown that salt enhances IL-17A production from CD4+ T cells via a serum- and glucocorticoid-regulated kinase 1-dependent (SGK1-dependent) pathway. Thus, we tested the hypothesis that SGK1 signaling in T cells promotes hypertension and contributes to end-organ damage. We show that loss of T cell SGK1 results in a blunted hypertensive response to Ang II infusion by 25 mmHg. Importantly, renal and vascular inflammation is abrogated in these mice compared with control mice. Furthermore, mice lacking T cell SGK1 are protected from Ang II-induced endothelial dysfunction and renal injury. Loss of T cell SGK1 also blunts blood pressure and vascular inflammation in response to deoxycorticosterone acetate-salt (DOCA-salt) hypertension. Finally, we demonstrate that the Na+-K+-2Cl- cotransporter 1 (NKCC1) is upregulated in Th17 cells and is necessary for the salt-induced increase in SGK1 and the IL-23 receptor. These studies demonstrate that T cell SGK1 and NKCC1 may be novel therapeutic targets for the treatment of hypertension and identify a potentially new mechanism by which salt contributes to hypertension.

Renal vascular and glomerular pathologies associated with spontaneous hypertension in the nonhuman primate Chlorocebus aethiops sabaeus

Hypertension is a complex, multifactorial disease affecting an estimated 78 million adults in the United States. Despite scientific gains, the etiology of human essential hypertension is unknown and current experimental models do not recapitulate all the behavioral and physiological characteristics of the pathology. Researchers should assess the translational capacity of these models and look to other animal models for the discovery of new therapies. Chlorocebus aethiops sabaeus, the African Green Monkey (AGM), is a nonhuman primate that develops spontaneous hypertension and may provide a novel translational model for the study of hypertension and associated diseases. In a randomly selected group of 424 adult AGMs, 37% (157/424) exhibited systolic blood pressures (SBP) >140 mmHg (SBP: 172.0 ± 2.2 mmHg) and were characterized as hypertensive (HT). 44% (187/424) were characterized as normotensive with SBP <120 mmHg (NT, SBP: 99.6 ± 1.0 mmHg) and the remaining 18% (80/424) as borderline hypertensive (BHT, SBP: 130.6 ± 0.6 mmHg). When compared with NT animals, HT AGMs are older (8.7 ± 0.6 vs. 12.4 ± 0.7 yr, P < 0.05) with elevated heart rates (125.7 ± 2.0 vs. 137.7 ± 2.2 beats/min, P < 0.05). BHT animals had average heart rates of 138.2 ± 3.1 beats/min (P < 0.05 compared with NT) and were 11.00 ± 0.9 yr old. NT and HT animals had similar levels of angiotensinogen gene expression, plasma renin activity, and renal cortical renin content (P > 0.05). HT monkeys exhibit renal vascular remodeling (wall-to-lumen ratio NT 0.11 ± 0.01 vs. HT 0.15 ± 0.02, P < 0.05) and altered glomerular morphology (Bowman's capsular space: NT 30.9 ± 1.9% vs. HT 44.4 ± 3.1%, P < 0.05). The hypertensive AGM provides a large animal model that is highly similar to humans and should be studied to identify novel, more effective targets for the treatment of hypertension.

Abstract P082: Mitochondrial Gene Expression is Decreased in Renal Inner Medulla of Non-human Primates with Spontaneous Hypertension

Mitochondrial gene expression may influence renal function and consequently, long-term blood pressure control. The African Green Monkey (Chlorocebus aethiops sabaeus; AGM) exhibits heritable, spontaneous hypertension and thus is a translational model for the study of human essential hypertension. We hypothesized that renal mitochondrial gene expression in hypertensive AGMs is decreased and may contribute to renal mitochondrial dysfunction in specific kidney regions. AGMs were phenotyped as normotensive (NT, systolic blood pressure; SBP <120 mmHg) or hypertensive (HT, SBP > 140 mmHg) by forearm plethysmography. Gene expression was determined using qRT-PCR with RNA extracted from renal cortex, outer medulla (OM), inner medulla (IM), and liver of 18 HT (mean SBP 166±7 mmHg) and 18 NT (mean SBP 98±3 mmHg) animals. In renal cortical and OM tissue of HT vervets, COX 3 (Complex IV), Cyt B (Complex III), NADH4 (Complex I) and ATP8 (Complex V) expression were similar in NT and HT AGMs. In IM of HT AGMs, COX3, NADH4, ATP8, and CYTB were downregulated by 4.7-fold (p=0.007), 4-fold (p=0.002), 4.1-fold (p=0.006), and 3-fold (p=0.018), respectively. In the liver, COX 3 expression was decreased 1.9-fold (p=0.04), 1.6-fold for ATP8 (p=0.03), and 2.0-fold for CYTB (p=0.01). Expression of SDH (Complex II) and COX4 (Complex IV), nuclear encoded subunits of the OXPHOS chain, was also assessed. SDH expression was up-regulated 8-fold in renal OM (p=0.005) but unchanged in liver, IM, and cortex. COX4 expression however, was down-regulated in OM by 8-fold (p=0.05), in IM by 5-fold (p=0.011) but unchanged in cortex and liver. Gene expression of the mitochondrial transcription factor TFAM was downregulated by 4-fold in renal OM, and unchanged in renal cortex and liver. Citrate synthase activity showed no difference in mitochondrial number between NT and HT AGMs (p=0.73). We suggest that reduced expression of mitochondrially encoded OXPHOS subunits in the renal IM may contribute to the development of hypertension in the AGM. Mitochondrial gene down-regulation in the liver may be a consequence of the systemic hypertension. We conclude that altered mitochondrial gene expression may be a cellular response to increased oxidative stress in the renal inner medulla of HT AGM.

Abstract 219: Altered Renal and Cardiac Morphology in a Non-Human Primate Model of Spontaneous Essential Hypertension

Current models of essential hypertension are limited in comparisons to human pathophysiology due to distinct differences in postural positioning, circadian rhythmicity, and genetic variation. We have identified a unique non-human primate model of spontaneous, essential hypertension in Chlorocebus aethiops sabaeus , the African Green, or vervet, monkey. Using forearm plethysmography under light ketamine sedation (~15mg/kg), adult vervets were categorized as hypertensive (HT; SBP >140mmHg), borderline hypertensive (BHT; 120mmHg < SBP < 140mmHg), or normotensive (NT; SBP < 120 mmHg). Of the 168 males phenotyped, 32% (53 of 168) were HT (average SBP = 168.79±3.29mmHg), 27% (45 of 168) were BHT (average SBP = 129.58±0.88mmHg), and 41% (70 of 168) were NT (average SBP = 100.21±1.65mmHg). In a smaller study of 18 adult females, 1 was found to be HT (SBP = 156mmHg), 2 were BHT (average SBP = 134±2.0mmHg), and 15 were NT (average SBP = 98.06±4.32mmHg). Fixed paraffin-embedded histological sections of vervet cardiac (H&E) and renal (periodic acid Schiff) tissues were digitized and analyzed. Wall/lumen ratios of renal arterioles (average diameter 33.01±0.5μm) were greater in HT (0.33±0.03) compared to NT (0.24±0.01) male monkeys. Bowman’s capsule space was increased in HT (43.14±3.98% area) glomeruli compared to NT (29.92±2.42% area) glomeruli (p<0.05). Indicative of left ventricular hypertrophy, left ventricular cells were decreased in HT (1833.33±65.48 cells/mm 2 ) compared to NT monkeys (2223.61±59.72 cells/mm 2 ; p<0.05). The increased wall/lumen ratios and Bowman’s capsular spaces indicate that renal vascular hypertrophy and glomerular damage occur in this non-human primate model of essential hypertension. Left ventricular hypertrophy combined with renal glomerular and vascular dysfunction provides strong evidence for using this non-human primate model as a translational experimental species.