Nucleic Acids Research - recent issues

A dual-tag microarray platform for high-performance nucleic acid and protein analyses

2008-05-06

DNA microarrays serve to monitor a wide range of molecular events, but emerging applications like measurements of weakly expressed genes or of proteins and their interaction patterns will require enhanced performance to improve specificity of detection and dynamic range. To further extend the utility of DNA microarray-based approaches we present a high-performance tag microarray procedure that enables probe-based analysis of as little as 100 target cDNA molecules, and with a linear dynamic range close to 10⁵. Furthermore, the protocol radically decreases the risk of cross-hybridization on microarrays compared to current approaches, and it also allows for quantification by single-molecule analysis and real-time on-chip monitoring of rolling-circle amplification. We provide proof of concept for microarray-based measurement of both mRNA molecules and of proteins, converted to tag DNA sequences by padlock and proximity probe ligation, respectively.

Detecting cis-regulatory binding sites for cooperatively binding proteins

2008-05-06

Several methods are available to predict cis-regulatory modules in DNA based on position weight matrices. However, the performance of these methods generally depends on a number of additional parameters that cannot be derived from sequences and are difficult to estimate because they have no physical meaning. As the best way to detect cis-regulatory modules is the way in which the proteins recognize them, we developed a new scoring method that utilizes the underlying physical binding model. This method requires no additional parameter to account for multiple binding sites; and the only necessary parameters to model homotypic cooperative interactions are the distances between adjacent protein binding sites in basepairs, and the corresponding cooperative binding constants. The heterotypic cooperative binding model requires one more parameter per cooperatively binding protein, which is the concentration multiplied by the partition function of this protein. In a case study on the bacterial ferric uptake regulator, we show that our scoring method for homotypic cooperatively binding proteins significantly outperforms other PWM-based methods where biophysical cooperativity is not taken into account.

Exalign: a new method for comparative analysis of exon-intron gene structures

Pavesi, G., Zambelli, F., Caggese, C., Pesole, G. — 2008-05-06

The evolution of genes is usually studied and reconstructed at the sequence level, that is, by comparing and aligning their genomic, transcript or protein sequences. However, including the exon–intron structure of genes in the analysis can provide further and useful information, for example to draw reliable phylogenetic relationships left unsolved by traditional sequence-based evolutionary studies, or to shed further light on patterns of intron gain and loss. In spite of this, no tool especially devised for this task is currently available. In this work we present Exalign, an algorithm designed to retrieve, compare and search for the exon-intron structure of existing gene annotations, that has been implemented in a software tool freely accessible through a web interface as well as available for download. We present different applications of our method, from the reconstruction of the evolutionary history of homologous gene families to the detection of as of today unknown cases of intron loss in human and rodents, and, remarkably, two never reported intron gain events in human and mouse. The web interface for accessing Exalign is available at http://www.pesolelab.it/exalign/ or http://www.beacon.unimi.it/exalign/

Combined in silico and experimental identification of the Pyrococcus abyssi H/ACA sRNAs and their target sites in ribosomal RNAs

2008-05-06

How far do H/ACA sRNPs contribute to rRNA pseudouridylation in Archaea was still an open question. Hence here, by computational search in three Pyrococcus genomes, we identified seven H/ACA sRNAs and predicted their target sites in rRNAs. In parallel, we experimentally identified 17 residues in P. abyssi rRNAs. By in vitro reconstitution of H/ACA sRNPs, we assigned 15 out of the 17 residues to the 7 identified H/ACA sRNAs: one H/ACA motif can guide up to three distinct pseudouridylations. Interestingly, by using a 23S rRNA fragment as the substrate, one of the two remaining residues could be formed in vitro by the aCBF5/aNOP10/aGAR1 complex without guide sRNA. Our results shed light on structural constraints in archaeal H/ACA sRNPs: the length of helix H2 is of 5 or 6 bps, the distance between the ANA motif and the targeted U residue is of 14 or 15 nts, and the stability of the interaction formed by the substrate rRNA and the 3'-guide sequence is more important than that formed with the 5'-guide sequence. Surprisingly, we showed that a sRNA–rRNA interaction with the targeted uridine in a single-stranded 5'-UNN-3' trinucleotide instead of the canonical 5'-UN-3' dinucleotide is functional.

Step-wise formation of eukaryotic double-row polyribosomes and circular translation of polysomal mRNA

2008-05-06

The time course of polysome formation was studied in a long-term wheat germ cell-free translation system using sedimentation and electron microscopy techniques. The polysomes were formed on uncapped luciferase mRNA with translation-enhancing 5' and 3' UTRs. The formation of fully loaded polysomes was found to be a long process that required many rounds of translation and proceeded via several phases. First, short linear polysomes containing no more than six ribosomes were formed. Next, folding of these polysomes into short double-row clusters occurred. Subsequent gradual elongation of the clusters gave rise to heavy-loaded double-row strings containing up to 30–40 ribosomes. The formation of the double-row polysomes was considered to be equivalent to circularization of polysomes, with antiparallel halves of the circle being laterally stuck together by ribosome interactions. A slow exchange with free ribosomes and free mRNA observed in the double-row type polysomes, as well as the resistance of translation in them to AMP-PNP, provided evidence that most polysomal ribosomes reinitiate translation within the circularized polysomes without scanning of 5' UTR, while de novo initiation including 5' UTR scanning proceeds at a much slower rate. Removal or replacements of 5' and 3' UTRs affected the initial phase of translation, but did not prevent the formation of the double-row polysomes during translation.

Two conserved modules of Schizosaccharomyces pombe Mediator regulate distinct cellular pathways

2008-05-06

Mediator is an evolutionary conserved coregulator complex required for transcription of almost all RNA polymerase II-dependent genes. The Schizosaccharomyces pombe Mediator consists of two dissociable components—a core complex organized into a head and middle domain as well as the Cdk8 regulatory subcomplex. In this work we describe a functional characterization of the S. pombe Mediator. We report the identification of the S. pombe Med20 head subunit and the isolation of ts alleles of the core head subunit encoding med17⁺. Biochemical analysis of med8^ts, med17^ts, med18, med20 and med27 alleles revealed a stepwise head domain molecular architecture. Phenotypical analysis of Cdk8 and head module alleles including expression profiling classified the Mediator mutant alleles into one of two groups. Cdk8 module mutants flocculate due to overexpression of adhesive cell-surface proteins. Head domain-associated mutants display a hyphal growth phenotype due to defective expression of factors required for cell separation regulated by transcription factor Ace2. Comparison with Saccharomyces cerevisiae Mediator expression data reveals that these functionally distinct modules are conserved between S. pombe and S. cerevisiae.

Single-molecule manipulation reveals supercoiling-dependent modulation of lac repressor-mediated DNA looping

Normanno, D., Vanzi, F., Pavone, F. S. — 2008-05-06

Gene expression regulation is a fundamental biological process which deploys specific sets of genomic information depending on physiological or environmental conditions. Several transcription factors (including lac repressor, LacI) are present in the cell at very low copy number and increase their local concentration by binding to multiple sites on DNA and looping the intervening sequence. In this work, we employ single-molecule manipulation to experimentally address the role of DNA supercoiling in the dynamics and stability of LacI-mediated DNA looping. We performed measurements over a range of degrees of supercoiling between –0.026 and +0.026, in the absence of axial stretching forces. A supercoiling-dependent modulation of the lifetimes of both the looped and unlooped states was observed. Our experiments also provide evidence for multiple structural conformations of the LacI–DNA complex, depending on torsional constraints. The supercoiling-dependent modulation demonstrated here adds an important element to the model of the lac operon. In fact, the complex network of proteins acting on the DNA in a living cell constantly modifies its topological and mechanical properties: our observations demonstrate the possibility of establishing a signaling pathway from factors affecting DNA supercoiling to transcription factors responsible for the regulation of specific sets of genes.

Evolution of acceptor stem tRNA recognition by class II prolyl-tRNA synthetase

An, S., Barany, G., Musier-Forsyth, K. — 2008-05-06

Aminoacyl-tRNA synthetases (AARS) are an essential family of enzymes that catalyze the attachment of amino acids to specific tRNAs during translation. Previously, we showed that base-specific recognition of the tRNA^Pro acceptor stem is critical for recognition by Escherichia coli prolyl-tRNA synthetase (ProRS), but not for human ProRS. To further delineate species-specific differences in acceptor stem recognition, atomic group mutagenesis was used to probe the role of sugar–phosphate backbone interactions in recognition of human tRNA^Pro. Incorporation of site-specific 2'-deoxynucleotides, as well as phosphorothioate and methylphosphonate modifications within the tRNA acceptor stem revealed an extensive network of interactions with specific functional groups proximal to the first base pair and the discriminator base. Backbone functional groups located at the base of the acceptor stem, especially the 2'-hydroxyl of A66, are also critical for aminoacylation catalytic efficiency by human ProRS. Therefore, in contrast to the bacterial system, backbone-specific interactions contribute significantly more to tRNA recognition by the human enzyme than base-specific interactions. Taken together with previous studies, these data show that ProRS-tRNA acceptor stem interactions have co-adapted through evolution from a mechanism involving ‘direct readout’ of nucleotide bases to one relying primarily on backbone-specific ‘indirect readout’.

Molecular dissection of Penelope transposable element regulatory machinery

2008-05-06

Penelope-like elements (PLEs) represent a new class of retroelements identified in more than 80 species belonging to at least 10 animal phyla. Penelope isolated from Drosophila virilis is the only known transpositionally active representative of this class. Although the size and structure of the Penelope major transcript has been previously described in both D. virilis and D. melanogaster transgenic strains, the architecture of the Penelope regulatory region remains unknown. In order to determine the localization of presumptive Penelope promoter and enhancer-like elements, segments of the putative Penelope regulatory region were linked to a CAT reporter gene and introduced into D. melanogaster by P-element-mediated transformation. The results obtained using ELISA to measure CAT expression levels and RNA studies, including RT–PCR, suggest that the active Penelope transposon contains an internal promoter similar to the TATA-less promoters of LINEs. The results also suggest that some of the Penelope regulatory sequences control the preferential expression in the ovaries of the adult flies by enhancing expression in the ovary and reducing expression in the carcass. The possible significance of the intron within Penelope for the function and evolution of PLEs, and the effect of Penelope insertions on adjacent genes, are discussed.

Bioinformatic and functional analysis of RNA secondary structure elements among different genera of human and animal caliciviruses

2008-05-06

The mechanism and role of RNA structure elements in the replication and translation of Caliciviridae remains poorly understood. Several algorithmically independent methods were used to predict secondary structures within the Norovirus, Sapovirus, Vesivirus and Lagovirus genera. All showed profound suppression of synonymous site variability (SSSV) at genomic 5' ends and the start of the sub-genomic (sg) transcript, consistent with evolutionary constraints from underlying RNA structure. A newly developed thermodynamic scanning method predicted RNA folding mapping precisely to regions of SSSV and at the genomic 3' end. These regions contained several evolutionarily conserved RNA secondary structures, of variable size and positions. However, all caliciviruses contained 3' terminal hairpins, and stem–loops in the anti-genomic strand invariably six bases upstream of the sg transcript, indicating putative roles as sg promoters. Using the murine norovirus (MNV) reverse-genetics system, disruption of 5' end stem–loops produced ~15- to 20-fold infectivity reductions, while disruption of the RNA structure in the sg promoter region and at the 3' end entirely destroyed replication ability. Restoration of infectivity by repair mutations in the sg promoter region confirmed a functional role for the RNA secondary structure, not the sequence. This study provides comprehensive bioinformatic resources for future functional studies of MNV and other caliciviruses.

A systematic characterization of factors that regulate Drosophila segmentation via a bacterial one-hybrid system

Noyes, M. B., Meng, X., Wakabayashi, A., Sinha, S., Brodsky, M. H., Wolfe, S. A. — 2008-05-06

Specificity data for groups of transcription factors (TFs) in a common regulatory network can be used to computationally identify the location of cis-regulatory modules in a genome. The primary limitation for this type of analysis is the paucity of specificity data that is available for the majority of TFs. We describe an omega-based bacterial one-hybrid system that provides a rapid method for characterizing DNA-binding specificities on a genome-wide scale. Using this system, 35 members of the Drosophila melanogaster segmentation network have been characterized, including representative members of all of the major classes of DNA-binding domains. A suite of web-based tools was created that uses this binding site dataset and phylogenetic comparisons to identify cis-regulatory modules throughout the fly genome. These tools allow specificities for any combination of factors to be used to perform rapid local or genome-wide searches for cis-regulatory modules. The utility of these factor specificities and tools is demonstrated on the well-characterized segmentation network. By incorporating specificity data on an additional 66 factors that we have characterized, our tools utilize ~14% of the predicted factors within the fly genome and provide an important new community resource for the identification of cis-regulatory modules.

Distinct roles of XRCC4 and Ku80 in non-homologous end-joining of endonuclease- and ionizing radiation-induced DNA double-strand breaks

2008-05-06

Non-homologous end-joining (NHEJ) of DNA double-strand breaks (DSBs) is mediated by two protein complexes comprising Ku80/Ku70/DNA-PKcs/Artemis and XRCC4/LigaseIV/XLF. Loss of Ku or XRCC4/LigaseIV function compromises the rejoining of radiation-induced DSBs and leads to defective V(D)J recombination. In this study, we sought to define how XRCC4 and Ku80 affect NHEJ of site-directed chromosomal DSBs in murine fibroblasts. We employed a recently developed reporter system based on the rejoining of I-SceI endonuclease-induced DSBs. We found that the frequency of NHEJ was reduced by more than 20-fold in XRCC4–/– compared to XRCC4+/+ cells, while a Ku80 knock-out reduced the rejoining efficiency by only 1.4-fold. In contrast, lack of either XRCC4 or Ku80 increased end degradation and shifted repair towards a mode that used longer terminal microhomologies for rejoining. However, both proteins proved to be essential for the repair of radiation-induced DSBs. The remarkably different phenotype of XRCC4- and Ku80-deficient cells with regard to the repair of enzyme-induced DSBs mirrors the embryonic lethality of XRCC4 knock-out mice as opposed to the viability of the Ku80 knock-out. Thus, I-SceI-induced breaks may resemble DSBs arising during normal DNA metabolism and mouse development. The removal of these breaks likely has different genetic requirements than the repair of radiation-induced DSBs.

The small RNA GlmY acts upstream of the sRNA GlmZ in the activation of glmS expression and is subject to regulation by polyadenylation in Escherichia coli

Reichenbach, B., Maes, A., Kalamorz, F., Hajnsdorf, E., Gorke, B. — 2008-05-06

In Escherichia coli the glmS gene encoding glucosamine 6-phosphate (GlcN-6-P) synthase GlmS is feedback regulated by GlcN-6-P in a pathway that involves the small RNA GlmZ. Expression of glmS is activated by the unprocessed form of GlmZ, which accumulates when the intracellular GlcN-6-P concentration decreases. GlmZ stabilizes a glmS transcript that derives from processing. Overexpression of a second sRNA, GlmY, also activates glmS expression in an unknown way. Furthermore, mutations in two genes, yhbJ and pcnB, cause accumulation of full-length GlmZ and thereby activate glmS expression. The function of yhbJ is unknown and pcnB encodes poly(A) polymerase PAP-I known to polyadenylate and destabilize RNAs. Here we show that GlmY acts indirectly in a way that depends on GlmZ. When the intracellular GlcN-6-P concentration decreases, GlmY accumulates and causes in turn accumulation of full-length GlmZ, which finally activates glmS expression. In glmZ mutants, GlmY has no effect on glmS, whereas artificially expressed GlmZ can activate glmS expression also in the absence of GlmY. Furthermore, we show that PAP-I acts at the top of this regulatory pathway by polyadenylating and destabilizing GlmY. In pcnB mutants, GlmY accumulates and induces glmS expression by stabilizing full-length GlmZ. Hence, the data reveal a regulatory cascade composed of two sRNAs, which responds to GlcN-6-P and is controlled by polyadenylation.

Real-time kinetics of restriction-modification gene expression after entry into a new host cell

Mruk, I., Blumenthal, R. M. — 2008-05-06

Most type II restriction–modification (R–M) systems produce separate restriction endonuclease (REase) and methyltransferase (MTase) proteins. After R–M system genes enter a new cell, protective MTase must appear before REase to avoid host chromosome cleavage. The basis for this apparent temporal regulation is not well understood. PvuII and some other R–M systems appear to achieve this delay by cotranscribing the REase gene with the gene for an autogenous transcription activator/repressor (the ‘C’ protein C.PvuII). To test this model, bacteriophage M13 was used to introduce the PvuII genes into a bacterial population in a relatively synchronous manner. REase mRNA and activity appeared ~10 min after those of the MTase, but never rose if there was an inactivating pvuIIC mutation. Infection with recombinant M13pvuII phage had little effect on cell growth, relative to infection with parental M13. However, infection of cells pre-expressing C.PvuII led to cessation of growth. This study presents the first direct demonstration of delayed REase expression, relative to MTase, when type II R–M genes enter a new host cell. Surprisingly, though the C and REase genes are cotranscribed, the pvuIIC portion of the mRNA was more abundant than the pvuIIR portion after stable establishment of the R–M system.

Mitogen-induced recruitment of ERK and MSK to SRE promoter complexes by ternary complex factor Elk-1

2008-05-06

Many eukaryotic genes are acutely regulated by extra-cellular signals. The c-fos serum response element (SRE) mediates transcriptional activation in response to mitogens through serum response factor (SRF)-dependent recruitment of Elk-1, a mitogen-activated protein kinase (MAPK)-responsive transcription factor. How subsequent events at SRE promoters stimulate initiation of transcription has yet to be fully resolved. Here we show that extra-cellular signal-regulated kinase (ERK) and mitogen and stress-activated kinase (MSK) are recruited to SRE promoter complexes in vitro and in vivo. Their recruitment in vitro correlates with Elk-1 binding and for ERK the D domain/KIM of Elk-1 is specifically involved. In vivo, recruitment of ERK and MSK is stimulated by mitogens, correlates with histone H3 phosphorylation and is impaired by Elk-1 knockdown. Immunocytochemistry and confocal microscopy reveal that ERK appears to associate to some extent with initiating rather than elongating RNA polymerase II. Taken together, our data add to the body of evidence implying that ERK and related MAPKs may fulfil a generic role at the promoters of acutely regulated genes.

Sequence homology and microhomology dominate chromosomal double-strand break repair in African trypanosomes

Glover, L., McCulloch, R., Horn, D. — 2008-05-06

Genetic diversity in fungi and mammals is generated through mitotic double-strand break-repair (DSBR), typically involving homologous recombination (HR) or non-homologous end joining (NHEJ). Microhomology-mediated joining appears to serve a subsidiary function. The African trypanosome, a divergent protozoan parasite, relies upon rearrangement of subtelomeric variant surface glycoprotein (VSG) genes to achieve antigenic variation. Evidence suggests an absence of NHEJ but chromosomal repair remains largely unexplored. We used a system based on I-SceI meganuclease and monitored temporally constrained DSBR at a specific chromosomal site in bloodstream form Trypanosoma brucei. In response to the lesion, adjacent single-stranded DNA was generated; the homologous strand-exchange factor, Rad51, accumulated into foci; a G₂M checkpoint was activated and >50% of cells displayed successful repair. Quantitative analysis of DSBR pathways employed indicated that inter-chromosomal HR dominated. HR displayed a strong preference for the allelic template but also the capacity to interact with homologous sequence on heterologous chromosomes. Intra-chromosomal joining was predominantly, and possibly exclusively, microhomology mediated, a situation unique among organisms examined to date. These DSBR pathways available to T. brucei likely underlie patterns of antigenic variation and the evolution of the vast VSG gene family.

A new kinetic model reveals the synergistic effect of E-, P- and A-sites on +1 ribosomal frameshifting

Liao, P.-Y., Gupta, P., Petrov, A. N., Dinman, J. D., Lee, K. H. — 2008-05-06

Programmed ribosomal frameshifting (PRF) is a process by which ribosomes produce two different polypeptides from the same mRNA. In this study, we propose three different kinetic models of +1 PRF, incorporating the effects of the ribosomal E-, P- and A-sites toward promoting efficient +1 frameshifting in Escherichia coli. Specifically, the timing of E-site tRNA dissociation is discussed within the context of the kinetic proofreading mechanism of aminoacylated tRNA (aa-tRNA) selection. Mathematical modeling using previously determined kinetic rate constants reveals that destabilization of deacylated tRNA in the E-site, rearrangement of peptidyl-tRNA in the P-site, and availability of cognate aa-tRNA corresponding to the A-site act synergistically to promote efficient +1 PRF. The effect of E-site codon:anticodon interactions on +1 PRF was also experimentally examined with a dual fluorescence reporter construct. The combination of predictive modeling and empirical testing allowed the rate constant for P-site tRNA slippage (k_s) to be estimated as k_s 1.9 s^–1 for the release factor 2 (RF2) frameshifting sequence. These analyses suggest that P-site tRNA slippage is the driving force for +1 ribosomal frameshifting while the presence of a ‘hungry codon’ in the A-site and destabilization in the E-site further enhance +1 PRF in E. coli.

A space-efficient and accurate method for mapping and aligning cDNA sequences onto genomic sequence

Gotoh, O. — 2008-05-06

The mapping and alignment of transcripts (cDNA, expressed sequence tag or amino acid sequences) onto a genomic sequence is a fundamental step for genome annotation, including gene finding and analyses of transcriptional activity, alternative splicing and nucleotide polymorphisms. As DNA sequence data of genomes and transcripts are accumulating at an unprecedented rate, steady improvement in accuracy, speed and space requirement in the computational tools for mapping/alignment is desired. We devised a multi-phase heuristic algorithm and implemented it in the development of the stand-alone computer program Spaln (space-efficient spliced alignment). Spaln is reasonably fast and space efficient; it requires <1 Gb of memory to map and align >120 000 Unigene sequences onto the unmasked whole human genome with a conventional computer, finishing the job in <6 h. With artificially introduced noise of various levels, Spaln significantly outperforms other leading alignment programs currently available with respect to the accuracy of mapped exon–intron structures. This performance is achieved without extensive learning procedures to adjust parameter values to a particular organism. According to the handiness and accuracy, Spaln may be used for studies on a wide area of genome analyses.

Determinants of a transcriptionally competent environment at the GM-CSF promoter

2008-05-06

Granulocyte macrophage-colony stimulating factor (GM-CSF) is produced by T cells, but not B cells, in response to immune signals. GM-CSF gene activation in response to T-cell stimulation requires remodelling of chromatin associated with the gene promoter, and these changes do not occur in B cells. While the CpG methylation status of the murine GM-CSF promoter shows no correlation with the ability of the gene to respond to activation, we find that the basal chromatin environment of the gene promoter influences its ability to respond to immune signals. In unstimulated T cells but not B cells, the GM-CSF promoter is selectively marked by enrichment of histone acetylation, and association of the chromatin-remodelling protein BRG1. BRG1 is removed from the promoter upon activation concomitant with histone depletion and BRG1 is required for efficient chromatin remodelling and transcription. Increasing histone acetylation at the promoter in T cells is paralleled by increased BRG1 recruitment, resulting in more rapid chromatin remodelling, and an associated increase in GM-CSF mRNA levels. Furthermore, increasing histone acetylation in B cells removes the block in chromatin remodelling and transcriptional activation of the GM-CSF gene. These data are consistent with a model in which histone hyperacetylation and BRG1 enrichment at the GM-CSF promoter, generate a chromatin environment competent to respond to immune signals resulting in gene activation.

The bacterial and mitochondrial ribosomal A-site molecular switches possess different conformational substates

Kondo, J., Westhof, E. — 2008-05-06

The A site of the small ribosomal subunit participates in the fidelity of decoding by switching between two states, a resting ‘off’ state and an active decoding ‘on’ state. Eight crystal structures of RNA duplexes containing two minimal decoding A sites of the Homo sapiens mitochondrial wild-type, the A1555G mutant or bacteria have been solved. The resting ‘off’ state of the mitochondrial wild-type A site is surprisingly different from that of the bacterial A site. The mitochondrial A1555G mutant has two types of the ‘off’ states; one is similar to the mitochondrial wild-type ‘off’ state and the other is similar to the bacterial ‘off’ state. Our present results indicate that the dynamics of the A site in bacteria and mitochondria are different, a property probably related to the small number of tRNAs used for decoding in mitochondria. Based on these structures, we propose a hypothesis for the molecular mechanism of non-syndromic hearing loss due to the mitochondrial A1555G mutation.

Autoregulation of the Escherichia coli melR promoter: repression involves four molecules of MelR

2008-05-06

The Escherichia coli MelR protein is a transcription activator that autoregulates its own promoter by repressing transcription initiation. Optimal repression requires MelR binding to a site that overlaps the melR transcription start point and to upstream sites. In this work, we have investigated the different determinants needed for optimal repression and their spatial requirements. We show that repression requires a complex involving four DNA-bound MelR molecules, and that the global CRP regulator plays little or no role.

Small ncRNA transcriptome analysis from Aspergillus fumigatus suggests a novel mechanism for regulation of protein synthesis

2008-05-06

Small non-protein-coding RNAs (ncRNAs) have systematically been studied in various model organisms from Escherichia coli to Homo sapiens. Here, we analyse the small ncRNA transcriptome from the pathogenic filamentous fungus Aspergillus fumigatus. To that aim, we experimentally screened for ncRNAs, expressed under various growth conditions or during specific developmental stages, by generating a specialized cDNA library from size-selected small RNA species. Our screen revealed 30 novel ncRNA candidates from known ncRNA classes such as small nuclear RNAs (snRNAs) and C/D box-type small nucleolar RNAs (C/D box snoRNAs). Additionally, several candidates for H/ACA box snoRNAs could be predicted by a bioinformatical screen. We also identified 15 candidates for ncRNAs, which could not be assigned to any known ncRNA class. Some of these ncRNA species are developmentally regulated implying a possible novel function in A. fumigatus development. Surprisingly, in addition to full-length tRNAs, we also identified 5'- or 3'-halves of tRNAs, only, which are likely generated by tRNA cleavage within the anti-codon loop. We show that conidiation induces tRNA cleavage resulting in tRNA depletion within conidia. Since conidia represent the resting state of A. fumigatus we propose that conidial tRNA depletion might be a novel mechanism to down-regulate protein synthesis in a filamentous fungus.

Transforming growth factor-{beta}-regulated miR-24 promotes skeletal muscle differentiation

2008-05-06

MicroRNAs (miRNAs) have recently been proposed as a versatile class of molecules involved in regulation of a variety of biological processes. However, the role of miRNAs in TGF-β-regulated biological processes is poorly addressed. In this study, we found that miR-24 was upregulated during myoblast differentiation and could be inhibited by TGF-β1. Using both a reporter assay and Northern blot analysis, we showed that TGF-β1 repressed miR-24 transcription which was dependent on the presence of Smad3 and a Smads binding site in the promoter region of miR-24. TGF-β1 was unable to inhibit miR-24 expression in Smad3-deficient myoblasts, which exhibited accelerated myogenesis. Knockdown of miR-24 led to reduced expression of myogenic differentiation markers in C2C12 cells, while ectopic expression of miR-24 enhanced differentiation, and partially rescued inhibited myogenesis by TGF-β1. This is the first study demonstrating a critical role for miRNAs in modulating TGF-β-dependent inhibition of myogenesis, and provides a novel mechanism of the genetic regulation of TGF-β signaling during skeletal muscle differentiation.

The relationship of potential G-quadruplex sequences in cis-upstream regions of the human genome to SP1-binding elements

Todd, A. K., Neidle, S. — 2008-05-06

We have carried out a survey of potential quadruplex structure sequences (PQSS), which occur in the immediate upstream region (500 bp) of human genes. By examining the number and distribution of these we have established that there is a clear link between them and the occurrence of the SP1-binding element ‘GGGCGG’, such that a large number of upstream PQSS incorporate the SP1-binding element.

Dissecting protein-RNA recognition sites

Bahadur, R. P., Zacharias, M., Janin, J. — 2008-05-06

We analyze the protein–RNA interfaces in 81 transient binary complexes taken from the Protein Data Bank. Those with tRNA or duplex RNA are larger than with single-stranded RNA, and comparable in size to protein–DNA interfaces. The protein side bears a strong positive electrostatic potential and resembles protein–DNA interfaces in its amino acid composition. On the RNA side, the phosphate contributes less, and the sugar much more, to the interaction than in protein–DNA complexes. On average, protein–RNA interfaces contain 20 hydrogen bonds, 7 that involve the phosphates, 5 the sugar 2'OH, and 6 the bases, and 32 water molecules. The average H-bond density per unit buried surface area is less with tRNA or single-stranded RNA than with duplex RNA. The atomic packing is also less compact in interfaces with tRNA. On the protein side, the main chain NH and Arg/Lys side chains account for nearly half of all H-bonds to RNA; the main chain CO and side chain acceptor groups, for a quarter. The 2'OH is a major player in protein–RNA recognition, and shape complementarity an important determinant, whereas electrostatics and direct base–protein interactions play a lesser part than in protein–DNA recognition.

Human abasic endonuclease action on multilesion abasic clusters: implications for radiation-induced biological damage

Paap, B., Wilson, D. M., Sutherland, B. M. — 2008-05-06

Clustered damages—two or more closely opposed abasic sites, oxidized bases or strand breaks—are induced in DNA by ionizing radiation and by some radiomimetic drugs. They are potentially mutagenic or lethal. High complexity, multilesion clusters (three or more lesions) are hypothesized as repair-resistant and responsible for the greater biological damage induced by high linear energy transfer radiation (e.g. charged particles) than by low linear energy transfer X- or -rays. We tested this hypothesis by assessing human abasic endonuclease Ape1 activity on two- and multiple-lesion abasic clusters. We constructed cluster-containing oligonucleotides using a central variable cassette with abasic site(s) at specific locations, and 5' and 3' terminal segments tagged with visually distinctive fluorophores. The results indicate that in two- or multiple-lesion clusters, the spatial arrangement of uni-sided positive [in which the opposing strand lesion(s) is 3' to the base opposite the reference lesion)] or negative polarity [opposing strand lesion(s) 5' to the base opposite the reference lesion] abasic clusters is key in determining Ape1 cleavage efficiency. However, no bipolar clusters (minimally three-lesions) were good Ape1 substrates. The data suggest an underlying molecular mechanism for the higher levels of biological damage associated with agents producing complex clusters: the induction of highly repair-resistant bipolar clusters.

Recognition of tRNALeu by Aquifex aeolicus leucyl-tRNA synthetase during the aminoacylation and editing steps

Yao, P., Zhu, B., Jaeger, S., Eriani, G., Wang, E.-D. — 2008-05-06

Recognition of tRNA by the cognate aminoacyl-tRNA synthetase during translation is crucial to ensure the correct expression of the genetic code. To understand tRNA^Leu recognition sets and their evolution, the recognition of tRNA^Leu by the leucyl-tRNA synthetase (LeuRS) from the primitive hyperthermophilic bacterium Aquifex aeolicus was studied by RNA probing and mutagenesis. The results show that the base A73; the core structure of tRNA formed by the tertiary interactions U8–A14, G18–U55 and G19–C56; and the orientation of the variable arm are critical elements for tRNA^Leu aminoacylation. Although dispensable for aminoacylation, the anticodon arm carries discrete editing determinants that are required for stabilizing the conformation of the post-transfer editing state and for promoting translocation of the tRNA acceptor arm from the synthetic to the editing site.

Structure of the DNA-binding domain of NgTRF1 reveals unique features of plant telomere-binding proteins

2008-05-06

Telomeres are protein–DNA elements that are located at the ends of linear eukaryotic chromosomes. In concert with various telomere-binding proteins, they play an essential role in genome stability. We determined the structure of the DNA-binding domain of NgTRF1, a double-stranded telomere-binding protein of tobacco, using multidimensional NMR spectroscopy and X-ray crystallography. The DNA-binding domain of NgTRF1 contained the Myb-like domain and C-terminal Myb-extension that is characteristic of plant double-stranded telomere-binding proteins. It encompassed amino acids 561–681 (NgTRF1^561–681), and was composed of 4 -helices. We also determined the structure of NgTRF1^561–681 bound to plant telomeric DNA. We identified several amino acid residues that interacted directly with DNA, and confirmed their role in the binding of NgTRF1 to telomere using site-directed mutagenesis. Based on a structural comparison of the DNA-binding domains of NgTRF1 and human TRF1 (hTRF1), NgTRF1 has both common and unique DNA-binding properties. Interaction of Myb-like domain with telomeric sequences is almost identical in NgTRF1^561–681 with the DNA-binding domain of hTRF1. The interaction of Arg-638 with the telomeric DNA, which is unique in NgTRF1^561–681, may provide the structural explanation for the specificity of NgTRF1 to the plant telomere sequences, (TTTAGGG)_n.

Spontaneous symmetry breaking in genome evolution

Ryabov, Y., Gribskov, M. — 2008-05-06

The quest for evolutionary mechanisms providing separation between the coding (exons) and noncoding (introns) parts of genomic DNA remains an important focus of genetics. This work combines an analysis of the most recent achievements of genomics and fundamental concepts of random processes to provide a novel point of view on genome evolution. Exon sizes in sequenced genomes show a lognormal distribution typical of a random Kolmogoroff fractioning process. This implies that the process of intron incretion may be independent of exon size, and therefore could be dependent on intron–exon boundaries. All genomes examined have two distinctive classes of exons, each with different evolutionary histories. In the framework proposed in this article, these two classes of exons can be derived from a hypothetical ancestral genome by (spontaneous) symmetry breaking. We note that one of these exon classes comprises mostly alternatively spliced exons.

Intracellular delivery of an anionic antisense oligonucleotide via receptor-mediated endocytosis

2008-05-06

We describe the synthesis and characterization of a 5' conjugate between a 2'-O-Me phosphorothioate antisense oligonucleotide and a bivalent RGD (arginine–glycine–aspartic acid) peptide that is a high-affinity ligand for the vβ3 integrin. We used vβ3-positive melanoma cells transfected with a reporter comprised of the firefly luciferase gene interrupted by an abnormally spliced intron. Intranuclear delivery of a specific antisense oligonucleotide (termed 623) corrects splicing and allows luciferase expression in these cells. The RGD–623 conjugate or a cationic lipid-623 complex produced significant increases in luciferase expression, while ‘free’ 623 did not. However, the kinetics of luciferase expression was distinct; the RGD–623 conjugate produced a gradual increase followed by a gradual decline, while the cationic lipid-623 complex caused a rapid increase followed by a monotonic decline. The subcellular distribution of the oligonucleotide delivered using cationic lipids included both cytoplasmic vesicles and the nucleus, while the RGD–623 conjugate was primarily found in cytoplasmic vesicles that partially co-localized with a marker for caveolae. Both the cellular uptake and the biological effect of the RGD–623 conjugate were blocked by excess RGD peptide. These observations suggest that the bivalent RGD peptide–oligonucleotide conjugate enters cells via a process of receptor-mediated endocytosis mediated by the vβ3 integrin.

The biological function of some human transcription factor binding motifs varies with position relative to the transcription start site

2008-05-06

A number of previous studies have predicted transcription factor binding sites (TFBSs) by exploiting the position of genomic landmarks like the transcriptional start site (TSS). The studies’ methods are generally too computationally intensive for genome-scale investigation, so the full potential of ‘positional regulomics’ to discover TFBSs and determine their function remains unknown. Because databases often annotate the genomic landmarks in DNA sequences, the methodical exploitation of positional regulomics has become increasingly urgent. Accordingly, we examined a set of 7914 human putative promoter regions (PPRs) with a known TSS. Our methods identified 1226 eight-letter DNA words with significant positional preferences with respect to the TSS, of which only 608 of the 1226 words matched known TFBSs. Many groups of genes whose PPRs contained a common word displayed similar expression profiles and related biological functions, however. Most interestingly, our results included 78 words, each of which clustered significantly in two or three different positions relative to the TSS. Often, the gene groups corresponding to different positional clusters of the same word corresponded to diverse functions, e.g. activation or repression in different tissues. Thus, different clusters of the same word likely reflect the phenomenon of ‘positional regulation’, i.e. a word's regulatory function can vary with its position relative to a genomic landmark, a conclusion inaccessible to methods based purely on sequence. Further integrative analysis of words co-occurring in PPRs also yielded 24 different groups of genes, likely identifying cis-regulatory modules de novo. Whereas comparative genomics requires precise sequence alignments, positional regulomics exploits genomic landmarks to provide a ‘poor man's alignment’. By exploiting the phenomenon of positional regulation, it uses position to differentiate the biological functions of subsets of TFBSs sharing a common sequence motif.

A BBP-Mud2p heterodimer mediates branchpoint recognition and influences splicing substrate abundance in budding yeast

Wang, Q., Zhang, L., Lynn, B., Rymond, B. C. — 2008-05-06

The 3' end of mammalian introns is marked by the branchpoint binding protein, SF1, and the U2AF65-U2AF35 heterodimer bound at an adjacent sequence. Baker's yeast has equivalent proteins, branchpoint binding protein (BBP) (SF1) and Mud2p (U2AF65), but lacks an obvious U2AF35 homolog, leaving open the question of whether another protein substitutes during spliceosome assembly. Gel filtration, affinity selection and mass spectrometry were used to show that rather than a U2AF65/U2AF35-like heterodimer, Mud2p forms a complex with BBP without a third (U2AF35-like) factor. Using mutants of MUD2 and BBP, we show that the BBP–Mud2p complex bridges partner-specific Prp39p, Mer1p, Clf1p and Smy2p two-hybrid interactions. In addition to inhibiting Mud2p association, the bbp56 mutation impairs splicing, enhances pre-mRNA release from the nucleus, and similar to a mud2::KAN knockout, suppresses a lethal sub2::KAN mutation. Unexpectedly, rather than exacerbating bbp56, the mud2::KAN mutation partially suppresses a pre-mRNA accumulation defect observed with bbp56. We propose that a BBP–Mud2p heterodimer binds as a unit to the branchpoint in vivo and serves as a target for the Sub2p-DExD/H-box ATPase and for other splicing factors during spliceosome assembly. In addition, our results suggest the possibility that the Mud2p may enhance the turnover of pre-mRNA with impaired BBP-branchpoint association.

Structural probing of the HIV-1 polypurine tract RNA:DNA hybrid using classic nucleic acid ligands

2008-05-06

The interactions of archetypical nucleic acid ligands with the HIV-1 polypurine tract (PPT) RNA:DNA hybrid, as well as analogous DNA:DNA, RNA:RNA and swapped hybrid substrates, were used to probe structural features of the PPT that contribute to its specific recognition and processing by reverse transcriptase (RT). Results from intercalative and groove-binding ligands indicate that the wild-type PPT hybrid does not contain any strikingly unique groove geometries and/or stacking arrangements that might contribute to the specificity of its interaction with RT. In contrast, neomycin bound preferentially and selectively to the PPT near the 5'(rA)₄:(dT)₄ tract and the 3' PPT-U3 junction. Nuclear magnetic resonance data from a complex between HIV-1 RT and the PPT indicate RT contacts within the same regions highlighted on the PPT by neomycin. These observations, together with the fact that the sites are correctly spaced to allow interaction with residues in the ribonuclease H (RNase H) active site and thumb subdomain of the p66 RT subunit, suggest that despite the long cleft employed by RT to make contact with nucleic acids substrates, these sites provide discrete binding units working in concert to determine not only specific PPT recognition, but also its orientation on the hybrid structure.

Elevated {alpha}-synuclein mRNA levels in individual UV-laser-microdissected dopaminergic substantia nigra neurons in idiopathic Parkinson's disease

Grundemann, J., Schlaudraff, F., Haeckel, O., Liss, B. — 2008-04-17

The presynaptic protein -synuclein is involved in several neurodegenerative diseases, including Parkinson's disease (PD). In rare familial forms of PD, causal mutations (PARK1) as well as multiplications (PARK4) of the -synuclein gene have been identified. In sporadic, idiopathic PD, abnormal accumulation and deposition of -synuclein might also cause degeneration of dopaminergic midbrain neurons, the clinically most relevant neuronal population in PD. Thus, cell-specific quantification of -synuclein expression-levels in dopaminergic neurons from idiopathic PD patients in comparison to controls would provide essential information about contributions of -synuclein to the etiology of PD. However, a number of previous studies addressing this question at the tissue-level yielded varying results regarding -synuclein expression. To increase specificity, we developed a cell-specific approach for mRNA quantification that also took into account the important issue of variable RNA integrities of the individual human postmortem brain samples. We demonstrate that PCR –amplicon size can confound quantitative gene-expression analysis, in particular of partly degraded RNA. By combining optimized UV-laser microdissection- and quantitative RT–PCR-techniques with suitable PCR assays, we detected significantly elevated -synuclein mRNA levels in individual, surviving neuromelanin- and tyrosine hydroxylase-positive substantia nigra dopaminergic neurons from idiopathic PD brains compared to controls. These results strengthen the pathophysiologic role of transcriptional dysregulation of the -synuclein gene in sporadic PD.

High-resolution array comparative genomic hybridization of single micrometastatic tumor cells

2008-04-17

Only few selected cancer cells drive tumor progression and are responsible for therapy resistance. Their specific genomic characteristics, however, are largely unknown because high-resolution genome analysis is currently limited to DNA pooled from many cells. Here, we describe a protocol for array comparative genomic hybridization (array CGH), which enables the detection of DNA copy number changes in single cells. Combining a PCR-based whole genome amplification method with arrays of highly purified BAC clones we could accurately determine known chromosomal changes such as trisomy 21 in single leukocytes as well as complex genomic imbalances of single cell line cells. In single T47D cells aberrant regions as small as 1–2 Mb were identified in most cases when compared to non-amplified DNA from 10⁶ cells. Most importantly, in single micrometastatic cancer cells isolated from bone marrow of breast cancer patients, we retrieved and confirmed amplifications as small as 4.4 and 5 Mb. Thus, high-resolution genome analysis of single metastatic precursor cells is now possible and may be used for the identification of novel therapy target genes.

Labeling of unique sequences in double-stranded DNA at sites of vicinal nicks generated by nicking endonucleases

Kuhn, H., Frank-Kamenetskii, M. D. — 2008-04-17

We describe a new approach for labeling of unique sequences within dsDNA under nondenaturing conditions. The method is based on the site-specific formation of vicinal nicks, which are created by nicking endonucleases (NEases) at specified DNA sites on the same strand within dsDNA. The oligomeric segment flanked by both nicks is then substituted, in a strand displacement reaction, by an oligonucleotide probe that becomes covalently attached to the target site upon subsequent ligation. Monitoring probe hybridization and ligation reactions by electrophoretic mobility retardation assay, we show that selected target sites can be quantitatively labeled with excellent sequence specificity. In these experiments, predominantly probes carrying a target-independent 3' terminal sequence were employed. At target labeling, thus a branched DNA structure known as 3'-flap DNA is obtained. The single-stranded terminus in 3'-flap DNA is then utilized to prime the replication of an externally supplied ssDNA circle in a rolling circle amplification (RCA) reaction. In model experiments with samples comprised of genomic -DNA and human herpes virus 6 type B (HHV-6B) DNA, we have used our labeling method in combination with surface RCA as reporter system to achieve both high sequence specificity of dsDNA targeting and high sensitivity of detection. The method can find applications in sensitive and specific detection of viral duplex DNA.

wuHMM: a robust algorithm to detect DNA copy number variation using long oligonucleotide microarray data

2008-04-17

Copy number variants (CNVs) are currently defined as genomic sequences that are polymorphic in copy number and range in length from 1000 to several million base pairs. Among current array-based CNV detection platforms, long-oligonucleotide arrays promise the highest resolution. However, the performance of currently available analytical tools suffers when applied to these data because of the lower signal:noise ratio inherent in oligonucleotide-based hybridization assays. We have developed wuHMM, an algorithm for mapping CNVs from array comparative genomic hybridization (aCGH) platforms comprised of 385 000 to more than 3 million probes. wuHMM is unique in that it can utilize sequence divergence information to reduce the false positive rate (FPR). We apply wuHMM to 385K-aCGH, 2.1M-aCGH and 3.1M-aCGH experiments comparing the 129X1/SvJ and C57BL/6J inbred mouse genomes. We assess wuHMM's performance on the 385K platform by comparison to the higher resolution platforms and we independently validate 10 CNVs. The method requires no training data and is robust with respect to changes in algorithm parameters. At a FPR of <10%, the algorithm can detect CNVs with five probes on the 385K platform and three on the 2.1M and 3.1M platforms, resulting in effective resolutions of 24 kb, 2–5 kb and 1 kb, respectively.

Sensitive Melting Analysis after Real Time- Methylation Specific PCR (SMART-MSP): high-throughput and probe-free quantitative DNA methylation detection

Kristensen, L. S., Mikeska, T., Krypuy, M., Dobrovic, A. — 2008-04-17

DNA methylation changes that are recurrent in cancer have generated great interest as potential biomarkers for the early detection and monitoring of cancer. In such situations, essential information is missed if the methylation detection is purely qualitative. We describe a new probe-free quantitative methylation-specific PCR (MSP) assay that incorporates evaluation of the amplicon by high-resolution melting (HRM) analysis. Depending on amplicon design, different types of information can be obtained from the HRM analysis. Much of this information cannot be obtained by electrophoretic analysis. In particular, identification of false positives due to incomplete bisulphite conversion or false priming is possible. Heterogeneous methylation can also be distinguished from homogeneous methylation. As proof of principle, we have developed assays for the promoter regions of the CDH1, DAPK1, CDKN2A (p16^INK4a) and RARB genes. We show that highly accurate quantification is possible in the range from 100% to 0.1% methylated template when 25 ng of bisulphite-modified DNA is used as a template for PCR. We have named this new approach to quantitative methylation detection, Sensitive Melting Analysis after Real Time (SMART)-MSP.

Positional gene enrichment analysis of gene sets for high-resolution identification of overrepresented chromosomal regions

De Preter, K., Barriot, R., Speleman, F., Vandesompele, J., Moreau, Y. — 2008-04-17

The search for feature enrichment is a widely used method to characterize a set of genes. While several tools have been designed for nominal features such as Gene Ontology annotations or KEGG Pathways, very little has been proposed to tackle numerical features such as the chromosomal positions of genes. For instance, microarray studies typically generate gene lists that are differentially expressed in the sample subgroups under investigation, and when studying diseases caused by genome alterations, it is of great interest to delineate the chromosomal regions that are significantly enriched in these lists. In this article, we present a positional gene enrichment analysis method (PGE) for the identification of chromosomal regions that are significantly enriched in a given set of genes. The strength of our method relies on an original query optimization approach that allows to virtually consider all the possible chromosomal regions for enrichment, and on the multiple testing correction which discriminates truly enriched regions versus those that can occur by chance. We have developed a Web tool implementing this method applied to the human genome (http://www.esat.kuleuven.be/~bioiuser/pge). We validated PGE on published lists of differentially expressed genes. These analyses showed significant overrepresentation of known aberrant chromosomal regions.

Generation of {rho}0 cells utilizing a mitochondrially targeted restriction endonuclease and comparative analyses

2008-04-17

Eukaryotic cells devoid of mitochondrial DNA (⁰ cells) were originally generated under artificial growth conditions utilizing ethidium bromide. The chemical is known to intercalate preferentially with the mitochondrial double-stranded DNA thereby interfering with enzymes of the replication machinery. ⁰ cell lines are highly valuable tools to study human mitochondrial disorders because they can be utilized in cytoplasmic transfer experiments. However, mutagenic effects of ethidium bromide onto the nuclear DNA cannot be excluded. To foreclose this mutagenic character during the development of ⁰ cell lines, we developed an extremely mild, reliable and timesaving method to generate ⁰ cell lines within 3–5 days based on an enzymatic approach. Utilizing the genes for the restriction endonuclease EcoRI and the fluorescent protein EGFP that were fused to a mitochondrial targeting sequence, we developed a CMV-driven expression vector that allowed the temporal expression of the resulting fusion enzyme in eukaryotic cells. Applied on the human cell line 143B.TK^– the active protein localized to mitochondria and induced the complete destruction of endogenous mtDNA. Mouse and rat ⁰ cell lines were also successfully created with this approach. Furthermore, the newly established 143B.TK^– ⁰ cell line was characterized in great detail thereby releasing interesting insights into the morphology and ultra structure of human ⁰ mitochondria.

PDZ domain-mediated dimerization and homeodomain-directed specificity are required for high-affinity DNA binding by SATB1

2008-04-17

To better understand DNA recognition and transcription activity by SATB1, the T-lineage-enriched chromatin organizer and transcription factor, we have determined its optimal DNA-binding sequence by random oligonucleotide selection. The consensus SATB1-binding sequence (CSBS) comprises a palindromic sequence in which two identical AT-rich half-sites are arranged as inverted repeats flanking a central cytosine or guanine. Strikingly, the CSBS half-site is identical to the conserved element ‘TAATA’ bound by the known homeodomains (HDs). Furthermore, we show that the high-affinity binding of SATB1 to DNA is dimerization-dependent and the HD also binds in similar fashion. Binding studies using HD-lacking SATB1 and binding target with increased spacer between the two half-sites led us to propose a model for SATB1–DNA complex in which the HDs bind in an antiparallel fashion to the palindromic consensus element via minor groove, bridged by the PDZ-like dimerization domain. CSBS-driven in vivo reporter analysis indicated that SATB1 acts as a repressor upon binding to the CSBS and most of its derivatives and the extent of repression is proportional to SATB1's binding affinity to these sequences. These studies provide mechanistic insights into the mode of DNA binding and its effect on the regulation of transcription by SATB1.

Lsr2 of Mycobacterium tuberculosis is a DNA-bridging protein

2008-04-17

Lsr2 is a small, basic protein present in Mycobacterium and related actinomycetes. Recent studies suggest that Lsr2 is a regulatory protein involved in multiple cellular processes including cell wall biosynthesis and antibiotic resistance. However, the underlying molecular mechanisms remain unknown. In this article, we performed biochemical studies of Lsr2–DNA interactions and structure–function analysis of Lsr2. Analysis by atomic force microscopy revealed that Lsr2 has the ability to bridge distant DNA segments, suggesting that Lsr2 plays a role in the overall organization and compactness of the nucleoid. Mutational analysis identified critical residues and selection of dominant negative mutants demonstrated that both DNA binding and protein oligomerization are essential for the normal functions of Lsr2 in vivo. These results provide strong evidence that Lsr2 is a DNA bridging protein, which represents the first identification of such proteins in bacteria phylogenetically distant from the Enterobacteriaceae. DNA bridging by Lsr2 also provides a mechanism of transcriptional regulation by Lsr2.

Functional dissection of siRNA sequence by systematic DNA substitution: modified siRNA with a DNA seed arm is a powerful tool for mammalian gene silencing with significantly reduced off-target effect

2008-04-17

Short interfering RNA (siRNA)-based RNA interference (RNAi) is widely used for target gene knockdown in mammalian cells. To clarify the position-dependent functions of ribonucleotides in siRNA, siRNAs with various DNA substitutions were constructed. The following could be simultaneously replaced with DNA without substantial loss of gene-silencing activity: the seed arm, which occupies positions 2–8 from the 5'end of the guide strand; its complementary sequence; the 5'end of the guide strand and the 3'overhang of the passenger strand. However, most part of the 3' two-thirds of the guide strand could not be replaced with DNA, possibly due to binding of RNA-recognition proteins such as TRBP2 and Ago2. The passenger strand with DNA in the 3'end proximal region was incapable of inducing off-target effect. Owing to lesser stability of DNA–RNA hybrid than RNA duplex, modified siRNAs with DNA substitution in the seed region were, in most cases, incapable to exert unintended gene silencing due to seed sequence homology. Thus, it may be possible to design DNA–RNA chimeras which effectively silence mammalian target genes without silencing unintended genes.

Eukaryotic Y-family polymerases bypass a 3-methyl-2'-deoxyadenosine analog in vitro and methyl methanesulfonate-induced DNA damage in vivo

2008-04-17

N3-methyl-adenine (3MeA) is the major cytotoxic lesion formed in DNA by S_N2 methylating agents. The lesion presumably blocks progression of cellular replicases because the N3-methyl group hinders interactions between the polymerase and the minor groove of DNA. However, this hypothesis has yet to be rigorously proven, as 3MeA is intrinsically unstable and is converted to an abasic site, which itself is a blocking lesion. To circumvent these problems, we have chemically synthesized a 3-deaza analog of 3MeA (3dMeA) as a stable phosphoramidite and have incorporated the analog into synthetic oligonucleotides that have been used in vitro as templates for DNA replication. As expected, the 3dMeA lesion blocked both human DNA polymerases and . In contrast, human polymerases , and , as well as Saccharomyces cerevisiae pol were able to bypass the lesion, albeit with varying efficiencies and accuracy. To confirm the physiological relevance of our findings, we show that in S. cerevisiae lacking Mag1-dependent 3MeA repair, pol (Rad30) contributes to the survival of cells exposed to methyl methanesulfonate (MMS) and in the absence of Mag1, Rad30 and Rev3, human polymerases , and are capable of restoring MMS-resistance to the normally MMS-sensitive strain.

Computer design of obligate heterodimer meganucleases allows efficient cutting of custom DNA sequences

Fajardo-Sanchez, E., Stricher, F., Paques, F., Isalan, M., Serrano, L. — 2008-04-17

Meganucleases cut long (>12 bp) unique sequences in genomes and can be used to induce targeted genome engineering by homologous recombination in the vicinity of their cleavage site. However, the use of natural meganucleases is limited by the repertoire of their target sequences, and considerable efforts have been made to engineer redesigned meganucleases cleaving chosen targets. Homodimeric meganucleases such as I-CreI have provided a scaffold, but can only be modified to recognize new quasi-palindromic DNA sequences, limiting their general applicability. Other groups have used dimer-interface redesign and peptide linkage to control heterodimerization between related meganucleases such as I-DmoI and I-CreI, but until now there has been no application of this aimed specifically at the scaffolds from existing combinatorial libraries of I-CreI. Here, we show that engineering meganucleases to form obligate heterodimers results in functional endonucleases that cut non-palindromic sequences. The protein design algorithm (FoldX v2.7) was used to design specific heterodimer interfaces between two meganuclease monomers, which were themselves engineered to recognize different DNA sequences. The new monomers favour functional heterodimer formation and prevent homodimer site recognition. This design massively increases the potential repertoire of DNA sequences that can be specifically targeted by designed I-CreI meganucleases and opens the way to safer targeted genome engineering.

Trace amounts of 8-oxo-dGTP in mitochondrial dNTP pools reduce DNA polymerase {gamma} replication fidelity

Pursell, Z. F., McDonald, J. T., Mathews, C. K., Kunkel, T. A. — 2008-04-17

Replication of the mitochondrial genome by DNA polymerase requires dNTP precursors that are subject to oxidation by reactive oxygen species generated by the mitochondrial respiratory chain. One such oxidation product is 8-oxo-dGTP, which can compete with dTTP for incorporation opposite template adenine to yield A-T to C-G transversions. Recent reports indicate that the ratio of undamaged dGTP to dTTP in mitochondrial dNTP pools from rodent tissues varies from ~1:1 to >100:1. Within this wide range, we report here the proportion of 8-oxo-dGTP in the dNTP pool that would be needed to reduce the replication fidelity of human DNA polymerase . When various in vivo mitochondrial dNTP pools reported previously were used here in reactions performed in vitro, 8-oxo-dGTP was readily incorporated opposite template A and the resulting 8-oxo-G-A mismatch was not proofread efficiently by the intrinsic 3' exonuclease activity of pol . At the dNTP ratios reported in rodent tissues, whether highly imbalanced or relatively balanced, the amount of 8-oxo-dGTP needed to reduce fidelity was <1% of dGTP. Moreover, direct measurements reveal that 8-oxo-dGTP is present at such concentrations in the mitochondrial dNTP pools of several rat tissues. The results suggest that oxidized dNTP precursors may contribute to mitochondrial mutagenesis in vivo, which could contribute to mitochondrial dysfunction and disease.

Saccharomyces cerevisiae Mus81-Mms4 is a catalytic, DNA structure-selective endonuclease

Ehmsen, K. T., Heyer, W.-D. — 2008-04-17

The DNA structure-selective endonuclease Mus81-Mms4/Eme1 is a context-specific recombination factor that supports DNA replication, but is not essential for DSB repair in Saccharomyces cerevisiae. We overexpressed Mus81-Mms4 in S. cerevisiae, purified the heterodimer to apparent homogeneity, and performed a classical enzymological characterization. Kinetic analysis (k_cat, K_M) demonstrated that Mus81-Mms4 is catalytically active and identified three substrate classes in vitro. Class I substrates reflect low K_M (3–7 nM) and high k_cat (~1 min^–1) and include the nicked Holliday junction, 3'-flapped and replication fork-like structures. Class II substrates share low K_M (1–6 nM) but low k_cat (≤0.3 min^–1) relative to Class I substrates and include the D-loop and partial Holliday junction. The splayed Y junction defines a class III substrate having high K_M (~30 nM) and low k_cat (0.26 min^–1). Holliday junctions assembled from oligonucleotides with or without a branch migratable core were negligibly cut in vitro. We found that Mus81 and Mms4 are phosphorylated constitutively and in the presence of the genotoxin MMS. The endogenous complex purified in either modification state is negligibly active on Holliday junctions. Hence, Holliday junction incision activity in vitro cannot be attributed to the Mus81-Mms4 heterodimer in isolation.

Epigenetics of a tandem DNA repeat: chromatin DNaseI sensitivity and opposite methylation changes in cancers

2008-04-17

DNA methylation and chromatin DNaseI sensitivity were analyzed in and adjacent to D4Z4 repeat arrays, which consist of 1 to ~100 tandem 3.3-kb units at subtelomeric 4q and 10q. D4Z4 displayed hypomethylation in some cancers and hypermethylation in others relative to normal tissues. Surprisingly, in cancers with extensive D4Z4 methylation there was a barrier to hypermethylation spreading to the beginning of this disease-associated array (facioscapulohumeral muscular dystrophy, FSHD) despite sequence conservation in repeat units throughout the array. We infer a different chromatin structure at the proximal end of the array than at interior repeats, consistent with results from chromatin DNaseI sensitivity assays indicating a boundary element near the beginning of the array. The relative chromatin DNaseI sensitivity in FSHD and control myoblasts and lymphoblasts was as follows: a non-genic D4Z4-adjacent sequence (p13E-11, array-proximal)> untranscribed gene standards > D4Z4 arrays> constitutive heterochromatin (satellite 2; P < 10^–4 for all comparisons). Cancers displaying D4Z4 hypermethylation also exhibited a hypermethylation-resistant subregion within the 3.3-kb D4Z4 repeat units. This subregion contains runs of G that form G-quadruplexes in vitro. Unusual DNA structures might contribute to topological constraints that link short 4q D4Z4 arrays to FSHD and make long ones phenotypically neutral.

Nucleosome formation with the testis-specific histone H3 variant, H3t, by human nucleosome assembly proteins in vitro

Tachiwana, H., Osakabe, A., Kimura, H., Kurumizaka, H. — 2008-04-17

Five non-allelic histone H3 variants, H3.1, H3.2, H3.3, H3t and CENP-A, have been identified in mammals. H3t is robustly expressed in the testis, and thus was assigned as the testis-specific H3 variant. However, recent proteomics and tissue-specific RT-PCR experiments revealed a small amount of H3t expression in somatic cells. In the present study, we purified human H3t as a recombinant protein, and showed that H3t/H4 forms nucleosomes with H2A/H2B by the salt-dialysis method, like the conventional H3.1/H4. We found that H3t/H4 is not efficiently incorporated into the nucleosome by human Nap1 (hNap1), due to its defective H3t/H4 deposition on DNA. In contrast, human Nap2 (hNap2), a paralog of hNap1, promotes nucleosome assembly with H3t/H4. Mutational analyses revealed that the Ala111 residue, which is conserved among H3.1, H3.2 and H3.3, but not in H3t, is the essential residue for the hNap1-mediated nucleosome assembly. These results suggest that H3t may be incorporated into chromatin by a specific chaperone-mediated pathway.

LARP7 is a stable component of the 7SK snRNP while P-TEFb, HEXIM1 and hnRNP A1 are reversibly associated

2008-04-17

Regulation of the elongation phase of RNA polymerase II transcription by P-TEFb is a critical control point for gene expression. The activity of P-TEFb is regulated, in part, by reversible association with one of two HEXIMs and the 7SK snRNP. A recent proteomics survey revealed that P-TEFb and the HEXIMs are tightly connected to two previously-uncharacterized proteins, the methyphosphate capping enzyme, MEPCE, and a La-related protein, LARP7. Glycerol gradient sedimentation analysis of lysates from cells treated with P-TEFb inhibitors, suggested that the 7SK snRNP reorganized such that LARP7 and 7SK remained associated after P-TEFb and HEXIM1 were released. Immunodepletion of LARP7 also depleted most of the 7SK regardless of the presence of P-TEFb, HEXIM or hnRNP A1 in the complex. Small interfering RNA knockdown of LARP7 in human cells decreased the steady-state level of 7SK, led to an initial increase in free P-TEFb and increased Tat transactivation of the HIV-1 LTR. Knockdown of LARP7 or 7SK ultimately caused a decrease in total P-TEFb protein levels. Our studies have identified LARP7 as a 7SK-binding protein and suggest that free P-TEFb levels are determined by a balance between release from the large form and reduction of total P-TEFb.

Phylogenetic classification of short environmental DNA fragments

2008-04-17

Metagenomics is providing striking insights into the ecology of microbial communities. The recently developed massively parallel 454 pyrosequencing technique gives the opportunity to rapidly obtain metagenomic sequences at a low cost and without cloning bias. However, the phylogenetic analysis of the short reads produced represents a significant computational challenge. The phylogenetic algorithm CARMA for predicting the source organisms of environmental 454 reads is described. The algorithm searches for conserved Pfam domain and protein families in the unassembled reads of a sample. These gene fragments (environmental gene tags, EGTs), are classified into a higher-order taxonomy based on the reconstruction of a phylogenetic tree of each matching Pfam family. The method exhibits high accuracy for a wide range of taxonomic groups, and EGTs as short as 27 amino acids can be phylogenetically classified up to the rank of genus. The algorithm was applied in a comparative study of three aquatic microbial samples obtained by 454 pyrosequencing. Profound differences in the taxonomic composition of these samples could be clearly revealed.