A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species

Advances in next generation technologies have driven the costs of DNA sequencing down to the point that genotyping-by-sequencing (GBS) is now feasible for high diversity, large genome species. Here, we report a procedure for constructing GBS libraries based on reducing genome complexity with restriction enzymes (REs). This approach is simple, quick, extremely specific, highly reproducible, and may reach important regions of the genome that are inaccessible to sequence capture approaches. By using methylation-sensitive REs, repetitive regions of genomes can be avoided and lower copy regions targeted with two to three fold higher efficiency. This tremendously simplifies computationally challenging alignment problems in species with high levels of genetic diversity. The GBS procedure is demonstrated with maize (IBM) and barley (Oregon Wolfe Barley) recombinant inbred populations where roughly 200,000 and 25,000 sequence tags were mapped, respectively. An advantage in species like barley that lack a complete genome sequence is that a reference map need only be developed around the restriction sites, and this can be done in the process of sample genotyping. In such cases, the consensus of the read clusters across the sequence tagged sites becomes the reference. Alternatively, for kinship analyses in the absence of a reference genome, the sequence tags can simply be treated as dominant markers. Future application of GBS to breeding, conservation, and global species and population surveys may allow plant breeders to conduct genomic selection on a novel germplasm or species without first having to develop any prior molecular tools, or conservation biologists to determine population structure without prior knowledge of the genome or diversity in the species.

The limits of selection during maize domestication

The domestication of all major crop plants occurred during a brief period in human history about 10,000 years ago. During this time, ancient agriculturalists selected seed of preferred forms and culled out seed of undesirable types to produce each subsequent generation. Consequently, favoured alleles at genes controlling traits of interest increased in frequency, ultimately reaching fixation. When selection is strong, domestication has the potential to drastically reduce genetic diversity in a crop. To understand the impact of selection during maize domestication, we examined nucleotide polymorphism in teosinte branched1, a gene involved in maize evolution. Here we show that the effects of selection were limited to the gene's regulatory region and cannot be detected in the protein-coding region. Although selection was apparently strong, high rates of recombination and a prolonged domestication period probably limited its effects. Our results help to explain why maize is such a variable crop. They also suggest that maize domestication required hundreds of years, and confirm previous evidence that maize was domesticated from Balsas teosinte of southwestern Mexico.

The origin of the naked grains of maize

The most critical step in maize (Zea mays ssp. mays) domestication was the liberation of the kernel from the hardened, protective casing that envelops the kernel in the maize progenitor, teosinte. This evolutionary step exposed the kernel on the surface of the ear, such that it could readily be used by humans as a food source. Here we show that this key event in maize domestication is controlled by a single gene (teosinte glume architecture or tga1), belonging to the SBP-domain family of transcriptional regulators. The factor controlling the phenotypic difference between maize and teosinte maps to a 1-kilobase region, within which maize and teosinte show only seven fixed differences in their DNA sequences. One of these differences encodes a non-conservative amino acid substitution and may affect protein function, and the other six differences potentially affect gene regulation. Molecular evolution analyses show that this region was the target of selection during maize domestication. Our results demonstrate that modest genetic changes in single genes can induce dramatic changes in phenotype during domestication and evolution.

Investigation of the bottleneck leading to the domestication of maize

Maize (Zea mays ssp. mays) is genetically diverse, yet it is also morphologically distinct from its wild relatives. These two observations are somewhat contradictory: the first observation is consistent with a large historical population size for maize, but the latter observation is consistent with strong, diversity-limiting selection during maize domestication. In this study, we sampled sequence diversity, coupled with simulations of the coalescent process, to study the dynamics of a population bottleneck during the domestication of maize. To do this, we determined the DNA sequence of a 1,400-bp region of the Adh1 locus from 19 individuals representing maize, its presumed progenitor (Z. mays ssp. parviglumis), and a more distant relative (Zea luxurians). The sequence data were used to guide coalescent simulations of population bottlenecks associated with domestication. Our study confirms high genetic diversity in maize-maize contains 75% of the variation found in its progenitor and is more diverse than its wild relative, Z. luxurians-but it also suggests that sequence diversity in maize can be explained by a bottleneck of short duration and very small size. For example, the breadth of genetic diversity in maize is consistent with a founding population of only 20 individuals when the domestication event is 10 generations in length.

A genomics approach to the comprehensive analysis of the glutathione S-transferase gene family in soybean and maize

By BLAST searching a large expressed sequence tag database for glutathione S-transferase (GST) sequences we have identified 25 soybean (Glycine max) and 42 maize (Zea mays) clones and obtained accurate full-length GST sequences. These clones probably represent the majority of members of the GST multigene family in these species. Plant GSTs are divided according to sequence similarity into three categories: types I, II, and III. Among these GSTs only the active site serine, as well as another serine and arginine in or near the "G-site" are conserved throughout. Type III GSTs have four conserved sequence patches mapping to distinct structural features. Expression analysis reveals the distribution of GSTs in different tissues and treatments: Maize GSTI is overall the most highly expressed in maize, whereas the previously unknown GmGST 8 is most abundant in soybean. Using DNA microarray analysis we observed increased expression among the type III GSTs after inducer treatment of maize shoots, with different genes responding to different treatments. Protein activity for a subset of GSTs varied widely with seven substrates, and any GST exhibiting greater than marginal activity with chloro-2,4 dinitrobenzene activity also exhibited significant activity with all other substrates, suggesting broad individual enzyme substrate specificity.

Maize WRKY Transcription Factor ZmWRKY106 Confers Drought and Heat Tolerance in Transgenic Plants

WRKY transcription factors constitute one of the largest transcription factor families in plants, and play crucial roles in plant growth and development, defense regulation and stress responses. However, knowledge about this family in maize is limited. In the present study, we identified a drought-induced WRKY gene, ZmWRKY106, based on the maize drought de novo transcriptome sequencing data. ZmWRKY106 was identified as part of the WRKYII group, and a phylogenetic tree analysis showed that ZmWRKY106 was closer to OsWRKY13. The subcellular localization of ZmWRKY106 was only observed in the nucleus. The promoter region of ZmWRKY106 included the C-repeat/dehydration responsive element (DRE), low-temperature responsive element (LTR), MBS, and TCA-elements, which possibly participate in drought, cold, and salicylic acid (SA) stress responses. The expression of ZmWRKY106 was induced significantly by drought, high temperature, and exogenous abscisic acid (ABA), but was weakly induced by salt. Overexpression of ZmWRKY106 improved the tolerance to drought and heat in transgenic Arabidopsis by regulating stress-related genes through the ABA-signaling pathway, and the reactive oxygen species (ROS) content in transgenic lines was reduced by enhancing the activities of superoxide dismutase (SOD), peroxide dismutase (POD), and catalase (CAT) under drought stress. This suggested that ZmWRKY106 was involved in multiple abiotic stress response pathways and acted as a positive factor under drought and heat stress.

Comparative analysis of SET domain proteins in maize and Arabidopsis reveals multiple duplications preceding the divergence of monocots and dicots

Histone proteins play a central role in chromatin packaging, and modification of histones is associated with chromatin accessibility. SET domain [Su(var)3-9, Enhancer-of-zeste, Trithorax] proteins are one class of proteins that have been implicated in regulating gene expression through histone methylation. The relationships of 22 SET domain proteins from maize (Zea mays) and 32 SET domain proteins from Arabidopsis were evaluated by phylogenetic analysis and domain organization. Our analysis reveals five classes of SET domain proteins in plants that can be further divided into 19 orthology groups. In some cases, such as the Enhancer of zeste-like and trithorax-like proteins, plants and animals contain homologous proteins with a similar organization of domains outside of the SET domain. However, a majority of plant SET domain proteins do not have an animal homolog with similar domain organization, suggesting that plants have unique mechanisms to establish and maintain chromatin states. Although the domains present in plant and animal SET domain proteins often differ, the domains found in the plant proteins have been generally implicated in protein-protein interactions, indicating that most SET domain proteins operate in complexes. Combined analysis of the maize and Arabidopsis SET domain proteins reveals that duplication of SET domain proteins in plants is extensive and has occurred via multiple mechanisms that preceded the divergence of monocots and dicots.

Relationship, evolutionary fate and function of two maize co-orthologs of rice GW2 associated with kernel size and weight

BACKGROUND

In rice, the GW2 gene, found on chromosome 2, controls grain width and weight. Two homologs of this gene, ZmGW2-CHR4 and ZmGW2-CHR5, have been found in maize. In this study, we investigated the relationship, evolutionary fate and putative function of these two maize genes.

RESULTS

The two genes are located on duplicated maize chromosomal regions that show co-orthologous relationships with the rice region containing GW2. ZmGW2-CHR5 is more closely related to the sorghum counterpart than to ZmGW2-CHR4. Sequence comparisons between the two genes in eight diverse maize inbred lines revealed that the functional protein domain of both genes is completely conserved, with no non-synonymous polymorphisms identified. This suggests that both genes may have conserved functions, a hypothesis that was further confirmed through linkage, association, and expression analyses. Linkage analysis showed that ZmGW2-CHR4 is located within a consistent quantitative trait locus (QTL) for one-hundred kernel weight (HKW). Association analysis with a diverse panel of 121 maize inbred lines identified one single nucleotide polymorphism (SNP) in the promoter region of ZmGW2-CHR4 that was significantly associated with kernel width (KW) and HKW across all three field experiments examined in this study. SNPs or insertion/deletion polymorphisms (InDels) in other regions of ZmGW2-CHR4 and ZmGW2-CHR5 were also found to be significantly associated with at least one of the four yield-related traits (kernel length (KL), kernel thickness (KT), KW and HKW). None of the polymorphisms in either maize gene are similar to each other or to the 1 bp InDel causing phenotypic variation in rice. Expression levels of both maize genes vary over ear and kernel developmental stages, and the expression level of ZmGW2-CHR4 is significantly negatively correlated with KW.

CONCLUSIONS

The sequence, linkage, association and expression analyses collectively showed that the two maize genes represent chromosomal duplicates, both of which function to control some of the phenotypic variation for kernel size and weight in maize, as does their counterpart in rice. However, the different polymorphisms identified in the two maize genes and in the rice gene indicate that they may cause phenotypic variation through different mechanisms.

Maize (Zea mays L.) genome diversity as revealed by RNA-sequencing

Maize is rich in genetic and phenotypic diversity. Understanding the sequence, structural, and expression variation that contributes to phenotypic diversity would facilitate more efficient varietal improvement. RNA based sequencing (RNA-seq) is a powerful approach for transcriptional analysis, assessing sequence variation, and identifying novel transcript sequences, particularly in large, complex, repetitive genomes such as maize. In this study, we sequenced RNA from whole seedlings of 21 maize inbred lines representing diverse North American and exotic germplasm. Single nucleotide polymorphism (SNP) detection identified 351,710 polymorphic loci distributed throughout the genome covering 22,830 annotated genes. Tight clustering of two distinct heterotic groups and exotic lines was evident using these SNPs as genetic markers. Transcript abundance analysis revealed minimal variation in the total number of genes expressed across these 21 lines (57.1% to 66.0%). However, the transcribed gene set among the 21 lines varied, with 48.7% expressed in all of the lines, 27.9% expressed in one to 20 lines, and 23.4% expressed in none of the lines. De novo assembly of RNA-seq reads that did not map to the reference B73 genome sequence revealed 1,321 high confidence novel transcripts, of which, 564 loci were present in all 21 lines, including B73, and 757 loci were restricted to a subset of the lines. RT-PCR validation demonstrated 87.5% concordance with the computational prediction of these expressed novel transcripts. Intriguingly, 145 of the novel de novo assembled loci were present in lines from only one of the two heterotic groups consistent with the hypothesis that, in addition to sequence polymorphisms and transcript abundance, transcript presence/absence variation is present and, thereby, may be a mechanism contributing to the genetic basis of heterosis.

Metabolite Profiles of Maize Leaves in Drought, Heat, and Combined Stress Field Trials Reveal the Relationship between Metabolism and Grain Yield

The development of abiotic stress-resistant cultivars is of premium importance for the agriculture of developing countries. Further progress in maize (Zea mays) performance under stresses is expected by combining marker-assisted breeding with metabolite markers. In order to dissect metabolic responses and to identify promising metabolite marker candidates, metabolite profiles of maize leaves were analyzed and compared with grain yield in field trials. Plants were grown under well-watered conditions (control) or exposed to drought, heat, and both stresses simultaneously. Trials were conducted in 2010 and 2011 using 10 tropical hybrids selected to exhibit diverse abiotic stress tolerance. Drought stress evoked the accumulation of many amino acids, including isoleucine, valine, threonine, and 4-aminobutanoate, which has been commonly reported in both field and greenhouse experiments in many plant species. Two photorespiratory amino acids, glycine and serine, and myoinositol also accumulated under drought. The combination of drought and heat evoked relatively few specific responses, and most of the metabolic changes were predictable from the sum of the responses to individual stresses. Statistical analysis revealed significant correlation between levels of glycine and myoinositol and grain yield under drought. Levels of myoinositol in control conditions were also related to grain yield under drought. Furthermore, multiple linear regression models very well explained the variation of grain yield via the combination of several metabolites. These results indicate the importance of photorespiration and raffinose family oligosaccharide metabolism in grain yield under drought and suggest single or multiple metabolites as potential metabolic markers for the breeding of abiotic stress-tolerant maize.

Identification and characterization of the zinc-regulated transporters, iron-regulated transporter-like protein (ZIP) gene family in maize

BACKGROUND

Zinc (Zn) and iron (Fe) are essential micronutrients for plant growth and development, their deficiency or excess severely impaired physiological and biochemical reactions of plants. Therefore, a tightly controlled zinc and iron uptake and homeostasis network has been evolved in plants. The Zinc-regulated transporters, Iron-regulated transporter-like Proteins (ZIP) are capable of uptaking and transporting divalent metal ion and are suggested to play critical roles in balancing metal uptake and homeostasis, though a detailed analysis of ZIP gene family in maize is still lacking.

RESULTS

Nine ZIP-coding genes were identified in maize genome. It was revealed that the ZmZIP proteins share a conserved transmembrane domain and a variable region between TM-3 and TM-4. Transiently expression in onion epidermal cells revealed that all ZmZIP proteins were localized to the endoplasmic reticulum and plasma membrane. The yeast complementation analysis was performed to test the Zn or Fe transporter activity of ZmZIP proteins. Expression analysis showed that the ZmIRT1 transcripts were dramatically induced in response to Zn- and Fe-deficiency, though the expression profiles of other ZmZIP changed variously. The expression patterns of ZmZIP genes were observed in different stages of embryo and endosperm development. The accumulations of ZmIRT1 and ZmZIP6 were increased in the late developmental stages of embryo, while ZmZIP4 was up-regulated during the early development of embryo. In addition, the expression of ZmZIP5 was dramatically induced associated with middle stage development of embryo and endosperm.

CONCLUSIONS

These results suggest that ZmZIP genes encode functional Zn or Fe transporters that may be responsible for the uptake, translocation, detoxification and storage of divalent metal ion in plant cells. The various expression patterns of ZmZIP genes in embryo and endosperm indicates that they may be essential for ion translocation and storage during differential stages of embryo and endosperm development. The present study provides new insights into the evolutionary relationship and putative functional divergence of the ZmZIP gene family during the growth and development of maize.

CCCH-type zinc finger family in maize: genome-wide identification, classification and expression profiling under abscisic acid and drought treatments

Background

CCCH-type zinc finger proteins comprise a large protein family. Increasing evidence suggests that members of this family are RNA-binding proteins with regulatory functions in mRNA processing. Compared with those in animals, functions of CCCH-type zinc finger proteins involved in plant growth and development are poorly understood.

Methodology/Principal Findings

Here, we performed a genome-wide survey of CCCH-type zinc finger genes in maize (Zea mays L.) by describing the gene structure, phylogenetic relationships and chromosomal location of each family member. Promoter sequences and expression profiles of putative stress-responsive members were also investigated. A total of 68 CCCH genes (ZmC3H1-68) were identified in maize and divided into seven groups by phylogenetic analysis. These 68 genes were found to be unevenly distributed on 10 chromosomes with 15 segmental duplication events, suggesting that segmental duplication played a major role in expansion of the maize CCCH family. The Ka/Ks ratios suggested that the duplicated genes of the CCCH family mainly experienced purifying selection with limited functional divergence after duplication events. Twelve maize CCCH genes grouped with other known stress-responsive genes from Arabidopsis were found to contain putative stress-responsive cis-elements in their promoter regions. Seven of these genes chosen for further quantitative real-time PCR analysis showed differential expression patterns among five representative maize tissues and over time in response to abscisic acid and drought treatments.

Conclusions

The results presented in this study provide basic information on maize CCCH proteins and form the foundation for future functional studies of these proteins, especially for those members of which may play important roles in response to abiotic stresses.

Cloning and characterization of GLOSSY1, a maize gene involved in cuticle membrane and wax production

The cuticle covering the aerial organs of land plants plays a protective role against several biotic and abiotic stresses and, in addition, participates in a variety of plant-insect interactions. Here, we describe the molecular cloning and characterization of the maize (Zea mays) GLOSSY1 (GL1) gene, a component of the pathway leading to cuticular wax biosynthesis in seedling leaves. The genomic and cDNA sequences we isolated differ significantly in length and in most of the coding region from those previously identified. The predicted GL1 protein includes three histidine-rich domains, the landmark of a family of membrane-bound desaturases/hydroxylases, including fatty acid-modifying enzymes. GL1 expression is not restricted to the juvenile developmental stage of the maize plant, pointing to a broader function of the gene product than anticipated on the basis of the mutant phenotype. Indeed, in addition to affecting cuticular wax biosynthesis, gl1 mutations have a pleiotropic effect on epidermis development, altering trichome size and impairing cutin structure. Of the many wax biosynthetic genes identified so far, only a few from Arabidopsis (Arabidopsis thaliana) were found to be essential for normal cutin formation. Among these is WAX2, which shares 62% identity with GL1 at the protein level. In wax2-defective plants, cutin alterations induce postgenital organ fusion. This trait is not displayed by gl1 mutants, suggesting a different role of the maize and Arabidopsis cuticle in plant development.

The future of C4 research--maize, Flaveria or Cleome?

C4 photosynthesis has evolved multiple times among the angiosperms: the spatial rearrangement of the photosynthetic apparatus, combined with alterations to the leaf structure, allows CO2 to be concentrated around Rubisco. Higher CO2 concentrations at Rubisco decrease the rate of oxygenation and therefore reduce the amount of energy lost through photorespiration. C4 plants are particularly prevalent in tropical and subtropical regions because they can sustain higher rates of net photosynthesis; they also represent some of our most productive crops. To date, most progress in identifying genes crucial for C4 photosynthesis has been made using maize and Flaveria. We propose that Cleome, the most closely related genus containing C4 species to the C3 model Arabidopsis, be used together with Arabidopsis resources to accelerate our progress in elucidating the genetic basis of C4 photosynthesis.

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.

Evolution and function of the sucrose-phosphate synthase gene families in wheat and other grasses

Suc-phosphate synthase (SPS) is a key regulatory enzyme in the pathway of Suc biosynthesis and has been linked to quantitative trait loci controlling plant growth and yield. In dicotyledonous plants there are three SPS gene families: A, B, and C. Here we report the finding of five families of SPS genes in wheat (Triticum aestivum) and other monocotyledonous plants from the family Poaceae (grasses). Three of these form separate subfamilies within the previously described A, B, and C gene families, but the other two form a novel and distinctive D family, which on present evidence is only found in the Poaceae. The D-type SPS proteins lack the phosphorylation sites associated with 14-3-3 protein binding and osmotic stress activation, and the linker region between the N-terminal catalytic glucosyltransferase domain and the C-terminal Suc-phosphatase-like domain is 80 to 90 amino acid residues shorter than in the A, B, or C types. The D family appears to have arisen after the divergence of mono- and dicotyledonous plants, with a later duplication event resulting in the two D-type subfamilies. Each of the SPS gene families in wheat showed different, but overlapping, spatial and temporal expression patterns, and in most organs at least two different SPS genes are expressed. Analysis of expressed sequence tags indicated similar expression patterns to wheat for each SPS gene family in barley (Hordeum vulgare) but not in more distantly related grasses. We identified an expressed sequence tag from rice (Oryza sativa) that appears to be derived from an endogenous antisense SPS gene, and this might account for the apparently low level of expression of the related OsSPS11 sense gene, adding to the already extensive list of mechanisms for regulating the activity of SPS in plants.

Brachytic2/ZmABCB1 functions in IAA export from intercalary meristems

Dwarfism traits in Zea mays are regulated by multiple factors including the hormone auxin. Dwarf brachytic2 (br2) mutants harbour lesions in the gene encoding an orthologue of Arabidopsis thaliana ABCB1 which functions in auxin efflux out of meristematic regions in the shoot and root. br2 mesocotyls and coleoptiles exhibit reduced auxin transport. However, the dwarf stature of br2 derives from shortened lower internodes whilst the upper portion of the plant is completely normal. As such, it is counter-intuitive to attribute br2 dwarfism exclusively to reduced auxin export out of the shoot apex. Arabidopsis abcb1 mutants exhibit only minor reductions in auxin transport and plant height unless combined with mutations in the ABCB19 auxin transporter. Phylogenetic modelling analysis excludes the possibility that BR2 is more closely related to ABCB19 which has three more closely related orthologues in maize. BR2 is expressed in nodal meristems, and analyses of auxin transport and content indicate that BR2 function in these grass-specific tissues is analogous to ABCB1 function in the shoot and root apex of Arabidopsis. These results indicate that ABCB1/BR2 function is conserved between dicots and monocots, but also suggests that this function must be understood in the context of the segmental organization of grass plants.

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

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

Population structure and strong divergent selection shape phenotypic diversification in maize landraces

To conserve the long-term selection potential of maize, it is necessary to investigate past and present evolutionary processes that have shaped quantitative trait variation. Understanding the dynamics of quantitative trait evolution is crucial to future crop breeding. We characterized population differentiation of maize landraces from the State of Oaxaca, Mexico for quantitative traits and molecular markers. Qst values were much higher than Fst values obtained for molecular markers. While low values of Fst (0.011 within-village and 0.003 among-villages) suggest that considerable gene flow occurred among the studied populations, high levels of population differentiation for quantitative traits were observed (ie an among-village Qst value of 0.535 for kernel weight). Our results suggest that although quantitative traits appear to be under strong divergent selection, a considerable amount of gene flow occurs among populations. Furthermore, we characterized nonproportional changes in the G matrix structure both within and among villages that are consequences of farmer selection. As a consequence of these differences in the G matrix structure, the response to multivariate selection will be different from one population to another. Large changes in the G matrix structure could indicate that farmers select for genes of major and pleiotropic effect. Farmers' decision and selection strategies have a great impact on phenotypic diversification in maize landraces.

Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays

While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes.

Genome size in plants can vary by orders of magnitude, but this variation has long been considered to be of little functional consequence. Studying three independent adaptations to high altitude in Zea mays, we find that genome size experiences parallel pressures from natural selection, causing a reduction in genome size with increasing altitude. Though reductions in overall repetitive content are responsible for the genome size change, we find that only those individual loci contributing most to the variation in genome size are individually targeted by selection. To identify the phenotype influenced by genome size, we study how variation in genome size within a single wild population impacts leaf growth and cell division. We find that genome size variation correlates negatively with the rate of cell division, suggesting that individuals with larger genomes require longer to complete a mitotic cycle. Finally, we reanalyze data from maize inbreds to show that faster cell division is correlated with earlier flowering, connecting observed variation in genome size to an important adaptive phenotype.

Genetic association mapping identifies single nucleotide polymorphisms in genes that affect abscisic acid levels in maize floral tissues during drought

In maize, water stress at flowering causes loss of kernel set and productivity. While changes in the levels of sugars and abscisic acid (ABA) are thought to play a role in this stress response, the mechanistic basis and genes involved are not known. A candidate gene approach was used with association mapping to identify loci involved in accumulation of carbohydrates and ABA metabolites during stress. A panel of single nucleotide polymorphisms (SNPs) in genes from these metabolic pathways and in genes for reproductive development and stress response was used to genotype 350 tropical and subtropical maize inbred lines that were well watered or water stressed at flowering. Pre-pollination ears, silks, and leaves were analysed for sugars, starch, proline, ABA, ABA-glucose ester, and phaseic acid. ABA and sugar levels in silks and ears were negatively correlated with their growth. Association mapping with 1229 SNPs in 540 candidate genes identified an SNP in the maize homologue of the Arabidopsis MADS-box gene, PISTILLATA, which was significantly associated with phaseic acid in ears of well-watered plants, and an SNP in pyruvate dehydrogenase kinase, a key regulator of carbon flux into respiration, that was associated with silk sugar concentration. An SNP in an aldehyde oxidase gene was significantly associated with ABA levels in silks of water-stressed plants. Given the short range over which decay of linkage disequilibrium occurs in maize, the results indicate that allelic variation in these genes affects ABA and carbohydrate metabolism in floral tissues during drought.

Allelic and haplotypic diversity at the rp1 rust resistance locus of maize

The maize Rp1 rust resistance locus is a complex consisting of a family of closely related resistance genes. The number of Rp1 paralogs in different maize lines (haplotypes) varied from a single gene in some stocks of the inbred A188 to >50 genes in haplotypes carrying the Rp1-A and Rp1-H specificities. The sequences of paralogs in unrelated haplotypes differ, indicating that the genetic diversity of Rp1-related genes is extremely broad in maize. Two unrelated haplotypes with five or nine paralogs had identical resistance phenotypes (Rp1-D) encoded in genes that differed by three nucleotides resulting in a single amino acid substitution. Genes in some haplotypes are more similar to each other than to any of the genes in other haplotypes indicating that they are evolving in a concerted fashion.

Evolutionary Metabolomics Identifies Substantial Metabolic Divergence between Maize and Its Wild Ancestor, Teosinte

Maize (Zea mays subsp mays) was domesticated from its wild ancestor, teosinte (Zea mays subsp parviglumis). Maize's distinct morphology and adaptation to diverse environments required coordinated changes in various metabolic pathways. However, how the metabolome was reshaped since domestication remains poorly understood. Here, we report a comprehensive assessment of divergence in the seedling metabolome between maize and teosinte. In total, 461 metabolites exhibited significant divergence due to selection. Interestingly, teosinte and tropical and temperate maize, representing major stages of maize evolution, targeted distinct sets of metabolites. Alkaloids, terpenoids, and lipids were specifically targeted in the divergence between teosinte and tropical maize, while benzoxazinoids were specifically targeted in the divergence between tropical and temperate maize. To identify genetic factors controlling metabolic divergence, we assayed the seedling metabolome of a large maize-by-teosinte cross population. We show that the recent metabolic divergence between tropical and temperate maize tended to have simpler genetic architecture than the divergence between teosinte and tropical maize. Through integrating transcriptome data, we identified candidate genes contributing to metabolic divergence, many of which were under selection at the nucleotide and transcript levels. Through overexpression or mutant analysis, we verified the roles of Flavanone 3-hydroxylase1, Purple aleurone1, and maize terpene synthase1 in the divergence of their related biosynthesis pathways. Our findings not only provide important insights into domestication-associated changes in the metabolism but also highlight the power of combining omics data for trait dissection.

Comparative profiling of the sense and antisense transcriptome of maize lines

Comparative transcriptome profiling of inbred maize lines demonstrates remarkable similarities and a large number of antisense transcripts.

Background

There are thousands of maize lines with distinctive normal as well as mutant phenotypes. To determine the validity of comparisons among mutants in different lines, we first address the question of how similar the transcriptomes are in three standard lines at four developmental stages.

Results

Four tissues (leaves, 1 mm anthers, 1.5 mm anthers, pollen) from one hybrid and one inbred maize line were hybridized with the W23 inbred on Agilent oligonucleotide microarrays with 21,000 elements. Tissue-specific gene expression patterns were documented, with leaves having the most tissue-specific transcripts. Haploid pollen expresses about half as many genes as the other samples. High overlap of gene expression was found between leaves and anthers. Anther and pollen transcript expression showed high conservation among the three lines while leaves had more divergence. Antisense transcripts represented about 6 to 14 percent of total transcriptome by tissue type but were similar across lines. Gene Ontology (GO) annotations were assigned and tabulated. Enrichment in GO terms related to cell-cycle functions was found for the identified antisense transcripts. Microarray results were validated via quantitative real-time PCR and by hybridization to a second oligonucleotide microarray platform.

Conclusion

Despite high polymorphisms and structural differences among maize inbred lines, the transcriptomes of the three lines displayed remarkable similarities, especially in both reproductive samples (anther and pollen). We also identified potential stage markers for maize anther development. A large number of antisense transcripts were detected and implicated in important biological functions given the enrichment of particular GO classes.