Hepatocyte nuclear factor-4 alpha regulates the human apolipoprotein AV gene: Identification of a novel response element and involvement in the control by peroxisome proliferator-activated receptor-gamma coactivator-1 alpha, AMP-activated protein kinase, and mitogen-activated protein kinase pathway

The recently discovered apolipoprotein AV (apoAV) gene has been reported to be a key player in modulating plasma triglyceride levels. Here we identify the hepatocyte nuclear factor-4 (HNF-4 ) as a novel regulator of human apoAV gene. Inhibition of HNF-4 expression by small interfering RNA resulted in down-regulation of apoAV. Deletion, mutagenesis, and binding assays revealed that HNF-4 directly regulates human apoAV promoter through DR1 [a direct repeat separated by one nucleotide (nt)], and via a novel element for HNF-4 consisting of an inverted repeat separated by 8 nt (IR8). In addition, we show that the coactivator peroxisome proliferator-activated receptor- coactivator-1 was capable of stimulating the HNF-4 -dependent transactivation of apoAV promoter. Furthermore, analyses in human hepatic cells demonstrated that AMP-activated protein kinase (AMPK) and the MAPK signaling pathway regulate human apoAV expression and suggested that this regulation may be mediated, at least in part, by changes in HNF-4 . Intriguingly, EMSAs and mice with a liver-specific disruption of the HNF-4 gene revealed a species-distinct regulation of apoAV by HNF-4 , which resembles that of a subset of HNF-4 target genes. Taken together, our data provide new insights into the binding properties and the modulation of HNF-4 and underscore the role of HNF-4 in regulating triglyceride metabolism.

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

Association of NTRK3 and its interaction with NGF suggest an altered cross-regulation of the neurotrophin signaling pathway in eating disorders

Eating disorders (EDs) are complex psychiatric diseases that include anorexia nervosa and bulimia nervosa, and have higher than 50% heritability. Previous studies have found association of BDNF and NTRK2 to ED, while animal models suggest that other neurotrophin genes might also be involved in eating behavior. We have performed a family-based association study with 151 TagSNPs covering 10 neurotrophin signaling genes: NGFB, BDNF, NTRK1, NGFR/p75, NTF4/5, NTRK2, NTF3, NTRK3, CNTF and CNTFR in 371 ED trios of Spanish, French and German origin. Besides several nominal associations, we found a strong significant association after correcting for multiple testing (P = 1.04 × 10−4) between ED and rs7180942, located in the NTRK3 gene, which followed an overdominant model of inheritance. Interestingly, HapMap unrelated individuals carrying the rs7180942 risk genotypes for ED showed higher levels of expression of NTRK3 in lymphoblastoid cell lines. Furthermore, higher expression of the orthologous murine Ntrk3 gene was also detected in the hypothalamus of the anx/anx mouse model of anorexia. Finally, variants in NGFB gene appear to modify the risk conferred by the NTRK3 rs7180942 risk genotypes (P = 4.0 × 10−5) showing a synergistic epistatic interaction. The reported data, in addition to the previous reported findings for BDNF and NTRK2, point neurotrophin signaling genes as key regulators of eating behavior and their altered cross-regulation as susceptibility factors for EDs.

Polymorphisms of pyrimidine pathway enzymes encoding genes and HLA-B*4001 carriage in stavudine-associated lipodystrophy in HIV-infected patients.

: To assess in a cohort of Caucasian patients exposed to stavudine (d4T) the association of polymorphisms in pyrimidine pathway enzymes and HLA-B*4001 carriage with HIV lipodystrophy syndrome (HALS). 336 patients, 187 with HALS and 149 without HALS, and 72 controls were recruited. HALS was associated with the presence of a low expression, thymidylate synthase (TS) genotype polymorphism. Methylene-tetrahydrofolate reductase (MTHFR) gene polymorphisms and HLA-B*4001 carriage were not associated with HALS or d4T-TP intracellular levels. In conclusion HALS is associated with combined low-expression TS and MTHFR associated with high activity polymorphisms but not with HLA-B*4001 carriage.

Fsp27/CIDEC is a CREB target gene induced during early fasting in liver and regulated by FA oxidation rate

FSP27 (CIDEC in humans) is a protein associated with lipid droplets that downregulates the fatty acid oxidation (FAO) rate when it is overexpressed. However, little is known about its physiological role in liver. Here, we show that fasting regulates liver expression of Fsp27 in a time-dependent manner. Thus, during the initial stages of fasting a maximal induction of 800-fold was achieved, while during the later phase of fasting, Fsp27 expression decreased. The early response to fasting can be explained by a canonical PKA-CREB-CRTC2 signaling pathway since: i) CIDEC expression was induced by forskolin, ii) Fsp27 promoter activity was increased by CREB, and iii) Fsp27 expression was upregulated in the liver of Sirt1 knockout animals. Interestingly, pharmacological (etomoxir) or genetic (Hmgcs2 interference) inhibition of the FAO rate increases the in vivo expression of Fsp27 during fasting. Similarly, CIDEC expression was upregulated in HepG2 cells by either etomoxir or HMGCS2 interference. Our data indicate that there is a kinetic mechanism of auto-regulation between short- and long-term fasting, by which free fatty acids delivered to the liver during early fasting are accumulated/exported by FSP27/CIDEC, while over longer periods of fasting they are degraded in the mitochondria through the carnitine palmitoyl transferase (CPT) system.

Comparative genomics of the vertebrate insulin/TOR signal transduction pathway: A network-level analysis of selective pressures

Complexity of biological function relies on large networks of interacting molecules. However, the evolutionary properties of these networks are not fully understood. It has been shown that selective pressures depend on the position of genes in the network. We have previously shown that in the Drosophila insulin/target of rapamycin (TOR) signal transduction pathway there is a correlation between the pathway position and the strength of purifying selection, with the downstream genes being most constrained. In this study, we investigated the evolutionary dynamics of this well-characterized pathway in vertebrates. More specifically, we determined the impact of natural selection on the evolution of 72 genes of this pathway. We found that in vertebrates there is a similar gradient of selective constraint in the insulin/TOR pathway to that found in Drosophila. This feature is neither the result of a polarity in the impact of positive selection nor of a series of factors affecting selective constraint levels (gene expression level and breadth, codon bias, protein length, and connectivity). We also found that pathway genes encoding physically interacting proteins tend to evolve under similar selective constraints. The results indicate that the architecture of the vertebrate insulin/TOR pathway constrains the molecular evolution of its components. Therefore, the polarity detected in Drosophila is neither specific nor incidental of this genus. Hence, although the underlying biological mechanisms remain unclear, these may be similar in both vertebrates and Drosophila.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Use of physiological constraints to identify quantitative design principles for gene expression in yeast adaptation to heat shock

Background: Understanding the relationship between gene expression changes, enzyme activity shifts, and the corresponding physiological adaptive response of organisms to environmental cues is crucial in explaining how cells cope with stress. For example, adaptation of yeast to heat shock involves a characteristic profile of changes to the expression levels of genes coding for enzymes of the glycolytic pathway and some of its branches. The experimental determination of changes in gene expression profiles provides a descriptive picture of the adaptive response to stress. However, it does not explain why a particular profile is selected for any given response. Results: We used mathematical models and analysis of in silico gene expression profiles (GEPs) to understand how changes in gene expression correlate to an efficient response of yeast cells to heat shock. An exhaustive set of GEPs, matched with the corresponding set of enzyme activities, was simulated and analyzed. The effectiveness of each profile in the response to heat shock was evaluated according to relevant physiological and functional criteria. The small subset of GEPs that lead to effective physiological responses after heat shock was identified as the result of the tuning of several evolutionary criteria. The experimentally observed transcriptional changes in response to heat shock belong to this set and can be explained by quantitative design principles at the physiological level that ultimately constrain changes in gene expression. Conclusion: Our theoretical approach suggests a method for understanding the combined effect of changes in the expression of multiple genes on the activity of metabolic pathways, and consequently on the adaptation of cellular metabolism to heat shock. This method identifies quantitative design principles that facilitate understating the response of the cell to stress.

The Saccharomyces cerevisiae Hot1p regulated gene YHR087W (HGI1) has a role in translation upon high glucose concentration stress.

Background While growing in natural environments yeasts can be affected by osmotic stress provoked by high glucose concentrations. The response to this adverse condition requires the HOG pathway and involves transcriptional and posttranscriptional mechanisms initiated by the phosphorylation of this protein, its translocation to the nucleus and activation of transcription factors. One of the genes induced to respond to this injury is YHR087W. It encodes for a protein structurally similar to the N-terminal region of human SBDS whose expression is also induced under other forms of stress and whose deletion determines growth defects at high glucose concentrations. Results In this work we show that YHR087W expression is regulated by several transcription factors depending on the particular stress condition, and Hot1p is particularly relevant for the induction at high glucose concentrations. In this situation, Hot1p, together to Sko1p, binds to YHR087W promoter in a Hog1p-dependent manner. Several evidences obtained indicate Yhr087wp"s role in translation. Firstly, and according to TAP purification experiments, it interacts with proteins involved in translation initiation. Besides, its deletion mutant shows growth defects in the presence of translation inhibitors and displays a slightly slower translation recovery after applying high glucose stress than the wild type strain. Analyses of the association of mRNAs to polysome fractions reveals a lower translation in the mutant strain of the mRNAs corresponding to genes GPD1, HSP78 and HSP104. Conclusions The data demonstrates that expression of Yhr087wp under high glucose concentration is controlled by Hot1p and Sko1p transcription factors, which bind to its promoter. Yhr087wp has a role in translation, maybe in the control of the synthesis of several stress response proteins, which could explain the lower levels of some of these proteins found in previous proteomic analyses and the growth defects of the deletion strain. Keywords: Saccharomyces cerevisiae; High glucose osmotic stress; Gene YHR087W; Gene expression; Translation; Hot1p; Hog1p; Polysomes

GD3 Synthase Overexpression Sensitizes Hepatocarcinoma Cells to Hypoxia and Reduces Tumor Growth by Suppressing the cSrc/NF-κB Survival Pathway

Abstract Background: Hypoxia-mediated HIF-1a stabilization and NF-kB activation play a key role in carcinogenesis by fostering cancer cell survival, angiogenesis and tumor invasion. Gangliosides are integral components of biological membranes with an increasingly recognized role as signaling intermediates. In particular, ganglioside GD3 has been characterized as a proapoptotic lipid effector by promoting cell death signaling and suppression of survival pathways. Thus, our aim was to analyze the role of GD3 in hypoxia susceptibility of hepatocarcinoma cells and in vivo tumor growth. Methodology/Principal Findings: We generated and characterized a human hepatocarcinoma cell line stably expressing GD3 synthase (Hep3B-GD3), which catalyzes the synthesis of GD3 from GM3. Despite increased GD3 levels (2-3 fold), no significant changes in cell morphology or growth were observed in Hep3B-GD3 cells compared to wild type Hep3B cells under normoxia. However, exposure of Hep3B-GD3 cells to hypoxia (2% O2) enhanced reactive oxygen species (ROS) generation, resulting in decreased cell survival, with similar findings observed in Hep3B cells exposed to increasing doses of exogenous GD3. In addition, hypoxia-induced c-Src phosphorylation at tyrosine residues, NF-kB activation and subsequent expression of Mn-SOD were observed in Hep3B cells but not in Hep3B-GD3 cells. Moreover, MnTBAP, an antioxidant with predominant SOD mimetic activity, reduced ROS generation, protecting Hep3B-GD3 cells from hypoxia-induced death. Finally, lower tumor growth, higher cell death and reduced Mn-SOD expression were observed in Hep3B-GD3 compared to Hep3B tumor xenografts. Conclusion: These findings underscore a role for GD3 in hypoxia susceptibility by disabling the c-Src/NF-kB survival pathway resulting in lower Mn-SOD expression, which may be of relevance in hepatocellular carcinoma therapy.

The SUD1 Gene Encodes a Putative E3 Ubiquitin Ligase and Is a Positive Regulator of 3-Hydroxy-3-Methylglutaryl Coenzyme A Reductase Activity in Arabidopsis

The 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) enzyme catalyzes the major rate-limiting step of the mevalonic acid (MVA) pathway from which sterols and other isoprenoids are synthesized. In contrast with our extensive knowledge of the regulation of HMGR in yeast and animals, little is known about this process in plants. To identify regulatory components of the MVA pathway in plants, we performed a genetic screen for second-site suppressor mutations of the Arabidopsis thaliana highly drought-sensitive drought hypersensitive2 (dry2) mutant that shows decreased squalene epoxidase activity. We show that mutations in SUPPRESSOR OF DRY2 DEFECTS1 (SUD1) gene recover most developmental defects in dry2 through changes in HMGR activity. SUD1 encodes a putative E3 ubiquitin ligase that shows sequence and structural similarity to yeast Degradation of a factor (Doa10) and human TEB4, components of the endoplasmic reticulum-associated degradation C (ERAD-C) pathway. While in yeast and animals, the alternative ERAD-L/ERAD-M pathway regulates HMGR activity by controlling protein stability, SUD1 regulates HMGR activity without apparent changes in protein content. These results highlight similarities, as well as important mechanistic differences, among the components involved in HMGR regulation in plants, yeast, and animals.

Three-step pathway engineering results in more incidence rate and higher emission of nerolidol and improved attraction of Diadegma semiclausum

The concentration and ratio of terpenoids in the headspace volatile blend of plants have a fundamental role in the communication of plants and insects. The sesquiterpene (E)-nerolidol is one of the important volatiles with effect on beneficial carnivores for biologic pest management in the field. To optimize de novo biosynthesis and reliable and uniform emission of (E)-nerolidol, we engineered different steps of the (E)-nerolidol biosynthesis pathway in Arabidopsis thaliana. Introduction of a mitochondrial nerolidol synthase gene mediates de novo emission of (E)-nerolidol and linalool. Co-expression of the mitochondrial FPS1 and cytosolic HMGR1 increased the number of emitting transgenic plants (incidence rate) and the emission rate of both volatiles. No association between the emission rate of transgenic volatiles and their growth inhibitory effect could be established. (E)-Nerolidol was to a large extent metabolized to non-volatile conjugates.

Gene Set Analysis for improving genetic association studies

Introduction. Genetic epidemiology is focused on the study of the genetic causes that determine health and diseases in populations. To achieve this goal a common strategy is to explore differences in genetic variability between diseased and nondiseased individuals. Usual markers of genetic variability are single nucleotide polymorphisms (SNPs) which are changes in just one base in the genome. The usual statistical approach in genetic epidemiology study is a marginal analysis, where each SNP is analyzed separately for association with the phenotype. Motivation. It has been observed, that for common diseases the single-SNP analysis is not very powerful for detecting genetic causing variants. In this work, we consider Gene Set Analysis (GSA) as an alternative to standard marginal association approaches. GSA aims to assess the overall association of a set of genetic variants with a phenotype and has the potential to detect subtle effects of variants in a gene or a pathway that might be missed when assessed individually. Objective. We present a new optimized implementation of a pair of gene set analysis methodologies for analyze the individual evidence of SNPs in biological pathways. We perform a simulation study for exploring the power of the proposed methodologies in a set of scenarios with different number of causal SNPs under different effect sizes. In addition, we compare the results with the usual single-SNP analysis method. Moreover, we show the advantage of using the proposed gene set approaches in the context of an Alzheimer disease case-control study where we explore the Reelin signal pathway.

The Plesiomonas shigelloides wbO1 gene cluster and the role of O1-antigen LPS in pathogenicity

The Plesiomonas shigelloides 302-73 strain (serotype O1) wb gene cluster encodes 15 proteins which are consistent with the chemical structure of the O1-antigen lypopolysaccharide (LPS) previously described for this strain. The P. shigelloides O1-antigen LPS export uses the Wzy-dependent pathway as correspond to heteropolysaccharides structures. By the isolation of two mutants lacking this O1-antigen LPS, we could establish that the presence of the O1-antigen LPS is crucial for to survive in serum mainly to become resistant to complement. Also, it is an important factor in the bacterial adhesion and invasion to some eukaryotic cells, and in the ability to form biofilms. This is the first report on the genetics from a P. shigelloides O-antigen LPS cluster (wb) not shared by Shigella like P. shigelloides O17, the only one reported until now.

Divergent evolution and purifying selection of the flaA gene sequences in Aeromonas

BACKGROUND: The bacterial flagellum is the most important organelle of motility in bacteria and plays a key role in many bacterial lifestyles, including virulence. The flagellum also provides a paradigm of how hierarchical gene regulation, intricate protein-protein interactions and controlled protein secretion can result in the assembly of a complex multi-protein structure tightly orchestrated in time and space. As if to stress its importance, plants and animals produce receptors specifically dedicated to the recognition of flagella. Aside from motility, the flagellum also moonlights as an adhesion and has been adapted by humans as a tool for peptide display. Flagellar sequence variation constitutes a marker with widespread potential uses for studies of population genetics and phylogeny of bacterial species. RESULTS: We sequenced the complete flagellin gene (flaA) in 18 different species and subspecies of Aeromonas. Sequences ranged in size from 870 (A. allosaccharophila) to 921 nucleotides (A. popoffii). The multiple alignment displayed 924 sites, 66 of which presented alignment gaps. The phylogenetic tree revealed the existence of two groups of species exhibiting different FlaA flagellins (FlaA1 and FlaA2). Maximum likelihood models of codon substitution were used to analyze flaA sequences. Likelihood ratio tests suggested a low variation in selective pressure among lineages, with an omega ratio of less than 1 indicating the presence of purifying selection in almost all cases. Only one site under potential diversifying selection was identified (isoleucine in position 179). However, 17 amino acid positions were inferred as sites that are likely to be under positive selection using the branch-site model. Ancestral reconstruction revealed that these 17 amino acids were among the amino acid changes detected in the ancestral sequence. CONCLUSION: The models applied to our set of sequences allowed us to determine the possible evolutionary pathway followed by the flaA gene in Aeromonas, suggesting that this gene have probably been evolving independently in the two groups of Aeromonas species since the divergence of a distant common ancestor after one or several episodes of positive selection. REVIEWERS: This article was reviewed by Alexey Kondrashov, John Logsdon and Olivier Tenaillon (nominated by Laurence D Hurst).