Nucleotide Context Research Articles

Computational analysis of auxin responsive elements in the Arabidopsis thaliana L. genome.

Auxin responsive elements (AuxRE) were found in upstream regions of target genes for ARFs (Auxin response factors). While Chip-seq data for most of ARFs are still unavailable, prediction of potential AuxRE is restricted by consensus models that detect too many false positive sites. Using sequence analysis of experimentally proven AuxREs, we revealed both an extended nucleotide context pattern for AuxRE itself and three distinct types of its coupling motifs (Y-patch, AuxRE-like, and ABRE-like), which together with AuxRE may form the composite elements. Computational analysis of the genome-wide distribution of the predicted AuxREs and their impact on auxin responsive gene expression allowed us to conclude that: (1) AuxREs are enriched around the transcription start site with the maximum density in 5'UTR; (2) AuxREs mediate auxin responsive up-regulation, not down-regulation. (3) Directly oriented single AuxREs and reverse multiple AuxREs are mostly associated with auxin responsiveness. In the composite AuxRE elements associated with auxin response, ABRE-like and Y-patch are 5'-flanking or overlapping AuxRE, whereas AuxRE-like motif is 3'-flanking. The specificity in location and orientation of the coupling elements suggests them as potential binding sites for ARFs partners.

Computational approach for calculating the probability of eukaryotic translation initiation from ribo-seq data that takes into account leaky scanning.

BackgroundRibosome profiling (ribo-seq) provides experimental data on the density of elongating or initiating ribosomes at the whole transcriptome level that can be potentially used for estimating absolute levels of translation initiation at individual Translation Initiation Sites (TISs). These absolute levels depend on the mutual organisation of TISs within individual mRNAs. For example, according to the leaky scanning model of translation initiation in eukaryotes, a strong TIS downstream of another strong TIS is unlikely to be productive, since only a few scanning ribosomes would be able to reach the downstream TIS. In order to understand the dependence of translation initiation efficiency on the surrounding nucleotide context, it is important to estimate the strength of TISs independently of their mutual organisation, i.e. to estimate with what probability a ribosome would initiate at a particular TIS.ResultsWe designed a simple computational approach for estimating the probabilities of ribosomes initiating at individual start codons using ribosome profiling data. The method is based on the widely accepted leaky scanning model of translation initiation in eukaryotes which postulates that scanning ribosomes may skip a start codon if the initiation context is unfavourable and continue on scanning. We tested our approach on three independent ribo-seq datasets obtained in mammalian cultured cells.ConclusionsOur results suggested that the method successfully discriminates between weak and strong TISs and that the majority of numerous non-AUG TISs reported recently are very weak. Therefore the high frequency of non-AUG TISs observed in ribosome profiling experiments is due to their proximity to mRNA 5′-ends rather than their strength. Detectable translation initiation at non-AUG codons downstream of AUG codons is comparatively infrequent. The leaky scanning method will be useful for the characterization of differences in start codon selection between tissues, developmental stages and in response to stress conditions.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-014-0380-4) contains supplementary material, which is available to authorized users.

Abstract 2401: Identification of new target genes in microsatellite unstable colorectal cancer by exome sequencing

Abstract Colorectal cancer (CRC) is the third most common cancer in Western countries. Approximately 15% of CRCs display microsatellite instability (MSI) caused by defective cellular mismatch repair. Most of the MSI CRCs are sporadic and result from biallelic inactivation of the MLH1 gene, which is most often due to hypermethylation of its promoter. Cells displaying MSI accumulate a high number of mutations throughout the genome, especially in short repeat areas called microsatellites. These mutations typically lead to a premature stop codon resulting in a truncated protein product that may inactivate the gene, a mechanism known typical of tumor suppressor genes. The aim of this study is to identify new target genes in MSI CRC. To date, several genes have been proposed as MSI target genes based on high mutation frequency in targeted searches of microsatellites. The availability of new sequencing and bioinformatic technologies has enabled genome-wide investigation of mutations in human cancers. This has led to identification of a vast amount of mutations in cancer, most of which are passenger mutations that do not confer selective growth advantage. The challenge is therefore to distinguish true driver genes from passengers. The high number of mutations in MSI tumors augments this challenge. In this study, 25 sporadic MSI CRCs and their corresponding normal samples have been exome sequenced to identify changes of somatic origin. An analysis and visualization program developed in our group, “RikuRator” (unpublished), will be utilized along with a script to estimate the accurate somatic mutation frequency of each coding mononucleotide repeat. A statistical model based on exome sequencing data will be developed that takes into account the background mutation frequency of microsatellites depending on the repeat length and nucleotide context. The most mutated repeats will be validated by Sanger sequencing in a set of additional MSI CRCs. Functional studies will be carried out to further investigate the pathogenic effect of the new target genes identified. The systematic screening of mutations in microsatellites will improve our understanding on the mutation profile typical of these tumors. With our approach that takes into account the background mutation frequency of microsatellites we aim to identify true driver genes in MSI CRC. A detailed understanding of the molecular background of this tumor type is important for the development of more efficient screening methods and personalized treatments to improve the prognosis of patients diagnosed with MSI type CRC. Citation Format: Johanna Kondelin, Esa Pitkänen, Alexandra E. Gylfe, Kimmo Palin, Heikki Ristolainen, Riku Katainen, Eevi Kaasinen, Minna Taipale, Jussi Taipale, Laura Renkonen-Sinisalo, Heikki Järvinen, Jan Böhm, Jukka-Pekka Mecklin, Pia Vahteristo, Sari Tuupanen, Lauri A. Aaltonen. Identification of new target genes in microsatellite unstable colorectal cancer by exome sequencing. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 2401. doi:10.1158/1538-7445.AM2014-2401

New insights into the incorporation of natural suppressor tRNAs at stop codons in Saccharomyces cerevisiae.

Stop codon readthrough may be promoted by the nucleotide environment or drugs. In such cases, ribosomes incorporate a natural suppressor tRNA at the stop codon, leading to the continuation of translation in the same reading frame until the next stop codon and resulting in the expression of a protein with a new potential function. However, the identity of the natural suppressor tRNAs involved in stop codon readthrough remains unclear, precluding identification of the amino acids incorporated at the stop position. We established an in vivo reporter system for identifying the amino acids incorporated at the stop codon, by mass spectrometry in the yeast Saccharomyces cerevisiae. We found that glutamine, tyrosine and lysine were inserted at UAA and UAG codons, whereas tryptophan, cysteine and arginine were inserted at UGA codon. The 5′ nucleotide context of the stop codon had no impact on the identity or proportion of amino acids incorporated by readthrough. We also found that two different glutamine tRNAGln were used to insert glutamine at UAA and UAG codons. This work constitutes the first systematic analysis of the amino acids incorporated at stop codons, providing important new insights into the decoding rules used by the ribosome to read the genetic code.

Genome evolution by matrix algorithms: cellular automata approach to population genetics.

Mammalian genomes are replete with millions of polymorphic sites, among which those genetic variants that are colocated on the same chromosome and exist close to one another form blocks of closely linked mutations known as haplotypes. The linkage within haplotypes is constantly disrupted due to meiotic recombination events. Whole ensembles of such numerous haplotypes are subjected to evolutionary pressure, where mutations influence each other and should be considered as a whole entity—a gigantic matrix, unique for each individual genome. This idea was implemented into a computational approach, named Genome Evolution by Matrix Algorithms (GEMA) to model genomic changes taking into account all mutations in a population. GEMA has been tested for modeling of entire human chromosomes. The program can precisely mimic real biological processes that have influence on genome evolution such as: 1) Authentic arrangements of genes and functional genomic elements, 2) frequencies of various types of mutations in different nucleotide contexts, and 3) nonrandom distribution of meiotic recombination events along chromosomes. Computer modeling with GEMA has demonstrated that the number of meiotic recombination events per gamete is among the most crucial factors influencing population fitness. In humans, these recombinations create a gamete genome consisting on an average of 48 pieces of corresponding parental chromosomes. Such highly mosaic gamete structure allows preserving fitness of population under the intense influx of novel mutations (40 per individual) even when the number of mutations with deleterious effects is up to ten times more abundant than those with beneficial effects.

Comparative Analysis of Context-Dependent Mutagenesis Using Human and Mouse Models

Substitution rates strongly depend on their nucleotide context. One of the most studied examples is the excess of C > T mutations in the CG context in various groups of organisms, including vertebrates. Studies on the molecular mechanisms underlying this mutation regularity have provided insights into evolution, mutagenesis, and cancer development. Recently several other hypermutable motifs were identified in the human genome. There is an increased frequency of T > C mutations in the second position of the words ATTG and ATAG and an increased frequency of A > C mutations in the first position of the word ACAA. For a better understanding of evolution, it is of interest whether these mutation regularities are human specific or present in other vertebrates, as their presence might affect the validity of currently used substitution models and molecular clocks. A comprehensive analysis of mutagenesis in 4 bp mutation contexts requires a vast amount of mutation data. Such data may be derived from the comparisons of individual genomes or from single nucleotide polymorphism (SNP) databases. Using this approach, we performed a systematical comparison of mutation regularities within 2–4 bp contexts in Mus musculus and Homo sapiens and uncovered that even closely related organisms may have notable differences in context-dependent mutation regularities.

New Universal Rules of Eukaryotic Translation Initiation Fidelity

The accepted model of eukaryotic translation initiation begins with the scanning of the transcript by the pre-initiation complex from the 5′end until an ATG codon with a specific nucleotide (nt) context surrounding it is recognized (Kozak rule). According to this model, ATG codons upstream to the beginning of the ORF should affect translation. We perform for the first time, a genome-wide statistical analysis, uncovering a new, more comprehensive and quantitative, set of initiation rules for improving the cost of translation and its efficiency. Analyzing dozens of eukaryotic genomes, we find that in all frames there is a universal trend of selection for low numbers of ATG codons; specifically, 16–27 codons upstream, but also 5–11 codons downstream of the START ATG, include less ATG codons than expected. We further suggest that there is selection for anti optimal ATG contexts in the vicinity of the START ATG. Thus, the efficiency and fidelity of translation initiation is encoded in the 5′UTR as required by the scanning model, but also at the beginning of the ORF.The observed nt patterns suggest that in all the analyzed organisms the pre-initiation complex often misses the START ATG of the ORF, and may start translation from an alternative initiation start-site. Thus, to prevent the translation of undesired proteins, there is selection for nucleotide sequences with low affinity to the pre-initiation complex near the beginning of the ORF. With the new suggested rules we were able to obtain a twice higher correlation with ribosomal density and protein levels in comparison to the Kozak rule alone (e.g. for protein levels r = 0.7 vs. r = 0.31; p<10−12).

Abstract 1112: Whole exome and whole genome sequence analysis of lung adenocarcinoma.

Abstract Lung adenocarcinoma is the leading cause of cancer-related death worldwide. Recent molecular characterization of this disease through large-scale sequencing efforts has identified distinct subsets driven by mutant oncogenes or kinase fusion proteins, many of which can be inhibited with targeted therapies. Despite these advances, almost half of all lung cancers still lack an identifiable driver. Here, we describe the genomic profiling of 230 normal-paired lung adenocarcinoma samples included as part of The Cancer Genome Atlas (TCGA) effort. All samples were subjected to whole exome analysis, copy number profiling and a subset were also subjected to whole genome sequencing. Mutation calling was performed with the MuTect algorithm. To identify significantly mutated genes, we used the MutSig CV algorithm, a statistically rigorous analysis that takes into account nucleotide context, gene-expression, replication time, and somatic background mutation rate. Mutation rate in lung adenocarcinoma was quite high with an average of 242 mutations/tumor observed (median: 161, range: 11-1328). In total, we identified mutations in over 13,500 genes of which 10 genes reached statistical significance (q&amp;lt;0.1). One significant gene was excluded from further analyses as it was not expressed in RNA-seq data. In addition to mutant genes with established roles in lung adenocarcinoma (e.g. TP53, KRAS, STK11, EGFR, RB1, KEAP1, and BRAF), we also identified other statistically significant mutant genes whose role in lung tumorigenesis is presently unclear. These included mutations in the RNA-binding protein RBM10, and the integrin protein ITGAL. Although statistically insignificant by a small degree, we also identified mutations in the splicing factor U2AF1, and the SWI/SNF complex proteins SMARCA4 and ARID1A. We are currently analyzing whole genome sequences to confirm these events, and identify known and novel fusion events that may be contributing to tumorigenesis. In conclusion, we have analyzed the exomes of 230 lung adenocarcinoma samples and identified known and unknown mutations in this disease. Ultimately, these data will be integrated with ongoing expression, methylation, pathway, miRNA, and proteomic analyses. At its conclusion, this effort will represent the most comprehensive profiling of lung adenocarcinoma samples to date, and will provide a detailed integrative picture of the molecular mechanisms contributing to this disease. Citation Format: Juliann Chmielecki, Mara Rosenberg, Marcin Imielinski, Bryan Hernandez, Michael Lawrence, Andrey Sivachenko, Kristian Cibulskis, Douglas Voet, Carrie Sougnez, Stacey Gabriel, Gad Getz, Matthew Meyerson, The Cancer Genome Atlas Research Network. Whole exome and whole genome sequence analysis of lung adenocarcinoma. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1112. doi:10.1158/1538-7445.AM2013-1112

The roles of the monomer length and nucleotide context of plant tandem repeats in nucleosome positioning

Similar to regularly spaced nucleosomes in chromatin, long tandem DNA arrays are composed of regularly alternating monomers that have almost identical primary DNA structures. Such a similarity in the structural organization makes these arrays especially interesting for studying the role of intrinsic DNA preferences in nucleosome positioning. We have studied the nucleosome formation potential of DNA tandem repeat families with different monomer lengths (ML). In total, 165 plant tandem repeat families from the PlantSat database (http://w3lamc.umbr.cas.cz/PlantSat/) were divided into two classes based on the number of nucleosome repeats in one DNA monomer. For predicting nucleosome formation potential, we developed the Phase method, which combines the advantages of multiple bioinformatics models. The Phase method was able to distinguish interfamily differences and intrafamily monomer variation and identify the influence of nucleotide context on nucleosome formation potential. Three main types of nucleosome arrangement in DNA tandem repeat arrays – regular, partially regular (partial), and flexible – were distinguished among a great variety of Phase profiles. The regular type, in which all nucleosomes of the monomer array are positioned in a context-dependent manner, is the most representative type of the class 1 families, with ML equal to or a multiple of the nucleosome repeat length (NRL). In the partially regular type, nucleotide context influences the positioning of only a subset of nucleosomes. The influence of the nucleotide context on nucleosome positioning has the least effect in the flexible type, which contains the greatest number of families (65). The majority of these families belong to class 2 and have nonmultiple ML to NRL ratios.

114 Kinetic and thermodynamic basis for damaged bases excision by 8-oxoguanine-DNA glycosylases

DNA glycosylases play the opening act in a highly conserved process for excision of damaged bases from DNA called the base excision repair pathway. DNA glycosylases attend to a wide variety of lesions arising from both endogenous and exogenous factors. The types of damage include alkylation, oxidation, and hydrolysis. A major DNA oxidation product is 8-oxoguanine (8-oxoG), a base with a high mutagenic potential. In bacteria, this lesion is repaired by formamidopyrimidine-DNA glycosylase (Fpg), while in the case of humans this function belongs to 8-oxoG-DNA glycosylase (OGG1). We have attempted a comprehensive characterization of 8-oxoG recognition by DNA glycosylases. First, we have obtained thermodynamic parameters for melting of DNA duplexes containing 8-oxoG in all possible nucleotide contexts. The energy of stacking interactions of 8-oxoG was in strict dependence on 8-oxoG nucleotide environment, which may affect the recognition of damage and the efficiency of eversion of 8-oxoG from DNA helix by glycosylases. Next, we established how the flexibility of DNA context affects damage recognition by these enzymes (Kirpota et al., 2011). Then, we have found that DNA containing 8-oxoG next to a single-strand break provides a good substrate for Fpg, as soon as all structural phosphate residues are maintained. Using site-directed mutagenesis, we have addressed the functions of many previously unstudied amino acid residuess that were predicted to be important for Fpg activity by molecular dynamics simulation and phylogenetic analysis. Of note, many substitutions abolished the excision of 8-oxoG, but did not affect the cleavage efficiency of abasic substrates. Finally, we investigated the contribution of separated structural domains of Fpg to specific enzyme-substrate interaction. Surprisingly, despite the absence of the catalytic domain, C-terminal domain of Fpg possessed a low- residual ability to recognize and cleave abasic substrates. Our study sheds light on mechanism details of Fpg and OGG1 activity, with the ultimate goal of understanding how binding energy can be spent by these enzymes for catalysis.

Translation efficiency in yeasts correlates with nucleosome formation in promoters

Elongation efficiency index (EEI) was suggested earlier to estimate gene expression efficiency by nucleotide context of coding sequence in unicellular organisms. We have analyzed association between EEI and nucleosome formation potential (NFP) in 5’ regulatory regions upstream translation initiation site (TIS) from two yeast species. Theoretical estimations of NFP based on DNA sequence were obtained by Recon method. Experimental estimation of nucleosome occupancy was obtained by high-throughput sequencing data of nucleosomal DNA in Saccharomyces cerevisiae . For the sample of all genes correlation coefficient was calculated between two vectors: vector of NFP values for fixed position relative to TIS and vector of EEI values. Profiles of correlation coefficients of NFP and EEI were counted in (−600; +600) regions relative to TIS for gene sequences extracted from GenBank. We found regions of strong negative dependence between NFP and EEI for all genes as well as for 10% highly expressed genes in Schizosaccharomyces pombe (10% of EEI-highest genes). At the same time, we found positive dependence between NFP and EEI for all genes and for low expressed genes in S. cerevisiae (10% of EEI-lowest genes). The association between NFP and EEI could be explained by evolutionary selection of context characteristics of nucleotide sequences for gene expression optimization.

Hidden coding potential of eukaryotic genomes: nonAUG started ORFs

It is widely considered that the vast majority of eukaryotic mRNAs contain only one open reading frame (ORF) and encode single protein. However, eukaryotic ribosomes can initiate translation at alternative start codons due to leaky scanning or reinitiation mechanisms that provides an opportunity to synthesize several protein products. Recent investigations also demonstrated that alternative translation from nonAUG start codons and AUG codons in a weak nucleotide context could make an important contribution to eukaryotic proteomes. However, accurate prediction of alternative start codons demands detailed investigation of mRNA features influencing their recognition by eukaryotic ribosomes. In this work, we present the results of computational analysis of characteristics of yeast and mammalian mRNAs potentially involved in the recognition of nonAUG start codons. It was found that sequence features of nonAUG started Saccharomyces cerevisiae upstream ORFs (uORFs) were adjusted for efficient translation and these uORFs could frequently encode functional polypeptides. In particular, our initial studies revealed that predicted tertiary structures downstream of nonAUG start sites in mammalian mRNAs were energetically more stable than those predicted for AUG start sites with strong Kozak context. We hypothesize that presence of such stable tertiary structure downstream of nonAUG start sites could be an important factor for the ribosome to recognize noncanonical start codons.

HIV-1 genome is often defective in PBMCs and rectal tissues after long-term HAART as a result of APOBEC3 editing and correlates with the size of reservoirs

Precise characterization of viruses present in reservoirs in long-term pretreated patients will be a major issue to consider in the context of viral eradication. We assessed the frequency of defective viruses present in cellular reservoirs. Peripheral blood mononuclear cells (PBMCs) and rectal biopsy samples were compared between five patients on successful long-term highly active antiretroviral therapy (HAART) (>7 years without blips) and five untreated patients. Molecular cloning and sequencing of the reverse transcriptase region were used to detect the presence of and quantify in-frame stop codons in HIV quasi-species. The relationship between the size of the reservoir and the frequency of defective genomes was assessed. Defective genomes were systematically detected in all patients on long-term HAART in both compartments (PBMCs and rectal tissues), with a higher level of defective genomes per sample compared with PBMCs of untreated patients. A high level of defective genomes was correlated with a small size of HIV proviral DNA. Regarding the nucleotide context, guanine (G) to adenine (A) substitution at tryptophan positions was responsible for the appearance of 89% of all in-frame stop codons in the context of G-to-A hypermutation, likely reflecting APOBEC3 footprints on the viral genome. We propose a scenario whereby defective genomes accumulate during HAART treatment, eventually reaching a viral extinction threshold. In the context of viral eradication, measurement of the relative amounts of defective and non-defective viruses (by molecular cloning and ultradeep sequencing) should be used as a new criterion for eradicating HIV.

Complexities of 5'splice site definition: Implications in clinical analyses

In higher eukaryotes, the 5' splice site (5'ss) is initially recognized through an RNA-RNA interaction by U1 small nuclear ribonucleoprotein (U1 snRNP). This event represents one of the key steps in initial spliceosomal assembly and many disease-associated mutations in humans often disrupt this process. Beside base pair complementarity, 5'ss recognition can also be modified by additional factors such as RNA secondary structures or the specific binding of other nuclear proteins. In this work, we have focused on investigating a few examples of changes detected within the 5'ss in patients, that would not be immediately considered "disease causing mutations". We show that the splicing outcome of very similar mutations can be very different due to variations in trans-acting factor(s) interactions and specific context influences. Using several NF1 donor sites and SELEX approaches as experimental models, we have examined the binding properties of particular sequence motifs such as GGGU found in donor sites, and how the sequence context can change their interaction with hnRNPs such as H/F and A1/A2. Our results clearly show that even minor differences in local nucleotide context can differentially affect the binding ability of these factors to the GGGU core. Finally, using a previously identified mutation in KCNH2 that resulted in intron retention we show how very similar 5'ss mutations found in patients can have a very different splicing outcome due to the neighbouring sequence context, thus highlighting the general need to approach splicing problems with suitable experimental approaches.

ANALYSIS OF CONTEXT-DEPENDENT ERRORS FOR ILLUMINA SEQUENCING

The new generation of short-read sequencing technologies requires reliable measures of data quality. Such measures are especially important for variant calling. However, in the particular case of SNP calling, a great number of false-positive SNPs may be obtained. One needs to distinguish putative SNPs from sequencing or other errors. We found that not only the probability of sequencing errors (i.e. the quality value) is important to distinguish an FP-SNP but also the conditional probability of "correcting" this error (the "second best call" probability, conditional on that of the first call). Surprisingly, around 80% of mismatches can be "corrected" with this second call. Another way to reduce the rate of FP-SNPs is to retrieve DNA motifs that seem to be prone to sequencing errors, and to attach a corresponding conditional quality value to these motifs. We have developed several measures to distinguish between sequence errors and candidate SNPs, based on a base call's nucleotide context and its mismatch type. In addition, we suggested a simple method to correct the majority of mismatches, based on conditional probability of their "second" best intensity call. We attach a corresponding second call confidence (quality value) of being corrected to each mismatch.

Statistical Analysis of Readthrough Levels for Nonsense Mutations in Mammalian Cells Reveals a Major Determinant of Response to Gentamicin

The efficiency of translation termination depends on the nature of the stop codon and the surrounding nucleotides. Some molecules, such as aminoglycoside antibiotics (gentamicin), decrease termination efficiency and are currently being evaluated for diseases caused by premature termination codons. However, the readthrough response to treatment is highly variable and little is known about the rules governing readthrough level and response to aminoglycosides. In this study, we carried out in-depth statistical analysis on a very large set of nonsense mutations to decipher the elements of nucleotide context responsible for modulating readthrough levels and gentamicin response. We quantified readthrough for 66 sequences containing a stop codon, in the presence and absence of gentamicin, in cultured mammalian cells. We demonstrated that the efficiency of readthrough after treatment is determined by the complex interplay between the stop codon and a larger sequence context. There was a strong positive correlation between basal and induced readthrough levels, and a weak negative correlation between basal readthrough level and gentamicin response (i.e. the factor of increase from basal to induced readthrough levels). The identity of the stop codon did not affect the response to gentamicin treatment. In agreement with a previous report, we confirm that the presence of a cytosine in +4 position promotes higher basal and gentamicin-induced readthrough than other nucleotides. We highlight for the first time that the presence of a uracil residue immediately upstream from the stop codon is a major determinant of the response to gentamicin. Moreover, this effect was mediated by the nucleotide itself, rather than by the amino-acid or tRNA corresponding to the −1 codon. Finally, we point out that a uracil at this position associated with a cytosine at +4 results in an optimal gentamicin-induced readthrough, which is the therapeutically relevant variable.

Heterogeneity of the Transition/Transversion Ratio in Drosophila and Hominidae Genomes

Mutation rate varies between sites in the genome. Part of this variation can be explained by well-recognized short nucleotide contexts, but a large component of this variation remains cryptic. We used data on interspecies divergence and intraspecies polymorphism in Drosophila and Hominidae to analyze variation of the average rate of the 12 possible kinds of single-nucleotide mutations and in the transition/transversion ratio κ at single-nucleotide resolution. Both the average mutation rate and κ vary by a factor of ~3 between nucleotide sites. The characteristic scale of variation in κ is up to at least ~30 nucleotides in Drosophila and ~5 nucleotides in Hominidae. Genome segments with locally elevated mutation rates possess lower values of κ; however, a substantial fraction of variation in κ cannot be directly explained by the local mutation rates.

Inferring Divergence of Context-Dependent Substitution Rates in Drosophila Genomes with Applications to Comparative Genomics

Nucleotide substitution is a major evolutionary driving force that can incrementally and stochastically give rise to broad divergence patterns among species. The substitution process at each genomic position is frequently modeled independently of the other positions, although complex interactions between nearby bases are known to significantly affect mutation rates. Here, we study the evolution of 12 fly genomes using new algorithms for accurate inference of parameter-rich substitution models. By comparing models between lineages, we reveal the evolutionary histories of substitution rates at different flanking nucleotide contexts. We demonstrate these driving forces of molecular evolution to be constantly changing, suggesting that neutral drift of mutation rates is an important factor in the evolution of genomes and their sequence composition. This observation is used to develop a scalable approach for parameter-rich comparative genomics. By screening short DNA sequences, we demonstrate how homeoboxes and other transcription factor binding motifs are highly conserved based on our parameter-rich models but not according to standard conservation assays. With the increasing availability of genome sequences, rich substitution models become an attractive and practical approach for evolutionary analysis in general and comparative genomics in particular.

Patterns of Repeat-Induced Point Mutation in Transposable Elements of Basidiomycete Fungi

Transposable elements (TEs) are ubiquitous genomic parasites that have prompted the evolution of genome defense systems that restrict their activity. Repeat-induced point mutation (RIP) is a homology-dependent genome defense that introduces C-to-T transition mutations in duplicated DNA sequences and is thought to control the proliferation of selfish repetitive DNA. Here, we determine the taxonomic distribution of hypermutation patterns indicative of RIP among basidiomycetes. We quantify C-to-T transition mutations in particular di- and trinucleotide target sites for TE-like sequences from nine fungal genomes. We find evidence of RIP-like patterns of hypermutation at TpCpG trinucleotide sites in repetitive sequences from all species of the Pucciniomycotina subphylum of the Basidiomycota, Microbotryum lychnidis-dioicae, Puccinia graminis, Melampsora laricis-populina, and Rhodotorula graminis. In contrast, we do not find evidence for RIP-like hypermutation in four species of the Agaricomycotina and Ustilaginomycotina subphyla of the Basidiomycota. Our results suggest that a RIP-like process and the specific nucleotide context for mutations are conserved within the Pucciniomycotina subphylum. These findings imply that coevolutionary interactions between TEs and a hypermutating genome defense are stable over long evolutionary timescales.

Stringency of start codon selection modulates autoregulation of translation initiation factor eIF5

An AUG in an optimal nucleotide context is the preferred translation initiation site in eukaryotic cells. Interactions among translation initiation factors, including eIF1 and eIF5, govern start codon selection. Experiments described here showed that high intracellular eIF5 levels reduced the stringency of start codon selection in human cells. In contrast, high intracellular eIF1 levels increased stringency. High levels of eIF5 induced translation of inhibitory upstream open reading frames (uORFs) in eIF5 mRNA that initiate with AUG codons in conserved poor contexts. This resulted in reduced translation from the downstream eIF5 start codon, indicating that eIF5 autoregulates its own synthesis. As with eIF1, which is also autoregulated through translation initiation, features contributing to eIF5 autoregulation show deep evolutionary conservation. The results obtained provide the basis for a model in which auto- and cross-regulation of eIF5 and eIF1 translation establish a regulatory feedback loop that would stabilize the stringency of start codon selection.