Nucleotide divergence of the rp49 gene region between Drosophila melanogaster and two species of the Obscura group of Drosophila

Summary A 2.1-kb SStI fragment including the rp49 gene and the 3 end of the -serendipity gene has been cloned and sequenced in Drosophila pseudoobscura. rp49 maps at region 62 on the tip of chromosome II of this species. Both the coding and flanking regions have been aligned and compared with those of D. subobscura. There is no evidence for heterogeneity in the rate of silent substitution between the rp49 coding region and the rate of substitutions in flanking regions, the overall silent divergence per site being 0.19. Noncoding regions also differ between both species by different insertions/deletions, some of which are related to repeated sequences. The rp49 region of D. pseudoobscura shows a strong codon bias similar to those of D. subobscura and D. melanogaster. Comparison of the rates of silent (K _S ) and nonsilent (K _a ) substitutions of the rp49 gene and other genes completely sequenced in D. pseudoobscura and D. melanogaster confirms previous results indicating that rp49 is evolving slowly both at silent and nonsilent sites. According to the data for the rp49 region, D. pseudoobscura and D. subobscura lineages would have diverged some 9 Myr ago, if one assumes a divergence time of 30 Myr for the melanogaster and obscura groups.Offprint requests to: C. Segarra     

The factors dictating the codon usage variation among the genes in the genome of Burkholderia pseudomallei

Burkholderia pseudomallei is a recognized biothreat agent and the causative agent of melioidosis. Codon usage biases of all protein-coding genes (length greater than or equal to 300 bp) from the complete genome of B. pseudomallei K96243 have been analyzed. As B. pseudomallei is a GC-rich organism (68.5%), overall codon usage data analysis indicates that indeed codons ending in G and/or C are predominant in this organism. But multivariate statistical analysis indicates that there is a single major trend in the codon usage variation among the genes in this organism, which has a strong positively correlation with the expressivities of the genes. The majority of the lowly expressed genes are scattered towards the negative end of the major axis whereas the highly expressed genes are clustered towards the positive end. At the same time, from the results that there were two significant correlations between axis 1 coordinates and the GC, GC3s content at silent sites of each sequence, and clearly significant negatively correlations between the ‘Effective Number of Codons’ values and GC, GC3s content, we inferred that codon usage bias was affected by gene nucleotide composition also. In addition, some other factors such as the lengths of the genes as well as the hydrophobicity of genes also influence the codon usage variation among the genes in this organism in a minor way. At the same time, notably, 21 codons have been defined as ‘optimal codons’ of the B. pseudomallei. In summary, our work have provided a basic understanding of the mechanisms for codon usage bias and some more useful information for improving the expression of target genes in vivo and in vitro. Sheng Zhao and Qin Zhang contributed equally to this work.     

Nonneutral GC3 and Retroelement Codon Mimicry in Phytophthora

Phytophthora is a genus entirely comprised of destructive plant pathogens. It belongs to the Stramenopila, a unique branch of eukaryotes, phylogenetically distinct from plants, animals, or fungi. Phytophthora genes show a strong preference for usage of codons ending with G or C (high GC3). The presence of high GC3 in genes can be utilized to differentiate coding regions from noncoding regions in the genome. We found that both selective pressure and mutation bias drive codon bias in Phytophthora. Indicative for selection pressure is the higher GC3 value of highly expressed genes in different Phytophthora species. Lineage specific GC increase of noncoding regions is reminiscent of whole-genome mutation bias, whereas the elevated Phytophthora GC3 is primarily a result of translation efficiency-driven selection. Heterogeneous retrotransposons exist in Phytophthora genomes and many of them vary in their GC content. Interestingly, the most widespread groups of retroelements in Phytophthora show high GC3 and a codon bias that is similar to host genes. Apparently, selection pressure has been exerted on the retroelement’s codon usage, and such mimicry of host codon bias might be beneficial for the propagation of retrotransposons. Reviewing Editor: Dr. Yves van de Peer     

栽培大豆和野生大豆线粒体基因组密码子使用偏性的比较分析

为分析栽培大豆和野生大豆线粒体基因组的密码子使用特征差异,该文以其线粒体基因组编码序列为研究对象,比较其密码子偏性形成的影响因素和演化过程。结果表明:(1)栽培大豆和野生大豆线粒体基因组编码区的GC含量分别为44.56%和44.58%,说明栽培大豆和野生大豆线粒体编码基因均富含A/T碱基。(2)栽培大豆和野生大豆线粒体基因组密码子第1位、第2位GC含量平均值与第3位GC含量的相关性均呈极显著水平,说明突变在其密码子偏性形成中的作用不可忽略; PR2-plot分析显示,在同义密码子第3位碱基的使用频率上,嘌呤低于嘧啶; Nc-plot分析中Nc比值位于-0.1～0.2区间的基因数占总基因数的95%以上;突变和选择等多重因素共同作用影响了大豆线粒体基因组编码序列密码子使用偏性的形成。(3)有20、21个密码子分别被确定为栽培大豆和野生大豆线粒体基因组编码序列的最优密码子,其中除丝氨酸TCC密码子外均以A或T结尾。综上结果认为,栽培大豆线粒体密码子偏性的形成受选择的影响要高于野生大豆,这可能是栽培大豆由野生大豆经长期人工栽培驯化的结果。     

紫花苜蓿叶绿体基因组密码子偏好性分析

为分析紫花苜蓿叶绿体基因组密码子偏好性的使用模式,该文以紫花苜蓿叶绿体基因组中筛选到的49条蛋白质编码序列为研究对象,利用CodonW、CUSP、CHIPS、SPSS等软件对其密码子的使用模式和偏好性进行研究。结果表明:(1)紫花苜蓿叶绿体基因的第3位密码子的平均GC含量为26.44%,有效密码子数(ENC)在40.6～51.41之间,多数密码子的偏好性较弱。(2)相对同义密码子使用度(RSCU)分析发现,RSCU>1 的密码子数目有30个,以A、U结尾的有29个,说明了紫花苜蓿叶绿体基因组A或U出现的频率较高。(3)中性分析发现,GC3与 GC12的相关性不显著,表明密码子偏性主要受自然选择的影响; ENC-plot 分析发现一部分基因落在曲线的下方及周围,表明突变也影响了部分密码子偏性的形成。此外,有17个密码子被鉴定为紫花苜蓿叶绿体基因组的最优密码子。紫花苜蓿叶绿体基因组的密码子偏好性可能受自然选择和突变的共同作用。该研究将为紫花苜蓿叶绿体基因工程的开展和目标性状的遗传改良奠定基础。     

On the rate of DNA sequence evolution inDrosophila

Summary Analysis of the rate of nucleotide substitution at silent sites inDrosophila genes reveals three main points. First, the silent rate varies (by a factor of two) among nuclear genes; it is inversely related to the degree of codon usage bias, and so selection among synonymous codons appears to constrain the rate of silent substitution in some genes. Second, mitochondrial genes may have evolved only as fast as nuclear genes with weak codon usage bias (and two times faster than nuclear genes with high codon usage bias); this is quite different from the situation in mammals where mitochondrial genes evolve approximately 5–10 times faster than nuclear genes. Third, the absolute rate of substitution at silent sites in nuclear genes inDrosophila is about three times hihger than the average silent rate in mammals.     

沙枣叶绿体全基因组序列及其使用密码子偏性分析

该研究以2株野生沙枣(Elaeagnus angustifolia Linn.)嫩枝经温室水培后的嫩叶为材料,采用CTAB法分别提取总DNA,并利用第二代测序技术进行总DNA从头测序,组装后得到2株沙枣叶绿体基因组全序列,并详细分析了其蛋白质编码基因密码子使用的偏好性及其原因,为沙枣叶绿体基因工程和分子系统进化等研究奠定基础。结果显示:(1)组装得到沙枣叶绿体基因组序列全长150 546 bp,由长度为81 113 bp的长单拷贝(LSC)区域和25 494 bp的短单拷贝(SSC)区域,以及1对分隔开它们的长18 445 bp的反向重复序列(IRS)组成;注释共得到132个基因,包括86个蛋白编码基因、38个tRNA基因和8个rRNA基因。(2)沙枣叶绿体基因组蛋白编码基因密码子的第三位碱基GC含量(GC_3)为28.47%,明显低于整个叶绿体基因组GC含量(37%),也低于第一位(GC_1)和第二位(GC_2)碱基的GC含量,说明密码子对AT碱基结尾有偏好性;其中, UCU、CCU、UGU、GCU、CUU、GAU、UCA和UAA为最优密码子。(3)同义密码子相对使用频率(RSCU)分析发现,影响密码子使用模式的因素并不单一,密码子的偏好性受到突变、选择及其他因素的共同影响,并且自然选择表达引起的序列差异比突变对密码子偏好性的影响要显著;中性绘图分析、有效密码子数(ENC-plot)分析和奇偶偏好性(PR2-plot)分析表明,沙枣叶绿体基因组使用密码子的偏性受选择的影响更大。(4)通过最大似然法、最大简约法和贝叶斯方法对胡颓子科6个物种和1个枣的叶绿体基因序列构建系统发育树,与它们使用密码子偏性聚类的结果一致,表明叶绿体基因组使用密码子偏性与物种的亲缘关系相关。     

The Effect of Local Nucleotides on Synonymous Codon Usage in the Honeybee (<Emphasis Type="Italic">Apis mellifera</Emphasis> L.) Genome

Using all currently predicted coding regions in the honeybee genome, a novel form of synonymous codon bias is presented that affects the usage of particular codons dependent on the surrounding nucleotides in the coding region. Nucleotides at the third codon site are correlated, dependent on their weak (adenine [A] or thyamine [T]) versus strong (guanine [G] or cytosine [C]) status, to nucleotides on the first codon site which are dependent on their purine (A/G) versus pyrimidine (C/T) status. In particular, for adjacent third and first site nucleotides, weak–pyrimidine and strong–purine nucleotide combinations occur much more frequently than the underabundant weak–purine and strong–pyrimidine nucleotide combinations. Since a similar effect is also found in the noncoding regions, but is present for all adjacent nucleotides, this coding effect is most likely due to a genome-wide context-dependent mutation error correcting mechanism in combination with selective constraints on adjacent first and second nucleotide pairs within codons. The position-dependent relationship of synonymous codon usage is evidence for a novel form of codon position bias which utilizes the redundancy in the genetic code to minimize the effect of nucleotide mutations within coding regions. [Reviewing Editor: Dr. Brian Morton]     

A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences—exons in vertebrates and genes in prokaryotes—are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content. Correspondence to: J. L. Oliver     

Intralineage variation in the pattern ofrbcL nucleotide substitution

Variation in chloroplastrbcL sequences was studied in representative species of four different lineages: the tribeRubieae (Rubiaceae), and the generaDrosera (Droseraceae),Nothofagus (Nothofagaceae) andIlex (Aquifoliaceae). Each lineage has its particular non-overlapping set ofrbcL polymorphic sites, indicating that common unconstrainedrbcL sites are not shared. Large differences in the rate and pattern of nucleotide substitution are observed among the four lineages. The genusIlex has the lowest rate of substitution, the lowest transition/transversion ratio, the lowest synonymous/replacement ratio and the lowest number of substitutions at the third codon position. An apparent relationship of these measures to the age of the lineages is observed. The A + T content and codon use among the four lineages are very similar and, apparently, cannot account for the observed differences in patterns of nucleotide substitution. However, the A + T content of the two bases immediately flanking the polymorphic sites is higher inIlex than in the other lineages. This could be correlated with the transversion/transition bias observed inIlex. The particularly low synonymous/replacement ratio found inIlex could also be explained by the small population sizes of species in this genus.     

杜仲基因密码子使用模式分析

为了解杜仲基因密码子使用模式,该文以杜仲基因组密码子为研究对象,运用CodonW软件对杜仲的320个蛋白编码基因进行同义密码子相对使用频率(RSCU)分析、ENC-GC3s关联分析编码基因的密码子ENC值、PR2-plot偏倚分析编码基因的密码子碱基使用频率,并运用CUSP软件与Codon Usage Database...     

The Role of Context-Dependent Mutations in Generating Compositional and Codon Usage Bias in Grass Chloroplast DNA

Abstract The influence of local base composition on mutations in chloroplast DNA (cpDNA) is studied in detail and the resulting, empirically derived, mutation dynamics are used to analyze both base composition and codon usage bias. A 4 × 4 substitution matrix is generated for each of the 16 possible flanking base combinations (contexts) using 17,253 noncoding sites, 1309 of which are variable, from an alignment of three complete grass chloroplast genome sequences. It is shown that substitution bias at these sites is correlated with flanking base composition and that the A+T content of these flanking sites as well as the number of flanking pyrimidines on the same strand appears to have general influences on substitution properties. The context-dependent equilibrium base frequencies predicted from these matrices are then applied to two analyses. The first examines whether or not context dependency of mutations is sufficient to generate average compositional differences between noncoding cpDNA and silent sites of coding sequences. It is found that these two classes of sites exist, on average, in very different contexts and that the observed mutation dynamics are expected to generate significant differences in overall composition bias that are similar to the differences observed in cpDNA. Context dependency, however, cannot account for all of the observed differences: although silent sites in coding regions appear to be at the equilibrium predicted, noncoding cpDNA has a significantly lower A+T content than expected from its own substitution dynamics, possibly due to the influence of indels. The second study examines the codon usage of low-expression chloroplast genes. When context is accounted for, codon usage is very similar to what is predicted by the substitution dynamics of noncoding cpDNA. However, certain codon groups show significant deviation when followed by a purine in a manner suggesting some form of weak selection other than translation efficiency. Overall, the findings indicate that a full understanding of mutational dynamics is critical to understanding the role selection plays in generating composition bias and sequence structure.     

Recessive allelic variations of three microsatellite sites within the<Emphasis Type="Italic">O2</Emphasis> gene in maize

Recessive allelic variations were investigated at 3 microsatellite (SSR) sites within theO2 gene by using 14 inbredo2 lines and a wild-type line in maize. Among the 15 lines, allelic variations were observed at umc1066, phi057, and phi112 sites. Two alleles were found at the umc1066 site—a recessive allele with 2 perfect GCCAGA repeats and a dominant allele with 3 perfect repeats. Three alleles were found at the phi057 site—2 recessive alleles with 3 and 5 perfect GCC repeats, respectively, and another with 4 perfect repeats consistent with a dominant allele. At least 4 alleles exist at the phi112 site—among which 1 recessive allele has a 1-bp deletion, another has a 15-bp deletion, and other has no PCR products compared to the dominant allele; all the alleles have unchanged AG repeats. The phi057 site in exon 6 was identified to be a hypervariable region in the coding sequence of the02 gene, in addition to the 2 hypervariable regions in exon 1 previously reported. The primary mechanisms underlying the variations in repeat numbers and regions flanking the SSR within theO2 gene appear to be unequal crossing over and replication slippage. Furthermore, base substitution of SSR motif can create heteroalleles and modify the repeat number of SSR. The lysine content of kernel in theO2 ando2 lines correlates to a considerable extent with nucleotide variations at the umc1066, phi057, and phi112 sites. Our study suggests that it is best to use the 3 markers together in molecular marker-assisted selection for high-lysine maize materials.     

Comparison of dinucleotide frequency and codon usage inToxoplasma andPlasmodium: Evolutionary implications

Summary The weight-averaged observed/expected dinucleotide frequencies for the sum total of the coding regions of fiveToxoplasma genes were compared with the same parameters previously determined for the coding regions of 21Plasmodium genes. In addition, codon usage in the fiveToxoplasma genes was compared with that in the 21Plasmodium genes, and the percent distribution of amino acids in theToxoplasma protein pool and thePlasmodium protein pool were compared with that in a general protein pool of 314 proteins. The results are consistent with the hypothesis that, contrary to currently held opinion, the generaToxoplasma andPlasmodium are not especially closely related.     

The mitochondrial genome of <Emphasis Type="Italic">Protohermes concolorus</Emphasis> Yang et Yang 1988 (Insecta: Megaloptera: Corydalidae)

The first complete mitochondrial genome of dobsonfly Protohermes concolorus Yang et Yang, 1988 (Megaloptera: Corydalidae) was sequenced in this study. The genome was a circular molecule of 15,851 bp containing the typical 37 genes that arranged in the same order as that of the putative ancestor of hexapods. Sequences overlaps were observed between several neighbor genes, which made the genome relatively compact. The tRNA-Ser (GCT) could not be folded into typical secondary structure because its DHU arm was replaced with a simple loop. Six of the 13 protein genes were terminated with a single T adjacent to a downstream tRNA gene in the same strand. The variation of GC content caused the different nucleotide substitution patterns of the protein genes. The genome was AT-biased with a total A + T content of 75.83% which was also demonstrated by the codon usage. The control region was the most AT-rich region with a sub-region of even higher A + T content. Protein genes of two strands presented opposite CG-skew trends which was also reflected by the codon usage. For most of the amino acids, the protein coding sequences did not prefer to use the cognate codons of corresponding tRNAs and the codon usage of the protein genes was not random. The variation of nucleotide substitution patterns of protein genes was significantly correlated with the GC content. The phylogenetic analyses based on all the 13 protein genes showed that Megaloptera was the sister group of other holometabolous insects except Coleoptera.     

Codon usage in Aspergillus nidulans.

Summary Synonymous codon usage in genes from the ascomycete (filamentous) fungus Aspergillus nidulans has been investigated. A total of 45 gene sequences has been analysed. Multivariate statistical analysis has been used to identify a single major trend among genes. At one end of this trend are lowly expressed genes, whereas at the other extreme lie genes known or expected to be highly expressed. The major trend is from nearly random codon usage (in the lowly expressed genes) to codon usage that is highly biased towards a set of 19–20 optimal codons. The G+C content of the A. nidulans genome is close to 50%, indicating little overall mutational bias, and so the codon usage of lowly expressed genes is as expected in the absence of selection pressure at silent sites. Most of the optimal codons are C- or G-ending, making highly expressed genes more G+C-rich at silent sites.     

Codon usage and base composition inRickettsia prowazekii

Codon usage and base composition in sequences from the A + T-rich genome ofRickettsia prowazekii, a member of the alpha Proteobacteria, have been investigated. Synonymous codon usage patterns are roughly similar among genes, even though the data set includes genes expected to be expressed at very different levels, indicating that translational selection has been ineffective in this species. However, multivariate statistical analysis differentiates genes according to their G + C contents at the first two codon positions. To study this variation, we have compared the amino acid composition patterns of 21R. prowazekii proteins with that of a homologous set of proteins fromEscherichia coli. The analysis shows that individual genes have been affected by biased mutation rates to very different extents: genes encoding proteins highly conserved among other species being the least affected. Overall, protein coding and intergenic spacer regions have G + C content values of 32.5% and 21.4%, respectively. Extrapolation from these values suggests thatR. prowazekii has around 800 genes and that 60–70% of the genome may be coding. Correspondence to: S.G.E. Andersson     

Compositional Correlations in the Chicken Genome

This paper analyses the compositional correlations that hold in the chicken genome. Significant linear correlations were found among the regions studied—coding sequences (and their first, second, and third codon positions), flanking regions (5′ and 3′), and introns—as is the case in the human genome. We found that these compositional correlations are not limited to global GC levels but even extend to individual bases. Furthermore, an analysis of 1037 coding sequences has confirmed a correlation among GC₃, GC₂, and GC₁. The implications of these results are discussed. Received: 9 December 1998 / Accepted: 18 April 1999     

The comparative amino acid sequences, substrate specificities and gene or cDNA nucleotide sequences of some prokaryote and eukaryote amidinotransferases: implications for evolution

The amino acid sequences of the amidinotransferases and the nucleotide sequences of their genes or cDNA from four Streptomyces species (seven genes) and from the kidneys of rat, pig, human and human pancreas were compared. The overall amino acid and nucleotide sequences of the prokaryotes and eukaryotes were very similar and further, three regions were identified that were highly identical. Evidence is presented that there is virtually zero chance that the overall and high identity regions of the amino acid sequence similarities and the overall nucleotide sequence similarities between Streptomyces and mammals represent random match. Both rat and lamprey amidinotransferases were able to use inosamine phosphate, the amidine group acceptor of Streptomyces. We have concluded that the structure and function of the amidinotransferases and their genes has been highly conserved through evolution from prokaryotes to eukaryotes. The evolution has occurred with: (1) a high degree of retention of nucleotide and amino acid sequences; (2) a high degree of retention of the primitive Streptomyces guanine+cytosine (G+C) third codon position composition in certain high identity regions of the eukaryote cDNA; (3) a decrease in the specificities for the amidine group acceptors; and (4) most of the mutations silent in the regions suggested to code for active sites in the enzymes.     

香樟转录组基因密码子偏好性分析

为了解香樟基因密码子偏好性,该文以NCBI网站中香樟转录组数据为材料,利用生物信息学手段评价转录组数据质量,选取高质量数据的转录组,去除低质量序列,组装转录组,预测基因结构,再利用自编perl脚本提取以AUG开头的基因序列37 Mb序列34 931个基因,进一步利用CodonW分析基因密码子偏好性。结果表明:GC含量的变化范围为0.273～0.742,均值为0.452; ENC的范围为26.29～61.00,均值为52.76; CAI的范围为0.064～0.401,均值为0.199; RSCU值大于1的密码子数目为27个,其中以U或A结尾的有22个; 中性分析表明,小部分基因在对角线上,大多数基因偏离对角线; ENC-plot分析表明小部分基因在标准曲线上,大多数基因偏离标准曲线。上述研究结果表明,香樟基因的密码子偏好性比较弱,密码子常以A/U结尾; 突变和选择两者都在密码子偏好中起作用,而选择作用更大; 最终确定了GUU、CAG、GAA、UCU、GCU、GGU为最优密码子,通过对目标基因密码子的校正,提高表达效率,从而为利用基因工程技术改良香樟重要性状奠定了基础。