Domestication and the storage starch biosynthesis pathway: Signatures of selection from a whole sorghum genome sequencing strategy

Next-generation sequencing of complete genomes has given researchers unprecedented levels of information to study the multifaceted evolutionary changes that have shaped elite plant germplasm. In conjunction with population genetic analytical techniques and detailed online databases, we can more accurately capture the effects of domestication on entire biological pathways of agronomic importance. In this study, we explore the genetic diversity and signatures of selection in all predicted gene models of the storage starch synthesis pathway of Sorghum bicolor, utilizing a diversity panel containing lines categorized as either ‘Landraces’ or ‘Wild and Weedy’ genotypes. Amongst a total of 114 genes involved in starch synthesis, 71 had at least a single signal of purifying selection and 62 a signal of balancing selection and others a mix of both. This included key genes such as STARCH PHOSPHORYLASE 2 (SbPHO2, under balancing selection), PULLULANASE (SbPUL, under balancing selection) and ADP-glucose pyrophosphorylases (SHRUNKEN2, SbSH2 under purifying selection). Effectively, many genes within the primary starch synthesis pathway had a clear reduction in nucleotide diversity between the Landraces and wild and weedy lines indicating that the ancestral effects of domestication are still clearly identifiable. There was evidence of the positional rate variation within the well-characterized primary starch synthesis pathway of sorghum, particularly in the Landraces, whereby low evolutionary rates upstream and high rates downstream in the metabolic pathway were expected. This observation did not extend to the wild and weedy lines or the minor starch synthesis pathways.

Domestication and the storage starch biosynthesis pathway: signatures of selection from a whole sorghum genome sequencing strategy

Next-generation sequencing of complete genomes has given researchers unprecedented levels of information to study the multifaceted evolutionary changes that have shaped elite plant germplasm. In conjunction with population genetic analytical techniques and detailed online databases, we can more accurately capture the effects of domestication on entire biological pathways of agronomic importance. In this study, we explore the genetic diversity and signatures of selection in all predicted gene models of the storage starch synthesis pathway of Sorghum bicolor, utilizing a diversity panel containing lines categorized as either ‘Landraces’ or ‘Wild and Weedy’ genotypes. Amongst a total of 114 genes involved in starch synthesis, 71 had at least a single signal of purifying selection and 62 a signal of balancing selection and others a mix of both. This included key genes such as STARCH PHOSPHORYLASE 2 (SbPHO2, under balancing selection), PULLULANASE (SbPUL, under balancing selection) and ADP-glucose pyrophosphorylases (SHRUNKEN2, SbSH2 under purifying selection). Effectively, many genes within the primary starch synthesis pathway had a clear reduction in nucleotide diversity between the Landraces and wild and weedy lines indicating that the ancestral effects of domestication are still clearly identifiable. There was evidence of the positional rate variation within the well-characterized primary starch synthesis pathway of sorghum, particularly in the Landraces, whereby low evolutionary rates upstream and high rates downstream in the metabolic pathway were expected. This observation did not extend to the wild and weedy lines or the minor starch synthesis pathways.

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.

Correct biological timing in Arabidopsis requires multiple light-signaling pathways.

Circadian oscillators provide rhythmic temporal cues for a range of biological processes in plants and animals, enabling anticipation of the day/night cycle and enhancing fitness-associated traits. We have used engineering models to understand the control principles of a plant's response to seasonal variation. We show that the seasonal changes in the timing of circadian outputs require light regulation via feed-forward loops, combining rapid light-signaling pathways with entrained circadian oscillators. Linear time-invariant models of circadian rhythms were computed for 3,503 circadian-regulated genes and for the concentration of cytosolic-free calcium to quantify the magnitude and timing of regulation by circadian oscillators and light-signaling pathways. Bioinformatic and experimental analysis show that rapid light-induced regulation of circadian outputs is associated with seasonal rephasing of the output rhythm. We identify that external coincidence is required for rephasing of multiple output rhythms, and is therefore important in general phase control in addition to specific photoperiod-dependent processes such as flowering and hypocotyl elongation. Our findings uncover a fundamental design principle of circadian regulation, and identify the importance of rapid light-signaling pathways in temporal control.

IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Biological response to phorbol ester determined by alternative G1 pathways.

A plethora of extracellular signals is known to induce a common set of immediate early genes. The immediate early response, therefore, must not be sufficient to determine the biological outcome. An example of this is found with the phorbol ester 12-O-tetradecanoylphorbol 13-acetate (TPA). A potent activator of protein kinase C, TPA can either stimulate or inhibit cell proliferation, depending on the cell type. This cell context-dependent response to TPA is observed with two subclones of NIH 3T3 cells, the P- and the N-3T3 clones. TPA is a mitogen for the P-3T3 but an antimitogen for the N-3T3 cells. The immediate early pathway is activated by TPA in both cell types, indicating that this pathway alone does not activate DNA synthesis. The delayed induction of cyclin D1 expression by TPA is observed only in the P-3T3 cells, correlating with mitogenesis. N-Acetylcysteine does not affect the immediate early pathway but can inhibit the TPA-mediated induction of cyclin D1 and DNA synthesis. In the N-3T3 cells, TPA causes an inhibition of the cyclin E-associated kinase at the G1/S transition, correlating with growth inhibition. The growth-inhibitory activity of TPA is not affected by N-acetylcysteine. Thus, the two TPA-regulated G1 pathways can be distinguished by their sensitivity to N-acetylcysteine. These results demonstrate that TPA can activate alternative G1 pathways. Moreover, the selection of the alternative G1 pathways is determined by the cell context, which, in turn, dictates the biological response to TPA.

Biological and environmental determinants

- describe the complex web of determinants as part of broad causal pathways that affect health - identify and discuss the range of physical, biological and environmental determinants that impact on health - suggest why it is important to the practice of public health that you understand how determinants contribute to health - understand the complexity of health and illness and the multifaceted role of health determinants - relate determinants of health to public health activity and realise the need for multisectoral action and multiple approaches when working to improve health

The effect of silicate ions on proliferation, osteogenic differentiation and cell signalling pathways (WNT and SHH) of bone marrow stromal cells

Silicon (Si) is a trace element, which plays an important role in human bone growth. Si has been incorporated into biomaterials for bone regeneration in order to improve their osteogenic potential, both in vitro and in vivo. Little is known, however, as to how Si ions elicit their biological response on bone-forming cells. The aim of this study was to investigate the effect of Si ions on the proliferation, differentiation, bone-related gene expression and cell signalling pathways of bone marrow stromal cells (BMSCs) by comparing the BMSC responses to different concentrations of NaCl and Na2SiO3, while taking into account and excluding the effect of Na ions. Our study showed that Si ions at a concentration of 0.625 mM significantly enhanced the proliferation, mineralization nodule formation, bone-related gene expression (OCN, OPN and ALP) and bone matrix proteins (ALP and OPN) of BMSCs. Furthermore, Si ions at 0.625 mM could counteract the effect of the WNT inhibitor (W.I.) cardamonin on the osteogenic genes expression, (OPN, OCN and ALP), WNT and SHH signalling pathway-related genes in BMSCs. These results suggest that Si ions by themselves play an important role in regulating the proliferation and osteogenic differentiation of BMSCs, with the involvement of WNT and SHH signalling pathways. Our study provides evidence to explain possible molecular mechanisms whereby Si ions released from Si-containing biomaterials can acquire enhanced bioactivity at desired concentration.

Cellular model systems to study the tumor biological role of kallikrein-related peptidases in ovarian and prostate cancer

Since the identification of the gene family of kallikrein related peptidases (KLKs), their function has been robustly studied at the biochemical level. In vitro biochemical studies have shown that KLK proteases are involved in a number of extracellular processes that initiate intracellular signaling pathways by hydrolysis, as reviewed in Chapters 8, 9, and 15, Volume 1. These events have been associated with more invasive phenotypes of ovarian, prostate, and other cancers. Concomitantly, aberrant expression of KLKs has been associated with poor prognosis of patients with ovarian and prostate cancer (Borgoño and Diamandis, 2004; Clements et al., 2004; Yousef and Diamandis, 2009), with prostate-specific antigen (PSA, KLK3) being a long standing, clinically employed biomarker for prostate cancer (Lilja et al., 2008). Data generated from patient samples in clinical studies, alongwith biochemical activity, suggests that KLKs function in the development and progression of these diseases. To bridge the gap between their function at the molecular level and the clinical need for efficacious treatment and prognostic biomarkers, functional assessment at the in vitro cellular level, using various culture models, is increasing, particularly in a three-dimensional (3D) context (Abbott, 2003; Bissell and Radisky, 2001; Pampaloni et al., 2007; Yamada and Cukierman, 2007).

MS-based metabolomics facilitates the discovery of in vivo functional small molecules with a diversity of biological contexts

In vivo small molecules as necessary intermediates are involved in numerous critical metabolic pathways and biological processes associated with many essential biological functions and events. There is growing evidence that MS-based metabolomics is emerging as a powerful tool to facilitate the discovery of functional small molecules that can better our understanding of development, infection, nutrition, disease, toxicity, drug therapeutics, gene modifications and host-pathogen interaction from metabolic perspectives. However, further progress must still be made in MS-based metabolomics because of the shortcomings in the current technologies and knowledge. This technique-driven review aims to explore the discovery of in vivo functional small molecules facilitated by MS-based metabolomics and to highlight the analytic capabilities and promising applications of this discovery strategy. Moreover, the biological significance of the discovery of in vivo functional small molecules with different biological contexts is also interrogated at a metabolic perspective.

Osteocyte-induced angiogenesis via VEGF–MAPK-dependent pathways in endothelial cells

Recently, it has been suggested osteocytes control the activities of bone formation (osteoblasts) and resorption (osteoclast), indicating their important regulatory role in bone remodelling. However, to date, the role of osteocytes in controlling bone vascularisation remains unknown. Our aim was to investigate the interaction between endothelial cells and osteocytes and to explore the possible molecular mechanisms during angiogenesis. To model osteocyte/endothelial cell interactions, we co-cultured osteocyte cell line (MLOY4) with endothelial cell line (HUVECs). Co-cultures were performed in 1:1 mixture of osteocytes and endothelial cells or by using the conditioned media (CM) transfer method. Real-time cell migration of HUVECs was measured with the transwell migration assay and xCELLigence system. Expression levels of angiogenesis- related genes were measured by quantitative real-time polymerase chain reaction (qRT-PCR). The effect of vascular endothelial growth factor (VEGF) and mitogen-activated phosphorylated kinase (MAPK) signaling were monitored by western blotting using relevant antibodies and inhibitors. During the bone formation, it was noted that osteocyte dendritic processes were closely connected to the blood vessels. The CM generated from MLOY4 cells-activated proliferation, migration, tube-like structure formation, and upregulation of angiogenic genes in endothelial cells suggesting that secretory factor(s) from osteocytes could be responsible for angiogenesis. Furthermore, we identified that VEGF secreted from MLOY4-activated VEGFR2–MAPK–ERK-signaling pathways in HUVECs. Inhibiting VEGF and/or MAPK–ERK pathways abrogated osteocyte-mediated angiogenesis in HUVEC cells. Our data suggest an important role of osteocytes in regulating angiogenesis.