Epidemiology and genetic architecture of blood pressure: a family based study of Generation Scotland

Hypertension is a major risk factor for cardiovascular disease and mortality, and a growing global public health concern, with up to one-third of the world’s population affected. Despite the vast amount of evidence for the benefits of blood pressure (BP) lowering accumulated to date, elevated BP is still the leading risk factor for disease and disability worldwide. It is well established that hypertension and BP are common complex traits, where multiple genetic and environmental factors contribute to BP variation. Furthermore, family and twin studies confirmed the genetic component of BP, with a heritability estimate in the range of 30-50%. Contemporary genomic tools enabling the genotyping of millions of genetic variants across the human genome in an efficient, reliable, and cost-effective manner, has transformed hypertension genetics research. This is accompanied by the presence of international consortia that have offered unprecedentedly large sample sizes for genome-wide association studies (GWASs). While GWAS for hypertension and BP have identified more than 60 loci, variants in these loci are associated with modest effects on BP and in aggregate can explain less than 3% of the variance in BP. The aims of this thesis are to study the genetic and environmental factors that influence BP and hypertension traits in the Scottish population, by performing several genetic epidemiological analyses. In the first part of this thesis, it aims to study the burden of hypertension in the Scottish population, along with assessing the familial aggregation and heritialbity of BP and hypertension traits. In the second part, it aims to validate the association of common SNPs reported in the large GWAS and to estimate the variance explained by these variants. In this thesis, comprehensive genetic epidemiology analyses were performed on Generation Scotland: Scottish Family Health Study (GS:SFHS), one of the largest population-based family design studies. The availability of clinical, biological samples, self-reported information, and medical records for study participants has allowed several assessments to be performed to evaluate factors that influence BP variation in the Scottish population. Of the 20,753 subjects genotyped in the study, a total of 18,470 individuals (grouped into 7,025 extended families) passed the stringent quality control (QC) criteria and were available for all subsequent analysis. Based on the BP-lowering treatment exposure sources, subjects were further classified into two groups. First, subjects with both a self-reported medications (SRMs) history and electronic-prescription records (EPRs; n =12,347); second, all the subjects with at least one medication history source (n =18,470). In the first group, the analysis showed a good concordance between SRMs and EPRs (kappa =71%), indicating that SRMs can be used as a surrogate to assess the exposure to BP-lowering medication in GS:SFHS participants. Although both sources suffer from some limitations, SRMs can be considered the best available source to estimate the drug exposure history in those without EPRs. The prevalence of hypertension was 40.8% with higher prevalence in men (46.3%) compared to women (35.8%). The prevalence of awareness, treatment and controlled hypertension as defined by the study definition were 25.3%, 31.2%, and 54.3%, respectively. These findings are lower than similar reported studies in other populations, with the exception of controlled hypertension prevalence, which can be considered better than other populations. Odds of hypertension were higher in men, obese or overweight individuals, people with a parental history of hypertension, and those living in the most deprived area of Scotland. On the other hand, deprivation was associated with higher odds of treatment, awareness and controlled hypertension, suggesting that people living in the most deprived area may have been receiving better quality of care, or have higher comorbidity levels requiring greater engagement with doctors. These findings highlight the need for further work to improve hypertension management in Scotland. The family design of GS:SFHS has allowed family-based analysis to be performed to assess the familial aggregation and heritability of BP and hypertension traits. The familial correlation of BP traits ranged from 0.07 to 0.20, and from 0.18 to 0.34 for parent-offspring pairs and sibling pairs, respectively. A higher correlation of BP traits was observed among first-degree relatives than other types of relative pairs. A variance-component model that was adjusted for sex, body mass index (BMI), age, and age-squared was used to estimate heritability of BP traits, which ranged from 24% to 32% with pulse pressure (PP) having the lowest estimates. The genetic correlation between BP traits showed a high correlation between systolic (SBP), diastolic (DBP) and mean arterial pressure (MAP) (G: 81% to 94%), but lower correlations with PP (G: 22% to 78%). The sibling recurrence risk ratio (λS) for hypertension and treatment were calculated as 1.60 and 2.04 respectively. These findings confirm the genetic components of BP traits in GS:SFHS, and justify further work to investigate genetic determinants of BP. Genetic variants reported in the recent large GWAS of BP traits were selected for genotyping in GS:SFHS using a custom designed TaqMan® OpenArray®. The genotyping plate included 44 single nucleotide polymorphisms (SNPs) that have been previously reported to be associated with BP or hypertension at genome-wide significance level. A linear mixed model that is adjusted for age, age-squared, sex, and BMI was used to test for the association between the genetic variants and BP traits. Of the 43 variants that passed the QC, 11 variants showed statistically significant association with at least one BP trait. The phenotypic variance explained by these variant for the four BP traits were 1.4%, 1.5%, 1.6%, and 0.8% for SBP, DBP, MAP, and PP, respectively. The association of genetic risk score (GRS) that were constructed from selected variants has showed a positive association with BP level and hypertension prevalence, with an average effect of one mmHg increase with each 0.80 unit increases in the GRS across the different BP traits. The impact of BP-lowering medication on the genetic association study for BP traits has been established, with typical practice of adding a fixed value (i.e. 15/10 mmHg) to the measured BP values to adjust for BP treatment. Using the subset of participants with the two treatment exposure sources (i.e. SRMs and EPRs), the influence of using either source to justify the addition of fixed values in SNP association signal was analysed. BP phenotypes derived from EPRs were considered the true phenotypes, and those derived from SRMs were considered less accurate, with some phenotypic noise. Comparing SNPs association signals between the four BP traits in the two model derived from the different adjustments showed that MAP was the least impacted by the phenotypic noise. This was suggested by identifying the same overlapped significant SNPs for the two models in the case of MAP, while other BP traits had some discrepancy between the two sources

The role of inflammation in vascular function and pregnancy outcome in the pregnant stroke prone spontaneously hypertensive rat

During pregnancy, the maternal cardiovascular system undergoes major adaptation. One of these changes is a 40-50 % increase in circulating blood volume which requires a systemic remodelling of the vasculature in order to regulate maternal blood pressure and maximise blood supply to the developing placenta and fetus. These changes are broadly conserved between humans and rats making them an appropriate pre-clinical model in which to study the underlying mechanisms of pregnancy-dependent cardiovascular remodelling. Whilst women are normally protected against cardiovascular disease; pregnancy marks a period of time where women are susceptible to cardiovascular complications. Cardiovascular disease is the leading cause of maternal mortality in the United Kingdom; in particular hypertensive conditions are among the most common complications of pregnancy. One of the main underlying pathologies of these pregnancy complications is thought to be a failure of the maternal cardiovascular system to adapt. The remodelling of the uterine arteries, which directly supply the maternal-fetal interface, is paramount to a healthy pregnancy. Failure of the uterine arteries to remodel sufficiently can result in a number of obstetric complications such as preeclampsia, fetal growth restriction and spontaneous pregnancy loss. At present, it is poorly understood whether this deficient vascular response is due to a predisposition from existing maternal cardiovascular risk factors, the physiological changes that occur during pregnancy or a combination of both. Previous work in our group employed the stroke prone spontaneously hypertensive rat (SHRSP) as a model to investigate pregnancy-dependent remodelling of the uterine arteries. The SHRSP develops hypertension from 6 weeks of age and can be contrasted with the control strain, the Wistar Kyoto (WKY) rat. The phenotype of the SHRSP is therefore reflective of the clinical situation of maternal chronic hypertension during pregnancy. We showed that the SHRSP exhibited a deficient uterine artery remodelling response with respect to both structure and function accompanied by a reduction in litter size relative to the WKY at gestational day (GD) 18. A previous intervention study using nifedipine in the SHRSP achieved successful blood pressure reduction from 6 weeks of age and throughout pregnancy; however uterine artery remodelling and litter size at GD18 was not improved. We concluded that the abnormal uterine artery remodelling present in the SHRSP was independent of chronic hypertension. From these findings, we hypothesised that the SHRSP could be a novel model of spontaneously deficient uterine artery remodelling in response to pregnancy which was underpinned by other as yet unidentified cardiovascular risk factors. In Chapter 1 of this thesis, I have characterised the maternal, placental and fetal phenotype in pregnant (GD18) SHRSP and WKY. The pregnant SHRSP exhibit features of left ventricular hypertrophy in response to pregnancy and altered expression of maternal plasma biomarkers which have been previously associated with hypertension in human pregnancy. I developed a protocol for accurate dissection of the rat uteroplacental unit using qPCR probes specific for each layer. This allowed me to make an accurate and specific statement about gene expression in the SHRSP GD18 placenta; where oxidative stress related gene markers were increased in the vascular compartments. The majority of SHRSP placenta presented at GD18 with a blackened ring which encircled the tissue. Further investigation of the placenta using western blot for caspase 3 cleavage determined that this was likely due to increased cell death in the SHRSP placenta. The SHRSP also presented with a loss of one particular placental cell type at GD18: the glycogen cells. These cells could have been the target of cell death in the SHRSP placenta or were utilised early in pregnancy as a source of energy due to the deficient uterine artery blood supply. Blastocyst implantation was not altered but resorption rate was increased between SHRSP and WKY; indicating that the reduction in litter size in the SHRSP was primarily due to late (>GD14) pregnancy loss. Fetal growth was not restricted in SHRSP which led to the conclusion that SHRSP sacrifice part of their litter to deliver a smaller number of healthier pups. Activation of the immune system is a common pathway that has been implicated in the development of both hypertension and adverse pregnancy outcome. In Chapter 2, I proposed that this may be a mechanism of interest in SHRSP pregnancy and measured the pro-inflammatory cytokine, TNFα, as a marker of inflammation in pregnant SHRSP and WKY and in the placentas from these animals. TNFα was up-regulated in maternal plasma and urine from the GD18 SHRSP. In addition, TNFα release was increased from the GD18 SHRSP placenta as was the expression of the pro-inflammatory TNFα receptor 1 (Tnfr1). In order to investigate whether this excess TNFα was detrimental to SHRSP pregnancy, a vehicle-controlled intervention study using etanercept (a monoclonal antibody which works as a TNFα antagonist) was carried out. Etanercept treatment at GD0, 6, 12 and 18 resulted in an improvement in pregnancy outcome in the SHRSP with an increased litter size and reduced resorption rate. Furthermore, there was an improved uterine artery function in GD18 SHRSP treated with etanercept which was associated with an improved uterine artery blood flow over the course of gestation. In Chapter 3, I sought to identify the source of this detrimental excess of TNFα by designing a panel for maternal leukocytes in the blood and placenta at GD18. A population of CD3- CD161+ cells, which are defined as rat natural killer (NK) cells, were increased in number in the SHRSP. Intracellular flow cytometry also identified this cell type as a source of excess TNFα in blood and placenta from pregnant SHRSP. I then went on to evaluate the effects of etanercept treatment on these CD3- CD161+ cells and showed that etanercept reduced the expression of CD161 and the cytotoxic molecule, granzyme B, in the NK cells. Thus, etanercept limits the cytotoxicity and potential damaging effect of these NK cells in the SHRSP placenta. Analysing the urinary peptidome has clinical potential to identify novel pathways involved with disease and/or to develop biomarker panels to aid and stratify diagnosis. In Chapter 4, I utilised the SHRSP as a pre-clinical model to identify novel urinary peptides associated with hypertensive pregnancy. Firstly, a characterisation study was carried out in the kidney of the WKY and SHRSP. Urine samples from WKY and SHRSP taken at pre-pregnancy, mid-pregnancy (GD12) and late pregnancy (GD18) were used in the peptidomic screen. In order to capture peptides which were markers of hypertensive pregnancy from the urinary peptidomic data, I focussed on those that were only changed in a strain dependent manner at GD12 and 18 and not pre-pregnancy. Peptide fragments from the uromodulin protein were identified from this analysis to be increased in pregnant SHRSP relative to pregnant WKY. This increase in uromodulin was validated at the SHRSP kidney level using qPCR. Uromodulin has previously been identified to be a candidate molecule involved in systemic arterial hypertension but not in hypertensive pregnancy thus is a promising target for further study. In summary, we have characterised the SHRSP as the first model of maternal chronic hypertension during pregnancy and identified that inflammation mediated by TNFα and NK cells plays a key role in the pathology. The evidence presented in this thesis establishes the SHRSP as a pre-clinical model for pregnancy research and can be continued into clinical studies in pregnant women with chronic hypertension which remains an area of unmet research need.

Filling in gaps of Drosophila melanogaster urate degradation metabolic pathway using metabolomics approaches: towards the core metabolome of the fruit fly

The primary goal of systems biology is to integrate complex omics data, and data obtained from traditional experimental studies in order to provide a holistic understanding of organismal function. One way of achieving this aim is to generate genome-scale metabolic models (GEMs), which contain information on all metabolites, enzyme-coding genes, and biochemical reactions in a biological system. Drosophila melanogaster GEM has not been reconstructed to date. Constraint-free genome-wide metabolic model of the fruit fly has been reconstructed in our lab, identifying gaps, where no enzyme was identified and metabolites were either only produced or consume. The main focus of the work presented in this thesis was to develop a pipeline for efficient gap filling using metabolomics approaches combined with standard reverse genetics methods, using 5-hydroxyisourate hydrolase (5-HIUH) as an example. 5-HIUH plays a role in urate degradation pathway. Inability to degrade urate can lead to inborn errors of metabolism (IEMs) in humans, including hyperuricemia. Based on sequence analysis Drosophila CG30016 gene was hypothesised to encode 5- HIUH. CG30016 knockout flies were examined to identify Malpighian tubules phenotype, and shortened lifespan might reflect kidney disorders in hyperuricemia in humans. Moreover, LC-MS analysis of mutant tubules revealed that CG30016 is involved in purine metabolism, and specifically urate degradation pathway. However, the exact role of the gene has not been identified, and the complete method for gap filling has not been developed. Nevertheless, thanks to the work presented here, we are a step closer towards the development of a gap-filling pipeline in Drosophila melanogaster GEM. Importantly, the areas that require further optimisation were identified and are the focus of future research. Moreover, LC-MS analysis confirmed that tubules rather than the whole fly were more suitable for metabolomics analysis of purine metabolism. Previously, Dow/Davies lab has generated the most complete tissue-specific transcriptomic atlas for Drosophila – FlyAtlas.org, which provides data on gene expression across multiple tissues of adult fly and larva. FlyAtlas revealed that transcripts of many genes are enriched in specific Drosophila tissues, and that it is possible to deduce the functions of individual tissues within the fly. Based on FlyAtlas data, it has become clear that the fly (like other metazoan species) must be considered as a set of tissues, each 2 with its own distinct transcriptional and functional profile. Moreover, it revealed that for about 30% of the genome, reverse genetic methods (i.e. mutation in an unknown gene followed by observation of phenotype) are only useful if specific tissues are investigated. Based on the FlyAtlas findings, we aimed to build a primary tissue-specific metabolome of the fruit fly, in order to establish whether different Drosophila tissues have different metabolomes and if they correspond to tissue-specific transcriptome of the fruit fly (FlyAtlas.org). Different fly tissues have been dissected and their metabolome elucidated using LC-MS. The results confirmed that tissue metabolomes differ significantly from each other and from the whole fly, and that some of these differences can be correlated to the tissue function. The results illustrate the need to study individual tissues as well as the whole organism. It is clear that some metabolites that play an important role in a given tissue might not be detected in the whole fly sample because their abundance is much lower in comparison to other metabolites present in all tissues, which prevent the detection of the tissue-specific compound.