MicroRNA modulation of aldosterone production in the adrenal gland

Hypertension is the major risk factor for coronary disease worldwide. Primary hypertension is idiopathic in origin but is thought to arise from multiple risk factors including genetic, lifestyle and environmental influences. Secondary hypertension has a more definite aetiology; its major single cause is primary aldosteronism (PA), the greatest proportion of which is caused by aldosteroneproducing adenoma (APA), where aldosterone is synthesized at high levels by an adenoma of the adrenal gland. There is strong evidence to show that high aldosterone levels cause adverse effects on cardiovascular, cerebrovascular, renal and other systems. Extensive studies have been conducted to analyse the role that regulation of CYP11B2, the gene encoding the aldosterone synthase enzyme plays in determining aldosterone production and the development of hypertension. One significant regulatory factor that has only recently emerged is microRNA (miRNA). miRNAs are small non-coding RNAs, synthesized by a series of enzymatic processes, that negatively regulate gene expression at the posttranscriptional level. Detection and manipulation of miRNA is now known to be a viable method in the treatment, prevention and prognosis of certain diseases. The aim of the present study was to identify miRNAs likely to have a role in the regulation of corticosteroid biosynthesis. To achieve this, the miRNA profile of APA and normal human adrenal tissue was compared, as was the H295R adrenocortical cell line model of adrenocortical function, under both basal conditions and following stimulation of aldosterone production. Key differentially-expressed miRNAs were then identified and bioinformatic tools used to identify likely mRNA targets and pathways for these miRNAs, several of which were investigated and validated using in vitro methods. The background to this study is set out in Chapter 1 of this thesis, followed by a description of the major technical methods employed in Chapter 2. Chapter 3 presents the first of the study results, analysing differences in miRNA profile between APA and normal human adrenal tissue. Microarray was implemented to detect the expression of miRNAs in these two tissue types and several miRNAs were found to vary significantly and consistently between them. Furthermore, members of several miRNA clusters exhibited similar changes in expression pattern between the two tissues e.g. members of cluster miR-29b-1 (miR-29a-3p and miR-29b-3p) and of cluster miR-29b-2 (miR-29b-3p and miR-29c- 3p) are downregulated in APA, while members of cluster let-7a-1 (let-7a-5p and let-7d-5p), cluster let-7a-3 (let-7a-5p and let-7b-5p) and cluster miR-134 (miR- 134 and miR-382) are upregulated. Further bioinformatic analysis explored the possible biological function of these miRNAs using Ingenuity® Systems Pathway Analysis software. This led to the identification of validated mRNAs already known to be targeted by these miRNAs, as well as the prediction of other mRNAs that are likely targets and which are involved in processes relevant to APA pathology including cholesterol synthesis (HMGCR) and corticosteroidogenesis (CYP11B2). It was therefore hypothesised that increases in miR-125a-5p or miR- 335-5p would reduce HMGCR and CYP11B2 expression. Chapter 4 describes the characterisation of H295R cells of different strains and sources (H295R Strain 1, 2, 3 and HAC 15). Expression of CYP11B2 was assessed following application of 3 different stimulants: Angio II, dbcAMP and KCl. The most responsive strain to stimulation was Strain 1 at lower passage numbers. Furthermore, H295R proliferation increased following Angio II stimulation. In Chapter 5, the hypothesis that increases in miR-125a-5p or miR-335-5p reduces HMGCR and CYP11B2 expression was tested using realtime quantitative RT-PCR and transfection of miRNA mimics and inhibitors into the H295R cell line model of adrenocortical function. In this way, miR-125a-5p and miR-335-5p were shown to downregulate CYP11B2 and HMGCR expression, thereby validating certain of the bioinformatic predictions generated in Chapter 3. The study of miRNA profile in the H295R cell lines was conducted in Chapter 6, analysing how it changes under conditions that increase aldosterone secretion, including stimulation Angiotensin II, potassium chloride or dibutyryl cAMP (as a substitute for adrenocorticotropic hormone). miRNA profiling identified 7 miRNAs that are consistently downregulated by all three stimuli relative to basal cells: miR-106a-5p, miR-154-3p, miR-17-5p, miR-196b-5p, miR-19a-3p, miR-20b- 5p and miR-766-3p. These miRNAs include those derived from cluster miR-106a- 5p/miR-20b-5p and cluster miR-17-5p/miR-19a-3p, each producing a single polycistronic transcript. IPA bioinformatic analysis was again applied to identify experimentally validated and predicted mRNA targets of these miRNAs and the key biological pathways likely to be affected. This predicted several interactions between miRNAs derived from cluster miR-17-5p/miR-19a-3p and important mRNAs involved in cholesterol biosynthesis: LDLR and ABCA1. These predictions were investigated by in vitro experiment. miR-17-5p/miR-106a-p and miR-20b-5p were found to be consistently downregulated by stimulation of aldosterone biosynthesis. Moreover, miR-766-3p was upregulation throughout. Furthermore, I was able to validate the downregulation of LDLR by miR-17 transfection, as predicted by IPA. In summary, this study identified key miRNAs that are differentially-expressed in vivo in cases of APA or in vitro following stimulation of aldosterone biosynthesis. The many possible biological actions these miRNAs could have were filtered by bioinformatic analysis and selected interactions validated in vitro. While direct actions of these miRNAs on steroidogenic enzymes were identified, cholesterol handling also emerged as an important target and may represent a useful point of intervention in future therapies designed to modulate aldosterone biosynthesis and reduce its harmful effects.

Reconstruction and analysis of the NF-kB pathway interactome: a systems biology approach

Background. One of the phenomena observed in human aging is the progressive increase of a systemic inflammatory state, a condition referred to as “inflammaging”, negatively correlated with longevity. A prominent mediator of inflammation is the transcription factor NF-kB, that acts as key transcriptional regulator of many genes coding for pro-inflammatory cytokines. Many different signaling pathways activated by very diverse stimuli converge on NF-kB, resulting in a regulatory network characterized by high complexity. NF-kB signaling has been proposed to be responsible of inflammaging. Scope of this analysis is to provide a wider, systemic picture of such intricate signaling and interaction network: the NF-kB pathway interactome. Methods. The study has been carried out following a workflow for gathering information from literature as well as from several pathway and protein interactions databases, and for integrating and analyzing existing data and the relative reconstructed representations by using the available computational tools. Strong manual intervention has been necessarily used to integrate data from multiple sources into mathematically analyzable networks. The reconstruction of the NF-kB interactome pursued with this approach provides a starting point for a general view of the architecture and for a deeper analysis and understanding of this complex regulatory system. Results. A “core” and a “wider” NF-kB pathway interactome, consisting of 140 and 3146 proteins respectively, were reconstructed and analyzed through a mathematical, graph-theoretical approach. Among other interesting features, the topological characterization of the interactomes shows that a relevant number of interacting proteins are in turn products of genes that are controlled and regulated in their expression exactly by NF-kB transcription factors. These “feedback loops”, not always well-known, deserve deeper investigation since they may have a role in tuning the response and the output consequent to NF-kB pathway initiation, in regulating the intensity of the response, or its homeostasis and balance in order to make the functioning of such critical system more robust and reliable. This integrated view allows to shed light on the functional structure and on some of the crucial nodes of thet NF-kB transcription factors interactome. Conclusion. Framing structure and dynamics of the NF-kB interactome into a wider, systemic picture would be a significant step toward a better understanding of how NF-kB globally regulates diverse gene programs and phenotypes. This study represents a step towards a more complete and integrated view of the NF-kB signaling system.

How do our 'quicks' generate pace? A cross sectional analysis of the Cricket Australia pace pathway

Technique and physical contributions to ball delivery speed in fast bowling have been popular research topics in sports science. However, a common limiting factor of this work is the level of expertise of participants and lack of within bowler investigations (Salter et al., 2007). The relationship between technique, anthropometry and ball speed has not been comprehensively investigated among elite fast bowlers. The purpose of this study was to examine the relationship between technique, anthropometric variables and ball speed using both within- and betweenbowler analyses in a cross section of the Cricket Australia high performance pace pathway.

Genomic organisation, activity and distribution analysis of the microbial putrescine oxidase degradation pathway

The catalytic action of putrescine specific amine oxidases acting in tandem with 4-aminobutyraldehyde dehydrogenase is explored as a degradative pathway in Rhodococcus opacus. By limiting the nitrogen source, increased catalytic activity was induced leading to a coordinated response in the oxidative deamination of putrescine to 4-aminobutyraldehyde and subsequent dehydrogenation to 4-aminobutyrate. Isolating the dehydrogenase by ion exchange chromatography and gel filtration revealed that the enzyme acts principally on linear aliphatic aldehydes possessing an amino moiety. Michaelis-Menten kinetic analysis delivered a Michaelis constant (KM=0.014mM) and maximum rate (Vmax=11.2μmol/min/mg) for the conversion of 4-aminobutyraldehyde to 4-aminobutyrate. The dehydrogenase identified by MALDI-TOF mass spectrometric analysis (E value=0.031, 23% coverage) belongs to a functionally related genomic cluster that includes the amine oxidase, suggesting their association in a directed cell response. Key regulatory, stress and transport encoding genes have been identified, along with candidate dehydrogenases and transaminases for the further conversion of 4-aminobutyrate to succinate. Genomic analysis has revealed highly similar metabolic gene clustering among members of Actinobacteria, providing insight into putrescine degradation notably among Micrococcaceae, Rhodococci and Corynebacterium by a pathway that was previously uncharacterised in bacteria.

Signaling pathway genes for blood pressure, folate and cholesterol levels among hypertensives : an epistasis analysis

Irregular atrial pressure, defective folate and cholesterol metabolism contribute to the pathogenesis of hypertension. However, little is known about the combined roles of the methylenetetrahydrofolate reductase (MTHFR), apolipoprotein-E (ApoE) and angiotensin-converting enzyme (ACE) genes, which are involved in metabolism and homeostasis. The objective of this study is to investigate the association of the MTHFR 677 C>T and 1298A>C, ACE insertion–deletion (I/D) and ApoE genetic polymorphisms with hypertension and to further explore the epistasis interactions that are involved in these mechanisms. A total of 594 subjects, including 348 normotensive and 246 hypertensive ischemic stroke subjects were recruited. The MTHFR 677 C>T and 1298A>C, ACE I/D and ApoEpolymorphisms were genotyped and the epistasis interaction were analyzed. The MTHFR 677 C>T and ApoE polymorphisms demonstrated significant associations with susceptibility to hypertension in multiple logistic regression models, multifactor dimensionality reduction and a classification and regression tree. In addition, the logistic regression model demonstrated that significant interactions between the ApoE E3E3, E2E4, E2E2 and MTHFR 677 C>T polymorphisms existed. In conclusion, the results of this epistasis study indicated significant association between the ApoE and MTHFR polymorphisms and hypertension.

targetTB: A target identification pipeline for Mycobacterium tuberculosis through an interactome, reactome and genome-scale structural analysis

Background: Tuberculosis still remains one of the largest killer infectious diseases, warranting the identification of newer targets and drugs. Identification and validation of appropriate targets for designing drugs are critical steps in drug discovery, which are at present major bottle-necks. A majority of drugs in current clinical use for many diseases have been designed without the knowledge of the targets, perhaps because standard methodologies to identify such targets in a high-throughput fashion do not really exist. With different kinds of 'omics' data that are now available, computational approaches can be powerful means of obtaining short-lists of possible targets for further experimental validation. Results: We report a comprehensive in silico target identification pipeline, targetTB, for Mycobacterium tuberculosis. The pipeline incorporates a network analysis of the protein-protein interactome, a flux balance analysis of the reactome, experimentally derived phenotype essentiality data, sequence analyses and a structural assessment of targetability, using novel algorithms recently developed by us. Using flux balance analysis and network analysis, proteins critical for survival of M. tuberculosis are first identified, followed by comparative genomics with the host, finally incorporating a novel structural analysis of the binding sites to assess the feasibility of a protein as a target. Further analyses include correlation with expression data and non-similarity to gut flora proteins as well as 'anti-targets' in the host, leading to the identification of 451 high-confidence targets. Through phylogenetic profiling against 228 pathogen genomes, shortlisted targets have been further explored to identify broad-spectrum antibiotic targets, while also identifying those specific to tuberculosis. Targets that address mycobacterial persistence and drug resistance mechanisms are also analysed. Conclusion: The pipeline developed provides rational schema for drug target identification that are likely to have high rates of success, which is expected to save enormous amounts of money, resources and time in the drug discovery process. A thorough comparison with previously suggested targets in the literature demonstrates the usefulness of the integrated approach used in our study, highlighting the importance of systems-level analyses in particular. The method has the potential to be used as a general strategy for target identification and validation and hence significantly impact most drug discovery programmes.

Mutation analysis of TGF-β signaling pathway genes among Finnish patients with primary pulmonary hypertension and hereditary hemorrhagic telangiectasia

Primary pulmonary hypertension (PPH), or according to the recent classification idiopathic pulmonary hypertension (IPAH), is a rare, progressive disease of pulmonary vasculature leading to pulmonary hypertension and right heart failure. Most of the patients are sporadic but in about 6% of cases the disease is familial (FPPH). In 2000 two different groups identified the gene predisposing to PPH. This gene, Bone morphogenetic protein receptor type 2 (BMPR2), encodes a subunit of transforming growth factor β (TGF-β) receptor complex. There is a genetic connection between PPH and hereditary hemorrhagic telangiectasia (HHT), a bleeding disorder characterized by local telangiectasias and sometimes with pulmonary hypertension. In HHT, mutations in ALK1 (activin like kinase type 1) and Endoglin, another members of the TGF-β signaling pathway are found. In this study we identified all of the Finnish PPH patients for the years 1986-1999 using the hospital discharge registries of Finnish university hospitals. During this period we found a total of 59 confirmed PPH patients: 55 sporadic and 4 familial representing 3 different families. In 1999 the prevalence of PPH was 5.8 per million and the annual incidence varied between 0.2-1.3 per million. Among 28 PPH patients studied, heterozygous BMPR2 mutations were found in 12% (3/26) of sporadic patients and in 33% of the PPH families (1/3). All the mutations found were different. Large deletions of BMPR2 were excluded by single-stranded chain polymomorphism analysis. As a candidate gene approach we also studied ALK1, Endoglin, Bone Morphogenetic Receptor Type IA (BMPR1A or ALK3), Mothers Against Decapentaplegic Homolog 4 (SMAD4) and Serotonine Transporter Gene (SLC6A4) using single-strand conformational polymorphism (SSCP) analysis and direct sequencing. Among patients and family members studied, we found two mutations in ALK1 in two unrelated samples. We also identified all the HHT patients treated at the Department of Otorhinolaryngology at Helsinki University Central Hospital between the years of 1990-2005 and 8 of the patients were studied for Endoglin and ALK1 mutations using direct sequencing. A total of seven mutations were found and all the mutations were different. The absence of a founder mutation in the Finnish population in both PPH and HHT was somewhat surprising. This suggests that the mutations of BMPR2, ALK1 and Endoglin are quite young and the older mutations have been lost due to repetitive genetic bottlenecks and/or negative selection. Also, other genes than BMPR2 may be involved in the pathogenesis of PPH. No founder mutations were found in PPH or HHT and thus no simple genetic test is available for diagnostics.

Analysis of the chlorophyll catabolism pathway in leaves of an introgression senescence mutant of Lolium temulentum

Maria Roca, Caron James, Adriana Pruzinsk?, Stefan H?rtensteiner, Howard Thomas and Helen Ougham. Analysis of the chlorophyll catabolism pathway in leaves of an introgression senescence mutant of Lolium temulentum. Phytochemistry, 65 (9), 1231-1238. Sponsorship: BBSRC RAE2008

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Assessment of LD matrix measures for the analysis of biological pathway association.

Complex diseases will have multiple functional sites, and it will be invaluable to understand the cross-locus interaction in terms of linkage disequilibrium (LD) between those sites (epistasis) in addition to the haplotype-LD effects. We investigated the statistical properties of a class of matrix-based statistics to assess this epistasis. These statistical methods include two LD contrast tests (Zaykin et al., 2006) and partial least squares regression (Wang et al., 2008). To estimate Type 1 error rates and power, we simulated multiple two-variant disease models using the SIMLA software package. SIMLA allows for the joint action of up to two disease genes in the simulated data with all possible multiplicative interaction effects between them. Our goal was to detect an interaction between multiple disease-causing variants by means of their linkage disequilibrium (LD) patterns with other markers. We measured the effects of marginal disease effect size, haplotype LD, disease prevalence and minor allele frequency have on cross-locus interaction (epistasis). In the setting of strong allele effects and strong interaction, the correlation between the two disease genes was weak (r=0.2). In a complex system with multiple correlations (both marginal and interaction), it was difficult to determine the source of a significant result. Despite these complications, the partial least squares and modified LD contrast methods maintained adequate power to detect the epistatic effects; however, for many of the analyses we often could not separate interaction from a strong marginal effect. While we did not exhaust the entire parameter space of possible models, we do provide guidance on the effects that population parameters have on cross-locus interaction.

Association analysis of Notch pathway signalling genes in diabetic nephropathy

Several studies have provided compelling evidence implicating the Notch signalling pathway in diabetic nephropathy. Co-regulation of Notch signalling pathway genes with GREM1 has recently been demonstrated and several genes involved in the Notch pathway are differentially expressed in kidney biopsies from individuals with diabetic nephropathy. We assessed single-nucleotide polymorphisms (SNPs; n = 42) in four of these key genes (JAG1, HES1, NOTCH3 and ADAM10) for association with diabetic nephropathy using a case-control design.
Tag SNPs and potentially functional SNPs were genotyped using Sequenom or Taqman technologies in a total of 1371 individuals with type 1 diabetes (668 patients with nephropathy and 703 controls without nephropathy). Patients and controls were white and recruited from the UK and Ireland. Association analyses were performed using PLINK (http://pngu.mgh.harvard.edu/similar to purcell/plink/) and haplotype frequencies in patients and controls were compared. Adjustment for multiple testing was performed by permutation testing.
In analyses stratified by centre, we identified six SNPs, rs8708 and rs11699674 (JAG1), rs10423702 and rs1548555 (NOTCH3), rs2054096 and rs8027998 (ADAM10) as being associated with diabetic nephropathy before, but not after, adjustment for multiple testing. Haplotype and subgroup analysis according to duration of diabetes also failed to find an association with diabetic nephropathy.
Our results suggest that common variants in JAG1, HES1, NOTCH3 and ADAM10 are not strongly associated with diabetic nephropathy in type 1 diabetes among white individuals. Our findings, however, cannot entirely exclude these genes from involvement in the pathogenesis of diabetic nephropathy.

Functional analysis of the Campylobacter jejuni N-linked protein glycosylation pathway

We describe in this report the characterization of the recently discovered N-linked glycosylation locus of the human bacterial pathogen Campylobacter jejuni, the first such system found in a species from the domain Bacteria. We exploited the ability of this locus to function in Escherichia coli to demonstrate through mutational and structural analyses that variant glycan structures can be transferred onto protein indicating the relaxed specificity of the putative oligosaccharyltransferase PglB. Structural data derived from these variant glycans allowed us to infer the role of five individual glycosyltransferases in the biosynthesis of the N-linked heptasaccharide. Furthermore, we show that C. jejuni- and E. coli-derived pathways can interact in the biosynthesis of N-linked glycoproteins. In particular, the E. coli encoded WecA protein, a UDP-GlcNAc: undecaprenylphosphate GlcNAc-1-phosphate transferase involved in glycolipid biosynthesis, provides for an alternative N-linked heptasaccharide biosynthetic pathway bypassing the requirement for the C. jejuni-derived glycosyltransferase PglC. This is the first experimental evidence that biosynthesis of the N-linked glycan occurs on a lipid-linked precursor prior to transfer onto protein. These findings provide a framework for understanding the process of N-linked protein glycosylation in Bacteria and for devising strategies to exploit this system for glycoengineering.

Rapid acidification and alkylation: Redox analysis of the MHC class I pathway

The technique of rapid acidification and alkylation can be used to characterise the redox status of oxidoreductases, and to determine numbers of free cysteine residues within substrate proteins. We have previously used this method to analyse interacting components of the MHC class I pathway, namely ERp57 and tapasin. Here, we have applied rapid acidification alkylation as a novel approach to analysing the redox status of MHC class I molecules. This analysis of the redox status of the MHC class I molecules HLA-A2 and HLA-B27, which is strongly associated with a group of inflammatory arthritic disorders referred to as Spondyloarthropathies, revealed structural and conformational information. We propose that this assay provides a useful tool in the study of in vivo MHC class I structure. (c) 2008 Elsevier B.V. All rights reserved.

Genetic association analysis of vitamin D pathway with obesity traits

Objective:Observational studies have examined the link between vitamin D deficiency and obesity traits. Some studies have reported associations between vitamin D pathway genes such as VDR, GC and CYP27B1 with body mass index (BMI) and waist circumference (WC); however, the findings have been inconsistent. Therefore, we investigated the involvement of vitamin D metabolic pathway genes in obesity-related traits in a large population-based study.Methods:We undertook a comprehensive analysis between 100 tagging single nucleotide polymorphisms (tagSNPs) in genes encoding for DHCR7, CYP2R1, VDBP, CYP27B1, CYP27A1, CYP24A1, VDR and RXRG, and obesity traits in 5224 participants (aged 45 years) in the 1958 British birth cohort (1958BC). We further extended our analyses to investigate the associations between SNPs and obesity traits using the summary statistics from the GIANT (Genetic Investigation of Anthropometric Traits) consortium (n=123 865).Results:In the 1958BC (n=5224), after Bonferroni correction, none of the tagSNPs were associated with obesity traits except for one tagSNP from CYP24A1 that was associated with waist-hip ratio (WHR) (rs2296239, P=0.001). However, the CYP24A1 SNP was not associated with BMI-adjusted WHR (WHRadj) in the 1958BC (rs2296239, P=1.00) and GIANT results (n=123 865, P=0.18). There was also no evidence for an interaction between the tagSNPs and obesity on BMI, WC, WHR and WHRadj in the 1958BC. In the GIANT consortium, none of the tagSNPs were associated with obesity traits.Conclusions:Despite a very large study, our findings suggest that the vitamin D pathway genes are unlikely to have a major role in obesity-related traits in the general population.