EGASP: the human ENCODE Genome Annotation Assessment Project

Background: Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes. Results: Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families. Conclusions: The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

SA2RAGE: a new sequence for fast B1+ -mapping.

At high magnetic field strengths (≥ 3T), the radiofrequency wavelength used in MRI is of the same order of magnitude of (or smaller than) the typical sample size, making transmit magnetic field (B1+) inhomogeneities more prominent. Methods such as radiofrequency-shimming and transmit SENSE have been proposed to mitigate these undesirable effects. A prerequisite for such approaches is an accurate and rapid characterization of the B1+ field in the organ of interest. In this work, a new phase-sensitive three-dimensional B1+-mapping technique is introduced that allows the acquisition of a 64 × 64 × 8 B1+-map in ≈ 20 s, yielding an accurate mapping of the relative B1+ with a 10-fold dynamic range (0.2-2 times the nominal B1+). Moreover, the predominant use of low flip angle excitations in the presented sequence minimizes specific absorption rate, which is an important asset for in vivo B1+-shimming procedures at high magnetic fields. The proposed methodology was validated in phantom experiments and demonstrated good results in phantom and human B1+-shimming using an 8-channel transmit-receive array.

Similarities and differences of polyadenylation signals in human and fly.

BACKGROUND: Cleavage of messenger RNA (mRNA) precursors is an essential step in mRNA maturation. The signal recognized by the cleavage enzyme complex has been characterized as an A rich region upstream of the cleavage site containing a motif with consensus AAUAAA, followed by a U or UG rich region downstream of the cleavage site. RESULTS: We studied these signals using exhaustive databases of cleavage sites obtained from aligning raw expressed sequence tags (EST) sequences to genomic sequences in Homo sapiens and Drosophila melanogaster. These data show that the polyadenylation signal is highly conserved in human and fly. In addition, de novo motif searches generated a refined description of the U-rich downstream sequence (DSE) element, which shows more divergence between the two species. These refined motifs are applied, within a Hidden Markov Model (HMM) framework, to predict mRNA cleavage sites. CONCLUSION: We demonstrate that the DSE is a specific motif in both human and Drosophila. These findings shed light on the sequence correlates of a highly conserved biological process, and improve in silico prediction of 3' mRNA cleavage and polyadenylation sites.

Short-term overexpression of a constitutively active form of AMP-activated protein kinase in the liver leads to mild hypoglycemia and fatty liver.

AMP-activated protein kinase (AMPK) is a major therapeutic target for the treatment of diabetes. We investigated the effect of a short-term overexpression of AMPK specifically in the liver by adenovirus-mediated transfer of a gene encoding a constitutively active form of AMPKalpha2 (AMPKalpha2-CA). Hepatic AMPKalpha2-CA expression significantly decreased blood glucose levels and gluconeogenic gene expression. Hepatic expression of AMPKalpha2-CA in streptozotocin-induced and ob/ob diabetic mice abolished hyperglycemia and decreased gluconeogenic gene expression. In normal mouse liver, AMPKalpha2-CA considerably decreased the refeeding-induced transcriptional activation of genes encoding proteins involved in glycolysis and lipogenesis and their upstream regulators, SREBP-1 (sterol regulatory element-binding protein-1) and ChREBP (carbohydrate response element-binding protein). This resulted in decreases in hepatic glycogen synthesis and circulating lipid levels. Surprisingly, despite the inhibition of hepatic lipogenesis, expression of AMPKalpha2-CA led to fatty liver due to the accumulation of lipids released from adipose tissue. The relative scarcity of glucose due to AMPKalpha2-CA expression led to an increase in hepatic fatty acid oxidation and ketone bodies production as an alternative source of energy for peripheral tissues. Thus, short-term AMPK activation in the liver reduces blood glucose levels and results in a switch from glucose to fatty acid utilization to supply energy needs.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control.

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10(-12)). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the 'intermediate phenotype' nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host-pathogen interaction. DOI:http://dx.doi.org/10.7554/eLife.01123.001.

Mapping the transcription start points of the Staphylococcus aureus eap, emp, and vwb promoters reveals a conserved octanucleotide sequence that is essential for expression of these genes.

Mapping the transcription start points of the eap, emp, and vwb promoters revealed a conserved octanucleotide sequence (COS). Deleting this sequence abolished the expression of eap, emp, and vwb. However, electrophoretic mobility shift assays gave no evidence that this sequence was a binding site for SarA or SaeR, known regulators of eap and emp.

Induction of the acyl-coenzyme A synthetase gene by fibrates and fatty acids is mediated by a peroxisome proliferator response element in the C promoter.

The long-chain acyl-coenzyme A synthetase (ACS) gene gives rise to three transcripts containing different first exons preceded by specific regulatory regions A, B, and C. Exon-specific oligonucleotide hybridization indicated that only A-ACS mRNA is expressed in rat liver. Fibrate administration induced liver C-ACS strongly and A-ACS mRNA to a lesser extent. B-ACS mRNA remained undetectable. In primary rat hepatocytes and Fa-32 hepatoma cells C-ACS mRNA increased after treatment with fenofibric acid, alpha-bromopalmitate, tetradecylthioacetic acid, or alpha-linolenic acid. Nuclear run-on experiments indicated that fenofibric acid and alpha-bromopalmitate act at the transcriptional level. Transient transfections showed a 3.4-, 2.3-, and 2.2-fold induction of C-ACS promoter activity after fenofibric acid, alpha-bromopalmitate, and tetradecylthioacetic acid, respectively. Unilateral deletion and site-directed mutagenesis identified a peroxisome proliferator activator receptor (PPAR)-responsive element (PPRE) mediating the responsiveness to fibrates and fatty acids. This ACS PPRE contains three imperfect half sites spaced by 1 and 3 oligonucleotides and binds PPAR.retinoid X receptor heterodimers in gel retardation assays. In conclusion, the regulation of C-ACS mRNA expression by fibrates and fatty acids is mediated by PPAR.retinoid X receptor heterodimers interacting through a PPRE in the C-ACS promoters. PPAR therefore occupies a key position in the transcriptional control of a pivotal enzyme controlling the channeling of fatty acids into various metabolic pathways.

Changes of the mucosal N3 and N6 fatty acid status occur early in the colorectal adenoma-carcinoma sequence

Despite data favouring a role of dietary fat in colonic carcinogenesis, no study has focused on tissue n3 and n6 fatty acid (FA) status in human colon adenoma-carcinoma sequence. Thus, FA profile was measured in plasma phospholipids of patients with colorectal cancer (n = 22), sporadic adenoma (n = 27), and normal colon (n = 12) (control group). Additionally, mucosal FAs were assessed in both diseased and normal mucosa of cancer (n = 15) and adenoma (n = 21) patients, and from normal mucosa of controls (n = 8). There were no differences in FA profile of both plasma phospholipids and normal mucosa, between adenoma and control patients. There were considerable differences, however, in FAs between diseased and paired normal mucosa of adenoma patients, with increases of linoleic (p = 0.02), dihomogammalinolenic (p = 0.014), and eicosapentaenoic (p = 0.012) acids, and decreases of alpha linolenic (p = 0.001) and arachidonic (p = 0.02) acids in diseased mucosa. A stepwise reduction of eicosapentaenoic acid concentrations in diseased mucosa from benign adenoma to the most advanced colon cancer was seen (p = 0.009). Cancer patients showed lower alpha linolenate (p = 0.002) and higher dihomogammalinolenate (p = 0.003) in diseased than in paired normal mucosa. In conclusion changes in tissue n3 and n6 FA status might participate in the early phases of the human colorectal carcinogenesis.

Epigenetic regulation of histone H3 serine 10 phosphorylation status by HCF-1 proteins in C. elegans and mammalian cells.

BACKGROUND: The human herpes simplex virus (HSV) host cell factor HCF-1 is a transcriptional coregulator that associates with both histone methyl- and acetyltransferases, and a histone deacetylase and regulates cell proliferation and division. In HSV-infected cells, HCF-1 associates with the viral protein VP16 to promote formation of a multiprotein-DNA transcriptional activator complex. The ability of HCF proteins to stabilize this VP16-induced complex has been conserved in diverse animal species including Drosophila melanogaster and Caenorhabditis elegans suggesting that VP16 targets a conserved cellular function of HCF-1. METHODOLOGY/PRINCIPAL FINDINGS: To investigate the role of HCF proteins in animal development, we have characterized the effects of loss of the HCF-1 homolog in C. elegans, called Ce HCF-1. Two large hcf-1 deletion mutants (pk924 and ok559) are viable but display reduced fertility. Loss of Ce HCF-1 protein at reduced temperatures (e.g., 12 degrees C), however, leads to a high incidence of embryonic lethality and early embryonic mitotic and cytokinetic defects reminiscent of mammalian cell-division defects upon loss of HCF-1 function. Even when viable, however, at normal temperature, mutant embryos display reduced levels of phospho-histone H3 serine 10 (H3S10P), a modification implicated in both transcriptional and mitotic regulation. Mammalian cells with defective HCF-1 also display defects in mitotic H3S10P status. CONCLUSIONS/SIGNIFICANCE: These results suggest that HCF-1 proteins possess conserved roles in the regulation of cell division and mitotic histone phosphorylation.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Sequence homologies within the 5' end region of the estrogen-controlled vitellogenin gene in Xenopus and chicken.

In oviparous vertebrates vitellogenin, the precursor of the major yolk proteins, is synthesized in the liver of mature females under the control of estrogen. We have established the organization and primary structure of the 5' end region of the Xenopus laevis vitellogenin A2 gene and of the major chicken vitellogenin gene. The first three homologous exons have exactly the same length in both species, namely 53, 21 and 152 nucleotides, and present an overall sequence homology of 60%. In both species, the 5'-non-coding region of the vitellogenin mRNA measures only 13 nucleotides, nine of which are conserved. In contrast, the corresponding introns of the Xenopus and the chicken vitellogenin gene show no significant sequence homology. Within the 500 nucleotides preceding the 5' end of the genes, at least six blocks with sequence homologies of greater than 70% were detected. It remains to be demonstrated which of these conserved sequences, if any, are involved in the hormone-regulated expression of the vitellogenin genes.

Noise-enhanced excitability in bistable activator-inhibitor media

We show that external fluctuations induce excitable behavior in a bistable spatially extended system with activator-inhibitor dynamics of the FitzHugh-Nagumo type. This can be understood as a mechanism for sustained signal propagation in bistable media. The phase diagram of the stochastic system is analytically obtained and numerically verified. For small-noise intensities, front propagation becomes unstable, and excitable pulses arise as the only possible spatiotemporal behavior of the system. For large-noise intensities, on the other hand, the system enters an effective regime of oscillatory behavior, where it exhibits spontaneous nucleation of pulses and synchronized firing.

aIr-1, a newly found locus on mouse chromosome 16 encoding a trans-acting activator factor for MHC class II gene expression.

RJ 2.2.5 is a human B cell line that has lost the capacity to express MHC class II genes. The human class II-positive phenotype is restored in somatic cell hybrids between RJ 2.2.5 and mouse spleen cells. By karyotype and molecular studies of an informative family of hybrids we have now shown that the reexpression of human class II gene products, as well as the maintenance of the mouse class II-positive phenotype, correlates with the presence of mouse chromosome 16. Thus, the existence on this mouse chromosome of a newly found locus, designated by us aIr-1, that determines a trans-acting activator function for class II gene expression, is established. Possible implications of this finding are discussed.