Transcriptome analysis of a respiratory Saccharomyces cerevisiae strain suggests the expression of its phenotype is glucose insensitive and predominantly controlled by Hap4, Cat8 and Mig1

BACKGROUND: We previously described the first respiratory Saccharomyces cerevisiae strain, KOY.TM6*P, by integrating the gene encoding a chimeric hexose transporter, Tm6*, into the genome of an hxt null yeast. Subsequently we transferred this respiratory phenotype in the presence of up to 50 g/L glucose to a yeast strain, V5 hxt1-7Delta, in which only HXT1-7 had been deleted. In this study, we compared the transcriptome of the resultant strain, V5.TM6*P, with that of its wild-type parent, V5, at different glucose concentrations. RESULTS: cDNA array analyses revealed that alterations in gene expression that occur when transitioning from a respiro-fermentative (V5) to a respiratory (V5.TM6*P) strain, are very similar to those in cells undergoing a diauxic shift. We also undertook an analysis of transcription factor binding sites in our dataset by examining previously-published biological data for Hap4 (in complex with Hap2, 3, 5), Cat8 and Mig1, and used this in combination with verified binding consensus sequences to identify genes likely to be regulated by one or more of these. Of the induced genes in our dataset, 77% had binding sites for the Hap complex, with 72% having at least two. In addition, 13% were found to have a binding site for Cat8 and 21% had a binding site for Mig1. Unexpectedly, both the up- and down-regulation of many of the genes in our dataset had a clear glucose dependence in the parent V5 strain that was not present in V5.TM6*P. This indicates that the relief of glucose repression is already operable at much higher glucose concentrations than is widely accepted and suggests that glucose sensing might occur inside the cell. CONCLUSION: Our dataset gives a remarkably complete view of the involvement of genes in the TCA cycle, glyoxylate cycle and respiratory chain in the expression of the phenotype of V5.TM6*P. Furthermore, 88% of the transcriptional response of the induced genes in our dataset can be related to the potential activities of just three proteins: Hap4, Cat8 and Mig1. Overall, our data support genetic remodelling in V5.TM6*P consistent with a respiratory metabolism which is insensitive to external glucose concentrations.

A Case-Only Genome-wide Association Study of Gender- and Age-specific Risk Markers for Childhood Leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.

Analysis of expression quantitative trait loci in D. melanogaster

Expression Quantitative Trait Loci (eQTL) analysis allows for the identification of genetic variation associated with variation in gene expression. It is often unclear however, which of the associated variants are causal, and by what mechanism. Integrating functional genomic data with eQTL data can provide insight into the impact of natural variation in the population, and the nature of the transcriptional machinery itself. In this thesis, I integrate functional genomic data with eQTL data derived from both 5’ CAGE and 3’ TagXseq expression assays, in developing embryos. I first use both datasets to analyse the transcription landscape in embryonic D., melanogaster, and then carry out an analysis of sequence motifs associated with transcription factor binding sites, promoters, and 3’ polyadenylation sites. Finally, I integrate functional genomic data, including these novel sequence motifs, to shed light on the mechanisms of gene expression variation in D.,melanogaster. I am able to demonstrate that some variants effecting gene regulation in Drosophila are found within haplotypes which buffer their effects.

Flexible strategies for association analysis with genomic phenotypes

Thesis (Ph.D.)--University of Washington, 2016-08

Investigating the role of glutathione and glutathione biosynthetic genes in the adaptation of Anopheles gambiae to insecticides

Malaria remains a serious public health challenge in the tropical world, with 584,000 deaths globally in 2013, of which 90% occurred in Africa, and mostly in pregnant women and children under the age of five. Anopheles gambiae (An. gambiae) is the principal malaria vector in Africa, where vector control measures involve the use of insecticides in the forms of long-lasting insecticide-treated nets (LLINs) and indoor residual spraying (IRS). The development of insecticides resistance mitigates these approaches. Glutathione (GSH) is widely distributed among all living organisms, and is associated with detoxification pathways, especially the Glutathione S-transferases (GSTs). Its direct involvement and relevance in insecticide resistance in An. gambiae has not been determined. Thus, this work examines the contribution of GSH, its biosynthetic genes (GCLM, GCLC) and their possible transcriptional regulator Nrf2 in insecticide resistance in An. gambiae sampled from agricultural setting (areas of intensive agriculture) and residential setting (domestic area). Bioinformatics analysis, W.H.O. adult susceptibility bioassays and molecular techniques were employed to investigate. Total RNA was first isolated from the adults An. gambiae mosquitoes raised from agricultural and residential field-caught larvae which had been either challenged or unchallenged with insecticides. Semi-quantitative RT-PCR using gel image densitometry was used to determine the expression levels of GCLM, GCLC genes and Nrf2. Bioinformatics’ results established the presence of putative AGAP010259 (AhR) and AGAP005300 (Nf2e1) transcription factor binding sites in An. gambiae GCLC and GCLM promoters in silico. An. gambiae s.l. studied here were highly resistant to DDT and permethrin but less resistant to bendiocarb. Both knockdown resistance (kdr) mutation variants L1014S and L1014F that confers resistance to pyrethroid insecticides were identified in both An. coluzzii and An. arabiensis sampled from northern Nigeria. The L1014F was much associated with An. coluzzii. A significant positive correlation (P=0.04) between the frequency of the L1014F point mutation and resistance to DDT and permethrin was observed. However, a weak or non-significant correlation (P=0.772) between the frequency of the L1014S point mutation and resistance was also found. L1014S and L1014F mutations co-occurred in both agricultural and residential settings with high frequencies. However, the frequencies of the two mutations were greater in the agricultural settings than in the residential settings. The levels of total, reduced and oxidized GSH were significantly higher in mosquitoes from agricultural sites than those from residential sites. Increased oxidized GSH levels appears to correlate with higher DDT resistance. The expression levels of GCLM, GCLC and Nrf2 were also significantly up-regulated in adults An. gambiae raised from agricultural and residential field-caught larvae when challenged with insecticide. However, there was higher constitutive expression of GCLM, GCLC and Nrf2 in mosquitoes from agricultural setting. The increased expression levels of these genes and also GSH levels in this population suggest their roles in the response and adaptation of An. gambiae to insecticide challenges. There exists the feasibility of using GSH status in An. gambiae to monitor adaptation and resistance to insecticides.

A case-only genome-wide association study of gender- and age-specific risk markers for childhood leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. ^ Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5 ) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. ^ The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.^

Identification of the N-Linked Glycosylation Sites of the Transcription Factor Rest and Effect of Glycosylation on DNA Binding and Transcriptional Activity

REST is a zinc-finger transcription factor implicated in several processes such as maintenance of embryonic stem cell pluripotency and regulation of mitotic fidelity in non-neuronal cells [Chong et al., 1995]. The gene encodes for a 116-kDa protein that acts as a molecular platform for co-repressors recruitment and promotes modifications of DNA and histones [Ballas, 2005]. REST showed different apparent molecular weights, consistent with the possible presence of post-translational modifications [Lee et al., 2000]. Among these the most common is glycosylation, the covalent attachment of carbohydrates during or after protein synthesis [Apweiler et al., 1999] My thesis has ascertained, for the first time, the presence of glycan chians in the transcription factor REST. Through enzymatic deglycosylation and MS, oligosaccharide composition of glycan chains was evaluated: a complex mixture of glycans, composed of N-acetylgalactosamine, galactose and mannose, was observed thus confirming the presence of O- and N-linked glycan chains. Glycosylation site mapping was done using a 18O-labeling method and MS/MS and twelve potential N-glycosylation sites were identified. The most probable glycosylation target residues were mutated through site-directed mutagenesis and REST mutants were expressed in different cell lines. Variations in the protein molecular weight and mutant REST ability to bind the RE-1 sequence were analyzed. Gene reporter assays showed that, altogether, removal of N-linked glycan chains causes loss of transcriptional repressor function, except for mutant N59 which showed a slight residual repressor activity in presence of IGF-I. Taken togheter these results demonstrate the presence of complex glycan chians in the transcription factor REST: I have depicted their composition, started defining their position on the protein backbone and identified their possible role in the transcription factor functioning. Considering the crucial role of glycosylation and transcription factors activity in the aetiology of many diseases, any further knowledge could find important and interesting pharmacological application.

Identification of novel functional TBP-binding sites and general factor repertoires.

Our current knowledge of the general factor requirement in transcription by the three mammalian RNA polymerases is based on a small number of model promoters. Here, we present a comprehensive chromatin immunoprecipitation (ChIP)-on-chip analysis for 28 transcription factors on a large set of known and novel TATA-binding protein (TBP)-binding sites experimentally identified via ChIP cloning. A large fraction of identified TBP-binding sites is located in introns or lacks a gene/mRNA annotation and is found to direct transcription. Integrated analysis of the ChIP-on-chip data and functional studies revealed that TAF12 hitherto regarded as RNA polymerase II (RNAP II)-specific was found to be also involved in RNAP I transcription. Distinct profiles for general transcription factors and TAF-containing complexes were uncovered for RNAP II promoters located in CpG and non-CpG islands suggesting distinct transcription initiation pathways. Our study broadens the spectrum of general transcription factor function and uncovers a plethora of novel, functional TBP-binding sites in the human genome.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Artificial zinc finger fusions targeting Sp1-binding sites and the trans-activator-responsive element potently repress transcription and replication of HIV-1.

Tat activates transcription by interacting with Sp1, NF-kappaB, positive transcription elongation factor b, and trans-activator-responsive element (TAR). Tat and Sp1 play major roles in transcription by protein-protein interactions at human immunodeficiency virus, type 1 (HIV-1) long terminal repeat. Sp1 activates transcription by interacting with cyclin T1 in the absence of Tat. To disrupt the transcription activation by Tat and Sp1, we fused Sp1-inhibiting polypeptides, zinc finger polypeptide, and the TAR-binding mutant Tat (TatdMt) together. A designed or natural zinc finger and Tat mutant fusion was used to target the fusion to the key regulatory sites (GC box and TAR) on the long terminal repeat and nascent short transcripts to disrupt the molecular interaction that normally result in robust transcription. The designed zinc finger and TatdMt fusions were targeted to the TAR, and they potently repressed both transcription and replication of HIV-1. The Sp1-inhibiting POZ domain, TatdMt, and zinc fingers are key functional domains important in repression of transcription and replication. The designed artificial zinc fingers were targeted to the high affinity Sp1-binding site, and by being fused with TatdMt and POZ domain, they strongly block both Sp1-cyclin T1-dependent transcription and Tat-dependent transcription, even in the presence of excess expressed Tat.

Phosphorylation of the rRNA transcription factor upstream binding factor promotes its association with TATA binding protein

rRNA synthesis by RNA polymerase I requires both the promoter selectivity factor 1, which is composed of TATA binding protein (TBP) and three TBP-associated factors, and the activator upstream binding factor (UBF). Whereas there is strong evidence implicating a role for phosphorylation of UBF in the control of growth-induced increases in rRNA transcription, the mechanism of this effect is not known. Results of immunoprecipitation studies with TBP antibodies showed increased recovery of phosphorylated UBF from growth-stimulated smooth muscle cells. Moreover, using an immobilized protein-binding assay, we found that phosphorylation of UBF in vivo in response to stimulation with different growth factors or in vitro with smooth muscle cell nuclear extract increased its binding to TBP. Finally, we demonstrated that UBF–TBP binding depended on the C-terminal ‘acidic tail’ of UBF that was hyperphosphorylated at multiple serine sites after growth factor stimulation. Results of these studies suggest that phosphorylation of UBF and subsequent binding to TBP represent a key regulatory step in control of growth-induced increases in rRNA synthesis.

E2F4-RB and E2F4-p107 complexes suppress gene expression by transforming growth factor β through E2F binding sites

Transforming growth factor β (TGF-β) causes growth arrest in most cell types. TGF-β induces hypophosphorylation of retinoblastoma susceptibility gene 1 product (RB), which sequesters E2F factors needed for progression into S phase of the cell cycle, thereby leading to cell cycle arrest at G1. It is possible, however, that the E2F-RB complex induced by TGF-β may bind to E2F sites and suppress expression of specific genes whose promoters contain E2F binding sites. We show here that TGF-β treatment of HaCaT cells induced the formation of E2F4-RB and E2F4-p107 complexes, which are capable of binding to E2F sites. Disruption of their binding to DNA with mutation in the E2F sites did not change the expression from promoters of E2F1, B-myb, or HsORC1 genes in cycling HaCaT cells. However, the same mutation stimulated 5- to 6-fold higher expression from all three promoters in cells treated with TGF-β. These results suggest that E2F binding sites play an essential role in the transcription repression of these genes under TGF-β treatment. Consistent with their repression of TGF-β-induced gene expression, introduction of E2F sites into the promoter of cyclin-dependent kinase inhibitor p15INK4B gene effectively inhibited its induction by TGF-β. Experiments utilizing Gal4-RB and Gal4-p107 chimeric constructs demonstrated that either RB or p107 could directly repress TGF-β induction of p15INK4B gene when tethered to p15INK4B promoter through Gal4 DNA binding sites. Therefore, E2F functions to bring RB and p107 to E2F sites and represses gene expression by TGF-β. These results define a specific function for E2F4-RB and E2F4-p107 complexes in gene repression under TGF-β treatment, which may constitute an integral part of the TGF-β-induced growth arrest program.

Domains of transcription factor Sp1 required for synergistic activation with sterol regulatory element binding protein 1 of low density lipoprotein receptor promoter.

Feedback regulation of transcription from the low density lipoprotein (LDL) receptor gene is fundamentally important in the maintenance of intracellular sterol balance. The region of the LDL receptor promoter responsible for normal sterol regulation contains adjacent binding sites for the ubiquitous transcription factor Sp1 and the cholesterol-sensitive sterol regulatory element-binding proteins (SREBPs). Interestingly, both are essential for normal sterolmediated regulation of the promoter. The cooperation by Sp1 and SREBP-1 occurs at two steps in the activation process. SREBP-1 stimulates the binding of Sp1 to its adjacent recognition site in the promoter followed by enhanced stimulation of transcription after both proteins are bound to DNA. In the present report, we have defined the protein domains of Sp1 that are required for both synergistic DNA binding and transcriptional activation. The major activation domains of Sp1 that have previously been shown to be essential to activation of promoters containing multiple Sp1 sites are required for activation of the LDL receptor promoter. Additionally, the C domain is also crucial. This slightly acidic approximately 120-amino acid region is not required for efficient synergistic activation by multiple Sp1 sites or in combination with other recently characterized transcriptional regulators. We also show that Sp1 domain C is essential for full, enhanced DNA binding by SREBP-1. Taken together with other recent studies on the role of Sp1 in promoter activation, the current experiments suggest a unique combinatorial mechanism for promoter activation by two distinct transcription factors that are both essential to intracellular cholesterol homeostasis.

neptune, a Kruppel-like transcription factor that participates in primitive erythropoiesis in Xenopus

The specification of the erythroid lineage from hematopoietic stem cells requires the expression and activity of lineage-specific transcription factors. One transcription factor family that has several members involved in hematopoiesis is the Kruppel-like factor (KLF) family [1]. For example, erythroid KLF (EKLF) regulates beta -globin expression during erythroid differentiation [2-6]. KLFs share a highly conserved zinc finger-based DNA binding domain (DBD) that mediates binding to CACCC-box and GC-rich sites, both of which are frequently found in the promoters of hematopoietic genes. Here, we identified a novel Xenopus KLF gene, neptune, which is highly expressed in the ventral blood island (VBI), cranial ganglia, and hatching and cement glands. neptune expression is induced in response to components of the BMP-4 signaling pathway in injected animal cap explants. Similar to its family member, EKLF, Neptune can bind CACCC-box and GC-rich DNA elements. We show that Neptune cooperates with the hematopoietic transcription factor XGATA-1 to enhance globin induction in animal cap explants. A fusion protein comprised of Neptune's DBD and the Drosophila engrailed repressor domain suppresses the induction of globin in ventral marginal zones and in animal caps. These studies demonstrate that Neptune is a positive regulator of primitive erythropoiesis in Xenopus.