Sequence and structural aspects of the functional diversification of plant alcohol dehydrogenases

The glycolytic proteins in plants are coded by small multigene families, which provide an interesting contrast to the high copy number of gene families studied to date. The alcohol dehydrogenase (Adh) genes encode glycolytic enzymes that have been characterized in some plant families. Although the amino acid sequences of zinc-containing long-chain ADHs are highly conserved, the metabolic function of this enzyme is variable. They also have different patterns of expression and are submitted to differences in nonsynonymous substitution rates between gene copies. It is possible that the Adh copies have been retained as a consequence of adaptative amino acid replacements which have conferred subtle changes in function. Phylogenetic analysis indicates that there have been a number of separate duplication events within angiosperms, and that genes labeled Adh1, Adh2 and Adh3 in different groups may not be homologous. Nonsynonymous/synonymous ratios yielded no signs of positive selection. However, the coefficients of functional divergence (theta) estimated between the Adh1 and Adh2 gene groups indicate statistically significant site-specific shift of evolutionary rates between them, as well as between those of different botanical families, suggesting that altered functional constraints may have taken place at some amino acid residues after their diversification. The theoretical three-dimensional structure of the alcohol dehydrogenase from Arabis blepharophylla was constructed and verified to be stereochemically valid.     

Adaptive evolution in the rat olfactory receptor gene family

Summary Comparison of DNA sequences of the rat (Rattus norvegicus) olfactory receptor gene family revealed an unusual pattern of nucleotide substitution in the gene region encoding the second extracellular domain (E2) of the protein. In this domain, nonsynonymous nucleotide differences between members of this subfamily that caused a change in amino acid residue polarity were over four times more frequent than nonsynonymous differences that did not cause a polarity change. This nonrandom pattern of nucleotide substitution is evidence of past directional selection favoring diversification of the E2 domain among members of this subfamily. This in turn suggests that E2 may play some important role in the functions unique to each member of the olfactory receptor family, and that it may perhaps be an odorant binding domain.Offprint requests to: A.L. Hughes     

Whole chloroplast genome comparison of rice,maize, and wheat: implications for chloroplast gene diversification and phylogeny of cereals

The fully sequenced chloroplast genomes of maize (subfamily Panicoideae), rice (subfamily Bambusoideae), and wheat (subfamily Pooideae) provide the unique opportunity to investigate the evolution of chloroplast genes and genomes in the grass family (Poaceae) by whole-genome comparison. Analyses of nucleotide sequence variations in 106 cereal chloroplast genes with tobacco sequences as the outgroup suggested that (1) most of the genic regions of the chloroplast genomes of maize, rice, and wheat have evolved at similar rates; (2) RNA genes have highly conservative evolutionary rates relative to the other genes; (3) photosynthetic genes have been under strong purifying selection; (4) between the three cereals, 14 genes which account for about 28% of the genic region have evolved with heterogeneous nucleotide substitution rates; and (5) rice genes tend to have evolved more slowly than the others at loci where rate heterogeneity exists. Although the mechanism that underlies chloroplast gene diversification is complex, our analyses identified variation in nonsynonymous substitution rates as a genetic force that generates heterogeneity, which is evidence of selection in chloroplast gene diversification at the intrafamilial level. Phylogenetic trees constructed with the variable nucleotide sites of the chloroplast genes place maize basal to the rice-wheat clade, revealing a close relationship between the Bambusoideae and Pooideae.     

Identification of genes with fast-evolving regions in microbial genomes

Complete sequences of multiple strains of the same microbial species provide an invaluable source for studying the evolutionary dynamics between orthologous genes over a relatively short time scale. Usually the intensity of the selection pressure is inferred from a comparison between the nonsynonymous substitution rate and the synonymous substitution rate. In this paper, we propose an alternative method for detecting genes with one or more fast-evolving regions from pairwise comparisons of orthologous genes. Our method looks for regions with overrepresented nonsynonymous mutations along the alignment, and requires a higher nonsynonymous evolution rate in those regions than the neutral evolution rate. It identifies gene targets under intensive selection pressure that are not detected from the conventional rate comparison analysis. For those identified genes with known annotations, most of them have a clear role in processes such as bacterial defense and host–pathogen interactions. Gene sets reported from our method provide a measure of the phenotypic divergence between two closely related genomes.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Molecular evolution in the hypervariable regions of fetuin: comparison between human and African green monkey fetuin

Sequences of fetuin cDNA and its deduced amino acid residues from the African green monkey cell line Vero were found to differ by 7.3% and 12.9%, respectively, from the corresponding human sequences. Most amino acid substitutions were clustered within a small segment of the third domain (D3). Calculations of nonsynonymous and synonymous nucleotide substitution rates suggest that this small segment was mutated under positive selection. cDNAs encoding alpha1-antitrypsin, beta-actin and the sequences of intron 4 of alpha1-antitrypsin gene in human liver and Vero cells were also investigated. The results substantiated the positive selection imposed on the D3 segment.     

Runaway evolution of the male-specific exon of the doublesex gene in Diptera

In Diptera (Insecta), alternatively spliced male-specific and female-specific products of the doublesex (dsx) gene play a key role in regulating development of the adult genital structures from the genital disc. Analysis of the pattern of nucleotide substitution of different domains of the dsx gene in 29 dipteran species showed that, over short evolutionary times, purifying selection predominated on the domain common to both sexes, the female-specific exons, and the and male-specific exon. However, over longer the evolutionary time frames represented by between-family comparisons, the male-specific exon accumulated nonsynonymous substitutions at a much more rapid rate than either the common domain or the female-specific exon. Overall, the accumulation of nonsynonymous substitutions in the male-specific exon occurred at a significantly greater than linear rate relative to the common domain, whereas the accumulation of nonsynonymous substitutions in the female-specific exon occurred at less than linear rate relative to the common domain. The evolution of the male-specific exon of dsx thus shows a pattern reminiscent of that seen in the "runaway" evolution of male secondary sexual characters at the morphological level, consistent with the hypothesis that female choice is an important factor in the morphological diversification of insect male genitalia.     

Evolutionary lineages of RT1.Ba in the Australian Rattus

In this study, the evolutionary history of the variable second exon of RT1.Ba and its adjoining intron b are compared across a number of species and subspecies of the Australian RATTUS: Three lineages are identified in the second intron across a range of Rattus species. Two of these lineages, separated by the insertion of a probable rodent short interspersed nucleotide element and by point mutations outside the indel region, are both found in each of the major clades of the endemic Australian RATTUS: This pattern of ancestral polymorphism is reflected in the adjoining exon 2 sequences, although phylogenetic constraints confirm that the clustering is not identical to that of the associated intron sequences. In addition, the coding sequences show evidence of the retention of ancestral polymorphism, with identical exon sequences found in two divergent species, and some indication of gene conversion detected for the exon sequences.     

The alpha-mannosidases: phylogeny and adaptive diversification

alpha-Mannosidase enzymes comprise a class of gylcoside hydrolases involved in the maturation and degradaton of glycoprotein-linked oligosaccharides. Various alpha-mannosidase enzymatic activities are encoded by an ancient and ubiquitous gene superfamily. A comparative sequence analysis was employed here to characterize the evolutionary relationships and dynamics of the alpha-mannosidase superfamily. A series of lineage-specific BLAST searches recovered the first ever recognized archaean and eubacterial alpha-mannosidase sequences, in addition to numerous eukaryotic sequences. Motif-based alignment and subsequent phylogenetic analysis of the entire superfamily revealed the presence of three well-supported monophyletic clades that represent discrete alpha-mannosidase families. The comparative method was used to evaluate the phylogenetic distribution of alpha-mannosidase functional variants within families. Results of this analysis demonstrate a pattern of functional diversification of alpha-mannosidase paralogs followed by conservation of function among orthologs. Nucleotide polymorphism among the most closely related pair of duplicated genes was analyzed to evaluate the role of natural selection in the functional diversification of alpha-mannosidase paralogs. Ratios of nonsynonymous and synonymous variation show an increase in the rate of nonsynonymous change after duplication and a relative excess of fixed nonsynonymous changes between the two groups of paralogs. These data point to a possible role for positive Darwinian selection in the evolution of alpha-mannosidase functional diversification following gene duplication.     

The Molecular Evolution of the Small Heat-Shock Proteins in Plants

The small heat-shock proteins have undergone a tremendous diversification in plants; whereas only a single small heat-shock protein is found in fungi and many animals, over 20 different small heat-shock proteins are found in higher plants. The small heat-shock proteins in plants have diversified in both sequence and cellular localization and are encoded by at least five gene families. In this study, 44 small heat-shock protein DNA and amino acid sequences were examined, using both phylogenetic analysis and analysis of nucleotide substitution patterns to elucidate the evolutionary history of the small heat-shock proteins. The phylogenetic relationships of the small heat-shock proteins, estimated using parsimony and distance methods, reveal that gene duplication, sequence divergence and gene conversion have all played a role in the evolution of the small heat-shock proteins. Analysis of nonsynonymous substitutions and conservative and radical replacement substitutions (in relation to hydrophobicity) indicates that the small heat-shock protein gene families are evolving at different rates. This suggests that the small heat-shock proteins may have diversified in function as well as in sequence and cellular localization.     

Evolution of alcohol dehydrogenase genes in peonies (Paeonia): phylogenetic relationships of putative nonhybrid species

Alcohol dehydrogenase genes were amplified by PCR, cloned, and sequenced from 11 putative nonhybrid species of the angiosperm genus Paeonia. Sequences of five exons and six intron regions of the Adh gene were used to reconstruct the phylogeny of these species. Two paralogous genes, Adh1A, and Adh2, were found; an additional gene, Adh1B, is also present in section Moutan. Phylogenetic analyses of exon sequences of the Adh genes of Paeonia and a variety of other angiosperms imply that duplication of Adh1 and Adh2 occurred prior to the divergence of Paeonia species and was followed by a duplication resulting in Adh1A and Adh1B. Concerted evolution appears to be absent between these paralogous loci. Phylogenetic analysis of only the Paeonia Adh exon sequences, positioning the root of the tree between the paralogous genes Adh1 and Adh2, suggests that the first evolutionary split within the genus occurred between the shrubby section Moutan and the other two herbaceous sections Oneapia and Paeonia. Restriction of Adh1B genes to section Moutan may have resulted from deletion of Adh1B from the common ancestor of sections Oneapia and Paeonia. A relative-rate test was designed to compare rates of molecular change among lineages based on the divergence of paralogous genes, and the results indicate a slower rate of evolution within the shrubby section Moutan than in section Oneapia. This may be responsible for the relatively long branch length of section Oneapia and the short branch length between section Moutan and the other two sections found on the Adh, ITS (nrDNA), and matK (cpDNA) phylogenies of the genus. Adh1 and Adh2 intron sequences cannot be aligned, and we therefore carried out separate analyses of Adh1A and Adh2 genes using exon and intron sequences together. The Templeton test suggested that there is not significant incongruence among Adh1A, ITS, and matK data sets, but that these three data sets conflict significantly with Adh2 sequence data. A combined analysis of Adh1A, ITS, and matK sequences produced a tree that is better resolved than that of any individual gene, and congruent with morphology and the results of artificial hybridization. It is therefore considered to be the current best estimate of the species phylogeny. Paraphyly of section Paeonia in the Adh2 gene tree may be caused by longer coalescence times and random sorting of ancestral alleles.      

The Evolutionary Dynamics of Human Influenza B Virus

Despite their close phylogenetic relationship, type A and B influenza viruses exhibit major epidemiological differences in humans, with the latter both less common and less often associated with severe disease. However, it is unclear what processes determine the evolutionary dynamics of influenza B virus, and how influenza viruses A and B interact at the evolutionary scale. To address these questions we inferred the phylogenetic history of human influenza B virus using complete genome sequences for which the date (day) of isolation was available. By comparing the phylogenetic patterns of all eight viral segments we determined the occurrence of segment reassortment over a 30-year sampling period. An analysis of rates of nucleotide substitution and selection pressures revealed sporadic occurrences of adaptive evolution, most notably in the viral hemagglutinin and compatible with the action of antigenic drift, yet lower rates of overall and nonsynonymous nucleotide substitution compared to influenza A virus. Overall, these results led us to propose a model in which evolutionary changes within and between the antigenically distinct 'Yam88' and 'Vic87' lineages of influenza B virus are the result of changes in herd immunity, with reassortment continuously generating novel genetic variation. Additionally, we suggest that the interaction with influenza A virus may be central in shaping the evolutionary dynamics of influenza B virus, facilitating the shift of dominance between the Vic87 and the Yam88 lineages.     

A deletion adjacent to the maize transposable element Mu-1 accompanies loss of Adh1 expression.

Insertion of the maize transposable element Mu-1 into the first intron of the alcohol dehydrogenase locus (Adh1) of maize produced mutant Adh1-S3034 with 40% of the wild-type level of protein and mRNA. Continued instability at this locus resulted in secondary mutations with lower levels of protein expression. One of these, Adh1-S3034a, has no detectable ADH1 expression. This paper describes the precise nature of the changes in the Adh1 gene that gave rise to the S3034a allele. The Mu-1 element is still present in the mutant, but Adh1 sequences immediately adjacent to the element are deleted. The deletion starts precisely at the Mu-1 insertion site and extends 74 bp leftward removing part of the first intron, the intron:exon junction and 2 bp of the eleventh amino acid codon in the first exon of the gene. Tests for reversion within the somatic tissue of plants show that mutant S3034a, unlike its progenitor, is stably null for ADH1 activity.     

Evolution of immunoglobulin VH pseudogenes in chickens

In chickens, there is a single functional gene (VH1) coding for the heavy chain variable region of immunoglobulins, and immunoglobulin diversity is generated by gene conversion of the VH1 gene by many variable region pseudogenes (psi VH's) that exist on the 5' side of the VH1 gene. To understand the evolution of this unique genetic system, we conducted statistical analyses of VH1 and psi VH genes together with functional VH genes from other higher vertebrate species. The results indicate, first, that chicken VH genes are all closely related to one another and were derived relatively recently from an ancestral gene belonging to one of the three major groups of VH genes in higher vertebrates. Second, the rate of nonsynonymous substitution is slightly higher than that of synonymous substitution in the complementarity- determining regions (CDRs), which suggests that diversity-enhancing selection has operated in the CDRs even for pseudogenes. However, both the rates of synonymous and nonsynonymous substitution are higher in the CDRs than in the framework regions (FRs), apparently because of an interaction between positive selection and meiotic gene conversion in the CDRs. Third, a dot matrix analysis of the psi VH genes and genomic diversity (D) genes has indicated that the 3' end of psi VH genes is attached by D-gene-like sequences, and this region of psi VH genes has high similarity with D gene sequences. This suggests that V and D genes were fused at some point of evolutionary time and this fused element multiplied by gene duplication. Finally, two alternative hypotheses of explaining the evolution of the chicken VH gene system are presented.      

Evolutionary rates and expression level in Chlamydomonas

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.     

Analysis of substitution events in HIV-1 vif gene of the Korean clade

Nucleotide and amino acid substitution pattern in vif gene of the Korean clade of HIV-1 isolated from Koreans were analyzed using consensus sequences. At nucleotide level, transition/transversion substitution ratio was 1.88, and nonsynonymous/synonymous substitution ratio was 2.67, suggesting a divergent evolution in the Korean clade. At amino acid level, there were 17 substitutions and G-->E substitution at position 37 may be responsible for change in predicted secondary structure.     

Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

蕨类植物rpoC1内含子缺失及其分子进化速率

rpoC1基因编码RNA聚合酶β°亚基蛋白, 在转录过程中与DNA模板结合, 与β亚基形成的β-β°亚基复合体构成RNA合成的催化中心。以rpoC1基因为研究对象, 在贝叶斯因子大于20的条件下, 用HyPhy软件位点模型检测到3个正选择位点和541个负选择位点; 用PAML软件位点模型检测到10个正选择位点, 其中3个位点的后验概率超过99%。此外, 基于最大似然法构建64种蕨类植物的系统发育树, 结合HyPhy软件分析rpoC1基因的转换率、颠换率、转换率/颠换率、同义替换率、非同义替换率以及同义替换率/非同义替换率, 探讨rpoC1基因内含子丢失与分子进化速率的关系。结果表明, rpoC1基因内含子缺失对转换率、颠换率以及非同义替换率有一定影响。