Eukaryotic noncoding DNA sequences provide genes with an additional protection against chemical mutagens

A quantitative model was developed that identified a new function of noncoding sequences in the eukaryotic genome, namely, the protection of coding sequences against chemical (mainly endogenous) mutagens. It was shown that, under common ecological conditions, the number of nucleotides damaged by endogenous and exogenous reactive species in coding sequences of the genome is inversely proportional to the size of their noncoding parts. Noncoding sequences can differently protect single genetic loci from chemical modification by the formation of specific spatial structures of the protected loci in the interphase nuclei. The significant differences in genome sizes between species (C-value paradox) can be explained by different contributions of noncoding sequences to the total effect of genome protection from endogenous chemical mutagens.     

Genomic noncoding sequences and the size of eukaryotic cell nucleus as important factors of gene protection from chemical mutagens

An improved quantitative model describing a protective function of eukaryotic genomic noncoding sequences was developed. In this new model, two factors affecting gene protection from chemical mutagens are considered: (1) the ratio of the total lengths of coding and noncoding genomic sequences and (2) the volume of the cell nucleus. An increase in the noncoding DNA in the genome reduces the number of mutagen-damaged nucleotides in the coding region, whereas an increase in the volume of the nucleus decreases the flow of mutagens per unit of nuclear volume that attacks its surface.     

Genomic noncoding sequences and the size of eukaryotic cell nucleus as important factors of gene protection from chemical mutagens

An improved quantitative model describing a protective function of eukaryotic genomic noncoding sequences was developed. In this new model, two factors affecting gene protection from chemical mutagensare considered: (1) the ratio of the total lengths of coding and noncoding genomic sequences and (2) the volume of the cell nucleus. An increase in the noncoding DNA in the genome reduces the number of mutagen-damaged nucleotides in the coding region, whereas an increase in the volume of the nucleus decreases the flow of mutagens per unit of nuclear volume that attacks its surface.     

The problem of the eukaryotic genome size

The current state of knowledge concerning the unsolved problem of the huge interspecific eukaryotic genome size variations not correlating with the species phenotypic complexity (C-value enigma also known as C-value paradox) is reviewed. Characteristic features of eukaryotic genome structure and molecular mechanisms that are the basis of genome size changes are examined in connection with the C-value enigma. It is emphasized that endogenous mutagens, including reactive oxygen species, create a constant nuclear environment where any genome evolves. An original quantitative model and general conception are proposed to explain the C-value enigma. In accordance with the theory, the noncoding sequences of the eukaryotic genome provide genes with global and differential protection against chemical mutagens and (in addition to the anti-mutagenesis and DNA repair systems) form a new, third system that protects eukaryotic genetic information. The joint action of these systems controls the spontaneous mutation rate in coding sequences of the eukaryotic genome. It is hypothesized that the genome size is inversely proportional to functional efficiency of the anti-mutagenesis and/or DNA repair systems in a particular biological species. In this connection, a model of eukaryotic genome evolution is proposed.     

Construction of viable deletion and insertion mutants of the Sabin strain of type 1 poliovirus: function of the 5'' noncoding sequence in viral replication.

A number of deletion and insertion sequences were introduced into the 5' noncoding sequence (742 nucleotides long) of the genome of the Sabin strain of type 1 poliovirus by using an infectious cDNA clone of the virus strain. The genomes of all three poliovirus serotypes contained highly homologous sequences (nucleotide positions 509 to 639) as well as highly variable sequences (positions 640 to 742) in the 5' noncoding region. The viability of mutant viruses was tested by transfecting mutant cDNA clones into African green monkey kidney cells and then estimating the plaque sizes displayed on the cells. The results suggested that the highly variable sequence next to the VP4 coding region did not play an important role, at least in the in vitro culture system used, that the loci of highly conserved nucleotide sequences were not always expected to be the genome regions essential for viral replication, that the sequence between positions 564 and 599 carried genetic information to maintain the efficiency of certain steps in viral replication, and that the sequence between positions 551 to 563 might play an essential role in viral replication. Four-base deletion or insertion mutations were introduced into relatively variable sequences in the genome region upstream of position 509. The results suggest that variable sequences do not always indicate that the corresponding genome regions are less important. Apparent revertants (large-plaque variants) were easily generated from one of the viable mutants with the small-plaque phenotype. The determination of nucleotide sequences of the revertant genomes revealed the second mutation site. The results suggested that the different loci at around positions 200 and 500 might specifically interact with each other. This interaction may result in the formation of a functional structure that influences the efficiency of certain steps in the viral replication.     

The contribution of slippage-like processes to genome evolution

Simple sequences present in long (>30 kb) sequences representative of the single-copy genome of five species (Homo sapiens, Caenorhabditis elegans Saccharomyces cerevisiae, E. coli, and Mycobacterium leprae) have been analyzed. A close relationship was observed between genome size and the overall level of sequence repetition. This suggested that the incorporation of simple sequences had accompanied increases of genome size during evolution. Densities of simple sequence motifs were higher in noncoding regions than in coding regions in eukaryotes but not in eubacteria. All five genomes showed very biased frequency distributions of simple sequence motifs in all species, particularly in eukaryotes where AAA and TTT predominated. Interspecific comparisons showed that noncoding sequences in eukaryotes showed highly significantly similar frequency distributions of simple sequence motifs but this was not true of coding sequences. ANOVA of the frequency distributions of simple sequence motifs indicated strong contributions from motif base composition and repeat unit length, but much of the variation remained unexplained by these parameters. The sequence composition of simple sequences therefore appears to reflect both underlying sequence biases in slippage-like processes and the action of selection. Frequency distributions of simple sequence motifs in coding sequences correlated weakly or not at all with those in noncoding sequences. Selection on coding sequences to eliminate undesirable sequences may therefore have been strong, particularly in the human lineage.     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Hypervariable noncoding sequences in Saccharomyces cerevisiae

Compared to protein-coding sequences, the evolution of noncoding sequences and the selective constraints placed on these sequences is not well characterized. To compare the evolution of coding and noncoding sequences, we have conducted a survey for DNA polymorphism at five randomly chosen loci among a diverse collection of 81 strains of Saccharomyces cerevisiae. Average rates of both polymorphism and divergence are 40% lower at noncoding sites and 90% lower at nonsynonymous sites in comparison to synonymous sites. Although noncoding and coding sequences show substantial variability in ratios of polymorphism to divergence, two of the loci, MLS1 and PDR10, show a higher rate of polymorphism at noncoding compared to synonymous sites. The high rate of polymorphism is not accompanied by a high rate of divergence and is limited to a few small regions. These hypervariable regions include sites with three segregating bases at a single site and adjacent polymorphic sites. We show that this clustering of polymorphic sites is significantly greater than one would expect on the basis of the spacing between polymorphic fourfold degenerate sites. Although hypervariable noncoding sequences could result from selection on regulatory mutations, they could also result from transient mutational hotspots.     

A catalog of neutral and deleterious polymorphism in yeast

The abundance and identity of functional variation segregating in natural populations is paramount to dissecting the molecular basis of quantitative traits as well as human genetic diseases. Genome sequencing of multiple organisms of the same species provides an efficient means of cataloging rearrangements, insertion, or deletion polymorphisms (InDels) and single-nucleotide polymorphisms (SNPs). While inbreeding depression and heterosis imply that a substantial amount of polymorphism is deleterious, distinguishing deleterious from neutral polymorphism remains a significant challenge. To identify deleterious and neutral DNA sequence variation within Saccharomyces cerevisiae, we sequenced the genome of a vineyard and oak tree strain and compared them to a reference genome. Among these three strains, 6% of the genome is variable, mostly attributable to variation in genome content that results from large InDels. Out of the 88,000 polymorphisms identified, 93% are SNPs and a small but significant fraction can be attributed to recent interspecific introgression and ectopic gene conversion. In comparison to the reference genome, there is substantial evidence for functional variation in gene content and structure that results from large InDels, frame-shifts, and polymorphic start and stop codons. Comparison of polymorphism to divergence reveals scant evidence for positive selection but an abundance of evidence for deleterious SNPs. We estimate that 12% of coding and 7% of noncoding SNPs are deleterious. Based on divergence among 11 yeast species, we identified 1,666 nonsynonymous SNPs that disrupt conserved amino acids and 1,863 noncoding SNPs that disrupt conserved noncoding motifs. The deleterious coding SNPs include those known to affect quantitative traits, and a subset of the deleterious noncoding SNPs occurs in the promoters of genes that show allele-specific expression, implying that some cis-regulatory SNPs are deleterious. Our results show that the genome sequences of both closely and distantly related species provide a means of identifying deleterious polymorphisms that disrupt functionally conserved coding and noncoding sequences.     

Base composition of Coffea AFLP sequences and their conservation within the genus

Amplified fragment length polymorphism (AFLP) is often used for genetic mapping and diversity analysis, but very little information is currently available on their sequence characteristics. Species-specific sequences were analyzed from a single Coffea genome (Coffea pseudozanguebariae) associated with clustered or nonclustered AFLP loci of known genetic position. Compared with the expressed sequence tag (EST) sequence composition, their AT content exhibited a bimodal distribution with AT-poor sequences corresponding mainly to putative coding sequences. AT-rich sequences, apart from the EST distribution, were usually clustered on the genetic map and might correspond to noncoding sequences. Conversion of these AFLP markers into sequence-characterized amplified region (SCAR) anchor markers allowed us to assess sequence conservation within Coffea species with respect to species relatedness.     

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.     

Genomic heterogeneity in the density of noncoding single-nucleotide and microsatellite polymorphisms in Plasmodium falciparum

The density and distribution of single-nucleotide polymorphisms (SNPs) across the genome has important implications for linkage disequilibrium mapping and association studies, and the level of simple-sequence microsatellite polymorphisms has important implications for the use of oligonucleotide hybridization methods to genotype SNPs. To assess the density of these types of polymorphisms in P. falciparum, we sampled introns and noncoding DNA upstream and downstream of coding regions among a variety of geographically diverse parasites. Across 36,229 base pairs of noncoding sequence representing 41 genetic loci, a total of 307 polymorphisms including 248 polymorphic microsatellites and 39 SNPs were identified. We found a significant excess of microsatellite polymorphisms having a repeat unit length of one or two, compared to those with longer repeat lengths, as well as a nonrandom distribution of SNP polymorphisms. Almost half of the SNPs localized to only three of the 41 genetic loci sampled. Furthermore, we find significant differences in the frequency of polymorphisms across the two chromosomes (2 and 3) examined most extensively, with an excess of SNPs and a surplus of polymorphic microsatellites on chromosome 3 as compared to chromosome 2 (P=0.0001). Furthermore, at some individual genetic loci we also find a nonrandom distribution of polymorphisms between coding and flanking noncoding sequences, where completely monomorphic regions may flank highly polymorphic genes. These data, combined with our previous findings of nonrandom distribution of SNPs across chromosome 2, suggest that the Plasmodium falciparum genome may be a mosaic with regard to genetic diversity, containing chromosomal regions that are highly polymorphic interspersed with regions that are much less polymorphic.     

Correlation approach to identify coding regions in DNA sequences.

Recently, it was observed that noncoding regions of DNA sequences possess long-range power-law correlations, whereas coding regions typically display only short-range correlations. We develop an algorithm based on this finding that enables investigators to perform a statistical analysis on long DNA sequences to locate possible coding regions. The algorithm is particularly successful in predicting the location of lengthy coding regions. For example, for the complete genome of yeast chromosome III (315,344 nucleotides), at least 82% of the predictions correspond to putative coding regions; the algorithm correctly identified all coding regions larger than 3000 nucleotides, 92% of coding regions between 2000 and 3000 nucleotides long, and 79% of coding regions between 1000 and 2000 nucleotides. The predictive ability of this new algorithm supports the claim that there is a fundamental difference in the correlation property between coding and noncoding sequences. This algorithm, which is not species-dependent, can be implemented with other techniques for rapidly and accurately locating relatively long coding regions in genomic sequences.     

Evolution of programmed DNA rearrangements in a scrambled gene

Gene unscrambling in spirotrichous ciliates involves massive genome-wide DNA deletion and rearrangement events during development. During each sexual cycle, the somatic nucleus (macronucleus) regenerates from the germ line nucleus (micronucleus). Development of the polyploid somatic genome requires programmed DNA deletion of micronuclear-limited intragenic noncoding sequences and permutation and amplification of the protein-coding regions. Recent studies suggest that, despite novel insertions of endogenous transposon or foreign DNA into the germ line genome, ciliates possess a whole-genome surveillance system that guides the recapitulation of a functional somatic genome. This renders the germ line genome an extremely dynamic structure over evolutionary time. Here we describe the germ line and somatic architectures of the gene encoding alpha-telomere-binding protein in three early-diverging species (Holosticha sp., Uroleptus sp., and Paraurostyla weissei) and trace the natural history of DNA rearrangements in this gene in six species, including three previously studied oxytrichids. Comparisons of homologous coding regions between earlier and later diverging species provide evidence for fusion of scrambled germ line fragments as small as 24 bp during evolution, as well as simultaneous fragmentation and scrambling of the germ line locus and shifting of the boundaries between coding and noncoding DNA, leading to distinct gene architectures in each species. We infer an evolutionary recombination pathway that passes through identified intermediate species and gives rise to the observed patterns in all known species, capitalizing on their unique DNA rearrangement machinery and germ line flexibility.