DNA structural variability as a factor in gene expression and evolution

Redundant DNA can buffer sequence dependent structural deviations from an ideal double helix. Buffering serves a mechanistic function by reducing extraneous conformational effects which could interfere with readout or which would impose energetic constraints on evolution. It also serves an evolutionary function by allowing for gradual variations in conformation-dependent regulation of gene expression. Such gradualism is critical for the rate of evolution. The buffer structure concept provides a new interpretation for repetitive DNA and for exons and introns.

Regulation of cytochrome P-450 gene expression. Studies with a cloned probe

A cDNA clone for cytochrome P-450e, a phenobarbitone-inducible species in rat liver, has been isolated and characterized. With the use of this cloned DNA, an attempt has been initiated to elucidate the factors regulating the cytochrome P-450 gene expression. Inhibitors of heme synthesis such as cobalt chloride and 3-amino-1,2,4-triazole block the induction of cytochrome P-450e by phenobarbitone at the level of transcription. This is evident from the decrease in the rate of synthesis of cytochrome P-450e, a decrease in the levels of specific translatable messenger RNA, a decrease in the specific cytoplasmic and nuclear messenger RNA contents, and nuclear transcription of cytochrome P-450e gene, as revealed by hybridization to the cloned probe, under these conditions. It is proposed that heme is a general regulator of cytochrome P-450 gene expression at the level of transcription, whereas the drug or its metabolite would impart the specificity needed for the induction of a particular species.

A comparative study of tbe regulation of cytochrome P-450 and glutathione transferase gene expression in rat liver

A cDNA clone for the Ya subunit of glutathione transferase from rat liver was constructed in E.coli. The clone hybridized to Ya and Yc subunit messenger RNAs. On the basis of experiments involving cell-free translation and hybridization to the cloned probe, it was shown that prototype inducers of cytochrome P-450 such as phenobarbitone and 3-methylcholanthrene as well as inhibitors such as CoCl2 and 3-amino-l,2,4-triazole enhanced the glutathione transferase (Ya+Yc) messenger RNA contents in rat liver. A comparative study with the induction of cytochrome P-450 (b+e) by phenobarbitone revealed that the drug manifested a striking increase in the nuclear pre-messenger RNAs for the cytochrome at 12 hr, but did not significantly affect the same in the case of glutathione transferase (Ya+Yc). 3-Amino-l, 2,4-tnazole and CoCl- blocked the phenobarbitone mediated increase in cytochrome P-450 (b+e) nuclear pre-messenger RNAs. These compounds did not significantly affect the glutathione transferase (Ya+Yc) nuclear pre-messenger RNA levels. The polysomal, poly (A)- containing messenger RNAs for cytochrome P-450 (b+e) increased by 12–15 fold after phenobarbitone administration, reached a maximum around 16hr and then decreased sharply. In comparison, the increase in the case of glutathione transferase (Ya+Yc) mesenger RNAs was sluggish and steady and a value of 3–4 fold was reached around 24 hr. Run-off transcription rates for cytochrome P-450 (b+e) increased by nearly 15 fold in 4 hr after phenobarbitone administration, whereas the increase for glutathione transferase (Ya+Yc) was only 2.0 fold. At 12 hr after the drug administration, the glutathione transferase (Ya+Yc) transcription rates were near normal. Administration of 3-amino-l,2,4-triazole and CoCl2 blocked the phenobarbitone-mediated increase in the transcription of cytochrome P-450 (b+e) messenger RNAs. These compounds at best had only marginal effects on the transcription of glutathione transferase (Ya+Yc) messenger RNAs. The half-life of cytochrome P-450 (b+e) messenger RNA was estimated to be 3–4 hr, whereas that for glutathione transferase (Ya+Yc) was found to be 8-9 hr. Administration of phenobarbitone enhanced the half-life of glutathione transferase (Ya+Yc) messenger RNA by nearly two fold. It is suggested that while transcription activation may play a primary role in the induction of cytochrome P-450 (b+e), the induction of glutathione transferase (Ya+Yc) may essentially involve stabilization of the messenger RNAs.

Gene expression evidence for off-target effects caused by RNA interference-mediated gene silencing of Ubiquitin-63E in the cattle tick Rhipicephalus microplus.

Knowledge of cattle tick (Rhipicephalus (Boophilus) microplus; Acari: Ixodidae) molecular and cellular pathways has been hampered by the lack of an annotated genome. In addition, most of the tick expressed sequence tags (ESTs) available to date consist of similar to 50% unassigned sequences without predicted functions. The most common approach to address this has been the application of RNA interference (RNAi) methods to investigate genes and their pathways. This approach has been widely adopted in tick research despite minimal knowledge of the tick RNAi pathway and double-stranded RNA (dsRNA) uptake mechanisms. A strong knockdown phenotype of adult female ticks had previously been observed using a 594 bp dsRNA targeting the cattle tick homologue for the Drosophila Ubiquitin-63E gene leading to nil or deformed eggs. A NimbleGen cattle tick custom microarray based on the BmiGI.V2 database of R. microplus ESTs was used to evaluate the expression of mRNAs harvested from ticks treated with the tick Ubiquitin-63E 594 bp dsRNA compared with controls. A total of 144 ESTs including TC6372 (Ubiquitin-63E) were down-regulated with 136 ESTs up-regulated following treatment. The results obtained substantiated the knockdown phenotype with ESTs identified as being associated with ubiquitin proteolysis as well as oogenesis, embryogenesis, fatty acid synthesis and stress responses. A bioinformatics analysis was undertaken to predict off-target effects (OTE) resulting from the in silico dicing of the 594 bp Ubiquitin-63E dsRNA which identified 10 down-regulated ESTs (including TC6372) within the list of differentially expressed probes on the microarrays. Subsequent knockdown experiments utilising 196 and 109 bp dsRNAs, and a cocktail of short hairpin RNAs (shRNA) targeting Ubiquitin-63E, demonstrated similar phenotypes for the dsRNAs but nil effect following shRNA treatment. Quantitative reverse transcriptase PCR analysis confirmed differential expression of TC6372 and selected ESTs. Our study demonstrated the minimisation of predicted OTEs in the shorter dsRNA treatments (similar to 100-200 bp) and the usefulness of microarrays to study knockdown phenotypes.

Molecular Pathways of Disturbed Sleep and Depression: Studies on Adenosine and Gene Expression Patterns

Background: Adenosine is a potent sleep-promoting substance, and one of its targets is the basal forebrain. Fairly little is known about its mechanism of action in the basal forebrain and about the receptor subtype mediating its regulating effects on sleep homeostasis. Homeostatic deficiency might be one of the causes of the profoundly disturbed sleep pattern in major depressive disorder, which could explain the reduced amounts of delta-activity-rich stages 3 and 4. Since major depression has a relatively high heritability, and on the other hand adenosine regulates sleep homeostasis and might also be involved in mood modulation, adenosine-related genes should be considered for their possible contribution to a predisposition for depression and disturbed sleep in humans. Depression is a complex disorder likely involving the abnormal functioning of several genes. Novel target genes which could serve as the possible common substrates for depression and comorbid disturbed sleep should be identified. In this way specific brain areas related to sleep regulation should be studied by using animal model of depression which represents more homogenous phenotype as compared to humans. It is also important to study these brain areas during the development of depressive-like features to understand how early changes could facilitate pathophysiological changes in depression. Aims and methods: We aimed to find out whether, in the basal forebrain, adenosine induces recovery non-rapid eye movement (NREM) sleep after prolonged waking through the A1 or/and A2A receptor subtype. A1 and A2A receptor antagonists were perfused into the rat basal forebrain during 3 h of sleep deprivation, and the amount of NREM sleep and delta power during recovery NREM sleep were analyzed. We then explored whether polymorphisms in genes related to the metabolism, transport and signaling of adenosine could predispose to depression accompanied by signs of disturbed sleep. DNA from 1423 individuals representative of the Finnish population and including controls and cases with depression, depression accompanied by early morning awakenings and depression accompanied by fatigue, was used in the study to investigate the possible association between polymorphisms from adenosine-related genes and cases. Finally to find common molecular substrates of depression and disturbed sleep, gene expression changes were investigated in specific brain areas in the rat clomipramine model of depression. We focused on the basal forebrain of 3-week old clomipramine-treated rats which develop depressive-like symptoms later in adulthood and on the hypothalamus of adult female clomipramine-treated rats. Results: Blocking of the A1 receptor during sleep deprivation resulted in a reduction of the recovery NREM sleep amount and delta power, whereas A2A receptor antagonism had no effect. Polymorphisms in adenosine-related genes SLC29A3 (equilibrative nucleoside transporter type 3) in women and SLC28A1 (concentrative nucleoside transporter type 1) in men associated with depression alone as well as when accompanied by early morning awakenings and fatigue. In Study III the basal forebrain of postnatal rats treated with clomipramine displayed disturbances in gamma-aminobutyric acid (GABA) receptor type A signaling, in synaptic transmission and possible epigenetic changes. CREB1 was identified as a common transcription denominator which also mediates epigenetic regulation. In the hypothalamus the major changes included the expression of genes in GABA-A receptor pathway, K+ channel-related, glutamatergic and mitochondrial genes, as well as an overexpression of genes related to RNA and mRNA processing. Conclusions: Adenosine plays an important role in sleep homeostasis by promoting recovery NREM sleep via the A1 receptor subtype in the basal forebrain. Also adenosine levels might contribute to the risk of depression with disturbed sleep, since the genes encoding nucleoside transporters showed the strongest associations with depression alone and when accompanied by signs of disturbed sleep in both women and men. Sleep and mood abnormalities in major depressive disorder could be a consequence of multiple changes at the transcriptional level, GABA-A receptor signaling and synaptic transmission in sleep-related basal forebrain and the hypothalamus.

Gene Expression Profiling of the Cephalothorax and Eyestalk in Penaeus Monodon during Ovarian Maturation.

In crustaceans, a range of physiological processes involved in ovarian maturation occurs in organs of the cephalothorax including the hepatopancrease, mandibular and Y-organ. Additionally, reproduction is regulated by neuropeptide hormones and other proteins released from secretory sites within the eyestalk. Reproductive dysfunction in captive-reared prawns, Penaeus monodon, is believed to be due to deficiencies in these factors. In this study, we investigated the expression of gene transcripts in the cephalothorax and eyestalk from wild-caught and captive-reared animals throughout ovarian maturation using custom oligonucleotide microarray screening. We have isolated numerous transcripts that appear to be differentially expressed throughout ovarian maturation and between wild-caught and captive-reared animals. In the cephalothorax, differentially expressed genes included the 1,3-beta-D-glucan-binding high-density lipoprotein, 2/3-oxoacyl-CoA thiolase and vitellogenin. In the eyestalk, these include gene transcripts that encode a protein that modulates G-protein coupled receptor activity and another that encodes an architectural transcription factor. Each may regulate the expression of reproductive neuropeptides, such as the crustacean hyperglycaemic hormone and molt-inhibiting hormone. We could not identify differentially expressed transcripts encoding known reproductive neuropeptides in the eyestalk of either wild-caught or captive-reared prawns at any ovarian maturation stage, however, this result may be attributed to low relative expression levels of these transcripts. In summary, this study provides a foundation for the study of target genes involved in regulating penaeid reproduction.

Resistance gene analogues in sugarcane and sorghum and their association with quantitative trait loci for rust resistance

Fifty-four different sugarcane resistance gene analogue (RGA) sequences were isolated, characterized, and used to identify molecular markers linked to major disease-resistance loci in sugarcane. Ten RGAs were identified from a sugarcane stem expressed sequence tag (EST) library; the remaining 44 were isolated from sugarcane stem, leaf, and root tissue using primers designed to conserved RGA motifs. The map location of 31 of the RGAs was determined in sugarcane and compared with the location of quantitative trait loci (QTL) for brown rust resistance. After 2 years of phenotyping, 3 RGAs were shown to generate markers that were significantly associated with resistance to this disease. To assist in the understanding of the complex genetic structure of sugarcane, 17 of the 31 RGAs were also mapped in sorghum. Comparative mapping between sugarcane and sorghum revealed syntenic localization of several RGA clusters. The 3 brown rust associated RGAs were shown to map to the same linkage group (LG) in sorghum with 2 mapping to one region and the third to a region previously shown to contain a major rust-resistance QTL in sorghum. These results illustrate the value of using RGAs for the identification of markers linked to disease resistance loci and the value of simultaneous mapping in sugarcane and sorghum.

Construction of a global map of human gene expression : the process, tools and analysis

This thesis studies human gene expression space using high throughput gene expression data from DNA microarrays. In molecular biology, high throughput techniques allow numerical measurements of expression of tens of thousands of genes simultaneously. In a single study, this data is traditionally obtained from a limited number of sample types with a small number of replicates. For organism-wide analysis, this data has been largely unavailable and the global structure of human transcriptome has remained unknown. This thesis introduces a human transcriptome map of different biological entities and analysis of its general structure. The map is constructed from gene expression data from the two largest public microarray data repositories, GEO and ArrayExpress. The creation of this map contributed to the development of ArrayExpress by identifying and retrofitting the previously unusable and missing data and by improving the access to its data. It also contributed to creation of several new tools for microarray data manipulation and establishment of data exchange between GEO and ArrayExpress. The data integration for the global map required creation of a new large ontology of human cell types, disease states, organism parts and cell lines. The ontology was used in a new text mining and decision tree based method for automatic conversion of human readable free text microarray data annotations into categorised format. The data comparability and minimisation of the systematic measurement errors that are characteristic to each lab- oratory in this large cross-laboratories integrated dataset, was ensured by computation of a range of microarray data quality metrics and exclusion of incomparable data. The structure of a global map of human gene expression was then explored by principal component analysis and hierarchical clustering using heuristics and help from another purpose built sample ontology. A preface and motivation to the construction and analysis of a global map of human gene expression is given by analysis of two microarray datasets of human malignant melanoma. The analysis of these sets incorporate indirect comparison of statistical methods for finding differentially expressed genes and point to the need to study gene expression on a global level.

U12-type Spliceosome: Localization and Effects of Splicing Efficiency on Gene Expression

The removal of non-coding sequences, introns, is an essential part of messenger RNA processing. In most metazoan organisms, the U12-type spliceosome processes a subset of introns containing highly conserved recognition sequences. U12-type introns constitute less than 0,5% of all introns and reside preferentially in genes related to information processing functions, as opposed to genes encoding for metabolic enzymes. It has previously been shown that the excision of U12-type introns is inefficient compared to that of U2-type introns, supporting the model that these introns could provide a rate-limiting control for gene expression. The low efficiency of U12-type splicing is believed to have important consequences to gene expression by limiting the production of mature mRNAs from genes containing U12-type introns. The inefficiency of U12-type splicing has been attributed to the low abundance of the components of the U12-type spliceosome in cells, but this hypothesis has not been proven. The aim of the first part of this work was to study the effect of the abundance of the spliceosomal snRNA components on splicing. Cells with a low abundance of the U12-type spliceosome were found to inefficiently process U12-type introns encoded by a transfected construct, but the expression levels of endogenous genes were not found to be affected by the abundance of the U12-type spliceosome. However, significant levels of endogenous unspliced U12-type intron-containing pre-mRNAs were detected in cells. Together these results support the idea that U12-type splicing may limit gene expression in some situations. The inefficiency of U12-type splicing has also promoted the idea that the U12-type spliceosome may control gene expression, limiting the mRNA levels of some U12-type intron-containing genes. While the identities of the primary target genes that contain U12-type introns are relatively well known, little has previously been known about the downstream genes and pathways potentially affected by the efficiency of U12-type intron processing. Here, the effects of U12-type splicing efficiency on a whole organism were studied in a Drosophila line with a mutation in an essential U12-type spliceosome component. Genes containing U12-type introns showed variable gene-specific responses to the splicing defect, which points to variation in the susceptibility of different genes to changes in splicing efficiency. Surprisingly, microarray screening revealed that metabolic genes were enriched among downstream effects, and that the phenotype could largely be attributed to one U12-type intron-containing mitochondrial gene. Gene expression control by the U12-type spliceosome could thus have widespread effects on metabolic functions in the organism. The subcellular localization of the U12-type spliceosome components was studied as a response to a recent dispute on the localization of the U12-type spliceosome. All components studied were found to be nuclear indicating that the processing of U12-type introns occurs within the nucleus, thus clarifying a question central to the field. The results suggest that the U12-type spliceosome can limit the expression of genes that contain U12-type introns in a gene-specific manner. Through its limiting role in pre-mRNA processing, the U12-type splicing activity can affect specific genetic pathways, which in the case of Drosophila are involved in metabolic functions.

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions.

Gene Expression : From Microarrays to Functional Genomics

The time of the large sequencing projects has enabled unprecedented possibilities of investigating more complex aspects of living organisms. Among the high-throughput technologies based on the genomic sequences, the DNA microarrays are widely used for many purposes, including the measurement of the relative quantity of the messenger RNAs. However, the reliability of microarrays has been strongly doubted as robust analysis of the complex microarray output data has been developed only after the technology had already been spread in the community. An objective of this study consisted of increasing the performance of microarrays, and was measured by the successful validation of the results by independent techniques. To this end, emphasis has been given to the possibility of selecting candidate genes with remarkable biological significance within specific experimental design. Along with literature evidence, the re-annotation of the probes and model-based normalization algorithms were found to be beneficial when analyzing Affymetrix GeneChip data. Typically, the analysis of microarrays aims at selecting genes whose expression is significantly different in different conditions followed by grouping them in functional categories, enabling a biological interpretation of the results. Another approach investigates the global differences in the expression of functionally related groups of genes. Here, this technique has been effective in discovering patterns related to temporal changes during infection of human cells. Another aspect explored in this thesis is related to the possibility of combining independent gene expression data for creating a catalog of genes that are selectively expressed in healthy human tissues. Not all the genes present in human cells are active; some involved in basic activities (named housekeeping genes) are expressed ubiquitously. Other genes (named tissue-selective genes) provide more specific functions and they are expressed preferably in certain cell types or tissues. Defining the tissue-selective genes is also important as these genes can cause disease with phenotype in the tissues where they are expressed. The hypothesis that gene expression could be used as a measure of the relatedness of the tissues has been also proved. Microarray experiments provide long lists of candidate genes that are often difficult to interpret and prioritize. Extending the power of microarray results is possible by inferring the relationships of genes under certain conditions. Gene transcription is constantly regulated by the coordinated binding of proteins, named transcription factors, to specific portions of the its promoter sequence. In this study, the analysis of promoters from groups of candidate genes has been utilized for predicting gene networks and highlighting modules of transcription factors playing a central role in the regulation of their transcription. Specific modules have been found regulating the expression of genes selectively expressed in the hippocampus, an area of the brain having a central role in the Major Depression Disorder. Similarly, gene networks derived from microarray results have elucidated aspects of the development of the mesencephalon, another region of the brain involved in Parkinson Disease.

FGFR1 regulated gene-expression, cell proliferation and differentiation in the developing midbrain and hindbrain

The neuroectodermal tissue close to the midbrain hindbrain boundary (MHB) is an important secondary organizer in the developing neural tube. This so-called isthmic organizer (IsO) regulates cellular survival, patterning and proliferation in the midbrain (Mb) and rhombomere 1 (R1) of the hindbrain. Signaling molecules of the IsO, such as fibroblast growth factor 8 (FGF8) and WNT1 are expressed in distinct bands of cells around the MHB. It has been previously shown that FGF-receptor 1 (FGFR1) is required for the normal development of this brain region in the mouse embryo. In the present study, we have compared the gene expression profiles of wild-type and Fgfr1 mutant embryos. We show that the loss of Fgfr1 results in the downregulation of several genes expressed close to the MHB and in the disappearance of gene expression gradients in the midbrain and R1. Our microarray screen identified several previously uncharacterized genes which may participate in the development of midbrain R1 region. Our results also show altered neurogenesis in the midbrain and R1 of the Fgfr1 mutants. Interestingly, the neuronal progenitors in midbrain and R1 show different responses to the loss of signaling through FGFR1. As Wnt1 expression at the MHB region requires the FGF signaling pathway, WNT target genes, including Drapc1, were also identified in our screen. The microarray data analysis also suggested that the cells next to the midbrain hindbrain boundary express distinct cell cycle regulators. We showed that the cells close to the border appeared to have unique features. These cells proliferate less rapidly than the surrounding cells. Unlike the cells further away from the boundary, these cells express Fgfr1 but not the other FGF receptors. The slowly proliferating boundary cells are necessary for development of the characteristic isthmic constriction. They may also contribute to compartmentalization of this brain region.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.