Did gene family expansions during the Eocene–Oligocene boundary climate cooling play a role in Pooideae adaptation to cool climates?

Adaptation to cool environments is a common feature in the core group of the grass subfamily Pooideae (Triticeae and Poeae). This suggest an ancient evolutionary origin of low temperature stress tolerance dating back prior to the initiation of taxonomic divergence of core Pooideae species. Viewing the Pooideae evolution in a palaeo‐climatic perspective reveals that taxonomic divergence of the core Pooideae group initiated shortly after a global super‐cooling period at the Eocene–Oligocene boundary (～33.5–26 Ma). This global climate cooling altered distributions of plants and animals and must have imposed selection pressure for improved low temperature stress responses. Lineage‐specific gene family expansions are known to be involved in adaptation to new environmental stresses. In Pooideae, two gene families involved in low temperature stress response, the C‐repeat binding factor (CBF) and fructosyl transferase (FT) gene families, has undergone lineage‐specific expansions. We investigated the timing of these gene family expansions by molecular dating and found that Pooideae‐specific expansion events in CBF and FT gene families took place during Eocene–Oligocene super‐cooling period. We hypothesize that the E–O super‐cooling exerted selection pressure for improved low temperature stress response and frost tolerance in a core Pooideae ancestor, and that those individuals with multiple copies of CBF and FT genes were favoured.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Evolution and expression analysis of the β-glucosidase (GLU) encoding gene subfamily in maize

The maize ??-glucosidase (ZmGLU1) hydrolyzes cytokinin-conjugates for releasing active cytokinins and thus plays important roles in cytokinin regulatory processes. ZmGLU1 belongs to glycosyl hydrolases 1 (GH1) gene family with a large number of members, and the gene function of other homologs remains to be investigated. In this study, 47 Arabidopsis, 34 rice, 31 brachypodium, 28 sorghum and 26 maize GH1 protein sequences were collected and subsequently used to construct a phylogenetic tree by Neighbor-Joining method. ZmGLU1 together with its 7 paralogs and 4 sorghum homologs were assigned into a distinct group (named GLU subfamily) with far evolutionary distance to other GH1 members. None of the Arabidopsis, rice and brachypodium gene falling into this group indicated a recent evolutionary emergence of GLU subfamily in some Poaceae plants after the divergence of Poaceae species. Phylogeny and comparative genome analysis revealed that GLU subfamily members of maize and sorghum evolved from a common ancestor, and expanded independently in each species by several duplications after maize-sorghum split. Ka/Ks analysis showed that purifying selection played important roles in maintenance of similar functions among the maize GLU paralogs. In addition, the similar protein properties and cytokinin-dependent gene expressions further suggested the similar functions of ZmGLUs in cytokinin activation. However, the organ-dependent expression of ZmGLUs exhibited diverse patterns, which might contribute to their diverse roles in cytokinin homeostasis. Taken together, this work put new insights into the evolution and expression of ZmGLU genes, and provided the foundation for future functional investigations.     

Sequencing chromosome 5D of Aegilops tauschii and comparison with its allopolyploid descendant bread wheat (Triticum aestivum)

Flow cytometric sorting of individual chromosomes and chromosome‐based sequencing reduces the complexity of large, repetitive Triticeae genomes. We flow‐sorted chromosome 5D of Aegilops tauschii, the D genome donor of bread wheat and sequenced it by Roche 454 GS FLX platform to approximately 2.2x coverage. Repetitive sequences represent 81.09% of the survey sequences of this chromosome, and Class I retroelements are the prominent type, with a particular abundance of LTR/Gypsy superfamily. Nonrepetitive sequences were assembled to cover 17.76% of the total chromosome regions. Up to 6188 nonrepetitive gene loci were predicted to be encoded by the 5D chromosome. The numbers and chromosomal distribution patterns of tRNA genes suggest abundance in tRNA^Lys and tRNA^Met species, while the nonrepetitive assembly reveals tRNA^Ala species as the most abundant type. A comparative analysis of the genomic sequences of bread wheat and Aegilops chromosome 5D indicates conservation of gene content. Orthologous unique genes, matching Aegilops 5D sequences, numbered 3730 in barley, 5063 in Brachypodium, 4872 in sorghum and 4209 in rice. In this study, we provide a chromosome‐specific view into the structure and organization of the 5D chromosome of Ae. tauschii, the D genome ancestor of bread wheat. This study contributes to our understanding of the chromosome‐level evolution of the wheat genome and presents a valuable resource in wheat genomics due to the recent hybridization of Ae. tauschii genome with its tetraploid ancestor.     

Comparative mapping in grasses. Oat relationships

The development of RFLP linkage maps in hexaploid and diploid oat allows us to study genetic relationships of these species at the DNA level. In this report, we present the extension of a previously developed diploid oat map (Avena atlantica x A. hirtula) and its molecular-genetic relationships with wheat, rice and maize. Examination of 92–99% of the length of the oat genome map with probes common to Triticeae species, rice or maize showed that 84, 79 and 71%, respectively, was conserved between these species and oat. Generally, the orders of loci among chromosomes homoeologous to oat chromosomes A and D were the most conserved and those of chromosomes homoeologous to oat chromosome G were the least conserved. Conservation was observed for blocks ranging from whole chromosomes 101 cM long to small segments 2.5 cM long containing two loci. Comparison of the homoeologous segments of Triticeae, rice and maize relative to oat indicated that certain regions have been maintained in all four species. The relative positions of major genes governing traits such as seed storage proteins and resistance to leaf rusts have been conserved between cultivated oat and Triticeae species. Also, the locations of three vernalization/or photo-period response genes identified in hexaploid oat correspond to the locations of similar genes in homoeologous chromosomes of wheat, rice or maize. The locations of the centromeres for six of the seven oat chromosomes were estimated based on the homoeologous segments between oat and Triticeae chromosomes.     

Colinearity and gene density in grass genomes

Grasses are the single most important plant family in agriculture. In the past years, comparative genetic mapping has revealed conserved gene order (colinearity) among many grass species. Recently, the first studies at gene level have demonstrated that microcolinearity of genes is less conserved: small scale rearrangements and deletions complicate the microcolinearity between closely related species, such as sorghum and maize, but also between rice and other crop plants. In spite of these problems, rice remains the model plant for grasses as there is limited useful colinearity between Arabidopsis and grasses. However, studies in rice have to be complemented by more intensive genetic work on grass species with large genomes (maize, Triticeae). Gene-rich chromosomal regions in species with large genomes, such as wheat, have a high gene density and are ideal targets for partial genome sequencing.     

Uncovering the evolutionary origin of blue anthocyanins in cereal grains

Functional divergence after gene duplication plays a central role in plant evolution. Among cereals, only Hordeum vulgare (barley), Triticum aestivum (wheat) and Secale cereale (rye) accumulate delphinidin‐derived (blue) anthocyanins in the aleurone layer of grains, whereas Oryza sativa (rice), Zea mays (maize) and Sorghum bicolor (sorghum) do not. The underlying genetic basis for this natural occurrence remains elusive. Here, we mapped the barley Blx1 locus involved in blue aleurone to an approximately 1.13 Mb genetic interval on chromosome 4HL, thus identifying a trigenic cluster named MbHF35 (containing HvMYB4H, HvMYC4H and HvF35H). Sequence and expression data supported the role of these genes in conferring blue‐coloured (blue aleurone) grains. Synteny analyses across monocot species showed that MbHF35 has only evolved within distinct Triticeae lineages, as a result of dispersed gene duplication. Phylogeny analyses revealed a shared evolution pattern for MbHF35 in Triticeae, suggesting that these genes have co‐evolved together. We also identified a Pooideae‐specific flavonoid 3′,5′‐hydroxylase (F3′5′H) lineage, termed here Mo_F35H2, which has a higher amino acid similarity with eudicot F3′5′Hs, demonstrating a scenario of convergent evolution. Indeed, selection tests identified 13 amino acid residues in Mo_F35H2 that underwent positive selection, possibly driven by protein thermostablility selection. Furthermore, through the interrogation of barley germplasm there is evidence that HvMYB4H and HvMYC4H have undergone human selection. Collectively, our study favours blue aleurone as a recently evolved trait resulting from environmental adaptation. Our findings provide an evolutionary explanation for the absence of blue anthocyanins in other cereals and highlight the importance of gene functional divergence for plant diversity and environmental adaptation.     

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.     

The colinearity of the Sh2/A1 orthologous region in rice, sorghum and maize is interrupted and accompanied by genome expansion in the triticeae

The Sh2/A1 orthologous region of maize, rice, and sorghum contains five genes in the order Sh2, X1, X2, and two A1 homologs in tandem duplication. The Sh2 and A1 homologs are separated by approximately 20 kb in rice and sorghum and by approximately 140 kb in maize. We analyzed the fate of the Sh2/A1 region in large-genome species of the Triticeae (wheat, barley, and rye). In the Triticeae, synteny in the Sh2/A1 region was interrupted by a break between the X1 and X2 genes. The A1 and X2 genes remained colinear in homeologous chromosomes as in other grasses. The Sh2 and X1 orthologs also remained colinear but were translocated to a nonhomeologous chromosome. Gene X1 was duplicated on two nonhomeologous chromosomes, and surprisingly, a paralog shared homology much higher than that of the orthologous copy to the X1 gene of other grasses. No tandem duplication of A1 homologs was detected but duplication of A1 on a nonhomeologous barley chromosome 6H was observed. Intergenic distances expanded greatly in wheat compared to rice. Wheat and barley diverged from each other 12 million years ago and both show similar changes in the Sh2/A1 region, suggesting that the break in colinearity as well as X1 duplications and genome expansion occurred in a common ancestor of the Triticeae species.     

K-Box Sequence of the Apetala1-Class MADS-Box Genes in Hexaploid Wheat

Direct amplification of genomic DNA from four wheat species produced DNA fragments corresponding to the K-box sequence of the apetala1/squamosa class of the MADS-box genes. Exons 3 to 5 were highly conserved within the tribe Triticeae and very similar to the apetala1 genes of darnel ryegrass, rice, and maize. Most of the variations observed were due to synonymous substitutions: the deduced amino acid sequences were 89–99% similar within the Triticeae and 88–94% within the entire family Poaceae. Introns 3 and 4 of the apetala1 class genes were similar in wheat and rye and differed from those in other MADS-box genes presently known.     

Targeted mutagenesis of a conserved anther‐expressed P450 gene confers male sterility in monocots

Targeted mutagenesis using programmable DNA endonucleases has broad applications for studying gene function in planta and developing approaches to improve crop yields. Recently, a genetic method that eliminates the need to emasculate the female inbred during hybrid seed production, referred to as Seed Production Technology, has been described. The foundation of this genetic system relied on classical methods to identify genes critical to anther and pollen development. One of these genes is a P450 gene which is expressed in the tapetum of anthers. Homozygous recessive mutants in this gene render maize and rice plants male sterile. While this P450 in maize corresponds to the male fertility gene Ms26, male fertility mutants have not been isolated in other monocots such as sorghum and wheat. In this report, a custom designed homing endonuclease, Ems26+, was used to generate in planta mutations in the rice, sorghum and wheat orthologs of maize Ms26. Similar to maize, homozygous mutations in this P450 gene in rice and sorghum prevent pollen formation resulting in male sterile plants and fertility was restored in sorghum using a transformed copy of maize Ms26. In contrast, allohexaploid wheat plants that carry similar homozygous nuclear mutations in only one, but not all three, of their single genomes were male fertile. Targeted mutagenesis and subsequent characterization of male fertility genes in sorghum and wheat is an important step for capturing heterosis and improving crop yields through hybrid seed.     

Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.     

Comparative Genome Analysis between Agrostis stolonifera and Members of the Pooideae Subfamily,including Brachypodium distachyon

Creeping bentgrass (Agrostis stolonifera, allotetraploid 2n = 4x = 28) is one of the major cool-season turfgrasses. It is widely used on golf courses due to its tolerance to low mowing and aggressive growth habit. In this study, we investigated genome relationships of creeping bentgrass relative to the Triticeae (a consensus map of Triticum aestivum, T. tauschii, Hordeum vulgare, and H. spontaneum), oat, rice, and ryegrass maps using a common set of 229 EST-RFLP markers. The genome comparisons based on the RFLP markers revealed large-scale chromosomal rearrangements on different numbers of linkage groups (LGs) of creeping bentgrass relative to the Triticeae (3 LGs), oat (4 LGs), and rice (8 LGs). However, we detected no chromosomal rearrangement between creeping bentgrass and ryegrass, suggesting that these recently domesticated species might be closely related, despite their memberships to different Pooideae tribes. In addition, the genome of creeping bentgrass was compared with the complete genome sequence of Brachypodium distachyon in Pooideae subfamily using both sequences of the above-mentioned mapped EST-RFLP markers and sequences of 8,470 publicly available A. stolonifera ESTs (AgEST). We discovered large-scale chromosomal rearrangements on six LGs of creeping bentgrass relative to B. distachyon. Also, a total of 24 syntenic blocks based on 678 orthologus loci were identified between these two grass species. The EST orthologs can be utilized in further comparative mapping of Pooideae species. These results will be useful for genetic improvement of Agrostis species and will provide a better understanding of evolution within Pooideae species.     

The wheat D-genome HMW-glutenin locus: BAC sequencing,gene distribution,and retrotransposon clusters

A bacterial-artificial-chromosome (BAC) clone from the genome of Triticum tauschii, the D-genome ancestor of hexaploid bread wheat, was sequenced and the presence of the two paralogous x- and y-type high-molecular-weight (HMW) glutenin genes of the Glu-D1 locus was confirmed. These two genes occur in the same orientation, are 51,893 bp apart, and the separating DNA includes a 31,000-bp cluster of retrotransposons. A second retrotransposon cluster of 32,000 bp follows the x-type HMW-glutenin gene region. Each HMW-glutenin gene is found within a region of mainly unique DNA sequence which includes multiple additional genes including an active endosperm globulin gene not previously reported in the Triticeae family, a leucine-rich-repeat (LRR) type gene truncated at the 5′ end of the BAC, a kinase gene of unknown activity, remnants of a paralogous second globulin gene, and genes similar to two hypothetical rice genes. The newly identified globulin genes are assigned to a locus designated Glo-2. Comparison to available orthologous regions of the wheat A and B genomes show rapid sequence divergences flanking the HMW-glutenin genes, and the absence of two hypothetical and unknown genes found 5′ to the B-genome x-type ortholog. The region surrounding the Glu-D1 locus is similar to other reported Triticeae BAC sequences; i.e. small gene islands separated by retrotransposon clusters. Electronic Publication     

Molecular phylogeny of the genus Triticum L

The genus Triticum L. includes the major cereal crop, common or bread wheat (hexaploid Triticum aestivum L.), and other important cultivated species. Here, we conducted a phylogenetic analysis of all known wheat species and the closely related Aegilops species. This analysis was based on chloroplast matK gene comparison along with trnL intron sequences of some species. Polyploid wheat species are successfully divided only into two groups – Emmer (sections Dicoccoides and Triticum) and Timopheevii (section Timopheevii). Results reveal strictly maternal plastid inheritance of synthetic wheat amphiploids included in the study. A concordance of chloroplast origin with the definite nuclear genomes of polyploid species that were inherited at the last hybridization events was found. Our analysis suggests that there were two ancestral representatives of Aegilops speltoides Tausch that participated in the speciation of polyploid wheats with B and G genome in their genome composition. However, G genome species are younger in evolution than ones with B genome. B genome-specific PCR primers were developed for amplification of Acc-1 gene.     

BAC library development for allotetraploid Leymus (Triticeae) wildryes enable comparative genetic analysis of lax-barrenstalk1 orthogene sequences and growth habit QTLs

Tall-caespitose basin wildrye (Leymus cinereus) and rhizomatous creeping wildrye (Leymus triticoides) are perennial Triticeae relatives of wheat and barley. Quantitative trait loci (QTLs) controlling rhizome proliferation were previously detected on homoeologous regions of LG3a and LG3b in two full-sib families derived from allotetraploid hybrids. Triticeae homoeologous group 3 aligns to rice chromosome 1, which contains the rice lax panicle and maize barrenstalk1 orthogene responsible for induction of axillary branch meristems, but this gene has not been mapped or sequenced in Triticeae. We developed bacterial artificial chromosome (BAC) libraries representing 6.1 haploid equivalents of the tetraploid Leymus genome (10.7 Mb). Overgo probes designed from the lax-barrenstalk1 orthogene hybridized to 12 Leymus BAC clones. Deduced amino-acid sequences from seven BAC clones were highly conserved with the rice, maize, and sorghum lax-barrenstalk1orthogenes. Gene specific primers designed from two of the most divergent BAC clones map to homoeologous regions of Leymus LG3a and LG3b and align with the lax-barrenstalk1 orthogene on rice 1L. Comparisons of genomic DNA sequences revealed two other conserved regions surrounding the Leymus LG3a, rice, and sorghum lax-barrenstalk1 ortholoci, and one of these regions was also present in maize and Leymus LG3b sequences. Comparisons of Leymus LG3a and LG3b lax-barrenstalk1 coding sequences and flanking genomic regions elucidate molecular differences between subgenomes.     

TOOLS FOR STUDYING THE CHLOROPLAST GENOME IN THE TRITICEAE (GRAMINEAE): AN ECORI MAP,A DIAGNOSTIC DELETION,AND SUPPORT FOR BROMUS AS AN OUTGROUP

Data on variation in the chloroplast genome can be used to estimate the phylogeny of the wheat tribe (Triticeae), but two tools have been needed—a restriction site map for at least one frequent-cutting enzyme, and a well-supported outgroup. This paper presents an EcoRI map, as well as maps for six other restriction enzymes, for two members of the tribe (Elytrigia repens and Secale cereale). These are compared to the maps for Bromus tectorum, a possible outgroup, and for Briza maxima, Lolium perenne, and Festuca rubra, other members of the same subfamily (Pooideae). A unique deletion of ca. 700 base pairs was discovered that provides a distinctive marker for the cytoplasm of plants with the S genome. The deletion also appears in Dasypyrum villosum (V genome). Cladograms based on restriction site data are unequivocal in the support for Bromus tectorum as a near relative of the Triticeae; the chloroplast phylogeny is very similar to the morphological cladogram for a comparable set of taxa.     

Updating of transposable element annotations from large wheat genomic sequences reveals diverse activities and gene associations

Triticeae species (including wheat, barley and rye) have huge and complex genomes due to polyploidization and a high content of transposable elements (TEs). TEs are known to play a major role in the structure and evolutionary dynamics of Triticeae genomes. During the last 5 years, substantial stretches of contiguous genomic sequence from various species of Triticeae have been generated, making it necessary to update and standardize TE annotations and nomenclature. In this study we propose standard procedures for these tasks, based on structure, nucleic acid and protein sequence homologies. We report statistical analyses of TE composition and distribution in large blocks of genomic sequences from wheat and barley. Altogether, 3.8 Mb of wheat sequence available in the databases was analyzed or re-analyzed, and compared with 1.3 Mb of re-annotated genomic sequences from barley. The wheat sequences were relatively gene-rich (one gene per 23.9 kb), although wheat gene-derived sequences represented only 7.8% (159 elements) of the total, while the remainder mainly comprised coding sequences found in TEs (54.7%, 751 elements). Class I elements [mainly long terminal repeat (LTR) retrotransposons] accounted for the major proportion of TEs, in terms of sequence length as well as element number (83.6% and 498, respectively). In addition, we show that the gene-rich sequences of wheat genome A seem to have a higher TE content than those of genomes B and D, or of barley gene-rich sequences. Moreover, among the various TE groups, MITEs were most often associated with genes: 43.1% of MITEs fell into this category. Finally, the TRIM and copia elements were shown to be the most active TEs in the wheat genome. The implications of these results for the evolution of diploid and polyploid wheat species are discussed. Electronic Supplementary Material Supplementary material is available for this article at