Predicting Transcript Production Rates in Yeast With Sparse Linear Models

To provide biological insights into transcriptional regulation, a couple of groups have recently presented models relating the promoter DNA-bound transcription factors (TFs) to downstream gene’s mean transcript level or transcript production rates over time. However, transcript production is dynamic in response to changes of TF concentrations over time. Also, TFs are not the only factors binding to promoters; other DNA binding factors (DBFs) bind as well, especially nucleosomes, resulting in competition between DBFs for binding at same genomic location. Additionally, not only TFs, but also some other elements regulate transcription. Within core promoter, various regulatory elements influence RNAPII recruitment, PIC formation, RNAPII searching for TSS, and RNAPII initiating transcription. Moreover, it is proposed that downstream from TSS, nucleosomes resist RNAPII elongation.

Here, we provide a machine learning framework to predict transcript production rates from DNA sequences. We applied this framework in the S. cerevisiae yeast for two scenarios: a) to predict the dynamic transcript production rate during the cell cycle for native promoters; b) to predict the mean transcript production rate over time for synthetic promoters. As far as we know, our framework is the first successful attempt to have a model that can predict dynamic transcript production rates from DNA sequences only: with cell cycle data set, we got Pearson correlation coefficient Cp = 0.751 and coefficient of determination r2 = 0.564 on test set for predicting dynamic transcript production rate over time. Also, for DREAM6 Gene Promoter Expression Prediction challenge, our fitted model outperformed all participant teams, best of all teams, and a model combining best team’s k-mer based sequence features and another paper’s biologically mechanistic features, in terms of all scoring metrics.

Moreover, our framework shows its capability of identifying generalizable fea- tures by interpreting the highly predictive models, and thereby provide support for associated hypothesized mechanisms about transcriptional regulation. With the learned sparse linear models, we got results supporting the following biological insights: a) TFs govern the probability of RNAPII recruitment and initiation possibly through interactions with PIC components and transcription cofactors; b) the core promoter amplifies the transcript production probably by influencing PIC formation, RNAPII recruitment, DNA melting, RNAPII searching for and selecting TSS, releasing RNAPII from general transcription factors, and thereby initiation; c) there is strong transcriptional synergy between TFs and core promoter elements; d) the regulatory elements within core promoter region are more than TATA box and nucleosome free region, suggesting the existence of still unidentified TAF-dependent and cofactor-dependent core promoter elements in yeast S. cerevisiae; e) nucleosome occupancy is helpful for representing +1 and -1 nucleosomes’ regulatory roles on transcription.

A functional variant in the Stearoyl-CoA desaturase gene promoter enhances fatty acid desaturation in pork

There is growing public concern about reducing saturated fat intake. Stearoyl-CoA desaturase (SCD) is the lipogenic enzyme responsible for the biosynthesis of oleic acid (18:1) by desaturating stearic acid (18:0). Here we describe a total of 18 mutations in the promoter and 3′ non-coding region of the pig SCD gene and provide evidence that allele T at AY487830:g.2228T>C in the promoter region enhances fat desaturation (the ratio 18:1/18:0 in muscle increases from 3.78 to 4.43 in opposite homozygotes) without affecting fat content (18:0+18:1, intramuscular fat content, and backfat thickness). No mutations that could affect the functionality of the protein were found in the coding region. First, we proved in a purebred Duroc line that the C-T-A haplotype of the 3 single nucleotide polymorphisms (SNPs) (g.2108C>T; g.2228T>C; g.2281A>G) of the promoter region was additively associated to enhanced 18:1/18:0 both in muscle and subcutaneous fat, but not in liver. We show that this association was consistent over a 10-year period of overlapping generations and, in line with these results, that the C-T-A haplotype displayed greater SCD mRNA expression in muscle. The effect of this haplotype was validated both internally, by comparing opposite homozygote siblings, and externally, by using experimental Duroc-based crossbreds. Second, the g.2281A>G and the g.2108C>T SNPs were excluded as causative mutations using new and previously published data, restricting the causality to g.2228T>C SNP, the last source of genetic variation within the haplotype. This mutation is positioned in the core sequence of several putative transcription factor binding sites, so that there are several plausible mechanisms by which allele T enhances 18:1/18:0 and, consequently, the proportion of monounsaturated to saturated fat.

Evolution of the leukotoxin promoter in genus Mannheimia

BACKGROUND: The Mannheimia species encompass a wide variety of bacterial lifestyles, including opportunistic pathogens and commensals of the ruminant respiratory tract, commensals of the ovine rumen, and pathogens of the ruminant integument. Here we present a scenario for the evolution of the leukotoxin promoter among representatives of the five species within genus Mannheimia. We also consider how the evolution of the leukotoxin operon fits with the evolution and maintenance of virulence. RESULTS: The alignment of the intergenic regions upstream of the leukotoxin genes showed significant sequence and positional conservation over a 225-bp stretch immediately proximal to the transcriptional start site of the lktC gene among all Mannheimia strains. However, in the course of the Mannheimia genome evolution, the acquisition of individual noncoding regions upstream of the conserved promoter region has occurred. The rate of evolution estimated branch by branch suggests that the conserved promoter may be affected to different extents by the types of natural selection that potentially operate in regulatory regions. Tandem repeats upstream of the core promoter were confined to M. haemolytica with a strong association between the sequence of the repeat units, the number of repeat units per promoter, and the phylogenetic history of this species. CONCLUSION: The mode of evolution of the intergenic regions upstream of the leukotoxin genes appears to be highly dependent on the lifestyle of the bacterium. Transition from avirulence to virulence has occurred at least once in M. haemolytica with some evolutionary success of bovine serotype A1/A6 strains. Our analysis suggests that changes in cis-regulatory systems have contributed to the derived virulence phenotype by allowing phase-variable expression of the leukotoxin protein. We propose models for how phase shifting and the associated virulence could facilitate transmission to the nasopharynx of new hosts.

RAG2 is regulated differentially in B and T cells by elements 5′ of the promoter

To study RAG2 gene regulation in vivo, we developed a blastocyst complementation method in which RAG2-deficient embryonic stem cells were transfected with genomic clones containing RAG2 and then assessed for their ability to generate lymphocytes. A RAG2 genomic clone that contained only the RAG2 promoter sequences rescued V(D)J recombination in RAG2-deficient pro-B cell lines, but did not rescue development of RAG2-deficient lymphocytes in vivo. However, inclusion of varying lengths of sequences 5′ of the RAG2 promoter generated constructs capable of rescuing only in vivo B cell development, as well as other constructs that rescued both B and T cell development. In particular, the 2-kb 5′ region starting just upstream of the RAG2 promoter, as well as the region from 2–7 kb 5′, could independently drive B cell development, but not efficient T cell development. Deletion of the 2-kb 5′ region from the murine germ line demonstrated that this region was not required for RAG expression sufficient to generate normal B or T cell numbers, implying redundancy among 5′ elements. We conclude that RAG2 expression in vivo requires elements beyond the core promoter, that such elements contribute to differential regulation in the B vs. T lineages, and that sequences sufficient to direct B cell expression are located in the promoter-proximal 5′ region.

Nutritional regulation of the fatty acid synthase promoter in vivo: Sterol regulatory element binding protein functions through an upstream region containing a sterol regulatory element

The transcription of fatty acid synthase (FAS), a central enzyme in de novo lipogenesis, is dramatically induced by fasting/refeeding and insulin. We reported that upstream stimulatory factor binding to the −65 E-box is required for induction of the FAS transcription by insulin in 3T3-L1 adipocytes. On the other hand, we recently found that two upstream 5′ regions are required for induction in vivo by fasting/refeeding and insulin; one at −278 to −131 albeit at a low level, and the other at −444 to −278 with an E-box at −332 where upstream stimulatory factor functions for maximal induction. Here, we generated double transgenic mice carrying the chloramphenicol acetyltransferase reporter driven by the various 5′ deletions of the FAS promoter region and a truncated active form of the sterol regulatory element (SRE) binding protein (SREBP)-1a. We found that SREBP participates in the nutritional regulation of the FAS promoter and that the region between −278 and −131 bp is required for SREBP function. We demonstrate that SREBP binds the −150 canonical SRE present between −278 and −131, and SREBP can function through the −150 SRE in cultured cells. These in vivo and in vitro results indicate that SREBP is involved in the nutritional induction of the FAS promoter via the −278/−131 region and that the −150 SRE is the target sequence.

Circadian-Regulated Transcription of the psbD Light-Responsive Promoter in Wheat Chloroplasts1

The level of mRNAs derived from the plastid-encoded psbD light-responsive promoter (LRP) is controlled by a circadian clock(s) in wheat (Triticum aestivum). The circadian oscillations in the psbD LRP mRNA level persisted for at least three cycles in continuous light and for one cycle in continuous dark, with maxima in subjective morning and minima in subjective early night. In vitro transcription in chloroplast extracts revealed that the circadian cycles in the psbD LRP mRNA level were dominantly attributed to the circadian-regulated transcription of the psbD LRP. The effects of various mutations introduced into the promoter region on the psbD LRP activity in vitro suggest the existence of two positive elements located between −54 and −36, which generally enhance the transcription activity, and an anomalous core promoter structure lacking the functional “−35” element, which plays a crucial role in the circadian fluctuation and light dependency of psbD LRP transcription activity.

Genetic polymorphisms in TLR4, CR1 and Duffy genes are not associated with malaria resistance in patients from Baixo Amazonas region, Brazil

The main purpose of this research was to analyze the relation of the genetic polymorphisms frequently expressed by antigen-presenting cells, erythrocytes and malaria susceptibility/resistance with the human malaria infection cases. The sample used consisted of 23 Plasmodium vivax ( Pv)- and P. falciparum ( Pf)-infected patients, and 21 healthy individuals as a control group, from the Baixo Amazonas population in Para, Brazil. The Asp299Gly polymorphisms in the Toll-like receptor 4 ( TLR4), and Gly42Asp, Arg89Cys, Ala100Thr, and T-33C in the Duffy gene ( FY) were analyzed by restriction fragment length polymorphism-polymerase chain reaction. The Lys1590Glu and Arg1601Gly polymorphisms in the complement receptor type 1 (CR1) were analyzed by DNA sequencing. According to the results obtained and statistical analysis considering a significance level or alpha = 0.01, we conclude that the low heterozygote frequency (2.27%) for the Asp299Gly mutation, detected in the TLR4 gene, is not related to the Pv and Pf infections in the patients analyzed. Also, the promoter region GATA-1 analysis of the FY gene in the Pv-infected patients showed that the heterozygote frequency for the T-33C mutation (11.36% of the infected patients and 20.45% of the control patients) is not related to infection resistance. Regarding the CR1 gene, the observed heterozygote frequency (9.09%) for the Arg1601Gly mutation in Pf-infected patients when compared to heterozygote frequency in the control group (18.18%) suggests that there is no correlation with infection resistance.

The effects of dietary supplementation of methionine on genomic stability and p53 gene promoter methylation in rats

Methionine is a component of one-carbon metabolism and a precursor of S-adenosylmethionine (SAM), the methyl donor for DNA methylation. When methionine intake is high, an increase of S-adenosylmethionine (SAM) is expected. DNA methyltransferases convert SAM to S-adenosylhomocysteine (SAH). A high intracellular SAH concentration could inhibit the activity of DNA methyltransferases. Therefore, high methionine ingestion could induce DNA damage and change the methylation pattern of tumor suppressor genes. This study investigated the genotoxicity of a methionine-supplemented diet. It also investigated the diet`s effects on glutathione levels, SAM and SAH concentrations and the gene methylation pattern of p53. Wistar rats received either a methionine-supplemented diet (2% methionine) or a control diet (0.3% methionine) for six weeks. The methionine-supplemented diet was neither genotoxic nor antigenotoxic to kidney cells, as assessed by the comet assay. However, the methionine-supplemented diet restored the renal glutathione depletion induced by doxorubicin. This fact may be explained by the transsulfuration pathway, which converts methionine to glutathione in the kidney. Methionine supplementation increased the renal concentration of SAH without changing the SAM/SAH ratio. This unchanged profile was also observed for DNA methylation at the promoter region of the p53 gene. Further studies are necessary to elucidate this diet`s effects on genomic stability and DNA methylation. (C) 2011 Elsevier ay. All rights reserved.

Promoter studies of a polydnavirus gene from Cotesia rubecula (Hym : Braconidae)

The Cotesia rubecula polydnavirus gene, CrV1, is expressed in a highly transient fashion. Within four hours after egg deposition and virus infection, tissues of the host caterpillar, Pieris rapae, express high levels of the transcript. Twelve hours after infection no transcripts are visible. We have previously shown that the CrV1 secreted protein is mainly produced in host haemocytes. In haemocytes, immune functions such as phagocytosis and cell spreading are abolished by destabilization of the cell cytoskeleton. To test whether the observed down-regulation of CrV1 transcripts is mediated by transcriptional control or by other factors, such as the disruption of cytoskeleton in CrV1-inactivated cells, we cloned the promoter and the 3' untranslated region of the CrV1 gene to study CrV1 expression. The promoter region of the CrV1 gene was cloned into baculovirus expression systems along with the CAT reporter gene. Molecular analyses showed that the CAT gene under the control of CrV1 promoter is expressed as early as 2 h post infection and continues until late phase of infection suggesting that down-regulation of CrV1 expression in host haemocytes is perhaps mediated by post-transcriptional mechanisms.

High Frequency of Lamivudine Resistance Mutations in Brazilian Patients Co-Infected With HIV and Hepatitis B

This study analyzed the genotype distribution and frequency of lamivudine (LAM) and tenofovir (TDF) resistance mutations in a group of patients co-infected with HIV and hepatitis B virus (HBV). A cross-sectional study of 847 patients with HIV was conducted. Patients provided blood samples for HBsAg detection. The load of HBV was determined using an ""in-house"" real-time polymerase chain reaction. HBV genotypes/subgenotypes, antiviral resistance, basal core promoter (BCP), and precore mutations were detected by DNA sequencing. Twenty-eight patients with co-infection were identified. The distribution of HBV genotypes among these patients was A (n = 9; 50%), D (n = 4; 22.2%), G (n = 3; 16.7%), and F (n = 2; 11.1%). Eighteen patients were treated with LAM and six patients were treated with LAM plus TDF. The length of exposure to LAM and TDF varied from 4 to 216 months. LAM resistance substitutions (rtL180M + rtM204V) were detected in 10 (50%) of the 20 patients with viremia. This pattern and an accompanying rtV173L mutation was found in four patients. Three patients with the triple polymerase substitution pattern (rtV173L+ rtL180M + rtM204V) had associated changes in the envelope gene (sE164D + sl195M). Mutations in the BCP region (A1762T, G1764A) and in the precore region (G1896A, G1899A) were also found. No putative TDF resistance substitution was detected. The data suggest that prolonged LAM use is associated with the emergence of particular changes in the HBV genome, including substitutions that may elicit a vaccine escape phenotype. No putative TDF resistance change was detected after prolonged use of TDF. J. Med. Virol. 82:1481-1488, 2010. (C) 2010 Wiley-Liss, Inc.

Silencing of LRRC49 and THAP10 genes by bidirectional promoter hypermethylation is a frequent event in breast cancer

Previously we found that levels of LRRC49 (leucine rich repeat containing 49; FLJ20156) transcripts were elevated in ER-positive breast tumors compared with ER-negative breast tumors. The LRRC49 gene is located on chromosome 15q23 in close proximity to the THAP10 (THAP domain containing 10) gene. These two genes have a bidirectional organization being arranged head-to-head on opposite strands, possibly sharing the same promoter region. Analysis of the promoter region of this gene pair revealed the presence of potential estrogen response elements (EREs), suggesting the potential of this promoter to be under the control of estrogen. We used quantitative real-time PCR (qPCR) to evaluate the expression of LRRC49 and THAP10 in a series of 72 primary breast tumors, and found reduced LRRC49 and THAP10 expression in 61 and 46% of the primary breast tumors analyzed, respectively. In addition, the occurrence of LRRC49/THAP10 promoter hypermethylation was examined by methylation specific PCR (MSP) in a sub-group of the breast tumors. Hypermethylation was observed in 57.5% of the breast tumors analyzed, and the levels of mRNA expression of both genes were inversely correlated with promoter hypermethylation. We investigated the effects of 17 beta-estradiol on LRRC49 and THAP10 expression in MCF-7 breast cancer cells and found both transcripts to be up-regulated 2- to 3-fold upon 17 beta-estradiol treatment. Our results show that the transcripts of LRRC49/THAP10 bidirectional gene pair are co-regulated by estrogen and that hypermethylation of the bidirectional promoter region simultaneously silences both genes. Further studies will be necessary to elucidate the role of LRRC49/THAP10 down-regulation in breast cancer.

DC-SIGN (CD209) gene promoter polymorphisms in a Brazilian population and their association with human T-cell lymphotropic virus type 1 infection

This study evaluated four polymorphisms located in the DC-SIGN (CD209) gene promoter region (positions -336, -332 -201 and -139) in DNA samples from four Brazilian ethnic groups (Caucasians, Afro-Brazilian, Asians and Amerindians) to establish the population distribution of these single-nucleotide polymorphisms (SNPs) and correlated DC-SIGN polymorphisms and infection in samples from human T-cell lymphotropic virus type 1 (HTLV-1)-infected individuals. To identify CD209 SNPs, 452 bp of the CD209 promoter region were sequenced and the genotype and allelic frequencies were evaluated. This is the first study to show genetic polymorphism in the CD209 gene in distinct Brazilian ethnic groups with the distribution of allelic and genotypic frequency. The results showed that -336A and -139A SNPs were quite common in Asians and that the -201T allele was not observed in Caucasians, Asians or Amerindians. No significant differences were observed between individuals with HTLV-1 disease and asymptomatic patients. However, the -336A variant was more frequent in HTLV-1 -infected patients [HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP), 80%; healthy asymptomatic HTLV-1 carriers, 90 %] than in the control group (70 %) [P=0.0197, odds ratio (OR)=2.511, 95 % confidence interval (CI)=1.218-5.179). In addition, the -139A allele was found to be associated with protection against HTLV-1 infection (P=0.0037, OR=0.3758, 95% CI=0.1954-0.7229) when the HTLV-1 -infected patients as a whole were compared with the healthy-control group. These observations suggest that the -139A allele may be associated with HTLV-1 infection, although no significant association was observed among asymptomatic and HAM/TSP patients. In conclusion, the variation observed in SNPs -336 and -139 indicates that this lectin may be of crucial importance in the susceptibility/transmission of HTLV-1 infections.

Absence of the HLA-G*0113N allele in Amerindian populations from the Brazilian Amazon region

The HLA-G gene is predominantly expressed at the maternal-fetal interface and has been associated with maternal-fetal tolerance. The HLA-G*0113N is a null allele defined by the insertion of a premature stop codon at exon 2, observed in a single Ghanaian individual. Likewise the G*0105N allele, the occurrence of the HL4-G*0113N in a population from an area with high pathogen load suggests that the reduced HLA-G expression in G*0113N heterozygous placentas could improve the intrauterine defense against infections. The presence of the G*0113N allele here was investigated in 150 Amerindians from five isolated tribes that inhabit the Central Amazon and in 295 admixed individuals from the State of Sao Paulo, Southeastern Brazil, previously genotyped for HLA-G. No copy of the G*0113N null allele was found in both population samples by exon 2 sequence-based analysis, reinforcing its restricted occurrence in Africa. (C) 2010 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.

The-202 A Allele of Insulin-Like Growth Factor Binding Protein-3 (IGFBP3) Promoter Polymorphism Is Associated with Higher IGFBP-3 Serum Levels and Better Growth Response to Growth Hormone Treatment in Patients with Severe Growth Hormone Deficiency

Context: Genetic factors that influence the response to recombinant human GH (rhGH) therapy remain mostly unknown. To date, only the GH receptor gene has been investigated. Objective: The aim of the study was to assess the influence of a polymorphism in the IGF-binding protein-3 (IGFBP-3) promoter region (-202 A/C) on circulating IGFBP-3 levels and growth response to rhGH therapy in children with GH deficiency (GHD). Design and Patients: -202 A/C IGFBP3 genotyping (rs2854744) was correlated with data of 71 children with severe GHD who remained prepubertal during the first year of rhGH treatment. Main Outcome Measures: We measured IGFBP-3 levels and first year growth velocity (GV) during rhGH treatment. Results: Clinical and laboratory data at the start of treatment were indistinguishable among patients with different -202 A/C IGFBP3 genotypes. Despite similar rhGH doses, patients homozygous for the A allele presented higher IGFBP-3 SD score levels and higher mean GV in the first year of rhGH treatment than patients with AC or CC genotypes (first year GV, AA = 13.0 +/- 2.1 cm/yr, AC = 11.4 +/- 2.5 cm/yr, and CC = 10.8 +/- 1.9 cm/yr; P = 0.016). Multiple linear regression analyses demonstrated that the influence of -202 A/C IGFBP3 genotype on IGFBP-3 levels and GV during the first year of rhGH treatment was independent of other variables. Conclusion: The -202 A allele of IGFBP3 promoter region is associated with increased IGFBP-3 levels and GV during rhGH treatment in prepubertal GHD children. (J Clin Endocrinol Metab 94: 588-595, 2009)

The genetic structure of 3 ` untranslated region of the HLA-G gene: polymorphisms and haplotypes

The HLA-G gene is predominantly expressed at the maternal-fetal interface. It has been associated with maternal-fetal tolerance and in the inhibition of cytotoxic T lymphocyte and natural killer cytolytic functions. At least two variations in the 3` untranslated region (UTR) of HLA-G locus are associated with HLA-G expression levels, the 14-bp deletion/insertion polymorphism and the +3142 single-nucleotide polymorphism (SNP). However, this region has not been completely characterized yet. The variability of the 3`UTR of HLA-G gene and its haplotype structure were characterized in 155 individuals from Brazil, as well as HLA-G alleles associated with each of the 3`UTR haplotype. The following eight variation sites were detected: the 14-bp polymorphism and SNPs at the positions +3003T/C, +3010C/G, +3027A/C, +3035C/T, +3142G/C, +3187A/G and +3196C/G. Similarly, 11 different 3`UTR haplotypes were identified and several HLA-G alleles presented only one 3`UTR haplotype. In addition, a high linkage disequilibrium among the variation sites was detected, especially among the 14-bp insertion and the alleles +3142G and +3187A, all previously associated with low mRNA availability, demonstrating that their effects are not independent. The detailed analyses of 3`UTR of the HLA-G locus may shed some light into mechanisms underlying the regulation of HLA-G expression. Genes and Immunity (2010) 11, 134-141; doi: 10.1038/gene.2009.74; published online 1 October 2009