Sequence diversity in the tetraploid Zea perennis and the closely related diploid Z. diploperennis: insights from four nuclear loci

Mysterious meiotic behavior of autopolyploid and allopolyploid maize

This study was aimed to investigate the stability of chromosomes during meiosis in autopolyploid and allopolyploid maize, as well as to determine an association of chromosomes between maize (Zea mays ssp. mays Linnaeus, 1753) and Z. perennis (Hitchcock, 1922) Reeves & Mangelsdor, 1942, by producing a series of autopolyploid and allopolyploid maize hybrids. The intra-genomic and inter-genomic meiotic pairings in these polyploids were quantified and compared using dual-color genomic in-situ hybridization. The results demonstrated higher level of chromosome stability in allopolyploid maize during meiosis as compared to autopolyploid maize. In addition, the meiotic behavior of Z. perennis was relatively more stable as compared to the allopolyploid maize. Moreover, ten chromosomes of "A" subgenome in maize were homologous to twenty chromosomes of Z. perennis genome with a higher pairing frequency and little evolutionary differentiation. At the same time, little evolutionary differentiation has been shown by chromosomes of "A" subgenome in maize, while chromosomes of "B" subgenome, had a lower pairing frequency and higher evolutionary differentiation. Furthermore, 5IM + 5IIPP + 5IIIMPP and 5IIMM + 5IIPP + 5IVMMPP were observed in allotriploids and allotetraploids respectively, whereas homoeologous chromosomes were found between the "A" and "B" genome of maize and Z. perennis.

Cytogeography of Callisia section Cuthbertia (Commelinaceae)

Determining the distribution of cytotypes across the geographic distribution of polyploid complexes can provide valuable information about the evolution of biodiversity. Here, the phytogeography of cytotypes in section Cuthbertia (Small, 1903) Hunt, 1986 is investigated. A total of 436 voucher specimens was georeferenced; 133 new specimens were collected. Based on flow cytometry data, DNA content of all cytotypes in section Cuthbertia was estimated. Utilizing chromosome counts and flow cytometric analysis, cytotype distribution maps were generated. Two disjunct groups of populations of diploid Callisia graminea (Small, 1903) Tucker, 1989 were discovered; tetraploid C. graminea ranges broadly from the coastal plain of North Carolina through central Florida. One hexaploid C. graminea individual was recorded in South Carolina, and numerous individuals of hexaploid C. graminea were found in central Florida. Diploid C. ornata (Small, 1933) Tucker, 1989 occurs in eastern Florida; previously unknown tetraploid and hexaploid populations of C. ornata were discovered in western and central Florida, respectively. Diploid C. rosea (Ventenat, 1800) Hunt, 1986 occurs in Georgia and the Carolinas, with populations occurring on both sides of the Fall Line. The cytotype and species distributions in Callisia are complex, and these results provide hypotheses, to be tested with morphological and molecular data, about the origins of the polyploid cytotypes.

Patterns of Selection in Plant Genomes

Plants show a wide range of variation in mating system, ploidy level, and demographic history, allowing for unique opportunities to investigate the evolutionary and genetic factors affecting genome-wide patterns of positive and negative selection. In this review, we highlight recent progress in our understanding of the extent and nature of selection on plant genomes. We discuss differences in selection as they relate to variation in demography, recombination, mating system, and ploidy. We focus on the population genetic consequences of these factors and argue that, although variation in the magnitude of purifying selection is well documented, quantifying rates of positive selection and disentangling the relative importance of recombination, demography, and ploidy are ongoing challenges. Large-scale comparative studies that examine the relative and joint importance of these processes, combined with explicit models of population history and selection, are key and feasible goals for future work.

Genetic diversity and molecular evolution of Chinese waxy maize germplasm

Waxy maize (Zea mays L. var. certaina Kulesh), with many excellent characters in terms of starch composition and economic value, has grown in China for a long history and its production has increased dramatically in recent decades. However, the evolution and origin of waxy maize still remains unclear. We studied the genetic diversity of Chinese waxy maize including typical landraces and inbred lines by SSR analysis and the results showed a wide genetic diversity in the Chinese waxy maize germplasm. We analyzed the origin and evolution of waxy maize by sequencing 108 samples, and downloading 52 sequences from GenBank for the waxy locus in a number of accessions from genus Zea. A sharp reduction of nucleotide diversity and significant neutrality tests (Tajima's D and Fu and Li's F*) were observed at the waxy locus in Chinese waxy maize but not in nonglutinous maize. Phylogenetic analysis indicated that Chinese waxy maize originated from the cultivated flint maize and most of the modern waxy maize inbred lines showed a distinct independent origin and evolution process compared with the germplasm from Southwest China. The results indicated that an agronomic trait can be quickly improved to meet production demand by selection.

Extending coalescent theory to autotetraploids

We develop coalescent models for autotetraploid species with tetrasomic inheritance. We show that the ancestral genetic process in a large population without recombination may be approximated using Kingman's standard coalescent, with a coalescent effective population size 4N. Numerical results suggest that this approximation is accurate for population sizes on the order of hundreds of individuals. Therefore, existing coalescent simulation programs can be adapted to study population history in autotetraploids simply by interpreting the timescale in units of 4N generations. We also consider the possibility of double reduction, a phenomenon unique to polysomic inheritance, and show that its effects on gene genealogies are similar to partial self-fertilization.

Three new teosintes (Zea spp., Poaceae) from México

PREMISE OF THE STUDY

Teosinte species are the closest relatives of maize and represent an important but increasingly rare genetic resource for maize improvement and the study of evolution by domestication. Three morphologically and ecologically distinct teosinte populations were recently discovered in México. The taxonomic status of these rare and endangered populations was investigated by detailed comparisons to previously characterized wild Zea species. •

METHODS

Three new teosinte populations were compared to known teosinte taxa on the basis of morphological, ecogeographic, cytological, and molecular characteristics. Phenetic and phylogenetic analyses were performed using morphological and molecular data, respectively. •

KEY RESULTS

The newly discovered populations are distinct from each other and from other Zea species to represent three new entities based on their unique combinations of morphological, ecological, ploidy, and DNA markers. A perennial diploid population from Nayarit is distinguished by early maturing plants, and having male inflorescences with few tassel branches and long spikelets. A perennial tetraploid population from Michoacán is characterized by tall and late maturing plants, and having male inflorescences with many branches. An annual diploid population from Oaxaca is characterized by having male inflorescences with fewer branches and longer spikelets than those found in the sister taxa Z. luxurians and Z. nicaraguensis, plants with high thermal requirements, and very long seed dormancy. •

CONCLUSIONS

Evidence from multiple independent sources suggests placement of the three new populations of teosinte as distinct entities within section Luxuriantes of the genus Zea. However, more extensive DNA marker or sequence data are required to resolve the taxonomy of this genus.

Multilocus analysis of genetic divergence between outcrossing Arabidopsis species: evidence of genome-wide admixture

• Outcrossing Arabidopsis species that diverged from their inbreeding relative Arabidopsis thaliana 5 million yr ago and display a biogeographical pattern of interspecific sympatry vs intraspecific allopatry provides an ideal model for studying impacts of gene introgression and polyploidization on species diversification. • Flow cytometry analyses detected ploidy polymorphisms of 2× and 4× in Arabidopsis lyrata ssp. kamchatica of Taiwan. Genomic divergence between species/subspecies was estimated based on 98 randomly chosen nuclear genes. Multilocus analyses revealed a mosaic genome in diploid A. l. kamchatica composed of Arabidopsis halleri-like and A. lyrata-like alleles. • Coalescent analyses suggest that the segregation of ancestral polymorphisms alone cannot explain the high inconsistency between gene trees across loci, and that gene introgression via diploid A. l. kamchatica likely distorts the molecular phylogenies of Arabidopsis species. However, not all genes migrated across species freely. Gene ontology analyses suggested that some nonmigrating genes were constrained by natural selection. • High levels of estimated ancestral polymorphisms between A. halleri and A. lyrata suggest that gene flow between these species has not completely ceased since their initial isolation. Polymorphism data of extant populations also imply recent gene flow between the species. Our study reveals that interspecific gene flow affects the genome evolution in Arabidopsis.

Nucleotide diversity and molecular evolution of the PSY1 gene in Zea mays compared to some other grass species

Phytoene synthase (PSY), which is encoded by the phytoene synthase 1 (PSY1) gene, is the first rate-limiting enzyme in the plant carotenoid biosynthetic pathway. In order to examine the genetic diversity and evolution pattern of PSY1 within the Andropogoneae, sequences of 76 accessions from 5 species (maize, teosinte, tripsacum, coix, and sorghum) of the Andropogoneae were tested, along with 4 accessions of rice (Oryza sativa L.) included as outliers. Both the number and the order of exons and introns were relatively conserved across the species tested. Three domains were identified in the coding sequence, including signal peptide (SP), PSY, and highly conserved squalene synthase (SQS) domain. Although no positive selection signal was detected at an overall coding level among all species tested, the SP domain and the region upstream of the SQS-PSY domain appear to have undergone rapid evolution, as evidenced by a high d (N)/d (S) ratio (>1.0). At the nucleotide level, positive selection and balancing selection were detected only among the yellow maize germplasm and the white maize germplasm, respectively. The phylogenetic tree based on full-length sequences of PSY1-like regions supported the monophyletic theory of the Andropogoneae and the closest relationship between Zea and Tripsacum among the Andropogoneae. Coix, which was theorized to have a closer relationship with maize due to similarities in morphology and chromosome number, has been shown in this study to have diverged relatively early from the other Andropogoneae, including maize.

Post-domestication selection in the maize starch pathway

Modern crops have usually experienced domestication selection and subsequent genetic improvement (post-domestication selection). Chinese waxy maize, which originated from non-glutinous domesticated maize (Zea mays ssp. mays), provides a unique model for investigating the post-domestication selection of maize. In this study, the genetic diversity of six key genes in the starch pathway was investigated in a glutinous population that included 55 Chinese waxy accessions, and a selective bottleneck that resulted in apparent reductions in diversity in Chinese waxy maize was observed. Significant positive selection in waxy (wx) but not amylose extender1 (ae1) was detected in the glutinous population, in complete contrast to the findings in non-glutinous maize, which indicated a shift in the selection target from ae1 to wx during the improvement of Chinese waxy maize. Our results suggest that an agronomic trait can be quickly improved into a target trait with changes in the selection target among genes in a crop pathway.

Sox genes evolution in closely related young tetraploid cyprinid fishes and their diploid relative

Previous study and analysis of cytochrome b suggested that polyploidization event in the genus Tor occurred about 10 Mya ago. In order to understand evolutionary fates of Sox gene in the early stage of genome duplication at the nucleotide level, PCR surveys for Sox genes in three closely related cyprinid fishes T. douronensis (2n=100), T. qiaojiensis (2n=?), T. sinensis (2n=100) and their relative T. brevifilis (2n=50) were performed. Totally, 52 distinct Sox genes were obtained in these four species, representing SoxB, SoxC, and SoxE group. As expected, isoforms of some Sox genes correspond with the ploidy of species, such as two copies of Sox9a exist in tetraploid species. Analysis indicated that duplicated Sox gene pairs caused by polyploidization evolved independently of each other within polyploid species. Results of substitution rate showed nearly equal rate of nonsynonymous substitution of duplicated Sox orthologs among different polyploid species and their diploid relative orthologs, suggesting at the early stage of genome duplicated Sox orthologs are under similar selective constraints in different polyploidy species and their diploid relative at the amino acid level. All PCR fragments of Sox genes obtained in this study are not accompanied by obvious increase in mutations and pseudogene formation which means that they are under strong purifying selection, suggesting that they are functional at the DNA level. Genealogical analysis revealed that T. qiaojiensis was tetraploid, and T. douronensis, T. qiaojiensis as well as T. sinensis had an allotetraploid ancestor.

Historical divergence and gene flow in the genus Zea

Gene flow plays a fundamental role in plant evolutionary history, yet its role in population divergence--and ultimately speciation--remains poorly understood. We investigated gene flow and the modalities of divergence in the domesticated Zea mays ssp. mays and three wild Zea taxa using sequence polymorphism data from 26 nuclear loci. We described diversity across loci and assessed evidence for adaptive and purifying selection at nonsynonymous sites. For each of three divergence events in the history of these taxa, we used approximate Bayesian simulation to estimate population sizes and divergence times and explicitly compare among alternative models of divergence. Our estimates of divergence times are surprisingly consistent with previous data from other markers and suggest rapid diversification of lineages within Zea in the last approximately 150,000 years. We found widespread evidence of historical gene flow, including evidence for divergence in the face of gene flow. We speculate that cultivated maize may serve as a bridge for gene flow among otherwise allopatric wild taxa.

Using multilocus sequence data to assess population structure, natural selection, and linkage disequilibrium in wild tomatoes

We employed a multilocus approach to examine the effects of population subdivision and natural selection on DNA polymorphism in 2 closely related wild tomato species (Solanum peruvianum and Solanum chilense), using sequence data for 8 nuclear loci from populations across much of the species' range. Both species exhibit substantial levels of nucleotide variation. The species-wide level of silent nucleotide diversity is 18% higher in S. peruvianum (pi(sil) approximately 2.50%) than in S. chilense (pi(sil) approximately 2.12%). One of the loci deviates from neutral expectations, showing a clinal pattern of nucleotide diversity and haplotype structure in S. chilense. This geographic pattern of variation is suggestive of an incomplete (ongoing) selective sweep, but neutral explanations cannot be entirely dismissed. Both wild tomato species exhibit moderate levels of population differentiation (average F(ST) approximately 0.20). Interestingly, the pooled samples (across different demes) exhibit more negative Tajima's D and Fu and Li's D values; this marked excess of low-frequency polymorphism can only be explained by population (or range) expansion and is unlikely to be due to population structure per se. We thus propose that population structure and population/range expansion are among the most important evolutionary forces shaping patterns of nucleotide diversity within and among demes in these wild tomatoes. Patterns of population differentiation may also be impacted by soil seed banks and historical associations mediated by climatic cycles. Intragenic linkage disequilibrium (LD) decays very rapidly with physical distance, suggesting high recombination rates and effective population sizes in both species. The rapid decline of LD seems very promising for future association studies with the purpose of mapping functional variation in wild tomatoes.

Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice

Varying degrees of reduction of genetic diversity in crops relative to their wild progenitors occurred during the process of domestication. Such information, however, has not been available for the Asian cultivated rice (Oryza sativa) despite its importance as a staple food and a model organism. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationship and demographic history of O. sativa and its close relatives (Oryza rufipogon and Oryza nivara), we investigated nucleotide diversity data from 10 unlinked nuclear loci in species-wide samples of these species. The results indicated that O. rufipogon and O. nivara possessed comparable levels of nucleotide variation ((sil) = 0.0077 approximately 0.0095) compared with the relatives of other crops. In contrast, nucleotide diversity of O. sativa was as low as (sil) = 0.0024 and even lower ((sil) = 0.0021 for indica and 0.0011 for japonica), if we consider the 2 subspecies separately. Overall, only 20-10% of the diversity in the wild species was retained in 2 subspecies of the cultivated rice (indica and japonica), respectively. Because statistic tests did not reject the assumption of neutrality for all 10 loci, we further used coalescent to simulate bottlenecks under various lengths and population sizes to better understand the domestication process. Consistent with the dramatic reduction in nucleotide diversity, we detected a severe domestication bottleneck and demonstrated that the sequence diversity currently found in the rice genome could be explained by a founding population of 1,500 individuals if the initial domestication event occurred over a 3,000-year period. Phylogenetic analyses revealed close genetic relationships and ambiguous species boundary of O. rufipogon and O. nivara, providing additional evidence to treat them as 2 ecotypes of a single species. Lowest linkage disequilibrium (LD) was found in the perennial O. rufipogon where the r(2) value dropped to a negligible level within 400 bp, and the highest in the japonica rice where LD extended to the entirely sequenced region ( approximately 900 bp), implying that LD mapping by genome scans may not be feasible in wild rice due to the high density of markers needed.

Multilocus patterns of nucleotide diversity, linkage disequilibrium and demographic history of Norway spruce [Picea abies (L.) Karst]

DNA polymorphism at 22 loci was studied in an average of 47 Norway spruce [Picea abies (L.) Karst.] haplotypes sampled in seven populations representative of the natural range. The overall nucleotide variation was limited, being lower than that observed in most plant species so far studied. Linkage disequilibrium was also restricted and did not extend beyond a few hundred base pairs. All populations, with the exception of the Romanian population, could be divided into two main domains, a Baltico-Nordic and an Alpine one. Mean Tajima's D and Fay and Wu's H across loci were both negative, indicating the presence of an excess of both rare and high-frequency-derived variants compared to the expected frequency spectrum in a standard neutral model. Multilocus neutrality tests based on D and H led to the rejection of the standard neutral model and exponential growth in the whole population as well as in the two main domains. On the other hand, in all three cases the data are compatible with a severe bottleneck occurring some hundreds of thousands of years ago. Hence, demographic departures from equilibrium expectations and population structure will have to be accounted for when detecting selection at candidate genes and in association mapping studies, respectively.

Domestication history in the Medicago sativa species complex: inferences from nuclear sequence polymorphism

DNA sequence polymorphism carries genealogical information and allows for testing hypotheses on selection and population history, especially through coalescent-based analysis. Understanding the evolutionary forces at work in plant domestication and subsequent selection is of critical importance for the management of genetic resources. In this study, we surveyed DNA sequence diversity at two assumed neutral nuclear loci in the wild-domesticated species complex of alfalfa (Medicago sativa L.). A high level of polymorphism was detected. The domesticated pool contains on average 31% less diversity than the wild pool, but with a high heterogeneity among loci. Coalescent simulations of the domestication process showed that this result cannot be explained by assuming a constant population size but is rather consistent with a demographic bottleneck during domestication. A very low level of divergence was detected between the wild and the domesticated forms as well as between the related subspecies of the M. sativa species complex. However, the originality of the Spanish wild populations, already observed based on mitochondrial DNA polymorphism, was confirmed. These results, together with patterns of intrapopulation polymorphism, suggest that nuclear sequence polymorphism could be a promising tool, complementary to mitochondrial DNA and phenotypic evaluations, to investigate historical demographic and evolutionary processes.

Patterns of nucleotide diversity in wild and cultivated sunflower

Interest in the level and organization of nucleotide diversity in domesticated plant lineages has recently been motivated by the potential for using association-based mapping techniques as a means for identifying the genes underlying complex traits. To date, however, such data have been available only for a relatively small number of well-characterized plant taxa. Here we provide the first detailed description of patterns of nucleotide polymorphism in wild and cultivated sunflower (Helianthus annuus), using sequence data from nine nuclear genes. The results of this study indicate that wild sunflower harbors at least as much nucleotide diversity as has been reported in other wild plant taxa, with randomly selected sequence pairs being expected to differ at 1 of every 70 bp. In contrast, cultivated sunflower has retained only 40-50% of the diversity present in the wild. Consistent with this dramatic reduction in polymorphism, a phylogenetic analysis of our data revealed that the cultivars form a monophyletic clade, adding to the growing body of evidence that sunflower is the product of a single domestication. Eight of the nine loci surveyed appeared to be evolving primarily under purifying selection, while the remaining locus may have been the subject of positive selection. Linkage disequilibrium (LD) decayed very rapidly in the self-incompatible wild sunflower, with the expected LD falling to negligible levels within 200 bp. The cultivars, on the other hand, exhibited somewhat higher levels of LD, with nonrandom associations persisting up to approximately 1100 bp. Taken together, these results suggest that association-based approaches will provide a high degree of resolution for the mapping of functional variation in sunflower.

Genetic diversity and the evolutionary history of plant immunity genes in two species of Zea

Plant pathogenesis-related genes (PR genes) code for enzymes, enzyme inhibitors, and other peptides that confer resistance to pathogens and herbivores. Although several PR genes have been the subject of molecular population genetic analyses, a general understanding of their long-term evolutionary dynamics remains incomplete. Here we analyze sequence data from 17 PR genes from two closely related teosinte species of central Mexico. In addition to testing whether patterns of diversity at individual loci depart from expectations under a neutral model, we compared patterns of diversity at defense genes, as a class, to nondefense genes. In Zea diploperennis, the majority of defense genes have patterns of diversity consistent with neutral expectations while at least two genes showed evidence of recent positive selection consistent with arms-race models of antagonistic coevolution. In Zea mays ssp. parviglumis, by contrast, analyses of both defense and nondefense genes revealed strong and consistent departures from the neutral model, suggestive of nonequilibrium population dynamics or population structure. Nevertheless, we found a significant excess of replacement polymorphism in defense genes compared to nondefense genes. Although we cannot exclude relaxed selective constraint as an explanation, our results are consistent with temporally variable (transient or episodic) selection or geographically variable selection acting on parviglumis defense genes. The different patterns of diversity found in the two Zea species may be explained by parviglumis' greater distribution and population structure together with geographic variation in selection.

Sequence polymorphism in polyploid wheat and their d-genome diploid ancestor

Sequencing was used to investigate the origin of the D genome of the allopolyploid species Triticum aestivum and Aegilops cylindrica. A 247-bp region of the wheat D-genome Xwye838 locus, encoding ADP-glucopyrophosphorylase, and a 326-bp region of the wheat D-genome Gss locus, encoding granule-bound starch synthase, were sequenced in a total 564 lines of hexaploid wheat (T. aestivum, genome AABBDD) involving all its subspecies and 203 lines of Aegilops tauschii, the diploid source of the wheat D genome. In Ae. tauschii, two SNP variants were detected at the Xwye838 locus and 11 haplotypes at the Gss locus. Two haplotypes with contrasting frequencies were found at each locus in wheat. Both wheat Xwye838 variants, but only one of the Gss haplotypes seen in wheat, were found among the Ae. tauschii lines. The other wheat Gss haplotype was not found in either Ae. tauschii or 70 lines of tetraploid Ae. cylindrica (genomes CCDD), which is known to hybridize with wheat. It is concluded that both T. aestivum and Ae. cylindrica originated recurrently, with at least two genetically distinct progenitors contributing to the formation of the D genome in both species.

Population genetic evidence for rapid changes in intraspecific diversity and allelic cycling of a specialist defense gene in Zea

Two patterns of plant defense gene evolution are emerging from molecular population genetic surveys. One is that specialist defenses experience stronger selection than generalist defenses. The second is that specialist defenses are more likely to be subject to balancing selection, i.e., evolve in a manner consistent with balanced-polymorphism or trench-warfare models of host-parasite coevolution. Because most of the data of specialist defenses come from Arabidopsis thaliana, we examined the genetic diversity and evolutionary history of three defense genes in two outcrossing species, the autotetraploid Zea perennis and its most closely related extant relative the diploid Z. diploperennis. Intraspecific diversity at two generalist defenses, the protease inhibitors wip1 and mpi, were consistent with a neutral model. Like previously studied genes in these taxa, wip1 and mpi harbored similar levels of diversity in Z. diploperennis and Z. perennis. In contrast, the specialist defense hm2 showed strong although distinctly different departures from a neutral model in the two species. Z. diploperennis appears to have experienced a strong and recent selective sweep. Using a rejection-sampling coalescent method, we estimate the strength of selection on Z. diploperennis hm2 to be approximately 3.0%, which is approximately equal to the strength of selection on tb1 during maize domestication. Z. perennis hm2 harbors three highly diverged alleles, two of which are found at high frequency. The distinctly different patterns of diversity may be due to differences in the phase of host-parasite coevolutionary cycles, although higher hm2 diversity in Z. perennis may also reflect reduced efficacy of selection in the autotetraploid relative to its diploid relative.

Comparative evolutionary histories of chitinase genes in the Genus zea and Family poaceae

Patterns of DNA sequence diversity vary widely among genes encoding proteins that protect plants against pathogens and herbivores. Comparative studies may help determine whether these differences are due to the strength of selection acting on different types of defense, in different evolutionary lineages, or both. I analyzed sequence diversity at three chitinases, a well-studied component of defense, in two species of Zea and several Poaceae taxa. Although the Zea species are closely related and these genes code for proteins with similar biochemical function, patterns of diversity vary widely within and among species. Intraspecific diversity at chiB, chiI, and Z. mays ssp. parviglumis chiA are consistent with a neutral-equilibrium model whereas chiA had no segregating sites within Z. diploperennis--consistent with a recent and strong selective sweep. Codons identified as having diverged among Poaceae taxa in response to positive selection were significantly overrepresented among targets of selection in Arabis, suggesting common responses to selection in distantly related plant taxa. Divergence of the recent duplicates chiA and chiB is consistent with positive selection but relaxed constraint cannot be rejected. Weak evidence for adaptive divergence of these duplicated downstream components of defense contrasts with strong evidence for adaptive divergence of genes involved in pathogen recognition.

Multilocus analysis of variation and speciation in the closely related species Arabidopsis halleri and A. lyrata

Nucleotide variation in eight effectively unlinked genes was surveyed in species-wide samples of the closely related outbreeding species Arabidopsis halleri and A. lyrata ssp. petraea and in three of these genes in A. lyrata ssp. lyrata and A. thaliana. Significant genetic differentiation was observed more frequently in A. l. petraea than in A. halleri. Average estimates of nucleotide variation were highest in A. l. petraea and lowest in A. l. lyrata, reflecting differences among species in effective population size. The low level of variation in A. l. lyrata is concordant with a bottleneck effect associated with its origin. The A. halleri/A. l. petraea speciation process was studied, considering the orthologous sequences of an outgroup species (A. thaliana). The high number of ancestral mutations relative to exclusive polymorphisms detected in A. halleri and A. l. petraea, the significant results of the multilocus Fay and Wu H tests, and haplotype sharing between the species indicate introgression subsequent to speciation. Average among-population variation in A. halleri and A. l. petraea was approximately 1.5- and 3-fold higher than that in the inbreeder A. thaliana. The detected reduction of variation in A. thaliana is less than that expected from differences in mating system alone, and therefore from selective processes related to differences in the effective recombination rate, but could be explained by differences in population structure.

Nucleotide polymorphism and linkage disequilibrium within and among natural populations of European aspen (Populus tremula L., Salicaceae)

Populus is an important model organism in forest biology, but levels of nucleotide polymorphisms and linkage disequilibrium have never been investigated in natural populations. Here I present a study on levels of nucleotide polymorphism, haplotype structure, and population subdivision in five nuclear genes in the European aspen Populus tremula. Results show substantial levels of genetic variation. Levels of silent site polymorphisms, pi(s), averaged 0.016 across the five genes. Linkage disequilibrium was generally low, extending only a few hundred base pairs, suggesting that rates of recombination are high in this obligate outcrossing species. Significant genetic differentiation was found at all five genes, with an average estimate of F(ST) = 0.116. Levels of polymorphism in P. tremula are 2- to 10-fold higher than those in other woody, long-lived perennial plants, such as Pinus and Cryptomeria. The high levels of nucleotide polymorphism and low linkage disequilibrium suggest that it may be possible to map functional variation to very fine scales in P. tremula using association-mapping approaches.

Pattern of diversity in the genomic region near the maize domestication gene tb1

Domesticated maize and its wild ancestor (teosinte) differ strikingly in morphology and afford an opportunity to examine the connection between strong selection and diversity in a major crop species. The tb1 gene largely controls the increase in apical dominance in maize relative to teosinte, and a region of the tb1 locus 5' to the transcript sequence was a target of selection during maize domestication. To better characterize the impact of selection at a major "domestication" locus, we have sequenced the upstream tb1 genomic region and systematically sampled nucleotide diversity for sites located as far as 163 kb upstream to tb1. Our analyses define a selective sweep of approximately 60-90 kb 5' to the tb1 transcribed sequence. The selected region harbors a mixture of unique sequences and large repetitive elements, but it contains no predicted genes. Diversity at the nearest 5' gene to tb1 is typical of that for neutral maize loci, indicating that selection at tb1 has had a minimal impact on the surrounding chromosomal region. Our data also show low intergenic linkage disequilibrium in the region and suggest that selection has had a minor role in shaping the pattern of linkage disequilibrium that is observed. Finally, our data raise the possibility that maize-like tb1 haplotypes are present in extant teosinte populations, and our findings also suggest a model of tb1 gene regulation that differs from traditional views of how plant gene expression is controlled.

Molecular evolution of the wound-induced serine protease inhibitor wip1 in Zea and related genera

Plant defense mechanisms have been the subject of intensive investigation. However, little is known about their long-term evolutionary dynamics. We investigated the molecular diversity of a wound-induced serine protease inhibitor, wip1, in the genus Zea, as well as the divergence of wip1 among four genera, Zea, Tripsacum, Sorghum, and Oryza, in order to gain insight into the long-term evolution of plant defense. The specific objectives of this study were to determine (1) whether wip1 has a history of positive or balancing selection, as has been shown for genes involved in plant defense against pathogens, and (2) if the evolutionary histories of wip1 inhibitory loops, which come into closest contact with proteases, differ from the evolutionary history of other parts of this gene. The Zea polymorphism data are consistent with a neutral evolutionary history. In contrast, relative-rate tests suggest a nonneutral evolutionary history. This inconsistency may indicate that selection acting on wip1 is episodic or that wip1 evolves in response to selection favoring novel alleles. We also detected significant heterogeneity in the evolutionary rates of the two inhibitory loops of wip1-one inhibitory loop is highly conserved, whereas the second has diverged rapidly. Because these two inhibitory loops are predicted to have very similar biochemical functions, the significantly different evolutionary histories suggest that these loops have different ecological functions.