Discovering transcriptional modules by Bayesian data integration.

MOTIVATION: We present a method for directly inferring transcriptional modules (TMs) by integrating gene expression and transcription factor binding (ChIP-chip) data. Our model extends a hierarchical Dirichlet process mixture model to allow data fusion on a gene-by-gene basis. This encodes the intuition that co-expression and co-regulation are not necessarily equivalent and hence we do not expect all genes to group similarly in both datasets. In particular, it allows us to identify the subset of genes that share the same structure of transcriptional modules in both datasets. RESULTS: We find that by working on a gene-by-gene basis, our model is able to extract clusters with greater functional coherence than existing methods. By combining gene expression and transcription factor binding (ChIP-chip) data in this way, we are better able to determine the groups of genes that are most likely to represent underlying TMs. AVAILABILITY: If interested in the code for the work presented in this article, please contact the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Global transcriptional responses of the toxic cyanobacterium, Microcystis aeruginosa, to nitrogen stress, phosphorus stress, and growth on organic matter

Whole transcriptome shotgun sequencing (RNA-seq) was used to assess the transcriptomic response of the toxic cyanobacterium Microcystis aeruginosa during growth with low levels of dissolved inorganic nitrogen (low N), low levels of dissolved inorganic phosphorus (low P), and in the presence of high levels of high molecular weight dissolved organic matter (HMWDOM). Under low N, one third of the genome was differentially expressed, with significant increases in transcripts observed among genes within the nir operon, urea transport genes (urtBCDE), and amino acid transporters while significant decreases in transcripts were observed in genes related to photosynthesis. There was also a significant decrease in the transcription of the microcystin synthetase gene set under low N and a significant decrease in microcystin content per Microcystis cell demonstrating that N supply influences cellular toxicity. Under low P, 27% of the genome was differentially expressed. The Pho regulon was induced leading to large increases in transcript levels of the alkaline phosphatase phoX, the Pst transport system (pstABC), and the sphX gene, and transcripts of multiple sulfate transporter were also significantly more abundant. While the transcriptional response to growth on HMWDOM was smaller (5–22% of genes differentially expressed), transcripts of multiple genes specifically associated with the transport and degradation of organic compounds were significantly more abundant within HMWDOM treatments and thus may be recruited by Microcystis to utilize these substrates. Collectively, these findings provide a comprehensive understanding of the nutritional physiology of this toxic, bloom-forming cyanobacterium and the role of N in controlling microcystin synthesis.

Evolution of conserved secondary structures and their function in transcriptional regulation networks

Background: Many conserved secondary structures have been identified within conserved elements in the human genome, but only a small fraction of them are known to be functional RNAs. The evolutionary variations of these conserved secondary structures in h

Statistical analysis of sequence features of introns with positive transcriptional regulation in yeast genes

A great deal of experimental studies have shown that many introns of eukaryotic genes function as regulators of transcription. However, comprehensive studies of this problem have not yet been conducted. After checking the transcription frequencies of some Saccharomyces cerevisiae (yeast), genes and their introns, a remarkable phenomenon was discovered that generally the introns of the genes with higher transcription frequencies are longer, and the introns of the genes with lower transcription frequencies are shorter. This suggests that the longer introns of genes with higher transcription frequencies may contain some characteristic sequence structures, which could enhance the transcription of genes. Therefore, two sets of introns of yeast genes were chosen for further study. The transcription frequencies of the first set of genes are higher (>30), and those of the second set of genes are lower (less than or equal to10). Some oligonucleotides are detected by statistically comparative analyses of the occurrence frequencies of oligonucleotides (mainly tetranucleotides and pentanucleotides), whose occurrence frequencies in the first set of introns; are significantly higher than those in the second set of introns, and are also significantly higher than those in the exons flanking the introns of the first set. Some of these extracted oligonucleotides are the same as the regulatory elements of transcription revealed by experimental analyses. Besides, the distributions of these extracted oligonucleotides in the two sets of introns and the exons show that the sequence structures of the first set of introns are favorable for transcription of genes.

The identification and quantification of highly stable 'common hairpin' in the dynamic process of co-transcriptional mRNA folding

In our studies, 88 human mRNA samples were collected from the Integrated Sequence-Structure database and then the dynamic process in co-transcriptional mRNA folding was simulated using the RNAstructure version 4.1 program. Through statistical analyses of the frequencies of occurrence of hairpins, a group of special folding structures-the 'common hairpins'-were identified. These 'common hairpins' have lower energies and occur in all the subsequent folding units that formed in the dynamic folding process. By applying the formulas (1)-(4) of the 'common hairpins' statistical model, 163 'common hairpins' were found, to make up about 7% of the total of 2286 hairpins. Classified studies further show that the 'common hairpins' that were studied may oscillate in the dynamic folding process. However, the hairpin loops of the 'common hairpins' and stems proximal to those 'common hairpins' loops maintain topologically stable structures, while other loops and stems distal to the 'common hairpins' loops are shown to be alterable structures. Strikingly, further studies indicate that the stable structures of these 'common hairpins' may have unbeknown effects on controlling the formation of protein structures in the translation process (unpublished results). (c) 2005 Elsevier B.V. All rights reserved.

Distribution of IgM, IgD and IgZ in mandarin fish, Siniperca chuatsi lymphoid tissues and their transcriptional changes after Flavobacterium columnare stimulation

The gene sequences of three different immunoglobulin (Ig) heavy chains, namely IgM, IgD and IgZ, were cloned from mandarin fish (Siniperca chuatsi) recently. In this study the distribution of these three kinds of Ig-producing cells in lymphoid-related tissues as head kidney, spleen, gill and intestine were investigated by using in situ hybridization, and their transcriptional changes were also analyzed by quantitative real-time PCR during 8 weeks after immunization. IgM-producing cells could be detected obviously and abundantly in all the tissues examined. A few numbers of IgD and IgZ positive cells were both detected in head kidney and spleen. IgZ positive cells could be detected in gill moderately while IgD showed negative results, otherwise no IgD or IgZ positive cells could be detected in intestine. After stimulated with bacterial pathogen Flavobacterium columnare G(4), the transcripts of these three Ig genes exhibited quite different kinetics. Significantly increased transcription of IgM gene was observed in almost all the tissues examined especially in boosted group. In contrast with IgM, seldom strong increase was examined for IgD and IgZ genes. For IgD, it seemed that the first injection could stimulate the immune response easier, since in almost all the tissues significant increase was detected at 1 or 2 weeks after injection. For IgZ, boosted injection could not enlarge the up-regulation of gene expression of first injection. This is the first case to report the transcriptional kinetics of three Ig genes in teleost after bacterin immunization. (C) 2008 Elsevier B.V. All rights reserved.

Genomic sequence of mandarin fish rhabdovirus with an unusual small non-transcriptional ORF

The complete genome of mandarin fish Siniperca chuatsi rhabdovirus (SCRV) was cloned and sequenced. It comprises 11,545 nucleotides and contains five genes encoding the nucleoprotein N, the phosphoprotein P, the matrix protein M, the glycoprotein G, and the RNA-dependent RNA polymerase protein L. At the 3' and 5' termini of SCRV genome, leader and trailer sequences show inverse complementarity. The N, P, M and G proteins share the highest sequence identities (ranging from 14.8 to 41.5%) with the respective proteins of rhabdovirus 903/87, the L protein has the highest identity with those of vesiculoviruses, especially with Chandipura virus (44.7%). Phylogenetic analysis of L proteins showed that SCRV clustered with spring vireamia of carp virus (SVCV) and was most closely related to viruses in the genus Vesiculovirus. In addition, an overlapping open reading frame (ORF) predicted to encode a protein similar to vesicular stomatitis virus C protein is present within the P gene of SCRV. Furthermore, an unoverlapping small ORF downstream of M ORF within M gene is predicted (tentatively called orf4). Therefore, the genomic organization of SCRV can be proposed as 3' leader-N-P/C-M-(orf4)-G-L-trailer 5'. Orf4 transcription or translation products could not be detected by northern or Western blot, respectively, though one similar mRNA band to M mRNA was found. This is the first report on one small unoverlapping ORF in M gene of a fish rhabdovirus. (c) 2007 Elsevier B.V. All rights reserved.

FB1, an E2A fusion partner in childhood leukemia, interacts with U19/EAF2 and inhibits its transcriptional activity

Background: U19/EAF2 is a potential tumor suppressor exhibiting frequent down-regulation and allelic loss in advanced human prostate cancer specimens. U 19/EAF2 has also been identified as ELL-associated factor 2 (EAF2) based on its binding to ELL, a fusion partner of MLL in acute myeloid leukemia. U19/EAF2 is a putative transcription factor with a transactivation domain and capability of sequence-specific DNA binding. Methods: Yeast-two-hybrid-screening was used to identify U19/EAF2-binding partners. Co-immunoprecipitation and mammalian 1-hybrid assay were used to characterize a U19/EAF2-binding partner. Results: FB1, an E2A fusion partner in childhood leukemia, was identified as a binding-partner of U19/EAF2. FB1 also binds to EAF1, the only homologue of U19/EAF2. FB1 also interacts and co-localizes with ELL in the nucleus. Interestingly, FB1 inhibited the transcriptional activity of U19/EAF2 but not EAF1. Conclusions: FB1 is an important binding partner and a functional regulator of U19/EAF2, EAF1, and/or ELL. (c) 2007 Elsevier Ireland Ltd. All rights reserved.

A host transcriptional signature for presymptomatic detection of infection in humans exposed to influenza H1N1 or H3N2.

There is great potential for host-based gene expression analysis to impact the early diagnosis of infectious diseases. In particular, the influenza pandemic of 2009 highlighted the challenges and limitations of traditional pathogen-based testing for suspected upper respiratory viral infection. We inoculated human volunteers with either influenza A (A/Brisbane/59/2007 (H1N1) or A/Wisconsin/67/2005 (H3N2)), and assayed the peripheral blood transcriptome every 8 hours for 7 days. Of 41 inoculated volunteers, 18 (44%) developed symptomatic infection. Using unbiased sparse latent factor regression analysis, we generated a gene signature (or factor) for symptomatic influenza capable of detecting 94% of infected cases. This gene signature is detectable as early as 29 hours post-exposure and achieves maximal accuracy on average 43 hours (p = 0.003, H1N1) and 38 hours (p-value = 0.005, H3N2) before peak clinical symptoms. In order to test the relevance of these findings in naturally acquired disease, a composite influenza A signature built from these challenge studies was applied to Emergency Department patients where it discriminates between swine-origin influenza A/H1N1 (2009) infected and non-infected individuals with 92% accuracy. The host genomic response to Influenza infection is robust and may provide the means for detection before typical clinical symptoms are apparent.

The Ets transcription factors of the PEA3 group: Transcriptional regulators in metastasis

The PEA3 group is composed of three highly conserved Ets transcription factors: Erm, Er81, and Pea3. These proteins regulate transcription of multiple genes, and their transactivating potential is affected by post-translational modifications. Among their target genes are several matrix metalloproteases (MMPs), which are enzymes degrading the extracellular matrix during normal remodelling events and cancer metastasis. In fact, PEA3-group genes are often over-expressed in different types of cancers that also over-express these MMPs and display a disseminating phenotype. Experimental regulation of the synthesis of PEA3 group members influences the metastatic process. This suggests that these factors play a key role in metastasis. © 2006 Elsevier B.V. All rights reserved.

Ets-1 p27: A novel Ets-1 isoform with dominant-negative effects on the transcriptional properties and the subcellular localization of Ets-1 p51

The transcription factor Ets-1 is implicated in various physiological processes and invasive pathologies. We identified a novel variant of ets-1, ets-1Delta(III-VI), resulting from the alternative splicing of exons III to VI. This variant encodes a 27 kDa isoform, named Ets-1 p27. Ets-1 p27 lacks the threonine-38 residue, the Pointed domain and the transactivation domain, all of which are required for the transactivation of Ets-1 target genes. Both inhibitory domains surrounding the DNA-binding domain are conserved, suggesting that Ets-1 p27, like the full-length Ets-1 p51 isoform, is autoinhibited for DNA binding. We showed that Ets-1 p27 binds DNA in the same way as Ets-1 p51 does and that it acts both at a transcriptional and a subcellular localization level, thereby constituting a dual-acting dominant negative of Ets-1 p51. Ets-1 p27 blocks Ets-1 p51-mediated transactivation of target genes and induces the translocation of Ets-1 p51 from the nucleus to the cytoplasm. Furthermore, Ets-1 p27 overexpression represses the tumor properties of MDA-MB-231 mammary carcinoma cells in correlation with the known implication of Ets-1 in various cellular mechanisms. Thus the dual-acting dominant-negative function of Ets-1 p27 gives to the Ets-1 p27/Ets-1 p51 ratio a determining effect on cell fate.