AtHVA22 gene family in Arabidopsis: phylogenetic relationship,ABA and stress regulation,and tissue-specific expression

HVA22 is an ABA- and stress-inducible gene first isolated from barley (Hordeum vulgare L.). Homologues of HVA22 have been found in plants, animals, fungi and protozoa, but not in prokaryotes, suggesting that HVA22 plays a unique role in eukaryotes. Five HVA22 homologues, designated AtHVA22a, b, c, d and e, have been identified in Arabidopsis. These five AtHVA22 homologues can be separated into two subfamilies, with AtHVA22a, b and c grouped in one subfamily and AtHVA22d and e in the other. Phylogenetic analyses show that AtHVA22d and e are closer to barley HVA22 than to AtHVA22a, b and c, suggesting that the two subfamilies had diverged before the divergence of monocots and dicots. The distribution and size of exons of AtHVA22 homologues and barley HVA22 are similar, suggesting that these genes are descendents of a common ancestor. AtHVA22 homologues are differentially regulated by ABA, cold, dehydration and salt stresses. These four treatments enhance AtHVA22a, d and e expression, but have little or even suppressive effect on AtHVA22c expression. ABA and salt stress induce AtHVA22b expression, but cold stress suppresses ABA induction of this gene. Expression of AtHVA22d is the most tightly regulated by these four treatments among the five homologues. In general, AtHVA22 homologues are expressed at a higher level in flower buds and inflorescence stems than in rosette and cauline leaves. The expression level of these homologues in immature siliques is the lowest among all tissues analyzed. It is suggested that some of these AtHVA22 family members may play a role in stress tolerance, and others are involved in plant reproductive development.     

Gene conversion as a source of nucleotide diversity in Plasmodium falciparum

Examination of polymorphisms in the Plasmodium falciparum gene for falcipain 2 revealed that this gene is one of two paralogs separated by 10.8 kb in chromosome 11. We designate the annotated gene denoted chr11.gen_424 as encoding falcipain 2A and the annotated gene denoted chr11.gen_427 as encoding falcipain 2B. The paralogs are 96% identical at the nucleotide level and 93% identical at the amino acid level. The consensus sequences differ in 31/309 synonymous sites and 45/1140 nonsynonymous sites, including three amino acid replacements (V393I, A400P, and Q414E) that are near the catalytic site and that may affect substrate affinity or specificity. In six reference isolates, among 36 synonymous sites and 46 nonsynonymous sites that are polymorphic in the gene for falcipain 2A, falcipain 2B, or both, significant spatial clustering is observed. All but one of the polymorphisms appear to result from gene conversion between the paralogs. The estimated rate of gene conversion between the paralogs may be as many as 1,400 to 1,700 times greater than the rate of mutation. Owing to gene conversion, one of the falcipain 2A alleles is more similar to the falcipain 2B alleles than it is to other falcipain 2A alleles. Divergence among the synonymous sites suggests that the paralogous genes last shared a common ancestor 15.2 MYA, with a range of 8.8 to 20.6 MYA. During this period, the paralogs have acquired 0.10 synonymous substitutions per synonymous site in the coding region. The 5' and 3' flanking regions differ in 47.7% and 39.8% of the nucleotide sites, respectively. Hence synonymous sites and flanking regions are not conserved in sequence in spite of their high AT content and T skew.     

Micro-scale signature of purifying selection in Marburg virus genomes

In the seven protein-coding genes in the Marburg virus (MARV) genome, the synonymous nucleotide diversity substantially exceeded the nonsynonymous nucleotide diversity, indicating strong purifying selection. Likewise, there was evidence of purifying selection on 5'UTR and 3'UTR, where nucleotide diversity (pi) was significantly less than piS in the coding regions. Nonsynonymous polymorphic sites showed significantly reduced mean gene diversity in comparison to other polymorphic sites, indicating that purifying selection at certain slightly deleterious nonsynonymous polymorphisms is ongoing. Moreover, nonsynonymous polymorphic sites showed significantly reduced gene diversity in comparison to adjacent synonymous sites, even though the vast majority of such adjacent synonymous sites were in the same codon or an adjacent codon. Thus purifying selection, in conjunction with recombination and/or backward mutation, can act to break up linkage relationships at a micro-scale in the MARV genome. The ability of purifying selection to break up linkage between synonymous and nonsynonymous polymorphisms on such a fine scale has not been reported in any other genome.     

AtHVA22 gene family in Arabidopsis: phylogenetic relationship,ABA and stress regulation,and tissue-specific expression

HVA22 is an ABA- and stress-inducible gene first isolated from barley (Hordeum vulgare L.). Homologues of HVA22 have been found in plants, animals, fungi and protozoa, but not in prokaryotes, suggesting that HVA22 plays a unique role in eukaryotes. Five HVA22 homologues, designated AtHVA22a, b, c, d and e, have been identified in Arabidopsis. These five AtHVA22 homologues can be separated into two subfamilies, with AtHVA22a, b and c grouped in one subfamily and AtHVA22d and e in the other. Phylogenetic analyses show that AtHVA22d and e are closer to barley HVA22 than to AtHVA22a, bandc, suggesting that the two subfamilies had diverged before the divergence of monocots and dicots. The distribution and size of exons of AtHVA22 homologues and barley HVA22 are similar, suggesting that these genes are descendents of a common ancestor. AtHVA22 homologues are differentially regulated by ABA, cold, dehydration and salt stresses. These four treatments enhance AtHVA22a, d and e expression, but have little or even suppressive effect on AtHVA22c expression. ABA and salt stress induce AtHVA22b expression, but cold stress suppresses ABA induction of this gene. Expression of AtHVA22d is the most tightly regulated by these four treatments among the five homologues. In general, AtHVA22 homologues are expressed at a higher level in flower buds and inflorescence stems than in rosette and cauline leaves. The expression level of these homologues in immature siliques is the lowest among all tissues analyzed. It is suggested that some of these AtHVA22 family members may play a role in stress tolerance, and others are involved in plant reproductive development.     

Genetic diversity and differentiation among Lagopus lagopus populations in Scandinavia and Scotland: evolutionary significant units confirmed by SNP markers

Single Nucleotide Polymorphism in four Scandinavian populations of willow grouse (Lagopus lagopus) and two Scottish populations of red grouse (Lagopus lagopus scoticus) were assessed at 13 protein‐coding loci. We found high levels of diversity, with one substitution every 55 bp as an average and a total of 76 unlinked parsimony informative SNPs. Different estimators of genetic diversity such as: number of synonymous and non‐synonymous sites, average number of alleles, number and percentage of polymorphic loci, mean nucleotide diversity (π_s, π_a) and gene diversity at synonymous and non‐synonymous sites showed higher diversity in the northern populations compared to southern ones. Strong levels of purifying selection found in all the populations together with neutrality tests conforming to neutral expectations agree with large effective population sizes. Assignment tests reported a clear distinction between Scandinavian and Scottish grouse suggesting the existence of two different evolutionary significant units. The divergence time between willow and red grouse ranging between 12 500 and 125 000 years, in conjunction with the presence of ‘specific’ markers for each subspecies prompt a reassessment of the taxonomical status of the Scottish red grouse.     

Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

Low Nucleotide Diversity in Man

The nucleotide diversity (pi) in humans is studied by using published cDNA and genomic sequences that have been carefully checked for sequencing accuracy. This measure of genetic variability is defined as the number of nucleotide differences per site between two randomly chosen sequences from a population. A total of more than 75,000 base pairs from 49 loci are compared. The DNA regions studied are the 5' and 3' untranslated regions and the amino acid coding regions. The coding regions are divided into nondegenerate sites (i.e., sites at which all possible changes are nonsynonymous), twofold degenerate sites (i.e., sites at each of which one of the three possible changes is synonymous) and fourfold degenerate sites (i.e., sites at which all three possible changes are synonymous). The pi values estimated are, respectively, 0.03 and 0.04% for the 5' and 3' UT regions, and 0.03, 0.06 and 0.11% for nondegenerate, twofold degenerate and fourfold degenerate sites. Since the highest pi value is only 0.11%, which is about one order of magnitude lower than those in Drosophila populations, the nucleotide diversity in humans is very low. The low diversity is probably due to a relatively small long-term effective population size rather than any severe bottleneck during human evolution.     

A conservative test of genetic drift in the endosymbiotic bacterium Buchnera: slightly deleterious mutations in the chaperonin groEL

The obligate endosymbiotic bacterium Buchnera aphidicola shows elevated rates of sequence evolution compared to free-living relatives, particularly at nonsynonymous sites. Because Buchnera experiences population bottlenecks during transmission to the offspring of its aphid host, it is hypothesized that genetic drift and the accumulation of slightly deleterious mutations can explain this rate increase. Recent studies of intraspecific variation in Buchnera reveal patterns consistent with this hypothesis. In this study, we examine inter- and intraspecific nucleotide variation in groEL, a highly conserved chaperonin gene that is constitutively overexpressed in Buchnera. Maximum-likelihood estimates of nonsynonymous substitution rates across Buchnera species are strikingly low at groEL compared to other loci. Despite this evidence for strong purifying selection on groEL, our intraspecific analysis of this gene documents reduced synonymous polymorphism, elevated nonsynonymous polymorphism, and an excess of rare alleles relative to the neutral expectation, as found in recent studies of other Buchnera loci. Comparisons with Escherichia coli generally show patterns predicted by their differences in N(e). The sum of these observations is not expected under relaxed or balancing selection, selective sweeps, or increased mutation rate. Rather, they further support the hypothesis that drift is an important force driving accelerated protein evolution in this obligate mutualist.     

On the effects of background selection in small populations on comparisons of molecular variation

Deleterious mutations affect genetic variation at linked neutral loci. Neutral variation can be reduced due to background selection, but in small population and with tight linkage such variation may increase due to associative overdominance. Here I report the results of computer simulations of diploid genotypes in small populations, where I look at the effect of deleterious mutations and linkage on comparisons of intra- and interspecific variation. Each chromosome consisted of 2000 loci where deleterious and neutral mutations occurred. The ratio of nonsynonymous to synonymous substitution rates (Ka/Ks) either increases with tight linkage or is unaffected, depending on the strength of selection. The ratio of the numbers of segregating mutations to the number of fixed mutations decreases under the conditions leading to background selection but can increase at tight linkage. Numbers of segregating sites (Sn) are less affected than nucleotide site diversity (pi), pi reduces more than Sn at intermediate linkage, but pi increases more than Sn when linkage is tight. Similar effects as found for Sn and pi are observed for heterozygosity and variance in allele size of tandem repeat loci.     

The PI/GLO-like locus in orchids: duplication and purifying selection at synonymous sites within Orchidinae (Orchidaceae)

Positive selection and relaxation of purifying constraints after duplication events have driven the functional diversification of gene families involved in development. One example of this occurred within the plant MADS-box genes. The evolution of the orchid flower was driven by duplication events followed by sub- and neo-functionalization of class B DEF-like MADS-box genes, which are present at three to four copies in the orchid genome. In contrast, the orchid PI/GLO-like class B MADS-box genes have been reported thus far as single-copy loci, with the only exception of Habenaria radiata.We isolated a novel PI/GLO-like gene (OrcPI2) in Orchis italica, which is different than the previously characterized OrcPI locus. The presence of two functional paralogs of PI/GLO-like genes in orchids is detectable only within the tribe Orchidinae. Evolutionary analyses revealed an apparent relaxation of purifying selection acting on the two PI/GLO-like paralogs of the Orchidinae when compared to the single-copy PI/GLO-like genes found in other orchid species. Furthermore, by measuring dN/dS (ω) ratios, we show that a high percentage of sites between the two PI/GLO-like paralogs have different evolutionary pressures. Interestingly, the apparent relaxation of selective constraints on the two PI/GLO-like paralogs is due to strong purifying selection at synonymous sites rather than to a high value of nonsynonymous substitution rate. This peculiar evolutionary pattern might be related to molecular processes such as mRNA folding and/or translational efficiency control. These processes could potentially be involved in or predate the functional diversification of the two PI/GLO-like paralogs within Orchidinae.     

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.     

Evolution of duplicate genes in a tetraploid animal, Xenopus laevis

To understand the evolution of duplicate genes, we compared rates of nucleotide substitution between 17 pairs of nonallelic duplicated genes in the tetraploid frog Xenopus laevis with rates between the orthologous loci of human and rodent. For all duplicated X. laevis genes, the number of synonymous substitutions per site (dS) was greater than the number of nonsynonymous substitutions per site (dN), indicating that these genes are subject to purifying selection. There was also a significant positive correlation (r = 0.915) between dN for the X. laevis genes and dN for the mammalian genes, suggesting that, at the amino acid level, the X. laevis genes and the mammalian genes are under similar constraints. Results of relative-rate tests showed nearly equal rates of nonsynonymous substitution in each copy of the X. laevis genes; apparently there are similar constraints on both copies. No correlation was found between dS for the X. laevis genes and dS for the mammalian genes. There was a significant positive correlation both between members of pairs of duplicated X. laevis genes (r = 0.951) and between human and rodent orthologues (r = 0.854) with respect to third- position G+C content but no such relationship between the X. laevis genes and either of their mammalian orthologues. The results indicate that both copies of a duplicate gene can be subject to purifying selection and thus support the hypothesis of selection against all genotypes containing a null allele at either of two duplicate loci.      

A maximum likelihood method for analyzing pseudogene evolution: implications for silent site evolution in humans and rodents.

We present a new likelihood method for detecting constrained evolution at synonymous sites and other forms of nonneutral evolution in putative pseudogenes. The model is applicable whenever the DNA sequence is available from a protein-coding functional gene, a pseudogene derived from the protein-coding gene, and an orthologous functional copy of the gene. Two nested likelihood ratio tests are developed to test the hypotheses that (1) the putative pseudogene has equal rates of silent and replacement substitutions; and (2) the rate of synonymous substitution in the functional gene equals the rate of substitution in the pseudogene. The method is applied to a data set containing 74 human processed-pseudogene loci, 25 mouse processed-pseudogene loci, and 22 rat processed-pseudogene loci. Using the informatics resources of the Human Genome Project, we localized 67 of the human-pseudogene pairs in the genome and estimated the GC content of a large surrounding genomic region for each. We find that, for pseudogenes deposited in GC regions similar to those of their paralogs, the assumption of equal rates of silent and replacement site evolution in the pseudogene is upheld; in these cases, the rate of silent site evolution in the functional genes is approximately 70% the rate of evolution in the pseudogene. On the other hand, for pseudogenes located in genomic regions of much lower GC than their functional gene, we see a sharp increase in the rate of silent site substitutions, leading to a large rate of rejection for the pseudogene equality likelihood ratio test.     

Nucleotide sequence diversity at the methionine synthase locus in endangered Dunnia sinensis (Rubiaceae): an evaluation of the positive selection hypothesis

Methionine synthase is a key enzyme for the synthesis of the aspartate-derived methionine, the immediate precursor of S-adenosyl-methionine, which has been illustrated to be associated with plant growth and pathogen interactions. In this study we tested the positive selection hypothesis of molecular evolution of the methionine synthase gene in Dunnia sinensis. In the entire sample of 87 sequences, 22 haplotypes of introns and 16 haplotypes of exons were identified. An excess of polymorphism over the neutral expectation for the class of unique nucleotide polymorphisms was observed in both exon and intron sequences. Ten replacement substitutions versus six synonymous substitutions among lineages, although nonsignificant, revealed that some advantageous mutants might have been favored. The distribution of d(N)/d(S) > 1 at nodes between closely related haplotypes in the gene network also indicated weak and variable positive selection. Nevertheless, low levels of genetic diversity in exons (theta; = 0.0052) and introns (theta; = 0.0070) of the methionine synthase gene of the outcrossing Dunnia were also attributed to the endangered status of the species. The atpB-rbcL intergenic spacer of cpDNA and the ribosomal internal transcribed spacer of mtDNA were used to discern the relative effectiveness of natural selection from intrinsic evolutionary forces. The low levels of nucleotide polymorphisms in both organelle spacers and the significant population differentiation reflected the effect of population-species history and demography. Two major lineages of the methionine synthase gene genealogy were recovered corresponding to two geographic regions, a result that was consistent with organelle phylogenies. Both past fragmentation and recent habitat disturbance causing complete bottlenecks may have resulted in population decline and geographic isolation and may have led to the depletion of genetic variation at loci in nuclear and organelle genomes.     

Extensive concerted evolution of rice paralogs and the road to regaining independence

Many genes duplicated by whole-genome duplications (WGDs) are more similar to one another than expected. We investigated whether concerted evolution through conversion and crossing over, well-known to affect tandem gene clusters, also affects dispersed paralogs. Genome sequences for two Oryza subspecies reveal appreciable gene conversion in the approximately 0.4 MY since their divergence, with a gradual progression toward independent evolution of older paralogs. Since divergence from subspecies indica, approximately 8% of japonica paralogs produced 5-7 MYA on chromosomes 11 and 12 have been affected by gene conversion and several reciprocal exchanges of chromosomal segments, while approximately 70-MY-old "paleologs" resulting from a genome duplication (GD) show much less conversion. Sequence similarity analysis in proximal gene clusters also suggests more conversion between younger paralogs. About 8% of paleologs may have been converted since rice-sorghum divergence approximately 41 MYA. Domain-encoding sequences are more frequently converted than nondomain sequences, suggesting a sort of circularity--that sequences conserved by selection may be further conserved by relatively frequent conversion. The higher level of concerted evolution in the 5-7 MY-old segmental duplication may reflect the behavior of many genomes within the first few million years after duplication or polyploidization.     

Contrasting selection modes at the Adh1 locus in outcrossing Miscanthus sinensis vs. inbreeding Miscanthus condensatus (Poaceae)

We estimated DNA sequence variation of the Adh1 locus in the outcrossing Miscanthus sinensis (Poaceae) and its close selfing relative, M. condensatus. Tajima's test of selection is significantly negative for both overall exons and replacement sites in M. sinensis. Among its entire sample, nucleotide diversity of nonsynonymous sites is higher than that of synonymous sites. A McDonald and Kreitman test of neutrality indicates an excess of intraspecific replacement polymorphisms, suggesting possible directional selection toward advantageous mutants. However, frequent intragenic recombination suggests both purifying and positive selection is unlikely. Recent demographic expansions coupled with relaxation of purifying selection may have resulted in elevated genetic diversity at the Adh1 locus as well as the trnL-trnF intergenic spacer of cpDNA in this outcrossing species. In contrast, low levels of genetic diversity were detected at both the Adh1 locus and the cpDNA spacer in M. condensatus, consistent with bottlenecks associated with selfing in all populations. While Tajima's D and Fu and Li's F statistics did not reveal deviation from neutrality at the Adh1 locus in M. condensatus, 12 replacements vs. 10 synonymous changes were detected. Based on pairwise comparisons of the d(N)/d(S) ratio, lineages of closely related populations of the species distributed along saline habitats appeared to be under directional selection.     

RECENT GENE‐CAPTURE ON THE UV SEX CHROMOSOMES OF THE MOSS CERATODON PURPUREUS

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex‐linked loci. In this system, recombination is suppressed on both the female‐transmitted (U) sex chromosome and the male‐transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex‐limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex‐linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (～0.6–1.3 MYA and ～2.8–3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U‐ and V‐linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.     

Molecular evolution of the phytochrome gene family in sorghum: changing rates of synonymous and replacement evolution

The photoreceptor phytochromes, encoded by a small gene family, are responsible for controlling the expression of a number of light-responsive genes and photomorphogenic events, including agronomically important phenotypes such as flowering time and shade-avoidance behavior. The understanding and control of flowering time are particularly important goals in sorghum cultivar development for diverse environments, and naturally occurring variation in the phytochrome genes might prove useful in breeding programs. Also of interest is whether variation observed at the phytochrome loci in domesticated sorghum, or in particular races, is a result of human selection. Population genetic studies can reveal evidence of such selection in patterns of polymorphism and divergence. In this study we report a population genetic analysis of the PHY gene family in Sorghum bicolor (L.) Moench in a diverse panel including both cultivated and wild accessions. We show that the level of nucleotide variation in all gene family members is about half the average for this species, consistent with purifying selection acting on these loci. However, the rate of amino acid substitution is accelerated at PHYC compared to the other two loci. In comparisons to a closely related sorghum species, PHYC shows a pattern of intermediate frequency amino acid changes that differ from the patterns observed in comparisons across longer evolutionary distances. There is also a departure from expected patterns of polymorphism and divergence at synonymous sites in PHYC, although the data do not fit a simple model of directional or diversifying selection. Cultivated sorghum has a level of variation similar to that of wild relatives (ssp. verticilliflorum), but many polymorphisms are subspecies-specific, including several amino acid variants.