On the tetraploid origin of the maize genome

Polyploidy in maize: from evolution to breeding

Polyploidy played an important role in the evolution of the three most important crops: wheat, maize and rice, each of them providing a unique model for studying allopolyploidy, segmental alloploidy or paleopolyploidy. However, its genetic and evolutionary role is still vague. The undelying mechanisms and consequences of polyploidy remain fundamental objectives in the study of eukaryotes. Maize is one of the underutilized crops at the polyploid level. This species has no stable natural polyploids, the existing ones being artificially obtained. From the experimental polyploid series of maize, only the tetraploid forms (4n = 40) are of interest. They are characterized by some valuable morphological, physiological and biochemical features, superior to the diploid forms from which they originated, but also by some drawbacks such as: reduced fertility, slower development, longer vegetation period, low productivity and adaptedness. Due to these barriers to using tetraploids in field production, maize tetraploids primarily found utility in scientific studies regarding genetic variability, inbreeding, heterosis and gene dosage effect. Since the first mention of a triploid maize plant to present, many scientists and schools, devoted their efforts to capitalize on the use of polyploidy in maize. Despite its common disadvantages as a crop, significant progress in developing tetraploid maize with good agronomic performance was achieved leading to registered tetraploid maize varieties. In this review we summarize and discuss the different aspects of polyploidy in maize, such as evolutionary context, methods of induction, morphology, fertility issue, inheritance patterns, gene expression and potential use.

Polyploid speciation in Zea (Poaceae): cytogenetic insights

MAIN CONCLUSION

The analysis of meiotic pairing affinities and genomic formulae in species and hybrids of Zea allowed us to speculate an evolutionary model to recreate the ancient polyploidization of maize and allied species. The meiotic pairing affinities and the genomic formulae analysis in Zea species and hybrids obtained in new and previous crosses, together with the molecular data known in the genus, allowed us to speculate an evolutionary model to attempt to recreate the ancient polyploidization process of Zea species. We propose that x = 5 semispecies are the ancestors of all modern species of the genus. The complex evolutionary process that originated the different taxa could be included hybridization between sympatric diploid ancestral semispecies (2n = 10) and recurrent duplication of the hybrid chromosome number, resulting in distinct auto- and allopolyploids. After the merger and doubling of independent genomes would have undergone cytological and genetical diploidization, implying revolutionary changes in genome organization and genic balance processes. Based on the meiotic behaviour of the 2n = 30 hybrids, that showed homoeology between the A subgenomes of all parental species, we propose that this subgenome A would be pivotal in all the species and would have conserved the rDNA sequences and the pairing regulator locus (PrZ). In the hypothetical model postulated here, the ancestral semispecies with the pivotal subgenome A would have had a wide geographic distribution, co-occurring and hybridizing with the semispecies harbouring B subgenomes, thus enabling sympatric speciation.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

A pan-grass transcriptome reveals patterns of cellular divergence in crops

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

Mutations in NAKED-ENDOSPERM IDD genes reveal functional interactions with SCARECROW during leaf patterning in C4 grasses

Leaves comprise a number of different cell-types that are patterned in the context of either the epidermal or inner cell layers. In grass leaves, two distinct anatomies develop in the inner leaf tissues depending on whether the leaf carries out C3 or C4 photosynthesis. In both cases a series of parallel veins develops that extends from the leaf base to the tip but in ancestral C3 species veins are separated by a greater number of intervening mesophyll cells than in derived C4 species. We have previously demonstrated that the GRAS transcription factor SCARECROW (SCR) regulates the number of photosynthetic mesophyll cells that form between veins in the leaves of the C4 species maize, whereas it regulates the formation of stomata in the epidermal leaf layer in the C3 species rice. Here we show that SCR is required for inner leaf patterning in the C4 species Setaria viridis but in this species the presumed ancestral stomatal patterning role is also retained. Through a comparative mutant analysis between maize, setaria and rice we further demonstrate that loss of NAKED-ENDOSPERM (NKD) INDETERMINATE DOMAIN (IDD) protein function exacerbates loss of function scr phenotypes in the inner leaf tissues of maize and setaria but not rice. Specifically, in both setaria and maize, scr;nkd mutants exhibit an increased proportion of fused veins with no intervening mesophyll cells. Thus, combined action of SCR and NKD may control how many mesophyll cells are specified between veins in the leaves of C4 but not C3 grasses. Together our results provide insight into the evolution of cell patterning in grass leaves and demonstrate a novel patterning role for IDD genes in C4 leaves.

Elucidation of Triacylglycerol Overproduction in the C4 Bioenergy Crop Sorghum bicolor by Constraint-Based Analysis

Upregulation of triacylglycerols (TAGs) in vegetative plant tissues such as leaves has the potential to drastically increase the energy density and biomass yield of bioenergy crops. In this context, constraint-based analysis has the promise to improve metabolic engineering strategies. Here we present a core metabolism model for the C4 biomass crop Sorghum bicolor (iTJC1414) along with a minimal model for photosynthetic CO2 assimilation, sucrose and TAG biosynthesis in C3 plants. Extending iTJC1414 to a four-cell diel model we simulate C4 photosynthesis in mature leaves with the principal photo-assimilatory product being replaced by TAG produced at different levels. Independent of specific pathways and per unit carbon assimilated, energy content and biosynthetic demands in reducing equivalents are about 1.3 to 1.4 times higher for TAG than for sucrose. For plant generic pathways, ATP- and NADPH-demands per CO2 assimilated are higher by 1.3- and 1.5-fold, respectively. If the photosynthetic supply in ATP and NADPH in iTJC1414 is adjusted to be balanced for sucrose as the sole photo-assimilatory product, overproduction of TAG is predicted to cause a substantial surplus in photosynthetic ATP. This means that if TAG synthesis was the sole photo-assimilatory process, there could be an energy imbalance that might impede the process. Adjusting iTJC1414 to a photo-assimilatory rate that approximates field conditions, we predict possible daily rates of TAG accumulation, dependent on varying ratios of carbon partitioning between exported assimilates and accumulated oil droplets (TAG, oleosin) and in dependence of activation of futile cycles of TAG synthesis and degradation. We find that, based on the capacity of leaves for photosynthetic synthesis of exported assimilates, mature leaves should be able to reach a 20% level of TAG per dry weight within one month if only 5% of the photosynthetic net assimilation can be allocated into oil droplets. From this we conclude that high TAG levels should be achievable if TAG synthesis is induced only during a final phase of the plant life cycle.

Genome wide association study reveals plant loci controlling heritability of the rhizosphere microbiome

Host genetics has recently been shown to be a driver of plant microbiome composition. However, identifying the underlying genetic loci controlling microbial selection remains challenging. Genome-wide association studies (GWAS) represent a potentially powerful, unbiased method to identify microbes sensitive to the host genotype and to connect them with the genetic loci that influence their colonization. Here, we conducted a population-level microbiome analysis of the rhizospheres of 200 sorghum genotypes. Using 16S rRNA amplicon sequencing, we identify rhizosphere-associated bacteria exhibiting heritable associations with plant genotype, and identify significant overlap between these lineages and heritable taxa recently identified in maize. Furthermore, we demonstrate that GWAS can identify host loci that correlate with the abundance of specific subsets of the rhizosphere microbiome. Finally, we demonstrate that these results can be used to predict rhizosphere microbiome structure for an independent panel of sorghum genotypes based solely on knowledge of host genotypic information.

Tissue-specific gene expression and protein abundance patterns are associated with fractionation bias in maize

BACKGROUND

Maize experienced a whole-genome duplication event approximately 5 to 12 million years ago. Because this event occurred after speciation from sorghum, the pre-duplication subgenomes can be partially reconstructed by mapping syntenic regions to the sorghum chromosomes. During evolution, maize has had uneven gene loss between each ancient subgenome. Fractionation and divergence between these genomes continue today, constantly changing genetic make-up and phenotypes and influencing agronomic traits.

RESULTS

Here we regenerate the subgenome reconstructions for the most recent maize reference genome assembly. Based on both expression and abundance data for homeologous gene pairs across multiple tissues, we observed functional divergence of genes across subgenomes. Although the genes in the larger maize subgenome are often expressing more highly than their homeologs in the smaller subgenome, we observed cases where homeolog expression dominance switches in different tissues. We demonstrate for the first time that protein abundances are higher in the larger subgenome, but they also show tissue-specific dominance, a pattern similar to RNA expression dominance. We also find that pollen expression is uniquely decoupled from protein abundance.

CONCLUSION

Our study shows that the larger subgenome has a greater range of functional assignments and that there is a relative lack of overlap between the subgenomes in terms of gene functions than would be suggested by similar patterns of gene expression and protein abundance. Our study also revealed that some reactions are catalyzed uniquely by the larger and smaller subgenomes. The tissue-specific, nonequivalent expression-level dominance pattern observed here implies a change in regulatory control which favors differentiated selective pressure on the retained duplicates leading to eventual change in gene functions.

Evolution of the RNA N 6-Methyladenosine Methylome Mediated by Genomic Duplication

RNA N ⁶-methyladenosine (m⁶A) modification is the most abundant form of RNA epigenetic modification in eukaryotes. Given that m⁶A evolution is associated with the selective constraints of nucleotide sequences in mammalian genomes, we hypothesize that m⁶A evolution can be linked, at least in part, to genomic duplication events in complex polyploid plant genomes. To test this hypothesis, we presented the maize (Zea mays) m⁶A modification landscape in a transcriptome-wide manner and identified 11,968 m⁶A peaks carried by 5,893 and 3,811 genes from two subgenomes (maize1 and maize2, respectively). Each of these subgenomes covered over 2,200 duplicate genes. Within these duplicate genes, those carrying m⁶A peaks exhibited significant differences in retention rate. This biased subgenome fractionation of m⁶A-methylated genes is associated with multiple sequence features and is influenced by asymmetric evolutionary rates. We also characterized the coevolutionary patterns of m⁶A-methylated genes and transposable elements, which can be mediated by whole genome duplication and tandem duplication. We revealed the evolutionary conservation and divergence of duplicated m⁶A functional factors and the potential role of m⁶A modification in maize responses to drought stress. This study highlights complex interplays between m⁶A modification and gene duplication, providing a reference for understanding the mechanisms underlying m⁶A evolution mediated by genome duplication events.

Exploiting protein modification systems to boost crop productivity: SUMO proteases in focus

In recent years, post-translational modification (PTM) of proteins has emerged as a key process that integrates plant growth and response to a changing environment. During the processes of domestication and breeding, plants were selected for various yield and adaptational characteristics. The post-translational modifier small ubiquitin-like modifier (SUMO) protein is known to have a role in the regulation of a number of these characteristics. Using bioinformatics, we mined the genomes of cereal and Brassica crops and their non-crop relatives Arabidopsis thaliana and Brachypodium distachyon for ubiquitin-like protease (ULP) SUMO protease sequences. We discovered that the SUMO system in cereal crops is disproportionately elaborate in comparison with that in B. distachyon. We use these data to propose deSUMOylation as a mechanism for specificity in the SUMO system.

Cytological diploidization of paleopolyploid genus Zea: Divergence between homoeologous chromosomes or activity of pairing regulator genes?

Cytological diploidization process is different in autopolyploid and allopolyploid species. Colchicine applied at the onset of meiosis suppresses the effect of pairing regulator genes resulting multivalents formation in bivalent-forming species. Colchicine treated maizes (4x = 2n = 20, A_mA_mB_mB_m) showed up to 5IV, suggesting pairing between chromosomes from genomes homoeologous A_m and B_m. In untreated individuals of the alloautooctoploid Zea perennis (8x = 2n = 40, A_pA_pA_p´A_p´B_p1B_p1B_p2B_p2) the most frequent configuration was 5IV+10II (formed by A and B genomes, respectively). The colchicine treated Z. perennis show up to 10IV revealing higher affinity within genomes A and B, but any homology among them. These results suggest the presence of a paring regulator locus (PrZ) in maize and Z. perennis, whose expression is suppressed by colchicine. It could be postulated that in Z. perennis, PrZ would affect independently the genomes A and B, being relevant the threshold of homology, the fidelity of pairing in each genomes and the ploidy level. Cytological analysis of the treated hexaploid hybrids (6x = 2n = 30), with Z. perennis as a parental, strongly suggests that PrZ is less effective in only one doses. This conclusion was reinforced by the homoeologous pairing observed in untreated dihaploid maizes, which showed up to 5II. Meiotic behaviour of individuals treated with different doses of colchicine allowed to postulate that PrZ affect the homoeologous association by controlling entire genomes (A_m or B_m) rather than individual chromosomes. Based on cytological and statistical results it is possible to propose that the cytological diploidization in Zea species occurs by restriction of pairing between homoeologous chromosomes or by genetical divergence of the homoeologous chromosomes, as was observed in untreated Z. mays ssp. parviglumis. These are independent but complementary systems and could be acting jointly in the same nucleus.

Genomic affinities revealed by GISH suggests intergenomic restructuring between parental genomes of the paleopolyploid genus Zea

The present work compares the molecular affinities, revealed by GISH, with the analysis of meiotic pairing in intra- and interspecific hybrids between species of Zea obtained in previous works. The joint analysis of these data provided evidence about the evolutionary relationships among the species from the paleopolyploid genus Zea (maize and teosintes). GISH and meiotic pairing of intraspecific hybrids revealed high genomic affinity between maize (Zea mays subsp. mays) and both Zea mays subsp. parviglumis and Zea mays subsp. mexicana. On the other hand, when Zea mays subsp. huehuetenanguensis DNA was probed on maize chromosomes, a lower affinity was detected, and the pattern of hybridization suggested intergenomical restructuring between the parental genomes of maize. When DNA from Zea luxurians was used as probe, homogeneous hybridization signals were observed through all maize chromosomes. Lower genomic affinity was observed when DNA from Zea diploperennis was probed on maize chromosomes, especially at knob regions. Maize chromosomes hybridized with Zea perennis DNA showed hybridization signals on four chromosome pairs: two chromosome pairs presented hybridization signal in only one chromosomal arm, whereas four chromosome pairs did not show any hybridization. These results are in agreement with previous GISH studies, which have identified the genomic source of the chromosomes involved in the meiotic configurations of Z. perennis × maize hybrids. These findings allow postulating that maize has a parental genome not shared with Z. perennis, and the existence of intergenomic restructuring between the parental genomes of maize. Moreover, the absence of hybridization signals in all maize knobs indicate that these heterochromatic regions were lost during the Z. perennis genome evolution.

Genome Alignment Spanning Major Poaceae Lineages Reveals Heterogeneous Evolutionary Rates and Alters Inferred Dates for Key Evolutionary Events

Multiple comparisons among genomes can clarify their evolution, speciation, and functional innovations. To date, the genome sequences of eight grasses representing the most economically important Poaceae (grass) clades have been published, and their genomic-level comparison is an essential foundation for evolutionary, functional, and translational research. Using a formal and conservative approach, we aligned these genomes. Direct comparison of paralogous gene pairs all duplicated simultaneously reveal striking variation in evolutionary rates among whole genomes, with nucleotide substitution slowest in rice and up to 48% faster in other grasses, adding a new dimension to the value of rice as a grass model. We reconstructed ancestral genome contents for major evolutionary nodes, potentially contributing to understanding the divergence and speciation of grasses. Recent fossil evidence suggests revisions of the estimated dates of key evolutionary events, implying that the pan-grass polyploidization occurred ∼96 million years ago and could not be related to the Cretaceous-Tertiary mass extinction as previously inferred. Adjusted dating to reflect both updated fossil evidence and lineage-specific evolutionary rates suggested that maize subgenome divergence and maize-sorghum divergence were virtually simultaneous, a coincidence that would be explained if polyploidization directly contributed to speciation. This work lays a solid foundation for Poaceae translational genomics.

The impact of widespread regulatory neofunctionalization on homeolog gene evolution following whole-genome duplication in maize

Whole-genome duplications are a widespread feature of plant genome evolution, having been detected in all flowering plant lineages. Despite the prevalence of these events, the extent to which duplicated genes (homeolog gene pairs) functionally diverge (neofunctionalization) is unclear. We present a genome-wide analysis of molecular evolution and regulatory neofunctionalization in maize (Zea mays L.). We demonstrate that 13% of all homeolog gene pairs in maize are regulatory neofunctionalized in leaves, and that regulatory neofunctionalized genes experience enhanced purifying selection. We show that significantly more genes have been regulatory neofunctionalized in foliar leaves than in husk leaves and that both leaf types have experienced selection for distinct functional roles. Furthermore, we demonstrate that biased subgenome expression dominance occurs only in the presence of regulatory neofunctionalization and that in nonregulatory neofunctionalized genes subgenome dominance is progressively acquired during development. Taken together, our study reveals several novel insights into the evolution of maize, genes, and gene expression, and provides a general model for gene evolution following whole-genome duplication in plants.

A near complete snapshot of the Zea mays seedling transcriptome revealed from ultra-deep sequencing

RNA-sequencing (RNA-seq) enables in-depth exploration of transcriptomes, but typical sequencing depth often limits its comprehensiveness. In this study, we generated nearly 3 billion RNA-Seq reads, totaling 341 Gb of sequence, from a Zea mays seedling sample. At this depth, a near complete snapshot of the transcriptome was observed consisting of over 90% of the annotated transcripts, including lowly expressed transcription factors. A novel hybrid strategy combining de novo and reference-based assemblies yielded a transcriptome consisting of 126,708 transcripts with 88% of expressed known genes assembled to full-length. We improved current annotations by adding 4,842 previously unannotated transcript variants and many new features, including 212 maize transcripts, 201 genes, 10 genes with undocumented potential roles in seedlings as well as maize lineage specific gene fusion events. We demonstrated the power of deep sequencing for large transcriptome studies by generating a high quality transcriptome, which provides a rich resource for the research community.

Expression and functional analysis of genes encoding cytokinin receptor-like histidine kinase in maize (Zea mays L.)

Cytokinin signaling is vital for plant growth and development which function via the two-component system (TCS). As one of the key component of TCS, transmembrane histidine kinases (HK) are encoded by a small gene family in plants. In this study, we focused on expression and functional analysis of cytokinin receptor-like HK genes (ZmHK) in maize. Firstly, bioinformatics analysis revealed that seven cloned ZmHK genes have different expression patterns during maize development. Secondly, ectopic expression by CaMV35S promoter in Arabidopsis further revealed that functional differentiation exists among these seven members. Among them, the ZmHK1a2-OX transgenic line has the lowest germination rate in the dark, ZmHK1-OX and ZmHK2a2-OX can delay leaf senescence, and seed size of ZmHK1-OX, ZmHK1a2-OX, ZmHK2-OX, ZmHK3b-OX and ZmHK2a2-OX was obviously reduced as compared to wild type. Additionally, ZmHK genes play opposite roles in shoot and root development; all ZmHK-OX transgenic lines display obvious shorter root length and reduced number of lateral roots, but enhanced shoot development compared with the wild type. Most notably, Arabidopsis response regulator ARR5 gene was up-regulated in ZmHK1-OX, ZmHK1a2-OX, ZmHK2-OX, ZmHK3b-OX and ZmHK2a2-OX as compared to wild type. Although the causal link between ZmHK genes and cytokinin signaling pathway is still an area to be further elucidated, these findings reflected that the diversification of ZmHK genes expression patterns and functions occurred in the course of maize evolution, indicating that some ZmHK genes might play different roles during maize development.

Genetic analysis of recombinant inbred lines for Sorghum bicolor × Sorghum propinquum

We describe a recombinant inbred line (RIL) population of 161 F5 genotypes for the widest euploid cross that can be made to cultivated sorghum (Sorghum bicolor) using conventional techniques, S. bicolor × Sorghum propinquum, that segregates for many traits related to plant architecture, growth and development, reproduction, and life history. The genetic map of the S. bicolor × S. propinquum RILs contains 141 loci on 10 linkage groups collectively spanning 773.1 cM. Although the genetic map has DNA marker density well-suited to quantitative trait loci mapping and samples most of the genome, our previous observations that sorghum pericentromeric heterochromatin is recalcitrant to recombination is highlighted by the finding that the vast majority of recombination in sorghum is concentrated in small regions of euchromatin that are distal to most chromosomes. The advancement of the RIL population in an environment to which the S. bicolor parent was well adapted (indeed bred for) but the S. propinquum parent was not largely eliminated an allele for short-day flowering that confounded many other traits, for example, permitting us to map new quantitative trait loci for flowering that previously eluded detection. Additional recombination that has accrued in the development of this RIL population also may have improved resolution of apices of heterozygote excess, accounting for their greater abundance in the F5 than the F2 generation. The S. bicolor × S. propinquum RIL population offers advantages over early-generation populations that will shed new light on genetic, environmental, and physiological/biochemical factors that regulate plant growth and development.

High-throughput discovery of mutations in tef semi-dwarfing genes by next-generation sequencing analysis

Tef (Eragrostis tef) is a major cereal crop in Ethiopia. Lodging is the primary constraint to increasing productivity in this allotetraploid species, accounting for losses of ∼15-45% in yield each year. As a first step toward identifying semi-dwarf varieties that might have improved lodging resistance, an ∼6× fosmid library was constructed and used to identify both homeologues of the dw3 semi-dwarfing gene of Sorghum bicolor. An EMS mutagenized population, consisting of ∼21,210 tef plants, was planted and leaf materials were collected into 23 superpools. Two dwarfing candidate genes, homeologues of dw3 of sorghum and rht1 of wheat, were sequenced directly from each superpool with 454 technology, and 120 candidate mutations were identified. Out of 10 candidates tested, six independent mutations were validated by Sanger sequencing, including two predicted detrimental mutations in both dw3 homeologues with a potential to improve lodging resistance in tef through further breeding. This study demonstrates that high-throughput sequencing can identify potentially valuable mutations in under-studied plant species like tef and has provided mutant lines that can now be combined and tested in breeding programs for improved lodging resistance.

Tackling drought stress: receptor-like kinases present new approaches

Global climate change and a growing population require tackling the reduction in arable land and improving biomass production and seed yield per area under varying conditions. One of these conditions is suboptimal water availability. Here, we review some of the classical approaches to dealing with plant response to drought stress and we evaluate how research on RECEPTOR-LIKE KINASES (RLKs) can contribute to improving plant performance under drought stress. RLKs are considered as key regulators of plant architecture and growth behavior, but they also function in defense and stress responses. The available literature and analyses of available transcript profiling data indeed suggest that RLKs can play an important role in optimizing plant responses to drought stress. In addition, RLK pathways are ideal targets for nontransgenic approaches, such as synthetic molecules, providing a novel strategy to manipulate their activity and supporting translational studies from model species, such as Arabidopsis thaliana, to economically useful crops.

RNA-seq in grain unveils fate of neo- and paleopolyploidization events in bread wheat (Triticum aestivum L.)

Background

Whole genome duplication is a common evolutionary event in plants. Bread wheat (Triticum aestivum L.) is a good model to investigate the impact of paleo- and neoduplications on the organization and function of modern plant genomes.

Results

We performed an RNA sequencing-based inference of the grain filling gene network in bread wheat and identified a set of 37,695 non-redundant sequence clusters, which is an unprecedented resolution corresponding to an estimated half of the wheat genome unigene repertoire. Using the Brachypodium distachyon genome as a reference for the Triticeae, we classified gene clusters into orthologous, paralogous, and homoeologous relationships. Based on this wheat gene evolutionary classification, older duplicated copies (dating back 50 to 70 million years) exhibit more than 80% gene loss and expression divergence while recent duplicates (dating back 1.5 to 3 million years) show only 54% gene loss and 36 to 49% expression divergence.

Conclusions

We suggest that structural shuffling due to duplicated gene loss is a rapid process, whereas functional shuffling due to neo- and/or subfunctionalization of duplicates is a longer process, and that both shuffling mechanisms drive functional redundancy erosion. We conclude that, as a result of these mechanisms, half the gene duplicates in plants are structurally and functionally altered within 10 million years of evolution, and the diploidization process is completed after 45 to 50 million years following polyploidization.

Dispersed benzoxazinone gene cluster: molecular characterization and chromosomal localization of glucosyltransferase and glucosidase genes in wheat and rye

Benzoxazinones (Bxs) are major defensive secondary metabolites in wheat (Triticum aestivum), rye (Secale cereale), and maize (Zea mays). Here, we identified full sets of homeologous and paralogous genes encoding Bx glucosyltransferase (GT) and Bx-glucoside glucosidase (Glu) in hexaploid wheat (2n = 6x = 42; AABBDD). Four GT loci (TaGTa-TaGTd) were mapped on chromosomes 7A, 7B (two loci), and 7D, whereas four glu1 loci (Taglu1a-Taglu1d) were on chromosomes 2A, 2B (two loci), and 2D. Transcript levels differed greatly among the four loci; B-genome loci of both TaGT and Taglu1 genes were preferentially transcribed. Catalytic properties of the enzyme encoded by each homeolog/paralog also differed despite high levels of identity among amino acid sequences. The predominant contribution of the B genome to GT and Glu reactions was revealed, as observed previously for the five Bx biosynthetic genes, TaBx1 to TaBx5, which are separately located on homeologous groups 4 and 5 chromosomes. In rye, where the ScBx1 to ScBx5 genes are dispersed to chromosomes 7R and 5R, ScGT and Scglu were located separately on chromosomes 4R and 2R, respectively. The dispersal of Bx-pathway loci to four distinct chromosomes in hexaploid wheat and rye suggests that the clustering of Bx-pathway genes, as found in maize, is not essential for coordinated transcription. On the other hand, barley (Hordeum vulgare) was found to lack the orthologous GT and glu loci like the Bx1 to Bx5 loci despite its close phylogenetic relationship with wheat and rye. These results contribute to our understanding of the evolutionary processes that the Bx-pathway loci have undergone in grasses.

Gene space dynamics during the evolution of Aegilops tauschii, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor genomes

Nine different regions totaling 9.7 Mb of the 4.02 Gb Aegilops tauschii genome were sequenced using the Sanger sequencing technology and compared with orthologous Brachypodium distachyon, Oryza sativa (rice), and Sorghum bicolor (sorghum) genomic sequences. The ancestral gene content in these regions was inferred and used to estimate gene deletion and gene duplication rates along each branch of the phylogenetic tree relating the four species. The total gene number in the extant Ae. tauschii genome was estimated to be 36,371. The gene deletion and gene duplication rates and total gene numbers in the four genomes were used to estimate the total gene number in each node of the phylogenetic tree. The common ancestor of the Brachypodieae and Triticeae lineages was estimated to have had 28,558 genes, and the common ancestor of the Panicoideae, Ehrhartoideae, and Pooideae subfamilies was estimated to have had 27,152 or 28,350 genes, depending on the ancestral gene scenario. Relative to the Brachypodieae and Triticeae common ancestor, the gene number was reduced in B. distachyon by 3,026 genes and increased in Ae. tauschii by 7,813 genes. The sum of gene deletion and gene duplication rates, which reflects the rate of gene synteny loss, was correlated with the rate of structural chromosome rearrangements and was highest in the Ae. tauschii lineage and lowest in the rice lineage. The high rate of gene space evolution in the Ae. tauschii lineage accounts for the fact that, contrary to the expectations, the level of synteny between the phylogenetically more related Ae. tauschii and B. distachyon genomes is similar to the level of synteny between the Ae. tauschii genome and the genomes of the less related rice and sorghum. The ratio of gene duplication to gene deletion rates in these four grass species closely parallels both the total number of genes in a species and the overall genome size. Because the overall genome size is to a large extent a function of the repeated sequence content in a genome, we suggest that the amount and activity of repeated sequences are important factors determining the number of genes in a genome.

Characterization of CENH3 proteins and centromere-associated DNA sequences in diploid and allotetraploid Brassica species

CENH3 is a centromere-specific histone H3 variant and has been used as a marker to identify active centromeres and DNA sequences associated with functional centromere/kinetochore complexes. In this study, up to four distinct CENH3 (BrCENH3) cDNAs were identified in individuals of each of three diploid species of Brassica. Comparison of the BrCENH3 cDNAs implied three related gene families: BrCENH3-A in Brassica rapa (AA), BrCENH3-B in B. nigra (BB), and BrCENH3-C in B. oleracea (CC). Each family encoded a histone fold domain and N-terminal histone tails that vary in length in all three families. The BrCENH3-B cDNAs have a deletion of two exons relative to BrCENH3-A and BrCENH3-C, consistent with the more ancient divergence of the BB genome. Chromatin immunoprecipitation and immunolabeling tests with anti-BrCENH3 antibodies indicated that both centromeric tandem repeats and the centromere-specific retrotransposons of Brassica are directly associated with BrCENH3 proteins. In three allotetraploid species, we find either co-transcription of the BrCENH3 genes of the ancestral diploid species or gene suppression of the BrCENH3 from one ancestor. Although B genome centromeres are occupied by BrCENH3-B in the ancestral species B. nigra, in allotetraploids both BrCENH3-A and BrCENH3-C proteins appear to assemble at these centromeres.

Bifurcation and enhancement of autonomous-nonautonomous retrotransposon partnership through LTR Swapping in soybean

Long terminal repeat (LTR) retrotransposons, the most abundant genomic components in flowering plants, are classifiable into autonomous and nonautonomous elements based on their structural completeness and transposition capacity. It has been proposed that selection is the major force for maintaining sequence (e.g., LTR) conservation between nonautonomous elements and their autonomous counterparts. Here, we report the structural, evolutionary, and expression characterization of a giant retrovirus-like soybean (Glycine max) LTR retrotransposon family, SNARE. This family contains two autonomous subfamilies, SARE(A) and SARE(B), that appear to have evolved independently since the soybean genome tetraploidization event approximately 13 million years ago, and a nonautonomous subfamily, SNRE, that originated from SARE(A). Unexpectedly, a subset of the SNRE elements, which amplified from a single founding SNRE element within the last approximately 3 million years, have been dramatically homogenized with either SARE(A) or SARE(B) primarily in the LTR regions and bifurcated into distinct subgroups corresponding to the two autonomous subfamilies. We uncovered evidence of region-specific swapping of nonautonomous elements with autonomous elements that primarily generated various nonautonomous recombinants with LTR sequences from autonomous elements of different evolutionary lineages, thus revealing a molecular mechanism for the enhancement of preexisting partnership and the establishment of new partnership between autonomous and nonautonomous elements.

Insights from the comparison of plant genome sequences

The next decade will see essentially completed sequences for multiple branches of virtually all angiosperm clades that include major crops and/or botanical models. These sequences will provide a powerful framework for relating genome-level events to aspects of morphological and physiological variation that have contributed to the colonization of much of the planet by angiosperms. Clarification of the fundamental angiosperm gene set, its arrangement, lineage-specific variations in gene repertoire and arrangement, and the fates of duplicated gene pairs will advance knowledge of functional and regulatory diversity and perhaps shed light on adaptation by lineages to whole-genome duplication, which is a distinguishing feature of angiosperm evolution. Better understanding of the relationships among angiosperm genomes promises to provide a firm foundation upon which to base translational genomics: the leveraging of hard-won structural and functional genomic information from crown botanical models to dissect novel and, in some cases, economically important features in many additional organisms.

A single molecule scaffold for the maize genome

About 85% of the maize genome consists of highly repetitive sequences that are interspersed by low-copy, gene-coding sequences. The maize community has dealt with this genomic complexity by the construction of an integrated genetic and physical map (iMap), but this resource alone was not sufficient for ensuring the quality of the current sequence build. For this purpose, we constructed a genome-wide, high-resolution optical map of the maize inbred line B73 genome containing >91,000 restriction sites (averaging 1 site/∼23 kb) accrued from mapping genomic DNA molecules. Our optical map comprises 66 contigs, averaging 31.88 Mb in size and spanning 91.5% (2,103.93 Mb/∼2,300 Mb) of the maize genome. A new algorithm was created that considered both optical map and unfinished BAC sequence data for placing 60/66 (2,032.42 Mb) optical map contigs onto the maize iMap. The alignment of optical maps against numerous data sources yielded comprehensive results that proved revealing and productive. For example, gaps were uncovered and characterized within the iMap, the FPC (fingerprinted contigs) map, and the chromosome-wide pseudomolecules. Such alignments also suggested amended placements of FPC contigs on the maize genetic map and proactively guided the assembly of chromosome-wide pseudomolecules, especially within complex genomic regions. Lastly, we think that the full integration of B73 optical maps with the maize iMap would greatly facilitate maize sequence finishing efforts that would make it a valuable reference for comparative studies among cereals, or other maize inbred lines and cultivars.

The maize genome contains abundant repeats interspersed by low-copy, gene-coding sequences that make it a challenge to sequence; consequently, current BAC sequence assemblies average 11 contigs per clone. The iMap deals with such complexity by the judicious integration of IBM genetic and B73 physical maps, but the B73 genome structure could differ from the IBM population because of genetic recombination and subsequent rearrangements. Accordingly, we report a genome-wide, high-resolution optical map of maize B73 genome that was constructed from the direct analysis of genomic DNA molecules without using genetic markers. The integration of optical and iMap resources with comparisons to FPC maps enabled a uniquely comprehensive and scalable assessment of a given BAC's sequence assembly, its placement within a FPC contig, and the location of this FPC contig within a chromosome-wide pseudomolecule. As such, the overall utility of the maize optical map for the validation of sequence assemblies has been significant and demonstrates the inherent advantages of single molecule platforms. Construction of the maize optical map represents the first physical map of a eukaryotic genome larger than 400 Mb that was created de novo from individual genomic DNA molecules.

Molecular and chromosomal evidence for allopolyploidy in soybean

Recent studies have documented that the soybean (Glycine max) genome has undergone two rounds of large-scale genome and/or segmental duplication. To shed light on the timing and nature of these duplication events, we characterized and analyzed two subfamilies of high-copy centromeric satellite repeats, CentGm-1 and CentGm-2, using a combination of computational and molecular cytogenetic approaches. These two subfamilies of satellite repeats mark distinct subsets of soybean centromeres and, in at least one case, a pair of homologs, suggesting their origins from an allopolyploid event. The satellite monomers of each subfamily are arranged in large tandem arrays, and intermingled monomers of the two subfamilies were not detected by fluorescence in situ hybridization on extended DNA fibers nor at the sequence level. This indicates that there has been little recombination and homogenization of satellite DNA between these two sets of centromeres. These satellite repeats are also present in Glycine soja, the proposed wild progenitor of soybean, but could not be detected in any other relatives of soybean examined in this study, suggesting the rapid divergence of the centromeric satellite DNA within the Glycine genus. Together, these observations provide direct evidence, at molecular and chromosomal levels, in support of the hypothesis that the soybean genome has experienced a recent allopolyploidization event.

Duplication of a well-conserved homeodomain-leucine zipper transcription factor gene in barley generates a copy with more specific functions

Three spikelets are formed at each rachis node of the cultivated barley (Hordeum vulgare ssp. vulgare) spike. In two-rowed barley, the central one is fertile and the two lateral ones are sterile, whereas in the six-rowed type, all three are fertile. This characteristic is determined by the allelic constitution at the six-rowed spike 1 (vrs1) locus on the long arm of chromosome 2H, with the recessive allele (vrs1) being responsible for the six-rowed phenotype. The Vrs1 (HvHox1) gene encodes a homeodomain-leucine zipper (HD-Zip) transcription factor. Here, we show that the Vrs1 gene evolved in the Poaceae via a duplication, with a second copy of the gene, HvHox2, present on the short arm of chromosome 2H. Micro-collinearity and polypeptide sequences were both well conserved between HvHox2 and its Poaceae orthologs, but Vrs1 is unique to the barley tribe. The Vrs1 gene product lacks a motif which is conserved among the HvHox2 orthologs. A phylogenetic analysis demonstrated that Vrs1 and HvHox2 must have diverged after the separation of Brachypodium distachyon from the Pooideae and suggests that Vrs1 arose following the duplication of HvHox2, and acquired its new function during the evolution of the barley tribe. HvHox2 was expressed in all organs examined but Vrs1 was predominantly expressed in immature inflorescence.

The Sorghum bicolor genome and the diversification of grasses

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the approximately 730-megabase Sorghum bicolor (L.) Moench genome, placing approximately 98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the approximately 75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization approximately 70 million years ago, most duplicated gene sets lost one member before the sorghum-rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum's drought tolerance.

The evolution of the starch biosynthetic pathway in cereals and other grasses

In most species, the precursor for starch synthesis, ADPglucose, is made exclusively in the plastids by the enzyme ADPglucose pyrophosphorylase (AGPase). However, in the endosperm of grasses, including the economically important cereals, ADPglucose is also made in the cytosol via a cytosolic form of AGPase. Cytosolic ADPglucose is imported into plastids for starch synthesis via an ADPglucose/ADP antiporter (ADPglucose transporter) in the plastid envelope. The genes encoding the two subunits of cytosolic AGPase and the ADPglucose transporter are unique to grasses. In this review, the evolutionary origins of this unique endosperm pathway of ADPglucose synthesis and its functional significance are discussed. It is proposed that the genes encoding the pathway originated from a whole-genome-duplication event in an early ancestor of the grasses.

Mutations in the rice liguleless gene result in a complete loss of the auricle, ligule, and laminar joint

The area between the upper part of the leaf sheath and the basal portion of the leaf blade contains several specialized organs, such as the laminar joint, auricle and ligule. Here we report the identification of T-DNA insertional mutant lines that lack all of these organs. The gene knocked out in the mutant lines encodes a protein that contains a SBP (SQUAMOSA promoter Binding Protein)-domain and is highly homologous to the maize LIGULELESS1 (LG1) gene. At the amino acid sequence level, the OsLG1 protein is 69% identical to maize LG1 and 78% identical to barley LG1. We named the rice gene OsLIGULELESS1 (OsLG1). Transient expression of an OsLG1:RFP (Red Fluorescent Protein) fusion protein indicated that the protein is localized to the nucleus. Transgenic plants harboring the OsLG1 promoter:GUS (beta-glucuronidase) reporter gene construct display preferential expression in developing laminar joint regions and meristemic regions. The gene is also weakly expressed in the ligule, auricles, and leaf sheaths at the basal region. These results indicate that OsLG1 is a transcriptional factor that plays an important role in building the laminar joint between leaf blade and leaf sheath boundary, thereby controlling ligule and auricle development.

Two paralogous genes encoding small subunits of ADP-glucose pyrophosphorylase in maize, Bt2 and L2, replace the single alternatively spliced gene found in other cereal species

Two types of gene encoding small subunits (SSU) of ADP-glucose pyrophosphorylase, a starch-biosynthetic enzyme, have been found in cereals and other grasses. One of these genes encodes two SSU proteins. These are targeted to different subcellular compartments and expressed in different organs of the plant: the endosperm cytosol and the leaf plastids. The SSU gene encoding two proteins evolved from an ancestral gene encoding a single protein by the acquisition of an alternative first exon. Prior to the work reported here, this type of SSU gene had been found in all grasses examined except maize. In maize, two separate genes, Bt2 and L2, were known to have the same roles as the alternatively spliced gene found in other grasses. The evolutionary origin of these maize genes and their relationship to the SSU genes in other grasses were unclear. Here we show that Bt2 and L2 are paralogous genes that arose as a result of the tetraploidization of the maize genome. Both genes derive from an ancestral alternatively spliced SSU gene orthologous to that found in other grasses. Following duplication, the Bt2 and L2 genes diverged in function. Each took a different one of the two functions of the ancestral gene. Now Bt2 encodes the endosperm cytosolic SSU but does not contribute significantly to leaf AGPase activity. Similarly, L2 has lost the use of one of its two alternative first exons. It can no longer contribute to the endosperm cytosolic SSU but is probably responsible for the bulk of the leaf AGPase SSU.

Comparative genomic analyses in Asparagus

Garden asparagus (Asparagus officinalis L.) belongs to the monocot family Asparagaceae in the order Asparagales. Onion (Allium cepa L.) and Asparagus officinalis are 2 of the most economically important plants of the core Asparagales, a well supported monophyletic group within the Asparagales. Coding regions in onion have lower GC contents than the grasses. We compared the GC content of 3374 unique expressed sequence tags (ESTs) from A. officinalis with Lycoris longituba and onion (both members of the core Asparagales), Acorus americanus (sister to all other monocots), the grasses, and Arabidopsis. Although ESTs in A. officinalis and Acorus had a higher average GC content than Arabidopsis, Lycoris, and onion, all were clearly lower than the grasses. The Asparagaceae have the smallest nuclear genomes among all plants in the core Asparagales, which typically have huge genomes. Within the Asparagaceae, European Asparagus species have approximately twice the nuclear DNA of that of southern African Asparagus species. We cloned and sequenced 20 genomic amplicons from European A. officinalis and the southern African species Asparagus plumosus and observed no clear evidence for a recent genome doubling in A. officinalis relative to A. plumosus. These results indicate that members of the genus Asparagus with smaller genomes may be useful genomic models for plants in the core Asparagales.

Molecular cloning and heterologous expression of beta1,2-xylosyltransferase and core alpha1,3-fucosyltransferase from maize

Maize is considered a promising alternative production system for pharmaceutically relevant proteins. However, like in all other plant species asparagine-linked oligosaccharides of maize glycoproteins are modified with beta1,2-xylose and core alpha1,3-fucose sugar residues, which are considered to be immunogenic in mammals. This altered N-glycosylation when compared to mammalian cells may reduce the potential of maize as a production system for heterologous glycoproteins. Here we report the cloning and characterization of the cDNA sequences coding for the maize enzymes beta1,2-xylosyltransferase (XylT) and core alpha1,3-fucosyltransferase (FucT). The cloned XylT and FucT cDNAs were shown to encode enzymatically active proteins, which were independently able to convert a mammalian acceptor glycoprotein into an antigen binding anti-plant N-glycan antibodies. The complete sequence of the XylT gene was determined. Evidence for the presence of at least three XylT and FucT gene loci in the maize genome was obtained. The identification of the two enzymes and their genes will allow the targeted downregulation or even elimination of beta1,2-xylose and core alpha1,3-fucose addition to recombinant glycoproteins produced in maize.

Differential regulation of five Pht1 phosphate transporters from maize (Zea mays L.)

Maize is one of the most important crops in the developing world, where adverse soil conditions and low fertilizer input are the two main constraints for stable food supply. Understanding the molecular and biochemical mechanisms involved in nutrient uptake is expected to support the development of future breeding strategies aimed at improving maize productivity on infertile soils. Phosphorus is the least mobile macronutrient in the soils and it is often limiting plant growth. In this work, five genes encoding Pht1 phosphate transporters which contribute to phosphate uptake and allocation in maize were identified. In phosphate-starved plants, transcripts of most of the five transporters were present in roots and leaves. Independent of the phosphate supply, expression of two genes was predominant in pollen or in roots colonized by symbiotic mycorrhizal fungi, respectively. Interestingly, high transcript levels of the mycorrhiza-inducible gene were also detectable in leaves of phosphate-starved plants. Thus, differential expression of Pht1 phosphate transporters in maize suggests involvement of the encoded proteins in diverse processes, including phosphate uptake from soil and transport at the symbiotic interface in mycorrhizas, phosphate (re)translocation in the shoot, and phosphate uptake during pollen tube growth.

Comparative physical mapping links conservation of microsynteny to chromosome structure and recombination in grasses

Nearly finished sequences for model organisms provide a foundation from which to explore genomic diversity among other taxonomic groups. We explore genome-wide microsynteny patterns between the rice sequence and two sorghum physical maps that integrate genetic markers, bacterial artificial chromosome (BAC) fingerprints, and BAC hybridization data. The sorghum maps largely tile a genomic component containing 41% of BACs but 80% of single-copy genes that shows conserved microsynteny with rice and partially tile a nonsyntenic component containing 46% of BACs but only 13% of single-copy genes. The remaining BACs are centromeric (4%) or unassigned (8%). The two genomic components correspond to cytologically discernible "euchromatin" and "heterochromatin." Gene and repetitive DNA distributions support this classification. Greater microcolinearity in recombinogenic (euchromatic) than nonrecombinogenic (heterochromatic) regions is consistent with the hypothesis that genomic rearrangements are usually deleterious, thus more likely to persist in nonrecombinogenic regions by virtue of Muller's ratchet. Interchromosomal centromeric rearrangements may have fostered diploidization of a polyploid cereal progenitor. Model plant sequences better guide studies of related genomes in recombinogenic than nonrecombinogenic regions. Bridging of 35 physical gaps in the rice sequence by sorghum BAC contigs illustrates reciprocal benefits of comparative approaches that extend at least across the cereals and perhaps beyond.

Identification and characterization of endoplasmic reticulum-associated degradation proteins differentially affected by endoplasmic reticulum stress

The disposal of misfolded proteins from the lumen of the endoplasmic reticulum (ER) is one of the quality control mechanisms present in the protein secretory pathway. Through ER-associated degradation, misfolded substrates are targeted to the cytosol where they are degraded by the proteasome. We have identified four maize (Zea mays) Der1-like genes (Zm Derlins) that encode homologs of Der1p, a yeast (Saccharomyces cerevisiae) protein implicated in ER-associated degradation. Zm Derlins are capable of functionally complementing a yeast Der1 deletion mutant. Such complementation indicates that the Der1p function is conserved among species. Zm Derlin genes are expressed at low levels throughout the plant, but appear prevalent in tissues with high activity of secretory protein accumulation, including developing endosperm cells. Expression of three of the four Zm Derlin genes increases during ER stress, with Zm Derlin1-1 showing the strongest induction. Subcellular fractionation experiments localized Zm Derlin proteins to the membrane fraction of microsomes. In maize endosperm, Zm Derlin proteins were found primarily associated with ER-derived protein bodies regardless of the presence of an ER stress response.