Molecular evolution of coding and non-coding regions in Plasmodium

Recurrence analysis provides a useful tool for the characterisation of oligonucleotide usage along genomic tracts. While coding regions are characterised by a low-recurrence regimen (except in the case of intragenic repeats) introns and intergenic regions exhibit a high density of recurring oligos, and appear to be correlated from the point of view of oligonucleotide preference. By comparing homologous loci in Plasmodium falciparum and P. berghei, it can be seen that introns and intergenic regions, though exhibiting very low sequence similarity, do not drift without constraints, but maintain a consistent use of the same oligos in the two species.     

酵母基因上游序列中潜在的转录正调控位点分析

前期研究表明,高效转录酵母基因内含子在序列长度、寡核苷酸使用、以及位置分布等方面都有着区别于低转录内含子的特征 . 进一步观察发现：上游基因间区域的序列长度与基因转录频率也有与内含子序列相同的现象,转录频率高的上游基因间序列一般都比转录频率低的长 . 对高效转录和低效转录上游基因间序列的寡核苷酸使用频率进行统计比较分析,抽提出高转录基因上游区可能的转录正调控元件 . 与酵母的所有非编码序列比较,这些可能的正调控元件基本上也是过表达的 (over-represented) ,其中多数和实验所得的一些位点特征相吻合 . 这些元件富含 G 、 C ,这与内含子中可能的正调控元件在碱基组成上有一定的互补性 . 从这些特征看,高效转录基因上游的序列结构确实有利于基因的转录 .     

Noncoding DNA, isochores and gene expression: nucleosome formation potential

The nucleosome formation potential of introns, intergenic spacers and exons of human genes is shown here to negatively correlate with among-tissues breadth of gene expression. The nucleosome formation potential is also found to negatively correlate with the GC content of genomic sequences; the slope of regression line is steeper in exons compared with noncoding DNA (introns and intergenic spacers). The correlation with GC content is independent of sequence length; in turn, the nucleosome formation potential of introns and intergenic spacers positively (albeit weakly) correlates with sequence length independently of GC content. These findings help explain the functional significance of the isochores (regions differing in GC content) in the human genome as a result of optimization of genomic structure for epigenetic complexity and support the notion that noncoding DNA is important for orderly chromatin condensation and chromatin-mediated suppression of tissue-specific genes.     

Insertion/deletion and nucleotide polymorphism data reveal constraints in Drosophila melanogaster introns and intergenic regions

Our study of nucleotide sequence and insertion/deletion polymorphism in Drosophila melanogaster noncoding DNA provides evidence for selective pressures in both intergenic regions and introns (of the large size class). Intronic and intergenic sequences show a similar polymorphic deletion bias. Insertions have smaller sizes and higher frequencies than deletions, supporting the hypothesis that insertions are selected to compensate for the loss of DNA caused by deletion bias. Analysis of a simple model of selective constraints suggests that the blocks of functional elements located in intergenic sequences are on average larger than those in introns, while the length distribution of relatively unconstrained sequences interspaced between these blocks is similar in intronic and intergenic regions.     

真核生物DNA非编码区的组分分析

在全基因组水平上,用直方图、混沌表示灰度图、距离差异度和信息熵差异度四种方法,研究了拟南芥、线虫、果蝇的DNA内含子、基因间隔区DNA、外显子三种区域的核苷酸短序列组分及组分复杂度.结果表明：a.不同基因组之间,不管基因数目多少,用4种方法得到的外显子部分其组分复杂度都比较接近,而非编码区部分的组分复杂度却很大.这一点定量地说明了物种之间的复杂程度,主要不体现在编码区部分,而体现在非编码区部分.b.同一基因组中,内含子的核苷酸短序列组分复杂度都是相似的,外显子和intergenic DNA部分的组分复杂度也是相似的.c.内含子和intergenic DNA在转录、剪切、二级结构等方面有很大的不同,但它们在核苷酸短序列组分上的差异却很小,说明内含子和intergenic DNA在转录、剪切、二级结构上的不同并不通过核苷酸短序列组分来进行限制.     

Evolution of noncoding and silent coding sites in the Plasmodium falciparum and Plasmodium reichenowi genomes

We compared levels of sequence divergence between fourfold synonymous coding sites and noncoding sites from the intergenic and intronic regions of the Plasmodium falciparum and Plasmodium reichenowi genomes. We observed significant differences in the level of divergence between these classes of silent sites. Fourfold synonymous coding sites exhibited the highest level of sequence divergence, followed by introns, and then intergenic sequences. This pattern of relative divergence rates has been observed in primate genomes but was unexpected in Plasmodium due to a paucity of variation at silent sites in P. falciparum and the corollary hypothesis that silent sites in this genome may be subject to atypical selective constraints. Exclusion of hypermutable CpG dinucleotides reduces the divergence level of synonymous coding sites to that of intergenic sites but does not diminish the significantly higher divergence level of introns relative to intergenic sites. A greater than expected incidence of CpG dinucleotides in intergenic regions less than 500 bp from genes may indicate selective maintenance of regulatory motifs containing CpGs. Divergence rates of different classes of silent sites in these Plasmodium genomes are determined by a combination of mutational and selective pressures.     

The Correlation Between Recombination Rate and Dinucleotide Bias in Drosophila melanogaster

Revealing how recombination affects genomic sequence is of great significance to our understanding of genome evolution. The present paper focuses on the correlation between recombination rate and dinucleotide bias in Drosophila melanogaster genome. Our results show that the overall dinucleotide bias is positively correlated with recombination rate for genomic sequences including untranslated regions, introns, intergenic regions, and coding sequences. The correlation patterns of individual dinucleotide biases with recombination rate are presented. Possible mechanisms of interaction between recombination and dinucleotide bias are discussed. Our data indicate that there may be a genome-wide universal mechanism acting between recombination rate and dinucleotide bias, which is likely to be neighbor-dependent biased gene conversion.     

Contrasting patterns of sequence divergence and base composition between Drosophila introns and intergenic regions

Two non-coding DNA classes, introns and intergenic regions, of Drosophila melanogaster exhibit contrasting evolutionary patterns. GC content is significantly higher in intergenic regions and affects their degree of nucleotide variability. Divergence is positively correlated with recombination rate in intergenic regions, but not in introns. We argue that these differences are due to different selective constraints rather than mutational or recombinational mechanisms.     

Background selection in single genes may explain patterns of codon bias

Background selection involves the reduction in effective population size caused by the removal of recurrent deleterious mutations from a population. Previous work has examined this process for large genomic regions. Here we focus on the level of a single gene or small group of genes and investigate how the effects of background selection caused by nonsynonymous mutations are influenced by the lengths of coding sequences, the number and length of introns, intergenic distances, neighboring genes, mutation rate, and recombination rate. We generate our predictions from estimates of the distribution of the fitness effects of nonsynonymous mutations, obtained from DNA sequence diversity data in Drosophila. Results for genes in regions with typical frequencies of crossing over in Drosophila melanogaster suggest that background selection may influence the effective population sizes of different regions of the same gene, consistent with observed differences in codon usage bias along genes. It may also help to cause the observed effects of gene length and introns on codon usage. Gene conversion plays a crucial role in determining the sizes of these effects. The model overpredicts the effects of background selection with large groups of nonrecombining genes, because it ignores Hill-Robertson interference among the mutations involved.     

Localizing triplet periodicity in DNA and cDNA sequences

Background  

The protein-coding regions (coding exons) of a DNA sequence exhibit a triplet periodicity (TP) due to fact that coding exons contain a series of three nucleotide codons that encode specific amino acid residues. Such periodicity is usually not observed in introns and intergenic regions. If a DNA sequence is divided into small segments and a Fourier Transform is applied on each segment, a strong peak at frequency 1/3 is typically observed in the Fourier spectrum of coding segments, but not in non-coding regions. This property has been used in identifying the locations of protein-coding genes in unannotated sequence. The method is fast and requires no training. However, the need to compute the Fourier Transform across a segment (window) of arbitrary size affects the accuracy with which one can localize TP boundaries. Here, we report a technique that provides higher-resolution identification of these boundaries, and use the technique to explore the biological correlates of TP regions in the genome of the model organism C. elegans.     

Sequence studies on the soybean chloroplast 16S–23S rDNA spacer region

The sequence of the ribosomal spacer region of soybean chloroplast DNA including the 3 end of the 16S rRNA gene, the tRNA^Ala and tRNA^Ile genes (but not their introns), the three intergenic regions and the 5 end of the 23S rRNA gene, has been determined. This sequence has been compared to corresponding regions of other angiosperm chloroplast DNAs. Secondary structure models are proposed for the entirety of the intergenic regions a, b and c and for the flanking rRNA regions. A model for a common secondary structure of the ribosomal spacer intergenic regions from chloroplasts of higher plants is proposed, which is supported by comparative evidence.     

The D4 set: primers that target highly variable intron loops in plant chloroplast genomes

Chloroplast group II introns offer high-quality, rapidly evolving single-copy loci for comparative sequence analysis. These introns feature diagnostic secondary structures with loops that are among the least evolutionarily constrained sequence in plastomes. We exploited these structures to develop universal primers that amplify and sequence the large Domain IV (D4) loop in several angiosperm introns. With a single sequence read, we recover 300-600 nucleotides of highly variable sequence across angiosperms, with rates of change that are equal to or higher than many of the best known intergenic spacers in plant chloroplast genomes.     

Organization and variation of angiosperm mitochondrial genome

The mitochondrial genomes of angiosperms are the largest mitochondrial genomes so far reported and are highly variable in size among plant species. The comparative analysis of the angiosperm mitochondrial genomes at the nucleotide level has now become feasible for addressing long-standing questions, owing to the publication of five dicot and three monocot genomes. Whereas the identified genes and introns are rather well conserved, intergenic regions are highly variable in sequence, even between two close relatives. Promiscuous DNA and horizontally transferred sequence constitute part of the intergenic regions, but the origin of the majority of these regions is unknown. On the other hand, duplication and extensive rearrangement of preexisting sequences may be one of the explanations for the occurrence of unknown sequences. Functional aspects of the mitochondrial genome, such as RNA editing and expression of unique open reading frames (ORFs), can be changed under certain nuclear genotypes.     

人与小鼠核糖体蛋白基因内含子中的转录调控位点分析

前期对酵母和果蝇核糖体蛋白(Ribosomal protein, RP)基因内含子序列中的寡核苷酸分析表明, 内含子中含有潜在的转录因子结合位点。为进一步发掘核糖体蛋白基因内含子参与转录调控的证据, 文章首先基于频率分析方法抽提出人和小鼠核糖体蛋白基因第一内含子中高频(Over-represented)出现的寡核苷酸片段 (亦称模体, Motif), 这些寡核苷酸中超过85%与已知的转录因子结合位点吻合, 是潜在的转录调控元件。对抽提出的寡核苷酸进行碱基组成分析, 发现95%以上的寡核苷酸富含碱基C和G, 而较少富含A和T。从寡核苷酸在内含子中的分布情况看, 它们相对靠近第一内含子的5′端, 即距离基因转录起始位点和上游区域较近。推测这些特征可能与基因转录调控有关。     

Chloroplast genome variation in upland and lowland switchgrass

Switchgrass (Panicum virgatum L.) exists at multiple ploidies and two phenotypically distinct ecotypes. To facilitate interploidal comparisons and to understand the extent of sequence variation within existing breeding pools, two complete switchgrass chloroplast genomes were sequenced from individuals representative of the upland and lowland ecotypes. The results demonstrated a very high degree of conservation in gene content and order with other sequenced plastid genomes. The lowland ecotype reference sequence (Kanlow Lin1) was 139,677 base pairs while the upland sequence (Summer Lin2) was 139,619 base pairs. Alignments between the lowland reference sequence and short-read sequence data from existing sequence datasets identified as either upland or lowland confirmed known polymorphisms and indicated the presence of other differences. Insertions and deletions principally occurred near stretches of homopolymer simple sequence repeats in intergenic regions while most Single Nucleotide Polymorphisms (SNPs) occurred in intergenic regions and introns within the single copy portions of the genome. The polymorphism rate between upland and lowland switchgrass ecotypes was found to be similar to rates reported between chloroplast genomes of indica and japonica subspecies of rice which were believed to have diverged 0.2-0.4 million years ago.