Genome comparisons reveal a dominant mechanism of chromosome number reduction in grasses and accelerated genome evolution in Triticeae.

Single-nucleotide polymorphism was used in the construction of an expressed sequence tag map of Aegilops tauschii, the diploid source of the wheat D genome. Comparisons of the map with the rice and sorghum genome sequences revealed 50 inversions and translocations; 2, 8, and 40 were assigned respectively to the rice, sorghum, and Ae. tauschii lineages, showing greatly accelerated genome evolution in the large Triticeae genomes. The reduction of the basic chromosome number from 12 to 7 in the Triticeae has taken place by a process during which an entire chromosome is inserted by its telomeres into a break in the centromeric region of another chromosome. The original centromere-telomere polarity of the chromosome arms is maintained in the new chromosome. An intrachromosomal telomere-telomere fusion resulting in a pericentric translocation of a chromosome segment or an entire arm accompanied or preceded the chromosome insertion in some instances. Insertional dysploidy has been recorded in three grass subfamilies and appears to be the dominant mechanism of basic chromosome number reduction in grasses. A total of 64% and 66% of Ae. tauschii genes were syntenic with sorghum and rice genes, respectively. Synteny was reduced in the vicinity of the termini of modern Ae. tauschii chromosomes but not in the vicinity of the ancient termini embedded in the Ae. tauschii chromosomes, suggesting that the dependence of synteny erosion on gene location along the centromere-telomere axis either evolved recently in the Triticeae phylogenetic lineage or its evolution was recently accelerated.

Viability versus Fecundity Selection in the Slender Wild Oat, Avena barbata L.

Estimates of relative selective values of different genotypes in natural populations of predominantly self-pollinated plants have consistently shown a net heterozygote advantage. Heterozygote advantage is further analyzed in the present study by estimating components of selective values corresponding approximately to viability and fecundity. The results show that a higher proportion of heterozygotes than of homozygotes survive from seeding to adult stages and that this viability selection accounts for the overall excess of heterozygotes.

Monte carlo studies of plant mating system estimation models: the one-pollen parent and mixed mating models.

Estimation of mating system parameters in plant populations typically employs family-structured samples of progeny genotypes. These estimation models postulate a mixture of self-fertilization and random outcrossing. One assumption of such models concerns the distribution of pollen genotypes among eggs within single maternal families. Previous applications of the mixed mating model to mating system estimation have assumed that pollen genotypes are sampled randomly from the total population in forming outcrossed progeny within families. In contrast, the one-pollen parent model assumes that outcrossed progeny within a family share a single-pollen parent genotype. Monte Carlo simulations of family-structured sampling were carried out to examine the consequences of violations of the different assumptions of the two models regarding the distribution of pollen genotypes among eggs. When these assumptions are violated, estimates of mating system parameters may be significantly different from their true values and may exhibit distributions which depart from normality. Monte Carlo methods were also used to examine the utility of the bootstrap resampling algorithm for estimating the variances of mating system parameters. The bootstrap method gives variance estimates that approximate empirically determined values. When applied to data from two plant populations which differ in pollen genotype distributions within families, the two estimation procedures exhibit the same behavior as that seen with the simulated data.

Nucleotide polymorphism in the Adh1 locus of pearl millet (Pennisetum glaucum) (Poaceae).

We investigated nucleotide polymorphism in the Adh1 locus of pearl millet (Pennisetum glaucum) (Poaceae) by determining the DNA sequence of 20 alleles from 10 individuals. The individuals were sampled from throughout pearl millet's indigenous range and represent both wild and cultivated accessions. Our results indicated that there is little nucleotide polymorphism in the Adh1 locus. Estimates of per site nucleotide polymorphism did not differ significantly between cultivated and wild millet accessions. We compared nucleotide polymorphism in pearl millet Adh1 with nucleotide polymorphism in maize (Zea mays) Adh1 and conclude that the maize Adh1 sample is more polymorphic. Increased polymorphism in maize Adh1 may be attributable, in part, to faster substitution rates in the maize lineage. Analysis suggests that substitution rates in the maize Adh1 lineage are approximately 1.7 times faster than substitution rates in the millet Adh1 lineage.

Temporal Stability of Allozyme Frequencies in a Natural Population of DROSOPHILA MELANOGASTER.

Seasonal patterns of allozyme variation are examined for 12 polymorphic enzyme loci in Drosophila melanogaster. The data derive from a total of 56 samples taken from a natural population in the Summer and Fall of 1978 and 1979. Samples were obtained at approximately five-day intervals and assayed for 6-phosphogluconate dehydrogenase (6Pgd), phosphoglucomutase (Pgm) and glucose-6-phosphate dehydrogenase (G6pd). The remaining nine enzymes were assayed in an average of eight samples per season. None of the loci exhibit regular seasonal cycles of gene-frequency change, although 6Pgd does show significant, but irregular, frequency oscillations. There is also little evidence for gene-frequency differences between years, although 6Pgd is again exceptional in showing significant frequency changes between years. In addition, genotypic frequency distributions are usually consistent with random mating expectations. With the notable exception of 6Pgd, the data give a strong impression of gene-frequency homogeneity within and among years, despite obvious seasonal changes in climate and in the distribution of breeding sites.

Dynamics of the Linkage Disequilibrium Function under Models of Gene-Frequency Hitchhiking.

The dynamic behavior of the linkage disequilibrium (D) between a neutral and a selected locus is analyzed for a variety of deterministic selection models. The time-dependent behavior of D is governed by the gene frequency at the selected locus (p) and by the selection (s) and recombination (r) parameters. Thomson (1977) showed numerically that D may increase under certain initial conditions. We give exact conditions for D to increase in time, which require that the selection intensity exceed the recombination fraction (s > r) and that p be near zero or one. We conclude from this result that gene frequency hitchhiking is most likely to be important when a new favorable mutant enters a population. We also show that, for what can be a wide range of gene frequencies, D will decay at a faster rate than the neutral rate. Consequently, the hitchhiking effect may quickly diminish as the selected gene becomes more common.-The method of analysis allows a complete qualitative description of the dynamics of D as a function of s and r. Two major findings concern the range of gene frequencies at the selected locus for which D either increases over time or decays at a faster rate than under neutrality. For all models considered, the region where D increases (i) first enlarges then shrinks as selection intensifies, and (ii) steadily shrinks as r increases. In contrast, the region of accelerated decay constantly enlarges as the selection intensity increases. This region will either shrink or enlarge as r increases, depending upon the form of selection in force.

Dynamics of correlated genetic systems. IV. Multilocus effects of ethanol stress environments.

Four replicate populations of Drosophila melanogaster, two reared on medium supplemented with ethanol and two reared on standard medium, were electrophoretically monitored for 28 generations. During the first 12 generations, allelic, genotypic and gametic frequencies were determined for eight polymorphic enzymes: GOT, alpha-GPDH, MDH, ADH, TO, E6, Ec and ODH. Samples from generation 18 and 28 were electrophoretically typed for ADH and alpha-GPDH. In addition, samples from generation 27 were analyzed for the presence of inversion heterozygotes. The experimental results showed rapid gene-frequency divergence between control and treatment populations at the Adh locus in a direction consistent with the activity hierarchy of Adh genotypes. Gene-frequency divergence between control and treatment populations also occurred at the alpha-Gpdh locus, although the agreement among replicates appeared to have broken down by generation 28. No differential gene-frequency change occurred at any of the six remaining marker loci. Furthermore, values of linkage disequilibria among all linked pairs of genes were initially small and remained small throughout the course of the experiment. Taking these facts into account, it is argued that the gene-frequency response observed at ADH is most probably caused by selection at the Adh locus. The gene frequency response at alpha-Gpdh can also be be accounted for in terms of the effect of ethanol on energy metabolism, although other explanations cannot be excluded.

Chloroplast DNA variation in pearl millet and related species.

The evolution of specific regions of the chloroplast genome was studied in five grass species in the genus Pennisetum, including pearl millet, and one species from a related genus (Cenchrus). Three different regions of the chloroplast DNA were investigated. The first region included a 12-kilobase pair (kbp) EcoRI fragment containing the 23S, 16S and 5S ribosomal RNA genes, which is part of a larger duplicated region of reverse orientation. The second region was contained in a 21-kbp Sa/I fragment, which spans the short single-copy sequence separating the two reverse repeat structures and which overlaps the duplicated copies of the 12-kbp Eco RI fragment. The third region was a 6-kbp EcoRI fragment located in the large single-copy region of the chloroplast genome. Together these regions account for slightly less than 25% of the chloroplast genome. Each of these DNA fragments was cloned and used as hybridization probes to determine the distribution of homologous DNA fragments generated by various restriction endonuclease digests.-A survey of 12 geographically diverse collections of pearl millet showed no indication of chloroplast DNA sequence polymorphism, despite moderate levels of nuclear-encoded enzyme polymorphism. Interspecific and intergeneric differences were found for restriction endonuclease sites in both the small and the large single-copy regions of the chloroplast genome. The reverse repeat structure showed identical restriction site distributions in all materials surveyed. These results suggest that the reverse repeat region is differentially conserved during the evolution of the chloroplast genome.

The Nature of Nucleotide Sequence Divergence between Barley and Maize Chloroplast DNA.

Analysis of a 2175-base pair (bp) SmaI-HindIII fragment of barley chloroplast DNA revealed that rbcL (the gene for the large subunit of ribulose 1,5-bisphosphate carboxylase) and atpB (the gene for the beta subunit of ATPase) are transcribed divergently and are separated by an untranscribed region of 155-166 bp. The rbcL mRNA has a 320-residue untranslated leader region, whereas the atpB mRNA has a 296- to 309-residue leader region. The sequence of these regions, together with the initial 113 bp of the atpB-coding region and the initial 1279 bp of the rbcL-coding region, is compared with the analogous maize chloroplast DNA sequences. Two classes of nucleotide differences are present, substitutions and insertions/deletions. Nucleotide substitutions show a 1.9-fold bias toward transitions in the rbcL-coding region and a 1.5-fold bias toward transitions in the noncoding region. The level of nucleotide substitutions between the barley and maize sequences is about 0.065/bp. Seventy-one percent of the substitutions in the rbcL-coding region are at the third codon position, and 95% of these are synonymous changes. Insertion/deletion events, which are confined to the noncoding regions, are not randomly distributed in these regions and are often associated with short repeated sequences. The extent of change for the noncoding regions (about 0.093 events/bp) is less than the extent of change at the third codon positions in the rbcL-coding region (about 0.135 events/bp), including insertion/delection events. Limited sequence analysis of the analogous DNA from a wild line ( Hordeum spontaneum) and a primitive Iranian barley (H. vulgare) suggested a low rate of chloroplast DNA evolution. Compared to spinach chloroplast DNA, the barley rbcL-atpB untranslated region is extremely diverged, with only the putative rbcL promoters and ribosome-binding site being extensively conserved.

Molecular evolution of the Ac/Ds transposable-element family in pearl millet and other grasses.

We report an Ac-like sequence from pearl millet (Pennisetum glaucum) and deletion derivative Ac-like sequences from pearl millet and another grass species, Bambusa multiplex. Sequence relationships between the pearl millet and maize Ac elements suggest that Ac/Ds transposable-element family is ancient. Further, the sequence identity between the Bambusa Ac-like sequence and maize Ac implies that the Ac/Ds transposable-element family has been in the grass family since its inception. The Ac-like sequences reported from pearl millet and maize Ac are statistically heterogeneous in pair-wise distance comparisons to each other. Yet, we are unable to discriminate between differential selection or ectopic exchange (recombination and conversion) between nonidentical transposable element homologues, as the cause of the heterogeneity. However, the more extreme heterogeneity exhibited between the previously described pearl millet element and maize Ac seems likely to derive from ectopic exchange between elements with different levels of divergence.

Dynamics of Correlated Genetic Systems. V. Rates of Decay of Linkage Disequilibria in Experimental Populations of DROSOPHILA MELANOGASTER.

The dynamic behavior of four-locus gametic frequency distributions was studied in five replicate cage populations of Drosophila melanogaster for up to 50 generations. The joint frequency distributions were resolved into gene frequencies and various disequilibrium measures. In addition, F statistics for marginal single-locus genotypic frequency distributions were followed through time. The gene frequency, disequilibrium and F statistics were obtained for four chromosome 3 enzyme marker loci [isocitrate dehydrogenase (3-27.1), esterase-6 (3-36.8), phosphoglucomutase (3-43.4) and esterase-C (3-49.0)]. The initial structure of the experimental populations featured random mating proportions, and two complementary gametic types with respect to the marker loci, thus assuring complete pairwise linkage disequilibrium among the markers.--The experimental results indicate: (1) the between-replicate variance in gene frequency varied substantially among loci, with isocitrate dehydrogenase showing the greatest between-replicate variance, and esterase-C the least. (2) The F statistics initially were strongly negative but decayed to the neighborhood of zero for all marker loci except esterase-C. The rate at which the F statistics approached zero varied among the marker loci, indicating substantial differences in the distribution of selective effects along the chromosome. The centromeric region, marked by esterase-C, shows the strongest selective effects. (3) The rate of decay of linkage disequilibrium was much faster than expected for pairs of neutral loci, averaging 1.82 times the neutral rate over all replicates and pairs of loci. This acceleration, which was observed for all six pairwise combinations of loci, was interpreted as resulting from the interaction between selection and recombination. Our experimental results are consistent with many investigations of linkage disequilibrium in natural populations of Drosophila melanogaster that show little or no disequilibrium among enzyme loci. (4) A fortuitous contamination of two cages revealed an apparent regulatory interaction between the migrant and nonmigrant chromosomes at the esterase-C locus. The migrant chromosomes were very rapidly absorbed into the recipient populations, despite this interaction. This result suggests that the dynamics of migration in populations may be phenomenologically richer than anticipated by simple theory.

Estimation of life cycle components of selection in an experimental plant population.

Viability and fertility components of selection associated with linkage blocks marked by four electrophoretically detectable loci were estimated in an experimental population of barley [Composite Cross V (CCV)]. The intensity of selection affecting the distribution of pollen types in the outcross pool was also estimated and comparisons were made between the selective values of genes in the pools of uniting ovules and pollen. The estimates show that selection was intense at various stages of the life cycle and that viability and fertility components often opposed one another. Estimates of viability and fertility components of selection were also extended to the three-locus level. The multilocus estimates reveal large differences in viability and fertility among homozygous genotypes. It is concluded that strong selection operates at all life cycle stages in CCV, although often in differing directions.

Regions of stable equilibria for models of differential selection in the two sexes under random mating.

The equilibrium structure of models of differential selection in the sexes is investigated. It is shown that opposing additive selection leads to stable polymorphic equilibria for only a restricted set of selection intensities, and that for weak selection intensities must be of approximately the same magnitude in the sexes. General models of opposing directional selection, with arbitrary dominance, are investigated by considering simultaneously the stability properties of the trivial equilibria and the curve along which multiple roots appear. Numerical calculations lead us to infer that the average degree of dominance determines the equilibrium characteristics of models of opposing selection. It appears that if the favored alleles are, on the average, recessive, there may be multiple polymorphic equilibria, whereas only a single polymorphic equilibrium can occur when the favored alleles are, on the average, dominant. The principle that the average degree of dominance controls equilibrium behavior is then extended to models allowing directional selection in one sex with overdominance in the other sex, by showing that polymorphism is maintained if and only if the average fitness in heterozygotes exceeds one.

Dynamics of correlated genetic systems. I. Selection in the region of the Glued locus of Drosophila melanogaster.

The dynamical behavior of chromosomal segments undergoing strong selection was investigated in four replicate populations of Drosophila melanogaster. This was accomplished by following the joint behavior of allozyme markers at the loci phosphoglucomutase and esterase C, adjacent to the recessive lethal locus Glued, during and following the course of selection against Glued. The results show strong selection at other loci in the region of the marked segment. Examination of the joint dynamic of the two markers indicates that there must be more than one, and probably several, selected loci in the region under observation, with large epistatic effects. The mode of selection on the segment often results in excess heterozygosity at the markers, but does not appear to be constant in time. It is concluded that the density of selective effects in the region under study is substantial.

A restriction map of the ribosomal RNA genes and the short single-copy DNA sequence of the pearl millet chloroplast genome.

The chloroplast rDNA genes of pearl millet (Pennisetum americanum) have been cloned and physically mapped. The chloroplast genome of the pearl millet contains two identical rRNA genes located on DNA sequences that are inverted with respect to one another and separated by 12 kb of single-copy DNA. The rRNA genes were positioned on a restriction endonuclease map by using as hybridization probes specific cloned rDNA sequences from the chloroplast DNA of the alga Euglena gracilis. The 16S and 23S rRNA genes were shown to be approx. 2 kb from one another, and the 5S RNA gene is immediately adjacent to the 23S tRNA gene.

Comparative genomics of chalcone synthase and Myb genes in the grass family.

Most plant genes occur as members of multigene families where new copies arise through duplication. Duplicate genes that do not confer an adaptive advantage to the plant are expected to rapidly erode into pseudogenes owing to the accumulation of transpositions, insertion/deletion mutations and nucleotide changes. Nonfunctional copies will drift to fixation within a few million years and ultimately erode beyond recognition. Duplicate genes that are retained over longer periods of evolutionary time must be positively selected based on some adaptive advantage conferred on the plant species. We explore the dynamics of the recruitment of new duplicate genes for chalcone synthase, the enzyme that catalyzes the first committed step of flavonoid biosynthesis, and for the myb family of transcriptional activators. Our analyses show that new chs genes are recruited into the genome of grasses at a rate of one new copy every 15 to 25 million years. In contrast, the myb gene family is much older and many duplicate copies appear to predate the separation of the angiosperm lineage from other seed plants. The general pattern suggests a rapid adaptive proliferation of new chs genes but a more ancient elaboration of regulatory gene functions. Our analyses also reveal accelerated rates of protein evolution following gene duplication and evidence is presented for interlocus exchange among duplicate gene loci.