DNA Sequence Motifs Research Articles

Identification of genomic features in environmentally induced epigenetic transgenerational inherited sperm epimutations.

A variety of environmental toxicants have been shown to induce the epigenetic transgenerational inheritance of disease and phenotypic variation. The process involves exposure of a gestating female and the developing fetus to environmental factors that promote permanent alterations in the epigenetic programming of the germline. The molecular aspects of the phenomenon involve epigenetic modifications (epimutations) in the germline (e.g. sperm) that are transmitted to subsequent generations. The current study integrates previously described experimental epigenomic transgenerational data and web-based bioinformatic analyses to identify genomic features associated with these transgenerationally transmitted epimutations. A previously identified genomic feature associated with these epimutations is a low CpG density (<12/100bp). The current observations suggest the transgenerational differential DNA methylation regions (DMR) in sperm contain unique consensus DNA sequence motifs, zinc finger motifs and G-quadruplex sequences. Interaction of molecular factors with these sequences could alter chromatin structure and accessibility of proteins with DNA methyltransferases to alter de novo DNA methylation patterns. G-quadruplex regions can promote the opening of the chromatin that may influence the action of DNA methyltransferases, or factors interacting with them, for the establishment of epigenetic marks. Zinc finger binding factors can also promote this chromatin remodeling and influence the expression of non-coding RNA. The current study identified genomic features associated with sperm epimutations that may explain in part how these sites become susceptible for transgenerational programming.

Application of the Codon-based Scoring Method in Motif Detection

A RNA or DNA sequence motif is a short sequence found within a particular nucleic acid sequence families. Most amino acid and nucleic acid sequences have some level of functional or structural similarities. These similarities are mostly represented by short, contiguous sequences called motif. Motif discovery is an important aspect of molecular biology. This is because the knowledge of these sequences helps determine their structural properties, signal sites and/or ligand-binding sites. In most cases, depending on the function of the motif, these contiguous regions can be highly conserved with an homology of nearly 100%. Several algorithms have been proposed for the discovery of motifs. In this paper, the codon-based scoring method is employed to detect motifs and with their invariants. The result obtained shows the reliability and robustness of the method as motifs are discovered irrespective of their length and position in a sequence.

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.

Gene expression signatures affected by alcohol-induced DNA methylomic deregulation in human embryonic stem cells

Stem cells, especially human embryonic stem cells (hESCs), are useful models to study molecular mechanisms of human disorders that originate during gestation. Alcohol (ethanol, EtOH) consumption during pregnancy causes a variety of prenatal and postnatal disorders collectively referred to as fetal alcohol spectrum disorders (FASDs). To better understand the molecular events leading to FASDs, we performed a genome-wide analysis of EtOH's effects on the maintenance and differentiation of hESCs in culture. Gene Co-expression Network Analysis showed significant alterations in gene profiles of EtOH-treated differentiated or undifferentiated hESCs, particularly those associated with molecular pathways for metabolic processes, oxidative stress, and neuronal properties of stem cells. A genome-wide DNA methylome analysis revealed widespread EtOH-induced alterations with significant hypermethylation of many regions of chromosomes. Undifferentiated hESCs were more vulnerable to EtOH's effect than their differentiated counterparts, with methylation on the promoter regions of chromosomes 2, 16 and 18 in undifferentiated hESCs most affected by EtOH exposure. Combined transcriptomic and DNA methylomic analysis produced a list of differentiation-related genes dysregulated by EtOH-induced DNA methylation changes, which likely play a role in EtOH-induced decreases in hESC pluripotency. DNA sequence motif analysis of genes epigenetically altered by EtOH identified major motifs representing potential binding sites for transcription factors. These findings should help in deciphering the precise mechanisms of alcohol-induced teratogenesis.

Methylome diversification through changes in DNA methyltransferase sequence specificity.

Epigenetic modifications such as DNA methylation have large effects on gene expression and genome maintenance. Helicobacter pylori, a human gastric pathogen, has a large number of DNA methyltransferase genes, with different strains having unique repertoires. Previous genome comparisons suggested that these methyltransferases often change DNA sequence specificity through domain movement—the movement between and within genes of coding sequences of target recognition domains. Using single-molecule real-time sequencing technology, which detects N6-methyladenines and N4-methylcytosines with single-base resolution, we studied methylated DNA sites throughout the H. pylori genome for several closely related strains. Overall, the methylome was highly variable among closely related strains. Hypermethylated regions were found, for example, in rpoB gene for RNA polymerase. We identified DNA sequence motifs for methylation and then assigned each of them to a specific homology group of the target recognition domains in the specificity-determining genes for Type I and other restriction-modification systems. These results supported proposed mechanisms for sequence-specificity changes in DNA methyltransferases. Knocking out one of the Type I specificity genes led to transcriptome changes, which suggested its role in gene expression. These results are consistent with the concept of evolution driven by DNA methylation, in which changes in the methylome lead to changes in the transcriptome and potentially to changes in phenotype, providing targets for natural or artificial selection.

Sequence, distribution and chromosomal context of class I and class II pilin genes of Neisseria meningitidis identified in whole genome sequences

BackgroundNeisseria meningitidis expresses type four pili (Tfp) which are important for colonisation and virulence. Tfp have been considered as one of the most variable structures on the bacterial surface due to high frequency gene conversion, resulting in amino acid sequence variation of the major pilin subunit (PilE). Meningococci express either a class I or a class II pilE gene and recent work has indicated that class II pilins do not undergo antigenic variation, as class II pilE genes encode conserved pilin subunits. The purpose of this work was to use whole genome sequences to further investigate the frequency and variability of the class II pilE genes in meningococcal isolate collections.ResultsWe analysed over 600 publically available whole genome sequences of N. meningitidis isolates to determine the sequence and genomic organization of pilE. We confirmed that meningococcal strains belonging to a limited number of clonal complexes (ccs, namely cc1, cc5, cc8, cc11 and cc174) harbour a class II pilE gene which is conserved in terms of sequence and chromosomal context. We also identified pilS cassettes in all isolates with class II pilE, however, our analysis indicates that these do not serve as donor sequences for pilE/pilS recombination. Furthermore, our work reveals that the class II pilE locus lacks the DNA sequence motifs that enable (G4) or enhance (Sma/Cla repeat) pilin antigenic variation. Finally, through analysis of pilin genes in commensal Neisseria species we found that meningococcal class II pilE genes are closely related to pilE from Neisseria lactamica and Neisseria polysaccharea, suggesting horizontal transfer among these species.ConclusionsClass II pilins can be defined by their amino acid sequence and genomic context and are present in meningococcal isolates which have persisted and spread globally. The absence of G4 and Sma/Cla sequences adjacent to the class II pilE genes is consistent with the lack of pilin subunit variation in these isolates, although horizontal transfer may generate class II pilin diversity. This study supports the suggestion that high frequency antigenic variation of pilin is not universal in pathogenic Neisseria.

High-Density Nucleosome Occupancy Map of Human Chromosome 9p21–22 Reveals Chromatin Organization of the Type I Interferon Gene Cluster

Genome-wide investigations have dramatically increased our understanding of nucleosome positioning and the role of chromatin in gene regulation, yet some genomic regions have been poorly represented in human nucleosome maps. One such region is represented by human chromosome 9p21-22, which contains the type I interferon gene cluster that includes 16 interferon alpha genes and the single interferon beta, interferon epsilon, and interferon omega genes. A high-density nucleosome mapping strategy was used to generate locus-wide maps of the nucleosome organization of this biomedically important locus at a steady state and during a time course of infection with Sendai virus, an inducer of interferon gene expression. Detailed statistical and computational analysis illustrates that nucleosomes in this locus exhibit preferences for particular dinucleotide and oligomer DNA sequence motifs in vivo, which are similar to those reported for lower eukaryotic nucleosome-DNA interactions. These data were used to visualize the region's chromatin architecture and reveal features that are common to the organization of all the type I interferon genes, indicating a common nucleosome-mediated gene regulatory paradigm. Additionally, this study clarifies aspects of the dynamic changes that occur with the nucleosome occupying the transcriptional start site of the interferon beta gene after virus infection.

A step‐by‐step protocol for formaldehyde‐assisted isolation of regulatory elements from Arabidopsis thaliana

The control of gene expression by transcriptional regulators and other types of functionally relevant DNA transactions such as chromatin remodeling and replication underlie a vast spectrum of biological processes in all organisms. DNA transactions require the controlled interaction of proteins with DNA sequence motifs which are often located in nucleosome-depleted regions (NDRs) of the chromatin. Formaldehyde-assisted isolation of regulatory elements (FAIRE) has been established as an easy-to-implement method for the isolation of NDRs from a number of eukaryotic organisms, and it has been successfully employed for the discovery of new regulatory segments in genomic DNA from, for example, yeast, Drosophila, and humans. Until today, however, FAIRE has only rarely been employed in plant research and currently no detailed FAIRE protocol for plants has been published. Here, we provide a step-by-step FAIRE protocol for NDR discovery in Arabidopsis thaliana. We demonstrate that NDRs isolated from plant chromatin are readily amenable to quantitative polymerase chain reaction and next-generation sequencing. Only minor modification of the FAIRE protocol will be needed to adapt it to other plants, thus facilitating the global inventory of regulatory regions across species.

Finding Motifs in DNA Sequences Using Low-Dispersion Sequences

Motif finding problems, abstracted as the planted (l, d)-motif finding problem, are a major task in molecular biology--finding functioning units and genes. In 2002, the random projection algorithm was introduced to solve the challenging (15, 4)-motif finding problem by using randomly chosen templates. Two years later, a so-called uniform projection algorithm was developed to improve the random projection algorithm by means of low-dispersion sequences generated by coverings. In this article, we introduce an improved projection algorithm called the low-dispersion projection algorithm, which uses low-dispersion sequences generated by developed almost difference families. Compared with the random projection algorithm, the low-dispersion projection algorithm can solve the (l, d)-motif finding problem with fewer templates without decreasing the success rate.

LDsplit: screening for cis-regulatory motifs stimulating meiotic recombination hotspots by analysis of DNA sequence polymorphisms

BackgroundAs a fundamental genomic element, meiotic recombination hotspot plays important roles in life sciences. Thus uncovering its regulatory mechanisms has broad impact on biomedical research. Despite the recent identification of the zinc finger protein PRDM9 and its 13-mer binding motif as major regulators for meiotic recombination hotspots, other regulators remain to be discovered. Existing methods for finding DNA sequence motifs of recombination hotspots often rely on the enrichment of co-localizations between hotspots and short DNA patterns, which ignore the cross-individual variation of recombination rates and sequence polymorphisms in the population. Our objective in this paper is to capture signals encoded in genetic variations for the discovery of recombination-associated DNA motifs.ResultsRecently, an algorithm called “LDsplit” has been designed to detect the association between single nucleotide polymorphisms (SNPs) and proximal meiotic recombination hotspots. The association is measured by the difference of population recombination rates at a hotspot between two alleles of a candidate SNP. Here we present an open source software tool of LDsplit, with integrative data visualization for recombination hotspots and their proximal SNPs. Applying LDsplit on SNPs inside an established 7-mer motif bound by PRDM9 we observed that SNP alleles preserving the original motif tend to have higher recombination rates than the opposite alleles that disrupt the motif. Running on SNP windows around hotspots each containing an occurrence of the 7-mer motif, LDsplit is able to guide the established motif finding algorithm of MEME to recover the 7-mer motif. In contrast, without LDsplit the 7-mer motif could not be identified.ConclusionsLDsplit is a software tool for the discovery of cis-regulatory DNA sequence motifs stimulating meiotic recombination hotspots by screening and narrowing down to hotspot associated SNPs. It is the first computational method that utilizes the genetic variation of recombination hotspots among individuals, opening a new avenue for motif finding. Tested on an established motif and simulated datasets, LDsplit shows promise to discover novel DNA motifs for meiotic recombination hotspots.

Erk1/2 Activity Promotes Chromatin Features and RNAPII Phosphorylation at Developmental Promoters in Mouse ESCs

Erk1/2 activation contributes to mouse ES cell pluripotency. We found a direct role of Erk1/2 in modulating chromatin features required for regulated developmental gene expression. Erk2 binds to specific DNA sequence motifs typically accessed by Jarid2 and PRC2. Negating Erk1/2 activation leads to increased nucleosome occupancy and decreased occupancy of PRC2 and poised RNAPII at Erk2-PRC2-targeted developmental genes. Surprisingly, Erk2-PRC2-targeted genes are specifically devoid of TFIIH, known to phosphorylate RNA polymerase II (RNAPII) at serine-5, giving rise to its initiated form. Erk2 interacts with and phosphorylates RNAPII at its serine 5 residue, which is consistent with the presence of poised RNAPII as a function of Erk1/2 activation. These findings underscore a key role for Erk1/2 activation in promoting the primed status of developmental genes in mouse ES cells and suggest that the transcription complex at developmental genes is different than the complexes formed at other genes, offering alternative pathways of regulation.

Transcription Factor Ets1, but Not the Closely Related Factor Ets2, Inhibits Antibody-Secreting Cell Differentiation

B cell differentiation into antibody-secreting cells (ASCs) is a tightly regulated process under the control of multiple transcription factors. One such transcription factor, Ets1, blocks the transition of B cells to ASCs via two separate activities: (i) stimulating the expression of target genes that promote B cell identity and (ii) interfering with the functional activity of the transcription factor Blimp1. Ets1 is a member of a multigene family, several members of which are expressed within the B cell lineage, including the closely related protein Ets2. In this report, we demonstrate that Ets1, but not Ets2, can block ASC formation despite the fact that Ets1 and Ets2 bind to apparently identical DNA sequence motifs and are thought to regulate overlapping sets of target genes. The DNA binding domain of Ets1 is required, but not sufficient by itself, to block ASC formation. In addition, less conserved regions within the N terminus of Ets1 play an important role in inhibiting B cell differentiation. Differences between the N termini of Ets1 and Ets2, rather than differences in the DNA binding domains, determine whether the proteins are capable of blocking ASC formation or not.

Transcriptional circuitry underlying seed coat development in Arabidopsis

We analyzed two sub-regions of the maternal seed coat, chalazal (CZSC) and distal (SC), using transcriptomic and histological analyses in the model plant Arabidopsis thaliana. Hierarchical clustering analysis showed that the CZSC and SC are transcriptionally distinct, though the two sub-regions are more similar during early stages of seed development. Robust statistical and network analysis revealed novel roles for both sub-regions during the course of the seed lifecycle and provides insight into the regulatory circuitry underlying these poorly studied sub-regions of the seed. Data show many of the processes that characterize the SC including starch deposition during the morphogenesis phase, and mucilage deposition and cell wall thickening during the maturation phase, are either absent or expressed to a much lesser extent in the CZSC. We further analyzed the CZSC in detail and show that this sub-region is likely involved in the control of information into the seed from the maternal plant and that some of these processes are predicted to operate through the activity of bZIP transcription factors through the G-box DNA sequence motif.

A new exhaustive method and strategy for finding motifs in ChIP-enriched regions.

ChIP-seq, which combines chromatin immunoprecipitation (ChIP) with next-generation parallel sequencing, allows for the genome-wide identification of protein-DNA interactions. This technology poses new challenges for the development of novel motif-finding algorithms and methods for determining exact protein-DNA binding sites from ChIP-enriched sequencing data. State-of-the-art heuristic, exhaustive search algorithms have limited application for the identification of short (, ) motifs (, ) contained in ChIP-enriched regions. In this work we have developed a more powerful exhaustive method (FMotif) for finding long (, ) motifs in DNA sequences. In conjunction with our method, we have adopted a simple ChIP-enriched sampling strategy for finding these motifs in large-scale ChIP-enriched regions. Empirical studies on synthetic samples and applications using several ChIP data sets including 16 TF (transcription factor) ChIP-seq data sets and five TF ChIP-exo data sets have demonstrated that our proposed method is capable of finding these motifs with high efficiency and accuracy. The source code for FMotif is available at http://211.71.76.45/FMotif/.

Preferred binding of gain-of-function mutant p53 to bidirectional promoters with coordinated binding of ETS1 and GABPA to multiple binding sites

Gain-of-function mutant p53 is thought to induce gene expression in part by binding transcription factors bound to promoters for genes that mediate oncogenesis. We investigated the mechanism of mutant p53 binding by mapping the human genomic binding sites for p53 R273H using ChIP-Seq and showed them to localize to ETS DNA sequence motifs and locations with ETS1 and GABPA binding, both within promoters and distal to promoters. Strikingly, p53 R273H showed statistically significant and substantial binding to bidirectional promoters, which are enriched for inverted repeated ETS DNA sequence motifs. p53 R273H exhibited an exponential increase in probability of binding promoters with a higher number of ETS motifs. Both ETS1 and GABPA also showed an increase in the probability of binding to promoters with a higher number of ETS motifs. However, despite this increase in probability of binding by p53 R273H and ETS1, there was no increase in the binding signal, suggesting that the number of ETS1 and p53 R273H proteins bound per promoter is being limited. In contrast, GABPA did exhibit an increase in binding signal with higher numbers of ETS motifs per promoter. Analysis of the distance between inverted pairs of ETS motifs within promoters and binding by p53 R273H, ETS1 and GABPA, showed a novel coordination of binding for the three proteins. Both ETS1 and p53 R273H exhibited preference for binding promoters with distantly spaced ETS motifs in face-to-face and back-to-back orientations, and low binding preference to promoters with closely spaced ETS motifs. GABPA exhibited the inverse pattern of binding by preferring to bind promoters with closely spaced ETS motifs. Analysis of the helical phase between ETS motifs showed that ETS1 and p53 R273H exhibited a low preference for binding promoters with ETS motifs on the same face of the DNA helix. We propose a model for the binding of ETS1 and p53 R273H in which two inverted ETS motifs on a looped DNA helix are juxtaposed for ETS1 binding as a homodimer, with p53 R273H bound to ETS1. We propose that the formation of this DNA loop and protein-bound complex prevents additional binding of ETS1 and p53 R273H proteins to other proximal binding sites.

C-Jun/c-Fos heterodimers regulate cellular genes via a newly identified class of methylated DNA sequence motifs

CpG methylation in mammalian DNA is known to interfere with gene expression by inhibiting the binding of transactivators to their cognate sequence motifs or recruiting proteins involved in gene repression. An Epstein–Barr virus-encoded transcription factor, Zta, was the first example of a sequence-specific transcription factor that preferentially recognizes and selectively binds DNA sequence motifs with methylated CpG residues, reverses epigenetic silencing and activates gene transcription. The DNA binding domain of Zta is homologous to c-Fos, a member of the cellular AP-1 (activator protein 1) transcription factor family, which regulates cell proliferation and survival, apoptosis, transformation and oncogenesis. We have identified a novel AP-1 binding site termed meAP-1, which contains a CpG dinucleotide. If methylated, meAP-1 sites are preferentially bound by the AP-1 heterodimer c-Jun/c-Fos in vitro and in cellular chromatin in vivo. In activated human primary B cells, c-Jun/c-Fos locates to these methylated elements in promoter regions of transcriptionally activated genes. Reminiscent of the viral Zta protein, c-Jun/c-Fos is the first identified cellular member of the AP-1 family of transactivators that can induce expression of genes with methylated, hence repressed promoters, reversing epigenetic silencing.

Education in Computational Biology Today and Tomorrow

Education in Computational Biology Today and Tomorrow

WhiB7, an Fe-S-dependent Transcription Factor That Activates Species-specific Repertoires of Drug Resistance Determinants in Actinobacteria

WhiB-like (Wbl) proteins are well known for their diverse roles in actinobacterial morphogenesis, cell division, virulence, primary and secondary metabolism, and intrinsic antibiotic resistance. Gene disruption experiments showed that three different Actinobacteria (Mycobacterium smegmatis, Streptomyces lividans, and Rhodococcus jostii) each exhibited a different whiB7-dependent resistance profile. Heterologous expression of whiB7 genes showed these resistance profiles reflected the host's repertoire of endogenous whiB7-dependent genes. Transcriptional activation of two resistance genes in the whiB7 regulon, tap (a multidrug transporter) and erm(37) (a ribosomal methyltransferase), required interaction of WhiB7 with their promoters. Furthermore, heterologous expression of tap genes isolated from Mycobacterium species demonstrated that divergencies in drug specificity of homologous structural proteins contribute to the variation of WhiB7-dependent drug resistance. WhiB7 has a specific tryptophan/glycine-rich region and four conserved cysteine residues; it also has a peptide sequence (AT-hook) at its C terminus that binds AT-rich DNA sequence motifs upstream of the promoters it activates. Targeted mutagenesis showed that these motifs were required to provide antibiotic resistance in vivo. Anaerobically purified WhiB7 from S. lividans was dimeric and contained 2.1 ± 0.3 and 2.2 ± 0.3 mol of iron and sulfur, respectively, per protomer (consistent with the presence of a 2Fe-2S cluster). However, the properties of the dimer's absorption spectrum were most consistent with the presence of an oxygen-labile 4Fe-4S cluster, suggesting 50% occupancy. These data provide the first insights into WhiB7 iron-sulfur clusters as they exist in vivo, a major unresolved issue in studies of Wbl proteins.

Detecting Motifs in DNA Sequences by Branching from Neighbors of Qualified Potential Motifs

Detecting Motifs in DNA Sequences by Branching from Neighbors of Qualified Potential Motifs

Meganuclease-mediated Virus Self-cleavage Facilitates Tumor-specific Virus Replication

Meganucleases can specifically cleave long DNA sequence motifs, a feature that makes them an ideal tool for gene engineering in living cells. In a proof-of-concept study, we investigated the use of the meganuclease I-Sce I for targeted virus self-disruption to generate high-specific oncolytic viruses. For this purpose, we provided oncolytic adenoviruses with a molecular circuit that selectively responds to p53 activation by expression of I-Sce I subsequently leading to self-disruption of the viral DNA via heterologous I-Sce I recognition sites within the virus genome. We observed that virus replication and cell lysis was effectively impaired in p53-normal cells, but not in p53-dysfunctional tumor cells. I-Sce I activity led to effective intracellular processing of viral DNA as confirmed by detection of specific cleavage products. Virus disruption did not interfere with E1A levels indicating that reduction of functional virus genomes was the predominant cause for conditional replication. Consequently, tumor-specific replication was further enhanced when E1A expression was additionally inhibited by targeted transcriptional repression. Finally, we demonstrated p53-dependent oncolysis by I-Sce I-expressing viruses in vitro and in vivo, and demonstrated effective inhibition of tumor growth. In summary, meganuclease-mediated virus cleavage represents a promising approach to provide oncolytic viruses with attractive safety profiles.