Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.      

Sense codons are found in specific contexts

The sequence environment of codons in structural genes has been investigated statistically, using computer methods. A set of Escherichia coli genes with abundant products was compared with a set having low gene product levels, in order to detect potential differences associated with expression. The results show striking non-randomness in the nucleotides occurring near codons. These effects are, unexpectedly, very much larger and more homogeneous among the genes with rare products. The intensity of effects in weakly expressed genes suggests that such non-random sequence environments decrease expression. In the weakly expressed set of genes, the 5' neighbor of a codon, and all positions of the 3' neighbor codon are biased. In the highly expressed genes, the first nucleotide of the next codon is a uniquely affected site. The distribution of non-randomness in weakly expressed genes suggests that sequence bias is primarily due to a constraint acting directly on the secondary or tertiary structure of the codon/anticodon. In highly expressed genes, the observed bias suggests an interaction between the codon/anticodon and a site outside the codon/anticodon. Much of the tendency to non-random near-neighbor sequences in weakly expressed genes can be ascribed to a correlation between nearby nucleotides and the wobble nucleotide of the codon, despite the fact that selection of such correlations will alter the amino acid sequence. The favored pattern, in genes expressed at low level, is R YYR or Y RRY. R indicates purine, Y indicates pyrimidine; the space is the boundary between codons. It seems likely that this preference for nearby sequences is the physical basis of the genetic context effect. Under this assumption such sequence biases will affect expression. On this basis, we predict new sites for contextual mutations which decrease expression, and suggest strategy for the design of messages having optimal translational activity.     

Isolation and characterization of the genomic region from Drosophila kuntzei containing the Adh and Adhr genes

The nucleotide sequences of the Adh and Adhr genes of Drosophila kuntzei were derived from combined overlapping sequences of clones isolated from a genomic library and from cloned PCR and inverse-PCR fragments. Only a proximal promoter was detected upstream of the Adh gene, indicating that D. kuntzei Adh is regulated by a one-promoter system. Further upstream of the Adh structural gene, an adult enhancer region (AAE) was found that contains most of the regulatory sequences described for AAEs of other Drosophila species. Analysis of the ADH protein showed an amino acid change from valine to threonine in the active site at position 189 which is also found in D. funebris but is otherwise unique among DROSOPHILA: This difference alone may be responsible for the very low ADH activity found in this species and may cause a difference in substrate usage pattern. Codon bias in Adh and Adhr was comparable and found to be very low compared with other species. Phylogenetic analysis showed that D. kuntzei is closest related to D. funebris and D. immigrans. The time of divergence between D. kuntzei and D. funebris was estimated to be 14.2-20.2 Myr and that between D. kuntzei-D. funebris and D. immigrans to be 30.8-44.0 Myr. An analysis of the genetic variation in the Adh gene and upstream sequences of four European strains showed that this gene was highly variable. Overall nucleotide diversity (pi) was 0.0139, which is two times higher than that in D. melanogaster.     

Influence of intercodon and base frequencies on codon usage in filarial parasites

Base frequency, codon usage, and intercodon identity were analyzed in five filarial parasite species representing five Onchocercidae genera. Wucheria bancrofti, Brugia malayi, Onchocerca volvulus, Acanthocheilonema viteae, and Dirofilaria immitis gene sequences were downloaded from NCBI, and analysis was performed using locally designed computer programs and other freely available applications. A clear sequence bias was observed among the nematode species examined. At the nucleotide level, AT basepairs were present in gene sequences at higher frequencies than GC. In addition, codons ending in A or T were used proportionately more than those with G or C in the third-codon position. In addition, the amino acids used most often corresponded to codons ending in AT basepairs. Intercodon base proportion was biased in that A was found most often at N4, second only to T in certain specific cases. Since all of these sequence biases were observed in a relatively consistent fashion among all of the organisms studied, we conclude that sequence bias is a genetic characteristic, which is associated with multiple filarial genera.     

Codon contexts in enterobacterial and coliphage genes

This investigation of the codon context of enterobacteria, plasmid, and phage protein genes was based on a search for correlations between the presence of one base type at codon position III and the presence of another base type at some other position in adjacent codons. Enterobacterial genes were compared with eukaryotic sequences for codon context effects. In enterobacterial genes, base usage at codon position III is correlated with the third position of the upstream adjacent codon and with all three positions of the downstream codon. Plasmid genes are free of context biases. Phage genes are heterogeneous: MS2 codons have no biased context, whereas lambda genes partly follow the trends of the host bacterium, and T7 genes have biased codon contexts that differ from those of the host. It has been reported that two successive third-codon positions tend to be occupied by two purines or two pyrimidines in Escherichia coli genes of low expression level. Here, the extent to which highly expressed protein genes can modulate base usage at two successive codon positions III, given the constraints on codon usage and protein sequence that act on them, was quantified. This demonstrates that the above-mentioned favored patterns are not a characteristic of weakly expressed genes but occur in all genes in which codon context can vary appreciably. The correlation between successive third-codon positions is a distinct feature of enterobacteria and of some phages, one that may result from adaptation of gene structure to translational efficiency. Conversely, codon context in yeast and human genes is biased--but for reasons unrelated to translation.      

Patterns of genetic diversification within the Adh gene family in the grasses (Poaceae).

We investigated the evolutionary dynamics of the Adh gene family within the grasses (Poaceae), with the goal of using molecular evolutionary tools to understand the process of gene family diversification. We analyzed 21 Adh sequences representing a broad array of grasses. Phylogenetic analyses suggested that Adh duplicated into Adh1 and Adh2 before the radiation of the grasses roughly 65 MYA. Gene structure, including intron length, has varied little over this period. Conservation of intron length prompted investigation into the dynamics of intron evolution, particularly the ability of intron sequences to form secondary structures. Intron sequences did not have an extremely high or low minimum free energy of folding relative to permuted sequences, suggesting that individual Adh introns do not evolve under secondary structural constraints. For coding sequences, the diversification of Adh1 and Adh2 was marked by a shift in third-position G + C content. This shift may reflect differential selection for codon use. Diversification between Adh1 and Adh2 was also typified by a shift in nonsynonymous nucleotide substitution rates, but there was no evidence that relatively fast nonsynonymous nucleotide substitution rates in the Adh2 clade were a product of diversifying selection. Gene conversion may have played a role in retarding diversification of Adh1 and Adh2 in rice, but there is no evidence of gene conversion between paralogs in other taxa. Although the reasons for retention of two functional Adh genes remain obscure, we propose that a shift in gene expression was important for the retention of the two Adh gene copies within the grasses.     

Selection intensity on preferred codons correlates with overall codon usage bias in Caenorhabditis remanei

Adaptive codon usage provides evidence of natural selection in one of its most subtle forms: a fitness benefit of one synonymous codon relative to another. Codon usage bias is evident in the coding sequences of a broad array of taxa, reflecting selection for translational efficiency and/or accuracy as well as mutational biases. Here, we quantify the magnitude of selection acting on alternative codons in genes of the nematode Caenorhabditis remanei, an outcrossing relative of the model organism C. elegans, by fitting the expected mutation-selection-drift equilibrium frequency distribution of preferred and unpreferred codon variants to the empirical distribution. This method estimates the intensity of selection on synonymous codons in genes with high codon bias as N(e)s = 0.17, a value significantly greater than zero. In addition, we demonstrate for the first time that estimates of ongoing selection on codon usage among genes, inferred from nucleotide polymorphism data, correlate strongly with long-term patterns of codon usage bias, as measured by the frequency of optimal codons in a gene. From the pattern of polymorphisms in introns, we also infer that these findings do not result from the operation of biased gene conversion toward G or C nucleotides. We therefore conclude that coincident patterns of current and ancient selection are responsible for shaping biased codon usage in the C. remanei genome.     

Relationships Among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

It is well known that stop codons play a critical role in the process of protein synthesis. However, little effort has been made to investigate whether stop codon usage exhibits biases, such as widely seen for synonymous codon usage. Here we systematically investigate stop codon usage bias in various eukaryotes as well as its relationships with its context, GC3 content, gene expression level, and secondary structure. The results show that there is a strong bias for stop codon usage in different eukaryotes, i.e., UAA is overrepresented in the lower eukaryotes, UGA is overrepresented in the higher eukaryotes, and UAG is least used in all eukaryotes. Different conserved patterns for each stop codon in different eukaryotic classes are found based on information content and logo analysis. GC3 contents increase with increasing complexity of organisms. Secondary structure prediction revealed that UAA is generally associated with loop structures, whereas UGA is more uniformly present in loop and stem structures, i.e., UGA is less biased toward having a particular structure. The stop codon usage bias, however, shows no significant relationship with GC3 content and gene expression level in individual eukaryotes. The results indicate that genomic complexity and GC3 content might contribute to stop codon usage bias in different eukaryotes. Our results indicate that stop codons, like synonymous codons, exhibit biases in usage. Additional work will be needed to understand the causes of these biases and their relationship to the mechanism of protein termination. [Reviewing Editor: Dr. Manyuan Long]     

Nucleotide distribution and the recognition of coding regions in DNA sequences: an information theory approach

A method for the recognition of coding-regions along DNA sequences is described. The method is based on the observation, made in several cases, that nucleotide distribution at the third position of the codon is more biased (less random) than that in the other two positions. It is suggested that since nucleotide distribution at the third position is only weakly influenced by the amino acid distribution in the coded protein, there must be some constraints at the DNA level which bias the nucleotide distribution at the third position. The distinction between DNA-level constraints and protein-level constraints is discussed in the frame of Information Theory, and the analysis of the Mitochondrial gene coding for subunit-1 of the yeast cytochrome oxidase is presented.     

Analysis of codon usage pattern in the radioresistant bacterium Deinococcus radiodurans

The main factors shaping codon usage bias in the Deinococcus radiodurans genome were reported. Correspondence analysis (COA) was carried out to analyze synonymous codon usage bias. The results showed that the main trend was strongly correlated with gene expression level assessed by the "Codon Adaptation Index" (CAI) values, a result that was confirmed by the distribution of genes along the first axis. The results of correlation analysis, variance analysis and neutrality plot indicated that gene nucleotide composition was clearly contributed to codon bias. CDS length was also key factor in dictating codon usage variation. A general tendency of more biased codon usage of genes with longer CDS length to higher expression level was found. Further, the hydrophobicity of each protein also played a role in shaping codon usage in this organism, which could be confirmed by the significant correlation between the positions of genes placed on the first axis and the hydrophobicity values (r=-0.100, P<0.01). In summary, gene expression level played a crucial role, nucleotide mutational bias, CDS length and the hydrophobicity of each protein just in a minor way in shaping the codon usage pattern of D. radiodurans. Notably, 19 codons firstly defined as "optimal codons" may provide useful clues for molecular genetic engineering and evolutionary studying.     

Evolution of the Adh locus in the Drosophila willistoni group: the loss of an intron, and shift in codon usage

We report here the DNA sequence of the alcohol dehydrogenase gene (Adh) cloned from Drosophila willistoni. The three major findings are as follows: (1) Relative to all other Adh genes known from Drosophila, D. willistoni Adh has the last intron precisely deleted; PCR directly from total genomic DNA indicates that the deletion exists in all members of the willistoni group but not in any other group, including the closely related saltans group. Otherwise the structure and predicted protein are very similar to those of other species. (2) There is a significant shift in codon usage, especially compared with that in D. melanogaster Adh. The most striking shift is from C to U in the wobble position (both third and first position). Unlike the codon-usage-bias pattern typical of highly biased genes in D. melanogaster, including Adh, D. willistoni has nearly 50% G + C in the third position. (3) The phylogenetic information provided by this new sequence is in agreement with almost all other molecular and morphological data, in placing the obscura group closer to the melanogaster group, with the willistoni group farther distant but still clearly within the subgenus Sophophora.      

Molecular population genetics of the alcohol dehydrogenase locus in the Hawaiian drosophilid D. mimica

Sequence variation among 10 alleles of the alcohol dehydrogenase (Adh) gene of the Hawaiian drosophilid D. mimica was analyzed with reference to the evolutionary history of the Hawaiian subgroup as well as to levels and patterns of polymorphism of the Adh gene in continental drosophilid species. The Adh gene of D. mimica is less polymorphic than that of other drosophilid species, and no replacement substitutions were found. Statistical analyses of the Adh alleles suggested the action of balancing selection and revealed significant linkage disequilibrium among three of the variable sites. The effective population size was estimated to be only slightly smaller than that of continental species and, surprisingly, on the same order of magnitude as the actual size.      

Cloning and characterization of Vine-1, a LTR-retrotransposon-like element in Vitis vinifera L., and other Vitis species.

We report the organization of a grapevine chimeric gene Adhr-Vine-1, composed by an Adhr gene, into which a retroelement, Vine-1, was inserted. Sequence analysis revealed that Adhr is a member of the Adh multigene family, but does not correspond to any other grapevine Adh described to date. Vine-1, albeit defective, is the most complete LTR (long terminal repeat)-retrotransposon-like element described in Vitis vinifera L. It is 2392 bp long, with two almost identical LTRs (287 bp) in the same orientation, and flanked by direct repeats of a 5 bp host DNA. This element presents other features, characteristic of retroviruses and retrotransposons including inverted repeats, a primer binding site, and a polypurine tract. It has a single open reading frame (ORF) of 581 amino acids, potentially encoding for a gag protein and parts of the protease and integrase proteins. Vine-1 is most likely related to the copia-like type family, but with no significant similarity to any previously described plant retrotransposon or inserted element, nor to any eukaryotic element described to date. Vine-1 element has been found in Adhr at the same location in different V. vinifera cultivars, but not in some other analyzed Vitis species. These data suggest that Vine-1 insertion in Adhr is specific to V. vinifera, and has occurred after the Adh isogene separation, but prior to cultivar development. Sequences related to Vine-1 were revealed in multiple copies in the V. vinifera genome and, to a lesser extent, in other analyzed Vitis species. The polymorphism observed prompts us to question the role played by transposition in the evolution of the Vitis genus.     

Determining the relationship of gene expression and global mRNA stability in Drosophila melanogaster and Escherichia coli using linear models

There are several sequence-dependent factors regulating gene expression. Some of them have been extensively studied, among the most prominent are GC content and codon usage bias. Other factors hypothesized to have an impact on gene expression are gene length and the thermodynamic stability of mRNA secondary structure. In this work, we analyzed two different microarray datasets of Drosophila melanogaster gene expression and one dataset of Escherichia coli. To investigate the relationship between gene expression, codon usage bias and GC content of first, second and third codon position, gene length and mRNA stability we employed a multiple regression analysis using a comprehensive linear model. It is shown that codon usage bias and GC content of the first, second and third codon position show a significant influence on gene expression, whereas no significant effect of mRNA secondary structure stability is observed.     

Translational effects of differential codon usage among intragenic domains of new genes in Drosophila

Evolved codon usages often pose a technical challenge over the expressing of eukaryotic genes in microbial systems because of changed translational machinery. In the present study, we investigated the translational effects of intragenic differential codon usage on the expression of the new Drosophila gene, jingwei (jgw), a chimera derived from two unrelated parental genes: Ymp and Adh. We found that jgw possesses a strong intragenic differential usage of synonymous codons, i.e. the Adh-derived C-domain has a significantly higher codon bias than that of the Ymp-derived N-domain (P=0.0023 by t-test). Additional evolutionary analysis revealed the heterogeneous distribution of rare codons, implicating its role in gene regulation and protein translation. The in vitro expression of jgw further demonstrated that the heterogeneous distribution of rare codons has played a role in regulating gene expression, particularly, affecting the quality of protein translation.