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.     

Pattern of selective constraint in C. elegans and C. briggsae genomes.

Similarity between related genomes may carry information on selective constraint in each of them. We analysed patterns of similarity between several homologous regions of Caenorhabditis elegans and C. briggsae genomes. All homologous exons are quite similar. Alignments of introns and of intergenic sequences contain long gaps, segments where similarity is low and close to that between random sequences aligned using the same parameters, and segments of high similarity. Conservative estimates of the fractions of selectively constrained nucleotides are 72%, 17% and 18% for exons, introns and intergenic sequences, respectively. This implies that the total number of constrained nucleotides within non-coding sequences is comparable to that within coding sequences, so that at least one-third of nucleotides in C. elegans and C. briggsae genomes are under strong stabilizing selection.     

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.     

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

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

Haplotype block partition with limited resources and applications to human chromosome 21 haplotype data

Recent studies have shown that the human genome has a haplotype block structure such that it can be decomposed into large blocks with high linkage disequilibrium (LD) and relatively limited haplotype diversity, separated by short regions of low LD. One of the practical implications of this observation is that only a small fraction of all the single-nucleotide polymorphisms (SNPs) (referred as "tag SNPs") can be chosen for mapping genes responsible for human complex diseases, which can significantly reduce genotyping effort, without much loss of power. Algorithms have been developed to partition haplotypes into blocks with the minimum number of tag SNPs for an entire chromosome. In practice, investigators may have limited resources, and only a certain number of SNPs can be genotyped. In the present article, we first formulate this problem as finding a block partition with a fixed number of tag SNPs that can cover the maximal percentage of the whole genome, and we then develop two dynamic programming algorithms to solve this problem. The algorithms are sufficiently flexible to permit knowledge of functional polymorphisms to be considered. We apply the algorithms to a data set of SNPs on human chromosome 21, combining the information of coding and noncoding regions. We study the density of SNPs in intergenic regions, introns, and exons, and we find that the SNP density in intergenic regions is similar to that in introns and is higher than that in exons, results that are consistent with previous studies. We also calculate the distribution of block break points in intergenic regions, genes, exons, and coding regions and do not find any significant differences.     

A statistical analysis of nucleotide sequences of introns and exons in human genes

DNA sequences of 56 human genes for which information on both exons and introns was available were examined. The variance in G+C content among genes is estimated and shown to be substantial. There is a high correlation in G+C content between exons and introns within the same gene. The dinucleotide frequencies of introns are similar to those of intergenic spacer regions and are in reasonable agreement with predictions from substitution rates estimated from pseudogenes, except that the observed deficiency of TA doublets is not predicted. Duplicated bases also show a frequency greater than the expectation under independence. There is marked variability among genes in the frequency of the doublet CG relative to its expectation under independence. This variation is evolutionarily conserved and is correlated with the G+C content. Pseudogenes behave as if they are in a low -G+C, CG-deficient part of the genome, although the genes from which they arose are variable in these respects.      

Population,evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.     

The Euglena gracilis chloroplast ribulose-1,5-bisphosphate carboxylase gene. I. Complete DNA sequence and analysis of the nine intervening sequences

The nucleotide sequence of 6225 base pairs (bp) of Euglena gracilis chloroplast DNA including the complete DNA sequence of the chloroplast-encoded ribulose-1,5-bisphosphate carboxylase large subunit gene along with the flanking DNA sequences is presented. The gene is greater than 5.5 kilobase pairs in length and is organized as 10 exons coding for 475 amino acids, separated by 9 introns. The exons range in size from 45 to 438 bp, while the introns range in size from 382 to 568 bp. The introns have highly conserved boundary sequences with the consensus, 5'-N GTGTGGATTT...(intron)...TTAATTTTAT N-3'. The introns are 82-85 mol% AT, with a pronounced T greater than A greater than G greater than C base bias in the RNA-like strand. They do not appear to encode any polypeptides. In addition, the introns have a conserved sequence 30-50 bp from their 3'-ends with the consensus, 5'-TACAGTTTGAAAATGA-3'. The 5'-TACA sequence bears some homology to the 5'-end of the TACTAACA sequence found in a similar location in yeast nuclear mRNA introns. The conserved sequences of the Euglena rbcL introns may be indicative of a splicing mechanism similar to that of eucaryotic nuclear mRNA introns and group II mitochondrial introns.     

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

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

Genomic MRI - a Public Resource for Studying Sequence Patterns within Genomic DNA

Non-coding genomic regions in complex eukaryotes, including intergenic areas, introns, and untranslated segments of exons, are profoundly non-random in their nucleotide composition and consist of a complex mosaic of sequence patterns. These patterns include so-called Mid-Range Inhomogeneity (MRI) regions -- sequences 30-10000 nucleotides in length that are enriched by a particular base or combination of bases (e.g. (G+T)-rich, purine-rich, etc.). MRI regions are associated with unusual (non-B-form) DNA structures that are often involved in regulation of gene expression, recombination, and other genetic processes (Fedorova & Fedorov 2010). The existence of a strong fixation bias within MRI regions against mutations that tend to reduce their sequence inhomogeneity additionally supports the functionality and importance of these genomic sequences (Prakash et al. 2009).Here we demonstrate a freely available Internet resource -- the Genomic MRI program package -- designed for computational analysis of genomic sequences in order to find and characterize various MRI patterns within them (Bechtel et al. 2008). This package also allows generation of randomized sequences with various properties and level of correspondence to the natural input DNA sequences. The main goal of this resource is to facilitate examination of vast regions of non-coding DNA that are still scarcely investigated and await thorough exploration and recognition.     

The evolution of eukaryotic ribosomal DNA

Mutations occur randomly throughout the ribosomal DNA (rDNA) sequence. Molecular drive (unequal crossing-over, gene conversion, and transposition) spreads these variations through the multiple copies of rDNA. Forces of selection act upon the variants to favor and fix them or disfavor and eliminate them. Selection has not permitted changes in regions within rRNA vital for its function; these sequences are evolutionarily conserved between diverse species. Possible functions for some of these conserved sequences are discussed. The secondary structure of rRNA is also highly conserved during evolution. However, eukaryotic rRNA is larger than prokaryotic rRNA due to blocks of "expansion segments". Arguments are put forward that expansion segments might not play any functional role. Other examples are reviewed of rDNA sequence insertion or deletion, including introns and the internal transcribed spacer 2.     

Transmission of the yeast mitochondrial genome to progeny: The impact of intergenic sequences

Summary In a previous publication it was shown that the output of yeast mitochondrial loci lacking nearby intergenic sequences (encompassing ori/rep elements) was reduced in crosses to strains with wild-type mtDNAs. In the present work, mitochondrial genomes carrying the intergenic deletions were marked at unlinked, loci by introducing specific antibiotic resistance mutations against erythromycin, oligomycin and paromomycin. These marked genomes were used to follow the output of unlinked regions of the genome from crosses between the intergenic deletion mutants and wild-type strains. Transmission of genetically unlinked markers in coding regions was substantially reduced when an intergenic deletion was present on the same genome. In general the transmission of the antibiotic markers was the same as or slightly higher than the corresponding intergenic marker. These results indicate that the presence of an intergenic deletion in the regions studied impairs the transmission to progeny of a mitochondrial genome as a whole. More specifically, the results suggest that ori/rep sequences, present in the regions that have been deleted, confer a competitive advantage over genomes lacking a full complement of such sequences. These results support the hypothesis that intergenic sequences, and specifically ori/rep elements, have a biological role in the mitochondrial genome. However, because of the exclusive presence of ori/rep sequences in the genus Saccharomyces, it may be that these sequences evolved in (or invaded) the mitochondrial genome relatively late in the evolution of the yeasts. Therefore, in a more general sense, variations in the amount and structure of intergenic sequences in various yeasts may reflect processes that have been of selective advantage in the metabolism of individual mitochondrial DNA in a particular environment and that have not drastically interrupted the respiratory phenotype.     

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 compositional properties of human genes

Summary The present work represents the first attempt to study in greater detail previously proposed compositional correlations in genomes, based on a body of additional data relating to gene localizations as well as to extended flanking sequences extracted from gene banks. We have investigated the correlations that exist between (1) the GC levels of exons of human genes, and (2) the GC levels of either intergenic sequences or introns associated with the genes under consideration. In both cases, linear relationships with slopes close to unity were found. The similarity of the linear relationships indicates similar GC levels in intergenic sequences and introns located in the same isochores. Moreover, both intergenic sequences and introns showed GC levels 5–10% lower than the corresponding exons. The above findings considerably strengthen the previously drawn conclusion that coding and noncoding sequences (both inter- and intragenic) from the same isochores of the human genome are compositionally correlated. In addition, we find linear correlations between the GC levels of codon positions and of the intergenic sequences or introns associated with the corresponding genes, as well as among the GC levels of codon positions of genes.     

Structure and evolution of the Xenopus laevis albumin genes

The 68K and 74K albumin genes of Xenopus laevis arose by duplication approximately 30 million years ago. Electron microscopic analysis showed that both genes contain 15 coding sequences. The lengths of corresponding coding sequences are almost identical and are extremely similar to those of mammalian albumin genes. A block of four coding sequences, which in mammals codes for one protein domain, is repeated three times. The corresponding introns are usually different in length and have therefore diverged as a result of insertion/deletion events. The extensive homology between these gene sequences is neither confined to nor most extensive in the coding sequences and similar amounts of homologous sequences are found in the flanking DNAs as in the gene regions. Various structures were formed in the 5'-flanking DNA by mutually exclusive pairing of different homology regions. Analysis of the two 74K albumin gene sequences isolated suggests that the X. laevis genome may contain one 68K albumin gene and two very closely related 74K albumin genes.     

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.