Local Sequence Context Research Articles

Target Cell-Mediated Editing of HIV-1 cDNA by APOBEC3 Proteins in Human Macrophages

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) proteins are encapsidated by assembling HIV-1 virions and edit viral cDNA in the next round of infection. Using alpha interferon (IFN-α)-treated monocyte-derived macrophages, we show that infrequent editing of HIV-1 reverse transcripts can also be mediated by APOBEC3 proteins supplied by the targets of infection. Based on the local sequence contexts of these mutations and the established characteristics of APOBEC3 protein expression in myeloid cells, we speculate that APOBEC3A may be responsible for a substantial proportion of this activity.

The influence of the local sequence environment on RNA loop structures

RNA folding is assumed to be a hierarchical process. The secondary structure of an RNA molecule, signified by base-pairing and stacking interactions between the paired bases, is formed first. Subsequently, the RNA molecule adopts an energetically favorable three-dimensional conformation in the structural space determined mainly by the rotational degrees of freedom associated with the backbone of regions of unpaired nucleotides (loops). To what extent the backbone conformation of RNA loops also results from interactions within the local sequence context or rather follows global optimization constraints alone has not been addressed yet. Because the majority of base stacking interactions are exerted locally, a critical influence of local sequence on local structure appears plausible. Thus, local loop structure ought to be predictable, at least in part, from the local sequence context alone. To test this hypothesis, we used Random Forests on a nonredundant data set of unpaired nucleotides extracted from 97 X-ray structures from the Protein Data Bank (PDB) to predict discrete backbone angle conformations given by the discretized η/θ-pseudo-torsional space. Predictions on balanced sets with four to six conformational classes using local sequence information yielded average accuracies of up to 55%, thus significantly better than expected by chance (17%-25%). Bases close to the central nucleotide appear to be most tightly linked to its conformation. Our results suggest that RNA loop structure does not only depend on long-range base-pairing interactions; instead, it appears that local sequence context exerts a significant influence on the formation of the local loop structure.

Intrinsic Characteristics of Neighboring DNA Modulate Transposable Element Activity in Drosophila melanogaster

Identifying factors influencing transposable element activity is essential for understanding how these elements impact genomes and their evolution as well as for fully exploiting them as functional genomics tools and gene-therapy vectors. Using a genetics-based approach, the influence of genomic position on piggyBac mobility in Drosophila melanogaster was assessed while controlling for element structure, genetic background, and transposase concentration. The mobility of piggyBac elements varied over more than two orders of magnitude solely as a result of their locations within the genome. The influence of genomic position on element activities was independent of factors resulting in position-dependent transgene expression ("position effects"). Elements could be relocated to new genomic locations without altering their activity if ≥ 500 bp of genomic DNA originally flanking the element was also relocated. Local intrinsic factors within the neighboring DNA that determined the activity of piggyBac elements were portable not only within the genome but also when elements were moved to plasmids. The predicted bendability of the first 50 bp flanking the 5' and 3' termini of piggyBac elements could account for 60% of the variance in position-dependent activity observed among elements. These results are significant because positional influences on transposable element activities will impact patterns of accumulation of elements within genomes. Manipulating and controlling the local sequence context of piggyBac elements could be a powerful, novel way of optimizing gene vector activity.

SH2 Domains Recognize Contextual Peptide Sequence Information to Determine Selectivity

Selective ligand recognition by modular protein interaction domains is a primary determinant of specificity in signaling pathways. Src homology 2 (SH2) domains fulfill this capacity immediately downstream of tyrosine kinases, acting to recruit their host polypeptides to ligand proteins harboring phosphorylated tyrosine residues. The degree to which SH2 domains are selective and the mechanisms underlying selectivity are fundamental to understanding phosphotyrosine signaling networks. An examination of interactions between 50 SH2 domains and a set of 192 phosphotyrosine peptides corresponding to physiological motifs within FGF, insulin, and IGF-1 receptor pathways indicates that individual SH2 domains have distinct recognition properties and exhibit a remarkable degree of selectivity beyond that predicted by previously described binding motifs. The underlying basis for such selectivity is the ability of SH2 domains to recognize both permissive amino acid residues that enhance binding and non-permissive amino acid residues that oppose binding in the vicinity of the essential phosphotyrosine. Neighboring positions affect one another so local sequence context matters to SH2 domains. This complex linguistics allows SH2 domains to distinguish subtle differences in peptide ligands. This newly appreciated contextual dependence substantially increases the accessible information content embedded in the peptide ligands that can be effectively integrated to determine binding. This concept may serve more broadly as a paradigm for subtle recognition of physiological ligands by protein interaction domains.

Deep sequencing analysis of mutations resulting from the incorporation of dNTP analogs

Next-generation DNA sequencing technology was used to score >100 000 mutations resulting from exposure of a nucleic acid template to a mutagenic dNTP analog during a single pass of a DNA polymerase. An RNA template of known secondary structure was reverse transcribed in the presence of 8-oxo-dGTP, dPTP or both, followed by forward transcription in the presence of standard NTPs. Each mutagen, whether used alone or in combination, resulted in a highly characteristic mutation profile. Mutations were generated at a mean frequency of 1–2% per eligible nucleotide position, but there was substantial variation in the frequency of mutation at different positions, with a SD close to the mean. This variation was partly due to the identity of the immediately surrounding nucleotides and was not significantly influenced by the secondary structure of the RNA template. Most of the variation appears to result from idiosyncratic features that derive from local sequence context, demonstrating how different genetic sequences have different chemical phenotypes.

Comparison of the in Vitro Replication of the 7-(2-Oxoheptyl)-1,N2-etheno-2′-deoxyguanosine and 1,N2-Etheno-2′-deoxyguanosine Lesions by Sulfolobus solfataricus P2 DNA Polymerase IV (Dpo4)

Oligonucleotides were synthesized containing the 7-(2-oxoheptyl)-etheno-dGuo adduct, which is derived from the reaction of dGuo and the lipid peroxidation product 4-oxo-2-nonenal. The in vitro replication of 7-(2-oxoheptyl)-etheno-dGuo by the model Y-family polymerase Sulfolobus solfataricus P2 DNA Polymerase IV (Dpo4) was examined in two sequences. The extension products were sequenced using an improved LC-ESI-MS/MS protocol developed in our laboratories, and the results were compared to that of the 1,N(2)-etheno-dGuo adduct in the same sequence contexts. Both etheno adducts were highly miscoding when situated in 5'-TXG-3' local sequence contexts with <4% of the extension products being derived from error-free bypass. The major extension products resulted from the misinsertion of Ade opposite the adduct and a one-base deletion. The major extension products from replication of the etheno lesions in a 5'-CXG-3' local sequence context were the result of misinsertion of Ade, a one-base deletion, and error-free bypass. Other minor extension products were also identified. The 7-(2-oxoheptyl)-etheno-dGuo lesion resulted in a larger frequency of misinsertion of Ade, whereas the 1,N(2)-etheno-dGuo gave more of the one-base deletion product. Conformational studies of duplex DNA containing the 7-(2-oxoheptyl)-etheno-dGuo in a 5'-TXG-3' sequence context by NMR indicated the presence of a pH-dependent conformational transition, likely involving the glycosyl bond at the adducted guanosine; the pK(a) for this transition was lower than that observed for the 1,N(2)-epsilon-dGuo lesion. However, the 7-(2-oxoheptyl)-etheno-dGuo lesion, the complementary Cyt, and both flanking base pairs remained disordered at all pH values, which is attributed to the presence of the hydrophobic heptyl group of the 7-(2-oxoheptyl)-etheno-dGuo lesion. The altered pK(a) value and the structural disorder at the 7-(2-oxoheptyl)-etheno-dGuo lesion site, as compared to the same sequence containing the 1,N(2)-etheno-dGuo, may contribute to higher frequency of misinsertion of Ade.

The role of local secondary structure in the function of the trans-splicing motif of Brugia malayi

A 7-nt motif (the trans-splicing motif or TSM) was previously shown to be necessary and sufficient to direct trans-splicing of transgenic mRNAs in transgenic Brugia malayi embryos. Insertion of the TSM into two genes lacking a TSM homologue resulted in trans-splicing of transgenic mRNAs from one transgene but not the other, suggesting that local sequence context might affect TSM function. To test this hypothesis, constructs inserting the TSM into different positions of two B. malayi genes were tested for their ability to support trans-splicing of transgenic mRNAs. Transgenic mRNAs derived from constructs in which the insertion of the TSM did not result in a perturbation of the local predicted secondary structure were trans-spliced, while those in which the TSM perturbed the local secondary structure were not. These data suggest that local secondary structure plays a role in the ability of the TSM to direct trans-splicing.

Influence of local sequence context on damaged base conformation in human DNA polymerase ι: molecular dynamics studies of nucleotide incorporation opposite a benzo[a]pyrene-derived adenine lesion

Human DNA polymerase ι is a lesion bypass polymerase of the Y family, capable of incorporating nucleotides opposite a variety of lesions in both near error-free and error-prone bypass. With undamaged templating purines polymerase ι normally favors Hoogsteen base pairing. Polymerase ι can incorporate nucleotides opposite a benzo[a]pyrene-derived adenine lesion (dA*); while mainly error-free, the identity of misincorporated bases is influenced by local sequence context. We performed molecular modeling and molecular dynamics simulations to elucidate the structural basis for lesion bypass. Our results suggest that hydrogen bonds between the benzo[a]pyrenyl moiety and nearby bases limit the movement of the templating base to maintain the anti glycosidic bond conformation in the binary complex in a 5′-CAGA*TT-3′ sequence. This facilitates correct incorporation of dT via a Watson−Crick pair. In a 5′-TTTA*GA-3′ sequence the lesion does not form these hydrogen bonds, permitting dA* to rotate around the glycosidic bond to syn and incorporate dT via a Hoogsteen pair. With syn dA*, there is also an opportunity for increased misincorporation of dGTP. These results expand our understanding of the versatility and flexibility of polymerase ι and its lesion bypass functions in humans.

Ride the wavelet: A multiscale analysis of genomic contexts flanking small insertions and deletions

Recent studies have revealed that insertions and deletions (indels) are more different in their formation than previously assumed. What remains enigmatic is how the local DNA sequence context contributes to these differences. To investigate the relative impact of various molecular mechanisms to indel formation, we analyzed sequence contexts of indels in the non protein- or RNA-coding, nonrepetitive (NCNR) portion of the human genome. We considered small (<or=30-bp) indels occurring in the human lineage since its divergence from chimpanzee and used wavelet techniques to study, simultaneously for multiple scales, the spatial patterns of short sequence motifs associated with indel mutagenesis. In particular, we focused on motifs associated with DNA polymerase activity, topoisomerase cleavage, double-strand breaks (DSBs), and their repair. We came to the following conclusions. First, many motifs are characterized by unique enrichment profiles in the vicinity of indels vs. indel-free portions of the genome, verifying the importance of sequence context in indel mutagenesis. Second, only limited similarity in motif frequency profiles is evident flanking insertions vs. deletions, confirming differences in their mutagenesis. Third, substantial similarity in frequency profiles exists between pairs of individual motifs flanking insertions (and separately deletions), suggesting "cooperation" among motifs, and thus molecular mechanisms, during indel formation. Fourth, the wavelet analyses demonstrate that all these patterns are highly dependent on scale (the size of an interval considered). Finally, our results depict a model of indel mutagenesis comprising both replication and recombination (via repair of paused replication forks and site-specific recombination).

A CA+ Pair Adjacent to a Sheared GA or AA Pair Stabilizes Size-Symmetric RNA Internal Loops

RNA internal loops are often important sites for folding and function. Residues in internal loops can have pKa values shifted close to neutral pH because of the local structural environment. A series of RNA internal loops were studied at different pH by UV absorbance versus temperature melting experiments and imino proton nuclear magnetic resonance (NMR). A stabilizing CA pair forms at pH 7 in the and nearest neighbors when the CA pair is the first noncanonical pair (loop-terminal pair) in 3 × 3 nucleotide and larger size-symmetric internal loops. These and nearest neighbors, with CA adjacent to a closing Watson−Crick pair, are further stabilized when the pH is lowered from 7 to 5.5. The results are consistent with a significantly larger fraction (from ∼20% at pH 7 to ∼90% at pH 5.5) of adenines being protonated at the N1 position to form stabilizing wobble CA+ pairs adjacent to a sheared GA or AA pair. The noncanonical pair adjacent to the GA pair in can either stabilize or destabilize the loop, consistent with the sequence-dependent thermodynamics of GA pairs. No significant pH-dependent stabilization is found for most of the other nearest neighbor combinations involving CA pairs (e.g., and ), which is consistent with the formation of various nonwobble pairs observed in different local sequence contexts in crystal and NMR structures. A revised free-energy model, including stabilization by wobble CA+ pairs, is derived for predicting stabilities of medium-size RNA internal loops.

Sequence context-specific profiles for homology searching

Sequence alignment and database searching are essential tools in biology because a protein's function can often be inferred from homologous proteins. Standard sequence comparison methods use substitution matrices to find the alignment with the best sum of similarity scores between aligned residues. These similarity scores do not take the local sequence context into account. Here, we present an approach that derives context-specific amino acid similarities from short windows centered on each query sequence residue. Our results demonstrate that the sequence context contains much more information about the expected mutations than just the residue itself. By employing our context-specific similarities (CS-BLAST) in combination with NCBI BLAST, we increase the sensitivity more than 2-fold on a difficult benchmark set, without loss of speed. Alignment quality is likewise improved significantly. Furthermore, we demonstrate considerable improvements when applying this paradigm to sequence profiles: Two iterations of CSI-BLAST, our context-specific version of PSI-BLAST, are more sensitive than 5 iterations of PSI-BLAST. The paradigm for biological sequence comparison presented here is very general. It can replace substitution matrices in sequence- and profile-based alignment and search methods for both protein and nucleotide sequences.

Genome and gene alterations by insertions and deletions in the evolution of human and chimpanzee chromosome 22

BackgroundUnderstanding structure and function of human genome requires knowledge of genomes of our closest living relatives, the primates. Nucleotide insertions and deletions (indels) play a significant role in differentiation that underlies phenotypic differences between humans and chimpanzees. In this study, we evaluated distribution, evolutionary history, and function of indels found by comparing syntenic regions of the human and chimpanzee genomes.ResultsSpecifically, we identified 6,279 indels of 10 bp or greater in a ~33 Mb alignment between human and chimpanzee chromosome 22. After the exclusion of those in repetitive DNA, 1,429 or 23% of indels still remained. This group was characterized according to the local or genome-wide repetitive nature, size, location relative to genes, and other genomic features. We defined three major classes of these indels, using local structure analysis: (i) those indels found uniquely without additional copies of indel sequence in the surrounding (10 Kb) region, (ii) those with at least one exact copy found nearby, and (iii) those with similar but not identical copies found locally. Among these classes, we encountered a high number of exactly repeated indel sequences, most likely due to recent duplications. Many of these indels (683 of 1,429) were in proximity of known human genes. Coding sequences and splice sites contained significantly fewer of these indels than expected from random expectations, suggesting that selection is a factor in limiting their persistence. A subset of indels from coding regions was experimentally validated and their impacts were predicted based on direct sequencing in several human populations as well as chimpanzees, bonobos, gorillas, and two subspecies of orangutans.ConclusionOur analysis demonstrates that while indels are distributed essentially randomly in intergenic and intronic genomic regions, they are significantly under-represented in coding sequences. There are substantial differences in representation of indel classes among genomic elements, most likely caused by differences in their evolutionary histories. Using local sequence context, we predicted origins and phylogenetic relationships of gene-impacting indels in primate species. These results suggest that genome plasticity is a major force behind speciation events separating the great ape lineages.

The Sequence Dependence of Human Nucleotide Excision Repair Efficiencies of Benzo[ a]pyrene-derived DNA Lesions: Insights into the Structural Factors that Favor Dual Incisions

Nucleotide excision repair (NER) is a vital cellular defense system against carcinogen–DNA adducts, which, if not repaired, can initiate cancer development. The structural features of bulky DNA lesions that account for differences in NER efficiencies in mammalian cells are not well understood. In vivo, the predominant DNA adduct derived from metabolically activated benzo[ a]pyrene (BP), a prominent environmental carcinogen, is the 10 S (+)- trans- anti-[BP]- N 2-dG adduct ( G*), which resides in the B-DNA minor groove 5′-oriented along the modified strand. We have compared the structural distortions in double-stranded DNA, imposed by this adduct, in the different sequence contexts 5′-…CG G*C…, 5′-…C G*GC…, 5′-…CI G*C… (I is 2′-deoxyinosine), and 5′-…C G*C…. On the basis of electrophoretic mobilities, all duplexes manifest moderate bends, except the 5′-…CG G*C…duplex, which exhibits an anomalous, slow mobility attributed to a pronounced flexible kink at the site of the lesion. This kink, resulting from steric hindrance between the 5′-flanking guanine amino group and the BP aromatic rings, both positioned in the minor groove, is abolished in the 5′-…CI G*C…duplex (the 2′-deoxyinosine group, I, lacks this amino group). In contrast, the sequence-isomeric 5′-…C G*GC…duplex exhibits only a moderate bend, but displays a remarkably increased opening rate at the 5′-flanking base pair of G*, indicating a significant destabilization of Watson–Crick hydrogen bonding. The NER dual incision product yields were compared for these different sequences embedded in otherwise identical 135-mer duplexes in cell-free human HeLa extracts. The yields of excision products varied by a factor of as much as ∼ 4 in the order 5′-...C G*GC…> 5′...CG G*C…≥ 5′...CI G*C…≥ 5′-…C G*C…. Overall, destabilized Watson–Crick hydrogen bonding, manifested in the 5′-...C G*GC...duplex, elicits the most significant NER response, while the flexible kink displayed in the sequence-isomeric 5′-...CG G*C...duplex represents a less significant signal in this series of substrates. These results demonstrate that the identical lesion can be repaired with markedly variable efficiency in different local sequence contexts that differentially alter the structural features of the DNA duplex around the lesion site.

Evaluation of next generation sequencing platforms for population targeted sequencing studies

Human sequence generated from three next-generation sequencing platforms reveals systematic variability in sequence coverage due to local sequence characteristics.

A Multifactorial Signature of DNA Sequence and Polycomb Binding Predicts Aberrant CpG Island Methylation

Aberrant CpG island methylation is associated with transcriptional silencing of regulatory genes in human cancer. Although most CpG islands remain unmethylated, a subset accrues aberrant methylation in cancer via unknown mechanisms. Previously, we showed that CpG islands differ in their intrinsic propensity towards hypermethylation. We developed a classifier (PatMAn) based on the frequencies of seven DNA sequence patterns that discriminated methylation-prone (MP) and methylation-resistant (MR) CpG islands. Here, we report on the genome-wide application and direct testing of PatMAn in cancer. Although trained on data from a cell culture model of de novo methylation involving the overexpression of DNMT1, PatMAn accurately predicted CpG islands at increased risk of hypermethylation in cancer cell lines and primary tumors. Analysis of CpG islands predicted to be MP revealed a strong association with embryonic targets of polycomb-repressive complex 2 (PRC2), indicating that PatMAn predicts not only aberrant methylation, but also PRC2 binding. A second classifier (SUPER-PatMAn) that integrates the seven PatMAn DNA patterns with SUZ12 enriched regions as a marker of PRC2 occupancy showed improved performance (prediction accuracy, 81-88%). In addition to many non-PRC2 targets, SUPER-PatMAn identified a subset of PRC2 targets that were more likely to be hypermethylated in cancer. Genome-wide, CpG islands predicted to be MP were enriched in genes known to undergo hypermethylation in cancer, genes functioning in transcriptional regulation, and components of developmental pathways. These findings show that hypermethylation of certain gene loci is controlled in part by an underlying susceptibility influenced by both local sequence context and trans-acting factors.

Hypermutable non-synonymous sites are under stronger negative selection.

Mutation rate varies greatly between nucleotide sites of the human genome and depends both on the global genomic location and the local sequence context of a site. In particular, CpG context elevates the mutation rate by an order of magnitude. Mutations also vary widely in their effect on the molecular function, phenotype, and fitness. Independence of the probability of occurrence of a new mutation's effect has been a fundamental premise in genetics. However, highly mutable contexts may be preserved by negative selection at important sites but destroyed by mutation at sites under no selection. Thus, there may be a positive correlation between the rate of mutations at a nucleotide site and the magnitude of their effect on fitness. We studied the impact of CpG context on the rate of human–chimpanzee divergence and on intrahuman nucleotide diversity at non-synonymous coding sites. We compared nucleotides that occupy identical positions within codons of identical amino acids and only differ by being within versus outside CpG context. Nucleotides within CpG context are under a stronger negative selection, as revealed by their lower, proportionally to the mutation rate, rate of evolution and nucleotide diversity. In particular, the probability of fixation of a non-synonymous transition at a CpG site is two times lower than at a CpG site. Thus, sites with different mutation rates are not necessarily selectively equivalent. This suggests that the mutation rate may complement sequence conservation as a characteristic predictive of functional importance of nucleotide sites.

Polypyrimidine tract binding protein regulates alternative splicing of an aberrant pseudoexon in NF1

In disease-associated genes, understanding the functional significance of deep intronic nucleotide variants represents a difficult challenge. We previously reported that an NF1 intron 30 exonization event is triggered from a single correct nomenclature is 'c.293-279 A>G' mutation [Raponi M, Upadhyaya M & Baralle D (2006) Hum Mutat 27, 294-295]. In this paper, we investigate which characteristics play a role in regulating inclusion of the aberrant pseudoexon. Our investigation shows that pseudoexon inclusion levels are strongly downregulated by polypyrimidine tract binding protein and its homologue neuronal polypyrimidine tract binding protein. In particular, we provide evidence that the functional effect of polypyrimidine tract binding protein is proportional to its concentration, and map the cis-acting elements that are principally responsible for this negative regulation. These results highlight the importance of evaluating local sequence context for diagnostic purposes, and the utility of developing therapies to turn off activated pseudoexons.

Genome‐wide transposon tagging reveals location‐dependent effects on transcription and chromatin organization in Arabidopsis

The interphase nucleus exists as a highly dynamic system, the physical properties of which have functional importance in gene regulation. Not only can gene expression be influenced by the local sequence context, but also by the architecture of the nucleus in three-dimensions (3D), and by the interactions between these levels via chromatin modifications. A challenging task is to resolve the complex interplay between sequence- and genome structure-based control mechanisms. Here, we created a collection of 277 Arabidopsis lines that allow the visual tracking of individual loci in living plants while comparing gene expression potential at these locations, via an identical reporter cassette. Our studies revealed regional gene silencing near a heterochromatin island, via DNA methylation, that is correlated with mobility constraint and nucleolar association. We also found an example of nucleolar association that does not correlate with gene suppression, suggesting that distinct mechanisms exist that can mediate interactions between chromatin and the nucleolus. These studies demonstrate the utility of this novel resource in unifying structural and functional studies towards a more comprehensive model of how global chromatin organization may coordinate gene expression over large scales.

DNA polymerase bypass in vitro and in E. coli of a C-nucleotide analogue of Fapy·dG

Bypass of the configurationally stable analogue (β-C-Fapy·dG) of the formamidopyrimidine lesion derived from 2′-deoxyguanosine oxidation (Fapy·dG) was studied in vitro and in Escherichia coli. The exonuclease deficient Klenow fragment of E. coli DNA polymerase I (Klenow exo −) misincorporated dA most frequently opposite β-C-Fapy·dG, but its efficiency was <0.2% of dC insertion. Klenow exo − fidelity was enhanced by the enzyme’s high selectivity for extending duplexes only when dC was opposite β-C-Fapy·dG. The expectations raised by these in vitro data were realized when β-C-Fapy·dG replication was studied in E. coli by transfecting M13mp7(L2) bacteriophage DNA containing the nucleotide analogue within the lacZ gene in 4 local sequence contexts. The bypass efficiency of β-C-Fapy·dG varied between 45% and 70% compared to a genome containing only native nucleotides. Mutation frequencies at the site of the lesions in the originally transfected genomes were determined using the REAP assay [Delaney, J. C.; Essigmann, J. M. Methods Enzymol. 2006, 408, 1]. The levels of mutations could not be distinguished between those observed when genomes containing native nucleotides were replicated, indicating that the mutagenicity of β-C-Fapy·dG was <1%. These data and previous reports indicate that β-C-Fapy·dG is a good model of Fapy·dG in E. coli. In addition, these results and the previous report of β-C-Fapy·dG binding to the base excision repair protein formamidopyrimidine glycosylase suggest that this analogue could be useful as a DNA repair inhibitor.

Genome wide prediction of HNF4α functional binding sites by the use of local and global sequence context

An application of machine learning algorithms enables prediction of the functional context of transcription factor binding sites in the human genome.