Glucosylated hydroxymethyluracil, DNA base J, prevents transcriptional readthrough in Leishmania.

Some Ts in nuclear DNA of trypanosomes and Leishmania are hydroxylated and glucosylated to yield base J (ß-D-glucosyl-hydroxymethyluracil). In Leishmania, about 99% of J is located in telomeric repeats. We show here that most of the remaining J is located at chromosome-internal RNA polymerase II termination sites. This internal J and telomeric J can be reduced by a knockout of J-binding protein 2 (JBP2), an enzyme involved in the first step of J biosynthesis. J levels are further reduced by growing Leishmania JBP2 knockout cells in BrdU-containing medium, resulting in cell death. The loss of internal J in JBP2 knockout cells is accompanied by massive readthrough at RNA polymerase II termination sites. The readthrough varies between transcription units but may extend over 100 kb. We conclude that J is required for proper transcription termination and infer that the absence of internal J kills Leishmania by massive readthrough of transcriptional stops.

Genetic analysis of Leishmania donovani tropism using a naturally attenuated cutaneous strain.

A central question in Leishmania research is why most species cause cutaneous infections but others cause fatal visceral disease. Interestingly, L. donovani causes both visceral and cutaneous leishmaniasis in Sri Lanka. L. donovani clinical isolates were therefore obtained from cutaneous leishmaniasis (CL-SL) and visceral leishmaniasis (VL-SL) patients from Sri Lanka. The CL-SL isolate was severely attenuated compared to the VL-SL isolate for survival in visceral organs in BALB/c mice. Genomic and transcriptomic analysis argue that gene deletions or pseudogenes specific to CL-SL are not responsible for the difference in disease tropism and that single nucleotide polymorphisms (SNPs) and/or gene copy number variations play a major role in altered pathology. This is illustrated through the observations within showing that a decreased copy number of the A2 gene family and a mutation in the ras-like RagC GTPase enzyme in the mTOR pathway contribute to the attenuation of the CL-SL strain in visceral infection. Overall, this research provides a unique perspective on genetic differences associated with diverse pathologies caused by Leishmania infection.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Overlapping repressor binding sites regulate expression of the Methanococcus maripaludis glnK(1) operon.

The euryarchaeal transcriptional repressor NrpR regulates a variety of nitrogen assimilation genes by 2-oxoglutarate-reversible binding to conserved palindromic operators. The number and positioning of these operators varies among promoter regions of regulated genes, suggesting NrpR can bind in different patterns. Particularly intriguing is the contrast between the nif and glnK(1) promoter regions of Methanococcus maripaludis, where two operators are present but with different configurations. Here we study NrpR binding and regulation at the glnK(1) promoter, where the two operator sequences overlap and occur on opposite faces of the double helix. We find that both operators function in binding, with a dimer of NrpR binding simultaneously to each overlapping operator. We show in vivo that the first operator plays a primary role in regulation and the second operator plays an enhancing role. This is the first demonstration of overlapping operators functioning in Archaea.

RNA-seq approaches for determining mRNA abundance in Leishmania.

High-throughput sequencing of cDNA copies of mRNA (RNA-seq) provides a digital read-out of mRNA levels over several orders of magnitude, as well as mapping the transcripts to the nucleotide level. Here we describe an RNA-seq approach that exploits the 39-nucleotide mini-exon or spliced leader (SL) sequence found at the 5' end of all Leishmania (and other trypanosomatid) mRNAs.

Continuous culture of Methanococcus maripaludis under defined nutrient conditions.

To study global regulation in the methanogenic archaeon Methanococcus maripaludis, we devised a system for steady-state growth in chemostats. New Brunswick Bioflo 110 bioreactors were equipped with controlled delivery of hydrogen, nitrogen, carbon dioxide, hydrogen sulfide, and anaerobic medium. We determined conditions and media compositions for growth with three different limiting nutrients, hydrogen, phosphate, and leucine. To investigate leucine limitation we constructed and characterized a mutant in the leuA gene for 2-isopropylmalate synthase, demonstrating for the first time the function of this gene in the Archaea. Steady state specific growth rates in these studies ranged from 0.042 to 0.24 h(-1). Plots of culture density vs. growth rate for each condition showed the behavior predicted by growth modeling. The results show that growth behavior is normal and reproducible and validate the use of the chemostat system for metabolic and global regulation studies in M. maripaludis.

Iron uptake controls the generation of Leishmania infective forms through regulation of ROS levels.

During its life cycle, Leishmania undergoes extreme environmental changes, alternating between insect vectors and vertebrate hosts. Elevated temperature and decreased pH, conditions encountered after macrophage invasion, can induce axenic differentiation of avirulent promastigotes into virulent amastigotes. Here we show that iron uptake is a major trigger for the differentiation of Leishmania amazonensis amastigotes, independently of temperature and pH changes. We found that iron depletion from the culture medium triggered expression of the ferrous iron transporter LIT1 (Leishmania iron transporter 1), an increase in iron content of the parasites, growth arrest, and differentiation of wild-type (WT) promastigotes into infective amastigotes. In contrast, LIT1-null promastigotes showed reduced intracellular iron content and sustained growth in iron-poor media, followed by cell death. LIT1 up-regulation also increased iron superoxide dismutase (FeSOD) activity in WT but not in LIT1-null parasites. Notably, the superoxide-generating drug menadione or H(2)O(2) was sufficient to trigger differentiation of WT promastigotes into fully infective amastigotes. LIT1-null promastigotes accumulated superoxide radicals and initiated amastigote differentiation after exposure to H(2)O(2) but not to menadione. Our results reveal a novel role for FeSOD activity and reactive oxygen species in orchestrating the differentiation of virulent Leishmania amastigotes in a process regulated by iron availability.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells.

BACKGROUND: Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian. RESULTS: We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types. CONCLUSIONS: The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Functionally distinct genes regulated by hydrogen limitation and growth rate in methanogenic Archaea.

The use of molecular hydrogen as electron donor for energy generation is a defining characteristic of the hydrogenotrophic methanogens, an ancient group that dominates the phylum Euryarchaeota. We present here a global study of changes in mRNA abundance in response to hydrogen availability for a hydrogenotrophic methanogen. Cells of Methanococcus maripaludis were grown by using continuous culture to deconvolute the effects of hydrogen limitation and growth rate, and microarray analyses were conducted. Hydrogen limitation markedly increased mRNA levels for genes encoding enzymes of the methanogenic pathway that reduce or oxidize the electron-carrying deazaflavin, coenzyme F(420). F(420)-dependent redox functions in energy-generating metabolism are characteristic of the methanogenic Archaea, and the results show that their regulation is distinct from other redox processes in the cell. Rapid growth increased mRNA levels of the gene for an unusual hydrogenase, the hydrogen-dependent methylenetetrahydromethanopterin dehydrogenase.

Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays.

Localized accessibility of critical DNA sequences to the regulatory machinery is a key requirement for regulation of human genes. Here we describe a high-resolution, genome-scale approach for quantifying chromatin accessibility by measuring DNase I sensitivity as a continuous function of genome position using tiling DNA microarrays (DNase-array). We demonstrate this approach across 1% ( approximately 30 Mb) of the human genome, wherein we localized 2,690 classical DNase I hypersensitive sites with high sensitivity and specificity, and also mapped larger-scale patterns of chromatin architecture. DNase I hypersensitive sites exhibit marked aggregation around transcriptional start sites (TSSs), though the majority mark nonpromoter functional elements. We also developed a computational approach for visualizing higher-order features of chromatin structure. This revealed that human chromatin organization is dominated by large (100-500 kb) 'superclusters' of DNase I hypersensitive sites, which encompass both gene-rich and gene-poor regions. DNase-array is a powerful and straightforward approach for systematic exposition of the cis-regulatory architecture of complex genomes.

Function and regulation of the formate dehydrogenase genes of the methanogenic archaeon Methanococcus maripaludis.

Methanococcus maripaludis is a mesophilic species of Archaea capable of producing methane from two substrates: hydrogen plus carbon dioxide and formate. To study the latter, we identified the formate dehydrogenase genes of M. maripaludis and found that the genome contains two gene clusters important for formate utilization. Phylogenetic analysis suggested that the two formate dehydrogenase gene sets arose from duplication events within the methanococcal lineage. The first gene cluster encodes homologs of formate dehydrogenase alpha (FdhA) and beta (FdhB) subunits and a putative formate transporter (FdhC) as well as a carbonic anhydrase analog. The second gene cluster encodes only FdhA and FdhB homologs. Mutants lacking either fdhA gene exhibited a partial growth defect on formate, whereas a double mutant was completely unable to grow on formate as a sole methanogenic substrate. Investigation of fdh gene expression revealed that transcription of both gene clusters is controlled by the presence of H(2) and not by the presence of formate.