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 domesticate 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 ~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.     

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.     

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.     

Effective population size is positively correlated with levels of adaptive divergence among annual sunflowers

The role of adaptation in the divergence of lineages has long been a central question in evolutionary biology, and as multilocus sequence data sets have become available for a wide range of taxa, empirical estimates of levels of adaptive molecular evolution are increasingly common. Estimates vary widely among taxa, with high levels of adaptive evolution in Drosophila, bacteria, and viruses but very little evidence of widespread adaptive evolution in hominids. Although estimates in plants are more limited, some recent work has suggested that rates of adaptive evolution in a range of plant taxa are surprisingly low and that there is little association between adaptive evolution and effective population size in contrast to patterns seen in other taxa. Here, we analyze data from 35 loci for six sunflower species that vary dramatically in effective population size. We find that rates of adaptive evolution are positively correlated with effective population size in these species, with a significant fraction of amino acid substitutions driven by positive selection in the species with the largest effective population sizes but little or no evidence of adaptive evolution in species with smaller effective population sizes. Although other factors likely contribute as well, in sunflowers effective population size appears to be an important determinant of rates of adaptive evolution.     

Speciation and domestication in maize and its wild relatives: evidence from the globulin-1 gene.

The grass genus Zea contains the domesticate maize and several wild taxa indigenous to Central and South America. Here we study the genetic consequences of speciation and domestication in this group by sampling DNA sequences from four taxa-maize (Zea mays ssp. mays), its wild progenitor (Z. mays ssp. parviglumis), a more distant species within the genus (Z. luxurians), and a representative of the sister genus (Tripsacum dactyloides). We sampled a total of 26 sequences from the glb1 locus, which encodes a nonessential seed storage protein. Within the Zea taxa sampled, the progenitor to maize contains the most sequence diversity. Maize contains 60% of the level of genetic diversity of its progenitor, and Z. luxurians contains even less diversity (32% of the level of diversity of Z. mays ssp. parviglumis). Sequence variation within the glb1 locus is consistent with neutral evolution in all four taxa. The glb1 data were combined with adh1 data from a previous study to make inferences about the population genetic histories of these taxa. Comparisons of sequence data between the two morphologically similar wild Zea taxa indicate that the species diverged approximately 700, 000 years ago from a common ancestor of intermediate size to their present populations. Conversely, the domestication of maize was a recent event that could have been based on a very small number of founding individuals. Maize retained a substantial proportion of the genetic variation of its progenitor through this founder event, but diverged rapidly in morphology.     

Allelic imbalance in Drosophila hybrid heads: exons, isoforms, and evolution

Unraveling how regulatory divergence contributes to species differences and adaptation requires identifying functional variants from among millions of genetic differences. Analysis of allelic imbalance (AI) reveals functional genetic differences in cis regulation and has demonstrated differences in cis regulation within and between species. Regulatory mechanisms are often highly conserved, yet differences between species in gene expression are extensive. What evolutionary forces explain widespread divergence in cis regulation? AI was assessed in Drosophila melanogaster-Drosophila simulans hybrid female heads using RNA-seq technology. Mapping bias was virtually eliminated by using genotype-specific references. Allele representation in DNA sequencing was used as a prior in a novel Bayesian model for the estimation of AI in RNA. Cis regulatory divergence was common in the organs and tissues of the head with 41% of genes analyzed showing significant AI. Using existing population genomic data, the relationship between AI and patterns of sequence evolution was examined. Evidence of positive selection was found in 30% of cis regulatory divergent genes. Genes involved in defense, RNAi/RISC complex genes, and those that are sex regulated are enriched among adaptively evolving cis regulatory divergent genes. For genes in these groups, adaptive evolution may play a role in regulatory divergence between species. However, there is no evidence that adaptive evolution drives most of the cis regulatory divergence that is observed. The majority of genes showed patterns consistent with stabilizing selection and neutral evolutionary processes.     

Zea systematics: ribosomal ITS evidence

Ribosomal internal transcribed spacer (ITS) sequences were used to evaluate the phylogenetics of Zea and Tripsacum. Maximum likelihood and polymorphism parsimony were used for phylogenetic reconstructions. Zea ITS nucleotide diversity was high compared to other plant species, but approximately equivalent to other maize loci. Coalescence of ITS alleles was rapid relative to other nuclear loci; however, there was still much diversity within populations. Zea and Tripsacum form a clade clearly differentiated from all other Poaceae. Four Zea ITS pseudogenes were identified by phylogenetic position and nucleotide composition. The phylogenetic position of Z. mays ssp. huehuetenangensis was clearly established as basal to the other Z. mays. The ITS phylogeny disfavored a Z. luxurians and Z. diploperennis clade, which conflicted with some previous studies. The introgression of Z. mays alleles into Z. perennis and Z. diploperennis was also established. The ITS data indicated a near contemporary divergence of domesticated maize and its two closest wild relatives.      

Cloning of the genes of the chitin utilization regulon of Serratia liquefaciens.

The set of genes that determine the expression of the enzymes involved in chitin degradation by Serratia liquefaciens was cloned. The role of each gene was investigated, and for the first time regulatory genes were identified in this system. The chiA and chiB genes coded for separate chitinase activities. The chiC region coded for a chitobiase activity, but it was not formally separated from chiB. Transposon mutagenesis and deletion analysis identified a region, chiD, whose absence led to higher expression of chiA, chiB, and chiC. chiD may therefore be a gene that codes for a repressor. Loss of function of another adjacent region, chiE, prevented induction unless a chiE+ strain was a near neighbor, suggesting that this gene may code for a protein that is involved in the synthesis of the inducer. chiB, chiC, chiD, and chiE are closely linked, while chiA is in a separate location on the chromosome.     

Contrasting modes of natural selection acting on pigmentation genes in the Drosophila dunni subgroup

Genes that encode for divergent adaptive traits may have genealogies that contrast with those from loci that are not functionally involved in differentiation. Here, we examine DNA sequence variation among the species of the eastern Caribbean Drosophila dunni subgroup at two loci, yellow and dopa decaboxylase (Ddc), which both play integral roles in pigmentation patterning of adult Drosophila. Phylogenetic analyses of these loci produce gene genealogies with topologies that mirror those described for other nuclear genes: the six morphologically distinct species within the subgroup are divided into only three lineages, with one lineage containing four species that share extensive ancestral polymorphism. At the Ddc locus these major lineages are delineated only by silent site variation. We observe a significantly higher rate of synonymous site divergence than non-synonymous divergence, consistent with strong purifying selection acting on the locus. In contrast, the yellow locus exhibits patterns of amino acid divergence and nucleotide diversity that are consistent with recent diversifying selection acting in two different lineages. This selection appears to be targeting amino acid variants in the signal sequence of the Yellow protein, a region which is tightly constrained among members of the larger D. cardini radiation. This result highlights not only the potential importance of yellow in the evolution of divergent pigmentation patterns among members of the D. dunni subgroup, but also hints that variation in signal peptide sequences may play a role in phenotypic diversification.     

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.     

Pollen- and anther-specific chi promoters from petunia: tandem promoter regulation of the chiA gene.

We have analyzed the spatial and temporal activities of chalcone flavanone isomerase (chi) A and B gene promoters from petunia. To study the tandem promoter regulation of chiA, various chiA promoter fragments were fused with the beta-glucuronidase (GUS) reporter gene. Analysis of transgenic plants containing these chimeric genes provided definitive proof that the chiA coding region is regulated by two distinct promoters (designated PA1 and PA2). We also showed that both promoters can function independently and that the chiA PA1 promoter is expressed in limb (epidermal and parenchyma cells), tube (inner epidermal and parenchyma cells), seed (seed coat, endosperm, and embryo), sepal, leaf, and stem. The use of chiA and chiB promoters in the regulation of anther- and pollen-specific gene expression has been studied. By analyzing transgenic plants containing chimeric genes consisting of chiA and B promoter fragments and the GUS reporter gene, we were able to identify a 0.44-kilobase chiA PA2 promoter fragment that drives pollen-specific gene expression and a 1.75-kilobase chiB PB promoter fragment that confers anther-specific (pollen and tapetum cells) expression to the GUS gene.     

Widespread adaptive evolution of Drosophila genes with sex-biased expression

Many genes in higher eukaryotes show sexually dimorphic expression, and these genes tend to be among the most divergent between species. In most cases, however, it is not known whether this rapid divergence is caused by positive selection or if it is due to a relaxation of selective constraint. To distinguish between these two possibilities, we surveyed DNA sequence polymorphism in 91 Drosophila melanogaster genes with male-, female-, or nonsex-biased expression and determined their divergence from the sister species D. simulans. Using several single- and multilocus statistical tests, we estimated the type and strength of selection influencing the evolution of the proteins encoded by genes of each expression class. Adaptive evolution, as indicated by a relative excess of nonsynonymous divergence between species, was common among the sex-biased genes (both male and female). Male-biased genes, in particular, showed a strong and consistent signal of positive selection, while female-biased genes showed more variation in the type of selection they experience. Genes expressed equally in the two sexes, in contrast, showed no evidence for adaptive evolution between D. melanogaster and D. simulans. This suggests that sexual selection and intersexual coevolution are the major forces driving genetic differentiation between species.     

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

Polyploidy has been an extremely common phenomenon in the evolutionary history of angiosperms. Despite this there are few data available to evaluate the effects of polyploidy on genetic diversity and to compare the relative effects of drift and selection in polyploids and related diploids. We investigated DNA sequence diversity at four nuclear loci (adh1, glb1, c1, and waxy) from the tetraploid Zea perennis and the closely related diploid Z. diploperennis. Contrary to expectations, we detected no strong evidence for greater genetic diversity in the tetraploid, or for consistent differences in the effects of either drift or selection between the tetraploid and the diploid. Our failure to find greater genetic diversity in Z. perennis may result from its relatively recent origin or demographic factors associated with its origin. In addition to comparing genetic diversity in the two species, we constructed genealogies to infer the evolutionary origin of Z. perennis. Although these genealogies are equivocal regarding the mode of origin, several aspects of these genealogies support an autotetraploid origin. Consistent with previous molecular data the genealogies do not, however, support the division of Zea into two sections, the section Zea and the section Luxuriantes.     

The molecular evolution of terminal ear1, a regulatory gene in the genus Zea.

Nucleotide diversity in the terminal ear1 (te1) gene, a regulatory locus hypothesized to be involved in the morphological evolution of maize (Zea mays ssp. mays), was investigated for evidence of past selection. Nucleotide polymorphism in a 1.4-kb region of te1 was analyzed for a sample of 26 sequences isolated from 12 maize lines, five populations of the maize progenitor, Z. mays ssp. parviglumis, six other Zea populations, and two Tripsacum species. Although nucleotide diversity in te1 in maize is reduced relative to ssp. parviglumis, phylogenetic and statistical analyses of the pattern of polymorphism among these sequences provided no evidence of past selection, indicating that the region of the gene studied was probably not involved in maize evolution. The level of reduction in genetic diversity in te1 in maize relative to its progenitor is comparable to that found in previous reports for isozymes and other neutrally evolving maize genes and is consistent with a genome-wide reduction of genetic diversity resulting from a domestication bottleneck. An estimate of the age (1.2-1.4 million yr) of the maize gene pool based on te1 is roughly consistent with previous estimates based on other neutral genes, but may be biased by the apparently slow synonymous substitution rate at te1.     

苏云金芽孢杆菌chiA,chiB全基因的克隆、表达及其序列分析

以苏云金芽孢杆菌科默尔亚种15A3菌株基因组DNA为模版,用touchdown PCR方法扩增几丁质酶ChiA和ChiB的全基因序列(GenBank登录号:EF103273和DQ512474)。将PCR产物连接pUCm-T克隆载体,获得重组质粒pUCm-chiA和pUCm-chiB,分别转化E.coliXL-Blue。克隆的几丁质酶基因可以利用本身的启动子异源表达各自的蛋白,不需要几丁质作为诱导物。表达的几丁质酶能够分泌到胞外。证明15A3菌株可组成型表达2种几丁质酶。经核苷酸及氨基酸序列分析证明,chiA基因全长1426bp,含有343bp的上游非编码区和1083bp的ORF,编码360个氨基酸。推测成熟蛋白分子量为36kD,只有一个几丁质酶催化域。chiB基因全长2279bp,含有248bp的上游非编码区和2031bp的ORF,编码676个氨基酸。推测成熟蛋白分子量约为70.6kD,具有三个功能域。核苷酸序列分析显示chiA和chiB的启动子所处的位置及转录起始碱基都不相同,-35区相同,而-10区有两个碱基不同,SD序列也不完全一致。     

Selection on rapidly evolving proteins in the Arabidopsis genome

Genes that have undergone positive or diversifying selection are likely to be associated with adaptive divergence between species. One indicator of adaptive selection at the molecular level is an excess of amino acid replacement fixed differences per replacement site relative to the number of synonymous fixed differences per synonymous site (omega = K(a)/K(s)). We used an evolutionary expressed sequence tag (EST) approach to estimate the distribution of omega among 304 orthologous loci between Arabidopsis thaliana and A. lyrata to identify genes potentially involved in the adaptive divergence between these two Brassicaceae species. We find that 14 of 304 genes (approximately 5%) have an estimated omega > 1 and are candidates for genes with increased selection intensities. Molecular population genetic analyses of 6 of these rapidly evolving protein loci indicate that, despite their high levels of between-species nonsynonymous divergence, these genes do not have elevated levels of intraspecific replacement polymorphisms compared to previously studied genes. A hierarchical Bayesian analysis of protein-coding region evolution within and between species also indicates that the selection intensities of these genes are elevated compared to previously studied A. thaliana nuclear loci.