TRANSFAC: an integrated system for gene expression regulation

Multiple evolutionary origins of bat papillomaviruses

Infection by papillomaviruses (PVs) has been linked to different types of neoplasias, in both human and non-human hosts. Knowledge about PV diversity is essential to reliably infer the evolutionary history of these pathogens and to elucidate the link between infection and disease. We cloned and sequenced the complete genome of a novel PV, EhelPV1, isolated from hair bulbs from a captive straw-colored fruit bat Eidolon helvum (Pteropodidae, Chiroptera). We also retrieved partial sequences of the E1 and L1 genes from hair bulbs from a captive Indian flying fox Pteropus giganteus (Pteropodidae, Chiroptera). The detected virus (PgigPV1) presumably corresponded to a novel type as well. Maximum likelihood phylogenetic analyses were conducted using a representative collection of 132 PVs. EhelPV1 belonged to the Lambda+Mu-PV crown group and was most closely related to another bat PV, MschPV2. Both fragments of PgigPV1 were placed alongside with EhelPV1. The novel PVs were phylogenetically distant from other previously described bat PVs, namely MrPV1, MschPV1 and RaPV1. We have further characterized the sequence patterns of the E2-binding sites occurring in the upstream regulatory region of Lambda+Mu-PVs. Common fingerprints within this region are shared by certain PVs. However, there is not a sharp correspondence between the repertoire of transcription factor binding sites in the viral regulatory region and host range, tissue tropism or viral life style. Our results reinforce the hypothesis that PVs have undergone an initial radiation prior to the divergence of the mammalian hosts, giving rise to the present-day PV crown groups.

The Transcriptional Response of Neurotrophins and Their Tyrosine Kinase Receptors in Lumbar Sensorimotor Circuits to Spinal Cord Contusion is Affected by Injury Severity and Survival Time

Traumatic spinal cord injury (SCI) results in changes to the anatomical, neurochemical, and physiological properties of cells in the central and peripheral nervous system. Neurotrophins, acting by binding to their cognate Trk receptors on target cell membranes, contribute to modulation of anatomical, neurochemical, and physiological properties of neurons in sensorimotor circuits in both the intact and injured spinal cord. Neurotrophin signaling is associated with many post-SCI changes including maladaptive plasticity leading to pain and autonomic dysreflexia, but also therapeutic approaches such as training-induced locomotor improvement. Here we characterize expression of mRNA for neurotrophins and Trk receptors in lumbar dorsal root ganglia (DRG) and spinal cord after two different severities of mid-thoracic injury and at 6 and 12 weeks post-SCI. There was complex regulation that differed with tissue, injury severity, and survival time, including reversals of regulation between 6 and 12 weeks, and the data suggest that natural regulation of neurotrophins in the spinal cord may continue for months after birth. Our assessments determined that a coordination of gene expression emerged at the 12-week post-SCI time point and bioinformatic analyses address possible mechanisms. These data can inform studies meant to determine the role of the neurotrophin signaling system in post-SCI function and plasticity, and studies using this signaling system as a therapeutic approach.

Forkhead box M1 is regulated by heat shock factor 1 and promotes glioma cells survival under heat shock stress

The forkhead box M1 (FoxM1) is a key transcription factor regulating multiple aspects of cell biology. Prior studies have shown that FoxM1 is overexpressed in a variety of human tumors, including brain tumor, and plays a critical role in cancer development and progression. In this study we found that FoxM1 was up-regulated by heat shock factor 1 (HSF1) under heat shock stress condition in multiple cell lines. Knockdown of HSF1 with HSF1 siRNA or inhibition of HSF1 with a HSF1 inhibitor abrogated heat shock-induced expression of FoxM1. Genetic deletion of HSF1 in mouse embryo fibroblast cells also abolished heat shock stress-induced FoxM1 expression. Moreover, we showed that HSF1 directly bound to FoxM1 promoter and increased FoxM1 promoter activity. Furthermore, we demonstrated that FoxM1 was required for the G(2)-M phase progression through regulating Cdc2, Cdc20, and Cdc25B under a mild heat shock stress but enhanced cell survival under lethal heat shock stress condition. Finally, in human glioblastoma specimens, FoxM1 overexpression correlated with elevated HSF1 expression. Our results indicate that FoxM1 is regulated by HSF1 and is critical for HSF1-mediated heat shock response. We demonstrated a novel mechanism of stress resistance controlled by HSF1 and a new HSF-FoxM1 connection that mediates cellular thermotolerance.

[The databases of transcription factors.]

The control of gene transcription is a critical level of gene expression regulation. The interactions between transcription factors (TF) and their DNA binding sites (TFBS) play a key role at this level. In order to decipher the molecular mechanism of the interactions of TFs with TFBSs and construct transcription regulatory network, it is necessary to systematically collect, save, and analyze the information of discovered TFs and their TFBSs. In recent years, multiple TF and TFBS-related databases have been established. These databeses significantly promoted the TF-related studies in the fields of molecular biology, bioinformatics, and system biology. This paper summarized the contents, characteristics, access, and advances of main TFs and TFBSs-related databases, including TRANSFAC, JASPAR, TFDB, TRRD, TRED, PAZAR, MAPPER and others.

An integrative computational systems biology approach identifies differentially regulated dynamic transcriptome signatures which drive the initiation of human T helper cell differentiation

Background

A proper balance between different T helper (Th) cell subsets is necessary for normal functioning of the adaptive immune system. Revealing key genes and pathways driving the differentiation to distinct Th cell lineages provides important insight into underlying molecular mechanisms and new opportunities for modulating the immune response. Previous computational methods to quantify and visualize kinetic differential expression data of three or more lineages to identify reciprocally regulated genes have relied on clustering approaches and regression methods which have time as a factor, but have lacked methods which explicitly model temporal behavior.

Results

We studied transcriptional dynamics of human umbilical cord blood T helper cells cultured in absence and presence of cytokines promoting Th1 or Th2 differentiation. To identify genes that exhibit distinct lineage commitment dynamics and are specific for initiating differentiation to different Th cell subsets, we developed a novel computational methodology (LIGAP) allowing integrative analysis and visualization of multiple lineages over whole time-course profiles. Applying LIGAP to time-course data from multiple Th cell lineages, we identified and experimentally validated several differentially regulated Th cell subset specific genes as well as reciprocally regulated genes. Combining differentially regulated transcriptional profiles with transcription factor binding site and pathway information, we identified previously known and new putative transcriptional mechanisms involved in Th cell subset differentiation. All differentially regulated genes among the lineages together with an implementation of LIGAP are provided as an open-source resource.

Conclusions

The LIGAP method is widely applicable to quantify differential time-course dynamics of many types of datasets and generalizes to any number of conditions. It summarizes all the time-course measurements together with the associated uncertainty for visualization and manual assessment purposes. Here we identified novel human Th subset specific transcripts as well as regulatory mechanisms important for the initiation of the Th cell subset differentiation.

Investigation of CD28 gene polymorphisms in patients with sporadic breast cancer in a Chinese Han population in Northeast China

Background

CD28 is one of a number of costimulatory molecules that play crucial roles in immune regulation and homeostasis. Accumulating evidence indicates that immune factors influence breast carcinogenesis. To clarify the relationships between polymorphisms in the CD28 gene and breast carcinogenesis, a case-control study was conducted in women from Heilongjiang Province in northeast of China.

Methodology/Principal Findings

Our research subjects consisted of 565 female patients with sporadic breast cancer and 605 age- and sex-matched healthy controls. In total, 12 single nucleotide polymorphisms (SNPs) in the CD28 gene were successfully determined using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. The relationship between the CD28 variants and clinical features, including histological grade, tumor size, lymph node metastasis, human epidermal growth factor receptor 2 (C-erbB2), estrogen receptor (ER), progesterone receptor (PR), and tumor protein 53 (P53) status were analyzed. A statistically significant association was observed between rs3116496 and breast cancer risk under different genetic models (additive P = 0.0164, dominant P = 0.0042). Different distributions of the rs3116496 ‘T’ allele were found in patients and controls, which remained significant after correcting the P value for multiple testing using Haploview with 10,000 permutations (corrected P = 0.0384). In addition, significant associations were observed between rs3116487/rs3116494 (D’ = 1, r² = 0.99) and clinicopathological features such as C-erbB2 and ER status, in breast cancer patients.

Conclusions/Significance

Our findings indicate that CD28 gene polymorphisms contribute to sporadic breast cancer risk and have a significant association with clinicopathological features in a northeast Chinese Han population.

MuMoD: a Bayesian approach to detect multiple modes of protein-DNA binding from genome-wide ChIP data

High-throughput chromatin immunoprecipitation has become the method of choice for identifying genomic regions bound by a protein. Such regions are then investigated for overrepresented sequence motifs, the assumption being that they must correspond to the binding specificity of the profiled protein. However this approach often fails: many bound regions do not contain the 'expected' motif. This is because binding DNA directly at its recognition site is not the only way the protein can cause the region to immunoprecipitate. Its binding specificity can change through association with different co-factors, it can bind DNA indirectly, through intermediaries, or even enforce its function through long-range chromosomal interactions. Conventional motif discovery methods, though largely capable of identifying overrepresented motifs from bound regions, lack the ability to characterize such diverse modes of protein-DNA binding and binding specificities. We present a novel Bayesian method that identifies distinct protein-DNA binding mechanisms without relying on any motif database. The method successfully identifies co-factors of proteins that do not bind DNA directly, such as mediator and p300. It also predicts literature-supported enhancer-promoter interactions. Even for well-studied direct-binding proteins, this method provides compelling evidence for previously uncharacterized dependencies within positions of binding sites, long-range chromosomal interactions and dimerization.

Functional classification of skeletal muscle networks. I. Normal physiology

Extensive measurements of the parts list of human skeletal muscle through transcriptomics and other phenotypic assays offer the opportunity to reconstruct detailed functional models. Through integration of vast amounts of data present in databases and extant knowledge of muscle function combined with robust analyses that include a clustering approach, we present both a protein parts list and network models for skeletal muscle function. The model comprises the four key functional family networks that coexist within a functional space; namely, excitation-activation family (forward pathways that transmit a motoneuronal command signal into the spatial volume of the cell and then use Ca(2+) fluxes to bind Ca(2+) to troponin C sites on F-actin filaments, plus transmembrane pumps that maintain transmission capacity); mechanical transmission family (a sophisticated three-dimensional mechanical apparatus that bidirectionally couples the millions of actin-myosin nanomotors with external axial tensile forces at insertion sites); metabolic and bioenergetics family (pathways that supply energy for the skeletal muscle function under widely varying demands and provide for other cellular processes); and signaling-production family (which represents various sensing, signal transduction, and nuclear infrastructure that controls the turn over and structural integrity and regulates the maintenance, regeneration, and remodeling of the muscle). Within each family, we identify subfamilies that function as a unit through analysis of large-scale transcription profiles of muscle and other tissues. This comprehensive network model provides a framework for exploring functional mechanisms of the skeletal muscle in normal and pathophysiology, as well as for quantitative modeling.

Conserved noncoding sequences highlight shared components of regulatory networks in dicotyledonous plants

Conserved noncoding sequences (CNSs) in DNA are reliable pointers to regulatory elements controlling gene expression. Using a comparative genomics approach with four dicotyledonous plant species (Arabidopsis thaliana, papaya [Carica papaya], poplar [Populus trichocarpa], and grape [Vitis vinifera]), we detected hundreds of CNSs upstream of Arabidopsis genes. Distinct positioning, length, and enrichment for transcription factor binding sites suggest these CNSs play a functional role in transcriptional regulation. The enrichment of transcription factors within the set of genes associated with CNS is consistent with the hypothesis that together they form part of a conserved transcriptional network whose function is to regulate other transcription factors and control development. We identified a set of promoters where regulatory mechanisms are likely to be shared between the model organism Arabidopsis and other dicots, providing areas of focus for further research.

Uncovering transcription factor modules using one- and three-dimensional analyses

Transcriptional regulation is a critical mediator of many normal cellular processes, as well as disease progression. Transcription factors (TFs) often co-localize at cis-regulatory elements on the DNA, form protein complexes, and collaboratively regulate gene expression. Machine learning and Bayesian approaches have been used to identify TF modules in a one-dimensional context. However, recent studies using high throughput technologies have shown that TF interactions should also be considered in three-dimensional nuclear space. Here, we describe methods for identifying TF modules and discuss how moving from a one-dimensional to a three-dimensional paradigm, along with integrated experimental and computational approaches, can lead to a better understanding of TF association networks.

Itga2b regulation at the onset of definitive hematopoiesis and commitment to differentiation

Product of the Itga2b gene, CD41 contributes to hematopoietic stem cell (HSC) and megakaryocyte/platelet functions. CD41 expression marks the onset of definitive hematopoiesis in the embryo where it participates in regulating the numbers of multipotential progenitors. Key to platelet aggregation, CD41 expression also characterises their precursor, the megakaryocyte, and is specifically up regulated during megakaryopoiesis. Though phenotypically unique, megakaryocytes and HSC share numerous features, including key transcription factors, which could indicate common sub-regulatory networks. In these respects, Itga2b can serve as a paradigm to study features of both developmental-stage and HSC- versus megakaryocyte-specific regulations. By comparing different cellular contexts, we highlight a mechanism by which internal promoters participate in Itga2b regulation. A developmental process connects epigenetic regulation and promoter switching leading to CD41 expression in HSC. Interestingly, a similar process can be observed at the Mpl locus, which codes for another receptor that defines both HSC and megakaryocyte identities. Our study shows that Itga2b expression is controlled by lineage-specific networks and associates with H4K8ac in megakaryocyte or H3K27me3 in the multipotential hematopoietic cell line HPC7. Correlating with the decrease in H3K27me3 at the Itga2b Iocus, we find that following commitment to megakaryocyte differentiation, the H3K27 demethylase Jmjd3 up-regulation influences both Itga2b and Mpl expression.

A metabolic prosurvival role for PML in breast cancer

Cancer cells exhibit an aberrant metabolism that facilitates more efficient production of biomass and hence tumor growth and progression. However, the genetic cues modulating this metabolic switch remain largely undetermined. We identified a metabolic function for the promyelocytic leukemia (PML) gene, uncovering an unexpected role for this bona fide tumor suppressor in breast cancer cell survival. We found that PML acted as both a negative regulator of PPARγ coactivator 1A (PGC1A) acetylation and a potent activator of PPAR signaling and fatty acid oxidation. We further showed that PML promoted ATP production and inhibited anoikis. Importantly, PML expression allowed luminal filling in 3D basement membrane breast culture models, an effect that was reverted by the pharmacological inhibition of fatty acid oxidation. Additionally, immunohistochemical analysis of breast cancer biopsies revealed that PML was overexpressed in a subset of breast cancers and enriched in triple-negative cases. Indeed, PML expression in breast cancer correlated strikingly with reduced time to recurrence, a gene signature of poor prognosis, and activated PPAR signaling. These findings have important therapeutic implications, as PML and its key role in fatty acid oxidation metabolism are amenable to pharmacological suppression, a potential future mode of cancer prevention and treatment.

Expression of complement component C7 and involvement in innate immune responses to bacteria in grass carp

Activation of the complement system, which forms a major part of the innate immune system, results in the formation of the terminal complement complex. The complement component, C7, plays an integral role in the assembly of this complex within target cell membranes. In this study, C7 was isolated and characterized from grass carp, an important cultured fish in China. The predicted amino acid sequence of C7 cDNA (2644 bp) exhibited 55.4 and 48.3% homology with trout C7-1 and zebrafish C7, respectively. The grass carp C7 gene was consisted of 18 exons and 17 introns. C7 gene expression was detected in the trunk kidney, liver, head kidney, skin, spleen, heart and intestine. Significant changes in C7 transcript expression (>20-fold) were detected following Aeromonas hydrophila infection, indicating C7 involvement in innate immune responses to bacteria in teleost fish.

Oct-1 cooperates with the TATA binding initiation complex to control rapid transcription of human iNOS

Expression of the human inducible nitric oxide synthase (hiNOS) is generally undetectable in resting cells, but stimulation by a variety of signals including cytokines induces transcription in most cell types. The tight transcriptional regulation of the enzyme is a complex mechanism many aspects of which remain unknown. Here, we describe an octamer (Oct) element in hiNOS proximal promoter, located close to the TATA box. This site constitutively binds Oct-1 and its deletion abrogates cytokine-induced transcription, showing that it is indispensable though not sufficient for transcription. Increasing the distance between Oct and the TATA box by inserting inert DNA sequence inhibits transcription, and footprinting of this region shows no other protein binding in resting cells, suggesting an interaction between the two complexes. Chromatin immunoprecipitation assays detect the presence of Oct-1, RNA polymerase II and trimethyl K4 histone H3 on the proximal promoter in resting cells, confirming that the gene is primed for transcription before stimulation. RT-PCR of various fragments along the hiNOS gene shows that transcription is initiated in resting cells and this is inhibited by interference with Oct-1 binding to the proximal site of the promoter. We propose that, through interaction with the initiation complex, Oct-1 regulates hiNOS transcription by priming the gene for the rapid response required in an immune response.

Network-enabled gene expression analysis

BACKGROUND

Although genome-scale expression experiments are performed routinely in biomedical research, methods of analysis remain simplistic and their interpretation challenging. The conventional approach is to compare the expression of each gene, one at a time, between treatment groups. This implicitly treats the gene expression levels as independent, but they are in fact highly interdependent, and exploiting this enables substantial power gains to be realized.

RESULTS

We assume that information on the dependence structure between the expression levels of a set of genes is available in the form of a Bayesian network (directed acyclic graph), derived from external resources. We show how to analyze gene expression data conditional on this network. Genes whose expression is directly affected by treatment may be identified using tests for the independence of each gene and treatment, conditional on the parents of the gene in the network. We apply this approach to two datasets: one from a hepatotoxicity study in rats using a PPAR pathway, and the other from a study of the effects of smoking on the epithelial transcriptome, using a global transcription factor network.

CONCLUSIONS

The proposed method is straightforward, simple to implement, gives rise to substantial power gains, and may assist in relating the experimental results to the underlying biology.

Effect of positional dependence and alignment strategy on modeling transcription factor binding sites

Background

Many consensus-based and Position Weight Matrix-based methods for recognizing transcription factor binding sites (TFBS) are not well suited to the variability in the lengths of binding sites. Besides, many methods discard known binding sites while building the model. Moreover, the impact of Information Content (IC) and the positional dependence of nucleotides within an aligned set of TFBS has not been well researched for modeling variable-length binding sites. In this paper, we propose ML-Consensus (Mixed-Length Consensus): a consensus model for variable-length TFBS which does not exclude any reported binding sites.

Methods

We consider Pairwise Score (PS) as a measure of positional dependence of nucleotides within an alignment of TFBS. We investigate how the prediction accuracy of ML-Consensus is affected by the incorporation of IC and PS with a particular binding site alignment strategy. We perform cross-validations for datasets of six species from the TRANSFAC public database, and analyze the results using ROC curves and the Wilcoxon matched-pair signed-ranks test.

Results

We observe that the incorporation of IC and PS in ML-Consensus results in statistically significant improvement in the prediction accuracy of the model. Moreover, the existence of a core region among the known binding sites (of any length) is witnessed by the pairwise coexistence of nucleotides within the core length.

Conclusions

These observations suggest the possibility of an efficient multiple sequence alignment algorithm for aligning TFBS, accommodating known binding sites of any length, for optimal (or near-optimal) TFBS prediction. However, designing such an algorithm is a matter of further investigation.

A defect in dystrophin causes a novel porcine stress syndrome

Background

Losses of slaughter-weight pigs due to transport stress are both welfare and economic concerns to pork producers. Historically, the HAL-1843 mutation in ryanodine receptor 1 was considered responsible for most of the losses; however, DNA testing has effectively eliminated this mutation from commercial herds. We identified two sibling barrows in the USMARC swine herd that died from apparent symptoms of a stress syndrome after transport at 12 weeks of age. The symptoms included open-mouth breathing, skin discoloration, vocalization and loss of mobility.

Results

We repeated the original mating along with sire-daughter matings to produce additional offspring. At 8 weeks of age, heart rate and electrocardiographs (ECG) were monitored during isoflurane anesthesia challenge (3% for 3 min). Four males from the original sire-dam mating and two males from a sire-daughter mating died after one minute of anesthesia. Animals from additional litters were identified as having a stress response, sometimes resulting in death, during regular processing and weighing. Affected animals had elevated plasma creatine phosphokinase (CPK) levels before and immediately after isoflurane challenge and cardiac arrhythmias. A pedigree containing 250 pigs, including 49 affected animals, was genotyped with the Illumina PorcineSNP60 Beadchip and only one chromosomal region, SSCX at 25.1-27.7 Mb over the dystrophin gene (DMD), was significantly associated with the syndrome. An arginine to tryptophan (R1958W) polymorphism in exon 41 of DMD was the most significant marker associated with stress susceptibility. Immunoblots of affected heart and skeletal muscle showed a dramatic reduction of dystrophin protein and histopathology of affected hearts indicated muscle fiber degeneration.

Conclusions

A novel stress syndrome was characterized in pigs and the causative genetic factor most likely resides within DMD that results in less dystrophin protein and cardiac abnormalities that can lead to death under stressful conditions. The identification of predictive markers will allow us to determine the prevalence of this disease in commercial swine populations. This defect also provides a unique biomedical model for human cardiomyopathy associated with muscular dystrophy that may be superior to those available because of the similarities in anatomy and physiology and allow advances in gene therapies for human disease.

The involvement of heat shock protein and cytochrome P450 genes in response to UV-A exposure in the beetle Tribolium castaneum

Sunlight is an important environmental factor that affects all living organisms on Earth. Ultraviolet A (UV-A) is one of the many frequency bands found in sunlight. Many animals use UV-A to attain visual cues, for example, in foraging and mate selection. However, UV-A can also induce damage, such as oxidative stress, DNA lesions and apoptosis. In the present study, we investigated the effects of UV-A on the survival, fecundity and expression profiles of several stress-responsive genes belonging to the heat shock protein (Hsp) and the cytochrome CYP6BQ families from the adult red flour beetle, Tribolium castaneum (Herbst) (Coleoptera: Tenebrionidae). The results showed that short-term UV-A exposure (365 nm, <4h) did not influence the survival or fecundity of the beetles; however, Hsp27, Hsp68, Hsp83, CYP6BQ4 and CYP6BQ8 mRNA levels significantly increased during the first 2h of UV-A exposure. Among them, Hsp68 was the most highly up-regulated, increasing by 8.9-fold. These results indicate that these genes may participate in the defense against harmful UV-A radiation. In addition, we investigated the potential transcription factor binding motifs (TFBMs) in the promoter sequences of genes induced in similar pattern from the Hsp and P450 gene families; the results indicated that, these motifs are highly homologous to environmental stress transcription factor binding sites in mammals. Our experiments revealed that UV-A irradiation could influence the expression profile of stress-responsive genes, such as Hsps and P450s, which have universal TFBMs, and that these genes may be involved in reducing the ecological challenges posed by irradiation.

GPUmotif: an ultra-fast and energy-efficient motif analysis program using graphics processing units

Computational detection of TF binding patterns has become an indispensable tool in functional genomics research. With the rapid advance of new sequencing technologies, large amounts of protein-DNA interaction data have been produced. Analyzing this data can provide substantial insight into the mechanisms of transcriptional regulation. However, the massive amount of sequence data presents daunting challenges. In our previous work, we have developed a novel algorithm called Hybrid Motif Sampler (HMS) that enables more scalable and accurate motif analysis. Despite much improvement, HMS is still time-consuming due to the requirement to calculate matching probabilities position-by-position. Using the NVIDIA CUDA toolkit, we developed a graphics processing unit (GPU)-accelerated motif analysis program named GPUmotif. We proposed a "fragmentation" technique to hide data transfer time between memories. Performance comparison studies showed that commonly-used model-based motif scan and de novo motif finding procedures such as HMS can be dramatically accelerated when running GPUmotif on NVIDIA graphics cards. As a result, energy consumption can also be greatly reduced when running motif analysis using GPUmotif. The GPUmotif program is freely available at http://sourceforge.net/projects/gpumotif/

Integration of MicroRNA databases to study MicroRNAs associated with multiple sclerosis

MicroRNAs (miRNAs) are small non-coding RNAs which regulate many genes post-transcriptionally. In various contexts of medical science, miRNAs gained increasing attention over the last few years. Analyzing the functions, interactions and cellular effects of miRNAs is a very complex and challenging task. Many miRNA databases with diverse data contents have been developed. Here, we demonstrate how to integrate their information in a reasonable way on a set of miRNAs that were found to be dysregulated in the blood of patients with multiple sclerosis (MS). Using the miR2Disease database, we retrieved 16 miRNAs associated with MS according to four different studies. We studied the predicted and experimentally validated target genes of these miRNAs, their expression profiles in different blood cell types and brain tissues, the pathways and biological processes affected by these miRNAs as well as their regulation by transcription factors. Only miRNA-mRNA interactions that were predicted by at least seven different prediction algorithms were considered. This resulted in a network of 1,498 target genes. In this network, the MS-associated miRNAs hsa-miR-20a-5p and hsa-miR-20b-5p occurred as central hubs regulating about 500 genes each. Strikingly, many of the putative target genes play a role in T cell activation and signaling, and many have transcription factor activity. The latter suggests that miRNAs often act as regulators of regulators with many secondary effects on gene expression. Our present work provides a guideline on how information of different databases can be integrated in the analysis of miRNAs. Future investigations of miRNAs shall help to better understand the mechanisms underlying different diseases and their treatments.

Role of pharmacogenomics in the management of traditional and novel oral anticoagulants

Warfarin is the most commonly prescribed oral anticoagulant. However, it remains a difficult drug to manage mostly because of its narrow therapeutic index and wide interpatient variability in anticoagulant effects. Over the past decade, there has been substantial progress in our understanding of genetic contributions to variable warfarin response, particularly with regard to warfarin dose requirements. The genes encoding for cytochrome P450 (CYP) 2C9 (CYP2C9) and vitamin K epoxide reductase complex subunit 1 (VKORC1) are the major genetic determinants of warfarin pharmacokinetics and pharmacodynamics, respectively. Numerous studies have demonstrated significant contributions of these genes to warfarin dose requirements. The CYP2C9 gene has also been associated with bleeding risk with warfarin. The CYP4F2 gene influences vitamin K availability and makes minor contributions to warfarin dose requirements. Less is known about genes influencing warfarin response in African-American patients compared with other racial groups, but this is the focus of ongoing research. Several warfarin pharmacogenetic dosing algorithms and United States Food and Drug Administration-cleared genotyping tests are available for clinical use. Clinical trials are ongoing to determine the clinical utility and cost-effectiveness of genotypeguided warfarin dosing. Results from these trials will likely influence clinical uptake and third party payer reimbursement for genotype-guided warfarin therapy. There is still a lack of pharmacogenetic data for the newly approved oral anticoagulants, dabigatran and rivaroxaban, and with other oral anticoagulants in the research and development pipeline. These data, once known, could be of great importance as routine monitoring parameters for these agents are not available.

Population differences in transcript-regulator expression quantitative trait loci

Gene expression quantitative trait loci (eQTL) are useful for identifying single nucleotide polymorphisms (SNPs) associated with diseases. At times, a genetic variant may be associated with a master regulator involved in the manifestation of a disease. The downstream target genes of the master regulator are typically co-expressed and share biological function. Therefore, it is practical to screen for eQTLs by identifying SNPs associated with the targets of a transcript-regulator (TR). We used a multivariate regression with the gene expression of known targets of TRs and SNPs to identify TReQTLs in European (CEU) and African (YRI) HapMap populations. A nominal p-value of <1×10(-6) revealed 234 SNPs in CEU and 154 in YRI as TReQTLs. These represent 36 independent (tag) SNPs in CEU and 39 in YRI affecting the downstream targets of 25 and 36 TRs respectively. At a false discovery rate (FDR) = 45%, one cis-acting tag SNP (within 1 kb of a gene) in each population was identified as a TReQTL. In CEU, the SNP (rs16858621) in Pcnxl2 was found to be associated with the genes regulated by CREM whereas in YRI, the SNP (rs16909324) was linked to the targets of miRNA hsa-miR-125a. To infer the pathways that regulate expression, we ranked TReQTLs by connectivity within the structure of biological process subtrees. One TReQTL SNP (rs3790904) in CEU maps to Lphn2 and is associated (nominal p-value = 8.1×10(-7)) with the targets of the X-linked breast cancer suppressor Foxp3. The structure of the biological process subtree and a gene interaction network of the TReQTL revealed that tumor necrosis factor, NF-kappaB and variants in G-protein coupled receptors signaling may play a central role as communicators in Foxp3 functional regulation. The potential pleiotropic effect of the Foxp3 TReQTLs was gleaned from integrating mRNA-Seq data and SNP-set enrichment into the analysis.

Unsupervised pattern discovery in human chromatin structure through genomic segmentation

We trained Segway, a dynamic Bayesian network method, simultaneously on chromatin data from multiple experiments, including positions of histone modifications, transcription-factor binding and open chromatin, all derived from a human chronic myeloid leukemia cell line. In an unsupervised fashion, we identified patterns associated with transcription start sites, gene ends, enhancers, transcriptional regulator CTCF-binding regions and repressed regions. Software and genome browser tracks are at http://noble.gs.washington.edu/proj/segway/.

Genome-wide association scan identifies a risk locus for preeclampsia on 2q14, near the inhibin, beta B gene

Elucidating the genetic architecture of preeclampsia is a major goal in obstetric medicine. We have performed a genome-wide association study (GWAS) for preeclampsia in unrelated Australian individuals of Caucasian ancestry using the Illumina OmniExpress-12 BeadChip to successfully genotype 648,175 SNPs in 538 preeclampsia cases and 540 normal pregnancy controls. Two SNP associations (rs7579169, p = 3.58×10⁻⁷, OR = 1.57; rs12711941, p = 4.26×10⁻⁷, OR = 1.56) satisfied our genome-wide significance threshold (modified Bonferroni p<5.11×10⁻⁷). These SNPs reside in an intergenic region less than 15 kb downstream from the 3′ terminus of the Inhibin, beta B (INHBB) gene on 2q14.2. They are in linkage disequilibrium (LD) with each other (r² = 0.92), but not (r²<0.80) with any other genotyped SNP ±250 kb. DNA re-sequencing in and around the INHBB structural gene identified an additional 25 variants. Of the 21 variants that we successfully genotyped back in the case-control cohort the most significant association observed was for a third intergenic SNP (rs7576192, p = 1.48×10⁻⁷, OR = 1.59) in strong LD with the two significant GWAS SNPs (r²>0.92). We attempted to provide evidence of a putative regulatory role for these SNPs using bioinformatic analyses and found that they all reside within regions of low sequence conservation and/or low complexity, suggesting functional importance is low. We also explored the mRNA expression in decidua of genes ±500 kb of INHBB and found a nominally significant correlation between a transcript encoded by the EPB41L5 gene, ∼250 kb centromeric to INHBB, and preeclampsia (p = 0.03). We were unable to replicate the associations shown by the significant GWAS SNPs in case-control cohorts from Norway and Finland, leading us to conclude that it is more likely that these SNPs are in LD with as yet unidentified causal variant(s).

pcaGoPromoter--an R package for biological and regulatory interpretation of principal components in genome-wide gene expression data

Analyzing data obtained from genome-wide gene expression experiments is challenging due to the quantity of variables, the need for multivariate analyses, and the demands of managing large amounts of data. Here we present the R package pcaGoPromoter, which facilitates the interpretation of genome-wide expression data and overcomes the aforementioned problems. In the first step, principal component analysis (PCA) is applied to survey any differences between experiments and possible groupings. The next step is the interpretation of the principal components with respect to both biological function and regulation by predicted transcription factor binding sites. The robustness of the results is evaluated using cross-validation, and illustrative plots of PCA scores and gene ontology terms are available. pcaGoPromoter works with any platform that uses gene symbols or Entrez IDs as probe identifiers. In addition, support for several popular Affymetrix GeneChip platforms is provided. To illustrate the features of the pcaGoPromoter package a serum stimulation experiment was performed and the genome-wide gene expression in the resulting samples was profiled using the Affymetrix Human Genome U133 Plus 2.0 chip. Array data were analyzed using pcaGoPromoter package tools, resulting in a clear separation of the experiments into three groups: controls, serum only and serum with inhibitor. Functional annotation of the axes in the PCA score plot showed the expected serum-promoted biological processes, e.g., cell cycle progression and the predicted involvement of expected transcription factors, including E2F. In addition, unexpected results, e.g., cholesterol synthesis in serum-depleted cells and NF-κB activation in inhibitor treated cells, were noted. In summary, the pcaGoPromoter R package provides a collection of tools for analyzing gene expression data. These tools give an overview of the input data via PCA, functional interpretation by gene ontology terms (biological processes), and an indication of the involvement of possible transcription factors.

Tcra enhancer activation by inducible transcription factors downstream of pre-TCR signaling

The Tcra enhancer (Eα) is essential for pre-TCR-mediated activation of germline transcription and V(D)J recombination. Eα is considered an archetypical enhanceosome that acts through the functional synergy and cooperative binding of multiple transcription factors. Based on dimethylsulfate genomic footprinting experiments, there has been a long-standing paradox regarding Eα activation in the absence of differences in enhancer occupancy. Our data provide the molecular mechanism of Eα activation and an explanation of this paradox. We found that germline transcriptional activation of Tcra is dependent on constant phospholipase Cγ, as well as calcineurin- and MAPK/ERK-mediated signaling, indicating that inducible transcription factors are crucially involved. NFAT, AP-1, and early growth response factor 1, together with CREB-binding protein/p300 coactivators, bind to Eα as part of an active enhanceosome assembled during pre-TCR signaling. We favor a scenario in which the binding of lymphoid-restricted and constitutive transcription factors to Eα prior to its activation forms a regulatory scaffold to recruit factors induced by pre-TCR signaling. Thus, the combinatorial assembly of tissue- and signal-specific transcription factors dictates the Eα function. This mechanism for enhancer activation may represent a general paradigm in tissue-restricted and stimulus-responsive gene regulation.

Meta-analysis identifies common variants associated with body mass index in east Asians

Multiple genetic loci associated with obesity or body mass index (BMI) have been identified through genome-wide association studies conducted predominantly in populations of European ancestry. We conducted a meta-analysis of associations between BMI and approximately 2.4 million SNPs in 27,715 East Asians, followed by in silico and de novo replication in 37,691 and 17,642 additional East Asians, respectively. We identified ten BMI-associated loci at the genome-wide significance level (P<5.0×10⁻⁸), including seven previously identified loci (FTO, SEC16B, MC4R, GIPR/QPCTL, ADCY3/RBJ, BDNF, and MAP2K5) and three novel loci in or near the CDKAL1,PCSK1, and GP2 genes. Three additional loci nearly reached the genome-wide significance threshold, including two previously identified loci in the GNPDA2 and TFAP2B genes and a new locus near PAX6, which all had P<5.0×10⁻⁷. Findings from this study may shed light on new pathways involved in obesity and demonstrate the value of conducting genetic studies in non-European populations.

Promoter variants determine γ-aminobutyric acid homeostasis-related gene transcription in human epileptic hippocampi

The functional consequences of single nucleotide polymorphisms associated with episodic brain disorders such as epilepsy and depression are unclear. Allelic associations with generalized epilepsies have been reported for single nucleotide polymorphisms rs1883415 (ALDH5A1; succinic semialdehyde dehydrogenase) and rs4906902 (GABRB3; GABAA β3), both of which are present in the 5' regulatory region of genes involved in γ-aminobutyric acid (GABA) homeostasis. To address their allelic association with episodic brain disorders and allele-specific impact on the transcriptional regulation of these genes in human brain tissue, DNA and messenger RNA (mRNA) isolated from hippocampi were obtained at epilepsy surgery of 146 pharmacoresistant mesial temporal lobe epilepsy (mTLE) patients and from 651 healthy controls. We found that the C allele of rs1883415 is accumulated to a greater extentin mTLE versus controls. By real-time quantitative reverse transcription-polymerase chain reaction analyses, individuals homozygous for the C allele showed higher ALDH5A1 mRNA expression. The rs4906902 G allele of the GABRB3 gene was overrepresented in mTLE patients with depression; individuals homozygous for the G allele showed reduced GABRB3 mRNA expression. Bioinformatic analyses suggest that rs1883415 and rs4906902 alter the DNA binding affinity of the transcription factors Egr-3 in ALDH5A1 and MEF-2 in GABRB3 promoters, respectively. Using in vitro luciferase transfection assays, we observed that, in both cases, the transcription factors regulate gene expression depending on the allelic variant in the same direction as in the human hippocampi. Our data suggest that distinct promoter variants may sensitize individuals for differential, potentially stimulus-induced alterations of GABA homeostasis-relevant gene expression. This might contribute to the episodic onset of symptoms and point to new targets for pharmacotherapies.

GPMiner: an integrated system for mining combinatorial cis-regulatory elements in mammalian gene group

Background

Sequence features in promoter regions are involved in regulating gene transcription initiation. Although numerous computational methods have been developed for predicting transcriptional start sites (TSSs) or transcription factor (TF) binding sites (TFBSs), they lack annotations for do not consider some important regulatory features such as CpG islands, tandem repeats, the TATA box, CCAAT box, GC box, over-represented oligonucleotides, DNA stability, and GC content. Additionally, the combinatorial interaction of TFs regulates the gene group that is associated with same expression pattern. To investigate gene transcriptional regulation, an integrated system that annotates regulatory features in a promoter sequence and detects co-regulation of TFs in a group of genes is needed.

Results

This work identifies TSSs and regulatory features in a promoter sequence, and recognizes co-occurrence of cis-regulatory elements in co-expressed genes using a novel system. Three well-known TSS prediction tools are incorporated with orthologous conserved features, such as CpG islands, nucleotide composition, over-represented hexamer nucleotides, and DNA stability, to construct the novel Gene Promoter Miner (GPMiner) using a support vector machine (SVM). According to five-fold cross-validation results, the predictive sensitivity and specificity are both roughly 80%. The proposed system allows users to input a group of gene names/symbols, enabling the co-occurrence of TFBSs to be determined. Additionally, an input sequence can also be analyzed for homogeneity of experimental mammalian promoter sequences, and conserved regulatory features between homologous promoters can be observed through cross-species analysis. After identifying promoter regions, regulatory features are visualized graphically to facilitate gene promoter observations.

Conclusions

The GPMiner, which has a user-friendly input/output interface, has numerous benefits in analyzing human and mouse promoters. The proposed system is freely available at http://GPMiner.mbc.nctu.edu.tw/.

Using ChIPMotifs for de novo motif discovery of OCT4 and ZNF263 based on ChIP-based high-throughput experiments

DNA motifs are short sequences varying from 6 to 25 bp and can be highly variable and degenerated. One major approach for predicting transcription factor (TF) binding is using position weight matrix (PWM) to represent information content of regulatory sites; however, when used as the sole means of identifying binding sites suffers from the limited amount of training data available and a high rate of false-positive predictions. ChIPMotifs program is a de novo motif finding tool developed for ChIP-based high-throughput data, and W-ChIPMotifs is a Web application tool for ChIPMotifs. It composes various ab initio motif discovery tools such as MEME, MaMF, Weeder and optimizes the significance of the detected motifs by using bootstrap re-sampling error estimation and a Fisher test. Using these techniques, we determined a PWM for OCT4 which is similar to canonical OCT4 consensus sequence. In a separate study, we also use de novo motif discovery to suggest that ZNF263 binds to a 24-nt site that differs from the motif predicted by the zinc finger code in several positions.

Modeling community-wide molecular networks of multicellular systems

MOTIVATION

Multicellular systems, such as tissues, are composed of different cell types that form a heterogeneous community. Behavior of these systems is determined by complex regulatory networks within (intracellular networks) and between (intercellular networks) cells. Increasingly more studies are applying genome-wide experimental approaches to delineate the contributions of individual cell types (e.g. stromal, epithelial, vascular cells) to collective behavior of heterogeneous cell communities (e.g. tumors). Although many computational methods have been developed for analyses of intracellular networks based on genome-scale data, these efforts have not been extended toward analyzing genomic data from heterogeneous cell communities.

RESULTS

Here, we propose a network-based approach for analyses of genome-scale data from multiple cell types to extract community-wide molecular networks comprised of intra- and intercellular interactions. Intercellular interactions in this model can be physical interactions between proteins or indirect interactions mediated by secreted metabolites of neighboring cells. Applying this method on data from a recent study on xenograft mouse models of human lung adenocarcinoma, we uncover an extensive network of intra- and intercellular interactions involved in the acquired resistance to angiogenesis inhibitors.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

BINDING MATRIX: A NOVEL APPROACH FOR BINDING SITE RECOGNITION

Recognition of protein-DNA binding sites in genomic sequences is a crucial step for discovering biological functions of genomic sequences. Explosive growth in availability of sequence information has resulted in a demand for binding site detection methods with high specificity. The motivation of the work presented here is to address this demand by a systematic approach based on Maximum Likelihood Estimation.

A general framework is developed in which a large class of binding site detection methods can be described in a uniform and consistent way. Protein-DNA binding is determined by binding energy, which is an approximately linear function within the space of sequence words. All matrix based binding word detectors can be regarded as different linear classifiers which attempt to estimate the linear separation implied by the binding energy function. The standard approaches of consensus sequences and profile matrices are described using this framework.

A maximum likelihood approach for determining this linear separation leads to a novel matrix type, called the binding matrix. The binding matrix is the most specific matrix based classifier which is consistent with the input set of known binding words. It achieves significant improvements in specificity compared to other matrices. This is demonstrated using 95 sets of experimentally determined binding words provided by the TRANSFAC database.

Left ventricular remodeling in swine after myocardial infarction: a transcriptional genomics approach

Despite the apparent appropriateness of left ventricular (LV) remodeling following myocardial infarction (MI), it poses an independent risk factor for development of heart failure. There is a paucity of studies into the molecular mechanisms of LV remodeling in large animal species. We took an unbiased molecular approach to identify candidate transcription factors (TFs) mediating the genetic reprogramming involved in post-MI LV remodeling in swine. Left ventricular tissue was collected from remote, non-infarcted myocardium, 3 weeks after MI-induction or sham-surgery. Microarray analysis identified 285 upregulated and 278 downregulated genes (FDR < 0.05). Of these differentially expressed genes, the promoter regions of the human homologs were searched for common TF binding sites (TFBS). Eighteen TFBS were overrepresented >two-fold (p < 0.01) in upregulated and 13 in downregulated genes. Left ventricular nuclear protein extracts were assayed for DNA-binding activity by protein/DNA array. Out of 345 DNA probes, 30 showed signal intensity changes >two-fold. Five TFs were identified in both TFBS and protein/DNA array analyses, which showed matching changes for COUP-TFII and glucocorticoid receptor (GR) only. Treatment of swine with the GR antagonist mifepristone after MI reduced the post-MI increase in LV mass, but LV dilation remained unaffected. Thus, using an unbiased approach to study post-MI LV remodeling in a physiologically relevant large animal model, we identified COUP-TFII and GR as potential key mediators of post-MI remodeling.

Epigenetic regulation of nuclear PI-PLCbeta1 signaling pathway in low-risk MDS patients during azacitidine treatment

Phosphoinositide-phospholipase C (PI-PLC) beta1 can be considered a specific target for demethylating therapy in high-risk myelodysplastic syndrome (MDS) patients, as azacitidine treatment has been associated with a PI-PLCbeta1-specific promoter demethylation, and induction of PI-PLCbeta1 gene and protein expression. However, little is known about the molecular effect of azacitidine in low-risk MDS or the functional mechanisms linked with azacitidine effect on PI-PLCbeta1 promoter. In the present study, we further investigated the role of epigenetic regulation of PI-PLCbeta1, mainly focusing on the structure of the PI-PLCbeta1 promoter. We first examined the effect of azacitidine on PI-PLCbeta1 promoter methylation and gene expression in low-risk MDS. Moreover, we studied the expression of key molecules associated with the nuclear inositide signaling pathways, such as cyclin D3. By applying a chromatin immunoprecipitation method, we also studied the correlation between the demethylating effect of azacitidine and the degree of recruitment to PI-PLCbeta1 promoter of some transcription factors implicated in hematopoietic stem cell proliferation and differentiation, as well as of the methyl-CpG-binding domain proteins, which specifically interact with methylated DNA. Taken together, our results hint at a specific involvement of PI-PLCbeta1 in epigenetic mechanisms, and are particularly consistent with the hypothesis of a role for PI-PLCbeta1 in azacitidine-induced myeloid differentiation.

A new methodology to associate SNPs with human diseases according to their pathway related context

Genome-wide association studies (GWAS) with hundreds of żthousands of single nucleotide polymorphisms (SNPs) are popular strategies to reveal the genetic basis of human complex diseases. Despite many successes of GWAS, it is well recognized that new analytical approaches have to be integrated to achieve their full potential. Starting with a list of SNPs, found to be associated with disease in GWAS, here we propose a novel methodology to devise functionally important KEGG pathways through the identification of genes within these pathways, where these genes are obtained from SNP analysis. Our methodology is based on functionalization of important SNPs to identify effected genes and disease related pathways. We have tested our methodology on WTCCC Rheumatoid Arthritis (RA) dataset and identified: i) previously known RA related KEGG pathways (e.g., Toll-like receptor signaling, Jak-STAT signaling, Antigen processing, Leukocyte transendothelial migration and MAPK signaling pathways); ii) additional KEGG pathways (e.g., Pathways in cancer, Neurotrophin signaling, Chemokine signaling pathways) as associated with RA. Furthermore, these newly found pathways included genes which are targets of RA-specific drugs. Even though GWAS analysis identifies 14 out of 83 of those drug target genes; newly found functionally important KEGG pathways led to the discovery of 25 out of 83 genes, known to be used as drug targets for the treatment of RA. Among the previously known pathways, we identified additional genes associated with RA (e.g. Antigen processing and presentation, Tight junction). Importantly, within these pathways, the associations between some of these additionally found genes, such as HLA-C, HLA-G, PRKCQ, PRKCZ, TAP1, TAP2 and RA were verified by either OMIM database or by literature retrieved from the NCBI PubMed module. With the whole-genome sequencing on the horizon, we show that the full potential of GWAS can be achieved by integrating pathway and network-oriented analysis and prior knowledge from functional properties of a SNP.

Global mapping of cell type-specific open chromatin by FAIRE-seq reveals the regulatory role of the NFI family in adipocyte differentiation

Identification of regulatory elements within the genome is crucial for understanding the mechanisms that govern cell type-specific gene expression. We generated genome-wide maps of open chromatin sites in 3T3-L1 adipocytes (on day 0 and day 8 of differentiation) and NIH-3T3 fibroblasts using formaldehyde-assisted isolation of regulatory elements coupled with high-throughput sequencing (FAIRE-seq). FAIRE peaks at the promoter were associated with active transcription and histone modifications of H3K4me3 and H3K27ac. Non-promoter FAIRE peaks were characterized by H3K4me1+/me3-, the signature of enhancers, and were largely located in distal regions. The non-promoter FAIRE peaks showed dynamic change during differentiation, while the promoter FAIRE peaks were relatively constant. Functionally, the adipocyte- and preadipocyte-specific non-promoter FAIRE peaks were, respectively, associated with genes up-regulated and down-regulated by differentiation. Genes highly up-regulated during differentiation were associated with multiple clustered adipocyte-specific FAIRE peaks. Among the adipocyte-specific FAIRE peaks, 45.3% and 11.7% overlapped binding sites for, respectively, PPARγ and C/EBPα, the master regulators of adipocyte differentiation. Computational motif analyses of the adipocyte-specific FAIRE peaks revealed enrichment of a binding motif for nuclear family I (NFI) transcription factors. Indeed, ChIP assay showed that NFI occupy the adipocyte-specific FAIRE peaks and/or the PPARγ binding sites near PPARγ, C/EBPα, and aP2 genes. Overexpression of NFIA in 3T3-L1 cells resulted in robust induction of these genes and lipid droplet formation without differentiation stimulus. Overexpression of dominant-negative NFIA or siRNA-mediated knockdown of NFIA or NFIB significantly suppressed both induction of genes and lipid accumulation during differentiation, suggesting a physiological function of these factors in the adipogenic program. Together, our study demonstrates the utility of FAIRE-seq in providing a global view of cell type-specific regulatory elements in the genome and in identifying transcriptional regulators of adipocyte differentiation.

Searching ChIP-seq genomic islands for combinatorial regulatory codes in mouse embryonic stem cells

Background

To facilitate deciphering underlying transcriptional regulatory circuits in mouse embryonic stem (ES) cells, recent ChIP-seq data provided genome-wide binding locations of several key transcription factors (TFs); meanwhile, existing efforts profiled gene expression in ES cells and in their early differentiated state. It has been shown that the gene expression profiles are correlated with the binding of these TFs. However, it remains unclear whether other TFs, referred to as cofactors, participate the gene regulation by collaborating with the ChIP-seq TFs.

Results

Based on our analyses of the ES gene expression profiles and binding sites of potential cofactors in vicinity of the ChIP-seq TF binding locations, we identified a list of co-binding features that show significantly different characteristics between different gene expression patterns (activated or repressed gene expression in ES cells) at a false discovery rate of 10%. Gene classification with a subset of the identified features achieved up to 20% improvement over classification only based on the ChIP-seq TFs. More than 1/3 of reasoned regulatory roles of cofactor candidates involved in these features are supported by existing literatures. Finally, the predicted target genes of the majority candidates present expected expression change in another independent data set, which serves as a supplementary validation of these candidates.

Conclusions

Our results revealed a list of combinatorial genomic features that are significantly associated with gene expression in ES cells, suggesting potential cofactors of the ChIP-seq TFs for gene regulation.

Endogenous myoglobin in breast cancer is hypoxia-inducible by alternative transcription and functions to impair mitochondrial activity: a role in tumor suppression?

Recently, immunohistochemical analysis of myoglobin (MB) in human breast cancer specimens has revealed a surprisingly widespread expression of MB in this nonmuscle context. The positive correlation with hypoxia-inducible factor 2α (HIF-2α) and carbonic anhydrase IX suggested that oxygen regulates myoglobin expression in breast carcinomas. Here, we report that MB mRNA and protein levels are robustly induced by prolonged hypoxia in breast cancer cell lines, in part via HIF-1/2-dependent transactivation. The hypoxia-induced MB mRNA originated from a novel alternative transcription start site 6 kb upstream of the ATG codon. MB regulation in normal and tumor tissue may thus be fundamentally different. Functionally, the knockdown of MB in MDA-MB468 breast cancer cells resulted in an unexpected increase of O(2) uptake and elevated activities of mitochondrial enzymes during hypoxia. Silencing of MB transcription attenuated proliferation rates and motility capacities of hypoxic cancer cells and, surprisingly, also fully oxygenated breast cancer cells. Endogenous MB in cancer cells is apparently involved in controlling oxidative cell energy metabolism, contrary to earlier findings on mouse heart, where the targeted disruption of the Mb gene did not effect myocardial energetics and O(2) consumption. This control function of MB seemingly impacts mitochondria and influences cell proliferation and motility, but it does so in ways not directly related to the facilitated diffusion or storage of O(2). Hypothetically, the mitochondrion-impairing role of MB in hypoxic cancer cells is part of a novel tumor-suppressive function.

A variant in the heart-specific fatty acid transport protein 6 is associated with lower fasting and postprandial TAG, blood pressure and left ventricular hypertrophy

Fatty acid transport protein 6 (FATP6) is primarily expressed in the heart and seems to be involved in cardiac fatty acid uptake. Therefore, we investigated whether a variation in the 5'-untranslated region of the FATP6 gene is associated with features of the metabolic syndrome and signs of myocardial alteration or heart failure. A total of 755 male participants from a Metabolic Intervention Cohort Kiel were genotyped for the FATP6-7T>A polymorphism (rs2526246) and phenotyped for features of the metabolic syndrome. Participants underwent a glucose tolerance test and the postprandial assessment of metabolic variables after a standardised mixed meal. Left ventricular heart function was evaluated in fifty-four participants. Fasting (P = 0·01) and postprandial (P = 0·02) TAG concentrations were significantly lower in AA homozygotes when compared with wild-type carriers. Homozygosity of allele A was associated with significantly lower postprandial insulin concentrations after a glucose load and significantly lower systolic (P = 0·01) and diastolic (P = 0·01) blood pressure values compared with wild-type carriers. Accordingly, left ventricular heart mass was significantly lower in twenty-seven AA homozygotes in comparison with twenty-seven TT homozygotes, matched for BMI (P = 0·04). In conclusion, the effects of the FATP6 polymorphism on TAG are mediated by affluent dietary fat. The FATP6-7T>A polymorphism may protect from traits of the metabolic syndrome and CVD.

Cloning and expression pattern of alkaline phosphatase during the development of Paralichthys olivaceus

Alkaline phosphatases are ubiquitous enzymes involved in many important biological processes. Mammalian tissue-nonspecific alkaline phosphatase has long been thought to feature in embryonic development and bone formation. In this study, an alkaline phosphatase (ALP) gene from Paralichthys olivaceus was identified by rapid amplification of cDNA ends and genome-walking PCR. The ALP gene extends 10,141 bp and contains 11 exons and 10 introns. The open reading frame of the ALP transcript consists of 1,431 bp, which encodes 476 amino acids products named as POALP. An analysis of its secondary and tertiary structure revealed that the POALP was conserved in different species, but one disulfide linkage made it possible to adapt to low-temperature environment. The ALP activity was found to be first detectable in the embryo before hatching. The POALP was distributed ubiquitously in the body of P. olivaceus and was particularly high in the digestive tract. These findings suggest the potential role of POALP in nutrient absorption and transportation. During the pre-metamorphosis (F stage), ALP gene expression is 2.5-folds of that in the pro-metamorphosis (E stage); but in the post-metamorphosis (I stage), it was 1.8-folds of that of pro-metamorphosis. Exogenetic thyroxine (T4) and thiourea (TU) influenced the ALP gene expression significantly during the metamorphosis. Bioinformatics analysis showed that Japanese flounder ALP promoter region contained promoter sequence and putative recognition site for several transcriptional factors, including SREBP-1, SYR, and CdxA. In vitro promoter assays employing EGFP reporter system demonstrated that the promoter of ALP was active.

Activation of the glucocorticoid receptor is associated with poor prognosis in estrogen receptor-negative breast cancer

Estrogen receptor-negative (ER(-)) breast cancers have limited treatment options and are associated with earlier relapses. Because glucocorticoid receptor (GR) signaling initiates antiapoptotic pathways in ER(-) breast cancer cells, we hypothesized that activation of these pathways might be associated with poor prognosis in ER(-) disease. Here we report findings from a genome-wide study of GR transcriptional targets in a premalignant ER(-) cell line model of early breast cancer (MCF10A-Myc) and in primary early-stage ER(-) human tumors. Chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) coupled to time-course expression profiling led us to identify epithelial-to-mesenchymal transition (EMT) pathways as an important aspect associated with GR activation. We validated these findings by carrying out a meta-analysis of primary breast tumor gene expression from 1,378 early-stage breast cancer patients with long-term clinical follow-up, confirming that high levels of GR expression significantly correlated with shorter relapse-free survival in ER(-) patients who were treated or untreated with adjuvant chemotherapy. Notably, in ER(+) breast cancer patients, high levels of GR expression in tumors were significantly associated with better outcome relative to low levels of GR expression. Gene expression analysis revealed that ER(-) tumors expressing high GR levels exhibited differential activation of EMT, cell adhesion, and inflammation pathways. Our findings suggest a direct transcriptional role for GR in determining the outcome of poor-prognosis ER(-) breast cancers.

Characterization of two common 5' polymorphisms in PEX1 and correlation to survival in PEX1 peroxisome biogenesis disorder patients

Background

Mutations in PEX1 are the most common primary cause of Zellweger syndrome. In addition to exonic mutations, deletions and splice site mutations two 5' polymorphisms at c.-137 and c.-53 with a potential influence on PEX1 protein levels have been described in the 5' untranslated region (UTR) of the PEX1 gene.

Methods

We used RACE and in silico promoter prediction analysis to study the 5' UTR of PEX1. We determined the distribution of PEX1 5' polymorphisms in a cohort of 30 Zellweger syndrome patients by standard DNA sequencing. 5' polymorphisms were analysed in relation to the two most common mutations in PEX1 and were incorporated into a novel genotype-phenotype analysis by correlation of three classes of PEX1 mutations with patient survival.

Results

We provide evidence that the polymorphism 137 bp upstream of the ATG codon is not part of the UTR, rendering it a promoter polymorphism. We show that the first, but not the second most common PEX1 mutation arose independently of a specific upstream polymorphic constellation. By genotype-phenotype analysis we identified patients with identical exonic mutation and identical 5' polymorphisms, but strongly differing survival.

Conclusions

Our study suggests that two different types of PEX1 5' polymorphisms have to be distinguished: a 5' UTR polymorphism at position c.-53 and a promoter polymorphism 137 bp upstream of the PEX1 start codon. Our results indicate that the exonic PEX1 mutation correlates with patient survival, but the two 5' polymorphisms analysed in this study do not have to be considered for diagnostic and/or prognostic purposes.

Functional analysis of the DAT gene promoter using transient Catharanthus roseus and stable Nicotiana tabacum transformation systems

The Catharanthus roseus DAT gene encodes the enzyme acetyl-CoA:deacetylvindoline-4-O-acetyltransferase involved in the last step of the indole alkaloid pathway leading to vindoline. This gene is characterized by specific cell type expression in idioblasts and laticifers. To understand the specific transcriptional regulation mechanism(s) of DAT, several DAT promoter GUS constructs were cloned into pCAMBIA1305.1. Agroinfiltration of different explant types of C. roseus resulted in organ-specific accumulation of GUS, albeit at various levels. Heterologous accumulation of GUS in transgenic tobacco revealed both general and non-specific expression with the exception of a stomata-specific expression when 2.3 kb of the DAT promoter was coupled with a portion of the DAT ORF. These results suggest that in addition to the 2.3 kb upstream of the DAT transcriptional start site, additional cis-acting elements may be responsible for the specific spatial expression of DAT in vivo. Furthermore, hairy roots transformed with DAT promoter GUS constructs demonstrated GUS expression in root tissues (visualized through GUS enzyme activity), even though DAT is repressed in non-transformed roots.

Polymorphisms in promoter sequences of MDM2, p53, and p16 genes in normal Japanese individuals

Research has been conducted to identify sequence polymorphisms of gene promoter regions in patients and control subjects, including normal individuals, and to determine the influence of these polymorphisms on transcriptional regulation in cells that express wild-type or mutant p53. In this study we isolated genomic DNA from whole blood of healthy Japanese individuals and sequenced the promoter regions of the MDM2, p53, and p16(INK4a) genes. We identified polymorphisms comprising 3 nucleotide substitutions at exon 1 and intron 1 regions of the MDM2 gene and 1 nucleotide insertion at a poly(C) nucleotide position in the p53 gene. The Japanese individuals also exhibited p16(INK4a) polymorphisms at several positions, including position -191. Reporter gene analysis by using luciferase revealed that the polymorphisms of MDM2, p53, and p16(INK4a) differentially altered luciferase activities in several cell lines, including the Colo320DM, U251, and T98G cell lines expressing mutant p53. Our results indicate that the promoter sequences of these genes differ among normal Japanese individuals and that polymorphisms can alter gene transcription activity.